Integrative analysis of DNA methylation-driven genes for the prognosis of lung squamous cell carcinoma using MethylMix

Immunogenic cell death-related genes predict prognosis and response to immunotherapy in lung squamous cell carcinoma

Lung squamous cell carcinoma (LUSC) is a malignancy with limited therapeutic options. Immunogenic cell death (ICD) has the potential to enhance the efficacy of cancer therapy by triggering immune responses. We aimed to explore the potential of ICD-based classification in predicting prognosis and response to immunotherapy for LUSC. RNA-seq information and clinical data of LUSC patients were obtained from The Cancer Genome Atlas (TCGA) dataset. ICD-related gene expressions in LUSC samples were analyzed by consensus clustering. Subsequently, differentially expressed genes (DEGs) between different ICD-related subsets were analyzed. Tumor mutation burden, immune cell infiltration, and survival analyses were conducted between different ICD subsets. Finally, an ICD-related risk signature was constructed and evaluated in LUSC patients, and the immunotherapy responses based on the gene expressions were also forecasted. ICD-high and ICD-low groups were defined, and 1466 DEGs were identified between the two subtypes. These DEGs were mainly enriched in collagen-containing extracellular matrix, cytokine-cytokine receptor interaction, the PI3K-Akt signaling pathway, and neuroactive ligand-receptor interaction. Furthermore, the ICD-low group exhibited a favorable prognosis, enhanced TTN and MUC16 mutation frequencies, increased infiltrating immune cells, and downregulated immune checkpoint expressions. Furthermore, we demonstrated that an ICD-related model (based on CD4, NLRP3, NT5E, and TLR4 genes) could forecast the prognosis of LUSC, and ICD risk scores were lower in the responder group. In summary, the predicted values of ICD-related genes (CD4, NLRP3, NT5E, and TLR4) for the prognosis and response to immunotherapy in LUSC were verified in the study, which benefits immunotherapy-based interventions for LUSC patients.

Identification of Key DNA methylation sites related to differentially expressed genes in Lung squamous cell carcinoma

Changes in DNA methylation level at some CpG locus are closely associated with the occurrence of lung squamous cell carcinoma (LUSC). However, its specific regulatory mechanism is still unclear. Therefore, it is necessary to systematically identify and analyze those key CpG sites whose DNA methylation levels are closely related to the differential expression of up- and down-regulated genes in LUSC. Due to the dispersion of DNA methylation sites in different regions of genome, to study the correlation between gene expression level and DNA methylation, we divided gene into 6 non-overlapping functional regions and proposed a two-step correlation analysis method to identify differential DNA methylation sites and matched differential expression genes. As a results, we obtained 39 key CpG sites scattered in 27 genes. Through comparative analysis of LUSC-normal sample pairs, we found that these sites and genes can accurately cluster LUSC samples and normal samples. Finally, we used these sites and genes to distinguish LUSC from normal samples. The results suggest that they can be used as effective biomarkers for identifying LUSC. In addition, the proposed two-step correlation analysis method can also be extended to the identification of biomarkers of other cancers and diseases.

Epithelial-mesenchymal transition-related gene signature for prognosis of lung squamous cell carcinoma

Epithelial-mesenchymal transition (EMT) is associated with tumor invasion and progression, and is regulated by DNA methylation. A prognostic signature of lung squamous cell carcinoma (LUSC) with EMT-related gene data has not yet been established. In our study, we constructed a co-expression network using differentially expressed genes (DEGs) obtained from The Cancer Genome Atlas (TCGA) to identify hub genes. We conducted a correlation analysis between the differentially methylated hub genes and differentially expressed EMT-related genes to screen EMT-related differentially methylated genes (ERDMGs). Functional enrichment was performed to annotate the ERDMGs. The least absolute shrinkage and selection operator (LASSO) and stepwise Cox regression analyses were performed to build a survival prognosis prediction model. Additionally, druggability analysis was performed to predict the potential drug targets of ERDMGs. We screened 11 ERDMGs that were enriched in cell adhesion molecules and other signaling pathways. Finally, we constructed a 4-ERDMG model, which showed good ability to predict survival prognosis in the training and validation sets. The model could serve as an independent predictive factor for patients with LUSC. Additionally, our druggability analysis predicted that CC chemokine ligand 23 (CCL23) and Hepatocyte nuclear factor 1b (HNF1B) may be the underlying drug targets of LUSC. We established a new risk score (RS) system as a prognostic indicator to predict the outcome of patients with LUSC, which will help in the improvement of treatment strategies.

Immune-related gene-based prognostic index for predicting survival and immunotherapy outcomes in colorectal carcinoma

Introduction

Colorectal cancer shows high incidence and mortality rates. Immune checkpoint blockade can be used to treat colorectal carcinoma (CRC); however, it shows limited effectiveness in most patients.

Methods

To identify patients who may benefit from immunotherapy using immune checkpoint inhibitors, we constructed an immune-related gene prognostic index (IRGPI) for predicting the efficacy of immunotherapy in patients with CRC. Transcriptome datasets and clinical information of patients with CRC were used to identify differential immune-related genes between tumor and para-carcinoma tissue. Using weighted correlation network analysis and Cox regression analysis, the IRGPI was constructed, and Kaplan-Meier analysis was used to evaluate its predictive ability. We also analyzed the molecular and immune characteristics between IRGPI high-and low-risk subgroups, performed sensitivity analysis of ICI treatment, and constructed overall survival-related receiver operating characteristic curves to validate the IRGPI. Finally, IRGPI genes and tumor immune cell infiltration in CRC model mice with orthotopic metastases were analyzed to verify the results.

Results

The IRGPI was constructed based on the following 11 hub genes: ADIPOQ, CD36, CCL24, INHBE, UCN, IL1RL2, TRIM58, RBCK1, MC1R, PPARGC1A, and LGALS2. Patients with CRC in the high-risk subgroup showed longer overall survival than those in the low-risk subgroup, which was confirmed by GEO database. Clinicopathological features associated with cancer progression significantly differed between the high- and low-risk subgroups. Furthermore, Kaplan-Meier analysis of immune infiltration showed that the increased infiltration of naïve B cells, macrophages M1, and regulatory T cells and reduced infiltration of resting dendritic cells and mast cells led to a worse overall survival in patients with CRC. The ORC curves revealed that IRGPI predicted patient survival more sensitive than the published tumor immune dysfunction and rejection and tumor inflammatory signature.

Discussion

Thus, the low-risk subgroup is more likely to benefit from ICIs than the high-risk subgroup. CRC model mice showed higher proportions of Tregs, M1 macrophages, M2 macrophages and lower proportions of B cells, memory B cell immune cell infiltration, which is consistent with the IRGPI results. The IRGPI can predict the prognosis of patients with CRC, reflect the CRC immune microenvironment, and distinguish patients who are likely to benefit from ICI therapy.

Regulatory pattern of abnormal promoter CpG island methylation in the glioblastoma multiforme classification

Glioblastoma (GBM) is characterized by extensive genetic and phenotypic heterogeneity. However, it remains unexplored primarily how CpG island methylation abnormalities in promoter mediate glioblastoma typing. First, we presented a multi-omics scale map between glioblastoma sample clusters constructed based on promoter CpG island (PCGI) methylation-driven genes, using datasets including methylation profiles, expression profiles, and single-cell sequencing data from multiple highly annotated public clinical cohorts. Second, we identified differences in the tumor microenvironment between the two glioblastoma sample clusters and resolved key signaling pathways between cell clusters at the single-cell level based on comprehensive comparative analyses to investigate the reasons for survival differences between two of these clusters. Finally, we developed a diagnostic map and a prediction model for glioblastoma, and compared theoretical differences of drug sensitivity between two glioblastoma sample clusters. In summary, this study established a classification system for dissecting promoter CpG island methylation heterogeneity in glioblastoma and provides a new perspective for the diagnosis and treatment of glioblastoma.

Pan-Cancer Methylated Dysregulation of Long Non-coding RNAs Reveals Epigenetic Biomarkers

Different cancer types not only have common characteristics but also have their own characteristics respectively. The mechanism of these specific and common characteristics is still unclear. Pan-cancer analysis can help understand the similarities and differences among cancer types by systematically describing different patterns in cancers and identifying cancer-specific and cancer-common molecular biomarkers. While long non-coding RNAs (lncRNAs) are key cancer modulators, there is still a lack of pan-cancer analysis for lncRNA methylation dysregulation. In this study, we integrated lncRNA methylation, lncRNA expression and mRNA expression data to illuminate specific and common lncRNA methylation patterns in 23 cancer types. Then, we screened aberrantly methylated lncRNAs that negatively regulated lncRNA expression and mapped them to the ceRNA relationship for further validation. 29 lncRNAs were identified as diagnostic biomarkers for their corresponding cancer types, with lncRNA AC027601 was identified as a new KIRC-associated biomarker, and lncRNA ACTA2-AS1 was regarded as a carcinogenic factor of KIRP. Two lncRNAs HOXA-AS2 and AC007228 were identified as pan-cancer biomarkers. In general, the cancer-specific and cancer-common lncRNA biomarkers identified in this study may aid in cancer diagnosis and treatment.

A methylation-driven genes prognostic signature and the immune microenvironment in epithelial ovarian cancer

Aberrant gene methylation has been implicated in the development and progression of tumors. In this study, we aimed to identity methylation driven genes involved in epithelial ovarian cancer (EOC) to establish a prognostic signature for patients with EOC. We identified and verified 6 MDGs that are closely related to the prognosis of ovarian cancer. A prognostic risk score model and nomogram for predicting the prognosis of ovarian cancer were constructed based on the six MDGs. It can also effectively reflect the immune environment and immunotherapy response of ovarian cancer. These MDGs have great significance to the implementation of individualized treatment and disease monitoring of ovarian cancer patients.

Clinicopathological implications of lncRNAs, immunotherapy and DNA methylation in lung squamous cell carcinoma: a narrative review

Objective

To explore the clinicopathological impact of lncRNAs, immunotherapy, and DNA methylation in lung squamous cell carcinoma (LUSC), emphasizing their exact roles in carcinogenesis and modes of action.

Background

LUSC is the second most prevalent form, accounting for around 30% of non-small cell lung cancer (NSCLC). To date, molecular-targeted treatments have significantly improved overall survival in lung adenocarcinoma patients but have had little effect on LUSC therapy. As a result, there is an urgent need to discover new treatments for LUSC that are based on existing genomic methods.

Methods

In this review, we summarized and analyzed recent research on the biological activities and processes of lncRNA, immunotherapy, and DNA methylation in the formation of LUSC. The relevant studies were retrieved using a thorough search of Pubmed, Web of Science, Science Direct, Google Scholar, and the university's online library, among other sources.

Conclusions

LncRNAs are the primary components of the mammalian transcriptome and are emerging as master regulators of a number of cellular processes, including the cell cycle, differentiation, apoptosis, and growth, and are implicated in the pathogenesis of a variety of cancers, including LUSC. Understanding their role in LUSC in detail may help develop innovative treatment methods and tactics for LUSC. Meanwhile, immunotherapy has transformed the LUSC treatment and is now considered the new standard of care. To get a better knowledge of LUSC biology, it is critical to develop superior modeling systems. Preclinical models, particularly those that resemble human illness by preserving the tumor immune environment, are essential for studying cancer progression and evaluating novel treatment targets. DNA methylation, similarly, is a component of epigenetic alterations that regulate cellular function and contribute to cancer development. By methylating the promoter regions of tumor suppressor genes, abnormal DNA methylation silences their expression. DNA methylation indicators are critical in the early detection of lung cancer, predicting therapy efficacy, and tracking treatment resistance. As such, this review seeks to explore the clinicopathological impact of lncRNAs, immunotherapy, and DNA methylation in LUSC, emphasizing their exact roles in carcinogenesis and modes of action.

Specific Lung Squamous Cell Carcinoma Prognosis-Subtype Distinctions Based on DNA Methylation Patterns

Background

Lung squamous cell carcinoma (LUSC) is one of the major types of non-small-cell lung cancer. Epigenetic alterations, such as DNA methylation, have been recognized to be closely associated with the tumorigenesis and progression.

Material/Methods

In this study, we investigated the prognosis subgroups and assessed their correlation with clinical characteristics in LUSC using a methylation array acquired from The Cancer Genome Atlas (TCGA) database.

Results

A total of 196 DNA methylation sites exhibited a significant association with patient prognosis, and patients were further stratified into 7 prognosis subgroups based upon the consensus clustering. The patients in every subgroup were different in terms of prognosis and TNM stage. In addition, we found these 196 significant methylation sites corresponded to 258 genes. The function enrichment analysis revealed that these 258 genes enriched in biological pathways were closely related to cancers, such as DNA methylation and demethylation, cell cycle DNA replication, regulation of signal transduction by p53 class mediator, and genetic imprinting. Subsequently, we determined the levels of methylation sites in 7 subgroups, and found 24 intra-subgroup-specific methylation sites. Meanwhile, we selected 3 subgroups-specific methylation sites to construct the prognosis model for LUSC patients using multivariate Cox proportional risk regression model analysis. This model can effectively predict the prognosis of LUSC patients.

Conclusions

Our study identified a new classification of LUSC into 7 prognosis subgroups on the basis of DNA methylation data in TCGA, which demonstrated that molecular subtypes are independent factor for prognosis in LUSC. This may provide a more detailed explanation for LUSC heterogeneity. Additionally, this classification will contribute to discovery of new biomarkers of LUSC and provide more accurate subdivision of LUSC. Furthermore, these specific DNA methylation sites and corresponding genes can serve as biomarkers for early diagnosis, accurate therapy, and prognosis prediction.

Epigenetic Regulatory Enzymes: mutation Prevalence and Coexistence in Cancers

Epigenetic regulation is an important layer of transcriptional control with the particularity to affect the broad spectrum of genome. Over the years, largely due to the substantial number of recurrent mutations, there have been hundreds of novel driver genes characterized in various cancers. Additionally, the relative contribution of two dysregulated epigenomic entities (DNA methylation and histone modifications) that gradually drive the cancer phenotype remains in the research focus. However, a complex scenario arises when the disease phenotype does not harbor any relevant mutation or an abnormal transcription level. Although the cancer landscape involves the contribution of multiple genetic and non-genetic factors, herein, we discuss specifically the mutation spectrum of epigenetically-related enzymes in cancer. In addition, we address the coexistence of these two epigenetic entities in malignant human diseases, especially cancer. We suggest that the study of epigenetically-related somatic mutations in the early cellular differentiation stage of embryonic development might help to understand their later-staged footprints in the cancer genome. Furthermore, understanding the co-occurrence and/or inverse association of different disease types and redefining the general definition of "healthy" controls could provide insights into the genome reorganization.

Integral Analysis of the RNA Binding Protein-associated Prognostic Model for Renal Cell Carcinoma

RNA binding protein (RBPs) dysregulation has been reported in various malignant tumors and plays a pivotal role in tumor carcinogenesis and progression. However, the underlying mechanisms in renal cell carcinoma (RCC) are still unknown. In the present study, we performed a bioinformatics analysis using data from TCGA database to explore the expression and prognostic value of RBPs. We identified 125 differently expressed RBPs between tumor and normal tissue in RCC patients, including 87 upregulated and 38 downregulated RBPs. Eight RBPs (RPL22L1, RNASE2, RNASE3, EZH2, DDX25, DQX1, EXOSC5, DDX47) were selected as prognosis-related RBPs and used to construct a risk score model. In the risk score model, the high-risk subgroup had a poorer overall survival (OS) than the low-risk subgroup, and we divided the 539 RCC patients into two groups and conducted a time-dependent receiver operating characteristic (ROC) analysis to further test the prognostic ability of the eight hub RBPs. The area under the curve (AUC) of the ROC curve was 0.728 in train-group and 0.688 in test-group, indicating a good prognostic model. More importantly, we established a nomogram based on the selected eight RBPs. The eight selected RBPS have predictive value for RCC patients, with potential applications in clinical decision-making and individualized treatment.

Potential Genes Associated with the Survival of Lung Adenocarcinoma Were Identified by Methylation

Background

Lung adenocarcinoma (LUAD) is the most common pathological type of lung cancer. The purpose of this study is to search for genes related to the prognosis of LUAD through methylation based on a linear mixed model (LMM).

Methods

Gene expression, methylation, and survival data of LUAD patients were downloaded from the TCGA database. Based on the LMM model, the GEMMA algorithm was used to screen the predictive genes related to LUAD survival. The Cox model was used to further screen the predicted genes, and then, protein-protein interaction (PPI) network was constructed. Through the software plugin Cytoscape MCODE 3.8.0, the most closely related genes in the PPI network module were selected for in-depth biological function analysis to further explore the interaction and correlation between genes.

Results

We screened out 97 predictive genes from 18,834 genes and eliminated one gene associated with lung squamous cell carcinoma from previous studies, leaving 96 genes. The MCODE and the Kaplan-Meier curve analysis were used to finally identify two genes ASB16 and NEDD4 that are related to the prognosis of LUAD.

Conclusions

The newly identified two genes associated with the prognosis of LUAD may provide a basis for the treatment of patients.