張 旭, 蘇世平, 師微檸, Alizet Didi Dom, 馬 強

(甘肅農(nóng)業(yè)大學， 甘肅 蘭州 730070)

紅砂(Reaumuriasoongorica)隸屬于檉柳科(Tamari-caceae)石竹亞綱(Caryophyllidae)[1]，主要分布于亞歐大陸(中亞)、南歐和北非，共有21種，亞洲有18種，我國分布有4種，多產(chǎn)自內(nèi)蒙古和西北地區(qū)[2]。紅砂是我國干旱荒漠區(qū)分布最廣的植物種之一，也是荒漠化、半荒漠化地區(qū)防風固沙和水土保持的優(yōu)良灌木樹種，和荒漠綠洲與荒漠過渡帶的重要植被類型，其具有抗旱、耐旱、耐鹽堿及很強的集沙固沙能力[3-4]。

近年來，由于干旱和高溫等極端氣候的頻發(fā)，導致紅砂種子繁殖困難，長期的無性繁殖導致遺傳多樣性降低[5]，因此探討紅砂種群的遺傳變異規(guī)律尤為重要，為其遺傳資源的保護和評價提供理論依據(jù)。李秀玲等[6]運用空間自相關(guān)對我國西北地區(qū)的28個紅砂種群進行空間分布的研究，結(jié)果表明紅砂種群的空間分布模式為隨機分布。徐莉等[7]應用了RAPD分子標記，15對引物對分布于新疆阜康荒漠7個亞群的136株紅砂進行擴增。Shannon多樣性指數(shù)(0.307)和居群間遺傳分化系數(shù)(0.312)揭示了紅砂的遺傳變異多存在于亞種群內(nèi)，亞種群間的遺傳分化較小。殷恒霞等[8]采用磁珠富集法構(gòu)建了紅砂基因組微衛(wèi)星富集文庫，分析得出微衛(wèi)星重復基元顯示為5種類型，其中三堿基基元(GGT)n,(TGG)n,(ACC)n最多，占總數(shù)的60.7%。楊九艷等[9]運用ISSR對內(nèi)蒙古高原荒漠區(qū)紅砂居群的遺傳多樣性進行研究，結(jié)果表明紅砂居群間有高度的遺傳分化。張穎娟等[10]應用RAPD分子標記對內(nèi)蒙古西部紅砂種群進行遺傳多樣性分析，結(jié)果表明不同生境下紅砂種群內(nèi)存在較高的遺傳多樣性，聚類分析表明遺傳距離與地理距離之間存在一定相關(guān)性。

本研究基于前人的研究基礎(chǔ)，應用SSR分子標記對來源于不同區(qū)域的16個紅砂居群進行多樣性水平評估及群體遺傳結(jié)構(gòu)分析，從分子水平上探究紅砂居群的遺傳規(guī)律，以期為紅砂的合理利用及有效保護提供理論基礎(chǔ)。盡管在紅砂遺傳多樣性研究進展中有相關(guān)研究，但是都基于小地域，而對全國大地域分布區(qū)的研究報道較少[11-14]。應用的SSR分子標記具有高度重復性，豐富多態(tài)性及共顯性的優(yōu)點，為目前應用比較廣泛的分子標記技術(shù)[15]。

1 材料與方法

1.1 試驗材料

2011年11月在我國甘肅、寧夏、內(nèi)蒙古的紅砂天然分布區(qū)的各居群中，隨機選取紅砂植株并采集種子，標記、分裝；2012年1月在武威市林業(yè)技術(shù)服務中心苗圃利用所采集的種子進行穴盤育苗；2013年4月土壤解凍后進行移栽，并按常規(guī)方法進行育苗期管理，建立了位于武威市涼州區(qū)下雙鎮(zhèn)國家樟子松基地向北1 000米的紅砂試驗地。本研究實驗材料均取自該試驗基地。

表1 供試材料紅砂不同居群信息Table 1 Information of different populations of R. Soongorica

1.2 DNA提取及擴增

用磁珠法基因組提取試劑盒(天根生化科技有限公司)根據(jù)產(chǎn)品說明按步驟進行總DNA的提取。從96對引物中篩選出11對多態(tài)性較高的引物(SSR引物由武漢天一輝遠生物科技有限公司合成)用于PCR擴增。反應體系為15 μL：2×Taq PCR Master Mix7.5 μL，Mix primer(引物混合液)2.0 μL,DNA template 1 μL,ddH2O 4.5 μL。反應程序:95℃預變性5 min；95℃變性30 s，52℃～62℃退火30 s，72℃延伸30 s,共10個循環(huán)；95℃變性30 s,52℃退火30 s,72℃延伸30 s,共22～25個循環(huán)；72℃延伸2 min;4℃保存[11-12]。

1.3 數(shù)據(jù)分析

通過群體遺傳學軟件GeneAlex6.41計算每個位點的等位基因數(shù)(Na)、有效等位基因數(shù)(Ne)、期望雜合度(Expected heterozygosity,He)、觀察雜合度(Observed heterozygosity,Ho)、基因分化系數(shù)(Fst)、基因流(Gene Flow,Nm)、香農(nóng)信息指數(shù)(I)等[16]；并對群體間進行分子方差分析(AMOVA)[17-18]、居群間的遺傳距離(Nei′s Genetic Distance,D)和遺傳一致度(Nei′s Genetic Identity,GI)[19]。基于Nei′s遺傳距離，用Phylip軟件構(gòu)建UPGMA聚類圖[20]。

有效等位基因數(shù)為純合度的倒數(shù)，是反應群體遺傳變異大小的指標之一，其值越接近等位基因數(shù)的絕對值，表明等位基因在群體中分布越均勻[21]。式中，Pi為第i個位點上等位基因出現(xiàn)的頻率。

Ne=1/1-He

香農(nóng)信息指數(shù)用于估算居群內(nèi)的遺傳分化,指數(shù)越大居群間或居群內(nèi)的分化程度越高，且遺傳多樣性越高[22]。

I=-∑pilnpi

期望雜合度(Expected heterozygosity,He)，根據(jù)每個等位基因的基因頻率計算，期望雜合度的值越高，表明物種多態(tài)性越高[23]。

多態(tài)信息含量指數(shù)用來衡量基因座多態(tài)性高低的程度[24]。式中，n代表的是樣本中序列的數(shù)量。

自然群體并非是理想群體，如果僅考慮單基因座上一對等位基因，對于由若干個地方群體組成的生物大群體，各層次現(xiàn)實群體中基因型頻率期望值偏差一般用Fst度量。計算公式如下[25]：

式中，群體內(nèi)平均雜合度(The mean Heterozygosity within Populations,Hs),為各個居群的期望雜合度的算數(shù)平方根。

整體期望雜合度(Total Expected Heterozygosity,HT)的計算公式如下：

基因流代表等位基因從一個居群到另一個居群間的轉(zhuǎn)移?；蛄鞔笥?，就能發(fā)揮均質(zhì)化作用，能有效抑制有遺傳漂變而引起的遺傳分化[26]。

Nm=[(1/FST)-1]/4

2 結(jié)果與分析

2.1 位點多態(tài)性評估

11對引物擴增出67個等位基因，其中多態(tài)等位基因為56個，多態(tài)率為83.58%。擴增長度為126～396 bp，位點ReS078上擴增出的等位基因數(shù)最多，位點ReS008最少。多態(tài)含量指數(shù)為0.061～0.910，均值為0.396，位點ReS078的多態(tài)含量指數(shù)(0.910)為最高，位點ReS008的PIC指數(shù)(0.061)為最低。紅砂基因位點的多態(tài)性水平為中度多態(tài)性(表2)。

2.2 居群多樣性分析

所測紅砂居群的等位基因數(shù)為2.364～4.455，均值為3.735，有效等位基因數(shù)為1.504～2.585，均值為2.221，結(jié)果顯示居群的觀測等位基因數(shù)均高于有效等位基因數(shù)(表3)。其中居群THCY的Na(4.455)和Ne(2.585)為最高，居群MQXX的Na(2.364)和Ne(1.504)為最低。觀測雜合度為0.295～0.418，均值為0.358。居群LZ的Ho(0.418)為最高，表明該居群內(nèi)雜合子占比為41.8%；MQXX的Ho(0.286)為最低，說明該居群內(nèi)純合子占比較多。居群觀測雜合度和期望雜合度均小于0.5。僅有MQXX觀測雜合度大于期望雜合度，在該居群內(nèi)存在雜合度缺失的現(xiàn)象。香農(nóng)信息指數(shù)為0.420～0.873，均值為0.763，說明居群內(nèi)的遺傳分化程度低，按照香農(nóng)多樣性指數(shù)評估紅砂居群的遺傳多樣性水平：THCY,ZZG>LZ>NMYH>QYS>HSZ>WWYH>MJW>SSC>BYHT>HG>XGG>SGK>WD>SMY>MQXX。

表2 11對引物多樣性評估Table 2 Diversity Evaluation of 11 pairs of primers

表3 居群各項多樣性參數(shù)(均值)Table 3 Values of population diversity parameters(Mean value)

2.3 遺傳分化及基因流

紅砂整體近交系數(shù)為-0.026～0.548，整體表現(xiàn)為近親交配，但在地方亞群內(nèi)存在異交現(xiàn)象(表4)。居群間遺傳分化系數(shù)為0.041～0.21，平均值為0.088，說明居群間遺傳分化程度較低。AMOVA結(jié)果表明，居群間的遺傳差異占變異總來源的5.8%，個體間占14.8%，個體內(nèi)的遺傳差異占79.4%，得出紅砂的遺傳變異材料主要來源于紅砂個體內(nèi)。各居群基因交流為0.940～5.807，平均值為3.163，該結(jié)果表明所測紅砂居群間有較強的基因交流。

表4 各位點F-Statistics統(tǒng)計值和基因流Table 4 Each point F-statistics Statistics and gene flow

2.4 遺傳結(jié)構(gòu)及Mantel檢測

紅砂居群間的遺傳距離為0.009～0.229，均值為0.062，遺傳距離最近的居群是WD和SMY，最遠的居群是MQXX和HG(表5)。遺傳距離小于0.1的居群占總數(shù)的76.7%，大于0.1的占總數(shù)的23.3%，表明所測紅砂居群間的遺傳距離相近。居群間遺傳一致度為0.795～0.983，均值為0.933，其中遺傳相似性最大的居群是WD和SMY，最小的是MQXX和HG。Mantel檢測(圖1)得出紅砂居群的遺傳距離和地理距離顯著性相關(guān)(R2=0.3254,P=0.02<0.05)。

表5 遺傳距離(左下)和遺傳一致度(右上)Table 5 Enetic Distance (Lower left) and Genetic Identity (Upper right)

圖1 遺傳距離與地理距離的Mantel檢測Fig.1 Mantel test of Genetic distance and geographical distance

2.5 聚類分析及PCoA主成分分析

基于紅砂居群間的遺傳距離并采用UPGMA進行聚類(圖2)。SSC,SGK,WD,SMY,THCY,BYHT,NMYH,MJW,QYS,WWYH,ZZG被聚為一類；然后和MQXX聚成第二類；和HSZ聚為第三類；再和XGG,HG聚為第四類；最后和遺傳距離較遠的LZ聚為第五類。PCoA分析結(jié)果(圖3)將紅砂居群被分為5個部分。第一部分為HG,XGG；第二部分為NMYH；第三部分為LZ；第四部分為MQXX；其余居群為第五部分。PCoA與UPGMA聚類結(jié)果基本一致。

圖2 不同居群紅砂UPGMA聚類圖Fig.2 UPGMA clustering diagram of R. soongorica in different cohabitation groups

圖3 不同居群的PCoA分析Fig.3 PCoA analysis of different populations

2.6 STRUCTURE結(jié)構(gòu)分析

當K=2時，299份材料被分為2個亞群；當K=3時，能將MQXX,HSZ和XGG,HG區(qū)分開；模擬k值2～20，得出的最大似然值為5，說明299份紅砂材料被分為5個亞群(圖4)。第Ⅰ亞群(綠色)共75份材料，占總數(shù)的25%，第Ⅱ亞群(藍色)共75份紅砂，占總體的25%，第Ⅲ亞群(粉色)包含66份材料，占總數(shù)的22.1%，第Ⅳ亞群(橙色)包含材料56份，占比為18.7%。第Ⅴ亞群(黃色)共27份，占比為9%。聚類分析將299份材料分為5個亞群(圖5)，第Ⅰ類中LZ為4份，THCY,HSZ各3份，ZZG,MJW,SSC,QYS,SGK各2份。第Ⅱ類中與XGG,HG(占總數(shù)37.8%)分在一起的有46份材料，各占比為NMYH(9.2%),WWYH(7.8%),THCY(6.5%)。第Ⅲ類以MJW為主共計聚類23份材料。第Ⅳ類共有76份紅砂材料，其中ZZG(13.1%)，HSZ(15.7%)，NMYH(9.2%)，SSC(9.2%)聚在一起。第Ⅴ類共有106份紅砂材料，占比為MQXX(15%)，SMY(13.2%)，WD(11.3%)，BYHT和SGK(各為8.4%)，THCY,SSC和QYS(各占6.6%),NMYH,ZZG(各占5.6%),MJW(4.7%),WWYH(2.8%),HSZ(1.8%),XGG(0.9%)。利用聚類分析校對了紅砂的亞群分類，兩者結(jié)果顯示299份紅砂材料被分為5個亞群。

圖4 不同k值的遺傳結(jié)構(gòu)與Δk值的折線圖Fig.4 Genetic structure of different K values and the magnitude of Δk

圖5 299份紅砂材料的聚類分析Fig.5 Clustering results of 299 R. Soongorica materials

3 討論

3.1 居群遺傳多樣性

本研究的11對引物能成功擴增出3～9個等位基因，與張穎娟[10]和馮亮亮[14]的研究結(jié)果大致相同。等位基因多態(tài)率為83.58%，多態(tài)信息含量指數(shù)為0.395，得出紅砂的基因位點具有較高的多態(tài)性，在今后的育種工作中可以利用該物種豐富的遺傳材料，通過雜交、篩選等方法獲得生態(tài)效益更高的紅砂品種[27]。紅砂居群的等位基因數(shù)均高于有效等位基因數(shù)，表明紅砂的等位基因在染色體上分布不均勻，紅砂居群在雜交過程中有明顯的、通過染色體較大片段而進行遺傳傳遞的現(xiàn)象。該現(xiàn)象可能會導致在遺傳時發(fā)生后代遺傳信息偏親遺傳的現(xiàn)象[28]?；螂s合度是衡量居群遺傳變異的參數(shù)之一[29]。經(jīng)研究表明紅砂居群的雜合度均較低，其中居群MQXX的觀測雜合度(0.286)最低，表明該居群內(nèi)雜合率低，在雜交過程中產(chǎn)生的雜合度丟失的現(xiàn)象[30]，推測產(chǎn)生該現(xiàn)象的原因與紅砂的繁殖特性有關(guān)，位于干旱、沙化生境的成熟階段的紅砂居群主要以無性繁殖為主[3]，而該紅砂分布區(qū)常年高溫少雨、干旱[31]，因此居群內(nèi)主要以根基劈裂和不定根的方式進行繁殖，使紅砂居群內(nèi)雜合度缺失及分化程度降低，進而導致對生存環(huán)境的適應性降低。本研究表明地方內(nèi)紅砂居群有一定程度的異交，這種交配所產(chǎn)生的后代具有雜合基因型。雜合基因型是產(chǎn)生基因交換、重組和產(chǎn)生新基因型的基礎(chǔ)，這種交配方式的存在，提高了紅砂居群內(nèi)的基因雜合度，在紅砂的進化過程中為其創(chuàng)造更廣泛的遺傳變異材料[32-33]。分化系數(shù)是從兩個地方群體中任意抽取兩個配子是同源的概率，可以度量紅砂居群間的遺傳分化程度[34]。結(jié)果表明紅砂居群間的遺傳分化主要來源于個體內(nèi)DNA堿基序列的差異，AMOVA分析表明居群間遺傳變異僅占遺傳總材料的5.8%,這種結(jié)果與前人研究結(jié)果一致[10,35]。

3.2 遺傳結(jié)構(gòu)

本研究表明WD和SMY的遺傳距離最近(0.007)，MQXX和HG的遺傳距離最遠(0.229)，Mantel檢測得出紅砂居群的遺傳距離和地理距離顯著相關(guān)(R2=0.3254,P=0.02<0.05)，由此類現(xiàn)象推測出可能隨著物種的進化所產(chǎn)生的地理隔離現(xiàn)象，經(jīng)歷不同的自然選擇、不同的進化方向、進而產(chǎn)生種內(nèi)差異而形成的地理小種[36-37]。聚類分析表明，來自騰格里沙漠東南邊沿的SSC,SGK,WD,SMY,THCY,BYHT,NMYH,MJW,QYS被聚在第一類，這些居群在地理位置上較為集中，并且在紅砂散粉期[38](7-8月)受中國夏季風的影響[39](風向由南向北，為夏季風盛行期)，導致這些居群間有較頻繁的基因交流，這種現(xiàn)象表明臨近紅砂居群間可以通過風媒促進居群間的基因交流，從而改變居群內(nèi)遺傳結(jié)構(gòu)及富集基因庫。來源于祁連山中段的LZ居群在遺傳結(jié)構(gòu)上其他居群分化明顯，推測該現(xiàn)象是由于地理分布不同引起的種間差異[40]，郝曉莉等[41]的研究結(jié)果表明五柱紅砂(R.kaschgaricaRupr.)分布延伸至祁連山中段，因此LZ和其他居群的差異極有可能是由屬內(nèi)的不同種間引起的。居群QYS和MJW間的基因交流較頻繁、遺傳分化較小，極有可能是居群間能夠自由交配的兩個亞群，基因流在居群間是隨機的、均一的，使得居群間相似性變大，進而導致居群間表現(xiàn)為基因頻率和基因型頻率的哈迪-溫伯格平衡[42]?；蛄魅趸俗匀贿x擇和遺傳漂變的作用，使居群間趨于一致，更有助于紅砂群體內(nèi)部能夠穩(wěn)定的生存和發(fā)育[43]。引物不能很好的區(qū)分THCY和BYHT，想要進一步了解紅砂居群的遺傳多樣性及親緣關(guān)系，需開發(fā)鑒別性更高的引物或在此基礎(chǔ)上增加其他引物[44]。通過種質(zhì)資源的遺傳多樣性分析，可有效降低資源利用過程中的盲目性?；诓煌尤哼z傳結(jié)構(gòu)表明，XGG,HSZ與其他居群遺傳距離相對較遠，居群內(nèi)遺傳多樣性豐富，在進行紅砂雜交組合選擇時，可利用這2個居群的材料與其他材料進行組配，從而為育種工作者親本選配提供科學依據(jù)[45]。

4 結(jié)論

本研究運用SSR分子標記技術(shù)，對16個不同地理來源的紅砂居群進行遺傳多樣性及遺傳結(jié)構(gòu)分析。研究結(jié)論得出紅砂的基因多態(tài)性為中度，等位基因不均勻的分布在染色體上。觀測雜合度和F指數(shù)說明居群內(nèi)雜合度缺失，純合體過量。按照香農(nóng)多樣性指數(shù)評估紅砂居群遺傳多樣性水平：THCY,ZZG>LZ>NMYH>QYS>HSZ>WWYH>MJW>SSC>BYHT>HG>XGG>SGK>WD>SMY>MQXX。Fst指數(shù)表明居群間分化程度較低，且紅砂的遺傳變異主要來自源于個體內(nèi)DNA序列的不同。各居群間基因交流為3.163，說明居群間有較頻繁的基因交流。Mantel檢測得出遺傳距離和地理距離顯著相關(guān)。UPGMA將16個紅砂居群分為5大類，聚類分析表明居群間的基因交流是影響居群聚類的重要因素。STRUCTURE結(jié)構(gòu)分析將材料分為5個亞群。各亞群中，第Ⅴ亞群內(nèi)的紅砂個體血統(tǒng)較為純正，第Ⅱ亞群內(nèi)的紅砂個體血統(tǒng)混雜。