The circular RNA circBCL2L1 regulates innate immune responses via microRNA-mediated downregulation of TRAF6 in teleost fish

Hsa_circ_0008085 acts as a miR-146a-5p sponge to suppress influenza a virus replication via modulating of TRAF6

Influenza A virus (IAV) has attracted considerable attention in recent years due to its increase in incidence and threat to human health. Circular RNA (circRNA) is a non-coding RNA with a closed structural pattern, which often acts as a sponge for microRNAs (miRNAs) in the cell. A growing body of evidence supports several crucial roles for circRNAs in viral infection, including regulation of viral replication, evasion of the host immune response and disease pathogenesis. However, due to the wide variety of circRNAs, their potential functions in IAV infection remain to be elucidated. In this study, we determined that the expression of hsa_circ_0008085 was induced by IAV infection. In addition, the JAK/STAT signaling pathway was associated with the induction of hsa_circ_0008085 in the context of IAV infection, as evidenced by the application of ruxolitinib, a compound that inhibits the activation of the JAK/STAT pathway. In A549 cells, enhanced expression of hsa_circ_0008085 was found to repress viral gene expression and diminish the production of infectious progeny. In contrast, when hsa_circ_0008085 was knocked down, a significant boost in IAV replication was observed. Following a pull-down of RNA utilizing biotin and a luciferase reporter assay, we established that hsa_circ_0008085 binds with miR-146a-5p, functioning as an endogenous sponge that suppresses the activity of miR-146a-5p. This interaction leads to an upregulation of TRAF6 expression, which in turn inhibits the replication of IAV. Conclusively, these data demonstrate that hsa_circ_0008085 inhibits viral replication via the miR-146a-5p/TRAF6 axis, highlighting the promising role of circRNAs in antiviral therapies.

circRNA/TLR interaction: key players in immune regulation and autoimmune diseases

Circular RNAs are a class of non-coding RNAs with covalently closed loops. They have been revealed to regulate immune responses by affecting gene expression. Although initially considered splicing byproducts, new studies have indicated their role in transcriptional and post-transcriptional control, especially with TLRs. TLRs start inflammatory signaling and let the innate immune system recognize PAMPs. circRNAs interact context-dependently with TLR pathways to influence immune homeostasis and inflammation in either pathogenic or protective roles. In autoimmune diseases, dysregulated circRNA expression can aggravate immune responses and damage tissue. CircRNAs can interact with RNA-binding proteins, function as molecular sponges for miRNAs, and change inflammatory pathways like the NF-κB signaling cascade, influencing immune responses. They control adaptive immunity, function of antigen-presenting cells, and cytokine generation. The stability and presence of circRNAs in many body fluids make them therapeutic targets and biomarkers for inflammatory and autoimmune diseases. The several immune control roles of circRNA-TLR interactions are discussed in this review, as well as their consequences for immunologically mediated disease diagnosis and treatment.

Identification and functional regulation of two alternative splicing isoforms of the Uhrf2 gene in Miichthysmiiuy

Alternative splicing can produce a variety of splicing isoforms to increase protein diversity, participate in the regulation of gene expression and the occurrence and development of diseases, and thus play an important role in innate immunity. Ubiquitin like with PHD and ring finger domains 2 (Uhrf2) protein is associated with cell proliferation, inflammation, tumors, and cancer, and is currently the focus of medical immunology research, but there is little research on alternative splicing of the Uhrf2 gene. In this study, we identified two different splicing isoforms of Uhrf2 in Miichthys miiuy through Sanger sequencing, dual-luciferase reporter gene assay, qRT-PCR, subcellular localization experiments, and named them Uhrf2-α and Uhrf2-β. Subcellular localization experiments found that Uhrf2-α was mainly located in the nucleus, while Uhrf2-β was mainly located in the cytoplasm. Although their localization was different, both could significantly inhibit the activation of IRF3 and NF-κB signaling pathways, and effectively inhibit the levels of inflammatory cytokines. These results indicate that Uhrf2-α and Uhrf2-β play important negative regulatory roles in innate immune responses in fish.

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.

Identification of immune-responsive circular RNAs in shrimp (Litopenaeus vannamei) upon yellow head virus infection

Circular RNAs (circRNAs) are a subclass of non-coding RNAs (ncRNAs) formed through a process known as back-splicing. They play a crucial role in the genetic regulation of various biological processes. Currently, circRNAs have been identified as participants in the antiviral response within mammalian cells. However, circRNAs in shrimp infected with the yellow head virus (YHV) remain largely unexplored. Therefore, this study aims to identify circRNAs in the hemocytes of Litopenaeus vannamei during YHV infection. We discovered 358 differentially expressed circRNAs (DECs), with 177 of them being up-regulated and 181 down-regulated. Subsequently, eight DECs, including circ_alpha-1-inhibitor 3, circ_CDC42 small effector protein 2, circ_hemicentin 2, circ_integrin alpha V, circ_kazal-type proteinase inhibitor, circ_phenoloxidase 3, circ_related protein rab-8B, and circ_protein toll-like, were randomly selected for analysis of their expression patterns during YHV infection using qRT-PCR. Furthermore, the circRNAs' characteristics were confirmed through PCR, RNase R treatment, and Sanger sequencing, all of which were consistent with the features of circRNAs. These findings contribute to a better understanding of circRNAs' involvement in the antiviral response in shrimp.

Hsa_circ_0007321 regulates Zika virus replication through miR-492/NFKBID/NF-κB signaling pathway

IMPORTANCE

Over the past decade, increasing evidence has shown that circular RNAs (circRNAs) play important regulatory roles in viral infection and host antiviral responses. However, reports on the role of circRNAs in Zika virus (ZIKV) infection are limited. In this study, we identified 45 differentially expressed circRNAs in ZIKV-infected A549 cells by RNA sequencing. We clarified that a downregulated circRNA, hsa_circ_0007321, regulates ZIKV replication through targeting of miR-492 and the downstream gene NFKBID. NFKBID is a negative regulator of nuclear factor-κB (NF-κB), and we found that inhibition of the NF-κB pathway promotes ZIKV replication. Therefore, this finding that hsa_circ_0007321 exerts its regulatory role on ZIKV replication through the miR-492/NFKBID/NF-κB signaling pathway has implications for the development of strategies to suppress ZIKV and possibly other viral infections.

CircRNA254 functions as the miR-375 sponge to inhibit coelomocyte apoptosis via targeting BAG2 in V. splendidus-challenged Apostichopus japonicus

Circular RNAs (circRNAs) function as immune regulators in many biological processes in mammals, while their function and underlying mechanisms in invertebrates are largely unexplored. In this study, the competing endogenous RNA (ceRNA) mechanism of circRNA that sponges miR-375 and thus regulates AjBAG2-mediated coelomocyte apoptosis was evaluated in Apostichopus japonicus. The results showed that circRNA254 (circ254) was significantly down-regulated in the intestines and coelomocytes after Vibrio splendidus challenge or Lipopolysaccharide exposure, which matched the RNA-seq results in A. japonicus within skin ulceration syndrome. Dual-luciferase and RNA FISH assays indicated that circ254 could directly combine with miR-375, in which circ254 possesses three binding sites of miR-375. Moreover, circ254 knockdown significantly promoted the coelomocyte apoptosis levels upon pathogen infection in vivo and in vitro. Furthermore, circ254 silencing could also down-regulate AjBAG2 expression and thereby promoting the levels of coelomocyte apoptosis levels and the expression of caspase 3, which the phenomenon could be reversed by treatment with miR-375 inhibitors. Taken together, our results confirmed that circ254 functions as a ceRNA of AjBAG2 by sponging miR-375, resulting in the inhibition of coelomocyte apoptosis in A. japonicus.

Functional Involvement of circRNAs in the Innate Immune Responses to Viral Infection

Effective viral clearance requires fine-tuned immune responses to minimize undesirable inflammatory responses. Circular RNAs (circRNAs) are a class of non-coding RNAs that are abundant and highly stable, formed by backsplicing pre-mRNAs, and expressed ubiquitously in eukaryotic cells, emerging as critical regulators of a plethora of signaling pathways. Recent progress in high-throughput sequencing has enabled a better understanding of the physiological and pathophysiological functions of circRNAs, overcoming the obstacle of the sequence overlap between circRNAs and their linear cognate mRNAs. Some viruses also encode circRNAs implicated in viral replication or disease progression. There is increasing evidence that viral infections dysregulate circRNA expression and that the altered expression of circRNAs is critical in regulating viral infection and replication. circRNAs were shown to regulate gene expression via microRNA and protein sponging or via encoding small polypeptides. Recent studies have also highlighted the potential role of circRNAs as promising diagnostic and prognostic biomarkers, RNA vaccines and antiviral therapy candidates due to their higher stability and lower immunogenicity. This review presents an up-to-date summary of the mechanistic involvement of circRNAs in innate immunity against viral infections, the current understanding of their regulatory roles, and the suggested applications.

Chitooligosaccharide boosts the immunity of immunosuppressed blunt snout bream against bacterial infections

The immunosuppression hazard of fish brought by intensive aquaculture needs to be addressed urgently, while chitooligosaccharide (COS) shows the potential application in the prevention the immunosuppression of fish due to its superior biological properties. In this study, COS reversed the cortisol-induced immunosuppression of macrophages and improved the immune activity of macrophages in vitro, promoting the expression of inflammatory genes (TNF-α, IL-1β, iNOS) and NO production, and increasing the phagocytic activity of macrophages. In vivo, the oral COS was absorbed directly through the intestine, significantly ameliorating the innate immunity of cortisol-induced immunosuppression of blunt snout bream (Megalobrama amblycephala). Such as facilitated the gene expression of inflammatory cytokines (TNF-α, IL-1β, IL-6) and pattern recognition receptors (TLR4, MR) and potentiated bacterial clearance, resulting in an effective improvement in survival and tissue damage. Altogether, this study demonstrates that COS offers potential strategies in the application of immunosuppression prevention and control in fish.

Transcriptome analysis of differentially expressed circular RNAs in the testis and ovary of golden pompano (Trachinotus blochii)

The artificial breeding of golden pompano (Trachinotus blochii) has expanded greatly in recent years, and after long-term breeding efforts, clear sexual dimorphisms have been observed in T. blochii growth traits, with females growing faster. As sponges of microRNA (miRNAs), circular RNAs (CircRNAs) can alleviate miRNA inhibition of target mRNA. However, few studies have examined sex-related CircRNAs and none of those have looked at T. blochii. To further understand the role of CircRNAs in sex differentiation and sexual size dimorphism in T. blochii, six CircRNA libraries were constructed from the testes and ovaries of T. blochii. A total of 1522 CircRNAs were found distributed over all 24 chromosomes of T. blochii. 135 differentially expressed CircRNAs (DECs) were identified by screening, These DECs were then subjected to GO enrichment, which found 47 enriched pathways. A number of CircRNAs were enriched in cellular processes and metabolic processes. According to the KEGG pathway analysis, a series of sex differentiation pathways were enriched, including the GnRH, calcium, and MAPK signaling pathways. Furthermore, we selected two CircRNAs from the DECs named circ-cacna1b and circ-octc. We found that the cacna1b gene is regulated by 7 miRNAs, 3 of which were regulated by circ-cacna1b, i.e., mmu-miR-138-5p, fru-miR-138, and pma-miR-138b. In addition, the miRNA named pma-miR-138b can regulate sex-related genes, such as sox9 and dmrt1, among others. The co-expression network of CircRNA-miRNA-mRNA showed circ-cacna1b may play a crucial role in T. blochii sex differentiation by regulating pma-miR-138b to affect the expression of sex differentiation genes. The circ-octc may be one of the largest contributors to sexual size dimorphism during growth through its effect on lipid metabolism. These findings could broaden our understanding of CircRNAs and provide new insight into their function in sex differentiation and growth.

Circular RNA circPlce1 regulates innate immune response in miiuy croaker, Miichthys miiuy

In recent years, more and more researchers have devoted to the study of circular RNAs (circRNAs) in noncoding RNAs. As an important regulator in a variety of biological processes, circRNAs are relatively abundant in the study of mammals, while research in lower vertebrates is still lacking. In this study, we found a circRNA, circPlce1, related to innate immune response in Miichthys miiuy (miiuy croaker). The experimental results confirmed that circPlce1 could promote the production of antiviral genes and inflammatory response under the stimulation of poly (I: C) and LPS. We also confirmed that circPlce1 can promote NF-κB and IRF3 pathways through luciferase reporter assay experiment. In addition, we also found that circPlce1 can promote cell proliferation and improve cell viability. In conclusion, our results showed that circPlce1 plays an active role in regulating inflammatory response, cell proliferation and cell viability, providing a foundation for the study of the biological function of circRNAs in the innate immune response in teleost fish.

Circ004463 promotes fibroblast proliferation and collagen I synthesis by sponging miR-23b and regulating CADM3/MAP4K4 via activation of AKT/ERK pathways

Skin thickness is closely related to the appearance of human skin, such as sagging and wrinkling, which primarily depends on the level of collagen I synthesized by fibroblasts in the dermal layer. To explore the underlying genetics of the development of skin thickness, we used the indigenous Chinese Chenghua pigs, considered to have superior skin thickness, as model animals. We first performed whole transcriptome sequencing analysis to identify significant skin morphological differences between Chenghua pigs and Large White pigs and obtained some differentially expressed coding RNAs (454 mRNAs) and noncoding RNAs (612 circRNAs, 188 miRNAs, and 19 lncRNAs); moreover, some competing endogenous RNA (ceRNA) networks were constructed. Interestingly, we then identified a circRNA, namely circ0044633, which plays an important role in promoting fibroblast proliferation along with myofibroblast transition and collagen I synthesis by sponging miR-23b and regulating CADM3 and MAP4K4 expression via activation of the downstream AKT and ERK pathways in vitro. Furthermore, overexpression of circ004463 increased the mouse skin thickness and collagen I content in vivo. These results revealed a whole transcript profile of skin tissue and identified an important circ0044633-miR-23b-CADM3/MAP4K4 axis related to fibroblast proliferation and collagen I synthesis during the development of skin thickness.

The immune-related circRNA-miRNA-mRNA ceRNA regulatory network in the liver of turbot (Scophthalmus maximus L.) induced by Vibrio anguillarum

Recently, Vibrio anguillarum, a Gram-negative pathogenic bacterium, has been becoming a major constraint on the development of the turbot aquaculture industry because of its characteristics of worldwide distribution, broad host range and potentially devastating impacts. Although the functions of protein-coding mRNAs in the immune response against bacterial infection have been reported, as well as several non-coding RNAs (ncRNAs), such as circular RNAs (circRNAs) and microRNAs (miRNAs), the relationships between mRNAs and ncRNAs in the immune system of turbot liver are still limited during bacterial infection. In present study, the comprehensive analyses of whole-transcriptome sequencing were conducted in turbot liver infected by V. anguillarum. The differential expression was analyzed in the data of circRNAs, miRNAs, and mRNAs. The interactions of miRNA-circRNA pairs and miRNA-mRNA pairs were predicted basing on the negative regulatory relationships between miRNAs and their target circRNAs\mRNAs. The circRNA-related ceRNA regulatory networks were constructed for the analyses of regulated mechanism in turbot immune system. Subsequently, the RT-qPCR was carried out to verify the results of sequencing. Finally, we identified 31 circRNAs, 53 miRNAs and 948 mRNAs with differential expression. Gene set enrichment analyses using Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways showed that innate immunity was principally activated at the early stages of infection, while adaptive immunity was activated after 24 h. Finally, 65 circRNA-miRNA-mRNA pathways were constructed, based on the hypothesis of ceRNA regulatory networks. In conclusion, our findings provide new insights on the underlying immune response to bacterial infection and identify novel target genes for the prevention and control of disease in turbot.

Non-coding RNAs targeting NF-κB pathways in aquatic animals: A review

Nuclear factor-kappa B (NF-κB) pathways have a close relationship with many diseases, especially in terms of the regulation of inflammation and the immune response. Non-coding RNAs (ncRNAs) are a heterogeneous subset of endogenous RNAs that directly affect cellular function in the absence of proteins or peptide products; these include microRNAs (miRNAs), long noncoding RNAs (lncRNAs), circular RNAs (circRNAs), etc. Studies on the roles of ncRNAs in targeting the NF-κB pathways in aquatic animals are scarce. A few research studies have confirmed detailed regulatory mechanisms among ncRNAs and the NF-κB pathways in aquatic animals. This comprehensive review is presented concerning ncRNAs targeting the NF-κB pathway in aquatic animals and provides new insights into NF-κB pathways regulatory mechanisms of aquatic animals. The review discusses new possibilities for developing non-coding-RNA-based antiviral applications in fisheries.

High Dietary Histamine Induces Digestive Tract Oxidative Damage in Juvenile Striped Catfish (Pangasianodon hypophthalmus)

A 56-day feeding trial investigated the effects of dietary histamine on the antioxidant capacity, gastric and intestinal barrier functions, and growth performance of striped catfish (Pangasianodon hypophthalmus). Seven isonitrogenous (34.0% crude protein) and isolipidic (10.5% crude lipid) diets were formulated with supplemental 0, 15, 30, 60, 120, 240, and 480 mg/kg of histamine, named H0, H15, H30, H60, H120, H240, and H480 group, respectively. Results showed that the weight gain rate, specific growth rate, relative intestinal length in the H240 and H480 groups, and the condition factors in the H480 group were significantly lower than those in the H0 group. Intestinal total antioxidant capacity, peroxidase, catalase, superoxide dismutase, glutathione peroxidase, and glutathione reductase activities in the H480 group were significantly lower than those in the H0 group, whereas intestinal malondialdehyde content exhibited the opposite trend. Intestinal complement 3, complement 4, immunoglobulin M, and Recombinant Mucin 2 in the H480 group were significantly lower than those in the H0 group, in contrast to intestinal lipopolysaccharide content. Intestinal IL-10 gene expression in the H480 group was significantly lower than that in the H0 group, whereas the TNF-α, IL-1, IL-6, and IL-8 gene expression exhibited opposite results. Scanning and transmission electron microscopic observation of the gastrointestinal tract revealed severe damage to the gastric mucosa and intestinal epithelium in the H480 group. The abundance of Treponema in the histamine groups was significantly higher than that in the H0 group. These results indicated that high dietary histamine decreases intestinal immunity and antioxidant capacity, inducing digestive tract oxidative damage and ultimately decreasing the growth of striped catfish.

Assessing effects of guar gum viscosity on the growth, intestinal flora, and intestinal health of Micropterus salmoides

A 56-day feeding trial was conducted to assess the effects of different viscous guar gum on the growth, intestinal flora, and intestinal health of Micropterus salmoides. Four practical diets with 42.5 % crude protein and 13.7 % crude lipid were formulated to contain 8 % cellulose and three different viscosities (2500, 5200, and 6000 mPa·s) of guar gum. Dietary guar gum inhibits fish growth and feed utilization, decreases the α-diversity of the intestinal flora, and negatively alters the intestinal flora structure and metabolite composition. High viscous guar gum down-regulated the intestinal tight junction, anti-inflammatory, and anti-apoptotic related gene's expression, decreased digesta butyrate/histamine ratio; and increased the abundance of Plesiomonas shigelloides. These results suggest that dietary guar gum adversely affects intestinal health by disrupting intestinal flora structure and metabolite composition, and that viscosity should be considered when using guar gum as a binder in aquafeeds.

Molecular characterization and expression analysis of TRIF, TRAF6, and TBK1 of golden pompano (Trachinotus ovatus)

Toll/IL-1R domain-containing adaptor-inducing IFN-β (TRIF), tumor necrosis factor receptor-associated factor 6 (TRAF6) and TANK-binding kinase 1 (TBK1) are critical signal transducers in toll-like receptors (TLRs) signaling pathway. In the present study, TRIF, TRAF6 and TBK1 were characterized from golden pompano (Trachinotus ovatus), named as TroTRIF, TroTRAF6 and TroTBK1, respectively. The full cDNA length of TroTRIF, TroTRAF6 and TroTBK1 was 2297 bp, 2293 bp, and 2482 bp, which respectively encoded 589, 573 and 723 amino acids. The deduced amino acids sequences of TroTRIF, TroTRAF6 and TroTBK1 contained conserved motifs, similar to their counterparts in other vertebrates. Phylogenetic tree analysis revealed that TroTRIF, TroTRAF6 and TroTBK1 were well clustered with their counterparts in other fish species. Quantitative Real-Time PCR (qPCR) analysis showed that TroTRIF, TroTBK1 and TroTRAF6 were detected in all examined tissues of healthy fish, but shared distinct transcript levels. Moreover, the expressions of TroTRIF, TroTBK1 and TroTRAF6 were generally induced by polyriboinosinic-polyribocytidylic acid (polyI:C), lipopolysaccharide (LPS), and Vibrio alginolyticus stimulation in vivo, indicating their critical roles in the immune defense of golden pompano against pathogen invasion. Our results provide valuable information for understanding the functions of these genes in golden pompano.

circMEF2D Negatively Regulated by HNRNPA1 Inhibits Proliferation and Differentiation of Myoblasts via miR-486-PI3K/AKT Axis

Circular RNA (circRNA) is a form of endogenous RNA that can regulate gene expression and participate in the regulation of myogenesis. However, the molecular mechanisms and potential roles of circRNAs in bovine muscle development remain largely unknown. Nevertheless, the RNA splicing factors regulating the biogenesis of bovine circRNA have not yet been characterized. In this study, we identified a novel circRNA, circMEF2D, formed by back-splicing of constitutive exons (exons 5-7) of the bovine MEF2D gene. Functional assays showed that circMEF2D inhibited the proliferation and differentiation of bovine myoblasts. Importantly, we showed that circMEF2D regulated the PI3K-AKT signaling pathway through direct and competitive binding to miR-486. Furthermore, to explore the formation mechanism of circMEF2D, we explored the MEF2D gene alternative splicing progress. Four alternative linear variants of MEF2D were found. Due to its role in alternative splicing, the RNA-binding protein HNRNPA1 was selected for further study and the modulation of HNRNPA1 levels showed that it negatively regulated both back-splicing and linear splicing of MEF2D gene. Overall, in addition to the characterization of bovine circRNAs, these findings revealed the crucial role of HNRNPA1 in MEF2D gene alternative splicing and demonstrated a regulatory circMEF2D-miR-486-PI3K-AKT axis.

RNA editing enzyme ADAR1 controls miR-381-3p-mediated expression of multidrug resistance protein MRP4 via regulation of circRNA in human renal cells

Multidrug resistance–associated protein 4 (MRP4), a member of the C subfamily of ATP-binding cassette transporters, is highly expressed in the kidneys of mammals and is responsible for renal elimination of numerous drugs. Adenosine deaminase acting on RNA 1 (ADAR1) has been reported to regulate gene expression by catalyzing adenosine-to-inosine RNA editing reactions; however, potential roles of ADAR1 in the regulation of MRP4 expression have not been investigated. In this study, we found that downregulation of ADAR1 increased the expression of MRP4 in human renal cells at the posttranscriptional level. Luciferase reporter assays and microarray analysis revealed that downregulation of ADAR1 reduced the levels of microRNA miR-381-3p, which led to the corresponding upregulation of MPR4 expression. Circular RNAs (circRNAs) are a type of closed-loop endogenous noncoding RNAs that play an essential role in gene expression by acting as miRNA sponges. We demonstrate that ADAR1 repressed the biogenesis of circRNA circHIPK3 through its adenosine-to-inosine RNA editing activity, which altered the secondary structure of the precursor of circHIPK3. Furthermore, in silico analysis suggested that circHIPK3 acts as a sponge of miR-381-3p. Indeed, we found overexpression of circHIPK3 induced the expression of MRP4 through its interference with miR-381-3p. Taken together, our study provides novel insights into regulation of the expression of xenobiotic transporters through circRNA expression by the RNA editing enzyme ADAR1.

Identification of mRNA-miRNA-lncRNA regulatory network associated with the immune response to Aeromonas salmonicides infection in the black rockfish (Sebastes schlegelii)

In aquaculture, Aeromonas salmonicides (A. salmonicida) is a main fish pathogen because of its nearly worldwide distribution, and broad host range. Recently, an increasing number of evidences have uncovered the roles of mRNA-miRNA-lncRNA network in fish diseases. In current study, RNA-seq was conducted in the black rockfish spleen following A. salmonicida infection at 0 h (Sp0 or control) and three different post-infection time-points (2 h: Sp2, 12 h: Sp12 and 24 h: Sp24, respectively) to comprehensively identify differentially expressed (DE) mRNAs, miRNAs and lncRNAs. Enrichment analysis and protein-protein interaction (PPI) analysis of DE mRNAs were performed. Then, expression and correlation analysis for mRNAs and their upstream miRNAs and lncRNAs were conducted. Finally, a total of 1364 mRNAs, 17 miRNAs and 1584 lncRNAs exhibited significantly differential expressions during bacterial infection in the black rockfish spleen. Functional enrichment analysis suggested that they were significantly enriched in several immune-related pathways, including Amino sugar and nucleotide sugar metabolism, Cell adhesion molecules (CAMs), Neuroactive ligand-receptor interaction, Nicotinate and nicotinamide metabolism, Pentose and glucuronate interconversions, Phagosome, Proteasome, etc. Subsequently, 1091 lncRNA-miRNA-mRNA pathways (323 in Sp2, 609 in Sp12 and 207 in Sp24) were constructed including 400 lncRNAs, 69 miRNAs, and 70 mRNAs. Meanwhile, NLRC3/novel-264/LNC_00116154 pathway demonstrated important immune modulating function in the black rockfish against A. salmonicida infection. Finally, the novel mRNA-miRNA-lncRNA sub-networks were established, among which all mRNAs and ncRNAs possessed significant predictive values for further studies for immune responses in the black rockfish.