Simple deterministic identity-by-descent coefficients and estimation of QTL allelic effects in full and half sibs

Accurate and rapid methods for the detection of quantitative trait loci (QTLs) and evaluation of consequent allelic effects are required to implement marker-assisted selection in outbred populations. In this study, we present a simple deterministic method for estimating identity-by-descent (IBD) coefficients in full- and half-sib families that can be used for the detection of QTLs via a variance-component approach. In a simulated dataset, IBD coefficients among sibs estimated by the simple deterministic and Markov chain Monte Carlo (MCMC) methods with three or four alleles at each marker locus exhibited a correlation of greater than 0.99. This high correlation was also found in QTL analyses of data from an outbred pig population. Variance component analysis used both the simple deterministic and MCMC methods to estimate IBD coefficients. Both procedures detected a QTL at the same position and gave similar test statistics and heritabilities. The MCMC method, however, required much longer computation than the simple method. The conversion of estimated QTL genotypic effects into allelic effects for use in marker-assisted selection is also demonstrated.     

QTL analysis in arbitrary pedigrees with incomplete marker information

Mapping quantitative trait loci (QTL) in arbitrary outbred pedigrees is complicated by the combinatorial possibilities of allele flow relationships and of the founder allelic configurations. Exact methods are only available for rather short and simple pedigrees. Stochastic simulation using Markov chain Monte Carlo (MCMC) integration offers more flexibility. MCMC methods are less natural in a frequentist than in a Bayesian context, which we therefore adopt. Among the MCMC algorithms for updating marker locus genotypes, we implement the descent-graph algorithm. It can be used to update marker locus allele flow relationships and can handle arbitrarily complex pedigrees and missing marker information. Compared with updating marker genotypic information, updating QTL parameters, such as position, effects, and the allele flow relationships is relatively easy with MCMC. We treat the effect of each diploid combination of founder alleles as a random variable and only estimate the variance of these effects, ie, we model diploid genotypic effects instead of the usual partition in additive and dominance effects. This is a variant of the random model approach. The number of QTL alleles is generally unknown. In the Bayesian context, the number of QTL present on a linkage group can be treated as variable. Computer simulations suggest that the algorithm can indeed handle complex pedigrees and detect two QTL on a linkage group, but that the number of individuals in a single extended family is limited to about 50 to 100 individuals.     

Detection of multiple QTL with epistatic effects under a mixed inheritance model in an outbred population

A quantitative trait depends on multiple quantitative trait loci (QTL) and on the interaction between two or more QTL, named epistasis. Several methods to detect multiple QTL in various types of design have been proposed, but most of these are based on the assumption that each QTL works independently and epistasis has not been explored sufficiently. The objective of the study was to propose an integrated method to detect multiple QTL with epistases using Bayesian inference via a Markov chain Monte Carlo (MCMC) algorithm. Since the mixed inheritance model is assumed and the deterministic algorithm to calculate the probabilities of QTL genotypes is incorporated in the method, this can be applied to an outbred population such as livestock. Additionally, we treated a pair of QTL as one variable in the Reversible jump Markov chain Monte Carlo (RJMCMC) algorithm so that two QTL were able to be simultaneously added into or deleted from a model. As a result, both of the QTL can be detected, not only in cases where either of the two QTL has main effects and they have epistatic effects between each other, but also in cases where neither of the two QTL has main effects but they have epistatic effects. The method will help ascertain the complicated structure of quantitative traits.     

Bayesian mapping of quantitative trait loci under complicated mating designs

Quantitative trait loci (QTL) are easily studied in a biallelic system. Such a system requires the cross of two inbred lines presumably fixed for alternative alleles of the QTL. However, development of inbred lines can be time consuming and cost ineffective for species with long generation intervals and severe inbreeding depression. In addition, restriction of the investigation to a biallelic system can sometimes be misleading because many potentially important allelic interactions do not have a chance to express and thus fail to be detected. A complicated mating design involving multiple alleles mimics the actual breeding system. However, it is difficult to develop the statistical model and algorithm using the classical maximum-likelihood method. In this study, we investigate the application of a Bayesian method implemented via the Markov chain Monte Carlo (MCMC) algorithm to QTL mapping under arbitrarily complicated mating designs. We develop the method under a mixed-model framework where the genetic values of founder alleles are treated as random and the nongenetic effects are treated as fixed. With the MCMC algorithm, we first draw the gene flows from the founders to the descendants for each QTL and then draw samples of the genetic parameters. Finally, we are able to simultaneously infer the posterior distribution of the number, the additive and dominance variances, and the chromosomal locations of all identified QTL.     

A random model approach to mapping quantitative trait loci for complex binary traits in outbred populations.

Mapping quantitative trait loci (QTL) for complex binary traits is more challenging than for normally distributed traits due to the nonlinear relationship between the observed phenotype and unobservable genetic effects, especially when the mapping population contains multiple outbred families. Because the number of alleles of a QTL depends on the number of founders in an outbred population, it is more appropriate to treat the effect of each allele as a random variable so that a single variance rather than individual allelic effects is estimated and tested. Such a method is called the random model approach. In this study, we develop the random model approach of QTL mapping for binary traits in outbred populations. An EM-algorithm with a Fisher-scoring algorithm embedded in each E-step is adopted here to estimate the genetic variances. A simple Monte Carlo integration technique is used here to calculate the likelihood-ratio test statistic. For the first time we show that QTL of complex binary traits in an outbred population can be scanned along a chromosome for their positions, estimated for their explained variances, and tested for their statistical significance. Application of the method is illustrated using a set of simulated data.     

Study on mapping Quantitative Trait Loci for animal complex binary traits using Bayesian-Markov chain Monte Carlo approach

It is a challenging issue to map Quantitative Trait Loci (QTL) underlying complex discrete traits, which usually show discontinuous distribution and less information, using conventional statistical methods. Bayesian-Markov chain Monte Carlo (Bayesian-MCMC) approach is the key procedure in mapping QTL for complex binary traits, which provides a complete posterior distribution for QTL parameters using all prior information. As a consequence, Bayesian estimates of all interested variables can be obtained straightforwardly basing on their posterior samples simulated by the MCMC algorithm. In our study, utilities of Bayesian-MCMC are demonstrated using simulated several animal outbred full-sib families with different family structures for a complex binary trait underlied by both a QTL and polygene. Under the Identity-by-Descent-Based variance component random model, three samplers basing on MCMC, including Gibbs sampling, Metropolis algorithm and reversible jump MCMC, were implemented to generate the joint posterior distribution of all unknowns so that the QTL parameters were obtained by Bayesian statistical inferring. The results showed that Bayesian-MCMC approach could work well and robust under different family structures and QTL effects. As family size increases and the number of family decreases, the accuracy of the parameter estimates will be improved. When the true QTL has a small effect, using outbred population experiment design with large family size is the optimal mapping strategy.     

Mapping quantitative trait loci in F2 incorporating phenotypes of F3 progeny

In plants and laboratory animals, QTL mapping is commonly performed using F(2) or BC individuals derived from the cross of two inbred lines. Typical QTL mapping statistics assume that each F(2) individual is genotyped for the markers and phenotyped for the trait. For plant traits with low heritability, it has been suggested to use the average phenotypic values of F(3) progeny derived from selfing F(2) plants in place of the F(2) phenotype itself. All F(3) progeny derived from the same F(2) plant belong to the same F(2:3) family, denoted by F(2:3). If the size of each F(2:3) family (the number of F(3) progeny) is sufficiently large, the average value of the family will represent the genotypic value of the F(2) plant, and thus the power of QTL mapping may be significantly increased. The strategy of using F(2) marker genotypes and F(3) average phenotypes for QTL mapping in plants is quite similar to the daughter design of QTL mapping in dairy cattle. We study the fundamental principle of the plant version of the daughter design and develop a new statistical method to map QTL under this F(2:3) strategy. We also propose to combine both the F(2) phenotypes and the F(2:3) average phenotypes to further increase the power of QTL mapping. The statistical method developed in this study differs from published ones in that the new method fully takes advantage of the mixture distribution for F(2:3) families of heterozygous F(2) plants. Incorporation of this new information has significantly increased the statistical power of QTL detection relative to the classical F(2) design, even if only a single F(3) progeny is collected from each F(2:3) family. The mixture model is developed on the basis of a single-QTL model and implemented via the EM algorithm. Substantial computer simulation was conducted to demonstrate the improved efficiency of the mixture model. Extension of the mixture model to multiple QTL analysis is developed using a Bayesian approach. The computer program performing the Bayesian analysis of the simulated data is available to users for real data analysis.     

Bayesian mapping of multiple quantitative trait loci from incomplete outbred offspring data

A general fine-scale Bayesian quantitative trait locus (QTL) mapping method for outcrossing species is presented. It is suitable for an analysis of complete and incomplete data from experimental designs of F2 families or backcrosses. The amount of genotyping of parents and grandparents is optional, as well as the assumption that the QTL alleles in the crossed lines are fixed. Grandparental origin indicators are used, but without forgetting the original genotype or allelic origin information. The method treats the number of QTL in the analyzed chromosome as a random variable and allows some QTL effects from other chromosomes to be taken into account in a composite interval mapping manner. A block-update of ordered genotypes (haplotypes) of the whole family is sampled once in each marker locus during every round of the Markov Chain Monte Carlo algorithm used in the numerical estimation. As a byproduct, the method gives the posterior distributions for linkage phases in the family and therefore it can also be used as a haplotyping algorithm. The Bayesian method is tested and compared with two frequentist methods using simulated data sets, considering two different parental crosses and three different levels of available parental information. The method is implemented as a software package and is freely available under the name Multimapper/outbred at URL http://www.rni.helsinki.fi/mjs/.     

Using dominance relationship coefficients based on linkage disequilibrium and linkage with a general complex pedigree to increase mapping resolution

Dominance (intralocus allelic interactions) plays often an important role in quantitative trait variation. However, few studies about dominance in QTL mapping have been reported in outbred animal or human populations. This is because common dominance effects can be predicted mainly for many full sibs, which do not often occur in outbred or natural populations with a general pedigree. Moreover, incomplete genotypes for such a pedigree make it infeasible to estimate dominance relationship coefficients between individuals. In this study, identity-by-descent (IBD) coefficients are estimated on the basis of population-wide linkage disequilibrium (LD), which makes it possible to track dominance relationships between unrelated founders. Therefore, it is possible to use dominance effects in QTL mapping without full sibs. Incomplete genotypes with a complex pedigree and many markers can be efficiently dealt with by a Markov chain Monte Carlo method for estimating IBD and dominance relationship matrices (D(RM)). It is shown by simulation that the use of D(RM) increases the likelihood ratio at the true QTL position and the mapping accuracy and power with complete dominance, overdominance, and recessive inheritance modes when using 200 genotyped and phenotyped individuals.     

Study on mapping Quantitative Trait Loci for animal complex binary traits using Bayesian-Markov chain Monte Carlo approach

It is a challenging issue to map Quantitative Trait Loci (QTL) underlying complex discrete traits,which usually show discontinuous distribution and less information,using conventional statisti-cal methods. Bayesian-Markov chain Monte Carlo (Bayesian-MCMC) approach is the key procedure in mapping QTL for complex binary traits,which provides a complete posterior distribution for QTL parameters using all prior information. As a consequence,Bayesian estimates of all interested vari-ables can be obtained straightforwardly basing on their posterior samples simulated by the MCMC algorithm. In our study,utilities of Bayesian-MCMC are demonstrated using simulated several ani-mal outbred full-sib families with different family structures for a complex binary trait underlied by both a QTL and polygene. Under the Identity-by-Descent-Based variance component random model,three samplers basing on MCMC,including Gibbs sampling,Metropolis algorithm and reversible jump MCMC,were implemented to generate the joint posterior distribution of all unknowns so that the QTL parameters were obtained by Bayesian statistical inferring. The results showed that Bayesian-MCMC approach could work well and robust under different family structures and QTL effects. As family size increases and the number of family decreases,the accuracy of the parameter estimates will be im-proved. When the true QTL has a small effect,using outbred population experiment design with large family size is the optimal mapping strategy.     

QTL detection and allelic effects for growth and fat traits in outbred pig populations

Quantitative trait loci (QTL) for growth and fatness traits have previously been identified on chromosomes 4 and 7 in several experimental pig populations. The segregation of these QTL in commercial pigs was studied in a sample of 2713 animals from five different populations. Variance component analysis (VCA) using a marker-based identity by descent (IBD) matrix was applied. The IBD coefficient was estimated with simple deterministic (SMD) and Markov chain Monte Carlo (MCMC) methods. Data for two growth traits, average daily gain on test and whole life daily gain, and back fat thickness were analysed. With both methods, seven out of 26 combinations of population, chromosome and trait, were significant. Additionally, QTL genotypic and allelic effects were estimated when the QTL effect was significant. The range of QTL genotypic effects in a population varied from 4.8% to 10.9% of the phenotypic mean for growth traits and 7.9% to 19.5% for back fat trait. Heritabilities of the QTL genotypic values ranged from 8.6% to 18.2% for growth traits, and 14.5% to 19.2% for back fat. Very similar results were obtained with both SMD and MCMC. However, the MCMC method required a large number of iterations, and hence computation time, especially when the QTL test position was close to the marker.     

A method for computing identity by descent probabilities and quantitative trait loci mapping with dominant (AFLP) markers

Amplified fragment length polymorphisms (AFLPs) are a widely used marker system: the technique is very cost-effective, easy and rapid, and reproducibly generates hundreds of markers. Unfortunately, AFLP alleles are typically scored as the presence or absence of a band and, thus, heterozygous and dominant homozygous genotypes cannot be distinguished. This results in a significant loss of information, especially as regards mapping of quantitative trait loci (QTLs). We present a Monte Carlo Markov Chain method that allows us to compute the identity by descent probabilities (IBD) in a general pedigree whose individuals have been typed for dominant markers. The method allows us to include the information provided by the fluorescent band intensities of the markers, the rationale being that homozygous individuals have on average higher band intensities than heterozygous individuals, as well as information from linked markers in each individual and its relatives. Once IBD probabilities are obtained, they can be combined into the QTL mapping strategy of choice. We illustrate the method with two simulated populations: an outbred population consisting of full sib families, and an F2 cross between inbred lines. Two marker spacings were considered, 5 or 20 cM, in the outbred population. There was almost no difference, for the practical purpose of QTL estimation, between AFLPs and biallelic codominant markers when the band density is taken into account, especially at the 5 cM spacing. The performance of AFLPs every 5 cM was also comparable to that of highly polymorphic markers (microsatellites) spaced every 20 cM. In economic terms, QTL mapping with a dense map of AFLPs is clearly better than microsatellite QTL mapping and little is lost in terms of accuracy of position. Nevertheless, at low marker densities, AFLPs or other biallelic markers result in very inaccurate estimates of QTL position.     

QTL mapping in outbred tetraploid (and diploid) diallel populations

Over the last decade, multiparental populations have become a mainstay of genetics research in diploid species. Our goal was to extend this paradigm to autotetraploids by developing software for quantitative trait locus (QTL) mapping in connected F1 populations derived from a set of shared parents. For QTL discovery, phenotypes are regressed on the dosage of parental haplotypes to estimate additive effects. Statistical properties of the model were explored by simulating half-diallel diploid and tetraploid populations with different population sizes and numbers of parents. Across scenarios, the number of progeny per parental haplotype (pph) largely determined the statistical power for QTL detection and accuracy of the estimated haplotype effects. Multiallelic QTL with heritability 0.2 were detected with 90% probability at 25 pph and genome-wide significance level 0.05, and the additive haplotype effects were estimated with over 90% accuracy. Following QTL discovery, the software enables a comparison of models with multiple QTL and nonadditive effects. To illustrate, we analyzed potato tuber shape in a half-diallel population with three tetraploid parents. A well-known QTL on chromosome 10 was detected, for which the inclusion of digenic dominance lowered the Deviance Information Criterion (DIC) by 17 points compared to the additive model. The final model also contained a minor QTL on chromosome 1, but higher-order dominance and epistatic effects were excluded based on the DIC. In terms of practical impacts, the software is already being used to select offspring based on the effect and dosage of particular haplotypes in breeding programs.     

Efficient Method for Analysis of QTL Using F1 Progenies in an Outcrossing Species

In the present paper, we proposed a statistical procedure based on composite interval mapping for accurate analysis of quantitative trait loci (QTL) for individuals sampled from an outcrossing population with two-generation families consisting of the sampled individuals and F1 progenies obtained by crossing them as parental individuals. In the proposed procedure, haplotypes of markers of parental individuals were reconstructed based on the genotypes of F1 progenies and QTL analyses with composite interval mapping were conducted separately for each of parents as well as jointly for both parents. A least squares method was applied to the composite interval mapping, where some of markers were selected as cofactors to absorb the variation induced by QTL located elsewhere in the genome. The procedure was evaluated for the efficiency in detecting QTL and the precision of estimates of locations and effects of QTL using simulations. It was shown that QTL with interaction between paternal and maternal alleles was effectively detected by joint analysis of both parents, while a QTL segregating only in one parent, closely linked to a QTL segregating only in the other parent, was successfully detected by analyzing separately each of the parents with inclusion of markers of both parents. The proposed procedure can provide detailed genetic information useful for marker assisted breeding in an outcrossing species such as forest trees.     

Detection of two QTL on chicken chromosome 14 for keyhole lymphet heamocyanin

A keyhole lymphet heamocyanin is an antigen which triggers Th1 type of immune response. A QTL for a primary immune response towards keyhole lymphet heamocyanin has been detected on chicken chromosome 14 in three populations. The results from the most recent population were inconsistent and varied depending on the applied QTL detection model. The major goal of the current study was the reanalysis of this data using a 2 QTL model. Additionally, in order to provide more accurate estimates of QTL effects and positions, epistasis between the QTL was considered as a potential important contributor to quantitative traits. Four statistical models were assumed: M1: A model assuming marginal additive effects of two QTL; M2: A model assuming marginal and epistatic additive effects of two QTL; M3: A model assuming marginal additive and dominance effects of two QTL; M4: A model assuming marginal additive and dominance effects of two QTL and all possible pairwise epistases. Two QTL with significant additive and dominance effects were detected on chicken chromosome 14 using model M3. One QTL was detected at 63 cM between MCW0123 and ROS0005, another at 76 cM between ROS0005 and MCW0225/NTN2Lsts1 (FDR = 0.0051). Modelling only additive effects resulted in a significantly worse fit. On the other hand, including epistatic effects did not improve fit significantly. The current study confirms previous reports of the QTL location on GGA14. A notable finding of this study is recognition of two closely related QTL for a keyhole lymphet heamocyanin response at the distal part of chicken chromosome 14.     

远交群体动态性状基因定位的似然分析Ⅰ.理论方法

受动物遗传育种中用来估计动态性状育种值的随机回归测定日模型思想的启发 ,将关于时间 (测定日期 )的Legendre多项式镶嵌在遗传模型的每个遗传效应中 ,以刻画QTL对动态性状变化过程的作用 ,从而建立起动态性状基因定位的数学模型。利用远交设计群体 ,阐述了动态性状基因定位的似然分析原理 ,推导了定位参数似然估计的EM法两步求解过程。结合动态性状遗传分析的特点和普通数量性状基因定位研究进展 ,还提出了有关动态性状基因定位进一步研究的设想     

杂交水稻中超显性效应的分析

本文提出亚群体分析法,在二个衍生于杂交稻推广组合的F2群体中,挑选均匀分布于连锁图谱的DNA标记作为固定因子,分别根据每个固定因子的基因型将F2群体分成三类亚群体：母本型（I型）、父本型（II型）和杂合型（III型）。在大量III型亚群体中,杂合度与产量和穗数呈显著正相关。在表现这种相关性的III型亚群体内分析产量QTL和穗数QTL,发现超显性作用是研究组合杂种优势的主要遗传基础。 Abstract:Two F2 populations were each derived from a commercial hybrid and used for identification of genetic factors contributing to heterosis of rice.DNA markers distributing evenly in the linkage maps were selected as fixing factors.The F2 populations were divided into three types of sub-populations based on the genotypes of each fixing factor:maternal type (Type I),paternal type (Type II) and heterozygous type (Type III).In a large number of Type III sub-populations,significant correlations were observed between heterozygosity and grain yield,and between heterozygosity and number of panicle.QTL analysis in Type III sub-populations showing such correlations indicated that over-dominance was the major genetic basis of heterosis in the two crosses.     

Mapping multiple QTL using linkage disequilibrium and linkage analysis information and multitrait data

A multi-locus QTL mapping method is presented, which combines linkage and linkage disequilibrium (LD) information and uses multitrait data. The method assumed a putative QTL at the midpoint of each marker bracket. Whether the putative QTL had an effect or not was sampled using Markov chain Monte Carlo (MCMC) methods. The method was tested in dairy cattle data on chromosome 14 where the DGAT1 gene was known to be segregating. The DGAT1 gene was mapped to a region of 0.04 cM, and the effects of the gene were accurately estimated. The fitting of multiple QTL gave a much sharper indication of the QTL position than a single QTL model using multitrait data, probably because the multi-locus QTL mapping reduced the carry over effect of the large DGAT1 gene to adjacent putative QTL positions. This suggests that the method could detect secondary QTL that would, in single point analyses, remain hidden under the broad peak of the dominant QTL. However, no indications for a second QTL affecting dairy traits were found on chromosome 14.     

A simulation study on the accuracy of position and effect estimates of linked QTL and their asymptotic standard deviations using multiple interval mapping in an F2 scheme

Approaches like multiple interval mapping using a multiple-QTL model for simultaneously mapping QTL can aid the identification of multiple QTL, improve the precision of estimating QTL positions and effects, and are able to identify patterns and individual elements of QTL epistasis. Because of the statistical problems in analytically deriving the standard errors and the distributional form of the estimates and because the use of resampling techniques is not feasible for several linked QTL, there is the need to perform large-scale simulation studies in order to evaluate the accuracy of multiple interval mapping for linked QTL and to assess confidence intervals based on the standard statistical theory. From our simulation study it can be concluded that in comparison with a monogenetic background a reliable and accurate estimation of QTL positions and QTL effects of multiple QTL in a linkage group requires much more information from the data. The reduction of the marker interval size from 10 cM to 5 cM led to a higher power in QTL detection and to a remarkable improvement of the QTL position as well as the QTL effect estimates. This is different from the findings for (single) interval mapping. The empirical standard deviations of the genetic effect estimates were generally large and they were the largest for the epistatic effects. These of the dominance effects were larger than those of the additive effects. The asymptotic standard deviation of the position estimates was not a good criterion for the accuracy of the position estimates and confidence intervals based on the standard statistical theory had a clearly smaller empirical coverage probability as compared to the nominal probability. Furthermore the asymptotic standard deviation of the additive, dominance and epistatic effects did not reflect the empirical standard deviations of the estimates very well, when the relative QTL variance was smaller/equal to 0.5. The implications of the above findings are discussed.     

Mapping quantitative trait loci in noninbred mosquito crosses

The identification of genes that affect quantitative traits has been of great interest to geneticists for many decades, and many statistical methods have been developed to map quantitative trait loci (QTL). Most QTL mapping studies in experimental organisms use purely inbred lines, where the two homologous chromosomes in each individual are identical. As a result, many existing QTL mapping methods developed for experimental organisms are applicable only to genetic crosses between inbred lines. However, it may be difficult to obtain inbred lines for certain organisms, e.g., mosquitoes. Although statistical methods for QTL mapping in outbred populations, e.g., humans, can be applied for such crosses, these methods may not fully take advantage of the uniqueness of these crosses. For example, we can generally assume that the two grandparental lines are homozygous at the QTL of interest, but such information is not be utilized through methods developed for outbred populations. In addition, mating types and phases can be relatively easy to establish through the analysis of adjacent markers due to the large number of offspring that can be collected, substantially simplifying the computational need. In this article, motivated by a mosquito intercross experiment involving two selected lines that are not genetically homozygous across the genome, we develop statistical methods for QTL mapping for genetic crosses involving noninbred lines. In our procedure, we first infer parental mating types and use likelihood-based methods to infer phases in each parent on the basis of genotypes of offspring and one parent. A hidden Markov model is then employed to estimate the number of high-risk alleles at marker positions and putative QTL positions between markers in each offspring, and QTL mapping is finally conducted through the inferred QTL configuration across all offspring in all crosses. The performance of the proposed methods is assessed through simulation studies, and the usefulness of this method is demonstrated through its application to a mosquito data set.