- Open Access
Genetic variance in Nitric Oxide Synthase and Endothelin Genes among children with and without Endothelial Dysfunction
Journal of Translational Medicine volume 11, Article number: 227 (2013)
The presence of endothelial dysfunction (ED) constitutes an early risk factor for cardiovascular disease (CVD) in children. Nitric oxide (NO) and endothelin (EDN) are generated in endothelial cells and are critical regulators of vascular function, with ED resulting from an imbalance between these two molecules. We hypothesized that genetic variants in NO synthase and EDN isoforms and its receptors (EDNRA and EDNRB) may account for a proportion of the risk for ED in developing children.
Consecutive children (ages 5–10 years) were prospectively recruited from the community. Time to peak post-occlusive reperfusion (Tmax) was considered as the indicator of either normal endothelial function (NEF; Tmax < 45 sec) or ED (Tmax ≥ 45 sec). Lipid profiles, high sensitivity C-reactive protein (hsCRP), fasting glucose and insulin were assayed using ELISA. Genomic DNA from peripheral blood was extracted and genotyped for NOS1 (209 SNPs), NOS2 (122 SNPs), NOS3 (50 SNPs), EDN1 (43 SNPs), EDN2 (48 SNPs), EDN3 (14 SNPs), EDNRA (27 SNPs), and EDNRB (23 SNPs) using a custom SNPs array. Linkage disequilibrium was analyzed using Haploview version 4.2 software.
The relative frequencies of SNPs were evaluated in 122 children, 84 with NEF and 38 with ED. The frequencies of NOS1 (11 SNPs), and EDN1 (2 SNPs) were differentially distributed between NEF vs. ED, and no significant differences emerged for all other genes. Significant SNPs for NOS1 and EDN1 SNPs were further validated with RT-PCR.
Genetic variants in the NOS1 and EDN1 genes appear to account for important components of the variance in endothelial function, particularly when concurrent risk factors such as obesity exist. Thus, analysis of genotype-phenotype interactions in children at risk for ED will be critical for more accurate formulation of categorical CVD risk estimates.
Obesity is one of the world’s greatest public health challenges affecting not only developed, but also developing countries. The presence of obesity has now been conclusively linked to heightened morbidity and mortality through increased risk for many chronic diseases, including type-2 diabetes, hypertension, dyslipidemia, and coronary artery disease [1–4]. The pathogenesis of obesity is multifactorial incorporating both genetics and lifestyle. In children, obesity is associated with increased risk for multifaceted derangements in metabolic and cardiovascular function, including endothelial dysfunction (ED) [5–8]. In general terms, overweight children are more likely to prematurely develop ED, hypertension and type 2 diabetes, an array of conditions that would normally be only found in older obese adults . Recently, we have shown that obesity in children is associated with an increased risk for the development of ED prior to the onset of hypertension . However, not every obese child will develop ED, suggesting that both genetic and environmental factors may play a role. Conversely, a small subset of otherwise healthy children who are not obese may manifest abnormal endothelial function, and such functional phenotype may be determined by genetic variance in endothelial function-related genes.
ED, an early risk marker of cardiovascular disease, refers to a loss of normal homeostatic function in the blood vessels and is characterized by altered vasodilatory and vasoconstrictive functions and inflammatory activity . ED is involved in the development of vascular complications related to dyslipidemia, and cardiovascular disease such as hypertension, coronary artery disease and chronic heart failure [10, 11]. A major contributor of ED involves reductions in the amount of bioavailable nitric oxide (NO) in the vasculature . Indeed, NO is an important factor in endothelial functional homeostasis, and also inhibits platelet aggregation, leukocyte adhesion, smooth muscle cell migration and proliferation , as well as oxidation of atherogenic low-density lipoprotein . NO is synthesized from the amino acid L-arginine in endothelial cells, as well as many other cell types by three nitric oxide synthase (NOS) isoforms; neuronal (nNOS, NOS1), inducible (iNOS, NOS2), and endothelial (eNOS, NOS3) [15–18]. NOS1 and NOS3 genes are constitutively expressed resulting in a low basal synthesis of NO [19, 20], whereas induction of NOS2 expression is regulated by transcription factors such as nuclear factor κB (NF-KB) . Several studies have examined the possibility that genetic variants in the genes encoding these enzymes could influence their expression and functional activity, potentially altering the predisposition to cardiovascular disease [22, 23]. Accordingly, single nucleotide polymorphisms (SNPs) have been identified in NOS genes, and their association with coronary artery disease, hypertension, and diabetes has been explored [24–26]. To the best of our knowledge, the potential associations between NOS polymorphisms and ED are currently unexplored in children.
The endothelin (EDN) system consists of three endothelin isoforms (EDN1, EDN2 and EDN3), and two receptors (endothelin type A and type B) linked to multiple signaling pathways . Endothelins (EDN1, EDN2 and EDN3) are 21-amino acid peptides that exert their effects through their cognate receptors EDNRA and EDNRB with different degrees of binding affinity. Various SNPs have been identified on EDN genes and EDNRA and EDNRB receptors genes, and have been associated with susceptibility and prognosis of diseases such as heart failure, dilated cardiomyopathy, diabetic retinopathy, and atherosclerosis [28–32]. Again, we are unaware of any published studies on the potential associations between ED and endothelin-related gene variants in children.
Based on aforementioned considerations, we aimed to identify single nucleotide polymorphisms (SNPs) in NOS family (3 isoforms), and EDN family (3 isoforms and their 2 cognate receptors) in order to identify potential associations of these SNPs with ED in children.
The study was approved by the University of Louisville Human Research Committee (Louisville, KY, USA), and informed consent was obtained from the legal caregiver of each participant. Consecutive healthy pre-pubertal children (ages 5–10 years) were recruited from the community to investigate endothelial function.
Children were weighed using the InBody 320 scale (Biospace; Cerritos, CA), and height (to 0.1 cm) was measured using a stadiometer (Holtain, Crosswell, UK). The BMI was calculated and the BMI z-score was computed using US Centers for Disease Control and Prevention 2000 growth standards (http://www.cdc.gov/growthcharts/) and online software (http://wwwn.cdc.gov/epiinfo/). A BMI z-score > 1.65 (> 95th percentile) was considered as fulfilling obesity criteria.
All children had arterial blood pressure measured noninvasively using an automated mercury sphygmomanometer (Welch Allyn; Skaneateles Falls, New York) at the brachial artery, using the appropriate cuff size on the nondominant arm . Systolic BP and diastolic BP indices were calculated by dividing the average systolic and diastolic pressure by the respective 95th percentile for BP using National Heart, Lung and Blood Institute guidelines http://www.nhlbi.nih.gov/guidelines/hypertension/child_tbl.htm), computed for age, sex, and height. Hypertension was defined when the SBPi or DBPi was > 1.
Endothelial function tests
Endothelial function was assessed using a modified hyperemic test after cuff-induced occlusion of the radial and ulnar arteries by placing the cuff over the wrist. All testing was performed in the morning after an overnight fast. A laser Doppler sensor (Periflux 5000 System integrated with the PF 5050 Pressure Unit, Perimed, Järfälla, Sweden) was applied over the volar aspect of the hand at the 1st finger distal metacarpal surface and the hand was gently immobilized. This site was chosen as an area in order to minimize the effects of motion artifact, and was also found to have a density of skin capillary blood flow that was of appropriate magnitude for detection. Children lay supine with the head of the bed elevated 45°. Once cutaneous blood flow over the area became stable, the pressure within an inflatable cuff placed at the forearm and connected to a computer- controlled manometer was raised to 200 mmHg for 60 seconds during which blood flow was reduced to undetectable levels. An occlusion time of 60 seconds was chosen in order to minimize discomfort for the child and thus prevent motion and invalidation of the test. Furthermore, we used a computer-controlled pressure release to allow for consistent deflation times. The cuff was rapidly deflated, and the laser Doppler measured hyperemic responses. The maneuver was repeated twice within 10 min with at least 2 min separating both trials to ensure a return to baseline perfusion. The average of both maneuvers was then computed for subsequent analyses, and demonstrated excellent reproducibility . Commercially available software (Perimed, Järfälla, Sweden) allowed for unbiased estimates of the time to peak regional blood flow response post-occlusion release (Tmax), which is considered representative of the hyperemic response and an index of endothelial function . As previously described, we defined abnormal endothelial function as a Tmax cutoff value of ≥45 seconds .
Fasting blood samples were drawn from the subjects and immediately centrifuged in a cooled centrifuge and plasma was frozen at −80°C until assay, usually within <2 weeks from collection. Plasma insulin levels were measured using a commercially available radioimmunoassay kit (Coat-A Count Insulin; Diagnostic Products Inc). This method has a detection level of 1.2 μIU/mL and exhibits linear behavior up to 350 μIU/mL, with intra-assay and inter-assay coefficients of variability of 3.1% and 4.9%, respectively. Plasma glucose level was measured using a commercial kit based on the hexokinase-glucose-6-phosphate dehydrogenase method (Flex Reagent Cartridges; Dade Behring, Newark, DE). Insulin resistance was assessed using the homeostasis model assessment (HOMA) equation (fasting insulin × fasting glucose/22.5) . Plasma hsCRP levels were measured within 2–3 hours after collection using the Flex reagent Cartridge (Date Behring, Newark, DE). This method has a detection level of 0.05 mg/dl, and exhibits linear behavior up to 255 mg/dl, with intra-assay and inter-assay coefficients of variability of 9% and 18% respectively for hsCRP. Serum lipids including total cholesterol, high-density lipoprotein (HDL) cholesterol, calculated low-density lipoprotein cholesterol (LDL), and triglycerides (TG) were also assessed using Flex Reagent Cartridges (Dade Behring). To ensure consistency and to prevent protein degradation, particular care was taken to standardize all steps of plasma sample processing and to minimize thawing more than once for each aliquot.
Peripheral blood samples were collected in vacutainer tubes containing EDTA (Becton Dickinson, Franklin Lakes, NJ, USA). Genomic DNAs of children were extracted using QIAmp DNA blood kit (Qiagen, Valencia, CA, USA). The concentration and quality of the DNA were determined using a ND-1000 Spectrophotometer (Nanodrop Technologies, Wilmington, DE, USA). The purity of the DNA were determined by calculating the ratio of absorbance at 260/280 nm, and all DNA samples had a ratio of 1.8–1.9. The precise length of genomic DNA was determined by gel electrophoresis using 1% agarose gel. All the purified samples were stored at −80°C until further analyses.
Custom cardiovascular gene SNP array
The IBC array was developed using SNP and linkage disequilibrium information from the HapMap as well as resequencing data from Seattle SNPs and National Institute of Environmental Health Sciences (NIEHS) SNPs . Briefly, the IBC array contains about 50,000 SNPS from genetic diversity across approximately 2100 genes related to cardiovascular, inflammatory, hemostasis/coagulation, and metabolic phenotypes and pathways. Among those genes, we selected the NOS genes which include NOS1 (209 SNPs), NOS2 (122 SNPs) and NOS3 (50 SNPs). Furthermore we selected the EDN and EDN receptor genes family which includes EDN1 (43 SNPs), EDN2 (48 SNPs), EDN3 (14SNPs), EDNRA (27 SNPs) and EDNRB (23 SNPs). SNPs were clustered into genotypes with the Illumina Beadstudio software and subjected to quality-control filters at the sample and SNP levels separately within each cohort. Samples were excluded for individual call rates <90%, gender mismatch, and duplicate discordance. SNPs were removed for call rates <95% or Hardy-Weinberg Equilibrium p < 10−7 in controls from each cohort (regardless of ethnicity). Because of the low-frequency SNPs included in the design and the aim to capture low-frequency variants of large effect across the large dataset, we filtered only on minor allele frequency (MAF) < 0.005.
Validation of SNP array using real-time PCR
SNPs that showed statistically significant differences between children with ED and without ED were further validated using TaqMan real-time PCR (Applied Biosystems, Inc.). Two fluorogenic minor groove binder probes were used for each locus using the dyes 6-carboxyfluorescein (FAM; excitation, 494 nm) and VIC (excitation, 538 nm) which are easily differentiated in RT-PCR system. The reaction mixture (25 μl total volume per single well reaction) containing 12.5 μl of TaqMan 2X universal master mix (Applied Biosystems, CA), 1.25 μl of each primer, 10.25 μl of RNase- and DNase-free water (Ambion, Austin, TX), and 1 μl of DNA template. Each sample was duplicated in two wells of a 96 well-plate (Applied Biosystems, CA). DNase-free water used as negative control was included in each assay run. Protocol consisted of 2 min at 50°C, 10 min at 95°C, and 42 cycles of 95°C for 15 sec and 60°C for 1 min. Initially, the SNP assay was set up using SDS, version 2.1, software (Applied Biosystems, CA) as an absolute quantification assay, but after assay completion the plate was read using the allelic discrimination settings. Post-assay analysis was performed using the SDS software. The goodness-of-fit test for Hardy–Weinberg equilibrium (HWE) was performed after analysis using the equation p2 + q2 + 2pq = 1, where p and q represent the wild-type and variant allele of a gene.
Total RNA isolation
Fasting peripheral blood samples were drawn from children within the first hour after awakening and collected in PAXgene Blood RNA tubes (Becton Dickinson, UK). Total RNA was isolated using PAXgene Blood RNA Kit and treated with DNase I (QIAGEN, CA), according to the manufacturer’s protocol. The RNA quantity and integrity were determined using a Nanodrop Spectrophotometer and Agilent 2100 Bioanalyzer Nano 6000 LabChip assay (Agilent Technologies).
Gene expression using qRT-PCR
Quantitative real time RT-PCR (qRT-PCR) analyses were performed using ABI 7500 (Applied Biosystems, Foster City, CA). cDNA synthesis was performed using a High-Capacity cDNA Archive Kit (Applied Biosystems, Foster City, CA). Five hundred nanograms (500 ng) of total RNA from NEF and ED samples were used to generate cDNA templates for RT-PCR. The TaqMan® Master Mix Reagent Kit (Applied Biosystems, Foster City, CA) was used to amplify and quantify each transcript of interest in 25 μl reactions. Various negative controls were included in the PCR reaction to ensure specific amplification. Triplicate PCR reactions were performed in 96-well plates for each gene in parallel with the 18S rRNA. The steps involved in the reaction program included: the initial step of 2 minutes at 50°C; denaturation at 95°C for 10 min, followed by 45 thermal cycles of denaturation (15 seconds at 95°C) and elongation (1 min at 60°C). The expression values were obtained from the cycle number (Ct value) using the Biosystems analysis software. The threshold cycle (CT) values were averaged from each reaction, and each gene was normalized to the 18S rRNA level. All the genes of interest and 18S rRNA were performed in triplicates to determine the Ct-diff. These Ct values were averaged and the difference between the 18S Ct (Avg) and the gene of interest Ct (Avg) was calculated (Ct-diff). The relative expression of the gene of interest was analyzed using the 2-ΔΔCT method . Quantitative results are expressed as the mean ± standard deviation (SD). Statistical significance was evaluated by the Student’s t-test.
All analyses were conducted using SPSS software (version 19.0; SPPS Inc., Chicago, Ill.), and data are presented as mean ± SD. The association analysis was assessed by using Pearson’s chi-square test implemented in SPSS. A P-value < 0.05 was considered statistically significant for all analyses. Odds ratio and 95% confidence interval were calculated for the minor allele of each SNP. We performed these analyses under the a priori assumption that since candidate-gene studies lack power to detect weak genetic risk effects of common variants, any findings from the current analyses will require cautious interpretation. Indeed, to achieve a power of >80% in the detection of a modest genetic risk (e.g., odds ratio = 1.2) for any SNP of interest with a known prevalence of 10% in the population, a sample size of more than 10,000 subjects would be needed [39, 40]. The Haploview version 4.2 software (http://www.borad.mit.edu/mpg/haploview) was used to analyze the linkage disequilibrium structure, calculating D’ to define haplotype block  and to estimate haplotype frequencies. Additionally, pair-wise linkage disequilibrium (LD) among the SNPs was examined using Lewontin’s standardized coefficient D’ and LD coefficient r2, and haplotype blocks were defined according to the method of Garbriel et al.  in Haploview 4.2 with default settings. Haplotypes within these blocks were estimated using the estimation of maximization algorithm .
Of a potential total of 850 subjects, > 600 subjects were recruited from the community, and of these, 122 children were randomly selected and their endothelial function was tested (Figure 1). Approximately 245 subjects were excluded from the study due to chronic medical conditions such as Down syndrome, craniofacial or known genetic syndromes, a known episode of infection in the eight weeks preceding the sleep study, asthma or allergies receiving specific therapy (desensitization, leukotriene inhibitors, steroids (topical or systemic). Eighty four children were found to have NEF and 38 had evidence of ED as defined by their individual Tmax values. The demographic characteristics of these subjects are shown in Table 1. The BMI- z score and the proportion of obese children were significantly higher (P-value ≤ 0.0002) in children with ED. However, there were no differences in either systolic or diastolic blood pressures between both groups. The metabolic data for the 2 groups are shown in Table 1. Only serum triglyceride levels (TG) were significantly higher in children with ED (P-value ≤ 0.04), with both groups exhibiting similar fasting serum TC, HDL and LDL cholesterol. Similarly, fasting glucose levels were similar, but plasma insulin concentrations were higher in ED children (P-value ≤0.05). Serum levels of high sensitivity C-reactive protein (hsCRP) were also similar (Table 1).
As a preamble, all the SNPs tested in this study were in Hardy-Weinberg equilibrium. From a total of 381 SNPs assayed for the 3 NOS genes, 11 SNPs in NOS1 exhibited statistically significant differences in their frequencies among children with and without ED both in the univariate analysis and after correction for multiple comparisons (Table 2). The P-value in Table 2 represents the significance of the allele frequency for each SNP. In addition, the odds ratio (OR) and the confidence intervals of 95% (CI) for all SNPs are shown in Table 2. For example, in addition, to the P-value for the allelic frequency, if a SNP presented in Table 2 is higher than 1.0 that indicates this SNPs is significantly associated with ED in the children. Furthermore, all SNPs in Table 2 have a higher OR and CI than 1.0 with the exception of one SNP (rs9658255), indicating that these SNPs are most likely to be associated with ED. No differences emerged for either NOS2 or NOS3 SNPs.
Linkage disequilibrium (LD) analysis of the 11 SNPs in the NOS1 gene was assessed for both NEF and ED subjects. In NEF subjects, two haplotype blocks emerged, and are outlined in black triangular regions in Figure 2 (Panel A). In ED subjects, the haplotype showed the presence of 2 blocks (Figure 2, Panel B). The haplotype of these blocks and their frequencies in NEF and ED are shown in Figure 3, Panels A and B, respectively. Taken together, the patterns of LD and haplotype frequencies differed between NEF and ED, suggesting that some of these SNPs may contribute to ED risk. Therefore, allelic data were analyzed in conjunction with Tmax values among ED and NEF subjects for the significant NOS1 SNPs as shown in Table 3. Of the 11 significant SNPs of the NOS1 gene, there were only 3 SNPs (rs6490121, rs483589, rs1879417) that showed significantly different Tmax based on the genotype. We found that the A allele in rs6490121 (A/G) and rs483589 (A/G) had significantly higher Tmax than the G allele (Figure 3, Panels A and B), while the C allele in rs1879417 (C/T) had significantly higher Tmax values compared to the T allele (Figure 4, Panel C).
Endothelin related genes
From a total of 155 available SNPs, there were only 2 SNPs in the EDN1 gene whose frequencies differed among children with and without ED. The allele and genotype frequencies of EDN1 SNPs are shown in Table 4. These 2 EDN1 SNPs were analyzed in conjunction with corresponding Tmax values among ED and NEF subjects. We found that both of these SNPs in EDN1 (rs1476046, rs4714384) showed significant differences in Tmax based on genotype (Table 5). The G allele in rs1476046 (A/G) had a significantly higher Tmax than the A allele (Figure 4, Panel D), and the T allele in rs4714384 (C/T) had a significantly higher Tmax than the C allele (Figure 4, Panel E).
qRT-PCR SNP assay validation
To confirm the SNP array-based findings, we performed qRT-PCR using TaqMan assays with specific primers for some of the SNPs of interest, i.e., those with significantly different frequencies in their allelic genotype frequencies (Additional file 1: Table S1), as well as some of those with significantly different genotype-phenotype allelic frequencies (Additional file 1: Table S1). The concordance findings between array and real-time PCR for allele and genotype frequencies of NOS1 and EDN1 SNPs in children with and without ED are shown in Additional file 2: Table S2. The concordance for both SNPs in the NOS1 gene (rs3825102, rs483589) was 53.21%, and 61.47%, respectively and 49% for EDN1 SNPS (rs1476046, rs4714384). The summary of the validated SNPs data for both NOS1 and EDN1 genes between SNPs array and qRT-PCR is shown in Additional file 3: Table S3.
qRT-PCR gene expression
Next we asked how the SNPs that putatively affected endothelial function, specifically, their association with gene expression. We selected a total of 20 matched subjects, 10 with endothelial dysfunction (ED), and 10 with normal endothelial function (NEF), and performed gene expression analysis using qRT-PCR. A significant increase in EDN1 gene expression in the peripheral blood of children with ED emerged compared to NEF as illustrated in Figure 5 (P-value 0.0004).
In this study, we report on the associations between the frequency of the 3 NOS genes, and 3 genes of the EDN family, along with its two cognate receptors in a carefully phenotyped cohort of pre-pubertal otherwise healthy community children with and without ED. Not surprisingly, an increased risk of obesity and of increased serum TG levels was present in children with ED. The latter may underlie previously reported lower nitric oxide levels in the context of elevated TG . In addition, the relative frequencies of only a selected few NOS1 and EDN1 gene polymorphisms differed among children with ED, and assessment of Tmax based on allelic variants further confirmed a higher risk of increased Tmax among those children harboring the polymorphic allele.
Before we discuss the potential significance of our findings, some of the strengths and limitations of the study should be mentioned. First, the use of the laser Doppler technique for assessment of vascular responses following cuff-induced arterial occlusion not only allows reproducible determinations of the kinetics of post-ischemic reperfusion, but also provides an accurate reporter of NO-mediated physiological recruitment of the microvasculature [34, 45]. The selection of the 60-seconds occlusion time, rather than the more extensively used 5 minutes in adults, was necessary based on a vast set of preliminary experiments aimed at preventing motion-induced artifacts, and thus preserving the validity of the tests. In addition, all Tmax assessments were conducted in the fasting state at the same time of the day to minimize potential confounders introduced by differences in meal content and timing relative to the testing, and also to ensure that circadian variation in endothelial function would not play a role. An additional strength of the study was the inclusion of a priori healthy children from the community. In addition, we excluded children with a variety of diagnoses that can be associated with endothelial dysfunction . Two important limitations of this study include the relatively small size of the cohort of children studied which could hamper statistical power, and the absence of endothelin plasma level measurements. However, these data may be incorporated into multicenter-based meta-analyses that will undoubtedly enhance the validity of the present putative findings. Finally, in the context of our qRT-pCR validation efforts using TaqMan SNP assays, we found substantial discrepancies between the two methods (i.e., SNP arrays and qRT-PCR), even if the significance of the findings remained unaltered. The concordance between the array data and qRT-PCR revolved around only 50-60%. This was surprising considering that the concordance between gene expression arrays and qRT-PCR is remarkably high. It therefore appears that this may not be case for SNP arrays, and the discordant findings necessarily mandate that future studies should require an additional validation step of the array findings using alternative methodologies, an approach that has not been pursued to date while using the same 50 k SNP array as used herein or other commercially available arrays [47–50].
Obesity is a prominent risk factor for the development of cardiovascular disease, diabetes mellitus, dyslipidemia, hypertension, and their constellation as the metabolic syndrome in both adults and children [51–55]. We and others have previously shown that obese children are at higher risk for ED and that both the vasodilatory and vasoconstrictive responses are important. Interestingly, in a recent study Tounian and colleagues reported their analyses of a cohort of 232 severely obese children in whom 12 gene polymorphisms representing selected variants previously associated with vascular ED in adults were examined . These investigators did not find any significant association between flow-mediated vasodilation and any of the 12 SNPs studied. Based on our current findings, we propose that the increased risk for ED in obesity may be explained at least in part by the presence of gene variants in NOS and EDN.
Several studies have investigated the NOS genes family SNPs in relation to different diseases, primarily involving the cardiovascular system [57–60]. In the present study, we found 11 SNPs in the NOS1 gene showing significant differences between children with and without ED. In addition, out of these 11 SNPs, we found 3 SNPs (rs6490121, rs483589, and rs1879417) that were significantly associated with altered Tmax values in children, suggesting that these SNPs carry biological significance. However, no significant differences or associations among SNPs in either NOS2 or NOS3 genes were found.
Similar to studies on NOS gene SNPs, several studies suggest that some of the EDN1 SNPs may be associated with hypertension. For example, EDN1 SNP (rs5370) Lys198Asn (G/T), interactions between BMI and blood pressure emerged in overweight/obese subjects [61–64]. This non-synonymous SNP has already been associated with heart failure and hypertension [61, 65]. Furthermore, in rheumatoid arthritis patients who carry the T-T haplotype in EDN1 SNPs (rs1800541 and rs5370), increased systolic blood pressure and high EDN1 plasma levels emerged . In this study, we identified 2 significant SNPs (rs1476046, rs4714384) in the EDN1 gene, and both of these SNPs showed significant differences in Tmax values, adding corroborative evidence to the assumption that these SNPs have biological significance. However, no significant differences emerged in the allelic frequencies in EDN2, EDN3, EDNRA, and EDNRB.
In summary, significant associations of polymorphisms in NOS and EDN genes are apparent among young children with ED. Of note, EDN1 gene expression appears to be increased in the context of ED. Accordingly; further functional studies are needed to clarify the role of variants of NOS and EDN genes in the pathophysiology of cardiovascular disease.
DG is supported by National Institutes of Health grants HL-065270 and HL-086662. SC is the recipient of the Royal Golden Jubilee Ph. D. Program Award from the Thailand Research Fund (PHD/0138/2552).
Kelishadi R: Childhood overweight, obesity, and the metabolic syndrome in developing countries. Epidemiol Rev. 2007, 29: 62-76. 10.1093/epirev/mxm003.
Haslam DW, James WP: Obesity. Lancet. 2005, 366 (9492): 1197-1209. 10.1016/S0140-6736(05)67483-1.
Wang Y, Mi J, Shan XY, Wang QJ, Ge KY: Is China facing an obesity epidemic and the consequences? The trends in obesity and chronic disease in China. International journal of obesity. 2007, 31 (1): 177-188. 10.1038/sj.ijo.0803354.
Reilly JJ, Kelly J: Long-term impact of overweight and obesity in childhood and adolescence on morbidity and premature mortality in adulthood: systematic review. International journal of obesity. 2011, 35 (7): 891-898. 10.1038/ijo.2010.222.
Ribeiro MM, Silva AG, Santos NS, Guazzelle I, Matos LN, Trombetta IC, Halpern A, Negrao CE, Villares SM: Diet and exercise training restore blood pressure and vasodilatory responses during physiological maneuvers in obese children. Circulation. 2005, 111 (15): 1915-1923. 10.1161/01.CIR.0000161959.04675.5A.
Aggoun Y, Farpour-Lambert NJ, Marchand LM, Golay E, Maggio AB, Beghetti M: Impaired endothelial and smooth muscle functions and arterial stiffness appear before puberty in obese children and are associated with elevated ambulatory blood pressure. Eur Heart J. 2008, 29 (6): 792-799. 10.1093/eurheartj/ehm633.
Bhattacharjee R, Alotaibi WH, Kheirandish-Gozal L, Capdevila OS, Gozal D: Endothelial dysfunction in obese non-hypertensive children without evidence of sleep disordered breathing. BMC pediatrics. 2010, 10: 8-10.1186/1471-2431-10-8.
Gozal D, Kheirandish-Gozal L: Cardiovascular morbidity in obstructive sleep apnea: oxidative stress, inflammation, and much more. Am J Respir Crit Care Med. 2008, 177 (4): 369-375. 10.1164/rccm.200608-1190PP.
Davignon J, Ganz P: Role of endothelial dysfunction in atherosclerosis. Circulation. 2004, 109 (23 Suppl 1): III27-III32.
Hirata Y, Nagata D, Suzuki E, Nishimatsu H, Suzuki J, Nagai R: Diagnosis and treatment of endothelial dysfunction in cardiovascular disease. International heart journal. 2010, 51 (1): 1-6. 10.1536/ihj.51.1.
Versari D, Daghini E, Virdis A, Ghiadoni L, Taddei S: Endothelial dysfunction as a target for prevention of cardiovascular disease. Diabetes Care. 2009, 32 (Suppl 2): S314-S321.
Hingorani AD: Polymorphisms in endothelial nitric oxide synthase and atherogenesis: John French Lecture 2000. Atherosclerosis. 2001, 154 (3): 521-527. 10.1016/S0021-9150(00)00699-7.
Kugiyama K, Yasue H, Okumura K, Ogawa H, Fujimoto K, Nakao K, Yoshimura M, Motoyama T, Inobe Y, Kawano H: Nitric oxide activity is deficient in spasm arteries of patients with coronary spastic angina. Circulation. 1996, 94 (3): 266-271. 10.1161/01.CIR.94.3.266.
Cooke JP, Dzau VJ: Nitric oxide synthase: role in the genesis of vascular disease. Annu Rev Med. 1997, 48: 489-509. 10.1146/annurev.med.48.1.489.
Lamas S, Marsden PA, Li GK, Tempst P, Michel T: Endothelial nitric oxide synthase: molecular cloning and characterization of a distinct constitutive enzyme isoform. Proc Natl Acad Sci USA. 1992, 89 (14): 6348-6352. 10.1073/pnas.89.14.6348.
Bogdan C: Nitric oxide and the immune response. Nature immunology. 2001, 2 (10): 907-916. 10.1038/ni1001-907.
Alderton WK, Cooper CE, Knowles RG: Nitric oxide synthases: structure, function and inhibition. Biochem J. 2001, 357 (Pt 3): 593-615.
Hansel TT, Kharitonov SA, Donnelly LE, Erin EM, Currie MG, Moore WM, Manning PT, Recker DP, Barnes PJ: A selective inhibitor of inducible nitric oxide synthase inhibits exhaled breath nitric oxide in healthy volunteers and asthmatics. FASEB journal : official publication of the Federation of American Societies for Experimental Biology. 2003, 17 (10): 1298-1300.
Cook S, Vollenweider P, Menard B, Egli M, Nicod P, Scherrer U: Increased eNO and pulmonary iNOS expression in eNOS null mice. The European respiratory journal : official journal of the European Society for Clinical Respiratory Physiology. 2003, 21 (5): 770-773. 10.1183/09031936.03.00121203.
Vaughan DJ, Brogan TV, Kerr ME, Deem S, Luchtel DL, Swenson ER: Contributions of nitric oxide synthase isozymes to exhaled nitric oxide and hypoxic pulmonary vasoconstriction in rabbit lungs. American journal of physiology Lung cellular and molecular physiology. 2003, 284 (5): L834-L843.
Silkoff PE, McClean PA, Slutsky AS, Furlott HG, Hoffstein E, Wakita S, Chapman KR, Szalai JP, Zamel N: Marked flow-dependence of exhaled nitric oxide using a new technique to exclude nasal nitric oxide. Am J Respir Crit Care Med. 1997, 155 (1): 260-267. 10.1164/ajrccm.155.1.9001322.
Alvarez R, Gonzalez P, Batalla A, Reguero JR, Iglesias-Cubero G, Hevia S, Cortina A, Merino E, Gonzalez I, Alvarez V: Association between the NOS3 (−786 T/C) and the ACE (I/D) DNA genotypes and early coronary artery disease. Nitric oxide : biology and chemistry/official journal of the Nitric Oxide Society. 2001, 5 (4): 343-348. 10.1006/niox.2001.0351.
Salazar LA, Hirata MH, Giannini SD, Forti N, Diament J, Lima TM, Hirata RD: Seven DNA polymorphisms at the candidate genes of atherosclerosis in Brazilian women with angiographically documented coronary artery disease. Clinica chimica acta; international journal of clinical chemistry. 2000, 300 (1–2): 139-149.
Zanchi A, Moczulski DK, Hanna LS, Wantman M, Warram JH, Krolewski AS: Risk of advanced diabetic nephropathy in type 1 diabetes is associated with endothelial nitric oxide synthase gene polymorphism. Kidney Int. 2000, 57 (2): 405-413. 10.1046/j.1523-1755.2000.00860.x.
Wang XL, Wang J: Endothelial nitric oxide synthase gene sequence variations and vascular disease. Mol Genet Metab. 2000, 70 (4): 241-251. 10.1006/mgme.2000.3033.
Zintzaras E, Kitsios G, Stefanidis I: Endothelial NO synthase gene polymorphisms and hypertension: a meta-analysis. Hypertension. 2006, 48 (4): 700-710. 10.1161/01.HYP.0000238124.91161.02.
Brunner F, Bras-Silva C, Cerdeira AS, Leite-Moreira AF: Cardiovascular endothelins: essential regulators of cardiovascular homeostasis. Pharmacology & therapeutics. 2006, 111 (2): 508-531. 10.1016/j.pharmthera.2005.11.001.
Colombo MG, Ciofini E, Paradossi U, Bevilacqua S, Biagini A: ET-1 Lys198Asn and ET(A) receptor H323H polymorphisms in heart failure. A case–control study. Cardiology. 2006, 105 (4): 246-252. 10.1159/000092374.
Herrmann S, Schmidt-Petersen K, Pfeifer J, Perrot A, Bit-Avragim N, Eichhorn C, Dietz R, Kreutz R, Paul M, Osterziel KJ: A polymorphism in the endothelin-A receptor gene predicts survival in patients with idiopathic dilated cardiomyopathy. Eur Heart J. 2001, 22 (20): 1948-1953. 10.1053/euhj.2001.2626.
Telgmann R, Harb BA, Ozcelik C, Perrot A, Schonfelder J, Nonnenmacher A, Brand M, Schmidt-Petersen K, Dietz R, Kreutz R: The G-231A polymorphism in the endothelin-A receptor gene is associated with lower aortic pressure in patients with dilated cardiomyopathy. American journal of hypertension. 2007, 20 (1): 32-37. 10.1016/j.amjhyper.2006.06.016.
Li H, Louey JW, Choy KW, Liu DT, Chan WM, Chan YM, Fung NS, Fan BJ, Baum L, Chan JC: EDN1 Lys198Asn is associated with diabetic retinopathy in type 2 diabetes. Mol Vis. 2008, 14: 1698-1704.
Yasuda H, Kamide K, Takiuchi S, Matayoshi T, Hanada H, Kada A, Yang J, Miwa Y, Yoshii M, Horio T: Association of single nucleotide polymorphisms in endothelin family genes with the progression of atherosclerosis in patients with essential hypertension. J Hum Hypertens. 2007, 21 (11): 883-892. 10.1038/sj.jhh.1002234.
National High Blood Pressure Education Program Working Group on High Blood Pressure in C, Adolescents: The fourth report on the diagnosis, evaluation, and treatment of high blood pressure in children and adolescents. Pediatrics. 2004, 114 (2 Suppl 4th Report): 555-576.
Bhattacharjee R, Kim J, Alotaibi WH, Kheirandish-Gozal L, Capdevila OS, Gozal D: Endothelial dysfunction in children without hypertension: potential contributions of obesity and obstructive sleep apnea. Chest. 2012, 141 (3): 682-691. 10.1378/chest.11-1777.
Wahlberg E, Olofsson P, Swendenborg J, Fagrell B: Changes in postocclusive reactive hyperaemic values as measured with laser Doppler fluxmetry after infrainguinal arterial reconstructions. Eur J Vasc Endovasc Surg. 1995, 9 (2): 197-203. 10.1016/S1078-5884(05)80090-5.
Matthews DR, Hosker JP, Rudenski AS, Naylor BA, Treacher DF, Turner RC: Homeostasis model assessment: insulin resistance and beta-cell function from fasting plasma glucose and insulin concentrations in man. Diabetologia. 1985, 28 (7): 412-419. 10.1007/BF00280883.
Keating BJ, Tischfield S, Murray SS, Bhangale T, Price TS, Glessner JT, Galver L, Barrett JC, Grant SF, Farlow DN: Concept, design and implementation of a cardiovascular gene-centric 50 k SNP array for large-scale genomic association studies. PLoS One. 2008, 3 (10): e3583-10.1371/journal.pone.0003583.
Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) Method. Methods. 2001, 25 (4): 402-408. 10.1006/meth.2001.1262.
Zintzaras E, Lau J: Synthesis of genetic association studies for pertinent gene-disease associations requires appropriate methodological and statistical approaches. JCE. 2008, 61 (7): 634-645.
Zintzaras E, Lau J: Trends in meta-analysis of genetic association studies. J Hum Genet. 2008, 53 (1): 1-9. 10.1007/s10038-007-0223-5.
Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, Blumenstiel B, Higgins J, DeFelice M, Lochner A, Faggart M: The structure of haplotype blocks in the human genome. Science. 2002, 296 (5576): 2225-2229. 10.1126/science.1069424.
Lewontin RC: On measures of gametic disequilibrium. Genetics. 1988, 120 (3): 849-852.
Schaid DJ, Rowland CM, Tines DE, Jacobson RM, Poland GA: Score tests for association between traits and haplotypes when linkage phase is ambiguous. Am J Hum Genet. 2002, 70 (2): 425-434. 10.1086/338688.
Sarwar N, Danesh J, Eiriksdottir G, Sigurdsson G, Wareham N, Bingham S, Boekholdt SM, Khaw KT, Gudnason V: Triglycerides and the risk of coronary heart disease: 10,158 incident cases among 262,525 participants in 29 Western prospective studies. Circulation. 2007, 115 (4): 450-458. 10.1161/CIRCULATIONAHA.106.637793.
Medow MS, Taneja I, Stewart JM: Cyclooxygenase and nitric oxide synthase dependence of cutaneous reactive hyperemia in humans. American journal of physiology Heart and circulatory physiology. 2007, 293 (1): H425-H432. 10.1152/ajpheart.01217.2006.
Valle Jimenez M, Estepa RM, Camacho RM, Estrada RC, Luna FG, Guitarte FB: Endothelial dysfunction is related to insulin resistance and inflammatory biomarker levels in obese prepubertal children. European journal of endocrinology/European Federation of Endocrine Societies. 2007, 156 (4): 497-502. 10.1530/EJE-06-0662.
Talmud PJ, Cooper JA, Gaunt T, Holmes MV, Shah S, Palmen J, Drenos F, Shah T, Kumari M, Kivimaki M: Variants of ADRA2A are associated with fasting glucose, blood pressure, body mass index and type 2 diabetes risk: meta-analysis of four prospective studies. Diabetologia. 2011, 54 (7): 1710-1719. 10.1007/s00125-011-2108-6.
Fox ER, Young JH, Li Y, Dreisbach AW, Keating BJ, Musani SK, Liu K, Morrison AC, Ganesh S, Kutlar A: Association of genetic variation with systolic and diastolic blood pressure among African Americans: the Candidate Gene Association Resource study. Hum Mol Genet. 2011, 20 (11): 2273-2284. 10.1093/hmg/ddr092.
Ronald J, Rajagopalan R, Cerrato F, Nord AS, Hatsukami T, Kohler T, Marcovina S, Heagerty P, Jarvik GP: Genetic variation in LPAL2, LPA, and PLG predicts plasma lipoprotein(a) level and carotid artery disease risk. Stroke; a journal of cerebral circulation. 2011, 42 (1): 2-9. 10.1161/STROKEAHA.110.591230.
Saleheen D, Alexander M, Rasheed A, Wormser D, Soranzo N, Hammond N, Butterworth A, Zaidi M, Haycock P, Bumpstead S: Association of the 9p21.3 locus with risk of first-ever myocardial infarction in Pakistanis: case–control study in South Asia and updated meta-analysis of Europeans. Arterioscler Thromb Vasc Biol. 2010, 30 (7): 1467-1473. 10.1161/ATVBAHA.109.197210.
Van Gaal LF, Mertens IL, De Block CE: Mechanisms linking obesity with cardiovascular disease. Nature. 2006, 444 (7121): 875-880. 10.1038/nature05487.
Freedman DS, Dietz WH, Srinivasan SR, Berenson GS: The relation of overweight to cardiovascular risk factors among children and adolescents: the Bogalusa Heart Study. Pediatrics. 1999, 103 (6 Pt 1): 1175-1182.
Morrison JA, Laskarzewski PM, Rauh JL, Brookman R, Mellies M, Frazer M, Khoury P, deGroot I, Kelly K, Glueck CJ: Lipids, lipoproteins, and sexual maturation during adolescence: the Princeton maturation study. Metabolism. 1979, 28 (6): 641-649. 10.1016/0026-0495(79)90017-9.
Pinhas-Hamiel O, Dolan LM, Daniels SR, Standiford D, Khoury PR, Zeitler P: Increased incidence of non-insulin-dependent diabetes mellitus among adolescents. J Pediatr. 1996, 128 (5 Pt 1): 608-615.
Sinha R, Fisch G, Teague B, Tamborlane WV, Banyas B, Allen K, Savoye M, Rieger V, Taksali S, Barbetta G: Prevalence of impaired glucose tolerance among children and adolescents with marked obesity. N Engl J Med. 2002, 346 (11): 802-810. 10.1056/NEJMoa012578.
Tounian A, Aggoun Y, Lacorte JM, Dubern B, Clement K, Bonnet D, Tounian P: Influence of polymorphisms in candidate genes on early vascular alterations in obese children. Archives of cardiovascular diseases. 2010, 103 (1): 10-18. 10.1016/j.acvd.2009.10.001.
Lillie EO, Mahata M, Khandrika S, Rao F, Bundey RA, Wen G, Chen Y, Taupenot L, Smith DW, Mahata SK: Heredity of endothelin secretion: human twin studies reveal the influence of polymorphism at the chromogranin A locus, a novel determinant of endothelial function. Circulation. 2007, 115 (17): 2282-2291. 10.1161/CIRCULATIONAHA.106.648345.
Gormley K, Bevan S, Hassan A, Markus HS: Polymorphisms in genes of the endothelin system and cerebral small-vessel disease. Stroke; a journal of cerebral circulation. 2005, 36 (8): 1656-1660. 10.1161/01.STR.0000173173.38289.69.
Dell’Omo G, Penno G, Pucci L, Fotino C, Lucchesi D, Del Prato S, Pedrinelli R: Lack of association between endothelial nitric oxide synthase gene polymorphisms, microalbuminuria and endothelial dysfunction in hypertensive men. J Hypertens. 2007, 25 (7): 1389-1395. 10.1097/HJH.0b013e3281268548.
Song XY, Lee SY, Ma RC, So WY, Cai JH, Tam C, Lam V, Ying W, Ng MC, Chan JC: Phenotype-genotype interactions on renal function in type 2 diabetes: an analysis using structural equation modelling. Diabetologia. 2009, 52 (8): 1543-1553. 10.1007/s00125-009-1400-1.
Tiret L, Poirier O, Hallet V, McDonagh TA, Morrison C, McMurray JJ, Dargie HJ, Arveiler D, Ruidavets JB, Luc G: The Lys198Asn polymorphism in the endothelin-1 gene is associated with blood pressure in overweight people. Hypertension. 1999, 33 (5): 1169-1174. 10.1161/01.HYP.33.5.1169.
Asai T, Ohkubo T, Katsuya T, Higaki J, Fu Y, Fukuda M, Hozawa A, Matsubara M, Kitaoka H, Tsuji I: Endothelin-1 gene variant associates with blood pressure in obese Japanese subjects: the Ohasama Study. Hypertension. 2001, 38 (6): 1321-1324. 10.1161/hy1101.095333.
Jin JJ, Nakura J, Wu Z, Yamamoto M, Abe M, Tabara Y, Yamamoto Y, Igase M, Kohara K, Miki T: Association of endothelin-1 gene variant with hypertension. Hypertension. 2003, 41 (1): 163-167. 10.1161/01.HYP.0000043680.75107.CF.
Barath A, Endreffy E, Bereczki C, Gellen B, Szucs B, Nemeth I, Turi S: Endothelin-1 gene and endothelial nitric oxide synthase gene polymorphisms in adolescents with juvenile and obesity-associated hypertension. Acta Physiologica. 2007, 94 (1–2): 49-66.
Pare G, Serre D, Brisson D, Anand SS, Montpetit A, Tremblay G, Engert JC, Hudson TJ, Gaudet D: Genetic analysis of 103 candidate genes for coronary artery disease and associated phenotypes in a founder population reveals a new association between endothelin-1 and high-density lipoprotein cholesterol. Am J Hum Genet. 2007, 80 (4): 673-682. 10.1086/513286.
Panoulas VF, Douglas KM, Smith JP, Taffe P, Stavropoulos-Kalinoglou A, Toms TE, Elisaf MS, Nightingale P, Kitas GD: Polymorphisms of the endothelin-1 gene associate with hypertension in patients with rheumatoid arthritis. Endothelium : journal of endothelial cell research. 2008, 15 (4): 203-212.
The authors declare that they have no competing interests.
SC performed data analysis, and qRT-PCR validation, and help in writing in method section, DG provided the conceptual design of the project, participated in the data analysis and editing final version of the manuscript, L K-G participated in recruited subjects, data analysis and sleep studies, R B participated in recruited subjects, AAK participated in genotypic-phenotypic polymorphisms for pairwise linkage disequilibrium structure and haplotype frequencies using Haploview software, YW reviewed data, BK and HH contributed to genomic polymorphism analyses, WS participated in general discussion and review data, and AK carried data analysis, overall the project and SNPs analysis, writing and editing manuscript. All authors read and approved the final manuscript.