长江中上游圆筒吻鮈群体线粒体Cyt b遗传多样性分析

长江中上游圆筒吻鮈群体线粒体Cyt b遗传多样性分析
蒲艳;田辉伍;陈大庆;段辛斌;刘绍平;汪登强
【摘要】为探明大坝建设及过度捕捞等可能会对圆筒吻鮈资源造成的影响,采用线粒体DNA Cyt b基因序列对长江中上游宜宾、泸州、合江、江津、巴南、涪陵和黄冈7个圆筒吻鮈(Rhinogobio cylindricus)群体进行研究,分析其遗传多样性及遗传结构.结果表明:在获得的998 bp序列长度中检测到14个变异位点,140尾样本定义15个单倍型,平均单倍型多样性(Hd)和平均核苷酸多样性(Pi)分别为
0.625、0.00077,说明圆筒吻鮈种群遗传多样性较贫乏.分子变异方差分析(AMOVA)表明圆筒吻鮈地理群体间未出现显著遗传分化.单倍型Network网络分析呈星型结构,邻接系统发育(NJ)树显示单倍型没有出现高支持率的分支,表明群体内部没有发生遗传分化.中性检验、错配分布及BSP分析均表明圆筒吻鮈发生历史种群扩张事件,种群扩张发生在0.06 Ma(百万年前).综上表明,圆筒吻鮈的群体遗传多样性较低,群体间遗传分化不显著,建议将圆筒吻鮈群体作为一个整体就地保护,并在今后应加强长期监测.
【期刊名称】《淡水渔业》
【年(卷),期】2019(049)001
【总页数】6页(P14-19)
【关键词】圆筒吻鮈(Rhinogobio cylindricus);线粒体Cytb;遗传多样性;长江中上游
【作者】蒲艳;田辉伍;陈大庆;段辛斌;刘绍平;汪登强
【作者单位】西南大学生命科学学院, 重庆400715;中国水产科学研究院长江水产研究所, 农业部长江中上游渔业资源环境科学观测实验站, 武汉430223;中国水产
科学研究院长江水产研究所, 农业部长江中上游渔业资源环境科学观测实验站, 武
汉430223;中国水产科学研究院长江水产研究所, 农业部长江中上游渔业资源环境科学观测实验站, 武汉430223;中国水产科学研究院长江水产研究所, 农业部长江
中上游渔业资源环境科学观测实验站, 武汉430223;中国水产科学研究院长江水产研究所, 农业部长江中上游渔业资源环境科学观测实验站, 武汉430223;中国水产
科学研究院长江水产研究所, 农业部长江中上游渔业资源环境科学观测实验站, 武
汉430223
【正文语种】中文
【中图分类】S917.4
圆筒吻鮈(Rhinogobio cylindricus)，隶属于鲤形目(Cypriniformes)鲤科(Cyprinidae)鮈亚科(Gobioninae)吻鮈属(Rhinogobio)，主要分布在长江上游干
流及其支流(如岷江、金沙江、嘉陵江等)[1-2]，近年有研究表明在长江中游也有分布[3]，是长江主要经济鱼类 [4]。

圆筒吻鮈喜生活于流水或缓流水环境，有洄游产卵行为[5]，主要以底栖生物(如摇蚊幼虫)为食[6]。

资源调查表明，目前长江上游
圆筒吻鮈资源开发率偏高，存在过度开发现象，资源量呈下降趋势[9]。

圆筒吻鮈
产漂流性鱼卵，受精卵随江水漂流发育，长江干流水利工程建设改变了河流自然水文状况、流水江段长度变短、阻隔了鱼类的洄游通道[4]，严重影响了其种质资源，大坝建设对圆筒吻鮈繁殖习性同样产生严重影响，进一步威胁其生存。

迄今，有关圆筒吻鮈研究报道主要在生物学[4]、物种分化[7]、组织学[8]、资源量[9]等方面，在遗传学方面，曾晓芸对其遗传结构分析也有了初步研究[10]，其研究区域为长江
上游宜宾至万州库区江段以及主要支流，属于长江上游及其支流，而长江中游却未见相关报道。

线粒体DNA(Mitochondrial DNA，mtDNA)因结构简单、呈母系遗传、拷贝数
量大、进化速率快、无组织特异性等特点，被广泛运用于物种分化、群体遗传多样性[10，19]、种属亲缘关系远近等研究中[11]。

细胞色素b(Cyt b)基因的结构和功能在mtDNA 13个蛋白质编码基因中了解的最为全面，且其进化速度适中、转换或颠换占很大比例等特点，被认为是研究种内或近缘中间系统发育、种质资源状况和群体遗传结构的理想工具，成为鱼类遗传多样性研究中常见的分子标记[14，26，27]。

本实验采用线粒体Cyt b序列分析长江中上游圆筒吻鮈群体的遗传多样性及
种群结构，以阐明其遗传多样性和群体分化状况，为圆筒吻鮈种质资源开发和保护提供科学依据。

1 材料与方法
1.1 材料
2017年7～10月在长江上游干流的宜宾(YB)、泸州(LZ)、合江(HJ)、江津(JJ)、巴南(BN)、涪陵(FL)及长江中游的黄冈(HG)7个采样点共采集圆筒吻鮈140尾(表1)，各采样点位置见图1。

所有样本根据其形态鉴定后剪取鳍条，保存于无水乙醇。

图1 采样点地理位置Fig.1 Sampling stations of R.cylindricus in the Yangtze River
1.2 基因组DNA提取、PCR扩增及测序
剪取适量鳍条，双蒸水洗去乙醇，用高盐法提取总DNA[12]。

用鱼类通用引物扩
增线粒体DNA Cyt b序列，引物序列为L14724：5′-GACTTGAAAAACCACCGTTG-3′和H15915：5′-CTCCGATCTCCGGATTACAAGAC-3′[13]。

PCR 扩增体系50 μL，包括：10× buffer(Mg2+)5 μL，dNTPs(10 md/L)0.5 μL，正反引物各2 μL，Taq DNA
Polymerase(2 U)0.4 μL，50 ng/μL DNA模板2 μL，灭菌去离子水补至50 μL。

反应程序如下：94 ℃预变性4 min；94 ℃变性30 s，52 ℃退火40 s，72 ℃延
伸1 min，35个循环；72 ℃再延伸8 min。

用1.0%琼脂糖凝胶电泳检测PCR 扩增产物，条带明亮、清晰、单一的样品，交由生物公司进行双向测序，测序引物与PCR扩增引物相同。

1.3 数据处理
采用Lasergene v7.1软件包的Seqman进行拼接，序列的对位排列用Clustal X
软件完成。

用MEGA v6.0分析软件统计序列的碱基组成、转换/颠换比率，并构
建邻接系统发育(NJ)树，树检验采用bootstrap法，设1 000次重复。

使用DnaSP v5.0获得单倍型分布，统计变异位点(s)、单倍型数目(h)、单倍型多样性
指数(Hd)和核苷酸多样性指数(Pi)，并进行中性检验和歧点分布分析及Nm值计算。

应用Arlequin v3.11软件进行AMOVA分析及群体的遗传分化指数FST。

使用Network v5.003程序构建单倍型网络结构图，对单倍型之间的进化关系进行分析。

根据公式τ=2 ut估算种群扩张时间，其中τ是扩张时间参数；u=2 μk，其中μ
为每个核苷酸的变异速率，k为序列长度；t表示自扩张以来所经历的代数；扩张
时间T=t×代时。

Beast 2.4.3.0进行Bayesian Skyline Plots(BSP)分析，采用鲤
科鱼类常用的进化速率每百万年1%[14]，MCMC长度为100 000 000，在TRACER ver.1.4构图，各项参数的ESS均大于200。

2 结果与分析
2.1 序列变异
经校正和比对后得到圆筒吻鮈线粒体Cyt b有效序列长度为998 bp，约为Cyt b
基因全长的87%。

序列中的转换数(Ts)明显高于颠换数(Tv)，Ts/Tv =6.01。

A、T、C、G的平均含量分别为27.64%，30.85%，27.11%和14.40%，A+T(58.49%)
大于C+G(41.51%)，并显示明显的反G含量倾向，与其他硬骨鱼类的情况相似。

在998 bp序列中发现14个变异位点，其中7个单一突变位点，另7个为简约信息位点，序列间平均核苷酸差异数(K)为0.772。

140个个体中共检测到15个单倍型(GenBank登录号：MH673879-MH673893)，其中单倍型H_4出现的频率最高，为58.6%，其次是H_2(13.6%)和H_3(12.1%)，有7个单倍型仅出现1次。

7个群体中，江津群体的单倍型数最多为13个，涪陵的单倍型数最少，仅有1个(表3)。

群体平均单倍型多样性指数和核苷酸多样性指数分别为0.625和0.000 77，最高的是泸州群体，其次是巴南群体，涪陵群体最低(表1)。

表1 圆筒吻鮈线粒体Cyt b遗传多样性Tab.1 Genetic diversity parameters within R.cylindricus based on mtDNA Cyt b sequences群体样本数(n)突变位
点数(s)单倍型数(h)单倍型多样性(Hd)核苷酸多样性(Pi)Tajima′s DFu′s
FsYB10230.5110.000 56-0.690 98-0.594LZ12560.8180.001 11-1.223 73-
3.022HJ23450.6050.000 71-0.951 63-1.830JJ6312130.6340.000 78-1.978
12-11.388BN20560.7580.001 01-0.865 88-2.229FL60100--
HG6220.3330.000 67-1.131 970.952合计14014150.6250.000 77-1.821 88-12.199
2.2 群体遗传结构
群体间分子变异方差分析(AMOVA)显示，群体主要变异来源为群体内，占
101.00%，群体间的遗传分化指数FST为 -0.998%(P>0.1)(表2)，表明圆筒吻鮈
7个群体间未出现显著遗传分化。

表2 圆筒吻鮈群体间分子变异分析Tab.2 Analysis of molecular
variance(AMOVA)among R.cylindricus populations变异来源自由度方差变异
组成变异百分比群体间61.939-0.003 85-1.00群体内13351.7400.389 02101.00总变异13953.6790.385 17遗传分化指数(FST)-0.009 98
对两两群体间的FST值进行计算(表4)，种群间FST值均低于0.05，说明群体之
间未出现遗传分化。

Nm值表示群体间基因交流程度，其Nm>1则表明群体间存在基因交流[15]，经计算总群体Nm值为115.80，两两群体Nm值见表4。

表3 圆筒吻鮈群体线粒体Cyt b序列单倍型数目及分布Tab.3 MtDNA Cyt b haplotype distribution across R.cylindricus populations单倍型YBLZHJJJBNFLHG总计
H_111215H_21249319H_32135617H_475143786582H_5112H_6112H_711H _8112H_911H_10224H_1111H_1211H_1311H_1411H_1511
表4 圆筒吻鮈群体间的FST值(对角线下方)和基因流Nm值(对角线上方)Tab.4 Pairwise FST(below diagonal)and Nm values(above diagonal)among
R.cylindricus populations群体YBLZHJJJBNFLHGYB-14.69-5.43-12.32-
5.923.13-3.44LZ-0.022 6-10.78-10.7731.604.5687.75HJ-0.053 3-0.015 6-12.93-48.132.50-21.79JJ-0.031 6-0.012 9-0.020 69.283.95-9
6.07BN-0.054 70.009 5-0.004 10.032 71.33-39.83FL0.006 6-0.025 2-0.017 3-0.050 90.040 0-HG-0.073 0-0.017 8-0.015 6-0.013 5-0.028 70
利用NETWORK里Nedian Joining方法构建圆筒吻鮈Cyt b序列的单倍型网络结构如图2所示。

圆筒吻鮈单倍型网络结构图以H_4为中心呈星状结构，相邻单倍型之间均是通过一个突变步骤连接。

以长鳍吻鮈(R.ventralis)为外类群，采用NJ法(neighbor-joining)构建的单倍型系统发育树如图3所示，系统树未出现高支持率(均低于60%)的分支，说明群体内部没有明显遗传分化。

图2 基于线粒体Cyt b序列构建圆筒吻鮈单倍型网络结构图Fig.2 Haplotype network of R.cylindricus based on the mtDNA Cyt b sequences
图3 圆筒吻鮈线粒体Cyt b 单倍型NJ树Fig.3 NJ phylogenetic tree of mtDNA Cyt b in R.cylindricus
2.3 历史种群动态
种群扩张历史动态利用Tajima′s D和Fu′s Fs中性检验和碱基错配来评估。

中性
检验的Tajima′s D=-1.821 88(P<0.05)及Fu′s Fs =-12.199(P<0.05)，均为显著
性负值。

单倍型错配分布分析显示圆筒吻鮈群体的歧点分布表现为单峰(图4)且BSP分析结果均表明圆筒吻鮈近期发生过经历过种群扩张。

鲤科鱼类常用的进化
速率每百万年1%[14]，根据τ=1.179，计算圆筒吻鮈在距今0.06 Ma(百万年)前
发生过种群扩张。

图4 圆筒吻鮈Cyt b序列单倍型错配分布图及BSP分析Fig.4 Mismatch distribution(A)and BSP analysis(B)of R.cylindricus based on the mtDNA Cyt
b sequences
3 讨论
3.1 遗传多样性
遗传多样性是生物多样性的总和，存在于生物种群间、种群内，是生物对环境适应性和进化的基础，物种遗传多样性的高低通常与该物种对环境的适应能力具有一定的线性关系，表现为遗传多样性越高，则表明对环境的适应能力越强[16-17]。

单
倍型多样性指数(Hd)和核苷酸多样性指数(Pi)是衡量物种遗传多样性水平高低的两个重要指标[18]，Pi值越大则表示群体多态性较高，反之亦然[16]。

与长江上游其它鱼类，如红唇薄鳅(L.rubrilaris)[19](Pi=0.006)、长鳍吻鮈[20](Pi=0.001)、异
鳔鳅鮀(X.boulengeri)[21](Pi=0.002)等比较，长江中上游圆筒吻鮈群体遗传多样
性处于较低水平。

本研究以前，曾晓芸[10]也采用Cyt b分析过长江上游圆筒吻鮈群体遗传结构，得到Pi=0.010 06，明显大于本研究的结果，暗示圆筒吻鮈群体遗传多样性存在下降趋势。

从单倍型来看，曾晓芸的研究结果中，频率最高的3个
单倍型(Hap_1：45.1%、Hap_2：20.5%和Hap_3：8.2%)占总样本量为73.8%，本研究的频率最高的3个单倍型比例为84.3%，提示近年来圆筒吻鮈群体同质化
有所提高，值得关注。

造成圆筒吻鮈遗传多样性下降的原因可能是过度捕捞造成种群大小下降，引起群体遗传漂变加剧。

研究表明，长江上游圆筒吻鮈资源开发率偏高[6，9]，资源量呈下降趋势。

研究还发现位于三峡库区的库尾涪陵群体没有检测到遗传变异，仅有一个单倍型(表3)。

圆筒吻鮈属喜流水性鱼类，涪陵江段在三峡水库低水位时江水存在一定的
流动性，高水位时，江水几乎不流动，因此不是圆筒吻鮈适宜的栖息地，但该江段的研究有助于分析大坝建设形成的水库对喜流水性鱼类的影响。

本研究还显示，长江中游的黄冈群体遗传多样性较低(表1)，有研究表明长江上游圆筒吻鮈资源量大
于中游[3，9]，说明长江上游较中游而言是圆筒吻鮈较适宜的栖息地，并且圆筒吻鮈在中游面临与更多大型鱼类竞争食物的压力等，这些可能是造成其群体遗传多样性降低的原因。

但由于这两个群体的样本量比较少，还需要加强采样，进一步研究。

3.2 遗传结构及种群历史
溯祖理论[22]认为分布最广泛的单倍型为原始单倍型。

基于Cyt b序列的单倍型网络结构图(图2)显示，H_4位于网络结构中心，为主要单倍型，在7个群体中均有分布，说明H_4单倍型的圆筒吻鮈对环境的适应能力较强。

AMOVA和基因流及
两两群体之间的遗传分化分析均表明长江中上游圆筒吻鮈群体遗传分化不显著(表
2和表4)，但基因流检测显示，涪陵群体与其它群体(泸州群体除外)的Nm均小于4，黄冈群体与其他群体(泸州群体除外)的Nm值也均小于4(表4)，说明这些群体间的基因交流可能存在一定障碍。

考虑到涪陵和黄冈群体样本量小，Cyt b序列的变异位点数量也较少，对分析的结果可能会造成一定的影响，因此有必要采集更多样本和采用多态性更高的标记，如mtDNA控制区[23]、微卫星DNA等进一步研究验证。

Grant等[24]根据海水鱼类单倍型多样性和核苷酸多样性的高低推测群体历史状况之间的关系，可分为四种情况：(1)低单倍型多样性(<0.5)和低核苷酸多样性
(<0.005)，说明种群近期发生瓶颈或奠定者由单个或少数mtDNA谱系组成；(2)
高单倍型多样性和低核苷酸多样性，说明群体种群发生瓶颈后迅速增长并积累突变；
(3)低单倍型多样性和高核苷酸多样性，说明地理隔离的亚群体出现分化；(4)高单
倍型多样性和高核苷酸多样性，说明具有长进化历史的稳定群体或有谱系差异的群体再次融合。

本研究结果显示，圆筒吻鮈群体遗传多样性与上述第二种类型相符，说明圆筒吻鮈群体发生遗传瓶颈和群体扩张历史，这与中性检验、碱基错配分布及BSP分析的结果吻合。

研究表明，长江上游的一些鱼类如红唇薄鳅、长薄鳅
(L.elongate)[19]、长鳍吻鮈[20] 在长江上游末次冰期，即大理亚冰期(1～11万)
发生过种群扩张，此时期长江上游气温回升，有利于鱼类的生存和扩张 [25]。

3.3 保护建议
本研究显示，长江中上游圆筒吻鮈群体没有发生遗传结构分化，与之前研究结果比较，其圆筒吻鮈群体遗传多样性存在下降趋势，需要开展长期遗传多样性监测。

近年来，渔民为了获得更多利益而改变了网具规格，不断缩小其网具的捕捞规格，资源量不断减少，有效繁殖群体难以恢复。

因此，应采取有效措施(如禁渔期、渔民
转业等)并加强对圆筒吻鮈的保护力度。

从本研究得知，7个群体之间及群体内部
都没有明显的遗传分化，因此，我们建议将圆筒吻鮈作为一个整体进行就地保护，并开展长期监测。

参考文献：
【相关文献】
[1]伍献文.中国鲤科鱼类志(上卷)[M].上海：上海科学技术出版社，1977.
[2]乐佩琦.中国动物志硬骨鱼纲、鲤形目(下卷)[M].北京：科学出版社，2000：273-390.
[3]林鹏程，刘焕章，曹文宣.长江宜昌江段三层流刺网对鱼类资源影响的分析[J].四川动物，2014，
33(5)：762-767.
[4]马惠钦.长江干流圆筒吻鮈生物学的初步研究[D].重庆：西南师范大学，2001.
[5]高天珩.长江上游鮈亚科鱼类资源及生境选择策略研究[D].重庆：西南大学，2016.
[6]熊星.圆筒吻鮈(Rhinogobio cylindricus Günther)基础生物学及资源量研究[D].重庆：重庆师范大学，2013.
[7]孙玉苗.吻鮈属鱼类物种分化的研究[D].济南：山东大学，2005.
[8]温龙岚，姚艳红，王志坚.吻鮈、圆筒吻鮈和福建紋胸鮡脾脏的组织学初步研究[J].遵义师范学院学报，2006，8(6)：49-51.
[9]熊飞，刘红艳，段辛斌，等.长江上游朱杨溪江段圆筒吻鮈种群参数和资源量[J].生态学报，2015，35(22)：7320-7327.
[10]曾晓芸.长江上游圆筒吻鮈遗传结构分析[D].重庆：西南大学，2015.
[11]杨子拓.珠江流域三种鲌亚科鱼类的遗传结构及遗传多样性分析[D].武汉：华中农业大学，2016.
[12] Aljanabi S M，Martinez I.Universal and rapid salt-extraction of high quality genomic DNA for PCR-based techniques [J].Nucleic Acids Res，1997，25(22)：4692-4693.
[13] Xiao W，Zhang Y，Liu H，et al.Molecular systematics of Xenocyprinae(Teleostei：Cyprinidae)：taxonomy，biogeography，and coevolution of a special group restricted in East Asia [J].Mol Phylogenet Evol，2001，18(2)：163-173.
[14]李小兵，唐琼英，俞丹，等.长江上游干流及赤水河蛇鮈遗传多样性与种群历史分析[J].动物学杂志，2016，51(5)：833-843.
[15]荣朝振，祖国掌，胡建华，等.泥鳅线粒体DNA控制区结构分析及遗传多样性研究[J].南方水产科学，2011，7(5)：55-62.
[16]胡静，侯新远，尹绍武，等.基于mtDNA COI和Cyt b基因序列对南中国海不同海域波纹唇鱼群体遗传多样性的研究[J].水生生物学报，2014，38(6)：1008-1016.
[17]田华.鲢鳙长江野生群体和养殖群体微卫星的遗传多样性分析[D].武汉：华中农业大学，2008.
[18]Vrijenhoek R C.Genetic diversity and fitness in small populations [J].Conserv Genet，1994，68：37-53.
[19]田辉伍.长江上游保护区长薄鳅和红唇薄鳅种群生态及遗传结构比较研究[D].重庆：西南大学，2013.
[20]程晓凤.长江上游特有鱼长鳍吻鮈遗传结构分析[D].重庆：西南大学，2013.
[21]董微微，汪登强，田辉伍，等.长江上游异鳔鳅鮀线粒体控制区遗传多样性[J].生态学杂志，2018，37(5)：1438-1443.
[22] Crandall K A，Templeton A R.Empirical tests of some predictions from coalescent theory with applications to intraspecific phylogeny reconstruction [J].Genetics，1993，134(3)：959-969.
[23]董志国，李晓英，王普力，等.基于线粒体D-loop基因的中国海三疣梭子蟹遗传多样性与遗传分化研究[J].水产学报，2013，37(9)：1304-1312.
[24] Grant W，Bowen B.Shallow population histories in deep evolutionary lineages of
marine fishes：insights from sardines and anchovies and lessons for conservation [J].J Hered，1998，89(5)：415-426.
[25] Bowen B W，Rocha L A，Grant D R，et al.Phylogeography of the
trumpetfishes(Aulostomus)：Ring species complex on a global scale [J].Evolution，2001，55(5)：1029-1039.
[26]张丽，滚双宝，雷天云，等.应用mtDNA Cyt b基因全序列分析中国5个马鹿群体的遗传多样性和系统发育[J].华北农学报，2010，25(4)：12-16.
[27]高天翔，毕潇潇，赵林林，等.基于线粒体Cyt b基因全序列的松江鲈群体遗传结构分析[J].水生生物学报，2013，37(2)：199-207.。