Table of Contents Author Guidelines Submit a Manuscript
The Scientific World Journal
Volume 2013, Article ID 541710, 7 pages
http://dx.doi.org/10.1155/2013/541710
Research Article

Human SLC26A1 Gene Variants: A Pilot Study

1Mater Research, Translational Research Institute, Woolloongabba, QLD 4102, Australia
2Department of Nephrology, University of Queensland at Princess Alexandra Hospital, Brisbane, QLD 4102, Australia
3Pathology Department, Mater Misercordiae Hospital, South Brisbane, QLD 4101, Australia

Received 5 August 2013; Accepted 4 September 2013

Academic Editors: L. Bao and S. Thamilselvan

Copyright © 2013 Paul A. Dawson et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Kidney stones are a global health problem, incurring massive health costs annually. Why stones recur in many patients remains unknown but likely involves environmental, physiological, and genetic factors. The solute linked carrier (SLC) 26A1 gene has previously been linked to kidney stones in mice. SLC26A1 encodes the sulfate anion transporter 1 (SAT1) protein, and its loss in mice leads to hyperoxaluria and calcium oxalate renal stones. To investigate the possible involvement of SAT1 in human urolithiasis, we screened the SLC26A1 gene in a cohort of 13 individuals with recurrent calcium oxalate urolithiasis, which is the commonest type. DNA sequence analyses showed missense mutations in seven patients: one individual was heterozygous R372H; 4 individuals were heterozygous Q556R; one patient was homozygous Q556R; and one patient with severe nephrocalcinosis (requiring nephrectomy) was homozygous Q556R and heterozygous M132T. The M132 amino acid in human SAT1 is conserved with 15 other species and is located within the third transmembrane domain of the predicted SAT1 protein structure, suggesting that this amino acid may be important for SAT1 function. These initial findings demonstrate genetic variants in SLC26A1 of recurrent stone formers and warrant wider independent studies of SLC26A1 in humans with recurrent calcium oxalate stones.

1. Introduction

The formation of insoluble calcium stones in the kidney and/or urinary tract is a major health problem worldwide with 6–9% of men and 3-4% of women affected [1]. This leads to significant cost in the form of expensive surgical procedures ( $2 billions annually in USA) and morbidity [2]. Most of these stones are composed of calcium and oxalate [3], the latter being a metabolic end product derived from the liver and diet [4], with high urinary oxalate levels (hyperoxaluria) being a major contributor to renal calcium oxalate deposition [5, 6]. We have previously shown normal levels of endogenous oxalate synthesis in recurrent stone formers [7] and individuals are generally compliant with dietary advice to reduce their oxalate intake, suggesting that most stones cannot be attributed to excess oxalate from the liver or diet. There is a genetic predisposition to renal stone formation with 15% of patients having an affected family member and 50% of patients form a second stone over the subsequent 10 years [8]. However, in most cases the aetiology of renal stone deposition is unknown [9].

Calcium oxalate stone formation is a multifactorial disease involving environmental, physiological, and genetic factors [12]. Numerous genetic studies have identified possible links between calcium oxalate stones and polymorphisms in genes that regulate calcium homeostasis [1215], mediate inflammatory responses [16, 17], and maintain the structure and function of renal tubules [1821]. Hyperoxaluria is the most frequent metabolic abnormality found in calcium oxalate stone formers and its basis can be idiopathic (most common), secondary (e.g., enteric hyperoxaluria due to intestinal dysfunction or resection), or due to primary hyperoxaluria (PH) [12]. The primary hyperoxalurias (PH 1 and PH 2) are rare forms of severe hyperoxaluria and calcium oxalate stones and are caused by mutations in the AGXT and GRHPR genes, respectively [22]. Monozygotic twin studies indicate that there is a genetic predisposition to idiopathic hyperoxaluria [12, 23], which may be relevant to the SAT1 sulfate anion transporter (protein is SAT1; gene is SLC26A1) that has been linked to calcium oxalate stone formation in mice [24].

SAT1 is expressed in the intestine, kidney, and liver where it exchanges sulfate ( ) with structurally similar anions including oxalate ( ) [2628]. This exchange process is proposed to be important for the removal of oxalate from the blood (via the intestine), as well as for maintaining circulating (via the kidneys) and hepatic levels of nutrient sulfate [4, 6, 24, 2931]. Loss of SAT1 in mice led to disturbances in sulfate homeostasis (hyposulfataemia and hypersulfaturia), as well as hyperoxalaemia and hyperoxaluria [24]. The finding of calcium oxalate stones in SAT1 null mice in this landmark paper was remarkable because rodents are extremely refractive to forming stones, even when mild-moderate hyperoxaluria occurs [32].

To date, the contribution of human SAT1 to calcium oxalate stones is unknown. As an initial approach, we compared the SLC26A1 gene sequence of 13 male patients with recurrent calcium oxalate stones to NCBI reference SLC26A1 sequences for human and 15 other species. Our findings of three nonsynonymous gene variants in SLC26A1, including a missense variant at a highly conserved residue in one patient with severe nephrocalcinosis, warrant wider independent studies of SLC26A1 in recurrent calcium oxalate stone formers.

2. Materials and Methods

2.1. Patient Cohort

The research protocol was approved by the Mater Human Research Ethics Committee. All persons gave their informed consent prior to their inclusion in the study. The eligibility criteria for our study were adult patients with recurrent calcium oxalate stones, based on laboratory stone analysis (calcium oxalate 75% of the stone), for whom the aetiology of calcium oxalate stones was unknown, and if the diagnosis of primary hyperoxaluria or enteric hyperoxaluria was able to be excluded by their treating clinician. Buccal swab samples for DNA analysis were collected but individual patient characteristics were not.

2.2. Genomic DNA Isolation and PCR and DNA Sequencing

Genomic DNA was prepared from buccal swabs according to the manufacturer’s protocol (Isohelix). Five overlapping fragments of SLC26A1 (Figure 1(a)) were PCR-amplified using 200 nM forward and reverse primers (Table 1) in a C1000 thermal cycler (BioRad). The thermal cycling protocol was 94°C for 1 min; 35 cycles of 94°C for 30 sec, 60°C for 30 sec, and 72°C for 2 min, followed by 1 cycle 72°C for 5 min. DNA sequencing of SLC26A1 was performed using 50 ng PCR product (purified with NucleoSpin Extract II, Macherey-Nagel) and 9.6 pmol forward and reverse primers (Table 1), ABI Prism BigDye Terminator Sequence chemistry and an AB3730x/96-capillary sequencer (Applied Biosystems). DNA sequence analyses included the SLC26A1 coding region (2 exons of 2103 bp total), exon-intron junctions, 5′-(378 bp) and 3′-(1189 bp) UTRs and 300 bp 5′-flanking region containing the minimal (135 bp) SLC26A1 promoter region [33]. Nucleotide changes in SLC26A1 were verified by DNA sequence analysis of 2 independent PCR-amplified products for both sense and antisense strands.

tab1
Table 1: Primers used for PCR and DNA sequencing.
fig1
Figure 1: Identification of variants in the human SLC26A1 gene and predicted amino acid substitutions. (a) Schematic of SLC26A1 showing 4 exons (boxes), transcription initiation site (+1, arrow), and the relative location of PCR primers (P1–P10) and PCR amplicon fragments (1–5) used for DNA sequence analysis. ((b)–(d)) Representative DNA sequence chromatographs showing variant (top panels) and control (bottom panels) sequences that predict (b) heterozygous M132T, (c) heterozygous R372H, and (d) homozygous Q556R. Reverse complement nucleotide sequences are shown. (e) Position of amino acid substitutions (•) in the predicted secondary structure model of the human SAT1 protein [25]. M132T within the third transmembrane domain; R372H in the fourth intracellular loop; and Q556R in the intracellular carboxy-terminal domain.
2.3. Bioinformatic Analyses

Prediction of protein stability changes for amino acid substitutions was determined as previously described [10] using SAT1 amino acid sequence NP_998778.1. Alignment of SAT1 amino acid sequences between human (NP_998778.1) and multiple species was generated using the Clustal W program [11]. Putative transcription factor binding sites in the human SLC26A1 5′-flanking region were determined using MatInspector [34] and NCBI reference sequence NT_037622.5.

3. Results and Discussion

To date, two variants (G368S and Q556R) in the human SLC26A1 gene have been reported, based on information from the NCBI single nucleotide polymorphism (SNP) database [25]. Whilst SLC26A1 has yet to be associated with any human disease, the mouse Slc26a1 gene has been linked to hyperoxaluria and calcium oxalate kidney stones [24]. Since human SAT1 and mouse Sat1 share high homology (78% amino acid identity) [33] and a similar tissue distribution with strong expression in the kidneys [27, 28], we propose a potential role for human SLC26A1 in renal stone pathology.

In this study, we screened the SLC26A1 gene in a cohort of recurrent calcium oxalate kidney stone formers and identified single nucleotide variants in the coding region of eight individuals (Table 2). Three of these variants predict amino acid substitutions (M132T, R372H, and Q556R), and four nucleotide variants are silent changes (T225T, S306S, R439R, and A625A). In addition, one variant was detected in the 5′-untranslated region of SLC26A1 exon 1 for three individuals, and four variants were found within the first 210 nucleotides of the SLC26A1 5′-flanking region for eight individuals (Table 2). The allele frequencies of 5 variants are shown in Table 3 and Figure 2. As yet, we do not know the allelic frequencies of SLC26A1 variants in recurrent calcium oxalate stone formers but that will be the next phase of our work.

tab2
Table 2: SLC26A1 sequence variants identified in a cohort of recurrent calcium oxalate stone formers.
tab3
Table 3: Predicted effect of missense gene variants on SAT1 protein structure stability, and allelic frequencies of these variants in the NCBI SNP database.
fig2
Figure 2: Sequence variants identified in this study and their location within putative transcription factor binding sites in the human SLC26A1 5′-flanking region. (a) Nucleotide sequence of the SLC26A1 5′-flanking region (NCBI accession number NT_037622.5) from −210 to is shown. Position denotes the putative transcription initiation site (NCBI accession number AK292747). Putative transcription factor binding motifs are boxed, and their position on the plus (+) or minus (−) strands is indicated: ZNF300: KRAB-containing zinc finger protein 300; PTX1: pituitary homeobox 1; HEN1: helix-loop-helix transcription factor; AP4: activator protein 4; BRN5: BM-5 POU domain factor; AP1: activator protein 1. (b) Allelic frequencies and accession numbers of nucleotide variants in the NCBI database: #accessed 7 November 2012. *Nucleotide numbering based on NCBI reference sequence NT_037622.

Destabilisation of protein structure is the most common mechanism by which an amino acid variant causes loss of protein function [35]. Using the Support Vector Machine program [10], we identified M132T and R372H to have a significant deleterious effect on the structural stability of human SAT1 protein (Table 3). In addition, the location of variants within membrane proteins can also be useful for predicting the functional consequences of amino acid changes. For example, nonconservative amino acid substitutions within the transmembrane domains (TMDs) of plasma membrane transporters frequently have a negative consequence on protein function [25]. The human SAT1 protein is predicted to contain 12 TMDs, with intracellular amino and carboxy-termini [28]. The nonconservative (hydrophobic to hydrophilic) M132T substitution is located within the third putative TMD of SAT1 (Figure 1(e)), suggesting that M132 may be important for maintaining the TMD structure of SAT1. R372H is located within the fourth intracellular loop, whereas Q556R is positioned in the intracellular carboxy-terminal sequence of the predicted secondary structure of SAT1 (Figure 1(e)), suggesting the possibility that these 2 amino acid substitutions may not disrupt the TMD structure of SAT1.

Conservation of amino acid sequences between multiple species can be useful to identify those residues that are important for protein function [36]. In order to determine if the nonsynonymous variants M132T, R372H, and Q556R occur at conserved positions of the SAT1 protein, we aligned the human SAT1 sequence with homologues from 15 different species (Table 4). Human M132 is conserved between all species, suggesting that this amino acid may play an important role in the structure and or function of SAT1. The R372H variant is a conservative amino acid change (i.e., both arginine and histidine are positively charged amino acids). Whilst human R372 is conserved with 11 other species (Table 4), this residue aligns with cysteine in rhesus monkey, serine in both rat and mouse, and glutamic acid in the African clawed frog (Table 4), implying that R372 may not be essential for SAT1 structure or function. Human Q556 aligns with an arginine residue in 10 species, including chimpanzee, rhesus monkey, dog, rat, and mouse, demonstrating that SAT1 sequence containing Q556R is conserved with other species. This finding implies that Q556 is a more recent amino acid in humans, due to genetic drift or selection, to become the major allele (0.65, Table 3). In addition, the presence of R556 in 10 species, together with the relatively high allelic frequency (0.35) of R556 in the human population (Table 3), suggests that the Q556R variant may not be pathogenetic in humans.

tab4
Table 4: Alignment of human SAT1 amino acid variants across species.

Previous studies identified the first 135 bp of the human SLC26A1 5′-flanking region to be sufficient for basal promoter activity [28]. In the SLC26A1 5′-flanking region of eight individuals, we identified four nucleotide variants (Table 2) which are located within putative transcription factor binding motifs (Figure 2(a)). One variant, g.977276A G, is located within an AP1 site, previously shown to be important for the transcriptional activity of human SLC26A1 [28]. However, this variant does not alter the essential core nucleotides of the AP1 binding motif sequence. Similarly, the other three variants found within the 5′-flanking region do not alter the essential core nucleotide sequences of the ZNF300, HEN1, or BRN5 transcription factor binding motifs, indicating that these variants may not be detrimental to the transcriptional activity of the SLC26A1 promoter. These findings, together with the high allelic frequencies (0.34 to 0.65, NCBI database) of g.977391G C, g.977343C A, and g.977276A G in the general population (Figure 2(b)), imply that these variants may not be pathogenetic.

4. Conclusions

This study provides preliminary and as yet unreported evidence of sequence variants in the human SLC26A1 gene of a small number of individuals with recurrent calcium oxalate kidney stones. Of great interest is the finding of a nonconservative amino acid substitution (M132T) within a predicted transmembrane domain of SAT1, in one patient with severe nephrocalcinosis. The novelty of these findings warrants further studies of SLC26A1 in humans with unexplained recurrent kidney stones.

Conflict of Interests

The authors declare that they have no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors thank S. Diaz-Guilas, S. Hancock, and the Mater Urology Unit for patient recruitment and sample collection. This study was supported by the Mater Medical Research Institute and the Mater Foundation.

References

  1. P. Hughes, “Kidney stones epidemiology,” Nephrology, vol. 12, no. 1, pp. S26–S30, 2007. View at Publisher · View at Google Scholar · View at Scopus
  2. Y. Lotan, “Economics and cost of care of stone disease,” Advances in Chronic Kidney Disease, vol. 16, no. 1, pp. 5–10, 2009. View at Publisher · View at Google Scholar · View at Scopus
  3. M. H. Gault and L. Chafe, “Relationship of frequency, age, sex, stone weight and composition in 15,624 stones: comparison of results for 1980 to 1983 and 1995 to 1998,” Journal of Urology, vol. 164, no. 2, pp. 302–307, 2000. View at Google Scholar · View at Scopus
  4. R. P. Holmes, H. O. Goodman, and D. G. Assimos, “Contribution of dietary oxalate to urinary oxalate excretion,” Kidney International, vol. 59, no. 1, pp. 270–276, 2001. View at Publisher · View at Google Scholar · View at Scopus
  5. S. R. Khan, P. A. Glenton, and K. J. Byer, “Dietary oxalate and calcium oxalate nephrolithiasis,” Journal of Urology, vol. 178, no. 5, pp. 2191–2196, 2007. View at Publisher · View at Google Scholar · View at Scopus
  6. J. M. Brown, A. H. Chalmers, D. M. Cowley, and B. C. McWhinney, “Enteric hyperoxaluria and urolithiasis,” The New England Journal of Medicine, vol. 315, no. 15, pp. 970–971, 1986. View at Google Scholar · View at Scopus
  7. B. C. McWhinney, S. L. Nagel, D. M. Cowley, J. M. Brown, and A. H. Chalmers, “Two-carbon oxalogenesis compared in recurrent calcium oxalate stone formers and normal subjects,” Clinical Chemistry, vol. 33, no. 7, pp. 1118–1120, 1987. View at Google Scholar · View at Scopus
  8. G. A. Rose, Urinary Stones: Clinical and Laboratory Aspects, MTP, Lancaster, UK, 1981.
  9. S. Park and M. S. Pearle, “Pathophysiology and management of calcium stones,” Urologic Clinics of North America, vol. 34, no. 3, pp. 323–334, 2007. View at Publisher · View at Google Scholar · View at Scopus
  10. J. Cheng, A. Randall, and P. Baldi, “Prediction of protein stability changes for single-site mutations using support vector machines,” Proteins, vol. 62, no. 4, pp. 1125–1132, 2006. View at Publisher · View at Google Scholar · View at Scopus
  11. D. G. Higgins, A. J. Bleasby, and R. Fuchs, “CLUSTAL V: improved software for multiple sequence alignment,” Computer Applications in the Biosciences, vol. 8, no. 2, pp. 189–191, 1992. View at Google Scholar · View at Scopus
  12. J. A. Sayer, “The genetics of nephrolithiasis,” Nephron, vol. 110, no. 2, pp. e37–e43, 2008. View at Publisher · View at Google Scholar · View at Scopus
  13. B. Y. Reed, W. L. Gitomer, H. J. Heller et al., “Identification and characterization of a gene with base substitutions associated with the absorptive hypercalciuria phenotype and low spinal bone density,” Journal of Clinical Endocrinology and Metabolism, vol. 87, no. 4, pp. 1476–1485, 2002. View at Publisher · View at Google Scholar · View at Scopus
  14. Y. Suzuki, A. Pasch, O. Bonny, M. G. Mohaupt, M. A. Hediger, and F. J. Frey, “Gain-of-function haplotype in the epithelial calcium channel TRPV6 is a risk factor for renal calcium stone formation,” Human Molecular Genetics, vol. 17, no. 11, pp. 1613–1618, 2008. View at Publisher · View at Google Scholar · View at Scopus
  15. G. Vezzoli, A. Terranegra, T. Arcidiacono et al., “R990G polymorphism of calcium-sensing receptor does produce a gain-of-function and predispose to primary hypercalciuria,” Kidney International, vol. 71, no. 11, pp. 1155–1162, 2007. View at Publisher · View at Google Scholar · View at Scopus
  16. R. D. Mittal, H. K. Bid, P. K. Manchanda, and R. Kapoor, “Association of interleukin-1β gene and receptor antagonist polymorphisms with calcium oxalate urolithiasis,” Journal of Endourology, vol. 21, no. 12, pp. 1565–1570, 2007. View at Publisher · View at Google Scholar · View at Scopus
  17. R. D. Mittal, H. K. Bid, and P. K. Manchanda, “Genotype and haplotype determination of IL1B (g. − 511C > T and g. + 3954C > T) and (IL1RN) in pediatric nephrolithiasis,” Clinica Chimica Acta, vol. 379, no. 1-2, pp. 42–47, 2007. View at Publisher · View at Google Scholar · View at Scopus
  18. M. Onaran, A. YIlmaz, I. Şen et al., “Heparan sulfate gene polymorphism in calcium oxalate nephrolithiasis,” Urological Research, vol. 37, no. 1, pp. 47–50, 2009. View at Publisher · View at Google Scholar · View at Scopus
  19. M. Onaran, A. Yilmaz, I. Şen et al., “A HindIII polymorphism of fibronectin gene is associated with nephrolithiasis,” Urology, vol. 74, no. 5, pp. 1004–1007, 2009. View at Publisher · View at Google Scholar · View at Scopus
  20. M. Ozturk, Y. Kordan, H. Cangul et al., “Association of urokinase gene 3′-UTR T/C polymorphism with calcium oxalate urolithiasis in children,” International Urology and Nephrology, vol. 40, no. 3, pp. 563–568, 2008. View at Publisher · View at Google Scholar · View at Scopus
  21. V. Tasic, D. Dervisov, S. Koceva, S. Weber, and M. Konrad, “Hypomagnesemia with hypercalciuria and nephrocalcinosis: case report and a family study,” Pediatric Nephrology, vol. 20, no. 7, pp. 1003–1006, 2005. View at Publisher · View at Google Scholar · View at Scopus
  22. B. Hoppe, B. B. Beck, and D. S. Milliner, “The primary hyperoxalurias,” Kidney International, vol. 75, no. 12, pp. 1264–1271, 2009. View at Publisher · View at Google Scholar · View at Scopus
  23. M. Monga, B. Macias, E. Groppo, and A. Hargens, “Genetic heritability of urinary stone risk in identical twins,” Journal of Urology, vol. 175, no. 6, pp. 2125–2128, 2006. View at Publisher · View at Google Scholar · View at Scopus
  24. P. A. Dawson, C. S. Russell, S. Lee et al., “Urolithiasis and hepatotoxicity are linked to the anion transporter Sat1 in mice,” The Journal of Clinical Investigation, vol. 120, no. 3, pp. 706–712, 2010. View at Publisher · View at Google Scholar
  25. P. A. Dawson and D. Markovich, “Genetic polymorphisms of human sulfate transporters,” Current Pharmacogenomics, vol. 5, no. 4, pp. 262–274, 2007. View at Publisher · View at Google Scholar · View at Scopus
  26. W. Krick, N. Schnedler, G. Burckhardt, and B. C. Burckhardt, “Ability of sat-1 to transport sulfate, bicarbonate, or oxalate under physiological conditions,” American Journal of Physiology: Renal Physiology, vol. 297, no. 1, pp. F145–F154, 2009. View at Publisher · View at Google Scholar · View at Scopus
  27. A. Lee, L. Beck, and D. Markovich, “The mouse sulfate anion transporter gene sat1 (Slc26a1): cloning, tissue distribution, gene structure, functional characterization, and transcriptional regulation by thyroid hormone,” DNA and Cell Biology, vol. 22, no. 1, pp. 19–31, 2003. View at Publisher · View at Google Scholar · View at Scopus
  28. R. R. Regeer, A. Lee, and D. Markovich, “Characterization of the human sulfate anion transporter (hsat-1) protein and gene (SAT1; SLC26A1),” DNA and Cell Biology, vol. 22, no. 2, pp. 107–117, 2003. View at Publisher · View at Google Scholar · View at Scopus
  29. J. M. Brown, G. Stratmann, and D. M. Cowley, “The variability and dietary dependence of urinary oxalate excretion in recurrent calcium stone formers,” Annals of Clinical Biochemistry, vol. 24, no. 4, pp. 385–390, 1987. View at Google Scholar · View at Scopus
  30. M. Hatch and R. W. Freel, “Intestinal transport of an obdurate anion: oxalate,” Urological Research, vol. 33, no. 1, pp. 1–16, 2005. View at Publisher · View at Google Scholar · View at Scopus
  31. S. R. Marengo and A. M. P. Romani, “Oxalate in renal stone disease: the terminal metabolite that just won't go away,” Nature Clinical Practice Nephrology, vol. 4, no. 7, pp. 368–377, 2008. View at Publisher · View at Google Scholar · View at Scopus
  32. Z. Jiang, J. R. Asplin, A. P. Evan et al., “Calcium oxalate urolithiasis in mice lacking anion transporter Slc26a6,” Nature Genetics, vol. 38, no. 4, pp. 474–478, 2006. View at Publisher · View at Google Scholar · View at Scopus
  33. A. Lee, P. A. Dawson, and D. Markovich, “NaSi-1 and Sat-1: structure, function and transcriptional regulation of two genes encoding renal proximal tubular sulfate transporters,” International Journal of Biochemistry and Cell Biology, vol. 37, no. 7, pp. 1350–1356, 2005. View at Publisher · View at Google Scholar · View at Scopus
  34. K. Cartharius, K. Frech, K. Grote et al., “MatInspector and beyond: promoter analysis based on transcription factor binding sites,” Bioinformatics, vol. 21, no. 13, pp. 2933–2942, 2005. View at Publisher · View at Google Scholar · View at Scopus
  35. Z. Wang and J. Moult, “SNPs, protein structure, and disease,” Human Mutation, vol. 17, no. 4, pp. 263–270, 2001. View at Publisher · View at Google Scholar · View at Scopus
  36. P. C. Ng and S. Henikoff, “SIFT: predicting amino acid changes that affect protein function,” Nucleic Acids Research, vol. 31, no. 13, pp. 3812–3814, 2003. View at Publisher · View at Google Scholar · View at Scopus