A DNA methylation-based definition of biologically distinct breast cancer subtypes

Obesity-associated epigenetic alterations and the obesity-breast cancer axis

Both breast cancer and obesity can regulate epigenetic changes or be regulated by epigenetic changes. Due to the well-established link between obesity and an increased risk of developing breast cancer, understanding how obesity-mediated epigenetic changes affect breast cancer pathogenesis is critical. Researchers have described how obesity and breast cancer modulate the epigenome individually and synergistically. In this review, the epigenetic alterations that occur in obesity, including DNA methylation, histone, and chromatin modification, accelerated epigenetic age, carcinogenesis, metastasis, and tumor microenvironment modulation, are discussed. Delineating the relationship between obesity and epigenetic regulation is vital to furthering our understanding of breast cancer pathogenesis.

Precise subtyping reveals immune heterogeneity for hormone receptor-positive breast cancer

A significant proportion of breast cancer cases are characterized by hormone receptor positivity (HR+). Clinically, the heterogeneity of HR+ breast cancer leads to different therapeutic effects on endocrine. Therefore, definition of subgroups in HR+ breast cancer is important for effective treatment. Here, we have developed a CMBR method utilizing computational functional networks based on DNA methylation to identify conserved subgroups in HR+ breast cancer. Calculated by CMBR, HR+ breast cancer was divided into five subgroups, of which HR+/negative epidermal growth factor receptor-2 (Her2-) was divided into two subgroups, and HR+/positive epidermal growth factor receptor-2 (Her2+) was divided into three subgroups. These subgroups had heterogeneity in the immune microenvironment, tumor infiltrating lymphocyte patterns, somatic mutation patterns and drug sensitivity. Specifically, CMBR identified two subgroups with the "Hot" tumor phenotype. In addition, these conserved subgroups were broadly validated on external validation datasets. CMBR identified the molecular signature of HR+ breast cancer subgroups, providing valuable insights into personalized treatment strategies and management options.

Anticancer Potential of Flavonoids: An Overview with an Emphasis on Tangeretin

Natural compounds with pharmacological activity, flavonoids have been the subject of an exponential increase in studies in the field of scientific research focused on therapeutic purposes due to their bioactive properties, such as antioxidant, anti-inflammatory, anti-aging, antibacterial, antiviral, neuroprotective, radioprotective, and antitumor activities. The biological potential of flavonoids, added to their bioavailability, cost-effectiveness, and minimal side effects, direct them as promising cytotoxic anticancer compounds in the optimization of therapies and the search for new drugs in the treatment of cancer, since some extensively antineoplastic therapeutic approaches have become less effective due to tumor resistance to drugs commonly used in chemotherapy. In this review, we emphasize the antitumor properties of tangeretin, a flavonoid found in citrus fruits that has shown activity against some hallmarks of cancer in several types of cancerous cell lines, such as antiproliferative, apoptotic, anti-inflammatory, anti-metastatic, anti-angiogenic, antioxidant, regulatory expression of tumor-suppressor genes, and epigenetic modulation.

Deep multi-omics integration by learning correlation-maximizing representation identifies prognostically stratified cancer subtypes

Motivation

Molecular subtyping by integrative modeling of multi-omics and clinical data can help the identification of robust and clinically actionable disease subgroups; an essential step in developing precision medicine approaches.

Results

We developed a novel outcome-guided molecular subgrouping framework, called Deep Multi-Omics Integrative Subtyping by Maximizing Correlation (DeepMOIS-MC), for integrative learning from multi-omics data by maximizing correlation between all input -omics views. DeepMOIS-MC consists of two parts: clustering and classification. In the clustering part, the preprocessed high-dimensional multi-omics views are input into two-layer fully connected neural networks. The outputs of individual networks are subjected to Generalized Canonical Correlation Analysis loss to learn the shared representation. Next, the learned representation is filtered by a regression model to select features that are related to a covariate clinical variable, for example, a survival/outcome. The filtered features are used for clustering to determine the optimal cluster assignments. In the classification stage, the original feature matrix of one of the -omics view is scaled and discretized based on equal frequency binning, and then subjected to feature selection using RandomForest. Using these selected features, classification models (for example, XGBoost model) are built to predict the molecular subgroups that were identified at clustering stage. We applied DeepMOIS-MC on lung and liver cancers, using TCGA datasets. In comparative analysis, we found that DeepMOIS-MC outperformed traditional approaches in patient stratification. Finally, we validated the robustness and generalizability of the classification models on independent datasets. We anticipate that the DeepMOIS-MC can be adopted to many multi-omics integrative analyses tasks.

Availability and implementation

Source codes for PyTorch implementation of DGCCA and other DeepMOIS-MC modules are available at GitHub (https://github.com/duttaprat/DeepMOIS-MC).

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

moBRCA-net: a breast cancer subtype classification framework based on multi-omics attention neural networks

BACKGROUND

Breast cancer is a highly heterogeneous disease that comprises multiple biological components. Owing its diversity, patients have different prognostic outcomes; hence, early diagnosis and accurate subtype prediction are critical for treatment. Standardized breast cancer subtyping systems, mainly based on single-omics datasets, have been developed to ensure proper treatment in a systematic manner. Recently, multi-omics data integration has attracted attention to provide a comprehensive view of patients but poses a challenge due to the high dimensionality. In recent years, deep learning-based approaches have been proposed, but they still present several limitations.

RESULTS

In this study, we describe moBRCA-net, an interpretable deep learning-based breast cancer subtype classification framework that uses multi-omics datasets. Three omics datasets comprising gene expression, DNA methylation and microRNA expression data were integrated while considering the biological relationships among them, and a self-attention module was applied to each omics dataset to capture the relative importance of each feature. The features were then transformed to new representations considering the respective learned importance, allowing moBRCA-net to predict the subtype.

CONCLUSIONS

Experimental results confirmed that moBRCA-net has a significantly enhanced performance compared with other methods, and the effectiveness of multi-omics integration and omics-level attention were identified. moBRCA-net is publicly available at https://github.com/cbi-bioinfo/moBRCA-net .

meth-SemiCancer: a cancer subtype classification framework via semi-supervised learning utilizing DNA methylation profiles

BACKGROUND

Identification of the cancer subtype plays a crucial role to provide an accurate diagnosis and proper treatment to improve the clinical outcomes of patients. Recent studies have shown that DNA methylation is one of the key factors for tumorigenesis and tumor growth, where the DNA methylation signatures have the potential to be utilized as cancer subtype-specific markers. However, due to the high dimensionality and the low number of DNA methylome cancer samples with the subtype information, still, to date, a cancer subtype classification method utilizing DNA methylome datasets has not been proposed.

RESULTS

In this paper, we present meth-SemiCancer, a semi-supervised cancer subtype classification framework based on DNA methylation profiles. The proposed model was first pre-trained based on the methylation datasets with the cancer subtype labels. After that, meth-SemiCancer generated the pseudo-subtypes for the cancer datasets without subtype information based on the model's prediction. Finally, fine-tuning was performed utilizing both the labeled and unlabeled datasets.

CONCLUSIONS

From the performance comparison with the standard machine learning-based classifiers, meth-SemiCancer achieved the highest average F1-score and Matthews correlation coefficient, outperforming other methods. Fine-tuning the model with the unlabeled patient samples by providing the proper pseudo-subtypes, encouraged meth-SemiCancer to generalize better than the supervised neural network-based subtype classification method. meth-SemiCancer is publicly available at https://github.com/cbi-bioinfo/meth-SemiCancer .

MIMESIS: minimal DNA-methylation signatures to quantify and classify tumor signals in tissue and cell-free DNA samples

DNA-methylation alterations are common in cancer and display unique characteristics that make them ideal markers for tumor quantification and classification. Here we present MIMESIS, a computational framework exploiting minimal DNA-methylation signatures composed by a few dozen informative DNA-methylation sites to quantify and classify tumor signals in tissue and cell-free DNA samples. Extensive analyses of multiple independent and heterogenous datasets including >7200 samples demonstrate the capability of MIMESIS to provide precise estimations of tumor content and to enable accurate classification of tumor type and molecular subtype. To assess our framework for clinical applications, we designed a MIMESIS-informed assay incorporating the minimal signatures for breast cancer. Using both artificial samples and clinical serial cell-free DNA samples from patients with metastatic breast cancer, we show that our approach provides accurate estimations of tumor content, sensitive detection of tumor signal and the ability to capture clinically relevant molecular subtype in patients' circulation. This study provides evidence that our extremely parsimonious approach can be used to develop cost-effective and highly scalable DNA-methylation assays that could support and facilitate the implementation of precision oncology in clinical practice.

Chromatin profile-based identification of a novel ER-positive breast cancer subgroup with reduced ER-responsive element accessibility

BACKGROUND

Oestrogen receptor (ER) signalling-dependent cancer cell growth is one of the major features of ER-positive breast cancer (BC). Inhibition of ER function is a standard and effective treatment for ER-positive tumours; however, ~20% of patients with ER-positive BC experience early or late recurrence. In this study, we examined intertumour heterogeneity from an epigenetic perspective based on the hypothesis that the intrinsic difference in epigenetic states around ER signalling pathway underlies endocrine therapy resistance.

METHODS

We performed transposase-accessible chromatin sequencing (ATAC-seq) analysis of 42 BC samples, including 35 ER-positive(+) human epidermal growth factor receptor 2 (HER2)-negative(-) and 7 triple-negative tumours. We also reanalysed ATAC-seq data of 45 ER + /HER2 - tumours in the Cancer Genome Atlas (TCGA) BC cohort to validate our observations.

RESULTS

We conducted a comprehensive analysis of cis-regulatory elements (CREs) using ATAC-seq, identifying three subgroups based on chromatin accessibility profiles. We identified a subgroup of ER-positive BCs with a distinctive chromatin accessibility pattern including reduced accessibility to ER-responsive elements (EREs). The same subgroup was also observed in TCGA BC cohort. Despite the reduced accessibility to EREs, the expression of ER and potential ER target genes were not decreased in these tumours.

CONCLUSION

Our findings highlight the existence of a subset of ER-positive BCs with unchanged ER expression but reduced EREs accessibility that cannot be distinguished by conventional immunostaining for ER. Future studies should determine whether these tumours are associated with resistance to endocrine therapy.

Genome-wide characterization of mitochondrial DNA methylation in human brain

Background

There is growing interest in the role of DNA methylation in regulating the transcription of mitochondrial genes, particularly in brain disorders characterized by mitochondrial dysfunction. Here, we present a novel approach to interrogate the mitochondrial DNA methylome at single base resolution using targeted bisulfite sequencing. We applied this method to investigate mitochondrial DNA methylation patterns in post-mortem superior temporal gyrus and cerebellum brain tissue from seven human donors.

Results

We show that mitochondrial DNA methylation patterns are relatively low but conserved, with peaks in DNA methylation at several sites, such as within the D-LOOP and the genes MT-ND2, MT-ATP6, MT-ND4, MT-ND5 and MT-ND6, predominantly in a non-CpG context. The elevated DNA methylation we observe in the D-LOOP we validate using pyrosequencing. We identify loci that show differential DNA methylation patterns associated with age, sex and brain region. Finally, we replicate previously reported differentially methylated regions between brain regions from a methylated DNA immunoprecipitation sequencing study.

Conclusions

We have annotated patterns of DNA methylation at single base resolution across the mitochondrial genome in human brain samples. Looking to the future this approach could be utilized to investigate the role of mitochondrial epigenetic mechanisms in disorders that display mitochondrial dysfunction.

Somatic Super-Enhancer Epigenetic Signature for Overall Survival Prediction in Patients with Breast Invasive Carcinoma

To analyze genome-wide super-enhancers (SEs) methylation signature of breast invasive carcinoma (BRCA) and its clinical value. Differential methylation sites (DMS) between BRCA and adjacent tissues from The Cancer Genome Atlas (TCGA) database were identified by using ChAMP package in R software. Super-enhancers were identified sing ROSE software. Overlap analysis was used to assess the potential DMS in SEs region. Feature selection was performed by Cox regression and least absolute shrinkage and selection operator (LASSO) algorithm based on TCGA training cohort. Prognosis model validation was performed in TCGA training cohort, TCGA validation cohort, and gene expression omnibus (GEO) test cohort. The gene ontology and KEGG analysis revealed that SEs target genes were significantly enriched in cell-migration-associated processes and pathways. A total of 83 654 DMS were identified between BRCA and adjacent tissues. Around 2397 DMS in SEs region were identified by overlap study and used to feature selection. By using Cox regression and LASSO algorithm, 42 features were selected to develop a clinical prediction model (CPM). Both training (TCGA) and validation cohorts (TCGA and GEO) show that the CPM has ideal discrimination and calibration. The CPM based on DMS at SE regions has ideal discrimination and calibration, which combined with tumor node metastasis (TNM) stage could improve prognostication, and thus contribute to individualized medicine.

Genome-wide differential DNA methylation analysis of MDA-MB-231 breast cancer cells treated with curcumin derivatives, ST08 and ST09

ST08 and ST09 are potent curcumin derivatives with antiproliferative, apoptotic, and migrastatic properties. Both ST08 and ST09 exhibit in vitro and in vivo anticancer properties. As reported earlier, these derivatives were highly cytotoxic towards MDA-MB-231 triple-negative breast cancer cells with IC50 values in the nanomolar (40-80nM) range.In this study,we performed whole-genome bisulfite sequencing(WGBS) of untreated (control), ST08 and ST09 (treated) triple-negative breast cancer cell line MDA-MB-231 to unravel epigenetic changes induced by the drug. We identified differentially methylated sites (DMSs) enriched in promoter regions across the genome. Analysis of the CpG island promoter methylation identified 12 genes common to both drugs, and 50% of them are known to be methylated in patient samples that were hypomethylated by drugs belonging to the homeobox family transcription factors.Methylation analysis of the gene body revealed 910 and 952 genes to be hypermethylatedin ST08 and ST09 treated MDA-MB-231 cells respectively. Correlation of the gene body hypermethylation with expression revealed CACNAH1 to be upregulated in ST08 treatment and CDH23 upregulation in ST09.Further, integrated analysis of the WGBS with RNA-seq identified uniquely altered pathways - ST08 altered ECM pathway, and ST09 cell cycle, indicating drug-specific signatures.

Classification of Subgroups with Immune Characteristics Based on DNA Methylation in Luminal Breast Cancer

Luminal breast cancer (BC) accounts for a large proportion of patients in BC, with high heterogeneity. Determining the precise subtype and optimal selection of treatment options for luminal BC is a challenge. In this study, we proposed an MSBR framework that integrate DNA methylation profiles and transcriptomes to identify immune subgroups of luminal BC. MSBR was implemented both on a key module scoring algorithm and “Boruta” feature selection method by DNA methylation. Luminal A was divided into two subgroups and luminal B was divided into three subgroups using the MSBR. Furthermore, these subgroups were defined as different immune subgroups in luminal A and B respectively. The subgroups showed significant differences in DNA methylation levels, immune microenvironment (immune cell infiltration, immune checkpoint PD1/PD-L1 expression, immune cell cracking activity (CYT)) and pathology features (texture, eccentricity, intensity and tumor-infiltrating lymphocytes (TILs)). The results also showed that there is a subgroup in both luminal A and B that has the benefit from immunotherapy. This study proposed a classification of luminal BC from the perspective of epigenetics and immune characteristics, which provided individualized treatment decisions.

Maternal obesity alters methylation level of cytosine in CpG island for epigenetic inheritance in fetal umbilical cord blood

BACKGROUND

Over the past few decades, global maternal obesity prevalence has rapidly increased. This condition may induce long-lasting pathophysiological effects on either fetal or infant health that could be attributable to unknown unique changes in the umbilical blood composition.

METHODS

A total of 34 overweight/obese and 32 normal-weight pregnant women were recruited. Fifteen umbilical blood samples including 8 overweight/obese subjects and 7 normal weight women were sequenced using Targeted Bisulfite Sequencing technology to detect the average methylation level of cytosine and identify the differentially methylated region (DMR). GO and KEGG analyses were then employed to perform pathway enrichment analysis of DMR-related genes and promoters. Moreover, the mRNA levels of methylation-related genes histone deacetylases (HDACs) and DNA methyltransferases (DNMTs) were characterized in the samples obtained from these two groups.

RESULTS

Average methylated cytosine levels in both the CpG islands (CGI) and promoter significantly decreased in overweight/obese groups. A total of 1669 DMRs exhibited differences in their DNA methylation status between the overweight/obese and control groups. GO and KEGG analyses revealed that DMR-related genes and promoters were enriched in the metabolism, cancer and cardiomyopathy signaling pathways. Furthermore, the HDACs and DNMTs mRNA levels trended to decline in overweight/obese groups.

CONCLUSIONS

Decreased methylated cytosine levels in overweight/obese women induce the gene expression activity at a higher level than in the control group. DMRs between these two groups in the fetal blood may contribute to the changes in gene transcription that underlie the increased risk of metabolic disorders, cancers and cardiomyopathy in their offspring.

DNA methylation-based epigenetic signatures predict somatic genomic alterations in gliomas

Molecular classification has improved diagnosis and treatment for patients with malignant gliomas. However, classification has relied on individual assays that are both costly and slow, leading to frequent delays in treatment. Here, we propose the use of DNA methylation, as an emerging clinical diagnostic platform, to classify gliomas based on major genomic alterations and provide insight into subtype characteristics. We show that using machine learning models, DNA methylation signatures can accurately predict somatic alterations and show improvement over existing classifiers. The established Unified Diagnostic Pipeline (UniD) we develop is rapid and cost-effective for genomic alterations and gene expression subtypes diagnostic at early clinical phase and improves over individual assays currently in clinical use. The significant relationship between genetic alteration and epigenetic signature indicates broad applicability of our approach to other malignancies.

No clinical assay currently exists to classify glioma tumours based on gene expression. Here, the authors develop a DNA methylation-based classifier, Unified Diagnostic Pipeline (UniD) that identifies genomic alterations and gene expression subtypes of gliomas.

DNA methylation changes in response to neoadjuvant chemotherapy are associated with breast cancer survival

BACKGROUND

Locally advanced breast cancer is a heterogeneous disease with respect to response to neoadjuvant chemotherapy (NACT) and survival. It is currently not possible to accurately predict who will benefit from the specific types of NACT. DNA methylation is an epigenetic mechanism known to play an important role in regulating gene expression and may serve as a biomarker for treatment response and survival. We investigated the potential role of DNA methylation as a prognostic marker for long-term survival (> 5 years) after NACT in breast cancer.

METHODS

DNA methylation profiles of pre-treatment (n = 55) and post-treatment (n = 75) biopsies from 83 women with locally advanced breast cancer were investigated using the Illumina HumanMethylation450 BeadChip. The patients received neoadjuvant treatment with epirubicin and/or paclitaxel. Linear mixed models were used to associate DNA methylation to treatment response and survival based on clinical response to NACT (partial response or stable disease) and 5-year survival, respectively. LASSO regression was performed to identify a risk score based on the statistically significant methylation sites and Kaplan-Meier curve analysis was used to estimate survival probabilities using ten years of survival follow-up data. The risk score developed in our discovery cohort was validated in an independent validation cohort consisting of paired pre-treatment and post-treatment biopsies from 85 women with locally advanced breast cancer. Patients included in the validation cohort were treated with either doxorubicin or 5-FU and mitomycin NACT.

RESULTS

DNA methylation patterns changed from before to after NACT in 5-year survivors, while no significant changes were observed in non-survivors or related to treatment response. DNA methylation changes included an overall loss of methylation at CpG islands and gain of methylation in non-CpG islands, and these changes affected genes linked to transcription factor activity, cell adhesion and immune functions. A risk score was developed based on four methylation sites which successfully predicted long-term survival in our cohort (p = 0.0034) and in an independent validation cohort (p = 0.049).

CONCLUSION

Our results demonstrate that DNA methylation patterns in breast tumors change in response to NACT. These changes in DNA methylation show potential as prognostic biomarkers for breast cancer survival.

CancerNet: a unified deep learning network for pan-cancer diagnostics

Background

Despite remarkable advances in cancer research, cancer remains one of the leading causes of death worldwide. Early detection of cancer and localization of the tissue of its origin are key to effective treatment. Here, we leverage technological advances in machine learning or artificial intelligence to design a novel framework for cancer diagnostics. Our proposed framework detects cancers and their tissues of origin using a unified model of cancers encompassing 33 cancers represented in The Cancer Genome Atlas (TCGA). Our model exploits the learned features of different cancers reflected in the respective dysregulated epigenomes, which arise early in carcinogenesis and differ remarkably between different cancer types or subtypes, thus holding a great promise in early cancer detection.

Results

Our comprehensive assessment of the proposed model on the 33 different tissues of origin demonstrates its ability to detect and classify cancers to a high accuracy (> 99% overall F-measure). Furthermore, our model distinguishes cancers from pre-cancerous lesions to metastatic tumors and discriminates between hypomethylation changes due to age related epigenetic drift and true cancer.

Conclusions

Beyond detection of primary cancers, our proposed computational model also robustly detects tissues of origin of secondary cancers, including metastatic cancers, second primary cancers, and cancers of unknown primaries. Our assessment revealed the ability of this model to characterize pre-cancer samples, a significant step forward in early cancer detection. Deployed broadly this model can deliver accurate diagnosis for a greatly expanded target patient population.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04783-y.

Development of an automated liquid biopsy assay for methylated markers in advanced breast cancer

Current molecular liquid biopsy assays to detect recurrence or monitor response to treatment require sophisticated technology, highly trained personnel, and a turnaround time of weeks. We describe the development and technical validation of an automated Liquid Biopsy for Breast Cancer Methylation (LBx-BCM) prototype, a DNA methylation detection cartridge assay that is simple to perform and quantitatively detects nine methylated markers within 4.5 h. LBx-BCM demonstrated high interassay reproducibility when analyzing exogenous methylated DNA (75-300 DNA copies) spiked into plasma (Coefficient of Variation, CV = 7.1 - 10.9%) and serum (CV = 19.1 - 36.1%). It also demonstrated high interuser reproducibility (Spearman r = 0.887, P < 0.0001) when samples of metastatic breast cancer (MBC, N = 11) and normal control (N = 4) were evaluated independently by two users. Analyses of interplatform reproducibility indicated very high concordance between LBx-BCM and the reference assay, cMethDNA, among 66 paired plasma samples (MBC N = 40, controls N = 26; Spearman r = 0.891; 95% CI = 0.825 - 0.933, P< 0.0001). LBx-BCM achieved a ROC AUC = 0.909 (95% CI = 0.836 - 0.982), 83% sensitivity and 92% specificity; cMethDNA achieved a ROC AUC = 0.896 (95% CI = 0.817 - 0.974), 83% sensitivity and 92% specificity in test set samples. The automated LBx-BCM cartridge prototype is fast, with performance levels equivalent to the highly sensitive, manual cMethDNA method. Future prospective clinical studies will evaluate LBx-BCM detection sensitivity and its ability to monitor therapeutic response during treatment for advanced breast cancer.

hTERT DNA Methylation Analysis Identifies a Biomarker for Retinoic Acid-Induced hTERT Repression in Breast Cancer Cell Lines

Telomerase reactivation is responsible for telomere preservation in about 90% of cancers, providing cancer cells an indefinite proliferating potential. Telomerase consists of at least two main subunits: a catalytic reverse transcriptase protein (hTERT) and an RNA template subunit. Strategies to inhibit hTERT expression seem promising for cancer treatment. Previous works showed that all-trans retinoic acid (ATRA) induces hTERT repression in acute promyelocytic leukemia cells, resulting in their death. Here, we investigated the effects of ATRA in a subset of breast cancer cell lines. The mutational status of hTERT promoter and the methylation patterns at a single CpG resolution were assessed. We observed an inverse relationship between hTERT expression after ATRA treatment and the methylation level of a specific CpG at chr5: 1,300,438 in a region of hTERT gene at −5 kb of the transcription initiation site. This observation highlighted the significance of this region, whose methylation profile could represent a promising biomarker to predict the sensitivity to ATRA-induced hTERT repression in specific breast cancer subtypes. As hTERT repression promotes drug-induced cell death, checking the methylation status of this unique region and the specific CpG included can help in decision-making to include ATRA in combination therapy and contributes to a better clinical outcome.

Epigenetic Inheritance From Normal Origin Cells Can Determine the Aggressive Biology of Tumor-Initiating Cells and Tumor Heterogeneity

The acquisition of genetic- and epigenetic-abnormalities during transformation has been recognized as the two fundamental factors that lead to tumorigenesis and determine the aggressive biology of tumor cells. However, there is a regularity that tumors derived from less-differentiated normal origin cells (NOCs) usually have a higher risk of vascular involvement, lymphatic and distant metastasis, which can be observed in both lymphohematopoietic malignancies and somatic cancers. Obviously, the hypothesis of genetic- and epigenetic-abnormalities is not sufficient to explain how the linear relationship between the cellular origin and the biological behavior of tumors is formed, because the cell origin of tumor is an independent factor related to tumor biology. In a given system, tumors can originate from multiple cell types, and tumor-initiating cells (TICs) can be mapped to different differentiation hierarchies of normal stem cells, suggesting that the heterogeneity of the origin of TICs is not completely chaotic. TIC’s epigenome includes not only genetic- and epigenetic-abnormalities, but also established epigenetic status of genes inherited from NOCs. In reviewing previous studies, we found much evidence supporting that the status of many tumor-related “epigenetic abnormalities” in TICs is consistent with that of the corresponding NOC of the same differentiation hierarchy, suggesting that they may not be true epigenetic abnormalities. So, we speculate that the established statuses of genes that control NOC’s migration, adhesion and colonization capabilities, cell-cycle quiescence, expression of drug transporters, induction of mesenchymal formation, overexpression of telomerase, and preference for glycolysis can be inherited to TICs through epigenetic memory and be manifested as their aggressive biology. TICs of different origins can maintain different degrees of innate stemness from NOC, which may explain why malignancies with stem cell phenotypes are usually more aggressive.

A Novel Gene Prognostic Signature Based on Differential DNA Methylation in Breast Cancer

Background: DNA methylation played essential roles in regulating gene expression. The impact of DNA methylation status on the occurrence and development of cancers has been well demonstrated. However, little is known about its prognostic role in breast cancer (BC).

Materials: The Illumina Human Methylation450 array (450k array) data of BC was downloaded from the UCSC xena database. Transcriptomic data of BC was downloaded from the Cancer Genome Atlas (TCGA) database. Firstly, we used univariate and multivariate Cox regression analysis to screen out independent prognostic CpGs, and then we identified methylation-associated prognosis subgroups by consensus clustering. Next, a methylation prognostic model was developed using multivariate Cox analysis and was validated with the Illumina Human Methylation27 array (27k array) dataset of BC. We then screened out differentially expressed genes (DEGs) between methylation high-risk and low-risk groups and constructed a methylation-based gene prognostic signature. Further, we validated the gene signature with three subgroups of the TCGA-BRCA dataset and an external dataset GSE146558 from the Gene Expression Omnibus (GEO) database.

Results: We established a methylation prognostic signature and a methylation-based gene prognostic signature, and there was a close positive correlation between them. The gene prognostic signature involved six genes: IRF2, KCNJ11, ZDHHC9, LRP11, PCMT1, and TMEM70. We verified their expression in mRNA and protein levels in BC. Both methylation and methylation-based gene prognostic signatures showed good prognostic stratification ability. The AUC values of 3-years, 5-years overall survival (OS) were 0.737, 0.744 in the methylation signature and 0.725, 0.715 in the gene signature, respectively. In the validation groups, high-risk patients were confirmed to have poorer OS. The AUC values of 3 years were 0.757, 0.735, 0.733 in the three subgroups of TCGA dataset and 0.635 in GSE146558 dataset.

Conclusion: This study revealed the DNA methylation landscape and established promising methylation and methylation-based gene prognostic signatures that could serve as potential prognostic biomarkers and therapeutic targets.

Targeted Bisulfite Sequencing Reveals DNA Methylation Changes in Zinc Finger Family Genes Associated With KRAS Mutated Colorectal Cancer

Background: Colorectal cancer (CRC) is a leading cause of cancer death, and early diagnosis of CRC could significantly reduce its mortality rate. Previous studies suggest that the DNA methylation status of zinc finger genes (ZFGs) could be of potential in CRC early diagnosis. However, the comprehensive evaluation of ZFGs in CRC is still lacking. Methods: We first collected 1,426 public samples on genome-wide DNA methylation, including 1,104 cases of CRC tumors, 54 adenomas, and 268 para-tumors. Next, the most differentially methylated ZFGs were identified and validated in two replication cohorts comprising 218 CRC patients. Finally, we compared the prediction capabilities between the ZFGs and the SEPT9 in all CRC patients and the KRAS + and KRAS- subgroup. Results: Five candidate ZFGs were selected: ESR1, ZNF132, ZNF229, ZNF542, and ZNF677. In particular, ESR1 [area under the curve (AUC) = 0.91] and ZNF132 (AUC = 0.93) showed equivalent or better diagnostic capability for CRC than SEPT9 (AUC = 0.91) in the validation dataset, suggesting that these two ZFGs might be of potential for CRC diagnosis in the future. Furthermore, we performed subgroup analysis and found a significantly higher diagnostic capability in KRAS + (AUC ranged from 0.97 to 1) than that in KRAS- patients (AUC ranged from 0.74 to 0.86) for all these five ZFGs, suggesting that these ZFGs could be ideal diagnostic markers for KRAS mutated CRC patients. Conclusion: The methylation profiles of the candidate ZFGs could be potential biomarkers for the early diagnosis of CRC, especially for patients carrying KRAS mutations.

Construction of Novel Methylation-Driven Gene Model and Investigation of PARVB Function in Glioblastoma

Background

Glioblastoma multiforme (GBM) is characterized by widespread genetic and transcriptional heterogeneity. Aberrant DNA methylation plays a vital role in GBM progression by regulating gene expression. However, little is known about the role of methylation and its association with prognosis in GBM. Our aim was to explore DNA methylation-driven genes (DMDGs) and provide evidence for survival prediction and individualized treatment of GBM patients.

Methods

Use of the MethylMix R package identified DMDGs in GBM. The prognostic signature of DMDGs based on the risk score was constructed by multivariate Cox regression analysis. Receiver operating characteristics (ROC) curve and C-index were applied to assess the predictive performance of the DMDG prognostic signature. The predictive ability of the multigene signature model was validated in TCGA and CGGA cohorts. Finally, the role of DMDG β-Parvin (PARVB) was explored in vitro.

Results

The prognostic signature of DMDGs was constructed based on six genes (MDK, NMNAT3, PDPN, PARVB, SERPINB1, and UPP1). The low-risk cohort had significantly better survival than the high-risk cohort (p < 0.001). The area under the curve of the ROC of the six-gene signature was 0.832, 0.927, and 0.980 within 1, 2, and 3 years, respectively. The C-index of 0.704 indicated superior specificity and sensitivity. The six-gene model has been demonstrated to be an independent prognostic factor for GBM. In addition, joint survival analysis indicated that the MDK, NMNAT3, PARVB, SERPINB1, and UPP1 genes were significantly associated with prognosis and therapeutic targets for GBM. Importantly, our DMDG prognostic model was more suitable and accurate for low-grade gliomas. Finally, we verified that PARVB induced epithelial-mesenchymal transition partially through the JAK2/STAT3 pathway, which in turn promoted GBM cell proliferation, migration, and invasion.

Conclusion

This study demonstrated the potential value of the prognostic signature of DMDGs and provided important bioinformatic and potential therapeutic target data to facilitate individualized treatment for GBM, and to elucidate the specific mechanism by which PARVB promotes GBM progression.

Gene-associated methylation status of ST14 as a predictor of survival and hormone receptor positivity in breast Cancer

Background

Genomic profiles of specific gene sets have been established to guide personalized treatment and prognosis for patients with breast cancer (BC). However, epigenomic information has not yet been applied in a clinical setting. ST14 encodes matriptase, a proteinase that is widely expressed in BC with reported prognostic value.

Methods

In this present study, we evaluated the effect of ST14 DNA methylation (DNAm) on overall survival (OS) of patients with BC as a representative example to promote the use of the epigenome in clinical decisions. We analyzed publicly available genomic and epigenomic data from 1361 BC patients. Methylation was characterized by the β-value from CpG probes based on sequencing with the Illumina Human 450 K platform.

Results

A high mean DNAm (β > 0.6779) across 34 CpG probes for ST14, as the gene-associated methylation (GAM) pattern, was associated with a longer OS after adjusting age, stage, histology and molecular features in Cox model (p value < 0.001). A high GAM status was also associated with a higher XBP1 expression level and higher proportion of hormone-positive BC (p value < 0.001). Pathway analysis revealed that altered GAM was related to matrisome-associated pathway.

Conclusions

Here we show the potential role of ST14 DNAm in BC prognosis and warrant further study.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-08645-3.

The Roles of DNA Demethylases in Triple-Negative Breast Cancer

Triple-negative breast cancers (TNBCs) are very heterogenous, molecularly diverse, and are characterized by a high propensity to relapse or metastasize. Clinically, TNBC remains a diagnosis of exclusion by the lack of hormone receptors (Estrogen Receptor (ER) and Progesterone Receptor (PR)) as well as the absence of overexpression and/or amplification of HER2. DNA methylation plays an important role in breast cancer carcinogenesis and TNBCs have a distinct DNA methylation profile characterized by marked hypomethylation and lower gains of methylations compared to all other subtypes. DNA methylation is regulated by the balance of DNA methylases (DNMTs) and DNA demethylases (TETs). Here, we review the roles of TETs as context-dependent tumor-suppressor genes and/or oncogenes in solid tumors, and we discuss the current understandings of the oncogenic role of TET1 and its therapeutic implications in TNBCs.

The epigenetics of breast cancer - Opportunities for diagnostics, risk stratification and therapy

The stage and molecular pathology-dependent prognosis of breast cancer, the limited treatment options for triple-negative carcinomas, as well as the development of resistance to therapies illustrate the need for improved early diagnosis and the development of new therapeutic approaches. Increasing data suggests that some answers to these challenges could be found in the area of epigenetics. In this study, we focus on the current research of the epigenetics of breast cancer, especially on the potential of epigenetics for clinical application in diagnostics, risk stratification and therapy. The differential DNA methylation status of specific gene regions has been used in the past to differentiate breast cancer cells from normal tissue. New technologies as detection of circulating nucleic acids including microRNAs to early detect breast cancer are emerging. Pattern of DNA methylation and expression of histone-modifying enzymes have been successfully used for risk stratification. However, all these epigenetic biomarkers should be validated in larger clinical studies. Recent preclinical and clinical studies show a therapeutic benefit of epigenetically active drugs for breast cancer entities that are still difficult to treat (triple negative, UICC stage IV). Remarkably, epigenetic therapies combined with chemotherapies or hormone-based therapies represent the most promising strategy. At the current stage, the integration of epigenetic substances into established breast cancer therapy protocols seems to hold the greatest potential for a clinical application of epigenetic research.

Integrated Analyses Reveal the Multi-Omics and Prognostic Characteristics of ATP5B in Breast Cancer

The beta subunit of F1Fo-ATP synthase (ATP5B) has been demonstrated to play an essential role in tumor progression and metastasis. However, there has been no comprehensive pan-cancer multi-omics analysis of ATP5B, while the clinical relevance of ATP5B and its potential mechanism in regulating breast cancer are still poorly understood. In this study, we demonstrated that ATP5B has a higher frequency of amplification than deletion in most cancer types, and the copy number variation (CNV) of ATP5B was significantly positively correlated with its mRNA expression level. DNA methylation analysis across pan-cancer also revealed a strong correlation between ATP5B expression and epigenetic changes. We identified 6 significant methylation sites involved in the regulation of ATP5B expression. Tissue microarrays (TMA) from 129 breast cancer samples, integrated with multiple additional breast cancer dataset, were used to evaluate the ATP5B expression and its correlation with prognosis. Higher levels of ATP5B expression were consistently associated with a worse OS in all datasets, and Cox regression analysis suggested that ATP5B expression was an independent prognostic factor. Gene enrichment analysis indicated that the gene signatures of DNA damage recognition, the E-cadherin nascent pathway and the PLK1 pathway were enriched in ATP5B-high patients. Moreover, somatic mutation analysis showed that a significant different mutation frequency of CDH1 and ADAMTSL3 could be observed between the ATP5B-high and ATP5B-low groups. In conclusion, this study reveals novel significance regarding the genetic characteristics and clinical value of ATP5B highlighted in predicting the outcome of breast cancer patients.

Accelerated aging in normal breast tissue of women with breast cancer

Background

DNA methylation alterations have similar patterns in normal aging tissue and in cancer. In this study, we investigated breast tissue-specific age-related DNA methylation alterations and used those methylation sites to identify individuals with outlier phenotypes. Outlier phenotype is identified by unsupervised anomaly detection algorithms and is defined by individuals who have normal tissue age-dependent DNA methylation levels that vary dramatically from the population mean.

Methods

We generated whole-genome DNA methylation profiles (GSE160233) on purified epithelial cells and used publicly available Infinium HumanMethylation 450K array datasets (TCGA, GSE88883, GSE69914, GSE101961, and GSE74214) for discovery and validation.

Results

We found that hypermethylation in normal breast tissue is the best predictor of hypermethylation in cancer. Using unsupervised anomaly detection approaches, we found that about 10% of the individuals (39/427) were outliers for DNA methylation from 6 DNA methylation datasets. We also found that there were significantly more outlier samples in normal-adjacent to cancer (24/139, 17.3%) than in normal samples (15/228, 5.2%). Additionally, we found significant differences between the predicted ages based on DNA methylation and the chronological ages among outliers and not-outliers. Additionally, we found that accelerated outliers (older predicted age) were more frequent in normal-adjacent to cancer (14/17, 82%) compared to normal samples from individuals without cancer (3/17, 18%). Furthermore, in matched samples, we found that the epigenome of the outliers in the pre-malignant tissue was as severely altered as in cancer.

Conclusions

A subset of patients with breast cancer has severely altered epigenomes which are characterized by accelerated aging in their normal-appearing tissue. In the future, these DNA methylation sites should be studied further such as in cell-free DNA to determine their potential use as biomarkers for early detection of malignant transformation and preventive intervention in breast cancer.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13058-021-01434-7.

Comprehensive analysis based on DNA methylation and RNA-seq reveals hypermethylation of the up-regulated WT1 gene with potential mechanisms in PAM50 subtypes of breast cancer

Background

Breast cancer (BC), one of the most widespread cancers worldwide, caused the deaths of more than 600,000 women in 2018, accounting for about 15% of all cancer-associated deaths in women that year. In this study, we aimed to discover potential prognostic biomarkers and explore their molecular mechanisms in different BC subtypes using DNA methylation and RNA-seq.

Methods

We downloaded the DNA methylation datasets and the RNA expression profiles of primary tissues of the four BC molecular subtypes (luminal A, luminal B, basal-like, and HER2-enriched), as well as the survival information from The Cancer Genome Atlas (TCGA). The highly expressed and hypermethylated genes across all the four subtypes were screened. We examined the methylation sites and the downstream co-expressed genes of the selected genes and validated their prognostic value using a different dataset (GSE20685). For selected transcription factors, the downstream genes were predicted based on the Gene Transcription Regulation Database (GTRD). The tumor microenvironment was also evaluated based on the TCGA dataset.

Results

We found that Wilms tumor gene 1 (WT1), a transcription factor, was highly expressed and hypermethylated in all the four BC subtypes. All the WT1 methylation sites exhibited hypermethylation. The methylation levels of the TSS200 and 1stExon regions were negatively correlated with WT1 expression in two BC subtypes, while that of the gene body region was positively associated with WT1 expression in three BC subtypes. Patients with low WT1 expression had better overall survival (OS). Five genes including COL11A1, GFAP, FGF5, CD300LG, and IGFL2 were predicted as the downstream genes of WT1. Those five genes were dysregulated in the four BC subtypes. Patients with a favorable 6-gene signature (low expression of WT1 and its five predicted downstream genes) exhibited better OS than that with an unfavorable 6-gene signature. We also found a correlation between WT1 and tamoxifen using STITCH. Higher infiltration rates of CD8 T cells, plasma cells, and monocytes were found in the lower quartile WT1 group and the favorable 6-gene signature group. In conclusion, we demonstrated that WT1 is hypermethylated and up-regulated in the four BC molecular subtypes and a 6-gene signature may predict BC prognosis.

Crosstalk between microRNA expression and DNA methylation drives the hormone-dependent phenotype of breast cancer

BACKGROUND

Abnormal DNA methylation is observed as an early event in breast carcinogenesis. However, how such alterations arise is still poorly understood. microRNAs (miRNAs) regulate gene expression at the post-transcriptional level and play key roles in various biological processes. Here, we integrate miRNA expression and DNA methylation at CpGs to study how miRNAs may affect the breast cancer methylome and how DNA methylation may regulate miRNA expression.

METHODS

miRNA expression and DNA methylation data from two breast cancer cohorts, Oslo2 (n = 297) and The Cancer Genome Atlas (n = 439), were integrated through a correlation approach that we term miRNA-methylation Quantitative Trait Loci (mimQTL) analysis. Hierarchical clustering was used to identify clusters of miRNAs and CpGs that were further characterized through analysis of mRNA/protein expression, clinicopathological features, in silico deconvolution, chromatin state and accessibility, transcription factor binding, and long-range interaction data.

RESULTS

Clustering of the significant mimQTLs identified distinct groups of miRNAs and CpGs that reflect important biological processes associated with breast cancer pathogenesis. Notably, two major miRNA clusters were related to immune or fibroblast infiltration, hence identifying miRNAs associated with cells of the tumor microenvironment, while another large cluster was related to estrogen receptor (ER) signaling. Studying the chromatin landscape surrounding CpGs associated with the estrogen signaling cluster, we found that miRNAs from this cluster are likely to be regulated through DNA methylation of enhancers bound by FOXA1, GATA2, and ER-alpha. Further, at the hub of the estrogen cluster, we identified hsa-miR-29c-5p as negatively correlated with the mRNA and protein expression of DNA methyltransferase DNMT3A, a key enzyme regulating DNA methylation. We found deregulation of hsa-miR-29c-5p already present in pre-invasive breast lesions and postulate that hsa-miR-29c-5p may trigger early event abnormal DNA methylation in ER-positive breast cancer.

CONCLUSIONS

We describe how miRNA expression and DNA methylation interact and associate with distinct breast cancer phenotypes.

Integrative analysis of DNA methylation and gene expression profiles identified potential breast cancer-specific diagnostic markers

Breast cancer is a common malignant tumor among women whose prognosis is largely determined by the period and accuracy of diagnosis. We here propose to identify a robust DNA methylation-based breast cancer-specific diagnostic signature. Genome-wide DNA methylation and gene expression profiles of breast cancer patients along with their adjacent normal tissues from the Cancer Genome Atlas (TCGA) were obtained as the training set. CpGs that with significantly elevated methylation level in breast cancer than not only their adjacent normal tissues and the other ten common cancers from TCGA but also the healthy breast tissues from the Gene Expression Omnibus (GEO) were finally remained for logistic regression analysis. Another independent breast cancer DNA methylation dataset from GEO was used as the testing set. Lots of CpGs were hyper-methylated in breast cancer samples compared with adjacent normal tissues, which tend to be negatively correlated with gene expressions. Eight CpGs located at RIIAD1, ENPP2, ESPN, and ETS1, were finally retained. The diagnostic model was reliable in separating BRCA from normal samples. Besides, chromatin accessibility status of RIIAD1, ENPP2, ESPN and ETS1 showed great differences between MCF-7 and MDA-MB-231 cell lines. In conclusion, the present study should be helpful for breast cancer early and accurate diagnosis.

Machine learning for deciphering cell heterogeneity and gene regulation

Epigenetics studies inheritable and reversible modifications of DNA that allow cells to control gene expression throughout their development and in response to environmental conditions. In computational epigenomics, machine learning is applied to study various epigenetic mechanisms genome wide. Its aim is to expand our understanding of cell differentiation, that is their specialization, in health and disease. Thus far, most efforts focus on understanding the functional encoding of the genome and on unraveling cell-type heterogeneity. Here, we provide an overview of state-of-the-art computational methods and their underlying statistical concepts, which range from matrix factorization and regularized linear regression to deep learning methods. We further show how the rise of single-cell technology leads to new computational challenges and creates opportunities to further our understanding of epigenetic regulation.

Genome-Wide Epigenetic Landscape of Lung Adenocarcinoma Links HOXB9 DNA Methylation to Intrinsic EGFR-TKI Resistance and Heterogeneous Responses

PURPOSE

Epidermal growth factor receptor (EGFR)-tyrosine kinase inhibitors (TKIs) show efficacy in treating patients with lung adenocarcinoma with EGFR-activating mutations. However, a significant subset of targeted patients fail to respond. Unlike acquired resistance (AR), intrinsic resistance (IR) remains poorly understood. We investigated whether epigenomic factors contribute to patient-to-patient heterogeneity in the EGFR-TKI response and aimed to characterize the IR subpopulation that obtains no benefit from EGFR-TKIs.

PATIENTS AND METHODS

We conducted genome-wide DNA methylation profiling of 79 tumors sampled from patients with advanced lung adenocarcinoma before they received EGFR-TKI treatment and analyzed the patient responses. Pyrosequencing was performed in a validation cohort of 163 patients with EGFR-activating mutations.

RESULTS

A DNA methylation landscape of 216 CpG sites with differential methylation was established to elucidate the association of DNA methylation with the characteristics and EGFR-TKI response status of the patients. Functional analysis of 37 transcription-repressive sites identified the enrichment of transcription factors, notably homeobox (HOX) genes. DNA methylation of HOXB9 (cg13643585) in the enhancer region yielded 88% sensitivity for predicting drug response (odds ratio [OR], 6.64; 95% CI, 1.98 to 25.23; P = .0009). Pyrosequencing validated that HOXB9 gained methylation in patients with a poor EGFR-TKI response (OR, 3.06; 95% CI, 1.13 to 8.19; P = .019).

CONCLUSION

Our data suggest that homeobox DNA methylation could be a novel tumor cellular state that can aid the precise categorization of tumor heterogeneity in the study of IR to EGFR-TKIs. We identified, for the first time, an epigenomic factor that can potentially complement DNA mutation status in discriminating patients with lung adenocarcinoma who are less likely to benefit from EGFR-TKI treatment, thereby leading to improved patient management in precision medicine.

Serum Long Non-Coding RNAs PVT1, HOTAIR, and NEAT1 as Potential Biomarkers in Egyptian Women with Breast Cancer

Long non-coding RNAs play an important role in tumor growth, angiogenesis, and metastasis in several types of cancer. However, the clinical significance of using lncRNAs as biomarkers for breast cancer diagnosis and prognosis is still poorly investigated. In this study, we analyzed the serum expression levels of lncRNAs PVT1, HOTAIR, NEAT1, and MALAT1, and their associated proteins, PAI-1, and OPN, in breast cancer patients compared to fibroadenoma patients and healthy subjects. Using quantitative real-time PCR (qRT-PCR), we compared the serum expression levels of the four circulating lncRNAs in patients with breast cancer (n = 50), fibroadenoma (n = 25), and healthy controls (n = 25). The serum levels of PAI-1 and OPN were measured using ELISA. Receiveroperating-characteristic (ROC) analysis and multivariate logistic regression were used to evaluate the diagnostic value of the selected parameters. The serum levels of HOTAIR, PAI-1, and OPN were significantly higher in breast cancer patients compared to controls and fibroadenoma patients. The serum level of PVT1 was significantly higher in breast cancer patients than in the controls, while that of NEAT1 was significantly lower in breast cancer patients compared to controls and fibroadenoma patients. Both ROC and multivariate logistic regression analyses revealed that PAI-1 has the greatest power in discriminating breast cancer from the control, whereas HOTAIR, PAI-1, and OPN have the greatest power in discriminating breast cancer from fibroadenoma patients. In conclusion, our data suggest that the serum levels of PVT1, HOTAIR, NEAT1, PAI-1, and OPN could serve as promising diagnostic biomarkers for breast cancer.

A Histone Acetylation Modulator Gene Signature for Classification and Prognosis of Breast Cancer

Regulators of histone acetylation are promising epigenetic targets for therapy in breast cancer. In this study, we comprehensively analyzed the expression of histone acetylation modulator genes in breast cancer using TCGA data sources. A gene signature composed of eight histone acetylation modulators (HAMs) was found to be effective for the classification and prognosis of breast cancers, especially in the HER2-enriched and basal-like molecular subtypes. The eight genes consist of two histone acetylation writers (GTF3C4 and CLOCK), two erasers (HDAC2 and SIRT7) and four readers (BRD4, BRD7, SP100, and BRWD3). Both histone acetylation writer genes and eraser genes were found to be differentially expressed between the two groups indicating a close relationship exists between overall histone acetylation level and prognosis of breast cancer in HER2-enriched and basal-like breast cancer.

Discovery of novel DNA methylation biomarkers for non-invasive sporadic breast cancer detection in the Latino population

Human diversity is one of the main pitfalls in the development of robust worldwide biomarkers in oncology. Epigenetic variability across human populations is associated with different genetic backgrounds, as well as variable lifestyles and environmental exposures, each of which should be investigated. To identify potential non-invasive biomarkers of sporadic breast cancer in the Uruguayan population, we studied genome-wide DNA methylation using Illumina methylation arrays in leukocytes of 22 women with sporadic breast cancer and 10 healthy women in a case-control study. We described a panel of 38 differentially methylated CpG positions that was able to cluster breast cancer patients (BCP) and controls, and that also recapitulated methylation differences in 12 primary breast tumors and their matched normal breast tissue. Moving forward, we simplified the detection method to improve its applicability in a clinical setting and used an independent well-characterized cohort of 80 leukocyte DNA samples from BCP and 80 healthy controls to validate methylation results at specific cancer-related genes. Our investigations identified methylation at CYFIP1 as a novel epigenetic biomarker candidate for sporadic breast cancer in the Uruguayan population. These results provide a proof-of-concept for the design of larger studies aimed at validating biomarker panels for the Latin American population.

Development and Validation of Epigenetic Signature Predict Survival for Patients with Laryngeal Squamous Cell Carcinoma

Establishing epigenetic signature to improve the accuracy of survival prediction and optimize therapeutic strategies for laryngeal squamous cell carcinoma (LSCC) by a genome-wide integrated analysis of methylation and the transcriptome. LSCC DNA methylation datasets and RNA sequencing datasets were acquired from the Cancer Genome Atlas (TCGA). MethylMix was applied to detect DNA methylation-driven genes (MDGs), which developed an epigenetic signature. The predictive accuracy and clinical value of the epigenetic signature were evaluated by receiver operating characteristic and decision curve analysis, and compared with tumor-node-metastasis (TNM) stage system. In addition, prognostic value of the epigenetic signature was validated by external Gene Expression Omnibus (GEO) database. According to five MDGs of epigenetic signature, the candidate small molecules for LSCC were screen out by the CMap database. A total of 88 DNA MDGs were identified, five of which (MAGEB2, SUSD1, ZNF382, ZNF418, and ZNF732) were chosen to construct an epigenetic signature. The epigenetic signature can effectively divide patients into high-risk and low-risk group, with the area under curve (AUC) of 0.8 (5-year overall survival [OS]) and AUC of 0.745 (3-year OS). Stratification analysis affirmed that the epigenetic signature was still a significant statistical prognostic model in subsets of patients with different clinical variables. Multivariate Cox regression analysis indicated that the efficacy of epigenetic signature appears independent of other clinicopathological characteristics. In terms of predictive capacity and clinical usefulness, the epigenetic signature was superior to traditional TNM stage. In addition, the epigenetic signature was confirmed in external LSCC cohorts from GEO. Finally, CMap matched the 10 most significant small molecules as promising therapeutic drugs to reverse the LSCC gene expression. An epigenetic signature, with five DNA MDGs, was identified and validated in LSCC patients by integrating multidimensional genomic data, which may offer novel research directions and prospects for individualized treatment of patients with LSCC.

Association ofSOD3 promoter DNA methylation with its down-regulation in breast carcinomas

Superoxide dismutase 3 (SOD3) is a secreted antioxidant enzyme that regulates reactive oxygen species in the microenvironment. It is also a potential tumour suppressor gene that is significantly downregulated in breast cancer. We have previously shown that its mRNA expression is inversely correlated with relapse free survival in breast cancer patients. This study aimed to investigate the correlation of SOD3 promoter DNA methylation with its expression in different molecular subtypes of breast carcinoma. We found that SOD3 expression was significantly reduced in breast carcinoma samples compared to normal tissues with the lowest levels observed in Luminal B subtype. Pyrosequencing analysis showed significant increase in methylation levels in the SOD3 promoter region (-108 and -19 from the TSS) in tumours vs normal tissues (53.6% vs 25.2%). The highest degree of correlation between methylation and SOD3 expression levels was observed in Luminal B subtype (Spearman's R = -0.540, P < 0.00093). In this subtype, the -78 CpG position is the most significantly methylated site. The Spearman's coefficient analysis also indicated the most significant correlation of DNA methylation at this site with SOD3 gene expression levels in tumours vs. normal tissues (R = -0.5816, P < 6.9E-12). Moreover, copy number variation analysis of TCGA database revealed that the more aggressive Triple Negative and Her2+ subtypes had higher levels of SOD3 gene deletion. The predominantly down-regulated expression pattern of SOD3 and the various genetic and epigenetic deregulations of its expression suggest that loss of this antioxidant promotes an advantageous tumour-promoting microenvironment in breast cancer.

Identification and validation of methylation-driven genes prognostic signature for recurrence of laryngeal squamous cell carcinoma by integrated bioinformatics analysis

Background

Recurrence remains a major obstacle to long-term survival of laryngeal squamous cell carcinoma (LSCC). We conducted a genome-wide integrated analysis of methylation and the transcriptome to establish methylation-driven genes prognostic signature (MDGPS) to precisely predict recurrence probability and optimize therapeutic strategies for LSCC.

Methods

LSCC DNA methylation datasets and RNA sequencing (RNA-seq) dataset were acquired from the Cancer Genome Atlas (TCGA). MethylMix was applied to detect DNA methylation-driven genes (MDGs). By univariate and multivariate Cox regression analyses, five genes of DNA MDGs was developed a recurrence-free survival (RFS)-related MDGPS. The predictive accuracy and clinical value of the MDGPS were evaluated by receiver operating characteristic (ROC) and decision curve analysis (DCA), and compared with TNM stage system. Additionally, prognostic value of MDGPS was validated by external Gene Expression Omnibus (GEO) database. According to 5 MDGs, the candidate small molecules for LSCC were screen out by the CMap database. To strengthen the bioinformatics analysis results, 30 pairs of clinical samples were evaluated by digoxigenin-labeled chromogenic in situ hybridization (CISH).

Results

A total of 88 DNA MDGs were identified, and five RFS-related MDGs (LINC01354, CCDC8, PHYHD1, MAGEB2 and ZNF732) were chosen to construct a MDGPS. The MDGPS can effectively divide patients into high-risk and low-risk group, with the area under curve (AUC) of 0.738 (5-year RFS) and AUC of 0.74 (3-year RFS). Stratification analysis affirmed that the MDGPS was still a significant statistical prognostic model in subsets of patients with different clinical variables. Multivariate Cox regression analysis indicated the efficacy of MDGPS appears independent of other clinicopathological characteristics. In terms of predictive capacity and clinical usefulness, the MDGPS was superior to traditional TNM stage. Additionally, the MDGPS was confirmed in external LSCC cohorts from GEO. CMap matched the 9 most significant small molecules as promising therapeutic drugs to reverse the LSCC gene expression. Finally, CISH analysis in 30 LSCC tissues and paired adjacent normal tissues revealed that MAGEB2 has significantly higher expression of LSCC compared to adjacent non-neoplastic tissues; LINC01354, CCDC8, PHYHD1, and ZNF732 have significantly lower expression of LSCC compared to adjacent non-neoplastic tissues, which were in line with bioinformatics analysis results.

Conclusion

A MDGPS, with five DNA MDGs, was identified and validated in LSCC patients by combining transcriptome and methylation datasets analysis. Compared TNM stage alone, it generates more accurate estimations of the recurrence prediction and maybe offer novel research directions and prospects for individualized treatment of patients with LSCC.

A comparative analysis of cell-type adjustment methods for epigenome-wide association studies based on simulated and real data sets

Technological advances and reduced costs of high-density methylation arrays have led to an increasing number of association studies on the possible relationship between human disease and epigenetic variability. DNA samples from peripheral blood or other tissue types are analyzed in epigenome-wide association studies (EWAS) to detect methylation differences related to a particular phenotype. Since information on the cell-type composition of the sample is generally not available and methylation profiles are cell-type specific, statistical methods have been developed for adjustment of cell-type heterogeneity in EWAS. In this study we systematically compared five popular adjustment methods: the factored spectrally transformed linear mixed model (FaST-LMM-EWASher), the sparse principal component analysis algorithm ReFACTor, surrogate variable analysis (SVA), independent SVA (ISVA) and an optimized version of SVA (SmartSVA). We used real data and applied a multilayered simulation framework to assess the type I error rate, the statistical power and the quality of estimated methylation differences according to major study characteristics. While all five adjustment methods improved false-positive rates compared with unadjusted analyses, FaST-LMM-EWASher resulted in the lowest type I error rate at the expense of low statistical power. SVA efficiently corrected for cell-type heterogeneity in EWAS up to 200 cases and 200 controls, but did not control type I error rates in larger studies. Results based on real data sets confirmed simulation findings with the strongest control of type I error rates by FaST-LMM-EWASher and SmartSVA. Overall, ReFACTor, ISVA and SmartSVA showed the best comparable statistical power, quality of estimated methylation differences and runtime.

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.

Unraveling the Genomic-Epigenomic Interaction Landscape in Triple Negative and Non-Triple Negative Breast Cancer

Background: The recent surge of next generation sequencing of breast cancer genomes has enabled development of comprehensive catalogues of somatic mutations and expanded the molecular classification of subtypes of breast cancer. However, somatic mutations and gene expression data have not been leveraged and integrated with epigenomic data to unravel the genomic-epigenomic interaction landscape of triple negative breast cancer (TNBC) and non-triple negative breast cancer (non-TNBC). Methods: We performed integrative data analysis combining somatic mutation, epigenomic and gene expression data from The Cancer Genome Atlas (TCGA) to unravel the possible oncogenic interactions between genomic and epigenomic variation in TNBC and non-TNBC. We hypothesized that within breast cancers, there are differences in somatic mutation, DNA methylation and gene expression signatures between TNBC and non-TNBC. We further hypothesized that genomic and epigenomic alterations affect gene regulatory networks and signaling pathways driving the two types of breast cancer. Results: The investigation revealed somatic mutated, epigenomic and gene expression signatures unique to TNBC and non-TNBC and signatures distinguishing the two types of breast cancer. In addition, the investigation revealed molecular networks and signaling pathways enriched for somatic mutations and epigenomic changes unique to each type of breast cancer. The most significant pathways for TNBC were: retinal biosynthesis, BAG2, LXR/RXR, EIF2 and P2Y purigenic receptor signaling pathways. The most significant pathways for non-TNBC were: UVB-induced MAPK, PCP, Apelin endothelial, Endoplasmatic reticulum stress and mechanisms of viral exit from host signaling Pathways. Conclusion: The investigation revealed integrated genomic, epigenomic and gene expression signatures and signing pathways unique to TNBC and non-TNBC, and a gene signature distinguishing the two types of breast cancer. The study demonstrates that integrative analysis of multi-omics data is a powerful approach for unravelling the genomic-epigenomic interaction landscape in TNBC and non-TNBC.

DNA methylation markers panel can improve prediction of response to neoadjuvant chemotherapy in luminal B breast cancer

Despite the advantages of neoadjuvant chemotherapy (NACT), associated toxicity is a serious complication that renders monitoring of the patients' response to NACT highly important. Thus, prediction of tumor response to treatment is imperative to avoid exposure of potential non-responders to deleterious complications. We have performed genome-wide analysis of DNA methylation by XmaI-RRBS and selected CpG dinucleotides differential methylation of which discriminates luminal B breast cancer samples with different sensitivity to NACT. With this data, we have developed multiplex methylation sensitive restriction enzyme PCR (MSRE-PCR) protocol for determining the methylation status of 10 genes (SLC9A3, C1QL2, DPYS, IRF4, ADCY8, KCNQ2, TERT, SYNDIG1, SKOR2 and GRIK1) that distinguish BC samples with different NACT response. Analysis of these 10 markers by MSRE-PCR in biopsy samples allowed us to reveal three top informative combinations of markers, (1) IRF4 and C1QL2; (2) IRF4, C1QL2, and ADCY8; (3) IRF4, C1QL2, and DPYS, with the areas under ROC curves (AUCs) of 0.75, 0.78 and 0.74, respectively. A classifier based on IRF4 and C1QL2 better meets the diagnostic panel simplicity requirements, as it consists of only two markers. Diagnostic accuracy of the panel of these two markers is 0.75, with the sensitivity of 75% and specificity of 75%.

DNMT1: A key drug target in triple-negative breast cancer

Triple-negative breast cancer (TNBC) is the most aggressive subtype of breast cancer. Altered epigenetics regulation including DNA hypermethylation by DNA methyltransferase 1 (DNMT1) has been implicated as one of the causes of TNBC tumorigenesis. In this review, the oncogenic functions rendered by DNMT1 in TNBCs, and DNMT1 inhibitors targeting TNBC cells are presented and discussed. In summary, DNMT1 expression is associated with poor breast cancer survival, and it is overexpressed in TNBC subtype. The oncogenic roles of DNMT1 in TNBCs include: (1) Repression of estrogen receptor (ER) expression; (2) Promotion of epithelial-mesenchymal transition (EMT) required for metastasis; (3) Induces cellular autophagy and; (4) Promotes the growth of cancer stem cells in TNBCs. DNMT1 confers these phenotypes by hypermethylating the promoter regions of ER, multiple tumor suppressor genes, microRNAs and epithelial markers involved in suppressing EMT. DNMT1 inhibitors exert anti-tumorigenic effects against TNBC cells. This includes the hypomethylating agents azacitidine, decitabine and guadecitabine that might sensitize TNBC patients to immune checkpoint blockade therapy. DNMT1 represents an epigenetic target for TNBC cells destruction as well as to derail their metastatic and aggressive phenotypes.

Estimating the magnitude of clinical benefit from (neo)adjuvant chemotherapy in patients with ER-positive/HER2-negative breast cancer

Gene-expression assays were originally validated retrospectively as tools of prognostication, with evidence emerging from more recent prospectively-conducted studies such as MINDACT and TAILORx supporting their clinical validity and utility as biomarkers in identifying patients with luminal breast cancer who might be spared chemotherapy. However, these assays still do not have the ability to identify all patients who may safely avoid chemotherapy, and may over-estimate the risk of relapse in some cases. Future studies should aim to prospectively integrate contemporary approaches that assume a theoretical risk of relapse (based on pathological and/or genomic evaluation of the primary tumour), with new tools that can detect signals of active micro-metastatic disease. Until current methods of estimating prognosis and predicting benefit from adjuvant chemotherapy are significantly refined, estimating and improving the true magnitude of benefit derived from chemotherapy remains a challenge for clinicians.

Aberrantly Methylated cfDNA in Body Fluids as a Promising Diagnostic Tool for Early Detection of Breast Cancer

Breast malignancies are the leading type of cancer among women. Its prevention and early detection, particularly in young women, remains challenging. To this end, cell-free DNA (cfDNA) detected in body fluids demonstrates great potential for early detection of tissue transformation and altered molecular setup, such as epigenetic profiles. Aberrantly methylated cfDNA in body fluids could therefore serve as a potential diagnostic and prognostic tool in breast cancer management. Abnormal methylation may lead to both an activation of oncogenes via hypomethylation and an inactivation of tumor suppressor genes by hypermethylation. We update the state of the art in the area of aberrant cfDNA methylation analyses as a diagnostic and prognostic tool in breast cancer, report on the main technological challenges, and provide an outlook for advancing the overall management of breast malignancies based on cfDNA as a target for diagnosis and tailored therapies.

Integrative cancer patient stratification via subspace merging

MOTIVATION

Technologies that generate high-throughput omics data are flourishing, creating enormous, publicly available repositories of multi-omics data. As many data repositories continue to grow, there is an urgent need for computational methods that can leverage these data to create comprehensive clusters of patients with a given disease.

RESULTS

Our proposed approach creates a patient-to-patient similarity graph for each data type as an intermediate representation of each omics data type and merges the graphs through subspace analysis on a Grassmann manifold. We hypothesize that this approach generates more informative clusters by preserving the complementary information from each level of omics data. We applied our approach to The Cancer Genome Atlas (TCGA) breast cancer dataset and show that by integrating gene expression, microRNA and DNA methylation data, our proposed method can produce clinically useful subtypes of breast cancer. We then investigate the molecular characteristics underlying these subtypes. We discover a highly expressed cluster of genes on chromosome 19p13 that strongly correlates with survival in TCGA breast cancer patients and validate these results in three additional breast cancer datasets. We also compare our approach with previous integrative clustering approaches and obtain comparable or superior results.

AVAILABILITY AND IMPLEMENTATION

https://github.com/michaelsharpnack/GrassmannCluster.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

5-Aza-2'-deoxycytidine advances the epithelial-mesenchymal transition of breast cancer cells by demethylating Sipa1 promoter-proximal elements

Human breast cancer cells exhibit considerable diversity in the methylation status of genomic DNA CpGs that regulate metastatic transcriptome networks. In this study, we identified human Sipa1 promoter-proximal elements that contained a CpG island and demonstrated that the methylation status of the CpG island was inversely correlated with SIPA1 protein expression in cancer cells. 5-Aza-2′-deoxycytidine (5-Aza-CdR), a DNA methyltransferase inhibitor, promoted the expression of Sipa1 in the MCF7 breast cancer cells with a low level of SIPA1 expression. On the contrary, in MDA-MB-231 breast cancer cells with high SIPA1 expression levels, hypermethylation of the CpG island negatively regulated the transcription of Sipa1. In addition, the epithelial–mesenchymal transition (EMT) was reversed after knocking down Sipa1 in MDA-MB-231 cells. However, the EMT was promoted in MCF7 cells with over-expression of SIPA1 or treated with 5-Aza-CdR. Taken together, hypomethylation of the CpG island in Sipa1 promoter-proximal elements could enhance SIPA1 expression in breast cancer cells, which could facilitate EMT of cancer cells, possibly increasing a risk of cancer cell metastasis in individuals treated with 5-Aza-CdR.

Summary: Hypomethylation by 5-Aza-CdR upregulates the SIPA1 expression and promotes epithelial–mesenchymal transition in breast cancer cells, possibly increasing the risk of cancer cell metastasis.

Intravenous 5-fluoro-2'-deoxycytidine administered with tetrahydrouridine increases the proportion of p16-expressing circulating tumor cells in patients with advanced solid tumors

PURPOSE

Following promising responses to the DNA methyltransferase (DNMT) inhibitor 5-fluoro-2'-deoxycytidine (FdCyd) combined with tetrahydrouridine (THU) in phase 1 testing, we initiated a non-randomized phase 2 study to assess response to this combination in patients with advanced solid tumor types for which tumor suppressor gene methylation is potentially prognostic. To obtain pharmacodynamic evidence for DNMT inhibition by FdCyd, we developed a novel method for detecting expression of tumor suppressor protein p16/INK4A in circulating tumor cells (CTCs).

METHODS

Patients in histology-specific strata (breast, head and neck [H&N], or non-small cell lung cancers [NSCLC] or urothelial transitional cell carcinoma) were administered FdCyd (100 mg/m2) and THU (350 mg/m2) intravenously 5 days/week for 2 weeks, in 28-day cycles, and progression-free survival (PFS) rate and objective response rate (ORR) were evaluated. Blood specimens were collected for CTC analysis.

RESULTS

Ninety-three eligible patients were enrolled (29 breast, 21 H&N, 25 NSCLC, and 18 urothelial). There were three partial responses. All strata were terminated early due to insufficient responses (H&N, NSCLC) or slow accrual (breast, urothelial). However, the preliminary 4-month PFS rate (42%) in the urothelial stratum exceeded the predefined goal-though the ORR (5.6%) did not. An increase in the proportion of p16-expressing cytokeratin-positive CTCs was detected in 69% of patients evaluable for clinical and CTC response, but was not significantly associated with clinical response.

CONCLUSION

Further study of FdCyd + THU is potentially warranted in urothelial carcinoma but not NSCLC or breast or H&N cancer. Increase in the proportion of p16-expressing cytokeratin-positive CTCs is a pharmacodynamic marker of FdCyd target engagement.

Targeting Chromatin Remodeling for Cancer Therapy

Background:

Epigenetic alterations comprise key regulatory events that dynamically alter gene expression and their deregulation is commonly linked to the pathogenesis of various diseases, including cancer. Unlike DNA mutations, epigenetic alterations involve modifications to proteins and nucleic acids that regulate chromatin structure without affecting the underlying DNA sequence, altering the accessibility of the transcriptional machinery to the DNA, thus modulating gene expression. In cancer cells, this often involves the silencing of tumor suppressor genes or the increased expression of genes involved in oncogenesis. Advances in laboratory medicine have made it possible to map critical epigenetic events, including histone modifications and DNA methylation, on a genome-wide scale. Like the identification of genetic mutations, mapping of changes to the epigenetic landscape has increased our understanding of cancer progression. However, in contrast to irreversible genetic mutations, epigenetic modifications are flexible and dynamic, thereby making them promising therapeutic targets. Ongoing studies are evaluating the use of epigenetic drugs in chemotherapy sensitization and immune system modulation. With the preclinical success of drugs that modify epigenetics, along with the FDA approval of epigenetic drugs including the DNA methyltransferase 1 (DNMT1) inhibitor 5-azacitidine and the histone deacetylase (HDAC) inhibitor vorinostat, there has been a rise in the number of drugs that target epigenetic modulators over recent years.

Conclusion:

We provide an overview of epigenetic modulations, particularly those involved in cancer, and discuss the recent advances in drug development that target these chromatin-modifying events, primarily focusing on novel strategies to regulate the epigenome.