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.     

Towards identifying genes underlying ecologically relevant traits in Arabidopsis thaliana

A major challenge in evolutionary biology and plant breeding is to identify the genetic basis of complex quantitative traits, including those that contribute to adaptive variation. Here we review the development of new methods and resources to fine-map intraspecific genetic variation that underlies natural phenotypic variation in plants. In particular, the analysis of 107 quantitative traits reported in the first genome-wide association mapping study in Arabidopsis thaliana sets the stage for an exciting time in our understanding of plant adaptation. We also argue for the need to place phenotype-genotype association studies in an ecological context if one is to predict the evolutionary trajectories of plant species.     

Association genetics of complex traits in plants

Association mapping is rapidly becoming the main method for dissecting the genetic architecture of complex traits in plants. Currently most association mapping studies in plants are preformed using sets of genes selected to be putative candidates for the trait of interest, but rapid developments in genomics will allow for genome-wide mapping in virtually any plant species in the near future. As the costs for genotyping are decreasing, the focus has shifted towards phenotyping. In plants, clonal replication and/or inbred lines allows for replicated phenotyping under many different environmental conditions. Reduced sequencing costs will increase the number of studies that use RNA sequencing data to perform expression quantitative trait locus (eQTL) mapping, which will increase our knowledge of how gene expression variation contributes to phenotypic variation. Current population sizes used in association mapping studies are modest in size and need to be greatly increased if mutations explaining less than a few per cent of the phenotypic variation are to be detected. Association mapping has started to yield insights into the genetic architecture of complex traits in plants, and future studies with greater genome coverage will help to elucidate how plants have managed to adapt to a wide variety of environmental conditions.     

Simulation of Finnish Population History,Guided by Empirical Genetic Data,to Assess Power of Rare-Variant Tests in Finland

Finnish samples have been extensively utilized in studying single-gene disorders, where the founder effect has clearly aided in discovery, and more recently in genome-wide association studies of complex traits, where the founder effect has had less obvious impacts. As the field starts to explore rare variants’ contribution to polygenic traits, it is of great importance to characterize and confirm the Finnish founder effect in sequencing data and to assess its implications for rare-variant association studies. Here, we employ forward simulation, guided by empirical deep resequencing data, to model the genetic architecture of quantitative polygenic traits in both the general European and the Finnish populations simultaneously. We demonstrate that power of rare-variant association tests is higher in the Finnish population, especially when variants’ phenotypic effects are tightly coupled with fitness effects and therefore reflect a greater contribution of rarer variants. SKAT-O, variable-threshold tests, and single-variant tests are more powerful than other rare-variant methods in the Finnish population across a range of genetic models. We also compare the relative power and efficiency of exome array genotyping to those of high-coverage exome sequencing. At a fixed cost, less expensive genotyping strategies have far greater power than sequencing; in a fixed number of samples, however, genotyping arrays miss a substantial portion of genetic signals detected in sequencing, even in the Finnish founder population. As genetic studies probe sequence variation at greater depth in more diverse populations, our simulation approach provides a framework for evaluating various study designs for gene discovery.     

Genome-wide prediction of functional gene-gene interactions inferred from patterns of genetic differentiation in mice and men

The human genome encodes a limited number of genes yet contributes to individual differences in a vast array of heritable traits. A possible explanation for the capacity our genome to generate this virtually unlimited range of phenotypic variation in complex traits is to assume functional interactions between genes. Therefore we searched two mammalian genomes to identify potential epistatic interactions by looking for co-adapted genes marked by excess two-locus genetic differentiation between populations/lineages using publicly available SNP genotype data. The practical motivation for this effort is to reduce the number of pair-wise tests that need to be performed in genome-wide association studies aimed at detecting GxG interactions, by focusing on pairs predicted to be more likely to jointly affect variation in complex traits. Hence, this approach generates a list of candidate interactions that can be empirically tested. In both the mouse and human data we observed two-locus genetic differentiation in excess of what can be expected from chance alone based on simulations. In an attempt to validate our hypothesis that pairs of genes showing excess genetic divergence represent potential functional interactions, we selected a small set of gene combinations postulated to be interacting based on our analyses and looked for a combined effect of the selected genes on variation in complex traits in both mice and man. In both cases the individual effect of the genes were not significant, instead we observed marginally significant interaction effects. These results show that genome wide searches for gene-gene interactions based on population genetic data are feasible and can generate interesting candidate gene pairs to be further tested for their contribution to phenotypic variation in complex traits.     

复杂疾病全基因组关联研究进展—— 研究设计和遗传标记

实现全基因组关联研究(Genome-wide association study, GWA)在数年前还是遗传学家们的梦想, 如今它已经变成了现实。自2005年Science杂志报道了第一项有关年龄相关性(视网膜)黄斑变性全基因组关联研究研究以来, 有关与复杂疾病的全基因组关联研究如雨后春笋般层出不穷。文中介绍了近两年来全基因组关联研究在复杂疾病研究领域内的主要发现、全基因组关联研究设计原理、遗传标记的选择、比较及相关商品信息。最后介绍了人类基因组拷贝数变异的研究进展, 总结了人类全基因组关联研究所取得成就和存在的问题, 并对全基因组关联研究未来的研究重点和要解决的问题进行了展望。     

Accounting for Population Structure in Gene-by-Environment Interactions in Genome-Wide Association Studies Using Mixed Models

Although genome-wide association studies (GWASs) have discovered numerous novel genetic variants associated with many complex traits and diseases, those genetic variants typically explain only a small fraction of phenotypic variance. Factors that account for phenotypic variance include environmental factors and gene-by-environment interactions (GEIs). Recently, several studies have conducted genome-wide gene-by-environment association analyses and demonstrated important roles of GEIs in complex traits. One of the main challenges in these association studies is to control effects of population structure that may cause spurious associations. Many studies have analyzed how population structure influences statistics of genetic variants and developed several statistical approaches to correct for population structure. However, the impact of population structure on GEI statistics in GWASs has not been extensively studied and nor have there been methods designed to correct for population structure on GEI statistics. In this paper, we show both analytically and empirically that population structure may cause spurious GEIs and use both simulation and two GWAS datasets to support our finding. We propose a statistical approach based on mixed models to account for population structure on GEI statistics. We find that our approach effectively controls population structure on statistics for GEIs as well as for genetic variants.     

Lessons learned from the dog genome

Extensive genetic resources and a high-quality genome sequence position the dog as an important model species for understanding genome evolution, population genetics and genes underlying complex phenotypic traits. Newly developed genomic resources have expanded our understanding of canine evolutionary history and dog origins. Domestication involved genetic contributions from multiple populations of gray wolves probably through backcrossing. More recently, the advent of controlled breeding practices has segregated genetic variability into distinct dog breeds that possess specific phenotypic traits. Consequently, genome-wide association and selective sweep scans now allow the discovery of genes underlying breed-specific characteristics. The dog is finally emerging as a novel resource for studying the genetic basis of complex traits, including behavior.     

On the Prospects of Whole-Genome Association Mapping in Saccharomyces cerevisiae

Advances in sequencing technology have enabled whole-genome sequences to be obtained from multiple individuals within species, particularly in model organisms with compact genomes. For example, 36 genome sequences of Saccharomyces cerevisiae are now publicly available, and SNP data are available for even larger collections of strains. One potential use of these resources is mapping the genetic basis of phenotypic variation through genome-wide association (GWA) studies, with the benefit that associated variants can be studied experimentally with greater ease than in outbred populations such as humans. Here, we evaluate the prospects of GWA studies in S. cerevisiae strains through extensive simulations and a GWA study of mitochondrial copy number. We demonstrate that the complex and heterogeneous patterns of population structure present in yeast populations can lead to a high type I error rate in GWA studies of quantitative traits, and that methods typically used to control for population stratification do not provide adequate control of the type I error rate. Moreover, we show that while GWA studies of quantitative traits in S. cerevisiae may be difficult depending on the particular set of strains studied, association studies to map cis-acting quantitative trait loci (QTL) and Mendelian phenotypes are more feasible. We also discuss sampling strategies that could enable GWA studies in yeast and illustrate the utility of this approach in Saccharomyces paradoxus. Thus, our results provide important practical insights into the design and interpretation of GWA studies in yeast, and other model organisms that possess complex patterns of population structure.     

Assumptions and Properties of Limiting Pathway Models for Analysis of Epistasis in Complex Traits

For most complex traits, results from genome-wide association studies show that the proportion of the phenotypic variance attributable to the additive effects of individual SNPs, that is, the heritability explained by the SNPs, is substantially less than the estimate of heritability obtained by standard methods using correlations between relatives. This difference has been called the “missing heritability”. One explanation is that heritability estimates from family (including twin) studies are biased upwards. Zuk et al. revisited overestimation of narrow sense heritability from twin studies as a result of confounding with non-additive genetic variance. They propose a limiting pathway (LP) model that generates significant epistatic variation and its simple parametrization provides a convenient way to explore implications of epistasis. They conclude that over-estimation of narrow sense heritability from family data (‘phantom heritability’) may explain an important proportion of missing heritability. We show that for highly heritable quantitative traits large phantom heritability estimates from twin studies are possible only if a large contribution of common environment is assumed. The LP model is underpinned by strong assumptions that are unlikely to hold, including that all contributing pathways have the same mean and variance and are uncorrelated. Here, we relax the assumptions that underlie the LP model to be more biologically plausible. Together with theoretical, empirical, and pragmatic arguments we conclude that in outbred populations the contribution of additive genetic variance is likely to be much more important than the contribution of non-additive variance.     

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.     

Tools for efficient epistasis detection in genome-wide association study

Background  

Genome-wide association study (GWAS) aims to find genetic factors underlying complex phenotypic traits, for which epistasis or gene-gene interaction detection is often preferred over single-locus approach. However, the computational burden has been a major hurdle to apply epistasis test in the genome-wide scale due to a large number of single nucleotide polymorphism (SNP) pairs to be tested.     

全基因组关联分析实现水稻粒型自然变异的分子解析

全基因组关联分析(GWAS)近年来被广泛应用于解析生物自然变异的遗传基础。但限于其遗传定位精度, 在水稻(Oryza sativa)遗传学研究中, 该方法尚无法取代传统的图位克隆法在克隆复杂性状调控基因中的作用。近期, 中国科学家在应用GWAS等大数据来克隆控制水稻粒长和粒重等复杂性状的QTL方面取得了新突破。     

Common garden experiments in the genomic era: new perspectives and opportunities

The study of local adaptation is rendered difficult by many evolutionary confounding phenomena (for example, genetic drift and demographic history). When complex traits are involved in local adaptation, phenomena such as phenotypic plasticity further hamper evolutionary biologists to study the complex relationships between phenotype, genotype and environment. In this perspective paper, we suggest that the common garden experiment, specifically designed to deal with phenotypic plasticity, has a clear role to play in the study of local adaptation, even (if not specifically) in the genomic era. After a quick review of some high-throughput genotyping protocols relevant in the context of a common garden, we explore how to improve common garden analyses with dense marker panel data and recent statistical methods. We then show how combining approaches from population genomics and genome-wide association studies with the settings of a common garden can yield to a very efficient, thorough and integrative study of local adaptation. Especially, evidence from genomic (for example, genome scan) and phenotypic origins constitute independent insights into the possibility of local adaptation scenarios, and genome-wide association studies in the context of a common garden experiment allow to decipher the genetic bases of adaptive traits.     

An Efficient Stepwise Statistical Test to Identify Multiple Linked Human Genetic Variants Associated with Specific Phenotypic Traits

Recent advances in genotyping methodologies have allowed genome-wide association studies (GWAS) to accurately identify genetic variants that associate with common or pathological complex traits. Although most GWAS have focused on associations with single genetic variants, joint identification of multiple genetic variants, and how they interact, is essential for understanding the genetic architecture of complex phenotypic traits. Here, we propose an efficient stepwise method based on the Cochran-Mantel-Haenszel test (for stratified categorical data) to identify causal joint multiple genetic variants in GWAS. This method combines the CMH statistic with a stepwise procedure to detect multiple genetic variants associated with specific categorical traits, using a series of associated I × J contingency tables and a null hypothesis of no phenotype association. Through a new stratification scheme based on the sum of minor allele count criteria, we make the method more feasible for GWAS data having sample sizes of several thousands. We also examine the properties of the proposed stepwise method via simulation studies, and show that the stepwise CMH test performs better than other existing methods (e.g., logistic regression and detection of associations by Markov blanket) for identifying multiple genetic variants. Finally, we apply the proposed approach to two genomic sequencing datasets to detect linked genetic variants associated with bipolar disorder and obesity, respectively.     

From QTL to QTN: Candidate Gene Set Approach and a Case Study in Porcine IGF1-FoxO Pathway

Unraveling the genetic background of economic traits is a major goal in modern animal genetics and breeding. Both candidate gene analysis and QTL mapping have previously been used for identifying genes and chromosome regions related to studied traits. However, most of these studies may be limited in their ability to fully consider how multiple genetic factors may influence a particular phenotype of interest. If possible, taking advantage of the combined effect of multiple genetic factors is expected to be more powerful than analyzing single sites, as the joint action of multiple loci within a gene or across multiple genes acting in the same gene set will likely have a greater influence on phenotypic variation. Thus, we proposed a pipeline of gene set analysis that utilized information from multiple loci to improve statistical power. We assessed the performance of this approach by both simulated and a real IGF1-FoxO pathway data set. The results showed that our new method can identify the association between genetic variation and phenotypic variation with higher statistical power and unravel the mechanisms of complex traits in a point of gene set. Additionally, the proposed pipeline is flexible to be extended to model complex genetic structures that include the interactions between different gene sets and between gene sets and environments.     

Genetic diversity and population structure of a diverse set of rice germplasm for association mapping

Germplasm diversity is the mainstay for crop improvement and genetic dissection of complex traits. Understanding genetic diversity, population structure, and the level and distribution of linkage disequilibrium (LD) in target populations is of great importance and a prerequisite for association mapping. In this study, 100 genome-wide simple sequence repeat (SSR) markers were used to assess genetic diversity, population structure, and LD of 416 rice accessions including landraces, cultivars and breeding lines collected mostly in China. A model-based population structure analysis divided the rice materials into seven subpopulations. 63% of the SSR pairs in these accessions were in LD, which was mostly due to an overall population structure, since the number of locus pairs in LD was reduced sharply within each subpopulation, with the SSR pairs in LD ranging from 5.9 to 22.9%. Among those SSR pairs showing significant LD, the intrachromosomal LD had an average of 25–50 cM in different subpopulations. Analysis of the phenotypic diversity of 25 traits showed that the population structure accounted for an average of 22.4% of phenotypic variation. An example association mapping for starch quality traits using both the candidate gene mapping and genome-wide mapping strategies based on the estimated population structure was conducted. Candidate gene mapping confirmed that the Wx and starch synthase IIa (SSIIa) genes could be identified as strongly associated with apparent amylose content (AAC) and pasting temperature (PT), respectively. More importantly, we revealed that the Wx gene was also strongly associated with PT. In addition to the major genes, we found five and seven SSRs were associated with AAC and PT, respectively, some of which have not been detected in previous linkage mapping studies. The results suggested that the population may be useful for the genome-wide marker–trait association mapping. This new association population has the potential to identify quantitative trait loci (QTL) with small effects, which will aid in dissecting complex traits and in exploiting the rich diversity present in rice germplasm.     

Searching for the regulators of human gene expression

Many common human traits are believed to be a composite reflection of multiple genetic and non-genetic factors and the genetic contribution is consequently often difficult to characterise. Recent advances suggest that subtle variation in the regulation of gene expression may contribute to complex human traits. In two reports, Cheung and colleagues scale up human genetics analysis to an impressive level in a genome-wide search for the regulators of gene expression. They perform linkage analysis on expression profiles for over 3,500 genes and then employ the HapMap resource to take positive findings through to association studies at the genome-wide level. This work gives new insights into the complexities of gene regulation and the plausibility of genome-wide study design.     

Classical phenotyping and deep learning concur on genetic control of stomatal density and area in sorghum

Stomatal density (SD) and stomatal complex area (SCA) are important traits that regulate gas exchange and abiotic stress response in plants. Despite sorghum (Sorghum bicolor) adaptation to arid conditions, the genetic potential of stomata-related traits remains unexplored due to challenges in available phenotyping methods. Hence, identifying loci that control stomatal traits is fundamental to designing strategies to breed sorghum with optimized stomatal regulation. We implemented both classical and deep learning methods to characterize genetic diversity in 311 grain sorghum accessions for stomatal traits at two different field environments. Nearly 12,000 images collected from abaxial (Ab) and adaxial (Ad) leaf surfaces revealed substantial variation in stomatal traits. Our study demonstrated significant accuracy between manual and deep learning methods in predicting SD and SCA. In sorghum, SD was 32%–39% greater on the Ab versus the Ad surface, while SCA on the Ab surface was 2%–5% smaller than on the Ad surface. Genome-Wide Association Study identified 71 genetic loci (38 were environment-specific) with significant genotype to phenotype associations for stomatal traits. Putative causal genes underlying the phenotypic variation were identified. Accessions with similar SCA but carrying contrasting haplotypes for SD were tested for stomatal conductance and carbon assimilation under field conditions. Our findings provide a foundation for further studies on the genetic and molecular mechanisms controlling stomata patterning and regulation in sorghum. An integrated physiological, deep learning, and genomic approach allowed us to unravel the genetic control of natural variation in stomata traits in sorghum, which can be applied to other plants.

High-throughput phenotyping using deep learning tools integrated with genome-wide association studies revealed genes that control SD and area in grain sorghum.