Detection and localization of a single binary trait locus in experimental populations

The advancements made in molecular technology coupled with statistical methodology have led to the successful detection and location of genomic regions (quantitative trait loci; QTL) associated with quantitative traits. Binary traits (e.g. susceptibility/resistance), while not quantitative in nature, are equally important for the purpose of detecting and locating significant associations with genomic regions. Existing interval regression methods used in binary trait analysis are adapted from quantitative trait analysis and the tests for regression coefficients are tests of effect, not detection. Additionally, estimates of recombination that fail to take into account varying penetrance perform poorly when penetrance is incomplete. In this work a complete probability model for binary trait data is developed allowing for unbiased estimation of both penetrance and recombination between a genetic marker locus and a binary trait locus for backcross and F2 experimental designs. The regression model is reparameterized allowing for tests of detection. Extensive simulations were conducted to assess the performance of estimation and testing in the proposed parameterization. The proposed parameterization was compared with interval regression via simulation. The results indicate that our parameterization shows equivalent estimation capabilities, requires less computational effort and works well with only a single marker.     

A correlation method for detecting and estimating linkage between a marker locus and a quantitative trait locus using inbred lines

The advent of molecular genetic markers has stimulated interest in detecting linkage between a marker locus and a quantitative trait locus (QTL) because the marker locus, even without direct effect on the quantitative trait, could be useful in increasing the response to selection. A correlation method for detecting and estimating linkage between a marker locus and a QTL is described using selfing and sib-mating populations. Computer simulations were performed to estimate the power of the method, the sample size (N) needed to detect linkage, and the recombination value (r). The power of this method was a function of the expected recombination value E(r), the standardized difference (d) between the QTL genotypic means, and N. The power was highest at complete linkage, decreased with an increase in E(r), and then increased at E(r)=0.5. A larger d and N led to a higher power. The sample size needed to detect linkage was dependent upon E(r) and d. The sample size had a minimum value at E(r)=0, increased with an increase in E(r) and a decrease in d. In general, the r was overestimated. With an increase in d, the r was closer to its expectation. Detection of linkage by the proposed method under incomplete linkage was more efficient than estimation of recombination values. The correlation method and the method of comparison of marker-genotype means have a similar power when there is linkage, but the former has a slightly higher power than the latter when there is no linkage.     

Linkage analysis in the presence of errors I: complex-valued recombination fractions and complex phenotypes

Linkage is a phenomenon that correlates the genotypes of loci, rather than the phenotypes of one locus to the genotypes of another. It is therefore necessary to convert the observed trait phenotypes into trait-locus genotypes, which can then be analyzed for coinheritance with marker-locus genotypes. However, if the mode of inheritance of the trait is not known accurately, this conversion can often result in errors in the inferred trait-locus genotypes, which, in turn, can lead to the misclassification of the recombination status of meioses. As a result, the recombination fraction can be overestimated in two-point analysis, and false exclusions of the true trait locus can occur in multipoint analysis. We propose a method that increases the robustness of multipoint analysis to errors in the mode of inheritance assumptions of the trait, by explicitly allowing for misclassification of trait-locus genotypes. To this end, the definition of the recombination fraction is extended to the complex plane, as Theta=straight theta+straightepsiloni; theta is the recombination fraction between actual ("real") genotypes of marker and trait loci, and straightepsilon is the probability of apparent but false ("imaginary") recombinations between the actual and inferred trait-locus genotypes. "Complex" multipoint LOD scores are proven to be stochastically equivalent to conventional two-point LOD scores. The greater robustness to modeling errors normally associated with two-point analysis can thus be extended to multiple two-point analysis and multipoint analysis. The use of complex-valued recombination fractions also allows the stochastic equivalence of "model-based" and "model-free" methods to be extended to multipoint analysis.     

Mapping quantitative trait loci using molecular marker linkage maps

Summary High-density restriction fragment length polymorphism (RFLP) and allozyme linkage maps have been developed in several plant species. These maps make it technically feasible to map quantitative trait loci (QTL) using methods based on flanking marker genetic models. In this paper, we describe flanking marker models for doubled haploid (DH), recombinant inbred (RI), backcross (BC), F₁ testcross (F₁TC), DH testcross (DHTC), recombinant inbred testcross (RITC), F₂, and F₃ progeny. These models are functions of the means of quantitative trait locus genotypes and recombination frequencies between marker and quantitative trait loci. In addition to the genetic models, we describe maximum likelihood methods for estimating these parameters using linear, nonlinear, and univariate or multivariate normal distribution mixture models. We defined recombination frequency estimators for backcross and F₂ progeny group genetic models using the parameters of linear models. In addition, we found a genetically unbiased estimator of the QTL heterozygote mean using a linear function of marker means. In nonlinear models, recombination frequencies are estimated less efficiently than the means of quantitative trait locus genotypes. Recombination frequency estimation efficiency decreases as the distance between markers decreases, because the number of progeny in recombinant marker classes decreases. Mean estimation efficiency is nearly equal for these methods.     

Estimation of segregation and linkage parameters in simulated data. II. Simultaneous estimation with one linked marker.

This study examined the method of simultaneous estimation of recombination frequency and parameters for a qualitative trait locus and compared the results with those of standard methods of linkage analysis. With both approaches we were able to detect linkage of an incompletely penetrant qualitative trait to highly polymorphic markers with recombination frequencies in the range of .00-.05. Our results suggest that detecting linkage at larger recombination frequencies may require larger data sets or large high-density families. When applied to all families without regard to informativeness of the family structure for linkage, analyses of simulated data could detect no advantage of simultaneous estimation over more traditional and much less time-consuming methods, either in detecting linkage, estimating frequency, refining estimates of parameters for the qualitative trait locus, or avoiding false evidence for linkage. However, the method of sampling affected results.     

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

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

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.     

Misclassification in Logistic Regression with Discrete Covariates

We study the effect of misclassification of a binary covariate on the parameters of a logistic regression model. In particular we consider 2 × 2 × 2 tables. We assume that a binary covariate is subject to misclassification that may depend on the observed outcome. This type of misclassification is known as (outcome dependent) differential misclassification. We examine the resulting asymptotic bias on the parameters of the model and derive formulas for the biases and their approximations as a function of the odds and misclassification probabilities. Conditions for unbiased estimation are also discussed. The implications are illustrated numerically using a case control study. For completeness we briefly examine the effect of covariate dependent misclassification of exposures and of outcomes.     

Linkage between quantitative trait loci and marker loci: resolution power of three statistical approaches in single marker analysis

This paper presents a comparison of three methods of parameter estimation in analysis of linkage between a quantitative trait locus (QTL) and a marker locus: maximum likelihood, mean square for trait cumulative distribution function, and method of moments, employing simulated backcross data. The sensitivity of estimates to violation of assumptions of normality and equal variances were also studied. Some measures of discrepancy between the trait distributions in the QTL groups are considered to evaluate the potential dependence of the resolution capacity of the QTL substitution effect with respect to trait mean value and variance.     

Use of Multiple Genetic Markers in Prediction of Breeding Values

Genotypes at a marker locus give information on transmission of genes from parents to offspring and that information can be used in predicting the individuals' additive genetic value at a linked quantitative trait locus (MQTL). In this paper a recursive method is presented to build the gametic relationship matrix for an autosomal MQTL which requires knowledge on recombination rate between the marker locus and the MQTL linked to it. A method is also presented to obtain the inverse of the gametic relationship matrix. This information can be used in a mixed linear model for simultaneous evaluation of fixed effects, gametic effects at the MQTL and additive genetic effects due to quantitative trait loci unlinked to the marker locus (polygenes). An equivalent model can be written at the animal level using the numerator relationship matrix for the MQTL and a method for obtaining the inverse of this matrix is presented. Information on several unlinked marker loci, each of them linked to a different locus affecting the trait of interest, can be used by including an effect for each MQTL. The number of equations per animal in this case is 2m + 1 where m is the number of MQTL. A method is presented to reduce the number of equations per animal to one by combining information on all MQTL and polygenes into one numerator relationship matrix. It is illustrated how the method can accommodate individuals with partial or no marker information. Numerical examples are given to illustrate the methods presented. Opportunities to use the presented model in constructing genetic maps are discussed.     

Genetic Control of the Leaf Angle and Leaf Orientation Value as Revealed by Ultra-High Density Maps in Three Connected Maize Populations

Plant architecture is a key factor for high productivity maize because ideal plant architecture with an erect leaf angle and optimum leaf orientation value allow for more efficient light capture during photosynthesis and better wind circulation under dense planting conditions. To extend our understanding of the genetic mechanisms involved in leaf-related traits, three connected recombination inbred line (RIL) populations including 538 RILs were genotyped by genotyping-by-sequencing (GBS) method and phenotyped for the leaf angle and related traits in six environments. We conducted single population quantitative trait locus (QTL) mapping and joint linkage analysis based on high-density recombination bin maps constructed from GBS genotype data. A total of 45 QTLs with phenotypic effects ranging from 1.2% to 29.2% were detected for four leaf architecture traits by using joint linkage mapping across the three populations. All the QTLs identified for each trait could explain approximately 60% of the phenotypic variance. Four QTLs were located on small genomic regions where candidate genes were found. Genomic predictions from a genomic best linear unbiased prediction (GBLUP) model explained 45±9% to 68±8% of the variation in the remaining RILs for the four traits. These results extend our understanding of the genetics of leaf traits and can be used in genomic prediction to accelerate plant architecture improvement.     

Multiple regression for molecular-marker,quantitative trait data from large F2 populations

Molecular marker-quantitative trait associations are important for breeders to recognize and understand to allow application in selection. This work was done to provide simple, intuitive explanations of trait-marker regression for large samples from an F₂ and to examine the properties of the regression estimators. Beginning with a(- 1,0,1) coding of marker classes and expected frequencies in the F₂, expected values, variances, and covariances of marker variables were calculated. Simple linear regression and regression of trait values on two markers were computed. The sum of coefficient estimates for the flanking-marker regression is asymptotically unbiased for an included additive effect with complete interference, and is only slightly biased with no interference and moderately close (15 cM) marker spacing. The variance of the sum of regression coefficients is much more stable for small recombination distances than variances of individual coefficients. Multiple regression of trait variables on coded marker variables can be interpreted as the product of the inverse of the marker correlation matrix R and the vector a of simple linear regression estimators for each marker. For no interference, elements of the correlation matrix R can be written as products of correlations between adjacent markers. The inverse of R is displayed and used to illustrate the solution vector. Only markers immediately flanking trait loci are expected to have non-zero values and, with at least two marker loci between each trait locus, the solution vector is expected to be the sum of solutions for each trait locus. Results of this work should allow breeders to test for intervals in which trait loci are located and to better interpret results of the trait-marker regression.     

人类混血群体连锁不平衡的遗传学模型与基因定位

人类混血群体可以说是混合群体的一种特例．在无选择、无突变、无限随机交配群体的假定前提下,研究了亲本群体的基因频率对混血群体及其衍生后代群体连锁不平衡结构的影响,导出了各群体连锁不平衡值的表达式,建立了一个估计基因间重组率的简便方法;同时, 采用估算分子标记与QTL之间连锁不平衡系数的统计分析方法,分析了人类混血群体及其衍生后代群体QTL检测与估计的关系,建立了该关系的系列理论公式．研究结果表明,本方法不仅适用于人类疾病(包括复杂遗传疾病)基因定位,而且适合于人类正常基因的定位,同时也适用于人类普通多基因性状的QTL分析．     

Impact of and Correction for Outcome Misclassification in Cumulative Incidence Estimation

Cohort studies and clinical trials may involve multiple events. When occurrence of one of these events prevents the observance of another, the situation is called “competing risks”. A useful measure in such studies is the cumulative incidence of an event, which is useful in evaluating interventions or assessing disease prognosis. When outcomes in such studies are subject to misclassification, the resulting cumulative incidence estimates may be biased. In this work, we study the mechanism of bias in cumulative incidence estimation due to outcome misclassification. We show that even moderate levels of misclassification can lead to seriously biased estimates in a frequently unpredictable manner. We propose an easy to use estimator for correcting this bias that is uniformly consistent. Extensive simulations suggest that this method leads to unbiased estimates in practical settings. The proposed method is useful, both in settings where misclassification probabilities are known by historical data or can be estimated by other means, and for performing sensitivity analyses when the misclassification probabilities are not precisely known.     

An improved procedure of mapping a quantitative trait locus via the EM algorithm using posterior probabilities

Mapping a locus controlling a quantitative genetic trait (e.g. blood pressure) to a specific genomic region is of considerable contemporary interest. Data on the quantitative trait under consideration and several codominant genetic markers with known genomic locations are collected from members of families and statistically analysed to estimate the recombination fraction, θ, between the putative quantitative trait locus and a genetic marker. One of the major complications in estimating θ for a quantitative trait in humans is the lack of haplotype information on members of families. We have devised a computationally simple two-stage method of estimation of θ in the absence of haplotypic information using the expectation-maximization (EM) algorithm. In the first stage, parameters of the quantitative trait locus (QTL) are estimated on the basis of data of a sample of unrelated individuals and a Bayes’s rule is used to classify each parent into a QTL genotypic class. In the second stage, we have proposed an EM algorithm for obtaining the maximum-likelihood estimate of θ based on data of informative families (which are identified upon inferring parental QTL genotypes performed in the first stage). The purpose of this paper is to investigate whether, instead of using genotypically ‘classified’ data of parents, the use of posterior probabilities of QT genotypes of parents at the second stage yields better estimators. We show, using simulated data, that the proposed procedure using posterior probabilities is statistically more efficient than our earlier classification procedure, although it is computationally heavier.     

Robust multipoint identical-by-descent mapping for affected relative pairs

The genetic mapping of complex traits has been challenging and has required new statistical methods that are robust to misspecified models. Liang et al. proposed a robust multipoint method that can be used to simultaneously estimate, on the basis of sib-pair linkage data, both the position of a trait locus on a chromosome and its effect on disease status. The advantage of their method is that it does not require specification of an underlying genetic model, so estimation of the position of a trait locus on a specified chromosome and of its standard error is robust to a wide variety of genetic mechanisms. If multiple loci influence the trait, the method models the marginal effect of a locus on a specified chromosome. The main critical assumption is that there is only one trait locus on the chromosome of interest. We extend this method to different types of affected relative pairs (ARPs) by two approaches. One approach is to estimate the position of a trait locus yet allow unconstrained trait-locus effects across different types of ARPs. This robust approach allows for differences in sharing alleles identical-by-descent across different types of ARPs. Some examples for which an unconstrained model would apply are differences due to secular changes in diagnostic methods that can change the frequency of phenocopies among different types of relative pairs, environmental factors that modify the genetic effect, epistasis, and variation in marker-information content. However, this unconstrained model requires a parameter for each type of relative pair. To reduce the number of parameters, we propose a second approach that models the marginal effect of a susceptibility locus. This constrained model is robust for a trait caused by either a single locus or by multiple loci without epistasis. To evaluate the adequacy of the constrained model, we developed a robust score statistic. These methods are applied to a prostate cancer-linkage study, which emphasizes their potential advantages and limitations.     

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.     

Linkage group correction using epistatic distorted markers in F2 and backcross populations

Epistasis has been frequently observed in all types of mapping populations. However, relatively little is known about the effect of epistatic distorted markers on linkage group construction. In this study, a new approach was proposed to correct the recombination fraction between epistatic distorted markers in backcross and F₂ populations under the framework of fitness and liability models. The information for three or four markers flanking with an epistatic segregation distortion locus was used to estimate the recombination fraction by the maximum likelihood method, implemented via an expectation–maximisation algorithm. A set of Monte Carlo simulation experiments along with a real data analysis in rice was performed to validate the new method. The results showed that the estimates from the new method are unbiased. In addition, five statistical properties for the new method in a backcross were summarised and confirmed by theoretical, simulated and real data analyses.     

Estimation of the contribution of quantitative trait loci (QTL) to the variance of a quantitative trait by means of genetic markers

The estimation of the contribution of an individual quantitative trait locus (QTL) to the variance of a quantitative trait is considered in the framework of an analysis of variance (ANOVA). ANOVA mean squares expectations which are appropriate to the specific case of QTL mapping experiments are derived. These expectations allow the specificities associated with the limited number of genotypes at a given locus to be taken into account. Discrepancies with classical expectations are particularly important for two-class experiments (backcross, recombinant inbred lines, doubled haploid populations) and F₂ populations. The result allows us firstly to reconsider the power of experiments (i.e. the probability of detecting a QTL with a given contribution to the variance of the trait). It illustrates that the use of classical formulae for mean squares expectations leads to a strong underestimation of the power of the experiments. Secondly, from the observed mean squares it is possible to estimate directly the variance associated with a locus and the fraction of the total variance associated to this locus (r _l ² ). When compared to other methods, the values estimated using this method are unbiased. Considering unbiased estimators increases in importance when (1) the experimental size is limited; (2) the number of genotypes at the locus of interest is large; and (3) the fraction of the variation associated with this locus is small. Finally, specific mean squares expectations allows us to propose a simple analytical method by which to estimate the confidence interval of r _l ² . This point is particularly important since results indicate that 95% confidence intervals for r _l ² can be rather wide:2–23% for a 10% estimate and 8–34% for a 20% estimate if 100 individuals are considered.     

A Method of Screening for Genes of Major Effect

This paper describes a method for screening animal populations on an index of calculated probabilities of genotype status at an unknown single locus. Animals selected by such a method might then be candidates in test matings and genetic marker analyses for major gene detection. The method relies on phenotypic measures for a continuous trait plus identification of sire and dam. Some missing phenotypes and missing pedigree information are permitted. The method is an iterative two-step procedure, the first step estimates genotype probabilities and the second step estimates genotypic effects by regressing phenotypes on genotype probabilities, modeled as true genotype status plus error. Prior knowledge or choice of major locus-free heritability for the trait of interest is required, plus initial starting estimates of the effect on phenotype of carrying one and two copies of the unknown gene. Gene frequency can be estimated by this method, but it is demonstrated that the consequences of using an incorrect fixed prior for gene frequency are not particularly adverse where true frequency of the allele with major effect is low. Simulations involving deterministic sampling from the normal distribution lead to convergence for estimates of genotype effects at the true values, for a reasonable range of starting values, illustrating that estimation of major gene effects has a rational basis. In the absence of polygenic effects, stochastic simulations of 600 animals in five generations resulted in estimates of genotypic effects close to the true values. However, stochastic simulations involving generation and fitting of both major genotype and animal polygenic effects showed upward bias in estimates of major genotype effects. This can be partially overcome by not using information from relatives when calculating genotype probabilities-a result which suggests a route to a modified method which is unbiased and yet does use this information.