Impact of tobacco, alcohol, and marijuana on genome-wide DNA methylation and its relationship with hypertension

Tobacco, alcohol, and marijuana consumption is an important public health problem because of their high use worldwide and their association with the risk of mortality and many health conditions, such as hypertension, which is the commonest risk factor for death throughout the world. A likely pathway of action of substance consumption leading to persistent hypertension is DNA methylation. Here, we evaluated the effects of tobacco, alcohol, and marijuana on DNA methylation in the same cohort (N = 3,424). Three epigenome-wide association studies (EWAS) were assessed in whole blood using the InfiniumHumanMethylationEPIC BeadChip. We also evaluated the mediation of the top CpG sites in the association between substance consumption and hypertension. Our analyses showed 2,569 CpG sites differentially methylated by alcohol drinking and 528 by tobacco smoking. We did not find significant associations with marijuana consumption after correcting for multiple comparisons. We found 61 genes overlapping between alcohol and tobacco that were enriched in biological processes involved in the nervous and cardiovascular systems. In the mediation analysis, we found 66 CpG sites that significantly mediated the effect of alcohol consumption on hypertension. The top alcohol-related CpG site (cg06690548, P-value = 5.9·10^-83) mapped to SLC7A11 strongly mediated 70.5% of the effect of alcohol consumption on hypertension (P-value = 0.006). Our findings suggest that DNA methylation should be considered for new targets in hypertension prevention and management, particularly concerning alcohol consumption. Our data also encourage further research into the use of methylation in blood to study the neurological and cardiovascular effects of substance consumption.

State-of-the-art methods for exposure-health studies: Results from the exposome data challenge event

The exposome recognizes that individuals are exposed simultaneously to a multitude of different environmental factors and takes a holistic approach to the discovery of etiological factors for disease. However, challenges arise when trying to quantify the health effects of complex exposure mixtures. Analytical challenges include dealing with high dimensionality, studying the combined effects of these exposures and their interactions, integrating causal pathways, and integrating high-throughput omics layers. To tackle these challenges, the Barcelona Institute for Global Health (ISGlobal) held a data challenge event open to researchers from all over the world and from all expertises. Analysts had a chance to compete and apply state-of-the-art methods on a common partially simulated exposome dataset (based on real case data from the HELIX project) with multiple correlated exposure variables (P > 100 exposure variables) arising from general and personal environments at different time points, biological molecular data (multi-omics: DNA methylation, gene expression, proteins, metabolomics) and multiple clinical phenotypes in 1301 mother-child pairs. Most of the methods presented included feature selection or feature reduction to deal with the high dimensionality of the exposome dataset. Several approaches explicitly searched for combined effects of exposures and/or their interactions using linear index models or response surface methods, including Bayesian methods. Other methods dealt with the multi-omics dataset in mediation analyses using multiple-step approaches. Here we discuss features of the statistical models used and provide the data and codes used, so that analysts have examples of implementation and can learn how to use these methods. Overall, the exposome data challenge presented a unique opportunity for researchers from different disciplines to create and share state-of-the-art analytical methods, setting a new standard for open science in the exposome and environmental health field.

Software Application Profile: exposomeShiny—a toolbox for exposome data analysis

Motivation

Studying the role of the exposome in human health and its impact on different omic layers requires advanced statistical methods. Many of these methods are implemented in different R and Bioconductor packages, but their use may require strong expertise in R, in writing pipelines and in using new R classes which may not be familiar to non-advanced users. ExposomeShiny provides a bridge between researchers and most of the state-of-the-art exposome analysis methodologies, without the need of advanced programming skills.

Implementation

ExposomeShiny is a standalone web application implemented in R. It is available as source files and can be installed in any server or computer avoiding problems with data confidentiality. It is executed in RStudio which opens a browser window with the web application.

General features

The presented implementation allows the conduct of: (i) data pre-processing: normalization and missing imputation (including limit of detection); (ii) descriptive analysis; (iii) exposome principal component analysis (PCA) and hierarchical clustering; (iv) exposome-wide association studies (ExWAS) and variable selection ExWAS; (v) omic data integration by single association and multi-omic analyses; and (vi) post-exposome data analyses to gain biological insight for the exposures, genes or using the Comparative Toxicogenomics Database (CTD) and pathway analysis.

Availability

The exposomeShiny source code is freely available on Github at [https://github.com/isglobal-brge/exposomeShiny], Git tag v1.4. The software is also available as a Docker image [https://hub.docker.com/r/brgelab/exposome-shiny], tag v1.4. A user guide with information about the analysis methodologies as well as information on how to use exposomeShiny is freely hosted at [https://isglobal-brge.github.io/exposome_bookdown/].

Advancing tools for human early lifecourse exposome research and translation (ATHLETE): Project overview

Early life stages are vulnerable to environmental hazards and present important windows of opportunity for lifelong disease prevention. This makes early life a relevant starting point for exposome studies. The Advancing Tools for Human Early Lifecourse Exposome Research and Translation (ATHLETE) project aims to develop a toolbox of exposome tools and a Europe-wide exposome cohort that will be used to systematically quantify the effects of a wide range of community- and individual-level environmental risk factors on mental, cardiometabolic, and respiratory health outcomes and associated biological pathways, longitudinally from early pregnancy through to adolescence. Exposome tool and data development include as follows: (1) a findable, accessible, interoperable, reusable (FAIR) data infrastructure for early life exposome cohort data, including 16 prospective birth cohorts in 11 European countries; (2) targeted and nontargeted approaches to measure a wide range of environmental exposures (urban, chemical, physical, behavioral, social); (3) advanced statistical and toxicological strategies to analyze complex multidimensional exposome data; (4) estimation of associations between the exposome and early organ development, health trajectories, and biological (metagenomic, metabolomic, epigenetic, aging, and stress) pathways; (5) intervention strategies to improve early life urban and chemical exposomes, co-produced with local communities; and (6) child health impacts and associated costs related to the exposome. Data, tools, and results will be assembled in an openly accessible toolbox, which will provide great opportunities for researchers, policymakers, and other stakeholders, beyond the duration of the project. ATHLETE's results will help to better understand and prevent health damage from environmental exposures and their mixtures from the earliest parts of the life course onward.

Performance of approaches relying on multidimensional intermediary data to decipher causal relationships between the exposome and health: A simulation study under various causal structures

Challenges in the assessment of the health effects of the exposome, defined as encompassing all environmental exposures from the prenatal period onwards, include a possibly high rate of false positive signals. It might be overcome using data dimension reduction techniques. Data from the biological layers lying between the exposome and its possible health consequences, such as the methylome, may help reducing exposome dimension. We aimed to quantify the performances of approaches relying on the incorporation of an intermediary biological layer to relate the exposome and health, and compare them with agnostic approaches ignoring the intermediary layer. We performed a Monte-Carlo simulation, in which we generated realistic exposome and intermediary layer data by sampling with replacement real data from the Helix exposome project. We generated a Gaussian outcome assuming linear relationships between the three data layers, in 2381 scenarios under five different causal structures, including mediation and reverse causality. We tested 3 agnostic methods considering only the exposome and the health outcome: ExWAS (for Exposome-Wide Association study), DSA, LASSO; and 3 methods relying on an intermediary layer: two implementations of our new oriented Meet-in-the-Middle (oMITM) design, using ExWAS and DSA, and a mediation analysis using ExWAS. Methods' performances were assessed through their sensitivity and FDP (False-Discovery Proportion). The oMITM-based methods generally had lower FDP than the other approaches, possibly at a cost in terms of sensitivity; FDP was in particular lower under a structure of reverse causality and in some mediation scenarios. The oMITM-DSA implementation showed better performances than oMITM-ExWAS, especially in terms of FDP. Among the agnostic approaches, DSA showed the highest performance. Integrating information from intermediary biological layers can help lowering FDP in studies of the exposome health effects; in particular, oMITM seems less sensitive to reverse causality than agnostic exposome-health association studies.

Early-life environmental exposure determinants of child behavior in Europe: A longitudinal, population-based study

BACKGROUND

Environmental exposures in early life influence the development of behavioral outcomes in children, but research has not considered multiple exposures. We therefore aimed to investigate the impact of a broad spectrum of pre- and postnatal environmental exposures on child behavior.

METHODS AND FINDINGS

We used data from the HELIX (Human Early Life Exposome) project, which was based on six longitudinal population-based birth cohorts in Europe. At 6-11 years, children underwent a follow-up to characterize their exposures and assess behavioral problems. We measured 88 prenatal and 123 childhood environmental factors, including outdoor, indoor, chemical, lifestyle and social exposures. Parent-reported behavioral problems included (1) internalizing, (2) externalizing scores, using the child behavior checklist (CBCL), and (3) the Conner's Attention Deficit Hyperactivity Disorder (ADHD) index, all outcomes being discrete raw counts. We applied LASSO penalized negative binomial regression models to identify which exposures were associated with the outcomes, while adjusting for co-exposures. In the 1287 children (mean age 8.0 years), 7.3% had a neuropsychiatric medical diagnosis according to parent's reports. During pregnancy, smoking and car traffic showing the strongest associations (e.g. smoking with ADHD index, aMR:1.31 [1.09; 1.59]) among the 13 exposures selected by LASSO, for at least one of the outcomes. During childhood, longer sleep duration, healthy diet and higher family social capital were associated with reduced scores whereas higher exposure to lead, copper, indoor air pollution, unhealthy diet were associated with increased scores. Unexpected decreases in behavioral scores were found with polychlorinated biphenyls (PCBs) and organophosphate (OP) pesticides.

CONCLUSIONS

Our systematic exposome approach identified several environmental contaminants and healthy lifestyle habits that may influence behavioral problems in children. Modifying environmental exposures early in life may limit lifetime mental health risk.

DNA methylation signatures of aggression and closely related constructs: A meta-analysis of epigenome-wide studies across the lifespan

DNA methylation profiles of aggressive behavior may capture lifetime cumulative effects of genetic, stochastic, and environmental influences associated with aggression. Here, we report the first large meta-analysis of epigenome-wide association studies (EWAS) of aggressive behavior (N = 15,324 participants). In peripheral blood samples of 14,434 participants from 18 cohorts with mean ages ranging from 7 to 68 years, 13 methylation sites were significantly associated with aggression (alpha = 1.2 × 10-7; Bonferroni correction). In cord blood samples of 2425 children from five cohorts with aggression assessed at mean ages ranging from 4 to 7 years, 83% of these sites showed the same direction of association with childhood aggression (r = 0.74, p = 0.006) but no epigenome-wide significant sites were found. Top-sites (48 at a false discovery rate of 5% in the peripheral blood meta-analysis or in a combined meta-analysis of peripheral blood and cord blood) have been associated with chemical exposures, smoking, cognition, metabolic traits, and genetic variation (mQTLs). Three genes whose expression levels were associated with top-sites were previously linked to schizophrenia and general risk tolerance. At six CpGs, DNA methylation variation in blood mirrors variation in the brain. On average 44% (range = 3-82%) of the aggression-methylation association was explained by current and former smoking and BMI. These findings point at loci that are sensitive to chemical exposures with potential implications for neuronal functions. We hope these results to be a starting point for studies leading to applications as peripheral biomarkers and to reveal causal relationships with aggression and related traits.

DNA methylation as the link between migration and the major noncommunicable diseases: the RODAM study

Aim: We assessed epigenome-wide DNA methylation (DNAm) differences between migrant and non-migrant Ghanaians. Materials & methods: We used the Illumina Infinium® HumanMethylation450 BeadChip to profile DNAm of 712 Ghanaians in whole blood. We used linear models to detect differentially methylated positions (DMPs) associated with migration. We performed multiple post hoc analyses to validate our findings. Results: We identified 13 DMPs associated with migration (delta-beta values: 0.2-4.5%). Seven DMPs in CPLX2, EIF4E3, MEF2D, TLX3, ST8SIA1, ANG and CHRM3 were independent of extrinsic genomic influences in public databases. Two DMPs in NLRC5 were associated with duration of stay in Europe among migrants. All DMPs were biologically linked to migration-related factors. Conclusion: Our findings provide the first insights into DNAm differences between migrants and non-migrants.

Association between the pregnancy exposome and fetal growth

BACKGROUND

Several environmental contaminants were shown to possibly influence fetal growth, generally from single exposure family studies, which are prone to publication bias and confounding by co-exposures. The exposome paradigm offers perspectives to avoid selective reporting of findings and to control for confounding by co-exposures. We aimed to characterize associations of fetal growth with the pregnancy chemical and external exposomes.

METHODS

Within the Human Early-Life Exposome project, 131 prenatal exposures were assessed using biomarkers and environmental models in 1287 mother-child pairs from six European cohorts. We investigated their associations with fetal growth using a deletion-substitution-addition (DSA) algorithm considering all exposures simultaneously, and an exposome-wide association study (ExWAS) considering each exposure independently. We corrected for exposure measurement error and tested for exposure-exposure and sex-exposure interactions.

RESULTS

The DSA model identified lead blood level, which was associated with a 97 g birth weight decrease for each doubling in lead concentration. No exposure passed the multiple testing-corrected significance threshold of ExWAS; without multiple testing correction, this model was in favour of negative associations of lead, fine particulate matter concentration and absorbance with birth weight, and of a positive sex-specific association of parabens with birth weight in boys. No two-way interaction between exposure variables was identified.

CONCLUSIONS

This first large-scale exposome study of fetal growth simultaneously considered >100 environmental exposures. Compared with single exposure studies, our approach allowed making all tests (usually reported in successive publications) explicit. Lead exposure is still a health concern in Europe and parabens health effects warrant further investigation.

Genomic analysis of diet composition finds novel loci and associations with health and lifestyle

We conducted genome-wide association studies (GWAS) of relative intake from the macronutrients fat, protein, carbohydrates, and sugar in over 235,000 individuals of European ancestries. We identified 21 unique, approximately independent lead SNPs. Fourteen lead SNPs are uniquely associated with one macronutrient at genome-wide significance (P < 5 × 10^-8), while five of the 21 lead SNPs reach suggestive significance (P < 1 × 10^-5) for at least one other macronutrient. While the phenotypes are genetically correlated, each phenotype carries a partially unique genetic architecture. Relative protein intake exhibits the strongest relationships with poor health, including positive genetic associations with obesity, type 2 diabetes, and heart disease (r_g ≈ 0.15-0.5). In contrast, relative carbohydrate and sugar intake have negative genetic correlations with waist circumference, waist-hip ratio, and neighborhood deprivation (|r_g| ≈ 0.1-0.3) and positive genetic correlations with physical activity (r_g ≈ 0.1 and 0.2). Relative fat intake has no consistent pattern of genetic correlations with poor health but has a negative genetic correlation with educational attainment (r_g ≈-0.1). Although our analyses do not allow us to draw causal conclusions, we find no evidence of negative health consequences associated with relative carbohydrate, sugar, or fat intake. However, our results are consistent with the hypothesis that relative protein intake plays a role in the etiology of metabolic dysfunction.

Comprehensive study of the exposome and omic data using rexposome Bioconductor Packages

SUMMARY

Genomics has dramatically improved our understanding of the molecular origins of certain human diseases. Nonetheless, our health is also influenced by the cumulative impact of exposures experienced across the life course (termed 'exposome'). The study of the high-dimensional exposome offers a new paradigm for investigating environmental contributions to disease etiology. However, there is a lack of bioinformatics tools for managing, visualizing and analyzing the exposome. The analysis data should include both association with health outcomes and integration with omic layers. We provide a generic framework called rexposome project, developed in the R/Bioconductor architecture that includes object-oriented classes and methods to leverage high-dimensional exposome data in disease association studies including its integration with a variety of high-throughput data types. The usefulness of the package is illustrated by analyzing a real dataset including exposome data, three health outcomes related to respiratory diseases and its integration with the transcriptome and methylome.

AVAILABILITY AND IMPLEMENTATION

rexposome project is available at https://isglobal-brge.github.io/rexposome/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Cadmium, Smoking, and Human Blood DNA Methylation Profiles in Adults from the Strong Heart Study

BACKGROUND

The epigenetic effects of individual environmental toxicants in tobacco remain largely unexplored. Cadmium (Cd) has been associated with smoking-related health effects, and its concentration in tobacco smoke is higher in comparison with other metals.

OBJECTIVES

We studied the association of Cd and smoking exposures with human blood DNA methylation (DNAm) profiles. We also evaluated the implication of findings to relevant methylation pathways and the potential contribution of Cd exposure from smoking to explain the association between smoking and site-specific DNAm.

METHODS

We conducted an epigenome-wide association study of urine Cd and self-reported smoking (current and former vs. never, and cumulative smoking dose) with blood DNAm in 790,026 CpGs (methylation sites) measured with the Illumina Infinium Human MethylationEPIC (Illumina Inc.) platform in 2,325 adults 45-74 years of age who participated in the Strong Heart Study in 1989-1991. In a mediation analysis, we estimated the amount of change in DNAm associated with smoking that can be independently attributed to increases in urine Cd concentrations from smoking. We also conducted enrichment analyses and in silico protein-protein interaction networks to explore the biological relevance of the findings.

RESULTS

At a false discovery rate (FDR)-corrected level of 0.05, we found 6 differentially methylated positions (DMPs) for Cd; 288 and 17, respectively, for current and former smoking status; and 77 for cigarette pack-years. Enrichment analyses of these DMPs displayed enrichment of 58 and 6 Gene Ontology and Kyoto Encyclopedia of Genes and Genomes gene sets, respectively, including biological pathways for cancer and cardiovascular disease. In in silico protein-to-protein networks, we observed key proteins in DNAm pathways directly and indirectly connected to Cd- and smoking-DMPs. Among DMPs that were significant for both Cd and current smoking (annotated to PRSS23, AHRR, F2RL3, RARA, and 2q37.1), we found statistically significant contributions of Cd to smoking-related DNAm.

CONCLUSIONS

Beyond replicating well-known smoking epigenetic signatures, we found novel DMPs related to smoking. Moreover, increases in smoking-related Cd exposure were associated with differential DNAm. Our integrative analysis supports a biological link for Cd and smoking-associated health effects, including the possibility that Cd is partly responsible for smoking toxicity through epigenetic changes. https://doi.org/10.1289/EHP6345.

Using methylome data to inform exposome-health association studies: An application to the identification of environmental drivers of child body mass index

BACKGROUND

The exposome is defined as encompassing all environmental exposures one undergoes from conception onwards. Challenges of the application of this concept to environmental-health association studies include a possibly high false-positive rate.

OBJECTIVES

We aimed to reduce the dimension of the exposome using information from DNA methylation as a way to more efficiently characterize the relation between exposome and child body mass index (BMI).

METHODS

Among 1,173 mother-child pairs from HELIX cohort, 216 exposures ("whole exposome") were characterized. BMI and DNA methylation from immune cells of peripheral blood were assessed in children at age 6-10 years. A priori reduction of the methylome to preselect BMI-relevant CpGs was performed using biological pathways. We then implemented a tailored Meet-in-the-Middle approach to identify from these CpGs candidate mediators in the exposome-BMI association, using univariate linear regression models corrected for multiple testing: this allowed to point out exposures most likely to be associated with BMI ("reduced exposome"). Associations of this reduced exposome with BMI were finally tested. The approach was compared to an agnostic exposome-wide association study (ExWAS) ignoring the methylome.

RESULTS

Among the 2284 preselected CpGs (0.6% of the assessed CpGs), 62 were associated with BMI. Four factors (3 postnatal and 1 prenatal) of the exposome were associated with at least one of these CpGs, among which postnatal blood level of copper and PFOS were directly associated with BMI, with respectively positive and negative estimated effects. The agnostic ExWAS identified 18 additional postnatal exposures, including many persistent pollutants, generally unexpectedly associated with decreased BMI.

DISCUSSION

Our approach incorporating a priori information identified fewer significant associations than an agnostic approach. We hypothesize that this smaller number corresponds to a higher specificity (and possibly lower sensitivity), compared to the agnostic approach. Indeed, the latter cannot distinguish causal relations from reverse causation, e.g. for persistent compounds stored in fat, whose circulating level is influenced by BMI.

Association of Fish Consumption and Mercury Exposure During Pregnancy With Metabolic Health and Inflammatory Biomarkers in Children

Importance

The balance of mercury risk and nutritional benefit from fish intake during pregnancy for the metabolic health of offspring to date is unknown.

Objective

To assess the associations of fish intake and mercury exposure during pregnancy with metabolic syndrome in children and alterations in biomarkers of inflammation in children.

Design, Setting, and Participants

This population-based prospective birth cohort study used data from studies performed in 5 European countries (France, Greece, Norway, Spain, and the UK) between April 1, 2003, and February 26, 2016, as part of the Human Early Life Exposome (HELIX) project. Mothers and their singleton offspring were followed up until the children were aged 6 to 12 years. Data were analyzed between March 1 and August 2, 2019.

Exposures

Maternal fish intake during pregnancy (measured in times per week) was assessed using validated food frequency questionnaires, and maternal mercury concentration (measured in micrograms per liter) was assessed using maternal whole blood and cord blood samples.

Main Outcomes and Measures

An aggregate metabolic syndrome score for children was calculated using the z scores of waist circumference, systolic and diastolic blood pressures, and levels of triglyceride, high-density lipoprotein cholesterol, and insulin. A higher metabolic syndrome score (score range, -4.9 to 7.5) indicated a poorer metabolic profile. Three protein panels were used to measure several cytokines and adipokines in the plasma of children.

Results

The study included 805 mothers and their singleton children. Among mothers, the mean (SD) age at cohort inclusion or delivery of their infant was 31.3 (4.6) years. A total of 400 women (49.7%) had a high educational level, and 432 women (53.7%) were multiparous. Among children, the mean (SD) age was 8.4 (1.5) years (age range, 6-12 years). A total of 453 children (56.3%) were boys, and 734 children (91.2%) were of white race/ethnicity. Fish intake consistent with health recommendations (1 to 3 times per week) during pregnancy was associated with a 1-U decrease in metabolic syndrome score in children (β = -0.96; 95% CI, -1.49 to -0.42) compared with low fish consumption (<1 time per week) after adjusting for maternal mercury levels and other covariates. No further benefit was observed with fish intake of more than 3 times per week. A higher maternal mercury concentration was independently associated with an increase in the metabolic syndrome score of their offspring (β per 2-fold increase in mercury concentration = 0.18; 95% CI, 0.01-0.34). Compared with low fish intake, moderate and high fish intake during pregnancy were associated with reduced levels of proinflammatory cytokines and adipokines in children. An integrated analysis identified a cluster of children with increased susceptibility to metabolic disease, which was characterized by low fish consumption during pregnancy, high maternal mercury levels, decreased levels of adiponectin in children, and increased levels of leptin, tumor necrosis factor α, and the cytokines interleukin 6 and interleukin 1β in children.

Conclusions and Relevance

Results of this study suggest that moderate fish intake consistent with current health recommendations during pregnancy was associated with improvements in the metabolic health of children, while high maternal mercury exposure was associated with an unfavorable metabolic profile in children.

CTDquerier: a bioconductor R package for Comparative Toxicogenomics DatabaseTM data extraction, visualization and enrichment of environmental and toxicological studies

Summary

Biomedical studies currently include a large volume of genomic and environmental factors for studying the etiology of human diseases. R/Bioconductor projects provide several tools for performing enrichment analysis at gene-pathway level, allowing researchers to develop novel hypotheses. However, there is a need to perform similar analyses at the chemicals-genes or chemicals-diseases levels to provide complementary knowledge of the causal path between chemicals and diseases. While the Comparative Toxicogenomics Database^TM (CTD) provides information about these relationships, there is no software for integrating it into R/Bioconductor analysis pipelines. CTDquerier helps users to easily download CTD data and integrate it in the R/Bioconductor framework. The package also contains functions for visualizing CTD data and performing enrichment analyses. We illustrate how to use the package with a real data analysis of asthma-related genes. CTDquerier is a flexible and easy-to-use Bioconductor package that provides novel hypothesis about the relationships between chemicals and diseases.

Availability and implementation

CTDquerier R package is available through Bioconductor and its development version at https://github.com/isglobal-brge/CTDquerier.

Supplementary information

Supplementary data are available at Bioinformatics online.

Genome-wide association analyses of risk tolerance and risky behaviors in over 1 million individuals identify hundreds of loci and shared genetic influences

Humans vary substantially in their willingness to take risks. In a combined sample of over 1 million individuals, we conducted genome-wide association studies (GWAS) of general risk tolerance, adventurousness, and risky behaviors in the driving, drinking, smoking, and sexual domains. Across all GWAS, we identified hundreds of associated loci, including 99 loci associated with general risk tolerance. We report evidence of substantial shared genetic influences across risk tolerance and the risky behaviors: 46 of the 99 general risk tolerance loci contain a lead SNP for at least one of our other GWAS, and general risk tolerance is genetically correlated ([Formula: see text] ~ 0.25 to 0.50) with a range of risky behaviors. Bioinformatics analyses imply that genes near SNPs associated with general risk tolerance are highly expressed in brain tissues and point to a role for glutamatergic and GABAergic neurotransmission. We found no evidence of enrichment for genes previously hypothesized to relate to risk tolerance.

A measure of agreement across numerous conditions: assessing when changes in network structures are tissue-specific

BACKGROUND

There is great interest to study how gene pathways change their structure across different tissues. The assessment of inter-study reliability of pathway changes across tissues can inform on the fraction of tissues with specific functional changes in network structure. However, there is a lack of agreement measures among studies that independently observe how a group of observations change across conditions. We, therefore, propose λ, a new inter-study reliability measure that determines the consistency to distinguish observations by condition.

RESULTS

We derived λ's distributional characteristics, determine its reliability properties and compared it with Cohen's κ. We studied the co-expression structure of 287 gene pathways across four brain regions in two transcriptomic studies and applied λ to assess the inter-study reliability of the pathways' brain-regional changes. Brain-related pathways showed highest λ; the top value was for the nicotine addiction pathway whose structure was reliably distinguishable among regions with dopaminergic projections.

CONCLUSION

Our results offer novel substantial evidence that changes in network structure across tissues can be inferred independently of samples, algorithms and experiments (RNA-sequencing or microarrays). Reliability measures, such as λ, can inform on the tissues where changes in a network's structure are likely functional. An R package is available at https://github.com/isglobal-brge/lambda .

Variability of urinary concentrations of non-persistent chemicals in pregnant women and school-aged children

BACKGROUND

Exposome studies are challenged by exposure misclassification for non-persistent chemicals, whose temporal variability contributes to bias in dose-response functions.

OBJECTIVES

We evaluated the variability of urinary concentrations of 24 non-persistent chemicals: 10 phthalate metabolites, 7 phenols, 6 organophosphate (OP) pesticide metabolites, and cotinine, between weeks from different pregnancy trimesters in pregnant women, and between days and between seasons in children.

METHODS

154 pregnant women and 152 children from six European countries were enrolled in 2014-2015. Pregnant women provided three urine samples over a day (morning, midday, and night), for one week in the 2nd and 3rd pregnancy trimesters. Children provided two urines a day (morning and night), over two one-week periods, six months apart. We pooled all samples for a given subject that were collected within a week. In children, we also made four daily pools (combining morning and night voids) during the last four days of the first follow-up week. Pools were analyzed for all 24 metabolites of interest. We calculated intraclass-correlation coefficients (ICC) and estimated the number of pools needed to obtain an ICC above 0.80.

RESULTS

All phthalate metabolites and phenols were detected in >90% of pools whereas certain OP pesticide metabolites and cotinine were detected in <43% of pools. We observed fair (ICC = 0.40-0.59) to good (0.60-0.74) between-day reliability of the pools of two samples in children for all chemicals. Reliability was poor (<0.40) to fair between trimesters in pregnant women and between seasons in children. For most chemicals, three daily pools of two urines each (for weekly exposure windows) and four weekly pools of 15-20 urines each would be necessary to obtain an ICC above 0.80.

CONCLUSIONS

This quantification of the variability of biomarker measurements of many non-persistent chemicals during several time windows shows that for many of these compounds a few dozen samples are required to accurately assess exposure over periods encompassing several trimesters or months.

Human Early Life Exposome (HELIX) study: a European population-based exposome cohort

PURPOSE

Essential to exposome research is the collection of data on many environmental exposures from different domains in the same subjects. The aim of the Human Early Life Exposome (HELIX) study was to measure and describe multiple environmental exposures during early life (pregnancy and childhood) in a prospective cohort and associate these exposures with molecular omics signatures and child health outcomes. Here, we describe recruitment, measurements available and baseline data of the HELIX study populations.

PARTICIPANTS

The HELIX study represents a collaborative project across six established and ongoing longitudinal population-based birth cohort studies in six European countries (France, Greece, Lithuania, Norway, Spain and the UK). HELIX used a multilevel study design with the entire study population totalling 31 472 mother-child pairs, recruited during pregnancy, in the six existing cohorts (first level); a subcohort of 1301 mother-child pairs where biomarkers, omics signatures and child health outcomes were measured at age 6-11 years (second level) and repeat-sampling panel studies with around 150 children and 150 pregnant women aimed at collecting personal exposure data (third level).

FINDINGS TO DATE

Cohort data include urban environment, hazardous substances and lifestyle-related exposures for women during pregnancy and their offspring from birth until 6-11 years. Common, standardised protocols were used to collect biological samples, measure exposure biomarkers and omics signatures and assess child health across the six cohorts. Baseline data of the cohort show substantial variation in health outcomes and determinants between the six countries, for example, in family affluence levels, tobacco smoking, physical activity, dietary habits and prevalence of childhood obesity, asthma, allergies and attention deficit hyperactivity disorder.

FUTURE PLANS

HELIX study results will inform on the early life exposome and its association with molecular omics signatures and child health outcomes. Cohort data are accessible for future research involving researchers external to the project.

Characterization of Plasmodium vivax Proteins in Plasma-Derived Exosomes From Malaria-Infected Liver-Chimeric Humanized Mice

Exosomes are extracellular vesicles of endocytic origin containing molecular signatures implying the cell of origin; thus, they offer a unique opportunity to discover biomarkers of disease. Plasmodium vivax, responsible for more than half of all malaria cases outside Africa, is a major obstacle in the goal of malaria elimination due to the presence of dormant liver stages (hypnozoites), which after the initial infection may reactivate to cause disease. Hypnozoite infection is asymptomatic and there are currently no diagnostic tools to detect their presence. The human liver-chimeric (FRG huHep) mouse is a robust P. vivax infection model for exo-erythrocytic development of liver stages, including hypnozoites. We studied the proteome of plasma-derived exosomes isolated from P. vivax infected FRG huHep mice with the objective of identifying liver-stage expressed parasite proteins indicative of infection. Proteomic analysis of these exosomes showed the presence of 290 and 234 proteins from mouse and human origin, respectively, including canonical exosomal markers. Human proteins include proteins previously detected in liver-derived exosomes, highlighting the potential of this chimeric mouse model to study plasma exosomes derived unequivocally from human hepatocytes. Noticeably, we identified 17 parasite proteins including enzymes, surface proteins, components of the endocytic pathway and translation machinery, as well as uncharacterized proteins. Western blot analysis validated the presence of human arginase-I and an uncharacterized P. vivax protein in plasma-derived exosomes. This study represents a proof-of-principle that plasma-derived exosomes from P. vivax infected FRG-huHep mice contain human hepatocyte and P. vivax proteins with the potential to unveil biological features of liver infection and identify biomarkers of hypnozoite infection.

Circulating miRNAs, isomiRs and small RNA clusters in human plasma and breast milk

Circulating small RNAs, including miRNAs but also isomiRs and other RNA species, have the potential to be used as non-invasive biomarkers for communicable and non-communicable diseases. This study aims to characterize and compare small RNA profiles in human biofluids. For this purpose, RNA was extracted from plasma and breast milk samples from 15 healthy postpartum mothers. Small RNA libraries were prepared with the NEBNext® small RNA library preparation kit and sequenced in an Illumina HiSeq2000 platform. miRNAs, isomiRs and clusters of small RNAs were annotated using seqBuster/seqCluster framework in 5 plasma and 10 milk samples that passed the initial quality control. The RNA yield was 81 ng/mL [standard deviation (SD): 41] and 3985 ng/mL (SD: 3767) for plasma and breast milk, respectively. Mean number of good quality reads was 4.04 million (M) (40.01% of the reads) in plasma and 12.5M (89.6%) in breast milk. One thousand one hundred eighty two miRNAs, 12,084 isomiRs and 1,053 small RNA clusters that included piwi-interfering RNAs (piRNAs), tRNAs, small nucleolar RNAs (snoRNA) and small nuclear RNAs (snRNAs) were detected. Samples grouped by biofluid, with 308 miRNAs, 1,790 isomiRs and 778 small RNA clusters differentially detected. In summary, plasma and milk showed a different small RNA profile. In both, miRNAs, piRNAs, tRNAs, snRNAs, and snoRNAs were identified, confirming the presence of non-miRNA species in plasma, and describing them for the first time in milk.

Impact of Healthy Vending Machine Options in a Large Community Health Organization

PURPOSE

To determine whether increasing the proportion of healthier options in vending machines decreases the amount of calories, fat, sugar, and sodium vended, while maintaining total sales revenue.

DESIGN

This study evaluated the impact of altering nutritious options to vending machines throughout the Banner Health organization by comparing vended items' sales and nutrition information over 6 months compared to the same 6 months of the previous year.

SETTING

Twenty-three locations including corporate and patient-care centers.

INTERVENTION

Changing vending machine composition toward more nutritious options.

MEASURES

Comparisons of monthly aggregates of sales, units vended, calories, fat, sodium, and sugar vended by site.

ANALYSIS

A pre-post analysis using paired t tests comparing 6 months before implementation to the equivalent 6 months postimplementation.

RESULTS

Significant average monthly decreases were seen for calories (16.7%, P = .002), fat (27.4%, P ≤ .0001), sodium (25.9%, P ≤ .0001), and sugar (11.8%, P = .045) vended from 2014 to 2015. Changes in revenue and units vended did not change from 2014 to 2015 ( P = .58 and P = .45, respectively).

CONCLUSION

Increasing the proportion of healthier options in vending machines from 20% to 80% significantly lowered the amount of calories, sodium, fat, and sugar vended, while not reducing units vended or having a negative financial impact.

Genes Interacting with Occupational Exposures to Low Molecular Weight Agents and Irritants on Adult-Onset Asthma in Three European Studies

BACKGROUND

The biological mechanisms by which cleaning products and disinfectants-an emerging risk factor-affect respiratory health remain incompletely evaluated. Studying genes by environment interactions (G × E) may help identify new genes related to adult-onset asthma.

OBJECTIVES

We identified interactions between genetic polymorphisms of a large set of genes involved in the response to oxidative stress and occupational exposures to low molecular weight (LMW) agents or irritants on adult-onset asthma.

METHODS

Our data came from three large European cohorts: Epidemiological Family-based Study of the Genetics and Environment of Asthma (EGEA), Swiss Cohort Study on Air Pollution and Lung and Heart Disease in Adults (SAPALDIA), and European Community Respiratory Health Survey in Adults (ECRHS). A candidate pathway-based strategy identified 163 genes involved in the response to oxidative stress and potentially related to exposures to LMW agents/irritants. Occupational exposures were evaluated using an asthma job-exposure matrix and job-specific questionnaires for cleaners and healthcare workers. Logistic regression models were used to detect G × E interactions, adjusted for age, sex, and population ancestry, in 2,599 adults (mean age, 47 years; 60% women, 36% exposed, 18% asthmatics). p-Values were corrected for multiple comparisons.

RESULTS

Ever exposure to LMW agents/irritants was associated with current adult-onset asthma [OR = 1.28 (95% CI: 1.04, 1.58)]. Eight single nucleotide polymorphism (SNP) by exposure interactions at five loci were found at p < 0.005: PLA2G4A (rs932476, chromosome 1), near PLA2R1 (rs2667026, chromosome 2), near RELA (rs931127, rs7949980, chromosome 11), PRKD1 (rs1958980, rs11847351, rs1958987, chromosome 14), and PRKCA (rs6504453, chromosome 17). Results were consistent across the three studies and after accounting for smoking.

CONCLUSIONS

Using a pathway-based selection process, we identified novel genes potentially involved in adult asthma by interaction with occupational exposure. These genes play a role in the NF-κB pathway, which is involved in inflammation. Citation: Rava M, Ahmed I, Kogevinas M, Le Moual N, Bouzigon E, Curjuric I, Dizier MH, Dumas O, Gonzalez JR, Imboden M, Mehta AJ, Tubert-Bitter P, Zock JP, Jarvis D, Probst-Hensch NM, Demenais F, Nadif R. 2017. Genes interacting with occupational exposures to low molecular weight agents and irritants on adult-onset asthma in three European studies. Environ Health Perspect 125:207-214; http://dx.doi.org/10.1289/EHP376.

Efficient and Powerful Method for Combining P-Values in Genome-Wide Association Studies

The goal of Genome-wide Association Studies (GWAS) is the identification of genetic variants, usually single nucleotide polymorphisms (SNPs), that are associated with disease risk. However, SNPs detected so far with GWAS for most common diseases only explain a small proportion of their total heritability. Gene set analysis (GSA) has been proposed as an alternative to single-SNP analysis with the aim of improving the power of genetic association studies. Nevertheless, most GSA methods rely on expensive computational procedures that make unfeasible their implementation in GWAS. We propose a new GSA method, referred as globalEVT, which uses the extreme value theory to derive gene-level p-values. GlobalEVT reduces dramatically the computational requirements compared to other GSA approaches. In addition, this new approach improves the power by allowing different inheritance models for each genetic variant as illustrated in the simulation study performed and allows the existence of correlation between the SNPs. Real data analysis of an Attention-deficit/hyperactivity disorder (ADHD) study illustrates the importance of using GSA approaches for exploring new susceptibility genes. Specifically, the globalEVT method is able to detect genes related to Cyclophilin A like domain proteins which is known to play an important role in the mechanisms of ADHD development.

Heritability and Genome-Wide Association Analyses of Sleep Duration in Children: The EAGLE Consortium

STUDY OBJECTIVES

Low or excessive sleep duration has been associated with multiple outcomes, but the biology behind these associations remains elusive. Specifically, genetic studies in children are scarce. In this study, we aimed to: (1) estimate the proportion of genetic variance of sleep duration in children attributed to common single nucleotide polymorphisms (SNPs), (2) identify novel SNPs associated with sleep duration in children, and (3) investigate the genetic overlap of sleep duration in children and related metabolic and psychiatric traits.

METHODS

We performed a population-based molecular genetic study, using data form the EArly Genetics and Life course Epidemiology (EAGLE) Consortium. 10,554 children of European ancestry were included in the discovery, and 1,250 children in the replication phase.

RESULTS

We found evidence of significant but modest SNP heritability of sleep duration in children (SNP h² 0.14, 95% CI [0.05, 0.23]) using the LD score regression method. A novel region at chromosome 11q13.4 (top SNP: rs74506765, P = 2.27e-08) was associated with sleep duration in children, but this was not replicated in independent studies. Nominally significant genetic overlap was only found (rG = 0.23, P = 0.05) between sleep duration in children and type 2 diabetes in adults, supporting the hypothesis of a common pathogenic mechanism.

CONCLUSIONS

The significant SNP heritability of sleep duration in children and the suggestive genetic overlap with type 2 diabetes support the search for genetic mechanisms linking sleep duration in children to multiple outcomes in health and disease.

Genetic variants associated with subjective well-being, depressive symptoms, and neuroticism identified through genome-wide analyses

Very few genetic variants have been associated with depression and neuroticism, likely because of limitations on sample size in previous studies. Subjective well-being, a phenotype that is genetically correlated with both of these traits, has not yet been studied with genome-wide data. We conducted genome-wide association studies of three phenotypes: subjective well-being (n = 298,420), depressive symptoms (n = 161,460), and neuroticism (n = 170,911). We identify 3 variants associated with subjective well-being, 2 variants associated with depressive symptoms, and 11 variants associated with neuroticism, including 2 inversion polymorphisms. The two loci associated with depressive symptoms replicate in an independent depression sample. Joint analyses that exploit the high genetic correlations between the phenotypes (|ρ^| ≈ 0.8) strengthen the overall credibility of the findings and allow us to identify additional variants. Across our phenotypes, loci regulating expression in central nervous system and adrenal or pancreas tissues are strongly enriched for association.

DNA Methylation in Newborns and Maternal Smoking in Pregnancy: Genome-wide Consortium Meta-analysis

Epigenetic modifications, including DNA methylation, represent a potential mechanism for environmental impacts on human disease. Maternal smoking in pregnancy remains an important public health problem that impacts child health in a myriad of ways and has potential lifelong consequences. The mechanisms are largely unknown, but epigenetics most likely plays a role. We formed the Pregnancy And Childhood Epigenetics (PACE) consortium and meta-analyzed, across 13 cohorts (n = 6,685), the association between maternal smoking in pregnancy and newborn blood DNA methylation at over 450,000 CpG sites (CpGs) by using the Illumina 450K BeadChip. Over 6,000 CpGs were differentially methylated in relation to maternal smoking at genome-wide statistical significance (false discovery rate, 5%), including 2,965 CpGs corresponding to 2,017 genes not previously related to smoking and methylation in either newborns or adults. Several genes are relevant to diseases that can be caused by maternal smoking (e.g., orofacial clefts and asthma) or adult smoking (e.g., certain cancers). A number of differentially methylated CpGs were associated with gene expression. We observed enrichment in pathways and processes critical to development. In older children (5 cohorts, n = 3,187), 100% of CpGs gave at least nominal levels of significance, far more than expected by chance (p value < 2.2 × 10(-16)). Results were robust to different normalization methods used across studies and cell type adjustment. In this large scale meta-analysis of methylation data, we identified numerous loci involved in response to maternal smoking in pregnancy with persistence into later childhood and provide insights into mechanisms underlying effects of this important exposure.

DNA methylation levels and long-term trihalomethane exposure in drinking water: an epigenome-wide association study

Trihalomethanes (THM) are undesired disinfection byproducts (DBPs) formed during water treatment. Mice exposed to DBPs showed global DNA hypomethylation and c-myc and c-jun gene-specific hypomethylation, while evidence of epigenetic effects in humans is scarce. We explored the association between lifetime THM exposure and DNA methylation through an epigenome-wide association study. We selected 138 population-based controls from a case-control study of colorectal cancer conducted in Barcelona, Spain, exposed to average lifetime THM levels ≤85 μg/L vs. >85 μg/L (N = 68 and N = 70, respectively). Mean age of participants was 70 years, and 54% were male. Average lifetime THM level in the exposure groups was 64 and 130 µg/L, respectively. DNA was extracted from whole blood and was bisulphite converted to measure DNA methylation levels using the Illumina HumanMethylation450 BeadChip. Data preprocessing was performed using RnBeads. Methylation was compared between exposure groups using empirical Bayes moderated linear regression for CpG sites and Gaussian kernel for CpG regions. ConsensusPathDB was used for gene set enrichment. Statistically significant differences in methylation between exposure groups was found in 140 CpG sites and 30 gene-related regions, after false discovery rate <0.05 and adjustment for age, sex, methylation first principal component, and blood cell proportion. The annotated genes were localized to several cancer pathways. Among them, 29 CpGs had methylation levels associated with THM levels (|Δβ|≥0.05) located in 11 genes associated with cancer in other studies. Our results suggest that THM exposure may affect DNA methylation in genes related to tumors, including colorectal and bladder cancers. Future confirmation studies are required.

Genome-wide association analysis identifies three new susceptibility loci for childhood body mass index

A large number of genetic loci are associated with adult body mass index. However, the genetics of childhood body mass index are largely unknown. We performed a meta-analysis of genome-wide association studies of childhood body mass index, using sex- and age-adjusted standard deviation scores. We included 35 668 children from 20 studies in the discovery phase and 11 873 children from 13 studies in the replication phase. In total, 15 loci reached genome-wide significance (P-value < 5 × 10(-8)) in the joint discovery and replication analysis, of which 12 are previously identified loci in or close to ADCY3, GNPDA2, TMEM18, SEC16B, FAIM2, FTO, TFAP2B, TNNI3K, MC4R, GPR61, LMX1B and OLFM4 associated with adult body mass index or childhood obesity. We identified three novel loci: rs13253111 near ELP3, rs8092503 near RAB27B and rs13387838 near ADAM23. Per additional risk allele, body mass index increased 0.04 Standard Deviation Score (SDS) [Standard Error (SE) 0.007], 0.05 SDS (SE 0.008) and 0.14 SDS (SE 0.025), for rs13253111, rs8092503 and rs13387838, respectively. A genetic risk score combining all 15 SNPs showed that each additional average risk allele was associated with a 0.073 SDS (SE 0.011, P-value = 3.12 × 10(-10)) increase in childhood body mass index in a population of 1955 children. This risk score explained 2% of the variance in childhood body mass index. This study highlights the shared genetic background between childhood and adult body mass index and adds three novel loci. These loci likely represent age-related differences in strength of the associations with body mass index.

Principal Component Analysis of the Effects of Environmental Enrichment and (-)-epigallocatechin-3-gallate on Age-Associated Learning Deficits in a Mouse Model of Down Syndrome

Down syndrome (DS) individuals present increased risk for Alzheimer's disease (AD) neuropathology and AD-type dementia. Here, we investigated the use of green tea extracts containing (-)-epigallocatechin-3-gallate (EGCG), as co-adjuvant to enhance the effects of environmental enrichment (EE) in Ts65Dn mice, a segmental trisomy model of DS that partially mimics DS/AD pathology, at the age of initiation of cognitive decline. Classical repeated measures ANOVA showed that combined EE-EGCG treatment was more efficient than EE or EGCG alone to improve specific spatial learning related variables. Using principal component analysis (PCA) we found that several spatial learning parameters contributed similarly to a first PC and explained a large proportion of the variance among groups, thus representing a composite learning measure. This PC1 revealed that EGCG or EE alone had no significant effect. However, combined EE-EGCG significantly ameliorated learning alterations of middle age Ts65Dn mice. Interestingly, PCA revealed an increased variability along learning sessions with good and poor learners in Ts65Dn, and this stratification did not disappear upon treatments. Our results suggest that combining EE and EGCG represents a viable therapeutic approach for amelioration of age-related cognitive decline in DS, although its efficacy may vary across individuals.

Multi-ancestry genome-wide association study of 21,000 cases and 95,000 controls identifies new risk loci for atopic dermatitis

Genetic association studies have identified 21 loci associated with atopic dermatitis risk predominantly in populations of European ancestry. To identify further susceptibility loci for this common complex skin disease, we performed a meta-analysis of >15 million genetic variants in 21,399 cases and 95,464 controls from populations of European, African, Japanese and Latino ancestry, followed by replication in 32,059 cases and 228,628 controls from 18 studies. We identified 10 novel risk loci, bringing the total number of known atopic dermatitis risk loci to 31 (with novel secondary signals at 4 of these). Notably, the new loci include candidate genes with roles in regulation of innate host defenses and T-cell function, underscoring the important contribution of (auto-)immune mechanisms to atopic dermatitis pathogenesis.

The Pregnancy Exposome: Multiple Environmental Exposures in the INMA-Sabadell Birth Cohort

The "exposome" is defined as "the totality of human environmental exposures from conception onward, complementing the genome" and its holistic approach may advance understanding of disease etiology. We aimed to describe the correlation structure of the exposome during pregnancy to better understand the relationships between and within families of exposure and to develop analytical tools appropriate to exposome data. Estimates on 81 environmental exposures of current health concern were obtained for 728 women enrolled in The INMA (INfancia y Medio Ambiente) birth cohort, in Sabadell, Spain, using biomonitoring, geospatial modeling, remote sensors, and questionnaires. Pair-wise Pearson's and polychoric correlations were calculated and principal components were derived. The median absolute correlation across all exposures was 0.06 (5th-95th centiles, 0.01-0.54). There were strong levels of correlation within families of exposure (median = 0.45, 5th-95th centiles, 0.07-0.85). Nine exposures (11%) had a correlation higher than 0.5 with at least one exposure outside their exposure family. Effectively all the variance in the data set (99.5%) was explained by 40 principal components. Future exposome studies should interpret exposure effects in light of their correlations to other exposures. The weak to moderate correlation observed between exposure families will permit adjustment for confounding in future exposome studies.

A novel common variant in DCST2 is associated with length in early life and height in adulthood

Common genetic variants have been identified for adult height, but not much is known about the genetics of skeletal growth in early life. To identify common genetic variants that influence fetal skeletal growth, we meta-analyzed 22 genome-wide association studies (Stage 1; N = 28 459). We identified seven independent top single nucleotide polymorphisms (SNPs) (P < 1 × 10(-6)) for birth length, of which three were novel and four were in or near loci known to be associated with adult height (LCORL, PTCH1, GPR126 and HMGA2). The three novel SNPs were followed-up in nine replication studies (Stage 2; N = 11 995), with rs905938 in DC-STAMP domain containing 2 (DCST2) genome-wide significantly associated with birth length in a joint analysis (Stages 1 + 2; β = 0.046, SE = 0.008, P = 2.46 × 10(-8), explained variance = 0.05%). Rs905938 was also associated with infant length (N = 28 228; P = 5.54 × 10(-4)) and adult height (N = 127 513; P = 1.45 × 10(-5)). DCST2 is a DC-STAMP-like protein family member and DC-STAMP is an osteoclast cell-fusion regulator. Polygenic scores based on 180 SNPs previously associated with human adult stature explained 0.13% of variance in birth length. The same SNPs explained 2.95% of the variance of infant length. Of the 180 known adult height loci, 11 were genome-wide significantly associated with infant length (SF3B4, LCORL, SPAG17, C6orf173, PTCH1, GDF5, ZNFX1, HHIP, ACAN, HLA locus and HMGA2). This study highlights that common variation in DCST2 influences variation in early growth and adult height.

The human early-life exposome (HELIX): project rationale and design

BACKGROUND

Developmental periods in early life may be particularly vulnerable to impacts of environmental exposures. Human research on this topic has generally focused on single exposure-health effect relationships. The "exposome" concept encompasses the totality of exposures from conception onward, complementing the genome.

OBJECTIVES

The Human Early-Life Exposome (HELIX) project is a new collaborative research project that aims to implement novel exposure assessment and biomarker methods to characterize early-life exposure to multiple environmental factors and associate these with omics biomarkers and child health outcomes, thus characterizing the "early-life exposome." Here we describe the general design of the project.

METHODS

In six existing birth cohort studies in Europe, HELIX will estimate prenatal and postnatal exposure to a broad range of chemical and physical exposures. Exposure models will be developed for the full cohorts totaling 32,000 mother-child pairs, and biomarkers will be measured in a subset of 1,200 mother-child pairs. Nested repeat-sampling panel studies (n = 150) will collect data on biomarker variability, use smartphones to assess mobility and physical activity, and perform personal exposure monitoring. Omics techniques will determine molecular profiles (metabolome, proteome, transcriptome, epigenome) associated with exposures. Statistical methods for multiple exposures will provide exposure-response estimates for fetal and child growth, obesity, neurodevelopment, and respiratory outcomes. A health impact assessment exercise will evaluate risks and benefits of combined exposures.

CONCLUSIONS

HELIX is one of the first attempts to describe the early-life exposome of European populations and unravel its relation to omics markers and health in childhood. As proof of concept, it will form an important first step toward the life-course exposome.

Genome-wide association study of body mass index in 23 000 individuals with and without asthma

BACKGROUND

Both asthma and obesity are complex disorders that are influenced by environmental and genetic factors. Shared genetic factors between asthma and obesity have been proposed to partly explain epidemiological findings of co-morbidity between these conditions.

OBJECTIVE

To identify genetic variants that are associated with body mass index (BMI) in asthmatic children and adults, and to evaluate if there are differences between the genetics of BMI in asthmatics and healthy individuals.

METHODS

In total, 19 studies contributed with genome-wide analysis study (GWAS) data from more than 23 000 individuals with predominantly European descent, of whom 8165 are asthmatics.

RESULTS

We report associations between several DENND1B variants (P = 2.2 × 10(-7) for rs4915551) on chromosome 1q31 and BMI from a meta-analysis of GWAS data using 2691 asthmatic children (screening data). The top DENND1B single nucleotide polymorphisms(SNPs) were next evaluated in seven independent replication data sets comprising 2014 asthmatics, and rs4915551 was nominally replicated (P < 0.05) in two of the seven studies and of borderline significance in one (P = 0.059). However, strong evidence of effect heterogeneity was observed and overall, the association between rs4915551 and BMI was not significant in the total replication data set, P = 0.71. Using a random effects model, BMI was overall estimated to increase by 0.30 kg/m(2) (P = 0.01 for combined screening and replication data sets, N = 4705) per additional G allele of this DENND1BSNP. FTO was confirmed as an important gene for adult and childhood BMI regardless of asthma status.

CONCLUSIONS AND CLINICAL RELEVANCE

DENND1B was recently identified as an asthma susceptibility gene in a GWAS on children, and here, we find evidence that DENND1B variants may also be associated with BMI in asthmatic children. However, the association was overall not replicated in the independent data sets and the heterogeneous effect of DENND1B points to complex associations with the studied diseases that deserve further study.

A flexible count data model to fit the wide diversity of expression profiles arising from extensively replicated RNA-seq experiments

Background

High-throughput RNA sequencing (RNA-seq) offers unprecedented power to capture the real dynamics of gene expression. Experimental designs with extensive biological replication present a unique opportunity to exploit this feature and distinguish expression profiles with higher resolution. RNA-seq data analysis methods so far have been mostly applied to data sets with few replicates and their default settings try to provide the best performance under this constraint. These methods are based on two well-known count data distributions: the Poisson and the negative binomial. The way to properly calibrate them with large RNA-seq data sets is not trivial for the non-expert bioinformatics user.

Results

Here we show that expression profiles produced by extensively-replicated RNA-seq experiments lead to a rich diversity of count data distributions beyond the Poisson and the negative binomial, such as Poisson-Inverse Gaussian or Pólya-Aeppli, which can be captured by a more general family of count data distributions called the Poisson-Tweedie. The flexibility of the Poisson-Tweedie family enables a direct fitting of emerging features of large expression profiles, such as heavy-tails or zero-inflation, without the need to alter a single configuration parameter. We provide a software package for R called tweeDEseq implementing a new test for differential expression based on the Poisson-Tweedie family. Using simulations on synthetic and real RNA-seq data we show that tweeDEseq yields P-values that are equally or more accurate than competing methods under different configuration parameters. By surveying the tiny fraction of sex-specific gene expression changes in human lymphoblastoid cell lines, we also show that tweeDEseq accurately detects differentially expressed genes in a real large RNA-seq data set with improved performance and reproducibility over the previously compared methodologies. Finally, we compared the results with those obtained from microarrays in order to check for reproducibility.

Conclusions

RNA-seq data with many replicates leads to a handful of count data distributions which can be accurately estimated with the statistical model illustrated in this paper. This method provides a better fit to the underlying biological variability; this may be critical when comparing groups of RNA-seq samples with markedly different count data distributions. The tweeDEseq package forms part of the Bioconductor project and it is available for download at http://www.bioconductor.org.

Genome-wide meta-analysis of common variant differences between men and women

The male-to-female sex ratio at birth is constant across world populations with an average of 1.06 (106 male to 100 female live births) for populations of European descent. The sex ratio is considered to be affected by numerous biological and environmental factors and to have a heritable component. The aim of this study was to investigate the presence of common allele modest effects at autosomal and chromosome X variants that could explain the observed sex ratio at birth. We conducted a large-scale genome-wide association scan (GWAS) meta-analysis across 51 studies, comprising overall 114 863 individuals (61 094 women and 53 769 men) of European ancestry and 2 623 828 common (minor allele frequency >0.05) single-nucleotide polymorphisms (SNPs). Allele frequencies were compared between men and women for directly-typed and imputed variants within each study. Forward-time simulations for unlinked, neutral, autosomal, common loci were performed under the demographic model for European populations with a fixed sex ratio and a random mating scheme to assess the probability of detecting significant allele frequency differences. We do not detect any genome-wide significant (P < 5 × 10(-8)) common SNP differences between men and women in this well-powered meta-analysis. The simulated data provided results entirely consistent with these findings. This large-scale investigation across ~115 000 individuals shows no detectable contribution from common genetic variants to the observed skew in the sex ratio. The absence of sex-specific differences is useful in guiding genetic association study design, for example when using mixed controls for sex-biased traits.

Influence of fetal glutathione S-transferase copy number variants on adverse reproductive outcomes

A nested case-control association study was designed to investigate the influence of maternal and fetal copy number variants (CNVs) on reproductive outcomes. Genotypes of ten CNVs encompassing GST and CYP genes were assessed. Significant associations were only found for child CNV genotypes. In particular, the child GSTM1 insertion allele was associated with prematurity protection (odds ratio, 95% CI: 0.67, 0.51-0.89; P < 0.01), whereas the child GSTT2B insertion allele was associated with an increased risk of being small for gestational age (odds ratio, 95% CI: 1.33, 1.07-1.67; P = 0.01). The study highlights the role of the fetal genome in prenatal development and also the need to analyse CNVs in a systematic manner.

Genome-wide association study of lung function decline in adults with and without asthma

BACKGROUND

Genome-wide association studies have identified determinants of chronic obstructive pulmonary disease, asthma, and lung function level; however, none have addressed decline in lung function.

OBJECTIVE

We conducted the first genome-wide association study on the age-related decrease in FEV(1) and its ratio to forced vital capacity (FVC) stratified a priori by asthma status.

METHODS

Discovery cohorts included adults of European ancestry (1,441 asthmatic and 2,677 nonasthmatic participants: the Epidemiological Study on the Genetics and Environment of Asthma, the Swiss Cohort Study on Air Pollution and Lung and Heart Disease in Adults, and the European Community Respiratory Health Survey). The associations of FEV(1) and FEV(1)/FVC ratio decrease with 2.5 million single nucleotide polymorphisms (SNPs) were estimated. Thirty loci were followed up by in silico replication (1,160 asthmatic and 10,858 nonasthmatic participants: Atherosclerosis Risk in Communities, the Framingham Heart Study, the British 1958 Birth Cohort, and the Dutch Asthma Study).

RESULTS

Main signals identified differed between asthmatic and nonasthmatic participants. None of the SNPs reached genome-wide significance. The association between the height-related gene DLEU7 and FEV(1) decrease suggested for nonasthmatic participants in the discovery phase was replicated (discovery, P = 4.8 × 10(-6); replication, P = .03), and additional sensitivity analyses point to a relation to growth. The top ranking signal, TUSC3, which is associated with FEV(1)/FVC ratio decrease in asthmatic participants (P = 5.3 × 10(-8)), did not replicate. SNPs previously associated with cross-sectional lung function were not prominently associated with decline.

CONCLUSIONS

Genetic heterogeneity of lung function might be extensive. Our results suggest that genetic determinants of longitudinal and cross-sectional lung function differ and vary by asthma status.

Cure frailty models for survival data: Application to recurrences for breast cancer and to hospital readmissions for colorectal cancer

Owing to the natural evolution of a disease, several events often arise after a first treatment for the same subject. For example, patients with a primary invasive breast cancer and treated with breast conserving surgery may experience breast cancer recurrences, metastases or death. A certain proportion of subjects in the population who are not expected to experience the events of interest are considered to be ‘cured’ or non-susceptible. To model correlated failure time data incorporating a surviving fraction, we compare several forms of cure rate frailty models. In the first model already proposed non-susceptible patients are those who are not expected to experience the event of interest over a sufficiently long period of time. The other proposed models account for the possibility of cure after each event. We illustrate the cure frailty models with two data sets. First to analyse time-dependent prognostic factors associated with breast cancer recurrences, metastases, new primary malignancy and death. Second to analyse successive rehospitalizations of patients diagnosed with colorectal cancer. Estimates were obtained by maximization of likelihood using SAS proc NLMIXED for a piecewise constant hazards model. As opposed to the simple frailty model, the proposed methods demonstrate great potential in modelling multivariate survival data with long-term survivors (‘cured’ individuals).

Multiple correspondence discriminant analysis: an application to detect stratification in copy number variation

We illustrate the use of multiple correspondence analysis (MCA) to correct for population stratification of copy number alteration data. In addition, we propose the use of multiple correspondence discriminant analysis (MCDA) to identify an optimal set of copy number variants (CNVs) that correctly infers the population stratification of a CNV map. Within MCDA, we highlight the novel use of correlation with class directions for variable ranking. We found a set of 20 CNVs with 98 per cent predictability in a CNV map of the HapMap populations. On this sample, the selection of variables based on centroid ranking outperformed the most common practice of ranking variables with their correlation to the principal axes.

CNVassoc: Association analysis of CNV data using R

Background

Copy number variants (CNV) are a potentially important component of the genetic contribution to risk of common complex diseases. Analysis of the association between CNVs and disease requires that uncertainty in CNV copy-number calls, which can be substantial, be taken into account; failure to consider this uncertainty can lead to biased results. Therefore, there is a need to develop and use appropriate statistical tools. To address this issue, we have developed CNVassoc, an R package for carrying out association analysis of common copy number variants in population-based studies. This package includes functions for testing for association with different classes of response variables (e.g. class status, censored data, counts) under a series of study designs (case-control, cohort, etc) and inheritance models, adjusting for covariates. The package includes functions for inferring copy number (CNV genotype calling), but can also accept copy number data generated by other algorithms (e.g. CANARY, CGHcall, IMPUTE).

Results

Here we present a new R package, CNVassoc, that can deal with different types of CNV arising from different platforms such as MLPA o aCGH. Through a real data example we illustrate that our method is able to incorporate uncertainty in the association process. We also show how our package can also be useful when analyzing imputed data when analyzing imputed SNPs. Through a simulation study we show that CNVassoc outperforms CNVtools in terms of computing time as well as in convergence failure rate.

Conclusions

We provide a package that outperforms the existing ones in terms of modelling flexibility, power, convergence rate, ease of covariate adjustment, and requirements for sample size and signal quality. Therefore, we offer CNVassoc as a method for routine use in CNV association studies.

Loss of function of transient receptor potential vanilloid 1 (TRPV1) genetic variant is associated with lower risk of active childhood asthma

Transient receptor potential cation channels of the vanilloid subfamily (TRPV) participate in the generation of Ca(2+) signals at different locations of the respiratory system, thereby controlling its correct functioning. TRPV1 expression and activity appear to be altered under pathophysiological conditions such as chronic cough and airway hypersensitivity, whereas TRPV4 single nucleotide polymorphisms (SNP) are associated with chronic obstructive pulmonary disease. However, to date, there is no information about the genetic impact of either TRPV1 or TRPV4 on asthma pathophysiology. We now report on the association of two functional SNPs, TRPV1-I585V and TRPV4-P19S, with childhood asthma. Both SNPs were genotyped in a population of 470 controls without respiratory symptoms and 301 asthmatics. Although none of the SNPs modified the risk of suffering from asthma, carriers of the TRPV1-I585V genetic variant showed a lower risk of current wheezing (odds ratio = 0.51; p = 0.01), a characteristic of active asthma, or cough (odds ratio = 0.57; p = 0.02). Functional analysis of TRPV1-I585V, using the Ca(2+)-sensitive dye fura-2 to measure intracellular [Ca(2+)] concentrations, revealed a decreased channel activity in response to two typical TRPV1 stimuli, heat and capsaicin. On the other hand, TRPV4-P19S, despite its loss-of-channel function, showed no significant association with asthma or the presence of wheezing. Our data suggest that genetically determined level of TRPV1 activity is relevant for asthma pathophysiology.

Positionally cloned genes and age-specific effects in asthma and atopy: an international population-based cohort study (ECRHS)

BACKGROUND

Several genes identified by positional cloning have been associated with asthma and atopy, but few findings have been replicated. Age at onset of asthma has been associated with different phenotypic characteristics, and with variants at chromosome 17q21 identified through genome-wide association. This study examined the associations and age-specific effects on asthma, atopy and bronchial hyper-responsiveness (BHR) of five candidate genes previously identified by positional cloning (ADAM33, PHF11, NPSR1, DPP10, SPINK5).

METHODS

51 polymorphisms from 2474 participants from 13 countries who took part in the European Community Respiratory Health Survey (1990-2000) were studied. Asthma and age at onset of asthma were assessed by questionnaire data, BHR by methacholine challenge and atopy by specific immunoglobulin E to four common allergens.

RESULTS

Significant associations with asthma, atopy and particularly for asthma with atopy were observed for a large region of 47 kb in the NPSR1 gene, even after Bonferroni correction for multiple comparisons (p<0.001). The associations with NPSR1 were stronger in those reporting a first attack of asthma before the age of 15, with statistically significant interactions with age of onset found for three SNPs. The evidence for ADAM33 and BHR and for an age-specific effect of two SNPs in DPP10 and asthma was weaker.

CONCLUSION

This study provides further evidence for an effect of NPSR1 on asthma, atopy and atopic asthma. In addition, this analysis suggests a role for NPSR1 in early-onset asthma driven by the strong effect of this gene on atopic asthma.

Joint effect of obesity and TNFA variability on asthma: two international cohort studies

Obesity is a risk factor for asthma. Adipose tissue expresses pro-inflammatory molecules including tumour necrosis factor (TNF), and levels of TNF are also related to polymorphisms in the TNF-alpha (TNFA) gene. The current authors examined the joint effect of obesity and TNFA variability on asthma in adults by combining two population-based studies. The European Community Respiratory Health Survey and the Swiss Cohort Study on Air Pollution and Lung and Heart Disease in Adults used comparable protocols, questionnaires and measures of lung function and atopy. DNA samples from 9,167 participants were genotyped for TNFA -308 and lymphotoxin-alpha (LTA) +252 gene variants. Obesity and TNFA were associated with asthma when mutually adjusting for their independent effects (odds ratio (OR) for obesity 2.4, 95% confidence interval (CI) 1.7-3.2; OR for TNFA -308 polymorphism 1.3, 95% CI 1.1-1.6). The association of obesity with asthma was stronger for subjects carrying the G/A and A/A TNFA -308 genotypes compared with the more common G/G genotype, particularly among nonatopics (OR for G/A and A/A genotypes 6.1, 95% CI 2.5-14.4; OR for G/G genotype 1.7, 95% CI 0.8-3.3). The present findings provide, for the first time, evidence for a complex pattern of interaction between obesity, a pro-inflammatory genetic factor and asthma.

Brain-derived neurotrophic factor modulates the severity of cognitive alterations induced by mutant huntingtin: involvement of phospholipaseCgamma activity and glutamate receptor expression

The involvement of brain-derived neurotrophic factor (BDNF) in cognitive processes and the decrease in its expression in Huntington's disease suggest that this neurotrophin may play a role in learning impairment during the disease progression. We therefore analyzed the onset and severity of cognitive deficits in two different mouse models with the same mutant huntingtin but with different levels of BDNF (R6/1 and R6/1:BDNF+/- mice). We observed that BDNF modulates cognitive function in different learning tasks, even before the onset of motor symptoms. R6/1:BDNF+/- mice showed earlier and more accentuated cognitive impairment than R6/1 mice at 5 weeks of age in discrimination learning; at 5 weeks of age in procedural learning; and at 9 weeks of age in alternation learning. At the earliest age at which cognitive impairment was detected, electrophysiological analysis was performed in the hippocampus. All mutant genotypes showed reduced hippocampal long term potentiation (LTP) with respect to wild type but did not show differences between them. Thus, we evaluated the involvement of BDNF-trkB signaling and glutamate receptor expression in the hippocampus of these mice. We observed a decrease in phospholipaseCgamma activity, but not ERK, in R61, BDNF+/- and R6/1:BDNF+/- hippocampus at the age when LTP was altered. However, a specific decrease in the expression of glutamate receptors NR1, NR2A and GluR1 was detected only in R6/1:BDNF+/- hippocampus. Therefore, these results show that BDNF modulates the learning and memory alterations induced by mutant huntingtin. This interaction leads to intracellular changes, such as specific changes in glutamate receptors and in BDNF-trkB signaling through phospholipaseCgamma.