Whole-genome QTL analysis for MAGIC

Key message

An efficient whole genome method of QTL analysis is presented for Multi-parent advanced generation integrated crosses.

Abstract

Multi-parent advanced generation inter-cross (MAGIC) populations have been developed for mice and several plant species and are useful for the genetic dissection of complex traits. The analysis of quantitative trait loci (QTL) in these populations presents some additional challenges compared with traditional mapping approaches. In particular, pedigree and marker information need to be integrated and founder genetic data needs to be incorporated into the analysis. Here, we present a method for QTL analysis that utilizes the probability of inheriting founder alleles across the whole genome simultaneously, either for intervals or markers. The probabilities can be found using three-point or Hidden Markov Model (HMM) methods. This whole-genome approach is evaluated in a simulation study and it is shown to be a powerful method of analysis. The HMM probabilities lead to low rates of false positives and low bias of estimated QTL effect sizes. An implementation of the approach is available as an R package. In addition, we illustrate the approach using a bread wheat MAGIC population.     

Joint linkage and linkage disequilibrium mapping of quantitative trait loci in natural populations

Linkage analysis and allelic association (also referred to as linkage disequilibrium) studies are two major approaches for mapping genes that control simple or complex traits in plants, animals, and humans. But these two approaches have limited utility when used alone, because they use only part of the information that is available for a mapping population. More recently, a new mapping strategy has been designed to integrate the advantages of linkage analysis and linkage disequilibrium analysis for genome mapping in outcrossing populations. The new strategy makes use of a random sample from a panmictic population and the open-pollinated progeny of the sample. In this article, we extend the new strategy to map quantitative trait loci (QTL), using molecular markers within the EM-implemented maximum-likelihood framework. The most significant advantage of this extension is that both linkage and linkage disequilibrium between a marker and QTL can be estimated simultaneously, thus increasing the efficiency and effectiveness of genome mapping for recalcitrant outcrossing species. Simulation studies are performed to test the statistical properties of the MLEs of genetic and genomic parameters including QTL allele frequency, QTL effects, QTL position, and the linkage disequilibrium of the QTL and a marker. The potential utility of our mapping strategy is discussed.     

Quantitative trait loci mapping in five new large recombinant inbred line populations of Arabidopsis thaliana genotyped with consensus single-nucleotide polymorphism markers

Quantitative approaches conducted in a single mapping population are limited by the extent of genetic variation distinguishing the parental genotypes. To overcome this limitation and allow a more complete dissection of the genetic architecture of complex traits, we built an integrated set of 15 new large Arabidopsis thaliana recombinant inbred line (RIL) populations optimized for quantitative trait loci (QTL) mapping, having Columbia as a common parent crossed to distant accessions. Here we present 5 of these populations that were validated by investigating three traits: flowering time, rosette size, and seed production as an estimate of fitness. The large number of RILs in each population (between 319 and 377 lines) and the high density of evenly spaced genetic markers scored ensure high power and precision in QTL mapping even under a minimal phenotyping framework. Moreover, the use of common markers across the different maps allows a direct comparison of the QTL detected within the different RIL sets. In addition, we show that following a selective phenotyping strategy by performing QTL analyses on genotypically chosen subsets of 164 RILs (core populations) does not impair the power of detection of QTL with phenotypic contributions >7%.     

Potential of a tomato MAGIC population to decipher the genetic control of quantitative traits and detect causal variants in the resequencing era

Identification of the polymorphisms controlling quantitative traits remains a challenge for plant geneticists. Multiparent advanced generation intercross (MAGIC) populations offer an alternative to traditional linkage or association mapping populations by increasing the precision of quantitative trait loci (QTL) mapping. Here, we present the first tomato MAGIC population and highlight its potential for the valorization of intraspecific variation, QTL mapping and causal polymorphism identification. The population was developed by crossing eight founder lines, selected to include a wide range of genetic diversity, whose genomes have been previously resequenced. We selected 1536 SNPs among the 4 million available to enhance haplotype prediction and recombination detection in the population. The linkage map obtained showed an 87% increase in recombination frequencies compared to biparental populations. The prediction of the haplotype origin was possible for 89% of the MAGIC line genomes, allowing QTL detection at the haplotype level. We grew the population in two greenhouse trials and detected QTLs for fruit weight. We mapped three stable QTLs and six specific of a location. Finally, we showed the potential of the MAGIC population when coupled with whole genome sequencing of founder lines to detect candidate SNPs underlying the QTLs. For a previously cloned QTL on chromosome 3, we used the predicted allelic effect of each founder and their genome sequences to select putative causal polymorphisms in the supporting interval. The number of candidate polymorphisms was reduced from 12 284 (in 800 genes) to 96 (in 54 genes), including the actual causal polymorphism. This population represents a new permanent resource for the tomato genetics community.     

Quantitative trait locus study design from an information perspective

We examine the efficiency of different genotyping and phenotyping strategies in inbred line crosses from an information perspective. This provides a mathematical framework for the statistical aspects of QTL experimental design, while guiding our intuition. Our central result is a simple formula that quantifies the fraction of missing information of any genotyping strategy in a backcross. It includes the special case of selectively genotyping only the phenotypic extreme individuals. The formula is a function of the square of the phenotype and the uncertainty in our knowledge of the genotypes at a locus. This result is used to answer a variety of questions. First, we examine the cost-information trade-off varying the density of markers and the proportion of extreme phenotypic individuals genotyped. Then we evaluate the information content of selective phenotyping designs and the impact of measurement error in phenotyping. A simple formula quantifies the information content of any combined phenotyping and genotyping design. We extend our results to cover multigenotype crosses, such as the F(2) intercross, and multiple QTL models. We find that when the QTL effect is small, any contrast in a multigenotype cross benefits from selective genotyping in the same manner as in a backcross. The benefit remains in the presence of a second unlinked QTL with small effect (explaining <20% of the variance), but diminishes if the second QTL has a large effect. Software for performing power calculations for backcross and F(2) intercross incorporating selective genotyping and marker spacing is available from http://www.biostat.ucsf.edu/sen.     

A Multiparent Advanced Generation Inter-Cross to Fine-Map Quantitative Traits in Arabidopsis thaliana

Identifying natural allelic variation that underlies quantitative trait variation remains a fundamental problem in genetics. Most studies have employed either simple synthetic populations with restricted allelic variation or performed association mapping on a sample of naturally occurring haplotypes. Both of these approaches have some limitations, therefore alternative resources for the genetic dissection of complex traits continue to be sought. Here we describe one such alternative, the Multiparent Advanced Generation Inter-Cross (MAGIC). This approach is expected to improve the precision with which QTL can be mapped, improving the outlook for QTL cloning. Here, we present the first panel of MAGIC lines developed: a set of 527 recombinant inbred lines (RILs) descended from a heterogeneous stock of 19 intermated accessions of the plant Arabidopsis thaliana. These lines and the 19 founders were genotyped with 1,260 single nucleotide polymorphisms and phenotyped for development-related traits. Analytical methods were developed to fine-map quantitative trait loci (QTL) in the MAGIC lines by reconstructing the genome of each line as a mosaic of the founders. We show by simulation that QTL explaining 10% of the phenotypic variance will be detected in most situations with an average mapping error of about 300 kb, and that if the number of lines were doubled the mapping error would be under 200 kb. We also show how the power to detect a QTL and the mapping accuracy vary, depending on QTL location. We demonstrate the utility of this new mapping population by mapping several known QTL with high precision and by finding novel QTL for germination data and bolting time. Our results provide strong support for similar ongoing efforts to produce MAGIC lines in other organisms.     

Selective phenotyping for increased efficiency in genetic mapping studies

The power of a genetic mapping study depends on the heritability of the trait, the number of individuals included in the analysis, and the genetic dissimilarity among them. In experiments that involve microarrays or other complex physiological assays, phenotyping can be expensive and time-consuming and may impose limits on the sample size. A random selection of individuals may not provide sufficient power to detect linkage until a large sample size is reached. We present an algorithm for selecting a subset of individuals solely on the basis of genotype data that can achieve substantial improvements in sensitivity compared to a random sample of the same size. The selective phenotyping method involves preferentially selecting individuals to maximize their genotypic dissimilarity. Selective phenotyping is most effective when prior knowledge of genetic architecture allows us to focus on specific genetic regions. However, it can also provide modest improvements in efficiency when applied on a whole-genome basis. Importantly, selective phenotyping does not reduce the efficiency of mapping as compared to a random sample in regions that are not considered in the selection process. In contrast to selective genotyping, inferences based solely on a selectively phenotyped population of individuals are representative of the whole population. The substantial improvement introduced by selective phenotyping is particularly useful when phenotyping is difficult or costly and thus limits the sample size in a genetic mapping study.     

Effect of population size on the estimation of QTL: a test using resistance to barley stripe rust

The limited population sizes used in many quantitative trait locus (QTL) detection experiments can lead to underestimation of QTL number, overestimation of QTL effects, and failure to quantify QTL interactions. We used the barley/barley stripe rust pathosystem to evaluate the effect of population size on the estimation of QTL parameters. We generated a large (n=409) population of doubled haploid lines derived from the cross of two inbred lines, BCD47 and Baronesse. This population was evaluated for barley stripe rust severity in the Toluca Valley, Mexico, and in Washington State, USA, under field conditions. BCD47 was the principal donor of resistance QTL alleles, but the susceptible parent also contributed some resistance alleles. The major QTL, located on the long arm of chromosome 4H, close to the Mlo gene, accounted for up to 34% of the phenotypic variance. Subpopulations of different sizes were generated using three methods—resampling, selective genotyping, and selective phenotyping—to evaluate the effect of population size on the estimation of QTL parameters. In all cases, the number of QTL detected increased with population size. QTL with large effects were detected even in small populations, but QTL with small effects were detected only by increasing population size. Selective genotyping and/or selective phenotyping approaches could be effective strategies for reducing the costs associated with conducting QTL analysis in large populations. The method of choice will depend on the relative costs of genotyping versus phenotyping. Electronic Supplementary Material Supplementary material is available for this article at     

Sampling strategy for wild soybean (Glycine soja) populations based on their genetic diversity and fine-scale spatial genetic structure

A total of 892 individuals sampled from a wild soybean population in a natural reserve near the Yellow River estuary located in Kenli of Shandong Province (China) were investigated.Seventeen SSR (simple sequence repeat) primer pairs from cultivated soybeans were used to estimate the genetic diversity of the population and its variation pattern versus changes of the sample size (sub-samples),in addition to investigating the fine-scale spatial genetic structure within the population.The results showed relatively high genetic diversity of the population with the mean value of allele number (A) being 2.88,expected heterozygosity (He) 0.431,Shannon diversity index (/) 0.699,and percentage of polymorphic loci (P) 100%.Sub-samples of different sizes (ten groups) were randomly drawn from the population and their genetic diversity was calculated by computer simulation.The regression model of the four diversity indexes with the change of sample sizes was computed.As a result,27-52 individuals can reach 95% of total genetic variability of the population.Spatial autocorrelation analysis revealed that the genetic patch size of this wild soybean population is about 18 m.The study provided a scientific basis for the sampling strategy of wild soybean populations.     

Selection strategies for linkage studies using twins.

Genetic linkage analysis for complex diseases offers a major challenge to geneticists. In these complex diseases multiple genetic loci are responsible for the disease and they may vary in the size of their contribution; the effect of any single one of them is likely to be small. In many situations, like in extensive twin registries, trait values have been recorded for a large number of individuals, and preliminary studies have revealed summary measures for those traits, like mean, variance and components of variance, including heritability. Given the small effect size, a random sample of twins will require a prohibitively large sample size. It is well known that selective sampling is far more efficient in terms of genotyping effort. In this paper we derive easy expressions for the information contributed by sib pairs for the detection of linkage to a quantitative trait locus (QTL). We consider random samples as well as samples of sib pairs selected on the basis of their trait values. These expressions can be rapidly computed and do not involve simulation. We extend our results for quantitative traits to dichotomous traits using the concept of a liability threshold model. We present tables with required sample sizes for height, insulin levels and migraine, three of the traits studied in the GenomEUtwin project.     

Where should wild species be sampled? New method based on isolation‐by‐distance objectively gives the answer

Random sampling is an important statistical assumption, but virtually impossible when sampling a wild species as we cannot know where all the individuals exist. While interpopulation or intrataxa sampling methods have been developed, there are currently few intrataxon sampling methods to objectively decide where to sample wild taxa. We suggest a new sampling method which computes appropriate sampling locations from coordinates, assuming geographical autocorrelation of phylogeny within a taxon (isolation‐by‐distance). The computed locations encompass the highest genetic diversity, providing a genetically representative sample. In addition, it can utilize presence/absence information during sampling to reoptimize sampling scheme. Comparing to the single existing method of the similar purpose, the merits of ours is unnecessity of environmental data resulting in easy application, and is theoretically deduced. We tested this method using published phylogeographical data. The test result was generally encouraging, but the method failed where species showed uniform genetic structure or recent distribution expansion which violate the assumption of geographical autocorrelation of phylogeny. Though simple, our method constructs a methodological and statistical foundation for sampling wild species, and is applicable to revising taxonomic study and conservation biology.     

A mixed-model approach to mapping quantitative trait loci in barley on the basis of multiple environment data

In this article, I propose a mixed-model method to detect QTL with significant mean effect across environments and to characterize the stability of effects across multiple environments. I demonstrate the method using the barley dataset by the North American Barley Genome Mapping Project. The analysis raises the need for mixed modeling in two different ways. First, it is reasonable to regard environments as a random sample from a population of target environments. Thus, environmental main effects and QTL-by-environment interaction effects are regarded as random. Second, I expect a genetic correlation among pairs of environments caused by undetected QTL. I show how random QTL-by-environment effects as well as genetic correlations are straightforwardly handled in a mixed-model framework. The main advantage of this method is the ability to assess the stability of QTL effects. Moreover, the method allows valid statistical inferences regarding average QTL effects.     

Optimum allocation of resources for QTL detection using a nested association mapping strategy in maize

In quantitative trait locus (QTL) mapping studies, it is mandatory that the available financial resources are spent in such a way that the power for detection of QTL is maximized. The objective of this study was to optimize for three different fixed budgets the power of QTL detection 1 − β* in recombinant inbred line (RIL) populations derived from a nested design by varying (1) the genetic complexity of the trait, (2) the costs for developing, genotyping, and phenotyping RILs, (3) the total number of RILs, and (4) the number of environments and replications per environment used for phenotyping. Our computer simulations were based on empirical data of 653 single nucleotide polymorphism markers of 26 diverse maize inbred lines which were selected on the basis of 100 simple sequence repeat markers out of a worldwide sample of 260 maize inbreds to capture the maximum genetic diversity. For the standard scenario of costs, the optimum number of test environments (E _opt) ranged across the examined total budgets from 7 to 19 in the scenarios with 25 QTL. In comparison, the E _opt values observed for the scenarios with 50 and 100 QTL were slightly higher. Our finding of differences in 1 − β* estimates between experiments with optimally and sub-optimally allocated resources illustrated the potential to improve the power for QTL detection without increasing the total resources necessary for a QTL mapping experiment. Furthermore, the results of our study indicated that also in studies using the latest genomics tools to dissect quantitative traits, it is required to evaluate the individuals of the mapping population in a high number of environments with a high number of replications per environment.     

A multiparent advanced generation inter-cross population for genetic analysis in wheat

We present the first results from a novel multiparent advanced generation inter-cross (MAGIC) population derived from four elite wheat cultivars. The large size of this MAGIC population (1579 progeny), its diverse genetic composition and high levels of recombination all contribute to its value as a genetic resource. Applications of this resource include interrogation of the wheat genome and the analysis of gene-trait association in agronomically important wheat phenotypes. Here, we report the utilization of a MAGIC population for the first time for linkage map construction. We have constructed a linkage map with 1162 DArT, single nucleotide polymorphism and simple sequence repeat markers distributed across all 21 chromosomes. We benchmark this map against a high-density DArT consensus map created by integrating more than 100 biparental populations. The linkage map forms the basis for further exploration of the genetic architecture within the population, including characterization of linkage disequilibrium, founder contribution and inclusion of an alien introgression into the genetic map. Finally, we demonstrate the application of the resource for quantitative trait loci mapping using the complex traits plant height and hectolitre weight as a proof of principle.     

Estimating sample sizes for DNA barcoding

Sample size has long been one of the basic issues since the start of the DNA barcoding initiative and the global biodiversity investigation. As a contribution to resolving this problem, we propose a simple resampling approach to estimate several key sampling sizes for a DNA barcoding project. We illustrate our approach using both structured populations simulated under coalescent and real species of skipper butterflies. We found that sample sizes widely used in DNA barcoding are insufficient to assess the genetic diversity of a species, population structure impacts the estimation of the sample sizes, and hence will bias the species identification potentially.     

Minimum sample sizes for invasion genomics: Empirical investigation in an invasive whitefly

Analysis of population genetics provides insights into the evolutionary processes, among which the sample size choice is per se a crucial issue in the analysis. Genome‐wide high‐throughput techniques based on RADseq have been increasingly used in studies on the population genomics of invasive species. However, there is little information available regarding optimal sample sizes for analyzing population genomics of invasive species. In this study, we first use type IIB endonucleases restriction site‐associated DNA (2b‐RAD) to mine thousands of single nucleotide polymorphisms (SNPs) for native and introduced populations in Q1 clade (SPB and 17JN) and Q2 clade (ISQ and UAS0601) of the whitefly, Bemisia tabaci (Gennadius) MED (also known as B. tabaci biotype Q). Then, we used resampling techniques to create simulated populations with a random subset of individuals and 3,000 SNPs to determine how many individuals should be sampled for accurate estimates of intra‐ and interpopulation genetic diversity. We calculated the intrapopulation genetic diversity parameters (unbiased expected heterozygosity, observed heterozygosity, and the number of effect alleles) and pairwise genetic differentiation F_ST; finally, an ad hoc statistic, ΔK, was used to determine the optimal value. Our results showed that a sample size greater than four individuals (n ≥ 4) has little impact on estimates of genetic diversity within whitefly populations; moreover, precise estimate of F_ST can be easily achieved at a very small simple size (n = 3 or 4). Our results will provide in‐depth understanding of the optimization of sampling schemes in population genomics of invasive species.     

Efficiency of selective genotyping for genetic analysis of complex traits and potential applications in crop improvement

Selective genotyping of individuals from the two tails of the phenotypic distribution of a population provides a cost efficient alternative to analysis of the entire population for genetic mapping. Past applications of this approach have been confounded by the small size of entire and tail populations, and insufficient marker density, which result in a high probability of false positives in the detection of quantitative trait loci (QTL). We studied the effect of these factors on the power of QTL detection by simulation of mapping experiments using population sizes of up to 3,000 individuals and tail population sizes of various proportions, and marker densities up to one marker per centiMorgan using complex genetic models including QTL linkage and epistasis. The results indicate that QTL mapping based on selective genotyping is more powerful than simple interval mapping but less powerful than inclusive composite interval mapping. Selective genotyping can be used, along with pooled DNA analysis, to replace genotyping the entire population, for mapping QTL with relatively small effects, as well as linked and interacting QTL. Using diverse germplasm including all available genetics and breeding materials, it is theoretically possible to develop an “all-in-one plate” approach where one 384-well plate could be designed to map almost all agronomic traits of importance in a crop species. Selective genotyping can also be used for genomewide association mapping where it can be integrated with selective phenotyping approaches. We also propose a breeding-to-genetics approach, which starts with identification of extreme phenotypes from segregating populations generated from multiple parental lines and is followed by rapid discovery of individual genes and combinations of gene effects together with simultaneous manipulation in breeding programs.     

The use of MapPop1.0 for choosing a QTL mapping sample from an advanced backcross population

QTL detection is a good way to assess the genetic basis of quantitative traits such as the plant response to its environment, but requires large mapping populations. Experimental constraints, however, may require a restriction of the population size, risking a decrease in the quality level of QTL mapping. The purpose of this paper was to test if an advanced backcross population sample chosen by MapPop 1.0 could limit the effect of size restriction and improve the QTL detection when compared to random samples. We used the genotypic and phenotypic data obtained for 280 genotypes, considered as the reference population. The “MapPop sample” of 100 genotypes was first compared to the reference population, and genetic maps, genotypic and phenotypic data and QTL results were analysed. Despite the increase in donor allele frequency in the MapPop sample, this did not lead to an increase of the genetic map length or a biased phenotypic distribution. Three QTL among the 10 QTL found in the reference population were also detected in the MapPop sample. Next, the MapPop sample results were compared to those from 500 random samples of the same size. The main conclusion was that the MapPop software avoided the selection of biased samples and the detection of false QTL and appears particularly interesting to select a sample from an unbalanced population.     

Simultaneous improvement and genetic dissection of grain yield and its related traits in a backbone parent of hybrid rice (Oryza sativa L.) using selective introgression

Three populations with a total of 125 BC2F3:4 introgression lines (ILs) selected for high yields from three BC2F2 populations were used for genetic dissection of rice yield and its related traits. The progeny testing in replicated phenotyping across two environments and genotyping with 140 polymorphic simple sequence repeat markers allowed the identification of 21 promising ILs that had significantly higher yields than the recurrent parent Shuhui527 (SH527). A total of 94 quantitative trait loci (QTL) were identified using the selective introgression method based on Chi-squared (χ ²) and multi-locus probability tests and the RSTEP-LRT method based on stepwise regression. These QTL were mostly mapped to 12 clusters on seven rice chromosomes. Several important properties of the QTL affecting grain yield (GY) and its related traits were revealed. The first one was the presence of strong and frequent non-random associations between or among QTL that affect low-heritability traits (GY and spikelet number per panicle, SN) in the ILs with high trait values. Second, beneficial alleles at 88.9 % GY and 75 % SN QTL for increased productivity were from the donors, suggesting that direct phenotypic selection for high yield in our introgression breeding program was a powerful way to transfer beneficial alleles at many loci from the donors into SH527. Third, most QTL were in clusters with large effects on multiple traits, which should be the focal points in further investigations and marker-assisted selection in rice. The majority of the QTL identified were expressed only in one of the environments, suggesting that differential expression of QTL in different environments is the primary genetic basis of genotype × environment interaction. Finally, a large variation in both the direction and magnitude of QTL effects was detected for different donor alleles at seven QTL in the same genetic background and environments. This finding suggests the possible presence of functional diversity among the donor alleles at these loci. The promising ILs and QTL identified provide valuable materials and genetic information for further improving the yield potential of SH527, which is a backbone restorer of hybrid rice in China.     

Using a reference population yardstick to calibrate and compare genetic diversity reported in different studies: an example from the brown bear

In species with large geographic ranges, genetic diversity of different populations may be well studied, but differences in loci and sample sizes can make the results of different studies difficult to compare. Yet, such comparisons are important for assessing the status of populations of conservation concern. We propose a simple approach of using a single well-studied reference population as a ‘yardstick'' to calibrate results of different studies to the same scale, enabling comparisons. We use a well-studied large carnivore, the brown bear (Ursus arctos), as a case study to demonstrate the approach. As a reference population, we genotyped 513 brown bears from Slovenia using 20 polymorphic microsatellite loci. We used this data set to calibrate and compare heterozygosity and allelic richness for 30 brown bear populations from 10 different studies across the global distribution of the species. The simplicity of the reference population approach makes it useful for other species, enabling comparisons of genetic diversity estimates between previously incompatible studies and improving our understanding of how genetic diversity is distributed throughout a species range.