Longitudinal data reveal strong genetic and weak non-genetic components of ethnicity-dependent blood DNA methylation levels

Mid-life anti-inflammatory metabolites are inversely associated with long-term cardiovascular disease events

BACKGROUND

Preclinical data have shown that low levels of metabolites with anti-inflammatory properties may impact metabolic disease processes. However, the association between mid-life levels of such metabolites and long-term ASCVD risk is not known.

METHODS

We characterised the plasma metabolomic profile (1228 metabolites) of 1852 participants (58.1 ± 7.5 years old, 69.6% female, 43.6% self-identified as Black) enrolled in the Heart Strategies Concentrating on Risk Evaluation (Heart SCORE) study. Logistic regression was used to assess the impact of metabolite levels on ASCVD risk (nonfatal MI, revascularisation, and cardiac mortality). We additionally explored the effect of genetic variants neighbouring ASCVD-related genes on the levels of metabolites predictive of ASCVD events. The Atherosclerosis Risk in Communities (ARIC) study (n = 4790; 75.5 ± 5.1 years old, 57.4% female, 19.5% self-identified as Black) was used as an independent validation cohort.

FINDINGS

In fully adjusted models, alpha-ketobutyrate [AKB] (OR 0.62 [95% CI, 0.49-0.80]; p < 0.001), and 1-palmitoyl-2-linoleoyl-GPI [OR, 0.62, 95% CI, 0.47-0.83; p < 0.001], two metabolites in amino acid and phosphatidylinositol lipid pathways, respectively, showed a significant protective association with incident ASCVD risk in both Heart SCORE and ARIC cohorts. Three plasmalogens and a bilirubin derivative, whose levels were regulated by genetic variants neighbouring FADS1 and UGT1A1, respectively, exhibited a significant protective association with ASCVD risk in the Heart SCORE only.

INTERPRETATION

Higher mid-life levels of AKB and 1-palmitoyl-2-linoleoyl-GPI metabolites may be associated with lower risk late-life ASCVD events. Further research can determine the causality and therapeutic potential of these metabolites in ASCVD.

FUNDING

This study was funded by the Pennsylvania Department of Health (ME-02-384). The department specifically disclaims responsibility for any analyses, interpretations, or conclusions. Additional funding was provided by National Institutes of Health (NIH) grant R01HL089292 and UL1 TR001857 (Steven Reis). Further, NIH funded R01HL141824 and R01HL168683 were used for the ARIC study validation (Bing Yu).

Map of epigenetic age acceleration: A worldwide analysis

We present a systematic analysis of epigenetic age acceleration based on by far the largest collection of publicly available DNA methylation data for healthy samples (93 datasets, 23 K samples), focusing on the geographic (25 countries) and ethnic (31 ethnicities) aspects around the world. We employed the most popular epigenetic tools for assessing age acceleration and examined their quality metrics and ability to extrapolate to epigenetic data from different tissue types and age ranges different from the training data of these models. In most cases, the models proved to be inconsistent with each other and showed different signs of age acceleration, with the PhenoAge model tending to systematically underestimate and different versions of the GrimAge model tending to systematically overestimate the age prediction of healthy subjects. Referring to data availability and consistency, most countries and populations are still not represented in GEO, moreover, different datasets use different criteria for determining healthy controls. Because of this, it is difficult to fully isolate the contribution of "geography/environment", "ethnicity" and "healthiness" to epigenetic age acceleration. Among the explored metrics, only the DunedinPACE, which measures aging rate, appears to adequately reflect the standard of living and socioeconomic indicators in countries, although it has a limited application to blood methylation data only. Invariably, by epigenetic age acceleration, males age faster than females in most of the studied countries and populations.

Epigenomic insights into common human disease pathology

The epigenome-the chemical modifications and chromatin-related packaging of the genome-enables the same genetic template to be activated or repressed in different cellular settings. This multi-layered mechanism facilitates cell-type specific function by setting the local sequence and 3D interactive activity level. Gene transcription is further modulated through the interplay with transcription factors and co-regulators. The human body requires this epigenomic apparatus to be precisely installed throughout development and then adequately maintained during the lifespan. The causal role of the epigenome in human pathology, beyond imprinting disorders and specific tumour suppressor genes, was further brought into the spotlight by large-scale sequencing projects identifying that mutations in epigenomic machinery genes could be critical drivers in both cancer and developmental disorders. Abrogation of this cellular mechanism is providing new molecular insights into pathogenesis. However, deciphering the full breadth and implications of these epigenomic changes remains challenging. Knowledge is accruing regarding disease mechanisms and clinical biomarkers, through pathogenically relevant and surrogate tissue analyses, respectively. Advances include consortia generated cell-type specific reference epigenomes, high-throughput DNA methylome association studies, as well as insights into ageing-related diseases from biological 'clocks' constructed by machine learning algorithms. Also, 3rd-generation sequencing is beginning to disentangle the complexity of genetic and DNA modification haplotypes. Cell-free DNA methylation as a cancer biomarker has clear clinical utility and further potential to assess organ damage across many disorders. Finally, molecular understanding of disease aetiology brings with it the opportunity for exact therapeutic alteration of the epigenome through CRISPR-activation or inhibition.

A consistent pattern of slide effects in Illumina DNA methylation BeadChip array data

Background: Recent studies have identified thousands of associations between DNA methylation CpGs and complex diseases/traits, emphasizing the critical role of epigenetics in understanding disease aetiology and identifying biomarkers. However, association analyses based on methylation array data are susceptible to batch/slide effects, which can lead to inflated false positive rates or reduced statistical powerResults: We use multiple DNA methylation datasets based on the popular Illumina Infinium MethylationEPIC BeadChip array to describe consistent patterns and the joint distribution of slide effects across CpGs, confirming and extending previous results. The susceptible CpGs overlap with the Illumina Infinium HumanMethylation450 BeadChip array content.Conclusions: Our findings reveal systematic patterns in slide effects. The observations provide further insights into the characteristics of these effects and can improve existing adjustment approaches.

Body Mass Index, Adverse Pregnancy Outcomes, and Cardiovascular Disease Risk

BACKGROUND

Obesity is a well-established risk factor for both adverse pregnancy outcomes (APOs) and cardiovascular disease (CVD). However, it is not known whether APOs are mediators or markers of the obesity-CVD relationship. This study examined the association between body mass index, APOs, and postpartum CVD risk factors.

METHODS

The sample included adults from the nuMoM2b (Nulliparous Pregnancy Outcomes Study: Monitoring Mothers-To-Be) Heart Health Study who were enrolled in their first trimester (6 weeks-13 weeks 6 days gestation) from 8 United States sites. Participants had a follow-up visit at 3.7 years postpartum. APOs, which included hypertensive disorders of pregnancy, preterm birth, small-for-gestational-age birth, and gestational diabetes, were centrally adjudicated. Mediation analyses estimated the association between early pregnancy body mass index and postpartum CVD risk factors (hypertension, hyperlipidemia, and diabetes) and the proportion mediated by each APO adjusted for demographics and baseline health behaviors, psychosocial stressors, and CVD risk factor levels.

RESULTS

Among 4216 participants enrolled, mean±SD maternal age was 27±6 years. Early pregnancy prevalence of overweight was 25%, and obesity was 22%. Hypertensive disorders of pregnancy occurred in 15%, preterm birth in 8%, small-for-gestational-age birth in 11%, and gestational diabetes in 4%. Early pregnancy obesity, compared with normal body mass index, was associated with significantly higher incidence of postpartum hypertension (adjusted odds ratio, 1.14 [95% CI, 1.10-1.18]), hyperlipidemia (1.11 [95% CI, 1.08-1.14]), and diabetes (1.03 [95% CI, 1.01-1.04]) even after adjustment for baseline CVD risk factor levels. APOs were associated with higher incidence of postpartum hypertension (1.97 [95% CI, 1.61-2.40]) and hyperlipidemia (1.31 [95% CI, 1.03-1.67]). Hypertensive disorders of pregnancy mediated a small proportion of the association between obesity and incident hypertension (13% [11%-15%]) and did not mediate associations with incident hyperlipidemia or diabetes. There was no significant mediation by preterm birth or small-for-gestational-age birth.

CONCLUSIONS

There was heterogeneity across APO subtypes in their association with postpartum CVD risk factors and mediation of the association between early pregnancy obesity and postpartum CVD risk factors. However, only a small or nonsignificant proportion of the association between obesity and CVD risk factors was mediated by any of the APOs, suggesting APOs are a marker of prepregnancy CVD risk and not a predominant cause of postpartum CVD risk.

An integrative framework and recommendations for the study of DNA methylation in the context of race and ethnicity

Human social epigenomics research is critical to elucidate the intersection of social and genetic influences underlying racial and ethnic differences in health and development. However, this field faces major challenges in both methodology and interpretation with regard to disentangling confounded social and biological aspects of race and ethnicity. To address these challenges, we discuss how these constructs have been approached in the past and how to move forward in studying DNA methylation (DNAm), one of the best-characterized epigenetic marks in humans, in a responsible and appropriately nuanced manner. We highlight self-reported racial and ethnic identity as the primary measure in this field, and discuss its implications in DNAm research. Racial and ethnic identity reflects the biological embedding of an individual's sociocultural experience and environmental exposures in combination with the underlying genetic architecture of the human population (i.e., genetic ancestry). Our integrative framework demonstrates how to examine DNAm in the context of race and ethnicity, while considering both intrinsic factors-including genetic ancestry-and extrinsic factors-including structural and sociocultural environment and developmental niches-when focusing on early-life experience. We reviewed DNAm research in relation to health disparities given its relevance to race and ethnicity as social constructs. Here, we provide recommendations for the study of DNAm addressing racial and ethnic differences, such as explicitly acknowledging the self-reported nature of racial and ethnic identity, empirically examining the effects of genetic variants and accounting for genetic ancestry, and investigating race-related and culturally regulated environmental exposures and experiences.

Supplementary Information

The online version contains supplementary material available at 10.1007/s44155-023-00039-z.

A longitudinal epigenome-wide association study of preeclamptic and normotensive pregnancy

Background

While preeclampsia (PE) is a leading cause of pregnancy-related morbidity/mortality, its underlying mechanisms are not fully understood. DNA methylation (DNAm) is a dynamic regulator of gene expression that may offer insight into PE pathophysiology and/or serve as a biomarker (e.g., risk, subtype, a therapeutic response). This study's purpose was to evaluate for differences in blood-based DNAm across all trimesters between individuals eventually diagnosed with PE (cases) and individuals who remained normotensive throughout pregnancy, did not develop proteinuria, and birthed a normally grown infant (controls).

Results

In the discovery phase, longitudinal, genome-wide DNAm data were generated across three trimesters of pregnancy in 56 participants (n=28 cases, n=28 controls) individually matched on self-identified race, pre-pregnancy body mass index, smoking, and gestational age at sample collection. An epigenome-wide association study (EWAS) was conducted, using surrogate variable analysis to account for unwanted sources of variation. No CpGs met the genome-wide significance p value threshold of 9×10-8, but 16 CpGs (trimester 1: 5; trimester 2: 1; trimester 3: 10) met the suggestive significance threshold of 1×10-5. DNAm data were also evaluated for differentially methylated regions (DMRs) by PE status. Three DMRs in each trimester were significant after Bonferonni-adjustment. Since only third-trimester samples were available from an independent replication sample (n=64 cases, n=50 controls), the top suggestive hits from trimester 3 (cg16155413 and cg21882990 associated with TRAF3IP2-AS1/TRAF3IP2 genes, which also made up the top DMR) were carried forward for replication. During replication, DNAm data were also generated for validation purposes from discovery phase third trimester samples. While significant associations between DNAm and PE status were observed at both sites in the validation sample, no associations between DNAm and PE status were observed in the independent replication sample.

Conclusions

The discovery phase findings for cg16155413/cg21882990 (TRAF3IP2-AS1/TRAF3IP2) were validated with a new platform but were not replicated in an independent sample. Given the differences in participant characteristics between the discovery and replication samples, we cannot rule out important signals for these CpGs. Additional research is warranted for cg16155413/cg21882990, as well as top hits in trimesters 1-2 and significant DMRs that were not examined in the replication phase.

Multi-omic association study identifies DNA methylation-mediated genotype and smoking exposure effects on lung function in children living in urban settings

Impaired lung function in early life is associated with the subsequent development of chronic respiratory disease. Most genetic associations with lung function have been identified in adults of European descent and therefore may not represent those most relevant to pediatric populations and populations of different ancestries. In this study, we performed genome-wide association analyses of lung function in a multiethnic cohort of children (n = 1,035) living in low-income urban neighborhoods. We identified one novel locus at the TDRD9 gene in chromosome 14q32.33 associated with percent predicted forced expiratory volume in one second (FEV1) (p = 2.4x10-9; βz = -0.31, 95% CI = -0.41- -0.21). Mendelian randomization and mediation analyses revealed that this genetic effect on FEV1 was partially mediated by DNA methylation levels at this locus in airway epithelial cells, which were also associated with environmental tobacco smoke exposure (p = 0.015). Promoter-enhancer interactions in airway epithelial cells revealed chromatin interaction loops between FEV1-associated variants in TDRD9 and the promoter region of the PPP1R13B gene, a stimulator of p53-mediated apoptosis. Expression of PPP1R13B in airway epithelial cells was significantly associated the FEV1 risk alleles (p = 1.3x10-5; β = 0.12, 95% CI = 0.06-0.17). These combined results highlight a potential novel mechanism for reduced lung function in urban youth resulting from both genetics and smoking exposure.

Evolutionary Implications of Environmental Toxicant Exposure

Homo sapiens have been exposed to various toxins and harmful compounds that change according to various phases of human evolution. Population genetics studies showed that such exposures lead to adaptive genetic changes; while observing present exposures to different toxicants, the first molecular mechanism that confers plasticity is epigenetic remodeling and, in particular, DNA methylation variation, a molecular mechanism proposed for medium-term adaptation. A large amount of scientific literature from clinical and medical studies revealed the high impact of such exposure on human biology; thus, in this review, we examine and infer the impact that different environmental toxicants may have in shaping human evolution. We first describe how environmental toxicants shape natural human variation in terms of genetic and epigenetic diversity, and then we describe how DNA methylation may influence mutation rate and, thus, genetic variability. We describe the impact of these substances on biological fitness in terms of reproduction and survival, and in conclusion, we focus on their effect on brain evolution and physiology.

Blood DNA methylation signatures of lifestyle exposures: tobacco and alcohol consumption

BACKGROUND

Smoking and alcohol consumption may compromise health by way of epigenetic modifications. Epigenetic signatures of alcohol and tobacco consumption could provide insights into the reversibility of phenotypic changes incurred with differing levels of lifestyle exposures. This study describes and validates two novel epigenetic signatures of tobacco (EpiTob) and alcohol (EpiAlc) consumption and investigates their association with disease outcomes.

METHODS

The epigenetic signatures, EpiTob and EpiAlc, were developed using data from the Swiss Kidney Project on Genes in Hypertension (SKIPOGH) (N = 689). Epigenetic and phenotypic data available from the 1921 (N = 550) and 1936 (N = 1091) Lothian Birth Cohort (LBC) studies, and two publicly available datasets on GEO Accession (GSE50660, N = 464; and GSE110043, N = 94) were used to validate the signatures. A multivariable logistic regression model, adjusting for age and sex, was used to assess the association between self-reported tobacco or alcohol consumption and the respective epigenetic signature, as well as to estimate the association between CVD and epigenetic signatures. A Cox proportional hazard model was used to estimate the risk of mortality in association with the EpiTob and EpiAlc signatures.

RESULTS

The EpiTob signature was positively associated with self-reported tobacco consumption for current or never smokers with explained variance ranging from 0.49 (LBC1921) to 0.72 (LBC1936) (pseudo-R2). In the SKIPOGH, LBC1921 and LBC1936 cohorts, the epigenetic signature for alcohol consumption explained limited variance in association with self-reported alcohol status [i.e., non-drinker, moderate drinker, and heavy drinker] (pseudo-R2 = 0.05, 0.03 and 0.03, respectively), although this improved considerably when measuring self-reported alcohol consumption with standardized units consumed per week (SKIPOGH R2 = 0.21; LBC1921 R2 = 0.31; LBC1936 R2 = 0.41). Both signatures were associated with history of CVD in SKIPOGH and LBC1936, but not in LBC1921. The EpiTob signature was associated with increased risk of all-cause and lung-cancer specific mortality in the 1936 and 1921 LBC cohorts.

CONCLUSIONS

This study found the EpiTob and EpiAlc signatures to be well-correlated with self-reported exposure status and associated with long-term health outcomes. Epigenetic signatures of lifestyle exposures may reduce measurement issues and biases and could aid in risk stratification for informing early-stage targeted interventions.

Dynamics of Methylation of CpG Sites Associated With Systemic Lupus Erythematosus Subtypes in a Longitudinal Cohort

OBJECTIVE

Findings from cross-sectional studies have revealed associations between DNA methylation and systemic lupus erythematosus (SLE) outcomes. This study was undertaken to investigate the dynamics of DNA methylation by examining participants from an SLE longitudinal cohort using samples collected at 2 time points.

METHODS

A total of 101 participants from the California Lupus Epidemiology Study were included in our analysis. DNA was extracted from blood samples collected at the time of enrolment in the cohort and samples collected after 2 years and was analyzed using Illumina EPIC BeadChip kit. Paired t-tests were used to identify genome-wide changes which included 256 CpG sites previously found to be associated with SLE subtypes. Linear mixed models were developed to understand the relationship between DNA methylation and disease activity, medication use, and sample cell-type proportions, adjusted for age, sex, and genetic principal components.

RESULTS

The majority of CpGs that were previously determined to be associated with SLE subtypes remained stable over 2 years (185 CpGs [72.3%]; t-test false discovery rate >0.05). Compared to background genome-wide methylation, there was an enrichment of SLE subtype-associated CpGs that changed over time (27.7% versus 0.34%). Changes in cell-type proportions were associated with changes at 67 CpGs (P < 2.70 × 10-5 ), and 15 CpGs had at least 1 significant association with immunosuppressant use.

CONCLUSION

In this longitudinal SLE cohort, we identified a subset of SLE subtype-associated CpGs that remained stable over time and may be useful as biomarkers of disease subtypes. Another subset of SLE subtype-associated CpGs changed at a higher proportion compared to the genome-wide methylome. Additional studies are needed to understand the etiology and impact of these changes on methylation of SLE-associated CpGs.

Sensitisation to Pollen Allergens in Children and Adolescents of Different Ancestry Born and Living in the Same Area

Background

Allergy can start at early ages, with genetic and environmental factors contributing to its development.

Aim

The study aimed to describe the pattern of sensitisation and allergy in children and adolescents of Spanish versus Moroccan ancestry but born in the same rural area of Spain.

Methods

Participants were children and adolescents (3–19 years) of Spanish or Moroccan descent, born in Blanca, Murcia (Spain). A detailed questionnaire was completed, and skin prick tests were performed to assess reactions to the most prevalent pollen allergens (O. europaea, P. pratense, S. kali, C. arizonica, P. acerifolia, A. vulgaris and P. judaica) plus molecular components Ole e 1 and Ole e 7. The association with ancestry was verified by studying participants’ parents.

Results

The study included 693 participants: 48% were aged 3–9 years and 52%, 10–19 years; 80% were of Spanish descent and 20% of Moroccan descent. Sensitisation to Olea europaea, Phleum pratense, Salsola kali and Cupressus arizonica were slightly higher in the Spanish group. The only significant differences were observed in sensitisation to Ole e 1 (p=0.02). Rhinitis, conjunctivitis, and rhinitis plus asthma were significantly higher in the Spanish group (p=0.03, p=0.02, p=0.007, respectively). The sensitisation pattern differed between Spanish and Moroccan parents, and between Moroccan parents and their children, but not between Spanish parents and their children.

Conclusion

Both environment and ancestry may influence sensitisation and symptoms. Although the environment seems to have a stronger influence, other factors may contribute to the differences in prevalence and in the clinical entities in people of Spanish versus Moroccan descent.

A biosocial return to race? A cautionary view for the postgenomic era

Recent studies demonstrating epigenetic and developmental sensitivity to early environments, as exemplified by fields like the Developmental Origins of Health and Disease (DOHaD) and environmental epigenetics, are bringing new data and models to bear on debates about race, genetics, and society. Here, we first survey the historical prominence of models of environmental determinism in early formulations of racial thinking to illustrate how notions of direct environmental effects on bodies have been used to naturalize racial hierarchy and inequalities in the past. Next, we conduct a scoping review of postgenomic work in environmental epigenetics and DOHaD that looks at the role of race/ethnicity in human health (2000-2021). Although there is substantial heterogeneity in how race is conceptualized and interpreted across studies, we observe practices that may unwittingly encourage typological thinking, including: using DNA methylation as a novel marker of racial classification; neglect of variation and reversibility within supposedly homogenous racial groups; and a tendency to label and reify whole groups as pathologized or impaired. Even in the very different politico-economic and epistemic context of contemporary postgenomic science, these trends echo deeply held beliefs in Western thinking which claimed that different environments shape different bodies and then used this logic to argue for essential differences between Europeans and non-Europeans. We conclude with a series of suggestions on interpreting and reporting findings in these fields that we feel will help researchers harness this work to benefit disadvantaged groups while avoiding the inadvertent dissemination of new and old forms of stigma or prejudice.

SHOX2 methylation in Vietnamese patients with lung cancer

BACKGROUND

DNA methylation on cytosine in the CpG dinucleotides is one of the most common epigenetic perturbations taking place during cancer initiation, progression, occurrence and resistance therapy. DNA methylation seems to be sufficiently stable epigenetic modification to be utilized as a cancer biomarker in in vitro diagnostic (IVD) settings. Nowadays, the SHOX2 methylation (mSHOX2) is one of the most valuable DNA methylation biomarkers of lung cancer that is the leading cause of cancer death. It is being continuously validated across ethnicities, lifestyles and lifespan. This study focused on characteristics of mSHOX2 in Vietnamese patients with lung cancer since a lack of investigation and evidence of its utility in this country.

METHODS

The probe and primer sets were designed according to the MethyLight method for quantitative assessment of the mSHOX2 in 214 formalin-fixed paraffin-embedded (FFPE) lung tissues and 57 plasma samples.

RESULTS

mSHOX2 in FFPE tissues allowed discriminating benign and malignant lung diseases with 60% (95% CI 50.7-68.8%) sensitivity and 90.4% (95% CI 82.6-95.5%) specificity. Importantly, based on mSHOX2 in plasma, lung cancer could be detected with 83.3% (95% CI 65.3-94.4%) sensitivity and 92.6% (95% CI 75.7-99.1%) specificity, respectively. There were insignificant associations between mSHOX2 with age, cancer stage, EGFR mutation and serum CEA, CYFRA21-1 concentrations except for that gender.

CONCLUSION

Our study indicated that mSHOX2 was satisfactory for distinguishing malignant from benign lung tissue and noninvasively detecting lung cancer.

Novel diagnostic DNA methylation episignatures expand and refine the epigenetic landscapes of Mendelian disorders

Overlapping clinical phenotypes and an expanding breadth and complexity of genomic associations are a growing challenge in the diagnosis and clinical management of Mendelian disorders. The functional consequences and clinical impacts of genomic variation may involve unique, disorder-specific, genomic DNA methylation episignatures. In this study, we describe 19 novel episignature disorders and compare the findings alongside 38 previously established episignatures for a total of 57 episignatures associated with 65 genetic syndromes. We demonstrate increasing resolution and specificity ranging from protein complex, gene, sub-gene, protein domain, and even single nucleotide-level Mendelian episignatures. We show the power of multiclass modeling to develop highly accurate and disease-specific diagnostic classifiers. This study significantly expands the number and spectrum of disorders with detectable DNA methylation episignatures, improves the clinical diagnostic capabilities through the resolution of unsolved cases and the reclassification of variants of unknown clinical significance, and provides further insight into the molecular etiology of Mendelian conditions.

Differential DNA Methylation by Hispanic Ethnicity Among Firefighters in the United States

Firefighters are exposed to a variety of environmental hazards and are at increased risk for multiple cancers. There is evidence that risks differ by ethnicity, yet the biological or environmental differences underlying these differences are not known. DNA methylation is one type of epigenetic regulation that is altered in cancers. In this pilot study, we profiled DNA methylation with the Infinium MethylationEPIC in blood leukocytes from 31 Hispanic white and 163 non-Hispanic white firefighters. We compared DNA methylation (1) at 12 xenobiotic metabolizing genes and (2) at all loci on the array (>740 000), adjusting for confounders. Five of the xenobiotic metabolizing genes were differentially methylated at a raw P-value <.05 when comparing the 2 ethnic groups, yet were not statistically significant at a 5% false discovery rate (q-value <.05). In the epigenome-wide analysis, 76 loci exhibited DNA methylation differences at q < .05. Among these, 3 CpG sites in the promoter region of the biotransformation gene SULT1C2 had lower methylation in Hispanic compared to non-Hispanic firefighters. Other differentially methylated loci included genes that have been implicated in carcinogenesis in published studies (FOXK2, GYLTL1B, ZBTB16, ARHGEF10, and more). In this pilot study, we report differential DNA methylation between Hispanic and non-Hispanic firefighters in xenobiotic metabolism genes and other genes with functions related to cancer. Epigenetic susceptibility by ethnicity merits further study as this may alter risk for cancers linked to toxic exposures.