Abstract

Cardiovascular development critical genes are key determinants in cardiovascular diseases. We hypothesize that SNPs in these genes may play critical roles in the development of hypertension. Therefore, we enrolled 516 paired hypertension patients and controls in a total of 2,742 subjects in a cross-sectional population study by the propensity score matching (PSM) method. Twenty-one SNPs from 5 cardiovascular developmental related genes were detected by the improved multiplex ligase detection reaction (iMLDR) method. Conditioned logistic regression under three different genetic models, namely, additive model, dominant model, and recessive model, was performed. The odds ratio (ORs) and 95% confidence intervals (95% CIs) were used to estimate the associations of SNPs with hypertension. We found that the distribution of genotypes at rs833061, rs3025010, and rs699947 within the VEGFA gene and the distribution of alleles at rs3025010 in hypertension subjects were different from those in controls. Both rs833061 and rs3025010 were associated with hypertension in crude models, but only rs3025010 remains associated with hypertension after adjusting with confounding factors in the additive model and the dominant model. We also found that hypertension subjects with C/T and C/C genotypes at rs3025010 had lower SBP and DBP levels. In addition, rs3025010 could interact with rs6784267 within the CCM3 gene in the association. In conclusion, our findings suggest that rs3025010 may play a role in the pathogenesis of hypertension, which may be a potential target for individualized prevention and treatment of hypertension.

1. Introduction

Hypertension is a common chronic disease and a major risk factor for stroke, myocardial infarction, and kidney failure [1, 2]. Single nucleotide polymorphisms (SNPs) are DNA sequence polymorphisms caused by a single nucleotide variation at the genome level and are the most common human genetic variants [3]. Genome-wide association studies (GWASs) have enabled the discovery of common genetic variation contributing to normal and pathological traits [4, 5], suggesting that the use of SNPs as biomarkers is useful for the screening of susceptible populations and is conducive to genetic diagnosis, gene therapy for hypertension [6, 7].

It is important to note that cardiovascular development critical genes are key determinants in cardiovascular diseases [8, 9]. Animal model studies showed that the deficiency of them is lethal. Taken CCM3 as an example, complete CCM3 knockout mice died early, whereas heterozygous deletion could survive and grow like wild-type, which suggests that it may act as a hidden primer for cardiovascular diseases in adulthood [10, 11]. However, there have been only limited significant loci of these genes which were found in GWASs [12]. Interestingly, we recently demonstrated that CCM3 gene polymorphism was associated with elevated blood pressure based on an arsenic-exposed population [13]. Therefore, we hypothesize that SNPs in cardiovascular developmental genes may play critical roles in the development of hypertension.

Propensity score matching (PSM) is a commonly used research method to solve the confounding bias in observational investigations [14]. This method was firstly proposed by Paul Rosenbaum and Donald Rubin in 1983 [15]. In recent years, it has been widely used in observational studies of nonrandomized large samples [1618]. PSM method displays power improvement in identifying causal SNPs over other state-of-the-art methods. We employed PSM method and conditional logistic regression analysis in a population-based cross-sectional study to identify whether the putative functional SNPs are associated with hypertension in cardiovascular developmental critical genes.

2. Methods

2.1. Populations

This study was approved by the Medical Ethics Committee of the School of Public Health, Sun Yat-sen University. The subjects in this study were adults who had physical examination in a physical examination center in Shenzhen city from July 2013 to January 2014. They were mainly employees of government departments, enterprises, and institutions. The individual data of the subjects, including age, gender, BMI, hypertension family history (more than one member of the immediate family has hypertension), smoking (current smoking or the cumulative smoking time ≥6 months), drinking (alcohol consumption ≥3 g/day), salt intake, and physical activity (moderate or heavy intensity activities lasting more than 10 minutes in work and household activities, ≥1 day/week) were collected by specially trained investigators through questionnaires. The clinical examination of blood pressure, ultrasound, and biochemical parameters was measured by professional medical doctors and nurses. Subjects with severe anemia, severe heart and kidney insufficiency, urinary tract infection, nephritis, endocrine disorders, and pregnant women were not included. All subjects were divided into hypertension and control groups according to blood pressure levels and clinical information. All individuals signed informed consent.

After questionnaire quality screening, clinical examination, medical history consultation, medical history review, and doctor’s judgment, 2,392 effective subjects were selected. Among them, 258 subjects with hypertension were finally included in this study, and 258 healthy controls were matched by the propensity score method in 1 : 1 manner. A flow chart of the design of this study is illustrated in Figure 1.

2.2. Blood Pressure Measurement

Blood pressure was measured by standard procedures. The classification and definition of blood pressure were referred to the WHO/ISH guidelines 2003. Systolic blood pressure ≥140 mmHg and (or) diastolic blood pressure ≥90 mmHg was defined as hypertension. In addition, individuals who reported taking antihypertensive drugs were also considered as hypertension patients. The blood pressure level to 120–139/80–89 mmHg was set as a normal high value, with systolic pressure <120 mmHg and diastolic pressure <80 mmHg as ideal blood pressure. Subjects with both normal high value and ideal blood pressure were defined as normal.

2.3. Propensity Score Matching (PSM)

The covariates used for the calculation of the propensity score were age, sex, BMI, waist circumference, smoking history, drinking history, physical activity level, salt intake, total cholesterol, triglycerides, and blood glucose. Subjects with hypertension and healthy controls were matched at a ratio of 1 : 1 by the propensity score. In brief, open PSM module in SPSS first, then select nearest neighbour matching, enter ID variable as the serial number, enter binary treatment indicator (0 = control and 1 = hypertension), enter covariates for scoring, and finally match individuals exactly based on selected variables. The matching efficiency was assessed by outputs such as relative multivariate imbalance L1 metric and propensity score distribution histograms. The L1 metric is theoretically between 0 and 1. The smaller the L1 metric, the better the matching result. And if the distribution histograms between the postmatching group and the control group is similar, it indicates that matching is good [18, 19]. Given the association with hypertension in the logistic work model, we postulated that quantitative blood pressure levels would be lower in individuals with the minor allele at expected SNPs. Thus, we further studied the effect of SNPs on blood pressure traits in hypertension patients. For these subjects, the PS was matched in a ratio of 1 : 1 : 1 among the three genotypes by hand.

2.4. Selection and Genotyping of SNPs

DNA extraction was referred to the blood genomic DNA extraction kit instructions (Tiangen Biotech, Beijing). The tag SNPs were verified in HapMap Project and selected by using Haploview software (version 4.2) with the following filters: minor allele frequency (MAF) > 0.1 and r2 > 0.8. Finally, there were totally 21 SNPs selected, which include 3 sites (rs9818496, rs3804610, and rs6784267) from CCM3, 3 sites (rs55805015, rs2277538, and rs3212278) from DLL4, 4 sites (rs7667298, rs2305948, rs13109660, and rs7671745) from KDR, 2 sites (rs3124591 and rs73668310) from Notch1, and 9 sites (rs833061, rs10434, rs833069, rs3025010, rs3025053, rs699947, rs2146323, rs3025035, and rs3025030) from VEGFA. The specific information of these SNPs is listed in Supplemental Table 1. Genotyping was done by the improved multiplex ligase detection reaction (iMLDR) method (Genesky Biotech, Shanghai) as previously described [20].

2.5. Statistical Analysis

Quantitative data were expressed as mean ± SD (standard deviation), which were compared using the t-test or one-way ANOVA using the χ2 test. The difference in the distribution of genotype frequencies between the two groups was tested using χ2 test. Conditioned logistic regression analysis was conducted to investigate the association of the differentially distributed SNPs with hypertension under three different genetic models including additive model, dominant model, and recessive model. In conditioned logistic regressions, the odds ratio (ORs) and 95% confidence intervals (95% CIs) were used to estimate the associations of SNPs with risk of hypertension. We also adjusted for SNPs that changed the matched regression coefficients by at least 10 percent. General linear regression was also performed to investigate the potential interaction among SNPs on blood pressure levels. The process of PSM was conducted by the PSM module in SPSS, and all statistical analyses were performed by SPSS version 20.0 software. was considered as the statistically significant level.

3. Results

3.1. Basic Characteristics

A total of 2,742 questionnaires were issued upon the physical examination, and 2,392 samples including complete questionnaires and physical examination data were collected (response rate, 87%). Among them, there were 258 subjects of hypertension and 2,134 subjects of normal blood pressure. The hypertension group and the control group (including normal high values and normal blood pressure) were 1 : 1 matched, and a total of 516 subjects were enrolled in this study. The matched L1 statistic is 0.050, less than 0.522 (value before matching), suggesting that the match was good (Supplemental Table 2). After matching, the distribution of hypertension and control groups was similar, further suggesting a good validity of matching (Supplemental Figure 1). Besides, we estimated the test power of our sample size when α level equals 0.05 and OR equals 1.5, and the results indicated that our study could provide 58.99% and 88.79% power for detection of genetic variation with a MAF of 0.1 and 0.5, respectively (Supplemental Doc 1). Systolic blood pressure, diastolic blood pressure, and pulse pressure in the hypertension group were significantly higher than those in the control group (Table 1).

3.2. The Distribution of Genotypes and Alleles

Among the polymorphism sites submitted, 19 SNP sites were analyzed since there was only 1 type of genotype in all subjects from rs73668310 (Notch1) and rs55805015 (DLL4) sites. After removing unreliable samples in the genotyping step, we finally obtained 225 pairs of good quality genotyping data. The quality control samples had the same reproducibility, and the negative controls had no band. All SNPs’ response rates were greater than 98%, and the minimum alleles were consistent with HapMap-HCB data (Supplemental Table 3). The results showed that the distribution of rs833061, rs3025010, and rs699947 genotypes in the VEGFA gene were significantly different between the hypertension group and the control group () (Table 2). The allele frequency distribution also showed that there were significant differences in rs3025010 loci () (Table 3).

3.3. The Association between SNPs and Hypertension

Because rs833061 and rs699947 were linearly correlated, we only studied rs833061. In additive model, the univariate analysis results showed that rs833061 and rs3025010 were associated with hypertension (), but only rs3025010 remained associated with hypertension in multivariate analysis after adjusting with confounding factors (OR = 0.103, 95% CI: 0.019–0.551, ). In the dominant model, the univariate analysis results showed that neither rs833061 nor rs3025010 was associated with hypertension (), but rs3025010 became associated with hypertension in multivariate analysis after adjusting with confounding factors (OR = 0.354, 95% CI: 0.147–0.854, ). In the recessive model, although rs833061 and rs3025010 were associated with hypertension, neither of the associations were significant in multivariate analysis after adjusting with confounding factors (). These results suggest that rs3025010 (CC genotype) may have an independent effect on hypertension (Table 4).

3.4. The Effect of SNPs on Blood Pressure Levels

In hypertension subjects with no hypertensive treatment, a total of 36 matched subjects were obtained, with 12 samples of each genotype at rs3025010 (TT, CT, and CC). We found subjects with CT and CC genotypes at rs3025010 had lower SBP and DBP levels although there were significant differences when compared with the T/T genotype (Figures 2(a) and 2(b)). In control subjects, a total of 75 matched subjects were obtained, with 25 samples of each genotype at rs3025010 (TT, CT, and CC). We found the levels of SBP and DBP were slightly lower in subjects with CT and CC genotypes at rs3025010, and only the difference of SBP was significant for the CT genotype when compared with the TT genotype (Figures 2(c) and 2(d)).

3.5. The Interactive Effect between SNPs in the CCM3 Gene and rs3025010 on Hypertension

Although the distribution of genotypes and alleles in SNPs from the CCM3 gene was not different between hypertension and control groups, our previous study found that CCM3 gene polymorphism is associated with hypertension in an occupational population exposed to low-level arsenic, and SNPs in the CCM3 gene could interact with it, suggesting that the CCM3 gene may assist with other factors. Therefore, we further examined the interaction between rs3025010 and SNPs of CCM3 genes. Among three SNPs in CCM3, we found that rs6784267 interacts with rs3025010 in the additive model and the dominant model but not in the recessive model. In particular, rs3025010 (CC) interacted with rs6784267 (C/T) in the additive model, and rs3025010 (C/T + C/C) interacted with rs6784267 (C/T + T/T) in the dominant model even after adjusting with confounding factors (Table 5).

4. Discussion

In this study, we found that rs3025010 within the VEGFA gene was associated with hypertension, and hypertension subjects with C/T and C/C genotypes at rs3025010 had lower SBP and DBP levels. In addition, rs3025010 could interact with rs6784267 within the CCM3 gene and further affecting hypertension.

Identifying putative functional SNPs associated with a disease can provide targets for gene therapy as early as possible, opening therapeutic strategies for precision and preventive medicine [21, 22]. Whole genome sequencing and screening for susceptible populations reported a suggestive association of SNPs and hypertension [2325]. In addition, emerging technologies have been developed to repair DNA by editing the wrong bases [26, 27]. However, less is known about putative functional SNPs in the critical genes that are associated with hypertension. According to our preliminary work and literature review, cardiovascular development critical genes such as CCM3, DLL4, KDR, Notch1, and VEGFA contain multiple SNP sites, and bioinformatics data suggest that certain SNP sites may affect protein expression and function. In this study, with the assistance of the PSM method to adjusting certain definite risk factors for hypertension, we found that three SNPs in the VEGFA gene named rs833061, rs3025010, and rs699947 were differential-distributed in hypertension, indicating that these three SNPs could be potential independent risk factors in hypertension.

VEGFA regulates vascular development during embryogenesis and organ formation under physiological condition [28]. Disruption of this gene in mice resulted in abnormal embryonic blood vessel formation. This gene is upregulated in many known tumors, and its expression is correlated with tumor stage and progression [29, 30]. A recent study showed that polymorphisms in the VEGFA gene may affect the antihypertensive responses to enalapril [31]. In particular, A/A genotype at the rs699947 site is associated with a more intense decrease in blood pressure in response to enalapril, underlying its protective effect in hypertension, consistent with our results observed among rs699947, rs833061, and rs3025010. Chen et al. found that rs3025010 is possibly associated with a reduced risk of human brain arteriovenous malformation [32]. Interestingly, in our study, we found that rs3025010 was negatively associated with hypertension, and C/T and C/C genotypes which contain minor allele had lower SBP and DBP levels. This consistency suggests that the rs3025010 variant may play similar roles in cardiovascular diseases. To the best of our knowledge, this study is the first direct report of the association of VEGFA gene polymorphism at rs3025010 with hypertension.

The vast majority of GWAS tag SNPs lie in intergenic or intronic regions (approximately 88%) and therefore are likely to influence gene regulations [33]. rs3025010 is located in the fifth intron. Based on searching results from ENCODE (data not shown), we found that there were some candidate cis-regulatory elements (ccREs) near this locus including H3K27ac, which is an enhancer marker, suggesting that this locus region has the function of regulating target gene expression. Our study shows that the polymorphism of this locus is associated with hypertension, and the blood pressure level varies among different genotypes, suggesting that the polymorphism of this locus may affect the expression of genes related to blood pressure regulation through regulating the activity of enhancers or binding transcription factors.

We consider that SNPs in CCM3 and other genes may affect blood pressure by interacting with the VEGFA gene. Furthermore, our previous study found that the CCM3 gene may represent a novel susceptibility gene for hypertension in a population with arsenic exposure [13]. Therefore, we step forward to analyze the SNP-SNP interactions in this study. As expected, we found rs6784267 in the CCM3 gene interacted with rs3025010 in the VEGFA gene while both of our studies did not found associations between CCM3 gene polymorphism and hypertension, which is in line with our previous hypothesis that the CCM3 gene could affect blood pressure in a manner of interacting with other factors.

One strength of our study is that we used propensity score matching to adjust many covariates to produce unbiased estimates of the treatment effects [14]. Propensity score matching can also be very cost-effective [34]. Besides, although the sample size was relatively small, we could provide 58.99%∼88.79% power for detection of genetic variation for the selected tag SNPs. Our study also has limitations since we mainly examined the association between SNPs with hypertension. SNP site in the genomic DNA loci will affect the expression and function of the target gene and then further lead to cardiovascular dysfunction. Thus, further study is needed to explore its function and molecular mechanisms.

In summary, our findings suggest that rs3025010 in VEGFA is associated with hypertension. We also observed lower levels of SBP and DBP for C/T and C/C genotypes at rs3025010 in hypertensive subjects. In addition, rs3025010 could interact with rs6784267 within the CCM3 gene to affect the risk of hypertension. In depth studies are warranted to determine the target gene and investigate its potential role in the development of hypertension. This may improve our understanding of the role of cardiovascular development-related genes in the maintenance of blood pressure homeostasis. These findings will shed new light on gene diagnosis, gene therapy, and individualized therapy in hypertension. The results and discussion may be presented separately, or in one combined section, and may optionally be divided into headed subsections.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Disclosure

Zhiqiang Zhao and Chunmei Gong are the co-first authors.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Authors’ Contributions

Zhiqiang Zhao and Chunmei Gong contributed equally and made the greatest contribution to this work.

Acknowledgments

Yun He was supported by the National Natural Science Foundation of China (nos. 81472998 and 81872661) and the Open Project of Key Laboratory of Prevention and Treatment of Cardiovascular and Cerebrovascular Diseases, Ministry of Education (no. XN201802). Xiaoli Liu was supported by National Key Technology Research and Development Program of China (#2012BAI02B02), and Chunmei Gong was supported by Research Program from Health and Family Planning Commission of Shenzhen Municipality (#SZGW2017015).

Supplementary Materials

Table 1 shows detailed SNP information. Figure 1 and Table 2 show the analysis results of propensity score matching. Doc 1 shows the output results of power tests. Table 3 shows the quality assessment result of genotyping. (Supplementary Materials)