SSR marker development, genetic diversity and population structure analysis in oil tree species Vernicia montana
-
摘要:目的 千年桐是大戟科油桐属的一种重要的工业油料树种，较同属的油桐具有更强的抗枯萎病能力，近年来受到了广泛关注。然而，因千年桐栽培历史短及种质管理的不足，目前对其分子遗传方面的研究还非常有限。本研究旨在开发千年桐的基因组SSR标记，进而开展其种质的鉴定、遗传多样性与群体结构分析。方法 采用RAD测序技术获得千年桐简化基因组以开发SSR标记，基于SSR标记利用分子变异分析（AMOVA）、非加权组平均法（UPGMA）聚类、主坐标（PCoA）与群体结构分析等方法对来自3个不同地理分布的共105份千年桐种质资源进行研究。结果 17个多态性的三核苷酸基因组SSR标记被开发，并能够很好地鉴别所有收集的105份千年桐种质。在其中62份种质中检测到85个私有等位基因，涉及15个SSR位点。AMOVA分析发现，千年桐群体间呈现出中等程度的遗传分化，但群体内的遗传变异远高于群体间的遗传变异。群体结构分析显示，3个来自不同地理分布的千年桐群体中存在4个不同的基因库，群体间既有进化独立性，又有较高程度的遗传混合，这一结果与UPGMA和PCoA分析的结果基本一致。结论 新开发的17个SSR标记有效鉴定了105份千年桐种质，揭示了其遗传多样性和群体遗传结构，对千年桐种质保存和育种计划具有非常重要的参考价值。Abstract:Objective Vernicia montana, belonging to the genus Vernicia of Euphorbiaceae family, is an important industrial oil tree species, which has received great attention in recent years because it is much more resistant to wilt disease than V. fordii of same genus. However, the current researches on the molecular genetics of V. montana are still very limited due to its short cultivation history and inadequate management of germplasm resources. This study aimed to develop V. montana’s genome SSR markers, and then to carry out its germplasm identification, genetic diversity and population structure analyses.Method Using RAD sequencing technology to obtain a simplified genome of V. montana to develop its SSR markers, based on the SSR markers, a series of methods including the analysis of molecular variance (AMOVA), the unweighted pair group method with arithmetic mean (UPGMA) clustering, principal coordinate analysis (PCoA) and population structure analysis were used to study 105 germplasm resources of V. montana from three different geographical distributions.Result 17 polymorphic trinucleotide genomic-SSR markers were developed, which well distinguished all the tested 105 V. montana germplasms. A total of 85 private alleles were detected in 62 germplasms at 15 SSR loci. The AMOVA analysis revealed a moderate degree of genetic differentiation among the populations of V. montana; however, the genetic variation within the populations was much higher than that among the populations. Population structure analysis showed that four different gene pools were present in the three V. montana populations from different geographical distributions, and there was both evolutionary independence and a relatively high degree of genetic admixture among the populations, which was basically consistent with the results of UPGMA and PCoA analyses.Conclusion The newly developed 17 SSR markers effectively identified 105 V. montana germplasms, and revealed the genetic diversity and population genetic structure of the germplasms, which would be very helpful for the conservation of germplasms and breeding program in V. montana.
-
Keywords:
- Vernicia montana /
- SSR marker /
- germplasm identification /
- genetic diversity /
- population structure
-
油桐（Vernicia fordii）和千年桐（V. montana）是大戟科（Euphorbiaceae）油桐属（Vernicia）中两个重要的工业木本油料树种，主要分布在中国南部和一些东南亚国家[1-2]。这两个树种的种子油富含共轭三烯十八碳不饱和脂肪酸即α-硬脂酸，是著名的工业干燥油[3]，已被广泛用于油漆、油墨等工业用料和生物柴油的生产[2, 4-5]。与油桐相比，千年桐的栽培和研究历史相对较短，但有研究表明其种子油产量高于油桐[1]。此外，千年桐对于土壤传播的镰刀菌枯萎病的抵抗力要显著强于油桐。在栽培中发现，枯萎病会严重阻碍油桐的生长与发育，而以千年桐作砧木，将油桐嫁接到千年桐上可以很好地抵抗枯萎病[6]。这些研究表明，千年桐具有更好的开发和利用价值。然而，千年桐栽培和发展历史较短，导致对其种质的管理不完善[1, 7]，相关种质鉴定和遗传多样性的研究还未充分开展。因此，本研究开展千年桐的种质鉴定和遗传多样性分析对于其种质保存、高产和抗病的新品种培育具有重要的实践意义。
近年来，分子标记已被证明是鉴定植物种质和揭示种质遗传多样性的有效工具[8-9]。许多类型的分子标记已被开发并广泛使用，包括限制性酶切片段长度多态性标记（restriction fragment length polymorphism，RFLP），扩增片段长度多态性标记（amplified fragment length polymorphism，AFLP），简单重复序列标记（simple sequence repeat，SSR），简单重复序列区间标记（inter-simple sequence repeat，ISSR）和单核苷酸多态性标记（single nucleotide polymorphism，SNP）[8, 10-11]。其中，SSR标记由于在基因组中的丰富性、高多态性、可再现性和共显性遗传等特点，已成为最优选的分子标记类型之一[12-13]，已成功应用于油桐、油茶（Camellia oleifera）和油橄榄（Olea europaea）等多种木本油料树种的种质鉴定、遗传多样性和群体结构分析等研究[14-18]。然而，SSR标记的应用在千年桐中非常有限。迄今为止，仅一篇关于千年桐SSR标记开发和应用的研究报道，且其开发的基因组SSR标记仅基于二核苷酸重复序列[3]。与较长的SSR重复序列类型相比，二核苷酸SSR标记在识别大小接近的等位基因时，不易检测到其差异，从而容易导致错误的判断或产生混淆[19-20]。因此，具有较长重复核心的SSR标记将更有利于千年桐的种质鉴定和遗传分析。
针对上述问题，本研究基于限制性位点关联DNA（restriction-site associated DNA，RAD）测序技术开发多态性的三核苷酸基因组SSR标记，对广西省收集的105份千年桐种质资源进行基因分型。开发的SSR标记在群体水平上进一步研究105份千年桐种质资源的遗传多样性和结构。
1. 材料与方法
1.1 试验材料
本研究从千年桐种质资源最丰富的中国广西壮族自治区收集到105份千年桐种质资源。根据其地理位置的分布情况，将这105份千年桐种质资源划分为3个群体。其中群体1分布于河池市凤山县（27份），群体2分布于百色市田林县（34份），群体3分布于南宁市西乡塘区（44份）。其中群体3中的“桂皱1号”“桂皱2号”“桂皱6号”和“桂皱27号”4个高产优良无性系是由广西林业科学研究院选育出的优良品种[21]。为了选择具有代表性的种质，所收集的千年桐种质资源之间的间隔至少为100 m，每一份种质由一棵树代表。
1.2 DNA提取和测序
利用植物DNA提取试剂盒（DP320-03，天根生化科技有限公司，北京）提取105份千年桐样品的总DNA，经质检合格后基于RAD测序技术构建测序文库，在Illumina平台进行高通量测序（ORI-GENE，中国北京）。高通量测序得到的原始数据（raw reads）过滤掉包含接头、poly N和poly A等低质量的序列后得到高质量的过滤后数据（cleaned reads），使用Trinity和FLASH软件将过滤后的短序列拼接得到重叠群（contigs），然后进一步组装成简化的基因组。
1.3 SSR分析
使用基于Perl脚本的MIcroSAtellite鉴定工具（MISA）扫描千年桐RAD测序文库以鉴别SSR位点，扫描的标准是三核苷酸基序至少5个重复，四和五核苷酸基序至少4个重复。SSR引物通过Primer Premier 5软件设计，由睿博兴科生物公司合成。PCR反应体系的总体积为20 μL，其中包含40 ~ 60 ng基因组DNA，25 μmol/L的dNTP，2.5 U的Taq DNA聚合酶，10 μmol/L的正、反向引物，10 μmol/L的荧光染料，及含有25 mmol/L Mg2+的10 × PCR缓冲液（天根生化科技有限公司，北京）。PCR反应程序采用两步法：94 ℃预变性5 min，然后94 ℃变性30 s，55 ℃退火40 s，72 ℃延伸30 s，共30个循环，接着进行94 ℃变性30 s，53 ℃退火40 s，72 ℃延伸30 s，共8个循环，最后72 ℃完全反应10 min。通过毛细管电泳检测具有荧光标记的PCR产物。每个试验进行3次重复。
1.4 数据分析
通过GeneMarker v2.2.0软件分析从毛细管电泳获得的SSR原始数据。使用POPGENE 32和Cervus v3.0.7软件计算每个SSR位点的遗传多样性信息参数，包括观察到的等位基因（Na）、有效等位基因（Ne），观察到的杂合度（Ho）、预期的杂合度（He）、香农信息指数（I）、基因流（Nm）、遗传分化系数（FST）、无效等位基因频率（FNull）和多态性信息含量（PIC）。使用GenALEx v6.503软件进行哈迪−温伯格平衡检测（HWE），同时检测105份千年桐种质资源中的私有等位基因。
GenALEx v6.503软件也用于AMOVA分子变异分析，同时AMOVA中的FST用于评估群体遗传分化。使用PowerMarker v3.25软件计算Nei’s遗传距离[22]，基于遗传距离，使用非加权组平均法（UPGMA）对105份千年桐种质进行聚类分析。为了获得更直观的遗传分类，利用iTOL在线工具进一步注释系统发育树。同时，在GenALEx v6.6.503软件中使用上述遗传距离数据矩阵进行主坐标分析（PCoA）。使用STRUCTURE v2.3软件对105份千年桐种质资源进行群体遗传结构分析。在标准STRUCTURE分析中，假定的种群数量用K值表示，设置其范围为1 ~ 8，对8个K值进行逐一测试，并针对每个K值执行5次独立运行。通过Structure Harvester程序使用Delta−K法评估最佳种群数量的K值。
2. 结果与分析
2.1 RAD测序与SSR位点的鉴定
本研究运用RAD测序技术共得到45.932 Mbp的原始数据，去除低质量的序列后获得30.996 Mbp的过滤后数据，占原始数据的67.5%（表1）。Q20和Q30分别为98.3%和95.4%，GC含量为40.51%（表1），这表明测序质量很高，数据可靠。过滤后的短序列经过拼接获得了139 988条重叠群，其平均长度为256 bp（表1），然后进一步组装成简化基因组。在千年桐简化基因组中共鉴定到3 425个SSR位点（表2），其中包括2 174个二核苷酸SSR位点，1 023个三核苷酸SSR位点，131个四核苷酸SSR位点，及97个五核苷酸SSR位点，分别占总SSR位点的63.47%、29.87%、3.83%和2.83%。除二核苷酸SSR位点外，剩余1 251个三、四、五核苷酸SSR位点均被用于基因组SSR标记的开发。
表 1 测序数据质量统计总表Table 1. Quality summary of sequencing data原始数据
Raw read/Mbp原始Q20
Raw Q20/Gbp原始Q30
Raw Q30/Gbp过滤后数据
Cleaned read/Mbp过滤后Q20
Cleaned Q20/Gbp过滤后Q30
Cleaned Q30/Gbp平均长度
Average length/bpGC含量
GC content/%45.932 6.580 (95.5%) 6.257 (90.8%) 30.996 (67.5%) 4.510 (98.3%) 4.374 (95.4%) 148.0 40.51 注：Q20. 质量值 ≥ 20的碱基所占百分比；Q30. 质量值 ≥ 30的碱基所占百分比。Notes: Q20, percentage of bases with the quanlity value ≥ 20; Q30, percentage of bases with the quanlity value ≥ 30. 表 2 具有不同重复基序的SSR位点在千年桐基因组中的分布Table 2. Distribution of SSR with different repeat motifs in the genome of Vernicia montanaSSR重复基序
SSR repeat motifSSR数量
Number of SSR频率
Frequency/%二核苷酸 Dinucleotide 2 174 63.47 三核苷酸 Trinucleotide 1 023 29.87 四核苷酸 Tetranucleotide 131 3.83 五核苷酸 Pentanucleotide 97 2.83 总计 Total 3 425 100.00 2.2 SSR标记的开发
针对上述1 251个三、四、五核苷酸SSR位点，共设计得到300对SSR引物，其中包括250对三核苷酸SSR引物、28对四核苷酸SSR引物和22对五核苷酸SSR引物。利用这300对引物对随机选取的一份千年桐种质资源（“桂皱1号”）的基因组DNA进行PCR扩增，毛细管电泳检测发现其中173对SSR引物扩增出了预期片段大小的PCR产物，具体包括162对、5对和6对三核苷酸、四核苷酸和五核苷酸SSR引物。这173对SSR引物被进一步在8份随机选取的千年桐种质资源中进行基因分型，筛选出56对多态性引物，包括52对三核苷酸SSR引物和4对五核苷酸SSR引物。随后，56对SSR引物被用于所有105份千年桐种质资源的基因分型，最终获得了17个能够完全鉴别105份千年桐种质资源的多态性SSR标记（表3）。
表 3 针对105份千年桐种质开发的17个多态性SSR标记的信息Table 3. Information of 17 polymorphic SSR markers developed for 105 V. montana germplasms引物编号
Primer No.重复基序
Repeat motif引物序列(5′—3′)
Primer sequence (5′−3′)退火温度
Annealing temperature/℃预期大小
Expected size/bpVm-BFU01 (TTG)5 AACCACTGCTACTTCACCATTTTC 55 167 ACCCAATGTTTTCTCCGACC Vm-BFU02 (AGA)5 GTGTTGAGCCAGAAACCCATTA 55 165 CAGAGAAGCCTCGGTCCCTA Vm-BFU03 (CAT)6 TCTACTCCCACACTTCCAAAACA 55 161 TAATCTCCTTCTGTCCTTCACGA Vm-BFU04 (CAG)9 GGGAGTGGTGGCAATGGC 55 181 GCTGGGAGGCATTGTTGAAG Vm-BFU05 (AGA)9 GAGCCAAGAGAAGACGAAAAGAG 55 146 ACCGTTTACAGTGTTTCGCTATG Vm-BFU06 (GGA)8 TGTGCCGCTTGTAACTGCC 55 164 TGCGGCTGTGTCAGGTGTAG Vm-BFU07 (TCT)7 AGCCTTTGCCACTGTTGAGC 55 127 GATGGGTCCGCCAAGTTCA Vm-BFU08 (TCT)5 ATTGTGAAGGATTTGCGATGG 55 161 CGGCGAAAACGAAACAGAG Vm-BFU09 (GTT)7 TCGCCTAAGGTGGTCTTGATG 55 178 GCCCAACGAAATCTAACTCTAATAA Vm-BFU10 (TAA)12 TTCCTCCTCTGGTGACGCTT 55 252 TTCCTTCCATCATCAACTTTTACC Vm-BFU11 (TCC)7 GCCGCCGCCTACTACTTACTT 55 131 TTTCTCAAACCAAACAGGAGTTG Vm-BFU12 (CTT)5 TATTTTTCTTGGGAGTAAAGTCACC 55 198 TATGTGAAATGGAGAGTTCGGAG Vm-BFU13 (CAG)6 TTGTCAACAAGCCTTCTCACCT 55 160 GCTCCAAGTCCCATCATCATTT Vm-BFU14 (GGC)5 CCAAAACCATCAATCTCTTCGC 55 204 TGATTTCGCACAAGTCCCAAG Vm-BFU15 (ACC)10 GCGTTCCTGACCCTACCTTT 55 184 AGAGAAACAAAAAAGCCACCAG Vm-BFU16 (TTA)12 TCCCCACAGCCATAAAACAAG 55 199 TTTCCAAAACTCTCAAACCACAA Vm-BFU17 (AAT)10 ACCCCATCTATGACATCCCACT 55 195 CCCCGTTCTTGCTCTCCC 注：重复基序括号外的数字为其重复次数。Note: the number of repeats is marked outside the parentheses of repeat motifs. 2.3 SSR标记的遗传多样性评价
对于所开发的17个SSR标记，PIC的值在0.309 ~ 0.824之间，平均值为0.552，其中10个为高度多态性标记（PIC ≥ 0.50）、7个为中度多态性标记（0.25 ≤ PIC < 0.50）。基于17个SSR标记，在105份千年桐种质中共检测到138个等位基因，每个SSR位点Na最少为2个，最多达14个，平均为8.117个（表4）。值得注意的是，遗传多样性评价的结果显示了私有等位基因（number of private allele，NPA）的存在。在17个多态性SSR位点中的15个位点上共检测到85个私有等位基因，平均每个位点有5个私有等位基因（表4和表5）。这85个私有等位基因存在于62份千年桐种质中，其中7份种质来自于群体1，22份种质来自于群体2，33份种质自于群体3（表5）。此外，其他遗传多样性评价参数，如Ho、I及Nm等与PIC呈现正相关的关系（表4）。
表 4 17个千年桐SSR位点的遗传多样性参数Table 4. Genetic diversity parameters of 17 SSR loci in V. montana位点 Locus Na Ne Ho He I FNull PIC FST Nm HWE NPA Vm-BFU01 2.000 1.618 0.190 0.384 0.570 0.335 0.309 0.147 1.446 *** 0 Vm-BFU02 6.000 1.491 0.086 0.331 0.704 0.590 0.310 0.012 20.486 *** 1 Vm-BFU03 10.000 2.400 0.365 0.586 1.223 0.218 0.539 0.055 4.286 *** 9 Vm-BFU04 11.000 3.262 0.476 0.697 1.620 0.178 0.670 0.064 3.629 *** 3 Vm-BFU05 4.000 2.187 0.286 0.546 0.923 0.306 0.466 0.018 13.956 *** 1 Vm-BFU06 5.000 2.595 0.417 0.618 1.075 0.179 0.547 0.018 13.645 *** 1 Vm-BFU07 6.000 3.135 0.419 0.684 1.357 0.232 0.638 0.035 6.850 *** 2 Vm-BFU08 8.000 2.191 0.583 0.546 1.160 0.054 0.510 0.034 7.124 *** 2 Vm-BFU09 4.000 2.237 0.390 0.556 0.908 0.162 0.453 0.072 3.214 *** 6 Vm-BFU10 9.000 2.078 0.200 0.521 1.149 0.450 0.491 0.063 3.695 *** 12 Vm-BFU11 13.000 6.289 0.235 0.845 2.101 0.559 0.824 0.057 4.163 *** 3 Vm-BFU12 2.000 1.873 0.320 0.469 0.659 0.186 0.358 0.050 4.742 ** 0 Vm-BFU13 14.000 3.157 0.410 0.687 1.617 0.252 0.657 0.093 2.439 *** 10 Vm-BFU14 10.000 1.886 0.248 0.472 1.144 0.304 0.457 0.097 2.338 *** 15 Vm-BFU15 11.000 3.155 0.352 0.686 1.544 0.333 0.648 0.062 3.808 *** 10 Vm-BFU16 9.000 5.253 0.548 0.814 1.787 0.193 0.782 0.056 4.256 *** 2 Vm-BFU17 14.000 3.913 0.538 0.748 1.852 0.158 0.725 0.031 7.705 *** 8 总数 Total 138.000 48.728 85 平均值 Mean 8.117 2.866 0.357 0.599 1.258 0.276 0.552 0.056 4.251 5 注：Na. 观测到的等位基因数；Ne. 有效等位基因数；Ho. 观测到的杂合度；He. 预期杂合度；I. 香农信息指数；Nm. 基因流；FNull. 无效等位基因频率；PIC. 多态信息含量；FST. 遗传分化指数；HWE. 哈迪−温伯格平衡；*. P < 0.05；**. P < 0.01； ***. P < 0.001；NPA. 私有等位基因的数量。下同。Notes: Na, observed number of alleles; Ne, effective number of alleles; Ho, observed heterozygosity; He, expected heterozygosity; I, Shannon’s information index; Nm, gene flow; FNull, null allele frequency; PIC, polymorphism information content; FST, genetic differentiation index; HWE, Hardy-Weinberg equilibrium; * means P < 0.05; ** means P < 0.01; *** means P < 0.001; NPA, number of private alleles. The same below. 表 5 具有一至多个私有等位基因的千年桐种质Table 5. V. montana germplasms with one or more private alleles种质
Germplasm群体
PopulationNPA 私有等位基因位点
Locus of private allele种质
Germplasm群体
PopulationNPA 私有等位基因位点
Locus of private allele1-3 1 1 Vm-BFU11 3-3 3 2 Vm-BFU03, Vm-BFU17 1-4 1 2 Vm-BFU03, Vm-BFU07 3-4 3 3 Vm-BFU03, Vm-BFU14, Vm-BFU16 1-7 1 2 Vm-BFU03, Vm-BFU07 3-5 3 1 Vm-BFU14 1-15 1 1 Vm-BFU06 3-6 3 1 Vm-BFU13 1-18 1 1 Vm-BFU11 3-9 3 3 Vm-BFU03, Vm-BFU11, Vm-BFU13 1-23 1 1 Vm-BFU13 3-10 3 1 Vm-BFU15 1-27 1 2 Vm-BFU10, Vm-BFU17 3-12 3 1 Vm-BFU14 2-2 2 1 Vm-BFU10 3-13 3 1 Vm-BFU15 2-3 2 1 Vm-BFU09 3-14 3 2 Vm-BFU14, Vm-BFU17 2-4 2 1 Vm-BFU10 3-15 3 1 Vm-BFU13 2-5 2 3 Vm-BFU09, Vm-BFU10, Vm-BFU13 3-16 3 1 Vm-BFU14 2-6 2 2 Vm-BFU03, Vm-BFU13 3-17 3 1 Vm-BFU14 2-11 2 1 Vm-BFU10 3-22 3 1 Vm-BFU03 2-13 2 1 Vm-BFU15 3-23 3 2 Vm-BFU04, Vm-BFU14 2-14 2 2 Vm-BFU10, Vm-BFU15 3-24 3 1 Vm-BFU15 2-15 2 1 Vm-BFU09 3-25 3 3 Vm-BFU04, Vm-BFU14, Vm-BFU15 2-16 2 1 Vm-BFU09 3-26 3 1 Vm-BFU14 2-17 2 1 Vm-BFU09 3-27 3 2 Vm-BFU15, Vm-BFU17 2-19 2 1 Vm-BFU15 3-28 3 2 Vm-BFU08, Vm-BFU17 2-20 2 1 Vm-BFU10 3-29 3 1 Vm-BFU14 2-23 2 1 Vm-BFU10 3-30 3 1 Vm-BFU14 2-24 2 1 Vm-BFU09 3-31 3 1 Vm-BFU14 2-26 2 1 Vm-BFU02 3-32 3 2 Vm-BFU03, Vm-BFU05 2-27 2 2 Vm-BFU10, Vm-BFU15 3-33 3 1 Vm-BFU16 2-28 2 1 Vm-BFU13 3-35 3 1 Vm-BFU14 2-29 2 1 Vm-BFU13 10-7 3 1 Vm-BFU17 2-30 2 2 Vm-BFU10, Vm-BFU13 10-13 3 2 Vm-BFU03, Vm-BFU15 2-31 2 1 Vm-BFU13 10-16 3 1 Vm-BFU14 2-33 2 1 Vm-BFU10 10-18 3 1 Vm-BFU10 3-1 3 1 Vm-BFU17 Guizhou-2 3 2 Vm-BFU08, Vm-BFU17 3-2 3 1 Vm-BFU14 Guizhou-27 3 1 Vm-BFU04 总计 Total 85 2.4 种质的遗传分化分析
遗传分化指数（FST）是评估群体间遗传分化程度的参数。本研究利用POPGENE软件分析得出：17个SSR标记的FST的范围在0.012 ~ 0.147之间，平均值为0.056（表4）。利用GenALEx软件计算得出的FST值为0.080（P < 0.001）（表6），与POPGENE的结果基本一致。基于所开发的17个SSR标记，AMOVA分子变异分析进一步显示，千年桐群体内的遗传变异远高于群体间的遗传变异，其群体内与群体间的遗传变异分别占总变异的92%和8%。
表 6 千年桐群体间与群体内的AMOVA分析Table 6. AMOVA analysis among and within V. montana populations变异来源
Source of variation自由度
Degree of freedom方差和
Sum of squares变异组分
Variance component占总变异的比例
Proportion in total variation/%FST 群体间 Among populations 2 110.841 1.209 8 0.080 (P < 0.001) 群体内 Within population 102 1 423.025 13.951 92 总计 Total 104 1 533.867 15.160 100 基于SSR位点的哈迪−温伯格平衡能够评估群体是否存在连锁不平衡，判断是否偏离遗传平衡，具体用显著性水平P−value来表示。当P > 0.05时，群体符合遗传平衡；而当P < 0.05时，群体显著偏离遗传平衡[23]。本研究中，17个多态性SSR位点经Bonferroni校正后显示，千年桐群体全部显著偏离哈迪−温伯格平衡（P < 0.05）（表4）。
2.5 种质的UPGMA聚类与PCoA分析
基于Nei’s遗传距离、利用UPGMA方法对来自3个群体的105份千年桐种质进行聚类分析，结果显示（图1）：群体1、2和3的部分种质（约占全部种质的三分之一）分别单独聚集在Ⅰ、Ⅱ和Ⅲ这3个分支中；来自3个群体的部分种质（约占全部种质的三分之一）共同聚集在较大的分支Ⅳ中；来自群体1和3的部分种质（不到种质总数的三分之一）共同聚集在分支Ⅴ中。
![]()
图 1 基于17个SSR标记的105份千年桐种质资源遗传关系的进化树红色、黄色和蓝色的线段分别代表群体1、2和3。The red, yellow and blue lines represent the population 1, 2 and 3, respectively.Figure 1. A dendrogram of genetic relationship among 105 V. montana germplasm resources based on 17 SSR markers105份千年桐种质的PCoA分析结果显示，第一个主坐标和第二个主坐标分别占所有千年桐种质总遗传变异的14.07%和9.01%（图2）。PCoA分析结果与UPGMA聚类分析的结果基本一致。
![]()
图 2 基于17个SSR标记的105份千年桐种质资源的主坐标分析(PCoA)红色、黄色和蓝色的点分别代表群体1、2和3。The red, yellow and blue dots represent population 1, 2 and 3, respectively.Figure 2. Principal coordinate analysis (PCoA) of the 105 V. montana germplsm resources based on 17 SSR markers2.6 种质的群体结构分析
为了揭示105份千年桐种质的群体结构，基于贝叶斯的方法分析了来自不同地理分布的3个群体的遗传结构，其中推定的种群数量（基因库）用K值表示。利用标准STRUCTURE分析，对8个K值（K = 1 ~ 8）逐一测试，根据Delta K值确定了最佳种群数量为2（K = 2），其次为4（K = 4）（图3a、b）。考虑到标准STRUCTURE方法中软件可能存在的偏差，因此分别在K = 2和K = 4的情况下分析了来自3个群体的105份千年桐种质的遗传结构（图3c、d）。结果表明，K = 4的群体结构比K = 2的群体结构更合理，因为来自不同地理分布的3个群体的所有种质都以各自独特的基因库出现在K = 4而不是K = 2的群体结构中。如图3所示，来自不同地理分布的3个千年桐群体中存在4个基因库。群体1、2和3分别由红色、黄色和蓝色片段代表的基因库主导，而绿色片段代表的基因库散布在这3个群体中，其在群体2和群体3中出现的频率比在群体1中更高。因此，在这105份千年桐种质资源中，基于群体结构推定的种群接近于地理划分的种群。
![]()
图 3 基于17个SSR标记的105份千年桐种质资源的群体结构分析a. K = 1 ~ 8时Delta K的变化图；b. 基于每个K值的5个重复，使用ln P(D)的平均值进行种群估计；c. 105份千年桐种质在K = 2时的种群结构；d. 105份千年桐种质在K = 4时的种群结构。a, plot of Delta K with K ranging from 1 to 8; b, estimation of population using the mean of ln P(D) based on five repetitions for each K-value; c, population structure of 105 V. montana germplasms at K = 2; d, population structure of 105 V. montana germplasms at K = 4.Figure 3. Population structure analysis of 105 V. montana germplasm resources based on 17 SSR markers3. 讨论与结论
3.1 RAD测序技术应用的可行性
RAD测序是一种基于全基因组酶切位点的简化基因组测序技术。它可以大规模快速鉴定基因组中的SSR位点，并具有高通量、高精度、低成本和短周期的优势[24-25]。对于没有参考基因组的物种，RAD测序技术是一种开发基因组SSR标记的替代方法。近年来，通过RAD测序技术已在大麦（Hordeum vulgare）、甘薯（Dioscorea esculenta）等多种植物中成功开发了SSR标记[25-27]，证实了RAD测序技术的可行性。因此，本研究首次利用RAD测序技术获得了无参考基因组的千年桐的简化基因组，开发SSR标记来鉴定种质，并分析千年桐种质的遗传多样性和群体结构。
3.2 千年桐SSR标记开发及其遗传多样性
SSR标记是植物基因组中以1 ~ 6个碱基为核心单位组成的串联重复DNA序列[28]。与单核苷酸和双核苷酸SSR标记相比，具有较长重复核心的三核苷酸、四核苷酸和五核苷酸SSR标记更容易检测到等位基因长度上的差异，因而更利于种质的鉴别[28-29]。在以往的千年桐研究中，只有一项研究开发了20对二核苷酸SSR标记来鉴别30份千年桐种质[3]。在本研究中，通过对300个包括三、四和五核苷酸SSR标记的筛选，最终发现只有17个多态性的三核苷酸SSR标记可以很好地鉴别所有105份千年桐种质（表2和表3），说明三核苷酸SSR标记对于鉴别千年桐种质是最有效的。三核苷酸SSR标记在种质鉴定中的高效性也已在其他植物中得到证实，例如大豆（Glycine max）、猕猴桃（Actinidia chinensis）、鹰嘴豆（Cicer arietinum）和油橄榄[18, 30-32]。
PIC作为一种重要的遗传多样性参数，代表了SSR位点的变异程度，可以用来评估SSR标记的鉴别能力，且不受种质数量、地理起源和所调查的基因位点的影响[33-34]。在本研究中，所开发的17个基因组SSR标记的平均PIC值为0.552（表3），高于油桐中开发的22个EST-SSR标记的平均PIC值（0.401）[1]。基因组SSR标记的多态性高于由转录组开发的EST-SSR标记的多态性，这在核桃（Juglans regia）的研究中也有报道[35-36]，这种现象很可能是由于基因组中含有不保守的内含子所致。
本研究所开发的基因组SSR位点，除Vm-BFU08外，观察到的杂合度（Ho）均低于预期的杂合度（He）（表4）。在核桃和白柳（Salix alba）中也观察到了类似的现象[37-38]，这被认为是由于纯合子过多或无效等位基因存在所致[37, 39]。在千年桐的一些SSR位点中，如Vm-BFU01、Vm-BFU02，Vm-BFU05、Vm-BFU10和Vm-BFU15，出现了过多的纯合子，而Vm-BFU11和Vm-BFU16位点显示出较高的无效等位基因频率（表4）。因此，本研究中SSR位点Ho低于He也是由过多纯合子和无效等位基因两方面原因导致。对于无效等位基因的存在，被认为是引物结合位点上插入或缺失突变导致的SSR位点侧翼序列变异引起的[40-41]。
私有等位基因的存在可以在一定程度上反映植物种质资源的遗传多样性[42-43]，这对种质的鉴定和后续的育种工作具有重要的指导意义。在本研究中，占总种质一半以上的千年桐种质资源（62份）中发现含私有等位基因（表4和表5），这揭示了105份千年桐种质资源中存在相对广泛的遗传变异。在植物栽培和驯化过程中，私有等位基因的存在被认为是由基因交换或可能的突变引起的[18, 44]。值得注意的是，在这62份具有私有等位基因的种质资源中，有一半以上来自群体3（表5），这表明群体3的种质比群体1和2具有更广泛的遗传变异。
3.3 千年桐种质的遗传分化和群体结构分析
遗传分化指数（FST）可以评估群体间遗传分化的程度[44]。FST < 0.05、0.05 < FST < 0.15、0.15 < FST < 0.25和FST > 0.25分别表示弱、中、高和极高的群体遗传分化水平[45]。在本研究中，由POPGENE和AMOVA分析计算的FST值在0.012至0.147之间（表4和表6），表明所测千年桐群体间存在中度程度的遗传分化。然而，千年桐群体内的遗传变异远高于群体间（表6），这一现象在多种木本植物中也被发现，例如核桃[35]、油棕（Elaeis guineensis）[46]、白柳[38]、油桐[7]和杉木（Cunninghamia lanceolata）[13]。造成这种现象的原因可能与不同群体种质资源之间的基因交流频率增加有关[13, 46]。在群体遗传学中，通常采用哈迪−温伯格平衡（HWE）来检测群体是否存在连锁不平衡和群体遗传混合的程度[47]。在本研究中，基于17个SSR位点的3个千年桐群体明显偏离了哈迪−温伯格平衡（表4）。一般说来，群体偏离哈迪−温伯格平衡被认为可能是SSR位点的杂合度或纯合度过高所致[38]。在本研究中发现，17个SSR位点中的大部分位点的杂合度过高（表4），这可能是千年桐群体偏离哈迪−温伯格平衡的主要原因，同时也反映了一定程度的群体遗传混合。在自然条件下，隔离、选择、杂交和漂移等过程会影响群体遗传结构的形成，而这些过程也会受到人为干预的影响[48]，从而可能导致SSR位点的杂合度过高。
群体结构分析显示，来自不同地理分布的3个千年桐群体既有进化上的相对独立性又有较高程度的遗传混合（图3），这与UPGMA和PCoA的结果基本一致（图1和图2）。植物群体之间的遗传混合通常与地理位置或地理环境密切相关[38, 44]。在土耳其的2个河流生态系统中，广泛人为介导的植物材料迁移显著影响了白柳群体的遗传结构[38]。在本研究中，3个群体的千年桐种质资源分别采自广西省中部红水河流域自西向东的3个不同地点。红水河流域具有丰富的水资源和植被资源，且该河流生态系统中人类活动较为普遍[49-50]。因此，千年桐3个群体之间较高程度的遗传混合与红水河流域及周边地区的人类活动密切相关。在一些天然河流生态系统中，河流下游的植物群体的遗传多样性高于上游，比如当归（Angelica archangelica）、疏花水柏枝（Myricaria laxiflora）和柳树（Salix hukaoana）[51-53]。植物种子或营养繁殖体通过河流系统自上游向下游流动可以解释这种现象，也能解释红水河下游的群体3比红水河上游的群体1和2遗传混合高的现象（图1 ~ 3）。群体3中含私有等位基因的种质资源也最多，与上述现象有较好的相关性。而对于红水河上游的2个群体，群体2的遗传混合高于群体1（图3），这很可能是受到了人类活动的干扰。人类活动对植物种质资源遗传混合的影响常见于一些国内外分布较广的植物物种中，包括油茶[15]、椰子（Cocos nucifera）[54]、核桃[16, 35-36]、油棕[17]和油桐[7]。例如，与千年桐同属的油桐，人类携带其种子在产地间迁移等活动对其种质资源的遗传混合产生了显著影响[7]。值得注意的是，本研究中利用群体结构分析推定的群体数量（基因库）略多于从不同地理起源收集的群体数量（图3），类似的现象也在葡萄（Vitis vinifera）、油橄榄和石榴（Punica granatum）的群体遗传结构分析中被发现[9, 41, 55]。因此，本研究中的3个千年桐群体可能由于历史起源确实存在4个基因库，但也不排除其由于群体间的基因交换、基因突变或其他原因创建了1个新的基因库。
综上所述，本研究首次在千年桐中成功开发了17个多态性的三核苷酸基因组SSR标记。这17个SSR标记能很好地鉴别来自不同地理分布的3个千年桐群体的105份种质资源，并揭示了种质的遗传多样性和群体遗传结构，为千年桐的种质保存和育种提供了有价值的科学依据。
-
![]()
图 1 基于17个SSR标记的105份千年桐种质资源遗传关系的进化树
红色、黄色和蓝色的线段分别代表群体1、2和3。The red, yellow and blue lines represent the population 1, 2 and 3, respectively.
Figure 1. A dendrogram of genetic relationship among 105 V. montana germplasm resources based on 17 SSR markers
![]()
图 2 基于17个SSR标记的105份千年桐种质资源的主坐标分析(PCoA)
红色、黄色和蓝色的点分别代表群体1、2和3。The red, yellow and blue dots represent population 1, 2 and 3, respectively.
Figure 2. Principal coordinate analysis (PCoA) of the 105 V. montana germplsm resources based on 17 SSR markers
![]()
图 3 基于17个SSR标记的105份千年桐种质资源的群体结构分析
a. K = 1 ~ 8时Delta K的变化图；b. 基于每个K值的5个重复，使用ln P(D)的平均值进行种群估计；c. 105份千年桐种质在K = 2时的种群结构；d. 105份千年桐种质在K = 4时的种群结构。a, plot of Delta K with K ranging from 1 to 8; b, estimation of population using the mean of ln P(D) based on five repetitions for each K-value; c, population structure of 105 V. montana germplasms at K = 2; d, population structure of 105 V. montana germplasms at K = 4.
Figure 3. Population structure analysis of 105 V. montana germplasm resources based on 17 SSR markers
表 1 测序数据质量统计总表
Table 1 Quality summary of sequencing data
原始数据
Raw read/Mbp原始Q20
Raw Q20/Gbp原始Q30
Raw Q30/Gbp过滤后数据
Cleaned read/Mbp过滤后Q20
Cleaned Q20/Gbp过滤后Q30
Cleaned Q30/Gbp平均长度
Average length/bpGC含量
GC content/%45.932 6.580 (95.5%) 6.257 (90.8%) 30.996 (67.5%) 4.510 (98.3%) 4.374 (95.4%) 148.0 40.51 注：Q20. 质量值 ≥ 20的碱基所占百分比；Q30. 质量值 ≥ 30的碱基所占百分比。Notes: Q20, percentage of bases with the quanlity value ≥ 20; Q30, percentage of bases with the quanlity value ≥ 30. 
下载: 导出CSV
表 2 具有不同重复基序的SSR位点在千年桐基因组中的分布
Table 2 Distribution of SSR with different repeat motifs in the genome of Vernicia montana
SSR重复基序
SSR repeat motifSSR数量
Number of SSR频率
Frequency/%二核苷酸 Dinucleotide 2 174 63.47 三核苷酸 Trinucleotide 1 023 29.87 四核苷酸 Tetranucleotide 131 3.83 五核苷酸 Pentanucleotide 97 2.83 总计 Total 3 425 100.00
下载: 导出CSV
表 3 针对105份千年桐种质开发的17个多态性SSR标记的信息
Table 3 Information of 17 polymorphic SSR markers developed for 105 V. montana germplasms
引物编号
Primer No.重复基序
Repeat motif引物序列(5′—3′)
Primer sequence (5′−3′)退火温度
Annealing temperature/℃预期大小
Expected size/bpVm-BFU01 (TTG)5 AACCACTGCTACTTCACCATTTTC 55 167 ACCCAATGTTTTCTCCGACC Vm-BFU02 (AGA)5 GTGTTGAGCCAGAAACCCATTA 55 165 CAGAGAAGCCTCGGTCCCTA Vm-BFU03 (CAT)6 TCTACTCCCACACTTCCAAAACA 55 161 TAATCTCCTTCTGTCCTTCACGA Vm-BFU04 (CAG)9 GGGAGTGGTGGCAATGGC 55 181 GCTGGGAGGCATTGTTGAAG Vm-BFU05 (AGA)9 GAGCCAAGAGAAGACGAAAAGAG 55 146 ACCGTTTACAGTGTTTCGCTATG Vm-BFU06 (GGA)8 TGTGCCGCTTGTAACTGCC 55 164 TGCGGCTGTGTCAGGTGTAG Vm-BFU07 (TCT)7 AGCCTTTGCCACTGTTGAGC 55 127 GATGGGTCCGCCAAGTTCA Vm-BFU08 (TCT)5 ATTGTGAAGGATTTGCGATGG 55 161 CGGCGAAAACGAAACAGAG Vm-BFU09 (GTT)7 TCGCCTAAGGTGGTCTTGATG 55 178 GCCCAACGAAATCTAACTCTAATAA Vm-BFU10 (TAA)12 TTCCTCCTCTGGTGACGCTT 55 252 TTCCTTCCATCATCAACTTTTACC Vm-BFU11 (TCC)7 GCCGCCGCCTACTACTTACTT 55 131 TTTCTCAAACCAAACAGGAGTTG Vm-BFU12 (CTT)5 TATTTTTCTTGGGAGTAAAGTCACC 55 198 TATGTGAAATGGAGAGTTCGGAG Vm-BFU13 (CAG)6 TTGTCAACAAGCCTTCTCACCT 55 160 GCTCCAAGTCCCATCATCATTT Vm-BFU14 (GGC)5 CCAAAACCATCAATCTCTTCGC 55 204 TGATTTCGCACAAGTCCCAAG Vm-BFU15 (ACC)10 GCGTTCCTGACCCTACCTTT 55 184 AGAGAAACAAAAAAGCCACCAG Vm-BFU16 (TTA)12 TCCCCACAGCCATAAAACAAG 55 199 TTTCCAAAACTCTCAAACCACAA Vm-BFU17 (AAT)10 ACCCCATCTATGACATCCCACT 55 195 CCCCGTTCTTGCTCTCCC 注：重复基序括号外的数字为其重复次数。Note: the number of repeats is marked outside the parentheses of repeat motifs.
下载: 导出CSV
表 4 17个千年桐SSR位点的遗传多样性参数
Table 4 Genetic diversity parameters of 17 SSR loci in V. montana
位点 Locus Na Ne Ho He I FNull PIC FST Nm HWE NPA Vm-BFU01 2.000 1.618 0.190 0.384 0.570 0.335 0.309 0.147 1.446 *** 0 Vm-BFU02 6.000 1.491 0.086 0.331 0.704 0.590 0.310 0.012 20.486 *** 1 Vm-BFU03 10.000 2.400 0.365 0.586 1.223 0.218 0.539 0.055 4.286 *** 9 Vm-BFU04 11.000 3.262 0.476 0.697 1.620 0.178 0.670 0.064 3.629 *** 3 Vm-BFU05 4.000 2.187 0.286 0.546 0.923 0.306 0.466 0.018 13.956 *** 1 Vm-BFU06 5.000 2.595 0.417 0.618 1.075 0.179 0.547 0.018 13.645 *** 1 Vm-BFU07 6.000 3.135 0.419 0.684 1.357 0.232 0.638 0.035 6.850 *** 2 Vm-BFU08 8.000 2.191 0.583 0.546 1.160 0.054 0.510 0.034 7.124 *** 2 Vm-BFU09 4.000 2.237 0.390 0.556 0.908 0.162 0.453 0.072 3.214 *** 6 Vm-BFU10 9.000 2.078 0.200 0.521 1.149 0.450 0.491 0.063 3.695 *** 12 Vm-BFU11 13.000 6.289 0.235 0.845 2.101 0.559 0.824 0.057 4.163 *** 3 Vm-BFU12 2.000 1.873 0.320 0.469 0.659 0.186 0.358 0.050 4.742 ** 0 Vm-BFU13 14.000 3.157 0.410 0.687 1.617 0.252 0.657 0.093 2.439 *** 10 Vm-BFU14 10.000 1.886 0.248 0.472 1.144 0.304 0.457 0.097 2.338 *** 15 Vm-BFU15 11.000 3.155 0.352 0.686 1.544 0.333 0.648 0.062 3.808 *** 10 Vm-BFU16 9.000 5.253 0.548 0.814 1.787 0.193 0.782 0.056 4.256 *** 2 Vm-BFU17 14.000 3.913 0.538 0.748 1.852 0.158 0.725 0.031 7.705 *** 8 总数 Total 138.000 48.728 85 平均值 Mean 8.117 2.866 0.357 0.599 1.258 0.276 0.552 0.056 4.251 5 注：Na. 观测到的等位基因数；Ne. 有效等位基因数；Ho. 观测到的杂合度；He. 预期杂合度；I. 香农信息指数；Nm. 基因流；FNull. 无效等位基因频率；PIC. 多态信息含量；FST. 遗传分化指数；HWE. 哈迪−温伯格平衡；*. P < 0.05；**. P < 0.01； ***. P < 0.001；NPA. 私有等位基因的数量。下同。Notes: Na, observed number of alleles; Ne, effective number of alleles; Ho, observed heterozygosity; He, expected heterozygosity; I, Shannon’s information index; Nm, gene flow; FNull, null allele frequency; PIC, polymorphism information content; FST, genetic differentiation index; HWE, Hardy-Weinberg equilibrium; * means P < 0.05; ** means P < 0.01; *** means P < 0.001; NPA, number of private alleles. The same below.
下载: 导出CSV
表 5 具有一至多个私有等位基因的千年桐种质
Table 5 V. montana germplasms with one or more private alleles
种质
Germplasm群体
PopulationNPA 私有等位基因位点
Locus of private allele种质
Germplasm群体
PopulationNPA 私有等位基因位点
Locus of private allele1-3 1 1 Vm-BFU11 3-3 3 2 Vm-BFU03, Vm-BFU17 1-4 1 2 Vm-BFU03, Vm-BFU07 3-4 3 3 Vm-BFU03, Vm-BFU14, Vm-BFU16 1-7 1 2 Vm-BFU03, Vm-BFU07 3-5 3 1 Vm-BFU14 1-15 1 1 Vm-BFU06 3-6 3 1 Vm-BFU13 1-18 1 1 Vm-BFU11 3-9 3 3 Vm-BFU03, Vm-BFU11, Vm-BFU13 1-23 1 1 Vm-BFU13 3-10 3 1 Vm-BFU15 1-27 1 2 Vm-BFU10, Vm-BFU17 3-12 3 1 Vm-BFU14 2-2 2 1 Vm-BFU10 3-13 3 1 Vm-BFU15 2-3 2 1 Vm-BFU09 3-14 3 2 Vm-BFU14, Vm-BFU17 2-4 2 1 Vm-BFU10 3-15 3 1 Vm-BFU13 2-5 2 3 Vm-BFU09, Vm-BFU10, Vm-BFU13 3-16 3 1 Vm-BFU14 2-6 2 2 Vm-BFU03, Vm-BFU13 3-17 3 1 Vm-BFU14 2-11 2 1 Vm-BFU10 3-22 3 1 Vm-BFU03 2-13 2 1 Vm-BFU15 3-23 3 2 Vm-BFU04, Vm-BFU14 2-14 2 2 Vm-BFU10, Vm-BFU15 3-24 3 1 Vm-BFU15 2-15 2 1 Vm-BFU09 3-25 3 3 Vm-BFU04, Vm-BFU14, Vm-BFU15 2-16 2 1 Vm-BFU09 3-26 3 1 Vm-BFU14 2-17 2 1 Vm-BFU09 3-27 3 2 Vm-BFU15, Vm-BFU17 2-19 2 1 Vm-BFU15 3-28 3 2 Vm-BFU08, Vm-BFU17 2-20 2 1 Vm-BFU10 3-29 3 1 Vm-BFU14 2-23 2 1 Vm-BFU10 3-30 3 1 Vm-BFU14 2-24 2 1 Vm-BFU09 3-31 3 1 Vm-BFU14 2-26 2 1 Vm-BFU02 3-32 3 2 Vm-BFU03, Vm-BFU05 2-27 2 2 Vm-BFU10, Vm-BFU15 3-33 3 1 Vm-BFU16 2-28 2 1 Vm-BFU13 3-35 3 1 Vm-BFU14 2-29 2 1 Vm-BFU13 10-7 3 1 Vm-BFU17 2-30 2 2 Vm-BFU10, Vm-BFU13 10-13 3 2 Vm-BFU03, Vm-BFU15 2-31 2 1 Vm-BFU13 10-16 3 1 Vm-BFU14 2-33 2 1 Vm-BFU10 10-18 3 1 Vm-BFU10 3-1 3 1 Vm-BFU17 Guizhou-2 3 2 Vm-BFU08, Vm-BFU17 3-2 3 1 Vm-BFU14 Guizhou-27 3 1 Vm-BFU04 总计 Total 85
下载: 导出CSV
表 6 千年桐群体间与群体内的AMOVA分析
Table 6 AMOVA analysis among and within V. montana populations
变异来源
Source of variation自由度
Degree of freedom方差和
Sum of squares变异组分
Variance component占总变异的比例
Proportion in total variation/%FST 群体间 Among populations 2 110.841 1.209 8 0.080 (P < 0.001) 群体内 Within population 102 1 423.025 13.951 92 总计 Total 104 1 533.867 15.160 100
下载: 导出CSV
-
[1] Xu W, Yang Q, Huai H Y, et al. Development of EST-SSR markers and investigation of genetic relatedness in tung tree[J]. Tree Genetics and Genomes, 2012, 8(4): 933−940. doi: 10.1007/s11295-012-0481-z
[2] Zhang Q Y, Gao M, Wu L W, et al. Expression network of transcription factors in resistant and susceptible tung trees responding to Fusarium wilt disease[J]. Industrial Crops and Products, 2018, 122: 716−725. doi: 10.1016/j.indcrop.2018.05.041
[3] Xu W, Yang Q, Huai H Y, et al. Microsatellite marker development in tung trees (Vernicia montana and V. fordii, Euphorbiaceae) [J]. American Journal of Botany, 2011, 98(8): e226−e228. doi: 10.3732/ajb.1100151
[4] Chen Y H, Chen J H, Chang C Y, et al. Biodiesel production from tung (Vernicia montana) oil and its blending properties in different fatty acid compositions[J]. Bioresource Technology, 2010, 101(24): 9521−9526. doi: 10.1016/j.biortech.2010.06.117
[5] He Z Q, Chapital D C, Cheng H N, et al. Application of tung oil to improve adhesion strength and water resistance of cottonseed meal and protein adhesives on maple veneer[J]. Industrial Crops and Products, 2014, 61: 398−402. doi: 10.1016/j.indcrop.2014.07.031
[6] Chen Y C, Yin H F, Gao M, et al. Comparative transcriptomics atlases reveals different gene expression pattern related to Fusarium wilt disease resistance and susceptibility in two Vernicia species[J/OL]. Frontiers in Plant Science, 2016, 7: 1974 [2020−12−25]. http://doi.org/10.3389/fpls.2016.01974.
[7] Zhang L L, Liu X L, Peng J H. Genetic diversity and geographic differentiation of tung tree, Vernicia fordii (Euphorbiaceae), a potential biodiesel plant species with low invasion risk[J/OL]. Agronomy, 2019, 9(7): 402 [2020−12−27]. http://doi.org/10.3390/agronomy9070402.
[8] Zhang L L, Luo M C, You F M, et al. Development of microsatellite markers in tung tree (Vernicia fordii) using cassava genomic sequences[J]. Plant Molecular Biology Reporter, 2015, 33(4): 893−904. doi: 10.1007/s11105-014-0804-3
[9] Luo X, Cao S Y, Hao Z X, et al. Analysis of genetic structure in a large sample of pomegranate (Punica granatum L.) using fluorescent SSR markers[J]. The Journal of Horticultural Science and Biotechnology, 2018, 93(6): 659−665. doi: 10.1080/14620316.2018.1432994
[10] Agarwal M, Shrivastava N, Padh H. Advances in molecular marker techniques and their applications in plant sciences[J]. Plant Cell Reports, 2008, 27(4): 617−631. doi: 10.1007/s00299-008-0507-z
[11] Francki M G, Walker E, Li D A, et al. High-density SNP mapping reveals closely linked QTL for resistance to Stagonospora nodorum blotch (SNB) in flag leaf and glume of hexaploid wheat[J]. Genome, 2018, 61(2): 145−149. doi: 10.1139/gen-2017-0203
[12] Kalia R K, Rai M K, Kalia S, et al. Microsatellite markers: an overview of the recent progress in plants[J]. Euphytica, 2011, 177(3): 309−334. doi: 10.1007/s10681-010-0286-9
[13] Lin E P, Zhuang H B, Yu J J, et al. Genome survey of Chinese fir (Cunninghamia lanceolata): identification of genomic SSRs and demonstration of their utility in genetic diversity analysis[J/OL]. Scientific Reports, 2020, 10(1): 4698 [2020−11−20]. http://doi.org/10.1038/s41598-020-61611-0.
[14] Jain N, Patil G B, Bhargava P, et al. In Silico mining of EST-SSRs in Jatropha curcas L. towards assessing genetic polymorphism and marker development for selection of high oil yielding clones[J]. American Journal of Plant Sciences, 2014, 5(11): 1521−1541. doi: 10.4236/ajps.2014.511167
[15] Chen Y N, Dai X G, Yin T M, et al. DNA fingerprinting of oil camellia cultivars with SSR markers[J/OL]. Tree Genetics and Genomes, 2016, 12(1): 7 [2021−01−02]. http://doi.org/10.1007/s11295-015-0966-7.
[16] Bernard A, Barreneche T, Lheureux F, et al. Analysis of genetic diversity and structure in a worldwide walnut (Juglans regia L.) germplasm using SSR markers[J/OL]. PLoS One, 2018, 13(11): e0208021 [2021−01−15]. https://doi.org/10.1371/journal.pone.0208021.
[17] Babu B K, Rani K L M, Sahu S, et al. Development and validation of whole genome-wide and genic microsatellite markers in oil palm (Elaeis guineensis Jacq.): first microsatellite database (OpSatdb)[J/OL]. Scientific Reports, 2019, 9: 1899 [2020−12−29]. http://doi.org/10.1038/s41598-018-37737-7.
[18] Li D Y, Long C, Pang X M, et al. The newly developed genomic-SSR markers uncover the genetic characteristics and relationships of olive accessions[J/OL]. Peer J, 2020, 8: e8573 [2020−12−18]. http://doi.org/10.7717/peerj.8573.
[19] Baldoni L, Cultrera N G, Mariotti R, et al. A consensus list of microsatellite markers for olive genotyping[J]. Molecular Breeding, 2009, 24(3): 213−231. doi: 10.1007/s11032-009-9285-8
[20] Beghè D, Molano J F G, Fabbri A, et al. Olive biodiversity in Colombia: a molecular study of local germplasm[J]. Scientia Horticulturae, 2015, 189: 122−131. doi: 10.1016/j.scienta.2015.04.003
[21] 广西林业科学研究所, 广西崇左县油桐试验站. “桂皱-27号”等四个千年桐高产优良无性系的选育[J]. 中国林业科学, 1977(1):41−45. Guangxi Forestry Research Institute, Guangxi Chongzuo Experimental Station of Tung tree. Breeding of high-yield and excellent clones of four V. montana including “Guizhou-27”[J]. Chinese Forestry Science, 1977(1): 41−45.
[22] Nei M, Tajima F, Tateno Y. Accuracy of estimated phylogenetic trees from molecular data(II): gene frequency data[J]. Journal of Molecular Evolution, 1983, 19(2): 153−170. doi: 10.1007/BF02300753
[23] Yu C, Zhang S G, Zhou C, et al. A likelihood ratio test of population Hardy-Weinberg equilibrium for case-control studies[J]. Genetic Epidemiology, 2009, 33(3): 275−280. doi: 10.1002/gepi.20381
[24] Barchi L, Lanteri S, Portis E, et al. Identification of SNP and SSR markers in eggplant using RAD tag sequencing[J/OL]. BMC Genomics, 2011, 12: 304 [2021−01−01]. http://doi.org/10.1186/1471-2164-12-304.
[25] Gao Y, Yin S, Liu C, et al. A rapid approach for SSR development in Amorphophallus paeoniifolius using RAD-seq[J]. Taiwania, 2018, 63(3): 281−285.
[26] Poland J A, Brown P J, Sorrells M E, et al. Development of high-density genetic maps for barley and wheat using a novel two-enzyme genotyping-by-sequencing approach[J/OL]. PLoS One, 2012, 7(2): e32253 [2021−01−22]. https://doi.org/10.1371/journal.pone.0032253.
[27] Feng J Y, Zhao S, Li M, et al. Genome-wide genetic diversity detection and population structure analysis in sweetpotato (Ipomoea batatas) using RAD-seq[J]. Genomics, 2020, 112(2): 1978−1987. doi: 10.1016/j.ygeno.2019.11.010
[28] Dettori M T, Micali S, Giovinazzi J, et al. Mining microsatellites in the peach genome: development of new long-core SSR markers for genetic analyses in five Prunus species[J/OL]. SpringerPlus, 2015, 4(1): 337 [2020−12−12]. http://doi.org/10.1186/s40064-015-1098-0.
[29] Song Q J, Fickus E W, Cregan P B. Characterization of trinucleotide SSR motifs in wheat[J]. Theoretical and Applied Genetics, 2002, 104(2−3): 286−293.
[30] Song Q J, Marek L F, Shoemaker R C, et al. A new integrated genetic linkage map of the soybean[J]. Theoretical and Applied Genetics, 2004, 109(1): 122−128. doi: 10.1007/s00122-004-1602-3
[31] Guo R, Landis J B, Moore M J, et al. Development and application of transcriptome-derived microsatellites in Actinidia eriantha (Actinidiaceae)[J/OL]. Frontiers in Plant Science, 2017, 8: 1383 [2021−01−12]. http://doi.org/10.3389/fpls.2017.01383.
[32] Asadi A, Ebrahimi A, Rashidi-Monfared S, et al. Comprehensive functional analysis and mapping of SSR markers in the chickpea genome (Cicer arietinum L.)[J/OL]. Computational Biology and Chemistry, 2020, 84: 107169 [2020−12−15]. https://doi.org/10.1016/j.compbiolchem.2019.107169.
[33] Delgado-Martinez F J, Amaya I, Sanchez-Sevilla J F, et al. Microsatellite marker-based identification and genetic relationships of olive cultivars from the Extremadura region of Spain[J]. Genetics and Molecular Research, 2012, 11(2): 918−932. doi: 10.4238/2012.April.10.7
[34] Nachimuthu V V, Muthurajan R, Duraialaguraja S, et al. Analysis of population structure and genetic diversity in rice germplasm using SSR markers: an initiative towards association mapping of agronomic traits in Oryza sativa[J/OL]. Rice, 2015, 8: 30 [2021−01−04]. http://doi.org/10.1186/s12284-015-0062-5.
[35] Vischi M, Chiabà C, Raranciuc S, et al. Genetic diversity of walnut (Juglans Regia L.) in the eastern Italian Alps[J/OL]. Forests, 2017, 8(3): 81 [2021−02−01]. http://doi.org/10.3390/f8030081.
[36] Torokeldiev N, Ziehe M, Gailing O, et al. Genetic diversity and structure of natural Juglans regia L. populations in the southern Kyrgyz Republic revealed by nuclear SSR and EST-SSR markers[J/OL]. Tree Genetics and Genomes, 2019, 15(1): 5 [2021−01−10]. https://doi.org/10.1007/s11295-018-1311-8.
[37] Balapanov I, Suprun I, Stepanov I, et al. Comparative analysis Crimean, Moldavian and Kuban Persian walnut collections genetic variability by SSR-markers[J]. Scientia Horticulturae, 2019, 253: 322−326. doi: 10.1016/j.scienta.2019.04.014
[38] Degirmenci F O, Acar P, Kaya Z. Consequences of habitat fragmentation on genetic diversity and structure of Salix alba L. populations in two major river systems of Turkey[J]. Tree Genetics and Genomes, 2019, 15(4): 59. doi: 10.1007/s11295-019-1365-2
[39] Hmmam I, Mariotti R, Ruperti B, et al. Venetian olive (Olea europaea) germplasm: disclosing the genetic identity of locally grown cultivars suited for typical extra virgin oil productions[J]. Genetic Resources and Crop Evolution, 2018, 65(6): 1733−1750. doi: 10.1007/s10722-018-0650-5
[40] Jones A G, Ardren W R. Methods of parentage analysis in natural populations[J]. Molecular Ecology, 2003, 12(10): 2511−2523. doi: 10.1046/j.1365-294X.2003.01928.x
[41] Noormohammadi Z, Trujillo I, Belaj A, et al. Genetic structure of Iranian olive cultivars and their relationship with Mediterranean’s cultivars revealed by SSR markers[J]. Scientia Horticulturae, 2014, 178: 175−183. doi: 10.1016/j.scienta.2014.08.002
[42] Mariotti R, Cultrera N G M, Mousavi S, et al. Development, evaluation, and validation of new EST-SSR markers in olive (Olea europaea L.)[J/OL]. Tree Genetics and Genomes, 2016, 12(6): 120 [2020−12−22]. http://doi.org/10.1007/s11295-016-1077-9.
[43] Boucheffa S, Miazzi M M, di Rienzo V, et al. The coexistence of oleaster and traditional varieties affects genetic diversity and population structure in Algerian olive (Olea europaea) germplasm[J]. Genetic Resources and Crop Evolution, 2017, 64(2): 379−390. doi: 10.1007/s10722-016-0365-4
[44] Li X, Li M, Hou L, et al. De Novo transcriptome assembly and population genetic analyses for an endangered Chinese endemic Acer miaotaiense (Aceraceae)[J/OL]. Genes, 2018, 9(8): 378 [2020−12−14]. https://doi.org/10.3390/genes9080378.
[45] Wright S. Evolution and the genetics of populations[J]. Experimental Results and Evolutionary Deductions, 1977, 59: 815−826.
[46] Ithnin M, Teh C K, Ratnam W. Genetic diversity of Elaeis oleifera (HBK) Cortes populations using cross species SSRs: implication’s for germplasm utilization and conservation[J/OL]. BMC Genetics, 2017, 18(1): 37 [2021−01−12]. http://doi.org/10.1186/s12863-017-0505-7.
[47] Deng H W, Chen W M, Recker R R. Population admixture: detection by Hardy-Weinberg test and its quantitative effects on Linkage-Disequilibrium methods for localizing genes underlying complex traits[J]. Genetics, 2001, 157(2): 885−897. doi: 10.1093/genetics/157.2.885
[48] Han H, Woeste K E, Hu Y H, et al. Genetic diversity and population structure of common walnut (Juglans regia) in China based on EST-SSRs and the nuclear gene phenylalanine ammonia-lyase (PAL)[J/OL]. Tree Genetics & Genomes, 2016, 12(6): 111 [2021−01−13]. http://doi.org/10.1007/s11295-016-1064-1.
[49] 周强, 周慧杰, 李史干, 等. 广西红水河流域生态经济区划初探[J]. 企业科技与发展, 2017(10):28−32. doi: 10.3969/j.issn.1674-0688.2017.10.010 Zhou Q, Zhou H J, Li S G, et al. A preliminary study on the ecological economic division of the Hongshui River Basin in Guangxi[J]. Technology and Development of Enterprise, 2017(10): 28−32. doi: 10.3969/j.issn.1674-0688.2017.10.010
[50] Huang Y, Wang H, Xiao W H, et al. Contributions of climate change and anthropogenic activities to runoff change in the Hongshui River, southwest China[J/OL]. IOP Conference Series:Earth and Environmental Science, 2018, 191(1): 012143 [2021−01−22]. http://doi.org/10.1088/1755-1315/191/1/012143.
[51] Lundqvist E, Andersson E. Genetic diversity in populations of plants with different breeding and dispersal strategies in a free-flowing boreal river system[J]. Hereditas, 2001, 135(1): 75−83.
[52] Liu Y F, Wang Y, Huang H W. High interpopulation genetic differentiation and unidirectional linear migration patterns in Myricaria laxiflora (Tamaricaceae) , an endemic riparian plant in the three gorges valley of the Yangtze River[J]. American Journal of Botany, 2006, 93(2): 206−215. doi: 10.3732/ajb.93.2.206
[53] Kikuchi S, Suzuki W, Sashimura N. Gene flow in an endangered willow Salix hukaoana (Salicaceae) in natural and fragmented riparian landscapes[J]. Conservation Genetics, 2011, 12(1): 79−89.
[54] Wu Y, Yang Y D, Qadri R, et al. Development of SSR markers for coconut (Cocos nucifera L.) by selectively amplified microsatellite (SAM) and its applications[J]. Tropical Plant Biology, 2019, 12(1): 32−43. doi: 10.1007/s12042-018-9215-1
[55] Emanuelli F, Lorenzi S, Grzeskowiak L, et al. Genetic diversity and population structure assessed by SSR and SNP markers in a large germplasm collection of grape[J/OL]. BMC Plant Biology, 2013, 13(1): 39 [2021−02−22]. http://doi.org/10.1186/1471-2229-13-39.
-
期刊类型引用(8)
1. 颜欢欢. 不同基质对千年桐容器苗光合生理特征的影响. 福建林业科技. 2024(02): 31-36 . 
百度学术
2. 赵永树，张世超，杨澜，谢佳鑫，刘明，杨汉波，辜云杰，彭建，刘闵豪. 基于SSR标记的四川楠木种源（家系）遗传多样性分析. 四川林业科技. 2024(03): 1-8 . 百度学术
3. 王芝懿，李振芳，彭婵，陈英，张新叶. 基于荧光SSR标记的紫薇遗传多样性分析. 南京林业大学学报(自然科学版). 2023(02): 61-69 . 百度学术
4. 赵慧琪，李栋梁，赵慧琳，吴芳芳，黄赛，张永豪，张浪，杨珺，任军方. 海南岛火焰兰SSR标记开发、遗传多样性与群体结构分析. 分子植物育种. 2023(07): 2300-2310 . 百度学术
5. 侯佳音，冯树香，代嵩华，闫淑芳. 7种观赏桃新种质的TP-M13-SSR分子标记鉴定及亲缘关系分析. 北京林业大学学报. 2023(08): 132-141 . 本站查看
6. 莫长明，谢文娟，郭文锋，张燕玲，黄江，张馨丹，李忠，唐其，马小军. 罗汉果遗传多样性与群体结构及核心种质研究. 中草药. 2023(18): 6040-6054 . 百度学术
7. 黄彬，黄建建，汤优令，周文才，汪雁楠，田仟仟，龚春，舒惠理，温强. 江西省浙江红山茶及近缘种群体遗传多样性的SSR分析. 江西农业大学学报. 2023(05): 1084-1095 . 百度学术
8. 秦英之，车佳航，尹跃，饶书培，安巍，戴国礼，陈金焕. 基于全长转录组信息的枸杞SSR标记开发. 植物遗传资源学报. 2022(06): 1816-1827 . 百度学术
其他类型引用(8)
计量
- 文章访问数: 1106
- HTML全文浏览量: 377
- PDF下载量: 81
- 被引次数: 16
