Atypical RNAs in the coelacanth transcriptome

Comprehensive Characterization of Circular RNAs in Ovary and Testis From Nile Tilapia

Circular RNA (circRNA) is an endogenous biomolecule in eukaryotes. It has tissue- and cell-specific expression patterns and can act as a microRNA sponge or competitive endogenous RNA. Although circRNA has been found in several species in recent years, the expression profiles in fish gonad are still not fully understood. We detected the expression of circRNA in the ovary, testis, and sex-changed gonad of tilapia by high-throughput deep sequencing, and circRNA-specific computing tools. A total of 20,607 circRNAs were obtained, of which 141 were differentially expressed in the testis and ovary. Among these circRNAs, 135 circRNAs were upregulated and 6 circRNAs were downregulated in female fish. In addition, GO annotation and KEGG pathway analysis of the host genes of circRNAs indicated that these host genes were mainly involved in adherens junction, androgen production, and reproductive development, such as ZP3, PLC, delta 4a, ARHGEF10, and HSD17b3. It is worth noting that we found that circRNAs in tilapia gonads have abundant miRNA-binding sites. Among them, 935 circRNAs have a regulatory effect on miR-212, 856 circRNAs have a regulatory effect on miR-200b-3p, and 529 circRNAs have a regulatory effect on miR-200b-5p. Thus, our findings provide a new evidence for circRNA–miRNA networks in the gonads in tilapia.

Emerging roles of circRNAs in regulating thermal and hypoxic stresses in Apostichopus japonicus (Echinodermata: Holothuroidea)

Some sea cucumbers are economically and ecologically important, but they are threatened by thermal and hypoxic stress in changing oceanographic conditions. We construct circRNAs profiles, reveal circRNAs characters, and illustrate the potential regulatory roles of circRNAs in one commercially important species of sea cucumber, Apostichopus japonicus. Reads are distributed in intergenic (44.14%), exonic (48.26%) and intronic (7.60%) regions of the genome. A total of 1684 circRNAs were identified, and the most common spliced length is 269 nt in the present study. In three treatments (HT [thermal stress], LO [hypoxic stress], and HL [combined thermal and hypoxic stress]), 24, 27 and 27 differentially expressed (DE) circRNAs were identified, respectively. Five novel DE-circRNAs commonly occur in these treatments (novel_circ_0003311, novel_circ_0000229, novel_circ_0003944, novel_circ_0001458 and novel_circ_0000707), and based on them, potential circRNA-miRNA binding pairs were predicted. Sanger sequencing, RNase R treatment experiment and qPCR validation identified the accuracy of the circRNAs. Key circRNAs identified in the present study were covalently closed and were more stable under RNase R treatment than linear RNAs. Based on function analysis, circRNAs could regulate metabolic process, signal transduction, and ion responses in A. japonicus when exposed to thermal and hypoxic stress, and 'regulation of response to stimulus' is a common gene ontology (GO) term that is significantly enriched in each treatment; GO terms for 'DNA' and 'stress' are commonly enriched in heat-related treatments (HT and HL); and GO terms for 'protein' are commonly enriched in hypoxia-related treatments (LO and HL). When environmentally stressed, 'metabolism,' 'transport and catabolism,' 'membrane transport,' and 'signal transduction' were significantly responded in sea cucumber based on KEGG analysis. We provide insights into circRNA functions in stress regulation and lay a foundation for invertebrate circRNA research.

Identification and Characterization of Novel circRNAs Involved in Muscle Growth of Blunt Snout Bream (Megalobrama amblycephala)

Circular RNAs (circRNAs), a novel class of endogenous RNAs, have been recognized to play important roles in the growth of animals. However, the regulatory mechanism of circRNAs on fish muscle growth is still unclear. In this study, we performed whole transcriptome analysis of skeletal muscles from two populations with different growth rates (fast-growing and slow-growing) of blunt snout bream (Megalobrama amblycephala), an important fish species for aquaculture. The selected circRNAs were validated by qPCR and Sanger sequencing. Pairs of circRNA–miRNA–mRNA networks were constructed with the predicted differentially expressed (DE) pairs, which revealed regulatory roles in muscle myogenesis and hypertrophy. As a result, a total of 445 circRNAs were identified, including 42 DE circRNAs between fast-growing (FG) and slow-growing (SG) groups. Many of these DE circRNAs were related with aminoglycan biosynthetic and metabolic processes, cytokinetic processes, and the adherens junction pathway. The functional prediction results showed that novel_circ_0001608 and novel_circ_0002886, competing to bind with dre-miR-153b-5p and dre-miR-124-6-5p, might act as competing endogenous RNAs (ceRNAs) to control MamblycephalaGene14755 (pik3r1) and MamblycephalaGene10444 (apip) level, respectively, thus playing an important regulatory role in muscle growth. Overall, these results will not only help us to further understand the novel RNA transcripts in M. amblycephala, but also provide new clues to investigate the potential mechanism of circRNAs regulating fish growth and muscle development.

seekCRIT: Detecting and characterizing differentially expressed circular RNAs using high-throughput sequencing data

Over the past two decades, researchers have discovered a special form of alternative splicing that produces a circular form of RNA. Although these circular RNAs (circRNAs) have garnered considerable attention in the scientific community for their biogenesis and functions, the focus of current studies has been on the tissue-specific circRNAs that exist only in one tissue but not in other tissues or on the disease-specific circRNAs that exist in certain disease conditions, such as cancer, but not under normal conditions. This approach was conducted in the relative absence of methods that analyze a group of common circRNAs that exist in both conditions, but are more abundant in one condition relative to another (differentially expressed). Studies of differentially expressed circRNAs (DECs) between two conditions would serve as a significant first step in filling this void. Here, we introduce a novel computational tool, seekCRIT (seek for differentially expressed CircRNAs In Transcriptome), that identifies the DECs between two conditions from high-throughput sequencing data. Using rat retina RNA-seq data from ischemic and normal conditions, we show that over 74% of identifiable circRNAs are expressed in both conditions and over 40 circRNAs are differentially expressed between two conditions. We also obtain a high qPCR validation rate of 90% for DECs with a FDR of < 5%. Our results demonstrate that seekCRIT is a novel and efficient approach to detect DECs using rRNA depleted RNA-seq data. seekCRIT is freely downloadable at https://github.com/UofLBioinformatics/seekCRIT. The source code is licensed under the MIT License. seekCRIT is developed and tested on Linux CentOS-7.

The focus of circRNA studies has been on condition-specific circRNAs, however, there are situations in which circRNAs exist in both conditions with different abundance. Here, we introduce a new and robust analytic software, seekCRIT (seek for differentially expressed CircRNAs In Transcriptome), that identifies the differentially expressed circRNAs (DECs) between two conditions from high-throughput sequencing data. seekCRIT provides a straightforward normalized quantification of circRNAs and statistical measures by adapting a junction-count-based estimation approach. Using publicly available ribosomal RNA depleted RNA-seq data and our own rat retina RNA-seq data, we show that seekCRIT can efficiently detect circRNAs and identify DECs. We also obtain a high qPCR validation rate of 90% for DECs with a FDR of < 5%. Our results demonstrate that seekCRIT is a novel and efficient software to detect DECs using rRNA depleted RNA-seq data.

Identification of circular RNAs in the ovarian follicles of Meishan and Duroc sows during the follicular phase

Circular RNAs (circRNAs) are a newly discovered class of endogenous non-coding RNAs that play an important role in growth and development by regulating gene expression and participating in a variety of biological processes. However, the role of circRNAs in porcine follicles remains unclear. Therefore, this study examined middle-sized ovarian follicles obtained from Meishan and Duroc sows at day 4 of the follicular phase. High-throughput RNA sequencing (RNA-seq) was utilized to construct circRNAs, and differential expression was identified. The findings were validated using reverse transcription PCR (RT-PCR) and DNA sequencing, GO and KEGG analyses were performed, and potential miRNA targets were identified. The RNA-seq identified a total of 15,866 circRNAs, with 244 differentially expressed in the Meishan relative to the Duroc (111 up-regulated and 133 down-regulated). The RT-PCR finding confirmed the RNA-seq results, and quantitative real-time PCR (qPCR) analysis examining a subset of the circRNAs showed that they are resistant to RNase R digestion. Bioinformatics analysis (GO and KEGG) showed that the host genes associated with the differentially expressed circRNAs are involved in reproduction and follicular development signaling pathways. Furthermore, many of the circRNAs were found to interact with miRNAs that are associated with follicular development. This study presents a new perspective for studying circRNAs and provides a valuable resource for further examination into the potential roles of circRNAs in porcine follicular development.

Single molecule real-time sequencing and RNA-seq unravel the role of long non-coding and circular RNA in the regulatory network during Nile tilapia (Oreochromis niloticus) infection with Streptococcus agalactiae

BACKGROUND

The tilapia aquaculture industry is facing heavy economic losses due to Streptococcus agalactiae (S. agalactiae) infections. While progress has been made in past years, the lack of a high-quality tilapia genome and transcript annotations makes systematic and comprehensive exploration for a non-coding RNA regulatory network associated with the infection process unfeasible, and it stunts further research focused on disease defense and treatment. Herein, single molecular real time sequencing (SMRT-Seq) and RNA-seq data were utilized to generate a high-quality transcript annotation. In addition, Changes in mRNA and non-coding RNA expression were also analyzed during a S. agalactiae infection in tilapia.

FINDINGS

In total, 16.79 Gb of clean data were obtained by sequencing on six SMRT cells, with 712,294 inserts (326,645 full-length non-chimeric reads and 354,188 non-full-length reads). A total of 197,952 consensus transcripts were obtained. Additionally, 55,857 transcript sequences were acquired, with 12,297 previously annotated and 43,560 newly identified transcripts. To further examine the immune response in Oreochromis niloticus following a S. agalactiae infection, a total of 470.62 Gb of clean data was generated by sequencing a library containing 18 S. agalactiae infected tilapia samples. Of the identified genes, 9911 were newly exploited, of which 7102 were functional annotated. Furthermore, 7874 mRNAs, 1281 long non-coding RNAs (out of 21,860 long non-coding RNAs), and 61 circular RNAs (out of 1026 circular RNAs) were found to be differentially expressed during infection, with the 1026 circRNAs not previously identified in tilapia. Moreover, k-means clustering and WGCNA analyses revealed that the immune response of tilapia to a S. agalactiae infection can be divided into three stages: cytokines driven rapid immune response, energy metabolism promotion, and the production of lysosomes and phagosomes. During this response, the head kidney and spleen have synergistic effects, while maintaining independent characteristics. Finally, lncRNA-mRNA (trans and cis), lncRNA-miRNA-mRNA, and circRNA-miRNA-mRNA regulatory networks were constructed and revealed that non-coding RNA is involved in the regulation of immune-related genes.

CONCLUSIONS

This study generated a greatly-improved transcript annotation for tilapia using long-read PacBio sequencing technology, and revealed the presence of a regulatory network comprised of non-coding RNAs in Nile tilapia infected with S. agalactiae.

Functions of Circular RNAs Involved in Animal Skeletal Muscle Development – A Review

Circular RNAs (circRNAs) have been identified in the skeletal muscle of numerous species of animals. Their abundance, diversity, and their dynamic expression patterns have been revealed in various developmental stages and physiological conditions in skeletal muscles. Recently, studies have made known that circRNAs widely participate in muscle cell proliferation and differentiation. They are also involved in other life processes such as functioning as microRNA (miRNA) sponges, regulators of splicing and transcription, and modifiers of parental gene expression with emerging pieces of evidence indicating a high chance of playing a vital role in several cells and tissues, especially the muscles. Other research has emphatically stated that the growth and development of skeletal muscle are regulated by proteins as well as non-coding RNAs, which involve circRNAs. Therefore, circRNAs have been considered significant biological regulators for understanding the molecular mechanisms of myoblasts. Here, we discuss how circRNAs are abundantly expressed in muscle (myoblast) and their critical roles in growth and development.

Identification of Potential Immune-Related circRNA-miRNA-mRNA Regulatory Network in Intestine of Paralichthys olivaceus During Edwardsiella tarda Infection

Olive flounder (Paralichthys olivaceus) is an important economical flatfish in Japan, Korea, and China, but its production has been greatly threatened by disease outbreaks. In this research, we aimed to explore the immune responsive mechanism of P. olivaceus against Edwardsiella tarda infection by profiling the expression of circRNA, miRNA, and mRNA by RNA-seq and constructing a regulatory circular circRNA–miRNA–mRNA network. Illumina sequencing of samples from normal control (H0), 2 h (H2), 8 h (H8), and 12 h (H12) post-challenge was conducted. Differentially expressed (DE) circRNA (DE–circRNAs), miRNAs (DE–miRNAs), and mRNAs [differential expression genes (DEGs)] between challenge and control groups were identified, resulting in a total of 62 DE–circRNAs, 39 DE–miRNAs, and 3,011 DEGs. Based on the differentially expressed gene results, miRNA target interactions (circRNA–miRNA pairs and miRNA–mRNA pairs) were predicted by MiRanda software. Once these paired were combined, a preliminary circRNA–miRNA–mRNA network was generated with 198 circRNA–miRNA edges and 3,873 miRNA–mRNA edges, including 44 DE–circRNAs, 32 DE–miRNAs, and 1,774 DEGs. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis was performed to evaluate the function of the DEGs in this network, and we focused and identified two important intestinal immune pathways (herpes simplex infection and intestinal immune network for IgA production) that showed statistical significance between the challenge and control groups. Furthermore, three critical DEGs (nectin2, MHC II α-chain, and MHC II β-chain) were identified, mapped into the preliminary circRNA–miRNA–mRNA network, and new circRNA–miRNA–mRNA regulatory networks were constructed. In conclusion, we, for the first time, identified circRNA–miRNA–mRNA network from P. olivaceus in the pathogenesis of E. tarda and provided valuable resources for further analyses of the molecular mechanisms and signaling networks.

Genome-Wide Identification of Circular RNAs Revealed the Dominant Intergenic Region Circularization Model in Apostichopus japonicus

Circular RNAs (circRNAs) were recently recognized to act as competing endogenous RNAs and play roles in gene expression regulation. Previous studies in humans and silkworms have shown that circRNAs take part in immune regulation. Here, we conducted coelomocyte circRNA sequencing to explore its immune functions in healthy and skin ulceration syndrome (SUS)-diseased sea cucumbers. A total of 3,592 circRNAs were identified in libraries with diversified circularization patterns compared with animal models. The common intron-pairing-driven circularization models are not popular in sea cucumber genome, which was replaced with intergenic region circularization. The accuracy of these identified circRNAs was further validated by Sanger sequencing and RNase R-treated assays. Expression profile analysis indicated that 117 circRNAs were upregulated and 144 circRNAs were downregulated in SUS-diseased condition, of which 71.6% were intergenic-type circRNAs. The interaction network of differentially expressed circRNAs and microRNAs (miRNAs) was constructed and showed that miR-2008 and miR-31, detected with significantly differential expression in SUS-affected samples in a previous study, were predicted to be regulated by 10 and 11 differentially expressed circRNAs with more than 10 binding sites, respectively. Moreover, seven circRNAs were further validated by quantitative real-time PCR, whose variation trends were consistent with circRNA sequencing. All our results supported that intergenic-type circRNAs might have a dominant function in Apostichopus japonicas immune response by acting as miRNA regulators.

Genome-Wide circRNA Profiling from RNA-seq Data

The genome-wide expression patterns of circular RNAs (circRNAs) are of increasing interest for their potential roles in normal cellular homeostasis, development, and disease. Thousands of circRNAs have been annotated from various species in recent years. Analysis of publically available or user-generated rRNA-depleted total RNA-seq data can be performed to uncover new circRNA expression trends. Here we provide a primer for profiling circRNAs from RNA-seq datasets. The description is tailored for the wet lab scientist with limited or no experience in analyzing RNA-seq data. We begin by describing how to access and interpret circRNA annotations. Next, we cover converting circRNA annotations into junction sequences that are used as scaffolds to align RNA-seq reads. Lastly, we visit quantifying circRNA expression trends from the alignment data.

From Glacier to Sauna: RNA-Seq of the Human Pathogen Black Fungus Exophiala dermatitidis under Varying Temperature Conditions Exhibits Common and Novel Fungal Response

Exophiala dermatitidis (Wangiella dermatitidis) belongs to the group of the so-called black yeasts. Thanks in part to its thick and strongly melanized cell walls, E. dermatitidis is extremely tolerant to various kinds of stress, including extreme pH, temperature and desiccation. E. dermatitidis is also the agent responsible for various severe illnesses in humans, such as pneumonia and keratitis, and might lead to fatal brain infections. Due to its association with the human environment, its poly-extremophilic lifestyle and its pathogenicity in humans, E. dermatitidis has become an important model organism. In this study we present the functional analysis of the transcriptional response of the fungus at 1°C and 45°C, in comparison with that at 37°C, for two different exposition times, i.e. 1 hour and 1 week. At 1°C, E. dermatitidis uses a large repertoire of tools to acclimatize, such as lipid membrane fluidization, trehalose production or cytoskeleton rearrangement, which allows the fungus to remain metabolically active. At 45°C, the fungus drifts into a replicative state and increases the activity of the Golgi apparatus. As a novel finding, our study provides evidence that, apart from the protein coding genes, non-coding RNAs, circular RNAs as well as fusion-transcripts are differentially regulated and that the function of the fusion-transcripts can be related to the corresponding temperature condition. This work establishes that E. dermatitidis adapts to its environment by modulating coding and non-coding gene transcription levels and through the regulation of chimeric and circular RNAs.

Purifying selection on splice-related motifs, not expression level nor RNA folding, explains nearly all constraint on human lincRNAs

There are two strong and equally important predictors of rates of human protein evolution: The amount the gene is expressed and the proportion of exonic sequence devoted to control splicing, mediated largely by selection on exonic splice enhancer (ESE) motifs. Is the same true for noncoding RNAs, known to be under very weak purifying selection? Prior evidence suggests that selection at splice sites in long intergenic noncoding RNAs (lincRNAs) is important. We now report multiple lines of evidence indicating that the great majority of purifying selection operating on lincRNAs in humans is splice related. Splice-related parameters explain much of the between-gene variation in evolutionary rate in humans. Expression rate is not a relevant predictor, although expression breadth is weakly so. In contrast to protein-coding RNAs, we observe no relationship between evolutionary rate and lincRNA stability. As in protein-coding genes, ESEs are especially abundant near splice junctions and evolve slower than non-ESE sequence equidistant from boundaries. Nearly all constraint in lincRNAs is at exon ends (N.B. the same is not witnessed in Drosophila). Although we cannot definitely answer the question as to why splice-related selection is so important, we find no evidence that splicing might enable the nonsense-mediated decay pathway to capture transcripts incorrectly processed by ribosomes. We find evidence consistent with the notion that splicing modifies the underlying chromatin through recruitment of splice-coupled chromatin modifiers, such as CHD1, which in turn might modulate neighbor gene activity. We conclude that most selection on human lincRNAs is splice mediated and suggest that the possibility of splice-chromatin coupling is worthy of further scrutiny.

circBase: a database for circular RNAs

Recently, several laboratories have reported thousands of circular RNAs (circRNAs) in animals. Numerous circRNAs are highly stable and have specific spatiotemporal expression patterns. Even though a function for circRNAs is unknown, these features make circRNAs an interesting class of RNAs as possible biomarkers and for further research. We developed a database and website, "circBase," where merged and unified data sets of circRNAs and the evidence supporting their expression can be accessed, downloaded, and browsed within the genomic context. circBase also provides scripts to identify known and novel circRNAs in sequencing data. The database is freely accessible through the web server at http://www.circbase.org/.