Volume 4 Supplement 1

## Genetic Analysis Workshop 13: Analysis of Longitudinal Family Data for Complex Diseases and Related Risk Factors

# Comparison of longitudinal variance components and regression-based approaches for linkage detection on chromosome 17 for systolic blood pressure

- Mariza de Andrade
^{1}Email author and - Curtis Olswold
^{1}

**4(Suppl 1)**:S17

**DOI: **10.1186/1471-2156-4-S1-S17

© Andrade and Olswold; licensee BioMed Central Ltd 2003

**Published: **31 December 2003

## Abstract

We compare two methods to detect genetic linkage by using serial observations of systolic blood pressure in pedigree data from the Framingham Heart Study focusing on chromosome 17. The first method is a variance components (VC) approach that incorporates longitudinal pedigree data, and the second method is a regression-based approach that summarizes all longitudinal measures in one single measure. No evidence of linkage was found either using the VC longitudinal approach or the regression-based approach, except when all time points were used from Cohorts 1 and 2 and only subjects aged 25 and 75 years were included.

## Background

Modeling becomes more complex when observations are recorded over time. Several authors have developed and reviewed statistical methods for longitudinal cohort studies [1]. For genetic analysis, Province and Rao [2, 3] used path analysis to estimate the genetic and environmental effects in families in the presence of temporal trends or a time effect, but did not include variance components (VC) to measure effects from specific genes. Models using structural equations have also been developed and widely applied in the field of behavioral genetics to twin studies [4–6], but these models, primarily directed to the estimation of polygenic and environmental effects, are difficult to use for studies of large families or extended pedigrees. Huggins and colleagues [7] applied cubic spline methods to analyze longitudinal data in twins, but their method considered only additive polygenic effects. Recently, de Andrade et al. [8] proposed an extension of the VC approach to incorporate longitudinal family data. Another approach to analyze longitudinal family data is to use regression methods ignoring the family structure, and then use the residuals as the quantitative traits. Several authors [9, 10] have proposed different variations of this approach. The regression-based approach summarizes all longitudinal measures in one single measure and uses this summarized measure that is the regression residuals, for linkage detection, i.e., each individual has one single value that represents his/her longitudinal measures of blood pressure. On the other hand, the longitudinal VC approach does not summarize the longitudinal measures and uses them all, i.e., each individual has a vector of measures that represents his/her measure of blood pressure levels at each time point. We focus our analyses primarily on chromosome 17 because it was the chromosome where Levy et al. [9] showed strong evidence of linkage for systolic blood pressure. Therefore, the purpose of this paper is to compare these two approaches, the longitudinal VC approach proposed by de Andrade et al. [8] and the residuals approach proposed by Levy et al. [9], using the Framingham Heart Study data set.

## Methods

### Longitudinal VCs

The longitudinal variance components (LVC) approach is an extension of the VC approach proposed by Amos [11]. For longitudinal familial data, let **Y**_{
i
} **= (Y**_{
i1
}**,...,Y**_{
iT
}**)'** be a vector of T time point trait values for k_{i} members of the i^{th} family, where **Y**_{
it
}**' = (Y**_{
i1t
}**, ..., Y**_{
ijt
}**,..., Y**_{
ikit
}**)'** for t = 1,..., T. Let **E(Y**_{
i
}**) = μ + X**_{
i
} **β** and **V**_{
i
} **= A**
**G**_{
i
} **+ B**
**Π**_{
i
} **+ C**
**I**_{
i
}, where
defines the direct product of two matrices; **G**_{
i
} is a k_{i} × k_{i} matrix of coefficient of relationship between pairs of relatives; **Π**_{
i
} is a k_{i} × k_{i} matrix of IBD values for the i^{th} family; **I**_{
i
} is a k_{i} × k_{i} identity matrix; and **A**, **B**, and **C** are, respectively, polygenic, major gene, and random environment variance-covariance matrices each of dimension T × T. These matrices are represented by **A** = (σ_{a.tt'}), σ_{a.tt} = σ^{2}_{a.t}, **B** = (σ_{g.tt'}), σ_{g.tt} = σ^{2}_{g.t}, **C** = (τ_{g.tt'}), τ_{g.tt'} = τ^{2}_{a.t}, with their typical elements in the parentheses. We assume **Y**_{
i
} follows a multivariate normal distribution. More details about this method can be found in de Andrade et al. [8].

To test for genetic linkage, we construct a likelihood ratio test. Under the null hypothesis, the major gene parameters, the σ_{g.tt} for all t and t', are restricted to be equal to **0**. The distribution of the longitudinal test is a mixture of χ^{2} values [12]. For example, for two time points linkage analysis of an additive genetic effect, the distribution of the bivariate test that the major-gene covariance components are zero is a mixture of 1/4 χ^{2}_{0}, 1/2 χ^{2}_{1}, and 1/4 χ^{2}_{3}.

The longitudinal feature was incorporated in the software ACT [13] within the module multic, which was used to run the analyses. The longitudinal multipoint linkage analysis was performed only in the concordant five time points from Cohorts 1 and 2. The trait of interest was systolic blood pressure (SBP), and the covariates were age, gender, and body mass index (BMI). Individuals with missing values were eliminated from the analysis.

### Residuals approach

We used the two-stage procedure described in Levy et al. [9]. This procedure first calculates the within-subject mean BP, and second, uses the sample-wide regressions adjusted for age and BMI, yielding a residual for each subject. Then, these residuals are the traits used in the quantitative linkage analysis. In our analyses, we consider two cases: 1) all subjects in both cohorts regardless of age and 2) only subjects between 25 and 75 years. Each of these analyses was then stratified by time points: 1) using all 21 time points from Cohort 1 and all five time points from Cohort 2 and 2) using only the concordant five time points from Cohorts 1 and 2. After the residuals were obtained, multipoint quantitative linkage analyses were performed for these four different scenarios using SOLAR [14]. We also performed a multipoint quantitative linkage analysis using the average SBP over all measurements for each subject for additional comparison. LOD score values were calculated by dividing the likelihood ratio statistic by 4.6 (2/log e).

## Results

## Discussion

Although evidence of a gene influencing blood pressure was reported on chromosome 17 using the Framingham Heart Study [9], the findings (LOD = 4.7, position = 67 cM) could not be replicated in our analyses. Only when SBP from all time points in Cohorts 1 and 2 were used and the subjects' age were restricted to be between 25 and 75 years was a comparable finding reached (LOD = 3.0, position = 68 cM). Slager and Iturria [10] also could not replicate Levy et al.'s findings when using two other approaches of phenotype definitions. We also expected that the longitudinal VC approach would perform much better because it takes into account a temporal trend affecting the genetic variability of the SBP. However, the maximum LOD found using this procedure was around 0.96 in the surrounding region of the 114 cM position for time points 1 and 4 (Figure 1C). There was no consistency among the pair-wise time points multipoint LOD values.

The reason why the results of the two approaches differ could be several. It could be due to missing values, since individuals with missing values were eliminated from the analysis in the longitudinal VC but not from the regression-based approach. The amount of missing SBP values among pair-wise time points is around 50%, where time point 5 has more missing values than the other time points. Another explanation could be the adjustment that was used for treatment effect in SBP in Levy et al. [9], i.e., the inferred values could be larger than expected, thus inflating the LOD values [15]. In our longitudinal VC analysis and in our regression-based analysis we did not adjust for treatment effect. We assume that if the treatment works it lowers the SBP. Therefore, the impact of not correcting for treatment makes the linkage results more conservative, which is preferable to false-positive results. Finally, the departure of normality could also inflate the LOD values, i.e., increase the false-positive rate [16].

In summary, the longitudinal VC and the regression approach showed similar results when all subjects were used regardless of their age, i.e., no evidence for linkage was found on chromosome 17. However, the evidence for linkage was only observed when all SBP time points were analyzed from both cohorts and when subjects' age were restricted to be between 25 and 75. Therefore, one possible conclusion is that the linkage results depend heavily on the phenotype. By using the residuals as a summarized measure of the longitudinal SPB, the time point variability and missing mechanism are not taken into account and could therefore lead to spurious linkage results.

## Declarations

### Acknowledgments

The authors thank Brooke Fridley, Beth Atkinson, and Kari Rabe for their help and two anonymous reviewers for their helpful comments.

## Authors’ Affiliations

## References

- Diggle PJ, Liang KY, Zeger SL: Analysis of Longitudinal Data. Oxford, Clarendon Press. 1995Google Scholar
- Province MA, Rao DC: A new model for the resolution of cultural and biological inheritance in the presence of temporal trends: application to systolic blood pressure. Genet Epidemiol. 1985, 2: 363-374. 10.1002/gepi.1370020405.View ArticlePubMedGoogle Scholar
- Province MA, Rao DC: Familial aggregation in the presence of temporal trends. Stat Med. 1988, 7: 185-198.View ArticlePubMedGoogle Scholar
- Eaves LJ, Long J, Heath AC: A theory of developmental change in quantitative phenotypes applied to cognitive development. Behav Genet. 1986, 16: 143-161. 10.1007/BF01065484.View ArticlePubMedGoogle Scholar
- Boomsma DI, Molenaar PCM: The genetic analysis of repeated measures. I. Simplex models. Behav Genet. 1987, 17: 111-123. 10.1007/BF01065991.View ArticlePubMedGoogle Scholar
- Neale MC, Cardon LR: Methodology for Genetic Studies of Twins and Families. Dordrecht, Kluwer Academic Publishers. 1992View ArticleGoogle Scholar
- Huggins RM, Hoang NH, Loesch NH: Analysis of longitudinal data in twins. Genet Epidemiol. 2000, 19: 345-353. 10.1002/1098-2272(200012)19:4<345::AID-GEPI6>3.0.CO;2-S.View ArticlePubMedGoogle Scholar
- de Andrade M, Gueguen R, Visvikis S, Sass C, Siest G, Amos CI: Extension of variance components approach to incorporate temporal trends and longitudinal pedigree data analysis. Genet Epidemiol. 2002, 22: 221-232. 10.1002/gepi.01118.View ArticlePubMedGoogle 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, 36: 477-483.View ArticlePubMedGoogle Scholar
- Slager SL, Iturria SJ: Genome-wide linkage analysis of systolic blood pressure: a comparison of two approaches to phenotype definition. BMC Genetics. 2003, 4 (suppl 1): S13-10.1186/1471-2156-4-S1-S13.PubMed CentralView ArticlePubMedGoogle Scholar
- Amos CI: Robust variance-components approaches for assessing genetic linkage in pedigrees. Am J Hum Genet. 1994, 54: 535-543.PubMed CentralPubMedGoogle 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
- de Andrade M, Krushkal J, Yu L, D Zhu D, Amos CI: ACT – A computer package for analysis of complex traits. Am J Hum Genet. 1998, 63: 287-Google Scholar
- Almasy L, Blangero J: Multipoint quantitative-trait linkage analysis in general pedigrees. Am J Hum Genet. 1998, 62: 1198-1211. 10.1086/301844.PubMed CentralView ArticlePubMedGoogle Scholar
- de Andrade M, Fridley B, Boerwinkle E, Turner ST: Diagnostic tools in linkage analysis for quantitative traits. Genet Epidemiol. 2003, 24: 302-308. 10.1002/gepi.10236.View ArticlePubMedGoogle Scholar
- Allison DB, Neale MC, Zannolli R, Schork NJ, Amos CI, Blangero J: Testing the robustness of the likelihood ratio test in a variance-component quantitative trait loci (QTL) mapping procedure. Am J Hum Genet. 1999, 65: 531-545. 10.1086/302487.PubMed CentralView ArticlePubMedGoogle Scholar

## Copyright

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.