Disrupted Ca2+ homeostasis and immunodeficiency in patients with functional Inositol 1,4,5-trisphosphate receptor subtype 3 defects

Julika Neumann, Erika Van Nieuwenhove, Lara E Terry, Frederik Staels, Taylor R Knebel, Kirsten Welkenhuyzen, Mariah R Baker, Margaux Gerbaux, Mathijs Willemsen, John S. Barber, Irina I Serysheva, Liesbeth De Waele, François Vermeulen, Isabelle Meyts, David I Yule, Geert Bultynck, Rik Schrijvers, Stephanie Humblet-Baron, Adrian Liston
doi: https://doi.org/10.1101/2021.05.29.21257775

Abstract

Calcium signaling is essential for lymphocyte activation, with genetic disruptions resulting in severe immunodeficiency. The inositol 1,4,5-trisphosphate receptor (IP3R), formed from homo- or hetero-tetramers of the IP3R isoforms 1-3, amplifies lymphocyte signaling by releasing Ca2+ from ER stores into the cytosol following antigen-stimulation. While knockout of all 3 IP3R isoforms results in immunodeficiency in mice, the seeming redundancy of subunits was thought to explain the absence of IP3R mutation as a cause of human immunodeficiency. Here, we identify compound heterozygous variants in ITPR3 in two unrelated Caucasian patients presenting with combined immunodeficiency, in one case requiring hematopoietic stem cell transplantation. We observed disrupted Calcium homeostasis in patient-derived fibroblasts and immune cells, with abnormal proliferation and activation responses following B and T cell receptor stimulation. Reconstitution of IP3R knockout cell lines identified the variants as functional hypomorphs with reduced discrimination between homeostatic and induced states, validating a link between genotype and phenotype. These results demonstrate a functional linkage between defective ER Ca2+ channels and immunodeficiency, and identify IP3Rs as diagnostic targets for patients with specific inborn errors of immunity.

Introduction

Spatiotemporally controlled changes in cytosolic Ca2+ concentrations serve as a key signaling mediator in a multitude of physiological processes, from neuron excitation through to cellular apoptosis and lymphocyte activation. In lymphocytes, elevation of the cytosolic Ca2+ concentration ([Ca2+]cyt) is a key event following the engagement of the B cell receptor (BCR) or T cell receptor (TCR). In these cells, Ca2+ signaling occurs via intracellular efflux and extracellular influx in multiple phases. The first step is mediated by the second messenger inositol 1,4,5-trisphosphate (IP3), generated upon activation of phospholipase C (PLC) γ. IP3 binds to and opens IP3 receptors (IP3Rs), thereby releasing Ca2+ from endoplasmic reticulum (ER) stores into the cytosol. While this event only slightly increases [Ca2+]cyt, in the second stage the ER transmembrane protein STIM1 senses lower ER [Ca2+], and via a conformational change leads to the opening of plasmalemmal ORAI1 channels. ORAI1 is a calcium-release activated Ca2+ channel (CRAC) that mediates influx of extracellular Ca2+, a process known as store-operated Ca2+ entry. This sustained increase in [Ca2+]cyt triggers downstream signaling, notably the NF-κB and calcineurin/NFAT pathways1, thereby activating the antigen-stimulated lymphocyte.

Due to the centrality of Ca2+ regulation, a large diversity of genetic disorders are associated with disrupted Ca2+ pathways2. Inadequate Ca2+ signaling has long been associated with severe immunodeficiencies, highlighting the importance of functional Ca2+ signaling in immune cell activation3,4. Genetic defects in extracellular Ca2+ influx are formally associated with primary immunodeficiency, with loss of ORAI1 or STIM1 causing severe combined immunodeficiency (SCID)5,6. By contrast, the primary ER Ca2+ release through IP3R channels, although associated with several non-immunological disorders7,8, has yet to be linked to primary immunodeficiencies. Variants in ITPR1, ITPR2 and ITPR3 have been reported to cause (spino-)cerebellar ataxia9,10, anhidrosis11, and Charcot-Marie-Tooth disease12, respectively. In mice, lymphocyte defects are not observed in individual knockouts, only arising in triple-knockout mice13–15. Based on the 60-80% sequence homology amongst the isoforms16 and the formation of either homo- or hetero-tetramers, subunit redundancy between subunits may occur in lymphocytes. However, whether lymphocytes are susceptible to genetic variants impairing the primary Ca2+ efflux from the ER in humans remained unknown.

Here we identified two patients from unrelated kindreds with immunodeficiency and immune dysregulation that harbor compound heterozygous variants in ITPR3. We show that the three variants identified differentially impact channel function and intracellular Calcium homeostasis, with proliferation and activation defects observed in stimulated lymphocytes. Together these results demonstrate that human lymphocytes are sensitive to genetic defects in the primary ER Ca2+ release system, with niche-filling ITPR3 variants altering IP3R sensitivity. This establishes genetic defects of IP3 receptors as a new class of inborn errors of immunity.

Results

Identification of compound heterozygous ITPR3 variants in two immunodeficient kindreds

To identify the cause of immunodeficiency in two patients without genetic diagnosis, comprehensive clinical, immunological and genetic analysis was performed. Both patients were born to non-consanguineous parents of European descent (Fig. 1A). P1 is male patient in his early adolescence who presented with a combined immunodeficiency with severely reduced numbers of B and T cells, requiring hematopoietic stem cell transplantation (HSCT) in the first decade of his life. He developed an EBV-induced leiomyoma after transplantation and showed abnormal tooth eruption and mineralization as well as thin hair since birth. In his early adolescence, P1 presented with peripheral neuropathy and was diagnosed with Charcot-Marie-Tooth. P2 is a male in his 30’s who presented with recurring immune thrombocytopenia (ITP) requiring splenectomy in his late adolescence. He subsequently suffered from autoimmune hemolytic anemia, susceptibility to infections, and enteropathy. Hypogammaglobulinemia and low switched memory B cells evoked diagnosis with CVID and led to treatment with monthly intravenous immunoglobulins (IVIG). Additional clinical details are available in the Supplementary Material (Supplementary Tables 1-3). None of the parents showed a clinically apparent immunological phenotype and all were self-reported healthy.

Figure 1:
Figure 1: Identification of compound heterozygous variants in ITPR3 in two kindreds with immunodeficiency or immune dysregulation.

(A) Pedigrees of two unrelated families, illustrating the genotype of the patients. Specific genotypes of the parents are not disclosed to protect their privacy. (B) Domain structure of IP3R3 with positions of the variants indicated by red lines. (C) PopViz plot showing the mean allelic frequency (MAF) and the Combined Annotation Dependent Deletion (CADD) score for the identified variants in red. The dashed line indicates the gene-specific mutation significance cutoff (MSC). (D) Cryo-EM structure of homo-tetrameric IP3R3 (6DQJ) depicted as ribbon diagram with different subunits color-coded. Two orthogonal views are shown: side view along the membrane plane (left) and view from the cytosol (right). (E) Illustration of the identified variants mapped onto the 3D structure of IP3R3 (6DQJ). Two opposing subunits are shown in side view. Unresolved sequences are indicated by dashed lines. Red dashed lines span known phosphorylation regions. Putative residues involving Ca2+ binding to the receptor are rendered as pink spheres. (F) Slices through the space-filled model of the IP3R3 TM region perpendicular to the 4-fold axis as seen from the cytosol. Surfaces are color-coded according to electrostatic charges calculated for the model (red: negative; blue: positive).

Whole-exome sequencing was performed on both patients, with filtering for rare damaging variants after excluding the presence of known variants of inborn errors of immunity. We identified compound heterozygote candidate variants in ITPR3, encoding the Inositol-1,4,5-trisphosphate receptor subtype 3 (IP3R3; NM_002224.4). Three distinct variants were identified across the two compound heterozygous patients. P1 harbors a de novo variant, not reported in any public database but recently suggested as a candidate for causing Charcot-Marie-Tooth syndrome12 (c.7570C>T:p.Arg2524Cys). Additionally, P1 inherited a c.5549G>A:p.Arg1850Gln substitution from one of his parents that is also present in P2 and one of his parents (referred to as RQ+/-)and reported with an allelic frequency of 6%. P2 inherited a second private variant from his other parent (referred to as FL+/-), resulting in a c.4882T>C:p.Phe1628Leu substitution. All variants were validated by Sanger sequencing (Supplementary Fig. 1A) and their position within the protein is depicted in Fig. 1B. All protein coding variants were predicted to be damaging by the Sorting Intolerant from Tolerant17 (SIFT) and Deleterious Annotation of Genetic Variants Using Neural Networks18 (DANN) algorithms (Supplementary Table 4). Additionally, all variants had a CADD score above the mutation significance cut-off (MSC 5.7) and are therefore predicted to have a high impact on protein structure and function (Fig. 1C)19. These results demonstrate the identification of two reportedly unrelated patients with immunodeficiency and compound heterozygosity of predicted damaging ITPR3 variants.

To further investigate the potential impact of the identified variants, we mapped them onto the IP3R3 structure solved by cryo-EM20,21. All altered residues are highly conserved across species, with residues F1628 and R2524 also being conserved among all isoforms of IP3Rs (Supplementary Fig. 1B,C). In vivo, the three isoforms IP3R1, IP3R2, and IP3R3 can assemble as homo- or hetero-tetramers to form the functional IP3-regulated Ca2+ channel in the ER membrane. Each IP3R monomer is made up of the cytosolic N-terminal IP3-binding domain followed by the large coupling/regulatory multi-domain region, the pore-forming region comprising six transmembrane (TM) helices and a C-terminal cytosolic tail16,20 (Fig. 1B, D). The F1628L and R1850Q variants are located in one of the regulatory alpha-helical domains, for which some flexible regions remain structurally unresolved. The R1850Q variant is close to the functionally important S1832 phosphorylation site22, while F1628L is in close proximity to several predicted phosphorylation and ubiquitination sites23 (Fig. 1E). The de novo R2524C variant is located in the sixth transmembrane domain helix (S6). The S6 helices from each subunit shape the channel’s ion conduction pathway and form the constriction gate that controls ion translocation across the ER membrane. R2524 is located in the cytosolic extension of the S6 helix just beyond the gate21 and its positive charge is likely neutralized through the coordination with neighboring Aspartate residues24. According to computational modeling, replacement of the positively charged Arginine by the uncharged Cysteine abolishes the charge coordination in the TM helix S6, thereby changing the pore charge in this region (Fig. 1F and Supplementary Fig. 1D). Due to the strong conservation of all identified variants and their location in or close to predicted functionally relevant sites we hypothesize that these variants alter channel function.

Functional impact of ITPR3 variants

To assess the functional impact of ITRP3 variants, we investigated expression levels and Calcium flux in patient cells. First, we analyzed the influence of the variants on mRNA and protein expression levels of all IP3R subtypes using isoform-specific primer pairs and antibodies selectively recognizing one of the isoforms or all of them simultaneously (panIP3R). The antibodies used were raised against immunogens N-terminal of the identified mutations, preserving recognition capacity of mutant IP3R3 proteins. mRNA expression was normal in fibroblasts (Fig. 2A), however ITPR3 was reduced in both patients and the parents of P2 in peripheral blood mononuclear cells (PBMCs, Fig. 2B). Immunoblot analysis yielded strong immunoreactive bands for IP3R1 and IP3R3, consistent with their higher prevalence in fibroblasts and immune cells compared to low levels of IP3R2 (Supplementary Fig. 2A,B). We detected reduced levels of IP3R3 in P1 fibroblasts (PBMCs were not available, due to HSCT) and P2 PBMCs (Fig. 2C,D). Parents of P2 showed intermediate expression levels (Fig. 2D). While expression of IP3R subtypes, and the impact of variants, varies across investigated cell types, these results suggest a reduction of ITPR3 in both patients at the mRNA and protein level in PBMCs.

Figure 2:
Figure 2: IP3R3 expression is reduced by variants in a cell type-dependent and variant-specific manner.

mRNA expression of all three receptor subtypes was measured in primary fibroblasts and PBMCs with isoform specific primers. Dots represent the mean of duplicates for independent runs normalized to the mean values of the housekeeping genes GAPDH and eIF4E. (A) Expression in fibroblasts (three independent runs for HD n=3, P1 n=2, P2 n=1 biological replicates), and (B) PBMCs (three independent runs for HD n=7, FL+/- n=2, RQ+/- n=2, P2 n=3 biological replicates). (C) Western Blot analysis of IP3R2 and IP3R3 isoform protein expression was performed using subtype-specific antibodies. Representative blots are shown (left) with the quantification of protein expression normalized to the housekeeping gene Vinculin (right). Quantification was performed in fibroblasts (HD n=2, P1 n=2, P2 n=2) and (D) PBMCs (HD n=4, FL+/- n=2, RQ+/- n=2, P2 n=2). Values are represented as mean + SEM. One-way ANOVA with multiple comparisons.

As IP3Rs are crucial for Calcium dynamics and mediating IP3-induced Ca2+ flux, we sought to investigate Ca2+ flux in primary patient cells to test IP3R functionality. As pre-transplantation PBMCs were limited from P1, we used fibroblasts cultured from skin biopsies and monitored the change in [Ca2+]cyt as a response to various stimuli using a ratiometric Ca2+ indicator (Fura2). We first assessed total intracellular Ca2+ response in fibroblasts stimulated with the Ca2+ ionophore Ionomycin in the presence of extracellular EGTA, a Ca2+-chelating agent. Both the integrated response as well as the peak amplitude were reduced in P1 (Fig. 3A), indicating impaired intracellular Calcium homeostasis. Likewise, P1 demonstrated reduced [Ca2+]cyt following treatment with Thapsigargin, an inhibitor of the sarco-/endoplasmic reticulum Ca2+-ATPase (Fig. 3B), indicating a reduced steady-state ER Ca2+ ER store content. This could be due to an increased leakiness of Ca2+ through the ER membrane, thereby resulting in a lower steady-state ER Ca2+ concentration and thus a reduced resting concentration gradient between the ER lumen and the cytosol. Finally, we monitored IP3R-mediated Ca2+ release in fibroblasts in response to the physiological GPCR agonist Bradykinin and observed a decreased response for P1 (Fig. 3C). These results indicate that overall Ca2+ homeostasis is altered in P1, with impaired IP3R-mediated Ca2+ release. In each of these assays P2 showed a comparable response to two different healthy controls.

Figure 3:
Figure 3: Variants in ITPR3 impede Calcium homeostasis and Calcium signaling.

Adherent fibroblasts were loaded with the ratiometric Ca2+ indicator Fura2/AM, and Calcium flux was monitored over time. The first dashed line indicates addition of EGTA for acquisition of background, the second dashed line indicates addition of (A) 10 µM Ionomycin, (B) 10 µM Thapsigargin or (C) 10 µM Bradykinin. Area under the curve (AUC) and maximal amplitude of the response were calculated respective to background values. HD n=12, P1 n=9, P2 n=4. (D) PBMCs were stained for CD4, loaded with the ratiometric Ca2+ indicator FuraRed and monitored for changes in cytosolic Ca2+ concentrations by flow cytometry. Mean of all experiments performed in duplicates or triplicates is depicted. Cells were treated with 10 µM Ionomycin or (E) 10 µM Thapsigargin. First dashed lines indicate addition of EGTA and second dashed lines indicate addition of the stimulant. HD n=7, FL+/- n=3, RQ+/- n=3, P2 n=5. Traces represent mean values of all experiments and response quantification is shown as mean + SEM. One-way ANOVA with multiple comparisons.

As PBMCs have a different Calcium signaling sensitivity than fibroblasts, we sought to test CD4+ T cells from P2 in stimulation assays (P1 was not available, due to successful HSCT). By labeling CD4+ T cells from different donors with different fluorochromes before combining the samples, we were able to reduce sample-to-sample variation. Cells were loaded with the ratiometric Ca2+ indicator FuraRed, and stimulated with Ionomycin or Thapsigargin. For both stimuli the increase in [Ca2+]cyt was reduced in individuals of kindred 2 when compared to healthy donors (Fig. 3D,E). These results suggest defects in IP3R-mediated Ca2+ release in both patients, with P1 potentially the more severe (also being detected in fibroblasts), although direct comparison in PBMCs was not possible due to HSCT. Defects were also observed in the heterozygous parents of P2.

To formally test the impact of the ITPR3 variants on Calcium signaling in a model allowing direct comparison in the absence of compensatory mechanisms, we used HEK cells that have been engineered to lack endogenous expression of all three IP3R subtypes25 (HEK-3KO) and reintroduced either wild-type (exogenous) or mutant IP3R3. Double-knockout of IP3R1 and IP3R2, resulting in cell lines only expressing endogenous IP3R3, were used as a control. We generated multiple monoclonal cell lines with transgenic expression of the mutant allele, confirmed localization to the ER by immunohistochemistry, and used Western Blot to benchmark protein expression levels against endogenous expression (Supplementary Fig. 3A,B). Using these knockout cell lines reconstituted with mutant or wild-type alleles, we investigated the Calcium flux induced by different concentrations of the muscarinic GPCR agonist Carbachol (CCH) and analyzed the magnitude of the response (Supplementary Fig. 3C-E). The inherited private F1628L variant of P2 and the inherited R1850Q variant shared by P1 and P2 resulted in reduced or absent responses under physiological protein levels, requiring more than a 40-fold over-expression to achieve similar sensitivity to the endogenous protein (Fig. 4 A,B, Supplementary Fig. 3F,G). The de novo R2524C variant of P1 showed even more severe defects, with absent Ca2+ mobilization across all clones (Fig. 4C, Supplementary Fig. 3H). Together, these results formally demonstrate IP3R-mediated Ca2+ release defects in the F1628L, R1850Q and R2524C mutants in an overexpression model, with the most striking defect being found in the mutation borne by P1, who presented with the more severe clinical manifestations of immunodeficiency.

Figure 4:
Figure 4: IP3-induced Ca2+ flux is reduced through F1628L, R1850Q, and R2524C mutant IP3R3.

HEK-3KO cell lines, not expressing any IP3R subtype, were stably transduced with the wild-type (exogenous) or mutated versions of the IP3R3 receptor and stimulated with different concentrations of the GPCR agonist Carbachol in Ca2+ imaging buffer. Cells only expressing endogenous IP3R3 (endogenous) and HEK-3KO cells were used as controls. Dose-response curves of the AUC of representative cell lines stably expressing IP3R3 with the (A) F1628L (n=3), (B) R1850Q (n=3) or (C) R2524C (n=3) variant at different concentrations of Carbachol (CCh) are shown on the left. IP3R3 protein expression normalized to the endogenous expression was analyzed by WB and is shown on the right. Values are represented as mean ± SEM.

Defective immunological responses in ITPR3 patients

Following our observations that Calcium homeostasis is disrupted in immune cells of our patients we investigated the downstream consequences on immune cell function. P1 was, of necessity, excluded due to successful transplantation of the hematopoietic compartment. First, we investigated the impact of variants in ITPR3 on early downstream activation via phosphorylation of Erk and PLCγ1 in response to stimulation. Although PLCγ1 signals upstream of IP3Rs, it is positively regulated by Ca2+ and thereby amplified following IP3R activation26. Following TCR stimulation, the percentage of pPLCγ1+ positive cells was reduced in members of kindred 2 at different time points (Fig. 5A,B). Cells that responded showed comparable phosphorylation of PLCγ1 (Fig. 5C), indicating that the defect was in triggering the pathway. Erk activation, while highly variable among experiments, suggested defects only in P2 but not the heterozygous parents (Fig. 5D-F). Stimulation of the BCR gave relatively normal pErk positivity in B cells (Fig. 5G-I). These results suggest a functional defect of IP3R3 in P2 with impaired lymphocyte receptor signaling.

Figure 5:
Figure 5: Antigen receptor downstream phosphorylation events are reduced in P2.

PBMCs were stimulated with anti-CD3 and anti-CD28 or immunoglobulins for 5, 15 or 30 min. (A) Representative histograms illustrating PLCγ1 phosphorylation following TCR stimulation after different stimulation intervals. (B) Percentage of T cells with phosphorylated PLCγ1 in (un)stimulated cells. (C) Fold-change of median fluorescence intensity in pPLCγ1+ cells compared to unstimulated HD cells. HD n=10, FL+/- n=2, RQ+/- n=2, P2 n=6. (D) Representative histograms illustrating Erk phosphorylation following TCR stimulation after different stimulation intervals. (E) Percentage of T cells with phosphorylated Erk in (un)stimulated cells. (F) Fold change of median fluorescence intensity in pErk+ cells compared to unstimulated HD cells. HD n=7, FL+/- n=1, RQ+/- n=1, P2 n=3. (G) Representative histograms illustrating Erk phosphorylation following BCR stimulation after different stimulation intervals. (H) Percentage of B cells with phosphorylated Erk in (un)stimulated cells. (I) Fold change of median fluorescence intensity in pErk+ cells compared to unstimulated HD cells. HD n=7, FL+/- n=2, RQ+/- n=2, P2 n=5. Values are represented as mean + SEM. One-way ANOVA with multiple comparisons.

Finally, we assessed the downstream consequences of lymphocyte activation by analyzing the proliferative potential of CD4+ and CD8+ T cells in response to TCR engagement. Total cell numbers after stimulation were reduced in all individuals from kindred 2 when compared to healthy donors (Supplementary Fig. 4A,B). TCR-induced cell division was substantially reduced in CD4+ and CD8+ T cells of P2, with intermediate phenotypes exhibited in the heterozygous parents (Fig. 6A,B, Supplementary Fig. 4C,D). Interestingly, bystander activation of CD19+ B cells27 was also sharply impeded in P2 cells (Fig. 6C, Supplementary Fig. 4E). The reduced proliferative capacity of P2 cells was even more pronounced when stimulating PBMCs with soluble anti-CD3 and soluble anti-CD28, mimicking suboptimal stimulation conditions (Fig. 6D-F, Supplementary Fig. 4E-G). Taken together, our results demonstrate impaired immunological responses to TCR stimulation in P2, with weak effects observed in the heterozygous patients, consistent with the defects in Calcium flux and phosphorylation of Erk and PLCγ1.

Figure 6:
Figure 6: Altered cellular responses with ITPR3 variants.

PBMCs were labelled with the CFSE cell proliferation dye and stimulated with anti-CD3 and anti-CD28 for up to 4 days. The mean number of proliferations of the whole cell population, including undivided cells, was calculated for CD4+, CD8+ and CD19+ subsets (Division Index). Stimulation was performed in plates coated with 10 µg/mL anti-CD3 and 5 µg/mL soluble anti-CD28. Exemplary histograms are shown on top and the proliferation index for HD n=9, FL+/- n=3, RQ+/- n=3, P2 n=5 after 3 and 4 days shown on the bottom for CD4+ T cells, (B) CD8+ T cells and (C) B cells, reading out on bystander activation rather than direct stimulation. (D) Stimulation was performed with 5 µg/mL soluble anti-CD3 and 5 µg/mL soluble anti-CD28. The Division Index for HD n=2, FL+/- n=1, RQ+/- n=1, P2 n=3 is shown for CD4+ T cells, (E) CD8+ T cells and (F) B cells. Values are represented as mean + SEM. One-way ANOVA with multiple comparisons.

Discussion

In this work we present two patients from unrelated families with compound heterozygous variants in ITPR3 that present with immunodeficiency and immune dysregulation. We show that both patients exhibit altered Calcium homeostasis and responses to in vitro stimulation, using cell proliferation and phosphorylation of major signaling proteins as readouts. Our results suggest that the R2524C variant results in the most severe impairment of channel function, a cellular phenotype in line with the more severe combined immunodeficiency observed clinically in P1. For P2, harboring the F1628L and R1850Q variants, we demonstrate defects in all in vitro assays, with a profound reduction of proliferative capacity following TCR stimulation. The partially reduced sensitivity of channel activation correlates with the weaker clinical phenotype. Together with the intermediate cellular responses in the clinically healthy parents of P2, these results suggest a spectrum of escalating clinical severity with progressive decline in channel function.

For the F1628L and R1850Q variants we propose a recessive mode of inheritance based on our observations. Both variants are located in the regulatory domain of the receptor that is known for its function as signaling hub, interacting with a plethora of proteins28,29. We expect that each variant impairs a specific aspect of regulation, which cannot entirely be compensated for during lymphocyte activation, thus leading to signaling defects and impaired Calcium homeostasis in vitro. As heterozygous individuals did not show a clinical phenotype we hypothesize that lymphocytes can partially compensate for individual losses in the heterozygous state. In this model, the heterozygous loss of two regulatory partners for IP3R3 is sufficient to induce combined immunodeficiency in P2. An alternative explanation for the additive effect of different variants is the reported cooperative nature of channel activation, which requires all subunits of the hetero-tetramer to bind IP3 and undergo the necessary changes for channel gating30. The compound heterozygosity increases the chance of an IP3R containing at least one defective subunit. The de novo generation of the R2524C variant in an individual inheriting the R1850Q allele does not allow us to discriminate between recessive and autosomal dominant effects of this allele. However, this allele was also observed with dominant inheritance for Charcot-Marie-Tooth in Rönkkö et al.12. The absence of a reported immunological phenotype in the previously described patient suggests that sole inheritance of R2524C creates sufficient impairment to lead to Charcot-Marie-Tooth, while preserving sufficient function for normal lymphocyte activation unless coupled with an additional defective allele, as in P1. This split inheritance pattern is consistent with data demonstrating cell type- and context-dependent partial redundancies of the IP3R isoforms observed in mice and humans 7,31,32. Alternatively, the unique genetic background in the Ashkenazi Jewish patient may alter the penetrance in lymphocytes, or subclinical effects may be present which have yet to develop into a clinical phenotype.

Despite advances in the accuracy of electron microscopy techniques revealing conformational changes and interactions in specific activation states of IP3Rs 20,21,24, the effect of post-translational modifications and protein interactions remains poorly studied and new interaction partners and signaling mechanisms are continuously being discovered. Genetic variants leading to human disease provide a unique way for better understanding protein regulation and function. Our results demonstrate that the F1628L and R1850Q variants each result in reduced sensitivity to channel activation while still maintaining its gating ability. Systematic investigation of mutated receptor variants could reveal the protein interactions and/or conformational changes that are hampered and provide better insight into channel regulation. The de novo R2524C variant in P1, however, completely abolishes gating of a homo-tetrameric receptor. The respective conserved residue in the related Ryanodine receptor (RyR) family has indeed been shown to stabilize both the closed and open state of the channel by alternating its interaction with two other highly conserved residues33. In line with the prediction of a drastic change in the pore charge introduced by this variant, its counterpart in the rat IP3R1 protein was suggested to have a crucial function for signal transmission as mutation into different amino acids drastically reduced channel function24,34.

With this study we add variants in ITPR3 as an underlying cause for disturbed Calcium signaling resulting in immunodeficiency. Although we observed a variable clinical presentation and effect on protein function depending on the location of the variants, the phenotype of the described patients resembles that of immune defects described in patients with ORAI1 and STIM1 deficiency. These patients show severely reduced or absent store-operated Ca2+ entry in vitro and present with normal lymphocyte counts but impaired T cell proliferation following in vitro stimulation 4,6. Comparable to the patients we report, they suffer from recurrent severe infections and can exhibit signs of lymphoproliferation and autoimmune disease, although the disease course is often more severe, requiring HSCT very early in life to avoid a lethal disease course 35,36. However, phenotypic presentation also varies greatly in STIM1 deficiency, with several STIM1 deficient patients described with late-onset or absent clinical immunodeficiency 37,38, and phenotypic expansion is commonly observed with increasing case recognition in primary immunodeficiency. In P1 we also observed non-immunological disease manifestations, such as mineralization defects of teeth with abnormalities of the hair and Charcot-Marie-Tooth disease, which is partly overlapping with dental enamel formation defects and myopathy - although congenital and non-progressive - in ORAI1 and STIM1 patients 6,37,39,40. We therefore propose that ITPR3 should be considered as a candidate gene in patients with B and T cell defects in association with impaired Calcium signaling function, adding variants in ITPR3 as a cause for inborn errors of immunity.

Material and Methods

Written informed consent to conduct and publish this study was obtained from all participants and the ethics committee of UZ/KU Leuven approved the study (S52653, S58466). For details on the clinical history of the patients please contact the clinical corresponding authors.

Whole exome sequencing and analysis

Whole exome sequencing for kindred 1 was performed using genomic DNA isolated from whole blood using the QIAmp DNA Blood Midi kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. The sequencing library was prepared using ligation-mediated PCR and hybridizing to the SureSelect Biotinylated RNA library (Agilent Technologies, Santa Clara, CA) for enrichment and performed by BGI (Shenzhen, China) before paired-end sequencing on a HiSeq2000 (Illumina, San Diego, CA) platform. The generated FastQ files were mapped to the human genome version 19 (hg19) using the Burrows Wheeler Aligner (BWA, v0.7.5a). Duplicates were then removed using Picard MarkDuplicates and realignment performed according to the Genome Analysis Toolkit 3.1 best practices guidelines. Variants were called using the HaplotypeCaller and filtered for coding, non-synonymous variants with a high CADD score and a mean allelic frequency of <0.005 in the gnomAD database.

For kindred 2, whole exome sequencing was performed by Macrogen (Seoul, Korea). Exome capture using SureSelect Human All Exon V7 (Agilent) preceded paired-end sequencing on a NovaSeq6000 (Illumina) platform. A computational pipeline was developed using bcbio-nextgen as backend (https://github.com/bcbio/bcbio-nextgenv1.1.5-b), to process the read data and perform tasks such as quality control (QC), variant discovery, annotation and filtering. Briefly, the sequencing reads in FASTQ format were aligned to the human reference genome (GRCh37) using BWA (v0.7.17). The resulting BAM files were further processed to remove duplicate reads using biobambam (v2.0.87). Resulting variants were annotated using snpEff (v4.3.1t) for function prediction and vcfanno (v0.3.1) for adding public databases like dbsnp (v151), dbNSFP (v3.5a), ExAC, gnomAD (v r2.1), 1000 genomes (v3) and clinvar (v2019-05-13). Likely disease-causing mutations were selected and prioritized based on quality score, allele frequency, functional impact, and probable inheritance model. Sanger confirmation for individuals from both kindreds was performed by amplification of variant-specific gene products by PCR using KOD polymerase (Sigma-Aldrich, St. Louis, MO, US) and primers listed in Supplementary Table 5. After gel purification with the NucleoSpin Gel and PCR Clean-up (Macherey-Nagel, Düren, Germany), according to the manufacturer’s instructions, sequencing was performed by Eurofins (Ebersberg, Germany). Variant effect prediction on protein structure

The cryo-EM structure of IP3R3 with the PDB accession number 6DQJ was used to analyze the identified variants. Computational mutagenesis was performed in COOT41 and subsequent molecular refinements of R2524C were performed in PHENIX42. Calculations of electrostatic potential and all molecular visualizations were generated using UCSF Chimera43.

RNA isolation and quantification

Total RNA was isolated from PBMCs and primary fibroblast cell lines generated from skin biopsies using TRIzol reagentList of Supplementary Materials (Ambion, Thermo Fisher, Waltham, MA), according to the manufacturer’s protocol. Complementary DNA was synthesized using the GoScript™ Reverse Transcription System (Promega, Madison, WI). Quantitative PCR was performed on a StepOnePlus real-time PCR system (ABI, Thermo Fisher) with Fast SYBR Green Master Mix (Applied Biosystems, Foster City, CA) supplemented with gene-specific primers (see Supplementary Table 6). Experiments were performed in duplicate and repeated thrice. Gene expression was normalised to the mean of two housekeeping genes before normalising to the mean expression of all healthy controls across experiments.

Western Blot

Primary fibroblasts and PBMCs were solubilized in lysis buffer containing 20 mM Tris-HCl pH 7.5, 150 mM NaCl, 1% Triton-X, 1.5 mM MgCl2, 0.5 mM DTT and protease inhibitor. Protein concentrations were determined by BCA and 20 µg of protein per sample separated on a 3-8% Tris-Acetate gel. Proteins were blotted on a PVDF membrane and incubated with protein-specific primary antibodies after blocking. Primary antibodies used at a 1/1000 dilution were in-house made rabbit anti-IP3R1 Rbt03 antibody44, rabbit IP3R2 (Abicode C2), mouse anti-IP3R3 (IP3R3, BD Biosciences) and in-house made rabbit pan-IP3R Rbt475 antibody raised against a peptide corresponding to amino acids 127-141 of human IP3R145,46. Vinculin (Sigma) was used as a loading control at a 1/2000 dilution. Secondary antibodies coupled to HRP were goat anti-mouse/rabbit IgG (Thermo Scientific). Membranes were developed using Pierce™ ECL (Thermo Fisher) in a Chemidoc MP Imaging system (BioRad, Hercules, CA). Quantification of protein expression was performed with ImageJ47,48. The fold-change in protein expression was normalized to the mean of all healthy controls within each independent experiment in Excel.

HEK-3KO cells, HEK293 cells expressing only endogenous IP3R3, and HEK-3KO cells stably expressing IP3R constructs were solubilized in RIPA lysis buffer supplemented with Halt™ protease inhibitor cocktail (Thermo Fisher). Protein concentrations were determined using Dc protein assay kit (BioRad) and 7.5 µg of protein per sample prepared in SDS loading buffer separated using 7.5% SDS-PAGE. Subsequently, proteins were transferred overnight to a nitrocellulose membrane (Pall Corporation, New York, NY) which was probed using the indicated primary antibodies (mouse monoclonal antibody recognizing residues 22-230 of human IP3R3 (BD Transduction Laboratories, San Jose, CA) and rabbit monoclonal antibody recognizing Calnexin (Cell Signaling, Danvers, MA)) and corresponding secondary antibodies. Membranes were imaged with an Odyssey infrared imaging system (LICOR Biosciences, Lincoln, NE). Resulting scans were exported to LICOR Image Studio where band intensity was calculated. In Excel, hR3 band intensity values were normalized to that of the calnexin loading control. Subsequently, these values were further normalized to that of the corresponding endogenous hIP3R3. Averages and statistical analysis were performed in GraphPad Prism.

Plasmid cloning

All plasmids used in this study were based on the primary sequence of human IP3R3 in a pcDNA3.1 backbone. Mutagenesis to introduce base changes identified in the reported patients was performed with the primers listed in Supplementary Table 7 that were synthesized by Integrated DNA Technologies. Suitable sequences were amplified by PCR and subcloned into pJet1.2/blunt using the CloneJET PCR Cloning kit (Thermo Fisher) according to the manufacturer’s instructions. Mutagenesis PCR was performed with outward-facing complementary primers introducing the desired base change and amplifying the full plasmid. DpnI digest in Tango Buffer (both Thermo Scientific) was performed overnight and the resulting circularized plasmids used to transform chemically competent DH5α bacteria. Following Ampicillin selection, plasmids were purified from liquid cultures of single colonies using the NucleoSpin Plasmid EasyPure kit (Macherey-Nagel) and base exchange verified by Sanger sequencing (Eurofins). Restriction-enzyme based cloning was then performed to introduce the mutated sequences into the original pcDNA3.1 construct and sequences verified again.

Generation of HEK-3KO cell lines stably expressing IP3R3 mutants

HEK-3KO cells are HEK293 cells engineered though CRISPR/Cas9 for the deletion of the three endogenous IP3R isoforms25, while cells expressing endogenous IP3R3 are HEK293 cells modified by CRISPR/Cas9 to only express endogenous IP3R3. The HEK-3KO cells, cells expressing endogenous IP3R3, and HEK-3KO cells stably expressing exogenous IP3R constructs were grown at 37 °C with 5 % CO2 in DMEM media supplemented with 10% fetal bovine serum, 100 U/ml penicillin, 100 µg/ml streptomycin. Cell lines were maintained by subculturing as necessary and any selection required was carried out using G418 (VWR, Radnor, PA). HEK-3KO transfection was also performed similarly to previously described 8,25. In brief, 1×106 cells were pelleted, washed once with PBS, and resuspended in a homemade transfection reagent containing ATP-disodium salt, MgCl2, KH2PO4, NaHCO3, glucose, and pH to 7.4. 2-5 µg of DNA were mixed with the resuspended cells and electroporated using the Amaxa cell nucleofector program Q-001. Cells recovered in fresh DMEM media, supplemented as listed above, for 48 h prior to passage into new 10 cm dishes containing DMEM supplemented with 1.5-2 mg/ml G418. About 7 days after start of selection individual colonies were picked by hand and transferred to 24-well plates containing fresh DMEM with G418. 10-14 days after transfection, wells that exhibited growth were expanded and screened using western blot to verify stable expression of the desired protein.

Immunocytochemistry and Confocal Microscopy

HEK-3KO cells stably expressing exogenous IP3R constructs were plated on poly-d-lysine coated coverslips. When roughly 50 % confluent, cells were fixed using 4 % PFA at room temperature for 15 minutes. Subsequently, coverslips were washed with PBS and cells were blocked in 10 % Bovine Serum Albumin (BSA) in PBS-T (PBS with 0.2 % Tween20) for 1.5 h. Following blocking, cells were incubated in anti-IP3R3 primary antibody in BSA overnight at 4 °C. The following day, primary antibody was removed, and coverslips were washed 3 times with PBS-T for 10 minutes with gentle rocking. Subsequently, secondary antibody conjugated to AlexaFluor488 was incubated for 1 h in BSA at RT with gentle rocking. After incubation, coverslips were washed with PBS-T and mounted on slides. After allowing slides to dry, coverslips were sealed onto slides and imaged using confocal microscopy.

Fluorescence measurement of Calcium fluxes

Fibroblasts were seeded in 96-well plates (Greiner, Kremsmünster, Austria) and assessed for [Ca2+]i concentrations using Fura-2/AM (Eurogentec, Seraing, Belgium). Briefly, cells were loaded with the ratiometric fluorescent Ca2+ indicator Fura2/AM (1 µM) at RT for 30 min in a modified Krebs solution (containing 150 mM NaCl, 5.9 mM KCl, 1.2 mM MgCl2, 11.6 mM HEPES (pH 7.3) and 1.5 mM CaCl2). Cells were rested for 30 min at RT in the absence of Fura-2 AM to allow complete dye de-esterification before proceeding to analysis on a FlexStation 3 microplate reader (Molecular Devices, Sunnyvale, CA, USA). The indicator was alternately excited at 340 and 380 nm and emission of fluorescence at 510 nm recorded. EGTA was added after 30 sec in all conditions at a final concentration of 3 mM, to chelate Ca2+ ions present in the buffer, and fluorescence recorded for 60 sec prior to stimulation. Responses to stimuli, prepared in Ca2+-free modified Krebs solution containing 3 mM EGTA, were acquired for 6 min. Ionomycin and the irreversible SERCA-inhibitor Thapsigargin were added at a final concentration of 10 µM and Bradykinin at a final concentration of 50 nM. All traces are shown as the ratio of both emission wavelengths F340/F380.

PBMCs from different donors were first stained with different CD4 antibodies for multiplexing (all RPA-T4, eBioscience, BD Biosciences, Invitrogen). After mixing together they were loaded with the ratiometric fluorescent Ca2+ indicator dye FuraRed AM (Thermo Fisher) as described above. Cells were analyzed on a BD FACSCanto (BD Biosciences, Franklin Lakes, NJ, US). The ratio of emission at 450/50 nm after excitation at 405 nm and emission at 670 nm after excitation at 488 nm was exported using the Kinetics platform of the FlowJo™ software (v10.7.1, Ashland, OR). Raw data from both approaches were smoothened using a 2nd order running average of 5 in GraphPad Prism (v9.0.0, San Diego, CA). A baseline value was calculated for each measurement as mean fluorescence after addition of EGTA and before addition of the stimulus. This baseline was used for calculation of the area under the curve (AUC) and the peak amplitude using GraphPad Prism with a minimum peak height of 10 % of the distance from minimum to maximum Y (v9.0.0). HEK-3KO cells and HEK-3KO cells stably expressing exogenous IP3R constructs were analyzed on a FlexStation 3 microplate reader (Molecular Devices) after loading with 4 µM Fura-2/AM (TEFLabs, Austin, TX) in complete DMEM. After 1 h, cells were harvested, washed, and resuspended in Ca2+ imaging buffer before dispensation into a black-walled flat bottom 96-well plate (Greiner). The plate was centrifuged (200xg for 2 min) and placed at 37 °C for 30 min prior to commencing the assay. Excitation was performed as described above and Ca2+ imaging buffer or varying concentrations of CCh were added to induce IP3R-mediated Ca2+ release. Readings were taken every 4 sec for a total of 200 sec and data was exported from SoftMax® Pro Microplate Data Acquisition and Analysis software to Excel where the F340/F380 ratio was calculated. Ratios were normalized to the average of the first five ratio values and the AUC and peak amplitude calculated in GraphPad Prism as described above. Data was averaged from at least 3 individual plates and logistic curve fits were calculated in GraphPad Prism.

Phosphorylation assay

Frozen PBMCs were thawed, plated and rested in complete RPMI at 37 °C prior to stimulation with 5 µg/mL aCD3 (UCHT1) and 5 µg/mL aCD28 (CD28.2, both eBioscience) or anti-human IgA/IgG/IgM (Jackson ImmunoResearch, West Grove, PA), for durations as indicated. Cells were fixed in 2 % paraformaldehyde and permeabilized with 100 % methanol, prior to staining for CD14 (M5E5, BioLegend), CD4 (RPA-T4, eBioscience), CD8 (SK1, eBioscience), CD19 (HIB19, eBioscience), phospho-Erk (Thr202/Tyr204, MILAN8R, eBioscience) and phospho-PLCγ1 (Tyr783, both Biolegend). Cells were acquired on a BD FACSCanto (BD Biosciences). After exclusion of doublets and monocytes via CD14, T cells were gated as being CD19- and B cells as CD19+. The fold change of the median fluorescence intensity was calculated based on the average of all healthy controls within one repeat of the experiment.

Proliferation assay

Frozen PBMCs were thawed in RPMI 1640 (Gibco™, ThermoFisher) supplemented with 10 % FBS (Tico Europe, Amstelveen, Netherlands) and 100 U/mL Penicillin/Streptomycin (Gibco™, ThermoFisher) prior to labelling with 1 µM CFSE (Invitrogen™, Carlsbad, CA) at 37 °C with regular shaking. The reaction was quenched with ice-cold FBS and cells plated after washing in complete medium containing 10 % FBS, 100 U/mL Penicillin, 25 mM HEPES and non-essential amino-acid. Cell stimulation was either performed in wells coated with 10 µg/mL anti-human CD3 (UCHT1) monoclonal antibodies and 5 µg/mL anti-human CD28 monoclonal antibodies (CD28.2, both eBioscience™, San Diego, CA) supplemented in the medium, or with both monoclonal antibodies supplemented in the medium at 5 µg/mL. Cells were stained with CD4 (RPA-T4, Invitrogen), CD8 (OKT-8, Invitrogen) and CD19 (HIB19, eBioscience) and analyzed by flow cytometry on a BD FACSCanto (BD Biosciences, Franklin Lakes, NJ, US) on day 3 and 4 after stimulation.

Statistics

Ordinary One-way ANOVA with Dunnett’s multiple comparisons test was performed using GraphPad Prism 9.0.0 (San Diego, CA).

Data Availability

All data will be made available upon reasonable request.

List of Supplementary Materials

  • ‐ Supplementary Table 1: Demographic, clinical, laboratory, technical and treatment data for investigated patients.

  • ‐ Supplementary Table 2: Laboratory values of P1 at different relevant time points.

  • ‐ Supplementary Table 3: Laboratory values of P2 at different relevant time points referred to in the manuscript.

  • ‐ Supplementary Table 4: Protein-coding variants in ITPR3 are predicted to be damaging.

  • ‐ Supplementary Table 5: Primers for Sanger sequencing of identified variants.

  • ‐ Supplementary Table 6: Primers used for subunit-specific qPCR analysis.

  • ‐ Supplementary Table 7: Primers used for cloning.

  • ‐ Supplementary Figure 1: Identified variants in ITPR3 are highly conserved.

  • ‐ Supplementary Figure 2: Western Blot analysis of IP3R2 and total IP3R expression is not reliably quantifiable.

  • ‐ Supplementary Figure 3: Reconstitution of HEK-3KO cells with mutated IP3R3 results in stable protein expression and reveals impaired channel function.

  • ‐ Supplementary Figure 4: Defective TCR-induced T cell proliferation in individuals from kindred 2.

Disclosure of Conflicts of Interest

IM holds the CSL Chair in Primary Immunodeficiencies, paid to KU Leuven. The other authors have declared that no conflict of interest exists.

Contributions

Patient selection, immunophenotypical characterization and ethics: IM, RS; Target identification: EVN, FS; Conceptualization: SHB, AL; Formal analysis: JN, EVN, LET; Investigation: JN, EVN, LET, FS, TK, KW, MG, MW, JB; Resources: IS, LDW, RS, IM, DIY, GB; Writing – original draft: JN, EVN, SHB, AL; Writing – review and editing: IM, RS, SHB, AL; Visualization: JN, LET; Supervision: SHB, AL; Funding acquisition: RS, SHB, AL

Funding and Acknowledgements

This work was supported by the VIB Grand Challenges Program, the KU Leuven C1 program, the European Union’s Horizon 2020 research and innovation programme under grant agreement No 779295 (to AL), and the Biotechnology and Biological Sciences Research Council (BBSRC) through Institute Strategic Program Grant funding BBS/E/B/000C0427 and BBS/E/B/000C0428 and the KU Leuven BOFZAP start-up grant (to SHB). GB is a recipient of a Hercules Foundation grant (AKUL/19/34) for the financing of equipment. GB, IS and DIY are partners of the FWO Scientific Research Network (CaSign W0.019.17N). IM and RS are FWO senior clinical investigator fellows. IM and RS are members of the European Reference Network for Rare Immunodeficiency, Autoinflammatory and Autoimmune Diseases (project ID No. 739543). JN, EVN and FS are FWO fellows. The authors acknowledge the important contributions of Tomas Luyten, Pier-Andrée Penttila and the KU Leuven FACS Core, and Stefanie Sente and Giorgia Bucciol for clinical trial support.

References

  1. 1.↵
    Feske S. Calcium signalling in lymphocyte activation and disease. Nat Rev Immunol. 2007;7(9):690–702. doi:10.1038/nri2152
    OpenUrlCrossRefPubMed
  2. 2.↵
    Berridge MJ. The Inositol Trisphosphate/Calcium Signaling Pathway in Health and Disease. Physiol Rev. 2016;96(4):1261–1296. doi:10.1152/physrev.00006.2016
    OpenUrlCrossRefPubMed
  3. 3.↵
    Partiseti M, Le Deist F, Hivroz C, Fischer A, Korn H, Choquet D. The calcium current activated by T cell receptor and store depletion in human lymphocytes is absent in a primary immunodeficiency. J Biol Chem. 1994;269(51):32327–32335. http://www.ncbi.nlm.nih.gov/pubmed/7798233
    OpenUrlAbstract/FREE Full Text
  4. 4.↵
    Le Deist F, Hivroz C, Partiseti M, et al. A primary T-cell immunodeficiency associated with defective transmembrane calcium influx. Blood. 1995;85(4):1053–1062. doi:10.1182/blood.V85.4.1053.bloodjournal8541053
    OpenUrlAbstract/FREE Full Text
  5. 5.↵
    Feske S, Gwack Y, Prakriya M, et al. A mutation in Orai1 causes immune deficiency by abrogating CRAC channel function. Nature. 2006;441(7090):179–185. doi:10.1038/nature04702
    OpenUrlCrossRefPubMedWeb of Science
  6. 6.↵
    Picard C, McCarl C-A, Papolos A, et al. STIM1 Mutation Associated with a Syndrome of Immunodeficiency and Autoimmunity. N Engl J Med. 2009;360(19):1971–1980. doi:10.1056/NEJMoa0900082
    OpenUrlCrossRefPubMedWeb of Science
  7. 7.↵
    Kerkhofs M, Seitaj B, Ivanova H, Monaco G, Bultynck G, Parys JB. Pathophysiological consequences of isoform-specific IP3 receptor mutations. Biochim Biophys Acta - Mol Cell Res. 2018;1865(11):1707–1717. doi:10.1016/j.bbamcr.2018.06.004
    OpenUrlCrossRefPubMed
  8. 8.↵
    Terry LE, Alzayady KJ, Furati E, Yule DI. Inositol 1,4,5-trisphosphate Receptor Mutations associated with Human Disease. Messenger (Los Angeles, Calif Print). 2018;6(1-2):29–44. http://www.ncbi.nlm.nih.gov/pubmed/30197841
    OpenUrl
  9. 9.↵
    Parolin Schnekenberg R, Perkins EM, Miller JW, et al. De novo point mutations in patients diagnosed with ataxic cerebral palsy. Brain. 2015;138(7):1817–1832. doi:10.1093/brain/awv117
    OpenUrlCrossRefPubMed
  10. 10.↵
    Tada M, Nishizawa M, Onodera O. Roles of inositol 1,4,5-trisphosphate receptors in spinocerebellar ataxias. Neurochem Int. 2016;94:1–8. doi:10.1016/j.neuint.2016.01.007
    OpenUrlCrossRef
  11. 11.↵
    Klar J, Hisatsune C, Baig SM, et al. Abolished InsP3R2 function inhibits sweat secretion in both humans and mice. J Clin Invest. 2014;124(11):4773–4780. doi:10.1172/JCI70720
    OpenUrlCrossRefPubMed
  12. 12.↵
    Rönkkö J, Molchanova S, Revah-Politi A, et al. Dominant mutations in ITPR3 cause Charcot-Marie-Tooth disease. Ann Clin Transl Neurol. 2020;7(10):1962–1972. doi:10.1002/acn3.51190
    OpenUrlCrossRef
  13. 13.↵
    Sugawara H. Genetic evidence for involvement of type 1,type 2and type 3inositol 1,4,5-trisphosphate receptors in signal transduction through the B-cell antigen receptor. EMBO J. 1997;16(11):3078–3088. doi:10.1093/emboj/16.11.3078
    OpenUrlAbstract/FREE Full Text
  14. 14.
    Ouyang K, Leandro Gomez-Amaro R, Stachura DL, et al. Loss of IP3R-dependent Ca2+ signalling in thymocytes leads to aberrant development and acute lymphoblastic leukemia. Nat Commun. 2014;5(1):4814. doi:10.1038/ncomms5814
    OpenUrlCrossRefPubMed
  15. 15.↵
    Tang H, Wang H, Lin Q, et al. Loss of IP 3 Receptor–Mediated Ca 2+ Release in Mouse B Cells Results in Abnormal B Cell Development and Function. J Immunol. 2017;199(2):570–580. doi:10.4049/jimmunol.1700109
    OpenUrlAbstract/FREE Full Text
  16. 16.↵
    Foskett JK, White C, Cheung K-H, Mak D-OD. Inositol trisphosphate receptor Ca2+ release channels. Physiol Rev. 2007;87(2):593–658. doi:10.1152/physrev.00035.2006
    OpenUrlCrossRefPubMed
  17. 17.↵
    Vaser R, Adusumalli S, Leng SN, Sikic M, Ng PC. SIFT missense predictions for genomes. Nat Protoc. 2016;11(1):1–9. doi:10.1038/nprot.2015.123
    OpenUrlCrossRefPubMed
  18. 18.↵
    Quang D, Chen Y, Xie X. DANN: a deep learning approach for annotating the pathogenicity of genetic variants. Bioinformatics. 2015;31(5):761–763. doi:10.1093/bioinformatics/btu703
    OpenUrlCrossRefPubMed
  19. 19.↵
    Zhang P, Bigio B, Rapaport F, et al. PopViz: a webserver for visualizing minor allele frequencies and damage prediction scores of human genetic variations. Stegle O, ed. Bioinformatics. 2018;34(24):4307–4309. doi:10.1093/bioinformatics/bty536
    OpenUrlCrossRef
  20. 20.↵
    Paknejad N, Hite RK. Structural basis for the regulation of inositol trisphosphate receptors by Ca2+ and IP3. Nat Struct Mol Biol. 2018;25(8):660–668. doi:10.1038/s41594-018-0089-6
    OpenUrlCrossRefPubMed
  21. 21.↵
    Azumaya CM, Linton EA, Risener CJ, Nakagawa T, Karakas E. Cryo-EM structure of human type-3 inositol triphosphate receptor reveals the presence of a self-binding peptide that acts as an antagonist. J Biol Chem. 2020;295(6):1743–1753. doi:10.1074/jbc.RA119.011570
    OpenUrlAbstract/FREE Full Text
  22. 22.↵
    Soulsby MD, Wojcikiewicz RJH. The type III inositol 1,4,5-trisphosphate receptor is phosphorylated by cAMP-dependent protein kinase at three sites. Biochem J. 2005;392(3):493–497. doi:10.1042/BJ20051325
    OpenUrlAbstract/FREE Full Text
  23. 23.↵
    Hornbeck P V., Zhang B, Murray B, Kornhauser JM, Latham V, Skrzypek E. PhosphoSitePlus, 2014: mutations, PTMs and recalibrations. Nucleic Acids Res. 2015;43(D1):D512–D520. doi:10.1093/nar/gku1267
    OpenUrlCrossRefPubMed
  24. 24.↵
    Fan G, Baker MR, Wang Z, et al. Cryo-EM reveals ligand induced allostery underlying InsP3R channel gating. Cell Res. 2018;28(12):1158–1170. doi:10.1038/s41422-018-0108-5
    OpenUrlCrossRefPubMed
  25. 25.↵
    Alzayady KJ, Wang L, Chandrasekhar R, Wagner LE, Van Petegem F, Yule DI. Defining the stoichiometry of inositol 1,4,5-trisphosphate binding required to initiate Ca 2+ release. Sci Signal. 2016;9(422):ra35–ra35. doi:10.1126/scisignal.aad6281
    OpenUrlAbstract/FREE Full Text
  26. 26.↵
    Wahl MI, Jones GA, Nishibe S, Rhee SG, Carpenter G. Growth factor stimulation of phospholipase C-gamma 1 activity. Comparative properties of control and activated enzymes. J Biol Chem. 1992;267(15):10447–10456. http://www.ncbi.nlm.nih.gov/pubmed/1316902
    OpenUrlAbstract/FREE Full Text
  27. 27.↵
    Jasiulewicz A, Lisowska KA, Pietruczuk K, Frąckowiak J, Fulop T, Witkowski JM. Homeostatic ‘bystander’ proliferation of human peripheral blood B cells in response to polyclonal T-cell stimulation in vitro. Int Immunol. 2015;27(11):579–588. doi:10.1093/intimm/dxv032
    OpenUrlCrossRefPubMed
  28. 28.↵
    Prole DL, Taylor CW. Inositol 1,4,5‐trisphosphate receptors and their protein partners as signalling hubs. J Physiol. 2016;594(11):2849–2866. doi:10.1113/JP271139
    OpenUrlCrossRef
  29. 29.↵
    Parys JB, Vervliet T. New Insights in the IP3 Receptor and Its Regulation. In: ; 2020:243–270. doi:10.1007/978-3-030-12457-1_10
    OpenUrlCrossRef
  30. 30.↵
    Marchant JS, Taylor CW. Cooperative activation of IP3 receptors by sequential binding of IP3 and Ca2+ safeguards against spontaneous activity. Curr Biol. 1997;7(7):510–518. doi:10.1016/S0960-9822(06)00222-3
    OpenUrlCrossRefPubMedWeb of Science
  31. 31.↵
    Futatsugi A. IP3 Receptor Types 2 and 3 Mediate Exocrine Secretion Underlying Energy Metabolism. Science (80-). 2005;309(5744):2232–2234. doi:10.1126/science.1114110
    OpenUrlAbstract/FREE Full Text
  32. 32.↵
    1. Clapham DE
    Uchida K, Aramaki M, Nakazawa M, et al. Gene Knock-Outs of Inositol 1,4,5-Trisphosphate Receptors Types 1 and 2 Result in Perturbation of Cardiogenesis. Clapham DE, ed. PLoS One. 2010;5(9):e12500. doi:10.1371/journal.pone.0012500
    OpenUrlCrossRefPubMed
  33. 33.↵
    des Georges A, Clarke OB, Zalk R, et al. Structural Basis for Gating and Activation of RyR1. Cell. 2016;167(1):145-157.e17. doi:10.1016/j.cell.2016.08.075
    OpenUrlCrossRefPubMed
  34. 34.↵
    Bhanumathy C, da Fonseca Pca, Morris EP, Joseph SK. Identification of Functionally Critical Residues in the Channel Domain of Inositol Trisphosphate Receptors. J Biol Chem. 2012;287(52):43674–43684. doi:10.1074/jbc.M112.415786
    OpenUrlAbstract/FREE Full Text
  35. 35.↵
    Feske S, Müller JM, Graf D, et al. Severe combined immunodeficiency due to defective binding of the nuclear factor of activated T cells in T lymphocytes of two male siblings. Eur J Immunol. 1996;26(9):2119–2126. doi:10.1002/eji.1830260924
    OpenUrlCrossRefPubMedWeb of Science
  36. 36.↵
    Lian J, Cuk M, Kahlfuss S, et al. ORAI1 mutations abolishing store-operated Ca2+ entry cause anhidrotic ectodermal dysplasia with immunodeficiency. J Allergy Clin Immunol. 2018;142(4):1297-1310.e11. doi:10.1016/j.jaci.2017.10.031
    OpenUrlCrossRef
  37. 37.↵
    Schaballie H, Rodriguez R, Martin E, et al. A novel hypomorphic mutation in STIM1 results in a late-onset immunodeficiency. J Allergy Clin Immunol. 2015;136(3):816-819.e4. doi:10.1016/j.jaci.2015.03.009
    OpenUrlCrossRef
  38. 38.↵
    Parry DA, Holmes TD, Gamper N, et al. A homozygous STIM1 mutation impairs store-operated calcium entry and natural killer cell effector function without clinical immunodeficiency. J Allergy Clin Immunol. 2016;137(3):955-957.e8. doi:10.1016/j.jaci.2015.08.051
    OpenUrlCrossRefPubMed
  39. 39.↵
    Fuchs S, Rensing-Ehl A, Speckmann C, et al. Antiviral and Regulatory T Cell Immunity in a Patient with Stromal Interaction Molecule 1 Deficiency. J Immunol. 2012;188(3):1523–1533. doi:10.4049/jimmunol.1102507
    OpenUrlAbstract/FREE Full Text
  40. 40.↵
    Rice L, Stockdale C, Berry I, et al. A Report of Novel STIM1 Deficiency and 6-Year Follow-Up of Two Previous Cases Associated with Mild Immunological Phenotype. J Clin Immunol. 2019;39(3):249–256. doi:10.1007/s10875-019-00618-3
    OpenUrlCrossRef
  41. 41.↵
    Emsley P, Lohkamp B, Scott WG, Cowtan K. Features and development of Coot. Acta Crystallogr D Biol Crystallogr. 2010;66(Pt 4):486–501. doi:10.1107/S0907444910007493
    OpenUrlCrossRefPubMedWeb of Science
  42. 42.↵
    Adams PD, Afonine P V, Bunkóczi G, et al. PHENIX: a comprehensive Python-based system for macromolecular structure solution. Acta Crystallogr D Biol Crystallogr. 2010;66(Pt 2):213–221. doi:10.1107/S0907444909052925
    OpenUrlCrossRefPubMedWeb of Science
  43. 43.↵
    Pettersen EF, Goddard TD, Huang CC, et al. UCSF Chimera?A visualization system for exploratory research and analysis. J Comput Chem. 2004;25(13):1605–1612. doi:10.1002/jcc.20084
    OpenUrlCrossRefPubMedWeb of Science
  44. 44.↵
    Parys JB, de Smedt H, Missiaen L, Bootman MD, Sienaert I, Casteels R. Rat basophilic leukemia cells as model system for inositol 1,4,5-trisphosphate receptor IV, a receptor of the type II family: functional comparison and immunological detection. Cell Calcium. 1995;17(4):239–249. doi:10.1016/0143-4160(95)90070-5
    OpenUrlCrossRefPubMedWeb of Science
  45. 45.↵
    Ma H-T, Venkatachalam K, Parys JB, Gill DL. Modification of store-operated channel coupling and inositol trisphosphate receptor function by 2-aminoethoxydiphenyl borate in DT40 lymphocytes. J Biol Chem. 2002;277(9):6915–6922. doi:10.1074/jbc.M107755200
    OpenUrlAbstract/FREE Full Text
  46. 46.↵
    Bultynck G, Szlufcik K, Kasri NN, et al. Thimerosal stimulates Ca2+ flux through inositol 1,4,5-trisphosphate receptor type 1, but not type 3, via modulation of an isoform-specific Ca2+-dependent intramolecular interaction. Biochem J. 2004;381(Pt 1):87–96. doi:10.1042/BJ20040072
    OpenUrlAbstract/FREE Full Text
  47. 47.↵
    Schindelin J, Arganda-Carreras I, Frise E, et al. Fiji: an open-source platform for biological-image analysis. Nat Methods. 2012;9(7):676–682. doi:10.1038/nmeth.2019
    OpenUrlCrossRefPubMedWeb of Science
  48. 48.↵
    Schneider CA, Rasband WS, Eliceiri KW. NIH Image to ImageJ: 25 years of image analysis. Nat Methods. 2012;9(7):671–675. doi:10.1038/nmeth.2089
    OpenUrlCrossRefPubMedWeb of Science
Back to top
Posted June 08, 2021.
Download PDF

Supplementary Material

Data/Code
Email
Disrupted Ca2+ homeostasis and immunodeficiency in patients with functional Inositol 1,4,5-trisphosphate receptor subtype 3 defects
Julika Neumann, Erika Van Nieuwenhove, Lara E Terry, Frederik Staels, Taylor R Knebel, Kirsten Welkenhuyzen, Mariah R Baker, Margaux Gerbaux, Mathijs Willemsen, John S. Barber, Irina I Serysheva, Liesbeth De Waele, François Vermeulen, Isabelle Meyts, David I Yule, Geert Bultynck, Rik Schrijvers, Stephanie Humblet-Baron, Adrian Liston
medRxiv 2021.05.29.21257775; doi: https://doi.org/10.1101/2021.05.29.21257775
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area