Person:

Christiani, David

GWAS has facilitated greatly the discovery of risk SNPs associated with complex diseases. Traditional methods analyze SNP individually and are limited by low power and reproducibility since correction for multiple comparisons is necessary. Several methods have been proposed based on grouping SNPs into SNP sets using biological knowledge and/or genomic features. In this article, we compare the linear kernel machine based test (LKM) and principal components analysis based approach (PCA) using simulated datasets under the scenarios of 0 to 3 causal SNPs, as well as simple and complex linkage disequilibrium (LD) structures of the simulated regions. Our simulation study demonstrates that both LKM and PCA can control the type I error at the significance level of 0.05. If the causal SNP is in strong LD with the genotyped SNPs, both the PCA with a small number of principal components (PCs) and the LKM with kernel of linear or identical-by-state function are valid tests. However, if the LD structure is complex, such as several LD blocks in the SNP set, or when the causal SNP is not in the LD block in which most of the genotyped SNPs reside, more PCs should be included to capture the information of the causal SNP. Simulation studies also demonstrate the ability of LKM and PCA to combine information from multiple causal SNPs and to provide increased power over individual SNP analysis. We also apply LKM and PCA to analyze two SNP sets extracted from an actual GWAS dataset on non-small cell lung cancer.

Background: In occupational settings, boilermakers are exposed to high levels of metallic fine particulate matter (PM2.5) generated during the welding process. The effect of welding PM2.5 on heart rate variability (HRV) has been described, but the relationship between PM2.5, DNA methylation, and HRV is not known. Methods: In this repeated-measures panel study, we recorded resting HRV and measured DNA methylation levels in transposable elements Alu and long interspersed nuclear element-1 (LINE-1) in peripheral blood leukocytes under ambient conditions (pre-shift) and right after a welding task (post-shift) among 66 welders. We also monitored personal PM2.5 level in the ambient environment and during the welding procedure. Results: The concentration of welding PM2.5 was significantly higher than background levels in the union hall (0.43 mg/m3 vs. 0.11 mg/m3, p < 0.0001). The natural log of transformed power in the high frequency range (ln HF) had a significantly negative association with PM2.5 exposure (β = -0.76, p = 0.035). pNN10 and pNN20 also had a negative association with PM2.5 exposure (β = -0.16%, p = 0.006 and β = -0.13%, p = 0.030, respectively). PM2.5 was positively associated with LINE-1 methylation [β = 0.79%, 5-methylcytosince (%mC), p = 0.013]; adjusted for covariates. LINE-1 methylation did not show an independent association with HRV. Conclusions: Acute decline of HRV was observed following exposure to welding PM2.5 and evidence for an epigenetic response of transposable elements to short-term exposure to high-level metal-rich particulates was reported. Electronic supplementary material The online version of this article (doi:10.1186/1471-2458-14-1279) contains supplementary material, which is available to authorized users.

Background: Previous studies in occupational exposure and lung function have focused only on the main effect of occupational exposure or genetics on lung function. Some disease-susceptible genes may be missed due to their low marginal effects, despite potential involvement in the disease process through interactions with the environment. Through comprehensive genome-wide gene-environment interaction studies, we can uncover these susceptibility genes. Our objective in this study was to explore gene by occupational exposure interaction effects on lung function using both the individual SNPs approach and the genetic network approach. Methods: The study population comprised the Offspring Cohort and the Third Generation from the Framingham Heart Study. We used forced expiratory volume in one second (FEV1) and ratio of FEV1 to forced vital capacity (FVC) as outcomes. Occupational exposures were classified using a population-specific job exposure matrix. We performed genome-wide gene-environment interaction analysis, using the Affymetrix 550 K mapping array for genotyping. A linear regression-based generalized estimating equation was applied to account for within-family relatedness. Network analysis was conducted using results from single-nucleotide polymorphism (SNP)-level analyses and from gene expression study results. Results: There were 4,785 participants in total. SNP-level analysis and network analysis identified SNP rs9931086 (Pinteraction =1.16 × 10-7) in gene SLC38A8, which may significantly modify the effects of occupational exposure on FEV1. Genes identified from the network analysis included CTLA-4, HDAC, and PPAR-alpha. Conclusions: Our study implies that SNP rs9931086 in SLC38A8 and genes CTLA-4, HDAC, and PPAR-alpha, which are related to inflammatory processes, may modify the effect of occupational exposure on lung function.

Background: Chronic inflammation from recurring trauma is an underlying pathophysiological basis of numerous diseases. Furthermore, it may result in cell death, scarring, fibrosis, and loss of tissue function. In states of inflammation, subsequent increases in oxidative stress and cellular division may lead to the accelerated erosion of telomeres, crucial genomic structures which protect chromosomes from decay. However, the association between plasma inflammatory marker concentrations and telomere length has been inconsistent in previous studies. Objective: The purpose of this study was to determine the longitudinal association between telomere length and plasma inflammatory biomarker concentrations including: CRP, SAA, sICAM-1, sVCAM-1, VEGF, TNF-α, IL-1β, IL-2, IL-6, IL-8, and IL-10. Methods: The longitudinal study population consisted of 87 subjects. The follow-up period was approximately 2 years. Plasma inflammatory biomarker concentrations were assessed using highly sensitive electrochemiluminescent assays. Leukocyte relative telomere length was assessed using Real-Time qPCR. Linear mixed effects regression models were used to analyze the association between repeated-measurements of relative telomere length as the outcome and each inflammatory biomarker concentration as continuous exposures separately. The analyses controlled for major potential confounders and white blood cell differentials. Results: At any follow-up time, each incremental ng/mL increase in plasma CRP concentration was associated with a decrease in telomere length of −2.6×10−2 (95%CI: −4.3×10−2, −8.2×10−3, p = 0.004) units. Similarly, the estimate for the negative linear association between SAA and telomere length was −2.6×10−2 (95%CI:−4.5×10−2, −6.1×10−3, p = 0.011). No statistically significant associations were observed between telomere length and plasma concentrations of pro-inflammatory interleukins, TNF-α, and VEGF. Conclusions: Findings from this study suggest that increased systemic inflammation, consistent with vascular injury, is associated with decreased leukocyte telomere length.

Background: Chronic exposure to arsenic in drinking water is associated with increased risk of type 2 diabetes mellitus (T2DM) but the underlying molecular mechanism remains unclear. Objectives: This study evaluated the interaction between single nucleotide polymorphisms (SNPs) in genes associated with diabetes and arsenic exposure in drinking water on the risk of developing T2DM. Methods: In 2009–2011, we conducted a follow up study of 957 Bangladeshi adults who participated in a case-control study of arsenic-induced skin lesions in 2001–2003. Logistic regression models were used to evaluate the association between 38 SNPs in 18 genes and risk of T2DM measured at follow up. T2DM was defined as having a blood hemoglobin A1C level greater than or equal to 6.5% at follow-up. Arsenic exposure was characterized by drinking water samples collected from participants' tubewells. False discovery rates were applied in the analysis to control for multiple comparisons. Results: Median arsenic levels in 2001–2003 were higher among diabetic participants compared with non-diabetic ones (71.6 µg/L vs. 12.5 µg/L, p-value <0.001). Three SNPs in ADAMTS9 were nominally associated with increased risk of T2DM (rs17070905, Odds Ratio (OR) = 2.30, 95% confidence interval (CI) 1.17–4.50; rs17070967, OR = 2.02, 95%CI 1.00–4.06; rs6766801, OR = 2.33, 95%CI 1.18–4.60), but these associations did not reach the statistical significance after adjusting for multiple comparisons. A significant interaction between arsenic and NOTCH2 (rs699780) was observed which significantly increased the risk of T2DM (p for interaction = 0.003; q-value = 0.021). Further restricted analysis among participants exposed to water arsenic of less than 148 µg/L showed consistent results for interaction between the NOTCH2 variant and arsenic exposure on T2DM (p for interaction = 0.048; q-value = 0.004). Conclusions: These findings suggest that genetic variation in NOTCH2 increased susceptibility to T2DM among people exposed to inorganic arsenic. Additionally, genetic variants in ADAMTS9 may increase the risk of T2DM.

Objective: Acceleration (AC) and deceleration (DC) capacities measure heart rate variability during speeding up and slowing down of the heart, respectively. We investigated associations between AC and DC with occupational short-term metal PM2.5 exposures. Methods: A panel of 48 male welders had particulate matter less than 2.5 microns in diameter (PM2.5) exposure measurements over 4–6 h repeated over 5 sampling periods between January 2010 and June 2012. We simultaneously obtained continuous recordings of digital ECG using a Holter monitor. We analysed ECG data in the time domain to obtain hourly AC and DC. Linear mixed models were used to assess the associations between hourly PM2.5 exposure and each of hourly AC and DC, controlling for age, smoking status, active smoking, exposure to secondhand smoke, season/time of day when ECG reading was obtained and baseline AC or DC. We also ran lagged exposure response models for each successive hour up to 3 h after onset of exposure. Results: Mean (SD) shift PM2.5 exposure during welding was 0.47 (0.43) mg/m3. Significant exposure–response associations were found for AC and DC with increased PM2.5 exposure. In our adjusted models without any lag between exposure and response, a 1 mg/m3 increase of PM2.5 was associated with a decrease of 1.46 (95% CI 1.00 to 1.92) ms in AC and a decrease of 1.00 (95% CI 0.53 to 1.46) ms in DC. The effect of PM2.5 on AC and DC was maximal immediately postexposure and lasted 1 h following exposure. Conclusions: There are short-term effects of metal particulates on AC and DC.

ABSTRACT DNA methylation is one of the potential epigenetic mechanisms associated with various adverse cardiovascular effects; however, its association with cardiac autonomic dysfunction, in particular, is unknown. In the current study, we aimed to identify epigenetic variants associated with alterations in cardiac autonomic responses. Cardiac autonomic responses were measured with two novel markers: acceleration capacity (AC) and deceleration capacity (DC). We examined DNA methylation levels at more than 472,506 CpG probes through the Illumina Infinium HumanMethylation450 BeadChip assay. We conducted separate linear mixed models to examine associations of DNA methylation levels at each CpG with AC and DC. One CpG (cg26829071) located in the GPR133 gene was negatively associated with DC values after multiple testing corrections through false discovery rate. Our study suggests the potential functional importance of methylation in cardiac autonomic responses. Findings from the current study need to be replicated in future studies in a larger population.

Background: Environmental and occupational exposure to metals is ubiquitous worldwide, and understanding the hazardous metal components in this complex mixture is essential for environmental and occupational regulations. Objective: To identify hazardous components from metal mixtures that are associated with alterations in cardiac autonomic responses. Methods: Urinary concentrations of 16 types of metals were examined and ‘acceleration capacity’ (AC) and ‘deceleration capacity’ (DC), indicators of cardiac autonomic effects, were quantified from ECG recordings among 54 welders. We fitted linear mixed-effects models with least absolute shrinkage and selection operator (LASSO) to identify metal components that are associated with AC and DC. The Bayesian Information Criterion was used as the criterion for model selection procedures. Results: Mercury and chromium were selected for DC analysis, whereas mercury, chromium and manganese were selected for AC analysis through the LASSO approach. When we fitted the linear mixed-effects models with ‘selected’ metal components only, the effect of mercury remained significant. Every 1 µg/L increase in urinary mercury was associated with −0.58 ms (−1.03, –0.13) changes in DC and 0.67 ms (0.25, 1.10) changes in AC. Conclusion: Our study suggests that exposure to several metals is associated with impaired cardiac autonomic functions. Our findings should be replicated in future studies with larger sample sizes.

Recent technological advancements have permitted high-throughput measurement of the human genome, epigenome, metabolome, transcriptome, and proteome at the population level. We hypothesized that subsets of genes identified from omic studies might have closely related biological functions and thus might interact directly at the network level. Therefore, we conducted an integrative analysis of multi-omic datasets of non-small cell lung cancer (NSCLC) to search for association patterns beyond the genome and transcriptome. A large, complex, and robust gene network containing well-known lung cancer-related genes, including EGFR and TERT, was identified from combined gene lists for lung adenocarcinoma. Members of the hypoxia-inducible factor (HIF) gene family were at the center of this network. Subsequent sequencing of network hub genes within a subset of samples from the Transdisciplinary Research in Cancer of the Lung-International Lung Cancer Consortium (TRICL-ILCCO) consortium revealed a SNP (rs12614710) in EPAS1 associated with NSCLC that reached genome-wide significance (OR = 1.50; 95% CI: 1.31–1.72; p = 7.75 × 10−9). Using imputed data, we found that this SNP remained significant in the entire TRICL-ILCCO consortium (p = .03). Additional functional studies are warranted to better understand interrelationships among genetic polymorphisms, DNA methylation status, and EPAS1 expression.

Background: Neurodevelopment is a complex process involving both genetic and environmental factors. Prenatal exposure to lead (Pb) has been associated with lower performance on neurodevelopmental tests. Adverse neurodevelopmental outcomes are more frequent and/or more severe when toxic exposures interact with genetic susceptibility. Methods: To explore possible loci associated with increased susceptibility to prenatal Pb exposure, we performed a genome-wide gene-environment interaction study (GWIS) in young children from Mexico (n = 390) and Bangladesh (n = 497). Prenatal Pb exposure was estimated by cord blood Pb concentration. Neurodevelopment was assessed using the Bayley Scales of Infant Development. Results: We identified a locus on chromosome 8, containing UNC5D, and demonstrated evidence of its genome-wide significance with mental composite scores (rs9642758, p meta = 4.35 × 10−6). Within this locus, the joint effects of two independent single nucleotide polymorphisms (SNPs, rs9642758 and rs10503970) had a p-value of 4.38 × 10−9 for mental composite scores. Correlating GWIS results with in vitro transcriptomic profiles identified one common gene, SLC1A5, which is involved in synaptic function, neuronal development, and excitotoxicity. Further analysis revealed interconnected interactions that formed a large network of 52 genes enriched with oxidative stress genes and neurodevelopmental genes. Conclusions: Our findings suggest that certain genetic polymorphisms within/near genes relevant to neurodevelopment might modify the toxic effects of Pb exposure via oxidative stress. Electronic supplementary material The online version of this article (doi:10.1186/s12940-017-0288-3) contains supplementary material, which is available to authorized users.