留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

抗松材线虫病马尾松种质资源遗传多样性分析及核心种质构建

邓莉丽 刘青华 周志春 高凯 骆定会

邓莉丽, 刘青华, 周志春, 高凯, 骆定会. 抗松材线虫病马尾松种质资源遗传多样性分析及核心种质构建[J]. 浙江农林大学学报, 2024, 41(1): 67-78. doi: 10.11833/j.issn.2095-0756.20230333
引用本文: 邓莉丽, 刘青华, 周志春, 高凯, 骆定会. 抗松材线虫病马尾松种质资源遗传多样性分析及核心种质构建[J]. 浙江农林大学学报, 2024, 41(1): 67-78. doi: 10.11833/j.issn.2095-0756.20230333
DENG Lili, LIU Qinghua, ZHOU Zhichun, GAO Kai, LUO Dinghui. Genetic diversity analysis and core collection of pinewood nematodiasis-resistant Pinus massoniana germplasm resources[J]. Journal of Zhejiang A&F University, 2024, 41(1): 67-78. doi: 10.11833/j.issn.2095-0756.20230333
Citation: DENG Lili, LIU Qinghua, ZHOU Zhichun, GAO Kai, LUO Dinghui. Genetic diversity analysis and core collection of pinewood nematodiasis-resistant Pinus massoniana germplasm resources[J]. Journal of Zhejiang A&F University, 2024, 41(1): 67-78. doi: 10.11833/j.issn.2095-0756.20230333

抗松材线虫病马尾松种质资源遗传多样性分析及核心种质构建

doi: 10.11833/j.issn.2095-0756.20230333
基金项目: 浙江省“十四五”育种专项林木协作组课题(2020C02007);江西省林业局林业科技创新专项(创新专项﹝2021﹞13号)
详细信息
    作者简介: 邓莉丽(ORCID: 0009-0002-3535-9972),从事林木遗传育种研究。E-mail: Lili_Deng1706@163.com
    通信作者: 刘青华(ORCID: 0000-0002-7960-6739),研究员,博士,从事抗性马尾松遗传改良研究。E-mail: liuqinghua642@163.com
  • 中图分类号: S722.3

Genetic diversity analysis and core collection of pinewood nematodiasis-resistant Pinus massoniana germplasm resources

  • 摘要:   目的  对来自安徽省林业科学研究院和异地保存在浙江省临海市林业技术推广和场圃旅游服务总站的114份抗松材线虫Bursaphelenchus xylophilus病马尾松Pinus massoniana种质进行遗传多样性与群体结构分析,构建抗性马尾松核心种质库。  方法  对114份抗性马尾松种质进行检测,使用分子生物学软件计算遗传多样性参数,并进行主坐标(PCoA)和群体结构分析。利用M策略和随机取样策略分别构建核心种质,分析不同核心种质的的遗传多样性参数,确定最适合的构建方法。  结果  114份抗性马尾松种质共检测到115个等位基因,平均有效等位基因数(Ne)为5.54,平均Shannon’s多样性指数(I)为1.51,平均多态信息含量(PIC)为0.90,结果表明其具有较高的遗传多样性水平。群体结构分析将114份抗性马尾松种质分为4个亚群,主坐标分析结果与上述基本一致。根据遗传多样性参数和抽样数量综合考虑,M策略构建的核心种质能以最小的种质数保留原有种质最大的遗传多样性,为最佳的取样策略。利用该策略得到了72份核心种质,其保留了原有种质100%的等位基因数,NeI、期望杂合度(He)、PIC等遗传参数的保留率分别为95.67%、94.96%、98.12%和100.00%,将构建的核心种质与原有种质进行t检验、PCoA分析和UPGMA聚类分析,结果表明两者间的遗传多样性无显著差异。  结论  114份抗性马尾松种质遗传多样性水平较高,构建的72份抗性马尾松核心种质,去除了遗传冗余,有利于抗性马尾松种质资源的有效保护和科学利用,可为优异基因发掘和新品种选育提供参考。图7表8参37
  • 图  1  引物多态性检测图

    Figure  1  Primer polymorphism detection diagram

    图  2  lnP(D)和∆KK变化的折线图

    Figure  2  Line plots of lnP(D) and ∆K as a function of K

    图  3  114份抗性马尾松样本的Structure聚类图

    Figure  3  Structure cluster diagram of 114 samples of resistant P. massoniana

    图  4  抗性马尾松114份样本主坐标分析

    Figure  4  Principal coordinate analysis of 114 samples

    图  5  72份核心种质与114份原有种质主坐标分析     

    Figure  5  Principal coordinate analysis of 72 core germplasm and 114 original germplasm

    图  6  72份抗性马尾松核心种质的UPGMA聚类图

    Figure  6  UPGMA clustering map of 72 samples of resistant P. massoniana

    表  1  材料编号及来源

    Table  1.   Test materials’ number and source

    编号来源编号来源编号来源编号来源
    1广西贵港35安徽黄山DAC12浙江临海HGS5浙江临海
    2广西桂平36安徽黄山DAC13浙江临海HGS6浙江临海
    3贵州黎平37安徽全椒DAC15浙江临海HGS8浙江临海
    4福建南平38安徽泾县DAC17浙江临海HGS9浙江临海
    6福建南平39安徽休宁DAC18浙江临海HGY1浙江临海
    7安徽休宁40安徽休宁DAC19浙江临海LS001浙江临海
    9安徽休宁41安徽休宁DAC20浙江临海LS003浙江临海
    10安徽祁县42安徽休宁DAC21浙江临海LS005浙江临海
    11福建南平43安徽休宁DAC22浙江临海LS007浙江临海
    12浙江淳安44安徽休宁DAC23浙江临海LS008浙江临海
    13浙江淳安45安徽休宁DAC24浙江临海LS009浙江临海
    14福建南平46安徽休宁DAC26浙江临海LS010浙江临海
    15浙江松阳47安徽休宁DAC34浙江临海XJ003江西峡江
    17福建南平48安徽休宁DLS27浙江临海XJ004江西峡江
    19江西婺源49安徽休宁DLS28浙江临海XJ005江西峡江
    20福建南平71江西分宜DLS29浙江临海XJ014江西峡江
    21安徽广德114江西分宜DLS30浙江临海AH1安徽广德
    22安徽广德161江西分宜DLS31浙江临海AH2安徽滁州
    23安徽广德231江西分宜DLS32浙江临海AH3安徽黄山
    24安徽广德242江西分宜DLS33浙江临海AH4安徽休宁
    25安徽广德305江西分宜DLS35浙江临海AH5安徽广德
    26安徽广德355江西分宜DLS36浙江临海AH6安徽广德
    27安徽广德398江西分宜DLS37浙江临海AH7安徽广德
    28安徽广德413江西分宜DLS38浙江临海AH8安徽泾县
    29安徽广德461江西分宜DLS40浙江临海AH9安徽广德
    30安徽和县473江西分宜HGS11浙江临海AH10安徽休宁
    31安徽和县485江西分宜HGS2浙江临海AH11安徽广德
    32安徽和县486江西分宜HGS3浙江临海
    34安徽和县490江西分宜HGS4浙江临海
      说明:AH1~AH11为 安徽省林业科学研究院审定的马尾松抗松材线虫病良种皖马抗1~11号。
    下载: 导出CSV

    表  2  14对SSR引物信息

    Table  2.   14 pairs of SSR primer information

    编号引物名称正向引物(5´→3´)反向引物(5´→3´)
    1 H03 CTCCAAAGGCGAGACTGC ACGAAAGCCAAGCTGAAC
    2 D04 AGGATGGTATGGTCGTGG CCCTTCTTCGCTCTGTGA
    3 H04 AAGAAGCGATCTGAGATGACTAA CTGTCATTGATTGTTTCCTTTTG
    4 B05 CCGTGCCTTCAGCATCTTCT CAGTGGATCTGTCACCTCCTCAT
    5 F05 AGAAGAAGAGCAGCAGTTTC GGTTTTCCATTGTTCTCACT
    6 H05 GTGCCTTCAGCATCTTCTAC ATCTGTCACCTCCTCATCTT
    7 B06 TATTAGCACCCCTCCAAAG TGGGGTAGAAGAATCGTAAGT
    8 D06 GATTCGGCTTCGTGACCTT CCCCCATAACCCCTGTC
    9 F06 AAGGACTTACAGAGGTTGGGTT GCTGCGACGAGCGTTTCT
    10 H06 TTCCTACCGCTGGGTTCTTG GACCTGACCTCGGGCATTAC
    11 D07 TCGCCTGGGCTTTGTCTG GCGGGTTGCATATTTGGTG
    12 F07 CAGCATCTTCTACATCTGAGTC CAATCAAAAAGACTACATCACT
    13 H07 CACCATCGGTTTCTCCATC CTCAATCAAAAAGACTACATCACT
    14 D08 TGCTTTCAGAAGGATAAGGGGT AAATACAATACTGGGTTTCGCC
    下载: 导出CSV

    表  3  不同SSR位点的遗传多样性统计

    Table  3.   Genetic diversity of different SSR microsatellite loci

    位点NaNeHoHePICIFisFNm
    H0353.540.000.590.831.181.001.000.68
    D0474.620.010.620.891.380.980.980.60
    H0495.510.040.650.931.550.940.950.70
    B0585.930.000.660.921.570.990.990.78
    F0563.570.000.510.821.150.990.990.43
    H0596.460.000.660.951.591.001.000.63
    B0695.870.010.650.921.540.990.990.67
    D0674.720.000.620.901.391.001.000.55
    F0675.090.020.640.931.470.980.980.62
    H06118.050.570.730.931.900.210.190.92
    D07106.350.010.660.931.620.980.980.64
    F07106.990.020.680.931.700.970.980.69
    H07117.450.140.700.931.790.790.760.76
    D0863.450.000.600.821.250.990.990.50
    平均8.175.540.060.640.901.510.920.910.66
      说明:Na. 等位基因数;Ne. 有效等位基因数;Ho. 观测杂合度;He. 期望杂合度;PIC. 多态信息含量;I. Shannon’s多样性指数;Fis. 居群近交系数;F. 固定指数;Nm. 基因流。
    下载: 导出CSV

    表  4  各亚群基因型情况

    Table  4.   Number of genotypes in each cluster

    亚群各亚群种质数量各亚群基因型数
    Q<0.6Q>0.8
    12012
    24912
    29517
    49045
      说明:Q为样本隶属不同亚群的比例。
    下载: 导出CSV

    表  5  马尾松各亚群遗传多样性水平

    Table  5.   Genetic diversity parameters for all populations of resistant P. massoniana

    亚群样本量/份NaNeIHoHeF
    126.074.871.640.060.780.93
    249.576.481.880.080.780.91
    297.504.881.590.030.700.96
    4913.437.142.160.060.850.93
      说明:Na. 等位基因数;Ne. 有效等位基因数;I. Shannon’s多样性指数;Ho. 观测杂合度;He. 期望杂合度;F. 固定指数。
    下载: 导出CSV

    表  6  基于微卫星数据的遗传变异分析

    Table  6.   Analysis of molecular variance from microsatellite data using GenALEx

    变异
    来源
    自由度方差
    总和
    平均
    方差
    方差
    分量
    变异百
    分比/%
    P
    亚群间3184.86561.6220.95714<0.01
    亚群内1101 216.32811.0585.32380<0.01
    FST0.143
      说明:FST为遗传分化系数。
    下载: 导出CSV

    表  7  不同策略构建的核心种质遗传多样性指标比较

    Table  7.   Comparisons of genetic diversity index of core collections constructed by different tactics

    构建方法种质数/份NaNeIHoHeFNmPIC
    原有种质  1141155.541.510.060.640.910.660.90
    M策略   721155.301.430.060.630.910.610.92
    随机取样策略791155.531.480.060.640.910.650.92
      说明:Na. 等位基因数;Ne. 有效等位基因数;Ho. 观测杂合度;He. 期望杂合度;PIC. 多态性信息含量;I. Shannon’s多样性指数;Fis. 居群近交系数;F. 固定指数;Nm. 基因流。* P<0.05, **P<0.01。
    下载: 导出CSV

    表  8  72份抗性马尾松核心种质与114份原有种质遗传多样性对比

    Table  8.   Comparison of genetic diversity between 72 core germplasm and 114 original germplasm of resistant P. massoniana

    种质种质数/份NaNeIHoHeFNmPIC
    原有种质 114 25.29 5.54 1.51 0.06 0.64 0.91 0.66 0.90
    核心种质 72 25.29 5.30 1.43 0.06 0.63 0.91 0.61 0.92
    保留比例/% 63.16 100.00 95.67 94.96 100.00 98.12 99.57 92.54 100.00
    t 0.000 −0.652 −1.428 0.059 −0.826 −0.052 0.014 −0.152
    P 1.000 0.525 0.177 0.953 0.424 0.960 0.989 0.882
      说明:Na. 等位基因数;Ne. 有效等位基因数;Ho. 观测杂合度;He. 期望杂合度;PIC. 多态信息含量;I. Shannon’s多样性指数;Fis. 居群近交系数;F. 固定指数;Nm. 基因流。
    下载: 导出CSV
  • [1] 叶建仁, 吴小芹. 松材线虫病研究进展[J]. 中国森林病虫, 2022, 41(3): 1 − 10.

    YE Jianren, WU Xiaoqin. Research progress of pine wilt disease [J]. Forest Pest and Disease, 2022, 41(3): 1 − 10.
    [2] 刘彬, 刘青华, 周志春, 等. 马尾松β-蒎烯合酶基因克隆以及对松材线虫侵染的响应[J]. 林业科学研究, 2020, 33(6): 1 − 12.

    LIU Bin, LIU Qinghua, ZHOU Zhichun, et al. Cloning of β-pinene synthase gene of Pinus massoniana and response to infection of pine wood nematode [J]. Forest Research, 2020, 33(6): 1 − 12.
    [3] 刘彬, 刘青华, 周志春, 等. 基于高通量转录组测序筛选马尾松抗松材线虫病相关基因[J]. 林业科学研究, 2019, 32(5): 1 − 10.

    LIU Bing, LIU Qinghua, ZHOU Zhichun, et al. Identification of candidate constitutive expressed resistant genes of pine wilt disease in Pinus massoniana based on high-throughput transcriptome sequencing [J]. Forest Research, 2019, 32(5): 1 − 10.
    [4] 徐六一, 章健, 高景斌, 等. 安徽省松材线虫病抗性育种研究进展[J]. 安徽林业科技, 2013, 39(2): 8 − 10, 14.

    XU Liuyi, ZHANG Jian, GAO Jingbin, et al. Research progress on breeding to pinewood nematodiasis in Anhui Province [J]. Anhui Forestry Science and Technology, 2013, 39(2): 8 − 10, 14.
    [5] 高景斌, 席启俊, 孙主义, 等. 松材线虫病抗性马尾松苗木的选育[J]. 林业科技开发, 2009, 23(1): 91 − 95.

    GAO Jingbin, XI Qijun, SUN Zhuyi, et al. Screening and breeding of Pinus massoniana seedlings for resistant to pine wood nematode [J]. Journal of Forestry Engineering, 2009, 23(1): 91 − 95.
    [6] 郝焰平, 徐六一, 姜春武, 等. 安徽省适生马尾松优良种源子代抗松材线虫病评价及生长性状研究[J]. 江西农业学报, 2019, 31(3): 41 − 45.

    HAO Yanping, XU Liuyi, JIANG Chunwu, et al. Nematode-resistance evaluation and growth trait study of fine offspring of suitable Pinus massoniana provenance in Anhui Province [J]. Acta Agriculturae Jiangxi, 2019, 31(3): 41 − 45.
    [7] FRANKEL O H, BROWN A H D. Plant genetic resources today: a critical appraisal [M]// HOLDON J H W. Crop Genetic Resources: Conservation and Evaluation. London: George Allen and Unwin Press, 1984: 249 − 257.
    [8] 徐晓美, 李颖, 孙启迪, 等. 辣椒种质材料疫病抗性鉴评及遗传多样性分析[J]. 广东农业科学, 2022, 49(10): 19 − 28.

    XU Xiaomei, LI Ying, SUN Qidi, et al. Resistance evaluation and genetic diversity analysis of Phytophthora disease in pepper germplasm materials [J]. Guangdong Agricultural Sciences, 2022, 49(10): 19 − 28.
    [9] 周延清. 遗传标记的发展[J]. 生物学通报, 2000, 35(5): 17 − 18.

    ZHOU Yanqing. Development of genetic markers [J]. Bulletin of Biology, 2000, 35(5): 17 − 18.
    [10] 杨海平, 李继生, 于国强, 等. 分子标记技术在林木育种中的应用[J]. 山东林业科技, 2017, 47(3): 111 − 114.

    YANG Haiping, LI Jisheng, YU Guoqiang, et al. Application of molecular marker technology in forest tree breeding [J]. Shandong Forestry Science and Technology, 2017, 47(3): 111 − 114.
    [11] 杨汉波, 张蕊, 王帮顺, 等. 基于SSR标记的木荷核心种质构建[J]. 林业科学, 2017, 53(6): 37 − 46.

    YANG Hanbo, ZHANG Rui, WANG Bangshun, et al. Construction of core collection of Schima superba based on SSR molecular markers [J]. Scientia Silvae Sinicae, 2017, 53(6): 37 − 46.
    [12] 陈存, 丁昌俊, 张静, 等. 美洲黑杨群体结构分析及核心种质库构建[J]. 林业科学, 2020, 56(9): 67 − 76.

    CHEN Cun, DING Changjun, ZHANG Jing, et al. Population structure analysis and core collection construction of Populus deltoides [J]. Scientia Silvae Sinicae, 2020, 56(9): 67 − 76.
    [13] 陈明堃, 陈璐, 孙维红, 等. 建兰种质资源遗传多样性分析及核心种质构建[J]. 园艺学报, 2022, 49(1): 175 − 186.

    CHEN Mingkun, CHEN Lu, SUN Weihong, et al. Genetic diversity analysis and core collection of Cymbidium ensifolium germplasm resources [J]. Acta Horticulturae Sinica, 2022, 49(1): 175 − 186.
    [14] RAJPUT S G, SANTRA D K. Evaluation of genetic diversity of Proso millet germplasm available in the United States using simple sequence repeat markers [J]. Crop Science, 2016, 56(5): 2401 − 2409.
    [15] LÜ Jiabin, LI Changrong, ZHOU Changpin, et al. Genetic diversity analysis of a breeding population of Eucalyptus cloeziana F. Muell. (Myrtaceae) and extraction of a core germplasm collection using microsatellite markers [J/OL]. Industrial Crops and Products, 2020, 145: 112157[2023-05-26]. doi: 10.1016/j.indcrop.2020.112157.
    [16] 唐玉娟, 罗世杏, 黄国弟, 等. 基于SSR荧光标记的杧果种质资源遗传多样性分析及分子身份证构建[J/OL]. 热带作物学报, 2023-01-17[2023-05-26]. http://kns.cnki.net/kcms/detail/46.1019.s.20221122.1956.004.html.

    TANG Yujuan, LUO Shixing, HUANG Guodi, et al. Genetic diversity analysis and molecular ID construction of mango germplasm based on SSR fluorescence marker [J/OL]. Chinese Journal of Tropical Crops: 2023-01-17[2023-05-26]. http://kns.cnki.net/kcms/detail/46.1019.s.20221122.1956.004.html.
    [17] HU Guangming, JIANG Quan, WANG Zhi, et al. Genetic diversity analysis and core collection construction of the Actinidia chinensis complex (Kiwi fruit) based on SSR Markers [J/OL]. Agronomy, 2022, 12(12): 3078[2023-05-26]. doi: 10.3390/agronomy12123078.
    [18] 董虹妤, 刘青华, 周志春, 等. 马尾松子代生长杂种优势与亲本配合力、遗传距离的相关性[J]. 林业科学, 2017, 53(2): 65 − 75.

    DONG Hongyu, LIU Qinghua, ZHOU Zhichun, et al. Correlation between heterosis in the growth of progeny and combining ability and genetic distance of the parents for Pinus massoniana [J] Scientia Silvae Sinicae, 2017, 53(2): 65 − 75.
    [19] YANG Hanbo, LIU Qinghua, ZHANG Rui, et al. Genetic diversity of second generation-parental germplasm of masson pine revealed by SSR markers and establishment of a core germplasm collection [J]. Scandinavian Journal of Forest Research, 2021, 36(7/8): 524 − 531.
    [20] ROD P, PETER E S. GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research an update [J]. Bioinformatics, 2012, 28(19): 2537 − 2539.
    [21] KIALINOWSK S, TAPER M L, MARSHALL T C. Revising how the computer program cervus accommodates genotyping error increases success in paternity assignment: CERVUS LIKELIHOOD MODEL [J]. Molecular Ecology, 2007, 16(5): 1099 − 1106.
    [22] PRITCHARD J K, STEPHENS M, DONNELLY P. Inference of population structure using multilocus genotype data [J]. Genetics, 2000, 155(2): 945 − 959.
    [23] KIM K W, CHUNG H K, CHO G T, et al. PowerCore: a program applying the advanced M strategy with a heuristic search for establishing core sets [J]. Bioinformatics, 2007, 23(16): 2155 − 2162.
    [24] CIPRIANI G, SPADOTTO A, JURMAN L, et al. The SSR-based molecular profile of 1005 grapevine (Vitis vinifera L. ) accessions uncovers new synonymy and parentages, and reveals a large admixture amongst varieties of different geographic origin [J]. TAG Theoretical and Applied Genetics, 2010, 121(8): 1569 − 1585.
    [25] LIU K, MUSE S V. PowerMarker: integrated analysis environment for genetic marker data [J]. Bioinformatics, 2005, 21(9): 2128 − 2129.
    [26] 杨梅, 张敏, 师守国, 等. 武当木兰种群遗传结构的ISSR分析[J]. 林业科学, 2014, 50(1): 76 − 81.

    YANG Mei, ZHANG Min, SHI Shouguo, et al. Analysis of genetic structure of Magonlia sprengeri populations based on ISSR marker [J]. Scientia Silvae Sinicae, 2014, 50(1): 76 − 81.
    [27] 沈敬理, 白天道, 陈亚斌, 等. 马尾松无性系种子园亲本指纹图谱构建[J]. 分子植物育种, 2015, 13(3): 646 − 652.

    SHEN Jingli, BAI Tiandao, CHEN Yabin, et al. The fingerprints construction of clonal in Pinus massoniana seed orchard [J]. Molecular Plant Breeding, 2015, 13(3): 646 − 652.
    [28] 高景斌, 徐六一, 叶建仁. 马尾松松材线虫病抗性无性系的筛选和遗传多样性分析[J]. 南京林业大学学报(自然科学版), 2021, 45(5): 109 − 118.

    GAO Jingbin, XU Liuyi, YE Jianren. Growth and genetic diversity analysis of clones screened by phenotypical resistance to pine wilt disease in Pinus massoniana [J]. Journal of Nanjing Forestry University (Natural Sciences Edition), 2021, 45(5): 109 − 118.
    [29] 谢立峰, 李宁, 李烨, 等. 茄子种质遗传多样性及群体结构的SRAP分析[J]. 植物学报, 2019, 54(1): 58 − 63.

    XIE Lifeng, LI Ning, LI Ye, et al. Genetic diversity and population structure of eggplant (Solanum melongena) germplasm resources based on SRAP method [J]. Chinese Bulletin of Botany, 2019, 54(1): 58 − 63.
    [30] 魏利斌, 苗红梅, 李春, 等. 芝麻SNP和InDel标记遗传多样性、群体结构及连锁不平衡分析[J]. 分子植物育种, 2017, 15(8): 3070 − 3079.

    WEI Libin, MIAO Hongmei, LI Chun, et al. Genetic diversity, population strcture and linkage disequilibrium analysis of Sesame using SNP and InDel markers [J]. Molecular Plant Breeding, 2017, 15(8): 3070 − 3079.
    [31] 方乐成, 夏慧敏, 麻文俊, 等. 基于SSR标记的楸树遗传多样性及核心种质构建[J]. 东北林业大学学报, 2017, 45(8): 1 − 5.

    FANG Lecheng, XIA Huimin, MA Wenjun, et al. Genetic diversity analysis and primary core collection of Catalpa bungei germplasm with SSR markers [J]. Journal of Northeast Forestry University, 2017, 45(8): 1 − 5.
    [32] WANG J C, HU J, XU H M, et al. A strategy on constructing core collections by least distance stepwise sampling [J]. Theoretical and Applied Genetics, 2007, 115(1): 1 − 8.
    [33] ESCRIBANO P, VIRUEL M A, HORMAZA J I. Comparison of different methods to construct a core germplasm collection in woody perennial species with simple sequence repeat markers. a case study in cherimoya (Annona cherimola, Annonaceae), an underutilised subtropical fruit tree species [J]. Annals of Applied Biology, 2008, 153(1): 25 − 32.
    [34] 李魁鹏, 陈仕昌, 程琳, 等. 基于SSR标记构建广西杉木核心种质[J]. 广西科学, 2021, 28(5): 511 − 519.

    LI Kuipeng, CHEN Shichang, CHENG Lin, et al. Construction of core germplasm of Cunninghamia lanceolata in Guangxi based on SSR marker [J]. Guangxi Sciences, 2021, 28(5): 511 − 519.
    [35] 黄小凤, 韦阳连, 袁叶, 等. 基于SNP分子标记的221份荔枝品种(品系)的遗传多样性分析及核心种质库构建[J]. 植物资源与环境学报, 2022, 31(4): 74 − 84.

    HUANG Xiaofeng, WEI Yanglian, YUAN Ye, et al. Genetic diversity analysis and core collection construction of 221 cultivars (strains) of Litchi chinensis based on SNP molecular markers [J]. Journal of Plant Resources and Environment, 2022, 31(4): 74 − 84.
    [36] 张馨方, 张树航, 李颖, 等. 基于SSR标记构建燕山板栗核心种质[J]. 华北农学报, 2021, 36(增刊1): 31 − 38.

    ZHANG Xinfang, ZHANG Shuhang, LI Ying, et al. Construction of core collection of Yanshan chestnut germplasm based on SSR marker [J]. Acta Agriculturae Boreali-Sinica, 2021, 36(suppl 1): 31 − 38.
    [37] 李自超, 张洪亮, 曾亚文, 等. 云南地方稻种资源核心种质取样方案研究[J]. 中国农业科学, 2000, 33(5): 1 − 7.

    LI Zichao, ZHANG Hongliang, ZENG Yawen, et al. Study on sampling schemes of core collection of local varieties of rice in Yunnan, China [J]. Scientia Agricultura Sinica, 2000, 33(5): 1 − 7.
  • [1] 黄锦春, 万思琦, 陈扬, 李丽红, 张自力, 朱建军, 吴梅, 邢丙聪, 邵清松, 陆晨飞.  利用ISSR与SRAP分子标记分析金线莲种质资源遗传多样性 . 浙江农林大学学报, 2023, 40(1): 22-29. doi: 10.11833/j.issn.2095-0756.20220473
    [2] 王洁, 朱锡彭, 王腾斐, 朱建军, 李文俊, 邢丙聪, 郑颖.  白及种质资源遗传多样性分析 . 浙江农林大学学报, 2023, 40(2): 321-329. doi: 10.11833/j.issn.2095-0756.20210690
    [3] 郑刘辉, 詹利云, 侯宇, 喻卫武, 曾燕如, 戴文圣.  雄性榧树遗传多样性的SSR荧光标记分析 . 浙江农林大学学报, 2022, 39(2): 329-337. doi: 10.11833/j.issn.2095-0756.20210279
    [4] 韩泽民, 李源, 王熊, 菅永峰, 周靖靖, 佃袁勇, 黄光体.  不同演替程度下马尾松人工林生物多样性对生物量的影响 . 浙江农林大学学报, 2021, 38(2): 246-252. doi: 10.11833/j.issn.2095-0756.20200334
    [5] 尹焕焕, 刘青华, 周志春, 万雪琴, 余启新, 丰忠平.  马尾松无性系木材基本密度和纤维形态的变异及选择 . 浙江农林大学学报, 2020, 37(6): 1186-1192. doi: 10.11833/j.issn.2095-0756.20190720
    [6] 王罗云, 何彩云, 罗红梅, 张建国, 段爱国, 曾艳飞.  沙棘7个亚种与26个重要品种的遗传多样性 . 浙江农林大学学报, 2019, 36(4): 670-677. doi: 10.11833/j.issn.2095-0756.2019.04.006
    [7] 柴晓娟, 郭玮龙, 陈慧, 王俊荣, 谢德志, 金水虎.  蔺草种质资源遗传多样性的ISSR分析 . 浙江农林大学学报, 2017, 34(3): 552-558. doi: 10.11833/j.issn.2095-0756.2017.03.022
    [8] 马中青, 张齐生.  温度对马尾松热解产物产率和特性的影响 . 浙江农林大学学报, 2016, 33(1): 109-115. doi: 10.11833/j.issn.2095-0756.2016.01.015
    [9] 吕常笑, 邓华锋, 王少杰, 陈振雄, 王雪军.  马尾松不同区域相容性立木材积和地上生物量模型 . 浙江农林大学学报, 2016, 33(5): 790-797. doi: 10.11833/j.issn.2095-0756.2016.05.010
    [10] 刘跃钧, 张媛, 蒋燕锋, 刘京晶.  黄精种质资源遗传多样性研究 . 浙江农林大学学报, 2016, 33(6): 1085-1091. doi: 10.11833/j.issn.2095-0756.2016.06.023
    [11] 刘青华, 周志春, 张开明, 兰永兆, 吴吉富, 聂国勤.  施用磷肥对马尾松种源生长和木材基本密度的影响 . 浙江农林大学学报, 2012, 29(2): 185-191. doi: 10.11833/j.issn.2095-0756.2012.02.006
    [12] 张振, 张含国, 张磊, 朱航勇, 李雪峰.  兴安落叶松基本群体与育种群体RAPD多样性分析 . 浙江农林大学学报, 2012, 29(1): 130-136. doi: 10.11833/j.issn.2095-0756.2012.01.022
    [13] 高燕会, 樊民亮, 骆文坚, 黄华宏, 童再康.  濒危树种金钱松RAPD体系的建立和遗传多样性分析 . 浙江农林大学学报, 2011, 28(5): 815-822. doi: 10.11833/j.issn.2095-0756.2011.05.023
    [14] 李永进, 丁贵杰.  不同家系马尾松插穗内源生根抑制物的分离、纯化及鉴定 . 浙江农林大学学报, 2010, 27(4): 507-512. doi: 10.11833/j.issn.2095-0756.2010.04.005
    [15] 何佩云, 丁贵杰.  猴樟、鹅掌楸对马尾松苗木生理活性的他感效应 . 浙江农林大学学报, 2008, 25(5): 604-608.
    [16] 唐娟娟, 范义荣, 朱睦元.  黄山松群体遗传多样性分析 . 浙江农林大学学报, 2003, 20(1): 23-26.
    [17] 郑仁华, 陈国金, 傅忠华, 俞白楠, 杨宗武, 傅玉狮, 潘琼蓉.  马尾松优树子代遗传评价及选择 . 浙江农林大学学报, 2001, 18(2): 144-149.
    [18] 季孔庶.  马尾松种源针叶内痕量元素数量变异与生长相关性* . 浙江农林大学学报, 1997, 14(3): 231-236.
    [19] 徐有明, 林汉, 万伏红.  马尾松纸浆材材性变异和采伐林龄的确定* . 浙江农林大学学报, 1997, 14(1): 8-15.
    [20] 林思祖, 黄青峥, 吴旺民.  杉木马尾松几个模型的预报有效性的计算机模拟试验 . 浙江农林大学学报, 1993, 10(2): 184-188.
  • 加载中
  • 链接本文:

    https://zlxb.zafu.edu.cn/article/doi/10.11833/j.issn.2095-0756.20230333

    https://zlxb.zafu.edu.cn/article/zjnldxxb/2024/1/67

图(6) / 表(8)
计量
  • 文章访问数:  417
  • HTML全文浏览量:  70
  • PDF下载量:  32
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-05-29
  • 修回日期:  2023-10-09
  • 录用日期:  2023-10-24
  • 网络出版日期:  2024-01-19
  • 刊出日期:  2024-02-20

抗松材线虫病马尾松种质资源遗传多样性分析及核心种质构建

doi: 10.11833/j.issn.2095-0756.20230333
    基金项目:  浙江省“十四五”育种专项林木协作组课题(2020C02007);江西省林业局林业科技创新专项(创新专项﹝2021﹞13号)
    作者简介:

    邓莉丽(ORCID: 0009-0002-3535-9972),从事林木遗传育种研究。E-mail: Lili_Deng1706@163.com

    通信作者: 刘青华(ORCID: 0000-0002-7960-6739),研究员,博士,从事抗性马尾松遗传改良研究。E-mail: liuqinghua642@163.com
  • 中图分类号: S722.3

摘要:   目的  对来自安徽省林业科学研究院和异地保存在浙江省临海市林业技术推广和场圃旅游服务总站的114份抗松材线虫Bursaphelenchus xylophilus病马尾松Pinus massoniana种质进行遗传多样性与群体结构分析,构建抗性马尾松核心种质库。  方法  对114份抗性马尾松种质进行检测,使用分子生物学软件计算遗传多样性参数,并进行主坐标(PCoA)和群体结构分析。利用M策略和随机取样策略分别构建核心种质,分析不同核心种质的的遗传多样性参数,确定最适合的构建方法。  结果  114份抗性马尾松种质共检测到115个等位基因,平均有效等位基因数(Ne)为5.54,平均Shannon’s多样性指数(I)为1.51,平均多态信息含量(PIC)为0.90,结果表明其具有较高的遗传多样性水平。群体结构分析将114份抗性马尾松种质分为4个亚群,主坐标分析结果与上述基本一致。根据遗传多样性参数和抽样数量综合考虑,M策略构建的核心种质能以最小的种质数保留原有种质最大的遗传多样性,为最佳的取样策略。利用该策略得到了72份核心种质,其保留了原有种质100%的等位基因数,NeI、期望杂合度(He)、PIC等遗传参数的保留率分别为95.67%、94.96%、98.12%和100.00%,将构建的核心种质与原有种质进行t检验、PCoA分析和UPGMA聚类分析,结果表明两者间的遗传多样性无显著差异。  结论  114份抗性马尾松种质遗传多样性水平较高,构建的72份抗性马尾松核心种质,去除了遗传冗余,有利于抗性马尾松种质资源的有效保护和科学利用,可为优异基因发掘和新品种选育提供参考。图7表8参37

English Abstract

邓莉丽, 刘青华, 周志春, 高凯, 骆定会. 抗松材线虫病马尾松种质资源遗传多样性分析及核心种质构建[J]. 浙江农林大学学报, 2024, 41(1): 67-78. doi: 10.11833/j.issn.2095-0756.20230333
引用本文: 邓莉丽, 刘青华, 周志春, 高凯, 骆定会. 抗松材线虫病马尾松种质资源遗传多样性分析及核心种质构建[J]. 浙江农林大学学报, 2024, 41(1): 67-78. doi: 10.11833/j.issn.2095-0756.20230333
DENG Lili, LIU Qinghua, ZHOU Zhichun, GAO Kai, LUO Dinghui. Genetic diversity analysis and core collection of pinewood nematodiasis-resistant Pinus massoniana germplasm resources[J]. Journal of Zhejiang A&F University, 2024, 41(1): 67-78. doi: 10.11833/j.issn.2095-0756.20230333
Citation: DENG Lili, LIU Qinghua, ZHOU Zhichun, GAO Kai, LUO Dinghui. Genetic diversity analysis and core collection of pinewood nematodiasis-resistant Pinus massoniana germplasm resources[J]. Journal of Zhejiang A&F University, 2024, 41(1): 67-78. doi: 10.11833/j.issn.2095-0756.20230333
  • 松材线虫Bursaphelenchus xylophilus病是20 世纪80 年代初传入中国的一种毁灭性森林病害,其传播速度快,防治难度大,危害中国松林面积超180 万 hm2,其中,马尾松Pinus massoniana最易受感染,松材线虫病对中国马尾松林业经济发展造成了极大破坏[1]。马尾松具有速生、耐旱等优良特点,松林面积达804 万 hm2,是中国南方荒山造林的先锋树种[2]。近40 a的林业生产实践和研究发现,松树种间和种内不同个体对松材线虫病的抗性存在较大差异,且马尾松群体中存在抗病基因型,因此选育抗性品种是当前治理松材线虫病害最经济有效的途径之一[3]。自2001年开始,中国开展了马尾松松材线虫病抗性育种研究,现已收集保存一批抗松材线虫病马尾松种质资源 [46]。种质资源是遗传改良的物质基础,种质资源收集保存的越多,其群体遗传多样性越高,但收集保存的成本也越高。FRANKEL等[7]提出核心种质的概念,即用最小的遗传资源数量和最小的遗传冗余最大限度地代表整个遗传资源的多样性。通过构建核心种质,能最大限度地去除重复和遗传关系较近的种质材料,加强现有抗性马尾松种质资源的有效保护利用,降低管理成本[8]。分子标记法是林木遗传多样性分析、种质资源评价和核心种质库构建常用的方法[913]。其中简单序列重复(SSR)分子标记具有多等位基因、共显性和高度多态性等优点,能够直观反映不同种质间遗传信息差异,在遗传育种中得到广泛应用[14]。LÜ等[15]利用 SSR 分子标记数据研究了桉树Eucalyptus cloeziana的遗传多样性并构建了核心种质;唐玉娟等[16]利用12对 SSR 引物对杧果Mangifera indica种质资源的遗传多样性进行分析并构建了分子身份证;HU等[17]利用 SSR 分子标记技术对猕猴桃Actinidia chinensis的遗传多样性进行研究,构建了猕猴桃核心种质。

    本研究利用 SSR 分子标记技术对前期筛选获得的103份抗松材线虫病马尾松种质和安徽省林业科学研究院审定的高抗良种进行遗传多样性及亲缘关系分析,并基于M策略和随机取样策略构建核心种质,通过分析不同构建策略的核心种质确定抗性马尾松核心种质的最适取样比例,为抗性马尾松种质资源的合理开发和利用以及挖掘优异基因提供理论支撑。

    • 2018年,在国家马尾松种质资源库中依据松脂和抗性的关系,从1 207株树中选择120份候选抗性马尾松无性系,同年在松材线虫病重灾区严重感染林分中收集健康树242份,并将所有无性系嫁接保存于浙江省临海市林业技术推广和场圃旅游服务总站(28°88′N,121°01′E)。2018—2020年连续3 a高强度对收集区内的候选抗性马尾松无性系进行人工接种松材线虫测定,筛选存活健康株系,共获得103份抗性马尾松无性系,另取安徽省林业科学研究院审定的抗松材线虫病良种(为目前中国所有的抗松材线虫病良种)作为对照。合计114份,其中浙江37份、江西26份、安徽42份、福建6份、贵州1份、广西2份。采集抗性样株幼嫩针叶,在−80 ℃冰箱中保存,各试验材料编号及来源见表1

      表 1  材料编号及来源

      Table 1.  Test materials’ number and source

      编号来源编号来源编号来源编号来源
      1广西贵港35安徽黄山DAC12浙江临海HGS5浙江临海
      2广西桂平36安徽黄山DAC13浙江临海HGS6浙江临海
      3贵州黎平37安徽全椒DAC15浙江临海HGS8浙江临海
      4福建南平38安徽泾县DAC17浙江临海HGS9浙江临海
      6福建南平39安徽休宁DAC18浙江临海HGY1浙江临海
      7安徽休宁40安徽休宁DAC19浙江临海LS001浙江临海
      9安徽休宁41安徽休宁DAC20浙江临海LS003浙江临海
      10安徽祁县42安徽休宁DAC21浙江临海LS005浙江临海
      11福建南平43安徽休宁DAC22浙江临海LS007浙江临海
      12浙江淳安44安徽休宁DAC23浙江临海LS008浙江临海
      13浙江淳安45安徽休宁DAC24浙江临海LS009浙江临海
      14福建南平46安徽休宁DAC26浙江临海LS010浙江临海
      15浙江松阳47安徽休宁DAC34浙江临海XJ003江西峡江
      17福建南平48安徽休宁DLS27浙江临海XJ004江西峡江
      19江西婺源49安徽休宁DLS28浙江临海XJ005江西峡江
      20福建南平71江西分宜DLS29浙江临海XJ014江西峡江
      21安徽广德114江西分宜DLS30浙江临海AH1安徽广德
      22安徽广德161江西分宜DLS31浙江临海AH2安徽滁州
      23安徽广德231江西分宜DLS32浙江临海AH3安徽黄山
      24安徽广德242江西分宜DLS33浙江临海AH4安徽休宁
      25安徽广德305江西分宜DLS35浙江临海AH5安徽广德
      26安徽广德355江西分宜DLS36浙江临海AH6安徽广德
      27安徽广德398江西分宜DLS37浙江临海AH7安徽广德
      28安徽广德413江西分宜DLS38浙江临海AH8安徽泾县
      29安徽广德461江西分宜DLS40浙江临海AH9安徽广德
      30安徽和县473江西分宜HGS11浙江临海AH10安徽休宁
      31安徽和县485江西分宜HGS2浙江临海AH11安徽广德
      32安徽和县486江西分宜HGS3浙江临海
      34安徽和县490江西分宜HGS4浙江临海
        说明:AH1~AH11为 安徽省林业科学研究院审定的马尾松抗松材线虫病良种皖马抗1~11号。
    • 马尾松针叶基因组DNA采用北京艾德莱生物科技有限公司生产的改良CTAB植物基因组DNA快速提取试剂盒提取。利用超微量分光光度计,根据在260 nm的吸光度D(260)和D(260)/D(280)确定DNA样品的纯度和质量浓度,经检验合格后的DNA用TE溶液或去离子水(灭菌后)稀释至20~50 mg·L−1,置于−20 ℃的冰箱保存备用。

      参考董虹妤等[18]和YANG等[19]使用的马尾松 SSR 引物,用质量浓度为1.2%的琼脂糖凝胶电泳检测,筛选条带清晰且多态性较高的SSR引物(图1),最终筛选出具有高多态性引物14对(表2),引物序列由浙江尚亚生物技术有限公司合成。

      图  1  引物多态性检测图

      Figure 1.  Primer polymorphism detection diagram

      表 2  14对SSR引物信息

      Table 2.  14 pairs of SSR primer information

      编号引物名称正向引物(5´→3´)反向引物(5´→3´)
      1 H03 CTCCAAAGGCGAGACTGC ACGAAAGCCAAGCTGAAC
      2 D04 AGGATGGTATGGTCGTGG CCCTTCTTCGCTCTGTGA
      3 H04 AAGAAGCGATCTGAGATGACTAA CTGTCATTGATTGTTTCCTTTTG
      4 B05 CCGTGCCTTCAGCATCTTCT CAGTGGATCTGTCACCTCCTCAT
      5 F05 AGAAGAAGAGCAGCAGTTTC GGTTTTCCATTGTTCTCACT
      6 H05 GTGCCTTCAGCATCTTCTAC ATCTGTCACCTCCTCATCTT
      7 B06 TATTAGCACCCCTCCAAAG TGGGGTAGAAGAATCGTAAGT
      8 D06 GATTCGGCTTCGTGACCTT CCCCCATAACCCCTGTC
      9 F06 AAGGACTTACAGAGGTTGGGTT GCTGCGACGAGCGTTTCT
      10 H06 TTCCTACCGCTGGGTTCTTG GACCTGACCTCGGGCATTAC
      11 D07 TCGCCTGGGCTTTGTCTG GCGGGTTGCATATTTGGTG
      12 F07 CAGCATCTTCTACATCTGAGTC CAATCAAAAAGACTACATCACT
      13 H07 CACCATCGGTTTCTCCATC CTCAATCAAAAAGACTACATCACT
      14 D08 TGCTTTCAGAAGGATAAGGGGT AAATACAATACTGGGTTTCGCC
    • 利用选取的14对 SSR 引物对114份马尾松DNA样本进行PCR扩增,设定本研究的反应体系为25.0 μL:2×Taq Plus Master Mix (Nazyme Code:P211-03)12.5 μL,引物F和R (10 μmol·L−1)各1.0 μL,DNA模板2.0 μL,ddH2O 8.5 μL。PCR扩增反应在Takara PCR Thermal cycler上进行,其程序为:95 ℃预变性3 min;95 ℃变性15 s,60 ℃退火15 s,72 ℃延伸30 s,35个循环;72 ℃延伸5 min,4 ℃保存。所得PCR产物采用Qsep100全自动毛细管电泳核酸分析仪进行检测和基因分型。

    • 利用GenAlex 6.5[20]计算各位点的遗传多样性参数,包括等位基因数(Na)、有效等位基因数(Ne)、观测杂合度(Ho)、期望杂合度(He)、Shannon’s多样性指数(I)、近交系数(Fis)、固定指数(F)和基因流(Nm),并对收集的抗松材线虫病马尾松种质进行主坐标分析(PCoA)。采用Cervus 2.0[21]获得微卫星位点的多态信息含量(PIC)。利用Structure 2.3.4[22]对抗性马尾松种质资源进行分类,推测最佳聚类组群数目,确定群体结构,并对各亚群的遗传多样性参数进行分子方差分析(AMOVA)。

    • 对抗松材线虫病马尾松种质进行群体结构分析后,利用Power Core[23]和Core Finder[24]对各亚群基于不同算法构建核心种质。其中,Power Core是利用随机取样策略和启发式算法,寻找从初始阶段到最终阶段的最优路径进行核心集合选择;Core Finder基于NP-完全覆盖问题,采用拉斯维加斯式(LasVegas-style)随机算法的M策略进行样本选择。在构建核心种质过程中,Power Core和Core Finder会自动生成合理的取样比例。最后利用SPSS 26.0和Excel中t检验对原有种质与核心种质各遗传多样性指标进行差异显著性分析,同时运用PCoA和Powermarker[25]对所构建的核心种质构建UPGMA进化树,确认构建的核心种质的代表性。

    • 表3可知:14对 SSR 引物在114份抗松材线虫病马尾松样本中共检测出115个等位位点,Na为5~11,平均为8.17。所有材料的Ne为3.45~8.05,平均为5.54。HoHe的分别为0~0.57和0.51~0.73,He均值(0.64)高于Ho均值(0.06),Fis平均为0.92,综合表明该抗性马尾松种质中存在杂合子缺失现象。PIC为0.82~0.95,平均为0.9,I为1.15~1.90,平均为1.51,表明所选引物多态性高。

      表 3  不同SSR位点的遗传多样性统计

      Table 3.  Genetic diversity of different SSR microsatellite loci

      位点NaNeHoHePICIFisFNm
      H0353.540.000.590.831.181.001.000.68
      D0474.620.010.620.891.380.980.980.60
      H0495.510.040.650.931.550.940.950.70
      B0585.930.000.660.921.570.990.990.78
      F0563.570.000.510.821.150.990.990.43
      H0596.460.000.660.951.591.001.000.63
      B0695.870.010.650.921.540.990.990.67
      D0674.720.000.620.901.391.001.000.55
      F0675.090.020.640.931.470.980.980.62
      H06118.050.570.730.931.900.210.190.92
      D07106.350.010.660.931.620.980.980.64
      F07106.990.020.680.931.700.970.980.69
      H07117.450.140.700.931.790.790.760.76
      D0863.450.000.600.821.250.990.990.50
      平均8.175.540.060.640.901.510.920.910.66
        说明:Na. 等位基因数;Ne. 有效等位基因数;Ho. 观测杂合度;He. 期望杂合度;PIC. 多态信息含量;I. Shannon’s多样性指数;Fis. 居群近交系数;F. 固定指数;Nm. 基因流。
    • Structure分析得出的结果经Structure Haevester处理后发现(图2A):lnP(D)(马尔科夫链不作数迭代后所得的对数似然值)持续增大,没有拐点,可用∆K [L(K)在连续K之间变化率的增量]来确定最优分组数(K)。由图2B可知:当K=4时,∆K取得最大值,表明参试的114份马尾松样本被划分为4个亚群。亚群Ⅰ(绿色)包含了12份种质资源,其中2份来自广西,1份来自贵州,3份来自福建,2份来自安徽,4份来自浙江;亚群Ⅱ(黄色)包含24份种质资源,其中8份来自安徽,7份来自江西,9份来自浙江,该亚群有12份抗性马尾松种质掺杂有来自红色部分的基因;亚群Ⅲ(红色)包含了29份种质资源,其中10份来自安徽(均为皖马抗良种),8份来自江西,11份来自浙江,共有12份抗性马尾松种质掺杂有来自黄色部分的基因;亚群Ⅳ(蓝色)包含49份种质资源,其中2份来自福建,21份来自安徽(其中1份为皖马抗良种),11份来自江西,15份来自浙江(图3)。由Structure输出的分组结果可知(表4):筛选获得的抗性马尾松种质中有75.4%的个体样本隶属不同亚群的比例(Q)>0.8,仅有12.3%的个体Q<0.6,表明抗性马尾松种质群体结构较强,不同来源的个体被聚在相同的遗传结构中表明不同省间的马尾松存在基因渗入现象。

      图  2  lnP(D)和∆KK变化的折线图

      Figure 2.  Line plots of lnP(D) and ∆K as a function of K

      图  3  114份抗性马尾松样本的Structure聚类图

      Figure 3.  Structure cluster diagram of 114 samples of resistant P. massoniana

      表 4  各亚群基因型情况

      Table 4.  Number of genotypes in each cluster

      亚群各亚群种质数量各亚群基因型数
      Q<0.6Q>0.8
      12012
      24912
      29517
      49045
        说明:Q为样本隶属不同亚群的比例。

      基于遗传距离矩阵对抗性马尾松材料进行主坐标分析(PCoA)。由图4可知:全部马尾松种质资源可大致分为3组,主坐标1和2分别解释了位点信息数据中10.45%和6.74%的变异。基于抗性马尾松114份样本主坐标分析图可知:大部分地理来源一致的抗性马尾松种质分布比较集中,但来自安徽、江西、浙江的部分抗性马尾松种质在3个组中均有分布。结合群体结构分析结果可知:114份抗性马尾松种质资源的分布情况基本符合群体结构分析。

      图  4  抗性马尾松114份样本主坐标分析

      Figure 4.  Principal coordinate analysis of 114 samples

    • 表5可见:抗性马尾松种质各亚群的遗传多样性整体水平较高。亚群Ⅳ的等位基因数在4个亚群中最大(13.43),亚群Ⅰ具有最低的等位基因数(6.07)。抗性马尾松4个亚群所包含的马尾松个体数不同,但亚群间遗传多样性水平差异不明显。亚群Ⅳ的Ne最高,为7.14,而最低的是亚群Ⅰ,为4.87。亚群Ⅳ的HoHe均最高,但与具有最小值的亚群Ⅲ差异不显著,分别高0.03和0.15。

      表 5  马尾松各亚群遗传多样性水平

      Table 5.  Genetic diversity parameters for all populations of resistant P. massoniana

      亚群样本量/份NaNeIHoHeF
      126.074.871.640.060.780.93
      249.576.481.880.080.780.91
      297.504.881.590.030.700.96
      4913.437.142.160.060.850.93
        说明:Na. 等位基因数;Ne. 有效等位基因数;I. Shannon’s多样性指数;Ho. 观测杂合度;He. 期望杂合度;F. 固定指数。

      AMOVA方差分析表明(表6):抗性马尾松遗传变异14%存在于亚群间,而亚群内变异为80%,表明亚群内的遗传变异是抗性马尾松变异的主要来源。4个亚群间FST为0.143,表明亚群间存在中等程度遗传分化。

      表 6  基于微卫星数据的遗传变异分析

      Table 6.  Analysis of molecular variance from microsatellite data using GenALEx

      变异
      来源
      自由度方差
      总和
      平均
      方差
      方差
      分量
      变异百
      分比/%
      P
      亚群间3184.86561.6220.95714<0.01
      亚群内1101 216.32811.0585.32380<0.01
      FST0.143
        说明:FST为遗传分化系数。
    • 利用M策略和随机取样策略构建核心种质的遗传多样性指标见表7。结果显示:基于随机取样策略的Power Core和基于M策略的Core Finder分别筛选得到的核心种质均保留了原有种质100%的等位基因数。但Core Finder核心种质包含的种质(72份)少于Power Core (79份),且两者间各遗传多样性参数差异不显著(P>0.05)。综合以上分析表明:选择M策略抽取72份抗性马尾松种质能够以最小的种质份数最大程度地代表收集区内抗性马尾松种质资源的遗传多样性。利用M策略构建的核心种质中包含广西2份(100.0%)、贵州1份(100.0%)、福建4份(66.7%)、安徽33份(其中4份为皖马抗良种)(78.6%)、江西8份(30.8%)和浙江24份(64.8%)种质材料。

      表 7  不同策略构建的核心种质遗传多样性指标比较

      Table 7.  Comparisons of genetic diversity index of core collections constructed by different tactics

      构建方法种质数/份NaNeIHoHeFNmPIC
      原有种质  1141155.541.510.060.640.910.660.90
      M策略   721155.301.430.060.630.910.610.92
      随机取样策略791155.531.480.060.640.910.650.92
        说明:Na. 等位基因数;Ne. 有效等位基因数;Ho. 观测杂合度;He. 期望杂合度;PIC. 多态性信息含量;I. Shannon’s多样性指数;Fis. 居群近交系数;F. 固定指数;Nm. 基因流。* P<0.05, **P<0.01。

      表8可见:构建的72份抗性马尾松核心种质占114份原有种质的63.16%,其Na均值、Ne均值、I均值分别占原有种质相应遗传参数的100.00%、95.67%和94.96%。t检验结果表明:列入研究的抗性马尾松种质与入选的核心种质的遗传多样性无显著差异(P>0.05),说明本研究基于M策略构建的抗性马尾松核心种质能充分代表列入研究的原抗性马尾松种质资源的遗传多样性,剔除了与114份原有种质重复或亲缘关系较近的材料。

      表 8  72份抗性马尾松核心种质与114份原有种质遗传多样性对比

      Table 8.  Comparison of genetic diversity between 72 core germplasm and 114 original germplasm of resistant P. massoniana

      种质种质数/份NaNeIHoHeFNmPIC
      原有种质 114 25.29 5.54 1.51 0.06 0.64 0.91 0.66 0.90
      核心种质 72 25.29 5.30 1.43 0.06 0.63 0.91 0.61 0.92
      保留比例/% 63.16 100.00 95.67 94.96 100.00 98.12 99.57 92.54 100.00
      t 0.000 −0.652 −1.428 0.059 −0.826 −0.052 0.014 −0.152
      P 1.000 0.525 0.177 0.953 0.424 0.960 0.989 0.882
        说明:Na. 等位基因数;Ne. 有效等位基因数;Ho. 观测杂合度;He. 期望杂合度;PIC. 多态信息含量;I. Shannon’s多样性指数;Fis. 居群近交系数;F. 固定指数;Nm. 基因流。

      PCoA分析(图5)显示:抗性马尾松核心种质均匀分布在整个马尾松种质资源中,表明本研究所构建的抗性马尾松核心种质具有很好的代表性。UPGMA进化树(图6)显示:72份抗性马尾松可分为3个类群,其中类群Ⅰ的11份抗性马尾松种质全部位于亚群Ⅰ(91.7%),类群Ⅱ的30份抗性马尾松种质全部位于亚群Ⅳ(61.2%),类群Ⅲ的抗性马尾松种质18份位于亚群Ⅱ(75.0%),其余13份位于亚群Ⅲ(44.8%)。结合结构分析可知:位于亚群Ⅱ中混杂有黄色和红色的核心种质与亚群Ⅲ中的核心种质聚为一个亚类,亚群Ⅱ中以黄色为主的核心种质聚为另一个亚类,但抗性马尾松核心种质的UPGMA聚类结果与列入研究的原有种质结构分析和PCoA分析结果基本一致,进一步对核心种质进行了确认。

      图  5  72份核心种质与114份原有种质主坐标分析     

      Figure 5.  Principal coordinate analysis of 72 core germplasm and 114 original germplasm

      图  6  72份抗性马尾松核心种质的UPGMA聚类图

      Figure 6.  UPGMA clustering map of 72 samples of resistant P. massoniana

    • 遗传多样性在生物育种工作中至关重要,是评价生物多样性的重要组成部分,物种的遗传多样性越高或遗传变异越丰富,其对环境变化的适应能力也就越强[26]。本研究利用14对 SSR 引物对经高强度人工接种松材线虫后存活下来的103份马尾松种质和安徽省林业科学研究院审定的11份高抗良种进行遗传多样性和群体结构分析。Na是衡量SSR位点多态性高低的重要指标,其中所有种质的平均Na为8.17,高于沈敬理等[27]利用SSR分子标记技术对131个马尾松无性系的研究(Na=3.67),可能是本研究选择的 SSR 引物具有较高多态性和基因分型检测手段不一样造成的。NeHeI等是评价遗传多样性的重要指标。高景斌等[28]利用 SCAR 分子标记对抗性马尾松材料进行遗传多样性研究,其NeI分别为1.51和0.46;YANG等[19]运用20对 SSR 引物分析马尾松二代亲本种质的遗传多样性,其Ne为3.01, I为1.26,PIC为0.89;董虹妤等[18]对选自马尾松第2代育种群体的12份马尾松材料运用13对 SSR 引物进行研究,其中NaI分别为3.17和0.47。与上述研究结果相比,本研究所选择的14对 SSR 引物能有效用于检测马尾松的遗传多样性(Ne=5.54、He=0.64、I=1.51、PIC=0.90),同时表明筛选获得的103份抗性马尾松种质和安徽省林业科学研究院审定的11份良种具有丰富的遗传多样性,在未来的抗性育种进程中具有很大的遗传增益潜力。

      Structure结构分析、UPGMA聚类分析和PCoA分析被广泛应用于推断植物的群体结构,但三者的统计原理不同。PCoA分析是一种非约束性的数据降维方法,通过对特征值和特征向量进行排序来寻找主坐标,从而有效地绘制所分析的样本;Structure结构分析基于贝叶斯模型的计算方法,根据基因型信息对多组样本进行分类;UPGMA聚类分析利用非加权组平均法对参试样本进行层次聚类。PCoA和UPGMA聚类分析可直观看到材料间的分类关系,但Structure结构分析能清晰地反映材料间的遗传背景和基因交流情况[2930]。总体而言,Structure结构分析、UPGMA聚类分析和PCoA分析相辅相成,两两使用有助于全面了解不同样本间的遗传差异。本研究Structure结构分析结果表明:114份抗性马尾松种质资源被分成4个亚群,而PCoA分析将所有抗性马尾松种质划分为3组,将结构分析所得4个亚群中的亚群Ⅱ和亚群Ⅲ合并为一组,可能是因为这2个亚群中的部分抗性马尾松种质Q<0.6,存在一定频率的基因交流和渗透。在结构分析中,每个亚群内都包含来自不同地区的个体,地理来源一致的抗性马尾松种质并没有完全聚在一起,这可能是由于不同地区育种者之间亲本交换或频繁的人工引种,导致新环境中种质的遗传信息发生了改变。安徽省林业科学研究院审定的11份抗松材线虫病良种是中国所有的抗松材线虫病马尾松优良种质,已在安徽省松材线虫病疫区进行了大面积的林相改造和造林,在松材线虫病防控中发挥了重要作用。这11份抗松材线虫病良种主要分布在亚群Ⅲ中,仅与19份抗性马尾松种质聚为一个亚群,这表明本研究筛选获得的大部分抗性马尾松优良品系与审定的良种遗传关系较远,具有较低的遗传相似性。因此,本研究筛选获得的103份抗性马尾松种质与审定的11份良种遗传信息同质性程度较低,且具有较高的遗传多样性,可保障今后营造的抗性林在生态系统上具有较强的稳定性。

    • 传统的核心种质构建方法是对比不同种质资源在植物学、形态学等性状上的差异来剔除亲缘关系相近的样本并且抽取核心种质,但表型性状容易受到外界因素影响而表现不稳定[31]。基于分子标记数据构建核心种质可减少温度、气候、环境等外界因素的影响,能够很好地作为植物形态和生理特性的有效补充。本研究根据SSR分子标记基因分型数据,利用Core Finder和Power Core,基于M策略和随机取样策略构建抗松材线虫病马尾松核心种质。综合对比分析不同构建策略,得到72份抗性马尾松核心种质,其保留了列入研究的原有种质全部的等位基因数,且各遗传多样性参数如NeIHePIC分别为5.30、1.43、0.63、0.92,均与114份原有种质的遗传多样性参数无显著差异(P>0.05)。本研究构建的72份抗性马尾松核心种质在广西、贵州、福建、安徽、江西和浙江6省中均有分布,但江西抗性马尾松核心种质抽取比例为30.8%,低于其他省份抗性马尾松种质的入选比例,表明该省抗性马尾松样本存在较高的遗传相似性。结合结构分析可知:72份抗性马尾松核心种质全面涵盖了列入研究的抗性马尾松种质的4个亚群,亚群Ⅲ抗性马尾松核心种质入选比例最低(44.8%),表明该亚群抗性马尾松种质遗传信息比较一致;亚群Ⅰ的抗性马尾松核心种质入选比例最高(91.7%),表明该亚群可能存在最丰富的遗传多样性。

      取样比例的大小是核心种质库能否构建成功的重要因素。WANG等[32]和ESCEIBANO等[33]认为核心种质取样比例为5%~40%,大多数为5%~15%。李魁鹏等[34]利用SSR分子标记从864份杉木Cunninghamia lanceolata中筛选获得了80份核心种质,占原有种质的9.3%;黄小凤等[35]利用SNP分子标记从221份荔枝Litchi chinensis品种(品系)中获得了44份核心种质,占原有种质的20%;张馨芳等[36]基于SSR分子标记数据,从161份板栗Yanshan chestnut中构建了含46份种质的核心种质,取样比例为28.6%。但本研究核心种质取样比例为63.13%,高于其他木本植物的取样比例。李自超等[37]认为核心种质取样比例通常由原有种质的数量、群体结构和遗传多样性水平等决定,对于遗传多样性丰富、原有种质数量较小的群体,可适当增大取样比例。

    • 利用14对 SSR 引物对114份抗松材线虫病马尾松种质的遗传多样性分析表明:所选抗性马尾松种质具有丰富的遗传多样性(Na=8.17,Ne=5.54,He=0.64,I=1.51,PIC=0.90),且本研究筛选的103份抗性马尾松材料遗传结构简单、遗传变异丰富,可用于开展马尾松重要性状的关联分析研究。根据SSR分子标记数据基于M策略构建了含有72份材料的抗性马尾松核心种质,保留了列入研究的原有种质100%的等位基因。经过t检验、PCoA分析和UPGMA聚类分析,构建的抗性马尾松核心种质与列入研究的原有种质之间无显著差异,剔除了遗传相近的材料,具有良好的代表性。

参考文献 (37)

目录

    /

    返回文章
    返回