Genetic structure analysis of Ochetobius elongatus between Yangtze River and Pearl River using multiple loci
-
摘要: 为保护极濒危物种鳤 (Ochetobius elongatus) 的遗传多样性,通过野外调查在长江获得了7尾鳤样本,利用Sanger测序技术测定了其2个线粒体基因 [细胞色素b (Cytb) 和酰胺腺嘌呤二核苷酸 (NADH) 脱氧酶亚单位2 (ND2)] 和2个核基因 [肌球蛋白重链6 (myh6) 和重组激活基因2 (RAG2)] 的序列,结合已公开发表的52尾珠江鳤的4个基因序列,采用系统发育分析、单倍型网状图构建、分化时间估算等方法,对珠江和长江鳤的遗传结构和分化历史展开研究。系统发育分析和单倍型网状图表明,珠江和长江鳤群体分别形成了高度分化的遗传谱系,并且在核基因水平上形成了特有的等位基因型,表明这两个水系的鳤群体独立进化且不存在基因交流。分化时间估算表明珠江和长江鳤两谱系在0.38~0.76 Ma (百万年前) 分化,中更新世青藏高原的快速隆升很可能是促进其分化的重要因素。鉴于珠江和长江的鳤群体在线粒体基因和核基因层面上存在严格的地理分化,建议将这两个群体看作不同的演化显著单位 (Evolutionary Significant Unit, ESU) 进行区别管理和保护。Abstract: To protect the genetic diversity of a critically endangered species Ochetobius elongatus, taking seven individuals of O. elongatus collected from the field investigation in the Yangtze River as samples, we sequenced two mitochondrial genes (Cytb and ND2) and two nuclear genes (mhy6 and RAG2) for the seven samples by Sanger sequencing technique. Combining with the published four gene sequences of 52 O. elongatus samples in the Pearl River, we explored the genetic structure of O. elongatus between the Pearl River and the Yangtze River, so as to provide scientific support for its conservation. We applied phylogenetic analyses, haplotype networks and divergence time estimation. Phylogenetic analyses and haplotype networks reveal that O. elongatus populations in the two rivers generated two deep and independent lineages, and formed private alleles at the nuclear gene level. The results suggest that O. elongatus populations in the two rivers evolved independently without agene flow. The rapid lifting of the Qinghai-Xizang Plateau during the middle Pleistocene might be an important factor that triggered the split of O. elongatus populations in the two rivers 0.38~0.76 million years ago (Ma). In view of the strictly geographical division of O. elongatus populations between the two rivers at both mitochondrial and nuclear gene levels, we suggest that these two populations should be regarded as two evolutionary significant units, and targeted strategies should be urgently put forward to manage and protect its resources.
-
鳤 (Ochetobius elongatus) 隶属于鲤形目、鲤科、鳤属,曾是广泛分布于我国南方长江、珠江等诸多水系的重要经济鱼类[1-3]。受到涉水工程、江河环境污染、栖息地破坏、过度捕捞等因素的影响,其资源量萎缩极其严重,在很多河流江段已多年未被监测到[4-5]。2016年修订的《中国脊椎动物红色名录》已将鳤的保护等级上升至极度濒危[6];近期颁布的《长江保护法》还特别针对鳤提出了重要指示,因此亟需给予重视和保护。然而,由于鳤样本较难获得,目前其相关研究鲜见报道,零星研究仅集中在大阶元的系统发育进化关系[7-8]、小样本的遗传变异分析[5]和局域分布鳤种群的遗传多样性[9]方面,这极大制约了有效保护策略的制定。
珠江和长江是我国的两大重要水系,在鱼类生物多样性和种质资源孕育方面起着重要作用[3,10]。南岭山脉作为珠江和长江两大重要水系的分水岭,是阻隔两个水系鱼类迁移扩散的天然屏障。已有大量遗传进化研究表明,南岭山脉导致珠江和长江两个水系共同分布的鱼类发生了遗传分化甚至形成高度分化的遗传谱系[11-13]。在此背景下,长江鳤种群与珠江鳤种群是否形成显著的遗传分化,遗传分化的程度如何,若存在遗传分化是在什么时间节点开始的,驱动分化的环境因素是什么,这些关键问题均尚未找到答案。
利用遗传进化方法揭示物种的遗传结构对于物种的管理保护具有重要意义[14]。本文使用珠江和长江两个水系鳤群体的2个线粒体基因 [细胞色素b (Cytb) 和酰胺腺嘌呤二核苷酸 (NADH) 脱氧酶亚单位2 (ND2)] 和2个核基因 [肌球蛋白重链6 (myh6) 和重组激活基因2 (RAG2)] 的序列,利用系统发育分析、单倍型网状图、分化时间估算等方法分析了两个水系鳤种群的遗传结构,估算了分化时间,探讨了导致其分化的可能环境因素,以期为鳤的管理和保护提供科学依据。
1. 材料与方法
1.1 样品采集
长江的7尾鳤样本 (体长21.0~24.2 cm,体质量74.5~112.0 g) 由中国水产科学研究院长江水产研究所于2020年在长江湖北公安江段常规渔业调查采集获得 (图1和表1)。剪取少量胸鳍妥善保存于无水乙醇中,用以基因组DNA的提取。处理后的鳤样本带回中国水产科学研究院长江水产研究所,用甲醛固定制作成标本用以科研和科普。另外,在珠江广东肇庆江段采集鳡 (Elopichthys bambusa) 1尾,作为系统发育分析的外类群。
表 1 采样信息和GenBank序列号Table 1. Sampling information and GenBank No.采样站位
Sampling site样本量
Sample number/尾GenBank序列号 GenBank No. Cytb ND2 myh6 RAG2 公安 Gongan 7 MW657590—MW657596 MW657598—MW657604 MW657606—MW657612 MW657614—MW657620 梧州 Wuzhou 1 MW657432 MW657484 MW657536 MW657588 桂平 Guiping 1 MW657431 MW657483 MW657535 MW657587 肇庆 Zhaoqing 50 MW657381—MW657430 MW657433—MW657482 MW657485—MW657534 MW657537—MW657586 1.2 基因组DNA提取、扩增与测序
采用高盐法提取基因组DNA,并利用质量分数为1.2%琼脂糖凝胶电泳检测DNA质量。本文分别选取了2个线粒体基因 (Cytb和ND2) 和2个核基因 (myh6和RAG2)。PCR反应体积为30 μL,其中Taq PCR Master Mix (上海生工) 15 μL,包括Taq DNA 聚合酶、dNTP、PCR Buffer、PCR stabilizers、Gel loading和核酸染液,模板DNA 1 μL,正反向引物各1 μL (浓度为10 μmol·L−1),加灭菌超纯水至总体积30 μL。鳤和鳡4个基因扩增与测序的引物信息和PCR反应体系参考杨计平等[9]。PCR产物经质量分数为1.2%琼脂糖凝胶电泳检测,对于扩增效果良好的产物交由测序公司双向测序。
1.3 数据处理与分析
测序所得序列使用Lasergene软件包 (DNASTAR Inc., Madison, WI) 来检测测序质量并生成一致序列。所得的一致序列使用MUSCLE[14]进行比对,然后剪掉首尾两端噪音序列至同样长度。在GenBank核苷酸序列数据库下载了珠江52尾鳤样本已公布的4个基因序列,具体采集信息参考杨计平等[9] (图1和表1)。这些公布序列将与新获得序列比对后一起用于后续分析。为了获得更多的变异信息,分别合并了2个线粒体基因和2个核基因用以后续分析。线粒体DNA单倍型和核基因等位基因均使用DnaSP 6.12[15]进行划分。利用DnaSP 6.12统计单倍型数、计算单倍型多样性和核苷酸多样性,并估算了两个水系鳤群体的遗传分化状况。使用Network 10.2 (http://www.fluxus-engineering.com) 构建了Cytb+ND2的单倍型中介网状图[16]和myh6+RAG2的等位基因中介网状图[16]。
为了揭示珠江和长江鳤种群的系统发育关系,本文通过MrBayes 3.2.7[17]构建贝叶斯树,并以鳡作为外类群。合并了4个基因进行系统发育树的构建。贝叶斯树使用马尔科夫链蒙特卡洛方法 (Markov chain Monte Carlo, MCMC) 分析运行108万代,每1 000代取样1次,并舍弃初始的25×105代。另外,本文使用同样的参数设置对Cytb+ND2和myh6+RAG2也分别进行了贝叶斯树的构建。不同基因组合的最优核苷酸替代模型分别利用MrModeltest 2.3[18]计算获得 (表2)。最后,基于Kimura双参数模型 (Kimura-2-parameter, K2P)[19],利用MEGA 7.0[20]计算两个水系群体之间的遗传距离。
表 2 系统发育树构建和分化时间估算的最优模型Table 2. Optimal model of phylogenetic tree construction and differentiation time estimation基因组合
Gene combination最优模型
Optimal model不变位点比例
Proportion of invariable siteGamma分布形状参数
Gamma distribution shape parameterCytb+ND2 GTR+I+G 0.474 0.986 myh6+RAG2 K80+I 0.881 — Cytb+ND2+myh6+RAG2 GTR+I+G 0.697 0.857 本文联合4个基因,使用BEAST v.1.8.1[21]软件中的贝叶斯马尔可夫链蒙特卡洛方法估算两个水系鳤群体的分化时间,采用严格分子钟模型、Yule先验模型和GTR+I+G最优核苷酸替代模型。由于缺乏有效的化石记录,本文采用固定进化速率来估算分化时间。联合基因的进化速率由Cytb的进化速率估算获得。由于联合基因的K2P平均遗传距离 (0.29%) 与Cytb的K2P平均遗传距离 (0.43%) 的比值为0.674,因此利用Cytb的进化速率乘以比值便可获得联合基因的进化速率。1%~2%是鲤科鱼类Cytb基因普遍接受的一个进化速率范围[22-24],因此联合基因的进化速率为0.674%~1.348%。分别利用0.674%和1.348%的进化速率进行运算,每次运算共运行1×109代,每1 000代进行1次抽样,舍弃初始5×108代。运行结束后,利用Tracer 1.6[25]评价每个参数的稳定性状况。最后使用TreeAnnotator v.1.8.1[21]总结出一致树,在Figtree v.1.3.1 (http://tree.bio.ed.ac.uk/software/figtree) 显示节点的分化时间和后验概率。
2. 结果
2.1 序列信息
本文分析数据包括新获得的长江7尾鳤的全部4个基因片段以及杨计平等[9] 52尾珠江样本的对应4个基因片段,共计59尾样本。序列比对后共获得895 bp的Cytb序列,包含了30个变异位点;935 bp的ND2序列,包含了26个变异位点;763 bp的myh6序列,包含了3个变异位点;857 bp的RAG2序列,包含了5个变异位点。在4个基因片段中均未发现缺失和插入。此外,成功获得了鳡的4个基因序列,GenBank序列号分别为MW657589、MW657597、MW657605和MW657613。
2.2 遗传多样性、系统发育分析与遗传结构
基于线粒体联合基因和核基因联合基因估算的单倍型数目、单倍型多样性和核苷酸多样性见表3。基于线粒体基因估算的两个水系鳤的单倍型多样性为0.950,核苷酸多样性为0.494%。基于核基因计算的单倍型多样性和核苷酸多样性明显低于线粒体联合基因,这主要是与核基因较慢的进化速率有关。基于线粒体基因组合发现两个水系之间的遗传分化系数 ФST为0.863 (P<0.001),基于核基因组合获得两个水系之间的遗传分化系数FST为0.518 (P<0.001)。
表 3 鳤种群的单倍型与遗传多样性指数Table 3. Haplotypes and genetic diversity indexes of O. elongatus populations基因组合
Gene combination水系
River system序列数
Sequence number单倍型数
Haplotype number单倍型多样性
Haplotype diversity核苷酸多样性
Nucleotide diversity/%Cytb+ND2 珠江 52 26 0.940±0.017 0.219±0.020 长江 7 3 0.714±0.127 0.187±0.003 合计 59 29 0.950±0.014 0.494±0.080 myh6+RAG2 珠江 104 5 0.546±0.050 0.040±0.006 长江 14 8 0.923±0.044 0.135±0.021 合计 118 13 0.643±0.036 0.071±0.009 线粒体基因组合共界定了29个单倍型,其中珠江群体26个 (H1—H26),长江群体3个 (H27—H29,图2-a)。对于核基因组合,共鉴定了13个等位基因型 (A1—A13),其中珠江和长江群体分别界定了5和8个等位基因型 (图2-b)。基于线粒体基因组合的单倍型中介网状图表明,珠江和长江的鳤群体形成显著的遗传分支结构,并且两个分支之间相差22个突变位点 (图2-a)。基于核基因组合的等位基因中介网状图虽然未形成明显的遗传分支结构,但两个水系群体之间却各自形成了特有的等位基因型 (图2-b)。
基于4个基因组合和线粒体基因组合获得的系统发育树拓扑结构一致,均支持珠江和长江形成了显著的支系结构 (图3-a和图3-c)。基于核基因组合构建的贝叶斯树并未形成可信的支系结构 (图3-b)。基于4个基因组合、线粒体基因组合和核基因组合估算的两个水系群体间的遗传距离分别为0.87%、1.52%和0.14%。
2.3 分化时间估算
Beast运算的有效取样大小 (Effect sample size, ESS) 参数均大于200,表明运算是收敛的。基于不同进化速率估算的分化时间显示,珠江和长江的鳤群体形成了2个进化支系,并获得100%的后验概率,分化时间大概在 0.38~0.76 Ma (百万年前) (图4)。
3. 讨论
3.1 谱系分化
基于线粒体基因组合的贝叶斯系统发育树、单倍型网状图和基于4个联合基因的系统发育树均发现,两个水系的鳤群体存在高度分化的遗传谱系,且形成了严格的珠江-长江分化格局。此外,两个水系群体之间的遗传距离在线粒体基因层面上也呈现较高水平 (1.52%)。这与范凤娟和章群[5]的研究结果相类似,其获得了4尾长江样本和2尾珠江样本,利用线粒体Cytb基因发现两个水系样本各自形成了两个独立遗传谱系,且谱系间的遗传距离为1.5%。这些研究结果表明两个水系鳤群体经历了长期的地理隔离。作为珠江和长江水系的分水岭,南岭山脉可能是导致其分化的重要推动力。类似的例子也发生在珠江和长江共存的其他鲤科鱼类上,比如䱗 (Hemiculter leucisculus)[26]、马口鱼 (Opsariichthys bidens)[27]。基于核基因组合获得的系统发育树和等位基因网状图发现,两个水系鳤群体虽未形成独立的分支结构,却形成了特有的等位基因型,说明两个水系群体之间不存在基因交流现象。然而,由于本文中涉及的样本量较少,后期需要使用更多样本量来验证这个观点。
草鱼 (Ctenopharyngodon idella)、青鱼 (Mylophayngodon piceus)、鲢 (Hypophthalmichthys molitrix) 和鳙 (H. nobilis) 是我国重要的经济鱼类,在珠江和长江均有广泛分布[1-3]。基于草鱼[28]、青鱼[29]、鲢[30]和鳙[30]的相关遗传分化研究发现,虽然珠江和长江野生群体形成了一定程度的群体分化,但均未形成独立的遗传谱系。本文发现两个水系的鳤群体形成了两个深度分化谱系,表明鳤的分化程度显著高于两个水系共同分布的这几种重要经济鱼类。鳤有限的分布范围和较小的种群数量可能是造成这种差异遗传模式的原因[1,3]。早期的地质学和鱼类区系研究认为,珠江和长江水系的部分河流相隔距离很近,甚至有少数河流直接连通[31-33]。四大家鱼分布广泛、种群数量庞大,促使两个水系群体之间可以通过一些相邻的河流、江段进行迁移扩散,降低了群间的遗传分化。另外,四大家鱼人工繁育种群数量庞大,不可避免的养殖逃逸和人工增殖放流也会导致两个水系四大家鱼的遗传种质趋于同质。相比之下,鳤主要分布在大型河流江段与湖泊,种群规模极其有限,加上人为干预极少,因此两个水系群体之间进行相互扩散的概率较小。
已有大量研究表明古地质历史事件是促进物种分化的重要驱动因子[26-27,34-35]。本文发现珠江和长江鳤群体的谱系分化大概发生在更新世中期。这一时期刚好与上新世后期开始的青藏高原快速隆升阶段第二幕的昆仑-黄河运动时间段吻合 (1.1~0.6 Ma)。研究表明,青藏高原在上新世后期进行快速抬升,共分为3个阶段[36-37]:第一阶段是青藏运动,时间介于3.6~1.7 Ma;第二阶段是昆仑-黄河运动,时间介于1.1~0.6 Ma;第三阶段是共和运动,时间大致在0.15 Ma。青藏高原的快速隆升重塑了东亚地区的景观特征,包括河流和山脉。早期研究也证明了青藏高原的快速抬升导致了一些鱼类的分化甚至成种[26,38-40]。因此,两个水系鳤的谱系分化很可能与青藏高原的快速隆升阶段相关。
3.2 保护建议
针对特定物种展开遗传结构解析能够有效辅助物种的管理和保护[26,41-42]。鳤的资源极其缺乏,近几年关于鳤的零星发现都能引起社会的广泛关注。由于鳤的野外资源稀少以及人工繁殖尚未成功,细化保护单元并展开就地保护是一个重要举措。因此,通过遗传进化分析解析鳤在不同水系的遗传结构对于细化其保护单元非常必要。针对濒危或稀有物种的保育,许多学者认为演化显著单位 (Evolutionary Significant Unit, ESU) 为适当的保护管理单元[43]。ESU是指一个种群在线粒体基因或者核基因水平已经形成独立的单系群 (Monophyly)[44]。本文发现珠江和长江鳤群体形成了独立的线粒体基因遗传谱系和特有的等位基因型,因此可以看作两个不同的ESU,在保护策略制定和实施上应该区别对待。然而,由于本研究样本量和分子标记有限,今后需要通过扩大研究范围和样本量,结合更多的标记 (如核基因组数据) 来全面评估鳤的种群遗传结构,以期为鳤的保育策略的制定提供更系统的支撑。
-
表 1 采样信息和GenBank序列号
Table 1 Sampling information and GenBank No.
采样站位
Sampling site样本量
Sample number/尾GenBank序列号 GenBank No. Cytb ND2 myh6 RAG2 公安 Gongan 7 MW657590—MW657596 MW657598—MW657604 MW657606—MW657612 MW657614—MW657620 梧州 Wuzhou 1 MW657432 MW657484 MW657536 MW657588 桂平 Guiping 1 MW657431 MW657483 MW657535 MW657587 肇庆 Zhaoqing 50 MW657381—MW657430 MW657433—MW657482 MW657485—MW657534 MW657537—MW657586 表 2 系统发育树构建和分化时间估算的最优模型
Table 2 Optimal model of phylogenetic tree construction and differentiation time estimation
基因组合
Gene combination最优模型
Optimal model不变位点比例
Proportion of invariable siteGamma分布形状参数
Gamma distribution shape parameterCytb+ND2 GTR+I+G 0.474 0.986 myh6+RAG2 K80+I 0.881 — Cytb+ND2+myh6+RAG2 GTR+I+G 0.697 0.857 表 3 鳤种群的单倍型与遗传多样性指数
Table 3 Haplotypes and genetic diversity indexes of O. elongatus populations
基因组合
Gene combination水系
River system序列数
Sequence number单倍型数
Haplotype number单倍型多样性
Haplotype diversity核苷酸多样性
Nucleotide diversity/%Cytb+ND2 珠江 52 26 0.940±0.017 0.219±0.020 长江 7 3 0.714±0.127 0.187±0.003 合计 59 29 0.950±0.014 0.494±0.080 myh6+RAG2 珠江 104 5 0.546±0.050 0.040±0.006 长江 14 8 0.923±0.044 0.135±0.021 合计 118 13 0.643±0.036 0.071±0.009 -
[1] 乐佩琦. 中国动物志. 硬骨鱼纲. 鲤形目 (中卷)[M]. 北京: 科学出版社, 2000: 106-108. [2] 周解, 张春光. 广西淡水鱼类志[M]. 南宁: 广西人民出版社, 2006: 140-142. [3] 郑慈英. 珠江鱼类志[M]. 北京: 科学出版社, 1989: 81-82. [4] 李捷, 李新辉, 贾晓平, 等. 西江鱼类群落多样性及其演变[J]. 中国水产科学, 2010, 17(2): 298-311. [5] 范凤娟, 章群. 汩罗江与柳江鳤鱼线粒体细胞色素b基因的遗传变异初探[J]. 生态科学, 2009, 28(5): 457-459. doi: 10.3969/j.issn.1008-8873.2009.05.015 [6] 蒋志刚, 江建平, 王跃招, 等. 中国脊椎动物红色名录[J]. 生物多样性, 2016, 24(5): 500-551. doi: 10.17520/biods.2016076 [7] WANG X, LI J, HE S. Molecular evidence for the monophyly of East Asian groups of Cyprinidae (Teleostei: Cypriniformes) derived from the nuclear recombination activating gene 2 sequences[J]. Mol Phylogenet Evol, 2007, 42(1): 157-170. doi: 10.1016/j.ympev.2006.06.014
[8] TAO W, ZOU M, WANG X, et al. Phylogenomic analysis resolves the formerly intractable adaptive diversification of the endemic clade of east Asian Cyprinidae (Cypriniformes)[J]. PLoS One, 2010, 5(10): e13508. doi: 10.1371/journal.pone.0013508
[9] 杨计平, 李策, 陈蔚涛, 等. 西江中下游鳤的遗传多样性与种群动态历史[J]. 生物多样性, 2018, 26(12): 1289-1295. doi: 10.17520/biods.2018121 [10] XING Y, ZHANG C, FAN E, et al. Freshwater fishes of China: species richness, endemism, threatened species and conservation[J]. Divers Distrib, 2016, 22(3): 358-370. doi: 10.1111/ddi.12399
[11] YANG L, MAYDEN R L, HE S. Population genetic structure and geographical differentiation of the Chinese catfish Hemibagrus macropterus (Siluriformes, Bagridae): evidence for altered drainage patterns[J]. Mol Phylogenet Evol, 2009, 51(2): 405-411. doi: 10.1016/j.ympev.2009.01.004
[12] PERDICES A, CUNHA C, COELHO M M. Phylogenetic structure of Zacco platypus (Teleostei, Cyprinidae) populations on the upper and middle Chang Jiang (=Yangtze) drainage inferred from cytochrome b sequences[J]. Mol Phylogenet Evol, 2004, 31(1): 192-203. doi: 10.1016/j.ympev.2003.07.001
[13] PERDICES A, SAYANDA D, COELHO M M. Mitochondrial diversity of Opsariichthys bidens (Teleostei, Cyprinidae) in three Chinese drainages[J]. Mol Phylogenet Evol, 2005, 37(3): 920-927. doi: 10.1016/j.ympev.2005.04.020
[14] EDGAR R C. MUSCLE: multiple sequence alignment with high accuracy and high throughput[J]. Nucleic Acids Res, 2004, 32(5): 1792-1797. doi: 10.1093/nar/gkh340
[15] ROZAS J, FERRER-MATA A, SÁNCHEZ-DELBARRIO J, et al. DnaSP 6: DNA sequence polymorphism analysis of large data sets[J]. Mol Biol Evol, 2017, 34(12): 3299-3302. doi: 10.1093/molbev/msx248
[16] BANDELT H J, FORSTER P, ROHL A. Median-joining networks for inferring intraspecific phylogenies[J]. Mol Biol Evol, 1999, 16(1): 37-48. doi: 10.1093/oxfordjournals.molbev.a026036
[17] RONQUIST F, TESLENKO M, van der MARK P, et al. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space[J]. Syst Biol, 2012, 61(3): 539-542. doi: 10.1093/sysbio/sys029
[18] NYLANDER J. MrModeltest v2. Program distributed by the author (ed. Nylander JAA)[M]. Uppsala: Evolutionary Biology Centre, Uppsala University, 2004: 1-2.
[19] KIMURA M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences[J]. J Mol Evol, 1980, 16(2): 111-120. doi: 10.1007/BF01731581
[20] KUMAR S, STECHER G, TAMURA K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets[J]. Mol Biol Evol, 2016, 33(7): 1870-1874. doi: 10.1093/molbev/msw054
[21] DRUMMOND A J, RAMBAUT A. BEAST: Bayesian evolutionary analysis by sampling trees[J]. BMC Evol Biol, 2007, 7: 214. doi: 10.1186/1471-2148-7-214
[22] DURAND J D, TSIGENOPOULOS C S, UNLU E, et al. Phylogeny and biogeography of the family Cyprinidae in the Middle East inferred from cytochrome b DNA: evolutionary significance of this region[J]. Mol Phylogenet Evol, 2002, 25(1): 218-218. doi: 10.1016/S1055-7903(02)00343-3
[23] KETMAIER V, BIANCO P G, COBOLLIA M et al. Molecular phylogeny of two lineages of Leuciscinae cyprinids (Telestes and Scardinius) from the peri-Mediterranean area based on cytochrome b data[J]. Mol Phylogenet Evol, 2004, 32(3): 1061-1071. doi: 10.1016/j.ympev.2004.04.008
[24] MEYER A. Evolution of mitochondrial DNA in fishes[M]. The Hague: Elsevier, 1993: 1-38.
[25] RAMBAUT A, DRUMMOND A J. Tracer v1. 4[CP/OL]. [2021-12-20]. http://beast.bio.ed.ac.uk/Tracer.
[26] CHEN W, ZHONG Z, DAI W, et al. Phylogeographic structure, cryptic speciation and demographic history of the sharpbelly (Hemiculter leucisculus), a freshwater habitat generalist from southern China[J]. BMC Evol Biol, 2017, 17(1): 216. doi: 10.1186/s12862-017-1058-0
[27] LI G Y, WANG X Z, ZHAO Y H, et al. Speciation and phylogeography of Opsariichthys bidens (Pisces: Cypriniformes: Cyprinidae) in China: analysis of the cytochrome b gene of mtDNA from diverse populations[J]. Zool Stud, 2009, 48(4): 569-583.
[28] 朱叶. 基于线粒体细胞色素b的3江水系草鱼群体遗传多样性分析[D]. 广州: 暨南大学, 2012: 1-54. [29] 付晓艳. 长江和珠江水系青鱼线粒体细胞色素b基因遗传多样性分析[D]. 广州: 暨南大学, 2011: 1-42. [30] LI C, WANG J, CHEN J, et al. Native bighead carp Hypophthalmichthys nobilis and silver carp Hypophthalmichthys molitrix populations in the Pearl River are threatened by Yangtze River introductions as revealed by mitochondrial DNA[J]. J Fish Biol, 2020, 96(3): 651-662. doi: 10.1111/jfb.14253
[31] 任美锷, 包浩生, 韩同春, 等. 云南西北部金沙江河谷地貌与河流袭夺问题[J]. 地理学报, 1959, 26(2): 135-155. doi: 10.11821/xb195902003 [32] LIU H Z. A preliminary analysis to biogeographical process of the eastern Asian freshwater fishes[J]. Acta Zootaxonomica Sinica, 1998, 23(Suppl): 49-55.
[33] YAP S Y. On the Distributional patterns of southeast-east Asian freshwater fish and their history[J]. J Biogeogr, 2002, 29(9): 1187-1199. doi: 10.1046/j.1365-2699.2002.00771.x
[34] DUFRESNES C, NICIEZA A G, LITVINCHUK S N, et al. Are glacial refugia hotspots of speciation and cytonuclear discordances? Answers from the genomic phylogeography of Spanish common frogs[J]. Mol Ecol, 2020, 29(5): 986-1000. doi: 10.1111/mec.15368
[35] RINCÓN-SANDOVAL M, BETANCUR R R, OCAMPO J A M. Comparative phylogeography of trans-Andean freshwater fishes based on genome-wide nuclear and mitochondrial markers[J]. Mol Ecol, 2019, 28(5): 1096-1115. doi: 10.1111/mec.15036
[36] 李吉均, 方小敏. 青藏高原隆起与环境变化研究[J]. 科学通报, 1998, 43(15): 1568-1574. [37] 李吉均, 方小敏, 潘保田, 等. 新生代晚期青藏高原强烈隆起及其对周边环境的影响[J]. 第四纪研究, 2001, 21(5): 381-391. doi: 10.3321/j.issn:1001-7410.2001.05.001 [38] HE D K, CHEN Y F. Biogeography and molecular phylogeny of the genus Schizothorax (Teleostei: Cyprinidae) in China inferred from cytochrome b sequences[J]. J Biogeogr, 2006, 33(8): 1448-1460. doi: 10.1111/j.1365-2699.2006.01510.x
[39] YANG J, YANG J X, CHEN X Y. A re-examination of the molecular phylogeny and biogeography of the genus Schizothorax (Teleostei: Cyprinidae) through enhanced sampling, with emphasis on the species in the Yunnan-Guizhou Plateau, China[J]. J Zool Syst Evol Res, 2012, 50(3): 184-191. doi: 10.1111/j.1439-0469.2012.00661.x
[40] 曹文宣, 陈宜瑜, 武云飞, 等. 裂腹鱼类的起源和演化及其与青藏高原隆起的关系[M]. 北京: 科学出版社, 1981: 118-130. [41] 李敏, 黄梓荣, 许友伟, 等. 基于线粒体cytb序列的花斑蛇鲻种群遗传结构研究[J]. 南方水产科学, 2019, 15(6): 41-48. doi: 10.12131/20190123 [42] 徐豪, 梁绪虹, 王丛丛, 等. 基于线粒体NADH脱氢酶亚基2标记的东南太平洋不同表型间茎柔鱼群体遗传学分析[J]. 南方水产科学, 2022, 18(1): 153-159. doi: 10.12131/20210119 [43] RYDER O A. Species conservation and systematics: the dilemma of subspecies[J]. Trends Ecol Evol, 1986, 1: 9-10. doi: 10.1016/0169-5347(86)90059-5
[44] MORITZ C. Defining 'Evolutionarily Significant Units' for conservation[J]. Trends Ecol Evol, 1994, 9(10): 373-375. doi: 10.1016/0169-5347(94)90057-4
-
期刊类型引用(1)
1. 翟东东,王东,何天琪,蔡方陶,段辛斌,胡兴坤,朱滨,谢仲桂,夏明,刘红艳,熊飞. 长江中游鳤的遗传多样性及种群历史动态分析. 水生生物学报. 2024(10): 1736-1744 . 百度学术
其他类型引用(1)