Characterization of lncRNAs involved in cold acclimation of zebrafish ZF4 cells

MSTRG3207 promotes apoptosis in zebrafish ZF4 cells via sponging dre-miR-736/bbc3/LOC101885512 axis during cold acclimation

Long non-coding RNAs (lncRNAs) play essential roles in a variety of biological processes. It has been recently reported that lncRNAs can regulate mRNA expression by binding to microRNAs (miRNAs) as competing endogenous RNAs (ceRNAs). However, the involvement of this regulatory mechanism during cold acclimation in fish remains unclear. In this study, we constructed a ceRNA network mediated by lncRNAs in cold-acclimated zebrafish ZF4 cells through bioinformatic analysis of the mRNA, miRNA, and lncRNA profiles obtained from ZF4 cells cultured at 18 °C for 30 days. A previously uncharacterized lncRNA, MSTRG3207, was selected for further analysis. MSTRG3207 was upregulated and dre-miR-736 was downregulated during cold acclimation. MSTRG3207 was cloned by rapid amplification of cDNA ends (RACE) and functionally characterized. The binding of MSTRG3207 to dre-miR-736 was validated by dual-luciferase reporter assay. Under cold acclimation, MSTRG3207 promoted apoptosis by sponging dre-miR-736 and upregulating bbc3 and LOC101885512, two apoptotic genes targeted by dre-miR-736. Taken together, our findings indicate that MSTRG3207 upregulation promotes apoptosis by sponging dre-miR-736 during cold acclimation in fish.

LIPRNAseq: a method to discover lipid interacting RNAs by sequencing

BACKGROUND

Current biological research extensively describes the interactions of molecules such as RNA with other nucleic acids or proteins. However, the relatively recent discovery of nuclear phospholipids playing biologically relevant processes outside membranes, as well as, RNA-lipid interactions shows the need for new methods to explore the identity of these RNAs.

METHODS AND RESULTS

In this study, we describe the method for LIPID-RNA isolation followed by sequencing and analysis of the RNA that has the ability to interact with the selected lipids. Here we utilized specific phospholipid coated beads for selective RNA binding. We tested RNA from organisms belonging to different realms (human, plant, and yeast), and tested their ability to bind a specific lipid.

CONCLUSIONS

The results show several RNAs differentially enriched in the pull-down of phosphatidyl Inositol 4,5 bisphosphate coated beads. This method is helpful to screen lipid-binding RNA, which may have relevant biological functions. The method can be used with different lipids and comparison of pull-downs and can narrow the selection of RNAs that interact with a particular lipid for further studies.

Transcriptome sequencing reveals the differentially expressed lncRNAs and mRNAs in response to cold acclimation and cold stress in Pomacea canaliculata

Background

Tolerance of low temperature has a significant impact on survival and expansion of invasive snail Pomacea canalicuata. Cold acclimation can enhance cold tolerance of Pomacea canalicuata. To elucidate the molecular mechanism of P. canaliculata’s responses to cold acclimation and cold stress, a high-throughput transcriptome analysis of P. canaliculata was performed, and gene expression following artificial cold acclimation and then cold stress at 0 °C for 24 h was compared using RNA sequencing.

Results

Using the Illumina platform, we obtained 151.59 G subreads. A total of 5,416 novel lncRNAs were identified, and 3166 differentially expressed mRNAs and 211 differentially expressed lncRNAs were screened with stringent thresholds. The potential antisense, cis and trans targets of lncRNAs were predicted. Kyoto Encyclopedia of Genes and Genomes enrichment analysis showed that many target genes were involved in proteasome, linoleic acid metabolism and retinol metabolism under cold acclimation. The lncRNA of P. canaliculata could participate in cold acclimation by regulating the expression of E3 ubiquitin protein ligase, 26S proteasome non-ATPase dependent regulation subunit, glutathione S-transferase, sodium/glucose cotransporter and cytochrome P450.

Conclusions

These results broaden our understanding of cold acclimation and cold stress associated lncRNAs and mRNAs, and provide new insights into lncRNA mediated regulation of P. canaliculata cold acclimation and cold stress response.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08622-5.

Integrative Analysis of lncRNA-mRNA Coexpression in Human Lung Epithelial Cells Exposed to Dimethyl Selenide-Derived Secondary Organic Aerosols

Dimethyl selenide (DMSe) is one of the major volatile organoselenium compounds released into the atmosphere through plant metabolism and microbial methylation. DMSe has been recently revealed as a precursor of secondary organic aerosol (SOA), and its resultant SOA possesses strong oxidizing capability toward thiol groups that can perturb several major biological pathways in human airway epithelial cells and is linked to genotoxicity, DNA damage, and p53-mediated stress responses. Mounting evidence has suggested that long noncoding RNAs (lncRNAs) are involved in stress responses to internal and environmental stimuli. However, the underlying molecular interactions remain to be elucidated. In this study, we performed integrative analyses of lncRNA-mRNA coexpression in the transformed human bronchial epithelial BEAS-2B cell line exposed to DMSe-derived SOA. We identified a total of 971 differentially expressed lncRNAs in BEAS-2B cells exposed to SOA derived from O₃ and OH oxidation of DMSe. Gene ontology (GO) network analysis of cis-targeted genes showed significant enrichment of DNA damage, apoptosis, and p53-mediated stress response pathways. trans-Acting lncRNAs, including PINCR, PICART1, DLGAP1-AS2, and LINC01629, known to be associated with human carcinogenesis, also showed altered expression in cell treated with DMSe-SOA. Overall, this study highlights the regulatory role of lncRNAs in altered gene expression induced by DMSe-SOA exposure.

Pan-tissue transcriptome analysis of long noncoding RNAs in the American beaver Castor canadensis

BACKGROUND

Long noncoding RNAs (lncRNAs) have roles in gene regulation, epigenetics, and molecular scaffolding and it is hypothesized that they underlie some mammalian evolutionary adaptations. However, for many mammalian species, the absence of a genome assembly precludes the comprehensive identification of lncRNAs. The genome of the American beaver (Castor canadensis) has recently been sequenced, setting the stage for the systematic identification of beaver lncRNAs and the characterization of their expression in various tissues. The objective of this study was to discover and profile polyadenylated lncRNAs in the beaver using high-throughput short-read sequencing of RNA from sixteen beaver tissues and to annotate the resulting lncRNAs based on their potential for orthology with known lncRNAs in other species.

RESULTS

Using de novo transcriptome assembly, we found 9528 potential lncRNA contigs and 187 high-confidence lncRNA contigs. Of the high-confidence lncRNA contigs, 147 have no known orthologs (and thus are putative novel lncRNAs) and 40 have mammalian orthologs. The novel lncRNAs mapped to the Oregon State University (OSU) reference beaver genome with greater than 90% sequence identity. While the novel lncRNAs were on average shorter than their annotated counterparts, they were similar to the annotated lncRNAs in terms of the relationships between contig length and minimum free energy (MFE) and between coverage and contig length. We identified beaver orthologs of known lncRNAs such as XIST, MEG3, TINCR, and NIPBL-DT. We profiled the expression of the 187 high-confidence lncRNAs across 16 beaver tissues (whole blood, brain, lung, liver, heart, stomach, intestine, skeletal muscle, kidney, spleen, ovary, placenta, castor gland, tail, toe-webbing, and tongue) and identified both tissue-specific and ubiquitous lncRNAs.

CONCLUSIONS

To our knowledge this is the first report of systematic identification of lncRNAs and their expression atlas in beaver. LncRNAs-both novel and those with known orthologs-are expressed in each of the beaver tissues that we analyzed. For some beaver lncRNAs with known orthologs, the tissue-specific expression patterns were phylogenetically conserved. The lncRNA sequence data files and raw sequence files are available via the web supplement and the NCBI Sequence Read Archive, respectively.

Identification and Evaluation of Long Noncoding RNAs in Response to Handling Stress in Red Cusk-Eel (Genypterus chilensis) via RNA-seq

The red cusk-eel (Genypterus chilensis) is a native species with strong potential to support Chilean aquaculture diversification. Under commercial conditions, fish are exposed to several stressors. To date, little is known about the mechanism involved in the stress response of red cusk-eel, and there is no information related to the regulation mediated by long noncoding RNAs (lncRNAs). The objective of this work was to identify for the first time the lncRNAs in the transcriptome of G. chilensis and to evaluate the differential expression levels of lncRNAs in the liver, head kidney, and skeletal muscle in response to handling stress. We used previously published transcriptome data to identify the lncRNAs by applying a series of filters based on annotation information in several databases to discard coding sequences. We identified a total of 14,614 putative lncRNAs in the transcriptome of red cusk-eel, providing a useful lncRNA reference resource to be used in future studies. We evaluated their differential expression in response to handling stress in the liver, head kidney, and skeletal muscle, identifying 112, 323, and 108 differentially expressed lncRNAs, respectively. The results suggest that handling stress in red cusk-eel generate an altered metabolic status in liver, altered immune response in head kidney, and skeletal muscle atrophy through an important coding and noncoding gene network. This is the first study that identifies lncRNAs in Genypterus genus and that evaluates the relation between handling stress and lncRNAs in teleost fish, thereby providing valuable information regarding noncoding responses to stress in Genypterus species.

Identification and characterization of miRNAs involved in cold acclimation of zebrafish ZF4 cells

MicroRNAs (miRNAs) play vital roles in various biological processes under multiple stress conditions by leading to mRNA cleavage or translational repression. However, the detailed roles of miRNAs in cold acclimation in fish are still unclear. In the present study, high-throughput sequencing was performed to identify miRNAs from 6 small RNA libraries from the zebrafish embryonic fibroblast ZF4 cells under control (28°C, 30 days) and cold-acclimation (18°C, 30 days) conditions. A total of 414 miRNAs, 349 known and 65 novel, were identified. Among those miRNAs, 24 (19 known and 5 novel) were up-regulated, and 23 (9 known and 14 novel) were down-regulated in cold acclimated cells. The Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) enrichment analyses indicated that the target genes of known differentially expressed miRNAs (DE-miRNA) are involved in cold acclimation by regulation of phosphorylation, cell junction, intracellular signal transduction, ECM-receptor interaction and so on. Moreover, both miR-100-3p inhibitor and miR-16b mimics could protect ZF4 cells under cold stress, indicating the involvement of miRNA in cold acclimation. Further study showed that miR-100-3p and miR-16b could regulate inversely the expression of their target gene (atad5a, cyp2ae1, lamp1, rilp, atxn7, tnika, btbd9), and that overexpression of miR-100-3p disturbed the early embryonic development of zebrafish. In summary, the present data show that miRNAs are closely involved in cold acclimation in zebrafish ZF4 cells and provide information for further understanding of the roles of miRNAs in cold acclimation in fish.

Comparative modulation of lncRNAs in wild-type and rag1-heterozygous mutant zebrafish exposed to immune challenge with spring viraemia of carp virus (SVCV)

Although the modulation of immune-related genes after viral infection has been widely described in vertebrates, the potential implications of non-coding RNAs (ncRNAs), especially long non-coding RNAs (lncRNAs), in immunity are still a nascent research field. The model species zebrafish could serve as a useful organism for studying the functionality of lncRNAs due to the numerous advantages of this teleost, including the existence of numerous mutant lines. In this work, we conducted a whole-transcriptome analysis of wild-type (WT) and heterozygous rag1 mutant (rag1^+/-) zebrafish after infection with the pathogen spring viraemia of carp virus (SVCV). WT and rag1^+/- zebrafish were infected with SVCV for 24 h. Kidney samples were sampled from infected and uninfected fish for transcriptome sequencing. From a total of 198,540 contigs, 12,165 putative lncRNAs were identified in zebrafish. Most of the putative lncRNAs were shared by the two zebrafish lines. However, by comparing the lncRNA profiles induced after SVCV infection in WT and rag1^+/- fish, most of the lncRNAs that were significantly induced after viral challenge were exclusive to each line, reflecting a highly differential response to the virus. Analysis of the neighboring genes of lncRNAs that were exclusively modulated in WT revealed high representation of metabolism-related terms, whereas those from rag1^+/- fish showed enrichment in terms related to the adaptive immune response, among others. On the other hand, genes involved in numerous antiviral processes surrounded commonly modulated lncRNAs, as expected. These results clearly indicate that after SVCV infection in zebrafish, the expression of an array of lncRNAs with functions in different aspects of immunity is induced.