Multilevel analysis of systolic blood pressure and ACE gene I/D polymorphism in 438 Swedish families – a public health perspective
© Merlo et al; licensee BioMed Central Ltd. 2006
Received: 12 September 2005
Accepted: 01 March 2006
Published: 01 March 2006
Individuals belonging to the same family share a number of genetic as well as environmental circumstances that may condition a common SBP level. Among the genetic factors, the angiotensin converting enzyme (ACE) gene I/D polymorphism appears as a possible candidate as it might influence both SBP and the pharmacological effect of ACE inhibitors. We aimed to combine genetic epidemiology with public health ideas concerning life-course and multilevel epidemiology in order to understand the role of familial factors regarding individual SBP.
We applied multilevel regression analysis on 1926 individuals nested within 438 families from South Sweden. Modelling familial SBP variance as a function of age and use of ACE inhibitors we calculates a variance partition coefficient and the proportional change in familial SBP variance attributable to differences in ACE gene I/D polymorphism
Our results suggest the existence of genetic or environmental circumstances that produce a considerable familial clustering of SBP, especially among individuals using ACE-inhibitors. However, ACE gene I/D polymorphism seems to play a minor role in this context. In addition, familial factors – genetic, environmental or their interaction – shape SBP among non-users of ACE inhibitors but their effect is expressed later in the life-course.
Strategies directed to prevent hypertension should be launched in younger rather than in older ages and both prevention of hypertension and its treatment with ACE inhibitors should be focused on families rather than on individuals.
Individuals belonging to the same family share a number of genetic as well as environmental factors like common physical environment, socioeconomic position, dietary habits, familial attitudes and health believes, etc. Such circumstances may condition a common level of systolic blood pressure (SBP) that may express itself as a clustering of SBP level within families.
It is known that SBP is a central risk factor for cardiovascular diseases . SBP level increases with age and also is conditioned by genetic factors . Moreover, in presence of established hypertension pharmacological treatment of blood pressure is as well a strong determinant of SBP level for obvious reasons.
One of the main mechanisms for blood pressure homeostasis is the renin angiotensin aldosterone system (RAAS) and a common effective group of antihypertensive drugs are angiotensin converting enzyme (ACE) inhibitors. By blocking the ACE these drugs reduce the levels of angiotensin II and thereby decrease vasoconstriction and renal sodium reabsorption which, in turn, leads to lower blood pressure . Among the genetic factors behind a possible familial clustering of SBP, the angiotensin converting enzyme (ACE) gene I/D polymorphism appears as a possible candidate. ACE levels are highly inherited and show a great interindividual variation . It has been shown that the ACE gene (chromosome 17q23) harbours an insertion/deletion (I/D) polymorphism that accounts for 47% of the variance in plasma ACE concentration . However, it is possible that the I/D polymorphism per se may not be responsible for the ACE levels  but rather an unknown nearby polymorphism that is in strong linkage disequilibrium with the ACE gene I/D locus . It is thus plausible that the pharmacological effect of ACE inhibitors drugs could be modified by familial factors, particularly by familial differences in the prevalence of the ACE gene I/D polymorphism. It has also been suggested that ACE gene polymorphism may modulate the physiological age related changes in blood pressure [8, 9].
Therefore, we expected to find a familial clustering for both ACE I/D polymorphism and for SBP level. If this were true, it might also be possible that familial differences in ACE I/D polymorphism could explain some of the familial differences in SBP and thus clustering of SBP. This possible phenomenon might differ by age and use of ACE inhibitors.
Considering all those factors (i.e. family environment, blood pressure, individual age, pharmacological treatment with ACE inhibitors and the ACE gene I/D polymorphism) simultaneously is a methodological challenge. In the present study we aim to apply a methodological approach based on multilevel regression analysis [10, 11] that may be useful in this context [12, 13]. By doing so, we aimed to combine genetic epidemiology with public health ideas concerning life-course and multilevel epidemiology in order to understand the role of familial factors regarding individual SBP.
In this study we pooled three population based family materials from South Sweden consisting of at least two siblings and, when available, their relatives and spouses. The information was collected during 1996–2000 in all the three samples. Two of them (Malmö-B and Skara) included siblings with established hypertension, whereas one (Malmö A) included siblings with a positive family history of established hypertension but presence of established hypertension was not a necessary inclusion criterion. In Malmö-A, probands were identified in the Malmö Diet and Cancer  and Malmö Preventive Project  databases. In Malmö B, probands were identified at five different primary health care centres in southern Sweden. In the Skara sample, probands and their families were ascertained from patients diagnosed with hypertension from two earlier population based surveys in Skara, one investigation of patients with hypertension and diabetes in primary health care , and one investigation of a random sample of Skara population .
The local ethics committee at Gothenburg and Lund University approved the study protocols and all participants gave their informed consent.
Assessment of blood pressure
All subjects were investigated in the morning after at least 10-hours fast. Blood pressure was measured in the right brachial artery with subjects in the supine position after at least 5–15 minutes rest by specially trained nurses using a using Tricuff® sphygmomanometer for adjustment of cuff size to arm circumference . The systolic and diastolic blood pressures were recorded to the nearest 2 mmHg and the mean of three measurements was recorded. Established hypertension was defined as having at least three consecutive blood pressure recordings at different occasions ≥ 160 mmHg systolic and/or ≥ 90 mmHg diastolic or being on medication for treatment of blood pressure.
High molecular DNA was extracted from blood leucocytes using standard techniques . The ACE gene I/D polymorphism was genotyped according to the methods described by Rigat et al  except that 1.5 mmol MgCl2 and 1% formamid were used in the PCR reaction. To exclude mistyping of the heterozygotes as DD homozygotes all the DD genotypes samples were confirmed with an insertion specific PCR .
Use of medication
In order to assess current pharmacological treatment the subjects were asked to bring receipts and pill boxes and the information was asserted by a face to face interview performed by specially trained nurses. ACE inhibitors were defined according the Anatomical Therapeutic Chemical classification system code C09A .
Statistical analysis: multilevel regression analysis
Individual SBP was analyzed by multilevel linear regressions models [10, 23] with individuals at the first level and families (i.e., pedigrees) at the second level. Multilevel regression analysis is a very appropriate technique for performing variance analysis in unbalanced clusters as families .
We fitted three models. Model I did not include any explanatory variables and only focused on decomposing the total SBP variance (VTotal) in its individual (VI) and familial (VF) components. In model II we included dummy variables for the original study samples (Skara, Malmö-A and B) individual age, and for use of ACE inhibitors. In the last model III we included the ACE gene I/D polymorphism in addition to the variables already included in model II. For the statistical analysis we constructed dummy variables for the presence or absence of homozygosis II, homozygosis DD, and heterozygosis ID, and used homozygosis II as reference in the comparisons. In model II and III we investigate family-individual interaction by allowing the slope of the association between on the one hand SBP and on the other age or ACE inhibitors to be random at the family level. That is, we investigate whether the family level modified those individual level associations.
We calculated the Variance Partition Coefficient (VPC). The VPC is a measure of the extent to which members of a family resemble each other more than they resemble individuals from other families in relation to SBP level. The VPC is the percentage of the total variance (VF + VI) in SBP that is attributed to the family level (VF) and is, therefore, a measure of clustering. In the simplest form (model i) the VPC corresponds to the intraclass correlation coefficient, being a measure of general clustering of individual SBP in the families.
In models II and III, because of the existence of random slopes at the familial level the variance of the slope of the association between SBP and age as well as the variance of the slope of the association between SBP and use of ACE inhibitors are also a coefficient in a function that describes how the between-families variation in SBP differs in individuals with different age and use of ACE inhibitors. This variance function includes also the intercept variance and the covariances. The output of this variance function was obtained from the software used for the multilevel regression analysis. [11, 24] By means of this variance function we calculated the VPC as a function of age and use of ACE inhibitors which provided information on the relevance of familial factors for understanding SBP differences in specific groups of individuals (i.e., defined by age and use of ACE inhibitors). A large VPC would indicate that familial factors are important in order to understand individual SBP differences. On the other hand, a VPC close to zero would indicate that the families exert only a small influence on individual SBP variance.
Using the VF from the models before (Model II) and after (Model III) adjustment for ACE I/D gene polymorphism we calculated the absolute [VF before - VF after] and the relative [((VF before - VF after)/VF before) × 100] change in variance.
The precision of the estimates was appraised by their SE. Parameters were tested in an approximate manner using a normal test (i.e., the ratio of an estimated variance to its standard error) and Wald statistic with a chi-square distribution , and calculated exact p-values. A p-value greater than 0.05 was considered non-significant (NS).
In addition to familial clustering (i.e., VPC) of SBP we also calculated the familial clustering of homozygosis II, homozygosis DD, and heterozygosis ID, as well as the familial clustering of ACE inhibitors use. Since these variables are binary we applied multilevel logistic regression analysis  and calculated the VPC using the latent variable threshold method recommended by Snijders .
The software MLwiN, Version 2  was used to perform the analyses. Regression parameters were estimated using Markov Chain Monte Carlo (MCMC). We use the Deviance Information Criterion (DIC) to compare the goodness of the fit of consecutive models. 
Characteristics of the 1926 individuals from 438 Swedish families with history of hypertension
Malmö Sample A
Number of people
Number of families
Age (mean years)
Systolic blood pressure (mean mmHg)
Use of ACE-inhibitors (%)
ACE gene I/D polymorphism
Family size in the study sample
Number of individuals
Number of families
Multilevel regression analysis of systolic blood pressure (SBP) in 1926 individuals from 438 Swedish families with hypertension
Model I (empty)
Regression coefficient (SE)
Sex (men vs. women)
ACE inhibitors use (yes vs. no) Cohort
ACE gene I/D polymorphism
ACE inhibitors-SBP slope
See figure 2
See figure 2
Goddness of the fit
The results of the complementary logistic regression analyses showed a strong familial genetic clustering for the different forms of the ACE gene I/D polymorphism, specially for homozygous forms DD (VPC = 43%) and II (VPC = 37%) and in a lower degree for the heterozygous form ID (VPC = 15%). Use if ACE inhibitors was also clustered within families (VPC = 9%).
The multilevel regression analysis allowed us to study the distribution of blood pressure in families along individual age and use of ACE inhibitors and simultaneously investigate the role of ACE gene I/D polymorphism in this context.
We observed that SBP level clustered within families to an appreciable degree, and that this clustering increased considerably with age in people that did not use ACE-inhibitors. This fact suggests the existence of specific familial genetic or environmental factors influencing individual SBP, and that the effect of these familial factors express later in the life course. In fact, familial differences in SBP were large in the elderly but very small in youth people.
It is possible that both familial genetic and environmental circumstances (e.g., common physical environment, socioeconomic position, dietary habits, familial attitudes and health believes, etc) affecting SBP act only after latency period. It is also probable that genetic susceptibility affects SBP level only when the effects of familial environmental factors become fully developed through the individual's life course. As a hypothetical example, learned familial behaviours may lead to obesity in adulthood, which in turn increases the risk of hypertension later in life in people with certain genetic predisposition. It is possible that these kinds of mechanisms could explain our results. However, the database used in this investigation was not primarily design to answer those question and information on socioeconomic and behavioural factors was not available. Moreover our study was cross-sectional and we cannot exclude – even if we dot believe it probable – the existence of a cohort effect confounding the age related familial differences in SBP.
In the present study we cannot differentiate between genetic effects (apart from the minor effect of the ACE I/D polymorphism) and effects of familial environment on SBP. However, it is reasonable to believe that, to a large extent, genetic and environmental factors determining SBP interact with each other throughout life. Thus, genetic susceptibility to high SBP is likely to be modified by reducing environmental risk factor load. Therefore, our results suggest that acting on familial environmental factors early in the life course, rather than in late adulthood, might be a more effective strategy to prevent high blood pressure in the population.  Future research will identify which those factors and how we can influence them.
The degree of familial clustering of SBP level was much higher in users of ACE-inhibitors than in non-users (Figure 2B) and this clustering was also considerably less age dependent. This observation indicates that ACE inhibitors produce a more intense blood pressure lowering effect in some families than in others. One reason for this familial heterogeneity might be genetic differences in the susceptibility to the effect of ACE inhibitor, and our hypothesis was that the ACE I/D polymorphism could be a candidate for explaining these familial differences. We also observed that – as expected – the different forms of the ACE I/D polymorphism presented a high degree of clustering within families. However, accounting for differences in the ACE I/D polymorphism explained only a minor part of the familial differences in SBP level among user of ACE inhibitors.
The role of the ACE I/D polymorphism in blood pressure regulation has been extensively studied. Although there is no doubt that the polymorphism is associated with plasma ACE activity, there is less evidence on the role of ACE I/D polymorphism in blood pressure regulation. [5, 27] Blood pressure is a complex trait that results from the effect of multiple genes and many environmental influences . There is a growing agreement from many studies in that single gene effects on antihypertensive drug responses are small, and even the combined effects of all presently known polymorphisms do not account for enough variation in response to be clinically useful. With this in mind, it is not surprising that the ACE I/D polymorphism only explained a very small part of the total familial variance in SBP level. Nevertheless our results suggest that the ACE I/D polymorphism might have some role in blood pressure control among users of ACE inhibitors.
ACE-inhibition was not a first-line antihypertensive therapy in Sweden at the time of patient-recruitment in the present study, indicating that these patients had responded less well on conventional treatment (thiazides and beta-blockers) or that they had concomitant diseases such as diabetes with microalbuminuria or heart failure. One possible explanation for the greater familial clustering of SBP in users than in non-users of ACE inhibitors could thus be that patients treated with ACE-inhibitors had more atherosclerosis and endothelial dysfunction than non-users of ACE inhibitors. In analogy with the effect of age on familial SBP clustering in non-users, it can be speculated that genetic and familial environmental factors affecting SBP have a greater impact when acting on the background of a relatively more pathological or "biologically aged" vascular system.
Most patients treated with ACE inhibitors were also on treatment with other antihypertensive agents, mostly beta-receptor blockers (i.e., 32%) and calcium channel blockers (i.e., 27%) introducing some uncertainty as to how much of the on-treatment blood pressure level was influenced by overall drug effects as opposed to ACE-inhibitor specific effects. However, adjustment for use of those medications in model III had only a minor effect on the slope variance, which was reduced from 131 to 128 (not show in tables).
We advocate a study design that may be useful for studying the influence of genetic and environmental familial factors in public health. This approach may facilitate the investigation of common processes that with multifaceted expression and development across time in different contexts. Multilevel regression analysis  emerges as an appropriate tool for investigating complex patterns of variance, an aspect that is especially useful in genetic epidemiology where the information has a natural hierarchical structure (e.g., pedigrees, generations, individuals) [12, 13].
Our study suggest that familial circumstances (genetic, environmental or their interaction) condition individual SBP level and seems to modify the effect of ACE inhibitors, a fact that deserves more investigation in order to improve the effectiveness of pharmacological treatment of hypertension. Familial circumstances also seem to play a relevant role for understanding SBP level in non-users of ACE inhibitors but those effects express later in the life course. Those facts suggests that strategies directed to prevent hypertension should be launched in younger rather than in older ages and both prevention of hypertension and its treatment with ACE inhibitors should be focused on families rather than on individuals.
This study was supported by grants from the Swedish Research Council (# 2004-6155, #2002-A0785), the Swedish Council for Working Life and Social Research (# 2002-054), Swedish Heart and Lung Foundation (#20020681), the Medical Faculty of Lund University, Malmö University Hospital, the Albert Påhlsson Research Foundation, the Crafoord foundation, the Ernhold Lundströms Research Foundation and the Region Skåne (Hälso- och sjukvårdsledning), Skaraborg Institute and the NEPI foundation. We are most grateful to Mrs Ann Carlsson, Ms Agneta Edman, Mrs Anita Mehner. Philippe Burri and Ms Lena Rosberg for excellent technical assistance.
- Mancia G, Seravalle G, Grassi G: Systolic blood pressure: an underestimated cardiovascular risk factor. J Hypertens. 2002, 20 (Suppl 5): S21-27.Google Scholar
- Turner ST, Boerwinkle E: Genetics of blood pressure, hypertensive complications, and antihypertensive drug responses. Pharmacogenomics. 2003, 4: 53-65. 10.1517/phgs.184.108.40.20687.View ArticlePubMedGoogle Scholar
- Kureya P, Beevers DG: The mechanisms of hypertension and the role of ACE inhibitors. J Hum Hypertens. 1991, 5 (Suppl 2): 3-8.PubMedGoogle Scholar
- Alhenc-Gelas F, Richard JL, Courbon D, Warnet JM, Corvol P: Distribution of plasma angiotensin I-converting enzyme levels in healthy men; relationship to environmental and hormonal parameters. J Lab Clin Med. 1991, 117: 33-39.PubMedGoogle Scholar
- Rigat B, Hubert C, Alhenc-Gelas F, Cambien F, Corvol P, Soubrier F: An insertion/deletion polymorphism in the angiotensin I-converting enzyme gene accounting for half the variance of serum enzyme levels. J Clin Invest. 1990, 86: 1343-1346.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhu X, Bouzekri N, Southam L, Cooper RS, Adeyemo A, McKenzie CA, Luke A, Chen G, Elston RC, Ward R: Linkage and association analysis of angiotensin I-converting enzyme (ACE)-gene polymorphisms with ACE concentration and blood pressure. Am J Hum Genet. 2001, 68: 1139-1148. 10.1086/320104.View ArticlePubMedPubMed CentralGoogle Scholar
- Villard E, Tiret L, Visvikis S, Rakotovao R, Cambien F, Soubrier F: Identification of new polymorphisms of the angiotensin l-converting enzyme (ACE) gene, and study of their relationship to plasma ACE levels by two-QTL segregation-linkage analysis. Am J Hum Genet. 1996, 58: 1268-1278.PubMedPubMed CentralGoogle Scholar
- Rudnichi A, Safar ME, Lajemi M, Benetos A: Gene polymorphisms of the renin-angiotensin system and age-related changes in systolic and diastolic blood pressure in subjects with hypertension. Am J Hypertens. 2004, 17: 321-327. 10.1016/j.amjhyper.2003.12.011.View ArticlePubMedGoogle Scholar
- Safar ME, Lajemi M, Rudnichi A, Asmar R, Benetos A: Angiotensin-converting enzyme D/I gene polymorphism and age-related changes in pulse pressure in subjects with hypertension. Arterioscler Thromb Vasc Biol. 2004, 24: 782-786. 10.1161/01.ATV.0000119354.41615.33.View ArticlePubMedGoogle Scholar
- Goldstein H: Multilevel Statistical Models. 2003, Hodder Arnold, London, 3Google Scholar
- Rasbash J, Browne W, Goldstein H, Yang M, Plewis I, Healy M, Woodhouse G, Draper D, Langford I, Lewis T: Modelling the variance as a function of explanatory variables. A user's guide to MLwiN. 2000, Multilevels Models Project, Institute of Eduaction, University of London, 77-88. [http://www.multilevelioeacuk/download/usermanpdf]version 2.1Google Scholar
- Briollais L, Tzontcheva A, Bull S: Multilevel modeling for the analysis of longitudinal blood pressure data in the Framingham Heart Study pedigrees. BMC Genet. 2003, 4 (Suppl 1): S19-10.1186/1471-2156-4-S1-S19.View ArticlePubMedPubMed CentralGoogle Scholar
- Burton P: Applications in genetic epidemiology: modellinh correlations within nuclear families using multilevel modelling. Multilevel Modelling Newsletter. 1995, 1: 5-8.Google Scholar
- Minisymposium: The Malmo Diet and Cancer Study. Design, biological bank and biomarker programme. 23 October 1991, Malmo, Sweden. J Intern Med. 1993, 233: 39-79.View ArticleGoogle Scholar
- Berglund G, Eriksson KF, Israelsson B, Kjellstrom T, Lindgarde F, Mattiasson I, Nilsson JA, Stavenow L: Cardiovascular risk groups and mortality in an urban swedish male population: the Malmo Preventive Project. J Intern Med. 1996, 239: 489-497. 10.1046/j.1365-2796.1996.483819000.x.View ArticlePubMedGoogle Scholar
- Bog-Hansen E, Lindblad U, Bengtsson K, Ranstam J, Melander A, Råstam L: Risk factor clustering in patients with hypertension and NIDDM. The Skaraborg Hypertension Project. J Int Med. 1998, 243: 223-232.View ArticleGoogle Scholar
- Bengtsson K, Orho-Melander M, Lindblad U, Melander O, Bög-Hansen E, Ranstam J, Råstam L, Groop L: Polymorphism in the angiotensin converting enzyme but not in the angiotensinogen gene is associated with hypertension and type 2 diabetes: the Skaraborg Hypertension and Diabetes Project. J Hypertens. 1999, 17: 1569-1575.View ArticlePubMedGoogle Scholar
- Rastam L, Sjonell G: Measurement of blood pressure with tricuff. Lancet. 1991, 337: 1162-10.1016/0140-6736(91)92828-P.View ArticlePubMedGoogle Scholar
- Sambrooks J, Fritsch EF, Maniatis T: Molecular cloning: A laboratory manual. 1989, NY: Cold Spring Harbour, 2Google Scholar
- Rigat B, Hubert C, Corvol P, Soubrier F: PCR detection of the insertion/deletion polymorphism of the human angiotensin converting enzyme gene (DCP1) (dipeptidyl carboxypeptidas 1). Nucl Acids Res. 1992, 20: 1433-View ArticlePubMedPubMed CentralGoogle Scholar
- Shanmugan V, Sell KW, Saha BK: Mistyping ACE heterozygotes. PCR Methods Applications. 1993, 3: 120-121.View ArticleGoogle Scholar
- Methodology WCCfDS: ATC INdex with DDDs. 1996, Oslo (Norway), [http://aje.oxfordjournals.org/cgi/content/full/161/1/81]Google Scholar
- Snijders TAB, Bosker RJ: Multilevel analysis: an introduction to basic and advanced multilevel modeling. 1999, Thousand Oaks, California: SAGE Publications, 1Google Scholar
- Rasbash J, Browne W, Goldstein H, Yang M, Plewis I, Healy M, Woodhouse G, Draper D, Langford I, Lewis T: A user's guide to MLwiN. 2000, Multilevels Models Project, Institute of Eduaction, University of London, version 2.1Google Scholar
- Larsen K, Merlo J: Appropriate assessment of neighborhood effects on individual health -integrating random and fixed effects in multilevel logistic regression. American Journal of Epidemiology. 2003.Google Scholar
- Ben-Shlomo Y, Kuh D: A life course approach to chronic disease epidemiology: conceptual models, empirical challenges and interdisciplinary perspectives. Int J Epidemiol. 2002, 31: 285-293. 10.1093/ije/31.2.285.View ArticlePubMedGoogle Scholar
- Staessen JA, Wang JG, Ginocchio G, Petrov V, Saavedra AP, Soubrier F, Vlietinck R, Fagard R: The deletion/ insertion polymorphism of the angiotensin converting enzyme gene and cardiovascular-renal disease. J Hypertens. 1997, 15: 1579-1592. 10.1097/00004872-199715120-00059.View ArticlePubMedGoogle Scholar
- Malats N, Calafell F: Basic glossary on genetic epidemiology. J Epidemiol Community Health. 2003, 57: 480-482. 10.1136/jech.57.7.480.View ArticlePubMedPubMed CentralGoogle Scholar
- Schwartz GL, Turner ST: Pharmacogenetics of antihypertensive drug responses. Am J Pharmacogenomics. 2004, 4: 151-160. 10.2165/00129785-200404030-00002.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/7/14/prepub