Harmonization of transcriptomic and methylomic analysis in environmental epidemiology studies for potential application in chemical risk assessment.

Recent efforts have posited the utility of transcriptomic-based approaches to understand chemical-related perturbations in the context of human health risk assessment. Epigenetic modification (e.g., DNA methylation) can influence gene expression changes and is known to occur as a molecular response to some chemical exposures. Characterization of these methylation events is critical to understand the molecular consequences of chemical exposures. In this context, a novel workflow was developed to interrogate publicly available epidemiological transcriptomic and methylomic data to identify relevant pathway level changes in response to chemical exposure, using inorganic arsenic as a case study. Gene Set Enrichment Analysis (GSEA) was used to identify causal methylation events that result in concomitant downstream transcriptional deregulation. This analysis demonstrated an unequal distribution of differentially methylated regions across the human genome. After mapping these events to known genes, significant enrichment of a subset of these pathways suggested that arsenic-mediated methylation may be both specific and non-specific. Parallel GSEA performed on matched transcriptomic samples determined that a substantially reduced subset of these pathways are enriched and that not all chemically-induced methylation results in a downstream alteration in gene expression. The resulting pathways were found to be representative of well-established molecular events known to occur in response to arsenic exposure. The harmonization of enriched transcriptional patterns with those identified from the methylomic platform promoted the characterization of plausibly causal molecular signaling events. The workflow described here enables significant gene and methylation-specific pathways to be identified from whole blood samples of individuals exposed to environmentally relevant chemical levels. As future efforts solidify specific causal relationships between these molecular events and relevant apical endpoints, this novel workflow could aid risk assessments by identifying molecular targets serving as biomarkers of hazard, informing mechanistic understanding, and characterizing dose ranges that promote relevant molecular/epigenetic signaling events occuring in response to chemical exposures.