Screening and functional identification of lncRNAs under β-diketone antibiotic exposure to zebrafish (Danio rerio) using high-throughput sequencing

Genome-Wide Identification and Involvement in Response to Biotic and Abiotic Stresses of lncRNAs in Turbot (Scophthalmus maximus)

Long non-coding RNAs (lncRNAs) play crucial roles in a variety of biological processes, including stress response. However, the number, characteristics and stress-related expression of lncRNAs in turbot are still largely unknown. In this study, a total of 12,999 lncRNAs were identified at the genome-wide level of turbot for the first time using 24 RNA-seq datasets. Sequence characteristic analyses of transcripts showed that lncRNA transcripts were shorter in average length, lower in average GC content and in average expression level as compared to the coding genes. Expression pattern analyses of lncRNAs in 12 distinct tissues showed that lncRNAs, especially lincRNA, exhibited stronger tissue-specific expression than coding genes. Moreover, 612, 1351, 1060, 875, 420 and 1689 differentially expressed (DE) lncRNAs under Vibrio anguillarum, Enteromyxum scophthalmi, and Megalocytivirus infection and heat, oxygen, and salinity stress conditions were identified, respectively. Among them, 151 and 62 lncRNAs showed differential expression under various abiotic and biotic stresses, respectively, and 11 lncRNAs differentially expressed under both abiotic and biotic stresses were selected as comprehensive stress-responsive lncRNA candidates. Furthermore, expression pattern analysis and qPCR validation both verified the comprehensive stress-responsive functions of these 11 lncRNAs. In addition, 497 significantly co-expressed target genes (correlation coefficient (R) > 0.7 and q-value < 0.05) for these 11 comprehensive stress-responsive lncRNA candidates were identified. Finally, GO and KEGG enrichment analyses indicated that these target genes were enriched mainly in molecular function, such as cytokine activity and active transmembrane transporter activity, in biological processes, such as response to stimulus and immune response, and in pathways, such as protein families: signaling and cellular processes, transporters and metabolism. These findings not only provide valuable reference resources for further research on the molecular basis and function of lncRNAs in turbot but also help to accelerate the progress of molecularly selective breeding of stress-resistant turbot strains or varieties.

Identification of a novel immune-related long noncoding RNA in carp primary macrophages associated with bisphenol A' s immunoregulatory effects

Increasing evidence suggests that long non-coding RNAs (lncRNAs) play pivotal roles in various biological processes. However, current studies on lncRNAs mostly focus on mammalian species, with little research on the functional roles of lncRNAs in teleost fish. Here, we identified a novel intergenic lncRNA (linc-93.2) in the head kidney primary macrophages of common carp (Cyprinus carpio) after exposure to a typical environmental endocrine disrupting chemical, bisphenol A (BPA). As a result, linc-93.2 was more than 3,619 bp in length and predominantly localized to the nucleus of primary macrophages other than cytoplasm, with the highest expression level in spleen followed by head kidney among different organs. Bioinformatic analysis predicted a cis-target gene, dennd1b, and 20 trans-target genes including hsp70, gna13 and rasgap, were potentially regulated by linc-93.2; NFκB and estrogen receptor (ERα) binding sites were located in the promoter region upstream of its transcription start site, which together suggested the involvement of linc-93.2 in immune and neurological functions in fish. Based on that, the expression level of linc-93.2 was determined in macrophages following acute lipopolysaccharide (LPS) and BPA treatments, both of which significantly induced linc-93.2 and IL-1β expression in cells. Moreover, a NF-κB inhibitor PDTC significantly reduced linc-93.2 expression in macrophages, but co-exposure of macrophages to PDTC with BPA or LPS could significantly rescue linc-93.2 expression, consistent with the observation on that LPS or BPA alone significantly induced both linc-93.2 and its target gene expression. Interestingly, linc-93.2 and its target gene expression was significantly suppressed by an ER antagonist ICI 182,780, however, the co-exposure of macrophages to ICI 182,780 with BPA failed to attenuate their declined expression. Overall, the current study demonstrated that linc-93.2, a novel immune-related lncRNA, may participate in the immune processes of common carp macrophages via the NF-κB and ER pathway. The results presented in this study enhance our understanding of the immunotoxin mechanisms of BPA in teleost fish.

LncRNA-adm2 targets adm2 via cid-miR-n3 and negatively regulates the inflammatory response in grass carp (Ctenopharyngodon idella)

Long non-coding RNAs (lncRNAs), which impact gene expression following pathogen infections, have garnered significant attention in recent years. Recent discoveries have revealed that lncRNAs play a crucial role in fish immune responses to pathogen infections. We investigated the influence of lncRNA-adm2 on the antibacterial immune response generated by Aeromonas hydrophila in grass carp (Ctenopharyngodon idella) through the adsorption of cid-miR-n3. Furthermore, we found that cid-miR-n3 interacts with lncRNA-adm2 and targets the 3' UTR of adm2. The upregulation of lncRNA-adm2 expression led to the suppression of pro-inflammatory cytokines (il-1β and il-6) in CIK cells, while anti-inflammatory cytokines (il-10) increased. Our research provides evidence that lncRNAs are involved in the antibacterial immune response of fish, expanding our understanding of the function of lncRNAs in teleosts.

Long non-coding RNAs as possible therapeutic targets in protozoa, and in Schistosoma and other helminths

Long non-coding RNAs (lncRNAs) emerged in the past 20 years due to massive amounts of scientific data regarding transcriptomic analyses. They have been implicated in a plethora of cellular processes in higher eukaryotes. However, little is known about lncRNA possible involvement in parasitic diseases, with most studies only detecting their presence in parasites of human medical importance. Here, we review the progress on lncRNA studies and their functions in protozoans and helminths. In addition, we show an example of knockdown of one lncRNA in Schistosoma mansoni, SmLINC156349, which led to in vitro parasite adhesion, motility, and pairing impairment, with a 20% decrease in parasite viability and 33% reduction in female oviposition. Other observed phenotypes were a decrease in the proliferation rate of both male and female worms and their gonads, and reduced female lipid and vitelline droplets that are markers for well-developed vitellaria. Impairment of female worms' vitellaria in SmLINC156349-silenced worms led to egg development deficiency. All those results demonstrate the great potential of the tools and methods to characterize lncRNAs as potential new therapeutic targets. Further, we discuss the challenges and limitations of current methods for studying lncRNAs in parasites and possible solutions to overcome them, and we highlight the future directions of this exciting field.

Integrated analysis of lncRNA and mRNA in liver of Megalobrama amblycephala post Aeromonas hydrophila infection

Background

As non-coding RNA molecules of more than 200 bp in length, long non-coding RNAs (lncRNAs) play a variety of roles in biological processes, including regulating the immune responses to bacterial infections. In recent years, there have been many in-depth studies on mammalian lncRNAs, but the relevant studies in fish are very limited. Meanwhile, since lncRNAs are not conserved among species, it is difficult to apply the existing results directly to unstudied species.

Results

To obtain the information of lncRNAs in Megalobrama amblycephala, one of the most economically important freshwater fish in China, also to better understand the biological significance of lncRNAs in the immunity system, the fish liver at 0, 4, 12, 24, and 72 h post Aeromonas hydrophila infection (hpi) were obtained for lncRNA-sequencing (lncRNA-seq). A total of 14,849 lncRNAs were identified, and 2196 lncRNAs showed significant differences at different time points post A. hydrophila infection. Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses showed that the target genes of the differentially expressed lncRNAs were enriched in several pathways related to immune such as apoptosis, inflammation, and immune response. Time-specific modules were then identified, using weighted correlation network analysis (WGCNA), and 28 modules significantly correlated with different time point after infection were found. Furthermore, four immune-related genes and six lncRNAs in the time-specific modules were subsequently verified by RT-qPCR.

Conclusions

The above findings reveal the discovery of widespread differentially expressed lncRNAs in the M. amblycephala liver post A. hydrophila infection, suggesting that lncRNAs might participate in the regulation of host response to bacterial infection, enriching the information of lncRNAs in teleost and providing a resources basis for further studies on the immune function of lncRNAs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07969-5.

Characterization and complexity of transcriptome in Gymnocypris przewalskii using single-molecule long-read sequencing and RNA-seq

The Tibetan Schizothoracinae fish Gymnocypris przewalskii has the ability to adapt to the extreme plateau environment, making it an ideal biological material for evolutionary biology research. However, the lack of well-annotated reference genomes has limited the study of the molecular genetics of G. przewalskii. To characterize its transcriptome features, we first used long-read sequencing technology in combination with RNA-seq for transcriptomic analysis. A total of 159,053 full-length (FL) transcripts were captured by Iso-Seq, having a mean length of 3,445 bp with N50 value of 4,348. Of all FL transcripts, 145,169 were well-annotated in the public database and 134,537 contained complete open reading frames. There were 4,149 pairs of alternative splicing events, of which three randomly selected were defined by RT–PCR and sequencing, and 13,293 long non-coding RNAs detected, based on all-vs.-all BLAST. A total of 118,185 perfect simple sequence repeats were identified from FL transcripts. The FL transcriptome might provide basis for further research of G. przewalskii.

Zebrafish as a Screening Model to Study the Single and Joint Effects of Antibiotics

The overuse of antibiotics combined with the limitation of wastewater facilities has resulted in drug residue accumulation in the natural environment. Thus, in recent years, the presence of antibiotic residues in the environment has raised concerns over the potential harmful effects on ecosystems and human health. The in vivo studies represent an essential step to study the potential impact induced by pharmaceutical exposure. Due to the limitations of traditional vertebrate model systems, zebrafish (Danio rerio) has recently emerged as a promising animal model to study the toxic effects of drugs and their therapeutic efficacy. The present review summarizes the recent advances made on the toxicity of seven representative classes of antibiotics, namely aminoglycosides, β-lactams, macrolides, quinolones, sulfonamides, tetracyclines and polyether antibiotics, in zebrafish, as well as the combined effects of antibiotic mixtures, to date. Despite a significant amount of the literature describing the impact of single antibiotic exposure, little information exists on the effects of antibiotic mixtures using zebrafish as an animal model. Most of the research papers on this topic have focused on antibiotic toxicity in zebrafish across different developmental stages rather than on their efficacy assessment.

Full-length transcriptome sequencing and identification of immune-related genes in the critically endangered Hucho bleekeri

Hucho bleekeri is a glacial relict and critically endangered fish restricted to the Yangtze River drainage in China. The lack of basic genomic information and immune characteristics will hinder the way toward protecting this species. In the present study, we conducted the first transcriptome analysis of H. bleekeri using the combination of SMRT and Illumina sequencing technology. Transcriptome sequencing generated a total of 93,330 non-redundant full-length unigenes with a mean length of 3072 bp. A total of 92,472 (99.08%) unigenes were annotated in at least one of the Nr protein, Swiss-Prot, KEGG, KOG, GO, Nt and Pfam databases. KEGG analysis showed that a total of 7240 unigenes belonging to 28 immune pathways were annotated to the immune system category. Meanwhile, differentially expressed genes between mucosa-associated tissues (skin, gill and hindgut) and systemic-immune tissues (spleen, head kidney and liver) were obtained. Importantly, genes participating in diverse immune signalling pathways and their expression profiles in H. bleekeri were discussed. In addition, a large number of long non-coding RNAs (lncRNAs) and simple sequence repeats (SSRs) were obtained in the H. bleekeri transcriptome. The present study will provide basic genomic information for H. bleekeri and for further research on analysing the characteristics of both the innate and adaptive immune systems of this critically endangered species.

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.

Transcriptomic Analysis of Potential "lncRNA-mRNA" Interactions in Liver of the Marine Teleost Cynoglossus semilaevis Fed Diets With Different DHA/EPA Ratios

Long non-coding RNAs (lncRNA) have emerged as important regulators of lipid metabolism and have been shown to play multifaceted roles in controlling transcriptional gene regulation, but very little relevant information has been available in fish, especially in non-model fish species. With a feeding trial on a typical marine teleost tongue sole C. semilaevis followed by transcriptomic analysis, the present study investigated the possible involvement of lncRNA in hepatic mRNA expression in response to different levels of dietary DHA and EPA, which are two most important fatty acids for marine fish. An 80-day feeding trial was conducted in a flow-through seawater system, and in this trial three experimental diets differing basically in DHA/EPA ratio, i.e., 0.61 (D/E-0.61), 1.46 (D/E-1.46), and 2.75 (D/E-2.75), were randomly assigned to 9 tanks of experimental fish. A total of 124.04 G high quality genome-wide clean data about coding and non-coding transcripts was obtained in the analysis of hepatic transcriptome. Compared to diet D/E-0.61, D/E-1.46 up-regulated expression of 178 lncRNAs and 2629 mRNAs, and down-regulated that of 47 lncRNAs and 3059 mRNAs, while D/E-2.75 resulted in much less change in gene expression. The co-expression and co-localization analysis of differentially expressed (DE) lncRNA and mRNA among dietary groups were then conducted. The co-expressed DE lncRNA and mRNA were primarily enriched in GO terms such as Metabolic process, Intracellular organelle, Catalytic activity, and Oxidoreductase activity, as well as in KEGG pathways such as Ribosome and Oxidative phosphorylation. Overlap of co-expression and co-localization analysis, i.e., lncRNA–mRNA matches “XR_523541.1–solute carrier family 16, member 5 (slc16a5)” and “LNC_000285–bromodomain adjacent to zinc finger domain 2A (baz2a),” were observed in all inter-group comparisons, indicating that they might crucially mediate the effects of dietary DHA and EPA on hepatic gene expression in tongue sole. In conclusion, this was the first time in marine teleost to investigate the possible lncRNA–mRNA interactions in response to dietary fatty acids. The results provided novel knowledge of lncRNAs in non-model marine teleost, and will serve as important resources for future studies that further investigate the roles of lncRNAs in lipid metabolism of marine teleost.

Transcriptome analysis of the effects of Cd and nanomaterial-loaded Cd on the liver in zebrafish

The wide application of engineered nanoparticles to remove heavy metals in aquatic environments has raised concerns over nanomaterial-adsorbed heavy metal toxicity. To ensure safe use of nanomaterial-heavy metal composites, understanding their biological effects at the molecular level is crucial. In the present study, we used the Illumina HiSeq technology to study the transcriptome changes induced by Cd²⁺ and nano-manganese dioxide- or nano-hydroxyapatite-adsorbed CdCl₂ composites (nMnO₂-Cd, nHAP₂₀-Cd, and nHAP₄₀-Cd) in zebrafish liver cells. We identified 545 differentially expressed genes (DEGs), 33 of which were in common between the nMnO₂-Cd, nHAP₂₀-Cd, and nHAP₄₀-Cd groups. The DEGs could be classified in four categories: hydrolases (enzymes involved in various physiological functions, including digestion, immune response, blood coagulation, and reproduction), biological binding (FMN-, actin- and metal ion-binding), metabolic enzymes (e.g., ceramidase, alpha-amylase, carboxylic ester hydrolase, and carboxypeptidase), and cell structure (cell surface, intermediate filament, and muscle myopen protein). The DEGs identified in this study are potentially useful markers to understand the physiological responses induced by Cd²⁺ and nano-Cd composites in zebrafish liver.

Genome-wide identification and differentially expression analysis of lncRNAs in tilapia

Background

Long noncoding RNAs (LncRNAs) play important roles in fundamental biological processes. However, knowledge about the genome-wide distribution and stress-related expression of lncRNAs in tilapia is still limited.

Results

Genome-wide identification of lncRNAs in the tilapia genome was carried out in this study using bioinformatics tools. 103 RNAseq datasets that generated in our laboratory or collected from NCBI database were analyzed. In total, 72,276 high-confidence lncRNAs were identified. The averaged positive correlation coefficient (r_mean = 0.286) between overlapped lncRNA and mRNA pairs showed significant differences with the values for all lncRNA-mRNA pairs (r_mean = 0.176, z statistics = − 2.45, p value = 0.00071) and mRNA-mRNA pairs (r_mean = 0.186, z statistics = − 2.23, p value = 0.0129). Weighted correlation network analysis of the lncRNA and mRNA datasets from 12 tissues identified 21 modules and many interesting mRNA genes that clustered with lncRNAs. Overrepresentation test indicated that these mRNAs enriched in many biological processes, such as meiosis (p = 0.00164), DNA replication (p = 0.00246), metabolic process (p = 0.000838) and in molecular function, e.g., helicase activity (p = 0.000102) and catalytic activity (p = 0.0000612). Differential expression (DE) analysis identified 99 stress-related lncRNA genes and 1955 tissue-specific DE lncRNA genes. MiRNA-lncRNA interaction analysis detected 72,267 lncRNAs containing motifs with sequence complementary to 458 miRNAs.

Conclusions

This study provides an invaluable resource for further studies on molecular bases of lncRNAs in tilapia genomes. Further function analysis of the lncRNAs will help to elucidate their roles in regulating stress-related adaptation in tilapia.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5115-x) contains supplementary material, which is available to authorized users.

Characterization of lncRNAs involved in cold acclimation of zebrafish ZF4 cells

Long non-coding RNAs (lncRNAs) are increasingly regarded as a key role in regulating diverse biological processes in various tissues and species. Although the cold responsive lncRNAs have been reported in plants, no data is available on screening and functional prediction of lncRNAs in cold acclimation in fish so far. Here we compared the expression profile of lncRNAs in cold acclimated zebrafish embryonic fibroblast cells (ZF4) cultured at 18°C for 30 days with that of cells cultured at 28°C as control by high-throughput sequencing. Totally 8,363 novel lncRNAs were identified. Including known and novel lncRNAs, there are 347 lncRNAs up-regulated and 342 lncRNAs down-regulated in cold acclimated cells. Among the differentially expressed lncRNAs, 74 and 61 were detected only in control cells or cold-acclimated cells, respectively. The Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) enrichment analyses of adjacent genes to the differentially expressed lncRNAs showed that the enriched genes are involved in electron transport, cell adhesion, oxidation-reduction process, and so on. We also predicted the target genes of the differentially expressed lncRNAs by looking for interactions between lncRNAs and mRNAs, and constructed an interaction network. In summary, our genome-wide systematic identification and functional prediction of cold responsive lncRNAs in zebrafish cells suggests a crucial role of lincRNAs in cold acclimation in fish.