Missing heritability problem

Last updated

The missing heritability problem [1] [2] [3] [4] [5] [6] is the fact that single genetic variations cannot account for much of the heritability of diseases, behaviors, and other phenotypes. This is a problem that has significant implications for medicine, since a person's susceptibility to disease may depend more on the combined effect of all the genes in the background than on the disease genes in the foreground, or the role of genes may have been severely overestimated.

Contents

Discovery

The missing heritability problem was named as such in 2008 (after the "missing baryon problem" in physics). The Human Genome Project led to optimistic forecasts that the large genetic contributions to many traits and diseases (which were identified by quantitative genetics and behavioral genetics in particular) would soon be mapped and pinned down to specific genes and their genetic variants by methods such as candidate-gene studies which used small samples with limited genetic sequencing to focus on specific genes believed to be involved, examining single-nucleotide polymorphisms (SNPs). While many hits were found, they often failed to replicate in other studies.

The exponential fall in genome sequencing costs led to the use of genome-wide association studies (GWASes) which could simultaneously examine all candidate-genes in larger samples than the original finding, where the candidate-gene hits were found to almost always be false positives and only 2-6% replicate; [7] [8] [9] [10] [11] [12] in the specific case of intelligence candidate-gene hits, only 1 candidate-gene hit replicated, [13] the top 25 schizophrenia candidate-genes were no more associated with schizophrenia than chance, [14] [15] and of 15 neuroimaging hits, none did. [16] In 2012, the editorial board of Behavior Genetics noted, in setting more stringent requirements for candidate-gene publications, that "the literature on candidate gene associations is full of reports that have not stood up to rigorous replication...it now seems likely that many of the published findings of the last decade are wrong or misleading and have not contributed to real advances in knowledge". [17] Other researchers have characterized the literature as having "yielded an infinitude of publications with very few consistent replications" and called for a phase out of candidate-gene studies in favor of polygenic scores. [18]

Dilemma

Standard genetics methods have long estimated large heritabilities such as 80% for traits such as height or intelligence, yet none of the genes had been found despite sample sizes that, while small, should have been able to detect variants of reasonable effect size such as 1 inch or 5 IQ points. If genes have such strong cumulative effects - where were they? Several resolutions have been proposed, that the missing heritability is some combination of:

  1. Twin studies and other methods were grossly biased by issues long raised by their critics; there was little genetic influence to be found. Therefore, it has been proposed that the genes that supposedly underlie behavior genetic estimates of heritability simply do not exist. [19]
  2. Genetic effects are actually epigenetics
  3. Genetic effects are generally non-additive and due to complex interactions. Among many proposals, a model has been introduced that takes into account epigenetic inheritance on the risk and recurrence risk of a complex disease. [4] The limiting pathway (LP) model has been introduced in which a trait depends on the value of k inputs that can have rate limitations due to stoichiometric ratios, reactants required in a biochemical pathway, or proteins required for transcription of a gene. Each of these k inputs is a strictly additive trait that depends on a set of common or rare variants. When k = 1, the LP model is simply a standard additive trait. [2]
  4. Genetic effects are not due to the common SNPs examined in the candidate-gene studies & GWASes, but due to very rare mutations, copy-number variations, and other exotic kinds of genetic variants. These variants tend to be harmful and kept at low frequencies by natural selection. Whole-genome sequencing would be required to track down specific rare variants.
  5. Traits are all misdiagnoses: one person's 'schizophrenia' is due to entirely different causes than another schizophrenic, and so while there may be a gene involved in one case, it will not be involved in another, rendering GWASes futile
  6. GWASes are unable to detect genes with moderate effects on phenotypes when those genes segregate at high frequencies [20]
  7. Traits are genuine but inconsistently diagnosed or genetically influenced from country to country and time to time, leading to measurement error, which combined with genetic heterogeneity, either due to race or environment, will bias meta-analyzed GWAS & GCTA results towards zero, [21] [22] [23] [24] [25] [26]
  8. Genetic effects are indeed through common SNPs acting additively, but are highly polygenic: dispersed over hundreds or thousands of variants each of small effect like a fraction of an inch or a fifth of an IQ point and with low prior probability: unexpected enough that a candidate-gene study is unlikely to select the right SNP out of hundreds of thousands of known SNPs, and GWASes up to 2010, with n<20000, would be unable to find hits which reach genome-wide statistical-significance thresholds. Much larger GWAS sample sizes, often n>100k, would be required to find any hits at all, and would steadily increase after that.
This resolution to the missing heritability problem was supported by the introduction of Genome-wide complex trait analysis (GCTA) in 2010, which demonstrated that trait similarity could be predicted by the genetic similarity of unrelated strangers on common SNPs treated additively, and for many traits the SNP heritability was indeed a substantial fraction of the overall heritability. The GCTA results were further supported by findings that a small percent of trait variance could be predicted in GWASes without any genome-wide statistically-significant hits by a linear model including all SNPs regardless of p-value; if there were no SNP contribution, this would be unlikely, but it would be what one expected from SNPs whose effects were very imprecisely estimated by a too-small sample. Combined with the upper bound on maximum effect sizes set by the GWASes up to then, this strongly implied that the highly polygenic theory was correct. Examples of complex traits where increasingly large-scale GWASes have yielded the initial hits and then increasing numbers of hits as sample sizes increased from n<20k to n>100k or n>300k include height, [27] educational attainment, [28] and schizophrenia.

Related Research Articles

<span class="mw-page-title-main">Heritability</span> Estimation of effect of genetic variation on phenotypic variation of a trait

Heritability is a statistic used in the fields of breeding and genetics that estimates the degree of variation in a phenotypic trait in a population that is due to genetic variation between individuals in that population. The concept of heritability can be expressed in the form of the following question: "What is the proportion of the variation in a given trait within a population that is not explained by the environment or random chance?"

<span class="mw-page-title-main">Single-nucleotide polymorphism</span> Single nucleotide in genomic DNA at which different sequence alternatives exist

In genetics and bioinformatics, a single-nucleotide polymorphism is a germline substitution of a single nucleotide at a specific position in the genome that is present in a sufficiently large fraction of considered population.

The candidate gene approach to conducting genetic association studies focuses on associations between genetic variation within pre-specified genes of interest, and phenotypes or disease states. This is in contrast to genome-wide association studies (GWAS), which is a hypothesis-free approach that scans the entire genome for associations between common genetic variants and traits of interest. Candidate genes are most often selected for study based on a priori knowledge of the gene's biological functional impact on the trait or disease in question. The rationale behind focusing on allelic variation in specific, biologically relevant regions of the genome is that certain alleles within a gene may directly impact the function of the gene in question and lead to variation in the phenotype or disease state being investigated. This approach often uses the case-control study design to try to answer the question, "Is one allele of a candidate gene more frequently seen in subjects with the disease than in subjects without the disease?" Candidate genes hypothesized to be associated with complex traits have generally not been replicated by subsequent GWASs or highly powered replication attempts. The failure of candidate gene studies to shed light on the specific genes underlying such traits has been ascribed to insufficient statistical power, low prior probability that scientists can correctly guess a specific allele within a specific gene that is related to a trait, poor methodological practices, and data dredging.

Peter McGuffin was a psychiatrist and geneticist from Belfast, Northern Ireland.

<span class="mw-page-title-main">Genome-wide association study</span> Study of genetic variants in different individuals

In genomics, a genome-wide association study, is an observational study of a genome-wide set of genetic variants in different individuals to see if any variant is associated with a trait. GWA studies typically focus on associations between single-nucleotide polymorphisms (SNPs) and traits like major human diseases, but can equally be applied to any other genetic variants and any other organisms.

In multivariate quantitative genetics, a genetic correlation is the proportion of variance that two traits share due to genetic causes, the correlation between the genetic influences on a trait and the genetic influences on a different trait estimating the degree of pleiotropy or causal overlap. A genetic correlation of 0 implies that the genetic effects on one trait are independent of the other, while a correlation of 1 implies that all of the genetic influences on the two traits are identical. The bivariate genetic correlation can be generalized to inferring genetic latent variable factors across > 2 traits using factor analysis. Genetic correlation models were introduced into behavioral genetics in the 1970s–1980s.

Behavioural genetics, also referred to as behaviour genetics, is a field of scientific research that uses genetic methods to investigate the nature and origins of individual differences in behaviour. While the name "behavioural genetics" connotes a focus on genetic influences, the field broadly investigates the extent to which genetic and environmental factors influence individual differences, and the development of research designs that can remove the confounding of genes and environment. Behavioural genetics was founded as a scientific discipline by Francis Galton in the late 19th century, only to be discredited through association with eugenics movements before and during World War II. In the latter half of the 20th century, the field saw renewed prominence with research on inheritance of behaviour and mental illness in humans, as well as research on genetically informative model organisms through selective breeding and crosses. In the late 20th and early 21st centuries, technological advances in molecular genetics made it possible to measure and modify the genome directly. This led to major advances in model organism research and in human studies, leading to new scientific discoveries.

Expression quantitative trait loci (eQTLs) are genomic loci that explain variation in expression levels of mRNAs.

<span class="mw-page-title-main">Phosphatase and actin regulator 1</span> Protein-coding gene in the species Homo sapiens

Phosphatase and actin regulator 1 (PHACTR1) is a protein that in humans is encoded by the PHACTR1 gene on chromosome 6. It is most significantly expressed in the globus pallidus of the brain. PHACTR1 is an actin and protein phosphatase 1 (PP1) binding protein that binds actin and regulates the reorganization of the actin cytoskeleton. This protein has been associated with coronary artery disease and migraines through genome-wide association studies. The PHACTR1 gene also contains one of 27 SNPs associated with increased risk of coronary artery disease.

Predictive genomics is at the intersection of multiple disciplines: predictive medicine, personal genomics and translational bioinformatics. Specifically, predictive genomics deals with the future phenotypic outcomes via prediction in areas such as complex multifactorial diseases in humans. To date, the success of predictive genomics has been dependent on the genetic framework underlying these applications, typically explored in genome-wide association (GWA) studies. The identification of associated single-nucleotide polymorphisms underpin GWA studies in complex diseases that have ranged from Type 2 Diabetes (T2D), Age-related macular degeneration (AMD) and Crohn's disease.

Genome-wide complex trait analysis (GCTA) Genome-based restricted maximum likelihood (GREML) is a statistical method for heritability estimation in genetics, which quantifies the total additive contribution of a set of genetic variants to a trait. GCTA is typically applied to common single nucleotide polymorphisms (SNPs) on a genotyping array and thus termed "chip" or "SNP" heritability.

A human disease modifier gene is a modifier gene that alters expression of a human gene at another locus that in turn causes a genetic disease. Whereas medical genetics has tended to distinguish between monogenic traits, governed by simple, Mendelian inheritance, and quantitative traits, with cumulative, multifactorial causes, increasing evidence suggests that human diseases exist on a continuous spectrum between the two.

<span class="mw-page-title-main">Polygenic score</span> Numerical score aimed at predicting a trait based on variation in multiple genetic loci

In genetics, a polygenic score (PGS) is a number that summarizes the estimated effect of many genetic variants on an individual's phenotype. The PGS is also called the polygenic index (PGI) or genome-wide score; in the context of disease risk, it is called a polygenic risk score or genetic risk score. The score reflects an individual's estimated genetic predisposition for a given trait and can be used as a predictor for that trait. It gives an estimate of how likely an individual is to have a given trait based only on genetics, without taking environmental factors into account; and it is typically calculated as a weighted sum of trait-associated alleles.

<span class="mw-page-title-main">Complex traits</span>

Complex traits, also known as quantitative traits, are traits that do not behave according to simple Mendelian inheritance laws. More specifically, their inheritance cannot be explained by the genetic segregation of a single gene. Such traits show a continuous range of variation and are influenced by both environmental and genetic factors. Compared to strictly Mendelian traits, complex traits are far more common, and because they can be hugely polygenic, they are studied using statistical techniques such as quantitative genetics and quantitative trait loci (QTL) mapping rather than classical genetics methods. Examples of complex traits include height, circadian rhythms, enzyme kinetics, and many diseases including diabetes and Parkinson's disease. One major goal of genetic research today is to better understand the molecular mechanisms through which genetic variants act to influence complex traits.

In statistical genetics, linkage disequilibrium score regression is a technique that aims to quantify the separate contributions of polygenic effects and various confounding factors, such as population stratification, based on summary statistics from genome-wide association studies (GWASs). The approach involves using regression analysis to examine the relationship between linkage disequilibrium scores and the test statistics of the single-nucleotide polymorphisms (SNPs) from the GWAS. Here, the "linkage disequilibrium score" for a SNP "is the sum of LD r2 measured with all other SNPs".

The Omnigenic Model, first proposed by Evan A. Boyle, Yang I. Li, and Jonathan K. Pritchard, describes a hypothesis regarding the heritability of complex traits. Expanding beyond polygenes, the authors propose that all genes expressed within a cell affect the expression of a given trait. In addition, the model states that the peripheral genes, ones that do not have a direct impact on expression, explain more heritability of traits than core genes, ones that have a direct impact on expression. The process that the authors propose that facilitates this effect is called “network pleiotropy”, in which peripheral genes can affect core genes, not by having a direct effect, but rather by virtue of being mediated within the same cell.

The GWAS catalog is a free online database that compiles data of genome-wide association studies (GWAS), summarizing unstructured data from different literature sources into accessible high quality data. It was created by the National Human Genome Research Institute (NHGRI) in 2008 and have become a collaborative project between the NHGRI and the European Bioinformatics Institute (EBI) since 2010. As of September 2018, it has included 71,673 SNP–trait associations in 3,567 publications.

Personality traits are patterns of thoughts, feelings and behaviors that reflect the tendency to respond in certain ways under certain circumstances.

Transcriptome-wide association study (TWAS) is a genetic methodology that can be used to compare the genetic components of gene expression and the genetic components of a trait to determine if an association is present between the two components. TWAS are useful for the identification and prioritization of candidate causal genes in candidate gene analysis following genome-wide association studies. TWAS looks at the RNA products of a specific tissue and gives researchers the abilities to look at the genes being expressed as well as gene expression levels, which varies by tissue type. TWAS are valuable and flexible bioinformatics tools that looks at the associations between the expressions of genes and complex traits and diseases. By looking at the association between gene expression and the trait expressed, genetic regulatory mechanisms can be investigated for the role that they play in the development of specific traits and diseases.

<span class="mw-page-title-main">Andre Franke</span> German geneticist

Andre Franke, born on 16 October 1978, is a geneticist, academic, and university professor. He is a Full W3 Professor of Molecular Medicine at the Christian-Albrechts-University of Kiel, and a managing director at the Institute of Clinical Molecular Biology.

References

  1. Manolio, T. A.; Collins, F. S.; Cox, N. J.; Goldstein, D. B.; Hindorff, L. A.; Hunter, D. J.; McCarthy, M. I.; Ramos, E. M.; Cardon, L. R.; Chakravarti, A.; Cho, J. H.; Guttmacher, A. E.; Kong, A.; Kruglyak, L.; Mardis, E.; Rotimi, C. N.; Slatkin, M.; Valle, D.; Whittemore, A. S.; Boehnke, M.; Clark, A. G.; Eichler, E. E.; Gibson, G.; Haines, J. L.; MacKay, T. F. C.; McCarroll, S. A.; Visscher, P. M. (2009). "Finding the missing heritability of complex diseases". Nature . 461 (7265): 747–753. Bibcode:2009Natur.461..747M. doi:10.1038/nature08494. PMC   2831613 . PMID   19812666.
  2. 1 2 Zuk, O.; Hechter, E.; Sunyaev, S. R.; Lander, E. S. (2012). "The mystery of missing heritability: Genetic interactions create phantom heritability". Proceedings of the National Academy of Sciences . 109 (4): 1193–1198. Bibcode:2012PNAS..109.1193Z. doi: 10.1073/pnas.1119675109 . PMC   3268279 . PMID   22223662.
  3. Lee, S. H.; Wray, N. R.; Goddard, M. E.; Visscher, P. M. (2011). "Estimating Missing Heritability for Disease from Genome-wide Association Studies". American Journal of Human Genetics . 88 (3): 294–305. doi:10.1016/j.ajhg.2011.02.002. PMC   3059431 . PMID   21376301.
  4. 1 2 Slatkin, M. (2009). "Epigenetic Inheritance and the Missing Heritability Problem". Genetics . 182 (3): 845–850. doi:10.1534/genetics.109.102798. PMC   2710163 . PMID   19416939.
  5. Eichler, E. E.; Flint, J.; Gibson, G.; Kong, A.; Leal, S. M.; Moore, J. H.; Nadeau, J. H. (2010). "Missing heritability and strategies for finding the underlying causes of complex disease". Nature Reviews Genetics . 11 (6): 446–450. doi:10.1038/nrg2809. PMC   2942068 . PMID   20479774.
  6. "Personal genomes: The case of the missing heritability", Maher 2008, Nature.
  7. "Replication Validity of Initial Association Studies: A Comparison between Psychiatry, Neurology and Four Somatic Diseases", Dumas-Mallet et al 2016
  8. "The False-positive to False-negative Ratio in Epidemiologic Studies", Ioannidis et al 2011
  9. "A Test-Replicate Approach to Candidate Gene Research on Addiction and Externalizing Disorders: A Collaboration Across Five Longitudinal Studies", Samek et al 2016
  10. Bevan et al 2012, "Genetic heritability of ischemic stroke and the contribution of previously reported candidate gene and genome-wide associations"
  11. Siontis et al 2010, "Replication of past candidate loci for common diseases and phenotypes in 100 genome-wide association studies"
  12. Duncan & Keller 2011, "A Critical Review of the First 10 Years of Candidate Gene-by-Environment Interaction Research in Psychiatry"
  13. Chabris, CF; Hebert, BM; Benjamin, DJ; Beauchamp, J; Cesarini, D; van der Loos, M; Johannesson, M; Magnusson, PK; Lichtenstein, P; Atwood, CS; Freese, J; Hauser, TS; Hauser, RM; Christakis, N; Laibson, D (2012). "Most reported genetic associations with general intelligence are probably false positives". Psychol Sci . 23 (11): 1314–23. doi:10.1177/0956797611435528. PMC   3498585 . PMID   23012269.
  14. Johnson et al 2017, "No evidence that schizophrenia candidate genes are more associated with schizophrenia than non-candidate genes"
  15. "Reproducibility in Imaging Genetics: The Case of Threat-Related Amygdala Reactivity", Avinun et al 2017
  16. Jahanshad et al 2017, "Do Candidate Genes Affect the Brain's White Matter Microstructure? Large-Scale Evaluation of 6,165 Diffusion MRI Scans"
  17. Hewitt 2012, "Editorial Policy on Candidate Gene Association and Candidate Gene-by-Environment Interaction Studies of Complex traits"
  18. Arango 2017, "Candidate gene associations studies in psychiatry: time to move forward"
  19. Chaufan, Claudia; Joseph, Jay (April 2013). "The 'Missing Heritability' of Common Disorders: Should Health Researchers Care?". International Journal of Health Services. 43 (2): 281–303. doi:10.2190/hs.43.2.f. ISSN   0020-7314. PMID   23821906. S2CID   25092977.
  20. Caballero, Armando; Tenesa, Albert; Keightley, Peter D. (December 2015). "The Nature of Genetic Variation for Complex Traits Revealed by GWAS and Regional Heritability Mapping Analyses". Genetics. 201 (4): 1601–1613. doi:10.1534/genetics.115.177220. ISSN   1943-2631. PMC   4676519 . PMID   26482794.
  21. "Meta-GWAS Accuracy and Power (MetaGAP) calculator shows that hiding heritability is partially due to imperfect genetic correlations across studies", de Vlaming et al 2016: MetaGAP
  22. Wray & Maier 2014, "Genetic basis of complex genetic disease: the contribution of disease heterogeneity to missing heritability"
  23. Wray et al 2012, "Impact of diagnostic misclassification on estimation of genetic correlations using genome-wide genotypes"
  24. Lee et al 2013a, "Genetic relationship between five psychiatric disorders estimated from genome-wide SNPs"
  25. Lee et al 2013b, "General framework for meta-analysis of rare variants in sequencing association studies"
  26. Sham & Purcell 2014, "Statistical power and significance testing in large-scale genetic studies"
  27. "Defining the role of common variation in the genomic and biological architecture of adult human height", Wood et al 2014
  28. Chabris et al 2012 reported only 1 possible hit using a few thousand; "GWAS of 126,559 Individuals Identifies Genetic Variants Associated with Educational Attainment", Rietveld et al 2013 with n=100k reported 3 hits; "Genome-wide association study identifies 74 loci associated with educational attainment", Okbay et al 2016 reported 74 hits using n=293k and ~160 when extended to n=404k