80份陸稻晚稻農(nóng)家種資源的遺傳多樣性綜合分析

徐志軍 　李亞波　歐陽紅軍　徐磊　安東升　劉洋

摘 ?要：為保護和利用陸稻農(nóng)家種資源，本研究利用SNP標(biāo)記和表型鑒定對80份陸稻晚稻農(nóng)家種資源進(jìn)行遺傳多樣性綜合分析。結(jié)果表明：利用水稻1KSNP芯片對80份資源進(jìn)行檢測后共獲得739個SNP的基因型數(shù)據(jù)，基于SNP的群體結(jié)構(gòu)分析、PCA分析和聚類分析，均將80份資源分為云南亞群（63份）、海南亞群（11份）和貴州亞群（6份）共3個亞群，3個亞群間的平均遺傳距離分別為0.83、0.63、0.80，其中海南亞群與其他2個亞群的遺傳距離較遠(yuǎn)。貴州的陸稻資源分別屬于貴州亞群和云南亞群，表明黔南地區(qū)的陸稻資源可能存在不同的來源。云南亞群可進(jìn)一步分為3個亞群，平均遺傳距離分別為0.30、0.28和0.22，表明細(xì)分亞群間的遺傳距離很近。表型鑒定結(jié)果表明，11個表型性狀的變異系數(shù)在4.09%～38.77%之間，其中有效分蘗數(shù)和穗重的變異系數(shù)較大，均超過25.00%，株高、種子長、種子寬、種子長寬比和種子圓度均值變異系數(shù)較小，均低于10%;遺傳多樣性指數(shù)范圍為1.76～2.84，種子長的遺傳多樣性最豐富。相關(guān)分析結(jié)果表明，11個表型性狀間存在著復(fù)雜的相關(guān)關(guān)系，其中種子長寬比與圓度均值的相關(guān)系數(shù)最高為0.98。在表型PCA分析中，前5個因子累計貢獻(xiàn)率為82.37%，可反映表型性狀的絕大部分信息，各成分的貢獻(xiàn)率分別為29.08%、23.42%、12.35%、9.77%和7.75%，主要反映了株型和產(chǎn)量相關(guān)信息。80份資源的綜合得分范圍為0.15～0.73，排名前10的種質(zhì)資源（ZRG36、ZRG101、ZRG58、ZRG31、ZRG1、ZRG102、ZRG104、ZRG63、ZRG29、ZRG11）均來自海拔較高地區(qū)，其中2份來源于海南，5份來源于貴州，3份來源于云南，且在穗長、種子長、種子寬、穗重和千粒重表現(xiàn)較好。該研究結(jié)果為陸稻農(nóng)家種資源的保護和利用提供了參考。

關(guān)鍵詞：陸稻農(nóng)家種;SNP;表型;遺傳多樣性;綜合評價中圖分類號：S511.6??????文獻(xiàn)標(biāo)識碼：A

Comprehensive Genetic Diversity Analysis of Eighty Late Upland Rice Landrace

Abstract： The comprehensive genetic diversity of 80 upland rice landrace?was evaluated based on SNP genotyping and phenotype identification， aiming to protect and utilize the upland rice landrace. The results showed that a total of 739 genotyping SNP data of the 80 landrace?was obtained using a 1K rice SNP array. Based on the SNP genotyping data， the upland rice landrace?was divided into Yunnan group， Hainan group and Guizhou group， through population structure analysis， PCA analysis and phylogenetic analysis. And there were 63， 11 and 6 landrace in the group， respectively. The genetic distance of the three groups was 0.83， 0.63 and 0.80， and the genetic distance between Hainan group and the other two groups was far relatively. The upland landrace， collected from Guizhou， were clustered into Guizhou group and Yunnan group respectively， indicating that there may be different sources of upland landrace in southern Guizhou. The coefficient variation of?11 agronomic and yield related traits ranged from 4.09% to 38.77%， and the?coefficient variation of panicle weight and effective tiller was larger than 25.00%， while the coefficient variation of plant height， grain length， grain width， grain length/width ratio and mean grain roundness was smaller than 10%. The traits diversity index ranged from 1.76 to 2.94， and seed length had the greatest genetic diversity among the traits. Correlation analysis revealed complex correlation ship between the traits， and the correlation coefficient between grain length/width ration and mean grain roundness was 0.98. The phenotype PCA analysis showed that the first five principal components?explained 29.08%， 23.42%， 12.35%， 9.77%?and 7.75% phenotypic variation respectively. And the first five principal components totally explained 82.37% of the phenotypic variation， reflecting most information of the 11 traits. The comprehensive score?of the landrace， extremely significant correlated with plant height， panicle length， flag leaf length， flag leaf width and effective tiller， was ranged from 0.15 to 0.73. The landrace ZRG36， ZRG101， ZRG58， ZRG31， ZRG1， ZRG102， ZRG104， ZRG63， ZRG29 and ZRG11 had the top comprehensive score all from higher elevations，. And the landrace， two from center region of Hainan Province， five from southern of Guizhou Province and three from Southwest and southeast of Yunnan Province， performed well in terms of panicle length， grain length， grain width， panicle weight and 1000-grain weight. The landrace were potential parent material for upland rice breeding or intermediate material for rice drought resistance breeding. The result will provide references for protecting and utilizing the upland rice landrace.AF2D856E-8F0A-401B-9F4A-53406B05FF5A

Keywords： upland rice landrace; SNP; phenotype; genetic variation; comprehensive evaluation

DOI： 10.3969/j.issn.1000-2561.2022.05.007

陸稻，也稱為旱稻，是栽培稻（Oryza sativa）的旱生類型，具有耐旱、耐貧瘠和抗逆性好的特點，是抗旱、氮高效栽培稻遺傳改良的寶貴資源。在全球氣溫升高、極端天氣頻繁和水資源短缺的背景下，我國稻作產(chǎn)業(yè)面臨的環(huán)境風(fēng)險不斷提高，因干旱或季節(jié)性干旱造成的水稻減產(chǎn)面積和范圍不斷擴大^[1]。因此，培育產(chǎn)量高、品質(zhì)優(yōu)、節(jié)水抗旱的水稻品種對于我國水稻產(chǎn)業(yè)發(fā)展具有重要意義^[2]。云南、貴州、廣西、海南等地是我國陸稻的傳統(tǒng)種植區(qū)，分布有豐富的陸稻資源，這些資源是水稻抗旱、抗逆等功能基因的重要來源^[3-5]。然而，隨著高產(chǎn)水稻品種的培育和應(yīng)用，傳統(tǒng)的農(nóng)家陸稻品系由于其產(chǎn)量水平較低導(dǎo)致種植面積急劇下降，陸稻農(nóng)家種流失嚴(yán)重^{[1， 6]}。此外，陸稻農(nóng)家種普遍具有植株較高、生育期較長等特點，也限制對這些資源的利用^[6]。表型是作物基因型和環(huán)境相互作用的綜合表現(xiàn)，可直觀反映作物的形態(tài)特征^{[3， 7]};遺傳多樣性可反映不同來源種質(zhì)資源的遺傳背景差異及變異水平^[8-9]。對作物種質(zhì)資源開展表型和遺傳多樣性評價，可以更全面地了解資源的特性，為育種利用提供重要參考依據(jù)。本課題組前期收集陸稻資源時，發(fā)現(xiàn)有80份陸稻晚稻農(nóng)家種資源?；诖?，本研究擬利用SNP標(biāo)記對從云南、貴州和海南等地收集的80份陸稻晚稻農(nóng)家種資源的基因型進(jìn)行檢測，結(jié)合農(nóng)藝性狀和產(chǎn)量性狀的鑒定，分析陸稻晚稻農(nóng)家種資源的遺傳多樣性特點，明確陸稻晚稻農(nóng)家種遺傳的特性和親緣關(guān)系，為陸稻晚稻農(nóng)家種的保護和利用提供參考。

1??材料與方法

1.1材料

以本課題組收集的來源于云南、貴州、海南和廣東的80份陸稻晚稻農(nóng)家種資源為供試材料（表1），上述材料在廣東湛江只能在6—11月種植（晚稻），春季種植均不能完成整個生育期。云南的陸稻晚稻資源主要來源于海拔1000～1400?m的滇東南（文山）和滇西南（臨滄、西雙版納和普洱），貴州的陸稻晚稻資源來源于平均海拔1000?m的黔南地區(qū)（平壩縣、冊亨縣和安龍縣），海南的陸稻晚稻資源來源于平均海拔1000?m的海南中部地區(qū)（瓊中縣和白沙縣）。

1.2 方法

1.2.1 ?SNP檢測??使用博瑞迪生物技術(shù)有限公司的水稻1KSNP芯片對80份資源進(jìn)行檢測。

1.2.2 ?表型鑒定??供試材料于2019—2020年連續(xù)兩年在廣東省湛江市中國熱帶農(nóng)業(yè)科學(xué)院湛江實驗站基地種植，每份材料種植3個小區(qū)，每小區(qū)

種5行，每行種10株，株距13.3 cm，行距25 cm，收獲前調(diào)查各材料的株高、穗長、劍葉長、劍葉寬、有效分蘗數(shù)共5個性狀，收獲后調(diào)查穗重、千粒種、種子長、種子寬、種子長寬比、種子圓度共6個性狀。所有性狀調(diào)查按照《水稻種質(zhì)資源描述規(guī)范和數(shù)據(jù)標(biāo)準(zhǔn)》執(zhí)行^[10]。

1.3數(shù)據(jù)處理

1.3.1 ?SNP數(shù)據(jù)分析??使用Plink（http：//www.?cog-genomics.org/plink2/）軟件對獲得的SNP芯片clean數(shù)據(jù)進(jìn)行質(zhì)控，去除缺失≥0.05、最小等位基因頻率（MAF）<0.05和=1的SNP位點。

1.3.2 ?群體結(jié)構(gòu)分析??使用fastSTRUCTURE軟件對群體結(jié)構(gòu)進(jìn)行分析，確定最佳分群^[11];使用R軟件包進(jìn)行PCA分析，并繪制群體結(jié)構(gòu)圖和PCA圖;使用MEGA 7.0軟件構(gòu)建進(jìn)化樹，并計算遺傳距離（P-distance）^[12]，同時使用iTol工具（http：//itol2.embl.de/）對進(jìn)化樹進(jìn)行美化。

1.3.3 ?多樣性分析??參照文獻(xiàn)[3， 13-14]中的方法對性狀劃分等級，即M±kσ（其中k=0、0.5、1.0、1.5、2.0，M為平均值，σ為標(biāo)準(zhǔn)差），將這些性狀表型值劃分10個等級，每1級中材料所占的比例（P_i）用于計算多樣性指數(shù)，遺傳多樣性指數(shù)即Shannon-Wiener index（H′）信息指數(shù)，其計算公式為：H′=-ΣP_i×lnP_i，式中P_i為某性狀第i級中包含的材料份數(shù)所占材料總份數(shù)的百分比。

1.3.4 ?綜合評價??使用模糊隸屬函數(shù)法將80份陸稻晚稻農(nóng)家種的11個表型性狀函數(shù)值進(jìn)行標(biāo)準(zhǔn)化^[13];參照綜合評價方法^[15]，使用R軟件包對11個表型性狀進(jìn)行PCA分析，計算主成分因子得分系數(shù)矩陣，結(jié)合權(quán)重系數(shù)來計算綜合評價F值。

2??結(jié)果與分析

2.1?陸稻晚稻農(nóng)家種資源的群體結(jié)構(gòu)分析

對80份陸稻晚稻農(nóng)家種資源進(jìn)行水稻1KSNP芯片檢測后，共獲得了1107個SNP位點，按照缺失和MAF進(jìn)行過濾后，共獲得739個符合要求的SNP位點。SNP的主基因頻率范圍為0.55～0.99，雜合度范圍為0～0.16，PIC值范圍為0～0.37。使用fastSTRUCTURE軟件對群體結(jié)構(gòu)進(jìn)行分析和最佳分組篩選，在K=3時，群體有最佳的分組（圖1A）。根據(jù)各品種的Q值，按照Q>0.5確定各品種的亞群，研究發(fā)現(xiàn)這些資源的遺傳背景比較清晰，不同亞群間的基因交流較少。按照地域可將這些資源分為云南亞群、海南亞群和貴州亞群。其中云南亞群共有63份資源，主要以云南的資源為主，包含10份來源于貴州和3份來源于廣東的陸稻晚稻農(nóng)家種資源;海南亞群包含11份海南的陸稻晚稻農(nóng)家種資源;貴州亞群包含6份來自于貴州的陸稻晚稻農(nóng)家種資源。分屬于不同亞群的貴州陸稻晚稻農(nóng)家種資源，其Q值差異較大。

2.2 陸稻晚稻農(nóng)家種資源的PCA分析和系統(tǒng)進(jìn)化分析

PCA分析結(jié)果表明，80份陸稻晚稻農(nóng)家種資源可分為云南亞群（Group1）、海南亞群（Group2）和貴州亞群（Group3）3個亞群（圖1B）。

進(jìn)化樹分析結(jié)果表明，80份陸稻晚稻農(nóng)家種資源可聚為云南亞群（Group1）、海南亞群（Group2）和貴州亞群（Group3）3類，與PCA和群體結(jié)構(gòu)的結(jié)果相一致，表明K=3為該陸稻晚稻資源的最佳分群數(shù)（圖1C）。3個亞群間的平均遺傳距離分別為0.83、0.63、0.80，表明云南亞群和貴州亞群的遺傳距離較近，海南亞群與其他兩個亞群的遺傳距離均較遠(yuǎn)。從圖1可看出，3個亞群處于同一水平，可進(jìn)一步將Group1所代表的云南亞群細(xì)分為3個亞群，其中Group1A為7份從貴州收集的陸稻晚稻農(nóng)家種資源，這些資源與Group3亞群中來自于同一地區(qū)的資源遺傳距離較遠(yuǎn)，而與來源于云南地區(qū)的資源遺傳距離較近，表明黔南地區(qū)的陸稻晚稻資源可能存在著不同的來源;Gruop1B和Gruop1C分別包含20份和36份陸稻晚稻農(nóng)家種資源，主要為云南西南地區(qū)和東南地區(qū)收集的陸稻晚稻農(nóng)家種資源，這2個亞群中均包含來自廣東的資源，Gruop1B中的2份廣東資源與云南西雙版納地區(qū)的2份資源聚在同一個分支，Gruop1C中的1份廣東資源與云

南臨滄的資源聚在同一個分支;3個細(xì)分亞群中均包含滇西南和滇東南地區(qū)的資源，西雙版納、臨滄、普洱和文山的資源均不能按地區(qū)分開，細(xì)分亞群兩兩的遺傳距離為0.30、0.28和0.22，表明云南亞群細(xì)分亞群間的遺傳距離很近，其中從貴州收集的Group1A與Gruop1B、Gruop1C的遺傳距離相對較遠(yuǎn)。

2.3陸稻晚稻農(nóng)家種資源的表型變異和多樣性分析

對80份陸稻晚稻農(nóng)家種資源連續(xù)2年進(jìn)行表型鑒定，結(jié)果表明，11個性狀的表型分布和變異差異明顯（圖2和表2）。11個性狀的變異系數(shù)為4.09%（種子長）～38.77%（穗重），其中有效分蘗數(shù)和穗重的變異系數(shù)較大，均超過25.00%;株高、種子長、種子寬、種子長寬比和種子圓度均值變異系數(shù)較小，均低于10%。11個性狀的遺傳多樣性指數(shù)范圍為1.76（劍葉寬）～2.94（種子長），遺傳多樣性豐富，其中種子長（2.94）和株高（2.85）的多樣性指數(shù)最高。

2.4 陸稻晚稻農(nóng)家種資源的表型相關(guān)分析

對陸稻晚稻農(nóng)家種11個性狀進(jìn)行相關(guān)分析，結(jié)果見表3，株高與穗長、劍葉長和劍葉寬呈極顯著正相關(guān)，穗長與劍葉長和劍葉寬呈極顯著正相關(guān)。種子長和種子長寬比呈極顯著正相關(guān)，與圓度均值呈極顯著負(fù)相關(guān);種子寬與種子長寬比呈極顯著負(fù)相關(guān)，與圓度均值和穗重呈極顯著正相關(guān);種子長寬比與圓度均值和穗重呈極顯著負(fù)相關(guān)。株高與有效分蘗數(shù)、劍葉長和劍葉寬、種子長和千粒重均呈顯著正相關(guān)，穗長與千粒重呈顯著負(fù)相關(guān)。其中種子長寬比與圓度均值的相關(guān)系數(shù)最高為0.98。

2.5 陸稻晚稻農(nóng)家種資源的表型綜合評價

對11個性狀進(jìn)行PCA分析，前5個主成分的累計貢獻(xiàn)率達(dá)82.37%，可以解釋表型性狀的絕大部分信息（表4）。第一主成分的特征值最大，為3.19，貢獻(xiàn)率為29.08%，其中種子長、寬、長

寬比和圓度均值的特征向量絕對值較大，主要反映產(chǎn)量相關(guān)信息;第二主成分的貢獻(xiàn)率為23.42%，其中株高、穗長、劍葉長和劍葉寬的特征向量較大，說明該成分主要反映株型相關(guān)的信息;第三主成分的貢獻(xiàn)率為12.35%，其中種子長、寬、千粒重、穗重的特征值較大，表明該特征值主要反映產(chǎn)量相關(guān)信息;第四主成分中有效分蘗數(shù)特征值最大為0.82，有效分蘗既與株型相關(guān)，又與產(chǎn)量相關(guān)，說明該成分主要反映株型和產(chǎn)量相關(guān)性狀;第五主成分中千粒重的特征值最大，其次為穗重和種子長，表明該成分主要反映產(chǎn)量相關(guān)性狀。

根據(jù)PCA分析獲得了5個因子得分公式：

F₁=0.04X₁+0.14X₂+0.19X₃+0.06X₄+0.11X₅–0.34X₆+?0.39X₇–0.53X₈+0.54X₉+0.27X₁₀–0.10X₁₁AF2D856E-8F0A-401B-9F4A-53406B05FF5A

F₂=0.52X₁+0.45X₂+0.45X₃+0.50X₄+0.05X₅+0.05X₆–0.11X₇+0.12X₈–0.10X₉–0.13X₁₀–0.14X₁₁

F₃=0.11X₁+0.02X₂+0.12X₃–0.01X₄+0.28X₅+0.55X₆+0.42X₇+0.08X₈–0.05X₉+0.34X₁₀+0.55X₁₁

F₄=0.17X₁+0.04X₂–0.22X₃–0.18X₄+0.82X₅–0.18X₆–0.29X₇+0.10X₈–0.06X₉+0.24X₁₀–0.18X₁₁

F₅=0.16X₁–0.14X₂+0.24X₃–0.14X₄+0.09X₅–0.41X₆–0.27X₇–0.08X₈+0.10X₉–0.33X₁₀+0.71X₁₁

將F₁～F₅的得分，根據(jù)各成分的權(quán)重（0.35、0.28、0.15、0.12、0.09），計算出各資源的綜合得分F，80份資源的得分范圍為0.15（ZRG7）～0.73（ZRG36），平均值為0.44（表5）。F值與11個性狀的相關(guān)分析表明，F值與株高、穗長、劍葉長、劍葉寬和有效分蘗數(shù)呈極顯著正相關(guān)，而與

其他性狀的相關(guān)關(guān)系不顯著（表6）。綜合排名前10位的農(nóng)家種為ZRG36、ZRG101、ZRG58、ZRG31、ZRG1、ZRG102、ZRG104、ZRG63、ZRG29、ZRG11，這些資源的穗長、種子長、種子寬、穗重和千粒重表型較好，且均分布于海拔較高地區(qū)，其中ZRG31和ZRG36來源于海南中部地區(qū)，ZRG1、ZRG63、ZRG101、ZRG102 和ZRG104來源于貴州黔南，ZRG11、ZRG29和ZRG58來源于滇西南和東南地區(qū)，是潛在的陸稻晚稻育種親本材料或水稻抗旱育種的中間材料。

3 ?討論

我國是栽培稻的起源地之一，稻作歷史悠久，稻種資源豐富。陸稻是栽培稻的旱作生態(tài)類型，是在自然環(huán)境與人工馴化共同作用下形成的一類具有耐旱、耐貧瘠等特性的特殊種質(zhì)資源，是水稻育種中最直接、最有效的抗旱基因供體來源^[16-18]。云南、海南和貴州具有豐富的陸稻資源^{[3， 19]}，本研究利用739個SNP標(biāo)記將從這些地區(qū)收集的80份陸稻晚稻資源按照來源劃分為云南亞群、海南亞群和貴州亞群，其中海南亞群與其他2個亞群的親緣關(guān)系均較遠(yuǎn)，可能與海南島與大陸隔離，與其他地區(qū)的資源交流較少有關(guān);云南亞群中包含部分貴州收集的晚稻資源，這些晚稻資源與云南的晚稻資源遺傳距離較近，而與同樣來源于貴州黔南的貴州亞群的晚稻資源遺傳距離較遠(yuǎn)，推測云南和貴州地區(qū)可能發(fā)生了晚稻農(nóng)家種的交流，且貴州地區(qū)的晚稻農(nóng)家種可能有不同的來源和馴化;此外，在云南亞群中，還存在3份來源于廣東的晚稻農(nóng)家種資源，分別與云南臨滄和西雙版納的晚稻資源遺傳距離較近，而與海南地區(qū)的晚稻資源遺傳距離較遠(yuǎn)，由于這3份資源（從國家水稻中期庫引進(jìn)）具體來源地不詳，且收集的資源中廣東地區(qū)的晚稻資源較少，初步推測可能與晚稻農(nóng)家種的引進(jìn)交流有關(guān)，還需進(jìn)一步結(jié)合本地區(qū)其他陸稻晚稻資源進(jìn)行綜合分析。前人研究發(fā)現(xiàn)，云南滇西南與滇東南的的陸稻品種之間的差異顯著^[6]，而本研中鑒定的滇東南和滇西南的陸稻晚稻資源差異并不明顯，可能與所鑒定的晚稻資源類型較為單一有關(guān)。本研究結(jié)果還表明，80份陸稻晚稻資源的11個表型性狀的表型分布和變異差異明顯，遺傳多樣性指數(shù)范圍為1.76～2.85，11個性狀間形成了復(fù)雜的相關(guān)關(guān)系，80個品種的得分范圍為0.15～0.73，這些結(jié)果對陸稻晚稻農(nóng)家種資源的利用和保護具有參考價值。AF2D856E-8F0A-401B-9F4A-53406B05FF5A

遺傳多樣性評價是認(rèn)識作物基因型和表型多樣性程度，以及挖掘和利用資源中優(yōu)異基因資源的理論和實踐基礎(chǔ)^[20-21]。表型多樣性有助于從整體上認(rèn)識作物性狀的多樣性程度，為優(yōu)異資源的挖掘和利用提供理論和實踐依據(jù)。本課題組在前期鑒定中發(fā)現(xiàn)，收集的陸稻資源中，有80份來源于云南、貴州、廣東和海南的農(nóng)家種晚稻資源，這些資源只有晚稻種植才能完成整個生育期。本研究對這80份陸稻晚稻資源進(jìn)行表型鑒定，發(fā)現(xiàn)所調(diào)查的11個性狀表現(xiàn)出豐富的表型變異。從地理分布來看，在海南和貴州收集的陸稻晚稻資源均分布在海拔約為1000 m的地區(qū);在云南滇西南和滇東南地區(qū)收集的陸稻晚稻資源，主要分布在海拔1000～1400 m的半山區(qū)^[6]。從氣候條件看，海南中部地區(qū)屬于熱帶季風(fēng)氣候，日照充足，陸稻晚稻種植期間平均溫度在26℃左右;黔南地區(qū)屬于亞熱帶濕潤季風(fēng)氣候，陰雨天多，日照少;云南滇南陸稻資源分布地區(qū)氣候與貴州類似，為半山云霧多濕區(qū)，日照不充分，陸稻晚稻生長期平均溫度在23℃左右。相似的地理和生態(tài)環(huán)境促進(jìn)了云南和貴州地區(qū)的資源引種和交流，可能是造成貴州黔南地區(qū)陸稻晚稻資源遺傳差異較大的原因。這些地理和氣候環(huán)境因素，是陸稻農(nóng)家種馴化的重要動力，造成了陸稻資源豐富的表型多樣性，也從生態(tài)方面解釋了為何海南的資源與貴州、云南的資源遺傳距離較遠(yuǎn)。

分子標(biāo)記是開展遺傳多樣性研究的重要工具，在陸稻研究中，RFLP、RAPD、SSR、SRAP等多種分子標(biāo)記已經(jīng)用于分析資源的遺傳背景、親緣關(guān)系和基因多態(tài)性^{[17-18， 22-23]}，如王英等^[22]利用22條多態(tài)性RAPD引物分析了49分陸稻的遺傳基礎(chǔ)，并構(gòu)建了指紋圖譜;張立娜等^[17]利用39對SSR引物分析了158份不同來源的陸稻資源的遺傳相似性，發(fā)現(xiàn)秈型旱稻地方品種的基因多樣性顯著高于粳型旱稻地方品種;徐建欣等^[6]利用24對SSR引物分析了131份云南陸稻品種的遺傳多樣性，引物的等位基因數(shù)平均為8.125個，基因多樣性指數(shù)平均為0.6543，研究發(fā)現(xiàn)滇西南和滇南地區(qū)存在豐富的遺傳變異。而本研究利用739個SNP分析了80份陸稻晚稻農(nóng)家種資源的親緣關(guān)系，經(jīng)群體結(jié)構(gòu)分析、PCA分析和聚類分析，均將80份陸稻晚稻農(nóng)家種資源分成了云南、海南和貴州3個亞群，說明了種群劃分的準(zhǔn)確性。SNP標(biāo)記與傳統(tǒng)的SSR和RAPD等標(biāo)記相比，由于只涉及單個堿基的變異，其檢測到的等位基因數(shù)一般是2～3個，但基于大量水稻資源測序構(gòu)建的水稻SNP芯片，能夠更均勻的檢測基因組上的變異，可以更加準(zhǔn)確對資源遺傳背景和親緣關(guān)系進(jìn)行判定。隨著芯片技術(shù)的成熟和測序成本的降低，基于基因芯片的SNP標(biāo)記展現(xiàn)出高通量、高效、快速的特點，已經(jīng)越來越多用于水稻^[24]、玉米^[25]、棉花^[26]、油菜^[27]等多種作物的遺傳多樣性研究和指紋圖譜構(gòu)建，將更加高效的用于種質(zhì)資源評價和品種鑒定，促進(jìn)優(yōu)異資源的挖掘和利用。

參考文獻(xiàn)

[1]?殷會德， 米鐵柱， 劉佳音， 于??萌， 顧曉振， 郭海鵬， 李繼明. 早稻研究現(xiàn)狀及發(fā)展前景[J]. 中國稻米， 2020（4）： 22-24.YIN H D， MI T Z， LIU J Y， YU M， GU X Z， GUO H P， LI J M. Research status and development prospects of upland rice[J]. China Rice， 2020（4）： 22-24. （in Chinese）

[2]?羅利軍， 梅捍衛(wèi)， 余新橋， 劉鴻艷， 馮芳君. 節(jié)水抗旱稻及其發(fā)展策略[J]. 科學(xué)通報， 2011， 56（11）： 804-811.LUO L J， MEI H W， YU X Q， LIU H Y， FENG F J. Water-saving and drought-resistance rice and its development strategy[J].?Chinese Science Bulletin， 2011， 56（11）： 804-811. （in Chinese）

[3]?宮彥龍， 雷??月， 閆志強， 劉雪薇， 張大雙， 吳健強， 朱速松. 不同生態(tài)區(qū)粳稻資源表型遺傳多樣性綜合評價[J]. 作物雜志， 2020（5）： 71-79.GONG Y L， LEI Y， YAN Z Q， LIU X W， ZHANG D S， WU J Q， ZHU S S. Comprehensive evaluation of phenotype genetic diversity in japonica rice germplasm resources in different ecological zones[J]. Crops， 2020（5）： 71-79. （in Chinese）

[4]?王一平. 旱稻種質(zhì)資源遺傳多樣性研究[D]. 北京： 中國農(nóng)業(yè)科學(xué)院， 2006.WANG Y P. Study on genetic diversity of upland rice resources[D]. Beijing： Chinese Academy of Agricultural Sciences， 2006. （in Chinese）

[5]?張立娜. 中國水旱稻地方品種遺傳多樣性分析[D]. 北京： 中國農(nóng)業(yè)科學(xué)院， 2010.ZHANG L N. Genetic diversity of lowland and upland rice resources in China[D]. Beijing： Chinese Academy of Agricultural Sciences， 2010. （in Chinese）AF2D856E-8F0A-401B-9F4A-53406B05FF5A

[6]?徐建欣， 王云月， 姚??春， 劉云霞， 湯??淼， 陸明興. 利用SSR分子標(biāo)記分析云南陸稻品種遺傳多樣性[J]. 中國水稻科學(xué)， 2012， 26（2）： 155-164.XU J X， WANG Y Y， YAO C， LIU Y X， TANG M， LU M X. Genetic diversity of upland rice varieties in Yunnan province revealed by SSR markers[J]. Chinese Journal of Rice Science， 2012， 26（2）： 155-164. （in Chinese）

[7]?都真真， 李錫香， 宋江萍， 武亞紅， 趙??青， 徐??婷， 張曉輝， HELLIER B， HU J G， 王海平.?228份引進(jìn)大蒜資源的表型多樣性分析及適應(yīng)性初步評價[J]. 植物遺傳資源學(xué)報， 2019， 20（5）： 1186-1196.DU Z Z， LI X X， SONG J P， WU Y H， ZHAO Q， XU T， ZHANG X H， HELLIER B， HU J G， WANG H P. Phenotypic diversity and adaptability analysis of 228 accessions of introduced garlic genetic resources[J]. Journal of Plant Genetic Resources， 2019， 20（5）： 1186-1196. （in Chinese）

[8]?萬映伶， 劉愛青， 張孔英， 劉??燕. 菏澤和洛陽芍藥品種資源表型多樣性研究[J]. 北京林業(yè)大學(xué)學(xué)報， 2018， 40（3）： 110-121.WAN Y L， LIU A Q， ZANG K Y， LIU Y. Phenotype diversity of herbaceous peony variety resources in Heze， Shandong of eastern China and Luoyang， Henan of central China[J]. Journal of Bjing Forestry University， 2018， 40（3）： 110-121. （in Chinese）

[9]?陳??越， 張敦宇， 丁明亮， 王玲仙， 肖素勤， 柯 ?學(xué)， 程在全. 多個省份水稻資源的表型多樣性與優(yōu)異資源的篩選[J]. 浙江農(nóng)業(yè)學(xué)報， 2019， 31（11）： 1779-1789.CHEN Y， ZHANG D Y， DING M L， WANG L X， XIAO S Q， KE X， CHENH Z Q. Phenotypic diversity of rice resources in multiple provinces and screening of excellent resources[J]. Acta Agriculturae Zhejiangensis， 2019， 31（11）： 1779-1789. （in Chinese）

[10]?韓龍植. 水稻種質(zhì)資源描述規(guī)范和數(shù)據(jù)標(biāo)準(zhǔn)[M]. 北京： 中國農(nóng)業(yè)出版社， 2006.HAN L Z. Descriptors and data standard for rice （Oryzaspp.）[M]. Beijing： China Agriculture Press， 2006. （in Chinese）

[11]?RAJ A， STEPHENS M， ?PRITCHARD J K. fastSTRUCTURE： variational inference of population structure in large SNP data sets[J]. Genetics， 2014， 197（2）： 573-589.

[12]?KUMAR S， STECHER G， TAMURA K. MEGA7： Molecular evolutionary genetics analysis version 7.0 for bigger datasets[J]. Molecular Biology and Evolution， 2016， 33（7）： 1870-1874.

[13]?徐志軍， 吳小麗， 胡小文， 劉??洋. 33份引進(jìn)花生資源表型遺傳多樣性分析及在粵西地區(qū)的適應(yīng)性初步評價[J]. 熱帶作物學(xué)報， 2021， 42（7）： 1885-1895.XU Z J， WU X L， HU X W， LIU Y. Phenotypic diversity of 33 introduced peanut germplasm accessions and preliminary adaptability evaluation in West Guangdong[J]. Chinese Journal of Tropical Crops， 2021， 42（7）： 1885-1895. （in Chinese）

[14]?孫珍珠， 李秋月， 王小柯， 趙婉彤， 薛??楊， 馮錦英， 劉小豐， 劉夢雨， 江??東. 寬皮柑橘種質(zhì)資源表型多樣性分析及綜合評價[J]. 中國農(nóng)業(yè)科學(xué)， 2017， 50（22）： 4362-4372.SUN Z Z， LI Q Y， WANG X K， ZHAO W T， XUE Y， FENG J Y， LIU X F， LIU M Y， JIANG D. Comprehensive evaluation and phenotypic diversity analysis of germplasm resources in mandarin[J]. Scientia Agricultura Sinica， 2017， 50（22）： 4362-4383.?（in Chinese）AF2D856E-8F0A-401B-9F4A-53406B05FF5A

[15]?胡標(biāo)林， 萬??勇， 李??霞， 雷建國， 羅向東， 嚴(yán)文貴， 謝建坤. 水稻核心種質(zhì)表型性狀遺傳多樣性分析及綜合評價[J]. 作物學(xué)報， 2012， 38（5）： 829-839.HU B L， WAN Y， LI X， LEI J G， LUO X D， YAN W G， XIE J K. Analysis on genetic diversity of phenotypic traits in rice （Oryza sativa） core collection and its comprehensive assessment[J]. Acta Agronomica Sinica， 2012， 38（5）： 829-839. （in Chinese）

[16]?王一平， 魏興華， 華??蕾， 袁筱萍， 余漢勇， 徐??群， 湯圣祥. 不同地理來源旱稻種質(zhì)資源的遺傳多樣性分析[J]. 作物學(xué)報， 2007， 33（12）： 2034-2040.WANG Y P， WEI X H， HUA L， YUAN X P， YU H Y， XU Q， TANG S X. Genetic diversity in upland rice germplasm from different geographic regions[J]. Acta Agronomica Sinica， 2007， 33（12）： 2034-2040. （in Chinese）

[17]?張立娜， 曹桂蘭， 韓龍植. 中國不同地理來源旱稻地方品種的遺傳相似性研究[J]. 中國農(nóng)業(yè)科學(xué)， 2010， 43（17）： 3481-3488.ZHANG L N， CAO G L， HAN L Z. Analysis of genetic similarity for upland landrace rices from different geographical origins in China[J]. Scientia Agricultura Sinica， 2010， 43（17）： 3481-3488. （in Chinese）

[18]?楊學(xué)輝， 袁??潔， 陳小均， 阮仁超， 何海永， 吳石平， 王云月. 貴州旱稻種質(zhì)資源的SSR遺傳多樣性分析[J]. 分子植物育種， 2009， 7（5）： 890-896.YANG X H， YUAN J， CHEN X J， RUAN R C， HE H Y， WU S P， WANG Y Y. Genetic diversity of upland rice germplasm resources from Guizhou based on SSR markers[J]. Molecular Plant Breeding， 2009， 7（5）： 890-896. （in Chinese）

[19]?游俊梅， 陳惠查， 金桃葉， 阮仁超. 貴州地方旱稻種質(zhì)資源遺傳多樣性評價[J]. 種子， 2005（4）： 79-81， 84.YOU J M， CHEN H C， JIN T Y， RUAN R C. Genetic diversity evaluation of upland rice germplasm resources in Guizhou province[J]. Seeds， 2005（4）： 79-81， 84. （in Chinese）

[20]?陳??巍， 王力榮， 朱更瑞， 方偉超， 陳昌文， 曹??珂， 張紹鈴. 基于SSR標(biāo)記和生物學(xué)性狀進(jìn)行桃遺傳多樣性的比較分析[J].?植物遺傳資源學(xué)報， 2009， 10（1）： 86-90， 96.CHEN W， WANG L R， ZHU G R， FANG W C， CHEN C W， CAO K， ZHANG S L. Comparison of genetic diversity among peach cultivars based on biological traits and SSR markers[J]. Journal of Plant Genetic Resources， 2009， 10（1）： 86-90， 96. （in Chinese）

[21]?CAMPBELL D R. Using phenotypic manipulations to study multivariate selection of floral trait associations[J]. Annals of Botany， 2009， 103（9）： 1557-1566.

[22]?王??英， 葉??通， 邱海燕， 莊南生. 49份旱稻種質(zhì)RAPD標(biāo)記遺傳多樣性分析[J]. 中國農(nóng)學(xué)通報， 2011， ?27（24）： 21-28.WANG Y， YE T， QIU H Y， ZHUANG N S. Analysis on genetic diversity of 49 upland rice germplasm using RAPD markers[J]. Chinese Agricultural Science Bulletin， 2011， 27（24）： 21-28. （in Chinese）

[23]?姚宇飛， 王??英， 莊南生， 高和瓊. 利用SRAP標(biāo)記分析海南旱稻地方品種的遺傳多樣性及其分子身份證的構(gòu)建[J]. 生物技術(shù)通報， 2014（11）： 97-106. （in Chinese）YAO Y F ， WANG Y， ZHUANG N S， GAO H Q. Analysis on genetic diversity of hainan upland rice local varietiesand establishment of molecular ID with SRAP markers[J]. Biotechnology Bulletin， 2014（11）： 97-106. （in Chinese）AF2D856E-8F0A-401B-9F4A-53406B05FF5A

[24]?李梓榕， 袁??雄， 陳??葉， 鄭興飛， 胡中立， 李蘭芝. 基于全基因組SNP高效鑒定水稻種質(zhì)資源并構(gòu)建指紋圖譜[J]. 分子植物育種， 2020， 18（18）： 6050-6057.LI Z R， YUAN X， CHEN Y， ZHENG X F， HU Z L， LI L Z. Effective identification for varieties by genome-wide SNPs and establishment of fingerprint for rice germplasms[J]. Molecular Plant Breeding， 2020， 18（18）： 6050-6057. （in Chinese）

[25]?田紅麗， 楊??揚， 王??璐， 王??蕊， 易紅梅， 許理文， 張云龍， 葛建镕， 王鳳格， 趙久然. 兼容型maizeSNP384標(biāo)記篩選與玉米雜交種DNA指紋圖譜構(gòu)建[J]. 作物學(xué)報， 2020， 46（7）： 1006-1015.TIAN H L， YANG Y， WANG L， WANG R， YI H M， XU L W， ZHANG Y L， GE J R， WANG F G， ZHAO J R. Screening of compatible maizeSNP384 markers and the construction of DNA fingerprints of maize varieties[J]. Acta Agronomica Sinica， 2020， 46（7）： 1006-1015. （in Chinese）

[26]?孫正文， 匡??猛， 馬峙英， 王省芬. 利用CottonSNP63K芯片構(gòu)建棉花品種的指紋圖譜[J]. 中國農(nóng)業(yè)科學(xué)， 2017， 50（24）： 4692-4704.SUN Z W， KUANG M， MA Z Y， WANG X F. Construction of cotton variety fingerprints using CottonSNP63K array[J]. Scientia Agricultura Sinica， 2017， 50（24）： 4692-4704. （in Chinese）

[27]?趙仁欣， 李森業(yè)， 郭瑞星， 曾新華， 文??靜， 馬朝芝， 沈金雄， 涂金星， 傅廷棟， 易??斌. 利用SNP芯片構(gòu)建我國冬油菜參試品種DNA指紋圖譜[J]. 作物學(xué)報， 2018， 44（7）： 956-965.ZHAO R X， LI S Y， GUO R X， ZENG X H， WEN J， MA C Z， SHEN J X， TU J X， FU T D， YI B. Construction of DNA fingerprinting forBrassica napusvarieties based on SNP chip[J]. Acta Agronomica Sinica， 2018， 44（7）： 956-965. （in Chinese）AF2D856E-8F0A-401B-9F4A-53406B05FF5A