Transcriptomic WGCNA analyses reveal endoplasmic reticulum response of Patinopecten yessoensis under acute heat stress

Abstract

Patinopecten yessoensis is one of the most important aquaculture shellfish worldwide. As a cold-water species, the increasing seawater temperatures due to global warming have been a great challenge for P. yessoensis. Understanding the response mechanism of P. yessoensis to high temperatures is essential to the genetic breeding of heat-tolerant varieties. In this study, histological, cytological and transcriptomic analyses were performed on the gill tissues of P. yessoensis under different acute heat stresses (temperature: 22 °C, 25 °C, 28 °C; time: 6 h, 36 h). A total of 10,601 differentially expressed genes (DEGs) were detected, among which, 141 DEGs were shared by all heat-stressed groups. Weighted gene co-expression network analysis (WGCNA) identified three modules (M1, M2 and M3) that were significantly correlated with heat stress. Functional enrichment analysis revealed that heat stress probably caused serious endoplasmic reticulum (ER) stress in gill cells, with many associated GO terms and KEGG pathways markedly enriched. Notably, the protein processing in ER pathway probably played a key role in the response to heat stress; this pathway was upregulated and was the most enriched in all heat-stressed groups, as well as in the M1 and M2 modules. DEGs from this pathway participated in different processes of the ER stress response, alleviating ER stress by upregulating the protein correct folding, ER-associated degradation (ERAD) and unfolded protein response (UPR) processes. Excessive ER stress induced apoptosis via the UPR, and excessive apoptosis resulted in tissue injury as observed by histological and TUNEL assays. Further, numerous Hsps were upregulated in all heat-stressed groups, and some members were identified as hub genes, including BiP, NEF, Hsp70 and Hsp90. These genes may perform crucial functions during heat stress as molecular chaperones. This study elucidates the key molecular components of the ER stress response to heat stress in scallops. The findings provide new insights into the adaptation mechanisms of molluscs to environmental stress.