Haplotype association analysis of genes within the WNT signalling pathways in diabetic nephropathy
BMC Nephrology volume 14, Article number: 126 (2013)
Renal interstitial fibrosis and glomerular sclerosis are hallmarks of diabetic nephropathy (DN) and several studies have implicated members of the WNT pathways in these pathological processes. This study comprehensively examined common genetic variation within the WNT pathway for association with DN.
Genes within the WNT pathways were selected on the basis of nominal significance and consistent direction of effect in the GENIE meta-analysis dataset. Common SNPs and common haplotypes were examined within the selected WNT pathway genes in a white population with type 1 diabetes, discordant for DN (cases: n = 718; controls: n = 749). SNPs were genotyped using Sequenom or Taqman assays. Association analyses were performed using PLINK, to compare allele and haplotype frequencies in cases and controls. Correction for multiple testing was performed by either permutation testing or using false discovery rate.
A logistic regression model including collection centre, duration of diabetes, and average HbA1c as covariates highlighted three SNPs in GSK3B (rs17810235, rs17471, rs334543), two in DAAM1 (rs1253192, rs1252906) and one in NFAT5 (rs17297207) as being significantly (P < 0.05) associated with DN, however these SNPs did not remain significant after correction for multiple testing. Logistic regression of haplotypes, with ESRD as the outcome, and pairwise interaction analyses did not yield any significant results after correction for multiple testing.
These results indicate that both common SNPs and common haplotypes of WNT pathway genes are not strongly associated with DN. However, this does not completely exclude these or the WNT pathways from association with DN, as unidentified rare genetic or copy number variants could still contribute towards the genetic architecture of DN.
Diabetic nephropathy (DN) is a microvascular complication of diabetes and is the most frequent cause of end-stage renal disease (ESRD) in western populations . Renal interstitial fibrosis and glomerular sclerosis are histological hallmarks of DN and several studies have implicated members of the WNT pathways in these pathological processes [2–9].
The WNT pathways can be subdivided into canonical β-catenin dependent and non-canonical β-catenin independent pathways (Figure 1). Canonical WNT signalling is central to numerous developmental processes and variants discovered within members of this pathway have been implicated in multiple complex diseases such as familial adenomatous polyposis coli, colorectal and hepatocellular cancers, type 2 diabetes and schizophrenia [10–14]. Non-canonical WNT signalling is less well defined, in part due to further subdivisions into the WNT/Ca2+ and the WNT planar cell polarity pathways. These pathways have been shown to modulate cytoskeletal reorganisation and activation of the JNK and MAPK signalling pathways [15–17], both of which can affect the motility and adherence of the mesangial cell, perturbing the cell's response to dynamic mechanical forces, which is the key function of the mesangial cell.
In vitro epithelial-to-mesenchymal transition (EMT) promotes renal fibrosis and can be induced by TGF-β1 , an integrin-linked kinase. Both the canonical WNT pathway and TGF-β1 require activation of β-catenin, in addition to the E-cadherin/β-catenin complexes localised to epithelial intercellular junctions, implicating both β-catenin and the WNT pathway in the regulation of the EMT . Furthermore, GSK3B, the protein responsible for phosphorylation of the β-catenin molecule and its subsequent proteosomal degradation, has been shown to prevent transition to a mesenchymal phenotype in human embryonic stem cells . Several WNT ligands, FZD receptors and β-catenin have all been reported to be differentially expressed in the unilateral ureteral obstructed (UUO) mouse model of renal injury . In addition, Dickkopf-1 (DKK1), a WNT signalling antagonist, was shown to promote hyperglycaemia-induced mesangial matrix expansion in rat mesangial cells .
Previously, common variants within four key genes in the WNT pathway have been investigated for association with diabetic nephropathy . In the present study, a more comprehensive assessment was undertaken of common variants in multiple genes within the WNT pathway. Due to the large number of WNT pathway genes (>65), eight potential candidate genes were chosen on the basis of single nucleotide polymorphisms (SNPs) reaching a nominal significance threshold of 0.05 from the meta-analysed Genetics of Nephropathy–an International Effort (GENIE) Consortium dataset (Table 1) . The chosen SNPs also showed a consistent direction of effect in each of the three case–control collections represented by the GENIE Consortium meta-analysed dataset, an international collaboration of three cohorts of type 1 diabetic patients discordant for DN totalling 2916 with nephropathy and 3315 without nephropathy . Three additional genes, CTNNB1, WNT5A and WNT6, were also included within the analysis despite failing to meet the inclusion criteria, on the basis of previous suggestion of their involvement in the pathogenesis of DN. Although the genotyping platforms used to determine the GENIE data provided reasonable coverage across the potential genes of interest, additional informative haplotype tagging SNPs identified through CEU participant data from HapMap offers a more comprehensive evaluation of any potential genetic effect.
Research ethics approval was obtained from the South and West Multicentre Research Ethics Committee (MREC/98/6/71) and Queens University Belfast Research Ethics Committee, and written informed consent obtained prior to participation. All recruited individuals were white, had type 1 diabetes mellitus (T1D) diagnosed before 32 years of age and were born in the UK or Ireland. Cases with nephropathy (n = 718) and controls without nephropathy (n = 749) were from the Warren 3/UK Genetics of Kidneys in Diabetes (GoKinD) and all-Ireland collections . The definition of DN in cases was based on development of persistent proteinuria (>0.5g protein/24h) at least 10 years after diagnosis of T1D, hypertension (blood pressure > 135/85 mmHg or treatment with antihypertensive agents) and associated diabetic retinopathy. Controls were individuals with T1D for at least 15 years with normal urinary albumin excretion rates and no evidence of microalbuminuria on repeated testing (at least 3 assays measuring albumin excretion over a minimum 12 month period, with each test separated by at least 3 months). In addition, control subjects had not been prescribed antihypertensive drug treatment avoiding possible misclassification of diabetic individuals with nephropathy as ‘control phenotypes’ when the use of antihypertensive treatment may have reduced urinary albumin excretion into the normal range. Individuals with micro-albuminuria were excluded from both case and control groups since it is not possible to confidently assign a case or control status to such individuals as their urinary albumin excretion may either regress or progress over time .
Haplotype definition, SNP selection and genotyping
A total of 11 genes (AXIN1, CALM3, CTNNB1, DAAM1, DKK2, GSK3B, NFAT5, WNT3, WNT5A, WNT6 and WNT16) were chosen for genotyping (Table 1). SNPs were selected from within these 11 genes to tag common haplotypes (>5% frequency within the HapMap CEU population). Haplotypes for each gene investigated were selected from Phase III, release 2 HapMap (http://www.hapmap.org) CEPH data (Utah residents with ancestry in northern and western Europe; CEU) using Haploview (http://www.broadinstitute.org/haploview) to visualise common haplotypes. Haplotypes were defined using the confidence interval method in Haploview as described in Gabriel et al. . Adjacent haplotypes that had a multi-allelic D-prime of greater 0.9 were combined in an iterative fashion. SNPs were selected using multi-marker tagging for their ability to tag unique haplotypes with r 2 > 0.8 (LOD threshold 3.0). All SNPs had a minor allele frequency (MAF) ≥5%, with quality control filters of genotype call rate ≥95%, and no deviation from Hardy–Weinberg equilibrium (HWE; P < 0.001).
Genotyping was performed by MassARRAY iPLEX (Sequenom, San Diego, CA, USA) or Taqman 5' nuclease (Applied Biosystems, Foster City, CA, USA) assays according to the manufacturers’ instructions. DNA samples were excluded if missing genotypes exceeded 10%. Other quality control measures included parent/offspring trio samples, duplicates on plates, random sample allocation to plates, independent scoring of problematic genotypes by two individuals and re-sequencing of selected DNAs to validate genotypes.
Clinical characteristics of cases and controls were compared using the z-test for large independent samples and the χ2 test. Association analyses were performed using PLINK . Initially a χ2 test for trend (1 df) was used with adjustment for collection centre. Logistic regression analysis was then performed on each SNP with terms for potential confounders (collection centre, gender, duration of T1D and HbA1c) included in the model. The level of statistical significance was set at 5% with correction for multiple testing performed by permutation test (n = 100,000). Pairwise interactions between SNPs were tested in the statistical programming package R, using logistic regression to compare models with and without the interaction terms to obtain a likelihood ratio test. The results of the interaction analysis were corrected for multiple testing by false discovery rate (FDR < 5%).
Results and discussion
A total of 90 SNPs were genotyped, 85 using MassARRAY iPLEX Gold technology (Sequenom, San Diego, CA, USA), and five using Taqman 5’ nuclease assay (Applied Biosystems, Foster City, CA, USA) in 719 cases and 748 controls. Quality criteria were applied to the data before association analysis. A total of 35 individuals with more than 10% missing genotype data were removed from the analysis. All SNPs passed the genotyping and Hardy-Weinberg thresholds of 95% and P < 0.001 respectively. No Mendelian errors or inconsistencies between duplicate samples were observed. The final average genotyping rate was 98.9% in 700 cases, and 732 controls.
The clinical characteristics of the DN cases (n = 700) and diabetic controls (n = 732) genotyped in this study, which met quality control filters, are listed in Table 2. There were more males, higher mean HbA1c and blood pressure values (despite the use of antihypertensive treatment) in the case group compared with the control group. All comparisons were significant at P < 0.001 with the exception of age at diagnosis which did not differ significantly between groups. Approximately one quarter of cases (26.6%) had ESRD.
SNPs chosen to tag common haplotypes across the 11 genes selected on the basis of their significant and common direction of effect across the GENIE cohorts (Table 1)  were assessed by logistic regression analysis with adjustment for collection centre, gender, duration of T1D and HbA1c (Table 3). Twenty-six putative linkage disequilibrium (LD) blocks were identified across the 11 genes, yielding 110 common haplotypes with an estimated frequency >5%. None of the haplotypes examined were significantly associated with DN at P < 0.01, however eight haplotypes were significantly associated with DN at P < 0.05. Of the eight haplotypes, three were in GSK3B, two in AXIN1, two in DAAM1, and one in NFAT5. However, no significant association between haplotype and DN remained after correction for multiple testing (data not shown).
In a single marker analysis, adjusted by collection centre, no SNPs were associated with DN at P < 0.01 (the significance threshold set was corrected for multiple testing), however five SNPs, rs17810235, rs11639947, rs11646942, rs17095819, and rs17510191 in GSK3B, NFAT5, AXIN1, DAAM1, DKK2 had P-values <0.05 as shown in Table 4a. Logistic regression analyses were performed with adjustment for collection centre, gender, duration of T1D, and average HbA1c as covariates in the model. The most significant association was reported for rs17810235 in GSK3B (Table 4b; P = 0.006). Five additional SNPs demonstrated a P <0.05, although they were not supported in the univariate analysis alone. Although limited in power, a subgroup analysis defined by comparison of ESRD as the primary phenotype versus non-ESRD, identified two significantly associated SNPs, rs1253192 and rs11079737 in DAAM1 and WNT3 respectively with P = 0.009, although neither association survived correction for multiple testing (Table 4c).
Assessment of 4005 pair-wise interactions between the 90 SNPs was performed by logistic regression analyses with adjustment for collection centre, gender, duration of T1D, and HbA1c. The χ2 likelihood ratio tests (Table 5) identified four interaction terms as nominally significant between SNPs in AXIN1, DAAM1, DKK2, WNT3 and WNT6. However, these interactions were not significant following correction for multiple testing by FDR P < 5%.
There is an increasing body of evidence to suggest that modulation of the WNT pathways may play a role in the development of DN. β-catenin and TCF/LEF have been shown to directly induce the expression of the cyclic nucleotide phosphodiesterase CNP, a regulator of fibroblast proliferation and activation . The expression of secreted frizzled-related protein 4 (SFRP4), an inhibitor of WNT signalling, is decreased following UUO. This decrease is concomitant with increased levels of WNT/β-catenin signalling, in tubular and interstitial cells, along with increased fibronectin and smooth muscle actin, both markers of fibrosis. Introduction of recombinant SFRP4 reduced the markers of fibrosis and WNT/β-catenin signalling. Furthermore E-cadherin expression was partially maintained by treatment with recombinant SFRP4, and the number of myofibroblasts decreased . DKK1 is shown to be increased in mesangial cells in response to stimulation with high concentrations of glucose . In addition high concentrations of glucose decreased WNT signalling and increased TGF-β1 and fibronectin expression in mesangial cells. Transfection of WNT4, WNT5a, GSK3β and β-catenin ameliorated the TGF-β1-induced fibrosis . Cultured podocytes with stabilised β-catenin are less motile and less adherent to the extracellular matrix whereas deletion of β-catenin rendered the cells more susceptible to apoptosis .
Gene-based assessments of association are increasingly been viewed as a useful complement to genome-wide association studies (GWAS) . The gene-based approach reduces the problems associated with multiple testing that inhibit GWAS by reducing the number of statistical tests under consideration. Our study has adopted a two stage approach to evaluate common variants in all WNT pathway members in relation to DN. SNPs located in genes implicated in the WNT pathways that failed to demonstrate significant association and direction of effect across all GENIE cohorts  were excluded at the first step. WNT pathway members that demonstrated significant association and direction of effect with DN across the three GENIE case–control collections were then evaluated more meticulously through refined genotyping of haplotype tagging SNPs. This approach offers a more comprehensive assessment of common variants across the WNT pathways in comparison to previously published studies. Univariate SNP analysis failed to identify any association with DN. Multivariate regression analyses of common haplotypic structure also failed to reveal any associations that remained significant after correction for multiple testing. All possible combinations of pair-wise SNP-SNP interactions were tested as an interaction term in a logistic regression model. Due to the large number of tests, and the unsuitability of permutations as a correction for multiple testing in interaction analyses, the false discovery rate method was used, although no associations remained significant after correction.
There are a number of inherent limitations associated with using a restricted number of SNPs across a chosen set of genes : (1) identification of association does not necessarily equate to functional significance given the concept of LD. (2) assessing one or two SNPs per gene may provide inadequate representation of the genetic architecture at that locus. (3) patterns of LD can vary significantly within and between different populations and therefore a significant association in one population may not necessarily translate across all populations. In addition, common genetic loci are likely to explain only a proportion of the variation contributing to the phenotype under consideration. Evidence in support of rare variants with potentially large individual effect size is currently under investigation in DN. Since our study focused only on common variants, untyped, highly penetrant rare variants in these genes could also contribute to DN. The size of our study was such that it had 90% power to detect variants with odds ratios of 1.69, 1.50, 1.44, and 1.42 and MAF of 10%, 20%, 30%, and 40%, respectively. Sample sizes of the magnitude required to detect variants present in the population with a lower frequency or with a smaller effect size were not available for analysis and would require international, multicentre collaborative efforts. Future amalgamation of independent cohorts with similar DN phenotypes will enable a more robust evaluation of such loci. In addition, other factors such as copy number variation or indeed epigenetic mechanisms (e.g. DNA methylation, histone modification and microRNAs) may also modify gene function and/or expression profiles affecting these pathways and modulating disease risk accordingly.
There was no association observed between DN and either common variants or haplotypes among any of the genes associated with the WNT pathway. It is unlikely that common variants located within WNT pathway genes have a major role in the underlying genetics of DN. Further investigation of rare variants, copy number variants or epigenetic mechanisms, in WNT pathway members may identify potential risk factors that contribute to the genetic susceptibility of DN, in addition to identifying potential therapeutic targets for this disease.
Gilg J, Castledine C, Fogarty D: Chapter 1 UK RRT incidence in 2010: national and centre-specific analyses. Nephron Clin Pract. 2012, 120 (Suppl 1): c1-27.
Surendran K, Simon TC: CNP gene expression is activated by Wnt signaling and correlates with Wnt4 expression during renal injury. Am J Physiol Renal Physiol. 2003, 284: F653-662.
Surendran K, Schiavi S, Hruska KA: Wnt-dependent beta-catenin signaling is activated after unilateral ureteral obstruction, and recombinant secreted frizzled-related protein 4 alters the progression of renal fibrosis. J Am Soc Nephrol. 2005, 16: 2373-2384. 10.1681/ASN.2004110949.
He W, Dai C, Li Y, Zeng G, Monga SP, Liu Y: Wnt/beta-catenin signaling promotes renal interstitial fibrosis. J Am Soc Nephrol. 2009, 20: 765-776. 10.1681/ASN.2008060566.
Lin CL, Wang JY, Ko JY, Huang YT, Kuo YH, Wang FS: Dickkopf-1 promotes hyperglycemia-induced accumulation of mesangial matrix and renal dysfunction. J Am Soc Nephrol. 2010, 21: 124-135. 10.1681/ASN.2008101059.
Kato H, Gruenwald A, Suh JH, et al: Wnt/β-catenin pathway in podocytes integrates cell adhesion, differentiation, and survival. J Biol Chem. 2011, 286: 26003-26015. 10.1074/jbc.M111.223164.
Rooney B, O’Donovan H, Gaffney A, et al: CTGF/CCN2 activates canonical Wnt signalling in mesangial cells through LRP6: implications for the pathogenesis of diabetic nephropathy. FEBS Lett. 2011, 585: 531-538. 10.1016/j.febslet.2011.01.004.
Ho C, Lee PH, Hsu YC, Wang FS, Huang YT, Lin CL: Sustained Wnt/β-Catenin Signaling Rescues High Glucose Induction of Transforming Growth Factor-β1-Mediated Renal Fibrosis. Am J Med Sci. 2012, 344: 374-382.
Zhou T, He X, Cheng R, et al: Implication of dysregulation of the canonical wingless-type MMTV integration site (WNT) pathway in diabetic nephropathy. Diabetologia. 2012, 55: 255-266. 10.1007/s00125-011-2314-2.
Hart M, Concordet JP, Lassot I, et al: The F-box protein beta-TrCP associates with phosphorylated beta-catenin and regulates its activity in the cell. Curr Biol. 1999, 9: 207-210. 10.1016/S0960-9822(99)80091-8.
Polakis P: Wnt signaling and cancer. Genes Dev. 2000, 14: 1837-1851.
Grant SFA, Thorleifsson G, Reynisdottir I, et al: Variant of transcription factor 7-like 2 (TCF7L2) gene confers risk of type 2 diabetes. Nat Genet. 2006, 38: 320-323. 10.1038/ng1732.
Saxena R, Gianniny L, Burtt NP, et al: Common single nucleotide polymorphisms in TCF7L2 are reproducibly associated with type 2 diabetes and reduce the insulin response to glucose in nondiabetic individuals. Diabetes. 2006, 55: 2890-2895. 10.2337/db06-0381.
Steinberg S, de Jong S, Andreassen OA, et al: Common variants at VRK2 and TCF4 conferring risk of schizophrenia. Human Mol Genet. 2011, 20: 4076-4081. 10.1093/hmg/ddr325.
Habas R, Dawid IB, He X: Coactivation of Rac and Rho by Wnt/Frizzled signaling is required for vertebrate gastrulation. Gen Dev. 2003, 17: 295-309. 10.1101/gad.1022203.
Veeman MT, Axelrod JD, Moon RT: A second canon. Functions and mechanisms of beta-catenin-independent Wnt signaling. Dev Cell. 2003, 5: 367-377. 10.1016/S1534-5807(03)00266-1.
Habas R, Dawid IB: Dishevelled and Wnt signaling: is the nucleus the final frontier?. J Biol. 2005, 4: 2-10.1186/jbiol22.
Liu Y: New insights into epithelial-mesenchymal transition in kidney fibrosis. J Am Soc Nephrol. 2010, 21: 212-222. 10.1681/ASN.2008121226.
Ivanova L, Butt MJ, Matsell DG: Mesenchymal transition in kidney collecting duct epithelial cells. Am J Physiol Renal Physiol. 2008, 294: F1238-1248. 10.1152/ajprenal.00326.2007.
Ullmann U, Gilles C, De Rycke M, Van de Velde H, Sermon K, Liebaers I: GSK-3-specific inhibitor-supplemented hESC medium prevents the epithelial-mesenchymal transition process and the up-regulation of matrix metalloproteinases in hESCs cultured in feeder-free conditions. Mol Hum Reprod. 2008, 14: 169-179. 10.1093/molehr/gan001.
Kavanagh DH, Savage DA, Patterson CC, et al: Association analysis of canonical Wnt signalling genes in diabetic nephropathy. PLoS One. 2011, 6: e23904-10.1371/journal.pone.0023904.
Sandholm N, Salem RM, McKnight AJ, et al: New Susceptibility Loci Associated with Kidney Disease in Type 1 Diabetes. PLoS Genet. 2012, 8: e1002921-10.1371/journal.pgen.1002921.
Kavanagh D, McKay GJ, Patterson CC, McKnight AJ, Maxwell AP, Savage DA: Association analysis of Notch pathway signalling genes in diabetic nephropathy. Diabetologia. 2011, 54: 334-338. 10.1007/s00125-010-1978-3.
Perkins BA, Ficociello LH, Silva KH, Finkelstein DM, Warram JH, Krolewski AS: Regression of microalbuminuria in type 1 diabetes. N Engl J Med. 2003, 348: 2285-2293. 10.1056/NEJMoa021835.
Gabriel SB, Schaffner SF, Nguyen H, et al: The structure of haplotype blocks in the human genome. Science. 2002, 296: 2225-2229. 10.1126/science.1069424.
Purcell S, Neale B, Todd-Brown K, et al: PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Human Genet. 2007, 81: 559-575. 10.1086/519795.
Liu JZ, McRae AF, Nyholt DR, et al: A versatile gene-based test for genome-wide association studies. Am J Hum Genet. 2010, 87: 139-145. 10.1016/j.ajhg.2010.06.009.
Cardon LR, Bell JI: Association study designs for complex diseases. Nat Rev Genet. 2001, 2: 91-99.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2369/14/126/prepub
David H Kavanagh is supported by a Northern Ireland Kidney Research Fund (NIKRF) PhD studentship. We thank the NIKRF and Northern Ireland Research and Development Office for supporting this work. We thank Dr Denise Sadlier, University College Dublin, for providing DNA samples from cases and controls from the Republic of Ireland. The Warren 3 / UK GoKinD Study Group was jointly funded by Diabetes UK and the Juvenile Diabetes Research Foundation and includes the following individuals: Belfast: Professor A. P. Maxwell, Dr A. J. McKnight, Dr D. A. Savage; Edinburgh: Dr J. Walker; London: Dr S. Thomas, Professor G. C. Viberti; Manchester: Professor A. J. M. Boulton; Newcastle: Professor S. Marshall; Plymouth: Professor A. G. Demaine and Dr B. A. Millward; Swansea: Professor S. C. Bain. We also thank the Genetics of Nephropathy – an International Effort (GENIE) study group for provision of data. GENIE was jointly funded by NIH NIDDK R01 DK081923, The Northern Ireland Research Development office and Science Foundation Ireland US-Ireland R&D partnership, and The Health Research Board Ireland.
The authors declare that they have no competing interests.
Conceived and designed the experiments: DHK DAS JKC APM GJM. Performed the experiments: DHK. Analyzed the data: DHK CCP GJM. Contributed reagents/materials/analysis tools: AJM DAS APM JKC. Wrote the paper: DHK DAS APM CCP GJM. All authors have approved the final version of the manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Kavanagh, D.H., Savage, D.A., Patterson, C.C. et al. Haplotype association analysis of genes within the WNT signalling pathways in diabetic nephropathy. BMC Nephrol 14, 126 (2013). https://doi.org/10.1186/1471-2369-14-126