PTPRG is an ischemia risk locus essential for HCO3–-dependent regulation of endothelial function and tissue perfusion
Abstract
Acid-base conditions modify artery tone and tissue perfusion but the involved vascular-sensing mechanisms and disease consequences remain unclear. We experimentally investigated transgenic mice and performed genetic studies in a UK-based human cohort. We show that endothelial cells express the putative HCO3–-sensor receptor-type tyrosine-protein phosphatase RPTPγ, which enhances endothelial intracellular Ca2+-responses in resistance arteries and facilitates endothelium-dependent vasorelaxation only when CO2/HCO3– is present. Consistent with waning RPTPγ-dependent vasorelaxation at low [HCO3–], RPTPγ limits increases in cerebral perfusion during neuronal activity and augments decreases in cerebral perfusion during hyperventilation. RPTPγ does not influence resting blood pressure but amplifies hyperventilation-induced blood pressure elevations. Loss-of-function variants in PTPRG, encoding RPTPγ, are associated with increased risk of cerebral infarction, heart attack, and reduced cardiac ejection fraction. We conclude that PTPRG is an ischemia susceptibility locus; and RPTPγ-dependent sensing of HCO3– adjusts endothelium-mediated vasorelaxation, microvascular perfusion, and blood pressure during acid-base disturbances and altered tissue metabolism.
eLife digest
Restricted blood flow in the heart or brain can deprive these vital organs of oxygen, thereby causing a heart attack or stroke. However, the body has compensatory mechanisms to mitigate damage: if the blood flow is reduced in one blood vessel, acidic waste accumulates locally. This causes nearby blood vessels to widen and increase the oxygen supply. Although scientists first observed this process 140 years ago, they have not yet devised a way to use it for treatment of heart attack or stroke.
Now, Hansen et al. discovered that a protein called RPTPγ, which is found on the lining of blood vessels, could be a good target for drugs intended to reduce the consequences of heart attacks and strokes. The protein RPTPγ has a similar structure to other proteins that bind bicarbonate, an important ion that buffers acids in the body. RPTPγ can also trigger signals to nearby cells, which suggests that the protein can monitor bicarbonate levels in the blood and tissue and alert blood vessels of the need to widen.
Hansen et al. found that the blood vessels of mice that lacked RPTPγ were unable to widen when needed. Moreover, mice without RPTPγ experienced abnormal changes in blood pressure and blood flow to the brain when oxygen consumption was elevated or pH was disrupted. Hansen et al. further analyzed genetic and health data from nearly 50,000 individuals in the UK Biobank. These analyses revealed that people with genetic changes that render RPTPγ ineffective are at higher risk of having a heart attack or stroke. People with these genetic variants also have reduced heart pumping ability.
The experiments suggest that a lack of functional RPTPγ affects an individual’s ability to adjust local blood flow in response to acid-base disturbances and oxygen deficits, increasing the risk of a heart attack or stroke. This information may help scientists develop new ways to prevent or treat heart attacks and strokes, or to treat other conditions like cancer, where pH is disturbed.
Introduction
Inadequate tissue perfusion relative to metabolic demand is a fundamental pathophysiological cause of acute (e.g. myocardial infarction and stroke) and chronic (e.g. heart failure and neurodegenerative disorders) cardiovascular disease. Newer treatment options improve the prognosis for some of these debilitating conditions but we need alternative therapeutic strategies in order to minimize ischemia-related morbidity and mortality.
Sensing of the chemical environment in the wall of arteries coordinates local blood flow to meet the oxidative metabolic demand and counteract ischemia. Despite their obvious clinical importance, the cellular and molecular mechanisms responsible for sensing metabolic disturbances and responding to inadequate perfusion are not well understood (Boedtkjer, 2018). Local acidification is an important signal of insufficient nutrient delivery and waste product elimination. Although functional effects of H+ in the vascular wall have long been appreciated (Boedtkjer, 2018; Gaskell, 1880; Boedtkjer et al., 2016a; Boedtkjer and Aalkjaer, 2012; Boedtkjer and Aalkjaer, 2013a; Boedtkjer and Aalkjaer, 2013b), the consequences of associated changes in CO2/HCO3– buffer composition have received far less attention. CO2 and HCO3– constitute an important buffer pair that minimizes acute changes in pH (Roos and Boron, 1981), improves spatial H+ mobility (Spitzer et al., 2002; Boedtkjer et al., 2016b), and provides substrate for acid-base transporters in cell membranes (Aalkjaer and Hughes, 1991; Boedtkjer et al., 2006). Evidence supports that cells also possess sensors that respond to changes in [HCO3–] in the intracellular and extracellular space (Zhou et al., 2016; Boedtkjer et al., 2016c; Chen et al., 2000). Interventions—such as NaHCO3 supplementation—that modify systemic acid-base status are available (Voss et al., 2020) but expected to carry considerable adverse effects and be too general for cardiovascular therapy. The identification of HCO3–-related proteins (e.g. carbonic anhydrases, Na+,HCO3–-cotransporters, Cl–/HCO3–-exchangers, and HCO3–-sensors) opens largely unexplored avenues for pharmacological treatment. Because the vasculature of vital organs—particularly the brain and heart—responds to metabolic deregulation and is sensitive to acid-base disturbances, it is likely that therapy targeting HCO3–-related proteins will be able to modify myocardial and cerebral perfusion preferentially in areas of unmet metabolic demand.
Through mechanisms that require the transmembrane Receptor Protein Tyrosine Phosphatase (RPTP)γ, isolated decreases in extracellular [HCO3–] enhance HCO3– reabsorption in the renal proximal tubule (Zhou et al., 2016) and contractions of basilar arteries (Boedtkjer et al., 2016c) when pCO2 and pH are kept constant using out-of-equilibrium technology (Zhao et al., 1995). The extracellular domain of RPTPγ resembles the active site of the carbonic anhydrases (Barnea et al., 1993) and therefore likely binds HCO3–. However, the carbonic anhydrase-like domain of RPTPγ lacks the histidine residues considered essential for catalyzing equilibration of the reaction: (Zhou et al., 2016). Instead, HCO3– may induce a dimerization-dependent auto-inhibitory response in RPTPγ and thereby regulate signaling via the intracellular tyrosine phosphatase domains (Barr et al., 2009). Apart from these structural observations, we do not currently understand the cellular and molecular mechanisms for impact of HCO3– and RPTPγ on resistance artery function, blood pressure, and control of tissue perfusion. Although CO2 has no direct net effect on cerebrovascular tone (Boedtkjer et al., 2016c), respiratory changes leading to hyper- or hypocapnia will lead to secondary rises or falls in [HCO3–] and pH that in turn influence artery contractions and tissue perfusion.
In the present study, we demonstrate—based on experimental investigations in transgenic mice—that RPTPγ (a) facilitates endothelium-dependent relaxation of resistance arteries through mechanisms regulated by extracellular HCO3– and (b) adjusts microvascular perfusion and blood pressure during increased tissue metabolism and acid-base disturbances. Through translational genetic studies in the UK Biobank cohort (Bycroft et al., 2018), we further substantiate the importance of the identified mechanisms for human cardiovascular health as we demonstrate that predicted loss-of-function variants in PTPRG, encoding RPTPγ, are associated with human ischemic vascular disease in the heart and brain.
Results
Ptprg is widely expressed in the vascular endothelium
We first evaluated promoter activity for Ptprg by β-galactosidase staining of mice with a promoterless LacZ insert under transcriptional control of the Ptprg promoter. We found signs of prominent Ptprg transcriptional activity in the endothelium of mouse basilar, middle cerebral, and coronary arteries (Figure 1A). We also saw signs of Ptprg transcriptional activity in pulmonary and skeletal muscle arteries, whereas transcription of Ptprg appeared lower in the aorta and mesenteric arteries (Figure 1A).
Since promoter activity is no direct measure of expression—and cytosolic dilution of the chromogenic reaction product lowers the sensitivity of the β-galactosidase reporter assay—we followed up with quantitative RT-PCR analyses that confirmed substantial steady-state levels of Ptprg mRNA in basilar, middle cerebral, coronary, and pulmonary arteries and identified clearly detectable levels in the aorta, mesenteric, and skeletal muscle arteries (Figure 1B).
We next explored the endothelial expression pattern of Ptprg based on single-cell transcriptomic data from healthy mice (Kalucka et al., 2020). These data verify Ptprg expression in the vascular endothelium of the brain, heart, lung, skeletal muscle, and small intestine (Figure 1C,D). The data further show that the endothelial Ptprg expression extends from the arterial vasculature, across capillaries, and into the veins (Figure 1D,E).
Taken together, we identify widespread endothelial Ptprg expression with prominent levels in the cerebral and coronary vasculature where metabolites, including acid-base equivalents, strongly influence arterial tone (Gaskell, 1880; Boedtkjer et al., 2016c; Wahl et al., 1970).
RPTPγ enhances endothelium-dependent vasorelaxation
Isolated small arteries pre-contracted with the thromboxane A2 analog U46619 relaxed concentration-dependently in response to the classical endothelium-dependent agonist acetylcholine (Figure 2). In order to evaluate a distinct additional endothelial signaling pathway (D'Andrea et al., 1998) capable of inducing potent endothelium-dependent vasorelaxation of cerebral and mesenteric small arteries (Boedtkjer et al., 2013a; McNeish et al., 2005; Bucci et al., 2005), we also tested the effect of the proteinase-activated receptor (PAR)2 agonist SLIGRL-amide (Figure 3).
Under physiological conditions with CO2/HCO3– present in the buffer solutions, acetylcholine- and SLIGRL-amide-induced vasorelaxation of basilar (Figures 2A–C and 3A–C), mesenteric (Figures 2I and 3I), and coronary (Figures 2M and 3L) arteries from RPTPγ knockout (KO) mice was reduced compared to arteries from wild type (WT) mice. We observed this difference between arteries from RPTPγ KO and WT mice whether experiments in presence of CO2/HCO3– were conducted with (Figures 2B and 3B) or without (Figures 2C and 3C) the non-physiological buffer HEPES. Addition of an artificial buffer is necessary for pH control in absence of CO2/HCO3–; and we included HEPES also in the CO2/HCO3–-containing solutions in order to separate the vascular effects of CO2/HCO3– omission from potential effects of adding HEPES (Altura et al., 1980). Acetylcholine-induced vasorelaxation of basilar arteries from RPTPγ KO mice was suppressed in presence of CO2/HCO3– mostly because of attenuated NO-dependent (i.e. N-nitro-L-arginine methyl ester (L-NAME)-sensitive) signaling (Figure 2D–G). In contrast, the diminished vasorelaxation of basilar arteries in response to SLIGRL-amide (Figure 3D–F) and of mesenteric and coronary arteries in response to acetylcholine (Figure 2J–L and Figure 2N–P) and SLIGRL-amide (Figure 3J,K,M,N) was mostly explained by smaller endothelium-dependent hyperpolarization (EDH)-type responses that were L-NAME-insensitive but inhibited by apamin and TRAM-34. Consistent with these findings, L-NAME-insensitive smooth muscle hyperpolarization was significantly attenuated in basilar arteries from RPTPγ KO mice when elicited by SLIGRL-amide (Figure 3G,H) but not when elicited by acetylcholine (Figure 2H).
The vasorelaxant influence of RPTPγ depends on CO2/HCO3–
In order to evaluate whether RPTPγ-dependent signaling is influenced by HCO3– in the vascular wall, we next repeated studies of endothelium-dependent vasorelaxation using nominally CO2/HCO3–-free, HEPES-buffered solutions (Figure 4). Under these conditions, we saw no RPTPγ-dependent differences in overall vasorelaxant function of basilar or mesenteric arteries in response to acetylcholine (Figure 4A,B,G) or SLIGRL-amide (Figure 4K,L,P). The relative contribution of NO-dependent signaling and EDH-type responses to vasorelaxation was also not significantly affected by RPTPγ KO in basilar (Figure 4C–F) or mesenteric (Figure 4H–J) arteries in response to acetylcholine. Whereas the dependency on NO-signaling vs EDH-type vasorelaxation in response to SLIGRL-amide was likewise not affected in mesenteric arteries (Figure 4Q,R), we saw indications of an apparent relative increase in NO-mediated vasorelaxation when basilar arteries from RPTPγ KO mice were stimulated with SLIGRL-amide in absence of CO2/HCO3– (Figure 4M–O). This unexpected change in the relative importance of the underlying vasorelaxant mechanisms in basilar arteries could be a compensation secondary to the predominant inhibition of NO-dependent vasorelaxation in these arteries in presence of CO2/HCO3– (Figure 2F) and hence presumably in vivo.
NO-synthase expression and responses to NO donors are largely independent of RPTPγ
We next explored the mechanism whereby RPTPγ influences resistance artery function (Figure 5).
We first evaluated whether the difference in vasorelaxation between arteries from RPTPγ KO and WT mice depends on a change in expression of the endothelial NO-synthase (eNOS). Based on immunoblotting of basilar (Figure 5A) and mesenteric (Figure 5D) arteries, we found equal eNOS expression between arteries from RPTPγ KO and WT mice.
We also evaluated whether the smooth muscle vasorelaxant responses to exogenous NO differ between arteries from WT and RPTPγ KO mice. In presence of CO2/HCO3–, spermine NONOate and S-nitroso-N-acetylpenicillamine (SNAP) produced equivalent vasorelaxation in basilar and mesenteric arteries from RPTPγ KO and WT mice (Figure 5B,C,E,F).
Together, our findings support that RPTPγ exerts acute influences on endothelial signaling pathways regulating the activity of endothelium-dependent vasorelaxation in response to [HCO3–]. This interpretation is in agreement with the near-complete normalization of vasorelaxation observed in arteries from RPTPγ KO mice upon omission of CO2/HCO3– (Figure 4).
RPTPγ enhances endothelial Ca2+ responses with no effect on intracellular pH
We then determined whether endothelial intracellular Ca2+ signals differ between arteries from RPTPγ KO and WT mice (Figure 5G–N). Loading of isolated arteries with Ca2+-sensitive fluorophores through luminal perfusion ensured endothelium-specific fluorescence signals (Figure 5G,K).
In the presence of CO2/HCO3–, application of acetylcholine elevated intracellular [Ca2+] substantially more in endothelial cells of basilar and mesenteric arteries from WT mice than in arteries from RPTPγ KO mice (Figure 5H,J,L,N). In contrast, we saw no significant difference between the intracellular Ca2+ responses of endothelial cells in arteries from WT and RPTPγ KO mice when they were investigated in the nominal absence of CO2/HCO3– (Figure 5I,J,M,N).
Reduced endothelial intracellular Ca2+ responses in arteries from RPTPγ KO mice in presence of CO2/HCO3– likely explain the observed mixed attenuation of NO-mediated and EDH-type vasorelaxation (Figures 2 and 3) as these responses rely on the Ca2+-sensitive NO-synthase and on intermediate- (IKCa) and small-conductance (SKCa) Ca2+-activated K+-channels, respectively (Figure 5O).
Intracellular H+ and Ca2+ can compete for buffer binding (Batlle et al., 1993), many Ca2+ handling proteins are pH-sensitive (Boedtkjer and Aalkjaer, 2012), and pH is an important regulator of NO-synthase activity (Boedtkjer et al., 2011). Therefore, we next evaluated whether RPTPγ KO influences intracellular pH—measured using the pH-sensitive fluorophore BCECF—in endothelial cells of small arteries. In basilar (Figure 5P) as well as mesenteric (Figure 5Q) arteries, intracellular pH of the endothelial cells was similar for RPTPγ KO and WT mice, irrespective of whether experiments were performed with or without CO2/HCO3–. Omission of CO2/HCO3– from the bath solutions had no net effect on endothelial steady-state intracellular pH (Figure 5P,Q), which is consistent with previous studies based on arteries from C57BL/6 mice (Boedtkjer et al., 2011; Boedtkjer et al., 2012).
RPTPγ amplifies hyperventilation-induced blood pressure elevations
If the vasomotor impact of RPTPγ is sufficiently widespread to affect total peripheral resistance, we predicted that KO of RPTPγ would influence blood pressure regulation. We therefore measured systemic blood pressure and heart rate by radiotelemetry (Figure 6A,B). Based on these recordings, we observed no difference in resting blood pressure or heart rate (Figure 6A) or in the blood pressure-elevating effect of L-NAME (Figure 6B) between RPTPγ KO and WT mice.
As the lack of difference in resting blood pressure between WT and RPTPγ mice could be due to compensatory adaptations, we next studied acute blood pressure responses of endotracheally intubated, mechanically ventilated mice under capnographic control. When mice were ventilated to a normal expiratory end-tidal CO2 fraction (FETCO2) of 3.5%, arterial blood gas parameters were very similar in RPTPγ KO and WT mice (Table 1) consistent with previous findings (Zhou et al., 2016). Hyperventilating mice until FETCO2 was lowered to 2% caused the expected decrease in PaCO2 and increase in pHa (Table 1). The shift of the chemical equilibrium also resulted in the anticipated decrease in [HCO3–]a (Table 1). The mean arterial blood pressure of WT mice increased around 20% during hyperventilation and this blood pressure response was reduced by approximately 1/3 in RPTPγ KO mice (Figure 6C). Although the observed hyperventilation-induced reduction in arterial [HCO3–] was of modest magnitude (1.4 mM, Table 1), the associated increase in blood pressure is consistent with attenuated RPTPγ-dependent vasodilator influences when extracellular [HCO3–] decreases (Boedtkjer et al., 2016c). The remaining hyperventilation-induced blood pressure elevation likely results from a direct vasocontractile effect of alkalosis (Boedtkjer et al., 2016c).
RPTPγ regulates cerebral perfusion
Metabolic waste products—that are generated at elevated rate and locally accumulate when neuronal activity increases—influence cerebral blood flow. We studied cerebral perfusion using laser speckle imaging (Figure 6D–F) in order to evaluate the involvement of RPTPγ in sensing of acid-base disturbances. Tissue perfusion in the somatosensory barrel cortex increased twice as much in RPTPγ KO mice compared to WT mice during whisker stimulation (Figure 6E). This finding is consistent with metabolically produced H+ leading to a local decrease in [HCO3–] that is sensed via RPTPγ in arteries from WT mice (Figure 5O). As RPTPγ-dependent vasorelaxation wanes at low [HCO3–] (Boedtkjer et al., 2016c), our data support that RPTPγ limits increases in cerebral perfusion induced by metabolic acidosis (Figure 6E).
Acid-base conditions during increased neuronal activity are complicated, as the increased metabolism will also induce some degree of CO2 accumulation that tends to increase [HCO3–]. In order to produce a simpler acid-base disturbance with predictable changes in extracellular [HCO3–], we next hyperventilated mice, which reduced pCO2, increased pH, and lowered [HCO3–] (Table 1). Under these conditions, cerebral vasoconstriction is expected in WT mice based on the combined effect of alkalosis and decreased [HCO3–] (Boedtkjer et al., 2016c). If RPTPγ is required for sensing of HCO3–, only the effect of the elevated pH is expected in the RPTPγ KO mice; and indeed, we saw a 30% smaller drop in cerebral perfusion in RPTPγ KO mice compared to WT mice during hyperventilation (Figure 6F). These recordings may in fact underestimate the difference in cerebrovascular resistance between RPTPγ KO and WT mice since the cerebral perfusion pressure—based on the blood pressure changes shown in (Figure 6C)—increased less in hyperventilated RPTPγ KO than WT mice. When whisker stimulation was performed during hyperventilation, we saw an even more prominent difference in cerebral perfusion between WT and RPTPγ KO mice (Figure 6F).
Together, our findings demonstrate that RPTPγ plays a substantial role for control of blood pressure (Figure 6C) and cerebral perfusion (Figure 6E,F) during acid-base deviations where [HCO3–] decreases below the normal level.
PTPRG is a susceptibility locus for ischemic vascular disease
Based on human exome sequencing data from the UK Biobank, we next explored whether the functional role of RPTPγ in regulating endothelial function and responding to metabolic disturbances in murine cerebral, mesenteric, and coronary arteries is corroborated by and translates to an altered risk for human ischemic vascular disease amongst carriers of predicted loss-of-function variants within PTPRG. We identified 334 missense variants and 75 predicted loss-of-function variants (categorized as low, moderate, and high impact) in PTPRG among the ~50,000 UK Biobank participants with available exome sequencing data.
The 72 carriers of predicted loss-of-function PTPRG variants with moderate or high impact showed substantially elevated risk (Figure 7A) of a combined vascular disease phenotype encompassing the three vascular beds functionally investigated in this study (Figures 2–6): cerebral infarct (ICD-10 diagnosis code I63), angina pectoris (ICD-10 diagnosis code I20), acute myocardial infarction (ICD-10 diagnosis code I21), and acute vascular disorders of the intestine (ICD-10 diagnosis code K55.0). In contrast, the 123 carriers of predicted loss-of-function PTPRG variants with low impact as well as the 29,975 carriers of one or more missense PTPRG variants had a risk similar to non-carriers for this aggregated phenotype of acute vascular disease (Figure 7A).
We followed-up by investigating the separate occurrences of diagnosed angina pectoris, acute myocardial infarction, and cerebral infarction amongst the individuals, who carried a predicted loss-of-function variant with moderate or high impact. The carriers of these most severe loss-of-function variants had dramatically increased risk (~7 fold) of cerebral infarction whereas the risk of acute myocardial infarction showed a tendency toward elevation that did not reach statistical significance (Figure 7B). The number of exome-sequenced UK Biobank participants with acute intestinal vascular disease was too low to provide statistical power for a separate analysis.
To complement the evidence based on ICD-10 diagnosis codes, we further evaluated whether carriers of missense and loss-of-function variants in PTPRG differed in self-reported diseases compared to non-carriers (Figure 7C). The risk of self-reported heart attacks was greatly elevated (~4 fold) amongst individuals carrying loss-of-function PTPRG variants of moderate and high impact (Figure 7C). The risk of self-reported stroke also showed a strong tendency toward elevation—with borderline statistical significance—amongst the carriers of moderate- and high-impact loss-of-function PTPRG variants (Figure 7C).
As cardiac pumping function often deteriorates following coronary ischemia, we next evaluated the cardiac contractile function in PTPRG loss-of-function carriers amongst the 12,851 exome-sequenced UK Biobank participants with available cardiac magnetic resonance imaging (MRI) data. We observed a significantly lower left ventricular ejection fraction in carriers of high-impact loss-of-function PTPRG variants compared to non-carriers (Figure 7D) consistent with the greater risk of ischemic heart disease (Figure 7C).
The statistical analyses linking PTPRG to risk of ischemic vascular disease (Figure 7A–C) and reduced left ventricular ejection fraction (Figure 7D) were all adjusted for sex, age, body mass index, genetic principal component, smoking status, dyslipidemia, diabetes, and hypertension between carriers and non-carriers of the evaluated PTPRG variants.
Taken together, we demonstrate that RPTPγ (a) CO2/HCO3–-dependently enhances endothelial intracellular Ca2+ responses and endothelium-dependent vasorelaxation, (b) regulates microvascular perfusion and blood pressure during acid-base disturbances, and (c) is associated with human ischemic vascular disease of the brain and heart.
Discussion
Local acid-base-dependent mechanisms of arterial tone regulation match tissue perfusion to the oxidative metabolic demand and thereby dynamically control cardiovascular function during cycles of intermittent rest and activity. Re-establishing proper tissue perfusion is of obvious clinical importance when managing patients with latent or fulminant ischemia. Responses of coronary and cerebral arteries to acid-base disturbances have been recognized for almost 140 years (Gaskell, 1880; Lassen, 1968) but this information has not yet been harnessed for therapeutic intervention for lack of identified molecular players and their interactions. In the current study, we identify PTPRG as a susceptibility locus for human ischemic vascular disease (Figure 7) and provide mechanistic evidence that RPTPγ regulates endothelial intracellular Ca2+ responses (Figure 5H,L), endothelium-dependent vasorelaxation (Figures 2 and 3), cerebral perfusion (Figure 6E,F), and blood pressure (Figure 6C) during acid-base disturbances and increased metabolic demand.
Variation in buffer composition accompanies changes in pH during acid-base disturbances. We show that the CO2/HCO3– buffer plays a hitherto unappreciated role for regulation of endothelial function and arterial tone through mechanisms that require RPTPγ (Figures 2–4). As schematically illustrated in Figure 5O, we demonstrate that RPTPγ enhances endothelial intracellular Ca2+ signals (Figure 5G–N) that in turn activate endothelial NO synthesis and EDH-type responses (Figures 2 and 3). These effects are in congruence with the Ca2+-sensitivity of eNOS and the SK and IK Ca2+-activated K+-channels. The vasorelaxant influence of RPTPγ wanes at low extracellular [HCO3–] (Boedtkjer et al., 2016c), which attenuates by half the elevation of cerebral perfusion during sensory input (Figure 6E) and amplifies by 30% the decrease in perfusion during hyperventilation (Figure 6F). We also observe that RPTPγ is necessary for approximately 1/3 of the blood pressure increase during hyperventilation (Table 1 and Figure 6C). Considering the prominent influence of RPTPγ on endothelial function in several vascular beds, somewhat surprisingly, we observe blood pressure consequences of RPTPγ KO only when we impose acute acid-base disturbances (Figure 6A–C). The unaltered resting blood pressure in RPTPγ KO mice (Figure 6A) most likely reflects that numerous discrete mechanisms are involved in blood pressure control and that compensation—for example, through nervous, hormonal or renal influences (Boedtkjer and Aalkjaer, 2013b)—can maintain blood pressure in the sustained phase and mask hemodynamic consequences of RPTPγ except when acute acid-base disturbances are imposed.
Previously recognized cellular acid-base sensors include G-protein coupled receptors (e.g. OGR1, GPR4, and TDAG8), ion channels (e.g. ASIC, BKCa), and enzymes (e.g. NO-synthase and rho-kinase) sensitive to H+ (Boedtkjer et al., 2011; Boedtkjer et al., 2012; Fleming et al., 1994; Schubert et al., 2001; Ludwig et al., 2003; Wemmie et al., 2006; Wenzel et al., 2020). In addition, we and others have described cellular functions responsive to changes in [HCO3–] and modified by RPTPγ (Zhou et al., 2016; Boedtkjer et al., 2016c) or the soluble adenylyl cyclase (Chen et al., 2000). Evidence connecting acid-base sensors to human disease has so far been scarce; but in the current study, we identify PTPRG as an ischemia susceptibility locus (Figure 7), which is supported by a recent meta-analysis of genome-wide association studies linking PTPRG to ischemic stroke in African Americans (Carty et al., 2015).
Using out-of-equilibrium technology—that permits separate control of pH and the individual CO2/HCO3– buffer components—we previously demonstrated that RPTPγ responds to changes in extracellular [HCO3–] independently of pCO2 and pH (Boedtkjer et al., 2016c). The homology between the extracellular domain of RPTPγ and the active site of the carbonic anhydrases supports the ability of RPTPγ to bind HCO3– even though it lacks histidine residues required for carbonic anhydrase activity (Zhou et al., 2016). The intracellular aspect of RPTPγ contains phosphatase domains with suggested auto-inhibitory activity; and HCO3– possibly alters the phosphatase activity by influencing the degree of RPTPγ dimerization (Barr et al., 2009). The separate signaling effects of H+ and HCO3– are further confirmed by the similar steady-state intracellular pH of endothelial cells in arteries from WT and RPTPγ KO mice (Figure 5P,Q).
The current study adds to our growing appreciation that acid-base equivalents, and HCO3– in particular, fulfill multifaceted functions (Boedtkjer et al., 2016a). In cerebral resistance arteries, HCO3– (a) contributes to buffering of acute acid loads (Rasmussen and Boedtkjer, 2018), (b) serves as substrate particularly for the Na+,HCO3–-cotransporter NBCn1 that protects against intracellular acidification (Thomsen et al., 2014), and (c) is sensed by RPTPγ to regulate cerebral perfusion (Figure 6). In migrating vascular smooth muscle cells from conduit arteries, the high spatial mobility of the CO2/HCO3– buffer system also contributes to dissipating local pH gradients in diffusion-restricted spaces of filopodia (Boedtkjer et al., 2016b).
Our current (Figures 1–4 and 6) and previous (Boedtkjer et al., 2016c) findings demonstrate that the endothelium of resistance arteries senses the local acid-base composition and has capacity to modify vascular resistance and perfusion in response to disturbances in extracellular [HCO3–]. Earlier studies of mouse arteries provide additional evidence that intracellular acidification inhibits endothelial NO-synthesis (Boedtkjer et al., 2011; Boedtkjer et al., 2012; Thomsen et al., 2014), intracellular alkalinisation interferes with myo-endothelial current transfer required for EDH-type responses (Boedtkjer et al., 2013a), and acid-base disturbances modify prostanoid-mediated endothelium-dependent vasocontraction (Baretella et al., 2014). Together, these studies highlight the sophisticated nature of vasomotor control that integrate the requirement for local blood flow and the necessity to minimize deviations in capillary pressure and fluid filtration during local metabolic disturbances, perturbed nutrient delivery, and restricted waste product elimination. Acting as a brake on vasodilation in regions of unmet metabolic demand, RPTPγ could reduce the degree of edema and consequent tissue damage caused by unopposed H+-induced vasodilation (Boedtkjer, 2018; Boedtkjer et al., 2016a).
As illustrated in Figure 1D and E, Ptprg mRNA is found in arteriolar, capillary, as well as venular endothelial cells (Kalucka et al., 2020; Vanlandewijck et al., 2018), and it remains a future task to evaluate physiological and pathophysiological roles of RPTPγ in these individual blood vessel segments. Transcript levels for Ptprg appear particularly high in a smaller subset of endothelial cells (up to 20–40%, depending on the vascular bed; Figure 1C–E), which raises the intriguing possibility that a distinct population of endothelial cells act as the primary sensors of the metabolic environment.
Hyperventilation-induced cerebral vasoconstriction can cause syncope but is also a therapeutic tool for rapidly reducing cerebral blood flow, capillary filtration, and intracranial pressure. Hypocapnia during hyperventilation leads to decreased [HCO3–] and alkalosis (Table 1); and the current study demonstrates that RPTPγ is necessary for approximately 30% of the associated decrease in cerebral perfusion (Figure 6F).
In conclusion, the single-pass transmembrane protein RPTPγ enhances endothelium-dependent vasorelaxation of resistance arteries—by amplifying endothelial intracellular Ca2+ responses, NO synthesis, and EDH-type responses—and adjusts cerebral perfusion during increased neuronal activity and acid-base disturbances. Although RPTPγ shows no major impact on resting blood pressure, it contributes markedly to the blood pressure increase observed during hyperventilation. Supporting the translational value of our findings, predicted loss-of-function variants in PTPRG are associated with human ischemic vascular disease in the brain and heart. The influence of RPTPγ on vascular resistance specifically in tissues with disturbed acid-base composition makes RPTPγ a promising focus for targeted therapeutic intervention against ischemia.
Materials and methods
RPTPγ KO mice (MGI:3663183) from Dr. Joseph Schlessinger, Yale University, USA (Lamprianou et al., 2006), were kindly provided by Dr. Walter F. Boron, Case Western Reserve University, USA (Zhou et al., 2016). The transgenic mice were originally generated by inserting an internal ribosome entry site–β-galactosidase–neomycin–poly(A) cassette immediately 3’ of exon 1 (Lamprianou et al., 2006). Homologous recombination in W4 129S6/SvEv mouse embryonic stem cells (RRID:CVCL_Y634) was achieved following electroporation, and correctly targeted embryonic stem cells were used for embryo aggregation to achieve germ line chimeras (Lamprianou et al., 2006). At the time of experiments, the transgene-carrying mice had been backcrossed to the C57BL/6 line for at least eight generations. In total, we used 125 WT and 112 RPTPγ KO mice that were kept on a 12 hr light/12 hr dark cycle with ad libitum access to tap water and chow (Altromin 1324; Brogaarden, Danmark). We used Eco-Pure Aspen Chips 2 Premium cage bedding (Datesand, UK). Genomic DNA used for genotyping was isolated from tail biopsies collected at weening and again after sacrifice. All animal experimental procedures were approved by the Danish Animal Experiments Inspectorate (2016-15-0201-00982).
Small vessel myography
Request a detailed protocolMale RPTPγ KO and matched WT mice (>8 weeks old) deeply anesthetized by intraperitoneal injection of pentobarbital were sacrificed by exsanguination. Approximately 2 mm long segments of basilar, mesenteric, and coronary septal arteries were dissected under a stereomicroscope, mounted in four-channel wire myographs (610M; DMT, Denmark) for isometric investigation, heated to 37°C, and normalized to 90% of the internal diameter corresponding to a transmural pressure of 100 mmHg (Mulvany and Halpern, 1977). The myograph chambers were aerated with 5% CO2/balance air (for CO2/HCO3–-containing physiological saline solutions) or nominally CO2-free air (for CO2/HCO3–-free solutions). The CO2/HCO3–-containing solution consisted of (in mM): 119 NaCl, 22 NaHCO3, 10 HEPES, 1.2 MgSO4, 2.82 KCl, 5.5 glucose, 1.18 KH2PO4, 0.03 EDTA, 1.6 CaCl2. The CO2/HCO3–-free solution was produced by substituting NaHCO3 with NaCl. HEPES is necessary for pH control in the CO2/HCO3–-free solution; and we included HEPES also in the CO2/HCO3–-containing solution in order to separate vascular effects of CO2/HCO3– omission from potential effects of adding HEPES (Altura et al., 1980). For the control experiments in Figures 2C and 3C, we used a CO2/HCO3–-containing solution without HEPES, i.e., where the 10 mM HEPES was substituted with 5 mM NaCl. After all contents were added, the solutions were vigorously bobbled with the appropriate CO2-containing or nominally CO2-free gas mixture at 37°C before pH was adjusted to 7.40 while allowing ample time for buffer equilibration. Arterial force was acquired using a PowerLab 4/25 recorder and LabChart 7 Pro software (RRID:SCR_001620; ADInstruments, New Zealand). As standard warm-up protocol, arteries were exposed twice to 3 µM of the thromboxane A2 analog U46619 for 2 min, and arteries that produced less than 1 mN force were excluded from further analysis. Vasorelaxation in response to acetylcholine, the PAR2 agonist SLIGRL-amide, and the NO-donors SNAP, and spermine NONOate was tested in U46619-contracted arteries that developed stable tension equivalent to ~70% of the maximal response during the standard warm-up procedure. Vasorelaxation, from the pre-contraction level, was quantified as the relative decrease in active tone during the last 30 s of each 2 min agonist application. We used two NO donors that differ, for instance, in their rate and location of NO release (Bradley and Steinert, 2015). Because acetylcholine-induced vasorelaxation shows substantial tachyphylaxis in arteries from C57BL/6 mice, we tested single doses of acetylcholine applied on separate U46619 pre-contractions. In contrast, SLIGRL-amide does not show this degree of tachyphylaxis, and we therefore tested its vasorelaxant properties through cumulative additions.
Endothelial Ca2+ measurements
Request a detailed protocolCa2+-sensitive fluorophores were loaded preferentially into endothelial cells by slow (~1 mL/hr for 45 min) perfusion of physiological saline solution containing 2.5 µM Fluo-4 and 3.5 µM Fura Red (Invitrogen) through the lumen of basilar and mesenteric arteries mounted in a pressure myograph (120CP; DMT) (Boedtkjer et al., 2011). The cell permeant acetoxymethyl ester (AM) forms of the fluorophores were first solubilized in a load mix containing dimethyl sulfoxide, Pluronic F127, and Cremophor EL. The loaded arteries were next mounted in a confocal wire myograph (360CW; DMT) and studied using a Zeiss Axiovert 200M confocal microscope equipped with an LSM Pascal exciter and a 40 × objective (LD C-Apochromat; N.A. 1.10; Zeiss, Germany). The arteries were excited at 488 nm and emission light collected at wavelengths in the range of 505–530 nm (F505-530; representing Fluo-4 signals that increase at elevated [Ca2+]) and longer than 600 nm (F>600; representing Fura Red signals that decrease at elevated [Ca2+]). Arteries with no or only solitary loaded endothelial cells were excluded from analysis. The F505-530/F>600 ratio normalized to the baseline ratio (the first 2 min of the recording) was used to evaluate relative intracellular Ca2+ dynamics.
Endothelial intracellular pH measurements
Request a detailed protocolInverted basilar and mesenteric arteries were mounted on a 40 µm wire in a confocal wire myograph (360CW, DMT) and loaded with the pH-sensitive fluorophore 2',7'-bis-(2-carboxyethyl)−5-(and-6)-carboxyfluorescein (BCECF) in physiological saline solution (1 µM BCECF-AM in 0.05‰ DMSO) to preferentially load endothelial cells, essentially as previously described (Boedtkjer et al., 2011; Boedtkjer and Aalkjaer, 2009). The loaded arteries were studied using an Olympus IX83 microscope equipped with a 20 × objective (LUCPlanFLN; N.A. 0.45; Olympus, Japan) and an ORCA-Flash 4.0 camera (Hamamatsu, Japan). The arteries were alternatingly excited at 490 and 436 nm and emission light collected at 530 nm. The F490/F436 fluorescence ratio was calibrated to intracellular pH using the high-[K+] nigericin method (Boedtkjer and Aalkjaer, 2009; Aalkjaer and Cragoe, 1988), which in arteries agrees well with the null-point technique (Danielsen et al., 2013). Steady-state intracellular pH was recorded in the presence and absence of CO2/HCO3–.
Membrane potential recordings
Request a detailed protocolMembrane potentials were measured in vascular smooth muscle cells of isolated basilar arteries using sharp electrodes as previously described (Boedtkjer et al., 2013a). The arteries were mounted in a wire myograph (420A; DMT) and microelectrodes that had resistances of 40–120 MΩ when backfilled with 3 M KCl were inserted into the vascular wall from the adventitial side. Cell impalement was observed as a sudden drop in voltage followed by sharp return to baseline upon retraction. Measurements with more than 10 mV difference between the baseline recording before and after impalement were excluded from analysis.
LacZ reporter studies
Request a detailed protocolThe genetic insert that disrupts RPTPγ expression in the employed KO mice contains a promotorless LacZ sequence allowing for β-galactosidase expression under control of the Ptprg promotor (Lamprianou et al., 2006). Four homozygous KO mice (PtprgLacZ/LacZ) and two WT mice were perfusion fixed with 4% (weight/volume) paraformaldehyde in phosphate-buffered saline (PBS, in mM: 137 NaCl, 2.5 KCl, 4.3 Na2HPO4, and 1 KH2PO4) and investigated for promoter activity essentially as previously described (Boedtkjer et al., 2008). Segments of basilar, middle cerebral, gracilis, coronary, mesenteric, and pulmonary arteries as well as thoracic aorta were dissected free from surrounding tissue and washed in PBS overnight at 4°C. The arteries were then placed in staining solution (in mM: 5 K4Fe(CN)6, 5 K3Fe(CN)6, 2 MgCl2, 0.1% (weight/volume) sodium dodecyl sulfate (SDS), 0.1% (volume/volume) TWEEN−20, and 0.1% (weight/volume) 5-bromo-4-chloro-indolyl-13-D-galactoside (X-Gal)) for 24 hr at room temperature (~21°C). Finally, samples were transferred to PBS containing 1% (weight/volume) ethylenediaminetetraacetic acid (EDTA) and 4% paraformaldehyde in order to stop the staining reaction. Whole mount micrographs of the arteries were captured using a Leica M165 C stereomicroscope equipped with a Leica M170 HD camera (Germany). In addition, we paraffin-embedded and cut basilar and middle cerebral arteries—that showed the strongest staining—to 4-µm-thick histological sections that were visualized on an upright Leica DM light microscope equipped with a Leica DM300 digital camera in order to identify the cellular expression pattern.
Reverse transcription and quantitative polymerase chain reaction
Request a detailed protocolArteries dissected free from surrounding tissue in cold physiological saline solution were stored in RNALater (Qiagen, Denmark) at 4°C. The isolated arteries were homogenized in RLT lysis buffer (Qiagen) using a TissueLyser II (Qiagen) at 30 Hz for 2 min. RNA was isolated using the RNeasy Micro Qiacube kit including carrier RNA (Qiagen). Samples were reverse transcribed using random decamer primers and Superscript III Reverse Transcriptase (Invitrogen, Fisher Scientific, Denmark). Reactions without reverse transcriptase were performed in order to test for genomic amplification. Quantitative PCR was performed on an MX3000P system (Agilent, USA) based on Maxima Hot Start Taq DNA polymerase (ThermoFisher, Denmark). Ptprg mRNA levels relative to the reference genes Rn18s (18S ribosomal subunit) and Actb (β-actin) were evaluated in the different vascular beds based on the 2–ΔΔCT method (Livak and Schmittgen, 2001). We used the following forward (F) and reverse (R) primers and probes (P) purchased from Eurofins Genomics (Germany): Ptprg (F: 5’ TGG TTA CAA CAA AGC GAA AGC CT 3’, R: 5’ ATA CTG ATC ACA CTT TCT CCT TCC 3’, P: 5’ ATC TGG GAA CAA AAC ACG GGA ATC ATC AT 3’), Rn18s (F: 5’ AAT AGC CTT CGC CAT CAC TGC 3', R: 5’ GTG AGG TCG ATG TCT GCT TTC C 3’, P: 5’ TGG GGC GGA GAT ATG CTC ATG TGG TGT T 3’), and Actb (F: 5’ TGA CGT TGA CAT CCG TAA AG 3’, R: 5’ CTG GAA GGT GGA CAG TGA GG 3’ and P: 5’ AGT GCT GTC TGG TGG TAC CAC CAT GTA CC 3’). Probes were modified with 5’ 6-FAM and 3’ TAMRA. Each reaction consisted of 10 min at 95°C followed by 50 cycles of 30 s at 95°C, 60 s at 55°C, and 60 s at 72°C.
Single-cell RNA sequencing data
Request a detailed protocolWe explored levels of Ptprg transcripts in individual endothelial cells based on the EC Atlas database from VIB-KU Leuven, which is accessible at https://endotheliomics.shinyapps.io/ec_atlas/ (downloaded on 11 June 2020). This online database includes data from a recent single-cell RNA sequencing study on endothelial cells from healthy mice (Kalucka et al., 2020).
Immunoblotting
Request a detailed protocolWe determined protein expression levels of eNOS in basilar and mesenteric arteries by immunoblotting using previously described antibodies (Boedtkjer et al., 2011; Voss et al., 2019). Arteries were snap frozen in liquid nitrogen and then homogenized using pellet pestles (Sigma-Aldrich, Denmark) in a lysis buffer at pH 7.5 containing (in mM) 20 Tris-HCl, 150 NaCl, 5 ethylene glycol tetraacetic acid (EGTA), 10 NaF, 20 β-glycerophosphate sodium salt, and HALT protease and phosphatase inhibitor cocktail (Thermo Scientific, Denmark). Samples were sonicated for 45 s and centrifuged at ~16,000 g for 10 min. Total protein concentrations in the supernatants were measured using a bicinchoninic acid (BCA) protein assay kit (Thermo Scientific); and 10 µg total protein diluted in Laemmlie sample buffer (Biorad, Denmark) was loaded in each lane of an SDS polyacrylamide gel (Biorad). Membranes were first probed with anti-eNOS (RRID:AB_304967; 0.2 µg/mL ab5589; Abcam, UK) or anti-pan-actin (RRID:AB_2313904; 40 ng/mL #4968; Cell Signaling Technology, USA) primary antibody and then with secondary goat anti-rabbit antibody (RRID:AB_2099233; 30 ng/mL #7074; Cell Signaling Technology) conjugated to horseradish peroxidase. Bound antibody was detected by enhanced chemiluminescence (ECL Plus; GE Healthcare, Denmark) using an ImageQuant LAS 4000 luminescent image analyzer (GE Healthcare). Densitometric analyses were performed using ImageJ software (RRID:SCR_003070; Rasband; NIH, USA). Band densities of eNOS relative to pan-actin were reported after normalization to the average WT level.
Blood pressure measurements
Request a detailed protocolMice were anesthetized by subcutaneous injection of ketamine and xylazine (80 mg/kg Ketaminol vet and 8 mg/kg Narcoxyl vet; Intervet International, The Netherlands) and placed on a thermostatically controlled heating platform.
For telemetry-based measurements of resting blood pressure and effects of L-NAME ingestion, the catheter of a telemetry transmitter (HD-X11; Data Sciences International, USA) was inserted in the common carotid artery through a midline incision in the neck, and the transmitter body placed in a subcutaneous pocket during stereomicroscopy. Pain relief was achieved through subcutaneous injection of buprenorphine (0.2 mL/kg, Temgesic, Schering-Plough, Europe). Telemetry measurements started one week after the operation. One 48 hr long registration at baseline was followed by two 72 hr long registrations during which first 0.5 mg/mL and then 3 mg/mL L-NAME was added to the drinking water. Telemetry signals were recorded for 10 s every minute using Dataquest A.R.T. 4.3 and analyzed with Ponemah 5.0 software (RRID:SCR_017107; Data Sciences International). We evaluated daytime and nighttime blood pressure from 11 AM to 1 PM and 11 PM to 1 AM, respectively.
For measurements of acute blood pressure responses to hyperventilation, a catheter was inserted in the common carotid artery through a midline incision and connected to a pressure transducer (MLT0699, ADInstruments). The mice were intubated and ventilated on a Minivent type 845 ventilator (Harvard Apparatus, USA) at a frequency of 125 min−1 and with tidal volume adjusted until capnography readings (Capnograph Type 240, Hugo-Sachs Electronics, Germany) showed an FETCO2 of 3.5%. During experiments, hypocapnia was induced by elevating the ventilation until FETCO2 decreased to 2%. Mean arterial blood pressure was derived from the pressure traces using the blood pressure add-on for LabChart 8 Pro (ADInstruments). A few mice with an initial systolic blood pressure below 70 mmHg were excluded from the analysis.
Laser speckle imaging
Request a detailed protocolMice were initially anesthetized by subcutaneous injection of ketamine (80 mg/kg) and xylazine (8 mg/kg) followed by injection with 1/4 the initial dose every 45 min. After endotracheal intubation, mice were ventilated under capnography control as described above. The head of the anesthetized mouse was fixed in an adaptor for a stereotaxic frame (World Precision Instruments, UK) while the rest of the body was kept warm on a heating pad (Fine Science Tools Inc, Canada). The skin covering the top of the skull was removed, the surface of the bone cleaned, and a coverslip mounted in agarose in order to minimize optical reflections. Whiskers were fixed to a metal pole mounted on a cylindrical solenoid controlled by an Arduino circuit (Funduino Kit, Germany). The whiskers on one side of the head were first stimulated by moving whiskers vertically with an amplitude of 8 mm and frequency of 4 Hz. Then, the stimulation protocol was repeated on the opposite side of the head, and the two responses were averaged for each mouse. Whisker stimulation was performed under control conditions (FETCO2 = 3.5%) and after FETCO2 was reduced to 2% by mechanical hyperventilation.
Speckle images of 1088 × 1088 pixels were captured with a Basler acA2000-165uc camera mounted on a VZM 200i Zoom Imaging Lens (Edmund Optics, USA) during transcranial illumination with near-infrared laser light (CLD 1011LP, Thorlabs Inc, USA). Speckle data were analyzed with MATLAB software (RRID:SCR_001622; MathWorks, USA). In each experiment, values were calculated for a region of interest (100 × 100 pixels) in the second posterior bifurcation of the middle cerebral artery corresponding to the somatosensory barrel cortex contralateral to the whisker stimulation (Aronoff and Petersen, 2008).
Arterial blood gas measurements
Request a detailed protocolArterial blood samples were collected through a catheter implanted in the common carotid artery of anesthetized, intubated mice ventilated to normocapnia (FETCO2 = 3.5%) or experimental hypocapnia (FETCO2 = 2%) on a Minivent type 845 ventilator (Harvard Apparatus). The blood was immediately analyzed using an ABL80 Flex blood gas analyzer (Radiometer, Denmark).
Statistics of mouse experiments
Request a detailed protocolData are expressed as mean ± SEM and n equals number of mice (i.e. biological replicates). Probability (p) values less than 0.05 were considered statistically significant. Sample sizes were selected based on previous experience (Boedtkjer et al., 2006; Boedtkjer et al., 2016c; Boedtkjer et al., 2011; Thomsen et al., 2014) to allow detection of biologically relevant differences. If distributions were approximately Gaussian and variances equal between groups, we compared (a) one variable between two groups using unpaired two-tailed Student’s t-test, (b) one variable between three or more groups using one-way ANOVA followed by Dunnett’s post-tests, and (c) effects of two variables on a third variable using two-way ANOVA followed by Sidak’s post-tests. If the distributions showed unequal variance (i.e. p<0.05 by F test or by Brown-Forsythe and Bartlett’s tests) or significant difference from normality (i.e. p<0.05 by D'Agostino and Pearson or Shapiro-Wilk normality tests) due to right-skewness, we performed square root- or log-transformation. If variances were still unequal, we compared one variable between two groups based on unpaired two-tailed t-tests with Welch’s correction. If data distributions still did not pass normality tests, we used the non-parametric Mann-Whitney test to compare one variable between two groups. Concentration-response relationships were fitted to sigmoidal functions using least-square regression analyses, and the derived log(EC50), Hill Slope, and bottom values were compared using extra sum-of-squares F-tests. Investigators were not blinded for genotype during experiments. Data processing and statistical analyses were performed using Microsoft Office Excel 2016 (RRID:SCR_016137) and GraphPad Prism 7.05 (RRID:SCR_002798) software.
Burden analysis of loss-of-function variants
Request a detailed protocolWe studied association between rare variation in PTPRG and human ischemic vascular disease based on exome sequencing data covering ~50,000 participants from the UK Biobank (RRID:SCR_012815), which is a cohort with deep genetic and phenotypic data collected from the United Kingdom (Bycroft et al., 2018; Van Hout et al., 2019).
Population-level variants (UK Biobank data field 23170) generated with the Functionally Equivalent pipeline (Regier et al., 2018) were used after ensuring that the genomic region containing PTPRG was unaffected by the current known error in the exome analysis protocol (UK Biobank Resource 3802). We predicted the consequences of alternative alleles within the full-length transcript for PTPRG (ENST00000474889.6) using the Ensembl Variant Effect Predictor (VEP) (McLaren et al., 2016); and included in our analyses, the 334 missense variants localized outside exon-intron boundary regions and 75 variants with predicted loss-of-function. Carriers of the identified variants were extracted using the R-package qgg (Rohde et al., 2020). The identified loss-of-function variants included 36 with low (intronic splice region variants, synonymous splice region variants), 17 with moderate (in-frame insertions, in-frame deletions, missense splice region variants), and 22 with high (splice acceptor variants, splice donor variants, stop-gain variants, frameshift variants) predicted impact. In total, we identified 29,975 individuals who carried at least one missense variant and 195 individuals, who carried a variant with predicted loss-of-function. Of the 195 loss-of-function carriers, 123, 32, and 40 individuals carried variants with low, moderate, and high predicted impact, respectively.
Known disease diagnoses for UK Biobank participants (data field 41270) were extracted from hospital inpatient records coded according to the International Classification of Disease, version 10 (ICD-10; RRID:SCR_010349). Amongst the 49,953 participants with exome data, 2529 individuals were diagnosed with angina pectoris (I20), 972 with acute myocardial infarction (I21), 392 with cerebral ischemia (I63), and 153 with vascular intestinal disorders (K55.0). Self-reported health status was extracted for angina (data field 1074; reported by 1621 individuals), heart attack (data field 1075; reported by 1158 individuals), and stroke (data field 1081; reported by 686 individuals). We further extracted MRI-based values for 12,851 exome-sequenced individuals for whom left ventricular ejection fraction (data field 22420) was available from a fully automated analysis approach that has shown good correlation with manual analysis results performed by trained readers (Suinesiaputra et al., 2018). The cardiac MRI data set included 6482 carriers of missense PTPRG variants, 32 carriers of low impact loss-of function variants, 13 carriers of moderate impact loss-of-function variants, and 13 carriers of high impact loss-of-function variants.
We performed logistic regression analyses to calculate odds ratios (± SEM) for carrying missense or predicted loss-of-function PTPRG variants when diagnosed with ischemic vascular disease. We fist evaluated acute ischemic vascular disease as a whole for the brain, heart, and intestine (aggregate of ICD-10 diagnosis codes I20, I21, I63, and K55.0). For the group of individuals carrying a predicted loss-of-function variant of moderate or high impact, we then calculated separate odds ratios for I20, I21, and I63. The number of exome-sequenced participants in the UK Biobank diagnosed with K55.0 was too low for a meaningful separate statistical analysis. We next evaluated association of the PTPRG loss-of-function variants with self-reported angina, heart attack, and stroke. Finally, we evaluated the relation between carrier status for PTPRG variants and left ventricular ejection fraction by multiple linear regression analysis. The logistic as well as multiple linear regression analyses were corrected for sex (data field 31), age (data field 21022), body mass index (data field 21001), the first four genetic principal components (data field 22009), smoking status (data field 20116), dyslipidemia (ICD-10 diagnosis code E78.5), diabetes mellitus (ICD-10 diagnosis code E11), and hypertension (ICD-10 diagnosis code I10).
Data availability
All data generated or analysed during this study are presented in the manuscript.
References
-
Intracellular pH regulation in resting and contracting segments of rat mesenteric resistance vesselsThe Journal of Physiology 402:391–410.https://doi.org/10.1113/jphysiol.1988.sp017211
-
Chloride and bicarbonate transport in rat resistance arteriesThe Journal of Physiology 436:57–73.https://doi.org/10.1113/jphysiol.1991.sp018539
-
Adverse effects of artificial buffers on contractile responses of arterial and venous smooth muscleBritish Journal of Pharmacology 69:207–214.https://doi.org/10.1111/j.1476-5381.1980.tb07892.x
-
Layer, column and cell-type specific genetic manipulation in mouse barrel cortexFrontiers in Neuroscience 2:64–71.https://doi.org/10.3389/neuro.01.001.2008
-
Acidosis prevents and alkalosis augments endothelium-dependent contractions in mouse arteriesPflügers Archiv - European Journal of Physiology 466:295–305.https://doi.org/10.1007/s00424-013-1323-z
-
Cytosolic free calcium regulation in response to acute changes in intracellular pH in vascular smooth muscleAmerican Journal of Physiology-Cell Physiology 264:C932–C943.https://doi.org/10.1152/ajpcell.1993.264.4.C932
-
Antibody-independent localization of the electroneutral Na+-HCO3– cotransporter NBCn1 (slc4a7) in miceAmerican Journal of Physiology-Cell Physiology 294:C591–C603.https://doi.org/10.1152/ajpcell.00281.2007
-
Extracellular HCO3– is sensed by mouse cerebral arteries: regulation of tone by receptor protein tyrosine phosphatase γJournal of Cerebral Blood Flow & Metabolism 36:965–980.https://doi.org/10.1177/0271678X15610787
-
Acid-base regulation and sensing: accelerators and brakes in metabolic regulation of cerebrovascular toneJournal of Cerebral Blood Flow & Metabolism 38:588–602.https://doi.org/10.1177/0271678X17733868
-
Insulin inhibits Na+/H+ exchange in vascular smooth muscle and endothelial cells in situ: involvement of H2O2 and tyrosine phosphatase SHP-2American Journal of Physiology-Heart and Circulatory Physiology 296:H247–H255.https://doi.org/10.1152/ajpheart.00725.2008
-
Intracellular pH in the resistance vasculature: regulation and functional implicationsJournal of Vascular Research 49:479–496.https://doi.org/10.1159/000341235
-
Disturbed acid-base transport: an emerging cause of hypertensionFrontiers in Physiology 4:388.https://doi.org/10.3389/fphys.2013.00388
-
Characterisation and comparison of temporal release profiles of nitric oxide generating donorsJournal of Neuroscience Methods 245:116–124.https://doi.org/10.1016/j.jneumeth.2015.02.024
-
Protease-activated receptor-2 (PAR2) in cardiovascular systemVascular Pharmacology 43:247–253.https://doi.org/10.1016/j.vph.2005.07.009
-
Characterization of protease-activated receptor-2 immunoreactivity in normal human tissuesJournal of Histochemistry & Cytochemistry 46:157–164.https://doi.org/10.1177/002215549804600204
-
On the tonicity of the heart and blood vesselsThe Journal of Physiology 3:48–92.https://doi.org/10.1113/jphysiol.1880.sp000083
-
Brain extracellular pH: the main factor controlling cerebral blood flowScandinavian Journal of Clinical and Laboratory Investigation 22:247–251.https://doi.org/10.3109/00365516809167060
-
The ensembl variant effect predictorGenome Biology 17:122.https://doi.org/10.1186/s13059-016-0974-4
-
Carbonic anhydrase inhibitors modify intracellular pH transients and contractions of rat middle cerebral arteries during CO2/HCO3-– fluctuationsJournal of Cerebral Blood Flow & Metabolism 38:492–505.https://doi.org/10.1177/0271678X17699224
-
Qgg: an R package for large-scale quantitative genetic analysesBioinformatics 36:2614–2615.https://doi.org/10.1093/bioinformatics/btz955
-
Protons inhibit the BKCa channel of rat small artery smooth muscle cellsJournal of Vascular Research 38:30–38.https://doi.org/10.1159/000051027
-
Fully-automated left ventricular mass and volume MRI analysis in the UK Biobank population cohort: evaluation of initial resultsThe International Journal of Cardiovascular Imaging 34:281–291.https://doi.org/10.1007/s10554-017-1225-9
-
Intracellular acidification alters myogenic responsiveness and vasomotion of mouse middle cerebral arteriesJournal of Cerebral Blood Flow & Metabolism 34:161–168.https://doi.org/10.1038/jcbfm.2013.192
-
Enhanced nitric oxide signaling amplifies vasorelaxation of human Colon cancer feed arteriesAmerican Journal of Physiology-Heart and Circulatory Physiology 316:H245–H254.https://doi.org/10.1152/ajpheart.00368.2018
-
Micropuncture evaluation of the importance of perivascular pH for the arteriolar diameter on the brain surfacePflügers Archiv European Journal of Physiology 316:152–163.https://doi.org/10.1007/BF00586483
-
Acid-sensing ion channels: advances, questions and therapeutic opportunitiesTrends in Neurosciences 29:578–586.https://doi.org/10.1016/j.tins.2006.06.014
-
Role of receptor protein tyrosine phosphatase γ in sensing extracellular CO2 and HCO3–Journal of the American Society of Nephrology : JASN 27:2616–2621.https://doi.org/10.1681/ASN.2015040439
Article and author information
Author details
Funding
Det Frie Forskningsråd (4183-00258B)
- Ebbe Boedtkjer
Det Frie Forskningsråd (7025-00050A)
- Ebbe Boedtkjer
Lundbeckfonden (R93-A8859)
- Ebbe Boedtkjer
Lundbeckfonden (R287-2018-735)
- Palle D Rohde
MEMBRANES Research Center
- Ebbe Boedtkjer
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
The authors wish to thank Drs. Joseph Schlessinger (Yale University) and Walter F Boron (Case Western Reserve University) for generously providing the RPTPγ KO mice. We are grateful to Dr. Joanna Kalucka (Aarhus University) for assistance with the single-cell RNA sequencing analyses. We thank Viola Larsen, Liisa Jacobsen, Jørgen Andresen, and Jane Rønn (Aarhus University) for expert technical assistance. Access to the UK Biobank resource was approved under application number 60032.
Ethics
Human subjects: The research based on the UK Biobank resource was conducted under Application Number 60032.
Animal experimentation: This study was performed in accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. Protocols were approved by the Danish Animal Experiments Inspectorate (2016-15-0201-00982). All surgery was performed under general anesthesia, and every effort was made to minimize suffering.
Copyright
© 2020, Hansen et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,199
- views
-
- 158
- downloads
-
- 17
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Genetics and Genomics
- Immunology and Inflammation
Systemic lupus erythematosus (SLE) is an autoimmune disease, the pathophysiology and genetic basis of which are incompletely understood. Using a forward genetic screen in multiplex families with SLE, we identified an association between SLE and compound heterozygous deleterious variants in the non-receptor tyrosine kinases (NRTKs) ACK1 and BRK. Experimental blockade of ACK1 or BRK increased circulating autoantibodies in vivo in mice and exacerbated glomerular IgG deposits in an SLE mouse model. Mechanistically, NRTKs regulate activation, migration, and proliferation of immune cells. We found that the patients’ ACK1 and BRK variants impair efferocytosis, the MERTK-mediated anti-inflammatory response to apoptotic cells, in human induced pluripotent stem cell (hiPSC)-derived macrophages, which may contribute to SLE pathogenesis. Overall, our data suggest that ACK1 and BRK deficiencies are associated with human SLE and impair efferocytosis in macrophages.
-
- Genetics and Genomics
- Microbiology and Infectious Disease
The sustained success of Mycobacterium tuberculosis as a pathogen arises from its ability to persist within macrophages for extended periods and its limited responsiveness to antibiotics. Furthermore, the high incidence of resistance to the few available antituberculosis drugs is a significant concern, especially since the driving forces of the emergence of drug resistance are not clear. Drug-resistant strains of Mycobacterium tuberculosis can emerge through de novo mutations, however, mycobacterial mutation rates are low. To unravel the effects of antibiotic pressure on genome stability, we determined the genetic variability, phenotypic tolerance, DNA repair system activation, and dNTP pool upon treatment with current antibiotics using Mycobacterium smegmatis. Whole-genome sequencing revealed no significant increase in mutation rates after prolonged exposure to first-line antibiotics. However, the phenotypic fluctuation assay indicated rapid adaptation to antibiotics mediated by non-genetic factors. The upregulation of DNA repair genes, measured using qPCR, suggests that genomic integrity may be maintained through the activation of specific DNA repair pathways. Our results, indicating that antibiotic exposure does not result in de novo adaptive mutagenesis under laboratory conditions, do not lend support to the model suggesting antibiotic resistance development through drug pressure-induced microevolution.