Identification of Methylation Markers and Differentially Expressed Genes with Prognostic Value in Breast Cancer

Comprehensive analysis of DNA methylation gene expression profiles in GEO dataset reveals biomarkers related to malignant transformation of sinonasal inverted papilloma

BACKGROUND

DNA methylation may be involved in the regulation of malignant transformation from sinonasal inverted papilloma (SNIP) to squamous cell carcinoma (SCC). The study of gene methylation changes and screening of differentially methylated loci (DMLs) are helpful to predict the possible key genes in the malignant transformation of SNIP-SCC.

MATERIALS AND METHODS

Microarray dataset GSE125399 was downloaded from the Gene Expression Omnibus (GEO) database and differentially methylated loci (DMLs) were analyzed using R language (Limma package). ClusterProfiler R package was used to perform Gene Ontology (GO) analysis on up-methylated genes and draw bubble maps. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway and its visualization analysis were analyzed to speculate the possible key Genes in SNIP-SCC malignant transformation. Subsequently, SNIP cases archived in our department were collected, tissue microarray was made, and immunohistochemical staining was performed to analyze the expression levels of UCKL1, GSTT1, HLA-G, MAML2 and NRGN in different grades of sinonasal papilloma tissues.

RESULTS

Analysis of dataset GSE125399 identified 56 DMLs, including 49 upregulated DMLs and 7 downregulated DMLs. Thirty-one genes containing upregulated DNA methylation loci and three genes containing downregulated DNA methylation loci were obtained by methylation microarray annotation analysis. In addition, KEGG pathway visualization analysis of 31 up-methylated genes showed that there were four significantly up-methylated genes including UCKL1, GSTT1, HLA-G and MAML2, and one significantly down-methylated gene NRGN. Subsequently, compared with non-neoplasia nasal epithelial tissues, the expression of HLA-G and NRGN was upregulated in grade I, II, III and IV tissues, while the expression of MAML2 was lost. The protein expression changes of MAML2 and NRGN were significantly negatively correlated with their gene methylation levels.

CONCLUSIONS

By analyzing the methylation dataset, we obtained four up-regulated methylation genes UCKL1, GSTT1, HLA-G and MAML2 and one down-regulated gene NRGN. MAML2, a tumor suppressor gene with high methylation modification but loss of protein expression, and NRGN, a tumor gene with low methylation modification but upregulated protein expression, can be used as biological indicators to judge the malignant transformation of SNIP-SCC.

A journey from omics to clinicomics in solid cancers: Success stories and challenges

The word 'cancer' encompasses a heterogenous group of distinct disease types characterized by a spectrum of pathological features, genetic alterations and response to therapies. According to the World Health Organization, cancer is the second leading cause of death worldwide, responsible for one in six deaths and hence imposes a significant burden on global healthcare systems. High-throughput omics technologies combined with advanced imaging tools, have revolutionized our ability to interrogate the molecular landscape of tumors and has provided unprecedented understanding of the disease. Yet, there is a gap between basic research discoveries and their translation into clinically meaningful therapies for improving patient care. To bridge this gap, there is a need to analyse the vast amounts of high dimensional datasets from multi-omics platforms. The integration of multi-omics data with clinical information like patient history, histological examination and imaging has led to the novel concept of clinicomics and may expedite the bench-to-bedside transition in cancer. The journey from omics to clinicomics has gained momentum with development of radiomics which involves extracting quantitative features from medical imaging data with the help of deep learning and artificial intelligence (AI) tools. These features capture detailed information about the tumor's shape, texture, intensity, and spatial distribution. Together, the related fields of multiomics, translational bioinformatics, radiomics and clinicomics may provide evidence-based recommendations tailored to the individual cancer patient's molecular profile and clinical characteristics. In this chapter, we summarize multiomics studies in solid cancers with a specific focus on breast cancer. We also review machine learning and AI based algorithms and their use in cancer diagnosis, subtyping, prognosis and predicting treatment resistance and relapse.

In silico analysis of differentially expressed-aberrantly methylated genes in breast cancer for prognostic and therapeutic targets

Breast cancer (BC) is the leading cause of death among women across the globe. Abnormal gene expression plays a crucial role in tumour progression, carcinogenesis and metastasis of BC. The alteration of gene expression may be through aberrant gene methylation. In the present study, differentially expressed genes which may be regulated by DNA methylation and their pathways associated with BC have been identified. Expression microarray datasets GSE10780, GSE10797, GSE21422, GSE42568, GSE61304, GSE61724 and one DNA methylation profile dataset GSE20713 were downloaded from Gene Expression Omnibus database (GEO). Differentially expressed-aberrantly methylated genes were identified using online Venn diagram tool. Based on fold change expression of differentially expressed-aberrantly methylated genes were chosen through heat map. Protein-protein interaction (PPI) network of the hub genes was constructed by Search Tool for the Retrieval of Interacting Genes (STRING). Gene expression and DNA methylation level of the hub genes were validated through UALCAN. Overall survival analysis of the hub genes was analysed through Kaplan-Meier plotter database for BC. A total of 72 upregulated-hypomethylated genes and 92 downregulated-hypermethylated genes were obtained from GSE10780, GSE10797, GSE21422, GSE42568, GSE61304, GSE61724, and GSE20713 datasets by GEO2R and Venn diagram tool. PPI network of the upregulated-hypomethylated hub genes (MRGBP, MANF, ARF3, HIST1H3D, GSK3B, HJURP, GPSM2, MATN3, KDELR2, CEP55, GSPT1, COL11A1, and COL1A1) and downregulated-hypermethylated hub genes were constructed (APOD, DMD, RBPMS, NR3C2, HOXA9, AMKY2, KCTD9, and EDN1). All the differentially expressed hub genes expression was validated in UALCAN database. 4 in 13 upregulated-hypomethylated and 5 in 8 downregulated-hypermethylated hub genes to be significantly hypomethylated or hypermethylated in BC were confirmed using UALCAN database (p < 0.05). MANF, HIST1H3D, HJURP, GSK3B, GPSM2, MATN3, KDELR2, CEP55, COL1A1, APOD, RBPMS, NR3C2, HOXA9, ANKMY2, and EDN1 were significantly (p < 0.05) associated with poor overall survival (OS). The identified aberrantly methylated-differentially expressed genes and their related pathways and function in BC can serve as novel diagnostic and prognostic biomarkers and therapeutic targets.Please confirm if the author names are presented accurately and in the correct sequence (given name, middle name/initial, family name). Author 4 Given name: [Jeewan Ram] Last name [Vishnoi]. Also, kindly confirm the details in the metadata are correct.It is correct.

Identification of a 5-gene-risk score model for predicting luminal A-invasive lobular breast cancer survival

Breast cancer is a devastating malignancy, among which the luminal A (LumA) breast cancer is the most common subtype. In the present study, we used a comprehensive bioinformatics approach in the hope of identifying novel prognostic biomarkers for LumA breast cancer patients. Transcriptomic profiling of 611 LumA breast cancer patients was downloaded from TCGA database. Differentially expressed genes (DEGs) between tumor samples and controls were first identified by differential expression analysis, before being used for the weighted gene co-expression network analysis. The subsequent univariate Cox regression and LASSO algorithm were used to uncover key prognostic genes for constructing multivariate Cox regression model. Patients were stratified into high-risk and low-risk groups according to the risk score, and subjected to multiple downstream analyses including survival analysis, gene set enrichment analysis (GSEA), inference on immune cell infiltration and analysis of mutation burden. Receiving operator curve analysis was also performed. A total of 7071 DEGs were first identified by edgeR package, pink module was found significantly associated with invasive lobular carcinoma (ILC). 105 prognostic genes and 9 predictors were identified, allowing the identification of a 5-key prognostic genes (LRRC77P, CA3, BAMBI, CABP1, ATP8A2) after intersection. These 5 genes, and the resulting Cox model, displayed good prognostic performance. Furthermore, distinct differences existed between two risk-score stratified groups at various levels. The identified 5-gene prognostic model will help deepen the understanding of the molecular and immunological mechanisms that affect the survival of LumA-ILC patients and guide and proper monitoring of these patients.

Comprehensive Analysis of Regulatory Factors and Immune-Associated Patterns to Decipher Common and BRCA1/2 Mutation-Type-Specific Critical Regulation in Breast Cancer

Background: BRCA1/2 mutations are closely related to high lifetime risk of breast cancer (BC). The objective of this study was to identify the genes, regulators, and immune-associated patterns underlying disease pathology in BC with BRCA1/2 somatic mutations and their associations with clinical traits. Methods: RNA sequencing data and clinical information from The Cancer Genome Atlas (TCGA; N = 36 BRCA1-mutant BC; N = 49 BRCA2-mutant BC; and N = 117 BRCA1/2-wild-type BC samples) were used for discovery, which included consensus network analysis, function enrichment, and analysis of hub genes; other TCGA data (N = 117 triple-negative BC) and two Gene Expression Omnibus database expression profiles were used as validation cohorts. Results: Consensus network analysis helped to identify specific co-expressed modules that showed positive correlations with tumor stage, number of positive lymph nodes, and margin status in BRCA1/2-mutant BC but lacking correlations in BRCA1/2-wild-type BC. Functional enrichment suggested potential mechanisms in BRCA1/2 carriers that could regulate the cell cycle, immune response, cellular metabolic processes, and cell migration, via enriched pathways including p53 and JAK-STAT signaling. Consensus network analysis identified the specific and common carcinogenic mechanisms involving BRCA mutations. Regulators cross-linking these modules include E2F or IRF transcription factor family, associated with cell cycle or immune response regulation module, respectively. Eight hub genes, including ISG15, BUB1, and TTK, were upregulated in several BRCA1/2-mutant BC datasets and showed prognostic value in BC. Furthermore, their genetic expression was related to higher levels of immune infiltration in BRCA1/2-mutant BC, which manifested as recruitment of T helper cells (Th1 cells), follicular helper T cells, and regulatory T cells, and T cell exhaustion. Moreover, important indicators for evaluation of BC immunotherapy, tumor mutational burden and neoantigen load also positively correlated with expression of some hub genes. Conclusion: We constructed a BRCA1/2 mutation-type-specific co-expressed gene network with related transcription factors and immune-associated patterns that could regulate and influence tumor metastasis and immune microenvironment, providing novel insights into the pathological process of this disease and the corresponding BRCA mutations.

Pan-cancer analysis of m5C regulator genes reveals consistent epigenetic landscape changes in multiple cancers

BACKGROUND

5-Methylcytosine (m5C) is a reversible modification to both DNA and various cellular RNAs. However, its roles in developing human cancers are poorly understood, including the effects of mutant m5C regulators and the outcomes of modified nucleobases in RNAs.

METHODS

Based on The Cancer Genome Atlas (TCGA) database, we uncovered that mutations and copy number variations (CNVs) of m5C regulatory genes were significantly correlated across many cancer types. We then assessed the correlation between the expression of individual m5C regulators and the activity of related hallmark pathways of cancers.

RESULTS

After validating m5C regulators' expression based on their contributions to cancer development and progression, we observed their upregulation within tumor-specific processes. Notably, our research connected aberrant alterations to m5C regulatory genes with poor clinical outcomes among various tumors that may drive cancer pathogenesis and/or survival.

CONCLUSION

Our results offered strong evidence and clinical implications for the involvement of m5C regulators.

Integrative Analysis of DNA Methylation and Gene Expression to Determine Specific Diagnostic Biomarkers and Prognostic Biomarkers of Breast Cancer

Background: DNA methylation is a common event in the early development of various tumors, including breast cancer (BRCA), which has been studies as potential tumor biomarkers. Although previous studies have reported a cluster of aberrant promoter methylation changes in BRCA, none of these research groups have proved the specificity of these DNA methylation changes. Here we aimed to identify specific DNA methylation signatures in BRCA which can be used as diagnostic and prognostic markers.

Methods: Differentially methylated sites were identified using the Cancer Genome Atlas (TCGA) BRCA data set. We screened for BRCA-differential methylation by comparing methylation profiles of BRCA patients, healthy breast biopsies and blood samples. These differential methylated sites were compared to nine main cancer samples to identify BRCA specific methylated sites. A BayesNet model was built to distinguish BRCA patients from healthy donors. The model was validated using three Gene Expression Omnibus (GEO) independent data sets. In addition, we also carried out the Cox regression analysis to identify DNA methylation markers which are significantly related to the overall survival (OS) rate of BRCA patients and verified them in the validation cohort.

Results: We identified seven differentially methylated sites (DMSs) that were highly correlated with cell cycle as potential specific diagnostic biomarkers for BRCA patients. The combination of 7 DMSs achieved ~94% sensitivity in predicting BRCA, ~95% specificity comparing healthy vs. cancer samples, and ~88% specificity in excluding other cancers. The 7 DMSs were highly correlated with cell cycle. We also identified 6 methylation sites that are highly correlated with the OS of BRCA patients and can be used to accurately predict the survival of BRCA patients (training cohort: likelihood ratio = 70.25, p = 3.633 × 10⁻¹³, area under the curve (AUC) = 0.784; validation cohort: AUC = 0.734). Stratification analysis by age, clinical stage, Tumor types, and chemotherapy retained statistical significance.

Conclusion: In summary, our study demonstrated the role of methylation profiles in the diagnosis and prognosis of BRCA. This signature is superior to currently published methylation markers for diagnosis and prognosis for BRCA patients. It can be used as promising biomarkers for early diagnosis and prognosis of BRCA.

Autoencoded DNA methylation data to predict breast cancer recurrence: Machine learning models and gene-weight significance

Breast cancer is the most frequent cancer in women and the second most frequent overall after lung cancer. Although the 5-year survival rate of breast cancer is relatively high, recurrence is also common which often involves metastasis with its consequent threat for patients. DNA methylation-derived databases have become an interesting primary source for supervised knowledge extraction regarding breast cancer. Unfortunately, the study of DNA methylation involves the processing of hundreds of thousands of features for every patient. DNA methylation is featured by High Dimension Low Sample Size which has shown well-known issues regarding feature selection and generation. Autoencoders (AEs) appear as a specific technique for conducting nonlinear feature fusion. Our main objective in this work is to design a procedure to summarize DNA methylation by taking advantage of AEs. Our proposal is able to generate new features from the values of CpG sites of patients with and without recurrence. Then, a limited set of relevant genes to characterize breast cancer recurrence is proposed by the application of survival analysis and a pondered ranking of genes according to the distribution of their CpG sites. To test our proposal we have selected a dataset from The Cancer Genome Atlas data portal and an AE with a single-hidden layer. The literature and enrichment analysis (based on genomic context and functional annotation) conducted regarding the genes obtained with our experiment confirmed that all of these genes were related to breast cancer recurrence.