Effect of linkage disequilibrium on the identification of functional variants

Heritability of complex traits in sub-populations experiencing bottlenecks and growth

Populations that have experienced a bottleneck are regularly used in Genome Wide Association Studies (GWAS) to investigate variants associated with complex traits. It is generally understood that these isolated sub-populations may experience high frequency of otherwise rare variants with large effect size, and therefore provide a unique opportunity to study said trait. However, the demographic history of the population under investigation affects all SNPs that determine the complex trait genome-wide, changing its heritability and genetic architecture. We use a simulation based approach to identify the impact of the demographic processes of drift, expansion, and migration on the heritability of complex trait. We show that demography has considerable impact on complex traits. We then investigate the power to resolve heritability of complex traits in GWAS studies subjected to demographic effects. We find that demography is an important component for interpreting inference of complex traits and has a nuanced impact on the power of GWAS. We conclude that demographic histories need to be explicitly modelled to properly quantify the history of selection on a complex trait.

Joint linkage and association mapping of complex traits in shrub willow (Salix purpurea L.)

BACKGROUND AND AIMS

Increasing energy demands and the necessity to reduce greenhouse gas emissions are key motivating factors driving the development of lignocellulosic crops as an alternative to non-renewable energy sources. The effects of global climate change will require a better understanding of the genetic basis of complex adaptive traits to breed more resilient bioenergy feedstocks, like willow (Salix spp.). Shrub willow is a sustainable and dedicated bioenergy crop, bred to be fast-growing and high-yielding on marginal land without competing with food crops. In a rapidly changing climate, genomic advances will be vital for the sustained improvement of willow and other non-model bioenergy crops. Here, joint genetic mapping was used to exploit genetic variation garnered from both recent and historical recombination events in S. purpurea.

METHODS

A panel of North American naturalized S. purpurea accessions and full-sib F2S. purpurea population were genotyped and phenotyped for a suite of morphological, physiological, pest and disease resistance, and wood chemical composition traits, collected from multi-environment and multi-year replicated field trials. Controlling for population stratification and kinship in the association panel and spatial variation in the F2, a comprehensive mixed model analysis was used to dissect the complex genetic architecture and plasticity of these important traits.

KEY RESULTS

Individually, genome-wide association (GWAS) models differed in terms of power, but the combined approach, which corrects for yearly and environmental co-factors across datasets, improved the overall detection and resolution of associated loci. Although there were few significant GWAS hits located within support intervals of QTL for corresponding traits in the F2, many large-effect QTL were identified, as well as QTL hotspots.

CONCLUSIONS

This study provides the first comparison of linkage analysis and linkage disequilibrium mapping approaches in Salix, and highlights the complementarity and limits of these two methods for elucidating the genetic architecture of complex bioenergy-related traits of a woody perennial breeding programme.

Population-based analysis of ocular Chlamydia trachomatis in trachoma-endemic West African communities identifies genomic markers of disease severity

BACKGROUND

Chlamydia trachomatis (Ct) is the most common infectious cause of blindness and bacterial sexually transmitted infection worldwide. Ct strain-specific differences in clinical trachoma suggest that genetic polymorphisms in Ct may contribute to the observed variability in severity of clinical disease.

METHODS

Using Ct whole genome sequences obtained directly from conjunctival swabs, we studied Ct genomic diversity and associations between Ct genetic polymorphisms with ocular localization and disease severity in a treatment-naïve trachoma-endemic population in Guinea-Bissau, West Africa.

RESULTS

All Ct sequences fall within the T2 ocular clade phylogenetically. This is consistent with the presence of the characteristic deletion in trpA resulting in a truncated non-functional protein and the ocular tyrosine repeat regions present in tarP associated with ocular tissue localization. We have identified 21 Ct non-synonymous single nucleotide polymorphisms (SNPs) associated with ocular localization, including SNPs within pmpD (odds ratio, OR = 4.07, p* = 0.001) and tarP (OR = 0.34, p* = 0.009). Eight synonymous SNPs associated with disease severity were found in yjfH (rlmB) (OR = 0.13, p* = 0.037), CTA0273 (OR = 0.12, p* = 0.027), trmD (OR = 0.12, p* = 0.032), CTA0744 (OR = 0.12, p* = 0.041), glgA (OR = 0.10, p* = 0.026), alaS (OR = 0.10, p* = 0.032), pmpE (OR = 0.08, p* = 0.001) and the intergenic region CTA0744-CTA0745 (OR = 0.13, p* = 0.043).

CONCLUSIONS

This study demonstrates the extent of genomic diversity within a naturally circulating population of ocular Ct and is the first to describe novel genomic associations with disease severity. These findings direct investigation of host-pathogen interactions that may be important in ocular Ct pathogenesis and disease transmission.

The value of extended pedigrees for next-generation analysis of complex disease in the rhesus macaque

Complex diseases (e.g., cardiovascular disease and type 2 diabetes, among many others) pose the biggest threat to human health worldwide and are among the most challenging to investigate. Susceptibility to complex disease may be caused by multiple genetic variants (GVs) and their interaction, by environmental factors, and by interaction between GVs and environment, and large study cohorts with substantial analytical power are typically required to elucidate these individual contributions. Here, we discuss the advantages of both power and feasibility afforded by the use of extended pedigrees of rhesus macaques (Macaca mulatta) for genetic studies of complex human disease based on next-generation sequence data. We present these advantages in the context of previous research conducted in rhesus macaques for several representative complex diseases. We also describe a single, multigeneration pedigree of Indian-origin rhesus macaques and a sample biobank we have developed for genetic analysis of complex disease, including power of this pedigree to detect causal GVs using either genetic linkage or association methods in a variance decomposition approach. Finally, we summarize findings of significant heritability for a number of quantitative traits that demonstrate that genetic contributions to risk factors for complex disease can be detected and measured in this pedigree. We conclude that the development and application of an extended pedigree to analysis of complex disease traits in the rhesus macaque have shown promising early success and that genome-wide genetic and higher order -omics studies in this pedigree are likely to yield useful insights into the architecture of complex human disease.

Genome-wide association mapping for wood characteristics in Populus identifies an array of candidate single nucleotide polymorphisms

Establishing links between phenotypes and molecular variants is of central importance to accelerate genetic improvement of economically important plant species. Our work represents the first genome-wide association study to the inherently complex and currently poorly understood genetic architecture of industrially relevant wood traits. Here, we employed an Illumina Infinium 34K single nucleotide polymorphism (SNP) genotyping array that generated 29,233 high-quality SNPs in c. 3500 broad-based candidate genes within a population of 334 unrelated Populus trichocarpa individuals to establish genome-wide associations. The analysis revealed 141 significant SNPs (α ≤ 0.05) associated with 16 wood chemistry/ultrastructure traits, individually explaining 3-7% of the phenotypic variance. A large set of associations (41% of all hits) occurred in candidate genes preselected for their suggested a priori involvement with secondary growth. For example, an allelic variant in the FRA8 ortholog explained 21% of the total genetic variance in fiber length, when the trait's heritability estimate was considered. The remaining associations identified SNPs in genes not previously implicated in wood or secondary wall formation. Our findings provide unique insights into wood trait architecture and support efforts for population improvement based on desirable allelic variants.

Lessons learned from Genetic Analysis Workshop 17: transitioning from genome-wide association studies to whole-genome statistical genetic analysis

Genetic Analysis Workshop 17 (GAW17) focused on the transition from genome-wide association study designs and methods to the study designs and statistical genetic methods that will be required for the analysis of next-generation sequence data including both common and rare sequence variants. In the 166 contributions to GAW17, a wide variety of statistical methods were applied to simulated traits in population- and family-based samples, and results from these analyses were compared to the known generating model. In general, many of the statistical genetic methods used in the population-based sample identified causal sequence variants (SVs) when the estimated locus-specific heritability, as measured in the population-based sample, was greater than about 0.08. However, SVs with locus-specific heritabilities less than 0.03 were rarely identified consistently. In the family-based samples, many of the methods detected SVs that were rarer than those detected in the population-based sample, but the estimated locus-specific heritabilities for these rare SVs, as measured in the family-based samples, were substantially higher (>0.2) than their corresponding heritabilities in the population-based samples. Substantial inflation of the type I error rate was observed across a wide variety of statistical methods. Although many of the contributions found little inflation in type I error for Q4, a trait with no causal SVs, type I error rates for Q1 and Q2 were well above their nominal levels with the inflation for Q1 being higher than that for Q2. It seems likely that this inflation in type I error is due to correlations among SVs.