CeRNA regulatory network-based analysis to study the roles of noncoding RNAs in the pathogenesis of intrahepatic cholangiocellular carcinoma

To explore and understand the competitive mechanism of ceRNAs in intrahepatic cholangiocarcinoma (ICC), we used bioinformatics analysis methods to construct an ICC-related ceRNA regulatory network (ceRNET), which contained 340 lncRNA-miRNA-mRNA regulatory relationships based on the RNA expression datasets in the NCBI GEO database. We identified the core regulatory pathway RP11-328K4.1-hsa-miR-27a-3p-PROS1, which is related to ICC, for further validation by molecular biology assays. GO analysis of 44 differentially expressed mRNAs in ceRNET revealed that they were mainly enriched in biological processes including “negative regulation of epithelial cell proliferation” and "positive regulation of activated T lymphocyte proliferation.” KEGG analysis showed that they were mainly enriched in the “complement and coagulation cascade” pathway. The molecular biology assay showed that lncRNA RP11-328K4.1 expression was significantly lower in the cancerous tissues and peripheral plasma of ICC patients than in normal controls (p<0.05). In addition, hsa-miR-27a-3p was found to be significantly upregulated in the cancer tissues and peripheral plasma of ICC patients (p<0.05). Compared to normal controls, the expression of PROS1 mRNA was significantly downregulated in ICC patient cancer tissues (p<0.05) but not in peripheral plasma (p>0.05). Furthermore, ROC analysis revealed that RP11-328K4.1, hsa-miR-27a-3p, and PROS1 had significant diagnostic value in ICC. We concluded that the upregulation of lncRNA RP11-328K4.1, which might act as a miRNA sponge, exerts an antitumor effect in ICC by eliminating the inhibition of PROS1 mRNA expression by oncogenic miRNA hsa-miR-27a.

ZEB1-AS1 initiates a miRNA-mediated ceRNA network to facilitate gastric cancer progression

Background

Currently, cancer-related competing endogenous RNA (ceRNA) networks are attracting significant interest. As long noncoding RNA ZEB1-AS1 has been reported to function as an oncogene due to sponging microRNAs (miRNAs) in several cancers, we hypothesized that it could interact with specific miRNAs to form regulatory networks and facilitate the growth of gastric cancer (GC).

Methods

MiRNAs interacting with ZEB1-AS1 were screened for and selected by bioinformatics analysis. Overexpression or repression of ZEB1-AS1 was performed to determine whether it could regulate selected miRNAs. Quantitative real-time polymerase chain reactions (qPCR) validated the expression profiles of ZEB1-AS1 and miR-149-3p in GC cell lines and tissue. Statistical analysis determined the clinical significance of ZEB1-AS1 in relation to miR-149-3p. Cell counting, wound healing and transwell assays were performed to assess cell proliferation, migration and invasion. A luciferase reporter assay was utilized to confirm the putative miR-149-3p-binding sites in ZEB1-AS1.

Results

Briefly, bioinformatics analysis inferred that ZEB1-AS1 interacts with miR-204, miR-610, and miR-149. Gain- or loss-of function assays suggested that ZEB1-AS1 negatively regulates miR-149-3p, miR-204-5p and miR-610 in GC cells. Validated by qPCR, ZEB1-AS1 was up-regulated and miR-149-3p down-regulated in GC cells and tissue. Data analyses indicated that ZEB1-AS1 and miR-149-3p are associated with the independent diagnosis and prognosis of GC. Functional assays support the theory that miR-149-3p hinders GC proliferation, migration and invasion, whereas its overexpression abrogates the corresponding effects induced by ZEB1-AS1. Lastly, dissection of the molecular mechanisms involved indicated that ZEB1-AS1 can regulate GC partly via a ZEB1-AS1/miR-149-3p axis.

Conclusions

ZEB1-AS1 can interact with specific miRNAs, forming a miRNA-mediated ceRNA network and promoting GC progress, partly through a ZEB1-AS1/miR-149-3p axis.

Electronic supplementary material

The online version of this article (10.1186/s12935-019-0742-0) contains supplementary material, which is available to authorized users.

Decoding ceRNA regulatory network in the pulmonary artery of hypoxia-induced pulmonary hypertension (HPH) rat model

Background

Hypoxia-induced pulmonary hypertension (HPH) is a lethal cardiovascular disease with the characteristic of severe remodeling of pulmonary vascular. Although a large number of dysregulated mRNAs, lncRNAs, circRNAs, and miRNAs related to HPH have been identified from extensive studies, the competitive endogenous RNA (ceRNA) regulatory network in the pulmonary artery that responds to hypoxia remains largely unknown.

Results

Transcriptomic profiles in the pulmonary arteries of HPH rats were characterized through high-throughput RNA sequencing in this study. Through relatively strict screening, a set of differentially expressed RNAs (DERNAs) including 19 DEmRNAs, 8 DElncRNAs, 19 DEcircRNAs, and 23 DEmiRNAs were identified between HPH and normal rats. The DEmRNAs were further found to be involved in cell adhesion, axon guidance, PPAR signaling pathway, and calcium signaling pathway, suggesting their crucial role in HPH. Moreover, a hypoxia-induced ceRNA regulatory network in the pulmonary arteries of HPH rats was constructed according to the ceRNA hypothesis. More specifically, the ceRNA network was composed of 10 miRNAs as hub nodes, which might be sponged by 6 circRNAs and 7 lncRNAs, and directed the expression of 18 downstream target genes that might play important role in the progression of HPH. The expression patterns of selected DERNAs in the ceRNA network were then validated to be consistent with sequencing results in another three independent batches of HPH and normal control rats. The diagnostic effectiveness of several hub mRNAs in ceRNA network was further evaluated through investigating their expression profiles in patients with pulmonary artery hypertension (PAH) recorded in the Gene Expression Omnibus (GEO) dataset GSE117261. Dysregulated POSTN, LTBP2, SPP1, and LSAMP were observed in both the pulmonary arteries of HPH rats and lung tissues of PAH patients.

Conclusions

A ceRNA regulatory network in the pulmonary arteries of HPH rats was constructed, 10 hub miRNAs and their corresponding interacting lncRNAs, circRNAs, and mRNAs were identified. The expression patterns of selected DERNAs were further validated to be consistent with the sequencing result. POSTN, LTBP2, SPP1, and LSAMP were suggested to be potential diagnostic biomarkers and therapeutic targets for PAH.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13578-022-00762-1.

Construction of Long Noncoding RNA-Associated ceRNA Networks Reveals Potential Biomarkers in Alzheimer's Disease

BACKGROUND

Alzheimer's disease (AD) is a chronic neurodegenerative disease that seriously impairs both cognitive and memory functions mainly in the elderly, and its incidence increases with age. Recent studies demonstrated that long noncoding RNAs (lncRNAs) play important roles in AD by acting as competing endogenous RNAs (ceRNAs).

OBJECTIVE

In this study, we aimed to construct lncRNA-associated ceRNA regulatory networks composed of potential biomarkers in AD based on the ceRNA hypothesis.

METHODS

A total of 20 genes (10 upregulated genes and 10 downregulated genes) were identified as the hub differentially expressed genes (DEGs). The functional enrichment analysis showed that the most significant pathways of DEGs involved include retrograde endocannabinoid signaling, synaptic vesicle circle, and AD. The upregulated hub genes were mainly enriched in the cytokine-cytokine receptor interaction pathway, whereas downregulated hub genes were involved in the neuroactive ligand-receptor interaction pathway. After convergent functional genomic (CFG) ranks and expression level analysis in different brain regions of hub genes, we found that CXCR4, GFAP, and GNG3 were significantly correlated with AD. We further identified crucial miRNAs and lncRNAs of targeted genes to construct lncRNA-associated ceRNA regulatory networks.

RESULTS

The results showed that two lncRNAs (NEAT1, MIAT), three miRNAs (hsa-miR-551a, hsa-miR-133b and hsa-miR-206), and two mRNA (CXCR4 and GNG3), which are highly related to AD, were preliminarily identified as potential AD biomarkers.

CONCLUSION

Our study provides new insights for understanding the pathogenic mechanism underlying AD, which may potentially contribute to the ceRNA mechanism in AD.

Integrated Comparative Transcriptome and circRNA-lncRNA-miRNA-mRNA ceRNA Regulatory Network Analyses Identify Molecular Mechanisms Associated with Intramuscular Fat Content in Beef Cattle

Intramuscular fat content (IMF), one of the most important carcass traits in beef cattle, is controlled by complex regulatory factors. At present, molecular mechanisms involved in regulating IMF and fat metabolism in beef cattle are not well understood. Our objective was to integrate comparative transcriptomic and competing endogenous RNA (ceRNA) network analyses to identify candidate messenger RNAs (mRNAs) and regulatory RNAs involved in molecular regulation of longissimus dorsi muscle (LDM) tissue for IMF and fat metabolism of 5 beef cattle breeds (Angus, Chinese Simmental, Luxi, Nanyang, and Shandong Black). In total, 34 circRNAs, 57 lncRNAs, 15 miRNAs, and 374 mRNAs were identified by integrating gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Furthermore, 7 key subnets with 16 circRNAs, 43 lncRNAs, 7 miRNAs, and 237 mRNAs were detected through clustering analyses, whereas GO enrichment analysis of identified RNAs revealed 48, 13, and 28 significantly enriched GO terms related to IMF in biological process, molecular function, and cellular component categories, respectively. The main metabolic-signaling pathways associated with IMF and fat metabolism that were enriched included metabolic, calcium, cGMP-PKG, thyroid hormone, and oxytocin signaling pathways. Moreover, MCU, CYB5R1, and BAG3 genes were common among the 10 comparative groups defined as important candidate marker genes for fat metabolism in beef cattle. Contributions of transcriptome profiles from various beef breeds and a competing endogenous RNA (ceRNA) regulatory network underlying phenotypic differences in IMF provided novel insights into molecular mechanisms associated with meat quality.

Identification of microRNA hsa-miR-30c-5p as an inhibitory factor in the progression of hepatocellular carcinoma and investigation of its regulatory network via comprehensive analysis

Hepatocellular carcinoma (HCC) is a primary liver cancer with high morbidity and mortality. An increasing number of abnormal gene expressions were identified to be associated with the progression of HCC. Previous studies showed that the hsa-miR-30 c-5p (miR-30 c), one of the miR-30 family members, might play a role in suppressing tumor progression in a variety of tumors. The present study aims to examine miR-30 c effects in the development of HCC. The role of miR-30 c in HCC was comprehensively investigated by using bioinformatics and experiments in vitro. The multiple databases were combined to predict and screen the target genes and upstream lncRNAs of miR-30 c, and then constructed a competitive endogenous RNA (ceRNA) regulatory network with miR-30 c as the central miRNA. The miR-30 c-related ceRNA regulatory network was also initially validated in vitro. The results showed that miR-30 c over-expression could inhibit proliferation, migration, invasion, induce apoptosis, and increase G0/G1 phase ratio of HCC cells. Three miR-30 c upstream lncRNAs and 12 miR-30 c target genes were expressed in HCC cells with increased expression and poor prognosis, and a miR-30 C-related ceRNA regulatory network was constructed. This study verified miR-30 c as an inhibitory factor in the progression of HCC and performed analyses on the miR-30 c regulatory network, which might provide potential target information for HCC prognoses and therapies. However, further experiments in vivo and studies including clinical trials will be conducted to validate our results.

lncRNA UCA1 inhibits mitochondrial dysfunction of skeletal muscle in type 2 diabetes mellitus by sequestering miR-143-3p to release FGF21

Increasing evidence suggests that insulin resistance in type 2 diabetes mellitus (T2DM) is associated with mitochondrial dysfunction in skeletal muscle, while the underlying molecular mechanisms remain elusive. This study aims to construct a ceRNA regulatory network that is involved in mitochondrial dysfunction of skeletal muscle in T2DM. Based on GEO database analysis, differentially expressed lncRNA and mRNA profiles were identified in skeletal muscle tissues of T2DM. Next, LASSO regression analysis was conducted to predict the key lncRNAs related to T2DM, which was validated by receiver operating characteristic (ROC) analysis. Moreover, the miRNAs related to skeletal muscle in T2DM were identified by WGCNA, followed by construction of gene-gene interaction network and GO and KEGG enrichment analyses. It was found that 12 lncRNAs and 6 miRNAs were related to skeletal muscle in T2DM. Moreover, the lncRNA-miRNA-mRNA ceRNA network involving UCA1, miR-143-3p, and FGF21 was constructed. UCA1, and FGF21 were downregulated, while miR-143-3p was upregulated in skeletal muscle cells (SkMCs) exposed to palmitic acid. Additionally, ectopic expression experiments were performed in SkMCs to confirm the effects of UCA1/miR-143-3p/FGF21 on mitochondrial dysfunction by determining mitochondrial ROS, oxygen consumption rate (OCR), membrane potential, and ATP level. Overexpression of miR-143-3p increased ROS accumulation and reduced the OCR, fluorescence ratio of JC-1, and ATP level, which were reversed by upregulation of UCA1 or FGF21. Collectively, lncRNA UCA1 inhibited mitochondrial dysfunction of skeletal muscle in T2DM by sequestering miR-143-3p away from FGF21, therefore providing a potential therapeutic target for alleviating mitochondrial dysfunction of skeletal muscle in T2DM.

Bioinformatics analysis of ceRNA regulatory network of baicalin in alleviating pathological joint alterations in CIA rats

Rheumatoid arthritis (RA) is an autoimmune disease characterized by chronic inffammation of synovium, leading to cartilage damage, bone erosion,even joint destruction and deformity. The conventional treatment modalities in RA are associated with side effects, emphasizing the need for alternative therapeutic remedies. Baicalin possesses multiple pharmacological effects and the advantage of low toxicity. This study aimed to reveal the potential gene regulatory mechanisms underlying the alleviating effects of baicalin in joint pathological alterations in Collagen-Induced Arthritis (CIA) rat models. At 28 days after the primary immunization, 60mg/kg/d of baicalin was administered via intraperitoneal injection once daily for 40 days, and the pathological alterations of hind paw joints were examined with X-ray imaging. Subsequently, the synovial tissue of knee joints was isolated, from which total RNA was extracted, and mRNA and miRNA sequencing libraries were established. Finally, High-throughput transcriptome sequencing(RNA-seq) technology was performed, and the lncRNAs/miRNAs/mRNAs competing endogenous RNA(ceRNA) regulatory network was analyzed. The CIA model was successfully established, and baicalin treatment significantly alleviated the destruction of distal joints of CIA rat models (p < 0.01). We found that 3 potential ceRNA regulatory networks of baicalin were established, including lncRNA ENSRNOT00000076420/miR-144-3p/Fosb, lncRNA MSTRG.1448.13/miR-144-3p/Atp2b2 and lncRNA MSTRG.1448.13/miR-144-3p/Shanks. The validation results from synovial tissue of CIA rats were consistent with the RNA-Seq results. Overall, this study revealed potentially important genes and ceRNA regulatory network that mediate the alleviating effects of baicalin on joint pathological alterations in CIA rats.

Identification and verification of m7G-Related genes as biomarkers for prognosis of sarcoma

Background: Increasing evidence indicates a crucial role for N7-methylguanosine (m7G) methylation modification in human disease development, particularly cancer, and aberrant m7G levels are closely associated with tumorigenesis and progression via regulation of the expression of multiple oncogenes and tumor suppressor genes. However, the role of m7G in sarcomas (SARC) has not been adequately evaluated. Materials and methods: Transcriptome and clinical data were gathered from the TCGA database for this study. Normal and SARC groups were compared for the expression of m7G-related genes (m7GRGs). The expression of m7GRGs was verified using real-time quantitative PCR (RT-qPCR) in SARC cell lines. Then, differentially expressed genes (DEGs) were identified between high and low m7GRGs expression groups in SARC samples, and GO enrichment and KEGG pathways were evaluated. Next, prognostic values of m7GRGs were evaluated by Cox regression analysis. Subsequently, a prognostic model was constructed using m7GRGs with good prognostic values by Lasso regression analysis. Besides, the relationships between prognostic m7GRGs and immune infiltration, clinical features, cuproptosis-related genes, and antitumor drugs were investigated in patients with SARC. Finally, a ceRNA regulatory network based on m7GRGs was constructed. Results: The expression of ten m7GRGs was higher in the SARC group than in the control group. DEGs across groups with high and low m7GRGs expression were enriched for adhesion sites and cGMP-PKG. Besides, we constructed a prognostic model that consists of EIF4A1, EIF4G3, NCBP1, and WDR4 m7GRGs for predicting the survival likelihood of sarcoma patients. And the elevated expression of these four prognostic m7GRGs was substantially associated with poor prognosis and elevated expression in SARC cell lines. Moreover, we discovered that these four m7GRGs expressions were negatively correlated with CD4⁺ T cell levels, dendritic cell level and tumor purity, and positively correlated with tumor mutational burden, microsatellite instability, drug sensitivity and cuproptosis-related genes in patients with sarcomas. Then, a triple regulatory network of mRNA, miRNA, and lncRNA was established. Conclusion: The current study identified EIF4A1, EIF4G3, NCBP1, and WDR4 as prognostic genes for SARC that are associated with m7G.These findings extend our knowledge of m7G methylation in SARC and may guide the development of innovative treatment options.

Identification of LINC02310 as an enhancer in lung adenocarcinoma and investigation of its regulatory network via comprehensive analyses

BACKGROUND

Lung adenocarcinoma (LADC) is a major subtype of non-small cell lung cancer and has one of the highest mortality rates. An increasing number of long non-coding RNAs (LncRNAs) were reported to be associated with the occurrence and progression of LADC. Thus, it is necessary and reasonable to find new prognostic biomarkers for LADC among LncRNAs.

METHODS

Differential expression analysis, survival analysis, PCR experiments and clinical feature analysis were performed to screen out the LncRNA which was significantly related to LADC. Its role in LADC was verified by CCK-8 assay and colony. Furthermore, competing endogenous RNA (ceRNA) regulatory network construction, enrichment analysis and protein-protein interaction (PPI) network construction were performed to investigate the downstream regulatory network of the selected LncRNA.

RESULTS

A total of 2431 differentially expressed LncRNAs (DELncRNAs) and 2227 differentially expressed mRNAs (DEmRNAs) were from The Cancer Genome Atlas database. Survival analysis results indicated that lnc-YARS2-5, lnc-NPR3-2 and LINC02310 were significantly related to overall survival. Their overexpression indicated poor prognostic. PCR experiments and clinical feature analysis suggested that LINC02310 was significantly correlated with TNM-stage and T-stage. CCK-8 assay and colony formation assay demonstrated that LINC02310 acted as an enhancer in LADC. In addition, 3 targeted miRNAs of LINC02310 and 414 downstream DEmRNAs were predicted. The downstream DEmRNAs were then enriched in 405 Gene Ontology terms and 11 Kyoto Encyclopedia of Genes and Genomes pathways, which revealed their potential functions and mechanisms. The PPI network showed the interactions among the downstream DEmRNAs.

CONCLUSIONS

This study verified LINC02310 as an enhancer in LADC and performed comprehensive analyses on its downstream regulatory network, which might benefit LADC prognoses and therapies.

Whole-transcriptome sequencing reveals a melanin-related ceRNA regulatory network in the breast muscle of Xichuan black-bone chicken

The economic losses incurred due to reduced muscle pigmentation highlight the crucial role of melanin-based coloration in the meat of black-bone chickens. Melanogenesis in the breast muscle of black-bone chickens is currently poorly understood in terms of molecular mechanisms. This study employed whole-transcriptome sequencing to analyze black and white breast muscle samples from black-bone chickens, leading to the identification of 367 differentially expressed (DE) mRNAs, 48 DElncRNAs, 104 DEcircRNAs, and 112 DEmiRNAs involved in melanin deposition. Based on these findings, a competitive endogenous RNA (ceRNA) network was developed to better understand the complex mechanisms of melanin deposition. Furthermore, our analysis revealed key DEmRNAs (TYR, DCT, EDNRB, MLPH and OCA2) regulated by DEmiRNAs (gga-miR-140-5p, gga-miR-1682, gga-miR-3529, gga-miR-499-3p, novel-m0012-3p, gga-miR-200b-5p, gga-miR-203a, gga-miR-6651-5p, gga-miR-7455-3p, gga-miR-31-5p, miR-140-x, miR-455-x, novel-m0065-3p, gga-miR-29b-1-5p, miR-455-y, novel-m0085-3p, and gga-miR-196-1-3p). These DEmiRNAs competitively interacted with DElncRNAs including MSTRG.2609.2, MSTRG.4185.1, LOC112530666, LOC112533366, LOC771030, LOC107054724, LOC121107411, LOC100859072, LOC101750037, LOC121108550, LOC121109224, LOC121110876, and LOC101749016, as well as DEcircRNAs, such as novel_circ_000158, novel_circ_000623, novel_001518, and novel_circ_003596. The findings from this study provide insight into the mechanisms that regulate lncRNA, circRNA, miRNA, and mRNA expression in chicken melanin deposition.

Preliminary verification of lncRNA ENST00000609755.1 potential ceRNA regulatory network in coronary heart disease

BACKGROUND

This study aims to explore the possible ceRNA regulatory network of lncRNA ENST00000609755.1 in CHD patients based on the population; reveal the possible regulatory mechanism of lncRNA ENST00000609755.1.

METHOD

Microarray analysis were used to identify differentially expressed miRNA, and mRNA profiles between 5 CHD and 5 healthy controls. The lncRNA ENST00000609755.1-miRNA-mRNA ceRNA regulatory network was constructed with lncRNA ENST00000609755.1 as the core based on microarray data and related prediction software (RNAhybird, miRanda, miRWalk 2.0). Furthermore, qRT-PCR was used to verify the expression levels of miRNA and mRNA. t-test and pearson correlation analysis were used to compare the expression differences and correlations of lncRNA, miRNA and mRNA. The receiver operating characteristic (ROC) curve was used to determine the discriminative ability of lncRNA ENST00000609755.1 and its downstream targets.

RESULTS

Totally 25 miRNAs and 953 mRNAs were differentially expressed between CHD and healthy control. The lncRNA ENST00000609755.1- miRNA- mRNA ceRNA regulatory network was constructed (5 miRNA and 58 mRNA). qRT-PCR results suggest that the expression of lncRNA ENST00000609755.1 and ELK1 were up-regulated in CHD group and positively correlated, the expression of miR-150 was down-regulated in CHD, which was negatively correlated with lncRNA ENST00000609755.1 and ELK1. The AUC was 0.777(95%CI, 0.659-0.895) when miRNA-150 and ELK1 was added, which was higher than that of lncRNA ENST00000609755.1 single indicator.

CONCLUSION

LncRNA ENST00000609755.1, miR-150 and ELK1 may have a potential ceRNA regulatory network relationship which could be considered to have a good combined diagnostic value for CHD. Also, preliminarily reveal the possible mechanism of lncRNA ENST00000609755.1 involved in CHD.

Constructing the ceRNA Regulatory Network and Combining Immune Cells to Evaluate Prognosis of Colon Cancer Patients

Objective: This study was conducted in order to construct a competitive endogenous RNA (ceRNA) network to screen RNA that plays an important role in colon cancer and to construct a model to predict the prognosis of patients.

Methods: The gene expression data of colon cancer were downloaded from the TCGA database. The difference was analyzed by the R software and the ceRNA network was constructed. The survival-related RNA was screened out by combining with clinical information, and the prognosis model was established by lasso regression. CIBERSORT was used to analyze the infiltration of immune cells in colon cancer, and the differential expression of immune cells related to survival was screened out by combining clinical information. The correlation between RNA and immune cells was analyzed by lasso regression. PCR was used to verify the expression of seven RNAs in colon cancer patients with different prognoses.

Results: Two hundred and fifteen lncRNAs, 357 miRNAs, and 2,955 mRNAs were differentially expressed in colon cancer. The constructed ceRNA network contains 18 lncRNAs, 42 miRNAs, and 168 mRNAs, of which 18 RNAs are significantly related to survival. Through lasso analysis, we selected seven optimal RNA construction models. The AUC value of the model was greater than 0.7, and there was a significant difference in the survival rate between the high- and low-risk groups. Two kinds of immune cells related to the prognosis of patients were screened out. The results showed that the expression of seven RNA markers in colon cancer patients with different prognoses was basically consistent with the model analysis.

Conclusion: We have established the regulatory network of ceRNA in colon cancer, screened out seven core RNAs and two kinds of immune cells, and constructed a comprehensive prognosis model of colon cancer patients.

Identification of Critical Modules and Biomarkers of Ulcerative Colitis by Using WGCNA

Background

Ulcerative colitis (UC) is a chronic inflammatory disease of the colon and rectum that has no exact cause and is characterized by relapsing and remitting episodes. We aimed to find biomarkers of UC and its causes.

Methods

We got GSE73661 from the GEO database and used WGCNA to find DEGs that were expressed in the same way in both normal and UC samples. To identify the co-expression modules, we used Weighted Gene Co-Expression Network Analysis. Next, we selected genes that were both DEGs and parts of main modules. Later, three datasets were used to find the hub genes, and qRT-PCR was utilized to confirm the in-silico findings. Additionally, we analyzed the connection between the hub genes and the filtration of immune cells in UC. Using the databases, we made predictions about the miRNAs and lncRNAs that regulate the hub genes and predicted possible therapeutic drugs.

Results

We found 822 DEGs and three main modules related to immunity, endoplasmic reticulum, and metabolism. Using another three datasets and human samples to confirm the mRNA expression of these genes in UC patients, XBP1 and PLPP1 were selected as hub genes, and had excellent diagnostic potential. According to the findings of the immune infiltration, patients with UC exhibited a larger proportion of immune cells. And hub genes, particularly XBP1, were closely linked to a number of immune cell infiltrations. Based on the databases and hub genes, a lncRNA-miRNA-mRNA network, including two miRNAs (miR-214-3p and miR-93-5p), two hub genes, and 124 lncRNAs, and potential therapeutic medicine were identified.

Conclusion

We found two new genes, XBP1 and PLPP1, that are involved in UC and can help diagnose and measure the disease. XBP1 also relates to clinical scores and immune cells. We suggested a gene network and possible drugs based on them.

Exploring biomarkers and therapeutic targets for pressure overload induced heart failure based on microarray data

Background

Heart failure (HF) is an end stage heart condition with poor prognosis which brings about tremendous social medical cost. Along decades, mechanism and treatments of HF have been under restless research.

Methods

In the present study, we first established pressure overload induced HF model using transaortic arch constriction (TAC) method in mice. The global expression profiles of long noncoding RNA (lncRNA), microRNA (miRNA) and messenger RNA (mRNA) were obtained by microarray probes, which were further confirmed by quantitative PCR (qPCR). Bioinformatics analysis was performed using multiple methods including volcano plotting, heatmapping and hierarchical clustering, Gene Ontology (GO) and pathway enrichment analysis, and competing endogenous RNA (ceRNA) regulatory network construction.

Results

Totally, 1,139 differentially expressed mRNAs (DEmRNAs), 3,830 lncRNAs (DElncRNAs) and 13 miRNAs (DEmiRNAs) were identified in HF group compared to control group, which could distinctly differentiate HF from normal control and were potential candidate biomarkers for HF. GO and pathway enrichment analysis revealed that multiple significant biological processes and pathways were involved in HF pathogenesis, such as extracellular matrix structural constituent, proteinaceous extracellular matrix, positive regulation of apoptotic process and integrin signaling pathway. Nine DElncRNAs, 3 DEmiRNAs and 25 DEmRNAs were filtrated out to construct a ceRNA network, which visually displayed their regulatory roles with therapeutic target potential.

Conclusions

The present study identified differentially expressed RNAs that might be involved in the pathogenesis and progression of HF. The outcomes shed lights into the underlying mechanisms for HF and provided candidate biomarkers and intervention targets for further research.

High throughput sequencing of whole transcriptome and construct of ceRNA regulatory network in RD cells infected with enterovirus D68

Background

With the advancement of sequencing technologies, a plethora of noncoding RNA (ncRNA) species have been widely discovered, including microRNAs (miRNAs), circular RNAs (circRNAs), and long ncRNAs (lncRNAs). However, the mechanism of these non-coding RNAs in diseases caused by enterovirus d68 (EV-D68) remains unclear. The goal of this research was to identify significantly altered circRNAs, lncRNAs, miRNAs, and mRNAs pathways in RD cells infected with EV-D68, analyze their target relationships, demonstrate the competing endogenous RNA (ceRNA) regulatory network, and evaluate their biological functions.

Methods

The total RNAs were sequenced by high-throughput sequencing technology, and differentially expressed genes between control and infection groups were screened using bioinformatics method. We discovered the targeting relationship between three ncRNAs and mRNA using bioinformatics methods, and then built a ceRNA regulatory network centered on miRNA. The biological functions of differentially expressed mRNAs (DEmRNAs) were discovered through GO and KEGG enrichment analysis. Create a protein interaction network (PPI) to seek for hub mRNAs and learn more about protein–protein interactions. The relative expression was verified using RT-qPCR. The effects of Fos and ARRDC3 on virus replication were confirmed using RT-qPCR, virus titer (TCID₅₀/ml), Western blotting.

Results

375 lncRNAs (154 upregulated and 221 downregulated), 33 circRNAs (32 upregulated and 1 downregulated), 96 miRNAs (49 upregulated and 47 downregulated), and 239 mRNAs (135 upregulated and 104 downregulated) were identified as differently in infected group compare to no-infected group. A single lncRNA or circRNA can be connected with numerous miRNAs, which subsequently coregulate additional mRNAs, according to the ceRNA regulatory network. The majority of DEmRNAs were shown to be connected to DNA binding, transcription regulation by RNA polymerase II, transcription factor, MAPK signaling pathways, Hippo signal pathway, and apoptosis pathway, according to GO and KEGG pathway enrichment analysis. The hub mRNAs with EGR1, Fos and Jun as the core were screened through PPI interaction network. We preliminarily demonstrated that the Fos and ARRDC3 genes can suppress EV-D68 viral replication in order to further verify the results of full transcriptome sequencing.

Conclusion

The results of whole transcriptome analysis after EV-D68 infection of RD cells were first reported in this study, and for the first time, a ceRNA regulation network containing miRNA at its center was established for the first time. The Fos and ARRDC3 genes were found to hinder viral in RD cells. This study establishes a novel insight host response during EV-D68 infection and further investigated potential drug targets.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12985-021-01686-x.

Identification and Functional Prediction of Poplar Root circRNAs Involved in Treatment With Different Forms of Nitrogen

Circular RNAs (circRNAs) are a class of noncoding RNA molecules with ring structures formed by covalent bonds and are commonly present in organisms, playing an important regulatory role in plant growth and development. However, the mechanism of circRNAs in poplar root responses to different forms of nitrogen (N) is still unclear. In this study, high-throughput sequencing was used to identify and predict the function of circRNAs in the roots of poplar exposed to three N forms [1 mM NO₃ ^- (T1), 0.5 mM NH₄NO₃ (T2, control) and 1 mM NH₄ ⁺ (T3)]. A total of 2,193 circRNAs were identified, and 37, 24 and 45 differentially expressed circRNAs (DECs) were screened in the T1-T2, T3-T2 and T1-T3 comparisons, respectively. In addition, 30 DECs could act as miRNA sponges, and several of them could bind miRNA family members that play key roles in response to different N forms, indicating their important functions in response to N and plant growth and development. Furthermore, we generated a competing endogenous RNA (ceRNA) regulatory network in poplar roots treated with three N forms. DECs could participate in responses to N in poplar roots through the ceRNA regulatory network, which mainly included N metabolism, amino acid metabolism and synthesis, response to NO₃ ^- or NH₄ ⁺ and remobilization of N. Together, these results provide new insights into the potential role of circRNAs in poplar root responses to different N forms.

An Integrated Immune-Related Bioinformatics Analysis in Glioma: Prognostic Signature's Identification and Multi-Omics Mechanisms' Exploration

As the traditional treatment for glioma, the most common central nervous system malignancy with poor prognosis, the efficacy of high-intensity surgery combined with radiotherapy and chemotherapy is not satisfactory. The development of individualized scientific treatment strategy urgently requires the guidance of signature with clinical predictive value. In this study, five prognosis-related differentially expressed immune-related genes (PR-DE-IRGs) (CCNA2, HMGB2, CASP3, APOBEC3C, and BMP2) highly associated with glioma were identified for a prognostic model through weighted gene co-expression network analysis, univariate Cox and lasso regression. Kaplan-Meier survival curves, receiver operating characteristic curves and other methods have shown that the model has good performance in predicting the glioma patients' prognosis. Further combined nomogram provided better predictive performance. The signature's guiding value in clinical treatment has also been verified by multiple analysis results. We also constructed a comprehensive competing endogenous RNA (ceRNA) regulatory network based on the protective factor BMP2 to further explore its potential role in glioma progression. Numerous immune-related biological functions and pathways were enriched in a high-risk population. Further multi-omics integrative analysis revealed a strong correlation between tumor immunosuppressive environment/IDH1 mutation and signature, suggesting that their cooperation plays an important role in glioma progression.

Long non-coding RNA crnde promotes deep vein thrombosis by sequestering miR-181a-5p away from thrombogenic Pcyox1l

BACKGROUND

Deep vein thrombosis (DVT) is an interplay of genetic and acquired risk factors, where functional interactions in lncRNA-miRNA-mRNA ceRNA networks contribute to disease pathogenesis. Based on the high-throughput transcriptome sequencing prediction, we have assessed the contribution of lncRNA Crnde/miR-181a-5p/Pcyox1l axis to thrombus formation.

METHODS

DVT was modeled in mice by inferior vena cava stenosis, and inferior vena cava tissues were harvested for high-throughput transcriptome sequencing to screen differentially expressed lncRNAs and mRNAs. The key miRNA binding to Crnde and Pcyox1l was obtained through searching the RNAInter and mirWalk databases. The binding affinity between Crnde, miR-181a-5p, and Pcyox1l was examined by FISH, dual luciferase reporter gene, RNA pull-down, and RIP assays. Functional experiments were conducted in DVT mouse models to assess thrombus formation and inflammatory injury in inferior vena cava.

RESULTS

It was noted that Crnde and Pcyox1l were upregulated in the blood of DVT mice. Crnde competitively bound to miR-181a-5p and inhibited miR-181a-5p expression, and Pcyox1l was the downstream target gene of miR-181a-5p. Silencing of Crnde or restoration of miR-181a-5p reduced inflammatory injury in the inferior vena cava, thus curtailing thrombus formation in mice. Ectopic expression of Pcyox1l counterweighed the inhibitory effect of Crnde silencing.

CONCLUSIONS

Therefore, Crnde sequesters miR-181a-5p to release Pcyox1l expression via ceRNA mechanism, thus aggravating thrombus formation in DVT.

Integrated Analysis of Transcriptome Profiles and lncRNA-miRNA-mRNA Competing Endogenous RNA Regulatory Network to Identify Biological Functional Effects of Genes and Pathways Associated with Johne's Disease in Dairy Cattle

Paratuberculosis or Johne's disease (JD), a chronic granulomatous gastroenteritis caused by Mycobacterium avium subsp. paratuberculosis (MAP), causes huge economic losses and reduces animal welfare in dairy cattle herds worldwide. At present, molecular mechanisms and biological functions involved in immune responses to MAP infection of dairy cattle are not clearly understood. Our purpose was to integrate transcriptomic profiles and competing endogenous RNA (ceRNA) network analyses to identify key messenger RNAs (mRNAs) and regulatory RNAs involved in molecular regulation of peripheral blood mononuclear cells (PBMCs) for MAP infection in dairy cattle. In total, 28 lncRNAs, 42 miRNAs, and 370 mRNAs were identified by integrating gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. In this regard, we identified 21 hub genes (CCL20, CCL5, CD40, CSF2, CXCL8, EIF2AK2, FOS, IL10, IL17A, IL1A, IL1B, IRF1, MX2, NFKB1, NFKBIA, PTGS2, SOCS3, TLR4, TNF, TNFAIP3, and VCAM1) involved in MAP infection. Furthermore, eight candidate subnets with eight lncRNAs, 29 miRNAs, and 237 mRNAs were detected through clustering analyses, whereas GO enrichment analysis of identified RNAs revealed 510, 22, and 11 significantly enriched GO terms related to MAP infection in biological process, molecular function, and cellular component categories, respectively. The main metabolic-signaling pathways related to MAP infection that were enriched included the immune system process, defense response, response to cytokine, leukocyte migration, regulation of T cell activation, defense response to bacterium, NOD-like receptor, B cell receptor, TNF, NF-kappa B, IL-17, and T cell receptor signaling pathways. Contributions of transcriptome profiles from MAP-positive and MAP-negative sample groups plus a ceRNA regulatory network underlying phenotypic differences in the intensity of pathogenicity of JD provided novel insights into molecular mechanisms associated with immune system responses to MAP infection in dairy cattle.

Characterization of the regulatory network and pathways in duodenum affecting chicken abdominal fat deposition

The excessive accumulation of abdominal fat in chickens has resulted in a reduction in both the feed conversion efficiency and the slaughter yield. To elucidate the regulatory mechanisms and metabolic pathways affecting abdominal fat deposition in the context of broiler breeding, a cohort of 400 Qingyuan partridge chickens with varying abdominal fat deposition was established. Whole transcriptome sequencing analyses were conducted on the duodenum of 20 representative chickens to ascertain the regulatory networks at this vital digestive and absorptive organ. Consequently, 116 differentially expressed genes were identified, exhibiting a trend of increasing or decreasing expression in correlation with the accumulation of abdominal fat. A total of 36 DEmRNAs, 170 DElncRNAs, 92 DEcircRNAs and 88 DEmiRNAs were identified as differentially expressed between chickens with extremely high and low abdominal fat deposition. The functional enrichment analyses demonstrated that the differentially expressed RNA in the duodenum were involved in the regulation of chicken abdominal fat deposition by mediating a series of metabolic pathways, including the Wnt signaling pathway, the PPAR signaling pathway, the Hippo signaling pathway, the FoxO signaling pathway, the MAPK signaling pathway and other signaling pathways that are involved in fatty acid metabolism and degradation. The construction of putative interaction pairs led to the suggestion of two lncRNA-miRNA-mRNA ceRNA networks comprising two mRNAs, two miRNAs, and 29 lncRNAs, as well as two circRNA-lncRNA-miRNA-mRNA ceRNA networks comprising 26 mRNAs, 12 miRNAs, 17 lncRNAs, and nine circRNAs, as core regulatory networks in the duodenum affecting chicken abdominal fat deposition. The aforementioned genes including TMEM150C, REXO1, PIK3C2G, ppp1cb, PARP12, SERPINE2, LRAT, CYP1A1, INSR and APOA4, were proposed as candidate genes, while the miRNAs, including miR-107-y, miR-22-y, miR-25-y, miR-2404-x and miR-16-x, as well as lncRNAs such as ENSGALT00000100291, TCONS_00063508, TCONS_00061201 and TCONS_00079402 were the candidate regulators associated with chicken abdominal fat deposition. The findings of this study provide a theoretical foundation for the molecular mechanisms of mRNAs and non-coding RNAs in duodenal tissues on abdominal fat deposition in chickens.

Whole-Transcriptome RNA Sequencing Uncovers the Global Expression Changes and RNA Regulatory Networks in Duck Embryonic Myogenesis

Duck meat is pivotal in providing high-quality protein for human nutrition, underscoring the importance of studying duck myogenesis. The regulatory mechanisms governing duck myogenesis involve both coding and non-coding RNAs, yet their specific expression patterns and molecular mechanisms remain elusive. To address this knowledge gap, we performed expression profiling analyses of mRNAs, lncRNAs, circRNAs, and miRNAs involved in duck myogenesis using whole-transcriptome RNA-seq. Our analysis identified 1733 differentially expressed (DE)-mRNAs, 1116 DE-lncRNAs, 54 DE-circRNAs, and 174 DE-miRNAs when comparing myoblasts and myotubes. A GO analysis highlighted the enrichment of DE molecules in the extracellular region, protein binding, and exocyst. A KEGG analysis pinpointed pathways related to ferroptosis, PPAR signaling, nitrogen metabolism, cell cycle, cardiac muscle contraction, glycerolipid metabolism, and actin cytoskeleton. A total of 51 trans-acting lncRNAs, including ENSAPLT00020002101 and ENSAPLT00020012069, were predicted to participate in regulating myoblast proliferation and differentiation. Based on the ceRNAs, we constructed lncRNA-miRNA-mRNA and circRNA-miRNA-mRNA ceRNA networks involving five miRNAs (miR-129-5p, miR-133a-5p, miR-22-3p, miR-27b-3p, and let-7b-5p) that are relevant to myogenesis. Furthermore, the GO and KEGG analyses of the DE-mRNAs within the ceRNA network underscored the significant enrichment of the glycerolipid metabolism pathway. We identified five different DE-mRNAs, specifically ENSAPLG00020001677, ENSAPLG00020002183, ENSAPLG00020005019, ENSAPLG00020010497, and ENSAPLG00020017682, as potential target genes that are crucial for myogenesis in the context of glycerolipid metabolism. These five mRNAs are integral to ceRNA networks, with miR-107_R-2 and miR-1260 emerging as key regulators. In summary, this study provides a valuable resource elucidating the intricate interplay of mRNA-lncRNA-circRNA-miRNA in duck myogenesis, shedding light on the molecular mechanisms that govern this critical biological process.

Construction of a circRNA- lincRNA-lncRNA-miRNA-mRNA ceRNA regulatory network identifies genes and pathways linked to goat fertility

Background: There is growing interest in the genetic improvement of fertility traits in female goats. With high-throughput genotyping, single-cell RNA sequencing (scRNA-seq) is a powerful tool for measuring gene expression profiles. The primary objective was to investigate comparative transcriptome profiling of granulosa cells (GCs) of high- and low-fertility goats, using scRNA-seq. Methods: Thirty samples from Ji'ning Gray goats (n = 15 for high fertility and n = 15 for low fertility) were retrieved from publicly available scRNA-seq data. Functional enrichment analysis and a literature mining approach were applied to explore modules and hub genes related to fertility. Then, interactions between types of RNAs identified were predicted, and the ceRNA regulatory network was constructed by integrating these interactions with other gene regulatory networks (GRNs). Results and discussion: Comparative transcriptomics-related analyses identified 150 differentially expressed genes (DEGs) between high- and low-fertility groups, based on the fold change (≥5 and ≤-5) and false discovery rate (FDR <0.05). Among these genes, 80 were upregulated and 70 were downregulated. In addition, 81 mRNAs, 58 circRNAs, 8 lincRNAs, 19 lncRNAs, and 55 miRNAs were identified by literature mining. Furthermore, we identified 18 hub genes (SMAD1, SMAD2, SMAD3, SMAD4, TIMP1, ERBB2, BMP15, TGFB1, MAPK3, CTNNB1, BMPR2, AMHR2, TGFBR2, BMP4, ESR1, BMPR1B, AR, and TGFB2) involved in goat fertility. Identified biological networks and modules were mainly associated with ovary signature pathways. In addition, KEGG enrichment analysis identified regulating pluripotency of stem cells, cytokine-cytokine receptor interactions, ovarian steroidogenesis, oocyte meiosis, progesterone-mediated oocyte maturation, parathyroid and growth hormone synthesis, cortisol synthesis and secretion, and signaling pathways for prolactin, TGF-beta, Hippo, MAPK, PI3K-Akt, and FoxO. Functional annotation of identified DEGs implicated important biological pathways. These findings provided insights into the genetic basis of fertility in female goats and are an impetus to elucidate molecular ceRNA regulatory networks and functions of DEGs underlying ovarian follicular development.

Integrative analysis of m7G methylation-associated genes prognostic signature with immunotherapy and identification of LARP1 as a key oncogene in head and neck squamous cell carcinoma

Background

N7-methylguanosine (m7G) methylation is an RNA modification associated with cancer progression, but its specific role in head and neck squamous cell carcinoma (HNSCC) remains unclear.

Methods

This study analyzed the differential expression of m7G-related genes (m7GRGs) in HNSCC using the TCGA-HNSCC dataset, identifying key pathways associated with the cell cycle, DNA replication, and focal adhesion. A LASSO-Cox regression model was constructed based on four m7GRGs (EIF3D, EIF1, LARP1, and METTL1) and validated with GEO datasets and clinical samples. Further validation of gene upregulation in HNSCC tissues was conducted using RT-qPCR and immunohistochemistry, while the role of LARP1 in HNSCC cells was assessed via knockout experiments.

Results

The constructed model demonstrated strong predictive performance, with the risk score significantly correlating with prognosis, immune infiltration, and drug sensitivity. An external dataset and clinical specimens further confirmed the model's predictive accuracy for immunotherapy response. Additionally, two regulatory axes-LINC00707/hsa-miR-30b-5p/LARP1 and SNHG16/hsa-miR-30b-5p/LARP1-were identified. LARP1 knockout experiments revealed that suppressing LARP1 markedly inhibited HNSCC cell proliferation, migration, and invasion.

Conclusion

The m7GRG-based prognostic model developed in this study holds strong clinical potential for predicting prognosis and therapeutic responses in HNSCC. The identification of LARP1 and its related regulatory pathways offers new avenues for targeted therapy in HNSCC.