There is no author summary for this article yet. Authors can add summaries to their articles on ScienceOpen to make them more accessible to a non-specialist audience.
Abstract
The identification of multiple signals at individual loci could explain additional
phenotypic variance (‘missing heritability’) of common traits, and help identify causal
genes. We examined gene expression levels as a model trait because of the large number
of strong genetic effects acting in
cis. Using expression profiles from 613 individuals, we performed genome-wide single
nucleotide polymorphism (SNP) analyses to identify
cis-expression quantitative trait loci (eQTLs), and conditional analysis to identify second
signals. We examined patterns of association when accounting for multiple SNPs at
a locus and when including additional SNPs from the 1000 Genomes Project. We identified
1298
cis-eQTLs at an approximate false discovery rate 0.01, of which 118 (9%) showed evidence
of a second independent signal. For this subset of 118 traits, accounting for two
signals resulted in an average 31% increase in phenotypic variance explained (Wilcoxon
P< 0.0001). The association of SNPs with
cis gene expression could increase, stay similar or decrease in significance when accounting
for linkage disequilibrium with second signals at the same locus. Pairs of SNPs increasing
in significance tended to have gene expression increasing alleles on opposite haplotypes,
whereas pairs of SNPs decreasing in significance tended to have gene expression increasing
alleles on the same haplotypes. Adding data from the 1000 Genomes Project showed that
apparently independent signals could be potentially explained by a single association
signal. Our results show that accounting for multiple variants at a locus will increase
the variance explained in a substantial fraction of loci, but that allelic heterogeneity
will be difficult to define without resequencing loci and functional work.
We have created a global map of the effects of polymorphism on gene expression in 400 children from families recruited through a proband with asthma. We genotyped 408,273 SNPs and identified expression quantitative trait loci from measurements of 54,675 transcripts representing 20,599 genes in Epstein-Barr virus-transformed lymphoblastoid cell lines. We found that 15,084 transcripts (28%) representing 6,660 genes had narrow-sense heritabilities (H2) > 0.3. We executed genome-wide association scans for these traits and found peak lod scores between 3.68 and 59.1. The most highly heritable traits were markedly enriched in Gene Ontology descriptors for response to unfolded protein (chaperonins and heat shock proteins), regulation of progression through the cell cycle, RNA processing, DNA repair, immune responses and apoptosis. SNPs that regulate expression of these genes are candidates in the study of degenerative diseases, malignancy, infection and inflammation. We have created a downloadable database to facilitate use of our findings in the mapping of complex disease loci.
Studies correlating genetic variation to gene expression facilitate the interpretation of common human phenotypes and disease. As functional variants may be operating in a tissue-dependent manner, we performed gene expression profiling and association with genetic variants (single-nucleotide polymorphisms) on three cell types of 75 individuals. We detected cell type-specific genetic effects, with 69 to 80% of regulatory variants operating in a cell type-specific manner, and identified multiple expressive quantitative trait loci (eQTLs) per gene, unique or shared among cell types and positively correlated with the number of transcripts per gene. Cell type-specific eQTLs were found at larger distances from genes and at lower effect size, similar to known enhancers. These data suggest that the complete regulatory variant repertoire can only be uncovered in the context of cell-type specificity.
Introduction Efforts to fine map the causal variants responsible for genome-wide association studies (GWAS) signals have been largely predicated on the common disease common variant theory, postulating a common variant as the culprit for observed associations. This has led to extensive resequencing efforts that have been largely unsuccessful [1]–[5]. Here, we explore the possibility that part of the reason for this may be that the disease class causing an observed association may consist of multiple low-frequency variants across large regions of the genome—a phenomenon we call synthetic association. For convenience, these less common variants will be referred to here as “rare,” but we emphasize that we use this term loosely, only to refer to variants less common than those routinely studied in GWAS. The basic idea of how synthetic associations emerge in this model is illustrated in Figure 1, which shows how rare variants, by chance, can occur disproportionately in some parts of a gene genealogy. Any variant “higher up in the genealogy” that partitions those parts of the genealogy containing more disease variants than average will be identified as disease-associated. It is well appreciated that a noncausal variant will show association with a causal variant if the two are in strong linkage disequilibrium (LD). We use the previously introduced term synthetic association [6], however, to describe how such indirect association can occur between a common variant and at least one and possibly many rarer causal variants. Using the term synthetic as opposed to indirect emphasizes that the properties of the association signal are very different when the responsible variant or variants are much less frequent than the marker that carries the signal, as we detail below. 10.1371/journal.pbio.1000294.g001 Figure 1 Example genealogies showing causal variants and the strongest association for a common variant. (A) A genealogy with 10,000 original haplotypes was generated with 3,000 cases and 3,000 controls, genotype relative risk (γ) = 4, and nine causal variants. The branches containing the strongest synthetic association are indicated in blue. The branches containing the rare causal variants are in red. (B) A second genealogy was generated using the same parameters. These genealogies demonstrate two scenarios with genome-wide significant synthetic associations: the first (upper genealogy) had a high risk allele frequency (RAF = 0.49), and the second (lower genealogy) had a low RAF (0.08). To assess the tendency of rare disease-causing variants to create synthetic signals of association that are credited to single polymorphisms that are much more common in the population than the causal variants, we have simulated 10,000 haplotypes based on a coalescent model in a region either with or without recombination (Materials and Methods). We assumed that gene variants that influence disease have an allele frequency between 0.005 and 0.02, which is generally below the range of reliable detection (either by inclusion or indirect representation) using the genome-wide association platforms currently in use. We assumed a baseline probability of disease of φ for individuals with none of the rare genetic risk factors. The presence of at least one rare risk allele at the locus increased the probability of disease from φ to γ. We considered two values of φ (0.01, 0.1) and chose values of the penetrance γ such that the genotypic relative risk (GRR) of the rare causal variants varied incrementally between 2 and 6, where GRR is the ratio γ/φ. These values were chosen to explore the space around a GRR of 4, a threshold above which consistent linkage signals would be expected [7]. We simulated scenarios with one, three, five, seven, and nine rare causal variants. Results Across the conditions we have studied, not only is it possible to achieve genome-wide significance for common variants when one or more rare variants are the only contributors to disease, it is often the likely outcome (Figure 2). Overall, 30% of the simulations were able to detect an association with a common SNP at genome-wide significance (p 5%, Hardy-Weinberg equilibrium p-value >1×10−6, SNP call rate >95%), using the PLINK software [40]. For the sickle cell anemia GWAS, we compared 194 cases and 7,407 controls of inferred African ancestry via multidimensional scaling, with a genomic control inflation factor of 1.01. For hearing loss, we performed a GWAS on 418 cases and 6,892 control subjects, all of whom were of genetically inferred European ancestry via multidimensional scaling, with a genomic control inflation factor of 1.02.
[1
]Genetics of Complex Traits, Peninsula College of Medicine and Dentistry, simpleUniversity of Exeter , Exeter EX1 2LU, UK,
[2
]Laboratory of Neurogenetics, simpleNational Institute on Aging, National Institutes of Health , 35 Lincoln Drive, Bethesda, MD, USA,
[3
]Department of Molecular Neuroscience and Reta Lila Laboratories, simpleInstitute of Neurology, UCL , Queen Square House, Queen Square, London WC1N 3BG, UK,
[4
]simpleInstitute of Biomedical and Clinical Sciences, Peninsula College of Medicine and Dentistry,
University of Exeter , Barrack Road, Exeter EX2 5DW, UK,
[5
]Department of Epidemiology and Public Health, simpleUniversity of Maryland School of Medicine , Baltimore, MD, USA,
[6
]Geriatric Unit, simpleAzienda Sanitaria di Firenze , Florence, Italy and
[7
]simpleClinical Research Branch, National Institute on Aging NIA-ASTRA Unit, Harbor Hospital , MD, USA
Author notes
[*
]To whom correspondence should be addressed. Tel: +44 1392722935; Fax: +44 1392722926;
Email:
tim.frayling@
123456pms.ac.uk
This is an Open Access article distributed under the terms of the Creative Commons
Attribution Non-Commercial License (
http://creativecommons.org/licenses/by-nc/2.5), which permits unrestricted non-commercial use, distribution, and reproduction in
any medium, provided the original work is properly cited.
scite shows how a scientific paper has been cited by providing the context of the citation, a classification describing whether it supports, mentions, or contrasts the cited claim, and a label indicating in which section the citation was made.