Joint linkage and linkage disequilibrium mapping of quantitative trait loci in natural populations

Linkage analysis and allelic association (also referred to as linkage disequilibrium) studies are two major approaches for mapping genes that control simple or complex traits in plants, animals, and humans. But these two approaches have limited utility when used alone, because they use only part of the information that is available for a mapping population. More recently, a new mapping strategy has been designed to integrate the advantages of linkage analysis and linkage disequilibrium analysis for genome mapping in outcrossing populations. The new strategy makes use of a random sample from a panmictic population and the open-pollinated progeny of the sample. In this article, we extend the new strategy to map quantitative trait loci (QTL), using molecular markers within the EM-implemented maximum-likelihood framework. The most significant advantage of this extension is that both linkage and linkage disequilibrium between a marker and QTL can be estimated simultaneously, thus increasing the efficiency and effectiveness of genome mapping for recalcitrant outcrossing species. Simulation studies are performed to test the statistical properties of the MLEs of genetic and genomic parameters including QTL allele frequency, QTL effects, QTL position, and the linkage disequilibrium of the QTL and a marker. The potential utility of our mapping strategy is discussed.     

A comparison between methods for linkage disequilibrium fine mapping of quantitative trait loci

We present a maximum likelihood method for mapping quantitative trait loci that uses linkage disequilibrium information from single and multiple markers. We made paired comparisons between analyses using a single marker, two markers and six markers. We also compared the method to single marker regression analysis under several scenarios using simulated data. In general, our method outperformed regression (smaller mean square error and confidence intervals of location estimate) for quantitative trait loci with dominance effects. In addition, the method provides estimates of the frequency and additive and dominance effects of the quantitative trait locus.     

Power and precision of alternate methods for linkage disequilibrium mapping of quantitative trait loci

Linkage disequilibrium (LD) analysis in outbred populations uses historical recombinations to detect and fine map quantitative trait loci (QTL). Our objective was to evaluate the effect of various factors on power and precision of QTL detection and to compare LD mapping methods on the basis of regression and identity by descent (IBD) in populations of limited effective population size (N(e)). An 11-cM region with 6-38 segregating single-nucleotide polymorphisms (SNPs) and a central QTL was simulated. After 100 generations of random mating with N(e) of 50, 100, or 200, SNP genotypes and phenotypes were generated on 200, 500, or 1000 individuals with the QTL explaining 2 or 5% of phenotypic variance. To detect and map the QTL, phenotypes were regressed on genotypes or (assumed known) haplotypes, in comparison with the IBD method. Power and precision to detect QTL increased with sample size, marker density, and QTL effect. Power decreased with N(e), but precision was affected little by N(e). Single-marker regression had similar or greater power and precision than other regression models, and was comparable to the IBD method. Thus, for rapid initial screening of samples of adequate size in populations in which drift is the primary force that has created LD, QTL can be detected and mapped by regression on SNP genotypes without recovering haplotypes.     

A haplotype-based algorithm for multilocus linkage disequilibrium mapping of quantitative trait loci with epistasis

For tightly linked loci, cosegregation may lead to nonrandom associations between alleles in a population. Because of its evolutionary relationship with linkage, this phenomenon is called linkage disequilibrium. Today, linkage disequilibrium-based mapping has become a major focus of recent genome research into mapping complex traits. In this article, we present a new statistical method for mapping quantitative trait loci (QTL) of additive, dominant, and epistatic effects in equilibrium natural populations. Our method is based on haplotype analysis of multilocus linkage disequilibrium and exhibits two significant advantages over current disequilibrium mapping methods. First, we have derived closed-form solutions for estimating the marker-QTL haplotype frequencies within the maximum-likelihood framework implemented by the EM algorithm. The allele frequencies of putative QTL and their linkage disequilibria with the markers are estimated by solving a system of regular equations. This procedure has significantly improved the computational efficiency and the precision of parameter estimation. Second, our method can detect marker-QTL disequilibria of different orders and QTL epistatic interactions of various kinds on the basis of a multilocus analysis. This can not only enhance the precision of parameter estimation, but also make it possible to perform whole-genome association studies. We carried out extensive simulation studies to examine the robustness and statistical performance of our method. The application of the new method was validated using a case study from humans, in which we successfully detected significant QTL affecting human body heights. Finally, we discuss the implications of our method for genome projects and its extension to a broader circumstance. The computer program for the method proposed in this article is available at the webpage http://www.ifasstat.ufl.edu/genome/~LD.     

Linkage disequilibrium mapping of quantitative trait loci under truncation selection

As a dense map of single nucleotide polymorphism (SNP) markers are available, population-based linkage disequilibrium (LD) mapping or association study is becoming one of the major tools for identifying quantitative trait loci (QTL) and for fine gene mapping. However, in many cases, LD between the marker and trait locus is not very strong. Approaches that maximize the potential of detecting LD will be essential for the success of LD mapping of QTL. In this paper, we propose two strategies for increasing the probability of detecting LD: (1) phenotypic selection and (2) haplotype LD mapping. To provide the foundations for LD mapping of QTL under selection, we develop analytic tools for assessing the impact of phenotypic selection on allele and haplotype frequencies, and LD under three trait models: single trait locus, two unlinked trait loci, and two linked trait loci with or without epistasis. In addition to a traditional chi(2) test, which compares the difference in allele or haplotype frequencies in the selected sample and population sample, we present multiple regression methods for LD mapping of QTL, and investigate which methods are effective in employing phenotypic selection for QTL mapping. We also develop a statistical framework for investigating and comparing the power of the single marker and multilocus haplotype test for LD mapping of QTL. Finally, the proposed methods are applied to mapping QTL influencing variation in systolic blood pressure in an isolated Chinese population.     

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

Causal mutations and their intra- and inter-locus interactions play a critical role in complex trait variation. It is often not easy to detect epistatic quantitative trait loci (QTL) due to complicated population structure requirements for detecting epistatic effects in linkage analysis studies and due to main effects often being hidden by interaction effects. Mapping their positions is even harder when they are closely linked. The data structure requirement may be overcome when information on linkage disequilibrium is used. We present an approach using a mixed linear model nested in an empirical Bayesian approach, which simultaneously takes into account additive, dominance and epistatic effects due to multiple QTL. The covariance structure used in the mixed linear model is based on combined linkage disequilibrium and linkage information. In a simulation study where there are complex epistatic interactions between QTL, it is possible to simultaneously map interacting QTL into a small region using the proposed approach. The estimated variance components are accurate and less biased with the proposed approach compared with traditional models.     

High-resolution mapping of quantitative trait loci by selective recombinant genotyping

Selective recombinant genotyping (SRG) is a three-stage procedure for high-resolution mapping of a QTL that has previously been mapped to a known confidence interval (target C.I.). In stage 1, a large mapping population is accessed and phenotyped, and a proportion, P, of the high and low tails is selected. In stage 2, the selected individuals are genotyped for a pair of markers flanking the target C.I., and a group of R individuals carrying recombinant chromosomes in the target interval are identified. In stage 3, the recombinant individuals are genotyped for a set of M markers spanning the target C.I. Extensive simulations showed that: (1) Standard error of QTL location (SEQTL) decreased when QTL effect (d) or population size (N) increased, but was constant for given "power factor" (PF = d(2)N); (2) increasing the proportion selected in the tails beyond 0.25 had only a negligible effect on SEQTL; and (3) marker spacing in the target interval had a remarkably powerful effect on SEQTL, yielding a reduction of up to 10-fold in going from highest (24 cM) to lowest (0.29 cM) spacing at given population size and QTL effect. At the densest marker spacing, SEQTL of 1.0-0.06 cM were obtained at PF = 500-16,000. Two new genotyping procedures, the half-section algorithm and the golden section/half-section algorithm, allow the equivalent of complete haplotyping of the target C.I. in the recombinant individuals to be achieved with many fewer data points than would be required by complete individual genotyping.     

The role of pedigree information in combined linkage disequilibrium and linkage mapping of quantitative trait loci in a general complex pedigree

Combined linkage disequilibrium and linkage (LDL) mapping can exploit historical as well as recent and observed recombinations in a recorded pedigree. We investigated the role of pedigree information in LDL mapping and the performance of LDL mapping in general complex pedigrees. We compared using complete and incomplete genotypic data, spanning 5 or 10 generations of known pedigree, and we used bi- or multiallelic markers that were positioned at 1- or 5-cM intervals. Analyses carried out with or without pedigree information were compared. Results were compared with linkage mapping in some of the data sets. Linkage mapping or LDL mapping with sparse marker spacing ( approximately 5 cM) gave a poorer mapping resolution without considering pedigree information compared to that with considering pedigree information. The difference was bigger in a pedigree of more generations. However, LDL mapping with closely linked markers ( approximately 1 cM) gave a much higher mapping resolution regardless of using pedigree information. This study shows that when marker spacing is dense and there is considerable linkage disequilibrium generated from historical recombinations between flanking markers and QTL, the loss of power due to ignoring pedigree information is negligible and mapping resolution is very high.     

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.     

Comparing linkage disequilibrium-based methods for fine mapping quantitative trait loci

Recently, a method for fine mapping quantitative trait loci (QTL) using linkage disequilibrium was proposed to map QTL by modeling covariance between individuals, due to identical-by-descent (IBD) QTL alleles, on the basis of the similarity of their marker haplotypes under an assumed population history. In the work presented here, the advantage of using marker haplotype information for fine mapping QTL was studied by comparing the IBD-based method with 10 markers to regression on a single marker, a pair of markers, or a two-locus haplotype under alternative population histories. When 10 markers were genotyped, the IBD-based method estimated the position of the QTL more accurately than did single-marker regression in all populations. When 20 markers were genotyped for regression, as single-marker methods do not require knowledge of haplotypes, the mapping accuracy of regression in all populations was similar to or greater than that of the IBD-based method using 10 markers. Thus for populations similar to those simulated here, the IBD-based method is comparable to single-marker regression analysis for fine mapping QTL.     

The efficiency of designs for fine-mapping of quantitative trait loci using combined linkage disequilibrium and linkage

In a simulation study, different designs were compared for efficiency of fine-mapping of QTL. The variance component method for fine-mapping of QTL was used to estimate QTL position and variance components. The design of many families with small size gave a higher mapping resolution than a design with few families of large size. However, the difference is small in half sib designs. The proportion of replicates with the QTL positioned within 3 cM of the true position is 0.71 in the best design, and 0.68 in the worst design applied to 128 animals with a phenotypic record and a QTL explaining 25% of the phenotypic variance. The design of two half sib families each of size 64 was further investigated for a hypothetical population with effective size of 1000 simulated for 6000 generations with a marker density of 0.25 cM and with marker mutation rate 4 × 10^-4 per generation. In mapping using bi-allelic markers, 42~55% of replicated simulations could position QTL within 0.75 cM of the true position whereas this was higher for multi allelic markers (48~76%). The accuracy was lowest (48%) when mutation age was 100 generations and increased to 68% and 76% for mutation ages of 200 and 500 generations, respectively, after which it was about 70% for mutation ages of 1000 generations and older. When effective size was linearly decreasing in the last 50 generations, the accuracy was decreased (56 to 70%). We show that half sib designs that have often been used for linkage mapping can have sufficient information for fine-mapping of QTL. It is suggested that the same design with the same animals for linkage mapping should be used for fine-mapping so gene mapping can be cost effective in livestock populations.     

High-resolution association mapping of quantitative trait loci: a population-based approach

In this article, population-based regression models are proposed for high-resolution linkage disequilibrium mapping of quantitative trait loci (QTL). Two regression models, the "genotype effect model" and the "additive effect model," are proposed to model the association between the markers and the trait locus. The marker can be either diallelic or multiallelic. If only one marker is used, the method is similar to a classical setting by Nielsen and Weir, and the additive effect model is equivalent to the haplotype trend regression (HTR) method by Zaykin et al. If two/multiple marker data with phase ambiguity are used in the analysis, the proposed models can be used to analyze the data directly. By analytical formulas, we show that the genotype effect model can be used to model the additive and dominance effects simultaneously; the additive effect model takes care of the additive effect only. On the basis of the two models, F-test statistics are proposed to test association between the QTL and markers. By a simulation study, we show that the two models have reasonable type I error rates for a data set of moderate sample size. The noncentrality parameter approximations of F-test statistics are derived to make power calculation and comparison. By a simulation study, it is found that the noncentrality parameter approximations of F-test statistics work very well. Using the noncentrality parameter approximations, we compare the power of the two models with that of the HTR. In addition, a simulation study is performed to make a comparison on the basis of the haplotype frequencies of 10 SNPs of angiotensin-1 converting enzyme (ACE) genes.     

Multitrait fine mapping of quantitative trait loci using combined linkage disequilibria and linkage analysis

A novel multitrait fine-mapping method is presented. The method is implemented by a model that treats QTL effects as random variables. The covariance matrix of allelic effects is proportional to the IBD matrix, where each element is the probability that a pair of alleles is identical by descent, given marker information and QTL position. These probabilities are calculated on the basis of similarities of marker haplotypes of individuals of the first generation of genotyped individuals, using "gene dropping" (linkage disequilibrium) and transmission of markers from genotyped parents to genotyped offspring (linkage). A small simulation study based on a granddaughter design was carried out to illustrate that the method provides accurate estimates of QTL position. Results from the simulation also indicate that it is possible to distinguish between a model postulating one pleiotropic QTL affecting two traits vs. one postulating two closely linked loci, each affecting one of the traits.     

Linkage disequilibrium fine mapping of quantitative trait loci: A simulation study

Recently, the use of linkage disequilibrium (LD) to locate genes which affect quantitative traits (QTL) has received an increasing interest, but the plausibility of fine mapping using linkage disequilibrium techniques for QTL has not been well studied. The main objectives of this work were to (1) measure the extent and pattern of LD between a putative QTL and nearby markers in finite populations and (2) investigate the usefulness of LD in fine mapping QTL in simulated populations using a dense map of multiallelic or biallelic marker loci. The test of association between a marker and QTL and the power of the test were calculated based on single-marker regression analysis. The results show the presence of substantial linkage disequilibrium with closely linked marker loci after 100 to 200 generations of random mating. Although the power to test the association with a frequent QTL of large effect was satisfactory, the power was low for the QTL with a small effect and/or low frequency. More powerful, multi-locus methods may be required to map low frequent QTL with small genetic effects, as well as combining both linkage and linkage disequilibrium information. The results also showed that multiallelic markers are more useful than biallelic markers to detect linkage disequilibrium and association at an equal distance.     

Combined linkage and association mapping of quantitative trait loci by multiple markers

Using multiple diallelic markers, variance component models are proposed for high-resolution combined linkage and association mapping of quantitative trait loci (QTL) based on nuclear families. The objective is to build a model that may fully use marker information for fine association mapping of QTL in the presence of prior linkage. The measures of linkage disequilibrium and the genetic effects are incorporated in the mean coefficients and are decomposed into orthogonal additive and dominance effects. The linkage information is modeled in variance-covariance matrices. Hence, the proposed methods model both association and linkage in a unified model. On the basis of marker information, a multipoint interval mapping method is provided to estimate the proportion of allele sharing identical by descent (IBD) and the probability of sharing two alleles IBD at a putative QTL for a sib-pair. To test the association between the trait locus and the markers, both likelihood-ratio tests and F-tests can be constructed on the basis of the proposed models. In addition, analytical formulas of noncentrality parameter approximations of the F-test statistics are provided. Type I error rates of the proposed test statistics are calculated to show their robustness. After comparing with the association between-family and association within-family (AbAw) approach by Abecasis and Fulker et al., it is found that the method proposed in this article is more powerful and advantageous based on simulation study and power calculation. By power and sample size comparison, it is shown that models that use more markers may have higher power than models that use fewer markers. The multiple-marker analysis can be more advantageous and has higher power in fine mapping QTL. As an application, the Genetic Analysis Workshop 12 German asthma data are analyzed using the proposed methods.     

Fine mapping of quantitative trait loci using linkage disequilibria with closely linked marker loci

A multimarker linkage disequilibrium mapping method was developed for the fine mapping of quantitative trait loci (QTL) using a dense marker map. The method compares the expected covariances between haplotype effects given a postulated QTL position to the covariances that are found in the data. The expected covariances between the haplotype effects are proportional to the probability that the QTL position is identical by descent (IBD) given the marker haplotype information, which is calculated using the genedropping method. Simulation results showed that a QTL was correctly positioned within a region of 3, 1.5, or 0.75 cM in 70, 62, and 68%, respectively, of the replicates using markers spaced at intervals of 1, 0.5, and 0.25 cM, respectively. These results were rather insensitive to the number of generations since the QTL occurred and to the effective population size, except that 10 generations yielded rather poor estimates of the QTL position. The position estimates of this multimarker disequilibrium mapping method were more accurate than those from a single marker transmission disequilibrium test. A general approach for identifying QTL is suggested, where several stages of disequilibrium mapping are used with increasingly dense marker spacing.     

Semiparametric methods for mapping quantitative trait loci with censored data

Statistical methods for the detection of genes influencing quantitative traits with the aid of genetic markers are well developed for normally distributed, fully observed phenotypes. Many experiments are concerned with failure-time phenotypes, which have skewed distributions and which are usually subject to censoring because of random loss to follow-up, failures from competing causes, or limited duration of the experiment. In this article, we develop semiparametric statistical methods for mapping quantitative trait loci (QTLs) based on censored failure-time phenotypes. We formulate the effects of the QTL genotype on the failure time through the Cox (1972, Journal of the Royal Statistical Society, Series B 34, 187-220) proportional hazards model and derive efficient likelihood-based inference procedures. In addition, we show how to assess statistical significance when searching several regions or the entire genome for QTLs. Extensive simulation studies demonstrate that the proposed methods perform well in practical situations. Applications to two animal studies are provided.     

High-resolution mapping of quantitative trait loci affecting increased life span in Drosophila melanogaster

Limited life span and senescence are near-universal characteristics of eukaryotic organisms, controlled by many interacting quantitative trait loci (QTL) with individually small effects, whose expression is sensitive to the environment. Analyses of mutations in model organisms have shown that genes affecting stress resistance and metabolism affect life span across diverse taxa. However, there is considerable segregating variation for life span in nature, and relatively little is known about the genetic basis of this variation. Replicated lines of Drosophila that have evolved increased longevity as a correlated response to selection for postponed senescence are valuable resources for identifying QTL affecting naturally occurring variation in life span. Here, we used deficiency complementation mapping to identify at least 11 QTL on chromosome 3 that affect variation in life span between five old (O) lines selected for postponed senescence and their five base (B) population control lines. Most QTL were sex specific, and all but one affected multiple O lines. The latter observation is consistent with alleles at intermediate frequency in the base population contributing to the response to selection for postponed senescence. The QTL were mapped with high resolution and contained from 12 to 170 positional candidate genes.