Genetic interactions affecting human gene expression identified by variance association mapping

  1. Andrew A Brown
  2. Alfonso Buil
  3. Ana Viñuela
  4. Tuuli Lappalainen
  5. Hou-Feng Zheng
  6. John B Richards
  7. Kerrin S Small
  8. Timothy D Spector
  9. Emmanouil T Dermitzakis
  10. Richard Durbin  Is a corresponding author
  1. Wellcome Trust Sanger Institute, United Kingdom
  2. University of Geneva, Switzerland
  3. King's College London, United Kingdom
  4. McGill University, Canada
  5. The Wellcome Trust Sanger Institute, United Kingdom

Abstract

Non-additive interaction between genetic variants, or epistasis, is a possible explanation for the gap between heritability of complex traits and the variation explained by identified genetic loci. Interactions give rise to genotype dependent variance, and therefore the identification of variance quantitative trait loci can be an intermediate step to discover both epistasis and gene by environment effects (GxE). Using RNA-sequence data from lymphoblastoid cell lines (LCLs) from the TwinsUK cohort, we identify a candidate set of 508 variance associated SNPs. Exploiting the twin design we show that GxE plays a role in ~70% of these associations. Further investigation of these loci reveals 57 epistatic interactions that replicated in a smaller dataset, explaining on average 4.3% of phenotypic variance. In 24 cases, more variance is explained by the interaction than their additive contributions. Using molecular phenotypes in this way may provide a route to uncovering genetic interactions underlying more complex traits.

Article and author information

Author details

  1. Andrew A Brown

    Wellcome Trust Sanger Institute, Cambridge, United Kingdom
    Competing interests
    No competing interests declared.
  2. Alfonso Buil

    University of Geneva, Geneva, Switzerland
    Competing interests
    No competing interests declared.
  3. Ana Viñuela

    King's College London, London, United Kingdom
    Competing interests
    No competing interests declared.
  4. Tuuli Lappalainen

    University of Geneva, Geneva, Switzerland
    Competing interests
    No competing interests declared.
  5. Hou-Feng Zheng

    McGill University, Montreal, Canada
    Competing interests
    No competing interests declared.
  6. John B Richards

    King's College London, London, United Kingdom
    Competing interests
    No competing interests declared.
  7. Kerrin S Small

    King's College London, London, United Kingdom
    Competing interests
    No competing interests declared.
  8. Timothy D Spector

    King's College London, London, United Kingdom
    Competing interests
    No competing interests declared.
  9. Emmanouil T Dermitzakis

    University of Geneva, Geneva, Switzerland
    Competing interests
    Emmanouil T Dermitzakis, Reviewing editor, eLife.
  10. Richard Durbin

    The Wellcome Trust Sanger Institute, Cambridge, United Kingdom
    For correspondence
    rd@sanger.ac.uk
    Competing interests
    No competing interests declared.

Reviewing Editor

  1. Philipp Khaitovich, Partner Institute for Computational Biology, China

Ethics

Human subjects: This project was approved by the ethics committee at St Thomas' Hospital London, where all the biopsies were carried out. Volunteers gave informed consent and signed an approved consent form prior to the biopsy procedure. Volunteers were supplied with an appropriate detailed information sheet regarding the research project and biopsy procedure by post prior to attending for the biopsy. The St. Thomas' Research Ethics Committee (REC) approved on 20th September 2007 the protocol for dissemination of data, including DNA, with the REC reference number RE04/015. On 12th of March of 2008, the St Thomas' REC confirmed this approval extends to expression data.

Version history

  1. Received: August 21, 2013
  2. Accepted: March 13, 2014
  3. Accepted Manuscript published: April 25, 2014 (version 1)
  4. Version of Record published: May 20, 2014 (version 2)

Copyright

© 2014, Brown et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 7,268
    views
  • 486
    downloads
  • 102
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Andrew A Brown
  2. Alfonso Buil
  3. Ana Viñuela
  4. Tuuli Lappalainen
  5. Hou-Feng Zheng
  6. John B Richards
  7. Kerrin S Small
  8. Timothy D Spector
  9. Emmanouil T Dermitzakis
  10. Richard Durbin
(2014)
Genetic interactions affecting human gene expression identified by variance association mapping
eLife 3:e01381.
https://doi.org/10.7554/eLife.01381

Share this article

https://doi.org/10.7554/eLife.01381

Further reading

    1. Chromosomes and Gene Expression
    Xinyang Bing, Wenfan Ke ... James B Jaynes
    Research Article

    Two different models have been proposed to explain how the endpoints of chromatin looped domains (‘TADs’) in eukaryotic chromosomes are determined. In the first, a cohesin complex extrudes a loop until it encounters a boundary element roadblock, generating a stem-loop. In this model, boundaries are functionally autonomous: they have an intrinsic ability to halt the movement of incoming cohesin complexes that is independent of the properties of neighboring boundaries. In the second, loops are generated by boundary:boundary pairing. In this model, boundaries are functionally non-autonomous, and their ability to form a loop depends upon how well they match with their neighbors. Moreover, unlike the loop-extrusion model, pairing interactions can generate both stem-loops and circle-loops. We have used a combination of MicroC to analyze how TADs are organized, and experimental manipulations of the even skipped TAD boundary, homie, to test the predictions of the ‘loop-extrusion’ and the ‘boundary-pairing’ models. Our findings are incompatible with the loop-extrusion model, and instead suggest that the endpoints of TADs in flies are determined by a mechanism in which boundary elements physically pair with their partners, either head-to-head or head-to-tail, with varying degrees of specificity. Although our experiments do not address how partners find each other, the mechanism is unlikely to require loop extrusion.

    1. Chromosomes and Gene Expression
    Wenfan Ke, Miki Fujioka ... James B Jaynes
    Research Article

    The chromosomes in multicellular eukaryotes are organized into a series of topologically independent loops called TADs. In flies, TADs are formed by physical interactions between neighboring boundaries. Fly boundaries exhibit distinct partner preferences, and pairing interactions between boundaries are typically orientation-dependent. Pairing can be head-to-tail or head-to-head. The former generates a stem-loop TAD, while the latter gives a circle-loop TAD. The TAD that encompasses the Drosophila even skipped (eve) gene is formed by the head-to-tail pairing of the nhomie and homie boundaries. To explore the relationship between loop topology and the physical and regulatory landscape, we flanked the nhomie boundary region with two attP sites. The attP sites were then used to generate four boundary replacements: λ DNA, nhomie forward (WT orientation), nhomie reverse (opposite of WT orientation), and homie forward (same orientation as WT homie). The nhomie forward replacement restores the WT physical and regulatory landscape: in MicroC experiments, the eve TAD is a ‘volcano’ triangle topped by a plume, and the eve gene and its regulatory elements are sequestered from interactions with neighbors. The λ DNA replacement lacks boundary function: the endpoint of the ‘new’ eve TAD on the nhomie side is ill-defined, and eve stripe enhancers activate a nearby gene, eIF3j. While nhomie reverse and homie forward restore the eve TAD, the topology is a circle-loop, and this changes the local physical and regulatory landscape. In MicroC experiments, the eve TAD interacts with its neighbors, and the plume at the top of the eve triangle peak is converted to a pair of ‘clouds’ of contacts with the next-door TADs. Consistent with the loss of isolation afforded by the stem-loop topology, the eve enhancers weakly activate genes in the neighboring TADs. Conversely, eve function is partially disrupted.