樣本數(shù)量對(duì)白樺群體遺傳參數(shù)估算的影響

畢志宏，魏敏靜，劉瑩瑩，楊傳平，魏志剛（東北林業(yè)大學(xué) 林木遺傳育種國(guó)家重點(diǎn)實(shí)驗(yàn)室，黑龍江 哈爾濱 150040）

樣本數(shù)量對(duì)白樺群體遺傳參數(shù)估算的影響

畢志宏，魏敏靜，劉瑩瑩，楊傳平，魏志剛
（東北林業(yè)大學(xué) 林木遺傳育種國(guó)家重點(diǎn)實(shí)驗(yàn)室，黑龍江 哈爾濱 150040）

群體遺傳參數(shù)是群體遺傳學(xué)研究的主要內(nèi)容之一，估算值的大小受群體樣本數(shù)量的影響。以帽兒山試驗(yàn)林場(chǎng)白樺Betula platyphylla種源試驗(yàn)林為對(duì)象，采用序列相關(guān)擴(kuò)增多態(tài)性（sequence-related amplified polymorphism，SRAP）技術(shù)分析了種源內(nèi)不同樣本數(shù)量對(duì)白樺群體遺傳多樣性的影響。結(jié)果表明：不同樣本數(shù)量對(duì)白樺的估算的各遺傳參數(shù)有一定影響，但影響各異。其中，對(duì)白樺群體多態(tài)位點(diǎn)百分率、群體內(nèi)遺傳變異和群體間遺傳參數(shù)影響顯著，對(duì)總位點(diǎn)數(shù)、多態(tài)位點(diǎn)數(shù)和基因流影響不顯著。樣本數(shù)量對(duì)于群體間與群體內(nèi)遺傳變異的大小的研究結(jié)果也產(chǎn)生影響，如樣本數(shù)量超過(guò)8個(gè)，顯示白樺種源內(nèi)的遺傳變異幅度大于種源間；而樣本數(shù)量為4個(gè)，顯示為種源間遺傳變異幅度大于種源內(nèi)。此外，不同樣本數(shù)量揭示出的種源間親緣關(guān)系不同，如樣本數(shù)量超過(guò)8個(gè)，15個(gè)種源的聚類(lèi)結(jié)果基本一致；而樣本數(shù)量為4個(gè)，聚類(lèi)結(jié)果無(wú)法反映白樺群體的親緣關(guān)系。最后，提出了白樺遺傳多樣性研究時(shí)樣本數(shù)量的計(jì)算公式，認(rèn)為白樺多態(tài)位點(diǎn)百分率達(dá)98.5%的樣本數(shù)量應(yīng)為11個(gè)。圖4表3參23

林木遺傳育種學(xué)；白樺；遺傳多樣性；樣本容量

林木種質(zhì)資源是陸地生物的避護(hù)所和主基因庫(kù)。林木種質(zhì)資源的丟失，將引起鄰近生物種及其種質(zhì)以4～13倍的速率丟失［1］。了解林木群體的遺傳背景、遺傳結(jié)構(gòu)及群體間的親緣關(guān)系，是林木遺傳多樣研究和采取科學(xué)有效措施保護(hù)其種質(zhì)資源的前提和基礎(chǔ)［2］，也為林木資源的開(kāi)發(fā)與利用提供信息。到目前為止，先后采用不同的標(biāo)記技術(shù)對(duì)眾多的林木遺傳多樣性進(jìn)行了研究［3-6］。群體遺傳參數(shù)是群體遺傳學(xué)研究的主要內(nèi)容之一，估算值的大小受群體樣本數(shù)量的影響。林木為異交物種，需要一定的數(shù)量才能代表1個(gè)群體，因此群體內(nèi)個(gè)體數(shù)量的多少可能對(duì)群體遺傳參數(shù)估算產(chǎn)生一定影響。然而，長(zhǎng)期以來(lái)，在樹(shù)木遺傳多樣性研究中，群體內(nèi)個(gè)體數(shù)量對(duì)遺傳參數(shù)研究結(jié)果的影響研究尚無(wú)報(bào)道。序列相關(guān)擴(kuò)增多態(tài)性（sequence-related amplified polymorphism，SRAP）標(biāo)記建立以來(lái)，由于其操作簡(jiǎn)便、成本低、可信度高、易于測(cè)序等特點(diǎn)倍受分子生物學(xué)家的青睞［7］，被迅速應(yīng)用到花生Arachis hypogaea，油菜Brassica napus等植物［8-14］的種質(zhì)資源鑒定評(píng)價(jià)、遺傳圖譜構(gòu)建、重要性狀標(biāo)記以及基因克隆等方面。白樺Betula platyphylla是一種分布范圍廣、適用性強(qiáng)、用途廣泛的闊葉速生樹(shù)種，然而，有關(guān)白樺遺傳多樣性研究中，同樣也沒(méi)有重視種源內(nèi)樣本數(shù)量對(duì)遺傳參數(shù)估算的影響問(wèn)題［15］。本研究以東北林業(yè)大學(xué)帽兒山試驗(yàn)林場(chǎng)15個(gè)白樺種源為對(duì)象，采用SRAP技術(shù)分析了白樺種源內(nèi)不同樣本數(shù)量對(duì)其遺傳多樣性參數(shù)估算的影響，以期確定白樺遺傳多樣性分析的種源內(nèi)個(gè)體最優(yōu)數(shù)量，也為其他樹(shù)種遺傳多樣性研究提供參考。

1　材料與方法

1.1植物材料

試驗(yàn)材料來(lái)自東北林業(yè)大學(xué)帽兒山實(shí)驗(yàn)林場(chǎng)內(nèi)白樺種源試驗(yàn)林，共15個(gè)種源，分別為涼水、汪清、新疆、露水河、小北湖、輝南、桓仁、莫爾道嘎、草河口、綽爾、寧夏、清源、東方紅、帽兒山、烏伊嶺。隨機(jī)選取個(gè)體20個(gè)·種源-1，共計(jì)300個(gè)樣本。2012年6月采取當(dāng)年生嫩葉，用密封袋封好，做好標(biāo)記，放在冰盒中帶回，置于－70℃冰箱中保存，以備提取DNA。

1.2試驗(yàn)方法

1.2.1白樺葉片基因組DNA提取白樺葉片基因組DNA提取參考文獻(xiàn)［16］。

1.2.2SRAP標(biāo)記分析SRAP反應(yīng)體系為20.00 μL：0.21×106μg·L－1DNA，1.50 μL；25 mmol·L－1Mg2+，1.40 μL；5×16.67 μkat·L－1Taq酶，0.25 μL；2.5 mmol·L－1三磷酸堿基脫氧核苷酸（dNTPs），2.00 μL；10 μmol·L－1引物，0.35 μL。聚合酶鏈?zhǔn)椒磻?yīng)（PCR）反應(yīng)程序?yàn)椋?4℃預(yù)變性5 min，94℃變性1 min，35℃復(fù)性1 min，72℃延伸1 min，5個(gè)循環(huán)；94℃變性1 min，50℃復(fù)性1 min，72℃延伸1 min，30個(gè)循環(huán)，72℃延伸7 min。根據(jù)文獻(xiàn)［17］中的引物序列，本研究中所用SRAP引物由上海生工合成。從親緣關(guān)系最近的2個(gè)種源（莫爾道嘎和桓仁）中分別選5個(gè)個(gè)體用于多態(tài)性引物的篩選。

1.2.3遺傳多樣性參數(shù)的估算15個(gè)白樺種源，每個(gè)種源按包含4，8，12和20個(gè)個(gè)體為研究對(duì)象，利用篩選出的17對(duì)引物對(duì)其擴(kuò)增后分析各項(xiàng)的遺傳參數(shù)。電泳圖譜中的每條帶均代表了引物與模板DNA互補(bǔ)的一對(duì)結(jié)合位點(diǎn)，可記為1個(gè)分子標(biāo)記。根據(jù)分子量標(biāo)準(zhǔn)對(duì)照反應(yīng)產(chǎn)物在膠上的對(duì)應(yīng)位置，估計(jì)擴(kuò)增產(chǎn)物的分子量大小，有帶的記為 “1”，無(wú)帶的記為 “0”。所得結(jié)果為一二元數(shù)據(jù)矩陣，利用Popgene 32軟件計(jì)算各項(xiàng)遺傳參數(shù)。利用NTsys 2.10e軟件分析各源間聚類(lèi)關(guān)系，利用Mintable 15軟件分析群體多態(tài)位點(diǎn)百分?jǐn)?shù)與群體內(nèi)個(gè)體數(shù)量之間的關(guān)系和作圖。

2　結(jié)果與分析

2.1基因組DNA的質(zhì)量檢測(cè)及多態(tài)性SRAP引物篩選

使用十六烷基三甲基溴化銨（CTAB）法提取的白樺基因組DNA在質(zhì)量濃度為10.0 g·L-1瓊脂糖凝膠電泳成像后如圖1所示。樣品DNA為1條清晰明亮的條帶，無(wú)明顯拖尾、彌散或其他異常現(xiàn)象。隨機(jī)抽取部分樣品經(jīng)紫外分光光度計(jì)檢測(cè)，D（260）/D（280）值約1.8。因此可以認(rèn)為提取的DNA純度較高，符合SRAP-PCR擴(kuò)增要求［18］。

前期基于表型研究的白樺種源實(shí)驗(yàn)結(jié)果表明：莫爾道嘎和桓仁親緣關(guān)系最近，因此，本研究中，從這2個(gè)種源內(nèi)各隨機(jī)選擇5個(gè)個(gè)體用于引物的篩選。以條帶清晰并且群體與個(gè)體間有差異條帶的引物為選擇標(biāo)準(zhǔn)，從合成的30對(duì)引物中進(jìn)行篩選。圖2中5個(gè)泳道為1組擴(kuò)增結(jié)果，從左開(kāi)始分別為莫爾道嘎和桓仁的種源及個(gè)體，圖2為引物me2em3組合擴(kuò)增結(jié)果。從圖2可看出：該引物擴(kuò)增時(shí)，不僅2個(gè)種源間存在差異性條帶，且同一種源內(nèi)不同個(gè)體間也有清晰差異條帶。這表明引物me2em3符合引物選擇的標(biāo)準(zhǔn)，可用于白樺群體遺傳多樣性的研究。按此標(biāo)準(zhǔn)，從合成的30對(duì)引物中，篩選出17個(gè)條帶清晰、多態(tài)性高的引物用于實(shí)驗(yàn)分析（表1）。

圖1　部分白樺樣本基因組DNA的瓊脂糖凝膠電泳檢測(cè)結(jié)果Figure 1　Detection of Betula platyphlla genomic DNA using agarose gel electrophoresis

圖2　白樺SRAP引物（me2em3）的篩選Figure 2　Primer selected for SRPA of B.platyphlla

表1　SRAP引物序列Table 1　primer sequences of SRAP

2.2白樺種源不同樣本數(shù)量估算的群體遺傳多樣性

表2　不同樣本數(shù)量估算的遺傳參數(shù)比較Table 2　Comparison of different sample size estimation for genetic parameters

利用17對(duì)SRAP引物分別對(duì)包含4，8，12和20個(gè)樣本個(gè)體的白樺共15個(gè)種源進(jìn)行擴(kuò)增后分析，并利用POP gene軟件計(jì)算各遺傳參數(shù)。表2為不同樣本數(shù)量時(shí)白樺群體遺傳多樣性各參數(shù)的平均值。結(jié)果發(fā)現(xiàn)：樣本數(shù)量不同時(shí)，各項(xiàng)遺傳參數(shù)平均數(shù)也存在一定差異。這說(shuō)明樣本數(shù)量對(duì)白樺群體的遺傳參數(shù)估算有一定的影響。如當(dāng)群體樣本數(shù)量分別為4，8，12，16和20時(shí)，各群體的多態(tài)位點(diǎn)分別為129，132，133，133和134個(gè)，多態(tài)位點(diǎn)百分率為95.56%～99.26%。平均等位基因數(shù)為1.955 6～1.985 2，有效等位基因數(shù)為1.420 7～1.427 0，Nei's指數(shù)的變動(dòng)范圍為0.256 9～0.259 3，Shannon指數(shù)的變動(dòng)范圍為0.400 8～0.404 1。進(jìn)一步分析發(fā)現(xiàn)，在不同樣本數(shù)量條件下，各項(xiàng)遺傳參數(shù)變異最大種源均為涼水，如當(dāng)群體樣本數(shù)量為 4時(shí)，共檢測(cè)到 135個(gè)位點(diǎn)，其中多態(tài)位點(diǎn)有 129個(gè)，多態(tài)位點(diǎn)百分率達(dá)95.56%。多態(tài)位點(diǎn)百分率為 17.78%～52.33%，平均等位基因數(shù)為1.177 8～1.533 3，有效等位基因數(shù)為1.112 6～1.356 6，Nei's指數(shù)的變動(dòng)范圍為0.066 0～0.205 9，Shannon指數(shù)的變動(dòng)范圍為0.098 3～0.304 2，上述的遺傳參數(shù)均表現(xiàn)為涼水種源各項(xiàng)遺傳參數(shù)最大，且在15個(gè)種源間變化趨勢(shì)一致。當(dāng)群體樣本數(shù)為8，12，16，20時(shí)涼水種源各項(xiàng)遺傳參數(shù)均最大。而各項(xiàng)遺傳參數(shù)最小的種源隨樣本數(shù)量的變異均不同，如樣本數(shù)量為12時(shí)，共檢測(cè)到135個(gè)位點(diǎn)，其中多態(tài)位點(diǎn)有133個(gè)，多態(tài)位點(diǎn)百分率達(dá)98.52%，多態(tài)位點(diǎn)百分率為34.07%～77.78%，平均等位基因數(shù)為1.340 7～1.777 8，有效等位基因數(shù)為1.166 7～1.420 9，Nei's指數(shù)為0.101 5～0.249 2，Shannon指數(shù)為0.156 6～0.378 5。上述遺傳參數(shù)中，汪清種源各項(xiàng)遺傳參數(shù)最?。欢鴺颖緮?shù)量為16時(shí)，共檢測(cè)到135個(gè)位點(diǎn)，其中多態(tài)位點(diǎn)有 133個(gè)，多態(tài)位點(diǎn)百分率達(dá)98.52%，種源的多態(tài)位點(diǎn)數(shù)為47～108，多態(tài)位點(diǎn)百分率為34.81%～80.00%，平均等位基因數(shù)為1.348 1～1.800 0，有效等位基因數(shù)為1.180 0～1.425 4，Nei's指數(shù)為0.108 1～0.251 1，Shannon指數(shù)為0.165 0～0.381 2，上述遺傳參數(shù)中，清源種源的各項(xiàng)遺傳參數(shù)最??；樣本數(shù)量為8和4時(shí)，各項(xiàng)遺傳參數(shù)最小的為淖爾種源。此外，4個(gè)樣本數(shù)量與其他樣本數(shù)量之間的平均等位基因數(shù)，多態(tài)位點(diǎn)數(shù)，多態(tài)位點(diǎn)百分率，Ht，Hs，Gst和Nm差異較大。

表2同時(shí)也表明：樣本數(shù)量對(duì)白樺群體內(nèi)與群體間遺傳變異大小的研究結(jié)果也有影響。如當(dāng)群體樣本數(shù)量為4時(shí)，47.78%的遺傳變異存在于種源內(nèi)，52.22%的遺傳變異存在于種源間。當(dāng)群體內(nèi)樣本數(shù)量超過(guò)8時(shí)，研究反映出白樺種群內(nèi)遺傳變異幅度超過(guò)群體間；如群體樣本數(shù)量為8時(shí)，58.82%的遺傳變異存在于種源內(nèi)，41.18%的遺傳變異存在于種源間；當(dāng)群體樣本數(shù)量為20時(shí)，65.88%的遺傳變異存在于種源內(nèi)，34.12%存在于種源間。

2.3不同樣本數(shù)量對(duì)群體間親緣關(guān)系分析結(jié)果的影響

為了進(jìn)一步分析種源內(nèi)樣本數(shù)量差異對(duì)于白樺遺傳結(jié)構(gòu)研究結(jié)果的影響，筆者又研究了不同樣本數(shù)量下的白樺種源遺傳進(jìn)化關(guān)系。圖3表明：樣本數(shù)量對(duì)于白樺種源遺傳進(jìn)化關(guān)系分析結(jié)果有一定的影響。如種源樣本數(shù)量超過(guò)12時(shí)，15個(gè)種源的聚類(lèi)結(jié)果基本一致（圖3 C，圖3D和圖3E），均可在遺傳距離為0.17時(shí)，分成4個(gè)類(lèi)群，而且每個(gè)類(lèi)群包括的種源基本一致。樣本數(shù)量為8時(shí)，在0.17的遺傳距離處，15個(gè)種源則劃分為5個(gè)類(lèi)群，其中涼水種源單獨(dú)聚成一類(lèi)（圖3B）。而種源樣本數(shù)量為4時(shí)，在0.17的遺傳距離處，15個(gè)種源被劃分成7個(gè)類(lèi)群，反映不出群體進(jìn)化的基本規(guī)律（圖3A）。此時(shí)，只有當(dāng)遺傳距離達(dá)到0.25左右時(shí)，15個(gè)種源才可劃分成4類(lèi)，但涼水種源單獨(dú)聚成一類(lèi)，帽兒山、東方紅和烏伊嶺種源聚成一類(lèi)，新疆和露水河種源聚成一類(lèi)，其余9個(gè)種源聚成一類(lèi)。由此可見(jiàn)：當(dāng)樣本數(shù)量?jī)H有4個(gè)時(shí)，其揭示的白樺群體間的親緣關(guān)系與樣本數(shù)量為8個(gè)和超過(guò)12時(shí)存在差異。而且，種源樣本數(shù)量為4時(shí)，15個(gè)種源遺傳距離變化范圍為0.050 6～0.280 0，其中清源和寧夏種源的遺傳距離最小，桓仁和烏伊嶺種源的遺傳距離最大。當(dāng)種源樣本數(shù)量為分別8，12，16和20，15個(gè)種源間在0.044 0～0.212 5范圍內(nèi)變化，而此時(shí)淖爾和涼水種源間的遺傳距離最大，莫爾道嘎和桓仁種源的遺傳距離最小。

2.4樣本數(shù)量和主要遺傳參數(shù)的相關(guān)分析

由于樣本數(shù)量對(duì)白樺遺傳參數(shù)與遺傳結(jié)構(gòu)的分析結(jié)果均產(chǎn)生一定的影響，因此為弄清樣本數(shù)量與各遺傳參數(shù)之間的具體關(guān)系，分析了不同樣本數(shù)量與其對(duì)應(yīng)的各項(xiàng)遺傳參數(shù)作相關(guān)性。結(jié)果（表3）表明：種源內(nèi)不同樣本數(shù)量與白樺群體多態(tài)位點(diǎn)百分率、群體內(nèi)遺傳變異和群體間遺傳變異參數(shù)相關(guān)性顯著，與總位點(diǎn)數(shù)、多態(tài)位點(diǎn)數(shù)和基因流的相關(guān)性不顯著。這進(jìn)一步表明，在研究異交群體遺傳多樣性時(shí)，應(yīng)科學(xué)估算出群體樣本的最佳數(shù)量。

2.5白樺遺傳多樣性研究時(shí)樣本數(shù)量的確定

表3　樣本容量與主要遺傳參數(shù)的相關(guān)性Table 3　Correlation between sample size and main genetic parameters

白樺種源樣本數(shù)量與一些遺傳參數(shù)之間存在顯著或極顯著相關(guān)性，這為確定白樺遺傳多樣性研究時(shí)，確定合理的樣本數(shù)量提供了可能。在群體遺傳多樣性研究中，多態(tài)位點(diǎn)百分率越高越有利于分析和研究群體遺傳多樣和遺傳結(jié)構(gòu)，因此對(duì)白樺種源不同樣本數(shù)量與多態(tài)位點(diǎn)百分率進(jìn)行擬合分析（圖4），在95%的置信區(qū)間內(nèi)，兩者之間呈一元二次曲線關(guān)系，決定系數(shù)為89%，表明兩者之間存在較好的相關(guān)性。我們?cè)诎讟暹z傳多樣性研究時(shí)，在綜合考慮成本、技術(shù)難度等基礎(chǔ)上，根據(jù)所需多態(tài)位點(diǎn)百分率大小，計(jì)算出所需樣本的數(shù)量。如當(dāng)白樺遺傳多樣性研究時(shí)，以最低所需多態(tài)位點(diǎn)百分率98.5%為標(biāo)準(zhǔn)，則每個(gè)種源所需樣本數(shù)量為10.673個(gè)，實(shí)際操作中可設(shè)為11個(gè)樣本個(gè)體。

圖3　不同樣本數(shù)量時(shí)15個(gè)白樺種源聚類(lèi)圖Figure 3　Different number of samples for 15 birch provenance clustering map

3　討論與結(jié)論

對(duì)于異交植物，種群內(nèi)和種群間遺傳變異較大，必須有一定數(shù)量個(gè)體才能代表該種群的基因庫(kù)，因此群體樣本數(shù)量會(huì)對(duì)研究結(jié)果產(chǎn)生影響。如HENSEN等［19］在研究白鮮Dictamnus albus時(shí)發(fā)現(xiàn)，其遺傳多樣性及多態(tài)位點(diǎn)隨群體樣本容量的不同而發(fā)生變化，而且群體大小、多態(tài)位點(diǎn)百分率和遺傳多樣性間存在顯著的相關(guān)性。本研究表明：白樺種源內(nèi)樣本數(shù)量對(duì)白樺群體多態(tài)位點(diǎn)百分率、群體內(nèi)遺傳變異和群體間遺傳參數(shù)影響顯著，對(duì)總位點(diǎn)數(shù)、多態(tài)位點(diǎn)數(shù)和基因流影響不顯著。ZHAO等［20］研究發(fā)現(xiàn)不同樣本容量的野生大豆Glycine soja種群間得到的遺傳多樣性指數(shù)存在較大差異。本研究同樣發(fā)現(xiàn)：白樺種源樣本數(shù)量不同時(shí)，其估算的遺傳多樣性指數(shù)也存在一定差異。如樣本數(shù)量為4時(shí)，多態(tài)位點(diǎn)數(shù)和百分率最大與最小的種源分別為涼水和淖爾，而樣本數(shù)量超過(guò)8個(gè)時(shí)，則為涼水和清源。此外，種源內(nèi)不同個(gè)體數(shù)量反映的種源內(nèi)與種源間的遺傳變異不同，如種源內(nèi)樣本數(shù)為4時(shí)，白樺種源間的遺傳變異大于種源內(nèi)。而當(dāng)群體樣本數(shù)量超過(guò)8時(shí)，白樺種源內(nèi)的遺傳變異幅度大于種源間。

本研究發(fā)現(xiàn)：群體樣本數(shù)量對(duì)遺傳群體間遺傳距離的分析結(jié)果產(chǎn)生一定的影響，如源樣本數(shù)量為4時(shí)，清源和寧夏種源的遺傳距離最小，桓仁和烏伊嶺種源的遺傳距離最大。而樣本數(shù)量超過(guò)8時(shí)，淖爾和涼水的遺傳距離最大，莫爾道嘎和桓仁的遺傳距離最小。這與樣本數(shù)量超過(guò)10時(shí)，通過(guò)表型分析得到的研究結(jié)果相同，如高玉池等［21］利用表型性狀對(duì)帽兒山地區(qū)10年生白樺種源間的親緣關(guān)系進(jìn)行了分析，結(jié)果表明淖爾和涼水的遺傳距離最大，莫爾道嘎和桓仁的遺傳距離相對(duì)較小。此外，樣本數(shù)量也對(duì)白樺種源間的遺傳進(jìn)化關(guān)系分析結(jié)果有一定的影響，如在相同的遺傳距離處（0.17），樣本數(shù)量為4時(shí)，15個(gè)種源被劃分成7個(gè)類(lèi)群；樣本數(shù)量為8時(shí)，15個(gè)種源被劃分成5個(gè)類(lèi)群；樣本數(shù)量超過(guò)12時(shí)，15個(gè)種源被劃分成4個(gè)類(lèi)群。上述分析表明：群體樣本數(shù)量不僅對(duì)白樺各群體的遺傳參數(shù)產(chǎn)生影響，也對(duì)群體間遺傳進(jìn)化關(guān)系的分析結(jié)果有一定的影響。

圖4　白樺種源樣本數(shù)量與多態(tài)位點(diǎn)百分率的關(guān)系Figure 4　Relationship between Betula platyphylla provenances samples and the percentage of polymorphic loci

異交植物多少個(gè)體才能代表一個(gè)群體，才能得到較為客觀真實(shí)的遺傳多樣性參數(shù)在草本植物或花卉植物中有所研究。如SINGH等［22］在研究野生二粒小麥Triticum dicoccoides的遺傳多樣性時(shí)發(fā)現(xiàn)，有的群體5～6個(gè)個(gè)體就可以使多樣性標(biāo)準(zhǔn)差不大于10%，個(gè)別一些則需12個(gè)以上的個(gè)體。JULIO等［23］發(fā)現(xiàn)，薔薇Rosaceae種群個(gè)體數(shù)超過(guò)10個(gè)以后，就不能提供更多的種群雜合度信息。本研究發(fā)現(xiàn)：白樺種源內(nèi)樣本的數(shù)量與其多態(tài)位點(diǎn)百分率呈一元二次曲線相關(guān)，且決定系數(shù)達(dá)89%。并以此為公式計(jì)算白樺遺傳多樣性研究時(shí)，如設(shè)定要探測(cè)到的多態(tài)位點(diǎn)百分率達(dá)98.5%以上時(shí)，每個(gè)種源至少必須含有11個(gè)樣本個(gè)體。

［1］魏志剛，高玉池，劉桂豐，等.白樺核心種質(zhì)初步構(gòu)建［J］．林業(yè)科學(xué)，2009，45（10）：74－80. WEI Zhigang，GAO Yuchi，LIU Guifeng，et al.Preliminary construction of core collection of Betula platyphylla germplasm［J］.Sci Silv Sin，2009，45（10）：74－80.

［2］顧萬(wàn)春．中國(guó)林木種質(zhì)資源保存研究與對(duì)策［C］//陳宜瑜.中國(guó)生物多樣性保護(hù)與研究進(jìn)展：第5屆全國(guó)生物多樣性保護(hù)與持續(xù)利用研討會(huì)論文集.北京：氣象出版社，2004：93－105.

［3］羅建勛，王國(guó)良，辜云杰，等．杉木優(yōu)樹(shù)半同胞子代測(cè)定及優(yōu)良家系選擇初步研究［C］//中國(guó)林學(xué)會(huì).第2屆中國(guó)林業(yè)學(xué)術(shù)大會(huì)：S2功能基因組時(shí)代的林木遺傳與改良論文集.南寧：［s.n.］.2009：671－673.

［4］李斌，顧萬(wàn)春.松屬植物遺傳多樣性研究進(jìn)展［J］．遺傳，2003，25（6）：740－748. LI Bin，GU Wanchun.Review on genetic diversity in Pinus［J］.Hereditas，2003，25（6）：740－748.

［5］張萍，周志春，金國(guó)慶，等.木荷種源遺傳多樣性和種源區(qū)初步劃分［J］．林業(yè)科學(xué)，2006，42（2）：38－42. ZHANG Ping，ZHOU Zhichun，JIN Guoqing，et al.Genetic diversity analysis and provenance zone allocation of Schima superba in China using RAPD markers［J］.Sci Silv Sin，2006，42（2）：38－42.

［6］莊振宏.運(yùn)用RAPD技術(shù)分析長(zhǎng)葉榧（Torreya jackii Chun）的遺傳多樣性［D］．福州：福建師范大學(xué)，2002. ZHUANG Zhenhong.Genetic Diversity Analysis by RAPD on Torreya jackii Chun［D］.Fuzhou：Fujian Normal University，2002.

［7］譚碧玥，王源秀，徐立安．分子標(biāo)記SRAP及其在林木研究中的應(yīng)用［J］．世界林業(yè)研究，2009，22（5）：45－50. TAN Biyue，WANG Yuanxiu，XU Li'an.Molecular marker SRAP and its application in forest research［J］.World For Res，2009，22（5）：45－50.

［8］劉月光，滕永勇，潘辰，等．應(yīng)用SRAP標(biāo)記對(duì)蓮藕資源的聚類(lèi)分析［J］．氨基酸和生物資源，2006，28（1）：29－32. LIU Yueguang，TENG Yongyong，PAN Chen，et al.Cluster analysis of Nelumbo based on SRAP markers［J］.AminoAcid Biol Resour，2006，28（1）：29-32.

［9］史紅麗，韓明玉，趙彩平.桃遺傳多樣性的SRAP和SSR標(biāo)記分析［J］.華北農(nóng)學(xué)報(bào)，2009，24（6）：187-192. SHI Hongli，HAN Mingyu，ZHAO Caiping.Genetic diversity of Prunus persica using SRAP and SSR markers［J］. Chin J North China，2009，24（6）：187-192.

［10］陳靜，胡曉輝，苗華榮，等.花生SRAP-PCR技術(shù)體系的優(yōu)化［J］.山東農(nóng)業(yè)科學(xué)，2008，10（5）：24-26. CHEN Jing，HU Xiaohui，MIAO Huarong，et al.Optimization of SRAP-PCR reaction system on peanut［J］.Shandong Agric Sci，2008，10（5）：24-26.

［11］張羽，李英，許偉，等.利用SRAP和SSR標(biāo)記對(duì)漢中地區(qū)主要油菜品種的遺傳多樣性分析［J］.四川農(nóng)業(yè)大學(xué)學(xué)報(bào)，2013，31（4）：370-376. ZHANG Yu，LI Ying，XU Wei，et al.Analysis of genetic diversity of Brassica napus mainly popularized in Hanzhong based on SRAP and SSR markers［J］.J Sichuan Agric Univ，2013，31（4）：370-376.

［12］王惠哲，李淑菊，管煒.應(yīng)用SRAP標(biāo)記分析黃瓜的遺傳差異［J］.生物技術(shù)通報(bào)，2009（2）：76-79. WANG Huizhe，LI Shuju，GUAN Wei.Analysis on the genetic relationship of four types of cucumber germphasms by SRAP markers［J］.Biotechnol Bull，2009（2）：76-79.

［13］張凱，羅小敏，蔣玉春，等.甘薯種質(zhì)資源的SRAP鑒定及遺傳多樣性分析［J］.核農(nóng)學(xué)報(bào)，2013，27（5）：568-575. ZHANG Kai，LUO Xiaomin，JIANG Yuchun，et al.Genetic diversity among main sweeppotato resources using SRAP markers［J］.J Nucl Agric Sci，2013，27（5）：568-575.

［14］李朋飛，霍秀愛(ài)，程永強(qiáng)，等.基于SRAP的西瓜種質(zhì)資源遺傳多樣性評(píng)價(jià)［J］.中國(guó)農(nóng)業(yè)科技導(dǎo)報(bào)，2013，15（2）：89-96. LI Pengfei，HUO Xiuai，CHENG Yongqiang，et al.Assement of genetic diversity in watermelon based on SRAP analysis［J］.J Agric Sci Technol，2013，15（2）：89-96.

［15］姜靜，楊傳平，劉桂豐，等.應(yīng)用RAPD技術(shù)對(duì)東北地區(qū)白樺種源遺傳變異的分析［J］.東北林業(yè)大學(xué)學(xué)報(bào)，2001，29（2）：30-34. JIANG Jing，YANG Chuanping，LIU Guifeng，et al.Analysis on genetic variation within and among Betula platyphylla provenances in northeast China using RAPD markers［J］.J Northeast For Univ，2001，29（2）：30-34.

［16］詹亞光，曾凡鎖.富含多糖的白樺成熟葉片DNA的提取方法［J］.東北林業(yè)大學(xué)學(xué)報(bào)，2005，33（3）：24-25. ZHAN Yaguang，ZENG Fansuo.A method of DNA extraction from mature birch leaves rich in polysaccharide［J］.J Northeast For Univ，2013，33（3）：24-25.

［17］任娜.紅豆杉內(nèi)生真菌遺傳多樣性分析及生物轉(zhuǎn)化研究［D］.長(zhǎng)沙：中南大學(xué)，2011. REN Na.Genetic Diversity Analysis and Biotransformation Research among Endophytic Fungi from Taxus［D］. Changsha：Central South University，2011.

［18］張景閣.金磁微粒純化細(xì)菌基因組DNA方法的研究［D］.西安：西北大學(xué)，2011. ZHANG Jingge.Study on the Genomic DNA Method by Goldmag Particles［D］.Xi'an：Northwestern University，2011.

［19］HENSEN I，OBERPRIELER C.Effects of population size on genetic diversity and seed production in the rare Dictamnus albus（Rutaceae）in central Germany［J］.Conserv Genet，2005，6（1）：63-73.

［20］ZHAO Ru，CHENG Zhou，LU Weifeng，et al.Estimating genetic diversity and sampling strategy for a wild soybean （Glycine soja）population based on different molecular markers［J］.Chin Sci Bull，2006，51（10）：1219-1227.

［21］高玉池，魏志剛，楊傳平，等.帽兒山地區(qū)10年生白樺種源試驗(yàn)［J］.浙江林學(xué)院學(xué)報(bào)，2009，26（6）：784-791. GAO Yuchi，WEI Zhigang，YANG Chuanping，et al.A 10-year-old provenance trial of Betula platyphylla in the Maoershan area of Heilongjiang Province［J］.J Zhejiang For Coll，2009，26（6）：784-791.

［22］SINGH M，CHABANE K，VALKOUN J，et al.Optimum sample size for estimating gene diversity in wild wheat using AFLP markers［J］.Genet Resour Crop Evol，2006，53（1）：23-33.

［23］JULIO N，SOBRAL A，BUE?AS J R，et al.RAPD and ISSR markers indicate diminished gene flow due to recent fragmentation of Polylepis australis woodlands in central Argentina［J］.Biochem Syst Ecol，2008，36（5/6）：329-335.

Genetic parameters of Betula platyphylla provenances due to sample size

BI Zhihong，WEI Minjing，LIU Yingying，YANG Chuanping，WEI Zhigang
（State Key Laboratory of Forest Genetics and Breeding，Northeast Forestry University，Harbin 150040，Heilongjiang，China）

Population genetic parameters，one of the main constituents of population genetics，and the size of estimated values are influenced by the number of population samples.In this study Betula platyphylla （birch）from Cap Mountain Experimental Forest Farm's Provenance Testing Forest was chosen as the experimental object and sequence related amplified polymorphism（SRAP）technology was used for analysis.Results showed that genetic parameters with different sample sizes had different effects on birch.Genetic variation，population of birch groups，and the percentage of polymorphic numbers within populations influenced genetic parameters.Influence from genetic variation and genetic parameters between groups was strong，but total numbers，polymorphic numbers，and gene flow were not affected.When the sample size was greater than eight，genetic variation amplitude of the birch source was greater than that of the provenances；when the sample size was four，genetic variation was greater than the source of genetic variation among provenances for different samples.For genetic relationships of different provenances with sample size greater than eight，15 kinds of source clustering results were basically the same.However，with a sample size of four，clustering results did not reflect the evolutionary relationship between birch groups.Also，11 birch samples were found to have a polymorphic number that was 98.5%.［Ch，4 fig.3 tab.23 ref.］

forest tree breeding；Betula platyphylla；genetic diversity；sample capacity

S722.3

A

2095-0756（2016）04-0564-07

10.11833/j.issn.2095-0756.2016.04.003

2015-06-10；

2015-10-09

國(guó)家高技術(shù)研究發(fā)展計(jì)劃（“863”計(jì)劃）項(xiàng)目（2011AA100202）

畢志宏，從事木材次生生長(zhǎng)等研究。E-mail：bizhihong9191@163.com。通信作者：魏志剛，教授，博士，從事林木遺傳育種等研究。E-mail：zhigangwei1973@163.com