Heritability in the genome-wide association era

Heritability, the fraction of phenotypic variation explained by genetic variation, has been estimated for many phenotypes in a range of populations, organisms, and time points. The recent development of efficient genotyping and sequencing technology has led researchers to attempt to identify the genetic variants responsible for the genetic component of phenotype directly via GWAS. The gap between the phenotypic variance explained by GWAS results and those estimated from classical heritability methods has been termed the "missing heritability problem". In this work, we examine modern methods for estimating heritability, which use the genotype and sequence data directly. We discuss them in the context of classical heritability methods, the missing heritability problem, and describe their implications for understanding the genetic architecture of complex phenotypes.     

使用多因素易患性阈值模型评估双相情感障碍的遗传度缺失

为了评估双向情感障碍的遗传度缺失,文章通过查询美国国家人类基因组研究所(National Human Genome Research Institute,NHGRI)的gwascatalog目录,检索出所有已发现的双相情感障碍易感变异,使用多因素易患性阈值模型计算每个易感变异对双相情感障碍遗传度的解释度。将所有易感变异遗传度解释度求和得到双相情感障碍已知易感变异对遗传度的总解释度,使用此总解释度评估双相情感障碍的遗传度缺失。结果显示,已知双相情感障碍易感变异对双相情感障碍遗传度的合计解释度为38.34%,尚有61.66%的遗传度无法被已有易感变异解释,属于遗传度缺失。双相情感障碍38.34%的遗传度解释度较早前国外同类研究大幅度提高,表明随着新的双相情感障碍易感变异被不断发现,双相情感障碍遗传度缺失得到大幅度减小。但双相情感障碍遗传度缺失依然存在且数目较大的事实也表明双相情感障碍尚存在许多未知的分子遗传学机制有待进一步阐明。     

Environment-sensitive epigenetics and the heritability of complex diseases

Genome-wide association studies have thus far failed to explain the observed heritability of complex human diseases. This is referred to as the "missing heritability" problem. However, these analyses have usually neglected to consider a role for epigenetic variation, which has been associated with many human diseases. We extend models of epigenetic inheritance to investigate whether environment-sensitive epigenetic modifications of DNA might explain observed patterns of familial aggregation. We find that variation in epigenetic state and environmental state can result in highly heritable phenotypes through a combination of epigenetic and environmental inheritance. These two inheritance processes together can produce familial covariances significantly higher than those predicted by models of purely epigenetic inheritance and similar to those expected from genetic effects. The results suggest that epigenetic variation, inherited both directly and through shared environmental effects, may make a key contribution to the missing heritability.     

同卵双胞胎在复杂性状DNA甲基化调控机制研究中的应用

同卵双胞胎来源于同一个受精卵,DNA序列基本一致,但在某些重要表型上如复杂疾病,并不完全一样。利用表型不一致的同卵双胞胎进行研究,能在遗传背景、母体效应、年龄性别效应等一致的基础上,深入研究分析复杂性状的表观调控机制。而DNA甲基化是最为稳定的一类表观遗传修饰。在人类中,利用同卵双胞胎对印记异常疾病、精神类疾病、自身免疫病及癌症等疾病的DNA甲基化调控研究已经揭示了多个致病基因,为研究疾病的表观调控以及表观遗传学药物的应用打下了基础。本文着重对同卵双胞胎DNA甲基化状态、DNA甲基化遗传力计算以及复杂性状DNA甲基化调控的研究应用及其进展展开综述,以期为复杂性状表观调控机制研究提供借鉴和参考。     

The Nature of Genetic Variation for Complex Traits Revealed by GWAS and Regional Heritability Mapping Analyses

We use computer simulations to investigate the amount of genetic variation for complex traits that can be revealed by single-SNP genome-wide association studies (GWAS) or regional heritability mapping (RHM) analyses based on full genome sequence data or SNP chips. We model a large population subject to mutation, recombination, selection, and drift, assuming a pleiotropic model of mutations sampled from a bivariate distribution of effects of mutations on a quantitative trait and fitness. The pleiotropic model investigated, in contrast to previous models, implies that common mutations of large effect are responsible for most of the genetic variation for quantitative traits, except when the trait is fitness itself. We show that GWAS applied to the full sequence increases the number of QTL detected by as much as 50% compared to the number found with SNP chips but only modestly increases the amount of additive genetic variance explained. Even with full sequence data, the total amount of additive variance explained is generally below 50%. Using RHM on the full sequence data, a slightly larger number of QTL are detected than by GWAS if the same probability threshold is assumed, but these QTL explain a slightly smaller amount of genetic variance. Our results also suggest that most of the missing heritability is due to the inability to detect variants of moderate effect (∼0.03–0.3 phenotypic SDs) segregating at substantial frequencies. Very rare variants, which are more difficult to detect by GWAS, are expected to contribute little genetic variation, so their eventual detection is less relevant for resolving the missing heritability problem.     

Laboratory evolution of adenylyl cyclase independent learning in Drosophila and missing heritability

Gene interactions are acknowledged to be a likely source of missing heritability in large‐scale genetic studies of complex neurological phenotypes. However, involvement of rare variants, de novo mutations, genetic lesions that are not easily detected with commonly used methods and epigenetic factors also are possible explanations. We used a laboratory evolution study to investigate the modulatory effects of background genetic variation on the phenotypic effect size of a null mutation with known impact on olfactory learning. To accomplish this, we first established a population that contained variation at just 23 loci and used selection to evolve suppression of the learning defect seen with null mutations in the rutabaga adenylyl cyclase. We thus biased the system to favor relatively simplified outcomes by choosing a Mendelian trait and by restricting the genetic variation segregating in the population. This experimental design also assures that the causal effects are among the known 23 segregating loci. We observe a robust response to selection that requires the presence of the 23 variants. Analyses of the underlying genotypes showed that interactions between more than two loci are likely to be involved in explaining the selection response, with implications for the missing heritability problem.     

Why genetic selection to reduce the prevalence of infectious diseases is way more promising than currently believed

Genetic selection for improved disease resistance is an important part of strategies to combat infectious diseases in agriculture. Quantitative genetic analyses of binary disease status, however, indicate low heritability for most diseases, which restricts the rate of genetic reduction in disease prevalence. Moreover, the common liability threshold model suggests that eradication of an infectious disease via genetic selection is impossible because the observed-scale heritability goes to zero when the prevalence approaches zero. From infectious disease epidemiology, however, we know that eradication of infectious diseases is possible, both in theory and practice, because of positive feedback mechanisms leading to the phenomenon known as herd immunity. The common quantitative genetic models, however, ignore these feedback mechanisms. Here, we integrate quantitative genetic analysis of binary disease status with epidemiological models of transmission, aiming to identify the potential response to selection for reducing the prevalence of endemic infectious diseases. The results show that typical heritability values of binary disease status correspond to a very substantial genetic variation in disease susceptibility among individuals. Moreover, our results show that eradication of infectious diseases by genetic selection is possible in principle. These findings strongly disagree with predictions based on common quantitative genetic models, which ignore the positive feedback effects that occur when reducing the transmission of infectious diseases. Those feedback effects are a specific kind of Indirect Genetic Effects; they contribute substantially to the response to selection and the development of herd immunity (i.e., an effective reproduction ratio less than one).     

A general model of the relation between phenotypic selection and genetic response

The phenotypic view of selection assumes that genetic responses can be predicted from selective forces and heritability — or in the classical quantitative genetic equation: R = h²S. However, data on selection in bird populations show that often no selection responses is found, despite consistent selective forces on phenotypes and significant heritable variation. Such discrepancies may arise due to the assumption that selection only acts on observed phenotypes. We derive a general selection equation that takes into account the possibility that some relevant (internal or external) traits are not measured. This equation shows that the classic equation applies if selection directly acts on the measured, phenotypic traits. This is not the case when, for instance, there are unknown internal genetic trade-offs, or unknown common environmental factors affecting both trait and fitness. In such cases, any relationship between phenotypic selection and genetic response is possible. Fortunately, the classical model can be tested by comparing phenotypic and genetic covariances between traits and fitness; an indication that important internal or external traits are missing can thus be obtained. Such an analysis was indeed found in the literature; for selection on fledging weight in Great Tits it yielded valuable extra information.     

Imaging genomics: data fusion in uncovering disease heritability

Sequencing of the human genome in the early 2000s enabled probing of the genetic basis of disease on a scale previously unimaginable. Now, two decades later, after interrogating millions of markers in thousands of individuals, a significant portion of disease heritability still remains hidden. Recent efforts to unravel this ‘missing heritability’ have focused on garnering new insight from merging different data types, including medical imaging. Imaging offers promising intermediate phenotypes to bridge the gap between genetic variation and disease pathology. In this review we outline this fusion and provide examples of imaging genomics in a range of diseases, from oncology to cardiovascular and neurodegenerative disease. Finally, we discuss how ongoing revolutions in data science and sharing are primed to advance the field.     

Selection,structure and the heritability of behaviour

Characters which are closely linked to fitness often have low heritabilities (V_A/V_P). Low heritabilities could be because of low additive genetic variation (V_A), that had been depleted by directional selection. Alternatively, low heritabilities may be caused by large residual variation (V_R=V_P – V_A) compounded at a disproportionately higher rate than V_A across integrated characters. Both hypotheses assume that each component of quantitative variation has an independent effect on heritability. However, V_A and V_R may also covary, in which case differences in heritability cannot be fully explained by the independent effects of elimination‐selection or compounded residual variation. We compared the central tendency of published behavioural heritabilities (mean=0.31, median=0.23) with morphological and life history data collected by 26 ). Average behavioural heritability was not significantly different from average life history heritability, but both were smaller than average morphological heritability. We cross‐classified behavioural traits to test whether variation in heritability was related to selection (dominance, domestic/wild) or variance compounding (integration level). There was a significant three‐way interaction between indices of selection and variance compounding, related to the absence of either effect at the highest integration level. At lower integration levels, high dominance variance indicated effects of selection. It was also indicated by the low CV_A of domestic species. At the same time CV_R increased disproportionately faster than CV_A across integration levels, demonstrating variance compounding. However, neither CV_R nor CV_A had a predominant effect on heritability. The partial regression coefficients of CV_R and CV_A on heritability were similar and a path analysis indicated that their (positive) correlation was also necessary to explain variation in heritability. These results suggest that relationships between additive genetic and residual components of quantitative genetic variation can constrain their independent direct effects on behavioural heritability.     

茄子果实相关农艺性状全基因组关联分析研究进展与展望

茄子是重要的园艺作物,也是茄科植物中种植最广泛的蔬菜之一。茄子果实相关农艺性状是一种复杂的数量性状,传统育种选育效率低、周期长。高通量测序技术与生物信息学技术的快速发展,使得全基因组关联分析(genome-wide association study, GWAS)在解析茄子果实相关复杂农艺性状的遗传规律方面展现出巨大的应用前景。本文对全基因组关联分析在茄子的果形、果色等果实相关农艺性状中的研究进展进行了综述;针对茄子数量性状遗传研究中普遍存在的“丢失遗传力”(missing heritability)问题,从4个GWAS策略在茄子果实相关农艺性状研究中的应用热点出发,提出了未来茄子GWAS的发展对策;并结合当前茄子遗传改良的实践需求,展望了GWAS策略在茄子分子育种领域的广阔应用前景。本文为今后利用GWAS解析各种茄子果实相关性状的遗传基础以及选育符合消费者需求的果实材料提供了理论依据和参考。     

Heritability of interpack aggression in a wild pedigreed population of North American grey wolves

Aggression is a quantitative trait deeply entwined with individual fitness. Mapping the genomic architecture underlying such traits is complicated by complex inheritance patterns, social structure, pedigree information and gene pleiotropy. Here, we leveraged the pedigree of a reintroduced population of grey wolves (Canis lupus) in Yellowstone National Park, Wyoming, USA, to examine the heritability of and the genetic variation associated with aggression. Since their reintroduction, many ecological and behavioural aspects have been documented, providing unmatched records of aggressive behaviour across multiple generations of a wild population of wolves. Using a linear mixed model, a robust genetic relationship matrix, 12,288 single nucleotide polymorphisms (SNPs) and 111 wolves, we estimated the SNP‐based heritability of aggression to be 37% and an additional 14% of the phenotypic variation explained by shared environmental exposures. We identified 598 SNP genotypes from 425 grey wolves to resolve a consensus pedigree that was included in a heritability analysis of 141 individuals with SNP genotype, metadata and aggression data. The pedigree‐based heritability estimate for aggression is 14%, and an additional 16% of the phenotypic variation was explained by shared environmental exposures. We find strong effects of breeding status and relative pack size on aggression. Through an integrative approach, these results provide a framework for understanding the genetic architecture of a complex trait that influences individual fitness, with linkages to reproduction, in a social carnivore. Along with a few other studies, we show here the incredible utility of a pedigreed natural population for dissecting a complex, fitness‐related behavioural trait.     

Statistical Power to Detect Genetic (Co)Variance of Complex Traits Using SNP Data in Unrelated Samples

We have recently developed analysis methods (GREML) to estimate the genetic variance of a complex trait/disease and the genetic correlation between two complex traits/diseases using genome-wide single nucleotide polymorphism (SNP) data in unrelated individuals. Here we use analytical derivations and simulations to quantify the sampling variance of the estimate of the proportion of phenotypic variance captured by all SNPs for quantitative traits and case-control studies. We also derive the approximate sampling variance of the estimate of a genetic correlation in a bivariate analysis, when two complex traits are either measured on the same or different individuals. We show that the sampling variance is inversely proportional to the number of pairwise contrasts in the analysis and to the variance in SNP-derived genetic relationships. For bivariate analysis, the sampling variance of the genetic correlation additionally depends on the harmonic mean of the proportion of variance explained by the SNPs for the two traits and the genetic correlation between the traits, and depends on the phenotypic correlation when the traits are measured on the same individuals. We provide an online tool for calculating the power of detecting genetic (co)variation using genome-wide SNP data. The new theory and online tool will be helpful to plan experimental designs to estimate the missing heritability that has not yet been fully revealed through genome-wide association studies, and to estimate the genetic overlap between complex traits (diseases) in particular when the traits (diseases) are not measured on the same samples.     

Progress and promise of genome-wide association studies for human complex trait genetics

Enormous progress in mapping complex traits in humans has been made in the last 5 yr. There has been early success for prevalent diseases with complex phenotypes. These studies have demonstrated clearly that, while complex traits differ in their underlying genetic architectures, for many common disorders the predominant pattern is that of many loci, individually with small effects on phenotype. For some traits, loci of large effect have been identified. For almost all complex traits studied in humans, the sum of the identified genetic effects comprises only a portion, generally less than half, of the estimated trait heritability. A variety of hypotheses have been proposed to explain why this might be the case, including untested rare variants, and gene-gene and gene-environment interaction. Effort is currently being directed toward implementation of novel analytic approaches and testing rare variants for association with complex traits using imputed variants from the publicly available 1000 Genomes Project resequencing data and from direct resequencing of clinical samples. Through integration with annotations and functional genomic data as well as by in vitro and in vivo experimentation, mapping studies continue to characterize functional variants associated with complex traits and address fundamental issues such as epistasis and pleiotropy. This review focuses primarily on the ways in which genome-wide association studies (GWASs) have revolutionized the field of human quantitative genetics.     

An assessment of gene-by-gene interactions as a tool to unfold missing heritability in dyslexia

Even if substantial heritability has been reported and candidate genes have been identified extensively, all known marker associations explain only a small proportion of the phenotypic variance of developmental dyslexia (DD) and related quantitative phenotypes. Gene-by-gene interaction (also known as “epistasis”—G × G) triggers a non-additive effect of genes at different loci and should be taken into account in explaining part of the missing heritability of this complex trait. We assessed potential G × G interactions among five DD candidate genes, i.e., DYX1C1, DCDC2, KIAA0319, ROBO1, and GRIN2B, upon DD-related neuropsychological phenotypes in 493 nuclear families with DD, by implementing two complementary regression-based approaches: (1) a general linear model equation whereby the trait is predicted by the main effect of the number of rare alleles of the two genes and by the effect of the interaction between them, and (2) a family-based association test to detect G × G interactions between two unlinked markers by splitting up the association effect into a between- and a within-family genetic orthogonal components. After applying 500,000 permutations and correcting for multiple testing, both methods show that G × G effects between markers within the DYX1C1, KIAA0319/TTRAP, and GRIN2B genes lower the memory letters composite z-score of on average 0.55 standard deviation. We provided initial evidence that the effects of familial transmission of synergistic interactions between genetic risk variants can be exploited in the study of the etiology of DD, explain part of its missing heritability, and assist in designing customized charts of individualized neurocognitive impairments in complex disorders, such as DD.     

A GWAS assessment of the contribution of genomic imprinting to the variation of body mass index in mice

Background

Genomic imprinting is an epigenetic mechanism that can lead to differential gene expression depending on the parent-of-origin of a received allele. While most studies on imprinting address its underlying molecular mechanisms or attempt at discovering genomic regions that might be subject to imprinting, few have focused on the amount of phenotypic variation contributed by such epigenetic process. In this report, we give a brief review of a one-locus imprinting model in a quantitative genetics framework, and provide a decomposition of the genetic variance according to this model. Analytical deductions from the proposed imprinting model indicated a non-negligible contribution of imprinting to genetic variation of complex traits. Also, we performed a whole-genome scan analysis on mouse body mass index (BMI) aiming at revealing potential consequences when existing imprinting effects are ignored in genetic analysis.

Results

10,021 SNP markers were used to perform a whole-genome single marker regression on mouse BMI using an additive and an imprinting model. Markers significant for imprinting indicated that BMI is subject to imprinting. Marked variance changed from 1.218 ×10⁻⁴ to 1.842 ×10⁻⁴ when imprinting was considered in the analysis, implying that one third of marked variance would be lost if existing imprinting effects were not accounted for. When both marker and pedigree information were used, estimated heritability increased from 0.176 to 0.195 when imprinting was considered.

Conclusions

When a complex trait is subject to imprinting, using an additive model that ignores this phenomenon may result in an underestimate of additive variability, potentially leading to wrong inferences about the underlying genetic architecture of that trait. This could be a possible factor explaining part of the missing heritability commonly observed in genome-wide association studies (GWAS).     

De novo mutations in human genetic disease

New mutations have long been known to cause genetic disease, but their true contribution to the disease burden can only now be determined using family-based whole-genome or whole-exome sequencing approaches. In this Review we discuss recent findings suggesting that de novo mutations play a prominent part in rare and common forms of neurodevelopmental diseases, including intellectual disability, autism and schizophrenia. De novo mutations provide a mechanism by which early-onset reproductively lethal diseases remain frequent in the population. These mutations, although individually rare, may capture a significant part of the heritability for complex genetic diseases that is not detectable by genome-wide association studies.     

Simultaneous Discovery,Estimation and Prediction Analysis of Complex Traits Using a Bayesian Mixture Model

Gene discovery, estimation of heritability captured by SNP arrays, inference on genetic architecture and prediction analyses of complex traits are usually performed using different statistical models and methods, leading to inefficiency and loss of power. Here we use a Bayesian mixture model that simultaneously allows variant discovery, estimation of genetic variance explained by all variants and prediction of unobserved phenotypes in new samples. We apply the method to simulated data of quantitative traits and Welcome Trust Case Control Consortium (WTCCC) data on disease and show that it provides accurate estimates of SNP-based heritability, produces unbiased estimators of risk in new samples, and that it can estimate genetic architecture by partitioning variation across hundreds to thousands of SNPs. We estimated that, depending on the trait, 2,633 to 9,411 SNPs explain all of the SNP-based heritability in the WTCCC diseases. The majority of those SNPs (>96%) had small effects, confirming a substantial polygenic component to common diseases. The proportion of the SNP-based variance explained by large effects (each SNP explaining 1% of the variance) varied markedly between diseases, ranging from almost zero for bipolar disorder to 72% for type 1 diabetes. Prediction analyses demonstrate that for diseases with major loci, such as type 1 diabetes and rheumatoid arthritis, Bayesian methods outperform profile scoring or mixed model approaches.     

An exact form of the breeder's equation for the evolution of a quantitative trait under natural selection

Starting with the Price equation, I show that the total evolutionary change in mean phenotype that occurs in the presence of fitness variation can be partitioned exactly into five components representing logically distinct processes. One component is the linear response to selection, as represented by the breeder's equation of quantitative genetics, but with heritability defined as the linear regression coefficient of mean offspring phenotype on parent phenotype. The other components are identified as constitutive transmission bias, two types of induced transmission bias, and a spurious response to selection caused by a covariance between parental fitness and offspring phenotype that cannot be predicted from parental phenotypes. The partitioning can be accomplished in two ways, one with heritability measured before (in the absence of) selection, and the other with heritability measured after (in the presence of) selection. Measuring heritability after selection, though unconventional, yields a representation for the linear response to selection that is most consistent with Darwinian evolution by natural selection because the response to selection is determined by the reproductive features of the selected group, not of the parent population as a whole. The analysis of an explicitly Mendelian model shows that the relative contributions of the five terms to the total evolutionary change depends on the level of organization (gene, individual, or mated pair) at which the parent population is divided into phenotypes, with each frame of reference providing unique insight. It is shown that all five components of phenotypic evolution will generally have nonzero values as a result of various combinations of the normal features of Mendelian populations, including biparental sex, allelic dominance, inbreeding, epistasis, linkage disequilibrium, and environmental covariances between traits. Additive genetic variance can be a poor predictor of the adaptive response to selection in these models. The narrow-sense heritability sigma2A/sigma2P should be viewed as an approximation to the offspring-parent linear regression rather than the other way around.     

Heritability of age at menarche in girls from the Fels Longitudinal Study

Menarche is the hallmark maturational event of female childhood. Many studies indicated a significant genetic contribution to the timing of the onset of menstruation, but most of these studies were limited by the use of retrospective data and by the use of data from only certain types of relatives (i.e., mothers and daughters, sisters, or twin sisters). The primary goal of this study was to use a modern maximum likelihood quantitative genetic method to estimate the heritability (h(2)) of age at menarche, using familial data collected over the course of the 74-year-old Fels Longitudinal Study. The secondary goal was to review earlier studies of the heritability of age at menarche. The study of the heritability of age at menarche presented here is unique for two reasons. First, because of the Fels Longitudinal Study's serial design, age-at-menarche data were collected prospectively from most participants. Second, because the Fels Longitudinal Study is a family study that has been conducted for decades, age-at-menarche data are available from many types of female relatives spanning multiple households and generations. The best-fitting and most parsimonious quantitative genetic model included provision for a secular decrease in age at menarche, and estimated the h(2) of age at menarche to be 0.49+/- 0.13 (95% confidence interval of h(2),=0.24-0.73). The results of this study are in general agreement with the findings of most previous studies of genetic influences on age at menarche, and suggest that it is reasonable to consider it well-established that approximately half the phenotypic variation among girls from developed nations in the timing of menarche is due to genetic factors.