Epigenetic and genetic burden measures are associated with tumor characteristics in invasive breast carcinoma

Associations in cell type-specific hydroxymethylation and transcriptional alterations of pediatric central nervous system tumors

Although intratumoral heterogeneity has been established in pediatric central nervous system tumors, epigenomic alterations at the cell type level have largely remained unresolved. To identify cell type-specific alterations to cytosine modifications in pediatric central nervous system tumors, we utilize a multi-omic approach that integrated bulk DNA cytosine modification data (methylation and hydroxymethylation) with both bulk and single-cell RNA-sequencing data. We demonstrate a large reduction in the scope of significantly differentially modified cytosines in tumors when accounting for tumor cell type composition. In the progenitor-like cell types of tumors, we identify a preponderance differential Cytosine-phosphate-Guanine site hydroxymethylation rather than methylation. Genes with differential hydroxymethylation, like histone deacetylase 4 and insulin-like growth factor 1 receptor, are associated with cell type-specific changes in gene expression in tumors. Our results highlight the importance of epigenomic alterations in the progenitor-like cell types and its role in cell type-specific transcriptional regulation in pediatric central nervous system tumors.

Alteration of DNA methyltransferases by eribulin elicits broad DNA methylation changes with potential therapeutic implications for triple-negative breast cancer

Background: Triple-negative breast cancer (TNBC) is an aggressive disease with limited treatment options. Eribulin, a chemotherapeutic drug, induces epigenetic changes in cancer cells, suggesting a unique mechanism of action. Materials & methods: MDA-MB 231 cells were treated with eribulin and paclitaxel, and the samples from 53 patients treated with neoadjuvant eribulin were compared with those from 14 patients who received the standard-of-care treatment using immunohistochemistry. Results: Eribulin treatment caused significant DNA methylation changes in drug-tolerant persister TNBC cells, and it also elicited changes in the expression levels of epigenetic modifiers (DNMT1, TET1, DNMT3A/B) in vitro and in primary TNBC tumors. Conclusion: These findings provide new insights into eribulin's mechanism of action and potential biomarkers for predicting TNBC treatment response.

Epigenetic Dysregulation and Its Correlation with the Steroidogenic Machinery Impacting Breast Pathogenesis: Data Mining and Molecular Insights into Therapeutics

Breast cancer (BC) is a heterogeneous condition and comprises molecularly distinct subtypes. An imbalance in the levels of epigenetic histone deacetylases (HDACs), modulating estrogen accumulation, especially 17β-estradiol (E2), promotes breast tumorigenesis. In the present study, analyses of The Cancer Genome Atlas (TCGA) pan-cancer normalized RNA-Seq datasets revealed the dysregulation of 16 epigenetic enzymes (among a total of 18 members) in luminal BC subtypes, in comparison to their non-cancerous counterparts. Explicitly, genomic profiling of these epigenetic enzymes displayed increases in HDAC1, 2, 8, 10, 11, and Sirtuins (SIRTs) 6 and 7, and decreases in HDAC4-7, -9, and SIRT1-4 levels, respectively, in TCGA breast tumors. Kaplan-Meier plot analyses showed that these HDACs, with the exception of HDAC2 and SIRT2, were not correlated with the overall survival of BC patients. Additionally, disruption of the epigenetic signaling in TCGA BC subtypes, as assessed using both heatmaps and boxplots, was associated with the genomic expression of factors that are instrumental for cholesterol trafficking/utilization for accelerating estrogen/E2 levels, in which steroidogenic acute regulatory protein (STAR) mediates the rate-limiting step in steroid biosynthesis. TCGA breast samples showed diverse expression patterns of a variety of key steroidogenic markers and hormone receptors, including LIPE, CYP27A1, STAR, STARD3, CYP11A1, CYP19A1, ER, PGR, and ERBB2. Moreover, regulation of STAR-governed steroidogenic machinery was found to be influenced by various transcription factors, i.e., CREB1, CREM, SF1, NR4A1, CEBPB, SREBF1, SREBF2, SP1, FOS, JUN, NR0B1, and YY1. Along these lines, ingenuity pathway analysis (IPA) recognized a number of new targets and downstream effectors influencing BCs. Of note, genomic, epigenomic, transcriptional, and hormonal anomalies observed in human primary breast tumors were qualitatively similar in pertinent BC cell lines. These findings identify the functional correlation between dysregulated epigenetic enzymes and estrogen/E2 accumulation in human breast tumors, providing the molecular insights into more targeted therapeutic approaches involving the inhibition of HDACs for combating this life-threatening disease.

Hydroxymethylation alterations in progenitor-like cell types of pediatric central nervous system tumors are associated with cell type-specific transcriptional changes

Although intratumoral heterogeneity has been established in pediatric central nervous system tumors, epigenomic alterations at the cell type level have largely remained unresolved. To identify cell type-specific alterations to cytosine modifications in pediatric central nervous system tumors we utilized a multi-omic approach that integrated bulk DNA cytosine modification data (methylation and hydroxymethylation) with both bulk and single-cell RNA-sequencing data. We demonstrate a large reduction in the scope of significantly differentially modified cytosines in tumors when accounting for tumor cell type composition. In the progenitor-like cell types of tumors, we identified a preponderance differential CpG hydroxymethylation rather than methylation. Genes with differential hydroxymethylation, like HDAC4 and IGF1R, were associated with cell type-specific changes in gene expression in tumors. Our results highlight the importance of epigenomic alterations in the progenitor-like cell types and its role in cell type-specific transcriptional regulation in pediatric CNS tumors.

Integrated Analysis of Multi-Omics Data to Establish a Hypoxia-Related Prognostic Model in Osteosarcoma

Background:

Osteosarcoma (OS) is the most common malignant bone tumor in clinical practice, and currently, the ability to predict prognosis in the diagnosis of OS is limited. There is an urgent need to find new diagnostic methods and treatment strategies for OS.

Material and methods:

We downloaded the multi-omics data for OS from the TARGET database. Prognosis-associated methylation sites were used to identify clustered subtypes of OS, and OS was classified into 3 subtypes (C1, C2, C3). Survival analysis showed significant differences between the C3 subtype and the other subtypes. Subsequently, differentially expressed genes (DEGs) across subtypes were screened and subjected to pathway enrichment analysis.

Results:

A total of 249 DEGs were screened from C3 subtype to other subtypes. Metabolic pathway enrichment analysis showed that DEGs were significantly enriched to the hypoxic pathway. Based on univariate and multivariate COX regression analysis, 12 genes from the hypoxia pathway were further screened and used to construct hypoxia-related prognostic model (HRPM). External validation of the HRPM was performed on the GSE21257 dataset. Finally, differences in survival and immune infiltration between high and low risk score groups were compared.

Conclusion:

In summary, we proposed a hypoxia-associated risk model based on a 12-gene expression signature, which is potentially valuable for prognostic diagnosis of patients with OS.

Genetic and Epigenetic Signatures Associated with the Divergence of Aquilegia Species

Widely grown in the Northern Hemisphere, the genus Aquilegia (columbine) is a model system in adaptive radiation research. While morphological variations between species have been associated with environmental factors, such as pollinators, how genetic and epigenetic factors are involved in the rapid divergence in this genus remains under investigated. In this study, we surveyed the genomes and DNA methylomes of ten Aquilegia species, representative of the Asian, European and North American lineages. Our analyses of the phylogeny and population structure revealed high genetic and DNA methylomic divergence across these three lineages. By multi-level genome-wide scanning, we identified candidate genes exhibiting lineage-specific genetic or epigenetic variation patterns that were signatures of inter-specific divergence. We demonstrated that these species-specific genetic variations and epigenetic variabilities are partially independent and are both functionally related to various biological processes vital to adaptation, including stress tolerance, cell reproduction and DNA repair. Our study provides an exploratory overview of how genetic and epigenetic signatures are associated with the diversification of the Aquilegia species.

Candidate methylation sites associated with endocrine therapy resistance in ER+/HER2- breast cancer

BACKGROUND

Estrogen receptor (ER) positive breast cancer is often effectively treated with drugs that inhibit ER signaling, i.e., tamoxifen (TAM) and aromatase inhibitors (AIs). However, 30% of ER+ breast cancer patients develop resistance to therapy leading to tumour recurrence. Changes in the methylation profile have been implicated as one of the mechanisms through which therapy resistance develops. Therefore, we aimed to identify methylation loci associated with endocrine therapy resistance.

METHODS

We used genome-wide DNA methylation profiles of primary ER+/HER2- tumours from The Cancer Genome Atlas in combination with curated data on survival and treatment to predict development of endocrine resistance. Association of individual DNA methylation markers with survival was assessed using Cox proportional hazards models in a cohort of ER+/HER2- tumours (N = 552) and two sub-cohorts corresponding to the endocrine treatment (AI or TAM) that patients received (N = 210 and N = 172, respectively). We also identified multivariable methylation signatures associated with survival using Cox proportional hazards models with elastic net regularization. Individual markers and multivariable signatures were compared with DNA methylation profiles generated in a time course experiment using the T47D ER+ breast cancer cell line treated with tamoxifen or deprived from estrogen.

RESULTS

We identified 134, 5 and 1 CpGs for which DNA methylation is significantly associated with survival in the ER+/HER2-, TAM and AI cohorts respectively. Multi-locus signatures consisted of 203, 36 and 178 CpGs and showed a large overlap with the corresponding single-locus signatures. The methylation signatures were associated with survival independently of tumour stage, age, AI treatment, and luminal status. The single-locus signature for the TAM cohort was conserved among the loci that were differentially methylated in endocrine-resistant T47D cells. Similarly, multi-locus signatures for the ER+/HER2- and AI cohorts were conserved in endocrine-resistant T47D cells. Also at the gene set level, several sets related to endocrine therapy and resistance were enriched in both survival and T47D signatures.

CONCLUSIONS

We identified individual and multivariable DNA methylation markers associated with therapy resistance independently of luminal status. Our results suggest that these markers identified from primary tumours prior to endocrine treatment are associated with development of endocrine resistance.

Prognostic significance of mutant-allele tumor heterogeneity in uterine corpus endometrial carcinoma

Background

Uterine corpus endometrial carcinoma (UCEC) is a clinically heterogeneous disease, and this heterogeneity is associated with tumor development, clinical characteristics, and prognostic outcomes. Mutant-allele tumor heterogeneity (MATH) is a novel, non-biased, quantitative measure to assess intra-tumor heterogeneity based on next-generation sequencing data. We aimed to explore the use of MATH as a measure for tumor heterogeneity and its prognostic role in UCEC patients.

Methods

We calculated MATH scores from the available data of 560 UCEC patients from The Cancer Genome Atlas (TCGA) and investigated their correlations with clinical characteristics, genetic alterations, and overall survival. Predictive accuracy was quantified using the area under the receiver operating characteristic curve (AUC) and the index of concordance (C-index).

Results

In total, 242 MATH scores were obtained from the UCEC cohort. MATH scores were significantly related to age, race, cancer type, clinical stage, histological grade, molecular type, targeted molecular therapy, and hormonal therapy. Furthermore, the genomic pattern on the basis of MATH scores showed that mutation rates of TP53 (tumor protein p53) and ARID1A (AT-rich interaction domain 1A) were independently associated with MATH scores. Correlation analysis revealed a significantly positive association of MATH scores with the fraction of somatic copy number alteration (SCNA). Importantly, a high MATH score was significantly associated with shorter overall survival [hazard ratio (HR), 2.342; 95% confidence interval (CI), 1.110-4.942]. Multivariate Cox regression combined with stratified analysis revealed that the MATH score is an independent prognostic factor in UCEC patients under 60 years old, and predictive quantification showed the MATH score had an AUC of 0.756 and a C-index of 0.845.

Conclusions

Our results suggest that MATH, a practical and useful way to measure intra-tumor heterogeneity, may serve as a significant biomarker for the prognosis of patients with UCEC, enabling more accurate prediction of clinical outcomes.

Deconvolution of the Genomic and Epigenomic Interaction Landscape of Triple-Negative Breast Cancer

Triple-negative breast cancer (TNBC) is the most aggressive form of breast cancer. Emerging evidenced suggests that both genetics and epigenetic factors play a role in the pathogenesis of TNBC. However, oncogenic interactions and cooperation between genomic and epigenomic variation have not been characterized. The objective of this study was to deconvolute the genomic and epigenomic interaction landscape in TNBC using an integrative genomics approach, which integrates information on germline, somatic, epigenomic and gene expression variation. We hypothesized that TNBC originates from a complex interplay between genomic (both germline and somatic variation) and epigenomic variation. We further hypothesized that these complex arrays of interacting genomic and epigenomic factors affect entire molecular networks and signaling pathways which, in turn, drive TNBC. We addressed these hypotheses using germline variation from genome-wide association studies and somatic, epigenomic and gene expression variation from The Cancer Genome Atlas (TCGA). The investigation revealed signatures of functionally related genes containing germline, somatic and epigenetic variations. DNA methylation had an effect on gene expression. Network and pathway analysis revealed molecule networks and signaling pathways enriched for germline, somatic and epigenomic variation, among them: Role of BRCA1 in DNA Damage Response, Hereditary Breast Cancer Signaling, Molecular Mechanisms of Cancer, Estrogen-Dependent Breast Cancer, p53, MYC Mediated Apoptosis, and PTEN Signaling pathways. The investigation revealed that integrative genomics is a powerful approach for deconvoluting the genomic-epigenomic interaction landscape in TNBC. Further studies are needed to understand the biological mechanisms underlying oncogenic interactions between genomic and epigenomic factors in TNBC.

Integrative epigenetic and genetic pan-cancer somatic alteration portraits

Genetic and epigenetic alterations are required for carcinogenesis and the mutation burden across tumor types has been investigated. Here, we investigate epigenetic alterations with a novel measure of global DNA methylation dysregulation, the methylation dysregulation index (MDI), across 14 cancer types in The Cancer Genome Atlas (TCGA) database. DNA methylation data-obtained using Illumina HumanMethylation450 BeadChip-was accessed from TCGA. We calculated the MDI in 14 tumor types (n = 5,592 tumors), using adjacent normal tissues (n = 701) from each tumor site. Copy number alteration, and mutation burden were retrieved from cBioportal (n = 5,152). We tested the relation of subject MDI across tumors and with age, gender, tumor stage, estimated tumor purity, and copy number alterations for both overall MDI and genomic-context-specific MDI. We also investigated the top most dysregulated loci shared across tumor types. There was a broad range of extent in methylation dysregulation across tumor types (P < 2.2E-16). However, a consistent pattern of methylation dysregulation stratified by genomic context was observed across tumor types where the highest dysregulation occurred at non-CpG island regions. Considering other summary measures of somatic alteration, MDI was correlated with copy number alterations but not with mutation burden. Using the top dysregulated CpG sites in common across tumors, 4 classes of cancer types were observed, and the functional consequences of these alterations to gene expression were confirmed. This work identified the global DNA methylation dysregulation patterns across 14 cancer types showing a higher impact for the non-CpG island areas. The most dysregulated loci across cancer types identified common clusters across cancer types that may have implications for future treatment and prevention measures.