Glioblastoma multiforme (GBM) is the most aggressive of all brain tumors, with a median survival of less than 1.5 years. Recently, epigenetic alterations were found to play key roles in both glioma genesis and clinical outcome, demonstrating the need to integrate genetic and epigenetic data in predictive models. To enhance current models through discovery of novel predictive biomarkers, we employed a genome-wide, agnostic strategy to specifically capture both methylation-directed changes in gene expression and alternative associations of DNA methylation with disease survival in glioma. Human GBM-associated DNA methylation, gene expression, IDH1 mutation status, and survival data were obtained from The Cancer Genome Atlas. DNA methylation loci and expression probes were paired by gene, and their subsequent association with survival was determined by applying an accelerated failure time model to previously published alternative and expression-based association equations. Significant associations were seen in 27 unique methylation/expression pairs with expression-based, alternative, and combinatorial associations observed (10, 13, and 4 pairs, respectively). The majority of the predictive DNA methylation loci were located within CpG islands, and all but three of the locus pairs were negatively correlated with survival. This finding suggests that for most loci, methylation/expression pairs are inversely related, consistent with methylation-associated gene regulatory action. Our results indicate that changes in DNA methylation are associated with altered survival outcome through both coordinated changes in gene expression and alternative mechanisms. Furthermore, our approach offers an alternative method of biomarker discovery using a priori gene pairing and precise targeting to identify novel sites for locus-specific therapeutic intervention.

Aim

We examined whether variation in blood-based epigenome-wide association studies could be more completely explained by augmenting existing reference DNA methylation libraries.

Materials & methods

We compared existing and enhanced libraries in predicting variability in three publicly available 450K methylation datasets that collected whole-blood samples. Models were fit separately to each CpG site and used to estimate the additional variability when adjustments for cell composition were made with each library.

Results

Calculation of the mean difference in the CpG-specific residual sums of squares error between models for an arthritis, aging and metabolic syndrome dataset, indicated that an enhanced library explained significantly more variation across all three datasets (p < 10(-3)).

Conclusion

Pathologically important immune cell subtypes can explain important variability in epigenome-wide association studies done in blood.

Natural killer (NK) cells are now recognized to exhibit characteristics akin to cells of the adaptive immune system. The generation of adaptive memory is linked to epigenetic reprogramming including alterations in DNA methylation. The study herein found reproducible genome wide DNA methylation changes associated with human NK cell activation. Activation led predominately to CpG hypomethylation (81% of significant loci). Bioinformatics analysis confirmed that non-coding and gene-associated differentially methylated sites (DMS) are enriched for immune related functions (i.e., immune cell activation). Known DNA methylation-regulated immune loci were also identified in activated NK cells (e.g., TNFA, LTA, IL13, CSF2). Twenty-one loci were designated high priority and further investigated as potential markers of NK activation. BHLHE40 was identified as a viable candidate for which a droplet digital PCR assay for demethylation was developed. The assay revealed high demethylation in activated NK cells and low demethylation in naïve NK, T- and B-cells. We conclude the NK cell methylome is plastic with potential for remodeling. The differentially methylated region signature of activated NKs revealed similarities with T cell activation, but also provided unique biomarker candidates of NK activation, which could be useful in epigenome-wide association studies to interrogate the role of NK subtypes in global methylation changes associated with exposures and/or disease states.

Quantitating the copy number of demethylated CpG promoter sites of the CD3Z gene can be used to estimate the numbers and proportions of T cells in human blood and tissue. Quantitative methylation specific PCR (qPCR) is useful for studying T cells but requires extensive calibration and is imprecise at low copy numbers. Here we compared the performance of a new digital PCR platform (droplet digital PCR or ddPCR) to qPCR using bisulfite converted DNA from 157 blood specimens obtained from ambulatory care controls and patients with primary glioma. We compared both ddPCR and qPCR with conventional flow cytometry (FACS) evaluation of CD3 positive T cells. Repeated measures on the same blood sample revealed ddPCR to be less variable than qPCR. Both qPCR and ddPCR correlated significantly with FACS evaluation of peripheral blood CD3 counts and CD3/total leukocyte values. However, statistical measures of agreement showed that linear concordance was stronger for ddPCR than for qPCR and the absolute values were closer to FACS for ddPCR. Both qPCR and ddPCR could distinguish clinically significant differences in T cell proportions and performed similarly to FACS. Given the higher precision, greater accuracy, and technical simplicity of ddPCR, this approach appears to be a superior DNA methylation based method than conventional qPCR for the assessment of T cells.

Background

Living near major roadways has been linked with increased risk of cardiovascular events and worse prognosis. Residential proximity to major roadways may also be associated with increased risk of hypertension, but few studies have evaluated this hypothesis.

Methods and results

We examined the cross-sectional association between residential proximity to major roadways and prevalent hypertension among 5401 postmenopausal women enrolled into the San Diego cohort of the Women's Health Initiative. We used modified Poisson regression with robust error variance to estimate the association between prevalence of hypertension and residential distance to nearest major roadway, adjusting for participant demographics, medical history, indicators of individual and neighborhood socioeconomic status, and for local supermarket/grocery and fast food/convenience store density. The adjusted prevalence ratios for hypertension were 1.22 (95% CI: 1.07, 1.39), 1.13 (1.00, 1.27), and 1.05 (0.99, 1.12) for women living ≤100, >100 to 200, and >200 to 1000 versus >1000 m from a major roadway (P for trend=0.006). In a model treating the natural log of distance to major roadway as a continuous variable, a shift in distance from 1000 to 100 m from a major roadway was associated with a 9% (3%, 16%) higher prevalence of hypertension.

Conclusions

In this cohort of postmenopausal women, residential proximity to major roadways was positively associated with the prevalence of hypertension. If causal, these results suggest that living close to major roadways may be an important novel risk factor for hypertension.

Living near major roadways has been associated with increased risk of cardiovascular morbidity and mortality, presumably from exposure to elevated levels of traffic-related air and/or noise pollution. This association may potentially be mediated through increased risk of incident hypertension, but results from prior studies are equivocal. Using Cox proportional hazards models we examined residential proximity to major roadways and incident hypertension among 38,360 participants of the Women's Health Initiative (WHI) Clinical Trial cohorts free of hypertension at enrollment and followed for a median of 7.9 years. Adjusting for participant demographics and lifestyle, trial participation, and markers of individual and neighborhood socioeconomic status, the hazard ratios for incident hypertension were 1.13 (95% CI: 1.00, 1.28), 1.03 (0.95, 1.11), 1.05 (0.99, 1.11), and 1.05 (1.00, 1.10) for participants living ≤50, >50-200, >200-400, and >400-1000 m vs >1000 m from the nearest major roadway, respectively (ptrend=0.013). This association varied substantially by WHI study region with hazard ratios for women living ≤50 m from a major roadway of 1.61 (1.18, 2.20) in the West, 1.51 (1.22, 1.87) in the Northeast, 0.89 (0.70, 1.14) in the South, and 0.94 (0.75, 1.19) in the Midwest. In this large, national cohort of post-menopausal women, residential proximity to major roadways was associated with incident hypertension in selected regions of the U.S. If causal, these results suggest residential proximity to major roadways, as a marker for air, noise and other traffic-related pollution, may be a risk factor for hypertension.

Linear regression is often used to estimate associations between chemical exposures and neurodevelopment at the mean of the outcome. However, the potential effect of chemicals may be greater among individuals at the 'tails' of outcome distributions. Here, we investigated distributional effects on the associations between gestational phthalate exposure and child Autism Spectrum Disorder (ASD)-related behaviors using quantile regression. We harmonized data from the Early Autism Risk Longitudinal Investigation (EARLI) (n = 140) Study, an enriched-risk cohort of mothers who had a child with ASD, and the Health Outcomes and Measures of the Environment (HOME) Study (n = 276), a general population cohort. We measured concentrations of 9 phthalate metabolites in urine samples collected twice during pregnancy. Caregivers reported children's ASD-related behaviors using the Social Responsiveness Scale (SRS) at age 3-8 years; higher scores indicate more ASD-related behaviors. In EARLI, associations between phthalate concentrations and SRS scores were predominately inverse or null across SRS score quantiles. In HOME, positive associations of mono-n-butyl phthalate, monobenzyl phthalate, mono-isobutyl phthalate, and di-2-ethylhexyl phthalate concentrations with SRS scores increased in strength from the median to 95th percentile of SRS scores. These results suggest associations between phthalate concentrations and SRS scores may be stronger in individuals with higher SRS scores.

Background

DNA methylation alterations may underlie associations between gestational perfluoroalkyl substances (PFAS) exposure and later-life health outcomes. To the best of our knowledge, no longitudinal studies have examined the associations between gestational PFAS and DNA methylation.

Objectives

We examined associations of gestational PFAS exposure with longitudinal DNA methylation measures at birth and in adolescence using the Health Outcomes and Measures of the Environment (HOME) Study (2003-2006; Cincinnati, Ohio).

Methods

We quantified serum concentrations of perfluorooctanoate (PFOA), perfluorooctane sulfonate (PFOS), perfluorononanoate (PFNA), and perfluorohexane sulfonate (PFHxS) in mothers during pregnancy. We measured DNA methylation in cord blood (n=266) and peripheral leukocytes at 12 years of age (n=160) using the Illumina HumanMethylation EPIC BeadChip. We analyzed associations between log2-transformed PFAS concentrations and repeated DNA methylation measures using linear regression with generalized estimating equations. We included interaction terms between children's age and gestational PFAS. We performed Gene Ontology enrichment analysis to identify molecular pathways. We used Project Viva (1999-2002; Boston, Massachusetts) to replicate significant associations.

Results

After adjusting for covariates, 435 cytosine-guanine dinucleotide (CpG) sites were associated with PFAS (false discovery rate, q<0.05). Specifically, we identified 2 CpGs for PFOS, 12 for PFOA, 8 for PFHxS, and 413 for PFNA; none overlapped. Among these, 2 CpGs for PFOA and 4 for PFNA were replicated in Project Viva. Some of the PFAS-associated CpG sites annotated to gene regions related to cancers, cognitive health, cardiovascular disease, and kidney function. We found little evidence that the associations between PFAS and DNA methylation differed by children's age.

Discussion

In these longitudinal data, PFAS biomarkers were associated with differences in several CpGs at birth and at 12 years of age in or near genes linked to some PFAS-associated health outcomes. Future studies should examine whether DNA methylation mediates associations between gestational PFAS exposure and health. https://doi.org/10.1289/EHP10118.

Maternal nutrition during gestation has been investigated for its role in child neurodevelopment. However, little is known about the potential impact of gestational caffeine exposure on child autistic behaviors. Here, we assess the relation between maternal caffeine intake during pregnancy and children's behavioral traits related to Autism Spectrum Disorder (ASD). We harmonized data from two pregnancy cohorts, Early Autism Risk Longitudinal Investigation (EARLI) (n = 120), an enriched-risk cohort of mothers who previously had a child with ASD, from Pennsylvania, Maryland, and Northern California (2009-2012), and the Health Outcomes and Measures of the Environment (HOME) Study (n = 269), a general population cohort from Cincinnati, Ohio (2003-2006). Mothers self-reported caffeine intake twice during pregnancy. Caregivers reported child behavioral traits related to ASD using the Social Responsiveness Scale (SRS) when children were aged 3-8 years. Higher scores indicate more ASD-related behaviors. We estimated covariate-adjusted differences in continuous SRS T-scores per interquartile range increase in caffeine intake. Self-reported caffeine intake during pregnancy was positively associated with SRS T-scores among children in EARLI (β: 2.0; 95% CI -0.1, 4.0), but to a lesser extent in HOME (β: 0.6; 95% CI -0.5, 1.6). In HOME, pre-pregnancy body mass index (BMI) modified the association between caffeine intake and SRS T-scores, where more positive associations were observed among women with higher BMIs. Our findings suggest gestational caffeine intake may represent a marker of vulnerability to childhood ASD-related behaviors. Additional studies are warranted to extend these findings.

Background

Myeloid-derived suppressor cells (MDSCs), which include monocytic (mMDSCs) and granulocytic (gMDSCs) cells, are an immunosuppressive, heterogeneous population of cells upregulated in cancer and other pathologic conditions, in addition to normal conditions of stress. The origin of MDSCs is debated, and the regulatory pattern responsible for gMDSC differentiation remains unknown. Since DNA methylation (DNAm) contributes to lineage differentiation, we have investigated whether it contributes to the acquisition of the gMDSC phenotype.

Results

Using the Illumina EPIC array to measure DNAm of gMDSCs and neutrophils from diverse neonatal and adult blood sources, we found 189 differentially methylated CpGs between gMDSCs and neutrophils with a core of ten differentially methylated CpGs that were consistent across both sources of cells. Genes associated with these loci that are involved in immune responses include VCL, FATS, YAP1, KREMEN2, UBTF, MCC-1, and EFCC1. In two cancer patient groups that reflected those used to develop the methylation markers (head and neck squamous cell carcinoma (HNSCC) and glioma), all of the CpG loci were differentially methylated, reaching statistical significance in glioma cases and controls, while one was significantly different in the smaller HNSCC group.

Conclusions

Our findings indicate that gMDSCs have a core of distinct DNAm alterations, informing future research on gMDSC differentiation and function.