基于微卫星标记对长江下游鲢遗传多样性现状的分析

罗宇婷, 方弟安, 周彦锋, 徐东坡, 彭云鑫, 彭飞, 张桂宁, 刘凯, 尤洋

罗宇婷, 方弟安, 周彦锋, 徐东坡, 彭云鑫, 彭飞, 张桂宁, 刘凯, 尤洋. 基于微卫星标记对长江下游鲢遗传多样性现状的分析[J]. 南方水产科学, 2021, 17(6): 48-57. DOI: 10.12131/20210067
引用本文: 罗宇婷, 方弟安, 周彦锋, 徐东坡, 彭云鑫, 彭飞, 张桂宁, 刘凯, 尤洋. 基于微卫星标记对长江下游鲢遗传多样性现状的分析[J]. 南方水产科学, 2021, 17(6): 48-57. DOI: 10.12131/20210067
LUO Yuting, FANG Di'an, ZHOU Yanfeng, XU Dongpo, PENG Yunxin, PENG Fei, ZHANG Guining, LIU Kai, YOU Yang. Genetic diversity of silver carp (Hypophthalmichthys molitrix) in lower reaches of Yangtze River based on microsatellite markers[J]. South China Fisheries Science, 2021, 17(6): 48-57. DOI: 10.12131/20210067
Citation: LUO Yuting, FANG Di'an, ZHOU Yanfeng, XU Dongpo, PENG Yunxin, PENG Fei, ZHANG Guining, LIU Kai, YOU Yang. Genetic diversity of silver carp (Hypophthalmichthys molitrix) in lower reaches of Yangtze River based on microsatellite markers[J]. South China Fisheries Science, 2021, 17(6): 48-57. DOI: 10.12131/20210067

基于微卫星标记对长江下游鲢遗传多样性现状的分析

基金项目: 国家重点研发计划项目 (2018YFD0900905);国家科技资源共享服务平台建设运行项目“国家淡水水产种质资源库” (2021DKA30470)
详细信息
    作者简介:

    罗宇婷 (1997—),女,硕士研究生,研究方向为渔业资源。E-mail: 996625504@qq.com

    通讯作者:

    方弟安 (1978—),男,博士,研究员,从事渔业资源研究。E-mail: fangdian@ffrc.cn

    尤 洋 (1967—),男,研究员,从事渔业资源研究。E-mail: youy@ffrc.cn

  • 中图分类号: S 932

Genetic diversity of silver carp (Hypophthalmichthys molitrix) in lower reaches of Yangtze River based on microsatellite markers

  • 摘要: 为了解长江下游鲢 (Hypophthalmichthys molitrix) 遗传多样性和遗传结构分布,该研究利用微卫星分子标记结合毛细管电泳技术,通过对2017—2019年长江湖口、安庆、芜湖、当涂、镇江、靖江、张家港和常熟8个江段鲢野生群体共327 份DNA进行基因分型与遗传多样性相关分析,比对了其遗传多样性、遗传分化和遗传结构情况。结果显示,长江下游8 个群体总体遗传多样性水平较高,但群体内缺乏足够的杂合子,8个群体等位基因数 (Na)介于6.0~12.3,有效等位基因数(Ne)介于3.94~6.10,观测杂合度 (Ho)和期望杂合度 (He)分别介于0.617~0.719和0.658~0.774,无偏期望杂合度(uHe)介于0.671~0.782,Shannon's信息指数 (I)介于1.38~1.86,群体内近交系数 (F)介于0.001~0.174;遗传分化分析显示Fst (F-statitics values) 介于0.006~0.068,基因流 (Nm)介于3.41~41.9,遗传距离 (D)介于0.001~0.106,湖口群体与其他群体表现出中等程度的遗传差异,其他群体间遗传差异很小,分子方差分析显示变异主要来自群体内 (97.6%);遗传结构分析显示湖口群体有不同于其他群体的基因库。结果表明江湖连通性可能对于鲢在湖口江段展现独特遗传分化特征具有积极意义,而形成这种种群分化究竟是由于地理距离的空间隔离还是生境差异的生理屏障有待进一步探究。
    Abstract: In order to understand the genetic diversity and distribution of silver carp (Hypophthalmichthys molitrix) in the lower reaches of the Yangtze River, we used simple sequence repeats (SSR) combining with capillary electrophoresis to analyze 327 DNA samples of silver carp from eight wild populations in Yangtze River sections (Hukou, Anqing, Wuhu, Dangtu, Zhenjiang, Jingjiang, Zhangjiagang and Changshu sections) from 2017 to 2019, so as to investigate their genetic diversity, genetic differentiation and genetic structure. The results show that the genetic diversity of the eight populations was high. However, heterozygotes were not enough in the populations. The number of alleles ranged from 6.00 to 12.3, while the number of effective alleles ranged from 3.94 to 6.10. The observed heterozygosity and expected heterozygosity ranged from 0.617 to 0.719 and from 0.658 to 0.774, respectively. The unbiased expected heterozygosity ranged from 0.671 to 0.782. The Shannon's information index ranged from 1.38 to 1.86, and the inbreeding coefficient ranged from 0.001 to 0.174. The genetic differentiation analysis shows that the Fst value ranged from 0.006 to 0.068, while the gene flow ranged from 3.41 to 41.9. The genetic distance ranged from 0.001 to 0.106. Hukou population showed a moderate degree of genetic difference with the other populations, while the genetic differences among the other populations were small. Analysis of molecular variance shows that the variation mainly was within the population (97.6%). Genetic structure analysis shows that Hukou population had different gene pool from the other populations. The river-lake connectivity may be of positive significance for silver carp to show unique genetic differentiation characteristics in Hukou section of the Yangtze River. Therefore, whether the formation of this population differentiation is caused by the spatial isolation of geographical distance or the physiological barrier of habitat difference needs to be further studied.
  • 鱼类生长激素(growth hormone,GH)是由脑垂体细胞合成和分泌的一种单链多肽,其分子量在22 kDa左右,可分为信号肽和成熟肽2个部分,信号肽一般由16~17个氨基酸残基组成,成熟蛋白由173~190个氨基酸组成。GH对鱼类最明显的作用是刺激鱼体重和体长的增长,促进发育,参与鱼体代谢调节,还可以调节海河洄游性鱼类的渗透压,并促进其对海水适应能力[1]。GH不但可以作为激素对鱼类的生长发育和生殖起重要作用,而且可以作为一种细胞分裂素,在免疫反应起作用,同时还参与调节蛋白质,糖和脂肪等的合成和释放,从而影响蛋白质、糖和脂肪的代谢[2-4]。深入研究鱼类生长激素的表达调控规律,弄清鱼类生长发育的调控机理,可为人工调控鱼类生长发育,促进水产养殖业的发展提供理论依据。

    鲮(Cirrhinus molitorella)隶属鲤形目(Cypriniformes)鲤科(Cyprinidae)中的野鲮亚科(Labeoninae),广泛分布于中国华南地区,东南亚和非洲的热带及亚热带地区[5]。鲮是中国华南地区重要的淡水养殖品种之一,其产量占广东、广西池塘养殖鱼类总产量的35%左右,具有产量高、抗病力强和肉质鲜美的特点,其适合制作罐装食品和鱼丸,是中国淡水鱼加工附加值最高的品种之一。但是鲮生长速度慢,且不耐低温,是目前制约鲮养殖产业发展的因素之一。与哺乳动物类似,鱼类生长发育主要受GH-IGF(胰岛素样生长因子)的调控,深入开展鲮生长激素表达调控机理的研究,以便为人工控制鲮的生长发育提供理论依据。目前该试验室已克隆了鲮生长激素cDNA并在大肠杆菌中进行了重组表达。为进一步研究在各种条件下鲮血液和脑垂体中生长激素蛋白的表达规律,必须要制备特异性强、灵敏度高的抗鲮生长激素抗体。因此,该研究在重组表达了鲮生长激素蛋白的基础上,制备了兔抗鲮生长激素多克隆抗体,并对其灵敏度和特异性进行了分析。

    重组工程菌M15(pQE30-GH)由笔者实验室构建,新西兰大白兔(Oryctolagus cuniculus)购于广州中医药大学实验动物中心,鲮鱼购于广州鹭江市场,镍琼脂糖凝胶FF填料购自北京韦氏博慧色谱科技有限公司,HRP标记的羊抗兔IgG购于Amersham公司,DAB显色试剂盒购自武汉博士德生物工程有限公司,其它试剂均为分析纯。

    取含有500 mL LB(100 μg · mL-1 Amp和25 μg · mL-1 Km)及种菌的烧瓶,最佳诱导条件参照江世贵等[6]对鲮生长激素cDNA的克隆及其重组表达产物的促生长活性的研究。在37℃加入浓度为100 mmol · L-1的IPTG溶液使其终浓度1 mmol ·L-1 250 rpm振荡培养诱导5 h,4℃保存。菌体经超声波破碎后,4℃,8 000 rpm离心20 min,离心收集到的包涵体经TE1(Tris-Cl 10 mmol · L-1,pH 8.0,EDTA 1 mmol · L-1),TE2[1%TritonX-100(V/V),Tris-Cl 10 mmol · L-1,pH 8.0,EDTA 1 mmol · L-1]和TE3(尿素2 mol · L-1,Tris-Cl 10 mmol · L-1 pH 8.0,EDTA 1 mmol · L-1)洗涤后,即可得到初步纯化的包涵体蛋白,然后将包涵体蛋白溶解在变性液(尿素8 mol · L-1,Tris-Cl 20 mmol ·L-1 pH 8.0,0.5 mol · L-1 NaCl,5 mmol · L-1 imidazole,1 mmol · L-1 2-mercaptoethanol)中。

    经过包涵体分离得到的包涵体蛋白,在4℃,8 000 rpm离心10 min后,上清液用0.2 μm的过滤器过滤,然后上样经镍琼脂糖凝胶FF(装填后柱床长5 cm,直径1.6 cm)亲合层析。上柱后以含20 mmol · L-1,500 mmol · L-1的咪唑以及ddH2O进行阶段洗脱,流速1.5 mL · min-1,手工收集洗脱峰。将收集后的样品进行SDS-PAGE分析,同时将纯化后的蛋白用稀释变性液(尿素8 mol · L-1,Tris-Cl 20 mmol · L-1 pH8.0,0.5 mol · L-1 NaCl)稀释到1 mg · mL-1装入处理好的透析袋中,采用逐步稀释的方法复性。具体操作如下,将装有蛋白的透析袋放到等体积的复性液(氧化型谷胱甘肽0.1 mmol · L-1,还原型谷胱甘肽0.9 mmol · L-1,Tris-Cl 20 mmol ·L-1,pH 8.0)中,4℃搅拌过夜,使尿素浓度由4~2.5 mol · L-1,再到0.8 mol · L-1逐步稀释,然后在ddH2O(pH 8.0)透析过夜以除去剩余的变性剂。最后,将复性后的蛋白用超滤离心管过滤浓缩。

    将浓缩后的重组鲮GH作为抗原,采用改进的方法[7]免疫新西兰大白兔。在第1天和第3天,将1 mL样品A(300 μg · 500 μL-1纯化浓缩蛋白和500 μL Freund′s完全佐剂)采用多点注射到兔子的背部和腿部;第28天将1 mL样品B(300 μg · 500 μL-1纯化浓缩蛋白和500 μL Freund′s不完全佐剂)加强免疫1次。第35天,颈动脉取血,37℃放置2 h,4℃过夜,离心收集血清。

    采用间接ELISA检测抗血清效价。用pH 9.6,0.05 mol · L-1的碳酸盐缓冲液将纯化的重组鲮生长激素稀释为20 μg · mL-1,包被酶标板,每孔100 μL,4℃包被过夜。次日每孔加入100 μL 5%的脱脂奶粉,37℃封闭1 h后,用pH 7.4的PBST洗涤3次。然后每孔加入200~51 200倍稀释的兔抗血清100 μL,阴性对照为1:50倍稀释的免疫前血清,空白对照为兔抗血清稀释液(PBS),每份样品均做1平行,37℃孵育1 h后,同上洗涤3次。再次每孔加入100 μL HRP标记的羊抗兔IgG,37℃孵育1 h后,同上洗涤3次。接下来每孔加入100 μL TMB显色液,37℃避光反应20 min后,每孔再加入50 μL 2 mol · L-1硫酸终止液。最后,在酶标仪上测450 nm下的OD值。

    将抗原稀释为0.5、1、2、4和10 ng,经15% SDS-PAGE电泳分离后转移到硝酸纤维素膜上。用稀释度1:500的抗血清作为一抗[8],比例1:2 000作为二抗稀释度,DAB显色。标准分子量蛋白质转膜后用丽春红染色液染色。

    鲮组织提取物中GH同源蛋白质的检测。取100 mg脑垂体组织加入到100 μL抽提液(4 M Urea,0.5% W/V SDS,10 mM EDTA和2 mM PMSF)中匀浆,94℃孵育10 min,10 000 g离心10 min后,用Bradford方法进行定量测定[9-10]。取鲮血4℃,6 500 rpm离心20 min分离血清,同样用Bradford方法进行定量测定。分别取1 μg脑垂体提取蛋白和1 μg血清蛋白进行Western blotting分析。

    图 1所示,IPTG可诱导M15(pQE30-GH)在大肠杆菌中大量表达。SDS-PAGE凝胶电泳结果显示,1条24.5 kDa的诱导表达重组鲮GH带。经过包涵体分离、纯化、浓缩得到纯度很高的蛋白。

    图  1  鲮鱼生长激素的重组融合表达,包涵体的分离,纯化和浓缩电泳图
    M. 低分子量蛋白标准;1. IPTG诱导生长激素的表达;2. 生长激素表达负对照;3. M15(pQE30)诱导表达;4. 包涵体分离后电泳图;5. 包涵体纯化后电泳图;6. 纯化蛋白浓缩后电泳图
    Fig. 1  The electrophoresis map of the recombinant expression and purification of mud carp growth hormone
    M. low molecular weight protein marker; 1. the growth hormone induced by IPTG; 2. growth hormone expressed negative control; 3. M15(pQE30)expressed; 4. inclusion body; 5. purified inclusion body; 6. condensed protein

    兔抗鲮生长激素的血清进行间接ELISA检测时,当(样品吸光值-空白吸光值)/(阴性吸光值-空白吸光值)>2.1即认为是阳性[11],经计算,此研究制备的兔抗血清效价为1 : 12 800(表 1)。

    表  1  ELISA检测结果
    Table  1  The results of ELISA
    免疫血清稀释倍数
    the antiserum dilution times
    免疫血清OD450
    the antiserum OD450
    阴性对照血清OD450
    normal rabbit serum OD450
    空白对照OD450
    black control OD450
    P/N
    1:200 1.186 0.336 0.175 +
    1:400 1.138 0.336 0.175 +
    1:800 1.180 0.336 0.175 +
    1:1 600 1.029 0.336 0.175 +
    1:3 200 0.872 0.336 0.175 +
    1:6 400 0.656 0.336 0.175 +
    1:12 800 0.533 0.336 0.175 +
    1:25 600 0.464 0.336 0.175
    1:51 200 0.379 0.336 0.175
    下载: 导出CSV 
    | 显示表格

    将抗原稀释为0.5、1、2、4和10 ng,经15% PAGE电泳分离后,转移到硝酸纤维素膜上。用稀释度1:500的抗血清作为一抗,比例1:2 000作为二抗稀释度,DAB显色。标准分子量蛋白质转膜后用丽春红染色液染色。结果显示,可检测到10 ng抗原量,表明此抗体灵敏度较高(图 2)。

    图  2  抗体灵敏度分析
    M. 低分子量蛋白标准;1. 10 ng抗原;2. 4 ng抗原;3. 2 ng抗原;4. 1 ng抗原;5. 0.5 ng抗原
    Fig. 2  The antibody sensitivety analysis
    M. low molecular weight protein marker; 1. 10 ng antigen; 2. 4 ng antigen; 3. 2 ng antigen; 4. 1 ng antigen; 5. 0.5 ng antigen

    为了解这一抗血清的特异性,取鲮脑垂体组织和血清进行Western blotting分析。结果显示,用此多抗可检测到1条大约为21.5 kDa的条带(图 3),这说明这一多抗具有较好的免疫特性。

    图  3  抗体特性性分析
    M. 低分子量蛋白标准;1. 血清;2. 脑组织
    Fig. 3  The antibody specificity analysis
    M. low molecular weight protein marker; 1. serum; 2. pituitary

    此试验所表达的重组鲮生长激素主要以包涵体形式存在。另外,包涵体中的杂蛋白较少、且50%以上为重组蛋白,给重组蛋白的纯化带来便利。经过3次TE缓冲液洗涤,包涵体达到了较好的纯化效果,洗涤后重组蛋白纯度达到85%以上。初步纯化后的包涵体重组蛋白经镍琼脂糖凝胶FF亲和层析后,由含20 mmol · L-1,500 mmol · L-1的咪唑以及ddH2O进行阶段洗脱,进一步纯化该重组蛋白,纯度可达到95%以上。由于GH基因内Cys含量较高,分子内可形成2个二硫键,在复性时容易形成蛋白质聚集体,所以笔者利用GSH/GSSG氧化交换系统促使二硫键的正确配对,同时采用稀释复性和透析复性相结合的方式对纯化蛋白进行复性。在复性过程中融合蛋白的再聚合是令人棘手的问题。姚燕等[12]中华绒鳌蟹(Eriocheir japonica sinensis)蜕皮抑制激素基因的表达及抗体制备研究表明在可在复性液中加入甘油增加粘度来减少蛋白分子间的相互作用,从而阻止透析复性过程中蛋白的再聚合,并且L-精氨酸作为小分子促溶剂,使整个透析过程中蛋白的聚合几率降到最低。张傅山等[13]研究表明在复性液中加入10%无水乙醇可提高复性效率。另外,蛋白质浓度也是影响蛋白质复性的重要因素。当蛋白浓度超过1 mg · mL-1时,即使复性液中加入甘油和L-精氨酸,复性后透析时仍会出现浑浊现象。笔者在试验中并没有加入促溶剂,直接将纯化后的蛋白用复性稀释液稀释到1 mg · mL-1再进行复性也得到了较好的复性效果。

    经过包涵体分离的蛋白通过镍琼脂糖凝胶FF填料来分离带His标签的重组蛋白,其纯度可达95%以上。分离纯化的蛋白质大小为24.5 kDa,这比GH成熟蛋白大了约3 kDa左右,这是因为6xHis标签的分子量约为3 kDa,而GH成熟蛋白的分子量约为22 kDa,所以该蛋白质带大小与预期大小一致。纯化的方法有很多,此试验采用的是亲和层析,该纯化方法简便,分离效果好,其纯化产物可达95%以上,且经过SDS-PAGE和考马斯亮蓝染色呈现均一的1条带。这说明该纯化产物可用于抗体的制备和后续的生化实验的需要。

    用经改良的方法进行免疫新西兰大白兔,经过35 d免疫注射后,用间接ELISA方法检测其抗体滴度为12 800,可以用来Western blotting分析。用改良后的方法来免疫新西兰大白兔在时间上由原来的3~4个月缩短到35 d,大大缩短了制备高效抗血清的时间,且所制备抗血清均能用于ELISA和Western blotting分析。胡志红等[14]用改良的方法制备了猴(Macaca mulatta)Esc-615基因多克隆抗体的制备,其抗血清滴度达到100 000。冯浩等[8]制备了青鱼(Mylopharyngodon piceus)生长激素多克隆抗血清也能用于ELISA和Western blotting分析。从而表明了上诉方法制备抗鲮GH的多克隆抗体不仅简便,而且行之有效,可为生长激素活性及功能的研究提供有效免疫学鉴定工具。

    此试验制备的兔抗鲮生长激素多克隆抗体,经ELISA检测其滴度为12 800,经Western blotting分析,可以检测到10 ng的抗原量,而且可以在脑垂体提取物和血清中检测到大小为21.5 kDa的蛋白质,充分说明该抗体具有较强的灵敏性和免疫原性,为进一步研究在不同条件下鲮血液和脑垂体中生长激素蛋白的表达规律奠定了基础。

  • 图  1   长江下游鲢采样位点

    Figure  1.   Sampling sites of H. molitrix in lower Yangtze River

    图  2   微卫星位点HLJBL168、HLJBL184、HLJBL220基因分型结果

    Figure  2.   Genotype of microsatellite loci HLJBL168, HLJBL184 and HLJBL220

    图  3   基于遗传距离的长江下游8个鲢群体UPGMA聚类图

    Figure  3.   UPGMA tree based on genetic distance among eight silver carp populations in lower Yangtze River

    图  4   Delta KK的变化曲线

    Figure  4.   Curve of change of Delta K with changing K value

    图  5   长江下游8个鲢群体聚类结果图 (K=4)

    Figure  5.   Clustering diagram of eight populations of H. molitrix in lower Yangtze River (K=4)

    表  1   长江下游鲢采样情况

    Table  1   Sampling of H. molitrix in lower Yangtze River

    地点
    Location
    简称
    Code
    坐标
    Coordinates
    采样时间段
    Sampling time
    数量
    Number
    湖口 Hukou HK 116°15'32''E, 29°47'02''N 2018.06—2019.07 26
    安庆 Anqing AQ 117°0'18''E, 30°29'28''N 2018.10—2019.12 53
    芜湖 Wuhu WH 118°19’59″E, 31°23’58″N 2017.03—2018.04 21
    当涂 Dangtu DT 118°27'3''E, 31°34'10''N 2017.03—2018.04 27
    镇江 Zhenjiang ZJ 119°20'45''E, 32°11'42''N 2017.03—2018.04 51
    靖江 Jingjiang JJ 120°17'31''E, 31°57'30''N 2017.03—2018.04 30
    张家港 Zhangjiagang ZJG 120°26'55''E, 31°59'20''N 2017.03—2018.04 67
    常熟 Changshu CS 120°54'12''E, 31°46'34''N 2017.03—2018.04 52
    下载: 导出CSV

    表  2   基于11个微卫星标记的鲢遗传多样性水平

    Table  2   Genetic diversity level of H. molitrix based on 11 simple sequence repeats

    基因座
    Locus
    等位
    基因数
    Na
    观测
    杂合度
    Ho
    期望
    杂合度
    He
    多态
    信息含量
    PIC
    哈迪-温伯格
    平衡
    HWE
    HLJBL16580.5940.6690.618ns
    HLJBL16870.6260.6910.630ns
    HLJBL169100.5290.5120.479ns
    HLJBL17090.5250.6700.625***
    HLJBL174210.8500.9160.909**
    HLJBL176170.7930.8820.869ns
    HLJBL184170.6700.7290.706*
    HLJBL203140.7420.7440.700ns
    HLJBL216170.6370.8860.874***
    HLJBL217250.7450.9360.930**
    HLJBL220190.6570.8650.851***
    平均值 Mean14.9090.6700.7730.744
    注:***. 极显著偏离哈迪-温伯格平衡;**. 中等显著偏离哈迪-温伯格平衡;*. 轻微显著偏离哈迪-温伯格平衡;ns. 未偏离哈迪-温伯格平衡。 Note: ***. Extremely significantly deviated from the Hardy-Weinberg equilibrium; **. Moderate significantly deviated from the Hardy-Weinberg equilibrium; *. Slightly significantly deviated from the Hardy-Weinberg equilibrium; ns. Non-significantly deviated from the Hardy-Weinberg equilibrium.
    下载: 导出CSV

    表  3   长江下游8个江段鲢的遗传多样性水平

    Table  3   Genetic diversity of H. molitrix in eight sections of lower Yangtze River

    群体
    Population
    等位
    基因
    Na
    有效等位
    基因
    Ne
    私有等位
    基因
    Ar
    观测
    杂合度
    Ho
    期望
    杂合度
    He
    无偏期望
    杂合度
    uHe
    Shannon's
    信息指数
    I
    群体内近交
    系数
    F
    湖口 HK 6.000±0.853 3.943±0.682 2 0.648±0.078 0.658±0.069 0.671±0.070 1.375±0.175 0.011±0.051
    安庆 AQ 12.273±2.273 5.867±1.272 17 0.718±0.039 0.751±0.042 0.758±0.042 1.832±0.202 0.040±0.019
    芜湖 WH 8.636±1.357 5.286±1.021 1 0.719±0.057 0.725±0.052 0.743±0.054 1.673±0.189 0.001±0.047
    当涂 DT 6.818±0.861 4.770±0.679 2 0.692±0.057 0.741±0.037 0.755±0.038 1.597±0.145 0.065±0.056
    镇江 ZJ 11.273±2.149 6.102±1.348 6 0.691±0.042 0.755±0.040 0.762±0.041 1.808±0.205 0.080±0.035
    靖江 JJ 8.909±1.665 5.206±1.192 3 0.653±0.046 0.717±0.046 0.730±0.047 1.645±0.194 0.087±0.034
    张家港 ZJG 11.273±1.711 5.927±1.009 3 0.617±0.036 0.772±0.036 0.778±0.037 1.848±0.177 0.174±0.072
    常熟 CS 11.000±1.673 5.962±1.012 7 0.655±0.038 0.774±0.037 0.782±0.037 1.859±0.174 0.137±0.060
    平均值 Mean 9.523±0.603 5.383±0.365 5.123 0.674±0.018 0.737±0.016 0.747±0.016 1.705±0.064 0.074±0.018
    下载: 导出CSV

    表  4   基于分子方差分析法 (AMOVA) 的8个鲢群体遗传变异结果

    Table  4   Molecular variance (AMOVA) results of eight H. molitrix populations

    变异来源
    Source of variation
    自由度
    df
    离差平方和
    Sum of squares
    方差分量
    Variance components
    方差分量比例
    Percentage of variation/%
    群体间 Among populations 7 82.294 0.097 2.37
    群体内 Within populations 646 2 581.102 3.996 97.63
    总变异 Total variation 653 2 663.396 4.093 100.00
    下载: 导出CSV

    表  5   基于微卫星标记的鲢群体Fst (对角线以下) 和基因流 (对角线以上) 配对估计值

    Table  5   Estimated pairwise Fst (below diagonal) and Nm (above diagonal) of H. molitrix based on simple sequence repeats

    湖口
    HK
    安庆
    AQ
    芜湖
    WH
    当涂
    DT
    镇江
    ZJ
    靖江
    JJ
    张家港
    ZJG
    常熟
    CS
    湖口 HK 4.265 4.553 3.748 3.538 3.674 3.885 3.409
    安庆 AQ 0.055 21.732 16.276 41.851 25.488 24.958 24.452
    芜湖 WH 0.052 0.011 11.802 18.199 17.587 15.180 15.284
    当涂 DT 0.063 0.015 0.021 17.548 17.505 12.762 14.626
    镇江 ZJ 0.066 0.006 0.014 0.014 28.015 27.298 33.488
    靖江 JJ 0.064 0.010 0.014 0.014 0.009 18.377 16.173
    张家港 ZJG 0.060 0.010 0.016 0.019 0.009 0.013 38.421
    常熟 CS 0.068 0.010 0.016 0.017 0.007 0.015 0.006
    下载: 导出CSV

    表  6   长江下游8个鲢群体遗传距离的配对矩阵

    Table  6   Genetic distance Pairwise matrix of eight H. molitrix populations in lower Yangtze River

    湖口
    HK
    安庆
    AQ
    芜湖
    WH
    当涂
    DT
    镇江
    ZJ
    靖江
    JJ
    张家口
    ZJG
    常熟
    CS
    湖口 HK
    安庆 AQ0.082
    芜湖 WH0.0780.005
    当涂 DT0.0920.0160.017
    镇江 ZJ0.0990.0010.0070.014
    靖江 JJ0.0940.0030.0030.0080.001
    张家港 ZJG0.0870.0100.0130.0250.0090.009
    常熟 CS0.1050.0090.0100.0200.0040.0110.003
    下载: 导出CSV
  • [1] 倪勇, 伍汉霖. 江苏鱼类志[M]. 北京: 中国农业出版社, 2006: 315-317.
    [2] 刘飞, 林鹏程, 黎明政, 等. 长江流域鱼类资源现状与保护对策[J]. 水生生物学学报, 2019, 43: 144-156.
    [3] 杨习文. 基于微卫星标记的长江江苏段鲢 (Hypophthalmichthys molitrix) 增殖放流效果评估[D]. 上海: 上海海洋大学, 2019: 50-52.
    [4] 张敏莹, 徐东坡, 刘凯, 等. 长江下游放流鲢群体遗传多样性的微卫星标记分析[J]. 江西农业大学学报, 2012, 34(1): 141-146. doi: 10.3969/j.issn.1000-2286.2012.01.026
    [5] 朱晓东, 耿波, 李娇, 等. 利用30个微卫星标记分析长江中下游鲢群体的遗传多样性[J]. 遗传, 2007, 29(6): 705-715. doi: 10.3321/j.issn:0253-9772.2007.06.010
    [6] 王长忠, 梁宏伟, 邹桂伟, 等. 长江中上游两个鲢群体遗传变异的微卫星分析[J]. 遗传, 2008, 30(10): 1341-1348. doi: 10.3321/j.issn:0253-9772.2008.10.016
    [7] 沙航, 罗相忠, 李忠, 等. 基于COI序列的长江中上游鲢6个地理群体遗传多样性分析[J]. 中国水产科学, 2018, 25(4): 783-792.
    [8] 陈会娟, 刘明典, 汪登强, 等. 长江中上游4个鲢群体遗传多样性分析[J]. 淡水渔业, 2018, 48(1): 20-25. doi: 10.3969/j.issn.1000-6907.2018.01.004
    [9] 陈会娟. 长江中游四大家鱼放流亲本对早期资源和遗传多样性的影响研究[D]. 重庆: 西南大学, 2019: 85.
    [10] 许承双, 艾志强, 肖鸣. 影响长江四大家鱼自然繁殖的因素研究现状[J]. 三峡大学学报 (自然科学版), 2017, 39(4): 27-30.
    [11] 刘熠, 杨习文, 任鹏, 等. 长江湖口段春夏季仔稚鱼群落结构研究[J]. 水生生物学报, 2019, 43(1): 141-154.
    [12] 杨习文, 刘熠, 薛向平, 等. 基于微卫星标记的长江江苏段鲢(Hypophthalmichthys molitrix)增殖放流资源贡献率的评估[J]. 湖泊科学, 2020, 32(4): 1154-1164. doi: 10.18307/2020.0422
    [13]

    KALINOWSKI S T, TAPER M L, MARSHALL T C. Revising how the computer program CERVUS accommodates genotyping error increases success in paternity assignment[J]. Mol Ecol, 2007, 16(5): 1099-1106. doi: 10.1111/j.1365-294X.2007.03089.x

    [14]

    PEAKALL R, SMOUSE P E. GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research[J]. Mol Ecol Notes, 2006, 6(1): 288-295. doi: 10.1111/j.1471-8286.2005.01155.x

    [15]

    EXCOFFIER L, LISCHER H. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows[J]. Mol Ecol Resour, 2010, 10(3): 564-567. doi: 10.1111/j.1755-0998.2010.02847.x

    [16]

    TAMURA K, PETERSON D, PETERSON N, et al. MEGA 5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods[J]. Mol Biol Evol, 2011, 28(10): 2731-2739. doi: 10.1093/molbev/msr121

    [17]

    SMOUSE P E, PEAKALL R. Spatial autocorrelation analysis of individual multiallele and multilocus genetic structure[J]. Heredity, 2010, 82(5): 561-573.

    [18]

    EVANNO G, REGNAUT S, GOUDET J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation studies[J]. Mol Ecol, 2005, 14(8): 2611-2620. doi: 10.1111/j.1365-294X.2005.02553.x

    [19]

    JAKOBSSON M, ROSENBERG N A. Clumpp: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure[J]. Bioinformatics, 2007, 23(14): 1801-1806. doi: 10.1093/bioinformatics/btm233

    [20]

    SWIFT M L. Graphpad prism, data analysis, and scientific graphing[J]. J Chem Inf Model, 1997, 37(2): 411-412.

    [21]

    BOTSTEIN D, WHITER R L, SKOLNICK M, et al. Construction of a genetic linkage map in man using restriction fragment length polymorphisms[J]. Am J Hum Genet, 1980, 32(3): 314-331.

    [22] 胡玉婷, 江河, 段国庆, 等. 安徽两水系黄颡鱼的微卫星遗传多样性分析[J]. 南方水产科学, 2020, 16(5): 33-41. doi: 10.12131/20200063
    [23]

    SHETE S, TIWARI H, ELSTON R C. On estimating the heterozygosity and polymorphism information content value[J]. Theor Popul Biol, 2000, 57(3): 265-271. doi: 10.1006/tpbi.2000.1452

    [24]

    WAPLES R S. A bias correction for estimates of effective population size based on linkage disequilibrium at unlinked gene loci[J]. Conserv Genet, 2006, 7(2): 167-184. doi: 10.1007/s10592-005-9100-y

    [25]

    KALINOWSKI S T. Counting alleles with rarefaction: private alleles and hierarchical sampling designs[J]. Conserv Genet, 2004, 5(4): 539-543. doi: 10.1023/B:COGE.0000041021.91777.1a

    [26] 姬长虹, 谷晶晶, 毛瑞鑫, 等. 长江、珠江、黑龙江水系野生鲢遗传多样性的微卫星分析[J]. 水产学报, 2009, 33(3): 364-371.
    [27]

    KALLINOWSKI S T. HP-RARE 1.0: a computer program for performing rarefaction on measures of allelic richness[J]. Mol Ecol Notes, 2005, 5(1): 187-189. doi: 10.1111/j.1471-8286.2004.00845.x

    [28]

    TAKEZAKI N, NEI M. Genetic distances and reconstruction of phylogenetic trees from microsatellite DNA[J]. Genetics, 1996, 144(1): 389-399. doi: 10.1093/genetics/144.1.389

    [29] 孙立元, 郭华阳, 朱彩艳, 等. 卵形鲳鲹育种群体遗传多样性分析[J]. 南方水产科学, 2014, 10(2): 67-71. doi: 10.3969/j.issn.2095-0780.2014.02.010
    [30] 李鹏飞, 刘萍, 柳学周, 等. 漠斑牙鲆引进种群同工酶的遗传多态性分析[J]. 中国水产科学, 2006, 13(1): 13-19. doi: 10.3321/j.issn:1005-8737.2006.01.003
    [31]

    MOSS D R, ARCES S M, OTOSHI C A, et al. Effects of inbreeding on survival and growth of Pacific white shrimp Penaeus (Litopenaeus) vannamei[J]. Aquaculture, 2007, 272(S1): S30-S37.

    [32]

    FANG D A, LUO Y T, XU D P, et al. Relationship between genetic risk and stock enhancement of the silver carp (Hypophthalmichthys molitrix) in the Yangtze River[J]. Fish Res, 2021, 235: 105829.

    [33]

    SHAKLEE J B, TAMARU C S, WAPLES R S. Speciation and evolution of marine fishes studied by the electrophoresis analysis of proteins[J]. Pac Sci, 1982, 36: 141-157.

    [34] 赵金良, 李思发. 长江中下游鲢、鳙、草鱼、青鱼种群分化的同工酶分析[J]. 水产学报, 1996, 20(2): 104-110.
    [35] 李思发, 吕国庆. 长江中下游鲢鳙草青四大家鱼线粒体DNA多样性分析[J]. 动物学报, 1998, 44(1): 82-93. doi: 10.3969/j.issn.1674-5507.1998.01.013
    [36] 张四明, 邓怀, 汪登强, 等. 长江水系鲢和草鱼遗传结构及变异性的RAPD研究[J]. 水生生物学报, 2001, 25(4): 325-330.
    [37]

    WEIGHT S. Evolution and the genetics of population variability within and among natural populations[M]. Chicago: University of Chicago Press, 1978: 4.

    [38]

    RAFAEL G A, BEATRIZ G, MARTÍNEZ S C G, et al. Extensive pollen flow but few pollen donors and high reproductive variance in an extremely fragmented landscape[J]. PLoS One, 2017, 7(11): 0049012.

    [39] 于悦. 长江、赣江、鄱阳湖鲢遗传多样性和群体遗传结构分析[D]. 武汉: 华中农业大学, 2016: 18.
  • 期刊类型引用(2)

    1. 王庚申,颜懿,李彤,董鹏生,谢建军,许文军,张德民,张化俊. 漂白粉消毒后对虾养殖源水细菌群落的响应特征研究. 海洋与湖沼. 2024(01): 202-212 . 百度学术
    2. 郑仕夫,徐慧敏,陈曦,裘丽萍,宋超,范立民,李丹丹,孟顺龙,徐跑. 水产养殖尾水处理技术的研究现状和发展趋势. 中国农学通报. 2024(12): 159-164 . 百度学术

    其他类型引用(1)

推荐阅读
免漂洗罗非鱼鱼糜凝胶品质改良及机理分析
罗颖莹 et al., 南方水产科学, 2025
南海岛礁鱼源产蛋白酶菌株热带芽孢杆菌btzb2的特性研究
胡晓娟 et al., 南方水产科学, 2024
草鱼subfatin分子鉴定及表达特性分析
杨博雅 et al., 南方水产科学, 2024
植物乳植杆菌通过抑制蛋白水解改善罗非鱼发酵鱼糜凝胶强度
崔巧燕 et al., 南方水产科学, 2024
环境微生物对水体中微囊藻毒素降解机制的研究进展
张赫 et al., 渔业研究, 2024
养殖环境及贝源溶藻弧菌mlst分型及其毒力基因、耐药性分析
ZHENG Yudong et al., PROGRESS IN FISHERY SCIENCES, 2024
Scrutinization of marangoni convective flow of dusty hybrid nanofluid with gyrotactic microorganisms and thermophoretic particle deposition
Abbas, Munawar et al., JOURNAL OF THERMAL ANALYSIS AND CALORIMETRY, 2024
Comprehensive investigation of cu2+ adsorption from wastewater using olive-waste-derived adsorbents: experimental and molecular insights
Elboughdiri, Noureddine et al., INTERNATIONAL JOURNAL OF MOLECULAR SCIENCES, 2024
Comparative laboratory wettability study of sandstone, tuff, and shale using 12-mhz nmr t1-t2 fluid typing: insight of shale
SPE JOURNAL
Isolation, identification, and characterization of an aspergillus niger bioflocculant-producing strain using potato starch wastewater as nutrilite and its application
PLOS ONE, 2018
Powered by
图(5)  /  表(6)
计量
  • 文章访问数: 
  • HTML全文浏览量: 
  • PDF下载量: 
  • 被引次数: 3
出版历程
  • 收稿日期:  2021-02-28
  • 修回日期:  2021-04-13
  • 录用日期:  2021-05-05
  • 网络出版日期:  2021-05-28
  • 刊出日期:  2021-12-05

目录

/

返回文章
返回