Introduction

The endothelin (ET) family consists of three members, ET-1, ET-2 and ET-3 [1], among which ET-1, being the most frequent peptide in blood vessels, is the most salient of all. Endothelin-1 contains 21 amino acids and two disulfide bonds and its biosynthesis is a multi-step process in which the mature form of ET-1 is produced from endothelin converting enzyme (ECE) activity. The ET-1 precursor has 212 amino acids [2]. Endothelin-1 can perform vasoconstriction by binding to two receptors: ETA and ETB [3]. The endothelin system has been known to participate in the pathophysiology of many diseases such as hypertension, atherosclerosis, coronary artery disease (CAD) and heart failure [1].

Endothelin-1 and hypertension

Hypertension is a health issue increasing constantly worldwide [4], and is a well-known factor leading to cardiovascular diseases [5]. Hypertensive individuals are believed to constitute a quarter of the global adult population, and it is estimated that this proportion will increase to 30% by 2025 [6]. It has been proposed that endothelin-1, as a potent vasoconstrictor, may have an effective role in some types of hypertension [7]. Ergul et al. found an association between ET-1 and hypertension in the Afro-Caribbean race [8]. This association was also reported in patients with specific conditions such as obesity and insulin resistance [9, 10]. Asia is the most populated region and there have been numerous reports regarding the increasing risk of cardiovascular diseases in this area [3]. Iran, as an Asian country, is facing the above-mentioned risk as it has been estimated that 46% (12 million) of Iranians aged between 25 to 64 years are currently suffering from hypertension [11]. Many studies confirm the role of ET-1 in regulating the vascular tone and hypertension [7, 12, 13]; yet, this is a controversial matter [14]. Several types of single nucleotide polymorphisms (SNPs) for the ET-1 gene have been reported so far [15], two of which have been considered to be most closely associated with the pathophysiology of cardiovascular diseases – rs5370 and rs10478694.

rs5370

The rs5370 SNP contains (G/T) substitution and (Lys/Asn) conversion in codon 198 of exon 5 of ET-1 and there are some presumptions regarding its role in hypertension [14]. Genotypes are GG for the wild type, TT for the mutant and GT for heterozygotes [16]. Previous studies have produced controversial results regarding the effect of this SNP on the development of hypertension. Several studies have shown the association of rs5370 with essential hypertension [17]; however, some researchers have reported no association between the distribution of the genotypes of this SNP and the development of hypertension [18]. Tiret along with some other researchers have claimed that this association can be dependent on the body mass index (BMI) [8, 19], while a cohort study in the UK did not confirm that [18].

rs10478694

According to NCBI this is a kind of indel variation on human chromosome 6 [15], while some studies have referred to it as an SNP of human ET-1 gene [20, 21]. This SNP contains an adenine insertion at position +138, 5′ untranslated region (UTR) and exon 1 of the ET-1 [22]. The genotypes are as follows: the mutant form (4A/4A), wild type (3A/3A) and the heterozygote (3A/4A) [17]. Some studies have reported increased plasma levels of ET-1 in individuals who have a mutant genotype [17]. Studies have reported different allele distributions among patients with pulmonary artery hypertension (PAH), idiopathic pulmonary artery hypertension (IPAH) and coronary heart disease (CHD), with the control group [23]. They showed a significant increase of alleles containing the 3A form in patients with hypertension [23]. Some researchers have shown that there is no significant association between this SNP and the development of hypertension [24, 25].

Endothelin-1 and coronary atherosclerosis

Coronary artery disease is the major risk factor for myocardial infarction (MI) and heart attack. People with a high cholesterol level and high blood pressure are more prone to develop CAD due to the intima stiffness caused by atherosclerosis [26, 27]. Atherosclerosis is the main cause of cardiovascular diseases [28] and almost 50% of the etiology could be explained by major risk factors of atherosclerosis [29]. Lerman et al. detected upregulation of ET-1 in coronary atherosclerosis (CA) [30, 31]. Endothelin-1 vasoconstriction and the myogenic effect on smooth muscle cells (SMCs) may have an effect on the development of atherosclerotic lesions [32, 33]. On the other hand, some studies have reported no significant association between ET-1 plasma level and the development of atherosclerosis [34]. Therefore, due to the considerable number of controversial reports, the present study aimed to investigate the association between ET-1 plasma level and the development of hypertension and CA in a population located in the south-eastern area of Iran. Also, rs5370 and rs10478694 polymorphisms of the ET-1 gene with hypertension and CA in the considered population.

Material and methods

This was a cross-sectional study that involved 260 individuals: 128 women and 132 men with the mean age of 56 ±10.05 years and mean body mass index (BMI) level of 25 ±4.9 with suspected cardiovascular diseases (CVD). Patients were candidates for coronary angiography at Shafa Hospital of Kerman, Iran. These participants were under CAD treatment before coronary angiography, they were symptomatic and had a history of hospitalization in coronary care unit (CCU), or some evidence of having myocardial ischemia.

The study has been approved by the research ethical committee of Kerman University of Medical Sciences (approval number: IR.KMU.REC/93/268).

All patients entered the study voluntarily, filling out the written consent form. All the participants underwent clinical and laboratory evaluation, including a complete medical history and physical examination, measurement of systolic and diastolic blood pressure, a complete blood count (CBC), fasting blood sugar (FBS) and a lipid profile. Clinical blood pressure recordings were performed according to the standard technique using a mercury sphygmomanometer [35]. Hypertension was defined as having a systolic blood pressure higher than 140 mm Hg and a diastolic blood pressure more than 90 mm Hg, or having records of taking any kind of anti-hypertension remedy [28, 36, 37]. Obesity was diagnosed based on the BMI level > 30 kg/m2 [38]. The presence of coronary atherosclerosis was identified by angiography and ascertained by a qualified cardiologist (angiography was done for all participants). Coronary artery disease was defined as ≥ 50% luminal diameters stenosis in ≥ 1 epicardial coronary artery and patients without any plaque and no intimal irregularity were placed in one group and regarded as normal according to angiographic data [39]. Then, all participants were divided into 4 main groups as follows: hypertensive without atherosclerosis (H+A–) (including 35 women and 26 men); hypertensive with atherosclerosis (H+A+) (including 39 women and 37 men); non-hypertensive with atherosclerosis (H–A+) (including 20 women and 41 men); and non-hypertensive without atherosclerosis (the control group) (including 34 women and 28 men). In order to check the allele and genotype distributions, we considered another grouping as follows: hypertensive (H+A+ and H+A–); non-hypertensive (H–A+ and H–A–); individuals with CAD (H–A+ and H+A+); and individuals without CAD (H–A– and H+A–). A 10 ml peripheral blood sample was collected in vials containing ethylenediaminetetraacetic acid (EDTA) from the participants when fasting, and after 7 min of centrifugation at a rate of 2500 revolutions per minute (rpm), plasma samples were collected in separated vials and kept at –20°C. DNA extraction was done using peripheral blood lymphocytes based on the “salting out” method [40].

rs5370: polymerase chain reaction – restriction fragment length polymorphism (PCR-RFLP)

Polymerase chain reaction (PCR) was processed (Corbett Research thermocycler) using the following primers based on a previous study on ET-1 [41]. Forward: 5′-TCTTGCTTTATTAGGTCGGAGACC-3′, reverse: 5′-TTTGAACGAGGACGCTGGTC-3′. Cycle conditions were as follows: initial denaturation 95°C for 5 min, 35 cycles of 95°C for 45 s, 59°C for 50 s, 72°C for 45 s, and finally 72°C for 10 min for a final extension. PCR product was a 262 base pair (bp) fragment which was incubated for 16 h in 37°C with restriction enzyme Cac8i (New England Biolabs) and the digested products were electrophoresed in 1.5% agarose gel.

rs10478694: PCR-RFLP

Primers were designed using Vector-NTI version 11.5 – forward: 5′-TCTTGCTTTATTAGGTCGGAGACC-3′, reverse: 5′-TTTGAACGAGGACGCTGGTC-3′. Cycle conditions were as follows: initial denaturation 95°C for 5 min, 35 cycles of 95°C for 45 s, 52°C for 50 s, 72°C for 50 s, and finally 72°C for 10 min for a final extension. The polymerase chain reaction product was a 262 bp fragment which was incubated for an hour at 55°C with the restriction enzyme Bsl1 (New England Biolabs) and the digested products were electrophoresed in 1.5% agarose gel.

Endothelin-1 measurement

Endothelin-1 of all participants was measured using enzyme linked immunosorbent assay (ELISA) Kit (R&D systems, Pan specific DuoSet DY1160), and measurement was carried out according to the manufacturer’s protocol.

Statistical analysis

All data for continuous variables are shown as mean ± standard deviation (M ± SD) and categorical variables are presented as frequency (percentages). The qualitative variables of the four groups were compared using the χ2 test (or Fisher’s exact test) and quantitative variables were compared using one-way ANOVA. The ET-1 plasma level of the four groups was compared by ANCOVA test. By use of Fisher’s exact test, the genotype and allele distributions between the hypertensive and non-hypertensive groups were compared and also the same comparison was made between people with coronary atherosclerosis and people without coronary atherosclerosis. Hardy-Weinberg equilibrium (HWE) for each polymorphism was assessed by χ2 test. The association of ET-1 plasma level and CAD was investigated by Pearson’s correlation coefficient. Statistical analysis was performed using SPSS version 22 (IBM Corp., Armonk, NY, US).

P-values less than 0.05 were considered statistically significant.

Results

Table I summarizes participants’ quantitative and qualitative characteristics. The average age among members of the H+A+ group was significantly higher than that of the other groups (p = 0.006). The frequency of women in groups H–A+ and H+A+ was significantly higher than in the other groups. Therefore, using the ANCOVA test, ET-1 was compared among the groups with controlled age and gender. The outlier data were identified and excluded using Boxplot. The averages of the other variables were not significantly different among the groups. The genotype and allele frequencies at both SNPs did not confirm the HWE (data not shown).

Table I

Comparison of qualitative and quantitative variables among the four groups under investigation

ParameterH–A– (N = 62)H+A– (N = 61)H–A+ (N = 61)H+A+ (N = 76)Statistical results
GenderFemale20 (32.5)*21 (35.0)44 (72.0)43 (57.5)χ2 = 17.5
p < 0.001
Male26 (42.5)16 (27.5)39 (65.0)51 (67.5)
Angiography informationNone62 (100)61 (100)13 (22.5)19 (25)χ2 = 100.2
p < 0.001
1 vessel0 (0)0 (0)20 (32.5)30 (40)
2 or more vessels0 (0)0 (0)27 (45)26 (35)
AddictionCigarettes7 (12.5)7 (12.5)13 (22.5)13.3 (17.5)χ2 = 2.02
p = 0.574
Opium7 (12.5)6 (10)15 (25)11 (15)χ2 = 3.93
p = 0.202
Others1 (2.5)1 (2.5)4 (7.5)2 (2.5)χ2 = 2.08
p = 0.551
Age [years]52.3 3 ±9.93**54.68 ±8.6857.30 ±9.0559.75 ±11.19F = 4.354
p =0.006
BMI [kg/m2]25.0 ±4.8726.68 ±5.6524.11 ±4.2524.37 ±4.48F = 2.27
p = 0.822
FBS [mg/dl]92.93 ±11.0594.38 ±10.2693.08 ±10.7192.23 ±9.74F = 0.29
p = 0.822
Cholesterol [mg/dl]155.28 ±41.24151.98 ±32.77149.20 ±37.39144.63 ±44.72F = 0.52
p = 0.662
TG [mg/dl]153.68 ±30.80170.63 ±44.18157.33 ±28.17175.20 ±69.59F = 2.00
p = 0.111
HDL [mg/dl]40.20 ±11.5938.25 ±10.2036.18 ±11.0837.60 ±9.51F = 0.99
p = 0.403
LDL [mg/dl]96.77 ±32.6289.51 ±24.5793.32 ±27.6093.19 ±29.83F = 0.42
p = 0.732
Urea [mg/dl]33.75 ±13.7135.98 ±10.1636.60 ±11.8735.28 ±10.04F = 0.45
p = 0.711
Creatinine [mg/dl]1.04 ±0.24.99 ±0.241.11 ±0.211.08 ±0.23F = 1.92
p = 0.131
Na [mg/dl]139.50 ±3.29139.70 ±3.44139.95 ±2.68139.90 ±3.24F = 0.16
p = 0.923
K [mg/dl]4.34 ±0.364.36 ±0.464.41 ±0.434.36 ±0.42F = 0.18
p = 0.911
SBP [mm Hg]118.88 ±15.38120.53 ±13.34116.13 ±12.58121.13 ±14.52F = 1.02
p = 0.381
DBP [mm Hg]74.00 ±8.4973.83 ±8.0670.63 ±8.1873.13 ±8.98F = 1.45
p = 0.232
Endothelin-1 [pg/ml]5.27 ±0.814.13 ±0.997.36 ±1.79.79 ±1.6F = 3.47
p = 0.018

[i] The average age among members of the H+A+ group was significantly higher than that of the other groups. Averages of the other quantitative variables were not significantly different among the groups. The number of women in groups H–A+ and H+A+ was significantly higher. Endothelin-1 plasma level was significantly higher in the H+A+ group in comparison with the other groups (p = 0.018). BMI – body mass index, DBP – diastolic blood pressure, FBS – fasting blood sugar, H–A– – no coronary artery disease and no hypertension, H+A– – hypertension and no coronary artery disease, H–A+ – coronary artery disease and no hypertension, H+A+ – coronary artery disease and hypertension, HDL – high-density lipoprotein, LDL – low-density lipoprotein, SBP – systolic blood pressure, TG – triglycerides.

Endothelin-1 plasma level

The comparison of plasma ET-1 level among the groups can be generally reported as follows: H+A– < H–A– < H–A+ < H+A+ (p = 0.018) (Figure 1). The Tukey results showed a significant difference between ET-1 plasma level in H+A+ (4.13 ±1.6 pg/ml) and H+A– (9.79 ±0.99 pg/ml) (p = 0.018).

Figure 1

Comparison of endothelin-1 plasma level among the four groups under investigation Endothelin-1 plasma level in the H+A+ group is significantly higher than in the other groups (p = 0.018). Only the difference of endothelin-1 plasma level observed between H+A+ and H+A– is significant (p = 0.017). H+A– – hypertension and no coronary artery disease, H–A+ – coronary artery disease and no hypertension, H+A+ – coronary artery disease and hypertension, H–A– – no coronary artery disease and no hypertension. Significant difference at < 0.05

https://www.archivesofmedicalscience.com/f/fulltexts/115937/AMS-17-3-91675-g001_min.jpg

Frequency of rs5370 and rs10478694 genotypes and alleles among groups

When analyzing rs5370, we observed that all four groups contained higher frequencies of GG genotype but comparing the allele and genotype frequencies it was found that the frequency of the T allele (p = 0.003) and TT genotype (p = 0.010) was significantly higher in people with atherosclerosis than in people without atherosclerosis (Table II).

Table II

Distribution of allele and genotype frequencies of rs5370

GroupSingle nucleotide polymorphismN
Allele frequencies (percent)Genotype frequencies (percent)
TGP-valueTTGTGGP-value
Hypertensive (H+A+ and H+A–)24 (9)122 (51)0.50115 (6)9 (3)113 (43)0.563137
Non-hypertensive (H–A+ and H–A–)25 (10)104 (40)19 (7.5)6 (2.5)98 (38)123
With coronary atherosclerosis (H–A+ and H+A+)37 (13)111 (41)0.00326 (9.6)11 (4)100 (38.5)0.011137
Without coronary atherosclerosis (H–A– and H+A–)14 (4)113 (42)10 (4)4 (1.5)109 (42.4)123

[i] Values are presented as number (percentage). Significant difference at < 0.05. The frequency of TT genotype is significantly higher in individuals with coronary atherosclerosis in comparison with those without coronary atherosclerosis (p = 0.011). The frequency of the T allele (Asn) is significantly higher in people with coronary atherosclerosis in comparison with those without coronary atherosclerosis (p = 0.003). There is no significant difference in allele and genotype frequencies of rs5370 between hypertensive and non-hypertensive individuals (p = 0.501, p = 0.563). N – indicates the total population of each section, H–A– – no coronary artery disease and no hypertension, H+A– – hypertension and no coronary artery disease, H–A+ – coronary artery disease and no hypertension, H+A+ – coronary artery disease and hypertension.

Analysis of rs10478694 showed significantly higher frequencies of 3A3A genotype (p = 0.007) and 3A allele (p = 0.004) in hypertensive people in comparison with non-hypertensive individuals (Table III).

Table III

Distribution of allele and genotype frequencies of rs10478694

GroupAllele frequencies (%)Genotype frequencies (%)N
4A3AP-value3A4A3A3A4A4AP-value
Hypertensive (H+A+ and H+A–)82 (31)93 (35)0.00438 (15)55 (21)44 (17)0.007137
Non-hypertensive (H–A+ and H–A–)98 (37)59 (22)34 (13)25 (9.5)64 (24.5)123
With coronary atherosclerosis (H–A+ and H+A+)90 (34)80 (30)0.86233 (12.7)47 (18)57 (22.2)0.241137
Without coronary atherosclerosis (H–A– and H+A–)90 (34)73 (28)40 (15.4)33 (12.7)50 (19)123

[i] Values are presented as number (percentage). Significant difference at < 0.05. The frequency of 3A3A genotype is significantly higher in hypertensive individuals in comparison with non-hypertensive ones (p = 0.007). The frequency of the 3A allele is significantly higher in hypertensive individuals in comparison with non-hypertensive ones (p = 0.004). There is no significant difference in allele and genotype frequencies of rs10478694 between individuals with coronary atherosclerosis and those without coronary atherosclerosis (p = 0.862, p = 0.241). N – indicates the total population of each section, H–A– – no coronary artery disease and no hypertension, H+A– – hypertension and no coronary artery disease, H–A+ – coronary artery disease and no hypertension, H+A+ – coronary artery disease and hypertension.

No association between ET-1 plasma level and CAD was found (p = 0.625). No association between ET-1 plasma level and genotypes in either of the two SNPs were seen (p = 0.332, p = 0.102).

Discussion

In this study, we provided evidence that ET-1 plasma level is positively correlated with the risk of CA development. The TT genotype and Asn (T) allele of rs5370 may be associated with the development of CA. The 3A3A genotype and 3A allele of the rs10478694 may be associated with the risk of hypertension.

According to the higher level of plasma ET-1 in the H+A+ group in comparison with the group H+A–, we can conclude that ET-1 is positively associated with the development of atherosclerosis. These findings also confirmed the results of several other studies [710]. Although Parinello et al. and Sarafidis and Bakris declared that the association of plasma level of ET-1 with the development of hypertension was dependent on specific conditions such as obesity and insulin resistance [9, 10], it has been reported in some studies that obesity can affect expression of some genes and cause metabolic imbalance which can lead to various diseases [42]. In this case there might be no connection between the increasing level of ET-1 in the blood and the development of hypertension; yet, in our study (as the BMI level did not exceed 30 kg/m2) obesity was not considered an influencing factor. Yoshibayashi et al. claimed that the effect of this peptide on the development of essential hypertension is clear due to the higher level of prepro ET-1 mRNA in the serum of people suffering from this disease; however, their serum level of ET-1 has not been reported [7]. Katona’s study shows that ET-1 plasma level was not significantly elevated in hypertension patients [43].

The ANCOVA and Tukey tests showed only a significant difference between H+A+ and H+A–; therefore, we concluded that ET-1 plasma level may be associated with the development of coronary atherosclerosis. In a healthy vessel, there is an equilibrium between ET-1 activities through ETA and ETB, which will help to maintain the endothelial function. If induced by the specific factors [44], ET-1 upregulation and elevated secretion can affect the endothelial function due to its known effect on the elevation of production of superoxides [36]. By disrupting nitric oxide synthase activity, superoxides can reduce the generation of nitric oxide and therefore its effects on the vascular balance and endothelial function [45, 46]. One of the main consequences of endothelial dysfunction is atherosclerosis [36, 46]. It is mentioned in Kawashima’s study that NO acts as an anti-atherosclerotic factor by controlling the vascular tone, preventing monocytes from sticking to the endothelial surface, preventing proliferation and immigration of smooth muscle cells [47].

The analysis of the distribution of rs5370 alleles and genotypes revealed significantly higher frequencies of the TT genotype and T allele (Asn) in people with CA in comparison with those without atherosclerosis (p = 0.011, p < 0.001). Yet, the present results showed no significant association between rs5370 allele and genotype distribution and the development of hypertension (p = 0.501, p = 0.563). Some researchers believe that this SNP is not located on a regulatory region of the ET-1 gene and therefore it is unable to have an effect on ET-1 gene expression [19]. According to the significantly higher frequency of Asn (T) and TT genotype in people with atherosclerosis in comparison with those without atherosclerosis (Table II), we concluded that this allele can be regarded as a risk factor causing the development of atherosclerosis.

As shown in Table II, the 3A allele and 3A3A genotype frequencies of rs10478694 were significantly higher in hypertensive individuals (p = 0.007, p = 0.004). Therefore, it can be deduced that the 3A allele and 3A3A genotype of this SNP can be regarded as risk factors that contribute to the development of hypertension. Although these findings are parallel to the findings of other studies [17, 23], some studies have shown contradictory results claiming that this SNP has no association with the development of hypertension [19, 24, 25]. Yet, the 3A allele and 3A3A genotype showed no association with the development of atherosclerosis. We found no association between ET-1 plasma level and the stage of atherosclerosis (the number of stiffened vessels). We also observed no association between ET-1 plasma level and genotypes in either of the two SNPs.

The strength of this study was the analysis of ET-1 plasma level and its two effective SNPs under different conditions, which revealed some important insights regarding ET-1 in hypertension and CA. However, this study was performed in one center, which may limit the external validity of current findings.

It is suggested to conduct similar studies on other races and populations worldwide to attain a more comprehensive understanding of the association of endothelin-1 and its polymorphisms with the development of hypertension and CA.

In conclusion, based on the present findings, ET-1 plasma level is positively correlated with the risk of CA development. The T allele can be regarded as a risk factor for CA development. The 3A allele (3A/3A genotype) of rs10478694 may be associated with the risk of hypertension in the south-eastern populations of Iran.