- Research article
- Open Access
- Open Peer Review
The impact of the rs8005161 polymorphism on G protein-coupled receptor GPR65 (TDAG8) pH-associated activation in intestinal inflammation
BMC Gastroenterology volume 19, Article number: 2 (2019)
Tissue inflammation in inflammatory bowel diseases (IBD) is associated with a decrease in local pH. The gene encoding G-protein-coupled receptor 65 (GPR65) has recently been reported to be a genetic risk factor for IBD. In response to extracellular acidification, proton activation of GPR65 stimulates cAMP and Rho signalling pathways. We aimed to analyse the clinical and functional relevance of the GPR65 associated single nucleotide polymorphism (SNP) rs8005161.
1138 individuals from a mixed cohort of IBD patients and healthy volunteers were genotyped for SNPs associated with GPR65 (rs8005161, rs3742704) and galactosylceramidase (rs1805078) by Taqman SNP assays. 2300 patients from the Swiss IBD Cohort Study (SIBDC) were genotyped for rs8005161 by mass spectrometry based SNP genotyping. IBD patients from the SIBDC carrying rs8005161 TT, CT, CC and non-IBD controls (CC) were recruited for functional studies. Human CD14+ cells were isolated from blood samples and subjected to an extracellular acidic pH shift, cAMP accumulation and RhoA activation were measured.
In our mixed cohort, but not in SIBDC patients, the minor variant rs8005161 was significantly associated with UC. In SIBDC patients, we observed a consistent trend in increased disease severity in patients carrying the rs8005161-TT and rs8005161-CT alleles. No significant differences were observed in the pH associated activation of cAMP production between IBD (TT, CT, WT/CC) and non-IBD (WT/CC) genotype carriers upon an acidic extracellular pH shift. However, we observed significantly impaired RhoA activation after an extracellular acidic pH shift in IBD patients, irrespective of the rs8005161 allele.
The T allele of rs8005161 might confer a more severe disease course in IBD patients. Human monocytes from IBD patients showed impaired pH associated RhoA activation upon an acidic pH shift.
Inflammatory bowel disease (IBD) is characterized by chronic inflammation of the intestinal tract. It is widely accepted that it is a multifactorial inflammatory disease determined by an interaction between genetic and environmental triggers. Inflammation is associated with an increase in local proton concentration and lactate production , with subsequent pro-inflammatory cytokine production. An acidic environment may affect the progression and resolution of inflammation [2,3,4]. To maintain pH homeostasis, cells are required to sense acidic changes in their environment and respond accordingly. Recently, three G protein-coupled receptors (GPCR), T cell death-associated gene 8 (TDAG8 also known as GPR65), ovarian cancer G protein-coupled receptor 1 (OGR1 also known as GPR68) and G protein coupled receptor 4 (GPR4), have been shown to sense extracellular protons and stimulate a variety of signalling pathways [5,6,7,8,9]. Accumulating evidence indicates that these particular proton-sensing receptors play a crucial role in pH homeostasis [6, 8, 10,11,12].
Studies in IBD genetics have identified more than 240 regions in the human genome that increase the risk of IBD [13,14,15,16]. The majority of IBD specific single nucleotide polymorphisms (SNPs) confer an increased risk for both CD and UC . Some of the CD specific genes are associated with bacterial response genes (e.g. NOD2 and autophagy), while for UC there are several specific genes that are associated with immune response and barrier function . More than 70% of the IBD-risk loci are shared with other immune-mediated inflammatory diseases . Genome wide association studies (GWAS) have identified a locus within the GPR65 (TDAG8) gene as one of the risk loci associated with CD and UC [13, 14]. GPR65 is highly expressed in spleen, thymus, lymph nodes and peripheral blood leukocytes, suggesting an important immune response function, which in turn plays a crucial role during the pathogenesis of IBD. In response to extracellular acidic pH, GPR65 activates the adenylyl cyclase (AC)/cAMP/Protein Kinase A (PKA) pathway through Gs proteins [5, 18] and the Rho signalling pathway via G12/13 [5, 7]. Inflammatory processes in the gut are frequently associated with a decrease in local pH, potentially explaining the contribution of GPR65 to intestinal inflammation in IBD.
Determining how genetic polymorphisms can affect functionality is currently a major challenge . In this study, we aimed to test the hypothesis that genetic polymorphisms lead to an altered activity of GPR65, which may result in a higher risk for gut inflammation and IBD. We found that the GPR65 rs8005161 polymorphism has a significant association with UC. Moreover, RhoA activation in human macrophages was significantly lower in acidic conditions for IBD patients vs. non-IBD group.
A second cohort was obtained from the Swiss Inflammatory Bowel Disease Cohort Study (SIBDCS), which includes patients with IBD from all regions of Switzerland since 2006 . The cohort goals and methodology are described elsewhere . We included 2300 adult IBD patients that were enrolled in the study and previously genotyped for the risk variant rs8005161 within the GPR65/GALC gene locus. Genotyping was performed as part of an analysis of the whole Swiss IBD cohort for all SNPs that were known to be associated with IBD at that point in time . Patients with IBD were recruited at the centres participating in SIBDCS . Genotyping of SIBDCS patients was performed using MALDI-TOFF mass spectrometry based SNP genotyping .
Eight SIBDC patients carrying rs8005161-CC, 9 SIBDC patients carrying rs8005161-CT and 9 SIBDC patients carrying rs8005161-TT provided blood samples for cAMP and RhoA assays. Demographic and clinical data were obtained at the time of the blood collection. Ten healthy volunteers were recruited as controls. All controls were rs8005161-CC. One individual in the control group used nonsteroidal anti-inflammatory drugs (NSAIDs) in low dosage on a regular basis. For this individual, cAMP values were in the medium range of the healthy volunteers and RhoA values could not be determined due to low cell numbers.
Isolation of CD14+ human peripheral blood monocytes
Human peripheral blood mononuclear cells (PBMCs) were isolated by density gradient centrifugation (Ficoll Histopaque 10,771 SIGMA, USA) and cryopreserved in foetal calf serum (FCS, Gibco, Thermo Fisher Scientific) supplemented with 10% dimethyl sulfoxide. Upon thawing, cell purification was then performed using EasySep™ Human Monocyte CD14 Enrichment Kit (17,858, Stemcell, Canada) according to the manufacturer’s instructions. Monocytes purity was > 85% as assessed by allophycocyanin (APC)-labelled anti-CD14 (#17–0149-42, eBioscience, USA) and Pacific Blue (PB)-labelled anti-CD45 (#304022, Biolegend, USA) by flow cytometry (Additional file 1: Figure S1).
Genomic DNA extraction
Genomic DNA was isolated from EDTA-blood or intestinal biopsies using the QIAamp DNA Mini Kit (QIAGEN, Hombrechtikon, Switzerland), or QIAzol (Qiagen), respectively, according to manufacturer’s instructions. The concentrations of genomic DNA were quantified using a NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies, Germany).
Genotyping of SNPs was performed using TaqMan allelic discrimination assays (TaqMan SNP Genotyping Assays C_1928636_10, C_1928640_1_ and C_11667238_10 for the SNPs rs8005161, rs3742704 and rs1805078, respectively, all from Applied Biosystems, Thermo Fisher Scientific, USA) on a 7900HT Fast Real-Time PCR instrument (Applied Biosystems, Rotkreuz, Switzerland) using the following cycling conditions: 10 min at 95 °C, 45 cycles of 95 °C for 15 s, and 60 °C for 1 min. Ten nanograms of each genomic DNA was used per PCR reaction in a volume of 5 μl.
The presence of either major or minor alleles of each subject participating in the cAMP and RhoA functional GPR65 assays was confirmed by Taqman genotyping (rs8005161, C_1928636_10 TaqMan SNP Genotyping Assay).
Monocytes were subjected to an extracellular acidic shift from pH 7.6, where pH-sensing GPR65 receptor is minimally active or almost silent, to pH 6.6, where it is maximally active. Cells treated at pH 7.6 served as negative controls. pH shift experiments measuring cAMP production were carried out in Hank’s Balanced Salt Solution (HBSS, 14065056, Gibco) with 25 mM HEPES (Gibco) in a 37 °C incubator without CO2. For the measurement of RhoA activation, pH shift experiments were carried out in serum free RPMI 1640 medium containing a bicarbonate buffer, supplemented with 2 mM Glutamax (35050–038, Gibco) and 25 mM HEPES. The pH of all solutions was adjusted using a calibrated pH meter (Metrohm, Herisau, Switzerland) by the addition of appropriate quantities of NaOH or HCl. Because the pH adjusted RPMI medium also contains bicarbonate buffer, we allowed the media to equilibrate in a 5% CO2 incubator for at least 36 h before use. Control experiments confirmed that the media pH was stable for at least one month under these conditions. The pH was checked before each experiment and found to be stable within a very narrow range (+/− 0.03). RhoA activity assays were incubated in a 5% CO2 humidified 37 °C incubator. All data presented in this paper are referenced to pH measured at room temperature.
cAMP accumulation following the activation of GPR65 by acidic pH was measured by a competitive cell-based sandwich immunoassay and quantified by homogenous time-resolved fluorescence (HTRF) technology (cAMP Dynamic 62AM4PEC, CisBio, France). Human CD14+ monocytes were seeded at non-activating pH 7.6 HBSS supplemented with 25 mM HEPES in 384-well plates (Cat. No. 781080, Greiner) at 10,000 cells/well in HBSS at pH 7.6 with or without the GPR65 antagonist (10uM) (provided by Novartis Institutes for Biomedical Research, Switzerland), and incubated for 15 min, followed by a 30 min pH shift, which was achieved by addition of the appropriate amount of HBSS buffer to obtain the desired final pH (pH 7.6 or pH 6.5). All incubations were carried out in a non-CO2 incubator at 37 °C. The optimal activating pH for GPR65/cAMP –mediated signalling was determined in extensive pH dose response experiments (Additional file 2: Figure S2), and described in detail in the following section. Phosphodiesterase inhibitors (1 mM IBMX, 10 μM Rolipram, 1 μM BAY) were used in all conditions. Samples and the cAMP standards were analysed using a sigmoidal dose response model with variable slope in using the software package GraphPad Prism, La Jolla California USA, www.graphpad.com (San Diego, CA, USA).
cAMP activation assay validation
Proton-activated GPR65/cAMP –mediated signalling was tested using the human monocytic cell line THP-1 and CD14+ primary human monocytes isolated from non-IBD subjects (WT/CC genotype). To confirm that pH was associated with GPR65/cAMP –mediated signalling, pH dose response experiments were performed. Cells were starved at pH 7.6 for 2 h and subsequently subjected to a pH shift for 10 min (pH 6.2 to 7.8 with 0.2 increments). The highest cAMP accumulation was observed at pH 6.4–6.8, whereas only low cAMP concentrations were observed at pH 7.6–7.8 (Additional file 2: Figure S2). Maximum activation and inactivation of GPR65 were achieved at pH 6.5 and pH 7.6, respectively.
RhoA GTPase activation assay
Human CD14+ cells were plated in RPMI plus 10% FCS and incubated for 1 h, followed by 2 h of starvation at pH 7.6 in serum free RPMI. The pH shift was performed for 10 min at pH 6.6, with the negative control (pH 7.6). All incubations were carried out in a 5% CO2 humidified 37 °C incubator. The pH range was established in pH dose response experiments (Additional file 3: Figure S3). 15 μg of protein was loaded per well and GTP-bound RhoA protein levels were measured in duplicates according to the manufacturer’s instructions (# BK124, Cytoskeleton, USA). Final absorbance (OD490) was measured in a Synergy 2 micro-plate reader (Biotek, Luzern, Switzerland). No internal standard for RhoA has been established, therefore baseline values for RhoA cannot be compared between experiments and only normalized values are presented.
RhoA activation assay validation
To confirm that RhoA activation is proton dependent, THP-1 cells and CD14+ monocytes were subjected to a pH shift (10 min) at different pH levels, with a preliminary starvation step (2 h) at non-activating pH (pH 7.6). As shown in Additional file 3: Figure S3, pH 6.6 elicited a significant increase in RhoA activation compared to pH 6.2, 7.4 and 7.6, which, in contrast, induced no significant activation in CD14+ monocytes and THP-1 cells. Since GPR65/G12/13/RhoA signalling was highest at pH 6.6, this pH was used in all RhoA activity assays.
Normalization of RhoA and cAMP levels in IBD patients and controls
cAMP and RhoA activity was measured at pH 6.5 and 6.6, respectively and all levels were normalized to RhoA and cAMP levels from the same participant tested at conditions with lowest RhoA and cAMP production. This residual activity (at non-activating pH and in the presence of the inhibitor) is unlikely to be due to GPR65. For cAMP, we normalized to levels at pH 7.6 in the presence of a GPR65 inhibitor. For RhoA activity we normalized to levels at pH 7.6. Due to a high demand of cells from human patients for the RhoA assay, the control at pH 7.6 with the GPR65 inhibitor was not feasible.
Clinical data were retrieved from the data centre of the SIBDCS at the Lausanne University Hospital. These data were entered into a database (Access 2000; Microsoft Switzerland Ltd., Liab., Co., Wallisellen, Switzerland). The Statistical Package for the Social Sciences (version 21; SPSS, Chicago, Ill., USA), GraphPad Prism, version 7, or Stata software (StataCorp., 2015. Stata Statistical Software: Release 14. College Station, TX, USA) was used for the statistical analysis.
The Chi-square test or Fisher’s exact test were used to determine associations between individual SNPs and subject phenotypes. Groups of data were compared using Kruskal-Wallis test, Fisher’s exact, student t test or one-way ANOVA. Data are presented as mean and interquartile range (IQR). Probabilities (p, two tailed) of p < 0.05 were considered statistically significant. Throughout this manuscript, asterisks denote significant differences at *, p < 0.05; **, p < 0.01; ***, p < 0.001; ****, p < 0.0001.
Genotyping for rs8005161, rs3742704 and rs1805078 in IBD vs. non-IBD patients by Taqman SNP assays
Genomic DNA from a group of 547 patients with IBD (mean age 55.7 years, range 20–81) and 591 non-IBD subjects (mean age 42.6 years, range 16–82) was used to genotype individuals for genetic variants rs8005161 [13, 14], rs3742704 [12, 13] (both assigned to the GPR65 gene), and an additional neighbouring SNP rs1805078 (galactosylceramidase, GALC). Frequencies of homozygous and heterozygous carriers of the major and minor alleles, respectively are provided in Table 1. For rs8005161 (GPR65), the frequency of individuals homozygous for the minor T allele was significantly higher in UC individuals compared to individuals without IBD (p < 0.05). When comparing the number of individuals carrying at least one T allele (i.e. CT + TT vs. CC) we noticed a trend for higher frequency of the T carriers in UC patients which did not reach significance. In contrast, in individuals with CD the frequency of rs8005161 alleles did not differ from non-IBD carriers. There were no significant differences in the genotype and allele frequencies distribution between the groups for the GPR65 variant rs3742704 and the GALC variant rs1805078. In a complementary analysis, we compared the total number of rs8005161 minor alleles which differed significantly between UC patients and healthy volunteers but not when CD patients or the other SNPs were evaluated (Additional file 4: Table S1).
When patients with UC and CD were compared, no difference between individuals carrying at least one T allele (i.e. CT + TT vs. CC) or number of T alleles was found (p = 0.52; p = 0.25).
Genotyping of Swiss IBD cohort study (SIBDC) patients for GPR65 IBD risk SNP rs8005161
The SIBDC comprises a data base of clinical and genetic data of patients with IBD. Genotype data for rs8005161 (GPR65) for 2300 SIBDCS patients are available; rs3742704 (GPR65) and rs1805078 (GALC) have not been genotyped in SIBDC. No control cohort is available for SIBDC. For the SNP rs8005161 (major allele: C, minor allele: T), 28/2300 patients (1.2%) were homozygous TT carriers, 430/2300 (18.7%) heterozygous CT carriers, and 1842/2300 (80.1%) were wildtype (WT) CC patients. Allele frequencies for T and C genotype were 10.6% (426) and 89.4% (3710), respectively (Table 2). Overall, carriers of the T allele tended to be diagnosed more often with CD (NS). Thus, the possible association of the T allele of rs8005161 with UC in our initial cohort could not be confirmed with SIBDC data.
Interestingly, we noted a consistent trend for a more severe disease course with an earlier age at diagnosis (NS), a lower body mass index (p = 0.02), a higher frequency of surgery (NS) and a higher frequency of therapy with biologics (p = 0.02). A similar trend for body mass index (BMI) and past or current therapy with biologics was confirmed when the CC, CT and TT genotypes were considered separately (Additional file 5: Table S2). However, for intestinal surgery this trend was not confirmed, and likely due to the small size of the TT group.
pH-dependent activation of cAMP production is not affected by the rs8005161 T-variant in CD14+ cells
Next, we sought to determine whether the GPR65 SNP variant rs8005161 has an impact on GPR65 protein functionality and downstream signalling cascades in vitro, by measuring production of cAMP and activation of RhoA, two well-established second massagers downstream of GPR65. To validate pH-dependent GPR65 activation, we performed pH dose response experiments in a human monocyte cell line (THP-1) and in CD14+ human monocytes, isolated from healthy volunteers (Additional file 1: Figure S1). Our data confirm maximum of cAMP and RhoA activity at pH 6.5 and 6.6, respectively (Additional file 2: Figure S2 and Additional file 3: Figure S3). For subsequent experiments, these pH values were used and normalized to levels from the same participant at conditions with the least activity.
We recruited patients from the SIBDC with TT, CT and CC genotype based on the results from SIBDCS sequencing, in addition to healthy CC volunteers, as described in the Materials and Methods section. Patient characteristics are shown in Table 3. Patients with different alleles of rs8005161 were of similar age and had a similar distribution of UC and CD; however, in our sample cohort, individuals with the TT genotype had lower disease activity with a lower frequency of TNF-inhibitor or Vedolizumab usage.
As shown in Fig. 1, normalized cAMP level increased in PBMCs from all groups upon an extracellular acidic pH shift from pH 7.6 to pH 6.5. In the presence of the GPR65 antagonist, no activation of cAMP production was observed upon pH shift, confirming that pH-stimulated cAMP production is mediated by GPR65. No difference between individuals with CC, CT and TT genotype could be detected (p > 0.05, data not shown), indicating that the genotype variant of GPR65 rs8995161 has no influence on the GPR65/Gs/cAMP pathway in human monocytes under the conditions tested.
We also note differences in the residual cAMP activity in the presence of the GPR65 inhibitor (Additional file 6: Figure S4). Since these differences were observed in the presence of the GPR65 inhibitor at the non-activating pH 7.6, they are unlikely to represent GPR65 activity and the underlying physiology remains unclear.
pH dependent activation of rho a production is significantly decreased in IBD patients
The GPR65 receptor is also coupled to G12/13 proteins, which are known to mediate small GTPase RhoA activation by phosphorylation. Consequently, we next investigated the effects of the GPR65 rs8005161 variant on this signalling pathway. An acidic pH shift from pH 7.6 to pH 6.6 resulted in activation of RhoA. Interestingly, this stimulation was lower for IBD patients compared to the healthy control subjects (p < 0.05, Fig. 2a). No differential RhoA GTPase activation in CD14+ monocytes according to the rs8005161 genotype could be detected (Fig. 2b). Experiments with the GPR65 antagonist were not feasible due to the limited amount of human material. Thus, CD14+ monocytes from IBD patients demonstrate reduced RhoA activation upon acidic pH shift compared to healthy control subjects.
In this study, we addressed the clinical relevance of the SNP rs8005161 GPR65 variant as a risk gene for IBD, and tested a potential functional consequence of this SNP variant in human CD14+ monocytes/ macrophages. We found several indications for a more severe clinical phenotype in IBD patients with the T allele of the GPR65 SNP rs8005161 in patients from the SIBDC. We examined proton-activated TDAG8 -mediated signalling pathways in CD14+ monocytes from IBD rs8005161 (WT/CC, CT, TT) and non-IBD (WT/CC) carriers. We observed lower RhoA GTPase activation upon an acidic pH shift in IBD patients compared to healthy volunteers. No differential activation of either RhoA or cAMP stimulated by acidosis was detected in individuals with different rs8005161 genotypes, thus rendering major effects of the CC, CT or TT variants on cAMP or RhoA activation under the conditions tested unlikely.
The GPR65 gene encodes a transmembrane receptor, which is reported to function as a proton sensor [5, 18, 24, 25]. Previously, psychosine (1-β-D-galactosylsphingosine) was proposed to activate GPR65 , however the study showed no data for the specific interaction of psychosine with GPR65, and subsequent reports have not supported this finding [5, 18, 24, 25]. Recently, BTB09089 (3-[(2,4-dichlorobenzyl)thio]-1,6-dimethyl-5,6-dihydro-1H-pyridazino[4,5 e][1,3,4]thiadiazin-5-one) has been reported to be an allosteric modulator for GPR65 . However, despite numerous studies, the physiological role of GPR65 under inflammatory conditions (acidic pH) is unclear. GPR65 is mainly expressed in immune cells [28, 29], suggesting an immunological role. Proton activation of GPR65 stimulates second messengers Gs/cAMP and G12/13/RhoA [5, 7, 30]. Second messenger cAMP is produced by the activation of adenylyl cyclase, which converts adenosine triphosphate (ATP) into the biologically active signalling mediator, with subsequent degradation by phosphodiesterases (PDEs). The role of cAMP in regulating inflammatory diseases has been well studied [31, 32] and there is continued interest in cAMP as a therapeutic target to treat inflammatory diseases . Increased cAMP levels can inhibit the secretion of proinflammatory cytokines and chemokines [27, 30, 34,35,36,37], inhibit inflammatory cell migration [38,39,40] and modulate epithelial barrier formation [41,42,43]. In the present study, we observed TDAG8–mediated increases in cAMP at reduced pH in peripheral monocytes. We could confirm TDAG8 signalling by using a specific TDAG8 inhibitor. However, no differences between IBD variants and non-IBD subjects were detected, suggesting that major effects of cAMP signalling at the conditions tested unlikely.
There is increasing interest in small GTPases of the Rho family as candidates for therapeutic intervention due to their involvement in a wide variety of diseases. Rho functions as a molecular switch, in the GTP-bound conformation, the proteins are able to interact with their downstream targets and transmit signals to the cell . Rho GTPases are critical regulators of many cellular functions including cytoskeletal remodelling, cell-cell adhesion, cell polarization, vesicle trafficking, morphogenesis, apoptosis, cell migration, organelle development, membrane transport pathways and tumour motility and proliferation [44,45,46,47,48,49,50,51]. Our study identified a lower RhoA activation in response to acidic pH from monocytes of IBD patients compared to healthy volunteers, but no differences between the rs8005161 GPR65 CC, CT or TT variants. These data are in line with a role of pH sensing for immune activation in IBD pathogenesis. The reasons for this differential activation remain unclear and possible explanations include a higher baseline activation in IBD patients with (subclinical) disease. Alternatively, cell migration and wound healing may be impaired in IBD patients. Thus, future studies should address this issue with more mechanistic studies.
Our analysis did not identify genotype dependent changes in either RhoA or cAMP activation. There might be several explanations for these findings: i) differential activation might only be relevant in subsets of macrophages or other cell types. ii) differential effects of the genotypes might not be mediated via cAMP or RhoA but other properties of GPR65. One example would be genotype dependent differential inhibition of Rac1. Inhibition of Rac1 via the via G12/13 pathway has as described previously  and Rac1 inhibition has been associated with remission in IBD . iii) finally, considering small effects of rs8005161 polymorphisms  (increase of CD risk by an odds ratio (OR) of 1.16 (1.09–1.22) and UC by an OR of 1.14 (1.08–1.21)) our study with a limited number of participants might be underpowered to detect small differences in the activation of secondary messengers.
Most IBD risk genes increase the risk of both, UC and CD in carriers . In a recent meta-analysis, OR for carriers of rs8005161 for CD was 1.156 (CI: 1.092–1.222), the OR for UC was 1.143 (CI: 1.076–1.213) . In agreement with a general increase in the risk for IBD, no differences between allelic frequencies of rs8005161 between CD and UC patients was found in our study.
Associations of IBD genotypes with subsequent disease course have been difficult to establish. However, disease location in CD (ileal vs. colonic) could be predicted by a compound genetic risk score . As expected, for rs8005161 no association with disease location was observed in agreement with highly similar OR for CD (predicting ileal disease) and UC (predicting colonic disease). Furthermore, genetic predisposition can predict onset of disease and individuals with a high genetic burden developed CD 5 years earlier than individuals with the lowest genetic risk . In agreement with this notion, carriers of the rs8005161 T allele developed disease 1.3 years before rs8005161-CC carriers; however, probably due to the wide distribution of the year of onset of disease, this difference failed to reach statistical significance. Associations of genetic risks and clinical course or treatment of IBD are insufficiently understood; most likely due to limited precision of data recorded in large IBD cohort data bases. Analysis of the SIBDC data evidenced an association of need for more intense treatment, i.e. biological therapy (p = 0.02) and a trend for more intestinal surgery (p = 0.13) in our patient cohort. However, even nominally significant associations would not remain significant after Bonferroni correction, pointing to the need of large, well-characterized cohorts for studying the important questions of IBD genetic risk and disease course and treatment.
Our study has various strengths and limitations. We were able to use a large prospective cohort of well-characterized IBD patients to test associations of rs8005161 genotypes with disease course. Moreover, we were able to recruit patients with IBD and rare genotypes, thus enabling functional studies with individuals of each genotype for RhoA and cAMP activation. Limitations include the small sample size for all experiments, and lack of non-wildtype carriers of rs8005161 in the control population. Limited amount of patient material prohibited inclusion of the GPR65 inhibitor for RhoA tests. Future studies should also address expression levels of GPR65 RNA and start gene sequencing to test for polymorphisms associated with specific variants of the SNP.
In summary, in the SIBDC patients, there is a trend for a more severe disease course for T allele carriers and a tendency towards a lower probability for anti-TNF treatment. In addition, we confirmed a significant association between the rare homozygote rs8005161 TT genotype and a diagnosis of UC. Our study did not identify biochemical changes in individuals with various genotypes of rs8005161, however we identified a lower activation of RhoA upon an acidic pH shift in IBD patients. Thus, pH sensors may be interesting new targets for pharmacological intervention in intestinal inflammation.
Cyclic adenosine monophosphate
G protein coupled receptor
G protein coupled receptor 65
Genome wide association studies
Hank’s Balanced Salt Solution
Inflammatory bowel disease
Phosphate buffered saline
Swiss IBD Cohort Study
Single nucleotide polymorphisms
T cell death associated gene 8
Lardner A. The effects of extracellular pH on immune function. J Leukoc Biol. 2001;69(4):522–30.
Brokelman WJ, Lensvelt M, Borel Rinkes IH, Klinkenbijl JH, Reijnen MM. Peritoneal changes due to laparoscopic surgery. Surg Endosc. 2011;25(1):1–9.
Hanly EJ, Aurora AA, Shih SP, Fuentes JM, Marohn MR, De Maio A, Talamini MA. Peritoneal acidosis mediates immunoprotection in laparoscopic surgery. Surgery. 2007;142(3):357–64.
Martinez D, Vermeulen M, von Euw E, Sabatte J, Maggini J, Ceballos A, Trevani A, Nahmod K, Salamone G, Barrio M, et al. Extracellular acidosis triggers the maturation of human dendritic cells and the production of IL-12. J Immunol. 2007;179(3):1950–9.
Ishii S, Kihara Y, Shimizu T. Identification of T cell death-associated gene 8 (TDAG8) as a novel acid sensing G-protein-coupled receptor. J Biol Chem. 2005;280(10):9083–7.
Ludwig MG, Vanek M, Guerini D, Gasser JA, Jones CE, Junker U, Hofstetter H, Wolf RM, Seuwen K. Proton-sensing G-protein-coupled receptors. Nature. 2003;425(6953):93–8.
Mogi C, Nakakura T, Okajima F. Role of extracellular proton-sensing OGR1 in regulation of insulin secretion and pancreatic beta-cell functions. Endocr J. 2014;61(2):101–10.
Okajima F. Regulation of inflammation by extracellular acidification and proton-sensing GPCRs. Cell Signal. 2013;25(11):2263–71.
Seuwen K, Ludwig MG, Wolf RM. Receptors for protons or lipid messengers or both? J Recept Signal Transduct Res. 2006;26(5–6):599–610.
Mohebbi N, Benabbas C, Vidal S, Daryadel A, Bourgeois S, Velic A, Ludwig MG, Seuwen K, Wagner CA. The proton-activated G protein coupled receptor OGR1 acutely regulates the activity of epithelial proton transport proteins. Cell Physiol Biochem. 2012;29(3–4):313–24.
Saxena H, Deshpande DA, Tiegs BC, Yan H, Battafarano RJ, Burrows WM, Damera G, Panettieri RA, Dubose TD Jr, An SS, et al. The GPCR OGR1 (GPR68) mediates diverse signalling and contraction of airway smooth muscle in response to small reductions in extracellular pH. Br J Pharmacol. 2012;166(3):981–90.
Lassen KG, McKenzie CI, Mari M, Murano T, Begun J, Baxt LA, Goel G, Villablanca EJ, Kuo SY, Huang H, et al. Genetic coding variant in GPR65 alters lysosomal pH and links lysosomal dysfunction with colitis risk. Immunity. 2016;44(6):1392–405.
Franke A, McGovern DP, Barrett JC, Wang K, Radford-Smith GL, Ahmad T, Lees CW, Balschun T, Lee J, Roberts R, et al. Genome-wide meta-analysis increases to 71 the number of confirmed Crohn's disease susceptibility loci. Nat Genet. 2010;42(12):1118–25.
Jostins L, Ripke S, Weersma RK, Duerr RH, McGovern DP, Hui KY, Lee JC, Schumm LP, Sharma Y, Anderson CA, et al. Host-microbe interactions have shaped the genetic architecture of inflammatory bowel disease. Nature. 2012;491(7422):119–24.
Liu JZ, van Sommeren S, Huang H, Ng SC, Alberts R, Takahashi A, Ripke S, Lee JC, Jostins L, Shah T, et al. Association analyses identify 38 susceptibility loci for inflammatory bowel disease and highlight shared genetic risk across populations. Nat Genet. 2015;47(9):979–86.
de Lange KM, Moutsianas L, Lee JC, Lamb CA, Luo Y, Kennedy NA, Jostins L, Rice DL, Gutierrez-Achury J, Ji SG, et al. Genome-wide association study implicates immune activation of multiple integrin genes in inflammatory bowel disease. Nat Genet. 2017;49(2):256–61.
Lee JC, Espeli M, Anderson CA, Linterman MA, Pocock JM, Williams NJ, Roberts R, Viatte S, Fu B, Peshu N, et al. Human SNP links differential outcomes in inflammatory and infectious disease to a FOXO3-regulated pathway. Cell. 2013;155(1):57–69.
Wang JQ, Kon J, Mogi C, Tobo M, Damirin A, Sato K, Komachi M, Malchinkhuu E, Murata N, Kimura T, et al. TDAG8 is a proton-sensing and psychosine-sensitive G-protein-coupled receptor. J Biol Chem. 2004;279(44):45626–33.
Altshuler D, Daly MJ, Lander ES. Genetic mapping in human disease. Science. 2008;322(5903):881–8.
Attinkara R, Mwinyi J, Truninger K, Regula J, Gaj P, Rogler G, Kullak-Ublick GA, Eloranta JJ, Swiss IBDCSG. Association of genetic variation in the NR1H4 gene, encoding the nuclear bile acid receptor FXR, with inflammatory bowel disease. BMC Res Notes. 2012;5:461.
Scharl M, Mwinyi J, Fischbeck A, Leucht K, Eloranta JJ, Arikkat J, Pesch T, Kellermeier S, Mair A, Kullak-Ublick GA, et al. Crohn's disease-associated polymorphism within the PTPN2 gene affects muramyl-dipeptide-induced cytokine secretion and autophagy. Inflamm Bowel Dis. 2012;18(5):900–12.
Pittet V, Juillerat P, Mottet C, Felley C, Ballabeni P, Burnand B, Michetti P, Vader JP, Swiss IBDCSG. Cohort profile: the Swiss inflammatory bowel disease cohort study (SIBDCS). Int J Epidemiol. 2009;38(4):922–31.
Storm N, Darnhofer-Patel B, van den Boom D, Rodi CP. MALDI-TOF mass spectrometry-based SNP genotyping. Methods Mol Biol. 2003;212:241–62.
Ihara Y, Kihara Y, Hamano F, Yanagida K, Morishita Y, Kunita A, Yamori T, Fukayama M, Aburatani H, Shimizu T, et al. The G protein-coupled receptor T-cell death-associated gene 8 (TDAG8) facilitates tumor development by serving as an extracellular pH sensor. Proc Natl Acad Sci U S A. 2010;107(40):17309–14.
Radu CG, Nijagal A, McLaughlin J, Wang L, Witte ON. Differential proton sensitivity of related G protein-coupled receptors T cell death-associated gene 8 and G2A expressed in immune cells. Proc Natl Acad Sci U S A. 2005;102(5):1632–7.
Im DS, Heise CE, Nguyen T, O'Dowd BF, Lynch KR. Identification of a molecular target of psychosine and its role in globoid cell formation. J Cell Biol. 2001;153(2):429–34.
Onozawa Y, Fujita Y, Kuwabara H, Nagasaki M, Komai T, Oda T. Activation of T cell death-associated gene 8 regulates the cytokine production of T cells and macrophages in vitro. Eur J Pharmacol. 2012;683(1–3):325–31.
Choi J-W, Lee SY, Choi Y. Identification of a putative G protein-coupled receptor induced during activation-induced apoptosis of T cells. Cell Immunol. 1996;168(1):78–84.
Kyaw H, Zeng Z, Su K, Fan P, Shell BK, Carter KC, Li Y. Cloning, characterization, and mapping of human homolog of mouse T-cell death-associated gene. DNA Cell Biol. 1998;17(6):493–500.
Mogi C, Tobo M, Tomura H, Murata N, He XD, Sato K, Kimura T, Ishizuka T, Sasaki T, Sato T, et al. Involvement of proton-sensing TDAG8 in extracellular acidification-induced inhibition of Proinflammatory cytokine production in peritoneal macrophages. J Immunol. 2009;182(5):3243–51.
Serezani CH, Ballinger MN, Aronoff DM, Peters-Golden M. Cyclic AMP: master regulator of innate immune cell function. Am J Respir Cell Mol Biol. 2008;39(2):127–32.
Yan K, Gao LN, Cui YL, Zhang Y, Zhou X. The cyclic AMP signaling pathway: exploring targets for successful drug discovery (review). Mol Med Rep. 2016;13(5):3715–23.
Raker VK, Becker C, Steinbrink K. The cAMP pathway as therapeutic target in autoimmune and inflammatory diseases. Front Immunol. 2016;7:123.
Aoki H, Mogi C, Okajima F. Ionotropic and metabotropic proton-sensing receptors involved in airway inflammation in allergic asthma. Mediat Inflamm. 2014;2014:712962.
Jin Y, Sato K, Tobo A, Mogi C, Tobo M, Murata N, Ishii S, Im DS, Okajima F. Inhibition of interleukin-1beta production by extracellular acidification through the TDAG8/cAMP pathway in mouse microglia. J Neurochem. 2014;129(4):683–95.
Nagasaka A, Mogi C, Ono H, Nishi T, Horii Y, Ohba Y, Sato K, Nakaya M, Okajima F, Kurose H. The proton-sensing G protein-coupled receptor T-cell death-associated gene 8 (TDAG8) shows cardioprotective effects against myocardial infarction. Sci Rep. 2017;7(1):7812.
Onozawa Y, Komai T, Oda T. Activation of T cell death-associated gene 8 attenuates inflammation by negatively regulating the function of inflammatory cells. Eur J Pharmacol. 2011;654(3):315–9.
Lang P, Gesbert F, Delespine-Carmagnat M, Stancou R, Pouchelet M, Bertoglio J. Protein kinase a phosphorylation of RhoA mediates the morphological and functional effects of cyclic AMP in cytotoxic lymphocytes. EMBO J. 1996;15(3):510–9.
Laudanna C, Campbell JJ, Butcher EC. Elevation of intracellular cAMP inhibits RhoA activation and integrin-dependent leukocyte adhesion induced by chemoattractants. J Biol Chem. 1997;272(39):24141–4.
Park PH, Huang H, McMullen MR, Bryan K, Nagy LE. Activation of cyclic-AMP response element binding protein contributes to adiponectin-stimulated interleukin-10 expression in RAW 264.7 macrophages. J Leukoc Biol. 2008;83(5):1258–66.
Hatzelmann A, Morcillo EJ, Lungarella G, Adnot S, Sanjar S, Beume R, Schudt C, Tenor H. The preclinical pharmacology of roflumilast--a selective, oral phosphodiesterase 4 inhibitor in development for chronic obstructive pulmonary disease. Pulm Pharmacol Ther. 2010;23(4):235–56.
Lorenowicz MJ, Fernandez-Borja M, van Stalborch AM, van Sterkenburg MA, Hiemstra PS, Hordijk PL. Microtubule dynamics and Rac-1 signaling independently regulate barrier function in lung epithelial cells. Am J Physiol Lung Cell Mol Physiol. 2007;293(5):L1321–31.
Zimmerman NP, Kumar SN, Turner JR, Dwinell MB. Cyclic AMP dysregulates intestinal epithelial cell restitution through PKA and RhoA. Inflamm Bowel Dis. 2012;18(6):1081–91.
Bustelo XR, Sauzeau V, Berenjeno IM. GTP-binding proteins of the rho/Rac family: regulation, effectors and functions in vivo. Bioessays. 2007;29(4):356–70.
Etienne-Manneville S, Hall A. Rho GTPases in cell biology. Nature. 2002;420(6916):629–35.
Hodge RG, Ridley AJ. Regulating rho GTPases and their regulators. Nat Rev Mol Cell Biol. 2016;17(8):496–510.
Konigs V, Jennings R, Vogl T, Horsthemke M, Bachg AC, Xu Y, Grobe K, Brakebusch C, Schwab A, Bahler M, et al. Mouse macrophages completely lacking rho subfamily GTPases (RhoA, RhoB, and RhoC) have severe lamellipodial retraction defects, but robust chemotactic navigation and altered motility. J Biol Chem. 2014;289(44):30772–84.
Kozasa T, Hajicek N, Chow CR, Suzuki N. Signalling mechanisms of RhoGTPase regulation by the heterotrimeric G proteins G12 and G13. J Biochem. 2011;150(4):357–69.
Vogt S, Grosse R, Schultz G, Offermanns S. Receptor-dependent RhoA activation in G12/G13-deficient cells: genetic evidence for an involvement of Gq/G11. J Biol Chem. 2003;278(31):28743–9.
Wettschureck N, Offermanns S. Mammalian G proteins and their cell type specific functions. Physiol Rev. 2005;85(4):1159–204.
Worzfeld T, Wettschureck N, Offermanns S. G(12)/G(13)-mediated signalling in mammalian physiology and disease. Trends Pharmacol Sci. 2008;29(11):582–9.
Sugimoto N, Takuwa N, Okamoto H, Sakurada S, Takuwa Y. Inhibitory and stimulatory regulation of Rac and cell motility by the G12/13-rho and Gi pathways integrated downstream of a single G protein-coupled sphingosine-1-phosphate receptor isoform. Mol Cell Biol. 2003;23(5):1534–45.
Parikh K, Zhou L, Somasundaram R, Fuhler GM, Deuring JJ, Blokzijl T, Regeling A, Kuipers EJ, Weersma RK, Nuij VJ, et al. Suppression of p21Rac signaling and increased innate immunity mediate remission in Crohn’s disease. Sci Transl Med. 2014;6(233):233ra253.
Cleynen I, Boucher G, Jostins L, Schumm LP, Zeissig S, Ahmad T, Andersen V, Andrews JM, Annese V, Brand S, et al. Inherited determinants of Crohn's disease and ulcerative colitis phenotypes: a genetic association study. Lancet. 2016;387(10014):156–67.
Ananthakrishnan AN, Huang H, Nguyen DD, Sauk J, Yajnik V, Xavier RJ. Differential effect of genetic burden on disease phenotypes in Crohn's disease and ulcerative colitis: analysis of a north American cohort. Am J Gastroenterol. 2014;109(3):395–400.
We thank Solange Vidal and Caterina Safina for assistance in cAMP assays and Christian Hiller for genomic DNA isolation. We thank all patients and the members of the Swiss IBD cohort for their participation in this study.
Members of the SIBDCS study group:
Karim Abdelrahman, Gentiana Ademi, Patrick Aepli, Amman Thomas, Claudia Anderegg, Anca-Teodora Antonino, Eva Archanioti, Eviano Arrigoni, Diana Bakker de Jong, Bruno Balsiger, Polat Bastürk, Peter Bauerfeind, Andrea Becocci, Dominique Belli, José M. Bengoa, Luc Biedermann, Janek Binek, Mirjam Blattmann, Stephan Boehm, Tujana Boldanova, Jan Borovicka, Christian P. Braegger, Stephan Brand, Lukas Brügger, Simon Brunner, Patrick Bühr, Bernard Burnand, Sabine Burk, Emanuel Burri, Sophie Buyse, Dahlia-Thao Cao, Ove Carstens, Dominique H. Criblez, Sophie Cunningham, Fabrizia D’Angelo, Philippe de Saussure, Lukas Degen, Joakim Delarive, Christopher Doerig, Barbara Dora, Susan Drerup, Mara Egger, Ali El-Wafa, Matthias Engelmann, Jessica Ezri, Christian Felley, Markus Fliegner, Nicolas Fournier, Montserrat Fraga, Yannick Franc, Pascal Frei, Remus Frei, Michael Fried, Florian Froehlich, Raoul Ivano Furlano, Luca Garzoni, Martin Geyer, Laurent Girard, Marc Girardin, Delphine Golay, Ignaz Good, Ulrike Graf Bigler, Beat Gysi, Johannes Haarer, Marcel Halama, Janine Haldemann, Pius Heer, Benjamin Heimgartner, Beat Helbling, Peter Hengstler, Denise Herzog, Cyrill Hess, Roxane Hessler, Klaas Heyland, Thomas Hinterleitner, Claudia Hirschi, Petr Hruz, Pascal Juillerat, Carolina Khalid-de Bakker, Stephan Kayser, Céline Keller, Christina Knellwolf (-Grieger), Christoph Knoblauch, Henrik Köhler, Rebekka Koller, Claudia Krieger(-Grübel), Patrizia Künzler, Rachel Kusche, Frank Serge Lehmann, Andrew Macpherson, Michel H. Maillard, Michael Manz, Astrid Marot, Rémy Meier, Christa Meyenberger, Pamela Meyer, Pierre Michetti, Benjamin Misselwitz, Patrick Mosler, Christian Mottet, Christoph Müller, Beat Müllhaupt, Leilla Musso, Michaela Neagu, Cristina Nichita, Jan Niess, Andreas Nydegger, Nicole Obialo, Diana Ollo, Cassandra Oropesa, Ulrich Peter, Daniel Peternac, Laetitia Marie Petit, Valérie Pittet, Daniel Pohl, Marc Porzner, Claudia Preissler, Nadia Raschle, Ronald Rentsch, Alexandre Restellini, Sophie Restellini, Jean-Pierre Richterich, Frederic Ris, Branislav Risti, Marc Alain Ritz, Gerhard Rogler, Nina Röhrich, Jean-Benoît Rossel, Vanessa Rueger, Monica Rusticeanu, Markus Sagmeister, Gaby Saner, Bernhard Sauter, Mikael Sawatzki, Michael Scharl, Martin Schelling, Susanne Schibli, Hugo Schlauri, Dominique Schluckebier, Daniela Schmid, Sybille Schmid (-Uebelhart), Jean-François Schnegg, Alain Schoepfer, Vivianne Seematter, Frank Seibold, Mariam Seirafi, Gian-Marco Semadeni, Arne Senning, Christiane Sokollik, Joachim Sommer, Johannes Spalinger, Holger Spangenberger, Philippe Stadler, Peter Staub, Dominic Staudenmann, Volker Stenz, Michael Steuerwald, Alex Straumann, Bruno Strebel, Andreas Stulz, Michael Sulz, Aurora Tatu, Michela Tempia-Caliera, Joël Thorens, Kaspar Truninger, Radu Tutuian, Patrick Urfer, Stephan Vavricka, Francesco Viani, Jürg Vögtlin, Roland Von Känel, Dominique Vouillamoz, Rachel Vulliamy, Paul Wiesel, Reiner Wiest, Stefanie Wöhrle, Samuel Zamora, Silvan Zander, Tina Wylie, Jonas Zeitz, Dorothee Zimmermann.
This work was supported by research grants to Gerhard Rogler from the Swiss National Science Foundation (SNF) awarded to Gerhard Rogler, Grant № 153380 (http://p3.snf.ch/Project-153380), Grant № 33CS30_148422, Swiss IBD Cohort Study (http://p3.snf.ch/project-148422) and to Benjamin Misselwitz (SNF Grant No. 32473B_156525). The funding institutions had no role in the study design, in the collection, analysis and interpretation of data and in the writing of the manuscript.
Availability of data and materials
Data from the Swiss IBD cohort study are not publicly available since Swiss law does not allow publication of individual data set. To obtain SIBDC data for further analyses, please contact the scientific board of SIBDC (http://www.ibdcohort.ch/index.php/informationen-fuer-forscher.html). The datasets from all other experiments and analyses of the current study are available from the corresponding author on reasonable request. The datasets analysed for this manuscript are not publicly available since Swiss law does not allow publication of individual data sets. However, further information is available from the corresponding author upon request.
Ethics approval and consent to participate
Ethical approvals were obtained from the local ethical committees of all study sites involved in the study: 1) SIBDCS was approved by the local ethical committees (IRB approval number: EK-1316, granted on February 5, 2007 by the Cantonal Ethics Committee of the Canton Zürich, Switzerland). 2) The Bioethical Committee at the Maria Sklodowska-Curie Memorial Cancer Centre and Institute of Oncology, Warsaw, Poland (ethical licenses 25/2006 and 25/2006/2007). All participants provided written informed consent.
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Quality control for human peripheral blood mononuclear cell enrichment. PBMCs were separated by Ficoll density gradient centrifugation and purified using the EasySep Human Monocyte CD14 Enrichment Kit. Flow cytometry data analysis using antibodies allophycocyanin (APC)-labelled anti-CD14 and Pacific Blue (PB)-labelled anti-CD45 was performed to check cell purity. CD14+ cell purity after enrichment was > 85%. (DOCX 183 kb)
Figure S2. Effect of pH on cAMP formation in (A) THP-1 cells and primary human CD14+ monocytes isolated from (B) healthy subjects. To confirm that the pH values were associated with the activation and inactivation of GPR65/cAMP G-protein mediated signalling a pH dose response curve was generated. THP-1 cells and primary CD14+ human monocytes (WT/CC) were starved at pH 7.6 for 2 h (to silence the receptor) and subsequently subjected to a pH shift for 10 min (pH 6.2 to 7.8 with 0.2 increments). The highest cAMP accumulation was observed at pH 6.4–6.8, whereas low cAMP concentrations were demonstrated at pH 7.6–7.8. (DOCX 674 kb)
Figure S3. Effect of pH on RhoA activation in (A) THP-1 cells and (B) primary human CD14+ monocytes. Description of data: To confirm pH dependent RhoA activity, THP-1 cells and CD14+ monocytes were subjected to different pH (10 min) after a preliminary starvation step (2 h) at non-activating pH (pH 7.6) to silence the receptor. pH 6.6 elicited a significant increase in RhoA activation compared to pH 6.2, 7.4 and 7.6. GPR65/G12/13/RhoA signalling exhibits the highest activity at pH 6.6. (DOCX 65 kb)
Table S1. Allele frequencies of SNP variants and allele association analysis within a population of IBD patients and healthy subjects. Allele frequencies of SNP variants GPR65, rs8005161 and rs3742704 and GALC rs1805078 for allele association analysis within a population of IBD patients and healthy subjects. (DOCX 43 kb)
Table S2. Allele frequencies and biological phenotypes GPR65 SNP rs8005161 for patients from the SIBDC. Allele frequencies and biological phenotypes GPR65 SNP rs8005161 for patients from the Swiss IBD cohort (SIBDC). (DOCX 41 kb)
Figure S4. Formation of cAMP in human CD14+ monocytes upon pH shift from pH 7.6 to pH 6.5. Description of data: (A) Baseline values pH 7.6 and (B) after 10 min at acidic pH (pH 6.5). Human CD14+ cells were obtained from IBD patients carrying either rs8005161 TT, CT or WT/CC genotype, and non-IBD control subjects - all WT/CC genotype. 10 μM of G protein-coupled receptor 65 (GPR65) antagonist was used (C, D). No significant differences between the genotypes were identified. These data are identical to Fig. 1 but presented without normalization. cAMP: cyclic adenosine monophosphate, IBD: inflammatory bowel disease, WT: Wild type. (DOCX 403 kb)