薜荔之传粉小蜂 ——我国东南沿海薜荔小蜂隐存种及分布格局

论文摘要

共生是最具影响力的种间相互作用之一，对生命各个层次都有深刻的且常常是决定性的影响；共生体系的长期维持引发了一系列的进化问题，也因此成为进化理论的一大挑战。榕和榕小蜂系统是目前所知的最严格的专性共生体系之一，两者异常丰富的物种多样性和宿主生活型为共生机理研究提供了理想的体系。榕-榕小蜂的共生关系古老且多样，早期的形态学研究显示1种榕往往只有1种对应的传粉小蜂，分子系统学研究也证实在较高的分类等级上两者间具有高度的协同枝系发生和协同进化关系，这种协同关系被应用到物种水平，“一对一”原则被认为在该系统中具有普适性。然而，近年来越来越多的案例证实许多榕属物种拥有1种以上的传粉小蜂，也有同一传粉小蜂为不同榕授粉的实例。特别是分子标记的应用，揭示出大量传粉小蜂隐存种，比例甚至达到所涉及榕属物种的60％；加之宿主转移现象和榕天然杂交后代的发现，冲击着榕-榕小蜂在物种水平高度对应的观点，Machado等也因此提出一个更为松散的进化模式，认为榕与榕小蜂的协同进化仅发生在两者群组之间，在榕群组内宿主转移现象可能经常发生。种间专一性屏障的松弛将极大地影响共生双方的动态，哪怕是低频率的例外事件也会导致严重的进化后果。Machado等的进化新模式主要依据巴拿马和墨西哥2个榕组及其对应的传粉小蜂，但有研究证实在某些地区和某些榕属类群中，榕-榕小蜂在物种水平高度对应，协同成种在其进化中占主导地位。因此要正确评价种间专一性在榕-榕小蜂共生体系中的真正地位，不同榕属类群和不同地区非一对一事件的发生频度及产生原因是不可或缺的基本信息。薜荔小蜂为薜荔及其变种爱玉子的专性传粉小蜂，至今形态学研究仍未发现它们在任何地区有差异，但薜荔和薜荔小蜂的分类地位飘浮不定，不同的研究者将它们归入不同的属，且在爱玉子人工栽培中薜荔小蜂能否为之有效传粉也有不同的结果。为探究薜荔传粉小蜂的物种组成、地域分布和系统发生关系，我们应用线粒体COI和Cytb、核核糖体ITS2基因对我国东南沿海岛屿及邻近大陆31个种群，331个个体进行了遗传组分、遗传多样性和种群扩张动态分析，并在GenBank中检索相关序列联合构建系统树，以期揭示薜荔传粉小蜂在该地区的物种组成、形成原因和可能的扩散路线；评价薜荔与其传粉小蜂间的物种专一性；评判岛屿和大陆生境对小蜂种群动态的影响；为印度-澳洲区系榕-榕小蜂专性对应程度提供案例；为协同进化模式研究提供参考。论文主要结论如下：1)薜荔传粉小蜂在我国东南沿海地区的遗传组成并不均一，可以划分为3个独立的枝系。3个枝系无论是线粒体基因还是核基因都超出了种内差异的范畴，达到物种分异水平；3个枝系没有明显的形态差异，为3个隐存种，分别为Sp．A、Sp．B、Sp．C。不同基因和不同方法均将31个种群的331个个体区分为3个枝系，枝系间的遗传分异明确。STUCTRUE软件模拟100万次将个体分配到各枝系的概率都在95％以上；系统树中3个单系的支持率都达99％；在巢式枝系结构中，3个枝系在95％的置信区间内完全分离。不同基因和不同方法区分得到的3个枝系在单倍型和个体水平都完全对应。比较3个枝系COI序列分化显示，平均枝系内K2P距离为0．57%，而枝系间分别达6．63％，10．94％和10．69％，任2个枝系间的序列分异都达到平均枝系内遗传差异的10倍以上。这不仅超过了区分动物物种的COI barcode标准距离域，也高于类似研究中榕小蜂隐存种对间的遗传距离，其线粒体基因达到物种分化水平。3个枝系的ITS2序列显示出致同进化的特征，枝系内部基本共享相同的ITS2单倍型，而不同枝系对应的单倍型间相差7-13个核苷酸。尽管枝系三有2个对应的ITS2单倍型，但它们仅相差1个碱基，且其中1个单倍型仅局限在1个南方种群。ITS2的单倍型格局从核基因水平证实3个枝系的物种性质。3个枝系达到物种分化水平的核基因证据还来源于Cytb的NUMTs序列，在Cytb单倍型和NUMTs的联合系统树中近期的NUMTs序列没有枝系的交叉，证实枝系间已完全隔离。2)薜荔传粉小蜂3个物种互为姊妹种。在薜荔小蜂COI和Cytb单倍型与其他小蜂相应序列联合构建的系统树中，薜荔小蜂序列聚为1个大单系的支持率为99%；2个物种中同时存在Cytb基因的NUMTs序列也说明它们为近缘种。联合Cytb单倍型和NUMTs序列进行系统重建，NUMTs古老单倍型聚在外类群位置，以之为根的系统树进一步证实了3个物种近缘的系统发生关系；而来源于不同物种的NUMTs序列在外类群中聚为单系更进一步确认它们互为姊妹种的性质。3)薜荔传粉小蜂3个物种在我们东南沿海地区有一定的地域分布和物候差异，生境偏好、地理和时间隔离可能是其物种分化的原因。薜荔传粉小蜂Sp．A完全分布于北部种群，Sp．B主要分布于南部种群，而Sp．C主要分布于北部的岛屿种群。地域分布差异反映了3个物种的生境偏好，纬度、地形地貌和岛屿大陆环境等差异导致的不同选择压力，山系和海洋导致的隔离都在薜荔小蜂祖先种群的分化中起到重要作用，小蜂先于薜荔的成种事件最终导致3种传粉小蜂在不同地域的共存。3个物种的物候有一定差异，Sp．C的出飞期较Sp．A晚10天左右，这种差异甚至逆纬度梯度存在，尽管这个差异在目前条件下尚不足以隔离物种，但在成种过程中起到一定作用。4)Sp．A非随机分布，呈现出明显的大陆-岛屿格局，岛屿群来源于大陆群的向外辐射。Sp．A广泛分布于北部的27个种群，STUCTURE和NCPA软件将之区分为大陆和岛屿2个下级组，X2置换检验表明2组非随机分布，AMOVA分析揭示大陆群与岛屿群间的差异达到显著水平，它们在单倍型组成、遗传多样性和种群动态上都有不同。与大陆群相比较，岛屿群表现出更明显的种群扩张表征，在巢式枝系图中呈现出更典型的星状结构、更多的单倍型末枝。失配分布分析和统计扩张检验都显示岛屿群经历过或正经历着扩张和瓶颈效应，而大陆群的Tajima’s D和Fs统计量都没有显著偏离0。比较遗传多样性可以看到，尽管岛屿群的样本量为大陆群的3倍，但其遗传多样性无论在整体水平还是种群平均水平均低于大陆群，特别是平均序列差异。COI岛屿特有单倍型间最多相差3个点突变，而大陆特有单倍型间可达12个。系统重建也显示大陆单倍型更为古老，它们在有根树中都深埋于内部节点。从COI单倍型格局看，所有岛屿特有单倍型都以1个点突变围绕着2个频率最高的大陆岛屿共享单倍型，且都分布在相应大陆种群的东北方向，这不仅符合小蜂的生活史特性，也与杭州湾风力场特征一致。岛屿生境剧烈波动和小种群效应导致的更频繁的岛屿种群灭绝事件、大陆种群的远距离救援、以及其后的定居、种群扩张和相对独立条件下的进化，构成了Sp．A典型的大陆-岛屿格局。5)宁德及其以南地区为Sp．B和Sp．C的潜在冰期避难所。Sp．B和Sp．C在南方和北方种群都有分布，但南方种群表现出更高的遗传多样性，拥有更多的古老特有单倍型。就物种C而言，宁德具有最多的特有单倍型，3个特有单倍型间相差8-11个点突变；而在其他11个种群中，仅桃花岛有1个特有单倍型，后者的特有单倍型与分布最广泛的共享单倍型间只相差1个点突变。从种群动态看，宁德为长期存在种群，失配分析和中性检验的各统计量都显著拒绝其经历过近期统计扩张的假设。宁德及其以南区域很可能为物种B、C潜在的冰期避难所，该地区的武夷山系为薜荔小蜂冰期种群提供了保护地，冰期后，小蜂沿山系从西南向东北方向扩散，直至薜荔分布区北缘。

论文目录

• 摘要
• Abstract
• 第一章 前言
• 1.1 榕传粉小蜂与其共生宿主的专一性对应关系
• 1.2 隐存种的揭示
• 1.3 薜荔及薜荔小蜂简介
• 1.4 论文的主要内容及研究路线
• 第二章 材料与方法
• 2.1 种群的选择及样品的采集
• 2.2 薜荔小蜂总DNA的提取
• 2.3 扩增测序
• 2.3.1 线粒体序列的直接测序
• 2.3.2 线粒体序列的克隆测序
• 2.3.3 核核糖体序列的直接测序
• 2.3.4 序列的拼接及比对
• 2.4 数据分析
• 2.4.1 序列特征的分析
• 2.4.2 遗传枝系的划分
• 2.4.3 枝系间遗传差异的计算
• 2.4.4 枝系间系统发生关系的确定
• 2.4.5 种群扩张的检验
• 2.4.6 遗传多样性的分析
• 第三章 结果
• 3.1 测序结果及序列特点
• 3.1.1 COI序列特征
• 3.1.2 Cytb序列特征
• 3.2 单倍型频率及分布
• 3.2.1 COI单倍型频率及分布
• 3.2.2 Cytb单倍型频率及分布
• 3.2.3 ITS2单倍型
• 3.3 遗传枝系
• 3.3.1 枝系的划分
• 3.3.2 枝系间的遗传差异
• 3.3.3 枝系的地域分布
• 3.3.4 枝系的系统发生关系
• 3.4 基于COI序列的枝系及种群扩张历史
• 3.4.1 枝系一的失配分析及中性检验
• 3.4.2 枝系二的失配分析及中性检验
• 3.4.3 枝系三的失配分析及中性检验
• 3.5.枝系及种群的COI序列遗传多样性
• 3.5.1 枝系遗传多样性
• 3.5.2 种群遗传多样性
• 第四章 讨论与结论
• 4.1 序列性质
• 4.1.1 薜荔传粉小蜂COI和Cytb清晰序列的mtDNA性质
• 4.1.2 薜荔小蜂Cytb克隆序列的核基因性质
• 4.2 薜荔传粉小蜂的物种组成
• 4.3 薜荔与薜荔传粉小蜂的种间专一性
• 4.4 薜荔小蜂3个物种的遗传多样性格局及其历史渊源推测
• 4.4.1 物种A的大陆岛屿格局
• 4.4.2 物种B、C的潜在冰期避难所
• 4.5 结论
• 附录
• 参考文献
• 发表论文
• 后记