Mapping gene-environment interactions at regulatory polymorphisms: insights into mechanisms of phenotypic variation

A computational investigation of cis-gene regulation in evolution

In biological processes involving gene networks, genes regulate other genes that determine the phenotypic traits. Gene regulation plays an important role in evolutionary dynamics. In a genetic algorithm, a trans-gene regulatory mechanism was shown to speed up adaptation and evolution. Here, we examine the effect of cis-gene regulation on an adaptive system. The model is haploid. A chromosome is partitioned into regulatory loci and structural loci. The regulatory genes regulate the expression and functioning of structural genes via the cis-elements in a probabilistic manner. In the simulation, the change in the allele frequency, the mean population fitness and the efficiency of phenotypic selection are monitored. Cis-gene regulation increases adaption and accelerates the evolutionary process in comparison with the case involving absence of gene regulation. Some special features of the simulation results are as follows. A low ratio of regulatory loci and structural loci gives higher adaptation for fixed total number of loci. Plasticity is advantageous beyond a threshold value. Adaptation is better for large number of total loci when the ratio of regulatory loci to structural loci is one. However, it reaches a saturation beyond which the increase in the total loci is not advantageous. Efficiency of the phenotypic selection is higher for larger value of the initial plasticity.

A signature of Neanderthal introgression on molecular mechanisms of environmental responses

Ancient human migrations led to the settlement of population groups in varied environmental contexts worldwide. The extent to which adaptation to local environments has shaped human genetic diversity is a longstanding question in human evolution. Recent studies have suggested that introgression of archaic alleles in the genome of modern humans may have contributed to adaptation to environmental pressures such as pathogen exposure. Functional genomic studies have demonstrated that variation in gene expression across individuals and in response to environmental perturbations is a main mechanism underlying complex trait variation. We considered gene expression response to in vitro treatments as a molecular phenotype to identify genes and regulatory variants that may have played an important role in adaptations to local environments. We investigated if Neanderthal introgression in the human genome may contribute to the transcriptional response to environmental perturbations. To this end we used eQTLs for genes differentially expressed in a panel of 52 cellular environments, resulting from 5 cell types and 26 treatments, including hormones, vitamins, drugs, and environmental contaminants. We found that SNPs with introgressed Neanderthal alleles (N-SNPs) disrupt binding of transcription factors important for environmental responses, including ionizing radiation and hypoxia, and for glucose metabolism. We identified an enrichment for N-SNPs among eQTLs for genes differentially expressed in response to 8 treatments, including glucocorticoids, caffeine, and vitamin D. Using Massively Parallel Reporter Assays (MPRA) data, we validated the regulatory function of 21 introgressed Neanderthal variants in the human genome, corresponding to 8 eQTLs regulating 15 genes that respond to environmental perturbations. These findings expand the set of environments where archaic introgression may have contributed to adaptations to local environments in modern humans and provide experimental validation for the regulatory function of introgressed variants.

Discovering in vivo cytokine-eQTL interactions from a lupus clinical trial

BACKGROUND

Cytokines are critical to human disease and are attractive therapeutic targets given their widespread influence on gene regulation and transcription. Defining the downstream regulatory mechanisms influenced by cytokines is central to defining drug and disease mechanisms. One promising strategy is to use interactions between expression quantitative trait loci (eQTLs) and cytokine levels to define target genes and mechanisms.

RESULTS

In a clinical trial for anti-IL-6 in patients with systemic lupus erythematosus, we measure interferon (IFN) status, anti-IL-6 drug exposure, and whole blood genome-wide gene expression at three time points. We show that repeat transcriptomic measurements increases the number of cis eQTLs identified compared to using a single time point. We observe a statistically significant enrichment of in vivo eQTL interactions with IFN status and anti-IL-6 drug exposure and find many novel interactions that have not been previously described. Finally, we find transcription factor binding motifs interrupted by eQTL interaction SNPs, which point to key regulatory mediators of these environmental stimuli and therefore potential therapeutic targets for autoimmune diseases. In particular, genes with IFN interactions are enriched for ISRE binding site motifs, while those with anti-IL-6 interactions are enriched for IRF4 motifs.

CONCLUSIONS

This study highlights the potential to exploit clinical trial data to discover in vivo eQTL interactions with therapeutically relevant environmental variables.

Determining the evolutionary forces shaping G × E

Contents Summary 31 I. Introduction 31 II. The maintenance of genetic variation for plasticity 32 III. Why is there environmental variation for genetic effects? 33 IV. Conclusions 35 Acknowledgements 35 References 35 SUMMARY: Phenotypic plasticity is common in nature, yet we lack a comprehensive understanding of the evolutionary forces that shape genetic variation for plasticity. This endeavor is especially important because variation for plasticity will result in genotype-by-environment interactions (G × E), a crucial component of variation in quantitative traits. Here, I review our understanding of the evolutionary forces shaping G × E, focusing specifically on: what evolutionary forces maintain variation for plasticity; and what forces maintain different genetic architectures across environments. My specific goal is to show that genomic data can be leveraged to explain the maintenance of G × E by contrasting patterns of genetic variation for plasticity with neutral expectations.

Interactions between genetic variation and cellular environment in skeletal muscle gene expression

From whole organisms to individual cells, responses to environmental conditions are influenced by genetic makeup, where the effect of genetic variation on a trait depends on the environmental context. RNA-sequencing quantifies gene expression as a molecular trait, and is capable of capturing both genetic and environmental effects. In this study, we explore opportunities of using allele-specific expression (ASE) to discover cis-acting genotype-environment interactions (GxE)—genetic effects on gene expression that depend on an environmental condition. Treating 17 common, clinical traits as approximations of the cellular environment of 267 skeletal muscle biopsies, we identify 10 candidate environmental response expression quantitative trait loci (reQTLs) across 6 traits (12 unique gene-environment trait pairs; 10% FDR per trait) including sex, systolic blood pressure, and low-density lipoprotein cholesterol. Although using ASE is in principle a promising approach to detect GxE effects, replication of such signals can be challenging as validation requires harmonization of environmental traits across cohorts and a sufficient sampling of heterozygotes for a transcribed SNP. Comprehensive discovery and replication will require large human transcriptome datasets, or the integration of multiple transcribed SNPs, coupled with standardized clinical phenotyping.

Combining genetic and nongenetic biomarkers to realize the promise of pharmacogenomics for inflammatory diseases

Many drugs used to treat inflammatory diseases are ineffective in a substantial proportion of patients. Identifying patients that are likely to respond to specific therapies would facilitate personalized treatment strategies that could improve outcomes while reducing costs and risks of adverse events. Despite these clear benefits, there are limited examples of predictive biomarkers of drug efficacy currently implemented into clinical practice for inflammatory diseases. We review efforts to identify genetic and nongenetic biomarkers of drug response in these diseases and consider potential benefits from combining multiple sources of biological data into multifeature predictive models.

The Red Queen's long race: human adaptation to pathogen pressure

Pathogens, and the infectious diseases they cause, have been paramount among the threats encountered by humans in their expansions throughout the globe. Numerous studies have identified immunity and host defence genes as being among the functions most strongly targeted by selection, most likely pathogen-driven. The dissection of the form and intensity of such selective pressures have increased our knowledge of the biological relevance of the underlying immunological mechanisms in host defence. Although the identities of the specific infectious agents imposing these selective pressures remain, in most cases, elusive, the impact of several pathogens, notably malaria and cholera, has been described. However, past selection against infectious diseases may have some fitness costs upon environmental changes, potentially leading to maladaptation and immunopathology.

Gene × smoking interactions on human brain gene expression: finding common mechanisms in adolescents and adults

BACKGROUND

Numerous studies have examined gene × environment interactions (G × E) in cognitive and behavioral domains. However, these studies have been limited in that they have not been able to directly assess differential patterns of gene expression in the human brain. Here, we assessed G × E interactions using two publically available datasets to assess if DNA variation is associated with post-mortem brain gene expression changes based on smoking behavior, a biobehavioral construct that is part of a complex system of genetic and environmental influences.

METHODS

We conducted an expression quantitative trait locus (eQTL) study on two independent human brain gene expression datasets assessing G × E for selected psychiatric genes and smoking status. We employed linear regression to model the significance of the Gene × Smoking interaction term, followed by meta-analysis across datasets.

RESULTS

Overall, we observed that the effect of DNA variation on gene expression is moderated by smoking status. Expression of 16 genes was significantly associated with single nucleotide polymorphisms that demonstrated G × E effects. The strongest finding (p = 1.9 × 10⁻¹¹) was neurexin 3-alpha (NRXN3), a synaptic cell-cell adhesion molecule involved in maintenance of neural connections (such as the maintenance of smoking behavior). Other significant G × E associations include four glutamate genes.

CONCLUSIONS

This is one of the first studies to demonstrate G × E effects within the human brain. In particular, this study implicated NRXN3 in the maintenance of smoking. The effect of smoking on NRXN3 expression and downstream behavior is different based upon SNP genotype, indicating that DNA profiles based on SNPs could be useful in understanding the effects of smoking behaviors. These results suggest that better measurement of psychiatric conditions, and the environment in post-mortem brain studies may yield an important avenue for understanding the biological mechanisms of G × E interactions in psychiatry.

Exploiting gene expression variation to capture gene-environment interactions for disease

Gene-environment interactions have long been recognized as a fundamental concept in evolutionary, quantitative, and medical genetics. In the genomics era, study of how environment and genome interact to shape gene expression variation is relevant to understanding the genetic architecture of complex phenotypes. While genetic analysis of gene expression variation focused on main effects, little is known about the extent of interaction effects implicating regulatory variants and their consequences on transcriptional variation. Here we survey the current state of the concept of transcriptional gene-environment interactions and discuss its utility for mapping disease phenotypes in light of the insights gained from genome-wide association studies of gene expression.