R/qtl2: Software for Mapping Quantitative Trait Loci with High-Dimensional Data and Multiparent Populations

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.

Natural variation in Arabidopsis shoot branching plasticity in response to nitrate supply affects fitness

The capacity of organisms to tune their development in response to environmental cues is pervasive in nature. This phenotypic plasticity is particularly striking in plants, enabled by their modular and continuous development. A good example is the activation of lateral shoot branches in Arabidopsis, which develop from axillary meristems at the base of leaves. The activity and elongation of lateral shoots depends on the integration of many signals both external (e.g. light, nutrient supply) and internal (e.g. the phytohormones auxin, strigolactone and cytokinin). Here, we characterise natural variation in plasticity of shoot branching in response to nitrate supply using two diverse panels of Arabidopsis lines. We find extensive variation in nitrate sensitivity across these lines, suggesting a genetic basis for variation in branching plasticity. High plasticity is associated with extreme branching phenotypes such that lines with the most branches on high nitrate have the fewest under nitrate deficient conditions. Conversely, low plasticity is associated with a constitutively moderate level of branching. Furthermore, variation in plasticity is associated with alternative life histories with the low plasticity lines flowering significantly earlier than high plasticity lines. In Arabidopsis, branching is highly correlated with fruit yield, and thus low plasticity lines produce more fruit than high plasticity lines under nitrate deficient conditions, whereas highly plastic lines produce more fruit under high nitrate conditions. Low and high plasticity, associated with early and late flowering respectively, can therefore be interpreted alternative escape vs mitigate strategies to low N environments. The genetic architecture of these traits appears to be highly complex, with only a small proportion of the estimated genetic variance detected in association mapping.

Genetic determinants of gut microbiota composition and bile acid profiles in mice

The microbial communities that inhabit the distal gut of humans and other mammals exhibit large inter-individual variation. While host genetics is a known factor that influences gut microbiota composition, the mechanisms underlying this variation remain largely unknown. Bile acids (BAs) are hormones that are produced by the host and chemically modified by gut bacteria. BAs serve as environmental cues and nutrients to microbes, but they can also have antibacterial effects. We hypothesized that host genetic variation in BA metabolism and homeostasis influence gut microbiota composition. To address this, we used the Diversity Outbred (DO) stock, a population of genetically distinct mice derived from eight founder strains. We characterized the fecal microbiota composition and plasma and cecal BA profiles from 400 DO mice maintained on a high-fat high-sucrose diet for ~22 weeks. Using quantitative trait locus (QTL) analysis, we identified several genomic regions associated with variations in both bacterial and BA profiles. Notably, we found overlapping QTL for Turicibacter sp. and plasma cholic acid, which mapped to a locus containing the gene for the ileal bile acid transporter, Slc10a2. Mediation analysis and subsequent follow-up validation experiments suggest that differences in Slc10a2 gene expression associated with the different strains influences levels of both traits and revealed novel interactions between Turicibacter and BAs. This work illustrates how systems genetics can be utilized to generate testable hypotheses and provide insight into host-microbe interactions.

Construction of Genetic Linkage Maps in Multiparental Populations

Construction of genetic linkage maps has become a routine step for mapping quantitative trait loci (QTL), particularly in animal and plant breeding populations. Many multiparental populations have recently been produced to increase genetic diversity and QTL mapping resolution. However, few software packages are available for map construction in these populations. In this paper, we build a general framework for the construction of genetic linkage maps from genotypic data in diploid populations, including bi- and multiparental populations, cross-pollinated (CP) populations, and breeding pedigrees. The framework is implemented as an automatic pipeline called magicMap, where the maximum multilocus likelihood approach utilizes genotypic information efficiently. We evaluate magicMap by extensive simulations and eight real datasets: one biparental, one CP, four multiparent advanced generation intercross (MAGIC), and two nested association mapping (NAM) populations, the number of markers ranging from a few hundred to tens of thousands. Not only is magicMap the only software capable of accommodating all of these designs, it is more accurate and robust to missing genotypes and genotyping errors than commonly used packages.

Gene loci associated with insulin secretion in islets from non-diabetic mice

Genetic susceptibility to type 2 diabetes is primarily due to β-cell dysfunction. However, a genetic study to directly interrogate β-cell function ex vivo has never been previously performed. We isolated 233,447 islets from 483 Diversity Outbred (DO) mice maintained on a Western-style diet, and measured insulin secretion in response to a variety of secretagogues. Insulin secretion from DO islets ranged >1,000-fold even though none of the mice were diabetic. The insulin secretory response to each secretagogue had a unique genetic architecture; some of the loci were specific for one condition, whereas others overlapped. Human loci that are syntenic to many of the insulin secretion QTL from mouse are associated with diabetes-related SNPs in human genome-wide association studies. We report on three genes, Ptpn18, Hunk and Zfp148, where the phenotype predictions from the genetic screen were fulfilled in our studies of transgenic mouse models. These three genes encode a non-receptor type protein tyrosine phosphatase, a serine/threonine protein kinase, and a Krϋppel-type zinc-finger transcription factor, respectively. Our results demonstrate that genetic variation in insulin secretion that can lead to type 2 diabetes is discoverable in non-diabetic individuals.

Testing Pleiotropy vs. Separate QTL in Multiparental Populations

The high mapping resolution of multiparental populations, combined with technology to measure tens of thousands of phenotypes, presents a need for quantitative methods to enhance understanding of the genetic architecture of complex traits. When multiple traits map to a common genomic region, knowledge of the number of distinct loci provides important insight into the underlying mechanism and can assist planning for subsequent experiments. We extend the method of Jiang and Zeng (1995), for testing pleiotropy with a pair of traits, to the case of more than two alleles. We also incorporate polygenic random effects to account for population structure. We use a parametric bootstrap to determine statistical significance. We apply our methods to a behavioral genetics data set from Diversity Outbred mice. Our methods have been incorporated into the R package qtl2pleio.

High-Diversity Mouse Populations for Complex Traits

Contemporary mouse genetic reference populations are a powerful platform to discover complex disease mechanisms. Advanced high-diversity mouse populations include the Collaborative Cross (CC) strains, Diversity Outbred (DO) stock, and their isogenic founder strains. When used in systems genetics and integrative genomics analyses, these populations efficiently harnesses known genetic variation for precise and contextualized identification of complex disease mechanisms. Extensive genetic, genomic, and phenotypic data are already available for these high-diversity mouse populations and a growing suite of data analysis tools have been developed to support research on diverse mice. This integrated resource can be used to discover and evaluate disease mechanisms relevant across species.

Determinants of QTL Mapping Power in the Realized Collaborative Cross

The Collaborative Cross (CC) is a mouse genetic reference population whose range of applications includes quantitative trait loci (QTL) mapping. The design of a CC QTL mapping study involves multiple decisions, including which and how many strains to use, and how many replicates per strain to phenotype, all viewed within the context of hypothesized QTL architecture. Until now, these decisions have been informed largely by early power analyses that were based on simulated, hypothetical CC genomes. Now that more than 50 CC strains are available and more than 70 CC genomes have been observed, it is possible to characterize power based on realized CC genomes. We report power analyses from extensive simulations and examine several key considerations: 1) the number of strains and biological replicates, 2) the QTL effect size, 3) the presence of population structure, and 4) the distribution of functionally distinct alleles among the founder strains at the QTL. We also provide general power estimates to aide in the design of future experiments. All analyses were conducted with our R package, SPARCC (Simulated Power Analysis in the Realized Collaborative Cross), developed for performing either large scale power analyses or those tailored to particular CC experiments.

Genetic Dissection of Resistance to the Three Fungal Plant Pathogens Blumeria graminis, Zymoseptoria tritici, and Pyrenophora tritici-repentis Using a Multiparental Winter Wheat Population

Bread wheat (Triticum aestivum L.) is one of the world’s most important crop species. The development of new varieties resistant to multiple pathogens is an ongoing task in wheat breeding, especially in times of increasing demand for sustainable agricultural practices. Despite this, little is known about the relations between various fungal disease resistances at the genetic level, and the possible consequences for wheat breeding strategies. As a first step to fill this gap, we analyzed the genetic relations of resistance to the three fungal diseases – powdery mildew (PM), septoria tritici blotch (STB), and tan spot (TS) – using a winter wheat multiparent advanced generation intercross population. Six, seven, and nine QTL for resistance to PM, STB, and TS, respectively, were genetically mapped. Additionally, 15 QTL were identified for the three agro-morphological traits plant height, ear emergence time, and leaf angle distribution. Our results suggest that resistance to STB and TS on chromosome 2B is conferred by the same genetic region. Furthermore, we identified two genetic regions on chromosome 1AS and 7AL, which are associated with all three diseases, but not always in a synchronal manner. Based on our results, we conclude that parallel marker-assisted breeding for resistance to the fungal diseases PM, STB, and TS appears feasible. Knowledge of the genetic co-localization of alleles with contrasting effects for different diseases, such as on chromosome 7AL, allows the trade-offs of selection of these regions to be better understood, and ultimately determined at the genic level.

Cleaning Genotype Data from Diversity Outbred Mice

Data cleaning is an important first step in most statistical analyses, including efforts to map the genetic loci that contribute to variation in quantitative traits. Here we illustrate approaches to quality control and cleaning of array-based genotyping data for multiparent populations (experimental crosses derived from more than two founder strains), using MegaMUGA array data from a set of 291 Diversity Outbred (DO) mice. Our approach employs data visualizations that can reveal problems at the level of individual mice or with individual SNP markers. We find that the proportion of missing genotypes for each mouse is an effective indicator of sample quality. We use microarray probe intensities for SNPs on the X and Y chromosomes to confirm the sex of each mouse, and we use the proportion of matching SNP genotypes between pairs of mice to detect sample duplicates. We use a hidden Markov model (HMM) reconstruction of the founder haplotype mosaic across each mouse genome to estimate the number of crossovers and to identify potential genotyping errors. To evaluate marker quality, we find that missing data and genotyping error rates are the most effective diagnostics. We also examine the SNP genotype frequencies with markers grouped according to their minor allele frequency in the founder strains. For markers with high apparent error rates, a scatterplot of the allele-specific probe intensities can reveal the underlying cause of incorrect genotype calls. The decision to include or exclude low-quality samples can have a significant impact on the mapping results for a given study. We find that the impact of low-quality markers on a given study is often minimal, but reporting problematic markers can improve the utility of the genotyping array across many studies.