欢迎访问《草业学报》官方网站,今天是 分享到:

草业学报 ›› 2024, Vol. 33 ›› Issue (11): 84-105.DOI: 10.11686/cyxb2023492

• 研究论文 • 上一篇    下一篇

青藏高原蕨麻的分子谱系地理学研究

白世俊1,2,3(), 李军乔1,2,3(), 刘欣1,2,3, 吕博文1,2,3   

  1. 1.青海民族大学生态环境与资源学院,青海 西宁 810007
    2.青海省特色经济植物高值化利用重点实验室,青海 西宁 810007
    3.青海民族大学青藏高原蕨麻产业研究院,青海 西宁 810007
  • 收稿日期:2023-12-19 修回日期:2024-03-21 出版日期:2024-11-20 发布日期:2024-09-09
  • 通讯作者: 李军乔
  • 作者简介:E-mail: jqlily2002@126.com
    白世俊(1970-),男,青海西宁人,副教授,本科。E-mail: bsj_1001@126.com
  • 基金资助:
    青海省重点研发与转化计划-蕨麻种质资源保存及品种提纯复壮关键技术研究与示范(2022-NK-120);青海省创新平台建设专项-青藏高原种质资源研究与利用实验室(2022-ZJ-Y01)

Phylogeography of Potentilla anserina distributed across the Qinghai-Tibet Plateau

Shi-jun BAI1,2,3(), Jun-qiao LI1,2,3(), Xin LIU1,2,3, Bo-wen LYU1,2,3   

  1. 1.College of Ecological Environment and Resources,Qinghai Nationalities University,Xining 810007,China
    2.Key Laboratory of High-Value Utilization of Characteristic Economic Plants of Qinghai Province,Xining 810007,China
    3.Tibetan Plateau Juema (Potentilla anserina) Research Centre,Qinghai Nationalities University,Xining 810007,China
  • Received:2023-12-19 Revised:2024-03-21 Online:2024-11-20 Published:2024-09-09
  • Contact: Jun-qiao LI

摘要:

青藏高原的隆升以及第四纪冰期气候的循环波动,对于青藏高原及其周边地区动植物的分布和遗传结构具有很大的影响。蕨麻是青藏高原极富营养、药用和生态价值的特有植物资源,蕨麻和鹅绒委陵菜的分类关系及分布在学术界存在争议,分子谱系地理学研究将为蕨麻遗传多样性形成机制及推断该物种迁移演化历史提供依据。以蕨麻为研究对象,对采集的30个居群的810个个体进行了叶绿体trnL-trnF序列和核基因ITS测序,揭示遗传变异在居群内和居群间的分布格局,结合群体遗传学和系统发生学,分析该物种的遗传结构与历史事件之间的关联,揭示物种及物种内不同种群形成现有分布格局的历史原因和演化过程。主要结论有:1)蕨麻具有较高水平的遗传多样性。cpDNA trnL-trnF片段共检测到40种单倍型,16个为共享单倍型,占比40%,24个为居群特有单倍型,占比60%,遗传多样性h=0.7078,单倍型多样性Hd=0.8217,核苷酸多样性π=0.010641,总遗传多样性HT=0.849;nrDNA ITS片段共检测到128种单倍型,共享单倍型42种,占比32.8%,居群特有单倍型86种,占比67.2%,遗传多样性h=0.7633,单倍型多样性Hd=0.8168,核苷酸多样性π=0.003584,总遗传多样性HT=0.844。2)居群内的遗传多样性大于居群间的遗传多样性。序列分析结果为居群内和居群间的遗传多样性都很高(cpDNA trnL-trnF:HT=0.849,居群内平均遗传多样性HS=0.640;nrDNA ITS序列:HT=0.844,HS=0.763,HT均大于HS)。蕨麻居群分为3个组:青海高原组、横断山脉组和藏南谷地组。蕨麻的遗传变异主要来源于居群内部。3)蕨麻种群具有明显的谱系地理结构。cpDNA trnL-trnF序列和nrDNA ITS遗传多态性分析及地理分布模式检验,cpDNA trnL-trnF序列:遗传分化系数GST=0.246,NST=0.417,nrDNA ITS序列:GST=0.096,NST=0.522,NST均显著大于GSTP<0.001),表明蕨麻所有居群单倍型存在显著的谱系地理结构,两种方法的结果高度一致。分子变异分析(AMOVA)表明,大部分遗传变异(59.69%)存在于居群内部,居群间分化水平很高(FST=0.40313)。4)共享单倍型和特有单倍型均由古老单倍型衍生而来。cpDNA trnL-trnF和nrDNA序列中央连接网状图呈以共享单倍型M4和H9位于中心,M1、M3和H2、H10、H11、H12位于主干位置的星状结构,其余共享单倍型和特有单倍型均由这些古老单倍型衍生而来,两者结果一致。5)蕨麻种群大小和范围发生过大规模扩张。利用cpDNA trnL-trnF和nrDNA序列进行歧点分析,前者歧点分布呈单峰曲线,表明近期群体大小和范围有大规模扩张发生;后者的歧点分布呈双峰曲线,反映基因谱系的高度复杂性,但Tajima’s D,Fu and Li’s D和Fu and Li’s F均为负值,且结果显著,且离差平方和(SSD)和扩张评估指数(HRag)的统计检验不显著,表明蕨麻居群近期有扩张的可能。6)蕨麻存在3个冰期避难所,即东喜马拉雅区域、青藏高原东南边缘及横断山脉区域。

关键词: 蕨麻, 叶绿体trnL-trnF序列, 核基因ITS, 青藏高原, 分子谱系地理学

Abstract:

The uplift of the Qinghai-Tibetan Plateau and the cyclic fluctuation of the climate during the Quaternary Ice Age have greatly affected the distribution and genetic structure of plants and animals in this area and neighboring areas. Juema is an endemic plant of the Tibetan Plateau with high nutritional, medicinal, and ecological values. The taxonomic relationship and distribution of juema and Potentilla anserina are controversial in academia. Therefore, we conducted molecular genealogical and geographical analyses to explore the formation of genetic diversity in juema and to infer its migratory evolutionary history. We collected 810 individual plants from 30 populations and determined the sequences of their chloroplast trnL-trnF and nuclear internal transcribed spacer (ITS) regions. The data were analyzed to reveal the distribution of genetic variation within and among populations, and combined population genetics and phylogenetic analyses were conducted to determine the correlation between the genetic structure of the species and historical events. The aim of these analyses was to determine the evolutionary history of the species and the reasons for the current distribution pattern of different populations within the species. The main results were as follows: 1) Juema shows a high level of genetic diversity. Forty haplotypesof the cpDNA trnL-trnF fragment were detected, including 16 shared haplotypes, accounting for 40% of the total, and 24 population-specific haplotypes, accounting for 60% of the total (genetic diversity h=0.7078, haplotype diversity Hd=0.8217, nucleotide diversity π=0.010641, total genetic diversity HT=0.849). For the ITS region, 128 haplotypes were detected, including 42 shared haplotypes, accounting for 32.8% of the total, and 86 population-specific haplotypes, accounting for 67.2% of the total (h=0.7633, Hd=0.8168, π=0.003584, and HT=0.844). 2) Genetic diversity is greater within than among populations. Sequence analyses revealed high genetic diversity both within and among populations (cpDNA trnL-trnF: HT=0.849, gene diversity within populations HS=0.640; nrDNA ITS: HT=0.844, HS=0.763; HT>HS for both sequences). The juema populations were divided into three groups: The Qinghai Plateau group, the Hengduan Mountains group, and the South Tibetan Valley group. The genetic variation of juema has mainly originated from within populations. 3) Juema populations show a distinct genealogical geographic structure. We conducted genetic polymorphism and geographic distribution pattern analyses based on the cpDNA trnL-trnF and nrDNA ITS data. For cpDNA trnL-trnFthe genetic differentiation coefficients were GST=0.246 and NST=0.417; for nrDNA ITS, GST=0.096 and NST=0.522. For both sequences, NST was significantly larger than GSTP<0.001). The results were highly consistent between the two methods, and indicated that a significant spectral geographic structure exists in all population haplotypes of juema. Analysis of molecular variance (AMOVA) showed that the majority of genetic variation (59.69%) is within populations, with a high level of differentiation among populations (FST=0.40313). 4) Both shared and endemic haplotypes are derived from ancient haplotypes. Centrally linked reticulation maps based on cpDNA trnL-trnF and nrDNA sequences showed a stellate structure with the shared haplotypes M4 and H9 at the center, and the shared haplotypes M1 and M3, and H2, H10, H11, and H12 at the trunk positions. The remaining shared and unique haplotypes were derived from these ancient haplotypes, and the results were consistent. 5) There has been a major expansion in the size and range of the plant population. When cpDNA trnL-trnF and nrDNA sequences were subjected to a manifold analysis, the former showed a single-peak curve in the manifold distribution, indicating that a large-scale expansion of the size and range of the population has occurred recently; and the latter showed a double-peaked curve, reflecting a high degree of complexity in its genetic genealogy. The results of Tajima’s D, Fu and Li’s D, and Fu and Li’s F tests were all negative and significant, and the statistical tests of sum of squares deviations and Harpending’s raggedness parameters were not significant, suggesting the possibility of recent expansion of juema populations. 6) Juema exists in three Ice Age refuges, i.e., the Eastern Himalaya region, the southeastern edge of the Tibetan Plateau, and the Hengduan Mountains region.

Key words: Potentilla anserina, chloroplast trnL-trnF fragment, nuclear gene ITS, Qinghai-Tibet Plateau, phylogeography