Research Article | Open Access
Novel Haplotype Indicator for End-Stage Renal Disease Progression among Saudi Patients
Background. End-stage renal disease (ESRD) is the result of hypertensive nephrosclerosis and chronic glomerular diseases and is associated with high morbidity and mortality. There are strong heritable components in the manifestation of the disease with a genetic predisposition to renal disorders, including focal segmental glomerulosclerosis and arterionephrosclerosis. Recent studies in genetics have examined modifiable risk factors that contribute to renal disease, and this has provided a deep insight into progressive kidney disease. Single-nucleotide polymorphisms at the proximity of SHROOM3, CST3, SLC7A9, and MYH9 genes have been associated with an increased risk of developing CKD and ESRD. Methods. A total of 160 CKD patients and 189 control subjects of Saudi origin participated in the study. Eight polymorphisms (SHROOM3-rs9992101, rs17319721; SLC7A9-rs4805834; MYH9-rs4821480, rs4821481, rs2032487, rs3752462; CST3-rs13038305) were genotyped using TaqMan assay, and the haplotype analysis was done using the HaploView 4.2 software. Results. Haplotype analysis revealed a novel haplotype “E6”-GTTT to be associated significantly with an increased risk for ESRD () and CKD (). Conclusion. CKD is often silent until symptomatic uremia during the advanced stages of the disease. The newly identified haplotype will help recognize patients at risk for a rapid progression of CKD to ESRD. Accurate detection and mapping of the genetic variants facilitates improved risk stratification and development of improved and targeted therapeutic management for CKD.
End-stage renal disease (ESRD) is commonly caused by hypertensive nephrosclerosis and chronic glomerular diseases  although a causative link between nephrosclerosis and hypertension is yet to be established . Studies have examined modifiable risk factors that contribute to CKD. However, there are strong heritable components in the manifestation of the disease [3, 4], with the vast majority of individuals suffering from comorbid conditions, such as hypertension or diabetes [5, 6].
Since the advent of genome-wide association studies (GWAS), many novel loci associated with common human diseases have been identified , including CKD [3, 8–10]. A genetic polymorphism in shroom family member 3 (SHROOM3) has been identified as a CKD susceptibility locus through GWAS. SHROOM3 is a regulator of epithelial cellular arrangement and planar remodeling , which contributes to glomerular filtration barrier integrity . Many studies have suggested that SHROOM3 plays an important role in mammalian kidney development and human kidney disease through estimated GFR (eGFR). One of these CKD-associated SHROOM3 variants, rs17319721, has been shown to impact cis-expression and renal allograft fibrosis . Genetic polymorphisms in solute carrier family 7-member 9 (SLC7A9) gene, an amino acid transporter in renal proximal tubule cells, cause cystinuria , showing an association with GFR [3, 15], and have been identified as a risk factor for CKD patients of European ancestry . Variants of cystatin C (CST3) have been also shown to impact altered eGFR and kidney disease . Polymorphisms in myosin heavy chain 9 (MYH9) gene on chromosome 22 have been shown to be associated with a risk for focal segmental glomerulosclerosis (FSGS), HIV-associated nephropathy CKD with admixed nondiabetic kidney disease, hypertension-associated ESRD, and nondiabetic etiologies of ESRD [17–20]. The MYH9 risk polymorphisms are common among African Americans, contributing approximately 40–45% of all ESRD and 70% of nondiabetic ESRD .
The United States Renal Data System 2016  reported the prevalence of ESRD to be 2,067 per million with an incidence rate of 370/million/year in the country. Genetic underpinnings of pediatric renal diseases, such as congenital and infantile nephrotic syndromes, are significantly higher in the Kingdom of Saudi Arabia (KSA) than in the Western world . The prevalence of ESRD in KSA has exhibited a rapid increase in the past decades resulting in a rate that exceeds those seen in European and American populations . The incidence and prevalence of CKD in the KSA is estimated to be approximately 1.72 million, equating to about 6% of the population. Out of these, only 7.1% are aware of their disease status, and this unawareness often results in poorer outcomes in such patients . Furthermore, there is a sharp annual increase in the rate of CKD patients who develop ESRD, and this accounts for 2.21% deaths annually . The prevalence of diabetic nephropathy among adult ESRD patients is 42.5% with a mortality rate of 18.6% compared to 6.9% of nondiabetic patients. While ethnicity is thought to play a large role in CKD genetics, very few genotyping studies of established CKD associations have been performed to date in Saudi Arabia and the surrounding regions. Our earlier study  highlighted the association of the eight SNPs in these four genetic regions and proposed a statistically significant method of predicting the CKD using FGF23, vitamin D3 level, and MYH9 genotypes . Haplotypes represent a combination of genetic determinants along a single chromosome that are either preserved intact or separated by recombination over time. Haplotypes are commonly used in research to identify a disease-conferring gene or locus. Presently, much interest surrounds the use of genetic association studies, which have an advisedly more powerful study design than the linkage- or family-based studies in localizing susceptibility loci for common diseases that confer moderate risk . Here, we present the haplotype association of these SNPs in 160 Saudi CKD and 189 non-CKD subjects from KSA towards the CKD progression risk assessment.
2. Materials and Methods
2.1. Study Population
The study included 160 Saudi CKD patients reporting to the Department of Nephrology at the King Fahd Hospital of the University, Al Khobar. Stringent inclusion criteria for both, cases and controls, were followed due to the high prevalence rate of CKD in the Saudi population. Patients with a history of any phosphate wasting disorder, including tumor-induced osteomalacia (TIO), X-linked hypophosphatemia (XLH), autosomal dominant hypophosphatemic rickets (ADHR), and untreated primary hyperthyroidism, and those undergoing renal replacement therapy were excluded from the study. Controls consisted of 189 healthy subjects of Saudi origin without evidence of renal disorders (serum creatinine <1.4 and <1.2 mg/dl in men and women, respectively). This study was approved by the Ethical Committee of the Imam Abdulrahman Bin Faisal University and conducted according to the Declaration of Helsinki. Signed written informed consent was obtained from all participants.
2.2. Genotyping Assay
DNA was isolated from the blood of 349 individuals using QIAamp DNA Mini Kit (Qiagen, USA) as per the manufacturer’s instructions. Genotyping analysis was performed by ABI TaqMan SNP genotyping assays on 160 CKD patients and 189 controls. Allele-specific TaqMan® PCR technology is a highly sensitive assay and was performed using ABI 7500 Fast real-time PCR. Disease-associated SNPs selected from published articles, namely, rs9992101, rs17319721 (SHROOM3), rs4805834 (SLC7A9), rs4821480, rs4821481, rs2032487, rs3752462 (MYH9), and rs13038305 (CST3), were included. An additional 100 coronary artery disease patients with hypertension were genotyped for the MYH9 (rs4821480, rs4821481, rs2032487, and rs3752462) SNPs.
2.3. E Haplotyping
MYH9 haplotypes spanning 12–23 introns were reconstructed by 4 tagging SNPs. The haplotypes comprised the four SNP loci of the MYH9 gene in the order rs4821480 (T/g), rs4821481 (T/c), rs2032487 (C/t), and rs3752462 (C/t). The 4 SNPs were genotyped, and the haplotype blocks were determined using the HaploView 4.2 software that assigned haplotypes into chromosome-specific blocks based on the partition-ligation approach through EM algorithm . To identify the nonrandom association of the eight SNPs by estimates of linkage disequilibrium (LD), each pair of SNPs was computed using the standard D-prime method. Patient population was further stratified into cohorts as CKD (eGFR > 15) and ESRD (eGFR < 15) to assess the significance of the haplotype variations.
2.4. Statistical Analysis
SPSS version 19 (Chicago, Illinois) was used to complete the statistical analysis. Mutation status was classified as positive or negative qualitatively. The genotype frequencies were tested for Hardy–Weinberg equilibrium among the control subjects.
The present study investigated 160 CKD cases, in which 85 were males and 75 were females, with a mean age of 47.68 ± 17.27 years. The baseline characteristics of the study participants are shown in Supplementary Table 1.
TaqMan analysis for genetic variants rs9992101, rs17319721 (SHROOM3), rs4805834 (SLC7A9), rs4821480, rs4821481, rs2032487, and rs3752462 (MYH9), and rs13038305 (CST3) was carried out for all samples. The control group was consistent with Hardy–Weinberg equilibrium for all SNPs. The distribution of the analyzed genotype polymorphisms was reported in our previous study .
The HaploView analysis results are represented in Figure 1. Tables 1 and 2 indicate the frequency of various MYH9 E and SHROOM3 haplotypes and their association, respectively. A new haplotype “E6”-GTTT is revealed to be significantly associated with the patient cohort (), stratified CKD cohort (), and ESRD cohort (). The common haplotype GTTC (E4) () was also found to be associated with CKD (Table 1). On further stratification, E6 haplotype was found to be strongly associated with ESRD () in ESRD vs. control and in CKD vs. ESRD () analysis.
To validate the new E6 haplotype, an additional 100 CAD patients with hypertension were genotyped for the MYH9 (rs4821480, rs4821481, rs2032487, and rs3752462) SNPs. Haplotype analysis revealed that the new haplotype was insignificant in the CAD patients. Also, the novel E6 haplotype was still found to be significantly associated () when ESRD was compared to the CAD cohort (Table 1).
CKD is becoming an important health issue worldwide and a major cause for morbidity and mortality. Genetic searches for vital markers of CKD are essential to identify individuals who are at risk for ESRD. Numerous genes have been shown to be associated with CKD. SHROOM3, a regulator of epithelial cellular arrangement contributes to glomerular filtration barrier integrity while, CST3 which is involved in creatinine and cystatin synthesis, is strongly associated with multiple kidney-related traits. The SLC7A9 gene, which is expressed in renal proximal tubule cells, was shown to be strongly associated with the markers of kidney function, creatinine, and eGFR. MYH9 gene encodes a protein which is expressed in the glomerular podocyte .
MYH9 gene mutations result in nephritis to varying degrees in Epstein and Fechtner syndromes . MYH9 SNPs and haplotypes are associated with a risk for T2D and non-T2D nephropathy, lupus nephritis, hypertensive nephropathy, and FSGS [17, 18, 20, 29–32]. The haplotype analysis of this study shows more combinatorial importance for the MYH9 gene region than SNPs. Oleksyk et al.  reported that the E1 haplotype of the MYH9 gene region, prominent in sub-Saharan Africa, is probably involved in the increased risk of developing CKD by increasing glomerulosclerosis and proteinuria through activation of nephritis by the deregulation of podocyte function and not by immunological mechanisms. In the present study, the E1 (GCCT) haplotype lacked any association towards CKD or ESRD, contradictory to the findings of Colares et al. . The rare E5 (GCTC) is the only other haplotype carrying the C allele at rs4821481, which was noted in only two of the human genome diversity project populations, with frequencies of 0.02 in Mandenka, a West African ethnic group, and 0.01 in Palestinians . Kopp et al.  reported the E2 haplotype prominent in the European, Middle Eastern, and South and Central Asian populations to be protective against renal disease, and Tavira et al.  reported the same protective effect in a Spanish cohort. The E2 (TTTC) and E3 (TTTT) haplotypes lacked an association in the present study. Interestingly, the haplotype frequencies of the present study were divergent from those reported from other populations . The E4 and E6 haplotypes were found to be significantly associated with an increased risk of CKD and ESRD, respectively. Even in the stratified groups, ESRD (eGFR < 15 mL/min/1.73 m2) and CKD (eGFR > 15 mL/min/1.73 m2), the E4 (GTTC) was associated with an increased risk of CKD, although the novel haplotype GTTT () was found to be only associated with ESRD.
The SNPs rs9992101 and rs17319721 located in the SHROOM3 gene on chromosome 4q21 are closely associated with CKD , and the former is in high linkage disequilibrium (LD) with rs17319721. The rs17319721 (A) allele is associated with increased SHROOM3 transcription and is therefore associated with an increased glomerular filtration rate, thereby increasing the risk for CKD . Both these SNPs and their haplotypes were not associated with a risk of CKD in the present study.
CKD is often silent until the advanced stages of the disorder. Thus, many individuals remain unaware until symptomatic uremia is detected, and they begin to suffer from renal ailments. The newly identified haplotype will help identify the patients at risk for quicker progression of CKD towards ESRD. This study will further our understanding of the biological mechanisms of kidney function by identifying loci which may potentially influence metabolic renal functions.
The datasets used and/or analyzed during the current study are available from the corresponding author.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
We thank the Deanship of Scientific Research of the University of Dammam (to Dr. Cyril Cyrus, Grant no. 2014115 and IRB # IRB-2014-08–046) and King Abdulaziz City for Science and Technology (KACST) (Grant no. 12-MED2799-46) for the funding. We also thank Mr. Geoffrey James Tam Moro and Mr. Mohammed Hassan Alshamlan for their technical and administrative support.
Supplementary Table 1: baseline characteristics of the total study population and stagewise distribution. (Supplementary Materials)
- US Renal Data System, “USRDS 2006 Annual Data Report: Atlas of End-Stage Renal Disease in the United States,” National Institutes of Health, National Institute of Diabetes and Digestive and Kidney Diseases, Bethesda, MD, USA, 2006.
- B. I. Freedman, S. S. Iskandar, and R. G. Appel, “The link between hypertension and nephrosclerosis,” American Journal of Kidney Diseases, vol. 25, no. 2, pp. 207–221, 1995.
- A. Köttgen, C. Pattaro, C. A. Böger et al., “Multiple new loci associated with kidney function and chronic kidney disease,” Nature Genetics, vol. 42, no. 5, pp. 376–384, 2010.
- V. M. Vehaskari, “Genetics and CKD,” Advances in Chronic Kidney Disease, vol. 18, no. 5, pp. 317–323, 2011.
- A. Ingsathit, A. Thakkinstian, A. Chaiprasert et al., “Prevalence and risk factors of chronic kidney disease in the Thai adult population: Thai SEEK study,” Nephrology Dialysis Transplantation, vol. 25, no. 5, pp. 1567–1575, 2010.
- R. A. Nugent, S. F. Fathima, A. B. Feigl, and D. Chyung, “The burden of chronic kidney disease on developing nations: a 21st century challenge in global health,” Nephron Clinical Practice, vol. 118, no. 3, pp. c269–c277, 2011.
- J. MacArthur, E. Bowler, M. Cerezo et al., “The new NHGRI-EBI catalog of published genome-wide association studies (GWAS catalog),” Nucleic Acids Research, vol. 45, no. D1, pp. D896–D901, 2017.
- J. C. Chambers, W. Zhang, G. M. Lord et al., “Genetic loci influencing kidney function and chronic kidney disease,” Nature Genetics, vol. 42, no. 5, pp. 373–375, 2010.
- M. Gorski, A. Tin, M. Garnaas et al., “Genome-wide association study of kidney function decline in individuals of European descent,” Kidney International, vol. 87, no. 5, pp. 1017–1029, 2015.
- C. Pattaro, A. Teumer, M. Gorski et al., “Genetic associations at 53 loci highlight cell types and biological pathways relevant for kidney function,” Nature Communications, vol. 7, no. 1, Article ID 10023, 2016.
- T. Nishimura and M. Takeichi, “Shroom3-mediated recruitment of Rho kinases to the apical cell junctions regulates epithelial and neuroepithelial planar remodeling,” Development, vol. 135, no. 8, pp. 1493–1502, 2008.
- N. C. Yeo, C. C. O’Meara, J. A. Bonomo et al., “Shroom3 contributes to the maintenance of the glomerular filtration barrier integrity,” Genome Research, vol. 25, no. 1, pp. 57–65, 2015.
- M. C. Menon, P. Y. Chuang, Z. Li et al., “Intronic locus determines SHROOM3 expression and potentiates renal allograft fibrosis,” Journal of Clinical Investigation, vol. 125, no. 1, pp. 208–221, 2015.
- A. Mattoo and D. S. Goldfarb, “Cystinuria,” Seminars in Nephrology, vol. 28, no. 2, pp. 181–191, 2008.
- C. A. Böger, M. Gorski, M. Li et al., “Association of eGFR-related loci identified by GWAS with incident CKD and ESRD,” PLoS Genetics, vol. 7, no. 9, Article ID e1002292, 2011.
- C. M. O’Seaghdha, A. Tin, Q. Yang et al., “Association of a cystatin C gene variant with cystatin C levels, CKD, and risk of incident cardiovascular disease and mortality,” American Journal of Kidney Diseases, vol. 63, no. 1, pp. 16–22, 2014.
- W. H. L. Kao, M. J. Klag, L. A. Meoni et al., “MYH9 is associated with nondiabetic end-stage renal disease in African Americans,” Nature Genetics, vol. 40, no. 10, pp. 1185–1192, 2008.
- J. B. Kopp, M. W. Smith, G. W. Nelson et al., “MYH9 is a major-effect risk gene for focal segmental glomerulosclerosis,” Nature Genetics, vol. 40, no. 10, pp. 1175–1184, 2008.
- B. I. Freedman, P. J. Hicks, M. A. Bostrom et al., “Polymorphisms in the non-muscle myosin heavy chain 9 gene (MYH9) are strongly associated with end-stage renal disease historically attributed to hypertension in African Americans,” Kidney International, vol. 75, no. 7, pp. 736–745, 2009.
- D. M. Behar, S. Rosset, S. Tzur et al., “African ancestry allelic variation at the MYH9 gene contributes to increased susceptibility to non-diabetic end-stage kidney disease in Hispanic Americans,” Human Molecular Genetics, vol. 19, no. 9, pp. 1816–1827, 2010.
- M. A. Bostrom and B. I. Freedman, “The spectrum of MYH9-associated nephropathy,” Clinical Journal of the American Society of Nephrology, vol. 5, no. 6, pp. 1107–1113, 2010.
- R. Saran, B. Robinson, K. C. Abbott et al., “United States Renal Data System. 2016 USRDS annual data report: epidemiology of kidney disease in the United States,” American Journal of Kidney Diseases, vol. 69, no. 3, p. A4, 2017.
- J. A. Kari, “Pediatric renal diseases in the kingdom of Saudi Arabia,” World Journal of Pediatrics, vol. 8, no. 3, pp. 217–221, 2012.
- A. A. Al-Sayyari and F. A. Shaheen, “End stage chronic kidney disease in Saudi Arabia. A rapidly changing scene,” Saudi Medical Journal, vol. 32, no. 4, pp. 339–346, 2011.
- A. O. Alsuwaida, Y. M. Farag, A. A. Al Sayyari et al., “Epidemiology of chronic kidney disease in the Kingdom of Saudi Arabia (SEEK-Saudi investigators)—a pilot study,” Saudi Journal of Kidney Diseases and Transplantation, vol. 21, no. 6, pp. 1066–1072, 2010.
- C. Cyrus, S. Al-Mueilo, C. Vatte et al., “Assessing known chronic kidney disease associated genetic variants in Saudi Arabian populations,” BMC Nephrology, vol. 19, no. 1, p. 88, 2018.
- D. C. Crawford and D. A. Nickerson, “Definition and clinical importance of haplotypes,” Annual Review of Medicine, vol. 56, no. 1, pp. 303–320, 2005.
- J. C. Barrett, B. Fry, J. Maller, and M. J. Daly, “Haploview: analysis and visualization of LD and haplotype maps,” Bioinformatics, vol. 21, no. 2, pp. 263–265, 2005.
- B. I. Freedman, J. B. Kopp, C. A. Winkler et al., “Polymorphisms in the nonmuscle myosin heavy chain 9 gene (MYH9) are associated with albuminuria in hypertensive African Americans: the HyperGEN study,” American Journal of Nephrology, vol. 29, no. 6, pp. 626–632, 2009.
- C. M. O’Seaghdha, R. S. Parekh, S.-J. Hwang et al., “The MYH9/APOL1 region and chronic kidney disease in European-Americans,” Human Molecular Genetics, vol. 20, no. 12, pp. 2450–2456, 2011.
- J. N. Cooke, M. A. Bostrom, P. J. Hicks et al., “Polymorphisms in MYH9 are associated with diabetic nephropathy in European Americans,” Nephrology Dialysis Transplantation, vol. 27, no. 4, pp. 1505–1511, 2012.
- N. Franceschini, V. S. Voruganti, K. Haack et al., “The association of the MYH9 gene and kidney outcomes in American Indians: the strong heart family study,” Human Genetics, vol. 127, no. 3, pp. 295–301, 2010.
- T. K. Oleksyk, G. W. Nelson, P. An, J. B. Kopp, and C. A. Winkler, “Worldwide distribution of the MYH9 kidney disease susceptibility alleles and haplotypes: evidence of historical selection in Africa,” PLoS One, vol. 5, no. 7, Article ID e11474, 2010.
- V. S. Colares, S. M. de Oliveira Titan, A. da Costa Pereira et al., “MYH9 and APOL1 gene polymorphisms and the risk of CKD in patients with lupus nephritis from an admixture population,” PLoS One, vol. 9, no. 3, Article ID e87716, 2014.
- B. Tavira, E. Coto, J. Gómez et al., “Association between a MYH9 polymorphism (rs3752462) and renal function in the Spanish RENASTUR cohort,” Gene, vol. 520, no. 1, pp. 73–76, 2013.
Copyright © 2019 Cyril Cyrus 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.