Prediction of hypertension based on the genetic analysis of longitudinal phenotypes: a comparison of different modeling approaches for the binary trait of hypertension
- Yun-Hee Choi^{1}Email author,
- Rafiqul Chowdhury^{1} and
- Balakumar Swaminathan^{1}
https://doi.org/10.1186/1753-6561-8-S1-S78
© Choi et al.; licensee BioMed Central Ltd. 2014
Published: 17 June 2014
Abstract
For the analysis of the longitudinal hypertension family data, we focused on modeling binary traits of hypertension measured repeatedly over time. Our primary objective is to examine predictive abilities of longitudinal models for genetic associations. We first identified single-nucleotide polymorphisms (SNPs) associated with any occurrence of hypertension over the study period to set up covariates for the longitudinal analysis. Then, we proceeded to the longitudinal analysis of the repeated measures of binary hypertension with covariates including SNPs by accounting for correlations arising from repeated outcomes and among family members.
We examined two popular models for longitudinal binary outcomes: (a) a marginal model based on the generalized estimating equations, and (b) a conditional model based on the logistic random effect model. The effects of risk factors associated with repeated hypertensions were compared for these two models and their prediction abilities were assessed with and without genetic information.
Based on both approaches, we found a significant interaction effect between age and gender where males were at higher risk of hypertension before age 35 years, but after age 35 years, women were at higher risk. Moreover, the SNPs were significantly associated with hypertension after adjusting for age, gender, and smoking status. The SNPs contributed more to predict hypertension in the marginal model than in the conditional model. There was substantial correlation among repeated measures of hypertension, implying that hypertension was considerably correlated with previous experience of hypertension. The conditional model performed better for predicting the future hypertension status of individuals.
Keywords
Background
Hypertension is a chronic condition resulting from high blood pressure in the arteries during circulation. Clinically, a person is said to be hypertensive if the individual's systolic blood pressure (SBP) is greater than 140 mm Hg or diastolic blood pressure (DBP) is greater than 90 mm Hg. With advances in genome-wide association studies, several researchers have also investigated the role of genes in this disease [1, 2]. It is essential to control hypertension in order to avoid consequences like cardiovascular diseases, stroke, and heart and kidney failure.
The San Antonio Family Study data for Genetic Analysis Workshop 18 (GAW18) contain up to 4 longitudinal measures of SBP and DBP along with their background and genetic information for a total of 932 individuals from 20 Mexican American families. In the analysis of the longitudinal hypertension family data, we focus on modeling longitudinal binary traits of hypertension, defined by SBP >140 mm Hg, DBP >90 mm Hg or use of antihypertensive medication, by taking into account correlations arising from repeated outcomes and among family members while controlling for covariates such as age, gender, smoking status, and genetic polymorphisms.
Our primary objective is to examine predictive abilities of longitudinal models with inclusion of genetic information. In the first step, we identify important single-nucleotide polymorphisms (SNPs) associated with any occurrence of hypertension over the study period in order to set up covariates for the longitudinal analysis. The selection of SNPs is based on chromosome 3 only. Then we proceed to the longitudinal analysis of repeated measures of hypertension with covariates, including SNPs that already identified. We examined two popular modeling approaches for longitudinal binary outcomes: the marginal model (population-average) and the conditional model (subject-specific). The effects of the risk factors associated with repeated hypertension from the two models were compared and their prediction abilities were assessed with and without genetic information using the areas under the receiver operating characteristic curve.
Methods
Selection of associated SNPs
where λ_{ij}(t) is the hazard function of individual i in family j, λ_{0}(t) is the baseline hazard function, b_{j} is the random effect for family j, and X_{ij} is the vector of covariates for the fixed effects. Each SNP entered into the model was assumed to be codominant (TT, TC, CC) as the inheritance mode is not known. Similarly, the logistic model was fitted for each SNP for the binary hypertension outcome, defined as 1 if one has ever experienced hypertension over the study period and 0 otherwise, with a random intercept to account for familial correlation. These two modeling approaches were employed to validate the choice of SNPs for the following longitudinal models.
Longitudinal models
We considered the two most commonly used modeling strategies for analyzing correlated data from longitudinal study: marginal and conditional models, also known as population-averaged and subject-specific approaches, respectively. To accommodate the dependence among longitudinal repeated outcomes, the marginal approach directly describes marginal means along with a prespecified working correlation structure as nuisance parameters, whereas the conditional approach describes an individual response conditioning on the unobserved values of random effects. The generalized estimating equations (GEEs) method by Liang and Zeger [5] was used to estimate population-average effects in the marginal model, whereas the logistic random effect model [6] was used to provide subject-specific effects in the conditional model.
Marginal model (population-averaged model)
where the risk factors of interest are age at examination, sex (1 = male, 2 = female), smoking status (1 = smoker, 0 = nonsmoker), and a vector of SNPs identified in the first step. The interaction between age and sex was also considered in the model.
The analysis of correlated traits was based on GEEs [5], using up to 3 visits per subject.
Conditional model (subject-specific model)
where ν_{jk}^{(2)} represents a random effect for individual j within family k following iid N(0, σ_{(2)}^{2}) and ν_{k}^{(3)} is a random effect for family k with iid N(0, σ_{(3)}^{2}). The two random intercepts describe the dependence at two levels: one at the individual level among repeated outcomes and another at the family level among family members.
Results
Our analysis focused on a binary trait of hypertension defined by SBP >140 mm Hg, DBP >90 mm Hg, or use of antihypertensive medication. A total of 932 participants from 20 families were considered. The genome-wide association analyses were performed for 65,519 SNPs on chromosome 3 using the software PLINK [7] with R-plugin, and statistical analysis was conducted using SAS [8] and R [9].
Data preprocessing
After the data were preprocessed for quality control, 850 individuals were kept with no missing phenotypes and 6261 SNPs were excluded with Hardy-Weinberg equilibrium test p value (among founders only) <0.001, or minor allele frequency <0.01, or with a missing genotyping rate >0.95 from our analyses.
SNP selection
Top 5 most significant SNPs associated with time to hypertension based on Cox PH model with frailty
Chr | SNP | Basepair position | pValue | In/near gene (within 60 kb) |
---|---|---|---|---|
3 | rs10510257 | 3346138 | 1.09080 × 10^{−5} | |
3 | rs1047115 | 186358366 | 1.38323 × 10^{−5} | FETUB |
3 | rs5024851 | 247473 | 1.67662 × 10^{−5} | CHL1 |
3 | rs7630698 | 189199930 | 2.19678 × 10^{−5} | |
3 | rs704903 | 43070847 | 4.09467 × 10^{−5} | FAM198A |
Top 5 most significant SNPs associated with any occurrence of hypertension over the study period based on logistic random effect model
Chr | SNP | Basepair position | pvalue | In/near gene (within 60 kb) |
---|---|---|---|---|
3 | rs10510257 | 3346138 | 6.95197 × 10^{−6} | |
3 | rs1047115 | 186358366 | 1.73097 × 10^{−5} | FETUB |
3 | rs719318 | 10474137 | 1.89300 × 10^{−5} | ATP2B2 |
3 | rs6807497 | 67015910 | 2.31716 × 10^{−5} | |
3 | rs1456217 | 66959472 | 3.16978 × 10^{−5} |
Marginal and conditional models for longitudinal data
Comparison of marginal and conditional models: estimated coefficients from the two models for the longitudinal binary hypertension traits
Marginal model | Conditional model | |||||
---|---|---|---|---|---|---|
Variables | Log OR | SE | pvalue | Log OR | SE | pvalue |
Intercept | −4.002 | 1.281 | 0.0018 | −5.012 | 1.618 | - |
Age (years) | 0.053 | 0.026 | 0.0429 | 0.072 | 0.032 | 0.0274 |
Gender | −1.209 | 0.811 | 0.1359 | −1.206 | 0.964 | 0.2112 |
Smoke | 0.203 | 0.253 | 0.4211 | 0.213 | 0.315 | 0.5001 |
Age × gender | 0.036 | 0.017 | 0.0331 | 0.036 | 0.019 | 0.0621 |
Rs10510257(AA) | −1.097 | 0.574 | 0.0558 | −1.399 | 0.747 | 0.0615 |
Rs10510257(AG) | −0.887 | 0.286 | 0.0019 | −1.119 | 0.313 | 0.0004 |
Rs1047115 (GT) | 0.712 | 0.431 | 0.0985 | 0.903 | 0.536 | 0.0925 |
Random effects | ||||||
σ_{(2)} for ID | 1.388 | |||||
σ_{(3)} for PEDNUM | 1.092 |
We observed that the estimates of the regression coefficients in the conditional analysis were slightly farther away from zero than those from the marginal model. For example, the log odds ratio (OR) of hypertension between smokers and nonsmokers was 0.213 in the conditional model compared to 0.203 in the marginal model. The log OR of rs10510257 genotype AA compared to genotype GG was −1.399 for the conditional model and −1.097 for the marginal model. In addition, the conditional model allowed us to measure intraclass correlations via the variances of random effects. The estimated variance of individual-level random effect was 1.388 and that of family-level random effect was 1.092. They yielded the estimated intraclass correlation across repeated measures of hypertension in the same individual equal to 0.43 and familial correlation equal to 0.19, which implies that hypertension is substantially correlated with previous experience of hypertension.
Prediction ability
Discussions
In the analysis of longitudinal family data for genetic association studies, the correlation structure of the data is often complex and requires adequate modeling. In addition, population-average effect estimates are smaller in magnitude than subject-specific effect estimates with the difference increasing with intraclass correlation, as our analyses revealed [11].
It is worth mentioning that in our conditional model, multiple random effects were included to address the nested structure of the data where individuals are repeatedly measured over time and nested within families; in particular, we assumed a common random effect shared within families to explain familial correlation induced by shared latent environmental and genetic risk factors. However, using a single random effect may not be sufficient if we want to describe more complex correlation structure within families. The relativeness within families can be explicitly modeled using multivariate random variables or using kinship coefficient.
In our modeling of longitudinal data, although the SNPs identified were significant risk factors in the model, the improvement of predicting capability of the model measured by the AUCs appeared to be very minor yet not negligible. Notice that our genome scan was only on chromosome 3. As one reviewer pointed out, a comprehensive predictive model would not be possible without whole-genome exploration.
Conclusion
Our study demonstrated that the genetic information plays an important role in predicting future hypertension event. In our modeling of longitudinal data, the SNPs identified were significant risk factors in the model whereas the improvement of predicting capability of the model measured by the AUCs appeared to be very minor yet not negligible. We found that the conditional modeling approach could make maximum use of the information provided by repeated phenotypes, which could lead to better prediction. Therefore, the longitudinal modeling approaches also can be helpful for identifying new genes and developing new treatments for repeated outcomes.
Declarations
Acknowledgements
This research was supported by Natural Sciences and Engineering Research Council of Canada. The GAW18 whole genome sequence data were provided by the T2D-GENES Consortium, which is supported by NIH grants U01 DK085524, U01 DK085584, U01 DK085501, U01 DK085526, and U01 DK085545. The other genetic and phenotypic data for GAW18 were provided by the San Antonio Family Heart Study and San Antonio Family Diabetes/Gallbladder Study, which are supported by NIH grants P01 HL045222, R01 DK047482, and R01 DK053889. The Genetic Analysis Workshop is supported by NIH grant R01 GM031575.
This article has been published as part of BMC Proceedings Volume 8 Supplement 1, 2014: Genetic Analysis Workshop 18. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcproc/supplements/8/S1. Publication charges for this supplement were funded by the Texas Biomedical Research Institute.
Authors’ Affiliations
References
- Ehret GB: Genome-wide association studies: contribution of genomics to understanding blood pressure and essential hypertension. Curr Hypertens Rep. 2010, 12: 17-25. 10.1007/s11906-009-0086-6.PubMed CentralView ArticlePubMedGoogle Scholar
- Barlassina C, Lanzani C, Manunta P, Bianchi G: Genetics of essential hypertension: from families to genes. J Am Soc Nephrol. 2002, 13: S155-S164. 10.1097/01.ASN.0000032524.13069.88.View ArticlePubMedGoogle Scholar
- Cox DR: Regression models and life tables (with discussion). J R Stat Soc Series B Stat Methodol. 1972, 34: 187-220.Google Scholar
- Hougaard P: Analysis of Multivariate Survival Data. New York, Springer. 2000Google Scholar
- Liang KY, Zeger SL: Longitudinal data analysis using generalized linear models. Biometrika. 1986, 73: 13-22. 10.1093/biomet/73.1.13.View ArticleGoogle Scholar
- Wong GY, Mason WM: The hierarchical logistic regression model for multilevel analysis. J Am Stat Assoc. 1985, 80: 513-524. 10.1080/01621459.1985.10478148.View ArticleGoogle Scholar
- Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, Maller J, Sklar P, de Bakker PIW, Daly MJ, et al: PLINK: a toolset for whole-genome association and population-based linkage analysis. Am J Hum Genet. 2007, 81: 559-575. 10.1086/519795.PubMed CentralView ArticlePubMedGoogle Scholar
- SAS Institute: Base SAS 9.2 Procedure Guide: Statistical Procedures. Cary, NC, SAS Institute. 2010Google Scholar
- R Development Core Team: R: A Language and Environment for Statistical Computing. 2011, Vienna, Austria, R Foundation for Statistical ComputingGoogle Scholar
- Pencina MJ, D'Agostino RB, D'Agostino RB, Vasan RS: Evaluating the added predictive ability of a new marker: from area under the ROC to reclassification and beyond. Stat Med. 2008, 27: 157-172. 10.1002/sim.2929.View ArticlePubMedGoogle Scholar
- Neuhaus JM: Estimation efficiency and tests of covariate effects with clustered binary data. Biometrics. 1993, 49: 989-996. 10.2307/2532241.View 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. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.