Joint tests of linkage and association for quantitative traits

In this paper we describe various study designs and analytic techniques for testing the joint hypothesis that a genetic marker is both linked to and associated with a quantitative phenotype. Issues of power and sampling are addressed. The distinction between methods that explicitly examine association and those that infer association by examining the distribution of allelic transmissions from a heterozygous parent is examined. Extensions to multivariate, multiallelic, and multilocus situations are addressed. Recent approaches that combine variance-components-based linkage analyses with joint tests of linkage in the presence of association for disentanglement of the linkage and association and the application of such methods to fine mapping are discussed. Finally, new classes of joint tests of linkage and association that do not require samples of related individuals are described.     

Sibling-based tests of linkage and association for quantitative traits.

The transmission/disequilibrium test (TDT) developed by Spielman et al. can be a powerful family-based test of linkage and, in some cases, a test of association as well as linkage. It has recently been extended in several ways; these include allowance for implementation with quantitative traits, allowance for multiple alleles, and, in the case of dichotomous traits, allowance for testing in the absence of parental data. In this article, these three extensions are combined, and two procedures are developed that offer valid joint tests of linkage and (in the case of certain sibling configurations) association with quantitative traits, with use of data from siblings only, and that can accommodate biallelic or multiallelic loci. The first procedure uses a mixed-effects (i.e., random and fixed effects) analysis of variance in which sibship is the random factor, marker genotype is the fixed factor, and the continuous phenotype is the dependent variable. Covariates can easily be accommodated, and the procedure can be implemented in commonly available statistical software. The second procedure is a permutation-based procedure. Selected power studies are conducted to illustrate the relative power of each test under a variety of circumstances.     

Combined linkage and association sib-pair analysis for quantitative traits.

An extension to current maximum-likelihood variance-components procedures for mapping quantitative-trait loci in sib pairs that allows a simultaneous test of allelic association is proposed. The method involves modeling of the allelic means for a test of association, with simultaneous modeling of the sib-pair covariance structure for a test of linkage. By partitioning of the mean effect of a locus into between- and within-sibship components, the method controls for spurious associations due to population stratification and admixture. The power and efficacy of the method are illustrated through simulation of various models of both real and spurious association.     

Testing quantitative traits for association and linkage in the presence or absence of parental data

Zhu and Elston developed a transmission disequilibrium test for quantitative traits by defining a linear transformation to condition out founder information. The method tests the null hypothesis of no linkage or association and can be applied to general pedigree structures. However, this method requires both genotype and phenotype parental information, which may be difficult to obtain. In this paper, we describe parametric and non-parametric methods to relax this requirement when only nuclear families are sampled. We show that neither method is affected by population stratification in the absence of linkage. The statistical power and validity of the tests are investigated by simulation. A simple simulation method to calculate the power of the nonparametric method is also discussed. In practice, the data may have some families with parental phenotype and genotype information available and some without. We briefly discuss how all the data may be analyzed jointly.     

Linkage analysis and association analysis in the presence of linkage using age at onset of COGA alcoholism data

Complex disease mapping usually involves a combination of linkage and association techniques. Linkage analysis can scan the entire genome in a few hundred tests. Association tests may involve an even greater number of tests. However, association tests can localize the susceptibility genes more accurately. Using a recently developed combined linkage and association strategy, we analyzed a subset of the Collaborative Study on the Genetics of Alcoholism (COGA) data for the Genetic Analysis Workshop 14 (GAW14). In this analysis, we first employed linkage analysis based on frailty models that take into account age of onset information to establish which regions along the chromosome are likely to harbor disease susceptibility genes for alcohol dependence. Second, we used an association analysis by exploiting linkage disequilibrium to narrow down the peak regions. We also compare the methods with mean identity-by-descent tests and transmission/disequilibrium tests that do not use age of onset information.     

Genome-wide linkage and peak-wide association study of obesity-related quantitative traits in Caribbean Hispanics

Although obesity is more prevalent in Hispanics than non-Hispanic whites in the United States, little is known about the genetic etiology of the related traits in this population. To identify genetic loci influencing obesity in non-Mexican Hispanics, we performed a genome-wide linkage scan in 1,390 subjects from 100 Caribbean Hispanic families on six obesity-related quantitative traits: body mass index (BMI), body weight, waist circumference, waist-to-hip ratio, abdominal and average triceps skinfold thickness after adjusting for significant demographic and lifestyle factors. We then carried out an association analysis of the linkage peaks and the FTO gene in an independent community-based Hispanic subcohort (N = 652, 64% Caribbean Hispanics) from the Northern Manhattan Study. Evidence of linkage was strongest on 1q43 with multipoint LOD score of 2.45 (p = 0.0004) for body weight. Suggestive linkage evidence of LOD > 2.0 was also identified on 1q43 for BMI (LOD = 2.03), 14q32 for abdominal skinfold thickness (LOD = 2.17), 16p12 for BMI (LOD = 2.27) and weight (LOD = 2.26), and 16q23–24 for average triceps skinfold thickness (LOD = 2.32). In the association analysis of 6,440 single nucleotide polymorphisms (SNPs) under 1-LOD unit down regions of our linkage peaks on chromosome 1q43 and 16p12 as well as in the FTO gene, we found that two SNPs (rs6665519 and rs669231) on 1q43 and one FTO SNP (rs12447427) were significantly associated with BMI or body weight after adjustment for multiple testing. Our results suggest that in addition to FTO, multiple genetic loci, particularly those on 1q43 region, may contribute to the variations in obesity-related quantitative traits in Caribbean Hispanics.     

Bivariate association analysis for quantitative traits using generalized estimation equation

Quantitative traits often underlie risk for complex diseases. Many studies collect multiple correlated quantitative phenotypes and perform univariate analyses on each of them respectively. However, this strategy may not be powerful and has limitations to detect plei- otropic genes that may underlie correlated quantitative traits. In addition, testing multiple traits individually will exacerbate perplexing problem of multiple testing. In this study, generalized estimating equation 2 （GEE2） is applied to association mapping of two correlated quantitative traits. We suppose that a quantitative trait locus is located in a chromosome region that exerts pleiotropic effects on multiple quantitative traits. In that region, multiple SNPs are genotyped. Genotypes of these SNPs and the two quantitative traits affected by a causal SNP were simulated under various parameter values： residual correlation coefficient between two traits, causal SNP heritability, minor allele frequency of the causal SNP, extent of linkage disequilibrium with the causal SNP, and the test sample size. By power ana- lytical analyses, it is showed that the bivariate method is generally more powerful than the univariate method. This method is robust and yields false-positive rates close to the pre-set nominal significance level. Our real data analyses attested to the usefulness of the method.     

Multivariate linkage analysis using the electrophysiological phenotypes in the COGA alcoholism data

Multivariate linkage analysis using several correlated traits may provide greater statistical power to detect susceptibility genes in loci whose effects are too small to be detected in univariate analysis. In this analysis, we apply a new approach and perform a linkage analysis of several electrophysiological phenotypes of the Collaborative Study on the Genetics of Alcoholism data of the Genetic Analysis Workshop 14. Our approach is based on a variance-component model to map candidate genes using repeated or longitudinal measurements. It can take into account covariate effects and time-dependent genetic effects in general pedigree data. We compare our results with the ones obtained by SOLAR using single measurement data. Our multivariate linkage analysis found linkage evidence on two regions on chromosome 4: around marker GABRB1 at 51.4 cM and marker FABP2 at 116.8 cM (unadjusted p-value = 0.00006).     

Comparison of marker types and map assumptions using Markov chain Monte Carlo-based linkage analysis of COGA data

We performed multipoint linkage analysis of the electrophysiological trait ECB21 on chromosome 4 in the full pedigrees provided by the Collaborative Study on the Genetics of Alcoholism (COGA). Three Markov chain Monte Carlo (MCMC)-based approaches were applied to the provided and re-estimated genetic maps and to five different marker panels consisting of microsatellite (STRP) and/or SNP markers at various densities. We found evidence of linkage near the GABRB1 STRP using all methods, maps, and marker panels. Difficulties encountered with SNP panels included convergence problems and demanding computations.     

Recursive partitioning models for linkage in COGA data

We have developed a recursive-partitioning (RP) algorithm for identifying phenotype and covariate groupings that interact with the evidence for linkage. This data-mining approach for detecting gene x environment interactions uses genotype and covariate data on affected relative pairs to find evidence for linkage heterogeneity across covariate-defined subgroups. We adapted a likelihood-ratio based test of linkage parameterized with relative risks to a recursive partitioning framework, including a cross-validation based deviance measurement for choosing optimal tree size and a bootstrap sampling procedure for choosing robust tree structure. ALDX2 category 5 individuals were considered affected, categories 1 and 3 unaffected, and all others unknown. We sampled non-overlapping affected relative pairs from each family; therefore, we used 144 affected pairs in the RP model. Twenty pair-level covariates were defined from smoking status, maximum drinks, ethnicity, sex, and age at onset. Using the all-pairs score in GENEHUNTER, the nonparametric linkage tests showed no regions with suggestive linkage evidence. However, using the RP model, several suggestive regions were found on chromosomes 2, 4, 6, 14, and 20, with detection of associated covariates such as sex and age at onset.     

A sib-pair regression model of linkage disequilibrium for quantitative traits

A multiple-regression model is described for the detection of linkage disequilibrium in quantitative trait loci. The model is developed for application to large numbers of single nucleotide polymorphism (SNP) markers genotyped on small nuclear families. Parental data are not required by the method, although it provides a direct means to test quantitative trait locus-marker allele association and to determine whether any such association is attributable to linkage disequilibrium or population admixture. Analytical expectations for the regression coefficients are derived, allowing direct interpretation of the parameter estimates. Simulation studies indicate a substantial improvement in power over classical linkage studies of sibling pairs and show the effects of population admixture on the model outcomes.     

Comparison of a multipoint identity-by-descent method with parametric multipoint linkage analysis for mapping quantitative traits.

We previously developed a method of partitioning genetic variance of a quantitative trait to loci in specific chromosomal regions. In this paper, we compare this method--multipoint IBD (identical by descent) method (MIM)--with parametric multipoint linkage analysis (MLINK). A simulation study was performed comparing the methods for the major-locus, mixed, and two-locus models. The criterion for comparisons between MIM and MLINK was the average lod score from multiple replicates of simulated data sets. The effect of gene frequency, dominance, model misspecification, marker spacing, and informativeness are also considered in a smaller set of simulations. Within the context of the models examined, the MIM approach was found to be comparable in power with parametric multipoint linkage analysis when (a) parental data are unknown, (b) the effect of the major locus is small and there is additional genetic variation, or (c) the parameters of the major-locus model are misspecified. The performance of the MIM method relative to MLINK was markedly lower when the allele frequency at the trait locus was .2 versus .5, particularly for the case when parental data were assumed to be known. Dominance at the trait major locus, as well as marker spacing and heterozygosity, did not appear to have a large effect on the ELOD comparisons.     

Accounting for expected linkage in biometric analysis of quantitative traits

Russian Journal of Genetics - The problem of accounting for a genetic estimation of expected linkage in the disposition of random loci was solved for the additive-dominant model. The...     

Genetic dissection of wheat panicle traits using linkage analysis and a genome-wide association study

Key message

Coincident regions on chromosome 4B for GW, on 5A for SD and TSS, and on 3A for SL and GNS were detected through an integration of a linkage analysis and a genome-wide association study (GWAS). In addition, six stable QTL clusters on chromosomes 2D, 3A, 4B, 5A and 6A were identified with high PVE% on a composite map.

Abstract

The panicle traits of wheat, such as grain number per spike and 1000-grain weight, are closely correlated with grain yield. Superior and effective alleles at loci related to panicles developments play a crucial role in the progress of molecular improvement in wheat yield breeding. Here, we revealed several notable allelic variations of seven panicle-related traits through an integration of genome-wide association mapping and a linkage analysis. The linkage analysis was performed using a recombinant inbred line (RIL) population (173 lines of F_8:9) with a high-density genetic map constructed with 90K SNP arrays, Diversity Arrays Technology (DArT) and simple sequence repeat (SSR) markers in five environments. Thirty-five additive quantitative trait loci (QTL) were discovered, including eleven stable QTLs on chromosomes 1A, 2D, 4B, 5B, 6B, and 6D. The marker interval between EX_C101685 and RAC875_C27536 on chromosome 4B exhibited pleiotropic effects for GW, SL, GNS, FSN, SSN, and TSS, with the phenotypic variation explained (PVE) ranging from 5.40 to 37.70%. In addition, an association analysis was conducted using a diverse panel of 205 elite wheat lines with a composite map (24,355 SNPs) based on the Illumina Infinium assay in four environments. A total of 73 significant marker-trait associations (MTAs) were detected for panicle traits, which were distributed across all wheat chromosomes except for 4D, 5D, and 6D. Consensus regions between RAC875_C27536_611 and Tdurum_contig4974_355 on chromosome 4B for GW in multiple environments, between QTSS5A.7-43 and BS00021805_51 on 5A for SD and TSS, and between QSD3A.2-164 and RAC875_c17479_359 on 3A for SL and GNS in multiple environments were detected through linkage analysis and a genome-wide association study (GWAS). In addition, six stable QTL clusters on chromosomes 2D, 3A, 4B, 5A, and 6A were identified with high PVE% on a composite map. This study provides potentially valuable information on the dissection of yield-component traits and valuable genetic alleles for molecular-design breeding or functional gene exploration.

     

A unified Haseman-Elston method for testing linkage with quantitative traits

The Haseman and Elston (H-E) method uses a simple linear regression to model the squared trait difference of sib pairs with the shared allele identical by descent (IBD) at marker locus for linkage testing. Under this setting, the squared mean-corrected trait sum is also linearly related to the IBD sharing. However, the resulting slope estimate for either model is not efficient. In this report, we propose a simple linkage test that optimally uses information from the estimates of both models. We also demonstrate that the new test is more powerful than both the traditional one and the recently revisited H-E methods.     

Analysis of alcohol dependence phenotype in the COGA families using covariates to detect linkage

Linkage analysis methods that incorporate etiological heterogeneity of complex diseases are likely to demonstrate greater power than traditional linkage analysis methods. Several such methods use covariates to discriminate between linked and unlinked pedigrees with respect to a certain disease locus. Here we apply several such methods including two mixture models, ordered subset analysis, and a conditional logistic model to genome scan data on the DSM-IV alcohol dependence phenotype on the Collaborative Studies on Genetics of Alcoholism families, and compare the results to traditional nonparametric linkage analysis. In general, there was little agreement among the various covariate-based linkage statistics. Linkage signals with empirical p-values less than 0.001 were detected on chromosomes 3, 4, 7, 10, and 12, with the highest peak occurring at the GABRB1 gene using the ecb21 covariate.     

Mapping quantitative traits and strategies to find quantitative trait genes

In 1999 a meeting took place at the Jackson Laboratory, a large mouse research centre in Bar Harbor, Maine, to consider the value of systematically collecting phenotypes on inbred strains of mice (Paigen and Eppig (2000) [1]). The group concluded that cataloguing the extensive phenotypic diversity present among laboratory mice, and in particular providing the research community with data from cohorts of animals, phenotyped according to standardized protocols, was essential if we were to take advantage of the possibilities of mouse genetics. Beginning with the collection of basic physiological, biochemical and behavioral data on nine commonly used inbred strains, the project has expanded so that by the beginning of 2010 data for 178 strains had been collected, with 105 phenotype projects yielding over 2000 different measurements (Bogue et al. (2007) [2].