药食同源植物薤白的谱系地理学研究

对我国14个地区的药食同源植物薤白（Allium macrostemon）的叶绿体基因片段（psbA-trnH、rps16和trnL-F）与核基因片段（ITS）进行测序分析,揭示薤白的遗传变异分布式样、单倍型地理分布格局,并推断其在第四纪冰期的避难所。结果表明：薤白叶绿体基因（cpDNA）遗传多样性低于核基因（nrDNA）遗传多样性（cpDNA：H_T=0.868;nrDNA：H_T=0.890）。cpDNA和nrDNA的分子变异分析（AMOVA）结果显示：薤白遗传变异主要发生在居群间（cpDNA：92.84%;nrDNA：98.40%）,存在遗传分化（cpDNA：N_st=0.918,G_st=0.866,F_st=0.928;nrDNA：N_st=0.984,G_st=0.855,F_st=0.984）,且N_st均大于G_st,表明该物种具有明显的谱系地理结构。在薤白居群中,共检测到11个叶绿体单倍型和14个nrDNA基因型;单倍型网络图及地理分布图表明,叶绿体单倍型H3、核DNA基因型H1频率最高,位于网络结构图的中心位置,可能为古老单倍型。此外,冰期避难所假说认为遗传多样性高、拥有古老单倍型和较多特有单倍型的区域可能是该物种的冰期避难所,因此推测薤白在第四纪冰期时可能在大盘山、天水和通化地区存在多个冰期避难所。这些分析可为类似草本植物的进化提供参考,丰富对东亚草本植物分子系统与生物地理学的认识。     

青藏高原岷县龙胆的遗传分化与种群动态历史分析

青藏高原是全球生物多样性中心之一,是研究物种形成和适应性进化的热点地区。岷县龙胆（Gentiana algida var. purdomii）是青藏高原地区特有植物,具有很高的药用价值。本研究以岷县龙胆种群为对象,基于不同遗传方式的分子标记,通过种群遗传结构、遗传分化、种群动态历史和物种分布模型等分析,探讨其在青藏高原的分化历史。结果表明：岷县龙胆具有高的遗传多样性,叶绿体数据表明种内遗传分化程度很高（F_ST=0.452）,但核基因数据表明种内遗传分化程度很低（F_ST=0.022）,揭示岷县龙胆种内存在强烈的基因流。种群动态历史分析表明岷县龙胆种群大小经历了近期扩张,物种分布模型的结果表明岷县龙胆在末次盛冰期以来分布范围略有扩张。分子钟模型的结果表明岷县龙胆种内的遗传分化主要发生在第四纪冰期以来。以上结果一致表明,第四纪气候波动是岷县龙胆分化的主要外部因素。本研究为龙胆属及青藏高原植物类群的物种分化、适应性研究以及资源开发与保护提供了参考。     

伏牛山区连翘遗传多样性研究

利用2个叶绿体及1个核基因片段,对伏牛山区3个连翘野生种群30个个体进行了遗传多样性分析.结果表明,叶绿体基因数据得到4个单倍型(H1～H4),单倍型多样性指数(h)为0.545±0.066,核酸多样性指数π为0.001 72±0.000 21.核基因数据得到10个单倍型(A～J),单倍型多样性指数为0.501±0.076,核酸多样性指数π为0.002 53±0.000 49.核基因种群间遗传分化(ФST)为0.167,小于叶绿体基因ФST(0.886),种群间基因流(Nm)为1.250,显著大于叶绿体基因Nm(0.030).研究发现,伏牛山区野生连翘种群目前仍具有较高的遗传多样性,但部分人类活动较多的地区已经造成了连翘遗传多样性的降低.连翘种群间地理距离会限制种子介导的基因流,但在小区域范围内对花粉所介导的基因流影响不明显.     

青藏高原东北部黄缨菊的谱系地理学

黄缨菊（Xanthopappus subacaulis）是青藏高原地区一种特有的高山草甸药用植物。为探讨第四纪冰期气候波动对黄缨菊居群遗传结构和空间分布格局的影响,对黄缨菊20个居群、123个个体的叶绿体DNA片段（psbA-trnH、rbcL和psbI-psbK）进行测序和数据分析。结果表明：黄缨菊居群共检测到6个单倍型,其中H1为古老单倍型,除居群P7外其余居群均具有单倍型H1,H3、H5和H6为特有单倍型,单倍型H3为居群P7的私有单倍型,单倍型H5和H6只存在于居群P18,单倍型H2和H4主要存在于青海湖流域的居群;总的遗传多样性（H_e）和核苷酸多样性（π）分别为3.101和0.008 903;居群间遗传变异（68.98%）大于居群内遗传变异（31.02%）,居群间遗传分化较高（F_ST=0.689 85,P<0.01）;居群遗传分化系数N_ST（0.727）大于G_ST（0.656）（P>0.05）,表明黄缨菊在分布区域内不存在明显的谱系地理结构;错配分布和中性检验结果显示,黄缨菊居群可能经历过近期扩张。据此,推测第四纪冰期黄缨菊可能在青海湖流域和甘肃临潭地区存在微型避难所,认为第四纪气候变迁及青藏高原隆升塑造了黄缨菊的现代地理分布格局。     

基于核基因和叶绿体基因的台湾水韭遗传多样性和遗传结构的分析

利用核基因和叶绿体基因对台湾水韭（Isoetes taiwanensis Devol）2个居群20个样本的遗传多样性和遗传结构进行了分析,并对台湾水韭的保育提出了建议。核基因的检测结果显示,2个居群共存在18个单倍型,单倍型多样性为0.9842,核苷酸多样性为0.00215,居群间基因流N_m为4.26,居群间分化系数G_st值为0.05543;叶绿体DNA的检测结果显示,单倍型数目为6,单倍型多样性为0.6211,核苷酸多样性为0.00326。其中台北居群的单倍型多样性和核苷酸多样性均比金门居群高。核基因和叶绿体基因的AMOVA分析显示,居群内的遗传差异远高于居群间。台湾水韭的核苷酸多样性相对其他水韭属植物较低,可能与台湾水韭的染色体倍型以及狭域分布有关。居群结构的形成可能和基因流以及奠基者效应有关。宜采用原地和迁地保护的方法对居群进行保护。     

基于叶绿体基因序列rbcL-accD研究莛子藨的分子遗传多样性和谱系地理结构

莛子藨是莛子藨属的一种多年生草本植物,主要分布在青藏高原东北部及邻近高山地区和日本,为东亚特有种。本文基于叶绿体DNA基因间隔区rbcL-accD序列,对莛子藨物种水平的序列特征、分子遗传多样性和谱系地理结构进行研究。分析结果发现该区域莛子藨种群的遗传变异主要来自于居群间（72.28%）,优势单倍型（H2）广泛分布,特有单倍型分散在多个居群。居群间遗传分化系数N_STST,没有检测到显著谱系地理结构,而且居群间基因流较低（N_m=0.096）,青藏高原及其邻近地区特殊地形和生态环境造成的地理隔离或者生境异质化可能是现有单倍型分布格局形成的原因。种群动态历史分析显示莛子藨经历过不显著的种群选择或者扩张事件,综合分析推测莛子藨居群应该是在间冰期或者冰期后从不同避难所扩散形成现在的分布格局,具体避难所还需进一步深入研究才能确认。     

全基因组重测序揭示青藏高原东缘及东南缘东方蜜蜂遗传多样性与适应性进化

【目的】对青藏高原东缘和东南缘及邻近地区东方蜜蜂Apis cerana样本进行群体遗传多样性与适应性进化研究,为进一步解析青藏高原东方蜜蜂的遗传资源多样性、种群扩散规律以及适应高原生境的分子进化机制提供参考。【方法】对青藏高原东缘和东南缘及邻近地区77群东方蜜蜂样本进行全基因组重测序;运用群体遗传学方法,基于群体结构、主成分分析、系统进化树、遗传分化指数、线粒体基因组单倍型以及选择信号分析对东方蜜蜂这77群及GenBank数据库下载的90群的全基因组重测序原始数据进行分析。【结果】这167群东方蜜蜂区分出来自川西高原、藏南高原、滇北高原和川北高原的4个高原型群,分属于两个进化支,平均遗传分化指数(Fst=0.1178)高于非高原型群的(Fst=0.0411)。基于种群间最小遗传距离分析,东南缘的藏南高原群、滇北高原群与滇南群具有更近的遗传关系;东缘的川西高原群、川北高原群分别与川西山地和秦巴群具有更近的遗传关系。结合线粒体基因组单倍型分析,初步推断出青藏高原东缘及东南缘高原型群体祖先单倍型及来源。选择性分析鉴定到了潜在的参与脂肪酸代谢、光转导、温度适应、卵巢发育等信号通路相关的潜在受选基因。在藏南高原和滇北高原群体中发现两个共同受选择基因ISL-1和FOXO,主要参与胰岛素分泌以应对细胞压力,暗示它们在东方蜜蜂适应青藏高原东南缘生境中发挥着重要作用。【结论】青藏高原东方蜜蜂遗传多样性丰富,东南缘的藏南高原和滇北高原群体以及东缘的川西高原和川北高原群体在遗传分化上明显区分,这4个高原群是由邻近地区非高原群扩散后的地理隔离形成了种群分化;初步筛选获得东方蜜蜂高原环境适应性潜在基因。本研究为进一步探析青藏高原东方蜜蜂适应高原生境下的分子进化机制奠定了基础。     

高山植物条纹狭蕊龙胆的分子亲缘地理学研究

草本植物由于较短的生活周期以及对环境变化的敏感性,可能会更好地揭示第四纪冰期以来植物居群变化的历史过程.分子亲缘地理学研究是揭示动植物居群历史的有力工具,但到目前为止对青藏高原草本植物的分子亲缘地理学研究几乎是空白.因此,本文选择在青藏高原及邻近地区生长的一年生高山草本植物条纹狭蕊龙胆为研究对象,进行了13个居群155个个体的叶绿体基因组(cpDNA)非编码片段trnH(GuG)-psbA基因间区序列变异检测,共发现7种单倍型,其中单倍型H印A是分布最广的,而单倍型Hap E、Hap F和Hap G是拥有的居群所特有的.青藏高原东北部、东部及邻近地区的每个居群拥有的单倍型非常单一,而高原东南部横断山区的单倍型分布很集中,遗传多样性也相对较高.分子变异分析(AMOVA)结果表明整个分布区条纹狭蕊龙胆的遗传变异主要存在于居群间(73.05%),且居群问的遗传分化很高(GsT=0.805,FsT=0.731,NsT=0.859),有着显著的亲缘地理学结构(NsT>Gsr,P<0.05)和较低的居群间的平均基因流(Nm=0.184).结合巢式支系法分析(NCA),根据本文的研究结果推测青藏高原东南部横断山区是该植物第四纪冰期时可能的避难所,而且在间冰期或冰期后.伴随着异域片段化和过去片段化从避难所发生范围扩张而形成当前单倍型及居群的分布格局.     

横断山区高山绣线菊的谱系地理学研究

横断山区作为青藏高原东南部主要的一个冰期避难所,第四纪冰期气候的变化对该地区的植物地理分布和居群遗传结构都产生了重要的影响。为了揭示该地区物种分布的分子系统地理学结构,选取在该地区广泛分布的一种高山灌木-高山绣线菊的叶绿体trnL-trnF序列进行研究。采集了15个居群182个个体进行测序,共发现7个单倍型。总的遗传多样性较高(HT=0.809),但居群内遗传多样性较低(HS=0.236)。分子变异分析(AMOVA)结果表明分布区内高山绣线菊的遗传变异主要存在于居群间(84.48%),且居群间的遗传分化很高(GST=0.708,FST=0.84476,NST=0.863),有着显著的谱系地理学结构(NST>GST,P<0.01)和较低的居群间平均基因流(Nm=0.09)。单倍型的系统进化树和进化分支网络分析得到了相似的拓扑结构,7种单倍型都按照地理分布聚为三支:横断山区西部、横断山区东部以及两者的交接地带。本研究推测该物种在横断山区存在多个冰期避难所,而没有表现出大规模的种群集体扩张和迁移的现象。青藏高原隆升、第四纪气候的反复波动以及横断山区特殊的地理环境使得原来连续的居群片段化,并发生范围扩张,从而塑造了高山绣线菊的现代生物地理分布格局。     

江孜沙棘和云南沙棘之间谱系分化和亲缘地理

研究具有邻近分布的近缘类群的遗传分化和谱系筛选对于进一步揭示物种形成机制具有重要意义。该文以沙棘属(Hippophae)临域分布的两个类群, 江孜沙棘(H. gyantsensis)和云南沙棘(H. rhamnoides subsp. yunnanensis)为研究对象, 进行群体水平上的母系分化研究。前一个种分布于西藏的中西部, 而后者主要分布于青藏高原东南部的云南西北部、四川西部和西藏东南部。两个类群在叶和果实性状存在明显区别。叶绿体在沙棘属植物中为母系遗传。共研究了两个类群14种群109个沙棘个体的叶绿体trnL-F、trnS-G序列; 序列排序后共发现11种单倍型, 江孜沙棘和云南沙棘分别有7种和6种, 两种单倍型为两个类群共享。分支分析和嵌套进化分析进一步表明, 两个类群之间的单倍型相互交错, 单倍型分化与形态上划分的两个类群不一致, 表明它们之间具有十分复杂的谱系筛选过程。这些发现明显不支持以前提出的有关江孜沙棘系统位置的假设。但是, 目前所获得的证据不能区分该物种究竟是通过异域分化还是同倍性杂交起源的。不同种群固定特有单倍型表明, 两个类群都在最后一次冰期可能存在多个避难所。     

Mitochondrial and chloroplast phylogeography of Picea crassifolia Kom. (Pinaceae) in the Qinghai-Tibetan Plateau and adjacent highlands

The disjunct distribution of forests in the Qinghai-Tibetan Plateau (QTP) and adjacent Helan Shan and Daqing Shan highlands provides an excellent model to examine vegetation shifts, glacial refugia and gene flow of key species in this complex landscape region in response to past climatic oscillations and human disturbance. In this study, we examined maternally inherited mitochondrial DNA (nad1 intron b/c and nad5 intron 1) and paternally inherited chloroplast DNA (trnC-trnD) sequence variation within a dominant forest species, Picea crassifolia Kom. We recovered nine mitotypes and two chlorotypes in a survey of 442 individuals from 32 populations sampled throughout the species' range. Significant mitochondrial DNA population subdivision was detected (G(ST) = 0.512; N(ST) = 0.679), suggesting low levels of recurrent gene flow through seeds among populations and significant phylogeographical structure (N(ST) > GST, P < 0.05). Plateau haplotypes differed in sequence from those in the adjacent highlands, suggesting a long period of allopatric fragmentation between the species in the two regions and the presence of independent refugia in each region during Quaternary glaciations. On the QTP platform, all but one of the disjunct populations surveyed were fixed for the same mitotype, while most populations at the plateau edge contained more than one haplotype with the mitotype that was fixed in plateau platform populations always present at high frequency. This distribution pattern suggests that present-day disjunct populations on the QTP platform experienced a common recolonization history. The same phylogeographical pattern, however, was not detected for paternally inherited chloroplast DNA haplotypes. Two chlorotypes were distributed throughout the range of the species with little geographical population differentiation (G(ST) = N(ST) = 0.093). This provides evidence for highly efficient pollen-mediated gene flow among isolated forest patches, both within and between the QTP and adjacent highland populations. A lack of isolation to pollen-mediated gene flow between forests on the QTP and adjacent highlands is surprising given that the Tengger Desert has been a geographical barrier between these two regions for approximately the last 1.8 million years.     

Responses of plants to changes in Qinghai–Tibetan Plateau and glaciations: Evidence from phylogeography of a Sibiraea (Rosaceae) complex

The genetic structure and phylogeographical history of the alpine shrubs Sibiraea angustata (Rosaceae) and Sibiraea laevigata from the Qinghai–Tibetan Plateau (QTP) were investigated to identify alpine plant responses to changes in the QTP and glaciations. Fifty-five populations were analyzed using four chloroplast DNA (cpDNA) regions and (nuclear ribosomal internal transcribed spacer) nrITS sequence data. In all, 21 cpDNA haplotypes and 13 nrITS sequence types were detected. Analyses of the genetic diversity and phylogenetic relationships detected two rarely reported glacial refugia. One was the Yushu–Nangqian area, and the other consisted of the area from the Songpan Plateau to the southeastern margin of the QTP. Sibiraea species populations experienced divergent evolution and founder effects when they recolonized the QTP platform and adjacent high-altitude regions following glaciations. The divergence times of the main lineages and haplotypes were in the range of 1.60–2.58 Ma. The population size of Sibiraea species in the QTP decreased approximately 23-fold during the last 0.12 Ma, indicating that Sibiraea species were significantly affected by environmental changes in the QTP. Therefore, the rapid uplift of the QTP and subsequent glaciations likely played an important role in driving genetic divergence and population size changes of Sibiraea species in the QTP.     

Migration patterns of Gentiana crassicaulis,an alpine gentian endemic to the Himalaya–Hengduan Mountains

The Himalaya–Hengduan Mountain region is one of the hotspots of biodiversity research. The uplift of the Qinghai–Tibetan Plateau (QTP) and the Quaternary glaciation caused great environmental changes in this region, and the responses of many species in the QTP to the Quaternary climate are still largely unknown. The genetic structure and phylogeographical history of Gentiana crassicaulis Duthie ex Burk, an endemic Chinese alpine species in this area, were investigated based on four chloroplast fragments and internal transcribed spacer region of the nuclear ribosomal DNA (nrITS) sequences of 11 populations. The populations with highly diverse chloroplast haplotypes were mainly found at the edge of the QTP. There were two main haplotypes of nrITS clones, one shared by the Yunnan and Guizhou populations, and the other by the remaining populations. The population with the highest diversity was the Gansu population, located at the edge of the plateau. Based on molecular dating, the diversification of G. crassicaulis at the edge of the plateau occurred before the Last Glacial Maximum (LGM), and the species may have completed its expansion from the edge to the platform. Ecological niche models were conducted to predict the distributional ranges of G. crassicaulis at present, during the LGM, and during the last interglacial (LIG) period. The results demonstrated that G. crassicaulis survived on the QTP platform and at the edge during the LGM but afterward retreated from the platform to the southern edge, followed by expansion to the platform.     

基于叶绿体DNA trnT-trnF序列研究祁连圆柏的谱系地理学

由于青藏高原的地貌效应,第四纪冰期气候的反复变化应对现今该地区生物的地理分布及其居群遗传结构产生重大影响。本文对这一地区特有分布物种祁连圆柏Juniperus przewalskii Kom.整个分布区内20居群392个个体的trnT-trnF序列变化进行了研究;共发现3种单倍型(haplotype),构成两种地理区域:高原台面上的居群主要固定Hap A,而Hap A、Hap B和Hap C在高原边缘居群均有分布。所有居群总的遗传多样性HT = 0.511,GST= 0.847。在低海拔的高原边缘,Hap A、Hap B和Hap C高频率固定在不同的居群中,表明可能存在多个不同的避难所,居群反复缩小和扩张的瓶颈效应造成了遗传多样性的丢失。而边缘的一个居群含有两种相似单倍型频率则可能是冰期后迁移融合而成或者该居群在冰期经受的瓶颈作用更弱。高原台面东部间断分布的居群只固定Hap A,表明它们可能经历了冰期后共同的回迁过程和由此产生的奠基者效应。我们的研究结果表明祁连圆柏在冰期可能存在多个避难所,瓶颈效应和奠基者效应造成了这些居群现在的遗传多样性分布式样。     

Multiple Pleistocene refugia and recent diversification for Streptocarpus ionanthus (Gesneriaceae) complex: Insights from multiple molecular sources

Streptocarpus ionanthus (Gesneriaceae) is endemic to Tanzania and Kenya, distributed in Tanga, Morogoro, and Kilifi regions. The species houses nine subspecies characterized by complex morphotypes and poorly understood evolutionary relationships, and thus is an ideal model for investigating evolutionary dynamics over time. Using multiple methods, we sought to test our hypothesis that the infraspecific taxa in Str. ionanthus are slightly variable and evolving populations. We first examined the genetic diversity, population differentiation, and phylogeographic structure among the populations of Str. ionanthus using both chloroplast and nuclear markers. We then estimated the divergence time of Str. ionanthus lineages and modeled past and future distribution. Despite Str. ionanthus exhibiting bottleneck events across its range, the populations maintain relatively high genetic diversity attributed to historical population admixture or local adaptation arising from habitat heterogeneity. The phylogeographic and genetic structure revealed a high connection among the Usambara mountains populations, while molecular dating suggested most diversification of haplotypes began ~1.32–0.18 million years ago and intensified toward the present, a conclusion of recent diversification. Phylogenetic relationship of Str. ionanthus cpDNA haplotypes revealed five main lineages with unique haplotypes that could be suggestive of past isolated refugia during the Pleistocene climate shifts. According to niche modeling, the stability of suitable areas during the Last Glacial Maximum (LGM) offered protective micro-habitats that have preserved the genetic diversity of Str. ionanthus to date. In conclusion, our findings suggest a complex Str. ionanthus with slightly variable lineages or populations attributed to multiple refugia and on the verge of divergence.     

Structure of Pinus sylvestris L. populations in Bulgaria revealed by chloroplast microsatellites and terpenes analysis: Provenance tests

In the present study we investigated the genetic structure and genetic diversity of Pinus sylvestris populations in Bulgaria using chloroplast microsatellite markers and terpene analysis. We were interested in addressing the following questions: (1) can population structure in Scots pine be detected via chloroplast microsatellites markers and terpenes; (2) are there differences in population differentiation between the two analyses; and (3) how are the patterns related to geographic distances. Twelve provenances were chosen throughout the species' range in Bulgaria. Following DNA extraction, chloroplast microsatellite (cpSSR) loci were surveyed using six primer pairs. Between 4 to 8 size variants were identified at each locus. A total of 35 size variants at the six loci were identified, 11 occurring at low frequencies (<1%). They were combined in 134 different haplotypes, of which seven represent 1/3 of the genetic structure. AMOVA analysis revealed that 10.99% of the variation was found among populations, while 89.01% was expressed within populations. The cpSSR analysis divided Scots pine populations into two groups, the first represented by populations located in the south-western part of the Rhodopes and Pirin mountains, while the second group is located in the northeast of Rhodopes and Rila mountains. Terpene analysis revealed that on average, 53% of the monoterpene pool in P. sylvestris was accounted for by -pinene (range 47–59%) followed by β-pinene (range 6–12%). The presence of two distinct groups is weekly consistent with physical distances between populations, similar significant correlation between genetic distance determined by chloroplast microsatellites analysis and chemotype distance (determined by terpenes) was observed. Our results suggest that the structural pattern of genetic diversity of cpDNA in Scots pine populations is the consequence of historical biogeographic processes.     

Phylogeography of Japanese water crowfoot based on chloroplast DNA haplotypes

Water crowfoot, Ranunculus subgenus Batrachium, are submerged macrophytes in the Ranunculaceae. We aimed to infer the Quaternary history of these aquatic macrophytes in the Japanese archipelago. We studied 212 individuals of three perennial and one annual species from 46 populations covering the entire geographic range and found eight haplotypes based on approximately 1800 bp of four spacers in chloroplast DNA. The relationships among haplotypes were resolved using maximum parsimony and parsimony network analyses. To identify the zones of clear genetic boundaries, Monmonier's algorithm was used. The eight haplotypes were distinguished from adjacent haplotypes by one substitution or indel. Each of the 46 populations was fixed for a single haplotype. Inconsistency between cpDNA haplotypes and Batrachium taxa was found, except in one annual species. The distribution of the haplotypes in perennial species was highly geographically structured. An abrupt genetic change was detected between the Tohoku region and more southerly regions. Since the perennial Batrachium are cold-adapted, this genetic differentiation may be due to historical changes in their distributions caused by Quaternary climatic oscillations: interglacial retreats in colder refugia and glacial range expansions. The single haplotype composition of each population may have been shaped by founder effects during colonization and/or by genetic drift. The interglacial refugial populations must have been small enough to deplete haplotype diversity. Inconsistency between cpDNA haplotypes and Batrachium taxa may be due to incomplete lineage sorting of ancestral polymorphic haplotypes.