利用微卫星标记分析7个凡纳滨对虾引进群体子一代的遗传多样性

黄小帅, 徐煜, 胡晓娟, 徐武杰, 苏浩昌, 文国樑, 杨铿, 曹煜成

黄小帅, 徐煜, 胡晓娟, 徐武杰, 苏浩昌, 文国樑, 杨铿, 曹煜成. 利用微卫星标记分析7个凡纳滨对虾引进群体子一代的遗传多样性[J]. 南方水产科学, 2019, 15(1): 54-62. DOI: 10.12131/20180135
引用本文: 黄小帅, 徐煜, 胡晓娟, 徐武杰, 苏浩昌, 文国樑, 杨铿, 曹煜成. 利用微卫星标记分析7个凡纳滨对虾引进群体子一代的遗传多样性[J]. 南方水产科学, 2019, 15(1): 54-62. DOI: 10.12131/20180135
HUANG Xiaoshuai, XU Yu, HU Xiaojuan, XU Wujie, SU Haochang, WEN Guoliang, YANG Keng, CAO Yucheng. Genetic diversity analysis of first filial generation of seven introduced Litopenaeus vannamei populations using microsatellite DNA markers[J]. South China Fisheries Science, 2019, 15(1): 54-62. DOI: 10.12131/20180135
Citation: HUANG Xiaoshuai, XU Yu, HU Xiaojuan, XU Wujie, SU Haochang, WEN Guoliang, YANG Keng, CAO Yucheng. Genetic diversity analysis of first filial generation of seven introduced Litopenaeus vannamei populations using microsatellite DNA markers[J]. South China Fisheries Science, 2019, 15(1): 54-62. DOI: 10.12131/20180135

利用微卫星标记分析7个凡纳滨对虾引进群体子一代的遗传多样性

基金项目: 中国水产科学研究院南海水产研究所中央级公益性科研院所基本科研业务费专项资金资助 (2017YB29);现代农业产业技术体系建设专项资金 (CARS-48);广东省渔港建设和渔业产业发展专项 (A201701B06);广东省促进经济发展专项资金 (现代渔业发展用途) 省级项目 (SDYY-2018-02)
详细信息
    作者简介:

    黄小帅(1988—),男,博士,助理研究员,从事虾蟹分子生理及健康养殖研究。E-mail: tianyi001@126.com

    通讯作者:

    曹煜成(1979 — ),男,博士,副研究员,从事虾蟹健康养殖研究。E-mail: cyc_169@163.com

  • 中图分类号: S 917.4

Genetic diversity analysis of first filial generation of seven introduced Litopenaeus vannamei populations using microsatellite DNA markers

  • 摘要:

    为明晰中国凡纳滨对虾(Litopenaeus vannamei)引进群体子一代的遗传多样性特征,于广东的3个对虾主产区采集7个养殖群体的种苗样品,其均为国外引进亲虾繁育的子一代。将之分别命名为TH-A1、TH-A2、TH-B、US-C1、US-C2、US-C3和US-C4,以微卫星标记检测其遗传多样性。结果显示,7个群体在12个位点呈现不同程度的多态性,其平均等位基因数(Na)、期望杂合度(He)、观测杂合度(Ho)和多态信息含量(PIC)分别为3.333~6.167、0.477~0.670、0.370~0.505和0.414~0.623;44个群体位点显著偏离哈迪-温伯格平衡,和文章中He大于Ho的结果对应。聚类分析显示7个群体共分为3支,其中TH-A1为一支,US-C1、US-C2和TH-A2为一支,其余的聚为一支。结果表明此次采集的养殖群体种苗样品存在一定的遗传差异。该结果可为后续挖掘种苗遗传背景与其养殖性能的关联性提供参考。

    Abstract:

    Microsatellite DNA markers were used to assess the genetic diversity of first filial generation of seven introduced Litopenaeus vannamei populations, named TH-A1, TH-A2, TH-B, US-C1, US-C2, US-C3 and US-C4, respectively. The results show that the seven populations had varying polymorphism at 12 microsatellite loci, with the mean allele number (Na) ranging from 3.333 to 6.167. The mean values of expected (He) and observed (Ho) heterozygosity were 0.477–0.670 and 0.370–0.505, respectively. The polymorphic information content (PIC) for each population varied from 0.414 to 0.623. For Hardy-Weinberg equilibrium test, 44 out of 88 exhibited significant deviation (P<0.05), which corresponded with the fact thatHe was slightly higher than Ho. Cluster analysis shows that the seven populations were clustered into three branches (TH-A1; US-C1, US-C2 and TH-A2; the rest). It is indicated that the genetic characteristics vary in different cultured populations of L. vannamei. The results provide references for further exploration of relevance between genetic characteristics and practical performance of L. vannamei.

  • 卵形鲳鲹 (Trachinotus ovatus) 具有生长快、肉质鲜美、经济价值高等优点,是我国南部沿海地区重要的海水养殖品种之一。其生存和生长容易受到各种环境因素的影响,其中水体中的溶解氧含量是关键因素之一。在自然条件下,海水通过风和浪的物理搅动使大量的氧气溶解在水中。水生生物光合作用产生的氧气也是水中溶解氧的重要来源。但当海水交换不及时,或赤潮爆发后,大量死亡的藻类在分解过程中消耗大量氧气[1],导致卵形鲳鲹缺氧。一些寄生虫例如刺激隐核虫 (Cryptocaryon irritans)[2]和眼点淀粉卵涡鞭虫 (Amylocdinium ocellatum)[3]的寄生会破坏鳃组织,使鱼无法进行气体交换,导致卵形鲳鲹缺氧死亡。

    溶解氧的减少会对鱼类产生诸多影响,包括食物摄入量减少、生长速度减慢、游泳能力下降和影响鱼类的形态学特征和生存策略[4]。研究发现,低氧环境会导致尼罗罗非鱼 (Oreochromis niloticus) 的食欲下降并影响其生长速度[5];相比正常溶解氧的生存环境,虹鳟 (Oncorhynchus mykiss) 在低氧条件下游得更慢[6];长期缺氧会降低大西洋鲑 (Salmo salar) 卵的孵化率,甚至导致胚胎发育畸形或死亡[7]。一些鱼类为适应低氧环境,鳃组织的形态结构发生了变化。例如,鲫 (Carassius carassius)[8]、金鱼 (C. auratus)[9]和团头鲂 (Megalobrama amblycephalala)[10]通过重建鳃组织、暴露鳃板增加鳃的表面积以应对缺氧胁迫。

    DNA标记是一种非常可靠的分子标记,随机扩增多态性DNA (Random Amplified Polymorphic DNA, RAPD)[11]、扩增片段长度多态性 (Amplified Fragment Length Polymorphism, AFLP)[12-13]、简单序列重复 (Simple Sequence Repeat, SSR)[14]、单核苷酸多态性 (Single Nucleotide Polymorphism, SNP)[15-17]和插入-缺失 (InDel)[18-19]等分子标记已广泛应用于遗传育种研究。在不同的个体中,由于基因的插入或缺失,相应的基因和基因组区域会有不同的序列长度。这些突变被称为InDel,可以通过插入转座元件、相似重复副本之间的不相等交叉事件或简单序列复制的滑移形成[20],并可能表现为功能丧失或无意义突变[21]。InDel和SNP是动物基因组中最丰富和分布最广泛的变异性来源,均非常适合于构建遗传连锁图谱、全基因组关联分析和遗传多样性分析等其他遗传学研究。然而,在分子标记辅助育种的过程中,InDel比SNP更方便可取,因为InDel的多态性更易通过设计特异性引物、PCR扩增和琼脂糖凝胶电泳来显示。以往有多项畜牧和农作物研究使用了InDel标记[22-23],但在水产动物中使用InDel标记的研究仅占少数[18-19,24]。因此,本研究利用全基因组重测序技术挖掘卵形鲳鲹与耐低氧性状显著关联的InDel标记,所筛选到的InDel标记可作为分子标记辅助育种,并有助于鉴定与卵形鲳鲹耐缺氧性状相关的候选基因,为进一步研究卵形鲳鲹耐缺氧性状的分子机制提供依据。

    在海南省陵水县中国水产科学研究院南海水产研究所热带水产研究开发中心随机抽取90对卵形鲳鲹亲本。人工催产孵化后,对幼鱼进行标准化养殖管理。随机选取500尾3月龄的幼鱼 [体质量 (22±2.9) g] 运至车间,以流水的养殖方式暂养2周,水体溶解氧水平维持在 (8.8±0.2) mg·L−1,水体pH为8.1。在暂养期间,每天投喂2次商品配合饲料,饲料投喂量为鱼体体质量的5%,实验开始前24 h停止投喂。

    将500尾卵形鲳鲹转移到体积为3 m3的养殖水桶中。在实验过程中,溶解氧含量由溶解氧测定仪 (WTW multi3620,德国) 监测。实验开始时溶解氧质量浓度为8.8 mg·L−1,实验开始后,关闭流水和氧气,以向水中注入氮气的方式降低溶解氧水平,1 h内将溶解氧水平降至窒息点 (1.3 mg·L−1),以1.3 mg·L−1水平一直维持到实验结束。当卵形鲳鲹在第一次触底后不能继续游超过10 s时,定义为死亡。快速捞起死鱼,记录死亡时间。收集鳍条保存在乙醇中,每个鳍条样品与死亡时间相对应。按照卵形鲳鲹的死亡时间分为2组,死亡顺序的前10%为缺氧敏感组,后10%为缺氧耐受组。对收集的共100尾卵形鲳鲹进行全基因组重测序。

    收集的鳍条组织使用MagPure tissue DNA试剂盒 (Magen,广州)并按照说明书步骤提取DNA。用1.0%琼脂糖凝胶电泳 (Biowest agrose, 西班牙) 和NanoDrop 2000 (Thermo Scientific, 美国) 分别测定提取的DNA的质量和浓度。检测合格的DNA送到天津诺禾致源生物公司进行全基因组重测序。

    基因组DNA被随机剪切成350 bp的片段。经过末端修复、磷酸化和添加poly-A尾构建DNA文库,DNA文库在Illumina HiSeq PE150平台进行重测序。100尾卵形鲳鲹重测序获得的原始数据使用软件fastq 0.23[25]进行质量过滤。过滤条件为:1) 删除包含接头 (adapter) 的reads;2) 删除N含量超过5%的reads;3) 删除低质量reads (Qphred≤20) 的碱基数占整个reads长度的50%以上的reads。

    利用BWA 0.7.17软件将过滤后的100尾卵形鲳鲹的reads与参考基因组进行比对定位,并输出BAM文件。使用gatk 4.2.6.0软件[26]检测每个样品的InDel并统计数目,检测得到包含每个个体InDel的VCF文件和gff3格式的基因组注释文件。通过snpEff 5.1和ANNOVAR软件[27-28]获得InDel在基因组的位置以及变异位点为同义突变或者为非同义突变的信息。

    比较低氧敏感组和耐受组两组的差异InDel,筛选低氧耐受组所特有的InDel位点,并定位到变异的基因。通过ClusterProfiler对变异基因进行GO、KEGG富集分析,筛选与耐低氧相关的通路。通过KEGG注释筛选耐低氧相关的候选基因,进行基因功能分析。

    利用PLINK 1.90b6.21软件[29]对100尾卵形鲳鲹InDel数据进行质控,去除次要等位基因频率小于5%和个体变异缺失率大于5%的位点,并进行主成分分析 (Principal component analysis, PCA),并使用ggplot2软件包可视化PCA图。使用GCTA 1.26.0软件[30]评估个体间的亲属关系系数。最后,利用GEMMA 0.98.3软件[31]的线性混合模型 (Linear mixed model, LMM) 将InDel位点和表型数据结合起来进行全基因组关联分析。计算模型为:

    $$ Y=Wα+Xβ+Zμ+ε $$

    式中:Y为个体在缺氧胁迫下的生存时间;W为协方差矩阵;α为群体均值和显著主成分的协变量;X为固定效应矩阵;β为每个InDel位点的效应值;Z为基于InDel的亲缘关系矩阵;μ为加性遗传效应;ε为残差向量。关联分析结果利用ggplot2包绘制曼哈顿图和QQ (Quantile-Quantile) 图。

    在显著相关InDel标记上下游50 kb的区域内筛选候选基因。利用BLAST在非冗余核酸数据库中进行序列比对,确定基因。根据基因的位置和功能,初步确定候选基因。

    通过使用Illumina平台对100尾卵形鲳鲹进行PE150测序,生成了693.48 Gb原始序列数据,并保存在NCBI数据库 (PRJNA658481)。统计100尾卵形鲳鲹的测序数据,测序读段在每个位置的GC含量约为41%,且在整个测序过程基本稳定不变,呈水平线,无偏好性。前几位碱基在核苷酸组成上有一定偏好性,产生波动是由于测序读段连接的接头所引起 (图1-a)。碱基的测序质量Q30平均为90.8% (Q30即碱基正确识别率为99.9%),表明碱基测序的出错率很低,保证了后续分析的准确性 (图1-b)。根据PE150测序技术的特点,测序片段末端的碱基质量一般会比前端低。

    图  1  卵形鲳鲹基因组 GC含量和测序质量分布图
    Fig. 1  GC content and sequence quality in T.ovatus

    通过GATK 4.2.6.0分析共检测到2 574 178个InDel位点。其中插入和缺失的个数分别为1103610和1 470 568个,缺失个数多于插入个数。大部分插入或缺失的片段大小为1~5 bp,分别占各自总数的75.8%和76.6%。1个碱基的插入和缺失的数量最多,分别为496 002和532 256个 (表1)。通过对InDel在基因组上的位置分析,发现有34.51%和35.64%的InDel位于基因间和内含子上,占据了InDel总数的一半以上,只有极少数 (0.68%) 的InDel位于外显子上 (图2)。本研究中采用了重测序的方法对InDel筛选,获得的InDel平均密度为372.1 InDel·Mb−1,24条染色体中,密度最大的是17号染色体 (4 886.6 InDel·Mb−1),密度最小的是4号染色体 (3 400.9 InDel·Mb−1),表明17号染色体发生最多的变异 (图3)。

    表  1  鉴定到的不同长度InDel的数量
    Table  1  Number of InDels of different sizes identified
    InDel类型
    InDel Type
    InDel大小
    InDel Size/bp
    数量
    Number
    百分比
    Percentage/%
    InDel类型
    InDel Type
    InDel大小
    InDel Size/bp
    数量
    Number
    百分比
    Percentage/%
    插入 Insertion 1 496002 44.94 缺失 Deletion 1 532256 37.7
    2 189179 17.14 2 292212 20.7
    3 58625 5.31 3 107471 7.6
    4 69620 6.31 4 106766 7.6
    5 23055 2.09 5 43919 3.1
    6 43192 3.91 6 62130 4.4
    7 13832 1.25 7 27167 1.9
    8 28729 2.60 8 39844 2.8
    9 11719 1.06 9 20564 1.5
    10 20236 1.83 10 27837 2.0
    ≥11 149421 13.54 ≥11 152545 10.8
    总计 Total    1103610 总计 Total    1412711
    下载: 导出CSV 
    | 显示表格
    图  2  InDel在基因组上的分布位置
    Fig. 2  Region of InDel on genome
    图  3  InDel在基因组上的分布密度
    Fig. 3  Distribution density of InDel on genome

    分别对敏感组和耐受组进行InDel筛选,敏感组和耐受组分别筛选到1 723 005和1 720 945个InDel,其中耐受组所特有的InDel 249 395个。使用snpEff 5.1和ANNOVAR软件对耐受组所特有的InDel进行基因分析,249 395个InDel中有2209个位于外显子上,涉及543个基因。通过GO注释分析,多数变异基因主要富集在染色体和顶体泡的组成、核酸结合和GTP结合的分子功能和细胞对DNA损伤刺激的反应,以及Wnt信号通路上 (表2)。通过KEGG富集分析,这些基因主要富集在核苷酸切除修复信号通路和细胞黏着分子上 (图4)。

    表  2  变异基因GO功能分类注释
    Table  2  GO classification of mutated genes
    类别
    Classification
    条目
    Term
    变异基因数量
    Number of
    variant genes
    细胞组分
    Cellular component
    染色体 7
    顶体泡 6
    线粒体外膜 4
    电压门控钾离子通道复合物 3
    肌间盘 3
    分子功能
    Molecular function
    核酸结合 11
    GTP结合 8
    电压门控离子通道活性 5
    镁离子结合 5
    转录辅抑制活性 5
    生物过程
    Biological process
    细胞对DNA损伤刺激的反应 7
    Wnt信号通路 6
    细胞器组织 6
    骨骼肌组织发育 5
    DNA整合 5
    下载: 导出CSV 
    | 显示表格
    图  4  KEGG富集分析
    Fig. 4  KEGG enrichment scatter plot

    利用InDel位点对实验群体进行主成分分析和亲缘关系分析。亲缘关系分析发现,大多数个体间的亲缘关系介于−0.1~0.3 (图5-a)。在主成分分析中,主成分因子1和2将实验样本分成5组 (图5-b)。采用GEMMA 0.98.3软件的线性混合模型对耐低氧性状进行全基因组关联分析,观察到的InDel位点的P值大于其预期值 (图5-c),表明线性混合模型处理样本数据是合理的,所得结果是可靠的。根据Bonferroni校正,阈值设置为P= (0.05/N),其中N为突变位点的数量。在本研究中,通过PLINK 1.90b6.21软件质控后,共有336 391个InDel位点被用于卵形鲳鲹耐低氧性状的关联分析。因此,本研究中的阈值设置为−lg (0.05/336 391)=6.83。由于Bonferroni校正过于严格会导致假阴性,所以需要根据曼哈顿图调整阈值[32-33]。阈值设置为−lg (1/336 391)=5.53作为建议显著性阈值,且达到该阈值的InDel位点可能与耐低氧性状关联。336 391个InDel位点中无InDel位点达到显著性阈值,包括调整后的阈值。但3个InDel位点的P值接近建议性显著性阈值。其中2个InDel位点位于3号染色体,根据InDel在染色体上的位置对其编号为InDel22883061和InDel24919481,对应的−lg P分别为5.26和5.11。另外一个InDel位于19号染色体,将其命名为InDel14451779,对应的−lg P 为5.25 (图5-d)。

    图  5  100尾卵形鲳鲹亲缘关系热图 (a)、主成分分析 (b)、耐低氧性状全基因组关联分析的QQ图 (c) 和曼哈顿图 (d)
    Fig. 5  Heat map of genetic relationship (a), principal component analysis (b), Q-Q plot (c) and Manhattan plot (d) of genome-wide association analysis for low oxygen tolerance of 100 individuals of T. ovatus

    对重要的3个InDel位点上下游50 kb的基因组序列进行注释后,探索可能与卵形鲳鲹耐缺氧相关的基因。3个位点共注释到9个耐低氧性状相关的候选基因。位于3号染色体上的2个InDel位点共注释到3个基因 (gpr153、acot7和lrfn2b)。另一位于19号染色体的InDel位点共注释到6个基因 (tmem237bmpp4、als2、LOC111232287、LOC111232276和LOC120797970)。

    InDel是广泛分布于基因组的结构变异的主要来源之一,是对其他基于序列的遗传标记 (如SSRs和SNPs) 的有价值补充,被公认为遗传分析的有效标记系统,且InDel标记具有密度大、可快速和经济有效地进行基因分型等优点[34]。以往有关鱼类耐低氧性状分子标记筛选研究大多集中在SNP上,很少有关于筛选InDel标记的报道[32-33]。本研究利用重测序技术对100尾卵形鲳鲹进行测序,获得了693.48 Gb的测序数据,其中Q30 (正确识别率大于99.9%的碱基占总体碱基的百分比) 平均值为90.8%,测序数据的质量主要分布在Q30 (≥80%),这样能够保证后续分析的正常进行。在此基础上,计算了InDel位点在基因组上的平均分布密度为3 972.1 InDel·Mb−1,高密度的InDel位点为挖掘与耐低氧性状关联的位点提供了基础。通过对InDel在基因组位置的注释发现,大部分InDel位点位于基因间和内含子上 (70.5%),只有极少的InDel位点位于外显子上 (0.7%),证明编码区的InDel位点数量较少,这与之前SNP在基因组上的分布情况一致[35]。这种分布模式不仅存在于卵形鲳鲹中,在狮头鹅 (Shitou goose)、陆地棉 (Gossypium hirsutum) 等生物中也表现出相同的分布模式[36-37],位于非编码区的InDel位点多于位于编码区的数量,这一结果符合生物学特性。

    通过对敏感组和耐受组的InDel位点进行分析,获得了249 395个耐受组特有的InDel位点,其中2 209个位于外显子上,涉及543个基因。在此基础上对这些基因进行功能注释,发现核苷酸切除修复信号通路和细胞黏着分子中均有InDel变异。核苷酸切除修复是4条DNA损伤修复的基本途径之一。本研究中发现rpa1、ercc3和pold3这3个基因富集到这条通路上,rpa1基因编码异三聚体复制蛋白A复合物的最大亚基,它与单链DNA结合,形成一个核蛋白复合物,在DNA代谢中发挥重要作用,参与DNA复制、修复[38]pold3基因编码DNA聚合酶delta的66-kD亚基,DNA聚合酶delta具有聚合酶和3'—5'外切酶活性,在DNA复制和修复中起着关键作用[39]。位于这些在DNA复制和修复过程中发挥重要作用的基因上的InDel位点,可能影响了卵形鲳鲹的耐低氧能力。另一较多基因富集的信号通路细胞黏着分子是一种表达于细胞表面的蛋白,参与细胞的活化和信号转导等生物过程[40],这些信号中的变异基因可能影响了低氧环境下卵形鲳鲹信号转导的能力。

    采用与敏感组比较、对耐受组所特有的位于外显子上的InDel进行基因注释的方法,获得的位点信息多,注释的基因数量大,仅依靠基因通路和基因功能等信息难以获得可靠的InDel位点用于分子标记的开发。InDel虽然在基因组上分布广泛,但在鱼类分子标记的开发中主要集中在SNP上,对InDel的关注很少,但也有利用InDel位点对斑节对虾 (Penaeus monodon) 早期性别鉴定的应用[41]。虽然在鱼类中仍无利用线性混合模型筛选与性状显著关联的InDel的研究,但已有在芝麻 (Sesamum indicum)、山羊 (Capra hircus) 等动植物中[42-43]的研究,利用全基因组关联分析 (Genome wide association study, GWAS) 的线性混合模型将InDel位点与表型信息关联,挖掘与表型信息显著关联的InDel位点,Yang等[44]也利用全基因组关联分析的方法筛选到了南丹瑶鸡与鸡冠性状显著相关的InDel位点。证明利用线性混合模型将InDel位点与性状关联的可行性。本研究为了筛选与耐低氧性状显著相关的InDel位点,首先利用InDel位点进行了主成分分析和亲缘关系分析。通过主成分分析,主成分因子1和2大致将实验群体分为5个群体,亲缘关系分析显示实验群体的亲缘关系大部分介于−0.1~0.3,这与San等[35]利用SNP分析的结果一致,保证了后续关联分析的准确性。在关联分析的结果中,InDel位点的P观察值和期望值相同,说明分析模型是合理的。但所有的P观测值均未明显超过期望值,说明分析结果未找到与耐低氧性状显著关联的位点,这可能是因为耐低氧性状是由微效多基因控制的,单个变异位点效应对性状的影响弱,不能达到显著性阈值;也可能是因为群体样本量小,本研究采用重测序的方法对InDel分型,成本虽然高但可以挖掘稀有变异,所以只对极端表型值的100尾个体进行了重测序,样本量偏小影响了检验效能。关联分析结果发现,没有位点达到显著性阈值,但有3个InDel位点接近显著性阈值,仍具有开发分子标记的价值。通过对显著性位点上下游50 kb的基因序列注释,共注释到9个候选基因,其中Acot7对中、长链脂酰辅酶a具有较高的酶活性。对小鼠 (Mus musculus) 的研究发现,Acot7具有保护神经元免受脂肪酸毒性的作用[45]Lrfn2b编码一个突触黏附样分子,在小鼠中Lrfn2b可能通过对n-甲基-d-天冬氨酸 (N-methyl-D- aspartate, NMDA) 受体的调节作用来调控红细胞生成[46]。由此推测这2个基因可能在卵形鲳鲹抵抗低氧环境的过程中发挥作用。Mpp4编码的蛋白是鸟苷酸激酶家族中的一种视网膜特异性支架蛋白,参与组织光感受器带状突触中的突触前蛋白复合物[47]Gpr135编码一个完整的膜蛋白,属于G蛋白偶联受体的A类视紫红质超家族,在大鼠 (Rattus norvegicus) 中敲除后其食物摄入量受到影响[48]Als2编码的蛋白可以激活GTPases的Ras超家族成员,研究发现als2基因被敲除的斑马鱼 (Danio rerio) 出现了严重的发育异常、游泳障碍和运动神经元紊乱[49]Tmem237b编码一种四跨膜蛋白,影响斑马鱼的胚胎发育[50]。这4个候选基因在低氧胁迫下发挥的功能仍需进一步研究。

    综上所述,卵形鲳鲹在低氧环境下,核苷酸切除修复信号通路和细胞黏着分子可能发挥了重要作用。通过InDel位点与耐低氧性状的关联分析,挖掘到了3个与选耐低氧能力相关的InDel位点,并对筛选到的InDel位点进行基因注释获得了9个候选基因。筛选到的InDel位点对后期分子标记选择育种的选择和鉴定有重要价值,注释到的候选基因为卵形鲳鲹耐低氧的机理研究提供了依据。

  • 图  1   凡纳滨对虾微卫星位点M1、HLJN-008、TUMXLv8.256、TUYFLvL16.1a (a),TUMXLv7.121、TUMXLv9.43、TUMXLv10.312、TUSWLvSU233 (b) 和Pvan1758、Pvan1815、HLJN-023、TUDGLv1-3.224、TUGAPv7-9.95 (c) 的毛细管电泳基因分型

    Figure  1.   Genotyping of microsatellite loci M1, HLJN-008, TUMXLv8.256, TUYFLvL16.1a (a), TUMXLv7.121, TUMXLv9.43, TUMXLv10.312, TUSWLvSU233 (b) and Pvan1758, Pvan1815, HLJN-023, TUDGLv1-3.224, TUGAPv7-9.95 (c) in L. vannamei by capillary electrophoresis

    图  2   凡纳滨对虾微卫星DNA中高频等位基因 (P>0.25) 在各养殖群体的频率分布情况

    Figure  2.   Frequency distribution of high frequency alleles (P>0.25) from microsatellite DNA in seven cultured populations of L. vannamei

    图  3   凡纳滨对虾7个养殖群体的UPGMA聚类图

    Figure  3.   UPGMA dendrogram of seven cultured populations of L. vannamei

    表  1   凡纳滨对虾7个养殖群体的取样信息

    Table  1   Sampling information of seven cultured populations of L. vannamei

    群体
    population
    亲虾来源
    parental origin
    采样地点
    sampling location
    种苗来源
    sample source
    采样时间
    sampling date
    TH-A1 CP 中山,正大神湾虾苗场育苗车间 正大神湾虾苗场 2017年8月
    TH-A2 CP 茂名,广东冠利海洋生物有限责任公司工厂化养殖池 中海水产种苗科技有限公司 2017年8月
    TH-B SyAqua 茂名,广东冠利海洋生物有限责任公司工厂化养殖池 2017年8月
    US-C1 SIS 茂名,广东冠利海洋生物有限责任公司工厂化养殖池 海威水产养殖有限公司 2017年9月
    US-C2 SIS 茂名,广东冠利海洋生物有限责任公司土塘养殖池 海威水产养殖有限公司 2017年9月
    US-C3 SIS 茂名,广东冠利海洋生物有限责任公司土塘养殖池 海威水产养殖有限公司 2017年9月
    US-C4 SIS 汕尾,红海湾养殖基地工厂化养殖池 海尚种苗培育基地 2017年8月
    下载: 导出CSV

    表  2   凡纳滨对虾微卫星位点扩增所用引物及其荧光标记

    Table  2   Primers and their fluorescence labeling for PCR amplification of L. vannamei microsatellite loci

    位点
    locus
    序列 (5'−3') 和荧光标记类型
    sequence (5'−3') and fluorescence type
    重复序列
    repeat motif
    参考文献
    Reference
    M1 F: TAMRA-CTAACCCAATATCGAATC
    R: GTGTGTTGCCGAATCGAA
    (TTTC)4...(TTTC)5...(TTTC)6 [10]
    Pvan1758 F: TAMRA-TATGCTCGTTCCCTTTGCTT
    R: TTGAAGGAAAAGTGTTGGGG
    (T)10(TC)7...(T)4(TC) [11]
    Pvan1815 F: ROX-GATCATTCGCCCCTCTTTTT
    R: ATCTACGGTTCGAGAGCAGA
    (T)7(CT)2(CTTT)4…(CT)6 [11]
    HLJN-008 F: HEX-CGTAAGTCCTGCAAAAGAAACT
    R: GCCGTTCAACTATATATCAGCA
    (TA)12(GA)(TAGA)4A(AT)(AG)49A(AG)10(AT)(AG)A(AG)10TG(GA)(AG)3A(AT)(AG)6 [12]
    HLJN-023 F: HEX-AAGAGAGATGGAAGGAGTAAGTGC
    R: GATCAATACCTTGCAGCGAAA
    (AG)3T(GA)24A(AG)15AAA(GA)2(AG)(GA)9 [12]
    TUMXLv7.121 F: 6-FAM-GGCACACTGTTTAGTCCTCG
    R: CGAACAGAATGGCAGAGGAG
    (GTT)3…(GA)3…(TC)3…(GT)3…(TC)3…(TC)3…(TC)3 [13]
    TUMXLv8.256 F: 6-FAM-GGACTCACACTTCTGGTTC
    R: GGCTGCACCTTGTAAGTC
    (AAT)4 [13]
    TUMXLv9.43 F: ROX-GAGAGCAAATAAGAAAGGGC
    R: AGGATGCAAATGATAACGAG
    (TG)4C(GT)4TC(GT)6(GC)3…(GT)3TTA(TG)4…(GGA)3(GA)3…(CT)3CC(CT)3CC(CCCT)6(CT)10…(TC)6…(TC)3 [13]
    TUMXLv10.312 F: HEX-ATACGAAACACCCCATCCC
    R: GTGGTCTTACCTCGTGGCTC
    (AG)30…(AG)5 [13]
    TUSWLvSU233 F: TAMRA-CCCGACTTGGCTTTTAGTTG
    R: GAGATTGCTATCCTCGGCTG
    (TGA)3ATG(GA)3…(GA)4AGTA(AGC)3 [14-15]
    TUYFLvL16.1a F: ROX-ATGGCACAAATAGGATCTTG
    R: GACTGGAAGAGCACTGATTC
    (AAT)4…(AC)3 [14-15]
    TUDGLv1-3.224 F: 6-FAM-ACTAGTGGATCTGTCTATTCAT
    R: ATACCCACCCATGCATGTTA
    (TAGA)3…(TAGA)3…(ACAG)4(AG)21A(AG)30 [16]
    TUGAPv7-9.95 F: ROX-GATCCTGCGAGTCACTTTATCTC
    R: TTTATTGCGTATCCCAGAAGC
    (TC)3…(AT)3…(AG)3…(TC)13…(AT)24G(TA)5TG(TA)4TG(TA)4…(AT)3…(AT)24 [16]
    下载: 导出CSV

    表  3   凡纳滨对虾7个养殖群体在12个微卫星位点的多态性信息

    Table  3   Polymorphic information at 12 microsatellite loci in seven cultured populations of L. vannamei

    群体
    population
    参数
    parameter
    位点 locus
    M1 Pvan1758 Pvan1815 HLJN-008 HLJN-023 TUMX Lv7.121 TUMX Lv8.256 TUMX Lv9.43 TUMX Lv10.312 TUYF LvL16.1a TUDG Lv1-3.224 TUGA Pv7-9.95 均值
    average
    TH-A1 Na 5 6 4 7 8 4 3 8 4 3 9 5 5.500
    He 0.719 0.730 0.624 0.817 0.799 0.504 0.385 0.769 0.496 0.254 0.774 0.678 0.629
    Ho 0.167 0.526 0.667 0.136 0.625 0.333 0.083 0.400 0.458 0.125 0.625 0.417 0.380
    PIC 0.675 0.686 0.572 0.794 0.771 0.458 0.325 0.743 0.442 0.231 0.754 0.631 0.590
    PHWE 0.000 0.056 0.260 0.000 0.000 0.119 0.000 0.000 0.269 0.011 0.001 0.000 0.060
    TH-A2 Na 7 7 5 5 7 4 3 11 3 3 7 3 5.417
    He 0.746 0.774 0.731 0.646 0.787 0.626 0.443 0.865 0.570 0.536 0.804 0.518 0.670
    Ho 0.304 0.375 0.500 0.417 0.625 0.667 0.391 0.625 0.583 0.250 0.875 0.292 0.492
    PIC 0.714 0.742 0.684 0.608 0.755 0.567 0.399 0.850 0.496 0.476 0.776 0.404 0.622
    PHWE 0.000 0.000 0.153 0.003 0.018 0.800 0.155 0.000 0.045 0.001 0.101 0.083 0.113
    TH-B Na 3 5 2 5 4 3 3 5 3 1 4 2 3.333
    He 0.508 0.719 0.080 0.559 0.555 0.556 0.503 0.583 0.633 0.000 0.532 0.500 0.477
    Ho 0.125 0.435 0.083 0.417 0.583 0.708 0.375 0.250 0.625 0.000 0.500 0.333 0.370
    PIC 0.428 0.673 0.077 0.461 0.482 0.496 0.435 0.503 0.556 0.000 0.488 0.375 0.414
    PHWE 0.000 0.000 0.882 0.000 0.868 0.090 0.033 0.002 0.166 0.062 0.082 0.199
    US-C1 Na 8 7 6 7 7 4 3 11 3 3 11 4 6.167
    He 0.837 0.768 0.743 0.796 0.790 0.539 0.348 0.861 0.555 0.344 0.797 0.463 0.653
    Ho 0.217 0.478 0.750 0.739 0.542 0.625 0.250 0.583 0.542 0.167 0.750 0.417 0.505
    PIC 0.816 0.732 0.703 0.767 0.759 0.501 0.316 0.846 0.456 0.307 0.772 0.401 0.615
    PHWE 0.000 0.001 0.835 0.294 0.011 0.908 0.061 0.204 0.939 0.004 0.042 0.643 0.329
    US-C2 Na 8 5 6 6 7 4 3 7 3 3 8 4 5.333
    He 0.785 0.747 0.700 0.819 0.816 0.679 0.430 0.835 0.570 0.292 0.760 0.541 0.664
    Ho 0.217 0.250 0.583 0.542 0.625 0.667 0.333 0.667 0.583 0.083 0.667 0.292 0.459
    PIC 0.760 0.704 0.657 0.794 0.791 0.622 0.388 0.813 0.496 0.272 0.727 0.450 0.623
    PHWE 0.000 0.000 0.317 0.006 0.001 0.750 0.159 0.000 0.617 0.000 0.337 0.000 0.182
    US-C3 Na 5 5 3 6 8 4 3 6 4 2 9 3 4.833
    He 0.674 0.739 0.284 0.735 0.790 0.669 0.260 0.695 0.598 0.041 0.701 0.617 0.567
    Ho 0.083 0.667 0.333 0.500 0.917 0.667 0.125 0.583 0.667 0.042 0.708 0.542 0.486
    PIC 0.632 0.703 0.254 0.689 0.760 0.601 0.244 0.666 0.519 0.040 0.676 0.544 0.527
    PHWE 0.000 0.000 0.843 0.007 0.589 0.895 0.000 0.551 0.551 1.000 0.825 0.019 0.440
    US-C4 Na 4 7 3 5 8 3 2 4 3 1 7 3 4.167
    He 0.738 0.785 0.403 0.733 0.773 0.562 0.153 0.600 0.594 0.000 0.745 0.539 0.552
    Ho 0.182 0.333 0.292 0.458 0.565 0.583 0.167 0.625 0.708 0.000 0.833 0.458 0.434
    PIC 0.689 0.757 0.363 0.694 0.740 0.482 0.141 0.533 0.511 0.000 0.708 0.432 0.504
    PHWE 0.000 0.000 0.146 0.000 0.000 0.788 0.703 0.942 0.681 / 0.669 0.380 0.392
    下载: 导出CSV

    表  4   凡纳滨对虾7个养殖群体的Nei's遗传相似度 (对角线上方) 和遗传距离 (对角线下方)

    Table  4   Nei's genetic identity (above diagonal) and genetic distance (below diagonal) among seven cultured populations of L. vannamei

    种群
    population
    TH-A1 TH-A2 TH-B US-C1 US-C2 US-C3 US-C4
    TH-A1 0.669 4 0.508 5 0.731 5 0.709 6 0.652 0 0.562 2
    TH-A2 0.401 4 0.703 1 0.921 0 0.921 9 0.758 3 0.669 2
    TH-B 0.676 4 0.352 2 0.695 0 0.769 6 0.858 7 0.767 5
    US-C1 0.312 7 0.082 3 0.363 9 0.944 3 0.788 4 0.661 8
    US-C2 0.343 1 0.081 3 0.261 9 0.057 3 0.847 0 0.738 7
    US-C3 0.427 7 0.276 7 0.152 3 0.237 8 0.166 1 0.883 1
    US-C4 0.575 9 0.401 7 0.264 7 0.412 8 0.302 8 0.124 3
    下载: 导出CSV
  • [1] 王兴强, 马甡, 董双林. 凡纳滨对虾生物学及养殖生态学研究进展[J]. 海洋湖沼通报, 2004(4): 94-100. doi: 10.3969/j.issn.1003-6482.2004.04.016
    [2] 农业部渔业渔政管理局. 2017中国渔业统计年鉴[M]. 北京: 中国农业出版社, 2017: 22-61.
    [3] 童馨, 龚世圆, 喻达辉, 等. 凡纳滨对虾(Litopenaeus vannanamei)不同世代养殖群体的遗传多样性分析[J]. 海洋与湖沼, 2009, 40(2): 214-220. doi: 10.3321/j.issn:0029-814X.2009.02.017
    [4] 颉晓勇, 苏天凤, 陈文. 凡纳滨对虾6个养殖群体遗传多样性的比较分析[J]. 南方水产, 2008, 4(6): 42-49. doi: 10.3969/j.issn.2095-0780.2008.06.006
    [5] 代平, 孔杰, 栾生. 我国凡纳滨对虾种质资源引进与分析[J]. 科学养鱼, 2018(1): 3-5.
    [6]

    ZHANG K, WANG W J, LI W Y, et al. Analysis of genetic diversity and differentiation of seven stocks of Litopenaeus vannamei using microsatellite markers[J]. J Ocean Univ China, 2014, 13(4): 647-656. doi: 10.1007/s11802-014-2208-2

    [7] 马春艳, 马洪雨, 马凌波, 等. 凡纳滨对虾引进群体和2个养殖群体遗传变异的微卫星分析[J]. 海洋渔业, 2011, 33(1): 1-8. doi: 10.3969/j.issn.1004-2490.2011.01.001
    [8] 包秀凤. 凡纳滨对虾选育群体遗传多样性分析[D]. 湛江: 广东海洋大学, 2014: 11-26.
    [9]

    HU X J, CAO Y C, WEN G L, et al. Effect of combined use of Bacillus and molasses on microbial communities in shrimp cultural enclosure systems[J]. Aquacult Res, 2017, 48(6): 2691-2705. doi: 10.1111/are.2017.48.issue-6

    [10]

    GARCIA D K, DHAR A K, ALCIVAR-WARREN A. Molecular analysis of a RAPD marker (B20) reveals two microsatellites and differential mRNA expression in Penaeus vannamei[J]. Mol Mar Biol Biotechnol, 1996, 5(1): 71-83. doi: 10.1007/BF02762417

    [11]

    CRUZ P, MEJIA-RUIZ C H, PEREZ-ENRIQUEZ R, et al. Isolation and characterization of microsatellites in Pacific white shrimp Penaeus (Litopenaeus) vannamei[J]. Mol Ecol Resour, 2002, 2(3): 239-241.

    [12]

    JIA Z, SUN X, LIANG L, et al. Isolation and characterization of microsatellite markers from Pacific white shrimp (Litopenaeus vannamei)[J]. Mol Ecol Resour, 2006, 6(4): 1282-1284.

    [13]

    MEEHAN D, XU Z, ZUNIGA G, et al. High frequency and large number of polymorphic microsatellites in cultured shrimp, Penaeus (Litopenaeus) vannamei[J]. Mar Biotechnol, 2003, 5(4): 311-330. doi: 10.1007/s10126-002-0092-z

    [14]

    ALCIVAR-WARREN A, MEEHAN-MEOLA D, PARK S W, et al. ShrimpMap: a low-density, microsatellite-based linkage map of the pacific whiteleg shrimp, Litopenaeus vannamei: identification of sex-linked markers in linkage group 4[J]. J Shellfish Res, 2007, 26(4): 1259-1277. doi: 10.2983/0730-8000(2007)26[1259:SALMLM]2.0.CO;2

    [15]

    ALCIVAR-WARREN A, SONG L, MEEHAN D, et al. Mapping simple sequence repeat markers identified in ESTs from a subtracted cDNA library of white spot virus-challenged shrimp Litopenaeus vannamei[J]. J Shellfish Res, 2007, 26(4): 1247-1258. doi: 10.2983/0730-8000(2007)26[1247:CAMOES]2.0.CO;2

    [16]

    GARCIA D K, ALCIVAR-WARREN A. Characterization of 35 new microsatellite genetic markers for the pacific whiteleg shrimp, Litopenaeus vannamei: their usefulness for studying genetic diversity of wild and cultured stocks, tracing pedigree in breeding programs, and linkage mapping[J]. J Shellfish Res, 2007, 26(4): 1203-1216. doi: 10.2983/0730-8000(2007)26[1203:CONMGM]2.0.CO;2

    [17]

    YEH F C, BOYLE T J B. Population genetic analysis of co-dominant and dominant markers and quantitative traits[J]. Belg J Botany, 1997, 129(2): 157.

    [18]

    NEI M. Estimation of average heterozygosity and genetic distance from a small number of individuals[J]. Genetics, 1978, 89(3): 583-590.

    [19]

    NAGY S, POCZAI P, CERNÁK I, et al. PICcalc: An online program to calculate polymorphic information content for molecular genetic studies[J]. Biochem Genet, 2012, 50(9/10): 670-672.

    [20]

    TAMURA K, PETERSON D, PETERSON N, et al. MEGA5: 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

    [21] 孙效文, 张晓锋, 赵莹莹, 等. 水产生物微卫星标记技术研究进展及其应用[J]. 中国水产科学, 2008, 15(4): 689-703. doi: 10.3321/j.issn:1005-8737.2008.04.022
    [22] 杨铭, 于洋, 张晓军, 等. 基于转录组数据的凡纳滨对虾微卫星标记开发[J]. 海洋科学, 2017, 41(2): 96-102.
    [23]

    YU Y, ZHANG X J, YUAN J B, et al. Genome survey and high-density genetic map construction provide genomic and genetic resources for the Pacific white shrimp Litopenaeus vannamei[J]. Sci Rep, 2015, 5: 15612. doi: 10.1038/srep15612

    [24] 谢丽, 陈国良, 叶富良, 等. 凡纳滨对虾4个选育群体遗传多样性的SSR分析[J]. 广东海洋大学学报, 2009, 29(4): 5-9. doi: 10.3969/j.issn.1673-9159.2009.04.002
    [25]

    LEBERG P L. Estimating allelic richness: effects of sample size and bottlenecks[J]. Mol Ecol, 2002, 11(11): 2445-2449.

    [26]

    CHANG Y M, LIANG L Q, MA H T, et al. Microsatellite analysis of genetic diversity and population structure of Chinese mitten crab (Eriocheir sinensis)[J]. J Genet Genomics, 2008, 35(3): 171-176. doi: 10.1016/S1673-8527(08)60023-5

    [27]

    QUAN Y C, SUN X W, LIANG L Q. Genetic popymorphism of microsatellite DNA in two populations of nothern sheatfish (Silurus soldatovi)[J]. Acta Genet Sin, 2006, 33: 908-916. doi: 10.1016/S0379-4172(06)60125-X

    [28] 孙成波, 陈国良, 童汉荣, 等. 美国4个凡纳滨对虾(Litopenaeus vannamei)种群形态差异与判别分析[J]. 海洋与湖沼, 2009, 40(1): 27-32. doi: 10.3321/j.issn:0029-814X.2009.01.005
  • 期刊类型引用(4)

    1. 海萨·艾也力汗,张钰,杨博文,咸玉兰,高攀,沈玉帮. 基于简化基因组测序筛选白斑狗鱼耐热性状关联的InDel标记. 水产学报. 2024(03): 58-65 . 百度学术
    2. 张鑫淼,赵春龙,谢鹏,冯广朋,陈朋,王立新,熊冬梅. 基于全基因组InDel分析我国细鳞鲑遗传分化及本地适应性. 水生生物学报. 2024(10): 1724-1735 . 百度学术
    3. 鲁翠云,罗鸣,郑先虎,刘龙龙,张国庆,刘天奇,陈有铭. 用微卫星标记分析卵形鲳鲹两个繁育群体的遗传结构. 水产学杂志. 2024(05): 28-36+44 . 百度学术
    4. 范嗣刚,黄皓,王鹏飞,赵超,闫路路,邱丽华. 基于微卫星标记的花鲈亲子鉴定技术. 广东海洋大学学报. 2023(05): 26-33 . 百度学术

    其他类型引用(7)

推荐阅读
低盐水体so4 2−/cl− 胁迫下凡纳滨对虾生长、肝胰腺与鳃组织结构及酶活力比较
贺铮 et al., 南方水产科学, 2025
Lh原油和0# 柴油乳化液对凡纳滨对虾肝胰腺抗氧化酶活性及相关功能基因表达的影响
沈楚焰 et al., 南方水产科学, 2025
Hctlr1通过myd88-nf-κb信号通路参与三角帆蚌抗菌免疫应答
路俊怡 et al., 南方水产科学, 2024
黄颡鱼源迟缓爱德华氏菌的分离鉴定及其药物敏感性研究
高宣 et al., 南方水产科学, 2024
基于简化基因组测序筛选白斑狗鱼耐热性状关联的indel标记
海萨·艾也力汗 et al., 水产学报, 2024
低压低氧环境暴露下睾丸单细胞转录组分析揭示生殖毒性机制
吉运华 et al., 协和医学, 2024
Sentinel node mapping in high-intermediate and high-risk endometrial cancer: analysis of 5-year oncologic outcomes
Cuccu, Ilaria et al., EJSO, 2024
Sub-2 nm irrunimoco high-entropy alloy with iridium-rich medium-entropy oxide shell to boost acidic oxygen evolution
Yao, Longping et al., ADVANCED MATERIALS, 2024
A novel method to identify influential nodes based on hybrid topology structure
PHYSICAL COMMUNICATION
Lightweight marine biodetection model based on improved yolov10
ALEXANDRIA ENGINEERING JOURNAL
Powered by
图(3)  /  表(4)
计量
  • 文章访问数:  3923
  • HTML全文浏览量:  1790
  • PDF下载量:  49
  • 被引次数: 11
出版历程
  • 收稿日期:  2018-06-19
  • 修回日期:  2018-09-04
  • 录用日期:  2018-10-09
  • 网络出版日期:  2018-12-04
  • 刊出日期:  2019-02-04

目录

/

返回文章
返回