Association of Glyphosate Exposure with Blood DNA Methylation in a Cross-Sectional Study of Postmenopausal Women.

BACKGROUND: Glyphosate is the most commonly used herbicide in the world and is purported to have a variety of health effects, including endocrine disruption and an elevated risk of several types of cancer. Blood DNA methylation has been shown to be associated with many other environmental exposures, but to our knowledge, no studies to date have examined the association between blood DNA methylation and glyphosate exposure.
OBJECTIVE: We conducted an epigenome-wide association study to identify DNA methylation loci associated with urinary glyphosate and its metabolite aminomethylphosphonic acid (AMPA) levels. Secondary goals were to determine the association of epigenetic age acceleration with glyphosate and AMPA and develop blood DNA methylation indices to predict urinary glyphosate and AMPA levels.
METHODS: For 392 postmenopausal women, white blood cell DNA methylation was measured using the Illumina Infinium MethylationEPIC BeadChip array. Glyphosate and AMPA were measured in two urine samples per participant using liquid chromatography-tandem mass spectrometry. Methylation differences at the probe and regional level associated with glyphosate and AMPA levels were assessed using a resampling-based approach. Probes and regions that had an false discovery rate q<0.1 in ≥90% of 1,000 subsamples of the study population were considered differentially methylated. Differentially methylated sites from the probe-specific analysis were combined into a methylation index. Epigenetic age acceleration from three epigenetic clocks and an epigenetic measure of pace of aging were examined for associations with glyphosate and AMPA.
RESULTS: We identified 24 CpG sites whose methylation level was associated with urinary glyphosate concentration and two associated with AMPA. Four regions, within the promoters of the MSH4, KCNA6, ABAT, and NDUFAF2/ERCC8 genes, were associated with glyphosate levels, along with an association between ESR1 promoter hypomethylation and AMPA. The methylation index accurately predicted glyphosate levels in an internal validation cohort. AMPA, but not glyphosate, was associated with greater epigenetic age acceleration. DISCUSSION: Glyphosate and AMPA exposure were associated with DNA methylation differences that could promote the development of cancer and other diseases. Further studies are warranted to replicate our results, determine the functional impact of glyphosate- and AMPA-associated differential DNA methylation, and further explore whether DNA methylation could serve as a biomarker of glyphosate exposure. https://doi.org/10.1289/EHP10174.

Introduction

Glyphosate is the most used pesticide in the world,[1] used in both agricultural (for weed control and preharvest desiccation) and nonagricultural settings.[1,2,3] Studies have detected glyphosate in the air, soil, drinking water, and food.[4] Use of glyphosate-based herbicides has increased dramatically since their introduction,[1] largely due to the growing use of genetically modified glyphosate-resistant crops starting in the late 1990s.[5] Glyphosate and/or its primary metabolite, aminomethylphosphonic acid (AMPA), are frequently detected in the food supply at different levels,[6,7,8,9,10] and recent studies in the United States have detected glyphosate in the urine of 70%–90% of participants.[11,12,13,14]

Concerns have been raised about the safety of glyphosate exposure in humans, and the topic remains controversial. In 2015 the International Agency for Research on Cancer (IARC) classified glyphosate as a probable human carcinogen.[15] This classification was supported by a meta-analysis that suggested an elevated risk of non-Hodgkin’s lymphoma associated with glyphosate exposure,[16] although another large cohort study did not find a relationship between glyphosate and cancer risk.[17] Epidemiological studies have also found associations with other health problems, including shortened gestational length,[14] birth defects,[18] endocrine disruption,[12] and thyroid dysfunction.[19] Animal and in vitro studies have shown associations with endocrine disruption,[20,21] fatty liver disease,[22] breast cancer cell proliferation,[23] changes to the microbiome,[24,25,26,27] and increased risk of antibiotic resistance in bacteria.[28] Although limited data exists about potential adverse effects of AMPA, some studies suggest that AMPA may have effects similar to those of glyphosate.[29,30,31] One recent study found elevated breast cancer risk among women with high urinary AMPA.[32]

Epigenetic markers, such as DNA methylation, may be a powerful tool for understanding the potential effects of glyphosate exposure in humans. White blood cell (WBC) DNA methylation has been associated with environmental exposures,[33,34] including endocrine-disrupting compounds,[35,36] other pesticides,[37] and air pollution.[38,39,40] In addition, DNA methylation indices combining information from multiple sites have been developed as biomarkers of some exposures, including smoking[41,42] and alcohol consumption.[43] Some smoking-associated DNA methylation markers have also been shown to be associated with risk for lung[44] and bladder[45] cancers.

Recent in vitro[29,46] and animal[47,48] evidence (reviewed in reference[49]) suggests that glyphosate and AMPA may influence DNA methylation, including in ways that could promote the development of breast cancer.[50] However, the association of glyphosate and AMPA with DNA methylation differences in humans has been rarely studied. One recent study of pesticide applicators identified a single CpG site (cg06950346) at which methylation was associated with history of self-reported glyphosate use.[51] To our knowledge, no studies to date have examined direct glyphosate and AMPA measurement in association with blood DNA methylation.

Another recent avenue of research is the development of “epigenetic clocks” that purport to measure biological age using DNA methylation at a selection of CpG sites.[52,53,54] Epigenetic age acceleration (biological age greater than chronological age) or a faster molecular pace of aging[55] have been associated with a host of adverse outcomes, including all-cause mortality,[56,57,58] risk for breast cancer,[59] B-cell lymphoma,[60] which have been associated with glyphosate[16] and AMPA[32] in epidemiological studies. Epigenetic age acceleration has also been associated with other environmental exposures,[61,62,63,64] but to our knowledge, the association between epigenetic age acceleration and glyphosate and AMPA measurements has not been previously explored.

We conducted an epigenome-wide association study to identify DNA methylation loci associated with urinary glyphosate and AMPA levels. Secondary goals were to determine the association of epigenetic age acceleration with glyphosate and AMPA and develop blood DNA methylation indices to predict urinary glyphosate and AMPA levels.

Methods

Study Population

The study population consisted of postmenopausal women residing in southern California between the ages of 45 and 66 y (). All participants provided informed consent to participate, and the study was approved by the University of California, Irvine institutional review board, HS #2016-3,127. Recruitment, questionnaires, and specimen collection were described in detail previously.[65] Briefly, women were recruited from 2017 to 2019 via two recruitment arms: targeted mailings to members of the Athena Breast Health Network (a cohort of women receiving screening mammograms at University of California, Irvine health facilities), and untargeted recruitment from the broader community via flyers, social media, events, and word of mouth. The study was designed to identify epigenetic markers for factors potentially related to breast cancer risk; thus, women with a personal history of breast cancer or mastectomy were excluded from the study.

Biospecimen Collection and Processing

Specimens were collected and processed as previously described.[65] Briefly, participants self-collected first-morning urine samples on 2 d within a 10-d period not known to them in advance to capture their typical behaviors and exposures. Urine samples were stored in a freezer () until they could be transported to the laboratory. At an appointment attempted to be scheduled within 10 d of the second urine sample collection (median 1 d, range 0–24 d), a certified phlebotomist collected peripheral blood in three glass Vacutainer blood collection tubes: one containing ethylenediaminetetraacetic acid (EDTA), one containing acid citrate dextrose (ACD), and one with no anticoagulants. DNA was extracted from the buffy coat (from EDTA and ACD tubes, pooled) within 6 h of collection using the QIAamp DNA Blood Maxi Kit (Cat. No. 51,194, QIAGEN) and stored at for later analysis.

Urinary glyphosate and AMPA were measured using liquid chromatography–tandem mass spectrometry (LC-MS/MS) at the Collaborative Center for Translational Mass Spectrometry (CCTMS) (Tgen) using a Vanquish UHPLC coupled to TSQ Altis triple quadrupole mass spectrometer (Thermo Scientific) as previously described.[12] Briefly, glyphosate and AMPA assay validation was performed in a commercially available urine pool from LeeBio Solutions prior to analysis. The calibration curves for glyphosate and AMPA were prepared over a linear range of (coefficient of determination ) by spiking variable concentrations of glyphosate and AMPA and their respective isotopically labeled internal standards, and (Sigma-Aldrich) at a fixed concentration of . Data acquisition and processing were performed using Xcalibur 4.1.50 and QuanBrowser, 4.1.50 (Thermo Scientific). Both assays were linear (coefficient of determination ) over a range (Figure S1). The limits of detection (LOD) for glyphosate and AMPA were 0.014 and , and the limits of quantitation (LOQ) were 0.041 and , respectively (Table S1). Creatinine was measured using the DetectX urinary creatinine detection kit (Arbor Assays, K002-H5) according to the manufacturer’s instructions. The data were tested for batch effects using the Kruskal-Wallis test, and batch correction was performed for glyphosate values using the removeBatchEffects function in limma.[66] Four samples with implausibly low creatinine () were excluded from subsequent analysis. Measurements were averaged for the two urine samples for each analyte after replacing values with .[67]

Genomic DNA was bisulfite converted using the Zymo EZ DNA methylation kit (Zymo Research), and then DNA methylation at over 850,000 CpG sites was measured using the Illumina Infinium MethylationEPIC BeadChip (Illumina) at the University of Southern California Molecular Genomics Core. Laboratory staff were blinded to glyphosate and AMPA measurements. Methylation data preprocessing was performed according to recommended steps,[68] including probe filtering, normalization, and batch correction. All data analysis was performed in R, (version 3.6.2; R Development Core Team). Probes with a detection were considered missing, and other low-quality probes were removed as follows: missing in at least 20% of samples (); had SNPs with global minor allele frequency within 5 base pairs of the target sequence or mapping problems with the probe sequence ();[69] hybridized to multiple locations ();[70] or located on the X or Y chromosome (). All samples passed quality control, including efficiency of bisulfite conversion, verification of reported sex and identity, and having failed probes (maximum value: 0.13%), and were included in the final data analysis. Values were normalized with noob normalization as implemented in the minfi package (version 1.32.0)[71] for background correction and dye bias adjustment followed by beta mixture quantile normalization (BMIQ) to correct for type II probe bias.[72] Except where noted, filtering and normalization was completed using the ChAMP package (version 2.12.4).[73] Finally, DNA methylation measurements were corrected for batch and position on chip using ComBat[74] as implemented in sva version 3.30.1.[75] We used a reference-based method[76] to estimate the proportions of six WBC types in our samples, which did not significantly differ with glyphosate or AMPA tertile (Table S2). We used the methylation M value (the base 2 logarithm of the ratio of methylated to unmethylated intensities) for all analyses; in some cases, we also reported values (percent methylation) for ease of interpretation.[77]

Covariates

Relevant covariates, including age, self-reported race/ethnicity, height, weight, smoking status (current, former, or never smoker), alcohol consumption, organic eating habits, physical activity, and recent herbicide use, were collected via questionnaires. Race and ethnicity were self-reported in response to two questions: “What is your racial background?” (response choices: “Black or African American,” “White,” “Asian,” “American Indian or Alaska Native,” “Native Hawaiian or other Pacific Islander,” “Some other race,” “Do not know,” or “Prefer not to answer”) and “Are you of Hispanic, Latino, or Spanish origin or ancestry?” (yes/no; “yes” responses were prompted to further self-identify as “Mexican, Mexican American, or Chicano,” “Puerto Rican,” “Cuban,” and/or “Other Hispanic, Latino, or Spanish origin”). Responses to these questions were pooled into a combined race/ethnicity measure with categories non-Hispanic White, non-Hispanic Asian, non-Hispanic Black, Hispanic, or Other. All those who answered “yes” to the Hispanic/Latino/Spanish origin question were placed into the Hispanic category, regardless of response to the race question. Multiple selections were allowed; individuals who selected more than one race () were placed into the “Other” category. Those who answered “American Indian or Alaska Native,” “Native Hawaiian or other Pacific Islander,” or “Some other race” were placed into the “Other” category due to small numbers in the study cohort (total including ). Race/ethnicity was included in the analysis due to known associations with DNA methylation.[78,79,80]

Self-reported height and weight were used to calculate body mass index (BMI) by dividing the weight in kilograms by the height in meters squared. Organic eating habits were self-reported as “seldom or never,” “sometimes,” or “often or always.”[81] Physical activity was reported as typical frequency and duration of moderate and vigorous physical activity. These responses were converted to minutes per week of moderate exercise, with each minute of vigorous exercise equivalent to 2 min of moderate exercise.[82] Diet quality was estimated by using up to three unannounced ASA24 24-h dietary recalls,[83] completed during the 10 d prior to blood collection, to calculate the Healthy Eating Index-2015 (HEI) averaged across all recalls. The HEI is scored on a 0–100 scale, with higher values indicating greater adherence to the Dietary Guidelines for Americans.[84] Some participants were asked whether they had used herbicides at home or work within the past 7 d. This question was added to the questionnaire partway through study recruitment, so data on recent herbicide use is missing for 37% () of participants.

Statistical Analyses

Variability and correlation of urinary glyphosate and AMPA levels between and within samples were characterized with the intraclass correlation (ICC) and Pearson’s correlation, respectively. The association between all covariates and the natural logarithm of glyphosate and AMPA concentration (adjusted for urinary creatinine) was assessed using linear regression. In all linear regression analyses, two-sided were considered statistically significant. In the DNA methylation analyses, missing data for covariates (from individuals) were substituted with the median (numeric variables) or mode (categorical variables).

DNA Methylation Analysis

We calculated epigenetic age for each sample according to three epigenetic clocks: Hannum,[52] Horvath,[53] and Levine,[54] and one epigenetic measure of pace of aging, DunedinPoAm.[55] The Hannum and Horvath clocks are designed to correlate with chronological age, whereas the Levine clock also incorporates phenotypic measures, including mortality and physical functioning. DunedinPoAm is intended to capture the current rate of change of biological age. Epigenetic age was calculated from the model coefficients provided by each author using normalized methylation beta values. Values for missing probes ( for Hannum, for Horvath, for Levine, for DunedinPoAm) were imputed using k-nearest neighbors imputation. DunedinPoAm was calculated using code provided by its developers.[55] Epigenetic age acceleration (residuals from a model regressing chronological age on epigenetic age from each clock) and epigenetic pace of aging from DunedinPoAm was examined for association with glyphosate and AMPA, first in univariate linear regression models adjusted for urinary creatinine concentration and then adjusted for age, race/ethnicity, BMI, smoking status, alcohol consumption, organic eating, and HEI, in addition to batch and position on chip. A third model also included estimated WBC type proportions predicted according to Houseman’s method.[76] These covariates were selected based on known relationships with either DNA methylation, glyphosate/AMPA, or both; creatinine was included to account for differences in urine concentration.

For probe- and region-specific analyses, the study samples were divided into training () and validation () sets stratified by glyphosate tertile. Demographic and dietary variables were compared for the training and validation sets (Table S3). All variables except smoking were not significantly associated with the randomly assigned set; there was a larger proportion of former smokers in the validation set (). Because all the former smokers in the validation set reported cessation prior to blood sample collection, we proceeded with the planned analysis. Using the training set, we identified differentially methylated probes (DMPs) and differentially methylated regions (DMRs) according to a resampling-based method.[85]

We selected a random subsample consisting of 90% of the training set (). Using this subsample, we identified candidate DMPs using limma (version 3.44.3),[66] to fit linear models and candidate DMRs using DMRcate (version 2.0.7),[86] with a kernel bandwidth of 1,000 base pairs and scaling factor of 2. Each model was adjusted for the same variables as in the final epigenetic aging model. Models without adjustment for dietary variables (organic eating and HEI) did not show substantially different results (Pearson’s for test statistics from the two versions of the model for both analytes), nor did models that used creatinine-standardized glyphosate and AMPA concentrations instead of including creatinine as a covariate in the model (Pearson’s for glyphosate and 0.985 for AMPA). A false discovery rate (FDR) q-value of was considered statistically significant. This process was repeated in a series of 1,000 subsamples of 90% of the training set and sampling distributions for all summary statistics constructed. Probes and regions that were selected as candidate DMPs or DMRs in of subsamples were considered differentially methylated and carried forward for further analysis. This approach results in a more stable list of DMPs and DMRs, because DNA methylation microarray results are known to be sensitive to small differences in the study cohort.[87,88,89,90] The median FDR q-value was used to rank the relative significance of DMPs. For DMRs, overlapping regions were combined, and then significance was ranked after combining median q-values for the probes within the region using Stouffer’s method. A traditional epigenome-wide association analysis with the entire training set[74] was also done to determine whether findings were robust to changes in the analysis approach.

DMPs and DMRs were annotated using the Illumina manifest (version B4) to identify associated genes and genomic context. Probes were also mapped to ChromHMM data from ENCODE[91] to determine predicted chromatin state; data from the GM12878 lymphoblastic cell line was used because it is the tissue most similar to leukocytes. Chromatin states were pooled into six categories as follows: promoter (active, weak, or poised promoters), enhancer (strong or weak enhancers), transcribed (transcriptional transition, transcriptional elongation, or weak transcribed), polycomb-repressed, inactive (heterochromatin and repetitive regions), and insulators. DMPs were examined for enrichment of certain locations relative to CpG islands or chromatin states, in comparison with the background of all probes included on the array, using Fisher’s exact test. We used the “gometh” function from the missMethyl package (version 1.16.0)[92] to examine enrichment of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and gene ontology terms in genes associated with DMPs and assessed enrichment of pathways using Ingenuity Pathway Analysis (IPA) software (version 62089861; release Date: 17 February 2021). Terms or pathways with and at least three genes associated with DMPs were considered enriched.

We developed methylation indices to predict urinary glyphosate and AMPA levels using blood DNA methylation at the glyphosate- and AMPA-associated DMPs. Methylation () values for each DMP were passed to an elastic net () regression model with 10-fold cross-validation repeated 100 times to estimate lambda (the penalty coefficient) and regression coefficients for each site. The final lambda was selected as the largest lambda within one standard error of the minimum prediction error. This approach was implemented using the glmnet package (version 2.0-18).[93] The performance of the methylation indices were assessed by calculating the Pearson correlation between predicted and actual glyphosate and AMPA concentration in the 60 validation samples. The indices’ relationship with glyphosate and AMPA tertile was also tested using analysis of variance (ANOVA), and the area under the receiver operating characteristic curve (AUC) was computed to characterize the discriminatory accuracy for the highest vs. lowest tertile.

Results

Glyphosate and AMPA Concentrations

A total of 96% of study participants had detectable () glyphosate, and 82% had detectable AMPA in at least one urine sample. Table 1 describes glyphosate and AMPA concentrations stratified by various demographic, dietary, and lifestyle factors. The median glyphosate concentration (averaged between two urine samples for each participant) was (range ), and the median AMPA concentration was (range ). Both glyphosate and AMPA concentrations were strongly right-skewed and thus transformed with the natural logarithm for all further analyses (Figure S2). Between-sample agreement was moderate for both glyphosate () and AMPA (), as was the within-sample correlation between glyphosate and AMPA measurements (Pearson’s , ) (Figure S3).

Table 1

Median urinary glyphosate and AMPA concentrations, averaged for two samples collected between 2017 and 2019 from 392 postmenopausal California women, stratified by cohort characteristics.

	n (%)	Glyphosate			AMPA
		ng/mLmedian (IQR)	μg/g creatininemedian (IQR)	p-Value	ng/mLmedian (IQR)	μg/g creatininemedian (IQR)	p-Value
Overall	392	0.12 (0.06, 0.22)	0.20 (0.11, 0.38)	NA	0.06 (0.02, 0.12)	0.10 (0.04, 0.22)	NA
Race/ethnicity
Asian	43 (11.1)	0.11 (0.07, 0.23)	0.19 (0.09, 0.34)	0.30	0.08 (0.02, 0.14)	0.13 (0.03, 0.24)	0.50
Hispanic	69 (17.9)	0.12 (0.07, 0.21)	0.20 (0.11, 0.37)	0.90	0.06 (0.02, 0.11)	0.08 (0.04, 0.17)	0.41
Other	20 (5.2)	0.07 (0.04, 0.13)	0.11 (0.06, 0.24)	0.004	0.02 (0.01, 0.05)	0.04 (0.02, 0.07)	0.005
White	254 (65.8)	0.12 (0.06, 0.24)	0.21 (0.12, 0.42)	Ref	0.05 (0.02, 0.12)	0.10 (0.04, 0.22)	Ref
Missing	6	—	—	—	—	—	—
Age (y)
45–49	27 (6.9)	0.18 (0.09, 0.24)	0.23 (0.11, 0.33)	0.72	0.11 (0.05, 0.19)	0.12 (0.05, 0.29)	0.12
50–54	85 (21.7)	0.14 (0.06, 0.19)	0.18 (0.12, 0.33)		0.07 (0.02, 0.11)	0.09 (0.04, 0.17)
55–59	155 (39.5)	0.11 (0.06, 0.21)	0.18 (0.10, 0.33)		0.05 (0.02, 0.12)	0.08 (0.03, 0.22)
60–66	125 (31.9)	0.11 (0.06, 0.29)	0.23 (0.11, 0.56)		0.05 (0.01, 0.12)	0.10 (0.04, 0.22)
Missing	0	—	—	—	—	—	—
BMI (kg/m^2)
<25	188 (48.0)	0.12 (0.06, 0.24)	0.21 (0.10, 0.44)	0.78	0.06 (0.02, 0.12)	0.10 (0.04, 0.21)	0.63
25–29.9	116 (29.6)	0.12 (0.06, 0.22)	0.20 (0.12, 0.33)		0.05 (0.02, 0.12)	0.09 (0.04, 0.20)
≥30	88 (22.4)	0.13 (0.07, 0.19)	0.19 (0.12, 0.31)		0.06 (0.01, 0.13)	0.09 (0.04, 0.22)
Missing	0	—	—	—	—	—	—
Smoking status
Never	285 (72.9)	0.12 (0.06, 0.22)	0.22 (0.11, 0.38)	Ref	0.06 (0.02, 0.13)	0.11 (0.04, 0.22)	Ref
Former	89 (22.8)	0.10 (0.05, 0.22)	0.17 (0.10, 0.38)	0.52	0.04 (0.01, 0.09)	0.07 (0.03, 0.15)	0.04
Current	17 (4.3)	0.11 (0.05, 0.21)	0.17 (0.09, 0.33)	0.60	0.05 (0.02, 0.14)	0.09 (0.03, 0.26)	0.59
Missing	1	—	—	—	—	—	—
Alcohol (drinks/wk)
None	107 (27.4)	0.12 (0.06, 0.22)	0.22 (0.11, 0.48)	Ref	0.05 (0.02, 0.11)	0.08 (0.04, 0.22)	Ref
1 or fewer	161 (41.2)	0.12 (0.05, 0.24)	0.19 (0.09, 0.37)	0.14	0.06 (0.01, 0.12)	0.08 (0.03, 0.18)	0.49
2–6	72 (18.4)	0.14 (0.07, 0.19)	0.23 (0.14, 0.41)	0.92	0.06 (0.03, 0.12)	0.13 (0.04, 0.24)	0.64
7 or more	51 (13.0)	0.10 (0.05, 0.21)	0.16 (0.12, 0.25)	0.17	0.07 (0.03, 0.15)	0.12 (0.06, 0.21)	0.16
Missing	1	—	—	—	—	—	—
Organic eating
Seldom/Never	124 (31.7)	0.14 (0.07, 0.23)	0.21 (0.13, 0.38)	Ref	0.07 (0.03, 0.13)	0.12 (0.05, 0.23)	Ref
Sometimes	114 (29.2)	0.12 (0.05, 0.24)	0.19 (0.11, 0.38)	0.47	0.06 (0.02, 0.13)	0.10 (0.04, 0.25)	0.44
Often/Always	153 (39.1)	0.10 (0.06, 0.20)	0.19 (0.09, 0.38)	0.05	0.04 (0.01, 0.11)	0.07 (0.03, 0.18)	0.01
Missing	1	—	—	—	—	—	—
HEI
Quartile 1	93 (25.1)	0.13 (0.07, 0.22)	0.22 (0.13, 0.37)	0.99	0.07 (0.02, 0.12)	0.12 (0.05, 0.22)	0.03
Quartile 2	93 (25.1)	0.10 (0.05, 0.19)	0.19 (0.09, 0.35)		0.07 (0.02, 0.13)	0.13 (0.04, 0.26)
Quartile 3	92 (24.8)	0.11 (0.06, 0.22)	0.17 (0.10, 0.41)		0.05 (0.02, 0.12)	0.08 (0.04, 0.16)
Quartile 4	93 (25.1)	0.12 (0.06, 0.25)	0.24 (0.11, 0.54)		0.04 (<LOD, 0.12)	0.06 (0.03, 0.21)
Missing	21	—	—	—	—	—	—
Physical activity
≥150 min/wk	153 (40.7)	0.10 (0.05, 0.24)	0.19 (0.09, 0.38)	0.03	0.05 (0.01, 0.11)	0.08 (0.04, 0.19)	0.08
<150 min/wk	223 (59.3)	0.13 (0.07, 0.22)	0.22 (0.12, 0.41)		0.06 (0.02, 0.13)	0.11 (0.04, 0.23)
Missing	16	—	—	—	—	—	—
Herbicide use (past week)
Yes	21 (8.5)	0.14 (0.08, 0.20)	0.22 (0.11, 0.44)	0.61	0.04 (<LOD, 0.11)	0.11 (0.04, 0.16)	0.39
No	227 (91.5)	0.12 (0.06, 0.24)	0.21 (0.11, 0.38)		0.06 (0.02, 0.12)	0.10 (0.03, 0.20)
Missing	144	—	—	—	—	—	—

Note: -Values are from linear regression with the natural logarithm of glyphosate or AMPA as outcome, adjusted for urinary creatinine. LOD: for glyphosate and for AMPA. Values of glyphosate or AMPA were substituted with . Age, BMI, and HEI were evaluated as continuous variables; results were comparable when evaluated as categorical variables. —, no data; AMPA, aminomethylphosphonic acid; BMI, body mass index; HEI, Healthy Eating Index; LOD, limit of detection; min, minutes; NA, not applicable; Ref, reference.

Women who reported “often or always” eating organic food had lower glyphosate [median (interquartile range, IQR) 0.10 (0.06, 0.20) ng/mL compared with 0.14 (0.07, 0.23) ng/mL for “seldom or never,” ] and AMPA [median (IQR) 0.04 (0.01, 0.11) ng/mL compared with 0.07 (0.03, 0.13) ng/mL for “seldom or never,” ]. Lower diet quality was associated with elevated AMPA [median (IQR) 0.07 (0.02, 0.12) ng/mL for highest quartile of diet quality vs. 0.04 (, 0.12) ng/mL for lowest quartile, ], but not glyphosate (). Those who met the Physical Activity Guidelines for Americans ( of moderate-intensity exercise per week) had lower glyphosate [median (IQR) 0.10 (0.05, 0.24) ng/mL for those meeting guidelines vs. 0.13 (0.07, 0.22) for those not meeting, ]. Former smokers had marginally lower AMPA concentrations in comparison with never smokers [median (IQR) 0.04 (0.01, 0.09) ng/mL vs. 0.06 (0.02, 0.13) ng/mL, ]. Women in the “Other” race/ethnicity category had lower concentrations of both glyphosate [median (IQR) 0.07 (0.04, 0.13) vs. 0.12 (0.06, 0.24) ng/mL, ] and AMPA [median (IQR) 0.02 (0.01, 0.05) vs. 0.05 (0.02, 0.12) ng/mL, ] in comparison with White women, although there were only 20 women in the “Other” category. None of the other variables examined (age, BMI, alcohol consumption, or herbicide use) were associated with urinary glyphosate or AMPA.

Each of the three epigenetic clocks (Hannum, Horvath, and Levine) was well-correlated with chronologic age (Pearson’s , 0.64, and 0.54, respectively), although the Hannum and Levine clocks were poorly calibrated in this sample (Figure S4). Glyphosate was not significantly associated with differences in epigenetic age acceleration according to any of the epigenetic clocks (Table 2). AMPA was associated with increased epigenetic age acceleration according to the Hannum clock [adjusted coefficient 0.36 y; 95% confidence interval (CI): 0.03, 0.70, ]. AMPA was also associated with increased epigenetic age acceleration according to the Horvath clock in the adjusted model (adjusted coefficient 0.41; 95% CI: 0.02, 0.80, ) but not the univariate model (); the opposite was true for the Levine clock (univariate coefficient 0.61; 95% CI: 0.04, 1.18, ; adjusted ). Epigenetic pace of aging according to DunedinPoAm was not significantly associated with glyphosate or AMPA concentration (Table 2).

Table 2

Association of epigenetic age acceleration from three epigenetic clocks with urinary glyphosate and AMPA concentration in 392 postmenopausal California women.

	Glyphosate			AMPA
	Coefficient (95% CI)	p-Value	Adjusted R^2	Coefficient (95% CI)	p-Value	Adjusted R^2
Univariate
Hannum	0.08 (−0.27, 0.44)	0.64	0.0058	0.35 (0.03, 0.67)	0.03	0.017
Horvath	−0.17 (−0.55 , 0.21)	0.39	−0.0025	0.21 (−0.14 , 0.56)	0.24	−0.00076
Levine	0.34 (−0.30, 0.97)	0.30	−0.0023	0.61 (0.04, 1.18)	0.04	0.0061
DunedinPoAm	0.08 (−0.02, 0.17)	0.11	0.016	0.04 (−0.05, 0.13)	0.36	0.011
Adjusted, no WBC types
Hannum	−0.24 (−0.64, 0.16)	0.24	0.056	0.43 (0.07, 0.80)	0.02	0.056
Horvath	−0.36 (−0.79, 0.07)	0.10	0.037	0.40 (0.001, 0.79)	0.0499	0.037
Levine	0.02 (−0.69, 0.72)	0.97	0.047	0.49 (−0.17, 1.14)	0.14	0.047
DunedinPoAm	0.08 (−0.01, 0.18)	0.09	0.25	0.02 (−0.07, 0.11)	0.64	0.25
Adjusted, with WBC types
Hannum	−0.22 (−0.58, 0.15)	0.24	0.22	0.36 (0.03, 0.70)	0.04	0.22
Horvath	−0.33 (−0.75, 0.10)	0.13	0.075	0.41 (0.02, 0.80)	0.04	0.075
Levine	0.03 (−0.60, 0.66)	0.93	0.26	0.30 (−0.29, 0.88)	0.32	0.26
DunedinPoAm	0.07 (−0.02, 0.16)	0.11	0.39	0.01 (−0.07, 0.09)	0.85	0.39

Note: Coefficients are from linear regression with epigenetic age acceleration (years) from the Hannum, Horvath, and Levine epigenetic clocks or DunedinPoAm (z-score) as the dependent variable and the natural logarithm of glyphosate or AMPA concentration as the independent variable, adjusted for urinary creatinine. The adjusted model is additionally adjusted for age, race/ethnicity, body mass index, smoking status, alcohol consumption, self-reported organic eating habits, diet quality (Healthy Eating Index), batch, and position on chip. The adjusted model with WBC types is additionally adjusted for WBC type proportions estimated via Houseman’s method. Glyphosate and AMPA concentrations are averaged for two urine samples collected within approximately 10 d of each other. Missing data for covariates (from individuals) were substituted with the median (numeric variables) or mode (categorical variables). AMPA, aminomethylphosphonic acid; CI, confidence interval; WBC, white blood cell.

Twenty-four probes associated with urinary glyphosate concentration and two with urinary AMPA were identified in of subsamples and considered differentially methylated (Table 3). Seventeen of the 24 probes (71%) were hypomethylated with higher glyphosate (Figure 1), whereas both AMPA-associated probes were hypermethylated (Table 3). Glyphosate-associated probes with the smallest median -values were located within the PEX26, SF3B2, and CHMP1A genes. The largest methylation difference between glyphosate tertiles were observed in probes within the VMO1, KCP, and QARS1 genes. Our findings were robust to changes in the analysis approach: All glyphosate- and AMPA-associated DMPs were statistically significant (FDR ) using a traditional epigenome-wide association analysis with the entire training set, except for one AMPA-associated DMP (cg23261070, ).

Table 3

Differentially methylated probes (DMPs) associated with urinary glyphosate and AMPA concentration in 332 postmenopausal California women.

Probe	Freq.	Median p-value	Median FDR q-value	Methylation (β) difference per unit glyphosate/AMPA	Chr.	Position (hg19)	Gene(s)	Gene summary	Location within gene(s)	CpG Island	Predicted chromatin state
Glyphosate
cg13499896	1,000	1.74×10^−8	0.006	−0.71	22	18,571,657	PEX26	Probably required for protein import into peroxisomes.	3’UTR	Open Sea	Strong enhancer
cg24576174	1,000	4.24×10^−8	0.010	−0.60	11	65,828,083	SF3B2	Involved in pre-mRNA splicing as a component of the splicing factor SF3B complex.	Body	Open Sea	Transcription elongation
cg00355690	996	5.64×10^−8	0.011	−1.89	16	89,722,417	CHMP1A	Involved in multivesicular body sorting of proteins to the interiors of lysosome. Involved in cytokinesis. May also be involved in chromosome condensation. May play a role in stable cell cycle progression and in PcG gene silencing.	Body	N. Shore	Weak enhancer
cg26833395	1,000	5.76×10^−8	0.012	−0.11	20	30,102,102	HM13	Plays a key role in generating signal sequence-derived human lymphocyte antigen-E epitopes that are recognized by the immune system.	Promoter	Island	Weak promoter
cg06993862	987	2.62×10^−7	0.023	−2.25	3	49,141,001	QARS1	Plays a critical role in brain development.	Body	N. Shore	Strong enhancer
cg13310154	978	4.56×10^−7	0.029	−0.11	17	55,910,843	—	—	Intergenic	Open Sea	Insulator
cg16601151	959	5.92×10^−7	0.031	−0.11	16	19,079,341	COQ7	Involved in ubiquinone biosynthesis. Involved in lifespan determination in a ubiquinone-independent manner.	Body	Island	Active promoter
cg26787244	965	6.96×10^−7	0.034	−0.55	11	644,011	—	—	Intergenic	Island	Weak transcribed
cg07261978	950	9.52×10^−7	0.038	−2.91	8	82,043,566	—	—	Intergenic	Open Sea	Weak enhancer
cg19029576	954	1.01×10^−6	0.040	−1.84	9	140,330,686	ENTPD8	Plays a central role in concentration of extracellular nucleotides.	Body	Island	Weak transcribed
cg25629796	971	1.12×10^−6	0.041	0.35	17	11,940,467	MAP2K4	Member of the mitogen-activated protein kinase (MAPK) family, which is involved in a wide variety of cellular processes such as proliferation, differentiation, transcription regulation, and development.	Body	Open Sea	Weak transcribed
cg01430385	954	1.20×10^−6	0.042	−1.43	10	1,033,349	GTPBP4	Involved in the biogenesis of the 60S ribosomal subunit.	Promoter	N. Shore	Active promoter
cg02519806	948	1.50×10^−6	0.045	−0.26	4	119,810,036	SYNPO2	Has an actin-binding and actin-bundling activity. Involved in regulation of cell migration; may be a tumor suppressor.	Promoter	Open Sea	Poised promoter
cg04915788	959	1.53×10^−6	0.044	−1.04	12	25,156,888	—	—	Intergenic	Open Sea	Weak enhancer
cg13170005	957	1.75×10^−6	0.048	−4.77	17	4,690,632	VMO1	Unknown function, homolog of a gene that encodes component of the outer membrane of the vitelline layer of the egg in chickens.	Promoter	N. Shore	Repressed
cg26483235	938	2.16×10^−6	0.050	−1.20	8	62,870,065	—	—	Intergenic	Open Sea	Heterochromatin
cg07318309	913	2.29×10^−6	0.051	0.19	5	156,537,633	HAVCR2	Cell surface receptor implicated in modulating innate and adaptive immune responses	Promoter	Open Sea	Weak enhancer
cg07509511	913	2.45×10^−6	0.051	−1.74	2	65,297,852	CEP68	Involved in maintenance of centrosome cohesion	Body	Open Sea	Transcription elongation
cg11062848	913	2.57×10^−6	0.052	2.08	6	13,616,537	NOL7	Maintains nucleolar structure and cell growth rates; also functions as a tumor suppressor and regulator of angiogenesis. Expression regulated by the RB tumor suppressor gene.	Body	S. Shore	Active promoter
cg18722557	917	2.80×10^−6	0.054	0.23	10	123,734,694	NSMCE4A	Component of a complex involved in DNA double-strand breaks by homologous recombination. Required for telomere maintenance via recombination in ALT (alternative lengthening of telomeres) cell lines.	Promoter	Island	Active promoter
cg20531550	923	2.81×10^−6	0.054	−3.11	7	128,550,902	KCP	Enhances bone morphogenetic protein (BMP) signaling in a paracrine manner.	Promoter	Island	Weak transcribed
cg05730283	933	2.96×10^−6	0.055	0.42	2	62,422,808	B3GNT2	Probably constitutes the main polylactosamine synthase	Promoter	Island	Active promoter
cg20540608	901	2.97×10^−6	0.053	1.53	5	60,241,003	NDUFAF2; ERCC8	NDUFAF2: Acts as a molecular chaperone for mitochondrial complex I assembly. ERCC8: Plays a role in DNA single-strand and double-strand breaks (DSSBs) repair; involved in repair of DSSBs by nonhomologous end joining (NHEJ)	Promoter	Island	Active promoter
cg25597976	929	3.00×10^−6	0.055	1.05	8	74,750,772	UBE2W	Promotes the ubiquitination of Fanconi anemia complementation group proteins and may be important in the repair of DNA damage	Body	Open Sea	Weak enhancer
AMPA
cg22171277	938	1.25×10^−7	0.038	0.32	2	70,363,119	—	—	Intergenic	Open Sea	Strong enhancer
cg23261070	946	1.55×10^−7	0.041	0.13	17	6,918,037	C17orf49	Component of chromatin complexes such as the MLL1/MLL and NURF complexes	Promoter	Island	Active promoter

Note: Linear models adjusted for urinary creatinine, age, race/ethnicity, BMI, smoking status, alcohol consumption, self-reported organic eating habits, diet quality (HEI), estimated WBC type proportions, batch, and position on chip, were fitted for all probes on the Illumina HumanMethylationEPIC array which remained after quality filtering. Results were pooled from 1,000 random subsamples of the training set (individuals per subsample). Probes that were statistically significant with FDR in 90% or more subsamples were considered differentially methylated. Missing data for covariates (from individuals in the training set) were substituted with the median (numeric variables) or mode (categorical variables). Chromatin state is predicted by ChromHMM from ENCODE data for the GM12878 cell line. Gene summaries are condensed from the RefSeq and UniProt gene summaries and provide a broad overview of the known function(s) of the gene product. —, no data; AMPA, aminomethylphosphonic acid; BMI, body mass index; Chr, Chromosome; FDR, false discovery rate; Freq., number of iterations in which DMP was identified as significant; HEI, Healthy Eating Index; WBC, white blood cell.

Figure 1.

Results from probe-level differential methylation analysis for urinary glyphosate in 332 postmenopausal California women. The volcano plot shows delta-M (difference in methylation M value for a 1-unit increase in the natural log of glyphosate) on the horizontal axis and on the vertical axis. Linear models adjusted for urinary creatinine, age, race/ethnicity, BMI, smoking status, alcohol consumption, self-reported organic eating habits, diet quality (Healthy Eating Index), estimated WBC type proportions, batch, and position on chip, were fitted for all probes on the Illumina HumanMethylationEPIC array, which remained after quality filtering. Results were pooled from 1,000 random subsamples of the training set ( individuals per subsample). Probes that were statistically significant with FDR in 90% or more subsamples are marked as hypermethylated or hypomethylated. Note: BMI, body mass index; FDR, false discovery rate; WBC, white blood cell.

Biological and functional analyses were performed on only glyphosate-associated probes because there were only two AMPA-associated probes (Figure 2). In comparison with all other probes on the array, a greater proportion of glyphosate-associated DMPs were located within CpG islands (33.3% vs. 19.2%), although this difference was not statistically significant (Fisher ). Glyphosate-associated DMPs were significantly enriched for enhancer regions (29.2% vs. 12.6%, Fisher ) and depleted for inactive/heterochromatin regions (4.2% vs. 36.8%, Fisher ). Three gene ontology terms were enriched in the 20 annotated genes containing glyphosate-associated probes: DNA metabolic process (biological process), intracellular organelle part (cellular component), and organelle part (cellular component), all with . No KEGG or IPA pathways were enriched in glyphosate-associated genes.

Figure 2.

Enrichment analysis for genomic context of glyphosate-associated differentially methylated probes (DMPs) identified in 332 postmenopausal California women. The proportions of glyphosate-associated DMPs compared to other probes on the array in each genomic context were compared with Fisher’s exact test; significant () -values are labeled.

Four regions for glyphosate and none for AMPA were significant in of subsamples and considered differentially methylated (Table 4). Three regions were significantly associated with AMPA at a relaxed threshold of of subsamples and are also included in Table 4. The glyphosate-associated regions were all located within gene promoters. Three were hypomethylated with greater glyphosate (MSH4, KCNA6, and ABAT), and the other was hypermethylated (NDUFAF2/ERCC8). The top AMPA-associated regions were all hypomethylated. Two were located within gene bodies (RNF39, TRIM31) and one within a gene promoter (ESR1).

Table 4

Differentially methylated regions (DMRs) associated with urinary glyphosate and AMPA concentration in 332 postmenopausal California women.

Chr.	Region start (hg19)	Region end (hg19)	Width (bases)	Number of CpGs in region	Freq.	Combined FDR q-value	Average methylation (β) difference per unit glyphosate/AMPA for probes in region	Gene(s)	Gene summary	Location within gene(s)
Glyphosate
1	76,262,302	76,262,984	682	9	979	0.023	−1.88	MSH4	Member of the DNA mismatch repair mutS family. Meiosis-specific protein that is required for reciprocal recombination and proper segregation of homologous chromosomes at meiosis I	Promoter
5	60,240,926	60,241,324	398	5	918	0.074	0.58	NDUFAF2; ERCC8	NDUFAF2: Acts as a molecular chaperone for mitochondrial complex I assembly. ERCC8: Plays a role in DNA single-strand and double-strand breaks (DSSBs) repair; involved in repair of DSSBs by nonhomologous end joining (NHEJ)	Promoter
12	4,918,169	4,919,230	1061	10	976	0.011	−5.60	KCNA6	Voltage-gated potassium channel that mediates transmembrane potassium transport in excitable membranes	Promoter
16	8,806,359	8,806,756	397	9	925	0.034	−4.98	ABAT	Responsible for catabolism of gamma-aminobutyric acid (GABA), an important, mostly inhibitory neurotransmitter in the central nervous system, into succinic semialdehyde	Promoter
AMPA
6	30,039,380	30,039,524	144	8	606	0.42	−13.02	RNF39	Lies within the major histocompatibility complex class I region on chromosome 6. Rat studies suggest that this gene encodes a protein that plays a role in an early phase of synaptic plasticity.	Body
6	30,071,496	30,071,612	116	5	636	0.025	−5.25	TRIM31	E3 ubiquitin-protein ligase. Shows altered expression in certain tumors and may be a negative regulator of cell growth.	Body
6	152,126,736	152,126,938	202	5	616	0.030	−3.63	ESR1	Estrogen receptor and ligand-activated transcription factor. Regulates the transcription of many estrogen-inducible genes that play a role in growth, metabolism, sexual development, gestation, and other reproductive functions and is expressed in many non-reproductive tissues. The receptor encoded by this gene plays a key role in breast cancer, endometrial cancer, and osteoporosis.	Promoter

Note: Linear models adjusted for urinary creatinine, age, race/ethnicity, BMI, smoking status, alcohol consumption, self-reported organic eating habits, diet quality (HEI), estimated WBC type proportions, batch, and position on chip, were fitted for all probes on the Illumina HumanMethylationEPIC array which remained after quality filtering. Results were pooled from 1,000 random subsamples of the training set ( individuals per subsample). Regions that were statistically significant with FDR in 90% or more subsamples were considered differentially methylated. Overlapping regions were combined, and FDR q-values were combined for all CpGs in region by Stouffer’s method. Missing data for covariates (from individuals in the training set) were substituted with the median (numeric variables) or mode (categorical variables). Gene summaries are condensed from the RefSeq and UniProt gene summaries and provide a broad overview of the known function(s) of the gene product. AMPA, aminomethylphosphonic acid; BMI, body mass index; Chr. Chromosome; FDR, false discovery rate; Freq., number of iterations out of 1,000 in which DMR was identified as significant; HEI, Healthy Eating Index; WBC, white blood cell.

Methylation Index

All 24 DMPs were selected by the elastic net for inclusion in the final glyphosate methylation index (Table S4). The index was significantly correlated with urinary glyphosate in the training (Pearson’s ; 95% CI: 0.55, 0.68, ; Figure 3A) and validation (; 95% CI: 0.10, 0.55, ; Figure 3B) sets. In the validation set, the methylation index was significantly associated with glyphosate tertile [median (IQR) ( to ) for lowest tertile vs. ( to ) for highest tertile, ANOVA ; Figure 3C] and showed excellent discrimination between the top and bottom tertiles of urinary glyphosate (, 95% CI: 0.57, 0.90; Figure 3D). A modified index trained using only probes included on the HumanMethylation450 BeadChip (the previous version of the Illumina methylation array), which selected 14 probes in the final model, achieved similar performance (Figure S5). The AMPA index was not significantly correlated with urinary AMPA in either the training or validation set (Figure S6).

Figure 3.

Performance of methylation index using 24 CpG sites to predict the natural logarithm of urinary glyphosate concentration in the training set (A) and the validation set (B, C). Panel C shows the methylation index was significantly associated with glyphosate tertile [median (IQR) ( to for lowest tertile vs. ( to ) for highest tertile, ANOVA ] in the validation set. Panel D shows the classification performance of the methylation index in the validation set for classifying the highest vs. the lowest tertile of urinary glyphosate. The index was developed using elastic net regression on methylation values of differentially methylated probes associated with glyphosate in the training set, a population of 332 postmenopausal California women. Note: ANOVA, analysis of variance; IQR, interquartile range.

Discussion

To our knowledge, this is the first study examining associations between urinary glyphosate and AMPA levels and DNA methylation. Higher AMPA, but not glyphosate, was associated with greater epigenetic age acceleration, a phenomenon which has previously been linked to the risk of all-cause mortality[56,57,58] and cancer.[59,60,94] We identified 24 CpG sites whose methylation was associated with glyphosate and two associated with AMPA. Four regions were associated with glyphosate, within the promoters of MSH4, KCNA6, ABAT, and NDUFAF2/ERCC8, and we found an association between ESR1 promoter hypomethylation and elevated AMPA. Finally, using 24 CpG sites, we developed a methylation index that was significantly associated with glyphosate concentration in an internal validation set. The AMPA index was not significantly associated with urinary AMPA in either set, likely due to the use of only two DMPs for prediction.

The significant and replicable differential DNA methylation associated with urinary glyphosate and AMPA informs the hypothesis that these compounds may have biological effects in humans, but the mechanisms by which glyphosate and AMPA could impact human health remain unclear. Although humans and other animals do not possess the shikimate pathway inhibited by glyphosate,[95,96,97] conversion of glyphosate to AMPA has been observed in multiple types of bacteria.[98] Previous studies have shown glyphosate-induced microbiome changes in honeybees,[99,100] birds,[101] and rats.[24,25,26,27] Because the microbiome can have significant and wide-ranging impacts on human health,[102] it is possible that glyphosate and/or AMPA could influence various aspects of human health via perturbations of the microbiome, changes that could be reflected in differential peripheral blood DNA methylation.[103,104]

In our study, AMPA was associated with epigenetic age acceleration, which has been previously associated with other environmental exposures[61,62,63,64] and risk of many diseases, including breast cancer,[59] B-cell lymphoma,[60] lung cancer,[96] obesity and metabolic syndrome,[105,106,107] and all-cause mortality.[56,57,58] The fact that this relationship was present only for AMPA supports the hypothesis that glyphosate and AMPA have distinct effects on the human body. This association was present for all epigenetic clocks considered (Hannum, Horvath, and Levine), although adjustment for age, race/ethnicity, BMI, smoking status, alcohol consumption, self-reported organic eating habits, and diet quality attenuated the relationship for the Levine clock, which incorporates phenotypic information and thus may be more influenced by these covariates. However, we observed no statistically significant relationship between epigenetic pace of aging and glyphosate or AMPA concentrations. Pace of aging aims to measure a distinct biological phenomenon in comparison with epigenetic age acceleration (rate of change of biological age vs. differences between chronological and biological age).

Furthermore, the genes whose methylation was associated with glyphosate and AMPA included those involved in various biological pathways related to cancer. SF3B2 is involved in RNA splicing and DNA repair[108] and associated with various types of cancer.[109,110,111] A germline mutation in the mismatch repair gene MSH4[112] was described in a family with a high incidence of nervous system tumors.[113] ERCC8 is also involved in DNA repair via transcription-coupled nucleotide excision repair and double-strand break repair.[114,115,116] Hypomethylation at TRIM31, which has been shown to promote progression in a variety of tumor types,[117,118,119] was associated with AMPA. The AMPA-associated hypomethylation at the ESR1 promoter is fascinating, given the potential link between AMPA and breast cancer risk[23,32] and the potential for glyphosate-induced endocrine disruption.[20] Perturbed expression of the ESR1 product in response to glyphosate exposure has been observed in in vitro studies of breast cancer cells.[23] Rat studies have also demonstrated altered expression accompanying changes in promoter methylation.[48,120] However, it should be noted that the relationship between AMPA and ESR1 promoter hypomethylation was only present after relaxing the original threshold for significance.

To our knowledge, this is also the largest report of urinary glyphosate and AMPA levels in the nonagricultural setting in the United States, complementing previous studies in the general population in other countries[121,122,123] and in limited settings in the United States.[11,12,124,125] The majority of women () had detectable glyphosate and AMPA in at least one urine sample, highlighting the near-ubiquitous exposure to these compounds, although the small convenience sample may not be representative of the general population. The levels observed among participants in our study (median of for glyphosate and for AMPA) were slightly lower than those previously reported, likely due to differences in the study population and detection assays. Our LC-MS/MS assay had a considerably higher analytic sensitivity (LOD/LOQ for glyphosate and for AMPA) than most assays previously described.[126,127]

This study used the most recent methylation array chip and best practices for methylation array data filtering, normalization, and analysis alongside a resampling-based approach intended to improve the stability and reproducibility of the results.[88,89,90,91] We previously demonstrated that this method is capable of identifying differentially methylated sites associated with smoking with a minimal number of false positives in comparison with a traditional epigenome-wide association approach.[86]

However, the lack of an external validation cohort is a limitation of the study. To our knowledge, no other studies have measured both urinary glyphosate and AMPA and DNA methylation, so an external data set is not currently available. The successful performance of the methylation index for glyphosate concentration in the internal validation set suggests that the glyphosate-associated differential methylation may be replicated in other cohorts with a similar participant profile. Our study cohort consisted of postmenopausal women age 45–66 y, so results may not be generalizable to other populations or those residing outside California. The power of these analyses may be limited by the relatively low exposure levels in this population. A further limitation of our study is the lack of gene expression data. Without this, we cannot state whether the differential DNA methylation identified in our study translates to differences in gene expression that could confer functional impacts in the body, nor can we infer the temporal relationship between DNA methylation differences and glyphosate exposure, given the cross-sectional nature of the study. Finally, although we did not observe a relationship between recent herbicide use and glyphosate and AMPA concentrations, those data were missing for 37% of participants, which should be considered another limitation of the study.

DNA methylation is highly dependent on tissue type.[128] Given the ambiguity of glyphosate’s impact on health and the systemic nature of exposure, blood likely represents the best available tissue at this time, especially considering its availability for epidemiological study. A major concern regarding the use of blood in epigenome-wide association studies is its heterogeneous cell composition, which may be impacted by disease states or exposure to pro-inflammatory compounds.[129] However, we saw no evidence of differences in WBC composition by glyphosate or AMPA levels (Table S2), and all analyses were adjusted for these proportions.

In humans, glyphosate and AMPA have short half-lives of approximately 5–10 h,[130,131] which makes accurate assessment of long-term exposure challenging. Our choice to use two spot urine samples from a 10-d period was intended to balance the need for a more complete picture of typical exposure with feasibility. The moderate between-sample ICCs of 0.53 for glyphosate and 0.34 for AMPA in our study participants highlight the pitfalls of relying on a single urine sample and suggest that more samples per individual may be needed to provide a better estimate of long-term exposure. However, because the dates of spot urine samples were not known in advance to the study participants, the glyphosate and AMPA measurements should reflect their usual habits.

This study identified differential DNA methylation associated with the herbicide glyphosate and its metabolite AMPA and developed a methylation index that accurately predicted urinary glyphosate concentration and tertile in an internal validation sample. Glyphosate- and AMPA-associated methylation occurred near genes associated with cancer (SF3B2,[109,110,111] MSH4,[113] and TRIM31[117,118,119]) and endocrine disruption (ESR1[132]), and AMPA was associated with greater epigenetic age acceleration. These results suggest that exposure to these common chemicals affects the epigenome, informing the hypothesis that glyphosate and/or AMPA exposure might elevate the risk for disease, including cancer. Further studies are warranted to replicate our results, determine the functional impact of glyphosate- and AMPA-associated differential DNA methylation, and explore whether DNA methylation could serve as a biomarker of long-term glyphosate exposure.

Click here for additional data file.

Click here for additional data file.