Localization of genes involved in the metabolic syndrome using multivariate linkage analysis
© Olswold and Andrade; licensee BioMed Central Ltd 2003
Published: 31 December 2003
There are no well accepted criteria for the diagnosis of the metabolic syndrome. However, the metabolic syndrome is identified clinically by the presence of three or more of these five variables: larger waist circumference, higher triglyceride levels, lower HDL-cholesterol concentrations, hypertension, and impaired fasting glucose. We use sets of two or three variables, which are available in the Framingham Heart Study data set, to localize genes responsible for this syndrome using multivariate quantitative linkage analysis. This analysis demonstrates the applicability of using multivariate linkage analysis and how its use increases the power to detect linkage when genes are involved in the same disease mechanism.
It has been shown that for correlated traits, multivariate approaches for genetic linkage analyses can increase the power and precision to identify genetic effects [1–4]. When correlated measures are considered, the composite score from joint consideration of all measures reflects a smaller level of measurement error than each of the univariate measures . Then, multivariate analysis provides a statistically efficient mechanism for controlling the analysis-wise significance level when there are multiple trait observations for each subject [3, 6]. Therefore, using methods that can analyze several traits jointly is likely to enhance the ability to identify genes influencing the metabolic syndrome. Although multivariate Haseman-Elston (H-E)  and variance-components (VC) methods  have been available for several years, only recently has the power of these methods been compared. Allison et al.  presented results from a large simulation study to assess the effectiveness of a bivariate H-E test for linkage versus the univariate H-E test . Their results showed that bivariate analyses can improve the power to detect linkage, with a greater gain in power when the genetic covariance due to a major locus linked to the marker studied is negative and the residual covariance among the traits is positive. Amos et al.  also showed that bivariate approaches are more powerful than univariate analyses except for traits with very high positive polygenic correlation. Evans  also reached similar conclusion.
Our approach is based on the assumption that it is easier to detect a quantitative trait locus (QTL) involved in the metabolic syndrome using multivariate linkage analysis. Our aim is to show that using combinations of traits related to the metabolic syndrome, and then using them in multivariate linkage analysis software, gives reliable results for linkage to genes associated with this syndrome.
The metabolic syndrome
Clinical identification of the metabolic syndrome
> 102 cm (> 40 in)
> 88 cm (> 35 in)
= 150 mg/dL
< 40 mg/dL
< 50 mg/dL
= 130/85 mm Hg
Multivariate linkage analysis
The multivariate variance-components (MVC) approach is an extension of the univariate approach described by Amos . For multivariate traits, let Y i = (Y11,...,Y1ki,...,Y mki )' be a vector of m multivariate trait values for k i members of the ith family. Let N be the total number of families, β a vector of dimension mp of the regression coefficients for the p covariates (including a vector of 1's corresponding to the overall mean), X i = I m X ki x m an mk i × mp known matrix of covariate values for the ith family, where is the Kronecker product, and V i a VC matrix of dimension mk i × mk i . Then, the variance-covariance matrix of the traits is V i = A G i + B Z i + C I i , where G i is the k i × k i matrix of the coefficients of relationship for the family i; Z i an k i × k i matrix of estimated proportion of alleles identical by decent (IBD) for pairs of related individuals for the ith pedigree; I i is the k i × k i identity matrix; and A, B, and C, are, respectively, polygenic, major-gene, and environment variance-covariance matrices each of dimension m × m. A more detailed description of these models was presented elsewhere [11, 12].
Multivariate VC test
To test for genetic linkage, we also construct a likelihood ratio test. Under the null hypothesis, the major gene parameter(s) are restricted to equal 0. The distribution of the multivariate test is a mixture of χ2 values . For trivariate linkage analysis of an additive genetic effect, the distribution of the trivariate test that the major-gene covariance components are zero is a mixture of 1/8 χ02, 3/8 χ12, 3/8 χ32 and 1/8 χ62. One-eighth of the time all the VCs are estimated to be positive with all the covariances different from 0 yielding 6 degrees of freedom. Three-eighths of the time, one of the VCs is estimated to be zero with two covariances fixed to zero (yielding 3 degrees of freedom). Another three-eighths of the time two VCs are fixed to zero with all covariances equal to zero yielding 1 degree of freedom. Finally, one-eighth of the time all the variances are fixed to zero resulting in a degenerate distribution of point mass at zero.
For the multivariate linkage analysis, we use the following four traits: triglycerides, HDL-cholesterol, systolic blood pressure (SBP), and fasting glucose. Since these variables, except for triglycerides, were measured at several time points, we applied a similar regression approach described in Levy et al.  for these four variables and then used their residuals as the quantitative traits in the multivariate genome-wide linkage analysis for quantitative traits. There are two packages that use the MVC approach: ACT  and EMVC . The analyses here presented were performed using the EMVC package using 330 families with 4692 individuals, of whom 1702 have genotype information.
Values of polygenic and QTL correlation between the variables involved in the metabolic syndrome at different locations
QTL Effects for Traits
Polygenic Effects for Traits
TraitsA (1, 2, 3)
1 and 2B
1 and 3
2 and 3
1 and 2C
1 and 3
2 and 3
The MVC approach appears to perform well in the identification of regions linked to genes associated with traits related to the metabolic syndrome, mainly on regions where the QTL effects were negatively correlated and there was a positively correlated polygenic effect as shown by Amos et al.  and Evans . Our results did identify a minor linkage peak to the same region of chromosome 17 described by Levy et al. . The only region on chromosome 17 using the trivariate VC approach that showed evidence for linkage was on the surrounding region of 10 cM, which was due primarily to the bivariate combination, SBP and triglycerides, (LOD 3.14, position 12 cM, results not shown). Furthermore, evidence for linkage was also found on chromosomes 2, 5, and 6. We also showed that the pair-wise combinations with evidence for linkage are the ones that have either small to moderate genetic correlation or negative genetic correlation. In summary, the use of multivariate quantitative trait loci linkage analysis can increase the power to detect a QTL. However, this procedure is computationally intensive, i.e., the CPU time increases exponentially as the number of traits increases additively.
The authors would like to thank Brooke Fridley and Beth Atkinson for their help and two anonymous reviewers for their helpful comments. This research was partially funded by NIH grant R01HL71917.
- Martin N, Boomsma D, Machin G: A twin-pronged attack on complex traits. Nat Genet. 1997, 17: 387-392. 10.1038/ng1297-387.View ArticlePubMedGoogle Scholar
- Boomsma DI, Dolan CV: A comparison of power to detect a QTL in sib-pair data using multivariate phenotypes, mean phenotypes, and factor scores. Behav Genet. 1998, 28: 329-340. 10.1023/A:1021665501312.View ArticlePubMedGoogle Scholar
- Amos CI, de Andrade M, Zhu D: Comparison of multivariate tests for genetic linkage. Hum Hered. 2001, 51: 133-144. 10.1159/000053334.View ArticlePubMedGoogle Scholar
- Evans DM: The power of multivariate quantitative-trait loci linkage analysis is influenced by the correlation between the variables. Am J Hum Genet. 2002, 70: 1599-1602. 10.1086/340850.PubMed CentralView ArticlePubMedGoogle Scholar
- Schmitz S, Cherny SS, Fulker DW: Increase in power through multivariate analyses. Behav Genet. 1998, 28: 357-363. 10.1023/A:1021669602220.View ArticlePubMedGoogle Scholar
- Allison DB, Thiel B, St Jean P, Elston RC, Infante MC, Schork NJ: Multiple phenotype modeling in gene-mapping studies of quantitative traits: power advantages. Am J Hum Genet. 1998, 63: 1190-1201. 10.1086/302038.PubMed CentralView ArticlePubMedGoogle Scholar
- Amos CI, Elston RC, Bonney GE, Keats BJB, Berenson GS: A multivariate method for detecting genetic linkage with application to the study of a pedigree with an adverse lipoprotein phenotype. Am J Hum Genet. 1990, 47: 247-54.PubMed CentralPubMedGoogle Scholar
- Amos CI: Robust variance-components approaches for assessing genetic linkage in pedigrees. Am J Hum Genet. 1994, 54: 535-543.PubMed CentralPubMedGoogle Scholar
- Haseman JK, Elston RC: The investigation of linkage between a quantitative trait and a marker locus. Behav Genet. 1972, 2: 3-19. 10.1007/BF01066731.View ArticlePubMedGoogle Scholar
- National Cholesterol Education Program: Executive Summary of the Third Report of The National Cholesterol Education Program (NCEP) Expert Panel on Detection, Evaluation, and Treatment of High Blood Cholesterol In Adults (Adult Treatment Panel III). JAMA. 2001, 285: 2486-2497. 10.1001/jama.285.19.2486.View ArticleGoogle Scholar
- Almasy L, Dyer TD, Blangero J: Bivariate quantitative trait linkage analysis: pleiotropy versus co-incident linkages. Genet Epidemiol. 1997, 14: 953-958. 10.1002/(SICI)1098-2272(1997)14:6<953::AID-GEPI65>3.0.CO;2-K.View ArticlePubMedGoogle Scholar
- de Andrade M, Thiel TJ, Yu L, Amos CI: Assessing linkage in chromosome 5 using components of variance approach: univariate versus multivariate. Genet Epidemiol. 1997, 14: 773-778. 10.1002/(SICI)1098-2272(1997)14:6<773::AID-GEPI35>3.0.CO;2-L.View ArticlePubMedGoogle Scholar
- Self SG, Liang K-L: Asymptotic properties of maximum likelihood estimators and likelihood ratio tests under nonstandard conditions. J Am Stat Assoc. 1987, 82: 605-610. 10.2307/2289471.View ArticleGoogle Scholar
- Levy D, DeStefano AL, Larson MG, O'Donnell CJ, Lifton RP, Gavras H, Cupples LA, Myers RH: Evidence for a gene influencing blood pressure on chromosome 17. Genome scan linkage results for longitudinal blood pressure phenotypes in subjects from the Framingham Heart Study. Hypertension. 2000, 6: 477-483.View ArticleGoogle Scholar
- de Andrade M, Krushkal J, Yu L, Zhu D, Amos CI: ACT – A computer package for analysis of complex traits. Denver, CO. 1998Google Scholar
- Iturria SJ, Blangero J: An EM algorithm for obtaining maximum likelihood estimates in the multi-phenotype variance components linkage model. Ann Hum Genet. 2000, 64: 349-369. 10.1046/j.1469-1809.2000.6440349.x.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.