Dynamic interplay between locus-specific DNA methylation and hydroxymethylation regulates distinct biological pathways in prostate carcinogenesis

Understanding the role of Pax5 in development of taxane-resistant neuroendocrine like prostate cancers

Resistance to the current Androgen Receptor Signaling Inhibitor (ARSI) therapies has led to higher incidences of therapy-induced neuroendocrine-like prostate cancer (t-NEPC). This highly aggressive subtype with predominant small cell-like characteristics is resistant to taxane chemotherapies and has a dismal overall survival. t-NEPCs are mostly treated with platinum-based drugs with a combination of etoposide or taxane and have less selectivity and high systemic toxicity, which often limit their clinical potential. During t-NEPC transformation, adenocarcinomas lose their luminal features and adopt neuro-basal characteristics. Whether the adaptive neuronal characteristics of t-NEPC are responsible for such taxane resistance remains unknown. Pathway analysis from patient gene-expression databases indicates that t-NEPC upregulates various neuronal pathways associated with enhanced cellular networks. To identify transcription factor(s) (TF) that could be important for promoting the gene expression for neuronal characters in t-NEPC, we performed ATAC-Seq, acetylated-histone ChIP-seq, and RNA-seq in our NE-like cell line models and analyzed the promoters of transcriptionally active and significantly enriched neuroendocrine-like (NE-like) cancer-specific genes. Our results indicate that Pax5 could be an important transcription factor for neuronal gene expression and specific to t-NEPC. Pathway analysis revealed that Pax5 expression is involved in axonal guidance, neurotransmitter regulation, and neuronal adhesion, which are critical for strong cellular communications. Further results suggest that depletion of Pax5 disrupts cellular interaction in NE-like cells and reduces surface growth factor receptor activation, thereby, sensitizing them to taxane therapies. Moreover, t-NEPC specific hydroxymethylation of Pax5 promoter CpG islands favors Pbx1 binding to induce Pax5 expression. Based on our study, we concluded that continuous exposure to ARSI therapies leads to epigenetic modifications and Pax5 activation in t-NEPC, which promotes the expression of genes necessary to adopt taxane-resistant NE-like cancer. Thus, targeting the Pax5 axis can be beneficial for reverting their taxane sensitivity.

Alterations in DNA 5-hydroxymethylation Patterns in the Hippocampus of an Experimental Model of Refractory Epilepsy

Temporal lobe epilepsy (TLE) is a type of focal epilepsy characterized by spontaneous recurrent seizures originating from the hippocampus. The epigenetic reprogramming hypothesis of epileptogenesis suggests that the development of TLE is associated with alterations in gene transcription changes resulting in a hyperexcitable network in TLE. DNA 5-methylcytosine (5-mC) is an epigenetic mechanism that has been associated with chronic epilepsy. However, the contribution of 5-hydroxymethylcytosine (5-hmC), a product of 5-mC demethylation by the Ten-Eleven Translocation (TET) family proteins in chronic TLE is poorly understood. 5-hmC is abundant in the brain and acts as a stable epigenetic mark altering gene expression through several mechanisms. Here, we found that the levels of bulk DNA 5-hmC but not 5-mC were significantly reduced in the hippocampus of human TLE patients and in the kainic acid (KA) TLE rat model. Using 5-hmC hMeDIP-sequencing, we characterized 5-hmC distribution across the genome and found bidirectional regulation of 5-hmC at intergenic regions within gene bodies. We found that hypohydroxymethylated 5-hmC intergenic regions were associated with several epilepsy-related genes, including Gal , SV2, and Kcnj11 and hyperdroxymethylation 5-hmC intergenic regions were associated with Gad65 , TLR4 , and Bdnf gene expression. Mechanistically, Tet1 knockdown in the hippocampus was sufficient to decrease 5-hmC levels and increase seizure susceptibility following KA administration. In contrast, Tet1 overexpression in the hippocampus resulted in increased 5-hmC levels associated with improved seizure resiliency in response to KA. These findings suggest an important role for 5-hmC as an epigenetic regulator of epilepsy that can be manipulated to influence seizure outcomes.

The 5-Hydroxymethylcytosine Landscape of Prostate Cancer

Analysis of DNA methylation is a valuable tool to understand disease progression and is increasingly being used to create diagnostic and prognostic clinical biomarkers. While conversion of cytosine to 5-methylcytosine (5mC) commonly results in transcriptional repression, further conversion to 5-hydroxymethylcytosine (5hmC) is associated with transcriptional activation. Here we perform the first study integrating whole-genome 5hmC with DNA, 5mC, and transcriptome sequencing in clinical samples of benign, localized, and advanced prostate cancer. 5hmC is shown to mark activation of cancer drivers and downstream targets. Furthermore, 5hmC sequencing revealed profoundly altered cell states throughout the disease course, characterized by increased proliferation, oncogenic signaling, dedifferentiation, and lineage plasticity to neuroendocrine and gastrointestinal lineages. Finally, 5hmC sequencing of cell-free DNA from patients with metastatic disease proved useful as a prognostic biomarker able to identify an aggressive subtype of prostate cancer using the genes TOP2A and EZH2, previously only detectable by transcriptomic analysis of solid tumor biopsies. Overall, these findings reveal that 5hmC marks epigenomic activation in prostate cancer and identify hallmarks of prostate cancer progression with potential as biomarkers of aggressive disease.

SIGNIFICANCE

In prostate cancer, 5-hydroxymethylcytosine delineates oncogene activation and stage-specific cell states and can be analyzed in liquid biopsies to detect cancer phenotypes. See related article by Wu and Attard, p. 3880.

BigMPI4py: Python Module for Parallelization of Big Data Objects Discloses Germ Layer Specific DNA Demethylation Motifs

Parallelization in Python integrates Message Passing Interface via the mpi4py module. Since mpi4py does not support parallelization of objects greater than 2³¹ bytes, we developed BigMPI4py, a Python module that wraps mpi4py, supporting object sizes beyond this boundary. BigMPI4py automatically determines the optimal object distribution strategy, and uses vectorized methods, achieving higher parallelization efficiency. BigMPI4py facilitates the implementation of Python for Big Data applications in multicore workstations and High Performance Computer systems. We use BigMPI4py to speed-up the search for germ line specific de novo DNA methylated/unmethylated motifs from the 59 whole genome bisulfite sequencing DNA methylation samples from 27 human tissues of the ENCODE project. We developed a parallel implementation of the Kruskall-Wallis test to find CpGs with differential methylation across germ layers. The parallel evaluation of the significance of 55 million CpG achieved a 22x speedup with 25 cores allowing us an efficient identification of a set of hypermethylated genes in ectoderm and mesoderm-related tissues, and another set in endoderm-related tissues and finally, the discovery of germ layer specific DNA demethylation motifs. Our results point out that DNA methylation signal provide a higher degree of information for the demethylated state than for the methylated state. BigMPI4py is available at https://https://www.arauzolab.org/tools/bigmpi4py and https://gitlab.com/alexmascension/bigmpi4py and the Jupyter Notebook with WGBS analysis at https://gitlab.com/alexmascension/wgbs-analysis.

Epigenetic reprogramming during prostate cancer progression: A perspective from development

Conrad Waddington's theory of epigenetic landscape epitomize the process of cell fate and cellular decision-making during development. Wherein the epigenetic code maintains patterns of gene expression in pluripotent and differentiated cellular states during embryonic development and differentiation. Over the years disruption or reprogramming of the epigenetic landscape has been extensively studied in the course of cancer progression. Cellular dedifferentiation being a key hallmark of cancer allow us to take cues from the biological processes involved during development. Here, we discuss the role of epigenetic landscape and its modifiers in cell-fate determination, differentiation and prostate cancer progression. Lately, the emergence of RNA-modifications has also furthered our understanding of epigenetics in cancer. The overview of the epigenetic code regulating androgen signalling, and progression to aggressive neuroendocrine stage of PCa reinforces its gene regulatory functions during the development of prostate gland as well as cancer progression. Additionally, we also highlight the clinical implications of cancer cell epigenome, and discuss the recent advancements in the therapeutic strategies targeting the advanced stage disease.

Holding on to Junk Bonds: Intron Retention in Cancer and Therapy

Intron retention (IR) in cancer was for a long time overlooked by the scientific community, as it was previously considered to be an artifact of a dysfunctional spliceosome. Technological advancements made in the last decade offer unique opportunities to explore the role of IR as a widespread phenomenon that contributes to the transcriptional diversity of many cancers. Numerous studies in cancer have shed light on dysregulation of cellular mechanisms that lead to aberrant and pathologic IR. IR is not merely a mechanism of gene regulation, but rather it can mediate cancer pathogenesis and therapeutic resistance in various human diseases. The burden of IR in cancer is governed by perturbations to mechanisms known to regulate this phenomenon and include epigenetic variation, mutations within the gene body, and splicing factor dysregulation. This review summarizes possible causes for aberrant IR and discusses the role of IR in therapy or as a consequence of disease treatment. As neoepitopes originating from retained introns can be presented on the cancer cell surface, the development of personalized cancer vaccines based on IR-derived neoepitopes should be considered. Ultimately, a deeper comprehension about the origins and consequences of aberrant IR may aid in the development of such personalized cancer vaccines.

Hydroxymethylation and tumors: can 5-hydroxymethylation be used as a marker for tumor diagnosis and treatment?

5-Methylcytosine (5mC) is considered as a common epigenetic modification that plays an important role in the regulation of gene expression. At the same time, 5-hydroxymethylcytosine (5hmC) has been found as an emerging modification of cytosine bases of recent years. Unlike 5mC, global 5hmC levels vary from tissues that have differential distribution both in mammalian tissues and in the genome. DNA hydroxymethylation is the process that 5mC oxidates into 5hmC with the catalysis of TET (ten-eleven translocation) enzymes. It is an essential option of DNA demethylation, which modulates gene expression by adjusting the DNA methylation level. Various factors can regulate the demethylation of DNA, such as environmental toxins and mental stress. In this review, we summarize the progress in the formation of 5hmC, and obtaining 5hmC in a cell-free DNA sample presents multiple advantages and challenges for the subject. Furthermore, the clinical potential for 5hmC modification in dealing with cancer early diagnosis, prognostic evaluation, and prediction of therapeutic effect is also mentioned.

DNA Methylation and Hydroxymethylation in Cervical Cancer: Diagnosis, Prognosis and Treatment

Recent discoveries have led to the development of novel ideas and techniques that have helped elucidate the correlation between epigenetics and tumor biology. Nowadays, the field of tumor genetics has evolved to include a new type of regulation by epigenetics. An increasing number of studies have demonstrated the importance of DNA methylation and hydroxymethylation in specific genes in the progression of cervical cancer. Determining the methylation and hydroxymethylation profiles of these genes will help in the early prevention and diagnosis, monitoring recurrence, prognosis, and treatment of patients with cervical cancer. In this review, we focus on the significance of aberrant DNA methylation and hydroxymethylation in cervical cancer and the use of these epigenetic signatures in clinical settings.

Exploring targets of TET2-mediated methylation reprogramming as potential discriminators of prostate cancer progression

Background

Global DNA methylation alterations are hallmarks of cancer. The tumor-suppressive TET enzymes, which are involved in DNA demethylation, are decreased in prostate cancer (PCa); in particular, TET2 is specifically targeted by androgen-dependent mechanisms of repression in PCa and may play a central role in carcinogenesis. Thus, the identification of key genes targeted by TET2 dysregulation may provide further insight into cancer biology.

Results

Using a CRISPR/Cas9-derived TET2-knockout prostate cell line, and through whole-transcriptome and whole-methylome sequencing, we identified seven candidate genes—ASB2, ETNK2, MEIS2, NRG1, NTN1, NUDT10, and SRPX—exhibiting reduced expression and increased promoter methylation, a pattern characteristic of tumor suppressors. Decreased expression of these genes significantly discriminates between recurrent and non-recurrent prostate tumors from the Cancer Genome Atlas (TCGA) cohort (n = 423), and ASB2, NUDT10, and SRPX were significantly correlated with lower recurrence-free survival in patients by Kaplan-Meier analysis. ASB2, MEIS2, and SRPX also showed significantly lower expression in high-risk Gleason score 8 tumors as compared to low or intermediate risk tumors, suggesting that these genes may be particularly useful as indicators of PCa progression. Furthermore, methylation array probes in the TCGA dataset, which were proximal to the highly conserved, differentially methylated sites identified in our TET2-knockout cells, were able to significantly distinguish between matched prostate tumor and normal prostate tissues (n = 50 pairs). Except ASB2, all genes exhibited significantly increased methylation at these probes, and methylation status of at least one probe for each of these genes showed association with measures of PCa progression such as recurrence, stage, or Gleason score. Since ASB2 did not have any probes within the TET2-knockout differentially methylated region, we validated ASB2 methylation in an independent series of matched tumor-normal samples (n = 19) by methylation-specific qPCR, which revealed concordant and significant increases in promoter methylation within the TET2-knockout site.

Conclusions

Our study identifies seven genes governed by TET2 loss in PCa which exhibit an association between their methylation and expression status and measures of PCa progression. As differential methylation profiles and TET2 expression are associated with advanced PCa, further investigation of these specialized TET2 targets may provide important insights into patterns of carcinogenic gene dysregulation.

Electronic supplementary material

The online version of this article (10.1186/s13148-019-0651-z) contains supplementary material, which is available to authorized users.

5hmC Level Predicts Biochemical Failure Following Radical Prostatectomy in Prostate Cancer Patients with ERG Negative Tumors

This study aimed to validate whether 5-hydroxymethylcytosine (5hmC) level in combination with ERG expression is a predictive biomarker for biochemical failure (BF) in men undergoing radical prostatectomy (RP) for prostate cancer (PCa). The study included 592 PCa patients from two consecutive Danish RP cohorts. 5hmC level and ERG expression were analyzed using immunohistochemistry in RP specimens. 5hmC was scored as low or high and ERG was scored as negative or positive. Risk of BF was analyzed using stratified cumulative incidences and multiple cause-specific Cox regression using competing risk assessment. Median follow-up was 10 years (95% CI: 9.5–10.2). In total, 246 patients (41.6%) had low and 346 patients (58.4%) had high 5hmC level. No significant association was found between 5hmC level or ERG expression and time to BF (p = 0.2 and p = 1.0, respectively). However, for men with ERG negative tumors, high 5hmC level was associated with increased risk of BF following RP (p = 0.01). In multiple cause-specific Cox regression analyses of ERG negative patients, high 5hmC expression was associated with time to BF (HR: 1.8; 95% CI: 1.2–2.7; p = 0.003). In conclusion, high 5hmC level was correlated with time to BF in men with ERG negative PCa, which is in accordance with previous results.

Pathway enrichment analysis and visualization of omics data using g:Profiler, GSEA, Cytoscape and EnrichmentMap

Pathway enrichment analysis helps researchers gain mechanistic insight into gene lists generated from genome-scale (omics) experiments. This method identifies biological pathways that are enriched in a gene list more than would be expected by chance. We explain the procedures of pathway enrichment analysis and present a practical step-by-step guide to help interpret gene lists resulting from RNA-seq and genome-sequencing experiments. The protocol comprises three major steps: definition of a gene list from omics data, determination of statistically enriched pathways, and visualization and interpretation of the results. We describe how to use this protocol with published examples of differentially expressed genes and mutated cancer genes; however, the principles can be applied to diverse types of omics data. The protocol describes innovative visualization techniques, provides comprehensive background and troubleshooting guidelines, and uses freely available and frequently updated software, including g:Profiler, Gene Set Enrichment Analysis (GSEA), Cytoscape and EnrichmentMap. The complete protocol can be performed in ~4.5 h and is designed for use by biologists with no prior bioinformatics training.

Global distribution of DNA hydroxymethylation and DNA methylation in chronic lymphocytic leukemia

BACKGROUND

Chronic lymphocytic leukemia (CLL) has been a good model system to understand the functional role of 5-methylcytosine (5-mC) in cancer progression. More recently, an oxidized form of 5-mC, 5-hydroxymethylcytosine (5-hmC) has gained lot of attention as a regulatory epigenetic modification with prognostic and diagnostic implications for several cancers. However, there is no global study exploring the role of 5-hydroxymethylcytosine (5-hmC) levels in CLL. Herein, using mass spectrometry and hMeDIP-sequencing, we analysed the dynamics of 5-hmC during B cell maturation and CLL pathogenesis.

RESULTS

We show that naïve B-cells had higher levels of 5-hmC and 5-mC compared to non-class switched and class-switched memory B-cells. We found a significant decrease in global 5-mC levels in CLL patients (n = 15) compared to naïve and memory B cells, with no changes detected between the CLL prognostic groups. On the other hand, global 5-hmC levels of CLL patients were similar to memory B cells and reduced compared to naïve B cells. Interestingly, 5-hmC levels were increased at regulatory regions such as gene-body, CpG island shores and shelves and 5-hmC distribution over the gene-body positively correlated with degree of transcriptional activity. Importantly, CLL samples showed aberrant 5-hmC and 5-mC pattern over gene-body compared to well-defined patterns in normal B-cells. Integrated analysis of 5-hmC and RNA-sequencing from CLL datasets identified three novel oncogenic drivers that could have potential roles in CLL development and progression.

CONCLUSIONS

Thus, our study suggests that the global loss of 5-hmC, accompanied by its significant increase at the gene regulatory regions, constitute a novel hallmark of CLL pathogenesis. Our combined analysis of 5-mC and 5-hmC sequencing provided insights into the potential role of 5-hmC in modulating gene expression changes during CLL pathogenesis.

Prostate Cancer Genomic Subtypes

Over the last decade, advancements in massively-parallel DNA sequencing and computational biology have allowed for unprecedented insights into the fundamental mutational processes that underlie virtually every major cancer type. Two major cancer genomics consortia-The Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium (ICGC)-have produced rich databases of mutational, pathological, and clinical data that can be mined through web-based portals, allowing for correlative studies and testing of novel hypotheses on well-powered patient cohorts.In this chapter, we will review the impact of these technological developments on the understanding of molecular subtypes that promote prostate cancer initiation, progression, metastasis, and clinical aggression. In particular, we will focus on molecular subtypes that define clinically-relevant patient cohorts and assess how a better understanding of how these subtypes-in both somatic and germline genomes-may influence the clinical course for individual men diagnosed with prostate cancer.

BS-virus-finder: virus integration calling using bisulfite sequencing data

Background

DNA methylation plays a key role in the regulation of gene expression and carcinogenesis. Bisulfite sequencing studies mainly focus on calling single nucleotide polymorphism, different methylation region, and find allele-specific DNA methylation. Until now, only a few software tools have focused on virus integration using bisulfite sequencing data.

Findings

We have developed a new and easy-to-use software tool, named BS-virus-finder (BSVF, RRID:SCR_015727), to detect viral integration breakpoints in whole human genomes. The tool is hosted at https://github.com/BGI-SZ/BSVF.

Conclusions

BS-virus-finder demonstrates high sensitivity and specificity. It is useful in epigenetic studies and to reveal the relationship between viral integration and DNA methylation. BS-virus-finder is the first software tool to detect virus integration loci by using bisulfite sequencing data.

GBX2 Methylation Is a Novel Prognostic Biomarker and Improves Prediction of Biochemical Recurrence Among Patients with Prostate Cancer Negative for Intraductal Carcinoma and Cribriform Architecture

BACKGROUND

Tumor intraductal carcinoma/cribriform architecture (IDC/C) is associated with an unfavorable prognosis and biochemical recurrence (BCR) in prostate cancer (PCa). Up to 70% of PCa patients are IDC/C-negative, but it is estimated that 20% of these cases still experience BCR. Thus, biomarkers for better detection of aggressive disease in IDC/C-negative patients are required.

OBJECTIVE

To investigate tumor-specific methylation of the transcription factor GBX2 as a novel prognosticator and predictor of BCR in PCa patients stratified by histopathologic features including IDC/C.

DESIGN, SETTING, AND PARTICIPANTS

Using genome-wide methylome profiling, we identified higher GBX2 methylation in grade group (GG) 4 tumors compared to GG1 (discovery cohort). The prognostic nature of GBX2 methylation was validated in silico using The Cancer Genome Atlas data (n=478) and a quantitative methylation assay for radical prostatectomy samples (n=254). Regulation of GBX2 methylation was investigated in prostate cells using methyl-CpG-binding domain sequencing and methylation analysis in functional knockouts of TET2, a key epigenetic player in prostate carcinogenesis.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS

The association of GBX2 methylation with Gleason score (GS), pathologic stage (pT), IDC/C, and BCR was analyzed using Kruskal-Wallis and Mann-Whitney tests. Univariate and multivariate Cox regression analyses were used to predict BCR.

RESULTS

GBX2 methylation was associated with GS (p<0.05), pT (p<0.01), and BCR (p<0.05). GBX2 methylation (p=0.004), GS (p<0.001), pT (p=0.012), and prostate-specific antigen (p=0.005) were independent predictors of BCR. Among IDC/C-negative patients, GBX2 methylation improved prediction of BCR (p=0.002). Loss of TET2 in prostate cells resulted in greater GBX2 methylation.

CONCLUSIONS

We identified GBX2 methylation as a novel prognostic factor in PCa and an independent predictor of BCR. We demonstrated the additive value of GBX2 methylation in predicting BCR among IDC/C-negative patients and elucidated a novel TET2-mediated upstream epigenetic regulatory mechanism of GBX2.

PATIENT SUMMARY

We identified GBX2 methylation as a promising prognostic biomarker that could improve the identification of prostate cancer patients at higher risk of biochemical recurrence.

Clinical Significance of Measuring Global Hydroxymethylation of White Blood Cell DNA in Prostate Cancer: Comparison to PSA in a Pilot Exploratory Study

This is the first study investigating the clinical relevance of 5-hydroxymethylcytosine (5hmC) in genomic DNA from white blood cells (WBC) in the context of prostate cancer (PCa) and other prostate pathologies. Using an enzyme-linked immunosorbent assay, we identified significantly different distributions of patients with low and elevated 5hmC content in WBC DNA across controls and patients with prostate cancer (PCa), atypical small acinar proliferation (ASAP), and benign prostatic hyperplasia (BPH). The measured values were within the normal range for most PCa patients, while the latter category was predominant for ASAP. We observed a wider heterogeneity in 5hmC content in all of the prostate pathologies analyzed when compared to the healthy age-matched controls. When compared to blood levels of prostate-specific antigen (PSA), this 5hmC-based biomarker had a lower performance in PCa detection than the use of a PSA cut-off of 2.5 nanograms per milliliter (ng/mL). Above this threshold, however, it delineated almost three quarters of PCa patients from controls and patients with other prostate pathologies. Overall, genome-wide 5hmC content of WBC DNA appears to be applicable for detecting non-cancerous prostate diseases, rather than PCa. Our results also suggest a potential clinical usefulness of complementing PSA as a PCa marker by the addition of a set of hydroxymethylation markers in the blood, but further studies are necessary to confirm these findings.

Prostate Cancer: Epigenetic Alterations, Risk Factors, and Therapy

Prostate cancer (PCa) is the most prevalent urological cancer that affects aging men in South Africa, and mechanisms underlying prostate tumorigenesis remain elusive. Research advancements in the field of PCa and epigenetics have allowed for the identification of specific alterations that occur beyond genetics but are still critically important in the pathogenesis of tumorigenesis. Anomalous epigenetic changes associated with PCa include histone modifications, DNA methylation, and noncoding miRNA. These mechanisms regulate and silence hundreds of target genes including some which are key components of cellular signalling pathways that, when perturbed, promote tumorigenesis. Elucidation of mechanisms underlying epigenetic alterations and the manner in which these mechanisms interact in regulating gene transcription in PCa are an unmet necessity that may lead to novel chemotherapeutic approaches. This will, therefore, aid in developing combination therapies that will target multiple epigenetic pathways, which can be used in conjunction with the current conventional PCa treatment.