Search

Simple Summary Breast cancer (BrCa) is characterized by aberrant DNA methylation. We leveraged high-throughput methylation data from BrCa and normal breast tissues and identified 11,176 to 27,786 differentially methylated genes (DMGs) against clinically relevant end-points. Innovative automated machine learning was employed to construct three highly performing signatures for (1) the discrimination of BrCa patients from healthy individuals, (2) the identification of BrCa metastatic disease and (3) the early diagnosis of BrCa. Furthermore, functional analysis revealed that most genes selected in the signatures showed associations to BrCa, with regulation of transcription being the main biological process, the nucleus being the main cellular component and transcription factor activity and sequence-specific DNA binding being the main molecular functions. Overall, revisiting methylome datasets led to three high-performance signatures that are readily available for improving BrCa precision management and significant knowledge mining related to disease pathophysiology. Abstract DNA methylation plays an important role in breast cancer (BrCa) pathogenesis and could contribute to driving its personalized management. We performed a complete bioinformatic analysis in BrCa whole methylome datasets, analyzed using the Illumina methylation 450 bead-chip array. Differential methylation analysis vs. clinical end-points resulted in 11,176 to 27,786 differentially methylated genes (DMGs). Innovative automated machine learning (AutoML) was employed to construct signatures with translational value. Three highly performing and low-feature-number signatures were built: (1) A 5-gene signature discriminating BrCa patients from healthy individuals (area under the curve (AUC): 0.994 (0.982–1.000)). (2) A 3-gene signature identifying BrCa metastatic disease (AUC: 0.986 (0.921–1.000)). (3) Six equivalent 5-gene signatures diagnosing early disease (AUC: 0.973 (0.920–1.000)). Validation in independent patient groups verified performance. Bioinformatic tools for functional analysis and protein interaction prediction were also employed. All protein encoding features included in the signatures were associated with BrCa-related pathways. Functional analysis of DMGs highlighted the regulation of transcription as the main biological process, the nucleus as the main cellular component and transcription factor activity and sequence-specific DNA binding as the main molecular functions. Overall, three high-performance diagnostic/prognostic signatures were built and are readily available for improving BrCa precision management upon prospective clinical validation. Revisiting archived methylomes through novel bioinformatic approaches revealed significant clarifying knowledge for the contribution of gene methylation events in breast carcinogenesis.

Autotaxin (ATX), encoded by the ctonucleotide pyrophosphatase/phosphodiesterase 2 (ENPP2) gene, is a key enzyme in lysophosphatidic acid (LPA) synthesis. We have recently described ENPP2 methylation profiles in health and multiple malignancies and demonstrated correlation to its aberrant expression. Here we focus on breast cancer (BrCa), analyzing in silico publicly available BrCa methylome datasets, to identify differentially methylated CpGs (DMCs) and correlate them with expression. Numerous DMCs were identified between BrCa and healthy breast tissues in the gene body and promoter-associated regions (PA). PA DMCs were upregulated in BrCa tissues in relation to normal, in metastatic BrCa in relation to primary, and in stage I BrCa in relation to normal, and this was correlated to decreased mRNA expression. The first exon DMC was also investigated in circulating cell free DNA (ccfDNA) isolated by BrCa patients; methylation was increased in BrCa in relation to ccfDNA from healthy individuals, confirming in silico results. It also differed between patient groups and was correlated to the presence of multiple metastatic sites. Our data indicate that promoter methylation of ENPP2 arrests its transcription in BrCa and introduce first exon methylation as a putative biomarker for diagnosis and monitoring which can be assessed in liquid biopsy.

Tissue-specific gene methylation events are key to the pathogenesis of several diseases and can be utilized for diagnosis and monitoring. Here, we established an in silico pipeline to analyze high-throughput methylome datasets to identify specific methylation fingerprints in three pathological entities of major burden, i.e., breast cancer (BrCa), osteoarthritis (OA) and diabetes mellitus (DM). Differential methylation analysis was conducted to compare tissues/cells related to the pathology and different types of healthy tissues, revealing Differentially Methylated Genes (DMGs). Highly performing and low feature number biosignatures were built with automated machine learning, including: (1) a five-gene biosignature discriminating BrCa tissue from healthy tissues (AUC 0.987 and precision 0.987), (2) three equivalent OA cartilage-specific biosignatures containing four genes each (AUC 0.978 and precision 0.986) and (3) a four-gene pancreatic β-cell-specific biosignature (AUC 0.984 and precision 0.995). Next, the BrCa biosignature was validated using an independent ccfDNA dataset showing an AUC and precision of 1.000, verifying the biosignature’s applicability in liquid biopsy. Functional and protein interaction prediction analysis revealed that most DMGs identified are involved in pathways known to be related to the studied diseases or pointed to new ones. Overall, our data-driven approach contributes to the maximum exploitation of high-throughput methylome readings, helping to establish specific disease profiles to be applied in clinical practice and to understand human pathology.

The corticotropin-releasing factor (CRF) system has been strongly associated with gastrointestinal pathophysiology, including colorectal cancer (CRC). We previously showed that altered expression of CRF receptors (CRFRs) in the colon critically affects CRC progression and aggressiveness through regulation of colonic inflammation. Here, we aimed to assess the potential of CRFR methylation levels as putative biomarkers in CRC. In silico methylation analysis of CRF receptor 1 (CRFR1) and CRF receptor 2 (CRFR2) was performed using methylome data derived by CRC and Crohn’s disease (CD) tissues and CRC-derived circulating cell-free DNAs (ccfDNAs). In total, 32 and 33 differentially methylated sites of CpGs (DMCs) emerged in CRFR1 and CRFR2, respectively, between healthy and diseased tissues. The methylation patterns were verified in patient-derived ccfDNA samples by qMSP and associated with clinicopathological characteristics. An automated machine learning (AutoML) technology was applied to ccfDNA samples for classification analysis. In silico analysis revealed increased methylation of both CRFRs in CRC tissue and ccfDNA-derived datasets. CRFR1 hypermethylation was also noticed in gene body DMCs of CD patients. CRFR1 hypermethylation was further validated in CRC adjuvant-derived ccfDNA samples, whereas CRFR1 hypomethylation, observed in metastasis-derived ccfDNAs, was correlated to disease aggressiveness and adverse prognostic characteristics. AutoML analysis based on CRFRs methylation status revealed a three-feature high-performing biosignature for CRC diagnosis with an estimated AUC of 0.929. Monitoring of CRFRs methylation-based signature in CRC tissues and ccfDNAs may be of high diagnostic and prognostic significance in CRC.