林业科学  2011, Vol. 47 Issue (7): 91-96   PDF    
0

文章信息

王玉山, 邢世岩, 唐海霞, 冯殿齐
Wang Yushan, Xing Shiyan, Tang Haixia, Feng Dianqi
侧柏种源遗传多样性分析
Genetic Diversity of Platycladus orientalis Provenances
林业科学, 2011, 47(7): 91-96.
Scientia Silvae Sinicae, 2011, 47(7): 91-96.

文章历史

收稿日期:2010-09-15
修回日期:2011-02-11

作者相关文章

王玉山
邢世岩
唐海霞
冯殿齐

侧柏种源遗传多样性分析
王玉山1,2, 邢世岩1, 唐海霞1, 冯殿齐2    
1. 山东农业大学 泰安 271018;
2. 泰安市泰山林业科学研究院 泰安 271000
摘要: 应用AFLP技术对17省市的18个侧柏种源进行遗传多样性分析。选用的8个引物共扩增出多态性带1 613条,占92.91%;平均有效等位基因数1.199 3,平均Nei's基因多样性指数0.123 9,Shannon's信息指数0.194 9,揭示侧柏具有丰富的遗传多样性,其中中部种源遗传多样性低于南、北部种源。AMOVA分析表明侧柏种源遗传分化大,74.86%的遗传变异主要存在种源内,11.12%存在区域间,14.02%存在于区域内种源间,其分布的间断性、地理隔离以及低水平基因流(Nm = 1.437 2)是导致其种源间分化的主要因素。应用Nei(1972)遗传距离进行非加权组平均法(UPGMA)聚类分析,结果显示纬度相近的种源聚在一起,把18个种源划分为北部、中部、南部和山东4个种源区。Mantel检验也证实种源间的遗传距离与地理距离呈正相关。
关键词:侧柏    种源    AFLP    遗传多样性    遗传结构    
Genetic Diversity of Platycladus orientalis Provenances
Wang Yushan1,2, Xing Shiyan1 , Tang Haixia1, Feng Dianqi2    
1. Shandong Agricultural University Tai'an 271018;
2. Taishan Forestry Institute of Tai'an City Tai'an 271000
Abstract: The genetic diversity of eighteen Platycladus orientalis provenances from seventeen provinces was analyzed by AFLP markers. Eight primers produced a total of 1 613 polymorphic bands(92.19%). The average effective number of alleles, Nei's gene diversity and Shannon's information index were 1.199 3, 0.123 9 and 0.194 9, respectively, indicating that there was plentiful genetic diversity among the provenances of P. orientalis. Genetic diversity of the middle provenances was less than that of provenances collected from the south and the north regions. Analysis of molecular variance(AMOVA) showed that variance component within provenances (74.86%) was higher than that among regions (11.12%) and among provenances within regions(14.02%). The distribution of discontinuity, geographical isolation and low levels of gene flow(Nm = 1.437 2) between the provenances were major factors that led to genetic differentiation. Nei's (1972) genetic distance was applied to unweighted pair-group method with arithmetic means (UPGMA) cluster analysis, and the results showed that the similar latitude provenances were clustered together, and eighteen provenances were clustered into four groups (north, south, middle and Shandong provenances). A Mantel test confirmed that there was a positive correlation between the genetic distance and geographical distance.
Key words: Platycladus orientalis    provenances    AFLP    genetic diversity    genetic structure    

侧柏(Platycladus orientalis)原产中国,分布远及日本和欧美,国内大部分省市均有栽培(石文玉等,1986)。侧柏耐干旱瘠薄、适应性强,是我国重要的造林绿化树种(曲绪奎,1987)。侧柏在长期生长发育过程中,受生态环境和异花授粉的影响,种内发生很大变异(李继华,1985),形成不同的形态、生理和生态类型(苏霞等,2006)。近年来,许多学者针对侧柏种源地理变异规律(吴夏明,1986; 陈晓阳等,1990; 施行博等,1992; 杨传强等,2005)、生理生态特性(吴夏明等,1988; 李良厚,1989; 王华芳,1996; 戴建良等,1999)、核型分析(李林初等,1984; 马颖敏等,2009)等进行了大量的研究,但未见侧柏种源遗传多样性的研究报道,仅见有关柏科(Cupressaceae)植物遗传多样性的研究。张国盛等(2005)红雨等(2006)对不同生境叉子圆柏(Sabina vulgaris)种群的RAPD分析表明:大部分的遗传变异存在于群体之内,种群的遗传分化与生境有着一定的相关性。张仁波等(2007)采用RAPD分子标记研究了崖柏(Thuja sutchuenensis)的遗传多样性,表明崖柏种源内和种源间均存在着广泛的表型及基因型的变异,空间距离和生态因子对崖柏种源进化的影响较小。张茜等(2005)通过对祁连圆柏(Sabina przewalskii)叶绿体DNA trnT-trnF序列研究,认为瓶颈效应和奠基者效应造成了祁连圆柏居群的遗传多样性分布式样。

作为显性分子标记,AFLP具有多态位点丰富、稳定重复性好、DNA用量少,可以在预先不知道基因组背景情况下进行研究等特点(白瑞霞,2008)。本文利用AFLP技术,对侧柏不同地理种源的遗传多样性进行研究,以探讨影响种源遗传多样性和遗传结构的因素,为侧柏种源选择、遗传改良、种源区划和遗传多样性保护提供依据。

1 材料与方法 1.1 材料

供试的18个侧柏种源来自于山东平阴大寨山和淄博原山林场的全国侧柏种源试验林。两试验林于1984年造林,大寨山试验林采取4株小区、9次重复的随机区组设计,淄博原山试验林为12株小区、4次重复设计,造林密度为2 m × 3 m。各种源分布于23°30'—43°57' N,81°20'—120°15' E之间。4月份采集当年生新叶放入自封袋中,每个种源随机选8~10个单株采样(表 1)。

表 1 侧柏种源的采样情况 Tab.1 Samples of Platycladus orientalis provenances
1.2 研究方法 1.2.1 DNA提取

采用改良的CTAB法从硅胶干燥的幼叶中提取总DNA。

1.2.2 PCR扩增

利用北京鼎国生物技术公司的试剂盒及其操作指南进行。用筛选出的8对引物(表 2)进行选择性扩增。PCR扩增总体积为25 μL,含预扩增稀释样品2 μL,10 × PCR buffer 2.5 μL,dNTPs 0.5 μL,EcoR Ⅰ引物1 μL,Mse Ⅰ引物1 μL,Taq酶0.5 μL,H2O 17.5 μL。PCR反应在Gene Amp PCR System 9600 (Perkin Elmer,USA)中进行。扩增程序为: 94 ℃预变性30 s; 65 ℃退火30 s,72 ℃延伸80 s,12个循环,每轮循环温度递减0.7 ℃; 接着94 ℃ 30 s,55 ℃ 30 s,72 ℃ 80 s,23个循环。

表 2 AFLP选择性扩增引物产生的条带多态性 Tab.2 Polymorphism of AFLP bands obtained by selective amplification based on the eight primer-combinations
1.2.3 电泳与数据处理

扩增后的样品在ABI PRISM 377测序仪上进行电泳分离检测,得到AFLP的指纹图谱,生成由“1”和“0”组成的原始矩阵。计算多态性位点百分率(p) : p = (k/n) × 100%,式中: k为多态位点数目,n为所测位点总数。采用POPGENE 32计算等位基因平均数Na、有效等位基因数Ne、Nei基因多样性指数H、Shannon信息指数I,总遗传多样性Ht、群体内遗传多样性Hs、遗传分化系数Gst和基因流Nm。采用Nei (1972)遗传距离,按非加权组平均法(UPGMA)进行聚类分析。应用IBD 1.52(Bohonak,2002)进行Mantel检验,应用AMOVA1.55 (Excoffier et al., 1992; Huff et al., 1993)进行分子方差分析。

2 结果与分析 2.1 AFLP引物筛选与扩增结果

采用筛选“3 + 3”核心引物组合策略,从64对引物组合中,筛选出8对银染引物用于本试验,分别对供试材料的基因组DNA进行扩增,在侧柏18个种源176个样品中共获得1 613条多态性条带,平均每对引物扩增的总带数为202条,多态带百分率为92.91% (表 2)。其中E-AGG/M-CAG是最有效率的1对引物,其产生位点的多态性条带比率达到98.62%;而引物E-ACG/M-CAG扩增出的条带的多态性比率最低,其多态性条带比率为89.40%。

2.2 种源的遗传多样性

表 3看出: 18个侧柏种源的多态位点百分率26.24%~55.94%,平均为46.04%。平均有效等位基因数为1.199 3,平均Nei's基因多样性指数0.123 9,平均Shannon信息指数0.194 9,而其种级水平分别达到1.26,0.168 4和0.277 7,结果说明侧柏的遗传多样性丰富。各种源的遗传多样性存在较大差别,云南文山Nei's基因多样性指数值最大,为0.153 1,山东平度的遗传多样性较低,仅为0.063 5; Shannon信息指数最大值为宁夏银川的0.237 1,最小值是山东平度的0.102 5。这些数据表明侧柏种源的遗传多样性比较丰富,其中南部的云南文山、福建南平,北部的新疆伊犁、内蒙古乌拉山等种源的遗传多样性程度高,山东种源的多样性程度较低。

表 3 侧柏种源的遗传多样性水平 Tab.3 Genetic variation within 18 Platycladus orientalis provenances based on AFLP markers
2.3 遗传分化程度

遗传多样性分析结果显示: 18个种源总的遗传多样性Ht,群体内遗传多样性Hs和群体间遗传多样性Dst分别为0.168 4,0.123 9和0.045 5 (表 4)。Gst是衡量种群遗传变异最常用的指标,表示在总的遗传变异中群体间变异所占的比例,侧柏种源间的遗传分化系数Gst为0.262 2,说明种源间的遗传变异占总变异的26.22%,种源内变异占73.78%。AMOVA分析结果(表 5)表明:种源间的遗传变异占总变异的14.02%,而种源内变异占74.86%,区域间变异占11.12%,且其变异均达极显著水平(P<0.001),表明侧柏种源出现了较高程度的遗传分化,遗传多样性主要存在于种源内。

表 4 侧柏不同种源间的遗传分化(由Nei指数估计) Tab.4 Genetic differentiations of Platycladus orientalis provenances by Nei's index
表 5 侧柏种源的分子方差分析结果 Tab.5 Analysis of molecular variance of Platycladus orientalis provenances
2.4 种源间的遗传一致度和遗传距离分析

遗传一致度和遗传距离分别是从相同和相反方面度量居群间遗传关系的指标,遗传距离反映居群亲缘关系的远近。从表 6中可以看出:遗传一致度最大的来自江苏铜山和陕西黄陵之间(0.988 8),其遗传距离最小(0.011 3);而四川西昌和山东泰安间的遗传一致度最小(0.829 9),遗传距离最大(0.186 5)。18个侧柏种源间的遗传一致度介于0.829 9~0.988 8之间,平均为0.938 9;遗传距离为0.011 3~0.186 5,平均值0.063 9,说明种源间保持了较高的基因相似度,遗传距离小。

表 6 18个侧柏种源的遗传一致度与遗传距离 Tab.6 Nei's genetic identity and genetic distance of 18 Platycladus orientalis provenances
2.5 聚类分析

基于侧柏样本间的遗传距离Nei (1972)进行UPGMA聚类(图 1),将所有种源分为4大类,大致按地域聚在一起:第1类为南部种源,包括云南文山、福建南平、贵州贵阳、四川西昌、湖北荆门和河南确山6个种源; 第2类为北部种源,包括北京密云、内蒙古乌拉山、辽宁朝阳和新疆伊犁4个种源; 第3类为中部的6个种源,为甘肃徽县、江苏铜山、陕西黄陵、山西石楼、宁夏银川及河北易县; 第4类为山东种源,包括泰安和平度2个种源。聚类图显示纬度(地理距离)相近的种源被聚到了一起,反映了侧柏的遗传分化与生境有着一定的相关性,经Mantel检验,种源的地理距离与遗传距离之间显示正相关,但未达到显著性水平(r = 0.288 6,P = 0.093 0)。

图 1 18种源聚类图 Figure 1 Dendrogram of 18 Platycladus orientalis provenances 种源号见表 1。Provenance code see Tab. 1
3 结论与讨论 3.1 侧柏种源遗传多样性

本次研究所选的18个侧柏种源,是根据侧柏的生境(经纬度)不同选取的,8个引物检测到1 613条多态带,多态带比率达到92.91%,表明侧柏遗传多样性高。不过多态带百分率对遗传多样性只是一个粗略的估计,它不能够确定各条带在频率上的均匀程度(红雨等,2006),且受到样本大小、条带多少、筛选的引物及筛选方法的影响。

侧柏在物种水平的多态位点百分率P为46.04%,Nei基因多样性指数He为0.168 4,Shannon信息指数I为0.277 7 (表 3),与岷江柏木(Cupressus chengiana) (P = 77.17%,He = 0.225 7,I = 0.347 2) (姚莉,2005)、叉子圆柏(P = 77.17%,90.70%,He = 0.301 8,I = 0.453 4) (张国盛等,2005)和崖柏(P = 72.09%,He = 0.301 5,I = 0.436 0) (张仁波等,2007)等柏科树种相比,侧柏的遗传多样性较低。

侧柏南部种源(纬度<31°14')和北部种源(纬度>36°59')的遗传多样性高于中部种源和山东种源。因南、北种源一般分布于地形复杂的山区,交通闭塞,天然地理隔离明显,缺乏与外界基因的交流,且侧柏生境的人为破坏少,得以保存较高的遗传多样性。此结论说明:在南部和北部种源内进行优良种源选择时,选种数量宜多,适当缩小距离; 而中部种源应加大选种距离。

3.2 侧柏种源遗传结构和分化

Nybom(2004)总结了27份AFLP标记的植物遗传多样性研究结果,得出的Gst平均值为0.21; Hamrick等(1990)认为以异交为主的物种,90%的遗传变异发生在种源内部。侧柏种源间的基因分化系数(Gst)为0.262 2,与以上结果相比,侧柏种源呈现出种源间遗传分化大的特点。由Gst估算的侧柏种源基因流值为1.437 2,说明侧柏种源之间存在一定程度的基因交流,但Hamrick (1987)报道的4种针叶树的Nm值一般在5.3~37.8之间,异交植物Nm平均值为3.5(陈小勇,1996),相比之下侧柏种源间的基因流明显处于较低水平,低水平的基因流可能造成侧柏种源对局部生态环境的适应,进而促使种源间的遗传分化。以往研究(Hamrick et al., 1990; Bartish et al., 1999)表明种子扩散机制会影响种群间的分化水平,靠风力散布的植物具有低水平的种群分化(Gst = 0.143 0),而靠重力散布的植物具有较高的种群分化水平(Gst = 0.277 0)。侧柏主要依靠重力进行散播,其遗传分化系数为0.262 2,与以上结论一致,种子散播方式在一定程度上限制了基因的流动,从而增加了侧柏种源间的遗传分化。由于侧柏为雌雄同株异花风媒传粉植物,同一地区中相互授粉的机会多,基因交换频繁,种源的遗传基础趋于相似,因而位置相近种源具有较高的遗传一致性和较近的亲缘关系,可见地理隔离及传粉方式对侧柏种源遗传分化具有显著影响。另外,分布区内多样性的生态条件也是侧柏产生遗传分化的因素之一(邓朝经等,1986),不同的基因型在不同生境上的适合度相异,导致具相同基因型的个体聚集在较适应的微生境中,从而产生遗传分化。因此,侧柏现有分布的间断性、地理隔离、环境因子造成的选择差、低水平基因流可能是导致其种源间分化的主要因素。

侧柏AMOVA分子变异分析显示:总遗传变异的14.02%来自种源间,74.86%来自种源内,11.12%来自区域间,因侧柏群体数目大并且为雌雄同株、风媒异花授粉植物,所以遗传变异的大部分驻留在种源内。此遗传结构说明在不同种源间进行选种具有很大的潜力,因此利用不同种源优异的性状与丰富的遗传多样性进行侧柏遗传改良具有广阔的前景; 开展侧柏良种选育时,在注重种源选择的同时,要加大种源内个体的选择力度,以选育出遗传多样性高的种源。

3.3 侧柏种源间的亲缘关系

本研究基于Nei(1972)遗传距离,利用UPGMA法进行聚类,将18个种源聚为4大类,湖北荆门以南种源聚为一类,山西石楼以北的聚为一类,两者之间的(中部)种源聚为一类,山东种源单独聚为一类,说明侧柏种源间的亲缘关系与地理位置相关,Mantel检验结果证实侧柏种源间遗传距离和地理距离间具有正相关性。不同种源聚类结果与梁一池(1990)孙仲序等(1992)用生长性状进行侧柏种源区划分的结果基本一致。聚类结果为侧柏种源区划提供了分子基础,同时说明在调拨种子时,应从临近地区购进。为更好地进行侧柏良种选育、种源区划以及种子调拨区划分,建议在分子标记基础上,结合形态标记以及核型分析对侧柏遗传多样性进行综合评价,同时进一步利用分子标记结果进行侧柏遗传结构分析。

参考文献(References)
白瑞霞. 2008. 枣种质资源遗传多样性的分子评价及其核心种质的构建. 河北农业大学博士学位论文. http: //cdmd. cnki. com. cn/article/cdmd-11920-2008074989. htm
陈小勇. 1996. 植物的基因流及其在濒危植物保护中的作用[J]. 生物多样性, 4(2): 97-102.
陈晓阳, 陈振丙, 吴检柱, 等. 1990. 侧柏种苗性状地理变异趋势及其气候生态基础[J]. .河南林业科技, (1): 5-10.
戴建良, 董源, 陈晓阳, 等. 1999. 不同种源侧柏鳞叶解剖构造及其与抗旱性的关系[J]. 北京林业大学学报, 21(1): 26-31.
邓朝经, 马军, 董素华. 1986. 侧柏种源苗期性状变异的初步研究[J]. 四川林业科技, 7(4): 27-34.
红雨, 王林和, 梁小荣. 2006. 不同生境臭柏种群的遗传多样性分析及其与环境因子的相关性[J]. 干旱区资源与环境, 20(3): 184-187.
李继华. 1985. 山东柏树的栽培历史和现状[J]. 山东林业科技, (1): 73-74.
李良厚. 1989. 侧柏抗旱特性与形态特征的关系[J]. 河南林业科技, (4): 40-42.
李林初, 徐炳生. 1984. 侧柏和福建柏染色体核型的研究[J]. 云南植物研究, 6(4): 447-451.
梁一池. 1990. 侧柏种源苗期生长与原产地气候因子间的典范相关分析[J]. 林业科学, 26(4): 308-315.
马颖敏, 邢世岩, 王玉山, 等. 2009. 中国侧柏地理种源核型分析与进化趋势[J]. 分子植物育种, 7(6): 1186-1192.
曲绪奎. 1987. 国内侧柏种源苗期研究综述[J]. 山东林业科技, (3): 50-53.
施行博, 郑吉联, 曲绪奎. 1992. 侧柏地理变异的研究[J]. 林业科学研究, 5(4): 402-408.
石文玉, 郭蓓. 1986. 侧柏种源子代试验[J]. 北京农学院学报, 11(1): 106-109.
苏霞, 卢春光, 汤安生, 等. 2006. 关于山东省侧柏采种基地建立技术的探析[J]. 防护林科技, (增刊): 98-99.
孙仲序, 丰镇, 曲绪奎, 等. 1992. 侧柏种源研究报告[J]. 山东林业科技, (3): 1-6.
王华芳. 1996. 三个种源侧柏一年生苗木对磷的吸收[J]. 林业科学, 32(3): 219-224.
吴夏明, 马骥. 1988. 不同种源侧柏的光合特性和种子发芽对水分胁迫的反应[J]. 林业科学, 24(4): 448-453.
吴夏明. 1986. 侧柏的地理变异[J]. 北京林业大学学报, (3): 1-15.
杨传强, 丰震, 孙仲序, 等. 2005. 21年生侧柏种源的地理变异及种源选择[J]. 山东农业大学学报:自然科学版, 36(2): l96-l98.
姚莉. 2005. 中国特有种———濒危植物岷江柏(Cupressus chenggiana S. Y. Hu)的遗传多样性研究. 四川大学硕士学位论文. http: //cdmd. cnki. com. cn/Article/CDMD-10610-2005127422. htm
张茜, 杨瑞, 王钦, 等. 2005. 基于叶绿体DNAtrnT-trnF序列研究祁连圆柏的谱系地理学[J]. 植物分类学报, 43(6): 503-512.
张国盛, 张小红, 王林和, 等. 2005. 内蒙古地区天然臭柏种群遗传多样性的RAPD分析[J]. 干旱区资源与环境, 19(7): 193-198.
张仁波, 窦全丽, 何平, 等. 2007. 濒危植物崖柏遗传多样性研究[J]. 广西植物, 27(5): 687-691.
Bartish I V, Jeppsson N, Nybow H. 1999. Population genetic structure in the dioecious pioneer plant species Hippophae rhamnoides investigated random amplified polymorphic DNA (RAPD) by markers[J]. Molecular Ecology, 8(5): 791-802. DOI:10.1046/j.1365-294X.1999.00631.x
Bohonak A J. 2002. IBD (Isolation By Distance) : a program for analyses of isolation by distance[J]. Journal of Heredity, 93(2): 153-154. DOI:10.1093/jhered/93.2.153
Excoffier L, Smouse P E, Quattro J M. 1992. Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data[J]. Genetics, 131: 479-491.
Hamrick J L. 1987. Gene flow and distribution of genetic variation in plant populations//Urbanska K M. Differetiation Patterns in Higher Plants. London: Academic Press Inc. https://www.researchgate.net/publication/313526586_Gene_flow_distribution_of_genetic_variation_in_plant_populations
Hamrick J L, Godt M J W. 1990. Allozyme diversity in plant species//Brown A H D, Clegg M T, Kahler A L, et al. Plant population genetics, breeding, and genetic resources. Sinauer Associates, Inc., Sunderland, MA, 43-63. https://www.researchgate.net/publication/245264868_Allozyme_Diversity_in_Plant_Species
Huff D R, Peakall R, Smouse P E. 1993. RAPD variation within and among natural populations of outcrossing buffalograss [Buchloe dactyloides (Nutt.) Engelm.][J]. Theoretical and Applied Genetics, 86(8): 927-934.
Nei M. 1972. Genetic distance between populations[J]. American Naturalist, 106: 283-292. DOI:10.1086/282771
Nybom H. 2004. Comparison of different nuclear DNA markers for estimating intraspecific genetic diversity in plants[J]. Molecular Ecology, 13: 1143-1155. DOI:10.1111/mec.2004.13.issue-5