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.     

Effectiveness of selective genotyping for detection of quantitative trait loci: an analysis of grain and malt quality traits in three barley populations.

Marker genotype data and grain and malt quality phenotype data from three barley (Hordeum vulgare L.) mapping populations were used to investigate the feasibility of selective genotyping for detection of quantitative trait loci (QTLs). With selective genotyping, only individuals with high and low phenotypic values for the trait of interest are genotyped. Here, genotyping of 10 to 70% of each population (i.e., 5 to 35% in each tail of the phenotypic distribution) was considered. Genomic positions detected by selective genotyping were compared to QTL position estimates from interval mapping analysis using marker genotype data from the entire population. Selective genotyping reliably detected almost all of the mapped QTLs, often with only 10% of the population genotyped. Selective genotyping also detected spurious QTLs in regions of the genome where no significant QTL had been mapped. Even with additional genotyping to verify putative QTLs, the total genotyping effort for detection of QTLs for a single trait by selective genotyping was usually less than 30% of that required for conventional interval mapping. Simultaneous investigation of two or more traits by selective genotyping would require additional genotyping effort, but could still be worthwhile.     

Association mapping in forest trees and fruit crops

Association mapping (AM), also known as linkage disequilibrium (LD) mapping, is a viable approach to overcome limitations of pedigree-based quantitative trait loci (QTL) mapping. In AM, genotypic and phenotypic correlations are investigated in unrelated individuals. Unlike QTL mapping, AM takes advantage of both LD and historical recombination present within the gene pool of an organism, thus utilizing a broader reference population. In plants, AM has been used in model species with available genomic resources. Pursuing AM in tree species requires both genotyping and phenotyping of large populations with unique architectures. Recently, genome sequences and genomic resources for forest and fruit crops have become available. Due to abundance of single nucleotide polymorphisms (SNPs) within a genome, along with availability of high-throughput resequencing methods, SNPs can be effectively used for genotyping trees. In addition to DNA polymorphisms, copy number variations (CNVs) in the form of deletions, duplications, and insertions also play major roles in control of expression of phenotypic traits. Thus, CNVs could provide yet another valuable resource, beyond those of microsatellite and SNP variations, for pursuing genomic studies. As genome-wide SNP data are generated from high-throughput sequencing efforts, these could be readily reanalysed to identify CNVs, and subsequently used for AM studies. However, forest and fruit crops possess unique architectural and biological features that ought to be taken into consideration when collecting genotyping and phenotyping data, as these will also dictate which AM strategies should be pursued. These unique features as well as their impact on undertaking AM studies are outlined and discussed.     

A New Simple Method for Improving QTL Mapping Under Selective Genotyping

The selective genotyping approach, where only individuals from the high and low extremes of the trait distribution are selected for genotyping and the remaining individuals are not genotyped, has been known as a cost-saving strategy to reduce genotyping work and can still maintain nearly equivalent efficiency to complete genotyping in QTL mapping. We propose a novel and simple statistical method based on the normal mixture model for selective genotyping when both genotyped and ungenotyped individuals are fitted in the model for QTL analysis. Compared to the existing methods, the main feature of our model is that we first provide a simple way for obtaining the distribution of QTL genotypes for the ungenotyped individuals and then use it, rather than the population distribution of QTL genotypes as in the existing methods, to fit the ungenotyped individuals in model construction. Another feature is that the proposed method is developed on the basis of a multiple-QTL model and has a simple estimation procedure similar to that for complete genotyping. As a result, the proposed method has the ability to provide better QTL resolution, analyze QTL epistasis, and tackle multiple QTL problem under selective genotyping. In addition, a truncated normal mixture model based on a multiple-QTL model is developed when only the genotyped individuals are considered in the analysis, so that the two different types of models can be compared and investigated in selective genotyping. The issue in determining threshold values for selective genotyping in QTL mapping is also discussed. Simulation studies are performed to evaluate the proposed methods, compare the different models, and study the QTL mapping properties in selective genotyping. The results show that the proposed method can provide greater QTL detection power and facilitate QTL mapping for selective genotyping. Also, selective genotyping using larger genotyping proportions may provide roughly equivalent power to complete genotyping and that using smaller genotyping proportions has difficulties doing so. The R code of our proposed method is available on http://www.stat.sinica.edu.tw/chkao/.     

Power of different sampling strategies to detect quantitative trait loci variance effects

Summary Many studies have shown that segregating quantitative trait loci (QTL) can be detected via linkage to genetic markers. Power to detect a QTL effect on the trait mean as a function of the number of individuals genotyped for the marker is increased by selectively genotyping individuals with extreme values for the quantitative trait. Computer simulations were employed to study the effect of various sampling strategies on the statistical power to detect QTL variance effects. If only individuals with extreme phenotypes for the quantitative trait are selected for genotyping, then power to detect a variance effect is less than by random sampling. If 0.2 of the total number of individuals genotyped are selected from the center of the distribution, then power to detect a variance effect is equal to that obtained with random selection. Power to detect a variance effect was maximum when 0.2 to 0.5 of the individuals selected for genotyping were selected from the tails of the distribution and the remainder from the center.     

Empirical bayes method for incorporating data from multiple genome scans

Individual genome scans tend to have low power and can produce markedly biased estimates of QTL effects. Further, the confidence interval for their location is often prohibitively large for subsequent fine mapping and positional cloning. Given that a large number of genome scans have been conducted, not to mention the large number of variables and subsets tested, it is difficult to confidently rule out type 1 error as an explanation for significant effects even when there is apparent replication in a separate data set. We adapted Empirical Bayes (EB) methods [1] to analyze data from multiple genome scans simultaneously and alleviate each of these problems while still allowing for different QTL population effects across studies. We investigated the effects of using the EB method to include data from background studies to update the results of a single study of interest via simulation and demonstrated that it has a stable confidence level over a wide range of parameters defining the background studies and increased the power to detect linkage, even when some of the background studies were null or had QTL effect at other markers. This EB method for incorporating data from multiple studies into genome scan analyses seems promising.     

Power of QTL mapping experiments in commercial Atlantic salmon populations, exploiting linkage and linkage disequilibrium and effect of limited recombination in males

Whereas detection and positioning of genes that affect quantitative traits (quantitative trait loci (QTL)) using linkage mapping uses only information from recombinants in the genotyped generations, linkage disequilibrium (LD) mapping uses historical recombinants. Thus, whereas linkage mapping requires large family sizes to detect and accurately position QTL, LD mapping is more dependent on the number of families sampled from the population. In commercial Atlantic salmon breeding programmes, only a small number of individuals per family are routinely phenotyped for traits such as disease resistance and meat colour. In this paper, we assess the power and accuracy of combined linkage disequilibrium linkage analysis (LDLA) to detect QTL in the commercial population using simulation. When 15 half-sib sire families (each sire mated to 30 dams, each dam with 10 progeny) were sampled from the population for genotyping, we were able to detect a QTL explaining 10% of the phenotypic variance in 85% of replicates and position this QTL within 3 cM of the true position in 70% of replicates. When recombination was absent in males, a feature of the salmon genome, power to detect QTL increased; however, the accuracy of positioning the QTL was decreased. By increasing the number of sire families sampled from the population to be genotyped to 30, we were able to increase both the proportion of QTL detected and correctly positioned (even with no recombination in males). QTL with much smaller effect could also be detected. The results suggest that even with the existing recording structure in commercial salmon breeding programmes, there is considerable power to detect and accurately position QTL using LDLA.     

Comparison of methods for analysis of selective genotyping survival data

Survival traits and selective genotyping datasets are typically not normally distributed, thus common models used to identify QTL may not be statistically appropriate for their analysis. The objective of the present study was to compare models for identification of QTL associated with survival traits, in particular when combined with selective genotyping. Data were simulated to model the survival distribution of a population of chickens challenged with Marek disease virus. Cox proportional hazards (CPH), linear regression (LR), and Weibull models were compared for their appropriateness to analyze the data, ability to identify associations of marker alleles with survival, and estimation of effects when all individuals were genotyped (full genotyping) and when selective genotyping was used. Little difference in power was found between the CPH and the LR model for low censoring cases for both full and selective genotyping. The simulated data were not transformed to follow a Weibull distribution and, as a result, the Weibull model generally resulted in less power than the other two models and overestimated effects. Effect estimates from LR and CPH were unbiased when all individuals were genotyped, but overestimated when selective genotyping was used. Thus, LR is preferred for analyzing survival data when the amount of censoring is low because of ease of implementation and interpretation. Including phenotypic data of non-genotyped individuals in selective genotyping analysis increased power, but resulted in LR having an inflated false positive rate, and therefore the CPH model is preferred for this scenario, although transformation of the data may also make the Weibull model appropriate for this case. The results from the research presented herein are directly applicable to interval mapping analyses.     

The effects of selective genotyping on estimates of proportion of recombination between linked quantitative trait loci

Selective genotyping is the marker assay of only the more extreme phenotypes for a quantitative trait and is intended to increase the efficiency of quantitative trait loci (QTL) mapping. We show that selective genotyping can bias estimates of the recombination frequency between linked QTLs — upwardly when QTLs are in repulsion phase, and downwardly when QTLs are in coupling phase. We examined these biases under simple models involving two QTLs segregating in a backcross or F₂ population, using both analytical models and computer simulations. We found that bias is a function of the proportion selected, the magnitude of QTL effects, distance between QTLs and the dominance of QTLs. Selective genotyping thus may decrease the power of mapping multiple linked QTLs and bias the construction of a marker map. We suggest a large proportion than previously suggested (50%) or the entire population be genotyped if linked QTLs of large effects (explain > 10% phenotypic variance) are evident. New models need to be developed to explicitly incorporate selection into QTL map construction.     

Novel linkage mapping approach using DNA pooling in human and animal genetics. II. Detection of quantitative traits loci in dairy cattle

Selective DNA pooling is an advanced methodology for linkage mapping of quantitative trait loci (QTL) in farm animals. The principle is based on densitometric estimates of marker allele frequency in pooled DNA samples of phenotypically extreme individuals from half-sib, backcross and F(2) experimental designs in farm animals. This methodology provides a rapid and efficient analysis of a large number of individuals with short tandem repeat markers that are essential to detect QTL through the genome - wide searching approach. Several strategies involving whole genome scanning with a high statistical power have been developed for systematic search to detect the quantitative traits loci and linked loci of complex traits. In recent studies, greater success has been achieved in mapping several QTLs in Israel-Holstein cattle using selective DNA pooling. This paper outlines the currently emerged novel strategies of linkage mapping to identify QTL based on selective DNA pooling with more emphasis on its theoretical pre-requisite to detect linked QTLs, applications, a general theory for experimental half-sib designs, the power of statistics and its feasibility to identify genetic markers linked QTL in dairy cattle. The study reveals that the application of selective DNA pooling in dairy cattle can be best exploited in the genome-wide detection of linked loci with small and large QTL effects and applied to a moderately sized half-sib family of about 500 animals.     

High-density genotyping: an overkill for QTL mapping? Lessons learned from a case study in maize and simulations

High-density genotyping is extensively exploited in genome-wide association mapping studies and genomic selection in maize. By contrast, linkage mapping studies were until now mostly based on low-density genetic maps and theoretical results suggested this to be sufficient. This raises the question, if an increase in marker density would be an overkill for linkage mapping in biparental populations, or if important QTL mapping parameters would benefit from it. In this study, we addressed this question using experimental data and a simulation based on linkage maps with marker densities of 1, 2, and 5 cM. QTL mapping was performed for six diverse traits in a biparental population with 204 doubled haploid maize lines and in a simulation study with varying QTL effects and closely linked QTL for different population sizes. Our results showed that high-density maps neither improved the QTL detection power nor the predictive power for the proportion of explained genotypic variance. By contrast, the precision of QTL localization, the precision of effect estimates of detected QTL, especially for small and medium sized QTL, as well as the power to resolve closely linked QTL profited from an increase in marker density from 5 to 1 cM. In conclusion, the higher costs for high-density genotyping are compensated for by more precise estimates of parameters relevant for knowledge-based breeding, thus making an increase in marker density for linkage mapping attractive.     

Optimal designs for linkage disequilibrium mapping and candidate gene association tests in livestock populations

Simulation was used to evaluate the performance of different selective genotyping strategies when using linkage disequilibrium across large half-sib families to position a QTL within a previously defined genomic region. Strategies examined included standard selective genotyping and different approaches of discordant and concordant sib selection applied to arbitrary or selected families. Strategies were compared as a function of effect and frequency of QTL alleles, heritability, and phenotypic expression of the trait. Large half-sib families were simulated for 100 generations and 2% of the population was genotyped in the final generation. Simple ANOVA was applied and the marker with the greatest F-value was considered the most likely QTL position. For traits with continuous phenotypes, genotyping the most divergent pairs of half-sibs from all families was the best strategy in general, but standard selective genotyping was somewhat more precise when heritability was low. When the phenotype was distributed in ordered categories, discordant sib selection was the optimal approach for positioning QTL for traits with high heritability and concordant sib selection was the best approach when genetic effects were small. Genotyping of a few selected sibs from many families was generally more efficient than genotyping many individuals from a few highly selected sires.     

Genetic mapping of quantitative phenotypic traits in Saccharomyces cerevisiae

Saccharomyces cerevisiae has become a favorite production organism in industrial biotechnology presenting new challenges to yeast engineers in terms of introducing advantageous traits such as stress tolerances. Exploring subspecies diversity of S. cerevisiae has identified strains that bear industrially relevant phenotypic traits. Provided that the genetic basis of such phenotypic traits can be identified inverse engineering allows the targeted modification of production strains. Most phenotypic traits of interest in S. cerevisiae strains are quantitative, meaning that they are controlled by multiple genetic loci referred to as quantitative trait loci (QTL). A straightforward approach to identify the genetic basis of quantitative traits is QTL mapping which aims at the allocation of the genetic determinants to regions in the genome. The application of high-density oligonucleotide arrays and whole-genome re-sequencing to detect genetic variations between strains has facilitated the detection of large numbers of molecular markers thus allowing high-resolution QTL mapping over the entire genome. This review focuses on the basic principle and state of the art of QTL mapping in S. cerevisiae. Furthermore we discuss several approaches developed during the last decade that allow down-scaling of the regions identified by QTL mapping to the gene level. We also emphasize the particular challenges of QTL mapping in nonlaboratory strains of S. cerevisiae.     

Quantitative trait loci mapping of phenotypic plasticity and genotype-environment interactions in plant and insect performance

Community genetic studies generally ignore the plasticity of the functional traits through which the effect is passed from individuals to the associated community. However, the ability of organisms to be phenotypically plastic allows them to rapidly adapt to changing environments and plasticity is commonly observed across all taxa. Owing to the fitness benefits of phenotypic plasticity, evolutionary biologists are interested in its genetic basis, which could explain how phenotypic plasticity is involved in the evolution of species interactions. Two current ideas exist: (i) phenotypic plasticity is caused by environmentally sensitive loci associated with a phenotype; (ii) phenotypic plasticity is caused by regulatory genes that simply influence the plasticity of a phenotype. Here, we designed a quantitative trait loci (QTL) mapping experiment to locate QTL on the barley genome associated with barley performance when the environment varies in the presence of aphids, and the composition of the rhizosphere. We simultaneously mapped aphid performance across variable rhizosphere environments. We mapped main effects, QTL × environment interaction (QTL×E), and phenotypic plasticity (measured as the difference in mean trait values) for barley and aphid performance onto the barley genome using an interval mapping procedure. We found that QTL associated with phenotypic plasticity were co-located with main effect QTL and QTL×E. We also located phenotypic plasticity QTL that were located separately from main effect QTL. These results support both of the current ideas of how phenotypic plasticity is genetically based and provide an initial insight into the functional genetic basis of how phenotypically plastic traits may still be important sources of community genetic effects.     

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.     

Extent and consistency across generations of linkage disequilibrium in commercial layer chicken breeding populations

Recent studies report a surprisingly high degree of marker-to-marker linkage disequilibrium (LD) in ruminant livestock populations. This has important implications for QTL mapping and marker-assisted selection. This study evaluated LD between microsatellite markers in a number of breeding populations of layer chickens using the standardized chi-square (chi(2')) measure. The results show appreciable LD among markers separated by up to 5 cM, decreasing rapidly with increased separation between markers. The LD within 5 cM was strongly conserved across generations and differed among chromosomal regions. Using marker-to-marker LD as an indication for marker-QTL LD, a genome scan of markers spaced 2 cM apart at moderate power would have good chances of uncovering most QTL segregating in these populations. However, of markers showing significant trait associations, only 57% are expected to be within 5 cM of the responsible QTL, and the remainder will be up to 20 cM away. Thus, high-resolution LD mapping of QTL will require dense marker genotyping across the region of interest to allow for interval mapping of the QTL.     

Selective genotyping for determination of linkage between a marker locus and a quantitative trait locus

Summary Selective genotyping is the term used when the determination of linkage between marker loci and quantitative trait loci (QTL) affecting some particular trait is carried out by genotyping only individuals from the high and low phenotypic tails of the entire sample population. Selective genotyping can markedly decrease the number of individuals genotyped for a given power at the expense of an increase in the number of individuals phenotyped. The optimum proportion of individuals genotyped from the point of view of minimizing costs for a given experimental power depends strongly on the cost of completely genotyping an individual for all of the markers included in the experiment (including the costs of obtaining a DNA sample) relative to the cost of rearing and trait evaluation of an individual. However, in single trait studies, it will almost never be useful to genotype more than the upper and lower 25% of a population. It is shown that the observed difference in quantitative trait values associated with alternative marker genotypes in the selected population can be much greater than the actual gene effect at the quantitative trait locus when the entire population is considered. An expression and a figure is provided for converting observed differences under selective genotyping to actual gene effects.     

On the generalized poisson regression mixture model for mapping quantitative trait loci with count data

Statistical methods for mapping quantitative trait loci (QTL) have been extensively studied. While most existing methods assume normal distribution of the phenotype, the normality assumption could be easily violated when phenotypes are measured in counts. One natural choice to deal with count traits is to apply the classical Poisson regression model. However, conditional on covariates, the Poisson assumption of mean-variance equality may not be valid when data are potentially under- or overdispersed. In this article, we propose an interval-mapping approach for phenotypes measured in counts. We model the effects of QTL through a generalized Poisson regression model and develop efficient likelihood-based inference procedures. This approach, implemented with the EM algorithm, allows for a genomewide scan for the existence of QTL throughout the entire genome. The performance of the proposed method is evaluated through extensive simulation studies along with comparisons with existing approaches such as the Poisson regression and the generalized estimating equation approach. An application to a rice tiller number data set is given. Our approach provides a standard procedure for mapping QTL involved in the genetic control of complex traits measured in counts.     

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.     

Quantitative Trait Loci Affecting Response to Crowding Stress in an F2 Generation of Rainbow Trout Produced Through Phenotypic Selection

Selective breeding programs for salmonids typically aim to improve traits associated with growth and disease resistance. It has been established that stressors common to production environments can adversely affect these and other traits which are important to producers and consumers. Previously, we employed phenotypic selection to create families that exhibit high or low plasma cortisol concentrations in response to crowding stress. Subsequent crosses of high × low phenotypes founded a multigenerational breeding scheme with the aim of dissecting the genetic basis for variation underlying stress response through the identification of quantitative trait loci (QTL). Multiple methods of QTL analyses differing in their assumptions of homozygosity of the causal alleles in the grandparental generation yielded similar results in the F₁ generation, and the analysis of two stress response phenotype measurement indexes were highly correlated. In the current study, we conducted a genome scan with microsatellites to detect QTL in the F₂ generation of two families created through phenotypic selection and having larger numbers of offspring than families screened in the previous generation. Seven suggestive and three significant QTL were detected, seven of which were not previously detected in the National Center for Cool and Cold Water Aquaculture germplasm, bringing the total number of chromosomes containing significant and suggestive stress response QTL to 4 and 15, respectively. One significant QTL which peaks at 7 cM on chromosome Omy12 spans 12 cM and explains 25 % of the phenotypic variance in family 2008052 particularly warrants further investigation. Five QTL with significant parent-of-origin effects were detected in family 2008052, including two QTL on Omy12. The 95 % confidence intervals for the remaining QTL we detected were broad, requiring validation and fine mapping with other genotyping approaches and mapping strategies. These results will facilitate identification of potential casual alleles that can be employed in strategies aimed at better understanding the genetic and physiological basis of stress responses to crowding in rainbow trout aquaculture production.