Epigenome-wide association studies of occupational exposure to benzene and formaldehyde.

Sufficient evidence supports a relationship between certain myeloid neoplasms and exposure to benzene or formaldehyde. DNA methylation could underlie benzene- and formaldehyde-induced health outcomes, but data in exposed human populations are limited. We conducted two cross-sectional epigenome-wide association studies (EWAS), one in workers exposed to benzene and another in workers exposed to formaldehyde. Using HumanMethylation450 BeadChips, we investigated differences in blood cell DNA methylation among 50 benzene-exposed subjects and 48 controls, and among 31 formaldehyde-exposed subjects and 40 controls. We performed CpG-level and regional-level analyses. In the benzene EWAS, we found genome-wide significant alterations, i.e., FWER-controlled P-values <0.05, in the mean and variance of methylation at 22 and 318 CpG sites, respectively, and in mean methylation of a large genomic region. Pathway analysis of genes corresponding to benzene-associated differential methylation sites revealed an impact on the AMPK signalling pathway. In formaldehyde-exposed subjects compared to controls, 9 CpGs in the DUSP22 gene promoter had genome-wide significant decreased methylation variability and a large region of the HOXA5 promoter with 44 CpGs was hypomethylated. Our findings suggest that DNA methylation may contribute to the pathogenesis of diseases related to benzene and formaldehyde exposure. Aberrant expression and methylation of HOXA5 previously has been shown to be clinically significant in myeloid leukaemias. The tumour suppressor gene DUSP22 is a potential biomarker of exposure to formaldehyde, and irregularities have been associated with multiple exposures and diseases.

Machine learning with routine electronic medical record data to identify people at high risk of disengagement from HIV care in Tanzania.

Machine learning methods for health care delivery optimization have the potential to improve retention in HIV care, a critical target of global efforts to end the epidemic. However, these methods have not been widely applied to medical record data in low- and middle-income countries. We used an ensemble decision tree approach to predict risk of disengagement from HIV care (missing an appointment by ≥28 days) in Tanzania. Our approach used routine electronic medical records (EMR) from the time of antiretroviral therapy (ART) initiation through 24 months of follow-up for 178 adults (63% female). We compared prediction accuracy when using EMR-based predictors alone and in combination with sociodemographic survey data collected by a research study. Models that included only EMR-based indicators and incorporated changes across past clinical visits achieved a mean accuracy of 75.2% for predicting risk of disengagement in the next 6 months, with a mean sensitivity of 54.7% for targeting the 30% highest-risk individuals. Additionally including survey-based predictors only modestly improved model performance. The most important variables for prediction were time-varying EMR indicators including changes in treatment status, body weight, and WHO clinical stage. Machine learning methods applied to existing EMR data in resource-constrained settings can predict individuals' future risk of disengagement from HIV care, potentially enabling better targeting and efficiency of interventions to promote retention in care.

Exposure to arsenic at different life-stages and DNA methylation meta-analysis in buccal cells and leukocytes.

BACKGROUND: Arsenic (As) exposure through drinking water is a global public health concern. Epigenetic dysregulation including changes in DNA methylation (DNAm), may be involved in arsenic toxicity. Epigenome-wide association studies (EWAS) of arsenic exposure have been restricted to single populations and comparison across EWAS has been limited by methodological differences. Leveraging data from epidemiological studies conducted in Chile and Bangladesh, we use a harmonized data processing and analysis pipeline and meta-analysis to combine results from four EWAS. METHODS: DNAm was measured among adults in Chile with and without prenatal and early-life As exposure in PBMCs and buccal cells (N = 40, 850K array) and among men in Bangladesh with high and low As exposure in PBMCs (N = 32, 850K array; N = 48, 450K array). Linear models were used to identify differentially methylated positions (DMPs) and differentially variable positions (DVPs) adjusting for age, smoking, cell type, and sex in the Chile cohort. Probes common across EWAS were meta-analyzed using METAL, and differentially methylated and variable regions (DMRs and DVRs, respectively) were identified using comb-p. KEGG pathway analysis was used to understand biological functions of DMPs and DVPs. RESULTS: In a meta-analysis restricted to PBMCs, we identified one DMP and 23 DVPs associated with arsenic exposure; including buccal cells, we identified 3 DMPs and 19 DVPs (FDR < 0.05). Using meta-analyzed results, we identified 11 DMRs and 11 DVRs in PBMC samples, and 16 DMRs and 19 DVRs in PBMC and buccal cell samples. One region annotated to LRRC27 was identified as a DMR and DVR. Arsenic-associated KEGG pathways included lysosome, autophagy, and mTOR signaling, AMPK signaling, and one carbon pool by folate. CONCLUSIONS: Using a two-step process of (1) harmonized data processing and analysis and (2) meta-analysis, we leverage four DNAm datasets from two continents of individuals exposed to high levels of As prenatally and during adulthood to identify DMPs and DVPs associated with arsenic exposure. Our approach suggests that standardizing analytical pipelines can aid in identifying biological meaningful signals.

Applying genome-wide CRISPR to identify known and novel genes and pathways that modulate formaldehyde toxicity.

Formaldehyde (FA), a ubiquitous environmental pollutant, is classified as a Group I human carcinogen by the International Agency for Research on Cancer. Previously, we reported that FA induced hematotoxicity and chromosomal aneuploidy in exposed workers and toxicity in bone marrow and hematopoietic stem cells of experimental animals. Using functional toxicogenomic profiling in yeast, we identified genes and cellular processes modulating eukaryotic FA cytotoxicity. Although we validated some of these findings in yeast, many specific genes, pathways and mechanisms of action of FA in human cells are not known. In the current study, we applied genome-wide, loss-of-function CRISPR screening to identify modulators of FA toxicity in the human hematopoietic K562 cell line. We assessed the cellular genetic determinants of susceptibility and resistance to FA at 40, 100 and 150 µM (IC10, IC20 and IC60, respectively) at two time points, day 8 and day 20. We identified multiple candidate genes that increase sensitivity (e.g. ADH5, ESD and FANC family) or resistance (e.g. FASN and KDM6A) to FA when disrupted. Pathway analysis revealed a major role for the FA metabolism and Fanconi anemia pathway in FA tolerance, consistent with findings from previous studies. Additional network analyses revealed potential new roles for one-carbon metabolism, fatty acid synthesis and mTOR signaling in modulating FA toxicity. Validation of these novel findings will further enhance our understanding of FA toxicity in human cells. Our findings support the utility of CRISPR-based functional genomics screening of environmental chemicals.

Ripples on spikes show increased phase-amplitude coupling in mesial temporal lobe epilepsy seizure-onset zones.

OBJECTIVE: Ripples (80-150 Hz) recorded from clinical macroelectrodes have been shown to be an accurate biomarker of epileptogenic brain tissue. We investigated coupling between epileptiform spike phase and ripple amplitude to better understand the mechanisms that generate this type of pathologic ripple (pRipple) event. METHODS: We quantified phase amplitude coupling (PAC) between epileptiform electroencephalography (EEG) spike phase and ripple amplitude recorded from intracranial depth macroelectrodes during episodes of sleep in 12 patients with mesial temporal lobe epilepsy. PAC was determined by (1) a phasor transform that corresponds to the strength and rate of ripples coupled with spikes, and a (2) ripple-triggered average to measure the strength, morphology, and spectral frequency of the modulating and modulated signals. Coupling strength was evaluated in relation to recording sites within and outside the seizure-onset zone (SOZ). RESULTS: Both the phasor transform and ripple-triggered averaging methods showed that ripple amplitude was often robustly coupled with epileptiform EEG spike phase. Coupling was found more regularly inside than outside the SOZ, and coupling strength correlated with the likelihood a macroelectrode's location was within the SOZ (p < 0.01). The ratio of the rate of ripples coupled with EEG spikes inside the SOZ to rates of coupled ripples in non-SOZ was greater than the ratio of rates of ripples on spikes detected irrespective of coupling (p < 0.05). Coupling strength correlated with an increase in mean normalized ripple amplitude (p < 0.01), and a decrease in mean ripple spectral frequency (p < 0.05). SIGNIFICANCE: Generation of low-frequency (80-150 Hz) pRipples in the SOZ involves coupling between epileptiform spike phase and ripple amplitude. The changes in excitability reflected as epileptiform spikes may also cause clusters of pathologically interconnected bursting neurons to grow and synchronize into aberrantly large neuronal assemblies.