Gills full-length transcriptomic analysis of osmoregulatory adaptive responses to salinity stress in Coilia nasus

Comprehensive analysis of hub mRNA, lncRNA and miRNA, and associated ceRNA networks implicated in cobia (Rachycentron canadum) scales under hypoosmotic adaption

Salinity plays a vital role in fish aquaculture, profoundly influencing the growth and development of fish. Scales, as the protective outer layer of fish, function as a critical defense against external factors. In this study, we employed transcriptome sequencing to analyze the ceRNA expression profiles to reveal the effect of salinity acclimation on transcriptional expression changes in the scales of cobia (Rachycentron canadum). The results revealed that after being exposed to a salinity level of 15 ‰ for just one day (1D), a total of 407 mRNAs/genes were significantly regulated; 66 miRNAs were respectively significantly regulated; and 109 target genes of the differentially expressed miRNAs were significantly regulated; a total of 185 differently expressed lncRNAs and 292 differently expressed target genes (DetGenes) of differently expressed lncRNAs were also identified. After 7 days (7D), a total of 2195 mRNAs/genes were found to be significantly regulated and 82 miRNAs were significantly regulated; among the target genes of the differentially expressed miRNAs, 245 were regulated. Moreover, 438 differently expressed lncRNAs and 681 DetGenes of these lncRNAs were identified. Subsequent analysis through GO, KEGG pathway, in 1D vs. CG (control group), DeGenes, which first respond to changes in salinity, are mainly involved in negative regulation of macrophage differentiation, negative regulation of granulocyte differentiation and negative regulation of phagocytosis, and are mainly related to biological processes related to the immune function of fish. After a 7-day process, DeGenes were enriched in the collagen fibril organization, regulation of nodal signaling pathway and cell recognition biology processes. These biological processes are not only related to the immune function of fish, but more importantly, to the physiological structure of fish. By analyzing the co down-regulated miRNAs of 1D vs. CG, as well as 7D vs. CG, the functions of these miRNAs are mainly related to bone differentiation and development. In addition，ceRNA network uncovered that the effect of salinity is temporal. The first competing lncRNAs mainly regulated genes related to physiological processes and biological development, while target genes related to immunity and body defense were less competitive. On the contrary, after a period of salinity treatment, the types of competing lncRNAs involved changed.

Identification of Genes and Long Non-Coding RNAs Putatively Related to Portunus trituberculatus Sex Determination and Differentiation Using Oxford Nanopore Technology Full-Length Transcriptome Sequencing

The swimming crab (Portunus trituberculatus) is an economically important species in China, and its growth traits show obvious sexual dimorphism. Thus, it is important to study the mechanism of sex determination and differentiation in this species. Herein, we identified 2138 differentially expressed genes and 132 differentially expressed long non-coding RNAs (lncRNAs) using Oxford Nanopore Technology full-length transcriptome sequencing. We predicted 561 target genes of the differentially expressed lncRNAs according to their location and base pair complimentary principles. Furthermore, pathways related to sex determination, differentiation, and reproduction were enriched for lncRNAs according to gene ontology and Kyoto Encyclopedia of Genes and Genomes analyses. This indicated that lncRNAs might play regulatory roles in these pathways. Our results could form the basis for future studies of sex determination and differentiation in P. trituberculatus.

Alternative splicing patterns of hnrnp genes in gill tissues of rainbow trout (Oncorhynchus mykiss) during salinity changes

Alternative splicing (AS) plays an important role in various physiological processes in eukaryotes, such as the stress response. However, patterns of AS events remain largely unexplored during salinity acclimation in fishes. In this study, we conducted AS analysis using RNA-seq datasets to explore splicing patterns in the gill tissues of rainbow trout exposed to altered salinity environments, ranging from 0 ‰ (T0) to 30 ‰ (T30). The results revealed 1441, 351, 483, 1051 and 1049 differentially alternatively spliced (DAS) events in 5 pairwise comparisons, including T6 vs. T0, T12 vs. T0, T18 vs. T0, T24 vs. T0, and T30 vs. T0, respectively. These DAS events were derived from 1290, 328, 444, 963 and 948 genes. Enrichment analysis indicated that these DAS genes were related to RNA splicing and processing. Among these, 14 DAS genes were identified as members of the large heterogeneous nuclear RNP (hnRNP) gene family. Alternative 3' splice site (A3SS), exon skipping (SE) and intron retention (RI) events resulted in the fragmentation or even loss of the functional RNA recognition motif (RRM) domains in hnrnpa0, hnrnp1a, hnrnp1b and hnrnpc genes. The incomplete RRM domains would hinder the interactions between hnRNP genes and pre-mRNAs. It would in turn influence the splicing patterns and mRNA stability of downstream target genes in response to salinity changes. The study provides insights into salinity acclimation in gill tissues of rainbow trout and serves as a significant reference on the osmoregulation mechanisms at post-transcription regulation levels in fish.

Transcriptome Responses to Different Salinity Conditions in Litoditis marina, Revealed by Long-Read Sequencing

The marine nematode Litoditis marina is widely distributed in intertidal zones around the globe, yet the mechanisms underlying its broad adaptation to salinity remain elusive. In this study, we applied ONT long-read sequencing technology to unravel the transcriptome responses to different salinity conditions in L. marina. Through ONT sequencing under 3‱, 30‱ and 60‱ salinity environments, we obtained 131.78 G clean data and 26,647 non-redundant long-read transcripts, including 6464 novel transcripts. The DEGs obtained from the current ONT lrRNA-seq were highly correlated with those identified in our previously reported Illumina short-read RNA sequencing data. When we compared the 30‱ to the 3‱ salinity condition, we found that GO terms such as oxidoreductase activity, cation transmembrane transport and ion transmembrane transport were shared between the ONT lrRNA-seq and Illumina data. Similarly, GO terms including extracellular space, structural constituents of cuticle, substrate-specific channel activity, ion transport and substrate-specific transmembrane transporter activity were shared between the ONT and Illumina data under 60‱ compared to 30‱ salinity. In addition, we found that 79 genes significantly increased, while 119 genes significantly decreased, as the salinity increased. Furthermore, through the GO enrichment analysis of 214 genes containing DAS, in 30‱ compared to 3‱ salinity, we found that GO terms such as cellular component assembly and coenzyme biosynthetic process were enriched. Additionally, we observed that GO terms such as cellular component assembly and coenzyme biosynthetic process were also enriched in 60‱ compared to 30‱ salinity. Moreover, we found that 86, 125, and 81 genes that contained DAS were also DEGs, in comparisons between 30‱ and 3‱, 60‱ and 30‱, and 60‱ and 3‱ salinity, respectively. In addition, we demonstrated the landscape of alternative polyadenylation in marine nematode under different salinity conditions This report provides several novel insights for the further study of the mechanisms by which euryhalinity formed and evolved, and it might also contribute to the investigation of salinity dynamics induced by global climate change.

Investigating the Impact of Disrupting the Glutamine Metabolism Pathway on Ammonia Excretion in Crucian Carp (Carassius auratus) under Carbonate Alkaline Stress Using Metabolomics Techniques

With the gradual decline in freshwater resources, the space available for freshwater aquaculture is diminishing and the need to maximize saline water for aquaculture is increasing. This study aimed to elucidate the impact mechanisms of the disruption of the glutamate pathway on serum metabolism and ammonia excretion in crucian carp (Carassius auratus) under carbonate alkaline stress. A freshwater control group (C group), a 20 mmol/L NaHCO3 stress group (L group), and a 40 mmol/L NaHCO3 stress group (H group) were established. After 30 days of exposure, methionine sulfoximine (MSO) was injected to block the glutamate pathway metabolism, and the groups post-blocking were labeled as MC, ML, and MH. Ultra-high-performance liquid chromatography coupled with the quadrupole time-of-flight mass spectrometry (UPLC-Q-TOF/MS) metabolomics technique was employed to detect changes in the composition and content of crucian carp serum metabolites. Significant differential metabolites were identified, and related metabolic pathways were analyzed. The results revealed that, following the glutamate pathway blockade, a total of 228 differential metabolites (DMs) were identified in the three treatment groups. An enrichment analysis indicated significant involvement in glycerophospholipid metabolism, arachidonic acid metabolism, sphingolipid metabolism, purine metabolism, arginine and proline biosynthesis, pantothenate and CoA biosynthesis, glutathione metabolism, and fatty acid degradation, among other metabolic pathways. The results showed that ROS imbalances and L-arginine accumulation in crucian carp after the glutamate pathway blockade led to an increase in oxidative stress and inflammatory responses in vivo, which may cause damage to the structure and function of cell membranes. Crucian carp improves the body's antioxidant capacity and regulates cellular homeostasis by activating glutathione metabolism and increasing the concentration of phosphatidylcholine (PC) analogs. Additionally, challenges such as aggravated ammonia excretion obstruction and disrupted energy metabolism were observed in crucian carp, with the upregulation of purine metabolism alleviating ammonia toxicity and maintaining energy homeostasis through pantothenate and CoA biosynthesis as well as fatty acid degradation. This study elucidated the metabolic changes in crucian carp under carbonate alkaline stress after a glutamate pathway blockade at the cellular metabolism level and screened out the key metabolic pathways, which provide a scientific basis for further in-depth studies on the ammonia excretion of freshwater scleractinian fishes under saline and alkaline habitats at a later stage.

Effects of low salinity stress on osmoregulation and gill transcriptome in different populations of mud crab Scylla paramamosain

Animals living in estuaries suffer from rapid and continuous salinity fluctuations, while the global warming and extreme precipitation aggravate this situation. Osmoregulation is important for estuarine animals adapt to salinity fluctuations. The present study investigated the effects of low salinity stress on osmoregulation and gill transcriptome in two populations of mud crab from Hangzhou Bay and Zhangzhou Bay of China, respectively. Crabs were transferred from salinity 25 ppt to 5 ppt for 96 h. Edematous swelling in gill filaments was caused by low salinity stress and was more serious in Zhangzhou Bay population. Gill Na⁺/K⁺-ATPase activity increased (p < 0.01) in both populations under the low salinity stress and was significantly higher (p < 0.01) in Hangzhou Bay population than in Zhangzhou Bay population. According to transcriptome analysis, there were 191 genes differentially expressed under the low salinity stress in gill tissue of both populations. Several ion transport and energy metabolism related pathways, as well as the arginine and proline metabolism pathway, were enriched by these genes. On the other hand, 272 genes were identified to differentially express between two populations under the low salinity stress, but not under the control salinity. The enrichment analysis showed that these genes were mainly related to ion transport, energy metabolism, osmolytes metabolism and methyltransferase activity. In conclusion, the present study suggested that mud crab exploited a combination of extracellular anisosmotic regulation and intracellular isosmotic regulation for osmoregulation under the low salinity stress. Hangzhou Bay population showed a greater osmoregulatory capacity, which is probably due to the enhanced ion transport, energy supply, and osmolytes regulation. Meanwhile, epigenetic modification might also contribute to an inherent osmoregulation ability for Hangzhou Bay population to response to salinity fluctuation rapidly.

Transcriptome and molecular regulatory mechanisms analysis of gills in the black tiger shrimp Penaeus monodon under chronic low-salinity stress

Background: Salinity is one of the main influencing factors in the culture environment and is extremely important for the survival, growth, development and reproduction of aquatic animals. Methods: In this study, a comparative transcriptome analysis (maintained for 45 days in three different salinities, 30 psu (HC group), 18 psu (MC group) and 3 psu (LC group)) was performed by high-throughput sequencing of economically cultured Penaeus monodon. P. monodon gill tissues from each treatment were collected for RNA-seq analysis to identify potential genes and pathways in response to low salinity stress. Results: A total of 64,475 unigenes were annotated in this study. There were 1,140 upregulated genes and 1,531 downregulated genes observed in the LC vs. HC group and 1,000 upregulated genes and 1,062 downregulated genes observed in the MC vs. HC group. In the LC vs. HC group, 583 DEGs significantly mapped to 37 signaling pathways, such as the NOD-like receptor signaling pathway, Toll-like receptor signaling pathway, and PI3K-Akt signaling pathway; in the MC vs. HC group, 444 DEGs significantly mapped to 28 signaling pathways, such as the MAPK signaling pathway, Hippo signaling pathway and calcium signaling pathway. These pathways were significantly associated mainly with signal transduction, immunity and metabolism. Conclusions: These results suggest that low salinity stress may affect regulatory mechanisms such as metabolism, immunity, and signal transduction in addition to osmolarity in P. monodon. The greater the difference in salinity, the more significant the difference in genes. This study provides some guidance for understanding the low-salt domestication culture of P. monodon.

Transcriptomic Analysis in Marine Medaka Gill Reveals That the Hypo-Osmotic Stress Could Alter the Immune Response via the IL17 Signaling Pathway

Fish gills are the major osmoregulatory tissue that contact the external water environment and have developed an effective osmoregulatory mechanism to maintain cellular function. Marine medaka (Oryzias melastigma) has the ability to live in both seawater and fresh water environments. The present study performed a seawater (SW) to 50% seawater (SFW) transfer, and the gill samples were used for comparative transcriptomic analysis to study the alteration of hypo-osmotic stress on immune responsive genes in this model organism. The result identified 518 differentiated expressed genes (DEGs) after the SW to SFW transfer. Various pathways such as p53 signaling, forkhead box O signaling, and the cell cycle were enriched. Moreover, the immune system was highlighted as one of the top altered biological processes in the enrichment analysis. Various cytokines, chemokines, and inflammatory genes that participate in the IL-17 signaling pathway were suppressed after the SW to SFW transfer. On the other hand, some immunoglobulin-related genes were up-regulated. The results were further validated by real-time qPCR. Taken together, our study provides additional gill transcriptome information in marine medaka; it also supports the notion that osmotic stress could influence the immune responses in fish gills.

PacBio Full-Length and Illumina Transcriptomes of the Gill Reveal the Molecular Response of Corbicula fluminea under Aerial Exposure

Air exposure is a common stress for Corbicula fluminea, an economically important freshwater shellfish consumed in China, during aquaculture and transportation. However, little is known about its molecular responses to air exposure. Therefore, this study used a combination of PacBio full-length and Illumina transcriptomes to investigate its molecular responses to air exposure. A total of 36,772 transcripts were obtained using PacBio sequencing. Structural analysis identified 32,069 coding sequences, 1906 transcription factors, 8873 simple sequence repeats, and 17,815 long non-coding RNAs. Subcellular localization analysis showed that most transcripts were located in the cytoplasm and nucleus. After 96-h of air exposure, 210 differentially expressed genes (DEGs) in the gill were obtained via Illumina sequencing. Among these DEGs, most of the genes related to glycolysis, tricarboxylic acid cycle, lipid metabolism, and amino acid metabolism were upregulated. Additionally, many DEGs associated with immunity, cytoskeleton reorganization, autophagy, and ferroptosis were identified. These findings indicated that metabolic strategy change, immune response, cytoskeleton reconstruction, autophagy, and ferroptosis might be the important mechanisms that C. fluminea use to cope with air exposure. This study will enrich the gene resources of C. fluminea and provide valuable data for studying the molecular mechanisms coping with air exposure in C. fluminea and other freshwater mollusks.