Genetic evaluation by best linear unbiased prediction using marker and trait information in a multibreed population.

Genetic evaluation by best linear unbiased prediction (BLUP) requires modeling genetic means, variances, and covariances. This paper presents theory to model means, variances, and covariances in a multibreed population, given marker and breed information, in the presence of gametic disequilibrium between the marker locus (ML) and linked quantitative trait locus (MQTL). Theory and algorithms are presented to construct the matrix of conditional covariances between relatives (Gv) for the MQTL effects in a multibreed population and to obtain the inverse of Gv efficiently. Theory presented here accounts for heterogeneity of variances among pure breeds and for segregation variances between pure breeds. A numerical example was used to illustrate how the theory and algorithms can be used for genetic evaluation by BLUP using marker and trait information in a multibreed population.     

The effect of using approximate gametic variance covariance matrices on marker assisted selection by BLUP

Under additive inheritance, the Henderson mixed model equations (HMME) provide an efficient approach to obtaining genetic evaluations by marker assisted best linear unbiased prediction (MABLUP) given pedigree relationships, trait and marker data. For large pedigrees with many missing markers, however, it is not feasible to calculate the exact gametic variance covariance matrix required to construct HMME. The objective of this study was to investigate the consequences of using approximate gametic variance covariance matrices on response to selection by MABLUP. Two methods were used to generate approximate variance covariance matrices. The first method (Method A) completely discards the marker information for individuals with an unknown linkage phase between two flanking markers. The second method (Method B) makes use of the marker information at only the most polymorphic marker locus for individuals with an unknown linkage phase. Data sets were simulated with and without missing marker data for flanking markers with 2, 4, 6, 8 or 12 alleles. Several missing marker data patterns were considered. The genetic variability explained by marked quantitative trait loci (MQTL) was modeled with one or two MQTL of equal effect. Response to selection by MABLUP using Method A or Method B were compared with that obtained by MABLUP using the exact genetic variance covariance matrix, which was estimated using 15 000 samples from the conditional distribution of genotypic values given the observed marker data. For the simulated conditions, the superiority of MABLUP over BLUP based only on pedigree relationships and trait data varied between 0.1% and 13.5% for Method A, between 1.7% and 23.8% for Method B, and between 7.6% and 28.9% for the exact method. The relative performance of the methods under investigation was not affected by the number of MQTL in the model.     

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.     

A model for marker-assisted selection among single crosses with multiple genetic markers

 Trait means of marker genotypes are often inconsistent across experiments, thereby hindering the use of regression techniques in marker-assisted selection. Best linear unbiased prediction based on trait and marker data (TM-BLUP) does not require prior information on the mean effects associated with specific marker genotypes and, consequently, may be useful in applied breeding programs. The objective of this paper is to present a flanking-marker, TM-BLUP model that is applicable to interpopulation single crosses that characterize maize (Zea mays L.) breeding programs. The performance of a single cross is modeled as the sum of testcross additive and dominance effects at unmarked quantitative trait loci (QTL) and at marked QTL (MQTL). The TM-BLUP model requires information on the recombination frequencies between flanking markers and the MQTL and on MQTL variances. A tabular method is presented for calculating the conditional probability that MQTL alleles in two inbreds are identical by descent given the observed marker genotypes (G ^k _obs) at the kth MQTL. Information on identity by descent of MQTL alleles can then be used to calculate the conditional covariance of MQTL effects between single crosses given G ^k _obs. The inverse of the covariance matrix for dominance effects at unmarked QTL and MQTL can be written directly from the inverse of the covariance matrices of the corresponding testcross additive effects. In practice, the computations required in TM-BLUP may be prohibitive. The computational requirements may be reduced with simplified TM-BLUP models wherein dominance effects at MQTL are excluded, only the single crosses that have been tested are included, or information is pooled across several MQTL. Received: 22 June 1997 / Accepted: 25 February 1998     

A recursive procedure to compute the gametic relationship matrix and its inverse for marked QTL clusters

Recently, a variety of mixed linear models have been proposed for marker-assisted prediction of the effects of quantitative trait loci (QTLs) in outbred populations of animals. One of them addresses the effects of a cluster of linked QTLs, or those of a particular chromosomal segment, marked by DNA marker(s) and requires that the inverse of the corresponding gametic relationship matrix whose elements are the conditional expected values of the identity-by-descent (IBD) proportions between gametes for individuals be evaluated. Here, for a model of this type, utilizing the property of the IBD set and using the information on the joint gametogenesis processes at the flanking marker loci, we present a recursive method to systematically calculate the elements of the gametic relationship matrix and its inverse. A numerical example is given to illustrate the proposed computing procedure.     

Prediction of additive genetic effects for the QTL-cluster on the basis of data on surrounding markers in outbred populations

This paper described a method for predicting additive effects of a cluster of tightly linked QTLs for outbred populations of animals in the situation where the QTLs are located on a chromosome segment surrounded by multiple linked DNA markers. We present a mixed model method for best linear unbiased prediction (conditional to the marker data) of the additive effects of the QTL-cluster and of the remaining QTLs unlinked to the marker linkage group. This method takes into consideration the identity-by-descent proportion (IBDP) for the particular chromosomal segment, in contrast to some other methods which use IBD probabilities at one specific location. In this method, fully informative data on different flanking markers is used to calculate the values of the expectations of the IBDPs (EIBDPs) between gametes for animals to be evaluated. Then the expected values are used as the elements of the gametic relationship matrix required in the best linear unbiased prediction. Giving a small numerical example, we illustrate how the present method can be used for the prediction of the QTL-cluster effects and for genetic evaluation of animals in outbred populations. A computational strategy is discussed on the basis of the calculation of the EIBDPs and the inverted gametic relationship matrix in complex pedigrees.     

Robust variance-components approach for assessing genetic linkage in pedigrees.

To assess evidence for genetic linkage from pedigrees, I developed a limited variance-components approach. In this method, variability among trait observations from individuals within pedigrees is expressed in terms of fixed effects from covariates and effects due to an unobservable trait-affecting major locus, random polygenic effects, and residual nongenetic variance. The effect attributable to a locus linked to a marker is a function of the additive and dominance components of variance of the locus, the recombination fraction, and the proportion of genes identical by descent at the marker locus for each pair of sibs. For unlinked loci, the polygenic variance component depends only on the relationship between the relative pair. Parameters can be estimated by either maximum-likelihood methods or quasi-likelihood methods. The forms of quasi-likelihood estimators are provided. Hypothesis tests derived from the maximum-likelihood approach are constructed by appeal to asymptotic theory. A simulation study showed that the size of likelihood-ratio tests was appropriate but that the monogenic component of variance was generally underestimated by the likelihood approach.     

Maximum likelihood techniques for the mapping and analysis of quantitative trait loci with the aid of genetic markers

A method is presented to estimate the biometric parameters of a quantitative trait locus linked to a genetic marker when both loci are segregating in the F-2 generation of a cross between two inbred lines. The method, which assumes underlying normal distributions, is a combination of maximum likelihood and moments methods and uses the statistics of the genetic marker genotype samples for the quantitative trait to estimate the recombination frequency between the two loci and the means and variances of the genotypes of the quantitative trait locus. With this method, the genetic parameters of a locus affecting plant height linked to an electrophoretic marker for esterase were accurately estimated from a sample of 1596 F-2 progeny of a cross between two species of Lycopersicon (tomato). Linkage distance between the two loci was 38 map units and the effect of the quantitative trait locus was 1.6 phenotypic standard deviation units. Accurate estimates of the genetic parameters and linkage distance for populations of 2000 individuals simulated with a segregating codominant locus with an effect of 1.63 standard deviations linked to a genetic marker with .2 recombination were also derived by this method. The method is not effective in distinguishing between complete and partial linkage in samples of only 500 individuals or for quantitative loci with effects less than a phenotypic standard deviation. The method is more effective for codominant than for dominant loci.     

Genetic evaluation methods for populations with dominance and inbreeding

Summary The effect of inbreeding on mean and genetic covariance matrix for a quantitative trait in a population with additive and dominance effects is shown. This genetic covariance matrix is a function of five relationship matrices and five genetic parameters describing the population. Elements of the relationship matrices are functions of Gillois (1964) identity coefficients for the four genes at a locus in two individuals. The equivalence of the path coefficient method (Jacquard 1966) and the tabular method (Smith and Mäki-Tanila 1990) to compute the covariance matrix of additive and dominance effects in a population with inbreeding is shown. The tabular method is modified to compute relationship matrices rather than the covariance matrix, which is trait dependent. Finally, approximate and exact Best Linear Unbiased Predictions (BLUP) of additive and dominance effects are compared using simulated data with inbreeding but no directional selection. The trait simulated was affected by 64 unlinked biallelic loci with equal effect and complete dominance. Simulated average inbreeding levels ranged from zero in generation one to 0.35 in generation five. The approximate method only accounted for the effect of inbreeding on mean and additive genetic covariance matrix, whereas the exact accounted for all of the changes in mean and genetic covariance matrix due to inbreeding. Approximate BLUP, which is computable for large populations where exact BLUP is not feasible, yielded unbiased predictions of additive and dominance effects in each generation with only slightly reduced accuracies relative to exact BLUP.     

A rapid method for computing the inverse of the gametic covariance matrix between relatives for a marked Quantitative Trait Locus

The inverse of the gametic covariance matrix between relatives, G^-1, for a marked quantitative trait locus (QTL) is required in best linear unbiased prediction (BLUP) of breeding values if marker data are available on a QTL. A rapid method for computing the inverse of a gametic relationship matrix for a marked QTL without building G itself is presented. The algorithm is particularly useful due to the approach taken in computing inbreeding coefficients by having to compute only few elements of G. Numerical techniques for determining, storing, and computing the required elements of G and the nonzero elements of the inverse are discussed. We show that the subset of G required for computing the inbreeding coefficients and hence the inverse is a tiny proportion of the whole matrix and can be easily stored in computer memory using sparse matrix storage techniques. We also introduce an algorithm to determine the maximum set of nonzero elements that can be found in G^-1 and a strategy to efficiently store and access them. Finally, we demonstrate that the inverse can be efficiently built using the present techniques for very large and inbred populations.     

Efficient computation of the inverse of gametic relationship matrix for a marked QTL

Best linear unbiased prediction of genetic merits for a marked quantitative trait locus (QTL) using mixed model methodology includes the inverse of conditional gametic relationship matrix (G^-1) for a marked QTL. When accounting for inbreeding, the conditional gametic relationships between two parents of individuals for a marked QTL are necessary to build G^-1 directly. Up to now, the tabular method and its adaptations have been used to compute these relationships. In the present paper, an indirect method was implemented at the gametic level to compute these few relationships. Simulation results showed that the indirect method can perform faster with significantly less storage requirements than adaptation of the tabular method. The efficiency of the indirect method was mainly due to the use of the sparseness of G^-1. The indirect method can also be applied to construct an approximate G^-1 for populations with incomplete marker data, providing approximate probabilities of descent for QTL alleles for individuals with incomplete marker data.     

Linkage Disequilibrium and Genetic Variances under Mutation-Selection Balance

I determine the contribution of linkage disequilibrium to genetic variances using results for two loci and for induced or marginal systems. The analysis allows epistasis and dominance, but assumes that mutation is weak relative to selection. The linkage disequilibrium component of genetic variance is shown to be unimportant for unlinked loci if the gametic mutation rate divided by the harmonic mean of the pairwise recombination rates is much less than one. For tightly linked loci, linkage disequilibrium is unimportant if the gametic mutation rate divided by the (induced) per locus selection is much less than one.     

Using molecular markers to map multiple quantitative trait loci: models for backcross,recombinant inbred,and doubled haploid progeny

Summary To maximize parameter estimation efficiency and statistical power and to estimate epistasis, the parameters of multiple quantitative trait loci (QTLs) must be simultaneously estimated. If multiple QTL affect a trait, then estimates of means of QTL genotypes from individual locus models are statistically biased. In this paper, I describe methods for estimating means of QTL genotypes and recombination frequencies between marker and quantitative trait loci using multilocus backcross, doubled haploid, recombinant inbred, and testcross progeny models. Expected values of marker genotype means were defined using no double or multiple crossover frequencies and flanking markers for linked and unlinked quantitative trait loci. The expected values for a particular model comprise a system of nonlinear equations that can be solved using an interative algorithm, e.g., the Gauss-Newton algorithm. The solutions are maximum likelihood estimates when the errors are normally distributed. A linear model for estimating the parameters of unlinked quantitative trait loci was found by transforming the nonlinear model. Recombination frequency estimators were defined using this linear model. Certain means of linked QTLs are less efficiently estimated than means of unlinked QTLs.     

The genetic variance for multiple linked quantitative trait loci conditional on marker information in a crossed population

In the prediction of genetic values and quantitative trait loci (QTLs) mapping via the mixed model method incorporating marker information in animal populations, it is important to model the genetic variance for individuals with an arbitrary pedigree structure. In this study, for a crossed population originated from different genetic groups such as breeds or outbred strains, the variance of additive genetic values for multiple linked QTLs that are contained in a chromosome segment, especially the segregation variance, is investigated assuming the use of marker data. The variance for a finite number of QTLs in one chromosomal segment is first examined for the crossed population with the general pedigree. Then, applying the concept of the expectation of identity-by-descent proportion, an approximation to the mean of the conditional probabilities for the linked QTLs over all loci is obtained, and using it an expression for the variance in the case of an infinite number of linked QTLs marked by flanking markers is derived. It appears that the approach presented can be useful in the segment mapping using, and in the genetic evaluation of, crosses with general pedigrees in the population of concern. The calculation of the segregation variance through the current approach is illustrated numerically, using a small data-set.     

Comparison of linkage-disequilibrium methods for localization of genes influencing quantitative traits in humans

Linkage disequilibrium has been used to help in the identification of genes predisposing to certain qualitative diseases. Although several linkage-disequilibrium tests have been developed for localization of genes influencing quantitative traits, these tests have not been thoroughly compared with one another. In this report we compare, under a variety of conditions, several different linkage-disequilibrium tests for identification of loci affecting quantitative traits. These tests use either single individuals or parent-child trios. When we compared tests with equal samples, we found that the truncated measured allele (TMA) test was the most powerful. The trait allele frequencies, the stringency of sample ascertainment, the number of marker alleles, and the linked genetic variance affected the power, but the presence of polygenes did not. When there were more than two trait alleles at a locus in the population, power to detect disequilibrium was greatly diminished. The presence of unlinked disequilibrium (D'*) increased the false-positive error rates of disequilibrium tests involving single individuals but did not affect the error rates of tests using family trios. The increase in error rates was affected by the stringency of selection, the trait allele frequency, and the linked genetic variance but not by polygenic factors. In an equilibrium population, the TMA test is most powerful, but, when adjusted for the presence of admixture, Allison test 3 becomes the most powerful whenever D'*>.15.     

Simulation of marker-assisted selection utilizing linkage disequilibrium: the effects of several additional factors

Use was made of our published model and methods to investigate the effects of several additional factors on marker-assisted selection (MAS) utilizing linkage disequilibrium. The additional factors were: size of the sample used to estimate the marker quantitative trait locus (MQTL) association effects, the method used to estimate the MQTL effects, use of the average of the top MQTL estimates in selection rather than individual estimates, size of the selection population, and the crossing of duplicate selection lines to generate further linkage disequilibrium and further selection response. The average map distance between the quantitative trait loci (QTLs) and their nearest marker was 0.15 Morgans. Use of estimates of MQTL effects derived by least squares yielded smaller selection responses than estimates derived by mixed-model methods. Selection responses were also reduced by using a smaller sample for estimating the associations because MQTL effects were less well estimated. This applied to selection on the MQTL effects themselves and to selection combining the MQTL with phenotypic information. Thus, poorly estimated MQTL effects added noise to the system and reduced selection response in combined selection. Using the average of the top MQTL estimates, rather than individual estimates, also reduced selection response. New linkage disequilibrium, generated by crossing two lines selected from the same population, did not lead to additional selection response in the cross line. These results show limitations to MAS using linkage disequilibrium until close linkages of markers and QTLs are available.     

Bayesian analysis of linkage between genetic markers and quantitative trait loci. I. Prior knowledge

Summary Prior information on gene effects at individual quantitative trait loci (QTL) and on recombination rates between marker loci and QTL is derived. The prior distribution of QTL gene effects is assumed to be exponential with major effects less likely than minor ones. The prior probability of linkage between a marker and another single locus is a function of the number and length of chromosomes, and of the map function relating recombination rate to genetic distance among loci. The prior probability of linkage between a marker locus and a quantitative trait depends additionally on the number of detectable QTL, which may be determined from total additive genetic variance and minimum detectable QTL effect. The use of this prior information should improve linkage tests and estimates of QTL effects.     

Effect of genetic heterogeneity and assortative mating on linkage analysis: a simulation study.

Linkage studies of complex genetic traits raise questions about the effects of genetic heterogeneity and assortative mating on linkage analysis. To further understand these problems, I have simulated and analyzed family data for a complex genetic disease in which disease phenotype is determined by two unlinked disease loci. Two models were studied, a two-locus threshold model and a two-locus heterogeneity model. Information was generated for a marker locus linked to one of the disease-defining loci. Random-mating and assortative-mating samples were generated. Linkage analysis was then carried out by use of standard methods, under the assumptions of a single-locus disease trait and a random-mating population. Results were compared with those from analysis of a single-locus homogeneous trait in samples with the same levels of assortative mating as those considered for the two-locus traits. The results show that (1) introduction of assortative mating does not, in itself, markedly affect the estimate of the recombination fraction; (2) the power of the analysis, reflected in the LOD scores, is somewhat lower with assortative rather than random mating. Loss of power is greater with increasing levels of assortative mating; and (3) for a heterogeneous genetic disease, regardless of mating type, heterogeneity analysis permits more accurate estimate of the recombination fraction but may be of limited use in distinguishing which families belong to each homogeneous subset. These simulations also confirmed earlier observations that linkage to a disease "locus" can be detected even if the disease is incorrectly defined as a single-locus (homogeneous) trait, although the estimated recombination fraction will be significantly greater than the true recombination fraction between the linked disease-defining locus and the marker locus.     

The effect of positive assortative mating at one locus on a second linked locus

Summary Considerations proceed from a model of positive assortative mating based on genotype at one locus, with an arbitrary number of alleles, assuming no selection, mutation, or migration, hypothetically infinite population size, and discrete non-overlapping generations. From these conditions, inferences are made about the genotypic structure at a linked locus, as well as about the corresponding 2-locus gametic structure.The following main results are presented: in the course of the generations, the genotypic structure at the second locus and the 2-locus gametic structure always tend to a limit responsive to the initial conditions concerning the joint genotypic structure at the two loci and the degree of assortativity and linkage. A complete, analytical representation of the limits is given. In particular, if assortative mating is only partial and at the same time linkage is not complete, a population is not able to maintain a permanent deviation of the gametic structure from linkage equilibrium, and thus the genotypic structure at the second locus tends to Hardy-Weinberg proportions. On the other hand, if initial linkage disequilibrium is combined with partial assortative mating and complete linkage (or with complete assortative mating and unlinked loci) the population maintains this disequilibrium and thus the genotypic structure at the second locus need not tend to Hardy-Weinberg proportions. It turns out that the conditions not only of complete linkage, but also of unlinked loci together with complete assortativity, imply no change in gametic structure from the initial structure.In order to demonstrate the influence of several parameters on the speed of convergence to and the magnitude of the respective limits, several graphs are included.     

Covariance between relatives in multibreed populations: additive model

Covariance between relatives in a multibreed population was derived for an additive model with multiple unlinked loci. An efficient algorithm to compute the inverse of the additive genetic covariance matrix is given. For an additive model, the variance for a crossbred individual is a function of the additive variances for the pure breeds, the covariance between parents, and segregation variances. Provided that the variance of a crossbred individual is computed as presented here, the covariance between crossbred relatives can be computed using formulae for purebred populations. For additive traits the inverse of the genotypic covariance matrix given here can be used both to obtain genetic evaluations by best linear unbiased prediction and to estimate genetic parameters by maximum likelihood in multibreed populations. For nonadditive traits, the procedure currently used to analyze multibreed data can be improved using the theory presented here to compute additive covariances together with a suitable approximation for nonadditive covariances.Supported in part by the Illinois Agricultural Experiment Station, Hatch Projects 35-0345 (RLF) and 35-0367 (MG)