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

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.

Sea cucumber physiological response to abiotic stress: Emergent contaminants and climate change

The ocean is facing a multitude of abiotic stresses due to factors such as climate change and pollution. Understanding how organisms in the ocean respond to these global changes is vital to better predicting consequences. Sea cucumbers are popular echinoderms with multiple ecological, nutritional, and pharmaceutical benefits. Here, we reviewed the effects of environmental change on an ecologically important echinoderm of the ocean, aiming to understand their response better, which could facilitate healthy culture programs under environmental changes and draw attention to knowledge gaps. After screening articles from the databases, 142 studies were included on the influence of emergent contaminants and climate variation on the early developmental stages and adults of sea cucumbers. We outlined the potential mechanism underlying the physiological response of sea cucumbers to emerging contaminants and climate change. It can be concluded that the physiological response of sea cucumbers to emergent contaminants differs from their response to climate change. Sea cucumbers could accumulate pollutants in their organs but are aestivated when exposed to extreme climate change. Research showed that the physiological response of sea cucumbers to pollutants indicates that these pollutants impair critical physiological processes, particularly during the more susceptible early phases of development compared to adults, and the accumulation of these pollutants in adults is often observed. For climate change, sea cucumbers showed gradual adaptation to the slight variation. However, sea cucumbers undergo aestivation under extreme conditions. Based on this review, critical suggestions for future research are presented, and we call for more efforts focusing on the co-occurrence of different stressors to extend the knowledge regarding the effects of environmental changes on these economically and ecologically important species.

Genome-wide analysis and expression profile of TCP gene family under drought and salinity stress condition in cowpea (Vigna unguiculata (L.) Walp.)

TCP transcription factors are known to regulate abiotic stress condition, but their role in V. unguiculata remains unexplored. So, in silico analysis and expression profile of the TCP gene family were performed in V. unguiculata to understand its role in response to heat and drought stress. A genome-wide search detected 28 TCPs (designated as VuTCPs) that were grouped into three subclasses by phylogenetic analysis. Gene structure, synteny, and phylogeny analyses of VuTCPs have shown a typical evolutionary path. One tandem and eight segmental duplication events were identified. Furthermore, identified duplicated, and orthologous VuTCP genes were under strong purifying selection pressure. A total of 15 SSRs were identified in the 12 VuTCPs, while 10 VuTCP genes were regulated by different miRNAs having a major role in abiotic stress tolerance. Analysed physicochemical properties, cis-acting elements, and gene ontology suggested that VuTCPs play various roles, including salinity and drought stress tolerance. qRT-PCR analysis showed that 11 and 15 VuTCPs were upregulated under drought and salinity stress conditions, respectively. Our findings provide comprehensive insights into the genomic characterization of the VuTCPs gene family in V. unguiculata, offering a foundation for understanding their structure, evolution, and role in abiotic stress tolerance.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-03976-x.

lncRNA-miRNA-mRNA network in kidney transcriptome of Labeo rohita under hypersaline environment

The present study describes the kidney transcriptome of Labeo rohita, a freshwater fish, exposed to gradually increased salinity concentrations (2, 4, 6 and 8ppt). A total of 10.25 Gbps data was generated, and a suite of bioinformatics tools, including FEELnc, CPC2 and BLASTn were employed for identification of long non-coding RNAs (lncRNAs) and micro RNAs (miRNAs). Our analysis revealed a total of 170, 118, 99, and 269 differentially expressed lncRNA and 120, 118, 99, and 124 differentially expressed miRNAs in 2, 4, 6 and 8 ppt treatment groups respectively. Two competing endogenous RNA (ceRNA) networks were constructed i.e. A* ceRNA network with up-regulated lncRNAs and mRNAs, down-regulated miRNAs; and B* ceRNA network vice versa. 2ppt group had 131 and 83 lncRNA-miRNA-mRNA pairs in A* and B* networks, respectively. 4ppt group featured 163 pairs in A* network and 191 in B* network, while the 6ppt had 103 and 105 pairs. 8ppt group included 192 and 174 pairs. These networks illuminate the intricate RNA interactions in freshwater fish to varying salinity conditions.

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.

Interplay of gene expression and regulators under salinity stress in gill of Labeo rohita

BACKGROUND

Labeo rohita is the most preferred freshwater carp species in India. The concern of increasing salinity concentration in freshwater bodies due to climate change may greatly impact the aquatic environment. Gills are one of the important osmoregulatory organs and have direct contact with external environment. Hence, the current study is conducted to understand the gill transcriptomic response of L. rohita under hypersalinity environment.

RESULTS

Comprehensive analysis of differentially expressed long non-coding RNAs (lncRNAs), microRNAs (miRNAs) and mRNAs was performed in gills of L. rohita treated with 2, 4, 6 and 8ppt salinity concentrations. Networks of lncRNA-miRNA-mRNA revealed involvement of 20, 33, 52 and 61 differentially expressed lncRNAs, 11, 13, 26 and 21 differentially expressed miRNAs in 2, 4, 6 and 8ppt groups between control and treatment respectively. These lncRNA-miRNA pairs were regulating 87, 214, 499 and 435 differentially expressed mRNAs (DE mRNAs) in 2, 4, 6 and 8ppt treatments respectively. Functional analysis of these genes showed enrichment in pathways related to ion transportation and osmolyte production to cope with induced osmotic pressure due to high salt concentration. Pathways related to signal transduction (MAPK, FOXO and phosphatidylinositol signaling), and environmental information processing were also upregulated under hypersalinity. Energy metabolism and innate immune response pathways also appear to be regulated. Protein turnover was high at 8ppt as evidenced by enrichment of the proteasome and aminoacyl tRNA synthesis pathways, along with other enriched KEGG terms such as apoptosis, cellular senescence and cell cycle.

CONCLUSION

Altogether, the RNA-seq analysis provided valuable insights into competitive endogenous (lncRNA-miRNA-mRNA) regulatory network of L. rohita under salinity stress. L. rohita is adapting to the salinity stress by means of upregulating protein turnover, osmolyte production and removing the damaged cells using apoptotic pathway and regulating the cell growth and hence diverting the essential energy for coping with salinity stress.

Transcriptome Analysis of Marbled Rockfish Sebastiscus marmoratus under Salinity Stress

The marbled rockfish, Sebastiscus marmoratus, belongs to the euryhaline fish and is an oviparous scleractinian fish. There are few studies on the adaptation mechanism, functional genes, and related pathways of S. marmoratus and salinity. The results showed that a total of 72.1 GB of clean reads were obtained and all clean reads annotated a total of 25,278 Unigenes, of which 2,160 were novel genes. Compared to 20‱, 479 and 520 differential genes were obtained for 35‱ and 10‱, respectively. Gene ontology (GO) enrichment analysis revealed significant enrichment in protein binding, ion binding, ATP binding, and catalytic activity. Kyoto Encyclopedia of Genes and Genomes (KEGG) showed that differentially expressed genes significantly expressed under salinity stress were mainly involved in the pathways of the cytochrome P450 metabolism of xenobiotics, tryptophan metabolism, cellular senescence, and calcium signaling pathways. Among them, pik3r6b, cPLA2γ-like, and WSB1 were differentially expressed in all three groups, and they were associated with apoptosis, inflammation, DNA damage, immune regulation, and other physiological processes. Six differentially expressed genes were randomly selected for qRT-PCR validation, and the results showed that the transcriptomic data were of high confidence.

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.

The Effect of Salinity Stress on Enzyme Activities, Histology, and Transcriptome of Silver Carp (Hypophthalmichthys molitrix)

A 56-day study was performed to examine the effect of freshwater (FW) and brackish water (BW 6‱ salinity) on the antioxidant ability, Na+/K+-ATPase (NKA) activities, histology, and transcriptome of the gill and kidney tissue in juvenile silver carp (Hypophthalmichthys molitrix). The results show that when juvenile silver carp were exposed to 6‱ salinity, the activities of superoxide dismutase (SOD) and catalase (CAT) were shown to be substantially increased (p < 0.05), while glutathione peroxidase (GSH-PX) activities in gill were not significantly affected (p < 0.05). In kidney tissue, SOD, CAT, and GSH-PX, enzyme activities peaked at 24, 8, and 4 h, respectively, but were not significantly different compared with the control group (p < 0.05). In addition, significant effects of salinity were observed for the NKA level in both the gills and kidney tissues (p < 0.05). The gill filaments of juvenile silver carp under the BW group all underwent adverse changes within 72 h, such as cracks and ruptures in the main part of the gill filaments, bending of the gill lamellae and enlargement of the gaps, and an increase in the number of mucus and chloride-secreting cells. Transcriptome sequencing showed 171 and 261 genes in the gill and kidney tissues of juvenile silver carp compared to the BW group, respectively. Based on their gene ontology annotations, transcripts were sorted into four functional gene groups, each of which may play a role in salt tolerance. Systems involved in these processes include metabolism, signal transduction, immunoinflammatory response, and ion transport. The above findings indicate that the regulation processes in juvenile silver carp under brackish water conditions are complex and multifaceted. These processes and mechanisms shed light on the regulatory mechanism of silver carp osmolarity and provide a theoretical foundation for future research into silver carp growth in brackish water aquaculture area.

MicroRNA-139-5p acts as a suppressor gene for depression by targeting nuclear receptor subfamily 3, group C, member 1

MicroRNA-139-5p (miR-139-5p) is one of the most differentially expressed miRNAs in the brain between healthy people and depressed patients. However, its function in depression is unclear. Therefore, we investigated the function of miR-139-5p in depression. Here, miR-139-5p expression was found to be upregulated in the model group. MiR-139-5p inhibition could increase sucrose preference and decrease mice immobility time after chronic corticosterone (CORT) injection. Furthermore, compared with the antago-NC group, 3 weeks of antagomiR-139-5p treatment significantly decreased miR-139-5p level in model group hippocampus, increased sucrose preference index, reduced neuron damages, and enhanced the levels of nuclear receptor subfamily 3 group C member 1 (NR3C1), brain-derived neurotrophic factor (BDNF), phosphorylated/total tyrosine kinase receptor B (p-TrkB/TrkB), phosphorylated/total cAMP-response element-binding protein (p-CREB/CREB) and phosphorylated/total extracellular regulated protein kinases (p-ERK/ERK). Moreover, as a potential target for miR-139-5p, NR3C1 level was reduced by miR-139-5p mimic. Altogether, by activating the BDNF-TrkB signaling pathway, miR-139-5p inhibition plays an antidepressant-like role and might serve as an effective depression target (Fig. graphical abstract).

Inhibition of microRNA-29b suppresses oxidative stress and reduces apoptosis in ischemic stroke

MicroRNAs (miRNAs) regulate protein expression by antagonizing the translation of mRNAs and are effective regulators of normal nervous system development, function, and disease. MicroRNA-29b (miR-29b) plays a broad and critical role in brain homeostasis. In this study, we tested the function of miR-29b in animal and cell models by inhibiting miR-29b expression. Mouse models of middle cerebral artery occlusion were established using the modified Zea-Longa suture method. Prior to modeling, 50 nmol/kg miR-29b antagomir was injected via the tail vein. MiR-29b expression was found to be abnormally increased in ischemic brain tissue. The inhibition of miR-29b expression decreased the neurological function score and reduced the cerebral infarction volume and cell apoptosis. In addition, the inhibition of miR-29b significantly decreased the malondialdehyde level, increased superoxide dismutase activity, and Bcl-2 expression, and inhibited Bax and Caspase3 expression. PC12 cells were treated with glutamate for 12 hours to establish in vitro cell models of ischemic stroke and then treated with the miR-29 antagomir for 48 hours. The results revealed that miR-29b inhibition in PC12 cells increased Bcl-2 expression and inhibited cell apoptosis and oxidative damage. These findings suggest that the inhibition of miR-29b inhibits oxidative stress and cell apoptosis in ischemic stroke, producing therapeutic effects in ischemic stroke. This study was approved by the Laboratory Animal Care and Use Committee of the First Affiliated Hospital of Zhengzhou University (approval No. 201709276S) on September 27, 2017.

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.

Clustering genomic organization of sea cucumber miRNAs impacts their evolution and expression

Echinoderms are marine deuterostomes with fascinating adaptation features such as aestivation and organ regeneration. However, post-transcriptional gene regulation by microRNAs (miRNAs) underlying these features are largely unexplored. Here, using homology-based and de novo approaches supported by expression data, we provided a comprehensive annotation of miRNA genes in the sea cucumber Apostichopus japonicus. By linkage and phylogenic analyses, we characterized miRNA genomic organization, evolutionary history and expression regulation. The results showed that sea cucumbers evolved a large number of new miRNAs, which tended to form polycistronic clusters via tandem duplication that had been especially active in the echinoderms. Most new miRNAs were weakly expressed, but miRNA clustering increased the expression level of clustered new miRNAs. The most abundantly expressed new miRNAs were organized in a single tandem cluster (cluster n2), which was activated during aestivation and intestine regeneration. Overall, our analyses suggest that clustering of miRNAs is important for their evolutionary origin, expression control, and functional cooperation.

Sea cucumbers in a high temperature and low dissolved oxygen world: Roles of miRNAs in the regulation of environmental stresses

The exacerbation of global warming has driven changes in environmental factors, including water temperature and oxygen concentration. The sea cucumber Apostichopus japonicus, an economically important aquatic animal, is constantly and directly challenged by heat and hypoxia. In this study, 12 small RNA libraries were constructed for this species, and a total of 21, 26 and 22 differentially expressed (DE) miRNAs were clarified in A. japonicus under thermal (26 °C), hypoxic (2 mg/L) and the combined stresses. Comparative miRNA sequencing analysis and real-time PCR were used to identify and validate the representative miRNAs, including Aja-miR-novel-299, Aja-let-7b-3p, Aja-miR-71b-5p, Aja-miR-novel-13218 and Aja-miR-2004 in response to high temperature, and Aja-miR-92b-3p, Aja-miR-210-5p and Aja-miR-novel-26331 in response to oxygen limitation. GO and KEGG pathway analysis revealed that the potential target genes of DE-miRNAs involved in biosynthesis, metabolism, immunity, cell growth and death, translation and signaling transduction. Key DE-miRNAs with potentially targeted genes associated with heat shock and hypoxia response were also determined. These results may help explaining the role of miRNA regulation in stress resistance, as well as the potential molecular regulation mechanism of the echinoderm A. japonicus in the context of global warming.

Integrated small RNA and Degradome sequencing provide insights into salt tolerance in sesame (Sesamum indicum L.)

BACKGROUND

MicroRNAs (miRNAs) exhibit important regulatory roles in the response to abiotic stresses by post-transcriptionally regulating the target gene expression in plants. However, their functions in sesame response to salt stress are poorly known. To dissect the complex mechanisms underlying salt stress response in sesame, miRNAs and their targets were identified from two contrasting sesame genotypes by a combined analysis of small RNAs and degradome sequencing.

RESULTS

A total of 351 previously known and 91 novel miRNAs were identified from 18 sesame libraries. Comparison of miRNA expressions between salt-treated and control groups revealed that 116 miRNAs were involved in salt stress response. Using degradome sequencing, potential target genes for some miRNAs were also identified. The combined analysis of all the differentially expressed miRNAs and their targets identified miRNA-mRNA regulatory networks and 21 miRNA-mRNA interaction pairs that exhibited contrasting expressions in sesame under salt stress.

CONCLUSIONS

This comprehensive integrated analysis may provide new insights into the genetic regulation mechanism of miRNAs underlying the adaptation of sesame to salt stress.