PNAS QnAs with NAS member and molecular biologist Steven A. Kliewer
New Research In
Physical Sciences
Featured Portals
Articles by Topic
Biological Sciences
Featured Portals
Articles by Topic
Contributed by Eric S. Lander, December 10, 2013 (sent for review September 23, 2013)
Significance
Discovering the genetic basis of common diseases, such as diabetes, heart disease, and schizophrenia, is a key goal in biomedicine. Genomic studies have revealed thousands of common genetic variants underlying disease, but these variants explain only a portion of the heritability. Rare variants are also likely to play an important role, but few examples are known thus far, and initial discovery efforts with small sample sizes have had only limited success. In this paper, we describe an analytical framework for the design of rare variant association studies of disease. It provides guidance with respect to sample size, as well as the roles of selection, disruptive and missense alleles, gene-specific allele frequency thresholds, isolated populations, gene sets, and coding vs. noncoding regions.
Abstract
Genetic studies have revealed thousands of loci predisposing to hundreds of human diseases and traits, revealing important biological pathways and defining novel therapeutic hypotheses. However, the genes discovered to date typically explain less than half of the apparent heritability. Because efforts have largely focused on common genetic variants, one hypothesis is that much of the missing heritability is due to rare genetic variants. Studies of common variants are typically referred to as genomewide association studies, whereas studies of rare variants are often simply called sequencing studies. Because they are actually closely related, we use the terms common variant association study (CVAS) and rare variant association study (RVAS). In this paper, we outline the similarities and differences between RVAS and CVAS and describe a conceptual framework for the design of RVAS. We apply the framework to address key questions about the sample sizes needed to detect association, the relative merits of testing disruptive alleles vs. missense alleles, frequency thresholds for filtering alleles, the value of predictors of the functional impact of missense alleles, the potential utility of isolated populations, the value of gene-set analysis, and the utility of de novo mutations. The optimal design depends critically on the selection coefficient against deleterious alleles and thus varies across genes. The analysis shows that common variant and rare variant studies require similarly large sample collections. In particular, a well-powered RVAS should involve discovery sets with at least 25,000 cases, together with a substantial replication set.
Footnotes
↵1Present address: Department of Statistics, The Hebrew University of Jerusalem, Mt. Scopus, Jerusalem 91905, Israel.
- ↵2To whom correspondence should be addressed. E-mail: lander{at}broadinstitute.org.
Author contributions: O.Z., E.H., M.J.D., B.M.N., S.R.S., and E.S.L. designed research; O.Z., S.F.S., K.S., R.D., E.H., B.M.N., S.R.S., and E.S.L. performed research; K.S., R.D., and S.K. contributed new reagents/analytic tools; O.Z., S.F.S., K.S., E.H., M.J.D., B.M.N., S.R.S., and E.S.L. analyzed data; and O.Z., M.J.D., B.M.N., S.R.S., and E.S.L. wrote the paper.
The authors declare no conflict of interest.
†We focus on risk to heterozygous carriers. Our calculations implicitly assume that the risk to individuals carrying two null alleles (λC*) is the same as the risk to heterozygous carriers (λC). Although such individuals may well have higher risk, they are much rarer than heterozygous carriers (because fC is small) and thus their impact on RVAS is typically negligible. [The effective relative risk is increased by fC(λC*/λC), which is ≪1 unless λC*/λC is huge; this case would essentially be a monogenic recessive trait.]
‡We focus on the number n of cases needed when the frequency fC in the population is known perfectly based on an (infinitely) large population survey. We make this assumption in the belief that very large datasets will become available in the coming years and that shared population controls can be used across studies. In the meanwhile, one can estimate fD within a study based on the frequency in either unaffected or random individuals. If a study involves cases and unaffecteds in proportion r and 1 − r, one requires approximately n/(1 − r) cases and n/r controls to detect association, where n is the number of cases given in Eq. 1 (SI Appendix, Section 3.3). For a balanced design, this corresponds to 2n cases and 2n controls.
§
¶The expected CAF for new neutral alleles born in a given generation is μC and thus the increase in the expected CAF over k generations cannot exceed kμC (and will be lower because many newborn alleles are lost). It follows that the collection of neutral alleles born since the onset of human population cannot increase the CAF of neutral alleles by more than 5% (= kμC/4μCNeq, where k =1,000 and Neq = 10,000). The actual increase is typically much smaller due to loss of newborn alleles.
||The relevant sampling frame for disease studies involves randomly selecting a chromosome and inspecting it for alleles. Alleles are thus sampled according to their frequency, yielding a frequency-weighted frequency distribution.
**One can also perform an analysis considering only those missense alleles with frequency ≤1% that are predicted to be probably damaging by PolyPhen-2 (that is, RVAS strategy 3). One observes 89 such missense alleles in cases vs. 32 in controls, which corresponds to an apparent excess relative risk λM(1%),PolyPhen2 = 2.5. This result agrees closely with the expectation of λM(1%),PolyPhen2 = 2.3, given the inferred value of s and the frequency threshold of 1%. For this type of analysis (assuming γnull = 80% and γneutral = 20%), the optimal threshold turns out to be T*= 0.15%, which would yield a much higher apparent effect size of λM(T*),PolyPhen2 = 8.2.
††The paper reports that the best nominal P value observed across the genome is 2.7 × 10−8 but does not address whether the value is significant. To correct for scanning the entire genome, one can apply extreme value theory for Orenstein-Uhlenbeck diffusions (SI Appendix, Section 11.1). The probability of such a P value arising by chance somewhere in the genome is ∼70%, and thus the observation is not statistically significant. Genomewide significance at the 5% level corresponds to P ∼ 2 × 10−9 (SI Appendix, Section 11.1).
This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1322563111/-/DCSupplemental.
Freely available online through the PNAS open access option.
Designing rare variant association studies
More Articles of This Classification
Related Content
Cited by...
- Taking Systems Medicine to Heart
- Statistical Validation of Rare Complement Variants Provides Insights into the Molecular Basis of Atypical Hemolytic Uremic Syndrome and C3 Glomerulopathy
- Evidence of directional and stabilizing selection in contemporary humans
- Identification of individuals by trait prediction using whole-genome sequencing data
- ANGPTL3 Deficiency and Protection Against Coronary Artery Disease
- The Effect of an Extreme and Prolonged Population Bottleneck on Patterns of Deleterious Variation: Insights from the Greenlandic Inuit
- A Genome-Wide Association Study of Resistance to Stripe Rust (Puccinia striiformis f. sp. tritici) in a Worldwide Collection of Hexaploid Spring Wheat (Triticum aestivum L.)
- Meiotic Interactors of a Mitotic Gene TAO3 Revealed by Functional Analysis of its Rare Variant
- Comprehensive population-based genome sequencing provides insight into hematopoietic regulatory mechanisms
- Mutation-Induced Functional Alterations of CCR6
- Distribution and clinical impact of functional variants in 50,726 whole-exome sequences from the DiscovEHR study
- Selection and explosive growth alter genetic architecture and hamper the detection of causal rare variants
- Quantitative Trait Loci and a Novel Genetic Candidate for Fear Learning
- Shroom3 contributes to the maintenance of the glomerular filtration barrier integrity
- Measuring missing heritability: Inferring the contribution of common variants
- Recent Developments in Cardiovascular Genetics and Genomics
- Identification of a Novel Gene for Diabetic Traits in Rats, Mice, and Humans
Similar Articles
You May Also be Interested in
Holographic acoustic tweezers, in which ultrasonic waves produced by arrays of sound emitters are used to individually manipulate up to 25 millimeter-sized particles in three dimensions, could be used to create 3D displays consisting of levitating physical voxels.
Human induced pluripotent stem cell-derived medial ganglionic eminence cells alleviate chronic epilepsy in rats, raising the possibility of autologous grafts, which could preclude the need for immune suppression after grafting and promote long-term graft–host integration.
By the end of the 21st century, annual mean residential electricity consumption in China’s Yangtze River Delta might increase by an estimated 9.2% for each 1 °C increase in global mean surface temperature relative to the late 20th century.
An analysis of journal articles, focusing on 11 writing-style components, suggests that good writing can noticeably enhance article impact.
Image courtesy of Tullio Rossi (artist).
