Increased NBCn1 expression, Na+/HCO3- co-transport and intracellular pH in human vascular smooth muscle cells with a risk allele for hypertension.

Genome-wide association studies have revealed an association between variation at the SLC4A7 locus and blood pressure. SLC4A7 encodes the electroneutral Na+/HCO3- co-transporter NBCn1 which regulates intracellular pH (pHi). We conducted a functional study of variants at this locus in primary cultures of vascular smooth muscle and endothelial cells. In both cell types, we found genotype-dependent differences for rs13082711 in DNA-nuclear protein interactions, where the risk allele is associated with increased SLC4A7 expression level, NBCn1 availability and function as reflected in elevated steady-state pHi and accelerated recovery from intracellular acidosis. However, in the presence of Na+/H+ exchange activity, the SLC4A7 genotypic effect on net base uptake and steady-state pHi persisted only in vascular smooth muscle cells but not endothelial cells. We found no discernable effect of the missense polymorphism resulting in the amino acid substitution Glu326Lys. The finding of a genotypic influence on SLC4A7 expression and pHi regulation in vascular smooth muscle cells provides an insight into the molecular mechanism underlying the association of variation at the SLC4A7 locus with blood pressure.

Introduction

Hypertension is a common and a major etiological factor in cardiovascular, cerebrovascular and renovascular disease, estimated to cause up to 12.8% of global mortality and 3.7% of morbidity (1), with its complications contributing to a major worldwide economic burden (2). To compound the population health issue, up to 30% of patients with hypertension are inadequately treated (3), indicating a need for novel therapeutic approaches. There is a strong contribution from genes and their interaction with the environment to blood pressure (BP) regulation and the pathogenesis of hypertension (4,5). Recently, genome-wide association studies (GWAS) have identified a number of genomic loci at which common genetic variants have influences on BP (6–14). The study presented here attempts to translate the genetic information on a gene relevant to hypertension into a biomedical context. The identification of a path linking genetic variation to protein availability, protein activity and subsequently cellular behaviours provides further confidence in pursuing NBCn1 as a potential therapeutic target in hypertension.

One of the BP-associated loci identified by GWAS is on chromosome 3p24.1 encompassing the SLC4A7 (solute carrier family 4 member 7, HGNC: 11033) gene (8). The lead GWAS single nucleotide polymorphism (SNP), rs13082711 (NC_000003.12:g.27496418T > C), at this locus showed an association with diastolic BP above the genome-wide significance threshold and a similar direction for systolic BP (P = 3.8 × 10−9 and 1.5 × 10−6, respectively) in a large-scale GWAS by the International Consortium for Blood Pressure (8). This SNP is in strong linkage disequilibrium with 92 other SNPs (LD, r2 > 0.8; Supplementary Material, Fig. S1A), with only one being a non-synonymous polymorphism (rs3755652; NC_000003.12:g.27431445C > T; Glu326Lys).

The SLC4A7 gene encodes NBCn1, an electroneutral Na+/ co-transporter (15). At typical intracellular and extracellular ion concentrations, NBCn1 mediates electroneutral symport of Na+ and into cells. Amongst cells expressing NBCn1 are those that play important roles in BP control, including vascular smooth muscle cells (VSMCs) (16–18), endothelial cells (VECs) (16,17,19) and epithelial cells of the medullary thick ascending limb of the loop of Henle (17,20,21). NBCn1 is known to have a role in controlling pH in VSMCs and VECs (18,19); and pH is a determinant of VSMC contractility (19,22–24) and endothelial function (19,23–26), both of which impact BP. In support, it has been shown that pH dysregulation in SLC4A7 knockout mice results in an altered BP phenotype where the knockout mice show resistance to hypertensive stimuli such as angiotensin II, as compared to wildtypes (19).

In this study, we investigated if the BP-associated genetic variants at the SLC4A7 locus identified by GWAS affect SLC4A7 expression, NBCn1 protein availability and pH regulation in VSMCs and VECs. This is important for understanding the functional effects of these genetic variants and the mechanisms associated with their influence on BP.

Results

Allelic difference in SLC4A7 expression level

The lead BP GWAS SNP rs13082711 is in high LD (r2 ≥ 0.8) with 92 other SNPs which together span a 134 kb genomic interval. SLC4A7 is the only gene located within this interval (Supplementary Material, Fig. S1A) and the genes closest on either side are NEK10 (NIMA-related kinase 10, HGNC: 18592) and EOMES (eomesodermin, HGNC: 3372) residing 255 and 232 kilobases, respectively, away from the lead BP SNP. To investigate if the BP-associated variants influence the expression of SLC4A7, NEK10 and/or EOMES, we first performed RT-PCR assays of these genes to determine if they were expressed in VSMCs and/or VECs. These assays showed that SLC4A7 was consistently expressed in both cell types, whereas neither NEK10 nor EOMES was expressed in either cell type (Supplementary Material, Fig. S1B). Having found that SLC4A7 was expressed, but both NEK10 and EOMES were not readily detected in VSMCs and VECs, we conducted subsequent experiments focusing on SLC4A7, starting with allelic expression imbalance analyses. Since this technique entailed the analysis of a SNP in the coding region, we analyzed the non-synonymous coding SNP rs13096477 (NC_000003.12:g.27448703T > C) which is in strong LD (r2 > 0.95) with the BP GWAS index SNP rs13082711 located upstream of the gene and therefore not directly analyzed in this assay. The analyses showed that in both VSMCs and VECs, SLC4A7 RNA expression level of the BP-raising (minor, C) allele at rs13096477 was higher as compared to that of the alternative (major, T) allele (Fig. 1A and B).

Figure 1

The BP risk (minor) allele at the SLC4A7 locus is associated with increased gene expression. (A) Allelic imbalance analyses using heterozygotes at an exonic SNP (rs13096477), which is in high LD with the GWAS lead SNP (rs13082711). Representative chromatographs of genomic DNA (gDNA) and complementary DNA (cDNA) for the same sample, with the C/T SNP rs13096477 highlighted in blue, which is quantified in (B), where the C (BP-risk) allele is preferentially observed in cDNA of VSMC and VEC samples. *P<0.05, **P<0.01, by paired t-test.

To investigate the potential molecular mechanism leading to the allelic difference in SLC4A7 RNA level described above, we ascertained whether there was an allelic effect on DNA-nuclear protein interaction, an important aspect in gene transcriptional regulation. To this end, we undertook formaldehyde-associated isolation of regulatory elements (FAIRE) studies coupled with allelic imbalance assay to investigate if there was an allelic-dependent difference in nuclear protein binding in the intact cellular environment. Since this technique required analysis of an intragenic SNP, as for the allelic expression imbalance assay described above, we again analyzed the SNP rs13096477. The analysis showed an allelic imbalance of rs13096477, with its minor (BP-elevating, T) allele being in preferentially open chromatin conformation in VSMCs (Supplementary Material, Fig. S2A). This allelic imbalance in DNA-protein binding was not detected in VECs.

Since rs13096477 and rs13082711 are in high LD (r2 ≥ 0.8) with 92 other SNPs, and the FAIRE allelic imbalance assay technique could not provide information about which of these SNPs in high LD was responsible for the allelic imbalance results described above, we carried out electrophoretic mobility shift assays (EMSA) on 10 selected SNPs among the 93 in high LD, prioritized based on their positions relative to the start of SLC4A7 transcription and transcription factor binding site predictions from various bioinformatics resources including ENCODE and RegulomeDB (Supplementary Material, Fig. S2B). EMSAs on the 10 tested SNPs showed VSMC nuclear protein binding preferentially to the major allele of rs13096477 but, the minor allele of rs2371065 (NC_000003.12:g.27453098A > C), and weakly to the minor allele of rs13077400 (NC_000003.12:g.27431575A > G) (Fig. 2A, Supplementary Material, Fig. S3A–C). In contrast, neither allele of the BP GWAS index SNP rs13082711 showed DNA-protein interaction in EMSAs. Consistent with the FAIRE results mentioned above, EMSAs with VEC nuclear protein extracts did not detect binding with any of the DNA probes for the 10 tested SNPs (Fig. 2B).

Figure 2

The SLC4A7 BP-associated locus shows allele-dependent DNA-nuclear protein interactions. (A) EMSA with labelled-probes corresponding to both alleles of the 10 highest priority SNPs at SLC4A7. The start codon NEK10 appears on the right or of the figure. Unbound probes as visualized at bottom of autoradiograph. Representative figure of n=3 replicates. (B) EMSA comparing VEC and VSMC nuclear extracts for three SNPs shows the shift present in VSMC nuclear extracts virtually absent with VEC nuclear extracts. Representative figure of n=2 replicates. All lanes had 10 femtomoles of labelled double-stranded oligonucleotide probe with 10 ng of corresponding nuclear extracts.

Since the above EMSAs showed allele-specific binding of a nuclear protein to the rs2371065 minor allele, we undertook further experiments in an attempt to determine the identity of this nuclear protein. We performed supershift assays with antibodies against the nuclear proteins PHOX2A, PLAG1 and TFAP2C, respectively, as bioinformatics analyses showed that their recognition DNA sequences had some similarities with the DNA sequence encompassing the rs2371065 site. However, these assays show no evidence to suggest that any of these proteins being the one interacting with the rs2371065 minor allele shown in EMSAs (Supplementary Material, Fig. S4). We next conducted a DNA pull-down assay using a biotin-labelled double-stranded oligonucleotide corresponding to the DNA sequence at and surrounding the site of the rs2371065 minor allele. An electrophoretic analysis of proteins pulled down by this oligonucleotide showed two bands consistent with the EMSA results described earlier, neither of which was present in the negative control (Supplementary Material, Fig. S5). However, N-terminal protein sequencing failed to reveal the identity of the protein(s) pulled down by the oligonucleotide.

Allelic difference in NBCn1 protein level

Having found an allelic difference in SLC4A7 RNA expression level as described above, we investigated if there was a corresponding difference in NBCn1 protein level. Immunoblot analyses of VSMC total cellular protein extracts showed that NBCn1 protein levels were the greatest in minor (BP-raising) allele homozygotes, intermediate in heterozygotes and lowest in major allele homozygotes (Fig. 3A and B). Despite the allelic imbalance in RNA/cDNA levels, there were no detected genotypic differences in VEC total cellular NBCn1 protein expression (Fig. 4A and B).

Figure 3

The BP risk (minor) allele at the SLC4A7 locus is associated with increased NBCn1 protein expression and availability at the plasma membrane in VSMCs. (A) Protein immunoblots for NBCn1 with VSMC samples stratified by genotype at rs13082711, where the two signals for NBCn1 reflecting the isoforms with and without Cassette II (band sizes approximately 150 kDa and 135 kDa, respectively), quantified in (B), where comparisons across the three genotypes show increased protein expression for risk allele carriers. *P<0.05, by one-way ANOVA after Bonferroni correction for multiple comparisons. (C) Subcellular fractionation samples from VSMCs, showing that the BP risk (minor) allele at the SLC4A7 locus is associated with increased NBCn1 signal in the membrane fraction. 2 μg protein of each cellular fraction was loaded in each lane. n=9 in each group *P<0.05 by Mann-Whitney U-test.

Figure 4

The BP risk (minor) allele at the SLC4A7 locus is associated with increased NBCn1 availability at the plasma membrane in VECs. (A) shows the protein immunoblots for NBCn1 with VEC samples stratified by genotype at rs13082711, and quantified in (B) P > 0.05 by Mann-Whitney U-test. (C) Subcellular fractionation samples from VECs, showing that the BP risk (minor) allele at the SLC4A7 locus is associated with increased NBCn1 signal in the membrane fraction. 2 μg protein of each cellular fraction was loaded in each lane. n=8 each group. *P<0.05 by Mann-Whitney U-test.

Following the above analysis of total cellular NBCn1, we examined NBCn1 in subcellular locations, more specifically, their presence at the cellular membrane. An immunoblot analysis of subcellular protein fractionations prepared by differential centrifugation showed that the cells carrying the BP-raising allele (C/T genotype for rs13082711) had higher levels of NBCn1 in the membrane fraction as compared to non-carriers (T/T genotype) in both VSMCs (Fig. 3C) and VECs (Fig. 4C). The vesicle and organelle fraction also contained NBCn1 but the difference between the genotypes was less pronounced. The supernatant (cytosolic) fraction had no detectable expression of NBCn1 or the plasma membrane markers N-cadherin and VE-cadherin.

Association of the BP-raising allele with increased Na+/ co-transport activity in VSMCs, which is not overcome by Na+/H+ exchange activity

To determine whether the increased SLC4A7 gene expression and NBCn1 protein availability at the plasma membrane in cells carrying the BP-raising allele gives rise to increased protein function, the capacity of these primary cultured cells for pH recovery from -prepulse-induced intracellular acidification was assessed. Comparing VSMCs from BP-raising allele carriers (C/T risk allele carriers at rs13082711, n = 5) and protective allele homozygotes (T/T genotype, n = 7), we observed a difference in pH recovery rate following intracellular acidosis in cells exposed to the Na+/H+ exchange inhibitor dimethylamiloride (DMA, 30 µM) in the presence of CO2/ (Fig. 5A). As the activities of acid-base transporters are regulated by pH, the rate of net base uptake was determined at regular pH intervals (Fig. 5B), demonstrating a genotype effect, where VSMCs from the rs13082711 BP-raising allele carriers had a higher rate of net base uptake at each pH compared to protective allele homozygotes. In addition to differences in net base uptake rates, there was also a genotype effect on final plateau pH that was revealed in the presence of DMA (ΔpH=0.087 ± 0.027, P < 0.01) (Fig. 5C). This observed genetic influence was confirmed to be CO2/-dependent as the difference in net base uptake and final resting pHdisappears in CO2/-free conditions (Supplementary Material, Fig. S6A–C).

Figure 5

The BP risk (minor) allele at the SLC4A7 locus is associated with increased Na+/ co-transport activity in VSMCs and is not overcome by Na+/H+ exchange activity. Intracellular pH recovery of VSMCs following -prepulse-induced intracellular acidosis. Experiments were performed with (A–C) CO2/ and 30 µM DMA, (D–F) CO2/ and 30 µM DMA, in cell cultures serum-starved for 48 h beforehand, and (G–I) in CO2/ conditions. The risk allele carriers (orange, C/T at rs13082711, n=5) had faster pH recovery and higher final plateau pH compared to protective allele homozygotes (green, T/T at rs13082711, n=7), even in the absence of DMA. (A,D,G) Intracellular pH traces of VSMCs during ammonium prepulse studies. Grey error bars denotes SEM. Periods of exposure to 20 mM NH4Cl, Na+-free buffers, or 30 µM DMA are marked above the traces. X-axis scale bar denotes 5 min. Groups compared by repeated-measures two-way ANOVA over the 280 time points between 60 s and 900 s after reintroduction of Na+-containing buffer. (B,E,H) Net base uptake for each group calculated at specified pH values. Slopes compared by least-squares linear regression analyses. (C,F,I) Final plateau pH after recovery from intracellular acidosis. *P<0.05, **P<0.01 by Mann–Whitney U-test.

To minimize the potential effects of variances in cell cycles, the VSMCs were also tested after 48 h of serum-starvation (27). In this synchronized state, the Na+- and CO2/-dependent, DMA-insensitive pH recovery and the final plateau pH were lower than under serum-stimulated conditions but still higher in BP-raising allele carriers compared to cells homozygous for the protective allele (Fig. 5D–F).

While these data show a genetic influence on Na+/ co-transport activity in the absence of Na+/H+ exchange activity (inhibited by DMA), it does not necessarily indicate an overall effect on pH under physiological conditions where Na+/H+ exchange may play a large role. It is therefore important to note that in the presence of Na+/H+ exchange activity (i.e. in the absence of DMA), the genotype-associated difference in net base uptake and final plateau pH persisted (Fig. 5G–I).

As expected, the buffering capacity of the VSMCs was higher in the presence of CO2/ than in its nominal absence. This was particularly evident at the pH ranges closer to physiological levels where [] is high and adds substantially to the buffering power (Supplementary Material, Fig. S7A). There were no differences in buffering capacity between cells of the two assessed genotypes (Supplementary Material, Fig. S7B) in the presence or absence of CO2/.

Association of the BP-raising allele with increased NBCn1 activity in VECs, which is overcome by Na+/H+ exchange activity

In a study of VEC samples from BP-raising allele carriers (C/T genotype at rs13082711, n = 10) and protective allele homozygotes (T/T genotype, n = 10), no significant difference in pH recovery following intracellular acidosis was detected by repeated measures two-way ANOVA (Fig. 6A) despite the apparent divergence of the curves in cells exposed to Na+/H+ exchange inhibition (30 µM DMA) in the presence of CO2/. When calculated and plotted as a function of pH, this revealed a genotype effect on the rate of net base uptake (Fig. 6B). Similar to VSMCs, there was a higher final plateau pHi that was revealed in the presence of DMA for C/T risk allele carriers compared to T/T protective allele homozygotes (ΔpH=0.082 ± 0.053, P = 0.10, Fig. 6C). This is a similar magnitude to that of VSMCs, but due to the larger confidence interval, did not reach statistical significance. Similarly, there was no significant difference in buffering capacities between the genotypes of VECs when compared in the pH range covered by both the C/T and T/T groups (Supplementary Material, Fig. S7C and D). Once again, the observed genetic influence on net base uptake was confirmed to be CO2/-dependent with its disappearance in CO2/-free conditions (Supplementary Material, Fig. S6D–F).

Figure 6

The BP risk (minor) allele at the SLC4A7 locus is associated with increased Na+/ co-transport activity in endothelial cells but it is masked by Na+/H+ exchange activity. Intracellular pH recovery of VECs following -prepulse-induced intracellular acidosis. Experiments were performed with (A–C) CO2/ and 30 µM DMA and, (D–F) in CO2/ conditions. The risk allele carriers (orange, C/T at rs13082711, n=10) had faster pH recovery and higher final plateau pH compared to protective allele homozygotes (green, T/T at rs13082711, n=10), but this difference was abolished in the absence of DMA. (A and D) Intracellular pH traces of VECs during ammonium prepulse studies. Grey error bars denotes SEM. Periods of exposure to 20 mM NH4Cl, Na+-free buffers or 30 µM DMA are marked above the traces. X-axis scale bar denotes 5 min. Groups compared by repeated-measures two-way ANOVA over the 280 time points between 60 s and 900 s after reintroduction of Na+-containing buffer. (B and E) Net base uptake for each group calculated at specified pH values. Slopes compared by least-squares linear regression analyses. (C and F) Final plateau pH after recovery from intracellular acidosis.

Similar to the analysis of VSMCs above, a more physiological condition with uninhibited Na+/H+ exchange needed to be considered. In this setting, unlike that of VSMCs, the genotype-associated difference in VEC net base uptake rate and final plateau pH was overcome by the presence of Na+/H+ exchange activity (Fig. 6D–F). Taken in combination, this may indicate that the BP-associated locus is more likely to exert its effects via VSMCs rather than VECs.

Bioinformatic tools predict the NBCn1 Glu326Lys variation to be well tolerated

Having identified that the BP-raising allele (C for rs13082711) is associated with higher Na+/ co-transport activity and steady-state resting pHi, we wondered whether the missense Glu326Lys amino acid change (SNP at rs3755652) in high LD with rs13082711 (r2 > 0.9) contributes to this effect. The main difference between glutamic acid and lysine amino acid residues is that the former has a negative charge; whist the latter has a positive charge. They are otherwise relatively similar - having similar molecular weights, both being hydrophilic, and neither is bulky enough to result in steric hindrance or changes in secondary structures. This similarity is reflected by the early analysis by Grantham (28), noting that the functional effect of this amino acid change is predicted to be relatively small. Utilizing seven other online prediction tools, the overall consensus was that the Glu326Lys is a well-tolerated amino acid change (Supplementary Material, Table S1).

The NBCn1 Glu326Lys variant does not affect NBCn1 activity

To assess the hypothesis that the Glu326Lys variant does not affect NBCn1 activity, we overexpressed the 326Glu and 326Lys variants and a shorter variant that lacks splice Cassette II (amino acids 251–374, thus lacking the Glu326Lys variation) of NBCn1 in A10 cells, and assessed their capacity for pH recovery following -prepulse-induced intracellular acidosis in the presence of CO2/ and DMA. Cells transfected with any of the three overexpression plasmids displayed faster pH recovery from intracellular acidification compared to cells transfected with the control vector, but there were no differences between the impact of the three overexpression plasmids (Fig. 7A and B). This was verified in the analysis accounting for net base uptake rate as a function of pH (Fig. 7C). In addition to the increased pH recovery rate, A10 cells transfected with overexpression plasmids also have a higher plateau pH as compared to cells transfected with the control vector, but again, there were no differences between cells transfected with the three overexpression plasmids (Fig. 7D). These findings were obtained with consistent and comparable overall overexpression across the three different NBCn1 variants (Fig. 7E and F), and similar to primary VSMCs and VECs observed before, there was no difference in buffering capacities (Supplementary Material, Fig. S7E). These findings suggest that the Glu326Lys amino acid variation does not alter the intrinsic acid-base transport activity of NBCn1, consistent with the aforementioned bioinformatics predictions.

Figure 7

The Glu326Lys amino acid variation at the BP risk locus does not influence NBCn1 activity. Intracellular pH recovery of A10 cells transfected with SLC4A7-expression plasmids (pcDNA) following -prepulse-induced intracellular acidosis. Experiments were performed with CO2/ and 30 µM DMA. Blue denotes Vector; Brown, 326K (lysine); Black, 326E (glutamic acid); Yellow, without Cassette II. Transfection with any of the three expression plasmids increased pH recovery and final plateau pH compared to vector control (n=6 sets). (A) Intracellular pH traces of A10 cells transfected with SLC4A7-expression plasmids during ammonium prepulse studies. Periods of exposure to 20 mM NH4Cl, Na+-free buffers or 30 µM DMA are marked above the traces. X-axis scale bar denotes 5 min. Groups compared by repeated-measures two-way ANOVA over the 280 time points between 60 and 900 s after reintroduction of Na+-containing buffer. (B) shows the zoomed in section between time points 0 and 600 s after reintroduction of Na+-containing buffer. Error bars denotes SEM. Groups compared by repeated-measures two-way ANOVA. (C) Net base uptake for each group calculated at specified pH values. Slopes compared by least-squares linear regression analyses. (D) Final plateau pH after recovery from intracellular acidosis. P<0.01, by repeated-measures one-way ANOVA; **P<0.01, by paired t-test after Bonferroni correction. (E) Representative immunoblot of samples after transfection with vector or overexpression plasmids, quantified in (F), where the cells transfected with overexpression plasmids had higher expression levels compared to vector control, *P<0.05, by paired t-test after Bonferroni correction. There were no differences when the three overexpression studies were compared to each other by repeated-measures one-way ANOVA.

Calcineurin inhibition did not influence Na+/ co-transport activity under the experimental conditions

Although there is evidence that NBCn1 activity is regulated by phosphorylation, perhaps via calcineurin (29,30), the absence of an effect of either the Glu326Lys amino acid variation, and the variant that lacks Cassette II (Fig. 7), led us to postulate that the genotypic effect was independent of calcineurin activity. Furthermore, the genotypic effect was observed in conditions with lower intracellular [Ca2+], where calcineurin is less likely to play a role (30). To assess this hypothesis, primary VSMCs were assessed for Na+- and CO2/-dependent, DMA-insensitive pH recovery rate following -prepulse-induced intracellular acidosis while being exposed to 10 µM of the calcineurin inhibitor FK506 (or equivolume DMSO vehicle). In these conditions, we did not observe an effect of 10 µM FK506 on DMA-insensitive Na+/ co-transport activity (Supplementary Fig. S8A and B) and it stands to follow that even in the presence of 10 µM FK506, the genotype-associated difference in DMA-insensitive Na+/ co-transport activity persists (Supplementary Material, Fig. S8C and D).

To confirm the above findings from primary VSMCs, A10 cells transfected with the three different overexpression plasmids (326Lys, 326Glu and the variant without splice cassette II) were also investigated under the same conditions. Consistent with the findings from primary VSMCs, A10 cells transfected with the different overexpression plasmids continues to show the lack of effect of calcineurin inhibition by 10 µM FK506 on Na+/ co-transport activity (Supplementary Material, Fig. S8E and F). Again these findings are consistent with previous reports that calcineurin inhibition only affects NBCn1 activity under conditions of elevated intracellular [Ca2+] (30).

Discussion

We here identify a potential mechanism for the impact of the SLC4A7 BP-associated locus on vascular cells. Mechanistically, our studies demonstrate allele-associated differences in DNA-nuclear protein interaction, gene expression and NBCn1 function, where the carriers of the SLC4A7 BP risk allele demonstrated increased NBCn1 protein expression, and in turn, upregulated net base uptake rate and higher steady-state pH. Although the increased expression and NBCn1 activity are also found in VECs, the effect on net acid extrusion rate and steady state pH was masked by a larger Na+/H+ exchange (DMA-sensitive) activity under our experimental conditions.

GWAS have identified multiple genetic loci associated with BP, each with a modest contribution to overall BP levels. One of these was found to encompass SLC4A7/NBCn1. NBCn1 regulates pH by electroneutral symport of Na+ and into cells. Reduction of NBCn1 activity in VSMCs and VECs, by knockdown (18) or knockout (19), abolishes Na+/ co-transport and markedly attenuates pH recovery from intracellular acidosis. As there appears to be opposing effects of altered NBCn1-function in VSMCs and VECs on BP demonstrated by the SLC4A7 knockout mouse (19), it is important to identify in which cell type, if any, NBCn1 allele variation has an effect.

The role of vascular pH regulation in human hypertension is supported by the finding that vascular segments from hypertensive patients are more resistant to noradrenaline-induced intracellular acidification than segments from normotensive controls (31). That study was conducted prior to in-depth characterization of Na+/ co-transporters, and more recent reports show that NBCn1-mediated Na+/ co-transport has subsequently been shown to protect VSMCs against intracellular acidification during contractions (18,30,32). As intracellular acidification of VSMCs lowers rho-kinase-dependent VSMC Ca2+ sensitivity (19,23), the improved ability to eliminate the contraction-induced intracellular acid load may contribute to the higher peripheral arterial resistance of hypertensive patients through increased rho-kinase activity. Consistent with the effect of sustained intracellular acidification on VSMC Ca2+ sensitivity, knockout of NBCn1 lowers noradrenaline-induced contractions of mesenteric arteries (19) and myogenic responses of pressurized middle cerebral arteries (24) after endothelial blockade without affecting VSMC membrane potential or [Ca2+]. There is also evidence that changes in acid-base transport function and/or pH impact in vitro VSMC proliferation (33), migration and viability (34), and medial wall thickness (23). NBCn1 plays a key role for VSMC migration and carotid artery remodeling most likely because it establishes local pH gradients and promotes filopodia, which can explain the decelerated directional migration of VSMCs from NBCn1 knockout mice (35). Altered NBCn1 activity could therefore modify vascular remodeling with long-term impact on peripheral resistance.

Endothelial function is impaired by intracellular acidosis (23,25); and endothelial NO production is reduced in SLC4A7 knockout mice without any change in endothelial [Ca2+] or NO synthase expression (19). Intracellular pH has also been shown to influence the generation of endothelial vasoactive substances such as the expression of endothelin (26) and the intrinsic activity of nitric oxide synthase (25). However, the disappearance of the genotypic effect in VECs when Na+/H+ exchange is present suggests that the GWAS identified genetic variance is unlikely to exert its effect through endothelial function. This is with the caveat of these studies being conducted in an in vitro system, whereas in vivo, endothelial cells would be exposed to shear stress as well as circulating hormonal factors.

We found that VECs have a higher level of SLC4A7 expression and faster Na+- and CO2/-dependent, DMA-insensitive recovery from intracellular acidosis as compared to VSMCs (Supplementary Material, Fig. S9). Importantly, the genotype-associated differences found in Na+/ co-transport activity are masked by Na+/H+ exchange in VECs, but not VSMCs. This may be related to the relative contributions of Na+/H+ exchange and Na+/ co-transport being similar to each other in VSMCs, especially in the pH ranges from around 6.7 upwards (Fig. 8A). This would enable subtle genotype-associated differences in NBCn1 activity to persist in VSMCs, but not in VECs where Na+/H+ exchange contribution far outweighs that of Na+/ co-transport (Fig. 8B). These results should be taken in relation to the expression of other pH regulators (Fig. 8C), where a panel of paired VSMC and VEC samples shows inter-sample variability with SLC4A7 short and long isoforms (confirming the previous qRT-PCR and immunoblot results), but also SLC4A4 (NBCe1) and SLC9A1 (NHE1). Although mRNA expression of multiple SLC4-family Na+/ co-transporters has also been identified in mouse carotid arteries, NBCn1 has been found to functionally dominate net acid extrusion (35). Notably, we found a minimal signal for SLC4A4 in VECs.

Figure 8

Na+/ co-transport and Na+/H+ exchange activities are similar in VSMCs, but Na+/H+ exchange dominates in VECs. Average net base uptake plotted against average pH for each corresponding time-point after -prepulse-induced intracellular acidosis of (A) VSMCs (n=12) and (B) VECs (n=20) in the presence and absence of either 30 µM DMA or CO2/. Open circles denote studies with CO2/; cross, CO2/ + DMA; plus, CO2/-free conditions. (C) shows result of end-point RT-PCR (10 ng of reverse-transcribed RNA) of selected pH regulators (long and short variants of SLC4A7, together with other pHi regulators SLC4A4 (NBCe1) and SLC9A1 (NHE1)) in paired samples of VSMCs and VECs, showing significant inter-sample variability. β-actin included as positive control.

The genotype-associated changes in overall VSMC pH regulation may be related not just to resting pH, but also to the rate from which pH recovers from intracellular acidosis; particularly as pH is not constant and pH recovery rates are important in response to a variety of in vivo vasoconstrictive stimuli such as angiotensin II, noradrenaline, endothelin-1 and cellular depolarization (31,36–38). It is hitherto unknown whether it is the resting steady-state pH, or the ability to return towards its set-point and avoid intracellular acidification during contractions that predominantly influences the overall in vivo phenotype but our results show that both the rate of recovery and the resting levels of pH are influenced by SLC4A7 genotype.

The direction of effect for the data presented here is congruent with that of the SLC4A7 knockout model (19); where the knockout mice are protected from the hypertensive stressor angiotensin II, we show that the BP-protective allele is associated with the reduced NBCn1 expression and slower recovery from intracellular acidosis, particularly in VSMCs. The clinical relevance of the described genotypic effects is supported by observational data from arterial segments from humans with hypertension: in human resistance arteries from hypertensive subjects, the VSMCs were protected from intracellular acidosis after noradrenaline contraction relative to those from normotensive subjects (31). Additionally, when intracellular acidosis was induced in mouse mesenteric artery studies, a fall of VSMC pH by approximately 0.1 was associated with a decreased ex vivo contractility (19). The differences in resting pH (95% CI, 0.08–0.12) and pH recovery rates (95% CI, 15–55%) that we describe between the SLC4A7 genotypes – with both parameters higher in rs13082711 risk allele carriers as compared to the protective allele homozygotes – may be clinically significant if they continue to persist in vivo.

A recently reported study shows that SNP rs820430, which is located in an intergenic region ∼23 kb from SLC4A7 and associated with BP in a Chinese Han population GWAS (39), has an effect on SLC4A7 transcription in peripheral blood mononuclear cells (40). In the present study, we investigated a separate genetic signal represented by rs13082711, the lead BP-associated SNP in GWASs in Europeans (8), which is not in high LD with rs820430 (r2 = 0.178 in CEU). Our study shows that rs13082711 also affects SLC4A7 expression and reveals that it leads to altered NBCn1 protein levels in VSMCs and VECs which are cell types that play important roles in controlling BP. Importantly, our study reveals, for the first time, that the BP-associated genetic variant impacts on pH regulation, which has a direct and important implication in BP control.

It should be acknowledged that apart from VECs and VSMCs, there are other tissues such as the medullary thick ascending limb of the loop of Henle where SLC4A7 may exert an effect on BP regulation. Additionally, data from the Genotype-Tissue Expression (GTEx) Project indicates that SLC4A7 is also expressed in other tissues such as mammary tissue, transformed fibroblasts, prostate, lymphocytes, nerves and subcutaneous adipose tissue. Furthermore, the same database shows an association between the SNP rs13082711 and the expression level of the neighbouring gene, NEK10, in left ventricles. Possible effects on these other tissues were not explored in this study, but remain a potential parallel mechanism of action for the BP-associated SNP.

In summary, the work presented here has revealed that the BP-raising allele of rs13082711 is associated with allele-dependent DNA-nuclear protein interactions, greater SLC4A7 transcript levels, higher NBCn1 protein levels and plasma membrane availability particularly in VSMCs, and increased rates of base uptake via Na+/ co-transport and higher resting pH, once again more apparent in VSMCs. Notably, due to the large contribution of Na+/H+ exchange to VEC pH regulation, the subtle genotypic effect is masked when Na+/H+ exchange is present. This is not the case for VSMCs, where the genotypic effect persists even in the presence of active Na+/H+ exchange. The missense variant Glu326Lys which is in high LD with rs13082711, had no significant effect on NBCn1 function, suggesting that the functional effect of the BP-associated variant is predominantly related to an influence on SLC4A7 expression levels. The identification of a pathophysiological path from BP-associated genetic variation, to gene expression, and subsequently to gene function that alters cellular behaviour supports these genetic studies as a route to the discovery of drug targets.

Materials and Methods

Study samples

This study had ethical approval from Queen Mary, University of London (Protocol No.: Plaque-WHRI-01; NRES ref: 08/H0704/140, and subsequent amendments). Umbilical cords for cell isolation were obtained from the Royal London Hospital which serves the east London population where the two largest self-reported ethnic groups were Bangladeshi (32%) and White British (31%), with 21% of households being multi-ethnic. All tissue samples were fully anonymous before distribution to the recipient analysis groups, as per ethical approval. We derived primary cultures of human umbilical artery vascular smooth muscle cells (VSMCs) based on cell explants cultured on 0.2% w/v gelatin-coated polystyrene as established previously (41). We derived primary cultures of human umbilical vein endothelial cells (VECs) based on endoluminal collagenase digestion, cultured on 0.2% w/v gelatin-coated polystyrene as established previously (42).

The A10 rat thoracic aorta vascular smooth muscle cell line was obtained commercially from ATCC (CRL-1476) for overexpression studies. All cells were used before a maximal passage number of 6. Where indicated, serum-free media was used to induce synchronised cell arrest of VSMCs (27).

Genotyping

DNA was extracted from cell culture samples (Qiagen, 69509) and genotyped in an array using the KASP™ method (Applied Bioscience). The primers used are listed in Supplementary Material, Table S2. The minor allele frequency of the study population (19.8%) approximates those in EUR (21.5%) and SAS (18.1%) in the 1000 Genomes Project samples.

Allelic imbalance analyses

Total cellular RNA was isolated from cell culture samples (Macherey-Nagel, 740955), and reverse transcribed (Promega, M170). End-point RT-PCR was conducted using the primers listed in Supplementary Material, Table S3. The PCR products were size-separated via agarose gel electrophoresis, excised and cleaned up (Promega, A9281) before Sanger sequencing was conducted by the commercial service at the Genome Centre, Queen Mary University of London. The PeakPicker software (McGill University and Génome Québec Innovation Centre, California) was used to determine relative allele ratios from heterozygotes of interest (43).

Formaldehyde-assisted isolation of regulatory elements

Formaldehyde-assisted isolation of regulatory elements (FAIRE) was conducted as adapted from the commercial kit from Millipore (#16-201). Briefly, formaldehyde was added to allow protein-DNA cross-linking, with glycine further added to quench the residual formaldehyde. The cells were lysed and the DNA sheared by a probe sonicator (Jencons Scientific, 690-024) for three pulses of 10 s, producing sheared DNA of lengths between 200 and 800 base pairs. The samples were de-crosslinked by heating, and the residual proteins were digested by proteinase K. The DNA fragments were then suitable for isolation by DNA extraction columns (Promega, A9281). The isolated chromatin-associated DNA was used for allelic imbalance analyses as described above.

Electrophoretic mobility assay (EMSA) and DNA pulldown

Commercial complementary single-stranded oligonucleotides (either unlabelled or biotin-labelled) were annealed to form double-stranded oligonucleotides (see Supplementary Material, Table S4 for sequences). Nuclear extracts were isolated from VSMCs and VECs (Millipore, #2900). Both the nuclear extracts and double-stranded oligonucleotides were used for both EMSAs and DNA pulldown assays.

For EMSAs, the reaction mixture was combined with study-specific quantities of nuclear proteins, labelled antibody and unlabelled competitor antibodies. Unless where specified, quantities of nuclear proteins were 10 µg (2 µg/µl) and labelled oligonucleotides were 10 femtomoles (10−14 moles). When required for supershift assays, 1 or 5 ng of the antibody (rabbit anti-AP2γ IgG antibody, Santa Cruz, sc-8977 X; goat anti-PHOX2A IgG antibody, Santa Cruz, sc-13229 X; or rabbit anti-ZAC1 IgG antibody, Santa Cruz, sc-22811 X) was pre-incubated with the nuclear proteins before adding to the reactions. To account for the variable volumes of reactants, the final reaction volume was made up to 20 µl with nuclease-free water. The reaction mixture was electrophoresed on a 4% non-denaturing acrylamide gel, electrotransferred onto a nylon membrane (Amersham, RPN119B) and crosslinked with oven drying followed by 254 nm ultraviolet light at a dose of 0.12 J/cm2 (Biolink, BLX-254E). The biotin-labelled double-stranded oligonucleotides were detected by the LightShift® Chemiluminescent EMSA Kit (ThermoScientific, #89880).

DNA pulldown assays were conducted using streptavidin-conjugated agarose beads mixture (Sigma, S1638) based on manufacturer’s protocol as summarised below. For every reaction, 1 nmol of biotin-labelled double-stranded oligonucleotides were immobilised onto 50 µl of streptavidin-conjugated agarose beads mixture and washed. The bead-oligonucleotide complex was incubated with 500 µg of nuclear extract, and subsequently, the bead-oligonucleotide-protein complex was washed. The pulled-down proteins were eluted with 100 µl of 2x Laemmli buffer at 70°C. The eluate underwent SDS-PAGE gel electrophoresis and Ponceau S staining.

Immunoblotting and subcellular fractionation

For immunoblotting, total cellular protein was isolated using RIPA buffer (50 mM Tris pH 7.4, 150 mM NaCl, 0.5% w/v sodium deoxycholate, 1% v/v NP-40, 0.1% w/v sodium dodecylsulfate) supplemented with protease inhibitors.

For subcellular fractionation by differential centrifugation, cells were homogenized with a homogenization buffer (10 mM Tris HCl pH 7.2, 1 mM EDTA pH 8.0, 250 mM sucrose) supplemented with protease inhibitor mixture. The cell homogenates were subjected to differential centrifugation with the supernatant removed at each step and resuspension in RIPA supplemented with protease inhibitors. The centrifugation sequences were: 900 g for 10 min (nuclear and cell debris fraction), 10,000 g for 5 min (mitochondrial/lysosomes/peroxisomes fraction) and 100,000 g for 60 min (membrane fraction). Each separate pellet was resuspended in RIPA buffer. The final supernatant was kept as the cytoplasmic fraction.

The samples derived from either method were electrophoresed using SDS-PAGE gels, electrotransferred to a PVDF membrane (Amersham, 10600021) and detected by enhanced chemiluminescence (Amersham, RPN2232). The antibodies used are listed in Supplementary Material, Table S5. Protease inhibitors reached the final concentration of 1 mM phenylmethanesulfonyl fluoride (Sigma, P7626), 2 µM leupeptin hemisulfate (Sigma, L2884), 1.5 µM pepstatin A (Sigma, P2465) and 0.15 µM aprotinin (Sigma, A1153) after addition to the samples.

Intracellular pH (pH) and Na+/-dependent pH recovery from intracellular acidosis

Experiments were conducted with cells cultured on 0.1% w/v poly-lysine (Sigma, P8920)-coated polystyrene flasks (BD-Falcon, #353107). Studies were conducted in a temperature-controlled stand and ambient air warmed to 37°C for 30 min before initiation of experiments. Ports were created to allow for aspiration or bubbling of gasses (5% CO2/95% air or 100% air) at a constant rate. Cells were incubated with 1 µM BCECF-AM (ThermoFisher, B-1170) for 20 min before being rinsed. A baseline F/F ratio was observed for 5 min before paired digital images at excitation wavelengths of 495 nm and 440 nm were obtained with exposure times of 800 and 1000 milliseconds respectively (Visitron Systems, Visiview Version 3), obtained every 3 s at periods of interest, otherwise at 15-s intervals. The initial resting baseline and final plateau was recorded for 1 min. A calibration curve was obtained at the end of the study where the cells are exposed to high-potassium buffers of varying pH with 5 mg/L nigericin (Sigma, N7143). These curves were highly reproducible, with coefficients of variation of <5%, and were near-linear in the pH section of interest (Supplementary Material, Fig. S10). For the constituents of buffers, see Supplementary Material, Table S6. The buffering capacity (mmol/L) was calculated based on the response to wash-out of NH4Cl using the formula β = Δ[]/ΔpH. The [] was calculated from the Henderson-Hasselbalch equation assuming equilibration of NH3 across the membrane and using a pKa of 9.2, and ΔpH calculated by the difference in pH before and after washout of NH4Cl (18) The rate of net base uptake (mmol/L/min) during recovery from intracellular acidosis was calculated as the product of the rate of pH gain (ΔpH per minute calculated by linear regression over 15 s) and the buffering capacity (18).

Cloning

The coding sequence of SLC4A7 (Homo sapiens solute carrier family 4 member 7, transcript variant 1, mRNA Sequence ID: ref|NM_003615.4) was subcloned to pcDNA3.1(+) vector. Site-directed in vitro mutagenesis was conducted using the QuikChange II Site-Directed Mutagenesis Kit (Agilent, #200523), using the primer pair, forward: CCTGACGCTGACTC TCT TGG GAACTTCTGGAGG; reverse: CCTCCAGAAGTTCCCAAGA GAGT CAGCGTCAGG. PCR-driven overlap extension to generate the splice variant lacking Cassette II (Homo sapiens solute carrier family 4 member 7, transcript variant 3, mRNA Sequence ID: ref|NM_001258380.1) was based on the protocol described by Heckman and Pease (44), using the primer pair, forward: CTTGAAAGGAATGGTATTTTGGCCTCT; reverse: AGAGGCC AAA ATACCATTCCTTTCAAG. All plasmids had their sequences confirmed by Sanger sequencing.

Cell culture and transfection

Transfections for A10 cells were conducted using liposome-based transfection with the X-tremeGENE reagent (Roche, 06366244001). For 10 cm2 surface area, ratios of 1 µg of plasmid, 200 µl of DMEM (Sigma, M4530) and 3 µl of the X-tremeGENE reagent were allowed to incubate for 20 min at room temperature before being added drop-wise into the culture surface already containing 2 ml of fresh media. The cells were then incubated in a humidified incubator kept at 37°C and 5% CO2 for 48 h without the need to replace the culture media prior to experiments.

Statistical analyses

Comparisons between two independent groups were conducted using unpaired, two-tailed Student’s t-test or the Mann-Whitney U-test for parametric and non-parametric distributions, respectively. Comparisons between two paired groups were conducted using two-tailed, one-sample t-test or the Wilcoxon sign-ranked for parametric and non-parametric distributions, respectively. Predicted linear relationships were analyzed by least-squares linear regression and the derived slopes and y-axis intercepts compared. Comparisons between multiple paired groups were conducted using repeated-measures analysis of variance (ANOVA). Examination of the influence of two different independent variables (e.g. genotype and dose) on one dependent variable was conducted using two-way ANOVA. Graphical presentation and statistical analyses were conducted using Prism v5 (GraphPad Software). Values were expressed as mean ± standard error of the mean. A p-value of less than 0.05 was considered statistically significant. Multiple testing was further adjusted for by Bonferroni correction.

Supplementary Material

Supplementary Material is available at HMG online.

Conflict of Interest statement. None declared.

Click here for additional data file.