Abstract
Aging leads to a number of disadvantageous changes in the cardiovascular system. Deterioration of vascular homoeostasis with increase in oxidative stress, chronic low-grade inflammation, and impaired nitric oxide bioavailability results in endothelial dysfunction, increased vascular stiffness, and compromised arterial-ventricular interactions. A chronic dietary supplementation with the micronutrient nitrate has been demonstrated to improve vascular function. Healthy dietary patterns may regulate gene expression profiles. However, the mechanisms are incompletely understood. The changes that occur at the gene expression level and transcriptional profile following a nutritional modification with nitrate have not been elucidated. To determine the changes of the vascular transcriptome, we conducted gene expression microarray experiments on aortas of old mice, which were treated with dietary nitrate. Our results highlight differentially expressed genes overrepresented in gene ontology categories. Molecular interaction and reaction pathways involved in the calcium-signaling pathway and the detoxification system were identified. Our results provide novel insight to an altered gene-expression profile in old mice following nitrate supplementation. This supports the general notion of nutritional approaches to modulate age-related changes of vascular functions and its detrimental consequences.
1. Introduction
Aging is considered to be the most important nonmodifiable cardiovascular risk factor, leading to increased cardiovascular morbidity and mortality [1]. It deteriorates vascular integrity and entails a number of disadvantageous structural and functional changes on the arterial system [2]. Structural alterations include increased luminal diameter, vascular wall remodeling with thickened intima and media, and altered extracellular matrix content [2]. Functional changes are characterized by increased aortic stiffness, endothelial dysfunction, and subsequent systolic hypertension. These alterations are regulated by signaling pathways activated by endogenous and exogenous factors [3–5]. On a cellular level, decreased protein synthesis, dysfunctional mitochondria, oxidative stress, impaired antioxidant defense mechanisms, disturbed calcium handling, and increased DNA and protein oxidation have been proposed [6–9]. On the genomic level, we recently identified differences in the transcriptional profile and the genes potentially responsible for the vascular aging process, with emphasis on the detoxification system [10].
Dietary interventions have emerged as novel options, counteracting age-related vascular dysfunction [11–14]. A definite composition of healthy diet is still elusive. Nevertheless, it is known, particularly since the DASH (dietary approaches to stop hypertension) trial, that certain dietary patterns, for example, the so-called Mediterranean diet, positively influence blood pressures and cardiovascular outcome [15–18].
Inorganic nitrate is abundant in our everyday diet and especially in leafy greens. The nutritional aspects of nitrate are intriguing since it represents a source for the endogenous generation of nitrite, nitric oxide (NO), and other reactive nitrogen oxides [19–21]. In experimental studies and first clinical observations, dietary inorganic nitrate interventions have emerged as novel nutrient options modulating cardiovascular functions. This relates to regenerative processes in ischemic tissues, enhancement of endothelial function, and limitation of mitochondrial reactive oxygen species in a cardiomyopathy model [22–24]. We have demonstrated that a chronic nitrate supplementation improves vascular stiffness and reduces systolic blood pressures and the level of proinflammatory cytokines in elderly volunteers [25, 26]. However, to comprehend the beneficial effects of a nitrate rich diet on the vasculature the underlying cellular and molecular changes have to be characterized.
Nitrate-derived nitrite and NO exert beneficial effects through second messengers, posttranslational modifications, and regulation of gene expression [27–29]. The genes potentially affected by a nitrate rich diet counteracting the vascular aging process have not been identified so far. We therefore set out to compare the transcriptional profile in aged aortas of old nitrate-treated mice with age-related endothelial dysfunction applying a high-throughput genomic screening approach. Using microarray analyses, identification of genes, and expression differences was simultaneously accomplished. Our study discloses alterations in gene expression following a dietary nitrate supplementation in the aged vascular system.
2. Material and Methods
2.1. Animals
Male C57BL/6 mice were obtained from Janvier (Saint Berthevin, France) and kept one week in the local animal house for acclimatization. Old (20 months) mice were housed in controlled temperature and light conditions with feeding (V1534-000, Sniff, Soest, Germany) and watering ad libitum. A group of old mice received 1 g/L sodium-nitrate (NaNO3, ~150 μmoles according to the measured daily drinking water consumption) in their drinking water for 8 weeks, as described previously [23, 26]. An equal concentration of sodium chloride served as control. For aortic vasoreactivity studies young (6 months) mice served as controls. For the gene-expression experiments different mice from the same supplier were used. The care and all experiments were performed according to the animal welfare regulations of the German local authorities conforming to NIH Guidelines.
2.2. Aortic Vascular Reactivity
Ex vivo endothelial functions were assessed using an aortic ring bioassay, as described [28, 30]. Briefly, three 2-3-mm wide aortic rings were prepared from each mouse, perivascular fat was removed and aortas were suspended in an organ bath containing 10 mL modified Krebs-Henseleit buffer and connected to an isometric force transducer (EMKA Technologies, Paris, France). Organ chambers were aerated continuously with 95% O2 and 5% CO2 to achieve pH of 7.4. Isometric tension was recorded continuously. The vascular rings with intact endothelium and similar dimensions equilibrated to a resting tension of 1.0 g. All rings were preconstricted with KCl (40 mM) and, after washout, phenylephrine (10 μM). Endothelium function was measured as the relaxation response to acetylcholine (from 10−9 to 10−4 M). Endothelium-independent relaxation was determined as the response to sodium nitroprusside (from 10−9 to 10−2 M).
2.3. Chemicals
All reagents were obtained from Sigma-Aldrich (Taufkirchen, Germany) unless indicated otherwise.
2.4. RNA Preparation and Microarray Assays
RNA preparation and microarray assays were conducted as described previously [10]. Briefly, the descending aorta of each mouse was used. The vessel was flushed thoroughly with ice-cold PBS, through the left ventricle of the heart, cleaned of periadventitial fat and connective tissues, snap-frozen in liquid nitrogen, and stored at −80°C. Extraction of total RNA from the tissue was performed (Qiagen RNeasy Mini Kit, Qiagen, Hilden, Germany). Removal of residual genomic DNA was achieved by treatment with DNase I. Total RNA preparations were checked for RNA integrity (Agilent 2100 Bioanalyzer, Agilent Technologies, Palo Alto, CA). All samples showed common high quality RNA Integrity Numbers (RIN 9.0–9.6) and RNA was quantified by photometric Nanodrop measurement.
Synthesis of cDNA and subsequent fluorescent labeling of cRNA was done on six replicates of each group (old and old nitrate treated mice) according to the manufacturers protocol (One-Color Microarray-Based Gene Expression Analysis/Low Input Quick Amp Labeling; Agilent Technologies, Palo Alto, CA). Briefly, 100 ng of total RNA was converted to cDNA, followed by in vitro transcription and incorporation of Cy3-CTP into nascent cRNA. After fragmentation labeled cRNA was hybridized to Agilent SurePrint G3 Mouse GE 8x60K Microarrays for 17 h at 65°C. Quality control parameters of cRNA labeling and hybridization performance were found within the manufacturers specifications. Arrays were scanned as described by the manufacturer. Signal intensities on 20 bit tiff images were calculated by Feature Extraction software (FE, Vers. 10.7.1.1; Agilent Technologies, Palo Alto, CA).
GeneSpring GX software was used to analyze data (Vers. 12.5; Agilent Technologies, Palo Alto, CA). Probe signal intensities were quantile normalized across all samples to reduce interarray variability [31]. Input data preprocessing was concluded by baseline transformation to the median of all samples. Biological replicates were grouped according to their respective experimental condition. A given transcript had to be expressed above background (“detected” by GeneSpring GX) in at least five of six replicates in any one of two, or both conditions to be additional analyzed in pairwise comparisons of conditions. Hierarchical clustering of gene expression patterns identified no obvious outlier samples. Differential gene expression was determined statistically by a moderated -test with genes showing considered to be significantly differentially expressed [32].
The microarray data set as been deposited in the NCBI Gene Expression Omnibus (GEO) public database in compliance with MIAME (Minimum Information About a Microarray Experiment) guidelines (GEO series accession number GSE66572).
2.5. Functional Analysis by Gene Ontology and Gene Pathways
GeneSpring GX and Database for Annotation, Visualization, and Integrated Discovery (DAVID Bioinformatics database) was used to identify functional related categories of genes.
Differentially expressed genes () with an additional fold change > 1.2 were considered relevant and were imported into DAVID Bioinformatics database [33, 34]. The DAVID software analysed significant enrichment of differentially expressed genes within gene ontology (GO) terms and additionally involves assessment of advanced pathway analysis [35, 36]. The background was chosen from the probes represented on the Agilent 8x60K mouse array, which were detected on either one of the three conditions to determine nitrate-associated genes in the aortic tissue. A modified Fisher’s exact test assessed the significance of the association between the observed data and the data in the GO and canonical pathway. Each identified pathway was assigned a significance score, defined by the number of differentially regulated focus genes in the data set. This score was the negative logarithm of the value, indicative of the likelihood that genes are found in a pathway randomly.
2.6. Statistical Analysis
The results of endothelial functions are presented as mean ± SEM. Data were analyzed by 1-way ANOVA and posthoc Bonferroni multiple comparison correction (for all pairwise tests) with GraphPad Prism 5 (GraphPad, California) software to compare differences between groups. A value of <0.05 was considered significant. Agilent microarray probe sets were analyzed as described above.
3. Results and Discussion
3.1. Dietary Nitrate Reverses Age-Related Endothelial Dysfunction
We studied old C57BL/6 mice and characterized age-related endothelial functions and the reversal through dietary nitrate. Using our model of aged mice, we were able to demonstrate an age-related vascular phenotype and associated alterations in vascular gene expression pattern [10]. Endothelial dysfunction is an important change that occurs in aged vessels, within increased cardiovascular risk and declining with advancing age [37, 38]. Decreases in endothelium-dependent dilator responses are believed to be a result of alterations in NO bioavailability, activity of the endothelial NO synthase, and increased formation of reactive oxygen species (ROS) and gene expression. We here show that aged mice have attenuated endothelial function as determined in isolated aortic bioassay studies compared to young mice (max vasodilator response to acetylcholine (ACh) young % versus old %, , Figures 1(a) and 1(b)). Importantly, following a chronic dietary nitrate supplementation a reversal of endothelial dysfunction was observed (max. vasodilator response to ACh old nitrate % versus %, , Figures 1(a) and 1(b)). No effect was observed for endothelium-independent function in young, old and old nitrate treated mice (Figure 1(c)). Our results are in line with previous findings showing decreased age-related endothelial functions and improvements through nitrate supplementation in elderly humans [25]. We have previously shown that a chronic nitrate supplementation regimen affects plasma nitrite and nitrate levels [23]. Although this dosage is relatively high as compared to what can be ingested from a diet, it yields the equivalent effect a diet rich in vegetables exerts in humans [25, 26]. In mice, a diet rich in nitrite led to a decrease in inflammatory cytokines, suggesting an enhancement in NO bioavailability through inhibition of inflammation and inactivation of ROS [11].

(a)

(b)

(c)
3.2. Microarrays, Gene Ontology, and Pathway Analysis
Dietary interventions are known to influence cardiovascular functions and impact gene expression patterns [39–42]. To explore the genes that underlie and are responsible for the beneficial effects of nitrate supplementation in aged mice we performed gene expression microarray experiments on whole thoracic aortas of old and old nitrate treated mice. Microarrays allow the simultaneous quantification of gene expression differences, obtaining insights in health and disease [43]. Aging as well as dietary interventions lead to subtle changes of the transcriptional profile and accordingly, we choose a statistical significance in differential gene expression with a -value set to <0.05, not performing control for false discovery rate or family-wise error rate, as previously reported [10, 40, 44, 45].
Overall probe sets on the chip before filtering were 55682 and after exclusion of undetectable transcripts 26621 expressed transcripts were found. About 2900 transcripts were expressed differentially in aortic tissue of old and old nitrate treated mice, with totally 1565 upregulated and 1335 downregulated transcripts (Figure 2), which correspond to 2231 genes. Since the absolute numbers in gene expression differences do not provide explanations and comprehension of changes that occur through a dietary modification in aged mice, we performed GO and pathway analysis to gain insight into the transcriptional changes following nitrate supplementation.

Pathway analysis was performed through DAVID Bioinformatics database to identify the biological processes associated to differential regulation in old nitrate treated mice. DAVIDs functional annotations tool identifies enriched GO annotations categorized in biological processes, molecular functions and cellular components, with highest GO terms displayed in Table 1.
To determine the targets altered by a dietary nitrate supplementation in old mice we conducted DAVID’s pathway analysis [33, 46]. The highest enriched databases of pathways that were in our study were KEGG, Panther, and Biocarta. Figure 3 depicts overrepresented pathways in aortas of old nitrate treated mice with the corresponding gene count. One of the highest enriched pathways was “calcium signaling” in the KEGG database and the “calcium/calmodulin dependent protein kinase activation” in Biocarta (Figure 3).

The calcium signaling pathways were followed by pathways related to neurological and brain function “neuroactive ligand-receptor interaction,” “Alzheimer disease,” and “GABA-B receptor II signaling.” This might be of particular interest given a potentially impaired NO bioavailability in neurological diseases like Alzheimer’s disease [47, 48]. The possible link gained from our study is of importance since dietary nitrate was shown to improve regional brain perfusion in older adults and to improve reaction time in diabetic patients [49, 50]. NO signaling was additionally linked to neuroprotective signaling pathways and to promote survival and inhibit apoptotic cell death in neuronal cell types [51]. On the contrary, a recent study showed no effects on cognitive performance in older adults after short-term beetroot juice ingestion [52]. However, since we used a chronic supplementation regimen in mice a direct conclusion remains unfeasible.
3.3. Calcium Signaling Pathway
Aging leads to alterations of vascular gene expression which ultimately lead to increased vascular tone perpetuating arterial stiffness and hypertension. Increased age-related stiffness has been demonstrated to be an independent predictor of CVD events [53]. Comparing old and old nitrate treated mice we found KEGGs “calcium signaling” and Biocartas “calcium/calmodulin dependent protein kinase activation” to be the highest enriched pathways in old mice following nitrate supplementation. Calcium signaling is fundamental to vascular functions and the vascular tone is maintained by a highly regulated calcium homeostasis. In smooth muscle cells calcium is the main trigger for the action potential, calcium-dependent phosphorylation of the myosin light chain leads to an increase in muscle cell tone and the changes in intracellular calcium parallel the changes in contractile forces. We and others have recently shown that a dietary nitrate supplementation reduces arterial stiffness and blood pressure in older adults and healthy volunteers [25, 54, 55]. While acute effects are exerted by nitrate to nitrite bioconversion and subsequent reduction to NO, the long-term effects are incompletely understood. NO is known to modulate calcium signaling in the cardiovascular system, regulating important calcium handling proteins, for example, SERCA (sarcoplamatic reticulum calcium pump), ryanodine receptor, and L-Type calcium hannel, which are responsible for smooth muscle cell relaxation [56–58]. Recently, it was shown that ingestion of nitrate over 7 days increased sarcoplasmatic reticulum calcium release and tetanic contractile force production via alterations in cellular calcium handling with altered calcium handling protein expressions in nitrate treated mice [59]. We now identified differentially expressed genes related to calcium pathway in old versus old nitrate treated mice (Table 2 and Figure 4). Major players of the cellular calcium handling process were found, namely the L-type calcium channel (Cacna1d and Ppapdc2), the sarcoplasmatic calcium release channel (Ryanodine Receptor 2, Ryr2), the calcium/calmodulin-dependent protein kinase type II (Calm2, Camk2g, Camk4), and the inositol triphosphate receptor (Itpr2, Itpr3, Itpka; for detailed regulation and fold change expression difference refer to Table 2).

3.4. Detoxification, Antioxidation, and Free Radical Removal
The vascular aging process is characterized by increased oxidative stress and impaired antioxidant defense mechanisms [6, 60]. We have recently shown a distinctive transcriptomic profile with altered expression patterns in the old aorta with emphasis on the detoxification system [10]. Oxidative stress progresses due to excessive generation of ROS, to reduced antioxidant capacity or due to uncoupled nitric oxide synthase, amongst others. We have demonstrated that nitrite-derived NO impacts ROS generation and additionally reduces protein damage in myocardial ischemia/reperfusion injury [30, 61]. It was evidenced that NO regulates cardiovascular cell signaling and modulates cellular energetics [62, 63]. While in the acute setting nitrite is reduced by nitrite reductases and exerts vascular effects, the consequences of a long-term administration are unclear [64–66]. Recently, decreased inflammatory cytokines, inhibition of inflammation and inactivation of ROS following a nitrite rich diet in mice and healthy volunteers with subsequent improved age-related vascular functions were observed [11, 26, 67]. We now determined the enriched genes in aortas of old nitrate treated mice associated with the biological process of detoxification, antioxidation, and free radical removal as defined by Panther GO terms (Table 3). Most of the genes altered by nitrate in the old aorta were related to glutathione metabolism, the major antioxidant system. Remarkably, a downregulation of inducible nitric oxide synthase was noted, which could suggest a reduced inflammation in nitrate treated animals. Of note, the diversity of transcriptional regulation is observed through the bilateral regulation of superoxide dismutases (Table 3).
3.5. Limitations
We used the whole thoracic aorta instead of a single cell type. Although a specific cell might provide fundamental information, cell-cell interactions are imperative for coordinated organ functions and we thus intended to examine the complete transcriptional vascular changes that occur following nitrate supplementation. As we used only male mice our results might not apply to female mice and this should be taken into account for further studies. Moreover, alternative splicing variants are inadequately reflected by our analysis. Verification of microarray data is needed in forthcoming studies to confirm the beneficial effects of dietary nitrate on vascular tissue. Important targets to confirm would be the family of superoxide dismutases (SOD 1 and 3) the nitric oxide synthase 2, the platelet-derived growth factor beta (PDGFR-beta), and the adrenergic receptor. Also, measurement of markers of oxidative stress or inflammation would be helpful to verify the presented findings. Finally, it has to be noted that gene expression data does not always correspond to protein expression and activity.
4. Conclusions
Aging is the most important nonmodifiable cardiovascular risk factor, which leads to disadvantageous changes of the vasculature. A dietary intervention with inorganic nitrate has been evidenced to improve vascular functions. We determined changes that occur at a gene expression level and transcriptional profile following a chronic nitrate supplementation in aged mice. Our results highlight differentially expressed genes overrepresented in GO categories, pathways related to the calcium-signaling, the detoxification and antioxidation system. Our results contribute to existing knowledge of nitrate’s beneficial effects. This supports the concept of nutritional approaches to modulate age-related changes of vascular functions and its detrimental consequences.
Conflict of Interests
The authors declare no conflict of interests.
Acknowledgments
The authors would like to thank Sebastian Wohlfromm, Ph.D., BMFZ Genomics and Transcriptomics Laboratory (GTL), University of Düsseldorf, for conducting the microarray experiments. This study was supported by grants from the Forschungskommission of the Heinrich-Hein-University (Matthias Totzeck) and the German Research Foundation (DFG) to Tienush Rassaf (RA969/4-2) and to Malte Kelm (Ke405/5-1). Tienush Rassaf is a Heisenberg professor funded by DFG (Ra969/7-2).