Specific breast cancer prognosis-subtype distinctions based on DNA methylation patterns

Tumor Microenvironment Characterization in Glioblastoma Identifies Prognostic and Immunotherapeutically Relevant Gene Signatures

Tumor microenvironment (TME) cells are important elements in tumor tissue. There is increasing evidence that they have important clinical pathological significance in predicting tumor clinical outcomes and therapeutic effects. However, no systematic analysis of TME cell interactions in glioblastoma (GBM) has been reported. We systematically analyzed the transcriptional sequencing data of GBM to find an immune gene marker to predict the clinical results of GBM. First, we downloaded the expression profiles and clinical follow-up information of GBM from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO). CIBERSORT was used to evaluate the infiltration mode of TME in 757 patients, systematically correlated TME phenotype with genomic characteristics and clinicopathological characteristics of GBM, defined four TME phenotypes, and TMEScore was constructed using algorithms such as random forest and principal component analysis. There is a significant correlation between TMEScore and age of onset. High TMEScore samples are characterized by immune activation, TGF pathway activation, and high expression of immune checkpoint genes, while low TMEScore samples are characterized by high-frequency IDH1 and MET mutations. Therefore, a comprehensive landscape depicting the TME characteristics of GBM may help explain GBM's response to immunotherapy and provide new strategies for cancer treatment. In this study, TMEScore can be used as a new prognostic marker to predict the survival of GBM patients, and as a potential predictor of immune checkpoint inhibitor response.

Identification of immune-enhanced molecular subtype associated with BRCA1 mutations, immune checkpoints and clinical outcome in ovarian carcinoma

Ovarian carcinoma has the highest mortality among the malignant tumours in gynaecology, and new treatment strategies are urgently needed to improve the clinical status of ovarian carcinoma patients. The Cancer Genome Atlas (TCGA) cohort were performed to explore the immune function of the internal environment of tumours and its clinical correlation with ovarian carcinoma. Finally, four molecular subtypes were obtained based on the global immune‐related genes. The correlation analysis and clinical characteristics showed that four subtypes were all significantly related to clinical stage; the immune scoring results indicated that most immune signatures were upregulated in C3 subtype, and the majority of tumour‐infiltrating immune cells were upregulated in both C3 and C4 subtypes. Compared with other subtypes, C3 subtype had a higher BRCA1 mutation, higher expression of immune checkpoints, and optimal survival prognosis. These findings of the immunological microenvironment in tumours may provide new ideas for developing immunotherapeutic strategies for ovarian carcinoma.

Gene expression and DNA methylation analyses suggest that immune process-related ADCY6 is a prognostic factor of luminal-like breast cancer

Breast cancer is a malignant tumor that seriously threatens women's health, and luminal-like cancer subtypes account for the majority of the cases. The purpose of this study was to investigate the relationships among DNA methylation, gene expression profile, and the tumor-immune microenvironment of luminal-like breast cancer, and to identify the potential key genes that regulate immune cell infiltration in luminal-like breast cancer. The ESTIMATE algorithm was applied to calculate immune scores and stromal scores of patients with breast cancer. Kaplan-Meier curves were generated for survival analysis. The clusterProfile package was used for Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analysis. The protein-protein interaction (PPI) network was constructed using the STRING database and Cytoscape software. Correlations between ADCY6 expression and immune cell infiltration-related pathways were analyzed by gene set variation analysis. R software was used for the statistical analysis and figure generation. Disease-free survival was higher in the immune score-high group than it was in the immune score-low group, while the stromal score had no correlation with prognosis. There were 515 genes that differed in both gene expression and DNA methylation levels, and these genes were mainly enriched in immune process-related pathways. ADCY6 was enriched in module A of the PPI network. Patients with downregulation and hypermethylation of ADCY6 associated with a better prognosis. ADCY6 expression was negatively correlated with the activation of immune process-related signaling pathways, immune checkpoint receptors, and ligands, except for CLEC4G. DNA methylation was found to be involved in the regulation of the key cellular pathways of luminal-like breast cancer immune cell infiltration. Additionally, ADCY6 was identified as a prognostic factor involved in the DNA methylation-regulated immune processes in luminal-like breast cancer.

Identification of diagnostic biomarkers in patients with gestational diabetes mellitus based on transcriptome gene expression and methylation correlation analysis

BACKGROUND

Gestational diabetes mellitus (GDM) has a high prevalence in the period of pregnancy. However, the lack of gold standards in current screening and diagnostic methods posed the biggest limitation. Regulation of gene expression caused by DNA methylation plays an important role in metabolic diseases. In this study, we aimed to screen GDM diagnostic markers, and establish a diagnostic model for predicting GDM.

METHODS

First, we acquired data of DNA methylation and gene expression in GDM samples (N = 41) and normal samples (N = 41) from the Gene Expression Omnibus (GEO) database. After pre-processing the data, linear models were used to identify differentially expressed genes (DEGs). Then we performed pathway enrichment analysis to extract relationships among genes from pathways, construct pathway networks, and further analyzed the relationship between gene expression and methylation of promoter regions. We screened for genes which are significantly negatively correlated with methylation and established mRNA-mRNA-CpGs network. The network topology was further analyzed to screen hub genes which were recognized as robust GDM biomarkers. Finally, the samples were randomly divided into training set (N = 28) and internal verification set (N = 27), and the support vector machine (SVM) ten-fold cross-validation method was used to establish a diagnostic classifier, which verified on internal and external data sets.

RESULTS

In this study, we identified 465 significant DEGs. Functional enrichment analysis revealed that these genes were associated with Type I diabetes mellitus and immunization. And we constructed an interactional network including 1091 genes by using the regulatory relationships of all 30 enriched pathways. 184 epigenetics regulated genes were screened by analyzing the relationship between gene expression and promoter regions' methylation in the network. Moreover, the accuracy rate in the training data set was increased up to 96.3, and 82.1% in the internal validation set, and 97.3% in external validation data sets after establishing diagnostic classifiers which were performed by analyzing the gene expression profiles of obtained 10 hub genes from this network, combined with SVM.

CONCLUSIONS

This study provided new features for the diagnosis of GDM and may contribute to the diagnosis and personalized treatment of GDM.

Differential DNA methylation following chemotherapy for breast cancer is associated with lack of memory improvement at one year

The biological basis underlying cognitive dysfunction in women with early-stage breast cancer (BC) remains unclear, but could reflect gene expression changes that arise from the acquisition and long-term retention of soma-wide alterations in DNA methylation in response to chemotherapy. In this longitudinal study, we identified differences in peripheral methylation patterns present in women prior to treatment (T1) and 1 year after receiving chemotherapy (T4) and evaluated relationships among the differential methylation (DM) ratios with changes in cognitive function. A total of 58 paired (T1 and T4) blood specimens were evaluated. Methylation values were determined for DNA isolated from whole blood using a genome-wide array . Cognitive function was measured using the validated, computerized CNS Vital Signs platform. Relationships between methylation patterns and cognitive domain scores were compared using a stepwise linear regression analysis, with demographic variables as covariates. The symptom comparison analysis was restricted to 2,199 CpG positions showing significant methylation ratio changes between T1 and T4. The positions with DM were enriched for genes involved in the modulation of cytokine concentrations. Significant DM ratios were associated with memory domain (56 CpGs). Eight of the ten largest DM ratio changes associated with lack of memory improvement were localized to genes involved in either neural function (ECE2, PPFIBP2) or signalling processes (USP6NL, RIPOR2, KLF5, UBE2V1, DGKA, RPS6KA1). These results suggest that epigenetic changes acquired and retained for at least one year in non-tumour cells following chemotherapy may be associated with a lack of memory improvement following treatment in BC survivors.

Integrative analysis of DNA methylation and gene expression to identify key epigenetic genes in glioblastoma

Glioblastoma (GBM) ranks the most common and aggressive primary brain malignant tumor worldwide. However, the survival rates of patients remain very poor. Therefore, molecular oncology of GBM are urgently needed. In this study, we performed an integrative analysis of DNA methylation and gene expression to identify key epigenetic genes in GBM. The methylation and gene expression of GBM patients in The Cancer Genome Atlas (TCGA) database were downloaded. After data preprocessing, we identified 4,881 differentially expressed genes (DEGs) between tumor and normal samples, including 1,111 upregulated and 3,770 downregulated genes. Then, we randomly separated all samples into training set (n = 69) and testing set (n = 69). We next obtained 11,269 survival-methylation sites by univariate and multivariate Cox regression analyses. In the correlation analysis, we defined 198 low promoter methylation with high gene expression as epigenetically induced (EI) genes and 111 high promoter methylation with low gene expression as epigenetically suppressed (ES) genes. Key markers including C1orf61 and FAM50B were selected with a Pearson correlation coefficient greater than 0.75. Further, we chose the 20 CpG methylation sites of above two genes in unsupervised clustering analysis using the Euclidean distance. We found that the prognosis of the hypomethylated group was significantly better than that in the hypermethylated group (log-rank test p-value = 0.011). Based on the validation in the TCGA testing set and GEO dataset, we validated the prognostic value of our signature (p-value = 0.02 in TCGA and 0.012 in GEO). In conclusion, our findings provided predictive and prognostic value as methylation-based biomarkers for the diagnosis and treatment of GBM.

DNA Methylation Module Network-Based Prognosis and Molecular Typing of Cancer

Achieving cancer prognosis and molecular typing is critical for cancer treatment. Previous studies have identified some gene signatures for the prognosis and typing of cancer based on gene expression data. Some studies have shown that DNA methylation is associated with cancer development, progression, and metastasis. In addition, DNA methylation data are more stable than gene expression data in cancer prognosis. Therefore, in this work, we focused on DNA methylation data. Some prior researches have shown that gene modules are more reliable in cancer prognosis than are gene signatures and that gene modules are not isolated. However, few studies have considered cross-talk among the gene modules, which may allow some important gene modules for cancer to be overlooked. Therefore, we constructed a gene co-methylation network based on the DNA methylation data of cancer patients, and detected the gene modules in the co-methylation network. Then, by permutation testing, cross-talk between every two modules was identified; thus, the module network was generated. Next, the core gene modules in the module network of cancer were identified using the K-shell method, and these core gene modules were used as features to study the prognosis and molecular typing of cancer. Our method was applied in three types of cancer (breast invasive carcinoma, skin cutaneous melanoma, and uterine corpus endometrial carcinoma). Based on the core gene modules identified by the constructed DNA methylation module networks, we can distinguish not only the prognosis of cancer patients but also use them for molecular typing of cancer. These results indicated that our method has important application value for the diagnosis of cancer and may reveal potential carcinogenic mechanisms.

Identification of core genes and clinical roles in pregnancy-associated breast cancer based on integrated analysis of different microarray profile datasets

More women are delaying child-birth. Thus, the diagnosis of pregnancy-associated breast cancer (PABC) will continue to increase. The aim of this study was to identify core candidate genes of PABC, and the relevance of the genes on the prognosis of PABC. GSE31192 and GSE53031 microarray profile datasets were downloaded from the Gene Expression Omnibus database and differentially expressed genes were analyzed using the R package and GEO2R tool. Then, Gene Ontology and Kyoto Encyclopedia of Gene and Genome pathway enrichment analyses were performed using the Database for Annotation, Visualization, and Integrated Discovery. Moreover, the Search Tool for the Retrieval of Interacting Genes and the Molecular Complex Detection Cytoscape software plug-in were utilized to visualize protein–protein interactions and to screen candidate genes. A total of 239 DEGs were identified in PABC, including 101 up-regulated genes mainly enriched in fatty acid activation and the fibroblast growth factor signaling pathway, while 138 down-regulated genes particularly involved in activation of DNA fragmentation factor and apoptosis-induced DNA fragmentation. Fourteen hub genes with a high degree of connectivity were selected, including CREB1, ARF3, UBA5, SIAH1, KLHL3, HECTD1, MMP9, TRIM69, MEX3C, ASB6, UBE2Q2, FBXO22, EIF4A3, and PXN. Overall survival (OS) analysis of core candidate genes was performed using the Gene Expression Profiling Interactive Analysis and UALCAN websites. High ASB6 expression was associated with worse OS of PABC patients. Molecular subtypes and menopause status were also associated with worse OS for PABC patients. In conclusion, ASB6 could be a potential predictor and therapeutic target in patient with PABC.