Skip to main content

Nidogen-1 could play a role in diabetic kidney disease development in type 2 diabetes: a genome-wide association meta-analysis

Abstract

Background

Diabetic kidney disease (DKD) affects about 40% of patients with diabetes. It is incurable and usually leads to end-stage renal disease (ESRD). The pathogenesis of DKD is still not fully understood, and the genetics of DKD have not yet been extensively studied. In this study, we investigate the genetic basis of DKD in type 2 diabetes (T2D) to provide more insights into the pathogenesis of the disease.

Results

Using the data provided by the UK Biobank (UKBB), we performed a DKD genome-wide association study (GWAS) in 13,123 individuals with T2D as well as two creatinine estimated glomerular filtration rate (eGFR) GWA studies: one in 26,786 individuals with T2D and the other in 339,080 non-diabetic individuals. We also conducted a DKD GWAS meta-analysis combining our results with those published by the surrogate markers for micro- and macro-vascular hard endpoints for Innovative diabetes Tools (SUMMIT) consortium. We confirm two loci previously reported to be associated with chronic kidney disease (CKD) and eGFR in T2D. The UMOD-PDILT locus is associated with DKD (P = 1.17E−09) as well as creatinine eGFR in both people with T2D (P = 1.31E−15) and people without diabetes (P = 3.95E−73). The PRKAG2 locus is associated with creatinine eGFR in people with (P = 2.78E−10) and without (P = 5.65E−72) T2D. Our meta-analysis reveals a novel association between DKD and variant rs72763500 (chr1:236116561) which is a splicing quantitative trait locus (sQTL) for nidogen-1 (NID1) gene.

Conclusion

Our data confirm two loci previously reported in association with CKD and creatinine eGFR in T2D. It also suggests that NID1, a major component of the renal tubular basement membrane, could play a role in DKD development in T2D. While our NID1 finding remains to be replicated, it is a step toward a more comprehensive understanding of DKD pathogenesis.

Background

Diabetes causes CKD in about 40% of patients with diabetes in the USA and is the leading cause of ESRD [1]. While DKD usually occurs after 10 years in type 1 diabetes (T1D), it can manifest at the time of diagnosis in T2D [2]. DKD is characterized by glomerular hyperfiltration at early stages and then progress to albuminuria, reduced eGFR, and subsequent ESRD [3]. These functional changes are accompanied by structural changes which include thickening of the glomerular basement membrane as well as capillary and tubular basement membranes [4].

Not all patients with diabetes develop DKD, and early heritability studies showed that DKD happens in family clusters [5, 6] which indicates genetic susceptibility. While the heritability of DKD is estimated to be 35% in T1D [7], it is only 8% in T2D [8] which could be due to the highly heterogeneous phenotype in T2D [9]. Although the genetics of T2D have been intensively explored by many GWA studies, the genetics of its complications have not received as much attention from researchers. In particular, the genetic of DKD is much less studied in T2D compared to DKD in T1D and most performed studies are of small sample sizes.

In 2003, the very first DKD GWAS in T2D was conducted using gene-based single nucleotide polymorphisms (SNPs) [10]. It initially identified the SLC12A3 gene as a candidate gene for DKD in a group of 188 (94 cases and 94 controls) Japanese individuals. The results were replicated in a larger group of 732 individuals (466 cases and 266 controls). In 2005, another GWAS using gene-based SNPs identified ELMO1 as a DKD-associated gene in a discovery cohort of 179 (87 cases and 92 controls) and a replication group of 701 (459 cases and 242 controls) Japanese individuals [11].

In 2009, a larger GWAS including 359,193 autosomal SNPs was performed in 820 cases and 885 controls of people with T1D identified four loci associated with DKD near FRMD3, CHN2-CPVL, CARS, and MYO16-IRS2 genes [12, 13]. These findings were replicated in African-American (965 cases and 1029 controls) [14], European (743 cases and 646 controls) [15], and in Japanese [16] cohorts with T2D.

In 2018, the SUMMIT consortium conducted the largest DKD GWAS [8, 17]. The study included seven phenotypes and a total of 26,827 individuals with T2D (12,710 with DKD). The main DKD phenotype included 5,717 individuals (3,345 DKD cases and 2,372 controls) with T2D. The GWAS identified a novel variant near GABRR1 gene associated with microalbuminuria in T2D.

One of the core advantages of conducting GWA studies is gaining insights into the pathophysiology of human diseases. Despite the effort to uncover the genetic basis of DKD, the heterogeneous phenotype and the complexity of the disease made the GWA studies challenging and limited their replicability [9, 18]. Up to now, these factors have restricted our full understanding of the pathogenesis of DKD.

Here, we conducted a GWAS for DKD in Europeans with T2D using data from the UKBB. To further explore the results of the DKD GWAS, we performed a GWAS for creatinine eGFR measurement in both people with T2D and people who have not been diagnosed with T1D or T2D. Furthermore, we used the summary statistics released by the SUMMIT consortium [8, 17] to perform a DKD meta-analysis study. While our study confirms the previously published findings in both CKD and creatinine eGFR in T2D, it also uncovers a novel association between NID1 and DKD.

Methods

DKD phenotype definition and study population

Clinically, a patient will be diagnosed with CKD if either reduced eGFR (GFR < 60 ml/min/1.73 m2) or increased albuminuria defined by urine albumin–creatinine ratio (ACR) ≥ 30 mg/g is present for more than 3 month [19].

These diagnostic criteria do not differentiate DKD from CKD unrelated to T2D. Despite being rarely performed in practice, a kidney biopsy is the gold standard for determining the cause of CKD [9, 20]. These factors make it challenging to define the phenotype of DKD, but considering the timing of a CKD diagnosis relative to the T2D diagnosis can help identify DKD. As described below, we adopted a CKD in T2D vs long duration T2D without CKD definition in this study leading to a diabetes duration that is significantly (P = 2.20E−16) shorter in DKD cases compared to complication-free T2D controls (Table 1).

Table 1 DKD GWAS sample characteristics

We used the International Classification of Diseases 10th revision (ICD-10) to identify the study’s cohort from the UKBB. We selected participants with any T2D ICD-10 code (E11.0, E11.1, E11.2, E11.3, E11.4, E11.5, E11.6, E11.7, E11.8, E11.9). Out of 46,308 participants with T2D, we included 13,123 (7982 controls and 5141 cases) unrelated White Europeans in our GWAS. Table 1 summarizes the characteristics of the study’s cohort.

We considered three criteria to define DKD: (1) participants who were assigned the ICD-10 code E11.2 (type 2 diabetes mellitus with kidney complications), (2) participants who were assigned any CKD ICD-10 code (N18.0, N18.1, N18.2, N18.3, N18.4, N18.5, N18.8, N18.9) after being diagnosed with T2D, and (3) participants with GFR < 60 ml/min/1.73 m2 and or ACR ≥ 30 mg/g after being diagnosed with T2D. Participants with T2D for more than 10 years and who have not developed any diabetes-related complications were included as controls.

eGFR calculation

eGFR was calculated using the Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) creatinine equation:

$$\begin{aligned}{\text{eGFRcr}} &= 142 \times \min \left( {{\text{Scr}}/\upkappa ,1} \right)\alpha \times \max \left( {{\text{Scr}}/\upkappa ,1} \right) \\ &\quad- 1.200 \times 0.9938\;{\text{Age}} \times 1.012\left[ {\text{if female}} \right] \end{aligned}$$

where Scr = serum creatinine in mg/dL, κ = 0.7 (females) or 0.9 (males), α =  − 0.241 (female) or − 0.302 (male), min(Scr/κ, 1) is the minimum of Scr/κ or 1.0, max(Scr/κ, 1) is the maximum of Scr/κ or 1.0, and age (years) [21].

GWAS quality control

SNPs and individuals with genotyping rate < 98% were excluded. SNPs with minor allele frequency < 5% were excluded. SNPs deviating from Hardy–Weinberg equilibrium test (P value< 1E−10 in cases and < 1E−6 in controls) were removed. Individuals who deviate > 3 SD from the samples mean heterozygosity rate were removed. Principal component analysis (PCA) was performed to correct for population substructure. Only White European participants were selected for further analysis. The quality control steps were performed using Plink V 1.9 [22]. The PCA for all studies was performed using bigsnpr [23]. Related individuals were identified using KING [24], and individuals with kinship coefficient ≥ 0.17 were excluded from the analysis.

Statistical analysis

Genome-wide association analyses were performed using an additive model. The DKD GWAS was corrected for age, sex, duration of diabetes, and the top 10 principal components. The same covariates were adjusted for in the eGFR GWA studies except for the duration of diabetes. The analyses were conducted using Plink V 1.9 [22].

DKD meta-analysis

We conducted an inverse-variance fixed effects meta-analysis including the summary statistics from our DKD GWAS and the above-mentioned SUMMIT study. The meta-analysis included a total of 18,840 (11,327 cases and 7513 controls) European individuals with T2D; 5717 (3345 cases and 2372 controls) from the SUMMIT study and 13,123 (5141 cases and 7982 controls) from our DKD GWAS. The meta-analysis was performed in Plink V 1.9 [22] combining the effect sizes from both studies. The meta-analysis results included only variants present in both summary statistics reports.

Results

Heritability estimation

Heritability (h2) was estimated using the LDAK model [25] from 241,883 directly genotyped SNPs. For DKD, h2 = 0.027 with standard deviation (SD) = 0.03. For eGFR in non-diabetic individuals, h2 = 0.036 with SD = 0.001. For eGFR in T2D, h2 = 0.1 with SD = 0.01.

DKD GWAS

In the DKD GWAS, only one locus reached genome-wide significance: UMOD-PDILT (Uromodulin-Protein disulfide-isomerase-like protein of the testis) which is located on chromosome 16 (Table 2; Fig. 1a). The same locus was previously reported to be associated with CKD and reduced renal function represented by both serum creatinine and cystatin c eGFR in diabetes [26, 27]. Genomic inflation factor was estimated based on median chi-square statistics, and λGC was 1.0 (See QQ plot; Additional file 1: Fig. S1A).

Table 2 Genome-wide significant associations with eGFR, DKD, and DKD meta-analysis
Fig. 1
figure 1

GWAS plots showing a diabetic kidney disease (DKD). b eGFR measurement in T2D. c eGFR measurement in non-diabetic (no T1D or T2D) individuals. The dotted line indicates the genome-wide significance threshold at P = 5 × 10−8

eGFR GWA studies

To further investigate the UMOD association with renal function, we performed two eGFR GWA studies in two unrelated White European groups selected from UKBB: one in 26,786 individuals (16,323 males and 10,463 females) with T2D and another 339,080 in individuals (186,027 females and 153,053 males) who have not developed T2D or T1D. Both GWA studies were corrected for sex and age at the time of creatinine measurement.

The calculation of eGFR is detailed in the methods section. Genomic inflation factors were calculated based on median chi-square statistics, and λGC values were 1.02, and 1.4 for eGFR in T2D GWAS, and eGFR in non-diabetic individuals, respectively (See QQ plots in Additional file 1: Fig. S1b and c). We applied genomic control correction to the latter, and the adjusted P values are presented in Table 2, the Manhattan plot (Fig. 1c).

In the T2D eGFR GWAS, 2 loci reached genome-wide significance: UMOD-PDILT and PRKAG2 (Table 2; Fig. 1b). These two loci were previously identified to be associated with CKD, creatinine eGFR, and cystatin c eGFR [27]. The UMOD-PDILT locus was reported in association with creatinine eGFR in T2D [8] and both T1D and T2D combined [8, 27, 28]. The PRKAG2 loci were identified in association with creatinine eGFR in T2D [8].

The same two loci, among many, showed genome-wide significance in the non-diabetic eGFR GWAS (Table 2; Fig. 1c). The difference in the number of the significant signals in the eGFR T2D vs non-diabetic GWA studies could be attributed to the much larger sample size of the non-diabetic GWAS.

The UMOD-PDILT loci have been reported to be associated with the rapid decline in eGFR [29] which could explain its significant association with both eGFR and DKD across all performed GWA studies.

DKD meta-analysis

As mentioned previously, the SUMMIT consortium released summary statistics for a meta-analysis of seven DKD phenotypes. The overall DKD was defined by microalbuminuria, late DKD (albumin excretion rate ≥ 200 µg/min), or ESRD (eGFR < 15 mL/min/1.73 m2 or kidney dialysis or transplant). The study comprised a total of 3345 cases and 2372 controls of European descendent and estimated the heritability of DKD in T2D to be 8%. We chose to include the SUMMIT consortium summary data in our meta-analysis because it is one of the largest DKD GWA studies to date. In addition, it matches both our study’s phenotype definition and population ancestry.

The UMOD locus remained significant in the meta-analysis as in the other three GWA studies. We also identified a novel variant, rs72763500 (chr1:236116561), that have not been reported to be associated with DKD before (Table 2; Fig. 2). This novel variant is unique to DKD and did not show a significant association with the eGFR measurement. We queried the NHGRI GWAS catalog (www.genome.gov) to investigate the association of variant rs72763500 with related phenotypes, but it has not been reported by any study. The variant is reported by the Genotype-Tissue Expression (GTEx) project (https://gtexportal.org/home/snp/rs72763500) as a splicing quantitative trait locus (sQTL) of NID1 gene.

Fig. 2
figure 2

GWAS plot showing: diabetic kidney disease (DKD) meta-analysis. The dotted line indicates the genome-wide significance threshold at P = 5 × 10−8

We used the Functional Mapping and Annotation of GWA Studies (FUMA GWAS) platform [30] for post-GWAS annotation of our results. To identify the enriched differentially expressed genes (DEG) in the provided GWAS summary statistics, FUMA employs two-sided t tests per gene per tissue against all other tissues. Although the tissue enrichment analysis using the GTEx v8 dataset showed no significant (P-value ≤ 0.05 after Bonferroni correction) DEG sets, it showed that kidney medulla is the tissue with highest enrichment in our data (Additional file 1: Fig. S2).

Additional file 1: Fig. S3 shows the genes in our data that have been reported by the GWAS Catalog. The related traits to our data are decline in kidney function traits and chronic kidney disease.

Discussion

For two decades, GWA studies have contributed to our understanding of the biological and physiological aspects of many human traits and diseases. They are powerful tools that have the potential to decipher the pathophysiology of many diseases as well as uncover new drug targets.

DKD is uncurable and usually leads to a debilitating ESRD which greatly reduces the patients’ quality of life. As we mentioned earlier, the genetics of DKD have not yet been extensively studied and, to date, the DKD GWA studies are still limited by the small sample sizes. While moderately extending the sample size of available studies, our study confirms the previously reported variants to be associated with CKD, eGFR in T2D, and reduced renal function. It also reports a novel association between NID1 gene that may be specific to DKD.

NID1 is a sulfated glycoprotein ubiquitous basement membrane protein and is involved in the glomerular basement membrane (GBM) development. We queried the Human Protein Atlas single-cell RNA expression database (https://www.proteinatlas.org/ENSG00000116962-NID1/single+cell+type) for NID1 tissue enrichment. NID1 is highly expressed in mesenchymal cells in adipose tissue, breast, esophagus, endometrium, liver, skeletal, and heart muscles. In the kidney, it is more abundant in the proximal tubular cells of the nephron in comparison to the distal tubular cells and collecting duct cells. On the contrary, UMOD is more abundant in the distal tubular cells and collecting duct cells compared to the proximal tubular cells.

The GBM is mainly composed of laminin, type IV collagen, nidogens (NID1 and NID2), and heparan sulfate proteoglycan [31, 32]. NID1 binds both laminin and type IV collagen [33]. Deletion of one of the nidogens in mice does not affect the development of basement membranes, but deletion of both genes leads to basement membrane abnormalities in the lungs and heart and subsequently perinatal death. The nidogen double mutant mice developed normal glomerular basement membrane, but the tubular basement membrane was sometimes absent or altered [34]. The effect of nidogen mutation on the human kidney is still largely unknown.

Our meta-analysis finding, although biologically relevant, is limited by the lack of a replication study. More studies are needed to substantiate and further investigate the role of NID1 in human DKD.

Availability of data and materials

The datasets generated during and/or analyzed during the current study are available from the corresponding author upon reasonable request.

Abbreviations

CKD:

Chronic kidney disease

DKD:

Diabetic kidney disease

eGFR:

Estimated glomerular filtration rate

ESKD:

End-stage kidney disease

FUMA GWAS:

Functional mapping and annotation of genome-wide association studies

GTEx:

Genotype-tissue expression

GWA:

Genome-wide association

GWAS:

Genome-wide association study

NID1:

Nidogen-1

PCA:

Principal component analysis

SNP:

Single nucleotide polymorphism

sQTL:

Splicing quantitative trait locus

SUMMIT:

Surrogate markers for micro- and macro-vascular hard endpoints for innovative diabetes tools.

T1D:

Type 1 diabetes

T2D:

Type 2 diabetes

UKBB:

UK Biobank

References

  1. Centers for Disease Control and Prevention. National Diabetes Statistics Report website. https://www.cdc.gov/diabetes/data/statistics-report/index.html. Accessed May 2022.

  2. Alicic RZ, Rooney MT, Tuttle KR. Diabetic kidney disease. Clin J Am Soc Nephrol. 2017;12:2032–45. https://doi.org/10.2215/cjn.11491116.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Mogensen CE, Christensen CK, Vittinghus E. The stages in diabetic renal disease: with emphasis on the stage of incipient diabetic nephropathy. Diabetes. 1983;32:64–78. https://doi.org/10.2337/diab.32.2.s64.

    Article  PubMed  Google Scholar 

  4. Reidy K, Kang HM, Hostetter T, Susztak K. Molecular mechanisms of diabetic kidney disease. J Clin Investig. 2014;124:2333–40. https://doi.org/10.1172/jci72271.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Seaquist ER, Goetz FC, Rich S, Barbosa J. Familial clustering of diabetic kidney disease. N Engl J Med. 1989;320(18):1161–5. https://doi.org/10.1056/NEJM198905043201801.

    Article  CAS  PubMed  Google Scholar 

  6. Freedman BI, Tuttle AB, Spray BJ. Familial predisposition to nephropathy in African-Americans with non-insulin-dependent diabetes mellitus. Am J Kidney Dis. 1995;25(5):710–3. https://doi.org/10.1016/0272-6386(95)90546-4.

    Article  CAS  PubMed  Google Scholar 

  7. Sandholm N, et al. The genetic landscape of renal complications in Type 1 diabetes. J Am Soc Nephrol. 2017;28:557–74. https://doi.org/10.1681/asn.2016020231.

    Article  PubMed  Google Scholar 

  8. Zuydam V, et al. A genome-wide association study of diabetic kidney disease in subjects with type 2 diabetes. Diabetes. 2018;67:1414–27. https://doi.org/10.2337/db17-0914.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Cole JB, Florez JC. Genetics of diabetes mellitus and diabetes complications. Nat Rev Nephrol. 2020;16:377–90. https://doi.org/10.1038/s41581-020-0278-5.

    Article  PubMed  Google Scholar 

  10. Tanaka N, Babazono T, Saito S, Sekine A, Tsunoda T, Haneda M, Tanaka Y, Fujioka T, Kaku K, Kawamori R, Kikkawa R, Iwamoto Y, Nakamura Y, Maeda S. Association of solute carrier family 12 (sodium/chloride) member 3 with diabetic nephropathy, identified by genome-wide analyses of single nucleotide polymorphisms. Diabetes. 2003;52:2848–53. https://doi.org/10.2337/diabetes.52.11.2848.

    Article  CAS  PubMed  Google Scholar 

  11. 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–8. https://doi.org/10.2337/diabetes.54.4.1171.

    Article  CAS  PubMed  Google Scholar 

  12. Pezzolesi MG, Poznik GD, Mychaleckyj JC, Paterson AD, Barati MT, Klein JB, Ng DPK, 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–10. https://doi.org/10.2337/db08-1514.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Pezzolesi MG, Skupien J, Mychaleckyj JC, Warram JH, Krolewski AS. Insights to the Genetics of Diabetic Nephropathy Through a Genome-Wide Association Study of the GoKinD Collection. Semin Nephrol. 2010;30(2):126–40. https://doi.org/10.1016/j.semnephrol.2010.01.004.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Palmer ND, Ng MCY, Hicks PJ, Mudgal P, Langefeld CD, Freedman BI, Bowden DW. Evaluation of candidate nephropathy susceptibility genes in a genome-wide association study of African American diabetic kidney disease. PLOS ONE. 2014;9:e88273. https://doi.org/10.1371/journal.pone.0088273.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Pezzolesi MG, Poznik GD, Skupien J, Smiles AM, Mychaleckyj JC, Rich SS, Warram JH, Krolewski AS. An intergenic region on chromosome 13q33.3 is associated with the susceptibility to kidney disease in type 1 and 2 diabetes. Kidney Int. 2011;80:105–11. https://doi.org/10.1038/ki.2011.64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Maeda S, Araki S-I, Babazono T, Toyoda M, Umezono T, Kawai K, Imanishi M, Uzu T, Watada H, Suzuki D, Kashiwagi A, Iwamoto Y, Kaku K, Kawamori R, Nakamura Y. Replication study for the association between four loci identified by a genome-wide association study on European American subjects with type 1 diabetes and susceptibility to diabetic nephropathy in Japanese subjects with type 2 diabetes. Diabetes. 2010;59:2075–9. https://doi.org/10.2337/db10-0067.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Data on diabetic kidney disease has been contributed by the SUMMIT Consortium and has been downloaded from mccarthy.well.ox.ac.uk

  18. Ahlqvist E, Van Zuydam NR, Groop LC, Mccarthy MI. The genetics of diabetic complications. Nat Rev Nephrol. 2015;11:277–87. https://doi.org/10.1038/nrneph.2015.37.

    Article  CAS  PubMed  Google Scholar 

  19. What is the criteria for CKD. National Kidney Foundation. (2015, April 20). Retrieved June 7, 2022, from https://www.kidney.org/professionals/explore-your-knowledge/what-is-the-criteria-for-ckd

  20. Anders H-J, Huber TB, Isermann B, Schiffer M. CKD in diabetes: diabetic kidney disease versus nondiabetic kidney disease. Nat Rev Nephrol. 2018;14:361–77. https://doi.org/10.1038/s41581-018-0001-y.

    Article  CAS  PubMed  Google Scholar 

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

    Article  Google Scholar 

  22. Chang CC, Chow CC, Tellier LC, Vattikuti S, Purcell SM, Lee JJ. Second-generation PLINK: rising to the challenge of larger and richer datasets. GigaScience. 2015. https://doi.org/10.1186/s13742-015-0047-8.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Privé F, Aschard H, Ziyatdinov A, Blum MGB. Efficient analysis of large-scale genome-wide data with two R packages: bigstatsr and bigsnpr. Bioinformatics. 2018;34:2781–7. https://doi.org/10.1093/bioinformatics/bty185.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Manichaikul A, Mychaleckyj JC, Rich SS, Daly K, Sale M, Chen W-M. Robust relationship inference in genome-wide association studies. Bioinformatics. 2010;26:2867–73. https://doi.org/10.1093/bioinformatics/btq559.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Speed D, Balding DJ. SumHer better estimates the SNP heritability of complex traits from summary statistics. Nat Genet. 2019;51:277–84. https://doi.org/10.1038/s41588-018-0279-5.

    Article  CAS  PubMed  Google Scholar 

  26. Köttgen A, et al. New loci associated with kidney function and chronic kidney disease. Nat Genet. 2010;42:376–84. https://doi.org/10.1038/ng.568.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Pattaro C, et al. Genetic associations at 53 loci highlight cell types and biological pathways relevant for kidney function. Nat Commun. 2016;7:10023. https://doi.org/10.1038/ncomms10023.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Hellwege JN, Velez Edwards DR, Giri A, Qiu C, Park J, Torstenson ES, Keaton JM, Wilson OD, Robinson-Cohen C, Chung CP, Roumie CL, Klarin D, Damrauer SM, Duvall SL, Siew E, Akwo EA, Wuttke M, Gorski M, Li M, Li Y, Gaziano JM, Wilson PWF, Tsao PS, O’Donnell CJ, Kovesdy CP, Pattaro C, Köttgen A, Susztak K, Edwards TL, Hung AM. Mapping eGFR loci to the renal transcriptome and phenome in the VA million veteran program. Nat Commun. 2019. https://doi.org/10.1038/s41467-019-11704-w.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Gorski M, Jung B, Li Y, Matias-Garcia PR, et al. Meta-analysis uncovers genome-wide significant variants for rapid kidney function decline. Kidney Int. 2021;99:926–39. https://doi.org/10.1016/j.kint.2020.09.030.

    Article  CAS  PubMed  Google Scholar 

  30. Watanabe K, Taskesen E, Van Bochoven A, Posthuma D. Functional mapping and annotation of genetic associations with FUMA. Nat Commun. 2017. https://doi.org/10.1038/s41467-017-01261-5.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Miner JH. The glomerular basement membrane. Exp Cell Res. 2012;318:973–8. https://doi.org/10.1016/j.yexcr.2012.02.031.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Miner JH. Renal basement membrane components. Kidney Int. 1999;56:2016–24. https://doi.org/10.1046/j.1523-1755.1999.00785.x.

    Article  CAS  PubMed  Google Scholar 

  33. Dziadek M. Role of laminin-Nidogen complexes in basement membrane formation during embryonic development. Experientia. 1995;51(9–10):901–13. https://doi.org/10.1007/bf01921740.

    Article  CAS  PubMed  Google Scholar 

  34. Bader BL, Smyth N, Nedbal S, Miosge N, Baranowsky A, Mokkapati S, Murshed M, Nischt R. Compound genetic ablation of Nidogen 1 and 2 causes basement membrane defects and perinatal lethality in mice. Mol Cell Biol. 2005;25:6846–56. https://doi.org/10.1128/mcb.25.15.6846-6856.2005.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Author information

Authors and Affiliations

Authors

Contributions

AK contributed to conceptualization, methodology, data curation, visualization, and writing—original draft preparation. AT was involved in supervision and writing—reviewing and editing. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Ali Torkamani.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

We report no conflicts of interest, for all listed authors, relevant to this article.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1

. Q–Q plots: (A) Diabetic Kidney Disease (DKD). (B) eGFR measurement in T2D. (C) eGFR measurement in non-diabetic (no T1D or T2D) individuals. λ = genomic inflation factor. Figure S2 Tissue enrichment analysis of the diabetic kidney disease (DKD) meta-analysis showing the levels of both up-regulated and down-regulated differentially expressed genes (DEG) in various tissues. Significant enrichment level is Bonferroni corrected P-value ≤ 0.05. Figure S3 GWAS Catalog reported genes: showing genes, from the diabetic kidney disease (DKD) meta-analysis, reported by the GWAS Catalog.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Khattab, A., Torkamani, A. Nidogen-1 could play a role in diabetic kidney disease development in type 2 diabetes: a genome-wide association meta-analysis. Hum Genomics 16, 47 (2022). https://doi.org/10.1186/s40246-022-00422-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40246-022-00422-y