基于基因组重测序的长江上游瓦氏黄颡鱼群体遗传结构_熊飞.pdf

2023,31(4):22391,pages 110 doi:10.17520/biods.2022391 收稿日期:2022-07-09;接受日期:2022-11-14 基金项目:国家自然科学基金(51979123)、湖北省高等学校优秀中青年科技创新团队计划项目(T2020034)和江汉大学科技创新专项项目(2021kjzx006)*共同通讯作者 Co-authors for correspondence.E-mail:; https:/www.biodiversity-基于基因组重测序的长江上游瓦氏黄颡鱼群体遗传结构 熊飞1,2,刘红艳1,2*,翟东东1,2,段辛斌3*,田辉伍3,陈大庆3 1.江汉大学生命科学学院,湖北省汉江流域特色生物资源保护开发与利用工程技术研究中心,武汉 430056;2.江汉大学持久性有毒污染物环境与健康危害湖北省重点实验室,武汉 430056;3.中国水产科学研究院长江水产研究所,武汉 430223 摘要:利用基因组重测序的方法获取高通量SNP标记,分析了长江上游三峡大坝-白鹤滩大坝之间8个不同江段(太平溪、巴南、合川、岷江口、宜宾、邵女坪、桧溪、冯家坪)共136尾瓦氏黄颡鱼(Pelteobagrus vachelli)的遗传多样性和遗传分化水平,阐明了长江上游瓦氏黄颡鱼群体遗传结构。结果显示:(1)三峡库区太平溪群体和巴南群体具有较高的SNP(single nucleotide polymorphism)数量和核苷酸多样性指数,遗传来源丰富,其遗传多样性高于其他群体;上游的岷江口、宜宾、邵女坪和冯家坪群体遗传来源单一。(2)瓦氏黄颡鱼存在3个不同的遗传分支,且不同遗传分支之间存在较大的遗传分化。(3)群体SNP数量和核苷酸多样性指数与河流坡降呈显著负相关,群体遗传分化指数与地理距离和隔离时间无显著相关性。研究结果表明,在三峡大坝-白鹤滩大坝江段,瓦氏黄颡鱼上游群体具有更低的遗传多样性,更易发生遗传漂变作用,在鱼类遗传多样性保护中需要特别关注;瓦氏黄颡鱼存在3种显著的遗传结构,应视为3个不同遗传单元进行种质资源管理。关键词:瓦氏黄颡鱼;长江上游;遗传分化;遗传多样性;基因组重测序 熊飞,刘红艳,翟东东,段辛斌,田辉伍,陈大庆(2023)基于基因组重测序的长江上游瓦氏黄颡鱼群体遗传结构.生物多样性,31,22391.doi:10.17520/biods.2022391.Xiong F,Liu HY,Zhai DD,Duan XB,Tian HW,Chen DQ(2023)Population genetic structure of Pelteobagrus vachelli in the upper Yangtze River based on genome re-sequencing.Biodiversity Science,31,22391.doi:10.17520/biods.2022391.Population genetic structure of Pelteobagrus vachelli in the upper Yangtze River based on genome re-sequencing Fei Xiong1,2,Hongyan Liu1,2*,Dongdong Zhai1,2,Xinbin Duan3*,Huiwu Tian3,Daqing Chen3 1 School of Life Sciences,Jianghan University,Hubei Engineering Research Center for Protection and Utilization of Special Biological Resources in the Hanjiang River Basin,Wuhan 430056 2 Hubei Key Laboratory of Environmental and Health Effects of Persistent Toxic Substances,Jianghan University,Wuhan 430056 3 Yangtze River Fisheries Research Institute,Chinese Academy of Fishery Sciences,Wuhan 430223 ABSTRACT Aims:Pelteobagrus vachelli is an important economic fish in the Yangtze River.Nevertheless,it shows a trend of younger age and miniaturization,and the number of reproductive population above the third age was small,which is not benefit to the long-term maintenance of its resources.Therefore,more attention and conservation measure are needed for this species.The aim of this study is to clarify genetic diversity and genetic structure of P.vachelli.Methods:We collected 136 samples of eight different river sections(Taipingxi,Banan,Hechuan,Minjiangkou,Yibin,Shaonping,Huixi,Fengjiaping)from the Three Gorges Dam to Baihetan Dam in the upper Yangtze River.High-throughput SNP(single nucleotide polymorphism)markers through genome re-sequencing method was used to study genetic diversity and genetic differentiation of P.vachelli populations.Results:The number of SNP and nucleotide diversity of Taipingxi and Banan populations in the Three Gorges Reservoir region were significantly higher than those in the other populations,and there were three genetic branches 研究报告 熊飞等:基于基因组重测序的长江上游瓦氏黄颡鱼群体遗传结构研究 2023年|31卷|4期|22391|第2页 and abundant genetic sources in the two populations.However,there was only one genetic branch and single genetic source in Minjiangkou,Yibin,Shaonping and Fengjiaping populations.The populations could be divided into three distinct genetic groups with great genetic differentiation.All individuals in Fengjiaping and Minjiangkou belonged to group 1,all individuals in Shaonping and Yibin belonged to group 2,Huixi and Hechuan populations had two genetic branches(group 1 and group 2),and Banan and Taipingxi populations had three genetic branches(group 1,group 2 and group 3).The SNP number and nucleotide diversity were significantly affected by channel slope.The genetic differentiation was not related to geographical distance and isolated time.Conclusion:The upstream populations with lower genetic diversity are more prone to genetic drift and thus need more attention in the fish protection.There are three distinct genetic structures in the upper Yangtze River,which can be regarded as three different genetic units for management of germplasm resource.Key words:Pelteobagrus vachelli;the upper Yangtze River;genetic differentiation;genetic diversity;genome re-sequencing 瓦氏黄颡鱼(Pelteobagrus vachelli)隶属于鲿科 黄颡鱼属,广泛分布于长江及其支流,是长江常见的底栖定居性鱼类(丁瑞华,1994)。过度捕捞等人为因素导致其资源衰退,瓦氏黄颡鱼被认为是中度超额开发物种(陈大庆等,2002)。另外,水电工程的建设导致许多河段水位抬高,淹没大量砾石和水草等附着物,影响粘性卵的附着,造成瓦氏黄颡鱼等产粘性卵鱼类的产卵场消失(姚维志,2018)。虽然渔业资源调查显示,瓦氏黄颡鱼目前在长江渔获物中的占比仍处于优势地位,但是其群体呈现出低龄化和小型化的现象,渔获物中大部分个体未达到性成熟,3龄以上繁殖群体较少,不利于其资源的可持续利用(杨家云,1994;熊飞等,2015)。遗传多样性是物种进化和适应的基础,维持群体遗传多样性有利于物种长期的进化和生存(Nakajima et al,2020)。瓦氏黄颡鱼是长江重要的经济鱼类,探明其群体遗传多样性和遗传结构,对于评估其种质状况和制定合理的管理政策至关重要。长江上游是水电大坝梯级开发集中、河流生境破碎化严重的区域,已建成三峡、向家坝、溪洛渡等多座大型水电站(Cheng et al,2015)。关于长江上游鱼类的相关研究主要集中在特有鱼类和洄游性鱼类等方面(Song et al,2016;Liu et al,2017),而对于广布种及定居性鱼类的关注较少。关于瓦氏黄颡鱼群体遗传结构的研究,涉及长江上游区域的较少,且所用的分子标记多为线粒体DNA标记(Wang et al,2004),获得的遗传信息有限。单核苷酸多态性(single nucleotide polymorphism,SNP)是基因组变异的主要方式之一,其遗传稳定性好、通量高。重测序技术将新获得的序列与已知的参考基因组序列进行比对,获取大量的SNP位点,为群体遗传研究提供了重要手段(Kang et al,2017)。本研究利用基因组重测序的方法获取高通量SNP标记,分析了长江上游三峡大坝-白鹤滩大坝之间8个不同江段(太平溪、巴南、合川、岷江口、宜宾、邵女坪、桧溪、冯家坪)瓦氏黄颡鱼的遗传多样性和遗传分化水平,探讨了长江上游瓦氏黄颡鱼群体遗传结构,以期为长江上游鱼类遗传种质资源保护和管理提供依据。1.1 研究区域和样本采集 2019年8月10月于长江上游采集瓦氏黄颡鱼样品,采集区域为三峡大坝至白鹤滩大坝约1,300 km的江段。共设置8个采样点,分别为太平溪(TP)、巴南(BN)、合川(HC)、岷江口(MJ)、宜宾(YB)、邵女坪(SN)、桧溪(HX)和冯家坪(FJ)(图1,表1)。样点太平溪、巴南和宜宾位于三峡大坝(2003年截流)和向家大坝(2008年截流)间,该河段长度约1,000 km,生境相对完好,存在从上游河流生境到库区生境的变化。其中,太平溪和巴南分别位于三峡水库的库首和库尾。样点邵女坪和桧溪位于向家大坝和溪洛渡大坝(2007年截流)之间。样点冯家坪位于溪洛渡大坝和白鹤滩大坝(2015年截流)之间。样点合川位于支流嘉陵江,与长江干流有草街水电站(2005年截流)阻隔。样点岷江口处于岷江支流的下游,与长江干流相隔约18 km的距离,其上游有犍为水电站(2019年截流)。取瓦氏黄颡鱼的尾鳍和背部肌肉组织,用95%乙醇固定,带回实验室保存备用。每个样点能获取的瓦氏黄颡鱼种群数量不同,当样点获取的样 1 材料与方法 熊飞等:基于基因组重测序的长江上游瓦氏黄颡鱼群体遗传结构研究 2023年|31卷|4期|22391|第3页 图1 长江上游瓦氏黄颡鱼的样点分布。S1:太平溪(TP);S2:巴南(BN);S3:合川(HC);S4:岷江口(MJ);S5:宜宾(YB);S6:邵女坪(SN);S7:桧溪(HX);S8:冯家坪(FJ)。Fig.1 Sampling sites of Pelteobagrus vachelli in the upper reaches of Yangtze River.S1,Taipingxi(TP);S2,Banan(BN);S3,Hechuan(HC);S4,Minjiangkou(MJ);S5,Yibin(YB);S6,Shaonping(SN);S7,Huixi(HX);S8,Fengjiaping(FJ).表1 瓦氏黄颡鱼8个群体的样点信息。样点缩写见图1。Table 1 Sampling information in eight populations of Pelteobagrus vachelli.Site abbreviations are shown in Fig.1.群体 Population 采样时间 Sampling time 样本数 No.of sample海拔 Elevation(m)坡降 Channel slope(m/km)生境类型 Habitat type 太平溪 TP 2019.8 14 110 0.12 水库 Reservoir 巴南 BN 2019.8 25 160 0.20 回水区 Backwater area 合川 HC 2019.9 4 187 1.85 水库 Reservoir 岷江口 MJ 2019.10 25 257 2.26 河流 River 宜宾 YB 2019.10 6 253 0.70 河流 River 邵女坪 SN 2019.10 12 296 0.90 水库 Reservoir 桧溪 HX 2019.10 25 356 1.40 河流 River 冯家坪 FJ 2019.10 25 578 1.50 河流 River 本量较大时,随机选取25个样本进行研究,当样点获取的样本量小于25个时,使用全部的样本进行研究(表1)。1.2 基因组DNA提取和重测序文库的构建 瓦氏黄颡鱼基因组DNA采用Foregene动物组织基因组试剂盒进行提取,获取的基因组DNA用TE缓冲液溶解。用1%的琼脂糖凝胶电泳检测基因组DNA样品是否有降解和杂质,用NanoDrop one超微量分光光度计(Thermo Scientific)检测基因组DNA纯度和浓度。取约1 g基因组DNA模板,根据TruSeq DNA Sample Preparation Guide(Illumina,15026486 Rev.C)的方法及流程进行重测序文库制备。具体如下:用超声破碎的方法将DNA随机打断,挑选350 bp的片段进行DNA片段的末端修复和3端加A反应,然后加上测序接头,通过PCR反应富集目的片段,形成测序文库。检测合格的文库在安诺优达基因科技有限公司的HiSeq X Ten测序平台上进行双端测序(PE150),得到150 bp的raw reads。对raw reads进行过滤,去除带接头的序列,过滤N含量超过10%的序列,去除质量值低于10的碱基数超过整条序列50%的reads,得到clean reads。1.3 序列比对及SNP变异检测 利用Bowtie2软件(Langmead et al,2012)将得到的瓦氏黄颡鱼clean reads与近缘种黄颡鱼参考基因组(Gong et al,2018)进行比对,把序列定位在参考熊飞等:基于基因组重测序的长江上游瓦氏黄颡鱼群体遗传结构研究 2023年|31卷|4期|22391|第4页 基因组上相应位置,统计测序深度、基因组覆盖度等信息。利用GATK(McKenna et al,2010)和Samtools(Li et al,2009)检测SNP位点。利用Vcftools软件(Danecek et al,2011)进行SNP位点过滤,过滤标准为:质量值(quality value)30,次等位基因频率(minor allele frequency,MAF)5%,哈迪-温伯格平衡检验(Hardy-Weinberg equilibrium test,HWE)0.01,所有个体中每个位点缺失基因型的个体不超过4个。1.4 遗传多样性和遗传结构分析 利用Vcftools软件(Danecek et al,2011)计算每个群体的SNP数量、核苷酸多样性指数()和遗传分化指数(Fst)。瓦氏黄颡鱼的群体遗传结构利用Admixture软件(Alexander et al,2009)进行分析,K值设置为110,对聚类结果进行交叉验证,采用交叉验证错误率最低值判断最优聚类数。在RAxML v8.1.17软件中(Stamatakis,2014)基于GTR+I+G模型 的 Bayesian 方 法 构 建 最 大 似 然 树(maximum likelihood,ML),设置100,000次重复。使用Eigen软件(Price et al,2006)进行主成分分析(principal component analysis,PCA),使用Plink软件(Purcell et al,2007)计算IBS(identity by state,IBS)值来评估群体的遗传相似度。1.5 环境因子对遗传结构的影响 环境因子主要考虑海拔高度、河流坡降、地理距离和隔离时间等因素。用GPS测量每个样点的海拔高度(精度1 m);河道坡降即单位河段长度的落差,按样点上游5 km和下游5 km的河段范围计算;地理距离(沿着河道测量)利用ArcGIS 10.5的大比例尺流域地形图进行测量;隔离时间从大坝截流时间算起。用R软件里的cor.test函数分析遗传多样性指数(SNP数量、核苷酸多样性指数)与河道坡降和海拔的相关性。用R软件vegan程序(Dixon,2003)进行Mantel检验,分析遗传分化指数与地理距离和隔离时 间 的 关 系,遗 传 分 化 指 数 采 用 线 性 化 的Fst/(1Fst)遗传分化指数。2.1 高通量测序数据分析 共获得了瓦氏黄颡鱼136个样本的基因组重测序序列,序列已上传至GenBank(Bioproject accession number为PRJNA803961,Biosample accession numbers为SAMN25692057SAMN25692192)。各群体的测序参数统计如表2,共得到了1,345.21 Gb的clean data,平均测序质量值Q30为92.78%,平均测序深度为8.95,平均测序比对率为50.73%,平均基因组覆盖率为69.43%。2.2 遗传多样性 通过与黄颡鱼参考基因组比对,共获得了瓦氏黄颡鱼7,341,959个SNP位点,SNP数量从419,733(岷江口)到4,172,876(太平溪),群体最大的SNP数量约为最小SNP数量的10倍,显示了SNP数量较大的变化范围。在三峡大坝-白鹤滩大坝江段,瓦氏黄颡鱼上游群体的SNP数量和核苷酸多样性指数低于下游的太平溪和巴南群体(表2)。2.3 群体遗传结构 当K=3时,交叉验证错误率最小(图2a),显示瓦氏黄颡鱼群体可分为3个不同遗传分支(图2b)。冯 表2 瓦氏黄颡鱼8个群体重测序数据及群体遗传多样性参数 Table 2 The genome re-sequencing data and genetic diversity parameters in eight populations of Pelteobagrus vachelli 群体 Population 样本数 No.of sample Clean data(Gb)Q30 (%)测序深度Sequencing depth 序列比对率 Sequence mapped(%)基因组覆盖度Coverage of genomic reference(%)SNP数目 Number of SNP 核苷酸多样性指数Nucleotide diversity 太平溪 TP 14 145.22 93.11 9.54 54.31 72.92 4,172,876 2.44 103 巴南 BN 25 236.41 92.82 8.71 53.43 71.68 3,543,176 2.81 103 合川 HC 4 35.77 92.81 8.38 49.94 66.28 534,504 8.59 104 岷江口 MJ 25 217.32 92.73 8.25 50.68 66.39 419,733 5.60 104 宜宾 YB 6 50.16 91.87 7.44 48.50 67.36 563,968 6.76 104 邵女坪 SN 12 128.29 92.79 9.37 49.09 69.64 457,603 5.75 104 桧溪 HX 25 258.76 92.79 9.20 49.22 68.87 597,465 6.60 104 冯家坪 FJ 25 273.27 92.78 9.57 49.06 69.73 441,092 5.84 104 总体 Total 136 1,345.21 92.78 8.95 50.73 69.43 7,341,959 1.19 103 2 结果 熊飞等:基于基因组重测序的长江上游瓦氏黄颡鱼群体遗传结构研究 2023年|31卷|4期|22391|第5页 图2 基于SNP标记分析8个瓦氏黄颡鱼群体的遗传结构。(a)由交叉验证错误率确定的最优分组数;(b)K=3时,所有个体的分组情况;(c)PCA聚类情况;(d)ML聚类树。样点缩写见图1。Fig.2 Genetic structure based on SNP markers in eight populations of Pelteobagrus vaclerii.(a)The optimal number of clusters determined by cross-validated error rate valley.(b)Population structure(K=3)of all individuals.(c)PCA plot of every individual.(d)ML phylogenetic tree.Site abbreviations are shown in Fig.1.家坪和岷江口的所有个体属于group 1,邵女坪和宜宾所有个体属于group 2,冯家坪、岷江口、邵女坪和宜宾群体的遗传背景较为单一;桧溪和合川的群体有2种遗传背景,桧溪绝大多数个体属于group 2,熊飞等:基于基因组重测序的长江上游瓦氏黄颡鱼群体遗传结构研究 2023年|31卷|4期|22391|第6页 少量个体属于group 1;合川群体里group 1和group 2遗传来源的个体数量大体相当;巴南和太平溪有3种遗传背景的个体,显示这两个地理位置的瓦氏黄颡鱼遗传来源较为丰富,其中group 1所占的比例最大,group 3所占的比例最小(图2b)。PCA分析(图2c)和ML聚类树(图2d)也进一步验证了这种结果。2.4 群体遗传分化 冯家坪与岷江口遗传分化指数较小(0.016),但这两个群体与其他群体间的遗传分化指数较大(0.0750.261)。桧溪、邵女坪和宜宾之间的遗传分化指数为负值或值较小(0.008至0.002),这3个群体的遗传背景非常相似,但桧溪和邵女坪群体与其他群体遗传分化指数较大(0.068至0.257)。合川与巴南和太平溪之间也为负值或者值较小(0.109至0.050),另外宜宾群体与巴南和太平溪之间的遗传分化指数也较小(0.020至0.050)(表3),结合structure结构图来看,巴南和太平溪的一部分遗传背景与合川相似,一部分遗传背景与宜宾相似。另外,我们比较3种遗传分支之间的遗传分化指数(structure图中的group 1、group 2、group 3),group 1与group 2间的遗传分化指数为0.260,group 2与group 3间的遗传分化指数为0.485,group 1与group 3间的遗传分化指数为0.510,可以看出,瓦氏黄颡鱼3种遗传分支之间存在较大的遗传分化。表3 长江上游瓦氏黄颡鱼两两群体之间的遗传分化指数 Table 3 Genetic differentiation between the pairwise populations of Pelteobagrus vaclerii in the upper Yangtze River 群体 Population 太平溪 TP 巴南 BN 合川 HC 岷江口 MJ 宜宾 YB 邵女坪SN 桧溪HX 巴南 BN 0.050 合川 HC 0.109 0.039 岷江口 MJ 0.112 0.143 0.075 宜宾 YB 0.020 0.043 0.152 0.233 邵女坪 SN 0.068 0.112 0.192 0.235 0.000 桧溪 HX 0.126 0.154 0.079 0.179 0.008 0.002 冯家坪 FJ 0.125 0.152 0.092 0.016 0.261 0.2570.193 IBS结果显示,3个遗传分支内个体间的遗传相似度值均较高,其中group 2的遗传相似度最高,group 1其次,group 3遗传相似度最低(图3)。2.5 环境因素对遗传结构的影响 相关性分析显示,SNP数量和核苷酸多样性指 图3 瓦氏黄颡鱼3个遗传分支内个体间的遗传相似度IBS值 Fig.3 Genetic similarity within three genetic groups based on identity by state(IBS)图4 瓦氏黄颡鱼SNP数量与河道海拔(a)和坡降(b)的相关性分析。P 0.05表示显著相关。Fig.4 Correlation analysis of the SNP numbers with elevation(a)and channel slope(b).P 0.05 indicates significant correlation.熊飞等:基于基因组重测序的长江上游瓦氏黄颡鱼群体遗传结构研究 2023年|31卷|4期|22391|第7页 数与海拔没有显著相关关系(图4a,图5a);SNP数量与河道坡降显著负相关(r=0.777,P=0.023)(图4b),核苷酸多样性指数也与河道坡降显著负相关(r=0.751,P=0.032)(图5b),说明瓦氏黄颡鱼的遗传多样性受河流坡降的影响。Mantel 检验显示,遗传分化指数与地理距离和隔离时间相关性不显著(图6),表明当前瓦氏黄颡鱼的群体遗传分化受地理距离和隔离时间的影响不明显。3.1 瓦氏黄颡鱼遗传多样性与遗传结构 在三峡大坝-白鹤滩大坝江段,瓦氏黄颡鱼上游群体SNP数量和核苷酸多样性指数明显小于下游的太平溪和巴南群体,从群体结构图中也可以看出,瓦氏黄颡鱼下游的太平溪和巴南群体有3种遗传背景,遗传来源较丰富,而上游的岷江口、宜宾、邵女坪和冯家坪群体仅有1种遗传背景,遗传来源较为单一。因此,瓦氏黄颡鱼在下游得到了较多的遗传补充,而在上游得到的遗传补充较少。类似的下游群体遗传多样性较上游丰富的现象在长江其他鱼类中也存在,如黑尾近红鲌(Ancherythroculter nigrocauda)(Zhai et al,2019)和 中 华 纹 胸 鮡(Glyptothorax sinensis)(张力文等,2020)。本研究中,瓦氏黄颡鱼上游群体具有更低的遗传多样性,更易发生遗传漂变作用,在鱼类遗传多样性保护中需要特别关注。长江上游个体之间的遗传相似度较高,这可能跟繁殖群体的减少相关,已有的研究表明瓦氏黄颡鱼资源呈现低龄化和小型化的趋势,3龄以上的繁 图5 瓦氏黄颡鱼核苷酸多样性指数与海拔(a)和坡降(b)的相关性分析。P 0.05表示显著相关。Fig.5 Correlation analysis of nucleotide diversity with elevation(a)and channel slope(b).P 0.05 indicates significant correlation.图6 遗传分化(a)与地理距离(b)和隔离时间的Mantel检验。Fst:遗传分化指数。Fig.6 Mantel test of genetic differentiation with geographical distance and isolated time.Fst,Index of genetic differentiation.3 讨论 熊飞等:基于基因组重测序的长江上游瓦氏黄颡鱼群体遗传结构研究 2023年|31卷|4期|22391|第8页 殖群体占比较少(杨家云,1994;熊飞等,2015)。虽然目前瓦氏黄颡鱼在渔获物中仍占较大的比例(吕浩等,2018),其高繁殖力能快速补偿损失的成熟个体,扩大群体规模,然而,局部种群的消失或大部分成熟个体的消失,会造成遗传多样性的丧失。较低的遗传多样性并不总是意味着种群数量的必然下降,长江上游大坝建设后,一些江段的水流变缓,瓦氏黄颡鱼的某些基因型可能更适应这种缓流的环境。但是,低遗传多样性对物种长期的生存和进化是不利的(Cook&Sullivan,2018)。与长江上游一些产漂流性卵鱼类遗传结构的同质性不同(Cheng et al,2013;Liu et al,2020),产粘性卵的瓦氏黄颡鱼存在明显的遗传结构。在被研究的长江上游群体中存在3个明显的遗传分支。宜宾与邵女坪和桧溪群体之间虽有向家坝阻隔,但为同一遗传背景;冯家坪与岷江口群体有向家坝和溪洛渡大坝阻隔,但这两个群体遗传背景也一致。瓦氏黄颡鱼的这种遗传结构可能是群体历史造成的,而不是大坝阻隔的影响。3.2 环境因素对瓦氏黄颡鱼群体结构的影响 河流坡降通常限制鱼类的迁移、扩散,从而影响鱼类的基因流和种群遗传结构。Kanno等(2011)研究表明,河流坡降是影响美洲红点鲑(Salvelinus fontinalis)种群遗传结构的主要因子。本研究结果显示坡降对瓦氏黄颡鱼SNP数量和核苷酸多样性指数造成了显著影响,水流的高度落差可能是阻碍瓦氏黄颡鱼扩散和基因交流的原因之一。瓦氏黄颡鱼不同江段的群体遗传分化较大,但是Mantel检验遗传分化指数与地理距离无显著相关性。宜宾与岷江口地理距离相隔较近,但是遗传分化值达到了0.233。冯家坪与岷江口相距较远,遗传分化值仅为0.016,在这些群体中,地理距离似乎并不是群体遗传分化的主要原因,奠基者效应和生境差异则更可能造成了目前的遗传分化格局。冯家坪到太平溪江段最大距离达1,300 km左右,海拔落差470 m左右,其间地形复杂,流域内环境时空变化较大,不同局域环境的适应性差异可能产生群体遗传分化(Xu et al,2017)。长江上游瓦氏黄颡鱼群体遗传分化与隔离时间相关关系不显著。一般认为,短期的时间尺度造成的遗传分化比较难检测出来。Deiner等(2007)发现大坝建成后的虹鳟(Oncorhynchus mykiss)遗传结构没有发生大的变化(繁殖了510代),Vera-Escalona等(2015)研究表明,大坝阻隔未引起普拉特南乳鱼(Galaxias platei)产生明显的遗传分化(繁殖了2030代)。本研究长江上游区域大坝截流时间约20年的时间或更短,按照瓦氏黄颡鱼最小性成熟年龄2龄一个世代计算,大约经历了10个或更少的世代,短时间内产生的遗传隔离效应还没有对瓦氏黄颡鱼遗传结构造成影响。高通量测序技术为鱼类遗传分析提供了大量的分子标记,为群体遗传模式的研究提供了强有力的工具。本研究鉴定了瓦氏黄颡鱼大量的SNP位点,阐明了长江上游瓦氏黄颡鱼群体遗传多样性和遗传结构。瓦氏黄颡鱼上游群体具有更低的遗传多样性,更易发生遗传漂变,在鱼类遗传多样性保护中需要特别关注。长江上游瓦氏黄颡鱼存在显著的遗传结构,可分为3种不同的遗传分支,不同遗传分支之间分化较大,可视为3个不同遗传单元进行种质资源管理。ORCID 熊飞 https:/orcid.org/0000-0002-0692-6862 参考文献 Alexander DH,Novembre J,Lange K(2009)Fast model-based estimation of ancestry in unrelated individuals.Genome Research,19,16551664.Chen DQ,Liu SP,Duan XB,Xiong F(2002)A preliminary study of the fisheries biology of main commercial fishes in the middle and upper reaches of the Yangtze River.Acta Hydrobiologica Sinica,26,618622.(in Chinese with English abstract)陈大庆,刘绍平,段辛斌,熊飞(2002)长江中上游主要经济鱼类的渔业生物学特征.水生生物学报,26,618622.Cheng F,Li W,Castello L,Murphy BR,Xie SG(2015)Potential effects of dam cascade on fish:Lessons from the Yangtze River.Reviews in Fish Biology and Fisheries,25,569585.Cheng F,Li W,Wu QJ,Hallerman E,Xie SG(2013)Microsatellite DNA variation among samples of bronze gudgeon,Coreius heterodon,in the mainstem of the Yangtze River,China.Ichthyological Research,60,165171.Cook DR,Sullivan SMP(2018)Associations between riffle development and aquatic biota following low head dam removal.Environmental Moni