Shichun Zheng

Aging and environmental exposures alter tissue-specific DNA methylation dependent upon CPG island context

Epigenetic control of gene transcription is critical for normal human development and cellular differentiation. While alterations of epigenetic marks such as DNA methylation have been linked to cancers and many other human diseases, interindividual epigenetic variations in normal tissues due to aging, environmental factors, or innate susceptibility are poorly characterized. The plasticity, tissue-specific nature, and variability of gene expression are related to epigenomic states that vary across individuals. Thus, population-based investigations are needed to further our understanding of the fundamental dynamics of normal individual epigenomes. We analyzed 217 non-pathologic human tissues from 10 anatomic sites at 1,413 autosomal CpG loci associated with 773 genes to investigate tissue-specific differences in DNA methylation and to discern how aging and exposures contribute to normal variation in methylation. Methylation profile classes derived from unsupervised modeling were significantly associated with age (P<0.0001) and were significant predictors of tissue origin (P<0.0001). In solid tissues (n = 119) we found striking, highly significant CpG island–dependent correlations between age and methylation; loci in CpG islands gained methylation with age, loci not in CpG islands lost methylation with age (P<0.001), and this pattern was consistent across tissues and in an analysis of blood-derived DNA. Our data clearly demonstrate age- and exposure-related differences in tissue-specific methylation and significant age-associated methylation patterns which are CpG island context-dependent. This work provides novel insight into the role of aging and the environment in susceptibility to diseases such as cancer and critically informs the field of epigenomics by providing evidence of epigenetic dysregulation by age-related methylation alterations. Collectively we reveal key issues to consider both in the construction of reference and disease-related epigenomes and in the interpretation of potentially pathologically important alterations.

Glioma Groups Based on 1p/19q, IDH, and TERT Promoter Mutations in Tumors

BACKGROUND The prediction of clinical behavior, response to therapy, and outcome of infiltrative glioma is challenging. On the basis of previous studies of tumor biology, we defined five glioma molecular groups with the use of three alterations: mutations in the TERT promoter, mutations in IDH, and codeletion of chromosome arms 1p and 19q (1p/19q codeletion). We tested the hypothesis that within groups based on these features, tumors would have similar clinical variables, acquired somatic alterations, and germline variants. METHODS We scored tumors as negative or positive for each of these markers in 1087 gliomas and compared acquired alterations and patient characteristics among the five primary molecular groups. Using 11,590 controls, we assessed associations between these groups and known glioma germline variants. RESULTS Among 615 grade II or III gliomas, 29% had all three alterations (i.e., were triple-positive), 5% had TERT and IDH mutations, 45% had only IDH mutations, 7% were triple-negative, and 10% had only TERT mutations; 5% had other combinations. Among 472 grade IV gliomas, less than 1% were triple-positive, 2% had TERT and IDH mutations, 7% had only IDH mutations, 17% were triple-negative, and 74% had only TERT mutations. The mean age at diagnosis was lowest (37 years) among patients who had gliomas with only IDH mutations and was highest (59 years) among patients who had gliomas with only TERT mutations. The molecular groups were independently associated with overall survival among patients with grade II or III gliomas but not among patients with grade IV gliomas. The molecular groups were associated with specific germline variants. CONCLUSIONS Gliomas were classified into five principal groups on the basis of three tumor markers. The groups had different ages at onset, overall survival, and associations with germline variants, which implies that they are characterized by distinct mechanisms of pathogenesis. (Funded by the National Institutes of Health and others.).

DNA Methylation, Isocitrate Dehydrogenase Mutation, and Survival in Glioma

BACKGROUND Although much is known about molecular and chromosomal characteristics that distinguish glioma histological subtypes, DNA methylation patterns of gliomas and their association with other tumor features such as mutation of isocitrate dehydrogenase (IDH) genes have only recently begun to be investigated. METHODS DNA methylation of glioblastomas, astrocytomas, oligodendrogliomas, oligoastrocytomas, ependymomas, and pilocytic astrocytomas (n = 131) from the Brain Tumor Research Center at the University of California San Francisco, as well as nontumor brain tissues (n = 7), was assessed with the Illumina GoldenGate methylation array. Methylation data were subjected to recursively partitioned mixture modeling (RPMM) to derive methylation classes. Differential DNA methylation between tumor and nontumor was also assessed. The association between methylation class and IDH mutation (IDH1 and IDH2) was tested using univariate and multivariable analysis for tumors (n = 95) with available substrate for sequencing. Survival of glioma patients carrying mutant IDH (n = 57) was compared with patients carrying wild-type IDH (n = 38) using a multivariable Cox proportional hazards model and Kaplan-Meier analysis. All statistical tests were two-sided. RESULTS We observed a statistically significant association between RPMM methylation class and glioma histological subtype (P < 2.2 × 10(-16)). Compared with nontumor brain tissues, across glioma tumor histological subtypes, the differential methylation ratios of CpG loci were statistically significantly different (permutation P < .0001). Methylation class was strongly associated with IDH mutation in gliomas (P = 3.0 × 10(-16)). Compared with glioma patients whose tumors harbored wild-type IDH, patients whose tumors harbored mutant IDH showed statistically significantly improved survival (hazard ratio of death = 0.27, 95% confidence interval = 0.10 to 0.72). CONCLUSION The homogeneity of methylation classes for gliomas with IDH mutation, despite their histological diversity, suggests that IDH mutation is associated with a distinct DNA methylation phenotype and an altered metabolic profile in glioma.

Model-based clustering of DNA methylation array data: a recursive-partitioning algorithm for high-dimensional data arising as a mixture of beta distributions.

BackgroundEpigenetics is the study of heritable changes in gene function that cannot be explained by changes in DNA sequence. One of the most commonly studied epigenetic alterations is cytosine methylation, which is a well recognized mechanism of epigenetic gene silencing and often occurs at tumor suppressor gene loci in human cancer. Arrays are now being used to study DNA methylation at a large number of loci; for example, the Illumina GoldenGate platform assesses DNA methylation at 1505 loci associated with over 800 cancer-related genes. Model-based cluster analysis is often used to identify DNA methylation subgroups in data, but it is unclear how to cluster DNA methylation data from arrays in a scalable and reliable manner.ResultsWe propose a novel model-based recursive-partitioning algorithm to navigate clusters in a beta mixture model. We present simulations that show that the method is more reliable than competing nonparametric clustering approaches, and is at least as reliable as conventional mixture model methods. We also show that our proposed method is more computationally efficient than conventional mixture model approaches. We demonstrate our method on the normal tissue samples and show that the clusters are associated with tissue type as well as age.ConclusionOur proposed recursively-partitioned mixture model is an effective and computationally efficient method for clustering DNA methylation data.

Methylation of the PTEN promoter defines low-grade gliomas and secondary glioblastoma

Glioblastoma multiforme (GBM) can present as either de novo or secondary tumors arising from previously diagnosed low-grade gliomas. Although these tumor types are phenotypically indistinguishable, de novo and secondary GBMs are associated with distinct genetic characteristics. PTEN mutations, which result in activation of the phosphoinositide 3-kinase (PI3K) signal transduction pathway, are frequent in de novo but not in secondary GBMs or their antecedent low-grade tumors. Results we present here show that grade II astrocytomas, oligodendrogliomas, and oligoastrocytomas commonly display methylation of the PTEN promoter, a finding that is absent in nontumor brain specimens and rare in de novo GBMs. Methylation of the PTEN promoter correlates with protein kinase B (PKB/Akt) phosphorylation, reflecting functional activation of the PI3K pathway. Our results also demonstrate frequent methylation of the PTEN promoter in grade III astrocytomas and secondary GBMs, consistent with the hypothesis that these tumors arise from lower grade precursors. PTEN methylation is rare in de novo GBMs and is mutually exclusive with PTEN mutations. We conclude that methylation of the PTEN promoter may represent an alternate mechanism by which PI3K signaling is increased in grade II and III gliomas as well as secondary GBMs, a finding that offers new therapeutic approaches in these patients.

Breast Cancer DNA Methylation Profiles Are Associated with Tumor Size and Alcohol and Folate Intake

Although tumor size and lymph node involvement are the current cornerstones of breast cancer prognosis, they have not been extensively explored in relation to tumor methylation attributes in conjunction with other tumor and patient dietary and hormonal characteristics. Using primary breast tumors from 162 (AJCC stage I–IV) women from the Kaiser Division of Research Pathways Study and the Illumina GoldenGate methylation bead-array platform, we measured 1,413 autosomal CpG loci associated with 773 cancer-related genes and validated select CpG loci with Sequenom EpiTYPER. Tumor grade, size, estrogen and progesterone receptor status, and triple negative status were significantly (Q-values <0.05) associated with altered methylation of 209, 74, 183, 69, and 130 loci, respectively. Unsupervised clustering, using a recursively partitioned mixture model (RPMM), of all autosomal CpG loci revealed eight distinct methylation classes. Methylation class membership was significantly associated with patient race (P<0.02) and tumor size (P<0.001) in univariate tests. Using multinomial logistic regression to adjust for potential confounders, patient age and tumor size, as well as known disease risk factors of alcohol intake and total dietary folate, were all significantly (P<0.0001) associated with methylation class membership. Breast cancer prognostic characteristics and risk-related exposures appear to be associated with gene-specific tumor methylation, as well as overall methylation patterns.

Epigenetic profiles distinguish pleural mesothelioma from normal pleura and predict lung asbestos burden and clinical outcome

Mechanisms of action of nonmutagenic carcinogens such as asbestos remain poorly characterized. As pleural mesothelioma is known to have limited numbers of genetic mutations, we aimed to characterize the relationships among gene-locus-specific methylation alterations, disease status, asbestos burden, and survival in this rapidly fatal asbestos-associated tumor. Methylation of 1505 CpG loci associated with 803 cancer-related genes were studied in 158 pleural mesotheliomas and 18 normal pleura. After false-discovery rate correction, 969 CpG loci were independently associated with disease status (Q < 0.05). Classifying samples based on CpG methylation profile with a mixture model approach, methylation classes discriminated tumor from normal pleura (permutation P < 0.0001). In a random forests classification, the overall misclassification error rate was 3.4%, with <1% (n = 1) of tumors misclassified as normal (P < 0.0001). Among tumors, methylation class membership was significantly associated with lung tissue asbestos body burden (P < 0.03), and significantly predicted survival (likelihood ratio P < 0.01). Consistent with prior work, asbestos burden was associated with an increased risk of death (hazard ratio, 1.4; 95% confidence interval, 1.1-1.8). Our results have shown that methylation profiles powerfully differentiate diseased pleura from nontumor pleura and that asbestos burden and methylation profiles are independent predictors of mesothelioma patient survival. We have added to the growing body of evidence that cellular epigenetic dysregulation is a critical mode of action for asbestos in the induction of pleural mesothelioma. Importantly, these findings hold great promise for using epigenetic profiling in the diagnosis and prognosis of human cancers.

Variants near TERT and TERC influencing telomere length are associated with high-grade glioma risk

Glioma, the most common central nervous system cancer in adults, has poor prognosis. Here we identify a new SNP associated with glioma risk, rs1920116 (near TERC), that reached genome-wide significance (Pcombined = 8.3 × 10−9) in a meta-analysis of genome-wide association studies (GWAS) of high-grade glioma and replication data (1,644 cases and 7,736 controls). This region has previously been associated with mean leukocyte telomere length (LTL). We therefore examined the relationship between LTL and both this new risk locus and other previously established risk loci for glioma using data from a recent GWAS of LTL (n = 37,684 individuals). Alleles associated with glioma risk near TERC and TERT were strongly associated with longer LTL (P = 5.5 × 10−20 and 4.4 × 10−19, respectively). In contrast, risk-associated alleles near RTEL1 were inconsistently associated with LTL, suggesting the presence of distinct causal alleles. No other risk loci for glioma were associated with LTL. The identification of risk alleles for glioma near TERC and TERT that also associate with telomere length implicates telomerase in gliomagenesis.

Epigenetic profiling reveals etiologically distinct patterns of DNA methylation in head and neck squamous cell carcinoma

Head and neck squamous cell carcinomas (HNSCCs) represent clinically and etiologically heterogeneous tumors affecting >40 000 patients per year in the USA. Previous research has identified individual epigenetic alterations and, in some cases, the relationship of these alterations with carcinogen exposure or patient outcomes, suggesting that specific exposures give rise to specific types of molecular alterations in HNSCCs. Here, we describe how different etiologic factors are reflected in the molecular character and clinical outcome of these tumors. In a case series of primary, incident HNSCC (n = 68), we examined the DNA methylation profile of 1413 autosomal CpG loci in 773 genes, in relation to exposures and etiologic factors. The overall pattern of epigenetic alteration could significantly distinguish tumor from normal head and neck epithelial tissues (P < 0.0001) more effectively than specific gene methylation events. Among tumors, there were significant associations between specific DNA methylation profile classes and tobacco smoking and alcohol exposures. Although there was a significant association between methylation profile and tumor stage (P < 0.01), we did not observe an association between these profiles and overall patient survival after adjustment for stage; although methylation of a number of specific loci falling in different cellular pathways was associated with overall patient survival. We found that the etiologic heterogeneity of HNSCC is reflected in specific patterns of molecular epigenetic alterations within the tumors and that the DNA methylation profiles may hold clinical promise worthy of further study.

Genetic Polymorphisms in MTHFR 677 and 1298, GSTM1 and T1, and Metabolism of Arsenic

Methylation is the primary route of metabolism of inorganic arsenic in humans, and previous studies showed that interindividual differences in arsenic methylation may have important impacts on susceptibility to arsenic-induced cancer. To date, the factors that regulate arsenic methylation in humans are mostly unknown. Urinary arsenic methylation patterns and genetic polymorphisms in methylenetetrahydrofolate reductase (MTHFR) and glutathione S-transferase (GST) were investigated in 170 subjects from an arsenic-exposed region in Argentina. Previous studies showed that subjects with the TT/AA polymorphisms at MTHFR 677 and 1298 have lower MTHFR activity than others. In this study, it was found that subjects with the TT/AA variant of MTHFR 677/1298 excreted a significantly higher proportion of ingested arsenic as inorganic arsenic and a lower proportion as dimethylarsinic acid. Women with the null genotype of GSTM1 excreted a significantly higher proportion of arsenic as monomethylarsonate than women with the active genotype. No associations were seen between polymorphisms in GSTT1 and arsenic methylation. This is the first study to report (1) associations between MTHFR and arsenic metabolism in humans, and (2) gender differences between genetic polymorphisms and urinary arsenic methylation patterns. Overall, this study provides evidence that MTHFR and GSTM1 are involved in arsenic metabolism in humans, and polymorphisms in the genes that encode these enzymes may play a role in susceptibility to arsenic-induced cancer.