新疆兔属物种的群体遗传学及系统发育关系

论文摘要

兔属物种由于适应性强,生存环境多样,其群体遗传学和系统发育研究对于揭示物种进化、探索物种与地理环境的适应性等热点问题具有重要意义。中国的兔属动物共有9种,新疆分布有三种:塔里木兔（Lepus yarkandensis Gunther 1875）,草兔（Lepus capensis Linnaeus, 1758）和雪兔（Lepus timidus Linnaeus, 1758）。目前,针对这三个物种的群体遗传学以及系统发育研究有限。塔里木兔和雪兔为国家二级保护动物,其中塔里木兔为我国特有种,仅分布于新疆塔里木盆地周围分散的绿洲,栖息地片段化程度日益严重。本研究选用两个线粒体基因（mtDNA）标记:控制区（D-loop）和细胞色素b（Cytb）基因,以检测塔里木兔种群的遗传多样性水平、基因交流状况;分析地理隔离对该物种群体遗传结构的作用;评价第四纪气候变化对塔里木兔遗传多样性及群体历史的影响。本研究测定了来自塔里木盆地周围二十个采集地224个塔里木兔样本的D-loop区553bp序列和Cytb全长1140bp的序列。通过对所得基因的序列特征分析,结合Mantel tests分析和多层次分子变异分析（AMOVA）,发现塔里木兔群体总体的mtDNA核苷酸多样度与其它哺乳动物的相比属于较低水平（D-loop 3.3%和Cytb 0.8%）,各地理种群具有丰富的特有单倍型,种群间已经出现显著分化,这种遗传分化和单倍型的地理分布模式是由于地理隔离（河流、绿洲、沙漠、戈壁）及栖息地片段化（绿洲退缩及分散分布）而导致的种群间有限的基因流造成的。对这两个多态性遗传标记的单独及合并数据分别用邻接法（Neighbor- joining;NJ）和贝叶斯法（Bayesian inference）构建系统发育树,结果显示:224个塔里木兔样品分为五个进化枝,每个进化枝包含的样品分布区稍显混杂,然而,当将这20个采样地划分为东西南北四个地理种群时,西南部群体与东北部群体便呈现出明显的遗传分化,合并数据的中介网络图更清晰地反映了这种系统地理分化格局。种群动力学及遗传多样性数据显示,塔里木盆地西南部可能是塔里木兔在最大冰期时的避难所,分别发生在0.21 MYA、0.090 MYA及0.054 MYA的三次冰期后扩张事件导致塔里木兔现今的分布格局。值得注意的是,不是冰期,而是水源,才是塔里木兔退缩至避难所的主要推动力。从系统树及实际的地理隔离看,塔里木兔可能经历了从塔里木盆地西南部到北部再到东部的扩散过程。种群历史过程及冰期避难所的研究也为塔里木兔生物多样性保护重点区域的划分提供了依据。基于本文线粒体基因的研究结果,塔里木兔五个进化枝具有较高的支持率和较多的突变步数,各进化枝间的遗传距离（0.7%-2.3%）达到亚种水平,西南组群与东北组群间分歧显著,均提示塔里木兔可能存在亚种的分化。但是,这种亚种划分仅基于我们线粒体基因的研究结果,是否有效还需要结合核基因、形态学以及生态学方面的证据。新疆草兔（Lepus capensis）的群体遗传结构至今无系统的研究报道,亚种水平的分类也长期存在争议。本文测定了形态分类上的新疆草兔三个亚种共74个个体的mtDNA控制区553bp和Cytb全长1140bp的序列。新疆草兔核苷酸多样度为0.021±0.010;单倍型多样度为0.974±0.007。Mantel tests和AMOVA分析显示,新疆草兔除东部与中部种群间没有明显分化外,其余种群间已产生十分明显的分化,种群间的基因交流有限,可能是由于山脉戈壁等地理隔离而造成的。新疆草兔群体的中性检验和核苷酸错配分析显示多峰分布,这主要是由新疆草兔群体存在较为明显的遗传分化而致。基于新疆草兔mtDNA D-loop和Cytb的合并数据所构的NJ树和贝叶斯树以较高的支持率聚为三枝,每枝包含了来自相同或相邻地理区域的草兔样品,显示出了一定的地理分布格局。中介网络图也得到了相同的结果。本研究的结果支持形态分类上草兔西域亚种（L.c. lehmanni）的分类地位;但中亚亚种（L.c. centrasiaticus）被分为两个独立的进化枝,提示可能存在两个亚种;帕米尔亚种（L.c. pamirensis）可能已达到种的分化水平。目前关于新疆兔属物种间的系统发育关系还不清楚。本研究采用线粒体两个基因（D-loop和Cytb）以及一个核基因（MGF）标记来探讨新疆兔属物种间的系统发育关系。线粒体系统发育树的构建采用邻接法和贝叶斯法,核基因采用邻接法和最小进化法（Minimum Evolution; ME）。我们共测定了353个新疆野兔样本的线粒体控制区553bp和Cytb全长1140bp的序列;由于存在高比例的杂合子,通过分子克隆共得到了304个个体的453条核基因MGF序列（592bp）。从线粒体的系统发育树上可以看出,新疆的草兔和塔里木兔序列分为3个世系:其中两个世系分别代表塔里木兔与草兔各自特有的世系,另一个世系为新发现的世系。新世系包含了来自塔里木兔及草兔采集地的样品。三个世系中草兔最先分歧出来,塔里木兔与新世系形成姐妹群的关系。唯一的一条新疆雪兔线粒体序列包含在草兔的世系中。核基因MGF的分析结果表明,草兔和塔里木兔序列在核基因进化树上相互混杂并大致分为两大类群。雪兔的核基因单独形成一枝,并且最先分歧出来。线粒体与核基因的结果相矛盾说明新疆草兔、塔里木兔以及雪兔之间可能存在杂交,但是进一步的证据还需要来自形态方面如外形及头骨的考查。

论文目录

• 摘要
• Abstract
• 第一章 文献综述
• 1.1 兔属动物分类及系统学研究进展
• 1.1.1 兔属动物的分类
• 1.1.2 兔属动物系统进化关系研究进展
• 1.1.3 中国野兔分类
• 1.1.4 中国野兔的地理区划
• 1.1.5 兔属起源
• 1.1.6 兔类研究工作的展望及存在的问题
• 1.2 分子系统地理学研究进展
• 1.2.1 发展简史及主要研究工具
• 1.2.2 研究内容和应用
• 1.2.3 前景
• 1.3 现代系统发育学研究常用的系统发育重建方法
• 1.3.1 最大简约法（Maximum Parsimony）
• 1.3.2 最大似然法（Maximum Likelihood）
• 1.3.3 距离法（Distance Methods）
• 1.3.4 系统树可靠性检验
• 1.3.5 贝叶斯法（Bayesian Methods）
• 1.3.6 小结
• 第二章 塔里木兔的群体遗传学研究
• 2.1 序言
• 2.2 材料和方法
• 2.2.1 样品收集
• 2.2.2 总DNA 提取
• 2.2.3 PCR 引物设计
• 2.2.4 PCR 扩增及测序测定
• 2.2.5 排序及数据统计分析
• 2.2.6 群体分析
• 2.2.7 系统分析及外群的选择
• 2.3 结果
• 2.3.1 序列特征
• 2.3.2 系统发育分析
• 2.3.3 中介网络图的构建
• 2.3.4 单倍型分布和群体遗传多样性
• 2.3.5 种群系统地理结构
• 2.4 讨论
• 2.4.1 塔里木兔群体遗传结构
• 2.4.2 塔里木兔群体的系统地理结构
• 2.4.3 冰期避难所及冰期后扩张
• 2.4.4 塔里木兔的亚种分化
• 2.4.5 保护管理建议
• 第三章 新疆草兔的群体遗传学研究
• 3.1 研究背景
• 3.2 材料和方法
• 3.2.1 样品收集
• 3.2.2 总DNA 提取
• 3.2.3 PCR 引物设计
• 3.2.4 PCR 扩增及测序测定
• 3.2.5 排序及数据统计分析
• 3.2.6 群体分析
• 3.2.7 系统分析及外群的选择
• 3.3 结果
• 3.3.1 序列特征
• 3.3.2 系统发育分析
• 3.3.3 中介网络图
• 3.3.4 草兔群体的遗传多样性及遗传结构
• 3.4 讨论
• 3.4.1 新疆草兔的遗传多样性
• 3.4.2 新疆草兔的遗传分化与遗传结构
• 3.4.3 亚种分化
• 第四章 新疆兔属动物的系统发育、分类及基因渗透
• 4.1 研究背景
• 4.2 材料和方法
• 4.2.1 样品收集
• 4.2.2 总DNA 提取
• 4.2.3 PCR 引物设计
• 4.2.4 PCR 扩增及测序测定
• 4.2.5 克隆
• 4.2.6 排序及数据统计分析
• 4.2.7 系统分析及外群的选择
• 4.3 结果
• 4.3.1 序列特征
• 4.3.2 基于D-loop 构建的系统发育关系
• 4.3.3 基于Cytb 构建的系统发育关系
• 4.3.4 基于D-loop 和Cytb 合并数据构建的系统发育关系
• 4.3.5 基于mtDNA 单独及合并数据的新疆野兔系统发育关系小结
• 4.3.6 基于核基因序列的系统发育
• 4.3.7 基于mtDNA 与核基因序列的新疆野兔系统发育关系小结
• 4.3.8 基因渗透（Introgression）
• 4.4 讨论
• 4.4.1 新疆野兔的亲缘关系与分类地位
• 4.4.2 基因渗透（Introgression）
• 第五章 结论与展望
• 参考文献
• 附录1
• 附录2
• 附录3
• 博士期间主持的课题和发表的论文
• 致谢