论文部分内容阅读
本论文采集了中华沙塘鳢(Odontobutis sinensis)的9个种群,用线粒体(Cyt b)及核基因(EPICs、微卫星)多种分子标记,全面分析了长江及珠江流域的中华沙塘鳢种群遗传结构,探讨了这种遗传差异的成因。主要结果如下:1.Cyt b基因对长江及珠江流域中华沙塘鳢的遗传多样性分析从宜都(YD)、当阳(DY)、荆州(JZ)、洪湖(HH)、鹰潭(YT)、桂林(GL)、鹰潭(DTH)、全州(QZ)、湘潭(XT)等9个地区采集108尾中华沙塘鳢个体,扩增线粒体Cyt b基因全长用于序列分析,并选取其中的1127 bp进行遗传学分析,结果如下:(1)从108尾中华沙塘鳢样本中,共获得29个Cyt b基因单倍型(序列长度1 127 bp),总的单倍型多样性(Hd)、核苷酸多样性(π)值分别为0.929 0、0.009 41,呈现出较高的遗传多样性和较低的核苷酸多样性的特点;(2)中华沙塘鳢9种群间的FST值在0.110 6~0.998 8间(P<0.01),K2-P遗传距离在0.002~0.022间,揭示各种群间存在显著的遗传分化;而多数种群间的基因交流值(Nm)小于1(P>0.05),表明这些中华沙塘鳢种群间的基因交流有限。(3)中性进化检测和网络亲缘关系分析表明,中华沙塘鳢的桂林(GL)和洞庭湖(DTH)种群经历过种群扩张事件,推测中华沙塘鳢基因交流的方向为由长江水系种群流向珠江水系种群。通过以上研究为中华沙塘鳢的资源保护、开发与利用提供科学的理论依据。2基于核基因EPICs分子标记的种群遗传多样性及遗传结构分析本研究克隆并测序了一个EPICs标记即55305E1的基因序列,用于中华沙塘鳢9个种群的遗传多样性和遗传结构分析。(1)序列分析。研究中,我们选取55305E1基因的1070bp,对长江及珠江流域共9个种群进行遗传多样性分析。结果表明,共有178个变异位点,变异比例为16.7%,其中单变异位点有105个,占总变异的59%,简约信息位点为63且占总变异的35%。碱基的平均含量为T32.6%,C 18.4%,A 31.7%,G17.7%。A+T(64.3%)>G+C(36.1%),表明中华沙塘鳢中线粒体与核基因均呈现AT偏好性。(2)遗传多样性分析。中华沙塘鳢9个种群的种群总的单倍型多样性(Hd)为0.99330,变化范围在0.6-1之间。核苷酸多样性(π)为0.01370,变化范围为在0.0035-0.0132之间。(3)种群亲缘关系及遗传结构分析。从9个种群的108个个体中,共检测到82个单倍型(H)。每个种群都有多个单倍型,其中DTH单倍型最多(19个),而GL单倍型最少(4个)较少。其中H3、H14为GL和QZ的共享单倍型,H13为GL、QZ、XT的共享单倍型。系统进化树显示长江及珠江流域中华沙塘鳢共分为两大支,即:YD的种群单独聚为一支,而其它的8个种群聚为另外一支,包GL、DTH、XT、QZ、YD、JZ、HH、YT种群。珠江流域的GL与QZ先聚类,然后再与XT和DTH聚类。单倍型网状结构关系显示该9个种群并没有按照地理位置进行聚类,表明中华沙塘鳢种群处于谱系重排的状态。AMOVA分析中华沙塘鳢最佳的分组方法,即YD种群单独分为一组而其它群体分为一组为最好的分组方法。(4)中华沙塘鳢种群历史动态检测。分析发现,Tajima’s D、Fu’s Fs小于0(P<0.05),加上单倍型最小网状图显示出星状散布的特点,表明长江及珠江流域的中华沙塘鳢可能经历了种群扩张事件。3、基于微卫星的中华沙塘鳢种群遗传多样性与遗传结构分析本章用13对微卫星引物,分析了以上9个种群的遗传多样性与遗传结构。主要结果如下:(1)种群遗传多样性。NA种群平均等位基因(Na)数为6,其中XT的等位基因数最多,而YD等位基因数最少。QZ的等位基因丰富度(AR)最高,而YD的最低。所有种群的观测杂合度(Ho)为0.347,期望杂合度(HE)为0.601。其中GL的Ho最高(0.471),而DY的Ho最低(0.212),而GL的HE最高(0.588),XT的最低(0.398)。各种群的近交系数Fis变化范围为0.713-0.745,都大于0.5,说明中华沙塘鳢的9个种群间分化较大,种群总的遗传多样性较高。(2)种群遗传差异分析。分析发现,种群间遗传距离较大,如DY和YT遗传距离最大为0.0458,而最小的为GL和QZ为0.101。不同种群的AMOVA分析显示,无论是将种群按地理距离分组还是系统进化树分支分组,种群间遗传差异都大于种群内的遗传差异,表明差异主要来源于种群间。NJ树显示,GL、QZ、XT、DTH聚为一类,而JZ、HH、YT、DY、YD聚为一类,显示出按照地理位置聚类的特征。(3)种群主成分分析。GL、QZ、XT、DTH聚为一组,而JZ、HH、YT、DY、YD聚为一组,且PC1为34.94%,PC2为20.87%。Structure显示K=2时为最佳分组。(4)种群的瓶颈效应。根据种群是否杂合度过剩来判断种群是否经历过瓶颈效应,(Sign test和Wilcoxon’s test)检测模型发现GL、DY、JZ种群经历显著的瓶颈效应。