论文部分内容阅读
八角莲属(Dysosma)植物隶属于小檗科(Berberidaceae)鬼臼亚科(Podophylloideae),其根状茎含有的生物活性成分鬼臼毒素(Podophylloxotin)等是合成抗癌药物的前体。该属植物分布于我国东南、中南至西南地区,其中,八角莲(D.versipellis)为广布种,其余种均为地方特有种。最西部的西藏八角莲(D.tsayuensis)分布于东喜马拉雅,最东部的六角莲(D.pleiantha)间断分布于我国东部与台湾。该属有着近200年的分类历史,但由于其形态性状的变异性以及采样不足,使得它与近缘属的分类问题和系统发育问题一直没有得到很好地解决。以前的系统发育和亲缘地理学研究仅涉及其中的少数种,因此有关属内7个种的物种形成与分化机制以及该属的整合分类系统是一个亟需解决的科学问题。基于本课题组对八角莲属近10年的材料积累和广泛的野外资源调查取样,在形态学研究的基础上,利用系统发育、生物地理以及群体遗传学原理与方法,对鬼臼亚科的系统发育以及八角莲属植物地理分布与物种形成的时间格局与机制进行了全面深入的研究,获得以下主要结果: 1)鬼臼亚科的系统发育及分化时间: 对鬼臼亚科4属12个种进行了广泛的野外取样,采集了中国、日本和北美的63个群体材料,采集群体涵盖了八角莲属的整个地理分布区域。以Wang et al.(2007)的小檗科主要类群的系统发育树为基础(matK、rbcL),补充八角莲属缺失的6个种,并对鬼臼亚科12个种增加了叶绿体trnL-trnF、atpH-atpI和ITS和ETS的扩增,联合6个片段,通过系统发育分析构建了鬼臼亚科12个种的分子系统进化树,并结合两个化石校正点(小檗科的crown节点分化时间91 Ma,Ranzania、Berberis和Mahonia三个属的crown节点时间45Ma)和鬼臼亚科二次校正点(8.8±1.4 Ma),利用BEAST进行分化时间估算,揭示了鬼臼亚科属间及属内种间的系统发育关系:鬼臼亚科是一个单系类群(PP=1),山荷叶属为一个单系(PP=1)位于鬼臼亚科的基部,八角莲属是一个很好的单系(PP=1),与东亚-北美间断分布的姐妹类群桃儿七和足叶草组成一个支系,但支持率较低(PP=0.67)。八角莲属内的广布种八角莲不是单系,而是与六角莲形成并系,分子系统树支持其它形态学种为单系类群。分歧时间估算结果表明鬼臼亚科的共祖时间为上新世[c.12.67 Ma(10.71-14.56 Ma)];山荷叶属的共祖时间为7.96(5.19-10.61) Ma;八角莲属、桃儿七属和足叶草属的分化(7.96-9.14 Ma)与中新世的全球气候变冷有关。BBM结果揭示了八角莲属起源于东喜马拉雅,并自上新世由西向东迁移。在八角莲属中,位于西部的云南八角莲、西藏八角莲和川八角莲分化时间较早(7.76-9.14 Ma),其分化时间与青藏高原早期快速隆升时间(8-13 Ma)一致,由此可知,早期地理隔离导致了西部类群异域成种。 2)八角莲属种水平亲缘地理学 选取八角莲属7个种46个群体进行了亲缘地理学研究,研究结果揭示了八角莲属内7个种的谱系关系及遗传多样性的时空分布格局。3个叶绿体片段(trnL-trnF、trnL-ndhJ、trnS-trnfM)序列变异的遗传多样性分析揭示了八角莲属共有66个叶绿体单倍型,八角莲属水平的遗传多样性相对较高(hT=0.929,πT=3.58×10-3),并具有显著的谱系地理结构(NST=0.866,GST=0.790,P<0.05)。广布种八角莲的单倍型多样性和核苷酸多样性(hT=0.314;πT=0.25×10-3)比其它窄域种都高;八角莲(NST=0.818,GST=0.611,P<0.05)、小八角莲(NST=0.865,GST=0.728,P<0.01)和川八角莲(NST=0.977,GST=0.870,P<0.05)的谱系地理结构明显。TCS单倍型谱系关系和系统发育分析结果均揭示,八角莲属共分为6个与地理分布对应的谱系(Lineage1-Lineage6),它包括分布于中国西部的云南八角莲(Lineage6)、西藏八角莲(Lineage5)和川八角莲(Lineage4)依次分化出来。分布于中国中部-东南部的小八角莲、贵州八角莲、六角莲及八角莲聚成一个单系[Lineage1+(Lineage2+Lineage3)];其中,广布种八角莲单倍型分为两个谱系(西北谱系Lineage2和东南谱系Lineage1),是一个并系类群;单系类群六角莲(Lineage3)与八角莲西北部谱系Lineage2为姐妹关系;小八角莲和贵州八角莲没有独立分支,小八角莲(H1,H52-56)镶嵌于八角莲的东南谱系中,并共享单倍型H1;贵州八角莲的所有个体和小八角莲共享一个单倍型H55。基于ITS序列分析表明,八角莲属内7个种没有共有单倍型。cpDNA与ITS单倍型的谱系关系在小八角莲、贵州八角莲、八角莲以及六角莲之间呈现不一致的现象,表明了这3个种在同域分布区可能存在杂交与基因渐渗现象。 失配分析表明,八角莲东南谱系(Lineage1)和六角莲谱系(Lineage3)都经历了显著的群体扩张,但八角莲西北谱系(Lineage2)和川八角莲谱系(Lineage4)长期居留原地。八角莲东南谱系的扩张时间为更新世早期[1.78 Ma(95% CI:1.01-3.03 Ma)],表明气候变冷导致暖温带落叶林下的八角莲表现出扩张信号。六角莲的扩张时间为更新世晚期[0.18 Ma(95% CI:0.00072-0.75 Ma)],与倒数第二个间冰期时间一致一致(c.0.25-0.135 Ma),表明气候变暖使暖温带常绿林的六角莲呈现扩张信号。 3) DNA条形码(DNA barcoding) 通过对9个候选的DNA条形码(ITS、matK、rbcL、trnL-trnF、trnL-ndhJ、trnS-trnfM、atpH-atpI、rpl32-trnL、rps18-clpp)在鬼臼亚科植物中的鉴定能力进行了评估,比较分析了各条形码的扩增成功率、种内和种间变异、barcodinggap,并运用BLAST、DISTANCE和构建系统发育树方法评估了不同序列的物种鉴定成功率,结果表明除rps18-clpp外其余8个条形码对鬼臼亚科12个种的通用性均为100%,其中通用引物ITS扩增的序列种内与种间变异的差异、barcodinggap较其他序列具有更明显的优势,其种鉴定成功率相对较高(83.3%)。因此,ITS是适合鬼臼亚科植物鉴别最好的DNA条形码。 4)云南八角莲和八角莲的比较转录组 通过Hiseq2000测序平台,对高海拔的云南八角莲和低海拔的八角莲进行高通量转录组测序,运用de novo从头组装的方法,分别获得了53929和44855条unigenes;基于序列相似性,我们找出两个种的直系同源基因(Orthologs)为4593对,其中3126对直系同源基因有Ka/Ks值,Ka/Ks>1的同源基因有202对,0.5< Ka/Ks<1的同源基因有535对,这737对同源基因只有一对具有显著P值(glycosyltransferase gene,Ka/Ks=0.606,P=0.02)。1273对表现出显著的净化选择(Ka/Ks<0.5,P<0.05)。另外,4593直系同源基因共含有SSR461个,其中,去除长度小于16bp的、侧翼序列不足够长的同源基因,最终成功设计了51对EST-SSR引物。我们检验了这些引物的通用性,其中26对(51%)和41对(80.4%)分别在鬼臼亚科和八角莲属具有通用性。最后,我们利用PCR直接扩增的方法得到了鬼臼亚科12个种33样本的19个CYP719A(鬼臼毒素合成路径的关键基因)的基因序列,并分析了它们的分子进化式样,结果阐明了CYP719A基因在鬼臼亚科经历了强烈的净化选择,在八角莲属内表现为分布海拔较高的物种具有松散的净化选择(relaxed purifying selection)或较强的正向选择(more positive selection)。比较转录组研究为后续八角莲属物种及其近缘物种的群体遗传学与进化生物学研究提供了大量的分子标记。 综上所述,鬼臼亚科是一个单系,山荷叶属为基部类群,北美山荷叶与日本山荷叶成姐妹关系,八角莲属是一个很好的单系。青藏高原快速隆升导致的地理隔离驱使了八角莲属西部类群的异域物种形成;上新世-更新世的气候动荡变化是八角莲属中东部物种分化和遗传多样性形成的主要原因,六角莲是近期分化的物种;ITS是鬼臼亚科鉴定种最好的DNA条形码。大量的属水平的EST-SSR引物和单拷贝核基因被开发出来,填补了八角莲属遗传资源缺失的空白,CYP719A在八角莲属内的分子进化模式表明该基因在高海拔的物种中经历了更强的正向选择或松散的净化选择。