Selection and subsequent analysis of sib pair data for QTL detection.

Haseman and Elston (1972) developed a robust regression method for the detection of linkage between a marker and a quantitative trait locus (QTL) using sib pair data. The principle underlying this method is that the difference in phenotypes between pairs of sibs becomes larger as they share a decreasing number of alleles at a particular QTL identical by descent (IBD) from their parents. In this case, phenotypically very different sibs will also on average share a proportion of alleles IBD at any marker linked to the QTL that is lower than the expected value of 0.5. Thus, the deviation of the proportion of marker alleles IBD from the expected value in pairs of sibs selected to be phenotypically different (i.e. discordant) can provide a test for the presence of a QTL. A simple regression method for QTL detection in sib pairs selected for high phenotypic differences is presented here. The power of the analytical method was found to be greater than the power obtained using the standard analysis when samples of sib pairs with high phenotypic differences were used. However, the use of discordant sib pairs was found to be less powerful for QTL detection than alternative selective genotyping schemes based on the phenotypic values of the sibs except with intense selection, when its advantage was only marginal. The most effective selection scheme overall was the use of sib pairs from entire families selected on the basis of high within-family variance for the trait in question. There is little effect of selection on QTL position estimates, which are in good agreement with the simulated values. However, QTL variance estimates are biased to a greater or lesser degree, depending on the selection method.     

Selective genotyping for QTL detection using sib pair analysis in outbred populations with hierarchical structures

A simulation study illustrates the effects of the inclusion of half-sib pairs as well as the effects of selective genotyping on the power of detection and the parameter estimates in a sib pair analysis of data from an outbred population. The power of QTL detection obtained from samples of sib pairs selected according to their within family variance or according to the mean within family variance within half sib family was compared and contrasted with the power obtained when only full sib pair analysis was used. There was an increase in power (4–16%) and decrease in the bias of parameter estimates with the use of half-sib information. These improvements in power and parameter estimates depended on the number of the half sib pairs (half sib family size). Almost the same power as that obtained using all the available sib pairs could be achieved by selecting only 50–60% the animals. The most effective method was to select both full and half sib pairs on the basis of high within full sib family variance for the trait in question. The QTL position estimates were in general slightly biased towards the center of the chromosome and the QTL variance estimates were biased upwards, there being quite large differences in bias depending on the selection method.     

The power to detect linkage disequilibrium with quantitative traits in selected samples

Results from power studies for linkage detection have led to many ongoing and planned collections of phenotypically extreme nuclear families. Given the great expense of collecting these families and the imminent availability of a dense diallelic marker map, the families are likely to be used in allelic-association as well as linkage studies. However, optimal selection strategies for linkage may not be equally powerful for association. We examine the power to detect linkage disequilibrium for quantitative traits after phenotypic selection. The results encompass six selection strategies that are in widespread use, including single selection (two designs), affected sib pairs, concordant and discordant pairs, and the extreme-concordant and -discordant design. Selection of sibships on the basis of one extreme proband with high or low trait scores provides as much power as discordant sib pairs but requires the screening and phenotyping of substantially fewer initial families from which to select. Analysis of the role of allele frequencies within each selection design indicates that common trait alleles generally offer the most power, but similarities between the marker- and trait-allele frequencies are much more important than the trait-locus frequency alone. Some of the most widespread selection designs, such as single selection, yield power gains only when both the marker and quantitative trait loci (QTL) are relatively rare in the population. In contrast, discordant pairs and the extreme-proband design provide power for the broadest range of QTL-marker-allele frequency differences. Overall, proband selection from either tail provides the best balance of power, robustness, and simplicity of ascertainment for family-based association analysis.     

The power of interval mapping of quantitative trait loci, using selected sib pairs.

The interval-mapping procedure of Fulker and Cardon for analysis of a quantitative-trait loci (QTL) is extended for application to selected samples of sib pairs. Phenotypic selection of sib pairs, which is known to yield striking increases in power when a single marker is used, provides further increases in power when the interval-mapping approach is used. The greatest benefits of the combined approach are apparent with coarse maps, where QTLs of relatively modest (15%-20%) heritability can be detected with widely spaced markers (40-60 cM apart) in reasonably sized sibling samples. Useful information concerning QTL location is afforded by interval mapping in both selected and unselected samples.     

Genetic Mapping of Quantitative Trait Loci Controlling Growth and Wood Quality Traits in Eucalyptus Grandis Using a Maternal Half-Sib Family and Rapd Markers

Quantitative trait loci (QTL) mapping of forest productivity traits was performed using an open pollinated half-sib family of Eucalyptus grandis. For volume growth, a sequential QTL mapping approach was applied using bulk segregant analysis (BSA), selective genotyping (SG) and cosegregation analysis (CSA). Despite the low heritability of this trait and the heterogeneous genetic background employed for mapping. BSA detected one putative QTL and SG two out of the three later found by CSA. The three putative QTL for volume growth were found to control 13.7% of the phenotypic variation, corresponding to an estimated 43.7% of the genetic variation. For wood specific gravity five QTL were identified controlling 24.7% of the phenotypic variation corresponding to 49% of the genetic variation. Overlapping QTL for CBH, WSG and percentage dry weight of bark were observed. A significant case of digenic epistasis was found, involving unlinked QTL for volume. Our results demonstrate the applicability of the within half-sib design for QTL mapping in forest trees and indicate the existence of major genes involved in the expression of economically important traits related to forest productivity in Eucalyptus grandis. These findings have important implications for marker-assisted tree breeding.     

Whole-genome scan to detect QTL for milk production, conformation, fertility and functional traits in two US Holstein families

A genome scan was conducted in two US Holstein half-sib families to identify quantitative trait loci (QTL) affecting milk production and conformation traits using the granddaughter design. The sires of the two studied families were related as sire and son and had 96 and 212 sons respectively. A total of 221 microsatellite loci were scored in both families. Statistical analysis was performed using two different analytical methods; half-sib least squares regression and Bayesian Monte Carlo Markov Chain. Traits analysed included five traditional milk production traits, somatic cell count, daughter pregnancy rate, male fertility and 20 conformation traits. A total of 47 tests achieved at least genome-wise significance. However, results from the two methods of analysis were only concordant for QTL location and level of significance in eight instances.     

Power of two and three-generation QTL mapping experiments in an outbred population containing full-sib or half-sib families

QTL mapping experiments involve many animals to be genotyped and performance tested. Consequently, experimental designs need to be optimized to minimize the costs of data collection and genotyping. The present study has analyzed the power and efficiency of experiments with two or three-generation family structures containing full-sib families, half-sib families, or both. The focus was on data from one outbred population because the main interest is to locate genes that can be used for within-line selection. For a two generation experiment more animals had to be typed for marker loci to obtain a certain power than for a three generation experiment. Fewer trait values, however, had to be obtained for a two-generation experiment than for a three-generation experiment. A two or three-generation family structure with full-sib offspring was more efficient than a two or three-generation family structure with half-sib offspring. A family structure with full-sib grand-offspring, however, was less efficient than a family structure with half-sib grand-offspring. For the most efficient family structure each pair of parents had full-sib offspring that were genotyped for the marker. For the most-efficient family structure each full-sib offspring had half-sib grand-offspring for which trait values were obtained. For equal power with a heritability of 0.1 and 100 grand-offspring per full-sib offspring, 30-times less marker typings were required for this most efficient family structure than for a two-generation half-sib structure in which marker genotypes and trait values were obtained for half-sib offspring. The effect of heritability and the type of analysis (single marker or interval analysis) on the efficiency of a family structure is described. The results of this study should help to design QTL mapping experiments in an outbred population.     

性状遗传力与QTL方差对标记辅助选择效果的影响

在采用动物模型标记辅助最佳线性无偏预测方法对个体育种值进行估计的基础上,模拟了在一个闭锁群体内连续对单个性状选择10个世代的情形,并系统地比较了性状遗传力和QTL方差对标记辅助选择所获得的遗传进展、QTL增效基因频率和群体近交系数变化的影响。结果表明：在对高遗传力和QTL方差较小的性状实施标记辅助选择时,可望获得更大的遗传进展;遗传力越高,QTL方差越大,则QTL增效基因频率的上升速度越快;遗传力较高时,群体近交系数上升的速度较为缓慢,而QTL方差对群体近交系数上升速度的影响则不甚明显。结合前人关于标记辅助选择相对效率的研究结果,可以认为：当选择性状的遗传力和QTL方差为中等水平时,标记辅助选择可望获得理想的效果。     

A unified sampling approach for multipoint analysis of qualitative and quantitative traits in sib pairs

Recent advances in molecular biology have enhanced the opportunity to conduct multipoint mapping for complex diseases. Concurrently, one sees a growing interest in the use of quantitative traits in linkage studies. Here, we present a multipoint sib-pair approach to locate the map position (tau) of a trait locus that controls the observed phenotype (qualitative or quantitative), along with a measure of statistical uncertainty. This method builds on a parametric representation for the expected identical-by-descent statistic at an arbitrary locus, conditional on an event reflecting the sampling scheme, such as affected sib pairs, for qualitative traits, or extreme discordant (ED) sib pairs, for quantitative traits. Our results suggest that the variance about tau&d4;, the estimator of tau, can be reduced by as much as 60%-70% by reducing the length of intervals between markers by one half. For quantitative traits, we examine the precision gain (measured by the variance reduction in tau&d4;) by genotyping extremely concordant (EC) sib pairs and including them along with ED sib pairs in the statistical analysis. The precision gain depends heavily on the residual correlation of the quantitative trait for sib pairs but considerably less on the allele frequency and exact genetic mechanism. Since complex traits involve multiple loci and, hence, the residual correlation cannot be ignored, our finding strongly suggests that one should incorporate EC sib pairs along with ED sib pairs, in both design and analysis. Finally, we empirically establish a simple linear relationship between the magnitude of precision gain and the ratio of the number of ED pairs to the number of EC pairs. This relationship allows investigators to address issues of cost effectiveness that are due to the need for phenotyping and genotyping subjects.     

Selection bias in quantitative trait loci mapping

A simulation study was performed to see whether selection affected quantitative trait loci (QTL) mapping. Populations under random selection, under selection among full-sib families, and under selection within a full-sib family were simulated each with heritability of 0.3, 0.5, and 0.7. They were analyzed with the marker spacing of 10 cM and 20 cM. The accuracy for QTL detection decreased for the populations under selection within full-sib family. Estimates of QTL effects and positions differed (P < .05) from their input values. The problems could be ignored when mapping a QTL for the populations under selection among full-sib families. A large heritability helped reduction of such problems. When the animals were selected within a full-sib family, the QTL was detected for the populations with heritability of 0.5 or larger using the marker spacing of 10 cM, and with heritability of 0.7 using the marker spacing of 20 cM. This study implied that when selection was introduced, the accuracy for QTL detection decreased and the estimates of QTL effects were biased. A caution was warranted on the decision of data (including selected animals to be genotyped) for QTL mapping.     

Genomic selection in a pig population including information from slaughtered full sibs of boars within a sib-testing program

Genomic selection is becoming a common practise in dairy cattle, but only few works have studied its introduction in pig selection programs. Results described for this species are highly dependent on the considered traits and the specific population structure. This paper aims to simulate the impact of genomic selection in a pig population with a training cohort of performance-tested and slaughtered full sibs. This population is selected for performance, carcass and meat quality traits by full-sib testing of boars. Data were simulated using a forward-in-time simulation process that modeled around 60K single nucleotide polymorphisms and several quantitative trait loci distributed across the 18 porcine autosomes. Data were edited to obtain, for each cycle, 200 sires mated with 800 dams to produce 800 litters of 4 piglets each, two males and two females (needed for the sib test), for a total of 3200 newborns. At each cycle, a subset of 200 litters were sib tested, and 60 boars and 160 sows were selected to replace the same number of culled male and female parents. Simulated selection of boars based on performance test data of their full sibs (one castrated brother and two sisters per boar in 200 litters) lasted for 15 cycles. Genotyping and phenotyping of the three tested sibs (training population) and genotyping of the candidate boars (prediction population) were assumed. Breeding values were calculated for traits with two heritability levels (h²=0.40, carcass traits, and h²=0.10, meat quality parameters) on simulated pedigrees, phenotypes and genotypes. Genomic breeding values, estimated by various models (GBLUP from raw phenotype or using breeding values and single-step models), were compared with the classical BLUP Animal Model predictions in terms of predictive ability. Results obtained for traits with moderate heritability (h²=0.40), similar to the heritability of traits commonly measured within a sib-testing program, did not show any benefit from the introduction of genomic selection. None of the considered genomic models provided improvements in prediction ability of pigs with no recorded phenotype. However, a few advantages were found for traits with low heritability (h²=0.10). These heritability levels are characteristic for meat quality traits recorded after slaughtering or for reproduction or health traits, typically recorded on field and not in performance stations. Other scenarios of data recording and genotyping should be evaluated before considering the implementation of genomic selection in a pig-selection scheme based on sib testing of boars.     

QTL detection with bidirectional and unidirectional selective genotyping: marker-based and trait-based analyses

Selective genotyping of one or both phenotypic extremes of a population can be used to detect linkage between markers and quantitative trait loci (QTL) in situations in which full-population genotyping is too costly or not feasible, or where the objective is to rapidly screen large numbers of potential donors for useful alleles with large effects. Data may be subjected to 'trait-based' analysis, in which marker allele frequencies are compared between classes of progeny defined based on trait values, or to 'marker-based' analysis, in which trait means are compared between progeny classes defined based on marker genotypes. Here, bidirectional and unidirectional selective genotyping were simulated, using population sizes and selection intensities relevant to cereal breeding. Control of Type I error was usually adequate with marker-based analysis of variance or trait-based testing using the normal approximation of the binomial distribution. Bidirectional selective genotyping was more powerful than unidirectional. Trait-based analysis and marker-based analysis of variance were about equally powerful. With genotyping of the best 30 out of 500 lines (6%), a QTL explaining 15% of the phenotypic variance could be detected with a power of 0.8 when tests were conducted at a marker 10 cM from the QTL. With bidirectional selective genotyping, QTL with smaller effects and (or) QTL farther from the nearest marker could be detected. Similar QTL detection approaches were applied to data from a population of 436 recombinant inbred rice lines segregating for a large-effect QTL affecting grain yield under drought stress. That QTL was reliably detected by genotyping as few as 20 selected lines (4.5%). In experimental populations, selective genotyping can reduce costs of QTL detection, allowing larger numbers of potential donors to be screened for useful alleles with effects across different backgrounds. In plant breeding programs, selective genotyping can make it possible to detect QTL using even a limited number of progeny that have been retained after selection.     

Identity-by-descent genomic selection using selective and sparse genotyping for binary traits

Background

Genomic selection (GS) allows estimation of the breeding value of individuals, even for non-phenotyped animals. The aim of the study was to examine the potential of identity-by-descent genomic selection (IBD-GS) in genomic selection for a binary, sib-evaluated trait, using different strategies of selective genotyping. This low-cost GS approach is based on linkage analysis of sparse genome-wide marker loci.

Findings

Lowly to highly heritable (h² = 0.15, 0.30 or 0.60) binary traits with varying incidences (10 to 90%) were simulated for an aquaculture-like population. Genotyping was restricted to the 30% best families according to phenotype, using three genotyping strategies for training sibs. IBD-GS increased genetic gain compared to classical pedigree-based selection; the differences were largest at incidences of 10 to 50% of the desired category (i.e. a relative increase in genetic gain greater than 20%). Furthermore, the relative advantage of IBD-GS increased as the heritability of the trait increased. Differences were small between genotyping strategies, and most of the improvement was achieved by restricting genotyping to sibs with the least common binary phenotype. Genetic gains of IBD-GS relative to pedigree-based models were highest at low to moderate (10 to 50%) incidences of the category selected for, but decreased substantially at higher incidences (80 to 90%).

Conclusions

The IBD-GS approach, combined with sparse and selective genotyping, is well suited for genetic evaluation of binary traits. Genetic gain increased considerably compared with classical pedigree-based selection. Most of the improvement was achieved by selective genotyping of the sibs with the least common (minor) binary category phenotype. Furthermore, IBD-GS had greater advantage over classical pedigree-based models at low to moderate incidences of the category selected for.     

Fine mapping and imprinting analysis for fatness trait QTLs in pigs

Quantitative trait loci (QTL) for fatness traits were reported recently in an experimental Meishan × Large White and Landrace F₂ cross. To further investigate the regions on pig Chr 2 (SSC2), SSC4, and SSC7, 25 additional markers from these regions were typed on 800 animals (619 F₂ animals, their F₁ parents, and F₀ grandfathers). Compared with the published maps, a modified order of markers was observed for SSC4 and SSC7. QTL analyses were performed both within the half-sib families as well as across families (line cross). Furthermore, a QTL model accounting for imprinting effects was tested. Information content could be increased considerably on all three chromosomes. Evidence for the backfat thickness QTL on SSC7 was increased, and the location could be reduced to a 33-cM confidence interval. The QTL for intramuscular fat on SSC4 could not be detected in this half-sib analysis, whereas under the line cross model a suggestive QTL on a different position on SSC4 was detected. For SSC2, in the half-sib analysis, a suggestive QTL for backfat thickness was detected with the best position at 26 cM. Imprinting analysis, however, revealed a genome-wise, significant, paternally expressed QTL on SSC2 with the best position at 63 cM. Our results suggest that this QTL is different from the previously reported paternally expressed QTL for muscle mass and fat deposition on the distal tip of SSC2p. Received: 15 October 1999 / Accepted: 21 February 2000     

Natural selection, evolvability and bias due to environmental covariance in the field in an annual plant

Estimates of the form and magnitude of natural selection based on phenotypic relationships between traits and fitness measures can be biased when environmental factors influence both relative fitness and phenotypic trait values. I quantified genetic variances and covariances, and estimated linear and quadratic selection coefficients, for seven traits of an annual plant grown in the field. For replicates of 50 paternal half-sib families, coefficients of selection were calculated both for individual phenotypic values of the traits and for half-sib family mean values. The potential for evolutionary response was supported by significant heritability and phenotypic directional selection for several traits but contradicted by the absence of significant genetic variation for fitness estimates and evidence of bias in phenotypic selection coefficients due to environmental covariance for at least two of the traits analysed. Only studies of a much wider range of organisms and traits will reveal the frequency and extent of such bias.     

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.     

The efficiency of designs for fine-mapping of quantitative trait loci using combined linkage disequilibrium and linkage

In a simulation study, different designs were compared for efficiency of fine-mapping of QTL. The variance component method for fine-mapping of QTL was used to estimate QTL position and variance components. The design of many families with small size gave a higher mapping resolution than a design with few families of large size. However, the difference is small in half sib designs. The proportion of replicates with the QTL positioned within 3 cM of the true position is 0.71 in the best design, and 0.68 in the worst design applied to 128 animals with a phenotypic record and a QTL explaining 25% of the phenotypic variance. The design of two half sib families each of size 64 was further investigated for a hypothetical population with effective size of 1000 simulated for 6000 generations with a marker density of 0.25 cM and with marker mutation rate 4 × 10^-4 per generation. In mapping using bi-allelic markers, 42~55% of replicated simulations could position QTL within 0.75 cM of the true position whereas this was higher for multi allelic markers (48~76%). The accuracy was lowest (48%) when mutation age was 100 generations and increased to 68% and 76% for mutation ages of 200 and 500 generations, respectively, after which it was about 70% for mutation ages of 1000 generations and older. When effective size was linearly decreasing in the last 50 generations, the accuracy was decreased (56 to 70%). We show that half sib designs that have often been used for linkage mapping can have sufficient information for fine-mapping of QTL. It is suggested that the same design with the same animals for linkage mapping should be used for fine-mapping so gene mapping can be cost effective in livestock populations.     

Identity-by-descent genomic selection using selective and sparse genotyping

Background

Genomic selection methods require dense and widespread genotyping data, posing a particular challenge if both sexes are subject to intense selection (e.g., aquaculture species). This study focuses on alternative low-cost genomic selection methods (IBD-GS) that use selective genotyping with sparse marker panels to estimate identity-by-descent relationships through linkage analysis. Our aim was to evaluate the potential of these methods in selection programs for continuous traits measured on sibs of selection candidates in a typical aquaculture breeding population.

Methods

Phenotypic and genomic data were generated by stochastic simulation, assuming low to moderate heritabilities (0.10 to 0.30) for a Gaussian trait measured on sibs of the selection candidates in a typical aquaculture breeding population that consisted of 100 families (100 training animals and 20 selection candidates per family). Low-density marker genotype data (~ 40 markers per Morgan) were used to trace genomic identity-by-descent relationships. Genotyping was restricted to selection candidates from 30 phenotypically top-ranking families and varying fractions of their phenotypically extreme training sibs. All phenotypes were included in the genetic analyses. Classical pedigree-based and IBD-GS models were compared based on realized genetic gain over one generation of selection.

Results

Genetic gain increased substantially (13 to 32%) with IBD-GS compared to classical selection and was greatest with higher heritability. Most of the extra gain from IBD-GS was obtained already by genotyping the 5% phenotypically most extreme sibs within the pre-selected families. Additional genotyping further increased genetic gains, but these were small when going from genotyping 20% of the extremes to all phenotyped sibs. The success of IBD-GS with sparse and selective genotyping can be explained by the fact that within-family haplotype blocks are accurately traced even with low-marker densities and that most of the within-family variance for normally distributed traits is captured by a small proportion of the phenotypically extreme sibs.

Conclusions

IBD-GS was substantially more effective than classical selection, even when based on very few markers and combined with selective genotyping of small fractions of the population. The study shows that low-cost GS programs can be successful by combining sparse and selective genotyping with pedigree and linkage information.     

Asymptotic variances of QTL estimators with selective DNA pooling

Investigation on QTL-marker linkage usually requires a great number of observed recombinations, inferred from combined analysis of phenotypes and genotypes. To avoid costly individual genotyping, inferences on QTL position and effects can instead make use of marker allele frequencies. DNA pooling of selected samples makes allele frequency estimation feasible for studies involving large sample sizes. Linkage studies in outbred populations have traditionally exploited half-sib family designs; within the animal production context, half-sibships provide large families that are highly suitable for DNA pooling. Estimators for QTL position and effect have been proposed that make use of information from flanking markers. We present formulas derived by the delta method for the asymptotic variance of these estimators.     

Association Mapping in Outbred Populations: Power and Efficiency When Genotyping Parents and Phenotyping Progeny

We develop expressions for the power to detect associations between parental genotypes and offspring phenotypes for quantitative traits. Three different “indirect” experimental designs are considered: full-sib, half-sib, and full-sib–half-sib families. We compare the power of these designs to detect genotype–phenotype associations relative to the common, “direct,” approach of genotyping and phenotyping the same individuals. When heritability is low, the indirect designs can outperform the direct method. However, the extra power comes at a cost due to an increased phenotyping effort. By developing expressions for optimal experimental designs given the cost of phenotyping relative to genotyping, we show how the extra costs associated with phenotyping a large number of individuals will influence experimental design decisions. Our results suggest that indirect association studies can be a powerful means of detecting allelic associations in outbred populations of species for which genotyping and phenotyping the same individuals is impractical and for life history and behavioral traits that are heavily influenced by environmental variance and therefore best measured on groups of individuals. Indirect association studies are likely to be favored only on purely economical grounds, however, when phenotyping is substantially less expensive than genotyping. A web-based application implementing our expressions has been developed to aid in the design of indirect association studies.