This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
The 24 potato cultivars grown in three regions with different cultivation environments were investigated genotype and environment interaction (GEI) using the Additive Main Effects and Multiplicative Interaction (AMMI), and Genotype main effect and genotype by environment interaction (GGE) biplot model analysis. Potato yield was highly related to environment and GEI factors but the number of tubers and dry matter content were related to the genotype factor. Potato yield was the highest in ‘Seohong’ in Gangneung, ‘Dejima’ in Jinbu area, and ‘Saebong’ in Daegwallyeong area. As a result of analysis of GEI effect between genotypes and environment, ‘Seohong’ in yield, ‘Daejima in tuber, and ‘Goun’ in dry matter content were superior and were stable. The ‘Seohong’, which has the highest yield and stability, showed promising results and can substitute the ‘Superior’ variety. The ‘Goun’ was selected as a chip potato processing grade to replace ‘Dejima’, which currently is the major processing grade. Therefore, as a result of analysis of AMMI and GGE biplot models of agricultural traits of potatoes showed that the cultivation environment has greatest influence on the yield. Therefore, proper site selection for the desired cultivar is necessary for best results. Though it is also important to select the best genotype with a high dry matter content.
서언
감자(Solanum tuberosum L.)는 전세계적으로 1.9천만 ha에서 재배되고 있는 중요한 식량작물로 연간 3.8억톤이 생산된다(FAOSTAT 2017). 우리나라에서는 주로 식용 및 가공용으로 연간 59만톤을 생산하고 있다(RDA 2016).
감자 품종은 용도에 따라 식용, 칩 가공용, 프렌치프라이 가공용 및 전분용으로 구분된다. 식용으로는 ‘수미’ 품종이 전체 면적의 60% 이상을 차지하고 있다. 칩 가공용으로는 건물함량이 높고 환원당 함량이 낮아 감자칩으로 제조하였을 때 칩의 색이 밝은 가공적성 품종으로 ‘대서’가 대표 품종이며, ‘고운’이 신품종으로 보급되고 있다(RDA 2016). 또한, 소비패턴의 변화에 따라 기능성이 강화된 보라색 감자, 병저항성 친환경 재배 감자 등 다양한 용도의 감자 품종이 활발히 개발되고 있다.
현재 국립종자원에 등록된 품종은 49품종이며, 이 중 농촌진흥청에서 개발한 품종은 17종이 있다(KSVS 2017). 각 국가별 주요 재배 품종으로는 우리나라의 경우 ‘수미(Superior)’, 미국은 가공용 품종인 ‘러셋버뱅크(Russet Burbank)’, 중동아시아 및 북아프리카에서는 고온에 강한 ‘데지레(Désirée)’ 품종이 50% 이상 재배되고 있다(van Dam 2002).
작물 품종은 서로 다른 환경에서 농업 형질이 현저한 차이를 나타낼 수 있다. 환경변화에 따른 작물 특성(performance)의 변화 양상을 유전형(Genotype)과 환경(Environment)의 상호작용(G×E interaction, GEI)이라고 하며(Vargas et al. 2001), GEI 분석을 통해 다양한 농업 환경에 적합한 품종을 선발 또는 추천할 수 있다(Nabavi et al. 2006). 또한, 품종개발단계에서 육성계통의 지역별 적응성 양상에 관한 정확한 정보를 제공함으로써 적응 지역 결정에 도움을 줄 수 있다(Van Oosterom et al. 1993, Yan 2001).
이처럼 GEI 효과에 대한 관심이 고조되면서 작물의 농업 형질에 대한 안정성 및 적응성을 정확하게 평가하기 위한 다양한 통계분석 방법이 개발되었다(Basford 1982). 그 중에서 ‘상가적인 주효과와 다양한 상호작용(Additive Main Effects and Multiplicative Interaction, AMMI)’모델과 ‘유전자형 주요 효과와 환경에 의한 유전자형(Genotype main effect and genotype by environment interaction, GGE) 행렬도 분석(biplot)’ 이 가장 많이 사용되고 있다(Gauch 1992).
AMMI 모델의 경우 수량성 등 양적 형질의 지역적응성 평가에 많이 사용되고 있으며, 상가적 요인뿐만 아니라 상승적 요인까지 해석이 가능하고 자유도가 작아서 GEI의 영향을 보다 정확하게 설명할 수 있다는 장점이 있다(Gauch 1988, Piepho 1997). 또한, 환경의 변화에 유전형의 반응이 충분하지 않은 경우에도 효과적인 분석이 가능하여(Yan & Hunt 1998, Zobel et al. 1988), 상가적인 주요 효과를 주성분 분석을 통해서 쉽게 분석할 수 있다.
All data values obtained from nearest Korea Meteorological Administration of Database. (http://www.kma.go.kr).
zEnvironment (masl) represented abbreviation of region and meter above the sea level.
yMean represented 2012-2013 climate average for cultivation period(E1, Mar. 21.-Jul. 11.; E2, Apr. 15.-Aug. 5.; E3, May 14.-Sep. 3.) including temperature and precipitation.
AMMI와 GGE biplot을 이용한 통계분석
모든 실험은 3반복으로 수행하였으며, 실험치는 평균값과 표준편차로 표시하였다. 조사된 수량성과 건물률 평가하기 위해 R analysis (R studio i386 3.3.2.)을 사용하였으며, AMMI와 GGE biplot 분석을 위해 ‘agricolae’와 ‘gge’ 라이브러리(library)를 호출하여 수량, 괴경수 및 건물률 데이터의 GEI를 시각적으로 표현하였으며 다음과 같이 계산하였다(Yan & Rajcan 2002, Park et al. 2016).
αij×βjδj=∑n=12λnξin×ηjn+εij=∑n=12ξin∓ηjn∓+εij
이 식에서 αij는 수량 등 농업 형질(trait) j에 대한 유전형(genotype) i의 평균값, βj는 형질 j에 대한 모든 유전형들의 평균값, δj는 유전형 평균값 사이의 형질 j의 표준편차, λn은 주성분의 특이치(PCn), ξin은 유전형 i에 대한 PCn 점수, ηjn은 형질 j에 대한 PCn 성적, εij는 형질 j에서 유전형 i와 연관된 잔차(residual effect)이다. 유전형 ξin와 형질 ηjn에 대한 특이벡터(singular vector)에 속하는 특이치(singular value) λn은 유전형과 형질 점수 사이에서 농업형질을 중심으로 스케일링 하여 만들었다. GT(Genotype by Trait) biplot은 각각의 유전형과 형질에 대한 주성분(Principal Component, PC) PC1과 PC2의 점수로 생성하였다(Yan & Rajcan 2002, Park et al. 2016).
각 유전형의 안정성 값인ASV(AMMI stability value)는 상호작용 주성분1(Interaction Principal Analysis 1, IPCA1)과 IPCA2 점수에 기초하여 계산하였다(Purchase 1997).
GGE biplot 분석은 유전형의 발현과 안정성을 가지는 유전형을 동시에 표현할 수 있으며, SVD(singular value decomposition) 형식으로 데이터를 분석함으로써, 각 유전형 값들이 유전형의 점수를 나타낸다(Yan & Rajcan 2002).
GGE biplot의 폴리곤 뷰(polygon view, which-won-where view)는 폴리곤면을 통과하는 수직선을 통해 biplot을 섹터로 구분할 수 있으며, 폴리곤은 biplot의 가장자리에 위치한 유전형들이 연결되어 그려진다. 환경이 포함된 섹터에서 유전형이 다른 섹터보다 유전형 발현이 높게 나타나면, GEI가 있음을 의미한다(Dia et al. 2016). 즉, 특정 환경에 적응한 유전형은 꼭지점에 위치한 유전형으로 표시된다.
결과 및 고찰
감자 품종의 재배 환경에 따른 농업적 형질 평가
감자 24품종을 대상으로 해발고도별 재배환경이 다른 3지역에서 2년간 농업적 특성을 평가한 결과는 다음과 같다(Table 3, Figs. 1-3). 재배환경별 10a당 수량성은 강릉(E1)에서 4,165 kg으로 서홍(G10), 진부(E2)에서 4,588 kg으로 대지(G13), 대관령(E3)에서 3,734 kg으로 새봉(G24)이 가장 많은 수량성을 보였다(Table 3). 식물체당 괴경수는 강릉(E1)에서 새봉(G24)이 12.6개, 진부(E2)와 대관령(E3)에서 홍영(G11)이 각각 12.2개와 11.0개로 가장 많은 괴경수를 보였다. 건물률은 강릉(E1), 진부(E2) 및 대관령(E3) 3지역 모두에서 고운(G23)이 각각 23.1%, 21.6%, 22.6%가장 높았다.
Table 3
Agronomic traits of the 24 potato genotypes across three environments during (2012-2013) main cropping season.
yE1, Gangneung; E2, Jinbu; E3, Daegwallyeong; TY, Tuber yield(kg/10a); NT, Number of tuber (ea/plant); DM, Dry matter of tuber (%).
xUnderlined values are the highest yield at each test environment.
wAverage data for three environment (E1, E2, and E3) ± standard deviation.
vC.V. Coefficient of variation
Fig. 3
Tuber dry matter (%) in the 24 potato genotypes across three environments during (2012-2013) main cropping season. G01, Irish Cobbler; G02, Superior; G03, Jopung; G04, Namsuh; G05, Jasim; G06, Jasuh; G07, Jowon; G08, Sinnamjag; G09, Haryeong; G10, Seohong; G11, Hongyoung; G12, Jayoung; G13, Dejima; G14, Chubaek; G15, Chudong; G16, Chugang; G17, Chuyoung; G18, Bangul; G19, Shepody; G20, Atlantic; G21, Gawon; G22 Gahwang; G23, Goun; G24 Saebong. Single, Single cultivation; Doble, Doble cultivation.
감자 품종의 3지역 평균 괴경 수량은 새봉(G24)이 4,118 kg/10a로 가장 많았으며, 방울(G18)이 2,257 kg/10a로 가장 적었다(Table 3, Fig. 1). 3지역 평균 괴경수는 홍영 (G11)이 식물체당 11.9개로 가장 많았으며, 자서(G06)가 5.7개로 가장 적었다(Table 3, Fig. 2). 3지역 평균 건물률은 고운(G23)이 22.4%로 가장 높았으며, 추백(G14)이 16.1%로 가장 낮았다.
3개 재배지역 중 강릉(E1)이 수량, 괴경수 및 건물률에서 표준편차(Standard Deviation, SD)와 변이계수(coefficient variation, CV)가 다른 2지역인 진부(E2), 대관령(E3)에 비해 작았다. Shim et al.(2015)은 재배환경과 GEI 효과가 미비할 경우 SD와 CV 값이 작아진다고 하였다. 본 연구에서는 수량성에서 남서(G014, 괴경수에서 홍영(G11), 건물률에서 새봉(G24)의 CV가 각각 0.171, 0.191, 0.037로 가장 낮아 재배환경에 따른 농업적 형질의 영향이 크지 않은 것은 것으로 나타났다.
농업적 형질에 대한 AMMI 모델 분석
해발고도별 재배환경에 따른 품종의 농업적 형질에 대한 GEI 효과를 구명하기 위해 분산분석과 AMMI 모델 분석을 수행하였다(Tables 4-6). 괴경 수량성, 식물체당 괴경수 및 건물률은 재배환경, 유전형 및 GEI에 고도의 유의성을 나타내었다(p<0.001). 수량성 변이는 재배환경이 전체 제곱합(Sum of Squares, SS)의 29.0%를 차지하여 가장 큰 영향을 미쳤으며, 그 다음으로 유전형과 GEI가 각각 19.2%, 10.4%를 차지하는 것으로 분석되었다(Table 4). 식물체당 괴경수의 경우 유전형이 SS의 33.6%를 차지하여 가장 큰 영향을 미치는 것으로 분석되었으며, 환경 요인은 19.7%, GEI는 10.1%인 것으로 나타났다(Table 5). 건물률은 유전형이 51.5%, 환경 요인은 14.7%, GEI는 5.8%로 유전형이 강한 영향을 미친 것으로 분석되었다(Table 6).
Table 4
AMMI analysis of variance for tuber yield (kg⋅10a-1) of 24 potato genotypes grown at three environments.
Abbreviations: df, degree of freedom; SS, sum of squares; MS, mean sum of squares.
zValues for the variables in each factor: nsnon-significant or
***Significant at p<0.001.
yValues in brackets indicate part of the E, G and G×E SS to the total yield variation.
기존 연구에서 수량성의 경우 환경요인이 51.9%, 유전형이 13.7%, GEI 효과가 31.4%으로 수량에서 환경요인이 가장 큰 영향을 미쳤다고 보고하였는데(Fufa 2013), 본 연구와 매우 유사한 경향이었다. 이러한 환경변이는 재배지역의 기후조건에 따라 달라질 수 있다고 하였다(Fufa 2013). 또한, 기존연구에서 AMMI 분석에서 환경에 대한 제곱합이 클 경우에는 GEI 효과가 있으며, 특정한 환경조건에서 최고 또는 최저의 수량성을 보일 수 있다고 하였다(Fufa 2013, Tolessa et al. 2013, Yan & Rajcan 2002). 본 연구에서도, 추백(G14)의 경우 강릉(E1)에서 수량성이 4,641 kg/10a로 높았으나, 대관령(E3)에서는 1,720 kg/10a로 낮아 GEI 효과가 큰 품종인 것으로 나타났다. 수량성과는 달리 괴경수와 건물률 경우에는 유전형이 GEI 효과보다 각각 3.3배, 8.9배 더 큰 영향을 미쳤는데, 이는 유전형이 강하게 작용하고 재배환경의 영향이 약했음을 의미한다.
감자 수량성에 대한 AMMI 주성분 분석을 수행하였을 때, 감자의 수량성에 영향을 미치는 유전형과 환경에 대한 상호작용의 정도는 주성분 1(PC1)이 64.5%, 주성분 2(PC2)가 35.8%이었다. 괴경수는 주성분 1(PC1)이 70.3%, 주성분 2(PC2)가 29.7%의 결과 값을 보였으며, 건물률은 주성분 1(PC1)이 68.4%, 주성분 2(PC2)가 31.6%인 것으로 분석되었다.
재배환경에 따른 농업형질의 특성인 수량성, 괴경수 및 건물율을 감자 24품종간 비교 검토하여 그 결과를 Table 7에 나타내었다. 수량성이 가장 높은(rank of tuber yield, rTY) 품종은 새봉(G24), 대지(G13), 서홍(G10) 순이었으나, 남작(G01)과 방울(G18)의 경우 안정성이 낮았다. 괴경수에서는 홍영(G11), 대지(G13), 자영(G12)이 괴경수가 많았으나, 자서(G06)는 괴경수가 적었다. 건물률의 경우 고운(G23), 자영(G12), 새봉(G24)이 건물율이 높았으나 추백(G14)이 건물율이 낮았다.
Table 7
Tuber yield, number, and dry matter for AMMI and rank of 24 potato genotypes tested for two years per location at three environments in Korea.
yrTY, rank of tuber yield (kg/10a); rASV_TY, rank of stability index in tuber yield.
xrNT, rank of tuber number (ea/plant); rASV_NT, rank of stability index in tuber number.
wrDM, rank of tuber dry matter (%); rASV_DM, rank of stability index in tuber dry matter.
재배환경에 따른 농업형질의 안정성(AMMI stability value, ASV)을 검토하고 그 결과를 Table 7에 나타내었다. ASV는 IPCA2 점수에 대한 IPCA1 점수의 2차원적 산점도를 바탕으로 좌표점에서 원점까지의 거리를 나타내며, 네거티브 또는 포지티브 점수에 상관없이 IPCA 점수가 클수록 특정 환경에 대한 유전형의 변이가 심하고, IPCA 점수가 작을수록 모든 환경에서 유전형 변이가 적음을 의미한다. rASV (rank of AMMI stability value) 높으면 광범위한 환경에 적응성을 가진 유전형으로 평가할 수 있으나, rASV 낮을 경우 재배환경에 따라 농업적 형질의 변이가 심한 유전형으로 분류할 수 있다고 하였다(Fufa 2013).
수량성에서 안정성이 높은 품종은 신남작(G08), 조원(G07), 가황(G22) 순이었으나, 하령(HR)과 추백의 경우 안정성이 낮았다. 괴경수에서는 홍영(G11), 대지(G13), 조원(G07)이 안정성이 높았으나, 자서(G06)는 안정성이 낮았다. 건물률의 경우 수미(G02), 대서(G20), 자영(G12)이 안정성이 높은 것으로 나타났다.
우리나라에서 식용으로 90% 이상 재배 되고 있는 수미(G02)의 경우 수량성에서 rTY가 6위로 높고 안정성에서도 rASV_TY 가 9위로 우수하였다(MFAFF 2013). 본 연구에서는 서홍(G10)이 수량이 많으면서(rTY, 3위), 안정성에서도 rASV_TY 가 6위로 우수하여 수미(G02)를 대체할 식용 품종으로 가능성이 매우 높은 것으로 판단되었다.
또한, 우리나라에서 가공용으로 대부분 이용되는 대서(G20)의 경우 건물률이 높고(rDM, 2위), 안정성에서도 rASV_DM 가 3위로 우수하여 가공용으로 오랫동안 이용되고 있다. 본 연구에서는 고운(23)이 건물률이 가장 높으면서(rDM 1위), rASV_DM 가 5위로 안정성도 높아 앞으로 대서(G20)를 대체할 가공용 품종으로 기대되었다.
농업적 형질의 AMMI biplot 분석
AMMI biplot은 농업적 형질에 대한 제 1 주성분(PCA 1)은 X축으로 제 2 주성분(PCA 2)은 Y축으로 GEI 효과를 시각적으로 보여준다(Fig. 4).
Fig. 4
AMMI biplot showing the mean (main effect) vs. stability (PC1) view of both genotypes and environments on tuber yield (A), tuber number (B), and tuber dry matter (C) in 24 potato genotypes. Abbreviations of genotypes and environments are as given in Table 1 and 2.
수량성은 강릉(E1)에서는 추백(G14), 진부(E2)는 고운(G23), 대관령(E3)은 하령(G09)이 높은 GEI 효과를 보이는 것으로 나타났다(Fig. 4A). 수량성 정도에 따라 적응성을 평가하였을 때, 수량이 낮은 그룹(2,100이상 2,700미만)에서는 신남작(G08), 홍영(G11), 세풍(G19), 가황(G22)이, 수량성이 중간인 그룹(2,700이상 3,500미만)은 조원(G07), 자심(G05)이, 수량성이 높은 그룹(3,500이상 4,300미만)은 서홍(G10)이 재배환경에 걸쳐 안정적인 적응성을 보이는 것으로 나타났다.
괴경수는 강릉(E1)에서는 새봉(G24), 대관령(E3)은 방울(G18)이 GEI 효과가 높은 것으로 나타났다(Fig. 4B). 괴경수 정도에 따른 적응성 평가에서 괴경수가 적은 그룹(4이상 7미만)에서는 자서(G06), 조원(G07), 조풍(G03)이, 괴경수가 중간 그룹(7이상 10미만)은 가원(G21), 고운(G23)이, 괴경수가 많은 그룹(10이상 13미만)에서는 대지(G13)가 재배환경에 안정적인 적응성을 보였다.
건물률은 강릉(E1)에서는 서홍(G10), 진부(E2)는 추백(G14)이 GEI 효과가 있는 것으로 나타났으나 그 효과는 앞에서 언급한 다른 두 농업적 형질에 비해 크지 않았다(Fig. 4B). 가공용 감자의 경우 건물률이 중요한 요소인데 건물률이 낮은 그룹(15.5이상 18미만)에서는 추동(G15), 서홍(G10)이, 건물률이 중간 그룹(18이상 20.5미만)은 수미(G02), 가원(G21), 추강(G16)이, 건물률이 높은 그룹(20.5이상 23미만)에서는 대지(G20), 자영(G12), 고운(G23)이 안정적인 적응성을 가지는 품종인 것으로 분석되었다.
이상의 결과로, 수량성에서는 서홍(G10), 괴경수에서 대지(G13), 건물률에서 고운(G23)이 농업적 형질이 우수하면서, 안정적인 적응성을 가지는 품종인 것으로 분석되었다.
농업적 형질의 GGE bioplot 모델 분석
GGE biplot(Genotype + Genotype and Environment)은 재배환경과 유전형 사이의 관계를 유전형 중심으로 설명하는 모델이며, 주성분 1(PC 1)과 주성분 2(PC 2)로 감자 품종별 재배환경에 대한 GEI 효과를 나타낸다. 수량성, 괴경수, 건물률을 대상으로 GGE biplot 분석을 수행하였을 때, PC 1과 PC 2 합이 각각 전체의 89, 93, 96%를 차지하여 유전형과 재배환경이 농업적 형질에 영향을 미친 것으로 나타났다(Fig. 5).
Fig. 5
Vector view of GGE (genotype main effects plus genotype x environment interaction effect) biplot of on potato tuber yield (A), tuber number (B), and tuber dry matter (C) in 24 potato genotypes across three environments.
재배환경을 가리키는 두 벡터 간의 각도가 좁을수록 재배환경간 연관성이 강하고, 재배환경 벡터와 인접한 위치에 있는 품종은 GEI 효과가 크다고 하였다(Kroonenburg 1995, Yan & Tinker 2006; Yan et al. 2007). 본 연구에서도 수량성과 괴경수의 경우 강릉(E1)과 진부(E2) 두 지역을 나타내는 벡터의 각도가 좁아 연관성이 강한 것으로 분석되었으나, 대관령(E3)은 다른 두 재배환경과 이루는 각이 매우 커서 독립적 환경인 것으로 나타났다. 그러나, 건물률은 수량성, 괴경수와 달리 3지역을 나타내는 벡터의 각도가 좁은 것으로 나타나 재배환경에 큰 영향을 받지 않음을 알 수 있었다. 또한, GGE biplot분석을 통해 지역적응성 평가시 여러 지역 중 유사 지역 중 한 지역만 품종특성 평가에 큰 문제가 없어 품종개발에 필요한 재배시험의 규모를 줄일 수 있다고 하였다(Choukan 2011). 본 연구에서도 수량성에서 강릉(E1)과 진부(E2) 지역이 유사성이 높은 것으로 분류되었다.
GGE biplot은 다양한 환경내에서 유전형의 발현 정도와 안정성을 동시에 표현함으로써 유전형간 비교가 용이하고 중심점에서 멀리 떨어져 있는 유전형은 농업형질의 발현에 미치는 특정 환경의 영향이 크다는 것을 의미한다(Yan et al. 2007, Yan et al. 2000). 본 연구에서는 수량성에서 새봉(G24), 괴경에서 홍영(G11), 건물률에서 신남작(08)의 유전형은 특정 환경에서 발현이 높았다(Fig. 5).
또한, GGE biplot 분석은 메가환경에서 안정적인 품종을 선발할 수 있는데 다각형 모양의 폴리곤 뷰(polygon view, which- won-where view)는 적지적품종 선정에 유용하게 활용될 수 있다(Gauch & Zobel 1997). 괴경수와 건물률의 경우 폴리곤뷰가 한 섹터의 환경을 가지고 있어 환경적인 영향보다 유전형이 큰 영향을 미치는 것으로 확인되었다(Fig. 6).
Fig. 6
Polygon (which-won-where) view of GGE (genotype main effects plus genotype x environment interaction effect) biplot of on potato tuber yield (A), tuber number (B), and tuber dry matter (C) in 24 potato genotypes across three environments.
폴리곤뷰는 biplot 의 가장 자리에 위치한 유전형들을 서로 선(segment)으로 연결하여 표시함으로 유전형이 다각형 내에 포함된다. 각각의 유전형을 잇는 선은 다각형의 측면 또는 그 연장 부분에 수직으로 형성되어 biplot을 여러 부분으로 나눌 수 있다(Dia et al. 2016). 수량성의 경우 총 5개의 섹터(sector)로 구분되었으며, 각 섹터에서 메가환경의 특성을 나타내는 유전형이 분포하였다. 첫 번째 섹터는 방울(G18)이 하나의 메가환경을 나타내고, 두 번째 섹터는 하령(G09), 세 번째 섹터는 새봉(G24), 네 번째 섹터는 대지(G13), 마지막 섹터는 추백(G14)이 독립적인 메가환경으로 구분되었다. 각각의 메가환경의 꼭지점에 위치한 유전형들은 각각의 환경에서 최고의 성능을 보이므로 대지(G13)는 강릉(E1), 새봉(G24)는 대관령(E3)에서 최고의 수량성을 나타내었다.
이상의 결과, AMMI와 GGE biplot는 유전형의 발현이나 안정성 분석에 효과적으로 활용 가능하였다.
이러한 결과를 바탕으로 감자의 수량성은 재배환경이 큰 영향을 미치므로 품종별 최적지 선정을 통해 재배하는 것이 바람직하고, 건물률은 유전적 요인이 크므로 가공용 품종육성을 위해서는 건물률이 높은 품종을 교배하기 위한 양친으로 활용하여 육성하는 것이 필요하다고 판단되었다.
적요
우리나라 주요 감자 품종 24품종을 대상으로, 해발고도별 재배환경이 다른 3개 지역에서 2년 동안 수량성, 괴경수, 건물률 등 농업 형질을 조사하고, 유전형(G)과 재배환경(E)과의 상호작용(G×E interaction, GEI) 효과를 구명하기 위해 AMMI와 GGE biplot 모델 분석을 수행하였다. 감자 수량성은 환경과 상호작용이 차지하는 비율이 높고, 괴경수와 건물률은 유전형 효과가 차지하는 비율이 높았다. 재배환경별 수량성은 강릉 지역에서 ‘서홍’, 진부 지역에서 ‘대지’, 대관령 지역에서 ‘새봉’이 높았다. 유전형과 환경과의 GEI 효과를 분석한 결과로서, 수량성에서 ‘서홍’, 괴경수에서 ‘대지’, 건물률에서 ‘고운’이 안정성이 높았다. 우리나라에서 ‘수미’ 대체 품종으로 수량성과 안정성이 높은 ‘서홍’이 유망하였고, 가공용 주요 품종인 ‘대서’를 대체할 품종으로 ‘고운’이 선발되었다. 따라서, 감자의 농업적 형질의 AMMI 과 GGE biplot 모델 분석 결과, 수량성은 재배환경이 큰 영향을 미치므로 품종별 재배적지 추천이 필요하고, 건물률은 유전적 요인이 커 가공용 품종육성을 위해서는 건물률이 높은 계통을 선발하는 것이 바람직한 것으로 판단되었다.
사사
본 논문은 농촌진흥청 연구사업(과제번호: PJ01024002)의 지원을 받아 수행되었습니다.
References
Basford KE. 1982. The use of multidimensional scaling in analyzing multi-attribute genotype response across environments. Aust J Agric Res 33: 473-480.
Choukan R. 2011. Genotype, environment and genotype ×environment interaction effects on the performance of maize (Zea mays L.). inbred lines. Crop Breed J. 1: 97-103. DOI: http://dx.doi.org/10.22092/cbj.2011.100349
Dia M, Wehner TC, Hassell R, Price DS, Boyhan GE, Olson S, King S, Davis AR, Tolla GE. 2016. Genotype ×environment interaction and stability analysis for watermelon fruit yield in the United States. Crop Sci 56: 1645-1661.
Kroonenberg PM. 1995. Introduction to biplots for GE tables. Res. Rep. 51: Dep. of Mathematics, Univ. of Queensland, Brisbane, Australia.
Ministry for Food, Agriculture, Forestry and Fisheries (MFAFF)2013. Production and supply plan of major crop seeds in summer for 2013. Ministry of Agriculture, Forestry, and Fisheries 57.
National Climate Data Service System (NCDSS) 2016. http://sts.kma.go.kr. (in Korean).
Navabi A, Yang RC, Helm CJ, Spawer DM. 2006. Can spring wheat growing mega-environments in the northern great plain be dissected for representative locations or niche-adapted genotypes? Crop Sci 46: 1107-1116.
Park MW. 2017. Study on effect of genotype and environment interaction on tomato fruit characteristics. Ph.D. thesis. Gyeongsang National Univ.
Park MW, Chung YS, Lee SH. 2016. Quality changes in tomato fruits caused by genotype and environment interactions. Horti Sci Technol 35: 361-372.
Purchase JL. 1997. Parametric analysis to described G x E interaction and yield stability in winter yield. Ph.D Thesis Department of Agronomy, Faculty of Agriculture, University of Orange Free State. Bloemfontein, South Africa, 4-83.
Rural Development Administration (RDA).2014. The Potato: Guidance of agricultural technology-31. Rural Development Administration, Suwon, Korea: pp. 24-120.
Rural Development Administration (RDA). 2016. Standard farming text -The potato. Rural Development Administration, Suwon, Korea: pp. 16-290.
Shim KB, Shin SH, Shon JY, Kang SG, Yang WH, Heu SG. 2015. Interpretation of genotype × environment interaction of sesame yield using GGE biplot analysis. Korean J Crop Sci 6: 349-354.
Tolessa TT, Keneni G, Sefera T, Jarso M, Bekele Y. 2013. Genotype × environment interaction and performance stability for grain yield in field pea (Pisum sativum L.). genotypes. Int J Plant Breed 7: 117-123.
van Dam J. 2002. Genetic characterisation of agronomic and morphological traits and the development of DNA markers associated with total glycoalkaloid content in the tubers of tetraploid potato (Solanum tuberosum L.). Ph. D Thesis Wageningen University.
van Oosterom EJ, Kleijn D, Ceccarelli S, Nachit MM. 1993. Genotype-by-Environment Interactions of Barley in the Mediterranean Region. Crop Sci 33: 669-674. doi: 10.2135/cropsci1993.0011183X003300040004x.
Yan W. 2001. GGE biplot : a window application for graphical analysis of multi-environmental data and other types of two-way data. Agron J 93: 1111-1118.
Yan W, Hunt LA. 1998. Genotype by environment interaction and crop yield. Plant Breed Rev 16: 135-178.
Interpretation of Genotype × Environment Interaction in Agronomic Traits in 24 Potato Cultivars Using AMMI and GGE Biplot Analysis
Fig. 1
Tuber yield (kg⋅10a-1) in the 24 potato genotypes across three environments during (2012-2013) main cropping season. G01, Irish Cobbler; G02, Superior; G03, Jopung; G04, Namsuh; G05, Jasim; G06, Jasuh; G07, Jowon; G08, Sinnamjag; G09, Haryeong; G10, Seohong; G11, Hongyoung; G12, Jayoung; G13, Dejima; G14, Chubaek; G15, Chudong; G16, Chugang; G17, Chuyoung; G18, Bangul; G19, Shepody; G20, Atlantic; G21, Gawon; G22 Gahwang; G23, Goun; G24 Saebong. Single, Single cultivation; Doble, Doble cultivation.
Fig. 2
Number of tuber (ea/plant) in the 24 potato genotypes across three environments during (2012-2013) main cropping season. G01, Irish Cobbler; G02, Superior; G03, Jopung; G04, Namsuh; G05, Jasim; G06, Jasuh; G07, Jowon; G08, Sinnamjag; G09, Haryeong; G10, Seohong; G11, Hongyoung; G12, Jayoung; G13, Dejima; G14, Chubaek; G15, Chudong; G16, Chugang; G17, Chuyoung; G18, Bangul; G19, Shepody; G20, Atlantic; G21, Gawon; G22 Gahwang; G23, Goun; G24 Saebong. Single, Single cultivation; Doble, Doble cultivation.
Fig. 3
Tuber dry matter (%) in the 24 potato genotypes across three environments during (2012-2013) main cropping season. G01, Irish Cobbler; G02, Superior; G03, Jopung; G04, Namsuh; G05, Jasim; G06, Jasuh; G07, Jowon; G08, Sinnamjag; G09, Haryeong; G10, Seohong; G11, Hongyoung; G12, Jayoung; G13, Dejima; G14, Chubaek; G15, Chudong; G16, Chugang; G17, Chuyoung; G18, Bangul; G19, Shepody; G20, Atlantic; G21, Gawon; G22 Gahwang; G23, Goun; G24 Saebong. Single, Single cultivation; Doble, Doble cultivation.
Fig. 4
AMMI biplot showing the mean (main effect) vs. stability (PC1) view of both genotypes and environments on tuber yield (A), tuber number (B), and tuber dry matter (C) in 24 potato genotypes. Abbreviations of genotypes and environments are as given in Table 1 and 2.
Fig. 5
Vector view of GGE (genotype main effects plus genotype x environment interaction effect) biplot of on potato tuber yield (A), tuber number (B), and tuber dry matter (C) in 24 potato genotypes across three environments.
Fig. 6
Polygon (which-won-where) view of GGE (genotype main effects plus genotype x environment interaction effect) biplot of on potato tuber yield (A), tuber number (B), and tuber dry matter (C) in 24 potato genotypes across three environments.
Fig. 1
Fig. 2
Fig. 3
Fig. 4
Fig. 5
Fig. 6
Interpretation of Genotype × Environment Interaction in Agronomic Traits in 24 Potato Cultivars Using AMMI and GGE Biplot Analysis
Description of the potato genotypes evaluated in this study.
Genotype
Cultivar
Released year
Utilization
Cropping systemz
Total yieldy(kg/10a)
Breeding nation
Cross combination
Maturity
G01
Irish Cobbler
1960
Table potato
Single
3,761
USA
bud-mutant of Early Rose
Early
G02
Superior
1978
Table potato, Potato chip
Single
3,191
USA
B96-56 × M59-44
Early
G03
Jopung
1988
Table potato
Single
3,510
Korea
Resy × Superior
Early
G04
Namsuh
1995
Table potato
Single
4,714
Korea
(Sroeken × Superior) × Wheeler
Medium
G05
Jasim
1999
Table potato
Single
2,733
Korea
Daekwan48 × B6603-6
Late
G06
Jasuh
2000
Table potato
Single
3,770
Korea
Daekwan49 × B6603-6
Early
G07
Jowon
2000
Table potato
Single
3,703
Korea
AG54071 × Katahdin
Early
G08
Sinnamjag
2001
Table potato
Single
4,444
Korea
Irish Cobbler × Katahdin
Early
G09
Haryeong
2005
Table potato
Single
3,899
Korea
Atlantic × Superior
Early-medium
G10
Seohong
2006
Table potato
Single
4,052
Korea
Jasim × 93K61-5
Early-medium
G11
Hongyoung
2007
Table potato, Potato chip
Single
3,112
Korea
Atlantic × AG34314
Medium-late
G12
Jayoung
2007
Table potato, Potato chip
Single
3,096
Korea
Atlantic × AG34314
Late
G13
Dejima
1978
Table potato
Double
S 2,956 A 2,359
Japan
Hokaido31 × Unzen
Medium-late
G14
Chubaek
1999
Table potato
Double
S 3,529 A 2,340
Korea
H83011-3 × Superior
Early
G15
Chudong
2001
Table potato
Double
S 3,243 A 2,181
Korea
H83520 × Superior
Medium
G16
Chugang
2002
Table potato
Double
S 3,462 A 2,746
Korea
H83005-2 × Superior
Medium
G17
Chuyoung
2004
Table potato
Double
S 2,952 A 2,057
Korea
Atlantic × (Russet Burbank ×Solanum phureja)
Medium
G18
Bangul
2010
Table potato
Double
S 4,040 A 3,053
Korea
Dejima × Superior
Early
G19
Shepody
1988
French fry
Single
3,684
Canada
F58050 × Wauseon
Medium
G20
Atlantic
1995
Potato chip
Single
2,934
USA
B5141-6 × Wauseon
Medium
G21
Gawon
1999
Potato chip
Single
3,173
Korea
Namsuh × Konahubuki
Early
G22
Gahwang
2002
Potato chip
Single
3,438
Korea
Atlantic × AG54701
Early
G23
Goun
2006
Potato chip
Double
S 3,322 A 2,622
Korea
Lemhi Russet × Chubaek
Early-medium
G24
Saebong
2010
Potato chip
Double
S 3,200 A 2,710
Korea
Daegwan2-8 × DH09-231
Early
zSingle, Single cultivation; Double, Double cultivation.
yS, Spring crop; A, Autumn crop.
Meteorological data of three environments of potato cultivation during 2012 and 2013 cropping season.
Environmentz
Factor
Mar.
Apr.
May
Jun.
Jul.
Aug.
Sep.
Meanz
E1(20 m)
Maximum temp. (°C)
12.2y
16.2
21.8
23.9
30.0
-
-
20.8
Minimum temp. (°C)
2.6
6.1
12.3
16.8
21.8
-
-
11.9
Average temp. (°C)
7.5
11.1
16.9
19.9
26.1
-
-
16.3
Accumulated precipitation (mm)
47.8
82.0
77.1
28.7
67.3
-
-
60.6
Sum of sunshine (hr)
151.9
231.2
245.7
155.3
43.8
-
-
165.6
Accumulated precipitation days
5.5
9.0
8.5
11.0
4.5
-
-
7.7
E2(560 m)
Maximum temp. (°C)
-
16.6
22.1
25.4
27.0
28.2
-
23.9
Minimum temp. (°C)
-
3.0
7.7
13.7
19.0
20.0
-
12.7
Average temp. (°C)
-
9.3
14.5
18.7
22.6
24.2
-
17.9
Accumulated precipitation (mm)
-
57.0
76.3
89.3
416.3
40.0
-
135.8
Sum of sunshine (hr)
-
-
-
- -
- -
-
-
-
Accumulated precipitation
-
4.0
7.5
9.0
15.5
3.0
-
7.8
E3(800 m)
Maximum temp. (°C)
-
-
20.2
21.7
25.2
25.8
23.1
23.2
Minimum temp. (°C)
-
-
8.0
12.6
18.3
18.0
14.1
14.2
Average temp. (°C)
-
-
14.2
16.8
21.4
21.5
17.6
18.3
Accumulated precipitation (mm)
-
-
143.3
94.3
330.6
188.6
4.0
152.1
Sum of sunshine (hr)
-
-
56.9
179.6
180.4
187.2
13.2
123.5
Accumulated precipitation
-
-
4.0
13.0
16.0
14.0
1.0
9.6
All data values obtained from nearest Korea Meteorological Administration of Database. (http://www.kma.go.kr).
zEnvironment (masl) represented abbreviation of region and meter above the sea level.
yMean represented 2012-2013 climate average for cultivation period(E1, Mar. 21.-Jul. 11.; E2, Apr. 15.-Aug. 5.; E3, May 14.-Sep. 3.) including temperature and precipitation.
Agronomic traits of the 24 potato genotypes across three environments during (2012-2013) main cropping season.
yrTY, rank of tuber yield (kg/10a); rASV_TY, rank of stability index in tuber yield.
xrNT, rank of tuber number (ea/plant); rASV_NT, rank of stability index in tuber number.
wrDM, rank of tuber dry matter (%); rASV_DM, rank of stability index in tuber dry matter.
Table 1 Description of the potato genotypes evaluated in this study.
Single, Single cultivation; Double, Double cultivation.
S, Spring crop; A, Autumn crop.
Table 2 Meteorological data of three environments of potato cultivation during 2012 and 2013 cropping season.
All data values obtained from nearest Korea Meteorological Administration of Database. (http://www.kma.go.kr).
Environment (masl) represented abbreviation of region and meter above the sea level.
Mean represented 2012-2013 climate average for cultivation period(E1, Mar. 21.-Jul. 11.; E2, Apr. 15.-Aug. 5.; E3, May 14.-Sep. 3.) including temperature and precipitation.
Table 3 Agronomic traits of the 24 potato genotypes across three environments during (2012-2013) main cropping season.
