- Proceedings
- Open access
- Published:

# Multivariate variance-components analysis of longitudinal blood pressure measurements from the Framingham Heart Study

*BMC Genetics*
**volume 4**, Article number: S55 (2003)

## Abstract

Multivariate variance-components analysis provides several advantages over univariate analysis when studying correlated traits. It can test for pleiotropy or (in the longitudinal context) gene × age interaction. It can also have more power than univariate analyses to detect a quantitative trait locus influencing several traits. We apply multivariate variance components to longitudinal systolic blood pressure data from the Framingham Heart Study. We find evidence for a polygenic influence on blood pressure (heritabilities at different ages range from 27% to 38%). Tests based on a factor-analytic parameterization of the polygenic variance find significant (*p* < 2 × 10^{-3}) evidence that different genes affect blood pressure at different ages. Still, estimates for the proportion of polygenic variance due to shared genes ran as high as 85% for some trait pairs. Univariate and multivariate linkage analyses replicate previous linkage results on chromosome 17 (maximum LOD scores of 2.2 and 2.4, respectively). In this study, multivariate analysis provides no increase in power; this is likely due to the strong positive correlation in systolic blood pressure measured at different ages.

## Background

High blood pressure is a complex disorder that results from environmental and genetic factors and their interactions. Levy et al. [1] found evidence for a gene influencing blood pressure on chromosome 17 using data from the Framingham Heart Study. However, this study analyzed average blood pressure over a 50-year period (ages 25 to 75), and may not have taken full advantage of the longitudinal nature of the Framingham study. Blood pressure increases with age; there may be genes that influence the rate of this increase. Similarly, there may be genes that influence blood pressure only at early or late ages. For example, a segregation analysis by Pérusse et al. [2] suggested that blood pressure is influenced by a major gene with age-dependent effects. Animal studies have also found that different genes can influence a trait at different ages [3]. Taking lifetime averages may mask such effects.

de Andrade et al. [4] recently analyzed longitudinal quantitative trait data using a multivariate variance components approach. This approach can be more powerful for correlated traits than a univariate approach [5]. It can also test for gene × time interaction, polygenic pleiotropy – defined in the longitudinal context as a trait being determined by the same set of genes at distinct time points – and distinguish between major gene pleiotropy and co-incident linkage [6–8].

de Andrade et al. [4] defined traits by calendar time. Thus, for a cohort study with age-staggered entry, each measurement will have been taken at approximately the same calendar time for all subjects, but at different biological ages. In the context of the Framingham study we propose to define traits by biological age, so as to distinguish genes involved in determining high blood pressure at young or old ages as opposed to uncovering a gene × calendar time (environment) interaction.

We apply univariate and multivariate variance components to systolic blood pressure measurements on subjects from the Framingham Heart Study taken in four different age ranges. We consider models with a polygenic component and both a polygenic and major gene component. We describe and apply a test of the null hypothesis of complete pleiotropy versus the alternative of incomplete pleiotropy based on a factor-analytic parameterization of the polygenic variance component. We reject the null hypothesis of complete pleiotropy, suggesting a different set of genes influence systolic blood pressure at different ages. We find linkage signals on chromosome 17 consistent with the earlier report of Levy et al. [1].

## Methods

### Subjects, trait definitions, marker data

Phenotype data were available for 2885 subjects from 330 pedigrees (with a total of 4692 members) from the Framingham Heart Study. The data included age, systolic blood pressure (mm Hg), hypertension treatment (yes/no), sex, height, and weight measured at 2- to 4-year intervals. Data were not available on every subject at the same set of ages because of staggered entry, drop-out, and intermittent missing data.

To ensure that we had phenotype data on comparable ages for as many subjects as possible, we averaged systolic blood pressure (SBP) and body mass index (BMI) over any measurements taken during four age intervals: younger than 35 years; between 35 and 50; between 50 and 55; and older than 65. On average, 1.6, 5.3, 6.3, and 5.9 SBP measurements were available on original cohort members in these four age intervals, respectively. The corresponding numbers for the offspring cohort were 1.8, 2.2, 2.2, and 1.6 (lower due to the longer interval between exams). Interval-specific hypertension treatment phenotypes were defined to be "yes" if the subject received any hypertension treatment during the interval and "no" otherwise.

We adjusted SBP for the effect of hypertension treatment using a procedure similar to that outlined by Levy et al. [1] for each age interval separately. The adjusted SBP values for each age interval were then regressed on sex and BMI. We used the residuals from this regression as quantitative trait(s) in (multivariate) variance components analyses described below.

Marker genotype data were available on 1702 subjects. We used data on 16 markers along chromosome 17; the markers were roughly equidistant, one per 10 cM.

### Variance components for multivariate traits

For each subject *j* = 1,...,*J*_{
i
}, in family *i* = 1,...,*I*, let **Y**_{
ij
}be the *p*-dimensional vector of SBP residuals. We take *p* = 1, 2, or 4 for (respectively) the four age intervals considered separately, the two intervals 35–50 and 50–65 considered simultaneously, and all four intervals considered simultaneously.

We use the variance components model **Y**_{
ij
}= **μ** + **X β** + **a**_{
ij
}+ **g**_{
ij
}+ **e**_{
ij
}[8–10]. Here μ and β are fixed effects, and **X** is a matrix of (possibly age-specific) covariates. For these analyses we fit no fixed effects as we regressed on relevant covariates when creating the age-interval-specific trait data. The **a**_{
ij
}, **g**_{
ij
}, and **e**_{
ij
}terms are multivariate normal random effects with mean **0**. Here **a**_{
ij
}is an additive polygenic effect, **g**_{
ij
}is the additive effect due to a specific locus, and **e**_{
ij
}is individual-specific error.

Writing **Y** = (**Y**_{11}' **Y**_{12}' ... **Y**_{
IJi
}')' as the *p* (Σ_{
i
}*J*_{
i
}) × 1 = *p n* × 1 concatenated vector of all subjects' trait vectors,

Cov(**Y**,**Y**) = **K****A** + **Π****G** + **I** **E**,

where K is the *n* × *n* matrix of subjects' kinship coefficients; Π is the *n* × *n* matrix of identity-by-descent (IBD) sharing probabilities for all possible pairs of subjects (calculated at a given location using marker data); **I** is the *n* × *n* identity matrix; **A** is Cov(**a**_{
ij
},**a**_{
ij
}) = {σ_{
akl
}^{2}}, with *k* and *l* indexing age intervals; **G** is Cov(**g**_{
ij
},**g**_{
ij
}) = {σ_{
gkl
}^{2}}; and **E** is Cov(**e**_{
ij
},**e**_{
ij
}) = {σ_{
ekl
}^{2}}. Particular parametric forms for **A**, **G** or **E** – such as autoregressive or exponential decay – can be adopted [6, 7]. In this case, since *p* is small, we used the general form for **E** (which can be parameterized in terms of the variances σ_{
ekk
}^{2} ≥ 0 and correlations -1 ≤ ρ_{
kl
}≤ 1, *k*, *l* = 1,...,*p*).

Polygenic heritabilities can be estimated by fitting the reduced model with **g**_{
ij
}≡ **0**. To test the null hypothesis of complete polygenic pleiotropy, we first fit the constrained model ρ_{
kl
}≡ 1 by writing **A** = Δ Δ', with Δ' = (δ_{1},...,δ_{
p
}). This model is compared to the general model (with the *q* = *p* (*p* - 1)/2 correlation parameters allowed to range between -1 and 1) via a likelihood ratio test. Asymptotically, this test statistic is distributed as the mixture of υ_{
r
}^{2} variables *r* = 0,...,*q*, with mixing probabilities Binom(*q*,*r*) 2^{-q}[11]. The proportion of polygenic variance for two traits due to shared genes can be estimated as ρ_{
kl
}^{2} or (σ_{
akl
}^{2}/σ_{akkσall})^{2} [6].

LOD scores for linkage can be calculated by taking the log_{10} of the likelihood ratio of the model, which estimates **G** to the polygenic model with **g**_{
ij
}≡ **0**. We use the general form for **A** and **G**. In principle, models with varying constraints on these variance matrices could be compared via Aikake's information criterion.

We used SIMWALK2 [12] to calculate IBD sharing and Fisher [13] to calculate maximum likelihood estimates for the variance components parameters. Fisher allows the user to specify that the polygenic (or major gene) variance have the structure **A** = Λ Λ ', with

Here Λ_{1} is a lower triangular *s* × *s* matrix (*s* <*p*) and Λ_{2} is a general (*p* - *s*) × *s* matrix. This amounts to parameterizing **A** in terms of *s* independent factors. Taking *s* = 1 leads to the parameterization **A** = Δ Δ' described above; taking *s* = *p* leads to a Cholesky decomposition of the general form for **A**.

## Results

Both the univariate and multivariate analyses find evidence for a polygenic component to SBP. The estimates of polygenic variance for residual SBP from all models are more than 3 standard errors away from 0 (Tables 1,2,3). Univariate heritabilities range from 27% to 37%. The multivariate analyses show strong correlation between the polygenic components for different age ranges; correlations range from 0.47 to 0.92 (Table 4).

Tests for incomplete pleiotropy are significant, however (test statistics of 8.4 and 21.8 with *p*-values of 2 × 10^{-3} and 2 × 10^{-4} for the two- and four-trait analyses, respectively). This suggests that the set of genes involved in regulating blood pressure differs with age. According to parameter estimates from the unconstrained four-trait analysis (Table 3), the proportion of variance due to shared genes for SBP between age intervals 35–50 and 50–65 is about 82%, but the proportion due to shared genes between age intervals 0–35 and 65+ is 22%.

Univariate LOD scores for the 50–65 age range achieve a peak of 2.2 on chromosome 17 at approximately 80 cM. Bivariate lod scores considering the residuals from the 35–50 and 50–65 age ranges peak just over 2.4 at the same location. This position corresponds roughly to that found by Levy et al. [1] with a lod score of 4.7.

## Discussion

Our results suggest that different sets of genes regulate blood pressure at different ages. Havill and Mahaney [14] find evidence for incomplete pleiotropy in the polygenic component of SBP in the fourth and sixth decades of life, although their estimates of the proportion of polygenic variance due to shared genes are lower than ours for comparable age ranges. These differences may be due to different sample sizes or trait definitions. Furthermore, Havill and Mahaney restrict their analysis to Framingham subjects measured at both age intervals; we did not.

We also find linkage signals on chromosome 17 consistent with the earlier report [1], although with smaller LOD scores. This reduction may be due to smaller sample size or differences in trait definition. Our analyses also do not take other known risk factors for high blood pressure such as smoking or cholesterol into account, although fixed effects for these factors can be included in the variance components model. Similarly, fixed effects for birth cohort or calendar time could be included, perhaps to capture differences in unmeasured risk factors. However, fitting age, birth cohort, and calendar time effects simultaneously could lead to problems with model identifiability [15].

The small increase in multivariate LOD scores over univariate LOD scores likely reflects the fact that when traits are strongly positively correlated, univariate analyses are more powerful [5]. This may also be why the approach of averaging SBP measures taken by Levy et al. [1] performs better than multivariate analysis in this case.

In contrast to our approach, which defines traits based on Framingham subjects' biological age, de Andrade and Olswold [16] define traits based on calendar time (exam number). They fail to replicate the linkage signal on chromosome 17. This is consistent with the results of Mathias et al. [17], who found that age-matched analyses of SBP from the Framingham Heart Study produced higher estimates of heritability than calendar-year-matched analyses. This suggests care should be taken when defining traits in a longitudinal analysis, depending on whether gene × age interaction or gene × calendar time (environment) interaction is more relevant.

A general drawback to the multivariate variance components approach sketched here is that missing observations on subjects with incomplete data are assumed to be missing at random [6] and hence ignorable. However, non-ignorable missingness is a concern in this case, where early withdrawal from study is likely correlated with elevated blood pressure (and hence any genes associated with high blood pressure). Multiple imputation methods can be applied in the case of non-ignorable missing data, as discussed by Allison [18] and several contributions to GAW13 [19, 20].

## References

Levy D, DeStefano AL, Larson MG, O'Donnell CJ, Lifton RP, Gavras H, Cupples LA, Myers RH: Evidence for a gene influencing blood pressure phenotypes in subjects from the Framingham study: genome scan linkage results for longitudinal blood pressure phenotypes in subjects from the Framingham Heart Study. Hypertension. 2000, 36: 477-483.

Pérusse L, Moll PP, Sing CF: Evidence that a single gene with gender- and age-dependent effects influences systolic blood pressure determination in a population-based sample. Am J Hum Genet. 1991, 49: 94-105.

Cheverud JM, Routman EJ, Duarte FAM, van Swinderen B, Cothran K, Perel C: Quantitative trait loci for murine growth. Genetics. 1996, 142: 1305-1319.

de Andrade M, Gueguen R, Visvikis S, Sass C, Siest G, Amos C: Extension of variance components approach to incorporate temporal trends and longitudinal pedigree data analysis. Genet Epidemiol. 2002, 22: 221-232. 10.1002/gepi.01118.

Amos C, de Andrade M, Zhu DK: Comparison of multivariate tests for genetic linkage. Hum Hered. 2001, 51: 133-144. 10.1159/000053334.

Blangero J: Statistical approaches to human adaptability. Hum Biol. 1993, 65: 941-966.

Almasy L, Towne B, Peterson C, Blangero J: Detecting genotype × age interaction. Genet Epidemiol. 2001, 21 (suppl 1): S819-S824.

Almasy L, Dyer TD, Blangero J: Bivariate quantitative trait linkage analysis: pleiotropy vs. 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.

Lange K, Boehnke M: Extensions to pedigree analysis. IV. Covariance components models for multivariate traits. Am J Med Genet. 1983, 14: 513-524. 10.1002/ajmg.1320140315.

Boehnke M, Moll MP: Univariate and bivariate analyses of cholesterol and triglyceride levels in pedigrees. Am J Med Genet. 1986, 23: 775-792. 10.1002/ajmg.1320230306.

Self SG, Liang K: Asymptotic properties of maximum likelihood estimators and likelihood ratio tests under non-standard conditions. J Am Stat Assoc. 1987, 82: 605-610. 10.2307/2289471.

Sobel E, Lange K: Descent graphs in pedigree analysis: applications to haplotyping, location scores, and marker-sharing statistics. Am J Hum Genet. 1996, 58: 1323-1337.

Lange K, Weeks D, Boehnke M: Programs for pedigree analysis: Mendel, Fisher and dGene. Genet Epidemiol. 1988, 5: 471-472. 10.1002/gepi.1370050611.

Havill LM, Mahaney MC: Pleiotropic effects on cardiovascular risk factors within and between fourth and sixth decades of life: implications for genotype-by-age interactions. BMC Genetics. 2003, 4 (suppl 1): S54-10.1186/1471-2156-4-S1-S54.

Kupper LL, Janis JM, Karmous A, Greenberg BG: Statistical age-period-cohort analyses: a review and critique. J Chronic Dis. 1985, 38: 811-830. 10.1016/0021-9681(85)90105-5.

de Andrade M, Olswold C: Comparison of longitudinal variance components and regression based approach for linkage detection on chromosome 17 for systolic blood pressure. BMC Genetics. 2003, 4 (suppl 1): S17-10.1186/1471-2156-4-S1-S17.

Mathias RA, Roy-Gagnon M, Justice CM, Papanicolaou GJ, Fan YT, Pugh EW, Wilson AF: Comparison of year-of-exam and age matched estimates of heritability in the Framingham Heart Study. BMC Genetics. 2003, 4 (suppl 1): S36-10.1186/1471-2156-4-S1-S36.

Allison PD: Missing Data. Thousand Oaks, Sage University Papers Series on Quantitative Applications in the Social Sciences. 2001

Fridley B, Rabe K, de Andrade M: Imputation methods for missing data for polygenic models. BMC Genetics. 2003, 4 (suppl 1): S42-10.1186/1471-2156-4-S1-S42.

Kang T, Kraft P, Gauderman WJ, Thomas DC: Multiple imputation methods for longitudinal blood pressure measurements from the Framingham Heart Study. BMC Genetics. 2003, 4 (suppl 1): S43-10.1186/1471-2156-4-S1-S43.

## Acknowledgements

Part of this work was funded by NIH/NIAD 5 T32 AI07370, Biostatistics Training in AIDS research.

## Author information

### Authors and Affiliations

### Corresponding author

## Rights and permissions

## About this article

### Cite this article

Kraft, P., Bauman, L., Yuan, J.Y. *et al.* Multivariate variance-components analysis of longitudinal blood pressure measurements from the Framingham Heart Study.
*BMC Genet* **4**
(Suppl 1), S55 (2003). https://doi.org/10.1186/1471-2156-4-S1-S55

Published:

DOI: https://doi.org/10.1186/1471-2156-4-S1-S55