Function of Circular RNAs in Fish and Their Potential Application as Biomarkers

Pathogen-induced apoptosis in echinoderms: A review

Echinoderms possess unique biological traits that make them valuable models in immunology, regeneration, and developmental biology studies. As a class rich in active substances with significant nutritional and medicinal value, echinoderms face threats from marine pathogens, including bacteria, viruses, fungi, protozoa, and parasites, which have caused substantial economic losses in echinoderm aquaculture. Echinoderms counteract pathogen invasion through innate immunity and programmed cell death, in particular, with apoptosis being essential for eliminating infected or damaged cells and maintaining homeostasis in many echinoderm cell types. Despite the importance of this process, there is a lack of comprehensive and updated reviews on this topic. This review underscores that echinoderm apoptotic pathways exhibit a complexity comparable to that of vertebrates, featuring proteins with unique domains that may indicate the presence of novel signaling mechanisms. We synthesize current knowledge on how echinoderms utilize diverse transcriptional and post-transcriptional mechanisms to regulate apoptosis in response to pathogen infections and explore how pathogens have evolved strategies to manipulate echinoderm apoptosis, either by inhibiting it to create survival niches or by inducing excessive apoptosis to weaken the host. By elucidating the primary apoptotic pathways in echinoderms and the host-pathogen interactions that modulate these pathways, this review aims to reveal new mechanisms of apoptosis in animal immune defense and provide insights into the evolutionary arms race between hosts and pathogens.

Endemic Radiation of African Moonfish, Selene dorsalis (Gill 1863), in the Eastern Atlantic: Mitogenomic Characterization and Phylogenetic Implications of Carangids (Teleostei: Carangiformes)

This study offers an in-depth analysis of the mitochondrial genome of Selene dorsalis (Gill 1863), a species native to the Eastern Atlantic Ocean. The circular mitochondrial DNA molecule measures 16,541 base pairs and comprises 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA genes, and a control region (CR). The nucleotide composition exhibits a notable adenine-thymine (AT) bias, accounting for 53.13%, which aligns with other species in the Carangidae family. Most PCGs initiate with the ATG codon, with the exception of Cytochrome C oxidase subunit I, which starts with GTG. Analysis of relative synonymous codon usage reveals that leucine and serine are the most prevalent amino acids in the mitochondrial genome of S. dorsalis and its congeners (S. vomer and S. setapinnis). All tRNAs display the typical cloverleaf structure, though tRNA Serine (S1) lacks a dihydrouracil arm. Pairwise comparisons of synonymous and nonsynonymous substitutions for all PCGs yielded values below '1', indicating strong purifying selection. The CR spans 847 bp, representing 5.12% of the mitochondrial genome, and is characterized by high AT content (62.81%). It is situated between tRNA-Pro (TGG) and tRNA-Phe (GAA). The CR contains conserved sequence blocks, with CSB-1 being the longest at 22 bp and CSB-D the shortest at 18 bp. Phylogenetic analysis, using Bayesian and Maximum-likelihood trees constructed from concatenated PCGs across 72 species, successfully differentiates S. dorsalis from other carangids. This study also explores how ocean currents and gyres might influence lineage diversification and parapatric speciation of Selene species between the Atlantic and Pacific Oceans. These results highlight the importance of the mitochondrial genome in elucidating the structural organization and evolutionary dynamics of S. dorsalis and its relatives within marine ecosystems.

Embryonic temperature has long-term effects on muscle circRNA expression and somatic growth in Nile tilapia

Embryonic temperature has a lasting impact on muscle phenotype in vertebrates, involving complex molecular mechanisms that encompass both protein-coding and non-coding genes. Circular RNAs (circRNAs) are a class of regulatory RNAs that play important roles in various biological processes, but the effect of variable thermal conditions on the circRNA transcriptome and its long-term impact on muscle growth plasticity remains largely unexplored. To fill this knowledge gap, we performed a transcriptomic analysis of circRNAs in fast muscle of Nile tilapia (Oreochromis niloticus) subjected to different embryonic temperatures (24°C, 28°C and 32°C) and then reared at a common temperature (28°C) for 4 months. Nile tilapia embryos exhibited faster development and subsequently higher long-term growth at 32°C compared to those reared at 28°C and 24°C. Next-generation sequencing data revealed a total of 5,141 unique circRNAs across all temperature groups, of which 1,604, 1,531, and 1,169 circRNAs were exclusively found in the 24°C, 28°C and 32°C groups, respectively. Among them, circNexn exhibited a 1.7-fold (log2) upregulation in the 24°C group and a 1.3-fold (log2) upregulation in the 32°C group when compared to the 28°C group. Conversely, circTTN and circTTN_b were downregulated in the 24°C groups compared to their 28°C and 32°C counterparts. Furthermore, these differentially expressed circRNAs were found to have multiple interactions with myomiRs, highlighting their potential as promising candidates for further investigation in the context of muscle growth plasticity. Taken together, our findings provide new insights into the molecular mechanisms that may underlie muscle growth plasticity in response to thermal variation in fish, with important implications in the context of climate change, fisheries and aquaculture.

Comprehensive analysis of circRNA-miRNA-mRNA networks in the kidney of snakehead (Channa argus) response to Nocardia seriolae challenge

Non-coding RNAs (ncRNAs) play vital roles in regulating the expression levels of genes that control essential biological functions, including immune response to bacterial infections in teleost. To dissect the roles of ncRNAs in the Channa argus (snakehead), a systematic analysis of the expression profiles of circRNA, miRNA and mRNA, as well as competing endogenous RNAs (ceRNA) regulatory networks in the kidney of snakehead following Nocardia seriolae infection were performed in the present study. A total of 111 differentially expressed circRNAs, 706 differentially expressed miRNAs, and 2548 differentially expressed mRNAs were identified in the N. seriolae infected snakehead. Based on these differently expressed RNAs, we identified 55 circRNA-mRNA pairs, 124 miRNA-mRNA pairs, and 35 circRNA-miRNA-mRNA regulatory networks, including dre-miR-103-CD302, dre-miR-27e-IGSF3, novel_circ_0005462/novel_403-IGKC, novel_circ_0001750/novel_circ_0002162-novel_477-OCLN, and novel_circ_0003847-novel_4-KCNAB3. In addition, luciferase reporter assay was employed to detect the target relationships of several circRNA-miRNA-mRNA pairs. Taken together, this study demonstrates that the genes associated with immunity and structures in the kidney of snakehead can be regulated by circRNAs and miRNAs at post-transcription levels, and provided theoretical guidance for ncRNAs studies for other teleost. However, further studies are still in great need to validate the regulatory mechanisms of ncRNAs in snakehead.

circWSB1 promotes tumor progression in ccRCC via circWSB1/miR-182-5p/WSB1 axis

Clear cell renal cell carcinoma (ccRCC) is one of the most prevalent urological carcinomas with a low overall 5-year survival rate, and its prognosis remains dismal. circular RNAs (circRNAs) has been discovered to be important regulators in ccRCC. However, the specific regulatory mechanisms of circRNAs and their impact on phenotypes require further in-depth research. circRNA microarray sequencing analysis was used in this study to explore the expression pattern of circRNAs in ccRCC. circWSB1 was discovered, and we evaluated its derivation, potential diagnostic efficacy, and prognostic significance in ccRCC tissues. We discovered that circWSB1 is highly expressed in ccRCC. We identified that circWSB1 interacts with miR-182-5p and upregulates the expression of its host gene, WSB1. Through models in vivo and in vitro models, we found that circWSB1 increases WSB1 expression via the circWSB1/miR-182-5p/WSB1 axis, which promotes ccRCC cell proliferation and migration. The high expression of circWSB1 and WSB1 is correlated with poorer clinical prognosis and pathological grading. circWSB1 diminishes the inhibitory impact of miR-182-5p on WSB1 and increases WSB1 expression, thereafter promoting ccRCC development. Our findings provide a promising predictive biomarker and therapeutic target for ccRCC.

Changes in the gut bacterial community affect miRNA profiles in Riptortus pedestris under different rearing conditions

Insects possess complex and dynamic gut microbial system, which contributes to host nutrient absorption, reproduction, energy metabolism, and protection against stress. However, there are limited data on interactions of host-gut bacterial microbiota through miRNA (microRNA) regulation in a significant pest, Riptortus pedestris. Here, we performed the 16S rRNA amplicon sequencing and small RNA sequencing from the R. pedestris gut under three environmental conditions and antibiotic treatment, suggesting that we obtained a large amount of reads by assembly, filtration and quality control. The 16S rRNA amplicon sequencing results showed that the abundance and diversity of gut bacterial microbiota were significantly changed between antibiotic treatment and other groups, and they are involved in metabolism and biosynthesis-related function based on functional prediction. Furthermore, we identified different numbers of differentially expressed unigenes (DEGs) and differentially expressed miRNAs (DEMs) based on high-quality mappable reads, which were enriched in various immune-related pathways, including Toll-like receptor, RIG-I-like receptor, NOD-like receptor, JAK/STAT, PI3K/Akt, NF-κB, MAPK signaling pathways, and so forth, using GO and KEGG enrichment analysis. Later on, the identified miRNAs and their target genes in the R. pedestris gut were predicted and randomly selected to construct an interaction network. Finally, our study indicated that alterations in the gut bacterial microbiota are significantly positively or negatively associated with DEMs of the Toll/Imd signaling pathway with Pearson correlation analysis. Taken together, the results of our study lay the foundation for further deeply understanding the interactions between the gut microbiota and immune responses in R. pedestris through miRNA regulation, and provide the new basis for pest management in hemipteran pests.

The emerging role of circRNAs on skeletal muscle development in economical animals

CircRNAs are a novel type of closed circular molecules formed through a covalent bond lacking a 5'cap and 3' end tail, which mainly arise from mRNA precursor. They are widely distributed in plants and animals and are characterized by stable structure, high conservativeness in cells or tissues, and showed the expression specificity at different stages of development in different tissues. CircRNAs have been gradually attracted wide attention with the development of RNA sequencing, which become a new research hotspot in the field of RNA. CircRNAs play an important role in gene expression regulation. Presently, the related circRNAs research in the regulation of animal muscle development is still at the initial stage. In this review, the formation, properties, biological functions of circRNAs were summarized. The recent research progresses of circRNAs in skeletal muscle growth and development from economic animals including livestock, poultry and fishes were introduced. Finally, we proposed a prospective for further studies of circRNAs in muscle development, and we hope our research could provide new ideas, some theoretical supports and helps for new molecular genetic markers exploitation and animal genetic breeding in future.

A circRNA therapy based on Rnf103 to inhibit Vibrio anguillarum infection

The losses caused by Vibrio infections in the aquaculture industry are challenging to quantify. In the face of antibiotic resistance, a natural and environmentally friendly alternative is urgently needed. In this study, we identify E3 ubiquitin-protein ligase RNF103 (rnf103) as a crucial target involved in immune evasion by Vibrio anguillarum. Our research demonstrates that Rnf103 promotes immune escape by inhibiting Traf6. Interestingly, we discover a circular RNA (circRNA), circRnf103, formed by reverse splicing of the Rnf103 gene. Predictive analysis and experimentation reveal that circRnf103 encodes Rnf103-177aa, a protein that competes with Rnf103 and binds to Traf6, preventing its degradation. Notably, circRnf103 therapy induces Rnf103-177aa protein production in zebrafish. In zebrafish models, circRnf103 exhibits significant effectiveness in treating V. anguillarum infections, reducing organ burden. These findings highlight the potential of circRNA therapy as a natural and innovative approach to combat infectious diseases sustainably, particularly in aquaculture and environmental management.

Circular RNAs are associated with the resistance to Newcastle disease virus infection in duck cells

Introduction

Newcastle disease virus (NDV) is prevalent worldwide with an extensive host range. Among birds infected with velogenic NDV strains, chickens experience high pathogenicity and mortality, whereas ducks mostly experience mild symptoms or are asymptomatic. Ducks have a unique, innate immune system hypothesized to induce antiviral responses. Circular RNAs (circRNAs) are among the most abundant and conserved eukaryotic transcripts. These participate in innate immunity and host antiviral response progression.

Methods

In this study, circRNA expression profile differences post-NDV infection in duck embryo fibroblast (DEF) cells were analyzed using circRNA transcriptome sequencing. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were used to reveal significant enrichment of differentially expressed (DE) circRNAs. The circRNA-miRNA-mRNA interaction networks were used to predict the related functions of circRNAs. Moreover, circ-FBXW7 was selected to determine its effect on NDV infection in DEFs.

Results

NDV infection altered circRNA expression profiles in DEF cells, and 57 significantly differentially expressed circRNAs were identified post-NDV infection. DEF responded to NDV by forming circRNAs to regulate apoptosis-, cell growth-, and protein degradation-related pathways via GO and KEGG enrichment analyses. circRNA-miRNA-mRNA interaction networks demonstrated that DEF cells combat NDV infection by regulating cellular pathways or apoptosis through circRNA-targeted mRNAs and miRNAs. circ-FBXW7 overexpression and knockdown inhibited and promoted viral replication, respectively. DEF cells mainly regulated cell cycle alterations or altered cellular sensing to combat NDV infection.

Conclusion

These results demonstrate that DEF cells exert antiviral responses by forming circRNAs, providing novel insights into waterfowl antiviral responses.

CircPrime: a web-based platform for design of specific circular RNA primers

BACKGROUND

Circular RNAs (circRNAs) are covalently closed-loop RNAs with critical regulatory roles in cells. Tens of thousands of circRNAs have been unveiled due to the recent advances in high throughput RNA sequencing technologies and bioinformatic tools development. At the same time, polymerase chain reaction (PCR) cross-validation for circRNAs predicted by bioinformatic tools remains an essential part of any circRNA study before publication.

RESULTS

Here, we present the CircPrime web-based platform, providing a user-friendly solution for DNA primer design and thermocycling conditions for circRNA identification with routine PCR methods.

CONCLUSIONS

User-friendly CircPrime web platform ( http://circprime.elgene.net/ ) works with outputs of the most popular bioinformatic predictors of circRNAs to design specific circular RNA primers. CircPrime works with circRNA coordinates and any reference genome from the National Center for Biotechnology Information database).

The novel circular RNA CircMef2c is positively associated with muscle growth in Nile tilapia

Muscle growth in teleosts is a complex biological process orchestrated by numerous protein-coding genes and non-coding RNAs. A few recent studies suggest that circRNAs are involved in teleost myogenesis, but the molecular networks involved remain poorly understood. In this study, an integrative omics approach was used to determine myogenic circRNAs in Nile tilapia by quantifying and comparing the expression profile of mRNAs, miRNAs, and circRNAs in fast muscle from full-sib fish with distinct growth rates. There were 1947 mRNAs, 9 miRNAs, and 4 circRNAs differentially expressed between fast- and slow-growing individuals. These miRNAs can regulate myogenic genes and have binding sites for the novel circRNA circMef2c. Our data indicate that circMef2c may interact with three miRNAs and 65 differentially expressed mRNAs to form multiple competing endogenous RNA networks that regulate growth, thus providing novel insights into the role of circRNAs in the regulation of muscle growth in teleosts.

Comprehensive RNA-Seq Analysis Pipeline for Non-Model Organisms and Its Application in Schmidtea mediterranea

RNA sequencing (RNA-seq) is a high-throughput technology that provides in-depth information on transcriptome. The advancement and dropping costs of RNA sequencing, accompanied by more available reference genomes for different species, make transcriptome analysis in non-model organisms possible. Current obstacles in analyzing RNA-seq data include a lack of functional annotation, which may complicate the process of linking genes to corresponding functions. Here, we provide a one-stop RNA-seq analysis pipeline, PipeOne-NM, for transcriptome functional annotation, non-coding RNA identification, and transcripts alternative splicing analysis of non-model organisms, intended for use with Illumina platform-based RNA-seq data. We performed PipeOne-NM on 237 Schmidtea mediterranea RNA-seq runs and assembled a transcriptome with 84,827 sequences from 49,320 genes, identifying 64,582 mRNA from 35,485 genes, 20,217 lncRNA from 17,084 genes, and 3481 circRNAs from 1103 genes. In addition, we performed a co-expression analysis of lncRNA and mRNA and identified that 1319 lncRNA co-express with at least one mRNA. Further analysis of samples from S. mediterranea sexual and asexual strains revealed the role of sexual reproduction in gene expression profiles. Samples from different parts of asexual S. mediterranea revealed that differential expression profiles of different body parts correlated with the function of conduction of nerve impulses. In conclusion, PipeOne-NM has the potential to provide comprehensive transcriptome information for non-model organisms on a single platform.

Integrative analysis of miRNA and mRNA expression associated with the immune response in the intestine of rainbow trout (Oncorhynchus mykiss) infected with infectious hematopoietic necrosis virus

Rainbow trout (Oncorhynchus mykiss), an economically important cold-water fish cultured worldwide, suffers from infectious hematopoietic necrosis virus (IHNV) infection, resulting in huge financial losses. In order to understand the immune response of rainbow trout during virus infection, we explored trout intestine transcriptome profiles following IHNV challenge, and identified 3355 differentially expressed genes (DEGs) and 80 differentially expressed miRNAs (DEMs). Transcriptome analysis revealed numerous DEGs involved in immune responses, such as toll-like receptor 3 (TLR3), toll-like receptor 7/8 (TLR7/8), tripartite motif-containing 25 (TRIM25), DExH-Box helicase 58 (DHX58), interferon-induced with helicase C domain 1 (IFIH1), interferon regulatory factor 3 (IRF3/7), signal transducer and activator of transcription 1 (STAT1) and heat shock protein 90-alpha 1 (HSP90A1). Integrated analysis identified five key miRNAs (miR-19-y, miR-181-z, miR-203-y, miR-143-z and miR-206-y) targeting at least two important immune genes (TRIM25, DHX58, STAT1, TLR7/8 and HSP90A1). Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses showed that DEGs and target genes were significantly enriched in various immune-related terms including immune system process, binding, cell part and pathways of Toll-like receptor signalling, RIG-I-like receptor signalling, NOD-like receptor signalling, JAK-STAT signalling, PI3K-Akt signalling, NF-kappa B signalling, IL-17 signalling and AGE-RAGE signalling. In addition, protein-protein interaction networks (PPI) was used to display highly interactive DEG networks involving eight immune-related pathways. The expression trends of 12 DEGs and 10 DEMs were further verified by quantitative real-time PCR, which confirmed the reliability of the transcriptome sequencing results. This study expands our understanding of the immune response of rainbow trout infected with IHNV, and provides valuable resources for future studies on the immune molecular mechanism and disease resistance breeding.

Functions and clinical significance of circular RNAs in acute myeloid leukemia

Circular RNAs (circRNAs) are a class of covalently closed single-stranded RNA molecules. Four types of circRNAs have been reported in animal cells, and they have typical characteristics in their biogenesis, nuclear export and degradation. Advances in our understanding of the molecular functions of circRNAs in sponging microRNAs, modulating transcription, regulating RNA-binding proteins, as well as encoding proteins have been made very recently. Dysregulated circRNAs are associated with human diseases such as acute myeloid leukemia (AML). In this review, we focus on the recently described mechanisms, role and clinical significance of circRNAs in AML. Although great progress of circRNAs in AML has been achieved, substantial efforts are still required to explore whether circRNAs exert their biological function by other mechanisms such as regulation of gene transcription or serving as translation template in AML. It is also urgent that researchers study the machineries regulating circRNAs fate, the downstream effectors of circRNAs modulatory networks, and the clinical application of circRNAs in AML.

CircRNA Identification and CircRNA-miRNA-mRNA Network in Cynoglossus semilaevis Sexual Size Dimorphism

Sexual size dimorphism (SSD), which is the sexual differences in body size, has been widely reported in various species including fishes. For Chinese tongue sole (Cynoglossus semilaevis), a flatfish exhibiting typically female-biased SSD, little is known for its epigenetic regulation mechanism, especially the role of circRNAs. Here, we identified the differently expressed abundances of circRNAs in females, males, and pseudo-males to explore the potential functions of circRNAs in Chinese tongue sole SSD. In total, 14,745 novel circRNAs were screened, among which 1461 DE circRNAs were identified from the brain, gonad, liver, and muscle in female, male, and pseudo-male individuals. The ceRNA network was subsequently constructed, including 10 circRNAs, 26 mRNAs, and 11 miRNAs. These DE mRNAs were mainly related to the mRNA surveillance pathway, metabolic pathways, and cellular senescence. Importantly, the ceRNA network has revealed that several circRNAs such as novel_circ_004374 and novel_circ_014597 may regulate homeodomain interacting protein kinase 2 (hipk2) expression by sponging miR-130-x. It is also worth exploring whether or how novel_circ_008696 regulates SET Domain Containing 2, histone lysine methyltransferase (setd2), which in turn affects the epigenetic patterns of different sexual individuals. The present study not only enriches the knowledge on the potential roles of circRNA in the physiological process, but also provides new clues for the explanation of fish SSD. In future studies, the precise function and involvement of circRNAs in female-biased SSD will require more efforts.

Identification and functional analysis of circRNAs in the skeletal muscle of juvenile and adult largemouth bass (Micropterus salmoides)

Circular RNA (circRNA) is a novel emerging type of endogenous regulatory non-coding RNA molecules with a covalent closed-loop configuration, which exerts important functions in multiple biological processes. CircRNAs are known to regulate gene expression as functional regulators interacting with miRNAs by sponge, which have been reported to regulate skeletal muscle development. Nevertheless, the information of circRNAs involved in regulating muscle growth and development in fish is largely unknown. Here, we first identified 312 and 511 circRNAs in skeletal muscle of juvenile and adult largemouth bass (LMB) using RNA sequencing, respectively. The differentially expressed circRNAs (DE-circRNAs) analysis showed that there are 44 DE-circRNAs at two different skeletal muscle growth stages. Six circRNAs were chosen randomly and their relative expression levels in juvenile and adult LMB were confirmed by real-time PCR, indicating that these circRNAs were existed authenticity. In addition, Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analysis showed that these hose genes (their linear mRNAs) of DE-circRNAs were mainly enriched in the regulation of actin cytoskeleton signaling pathways. The circRNA-miRNA interaction regulatory networks indicated that one circRNA can regulate one or more miRNA. For instance, more than 30 miRNAs were regulated by two circRNAs (circRNA389 and circRNA399). Of them, the muscle-related miRNAs including the let-7 family, miR-133 and miR-26 and so on were found acting as miRNAs sponge regulated by circRNAs, indicating the roles of circRNAs in regulating muscle growth-related genes expression. Overall, these findings will not only broaden our understanding of circRNAs regulation mechanisms underlying muscle growth and development in LMB but also provides a novel clue for further functional research in carnivorous fish genetic breeding.