Orphan receptor GPR37L1 contributes to the sexual dimorphism of central cardiovascular control
Biology of Sex Differences volume 9, Article number: 14 (2018)
Over 100 mammalian G protein-coupled receptors are yet to be matched with endogenous ligands; these so-called orphans are prospective drug targets for the treatment of disease. GPR37L1 is one such orphan, abundant in the brain and detectable as mRNA in the heart and kidney. GPR37L1 ablation was reported to cause hypertension and left ventricular hypertrophy, and thus, we sought to further define the role of GPR37L1 in blood pressure homeostasis.
We investigated the cardiovascular effects of GPR37L1 using wild-type (GPR37L1wt/wt) and null (GPR37L1KO/KO) mice established on a C57BL/6J background, both under baseline conditions and during AngII infusion. We profiled GPR37L1 tissue expression, examining the endogenous receptor by immunoblotting and a β-galactosidase reporter mouse by immunohistochemistry.
GPR37L1 protein was abundant in the brain but not detectable in the heart and kidney. We measured blood pressure in GPR37L1wt/wt and GPR37L1KO/KO mice and found that deletion of GPR37L1 causes a female-specific increase in systolic, diastolic, and mean arterial pressures. When challenged with short-term AngII infusion, only male GPR37L1KO/KO mice developed exacerbated left ventricular hypertrophy and evidence of heart failure, while the female GPR37L1KO/KO mice were protected from cardiac fibrosis.
Despite its absence in the heart and kidney, GPR37L1 regulates baseline blood pressure in female mice and is crucial for cardiovascular compensatory responses in males. The expression of GPR37L1 in the brain, yet absence from peripheral cardiovascular tissues, suggests this orphan receptor is a hitherto unknown contributor to central cardiovascular control.
High blood pressure (BP) is responsible for > 45% of heart disease deaths worldwide and is the single biggest contributing risk factor to global disease burden . The human cost of hypertension is partially explained by the disorder’s multifactorial pathogenesis and is accordingly difficult to treat . Approximately 30% of patients receiving treatment for hypertension do not have their BP adequately controlled by current pharmacotherapy . Clearly, a better understanding of the mechanisms governing BP is necessary to more effectively manage cardiovascular disease.
One such BP regulator may be the G protein-coupled receptor (GPCR) GPR37L1. This orphan receptor (orphan—without a known endogenous ligand) is predominantly expressed in the brain [4, 5] and contributes to cerebellar development and motor skills . Gpr37l1 also lies within eight quantitative trait loci linked to BP regulation in the rat . Moreover, a previous study  reported that GPR37L1 knockout mice have left ventricular hypertrophy (LVH) and ~ 60 mmHg higher systolic BP than that of cardiac-specific, GPR37L1-overexpressing transgenic mice. However, the mouse strain, sex, age, heart rate, and depth of anesthesia used for BP measurement were not indicated, and BP was not compared to appropriate wild-type controls. Here, we sought to definitively establish whether GPR37L1 regulates BP using GPR37L1 null mice that we generated and maintained on a pure C57BL/6J background [9, 10]. We report that GPR37L1 is indeed involved in BP regulation but in a sexually dimorphic manner.
Mouse line generation and maintenance
Male and female C57BL/6J GPR37L1wt/wt, GPR37L1flx/flx, GPR37L1wt/KO, GPR37L1KO/KO, and GPR37L1lacZ/wt were generated using the European Conditional Mouse Mutagenesis Program (EUCOMM) method for conditional gene inactivation using Cre/loxP and Flp1/FRT gene trapping  described as GOI#1 in a previous report . Mice were group-housed in a temperature-controlled environment with a 12-h light-dark cycle (lights on at 0700). Standard chow and water were available ad libitum. All animal work was performed in accordance with the Australia Code for the Care and Use of Animals for Scientific Purposes, 8th Edition (2013), and was approved by the relevant Animal Ethics Committees (Garvan Institute for Medical Research/St Vincent’s Hospital, project number AEC 13/30; University of Melbourne, project number 1212573). All animals were entered into the study in a randomized order.
Tissue was homogenized in TRIzol (15,596,026, Thermo Scientific, USA) using a PRO200 homogenizer (Bio-Gen, USA), and RNA was extracted using manufacturer’s protocol (Thermo Scientific). DNA contamination was removed from extracted RNA with TURBO DNAse kit treatment (AM1907, Invitrogen, USA) before cDNA synthesis was performed using SuperScript® III First-Strand Synthesis Supermix as per manufacturer’s protocol (11752050, Invitrogen). cDNA was stored at − 20 °C until use. Quantitative PCR (qPCR) was performed on cDNA using TaqMan Gene Expression Master Mix (4369016, Applied Biosystems, USA) and TaqMan probes (Applied Biosystems). All qPCR was performed using a LightCycler 480 (Roche, Switzerland) (cycling conditions: preincubation at 95 °C for 10 min; 55 amplification cycles of 95 °C for 10 s, 58 °C for 30 s, and 72 °C for 30 s; final cooling at 50 °C for 10 s). Probe Mm00661872_m1 was used to detect GPR37L1 in mouse. The following probes were used to detect components of the fetal gene program: ANP (Mm01255747_g1), BNP (Mm01255770_g1), β-MHC (Mm00600555_m1), and α-SkA (Mm00808218_g1). ΔCt values were calculated in reference to HPRT (mouse, Mm03024075_m1) or when GPR37L1 transcript abundance was determined in EUCOMM mice; ΔCt values were calculated in reference to β2M (mouse, Mm00437762_m1). Relative mRNA abundance was determined using the 2−ΔΔCt method .
Tissues for pilot assay of GPR37L1 transcription
In reference to the experiment shown in Additional file 1: Figure S1, surplus tissue (left and right ventricles, left and right atria, aorta, kidney, lung, liver, and brain) from sham-operated wild-type C57BL/6J mice (12 weeks old, male), as harvested for a previous study , was used to generate cDNA to quantify GPR37L1 by qPCR.
Tissue fixation for β-galactosidase imaging
Tissue fixation was performed as adapted from a previous method . At ≥ 12 weeks of age, GPR37L1wt/wt and GPR37L1lacZ/wt mice of either sex were anesthetized with 5% isofluorane in oxygen (flow rate 1 L/min) delivered via nose cone. The thoracic cavity was opened, and the LV was punctured with a needle that was then used to perfuse the heart with phosphate-buffered saline (PBS) (1 min; flow rate 20 mL/min), achieving sacrifice by exsanguination. Immediately after, tissues were fixed by further 2 min perfusion with 4% paraformaldehyde (PFA) in PBS (flow rate 20 mL/min). Following fixation, the brain, heart, and kidney were collected and further fixed in 4% PFA/PBS for 48 h at room temperature (RT). After post-fixation, tissue was stored in sucrose cryoprotectant buffer (5 mmol/L Trizma base, 77 mmol/L NaCl, 4 mmol/L Na2HPO4, 1.3 mmol/L NaH2PO4, 20% w/v sucrose) at 4 °C for ≥ 48 h. Fixed and cryoprotected tissue was embedded in Tissue Freezing Medium (TFM-5, Triangle Biomedical, USA) and frozen. Frozen tissue was cut into 12 μm sections using a CM1950 Cryostat (Leica Biosystems, Germany) and mounted on Superfrost Plus Slides (Thermofisher Scientific).
β-galactosidase immunofluorescent staining
The protocol was performed at RT unless otherwise specified. Slides containing mounted tissue sections from the above procedure were washed for 10 min, three times, in immunobuffer [TPBS (10 mmol/L Tris, 8 mmol/L Na2HPO4, 2.6 mmol/L NaH2PO4, 0.9% NaCl, 0.05% thimerosal) with 0.3% Triton X-100]. Non-specific antibody binding was blocked by incubating each section in immunobuffer containing 10% v/v normal donkey serum (NDS) (017-000-121, Jackson ImmunoResearch Laboratories, USA) for 30 min, in a humidifying chamber. Primary antibodies [1:250 chicken anti-β-galactosidase (AB9361, Abcam, Cambridge, UK) for all tissue, 1:500 rabbit anti-GFAP (Z0334, Agilent Technologies, USA) for brain only] were diluted in immunobuffer with 10% v/v NDS and applied to each section for 16 h, in a humidifying chamber at 4 °C. Slides were then washed for 10 min, three times, in TPBS. For the remainder of the protocol, slides were shielded from ambient light. Secondary antibodies [1:500 donkey anti-chicken Cy3 (703-165-155, Jackson ImmunoResearch Laboratories) for all tissue, 1:500 donkey anti-rabbit Alexa Fluor 647 (711-605-152, Jackson ImmunoResearch) for brain only], DAPI (1:5000 4′,6-diamidino-2-phenylindole dihydrochloride; D9542, Sigma Aldrich, USA) nuclear stain, and wheat germ agglutinin (WGA) membrane stain (for heart and kidney only, 1:500 Alexa Fluor 488 WGA conjugate; W11261, Invitrogen) were diluted in immunobuffer containing 10% v/v NDS and applied to each section for 1 h, in a humidifying chamber. Slides were then washed for 10 min in TPBS. To reduce autofluorescence, each section was then treated with 0.1% w/v Sudan Black B (S2380, Sigma Aldrich) in 70% v/v ethanol for 1 h, in a humidifying chamber. Excess dye was then removed by washing for 10 min, three times, in TPBS. Slides were then dried for 16 h in a humidifying chamber before being sealed with glass coverslips mounted using 2.5% w/v DABCO (1,4-diazabicyclo-[2,2,2]-octane; D2522, Sigma Aldrich), 10% w/v polyvinyl alcohol (P8136, Sigma Aldrich), 25% w/v glycerol (G5516, Sigma Aldrich), and 100 mmol/L Tris buffer (pH 8.7), as adapted from a previous method . Slides were imaged on a LSM 7 Duo confocal microscope (Zeiss, Germany).
Immunoblotting and densitometry
Tissue homogenates and crude membranes were prepared as previously described [15, 16]. SDS-PAGE and immunoblotting was performed as previously described . GPR37L1 was detected using goat anti-GPR37L1 (C-12) antibody (1:1000; sc-164532, Santa Cruz Biotechnology, USA) that was previously confirmed to be GPR37L1-specific using knockout tissue  and then rabbit anti-goat antibody (1:7500; 61-1620, Invitrogen). β-galactosidase was detected using rabbit anti-β-galactosidase antibody (1:5000; A-11132; Invitrogen) and then donkey anti-rabbit antibody (1:15000, NA934, GE Healthcare, Australia). GAPDH was detected using rabbit anti-GAPDH antibody (1:10000; 14C10, Cell Signaling, USA) and then donkey anti-rabbit antibody (1:20000, NA934, GE Healthcare). SuperSignal™ West Pico (34080, Invitrogen) or Clarity™ (1705060, Bio-Rad, USA) were used as chemiluminescent substrates. Densitometry was performed using ImageJ software (https://imagej.nih.gov/ij/) on non-saturated chemiluminescent exposures. Pixel density of background was measured and subtracted from subsequent measurements taken on the same image. Pixel density of GPR37L1 was measured in reference to the predominant, cleaved receptor species  and adjusted to equivalent protein abundance using each sample’s own GAPDH pixel density score.
Hemodynamic measurement and analysis
At age 10–12 weeks, mice were anesthetized with 1–2% isofluorane (Zoetis, USA) in oxygen (flow rate 0.5 L/min), delivered via nose cone regulated using a mechanical ventilator (150 strokes/min, 200 μl stroke volume). Hemodynamics in the aorta and LV were recorded with a high-fidelity pressure transducer (Millar Instruments, USA) following its insertion into the right carotid artery as described previously [17, 18]. Heart rate was maintained at approximately 500 beats per minute (bpm) by manual adjustment of the isofluorane concentration. Data were recorded and analyzed using AcqKnowledge software (version 3.9.0 for Windows, BIOPAC Systems, USA) by an observer blinded to genotype (and to treatment where applicable). Hemodynamics on osmotic mini pump animals were recorded using a Transonic Scisense (FTH-1211B-001, USA) high-fidelity pressure transducer.
Radiotelemetry blood pressure measurement
At 12–14 weeks of age, mice were anesthetized with 2% isofluorane (Abbott Laboratories, USA), with the rate adjusted to provide a deep surgical plane. PA-C10 telemeters (Data Sciences International, USA) were implanted as previously described , with the pressure-sensing catheter inserted into the left carotid artery. Analgesia was administered before and after surgery (carprofen, 0.5 mg/100 g, IP; Norbrook, UK). Hydration was maintained by administration of 0.5 mL 0.9% NaCl (IP), given at the time of operation and again 24 h later. Approximately 2 weeks after telemeter implantation, blood pressure monitoring was performed twice daily in the home cage, between 09:00 and 12:00 (light phase) and between 21:00 and 24:00 (dark phase) for 9 days. The average of the 9 days of measurements in each phase was taken to give a single value per animal per phase. [N.B. radiotelemetry blood pressure monitoring was only performed during initial phenotyping experiments in Fig. 2c, d, as indicated].
Angiotensin II infusion
Osmotic mini pumps (product 1007D, Alzet) were filled with either vehicle (0.15 mol/L NaCl, 1 mmol/L acetic acid) or AngII (A9525, Sigma Aldrich) at a concentration adjusted to allow 2 mg/kg/day release of AngII for 7 days. After filling, pumps were primed by incubation in 0.9% NaCl for 16 h at 37 °C. Mice (10–12 weeks old) were anesthetized [inhaled 3–4% isofluorane (Zoetis) in oxygen (flow rate 0.8–1.0 L/min) delivered by nose cone] and pumps implanted subcutaneously, as per the manufacturer's instructions. Following pump implantation, topical analgesia was administered in the form of bupivacaine (AstraZeneca, UK), and systemic analgesia was administered (SC) in the form of buprenorphine (0.075 mg/kg; Reckitt Benckiser, UK), before completion of surgery and recovery from anesthesia. Animals were allowed to recover from anesthesia on a heated pad, before being single-housed for the remainder of the study. [N.B. blood pressure following the 7 days of infusion was measured by anesthetized catheterization, as indicated].
Tissues were retrieved after sacrifice, briefly rinsed of blood in 0.9% NaCl, and then blotted dry. Tissue weights were measured immediately and expressed as a ratio of each animal’s own tibia length, measured with caliMAX fine calipers (WIHA, Germany).
Cardiomyocyte density was determined using an adaptation from a previous method . Briefly, LV tissue was fixed in 4% PFA/PBS for 24 h at RT before post-fixation, sectioning, and processing as described above for immunofluorescent staining and confocal microscopy. Cell membranes were visualized using 1:500 Alexa Fluor 488 WGA conjugate. Using ImageJ, an observer blinded to sex, genotype, and treatment counted cardiomyocytes from two representative areas of each LV and averaged them to give one value per animal. Cardiomyocyte counts were expressed as cells per square millimeter. Cardiomyocyte density was used instead of cardiomyocyte cross-sectional area as we found density to show less inter-observer variability (data not shown).
LV tissue was fixed as described for cardiomyocyte density. Cytoplasm (red) and collagen (blue) were stained using Masson’s trichrome stain as per the manufacturer’s protocol (HT15, Sigma Aldrich). Coverslips were mounted onto slides using DEPEX mounting media (VWR International, USA). Sections were visualized with bright-field microscopy on a DM6000 power mosaic microscope (Leica Biosystems), and images were retained in RGB format. The degree of fibrosis was quantified using ImageJ by an observer blinded to sex, genotype, and treatment and expressed as a percentage of area occupied by blue-hued pixels relative to total cross-sectional tissue area, as adapted from previous methods .
RNA sequencing analysis
Human postmortem brain RNA sequencing data was obtained from the Genotype-Tissue Expression (GTEx) Consortium . GTEx raw sequence counts were examined for differentially expressed genes using the edgeR statistical package , where samples were grouped by genders for each brain tissue sample.
All data are presented as mean ± s.e.m. The number of independent observations (n) in animal/human data represents a single animal/donor. Where multiple measurements were taken per animal, measurements were averaged to give a single observation (n) per animal. Where two untreated genotypes (per sex) were analyzed, a two-tailed Student’s t test was used, or a one-tailed Mann-Whitney test, as indicated. Where > 2 untreated genotypes (per sex) were analyzed, a one-way ANOVA was used with Dunnett’s multiple comparisons tests, or Kruskal-Wallis test with Dunn’s multiple comparisons test where ANOVA assumptions were violated. Where two genotypes (per sex), with and without treatment, were analyzed, a two-way ANOVA was used with Tukey’s multiple comparisons test. For all tests, p ≤ 0.05 was considered significant.
GPR37L1 protein is abundant in the central nervous system but not detectable in the heart or kidney
Given the reported role of GPR37L1 in the cardiovascular system , we investigated GPR37L1 transcription by qPCR in several mouse tissues of interest, including the heart and kidney. We noted that GPR37L1 mRNA appeared at late cycle numbers in both tissues, suggesting very low relative expression (Additional file 1: Figure S1). To clarify the expression profile of GPR37L1, we used EUCOMM ES cells  to generate GPR37L1lacZ/wt mice (expressing one allele of β-galactosidase in place of GPR37L1) as we have described . Consistent with previous reports of GPR37L1 localization [4,5,6, 23], β-galactosidase expression was most abundant in the brain, particularly in the cerebellar Bergmann glia (Fig. 1a, Additional file 1: Figure S2A). Astrocyte β-galactosidase co-localization was observed throughout the brain but was most evident in the hippocampus (Additional file 1: Figure S2B), while other regions such as the brain stem (Fig. 1a) and neocortex (Additional file 1: Figure S2C) contained GFAP-negative β-galactosidase-positive cells that were consistent with the reported expression of GPR37L1 in cells of oligodendrocyte lineage [23, 24]. In several key cardiovascular centers of the CNS including the caudal and rostral ventrolateral medulla (CVLM, RVLM), the nucleus of the solitary tract (NTS), and the A5 nucleus, diffuse β-galactosidase staining was observed close to, but not within, tyrosine hydroxylase-positive catecholaminergic neurons (Additional file 1: Figure S3). In contrast, we were unable to detect GPR37L1 β-galactosidase reporter immunofluorescence in peripheral tissues. Although we measured low levels of GPR37L1 mRNA in the mouse heart (Additional file 1: Figure S1), we could not detect β-galactosidase by immunofluorescence (Fig. 1b, Additional file 1: Figure S4A) or Western blot (Additional file 1: Figure S4B), nor wild-type GPR37L1 using a validated  anti-GPR37L1 antibody (Additional file 1: Figure S4C-E). Similarly, we did not detect β-galactosidase in the renal medulla (Fig. 1c) or cortex (Additional file 1: Figure S5A), nor could it be detected in GPR37L1lacZ/wt kidney homogenates (Additional file 1: Figure S5B). Endogenous GPR37L1 was not detected in GPR37L1wt/wt kidney (Additional file 1: Figure S5C). Because only GPR37L1 mRNA, but not protein, was found in the heart and kidney, we hypothesize that GPR37L1 exerts its cardiovascular effect via the CNS.
GPR37L1 regulates blood pressure in female mice
To investigate the role of GPR37L1 in the cardiovascular system, we further generated the GPR37L1flx/flx and GPR37L1KO/KO mice [9, 10]. Receptor ablation in both genders was verified in the cerebellum (mRNA and protein, respectively; Fig. 1d–f), and in the LV (mRNA; Additional file 1: Figure S4F). GPR37L1 expression was unchanged in GPR37L1flx/flx as compared to GPR37L1wt/wt mice (generated from GPR37L1wt/KO breeding pairs) (Fig. 1d, Additional file 1: Figure S6A and B). Thus, GPR37L1wt/wt mice were used as controls for the remainder of our studies.
Baseline cardiovascular measures were first assessed in anesthetized females by micromanometry (Fig. 2a) at a constant heart rate of ~ 500 bpm (Fig. 2b, Additional file 1: Table S1). Female GPR37L1wt/KO and GPR37L1KO/KO mice displayed significantly elevated aortic systolic pressure compared to controls (Fig. 2a, Additional file 1: Table S1). Furthermore, GPR37L1KO/KO females showed increased aortic diastolic, mean arterial, and pulse pressures (Fig. 2a, Additional file 1: Table S1). No genotype differences in LV systolic pressure, end-diastolic pressure, and dP/dT max or dP/dT min were observed. To confirm that this phenotype was preserved in conscious animals, radiotelemetry BP monitoring was performed on a separate cohort of conscious GPR37L1wt/wt and GPR37L1KO/KO females. Again, GPR37L1KO/KO females had elevated systolic and diastolic pressure, during both light and dark phases of the diurnal cycle, relative to GPR37L1wt/wt controls (Fig. 2c).
To investigate if GPR37L1 deletion is deleterious in the setting of hypertensive stress, female GPR37L1wt/wt and GPR37L1KO/KO mice were infused for 7 days with AngII (2 mg/kg/day) or vehicle, before anesthetized hemodynamic measurement by pressure-transducing catheter and tissue collection. AngII infusion caused greater increases in systolic and mean arterial pressure in the GPR37L1wt/wt females (Table 1), presumably because GPR37L1KO/KO females have higher baseline systolic pressure (one-tailed Mann-Whitney test, p ≤ 0.05), consistent with our previous, independent cohort (Fig. 2). However, BP after 7 days of AngII did not differ between genotypes (Table 1), and both developed equivalent gross, cellular, and molecular LVH (Fig. 3a, b, Additional file 1: Figure S7, Figure S8), without progression to overt heart failure (see relative lung and right ventricle weights; Fig. 3c, Table 1). Although GPR37L1wt/wt females developed cardiac fibrosis in response to AngII treatment, this was not observed in GPR37L1KO/KO females (Fig. 4).
GPR37L1 protects against AngII cardiovascular stress in male mice
Baseline cardiovascular measures were also assessed in males. GPR37L1wt/wt, GPR37L1wt/KO and GPR37L1KO/KO males showed no differences in any hemodynamic parameter by micromanometry (Fig. 2a, Additional file 1: Table S1). BP radiotelemetry was performed on a separate cohort of conscious male GPR37L1wt/wt and GPR37L1KO/KO animals. Consistent with previous reports , males had higher systolic and diastolic BP than females (Fig. 2c) and lower heart rates (Fig. 2d). The expected circadian dip was seen in both sexes (Fig. 2c).
In contrast to the females, hypertension was not sustained in male mice after 7 days of AngII infusion (Table 2), as measured by pressure-transducing catheter under anesthesia, although it was clear from the various physiological endpoints that the mice received pressure overload. While AngII treatment induced LVH for males of both genotypes, the extent of this enlargement was significantly greater in GPR37L1KO/KO mice (Fig. 3a, Table 2). Cardiomyocyte density reflected this difference, with only GPR37L1KO/KO LV sections indicating cardiomyocyte hypertrophy in response to AngII (Fig. 3b, Additional file 1: Figure S7). In concert with exacerbated LVH, male AngII-treated GPR37L1KO/KO mice had increased relative lung weight (Fig. 3c), indicative of pulmonary edema resulting from congestive heart failure [26,27,28] and suggesting that GPR37L1KO/KO males are susceptible to AngII-induced heart failure. Indeed, we observed a blunted contractile response in AngII-treated GPR37L1KO/KO males compared to GPR37L1wt/wt counterparts (Table 2). Consistent with our previous cohort (Fig. 2), the baseline aortic systolic pressure of male GPR37L1KO/KO did not differ from that of GPR37L1wt/wt controls (one-tailed Mann-Whitney test, p > 0.05) (Table 2).
As for the females, GPR37L1KO/KO males appeared to be protected from myocardial fibrosis (Fig. 4a, c). However, the more marked hypertrophy and therefore increased tissue area in the GPR37L1KO/KO males may have obscured the extent of absolute fibrosis which was not grossly different between genotypes (Fig. 4a).
Here, we evaluated the cardiovascular effects of GPR37L1 using a gene inactivation mouse model. We found that BP was increased only in female GPR37L1KO/KO mice, compared to GPR37L1wt/wt counterparts. In contrast, AngII-induced LVH was exacerbated only in male GPR37L1KO/KO mice and was associated with depressed contractility (reduced dP/dT max) and heart failure (increased relative lung weight). Using a range of approaches, we have observed that the orphan receptor GPR37L1 is predominantly expressed in glial cells of the brain, yet we could not detect the protein in peripheral tissues. Taken together, these findings indicate GPR37L1 is an important central regulator of BP homeostasis but that its cardiovascular effects differ in males and females.
We show that female GPR37L1KO/KO mice had ~ 12 mmHg greater systolic BP than GPR37L1wt/wt counterparts. Such an increase in BP is not trivial, as it would translate to a doubling in cardiovascular disease risk if it were observed in humans . Basal BP was unaltered in GPR37L1KO/KO males, but they failed to mount compensatory responses to AngII-induced hypertensive stress, exhibiting heightened LVH, impaired LV contractile function, and lung congestion. The most parsimonious interpretation of our findings is that GPR37L1KO/KO males had an exaggerated pressor response to AngII. However, we did not observe sustained AngII hypertension in males of either genotype, despite each having obvious LVH, and this is a potential limitation of the present study. Yet, given that LVH is known to increase in proportion to BP in both mouse and human hypertension [30, 31], and that males are more susceptible to AngII-induced hypertension than females , the data suggest that AngII induced pressure overload in males of both genotypes, causing systolic dysfunction, as seen in previous rodent pressure overload studies . Indeed, only males lost body weight in response to AngII, supportive of AngII being more detrimental to their health (Table 2). This broadly agrees with human heart failure epidemiology—low systolic pressure and LV function at hospital admission is more common in men . Further, systolic function is better preserved in women with aortic stenosis than men, despite equivalent hypertensive stress . Interestingly, GPR37L1KO/KO females appeared to be protected from cardiac fibrosis in response to AngII, despite their equivalent LVH and BP following infusion. The lack of detectable GPR37L1 protein in the heart suggests that, unlike the direct, pro-fibrotic effects of cardiac-expressed GPCRs (e.g., the β2 adrenoceptor, AT1 and ETA receptors ), GPR37L1 may indirectly influence fibrosis by central modulation of circulatory factors that regulate the extracellular matrix, though this remains to be investigated. Alternatively, it is possible that the reduced cardiac fibrosis is indirectly linked to LV GPR37L1 transcript, rather than to protein, through regulation of miRNA, which are known to play roles in cardiac fibrosis . Ultimately, detailed studies using radiotelemetry in concert with echocardiography are required to fully evaluate temporal changes in hemodynamic response to AngII, and the mechanism of decompensation in the GPR37L1KO/KO male mice, and whether GPR37L1 conveys cardioprotection in other cardiovascular disease models (i.e., thoracic aortic constriction, myocardial infarction). Such studies are currently underway.
Our observation of sexual dimorphism in the cardiovascular effects of GPR37L1 is intriguing, although sex differences in the cardiovascular control systems are known to contribute to differential disease etiology and outcomes [38, 39]. While GPR37L1 had previously been linked to hypertension , the animals’ sex was not reported. Meanwhile, previous studies of GPR37L1 in the brain investigated only male mice [6, 40] or did not disclose sex . Notably, sexual dimorphism has been reported for GPR37, the closest phylogenetic neighbor of GPR37L1 , and shown to be more abundant in human females in a variety of brain sub-regions at different developmental stages . Further, GPR37 deletion in mice leads to sex-dependent effects on neurotransmitter levels and anxiety . In our study, we did not see significant gender differences in cerebellar GPR37L1 mRNA or protein levels (Additional file 1: Figure S9A–C). Additionally, our analysis of human postmortem brain RNA-sequencing (collected by the Genotype-Tissue Expression Consortium ) shows the cortex to be the only sub-region with sex differences in GPR37L1 mRNA abundance, which is increased in females (Additional file 1: Figure S9D, Table S4). Taken together, this suggests that GPR37L1 may exert its gender-specific effects on BP via a mechanism that is independent of its expression level. Such effects might involve hormonal influences; estrogen, for example, is well known to sensitize central baroreflex and modulate sympathetic nerve activity .
Sex differences alone cannot reconcile the BP phenotype of GPR37L1 deletion in the present study (+ 12 mmHg, female-specific) with that previously reported (+ 60 mmHg, sex not disclosed) . Although both studies determined BP in anesthetized animals, our additional baseline BP measures by radiotelemetry in conscious, free-moving mice show that the ~ 12 mmHg increase in GPR37L1KO/KO females is independent of heart rate and anesthesia response. The magnitude of the previously reported GPR37L1 hypertensive phenotype  is also difficult to reconcile with the effects of deletion of other major cardiovascular regulators. For example, deletion of regulator of G protein signaling 2 , endothelial nitric oxide synthase , adenosine A2A receptor , and bradykinin B2 receptor  reportedly increased BP by 45, 28, 25, and 20 mmHg, respectively. Therefore, it is likely that the effect of GPR37L1 on BP was overestimated due to the lack of comparisons to wild-type controls. Moreover, we obtained the same reported GPR37L1 knockout mouse line  (reanimated from frozen sperm by Charles River Laboratories, Japan) and failed to detect a change in BP by tail cuff plethysmography in a small cohort, suggesting that any BP difference in these mice is also likely to be modest (data not shown). We additionally found these mice to be of mixed genetic background. Such confounding genetic issues were avoided in our study by the use of EUCOMM conditional-ready mice [9, 10], enabling comparison of wild-type and null mice on the same C57BL/6J background.
How then might GPR37L1 regulate cardiovascular function? While our laboratory has shown GPR37L1 to constitutively couple to the Gαs G protein , and to share antagonists with orexin and bombesin GPCRs , the lack of a high-affinity ligand means pharmacological probing of its physiology is intractable at present. GPR37L1 expression in the brain is well established, while its role and presence in the periphery remains unclear . For example, absence of GPR37L1 mRNA in human heart and kidney has been reported , while another study reported low levels of GPR37L1 mRNA in both tissues in the rat . Here, we also detected GPR37L1 mRNA at very low levels in the heart and kidney but could not detect GPR37L1 protein in either. Our data agree with that of the Human Protein Atlas , which reports no expression of GPR37L1 protein outside the brain. While vascular function remains to be explored in GPR37L1KO/KO mice, we did not observe any vascular β-galactosidase immunoreactivity in the heart and kidney of GPR37L1lacZ/wt mice, and the RNA abundance of GPR37L1 in the human artery and aorta is comparable to that of the heart and kidney , where we could not detect receptor protein. This raises the fascinating possibility that GPR37L1 is a previously undiscovered component of the CNS-cardiovascular axis. Widely expressed in astrocytes throughout the brain/brainstem, GPR37L1 may exert cardiovascular control through its proximity to autonomic control regions such as the CVLM, RVLM, and NTS , where astrocytes contribute to the brain renin-angiotensin system [53, 54]. Alternatively, GPR37L1 could be involved in cerebellar circuits known to influence the cardiovascular system [55, 56]. Indeed, the cerebellum has a long-suspected role in the modulation of baroreflexes and the development of essential hypertension [57, 58]. More speculatively, cerebellar control of respiration  may exert secondary influences on BP, with a recent study implicating excessive respiratory modulation by brainstem neurons in the development of rat and human hypertension .
The present study found that deletion of the orphan GPCR, GPR37L1, causes increased systolic and diastolic blood pressure in female mice. In males, however, deletion of GPR37L1 blunted a compensatory response to AngII challenge. Although we observed low levels of GPR37L1 mRNA in the mouse heart and kidney, we could not detect GPR37L1 protein in either organ. Thus, while the mechanism by which GPR37L1 regulates BP in a sexually dimorphic manner has yet to be determined, it is likely its role in cardiovascular homeostasis is confined to its CNS actions, where the protein is abundant. Thus, the present study suggests that GPR37L1 is a contributor to CNS cardiovascular control, and the sex differences therein, and may contribute to the role of the CNS in hypertension and heart failure.
Atrial natriuretic peptide
Brain natriuretic peptide
Central nervous system
Caudal ventrolateral medulla
European Conditional Mouse Mutagenesis Program
Glyceraldehyde 3-phosphate dehydrogenase
Glial fibrillary acidic protein
G protein-coupled receptor
Genotype-Tissue Expression Consortium
Hypoxanthine guanine phosphoribosyl transferase
Left ventricular hypertrophy
Normal donkey serum
Nucleus of the solitary tract
Quantitative polymerase chain reaction
Rostral ventrolateral medulla
Wheat germ agglutinin
β-myosin heavy chain
WHO. A global brief on hypertension: silent killer, global public health crisis. 2013.
Coffman TM. Under pressure: the search for the essential mechanisms of hypertension. Nat Med. 2011;17(11):1402–9.
Rimoldi SF, Messerli FH, Bangalore S, Scherrer U. Resistant hypertension: what the cardiologist needs to know. Eur Heart J. 2015;36(40):2686–95.
Valdenaire O, Giller T, Breu V, Ardati A, Schweizer A, Richards JG. A new family of orphan G protein-coupled receptors predominantly expressed in the brain. FEBS Lett. 1998;424(3):193–6.
Zeng Z, Su K, Kyaw H, Li Y. A novel endothelin receptor type-B-like gene enriched in the brain. Biochem Biophys Res Commun. 1997;233(2):559–67.
Marazziti D, Di Pietro C, Golini E, Mandillo S, La Sala G, Matteoni R, et al. Precocious cerebellum development and improved motor functions in mice lacking the astrocyte cilium-, patched 1-associated Gpr37l1 receptor. Proc Natl Acad Sci U S A. 2013;110(41):16486–91.
Kent WJ, Sugnet CW, Furey TS, Roskin KM, Pringle TH, Zahler AM, et al. The human genome browser at UCSC. Genome Res. 2002;12(6):996–1006.
Min KD, Asakura M, Liao Y, Nakamaru K, Okazaki H, Takahashi T, et al. Identification of genes related to heart failure using global gene expression profiling of human failing myocardium. Biochem Biophys Res Commun. 2010;393(1):55–60.
Skarnes WC, Rosen B, West AP, Koutsourakis M, Bushell W, Iyer V, et al. A conditional knockout resource for the genome-wide study of mouse gene function. Nature. 2011;474(7351):337–42.
Coleman JL, Brennan K, Ngo T, Balaji P, Graham RM, Smith NJ. Rapid knockout and reporter mouse line generation and breeding colony establishment using EUCOMM conditional-ready embryonic stem cells: a case study. Front Endocrinol. 2015;6:105.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods. 2001;25(4):402–8.
Soetanto R, Hynes CJ, Patel HR, Humphreys DT, Evers M, Duan G, et al. Role of miRNAs and alternative mRNA 3′-end cleavage and polyadenylation of their mRNA targets in cardiomyocyte hypertrophy. Biochim Biophys Acta. 2016;1859(5):744–56.
Chen D, Bassi JK, Walther T, Thomas WG, Allen AM. Expression of angiotensin type 1A receptors in C1 neurons restores the sympathoexcitation to angiotensin in the rostral ventrolateral medulla of angiotensin type 1A knockout mice. Hypertension. 2010;56(1):143–50.
Ono M, Murakami T, Kudo A, Isshiki M, Sawada H, Segawa A. Quantitative comparison of anti-fading mounting media for confocal laser scanning microscopy. The journal of histochemistry and cytochemistry : official journal of the Histochemistry Society. 2001;49(3):305–12.
Coleman JL, Ngo T, Schmidt J, Mrad N, Liew CK, Jones NM, et al. Metalloprotease cleavage of the N terminus of the orphan G protein-coupled receptor GPR37L1 reduces its constitutive activity. Sci Signal. 2016;9(423):ra36.
Smith NJ, Ward RJ, Stoddart LA, Hudson BD, Kostenis E, Ulven T, et al. Extracellular loop 2 of the free fatty acid receptor 2 mediates allosterism of a phenylacetamide ago-allosteric modulator. Mol Pharmacol. 2011;80(1):163–73.
Pacher P, Nagayama T, Mukhopadhyay P, Batkai S, Kass DA. Measurement of cardiac function using pressure-volume conductance catheter technique in mice and rats. Nat Protoc. 2008;3(9):1422–34.
Li M, Naqvi N, Yahiro E, Liu K, Powell PC, Bradley WE, et al. c-kit is required for cardiomyocyte terminal differentiation. Circ Res. 2008;102(6):677–85.
Butz GM, Davisson RL. Long-term telemetric measurement of cardiovascular parameters in awake mice: a physiological genomics tool. Physiol Genomics. 2001;5(2):89–97.
Tejada T, Tan L, Torres RA, Calvert JW, Lambert JP, Zaidi M, et al. IGF-1 degradation by mouse mast cell protease 4 promotes cell death and adverse cardiac remodeling days after a myocardial infarction. Proc Natl Acad Sci U S A. 2016;113(25):6949–54.
Consortium GT. Human genomics. The Genotype-Tissue Expression (GTEx) pilot analysis: multitissue gene regulation in humans. Science. 2015;348(6235):648–60.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.
Jolly S, Bazargani N, Quiroga AC, Pringle NP, Attwell D, Richardson WD, et al. G protein-coupled receptor 37-like 1 modulates astrocyte glutamate transporters and neuronal NMDA receptors and is neuroprotective in ischemia. Glia. 2017;
Reddy AS, O'Brien D, Pisat N, Weichselbaum CT, Sakers K, Lisci M, et al. A comprehensive analysis of cell type-specific nuclear RNA from neurons and glia of the brain. Biol Psychiatry. 2017;81(3):252–64.
Barsha G, Denton KM, Mirabito Colafella KM. Sex- and age-related differences in arterial pressure and albuminuria in mice. Biol Sex Differ. 2016;7:57.
Liao Y, Takashima S, Maeda N, Ouchi N, Komamura K, Shimomura I, et al. Exacerbation of heart failure in adiponectin-deficient mice due to impaired regulation of AMPK and glucose metabolism. Cardiovasc Res. 2005;67(4):705–13.
Adachi Y, Saito Y, Kishimoto I, Harada M, Kuwahara K, Takahashi N, et al. Angiotensin II type 2 receptor deficiency exacerbates heart failure and reduces survival after acute myocardial infarction in mice. Circulation. 2003;107(19):2406–8.
Kuroda J, Ago T, Matsushima S, Zhai P, Schneider MD, Sadoshima J. NADPH oxidase 4 (Nox4) is a major source of oxidative stress in the failing heart. Proc Natl Acad Sci U S A. 2010;107(35):15565–70.
Jones DW, Hall JE. Seventh report of the Joint National Committee on prevention, detection, evaluation, and treatment of high blood pressure and evidence from new hypertension trials. Hypertension. 2004;43(1):1–3.
Crowley SD, Gurley SB, Herrera MJ, Ruiz P, Griffiths R, Kumar AP, et al. Angiotensin II causes hypertension and cardiac hypertrophy through its receptors in the kidney. Proc Natl Acad Sci U S A. 2006;103(47):17985–90.
Rowlands DB, Ireland MA, Glover DR, McLeay RA, Stallard TJ, Littler WA. The relationship between ambulatory blood pressure and echocardiographically assessed left ventricular hypertrophy. Clin Sci. 1981;61 Suppl 7:101s–3s.
Xue B, Pamidimukkala J, Hay M. Sex differences in the development of angiotensin II-induced hypertension in conscious mice. Am J Phys Heart Circ Phys. 2005;288(5):H2177–84.
Norton GR, Woodiwiss AJ, Gaasch WH, Mela T, Chung ES, Aurigemma GP, et al. Heart failure in pressure overload hypertrophy. The relative roles of ventricular remodeling and myocardial dysfunction. J Am Coll Cardiol. 2002;39(4):664–71.
Gheorghiade M, Abraham WT, Albert NM, Greenberg BH, O'Connor CM, She L, et al. Systolic blood pressure at admission, clinical characteristics, and outcomes in patients hospitalized with acute heart failure. JAMA. 2006;296(18):2217–26.
Villari B, Campbell SE, Schneider J, Vassalli G, Chiariello M, Hess OM. Sex-dependent differences in left ventricular function and structure in chronic pressure overload. Eur Heart J. 1995;16(10):1410–9.
Travers JG, Kamal FA, Robbins J, Yutzey KE, Blaxall BC. Cardiac fibrosis: the fibroblast awakens. Circ Res. 2016;118(6):1021–40.
Small EM, Olson EN. Pervasive roles of microRNAs in cardiovascular biology. Nature. 2011;469(7330):336–42.
Fermin DR, Barac A, Lee S, Polster SP, Hannenhalli S, Bergemann TL, et al. Sex and age dimorphism of myocardial gene expression in nonischemic human heart failure. Circ Cardiovasc Genet. 2008;1(2):117–25.
Mouat MA, Coleman JLJ, Smith NJ. GPCRs in context: sexual dimorphism in the cardiovascular system. Brit J Pharmacol. 2018; In Press
Meyer RC, Giddens MM, Schaefer SA, Hall RA. GPR37 and GPR37L1 are receptors for the neuroprotective and glioprotective factors prosaptide and prosaposin. Proc Natl Acad Sci U S A. 2013;110(23):9529–34.
Smith NJ. Drug discovery opportunities at the endothelin B receptor-related orphan G protein-coupled receptors, GPR37 and GPR37L1. Front Pharmacol. 2015;6:275.
Shi L, Zhang Z, Su B. Sex biased gene expression profiling of human brains at major developmental stages. Sci Rep. 2016;6:21181.
Mandillo S, Golini E, Marazziti D, Di Pietro C, Matteoni R, Tocchini-Valentini GP. Mice lacking the Parkinson’s related GPR37/PAEL receptor show non-motor behavioral phenotypes: age and gender effect. Genes Brain Behav. 2013;12(4):465–77.
Saleh MC, Connell BJ, Saleh TM. Autonomic and cardiovascular reflex responses to central estrogen injection in ovariectomized female rats. Brain Res. 2000;879(1–2):105–14.
Heximer SP, Knutsen RH, Sun X, Kaltenbronn KM, Rhee MH, Peng N, et al. Hypertension and prolonged vasoconstrictor signaling in RGS2-deficient mice. J Clin Invest. 2003;111(8):1259.
Godecke A, Decking UK, Ding Z, Hirchenhain J, Bidmon HJ, Godecke S, et al. Coronary hemodynamics in endothelial NO synthase knockout mice. Circ Res. 1998;82(2):186–94.
Ledent C, Vaugeois JM, Schiffmann SN, Pedrazzini T, El Yacoubi M, Vanderhaeghen JJ, et al. Aggressiveness, hypoalgesia and high blood pressure in mice lacking the adenosine A2a receptor. Nature. 1997;388(6643):674–8.
Madeddu P, Milia AF, Salis MB, Gaspa L, Gross W, Lippoldt A, et al. Renovascular hypertension in bradykinin B2-receptor knockout mice. Hypertension. 1998;32(3):503–9.
Ngo T, Ilatovskiy AV, Stewart AG, Coleman JL, McRobb FM, Riek RP, et al. Orphan receptor ligand discovery by pickpocketing pharmacological neighbors. Nat Chem Biol. 2017;13(2):235–42.
Leng N, Gu G, Simerly RB, Spindel ER. Molecular cloning and characterization of two putative G protein-coupled receptors which are highly expressed in the central nervous system. Brain Res Mol Brain Res. 1999;69(1):73–83.
Uhlen M, Fagerberg L, Hallstrom BM, Lindskog C, Oksvold P, Mardinoglu A, et al. Proteomics. Tissue-based map of the human proteome. Science. 2015;347(6220):1260419.
Guyenet PG. The sympathetic control of blood pressure. Nat Rev Neurosci. 2006;7(5):335–46.
O'Callaghan EL, Bassi JK, Porrello ER, Delbridge LM, Thomas WG, Allen AM. Regulation of angiotensinogen by angiotensin II in mouse primary astrocyte cultures. J Neurochem. 2011;119(1):18–26.
Stornetta RL, Hawelu-Johnson CL, Guyenet PG, Lynch KR. Astrocytes synthesize angiotensinogen in brain. Science. 1988;242(4884):1444–6.
Bradley DJ, Pascoe JP, Paton JF, Spyer KM. Cardiovascular and respiratory responses evoked from the posterior cerebellar cortex and fastigial nucleus in the cat. J Physiol. 1987;393:107–21.
Bradley DJ, Ghelarducci B, Paton JF, Spyer KM. The cardiovascular responses elicited from the posterior cerebellar cortex in the anaesthetized and decerebrate rabbit. J Physiol. 1987;383:537–50.
Korner PI. The central nervous system and physiological mechanisms of “optimal” cardiovascular control. The Australian journal of experimental biology and medical science. 1971;49(4):319–43.
Korner PI. The phenotypic patterns of essential hypertension are the key to identifying “high blood pressure” genes. Physiol Res. 2010;59(6):841–57.
Xu F, Frazier DT. Role of the cerebellar deep nuclei in respiratory modulation. Cerebellum. 2002;1(1):35–40.
Menuet C, Le S, Dempsey B, Connelly AA, Kamar JL, Jancovski N, et al. Excessive respiratory modulation of blood pressure triggers hypertension. Cell Metab. 2017;25(3):739–48.
We thank the VCCRI BioCORE for their assistance with animal breeding and husbandry.
This work was supported in part by NHMRC & NHF CJ Martin and Future Leader Fellowships (NJS), a St Vincent’s Clinic Foundation grant-in-aid (NJS and RMG), NHMRC Program Grants 573732 and 1074386 (RMG), and Australian Postgraduate Scholarships to JLJC, MAM, and TN. JLJC is additionally supported by a Simon and Michal Wilkenfeld Scholarship. NJS thanks the Mostyn Family Foundation for philanthropic support.
Availability of data and materials
The datasets generated and analyzed during the current study are available from the corresponding author on reasonable request.
All animal work was performed in accordance with the Australia Code for the Care and Use of Animals for Scientific Purposes, 8th Edition (2013), and was approved by the relevant Animal Ethics Committees (Garvan Institute for Medical Research/St Vincent’s Hospital, project number AEC 13/30; University of Melbourne, project number 1212573).
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Coleman, J.L.J., Mouat, M.A., Wu, J. et al. Orphan receptor GPR37L1 contributes to the sexual dimorphism of central cardiovascular control. Biol Sex Differ 9, 14 (2018). https://doi.org/10.1186/s13293-018-0173-y