- Open Access
Sequencing the CaSR locus in Pakistani stone formers reveals a novel loss-of-function variant atypically associated with nephrolithiasis
BMC Medical Genomics volume 14, Article number: 266 (2021)
Nephrolithiasis (NL) affects 1 in 11 individuals worldwide and causes significant morbidity and cost. Common variants in the calcium sensing receptor gene (CaSR) have been associated with NL. Rare inactivating CaSR variants classically cause hyperparathyroidism, hypercalcemia and hypocalciuria. However, NL and familial hypercalciuria have been paradoxically associated with select inactivating CaSR variants in three kindreds from Europe and Australia.
To discover novel NL-associated CaSR variants from a geographically distinct cohort, 57 Pakistani families presenting with pediatric onset NL were recruited. The CaSR locus was analyzed by directed or exome sequencing.
We detected a heterozygous and likely pathogenic splice variant (GRCh37 Chr3:122000958A>G; GRCh38 Chr3:12228211A>G; NM_000388:c.1609-2A>G) in CaSR in one family with recurrent calcium oxalate stones. This variant would be predicted to cause exon skipping and premature termination (p.Val537Metfs*49). Moreover, a splice variant of unknown significance in an alternative CaSR transcript (GRCh37 Chr3:122000929G>C; GRCh38 Chr3:122282082G >C NM_000388:c.1609-31G >C NM_001178065:c.1609-1G >C) was identified in two additional families.
Sequencing of the CaSR locus in Pakistani stone formers reveals a novel loss-of-function variant, expanding the connection between the CaSR locus and nephrolithiasis.
Nephrolithiasis (NL) affects 1 in 11 individuals during their lifetime [1, 2]. NL is associated with significant patient morbidity, recurrence, and healthcare costs . The causes of NL are not well understood. Formerly, monogenic causes of NL were thought to be restricted to rare tubulopathies and genetic syndromes. However, we demonstrated that a causative monogenic mutation can be detected in one of 30 known NL causing genes in ~ 6 to 11% of adult and ~ 17 to 29% pediatric cases [4,5,6,7].
The calcium sensing receptor (CaSR) gene regulates calcium homeostasis , and rare variants in CaSR cause calcium disorders [9, 10]. Inactivating CaSR variants are associated with hyperparathyroidism, hypercalcemia, and hypocalciuria with dominant (OMIM: 145980) and recessive (OMIM: 239200) modes of inheritance . Heterozygous activating CaSR variants cause autosomal dominant hypoparathyroidism, hypocalcemia and hypercalciuria (OMIM: 601198) [9, 11].
NL has been associated with common CaSR variants, which cause either increased activity or decreased expression [12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27]. NL is less typically the presenting problem in patients with calcium disorders caused by rare CaSR variants. 3.5% of cases with activating CaSR variants have NL or nephrocalcinosis, but this is typically an iatrogenic complication of treatment . NL has also been observed in three atypical kindreds from Europe and Australia with dominant inactivating variants (L650P, F881L, K336del), hypercalcemia, and familial hypercalciuria [29, 30]. These kindreds suggest that the functions of the calcium sensing receptor in serum calcium homeostasis and parathyroid hormone regulation can be uncoupled from its role in renal calcium excretion. Nevertheless, these position-specific inactivating variants are exceedingly rare, given that a deleterious CaSR variant was not detected in a worldwide cohort of 697 NL families using gene panel or exome sequencing [4,5,6,7]. Therefore, discovery of novel NL-associated alleles can provide a deeper understanding of the role of CaSR in calcium handling and kidney stone disease.
To uncover novel NL-associated variants, 57 Pakistani families presenting with pediatric onset NL were recruited. The CaSR locus was analyzed by directed sequencing in 20 families (47 individuals) and from exome sequencing data in an additional 37 families (46 individuals). We detected a novel heterozygous and likely pathogenic splice variant (NM_000388:c.1609-2A>G) in CaSR in one family with recurrent calcium oxalate stones. This variant would be predicted to cause skipping of exon 6 and premature termination (p.Pro537Metfs*49). Moreover, a splice variant of unknown significance in an alternative CaSR transcript (NM_001178065:c.1609-1G>C) was identified in two additional families. In summary, sequencing of the CaSR locus in Pakistani stone formers reveals a novel loss-of-function variant atypically associated with nephrolithiasis.
Subject recruitment and DNA sample extraction
The period of recruitment was 01/2013 to 08/2017. Subjects were evaluated by a Urologist for an NL episode at the Lahore General Hospital, Services Hospital, and associated clinics in Punjab, Pakistan. For inclusion, subjects presented with evidence of NL on ultrasonography and had onset of NL before age 21 years. Results of stone analyses (performed by Fourier transform infrared spectroscopy) were recovered where available. Subjects were excluded if the primary physician determined that the NL episode was potentially caused by an underlying medical condition (gastrointestinal disease, diuretic use), primary non-calcium metabolic disorder (e.g. primary or secondary hyperoxaluria), or cystic kidney disease (e.g. the presence of at least one renal cyst). At enrollment, informed consent, clinical data, pedigree information, and DNA samples from subjects were collected.
Blood samples from affected families were stored in ethylenediaminetetraacetic acid (EDTA) vacutainers, and DNA was extracted from these samples using an organic extraction method. DNA was stored in double-distilled deionized water or tris(hydroxymethyl)aminomethane (Tris) EDTA (TE) buffer (10 mM Tris and 0.1 mM EDTA).
Directed CaSR gene amplification
All coding exons and flanking regions of the CaSR gene were amplified by PCR using 15 pairs of primers. The annealing temperature, amplification product size, and PCR optimization conditions for each primer are given in Additional file 5: Table S1 and previously described . These amplified products were sequenced using the Sanger sequencing method. The sequencing data was analyzed using the CLC Genomics Workbench (version 6.5.2) software (CLC Bio, Aarhus, Denmark). The data was aligned to the CaSR reference genome sequence (chr3:121,902,530–122,005,342 forward strand; GRCh37:CM000665.1). Variants from the reference genome were analyzed as shown in Additional file 1: Figure S1A and as described in the Variant Analysis section below.
Genomic DNA was isolated from blood lymphocytes and subjected to exome capture using Agilent SureSelect (Santa Clara, CA) and human exome capture arrays (Life Technologies, Carlsbad, CA) followed by next-generation sequencing on the Illumina HiSeq sequencing platform (San Diego, CA). Sequence reads were mapped to the human reference genome assembly (NCBI build 3/hg19) using CLC Genomics Workbench version 6.5.2 software (CLC bio, Aarhus, Denmark). Following alignment to the human reference genome, variants within CaSR were filtered as previously described , and as described in the Variant Analysis section below, and in Additional file 1: Figure S1B. In any families with rare potentially deleterious variants, the exome data was further analyzed for causative variants in thirty established known NL-disease genes as previously described .
Variants were filtered as followed (see Additional file 1: Figure S1). First, different population databases were queried to include only rare alleles (minor allele frequency < 1%) including the Exome Sequencing Project (http://evs.gs.washington.edu/EVS), Exome Aggregation Consortium (http://exac.broadinstitute.org), Genome Aggregation Database (gnomAD; http://gnomad.broadinstitute.org), and 1000 Genomes Project (http://www.internationalgenome.org/1000-genomes-browsers). Synonymous and intronic variants that were not located within the splice site regions were excluded. Surviving variants were then ranked based on their likelihood to cause disease, taking into consideration evolutionary conservation among orthologs using the ENSEMBL Genome Browser (http://www.ensembl.org), and assembled using Clustal Omega (https://www.ebi.ac.uk/Tools/msa/clus-talo/). Three web-based programs were used to establish pathogenicity prediction scores of missense variants: Mutation Taster (http://www.mutationtaster.org), SIFT (http://sift.bii.a-star.edu.sg/), and PolyPhen-2 (http://genetics.bwh.harvard.edu/pph2). Splice variants were evaluated using the web-based programs MaxEnt (http://hollywood.mit.edu/burgelab/maxent/Xmaxentscan_scoreseq.html), NNSPLICE (https://www.fruitfly.org/seq_tools/splice.html), and splice site finder (http://www.umd.be/searchSpliceSite.html). The remaining variants were further evaluated by reviewing the existing literature and determining the phenotypic match. We consulted clinician-scientists and geneticists who are knowledgeable about the clinical phenotypes and pedigree structure and are experienced in identifying exome level mutations. In addition, the ACMG guidelines for variant classification were applied, and variants were disease-causing if they were classified as pathogenic or likely pathogenic.
Pediatric nephrolithiasis cohort characteristics
We recruited 57 families with pediatric nephrolithiasis. 31 probands were female, while 26 were male. 23 probands (40.4%) were born from a consanguineous union based on self-reporting that parents were first-degree cousins.
Directed and exome sequencing of CaSR locus
We performed directed sequencing of CaSR coding regions and adjacent splice sites in the first 20 families (47 individuals) with pediatric NL (Fig. 1). We identified rare variants within coding or splice regions in two families (KS-72 and KS-6). We, subsequently, queried exome sequencing data from an additional 37 Pakistani families (46 families) with pediatric NL and detected a rare variant within a splice region within one family (KS-71) (Fig. 1).
In family KS-72, there were three affected family members (Fig. 2; Table 1). The deceased mother in this pedigree had a history of NL without further clinical or genetic data. Her two sons developed bilateral NL with onset at age 18 years (II-3) and 19 years (II-4). Both had two stone episodes thus far in their lifetime. Stone composition analysis for subjects II-3 and II-4 revealed predominantly calcium oxalate nephrolithiasis (Table 1). Targeted sequencing detected a deleterious essential splice site variant (GRCh37 Chr3:122000958A>G; GRCh38 Chr3:12228211A>G; NM_000388:c.1609-2A>G; NM_001178065:c.1637A>G) in CaSR (Fig. 2, Additional file 3: Figure S3, Additional file 4: Figure S4; Table 1). This variant was deemed to be deleterious because it is (1) predicted to impair splicing of the predominantly expressed transcript NM_000388 [32,33,34] and (2) absent from population genomic databases ExAC and gnomAD (Tables 1, Additional file 7: S3). This variant would be predicted to cause skipping of exon 6 and premature termination (p. Val537Metfs*49) prior to the receptor’s seven transmembrane spanning domain, similar to a previous example from the CaSR locus . Of note, this variant would be predicted to cause a missense variant in the alternative transcript NM_001178065, which does encode a functional calcium sensing receptor but has only the third highest mRNA expression in the kidney relative to other alternative CaSR transcripts [32,33,34]. This missense variant was deemed deleterious in only one of three in silico prediction tools (Table 1). Sanger sequencing confirmed segregation of this variant in both affected siblings and lack of segregation with the unaffected father (Fig. 2). This variant, therefore, met the ACMG criteria  for a likely pathogenic variant.
In family KS-6, four family members had a history of urinary stones: II-2 (age of onset 18 years), the identical twins II-4 and II-5 (age of onset 17 years), and their mother I-2 (Additional file 2: Figure S2A, Additional file 6: Table S2). Stone composition analysis revealed predominantly calcium oxalate stones in all three affected siblings (Additional file 6: Table S2). All four affected family members report multiple NL recurrences since recruitment. Targeted sequencing detected a rare splice site variant (GRCh37 Chr3:122000929G >C GRCh38 Chr3:122282082G >C NM_000388:c.1609-31G >C NM_001178065:c.1609-1G>C) in CaSR (Additional file 6: Table S2; Additional file 2: Figure S2A, Additional file 3: S3). This variant is not predicted to effect splicing of the canonical transcript NM_000388, but rather impact an essential acceptor splice site position of exon 6 for the non-canonical transcript NM_001178065. This latter effect would result in premature termination (p.Pro537Metfs*49). This variant is rare in genome population databases overall. However, there is one homozygote and five heterozygotes present in the South Asian populations of these databases with no alleles identified in other populations (Additional file 6: Tables S2, Additional file 7: S3). Segregation analysis demonstrated that the variant was detected in all three affected siblings and not an unaffected sibling II-3 (Additional file 2: Figure S2A). The variant was not detected in the unaffected father, suggesting it may have been inherited from the affected mother. Because the biological importance of the alternative transcript NM_001178065 remains unclear, this variant was deemed to have unknown significance based on ACMG criteria .
In family KS-71, there were three affected family members (Additional file 2: Figure S2B; Additional file 6: Table S2). Subject II-5 had calcium oxalate NL beginning at age 20 years and had four NL episodes (Additional file 6: Table S2). Exome sequence data from this subject revealed the same rare splice site variant of unknown significance (NM_000388:c.1609-31G >C NM_001178065:c.1609-1G>C) in CaSR as was detected in family KS-6 (Additional file 6: Table S2). Additional DNA samples were not available from other family members for segregation analysis.
Nephrolithiasis affects 1 in 11 individuals worldwide. While NL has been associated with common CaSR variants [12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27], only specific rare inactivating CaSR variants in European and Australian kindreds cause NL [28,29,30]. These rare inactivating NL-associated variants are critical to discover and investigate, as they suggest the functions of the calcium sensing receptor in serum calcium homeostasis and parathyroid hormone regulation can be uncoupled from its role in renal calcium excretion. Discovery of novel NL-associated alleles can, therefore, provide a deeper understanding of the role of CaSR in calcium handling and kidney stone disease. For these reasons, the CaSR locus was sequenced in a cohort of 57 Pakistani families with pediatric NL, revealing rare loss-of-function variants in Pakistani kindreds with pediatric nephrolithiasis.
First, we detected a novel heterozygous and likely pathogenic splice variant (NM_000388:c.1609-2A>G) in CaSR in one family (KS-72) with pediatric calcium oxalate stones. This variant would be predicted to cause loss-of-function. NL has been observed in three atypical adult cases with dominant inactivating variants (K336del, L650P, F881L), hypercalcemia, and paradoxically familial hypercalciuria [29, 30]. These variants are spatially distinct and are predicted to impact the extracellular receptor domain, the transmembrane domain, or the intracellular C-terminal tail of the calcium sensing receptor, respectively. The variant in KS-72 is predicted to cause premature truncation of the receptor prior to its seven transmembrane spanning domain, which would be potentially more deleterious than the previous three variants identified in adult patients. A limitation of our study was that serum and urinary laboratory studies (e.g. PTH levels, serum calcium or phosphorous levels, urinary calcium excretion) to further characterize the NL-associated disease in KS-72 and other subjects were not available, because the primary physicians did not obtain these routinely during their clinical evaluation and because patients were not readily available to return to the clinical laboratory. Nevertheless, this variant was reported in ClinVar (VCV000854168.2) in a patient with a calcium disorder while our study was in progress. Overall, the novel variant expands the connection between rare inactivating CaSR variants and nephrolithiasis.
Secondly, we identified a splice variant of unknown significance in an alternative transcript (NM_001178065:c.1609-1G>C) in two additional families. Intriguingly, this variant would also be predicted to cause skipping of exon 6 in the alternative transcript and premature truncation prior to the receptor’s transmembrane domain. Alternatively, this variant may lead to further preference of the 5’ splice site of the canonical transcript NM_000388 exon 6. While aspects of our variant analysis would suggest this is a deleterious variant, further functional studies are warranted to understand if this alternative transcript plays an important or different role in calcium sensing relative to the canonical transcript (e.g. isoform-specific mouse knockouts).
Overall, our study reveals a novel rare CaSR variant associated with nephrolithiasis and suggests that specific inactivating variants in CaSR may have distinct effects on its role in calcium homeostasis.
Availability of data and material
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
- CaSR :
Calcium sensing receptor gene
Scales CD, Smith AC, Hanley JM, Saigal CS. Prevalence of kidney stones in the United States. Eur Urol. 2012;62:160–5.
Tasian GE, et al. Annual incidence of nephrolithiasis among children and adults in South Carolina from 1997 to 2012. Clin J Am Soc Nephrol. 2016;11:488–96.
Rule AD, et al. Kidney stones and the risk for chronic kidney disease. Clin J Am Soc Nephrol. 2009;4:804–11.
Amar A, et al. Gene panel sequencing identifies a likely monogenic cause in 7% of 235 Pakistani families with nephrolithiasis. Hum Genet. 2019;138:211–9.
Braun DA, et al. Prevalence of monogenic causes in pediatric patients with nephrolithiasis or nephrocalcinosis. Clin J Am Soc Nephrol. 2016;11:664–72.
Daga A, et al. Whole exome sequencing frequently detects a monogenic cause in early onset nephrolithiasis and nephrocalcinosis. Kidney Int. 2017. https://doi.org/10.1016/j.kint.2017.06.025.
Halbritter J, et al. Fourteen monogenic genes account for 15% of nephrolithiasis/nephrocalcinosis. J Am Soc Nephrol. 2015;26:543–51.
Toka HR, Pollak MR. The role of the calcium-sensing receptor in disorders of abnormal calcium handling and cardiovascular disease. Curr. Opin. Nephrol. Hypertens. 23; 2014.
Pollak MR, et al. Autosomal dominant hypocalcaemia caused by a Ca2+-sensing receptor gene mutation. Nat Genet. 1994;8:303–7.
Pollak MR, et al. Mutations in the human Ca2+-sensing receptor gene cause familial hypocalciuric hypercalcemia and neonatal severe hyperparathyroidism. Cell. 1993;75:1297–303.
Watanabe S, et al. Association between activating mutations of calcium-sensing receptor and Bartter’s syndrome. The Lancet. 2002;360:692–4.
Chen W-C, et al. The rs1256328 (ALPL) and rs12654812 (RGS14) polymorphisms are associated with susceptibility to calcium nephrolithiasis in a Taiwanese population. Sci Rep. 2019;9:17296.
Chou Y-H, et al. A genetic polymorphism (rs17251221) in the calcium-sensing receptor gene (CASR) is associated with stone multiplicity in calcium nephrolithiasis. PLOS ONE. 2011;6:e25227.
Corbetta S, et al. R990G polymorphism of the calcium-sensing receptor and renal calcium excretion in patients with primary hyperparathyroidism. Eur J Endocrinol Eur J Endocrinol. 2006;155:687–92.
Ding Q, et al. Calcium-sensing receptor genetic polymorphisms and risk of developing nephrolithiasis in a Chinese population. Urol Int. 2017;99:331–7.
Grzegorzewska AE, Bednarski D, Świderska M, Mostowska A, Jagodziński PP. The calcium-sensing receptor gene polymorphism rs1801725 and calcium-related phenotypes in hemodialysis patients. Kidney Blood Press Res. 2018;43:719–34.
Grzegorzewska AE, et al. Associations of the calcium-sensing receptor gene CASR rs7652589 SNP with nephrolithiasis and secondary hyperparathyroidism in haemodialysis patients. Sci Rep. 2016;6:35188.
Guha M, et al. Polymorphisms in CaSR and CLDN14 genes associated with increased risk of kidney stone disease in patients from the eastern part of India. PLoS ONE. 2015;10:e0130790.
Li H, Zhang J, Long J, Shi J, Luo Y. Calcium-sensing receptor gene polymorphism (rs7652589) is associated with calcium nephrolithiasis in the population of Yi nationality in Southwestern China. Ann Hum Genet. 2018;82:265–71.
Oddsson A, et al. Common and rare variants associated with kidney stones and biochemical traits. Nat Commun. 2015;6:7975.
Scillitani A, et al. Primary hyperparathyroidism and the presence of kidney stones are associated with different haplotypes of the calcium-sensing receptor. J Clin Endocrinol Metab. 2007;92:277–83.
Shakhssalim N, et al. Association between calcium-sensing receptor gene polymorphisms and recurrent calcium kidney stone disease: a comprehensive gene analysis. Scand J Urol Nephrol. 2010;44:406–12.
Vezzoli G, et al. Risk of nephrolithiasis in primary hyperparathyroidism is associated with two polymorphisms of the calcium-sensing receptor gene. J Nephrol. 2015;28:67–72.
Vezzoli G, Terranegra A, Soldati L. Calcium-sensing receptor gene polymorphisms in patients with calcium nephrolithiasis. Curr Opin Nephrol Hypertens. 21; 2012.
Vezzoli G, et al. Polymorphisms at the regulatory regions of the CASR gene influence stone risk in primary hyperparathyroidism. Eur J Endocrinol. 2011;164:421–7.
Vezzoli G, et al. Calcium kidney stones are associated with a haplotype of the calcium-sensing receptor gene regulatory region. Nephrol Dial Transplant. 2010;25:2245–52.
Vezzoli G, et al. Influence of calcium-sensing receptor gene on urinary calcium excretion in stone-forming patients. J Am Soc Nephrol. 2002;13:2517.
Pearce SHS, et al. A familial syndrome of hypocalcemia with hypercalciuria due to mutations in the calcium-sensing receptor. N Engl J Med. 1996;335:1115–22.
Carling T, et al. Familial hypercalcemia and hypercalciuria caused by a novel mutation in the cytoplasmic tail of the calcium receptor*. J Clin Endocrinol Metab. 2000;85:2042–7.
Warner J, et al. Genetic testing in familial isolated hyperparathyroidism: unexpected results and their implications. J Med Genet. 2004;41:155–60.
Álvarez-Hernández D, et al. A novel mutation in the calcium-sensing receptor responsible for autosomal dominant hypocalcemia in a family with two uncommon parathyroid hormone polymorphisms. J Mol Endocrinol. 2003;31(2):255–62.
Aida K, Koishi S, Tawata M, Onaya T. Molecular cloning of a putative Ca2+-sensing receptor cDNA from human kidney. Biochem Biophys Res Commun. 1995;214:524–9.
Garrett JE, et al. Molecular cloning and functional expression of human parathyroid calcium receptor cDNAs (∗). J Biol Chem. 1995;270:12919–25.
GTEx Consortium. The GTEx Consortium atlas of genetic regulatory effects across human tissues. Science 2020;369:1318.
D’Souza-Li L, Canaff L, Janicic N, Cole DEC, Hendy GN. An acceptor splice site mutation in the calcium-sensing receptor (CASR) gene in familial hypocalciuric hypercalcemia and neonatal severe hyperparathyroidism. Hum Mutat. 2001;18:411–21.
Richards S, et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med Off J Am Coll Med Genet. 2015;17:405–24.
We thank to patients and families who participated in this study and clinicians for their support in enrolment.
F.H. is the William E. Harmon Professor of Pediatrics. This research is supported by a grant from the National Institutes of Health to F.H. (DK-076683–13). A.J.M. was supported by NIH Institutional K12 Child Health Research Center Career Development Award (5K12HD052896-13), ASN Foundation for Kidney Research (FP01025169), and Boston Children's Hospital Manton Center for Orphan Disease Research Junior Faculty Award. The Yale Center for Mendelian Genomics (UM1HG006504) is funded by the National Human Genome Research Institute. The GSP Coordinating Center (U24 HG008956) contributed to cross-program scientific initiatives and provided logistical and general study coordination. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
Ethics approval and consent to participate
Human subject’s research performed in this study was in accordance with the ethical standards of and approved by the Institutional Review Boards of Boston Children’s Hospital and University of Veterinary and Animal Sciences Lahore, Pakistan. Informed consent was obtained from all NL subjects at time of enrolment to obtain and report clinical data and pedigree information in publications.
Consent for publication
F.H. is a co-founder of Goldfinch Biopharma Inc. No other authors declare potential competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Flow diagram of genetic variant filtering analysis.
Pedigree structure of Pakistani NL families with CaSR VUS.
Comparison of the canonical and alternative CaSR transcripts.
Impact of CaSR variant in family KS-72 on the canonical and alternative CaSR transcripts.
Primers sequence used for CaSR screening in nephrolithiasis patients.
Genetic description and clinical data of families with calcium-sensing receptor (CaSR) variants of unknown significance.
Allele frequency of NL-associated rare CaSR variants in South Asian and other populations of public genome databases.
About this article
Cite this article
Ullah, I., Ottlewski, I., Shehzad, W. et al. Sequencing the CaSR locus in Pakistani stone formers reveals a novel loss-of-function variant atypically associated with nephrolithiasis. BMC Med Genomics 14, 266 (2021). https://doi.org/10.1186/s12920-021-01116-5
- Calcium sensing receptor
- Rare disease