- Regular Article
Whole-exome sequencing reveals genetic variants associated with chronic kidney disease characterized by tubulointerstitial damages in North Central Region, Sri Lanka
Environmental Health and Preventive Medicine volume 20, pages 354–359 (2015)
The familial clustering observed in chronic kidney disease of uncertain etiology (CKDu) characterized by tubulointerstitial damages in the North Central Region of Sri Lanka strongly suggests the involvement of genetic factors in its pathogenesis. The objective of the present study is to use whole-exome sequencing to identify the genetic variants associated with CKDu.
Whole-exome sequencing of eight CKDu cases and eight controls was performed, followed by direct sequencing of candidate loci in 301 CKDu cases and 276 controls.
Association study revealed rs34970857 (c.658G > A/p.V220M) located in the KCNA10 gene encoding a voltage-gated K channel as the most promising SNP with the highest odds ratio of 1.74. Four rare variants were identified in gene encoding Laminin beta2 (LAMB2) which is known to cause congenital nephrotic syndrome. Three out of four variants in LAMB2 were novel variants found exclusively in cases.
Genetic investigations provide strong evidence on the presence of genetic susceptibility for CKDu. Possibility of presence of several rare variants associated with CKDu in this population is also suggested.
Chronic kidney disease (CKD) is a global public health issue due to increasing prevalence and economic burden incurred [1–3]. This has been attributed to an increasing incidence of conventional risk factors such as diabetes and hypertension in populations around the world. However, in the recent past involvement of genetic factors in pathogenesis of CKD has been highlighted and a number of genetic polymorphisms have been identified .
Increasing incidence of CKD during the past two decades was reported in the North Central Province and the neighboring provinces (hereafter referred to as North Central Region NCR) in Sri Lanka . Aturaliya and co-workers reported that major risk factors for CKD could not be identified for 87 % of these CKD patients through the available diagnostic tests of the time. This entity characterized by tubulointerstitial damages has since been termed CKD of uncertain etiology (CKDu) .
Over the years, several hypotheses have been proposed on etiology and pathogenesis of CKDu, mostly on possible involvement of environmental factors. Involvement of heavy metals such as Cadmium and Arsenic has been suggested but, the results have been inconsistent [6–9]. Fluoride has also been one of the earliest suggestions as an etiological factor for CKDu . In addition, involvement of herbicides, mycotoxins, pesticides, and algal toxins in CKDu has been proposed and debated in the recent past [9, 11]. Despite continuous efforts, the investigations to date have been unable to provide adequate scientific evidence to identify the exact cause/s of CKDu and a complex involvement of both genetic and environmental factors have been highly suspected .
In one of the earlier epidemiological studies, it has been postulated that a family history of renal dysfunction as a risk factor for CKDu. This has been largely supported by the prevalence of the disease in families or familial clustering [4, 12]. The inheritance pattern in CKDu did not support a single-gene disease but, matched that of multifactorial diseases such as diabetes, coronary heart disease and cancer . Thus, it could be surmised that the pathogenesis occurs in genetically susceptible individuals when exposed to one or more environmental toxin/s. With this working hypothesis, first, we screened a population for genetic polymorphisms that may predispose to CKDu using Genome-Wide Association Study (GWAS), which identified a GWAS significant single-nucleotide polymorphism in SLC13A3 gene coding sodium dicarboxylate cotransporter 3 located in renal tubular cells with a population attributable risk of 50 % .
GWAS of complex traits have been proven successful in identifying common variant associations. However, such methods often overlook large contributions of ethnicity-specific minor alleles, which can also cause substantial risk of complex trait diseases [14, 15]. As the second step, in the present study, whole-exome sequencing was performed to identify probable minor genetic variants associated with increased genetic susceptibility of this cohort for CKDu.
Materials and methods
Ethical approval for the study was obtained from the ethical committees of Kyoto University, Japan, and Faculty of Medicine, University of Peradeniya, Sri Lanka. All specimens, clinical reports, and questionnaire data were obtained following informed written consent.
A cohort of 311 CKDu male cases and 286 male controls from two highly affected areas (Medawachchiya and Girandurukotte) was established in 2011 for the GWAS (methodology of the recruitment process and sample collection is mentioned in detail in our latest publication) . CKDu cases are defined by the absence of known risk factors for CKD and having typical pathological lesion of tubulointerstitial fibrosis. Hypertensive nephropathy was excluded by the absence of history of chronic or uncontrolled hypertension before the diagnosis of CKDu. Diabetic nephropathy was excluded by the absence of history of diabetes preceding the diagnosis of CKDu and having glycosylated hemoglobin (HbA1c) level <6.5 % at the time of recruitment. Other known risk factors such as autoimmune diseases were also excluded by clinical history, examination findings, antibody assays and immunohistochemistry. Controls were the unrelated apparently healthy males from the same region without diabetes (no history of diabetes and HbA1c <6.5 % at the time of recruitment), hypertension (no history of hypertension and blood pressure <140/90 mmHg at the time of recruitment) and no evidence of renal function impairment (no history of renal diseases, serum creatinine <1.2 mg/dL and alpha 1 microglobulin <15.5 mg/L at the time of recruitment). For the whole-exome sequencing eight individuals with CKDu and eight controls were randomly selected to equally represent the two endemic areas.
Following informed written consent, blood samples (10 mL) were collected from peripheral veins into K-EDTA tubes. Serum was separated immediately by centrifugation at 3000 rpm for 10 min. Blood samples were immediately stored at 4 °C and then transferred to −30 °C within 6 h of collection. The samples were shipped to Japan at −20 °C and stored at −80 °C until analysis at the human specimen bank in Kyoto University, Japan. DNA was extracted from the whole blood samples using a commercially available extraction kit (QIAamp DNA Blood Mini Kit, Qiagen, Chatsworth, CA, USA) as per the manufacturer’s instructions. HiSeq 2000 platform was used to perform exome sequencing. The target regions were captured using SureSelect Human All Exon 50 M Kit (Agilent Technologies, Santa Clara, CA, USA) according to the manufacturer’s protocols. Briefly, genomic DNA was extracted from peripheral blood and randomly fragmented by acoustic fragmentation (Covaris) and purified using a QIAquick PCR Purification (Qiagen) kit. Then, adapters were ligated to both ends of the fragments. The resulting DNA library was purified using a QIAquick PCR Purification kit, amplified by ligation-mediated PCR and captured by hybridization to the SureSelect Biotinylated RNA library “baits” (Agilent) for enrichment. Captured ligation-mediated PCR products were subjected to an Aglient 2100 Bioanalyzer to estimate the magnitude of enrichment. Each captured library was then loaded on a Hiseq 2000 platform (Aligent), and paired-end sequencing was performed with read lengths of 101 bp using two channels. Sequence reads were mapped to the reference human genome (Ghr37/hg19) (http://hgdownload.cse.ucsc.edu/goldenPath/gh19/chromosomes/) using Burrows-Wheeler Aligner 0.6.2 (BWA) software (http://bio-bwa.sourceforge.net/index.shtml) [16, 17]. Single-nucleotide variants and small insertion/deletion were detected using The Genome Analysis Toolkit (GATK) software http://www.broadinstitute.org/gsa/wiki/index.php/The_Genome_Analysis_Toolkit) .
Bioinformatics analyses were mainly performed using SVS version 7.6.10, GoldenHelix software (Golden Helix, Inc., Bozeman, MT, USA). Figure 1 illustrates the step wise filtering process used to identify the candidate variants. The obtained sequence reads (SNPs and INDEL variants) were aligned to the human reference genome GRCh 37. Non-coding/splicing variants with poor read depths (<25) were excluded from further analysis. Synonymous variants were also filtered out. Coding/splicing, non-synonymous variants were isolated, and the sorting intolerant from tolerant (SIFT) tool [http://sift.bii.a-star.edu.sg] was used to predict SNPs with deleterious amino acid changes, and all SNPs with SIFT scores of >0.05 were eliminated .
Common variants with MAFs of ≥0.2 and indels with MAFs ≥0.01 that were present in the NHLBI exome sequencing project (ESP) database (data release-December 2011 https://esp.gs.washington.edu) were eliminated. The finalized SNPs and indels were used to check the gene associations with phenotypes using the kernel-based adaptive cluster method (KBAC) with a permutation (1,000 times) testing method . All SNPs and indels with values of P < 0.2 were selected, and the Endeavour software [http://homes.esat.kuleuven.be/~bioiuser/endeavour/index.php] was used to identify the genes related to renal function. Seventeen known genes associated with kidney diseases (FGF23 , UMOD [22, 23], PLA2R1 , NPHS1 , NPHS2 , ACTN4 , CD2AP , TRPC6 , LAMB2 [29, 30], PRCP , COL4A3 , CLCN5 , OCRL , CTNS , SLC3A1 , SLC7A9 , and SLC4A4 ) available in published data were used to educate the Endeavour software.
Selected 301 CKDu cases and 276 controls from the cohort (10 CKDu cases and 10 controls were eliminated during quality control process) used for GWAS analysis were screened for all the shortlisted SNPs in whole-exome sequencing using TaqMan probes or the restriction fragment length polymorphism method.
After removing the non-coding/splicing variants and variants with poor read depths (<25), 48,176 coding variants were left for further analysis. Among the remaining variants there were 22,919 synonymous variants which were also removed. Out of the 4,532 deleterious variants identified by the SIFT tool only 3,940 variants had MAF <0.2 for SNPs and <0.01 for Indels. These 3,940 variants were from 3,041 genes. KBAC method identified 120 genes with P < 0.2 for association out of which 7 genes were identified as related to kidney function by Endeavour software. The analysis process is illustrated in Fig. 1. Seven kidney-related genes (containing 11 deleterious SNPs) filtered through the analysis process were: PRCP (prolylcarboxypeptidase angiotensinase C, Chr 11); LAMB2 (laminin, beta 2 (laminin S), Chr 3); KNG1 (kininogen 1, Chr 3); BANK1 (B cell scaffold protein with ankyrin repeats, Chr 4); KCNA10 (potassium voltage-gated channel, Chr 1); SLC7A13 (sodium-independent aspartate/glutamate transporter, Chr 8); and TJP1 (tight junction protein 1,Chr 15) (Table 1). Out of the 11 SNPs, three were novel rare variants in the LAMB2 gene exclusively found in CKDu cases.
An association study was further conducted for these isolated 11 SNPs. The observed MAFs through direct sequencing of all the cases and controls as well as the odds ratios are summarized in Table 1. rs34970857 (c.658G > A/p.V220M) located in the KCNA10 gene encoding a voltage-gated K channel was the most suspected SNP with a significant odds ratio of 1.74 (P = 0.005 for QTL analysis; P = 0.06 for dichotomous analysis) (Table 1). The bioinformatics analysis by SIFT Score (http://sift.jcvi.org/) and Polyphen (2.2.2) (http://genetics.bwh.harvard.edu/pph2/) for p.V220M of KCNA10 predicted that the variant is possibly damaging and the missense mutations in LAMB2 gene are probably damaging.
Studies have consistently demonstrated the importance of genetic contributions to renal diseases [3, 37, 38]. In our previous investigations using GWAS, we reported a single common genetic variant in SLC13A3 gene associated with CKDu . Rare and generally deleterious genetic variants are also known to cause strong impact on risk of diseases in individual patients . Rapid advances in next-generation sequencing techniques such as whole-exome sequencing have facilitated the identification of these rare variants. In the present study, we identified one genetic variant in KCNA10 gene and four rare variants in LAMB2 gene which were not identified in our previous GWAS. The three rare variants in LAMB2 (c.3838G > A/E1280K, c.3583G > T/G1195W, c.1931G > A/R644H) and the variant identified in KCNA10 (c.658G > A/V220 M) were observed exclusively among the CKDu cases.
KCNA10 encodes a recently discovered voltage-gated Potassium (K) channel located in the heart, aorta glomerular endothelium, and apical membrane of renal proximal tubular cells [39, 40]. Its presence in endothelial and vascular smooth muscle cells supports the notion that it may also regulate vascular tone . Dysfunction of these K channels is known to be related to cardiovascular diseases, such as hypertension . Previously reported SLC13A3 gene is also reported to have suggestive association with blood pressure . Thus, both SLC13A3 and KCNA10 may play roles in blood pressure regulation and the pathogenesis of hypertension. In our previous study, as in the risk factor analyses for disease progression and mortality in CKDu, hypertension was identified as a major risk factor for disease progression . This evidence suggests that although clinically evident hypertension is not known to precede CKDu, hypertension may have an unidentified albeit an important role in the pathophysiology and progression of CKDu.
LAMB2 encodes laminin β2 is a component of laminin-521, which is an important constituent of the glomerular basement membrane . Detrimental missense mutations in the LAMB2 gene cause Pierson syndrome, a severe congenital nephrotic syndrome with ocular and neurologic defects [29, 30, 44, 45]. In the present study, we could not find extra-renal manifestations ocular and neurologic defects associated with mutations in LAMB2 gene. Family history has been reported to be a significant predictor for CKDu, and our investigations so far supports that genetic susceptibility plays an important role in pathogenesis of CKDu [12, 13]. In the present study, we identified four different rare variants by the exome sequencing of only a limited number of cases and controls. Thus, this suggests that there is a possibility for the existence of more rare variants in this population.
Previous evaluations of the pedigrees of CKDu patients have provided evidence in favor of a multifactorial disease rather than a single-gene disease . The past and present genetic analysis also shows that pathogenesis of CKDu cannot be completely explained by genetic factors but, genetic susceptibility enforces a discernible risk for disease occurrence.
Levey AS, Atkins R, Coresh J, Cohen EP, Collins AJ, Eckardt KU, et al. Chronic kidney disease as a global public health problem: approaches and initiatives—a position statement from Kidney Disease Improving Global Outcomes. Kidney Int. 2007;72(3):247–59.
Schieppati A, Remuzzi G. Chronic renal diseases as a public health problem: epidemiology, social, and economic implications. Kidney Int Suppl. 2005;98:S7–10.
Köttgen A, Pattaro C, Böger CA, Fuchsberger C, Olden M, Glazer NL, et al. New loci associated with kidney function and chronic kidney disease. Nat Genet. 2010;42(5):376–84.
Nanayakkara S, Senevirathna ST, Karunaratne U, Chandrajith R, Harada KH, Hitomi T, et al. Evidence of tubular damage in the very early stage of chronic kidney disease of uncertain etiology in the North Central Province of Sri Lanka: a cross-sectional study. Environ Health Prev Med. 2012;17(2):109–17.
Athuraliya NT, Abeysekera TD, Amerasinghe PH, Kumarasiri R, Bandara P, Karunaratne U, et al. Uncertain etiologies of proteinuric-chronic kidney disease in rural Sri Lanka. Kidney Int. 2011;80(11):1212–21.
Jayatilake N, Mendis S, Maheepala P. Mehta FR; CKDu National Research Project Team. Chronic kidney disease of uncertain aetiology: prevalence and causative factors in a developing country. BMC Nephrol. 2013;14:180.
Bandara JM, Wijewardena HV, Liyanege J, Upul MA, Bandara JM. Chronic renal failure in Sri Lanka caused by elevated dietary cadmium: trojan horse of the green revolution. Toxicol Lett. 2012;198(1):33–9.
Chandrajith R, Nanayakkara S, Itai K, Aturaliya TN, Dissanayake CB, Abeysekera T, et al. Chronic kidney diseases of uncertain etiology (CKDue) in Sri Lanka: geographic distribution and environmental implications. Environ Geochem Health. 2011;33(3):267–78.
Wimalawansa SJ. Escalating chronic kidney diseases of multi-factorial origin in Sri Lanka: causes, solutions, and recommendations. Environ Health Prev Med. 2014;19(6):375–94.
Ileperuma OA, Dharmagunawardhane HA, Herath KPRP. Dissolution of aluminium from sub-standard utensils under high fluoride stress: a possible risk factor for chronic renal failure in the North Central Province. J Natl Sci Found Sri Lanka. 2009;37(3):219–22.
Wickremasinghe AR, Peiris-John RJ, Wanigasuriya KP. Chronic kidney disease of unknown aetiology in the North Central Province of Sri Lanka: trying to unravel the mystery. Ceylon Med J. 2011;56(4):143–6.
Wanigasuriya KP, Peiris-John RJ, Wickremasinghe R, Hittarage A. Chronic renal failure in North Central Province of Sri Lanka: an environmentally induced disease. Trans R Soc Trop Med Hyg. 2007;101:1013–7.
Nanayakkara S, Senevirathna S, Abeysekera T, Chandrajith R, Ratnatunga N, Gunarathne E, et al. An integrative study of the genetic, social and environmental determinants of chronic kidney disease characterized by tubulointerstitial damages in the north central region of Sri Lanka. J Occup Health. 2014;56(1):28–38.
Casto AM, Feldman MW. Genome-wide association study SNPs in the human genome diversity project populations: does selection affect unlinked SNPs with shared trait associations? PLoS Genet. 2011;7(1):e1001266.
Majewski J, Schwartzentruber J, Lalonde E, Montpetit A, Jabado N. What can exome sequencing do for you? J Med Genet. 2011;48(9):580–9.
Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler Transform. Bioinformatics. 2009;25:1754–60.
Li H, Durbin R. Fast and accurate long-read alignment with Burrows-Wheeler Transform. Bioinformatics. 2010;26:589–95.
McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297–303.
Ng PC, Henikoff S. Predicting deleterious amino acid substitutions. Genome Res. 2001;11(5):863–74.
Liu DJ, Leal SM. A novel adaptive method for the analysis of next-generation sequencing data to detect complex trait associations with rare variants due to gene main effects and interactions. PLoS Genet. 2010;6(10):e1001156.
Dai B, David V, Martin A, Huang J, Li H, Jiao Y, et al. Comparative Transcriptome Analysis Identifying FGF23 Regulated Genes in the Kidney of a Mouse CKD Model. PLoS One. 2012;7(9):e44161.
Liu CT, Garnaas MK, Tin A, Kottgen A, Franceschini N, Peralta CA, et al. Genetic association for renal traits among participants of African ancestry reveals new loci for renal function. PLoS Genet. 2011;7(9):e1002264.
Köttgen A, Glazer NL, Dehghan A, Hwang SJ, Katz R, Li M, et al. Multiple loci associated with indices of renal function and chronic kidney disease. Nat Genet. 2009;41(6):712–7.
Stanescu HC, Arcos-Burgos M, Medlar A, Bockenhauer D, Kottgen A, Dragomirescu L, et al. Risk HLA-DQA1 and PLA(2)R1 alleles in idiopathic membranous nephropathy. N Engl J Med. 2011;364(7):616–26.
Abid A, Khaliq S, Shahid S, Lanewala A, Mubarak M, Hashmi S, et al. A spectrum of novel NPHS1 and NPHS2 gene mutations in pediatric nephrotic syndrome patients from Pakistan. Gene. 2012;502(2):133–7.
Bostrom MA, Perlegas P, Lu L, Hicks PJ, Hawkins G, Ng MC, et al. Relevance of the ACTN4 Gene in African-Americans with non-diabetic end-stage renal disease. Am J Nephrol. 2012;36(3):252–60.
Putaala H, Soininen R, Kilpeläinen P, Wartiovaara J, Tryggvason K. The murine nephrin gene is specifically expressed in kidney, brain and pancreas: inactivation of the gene leads to massive proteinuria and neonatal death. Hum Mol Genet. 2001;10(1):1–8.
Thilo F, Liu Y, Loddenkemper C, Schuelein R, Schmidt A, Yan Z, et al. VEGF regulates TRPC6 channels in podocytes. Nephrol Dial Transplant. 2012;27(3):921–9.
Matejas V, Hinkes B, Alkandari F, Al-Gazali L, Annexstad E, Aytac MB, et al. Mutations in the human laminin beta2 (LAMB2) gene and the associated phenotypic spectrum. Hum Mutat. 2010;31(9):992–1002.
Chen YM, Kikkawa Y, Miner JH. A missense LAMB2 mutation causes congenital nephrotic syndrome by impairing laminin secretion. J Am Soc Nephrol. 2011;22(5):849–58.
Mallela J, Yang J, Shariat-Madar Z. Prolylcarboxypeptidase: a cardioprotective enzyme. Int J Biochem Cell Biol. 2009;41(3):477–81.
Grand T, L’Hoste S, Mordasini D, Defontaine N, Keck M, Pennaforte T, et al. Heterogeneity in the processing of CLCN5 mutants related to Dent disease. Hum Mutat. 2011;32(4):476–83.
Monnens L, Levtchenko E. Evaluation of the proximal tubular function in hereditary renal Fanconi syndrome. Nephrol Dial Transplant. 2008;23(9):2719–22.
Servais A, Morinière V, Grünfeld JP, Noël LH, Goujon JM, Chadefaux-Vekemans B, et al. Late-onset nephropathic cystinosis: clinical presentation, outcome, and genotyping. Clin J Am Soc Nephrol. 2008;3(1):27–35.
Bisceglia L, Fischetti L, Bonis PD, Palumbo O, Augello B, Stanziale P, et al. Large rearrangements detected by MLPA, point mutations, and survey of the frequency of mutations within the SLC3A1 and SLC7A9 genes in a cohort of 172 cystinuric Italian patients. Mol Genet Metab. 2010;99(1):42–52.
Igarashi T, Inatomi J, Sekine T, Cha SH, Kanai Y, Kunimi M, et al. Mutations in SLC4A4 cause permanent isolated proximal renal tubular acidosis with ocular abnormalities. Nat Genet. 1999;23(3):264–6.
Satko SG, Sedor JR, Iyengar SK, Freedman BI. Familial clustering of chronic kidney disease. Semin Dial. 2007;20(3):229–36.
Böger CA, Heid IM. Chronic kidney disease: novel insights from genome-wide association studies. Kidney Blood Press Res. 2011;34(4):225–34.
Lang R, Lee G, Liu W, Tian S, Rafi H, Orias M, et al. KCNA10: a novel ion channel functionally related to both voltage-gated potassium and CNG cation channels. Am J Physiol Renal Physiol. 2000;278(6):F1013–21.
Yao X, Tian S, Chan HY, Biemesderfer D, Desir GV. Expression of KCNA10, a voltage-gated K channel, in glomerular endothelium and at the apical membrane of the renal proximal tubule. J Am Soc Nephrol. 2002;13(12):2831–9.
Köhler R, Ruth P. Endothelial dysfunction and blood pressure alterations in K + -channel transgenic mice. Pflugers Arch. 2010;459(6):969–76.
Simino J, Shi G, Arnett D, Broeckel U, Hunt SC, Rao DC. Variants on Chromosome 6p22.3 Associated with Blood Pressure in the HyperGEN Study: follow-up of FBPP Quantitative Trait Loci. Am J Hypertens. 2011;24(11):1227–33.
Senevirathna L, Abeysekera T, Nanayakkara S, Chandrajith R, Ratnatunga N, Harada KH, et al. Risk factors associated with disease progression and mortality in chronic kidney disease of uncertain etiology: a cohort study in Medawachchiya. Sri Lanka. Environ Health Prev Med. 2011;17(3):191–8.
Chen YM, Zhou Y, Go G, Marmerstein JT, Kikkawa Y, Miner JH. Laminin β2 gene missense mutation produces endoplasmic reticulum stress in podocytes. J Am SocNephrol. 2013;24(8):1223–33.
Chen YM, Kikkawa Y, Miner JH. A missense LAMB2 mutation causes congenital nephrotic syndrome by impairing laminin secretion. J Am Soc Nephrol. 2009;22(5):849–58.
The authors would like to acknowledge the Special Coordination Funds for Promoting Science and Technology from the Ministry of Education, Culture, Sports, Science and Technology in Japan.
Conflict of interests
Authors declare that they do not have any competing interests. The funding agency had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.
About this article
Cite this article
Nanayakkara, S., Senevirathna, S., Parahitiyawa, N.B. et al. Whole-exome sequencing reveals genetic variants associated with chronic kidney disease characterized by tubulointerstitial damages in North Central Region, Sri Lanka. Environ Health Prev Med 20, 354–359 (2015). https://doi.org/10.1007/s12199-015-0475-1