Objectives To explore the effect of interaction between ACE genotype and salt intake on hypertension among Chinese Kazakhs, and to compare applications of interactions between logistic model and generalised partially linear tree-based regression (GPLTR) model.
Design Population-based cross-sectional study.
Setting Hong Dun, North Xinjiang, China.
Participants Non-consanguineous Chinese Kazakh participants (n=916, 342 men and 574 women) aged ≥30 years.
Main outcome measures Association between ACE genotype and hypertension, association between salt intake and hypertension, and interaction of ACE genotype and salt intake on hypertension in two models.
Results Associations between salt intake and hypertension were different in ACE genotype of II and ID+DD. Under the logistic models, main and interaction effects were not observed for men, but effects were present in opposite directions for women (main effect of ACE: OR=0.20, p=0.003; interaction effect: OR=1.07, p=0.027). Under the GPLTR model, Bayesian information criterion trees included both salt intake and ACE genotype as split variables. Individuals with a salt intake ≥19.5 g/day and ID+DD genotypes had a 3.99-fold (p=0.004) higher risk of hypertension compared with the II genotype for men, whereas salt intake <20.1 g/day and ID+DD genotypes had an OR=0.55 (p=0.014) compared with the II genotype for women.
Conclusions An interaction of ACE genotype and salt intake on hypertension was observed among Chinese Kazakhs but in different ways according to sex. The GPLTR model appears to be more suitable for an exploration of interactions in complex diseases.
Statistics from Altmetric.com
Strengths and limitations of this study
This is one of the few studies to explore the interaction of ACE genotype and salt intake on hypertension and is also the first such study in Chinese Kazakhs.
The design strengths of the study include the extraction of participants without consanguinity, the use of stratified analyses to test the existence of an interaction effect and a comparison of the interaction exploration results from two models.
The cross-sectional design cannot realise an investigation of the causal pathway of the association between the ACE genotype and salt intake interaction with hypertension, although possible mechanisms were suggested for the interaction effect on hypertension.
Salt intake was assessed by urinary sodium excretion and might be underestimated because salt can be removed from the body through other means.
Hypertension is a complex disorder resulting from an interaction of genetic and environmental factors.1 Blood pressure (BP) is regulated by the renin–angiotensin–aldosterone system (RAAS),2 which is heavily influenced by dietary salt intake. Indeed, dietary salt is recognised as the most important environmental contributor to hypertension.3 The most common genetic variation in RAAS is the insertion/deletion polymorphism of the ACE gene,4 and its association with BP has been reported to be dependent on sex and race/ethnicity.4–7 Although the presence of an interaction between ACE genotype and dietary salt intake has been supported by research findings,8 9 few reports have investigated the effect of this interaction on hypertension in different populations.10
Measurements and explanations of candidate gene–environment interaction (cG×E) in populations have long been explored,11 resulting in a general consensus that statistical concerns about their modelling are often overlooked.12 Popular logistic regression models including interaction terms can be problematic,13 14 but the newly introduced generalised partially linear tree-based regression (GPLTR) model, which combines tree structure and linear adjustment, may offer more flexibility and power to the exploration of interactions in complex diseases.15
The aims of this study were to explore the interaction of ACE genotype and salt intake on hypertension among Chinese Kazakhs, who have a high prevalence of hypertension and a high dietary salt intake, and to compare the logistic model and the GPLTR model in interaction applications.
The design scheme, sampling process and data collection of the Chinese Xinjiang Altay Kazakh Heart Study have been presented elsewhere.16 17 Hong Dun town is located in the urban rural fringe of Altay prefecture in North Xinjiang, China. From October 2012 to February 2013, a total of 1805 individuals aged ≥30 years were recruited, and 1668 participants from 601 families completed the baseline survey and examination under standard procedures. Participants were required to have at least three paternal/maternal biological generations living in the same region without a history of intermarriage. Non-consanguineous spouses without consanguinity and only one member of biological relatives from each family were extracted to ensure independence (figure 1). Approval for the study was given by the Ethics Committee of the Institute of Basic Medical Sciences of the Chinese Academy of Medical Sciences, and written informed consent was obtained from all participants prior to enrolment.
We complied with current Chinese Guidelines for Blood Pressure Measurement.18 A quiet room in a village clinic was selected for BP measurements, and all people in the room were required to keep quiet. Participants were comfortably seated, with legs uncrossed and the back and arm supported, such that the middle of the cuff on the upper arm was at the level of the right atrium (the mid-point of the sternum). BP was measured with an appropriate arm cuff and a mercury column sphygmomanometer on the left arm after a resting period of at least 10 min in the supine position. All measurements used the same standard procedure and calibrated mercury sphygmomanometers, and were conducted by trained and qualified staff. Currently, a large-size arm cuff for a mercury column sphygmomanometer is not available in the market of China, so radial artery auscultation on the upper forearm was used for obese individuals with large upper arms, as recommended by BP measurement guidelines.18 BP was defined as the mean of two readings, or three readings if there was a difference of more than 5 mm Hg between the initial readings. Participants were classified as hypertensive if they met at least one of the following three criteria: (1) a systolic BP ≥140 mm Hg; (2) a diastolic BP ≥90 mm Hg and (3) the use of antihypertensive medications within the last 2 weeks.
Venous blood samples were drawn after an overnight fast of at least 10 hours and centrifuged within 20 min, then divided into two. One sample was used for biochemistry assays performed within 4 hours, and the other was stored in a portable refrigerator at −20°C then transferred as soon as possible to Beijing, where it was kept at −80°C for genetic analysis. D and I ACE alleles were identified by PCR amplification of the respective fragments from ACE intron 16 and size fractionation and visualisation by electrophoresis. The ACE genotype was classified as II, ID or DD depending on the presence of each allele, with a 99.76% genotyping success rate.
The second urine sample after waking was collected, and urinary sodium was measured with a Caretium XI-921CT automatic electrolyte analyzer (Shenzhen, China). Salt intake was assessed by urinary sodium excretion per day from the amounts of urinary sodium, urine creatinine (Cr) and the 24 hours urinary Cr excreted as estimated from height, body weight and age,19 for which each 43 mmol of sodium is approximately equivalent to 1 g of sodium or 2.5 g of salt (sodium chloride).
Other covariates were also collected using standard procedures. Demographic information (including age, sex, occupation and education) and information about diet, alcohol consumption, cigarette smoking and medical history were obtained from a face-to-face questionnaire survey. Body weight and height were measured in light indoor clothes and without shoes. Body mass index (BMI) was calculated as body weight in kilograms divided by height in metres squared (kg/m2). Fasting blood glucose (FBG) was collected according to standardised protocols by a Beckman Coulter AU2700 clinical chemistry analyser (California, USA).
GPLTR model and selection criteria
was denoted as the outcome of interest, were the confounding variables to be modelled linearly and were the group of risk factors whose interaction effect was to be approximated by the tree-based structure. The main purpose of this GPLTR model was to test the combined effect of on while adjusting for the confounders ; the model was expressed as
, where is a specific link function (generalised linear model), such as the logit link function for a binary outcome, and is the vector of indicator variables signifying the leaves of the tree .
In this study, we used to represent the presence of hypertension, to be the ACE genotype and the estimated salt intake and to be the environmental risk factors having an apparent impact on the hypertension rate. First, we constructed a maximal tree with successive steps of initialization, iteration and ending conditions. The goodness of a candidate split was assessed for each node by the deviance of a generalised linear model, which equals the test deviance of the parent node minus the sum of the deviance of the two child nodes. Second, the sequence of nested candidate subtrees was created, and Akaike information criterion (AIC) and Bayesian information criterion (BIC) were used to determine the best tree in the sequence.
All analyses were conducted separately by sex. Basic characteristics were compared with regard to the status of hypertension (table 1). Variables significantly different between hypertension and normotension were chosen as confounders; therefore, smoking and family history of hypertension were not included as confounders in further analyses. ORs with 95% CIs were calculated to assess the associations of the ACE polymorphism with hypertension (table 2). To examine the presence of an interaction, associations of salt intake with hypertension were illustrated as stratified by ACE genotype as II and ID+DD (table 3). To compare the effect of the interaction, two models were applied (table 4), a logistic regression model and a GPLTR model adjusting for the same covariates to extract the high-risk group of hypertension using AIC and BIC (figure 2). Description and preliminary analyses were carried out by SAS V.9.3 software, and GPLTR models were conducted by R V.3.1.3 with the ‘GPLTR’ package. One-sided hypotheses were established for trend tests based on the logistic model, and all p values were two-tailed and considered significant at ≤0.05.
This study enrolled 916 independent participants (342 men and 574 women) who had complete data for all study variables. The prevalence of hypertension was 49.4% for men and 35.9% for women. Participant characteristics are shown in table 1. Briefly, participants with hypertension had increased age, BMI and FBG compared with normotensive participants for both sexes.
The independent associations of the ACE genotype and salt intake with hypertension are shown in table 2. The D allele had a protective effect on hypertension in women, with the DD genotype having an OR=0.47 (95% CI 0.25 to 0.88) compared with the II genotype. This protective effect was not seen in men, with the DD genotype having an OR=1.27 (95% CI 0.66 to 2.46) compared with the II genotype. Salt intake by quartiles did not correlate with hypertension for either sex (p for trend=0.099 for men and 0.178 for women).
Associations between salt intake and hypertension stratified by ACE II genotype and ID+DD genotypes were evaluated (table 3). No correlation was observed in the II genotype group for either men (p for trend=0.967) or women (p for trend=0.904). An increase in hypertension with a higher salt intake was observed in the ID+DD genotypes group for both men (p for trend=0.028) and women (p for trend=0.045).
Results from the logistic regression model adjusted for multiple confounders are shown in table 4. For men, no main effects of the ACE genotype (OR=1.33, p=0.426) or salt intake (OR=1.00, p=0.415) were observed, and there was no interaction effect (OR=0.99, p=0.464). For women, no main effect of salt intake (OR=0.99, p=0.523) was observed, but a main effect of the ACE genotype (OR=0.20, p=0.003) and an interaction effect (OR=1.07, p=0.027) were observed to be present in opposite directions. Recessive model was also built as sensitivity analysis, and ORs were turned as reciprocals as expected, which confirmed the robustness of the findings.
All iterative procedures converged after 15 iterations, and trees selected by AIC and BIC with different degrees of pruning resulted in 10 and 7 leaves for men, and 9 and 3 leaves for women, respectively (figure 2). BIC trees were chosen for the final interpretation based on their stability. Risk assessments of trees are shown in table 4. For men, leaves are identified from the left as
Node 1: 57 participants characterised by salt intake ≥19.5 g/day and ACE II genotype;
Node 2: 40 participants characterised by salt intake ≥19.5 g/day and ACE ID+DD genotype;
Node 3: 127 participants characterised by 15.3 g/day ≤ salt intake <19.5 g/day;
Node 4: 13 participants characterised by 14.6 g/day ≤ salt intake <15.3 g/day;
Node 5: 57 participants characterised by 12.2 g/day ≤ salt intake <14.6 g/day;
Node 6: 17 participants characterised by salt intake <9.6 g/day;
Node 7: 31 participants characterised by 9.6 g/day ≤ salt intake <12.2 g/day.
Split variables for the selected tree included both salt intake and ACE genotype as II and ID+DD. After adjusting for the same confounders in the multivariable logistic model, the final tree showed that Node 2 had 3.99-fold (p=0.004) higher risk for hypertension compared with Node 1, and Node 3 had an OR=0.93 (p=0.860). For women, leaves were identified from the left as
Node 1: 185 participants characterised by salt intake <20.1 g/day and ACE II genotype;
Node 2: 288 participants characterised by salt intake <20.1 g/day and ACE ID+DD genotype;
Node 3: 101 participants characterised by salt intake ≥20.1 g/day.
Split variables also included both salt intake and ACE genotype as II and ID+DD. Node 2 had an OR=0.55 (p=0.014) and Node 3 had an OR=1.26 (p=0.444) compared with Node 1.
There are three main findings from the present study. First, we observed that the association between ACE gene polymorphism and hypertension was sex specific. Second, stratified analysis indicated that an interaction between ACE genotype and salt intake existed for both sexes. Third, an application comparison of the two models suggested that the traditional logistic model was not appropriate in this situation, and that the GPLTR model could effectively determine the high-risk characteristics of hypertension.
The distribution of the ACE genotype observed in this study is similar to that reported for an isolated group of Kazakhs in 2003,20 with a DD genotype frequency of less than 20%, which differs from that observed in other studies on Kazakhs21 and in studies of other ethnicities.4–7 This difference may arise from the diverse structure of the population and racial heterogeneity. Our current results are consistent with previous studies that reported ACE genotype to be a sex-specific candidate gene for hypertension.5 6 22 23 No positive association was detected between the D allele and hypertension in men, but a protective effect of the D allele was observed for women, presumably through sexes, ethnicities and other factors including gene–gene or gene–environment interactions.6 A non-linear association between salt intake and risk of hypertension was confirmed by restricted cubic spline with automatically computed ‘knots’ (data not shown). With regard to the interaction of ACE genotype and salt intake on hypertension, a previous study that examined 284 Japanese men using logistic regression analysis10 reported a negative association for ACE (II+ID vs DD: OR=0.1, p=0.024) and a positive association for interaction [ACE (II+ID vs DD)× salt intake level: OR=3.6, p=0.047] in opposite directions, which was similar to our current results when the II genotype was set as the reference.
Strong evidence from association and linkage studies suggests that ACE has a profound effect on the variance of plasma ACE levels,24 which maintain BP homeostasis. Among the major factors that activate RAAS, dietary sodium/potassium balance is the most important and is mainly influenced by dietary salt intake.3 The sodium status has been reported to affect the ACE I/D polymorphism phenotype, and to enhance BP and renal function responses, while the conversion of aldosterone to angiotensin I with the DD genotype is blunted by low a sodium intake.8 These reports suggest that the ACE polymorphism interacts with dietary salt intake, with subsequent effects on the regulation of BP. Gene-targeting experiments that inactivate ACE in mice have provided evidence for the sex-specific effect of ACE gene on blood pressure.25 Population-based research also supports the sex-specific influence of ACE on BP variability.4 26 Assessment of the dietary habits of individuals in this study (data not shown) revealed average salt intake of 17.3 g/day, which is far greater than that of other populations and WHO recommendations.27 Given the high daily salt intake for this population, we suggest that the potential harmful effects of high salt intake on the development of hypertension should be emphasised to Chinese Kazakhs. Previous studies into the relationship between ACE and salt sensitivity indicated possession of the ACE D allele by men may lead to an increased response to a high salt intake4 28; thus, carriers of the D allele may trigger hypertension if they consume a high salt diet. For women, the mechanism behind the interaction is uncertain; however, the combination of the ACE D allele and a lower salt intake may have a synergistic effect on hypertension.28
Previous discussions into gene–environment interactions have historically progressed through three main stages, including eugenics, IQ and the serotonin transporter gene controversy.29 This has resulted in an apparent divide between variation-partitioning and mechanism-elucidation approaches. The former approach interpreted interactions as strictly statistical phenomena, as suggested by Fisher,11 while the latter emphasised the need to probe and understand the causal mechanisms behind interaction responsible for the developmental process, as suggested by Hogben.11 For complex diseases like hypertension, failure to pay proper attention to the methodological checks that must be undertaken for cG×E has resulted in some studies making dubious claims.30 31 Logistic regression requiring main effects and linear conditions is usually restricted to low-order interactions, and exploration relying on cross-product terms to measure the cause of the variation responsible for the interaction can be particularly problematic. We have built a full model using three main factors (ie, gender, ACE genotype and salt intake) in which the three main effects and the second-order interaction term are positive, but the three first-order interaction terms are all negative. Moreover, explanatory and confounding variables are entered equally into logistic regression models, which may disperse or conceal the actual main and interactive effects.32 For example, the main protective effects from the ACE genotype and the risk interaction effects of the ACE genotype and salt intake in women of the current study cannot be readily explained because the ACE genotype has no reported independent effect on hypertension and rarely shows the opposite effect after interacting with other factors.8 Therefore, new models suitable for complicated high-order interactions and appropriate confounder adjustments, designed to be distinguished from generalised linear regression models (GLMs), are warranted to reveal gene–environment interactions for complex diseases.33 The GPLTR model is semi-parametric and integrates the advantages of the tree structure and GLMs. Its purpose is to decompose a data space recursively into smaller areas defined by a set of explanatory variables with adjustment for confounding variables,15 34 35 where explanatory and confounding variables are not treated in the same way. The modified procedure15 can explicitly identify subgroups of individuals at high risk for hypertension with fewer computational demands, where highlighted features of subgroup provide clues for the mechanisms of disease development. Taking the results of the GPLTR model in the current study as an example, the ACE D allele was suggested to be the trigger for hypertension among men with a high salt intake, and female carriers of the D allele combined with a lower salt intake may offer help protection from the development of hypertension. Our calculation of the 10-fold cross-validation prediction error (data not shown) of the logistic regression model and GPLTR model for both men and women revealed the overall improved accuracy of GPLTR model.
The high-quality study design, a response rate of 92.4% and high aggregation of Chinese Kazakhs in this study provide good internal validity of the results. Sex-stratified analyses were performed to reveal potential differences in mechanisms behind the effect of the ACE genotype and salt intake interaction on hypertension. The high prevalence of hypertension (49.4% for men and 35.9% for women) allowed for robust comparisons between individuals with hypertension and non-hypertensive controls. Application comparisons between the logistic model and the GPLTR model was based on the same confounding variables by sex, which were cautiously chosen according to baseline results. The study should also be considered in light of its limitations. The cross-sectional design cannot realise investigation of the causal pathway of the association between the ACE genotype and salt intake interaction with hypertension, so confirmatory studies in other large population-based samples are warranted. The relatively small sample size may cause statistical instability. However, all enrolled participants did not have consanguinity, which is necessary for gene–environment interaction research. Salt intake was assessed by urinary sodium excretion, so it might be underestimated because salt can be removed from the body through other means, for example, sweat. Nevertheless, this sodium excretion estimation method was shown to be closely correlated with the value determined from 24 hours pooled urine.19 In addition, because of field conditions, this method was deemed superior to an estimation from a dietary questionnaire.
The interaction of ACE genotype and salt intake on hypertension was observed for both men and women among Chinese Kazakhs, although in different ways. Application comparison of the two models should be of interest to epidemiologists and geneticists for cG×E exploration, and revealed that the GPLTR model might be more suitable to help elucidate the aetiology of complex diseases. Further research is needed to determine the causes for the observed differences behind the interaction mechanisms for men and women.
We thank the government of Altay City and Hong Dun Town and the Health Bureau of Altay City who supported this project, the village head and staff of each sample village, clinical doctors and nurses of town hospital and people hospital in Altay region, Kazakh college students from Xinjiang Medical University who participated in the surveys, and all participants of this study.
Contributors YW and BZ abstracted data, performed data analysis and interpreted results with substantial contributions from the other co-authors, and drafted the manuscript, which was critically revised by all co-authors. JJ and LH conceived the study and assisted in data analysis, design of figures and tables. WH, FX and YW assisted in data analysis and interpretation, added important background knowledge and improved the manuscript by repeated readings and rephrasing as well as critical discussions of the intellectual content. MZ, YT, SL, WW, KA and CQ contributed to the data collection, quality control, general editing purposes and comments regarding medical terminology and critical review of the manuscript. ZW, YH and LW mainly performed statistical analyses, designed tables and figures and reviewed the manuscript. All authors have read and approved the final version of this manuscript for submission.
Funding This work was supported by National Natural Science Foundation of China (81273181).
Competing interests None declared.
Patient consent Obtained.
Ethics approval Committee of Institute of the Basic Medical Sciences Chinese Academy of Medical Sciences.
Provenance and peer review Not commissioned; externally peer reviewed.
Data sharing statement No additional data available.
If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.