Differentially Expressed circRNAs Research Articles

CircN4bp1 Facilitates Sepsis-Induced Acute Respiratory Distress Syndrome through Mediating Macrophage Polarization via the miR-138-5p/EZH2 Axis.

We recently reported the differential circRNA expression patterns of the pulmonary macrophages in sepsis-induced acute respiratory distress syndrome (ARDS) mice model by microarray analysis. However, their function and hidden molecular mechanism in regulation of macrophage activation and inflammation remain poorly understood. In this study, we found that circN4bp1was overexpressed in PBMC and monocytes, and its expression levels were correlated with a poor prognosis in sepsis induced ARDS patients induced by sepsis. Knockdown of circN4bp1 inhibited the lung injury and improved the long-time survival through blunting the M1 macrophage activation in cecal ligation and puncture- (CLP-) induced ARDS mice. Moreover, bioinformatics analysis predicated a circN4bp1/miR-138-5p ceRNA network, which was confirmed by luciferase reporter assay and RNA binding protein immunoprecipitation (RIP). CircN4bp1 affected macrophage differentiation by binding to miR-138-5p, thus regulating the expression of EZH2 in vivo and ex vivo. Lastly, the m6A level of circN4bp1was found to be elevated in ARDS mice; inhibition of m6A methyltransferase METTL3 blocked this response in vitro. Therefore, circN4bp1 can function as a miR-138-5p sponge for the modulation of macrophage polarization through regulation the expression of EZH2 and may serve as a potential target and/or prognostic marker for ARDS patients following sepsis.

The exosome-circ_0001359 derived from cigarette smoke exposed-prostate stromal cells promotes epithelial cells collagen deposition and primary ciliogenesis

Cigarettes consumption is continued to be popular. We found that cigarette smoke (CS) exposure promoted prostatic fibrosis. In this study, human prostate epithelial RWPE-1 cells were co-cultured with exosomes derived from CS exposed-WPMY-1 cells (CS-WPMY-1-exo). The collagen deposition, primary ciliogenesis, epithelial-mesenchymal transition (EMT) and transforming growth factor (TGF)-β1 level of RWPE-1 were evaluated. The circRNAs profiles of WPMY-1-exo were explored by high-throughput RNA sequencing. It was found that CS-WPMY-1-exo significantly promoted RWPE-1 collagen deposition, EMT and primary ciliogenesis. There were 17 differentially expressed (DE) circRNAs (including circ_0001359) between CS-WPMY-1-exo and the negative control. Functional enrichment analyses showed that the DE circRNAs played important roles in ciliary basal body, spindle microtubule and TGF-β signaling pathway. Circ_0001359 siRNA attenuated CS-WPMY-1 induced RWPE-1 cells collagen deposition, EMT and primary ciliogenesis, as well as inhibited the level of TGF-β1. The whole results showed that circ_0001359 derived from CS-WPMY-1-exo contributed to prostatic fibrosis via stimulating epithelial cells phenotypes changes and collagen deposition.

Expression Profiles of Circular RNA in Aortic Vascular Tissues of Spontaneously Hypertensive Rats.

Background: Circular RNAs (circRNAs), as a kind of endogenous non-coding RNA, have been implicated in ischemic heart diseases and vascular diseases. Based on theirs high stability with a closed loop structure, circRNAs function as a sponge and bind specific miRNAs to exert inhibitory effects in heart and vasculature, thereby regulating their target gene and protein expression, via competitive endogenous RNA (ceRNA) mechanism. However, the exact roles and underlying mechanisms of circRNAs in hypertension and related cardiovascular diseases remain largely unknown.Methods and Results: High-throughput RNA sequencing (RNA-seq) was used to analyze the differentially expressed (DE) circRNAs in aortic vascular tissues of spontaneously hypertensive rats (SHR). Compared with the Wistar-Kyoto (WKY) rats, there were marked increases in the levels of systolic blood pressure, diastolic blood pressure and mean blood pressure in SHR under awake conditions via the tail-cuff methodology. Totally, compared with WKY rats, 485 DE circRNAs were found in aortic vascular tissues of SHR with 279 up-regulated circRNAs and 206 down-regulated circRNAs. Furthermore, circRNA-target microRNAs (miRNAs) and the target messenger RNAs (mRNAs) of miRNAs were predicted by the miRanda and Targetscan softwares, respectively. Additionally, real-time RT-PCR analysis verified that downregulation of rno_circRNA_0009197, and upregulation of rno_circRNA_0005818, rno_circRNA_0005304, rno_circRNA_0005506, and rno_circRNA_0009301 were observed in aorta of SHR when compared with that of WKY rats. Then, the potential ceRNA regulatory mechanism was constructed via integrating 5 validated circRNAs, 31 predicted miRNAs, and 266 target mRNAs. More importantly, three hub genes (NOTCH1, FOXO3, and STAT3) were recognized according to PPI network and three promising circRNA-miRNA-mRNA regulatory axes were found in hypertensive rat aorta, including rno_circRNA_0005818/miR-615/NOTCH1, rno_circRNA_0009197/ miR-509-5p/FOXO3, and rno_circRNA_0005818/miR-10b-5p/STAT3, respectively.Conclusions: Our results demonstrated for the first time that circRNAs are expressed aberrantly in aortic vascular tissues of hypertensive rats and may serve as a sponge linking with relevant miRNAs participating in pathogenesis of hypertension and related ischemic heart diseases via the circRNA-miRNA-mRNA ceRNAnetwork mechanism.

Integral Analyses of Competing Endogenous RNA Mechanisms and DNA Methylation Reveal Regulatory Mechanisms in Osteosarcoma.

Background: Osteosarcoma (OS) is the most common primary malignant bone tumour in children and adolescents, with rapid growth, frequent metastasis, and a poor prognosis, but its pathogenesis has not been fully elucidated. Exploring the pathogenesis of OS is of great significance for improving diagnoses and finding new therapeutic targets. Methods: Differentially expressed circRNAs (DECs), miRNAs (DEMs), methylated DNA sites (DMSs), and mRNAs (DEGs) were identified between OS and control cell lines. GSEA of DEGs and functional enrichment analysis of methylated DEGs were carried out to further identify potential biological processes. Online tools were used to predict the miRNA binding sites of DECs and the mRNA binding sites of DEMs, and then construct a circRNA-miRNA-mRNA network. Next, an analysis of the interaction between methylated DEGs was performed with a protein-protein interaction (PPI) network, and hub gene identification and survival analysis were carried out. The expression pattern of circRNA-miRNA-mRNA was validated by real-time PCR. Results: GSEA and functional enrichment analysis indicated that DEGs and methylated DEGs are involved in important biological processes in cancer. Hsa_circ_0001753/has_miR_760/CD74 network was constructed and validated in cell lines. Low expression levels of CD74 are associated with poor overall survival times and show good diagnostic ability. Conclusion: Methylated DEGs may be involved in the development of OS, and the hsa_circ_0001753/has_miR_760/CD74 network may serve as a target for the early diagnosis of and targeted therapy for OS.

Construction and analysis of the ceRNA network hsa_circ_0031968/miR-3611/GCG in lung adenocarcinoma.

BackgroundA competitive endogenous RNA (ceRNA) network was constructed to examine the potential mechanisms of circular RNAs (circRNAs) in lung adenocarcinoma (LUAD).MethodsLUAD-related data sets were obtained from the Gene Expression Omnibus (GEO) database and screened for differentially expressed circRNAs (DECs) and differentially expressed microRNAs (DEMs). We identified the target microRNAs (miRNAs) regulated by the DECs and the potential target genes of the miRNA. The basic structure of the DECs were analyzed and enrichment analyses were conducted to determine the function of the circRNA. The Kaplan-Meier method for survival analysis was used to examine the clinical data from The Cancer Genome Atlas (TCGA) database. A protein-protein interaction (PPI) network was constructed to determine the hub genes. The relative expression of the RNA molecules on the ceRNA axis was verified by quantitative real-time polymerase chain reaction (qRT-PCR) and Western blot analysis.ResultsA total of 17 DECs and 237 DEMs were selected for analysis. After reviewing the cancer-specific circRNA database (CSCD), 10 circRNAs were identified. The 432 target miRNAs were screened by circRNA interactome (CRI) and cross-referenced with the DEMs to obtain 126 miRNAs of interest. The expression of miR-3611, which is regulated by hsa_circ_0031968, was found to significantly affect the survival and prognosis of patients with LUAD (P≤0.05). The target gene function of hsa_circ_0031968 was determined to be mainly enriched in SMAD binding, and the signaling pathway was primarily enriched in miRNAs related to cancer. The TCGA database screened out 2,484 differentially expressed mRNAs (DEmRNAs) and intersection analysis with the target gene of miR-3611 revealed 1 gene, namely the proglucagon gene (GCG). Therefore, we chose the hsa_circ_0031968/miR-3611/GCG axis for further research. The expression of GCG was determined to be associated with a poorer survival rate and higher T stage in LUAD patients. Finally, 17 hub genes that interact with GCG were identified.ConclusionsThe ceRNA regulatory network hsa_circ_0031968/miR-3611/GCG was successfully constructed and this provided novel insights into the identification of biomarkers and the pathogenesis of LUAD. This knowledge will contribute to the early diagnosis and development of potential treatment for patients with LUAD.

Analysis of the Differential Expression of circRNA in Acute Myeloid Leukemia by GEO Database

To investigate the difference expression of circular RNA (circRNA) in acute myeloid leukemia (AML) by using bioinformatics method. The microarray chip data of AML was searched and downloaded from the Gene Expression Omnibus (GEO) of the National Center for Bioinformatics (NCBI). The differences between AML samples and control samples were analyzed by R software. The interaction between deregulated circRNA, miRNA and mRNA were predicted by miranda software and miRTarBase software. The circRNA-miRNA-mRNA regulatory network was constructed by using the cytoHubba plugin based on the Cytoscape software. A total of 203 differential expression of circRNAs were finally collected, including down-regulated 161 circRNAs and up-regulated 42 circRNAs. CircRNA/miRNA/mRNA interaction network was constructed through software prediction. hsa_circ_0001080, hsa_circ_0004511, hsa_circ_0054211, hsa_circ_0001944 may be positively regulated the gene expression in AML. Abnormal expression of circRNA in AML may become a new target for AML treatment.

The circular RNA expression profile in ovarian serous cystadenocarcinoma reveals a complex circRNA\u2013miRNA regulatory network

BackgroundOvarian serous cystadenocarcinoma is one of the most serious gynecological malignancies. Circular RNA (circRNA) is a type of noncoding RNA with a covalently closed continuous loop structure. Abnormal circRNA expression might be associated with tumorigenesis because of its complex biological mechanisms by, for example, functioning as a microRNA (miRNA) sponge. However, the circRNA expression profile in ovarian serous cystadenocarcinoma and their associations with other RNAs have not yet been characterized. The main purpose of this study was to reveal the circRNA expression profile in ovarian serous cystadenocarcinoma.MethodsWe collected six specimens from three patients with ovarian serous cystadenocarcinoma and adjacent normal tissues. After RNA sequencing, we analyzed the expression of circRNAs with relevant mRNAs and miRNAs to characterize potential function.Results15,092 unique circRNAs were identified in six specimens. Approximately 46% of these circRNAs were not recorded in public databases. We then reported 353 differentially expressed circRNAs with oncogenes and tumor-suppressor genes. Furthermore, a conjoint analysis with relevant mRNAs revealed consistent changes between circRNAs and their homologous mRNAs. Overall, construction of a circRNA–miRNA network suggested that 4 special circRNAs could be used as potential biomarkers.ConclusionsOur study revealed the circRNA expression profile in the tissues of patients with ovarian serous cystadenocarcinoma. The differential expression of circRNAs was thought to be associated with ovarian serous cystadenocarcinoma in the enrichment analysis, and co-expression analysis with relevant mRNAs and miRNAs illustrated the latent regulatory network. We also constructed a complex circRNA–miRNA interaction network and then demonstrated the potential function of certain circRNAs to aid future diagnosis and treatment.

Analysis and validation of m6A regulatory network: a novel circBACH2/has-miR-944/HNRNPC axis in breast cancer progression

BackgroundN6-methyladenosine (m6A), the most abundant and reversible modification of mRNAs in eukaryotes, plays pivotal role in breast cancer (BC) tumorigenesis and progression. Circular RNAs (circRNAs) can act as tumor promoters or suppressors by microRNA (miRNA) sponges in BC. However, the underlying mechanism of circRNAs in BC progression via regulating m6A modulators remains unclear.MethodsPrognostic m6A RNA methylation regulators were identified in 1065 BC patients from The Cancer Genome Atlas (TCGA) project. Differentially expressed (DE) miRNAs and DE circRNAs were identified between BC and normal samples in TCGA and GSE101123, respectively. MiRNA-mRNA interactive pairs and circRNA-miRNA interactive pairs were verified by MiRDIP and Circular RNA Interactome. GSEA, KEGG, and ssGSEA were executed to explore the potential biological and immune functions between HNRNPC-high and HNRNPC-low expression groups. qRT-PCR and Western blot were used to quantify the expression of HNRNPC and circBACH2 in MCF-7 and MDA-MB-231 cells. The proliferation of BC cells was assessed by CCK-8 and EdU assay.Results2 m6A RNA methylation regulators with prognostic value, including HNRNPC and YTHDF3, were identified in BC patients. Then, the regulatory network of circRNA-miRNA-m6A modulators was constructed, which consisted of 2 DE m6A modulators (HNRNPC and YTHDF3), 12 DE miRNAs, and 11 DE circRNAs. Notably, BC patients with high expression of HNRNPC and low expression of hsa-miR-944 were correlated with late clinical stages and shorter survival times. Besides, the results from the KEGG inferred that the DE HNRNPC was associated with the MAPK signaling pathway in BC. Moreover, the circBACH2 (hsa_circ_0001625) was confirmed to act as hsa-miR-944 sponge to stimulate HNRNPC expression to promote BC cell proliferation via MAPK signaling pathway, thus constructing a circBACH2/hsa-miR-944/HNRNPC axis in BC.ConclusionsOur findings decipher a novel circRNA-based m6A regulatory mechanism involved in BC progression, thus providing attractive diagnostic and therapeutic strategies for combating BC.

Analysis on the expression profile of circRNAs in hypertrophic myocardium mice

Objective: To explore the differential expression of circRNAs and their potential impact on the pathophysiological process in cardiac hypertrophy. Methods: Six SPF C57BL/6J male mice, aged 8 to 10 weeks, were randomly divided into transverse aortic constriction (TAC) group (n=3) or sham operation(sham) group (n=3) according to random number table method. TAC mouse model was used to induce cardiac hypertrophy. Four weeks after surgery, high-throughput sequencing analysis was performed to detect differentially expressed circRNA in left myocardial tissues of mice between TAC group and sham group, and principal component analysis of circRNA was performed by R language software. Enrichment analysis was performed by GO and KEGG databases to predict the basic functions of differentially expressed circRNA-derived genes and their biological pathways. The differentially expressed circRNAs in the sequencing results were verified by real-time fluorescence quantitative polymerase chain reaction. Cytoscape software was used to construct circRNA-microRNA (miRNA) network maps to predict their interactions by combining differentially expressed circRNA and TargetScan predicted miRNA sites. Results: Principal component analysis was performed on 4 580 circRNAs detected from 6 samples of mice in TAC group and sham group. The results of R language software indicated that the variance contribution rate of the first 3 principal components, namely the first, second and third principal components, was 91.01%, 3.19% and 2.01%, respectively, and the cumulative variance contribution rate of the 3 components was 96.21%. Among the differentially expressed circRNAs, 6 (19%) were up-regulated and 25 (81%) were down-regulated in the TAC group. GO analysis showed that differentially expressed circRNA was closely related to the occurrence and development of cardiac hypertrophy, and KEGG pathway analysis suggested that downregulated circRNA expression was involved in the regulation of actin cytoskeleton. Fifteen out of the 31 differentially expressed circRNAs were selected for real-time fluorescence quantitative polymerase chain reaction verification, and the results showed that 8 circRNAs were consistent with sequencing results. circRNA-miRNA co-expression network analysis results showed that chr11:65218529-65233184-interacts with mmu-miRNA-30e-3p and mmu-miRNA-30a-3p. Conclusions The differential expression of circRNA in hypertrophic myocardium mice is evidenced in TAC mouse model. circRNA may interact with the corresponding miRNA to influence the occurrence and development of cardiac hypertrophy through autophagy-related cellular hypertrophy pathway or apoptosis-related pathological phenotypes.

Genome-Wide RNA-Sequencing Reveals Massive Circular RNA Expression Changes of the Neurotransmission Genes in the Rat Brain after Ischemia-Reperfusion.

Ischemic brain stroke is one of the most serious and socially significant diseases. In addition to messenger RNAs (mRNAs), encoding protein, the study of regulatory RNAs in ischemic has exceptional importance for the development of new strategies for neuroprotection. Circular RNAs (circRNAs) have a closed structure, predominantly brain-specific expression, and remain highly promising targets of research. They can interact with microRNAs (miRNAs), diminish their activity and thereby inhibit miRNA-mediated repression of mRNA. Genome-wide RNA-Seq analysis of the subcortical structures of the rat brain containing an ischemic damage focus and penumbra area revealed 395 circRNAs changed their expression significantly at 24 h after transient middle cerebral artery occlusion model (tMCAO) conditions. Furthermore, functional annotation revealed their association with neuroactive signaling pathways. It was found that about a third of the differentially expressed circRNAs (DECs) originate from genes whose mRNA levels also changed at 24 h after tMCAO. The other DECs originate from genes encoding non-regulated mRNAs under tMCAO conditions. In addition, bioinformatic analysis predicted a circRNA–miRNA–mRNA network which was associated with the neurotransmission signaling regulation. Our results show that such circRNAs can persist as potential miRNA sponges for the protection of mRNAs of neurotransmitter genes. The results expanded our views about the neurotransmission regulation in the rat brain after ischemia–reperfusion with circRNA action.

Identification of circular RNAs and corresponding regulatory networks reveals potential roles in the brains of honey bee workers exposed to dinotefuran

Honey bees are important and highly efficient pollinators of agricultural crops and have been negatively affected by insecticides in recent years. Circular RNA (circRNA) plays an important role in the regulation of multiple biological and pathological processes; however, its role in the honey bee brain after exposure to dinotefuran is not well understood. Here, the expression profiles and potential modulation networks of circRNAs in the brains of workers (Apis mellifera) were comprehensively investigated using RNA sequencing and bioinformatics. In total, 33, 144, and 211 differentially expressed (DE) circRNAs were identified on the 1st, 5th and 10th days after exposure to dinotefuran, respectively. Enrichment analyses revealed that the host genes of DE circRNAs were enriched in the Hippo signaling pathway-fly, Wnt signaling pathway, and neuroactive ligand-receptor interaction. circ_0002266, circ_0005080, circ_0010239 and circ_0005415 were found to have translational potential due to the presence of an internal ribosome entry site (IRES). An integrated analysis of the DE circRNA-miRNA-mRNA networks suggest that circ_0008898 and circ_0001829 may participate in the immune response to dinotefuran exposure by acting as miRNA sponges. Our results provide invaluable basic data on A. mellifera brain circRNA patterns and a molecular basis for further study of the biological function of circRNAs in the development and immune response of honey bees.

Circular RNA Expression and Regulation Profiling in Testicular Tissues of Immature and Mature Wandong Cattle (Bos taurus).

Wandong cattle are an autochthonous Chinese breed used extensively for beef production. The breed tolerates extreme weather conditions and raw feed and is resistant to tick-borne diseases. However, the genetic basis of testis development and sperm production as well as breeding management is not well established in local cattle. Therefore, improving the reproductive efficiency of bulls via genetic selection is crucial as a single bull can breed thousands of cows through artificial insemination (AI). Testis development and spermatogenesis are regulated by hundreds of genes and transcriptomes. However, circular RNAs (circRNAs) are the key players in many biological developmental processes that have not been methodically described and compared between immature and mature stages in Bovine testes. In this study, we performed total RNA-seq and comprehensively analyzed the circRNA expression profiling of the testis samples of six bulls at 3 years and 3 months of developmental age. In total, 17,013 circRNAs were identified, of which 681 circRNAs (p-adjust < 0.05) were differentially expressed (DE). Among these DE circRNAs, 579 were upregulated and 103 were downregulated in calf and bull testes. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses revealed that the identified target genes were classified into three broad functional categories, including biological process, cellular component, and molecular function, and were enriched in the lysine degradation, cell cycle, and cell adhesion molecule pathways. The binding interactions between DE circRNAs and microRNAs (miRNAs) were subsequently constructed using bioinformatics approaches. The source genes ATM, CCNA1, GSK3B, KMT2C, KMT2E, NSD2, SUCLG2, QKI, HOMER1, and SNAP91 were found to be actively associated with bull sexual maturity and spermatogenesis. In addition, a real-time quantitative polymerase chain reaction (RT-qPCR) analysis showed a strong correlation with the sequencing data. Moreover, the developed model of Bovine testes in the current study provides a suitable framework for understanding the mechanism of circRNAs in the development of testes and spermatogenesis.

Analysis of differentially expressed circRNAs in longissimus muscle between castrated and intact male pigs.

Castration can reduce odor and fights in boars, but the carcass yield is reduced, and the intramuscular fat content is increased. Understanding its molecular mechanism is of great significance for production. Recent studies have shown that circular RNAs (circRNAs) play an important role(s) in the regulation of muscle development. To explore the effects of circRNAs on the development of longissimus dorsi (LD) muscle after castration, six Huainan male pigs were selected and three of which were randomly castrated. Six pigs were slaughtered when their body weight reached around 130 kg, and the LD muscle samples were collected. The differentially expressed circRNAs (DECs) were screened by high-throughput sequencing and functionally analyzed using the KEGG databases. DECs-miRNAs network was constructed, and the expression profiles of candidate circRNAs and their interactions with miRNAs were verified in porcine skeletal muscle satellite cells. The results showed that a total of 5866 circRNAs were obtained, and 370 DECs were identified in LD muscle between the castrated and intact groups (| log2Foldchange | > 1, padj <0.8). KEGG enrichment indicated that the parental genes for the DECs were mainly enriched in the pathways associated with muscle development, muscle fiber type transformation, and energy metabolism. There were 8 miRNAs and 69 circRNAs enriched in the DECs-miRNA network. circRNA_2241 and circRNA_4237 were selected for verification, which showed that these two circRNAs really existed and their expression profiles were consistent with the sequencing results. Further, preliminary analysis showed that circRNA_2241 interacted with miR-1, and testosterone promoted circRNA_2241 but inhibited miR-1 expression. These results confirmed that circRNAs might participate in the regulation of LD muscle development after castration by interacting with miRNAs, thereby providing new materials and references for analyses on the molecular mechanisms of castration on the regulation of muscle development.

Comprehensive analysis of ceRNA networks to determine genes related to prognosis, overall survival, and immune infiltration in clear cell renal carcinoma

BackgroundClear cell renal cell carcinoma (ccRCC) is one of the common subtypes of kidney cancer. Circular RNAs (circRNAs) act as competing endogenous RNAs (ceRNAs) to affect the expression of microRNAs (miRNAs), and hence the expression of genes involved in the development and progression of ccRCC. However, these interactions have not been sufficiently explored. MethodsThe differential expression of circRNAs (DEC) was extracted from the GEO database, and the expression of circRNAs was analyzed by the Limma R package. The interaction of miRNAs with circRNAs was predicted using (cancer-specific circRNA database) CSCD and circinteractome database. The genes affected by the miRNAs were predicted by miRwalk version 3, and the differential expression was retrieved using TCGA. Functional enrichment was assessed and a PPI network was created using DAVID and Cytoscape, respectively. The genes with significant interactions (hub-genes) were screened, and the total survival rate of ccRCC patients was extracted from the Gene Expression Profiling Interactive Analysis (GEPIA) database. To confirm the expression of OS genes we used the Immunohistochemistry (IHC) data and TCGA database. The correlation between gene expression and immune cell infiltration was investigated using TIMER2.0. Finally, potential drug candidates were predicted by the cMAP database. ResultsFour DECs (hsa_circ_0003340, hsa_circ_0007836, hsa_circ_0020303, and hsa_circ_0001873) were identified, along with 11 interacting miRNAs (miR-1224-3p, miR-1294, miR-1205, miR-1231, miR-615-5p, miR-940, miR-1283, and miR-1305). These miRNAs were predicted to affect 1282 target genes, and function enrichment was used to identify the genes involved in cancer biology. 18 hub-genes (CCR1, VCAM1, NCF2, LAPTM5, NCKAP1L, CTSS, BTK, LILRB2, CD53, MPEG1, C3AR1, GPR183, C1QA, C1QC, P2RY8, LY86, CYBB, and IKZF1) were identified from a PPI network. VCAM1, NCF2, CTSS, LILRB2, MPEG1, C3AR1, P2RY8, and CYBB could affect the survival of ccRCC patients. The hub-gene expression was correlated with tumor immune cell infiltration and patient prognosis. Two potantial drug candidates, naphazoline and lithocholic acid could play a role in ccRCC therapy, as well other cancers. ConclusionThis bioinformatics analysis brings a new insight into the role of circRNA/miRNA/mRNA interactions in ccRCC pathogenesis, prognosis, and possible drug treatment or immunotherapy.

Abstract 10192: Circular RNA Formation in the Failing Heart is Enhanced by the Reduction of Adenosine-to-Inosine RNA Editing

Introduction: Adenosine-to-Inosine (A-to-I) RNA editing is a post-transcriptional modification process regulating RNA stability and alternative splicing. A-to-I RNA editing is conducted by the enzymes ADAR1 and ADAR2 and mainly targets Alu elements, primate-specific elements which have been associated with the formation of circular RNA (circRNA). Although differential expression of circRNAs has been studied in heart failure (HF), the extent of A-to-I RNA editing and consequences in the human heart remain largely unknown. Methods and Results: We analyzed RNA editing in human heart samples of HF (n=20) patients and controls (n=10) using RNA sequencing. We found a reduction of A-to-I RNA editing in intronic Alu elements of protein-coding genes in HF patients compared to controls. The majority (96%) of regulated circRNAs were upregulated. The predicted back-splice sites (BSS) of 20 circRNAs were validated by qPCR. The circRNA candidates correlated with RNA editing (R=0.47, P=0.02). Among the upregulated circRNAs, we identified two circular transcripts (circAKAP13) derived from the AKAP13 gene, which showed reduced A-to-I RNA editing in HF (-70.7%, n=20). In HF, ADAR2 was reduced (-68.2%) and ADAR1 was increased (7.41±0.13 -fold) on protein level (n=3-6). The knockdown of ADAR1 did not alter circRNA levels, whereas the knockdown of ADAR2 led to significantly upregulated levels of circAKAP13 (1.88±0.42 -fold, n=6). Consistently, ADAR2 overexpression reduced circAKAP13 expression (-41%, n=3). Using two mini-genes containing exons 15-19 of the AKAP13 gene and flanking Alu elements, we found convergent Alu elements enhancing circAKAP13 expression. Conclusion: In conclusion, these data describe the A-to-I RNA editome in the human heart for the first time. Reduced A-to-I RNA editing in HF patients is associated with elevated circRNA levels. We propose a primate-specific splicing mechanism mediated by A-to-I RNA editing in the human heart. These findings contribute to a better mechanistic understanding of A-to-I RNA editing in cardiac diseases.

Circular RNA screening identifies circMYLK4 as a regulator of fast/slow myofibers in porcine skeletal muscles

The type of myofiber is related to the quality of meat. The slow oxidized myofiber helps to increase the tenderness and juiciness of muscle. Numerous studies have shown that circRNA plays a key role in skeletal muscle development. However, the role of circRNA in porcine skeletal myofiber types is unclear. In this study, we performed high-throughput RNA sequencing to study the differential expression of circRNA in the longissimus dorsi and the soleus muscle. A total of 40,757 circRNAs were identified, of which 181 were significantly different. Interestingly, some circRNAs were involved in metabolism pathways, AMPK, FoxO, and PI3K-Akt signaling pathways. Besides, we focused on a novel circRNA-circMYLK4. By injecting circMYLK4-AAV into piglets, we found that circMYLK4 significantly increased the mRNA and protein levels of the slow muscle marker genes. In summary, our study laid an essential foundation for further research of circRNA in myofiber type conversion and higher meat quality.

Whole-transcriptome sequencing reveals a vernalization-related ceRNA regulatory network in chinese cabbage (Brassica campestris L. ssp. pekinensis)

BackgroundThe transition from vegetative growth to reproductive growth involves various pathways. Vernalization is a crucial process for floral organ formation and regulation of flowering time that is widely utilized in plant breeding. In this study, we aimed to identify the global landscape of mRNAs, microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs) related to vernalization in Chinese cabbage. These data were then used to construct a competitive endogenous RNA (ceRNA) network that provides valuable information to better understand the vernalization response.ResultsIn this study, seeds sampled from the Chinese cabbage doubled haploid (DH) line ‘FT’ with or without vernalization treatment were used for whole-transcriptome sequencing. A total of 2702 differentially expressed (DE) mRNAs, 151 DE lncRNAs, 16 DE circRNAs, and 233 DE miRNAs were identified in the vernalization-treated seeds. Various transcription factors, such as WRKY, MYB, NAC, bHLH, MADS-box, zinc finger protein CONSTANS-like gene, and B3 domain protein, and regulatory proteins that play important roles in the vernalization pathway were identified. Additionally, we constructed a vernalization-related ceRNA–miRNA–target gene network and obtained 199 pairs of ceRNA relationships, including 108 DEmiRNA‒DEmRNA, 67 DEmiRNA‒DElncRNA, and 12 DEmiRNA‒DEcircRNA interactions, in Chinese cabbage. Furthermore, several important vernalization-related genes and their interacting lncRNAs, circRNAs, and miRNAs, which are involved in the regulation of flowering time, floral organ formation, bolting, and flowering, were identified.ConclusionsOur results reveal the potential mRNA and non-coding RNAs involved in vernalization, providing a foundation for further studies on the molecular mechanisms underlying vernalization in Chinese cabbage.

Construction of the circRNA-miRNA-mRNA Regulatory Network of an Abdominal Aortic Aneurysm to Explore Its Potential Pathogenesis

Background Abdominal aortic aneurysm (AAA) is a progressive cardiovascular disease, which is a permanent and localized dilatation of the abdominal aorta with potentially fatal consequence of aortic rupture. Dysregulation of circRNAs is correlated with the development of various pathological events in cardiovascular diseases. However, the function of circRNAs in abdominal aortic aneurysm (AAA) is unknown and remains to be explored. This study is aimed at determining the regulatory mechanisms of circRNAs in AAAs. This study was aimed at exploring the underlying molecular mechanisms of abdominal aortic aneurysms based on the competing endogenous RNA (ceRNA) regulatory hypothesis of circRNA, miRNA, and mRNA. Methods The expression profiles of circRNAs (GSE144431), miRNAs (GSE62179), and mRNAs (GSE7084, GSE57691, and GSE47472) in human tissue sample from the aneurysm group and normal group were obtained from the Gene Expression Omnibus database, respectively. The circRNA-miRNA-mRNA network was constructed by using Cytoscape 3.7.2 software; then, the protein-protein interaction (PPI) network was constructed by using the STRING database, and the hub genes were identified by using the cytoHubba plug-in. The circRNA-miRNA-hub gene regulatory subnetwork was formed to understand the regulatory axis of hub genes in AAAs. Results The present study identified 40 differentially expressed circRNAs (DECs) in the GSE144431, 90 differentially expressed miRNAs (DEmiRs) in the GSE62179, and 168 differentially expressed mRNAs (DEGs) with the same direction regulation (130 downregulated and 38 upregulated) in the GSE7084, GSE57691, and GSE47472 datasets identified regarding AAAs. The miRNA response elements (MREs) of three DECs were then predicted. Four overlapping miRNAs were obtained by intersecting the predicted miRNA and DEmiRs. Then, 17 overlapping mRNAs were obtained by intersecting the predicted target mRNAs of 4 miRNAs with 168 DEGs. Furthermore, the circRNA-miRNA-mRNA network was constructed through 3 circRNAs, 4 miRNAs, and 17 mRNAs, and three hub genes (SOD2, CCR7, and PGRMC1) were identified. Simultaneously, functional enrichment and pathway analysis were performed within genes in the circRNA-miRNA-mRNA network. Three of them (SOD2, CCR7, and PGRMC1) were suggested to be crucial based on functional enrichment, protein-protein interaction, and ceRNA network analysis. Furthermore, the expression of SOD2 and CCR7 may be regulated by hsa_circ_0011449/hsa_circ_0081968/hsa-let-7f-5p; the expression of PGRMC1 may be regulated by hsa_circ_0011449/hsa_circ_0081968-hsa-let-7f-5p/hsa-let-7e-5p. Conclusion In conclusion, the ceRNA interaction axis we identified may be an important target for the treatment of abdominal aortic aneurysms. This study provided further understanding of the potential pathogenesis from the perspective of the circRNA-related competitive endogenous RNA network in AAAs.

Circular RNA Dysregulation Characterizes Symptomatic T-LGL Leukemia Patients with STAT3 Mutation

T-cell large granular lymphocyte leukemia (T-LGLL) is a rare disease characterized by clonal proliferation of CD3+ lymphocytes whose pathogenesis is not completely understood. LGL resistance to the activation-induced cell death is sustained by anti-apoptotic pathways, mainly due to constitutive activation of the JAK-STAT axis by activating somatic mutations of STAT3 or STAT5B genes or epigenetic disturbances. STAT3 mutations are peculiar to CD8+ T-LGLL and are associated with cytopenias, treatment requirement, and reduced survival. STAT5B lesions are mostly found in CD4+ T-LGLL, a disorder presenting with a mild clinical course. A similar behaviour has been observed also in both STAT wild type (WT) CD8+ and CD4+ forms. The link between STAT3-dependent miR-146b repression and neutropenia exemplifies the role that the pathogenetic effect of somatic mutations can exert through complex axes, also involving non-coding RNAs. Since in the most symptomatic form of LGLL an urgent need to discover new mechanisms in place is demanded, we focused on additional regulatory molecules such as circular RNAs (circRNAs), single-stranded covalently closed RNA molecules generated by backsplicing (Figure 1A). CircRNAs' oncogenic potential is linked to their activity as miRNA sponges, modulators of RNA-binding protein activity, and their ability to encode oncogenic peptides.The study group includes 20 patients recruited at the Hematology Unit of Padua University Hospital (Italy) and diagnosed with T-LGLL according to the WHO criteria, plus 5 healthy controls (CTR). Patient leukemic clones were equally distributed among CD8+ WT and STAT3-mutated and CD4+ WT and STAT5B-mutated. RNA-seq of ribodepleted RNA from immunomagnetically purified CD57+ leukemic clones and CD8+ cells from CTR was obtained to study circRNAs. CircRNAs were identified and quantified from RNA-seq data by CirComPara v0.6.3. Differential expression was assessed by edgeR with robust estimation of dispersion. The circular to linear proportion (CLP), indicating the relative abundance of each circRNA with respect to all the overlapping host gene transcripts, was calculated, and differential CLP across conditions was assessed with CircTest.We detected 68.619 circRNAs and focused on 5.948 with high expression. These were expressed from 2,940 loci, mostly from annotated genes (99.95%). We identified 28 circRNAs derived from genomic regions without known genes. About one-half of circRNA host genes expressed multiple, up to 20, circular isoforms. About 4% of circRNAs had a CLP of at least 0.25, with 95 circRNAs more expressed than the linear counterpart, including circGUSBP2, circSHOC1, circTCEANC, circZNF609, and circKLHL8.Unsupervised analysis of circRNA expression showed a clear separation of T-LGLL from CTR samples, pointing toward circRNA dysregulation in T-LGLL, particularly for the CD8+ STAT3 mutated group, which diverged from the other subtypes (Figure 1B). The direct comparison between the CD8+ STAT3-mutated, the three molecular subtypes as a single group (OTH), and the CTR samples identified differences in circRNA expression specifically linked to the most symptomatic group with STAT3 lesions and those commonly observed in T-LGLL. In STAT3-mutated, 500 circRNAs were dysregulated with an excess (71%) of upregulation compared to CTR (Figure 1C). Most circRNAs commonly altered in both LGLL groups had a more marked dysregulation in patients with STAT3 mutations, as observed for circIKZF2 and a newly discovered intergenic circRNA. Most circRNAs with a significantly varied CLP (96%) showed an increased expression of the circular than linear transcripts in T-LGLL compared to CTR, uncovering genes with an imbalance of circular to linear transcripts proportion in leukemic cells (Figure 1D).The ongoing validation of the differential circRNA expression in an extended cohort of 40 LGLL patients is prioritizing circRNAs for functional investigation. The integration of computational predictions and in vitro functional screening through specific circRNA silencing will attempt to understand the link between circRNAs dysregulated in T-LGLL and disease mechanisms, unveiling possible circRNA-involving axes governed by hyperactivating STAT3 somatic variants that contribute to the malignant LGL phenotypes. [Display omitted] DisclosuresNo relevant conflicts of interest to declare.

WGCNA Analysis Identifies Polycystic Ovary Syndrome-Associated Circular RNAs That Interact with RNA-Binding Proteins and Sponge miRNAs.

ObjectiveDysfunction of cumulus granulosa cells has been suggested as a contributor to abnormal folliculogenesis and the development of polycystic ovary syndrome (PCOS), but the underlying molecular mechanisms remain unclear. Recent studies indicate that circular RNAs (circRNAs) exert important roles for diseases. We aimed to screen crucial circRNAs of PCOS patients and predict their functions.MethodsThe high-throughput datasets of circRNAs (GSE145296), microRNAs (miRNAs; GSE72274) and messenger RNAs (mRNAs; GSE155489) in cumulus cells of PCOS patients and controls were collected from the Gene Expression Omnibus database. Differentially expressed circRNAs (DECs), miRNAs (DEMs) and protein-coding genes (DEGs) were identified by the limma method. The weighted correlation network analysis (WGCNA) was conducted using the DECs to mine PCOS-associated modules. Hub DECs in modules were defined as both of |gene significance| and |module membership| >0.8. The downstream effectors of hub DECs were predicted by constructing DEC-DEM-DEG ceRNA and DEC-RNA binding protein (RBP) networks. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses were performed to explore the functions of circRNAs.ResultsA total of 3614 DECs, 3544 DEGs and 1469 DEMs were identified between PCOS and controls. WGCNA analysis yielded five PCOS-related modules, of which 190 DECs were hub circRNAs. Seventeen hub DECs, nine DEMs, and 315 DEGs were identified to construct the ceRNA network, while 56 hub DECs and two DEGs (MBNL2, RBPMS) constituted the circRNA-RBP network. Five hub DECs (hsa_circ_0063309, hsa_circ_0054275, hsa_circ_0056196, hsa_circ_0018108 and hsa_circ_0070987) were overlapped between ceRNA and DEC-MBNL2 regulatory networks and thus they may be pivotal for PCOS. Furthermore, hsa_circ_0099109 could interact with the RBP gene RBPMS. Function analyses showed these circRNAs were inflammation-, apoptosis- or steroidogenesis-related.ConclusionAberrant expression of six circRNAs that function as RBP regulators or miRNA sponges may be possible mechanisms underlying the pathogenesis of PCOS by affecting apoptosis and steroidogenesis in cumulus cells.