Power of mixed-model QTL mapping from phenotypic, pedigree and marker data in self-pollinated crops

The power of QTL mapping by a mixed-model approach has been studied for hybrid crops but remains unknown in self-pollinated crops. Our objective was to evaluate the usefulness of mixed-model QTL mapping in the context of a breeding program for a self-pollinated crop. Specifically, we simulated a soybean (Glycine max L. Merr.) breeding program and applied a mixed-model approach that comprised three steps: variance component estimation, single-marker analyses, and multiple-marker analysis. Average power to detect QTL ranged from <1 to 47% depending on the significance level (0.01 or 0.0001), number of QTL (20 or 80), heritability of the trait (0.40 or 0.70), population size (600 or 1,200 inbreds), and number of markers (300 or 600). The corresponding false discovery rate ranged from 2 to 43%. Larger populations, higher heritability, and fewer QTL controlling the trait led to a substantial increase in power and to a reduction in the false discovery rate and bias. A stringent significance level reduced both the power and false discovery rate. There was greater power to detect major QTL than minor QTL. Power was higher and the false discovery rate was lower in hybrid crops than in self-pollinated crops. We conclude that mixed-model QTL mapping is useful for gene discovery in plant breeding programs of self-pollinated crops.     

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.     

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.     

On the use of mathematically-derived traits in QTL mapping

Mathematically-derived traits from two or more component traits, either by addition, subtraction, multiplication, or division, have been frequently used in genetics and breeding. When used in quantitative trait locus (QTL) mapping, derived traits sometimes show discrepancy with QTL identified for the component traits. We used three QTL distributions and three genetic effects models, and an actual maize mapping population, to investigate the efficiency of using derived traits in QTL mapping, and to understand the genetic and biological basis of derived-only QTL, i.e., QTL identified for a derived trait but not for any component trait. Results indicated that the detection power of the four putative QTL was consistently greater than 90% for component traits in simulated populations, each consisting of 200 recombinant inbred lines. Lower detection power and higher false discovery rate (FDR) were observed when derived traits were used. In an actual maize population, simulations were designed based on the observed QTL distributions and effects. When derived traits were used, QTL detected for both component and derived traits had comparable power, but those detected for component traits but not for derived traits had low detection power. The FDR from subtraction and division in the maize population were higher than the FDR from addition and multiplication. The use of derived traits increased the gene number, caused higher-order gene interactions than observed in component traits, and possibly complicated the linkage relationship between QTL as well. The increased complexity of the genetic architecture with derived traits may be responsible for the reduced detection power and the increased FDR. Derived-only QTL identified in practical genetic populations can be explained either as minor QTL that are not significant in QTL mapping of component traits, or as false positives.     

Genomewide markers as cofactors for precision mapping of quantitative trait loci

In composite interval mapping of quantitative trait loci (QTL), subsets of background markers are used to account for the effects of QTL outside the marker interval being tested. Here, I propose a QTL mapping approach (called G model) that utilizes genomewide markers as cofactors. The G model involves backward elimination on a given chromosome after correcting for genomewide marker effects, calculated under a random effects model, at all the other chromosomes. I simulated a trait controlled by 15 or 30 QTL, mapping populations of N = 96, 192, and 384 recombinant inbreds, and N _M = 192 and 384 evenly spaced markers. In the C model, which utilized subsets of background markers, the number of QTL detected and the number of false positives depended on the number of cofactors used, with five cofactors being too few with N = 384 and 20–40 cofactors being too many with N = 96. A window size of 0 cM for excluding cofactors maintained the number of true QTL detected while decreasing the number of false positives. The number of true QTL detected was generally higher with the G model than with the C model, and the G model led to good control of the type I error rate in simulations where the null hypothesis of no marker–QTL linkage was true. Overall, the results indicated that the G model is useful in QTL mapping because it is less subjective and has equal, if not better, performance when compared with the traditional approach of using subsets of markers to account for background QTL.     

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.     

Power of quantitative trait locus mapping for polygenic binary traits using generalized and regression interval mapping in multi-family half-sib designs

A generalized interval mapping (GIM) method to map quantitative trait loci (QTL) for binary polygenic traits in a multi-family half-sib design is developed based on threshold theory and implemented using a Newton-Raphson algorithm. Statistical power and bias of QTL mapping for binary traits by GIM is compared with linear regression interval mapping (RIM) using simulation. Data on 20 paternal half-sib families were simulated with two genetic markers that bracketed an additive QTL. Data simulated and analysed were: (1) data on the underlying normally distributed liability (NDL) scale, (2) binary data created by truncating NDL data based on three thresholds yielding data sets with three different incidences, and (3) NDL data with polygenic and QTL effects reduced by a proportion equal to the ratio of the heritabilities on the binary versus NDL scale (reduced-NDL). Binary data were simulated with and without systematic environmental (herd) effects in an unbalanced design. GIM and RIM gave similar power to detect the QTL and similar estimates of QTL location, effects and variances. Presence of fixed effects caused differences in bias between RIM and GIM, where GIM showed smaller bias which was affected less by incidence. The original NDL data had higher power and lower bias in QTL parameter estimates than binary and reduced-NDL data. RIM for reduced-NDL and binary data gave similar power and estimates of QTL parameters, indicating that the impact of the binary nature of data on QTL analysis is equivalent to its impact on heritability.     

远交群体动态性状基因定位的似然分析Ⅱ.模拟研究

以远交设计群体为例,在推导出动态性状基因定位的似然法分析过程的基础上,选择3阶Legendre多项式为子模型,采用Monte-Carlo方法模拟研究了不同个体数、测定日频数、标记密度和QTL遗传贡献率对两种分析方法检测QTL效率的影响。每个因素都取高、中和低3个水平,利用正交设计安排模拟因素试验组合。模拟试验结果表明：高QTL遗传贡献率要比低QTL遗传贡献率的QTL在检测时需要较少个体数和测定日抽样;但不论QTL遗传贡献率多大,300以上的群体大小和5％以上的测定日频数都可以保证足够高的检测效率。个体数和测定日频数对动态性状QTL的分析和检测具有几乎相同的作用,而且相同样本含量条件下两者呈现互补的关系。就某个动态点的QTL检测而言,模拟试验也同时证明：提出的这种以整个动态过程为定位目标的动态性状基因定位方法明显优于传统的逐个动态点的定位分析方法。     

Enhanced efficiency of quantitative trait loci mapping analysis based on multivariate complexes of quantitative traits

An approach to increase the efficiency of mapping quantitative trait loci (QTL) was proposed earlier by the authors on the basis of bivariate analysis of correlated traits. The power of QTL detection using the log-likelihood ratio (LOD scores) grows proportionally to the broad sense heritability. We found that this relationship holds also for correlated traits, so that an increased bivariate heritability implicates a higher LOD score, higher detection power, and better mapping resolution. However, the increased number of parameters to be estimated complicates the application of this approach when a large number of traits are considered simultaneously. Here we present a multivariate generalization of our previous two-trait QTL analysis. The proposed multivariate analogue of QTL contribution to the broad-sense heritability based on interval-specific calculation of eigenvalues and eigenvectors of the residual covariance matrix allows prediction of the expected QTL detection power and mapping resolution for any subset of the initial multivariate trait complex. Permutation technique allows chromosome-wise testing of significance for the whole trait complex and the significance of the contribution of individual traits owing to: (a) their correlation with other traits, (b) dependence on the chromosome in question, and (c) both a and b. An example of application of the proposed method on a real data set of 11 traits from an experiment performed on an F(2)/F(3) mapping population of tetraploid wheat (Triticum durum x T. dicoccoides) is provided.     

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.     

家畜抗性等级性状的QTL定位方法

在广义线性模型的框架内模拟研究了家畜抗性等级性状的QTL定位方法,QTL参数的估计采用最大似然方法,比较了阈模型方法与一般线性方法的QTL定位效率,并对影响等级性状QTL定位效率的主要因素（QTL效应、性状的遗传力）进行了模拟研究,实验设计为多个家系的女儿设计,资源群体大小为500头。研究结果表明：在QTL位置参数估计及检验功效方面,阈模型方法具有一定的优势,对抗性等级性状QTL定位的功效也高于线性方法。另外,性状遗传力和QTL效应的大小对QTL定位的准确度也有直接的影响,随着性状遗传力QTL效应的     

Power of QTL detection by either fixed or random models in half-sib designs

The aim of this study was to compare the variance component approach for QTL linkage mapping in half-sib designs to the simple regression method. Empirical power was determined by Monte Carlo simulation in granddaughter designs. The factors studied (base values in parentheses) included the number of sires (5) and sons per sire (80), ratio of QTL variance to total genetic variance (λ = 0.1), marker spacing (10 cM), and QTL allele frequency (0.5). A single bi-allelic QTL and six equally spaced markers with six alleles each were simulated. Empirical power using the regression method was 0.80, 0.92 and 0.98 for 5, 10, and 20 sires, respectively, versus 0.88, 0.98 and 0.99 using the variance component method. Power was 0.74, 0.80, 0.93, and 0.95 using regression versus 0.77, 0.88, 0.94, and 0.97 using the variance component method for QTL variance ratios (λ) of 0.05, 0.1, 0.2, and 0.3, respectively. Power was 0.79, 0.85, 0.80 and 0.87 using regression versus 0.80, 0.86, 0.88, and 0.85 using the variance component method for QTL allele frequencies of 0.1, 0.3, 0.5, and 0.8, respectively. The log₁₀ of type I error profiles were quite flat at close marker spacing (1 cM), confirming the inability to fine-map QTL by linkage analysis in half-sib designs. The variance component method showed slightly more potential than the regression method in QTL mapping.     

Efficiency of genome-wide association studies in random cross populations

Genome-wide association studies (GWAS) with plant species have employed inbred lines panels. We evaluated the efficiency of GWAS in non-inbred and inbred populations and assessed factors affecting GWAS. Fifty samples of 800 individuals from populations with linkage disequilibrium were simulated. Individuals were genotyped for 10,000 single nucleotide polymorphisms (SNPs) and phenotyped for traits controlled by ten quantitative trait loci (QTLs) and 90 minor genes, assuming different degrees of dominance and broad sense heritabilities of 40 and 80%. The average SNP density was 0.1 centiMorgan (cM) and the QTL heritabilities ranged from 3.2 to 11.8%. The results for random cross populations evidenced that to increase the QTL detection power, the additive-dominance model must be fitted for traits controlled by dominance effects but must not be fitted for traits showing no dominance. The power of detection was maximized by increasing the sample size to 400 and the false discovery rate (FDR) to 5%. The average power of detection for the low, intermediate, and high heritability QTLs achieved 52.4, 87.0, and 100.0%, respectively. Assuming sample sizes of 400 and 800, the observed FDR was equal to or lower than the specified level of significance. The association mapping was highly precise, since at least 97% of the declared QTLs were detected by the SNP inside it (average bias of 0.4 cM). Besides controlling the FDR, relatedness (and identity by state) efficiently controls the number of significant associations outside the QTL interval (not all false positive associations). The analysis of the inbred random cross population provided essentially the same results as the non-inbred populations.     

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.     

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.     

Power of in silico QTL mapping from phenotypic, pedigree, and marker data in a hybrid breeding program

Most quantitative trait locus (QTL) mapping studies in plants have used designed mapping populations. As an alternative to traditional QTL mapping, in silico mapping via a mixed-model approach simultaneously exploits phenotypic, genotypic, and pedigree data already available in breeding programs. The statistical power of this in silico mapping method, however, remains unknown. Our objective was to evaluate the power of in silico mapping via a mixed-model approach in hybrid crops. We used maize (Zea mays L.) as a model species to study, by computer simulation, the influence of number of QTLs (20 or 80), heritability (0.40 or 0.70), number of markers (200 or 400), and sample size (600 or 2,400 hybrids). We found that the average power to detect QTLs ranged from 0.11 to 0.59 for a significance level of =0.01, and from 0.01 to 0.47 for =0.0001. The false discovery rate ranged from 0.22 to 0.74 for =0.01, and from 0.05 to 0.46 for =0.0001. As with designed mapping experiments, a large sample size, high marker density, high heritability, and small number of QTLs led to the highest power for in silico mapping via a mixed-model approach. The power to detect QTLs with large effects was greater than the power to detect QTL with small effects. We conclude that gene discovery in hybrid crops can be initiated by in silico mapping. Finding an acceptable compromise, however, between the power to detect QTL and the proportion of false QTL would be necessary.     

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.     

孙女设计中标记密度对QTL定位精确性的影响

采用蒙特卡罗方法分析了在孙女设计中不同的嫩体结构、性状遗传力、ＱＴＬ效应大小和ＱＴＬ在染色体上的位置中个因素不同水平组合下４种标记密度（标记间隔５ｃＭ,１０ｃＭ,２０ｃＭ、５０ｃＭ对ＱＴＬ定位精确性（以均方误ＭＳＥ为衡量指标）的影响,并从经济学角度探讨了应用于标记辅助选（ＭＡＳ）的ＱＴＬ定位的最佳标记密度。结果表明,一般说来,在各因素水平都较低时,ＭＳＥ随标记密度加大而下降的相对幅度也较 小,反之     

Mapping of multiple quantitative trait loci affecting bovine spongiform encephalopathy

A whole-genome scan was conducted to map quantitative trait loci (QTL) for BSE resistance or susceptibility. Cows from four half-sib families were included and 173 microsatellite markers were used to construct a 2835-cM (Kosambi) linkage map covering 29 autosomes and the pseudoautosomal region of the sex chromosome. Interval mapping by linear regression was applied and extended to a multiple-QTL analysis approach that used identified QTL on other chromosomes as cofactors to increase mapping power. In the multiple-QTL analysis, two genome-wide significant QTL (BTA17 and X/Y(ps)) and four genome-wide suggestive QTL (BTA1, 6, 13, and 19) were revealed. The QTL identified here using linkage analysis do not overlap with regions previously identified using TDT analysis. One factor that may explain the disparity between the results is that a more extensive data set was used in the present study. Furthermore, methodological differences between TDT and linkage analyses may affect the power of these approaches.     

离散性状QTL区间定位的最大似然方法

采用最大似然区间定位法对阈模型与一般线性模型的QTL定位效率进行了比较,并对影响离散性状QTL检测效率的主要因素（QTL效应、性状的遗传力和表型发生率）进行了模拟研究,实验设计为多个家系的女儿设计．资源群体大小为500头。研究结果表明：在QTL参数估计及检验功效方面,阈模型方法具有较大的优势,对离散性状QTL定位的效率明显高于LM（Linear Model）方法,定位的准确性也较高。另外,性状遗传力、QTL效应的大小和性状表型发生率对QTL定位的准确度也有直接的影响,随着性状遗传力和表型发生率的提高,随着QTL效应的增大,QTL定位的效率也进一步提高。