基于SSR熒光標(biāo)記構(gòu)建板栗品種（系）核心種質(zhì)群

關(guān)鍵詞：板栗；品種；SSR；遺傳多樣性；核心種質(zhì)

中圖分類號：S664.2 文獻(xiàn)標(biāo)志碼：A 文章編號：1009-9980（2023）02-0230-12

隨著中國種業(yè)創(chuàng)新力度的不斷加大，人們對各類種質(zhì)資源的保護(hù)和利用愈發(fā)重視，但是植物資源相當(dāng)豐富，且植株普遍高大，全部保存成本會很高。因此，構(gòu)建植物資源核心種質(zhì)，是一種有效解決該問題的方法。核心種質(zhì)的概念最早由Frankel 等[1]提出，即采用一定的方法選取最少的種質(zhì)資源數(shù)量以最大程度地代表該類種質(zhì)資源的遺傳多樣性。核心種質(zhì)的構(gòu)建主要有兩種途徑，分別是基于表型性狀和基于分子標(biāo)記基因型數(shù)據(jù)的構(gòu)建方法?；诒硇托誀顦?gòu)建核心種質(zhì)極易受環(huán)境的影響，因此誤差較大，且需要對多種性狀進(jìn)行分析，工作量巨大。分子標(biāo)記技術(shù)以信息量大、效率高、不受環(huán)境影響等特點(diǎn)[2]，尤其是迅速發(fā)展的SSR標(biāo)記具有穩(wěn)定好、多態(tài)性高、共顯性遺傳、操作簡便等優(yōu)點(diǎn)[3]，已逐漸成為構(gòu)建核心種質(zhì)的主要方法。目前，SSR 分子標(biāo)記技術(shù)已經(jīng)成功應(yīng)用到多種作物的核心種質(zhì)構(gòu)建中，對資源保護(hù)和創(chuàng)新利用以及新品種的選育發(fā)揮了重要作用。如在蔬菜中，崔竣杰等[4]根據(jù)表型和SSR 標(biāo)記構(gòu)建了苦瓜的核心種質(zhì)；李金龍等[5]利用SSR 分子標(biāo)記技術(shù)研究甜蕎種質(zhì)資源的群體結(jié)構(gòu)，進(jìn)而構(gòu)建甜蕎的初級核心種質(zhì)庫；吳茵[6]在利用SRAP分子標(biāo)記構(gòu)建辣椒初級核心種質(zhì)的基礎(chǔ)上，利用SSR 分子標(biāo)記壓縮初級核心種質(zhì)進(jìn)而構(gòu)建了最終的核心種質(zhì)；果樹中魯敏等[7]采用位點(diǎn)優(yōu)先取樣策略對收集的102 份自然種質(zhì)構(gòu)建核心種質(zhì)；李冬波等[8]利用SSR分子標(biāo)記對廣西荔枝種質(zhì)資源進(jìn)行遺傳多樣性分析并構(gòu)建核心種質(zhì)，從而對資源保護(hù)發(fā)揮了重要作用；鄒梁峰[9]基于表型性狀和SSR 分子標(biāo)記構(gòu)建了毛花獼猴桃的核心種質(zhì)，從而對中國廣泛分布的這種特有的種質(zhì)資源進(jìn)行了很好的保護(hù)利用。劉娟等[10-11]分別基于表型性狀和ISSR 分子標(biāo)記構(gòu)建了野杏和南疆杏資源的核心種質(zhì)，花卉中陳明堃等[12]統(tǒng)計(jì)建蘭的表型并利用SSR 標(biāo)記篩選了核心種質(zhì)，能夠最大程度代表原始種質(zhì)的遺傳多樣性。此外，林木中的杉木[13]、毛白楊[14]、刺槐[15]、核桃[16]、白蠟[17]、楸樹[18]等都構(gòu)建了核心種質(zhì)，對于資源保護(hù)和創(chuàng)新利用有重要作用。

板栗（Castanea mollissima Bl.）為殼斗科栗屬植物，果實(shí)味道鮮美，具有很高的營養(yǎng)價(jià)值和藥用價(jià)值[19]，早在6000 年前，古人就已經(jīng)食用野生板栗[20]。中國板栗資源豐富，可分為華北、長江流域、西南、東南、西北與東北6 個品種群[21]，現(xiàn)階段約有300 個栽培品種（系）。同時，隨著自然變異和人工選擇，市場上逐漸出現(xiàn)新品種，同時淘汰一些過時的品種。因此，如果不及時發(fā)現(xiàn)和保護(hù)，就會造成優(yōu)異種質(zhì)資源的遺失和浪費(fèi)，從而對今后資源的創(chuàng)新利用造成困難。所以，構(gòu)建中國板栗品種的核心種質(zhì)，對其保護(hù)和利用具有重要的意義。前人對板栗品種核心種質(zhì)的構(gòu)建僅在某一地區(qū)或所采集的板栗品種資源不夠全面[22-23]，不足以代表全國板栗品種資源的遺傳多樣性。筆者以國家板栗資源圃和江蘇良種繁育基地收集的342 個品種（系）資源為材料，應(yīng)用熒光SSR 分子標(biāo)記技術(shù)獲取片段信息，基于分層取樣法、模擬退火算法和隨機(jī)搜索算法構(gòu)建全國板栗品種核心種質(zhì)，從而為今后板栗新品種的創(chuàng)制、品種的種質(zhì)鑒定和保存利用提供理論基礎(chǔ)。

1 材料和方法

1.1 試驗(yàn)材料

試驗(yàn)材料為從山東省泰安市國家板栗資源圃和位于江蘇板栗良種繁殖基地搜集的共342 個板栗品種，它們原生地分布在16 個省份。資源圃位于當(dāng)?shù)氐那鹆晟降氐貐^(qū)，降水量較大，土壤適合板栗生長和栽培。

1.2 試驗(yàn)方法

1.2.1 DNA提取 板栗葉片DNA提取采用美基生物公司生產(chǎn)的廣譜植物DNA 提取試劑盒。采用1.5%瓊脂糖凝膠電泳檢測DNA完整度，用超微量分光光度計(jì)測定DNA濃度和純度，將合格的DNA保存在-20 ℃冰箱中。

1.2.2 熒光毛細(xì)管電泳SSR-PCR 擴(kuò)增 SSR 引物采用本實(shí)驗(yàn)室之前所篩選得到的21 對高多態(tài)SSR標(biāo)記進(jìn)行PCR擴(kuò)增[22]。PCR擴(kuò)增采用20 μL反應(yīng)體系：其中包含1 μL模板DNA（20～50 ng· μL-1），10 μL2×Taq PCR Master Mixes（Takara， DaLian China），0.1 μL 正向引物（10 μmol · L- 1），0.3 μL 反向引物（10 μmol·L-1），0.2 μL 熒光標(biāo)記（FAM，HEX，ROX）的M13引物（10 μmol·L-1）（中國天津，擎科）和9.4 μLddH2O。反應(yīng)程序?yàn)椋?4 ℃預(yù)變性3 min，94 ℃變性30 s，56 ℃ 退火30 s，72 ℃延伸30 s，35 個循環(huán)，最后72 ℃ 延伸10 min，10 ℃ 保存（BIO- RAD PCRThermal Cycler T100， USA）。

將3 對熒光SSR引物的PCR擴(kuò)增產(chǎn)物合并在一起，然后使用ABI 3730XL DNA 測序儀（AppliedBiosystems， Foster City， CA， USA）通過毛細(xì)管電泳分析混合物中片段位置信息，并使用Gene Marker V2.2.0 軟件（Soft Genetics LLC， State College， PA，USA）讀取等位位點(diǎn)具體的片段大小信息。

1.2.3 數(shù)據(jù)分析 按照PowerMarker V3.25[24]軟件要求的格式，在Excel 表格中將讀取的SSR 原始“bp值型”數(shù)據(jù)轉(zhuǎn)換成powermarker 格式，缺失的數(shù)據(jù)賦值為“0”，并利用該軟件及GenAlEx6.51[25-26]計(jì)算遺傳多樣性指標(biāo)參數(shù)，包括等位基因數(shù)量（Na）、有效等位基因數(shù)（Ne）、觀測雜合度（Ho）、期望雜合度（He）、多態(tài)信息含量（PIC）。

1.2.4 核心種質(zhì)的構(gòu)建 （1）基于分層取樣法篩選核心種質(zhì)?；谧畲蟮任换驍?shù)原則，根據(jù)聚類圖采用逐步聚類留樣，每一輪聚類，隨機(jī)去除最低分類水平上兩個個體中的一個，剩下的一個入選提取比例，第二輪提取比例為從上次提取的樣品中聚類篩選，從而形成八個提取比例（100%、83.33%、66.67%、46.49%、31.58%、20.76%、13.74%、8.77%），并計(jì)算這8 個提取比例的等位位點(diǎn)數(shù)量。

（2）基于模擬退火算法篩選核心種質(zhì)。為構(gòu)建合適的核心種質(zhì)，再通過基于最大等位基因數(shù)（SAAN）的模擬退火算法來篩選和構(gòu)建中國板栗品種的核心種質(zhì)。在Power Marker v3.25 軟件中建立一個梯度樣本抓取數(shù)量方案來計(jì)算所選的不同容量的樣品所保留的等位位點(diǎn)的數(shù)量，從而確定最終篩選的核心種質(zhì)樣品數(shù)量。為保證盡可能的保留等位位點(diǎn)數(shù)，采用5，10，15，20，25，30，35，40，45，55，60，65，70，75，80，85，90，95，100，105，110，115，120，125，130，135，140，145，150，155，160，165，170，175，180，185，190，195，200，205，210，215，220，225，230，235，240，245，250，255，260，265，270，275，280，285，290，295，300，305，310，315，320，325，330，335，340的梯度抓取方案來計(jì)算不同取樣數(shù)量的等位位點(diǎn)數(shù)，從而確定等位位點(diǎn)最大時的抓取數(shù)量。

（3）基于隨機(jī)搜索算法篩選核心種質(zhì)?；谧畲蟮任换驍?shù)（RS， AN， AN）的隨機(jī)搜索算法來構(gòu)建板栗品種的核心種質(zhì)。該方法同樣利用上述基于模擬退火算法的梯度抓取方案，比較不同取樣數(shù)量的等位位點(diǎn)數(shù)量，從而確定抓取數(shù)量。

將上述3 種策略所獲得的取樣數(shù)量、等位位點(diǎn)數(shù)量等進(jìn)行比較以確定最終的取樣數(shù)量，進(jìn)而來評估出3 種核心種質(zhì)構(gòu)建策略中的最優(yōu)策略。

最后通過Prism 軟件[27]對核心種質(zhì)和原始種質(zhì)之間的Na、Ne、Ho、He 和PIC 值進(jìn)行T 檢驗(yàn)，從而驗(yàn)證所構(gòu)建的核心種質(zhì)的準(zhǔn)確性、科學(xué)性和可靠性。

2 結(jié)果與分析

2.1 原始板栗品種（系）種質(zhì)聚類分析

對所有資源進(jìn)行遺傳多樣性分析，結(jié)果表明342 份板栗品種資源共有212 個等位位點(diǎn)，平均每對引物所保留的等位基因數(shù)（Na）為10.095，有效等位基因數(shù)（Ne）為3.488，觀測雜合度為0.596，期望雜合度為0.658，多態(tài)信息含量為0.642。同時基于UPGMA方法進(jìn)行聚類分析，結(jié)果如圖1 所示。全部資源可以劃分為3 個組，藍(lán)色部分大多為北方品種，但混有31 份南方品種；紅色部分大多為南方品種，但混有39 份北方品種；紫色部分共有42 份品種，其中大部分為山東資源，占78.57%，在地理位置上正好也為南北兩大類群的交匯地帶。由此推測，全部資源大致分為南北兩個群體，但也存在較多的混雜，南北方板栗品種存在較多基因交流，山東地區(qū)為南北品種類群的混交地帶。群體結(jié)構(gòu)分析進(jìn)一步證明，全部板栗資源大致分為兩部分，但存在較多混雜。群體結(jié)構(gòu)中混雜系數(shù)大于0.8 被認(rèn)為存在較多基因交流，如北方品種群中混入的南方品種Y140 系數(shù)大于0.8，因此會聚類到北方品種群中，表明該品種存在頻繁的基因交流。

2.2 核心種質(zhì)構(gòu)建與評價(jià)

2.2.1 構(gòu)建必選種質(zhì)庫 種質(zhì)資源庫對于育種和基因挖掘具有重要作用。如圖2 所示，根據(jù)板栗品種獨(dú)特的成熟期、枝條長度和彎曲度、刺蓬色澤等，從全部資源中挑選出8 個具有特殊性狀的品種作為必選品種，而這些資源會在未來育種和遺傳研究等方面發(fā)揮重要作用。其中，早熟型品種有京暑紅、燕山早豐；短枝型品種有短枝1 號、短枝3 號和超短枝2號；垂枝型品種有垂枝2 號；雄花不育型品種有無花；紅刺型品種有山東紅栗。

2.2.2 核心種質(zhì)構(gòu)建 為了確定最合適的核心種質(zhì)構(gòu)建策略以及取樣數(shù)量，筆者基于最大等位基因數(shù)，利用分層抽樣法、模擬退火算法和隨機(jī)搜索算法3種方法構(gòu)建最優(yōu)核心種質(zhì)（圖3）。

利用分層抽樣策略，通過逐步聚類取樣形成了8 個取樣比例為100%、83.33%、66.67%、46.49%、31.58%、20.76%、13.74%、8.77%，并計(jì)算每個取樣比例下的等位位點(diǎn)數(shù)量分別為212、207、202、196、185、179、163、152。隨著取樣比例的升高，所保留的等位位點(diǎn)數(shù)量也逐漸增多，當(dāng)取樣比例達(dá)到100%時，等位位點(diǎn)數(shù)量達(dá)到最大為212。

RS，AN，AN為利用基于等位位點(diǎn)數(shù)量的隨機(jī)搜索算法來構(gòu)建核心種質(zhì)。以5 為梯度計(jì)算等位位點(diǎn)數(shù)量，結(jié)果表明，等位位點(diǎn)數(shù)量呈現(xiàn)波浪變化，當(dāng)取樣數(shù)量達(dá)到325 份時所保留的等位位點(diǎn)數(shù)量達(dá)到最大值212。

SA，AN，AN為利用基于等位位點(diǎn)數(shù)量的模擬退火算法來構(gòu)建核心種質(zhì)。結(jié)果表明當(dāng)以5 為梯度計(jì)算等位位點(diǎn)數(shù)量時，隨著取樣數(shù)量的增加等位位點(diǎn)數(shù)量也逐漸升高。當(dāng)取樣數(shù)量達(dá)到85 份資源時，所保留的等位位點(diǎn)數(shù)達(dá)到最大值212，并且隨著取樣數(shù)量繼續(xù)增加時，等位位點(diǎn)數(shù)量趨于平穩(wěn)，不再升高或降低。

綜上所述，通過對3 種方法在不同取樣數(shù)量時的等位位點(diǎn)數(shù)比較分析，結(jié)果表明SA，AN，AN策略為構(gòu)建核心種質(zhì)的最優(yōu)策略，且當(dāng)種質(zhì)資源數(shù)量達(dá)到85 份時能夠代表全部資源的等位基因數(shù)，即所保留的等位位點(diǎn)數(shù)為212。

2.2.3 核心種質(zhì)評價(jià)及遺傳多樣性分析 為了驗(yàn)證這85 份核心種質(zhì)資源的合理性、可靠性及科學(xué)性，首先，筆者對核心種質(zhì)在全部資源的分布情況進(jìn)行主坐標(biāo)分析（PCoA）以及對核心種質(zhì)進(jìn)行群體結(jié)構(gòu)分析。如圖4 所示，Pop2 部分為核心種質(zhì)，Pop1 部分為篩選之后的剩余種質(zhì)，主坐標(biāo)1 和2 分別解釋了位點(diǎn)信息數(shù)據(jù)中34.36%和15.71%的變異，初步顯示核心種質(zhì)在全部種質(zhì)資源中較為均勻的分布，從而表明所篩選結(jié)果的合理性。群體結(jié)構(gòu)分析表明，85個核心種質(zhì)資源與全部種質(zhì)資源資源相同，大致分為兩部分。

隨后筆者對所篩選的核心種質(zhì)和原始種質(zhì)之間的遺傳多樣性參數(shù)進(jìn)行分析。如表1 所示，原始種質(zhì)共342 份，核心種質(zhì)共85 份，保留率為24.85%；核心種質(zhì)的等位位點(diǎn)數(shù)量和原始種質(zhì)的等位位點(diǎn)數(shù)相同，保留率為100%；等位基因數(shù)（Na）、有效等位基因數(shù)（Ne）、觀測雜合度（Ho）、期望雜合度（He）、多態(tài)信息含量（PIC）的保留率分別為86.32%、85.46%、98.15%、90.12%和100.79%。

對核心種質(zhì)和原始種質(zhì)的遺傳多樣性參數(shù)進(jìn)行T 檢驗(yàn)，結(jié)果顯示所構(gòu)建的核心種質(zhì)各遺傳多樣性參數(shù)與原始種質(zhì)之間沒有顯著差異，且其原生地涵蓋11 個省份，其他省份樣品較少，分析中沒有入選，表明所構(gòu)建的85 個核心種質(zhì)資源能夠很好地保存原始種質(zhì)資源的遺傳多樣性（表2）。

3 討論

3.1 板栗品種遺傳多樣性分析

基于SSR 標(biāo)記研究種群遺傳多樣性已成為最方便快捷的方法，前人[28-33]均利用SSR 技術(shù)研究種群遺傳多樣性。而關(guān)于板栗遺傳多樣性的研究也已有報(bào)道。張馨方等[34]利用SSR 分子標(biāo)記技術(shù)對279份板栗資源進(jìn)行遺傳多樣性分析，平均每對引物檢測到3.38 個等位位點(diǎn)，有效等位基因數(shù)為1.53，聚類結(jié)果表明群體間遺傳關(guān)系和地理來源有一定聯(lián)系。此外，江錫兵等[35]、張馨方等[36]均采用SSR標(biāo)記對部分地區(qū)板栗資源進(jìn)行遺傳多樣性分析。筆者對全國板栗品種資源進(jìn)行廣泛收集，并且基于SSR 熒光分子標(biāo)記技術(shù)對其進(jìn)行遺傳多樣性分析，總等位位點(diǎn)為212 個，平均每對引物檢測到的等位基因數(shù)為10.095，有效等位基因數(shù)為3.488，明顯高于前人的研究結(jié)果。因此，筆者在本研究中的21 對SSR分子標(biāo)記能夠反映出板栗品種的遺傳多樣性。此外，對342 份板栗品種資源進(jìn)行聚類分析，結(jié)果表明中國板栗品種資源大致分為南北兩個群體。但是存在大部分南北混雜品種，對于其混雜原因還有待進(jìn)一步研究分析。另外，選取聚類分析中得到的南北方板栗品種的含水量、可溶性糖含量、直鏈淀粉含量、支鏈淀粉含量、總淀粉含量等19 個表型性狀進(jìn)行差異顯著性分析，結(jié)果顯示，南北方板栗品種的可溶性糖含量、直鏈淀粉含量、支鏈淀粉含量、刺苞厚、堅(jiān)果厚、堅(jiān)果寬呈顯著性差異；總苞質(zhì)量、苞質(zhì)量、刺苞寬、刺苞高、苞肉厚、單粒質(zhì)量呈極顯著差異；含水量、總淀粉含量、出實(shí)率等其他指標(biāo)無顯著差異。此外，上述指標(biāo)中除了直鏈淀粉含量、出實(shí)率和果形指數(shù)外，其余南方品種的表型指標(biāo)值均大于北方品種，如含水量、可溶性糖含量、單粒質(zhì)量等指標(biāo)值，推測可能與南北方的環(huán)境差異有關(guān)。通過表型分析，進(jìn)一步證明板栗品種（系）大致分為南北兩個群體。

3.2 板栗品種核心種質(zhì)構(gòu)建分析

種質(zhì)資源的收集和保存對物種的遺傳多樣性保護(hù)和育種具有重要作用，但是一個物種的種質(zhì)資源往往比較豐富，不方便進(jìn)行遺傳育種研究[12]。因此，選擇一定方法構(gòu)建核心種質(zhì)對深入挖掘優(yōu)異種質(zhì)資源，提高資源利用價(jià)值和效率具有重要意義[1]。構(gòu)建策略的選擇是核心種質(zhì)的重點(diǎn)，不同的取樣策略對核心種質(zhì)樣品的選擇具有直接的影響。此外，等位基因的數(shù)量是多樣性的反映[37]，保留更多的等位基因?qū)τN有重要作用[38]。李金龍等[5]基于聚類和表型，構(gòu)建了甜蕎530 份種質(zhì)；李洪果等[39]基于等位基因最大化原則構(gòu)建了杜仲189 份種質(zhì)；王春梅等[40]、陳天青等[41]、何建文等[42]都用了一種方法構(gòu)建了玉米、小麥、辣椒等的核心種質(zhì)。因此，為了使核心種質(zhì)更具有代表性和可靠性，筆者在本研究中采用3 種方法進(jìn)行核心種質(zhì)構(gòu)建，同時比較3 種方法所保留的等位位點(diǎn)數(shù)量，進(jìn)而篩選出最少的種質(zhì)資源以最大程度地代表原始種質(zhì)的遺傳多樣性。

中國板栗品種豐富，資源眾多，不易保存，因此構(gòu)建核心種質(zhì)尤為重要。前人對于板栗核心種質(zhì)也有研究，如張馨方等[23]采用多次聚類的方法構(gòu)建了燕山板栗的核心種質(zhì)，不足以代表全國范圍內(nèi)的板栗品種資源；Nie 等[22]基于最大遺傳多樣性的模擬退火算法構(gòu)建了45份板栗品種的核心種質(zhì)，但是只運(yùn)用一種取樣策略且研究的樣品較少，具有局限性。因此亟需構(gòu)建涵蓋全國板栗品種資源遺傳多樣性的核心種質(zhì)。

筆者通過SA，AN，AN、RS，AN，AN和分層抽樣3 種方法構(gòu)建核心種質(zhì)，最終選擇SA，AN，AN作為最優(yōu)策略并且構(gòu)建的85 份資源約占原始種質(zhì)庫的24.85%，能夠充分代表原始種質(zhì)的遺傳多樣性，符合前人關(guān)于核心種質(zhì)取樣比例在5%～40%的研究結(jié)論[43]?；诘任换驍?shù)的模擬退火算法相比另外兩種構(gòu)建策略而言，能夠以最少的資源代表整個原始種群的遺傳多樣性，且能夠更好地保存資源。隨著市場的需要及科技的發(fā)展，會不斷出現(xiàn)新的優(yōu)異板栗品種資源，因此構(gòu)建核心種質(zhì)是一個動態(tài)的過程，在今后的研究中需要不斷完善核心種質(zhì)。

4 結(jié)論

筆者在本研究中基于SSR 分子標(biāo)記，揭示了南北方板栗品種（系）存在較多基因交流，同時利用SA，AN，AN、RS，AN，AN和分層抽樣3 種策略構(gòu)建核心種質(zhì)，確定基于等位基因數(shù)的模擬退火算法為最優(yōu)取樣策略，最終構(gòu)建85 份核心種質(zhì)。