circRNA-miRNA Network Research Articles

Systematic identification and analysis of immune-related circRNAs of Pelteobagrus fulvidraco involved in Aeromonas veronii infection

Circular RNA (circRNA) represents a type of newly discovered non-coding RNA, distinguished by its closed loop structure formed through covalent bonds. Recent studies have revealed that circRNAs have crucial influences on host anti-pathogen responses. Yellow catfish (Pelteobagrus fulvidraco), an important aquaculture fish with great economic value, is susceptible to Aeromonas veronii, a common aquatic pathogen that can cause acute death. Here, we reported the first systematic investigation of circRNAs in yellow catfish, especially those associated with A. veronii infection at different time points. A total of 1205 circRNAs were identified, which were generated from 875 parental genes. After infection, 47 circRNAs exhibited differential expression patterns (named DEcirs). The parental genes of these DEcirs were functionally engaged in immune-related processes. Accordingly, seven DEcirs (novel_circ_000226, 278, 401, 522, 736, 843, and 975) and six corresponding parental genes (ADAMTS13, HAMP1, ANG3, APOA1, FGB, and RALGPS1) associated with immunity were obtained, and their expression was confirmed by RT-qPCR. Moreover, we found that these DEcir-gene pairs likely acted through pathways, such as platelet activation, antimicrobial humoral response, and regulation of Ral protein signal transduction, to influence host immune defenses. Additionally, integrated analysis showed that, of the 7 immune-related DEcirs, three targeted 16 miRNAs, which intertwined into circRNA-miRNA networks. These findings revealed that circRNAs, by targeting genes or miRNAs are highly involved in anti-bacterial responses in yellow catfish. Our study comprehensively illustrates the roles of circRNAs in yellow catfish immune defenses. The identified DEcirs and the circRNA-miRNA network will contribute to the further investigations on the molecular mechanisms underlying yellow catfish immune responses.

Identification of differentially expressed circRNAs in prostate cancer of different clinical stages by RNA sequencing

Circular RNAs (circRNAs) are linked to cancer, but it's still not clear what role they play in prostatic cancer. Through high-throughput sequencing, the goal of this study was to compare how circRNAs are expressed at different stages of prostate cancer. 12 patients attending the Department of Urology at the Second Affiliated Hospital of Nanchang University between June 2020 and October 2021 were used for RNA sequencing, and 14 patients were used for real-time fluorescent quantitative PCR (qRT-PCR). The expression profiles of prostate cancer circRNAs were constructed by sequencing with the help of next-generation high-throughput sequencing technology, and the differentially expressed circRNAs were analyzed by targeting microRNA (miRNA) loci and Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of the genes from which circRNAs originated. Finally, the expression of target circRNAs in two prostate tissues was verified by qRT-PCR. Following high-throughput sequencing, 13,047 circRNAs were identified, and 605 circRNAs with significant differential expression were identified, of which 361 circRNAs were up-regulated, and 244 circRNAs were down-regulated. Analysis of circRNA-originated genes using GO and the KEGG enrichment analysis showed that circRNA host genes can regulate and influence multiple signaling pathways in prostate cancer with important biological functions. And the circRNA–miRNA network was constructed. The highest number of differentially expressed circRNA-binding miRNAs were: hsa_circ_000 7582 (52), hsa_circ_000 6198 (37), hsa_circ_000 6759 (28), hsa_circ_000 5675 (25), and hsa_circ_000 2172 (22). Moreover, we further screened out the circRNA (hsa_circ_0005692) that was significantly differentially expressed and common to all groups and verified by qRT-PCR that the expression of the target circRNA (hsa_circ_0005692) was significantly downregulated in prostate cancer compared with benign prostatic hyperplasia (BPH) tissues.

LncRNA and circRNA expression profiles in the hippocampus of Aβ25‑35‑induced AD mice treated with Tripterygium glycoside.

Tripterygium glycosides (TG) have been reported to ameliorate Alzheimer's disease (AD), although the mechanism involved remains to be determined. In the present study, the lncRNA and circRNA expression profiles of an AD mouse model treated with TG were assessed using microarrays. lncRNAs, mRNAs, and circRNAs in the hippocampi of 3 AD+normal saline (NS) mice and 3 AD+TG mice were detected using microarrays. The most differentially expressed lncRNAs, mRNAs, and circRNAs were screened between the AD+NS and AD+TG groups. The differentially expressed lncRNAs and circRNAs were analyzed using GO enrichment and KEGG analyses. Co-expression analysis of lncRNAs, circRNAs, and mRNAs was performed by calculating the correlation coefficients. Protein-protein interaction (PPI) network analysis was performed on mRNAs using STRING. The lncRNA-target-transcription factor (TF) network was analyzed using the Network software. In total, 661 lncRNAs, 64 circRNAs, and 503 mRNAs were found to be differentially expressed in AD mice treated with TG. Pou4f1, Egr2, Mag, and Nr4a1 were the hub genes in the PPI network. The KEGG results showed that the mRNAs that were co-expressed with lncRNAs were enriched in the TNF, PI3K-Akt, and Wnt signaling pathways. LncRNA-target-TF network analysis indicated that TFs, including Cebpa, Zic2, and Rxra, were the most likely to regulate the detected lncRNAs. The circRNA-miRNA interaction network indicated that 275 miRNAs may bind to the 64 circRNAs. In conclusion, these findings provide a novel perspective on AD pathogenesis, and the detected lncRNAs, mRNAs, and circRNAs may serve as novel therapeutic targets for the management of AD.

Circular RNA differential expression profiles and bioinformatics analysis of hsa_circRNA_079422 in human endometrial carcinoma

The aim of our study was to explore circular RNA (circRNA) expression profiles associated with human endometrial carcinoma (EC) and to analyse the molecular mechanisms involved in cancer development and their potential clinical importance. Differential expression profiles were revealed by Arraystar human circRNA microarray analysis. The results of the circRNA microarray were confirmed by quantitative real-time PCR. Interactions between circRNAs and microRNAs (miRNAs) were predicted using Arraystar’s miRNA target prediction software. The functions of the circRNA-miRNA coexpression network were identified by KEGG pathway analysis and GO analysis. Compared with para-tumorous tissues, 14 genes were significantly upregulated and 12 genes were significantly downregulated in EC tissues (P < 0.05). The quantitative real-time PCR data demonstrated consistency with the results of the microarray profile analysis. We generated a circRNA-miRNA coexpression network. Hsa_circRNA_079422 expression was significantly lower and miR-136-5p expression was higher in EC tissues than in normal endometrial tissues. KEGG pathway analysis and GO analysis indicated that hsa_circRNA_079422 might play roles in different signalling pathways and biological functions. We confirmed the presence of different circRNA expression profiles and predicted the circRNA-miRNA coexpression network in human EC tissues. Hsa_circRNA_079422 might be involved in the pathogenesis and biological process of EC via interactions with miRNAs. IMPACT STATEMENT What is already known on this subject? EC is a common malignancy of the female reproductive system. CircRNAs were demonstrated to exert critical roles in cancers, including EC. What do the results of this study add? The results of this study add circRNAs expression profiles, the circRNA-miRNA coexpression network and cancer-related circRNA-miRNA target genes in EC. It was first found that hsa_circRNA_079422 was downregulated, while miR-136-5p was upregulated in EC tissues. What are the implications of these findings for clinical practice and/or further research? In clinical practice, early EC diagnosis lacks specific biomarkers, so most EC patients are diagnosed at an advanced stage. In the management of EC patients, we also lack personalised adjuvant treatment that combines the clinical pathological characteristics. For the existing literature, we identified a new EC differential expression biomarker, hsa_circ_079422. It can be used to verify the correlation with EC clinical severity or poor prognosis. Its targeting can also be used to stratify EC patients with different molecular types, including to guide adjuvant therapy. In addition, we can verify and analyse regulatory pathways associated with it for the design of regulating engineering circRNA.

Identification and bioinformatic analysis of CircRNAs in the plasma of patients with very severe chronic obstructive pulmonary disease

BackgroundThe differential expression of circular RNAs (circRNAs) in individuals with very severe chronic obstructive pulmonary disease (COPD) and healthy individuals was screened using microarray technology. The related functions and mechanisms were analyzed using bioinformatic methods to explore the potential of target circRNAs as biomarkers of COPD and provide insights for future pathogenesis.Patients and methodsThirty patients with very severe COPD and thirty healthy controls were diagnosed at The Second People’s Hospital of Hefei from September 2021 to September 2022. The differential expression of circRNAs was compared and analyzed using a gene microarray and verified using quantitative real-time polymerase chain reaction (qRT-PCR) technology.ResultsA total of 90 upregulated and 29 downregulated circRNAs were screened in patients with very severe COPD and compared with those in healthy controls. qRT-PCR analysis showed that hsa_circ_0062683 of patients with very severe COPD was significantly upregulated, and hsa_circ_0089763 and hsa_circ_0008882 were significantly downregulated. By constructing the circRNA-miRNA interaction network, it was found that hsa-miR-612, hsa-miR-593-5p, hsa-miR-765, and hsa-miR-103a-2-5p are the miRNAs regulated by more differentially expressed circRNAs (DEcircRNAs). DEcircRNAs may participate in the development of COPD through hypoxia or regulation of various immune cells.ConclusionPlasma circRNAs may play a helpful role in the diagnosis and assessment of COPD and be valuable disease biomarkers.

A circRNA-miRNA-mRNA network analysis underlying pathogenesis of human heart failure.

The molecular mechanisms of heart failure (HF) are still poorly understood. Circular RNA (circRNA) has been discovered in the heart in increasing numbers of studies. The goal of this research is to learn more about the potential roles of circRNAs in HF. We used RNA sequencing data to identify the characteristics of circRNAs expressed in the heart and discovered that the majority of circRNAs screened were less than 2000 nt. Additionally, chromosomes One and Y had the most and least number of circRNAs, respectively. After excluding duplicate host genes and intergenic circRNAs, a total of 238 differentially expressed circRNAs (DECs) and 203 host genes were discovered. However, only four of the 203 host genes of DECs were examined in HF differentially expressed genes. Another study used Gene Oncology analysis of DECs host genes to elucidate the underlying pathogenesis of HF, and it found that binding and catalytic activity accounted for a large portion of DECs. Immune system, metabolism, and signal transduction pathways were significantly enriched. Furthermore, 1052 potentially regulated miRNAs from the top 40 DECs were collected to build a circRNA-miRNA network, and it was discovered that 470 miRNAs can be regulated by multiple circRNAs, while others are regulated by a single circRNA. In addition, a comparison of the top 10 mRNAs in HF and their targeted miRNAs revealed that DDX3Y and UTY were regulated by the most and least circRNA, respectively. These findings demonstrated circRNAs have species and tissue specific expression patterns; while circRNA expression is independent on host genes, the same types of genes in DECs and DEGs worked in HF. Our findings would contribute to a better understanding of the critical roles of circRNAs and lay the groundwork for future studies of HF molecular functions.

The diagnostic potential of a circRNA-miRNA network in non-small cell lung cancer.

Increasing studies demonstrate the significant contributions of circRNA-related competitive endogenous RNA (ceRNA) regulatory networks to tumorigenesis and cancer progression. Here, we aimed to construct a non-small cell lung cancer (NSCLC)-specific circRNA-miRNA network and evaluate its diagnostic potential in NSCLC. MiRNA deep sequencing was performed to screen differentially-expressed serum miRNAs in NSCLC. Four bioinformatics databases (TargetScan, miRanda, starBase, and RNAhybrid) were used to analyze the integrated circRNA-miRNA interaction network. The circRNA-miRNA network, including hsa-miR-4482-3p, hsa-miR-146a-3p, hsa_circ_0008167 and hsa_circ_0003317 was constructed based on their interactions and preliminary testing in NSCLC cells. The relative levels of the selected non-coding RNAs (ncRNAs) were quantified using quantitative real-time polymerase chain reaction (qRT-PCR) in the healthy, pneumonia, benign lung tumor and NSCLC cohorts. The diagnostic power of the circRNA-miRNA network was evaluated using receiver operating characteristic (ROC) analyses. The serum levels of hsa-miR-4482-3p, hsa-miR-146a-3p, hsa_circ_0008167, and hsa_circ_0003317 were dysregulated in NSCLC. The combination of the four ncRNAs showed the highest diagnostic value to discriminate between benign lung tumors and NSCLC. Additionally, the upregulated levels of hsa_circ_0008167 were correlated to more aggressive features of NSCLC, such as lymph node metastasis, distant metastasis, and higher stage. Furthermore, the combination of hsa_circ_0008167 + hsa-miR-4482-3p, and hsa_circ_0008167 + hsa-miR-4482-3p + hsa-miR-146a-3p had the greatest diagnostic power to differentiate between lymph node +/- metastases and higher/lower stages, respectively, compared to circRNAs or miRNAs alone, and traditional tumor markers. In conclusion, we identified a specific circRNA-miRNA network with higher sensitivity and specificity to diagnose NSCLC, thereby providing a new strategy for further development of ceRNA-related tumor markers in other cancers. KEY MESSAGES: Serum miR-4482-3p, miR-146a-3p, circ_0008167 and circ_0003317 are dysregulated in NSCLC. Higher levels of serum circ_0008167 are associated with more malignant NSCLC. Multiple combinations of circRNAs and miRNAs show higher value to diagnose NSCLC.

Profiling of circRNA expressions in radiation-treated head and neck cancer cells and the potential role of circPVT1

ObjectiveRadiotherapy is an indispensable treatment modality for head and neck cancers (HNCs). Due to their stable structure, circular RNAs (circRNA) have been implicated as potential biomarkers for clinical use in cancers. The purpose of this study was profiling the circRNA in radiation-treated head and neck cancer cells to identify potential differentially expressed circRNAs. DesignThe effects of radiation on the expression level of circRNAs were investigated in HNCs cells, compared to healthy cell lines. To predict the potential role of circRNAs in HNC patients, tissue expression levels, survival analyses of circRNAs, and circRNA-miRNA network were evaluated using TCGA/CPTAC datasets. Based on expression level in irradiated cells, circPVT1 (plasmacytoma variant translocation 1) was further investigated by sequence analysis. ResultsThe study revealed the characterization of differentially expressed circRNAs in cancer cells and that irradiation made significant changes in the expression of circRNAs. These findings suggest that certain circRNAs, especially circPVT1, may be potential biomarkers to monitor radiotherapy effects in patients with HNCs. ConclusionsCircRNAs may be promising molecules for improving and understanding radiotherapy efficacy in HNCs.

Key circRNAs from goat: discovery, integrated regulatory network and their putative roles in the differentiation of intramuscular adipocytes

BackgroundThe procession of preadipocytes differentiation into mature adipocytes involves multiple cellular and signal transduction pathways. Recently. a seirces of noncoding RNAs (ncRNAs), including circular RNAs (circRNAs) were proved to play important roles in regulating differentiation of adipocytes.ResultIn this study, we aimed to identificate the potential circRNAs in the early and late stages of goat intramuscular adipocytes differentiation. Using bioinformatics methods to predict their biological functions and map the circRNA-miRNA interaction network. Over 104 million clean reads in goat intramuscular preadipocytes and adipocytes were mapped, of which16 circRNAs were differentially expressed (DE-circRNAs). Furthermore, we used real-time fluorescent quantitative PCR (qRT-PCR) technology to randomly detect the expression levels of 8 circRNAs among the DE-circRNAs, and our result verifies the accuracy of the RNA-seq data. From the Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of the DE-circRNAs, two circRNAs, circ_0005870 and circ_0000946, were found in Focal adhesion and PI3K-Akt signaling pathway. Then we draw the circRNA-miRNA interaction network and obtained the miRNAs that possibly interact with circ_0005870 and circ_0000946. Using TargetScan, miRTarBase and miR-TCDS online databases, we further obtained the mRNAs that may interact with the miRNAs, and generated the final circRNA-miRNA-mRNA interaction network. Combined with the following GO (Gene Ontology) and KEGG enrichment analysis, we obtained 5 key mRNAs related to adipocyte differentiation in our interaction network, which are FOXO3(forkhead box O3), PPP2CA (protein phosphatase 2 catalytic subunit alpha), EEIF4E (eukaryotic translation initiation factor 4), CDK6 (cyclin dependent kinase 6) and ACVR1 (activin A receptor type 1).ConclusionsBy using Illumina HiSeq and online databases, we generated the final circRNA-miRNA-mRNA interaction network that have valuable functions in adipocyte differentiation. Our work serves as a valuable genomic resource for in-depth exploration of the molecular mechanism of ncRNAs interaction network regulating adipocyte differentiation.

Identification and Characterization of Circular RNAs (circRNAs) Using RNA-Seq in Two Breeds of Cashmere Goats.

Circular RNA (circRNA) is a type of non-coding RNA generated from back-splicing the reactions of linear RNA. It plays an important role in various cellular and biological processes. However, there are few studies about the regulatory effect of circRNAs on cashmere fiber traits in cashmere goats. In this study, the expression profiles of circRNAs in skin tissue were compared between Liaoning cashmere (LC) goats and Ziwuling black (ZB) goats, with a significant difference in cashmere fiber yield, cashmere fiber diameter, and cashmere fiber color, using RNA-seq. A total of 11,613 circRNAs were expressed in the caprine skin tissue, and their type, chromosomal distribution, and length distribution were characterized. A total of 115 up-regulated circRNAs and 146 down-regulated circRNAs in LC goats were screened compared to ZB goats. The authenticity of 10 differentially expressed circRNAs was validated by detecting their expression levels and the head-to-tail splice junction using RT-PCR and DNA sequencing, respectively. The parent genes of differentially expressed circRNA were mainly enriched in some Gene Ontology (GO) terms and pathways related to cashmere fiber traits, such as the canonical Wnt signaling pathway, which is involved in the regulation of cell promotion, stem cell proliferation, Wnt signaling pathway regulation, epithelial morphogenesis, MAPK signaling pathway, and cell adhesion molecules pathway. Eight differentially expressed circRNAs were further selected to construct a circRNA-miRNA network, and some miRNAs that were previously reported as related to fiber traits were found in the network. This study provides a deep understanding of the roles of circRNAs in the regulation of cashmere fiber traits in cashmere goats and the involvement of differential splicing in phenotypic expression according to breed and region.

Comprehensive analysis of dysregulated circular RNAs and construction of a ceRNA network involved in the pathology of Alzheimer's disease in a 5 × FAD mouse model.

Alzheimer's disease (AD) causes a decline in cognitive function that poses a significant hazard to human health. However, the exact pathogenesis of AD and effective treatment have both proven elusive. Circular RNAs (circRNAs), which were initially deemed as meaningless non-coding RNAs, have been shown to participate in a variety of physiological and pathological processes. However, the variations and characteristics of circRNAs are not fairly well understood during the occurrence and development of AD. In this study, we performed RNA sequencing analyses, identified circRNA expression profiles, and explored the circRNA-associated competing endogenous RNA (ceRNA) relationship in the hippocampus of five familial AD (5 × FAD) mice with cognitive dysfunction. The RNA sequencing results identified 34 dysregulated circRNAs in the hippocampus of 5 × FAD mice, including 17 upregulated and 17 downregulated circRNAs. The circRNA-miRNA interaction network for the dysregulated circRNAs was generated, and it was found to include 34 circRNAs and 711 miRNAs. Next, 2067 mRNAs potentially modulated by upregulated circRNA-interacting miRNAs and 2297 mRNAs potentially modulated by downregulated circRNA-interacting miRNAs were identified. Pathway enrichment analyses revealed that the circRNA-miRNA-mRNA network modulated AD development via multiple pathways, such as axon guidance, mitogen-activated protein kinase, and neurotrophin. The associated biological processes were mainly related to neuron projection development, cell morphogenesis, and head development. Their corresponding distributions were especially high in the axon, postsynapse, and neuronal body. We constructed a ceRNA network that included five circRNAs, four miRNAs, and 188 mRNAs. In this network, the differential expressions of three circRNAs (circRNA04655, circRNA00723, and circRNA01891), two miRNAs (miR-3470b and miR-6240), and 13 mRNAs (Vgll3, Nhsl2, Rab7, Tardbp, Vps33b, Fam107a, Tacr1, Ankrd40, Creb1, Snap23, Csnk1a1, Bmi1, and Bfar) in the hippocampus of 5 × FAD mice using qRT-PCR analyses were consistent with the RNA sequencing results. Another one circRNAs (circRNA00747) and two mRNAs (Zfp37 and Polr1e) had similar expression trends to the sequencing data, while circRNA03723 and Mapk10 had deviated expression trends to the sequencing data. In conclusion, our study uncovered dysregulated circRNA expression profiles in the hippocampus of 5 × FAD mice, stretched comprehension of ceRNA biology, investigated the potential role of this ceRNA network in pathogenesis and progression, and identified potential biomarkers and therapeutic targets for AD.

Transcriptome Sequencing Analysis of circRNA in Skeletal Muscle between Fast- and Slow-Growing Chickens at Embryonic Stages.

Skeletal muscle growth has always been the focus of the broiler industry, and circRNAs play a significant role in this process. We collected leg muscles of slow- and fast-growing Bian chicken embryos in the study at 14 (S14 and F14) and 20 (S20 and F20) days for RNA-seq. Finally, 123 and 121 differentially expressed circRNAs (DECs) were identified in S14 vs. F14 and S20 vs. F20, respectively. GO enrichment analysis for DECs obtained important biological process (BP) terms including nicotinate nucleotide biosynthetic process, nicotinate nucleotide salvage, and NAD salvage in S20 vs. F20 and protein mannosylation in S14 vs. F14. KEGG pathway analysis showed Wnt signaling pathway, Tight junction, Ubiquitin mediated proteolysis, and Notch signaling pathway were enriched in the top 20. Based on the GO and KEGG analysis results, we found some significant host genes and circRNAs such as NAPRT and novel_circ_0004547, DVL1 and novel_circ_0003578, JAK2 and novel_circ_0010289, DERA and novel_circ_0003082, etc. Further analysis found 19 co-differentially expressed circRNAs between the two comparison groups. We next constructed a circRNA-miRNA network for them, and some candidate circRNA-miRNA pairs related to skeletal muscle were obtained, such as novel_circ_0002153-miR-12219-5p, novel_circ_0003578-miR-3064-3p, and novel_circ_0010661-miR-12260-3p. These results would help to reveal the mechanism for circRNAs in skeletal muscle and also provide some guidance for the breeding of broilers.

Identification and characterization of human cytomegalovirus-encoded circular RNAs

Circular RNA (circRNA) exists extensively and plays essential roles in serving as microRNA (miRNA) or protein sponges and protein scaffolding in many organisms. However, the profiles and potential functions of the virus-encoded circRNA, including human cytomegalovirus (HCMV)-encoded circular RNAs, remain unclear. In the present study, HCMV-encoded circRNAs profile in human embryonic lung fibroblasts (HELF) with lytic infection was investigated using RNA deep sequencing and bioinformatics analysis. In total, 629 HCMV-encoded circRNAs were identified with various expression patterns in our results. The full sequences and alternative splicings of circUS12, circUL55, and circUL89 were verified by reverse transcriptase-PCR (RT-PCR) with divergent primers followed and Sanger sequencing. Transcription of circUL89 was validated by Northern blot. The HCMV-encoded circRNA-miRNA network analyses revealed the potential function of HCMV-encoded circRNAs during HCMV infection in HELFs. Collectively, HCMV infection deduced abundant HCMV-associated circRNAs during infection, and the HCMV-encoded circRNAs might play important roles in benefiting HCMV infection.

Identification and characterization of hypothalamic circular RNAs associated with bovine residual feed intake

Residual feed intake (RFI) is crucial economic indicator used for calculating the feed efficiency of growing beef cattle. circRNA plays an important biological role in gene transcriptional regulation, but little is known about its potential functional regulation underlying RFI phenotypic variation. As the core center of regulation of animal feeding, the hypothalamus is closely associated with RFI. Therefore, the present study aimed to identify the key genes and functional pathways contributing to variance in cattle RFI phenotypes using RNA sequencing from hypothalamic tissue samples, in order to gain insight into the potential regulatory role of circRNAs in bovine RFI phenotypic variation. Differentially expressed genes were detected by RNA sequencing for beef cattle in the high and low RFI groups, followed by GO, KEGG enrichment, and circRNA-miRNA co-expression network analysis. A total of 257 circRNAs were differentially expressed between the two groups, with 128 significantly upregulated and 129 significantly downregulated genes in H group compared to L group. Among them, 9 unique circRNAs were present in group L and 4 unique circRNAs were present in group H. GO and KEGG enrichment analysis of the source genes of the differentially expressed circRNAs revealed that they were mainly involved in metabolic processes, such as cellular metabolic processes, cellular macromolecular metabolic processes, and regulatory pathways related to nutrient metabolism, including protein and amino acid metabolism, as well as vitamin metabolism and pancreatic secretion associated with the animal feeding behavior. The circRNAs detected in this study were mostly novel, and have not been investigated directly to be associated with the RFI phenotype. Interestingly, most miRNAs of differentially expressed circRNAs predicted based on the circRNA-miRNA co-expression network analysis by using top 50 differentially expressed circRNAs and 13 unique circRNAs, have been reported to be related to animal RFIs, implying that circRNAs in bovine hypothalamic tissue may regulate phenotypic variation in RFI through miRNAs. The study results illustrate the complex biological functions of the hypothalamus in regulating feed efficiency and showing the potential role of circRNAs in the feeding behavior regulation of livestock, which would contributing to expanding the understanding of circRNA.

The comprehensive detection of miRNA and circRNA in the regulation of intramuscular and subcutaneous adipose tissue of Laiwu pig

circRNAs, as miRNA sponges, participate in many important biological processes. However, it remains unclear whether circRNAs can regulate lipid metabolism. This study aimed to explore the competing endogenouse RNA (ceRNA) regulatory network that affects the difference between intramuscular fat (IMF) and subcutaneous fat (SCF) deposition, and to screen key circRNAs and their regulatory genes. In this experiment, we identified 265 differentially expressed circRNAs, of which 187 up-regulated circRNA and 78 down-regulated circRNA in IMF. Subsequently, we annotated the function of DEcircRNA's host genes, and found that DEcircRNA's host genes were mainly involved in GO terms (including cellular response to fatty acids, lysophosphatidic acid acyltransferase activity, R-SMAD binding, etc.) and signaling pathways (fatty acid biosynthesis, Citrate cycle, TGF- β Signal pathway) related to adipogenesis, differentiation and lipid metabolism. By constructing a circRNA-miRNA network, we screened out DEcircRNA that can competitively bind to more miRNAs as key circRNAs (circRNA_06424 and circRNA_08840). Through the functional annotation of indirect target genes and protein network analysis, we found that circRNA_06424 affects the expression of PPARD, MMP9, UBA7 and other indirect target genes by competitively binding to miRNAs such as ssc-miR-339-5p, ssc-miR-744 and ssc-miR-328, and participates in PPAR signaling pathway, Wnt signaling pathway, unsaturated fatty acid and other signaling pathways, resulting in the difference of fat deposition between IMF and SCF. This study provide a theoretical basis for further research investigating the differences of lipid metabolism in different adipose tissues, providing potential therapeutic targets for ectopic fat deposition and lipid metabolism diseases.

Analysis of circRNA expression in chicken HD11 cells in response to avian pathogenic E.coli.

Avian pathogenic E. coli (APEC), one of the widespread zoonotic-pathogen, can cause a series of diseases collectively known as colibacillosis. This disease can cause thousands of million dollars economic loss each year in poultry industry and threaten to human health via meat or egg contamination. However, the detailed molecular mechanism underlying APEC infection is still not fully understood. Circular RNAs, a new type of endogenous noncoding RNA, have been demonstrated to involve in various biological processes. However, it is still not clear whether the circRNAs participate in host response against APEC infection. Herein, we utilized the high-throughput sequence technology to identify the circRNA expression profiles in APEC infected HD11 cells. A total of 49 differentially expressed (DE) circRNAs were detected in the comparison of APEC infected HD11 cells vs. wild type HD11 cells, which were involved in MAPK signaling pathway, Endocytosis, Focal adhesion, mTOR signaling pathway, and VEGF signaling pathway. Specifically, the source genes (BRAF, PPP3CB, BCL2L13, RAB11A, and TSC2) and their corresponding DE circRNAs may play a significant role in APEC infection. Moreover, based on ceRNA regulation, we constructed the circRNA-miRNA network and identified a couple of important regulatory relationship pairs related to APEC infection, including circRAB11A-gga-miR-125b-3p, circRAB11A-gga-miR-1696, and circTSC2-gga-miR-1649-5p. Results indicate that the aforementioned specific circRNAs and circRNA-miRNA network might have important role in regulating host immune response against APEC infection. This study is the first time to investigate the circRNAs expression profile and the biological function of the source genes of the identified DE circRNAs after APEC infection of chicken HD11 cells. These results would contribute to a better understanding of the molecular mechanisms in host response against APEC infection.

Alterations in cirRNA Expression Profile during the Development of Acute-on-Chronic Liver Failure Induced by Hepatitis B

Background. Hundreds of millions of people worldwide suffer from chronic hepatitis B (CHB), making it one of the most serious public health problems. CHB can lead to serious acute-on-chronic liver failure (ACLF). Since ACLF has a high mortality rate, predicting the risk of ACLF is a critical issue in clinical practice. Methods. To investigate the occurrence of ACLF in CHB, we used high-throughput RNA sequencing to detect the expression of a novel endogenous noncoding RNA (circRNA) in healthy controls (HC) as well as CHB and ACLF patients. Differentially expressed circRNAs were selected and analyzed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) biological pathway analyses, circRNA-miRNA coexpression networks, and statistical analyses. Results. A total of 21,101 circRNAs were identified in HC, CHB, and ACLF subjects. Compared with the HCs, 4 circRNAs were upregulated and 14 circRNAs were downregulated in ACLF subjects, with a consistent trend in all three groups. GO analysis revealed that regulation of lymphocyte activation and macrophage tolerance induction were the most important biological processes in ACLF progression. KEGG pathway analysis revealed that primary immunodeficiency and NOD-like receptor signaling pathways were the most enriched terms. The circRNA-miRNA coexpression network and statistical analyses showed that circRNA_07734, circRNA_08533, and circRNA_16083 were the most important circRNAs in the process of ACLF. Conclusions. Immune dysfunction and deficiency may play a key role in the development of ACLF, especially circRNA_07734, circRNA_08533, and circRNA_16083.

Circ_0047339 promotes the activation of fibroblasts and affects the development of urethral stricture by targeting the miR-4691-5p/TSP-1 axis

Urethral stricture is related to scar tissue fibrosis, but its pathogenesis is still unclear. This study aims to explore the regulatory mechanism of circular RNA (circRNA) in the occurrence and development of urethral stricture. CircRNA microarray was employed to analyze circRNA expression profiles between human urethral scar tissue and normal urethral tissue. The results of circRNA microarray showed that there were 296 differentially expressed genes between urethral scar tissue and normal urethral tissue. The enrichment analysis of Kyoto encyclopedia of genes and genomes showed that these circRNAs were significantly correlated with ECM–receptor interaction. The first nine differentially expressed circRNA were selected to predict the circRNA–miRNA network. RT-qPCR results showed that circ_0047339 was upregulated considerably in urethral scar tissue. Urethral scar fibroblasts were isolated from human urethral scar tissue and cultured in vitro. After silencing circ_0047339, the proliferation of urethral scar cells decreased significantly, and the expressions of Collagen I (COL-1) and α-smooth muscle actin (α-SMA) also reduced. As a competing endogenous RNA, circ_0047339 could increase the expression of TSP-1 by competitively binding miR-4691-5p. In addition, miR-4691-5p mimic transfection could inhibit the proliferation of urethral scar fibroblasts and the presentation of thrombospondin-1 (TSP-1), α-SMA and COL-1, while circ_0047339 overexpression eliminated this inhibition. Our results showed that circ_0047339 might promote the growth and fibrosis of urethral scar fibroblasts through miR-4691-5p/TSP-1 axis, thus promoting the development of urethral stricture.

Differential expression and functional analysis of circRNA in the ovaries of Yili geese at different egg-laying stages.

Ovarian development is regulated by genes that are expressed dynamically and stage-specifically. Circular RNA (circRNA) has been proven to play a significant role in the regulation of animal reproduction. Studying the expression characteristics of circRNAs in goose ovaries at various egg-laying stages can provide a reference for the molecular regulation mechanism of ovary development in geese that is mediated by circRNAs. In this study, the expression profiles of circRNAs were compared in ovary tissues from Yili geese in three different breeding periods, namely the prelaying period (KL), laying period (CL), and ceased period (XL), and differentially expressed circRNAs related to ovarian development in Yili geese were screened. The potential biological functions of differential circRNAs were predicted by bioinformatics, and the differential circRNA-miRNA regulatory network was constructed. The results showed that a total of 4483 circRNAs were identified in 12 ovarian tissue samples from Yili geese at different laying stages. In the KL vs. CL, XL vs. CL, and XL vs. KL groups, 159, 455, and 383 differentially expressed circRNAs were identified, respectively. The host genes of the differential circRNAs were mostly enriched in the signal transduction, metabolism, and other related pathways, such as those for phototransduction, glycerophospholipid metabolism, aminoacyl-tRNA biosynthesis, and retinol metabolism. Finally, we constructed circRNA-miRNA regulation networks. Nine differential circRNAs were randomly selected for qRT-PCR verification, and the expression trends were consistent with the sequencing results. Our results indicated that significant differences in the expression profiles of circRNAs in the ovaries of Yili geese at different egg-laying stages. Meanwhile, through analyzing the differential circRNA-miRNA interaction network, core regulators such as circRNA NW_013186107.1:36835|52,574 and gga-miR-34b-5p were screened. This study provides a reference for the further analysis of the molecular regulatory mechanism of the circRNAs regulating goose ovary development and enriches the theory of genetic regulation during goose ovary development.

Expression profile screening and bioinformatics analysis of CircRNA, LncRNA, and mRNA in HeLa cells infected with Chlamydia muridarum.

We have previously shown that circRNAs in host cells are involved in the process of Chlamydia trachomatis infection. In this study we aimed to identify significantly altered circRNAs/lncRNAs/mRNAs in Chlamydia muridarum infected cells and investigate their biological functions in the interaction between Chlamydia muridarum and host cells. For this purpose, circRNA, lncRNA and mRNA expression profiles were screened and identified in HeLa cells with or without Chlamydia muridarum infection by microarray. Bioinformatics analyses including Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis and Gene Ontology (GO) analysis were then carried out and the circRNA-miRNA ceRNA network was constructed. The differentially expressed circRNAs and lncRNAs were selected for validation by RT-qPCR. The results shown that a total of 834 circRNAs, 2149 lncRNAs and 1283 mRNAs were found to be differentially expressed. Enrichment analysis of GO and KEGG showed that the dysregulated genes involved nuclear-transcribed mRNA catabolic process, protein binding, RNA catabolic process and translation, the MAPK signaling pathway, apoptosis, Toll-like receptor signaling pathway, cAMP signaling pathway and Notch signaling pathway may play important roles in Chlamydia infection. Our study provides a systematic outlook on the potential function of non-coding RNAs in the molecular basis of Chlamydia infection.