Roadmap for investigating epigenome deregulation and environmental origins of cancer.

The interaction between the (epi)genetic makeup of an individual and his/her environmental exposure record (exposome) is accepted as a determinant factor for a significant proportion of human malignancies. Recent evidence has highlighted the key role of epigenetic mechanisms in mediating gene-environment interactions and translating exposures into tumorigenesis. There is also growing evidence that epigenetic changes may be risk factor-specific ("fingerprints") that should prove instrumental in the discovery of new biomarkers in cancer. Here, we review the state of the science of epigenetics associated with environmental stimuli and cancer risk, highlighting key developments in the field. Critical knowledge gaps and research needs are discussed and advances in epigenomics that may help in understanding the functional relevance of epigenetic alterations. Key elements required for causality inferences linking epigenetic changes to exposure and cancer are discussed and how these alterations can be incorporated in carcinogen evaluation and in understanding mechanisms underlying epigenome deregulation by the environment.

Childhood Leukemia: A Preventable Disease.

In contrast to most pediatric cancers, there is a growing body of literature, nationally and internationally, that has implicated the role of several environmental indoor and outdoor hazards in the etiology of childhood leukemia. For example, exposures to solvents, traffic, pesticides, and tobacco smoke have consistently demonstrated positive associations with the risk of developing childhood leukemia. Intake of vitamins and folate supplementation during the preconception period or pregnancy has been demonstrated to have a protective effect. Despite the strength of these findings, the dissemination of this knowledge to clinicians has been limited. Some children may be more vulnerable than others as documented by the high and increasing incidence of childhood leukemia in Hispanics. To protect children's health, it is prudent to establish programs to alter exposure to those factors with well-established associations with leukemia risk rather than to suspend judgment until no uncertainty remains. This is particularly true because other serious health outcomes (both negative and positive) have been associated with the same exposures. We draw from historical examples to put in perspective the arguments of association versus causation, as well as to discuss benefits versus risks of immediate and long-term preventive actions.

Genetic variants in telomerase-related genes are associated with an older age at diagnosis in glioma patients: evidence for distinct pathways of gliomagenesis.

BACKGROUND: Genome-wide association studies have implicated single nucleotide polymorphisms (SNPs) in 7 genes as glioma risk factors, including 2 (TERT, RTEL1) involved in telomerase structure/function. We examined associations of these 7 established glioma risk loci with age at diagnosis among patients with glioma. METHODS: SNP genotype data were available for 2286 Caucasian glioma patients from the University of California, San Francisco (n = 1434) and the Mayo Clinic (n = 852). Regression analyses were performed to test for associations between "number of risk alleles" and "age at diagnosis," adjusted for sex and study site and stratified by tumor grade/histology where appropriate. RESULTS: Four SNPs were significantly associated with age at diagnosis. Carrying a greater number of risk alleles at rs55705857 (CCDC26) and at rs498872 (PHLDB1) was associated with younger age at diagnosis (P = 1.4 × 10(-22) and P = 9.5 × 10(-7), respectively). These SNPs are stronger risk factors for oligodendroglial tumors, which tend to occur in younger patients, and their association with age at diagnosis varied across tumor subtypes. In contrast, carrying more risk alleles at rs2736100 (TERT) and at rs6010620 (RTEL1) was associated with older age at diagnosis (P = 6.2 × 10(-4) and P = 2.5 × 10(-4), respectively). These SNPs are risk factors for all glioma grades/histologies, and their association with age at diagnosis was consistent across tumor subgroups. CONCLUSIONS: Carrying a greater number of risk alleles might be expected to decrease age at diagnosis. However, glioma susceptibility conferred by variation in telomerase-related genes did not follow this pattern. This supports the hypothesis that telomerase-related mechanisms of telomere maintenance are more associated with gliomas that develop later in life than those utilizing telomerase-independent mechanisms (ie, alternative lengthening of telomeres).

Inherited variant on chromosome 11q23 increases susceptibility to IDH-mutated but not IDH-normal gliomas regardless of grade or histology.

INTRODUCTION: Recent discoveries of inherited glioma risk loci and acquired IDH mutations are providing new insights into glioma etiology. IDH mutations are common in lower grade gliomas and secondary glioblastomas and uncommon in primary glioblastomas. Because the inherited variant in 11q23 has been associated with risk of lower grade glioma and not with glioblastomas, we hypothesized that this variant increases susceptibility to IDH-mutated gliomas, but not to IDH-wild-type gliomas. METHODS: We tested this hypothesis in patients with glioma and controls from the San Francisco Adult Glioma Study, the Mayo Clinic, and Illumina controls (1102 total patients, 5299 total controls). Case-control additive associations of 11q23 risk alleles (rs498872, T allele) were calculated using logistic regression, stratified by tumor IDH status (mutated or wild-type) and by histology and grade. We also adjusted for the recently discovered 8q24 glioma risk locus rs55705857 G allele. RESULTS: The 11q23 glioma risk locus was associated with increased risk of IDH-mutated gliomas of all histologies and grades (odds ratio [OR] = 1.50; 95% confidence interval [CI] = 1.29-1.74; P = 1.3X10(-7)) but not with IDH-wild-type gliomas of any histology or grade (OR = 0.91; 95% CI = 0.81-1.03; P = 0.14). The associations were independent of the rs55705857 G allele. CONCLUSION: A variant at the 11q23 locus increases risk for IDH-mutated but not IDH-wild-type gliomas, regardless of grade or histology.

Analysis of 60 reported glioma risk SNPs replicates published GWAS findings but fails to replicate associations from published candidate-gene studies.

Genomewide association studies (GWAS) and candidate-gene studies have implicated single-nucleotide polymorphisms (SNPs) in at least 45 different genes as putative glioma risk factors. Attempts to validate these associations have yielded variable results and few genetic risk factors have been consistently replicated. We conducted a case-control study of Caucasian glioma cases and controls from the University of California San Francisco (810 cases, 512 controls) and the Mayo Clinic (852 cases, 789 controls) in an attempt to replicate previously reported genetic risk factors for glioma. Sixty SNPs selected from the literature (eight from GWAS and 52 from candidate-gene studies) were successfully genotyped on an Illumina custom genotyping panel. Eight SNPs in/near seven different genes (TERT, EGFR, CCDC26, CDKN2A, PHLDB1, RTEL1, TP53) were significantly associated with glioma risk in the combined dataset (P < 0.05), with all associations in the same direction as in previous reports. Several SNP associations showed considerable differences across histologic subtype. All eight successfully replicated associations were first identified by GWAS, although none of the putative risk SNPs from candidate-gene studies was associated in the full case-control sample (all P values > 0.05). Although several confirmed associations are located near genes long known to be involved in gliomagenesis (e.g., EGFR, CDKN2A, TP53), these associations were first discovered by the GWAS approach and are in noncoding regions. These results highlight that the deficiencies of the candidate-gene approach lay in selecting both appropriate genes and relevant SNPs within these genes.

A low-frequency variant at 8q24.21 is strongly associated with risk of oligodendroglial tumors and astrocytomas with IDH1 or IDH2 mutation.

Variants at 8q24.21 have been shown to be associated with glioma development. By means of tag SNP genotyping and imputation, pooled next-generation sequencing using long-range PCR and subsequent validation SNP genotyping, we identified seven low-frequency SNPs at 8q24.21 that were strongly associated with glioma risk (P=1×10(-25) to 1×10(-14)). The most strongly associated SNP, rs55705857, remained highly significant after individual adjustment for the other top six SNPs and two previously published SNPs. After stratifying by histological and tumor genetic subtype, the most significant associations of rs55705857 were with oligodendroglial tumors and gliomas with mutant IDH1 or IDH2 (odds ratio (OR)=5.1, P=1.1×10(-31) and OR=4.8, P=6.6×10(-22), respectively). Strong associations were observed for astrocytomas with mutated IDH1 or IDH2 (grades 2-4) (OR=5.16-6.66, P=4.7×10(-12) to 2.2×10(-8)) but not for astrocytomas with wild-type IDH1 and IDH2 (smallest P=0.26). The conserved sequence block that includes rs55705857 is consistently modeled as a microRNA.

SSBP2 variants are associated with survival in glioblastoma patients.

PURPOSE: Glioblastoma is a devastating, incurable disease with few known prognostic factors. Here, we present the first genome-wide survival and validation study for glioblastoma. EXPERIMENTAL DESIGN: Cox regressions for survival with 314,635 inherited autosomal single-nucleotide polymorphisms (SNP) among 315 San Francisco Adult Glioma Study patients for discovery and three independent validation data sets [87 Mayo Clinic, 232 glioma patients recruited from several medical centers in Southeastern United States (GliomaSE), and 115 The Cancer Genome Atlas patients] were used to identify SNPs associated with overall survival for Caucasian glioblastoma patients treated with the current standard of care, resection, radiation, and temozolomide (total n = 749). Tumor expression of the gene that contained the identified prognostic SNP was examined in three separate data sets (total n = 619). Genotype imputation was used to estimate hazard ratios (HR) for SNPs that had not been directly genotyped. RESULTS: From the discovery and validation analyses, we identified a variant in single-stranded DNA-binding protein 2 (SSBP2) on 5q14.1 associated with overall survival in combined analyses (HR, 1.64; P = 1.3 × 10(-6)). Expression of SSBP2 in tumors from three independent data sets also was significantly related to patient survival (P = 5.3 × 10(-4)). Using genotype imputation, the SSBP2 SNP rs17296479 had the strongest statistically significant genome-wide association with poorer overall patient survival (HR, 1.79; 95% CI, 1.45-2.22; P = 1.0 × 10(-7)). CONCLUSION: The minor allele of SSBP2 SNP rs17296479 and the increased tumor expression of SSBP2 were statistically significantly associated with poorer overall survival among glioblastoma patients. With further confirmation, previously unrecognized inherited variations influencing survival may warrant inclusion in clinical trials to improve randomization. Unaccounted for genetic influence on survival could produce unwanted bias in such studies.

Inherited variation in immune genes and pathways and glioblastoma risk.

To determine whether inherited variations in immune function single-nucleotide polymorphisms (SNPs), genes or pathways affect glioblastoma risk, we analyzed data from recent genome-wide association studies in conjunction with predefined immune function genes and pathways. Gene and pathway analyses were conducted on two independent data sets using 6629 SNPs in 911 genes on 17 immune pathways from 525 glioblastoma cases and 602 controls from the University of California, San Francisco (UCSF) and a subset of 6029 SNPs in 893 genes from 531 cases and 1782 controls from MD Anderson (MDA). To further assess consistency of SNP-level associations, we also compared data from the UK (266 cases and 2482 controls) and the Mayo Clinic (114 cases and 111 controls). Although three correlated epidermal growth factor receptor (EGFR) SNPs were consistently associated with glioblastoma in all four data sets (Mantel-Haenzel P values = 1 × 10⁻5 to 4 × 10⁻³), independent replication is required as genome-wide significance was not attained. In gene-level analyses, eight immune function genes were significantly (minP < 0.05) associated with glioblastoma; the IL-2RA (CD25) cytokine gene had the smallest minP values in both UCSF (minP = 0.01) and MDA (minP = 0.001) data sets. The IL-2RA receptor is found on the surface of regulatory T cells potentially contributing to immunosuppression characteristic of the glioblastoma microenvironment. In pathway correlation analyses, cytokine signaling and adhesion-extravasation-migration pathways showed similar associations with glioblastoma risk in both MDA and UCSF data sets. Our findings represent the first systematic description of immune genes and pathways that characterize glioblastoma risk.

Variants in the CDKN2B and RTEL1 regions are associated with high-grade glioma susceptibility.

The causes of glioblastoma and other gliomas remain obscure. To discover new candidate genes influencing glioma susceptibility, we conducted a principal component-adjusted genome-wide association study (GWAS) of 275,895 autosomal variants among 692 adult high-grade glioma cases (622 from the San Francisco Adult Glioma Study (AGS) and 70 from the Cancer Genome Atlas (TCGA)) and 3,992 controls (602 from AGS and 3,390 from Illumina iControlDB (iControls)). For replication, we analyzed the 13 SNPs with P < 10(-6) using independent data from 176 high-grade glioma cases and 174 controls from the Mayo Clinic. On 9p21, rs1412829 near CDKN2B had discovery P = 3.4 x 10(-8), replication P = 0.0038 and combined P = 1.85 x 10(-10). On 20q13.3, rs6010620 intronic to RTEL1 had discovery P = 1.5 x 10(-7), replication P = 0.00035 and combined P = 3.40 x 10(-9). For both SNPs, the direction of association was the same in discovery and replication phases.

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.

Nonsynonymous coding single-nucleotide polymorphisms spanning the genome in relation to glioblastoma survival and age at diagnosis.

PURPOSE: Our aim was to discover possible inherited factors associated with glioblastoma age at diagnosis and survival. Although new genotyping technologies allow greatly expanded exploration of such factors, they pose many challenges. EXPERIMENTAL DESIGN: In this pilot study, we (a) genotyped 112 newly diagnosed glioblastoma patients ascertained through a population-based study (group 1) with the ParAllele assay panel of approximately 10,000 nonsynonymous coding single-nucleotide polymorphisms (SNP), (b) used several statistical and bioinformatic techniques to identify 17 SNPs potentially related to either glioblastoma age at diagnosis or survival, and (c) genotyped 16 of these SNPs using conventional PCR methods in an independent group of 195 glioblastoma patients (group 2). RESULTS: In group 2, only one of the 16 SNPs, rs8057643 (located on 16p13.2), was significantly associated with glioblastoma age at diagnosis (nominal P = 0.0017; Bonferroni corrected P = 0.054). Median ages at diagnosis for those with 0, 1, or 2 T alleles were 66, 57, and 59 years in group 1 and 64, 57, and 55 years in group 2 (combined P = 0.001). Furthermore, Cox regression analyses of time to death with number of T alleles adjusted for gender and patient group yielded a hazard ratio of 0.82 (95% confidence interval, 0.68-0.98; P = 0.03). CONCLUSIONS: Although limited by a relatively small sample size, this pilot study, using well-characterized, unambiguous disease characteristics, illustrates the necessity of independent replication owing to the likelihood of false positives. Several other challenges are discussed, including attempts to incorporate information on the potential functional importance of SNPs in genome-disease association studies.

Serum IgE, tumor epidermal growth factor receptor expression, and inherited polymorphisms associated with glioma survival.

In population-based glioma patients, we examined survival in relation to potentially pertinent constitutive polymorphisms, serologic factors, and tumor genetic and protein alterations in epidermal growth factor receptor (EGFR), MDM2, and TP53. Subjects were newly diagnosed adults residing in the San Francisco Bay Surveillance Epidemiology and End Results Area during 1991 to 1994 and 1997 to 1999 with central neuropathology review (n = 873). Subjects provided blood for serologic studies of IgE and IgG to four herpes viruses and constitutive specimens for genotyping 22 polymorphisms in 13 genes (n = 471). We obtained 595 of 697 astrocytic tumors for marker studies. We determined treatments, vital status, and other factors using registry, interview, medical record, and active follow-up data. Cox regressions for survival were adjusted for age, gender, ethnicity, study series, resection versus biopsy only, radiation, and chemotherapy. Using a stringent P < 0.001, glioma survival was associated with ERCC1 C8092A [hazard ratio (HR), 0.72; 95% confidence limits (95% CL), 0.60-0.86; P = 0.0004] and GSTT1 deletion (HR, 1.64; 95% CL, 1.25-2.16; P = 0.0004); glioblastoma patients with elevated IgE had 9 months longer survival than those with normal or borderline IgE levels (HR, 0.62; 95% CL, 0.47-0.82; P = 0.0007), and EGFR expression in anaplastic astrocytoma was associated with nearly 3-fold poorer survival (HR, 2.97; 95% CL, 1.70-5.19; P = 0.0001). Based on our and others' findings, we recommend further studies to (a) understand relationships of elevated IgE levels and other immunologic factors with improved glioblastoma survival potentially relevant to immunologic therapies and (b) determine which inherited ERCC1 variants or other variants in the 19q13.3 region influence survival. We also suggest that tumor EGFR expression be incorporated into clinical evaluation of anaplastic astrocytoma patients.

Molecular features of adult glioma associated with patient race/ethnicity, age, and a polymorphism in O6-methylguanine-DNA-methyltransferase.

BACKGROUND: Risk factors for adult glioma in the San Francisco Bay Area include well-known demographic features such as age and race/ethnicity, and our previous studies indicated that these characteristics are associated with the TP53 mutation status of patients' tumors. We enlarged our study to assess the relationships of risk factors with TP53 as well as epidermal growth factor receptor (EGFR) and murine double minute-2 (MDM2) gene amplification and expression and the germ line Leu84Phe polymorphism in the DNA repair protein O6-methylguanine-DNA-methyltransferase (MGMT). MGMT expression may depend on the TP53 status of cells. METHODS: Molecular analyses were carried out on 556 incident astrocytic tumors. MGMT genotype data were collected on germ line DNA from 260 of these cases. RESULTS: The tumor data confirm the inverse relationships between TP53 mutation and MDM2 (P = 0.04) or EGFR (P = 0.004) amplification and that patients whose tumors contain TP53 mutations are younger than those without (P < 0.001). Although there was little difference in age of patient by EGFR amplification or expression among glioblastoma multiforme cases, EGFR gene amplification was associated with much older age of onset of anaplastic astrocytoma; for example, EGFR-amplified anaplastic astrocytoma cases were on average 63 years old compared with 48 years for nonamplified cases (P = 0.005). An increased prevalence of TP53 mutation positive glioblastoma multiforme was noted among nonwhites (African American and Asian) compared with whites (Latino and non-Latino; P = 0.004). Carriers of the MGMT variant 84Phe allele were significantly less likely to have tumors with TP53 overexpression (odds ratio, 0.30; 95% confidence interval, 0.13-0.71) and somewhat less likely to have tumors with any TP53 mutation (odds ratio, 0.47; 95% confidence interval, 0.13-1.69) after adjusting for age, gender, and ethnicity. Interestingly, EGFR gene amplification and EGFR protein overexpression were also inversely associated with the MGMT 84Phe allele. CONCLUSIONS: Our results are consistent with ethnic variation in glioma pathogenesis. The data on MGMT show that an inherited factor involving the repair of methylation and other alkylation damage, specifically to the O6 position of guanine, may be associated with the development of tumors that proceed in their development without TP53 mutations or accumulation of TP53 protein and possibly also those that do not involve amplification of the EGFR locus.

Origins of chromosome translocations in childhood leukaemia.

Chromosome translocations are often early or initiating events in leukaemogenesis, occurring prenatally in most cases of childhood leukaemia. Although these genetic changes are necessary, they are usually not sufficient to cause leukaemia. How, when and where do translocations arise? And can these insights aid our understanding of the natural history, pathogenesis and causes of leukaemia?