Convergent alterations in the tumor microenvironment of MYC-driven human and murine prostate cancer.

The tissue microenvironment in prostate cancer is profoundly altered. While such alterations have been implicated in driving prostate cancer initiation and progression to aggressive disease, how prostate cancer cells and their precursors mediate those changes is unclear, in part due to the inability to longitudinally study the disease evolution in human tissues. To overcome this limitation, we performed extensive single-cell RNA-sequencing (scRNA-seq) and rigorous molecular pathology of the comparative biology between human prostate cancer and key time points in the disease evolution of a genetically engineered mouse model (GEMM) of prostate cancer. Our studies of human tissues, with validation in a large external data set, revealed that cancer cell-intrinsic activation of MYC signaling was the top up-regulated pathway in human cancers, representing a common denominator across the well-known molecular and pathological heterogeneity of human prostate cancer. Likewise, numerous non-malignant cell states in the tumor microenvironment (TME), including non-cancerous epithelial, immune, and fibroblast cell compartments, were conserved across individuals, raising the possibility that these cell types may be a sequelae of the convergent MYC activation in the cancer cells. To test this hypothesis, we employed a GEMM of prostate epithelial cell-specific MYC activation in two mouse strains. Cell communication network and pathway analyses suggested that MYC oncogene-expressing neoplastic cells, directly and indirectly, reprogrammed the TME during carcinogenesis, leading to the emergence of cascading cell state alterations in neighboring epithelial, immune, and fibroblast cell types that paralleled key findings in human prostate cancer. Importantly, among these changes, the progression from a precursor-enriched to invasive-cancer-enriched state was accompanied by a cell-intrinsic switch from pro-immunogenic to immunosuppressive transcriptional programs with coinciding enrichment of immunosuppressive myeloid and Treg cells in the immune microenvironment. These findings implicate activation of MYC signaling in reshaping convergent aspects of the TME of prostate cancer as a common denominator across the otherwise well-documented molecular heterogeneity of human prostate cancer.

Corrigendum: Convergent antibody responses are associated with broad neutralization of hepatitis C virus.

[This corrects the article DOI: 10.3389/fimmu.2023.1135841.].

Convergent antibody responses are associated with broad neutralization of hepatitis C virus.

Introduction: Early development of broadly neutralizing antibodies (bNAbs) targeting the hepatitis C virus (HCV) envelope glycoprotein E2 is associated with spontaneous clearance of infection, so induction of bNAbs is a major goal of HCV vaccine development. However, the molecular antibody features important for broad neutralization are not known. Methods: To identify B cell repertoire features associated with broad neutralization, we performed RNA sequencing of the B cell receptors (BCRs) of HCV E2-reactive B cells of HCV-infected individuals with either high or low plasma neutralizing breadth. We then produced a monoclonal antibody (mAb) expressed by pairing the most abundant heavy and light chains from public clonotypes identified among clearance, high neutralization subjects. Results: We found distinctive BCR features associated with broad neutralization of HCV, including long heavy chain complementarity determining region 3 (CDRH3) regions, specific VH gene usage, increased frequencies of somatic hypermutation, and particular VH gene mutations. Most intriguing, we identified many E2-reactive public BCR clonotypes (heavy and light chain clones with the same V and J-genes and identical CDR3 sequences) present only in subjects who produced highly neutralizing plasma. The majority of these public clonotypes were shared by two subjects who cleared infection. A mAb expressing the most abundant public heavy and light chains from these clearance, high neutralization subjects had features enriched in high neutralization clonotypes, such as increased somatic hypermutation frequency and usage of IGHV1-69, and was cross-neutralizing. Discussion: Together, these results demonstrate distinct BCR repertoires associated with high plasma neutralizing capacity. Further characterization of the molecular features and function of these antibodies can inform HCV vaccine development.

The endocannabinoid system promotes hepatocyte progenitor cell proliferation and maturation by modulating cellular energetics.

The proliferation and differentiation of hepatic progenitor cells (HPCs) drive the homeostatic renewal of the liver under diverse conditions. Liver regeneration is associated with an increase in Axin2+Cnr1+ HPCs, along with a marked increase in the levels of the endocannabinoid anandamide (AEA). But the molecular mechanism linking AEA signaling to HPC proliferation and/or differentiation has not been explored. Here, we show that in vitro exposure of HPCs to AEA triggers both cell cycling and differentiation along with increased expression of Cnr1, Krt19, and Axin2. Mechanistically, we found that AEA promotes the nuclear localization of the transcription factor ß-catenin, with subsequent induction of its downstream targets. Systemic analyses of cells after CRISPR-mediated knockout of the ß-catenin-regulated transcriptome revealed that AEA modulates ß-catenin-dependent cell cycling and differentiation, as well as interleukin pathways. Further, we found that AEA promotes OXPHOS in HPCs when amino acids and glucose are readily available as substrates, but AEA inhibits it when the cells rely primarily on fatty acid oxidation. Thus, the endocannabinoid system promotes hepatocyte renewal and maturation by stimulating the proliferation of Axin2+Cnr1+ HPCs via the ß-catenin pathways while modulating the metabolic activity of their precursor cells.

Transcriptional programs of neoantigen-specific TIL in anti-PD-1-treated lung cancers.

PD-1 blockade unleashes CD8 T cells1, including those specific for mutation-associated neoantigens (MANA), but factors in the tumour microenvironment can inhibit these T cell responses. Single-cell transcriptomics have revealed global T cell dysfunction programs in tumour-infiltrating lymphocytes (TIL). However, the majority of TIL do not recognize tumour antigens2, and little is known about transcriptional programs of MANA-specific TIL. Here, we identify MANA-specific T cell clones using the MANA functional expansion of specific T cells assay3 in neoadjuvant anti-PD-1-treated non-small cell lung cancers (NSCLC). We use their T cell receptors as a 'barcode' to track and analyse their transcriptional programs in the tumour microenvironment using coupled single-cell RNA sequencing and T cell receptor sequencing. We find both MANA- and virus-specific clones in TIL, regardless of response, and MANA-, influenza- and Epstein-Barr virus-specific TIL each have unique transcriptional programs. Despite exposure to cognate antigen, MANA-specific TIL express an incompletely activated cytolytic program. MANA-specific CD8 T cells have hallmark transcriptional programs of tissue-resident memory (TRM) cells, but low levels of interleukin-7 receptor (IL-7R) and are functionally less responsive to interleukin-7 (IL-7) compared with influenza-specific TRM cells. Compared with those from responding tumours, MANA-specific clones from non-responding tumours express T cell receptors with markedly lower ligand-dependent signalling, are largely confined to HOBIThigh TRM subsets, and coordinately upregulate checkpoints, killer inhibitory receptors and inhibitors of T cell activation. These findings provide important insights for overcoming resistance to PD-1 blockade.

Bile acid-activated macrophages promote biliary epithelial cell proliferation through integrin αvß6 upregulation following liver injury.

Cholangiopathies caused by biliary epithelial cell (BEC) injury represent a leading cause of liver failure. No effective pharmacologic therapies exist, and the underlying mechanisms remain obscure. We aimed to explore the mechanisms of bile duct repair after targeted BEC injury. Injection of intermedilysin into BEC-specific human CD59 (hCD59) transgenic mice induced acute and specific BEC death, representing a model to study the early signals that drive bile duct repair. Acute BEC injury induced cholestasis followed by CCR2+ monocyte recruitment and BEC proliferation. Using microdissection and next-generation RNA-Seq, we identified 5 genes, including Mapk8ip2, Cdkn1a, Itgb6, Rgs4, and Ccl2, that were most upregulated in proliferating BECs after acute injury. Immunohistochemical analyses confirmed robust upregulation of integrin αvß6 (ITGß6) expression in this BEC injury model, after bile duct ligation, and in patients with chronic cholangiopathies. Deletion of the Itgb6 gene attenuated BEC proliferation after acute bile duct injury. Macrophage depletion or Ccr2 deficiency impaired ITGß6 expression and BEC proliferation. In vitro experiments revealed that bile acid-activated monocytes promoted BEC proliferation through ITGß6. Our data suggest that BEC injury induces cholestasis, monocyte recruitment, and induction of ITGß6, which work together to promote BEC proliferation and therefore represent potential therapeutic targets for cholangiopathies.

Longitudinal gut microbiome changes in alcohol use disorder are influenced by abstinence and drinking quantity.

Many patients with alcohol use disorder (AUD) consume alcohol chronically and in large amounts that alter intestinal microbiota, damage the gastrointestinal tract, and thereby injure other organs via malabsorption and intestinal inflammation. We hypothesized that alcohol consumption and subsequent abstinence would change the gut microbiome in adults admitted to a treatment program. Stool and oral specimens, diet data, gastrointestinal assessment scores, anxiety, depression measures and drinking amounts were collected longitudinally for up to 4 weeks in 22 newly abstinent inpatients with AUD who were dichotomized as less heavy drinkers (LHD, <10 drinks/d) and very heavy drinkers (VHD, 10 or more drinks/d). Next-generation 16 S rRNA gene sequencing was performed to measure the gut and oral microbiome at up to ten time points/subject and LHD and VHD were compared for change in principal components, Shannon diversity index and specific genera. The first three principal components explained 46.7% of the variance in gut microbiome diversity across time and all study subjects, indicating the change in gut microbiome following abstinence. The first time point was an outlier in three-dimensional principal component space versus all other time points. The gut microbiota in LHD and VHD were significantly dissimilar in change from day 1 to day 5 (p = .03) and from day 1 to week 3 (p = .02). The VHD drinking group displayed greater change from baseline. The Shannon diversity index of the gut microbiome changed significantly during abstinence in five participants. In both groups, the Shannon diversity was lower in the oral microbiome than gut. Ten total genera were shared between oral and stool in the AUD participants. These data were compared with healthy controls from the Human Microbiome Project to investigate the concept of a core microbiome. Rapid changes in gut microbiome following abstinence from alcohol suggest resilience of the gut microbiome in AUD and reflects the benefits of refraining from the highest levels of alcohol and potential benefits of abstinence.

Ventral midbrain astrocytes display unique physiological features and sensitivity to dopamine D2 receptor signaling.

Astrocytes are ubiquitous CNS cells that support tissue homeostasis through ion buffering, neurotransmitter recycling, and regulation of CNS vasculature. Yet, despite the essential functional roles they fill, very little is known about the physiology of astrocytes in the ventral midbrain, a region that houses dopamine-releasing neurons and is critical for reward learning and motivated behaviors. Here, using a combination of whole-transcriptome sequencing, histology, slice electrophysiology, and calcium imaging, we performed the first functional and molecular profiling of ventral midbrain astrocytes and observed numerous differences between these cells and their telencephalic counterparts, both in their gene expression profile and in their physiological properties. Ventral midbrain astrocytes have very low membrane resistance and inward-rectifying potassium channel-mediated current, and are extensively coupled to surrounding oligodendrocytes through gap junctions. They exhibit calcium responses to glutamate but are relatively insensitive to norepinephrine. In addition, their calcium activity can be dynamically modulated by dopamine D2 receptor signaling. Taken together, these data indicate that ventral midbrain astrocytes are physiologically distinct from astrocytes in cortex and hippocampus. This work provides new insights into the extent of functional astrocyte heterogeneity within the adult brain and establishes the foundation for examining the impact of regional astrocyte differences on dopamine neuron function and susceptibility to degeneration.

Local Cues Establish and Maintain Region-Specific Phenotypes of Basal Ganglia Microglia.

Microglia play critical roles in tissue homeostasis and can also modulate neuronal function and synaptic connectivity. In contrast to astrocytes and oligodendrocytes, which arise from multiple progenitor pools, microglia arise from yolk sac progenitors and are widely considered to be equivalent throughout the CNS. However, little is known about basic properties of deep brain microglia, such as those within the basal ganglia (BG). Here, we show that microglial anatomical features, lysosome content, membrane properties, and transcriptomes differ significantly across BG nuclei. Region-specific phenotypes of BG microglia emerged during the second postnatal week and were re-established following genetic or pharmacological microglial ablation and repopulation in the adult, indicating that local cues play an ongoing role in shaping microglial diversity. These findings demonstrate that microglia in the healthy brain exhibit a spectrum of distinct functional states and provide a critical foundation for defining microglial contributions to BG circuit function.

A Four-Gene Promoter Methylation Marker Panel Consisting of GREM1, NEURL, LAD1, and NEFH Predicts Survival of Clear Cell Renal Cell Cancer Patients.

Purpose: The currently used prognostic models for patients with nonmetastatic clear cell renal cell carcinoma (ccRCC) are based on clinicopathologic features and might be improved by adding molecular markers. Epigenetic alterations occur frequently in ccRCC and are promising biomarkers. The aim of this study is to identify prognostic promoter methylation markers for ccRCC.Experimental Design: We integrated data generated by massive parallel sequencing of methyl-binding domain enriched DNA and microarray-based RNA expression profiling of 5-aza-2'-deoxycytidine-treated ccRCC cell lines to comprehensively characterize the ccRCC methylome. A selection of the identified methylation markers was evaluated in two independent series of primary ccRCC (n = 150 and n = 185) by methylation-specific PCR. Kaplan-Meier curves and log-rank tests were used to estimate cause-specific survival. HRs and corresponding 95% confidence intervals (CI) were assessed using Cox proportional hazard models. To assess the predictive capacity and fit of models combining several methylation markers, HarrellC statistic and the Akaike Information Criterion were used.Results: We identified four methylation markers, that is, GREM1, NEURL, LAD1, and NEFH, that individually predicted prognosis of patients with ccRCC. The four markers combined were associated with poorer survival in two independent patient series (HR, 3.64; 95% CI, 1.02-13.00 and HR, 7.54; 95% CI, 2.68-21.19). These findings were confirmed in a third series of ccRCC cases from The Cancer Genome Atlas (HR, 3.60; 95% CI, 2.02-6.40).Conclusions: A four-gene promoter methylation marker panel consisting of GREM1, NEURL, LAD1, and NEFH predicts outcome of patients with ccRCC and might be used to improve current prognostic models. Clin Cancer Res; 23(8); 2006-18. ©2016 AACR.

Making Sense of Epigenetics.

The gene-environment interactions that underlie development and progression of psychiatric illness are poorly understood. Despite a century of progress, genetic approaches have failed to identify new treatment modalities, perhaps because of the heterogeneity of the disorders and lack of understanding of mechanisms. Recent exploration into epigenetic mechanisms in health and disease has uncovered changes in DNA methylation and chromatin structure that may contribute to psychiatric disorders. Epigenetic changes suggest a variety of new therapeutic options due to their reversible chemistry. However, distinguishing causal links between epigenetic changes and disease from changes consequent to life experience has remained problematic. Here we define epigenetics and explore aspects of epigenetics relevant to causes and mechanisms of psychiatric disease, and speculate on future directions.

Epigenetic silencing of neurofilament genes promotes an aggressive phenotype in breast cancer.

Neurofilament heavy polypeptide (NEFH) has recently been identified as a candidate DNA hypermethylated gene within the functional breast cancer hypermethylome. NEFH exists in a complex with neurofilament medium polypeptide (NEFM) and neurofilament light polypeptide (NEFL) to form neurofilaments, which are structural components of the cytoskeleton in mature neurons. Recent studies reported the deregulation of these proteins in several malignancies, suggesting that neurofilaments may have a role in other cell types as well. Using a comprehensive approach, we studied the epigenetic inactivation of neurofilament genes in breast cancer and the functional significance of this event. We report that DNA methylation-associated silencing of NEFH, NEFL, and NEFM in breast cancer is frequent, cancer-specific, and correlates with clinical features of disease progression. DNA methylation-mediated inactivation of these genes occurs also in multiple other cancer histologies including pancreas, gastric, and colon. Restoration of NEFH function, the major subunit of the neurofilament complex, reduces proliferation and growth of breast cancer cells and arrests them in Go/G1 phase of the cell cycle along with a reduction in migration and invasion. These findings suggest that DNA methylation-mediated silencing of the neurofilament genes NEFH, NEFM, and NEFL are frequent events that may contribute to the progression of breast cancer and possibly other malignancies.

DNA methylation in the medial prefrontal cortex regulates alcohol-induced behavior and plasticity.

Recent studies have suggested an association between alcoholism and DNA methylation, a mechanism that can mediate long-lasting changes in gene transcription. Here, we examined the contribution of DNA methylation to the long-term behavioral and molecular changes induced by a history of alcohol dependence. In search of mechanisms underlying persistent rather than acute dependence-induced neuroadaptations, we studied the role of DNA methylation regulating medial prefrontal cortex (mPFC) gene expression and alcohol-related behaviors in rats 3 weeks into abstinence following alcohol dependence. Postdependent rats showed escalated alcohol intake, which was associated with increased DNA methylation as well as decreased expression of genes encoding synaptic proteins involved in neurotransmitter release in the mPFC. Infusion of the DNA methyltransferase inhibitor RG108 prevented both escalation of alcohol consumption and dependence-induced downregulation of 4 of the 7 transcripts modified in postdependent rats. Specifically, RG108 treatment directly reversed both downregulation of synaptotagmin 2 (Syt2) gene expression and hypermethylation on CpG#5 of its first exon. Lentiviral inhibition of Syt2 expression in the mPFC increased aversion-resistant alcohol drinking, supporting a mechanistic role of Syt2 in compulsive-like behavior. Our findings identified a functional role of DNA methylation in alcohol dependence-like behavioral phenotypes and a candidate gene network that may mediate its effects. Together, these data provide novel evidence for DNA methyltransferases as potential therapeutic targets in alcoholism.

Cannabinoid receptor 1 promotes hepatocellular carcinoma initiation and progression through multiple mechanisms.

UNLABELLED: Hepatocellular carcinoma (HCC) has high mortality and no adequate treatment. Endocannabinoids interact with hepatic cannabinoid 1 receptors (CB1Rs) to promote hepatocyte proliferation in liver regeneration by inducing cell cycle proteins involved in mitotic progression, including Forkhead Box M1. Because this protein is highly expressed in HCC and contributes to its genesis and progression, we analyzed the involvement of the endocannabinoid/CB1R system in murine and human HCC. Postnatal diethylnitrosamine treatment induced HCC within 8 months in wild-type mice but fewer and smaller tumors in CB1R(-/-) mice or in wild-type mice treated with the peripheral CB1R antagonist JD5037, as monitored in vivo by serial magnetic resonance imaging. Genome-wide transcriptome analysis revealed CB1R-dependent, tumor-induced up-regulation of the hepatic expression of CB1R, its endogenous ligand anandamide, and a number of tumor-promoting genes, including the GRB2 interactome as well as Forkhead Box M1 and its downstream target, the tryptophan-catalyzing enzyme indoleamine 2,3-dioxygenase. Increased indoleamine 2,3-dioxygenase activity and consequent induction of immunosuppressive T-regulatory cells in tumor tissue promote immune tolerance. CONCLUSION: The endocannabinoid/CB1R system is up-regulated in chemically induced HCC, resulting in the induction of various tumor-promoting genes, including indoleamine 2,3-dioxygenase; and attenuation of these changes by blockade or genetic ablation of CB1R suppresses the growth of HCC and highlights the therapeutic potential of peripheral CB1R blockade.

Spectrin repeat containing nuclear envelope 1 and forkhead box protein E1 are promising markers for the detection of colorectal cancer in blood.

Identifying biomarkers in body fluids may improve the noninvasive detection of colorectal cancer. Previously, we identified N-Myc downstream-regulated gene 4 (NDRG4) and GATA binding protein 5 (GATA5) methylation as promising biomarkers for colorectal cancer in stool DNA. Here, we examined the utility of NDRG4, GATA5, and two additional markers [Forkhead box protein E1 (FOXE1) and spectrin repeat containing nuclear envelope 1 (SYNE1)] promoter methylation as biomarkers in plasma DNA. Quantitative methylation-specific PCR was performed on plasma DNA from 220 patients with colorectal cancer and 684 noncancer controls, divided in a training set and a test set. Receiver operating characteristic analysis was performed to measure the area under the curve of GATA5, NDRG4, SYNE1, and FOXE1 methylation. Functional assays were performed in SYNE1 and FOXE1 stably transfected cell lines. The sensitivity of NDRG4, GATA5, FOXE1, and SYNE1 methylation in all stages of colorectal cancer (154 cases, 444 controls) was 27% [95% confidence interval (CI), 20%-34%), 18% (95% CI, 12%-24%), 46% (95% CI, 38%-54%), and 47% (95% CI, 39%-55%), with a specificity of 95% (95% CI, 93%-97%), 99% (95% CI, 98%-100%), 93% (95% CI, 91%-95%), and 96% (95% CI, 94%-98%), respectively. Combining SYNE1 and FOXE1, increased the sensitivity to 56% (95% CI, 48%-64%), while the specificity decreased to 90% (95% CI, 87%-93%) in the training set and to 58% sensitivity (95% CI, 46%-70%) and 91% specificity (95% CI, 80%-100%) in a test set (66 cases, 240 controls). SYNE1 overexpression showed no major differences in cell proliferation, migration, and invasion compared with controls. Overexpression of FOXE1 significantly decreased the number of colonies in SW480 and HCT116 cell lines. Overall, our data suggest that SYNE1 and FOXE1 are promising markers for colorectal cancer detection.

Functional identification of cancer-specific methylation of CDO1, HOXA9, and TAC1 for the diagnosis of lung cancer.

PURPOSE: Non-small cell lung cancer (NSCLC) is the leading cause of cancer mortality in the world. Novel diagnostic biomarkers may augment both existing NSCLC screening methods as well as molecular diagnostic tests of surgical specimens to more accurately stratify and stage candidates for adjuvant chemotherapy. Hypermethylation of CpG islands is a common and important alteration in the transition from normal tissue to cancer. EXPERIMENTAL DESIGN: Following previously validated methods for the discovery of cancer-specific hypermethylation changes, we treated eight NSCLC cell lines with the hypomethylating agent deoxyazacitidine or trichostatin A. We validated the findings using a large publicly available database and two independent cohorts of primary samples. RESULTS: We identified >300 candidate genes. Using The Cancer Genome Atlas (TCGA) and extensive filtering to refine our candidate genes for the greatest ability to distinguish tumor from normal, we define a three-gene panel, CDO1, HOXA9, and TAC1, which we subsequently validate in two independent cohorts of primary NSCLC samples. This three-gene panel is 100% specific, showing no methylation in 75 TCGA normal and seven primary normal samples and is 83% to 99% sensitive for NSCLC depending on the cohort. CONCLUSION: This degree of sensitivity and specificity may be of high value to diagnose the earliest stages of NSCLC. Addition of this three-gene panel to other previously validated methylation biomarkers holds great promise in both early diagnosis and molecular staging of NSCLC.

Frequent inactivation of cysteine dioxygenase type 1 contributes to survival of breast cancer cells and resistance to anthracyclines.

PURPOSE: Genome-wide DNA methylation analyses have identified hundreds of candidate DNA-hypermethylated genes in cancer. Comprehensive functional analyses provide an understanding of the biologic significance of this vast amount of DNA methylation data that may allow the determination of key epigenetic events associated with tumorigenesis. EXPERIMENTAL DESIGN: To study mechanisms of cysteine dioxygenase type 1 (CDO1) inactivation and its functional significance in breast cancer in a comprehensive manner, we screened for DNA methylation and gene mutations in primary breast cancers and analyzed growth, survival, and reactive oxygen species (ROS) production in breast cancer cells with restored CDO1 function in the context of anthracycline treatment. RESULTS: DNA methylation-associated silencing of CDO1 in breast cancer is frequent (60%), cancer specific, and correlates with disease progression and outcome. CDO1 function can alternatively be silenced by repressive chromatin, and we describe protein-damaging missense mutations in 7% of tumors without DNA methylation. Restoration of CDO1 function in breast cancer cells increases levels of ROS and leads to reduced viability and growth, as well as sensitization to anthracycline treatment. Priming with 5-azacytidine of breast cancer cells with epigenetically silenced CDO1 resulted in restored expression and increased sensitivity to anthracyclines. CONCLUSION: We report that silencing of CDO1 is a critical epigenetic event that contributes to the survival of oxidative-stressed breast cancer cells through increased detoxification of ROS and thus leads to the resistance to ROS-generating chemotherapeutics including anthracyclines. Our study shows the importance of CDO1 inactivation in breast cancer and its clinical potential as a biomarker and therapeutic target to overcome resistance to anthracyclines.

Biomarkers for detection and prognosis of breast cancer identified by a functional hypermethylome screen.

Breast cancer (BC) is a disease with diverse tumor heterogeneity, which challenges conventional approaches to develop biomarkers for early detection and prognosis. To identify effective biomarkers, we performed a genome-wide screen for functional methylation changes in BC, i.e., genes silenced by promoter hypermethylation, using a functionally proven gene expression approach. A subset of candidate hypermethylated genes were validated in primary BCs and tested as markers for detection and prognosis prediction of BC. We identified 33 cancer specific methylated genes and, among these, two categories of genes: (1) highly frequent methylated genes that detect early stages of BC. Within that category, we have identified the combination of NDRG2 and HOXD1 as the most sensitive (94%) and specific (90%) gene combination for detection of BC; (2) genes that show stage dependent methylation frequency pattern, which are candidates to help delineate BC prognostic signatures. For this category, we found that methylation of CDO1, CKM, CRIP1, KL and TAC1 correlated with clinical prognostic variables and was a significant prognosticator for poor overall survival in BC patients. CKM [Hazard ratio (HR) = 2.68] and TAC1 (HR = 7.73) were the strongest single markers and the combination of both (TAC1 and CKM) was associated with poor overall survival independent of age and stage in our training (HR = 1.92) and validation cohort (HR = 2.87). Our study demonstrates an efficient method to utilize functional methylation changes in BC for the development of effective biomarkers for detection and prognosis prediction of BC.