The identification and characteristics of salinity-related microRNAs in gills of Portunus trituberculatus

miRNA-seq provides novel insight into the response to hyper- and hypo- salinity acclimation in Crassostrea hongkongensis

The Hong Kong oyster, Crassostrea hongkongensis, is a significant bivalve species with economic importance. It primarily inhabits the estuarine intertidal zones in southern China, making it susceptible to salinity fluctuations. Consequently, investigating the molecular mechanisms governing salinity regulation in C. hongkongensis is essential. In this study, we conducted miRNA-seq on C. hongkongensis to compare miRNA expression differences under varying salinities (5‰, 25‰, and 35‰). The miRNA sequencing revealed 51 known miRNAs and 95 novel miRNAs across nine small RNA libraries (S5, S25, and S35). Among these miRNAs, we identified 6 down-regulated differentially expressed (DE) miRNAs in response to hypo-salinity stress (5‰), while 1 up-regulated DE miRNA and 5 down-regulated DE miRNAs were associated with hyper-salinity stress (35‰). Additionally, we predicted 931 and 768 potential target genes for hypo- and hyper-salinity stress, respectively. Functional gene annotation indicated that the target genes under hypo-salinity stress were linked to vesicle-mediated transport and metal ion binding. Conversely, those under hyper-salinity stress were primarily involved in signal transduction and metabolic processes. These findings have provided insights into the regulatory role of miRNAs, their potential target genes and associated pathways in oyster hypo- and hyper-salinity stress, which establish a foundation for future studies on the roles of miRNAs in salinity acclimation mechanisms in C. hongkongensis.

Identification and functional characterization of a novel TRPA1 gene from sea cucumber Apostichopus japonicus and interaction with miR-2013 in response to salt stress

Salinity is important abiotic factor influencing sea cucumber aquaculture. This study aimed to identify and functional study of a novel transient receptor potential cation channel subfamily A member 1 (TRPA1) involved in salinity stress through interaction with miR-2013 in the sea cucumber. The full-length cDNA sequence was 1369 bp in length and encoded 138 amino acids. The TRPA1 homolog protein was a hydrophilic protein without a signal peptide and was predicted to a spatial structure of seven helices and eight random coils and two major ANK functional domains. Bioinformatic analysis and luciferase reporter assays confirmed TRPA1 as a target gene of miR-2013. Quantitative PCR revealed that miR-2013 was induced upregulation after salinity stress, while TRPA1 showed upregulated expression with maximum expression at 24 h. The expression of miR-2013 and TRPA1 was negatively regulated. Transfection experiments were conducted to validate the role of miR-2013 and TRPA1 in salinity response. The results showed that miR-2013 was upregulated and TRPA1 was downregulated after transfection with miR-2013 mimics, while miR-2013 was downregulated and TRPA1 was upregulated after transfection with miR-2013 inhibitor. Transfection with si-TRPA1 homolog resulted in upregulation of miR-2013 and downregulation of TRPA1 homolog. These findings suggest that miR-2013 can regulate the expression of TRPA1 under salt stress, and highlight the importance of miR-2013 and TRPA1 in salt stress response. miR-2013 mimics improved the survival rate, while miR-2013 inhibitor and si-TRPA1 reduced it. These findings suggest that miR-2013 and TRPA1 play important roles in sea cucumbers adaptation to salinity changes.

Comprehensive Analysis of Whole-Transcriptome Profiles in Response to Acute Hypersaline Challenge in Chinese Razor Clam Sinonovacula constricta

The Chinese razor clam (Sinonovacula constricta) is an important for Chinese aquaculture marine bivalve that naturally occurs across intertidal and estuarine areas subjected to significant changes in salinity level. However, the information on the molecular mechanisms related to high salinity stress in the species remain limited. In this study, nine gill samples of S. constricta treated with 20, 30, and 40 ppt salinity for 24 h were used for whole-transcriptome RNA sequencing, and a regulatory network of competing endogenous RNAs (ceRNAs) was constructed to better understand the mechanisms responsible for adaptation of the species to high salinity. A total of 83,262 lncRNAs, 52,422 mRNAs, 2890 circRNAs, and 498 miRNAs were identified, and 4175 of them displayed differential expression pattern among the three groups examined. The KEGG analyses of differentially expressed RNAs evidenced that amino acid synthesis and membrane transport were the dominant factors involved in the adaptation of the Chinese razor clam to acute salinity increase, while lipid metabolism and signaling played only a supporting role. In addition, lncRNA/circRNA-miRNA-mRNA regulatory networks (ceRNA network) showed clearly regulatory relationships among different RNAs. Moreover, the expression of four candidate genes, including tyrosine aminotransferase (TAT), hyaluronidase 4 (HYAL4), cysteine sulfinic acid decarboxylase (CSAD), and ∆¹-pyrroline-5-carboxylate synthase (P5CS) at different challenge time were detected by qRT-PCR. The expression trend of TAT and HYAL4 was consistent with that of the ceRNA network, supporting the reliability of established network. The expression of TAT, CSAD, and P5CS were upregulated in response to increased salinity. This might be associated with increased amino acid synthesis rate, which seems to play an essential role in adaptation of the species to high salinity stress. In contrast, the expression level of HYAL4 gene decreased in response to elevated salinity level, which is associated with reduction Hyaluronan hydrolysis to help maintain water in the cell. Our findings provide a very rich reference for understanding the important role of ncRNAs in the salinity adaptation of shellfish. Moreover, the acquired information may be useful for optimization of the artificial breeding of the Chinese razor clam under aquaculture conditions.

Salinity-induced transcriptome profiles in marine and freshwater threespine stickleback after an abrupt 6-hour exposure

Saltwater and freshwater environments have opposing physiological challenges, yet, there are fish species that are able to enter both habitats during short time spans, and as individuals they must therefore adjust quickly to osmoregulatory contrasts. In this study, we conducted an experiment to test for plastic responses to abrupt salinity changes in two populations of threespine stickleback, Gasterosteus aculeatus, representing two ecotypes (freshwater and ancestral saltwater). We exposed both ecotypes to abrupt native (control treatment) and non-native salinities (0‰ and 30‰) and sampled gill tissue for transcriptomic analyses after 6 h of exposure. To investigate genomic responses to salinity, we analyzed four different comparisons; one for each ecotype (in their control and exposure salinity; (1) and (2), one between ecotypes in their control salinity (3), and the fourth comparison included all transcripts identified in (3) that did not show any expressional changes within ecotype in either the control or the exposed salinity (4)). Abrupt salinity transfer affected the expression of 10 and 1530 transcripts for the saltwater and freshwater ecotype, respectively, and 1314 were differentially expressed between the controls, including 502 that were not affected by salinity within ecotype (fixed expression). In total, these results indicate that factors other than genomic expressional plasticity are important for osmoregulation in stickleback, due to the need for opposite physiological pathways to survive the abrupt change in salinity.

Effects of polystyrene nanoplastics with different functional groups on rice (Oryza sativa L.) seedlings: Combined transcriptome, enzymology, and physiology

Plastic particles originate from the daily use of plastics and have become a new form of pollutant. However, the effects of nanoplastics (NPs) on higher plants remain largely unclear, thus requiring further research. In this study, rice plants were exposed to polystyrene (PS) NPs with different functional groups to determine their toxicity. The presence of NPs reduced the biomass and photosynthetic capacity of rice. Compared with control (CK), the heights of rice plants exposed to no-modified PS, carboxyl-modified PS (PS-COOH) and amino-modified PS (PS-NH₂) groups decreased by 13.59%, 26.61%, and 42.71%, while the dry shoot weight decreased by 47.46%, 50.09%, and 71.04%, respectively. All treatments activated the antioxidant levels of rice and reduced photosynthesis. Transcriptome analysis showed that NPs induced the expression of genes related to antioxidant enzyme activity in rice roots. Rice could partially reduce the xenobiotic toxicity caused by external sources by regulating phenylpropane biosynthesis and the processes involved in cell detoxification. PS mainly affected the process of RNA metabolism, while PS-COOH mainly affected ion transport, and PS-NH₂ mainly affected the synthesis of macromolecular protein, which had different effects on rice growth.

Novel lncRNA lncRNA001074 participates in the low salinity-induced response in the sea cucumber Apostichopus japonicus by targeting the let-7/NKAα axis

Salinity fluctuations have severe impacts on sea cucumbers and therefore important consequences in sea cucumber farming. The responses of sea cucumbers to salinity changes are reflected in the expression profiles of multiple genes and non-coding RNAs (ncRNAs). The microRNA (let-7) which is a developmental regulator, the ion transporter gene sodium potassium ATPase gene (NKAα), and the long ncRNA lncRNA001074 were previously shown to be involved in responses to salinity changes in various marine species. To better understand the relationship between ncRNAs and target genes, the let-7/NKAα/lncRNA001074 predicted interaction was investigated in this study using luciferase reporter assays and gene knockdowns in the sea cucumber Apostichopus japonicus. The results showed that NKAα was the target gene of let-7 and NKAα expression levels were inversely correlated with let-7 expression based on the luciferase reporter assays and western blots. The let-7 abundance was negatively regulated by lncRNA001074 and NKAα both in vitro and in vivo. Knockdown of lncRNA001074 led to let-7 overexpression. These results demonstrated that lncRNA001074 binds to the 3'-UTR binding site of let-7 in a regulatory manner. Furthermore, the expression profiles of let-7, NKAα, and lncRNA001074 were analyzed in sea cucumbers after the knockdown of each of these genes. The results found that lncRNA001074 competitively bound let-7 to suppress NKAα expression under low salinity conditions. The downregulation of let-7, in conjunction with the upregulation of lncRNA001074 and NKAα, may be essential for the response to low salinity change in sea cucumbers. Therefore, the dynamic balance of the lncRNA001074, NKAα, and let-7 network might be a potential response mechanism to salinity change in sea cucumbers.

Identification and Characterization of Osmoregulation Related MicroRNAs in Gills of Hybrid Tilapia Under Three Types of Osmotic Stress

Researchers have increasingly suggested that microRNAs (miRNAs) are small non-coding RNAs that post-transcriptionally regulate gene expression and protein translation in organs and respond to abiotic and biotic stressors. To understand the function of miRNAs in osmotic stress regulation of the gills of hybrid tilapia (Oreochromis mossambicus ♀ × Oreochromis urolepis hornorum ♂), high-throughput Illumina deep sequencing technology was used to investigate the expression profiles of miRNAs under salinity stress (S, 25‰), alkalinity stress (A, 4‰) and salinity-alkalinity stress (SA, S: 15‰, A: 4‰) challenges. The results showed that 31, 41, and 27 upregulated and 33, 42, and 40 downregulated miRNAs (P < 0.05) were identified in the salt stress, alkali stress, and saline-alkali stress group, respectively, which were compared with those in the control group (C). Fourteen significantly differently expressed miRNAs were selected randomly and then validated by a quantitative polymerase chain reaction. On the basis of Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis, genes related to osmoregulation and biosynthesis were enriched in the three types of osmotic stress. In addition, three miRNAs and three predicted target genes were chosen to conduct a quantitative polymerase chain reaction in the hybrid tilapia and its parents during 96-h osmotic stress. Differential expression patterns of miRNAs provided the basis for research data to further investigate the miRNA-modulating networks in osmoregulation of teleost.

MicroRNA-mediated stress response in bivalve species

Bivalve mollusks are important aquatic organisms, which are used for biological monitoring because of their abundance, ubiquitous nature, and abilities to adapt to different environments. MicroRNAs (miRNAs) are small noncoding RNAs, which typically silence the expression of target genes; however, certain miRNAs directly or indirectly upregulate their target genes. They are rapidly modulated and play an essential role in shaping the response of organisms to stresses. Based on the regulatory function and rapid alteration of miRNAs, they could act as biomarkers for biotic and abiotic stress, including environmental stresses and contaminations. Moreover, mollusk, particularly hemocytes, rapidly respond to environmental changes, such as pollution, salinity changes, and desiccation, which makes them an attractive model for this purpose. Thus, bivalve mollusks could be considered a good animal model to examine a system's response to different environmental conditions and stressors. miRNAs have been reported to adjust the adaptation and physiological functions of bivalves during endogenous and environmental stressors. In this review, we aimed to discuss the potential mechanisms underlying the response of bivalves to stressors and how miRNAs orchestrate this process; however, if necessary, other organisms' response is included to explain specific processes.

Quantitative Trait Loci Mapping and Marker Identification for Low Salinity Tolerance Trait in the Swimming Crab (Portunus trituberculatus)

Low salinity is one of the most important abiotic factors that directly affect the abundance of the swimming crab, Portunus trituberculatus. Quantitative trait loci (QTL) mapping could be helpful in identifying the markers and genes involved in low salinity tolerance. In this study, two QTLs of low salt tolerance were mapped on linkage group 17 (LG17, 2.6-5.2 cM) based on a high-density linkage map. Ninety-five markers related to low salinity tolerance were identified via association analysis, and seventy-nine low salt-related candidate genes (including ammonium transport, aldehyde dehydrogenase, and glucosyltransferase) were screened from draft genome of the species via these markers. This represents the first report of QTL mapping for low salinity tolerance in the swimming crab, which may be useful to elucidate salinity adaptation mechanisms.

Identification and profiling of microRNAs responsive to cadmium toxicity in hepatopancreas of the freshwater crab Sinopotamon henanense

Background

Cadmium (Cd) is a ubiquitous environmental toxicant for aquatic animals. The freshwater crab, Sinopotamon henanense (S. henanense), is a useful model for monitoring Cd exposure since it is widely distributed in sediments whereby it tends to accumulate several toxicants, including Cd. In the recent years, the toxic effects of Cd in the hepatopancreas of S. henanense have been demonstrated by a series of biochemical analysis and ultrastructural observations as well as the deep sequencing approaches and gene expression profile analysis. However, the post-transcriptional regulatory network underlying Cd toxicity in S.henanense is still largely unknown.

Results

The miRNA transcriptional profile of the hepatopancreas of S. henanense was used to investigate the expression levels of miRNAs in response to Cd toxicity. In total, 464 known miRNAs and 191 novel miRNAs were identified. Among these 656 miRNAs, 126 known miRNAs could be matched with the miRNAs of Portunus trituberculatus, Eriocheir sinensis and Scylla paramamosain. Furthermore, a total of 24 conserved miRNAs were detected in these four crab species. Fifty-one differentially expressed miRNAs were identified in the Cd-exposed group, with 31 up-regulated and 20 down-regulated. Eight of the differentially expressed miRNAs were randomly selected and verified by the quantitative real-time PCR (qRT-PCR), and there was a general consistency (87.25%) between the qRT-PCR and miRNA transcriptome data. A total of 5258 target genes were screened by bioinformatics prediction. GO term analysis showed that, 17 GO terms were significantly enriched, which were mainly related to the regulation of oxidoreductase activity. KEGG pathway analysis showed that 18 pathways were significantly enriched, which were mainly associated with the biosynthesis, modification and degradation of proteins.

Conclusion

In response to Cd toxicity, in the hepatopancreas of S. henanense, the expressions of significant amount of miRNAs were altered, which may be an adaptation to resist the oxidative stress induced by Cd. These results provide a basis for further studies of miRNA-mediated functional adaptation of the animal to combat Cd toxicity.

Salinity stress-induced differentially expressed miRNAs and target genes in sea cucumbers Apostichopus japonicus

Environmental salinity is an important abiotic factor influencing normal physiological functions and productive performance in the sea cucumber Apostichopus japonicus. It is therefore important to understand how changes in salinity affect sea cucumbers in the face of global climate change. In this study, we investigated the responses to salinity stress in sea cucumbers using mRNA and miRNA sequencing. The regulatory network of mRNAs and miRNAs involved in salinity stress was examined, and the metabolic pathways enriched for differentially expressed miRNAs and target mRNAs were identified. The top 20 pathways were involved in carbohydrate metabolism, fatty acid metabolism, degradation, and elongation, amino acid metabolism, genetic information processing, metabolism of cofactors and vitamins, transport and catabolism, and environmental information processing. A total of 22 miRNAs showed differential expression during salinity acclimation. The predicted 134 target genes were enriched in functions consistent with the results of gene enrichment based on transcriptome analysis. These results suggested that sea cucumbers deal with salinity stress via changes in amino acid metabolism, ion channels, transporters, and aquaporins, under stimulation by environmental signals, and that this process requires energy from carbohydrate and fatty acid metabolism. Salinity challenge also induced miRNA expression. These results provide a valuable genomic resource that extends our understanding of the unique biological characteristics of this economically important species under conditions of salinity stress.

Genome-wide identification of Chinese shrimp (Fenneropenaeus chinensis) microRNA responsive to low pH stress by deep sequencing

pH has a great impact on the distribution, growth, behavior, and physiology in many aquatic animals. Here, we analyzed miRNA expression profiles of Chinese shrimp (Fenneropenaeus chinensis) from control pH (8.2) and low pH (6.5)-treated shrimp. Expression analysis identified 6 known miRNAs and 23 novel miRNAs with significantly different expression between control pH 8.2 and low pH 6.5; the predicted target genes of differentially expressed miRNAs were significantly enriched in organic acid metabolic process, oxidoreductase activity, coenzyme binding, cofactor binding, and collagen trimer. Moreover, target genes were significantly enriched in several Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways including citrate cycle, pyruvate metabolism, cytokine-cytokine receptor interaction, tight junction, carbon metabolism, etc. Our survey expanded the number of known shrimp miRNAs and provided comprehensive information about miRNA response to low pH stress.

Potential role for microRNA in facilitating physiological adaptation to hypoxia in the Pacific whiteleg shrimp Litopenaeus vannamei

Hypoxia is one of the most common physiological stressors in shrimp farming. Post-transcriptional regulation by microRNAs has been recognized as a ubiquitous strategy to enable transient phenotypic plasticity and adaptation to stressful environment, but involvement of microRNAs in hypoxia stress response of penaeid shrimp remains elusive. In this study, small RNA sequencing and comparative transcriptomic analysis was conducted to construct a comprehensive microRNA dataset for the whiteleg shrimp Litopenaeus vannamei exposed to hypoxia challenge. A total of 3324 known miRNAs and 8 putative novel miRNAs were identified, providing a valuable resource for future investigation on the functional mechanism of miRNAs in shrimp. Upon hypoxia, 1213 miRNAs showed significant differential expression, and many well-known miRNAs involved in hypoxia tolerance such as miR-210, let-7, miR-143 and miR-101 were identified. Remarkably, the vast majority of these miRNAs were up-regulated, suggesting that up-regulation of miRNAs may represent an effective strategy to inhibit protein translation under stressful hypoxic condition. The differentially expressed miRNAs were potentially targeting a wide variety of genes, including those with essential roles in hypoxia tolerance such as HIF1a and p53. GO and KEGG enrichment analysis further revealed that a broad range of biological processes and metabolic pathways were over-represented. Several GO terms associated with gene transcription and translation and KEGG pathways related to cytoskeleton remodeling, immune defense and signaling transduction were enriched, highlighting the crucial roles of these cellular events in the adaptation to hypoxia. Taken together, our study revealed that the differentially expressed miRNAs may regulate host response to hypoxia by modulating the expression of stress response genes such as HIF1a and p53 and affecting key cellular events involved in hypoxia adaptation. The findings would expand our knowledge of the biochemical and molecular underpinnings of hypoxia response strategies used by penaeid shrimp, and contribute to a better understanding of the molecular mechanisms of hypoxia tolerance in decapod crustaceans.

Transcriptional responses to low-salinity stress in the gills of adult female Portunus trituberculatus

The swimming crab (Portunus trituberculatus, Portunus) can tolerate low salinity, but the mechanism of its varied salinity adaptation at the molecular level remains unclear. In this study, we prepared four mRNA and microRNA (miRNA) libraries using the gills collected from four salinity groups and performed RNA-sequencing (RNA-Seq) to identify the genes related to the low salinity. We set 25 ppt as the control group. A total of 659 genes were differentially expressed in at least one of the six comparison groups (25 ppt vs. 20 ppt, 25 ppt vs. 15 ppt, 25 ppt vs. 10 ppt, 20 ppt vs. 15 ppt, 20 ppt vs. 10 ppt and 15 ppt vs. 10 ppt). A total of 15 and 9 unigenes were downregulated and upregulated under low salinity compared with that in 25 ppt, respectively. Six genes, namely, aminopeptidase, centromere protein, cytochrome b5 reductase, bone morphogenetic protein, and two carbonic anhydrases, were selected for verification through quantitative real-time PCR. The results were consistent with the RNA-Seq results. Furthermore, 95 conserved miRNAs and 16 novel miRNAs were differentially expressed in at least one of the six comparison groups. Analysis of the miRNA-mRNA interaction showed that miR-2 and miR-317 regulated >50 mRNA targets. In addition, let-7c was downregulated in all groups under low salinity compared with that in the control group. This study helped elucidate the adaptation mechanism of the swimming crab in low-saline environments.

Identification and characterization of miRNAs in the gills of the mud crab (Scylla paramamosain) in response to a sudden drop in salinity

BACKGROUND

The mud crab (Scylla paramamosain) is a euryhaline and commercially important species. MiRNAs participate in the regulation of many physiological activities.

RESULTS

The miRNA transcriptome of the gills of S. paramamosain was used to investigate the expression profiles of miRNAs in response to a sudden drop in salinity. In total, seven known miRNAs and 43 novel miRNAs were identified, with 18 differentially expressed small RNAs. Fourteen thousand nine hundred fifty-one differentially expressed miRNAs target genes were screened by prediction. GO analysis of differentially expressed miRNAs target genes indicated that 578 genes associated with cellular processes, 523 associated with metabolic processes, and 422 associated with single-organism processes were the most strongly affected by a sudden drop in salinity from 23‰ to 3‰. KEGG pathway analysis showed 14 pathways were related to amino acid metabolism, which plays an important role in osmoregulation. Besides, several pathways were associated with starch and sucrose metabolism (ko00500), glycosaminoglycan degradation (ko00531), and galactose metabolism (ko00052).

CONCLUSIONS

S. paramamosain regulated osmotic pressure and energy balance by regulating target genes to adapt to a sudden changes in salinity. These results provided a basis for further investigations of miRNA-modulating networks underlying the osmoregulation of S. paramamosain.

Identification and profiling of growth-related microRNAs of the swimming crab Portunus trituberculatus by using Solexa deep sequencing

MicroRNAs (miRNAs) are a class of endogenous small non-coding RNAs that regulate gene expression by post-transcriptional repression of mRNAs. The swimming crab Portunus trituberculatus is one of the most important crustacean species for aquaculture in China. However, to date no miRNAs have been reported to for modulating growth in P. trituberculatus. To investigate miRNAs involved in the growth of this species, we constructed six small RNA libraries for big individuals (BIs) and small individuals (SIs) from a highly inbred family. Six mixed RNA pools of five tissues (eyestalk, gill, heart, hepatopancreas, and muscle) were obtained. By aligning sequencing data with those for known miRNAs, a total of 404 miRNAs, including 339 known and 65 novel miRNAs, were identified from the six libraries. MiR-100 and miR-276a-3p were among the most prominent miRNA species. We identified seven differentially expressed miRNAs between the BIs and SIs, which were validated using real-time PCR. Preliminary analyzes of their putative target genes and GO and KEGG pathway analyzes showed that these differentially expressed miRNAs could play important roles in global transcriptional depression and cell differentiation of P. trituberculatus. This study reveals the first miRNA profile related to the body growth of P. trituberculatus, which would be particularly useful for crab breeding programs.