基于卵形鲳鲹基因组重测序的InDel标记挖掘与耐低氧性状关联分析

伞利择, 刘宝锁, 张楠, 郭梁, 郭华阳, 朱克诚, 张殿昌

伞利择, 刘宝锁, 张楠, 郭梁, 郭华阳, 朱克诚, 张殿昌. 基于卵形鲳鲹基因组重测序的InDel标记挖掘与耐低氧性状关联分析[J]. 南方水产科学, 2022, 18(5): 100-109. DOI: 10.12131/20210347
引用本文: 伞利择, 刘宝锁, 张楠, 郭梁, 郭华阳, 朱克诚, 张殿昌. 基于卵形鲳鲹基因组重测序的InDel标记挖掘与耐低氧性状关联分析[J]. 南方水产科学, 2022, 18(5): 100-109. DOI: 10.12131/20210347
SAN Lize, LIU Baosuo, ZHANG Nan, GUO Liang, GUO Huayang, ZHU Kecheng, ZHANG Dianchang. Mining of InDel marker and association analysis of hypoxia tolerance traits in Trachinotus ovatus based on resequencing[J]. South China Fisheries Science, 2022, 18(5): 100-109. DOI: 10.12131/20210347
Citation: SAN Lize, LIU Baosuo, ZHANG Nan, GUO Liang, GUO Huayang, ZHU Kecheng, ZHANG Dianchang. Mining of InDel marker and association analysis of hypoxia tolerance traits in Trachinotus ovatus based on resequencing[J]. South China Fisheries Science, 2022, 18(5): 100-109. DOI: 10.12131/20210347

基于卵形鲳鲹基因组重测序的InDel标记挖掘与耐低氧性状关联分析

基金项目: 国家自然科学基金项目 (U20A2064);财政部和农业农村部: 国家海水鱼产业技术体系资助 (CARS-47);国家重点研发计划项目 (2018YFD0900301);中国水产科学研究院中央级公益性科研院所基本科研业务费专项资金资助 (2020TD29)
详细信息
    作者简介:

    伞利择 (1996—),男,硕士研究生,研究方向为水产动物遗传育种。E-mail: cdsanlz@163.com

    通讯作者:

    张殿昌 (1977—),男,研究员,博士,从事水产种质资源与遗传育种研究。E-mail: zhangdch@scsfri.ac.cn

  • 中图分类号: S 917.4

Mining of InDel marker and association analysis of hypoxia tolerance traits in Trachinotus ovatus based on resequencing

  • 摘要: 卵形鲳鲹 (Trachinotus ovatus) 是高耗氧率的鱼类,低氧极易导致其死亡。筛选其耐低氧性状InDel分子标记并挖掘影响耐低氧性状的功能基因,可为选育具有较强耐低氧能力的卵形鲳鲹提供有益指导。利用全基因组重测序技术分析了卵形鲳鲹基因InDel差异,挖掘卵形鲳鲹与耐低氧性状显著关联的InDel位点,并探讨了与耐低氧性状相关的候选基因。结果表明,测序共获得693.48 Gb数据,其中Q30的平均值为90.8%,注释分析发现了共2 574 178个InDel位点。敏感组50尾卵形鲳鲹所特有的InDel共249 395个,其中2 209个位于外显子上,在核苷酸切除修复信号通路和细胞黏着分子中存在变异基因。将InDel位点与耐低氧性状关联分析发现,3个InDel位点 (InDel 22883061、InDel 24919481和InDel 14451779) 接近显著性阈值,注释分析获得了9个候选基因。筛选到的InDel位点对后期分子标记选择育种的选择和鉴定有重要价值,注释到的候选基因为卵形鲳鲹耐低氧的机理研究提供了基础和依据。
    Abstract: Trachinotus ovatus is a fish with high oxygen consumption rate, and hypoxia can easily lead to its death. Screening for InDel molecular markers for hypoxia tolerance traits in T. ovatus and discovering functional genes affecting hypoxia tolerance traits can provide guidances for selecting and breeding T. ovatus with strong hypoxia tolerance. In this study, we applied the whole genome resequencing technology to analyze the InDel differences of T. ovatus gene, excavated the InDel sites that were significantly associated with hypoxia tolerance traits in T. ovatus, and explored the candidate genes related to hypoxia tolerance traits. Altogether 693.48 Gb was obtained by sequencing, and the average value of Q30 was 90.8%. A total of 2 574 178 InDel markers were found by annotation analysis. A total of 249 395 InDel markers were found in the 50 individuals in the tolerance group, of which 2 209 were located in exons. It is found that there were mutated genes in the nucleotide excise repair signaling pathway and cell adhesion molecules. Three InDel markers (InDel 22883061, InDel 24919481 and InDel 14451779) were found to be close to the significance threshold by association analysis between InDel markers and hypoxia tolerance trait. Nine candidate genes were obtained by annotation analysis. The selected InDel markers were of great value for selecting and identifying molecular marker selection breeding. The annotated candidate genes provide a basis and references for the study of the hypoxia tolerance mechanism of T. ovatus.
  • 卵形鲳鲹 (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   卵形鲳鲹基因组 GC含量和测序质量分布图

    Figure  1.   GC content and sequence quality in T.ovatus

    图  2   InDel在基因组上的分布位置

    Figure  2.   Region of InDel on genome

    图  3   InDel在基因组上的分布密度

    Figure  3.   Distribution density of InDel on genome

    图  4   KEGG富集分析

    Figure  4.   KEGG enrichment scatter plot

    图  5   100尾卵形鲳鲹亲缘关系热图 (a)、主成分分析 (b)、耐低氧性状全基因组关联分析的QQ图 (c) 和曼哈顿图 (d)

    Figure  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

    表  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   变异基因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
  • [1]

    PAERL H W, OTTEN T G. Harmful cyanobacterial blooms: causes, consequences, and controls[J]. Microb Ecol, 2013, 65(4): 995-1010. doi: 10.1007/s00248-012-0159-y

    [2]

    YIN F, GONG H, KE Q, et al. Stress, antioxidant defence and mucosal immune responses of the large yellow croaker Pseudosciaena crocea challenged with Cryptocaryon irritans[J]. Fish Shellfish Immunol, 2015, 47(1): 344-351. doi: 10.1016/j.fsi.2015.09.013

    [3]

    MOREIRA M, SCHRAMA D, SOARES F, et al. Physiological responses of reared sea bream (Sparus aurata Linnaeus, 1758) to an Amyloodinium ocellatum outbreak[J]. J Fish Dis, 2017, 40(11): 1545-1560. doi: 10.1111/jfd.12623

    [4]

    ONUKWUFOR J O, WOOD C M. The osmorespiratory compromise in rainbow trout (Oncorhynchus mykiss): the effects of fish size, hypoxia, temperature and strenuous exercise on gill diffusive water fluxes and sodium net loss rates[J]. Comp Biochem Physiol A, 2018, 219: 10-8.

    [5]

    OBIRIKORANG K A, ACHEAMPONG J N, DUODU C P, et al. Growth, metabolism and respiration in Nile tilapia (Oreochromis niloticus) exposed to chronic or periodic hypoxia[J]. Comp Biochem Physiol A, 2020, 248: 110768. doi: 10.1016/j.cbpa.2020.110768

    [6]

    POULSEN S B, JENSEN L F, NIELSEN K S, et al. Behaviour of rainbow trout Oncorhynchus mykiss presented with a choice of normoxia and stepwise progressive hypoxia[J]. J Fish Biol, 2011, 79(4): 969-979. doi: 10.1111/j.1095-8649.2011.03069.x

    [7]

    WOOD A T, CLARK T D, ELLIOTT N G, et al. The effects of constant and cyclical hypoxia on the survival, growth and metabolic physiology of incubating Atlantic salmon (Salmo salar)[J]. Aquaculture, 2020, 527: 735449. doi: 10.1016/j.aquaculture.2020.735449

    [8]

    SOLLID J, de ANGELIS P, GUNDERSEN K, et al. Hypoxia induces adaptive and reversible gross morphological changes in crucian carp gills[J]. J Exp Biol, 2003, 206(20): 3667-3673. doi: 10.1242/jeb.00594

    [9]

    MITROVIC D, DYMOWSKA A, NILSSON G E, et al. Physiological consequences of gill remodeling in goldfish (Carassius auratus) during exposure to long-term hypoxia[J]. Am J Physiol Regul Integr Comp Physiol, 2009, 297(1): R224-234. doi: 10.1152/ajpregu.00189.2009

    [10]

    WU C B, LIU Z Y, LI F G, et al. Gill remodeling in response to hypoxia and temperature occurs in the hypoxia sensitive blunt snout bream (Megalobrama amblycephala)[J]. Aquaculture, 2017, 479: 479-486. doi: 10.1016/j.aquaculture.2017.06.020

    [11]

    MUNEER P M A, SIVANANDAN R, GOPALAKRISHNAN A, et al. Development and characterization of RAPD and microsatellite markers for genetic variation analysis in the critically endangered yellow catfish Horabagrus nigricollaris (Teleostei: Horabagridae)[J]. Biochem Genet, 2011, 49(1/2): 83-95.

    [12]

    WATANABE T, YOSHIDA M, NAKAJIMA M, et al. Linkage mapping of AFLP and microsatellite DNA markers with the body color-and sex-determining loci in the guppy (Poecilia reticulata)[J]. Zool Sci, 2005, 22(8): 883-889. doi: 10.2108/zsj.22.883

    [13]

    POOMPUANG S, NA-NAKORN U. A preliminary genetic map of walking catfish (Clarias macrocephalus)[J]. Aquaculture, 2004, 232(1/2/3/4): 195-203.

    [14]

    LEI D J, ZHAO G, XIE P, et al. Analysis of genetic diversity of Leuciscus leuciscus baicalensis using novel microsatellite markers with cross-species transferability[J]. Genet Mol Res, 2017, 16(2): 16029376.

    [15]

    TSAI H Y, HAMILTON A, TINCH A E, et al. Genome wide association and genomic prediction for growth traits in juvenile farmed Atlantic salmon using a high density SNP array[J]. BMC Genom, 2015, 16: 969. doi: 10.1186/s12864-015-2117-9

    [16]

    TSAI H Y, HAMILTON A, TINCH A E, et al. Genomic prediction of host resistance to sea lice in farmed Atlantic salmon populations[J]. Genet Sel Evol, 2016, 48(1): 47. doi: 10.1186/s12711-016-0226-9

    [17]

    JIN Y, ZHOU T, GENG X, et al. A genome-wide association study of heat stress-associated SNPs in catfish[J]. Anim Genet, 2017, 48(2): 233-236. doi: 10.1111/age.12482

    [18]

    CHEN F, LAI F L, LUO M J, et al. The genome-wide landscape of small insertion and deletion mutations in Monopterus albus[J]. J Genet Genom, 2019, 46(2): 75-86. doi: 10.1016/j.jgg.2019.02.002

    [19]

    VASEMAGI A, GROSS R, PALM D, et al. Discovery and application of insertion-deletion (INDEL) polymorphisms for QTL mapping of early life-history traits in Atlantic salmon[J]. BMC Genom, 2010, 11: 156. doi: 10.1186/1471-2164-11-156

    [20]

    BRITTEN R J, ROWEN L, WILLIAMS J, et al. Majority of divergence between closely related DNA samples is due to InDels[J]. Proc Nat Acad Sci USA, 2003, 100(8): 4661-4665. doi: 10.1073/pnas.0330964100

    [21]

    ROCKAH-SHMUEL L, TOTH-PETROCZY A, SELA A, et al. Correlated occurrence and bypass of frame-shifting insertion-deletions (InDels) to give functional proteins[J]. PLOS Genet, 2013, 9(10): e1003882. doi: 10.1371/journal.pgen.1003882

    [22] 罗林. 奶牛乳脂性状候选基因SNP及InDel位点筛选及其与产奶性状关联分析[D]. 大庆: 黑龙江八一农垦大学, 2021: 4-8.
    [23] 岳晓鹏. 基于甘蓝型油菜基因组重测序开发InDel标记[D]. 武汉: 华中农业大学, 2014: 6.
    [24] 陈静, 何吉祥, 樊佳佳, 等. 草鱼MyoD基因SNP和InDel标记的筛选及其与生长性状的关联分析[J]. 江苏农业学报, 2018, 34(3): 612-616. doi: 10.3969/j.issn.1000-4440.2018.03.019
    [25]

    CHEN S, ZHOU Y, CHEN Y, et al. Fastp: an ultra-fast all-in-one FASTQ preprocessor[J]. Bioinformatics, 2018, 34(17): 884-890. doi: 10.1093/bioinformatics/bty560

    [26]

    MCKENNA A, HANNA M, BANKS E, et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data[J]. Genome Res, 2010, 20(9): 1297-1303. doi: 10.1101/gr.107524.110

    [27]

    CINGOLANI P, PATEL V M, COON M, et al. Using Drosophila melanogaster as a model for genotoxic chemical mutational studies with a new program, SnpSift[J]. Front Genet, 2012, 3: 35.

    [28]

    WANG K, LI M Y, HAKONARSON H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data[J]. Nucleic Acids Res, 2010, 38(16): e164. doi: 10.1093/nar/gkq603

    [29]

    PURCELL S, NEALE B, TODD-BROWN K, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses[J]. Am J Human Genet, 2007, 81(3): 559-575. doi: 10.1086/519795

    [30]

    YANG J, LEE S H, GODDARD M E, et al. GCTA: a tool for genome-wide complex trait analysis[J]. Am J Human Genet, 2011, 88(1): 76-82. doi: 10.1016/j.ajhg.2010.11.011

    [31]

    ZHOU X, STEPHENS M. Genome-wide efficient mixed-model analysis for association studies[J]. Nature Genet, 2012, 44(7): 821-824. doi: 10.1038/ng.2310

    [32]

    ZHONG X, WANG X, ZHOU T, et al. Genome-wide association study reveals multiple novel QTL associated with low oxygen tolerance in hybrid catfish[J]. Mar Biotechnol (NY), 2017, 19(4): 379-390. doi: 10.1007/s10126-017-9757-5

    [33]

    WANG X, LIU S, JIANG C, et al. Multiple across-strain and within-strain QTLs suggest highly complex genetic architecture for hypoxia tolerance in channel catfish[J]. Mol Genet Genom, 2017, 292(1): 63-76. doi: 10.1007/s00438-016-1256-2

    [34] 杨洁, 赫佳, 王丹碧, 等. InDel标记的研究和应用进展[J]. 生物多样性, 2016, 24(2): 237-243. doi: 10.17520/biods.2015205
    [35]

    SAN L Z, LIU B S, LIU B, et al. Genome-wide association study reveals multiple novel SNPs and putative candidate genes associated with low oxygen tolerance in golden pompano Trachinotus ovatus (Linnaeus 1758)[J]. Aquaculture, 2021, 544: 737098. doi: 10.1016/j.aquaculture.2021.737098

    [36] 王慧芳, 周光现, 孙永峰, 等. 基于全基因组重测序技术的狮头鹅InDel标记分析[J]. 畜牧兽医学报, 2021, 52(3): 662-675. doi: 10.11843/j.issn.0366-6964.2021.03.010
    [37] 吴迷, 汪念, 沈超, 等. 基于重测序的陆地棉InDel标记开发与评价[J]. 作物学报, 2019, 45(2): 196-203.
    [38] 常玉晓. 水稻T-DNA插入突变体侧翼序列的分离和水稻DNA复制蛋白RPA1a及RPA2-3基因的功能研究[D]. 武汉: 华中农业大学, 2010: 16-20.
    [39] 黄芳. 斑马鱼DNA聚合酶Delta四亚基和相关因子的制备鉴定及赤点石斑鱼IκBα基因的克隆与功能分析[D]. 镇江: 江苏大学, 2016: 3-8.
    [40]

    ZEBDA N, DUBROVSKYI O, BIRUKOV K G. Focal adhesion kinase regulation of mechanotransduction and its impact on endothelial cell functions[J]. Microvasc Res, 2012, 83(1): 71-81. doi: 10.1016/j.mvr.2011.06.007

    [41] 黄智康, 江世贵, 周发林, 等. 基于InDel标记的斑节对虾早期性别鉴定方法的建立[J]. 南方水产科学, 2020, 16(3): 113-118. doi: 10.12131/20190222
    [42] 牛姣姣. 芝麻枯萎病抗性关联分析[D]. 南京: 南京农业大学, 2014: 24-43.
    [43] 薛蕾. 山羊InDels筛选鉴定及其与生产性状关联分析[D]. 重庆: 西南大学, 2018: 30-33.
    [44]

    YANG Z L, ZOU L Q, SUN T T, et al. Genome-wide association study using whole-genome sequencing identifies a genomic region on chromosome 6 associated with comb traits in Nandan-Yao chicken[J]. Front Genet, 2021, 12: 682501. doi: 10.3389/fgene.2021.682501

    [45]

    ELLIS J M, WONG W, WOLFGANG M J. Acyl coenzyme a thioesterase 7 regulates neuronal fatty acid metabolism to prevent neurotoxicity[J]. Mol Cell Biol, 2013, 33(9): 1869-1882. doi: 10.1128/MCB.01548-12

    [46]

    MAEKAWA R, MUTO H, HATAYAMA M, et al. Dysregulation of erythropoiesis and altered erythroblastic NMDA receptor-mediated calcium influx in Lrfn2-deficient mice[J]. PLOS ONE, 2021, 16(1): e0245624. doi: 10.1371/journal.pone.0245624

    [47]

    FORSTER J R, LOCHNIT G, STOHR H. Proteomic analysis of the membrane palmitoylated protein-4 (MPP4)-associated protein complex in the retina[J]. Exp Eye Res, 2009, 88(1): 39-46. doi: 10.1016/j.exer.2008.09.016

    [48]

    RUIZ-HERNANDEZ A, ROMERO-NAVA R, HUANG F Y, et al. Altered function and expression of the orphan GPR135 at the cardiovascular level in diabetic Wistar rats[J]. J Recept Sig Transd, 2018, 38(5/6): 484-491.

    [49]

    GROS-LOUIS F, KRIZ J, KABASHI E, et al. Als2 mRNA splicing variants detected in KO mice rescue severe motor dysfunction phenotype in Als2 knock-down zebrafish[J]. Human Mol Genet, 2008, 17(17): 2691-2702. doi: 10.1093/hmg/ddn171

    [50]

    HUANG L J, SZYMANSKA K, JENSEN V L, et al. TMEM237 is mutated in individuals with a joubert syndrome related disorder and expands the role of the TMEM family at the ciliary transition zone[J]. Am J Human Genet, 2011, 89(6): 713-730. doi: 10.1016/j.ajhg.2011.11.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
图(5)  /  表(2)
计量
  • 文章访问数:  758
  • HTML全文浏览量:  201
  • PDF下载量:  46
  • 被引次数: 11
出版历程
  • 收稿日期:  2021-11-23
  • 修回日期:  2022-02-05
  • 录用日期:  2022-02-17
  • 网络出版日期:  2022-03-02
  • 刊出日期:  2022-10-04

目录

/

返回文章
返回