Benjamin G. Barwick

Age-associated DNA methylation in pediatric populations

DNA methylation (DNAm) plays diverse roles in human biology, but this dynamic epigenetic mark remains far from fully characterized. Although earlier studies uncovered loci that undergo age-associated DNAm changes in adults, little is known about such changes during childhood. Despite profound DNAm plasticity during embryogenesis, monozygotic twins show indistinguishable childhood methylation, suggesting that DNAm is highly coordinated throughout early development. Here we examine the methylation of 27,578 CpG dinucleotides in peripheral blood DNA from a cross-sectional study of 398 boys, aged 3-17 yr, and find significant age-associated changes in DNAm at 2078 loci. These findings correspond well with pyrosequencing data and replicate in a second pediatric population (N = 78). Moreover, we report a deficit of age-related loci on the X chromosome, a preference for specific nucleotides immediately surrounding the interrogated CpG dinucleotide, and a primary association with developmental and immune ontological functions. Meta-analysis (N = 1158) with two adult populations reveals that despite a significant overlap of age-associated loci, most methylation changes do not follow a lifelong linear pattern due to a threefold to fourfold higher rate of change in children compared with adults; consequently, the vast majority of changes are more accurately modeled as a function of logarithmic age. We therefore conclude that age-related DNAm changes in peripheral blood occur more rapidly during childhood and are imperfectly accounted for by statistical corrections that are linear in age, further suggesting that future DNAm studies should be matched closely for age.

A heterozygous IDH1R132H/WT mutation induces genome-wide alterations in DNA methylation

Monoallelic point mutations of the NADP(+)-dependent isocitrate dehydrogenases IDH1 and IDH2 occur frequently in gliomas, acute myeloid leukemias, and chondromas, and display robust association with specific DNA hypermethylation signatures. Here we show that heterozygous expression of the IDH1(R132H) allele is sufficient to induce the genome-wide alterations in DNA methylation characteristic of these tumors. Using a gene-targeting approach, we knocked-in a single copy of the most frequently observed IDH1 mutation, R132H, into a human cancer cell line and profiled changes in DNA methylation at over 27,000 CpG dinucleotides relative to wild-type parental cells. We find that IDH1(R132H/WT) mutation induces widespread alterations in DNA methylation, including hypermethylation of 2010 and hypomethylation of 842 CpG loci. We demonstrate that many of these alterations are consistent with those observed in IDH1-mutant and G-CIMP+ primary gliomas and can segregate IDH wild-type and mutated tumors as well as those exhibiting the G-CIMP phenotype in unsupervised analysis of two primary glioma cohorts. Further, we show that the direction of IDH1(R132H/WT)-mediated DNA methylation change is largely dependent upon preexisting DNA methylation levels, resulting in depletion of moderately methylated loci. Additionally, whereas the levels of multiple histone H3 and H4 methylation modifications were globally increased, consistent with broad inhibition of histone demethylation, hypermethylation at H3K9 in particular accompanied locus-specific DNA hypermethylation at several genes down-regulated in IDH1(R132H/WT) knock-in cells. These data provide insight on epigenetic alterations induced by IDH1 mutations and support a causal role for IDH1(R132H/WT) mutants in driving epigenetic instability in human cancer cells.

Wnt signaling in triple negative breast cancer is associated with metastasis

BackgroundTriple Negative subset of (TN) Breast Cancers (BC), a close associate of the basal-like subtype (with limited discordance) is an aggressive form of the disease which convey unpredictable, and poor prognosis due to limited treatment options and lack of proven effective targeted therapies.MethodsWe conducted an expression study of 240 formalin-fixed, paraffin-embedded (FFPE) primary biopsies from two cohorts, including 130 TN tumors, to identify molecular mechanisms of TN disease.ResultsThe annotation of differentially expressed genes in TN tumors contained an overrepresentation of canonical Wnt signaling components in our cohort and others. These observations were supported by upregulation of experimentally induced oncogenic Wnt/β-catenin genes in TN tumors, recapitulated using targets induced by Wnt3A. A functional blockade of Wnt/β-catenin pathway by either a pharmacological Wnt-antagonist, WntC59, sulidac sulfide, or β-catenin (functional read out of Wnt/β-catenin pathway) SiRNA mediated genetic manipulation demonstrated that a functional perturbation of the pathway is causal to the metastasis- associated phenotypes including fibronectin-directed migration, F-actin organization, and invasion in TNBC cells. A classifier, trained on microarray data from β-catenin transfected mammary cells, identified a disproportionate number of TNBC breast tumors as compared to other breast cancer subtypes in a meta-analysis of 11 studies and 1,878 breast cancer patients, including the two cohorts published here. Patients identified by the Wnt/β-catenin classifier had a greater risk of lung and brain, but not bone metastases.ConclusionThese data implicate transcriptional Wnt signaling as a hallmark of TNBC disease associated with specific metastatic pathways.

Global DNA methylation remodeling accompanies CD8 T cell effector function

The differentiation of CD8 T cells in response to acute infection results in the acquisition of hallmark phenotypic effector functions; however, the epigenetic mechanisms that program this differentiation process on a genome-wide scale are largely unknown. In this article, we report the DNA methylomes of Ag-specific naive and day-8 effector CD8 T cells following acute lymphocytic choriomeningitis virus infection. During effector CD8 T cell differentiation, DNA methylation was remodeled such that changes in DNA methylation at gene promoter regions correlated negatively with gene expression. Importantly, differentially methylated regions were enriched at cis-elements, including enhancers active in naive T cells. Differentially methylated regions were associated with cell type–specific transcription factor binding sites, and these transcription factors clustered into modules that define networks targeted by epigenetic regulation and control of effector CD8 T cell function. Changes in the DNA methylation profile following CD8 T cell activation revealed numerous cellular processes, cis-elements, and transcription factor networks targeted by DNA methylation. Together, the results demonstrated that DNA methylation remodeling accompanies the acquisition of the CD8 T cell effector phenotype and repression of the naive cell state. Therefore, these data provide the framework for an epigenetic mechanism that is required for effector CD8 T cell differentiation and adaptive immune responses.

Prostate cancer genes associated with TMPRSS2–ERG gene fusion and prognostic of biochemical recurrence in multiple cohorts

Background:Recent studies have indicated that prostate cancer patients with the TMPRSS2–ERG gene fusion have a higher risk of recurrence. To identify markers associated with TMPRSS2–ERG fusion and prognostic of biochemical recurrence, we analysed a cohort of 139 men with prostate cancer for 502 molecular markers.Methods:RNA from radical prostatectomy tumour specimens was analysed using cDNA-mediated, annealing, selection, extension and ligation (DASL) to determine mRNAs associated with TMPRSS2–ERG T1/E4 fusion and prognostic of biochemical recurrence. Differentially expressed mRNAs in T1/E4-positive tumours were determined using significance analysis of microarrays (false discovery rate (FDR) <5%). Univariate and multivariate Cox regression determined genes, gene signatures and clinical factors prognostic of recurrence (P-value <0.05, log–rank test). Analysis of two prostate microarray studies (GSE1065 and GSE8402) validated the findings.Results:In the 139 patients from this study and from a 455-patient Swedish cohort, 15 genes in common were differentially regulated in T1/E4 fusion-positive tumours (FDR <0.05). The most significant mRNAs in both cohorts coded ERG. Nine genes were found prognostic of recurrence in this study and in a 596-patient Minnesota cohort. A molecular recurrence score was significant in prognosticating recurrence (P-value 0.000167) and remained significant in multivariate analysis of a mixed clinical model considering Gleason score and TMPRSS2–ERG fusion status.Conclusions:TMPRSS2–ERG T1/E4 fusion-positive tumours had differentially regulated mRNAs observed in multiple studies, the most significant one coded for ERG. Several mRNAs were consistently associated with biochemical recurrence and have potential clinical utility but will require further validation for successful translation.

Heterozygosity for Pten Promotes Tumorigenesis in a Mouse Model of Medulloblastoma

Background Recent publications have described an important role for cross talk between PI-3 kinase and sonic hedgehog signaling pathways in the pathogenesis of medulloblastoma. Methodology/Principal Findings We crossed mice with constitutive activation of Smoothened, SmoA1, with Pten deficient mice. Both constitutive and conditional Pten deficiency doubled the incidence of mice with symptoms of medulloblastoma and resulted in decreased survival. Analysis revealed a clear separation of gene signatures, with up-regulation of genes in the PI-3 kinase signaling pathway, including downstream activation of angiogenesis in SmoA1+/−; Pten +/− medulloblastomas. Western blotting and immunohistochemistry confirmed reduced or absent Pten, Akt activation, and increased angiogenesis in Pten deficient tumors. Down-regulated genes included genes in the sonic hedgehog pathway and tumor suppressor genes. SmoA1+/−; Pten +/+ medulloblastomas appeared classic in histology with increased proliferation and diffuse staining for apoptosis. In contrast, Pten deficient tumors exhibited extensive nodularity with neuronal differentiation separated by focal areas of intense staining for proliferation and virtually absent apoptosis. Examination of human medulloblastomas revealed low to absent PTEN expression in over half of the tumors. Kaplan-Meier analysis confirmed worse overall survival in patients whose tumor exhibited low to absent PTEN expression. Conclusions/Significance This suggests that PTEN expression is a marker of favorable prognosis and mouse models with activation of PI-3 kinase pathways may be important tools for preclinical evaluation of promising agents for the treatment of medulloblastoma.

Differential activation of Wnt-β-catenin pathway in triple negative breast cancer increases MMP7 in a PTEN dependent manner.

Mutations of genes in tumor cells of Triple Negative subset of Breast Cancer (TNBC) deregulate pathways of signal transduction. The loss of tumor suppressor gene PTEN is the most common first event associated with basal-like subtype (Martins, De, Almendro, Gonen, and Park, 2012). Here we report for the first time that the functional upregulation of secreted-MMP7, a transcriptional target of Wnt-β-catenin signature pathway in TNBC is associated to the loss of PTEN. We identified differential expression of mRNAs in several key-components genes, and transcriptional target genes of the Wnt-β-catenin pathway (WP), including beta-catenin, FZD7, DVL1, MMP7, c-MYC, BIRC5, CD44, PPARD, c-MET, and NOTCH1 in FFPE tumors samples from TNBC patients of two independent cohorts. A similar differential upregulation of mRNA/protein for beta-catenin, the functional readout of WP, and for MMP7, a transcriptional target gene of beta-catenin was observed in TNBC cell line models. Genetic or pharmacological attenuation of beta-catenin by SiRNA or WP modulators (XAV939 and sulindac sulfide) and pharmacological mimicking of PTEN following LY294002 treatment downregulated MMP7 levels as well as enzymatic function of the secreted MMP7 in MMP7 positive PTEN-null TNBC cells. Patient data revealed that MMP7 mRNA was high in only a subpopulation of TNBC, and this subpopulation was characterized by a concurrent low expression of PTEN mRNA. In cell lines, a high expression of casein-zymograph-positive MMP7 was distinguished by an absence of functional PTEN. A similar inverse relationship between MMP7 and PTEN mRNA levels was observed in the PAM50 data set (a correlation coefficient of -0.54). The PAM50 subtype and outcome data revealed that the high MMP7 group had low pCR (25%) and High Rd (74%) in clinical stage T3 pathologic response in contrast to the high pCR (40%) and low residual disease (RD) (60%) of the low MMP7 group.

Plasma cell differentiation is coupled to division-dependent DNA hypomethylation and gene regulation

The epigenetic processes that regulate antibody-secreting plasma cells are not well understood. Here, analysis of plasma cell differentiation revealed DNA hypomethylation of 10% of CpG loci that were overrepresented at enhancers. Inhibition of DNA methylation enhanced plasma cell commitment in a cell-division-dependent manner. Analysis of B cells differentiating in vivo stratified by cell division revealed a fivefold increase in mRNA transcription coupled to DNA hypomethylation. Demethylation occurred first at binding motifs for the transcription factors NF-κB and AP-1 and later at those for the transcription factors IRF and Oct-2 and was coincident with activation and differentiation gene-expression programs in a cell-division-dependent manner. These data provide mechanistic insight into cell-division-coupled transcriptional and epigenetic reprogramming and suggest that DNA hypomethylation reflects the cis-regulatory history of plasma cell differentiation.

DNA Extraction from Formalin-Fixed, Paraffin-Embedded Tissue

Formalin-fixed, paraffin-embedded (FFPE) tissue is one of the most widely practiced methods for clinical sample preservation and archiving. It is estimated that, worldwide, over a billion tissue samples, most of them FFPE, are being stored in numerous hospitals, tissue banks, and research laboratories. These archived samples could potentially provide a wealth of information in retrospective molecular studies of diseased tissues. While standard for histopathology and microscopic investigation (e.g., hematoxylin and eosin [H&E] staining, immunohistochemistry [IHC], and tissue microarray [TMA]), FFPE samples pose a major challenge for molecular pathologists, because nucleic acids are heavily modified and trapped by extensive protein-nucleic acid and protein-protein cross linking. Historically, FFPE samples were not considered to be a viable source for molecular analyses. Recently, however, it has been discovered that with appropriate protease digestion, it is possible to release microgram amounts of DNA and RNA from FFPE samples. The purified nucleic acids, although highly fragmented, are suitable for a variety of downstream genomic and gene expression analyses, such as polymerase chain reaction (PCR), quantitative reverse transcription PCR (qRT-PCR), microarray, array comparative genomic hybridization (CGH), microRNA, and methylation profiling. Several commercial kits are currently available for FFPE extraction. The protocol reported here is adapted from the Ambion RecoverAll Total Nucleic Acid Isolation Kit, but includes several modifications. Although our protocol focuses on DNA isolation, the RecoverAll Kit can also be utilized to recover RNA, including microRNA.

MAX is an epigenetic sensor of 5-carboxylcytosine and is altered in multiple myeloma.

Abstract The oncogenic transcription factor MYC and its binding partner MAX regulate gene expression by binding to DNA at enhancer-box (E-box) elements 5΄-CACGTG-3΄. In mammalian genomes, the central E-box CpG has the potential to be methylated at the 5-position of cytosine (5mC), or to undergo further oxidation to the 5-hydroxymethyl (5hmC), 5-formyl (5fC), or 5-carboxyl (5caC) forms. We find that MAX exhibits the greatest affinity for a 5caC or unmodified C-containing E-box, and much reduced affinities for the corresponding 5mC, 5hmC or 5fC forms. Crystallization of MAX with a 5caC modified E-box oligonucleotide revealed that MAX Arg36 recognizes 5caC using a 5caC–Arg–Guanine triad, with the next nearest residue to the carboxylate group being Arg60. In an analysis of >800 primary multiple myelomas, MAX alterations occurred at a frequency of ∼3%, more than half of which were single nucleotide substitutions affecting a basic clamp-like interface important for DNA interaction. Among these, arginines 35, 36 and 60 were the most frequently altered. In vitro binding studies showed that whereas mutation of Arg36 (R36W) or Arg35 (R35H/L) completely abolished DNA binding, mutation of Arg60 (R60Q) significantly reduced DNA binding, but retained a preference for the 5caC modified E-box. Interestingly, MAX alterations define a subset of myeloma patients with lower MYC expression and a better overall prognosis. Together these data indicate that MAX can act as a direct epigenetic sensor of E-box cytosine modification states and that local CpG modification and MAX variants converge to modulate the MAX-MYC transcriptional network.