Towards linkage analysis with markers in linkage disequilibrium by graphical modelling

We review recent developments of MCMC integration methods for computations on graphical models for two applications in statistical genetics: modelling allelic association and pedigree based linkage analysis. We discuss and illustrate estimation of graphical models from haploid and diploid genotypes, and the importance of MCMC updating schemes beyond what is strictly necessary for irreducibility. We then outline an approach combining these methods to compute linkage statistics when alleles at the marker loci are in linkage disequilibrium. Other extensions suitable for analysis of SNP genotype data in pedigrees are also discussed and programs that implement these methods, and which are available from the author's web site, are described. We conclude with a discussion of how this still experimental approach might be further developed.     

Identifying quantitative trait locus by genetic background interactions in association studies

Association studies are designed to identify main effects of alleles across a potentially wide range of genetic backgrounds. To control for spurious associations, effects of the genetic background itself are often incorporated into the linear model, either in the form of subpopulation effects in the case of structure or in the form of genetic relationship matrices in the case of complex pedigrees. In this context epistatic interactions between loci can be captured as an interaction effect between the associated locus and the genetic background. In this study I developed genetic and statistical models to tie the locus by genetic background interaction idea back to more standard concepts of epistasis when genetic background is modeled using an additive relationship matrix. I also simulated epistatic interactions in four-generation randomly mating pedigrees and evaluated the ability of the statistical models to identify when a biallelic associated locus was epistatic to other loci. Under additive-by-additive epistasis, when interaction effects of the associated locus were quite large (explaining 20% of the phenotypic variance), epistasis was detected in 79% of pedigrees containing 320 individuals. The epistatic model also predicted the genotypic value of progeny better than a standard additive model in 78% of simulations. When interaction effects were smaller (although still fairly large, explaining 5% of the phenotypic variance), epistasis was detected in only 9% of pedigrees containing 320 individuals and the epistatic and additive models were equally effective at predicting the genotypic values of progeny. Epistasis was detected with the same power whether the overall epistatic effect was the result of a single pairwise interaction or the sum of nine pairwise interactions, each generating one ninth of the epistatic variance. The power to detect epistasis was highest (94%) at low QTL minor allele frequency, fell to a minimum (60%) at minor allele frequency of about 0.2, and then plateaued at about 80% as alleles reached intermediate frequencies. The power to detect epistasis declined when the linkage disequilibrium between the DNA marker and the functional polymorphism was not complete.     

Fine mapping of complex trait genes combining pedigree and linkage disequilibrium information: a Bayesian unified framework

We present a Bayesian method that combines linkage and linkage disequilibrium (LDL) information for quantitative trait locus (QTL) mapping. This method uses jointly all marker information (haplotypes) and all available pedigree information; i.e., it is not restricted to any specific experimental design and it is not required that phases are known. Infinitesimal genetic effects or environmental noise ("fixed") effects can equally be fitted. A diallelic QTL is assumed and both additive and dominant effects can be estimated. We have implemented a combined Gibbs/Metropolis-Hastings sampling to obtain the marginal posterior distributions of the parameters of interest. We have also implemented a Bayesian variant of usual disequilibrium measures like D' and r(2) between QTL and markers. We illustrate the method with simulated data in "simple" (two-generation full-sib families) and "complex" (four-generation) pedigrees. We compared the estimates with and without using linkage disequilibrium information. In general, using LDL resulted in estimates of QTL position that were much better than linkage-only estimates when there was complete disequilibrium between the mutant QTL allele and the marker. This advantage, however, decreased when the association was only partial. In all cases, additive and dominant effects were estimated accurately either with or without disequilibrium information.     

A Monte Carlo method for Bayesian analysis of linkage between single markers and quantitative trait loci. I. Methodology

A Bayesian approach to the statistical mapping of Quantitative Trait Loci (QTLs) using single markers was implemented via Markov Chain Monte Carlo (MCMC) algorithms for parameter estimation and hypothesis testing. Parameter estimators were marginal posterior means computed using a Gibbs sampler with data augmentation. Variables sampled included the augmented data (marker-QTL genotypes, polygenic effects), an indicator variable for linkage, and the parameters (allele frequency, QTL substitution effect, recombination rate, polygenic and residual variances). Several MCMC algorithms were derived for computing Bayesian tests of linkage, which consisted of the marginal posterior probability of linkage and the marginal likelihood of the QTL variance associated with the marker.     

Linkage analysis of quantitative trait loci in multiple line crosses

Simple line crosses, for example, backcross and F₂, are commonly used in mapping quantitative trait loci (QTL). However, these simple crosses are rarely used alone in commercial plant breeding; rather, crosses involving multiple inbred lines or several simple crosses but connected by shared inbred lines may be common in plant breeding. Mapping QTL using crosses of multiple lines is more relevant to plant breeding. Unfortunately, current statistical methods and computer programs of QTL mapping are all designed for simple line crosses or multiple line crosses but under a regular mating system. It is not straightforward to extend the existing methods to handle multiple line crosses under irregular and complicated mating designs. The major hurdle comes from irregular inbreeding, multiple generations, and multiple alleles. In this study, we develop a Bayesian method implemented via the Markov chain Monte Carlo (MCMC) algorithm for mapping QTL using complicated multiple line crosses. With the MCMC algorithm, we are able to draw a complete path of the gene flow from founder alleles to their descendents via a recursive process. This has greatly simplified the problem caused by irregular mating and inbreeding in the mapping population. Adopting the reversible jump MCMC algorithm, we are able to simultaneously search for multiple QTL along the genome. We can even infer the posterior distribution of the number of QTL, one of the most important parameters in QTL study. Application of the new MCMC based QTL mapping procedure is demonstrated using two different mating designs. Design I involves two inbred lines and their derived F₁, F₂, and BC populations. Design II is a half-diallel cross involving three inbred lines. The two designs appear different, but can be handled with the same robust computer program.     

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.     

Simultaneous fine mapping of multiple closely linked quantitative trait Loci using combined linkage disequilibrium and linkage with a general pedigree

Within a small region (e.g., <10 cM), there can be multiple quantitative trait loci (QTL) underlying phenotypes of a trait. Simultaneous fine mapping of closely linked QTL needs an efficient tool to remove confounded shade effects among QTL within such a small region. We propose a variance component method using combined linkage disequilibrium (LD) and linkage information and a reversible jump Markov chain Monte Carlo (MCMC) sampling for model selection. QTL identity-by-descent (IBD) coefficients between individuals are estimated by a hybrid MCMC combining the random walk and the meiosis Gibbs sampler. These coefficients are used in a mixed linear model and an empirical Bayesian procedure combines residual maximum likelihood (REML) to estimate QTL effects and a reversible jump MCMC that samples the number of QTL and the posterior QTL intensities across the tested region. Note that two MCMC processes are used, i.e., an (internal) MCMC for IBD estimation and an (external) MCMC for model selection. In a simulation study, the use of the multiple-QTL model clearly removes the shade effects between three closely linked QTL located at 1.125, 3.875, and 7.875 cM across the region of 10 cM, using 40 markers at 0.25-cM intervals. It is shown that the use of combined LD and linkage information gives much more useful information compared to using linkage information alone for both single- and multiple-QTL analyses. When using a lower marker density (11 markers at 1-cM intervals), the signal of the second QTL can disappear. Extreme values of past effective size (resulting in extreme levels of LD) decrease the mapping accuracy.     

Bayesian analysis of linkage between genetic markers and quantitative trait loci. II. Combining prior knowledge with experimental evidence

Summary A Bayesian method was developed for identifying genetic markers linked to quantitative trait loci (QTL) by analyzing data from daughter or granddaughter designs and single markers or marker pairs. Traditional methods may yield unrealistic results because linkage tests depend on number of markers and QTL gene effects associated with selected markers are overestimated. The Bayesian or posterior probability of linkage combines information from a daughter or granddaughter design with the prior probability of linkage between a marker locus and a QTL. If the posterior probability exceeds a certain quantity, linkage is declared. Upon linkage acceptance, Bayesian estimates of marker-QTL recombination rate and QTL gene effects and frequencies are obtained. The Bayesian estimates of QTL gene effects account for different amounts of information by shrinking information from data toward the mean or mode of a prior exponential distribution of gene effects. Computation of the Bayesian analysis is feasible. Exact results are given for biallelic QTL, and extensions to multiallelic QTL are suggested.     

Derivation of single-locus relationship coefficients conditional on marker information

The coefficient of relationship is defined as the correlation between the additive genetic values of two individuals. This coefficient can be defined specifically for a single quantitative trait locus (QTL) and may deviate considerably from the overall expectation if it is taken conditional on information from linked marker loci. Conditional halfsib correlations are derived under a simple genetic model with a biallelic QTL linked to a biallelic marker locus. The conditional relationship coefficients are shown to depend on the recombination rate between the marker and the QTL and the population frequency of the marker alleles, but not on parameters of the QTL, i.e. number and frequency of QTL alleles, degree of dominance etc., nor on the (usually unknown) QTL genotype of the sire. Extensions to less simplified cases (multiple alleles at the marker locus and the QTL, two marker loci flanking the QTL) are given. For arbitrary pedigrees, conditional relationship coefficients can also be derived from the conditional gametic covariance matrix suggested by Fernando and Grossman (1989). The connection of these two approaches is discussed. The conditional relationship coefficient can be used for marker-assisted genetic evaluation as well as for the detection of QTL and the estimation of their effects.     

Identification of Quantitative Trait Loci Influencing Wood Specific Gravity in an Outbred Pedigree of Loblolly Pine

We report the identification of quantitative trait loci (QTL) influencing wood specific gravity (WSG) in an outbred pedigree of loblolly pine (Pinus taeda L.). QTL mapping in an outcrossing species is complicated by the presence of multiple alleles (>2) at QTL and marker loci. Multiple alleles at QTL allow the examination of interaction among alleles at QTL (deviation from additive gene action). Restriction fragment length polymorphism (RFLP) marker genotypes and wood specific gravity phenotypes were determined for 177 progeny. Two RFLP linkage maps were constructed, representing maternal and paternal parent gamete segregations as inferred from diploid progeny RFLP genotypes. RFLP loci segregating for multiple alleles were vital for aligning the two maps. Each RFLP locus was assayed for cosegregation with WSG QTL using analysis of variance (ANOVA). Five regions of the genome contained one or more RFLP loci showing differences in mean WSG at or below the P = 0.05 level for progeny as grouped by RFLP genotype. One region contained a marker locus (S6a) whose QTL-associated effects were highly significant (P > 0.0002). Marker S6a segregated for multiple alleles, a prerequisite for determining the number of alleles segregating at the linked QTL and analyzing the interactions among QTL alleles. The QTL associated with marker S6a appeared to be segregating for multiple alleles which interacted with each other and with environments. No evidence for digenic epistasis was found among the five QTL.     

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.     

Mixed model approaches for the identification of QTLs within a maize hybrid breeding program

Two outlines for mixed model based approaches to quantitative trait locus (QTL) mapping in existing maize hybrid selection programs are presented: a restricted maximum likelihood (REML) and a Bayesian Markov Chain Monte Carlo (MCMC) approach. The methods use the in-silico-mapping procedure developed by Parisseaux and Bernardo (2004) as a starting point. The original single-point approach is extended to a multi-point approach that facilitates interval mapping procedures. For computational and conceptual reasons, we partition the full set of relationships from founders to parents of hybrids into two types of relations by defining so-called intermediate founders. QTL effects are defined in terms of those intermediate founders. Marker based identity by descent relationships between intermediate founders define structuring matrices for the QTL effects that change along the genome. The dimension of the vector of QTL effects is reduced by the fact that there are fewer intermediate founders than parents. Furthermore, additional reduction in the number of QTL effects follows from the identification of founder groups by various algorithms. As a result, we obtain a powerful mixed model based statistical framework to identify QTLs in genetic backgrounds relevant to the elite germplasm of a commercial breeding program. The identification of such QTLs will provide the foundation for effective marker assisted and genome wide selection strategies. Analyses of an example data set show that QTLs are primarily identified in different heterotic groups and point to complementation of additive QTL effects as an important factor in hybrid performance.     

Fine mapping of a quantitative trait locus for twinning rate using combined linkage and linkage disequilibrium mapping

A novel and robust method for the fine-scale mapping of genes affecting complex traits, which combines linkage and linkage-disequilibrium information, is proposed. Linkage information refers to recombinations within the marker-genotyped generations and linkage disequilibrium to historical recombinations before genotyping started. The identity-by-descent (IBD) probabilities at the quantitative trait locus (QTL) between first generation haplotypes were obtained from the similarity of the marker alleles surrounding the QTL, whereas IBD probabilities at the QTL between later generation haplotypes were obtained by using the markers to trace the inheritance of the QTL. The variance explained by the QTL is estimated by residual maximum likelihood using the correlation structure defined by the IBD probabilities. Unlinked background genes were accounted for by fitting a polygenic variance component. The method was used to fine map a QTL for twinning rate in cattle, previously mapped on chromosome 5 by linkage analysis. The data consisted of large half-sib families, but the method could also handle more complex pedigrees. The likelihood of the putative QTL was very small along most of the chromosome, except for a sharp likelihood peak in the ninth marker bracket, which positioned the QTL within a region <1 cM in the middle part of bovine chromosome 5. The method was expected to be robust against multiple genes affecting the trait, multiple mutations at the QTL, and relatively low marker density.     

Multiple QTL mapping in related plant populations via a pedigree-analysis approach

QTL mapping experiments in plant breeding may involve multiple populations or pedigrees that are related through their ancestors. These known relationships have often been ignored for the sake of statistical analysis, despite their potential increase in power of mapping. We describe here a Bayesian method for QTL mapping in complex plant populations and reported the results from its application to a (previously analysed) potato data set. This Bayesian method was originally developed for human genetics data, and we have proved that it is useful for complex plant populations as well, based on a sensitivity analysis that was performed here. The method accommodates robustness to complex structures in pedigree data, full flexibility in the estimation of the number of QTL across multiple chromosomes, thereby accounting for uncertainties in the transmission of QTL and marker alleles due to incomplete marker information, and the simultaneous inclusion of non-genetic factors affecting the quantitative trait.     

Multiple Alleles for Tuber Shape in Diploid Potato Detected by Qualitative and Quantitative Genetic Analysis Using Rflps

Tuber shape in potato is commonly regarded as displaying continuous variation, yet at the diploid level phenotypes can be discerned visually, having round or long tubers. Inheritance of qualitative tuber shape can be explained by a single locus Ro, round being dominant to long. With restriction fragment length polymorphisms (RFLPs) the Ro locus was mapped on chromosome 10. Tuber shape was also studied as a quantitative trait, using the length/width ratio as trait value. The estimated broad sense heritability was h(2) = 0.80. The morphologically mapped Ro locus explained 75% of the genetic variation, indicating the presence of a major quantitative trait locus (QTL) at the Ro locus and minor genetic factors. RFLP alleles linked with Ro alleles were used to divide the progeny into four genotypic classes: Ro( &) Ro( &) : Ro( &) ro : roRo( &) : roro = 1 : 1 : 1 : 1. The recessive ro allele is identical by descent in both parents. The significantly different effects (P = 0.0157) of the non-identical alleles Ro( &) and Ro( &) provided evidence for multiallelism at the Ro locus. Linkage mapping of the Ro locus was compared with QTL mapping. Only those markers which are polymorphic in both parents allow accurate QTL mapping when genetic factors segregate from both parents. This finding applies to QTL mapping in all outbreeders without homozygous inbred strains.     

Bayesian mapping of quantitative trait loci under the identity-by-descent-based variance component model

Variance component analysis of quantitative trait loci (QTL) is an important strategy of genetic mapping for complex traits in humans. The method is robust because it can handle an arbitrary number of alleles with arbitrary modes of gene actions. The variance component method is usually implemented using the proportion of alleles with identity-by-descent (IBD) shared by relatives. As a result, information about marker linkage phases in the parents is not required. The method has been studied extensively under either the maximum-likelihood framework or the sib-pair regression paradigm. However, virtually all investigations are limited to normally distributed traits under a single QTL model. In this study, we develop a Bayes method to map multiple QTL. We also extend the Bayesian mapping procedure to identify QTL responsible for the variation of complex binary diseases in humans under a threshold model. The method can also treat the number of QTL as a parameter and infer its posterior distribution. We use the reversible jump Markov chain Monte Carlo method to infer the posterior distributions of parameters of interest. The Bayesian mapping procedure ends with an estimation of the joint posterior distribution of the number of QTL and the locations and variances of the identified QTL. Utilities of the method are demonstrated using a simulated population consisting of multiple full-sib families.     

Mapping quantitative trait loci using the MCMC procedure in SAS

The MCMC procedure in SAS (called PROC MCMC) is particularly designed for Bayesian analysis using the Markov chain Monte Carlo (MCMC) algorithm. The program is sufficiently general to handle very complicated statistical models and arbitrary prior distributions. This study introduces the SAS/MCMC procedure and demonstrates the application of the program to quantitative trait locus (QTL) mapping. A real life QTL mapping experiment in wheat female fertility trait was used as an example for the demonstration. The fertility trait phenotypes were described under three different models: (1) the Poisson model, (2) the Bernoulli model and (3) the zero-truncated Poisson model. One QTL was identified on the second chromosome. This QTL appears to control the switch of seed-producing ability of female plants but does not affect the number of seeds produced once the switch is turned on.     

Simultaneous mining of linkage and linkage disequilibrium to fine map quantitative trait loci in outbred half-sib pedigrees: revisiting the location of a quantitative trait locus with major effect on milk production on bovine chromosome 14

A maximum-likelihood QTL mapping method that simultaneously exploits linkage and linkage disequilibrium and that is applicable in outbred half-sib pedigrees is described. The method is applied to fine map a QTL with major effect on milk fat content in a 3-cM marker interval on proximal BTA14. This proximal location is confirmed by applying a haplotype-based association method referred to as recombinant ancestral haplotype analysis. The origin of the discrepancy between the QTL position derived in this work and that of a previous analysis is examined and shown to be due to the existence of distinct marker haplotypes associated with QTL alleles having large substitution effects.     

The genetic architecture of selection response. Inferences from fine-scale mapping of bristle number quantitative trait loci in Drosophila melanogaster.

Quantitative trait loci (QTL) affecting responses and correlated responses to selection for abdominal and sternopleural bristle number have been mapped with high resolution to the X and third chromosomes. Advanced intercross recombinant isogenic chromosomes were constructed from high and low selection lines in an unselected inbred background, and QTL were detected using composite interval mapping and high density transposable element marker maps. We mapped a total of 26 bristle number QTL with large effects, which were in or immediately adjacent to intervals previously inferred to contain bristle number QTL on these chromosomes. The QTL contributing to response to selection for high bristle number were not the same as those contributing to response to selection for low bristle number, suggesting that distributions of allelic effects per locus may be asymmetrical. Correlated responses were more often attributable to loose linkage than pleiotropy or close linkage. Bristle number QTL mapping to the same locations have been inferred in studies with different parental strains. Of the 26 QTL, 20 mapped to locations consistent with candidate genes affecting peripheral nervous system development and/or bristle number. This facilitates determining the molecular basis of quantitative variation and allele frequencies by associating molecular variation at the candidate genes with phenotypic variation in bristle number in samples of alleles from nature.     

Comparison of methods used for recovering the line origin of alleles in a cross between outbred lines

Here, we introduce the idea of probabilities of line origins for alleles in general pedigrees as found in crosses between outbred lines. We also present software for calculating these probabilities. The proposed algorithm is based on the linear regression method of Haley, Knott and Elsen (1994) combined with the Markov chain Monte Carlo (MCMC) method for estimating quantitative trait locus coefficients used as regressors. We compared the relative precision of our method and the original method as proposed by Haley et al. (1994). The scenarios studied varied in the allelic distribution of marker alleles in parental lines and in the frequency of missing marker genotypes. We found that the MCMC method achieves a higher accuracy in all scenarios considered. The benefits of using MCMC approximation are substantial if the frequency of missing marker data is high or the number of marker alleles is low and the allelic frequency distribution is similar in both parental lines.