-
楸树(Catalpa bungei C. A. Mey.)原产中国,是紫葳科(Bignoniaceae)梓属(Catalpa Scop.)的高大落叶乔木,由于其对环境的适应性较好、生长速度迅速、木材材质优良,同时具用材和观赏两方面的特性,是理想的复合农林业建设树种[1]。金丝楸是楸树中的类型之一,在河南、山东、江苏、安徽等省有自然分布,由于它干型好、生长快、材性好、木材剖面色泽金黄,因此,在楸树适生区域备受青睐。河南是我国金丝楸的主产区之一,深入研究金丝楸优异种质资源,认识其遗传多样性水平、明确其亲缘关系,可以为金丝楸种质资源的保护和利用提供理论支持。
楸树种质资源研究始于20世纪70年代,早期主要是采用形态学的方法,姚庆渭等[2]研究了分布于亚洲东部和北美洲的梓树属不同种的形态学特征,建议设立梓树组、大梓树组和楸树组,其中,楸树组包括楸树和灰楸两个种。潘庆凯等[3]首次将楸树和灰楸2个种的种质资源从生产利用角度划分成了不同类型并进行了详尽的描述。河南省楸树研究组在1983年至1985年期间对楸树和灰楸资源开展了全国范围的调查,将我国的楸树种质资源在种以下划分成了11个类型,这是我国目前最为系统的楸树资源调查和分类研究[4]。分子标记技术的广泛应用使得楸树种质资源的研究已经深入到了分子水平。Li等[5]采用叶绿体ndhF序列和核糖体DNA 序列对紫葳科梓属的系统发育关系进行了探讨。石欣等[6]采集了156个楸树不同类型的单株材料,利用 ISSR分子标记分析发现,其遗传变异主要存在于类型内,且聚类分析发现,10个类型可划分为4类。Wang等[7]从楸树转录组中设计并筛选了30对EST-SSR引物,这些引物可应用于楸树及其近缘种的遗传分析。方乐成等[8]应用 SSR分子标记解析楸树种质资源遗传多样性和亲缘关系,构建了其核心种质。这些研究的开展为楸树种质资源的保存保护和开发利用提供了参考,但在分子水平上对其遗传多样性的研究相对较少,采用 AFLP技术进行金丝楸遗传多样性分析的研究还未见报道。
扩增片段长度多态性(AFLP)具有重复性好、可靠性强、灵敏度高等优点[9],是进行种质资源遗传多样性和遗传结构研究的有效分子标记方法,目前已在很多多年生木本植物中成熟应用。为了明确金丝楸优良无性系的遗传多样性,解析其与楸树种下不同类型及梓树属其它种之间的亲缘关系,本研究以75个金丝楸优良无性系单株和梓树属不同类型为研究对象,应用 AFLP技术分析它们的遗传多样性,建立DNA指纹图谱,揭示金丝楸优良无性系和楸树不同类型间的遗传关系,为金丝楸品种鉴定提供依据,同时也为确定金丝楸的分类地位、种质资源保护及可持续利用和新品种选育提供参考。
HTML
-
供试材料为定植于我国河南省洛宁县的75个金丝楸优良无性系、‘豫楸1号’、光叶楸、灰楸和梓树(图1、表1),采其健康、无病虫为害的嫩叶,低温运输、保存并立即提取其基因组DNA。
无性系号
Code母树来源地
Origin纬度(N)
Latitude经度(E)
Longitude无性系号
Code母树来源地
Origin纬度(N)
Latitude经度(E)
LongitudeJSQ-01 罗岭乡讲理村 34°23′40″ 111°22′46″ JSQ-41 长水乡连山村 34°24′33″ 111°24′07″ JSQ-02 罗岭乡讲理村 34°23′40″ 111°22′46″ JSQ-42 长水乡连山村 34°24′33″ 111°24′07″ JSQ-03 罗岭乡讲理村 34°23′40″ 111°22′46″ JSQ-43 长水乡连山村 34°24′33″ 111°24′07″ JSQ-04 罗岭乡讲理村 34°23′40″ 111°22′46″ JSQ-44 长水乡连山村 34°24′33″ 111°24′07″ JSQ-05 罗岭乡核桃沟村 34°23′59″ 111°21′51″ JSQ-45 长水乡连山村 34°24′33″ 111°24′07″ JSQ-06 罗岭乡核桃沟村 34°23′59″ 111°21′51″ JSQ-46 长水乡连山村 34°24′33″ 111°24′07″ JSQ-07 罗岭乡核桃沟村 34°23′59″ 111°21′51″ JSQ-47 长水乡连山村 34°24′33″ 111°24′07″ JSQ-08 罗岭乡核桃沟村 34°23′59″ 111°21′51″ JSQ-48 长水乡连山村 34°24′33″ 111°24′07″ JSQ-09 罗岭乡核桃沟村 34°23′59″ 111°21′51″ JSQ-49 长水乡连山村 34°24′33″ 111°24′07″ JSQ-10 罗岭乡核桃沟村 34°23′59″ 111°21′51″ JSQ-50 长水乡连山村 34°24′33″ 111°24′07″ JSQ-11 罗岭乡核桃沟村 34°23′59″ 111°21′51″ JSQ-51 长水乡连山村 34°24′33″ 111°24′07″ JSQ-12 罗岭乡贾村 34°24′36″ 111°20′03″ JSQ-52 长水乡连山村 34°24′33″ 111°24′07″ JSQ-13 罗岭乡贾村 34°24′36″ 111°20′03″ JSQ-53 长水乡连山村 34°24′33″ 111°24′07″ JSQ-14 罗岭乡贾村 34°24′36″ 111°20′03″ JSQ-54 长水乡连山村 34°24′33″ 111°24′07″ JSQ-15 罗岭乡贾村 34°24′36″ 111°20′03″ JSQ-55 长水乡连山村 34°24′33″ 111°24′07″ JSQ-16 罗岭乡贾村 34°24′36″ 111°20′03″ JSQ-56 长水乡连山村 34°24′33″ 111°24′07″ JSQ-17 罗岭乡贾村 34°24′36″ 111°20′03″ JSQ-57 长水乡连山村 34°24′33″ 111°24′07″ JSQ-18 罗岭乡贾村 34°24′36″ 111°20′03″ JSQ-58 长水乡连山村 34°24′33″ 111°24′07″ JSQ-19 罗岭乡贾村 34°24′36″ 111°20′03″ JSQ-59 兴华乡瓦寨庙村 34°15′41″ 111°27′52″ JSQ-20 罗岭乡后韩沟村 34°23′28″ 111°21′49″ JSQ-60 兴华乡瓦寨庙村 34°15′41″ 111°27′52″ JSQ-21 罗岭乡后韩沟村 34°23′28″ 111°21′49″ JSQ-61 兴华乡瓦寨庙村 34°15′41″ 111°27′52″ JSQ-22 罗岭乡后韩沟村 34°23′28″ 111°21′49″ JSQ-62 兴华乡瓦寨庙村 34°15′41″ 111°27′52″ JSQ-23 罗岭乡后韩沟村 34°23′28″ 111°21′49″ JSQ-63 兴华乡瓦寨庙村 34°15′41″ 111°27′52″ JSQ-24 罗岭乡后韩沟村 34°23′28″ 111°21′49″ JSQ-64 底张乡中沟村 34°17′44″ 111°28′10″ JSQ-25 罗岭乡后韩沟村 34°23′28″ 111°21′49″ JSQ-65 底张乡中沟村 34°17′44″ 111°28′10″ JSQ-26 罗岭乡后韩沟村 34°23′28″ 111°21′49″ JSQ-66 兴华乡烟洞村 34°17′27″ 111°27′28″ JSQ-27 罗岭乡窑沟村 34°23′22″ 111°21′05″ JSQ-67 兴华乡烟洞村 34°17′27″ 111°27′28″ JSQ-28 罗岭乡窑沟村 34°23′22″ 111°21′05″ JSQ-68 兴华乡北岭村 34°15′51″ 111°22′41″ JSQ-29 罗岭乡窑沟村 34°23′22″ 111°21′05″ JSQ-69 兴华乡北岭村 34°15′51″ 111°22′41″ JSQ-30 罗岭乡窑沟村 34°23′22″ 111°21′05″ JSQ-70 兴华乡北岭村 34°15′51″ 111°22′41″ JSQ-31 罗岭乡窑沟村 34°23′22″ 111°21′05″ JSQ-71 兴华乡北岭村 34°15′51″ 111°22′41″ JSQ-32 罗岭乡窑沟村 34°23′22″ 111°21′05″ JSQ-72 兴华乡北岭村 34°15′51″ 111°22′41″ JSQ-33 罗岭乡罗岭村 34°22′35″ 111°22′24″ JSQ-73 兴华乡北岭村 34°15′51″ 111°22′41″ JSQ-34 罗岭乡罗岭村 34°22′35″ 111°22′24″ JSQ-74 兴华乡北岭村 34°15′51″ 111°22′41″ JSQ-35 罗岭乡前阴坡村 34°22′15″ 111°23′51″ JSQ-76 兴华乡北岭村 34°15′51″ 111°22′41″ JSQ-36 罗岭乡前阴坡村 34°22′15″ 111°23′51″ Yuqiu1 南阳邓州市 32°44′56″ 112°12′49″ JSQ-37 罗岭乡前阴坡村 34°22′15″ 111°23′51″ GYQ 洛宁县城关 34°23′14″ 111°39′44″ JSQ-38 罗岭乡前阴坡村 34°22′15″ 111°23′51″ HuiQiu 洛宁县城关 34°23′14″ 111°39′44″ JSQ-39 罗岭乡前阴坡村 34°22′15″ 111°23′51″ ZiShu 洛宁县城关 34°23′14″ 111°39′44″ JSQ-40 罗岭乡前阴坡村 34°22′15″ 111°23′51″ 注:“JSQ-1”~“JSQ-76”——金丝楸优良单株采用根蘖繁殖产生的无性系;“Yuqiu1”——‘豫楸1号’;“GYQ”——光叶楸;“HuiQiu”——灰楸;“ZiShu”——梓树。
Notes: “JSQ-1”~“JSQ-76”- Superior clones of Jinsiqiu (Catalpa bungei C. A. Mey.) propagated from root suckers ; “Yuqiu1”- C. bungei var. 'Yuqiu 1'; “GYQ”- C. bungei Guangyeqiu; “HuiQiu”- C. fargesii Bur.; “ZiShu”- C. ovata G. Don.Table 1. Superior clones of Jinsiqiu used in the AFLP analysis and their origin
-
采用改良的CTAB法[10]提取其基因组DNA。
-
从64对EcoR Ⅰ+3/Mse Ⅰ+3引物组合(上海生工)中筛选出10对(表2)用于本研究。选取‘豫楸1号’、梓树等5个样品采用10对AFLP引物进行扩增,试验重复3次,验证其谱带质量和可重复性。
序号
No.引物组合
Primers总条带数/条 No. of total bands 多态性比例/% Percentage of polymorphic bands (PPBs) 总体 Total level 金丝楸 JSQ only 总体 Total level 金丝楸 JSQ only 1 E-AAC/M-CAA 78 71 94.87 83.09 2 E-AAC/M-CTT 99 83 91.92 73.49 3 E-AAG/M-CAA 149 128 94.63 75.78 4 E-ACT/M-CAA 139 125 97.84 76.80 5 E-ACT/M-CTA 143 124 95.80 72.58 6 E-ACT/M-CTG 127 116 99.21 88.79 7 E-ACT/M-CTT 109 105 97.25 62.86 8 E-GTA/M-CAC 97 86 91.75 75.58 9 E-TCC/M-CAG 125 110 96.80 76.36 10 E-TCC/M-CTT 111 93 100.00 93.55 合计 Total 1 177 1 041 − − 平均 Average 117.7 104.1 96.01 77.89 注:选择性引物的核心碱基序列: E-为5′–GACTGCGTACCAATTTC–3′; M-为5′–GATGAGTCCTGAGTAA–3′,每条引物的3′端包含3个选择性碱基,例如:E-AAC 包含核心碱基序列及其3′端的选择性碱基AAC。
Notes: Core sequence of selective primer combinations: E: 5′-GACTGCGTACCAATTTC-3′; M: 5′-GATGAGTCCTGAGTAA-3′. The 3′ end of each primer contains three selective bases, for example, E-AAC contains the core base sequence and its 3′ selective base AAC.Table 2. Primers used for AFLP and amplification results of Jinsiqiu
-
对供试样品采用EcoR Ⅰ/Mse Ⅰ内切酶(NEB公司)进行双酶切,酶切、连接、预扩增、选择性扩增程序参考李辛雷等[11]描述的方法进行。
-
将甲酰胺与分子量内标(100∶1)混合均匀后取9 μL进行上样,再将PCR产物稀释10×后吸取1 μL加入上样板,使用3730XL测序仪进行毛细管电泳检测。
-
基于采样地点的经纬度利用Google earth和Arcmap生成采样地点分布图。利用GeneMarker 2.2软件将79个材料10对荧光引物毛细管电泳的检测结果与分子量内标进行比较,得到片段大小,然后根据无带和有带的情况转换为“0/1”矩阵。采用SPSS 22.0软件的simple match模块检验AFLP试验的可重复性。用NTSYSpc 2.11F软件计算DICE遗传相似系数,并基于相似系数矩阵采用UPGMA法(非加权算术平均法)进行聚类分析,通过 Tree plot模块生成聚类图;用EIGEN模块进行基于Jaccard遗传相似系数的主坐标分析(PCA)[12],以进一步验证并分析聚类结果。使用POPGENE32 version 1.32软件[13]计算多态性比例PPBs (percentage of polymorphic bands),并在假定哈迪-温伯格平衡的前提下计算其有效等位基因数(Ne)、基因多样度(H)和Shannon信息指数(I)[14]。利用SPSS 22.0软件基于遗传距离和地理距离进行Mantel检验[15]。
1.1. 材料
1.2. 方法
1.2.1. DNA提取
1.2.2. 引物筛选和重复性验证
1.2.3. AFLP分析
1.2.4. 电泳检测
1.2.5. 数据分析
-
本研究筛选10对EcoR Ⅰ/Mse Ⅰ引物组合对金丝楸优良无性系基因组DNA进行AFLP分析,扩增取得了较好的效果(表2,图2),扩增谱带清晰、位置一致性良好(重复性>92%)。总体水平上,10对引物组合共扩增出1 177条带,其中,1 130条带是多态性带,平均多态性比例为96.01%;对75个金丝楸无性系而言,总的条带数是1 041条,其中,多态性带811条,平均的多态性比例是77.89%。其中,引物组合E-TCC/ M-CTT所检测的多态性比例最高为总体水平100%及金丝楸93.55%,总条带数分别为111条(总体)和93条(金丝楸);引物组合E-AAG/ M-CAA所获得的总条带数最多,为149条(总体)和128条(金丝楸),其多态性比例为94.63%(总体)和75.78%(金丝楸)。
-
利用10对选择性引物组合进行AFLP检测,79个样品的DICE遗传相似系数范围为0.413 4~0.905 9,平均为0.811 6。其中,JSQ-54和JSQ-76的遗传相似系数最大为0.905 9,JSQ-37和梓树间的遗传相似系数最小,为0.413 4。
-
基于等位基因频率计算有效等位基因数(Ne)、基因多样度(H)及Shannon信息指数(I)是进行遗传多样性水平评估的重要方式[12]。表3表明:10对引物所检测到的不同位点的有效等位基因数范围为1.00~1.99,平均为1.48±0.03(总体)和1.36±0.04(金丝楸);基因多样度为0.00~0.49,平均为0.28±0.01(总体)和0.21±0.02(金丝楸);Shannon信息指数为0.00~0.69,平均为0.43±0.02(总体)和0.32±0.04(金丝楸)。75个金丝楸无性系的有效等位基因数、基因多样度和Shannon信息指数高于对白榆无性系[16]、核用银杏[17]和古银杏资源[18]的研究结果(Ne=1.33,H= 0.20,I=0.33;Ne=1.33,H=0.21,I=0.34;Ne=1.26、H=0.17,I=0.27),低于利用ISSR分子标记对观赏石榴[19]的研究结果(Ne=1.94,H=0.38,I=0.55)和采用SSR标记对西藏光核桃[20]的研究结果,与利用ISSR技术对楸树种质资源[6]的分析(Ne=1.48,H= 0.29,I=0.45)及王庆军等[21]对对山东石榴品种(Ne=1.38,H=0.23,I=0.35)的研究较接近。
序号
No.引物组合
Primer有效等位基因数
Effective number of alleles (Ne)基因多样度
Nei's gene diversity (H)Shannon信息指数
Shannon's Information index (I)1 E-AAC/M-CAA 1.43±0.37 0.25±0.19 0.38±0.27 2 E-AAC/M-CTT 1.33±0.36 0.20±0.20 0.30±0.28 3 E-AAG/M-CAA 1.34±0.36 0.20±0.19 0.31±0.28 4 E-ACT/M-CAA 1.36±0.37 0.21±0.20 0.32±0.28 5 E-ACT/M-CTA 1.31±0.36 0.19±0.19 0.28±0.28 6 E-ACT/M-CTG 1.43±0.37 0.25±0.19 0.38±0.26 7 E-ACT/M-CTT 1.31±0.37 0.18±0.20 0.28±0.28 8 E-GTA/M-CAC 1.37±0.39 0.22±0.20 0.32±0.29 9 E-TCC/M-CAG 1.34±0.36 0.20±0.19 0.31±0.27 10 E-TCC/M-CTT 1.38±0.34 0.23±0.18 0.35±0.25 平均 Average 1.36±0.04 0.21±0.02 0.32±0.04 注:表中数值为平均值±标准差。
Note: The values in the table are presented as mean ± SD.Table 3. Genetic diversity of 75 superior clones of Jinsiqiu
-
分析聚类图(图3)发现:当阈值设定为0.46时,可将79个供试样品分成2个大组(A和B):A组中仅有梓树(ZiShu),B组包含了其他78个单株。分析B组成员发现:75个金丝楸无性系首先聚在一起,然后与光叶楸(GYQ)和灰楸(HuiQiu)聚在一起,‘豫楸1号’则单独聚为一组。
当以DICE遗传相似系数平均值(0.8277)划分时,可将79个供试材料分为8个组(分别标记为Ⅰ~Ⅷ),其中,Ⅰ组成员最多,包含71个金丝楸无性系,其遗传相似系数为0.806 6~0.905 9,JSQ-30和JSQ-58间的遗传相似系数最小为0.806 6,JSQ-54和JSQ-76间的遗传相似系数最大为0.905 9,这些无性系的亲缘关系较密切。Ⅲ组包含2个无性系,即JSQ-7和JSQ-8(DICE遗传相似系数0.837 3),它们的母树均采自罗岭乡核桃沟村,亲缘关系较近。JSQ-37和JSQ-63分别构成了单独的组(Ⅱ、Ⅳ)。JSQ-37的母树采自罗岭乡前阴坡村,JSQ-63的母树采自兴华乡瓦寨庙村,这2个金丝楸无性系表现出了与其它无性系较远的亲缘关系。GYQ、HuiQiu、Yuqiu1和ZiShu作为种内和种间对照各自聚为一组(Ⅴ、Ⅵ、Ⅶ和Ⅷ),与金丝楸亲缘关系较远。GYQ是楸树中不同于金丝楸的一个类型[4],它单独聚为一组。值得注意的是,在本研究中与近缘种灰楸相比,普通楸树品种‘豫楸1号’与金丝楸的亲缘关系更远。
通过NTSYSpc 2.11F软件EIGEN程序进行的PCA分析表明:前2个主坐标Dim-1和Dim-2所能解释的相关性为84.94%,PCA分析结果与聚类结果基本一致(图4),这2种分析相结合可以更全面地反映金丝楸优良无性系的遗传多样性和亲缘关系。进一步对各无性系的遗传距离和地理距离进行相关分析,结果呈极显著正相关(图5),表明地理隔离对楸树遗传多样性有显著的影响。结合聚类分析结果,金丝楸无性系当中的63号与其他遗传距离最远,值得重点推广。