Identification of a ferroptosis-related gene pair biomarker with immune infiltration landscapes in ischemic stroke: a bioinformatics-based comprehensive study

Background Ischemic stroke (IS) is a principal contributor to long-term disability in adults. A new cell death mediated by iron is ferroptosis, characterized by lethal aggregation of lipid peroxidation. However, a paucity of ferroptosis-related biomarkers early identify IS until now. This study investigated potential ferroptosis-related gene pair biomarkers in IS and explored their roles in immune infiltration. Results In total, we identified 6 differentially expressed ferroptosis-related genes (DEFRGs) in the metadata cohort. Of these genes, 4 DEFRGs were incorporated into the competitive endogenous RNA (ceRNA) network, including 78 lncRNA-miRNA and 16 miRNA-mRNA interactions. Based on relative expression values of DEFRGs, we constructed gene pairs. An integrated scheme consisting of machine learning algorithms, ceRNA network, and gene pair was proposed to screen the key DEFRG biomarkers. The receiver operating characteristic (ROC) curve witnessed that the diagnostic performance of DEFRG pair CDKN1A/JUN was superior to that of single gene. Moreover, the CIBERSORT algorithm exhibited immune infiltration landscapes: plasma cells, resting NK cells, and resting mast cells infiltrated less in IS samples than controls. Spearman correlation analysis confirmed a significant correlation between plasma cells and CDKN1A/JUN (CDKN1A: r = − 0.503, P < 0.001, JUN: r = − 0.330, P = 0.025). Conclusions Our findings suggested that CDKN1A/JUN could be a robust and promising gene-pair diagnostic biomarker for IS, regulating ferroptosis during IS progression via C9orf106/C9orf139-miR-22-3p-CDKN1A and GAS5-miR-139-5p/miR-429-JUN axes. Meanwhile, plasma cells might exert a vital interplay in IS immune microenvironment, providing an innovative insight for IS therapeutic target.

Identification of Characteristic Genes in Whole Blood of Intervertebral Disc Degeneration Patients by Weighted Gene Coexpression Network Analysis (WGCNA)

Intervertebral disc degeneration (IDD) is a major cause of lower back pain. However, to date, the molecular mechanism of the IDD remains unclear. Gene expression profiles and clinical traits were downloaded from the Gene Expression Omnibus (GEO) database. Firstly, weighted gene coexpression network analysis (WGCNA) was used to screen IDD-related genes. Moreover, least absolute shrinkage and selection operator (LASSO) logistic regression and support vector machine (SVM) algorithms were used to identify characteristic genes. Furthermore, we further investigated the immune landscape by the Cell-type Identification By Estimating Relative Subsets Of RNA Transcripts (CIBERSORT) algorithm and the correlations between key characteristic genes and infiltrating immune cells. Finally, a competing endogenous RNA (ceRNA) network was established to show the regulatory mechanisms of characteristic genes. A total of 2458 genes were identified by WGCNA, and 48 of them were disordered. After overlapping the genes obtained by LASSO and SVM-RFE algorithms, genes including LINC01347, ASAP1-IT1, lnc-SEPT7L-1, B3GNT8, CHRNB3, CLEC4F, LOC102724000, SERINC2, and LOC102723649 were identified as characteristic genes of IDD. Moreover, differential analysis further identified ASAP1-IT1 and SERINC2 as key characteristic genes. Furthermore, we found that the expression of both ASAP1-IT1 and SERINC2 was related to the proportions of T cells gamma delta and Neutrophils. Finally, a ceRNA network was established to show the regulatory mechanisms of ASAP1-IT1 and SERINC2. In conclusion, the present study identified ASAP1-IT1 and SERINC2 as the key characteristic genes of IDD through integrative bioinformatic analyses, which may contribute to the diagnosis and treatment of IDD.

LINC00467 Is Upregulated by DNA Copy Number Amplification and Hypomethylation and Shows ceRNA Potential in Lung Adenocarcinoma

Lung adenocarcinoma (LUAD) is the most common histological lung cancer, and it is the leading cause of cancer-related deaths worldwide. Long noncoding RNAs (lncRNAs) have been implicated in tumorigenesis. LINC00467 is a novel lncRNA that is abnormally expressed in several cancer types including LUAD. However, its function and regulatory mechanism in LUAD progression remain unclear. In this study, based on The Cancer Genome Atlas data mining, we demonstrated that DNA copy number amplification and hypomethylation was positively correlated with LINC00467 expression in LUAD. In addition, DNA copy number amplification was significantly associated with distant metastasis, immune infiltration and poor survival. Microarray analysis demonstrated that LINC00467 knockdown in the LUAD A549 cell line led to a distinct microRNA expression profile that impacted various target genes involved in multiple biological processes. This finding suggests that LINC00467 may regulate LUAD progression by functioning as a competing endogenous RNA (ceRNA). Finally, we constructed a ceRNA network that included two microRNAs (hsa-miR-1225-5p, hsa-miR-575) and five mRNAs (BARX2, BCL9, KCNK1, KIAA1324, TMEM182) specific to LINC00467 in LUAD. Subsequent Kaplan-Meier survival analysis in both The Cancer Genome Atlas and Gene Expression Omnibus databases revealed that two genes, BARX2 and BCL9, were potential prognostic biomarkers for LUAD patients. In conclusion, our data provide possible mechanisms underlying the abnormal upregulation of LINC00467 as well as a comprehensive view of the LINC00467-mediated ceRNA network in LUAD, thereby highlighting its potential role in diagnosis and therapy.

Comprehensive Analysis of Hub Genes Associated With Competing Endogenous RNA Networks in Stroke Using Bioinformatics Analysis

Background: Acute ischemic stroke (AIS) is the second leading cause of death and the third leading cause of disability worldwide. Long noncoding RNAs (lncRNAs) are promising biomarkers for the early diagnosis of AIS and closely participate in the mechanism of stroke onset. However, studies focusing on lncRNAs functioning as microRNA (miRNA) sponges to regulate the mRNA expression are rare and superficial.Methods: In this study, we systematically analyzed the expression profiles of lncRNA, mRNA (GSE58294), and miRNA (GSE110993) from the GEO database. Gene ontology (GO) analysis was performed to reveal the functions of differentially expressed genes (DEGs), and we used weighted gene co-expression network analysis (WGCNA) to investigate the relationships between clinical features and expression profiles and the co-expression of miRNA and lncRNA. Finally, we constructed a lncRNA–miRNA–mRNA competing endogenous RNA (ceRNA) network with selected DEGs using bioinformatics methods and obtained ROC curves to assess the diagnostic efficacy of differentially expressed lncRNAs (DElncRNAs) and differentially expressed mRNAs (DEmRNAs) in our network. The GSE22255 dataset was used to confirm the diagnostic value of candidate genes.Results: In total, 199 DElncRNAs, 2068 DEmRNAs, and 96 differentially expressed miRNAs were detected. The GO analysis revealed that DEmRNAs primarily participate in neutrophil activation, neutrophil degranulation, vacuolar transport, and lysosomal transport. WGCNA screened out 16 lncRNAs and 195 mRNAs from DEGs, and only eight DElncRNAs maintained an area under the curve higher than 0.9. By investigating the relationships between lncRNAs and mRNAs, a ceRNA network containing three lncRNAs, three miRNAs, and seven mRNAs was constructed. GSE22255 confirmed that RP1-193H18.2 is more advantageous for diagnosing stroke, whereas no mRNA showed realistic diagnostic efficacy.Conclusion: The ceRNA network may broaden our understanding of AIS pathology, and the candidate lncRNA from the ceRNA network is assumed to be a promising therapeutic target and diagnostic biomarker for AIS.

Analysis of CircRNAs and CircRNA-Associated Competing Endogenous RNA Networks in β-Thalassemia

The involvement of circRNAs in β-thalassemia and their actions on fetal hemoglobin (HbF) is unclear. Here, the circRNAs in β-thalassemia carriers with high HbF levels were comprehensively analyzed in comparison with healthy individuals. Differential expression of 2183 circRNAs was observed and their correlations with hematological parameters were investigated. Down-regulated has-circRNA-100466 had a strong negative correlation with HbF and HbA2. Bioinformatics was employed to construct a has-circRNA-100466‑associated competing endogenous RNA (ceRNA) network with the determination of hub genes and associated miRNAs. In combination with previous reports, the has-circRNA-100466▁miR-19b-3p▁SOX6 pathway was identified. The ceRNA network was verified by qRT-PCR on β-thalassemia samples and RNA immunoprecipitation of K562 cell lysates. Has-circRNA-100466, miR-19b-3p, and SOX6 were present together in anti-argonaute 2 immunoprecipitates, indicating involvement with HbF induction. Furthermore, spearman correlation coefficients revealed their significant correlations with HbF. In conclusion, a novel has-circRNA-100466▁miR-19b-3p▁SOX6 pathway was identified, providing insight into HbF induction and suggesting targets β-thalassemia treatment.

Identification of a Pyroptosis-Related lncRNA Risk Model for Predicting Prognosis and Immune Response in Colon Adenocarcinoma

Background: Colon adenocarcinoma (COAD) is one of the most common malignant tumors and diagnosed at an advanced stage with poor prognosis in the world. Pyroptosis is involved in the initiation and progression of tumors. This research focused on constructing a pyroptosis-related ceRNA network to generate a reliable risk model for risk prediction and immune infiltration analysis of COAD.Methods: Transcriptome data, miRNA-sequencing data and clinical information were downloaded from the TCGA database. Firstly, differentially expressed mRNAs (DEmRNAs), miRNAs (DEmiRNAs), and lncRNAs (DElncRNAs) were identified to construct a pyroptosis-related ceRNA network. Secondly, a pyroptosis-related lncRNA risk model was developed applying univariate Cox regression analysis and least absolute shrinkage and selection operator method (LASSO) regression analysis. The Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment analyses were utilized to functionally annotate RNAs contained in the ceRNA network. In addition, Kaplan-Meier analysis, receiver operating characteristic (ROC) curves, univariate and multivariate Cox regression, and nomogram were applied to validate this risk model. Finally, the relationship of this risk model with immune cells and immune checkpoint blockade (ICB) related genes were analyzed.Results: Totally 5373 DEmRNAs, 1159 DElncRNAs and 355 DEmiRNAs were identified. A pyroptosis-related ceRNA regulatory network containing 132 lncRNAs, 7miRNAs and 5 mRNAs was constructed and a ceRNA-based pyroptosis-related risk model including 11 lncRNAs was built. Tumor tissues were classified into high- and low- risk groups according to the median risk score. Kaplan-Meier analysis showed that the high-risk group had a shorter survival time; ROC analysis, independent prognostic analysis and nomogram further indicated the risk model was a significant independent prognostic factor which had excellent ability to predict patients’ risk. Moreover, immune infiltration analysis indicated that the risk model was related to immune infiltration cells (i.e., B cells naïve, T cells follicular helper, Macrophages M1, etc.) and ICB-related genes (i.e., PD-1, CTLA4, HAVCR2, etc).Conclusions: This pyroptosis-related lncRNA risk model possessed good prognostic value and the ability to predict the outcome of ICB immunotherapy in COAD.

Hsa_circ_0060927 Is a Novel Tumor Biomarker by Sponging miR-195-5p in the Malignant Transformation of OLK to OSCC

ObjectiveTo investigate the clinical significance of differentially expressed circRNAs and candidate circRNAs in the transformation of oral leukoplakia (OLK) to oral squamous cell carcinoma (OSCC).MethodsWe performed high-throughput circRNA sequencing in six cases of normal oral mucosal (NOM) tissues, six cases of OLK tissues, and six cases of OSCC tissues. Ten circRNAs with significant differential expression were verified by qRT-PCR. Enzyme tolerance assay and Sanger sequencing were performed on the screened target circRNA hsa_circ_0060927, and a qRT-PCR assay of hsa_circ_0060927 was performed in three tissues (24 cases in each group); this was followed by an ROC analysis. The ceRNA network was predicted using TargetScan and miRanda. MiR-195-5p and TRIM14 were selected as the downstream research objects of hsa_circ_0060927. The sponge mechanism of hsa_circ_0060927 was detected by AGO2 RIP. The interaction between hsa_circ_0060927 and miR-195-5p was verified by RNA pull-down assay and dual luciferase reporter gene assay. The expressions of hsa_circ_0060927, miR-195-5p, and TRIM14 were verified by normal oral epithelial primary cells and cell lines of LEUK1, SCC9, and SCC25. The hsa_circ_0060927 overexpressed plasmid and miR-195-5p mimics were constructed to transfection LEUK1 to detect the changes in cell proliferation, apoptosis, and migration.ResultsThe results of qRT-PCR validation were consistent with the sequencing results. Hsa_circ_0060927 is a true circRNA with trans-splicing sites. The expression of hsa_circ_0060927 increased in NOM, OLK, and OSCC. Overexpression of hsa_circ_0060927 enhanced the ability of cell proliferation and migration, and decreased cell apoptosis capacity. The prediction of ceRNA network suggested that hsa_circ_0060927 could regulate the target gene TRIM14 through sponging miR-195-5p. AGO2 RIP indicated that hsa_circ_0060927 had a sponge mechanism. RNA pull-down and dual luciferase reporter gene assay suggested that hsa_circ_0060927 interacted with miR-195-5p. Hsa_circ_0060927 was positively correlated with the expression of TRIM14, and could relieve the inhibition of miR-195-5p on TRIM14 to regulate cell proliferation, apoptosis, and migration of LEUK1 cells.ConclusionHsa_circ_0060927 acted as a potential key ceRNA to sponge downstream miR-195-5p and promote OLK carcinogenesis by upregulating TRIM14. Hsa_circ_0060927 was expected to be a molecular marker for the prevention and treatment of OLK carcinogenesis through the hsa_circ_0060927/miR-195-5p/TRIM14 axis.

Multivariate competing endogenous RNA network characterization for cancer MicroRNA biomarker discovery: a novel bioinformatics model with application to prostate cancer metastasis

Background MicroRNAs (miRNAs) are post-transcriptional regulators with the potential as biomarkers for cancer management. Data-driven competing endogenous RNA (ceRNA) network modeling is an effective way to decipher the complex interplay between miRNAs and spongers. However, no general rules are discovered for ceRNA network-based biomarker prioritization. Methods and Results In this study, a novel bioinformatics model was developed by integrating gene expression with multivariate miRNA-target data for ceRNA network-based biomarker discovery. Compared with traditional methods, the structural vulnerability in human lncRNA-miRNA-mRNA network was comprehensively analyzed, and the single-line regulatory or competing mode among miRNAs, lncRNAs and mRNAs was characterized and quantified as statistical evidence for miRNA biomarker identification. The application of this model to prostate cancer (PCa) metastasis identified a total of 12 miRNAs as putative biomarkers from metastatic PCa-specific lncRNA-miRNA-mRNA network and nine of them have been previously reported as biomarkers for PCa metastasis. The receiver operating characteristic curve and cell line qRT-PCR experiments demonstrated the power of miR-26b-5p, miR-130a-3p, and miR-363-3p as novel candidates for predicting PCa metastasis. Moreover, PCa-associated pathways such as prostate cancer signaling, ERK/MAPK signaling, and TGF-β signaling were significantly enriched by targets of identified miRNAs, indicating the underlying mechanisms of miRNAs in PCa carcinogenesis. Conclusions A novel ceRNA-based bioinformatics model was proposed and applied to screen candidate miRNA biomarkers for PCa metastasis. Functional validations using human samples and clinical data will be performed for future translational studies on the identified miRNAs.

ceRNA Network and Functional Enrichment Analysis of Preeclampsia by Weighted Gene Coexpression Network Analysis

Background. Preeclampsia (PE) is a multisystemic syndrome which has short- and long-term risk to mothers and children and has pluralistic etiology. Objective. This study is aimed at constructing a competitive endogenous RNA (ceRNA) network for pathways most related to PE using a data mining strategy based on weighted gene coexpression network analysis (WGCNA). Methods. We focused on pathways involving hypoxia, angiogenesis, and epithelial mesenchymal transition according to the gene set variation analysis (GSVA) scores. The gene sets of these three pathways were enriched by gene set enrichment analysis (GSEA). WGCNA was used to study the underlying molecular mechanisms of the three pathways in the pathogenesis of PE by analyzing the relationship among pathways and genes. The soft threshold power (β) and topological overlap matrix allowed us to obtain 15 modules, among which the red module was chosen for the downstream analysis. We chose 10 hub genes that satisfied ∣ log 2 Fold Change ∣ > 2 and had a higher degree of connectivity within the module. These candidate genes were subsequently confirmed to have higher gene significance and module membership in the red module. Coexpression networks were established for the hub genes to unfold the connection between the genes in the red module and PE. Finally, ceRNA networks were constructed to further clarify the underlying molecular mechanism involved in the occurrence of PE. 56 circRNAs, 17 lncRNAs, and 20 miRNAs participated in the regulation of the hub genes. Coagulation factor II thrombin receptor (F2R) and lumican (LUM) were considered the most relevant genes, and ceRNA networks of them were constructed. Conclusion. The microarray data mining process based on bioinformatics methods constructed lncRNA and miRNA networks for ten hub genes that were closely related to PE and focused on ceRNAs of F2R and LUM finally. The results of our study may provide insight into the mechanisms underlying PE occurrence.

Aberrant Expression of miR-1301 in Human Cancer

miR-1301 is a newly discovered miRNA, which is abnormally expressed in 14 types of tumors. miR-1301 inhibits 23 target genes, forms a ceRNA network with 2 circRNAs and 8 lncRNAs, and participates in 6 signaling pathways, thereby affecting tumor cell proliferation, invasion, metastasis, apoptosis, angiogenesis, etc. Abnormal expression of miR-1301 is often associated with poor prognosis of cancer patients. In addition, miR-1301 is related to the anti-tumor effect of epirubicin on osteosarcoma and imatinib on chronic myeloid leukemia(CML) and can enhance the cisplatin sensitivity of ovarian cancer. This work systematically summarizes the abnormal expression and prognostic value of miR-1301 in a variety of cancers, depicts the miR-1301-related signaling pathways and ceRNA network, and provides potential clues for future miR-1301 research.