Laura S. Rozek

Genome-wide association scan identifies a colorectal cancer susceptibility locus on 11q23 and replicates risk loci at 8q24 and 18q21.

In a genome-wide association study to identify loci associated with colorectal cancer (CRC) risk, we genotyped 555,510 SNPs in 1,012 early-onset Scottish CRC cases and 1,012 controls (phase 1). In phase 2, we genotyped the 15,008 highest-ranked SNPs in 2,057 Scottish cases and 2,111 controls. We then genotyped the five highest-ranked SNPs from the joint phase 1 and 2 analysis in 14,500 cases and 13,294 controls from seven populations, and identified a previously unreported association, rs3802842 on 11q23 (OR = 1.1; P = 5.8 × 10−10), showing population differences in risk. We also replicated and fine-mapped associations at 8q24 (rs7014346; OR = 1.19; P = 8.6 × 10−26) and 18q21 (rs4939827; OR = 1.2; P = 7.8 × 10−28). Risk was greater for rectal than for colon cancer for rs3802842 (P < 0.008) and rs4939827 (P < 0.009). Carrying all six possible risk alleles yielded OR = 2.6 (95% CI = 1.75–3.89) for CRC. These findings extend our understanding of the role of common genetic variation in CRC etiology.

Genome-Wide DNA Methylation Differences Between Late-Onset Alzheimer's Disease and Cognitively Normal Controls in Human Frontal Cortex

Evidence supports a role for epigenetic mechanisms in the pathogenesis of late-onset Alzheimers disease (LOAD), but little has been done on a genome-wide scale to identify potential sites involved in disease. This study investigates human postmortem frontal cortex genome-wide DNA methylation profiles between 12 LOAD and 12 cognitively normal age- and gender-matched subjects. Quantitative DNA methylation is determined at 27,578 CpG sites spanning 14,475 genes via the Illumina Infinium HumanMethylation27 BeadArray. Data are analyzed using parallel linear models adjusting for age and gender with empirical Bayes standard error methods. Gene-specific technical and functional validation is performed on an additional 13 matched pair samples, encompassing a wider age range. Analysis reveals 948 CpG sites representing 918 unique genes as potentially associated with LOAD disease status pending confirmation in additional study populations. Across these 948 sites the subtle mean methylation difference between cases and controls is 2.9%. The CpG site with a minimum false discovery rate located in the promoter of the gene Transmembrane Protein 59 (TMEM59) is 7.3% hypomethylated in cases. Methylation at this site is functionally associated with tissue RNA and protein levels of the TMEM59 gene product. The TMEM59 gene identified from our discovery approach was recently implicated in amyloid-β protein precursor post-translational processing, supporting a role for epigenetic change in LOAD pathology. This study demonstrates widespread, modest discordant DNA methylation in LOAD-diseased tissue independent from DNA methylation changes with age. Identification of epigenetic biomarkers of LOAD risk may allow for the development of novel diagnostic and therapeutic targets.

Genome-wide methylation and expression differences in HPV(+) and HPV(-) squamous cell carcinoma cell lines are consistent with divergent mechanisms of carcinogenesis

Oncogenic human papillomaviruses (HPV) are associated with nearly all cervical cancers and are increasingly important in the etiology of oropharyngeal tumors. HPV-associated head and neck squamous cell carcinomas (HNSCC) have distinct risk profiles and appreciate a prognostic advantage compared to HPV-negative HNSCC. Promoter hypermethylation is widely recognized as a mechanism in the progression of HNSCC, but the extent to which this mechanism is consistent between HPV(+) and HPV(-) tumors is unknown. To investigate the epigenetic regulation of gene expression in HPV-induced and carcinogen-induced cancers, we examined genome-wide DNA methylation and gene expression in HPV(+) and HPV(-) SCC cell lines. We used two platforms: the Illumina Infinium Methylation BeadArray and tiling arrays, and confirmed illustrative examples with pyrosequencing and quantitative PCR. These analyses indicate that HPV(+) cell lines have higher DNA methylation in genic and LINE-1 regions than HPV(-) cell lines. Differentially methylated loci between HPV(+) and HPV(-) cell lines significantly correlated with HPV-typed HNSCC primary tumor DNA methylation levels. Novel findings include higher promoter methylation of polycomb repressive complex 2 target genes in HPV(+) cells compared to HPV(-) cells and increased expression of DNMT3A in HPV(+) cells. Additionally, CDKN2A and KRT8 were identified as interaction hubs among genes with higher methylation and lower expression in HPV(-) cells. Conversely, RUNX2, IRS-1 and CCNA1 were major hubs with higher methylation and lower expression in HPV(+) cells. Distinct HPV(+) and HPV(-) epigenetic profiles should provide clues to novel targets for development of individualized therapeutic strategies.

Epigenetic Responses Following Maternal Dietary Exposure to Physiologically Relevant Levels of Bisphenol A

Animal studies have linked perinatal bisphenol A (BPA) exposure to altered DNA methylation, but little attention is given to analyzing multiple physiologically relevant doses. Utilizing the viable yellow agouti (Avy) mouse, we examine the effects of developmental exposure through maternal diet to 50 ng BPA/kg (n = 14 litters), 50 μg BPA/kg (n = 9 litters), or 50 mg BPA/kg (n = 13 litters) on global and candidate gene methylation at postnatal day 22. Global methylation analysis reveals hypermethylation in tail tissue of a/a and Avy/a offspring across all dose groups compared with controls (n = 11 litters; P < 0.02). Analysis of coat color phenotype replicates previous work showing that the distribution of 50 mg BPA/kg Avy/a offspring shifts toward yellow (P = 0.006) by decreasing DNA methylation in the retrotransposon upstream of the Agouti gene (P = 0.03). Maternal exposure to 50 μg or 50 ng BPA/kg, however, results in altered coat color distributions in comparison with control (P = 0.04 and 0.02), but no DNA methylation effects at the Agouti gene are noted. DNA methylation at the CDK5 activator‐binding protein (CabpIAP) metastable epiallele shows hypermethylation in the 50 μg BPA/kg offspring, compared with controls (P = 0.02). Comparison of exposed mouse liver BPA levels to human fetal liver BPA levels indicates that the three experimental exposures are physiologically relevant. Thus, perinatal BPA exposure affects offspring phenotype and epigenetic regulation across multiple doses, indicating the need to evaluate dose effects in human clinical and population studies. Environ. Mol. Mutagen. 2012.

The MLH1 D132H variant is associated with susceptibility to sporadic colorectal cancer

Most susceptibility to colorectal cancer (CRC) is not accounted for by known risk factors. Because MLH1, MSH2 and MSH6 mutations underlie high-penetrance CRC susceptibility in hereditary nonpolyposis colon cancer (HNPCC), we hypothesized that attenuated alleles might also underlie susceptibility to sporadic CRC. We looked for gene variants associated with HNPCC in Israeli probands with familial CRC unstratified with respect to the microsatellite instability (MSI) phenotype. Association studies identified a new MLH1 variant (415G→C, resulting in the amino acid substitution D132H) in ∼1.3% of Israeli individuals with CRC self-described as Jewish, Christian and Muslim. MLH1 415C confers clinically significant susceptibility to CRC. In contrast to classic HNPCC, CRCs associated with MLH1 415C usually do not have the MSI defect, which is important for clinical mutation screening. Structural and functional analyses showed that the normal ATPase function of MLH1 is attenuated, but not eliminated, by the MLH1 415G→C mutation. The new MLH1 variant confers a high risk of CRC and identifies a previously unrecognized mechanism in microsatellite-stable tumors. These studies suggest that variants of mismatch repair proteins with attenuated function may account for a higher proportion of susceptibility to sporadic microsatellite-stable CRC than previously assumed.

Accurate classification of MLH1/MSH2 missense variants with multivariate analysis of protein polymorphisms-mismatch repair (MAPP-MMR).

Lynch syndrome, also known as hereditary nonpolyposis colon cancer (HNPCC), is the most common known genetic syndrome for colorectal cancer (CRC). MLH1/MSH2 mutations underlie approximately 90% of Lynch syndrome families. A total of 24% of these mutations are missense. Interpreting missense variation is extremely challenging. We have therefore developed multivariate analysis of protein polymorphisms–mismatch repair (MAPP‐MMR), a bioinformatic algorithm that effectively classifies MLH1/MSH2 deleterious and neutral missense variants. We compiled a large database (n>300) of MLH1/MSH2 missense variants with associated clinical and molecular characteristics. We divided this database into nonoverlapping training and validation sets and tested MAPP‐MMR. MAPP‐MMR significantly outperformed other missense variant classification algorithms (sensitivity, 94%; specificity, 96%; positive predictive value [PPV] 98%; negative predictive value [NPV], 89%), such as SIFT and PolyPhen. MAPP‐MMR is an effective bioinformatic tool for missense variant interpretation that accurately distinguishes MLH1/MSH2 deleterious variants from neutral variants. Hum Mutat 29(6), 852–860, 2008.

Aberrant DNA methylation in hereditary nonpolyposis colorectal cancer without mismatch repair deficiency.

BACKGROUND & AIMS Approximately half of the families that fulfill Amsterdam criteria for Lynch syndrome or hereditary nonpolyposis colorectal cancer (HNPCC) do not have evidence of the germline mismatch repair gene mutations that define this syndrome and result in microsatellite instability (MSI). The carcinogenic pathways and the best diagnostic approaches to detect microsatellite stable (MSS) HNPCC tumors are unclear. We investigated the contribution of epigenetic alterations to the development of MSS HNPCC tumors. METHODS Colorectal cancers were divided into 4 groups: (1) microsatellite stable, Amsterdam-positive (MSS HNPCC) (N = 22); (2) Lynch syndrome cancers (identified mismatch repair mutations) (N = 21); (3) sporadic MSS (N = 92); and (4) sporadic MSI (N = 46). Methylation status was evaluated for CACNAG1, SOCS1, RUNX3, NEUROG1, MLH1, and long interspersed nucleotide element-1 (LINE-1). KRAS and BRAF mutation status was analyzed. RESULTS MSS HNPCC tumors displayed a significantly lower degree of LINE-1 methylation, a marker for global methylation, than any other group. Although most MSS HNPCC tumors had some degree of CpG island methylation, none presented a high index of methylation. MSS HNPCC tumors had KRAS mutations exclusively in codon 12, but none harbored V600E BRAF mutations. CONCLUSIONS Tumors from Amsterdam-positive patients without mismatch repair deficiency (MSS HNPCC) have certain molecular features, including global hypomethylation, that distinguish them from all other colorectal cancers. These characteristics could have an important impact on tumor behavior or treatment response. Studies are underway to further assess the cause and effects of these features.

MethylSig: a whole genome DNA methylation analysis pipeline

MOTIVATION DNA methylation plays critical roles in gene regulation and cellular specification without altering DNA sequences. The wide application of reduced representation bisulfite sequencing (RRBS) and whole genome bisulfite sequencing (bis-seq) opens the door to study DNA methylation at single CpG site resolution. One challenging question is how best to test for significant methylation differences between groups of biological samples in order to minimize false positive findings. RESULTS We present a statistical analysis package, methylSig, to analyse genome-wide methylation differences between samples from different treatments or disease groups. MethylSig takes into account both read coverage and biological variation by utilizing a beta-binomial approach across biological samples for a CpG site or region, and identifies relevant differences in CpG methylation. It can also incorporate local information to improve group methylation level and/or variance estimation for experiments with small sample size. A permutation study based on data from enhanced RRBS samples shows that methylSig maintains a well-calibrated type-I error when the number of samples is three or more per group. Our simulations show that methylSig has higher sensitivity compared with several alternative methods. The use of methylSig is illustrated with a comparison of different subtypes of acute leukemia and normal bone marrow samples. AVAILABILITY methylSig is available as an R package at http://sartorlab.ccmb.med.umich.edu/software. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

Alzheimer’s Disease and Environmental Exposure to Lead: The Epidemiologic Evidence and Potential Role of Epigenetics

Several lines of evidence indicate that the etiology of late-onset Alzheimers disease (LOAD) is complex, with significant contributions from both genes and environmental factors. Recent research suggests the importance of epigenetic mechanisms in defining the relationship between environmental exposures and LOAD. In epidemiologic studies of adults, cumulative lifetime lead (Pb) exposure has been associated with accelerated declines in cognition. In addition, research in animal models suggests a causal association between Pb exposure during early life, epigenetics, and LOAD. There are multiple challenges to human epidemiologic research evaluating the relationship between epigenetics, LOAD, and Pb exposure. Epidemiologic studies are not well-suited to accommodate the long latency period between exposures during early life and onset of Alzheimers disease. There is also a lack of validated circulating epigenetics biomarkers and retrospective biomarkers of Pb exposure. Members of our research group have shown bone Pb is an accurate measurement of historical Pb exposure in adults, offering an avenue for future epidemiologic studies. However, this would not address the risk of LOAD attributable to early-life Pb exposures. Future studies that use a cohort design to measure both Pb exposure and validated epigenetic biomarkers of LOAD will be useful to clarify this important relationship.

Genetic variation in 8q24 associated with risk of colorectal cancer

Chromosome 8q24 harbors oncogenes known to be involved in pathogenesis of colorectal cancer (CRC) as well as uncharacterized genetic variants that have recently been shown to influence inherited risk of prostate cancer. In a population-based case-control study of colorectal cancer in northern Israel, we investigated the association between variation in 8q24 and risk of CRC. Among 1,861 incident cases and 1,937 population-based controls matched on age, gender, ethnicity, and clinic, rs10505477 was associated with risk of CRC in a dominant model, with an odds ratio = 1.23, 95% confidence interval = 1.05 -1.43, (p=0.008). This association was independently validated in an analysis of cancer among relatives of carriers of the risk allele, with a hazard ratio of 3.2 (95% bootstrap CI = 1.16 – 17.8). Genetic variation at rs10505477 on 8q24 potentially accounts for 14% of CRC in this population and should be replicated in other studies.