Skip to main content


Identified single-nucleotide polymorphisms and haplotypes at 16q22.1 increase diabetic nephropathy risk in Han Chinese population



Diabetic nephropathy (DN) has become one of the most common causes of end-stage renal disease (ESRD) in many countries, such as 44.5% in Taiwan. Previous studies have shown that there is a genetic component to ESRD. Studies attempting to determine which genetic variants are related to DN in Han Chinese are limited. A case-control study was conducted to identify DN susceptibility variants in Han Chinese patients with type 2 diabetes.


We included 574 unrelated type 2 diabetes patients (217 DN cases and 357 controls), who were genotyped using Illumina HumanHap550-Duo BeadChip. In single-SNP association tests, the SNPs rs11647932, rs11645214, and rs6499323 located at 16q22.1 under the additive-effect disease model were significantly associated with an approximately 2-fold increased risk of DN. In haplotype association tests, identified haplotypes located in the chromosome 16q22.1 region (containing ST3GAL2, COG4, SF3B3, and IL34 genes) raised DN risk. The strongest association was found with haplotype rs2288491-rs4985534-rs11645214 (C-C-G) (adjusted odds ratio [AOR] 1.93, 95% confidence interval [CI] 1.83-2.03, p = 6.25 × 10-7), followed by haplotype rs8052125-rs2288491-rs4985534-rs11645214 (G-C-C-G) (AOR 1.92, 95% CI 1.82-2.02, p = 6.56 × 10-7), and haplotype rs2303792-rs8052125-rs2288491-rs4985534-rs11645214 (A-G-C-C-G) (AOR 1.91, 95% CI 1.81-2.01, p = 1.15 × 10-6).


Our results demonstrate that the novel SNPs and haplotypes located at the 16q22.1 region may involve in the biological pathways of DN in Han Chinese patients with type 2 diabetes. This study can provide new insights into the etiology of DN.


Diabetic nephropathy (DN) has become one of the most common causes of end-stage renal disease (ESRD) in many countries. In Taiwan, diabetes accounted for 44.5% of all new cases of ESRD [1], and more than 99% diabetes patients were with type 2 diabetes [2]. Diabetic ESRD patients had worse survival than non-diabetic ESRD patients [3].

It has been reported that genetic predisposition is one of the main risk factors for the development of DN [4]. Numerous familial aggregation studies have suggested that genetic susceptibility plays an important role in the development and progression of DN [5],[6]. Epidemiologic studies have shown that 35% of patients with diabetes develop nephropathy, irrespective of glycemic control [7],[8].

Mooyaart et al. performed a meta-analysis to evaluate the pooled effect of each genetic variant reproducibly associated with DN [9]. They reported that 21 of 34 replicated genetic variants remained significantly associated with DN. These 34 variants were in or near the following genes: ACE, ELMO1, PPARG, etc. Recently, several genome-wide association studies (GWASs) have attempted to detect genetic variants associated with the risk of DN or diabetic ESRD in those of Japanese [10], Pima Indian [11], and African American [12] with type 2 diabetes, as well as European ancestry [13],[14] with type 1 diabetes.

Studies attempting to determine which genetic variants are related to DN in Han Chinese patients with type 2 diabetes are limited. In the present genetic association study, a case-control study was carried out to identify DN susceptibility variants in Han Chinese patients with type 2 diabetes, which can provide new insights into the etiology of DN.


Study subjects

Individuals with type 2 diabetes and aged over 20 years were recruited using the American Diabetes Association (ICD-9-CM, Diagnosis code 250) criteria for diagnosis of type 2 diabetes. Individuals with type 1 diabetes, gestational diabetes, and maturity-onset diabetes of the young were excluded. The dataset used is part of the whole dataset for the published paper titled “A Genome-wide Association Study Identifies Susceptibility Variants for Type 2 Diabetes in Han Chinese” [15]. A total of 995 type 2 diabetes subjects recruited from China Medical University Hospital were included in the current study. All patients with type 2 diabetes were of Han Chinese origin, including Minnan, Hakka, and Mainland Chinese. Individuals with significant aboriginal ancestry were excluded by using self-administered questionnaires with six items regarding the ancestral origin of their parents and grandparents. Chronic kidney disease (CKD) was determined by estimated glomerular filtration rate (eGFR) and urine protein. Diabetic patients with eGFR <60 mL/min/1.73 m2 or with proteinuria determined by a spot urine dipstick of >1+ were defined as DN cases. A total of 217 DN cases were eligible for the study. To increase comparability between DN cases and controls, 357 controls were randomly selected based on frequency-matching of their age and durations of diabetes. This study was approved by the Human Research Committee of China Medical University Hospital. All patients signed informed consent forms.


Blood samples were collected in the morning after a 12-h overnight fast and were sent for analysis within 4-h of collection. Spot morning urine samples were collected. Triglycerides, total cholesterol, high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), creatinine, blood urea nitrogen (BUN), and uric acid were measured by a biochemical autoanalyser (Beckman Coulter, Synchron LX20, Fullerton, CA, USA). Renal function was evaluated by eGFR, which was estimated by using the Modification of Diet in Renal Disease Study equation for Taiwanese: eGFR (ml/min/1.73 m2) = 175 × (serum creatinine (mg/dL)-1.154 × (age)-0.203 × (0.742 if female) × 0.945 [16]. We also used the Chronic Kidney Disease Epidemiology Collaboration equation to obtain eGFR [17]. The proportions of eGFR <60 mL/min/1.73 m2 by using above two equations were similar. A spot urine dipstick test was used for detecting proteinuria. Patients with a positive dipstick test (>1+) were classified as with proteinuria [18]. The sociodemographic and lifestyle characteristics and the self-reported health status for each subject were recorded using self-administered questionnaires.

Genotyping and quality control

Genomic DNA from peripheral blood was prepared using the Puregene DNA isolation kit (Gentra Systems, Minneapolis, MN, USA). Our samples were genotyped using Illumina HumanHap550-Duo BeadChip, which was performed by deCODE Genetics (Reykjavík, Iceland). Genotypes were called using the standard procedure in BeadStudio (Illumina, Inc., San Diego, CA, USA), with the default parameters recommended by the platform manufacturer. The genotyping quality control procedures used to identify and remove poor-quality data were described previously [15]. Individual SNPs were excluded if they had a total call rate <95% across all individuals, a minor allele frequency <5% and a total call rate <99%, or had significant deviation from Hardy-Weinberg disequilibrium (p-value <10-7) in these subjects. Further details of genotyping quality control procedures are available in the published study [15]. After SNPs quality control, 429,018 SNPs were used and their total call rate was 99.9%. We also examined population stratification by using multidimensional scaling (MDS) analysis as implemented in PLINK. The results of MDS analysis showed that there was no evidence for population stratification.

Statistical analysis

Demographic and clinical characteristics of study subjects were examined, including sex, age, diabetes-related variables, behaviors, biochemical variables, and history of diseases. Continuous variables are reported as mean × standard deviation (SD), and categorical variables are reported as number and percentage. Two-sample t tests and Chi-square tests were used for the bivariate analyses. Because the distribution of the triglycerides was skewed, the data were normalized using a natural log-transformation, and the geometric mean × SD was calculated. To identify the DN susceptibility variants, single-SNP association tests using Cochran-Armitage trend test were performed. Then, multiple logistic regression analysis using an additive-effect disease model (an ordinal genotype model: 0, 1 and 2 of a minor allele) was performed for each SNP with adjustment of the subject's sex, age, BMI, and durations of diabetes. The Bonferroni correction was used to adjust for multiple comparisons. P-value less than 10-4, association was considered to be statistically significant. According to previously published studies, as well as the results of single-SNP association tests and the Manhattan plot in our current study, potential susceptibility regions for DN were selected to perform haplotype analysis. In haplotype-based association analysis, the sliding window approach was adopted to detect haplotype effects. The window sizes of 3-SNP, 4-SNP, and 5-SNP haplotypes were used. Each haplotype with a frequency of >0.05 in this population was analyzed. Odds ratios (ORs) and their corresponding 95% confidence intervals (CIs) were calculated to estimate the effect sizes of the identified SNPs and haplotypes. In addition, the linkage disequilibrium (LD) structures of the identified contiguous SNPs were examined. Pairwise LD was measured by the r2 statistic. For power calculation of our case-control study (217 DN cases and 357 controls), Quanto software [19] was used. Under an additive effect disease model with a prevalence of 10% for DN (from our dataset), given a genetic relative risk of 1.85 and a disease allele frequency of 0.25-0.45, the power of our study was 0.76-0.86 at an alpha level of 10-4. All analyses were carried out using Haploview (v4.2) [20], PLINK (v1.07) ( [21], and SAS (v9.3, SAS Institute Inc, Cary, NC, USA) software. The regional plot was plotted from the LocusZoom, a web-based plotting tool ( [22]. The in silico prediction tool is-rSNP was used to predict potential regulatory SNPs (rSNPs) [23].


Demographic and clinical characteristics of study samples

A total of 574 type 2 diabetes study participants, comprising 217 cases of DN and 357 type 2 diabetes patients without DN controls, were included in the analysis. Table 1 shows the sociodemographic factors, diabetes-related variables, lifestyle behaviors, biochemical variables, and history of diseases according to DN status. The mean age of the participants was about 62 years for both groups and their mean HbA1c at enrollment was 8.2 × 1.7% for the cases and 7.9 × 1.4% for the controls. The mean durations of diabetes were 11.1 × 7.7 years for the cases and 10.2 × 7.2 years for the controls. Compared with the controls, DN cases had higher mean BMI, triglycerides, creatinine, uric acid, BUN, and prevalence of hypertension and heart disease.

Table 1 Demographic and clinical characteristics of study samples

The graphical summary plot and single-SNP association tests

The graphical summary plot of DN in single-SNP association analysis under additive-effect disease model is shown in Figure 1. There were four SNPs with p-value <10-5 (above the red line), including rs10963767 (intron; ADAMTSL1 gene; 9p21.3), rs2058289 (intergenic; 12q24.3), rs11645214 (3′ UTR; SF3B3 gene; 16q22.1), and rs6499323 (intron; IL34; 16q22.1). Furthermore, two signals with p-value <10-4 were observed, indicating the SNPs around this areas may be potential susceptibility variants for DN; one signal (rs876142, rs9928626, rs11647932, rs11645214, and rs6499323) was on chromosome 16 (16q21-16q22.1 regions) and the other (rs1028555, rs4815800, rs6065925, rs6131015, rs4812997, rs6074024, rs6127983, rs182784, rs1885580, rs7273764, rs6127999, rs6014975, rs4811839, rs6025517, and rs2426712) was on chromosome 20 (20p12.3 and 20q13.1-20q13.3 regions). Among these 15 SNPs on chromosome 20, only SNPs rs6127983 and rs182784 (intron; BMP7 gene; 20q13.3) and SNPs rs4811839 and rs6025517 (intron; RAE1 gene; 20q13.3) were inside of genes. A previous study reviewed in a meta-analysis by Mooyaart et al. showed that the Hp 1/2 variant in the HP gene (located at 16q22.2) was associated with DN. Therefore, a total of 383 genotyped SNPs at the 16q22.1-16q22.2 regions were selected for this study. Significant SNPs (p-value <10-4) identified from single-SNP association analysis, along with the other 8 SNPs identified from haplotype association analysis, are presented in Table 2. The SNP with the lowest p-value (6.25 × 10-7) was rs11645214 located on chromosome 16q22.1. This SNP, with a higher risk of DN (adjusted odds ratio [AOR] 1.94, 95% confidence interval [CI] 1.49-2.51), is inside the SF3B3 gene (Additional file 1: Figure S1). Two other SNPs that showed a significant or borderline significant association with DN, each with an additive-effect, were rs6499323 (IL34) and rs11647932 (ST3GAL2) on chromosome 16q22.1 (Additional file 1: Figure S1). No significant variants were detected for DN at the 16q22.2 region.

Figure 1

Graphical summary of DN in single-SNP association analysis under additive-effect disease model.

Table 2 Single-SNP association analysis for DN at the 16q22.1-16q22.2 regions under the additive model a

The LD structures and haplotypic association tests

Pairwise LDs among the 11 SNPs located at 16q22.1 region are shown in Figure 2. SNPs rs8062183 in ST3GAL2, rs2303792 and rs8052125 in COG4, and rs2288491 in SF3B3 were in strong LD with each other (r2 > 0.94). SNP rs11645214 in SF3B3 was highly in LD with rs6499323 in IL34 (r2 = 0.92). Nevertheless, SNP rs11647932 and other SNPs were not in LD with r2 values ranging from 0.04 to 0.30.

Figure 2

Pairwise linkage disequilibrium (LD) among the 11 SNPs located at 16q22.1 region. Pairwise LD was measured by the r2 statistic. An r2 value of 1 indicates a complete LD between two SNPs, and an r2 value of 0 indicates a perfect equilibrium. Lower r2 value shows a lower degree of LD. Pairwise D′ values along with LOD scores were displayed in different colors. White indicates D′ <1 and LOD <2; blue indicate D′ = 1 and LOD <2; bright red indicate D′ = 1 and LOD <2; shades of pink/red indicate D′ <1 and LOD ≥2.

Genotyped SNPs at the 16q22.1-16q22.2 regions were selected for the haplotypic association tests in this study. The window sizes of 3-SNP, 4-SNP, and 5-SNP haplotypes were used to detect haplotype effects on DN, and a frequency of >0.05 was analyzed. The haplotypes with a p-value <10-4 are listed in Table 3. There were 4, 3, and 4 haplotypes with p-value <10-5 in window size of 3-SNP, 4-SNP, and 5-SNP significantly associated with DN in this sample, respectively, after adjusted for sex, age, BMI, and durations of diabetes. In particular, haplotype rs2288491-rs4985534-rs11645214 (C-C-G) (AOR 1.93, 95% CI 1.83-2.03, p = 6.25 × 10-7), haplotype rs8052125-rs2288491-rs4985534-rs11645214 (G-C-C-G) (AOR 1.92, 95% CI 1.82-2.02, p = 6.56 × 10-7), and haplotype rs2303792-rs8052125-rs2288491-rs4985534-rs11645214 (A-G-C-C-G) (AOR 1.91, 95% CI 1.81-2.01, p = 1.15 × 10-6) seem to be associated with a higher risk of DN. All haplotypes with significant associations with DN consisted of SNPs that exerted significant associations with DN in single-SNP association analysis.

Table 3 Haplotypic association tests for DN at the 16q22.1-16q22.2 regions a


There is very limited evidence that any of genetic variants contribute to DN in Han Chinese populations. In this study, we expand on the previous investigation by Tsai et al. [15], which described compelling evidence of an association of type 2 diabetes with polymorphisms. We found that p-values of four SNPs were less than 10-5 in the Figure 1, and two signals, i.e. p-value <10-4, were observed. We identified 11 SNPs in four loci located in 16q22.1 that were associated with DN in a Han Chinese population. As a consequence, we were able to identify novel loci specifically associated with DN. Two SNPs of two loci with p-values in the 10-5 to 10-7 range were identified for DN. These include SNPs rs11645214 (SF3B3) and rs6499323 (IL34), with about a 2-fold increased risk of DN. In the haplotype association analysis, 15 haplotypes, including these three SNPs (rs11647932, rs11645214, and rs6499323), carried an approximately 2-fold increased risk of DN, which further confirmed the findings of the single-SNP association tests.

Previous studies using GWAS approach on DN in different racial groups have been conducted [10]-[14]. These studies showed that different polymorphisms in or near genes may be related to DN susceptibility, for example, intron 18 + 9170 (ELMO1) in Japanese with type 2 diabetes [10]; rs2648875 (PVT1) and rs2720709 (PVT1) in Pima Indians with type 2 diabetes [11]; rs6930576 (SASH1), rs7769051 (near RPS12), rs2358944 (MSRB3 ~ HMGA2), rs2106294 (LIMK2), rs4820043 (LIMK2), rs7735506 (AUH), and rs5749286 (SFI1) in African Americans with type 2 diabetes [12]; rs10868025 (near FRMD3), rs1888747 (near FRMD3), rs451041 (CARS), and rs739401 (CARS) in European ancestry with type 1 diabetes [13]; and rs12437854 (RGMA ~ MCTP2), rs7583877 (AFF3), and rs7588550 (ERBB4) in European ancestry with type 1 diabetes [14] (Additional file 1: Table S1). The susceptibility variants for DN varied widely among different population groups, which may be due to ethnic differences and genetic backgrounds. In this study, we found that polymorphisms in the ST3GAL2, COG4, SF3B3, IL34, and BMP7 genes may contribute to DN susceptibility in Han Chinese residing in Taiwan. In previous candidate gene studies, these genes have not been reported for DN in Han Chinese (Additional file 1: Table S2). The proportions of minor alleles of SNPs rs11647932 and rs6499323 among Taiwanese (14% and 34%, respectively) were similar to those of Han Chinese in China (CHB: 14% and 36%), but higher than those of Japanese (JPT: 12% and 28%), Europeans (CEU: 12% and 20%), and Africans (YRI: 0% and 12%) [24].

Many plausible mechanisms have been hypothesized to be involved in the association between DN and identified SNPs. First, the protein encoded by the ST3GAL2 gene mapped on chromosome 16q22.1 catalyzes the transfer of sialic acid from CMP-sialic acid to galactose-containing substrates [25]. Sialic acid is one of the inflammatory biomarkers of the acute-phase response and is a possible risk factor for cardiovascular disease [26]. Sialic acid was synthesized as a result of cytokines released by inflammatory cells and damaged endothelia due to angiopathy-associated tissue injury [26]. This cytokine response could directly result in vasculopathy in diabetes through atherosclerosis [26]. Moreover, cross-sectional and longitudinal studies have shown that elevated serum sialic acid concentrations are related to diabetic microvascular complications (including DN) in type 2 diabetes [27],[28] or type 1 diabetes [29],[30]. A strong association between increased serum sialic acid concentrations and macrovascular complications (coronary heart disease; CHD) in type 2 diabetes was also observed [31]. In addition, a 7-year prospective cohort study revealed that a raised serum sialic acid concentration is associated with CHD in type 1 diabetes [32].

Second, the protein encoded by the COG4 (a component of oligomeric golgi complex 4) gene, including COG4, is critical for the structure and function of the Golgi apparatus and can influence intracellular membrane trafficking, including the addition and processing of carbohydrates (glycosylation) in the rough endoplasmic reticulum and Golgi [33],[34]. Mutations in COG4 may lead to abnormal functioning in the Golgi apparatus, resulting in an excessive amount of glucose circulating in the blood plasma. Diabetes mellitus is characterized by chronic hyperglycaemia. A number of studies have examined the relationship of hyperglycaemia, classified according to mean levels of HbA1c [35]-[37] or fasting plasma glucose [38],[39], with diabetic complications.

Third, the IL34 (interleukin 34) gene maps on chromosome 16q22.1. Interleukin 34 is a cytokine that promotes the differentiation and viability of monocytes and macrophages through the colony-stimulating factor-1 receptor (CSF1R) [40], and the causative role of macrophages and/or monocytes in DN has been demonstrated by cell depletion studies in animal models [41],[42]. Many studies have shown that inflammatory cytokines such as IL-1, IL-6, IL-18, and TNF-α are important in the pathogenesis of diabetic microvascular complications (including DN) [43]. Furthermore, previous studies have revealed that cytokine gene polymorphisms such as IL-1, IL-6, and TNF-α are related to DN [44]. Variants in IL34 may play a role in the pathogenesis of DN because of the highly correlated cytokine network, although there is no direct evidence for the association between IL34 gene polymorphisms and DN.

Fourth, the protein encoded by the BMP7 (bone morphogenetic protein 7) gene mapped on 20q13.3 is a member of the transforming growth factor-beta (TGF-α) superfamily, and is highly expressed in the kidney tubules and glomeruli [45]. Previous studies showed that TGF-α cellular signaling is critical to the induction of glomerular and tubulointerstitial fibrosis in DN [46]. Animal studies demonstrated that the renal BMP7 protected against DN [47], and the BMP7 partially reversed kidney hypertrophy induced by diabetes, urine albumin excretion, restoring GFR, as well as glomerular histology toward normal [48]. McKnight et al. investigated the association of BMP gene variants with DN in White individuals with type 1 diabetes, and they reported that the common variants in BMP2, BMP4 and BMP7 genes did not strongly influence genetic susceptibility to DN [49]. In our study, we found that two SNPs in BMP7 gene may be susceptibility variants of DN in Han Chinese patients with type 2 diabetes.

For predicting potential rSNPs, we used the in silico prediction tool is-rSNP [23]. The is-rSNP predicted that rs2058289, rs9928626, rs11647932, rs1028555, rs6127999, rs6131015, rs7273764, and rs4815800 are potential rSNPs; that is, they affect the binding affinity of a transcription factor to the DNA (Additional file 1: Table S3). For example, alleles of rs2058289 alter the binding affinity of Sp4_2 (p = 0.021) and Eomes_1 (p = 0.036); alleles of rs9928626 alter the binding affinity of LYS14 (p = 0.018), YKL222C (p = 0.043), and YDR520C (p = 0.043); and alleles of rs11647932 alter the binding affinity of LM226 (p = 0.050). Furthermore, from the online Nephromine database (, there were some evidences in the literature that the expression of candidate genes in our study were associated with DN in human population (Additional file 1: Table S4) or mice (Additional file 1: Table S5).

The subjects in this study were a portion of the cohort studied in a published GWAS paper [15]. Although the sample size in this study was limited, the DN cases and the controls were much more homogeneous in terms of gender, age, durations of diabetes, and HbA1c level after excluding subjects with extreme characteristics and matching their age and durations of diabetes. A previous study found that the Han Chinese residing in Taiwan were relatively homogenous in genetic background, spread among the three major ethnic groups of Minnan, Hakka, and Mainland Chinese [50], which is consistent with the results of our population stratification analysis. Therefore, the impact of population stratification on our genetic association study is small.


We identified 3 novel SNPs (rs11647932, rs11645214, and rs6499323) and 11 haplotypes (4, 3, and 4 in window size of 3-SNP, 4-SNP, and 5-SNP, respectively), located at the 16q22.1 region, that are susceptibility variants of DN in a Han Chinese population in Taiwan. The findings reveal that the susceptibility variants located at the 16q22.1 region may involve in the biological pathways of DN in Han Chinese patients with type 2 diabetes. Our results can provide new insights into the etiology of DN.

Additional file



Diabetic nephropathy


End-stage renal disease


Chronic kidney disease


Estimated glomerular filtration rate


Hhigh-density lipoprotein cholesteroll


Low-density lipoprotein cholesterol


Blood urea nitrogen


Multidimensional scaling


Linkage disequilibrium


  1. 1.

    Collins AJ, Foley RN, Chavers B, Gilbertson D, Herzog C, Johansen K, Kasiske B, Kutner N, Liu J, St Peter W, Guo H, Gustafson S, Heubner B, Lamb K, Li S, Li S, Peng Y, Qiu Y, Roberts T, Skeans M, Snyder J, Solid C, Thompson B, Wang C, Weinhandl E, Zaun D, Arko C, Chen SC, Daniels F, Ebben J: United States Renal Data System 2011 Annual Data Report: Atlas of chronic kidney disease & end-stage renal disease in the United States. Am J Kidney Dis. 2012, 59 (A7): e1-e420.

  2. 2.

    Jiang YD, Chang CH, Tai TY, Chen JF, Chuang LM: Incidence and prevalence rates of diabetes mellitus in Taiwan: analysis of the 2000-2009 Nationwide Health Insurance database. J Formos Med Assoc. 2012, 111: 599-604. 10.1016/j.jfma.2012.09.014.

  3. 3.

    Yang WC, Hwang SJ: Incidence, prevalence and mortality trends of dialysis end-stage renal disease in Taiwan from 1990 to 2001: the impact of national health insurance. Nephrol Dial Transplant. 2008, 23: 3977-3982. 10.1093/ndt/gfn406.

  4. 4.

    Gross JL, de Azevedo MJ, Silveiro SP, Canani LH, Caramori ML, Zelmanovitz T: Diabetic nephropathy: diagnosis, prevention, and treatment. Diabetes Care. 2005, 28: 164-176. 10.2337/diacare.28.1.164.

  5. 5.

    Seaquist ER, Goetz FC, Rich S, Barbosa J: Familial clustering of diabetic kidney disease. Evidence for genetic susceptibility to diabetic nephropathy. N Engl J Med. 1989, 320: 1161-1165. 10.1056/NEJM198905043201801.

  6. 6.

    Imperatore G, Knowler WC, Pettitt DJ, Kobes S, Bennett PH, Hanson RL: Segregation analysis of diabetic nephropathy in Pima Indians. Diabetes. 2000, 49: 1049-1056. 10.2337/diabetes.49.6.1049.

  7. 7.

    Andersen AR, Christiansen JS, Andersen JK, Kreiner S, Deckert T: Diabetic nephropathy in Type 1 (insulin-dependent) diabetes: an epidemiological study. Diabetologia. 1983, 25: 496-501. 10.1007/BF00284458.

  8. 8.

    Krolewski AS, Warram JH, Christlieb AR, Busick EJ, Kahn CR: The changing natural history of nephropathy in type I diabetes. Am J Med. 1985, 78: 785-794. 10.1016/0002-9343(85)90284-0.

  9. 9.

    Mooyaart AL, Valk EJ, van Es LA, Bruijn JA, de Heer E, Freedman BI, Dekkers OM, Baelde HJ: Genetic associations in diabetic nephropathy: a meta-analysis. Diabetologia. 2011, 54: 544-553. 10.1007/s00125-010-1996-1.

  10. 10.

    Shimazaki A, Kawamura Y, Kanazawa A, Sekine A, Saito S, Tsunoda T, Koya D, Babazono T, Tanaka Y, Matsuda M, Kawai K, Iiizumi T, Imanishi M, Shinosaki T, Yanagimoto T, Ikeda M, Omachi S, Kashiwagi A, Kaku K, Iwamoto Y, Kawamori R, Kikkawa R, Nakajima M, Nakamura Y, Maeda S: Genetic variations in the gene encoding ELMO1 are associated with susceptibility to diabetic nephropathy. Diabetes. 2005, 54: 1171-1178. 10.2337/diabetes.54.4.1171.

  11. 11.

    Hanson RL, Craig DW, Millis MP, Yeatts KA, Kobes S, Pearson JV, Lee AM, Knowler WC, Nelson RG, Wolford JK: Identification of PVT1 as a candidate gene for end-stage renal disease in type 2 diabetes using a pooling-based genome-wide single nucleotide polymorphism association study. Diabetes. 2007, 56: 975-983. 10.2337/db06-1072.

  12. 12.

    McDonough CW, Palmer ND, Hicks PJ, Roh BH, An SS, Cooke JN, Hester JM, Wing MR, Bostrom MA, Rudock ME, Lewis JP, Talbert ME, Blevins RA, Lu L, Ng MC, Sale MM, Divers J, Langefeld CD, Freedman BI, Bowden DW: A genome-wide association study for diabetic nephropathy genes in African Americans. Kidney Int. 2011, 79: 563-572. 10.1038/ki.2010.467.

  13. 13.

    Pezzolesi MG, Poznik GD, Mychaleckyj JC, Paterson AD, Barati MT, Klein JB, Ng DP, Placha G, Canani LH, Bochenski J, Waggott D, Merchant ML, Krolewski B, Mirea L, Wanic K, Katavetin P, Kure M, Wolkow P, Dunn JS, Smiles A, Walker WH, Boright AP, Bull SB, Doria A, Rogus JJ, Rich SS, Warram JH, Krolewski AS: Genome-wide association scan for diabetic nephropathy susceptibility genes in type 1 diabetes. Diabetes. 2009, 58: 1403-1410. 10.2337/db08-1514.

  14. 14.

    Sandholm N, Salem RM, McKnight AJ, Brennan EP, Forsblom C, Isakova T, McKay GJ, Williams WW, Sadlier DM, Makinen VP, Swan EJ, Palmer C, Boright AP, Ahlqvist E, Deshmukh HA, Keller BJ, Huang H, Ahola AJ, Fagerholm E, Gordin D, Harjutsalo V, He B, Heikkilä O, Hietala K, Kytö J, Lahermo P, Lehto M, Lithovius R, Osterholm AM, Parkkonen M, et al: New susceptibility loci associated with kidney disease in type 1 diabetes. PLoS Genet. 2012, 8: e1002921-10.1371/journal.pgen.1002921.

  15. 15.

    Tsai FJ, Yang CF, Chen CC, Chuang LM, Lu CH, Chang CT, Wang TY, Chen RH, Shiu CF, Liu YM, Chang CC, Chen P, Chen CH, Fann CS, Chen YT, Wu JY: A genome-wide association study identifies susceptibility variants for type 2 diabetes in Han Chinese. PLoS Genet. 2010, 6: e1000847-10.1371/journal.pgen.1000847.

  16. 16.

    The Report of Chronic Kidney Disease Prevention, 2011 [in Chinese]. [], []

  17. 17.

    Levey AS, Stevens LA, Schmid CH, Zhang YL, Castro AF, Feldman HI, Kusek JW, Eggers P, Van Lente F, Greene T, Van Lente F, Greene T, Coresh J: A new equation to estimate glomerular filtration rate. Ann Intern Med. 2009, 150: 604-612. 10.7326/0003-4819-150-9-200905050-00006.

  18. 18.

    K/DOQI clinical practice guidelines for chronic kidney disease: evaluation, classification, and stratification. Am J Kidney Dis. 2002, 39: S1-S266.

  19. 19.

    Gauderman WJ, Morrison JM: QUANTO 1.1: A computer program for power and sample size calculations for genetic-epidemiology studies, 2006

  20. 20.

    Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21: 263-265. 10.1093/bioinformatics/bth457.

  21. 21.

    Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, Sham PC: PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007, 81: 559-575. 10.1086/519795.

  22. 22.

    Pruim RJ, Welch RP, Sanna S, Teslovich TM, Chines PS, Gliedt TP, Boehnke M, Abecasis GR, Willer CJ: LocusZoom: regional visualization of genome-wide association scan results. Bioinformatics. 2010, 26: 2336-2337. 10.1093/bioinformatics/btq419.

  23. 23.

    Macintyre G, Bailey J, Haviv I, Kowalczyk A: is-rSNP: a novel technique for in silico regulatory SNP detection. Bioinformatics. 2010, 26: i524-i530. 10.1093/bioinformatics/btq378.

  24. 24.

    Data Rel 27 Phase II+III, Feb 09. International HapMap Consortium.. Accessed May 21, 2013., []

  25. 25.

    Taniguchi A, Morishima T, Tsujita Y, Matsumoto Y, Matsumoto K: Genomic structure, expression, and transcriptional regulation of human Gal beta 1,3 GalNAc alpha 2,3-sialyltransferase gene. Biochem Biophys Res Commun. 2003, 300: 570-576. 10.1016/S0006-291X(02)02899-1.

  26. 26.

    Gopaul KP, Crook MA: Sialic acid: a novel marker of cardiovascular disease?. Clin Biochem. 2006, 39: 667-681. 10.1016/j.clinbiochem.2006.02.010.

  27. 27.

    Chen J, Gall MA, Yokoyama H, Jensen JS, Deckert M, Parving HH: Raised serum sialic acid concentration in NIDDM patients with and without diabetic nephropathy. Diabetes Care. 1996, 19: 130-134. 10.2337/diacare.19.2.130.

  28. 28.

    Nayak BS, Roberts L: Relationship between inflammatory markers, metabolic and anthropometric variables in the Caribbean type 2 diabetic patients with and without microvascular complications. J Inflamm. 2006, 3: 17-10.1186/1476-9255-3-17.

  29. 29.

    Crook MA, Pickup JC, Lumb PJ, Giorgino F, Webb DJ, Fuller JH, Group EICS: Relationship between plasma sialic acid concentration and microvascular and macrovascular complications in type 1 diabetes: the EURODIAB Complications Study. Diabetes Care. 2001, 24: 316-322. 10.2337/diacare.24.2.316.

  30. 30.

    Yokoyama H, Jensen JS, Myrup B, Mathiesen ER, Ronn B, Deckert T: Raised serum sialic acid concentration precedes onset of microalbuminuria in IDDM. A 10-year follow-up study. Diabetes Care. 1996, 19: 435-440. 10.2337/diacare.19.5.435.

  31. 31.

    Pickup JC, Mattock MB, Crook MA, Chusney GD, Burt D, Fitzgerald AP: Serum sialic acid concentration and coronary heart disease in NIDDM. Diabetes Care. 1995, 18: 1100-1103. 10.2337/diacare.18.8.1100.

  32. 32.

    Soedamah-Muthu SS, Chaturvedi N, Pickup JC, Fuller JH, Group EPCS: Relationship between plasma sialic acid and fibrinogen concentration and incident micro- and macrovascular complications in type 1 diabetes. The EURODIAB Prospective Complications Study (PCS). Diabetologia. 2008, 51: 493-501. 10.1007/s00125-007-0905-8.

  33. 33.

    Lodish H, Berk A, Matsudaira P, Kaiser CA, Krieger M, Scott MP, Zipursky L, Darnell J: Molecular Cell Biology. 2004, W.H. Freeman and Company, New York

  34. 34.

    Ungar D, Oka T, Brittle EE, Vasile E, Lupashin VV, Chatterton JE, Heuser JE, Krieger M, Waters MG: Characterization of a mammalian Golgi-localized protein complex, COG, that is required for normal Golgi morphology and function. J Cell Biol. 2002, 157: 405-415. 10.1083/jcb.200202016.

  35. 35.

    Hemmingsen B, Lund SS, Gluud C, Vaag A, Almdal T, Hemmingsen C, Wetterslev J: Intensive glycaemic control for patients with type 2 diabetes: systematic review with meta-analysis and trial sequential analysis of randomised clinical trials. BMJ. 2011, 343: d6898-10.1136/bmj.d6898.

  36. 36.

    Tu ST, Chang SJ, Chen JF, Tien KJ, Hsiao JY, Chen HC, Hsieh MC: Prevention of diabetic nephropathy by tight target control in an asian population with type 2 diabetes mellitus: a 4-year prospective analysis. Arch Intern Med. 2010, 170: 155-161. 10.1001/archinternmed.2009.471.

  37. 37.

    Intensive blood-glucose control with sulphonylureas or insulin compared with conventional treatment and risk of complications in patients with type 2 diabetes (UKPDS 33). UK Prospective Diabetes Study (UKPDS) Group. Lancet. 1998, 352: 837-853. 10.1016/S0140-6736(98)07019-6.

  38. 38.

    Lawes CM, Parag V, Bennett DA, Suh I, Lam TH, Whitlock G, Barzi F, Woodward M: Blood glucose and risk of cardiovascular disease in the Asia Pacific region. Diabetes Care. 2004, 27: 2836-2842. 10.2337/diacare.27.12.2836.

  39. 39.

    Sorkin JD, Muller DC, Fleg JL, Andres R: The relation of fasting and 2-h postchallenge plasma glucose concentrations to mortality: data from the Baltimore Longitudinal Study of Aging with a critical review of the literature. Diabetes Care. 2005, 28: 2626-2632. 10.2337/diacare.28.11.2626.

  40. 40.

    Lin H, Lee E, Hestir K, Leo C, Huang M, Bosch E, Halenbeck R, Wu G, Zhou A, Behrens D, Hollenbaugh D, Linnemann T, Qin M, Wong J, Chu K, Doberstein SK, Williams LT: Discovery of a cytokine and its receptor by functional screening of the extracellular proteome. Science. 2008, 320: 807-811. 10.1126/science.1154370.

  41. 41.

    Shikata K, Makino H: Role of macrophages in the pathogenesis of diabetic nephropathy. Contrib Nephrol. 2001, 134: 46-54. 10.1159/000060147.

  42. 42.

    Tesch GH: Role of macrophages in complications of type 2 diabetes. Clin Exp Pharmacol Physiol. 2007, 34: 1016-1019. 10.1111/j.1440-1681.2007.04729.x.

  43. 43.

    Navarro-Gonzalez JF, Mora-Fernandez C, Muros De Fuentes M, Garcia-Perez J: Inflammatory molecules and pathways in the pathogenesis of diabetic nephropathy. Nat Rev Nephrol. 2011, 7: 327-340. 10.1038/nrneph.2011.51.

  44. 44.

    Navarro-Gonzalez JF, Mora-Fernandez C: The role of inflammatory cytokines in diabetic nephropathy. J Am Soc Nephrol. 2008, 19: 433-442. 10.1681/ASN.2007091048.

  45. 45.

    Zhang Y, Zhang Q: Bone morphogenetic protein-7 and gremlin: new emerging therapeutic targets for diabetic nephropathy. Biochem Biophys Res Commun. 2009, 383: 1-3. 10.1016/j.bbrc.2009.03.086.

  46. 46.

    Kanwar YS, Sun L, Xie P, Liu FY, Chen S: A glimpse of various pathogenetic mechanisms of diabetic nephropathy. Annu Rev Pathol. 2011, 6: 395-423. 10.1146/annurev.pathol.4.110807.092150.

  47. 47.

    Wang S, de Caestecker M, Kopp J, Mitu G, Lapage J, Hirschberg R: Renal bone morphogenetic protein-7 protects against diabetic nephropathy. J Am Soc Nephrol. 2006, 17: 2504-2512. 10.1681/ASN.2006030278.

  48. 48.

    Wang S, Chen Q, Simon TC, Strebeck F, Chaudhary L, Morrissey J, Liapis H, Klahr S, Hruska KA: Bone morphogenic protein-7 (BMP-7), a novel therapy for diabetic nephropathy. Kidney Int. 2003, 63: 2037-2049. 10.1046/j.1523-1755.2003.00035.x.

  49. 49.

    McKnight AJ, Pettigrew KA, Patterson CC, Kilner J, Sadlier DM, Maxwell AP, Warren UKGSG: Investigation of the association of BMP gene variants with nephropathy in Type 1 diabetes mellitus. Diabet Med. 2010, 27: 624-630. 10.1111/j.1464-5491.2010.02976.x.

  50. 50.

    Yang HC, Lin CH, Hsu CL, Hung SI, Wu JY, Pan WH, Chen YT, Fann CS: A comparison of major histocompatibility complex SNPs in Han Chinese residing in Taiwan and Caucasians. J Biomed Sci. 2006, 13: 489-498. 10.1007/s11373-006-9077-7.

  51. 51.

    Liu ZH, Guan TJ, Chen ZH, Li LS: Glucose transporter (GLUT1) allele (XbaI-) associated with nephropathy in non-insulin-dependent diabetes mellitus. Kidney Int. 1999, 55: 1843-1848. 10.1046/j.1523-1755.1999.00449.x.

  52. 52.

    Hsieh MC, Lin SR, Hsieh TJ, Hsu CH, Chen HC, Shin SJ, Tsai JH: Increased frequency of angiotensin-converting enzyme DD genotype in patients with type 2 diabetes in Taiwan. Nephrol Dial Transplant. 2000, 15: 1008-1013. 10.1093/ndt/15.7.1008.

  53. 53.

    Huo P, Zhang D, Guan X, Mei Y, Zheng H, Feng X: Association between genetic polymorphisms of ACE & eNOS and diabetic nephropathy.Mol Biol Rep 2014, [Epub ahead of print] doi:10.1007/s11033-014-3736-y.,

  54. 54.

    Liu L, Zheng T, Wang N, Wang F, Li M, Jiang J, Zhao R, Li L, Zhao W, Zhu Q, Jia W: The manganese superoxide dismutase Val16Ala polymorphism is associated with decreased risk of diabetic nephropathy in Chinese patients with type 2 diabetes. Mol Cell Biochem. 2009, 322: 87-91. 10.1007/s11010-008-9943-x.

  55. 55.

    Maeda S, Kobayashi MA, Araki S, Babazono T, Freedman BI, Bostrom MA, Cooke JN, Toyoda M, Umezono T, Tarnow L, Hansen T, Gaede P, Jorsal A, Ng DP, Ikeda M, Yanagimoto T, Tsunoda T, Unoki H, Kawai K, Imanishi M, Suzuki D, Shin HD, Park KS, Kashiwagi A, Iwamoto Y, Kaku K, Kawamori R, Parving HH, Bowden DW, Pedersen O, et al: A single nucleotide polymorphism within the acetyl-coenzyme A carboxylase beta gene is associated with proteinuria in patients with type 2 diabetes. PLoS Genet. 2010, 6: e1000842-10.1371/journal.pgen.1000842.

  56. 56.

    Liu L, Zheng T, Wang F, Wang N, Song Y, Li M, Li L, Jiang J, Zhao W: Pro12Ala polymorphism in the PPARG gene contributes to the development of diabetic nephropathy in Chinese type 2 diabetic patients. Diabetes Care. 2010, 33: 144-149. 10.2337/dc09-1258.

  57. 57.

    Tang SC, Leung VT, Chan LY, Wong SS, Chu DW, Leung JC, Ho YW, Lai KN, Ma L, Elbein SC, Bowden DW, Hicks PJ, Comeau ME, Langefeld CD, Freedman BI: The acetyl-coenzyme A carboxylase beta (ACACB) gene is associated with nephropathy in Chinese patients with type 2 diabetes. Nephrol Dial Transplant. 2010, 25: 3931-3934. 10.1093/ndt/gfq303.

  58. 58.

    Zhang H, Wang J, Yi B, Zhao Y, Liu Y, Zhang K, Cai X, Sun J, Huang L, Liao Q: BsmI polymorphisms in vitamin D receptor gene are associated with diabetic nephropathy in type 2 diabetes in the Han Chinese population. Gene. 2012, 495: 183-188. 10.1016/j.gene.2011.12.049.

  59. 59.

    Zhou J, Peng R, Li T, Luo X, Peng H, Zha H, Yin P, Wen L, Zhang Z: A potentially functional polymorphism in the regulatory region of let-7a-2 is associated with an increased risk for diabetic nephropathy. Gene. 2013, 527: 456-461. 10.1016/j.gene.2013.06.088.

  60. 60.

    Wu HY, Wang Y, Chen M, Zhang X, Wang D, Pan Y, Li L, Liu D, Dai XM: Association of ELMO1 gene polymorphisms with diabetic nephropathy in Chinese population. J Endocrinol Invest. 2013, 36: 298-302.

  61. 61.

    Schmid H, Boucherot A, Yasuda Y, Henger A, Brunner B, Eichinger F, Nitsche A, Kiss E, Bleich M, Grone HJ, Nelson PJ, Schlöndorff D, Cohen CD, Kretzler M: Modular activation of nuclear factor-kappaB transcriptional programs in human diabetic nephropathy. Diabetes. 2006, 55: 2993-3003. 10.2337/db06-0477.

  62. 62.

    Woroniecka KI, Park AS, Mohtat D, Thomas DB, Pullman JM, Susztak K: Transcriptome analysis of human diabetic kidney disease. Diabetes. 2011, 60: 2354-2369. 10.2337/db10-1181.

  63. 63.

    Hodgin JB, Nair V, Zhang H, Randolph A, Harris RC, Nelson RG, Weil EJ, Cavalcoli JD, Patel JM, Brosius FC, Kretzler M: Identification of cross-species shared transcriptional networks of diabetic nephropathy in human and mouse glomeruli. Diabetes. 2013, 62: 299-308. 10.2337/db11-1667.

Download references


We would like to thank all of the patients and study personnel who participated in this study. We also thank the Academia Sinica for technical support. This study was supported by the Academia Sinica (BM102010130), Ministry of Science and National Science Council Technology of Taiwan (NSC101-2314-B-039-017-MY3) and Ministry of Health & Welfare Clinical Trial and Research Center of Excellence (MOHW103-TDU-B-212- 113002).

Author information

Correspondence to Tsai-Chung Li or Fuu-Jen Tsai.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

LNL, CCC, CCL, CTC, SLRK, JHH, TCL, and FJT carried out the conception and design, acquisition of data, and analysis and interpretation of data. LNL, FYW, and TCL drafted the article or revised it critically for important intellectual content. All authors read and approved the final manuscript.

Electronic supplementary material


Additional file 1: Figure S1.: Regional plot of these 11 SNPs located on chromosome 16q22.1. The -log10 (p-value) (left y-axis) was from the trend test using our data. The estimated r2 and recombination rate (right y-axis) based on the HapMap Phase II JPT + CHB populations were plotted to reflect the LD structure. The gene information was from the UCSC (Build hg18). The regional plot was plotted from the LocusZoom, a web-based plotting tool ( Table S1. Summary of identified polymorphisms through GWAS. Table S2. Summary of variants associated with nephropathy in Han Chinese patients with type 2 diabetes through candidate gene approach. Table S3. is-rSNP prediction resultsa. Table S4. Gene expression from human diabetic nephropathy datasets in the Nephromine database. Table S5. Gene expression from Hodgin Diabetes Mouse dataseta in the Nephromine database [51]-[63]. (DOCX 111 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Liao, L., Chen, C., Wu, F. et al. Identified single-nucleotide polymorphisms and haplotypes at 16q22.1 increase diabetic nephropathy risk in Han Chinese population. BMC Genet 15, 113 (2014).

Download citation


  • Diabetic nephropathy
  • Single-nucleotide polymorphism
  • Haplotype
  • Han Chinese