Molecular Mechanism of Pacific Oyster Response to Ocean Acidification ——Based on Transcriptome Analysis

doi:10.11924/j.issn.1000-6850.casb2022-0785

Abstract

Abstract:

The aim of this study was to explore the molecular mechanism of the response of Pacific oyster to ocean acidification through transcriptome analysis. In this study, the RNA-seq data (including pH 7.8, pH 7.4, pH 7.0 and pH 6.6 groups) of Pacific oyster under ocean acidification stress in NCBI were used to analyze the differential genes, GO and KEGG by Fastp, Hisat2, Samtools and R language analysis. The results showed that there were 61, 93 and 943 significantly differentially expressed genes in pH 7.4 vs pH 7.8, pH 7.0 vs pH 7.8 and pH 6.6 vs pH 7.8 groups. GO analysis of the significantly differentially expressed genes in pH 6.6 vs pH 7.8 group showed that these genes were significantly enriched in metabolism, genetic information processing, human disease and other molecular functions. KEGG enrichment analysis showed that these genes were significantly enriched in glucose metabolism, lipid metabolism and immune function. Expression analysis showed that with the decrease of pH, the amount of gene expression related to energy metabolism and immunity decreased significantly. For example, the expression of glycogen phosphorylase in saccharometabolism decreased from 41.5 to 0. The expression of arginine kinase 1 in amino acid metabolism decreased from 1438 to 0, and the expression of arginine kinase 2 decreased from 27 to 3. The expression of fatty acid synthetase in fat metabolism decreased from 28 to 2.5. The significant inhibition of these genes may lead to the disturbance of energy metabolism in oyster body. In addition, the expression of immune-related genes, such as CD151, decreased from 45 almost to 0, and was extremely significantly inhibited, which may alter osmotic pressure in oysters, lead to immune system disorders and reduce their ability to withstand changes in the external environment. In conclusion, ocean acidification is likely to lead to the disorder of metabolism and immune system of oyster body, and reduce the ability of Pacific oyster to resist environmental changes.

Key words: Pacific oyster, ocean acidification, transcriptome, differentially expressed genes, energy metabolism

LU Ziya, CHEN Xiaolin, HUANG Simin, DUAN Qianqian, PENG Boyan, XIAO Zijian, XIANG Zaiying, GUO Xiaomei, LIU Yaqi, LIN Siqing, TAN Karsoon, ZHANG Hongkuan, ZHENG Huaiping. Molecular Mechanism of Pacific Oyster Response to Ocean Acidification ——Based on Transcriptome Analysis[J]. Chinese Agricultural Science Bulletin, 2023, 39(26): 137-146.

Figures/Tables 7

References 27

[1]	张美昭. 海洋渔业产业发展现状与前景研究[M]. 广东: 广东经济出版社, 2018:154.
[2]	邹琰, 刘广斌, 张天文, 等. 我国太平洋牡蛎主要养殖模式[J]. 科学养鱼, 2021(11):63-64.
[3]	农业农村部渔业渔政管理局. 2021中国渔业统计年鉴[M]. 北京: 中国农业出版社, 2021.
[4]	李娜, 孙敏, 王春华, 等. 水产品保鲜贮藏期间品质评价方式的研究进展[J]. 食品安全质量检测学报, 2022, 13(4):1099-1105.
[5]	王婷, 郑佳慧, 胡梦红, 等. 海洋酸化对贝类的生理生态学影响研究进展[J]. 海洋科学, 2022, 46(1):192-202.
[6]	CALDEIRA K, WICKETT M E. Anthropogenic carbon and ocean pH[J]. Nature, 2003, 425(6956):365. doi: 10.1038/425365a
[7]	孟紫强. 生态毒理学[M]. 北京: 中国环境出版社, 2019:245.
[8]	章恒, 赵晟, 朱爱意. 厚壳贻贝的钙化和呼吸对海洋酸化的响应[J]. 安徽农业科学, 2015, 43(24):104-105.
[9]	莫知. 海洋酸化之牡蛎杀手[J]. 海洋世界, 2010(9):37.
[10]	COOLEY S R, DONEY S C. Anticipating ocean acidification's economic consequences for commercial fisheries[J]. Environmental research letters, 2009, 4(2):2-4.
[11]	PARKER L M, ROSS P M, O'CONNOR W A, et al. Predicting the response of molluscs to the impact of ocean acidification[J]. Biology, 2013, 2(2):676-679.
[12]	FEDERICA S, ASSAF M, PUTNAM H M, et al. Genetic and physiological traits conferring tolerance to ocean acidification in mesophotic corals[J]. Global change biology, 2021, 27(20):12-13.
[13]	KURIHARA H. Effects of CO₂-driven ocean acidification on the early developmental stages of invertebrates[J]. Marine ecology progress series, 2008, 373:276-282.
[14]	王晓芹. 海洋酸化胁迫对紫贻贝和长牡蛎生理活动的影响及其缓解途径[D]. 上海: 上海海洋大学, 2018:24-25.
[15]	程青丽, 刘文颖, 赵晓涵, 等. 水产品中过敏原精氨酸激酶研究进展[J]. 食品工业, 2021, 42(5):403-406.
[16]	王婷, 郑佳慧, 胡梦红, 等. 海洋酸化对贝类的生理生态学影响研究进展[J]. 海洋科学, 2022, 46(1):194.
[17]	TU Y J, SU Y J, LI G H, et al. Expression of lipid metabolism-associated genes in male and female white feather chicken[D]. Japan poultry science association, 2016.
[18]	王爱民. 饲料脂肪水平对吉富罗非鱼生长及脂肪代谢调节的研究[D]. 南京: 南京农业大学, 2011:19-20.
[19]	向枭, 曾本和, 王睿, 等. 胆汁酸对齐口裂腹鱼幼鱼脂肪沉积、脂肪代谢酶活性及相关基因表达的影响[J]. 水产学报, 2022, 46(6):1045-1052.
[20]	王晓芹, RASTRICK SPS, 吴亚林, 等. 海水酸化胁迫对紫贻贝能量分配的影响[J]. 渔业科学进展, 2019, 40(3):21-30.
[21]	刘世良, 麦康森. 贝类免疫系统和机理的研究进展[J]. 海洋学报(中文版), 2003(2):95-105.
[22]	YOO S H, LEE K, CHAE J Y, et al. CD151 expression can predict cancer progression in clear cell renal cell carcinoma[J]. Histopathology, 2011, 58(2):191-197. doi: 10.1111/his.2011.58.issue-2 URL
[23]	WANG H X, LI Q L, CHANDAN S, et al. Tetraspanin protein contributions to cancer[J]. Biochemical society transactions, 2011, 39(2):547-552 doi: 10.1042/BST0390547 URL
[24]	WANG Q, CAO R W, NING X X, et al. Effects of ocean acidification on immune responses of the Pacific oyster Crassostrea gigas[J]. Fish and shellfish immunology, 2016, 49:24-33. doi: 10.1016/j.fsi.2015.12.025 URL
[25]	魏磊. 长牡蛎(Crassostrea gigas)对海水酸化生理响应的组学研究[D]. 烟台: 中国科学院烟台海岸带研究所, 2015:63-69.
[26]	LIU S X, WEI X, CHENG G, et al. Ocean acidification weakens the immune response of blood clam through hampering the NF-kappa β and toll-like receptor pathways[J]. Fish and shellfish immunology, 2016, 54:322-327. doi: 10.1016/j.fsi.2016.04.030 URL
[27]	尹双, 代燕辉, 徐立娜, 等. 海洋酸化条件下重金属和金属纳米颗粒的环境行为及其对海洋生物的影响[J]. 中国科学:化学, 2018, 48(3):256-265.

samples	reads	overall alignment rate/%	Unique mapping/%	Multiple mapping/%
pHT_6.6_16A	32643672	78.91	74.41	4.50
pHT_6.6_16B	37993685	77.51	72.99	4.53
pHT_6.6_16C	28818302	77.30	72.76	4.53
pHT_6.6_16D	34403075	77.82	73.35	4.47
pHT_6.6_16E	28612020	76.54	71.96	4.58
pHT_7.0_11A	33906611	78.23	73.44	4.79
pHT_7.0_11B	44061225	76.43	71.74	4.70
pHT_7.0_11C	35516169	76.32	71.62	4.70
pHT_7.0_11D	25684693	78.75	74.44	4.31
pHT_7.0_11E	25482477	79.04	74.60	4.44
pHT_7.4_13A	45920961	79.36	74.31	5.05
pHT_7.4_13B	31837082	77.55	72.32	5.23
pHT_7.4_13C	33694916	77.99	73.14	4.85
pHT_7.4_13D	19610155	78.70	74.47	4.23
pHT_7.4_13E	34750182	78.55	73.98	4.57
pHT_7.8_3A	30350518	79.08	72.96	6.12
pHT_7.8_3B	36727612	76.26	71.31	4.94
pHT_7.8_3C	39098640	77.07	72.39	4.68
pHT_7.8_3D	29062777	78.62	74.06	4.56
pHT_7.8_3E	38749898	77.31	72.53	4.78

samples	reads	overall alignment rate/%	Unique mapping/%	Multiple mapping/%
pHT_6.6_16A	32643672	78.91	74.41	4.50
pHT_6.6_16B	37993685	77.51	72.99	4.53
pHT_6.6_16C	28818302	77.30	72.76	4.53
pHT_6.6_16D	34403075	77.82	73.35	4.47
pHT_6.6_16E	28612020	76.54	71.96	4.58
pHT_7.0_11A	33906611	78.23	73.44	4.79
pHT_7.0_11B	44061225	76.43	71.74	4.70
pHT_7.0_11C	35516169	76.32	71.62	4.70
pHT_7.0_11D	25684693	78.75	74.44	4.31
pHT_7.0_11E	25482477	79.04	74.60	4.44
pHT_7.4_13A	45920961	79.36	74.31	5.05
pHT_7.4_13B	31837082	77.55	72.32	5.23
pHT_7.4_13C	33694916	77.99	73.14	4.85
pHT_7.4_13D	19610155	78.70	74.47	4.23
pHT_7.4_13E	34750182	78.55	73.98	4.57
pHT_7.8_3A	30350518	79.08	72.96	6.12
pHT_7.8_3B	36727612	76.26	71.31	4.94
pHT_7.8_3C	39098640	77.07	72.39	4.68
pHT_7.8_3D	29062777	78.62	74.06	4.56
pHT_7.8_3E	38749898	77.31	72.53	4.78