PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of wtpaEurope PMCEurope PMC Funders GroupSubmit a Manuscript
 
Diabetologia. Author manuscript; available in PMC 2009 May 28.
Published in final edited form as:
PMCID: PMC2687720
EMSID: UKMS4502

Genetic association analyses of non-synonymous single nucleotide polymorphisms in diabetic nephropathy

D. A. Savagecorresponding author
Nephrology Research Laboratory, Queen's University Belfast, c/o Regional Genetics Centre, Floor A, Tower Block, Belfast City Hospital, Lisburn Road, Belfast BT9 7AB, Northern Ireland, UK
C. C. Patterson
Department of Epidemiology and Public Health, School of Medicine and Dentistry, Queen's University Belfast, Belfast, UK
P. Deloukas and P. Whittaker
Wellcome Trust Sanger Institute, Cambridge, UK
A. J. McKnight
Nephrology Research Laboratory, Queen's University Belfast, c/o Regional Genetics Centre, Floor A, Tower Block, Belfast City Hospital, Lisburn Road, Belfast BT9 7AB, Northern Ireland, UK
J. Morrison
Wellcome Trust Sanger Institute, Cambridge, UK
A. J. Boulton
Manchester Royal Infirmary, University Department of Medicine, Manchester, UK
A. G. Demaine
Molecular Medicine Research Group, Peninsula Medical School, Universities of Exeter and Plymouth, Plymouth, UK
S. M. Marshall
Department of Diabetes and Metabolism, University of Newcastle upon Tyne, Newcastle upon Tyne, UK
B. A. Millward
Molecular Medicine Research Group, Peninsula Medical School, Universities of Exeter and Plymouth, Plymouth, UK
S. M. Thomas
Department of Diabetes and Endocrinology, Guy's and St Thomas' Hospitals, London, UK
G. C. Viberti
Unit for Metabolic Medicine, KCL School of Medicine, University of London, London, UK
J. D. Walker
Medical Unit, St John's Hospital, Livingston, UK
D. Sadlier
Conway Institute, University College Dublin, Dublin, Ireland
A. P. Maxwell
Nephrology Research Laboratory, Queen's University Belfast, c/o Regional Genetics Centre, Floor A, Tower Block, Belfast City Hospital, Lisburn Road, Belfast BT9 7AB, Northern Ireland, UK

Abstract

Aims/hypothesis

Diabetic nephropathy, characterised by persistent proteinuria, hypertension and progressive kidney failure, affects a subset of susceptible individuals with diabetes. It is also a leading cause of end-stage renal disease (ESRD). Non-synonymous (ns) single nucleotide polymorphisms (SNPs) have been reported to contribute to genetic susceptibility in both monogenic disorders and common complex diseases. The objective of this study was to investigate whether nsSNPs are involved in susceptibility to diabetic nephropathy using a case-control design.

Methods

White type 1 diabetic patients with (cases) and without (controls) nephropathy from eight centres in the UK and Ireland were genotyped for a selected subset of nsSNPs using Illumina's GoldenGate BeadArray assay. A χ2 test for trend, stratified by centre, was used to assess differences in genotype distribution between cases and controls. Genomic control was used to adjust for possible inflation of test statistics, and the False Discovery Rate method was used to account for multiple testing.

Results

We assessed 1,111 nsSNPs for association with diabetic nephropathy in 1,711 individuals with type 1 diabetes (894 cases, 817 controls). A number of SNPs demonstrated a significant difference in genotype distribution between groups before but not after correction for multiple testing. Furthermore, neither subgroup analysis (diabetic nephropathy with ESRD or diabetic nephropathy without ESRD) nor stratification by duration of diabetes revealed any significant differences between groups.

Conclusions/interpretation

The nsSNPs investigated in this study do not appear to contribute significantly to the development of diabetic nephropathy in patients with type 1 diabetes.

Keywords: Diabetic nephropathy, Genetic association study, Genetics of Kidneys in Diabetes, Genetic susceptibility, Non-synonymous single nucleotide polymorphisms

Introduction

Diabetic nephropathy is caused by hyperglycaemia-induced renal injury in genetically predisposed individuals. Patients affected with diabetic nephropathy have persistent proteinuria, hypertension and a declining glomerular filtration rate. Diabetic nephropathy is now the leading cause of end-stage renal disease (ESRD) and affected patients have a significantly increased cardiovascular mortality rate.

Despite convincing support for the existence of susceptibility and protective genes in diabetic nephropathy, genetic association studies have not yet identified consistently replicated associations in well-powered studies with large sample sizes [1]. Non-synonymous (ns) single nucleotide polymorphisms (SNPs) resulting in amino acid changes may directly affect protein structure and/or function. Evidence of association between nsSNPs and common multifactorial diseases has recently been reported [2-4]. It is also noteworthy that missense and nonsense mutations are the predominant type of mutation (~60%) reported in monogenic disorders [5]. We therefore reasoned that nsSNPs could directly contribute to genetic susceptibility to diabetic nephropathy. In the present study we tested our hypothesis by using Illumina's GoldenGate BeadArray assay (http://www.illumina.com/pages.ilmn?ID=11; [6]) to assess the largest panel of nsSNPs to date for association with diabetic nephropathy in over 1,700 individuals with type 1 diabetes.

Methods

Participants

The type 1 diabetic participants with (cases) and without (controls) diabetic nephropathy who were recruited in this study were white and had parents and grandparents born in the UK or Ireland. All participants were diagnosed with type 1 diabetes before the age of 35 years and required insulin from time of diagnosis. Cases and controls were derived from the Warren 3 (http://www.diabetes.org.uk/en/About_us/News_Landing_Page/2865/) and Genetics of Kidneys in Diabetes (GoKinD; http://www.gokind.org) collections recruited from seven centres in the UK (Belfast, Birmingham, Edinburgh, London, Manchester, Newcastle and Plymouth), as well as additional cases and controls from both the Belfast centre and a single centre (Dublin) in the Republic of Ireland.

Cases with nephropathy were defined by development of persistent proteinuria (>0.5 g protein/24 h) at least 10 years after diagnosis of diabetes, hypertension (BP >135/85 mmHg and/or treatment with antihypertensive agents) and presence of diabetic retinopathy. In contrast, controls were defined as patients with type 1 diabetes of at least 15 years' duration, but with urinary albumin in the normal range and not receiving antihypertensive treatment. Patients with microalbuminuria were excluded from both groups. Ethical approval was obtained from the appropriate research ethics committees, and written informed consent was obtained from all participants prior to conducting this study.

Non-synonymous SNP selection and genotyping

The nsSNPs surveyed were derived from candidate genes (1) identified following extensive literature reviews, (2) from reports on transcription profiling, (3) inferred by function or (4) known to be expressed in the kidney.

Bioinformatic interrogation was applied to the gene list to identify nsSNPs with a minor allele frequency (MAF) ≥ 0.05 in the HapMap CEPH (Utah residents with ancestry in northern and western Europe) collection (http://www.hapmap.org) to be considered for inclusion in a custom Illumina GoldenGate array comprising nsSNPs with an assay design score >0.4. Genotyping of SNPs was performed on DNA samples (250 ng) in 96-well plates using the Illumina GoldenGate assay as described elsewhere [6]. Each plate contained a combination of case and control DNA samples, one non-template control and duplicates of three different samples in positions unique to a given plate. Genotyping data were then imported into Illumina's Gencall data analysis software, which employs algorithms for clustering, calling and scoring of genotypes [6]. Assignment of GenCall scores permits ranking of the data and where appropriate assists in the identification and subsequent exclusion of failed DNAs, genotypes and/or loci; only nsSNPs with GenCall scores >0.25 were considered acceptable in this study.

Statistical analyses

The clinical characteristics of the patient groups were compared by the independent samples t test. The Cochran–Armitage χ2 test for trend, stratified by centre, was used to identify differences in genotype distribution between cases and controls [7]. Possible inflation of the χ2 test statistics was assessed by a quantile–quantile plot of observed χ2 test statistics against expected order statistics [8]. The lowest 90% of test statistics were averaged and the result was divided by the average of the lowest 90% of expected order statistics from a χ2 distribution to derive an inflation factor, λ. The False Discovery Rate method was used to account for multiple testing [9].

Results

The age at diagnosis of diabetes was not significantly different between the cases and controls, with a mean duration of diabetes of at least 27 years for each group (Table 1). As expected the average blood pressure was higher in cases than in controls, consistent with the renal disease phenotype and despite the concurrent use of antihypertensive drug treatment by patients with diabetic nephropathy (Table 1). Within the diabetic nephropathy group, a subset of 30% (n=265) had ESRD that was being managed by dialysis or renal transplantation. As reported for most national registries of renal disease, men were more likely to have diabetic nephropathy than women [10].

Table 1
Clinical characteristics of cases and controls

The nsSNPs with individual GenCall scores >0.25 were filtered on the basis of our criteria for call rates of at least 95% for cases and controls (of those meeting these criteria, 88% had call rates of at least 99%) and a differential in data completeness no greater than 2.0% between cases and controls. Further criteria were then applied that excluded nsSNPs with more than one duplicate discrepancy (n=16), a MAF<1% (n=16; nine of which were monomorphic) and departures from Hardy–Weinberg equilibrium (p<0.001) in either cases or controls (n=17). This resulted in high-quality genotyping data for 1,111 nsSNPs in 894 cases and 817 controls (Electronic Supplementary Material [ESM] Table 1).

These 1,111 SNPs, covering 1,021 unique genes (some SNPs overlap more than one gene, see ESM Table 2), were first compared in cases and controls using the unstratified χ2 test for trend, and a factor λ=1.121 (12.1% inflation, 95% CI 4.9–19.3) was obtained. Stratifying the χ2 statistics by centre reduced this inflation factor to λ=1.079 (7.9% inflation, 95% CI 1.1–14.7) (ESM Fig. 1). This indicated that the observed test statistics were 7.9% greater than would have been anticipated from the expected distribution of test statistics; the test statistics were adjusted accordingly. Before making allowance for the large number of tests (n=1,111) performed, 49 nsSNPs generated results that were significant at the p<0.05 level (55.5 expected by chance alone) and 13 nsSNPs gave results that were significant at the p< 0.01 level (11.1 expected by chance alone) (Table 2). None of the p values for these nsSNPs fell below the False Discovery Rate cut-off for significance when taking into account the number of tests (Table 2). We then analysed the data in subgroups defined by the presence of ESRD, i.e. diabetic nephropathy cases with ESRD (n=265) vs controls (n=817) and diabetic nephropathy cases without ESRD (n=629) vs controls (n=817). These additional subgroup analyses also revealed no significant differences in allele frequencies between groups.

Table 2
Distribution of allele frequencies (%) for 13 top-ranked nsSNPs (uncorrected p<0.01) in 894 cases and 817 controls

Given the relevance of diabetes duration in identifying susceptibility factors of moderate effect [11], the data were also analysed following stratification according to duration of diabetes (<25, 25–34 and ≥35 years). This stratification also failed to identify an association between any of the nsSNPs and diabetic nephropathy irrespective of diabetes duration.

Discussion

Statistical analysis of the 1,111 nsSNPs in 894 cases and 817 controls revealed no significant associations with diabetic nephropathy after correction for multiple genotype comparisons. Further subgroup analyses examining the influence of diabetes duration or renal status (diabetic nephropathy with ESRD or diabetic nephropathy without ESRD) also showed no significant associations between genotypes and diabetic nephropathy. However, individually a number of SNPs had statistically significant associations before correction for multiple testing (Table 2) and may warrant further follow-up, even though a large number of false positive findings are likely.

Several explanations might account for the lack of statistically significant associations between the nsSNPs analysed in the present study and diabetic nephropathy. One possibility is that we failed to identify an association due to a false negative outcome(s) (i.e. type 2 errors). The sample sizes employed in our analysis (894 cases, 817 controls) gave at least 80% power to detect, at p=5.0×10−5 (=0.05/1,000), a disease effect at a typed locus with odds ratio equal to 1.5 for a MAF of no less than 20%; however, for the same MAF we only had 3% and 18% power to detect associations for odds ratios of 1.2 and 1.3 respectively.

Another possibility for the lack of association is technical errors in genotyping. However, this is extremely unlikely, as we employed the highly accurate (>99%) and reproducible (>99%) GoldenGate assay from Illumina [6], in conjunction with strict quality control criteria.

In conclusion, it is likely that the nsSNPs analysed in this case–control study do not contribute significantly to a predisposition to diabetic nephropathy. Nevertheless 1,111 nsSNPs represent only a small proportion of nsSNPs [3], and the hypothesis, that nsSNPs may be involved in susceptibility to diabetic nephropathy remains plausible. To comprehensively test this hypothesis, a much larger nsSNP panel would have to be fully evaluated for association with diabetic nephropathy.

Supplementary Material

S1

S2

Acknowledgements

We thank Diabetes UK for supporting this project, and the Juvenile Diabetes Research Foundation and Diabetes UK for access to DNA samples from the GoKinD UK and Warren 3 collections respectively. We also thank T. Dibling and C. Hind at the Wellcome Trust Sanger Institute for genotyping. The work of P. Whittaker, J. Morrison and P Deloukas was supported by the Wellcome Trust.

Abbreviations

ESRD
end-stage renal disease
GoKinD
Genetics of Kidneys in Diabetes
MAF
minor allele frequency
ns
non-synonymous
SNP
single nucleotide polymorphism

Footnotes

Electronic supplementary material The online version of this article (doi:10.1007/s00125-008-1142-5) contains supplementary material, which is available to authorised users.

Duality of interest The authors declare that there is no duality of interest associated with this manuscript.

References

1. Savage DA, Bain SC, McKnight AJ, Maxwell AP. Gene discovery in diabetic nephropathy. Curr Diab Rep. 2007;7:139–145. [PubMed]
2. Smyth DJ, Cooper JD, Bailey R, et al. A genome-wide association study of nonsynonymous SNPs identifies a type 1 diabetes locus in the interferon-induced helicase (IFIH1) region. Nat Genet. 2006;38:617–619. [PubMed]
3. Wellcome Trust Case Control Consortium. The Australo–Anglo–American Spondylitis Consortium Association scan of 14,500 nonsynonymous SNPs in four diseases identifies autoimmunity variants. Nat Genet. 2007;39:1329–1337. [PMC free article] [PubMed]
4. Hampe J, Franke A, Rosenstiel P, et al. A genome-wide association scan of nonsynonymous SNPs identifies a susceptibility variant for Crohn disease in ATG16L1. Nat Genet. 2007;39:207–211. [PubMed]
5. Botstein D, Risch N. Discovering genotypes underlying human phenotypes: past successes for Mendelian disease, future approaches for complex disease. Nat Genet. 2003;33:228–237. [PubMed]
6. Fan JB, Oliphant A, Shen R, et al. Highly parallel SNP genotyping. Cold Spring Harb Symp Quant Biol. 2003;68:69–78. [PubMed]
7. Schlesselman JJ. Basic methods of analysis. In: Schlesselman JJ, editor. Case-control studies. Oxford University Press; New York: 1982. pp. 171–226.
8. Clayton DG, Walker NM, Smyth DJ, et al. Population structure, differential bias and genomic control in a large-scale, case-control association study. Nat Genet. 2005;37:1243–1246. [PubMed]
9. Sabatti C, Service S, Freimer N. False discovery rate in linkage and association genome screens for complex disorders. Genetics. 2003;164:829–833. [PubMed]
10. Ansell D, Feehally J, Feest TG, Tomson C, Williams AJ, Warwick G. UK Renal Registry Report 2007. UK Renal Registry; Bristol, UK: 2007.
11. Rogus JJ, Warram JH, Krolewski AS. Genetic studies of late diabetic complications: the overlooked importance of diabetes duration before complication onset. Diabetes. 2002;51:1655–1662. [PubMed]