DNA fingerprint bands applied to linkage analysis with quantitative trait loci in chickens

Summary
An efficient approach to detect association between quantitative traits and bands of DNA fingerprint patterns uses intra-family tail analysis, which compares fingerprints of DNA mixes from individuals at the two tails of a phenotypic distribution. In analysis of 67 paternal half-sibs of a meat-type chicken family, of 57 sire bands generated by two probes, one sire-specific band (S6–6) was associated with abdominal fat deposition. The band effect was estimated by a linear model analysis to be 0–88 standard deviations, or about 30% of the family mean. The association between band S6–6 and abdominal fat was further examined by testing progeny of paternal half-sibs of the chickens which were used in the tail analysis, establishing genetic linkage between the DNA marker and a genetic locus affecting abdominal fat deposition.     

An efficient resampling method for assessing genome-wide statistical significance in mapping quantitative trait Loci

Assessing genome-wide statistical significance is an important and difficult problem in multipoint linkage analysis. Due to multiple tests on the same genome, the usual pointwise significance level based on the chi-square approximation is inappropriate. Permutation is widely used to determine genome-wide significance. Theoretical approximations are available for simple experimental crosses. In this article, we propose a resampling procedure to assess the significance of genome-wide QTL mapping for experimental crosses. The proposed method is computationally much less intensive than the permutation procedure (in the order of 10(2) or higher) and is applicable to complex breeding designs and sophisticated genetic models that cannot be handled by the permutation and theoretical methods. The usefulness of the proposed method is demonstrated through simulation studies and an application to a Drosophila backcross.     

Bayesian shrinkage analysis of quantitative trait Loci for dynamic traits

Many quantitative traits are measured repeatedly during the life of an organism. Such traits are called dynamic traits. The pattern of the changes of a dynamic trait is called the growth trajectory. Studying the growth trajectory may enhance our understanding of the genetic architecture of the growth trajectory. Recently, we developed an interval-mapping procedure to map QTL for dynamic traits under the maximum-likelihood framework. We fit the growth trajectory by Legendre polynomials. The method intended to map one QTL at a time and the entire QTL analysis involved scanning the entire genome by fitting multiple single-QTL models. In this study, we propose a Bayesian shrinkage analysis for estimating and mapping multiple QTL in a single model. The method is a combination between the shrinkage mapping for individual quantitative traits and the Legendre polynomial analysis for dynamic traits. The multiple-QTL model is implemented in two ways: (1) a fixed-interval approach where a QTL is placed in each marker interval and (2) a moving-interval approach where the position of a QTL can be searched in a range that covers many marker intervals. Simulation study shows that the Bayesian shrinkage method generates much better signals for QTL than the interval-mapping approach. We propose several alternative methods to present the results of the Bayesian shrinkage analysis. In particular, we found that the Wald test-statistic profile can serve as a mechanism to test the significance of a putative QTL.     

Bayesian model selection for genome-wide epistatic quantitative trait loci analysis

The problem of identifying complex epistatic quantitative trait loci (QTL) across the entire genome continues to be a formidable challenge for geneticists. The complexity of genome-wide epistatic analysis results mainly from the number of QTL being unknown and the number of possible epistatic effects being huge. In this article, we use a composite model space approach to develop a Bayesian model selection framework for identifying epistatic QTL for complex traits in experimental crosses from two inbred lines. By placing a liberal constraint on the upper bound of the number of detectable QTL we restrict attention to models of fixed dimension, greatly simplifying calculations. Indicators specify which main and epistatic effects of putative QTL are included. We detail how to use prior knowledge to bound the number of detectable QTL and to specify prior distributions for indicators of genetic effects. We develop a computationally efficient Markov chain Monte Carlo (MCMC) algorithm using the Gibbs sampler and Metropolis-Hastings algorithm to explore the posterior distribution. We illustrate the proposed method by detecting new epistatic QTL for obesity in a backcross of CAST/Ei mice onto M16i.     

Autosomal genome-wide linkage analysis to identify loci for gallbladder wall thickness in Mexican Americans

The significance of gallbladder wall thickness (GBWT) in regard to gallbladder disease (GBD) is not completely understood. Thickening of the gallbladder wall has been observed in patients with acute calculous and acalculous cholecystitis and chronic cholecystitis. However, various pathologic processes, such as gallbladder cancer and nonbiliary disorders such as liver cirrhosis and viral hepatitis, could also cause thickening of the gallbladder wall. To date, there is no report available on the genetic factors influencing GBWT. Therefore we sought to estimate the heritability (h2) of GBWT and to perform a genome-wide search to identify the susceptibility genes for GBWT, using data from the San Antonio Family Diabetes/Gallbladder Study (SAFDGS), a family study of Mexican Americans. GBWT was measured by ultrasound. After adjusting for the significant effects of age, sex, GBD (i.e., asymptomatic gallstones), metabolic syndrome, and duration of type 2 diabetes (T2DM), GBWT was found to be under significant and appreciable additive genetic influences (h2 +/- SE = 0.38 +/- 0.09, P < 0.0001). The strongest evidence for linkage occurred between markers D11S912 and D11S968 on chromosome 11q24-q25 (LOD = 2.7), where we have already shown suggestive evidence for linkage of GBD (LOD = 2.7) in a subset of our SAFDGS data. Potential evidence for linkage occurred at markers D1S1728 (1p31.1; LOD = 1.4) and D16S748 (16p13.1; LOD = 1.4), respectively. In conclusion, our study provides suggestive evidence for linkage of GBWT on chromosome 11q in Mexican Americans, and future tasks of mapping susceptibility gene(s) for GBD and its related traits, such as GBWT, in this chromosomal region can be fruitful.     

Models and tests of linkage and association studies of quantitative trait locus for multi-allele marker Loci

In this paper, we explore models and tests for association and linkage studies of a quantitative trait locus (QTL) linked to a multi-allele marker locus. Based on the difference between an offspring's conditional trait means of receiving and not receiving an allele from a parent at marker locus, we propose three statistics T(m), T(m,row) and T(m,col) to test association or linkage disequilibrium between the marker locus and the QTL. These tests are composite tests, and use the offspring marginal sample means including offspring data of both homozygous and heterozygous parents. For the linkage study, we calculate the offspring's conditional trait mean given the allele transmission status of a heterozygous parent at the marker locus. Based on the difference between the conditional means of a transmitted and a nontransmitted allele from a heterozygous parent, we propose statistics T(parsi), T(satur), T(gen) and T(m,het) to perform composite tests of linkage between the marker locus and the quantitative trait locus in the presence of association. These tests only use the offspring data that are related to the heterozygous parents at the marker locus. T(parsi) is a parsimonious or allele-wise statistic, T(satur) and T(gen )are satured or genotype-wise statistics, and T(m,het) compares the row and column sample means for offspring data of heterozygous parents. After comparing the powers and the sample sizes, we conclude that T(parsi) has higher power than those of the bi-allele tests, T(satur), T(gen), and T(m,het). If there is tight linkage between the marker and the trait locus, T(parsi) is powerful in detecting linkage between the marker and the trait locus in the presence of association. By investigating the goodness-of-fit of T(parsi), we find that T(satur) does not gain much power compared to that of T(parsi). Moreover, T(parsi) takes into account the pattern of the data that is consistent with linkage and linkage disequilibrium. As the number of alleles at the marker locus increases, T(parsi) is very conservative, and can be useful even for sparse data. To illustrate the usefulness and the power of the methods proposed in this paper, we analyze the chromosome 6 data of the Oxford asthma data, Genetic Analysis Workshop 12.     

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.     

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.     

Modeling quantitative trait Loci and interpretation of models

A quantitative genetic model relates the genotypic value of an individual to the alleles at the loci that contribute to the variation in a population in terms of additive, dominance, and epistatic effects. This partition of genetic effects is related to the partition of genetic variance. A number of models have been proposed to describe this relationship: some are based on the orthogonal partition of genetic variance in an equilibrium population. We compare a few representative models and discuss their utility and potential problems for analyzing quantitative trait loci (QTL) in a segregating population. An orthogonal model implies that estimates of the genetic effects are consistent in a full or reduced model in an equilibrium population and are directly related to the partition of the genetic variance in the population. Linkage disequilibrium does not affect the estimation of genetic effects in a full model, but would in a reduced model. Certainly linkage disequilibrium would complicate the detection of QTL and epistasis. Using different models does not influence the detection of QTL and epistasis. However, it does influence the estimation and interpretation of genetic effects.     

Sex-specific quantitative trait loci govern susceptibility to Theiler's murine encephalomyelitis virus-induced demyelination

Susceptibility to Theiler's murine encephalomyelitis virus-induced demyelination (TMEVD), a mouse model for multiple sclerosis (MS), is genetically controlled. Through a mouse-human comparative mapping approach, identification of candidate susceptibility loci for MS based on the location of TMEVD susceptibility loci may be possible. Composite interval mapping (CIM) identified quantitative trait loci (QTL) controlling TMEVD severity in male and female backcross populations derived from susceptible DBA/2J and resistant BALBc/ByJ mice. We report QTL on chromosomes 1, 5, 15, and 16 affecting male mice. In addition, we identified two QTL in female mice located on chromosome 1. Our results support the existence of three linked sex-specific QTL on chromosome 1 with opposing effects on the severity of the clinical signs of TMEV-induced disease in male and female mice.     

Segregation and linkage analysis for longitudinal measurements of a quantitative trait

We present a method for using slopes and intercepts from a linear regression of a quantitative trait as outcomes in segregation and linkage analyses. We apply the method to the analysis of longitudinal systolic blood pressure (SBP) data from the Framingham Heart Study. A first-stage linear model was fit to each subject's SBP measurements to estimate both their slope over time and an intercept, the latter scaled to represent the mean SBP at the average observed age (53.7 years). The subject-specific intercepts and slopes were then analyzed using segregation and linkage analysis. We describe a method for using the standard errors of the first-stage intercepts and slopes as weights in the genetic analyses. For the intercepts, we found significant evidence of a Mendelian gene in segregation analysis and suggestive linkage results (with LOD scores >or= 1.5) for specific markers on chromosomes 1, 3, 5, 9, 10, and 17. For the slopes, however, the data did not support a Mendelian model, and thus no formal linkage analyses were conducted.     

A genome-wide linkage scan for quantitative trait loci underlying obesity related phenotypes in 434 Caucasian families

To identify quantitative trait loci (QTLs) that contribute to obesity, we performed a large-scale whole genome linkage scan (WGS) involving 4,102 individuals from 434 Caucasian families. The most pronounced linkage evidence was found at the genomic region 20p11-12 for fat mass (LOD = 3.31) and percentage fat mass (PFM) (LOD = 2.92). We also identified several regions showing suggestive linkage signals (threshold LOD = 1.9) for obesity phenotypes, including 5q35, 8q13, 10p12, and 17q11.     

Optimal sibship selection for genotyping in quantitative trait locus linkage analysis

In this paper we present a novel method for selecting optimally informative sibships of any size for quantitative trait locus (QTL) linkage analysis. The method allocates a quantitative index of potential informativeness to each sibship on the basis of observed trait scores and an assumed true QTL model. Any sample of phenotypically screened sibships can therefore be easily rank-ordered for selective genotyping. The quantitative index is the sibship's expected contribution to the non-centrality parameter. This expectation represents the weighted sum of chi(2) test statistics that would be obtained given the observed trait values over all possible sibship genotypic configurations; each configuration is weighted by the likelihood of it occurring given the assumed true genetic model. The properties of this procedure are explored in relation to the accuracy of the assumed true genetic model and sibship size. In comparison to previous methods of selecting phenotypically extreme sibships for genotyping, the proposed method is considerably more efficient and is robust with regard to the specification of the genetic model.     

Theobroma cacao L.: a genetic linkage map and quantitative trait loci analysis

A genetic linkage map of Theobroma cacao (cocoa) has been constructed from 131 backcross trees derived from a cross between a single tree of the variety Catongo and an F1 tree from the cross of Catongo by Pound 12. The map comprises 138 markers: 104 RAPD loci, 32 RFLP loci and two morphologic loci. Ten linkage groups were found which cover 1068 centimorgans (cM). Only six (4%) molecular-marker loci show a significant deviation from the expected 11 segregation ratio.The average distance between two adjacent markers is 8.3 cM. The final genome-size estimates based on two-point linkage data ranged from 1078 to 1112 cM for the cocoa genome. This backcross progeny segregates for two apparently single gene loci controlling (1) anthocyanidin synthesis (Anth) in seeds, leaves and flowers and (2) self-compatibility (Autoc). The Anth locus was found to be 25 cM from Autoc and two molecular markers co-segregate with Anth. The genetic linkage map was used to localize QTLs for early flowering, trunk diameter, jorquette height and ovule number in the BC₁ generation using both single-point ANOVA and interval mapping. A minimum number of 2–4 QTLs (P<0.01) involved in the genetic expression of the traits studied was detected. Coincident map locations of a QTL for jorquette height and trunk diameter suggests the possibility of pleiotropic effects in cocoa for these traits. The combined estimated effects of the different mapped QTLs explained between 11.2% and 25.8% of the phenotypic variance observed in the BC₁ population.     

Joint mapping of quantitative trait Loci for multiple binary characters

Joint mapping for multiple quantitative traits has shed new light on genetic mapping by pinpointing pleiotropic effects and close linkage. Joint mapping also can improve statistical power of QTL detection. However, such a joint mapping procedure has not been available for discrete traits. Most disease resistance traits are measured as one or more discrete characters. These discrete characters are often correlated. Joint mapping for multiple binary disease traits may provide an opportunity to explore pleiotropic effects and increase the statistical power of detecting disease loci. We develop a maximum-likelihood method for mapping multiple binary traits. We postulate a set of multivariate normal disease liabilities, each contributing to the phenotypic variance of one disease trait. The underlying liabilities are linked to the binary phenotypes through some underlying thresholds. The new method actually maps loci for the variation of multivariate normal liabilities. As a result, we are able to take advantage of existing methods of joint mapping for quantitative traits. We treat the multivariate liabilities as missing values so that an expectation-maximization (EM) algorithm can be applied here. We also extend the method to joint mapping for both discrete and continuous traits. Efficiency of the method is demonstrated using simulated data. We also apply the new method to a set of real data and detect several loci responsible for blast resistance in rice.     

Genetic linkage analysis of a dichotomous trait incorporating a tightly linked quantitative trait in affected sib pairs

Many complex diseases are usually considered as dichotomous traits but are also associated with quantitative biological markers or quantitative risk factors. For such dichotomous traits, although their associated quantitative traits may not directly underly the diagnosis of the disease status, if the associated quantitative trait is also linked to the chromosomal regions linked to the dichotomous trait, then joint analysis of dichotomous and quantitative traits should be more efficient than consideration of them separately. Previous studies have focused on the situation when a dichotomous trait can be modeled by a threshold process acting on a single underlying normal liability distribution. However, for many complex disorders, including most psychiatric disorders, diagnosis is generally based on a set of binary or discrete criteria. These traits cannot be modeled on the basis of a threshold process acting on an underlying continuous trait. We propose a likelihood-based method that efficiently combines such a discrete trait and an associated quantitative trait in the analysis, using affected-sib-pair data. Our simulation studies suggest that joint analysis increases the power to detect linkage of dichotomous traits. We also apply the proposed new method to an asthma genome-scan data set and incorporate the total serum immunoglobulin E level in the analysis.     

GCTA: a tool for genome-wide complex trait analysis

For most human complex diseases and traits, SNPs identified by genome-wide association studies (GWAS) explain only a small fraction of the heritability. Here we report a user-friendly software tool called genome-wide complex trait analysis (GCTA), which was developed based on a method we recently developed to address the "missing heritability" problem. GCTA estimates the variance explained by all the SNPs on a chromosome or on the whole genome for a complex trait rather than testing the association of any particular SNP to the trait. We introduce GCTA's five main functions: data management, estimation of the genetic relationships from SNPs, mixed linear model analysis of variance explained by the SNPs, estimation of the linkage disequilibrium structure, and GWAS simulation. We focus on the function of estimating the variance explained by all the SNPs on the X chromosome and testing the hypotheses of dosage compensation. The GCTA software is a versatile tool to estimate and partition complex trait variation with large GWAS data sets.     

Variance components models for gene-environment interaction in quantitative trait locus linkage analysis.

Gene-environment interaction (G x E) is likely to be a common and important source of variation for complex behavioral traits. Gene-environment interaction, or genetic control of sensitivity to the environment, can be incorporated into variance components twin and sib-pair analyses by partitioning genetic effects into a mean part, which is independent of the environment, and a part that is a linear function of the environment. An approach described in a companion paper (Purcell, 2002) is applied to sib-pair variance components linkage analysis in two ways: allowing for quantitative trait locus by environment interaction and utilizing information on any residual interactions detected prior to analysis. As well as elucidating environmental pathways, consideration of G x E in quantitative and molecular studies will potentially direct and enhance gene-mapping efforts.     

Identification of plant quantitative trait Loci modulating a rhizobacteria-aphid indirect effect

Plants simultaneously interact with a plethora of species both belowground and aboveground, which can result in indirect effects mediated by plants. Studies incorporating plant genetic variation indicate that indirect effects mediated by plants may be a significant factor influencing the ecology and evolution of species within a community. Here, we present findings of a Quantitative Trait Locus (QTL) mapping study, where we mapped a rhizobacteria-aphid indirect effect onto the barley genome. We measured the size of aphid populations on barley when the barley rhizosphere either was or was not supplemented with a rhizobacterial species. Using a QTL mapping subset, we located five regions of the barley genome associated with the rhizobacteria-aphid indirect effect. Rhizobacterial supplementation led to an increase in aphid population size (mapped to three barley QTL), or a decrease in aphid population size (mapped to two barley QTL). One QTL associated with plant resistance to aphids was affected by a significant QTL-by-environment interaction, because it was not expressed when rhizobacteria was supplemented. Our results indicated that rhizobacterial supplementation of barley roots led to either increased or reduced aphid population size depending on plant genotype at five barley QTL. This indicates that the direction of a rhizobacteria-aphid indirect effect could influence the selection pressure on plants, when considering species that affect plant fitness. Further research may build on the findings presented here, to identify genes within QTL regions that are involved in the indirect interaction.     

Using progenitor strain information to identify quantitative trait nucleotides in outbred mice

We have developed a fast and economical strategy for dissecting the genetic architecture of quantitative trait loci at a molecular level. The method uses two pieces of information: mapping data from crosses that involve more than two inbred strains and sequence variants in the progenitor strains within the interval containing a quantitative trait locus (QTL). By testing whether the strain distribution pattern in the progenitor strains is consistent with the observed genetic effect of the QTL we can assign a probability that any sequence variant is a quantitative trait nucleotide (QTN). It is not necessary to genotype the animals except at a skeleton of markers; the genotypes at all other polymorphisms are estimated by a multipoint analysis. We apply the method to a 4.8-Mb region on mouse chromosome 1 that contains a QTL influencing anxiety segregating in a heterogeneous stock and show that, under the assumption that a single QTN is present and lies in a region conserved between the human and mouse genomes, it is possible to reduce the number of variants likely to be the quantitative trait nucleotide from many thousands to <20.