Study / S3Y6NAR272024-02-26

DNA methylation-environment interactions in the human genome

Rachel A. Johnston, Katherine A. Aracena, Luis B. Barreiro, Amanda J. Lea, Jenny Tung

About this study

Previously, we showed that a massively parallel reporter assay, mSTARR-seq, could be used to simultaneously test for both enhancer-like activity and DNA methylation-dependent enhancer activity for millions of loci in a single experiment (Lea et al., 2018). Here, we apply mSTARR-seq to query nearly the entire human genome, including almost all CpG sites profiled either on the commonly used Illumina Infinium MethylationEPIC array or via reduced representation bisulfite sequencing. We show that fragments containing these sites are enriched for regulatory capacity, and that methylation-dependent regulatory activity is in turn sensitive to the cellular environment. In particular, regulatory responses to interferon alpha (IFNA) stimulation are strongly attenuated by methyl marks, indicating widespread DNA methylation-environment interactions. In agreement, methylation-dependent responses to IFNA identified via mSTARR-seq predict methylation-dependent transcriptional responses to challenge with influenza virus in human macrophages. Our observations support the idea that pre-existing DNA methylation patterns can influence the response to subsequent environmental exposures-one of the tenets of biological embedding. However, we also find that, on average, sites previously associated with early life adversity are not more likely to functionally influence gene regulation than expected by chance.

Full author list & citation

Rachel A. Johnston, Katherine A. Aracena, Luis B. Barreiro, Amanda J. Lea, Jenny Tung. DNA methylation-environment interactions in the human genome. 2024-02-26. https://doi.org/10.7554/eLife.89371

Experiments 4

E3OOE73G8

K562 mSTARR-seq baseline: methylated versus sham-treated library

Captured human genomic DNA from GM12878 was cloned into the episomal pmSTARRseq1 reporter, treated with M.SssI or a sham reaction, and transfected into K562 cells. The table contains model statistics for 216,091 study-QC-passed windows under baseline/vehicle conditions.

Standard STARR-seqHumanhg38
Explore data
E4ONEXPBI

K562 mSTARR-seq after dexamethasone stimulation

The captured GM12878 genomic mSTARR-seq library was treated with M.SssI or sham-treated, transfected into K562, and exposed to 1 µM dexamethasone 42 hours after transfection. The table contains model statistics for 250,569 study-QC-passed windows harvested 6 hours after stimulation.

Standard STARR-seqHumanhg38
Explore data
E6E0ICW4Y

K562 mSTARR-seq after IFNA2b stimulation

The captured GM12878 genomic mSTARR-seq library was treated with M.SssI or sham-treated, transfected into K562, and exposed to 2,000 U/mL IFNA2b 42 hours after transfection. The table contains model statistics for 193,948 study-QC-passed windows harvested 6 hours after stimulation.

Standard STARR-seqHumanhg38
Explore data

Raw source data 9 files

Original supplemental and deposited inputs retained for this study. Download files individually or together as a ZIP; nested folders are preserved. Source reuse terms apply, and sequencing reads may be omitted.

Download all 9 files (ZIP)elife-89371-supp1-v1.xlsxelife-89371-supp10-v1.xlsxelife-89371-supp11-v1.xlsxelife-89371-supp2-v1.xlsxelife-89371-supp21-v1.xlsxelife-89371-supp3-v1.xlsxelife-89371-supp7-v1.xlsxelife-89371-supp9-v1.xlsxSOURCE_MANIFEST.txt

Cite OpenMPRA

Cite the OpenMPRA database. Include your access date because the collection changes over time.

Please also cite the source studies when using their data.