怒江州深紋核桃種質資源SSR遺傳多樣性分析

吳 濤, 陳少瑜,3, 寧德魯, 肖良俊, 賀 娜, 唐春云, 和玉德

(1.云南省林業科學院，云南 昆明650201；2.云南省森林植物培育與開發利用重點實驗室，云南 昆明650201；3.云南省木本油料工程技術研究中心，云南 昆明 650201；4.怒江州木本油料產業發展研究所， 云南 瀘水 673199；5.怒江州林業科學研究所，云南 瀘水 673199)

我國現有核桃屬(Juglans)植物8種，作為堅果栽培的有2種，包括核桃(Juglansregia)和深紋核桃(泡核桃、云南核桃)(Juglanssigillata).核桃適應性較強、分布廣泛，主要栽培于山西、河北、新疆等北方地區，是核桃屬中栽培最廣的樹種.深紋核桃是我國的特有種，主要分布在我國西南地區(云南、貴州、四川西部、西藏南部最為集中)[1]，也是云南的主要栽培種.怒江傈僳族自治州地處橫斷山三江并流核心區，屬青藏高原和橫斷山脈的過度地帶.特殊的自然地理環境和立體氣候，為動植物的生長提供了多樣的大環境與異常豐富的小環境，是全球生物資源最富集的地區之一[2].作為深紋核桃的自然分布區，該區域也蘊育了豐富的核桃種質資源，初步調查發現該區域現保存有一定規模的深紋核桃天然群體和數量可觀的實生群體，這些豐富和多樣的資源將為核桃育種以及資源開發和利用提供物質基礎.但對現有核桃種質資源遺傳背景、親緣關系缺乏深入系統的了解使得進一步的育種工作受到很大限制.對核桃種質資源，尤其是野生資源進行遺傳多樣性分析和評價，可為核桃資源的保護利用提供可靠的理論依據，也可為核桃遺傳育種及遺傳規律的研究奠定基礎.

DNA分子標記技術是當前植物資源遺傳多樣性評價的強有力工具，尤其是基于PCR基礎的分子標記已廣泛應用在植物遺傳多樣性分析[3].近年來，國內外采用隨機擴增多態性DNA(random amplification polymorphism DNA, RAPD)，簡單重復序列間擴增(inter-simple sequence repeats, ISSR)，簡單序列重復長度多態性(simple sequence repeats, SSR)，擴增片段長度多態性(amplified fragment length polymorphis, AFLP)等分子標記技術對核桃種質資源進行遺傳多樣性研究[4-9],為核桃種質資源利用及保存提供技術支撐.SSR分子標記技術是以PCR為基礎的DNA多態性檢測技術，具有穩定、共顯性和信息量大等優點，被廣泛地應用于核桃種質鑒定和居群遺傳特點的研究中[10].本課題組在深紋核桃SSR標記開發及云南省核桃種質資源調查基礎上，對怒江州豐富的核桃種質資源進行遺傳多樣性研究，揭示其遺傳多樣性水平及遺傳結構狀況，為資源的利用和保護提供科學依據.

1 材料與方法

1.1 研究區域概況及材料來源

怒江傈僳族自治州位于云南省西北部(25°33′—28°23′N, 98°09′—99°39′E)，怒江中游.總面積14 703 km2，州內最低海拔738 m，最高海拔5 128 m，年均溫度11～20 ℃，年均降雨量1 000～1 700 mm.怒江州屬低緯高原季風氣候，同時受地形地貌和緯度差異的影響，形成獨特的立體氣候.該區下轄貢山、福貢、瀘水和蘭坪4個縣，4個縣均有豐富的核桃資源分布.

根據當地林業部門的相關資料及走訪當地群眾獲取的深紋核桃的資源信息，確定用以分析的資源.共采集122份深紋核桃資源樣品，每份資源相距至少100 m,分別采自怒江州的4個縣，4個核桃群體地理信息、樣本數列于表1.每份資源以采集的單株(二倍體)穗條作為1個樣本，穗條采回后取韌皮提取DNA.

表1 地理信息及樣本容量Table 1 Geographical information of sampling sites and sample sizes

1.2 試驗方法

1.2.1 DNA提取 采用改進CTAB法提取核桃穗條韌皮基因組DNA[11]，提取的DNA用瓊脂糖凝膠電泳結合紫外可見分光光度計(Thermo Scientific)檢測濃度與純度，用Tris-EDTA緩沖液將樣品稀釋至20 ng·μL-1，-20 ℃保存備用.

1.2.2 SSR-PCR擴增及毛細管電泳 SSR-PCR反應體系及擴增程序參照陳少瑜等方法[12]，用20對多態性高且穩定性好的SSR引物進行核桃種質資源SSR-PCR擴增，20對引物信息見表2.毛細管電泳參照Chen et al的法[10].

1.3 數據統計及分析

用GenAIEx 6.5軟件計算等位基因數，用POPGENE version 1.32軟件分析群體遺傳多樣性參數[13]，包括多態位點百分率(percentage of polymorphic bands, PPB)、觀察等位基因數(observed number of alleles,Na)、有效等位基因數(effective number of alleles,Ne)、觀察雜合度(obseved heterozygosity,Ho)、期望雜合度(expected heterozygosity,He)、Shannon信息指數(Shannon′s information index,I)、Nei′s遺傳多樣性指數(Nei′s gene diversity,H)、群體間的遺傳距離(genetic distance,D)和遺傳一致度(genetic identity,I).

表2 20對SSR引物信息及擴增片段Table 2 Profiles of 20 SSR primers and amplified fragments

用POPGENE version 1.32軟件獲得群體遺傳分化狀況參數，包括群體內繁育系數(inbreeding coefficient at the population level,FIS)、總群體內繁育系數(inbreeding coefficient at the total sample level,FIT)、群體分化系數(proportion of differentiation among populations,FST)[14]以及反映基因流強度的群體每代遷移數(gene flow,Nm)[15].基于4個群體的遺傳一致度，運用Mega 6.06軟件進行UPGMA聚類并繪制各群體遺傳關系樹狀圖[16].用Mantel檢驗估算群體間遺傳距離和地理距離的相關性[17].

2 結果與分析

2.1 怒江州核桃種質資源的遺傳多樣性

將篩選出的穩定性強、多態性高的20對SSR引物用于怒江州4個群體122份核桃種質資源的SSR擴增，擴增片段大小在96～250 bp之間.共檢測到等位基因184個，平均為10.82個.擴增結果經POPGENE version 1.32軟件分析得到4個群體的遺傳多樣性參數(表3).

表3 怒江州4個核桃群體的遺傳多樣性參數1)Table 3 Genetic diversity parameters of 4 J.sigillata populations

1)Na：觀察等位基因數；Ne：有效等位基因數；Ho：觀察雜合度；He：期望雜合度；I：Shannon信息指數；H：Nei′基因多樣性指數.

由表3可得，4個群體的多樣性參數存在明顯差異，其中蘭坪群體的有效等位基因數、期望雜合度、Shannon信息指數和Nei′基因多樣性指數在所有群體中最高，分別是3.42、0.650、1.365和0.641；而瀘水群體的各參數最低，分別是2.41、0.503、0.937和0.485.因此，怒江州4個核桃群體中蘭坪群體的遺傳多樣性相對較高，瀘水群體的遺傳多樣性相對較低.

Nei′s基因多樣性指數及Shannon′s信息指數可以有效反映群體的遺傳多樣性[18]，期望雜合度是反映種質資源遺傳多樣性的重要參考指標[19]，其中期望雜合度不僅是衡量群體遺傳多樣性最常用的指標，同時也反映群體中等位基因的豐富度和均勻程度[7].本研究結果Nei′基因多樣性指數平均為0.576，Shannon信息指數平均為1.157，期望雜合度平均為0.591，與相關采用SSR分子標記的核桃相比(表4)，怒江州深紋核桃種質資源的遺傳多樣性屬于中等水平.

表4 核桃屬植物SSR標記的遺傳多樣性參數比較Table 4 Comparison on genetic parameters of walnut populations based on SSR

2.2 怒江州核桃種質資源的遺傳分化

怒江州核桃群體間遺傳分化狀況詳見表5.由表可知，群體間的近交系數介于-0.109 6～0.690 5之間，平均為0.244 6，其中3個位點(JSI-63，JSI-15和ZMZ11)的雜合子過剩.群體間的遺傳分化系數介于0.027 1～0.126 4之間，平均為0.060 9，表明遺傳變異主要存在于群體內，只有6.09%的遺傳變異來源于群體間.各位點的基因流(Nm)差異較大，變化范圍在1.727 2～8.992 0之間，平均為3.858 3.

2.3 群體的遺傳關系及聚類分析

4個核桃群體間Nei′s遺傳距離D值的變化范圍為0.083 6～0.164 6(表6).其中蘭坪群體和福貢群體的遺傳距離最小(0.083 6)；貢山群體和瀘水群體的遺傳距離最大(0.164 6).根據Nei′s遺傳一致度，構建群體遺傳關系UPGMA聚類圖(圖1).從聚類圖可以看出，當遺傳距離約為0.1時，4個群體聚成了3組，其中蘭坪和福貢聚為一組，貢山和瀘水分別為一組.

現階段我們國家對于稅源的管理不夠完善，征管措施比較落后，公民納稅意識比較薄弱，所以在政所的時候需要循序漸進，按照合理的征收過程，把現階段的分類稅制過渡到綜合分類稅制的形式，而且需要把綜合所得課稅為基礎，以分類所得稅當作輔助內容。進而將連續性較強以及經營性的收入增加到綜合所得的征收項目里面，而且需要選擇累進稅率，按照所得比例進行征收。

利用Mantel檢驗計算種群的地理距離和遺傳距離相關性，當P<0.05時，認為地理距離和遺傳距離具有顯著相關性.4個核桃群體Nei′s遺傳距離與地理距離的Mantel檢驗結果顯示群體間的遺傳距離與地理距離之間呈現顯著相關性(r=0.702 9，P=0.032 9).

3 小結與討論

3.1 群體的遺傳多樣性

研究表明，怒江州深紋核桃種質資源的遺傳多樣性屬于中等水平，其中蘭坪群體的遺傳多樣性較高，瀘水群體的遺傳多樣性較低.

表5 怒江州核桃種質資源的遺傳分化與基因流1)Table 5 Genetic differentiation and gene flow of 4 J.sigillata populations in Nujiang Prefecture

1)FIS：群體間的近交系數；FIT：總的群體間近交系數；FST：群體間的遺傳分化系數；Nm：基因流.

表6 怒江州4個核桃群體間Nei′s遺傳距離(D)與遺傳一致度(I)1)Table 6 Nei′s genetic distances and genetic identity among 4 populations

1)遺傳一致度(I)(對角線之上)，遺傳距離(D)(對角線之下).

圖1 基于Nei′s遺傳距離的UPGMA聚類圖Fig.1 Phylogenetic dendrogram based on genetic distance by UPGMA cluster

準確估算一個樹種在其分布區的遺傳多樣性對于基因保存和優樹選擇非常重要[23]，而某一樹種的遺傳多樣性水平及分布格局是地理分布、繁育系統、人為影響等諸多因素共同作用的結果.影響遺傳多樣性的因素有繁育系統、遺傳漂變、基因突變和基因流等內部因素，還包括由于環境變化和人為干擾引起的種群隔離、生境片斷化等外部因素[24].本研究中，怒江州復雜的地形地貌和多樣氣候環境是維持該區域核桃種質資源多樣性的重要原因之一，其次核桃為同株異花、雌雄異熟的風媒傳粉樹種，長期的自然演變和相互授粉，導致核桃群體高度雜合和較為復雜的遺傳背景，使核桃群體產生較高的遺傳多樣性.然而，另一方面，風媒傳粉距離的限制，地理分布和人為活動導致的片段化又會導致一定水平的遺傳漂變[8]；同時，核桃作為重要的經濟林樹種之一，人為干擾較大，諸如種質的交換、良種的種植推廣、野生資源及其生境的破壞等致使多樣性降低，這些因素綜合作用的結果使得此區域核桃資源的遺傳多樣性僅處于中等水平.其中蘭坪群體具有相對較高的多樣性是由于收集的樣品數量較多，另外調查也發現蘭坪的野生資源較為豐富.

3.2 群體間的遺傳分化

4個群體的遺傳分化系數FST=0.060 9，表明僅有6.09%的遺傳變異來源于群體間，遺傳變異主要來自群體內，群體內的分化是該區域核桃種質資源變異的主要途徑.此分化系數高于基于SSR標記的黑核桃居群(FST=0.017)[22]，低于西藏核桃居群(FST=0.103)和泡核桃居群(FST=0.111)[21]，也低于中國8個核桃居群(FST=0.196)[7]，說明該區域核桃群體間的遺傳分化水平較低.

基因流是影響群體內部和群體之間遺傳變異程度的重要因素[25]，而群體的遺傳分化狀況又直接反映了群體的遺傳結構.基因流的大小可以反映群體遺傳分化的狀況.Hamrick et al[26]認為當Nm>1時，基因流足以抵制遺傳漂變的作用，也可以防止群體間的分化；當Nm<1時，遺傳漂變可導致明顯的群體間遺傳分化.怒江州4個核桃群體的平均Nm為3.858 3(Nm>1)，說明群體間的基因交流阻止了群體間的分化，這與陳良華等對四川核桃資源的遺傳多樣性研究結果相似[27].Nm水平主要與物種的繁殖方式、花粉和種子傳播方式、地理分布有關，也與環境變化或人為活動密切相關.有研究表明，具有高度異交率的樹種Nm值較大，如冷杉、桉樹(Nm>5)[28].盡管核桃通過雌雄花的異熟性避免自交又是典型的風媒傳粉，但研究表明核桃存在一定的自交現象[29]，風媒花粉傳播距離有限，尤其對于分布于怒江山谷溝壑中的野生資源花粉的傳播受到地理隔離，因而核桃群體間的基因交流更可能是種子傳播造成的.人為馴化也是群體遺傳分化低下的重要原因[30].總之，怒江州為傳統的農業區，核桃又是具有特殊經濟價值的樹種，受人為影響較大，人為實生選種、種子的交換加強了群體間的基因交流，降低了群體間的遺傳分化.

基于遺傳距離的UPGMA聚類表明，聚類結果與群體的地理位置和距離有顯著關系.首先，同屬怒江中游，而且地理距離相對近的蘭坪和福貢群體聚為一組；然后與怒江下游地理距離相對較遠的瀘水群體聚類，最后再與怒江上游、地理距離最遠的貢山聚類.蘭坪、福貢和瀘水位于怒江中下游，三者與上游的貢山間分布有怒江大峽谷，一定程度上阻隔了基因流動.4個群體的地理分布特點導致了其特有的遺傳結構和聚類結果.

3.3 怒江州核桃種質資源的保護及利用

遺傳多樣性是資源利用的基礎，應重視對它的保護和利用.4個群體中，蘭坪群體的遺傳多樣性最高，其野生資源也相對豐富.野生資源是遺傳多樣性的重要來源，還可能蘊藏著珍貴的基因資源，應加強系統研究和科學保護；其次在根據怒江州深紋核桃資源群體遺傳分化特征進行選育和保護工作時，應在選擇遺傳多樣性高的群體基礎上，側重于群體內家系和單株的選擇，以盡可能地保存和利用好該區域的多樣性.