Devin C. Koestler

Review of processing and analysis methods for DNA methylation array data.

The promise of epigenome-wide association studies and cancer-specific somatic DNA methylation changes in improving our understanding of cancer, coupled with the decreasing cost and increasing coverage of DNA methylation microarrays, has brought about a surge in the use of these technologies. Here, we aim to provide both a review of issues encountered in the processing and analysis of array-based DNA methylation data and a summary of the advantages of recent approaches proposed for handling those issues, focusing on approaches publicly available in open-source environments such as R and Bioconductor. We hope that the processing tools and analysis flowchart described herein will facilitate researchers to effectively use these powerful DNA methylation array-based platforms, thereby advancing our understanding of human health and disease.

Associations between Gut Microbial Colonization in Early Life and Respiratory Outcomes in Cystic Fibrosis

OBJECTIVEnTo examine patterns of microbial colonization of the respiratory and intestinal tracts in early life in infants with cystic fibrosis (CF) and their associations with breastfeeding and clinical outcomes.nnnSTUDY DESIGNnA comprehensive, prospective longitudinal analysis of the upper respiratory and intestinal microbiota in a cohort of infants and young children with CF followed from birth was performed. Genus-level microbial community composition was characterized using 16S-targeted pyrosequencing, and relationships with exposures and outcomes were assessed using linear mixed-effects models, time-to-event analysis, and principal components analysis.nnnRESULTSnSequencing of 120 samples from 13 subjects collected from birth to 34 months revealed relationships between breastfeeding, microbial diversity in the respiratory and intestinal tracts, and the timing of onset of respiratory complications, including exacerbations and colonization with Pseudomonas aeruginosa. Fluctuations in the abundance of specific bacterial taxa preceded clinical outcomes, including a significant decrease in bacteria of the genus Parabacteroides within the intestinal tract prior to the onset of chronic P aeruginosa colonization. Specific assemblages of bacteria in intestinal samples, but not respiratory samples, were associated with CF exacerbation in early life, indicating that the intestinal microbiome may play a role in lung health.nnnCONCLUSIONSnOur findings relating breastfeeding to respiratory outcomes, gut diversity to prolonged periods of health, and specific bacterial communities in the gut prior to respiratory complications in CF highlight a connection between the intestinal microbiome and health and point to potential opportunities for antibiotic or probiotic interventions. Further studies in larger cohorts validating these findings are needed.

Differential DNA methylation in umbilical cord blood of infants exposed to mercury and arsenic in utero

Mercury and arsenic are known developmental toxicants. Prenatal exposures are associated with adverse childhood health outcomes that could be in part mediated by epigenetic alterations that may also contribute to altered immune profiles. In this study, we examined the association between prenatal mercury exposure on both DNA methylation and white blood cell composition of cord blood, and evaluated the interaction with prenatal arsenic exposure. A total of 138 mother-infant pairs with postpartum maternal toenail mercury, prenatal urinary arsenic concentrations, and newborn cord blood were assessed using the Illumina Infinium Methylation450 array. White blood cell composition was inferred from DNA methylation measurements. A doubling in toenail mercury concentration was associated with a 2.5% decrease (95% CI: 5.0%, 1.0%) in the estimated monocyte proportion. An increase of 3.5% (95% CI: 1.0, 7.0) in B-cell proportion was observed for females only. Among the top 100 CpGs associated with toenail mercury levels (ranked on P-value), there was a significant enrichment of loci located in North shore regions of CpG islands (P = 0.049), and the majority of these loci were hypermethylated (85%). Among the top 100 CpGs for the interaction between arsenic and mercury, there was a greater than expected proportion of loci located in CpG islands (P = 0.045) and in South shore regions (P = 0.009) and all of these loci were hypermethylated. This work supports the hypothesis that mercury may be contributing to epigenetic variability and immune cell proportion changes, and suggests that in utero exposure to mercury and arsenic, even at low levels, may interact to impact the epigenome.

DNA methylation profiles delineate etiologic heterogeneity and clinically important subgroups of bladder cancer

DNA methylation profiles can be used to define molecular cancer subtypes that may better inform disease etiology and clinical decision-making. This investigation aimed to create DNA methylation profiles of bladder cancer based on CpG methylation from almost 800 cancer-related genes and to then examine the relationship of those profiles with exposures related to risk and clinical characteristics. DNA, derived from formalin-fixed paraffin-embedded tumor samples obtained from incident cases involved in a population-based case-control study of bladder cancer in New Hampshire, was used for methylation profiling on the Illumina GoldenGate Methylation Bead Array. Unsupervised clustering of those loci with the greatest change in methylation between tumor and non-diseased tissue was performed to defined molecular subgroups of disease, and univariate tests of association followed by multinomial logistic regression was used to examine the association between these classes, bladder cancer risk factors and clinical phenotypes. Membership in the two most methylated classes was significantly associated with invasive disease (P < 0.001 for both class 3 and 4). Male gender (P = 0.04) and age >70 years (P = 0.05) was associated with membership in one of the most methylated classes. Finally, average water arsenic levels in the highest percentile predicted membership in an intermediately methylated class of tumors (P = 0.02 for both classes). Exposures and demographic associated with increased risk of bladder cancer specifically associate with particular subgroups of tumors defined by DNA methylation profiling and these subgroups may define more aggressive disease.

Correcting for cell-type heterogeneity in epigenome-wide association studies: Revisiting previous analyses

Correcting for cell-type heterogeneity in epigenome-wide association studies: revisiting previous analyses

Improving cell mixture deconvolution by identifying optimal DNA methylation libraries (IDOL).

BackgroundConfounding due to cellular heterogeneity represents one of the foremost challenges currently facing Epigenome-Wide Association Studies (EWAS). Statistical methods leveraging the tissue-specificity of DNA methylation for deconvoluting the cellular mixture of heterogenous biospecimens offer a promising solution, however the performance of such methods depends entirely on the library of methylation markers being used for deconvolution. Here, we introduce a novel algorithm for Identifying Optimal Libraries (IDOL) that dynamically scans a candidate set of cell-specific methylation markers to find libraries that optimize the accuracy of cell fraction estimates obtained from cell mixture deconvolution.ResultsApplication of IDOL to training set consisting of samples with both whole-blood DNA methylation data (Illumina HumanMethylation450 BeadArray (HM450)) and flow cytometry measurements of cell composition revealed an optimized library comprised of 300 CpG sites. When compared existing libraries, the library identified by IDOL demonstrated significantly better overall discrimination of the entire immune cell landscape (p = 0.038), and resulted in improved discrimination of 14 out of the 15 pairs of leukocyte subtypes. Estimates of cell composition across the samples in the training set using the IDOL library were highly correlated with their respective flow cytometry measurements, with all cell-specific R2>0.99 and root mean square errors (RMSEs) ranging from [0.97 % to 1.33 %] across leukocyte subtypes. Independent validation of the optimized IDOL library using two additional HM450 data sets showed similarly strong prediction performance, with all cell-specific R2>0.90 and RMSE<4.00 %. In simulation studies, adjustments for cell composition using the IDOL library resulted in uniformly lower false positive rates compared to competing libraries, while also demonstrating an improved capacity to explain epigenome-wide variation in DNA methylation within two large publicly available HM450 data sets.ConclusionsDespite consisting of half as many CpGs compared to existing libraries for whole blood mixture deconvolution, the optimized IDOL library identified herein resulted in outstanding prediction performance across all considered data sets and demonstrated potential to improve the operating characteristics of EWAS involving adjustments for cell distribution. In addition to providing the EWAS community with an optimized library for whole blood mixture deconvolution, our work establishes a systematic and generalizable framework for the assembly of libraries that improve the accuracy of cell mixture deconvolution.

Placental FKBP5 Genetic and Epigenetic Variation Is Associated with Infant Neurobehavioral Outcomes in the RICHS Cohort

Adverse maternal environments can lead to increased fetal exposure to maternal cortisol, which can cause infant neurobehavioral deficits. The placenta regulates fetal cortisol exposure and response, and placental DNA methylation can influence this function. FK506 binding protein (FKBP5) is a negative regulator of cortisol response, FKBP5 methylation has been linked to brain morphology and mental disorder risk, and genetic variation of FKBP5 was associated with post-traumatic stress disorder in adults. We hypothesized that placental FKBP5 methylation and genetic variation contribute to gene expression control, and are associated with infant neurodevelopmental outcomes assessed using the Neonatal Intensive Care Unit (NICU) Network Neurobehavioral Scales (NNNS). In 509 infants enrolled in the Rhode Island Child Health Study, placental FKBP5 methylation was measured at intron 7 using quantitative bisulfite pyrosequencing. Placental FKBP5 mRNA was measured in a subset of 61 infants by quantitative PCR, and the SNP rs1360780 was genotyped using a quantitative allelic discrimination assay. Relationships between methylation, expression and NNNS scores were examined using linear models adjusted for confounding variables, then logistic models were created to determine the influence of methylation on membership in high risk groups of infants. FKBP5 methylation was negatively associated with expression (Pu200a=u200a0.08, ru200a=u200a−0.22); infants with the TT genotype had higher expression than individuals with CC and CT genotypes (Pu200a=u200a0.06), and those with CC genotype displayed a negative relationship between methylation and expression (Pu200a=u200a0.06, ru200a=u200a−0.43). Infants in the highest quartile of FKBP5 methylation had increased risk of NNNS high arousal compared to infants in the lowest quartile (OR 2.22, CI 1.07–4.61). TT genotype infants had increased odds of high NNNS stress abstinence (OR 1.98, CI 0.92–4.26). Placental FKBP5 methylation reduces expression in a genotype specific fashion, and genetic variation supersedes this effect. These genetic and epigenetic differences in expression may alter the placenta’s ability to modulate cortisol response and exposure, leading to altered neurobehavioral outcomes.

Integrative genomic analysis identifies epigenetic marks that mediate genetic risk for epithelial ovarian cancer

BackgroundBoth genetic and epigenetic factors influence the development and progression of epithelial ovarian cancer (EOC). However, there is an incomplete understanding of the interrelationship between these factors and the extent to which they interact to impact disease risk. In the present study, we aimed to gain insight into this relationship by identifying DNA methylation marks that are candidate mediators of ovarian cancer genetic risk.MethodsWe used 214 cases and 214 age-matched controls from the Mayo Clinic Ovarian Cancer Study. Pretreatment, blood-derived DNA was profiled for genome-wide methylation (Illumina Infinium HumanMethylation27 BeadArray) and single nucleotide polymorphisms (SNPs, Illumina Infinium HD Human610-Quad BeadArray). The Causal Inference Test (CIT) was implemented to distinguish CpG sites that mediate genetic risk, from those that are consequential or independently acted on by genotype.ResultsControlling for the estimated distribution of immune cells and other key covariates, our initial epigenome-wide association analysis revealed 1,993 significantly differentially methylated CpGs that between cases and controls (FDR, qu2009<u20090.05). The relationship between methylation and case-control status for these 1,993 CpGs was found to be highly consistent with the results of previously published, independent study that consisted of peripheral blood DNA methylation signatures in 131 pretreatment cases and 274 controls. Implementation of the CIT test revealed 17 CpG/SNP pairs, comprising 13 unique CpGs and 17 unique SNPs, which represent potential methylation-mediated relationships between genotype and EOC risk. Of these 13 CpGs, several are associated with immune related genes and genes that have been previously shown to exhibit altered expression in the context of cancer.ConclusionsThese findings provide additional insight into EOC etiology and may serve as novel biomarkers for EOC susceptibility.

DNA methylation changes in the placenta are associated with fetal manganese exposure

Adequate micronutrient intake, including manganese (Mn), is important for fetal development. Both Mn deficiencies and excess exposures are associated with later-life disease, and Mn accumulates in the placenta. Placental functional alterations may alter fetal programming and lifelong health, and we hypothesized that prenatal exposures to Mn may alter placental function through epigenetic mechanisms. Using Illuminas HumanMethylation450 BeadArray, DNA methylation of >485,000 CpG loci genome-wide was interrogated in 61 placental samples and Mn associations assessed genome-wide via omnibus test (p=0.045). 713 loci were associated with Mn exposure (p<0.0001). Five significantly differentially-methylated (p<1.3×10(-7)) loci reside in neurodevelopmental, fetal growth and cancer-related genes. cg22284422, within the uncharacterized LOC284276 gene, was associated with birth weight; for every 10% increase in methylation, lower birth weights were observed, with an average decrease of 293.44g. Our observations suggest a link between prenatal micronutrient levels, placental epigenetic status and birth weight, although these preliminary results require validation.

DNA methylation in ductal carcinoma in situ related with future development of invasive breast cancer

BackgroundDuctal carcinoma in situ (DCIS) is a heterogeneous, pre-invasive lesion associated with an increased risk for future invasive ductal carcinoma. However, accurate risk stratification for development of invasive disease and appropriate treatment decisions remain clinical challenges. DNA methylation alterations are early events in the progression of cancer and represent emerging molecular markers that may predict invasive recurrence more accurately than traditional measures of DCIS prognosis.ResultsWe measured DNA methylation using the Illumina HumanMethylation450K array of estrogen-receptor positive DCIS (nu2009=u200940) and adjacent-normal (nu2009=u200915) tissues from subjects in the New Hampshire Mammography Network longitudinal breast imaging registry. We identified locus-specific methylation differences between DCIS and matched adjacent-normal tissue (95,609 CpGs, Qu2009<u20090.05). Among 40 DCIS cases, 13 later developed invasive disease and we identified 641 CpG sites that exhibited differential DNA methylation (Pu2009<u20090.01 and median |∆β|u2009>u20090.1) in these cases compared with age-matched subjects without invasive disease. The set of differentially methylated CpG loci associated with disease progression was enriched in homeobox-containing genes (Pu2009=u20091.3E-09) and genes involved with limb morphogenesis (Pu2009=u20091.0E-05). In an independent cohort, a subset of genes with progression-related differential methylation between DCIS and invasive breast cancer were confirmed. Further, the functional relevance of these genes’ regulation by methylation was demonstrated in early stage breast cancers from The Cancer Genome Atlas database.ConclusionsThis work contributes to the understanding of epigenetic alterations that occur in DCIS and illustrates the potential of DNA methylation as markers of DCIS progression.