油松天然次生林居群遺傳多樣性及與產(chǎn)地地理氣候因子的關(guān)聯(lián)分析

李 明，王樹香，高寶嘉

(1.河北農(nóng)業(yè)大學(xué)生命科學(xué)學(xué)院，保定 071000;2.河北農(nóng)業(yè)大學(xué)林學(xué)院，保定 071000)

油松(Pinus tabulaeformis Carr.)是我國特有樹種，分布全國14個省(自治區(qū))［1］，具有重要的生態(tài)、經(jīng)濟價值，是我國北方天然油松分布區(qū)的主要建群種，同時也是我國北方地區(qū)治理水土流失、改善生態(tài)環(huán)境的重要造林樹種［2］。近年來，人們對油松天然林進行了大量利用，加之氣候影響，大面積的油松天然林遭到破壞，甚至消亡。因而，迫切需要對油松天然林遺傳資源進行保護，以防止其優(yōu)良基因和遺傳多樣性喪失。在上世紀中后期徐化成等從形態(tài)特征、化學(xué)成分及分布特征等方面做了大量關(guān)于油松天然林的地理變異研究［3-5］，認為生長特性與緯度呈負相關(guān)，與年均溫呈正相關(guān)，同時形態(tài)指標、抗逆指標又存在地理變異的不連續(xù)性，并對油松種源的20個特征進行了遺傳變異研究，劃分了7個生態(tài)型。在油松的分子遺傳方面，張春曉、李悅等采用水平淀粉凝膠電泳方法研究了油松遺傳多樣性指數(shù)與海拔、年降雨量存在相關(guān)關(guān)系［6-7］;王磊等利用隨機微衛(wèi)星擴增多態(tài)性DNA(RMAPD)分子標記技術(shù)，揭示了油松各天然居群間存在較高遺傳多樣性差異［8］;李毳、柴寶峰等用酸性聚丙烯酰胺凝膠(APAGE)技術(shù)分析了華北地區(qū)油松居群在醇溶蛋白及分子水平的遺傳特征［9-10］。近年來，李明、張靜潔等采用ISSR標記技術(shù)對油松種質(zhì)資源保護、遺傳多樣性與生境環(huán)境作了一定研究，認為不同海拔生境及土壤因子影響下，其天然次生林存在明顯變異規(guī)律［11-12］;王孟本、Eckert等通過對不同居群油松與環(huán)境因子研究，認為不同油松居群間遺傳分化顯著，遺傳多樣性指數(shù)與年均溫、年降雨量不存在顯著相關(guān)性［13-15］。系統(tǒng)研究產(chǎn)地環(huán)境因子與油松天然次生林遺傳變異的關(guān)聯(lián)以及地理變異規(guī)律將對油松種質(zhì)資源的保護和利用具有重要的理論意義。本研究采用ISSR技術(shù)對油松天然群體遺傳多樣性變異進行研究，據(jù)此分析其與產(chǎn)地生態(tài)環(huán)境因子的關(guān)聯(lián)，從而更清楚地認識和預(yù)測天然油松次生林分子特征的地理變異規(guī)律，并為正確、合理利用與保護天然油松種質(zhì)資源提供基礎(chǔ)資料和科學(xué)依據(jù)。

1 研究地區(qū)與研究方法

1.1 野外取樣

2011年7底到8月初，于河南、山東、河北、寧夏、山西、陜西、甘肅、遼寧等8省10地區(qū)的油松天然次生林中取樣。每居群取25株個體，每兩個體間距大于50 m，共250個樣本。采集當年生新鮮針葉，將針葉放入保鮮盒中帶回，置于-20℃冰箱保存。其中，經(jīng)緯度、海拔由GPS測得，年平均氣溫、極端最高溫、極端最低溫、年降雨量、年日照時數(shù)等氣象因子由各居群所在地氣象站資料計算得到，各采樣點位置及氣象數(shù)據(jù)如表1。

1.2 研究方法

1.2.1 DNA的提取和反應(yīng)體系

采用UNIQ-10柱式新型植物基因組DNA抽提試劑盒(購自上海生物工程有限公司)提取DNA，并稀釋成終濃度20 ng/μL，放入-20℃冰箱儲存?zhèn)溆谩８鶕?jù)加拿大哥倫比亞大學(xué)公布的100條ISSR標記序列(http://www.biotech.ub.ca/services/naps/primers/Primers.pdf)，篩選出 13 個重復(fù)性好、特異性高的引物用于 ISSRPCR反應(yīng)［12］，ISSR引物由上海生物工程有限公司合成。擴增產(chǎn)物經(jīng)1.5%瓊脂糖凝膠電泳分離，溴化乙錠染色后，在UVP凝膠成像系統(tǒng)下觀察拍照。

1.2.3 數(shù)據(jù)統(tǒng)計和分析

以產(chǎn)生分子量為100—2000 bp的Marker DL2000(寶泰克)為標準譜帶，有帶記作“1”，無帶記作“0”。假設(shè)群體處于哈迪-韋伯平衡，利用POPGENE 1.32和GENAlEX6.3軟件對譜帶進行遺傳特征分析，分別計算:(1)觀測等位基因數(shù)Na;(2)有效等位基因數(shù)Ne;(3)期望雜合度He;(4)Shannon信息指數(shù)I;(5)多態(tài)位點百分率P;(6)總基因多樣度Ht;(7)居群內(nèi)基因多樣度Hs。

應(yīng)用SPSS 18.0軟件進行Pearson相關(guān)性分析，應(yīng)用POPGENE軟件以UPGMA法對居群間的親緣關(guān)系進行聚類分析，并繪制聚類圖。

2 結(jié)果與分析

2.1 油松天然居群遺傳多樣性

13個ISSR引物共擴增出137個多態(tài)位點，平均多態(tài)位點百分率60.72%，總的多態(tài)位點百分率100%，平均每條引物擴增10.54條譜帶(表2)，絕大多數(shù)譜帶集中在400—1600 bp(圖1)。不同油松居群的總基因多樣度(Ht)為0.3461，其中存在于居群內(nèi)的基因多樣度(Hs)為0.2162，居群間的基因多樣度(DST)為0.1299。不同居群間的遺傳分化系數(shù)(GST)為0.3753，表明大部分分子變異發(fā)生在居群內(nèi)(62.47%)。一般風媒傳粉的異交木本植物居群間變異占總變異的10%左右，油松與之相比，顯然是偏高的。根據(jù)Nm=0.5(1-Gst)/Gst計算出不同生境下的油松居群基因流(Nm)為0.8323，說明居群間存在有限的基因交流。

表2 ISSR擴增引物、序列及擴增結(jié)果Table 2 ISSR primer labels，primer sequences and amplification results

圖1 ISSR引物(826)對蘇峪口居群部分個體的擴增結(jié)果Fig.1 The profiles of SDC population using the ISSR primer(826)

由表3可知，各居群間 Shannon信息指數(shù)變化范圍在0.1923—0.2490，平均值(0.2165)與 Nybom［21］所統(tǒng)計的多種植物種群水平的遺傳多樣性(基于RAPD和ISSR等顯性標記)平均值(0.22)相似;各居群間期望雜合度He變化范圍在0.2824—0.3702，其中山西和順云龍公園居群最高，甘肅冶力關(guān)居群最低，與觀測等位基因數(shù)變化趨勢基本一致;多態(tài)位點百分率以河南寶天曼油松居群最高(71.53%)，山西和順云龍公園居群次之(70.7%)，甘肅冶力關(guān)居群最低(51.09%)。

表3 油松天然居群遺傳參數(shù)Table 3 Genetic parameters of Pinus tabulaeformis populations

2.2 油松居群ISSR差異帶數(shù)分析

10個油松天然居群中，13條引物擴增出差異帶數(shù)最少為84條(DWP)，最多為127條(YLGY)，各油松天然居群擴增譜帶及頻率存在較大的變異。另外，在差異帶數(shù)出現(xiàn)頻率很低(≤25%)時，少數(shù)居群如YLG、SYK、DWP三居群沒有擴增特異條帶，說明三居群變異程度低;而YLGY居群有3條，YWL和LKS居群各有2條，其它居群各1條(表4)，說明這些油松天然次生林居群具有更豐富的變異基礎(chǔ)。在使用的引物中，SDC居群有1條特異性擴增條帶，為該居群種質(zhì)資源的變異特性研究提供了重要前提。

表4 油松居群ISSR差異帶數(shù)分析Table 4 Different bands for ISSR data of Pinus tabulaeformis populations

2.3 油松天然居群間的遺傳距離和UPGMA分析

表5顯示，10個天然油松居群間的遺傳距離(D)變化范圍為0.083(LKS-CJC)—0.2939(SDC-DWP)，地理距離變化范圍為110.9252 km(LKS-SDC)—1741.0933 km(YWL-YLG)。經(jīng)Mantel檢驗，10個居群間的地理距離和遺傳距離不存在顯著相關(guān)性(r=0.069，P=0.360)，說明地理隔離不是遺傳變異的主要原因，而產(chǎn)地環(huán)境因子長期脅迫，可能使油松居群發(fā)生一定程度的地理適應(yīng)性分化。

表5 油松居群間遺傳距離(對角線下方)和地理距離(對角線上方)/kmTable 5 Nei's genetic distance(below diagonal)and geographic distance(above diagonal)/km of Pinus tabulaeformis populations

根據(jù)Nei無偏遺傳距離構(gòu)建居群間UPGMA聚類圖(圖2)顯示，10個油松天然居群共聚成3大類:河南寶天曼(BTM)、承德大窩鋪(DWP)、寧夏蘇峪口(SYK)和甘肅冶力關(guān)(YLG)居群聚為一類，遼寧醫(yī)巫閭(YWL)、山西沁源靈空山(LKS)、陜西蔡家川林場(CJC)、山西和順云龍公園(YLGY)和山西汶水三道川林場(SDC)居群聚為一類，山東蒙山(MS)居群獨立為一類。

圖2 10個天然油松居群基于Nei's(1972)遺傳距離聚類圖Fig.2 Dendrogram based Nei's(1972)genetic distance of 10 Pinus tabulaeformis populations

2.4 油松天然居群遺傳多樣性及與產(chǎn)地環(huán)境因子相關(guān)性分析

進一步對油松天然居群遺傳多樣性與產(chǎn)地環(huán)境因子進行相關(guān)性分析(表6)，油松天然居群的觀測等位基因數(shù)(Na)及多態(tài)條帶百分率(P)與年平均氣溫、1月平均氣溫、極端最低溫等產(chǎn)地氣候因子呈顯著正相關(guān)，與海拔呈顯著負相關(guān)。Shannon信息指數(shù)(I)與海拔(r=-0.5788)、極端最高溫(r=0.5289)及年降雨量(r=-0.5512)呈顯著相關(guān)，即隨著極端最高溫升高、海拔降低、年降雨量減少，Shannon信息指數(shù)升高。

表6 油松居群遺傳參數(shù)與環(huán)境因子相關(guān)性Table 6 The relationships between genetic diversity parameters of Pinus tabulaeformis populations and environmental factors

3 結(jié)論與討論

3.1 油松天然居群的遺傳多樣性與環(huán)境因子相關(guān)性

油松居群的平均多態(tài)位點百分率60.72%，總的多態(tài)位點比率為100%，變異幅度范圍為55.47%—71.53%，表明油松在種水平上具有較高的遺傳多樣性。Wright［17］認為居群間基因流大于1能發(fā)揮其均質(zhì)化作用，反之，則表明基因流成為遺傳分化的主要原因。本研究天然油松林的基因流為0.8323，說明不同產(chǎn)地環(huán)境的長期影響，可能造成天然油松居群間基因流較小，是居群間遺傳分化重要原因。居群間遺傳變異占總變異的37.53%，表明居群間已產(chǎn)生較大程度的分化。故表明居群間地理及環(huán)境因子的差異，可能使花粉的傳播受到較大程度的阻礙，從而限制了居群間的隨機交配，居群間遺傳分化越來越明顯。

本研究的目的之一是分析油松居群遺傳多樣性及與生態(tài)因子間相關(guān)性。13條引物擴增的多態(tài)條帶比率非常高(100%)，在137個ISSR擴增譜帶中，特有條帶1條(山西三道川林場)，頻率≤25%差異條帶數(shù)量1—3條，頻率≤50%差異條帶數(shù)量7—24條，出現(xiàn)頻率≥5%的差異帶數(shù)76—125條，從而進一步說明作為分布于我國北方各省的廣布種，天然油松次生林具有極豐富的基因多態(tài)性。一般認為一個物種的遺傳多樣性愈豐富，對環(huán)境變化的適應(yīng)性就愈大［18］，油松天然居群間遺傳結(jié)構(gòu)的差異有利于油松提高抵御不良環(huán)境的能力，擴展其分布范圍，增強其適應(yīng)性。很多學(xué)者對不同植物居群進行了環(huán)境因子與遺傳多樣性指數(shù)間的相關(guān)性研究，發(fā)現(xiàn)居群遺傳多樣性與產(chǎn)地海拔、溫度、降雨量、土壤礦物質(zhì)含量等因子間存在一定的相關(guān)關(guān)系［16，19-24］，說明環(huán)境脅迫影響下，導(dǎo)致譜帶頻率和遺傳多樣性發(fā)生一定程度的梯度變化，引起相應(yīng)的適應(yīng)性生態(tài)地理分化。有大量研究［25-27］指出，植物抗寒性并非由單一基因控，是由多種特異的數(shù)量性抗寒基因調(diào)控的。本研究結(jié)果顯示(表6)，隨著產(chǎn)地年平均氣溫、極端最低溫的降低，天然油松林居群的觀測等位基因數(shù)量及多態(tài)位點百分率顯著降低，Shannon信息指數(shù)與降雨量呈現(xiàn)顯著負相關(guān)，從而說明低溫和降雨量影響了抗逆相關(guān)基因數(shù)量的表達，對天然油松次生林ISSR多態(tài)性產(chǎn)生選擇壓力。這與借助RAPD、等位酶等標記技術(shù)對其它植物的研究結(jié)果相似［19-20，22-23］?？傊瑢?dǎo)致油松天然居群間遺傳分化的原因除了與自然選擇及居群間缺乏有效的基因流有關(guān)外，還可能與生物學(xué)特性及其它原因有關(guān)，故遺傳分化—自然選擇—生物學(xué)特性三者間關(guān)聯(lián)性有待進一步研究探討。

3.2 油松天然居群遺傳多樣性地理變異模式

20世紀七八十年代，徐化成等對油松生長、表型性狀與地理-氣候因素關(guān)系進行了系統(tǒng)研究，認為生長特性與年均溫具有關(guān)聯(lián)性，油松在形態(tài)學(xué)水平的地理變異具有隨機性［1-5］。本研究從DNA水平研究得到不同結(jié)論，即產(chǎn)地氣溫、海拔和降雨量等產(chǎn)地氣候因子是影響天然油松次生林遺傳變異的重要環(huán)境因子，而不同經(jīng)緯度產(chǎn)地的本質(zhì)是其氣候條件的差異，與產(chǎn)地氣候因子相關(guān)也就意味著存在一定的地理模式，進而說明油松天然居群基因水平的地理變異并不是完全隨機的，具有一定的規(guī)律性。基于UPGMA聚類分析顯示，作為天然油松林分布中心的太行山脈、呂梁山脈及黃龍山脈的天然油松次生林相似性高，而分布于其它地域天然油松林具有一定的獨立遺傳特征。Shannon多樣性指數(shù)(I)和多態(tài)位點百分率(P)可評價遺傳多樣性的高低，指數(shù)越大，表明居群的遺傳多樣性水平越高。山西和順云龍公園、河南寶天曼及遼寧醫(yī)巫閭居群的Shannon信息指數(shù)及多態(tài)位點百分率均較高，并地處中國地勢二級和三級地勢階梯分界線附近，即油松天然居群分布區(qū)南北向中軸區(qū)域，故推測三居群與其它居群受季風影響種子及花粉交流相對頻繁，基因多態(tài)性水平高;而山東蒙山、甘肅冶力關(guān)居群相對獨立，說明該居群環(huán)境限制了與其它居群的交流，遺傳特征獨特，Shannon信息指數(shù)及多態(tài)位點比率均較低。由此進一步說明，油松天然居群遺傳變異存在一定的地理變異模式。

森林樹種應(yīng)對環(huán)境脅迫表現(xiàn)出明顯的表型適應(yīng)性，而關(guān)于基因區(qū)域的中性標記是否作為植物適應(yīng)性的選擇始終具有爭論，Neale和Kremer［28］認為，森林基因組學(xué)將為適應(yīng)環(huán)境脅迫的功能基因的研究提供重要工具，研究結(jié)果將為適應(yīng)性遺傳變異的中性選擇理論提供科學(xué)依據(jù)［29-30］。本研究結(jié)果借助隨機性ISSR標記，說明不同環(huán)境脅迫下油松天然次生林群體多態(tài)性差異和遺傳分化較大，反映油松天然居群一定的地理變異趨勢，雖不能準確反映與當?shù)丨h(huán)境適應(yīng)的生長、抗性、形態(tài)等指標功能基因的變異趨勢，但為中性標記的適應(yīng)性選擇理論提供了重要支持。

［1］ Xu H C.Chinese pine.Beijing:China Forestry Publishing House，1990.

［2］ Xu H C，Guo G R，F(xiàn)en L，Yang L.The relativity between natural Pinus tabulaeformis forest growth and geographical climatic factors.Journal of Beijing Forestry University，1981，4:9-14.

［3］ Xu H C，Sun Z F，Guo G R，F(xiàn)en L.Geographic distribution of Pinrs tabulaeformis carr.and classification of provenance regions.Scientia Silvae Sinicae，1981，17(3):258-270.

［4］ Xu H C，Tang J，Zhang S J，Ma H.A study on climatic ecotypes in pinus tabulaeformis Carr..Scientia Silvae Sinicae，1986，22(1):10-20.

［5］ Xu H C，Tang Q.A preliminary study on geographic variation in pinus tabulaeformis Carr..Journal of Beijing forestry college，1984，2:57-72.

［6］ Zhang C X，Li Y，Shen X H.Genetic diversity of ten natural and artificial populations of Pinus tabulaeformis in two isozyme systems.China's biodiversity conservation toward the 21st century.Beijing:China Forestry Publishing House，2000:137-144.

［7］ Li Y，Zhang C X.A Study on genetic diversity and breeding value of a series of clonal groups in Pinus tabulaeformis.Journal of Beijing Forestry University，1998，20(4):12-18.

［8］ Wang L，F(xiàn)an J F，Liu Y H，Yang P H，Wang M B.Genetic polymorphism of Pinus tabulaeformis from main distribution areas in China.Journal of Northeast Forestry University，2009，37(12):3-7.

［9］ Li C，Chai B F，Wang M B.Genetic diversity of Pinus tabulaeformis populations in North China.Bulletin of Botanical Research，2006，26(1):98-102.

［10］ Li C，Chai B F.Genetic Conservation of Pinus tabulaeformis Population Implicated by Molecular Markers.Chinese Journal of Applied ＆Environmental Biology，2009，15(4):459-463.

［11］ Zhang J J，Li M，Gao B J.Genetic diversity of natural Pinus tabulaeformis populations at different altitudes and in four communities types in Dawopu of Chengde，Hebei Province of North China.Chinese Journal of Ecology，2011，30(11):2421-2426.

［12］ Li M，Gao B J，Zhang J J.Genetic diversity of Pinus tabulaeformis populations at different altitudes in Guangtushan Mountains and relationship to environmental changes.Journal of Plant Genetic Resources，2012，13(3):350-356.

［13］ Wang M B，Gao F Q.Genetic Variation in Chinese Pine(Pinus tabulaeformis)，a Woody Species Endemic to China.Biochem Genet，2009，47:154-164.

［14］ Chen K M，Abbott R J，Milne R I，Tian X M，Liu J Q.Phylogeography of Pinus tabulaeformis Carr.(Pinaceae)，a dominant species of coniferous forest in northern China.Molecular Ecology，2008，17:4276-4288.

［15］ Eckert C G，Samis K E，Lougheed S C.Genetic variation across species'geographical ranges:the central-marginal hypothesis and beyond.Molecular Ecology，2008，17:1170-1188.

［16］ Nybom H.Comparison of different nuclear DNA markers for estimating intraspecific genetic diversity in plants.Molecular Ecology，2004，13:1143-1155.

［17］ Wright S.The genetic structure of populations.Annals of Eugenics，1951，15:323-354.

［18］ Ma W L.Higher plants and diversity.Beijing:Higher Education Press，1998.

［19］ Zhao N X，Gao Y B，Wang J L，Ren A Z，Ruan W B，Chen L，Liu H F.An analysis of genetic differentiation between Stipa grandis populations in middle and eastern Inner Mongolia steppe.Acta Ecologica Sinica，2004，24(10):2178-2185.

［20］ Li Z B，Lin P.Relationships between environmental factors and genetic variation of Aegiceras corniculatum populations.Acta Ecologica Sinica，2002，22(11):1912-1916.

［21］ Turpeinen T，Tenhola T，Manninen O，Nevo E，Nissil E.Microsatellite diversity associated with ecological factors in Hordeum spontaneum populations in Israel.Molecular Ecology，2001，10:1577-1591.

［22］ Xie Y Q，Li Z Z，Huang R Z，Xiao X X，Huang Y.Genetic diversity of betula luminifera populations at different altitudes in Wuyi Mountains.Scientia Silvae Sinicae，2008，44(3):50-55.

［23］ Xie Y Q，Huang R Z，Li Z Z，Huang Y，Yang Z W.Genetic variation of natural populations of Betula luminifera in Fujian and its relationship with the habitat.Scientia Silvae Sinicae，2009，45(9):60-65.

［24］ Zhou Z Q，Hao Y，Liu T，Jiang J.Corelativity analysis between the main ecological factors and genetic diversity of Pinus sylvestris var· mongolica population in the north part of Great Xing'an Mountains.Journal of Beijing Forestry University，2006，28(6):22-27.

［25］ Howe G T，Aitken S N，Neale D B，Jermstad K D，Wheeler N C，Chen T H H.From genotype to phenotype:unraveling the complexities of cold adaptation in forest trees.Canadian Journal of Botany，2003，81:1247-1266.

［26］ Wang N，Kang J M，Yang Q C，Xiong J B，Jin H C.Research progress in proteomics researches on plant cold-stress.Pratacultural scienc，2008，25(12):93-98.

［27］ Eckert A J，Heerwaarden J，Wegrzyn J L，Nelson C D，Ross-Ibarra J，Gonza'lez-Martínez C，Neale D B.Patterns of Population Structure and Environmental Associations to Aridity Across the Range of Loblolly Pine(Pinus taeda L Pinaceae).Genetics，2010，185(3):969-982.

［28］ Neale D B，Kremer A.Forest tree genomics:growing resources and applications.Nature Reviews Genetics，2011，12:111-122.

［29］ Grivet D，Sebastiani F，Gonzalez-Martinez SC，Vendramin GG.Patterns of polymorphism resulting from long-range colonization in the Mediterranean conifer Aleppo pine.New Phytologist，2009，184:1016-1028.

［30］ Eckert A J，Bower A D，Gonzalez-Martinez S C，Wegrzyn J L，Coop G，Neale D B.Back to nature:ecological genomics of loblolly pine(Pinus taeda Pinaceae).Molecular Ecology，2010，19:3789-3805.

參考文獻:

［1］ 徐化成.油松.北京:中國林業(yè)出版社，1990.

［2］ 徐化成，郭廣榮，馮林，楊雷.油松天然林的生長與地理-氣候因素的關(guān)系.北京林學(xué)院學(xué)報，1981，4:9-14.

［3］ 徐化成，孫肇鳳，郭廣榮，馮林.油松天然林的地理分布和種源區(qū)的劃分.林業(yè)科學(xué)，1981，17(3):258-270.

［4］ 徐化成，唐謙，張淑娟，馬虹.油松氣候生態(tài)型的研究.林業(yè)科學(xué)，1986，22(1):10-20.

［5］ 徐化成，唐謙.油松地理變異的初步研究.北京林學(xué)院學(xué)報，1984，2:57-72.

［6］ 張春曉，李悅，沈熙環(huán).油松10個天然與人工群體兩個同工酶系統(tǒng)的遺傳多樣性研究.面向21世紀的中國林木遺傳育種.北京:中國林業(yè)出版社，2000:137-144.

［7］ 李悅，張春曉.油松無性系群體育種值與遺傳多樣性研究.北京林業(yè)大學(xué)學(xué)報，1998，20(4):12-18.

［8］ 王磊，樊軍鋒，劉永紅，楊培華，王孟昌.我國油松主要分布區(qū)種質(zhì)資源遺傳多樣性.東北林業(yè)大學(xué)學(xué)報，2009，37(12):3-7.

［9］ 李毳，柴寶峰，王孟本.華北地區(qū)油松種群遺傳多樣性分析.植物研究，2006，26(1):98-102.

［10］ 李毳，柴寶峰.基于分子標記的油松種群遺傳保護分析.應(yīng)用與環(huán)境生物學(xué)報，2009，15(4):459-463.

［11］ 張靜潔，李明，高寶嘉.承德大窩鋪不同海拔與群落類型天然油松種群遺傳多樣性.生態(tài)學(xué)雜志，2011，30(11):2421-2426.

［12］ 李明，高寶嘉，張靜潔.承德光禿山不同海拔油松居群遺傳多樣性與生境因子關(guān)聯(lián)研究.植物遺傳資源學(xué)報，2012，13(3):350-356.

［18］ 馬煒梁.高等植物及其多樣性.高等教育出版社，北京，1998.

［19］ 趙念席，高玉葆，王金龍，任安芝，阮維斌，陳磊，劉惠芬.內(nèi)蒙古中東部草原大針茅的種群遺傳分化.生態(tài)學(xué)報，2004，24(10):2178-2185.

［20］ 黎中寶，林鵬.桐花樹種群遺傳變異與環(huán)境變量的關(guān)系.生態(tài)學(xué)報，2002，22(11):1912-1916.

［22］ 謝一青，李志真，黃儒珠，肖祥希，黃勇.武夷山不同海拔光皮樺種群遺傳多樣性及其與生態(tài)因子的相關(guān)性.林業(yè)科學(xué)，2008，44(3):50-55.

［23］ 謝一青，黃儒珠，李志真，黃勇，楊宗武.福建光皮樺野生種群遺傳變異及其與生境的關(guān)系.林業(yè)科學(xué)，2009，45(9):60-65.

［24］ 周志強，郝雨，劉彤，姜靜.大興安嶺北段天然樟子松林遺傳多樣性與主要生態(tài)因子的相關(guān)性研究.北京林業(yè)大學(xué)學(xué)報，2006，28(6):22-27.

［26］ 王楠，康俊梅，楊青川，熊軍波，金后聰.植物低溫脅迫蛋白質(zhì)組學(xué)研究進展.草業(yè)科學(xué)，2008，25(12):93-98.