ceRNA Regulatory Network Research Articles

Identification of CFHR4 as a Potential Prognosis Biomarker Associated With lmmune Infiltrates in Hepatocellular Carcinoma.

BackgroundComplement factor H-related 4 (CFHR4) is a protein-coding gene that plays an essential role in multiple diseases. However, the prognostic value of CFHR4 in hepatocellular carcinoma (HCC) is unknown.MethodsUsing multiple databases, we investigated CFHR4 expression levels in HCC and multiple cancers. The relationship between CFHR4 expression levels and clinicopathological variables was further analyzed. Various potential biological functions and regulatory pathways of CFHR4 in HCC were identified by performing a Gene Ontology (GO) analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis and Gene Set Enrichment Analysis (GSEA). Single-sample gene set enrichment analysis (ssGSEA) was performed to confirm the correlation between CFHR4 expression and immune cell infiltration. The correlations between CFHR4 expression levels in HCC and N6-methyladenosine (m6A) modifications and the competing endogenous RNA (ceRNA) regulatory networks were confirmed in TCGA cohort.ResultsCFHR4 expression levels were significantly decreased in HCC tissues. Low CFHR4 expression in HCC tissues was significantly correlated with the patients’ sex, race, age, TNM stage, pathological stage, tumor status, residual tumor, histologic grade and alpha fetal protein (AFP) level. GO and KEGG analyses revealed that differentially expressed genes related to CFHR4 may be involved in the synaptic membrane, transmembrane transporter complex, gated channel activity, chemical carcinogenesis, retinol metabolism, calcium signaling pathway, PPAR signaling pathway, insulin and gastric acid secretion. GSEA revealed that the FCGR-activated reaction, PLK1 pathway, ATR pathway, MCM pathway, cascade reactions of PI3K and FGFR1, reactant-mediated MAPK activation and FOXM1 pathway were significantly enriched in HCC with low CFHR4 expression. Moreover, CFHR4 expression was inversely correlated the levels of infiltrating Th2 cells, NK CD56bright cells and Tfh cells. In contrast, we observed positive correlations with the levels of infiltrating DCs, neutrophils, Th17 cells and mast cells. CFHR4 expression showed a strong correlation with various immunomarker groups in HCC. In addition, high CFHR4 expression significantly prolonged the overall survival (OS), disease-specific survival (DSS) and progression-free interval (PFI). We observed a substantial correlation between the expression of CFHR4 and multiple N6-methyladenosine genes in HCC and constructed potential CFHR4-related ceRNA regulatory networks.ConclusionsCFHR4 might be a potential therapeutic target for improving the HCC prognosis and is closely related to immune cell infiltration.

Identification and Validation of Immune-Related Biomarker Gene and Construction of ceRNA Networks in Septic Cardiomyopathy.

Septic cardiomyopathy (SCM) is a cardiac dysfunction caused by severe sepsis, which greatly increases the risk of heart failure and death, and its molecular mechanism is unclear. The immune response has been reported to be an important process in septic cardiomyopathy and is present in the cardiac tissue of patients with sepsis, suggesting that the immune response may be an underlying mechanism of myocardial injury in SCM. Therefore, we explored the role of immune-related genes (IRGs) in SCM and aimed to identify pivotal immune-related targets with the aim of identifying key immune-related targets in SCM and potential therapeutic mechanisms involved in the pathological process of SCM. To explore the regulatory mechanisms of immune responses in SCM, we identified differentially expressed genes (DEGs) shared in the SCM datasets GSE179554 and GSE40180 by bioinformatics analysis and then obtained hub genes from the DEGs. Then, we obtained the immune-related hub genes (IRHGs) by intersecting the hub genes with IRGs and performed quantitative reverse transcription polymerase chain reaction to confirm the abnormal expression of IRHGs. Finally, we further constructed an immune-related lncRNA–miRNA–IRHG ceRNA regulatory network. In this study, we identified an IRHG that may be involved in the pathogenesis of SCM, which helps us to further elucidate the role of immune response in SCM and gain insights into the molecular mechanisms and potential therapeutic targets of SCM.

Integrated Analysis of Immune-Related circRNA-miRNA-mRNA Regulatory Network in Ischemic Stroke.

BackgroundStroke is the leading cause of death and disability worldwide, with ischemic stroke (IS) being the most prevalent type. Circular RNAs (circRNAs) are involved in the pathological process of IS and are promising biomarkers for the diagnosis of IS. However, studies focusing on circRNAs acting as microRNAs (miRNAs) sponges in regulating mRNA expression are currently scarce.MethodsIn this study, expression profiles of circRNAs (GSE195442), miRNAs (GSE117064), and mRNAs (GSE58294) from the Gene Expression Omnibus (GEO) database were analyzed. Differentially expressed circRNAs (DEcircRNAs), differentially expressed miRNAs (DEmiRNAs), and differentially expressed mRNAs (DEmRNAs) were identified by R software. The target miRNAs and target genes were predicted by several bioinformatics methods. Then, we performed Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of the DEmRNAs. Subsequently, the protein-protein interaction (PPI) network and the competing endogenous RNA (ceRNA) regulatory network were visualized by Cytoscape software. Finally, we further constructed an immune-related circRNA-miRNA-mRNA regulatory sub-network in IS.ResultsA total of 35 DEcircRNAs, 141 DEmiRNAs, and 356 DEmRNAs were identified. By comprehensive analysis of bioinformatics methods, we constructed a circRNA-miRNA-mRNA regulatory network, including 15 DEcircRNAs, eight DEmiRNAs, and 39 DEmRNAs. FGF9 was identified as an immune-related hub gene. Immune cell analysis indicated a significantly higher level of neutrophils in IS, and the expression of FGF9 was significantly negatively correlated with the level of neutrophils. Eventually, miR-767-5p was predicted as the upstream molecules of FGF9, and circ_0127785 and circ_0075008 were predicted as the upstream circRNAs of miR-767-5p.ConclusionOur study provides novel insights into the molecular mechanisms governing the progression of IS from the perspective of immune-related ceRNA networks.

Transcriptomics analysis reveals key lncRNAs and genes related to the infection of porcine lung macrophages by Glaesserella parasuis

Glaesserella parasuis (G. parasuis) is the pathogen of Glässer's disease in pig herds, which can cause severe inflammatory responses. However, at present, the pathogenic mechanism of G. parasuis is not very clear. LncRNAs can regulate the expression of mRNA in a variety of ways, thereby causing host cells to produce a variety of functional changes in response to bacterial infection. Here, we detected the changes in lncRNAs and mRNAs of 3D4/21 cells after G. parasuis CY1201 strain (serotype 13) infection. A total of 876 lncRNAs and 2166 mRNAs were differentially expression in 3D4/21 cells after G. parasuis infection. GO and KEGG enrichment analysis showed that the differentially up-regulated lncRNA target genes were mainly involved in the response to extracellular stimuli, cell receptor signaling pathways and chemokine signaling pathways. The differentially down-regulated lncRNA target genes were mainly involved in ERK1/ERK2 cascade reaction and adhesion junctions. 44 lncRNAs were screened that might be related in inflammation. CeRNA regulatory network of the top five difference inflammation-related lncRNAs showed that the up-regulated lncRNA group involved 5 lncRNAs, 50 miRNAs and 49 mRNAs. Meanwhile, there were 26 miRNAs and 36 mRNAs in the top five down-regulated lncRNA group. Our results contribute to understand the basic role of lncRNAs in 3D4/21 cells during G. parasuis infection, and lay the foundation for following research.

Construction and Investigation of circRNA-associated ceRNA Regulatory Network in Molecular Subtypes of Breast Cancer.

Circular RNAs (circRNAs) act as competing endogenous RNAs (ceRNAs) that indirectly regulate gene expression and function by binding microRNAs (miRNAs). A growing body of evidence indicates that the ceRNA networks can be used as an effective method to investigate cancer; however, the construction and analysis of ceRNA networks, especially circRNA-miRNA-mRNA regulatory network, in different subtypes of breast cancer have not been previously performed. In the present study, we generated a ceRNA network to explore their roles in two BC subtypes, namely Luminal A and triple negative breast cancer (TNBC). First, the expression profiles of circRNA, miRNA, and mRNA were downloaded from the GEO database, differentially expressed genes were obtained using GEO2R, and a ceRNA network was constructed based on circRNA-miRNA pairs and miRNA-mRNA pairs, consisted of 10 circRNAs, 25 miRNAs and 39 mRNAs. Further studies of BC subtypes based on TCGA datasets were also performed to validate the effect of a novel ceRNA network. Then, the related genes in the regulatory network were analyzed by GO functional annotation and KEGG pathway enrichment. The analysis showed that targeted genes were enriched in 97 GO terms and 25 KEGG pathways, involved in the molecular typing of breast cancer. Meanwhile, Kaplan-Meier analysis revealed that three key genes (MKI67, DEF8, and GFRA1) were significantly associated with BC tumor differentiation and prognosis. The current study provides a potential application of the ceRNA network within BC subtypes and may offer new targets for their diagnosis, therapy and prognosis.

Comprehensive Expression Profile Analysis of Neutrophil Extracellular Trap-Affected Genes in Gastric Cancer Cells and the Clinical Significance of lncRNA NEAT1-Related Signaling.

BackgroundGastric cancer (GC) is the fifth most common malignant tumor and the third leading cause of cancer-related deaths worldwide. Neutrophil extracellular traps (NETs) can enhance the invasion of GC cells and are associated with poor prognosis in patients. However, its mechanism of action is not completely understood.MethodsThe content of NETs in the peripheral blood of patients with GC was detected by enzyme-linked immunosorbent assay. GC AGS cells were treated with or without NETs for 24 h. High-throughput RNA sequencing was performed to screen differentially expressed long non-coding RNAs (lncRNAs), microRNAs (miRNAs), and messenger RNAs (mRNAs). Real-time polymerase chain reaction (PCR) was used to verify gene expression. A competing endogenous RNA (ceRNA) regulatory network was constructed. Modules were screened using the molecular complex detection (MCODE) plug-in. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses were performed using the genes in the network. The role and clinical significance of the lncRNA NEAT1-related signaling pathway were validated.ResultsThe content of NETs in the patients with GC was significantly higher than that in healthy controls and was also higher in patients with high-grade (stages III and IV) GC. NETs promoted the invasion of AGS cells. A total of 1,340 lncRNAs, 315 miRNAs, and 1,083 mRNAs were differentially expressed after NET treatment. The expression of five genes was validated using real-time PCR, which were in accordance with the RNA sequencing results. A ceRNA regulatory network was constructed with 1,239 lncRNAs, 310 miRNAs, and 1,009 mRNAs. Four genes (RAB3B, EPB41L4B, ABCB11, and CCDC88A) in the ceRNA network were associated with patient prognosis, with RAB3B being the most prominent and with signaling among the lncRNA NEAT1, the miRNA miR-3158-5p, and RAB3B. NEAT1 was upregulated in AGS cells after NET treatment. RNA interference of NEAT1 inhibited the invasion of AGS cells induced by NETs, inhibited miR-3158-5p expression, and promoted RAB3B expression. NEAT1 and RAB3B expression were positively correlated in patients with GC. Furthermore, RAB3B was upregulated and miR-3158-5p was downregulated in GC tissues compared with adjacent normal tissues, which was also associated with cancer stage.ConclusionThis study provides a comprehensive analysis of differentially expressed genes in NET-treated GC cells and validated the clinical significance of NEAT1-related signaling.

Integrated Analysis of lncRNA and circRNA Mediated ceRNA Regulatory Networks in Skin Reveals Innate Immunity Differences Between Wild-Type and Yellow Mutant Rainbow Trout (Oncorhynchus mykiss).

Fish skin is a vital immune organ that forms the first protective barrier preventing entry of external pathogens. Rainbow trout is an important aquaculture fish species that is farmed worldwide. However, our knowledge of innate immunity differences between wild-type (WR_S) and yellow mutant rainbow trout (YR_S) remains limited. In this study, we performed whole transcriptome analysis of skin from WR_S and YR_S cultured in a natural flowing water pond. A total of 2448 mRNAs, 1630 lncRNAs, 22 circRNAs and 50 miRNAs were found to be differentially expressed (DE). Among these DEmRNAs, numerous key immune-related genes, including ifih1, dhx58, trim25, atp6v1e1, tap1, tap2, cd209, hsp90a.1, nlrp3, nlrc3, and several other genes associated with metabolism (gstp1, nampt, naprt and cd38) were identified. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of DEmRNAs revealed that many were significantly enriched in innate immune-related GO terms and pathways, including NAD+ADP-ribosyltransferase activity, complement binding, immune response and response to bacterium GO terms, and RIG-I-like receptor signaling, NOD-like receptor signaling and phagosome KEGG pathways. Furthermore, the immune-related competing endogenous RNA networks were constructed, from which we found that lncRNAs MSTRG.11484.2, MSTRG.32014.1 and MSTRG.29012.1 regulated at least three immune-related genes (ifih1, dhx58 and irf3) through PC-5p-43254_34, PC-3p-28352_70 and bta-miR-11987_L-1R-1_1ss8TA, and tap2 was regulated by two circRNAs (circRNA5279 and circRNA5277) by oni-mir-124a-2-p5_1ss13GA. The findings expand our understanding of the innate immune system of rainbow trout, and lay the foundation for further study of immune mechanisms and disease resistance breeding.

Analysis of the lncRNA-Associated Competing Endogenous RNA (ceRNA) Network for Tendinopathy.

Background We aimed to construct the lncRNA-associated competing endogenous RNA (ceRNA) network and distinguish feature lncRNAs associated with tendinopathy. Methods We downloaded the gene profile of GSE26051 from the Gene Expression Omnibus (GEO), including 23 normal samples and 23 diseased tendons. Differentially expressed mRNAs (DEmRNAs) and differentially expressed lncRNAs (DElncRNAs) were identified, and functional and pathway enrichment analyses were performed. Protein-protein interaction (PPI) network was constructed and further analyzed by module mining. Moreover, a ceRNA regulatory network was constructed based on the identified lncRNA–mRNA coexpression relationship pairs and miRNA–mRNA regulation pairs. Results We identified 1117 DEmRNAs and 57 DElncRNAs from the GEO data. The downregulated DEmRNAs were particularly associated with muscle contraction and muscle filament, while the upregulated ones were linked to extracellular matrix organization and cell adhesion. From the PPI network, 11 modules were extracted. Genes in MCODE 2 (such as TPM4) were significantly involved in cardiomyopathy, and genes in MCODE 4 (such as COL4A3 and COL4A4) were involved in focal adhesion, ECM-receptor interaction, and PI3K-Akt signaling pathway. The ceRNA network contained 7 lncRNAs (MIR133A1HG, LINC01405, PRKCQ-AS1, C10orf71-AS1, MBNL1-AS1, HOTAIRM1, and DNM3OS), 63 mRNAs, and 41 miRNAs. Downregulated lncRNA MIR133A1HG could competitively bind with hsa-miR-659-3p and hsa-miR-218-1-3p to regulate the TPM3. Meanwhile, MIR133A1HG could competitively bind with hsa-miR-1179 to regulate the COL4A3. Downregulated C10orf71-AS1 could competitively bind with hsa-miR-130a-5p to regulate the COL4A4. Conclusions Seven important lncRNAs, particularly MIR133A1HG and C10orf71-AS1, were found associated with tendinopathy according to the lncRNA-associated ceRNA network.

MiR-29a-5p Inhibits Prenatal Hair Placode Formation Through Targeting EDAR by ceRNA Regulatory Network.

Hair placode formation is an important stage of hair follicle morphogenesis and it is a complex process facilitated by non-coding RNAs. In this study, we conducted whole transcriptome sequencing analysis of skin, heart, liver, lung, and kidney tissues of day 41 (E41) normal and hairless pig embryos, and respectively detected 15, 8, and 515 skin-specific differentially expressed (DE) lncRNAs, miRNAs, and mRNAs. Furthermore, 18 competing endogenous RNA (ceRNA) networks were constructed. Following weighted gene co-expression network analysis (WGCNA) of stages E39, E41, E45, E52, and E60, between normal and hairless pig embryos, only two ceRNAs (lncRNA2162.1/miR-29a-5p/BMPR1b and lncRNA627.1/miR-29a-5p/EDAR) that showed period-specific differential expression in E41 skin were retained. Dual-luciferase reporter assays further indicated that EDAR was a direct, functioning target of miR-29a-5p and that no binding site was found in BMPR1b. Moreover, miR-29a-5p overexpression inhibited the mRNA and protein expression of EDAR while no significant differential expression of BMPR1b was detected. In addition, over-expressed lncRNA627.1 reduces the expression of miR-29a-5p and increase EDAR expression while inhibits lncRNA627.1 resulted in a opposite expression trend. Cell proliferation result demonstrated that lower expression of EDAR and lncRNA627.1 inhibited hair placode precursor cells (HPPCs) proliferation in a manner similar to that shown by over-expressed miR-29a-5p. This study identified that miR-29a-5p inhibited HPPCs proliferation via the suppression of EDAR expression in the EDA/EDAR signaling pathway, while lncRNA627.1 rescues EDAR expression. Our study provides a basis for a better understanding of the mechanisms underlying the ceRNA complex, miR29a-5p/EDAR/lncRNA627.1, that could regulate hair placode formation, which may help decipher diseases affecting human hair.

Genome-Wide Identification and Characterization of Long Noncoding RNAs in Populus × canescens Roots Treated With Different Nitrogen Fertilizers.

Nitrate (NO3−) and ammonium (NH4+) are the primary forms of inorganic nitrogen acquired by plant roots. LncRNAs, as key regulators of gene expression, are a class of non-coding RNAs larger than 200 bp. However, knowledge about the regulatory role of lncRNAs in response to different nitrogen forms remains limited, particularly in woody plants. Here, we performed strand-specific RNA-sequencing of P. × canescens roots under three different nitrogen fertilization treatments. In total, 324 lncRNAs and 6,112 mRNAs were identified as showing significantly differential expression between the NO3− and NH4NO3 treatments. Moreover, 333 lncRNAs and 6,007 mRNAs showed significantly differential expression between the NH4+ and NH4NO3 treatments. Further analysis suggested that these lncRNAs and mRNAs have different response mechanisms for different nitrogen forms. In addition, functional annotation of cis and trans target mRNAs of differentially expressed lncRNAs indicated that 60 lncRNAs corresponding to 49 differentially expressed cis and trans target mRNAs were involved in plant nitrogen metabolism and amino acid biosynthesis and metabolism. Furthermore, 42 lncRNAs were identified as putative precursors of 63 miRNAs, and 28 differentially expressed lncRNAs were potential endogenous target mimics targeted by 96 miRNAs. Moreover, ceRNA regulation networks were constructed. MSTRG.6097.1, MSTRG.13550.1, MSTRG.2693.1, and MSTRG.12899.1, as hub lncRNAs in the ceRNA networks, are potential candidate lncRNAs for studying the regulatory mechanism in poplar roots under different nitrogen fertilization treatments. The results provide a basis for obtaining insight into the molecular mechanisms of lncRNA responses to different nitrogen forms in woody plants.

The LINC01929/miR-6875-5p/ADAMTS12 Axis in the ceRNA Network Regulates the Development of Advanced Bladder Cancer.

Considering its speedy development and extremely low 5-year overall survival rate worldwide, bladder cancer (BCa) is one of the most common and highly malignant tumors. Increasing evidence suggests that protein-coding mRNAs and non-coding RNAs, including long non-coding RNAs (lncRNAs) and micro RNAs (miRNAs), play an essential role in regulating the biological processes of cancer. To investigate the molecular regulation associated with poor prognosis during advanced BCa development, we constructed a competitive endogenous RNA (ceRNA) network. Using transcriptome profiles from The Cancer Genome Atlas and Gene Expression Omnibus databases, we performed differential expression (DE) analysis, weighted gene co-expression network analysis, functional enrichment analysis, survival analysis, prediction of miRNA targeting, and Pearson correlation analysis. Through layers of selection, 8 lncRNAs-28 mRNAs and 8 miRNAs-28 mRNAs pairs shared similar expression patterns, constituting a core ceRNA regulatory network related to the invasion, progression, and metastasis of advanced clinical stage (ACS) BCa. Subsequently, we conducted real time qPCR, western blotting, and immunohistochemistry to validate expression trend bioinformatics analysis on 3, 2, and 3 differentially expressed mRNAs, lncRNAs, and miRNAs, respectively. The most significantly differentially expressed LINC01929, miR-6875-5p and ADAMTS12 were selected for in vitro experiments to assess the functional role of the LINC01929/miR-6875-5p/ADAMTS12 axis. RNA pull-down, luciferase assays, and rescue assays were performed to examine the binding of LINC01929 and miR-6875-5p. Increasing trends in COL6A1, CDH11, ADAMTS12, LINC01705, and LINC01929 expression variation were verified as consistent with previous DE analysis results in ACS-BCa, compared with low clinical stage (LCS) BCa. Expression trends in parts of these RNAs, such as hsa-miR-6875-5p, hsa-miR-6784-5p, COL6A1, and CDH11, were measured in accordance with DE analysis in LCS-BCa, compared with normal bladder urothelium. Through experimental validation, the cancer-promoting molecule ADAMST12 was found to play a key role in the development of advanced BCa. Functionally, ADAMTS12 knockdown inhibited the progression of bladder cancer. Overexpression of LINC01929 promoted bladder cancer development, while overexpression of miR-6785-5p inhibited bladder cancer development. Mechanistically, LINC01929 acted as a sponge for miR-6785-5p and partially reversed the role of miR-6785-5p. Our findings provide an elucidation of the molecular mechanism by which advanced bladder cancer highly expressed LINC01929 upregulates ADAMTS12 expression through competitive adsorption of miR-6875-5p. It provides a new target for the prognosis and diagnosis of advanced bladder cancer.

The ceRNA Crosstalk between mRNAs and lncRNAs in Diabetes Myocardial Infarction

Competitive endogenous RNA regulation suggests an intricate network of all transcriptional RNAs that have the function of repressing miRNA function and regulating mRNA expression. Today, the specific ceRNA regulatory mechanisms of lncRNA–miRNA–mRNA in patients who have diabetes mellitus (DM) and myocardial infarction (MI) are still unknown. Two data sets, GSE34198 and GSE112690, were rooted in the Gene Expression Omnibus database to search for changes of lncRNA, miRNA, and mRNA in MI patients with diabetes. Weighted gene correlation network analysis (WGCNA) was used to identify the modules related to the development of diabetes in patients with MI. Target genes of miRNAs were predicted using miRWalk, TargetScan, mirDB, RNA22, and miRanda. Then, functional and enrichment analyses were performed to build the lncRNA–miRNA–mRNA interaction network. We built ceRNA regulatory networks with three lncRNAs, two miRNAs, and nine mRNAs. Differentially expressed genes enriched in biological process, including neutrophil activation, refer to immune response and positive system of defense feedback. Besides, there is significant enrichment in molecular function of calcium toll−like receptor binding, icosanoid binding, RAGE receptor binding, and arachidonic acid binding. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis enriched differentially expressed genes (DEGs) in pathways that were well known in MI, indicating inflammation and immune response. Pathways associated with diabetes were also significantly enriched. We confirmed significantly altered lncRNA, miRNA, and mRNA in MI patients with diabetes, which might serve as biomarkers for the progress and development of diabetic cardiovascular diseases. We constructed a ceRNA regulatory network of lncRNA–miRNA–mRNA, which will enable us to understand the novel molecular mechanisms included in the initiation, progression, and interaction between DM and MI, laying the foundation for clinical diagnosis and treatment.

Identification of Angiogenesis-Related Prognostic Biomarkers Associated With Immune Cell Infiltration in Breast Cancer.

Background: This study aimed to explore the prognostic value of angiogenesis-related genes (ARGs) and their association with immune cell infiltration (ICI) in breast cancer (BC). Methods: Transcriptome data of BC were obtained from the TCGA and GEO databases. Differentially expressed ARGs were identified by the limma package. The identification of key genes and construction of the risk score model were performed by univariate and multivariate Cox regression algorithms. The prognostic value of the risk score was assessed by ROC curves and nomogram. GO, KEGG pathway, and GSEA were used to investigate the biological functions of differentially expressed genes (DEGs), and CIBERSORT, ssGSEA, and xCell algorithms were performed to estimate the ICI in high-risk and low-risk groups. The correlations between prognostic biomarkers and differentially distributed immune cells were assessed. Moreover, a ceRNA regulatory network based on prognostic biomarkers was constructed and visualized by Cytoscape software. Results: A total of 18 differentially expressed ARGs were identified between tumor and adjacent normal tissue samples. TNFSF12, SCG2, COL4A3, and TNNI3 were identified as key prognostic genes by univariate and multivariate Cox regression analyses. The risk score model was further constructed based on the four-gene signature and validated in GSE7390 and GSE88770 datasets. ROC curves and nomogram indicated that the risk score had good accuracy for determining BC patient survival. Biological function analysis showed that DEGs in high- and low-risk groups had a high enrichment in immune-related biological processes and signaling pathways. Moreover, significantly different ICIs were found between high- and low-risk groups, such as memory B cells, CD8+ T cells, resting memory CD4+ T cells, follicular helper T cells, regulatory T cells, monocytes, M2 macrophages, and neutrophils, and each prognostic biomarker was significantly correlated with one or more immune cell types. Conclusion: The current study identified novel prognostic ARGs and developed a prognostic model for predicting survival in patients with BC. Furthermore, this study indicated that ICI may act as a bond between angiogenesis and BC. These findings enhance our understanding of angiogenesis in BC and provide novel guidance on developing therapeutic targets for BC patients.

Construction of lncRNA-miRNA-mRNA network based on ceRNA mechanism reveals the function of lncRNA in the pathogenesis of gout.

ObjectiveTo identify differentially expressed lncRNA, miRNA, and mRNA during the pathogenesis of gout, explore the ceRNA network regulatory mechanism of gout, and seek potential therapeutic targets.MethodFirst, gout‐related chips were retrieved by GEO database. Then, the analysis of differentially expressed lncRNAs and mRNAs was conducted by R language and other software. Besides, miRNA and its regulated mRNA were predicted based on public databases, the intersection of differentially expressed mRNA and predicated mRNA was taken, and the lncRNA‐miRNA‐mRNA regulatory relationships were obtained to construct the ceRNA regulatory network. Subsequently, hub genes were screened by the STRING database and Cytoscape software. Then the DAVID database was used to illustrate the gene functions and related pathways of hub genes and to mine key ceRNA networks.ResultsThree hundred and eighty‐eight lncRNAs and 758 mRNAs were identified with significant differential expression in gout patient, which regulates hub genes in the ceRNA network, such as JUN, FOS, PTGS2, NR4A2, and TNFAIP3. In the ceRNA network, lncRNA competes with mRNA for miRNA, thus affecting the IL‐17 signaling pathway, TNF signaling pathway, Oxytocin signaling pathway, and NF‐κB signaling pathway through regulating the cell's response to chemical stress. The research indicates that five miRNAs (miR‐429, miR‐137, miR‐139‐5p, miR‐217, miR‐23b‐3p) and five lncRNAs (SNHG1, FAM182A, SPAG5‐AS1, HNF1A‐AS1, UCA1) play an important role in the formation and development of gout.ConclusionThe interaction in the ceRNA network can affect the formation and development of gout by regulating the body's inflammatory response as well as proliferation, differentiation, and apoptosis of chondrocytes and osteoclasts. The identification of potential therapeutic targets and signaling pathways through ceRNA network can provide a reference for further research on the pathogenesis of gout.

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.

A ceRNA regulatory network in systemic lupus erythematosus and its molecular interplay with cancer.

BackgroundSystemic lupus erythematosus (SLE) is an autoimmune disease defined by the production of autoantibodies and involves multiple organs and systems. Although there are reports on SLE, data on its pathogenesis is limited.MethodsUsing R language software, we constructed a competing endogenous RNA (ceRNA) network. We then utilized the Search Tool for Recurring Instances of Neighbouring Genes (STRING) and cytoHubba databases to generate a protein-protein interaction (PPI) network, which led to the identification of hub genes. The top two hub genes with the highest Maximal Clique Centrality (MCC) score in the PPI network were further validated via quantitative real-time polymerase chain reaction (qRT-PCR) using in-house clinical samples. Also, weighted gene co-expression network analysis (WGCNA) with genes from the Gene Expression Omnibus Series (GSE)121239 dataset identified hub modules that were associated with clinical indicators. In addition, the genes contained in key modules as obtained by WGCNA were enriched and analyzed using the Database for Annotation, Visualization and Integrated Discovery (DAVID) online tool. The top hub gene, X-linked apoptosis inhibitory protein-associated factor (XAF1), was then identified by intersection of the PPI and WGCNA outcomes, and a pan-cancer analysis of this hub gene was subsequently performed.ResultsWe comprehensively profiled the expression of Circular RNAs (circRNAs), MicroRNAs (miRNAs), and messenger RNAs (mRNAs) in SLE. We identified a hub gene, XAF1, based on evidence from the ceRNA network, WGCNA key module genes, and PPI network analyses. Moreover, qRT-PCR analysis demonstrated that the expression of XAF1 was significantly upregulated in SLE. Through the pan-cancer analysis, we demonstrated the common molecular roles of XAF1 in the pathogenesis of SLE and tumors, especially cutaneous melanoma.ConclusionsXAF1 is a key molecular biomarker in SLE. The pan-cancer analysis in this study provided shared genomic characteristics in SLE and cancers, especially for skin cutaneous melanoma (SKCM).

An analysis of lncRNA-miRNA-mRNA networks to investigate the effects of HDAC4 inhibition on skeletal muscle atrophy caused by peripheral nerve injury.

BackgroundMuscle atrophy caused by peripheral nerve injury is a common clinical disease, with no effective treatments currently available. Our previous studies have found that denervation-induced muscle atrophy can be alleviated by inhibiting histone deacetylase 4 (HDAC4). An increasing amount of evidence shows that microRNA (miRNA) and long noncoding RNA (lncRNA) are involved in the occurrence of muscle atrophy. This study aimed to find the mechanism by which HDAC4 regulates denervation-induced muscle atrophy based on lncRNA-associated competing endogenous RNA (ceRNA) networks.MethodsWe analyzed the influence of short hairpin RNA (shRNA) knockdown of HDAC4 on lncRNAs and miRNAs after denervated muscle atrophy using RNA sequencing. A Pearson’s correlation heat map and principal component analysis were employed to analyze differentially expressed miRNAs and lncRNAs. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses of target genes were conducted. The ceRNA network of lncRNA-miRNA-mRNA was constructed, and the core regulatory molecules in the ceRNA network were analyzed.ResultsWe found 32 miRNAs and 111 lncRNAs related to denervated muscle atrophy regulated by HDAC4. Moreover, 15 downregulated lncRNAs, 14 upregulated miRNAs, and 61 downregulated mRNAs constituted a ceRNA regulatory network, participating in the biological processes including response to denervation involved in regulation of muscle adaptation, along with the signaling pathways including autophagy, FoxO signaling pathways, and Jak-STAT signaling pathways. Additionally, 6 upregulated lncRNAs, 8 downregulated miRNAs, and 66 upregulated mRNAs constituted another ceRNA regulatory network, which was mainly involved in cell cycle-related biological processes and pathways. Finally, 3 lncRNAs, 4 miRNAs, and 12 mRNAs constituted a ceRNA sub-network, and XR_377582.2 and ENSMUST00000143649 were considered to be the key lncRNAs.ConclusionsIn the ceRNA network, all nodes are directly or indirectly involved in the process by which HDAC4 regulates skeletal muscle atrophy caused by peripheral nerve injury. XR_377582.2 and ENSMUST00000143649 may be the key lncRNAs related to HDAC4 involved in the regulation of muscle atrophy.

Construction of PARPi Resistance-related Competing Endogenous RNA Network.

Objective: Ovarian cancer is a kind of common gynecological malignancy in women. PARP inhibitors (PARPi) have been approved for ovarian cancer treatment. However, the primary and acquired resistance have limited the application of PARPi. The mechanisms remain to be elucidated. Methods: In this study, we characterized the expression profiles of mRNA and nonconding RNAs (ncRNAs) and constructed the regulatory networks based on RNA sequencing in PARPi Olaparib-induced ovarian cancer cells. Results: We found that the functions of the differentially expressed genes were enriched in "PI3K/AKT signaling pathway," "MAPK signaling pathway" and "metabolic process". The functions of DELs (cis) were enriched in "Human papillomavirus infection""tight junction" "MAPK signaling pathway". As the central regulator of ceRNAs, the differentially expressed miRNAs were enriched in "Human papillomavirus infection" "MAPK signaling pathway" "Ras signaling pathway". According to the degree of interaction, we identified 3 lncRNAs, 2 circRNAs, 7 miRNAs, and 12 mRNA as the key regulatory ceRNA axis, in which miR-320b was the important mediator. Conclusion: Here, we revealed the key regulatory lncRNA (circRNA)-miRNA-mRNA axis and their involved pathways in the PARPi resistant ovarian cancer cells. These findings provide new insights into exploring the ceRNA regulatory networks and developing new targets for PARPi resistance.

Construction of a ceRNA Network and Comprehensive Analysis of lncRNA in Hepatocellular Carcinoma.

To explore the RNA biomolecular marker associated with hepatocellular carcinoma (HCC) prognosis, we constructed a regulatory network of competitive endogenous RNAs (ceRNAs), which provides favorable conditions for the early diagnosis, prognostic monitoring, and personalized treatment of HCC. In this study, the differentially expressed genes (DEGs) of patients with HCC were obtained from the Gene Expression Omnibus. We identified 574 upregulated genes and 274 downregulated genes relevant to HCC occurrence (p < 0.05). Subsequently, we constructed the protein–protein interaction (PPI) network using these DEGs and identified the hub genes from the PPI. We then determined the expression and prognostic values of the hub genes from the GEPIA and Kaplan–Meier plotter databases. After the upstream microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) were respectively identified by miRTarBase and miRNet, we validated the expression of the key miRNAs in the serum using qPCR experiments. Moreover, we identified a two-lncRNA (LINC01184 and ADORA2A-AS1) signature from the upstream lncRNA that effectively predicted overall survival and had promotive effects for HCC. To verify the clinical significance of the signature, we validated the expression of the lncRNA in HCC tissues. Finally, we discovered and identified four mRNAs, four miRNAs, and five lncRNAs associated with the prognosis of HCC and constructed a new ceRNA regulatory network, which will be beneficial for the accurate diagnosis and treatment of HCC.

Comprehensive analysis of GSEC/miR-101-3p/SNX16/PAPOLG axis in hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) is one of the most lethal malignancies. A growing number of studies have shown that competitive endogenous RNA (ceRNA) regulatory networks might play important roles during HCC process. The present study aimed to identify a regulatory axis of the ceRNA network associated with the development of HCC. The roles of SNX16 and PAPOLG in HCC were comprehensively analyzed using bioinformatics tools. Subsequently, the “mRNA-miRNA-lncRNA” model was then used to predict the upstream miRNAs and lncRNAs of SNX16 and PAPOLG using the miRNet database, and the miRNAs with low expression and good prognosis in HCC and the lncRNAs with high expression and poor prognosis in HCC were screened by differential expression and survival analysis. Finally, the risk-prognosis models of ceRNA network axes were constructed by univariate and multifactorial Cox proportional risk analysis, and the immune correlations of ceRNA network axes were analyzed using the TIMER and GEPIA database. In this study, the relevant ceRNA network axis GSEC/miR-101-3p/SNX16/PAPOLG with HCC prognosis was constructed, in which GSEC, SNX16, and PAPOLG were highly expressed in HCC with poor prognosis, while miR-101-3p was lowly expressed in HCC with good prognosis. The risk-prognosis model predicted AUC of 0.691, 0.623, and 0.626 for patient survival at 1, 3, and 5 years, respectively. Immuno-infiltration analysis suggested that the GSEC/miR-101-3p/SNX16/PAPOLG axis might affect macrophage polarization. The GSEC/miR-101-3p/SNX16/PAPOLG axis of the ceRNA network axis might be an important factor associated with HCC prognosis and immune infiltration.