Comparison of three multitrait methods for QTL detection

A comparison of power and accuracy of estimation of position and QTL effects of three multitrait methods and one single trait method for QTL detection was carried out on simulated data, taking into account the mixture of full and half-sib families. One multitrait method was based on a multivariate function as the penetrance function (MV). The two other multitrait methods were based on univariate analysis of linear combination(s) (LC) of the traits. One was obtained by a principal component analysis (PCA) performed on the phenotypic data. The second was based on a discriminate analysis (DA). It calculates a LC of the traits at each position, maximising the ratio between the genetic and the residual variabilities due to the putative QTL. Due to its number of parameters, MV was less powerful and accurate than the other methods. In general, DA better detected QTL, but it had lower accuracy for the QTL effect estimation when the detection power was low, due to higher bias than the other methods. In this case, PCA was better. Otherwise, PCA was slightly less powerful and accurate than DA. Compared to the single trait method, power can be improved by 30% to 100% with multitrait methods.     

Mapping of quantitative trait loci for flesh colour and growth traits in Atlantic salmon (Salmo salar)

Background

Flesh colour and growth related traits in salmonids are both commercially important and of great interest from a physiological and evolutionary perspective. The aim of this study was to identify quantitative trait loci (QTL) affecting flesh colour and growth related traits in an F2 population derived from an isolated, landlocked wild population in Norway (Byglands Bleke) and a commercial production population.

Methods

One hundred and twenty-eight informative microsatellite loci distributed across all 29 linkage groups in Atlantic salmon were genotyped in individuals from four F2 families that were selected from the ends of the flesh colour distribution. Genotyping of 23 additional loci and two additional families was performed on a number of linkage groups harbouring putative QTL. QTL analysis was performed using a line-cross model assuming fixation of alternate QTL alleles and a half-sib model with no assumptions about the number and frequency of QTL alleles in the founder populations.

Results

A moderate to strong phenotypic correlation was found between colour, length and weight traits. In total, 13 genome-wide significant QTL were detected for all traits using the line-cross model, including three genome-wide significant QTL for flesh colour (Chr 6, Chr 26 and Chr 4). In addition, 32 suggestive QTL were detected (chromosome-wide P < 0.05). Using the half-sib model, six genome-wide significant QTL were detected for all traits, including two for flesh colour (Chr 26 and Chr 4) and 41 suggestive QTL were detected (chromosome-wide P < 0.05). Based on the half-sib analysis, these two genome-wide significant QTL for flesh colour explained 24% of the phenotypic variance for this trait.

Conclusions

A large number of significant and suggestive QTL for flesh colour and growth traits were found in an F2 population of Atlantic salmon. Chr 26 and Chr 4 presented the strongest evidence for significant QTL affecting flesh colour, while Chr 10, Chr 5, and Chr 4 presented the strongest evidence for significant QTL affecting growth traits (length and weight). These QTL could be strong candidates for use in marker-assisted selection and provide a starting point for further characterisation of the genetic components underlying flesh colour and growth.     

Validation of growth-related quantitative trait loci markers in different Exopalaemon carinicauda families for marker-assisted selection

We detected growth-related QTL and associated markers from the backcross population of Exopalaemon carinicauda in the previous study. Based on our previous study, the 47 SNP markers associated with candidate growth trait QTL were selected to analyze the association between these markers and body weight (BW), body length and abdominal segment length traits in four different populations including wild population, a full-sib family, a half-sib family and a backcross population for evaluating their potential application of marker-assisted selection in E. carinicauda. The general linear model (GLM) and mixed linear model were applied and the associations between SNP loci and three growth-related traits verified. The results showed that the Marker79268 and Marker100644 were significantly associated with the BW trait in more than three populations by the GLM method. The Marker100644 was significantly associated with BW in the full-sib family, half-sib family and backcross populations by the GLM and mixed linear model methods. Our findings will provide useful SNP markers to go forward to improve growth performance through more refined marker-assisted selection in E. carinicauda.     

Strategies for genetic mapping of categorical traits

The search for efficient and powerful statistical methods and optimal mapping strategies for categorical traits under various experimental designs continues to be one of the main tasks in genetic mapping studies. Methodologies for genetic mapping of categorical traits can generally be classified into two groups, linear and non-linear models. We develop a method based on a threshold model, termed mixture threshold model to handle ordinal (or binary) data from multiple families. Monte Carlo simulations are done to compare its statistical efficiencies and properties of the proposed non-linear model with a linear model for genetic mapping of categorical traits using multiple families. The mixture threshold model has notably higher statistical power than linear models. There may be an optimal sampling strategy (family size vs number of families) in which genetic mapping reaches its maximal power and minimal estimation errors. A single large-sibship family does not necessarily produce the maximal power for detection of quantitative trait loci (QTL) due to genetic sampling of QTL alleles. The QTL allelic model has a marked impact on efficiency of genetic mapping of categorical traits in terms of statistical power and QTL parameter estimation. Compared with a fixed number of QTL alleles (two or four), the model with an infinite number of QTL alleles and normally distributed allelic effects results in loss of statistical power. The results imply that inbred designs (e.g. F₂ or four-way crosses) with a few QTL alleles segregating or reducing number of QTL alleles (e.g. by selection) in outbred populations are desirable in genetic mapping of categorical traits using data from multiple families. This revised version was published online in July 2006 with corrections to the Cover Date.     

Selecting Informative Traits for Multivariate Quantitative Trait Locus Mapping Helps to Gain Optimal Power

A major consideration in multitrait analysis is which traits should be jointly analyzed. As a common strategy, multitrait analysis is performed either on pairs of traits or on all of traits. To fully exploit the power of multitrait analysis, we propose variable selection to choose a subset of informative traits for multitrait quantitative trait locus (QTL) mapping. The proposed method is very useful for achieving optimal statistical power for QTL identification and for disclosing the most relevant traits. It is also a practical strategy to effectively take advantage of multitrait analysis when the number of traits under consideration is too large, making the usual multivariate analysis of all traits challenging. We study the impact of selection bias and the usage of permutation tests in the context of variable selection and develop a powerful implementation procedure of variable selection for genome scanning. We demonstrate the proposed method and selection procedure in a backcross population, using both simulated and real data. The extension to other experimental mapping populations is straightforward.     

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.     

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.     

Combined line-cross and half-sib QTL analysis in Duroc–Pietrain population

A Duroc–Pietrain resource population was built to detect quantitative trait loci (QTL) that affect growth, carcass composition, and pork quality. The data were analyzed by applying three least-squares Mendelian models: a line-cross (LC) model, a half-sib (HS) model, and a combined LC and HS model (CB), which enabled the detection of QTL that had fixed, equal, and different allele frequencies for alternate breed alleles, respectively. Permutation tests were performed to determine 5% chromosome-wide and 5% genome-wide threshold values. A total of 40 (137) QTL were detected at the 5% genome-wide (chromosome-wide) level for the 35 traits analyzed. Of the 137 QTL detected, 62 were classified as the LC type (LC-QTL), 47 as the HS type (HS-QTL), and 28 as the CB type (CB-QTL). The results indicate that implementation of a series of model-based framework is not only beneficial to detect QTL, but also provides us with a new and more robust interpretation from which further methodology could be developed. G. Liu and J. J. Kim contributed equally to this work.     

Optimal sampling of a population to determine QTL location, variance, and allelic number

In a population intended for breeding and selection, questions of interest relative to a specific segregating QTL are the variance it generates in the population, and the number and effects of its alleles. One approach to address these questions is to extract several inbreds from the population and use them to generate multiple mapping families. Given random sampling of parents, sampling strategy may be an important factor determining the power of the analysis and its accuracy in estimating QTL variance and allelic number. We describe appropriate multiple-family QTL mapping methodology and apply it to simulated data sets to determine optimal sampling strategies in terms of family number versus family size. Genomes were simulated with seven chromosomes, on which 107 markers and six QTL were distributed. The total heritability was 0.60. Two to ten alleles were segregating at each QTL. Sampling strategies ranged from sampling two inbreds and generating a single family of 600 progeny to sampling 40 inbreds and generating 40 families of 15 progeny each. Strategies involving only one to five families were subject to variation due to the sampling of inbred parents. For QTL where more than two alleles were segregating, these strategies did not sample QTL alleles representative of the original population. Conversely, strategies involving 30 or more parents were subject to variation due to sampling of QTL genotypes within the small families obtained. Given these constraints, greatest QTL detection power was obtained for strategies involving five to ten mapping families. The most accurate estimation of the variance generated by the QTL, however, was obtained with strategies involving 20 or more families. Finally, strategies with an intermediate number of families best estimated the number of QTL alleles. We conclude that no overall optimal sampling strategy exists but that the strategy adopted must depend on the objective.Communicated by P. Langridge     

Combined linkage and linkage disequilibrium QTL mapping in multiple families of maize (Zea mays L.) line crosses highlights complementarities between models based on parental haplotype and single locus polymorphism

Advancements in genotyping are rapidly decreasing marker costs and increasing marker density. This opens new possibilities for mapping quantitative trait loci (QTL), in particular by combining linkage disequilibrium information and linkage analysis (LDLA). In this study, we compared different approaches to detect QTL for four traits of agronomical importance in two large multi-parental datasets of maize (Zea mays L.) of 895 and 928 testcross progenies composed of 7 and 21 biparental families, respectively, and genotyped with 491 markers. We compared to traditional linkage-based methods two LDLA models relying on the dense genotyping of parental lines with 17,728 SNP: one based on a clustering approach of parental line segments into ancestral alleles and one based on single marker information. The two LDLA models generally identified more QTL (60 and 52 QTL in total) than classical linkage models (49 and 44 QTL in total). However, they performed inconsistently over datasets and traits suggesting that a compromise must be found between the reduction of allele number for increasing statistical power and the adequacy of the model to potentially complex allelic variation. For some QTL, the model exclusively based on linkage analysis, which assumed that each parental line carried a different QTL allele, was able to capture remaining variation not explained by LDLA models. These complementarities between models clearly suggest that the different QTL mapping approaches must be considered to capture the different levels of allelic variation at QTL involved in complex traits.     

Simultaneous selection of major and minor genes: use of QTL to increase selection efficiency of coleoptile length of wheat (Triticum aestivum L.)

Plant breeders simultaneously select for qualitative traits controlled by one or a small number of major genes, as well as for polygenic traits controlled by multiple genes that may be detected as quantitative trait loci (QTL). In this study, we applied computer simulation to investigate simultaneous selection for alleles at both major and minor gene (as QTL) loci in breeding populations of two wheat parental lines, HM14BS and Sunstate. Loci targeted for selection included six major genes affecting plant height, disease resistance, and grain quality, plus 6 known and 11 “unidentified” QTL affecting coleoptile length (CL). Parental line HM14BS contributed the target alleles at two of the major gene loci, while parental line Sunstate contributed target alleles at four loci. The parents have similar plant height, but HM14BS has a longer coleoptile, a desirable attribute for deep sowing in rainfed environments. Including the wild-type allele at the major reduced-height locus Rht-D1, HM14BS was assumed to have 13 QTL for increased CL, and Sunstate four; these assumptions being derived from mapping studies and empirical data from an actual HM14BS/Sunstate population. Simulation indicated that compared to backcross populations, a single biparental F₁ cross produced the highest frequency of target genotypes (six desired alleles at major genes plus desired QTL alleles for long CL). From 1,000 simulation runs, an average of 2.4 individuals with the target genotype were present in unselected F₁-derived doubled haploid (DH) or recombinant inbred line (RIL) populations of size 200. A selection scheme for the six major genes increased the number of target individuals to 19.1, and additional marker-assisted selection (MAS) for CL increased the number to 23.0. Phenotypic selection (PS) of CL outperformed MAS in this study due to the high heritability of CL, incompletely linked markers for known QTL, and the existence of unidentified QTL. However, a selection scheme combining MAS and PS was equally as efficient as PS and would result in net savings in production and time to delivery of long coleoptile wheats containing the six favorable alleles.     

Quantitative trait loci variation for growth and obesity between and within lines of pigs (Sus scrofa)

The hypothesis that quantitative trait loci (QTL) that explain variation between divergent populations also account for genetic variation within populations was tested using pig populations. Two regions of the porcine genome that had previously been reported to harbor QTL with allelic effects that differed between the modern pig and its wild-type ancestor and between the modern pig and a more distantly related population of Asian pigs were studied. QTL for growth and obesity traits were mapped using selectively genotyped half-sib families from five domesticated modern populations. Strong support was found for at least one QTL segregating in each population. For all five populations there was evidence of a segregating QTL affecting fatness in a region on chromosome 7. These findings confirm that QTL can be detected in highly selected commercial populations and are consistent with the hypothesis that the same chromosome locations that account for variation between populations also explain genetic variation within populations.     

Use of the EM algorithm to detect QTL affecting multiple-traits in an across half-sib family analysis

QTL detection experiments in livestock species commonly use the half-sib design. Each male is mated to a number of females, each female producing a limited number of progeny. Analysis consists of attempting to detect associations between phenotype and genotype measured on the progeny. When family sizes are limiting experimenters may wish to incorporate as much information as possible into a single analysis. However, combining information across sires is problematic because of incomplete linkage disequilibrium between the markers and the QTL in the population. This study describes formulæ for obtaining MLEs via the expectation maximization (EM) algorithm for use in a multiple-trait, multiple-family analysis. A model specifying a QTL with only two alleles, and a common within sire error variance is assumed. Compared to single-family analyses, power can be improved up to fourfold with multi-family analyses. The accuracy and precision of QTL location estimates are also substantially improved. With small family sizes, the multi-family, multi-trait analyses reduce substantially, but not totally remove, biases in QTL effect estimates. In situations where multiple QTL alleles are segregating the multi-family analysis will average out the effects of the different QTL alleles.     

Combined line-cross and half-sib QTL analysis of crosses between outbred lines

Data from an F 2 cross between breeds of livestock are typically analysed by least squares line-cross or half-sib models to detect quantitative trait loci (QTL) that differ between or segregate within breeds. These models can also be combined to increase power to detect QTL, while maintaining the computational efficiency of least squares. Tests between models allow QTL to be characterized into those that are fixed (LC QTL), or segregating at similar (HS QTL) or different (CB QTL) frequencies in parental breeds. To evaluate power of the combined model, data wih various differences in QTL allele frequencies (FD) between parental breeds were simulated. Use of all models increased power to detect QTL. The line-cross model was the most powerful model to detect QTL for FD>0.6. The combined and half-sib models had similar power for FD<0.4. The proportion of detected QTL declared as LC QTL decreased with FD. The opposite was observed for HS QTL. The proportion of CB QTL decreased as FD deviated from 0.5. Accuracy of map position tended to be greatest for CB QTL. Models were applied to a cross of Berkshire and Yorkshire pig breeds and revealed 160 (40) QTL at the 5% chromosome (genome)-wise level for the 39 growth, carcass composition and quality traits, of which 72, 54, and 34 were declared as LC, HS and CB QTL. Fourteen CB QTL were detected only by the combined model. Thus, the combined model can increase power to detect QTL and mapping accuracy and enable characterization of QTL that segregate within breeds.     

Functional dissection of the Drosophila Kallmann's syndrome protein DmKal-1

Background

The proportion of muscle fibre types and their size affect muscularity as well as functional properties of the musculature and meat quality. We aimed to identify QTL for microstructural muscle properties including muscle fibre size, their numbers and fibre type proportions as well as biophysical parameters of meat quality and traits related to body composition, i.e. pH, conductivity, area of M. longissimus dorsi and lean meat content. A QTL scan was conducted in a porcine experimental population that is based on Duroc and Berlin Miniature Pig.

Results

Least square regression interval mapping revealed five significant and 42 suggestive QTL for traits related to muscle fibre composition under the line-cross model as well as eight significant and 40 suggestive QTL under the half-sib model. For traits related to body composition and biophysical parameters of meat quality five and twelve significant plus nine and 22 suggestive QTL were found under the line-cross and half-sib model, respectively. Regions with either significant QTL for muscle fibre traits or significant QTL for meat quality and muscularity or both were detected on SSC1, 2, 3, 4, 5, 13, 14, 15, and 16. QTL for microstructural properties explained a larger proportion of variance than did QTL for meat quality and body composition.

Conclusion

Microstructural properties of pig muscle and meat quality are governed by genetic variation at many loci distributed throughout the genome. QTL analysis under both, the line-cross and half-sib model, allows detecting QTL in case of fixation or segregation of the QTL alleles among the founder populations and thus provide comprehensive insight into the genetic variation of the traits under investigation. Genomic regions affecting complex traits of muscularity and meat quality as well as microstructural properties might point to QTL that in first instance affect muscle fibre traits and by this in second instance meat quality. Disentangling complex traits in their constituent phenotypes might facilitate the identification of QTL and the elucidation of the pleiotropic nature of QTL effects.     

Connected populations for detecting quantitative trait loci and testing for epistasis: an application in maize

Quantitative trait loci (QTL) detection experiments have often been restricted to large biallelic populations. Use of connected multiparental crosses has been proposed to increase the genetic variability addressed and to test for epistatic interactions between QTL and the genetic background. We present here the results of a QTL detection performed on six connected F₂ populations of 150 F_2:3 families each, derived from four maize inbreds and evaluated for three traits of agronomic interest. The QTL detection was carried out by composite interval mapping on each population separately, then on the global design either by taking into account the connections between populations or not. Epistatic interactions between loci and with the genetic background were tested. Taking into account the connections between populations increased the number of QTL detected and the accuracy of QTL position estimates. We detected many epistatic interactions, particularly for grain yield QTL (R ² increase of 9.6%). Use of connections for the QTL detection also allowed a global ranking of alleles at each QTL. Allelic relationships and epistasis both contribute to the lack of consistency for QTL positions observed among populations, in addition to the limited power of the tests. The potential benefit of assembling favorable alleles by marker-assisted selection are discussed.     

The number of loci that affect milk production traits in dairy cattle

We have used the results of an experiment mapping quantitative trait loci (QTL) affecting milk yield and composition to estimate the total number of QTL affecting these traits. We did this by estimating the number of segregating QTL within a half-sib daughter design using logic similar to that used to estimate the "false discovery rate" (FDR). In a half-sib daughter design with six sire families we estimate that the average sire was heterozygous for approximately 5 QTL per trait. Also, in most cases only one sire was heterozygous for any one QTL; therefore at least 30 QTL were likely to be segregating for these milk production traits in this Holstein population.     

Pleiotropic quantitative trait loci contribute to population divergence in traits associated with life-history variation in Mimulus guttatus

Evolutionary biologists seek to understand the genetic basis for multivariate phenotypic divergence. We constructed an F2 mapping population (N = 539) between two distinct populations of Mimulus guttatus. We measured 20 floral, vegetative, and life-history characters on parents and F1 and F2 hybrids in a common garden experiment. We employed multitrait composite interval mapping to determine the number, effect, and degree of pleiotropy in quantitative trait loci (QTL) affecting divergence in floral, vegetative, and life-history characters. We detected 16 QTL affecting floral traits; 7 affecting vegetative traits; and 5 affecting selected floral, vegetative, and life-history traits. Floral and vegetative traits are clearly polygenic. We detected a few major QTL, with all remaining QTL of small effect. Most detected QTL are pleiotropic, implying that the evolutionary shift between these annual and perennial populations is constrained. We also compared the genetic architecture controlling floral trait divergence both within (our intraspecific study) and between species, on the basis of a previously published analysis of M. guttatus and M. nasutus. Eleven of our 16 floral QTL map to approximately the same location in the interspecific map based on shared, collinear markers, implying that there may be a shared genetic basis for floral divergence within and among species of Mimulus.     

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.     

Methodology and Accuracy of Estimation of Quantitative Trait Loci Parameters in a Half-Sib Design Using Maximum Likelihood

Maximum likelihood methods were developed for estimation of the six parameters relating to a marker-linked quantitative trait locus (QTL) segregating in a half-sib design, namely the QTL additive effect, the QTL dominance effect, the population mean, recombination between the marker and the QTL, the population frequency of the QTL alleles, and the within-family residual variance. The method was tested on simulated stochastic data with various family structures under two genetic models. A method for predicting the expected value of the likelihood was also derived and used to predict the lower bound sampling errors of the parameter estimates and the correlations between them. It was found that standard errors and confidence intervals were smallest for the population mean and variance, intermediate for QTL effects and allele frequency, and highest for recombination rate. Correlations among standard errors of the parameter estimates were generally low except for a strong negative correlation (r = -0.9) between the QTL's dominance effect and the population mean, and medium positive and negative correlations between the QTL's additive effect and, respectively, recombination rate (r = 0.5) and residual variance (r = -0.6). The implications for experimental design and method of analysis on power and accuracy of marker-QTL linkage experiments were discussed.