Genome-wide meta-analysis and omics integration identifies novel genes associated with diabetic kidney disease
==============================================================================================================

* Niina Sandholm
* Joanne B Cole
* Viji Nair
* Xin Sheng
* Hongbo Liu
* Emma Ahlqvist
* Natalie van Zuydam
* Emma H Dahlström
* Damian Fermin
* Laura J Smyth
* Rany M Salem
* Carol Forsblom
* Erkka Valo
* Valma Harjutsalo
* Eoin P Brennan
* Gareth McKay
* Darrell Andrews
* Ross Doyle
* Helen C Looker
* Robert G Nelson
* Colin Palmer
* Amy Jayne McKnight
* Catherine Godson
* Alexander P Maxwell
* Leif Groop
* Mark I McCarthy
* Matthias Kretzler
* Katalin Susztak
* Joel N Hirschhorn
* Jose C Florez
* Per-Henrik Groop
* for the GENIE Consortium

## Abstract

**Background** Diabetes is the leading cause of kidney disease, and heritability studies demonstrate a substantial, yet poorly understood, contribution of genetics to kidney complications in people with diabetes.

**Methods** We performed genome-wide association study (GWAS) meta-analyses using ten different phenotypic definitions of diabetic kidney disease (DKD), including nearly 27,000 individuals with diabetes, and integrated the results with various kidney omics datasets.

**Results** The meta-analysis identified a novel low frequency intronic variant (rs72831309) in the *TENM2* gene encoding teneurin transmembrane protein 2 associated with a lower risk of the combined chronic kidney disease (CKD; eGFR<60 ml/min/1.73 m2) and DKD (microalbuminuria or worse) phenotype (“CKD-DKD”, odds ratio 2.08, *p*=9.8×10−9). Gene-level analysis identified ten genes associated with DKD (*COL20A1, DCLK1, EIF4E, PTPRN-RESP18, GPR158, INIP-SNX30, LSM14A*, and *MFF, p*<2.7×10−6). Integration of GWAS data with human glomerular and tubular expression data in a transcriptome-wide association study demonstrated higher tubular *AKIRIN2* gene expression in DKD versus non-DKD controls (*p*=1.1×10−6). The lead SNPs within the *DCLK1, AKIRIN2, SNX30* and three other gene regions significantly alterated the methylation at this region in kidneys (*p*<2.2×10−11). Expression of target genes in kidney tubules or glomeruli correlated with relevant pathological phenotypes. For example, tubular *TENM2* expression positively correlated with eGFR (*p*=2.3×10−9) and negatively with tubulointerstitial fibrosis (*p*=4.7×10−9), tubular *DCLK1* expression positively correlated with fibrosis (*p*=1.6×10−12), and *SNX30* level positively correlated with eGFR (*p*=7.6×10−13) and negatively with fibrosis (*p*<2×10−16).

**Conclusions** GWAS meta-analysis and integration with renal omics data points to novel genes contributing to pathogenesis of DKD.

Keywords
*   Diabetic kidney disease
*   genome-wide association study
*   genetics
*   transcriptomics
*   diabetes complications

## INTRODUCTION

Diabetes is the leading cause of kidney disease. Diabetic kidney disease (DKD) is associated with high cardiovascular risk1 and mortality2, and consequently, both diabetes and kidney disease are among the most important causes of death worldwide3. While environmental factors, and especially blood glucose control, have a major impact on the risk of developing DKD, genetic factors also contribute to disease4,5. Although more than 300 genetic loci have been associated with chronic kidney disease (CKD) in the general population, these loci show limited effect in DKD, especially in individuals with type 1 diabetes (T1D)6. Genome-wide association studies (GWAS) have previously identified a handful of genetic loci for DKD in individuals with T1D at the genome-wide significance level (*p*<5×10−8)7–10. Recently, a GWAS meta-analysis including up to 19,406 individuals with T1D from the Diabetic Nephropathy Collaborative Research Initiative (DNCRI) identified 16 loci associated with various DKD definitions. The strongest association was a common missense variant on the *COL4A3* gene which also showed evidence of association in individuals with type 2 diabetes (T2D).6 A GWAS meta-analysis from the The SUrrogate markers for Micro- and Macrovascular hard endpoints for Innovative diabetes Tools (SUMMIT) consortium, including 6,000 individuals with T2D from five different studies, identified three loci for DKD in T2D, including the *UMOD* and *PRKAG2* loci previously identified in the general population11. However, meta-analysis with SUMMIT T1D studies (N=5,156)4 did not yield any genome-wide significant findings11. In this study, we performed GWAS meta-analyses on ten different DKD case-control definitions, including nearly 27,000 individuals with T1D or T2D from the two large consortia (including DNCRI6, SUMMIT-T1D4 and SUMMIT-T2D11 studies) to increase the power to detect novel genetic risk factors, followed by integration with versatile biological data to improve our understanding of the underlying biological mechanisms and clinical correlations (**Figure 1**).

![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/29/2021.08.27.21262264/F1.medium.gif)

[Figure 1:](http://medrxiv.org/content/early/2021/08/29/2021.08.27.21262264/F1)

Figure 1: Schematic illustration of the study design from GWAS meta-analysis to integration with various omics data sets.
GWAS meta-analysis for ten different phenotypic definitions of DKD included up to 26,785 individuals with either T1D or T2D from the previous DNCRI and SUMMIT GWAS meta-analyses. TWAS: Transcriptome-wide association study, integrating the GWAS meta-analysis results with kidney expression quantitative trait locus (eQTL) data for tubular and glomerular compartments, identifying genes with differential expression in DKD. mQTL: methylation quantitative trait locus (mQTL) data, identifying SNPs associated with DNA methylation at CpG sites. snATACseq: single nucleus ATAC sequencing, informative of chromatin openness in various kidney cell types. RegulomeDB: database with extensive epigenetic annotation for SNPs. PCHiC: promoter capture HiC sequencing data for sequence interaction with gene promoters, proposing target genes. Kidney transcriptomics: gene expression in glomerular and tubular tissue in nephrectomy samples, or in Pima Indian biopsies, correlated with various renal parameters.

## METHODS

### Participating studies and phenotype definitions

A total of ten case – control definitions for DKD were included in DNCRI6, based on either urinary albumin excretion rate (AER; divided into controls with normal AER, and cases with microalbuminuria, macroalbuminuria, or ESRD), estimated glomerular filtration rate (eGFR), or both, and harmonized to match and include all seven phenotypic definitions assessed in SUMMIT-T1D4 and SUMMIT-T2D11 analyses (**Supplemental Table 1**). All individuals (both cases and controls) had diabetes (either T1D or T2D) and cases had some form of kidney disease. The phenotypic comparisons are as follows: controls with normal AER vs. DKD cases with microalbuminuria or worse (“All vs. Ctrl”), macroalbuminuria or worse (“Severe DKD”), microalbuminuria alone (“Micro”), or “ESRD”; ESRD cases vs. everyone else (“ESRD vs. All”); controls with normal eGFR defined as eGFR ≥ 60 ml/min/1.73 m2 vs. CKD defined as eGFR < 60 ml/min/1.73 m2 (“CKD”); and “CKD-DKD” based on both AER and eGFR, with controls with normal AER and eGFR vs. cases with microalbuminuria or worse and eGFR < 45 ml/min/1.73 m2. For the three phenotypic comparisons not initially part of the SUMMIT analysis (normal AER vs. macroalbuminuria [“Macro”], ESRD vs. macroalbuminuria [“ESRD vs. macro”], and controls with eGFR ≥ 60 ml/min/1.73m2 vs. CKD cases with eGFR < 15 ml/min/1.73m2 or ESRD [“CKD extremes”]), GWAS and meta-analysis were performed with three SUMMIT T2D studies (Genetics of Diabetes Audit and Research, Tayside and Scotland [GoDARTS] 1 and 2, and Scania Diabetes Registry [SDR] T2D cohort) and the SDR T1D cohort. Individuals from the Finnish Diabetic Nephropathy Study (FinnDiane) were included in both the original DNCRI (N=6,019) and SUMMIT-T1D analyses (N=3,415), but for the purpose of this meta-analysis, FinnDiane was only included in the DNCRI meta-analysis and therefore excluded from the SUMMIT-T1D data (**Supplemental Table 2**). All contributing studies were performed in accordance with the Declaration of Helsinki and Declaration of Istanbul.

### Statistical analysis

Genotyping and statistical analysis of the DNCRI6 and SUMMIT4,11 cohorts have been previously described. Analysis plans were similar in both consortia, and the main characteristics are described in **Supplemental Table 3**. Imputation was performed using 1000Genomes Phase 3 reference panel in DNCRI, and the older 1000Genomes phase 1 panel in the SUMMIT cohorts. In both consortia, analyses were performed in unrelated individuals using the SNPtest additive score test, adjusting for age, sex, diabetes duration, the genetic principal components, and study-specific covariates (e.g., site or genotyping batch). Variants were filtered for INFO imputation quality score ≥ 0.3 (DNCRI) or ≥ 0.4 (SUMMIT) and minor allele count ≥ 10 in both cases and in controls. In SUMMIT, variants were further filtered to those with minor allele frequency (MAF) ≥ 0.01. Within-consortium meta-analyses were performed with inverse-variance fixed effects meta-analysis based on the effect size estimates. Meta-analyses between DNCRI, SUMMIT-T1D, and SUMMIT-T2D were performed with inverse-variance fixed effect methods based on the effect size estimates from the summary statistics for each of the three datasets. Finally, variants were limited to those found in at least two studies and reported in the 1000Genomes phase 3 reference panel. Study-wise summary statistics were available for DNCRI and SUMMIT-T1D studies. Regional association plots were plotted with LocusZoom.12

### Conditional analysis

We performed conditional analysis of the *COL4A3* locus with apparent secondary association peak using GCTA v1.93β 13 and FinnDiane GWAS data as the reference panel.

### Gene-level analysis

SNP summary statistics from the GWAS meta-analysis were aggregated by gene-level regression analysis using two related software programs, MAGMA14 and PASCAL15, using default parameters. Gene-level significance thresholds were determined by a Bonferroni multiple-testing correction based on the number of genes tested for each of the ten phenotypes within each software program (number of genes ranged from 18,439-21,790; significance thresholds ranged from 2.7×10−6 to 2.3×10−6).

### Transcriptome-wide association study (TWAS)

MetaXcan16 was applied to integrate kidney eQTL datasets with the GWAS meta-analysis results and to map disease-associated genes. The *cis-* eQTL data for microdissected human glomerular (N=119) and tubular (N=121) samples were obtained from Susztaklab Kidney Biobank ([https://susztaklab.com/eQTLci/download.php](https://susztaklab.com/eQTLci/download.php))17, and were analyzed jointly to infer differential gene expression in cases vs. controls using MetaXcan software with default parameters. The GTEx Elastic-Net Model pipeline ([https://github.com/hakyimlab/PredictDBPipeline](https://github.com/hakyimlab/PredictDBPipeline)) was applied to prepare the model used for MetaXcan. The LD references were estimated based on genotypes of European individuals from the 1000 Genome Project. Using FDR < 0.05, the method indentified 5,990 coding genes with significant models for glomerular eQTL, and 5,371 coding genes for tubular eQTL. Significant association was defined as *p*<0.05/2/6050=4.1×10−6, i.e., corrected for two tissues and 6,050 genes found in either tubular or glomerular eQTL data.

### Gene-prioritization analysis

Gene prioritization at each of our top loci was performed using two complementary similarity-based gene prioritization approaches (PoPS v0.118 and MAGMA14), which integrate GWAS summary statistics with gene set enrichment analysis based on a variety of biological annotation datasets including gene expression, curated pathways, protein-protein interactions, and mouse gene knock-out studies.

For PoPS gene prioritization, MAGMA is first used to calculate gene-level association statistics for 18,383 protein-coding genes in the genome, which is used to assess feature enrichment. PoPS then calculates polygenic priority (PoP) scores for each gene based on its membership to enriched features. For each of our top loci, we annotated the PoPS prioritized gene as the one with the highest PoP score within a 500kb flanking window of each of our lead SNPs.

Of note, the *PRNCR1* gene annotated as the nearest gene to SNP rs185299109 was not included in the PoPS protein-coding gene dataset and the CKD-associated SNP rs185299109 located in an intergenic region was also excluded from this analysis. MAGMA gene prioritization was conducted using a recently developed extension to the method as described and implemented in Benchmarker software19, enabling the explicit derivation of gene prioritization results from gene set enrichment analysis. Like the Benchmarker approach, we classified genes as members of each gene set using the top 50, 100, and 200 ranked genes, and obtained similar results from all three. To identify the PoPS features that contributed to the prioritization of *COL4A3*, we limited it to the selected marginal gene features (PoPS step 1), multiplied the *COL4A3* beta hats (PoPS step 2) by the *COL4A3* feature’s scores, and ranked the features by the highest overall score.

### Expression quantitative trait loci (eQTL)

eQTL associations were sought from the eQTLGen database for eQTL in whole blood from >30,000 participants ([http://www.eqtlgen.org/](http://www.eqtlgen.org/))20. Kidney specific eQTL associations were queried from eQTL datasets for glomeruli, tubules17, and a meta-analysis of four eQTL studies with 451 kidney samples. The meta-analysis of four eQTLs datasets obtained from the Susztak lab, The Cancer Genome Atlas (TCGA), the Genotype-Tissue Expression (GTEx v8), and the Nephrotic Syndrome Study Network (NephQTL)17,21–23, was performed using METAL with fixed effects inverse-variance meta-analysis24.

### Methylation quantitative trait loci (mQTL)

mQTL associations were sought for the lead SNPs in 188 healthy kidney samples (eGFR > 60 and fibrosis < 10%), with Bonferroni threshold (*p*<1.5×10−11) considered genome-wide significant. DNA methylation of CpG sites were profiled in 188 healthy kidney samples by the Infinium MethylationEPIC Kit and BeadChips (Illumina, USA) and were transformed by an inverse-normal transformation after quality control using SeSAMe.25 Genotypes for these samples were profiled by from Axiom Tx and Axiom Biobank arrays and imputed using the multiethnic panel reference from 1000Genomes Phase 3 (NCBI build 37, released in October 2014). The association between CpG site and the SNPs within 1Mb were estimated by linear regression model using MatrixQTL,26 with covariates including collection site, age, sex, top five genotype PCs, degree of bisulfite conversion, sample plate, and sentrix position and PEER factors. For the significant CpG sites, we then sought for evidence of association between blood methylation levels and eGFR, or eGFR decline, in 500 individuals with diabetes27; we furthermore tested association with DKD in our epigenome-wide association study in 1304 UK-ROI and FinnDiane participants, analyzed using the Infinium MethylationEPIC Kit and BeadChips (Illumina, USA), following the QC and analysis procedures described earlier for UK-ROI28. Meta-analysis of the two data sets was performed with METAL software24, based on *p*-values and direction of effect.

### Human kidney gene expression

For the 29 lead genes or transcripts underlying or located near the lead SNPs, or based on gene-level analyses, TWAS, PoPS, or kidney eQTL data, we studied gene expression in kidneys in human transcriptomics data from nephrectomy samples (433 tubule and 335 glomerulus samples)29 and kidney biopsies from the Pima Indian cohort (67 glomerular and 47 tubulointerstitial tissues),30 and tested for correlation with relevant pathological phenotypes. The microdissected nephrectomy samples were from individuals with varying degree of diabetic and hypertensive kidney disease, and gene expression was defined with RNA sequencing. Pearson correlation *p*-values below 2.2×10−4 were considered significant, corrected for multiple testing for 29 genes, two tissue compartments, and four phenotypes (eGFR, fibrosis, glomerulosclerosis and group comparison). The study was approved by the institutional review board of the University of Pennsylvania.

In the Pima Indian cohort, gene expression profiling in the first biopsy was performed with Affymetrix gene chip arrays30, and with Illumina RNA-sequencing for the second biopsy, as described earlier6. Available phenotypes included progression to ESRD, measured GFR (mGFR), albumin-to-creatinine ratio (ACR), glycated hemoglobin (HbA1c) and six kidney morphological parameters for both biopsies, and change in the phenotypes between the first and the second study biopsies (27 phenotypes in total).31 Pearson correlation *p*-values below 3.2×10−5were considered significant, corrected for 29 genes, 2 tissues, and 27 phenotypes; *p*-values below 8.6×10−4 (i.e. without correction for 27 phenotypes) were considered suggestive. The study was approved by the Institutional Review Board of the National Institute of Diabetes and Digestive and Kidney Diseases.

### Further annotation of the lead variants

Chromatin 3D conformation interactions with gene transcription start sites (TSS) were queried for the most significant SNPs from the promoter capture Hi-C (PCHiC) data from the [www.chicp.org](http://www.chicp.org) web interface, including data for GM12878 lymphoblastoid cell line and CD34 cells32, hESC derived cardiomyocytes33, 16 primary blood cells34, and pancreatic islets35. Interactions with score ≥ 5 were considered significant. We queried chromatin accessibility in kidney single-nucleus ATAC-sequencing (snATACseq) data available at [https://susztaklab.com/human_kidney/igv/](https://susztaklab.com/human_kidney/igv/) (accessed 24 June 2021)36. Detailed gene expression in kidney single cell RNA sequencing (scRNAseq) data was queried in the Human Diabetic Kidney data set (23,980 nuclei) by Wilson *et al*.37, accessed through [http://humphreyslab.com/SingleCell](http://humphreyslab.com/SingleCell). Further epigenetic annotation was sought from the regulomeDB38, and differential renal gene expression in DKD versus healthy controls from the Ju CKD39 and Woroniecka40 data sets in the NephroSeq portal ([www.nephroseq.org](http://www.nephroseq.org)). Of note, samples in the queried Woroniecka40 data are a subset (N=22) of the more recent RNA sequencing based nephrectomy samples mentioned above (N=433)29.

### LD score regression (LDSR) and Mendelian Randomization (MR)

LDSR41 was performed at LDhub ([http://ldsc.broadinstitute.org/](http://ldsc.broadinstitute.org/)) for 78 glycemic, autoimmune, anthropometric, bone, smoking behavior, lipid, kidney, uric acid, cardiometabolic, and aging related traits, based on the GWAS summary statistics of the ten DKD phenotypes explored. Variants were filtered to those with MAF ≥1%. LDSR associations with *p*<6.4×10−4 were defined significant after Bonferroni correction for 78 traits. To identify causal relationships for significant traits in the LDSR against DKD, we performed summary-based two-sample MR implemented in the R package TwoSampleMR42. For the SNP-trait associations, we selected genetic variants as instrumental variables (IV) that were independently associated with the selected traits (*p*<5×10−8; r2 < 0.001 based on the 1000Genomes EUR panel; LD window=10,000 kb) from published GWAS. Palindromic SNPs with intermediate allele frequency (MAF close to 50%) were removed. Traits with less than five IVs were excluded from the MR analysis. Primarily, we used inverse variance-weighted (IVW) regression, but causality was further assessed using methods less sensitive to pleiotropy/heterogeneity (weighted median and MR-Egger regression)43. Heterogeneity of SNP estimates in MR was assessed with Cochran’s Q statistic *p*-value and the I2 statistic. The MR–Egger intercept test was used to detect unbalanced horizontal pleiotropy.

### Data sharing information

The GWAS meta-analysis results can be accessed via the T1D, T2D, and Cardiometabolic Disease (CMD) Knowledge Portals ([https://t1d.hugeamp.org/](https://t1d.hugeamp.org/); [https://t2d.hugeamp.org/](https://t2d.hugeamp.org/); [https://hugeamp.org/](https://hugeamp.org/)), and downloaded on their respective downloads pages.

## RESULTS

### GWAS meta-analysis

The GWAS meta-analysis of the DNCRI (T1D), and SUMMIT-T1D and SUMMIT-T2D cohorts included up to 26,785 individuals with diabetes from 25 studies; 11,380 individuals with any DKD (micro- or macroalbuminuria or ESRD) and 15,405 individuals with normal AER (**Supplemental Table 2**). QQ plots and genomic correlation λGC of the meta-analysis indicated no marked inflation of the results; similarly, also LDSR intercepts were close to 1.00 (1.00 – 1.05) for each phenotype, demonstrating no evidence of population stratification bias in the GWAS meta-analysis (**Supplemental Figure 1)**.

The meta-analysis identified a novel association between the combined CKD-DKD phenotype and rs72831309 (minor allele frequency [MAF] =4%, odds ratio [OR] = 2.08, 95% confidence interval [CI] 1.62 - 2.67, *p*=9.8×10−9; **Figure 2A**; **Table 1; Supplemental Table 4**); the variant is located in an intron of the *TENM2* gene encoding the teneurin transmembrane protein 2. Exploration of the RegulomeDB revealed that rs72831309 alters a predicted CREB1 transcription factor binding site (**Figure 2C**). Kidney eQTL data indicated that rs72831309 was nominally associated with expression of a *TENM2* antisense transcript *CTB-178M22*.*2* in kidneys (*p*=6.9×10−3). Furthermore, chromatin conformation data in the GM12878 cell line indicated that the rs72831309 containing DNA fragment interacts with the *TENM2* transcription start site (TSS), as well as with three antisense transcripts (*CTB-180C19*.*1, CTB-105L4*.*2*, and *CTB-78F1*.*1*) within the *TENM2* gene32. Finally, human kidney scRNAseq showed that *TENM2* is expressed in podocytes and in the parietal epithelial cells (**Figure 2D**)37.

![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/29/2021.08.27.21262264/F2.medium.gif)

[Figure 2:](http://medrxiv.org/content/early/2021/08/29/2021.08.27.21262264/F2)

Figure 2: *TENM2* gene rs72831309 is associated with CKD-DKD
**A:** Regional association plot of the meta-analysis results. **B:** Forest plot of association across the contributing cohorts. **C**: SNP rs72831309 overlaps a predicted CREB1 motif sequence; data from RegulomeDB.org. **D:** Human kidney single cell RNA expression of *TENM2* showing strongest expression in podocytes (PODO), parietal epithelial cells (PEC) and proximal convoluted tubules (PCD). **E** and **F**: Tubular TENM2 expression is correlated with higher eGFR (E) and less fibrosis (F).

View this table:
[Table 1:](http://medrxiv.org/content/early/2021/08/29/2021.08.27.21262264/T1)

Table 1: GWAS meta-analysis result summary: loci with *p*<5×10−8

In addition to rs72831309, ten previously identified, mostly low-frequency or rare variants were associated with various kidney phenotypes (**Table 1, Supplemental Figure 2**). Apart from the *COL4A3* locus, none of these variants were found in the SUMMIT cohorts (filtered to MAF≥1%), and thus, these associations represent the originally reported associations from the DNCRI6. Here, we identified a secondary signal with a higher frequency (MAF=0.44) in only partial LD (D’=0.51, r2=0.08) with the lead signal at the *COL4A3* locus (rs6436688 with *p*=1.8×10−7 for Severe DKD; **Supplemental Figure 3**), which after conditional analysis for the rs55703767 lead variant remained nominally significant (*p*=0.002).

Two variants, chr3:141792314:I and rs186434345, were associated with ESRD (N=940, *p*=4.6×10−10), and with the CKD-DKD phenotype (N=2,571, *p*=4.0×10−8) in the SUMMIT-T2D and SUMMIT-T1D cohorts.The variants were absent (or poor imputation quality, r2<0.3) in all DNCRI cohorts derived from the more recent 1000Genomes phase 3 reference. When the original SUMMIT-FinnDiane GWAS was included in the analysis, both associations were non-significant (chr3:141792314:I: *p*=0.056, N=3,207, and rs186434345: *p*=0.002, N=4,782), and thus, excluded from further consideration.

### Gene-level analysis

To improve power and jointly test all available common genetic markers within a gene, SNP summary statistics from the GWAS meta-analysis were aggregated by gene and tested jointly for association using two similar programs, MAGMA and PASCAL. In addition to *COL20A1* and *SNX30* identified in the previously published DNCRI gene-level analysis6, we identified eight novel gene associations including *GPR158, LSM14A*, and *MFF* associated with severe DKD; *INIP* associated with any DKD; *PTPRN* and *RESP18* associated with CKD; and *DCLK1* and *EIF4E* associated with ESRD vs. macroalbuminuria (*p*<2.7×10−6; **Table 2, Supplemental Figures 4A-J**). Kidney eQTL data for the lead SNPs in the *INIP-SNX30* region suggested *SNX30* as the target gene (rs786959 eQTL *p*=4.6×10−7; **Supplemental Table 5**).

View this table:
[Table 2:](http://medrxiv.org/content/early/2021/08/29/2021.08.27.21262264/T2)

Table 2: Significant Gene-level DKD association results from MAGMA and PASCAL.

### Genome-wide integration of GWAS results with kidney eQTL data

We performed TWAS for each of the ten DKD meta-analyses to predict differential gene expression between cases and controls in human glomeruli and tubules, based on eQTL data in glomerular and tubulointerstitial samples from histologically normal kidneys.17 Expression levels of *AKIRIN2*, encoding a nuclear protein involved in stimulating pro-inflammatory pathways such as NF-κB44, were predicted to be higher in the tubular tissue of cases with severe DKD as compared to controls with normal AER (*p*=1.1×10−6, **Supplemental Table 6**). The eQTL data with 39 SNPs explained 5% of the variance in tubular *AKIRIN2* expression (*p*=0.01), which was further correlated with fibrosis (**Figure 3B**).

![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/29/2021.08.27.21262264/F3.medium.gif)

[Figure 3:](http://medrxiv.org/content/early/2021/08/29/2021.08.27.21262264/F3)

Figure 3: TWAS indicates increased *AKIRIN2* expression in Severe DKD.
**A:** The GWAS SNP effect sizes for association with Severe DKD (normal AER vs. macroalbuminuria or ESRD) are correlated with TWAS eQTL weights to predict *AKIRIN2* expression, suggesting that elevated *AKIRIN2* levels in tubules are associated with Severe DKD (*p*=1.1×10−6). **B:** *AKIRIN2* expression is correlated with renal fibrosis, correlation = 0.252, *p*-value=2.83×10−7.

### Gene prioritization

To identify the underlying causal genes within each of our top loci, we used the PoPS18 method that leverages genome-wide enrichment of biological annotations in combination with GWAS summary statistics to prioritize candidate genes. To increase precision, we intersected the results with both the simple nearest gene approach, as well as MAGMA gene prioritization. Four genes were both the PoPS prioritized gene and the nearest gene to the lead associated SNP, including *COL4A3, PLEKHA7, CNTN6*, and *TENM2* (**Table 1**). Of note, the *CNTN6* locus contained only two genes and the *TENM2* locus only one, with a relatively low PoP Score (−0.25). When taking the intersect between PoPS genes and genes that were within MAGMA’s top 10% of prioritized genes genome-wide, *COL4A3* was the only prioritized gene **(Supplemental Figure 5)**.

Given the robust evidence for *COL4A3* as a causal gene for DKD6, we interrogated which gene features contributed to its prioritization in MAGMA, reasoning that additional genes that are key members of these same gene features are also candidate genes for DKD and may provide a deeper understanding of disease mechanism. The specific enriched gene set that prioritized *COL4A3* for Severe DKD using the MAGMA gene prioritization approach was the fibulin 2 protein-protein interaction network (“FBLN2 PPI subnetwork”), which together with 26 correlated reconstituted gene-sets makes up the “basement membrane” meta-gene set derived in Marouli *et al*.45 These gene-sets prioritized 86 genes in MAGMA, including three other genes specific to the “FBLN2 PPI subnetwork” **(Supplemental Table 7)**.

### Kidney mQTL

To see whether the identified loci might affect the risk of DKD through DNA methylation, kidney mQTL data were queried for the top three variants at each lead locus from the GWAS meta-analyses, gene-level analyses, and TWAS. Altogether, 17 variants were significantly associated with kidney DNA methylation levels at six CpG sites (*p*<1.5×10−11**; Supplemental Table 8**). These included SNPs in the *LSM14A* gene, associated with severe DKD and cg14143166 methylation levels (*p*=1.9×10−28). Interestingly, cg14143166 methylation in blood was nominally associated with DKD status in our epigenome-wide association study (EWAS) in the UK-ROI and FinnDiane cohorts (*p*=0.03), further supporting the hypothesis that the DKD association at *LSM14A* might be mediated through methylation changes. Similarly, blood methylation levels at significant kidney mQTL CpG sites (rs7664964–cg25974308 *p*=1.1×10−11) in *EIF4E* were nominally associated also with the eGFR slope in diabetes (*p*=0.04)27.

### Gene expression and pathological phenotypes

Altogether, we identified 29 lead genes or transcripts either underlying or located near the lead SNPs, or based on gene-level analyses, TWAS, PoPS, or kidney eQTL data. Among these, the expression levels of 14 genes were significantly correlated with eGFR, glomerulosclerosis, or fibrosis in transcriptomics data obtained from 433 tubular and 335 glomerular nephrectomy samples with varying degree of diabetic and hypertensive kidney disease (*p*<2.2×10−4; **Figure 4, Supplemental Table 9**).29 The strongest correlations were observed for the tubular expression of *DCLK1* and *COL4A3 (*positive correlation), and *TENM2, COLEC11, ALLC, PLEKHA7*, and *SNX30* (negative correlation) with the level of fibrosis; and tubular *TENM2, PLEKHA7, ALLC*, and *SNX30* expression positively correlated with eGFR (*p*<1×10−7, |*r*| 0.29-0.56). In the Pima Indian kidney biopsy data, tubular *DCLK1* expression levels at the first biopsies were suggestively correlated (*p*<8.6×10−4, corrected for 29 genes and two tissues) with higher level of fibrosis, and *LSM14A* negatively correlated with the change in mesangial volume between the two study biopsies (**Figure 4)**; however, these correlations did not remain significant after a conservative correction for altogether 27 correlated tested phenotypes from two different time points. Multiple genes were nominally (*p*<0.05) correlated with these renal parameters (**Supplemental Figure 6, Supplemental Table 10**).

![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/29/2021.08.27.21262264/F4.medium.gif)

[Figure 4:](http://medrxiv.org/content/early/2021/08/29/2021.08.27.21262264/F4)

Figure 4: Tubular and glomerular gene expression of the lead genes correlates with multiple morphological and pathological renal parameters and with DKD.
Golden rectangles indicate glomerular gene expression, green ellipses tubular gene expression, and grey circles the morphological phenotypes. Blue edges indicate negative correlation, red edges positive correlation. Correlation with fibrosis, glomerulosclerosis (GlomScl), and eGFR were measured in the nephrectomy samples29, shown are correlations with p<2.2×10−4 (corrected for 29 genes, 2 tissues, and 4 tests). For the biopsy data in Pima Indians, suggestive correlations with p<8.6×10−4 are shown (corrected only for 29 genes and 2 tissues), including BX1 Fibr (fibrosis at first biopsy) and BX1 ΔMesV, change in the mesangial volume between the first and the second biopsies. Association with DKD (diabetic nephropathy, DN) was queried in two data sets, DN Wor (Woroniecka et al.40) and DN Ju (Ju et al39) with p<4.3×10−4, or p<0.05 and fold change>1.5.

### Genetic correlation and Mendelian randomization (MR) with related traits

LDSR was used to assess the shared inheritance across the genome between the DKD phenotypes and the related metabolic traits. The analysis revealed significant genetic correlation (*p*<6.4×10−4) between DKD and 15 traits including multiple obesity-related traits, mother’s age at death, T2D, coronary artery disease, HDL cholesterol, urate, and two smoking-related traits (**Figure 5**; **Supplemental Figure 7**). Among these 15 traits, all but “mother’s age at death”, which had fewer than five genome-wide significant SNP associations (*p*<5×10−8), were used in subsequent MR analysis. MR suggested that the overweight and obesity related traits were causal risk factors for DKD (**Supplemental Table 11, Figure 5B**). The causal effects were directionally consistent across methods, with no evidence of heterogeneity (I2=0-42.9%, *p*>0.05, **Supplemental Table 11**) or unbalanced horizontal pleiotropy (**Supplemental Table 12**). The MR Egger method, more robust for pleiotropic effects, further supported a causal role of higher BMI, waist, and hip circumference on DKD risk (*p*<0.05, **Supplemental Table 11, Supplemental Figure 8**). The other significant LDSR associations that indicated no causal association in MR analysis may represent reverse causality or shared upstream effects.

![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/29/2021.08.27.21262264/F5.medium.gif)

[Figure 5:](http://medrxiv.org/content/early/2021/08/29/2021.08.27.21262264/F5)

Figure 5: Genetic correlation between DKD phenotypes and various traits based on LDSR (A), and estimates of causal associations based on Mendelian Randomization (B).
For LDSR only significant trait combinations are shown (p-value < 0.05/78 = 6.4×10−4). **A: B:** Mendelian Randomization results for DKD (All vs. Ctrl) with inverse variance weighted (IVW) method for the traits significant in LDSR, and with at least 5 genome-wide significant SNPs.

## DISCUSSION

We have performed the largest GWAS meta-analysis to date on kidney complications in diabetes, including ten different phenotypic definitions in up to 26,785 individuals with either T1D or T2D, and integrated the results with emerging kidney omics data (**Figure 1**). In the single variant analysis with the combined CKD-DKD phenotype (cases with eGFR<45 ml/min/1.73 m2 and microalbuminuria or worse, vs. controls with normal AER and eGFR ≥ 60 ml/min/1.73 m2), we identified one novel locus, rs72831309 intronic in the *TENM2* gene. *TENM2* encodes the teneurin transmembrane protein 2, which is involved in cell-cell adhesion. *TENM2* rs72831309 was one of the lead loci associated with DKD in the original SUMMIT T1D GWAS meta-analysis, but it failed to reach genome-wide significance at the time4. Supporting the functionality of rs72831309, the SNP alters a predicted CREB1 transcription factor binding site. Even though no direct eQTL evidence was found for rs72831309 affecting *TENM2* expression, it was nominally associated with expression of a *TENM2* antisense transcript *CTB-178M22*.*2* in kidneys (*p*=0.007), while chromatin conformation data in the GM12878 cell line also pointed towards *TENM2* as the effector gene. Whereby kidney scRNAseq indicated *TENM2* expression in podocytes and proximal tubular cells, lower tubular *TENM2* expression was associated with renal fibrosis (*p*=2.0×10−9) and lower eGFR (*p*=1.6×10−8) in the nephrectomy samples; it was also lower among individuals with DKD versus controls (*p*=6.6×10−4, not significant after correction for multiple testing). Other variants in the *TENM2* gene were associated with multiple traits including smoking status (*p*=1.6×10−17) and BMI (*p*=2.6×10−8) in the UK Biobank. Furthermore, it was previously reported that the DNA methylation status in *TENM2* is among the strongest predictors of incident type 2 diabetes46.

Gene level analysis identified ten genes associated with various definitions of DKD. The *DCLK1* gene encodes a doublecortin-like kinase and is a known cancer stem cell marker involved in epithelial to mesenchymal cell transition.47 The histone modification based ChromHMM 15-state model for fetal kidney indicated strong transcription overlapping one of the three lead SNPs in the *DCLK1* locus (rs61948262), and ChIP-seq data supported *ZSCAN4* binding to the locus in the HEK293 kidney epithelial cell line. The lead SNPs were kidney mQTLs for *DCLK1* CpG sites (*p*=6.8×10−22), further supporting their functional relevance. Furthermore, multiple lines of evidence highlight the importance of *DCLK1* in DKD **(Figure 6, Supplemental Table 10)**: The correlation between tubular *DCLK1* expression and fibrosis was amongst the strongest correlations both in the nephrectomy samples (*p*=7.4×10−16) and in the Pima Indian biopsies (*p*=3.0×10−4), and glomerular *DCLK1* expression was nominally associated with glomerular width, mesangial volume, and podocyte foot process width in the Pima Indian biopsies (*p*<0.05). Furthermore, both glomerular and tubular *DCLK1* expression were elevated in DKD in two independent data sets (fold change [FC] 1.98, *p*=1.2×10−4 for glomeruli;39 FC 2.1, *p*=0.003 for tubules40). Finally, we have previously identified a subset of transcripts, including *DCLK1*, targeted by the early growth response-1 (egr-1) transcription factor in a murine model of DKD. In this study, *DCLK1* expression was upregulated in diabetic versus non-diabetic ApoE-/- mouse kidneys.48 Taken together, these expression data in human and experimental DKD identify *DCLK1* as a novel target.

![Figure 6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/29/2021.08.27.21262264/F6.medium.gif)

[Figure 6:](http://medrxiv.org/content/early/2021/08/29/2021.08.27.21262264/F6)

Figure 6: *DCLK1* is associated with ESRD.
**A:** The *DCLK1* gene region was associated with ESRD vs. macroalbuminuria in the MAGMA gene-level analysis (p=1.39×10−6). **EB and C:** Tubular *DCLK1* expression is highest in DKD **(B)**, and correlated with the level of fibrosis **(C)** in the nephrectomy samples. **D:** Glomerular *DCLK1* expression is higher in DKD than in healthy controls (Ju CKD39: Fold change 1.98, p=1.2×10−4). **E:** Tubular *DCLK1* expression is higher in DKD than in healthy controls (Woroniecka data set40: Fold change 2.1, p=0.003). **F and G:** Kidney *DCLK1* expression is strongest in mesangial cells in human scRNAseq data of individuals with diabetes and healthy controls37.

Of note, one of the previously identified associations with ESRD at rs116216059, intronic in the *STAC* gene, is located <200kb downstream of *DCLK3* from the same protein family. The SNP is located on a chromatin accessibility peak border in podocytes (snATACseq peak value 1.1; **Supplemental Figure 2F**). PCHiC data in GM12878 cell line supports both *STAC* and *DCLK3* as the target gene, and the expression levels of both genes were nominally associated with glomerular pathology parameters in the transcriptomics data sets (**Supplemental Table 10**).

The *SNX30* gene, encoding the Sorting Nexin Family Member 30, was associated with DKD in the gene-level analysis. While the association signal spans also the neighboring *INIP* gene, kidney eQTL data for the top SNPs in the locus pointed towards *SNX30*, with the DKD risk associated rs786959 A allele associated with higher *SNX30* expression (*p*=4.6×10−7). On the contrary, the nephrectomy sample transcriptomics data indicated that higher tubular *SNX30* was correlated with higher eGFR (*p*=5.8×10−14) and lower level of fibrosis (*p*<2×10−16); glomerular expression was correlated with less glomerulosclerosis (*p*=8×10−5). Finally, kidney *SNX30* expression was associated with higher eGFR in the general population using TWAS based on kidney eQTL17 and GWAS on eGFR49 (*p*=0.046).

The TWAS analysis based on our GWAS results, integrated with microdissected tubular and glomerular eQTL data, predicted that the *AKIRIN2* gene expression is elevated in tubules in individuals with severe DKD compared with individuals with normal AER. Furthermore, the *AKIRIN2* gene expression was highly correlated with the level of fibrosis (*p*=2.8×10−7). *AKIRIN2* encodes the akirin-2 protein, a conserved nuclear factor required for the innate immune response. Akirin-2 is a downstream effector of the toll-like receptor, tumor necrosis factor and IL-1 beta signaling pathways. It binds to nuclear NF-κB complexes and is required for the transcription of a subset of NF-κB–dependent genes such as *IL-6, CXCL10*, and *CCL5;*50 NF-κB activation drives inflammatory responses and is activated in DKD.51

The strongest regulatory evidence in RegulomeDB was obtained for rs1260634 intronic in the *LSM14A* gene: rs1260634 is located in a *ZNF362* binding sequence in HEK293 cell line, exerts strong transcription in 125 tissue types including fetal kidney chromatin state model, and affects a predicted transcription factor binding motif for KLF4, KLF12, and SP8 (**Supplemental Figure 9)**. Furthermore, in our kidney mQTL data, rs1260634 showed strong association (*p*=2.1×10−28) with a CpG site cg14143166, where methylation in blood was associated with DKD in our EWAS data (*p*=0.03). Tubular *LSM14A* expression correlated with higher eGFR (*p*=2.9×10−6), and glomerular expression with the decrease in mesangial volume (6.5×10−4). *LSM14A* encodes a Sm-like protein, thought to participate in pre-mRNA splicing, and implicated in innate antiviral responses52. Other noteworthy novel genes include *EIF4E* and *PTPRN*: *EIF4E* encodes a common mRNA translation initiation factor; its activation and/or suppression are influenced by mTOR signaling cascades involved in DKD53 as well as high glucose and high insulin environments in renal epithelial cells.54 *PTPRN* (protein tyrosine phosphatase receptor type N) encodes IA-2, a major T1D autoantigen involved in glucose-stimulated insulin secretion55. In mice, IA-2 is required to maintain normal levels of renin expression in kidneys56. Finally, the *MFF* gene encoding the mitochondrial fission factor, which was identified in our gene level analysis, has been previously related to DKD57,58. However, the association may be driven by the neighboring *COL4A3* associations, as was first suggested in the DNCRI GWAS eQTL analysis.6

In the transcriptomics analyses, as expected, all 13 significant correlations with the level of fibrosis were observed for tubular gene expression, whereby the two observed correlations for glomerulosclerosis were for glomerular expression of *SNX30* and *COLEC11* (**Figure 3**). Interestingly, eight out of ten correlations with eGFR were obtained for tubular rather than glomerular gene expression, supporting the importance of tubular damage for the loss of renal function.

Previous genetic risk scores and LDSR analyses have supported a role of T2D genetic factors, obesity and smoking for the development of DKD4. While there is a strong epidemiological link between DKD and coronary artery disease in diabetes1, this is the first study to report also a genetic correlation between these major diabetic micro- and macrovascular complications. Among the lipid traits, significant correlation with DKD was found only for lower HDL cholesterol, despite previous MR studies within general population cohorts implicating HDL as a marker of dyslipidemia rather than a causal factor59. Indeed, our subsequent MR found no evidence of causality between HDL cholesterol and DKD; in concordance with our previous MR on BMI,60 only body weight related measurements were causal risk factors for DKD. Of note, our current MR was in line with our previous MR in T1D suggesting that serum urate levels are not a causal risk factor for DKD61; similar negative results were also reported for non-diabetic CKD62.

The majority of kidney disease in individuals with T1D is considered to occur due to diabetic nephropathy, histologically characterized by thickening of glomerular basement membrane and mesangial expansion, as well as renal tubular, interstitial and arteriolar lesions.63 In individuals with T2D, only a proportion of DKD is purely due to diabetic nephropathy, whereas aging, obesity, and hypertension also contribute to the development of kidney complications. Thus, including individuals with T2D in the meta-analysis increases the heterogeneity of the underlying disease, and may select findings for those directly related to hyperglycemia as the common endpoint shared by both forms of diabetes. However, as T2D represents 95% of all diabetes cases, including those individuals greatly increases statistical power for our current work and future GWAS meta-analyses integrating multiple subtypes of diabetes to identify shared genetic risk factors for DKD.

## Supporting information

Supplemental material [[supplements/262264_file02.pdf]](pending:yes)

## Data Availability

The GWAS meta-analysis results can be accessed via the T1D, T2D, and Cardiometabolic Disease (CMD) Knowledge Portals ([https://t1d.hugeamp.org/](https://t1d.hugeamp.org/); [https://t2d.hugeamp.org/](https://t2d.hugeamp.org/); [https://hugeamp.org/](https://hugeamp.org/)), and downloaded on their respective downloads pages.

## Author contributions

NS and JC performed the main analyses and drafted the manuscript. VN, XS, HL, EA, NvZ, EHD, DF, LS performed statistical analyses, contributed to interpretation of the results, and revised the manuscript. RS, CF, EV, VH, EB, GM, DA, RD, HCL, RGN, AJM, CP, CG, APM, LG, MIM, MK, KS, JNH, JCF, PHG contributed to the data acquisition and interpretation of the results, and critically revised the manuscript.

## Disclosures

P-H.G. reports receiving lecture honorariums from Astellas, Astra Zeneca, Boehringer Ingelheim, Eli Lilly, Medscape, MSD, Mundipharma, Novo Nordisk, PeerVoice, Sanofi, SCIARC and being an advisory board member of Astellas, Astra Zeneca, Bayer, Boehringer Ingelheim, Eli Lilly, Elo Water, Medscape, MSD, Mundipharma, Nestlé, Novo Nordisk, PeerVoice, Sanofi and Sciarc. J.C.F. has received speaking honoraria from AstraZeneca, Merck and Novo Nordisk for scientific talks over which he had full control of content.

## Funding

This study was funded by JDRF (grants S-SRA-2014-276-Q-R and 17-2013-7), and NIDDK R01 DK105154. Research in FinnDiane was supported by Folkhälsan Research Foundation, Wilhelm and Else Stockmann Foundation, “Liv och Hälsa” Society, Helsinki University Central Hospital Research Funds (EVO TYH2018207), Academy of Finland (299200 and 316664), and Novo Nordisk Foundation (NNF OC0013659). Acquisition of clinical data and samples from the Pima Indians was supported by the Intramural Research Program of the National Institute of Diabetes and Digestive and Kidney Diseases and by the American Diabetes Association (Clinical Science Award 1-08-CR-42). JBC was funded by the National Institute of Diabetes and Digestive and Kidney Diseases (K99DK127196). EA was funded by grants from the Swedish Research Council (2017-02688, 2020-02191) and the Novo Nordisk foundation (NNF18OC0034408).

## Supplemental Table of contents

Supplemental Table 1: A total of ten case – control definitions for diabetic kidney disease (DKD; used as a general term to describe renal complications in diabetes)

Supplemental Table 2: Number of individuals in each cohort for each phenotypic comparison

Supplemental Table 3: Key characteristics of the genotyping and statistical analyses in DNCRI and SUMMIT cohorts.

Supplemental Table 4: Association details for the novel genome-wide significant locus rs72831309 in *TENM2*.

Supplemental Table 5: Kidney eQTL associations with p<1×10−4 in tubular or glomerular eQTL data, or in the kidney eQTL meta-analysis for the lead SNPs.

Supplemental Table 6: TWAS results with p<1×10−4.

Supplemental Table 7: Highly correlated reconstituted gene-sets that make-up the “basement membrane” meta-gene set derived in Marouli et al.

Supplemental Table 8: Significant kidney mQTL associations (p<1.46×10−11 Bonferroni-adjusted genome-wide significance; 1×10−7 suggestive significance) for lead loci.

Supplemental Table 9: Correlation between glomerular and tubular gene expression and glomerulosclerosis, fibrosis, and eGFR in nephrectomy samples.

Supplemental Table 10: Gene centric summary of the lead genes

Supplemental Table 11: Mendelian Randomization (MR) results for DKD (All vs Ctrl)

Supplemental Table 12. Egger intercepts for Mendelian Randomization analyses on DKD (all vs. Ctrl phenotype).

Supplemental Figure 1: Manhattan and QQ-plots of the ten DKD GWAS meta-analysis results.

Supplemental Figure 2: Regional association plots for the GWAS meta-analysis lead loci (A-K).

Supplemental Figure 3: Regional association plot for the *COL4A3* gene region associated with Severe DKD, indicating a secondary association peak at chr2:228259302 (rs6436688, effect allele (A) frequency 56%, OR = 1.13 (95% confidence interval 1.08 – 1.19), p-value 1.79×10−7).

Supplemental Figure 4: Regional association plots for the gene-level analysis results from MAGMA and PASCAL analysis.

Supplemental Figure 5: Gene prioritization for the *COL4A3* gene at lead SNP rs5570367 associated with Severe DKD using multiple intersecting gene prioritization approaches (PoPS, nearest gene, and MAGMA).

Supplemental Figure 6: Tubular and glomerular gene expression of the lead genes correlate with multiple morphological and pathological renal parameters. Golden

Supplemental Figure 7: Genetic correlation between DKD and related traits based on LD score regression.

Supplemental Figure 8: Mendelian Randomization scatter plots for SNP effects for the metabolic traits vs. DKD (All vs. Ctrl).

Supplemental Figure 9: rs1260634 intronic in the *ALLC* gene affects the predicted binding motifs for KLF12, KLF4, and SP8 (top to bottom).

*   Received August 27, 2021.
*   Revision received August 27, 2021.
*   Accepted August 29, 2021.


*   © 2021, Posted by Cold Spring Harbor Laboratory

The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission.

## References

1.  1.Harjutsalo, V., Thomas, M. C., Forsblom, C., Groop, P. H., & FinnDiane Study Group. Risk of coronary artery disease and stroke according to sex and presence of diabetic nephropathy in type 1 diabetes. Diabetes Obes.Metab. 20, 2759–2767 (2018).
    
    

2.  2.Groop, P. H. et al. The presence and severity of chronic kidney disease predicts all-cause mortality in type 1 diabetes. Diabetes 58, 1651–1658 (2009).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZGlhYmV0ZXMiO3M6NToicmVzaWQiO3M6OToiNTgvNy8xNjUxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDgvMjkvMjAyMS4wOC4yNy4yMTI2MjI2NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

3.  3.World Health Organization. The top 10 causes of death. (2020).
    
    

4.  4.Sandholm, N. et al. The Genetic Landscape of Renal Complications in Type 1 Diabetes. J.Am.Soc.Nephrol. 28, 557–574 (2017).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiam5lcGhyb2wiO3M6NToicmVzaWQiO3M6ODoiMjgvMi81NTciO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wOC8yOS8yMDIxLjA4LjI3LjIxMjYyMjY0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

5.  5.Harjutsalo, V., Katoh, S., Sarti, C., Tajima, N. & Tuomilehto, J. Population-based assessment of familial clustering of diabetic nephropathy in type 1 diabetes. Diabetes 53, 2449–2454 (2004).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZGlhYmV0ZXMiO3M6NToicmVzaWQiO3M6OToiNTMvOS8yNDQ5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDgvMjkvMjAyMS4wOC4yNy4yMTI2MjI2NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

6.  6.Salem, R. M. et al. Genome-Wide Association Study of Diabetic Kidney Disease Highlights Biology Involved in Glomerular Basement Membrane Collagen. J.Am.Soc.Nephrol. 30, 2000–2016 (2019).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiam5lcGhyb2wiO3M6NToicmVzaWQiO3M6MTA6IjMwLzEwLzIwMDAiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wOC8yOS8yMDIxLjA4LjI3LjIxMjYyMjY0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

7.  7.Sandholm, N. et al. New susceptibility Loci associated with kidney disease in type 1 diabetes. PLoS Genet. 8, e1002921 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1002921&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23028342&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F29%2F2021.08.27.21262264.atom) 

8.  8.Sandholm, N. et al. Chromosome 2q31.1 associates with ESRD in women with type 1 diabetes. J.Am.Soc.Nephrol. 24, 1537–1543 (2013).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiam5lcGhyb2wiO3M6NToicmVzaWQiO3M6MTA6IjI0LzEwLzE1MzciO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wOC8yOS8yMDIxLjA4LjI3LjIxMjYyMjY0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

9.  9.Sandholm, N. et al. Genome-wide association study of urinary albumin excretion rate in patients with type 1 diabetes. Diabetologia 57, 1143–1153 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00125-014-3202-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24595857&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F29%2F2021.08.27.21262264.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000335922800007&link_type=ISI) 

10. 10.Sandholm, N. et al. Confirmation of GLRA3 as a susceptibility locus for albuminuria in Finnish patients with type 1 diabetes. Sci.Rep. 8, 12408–018-29211–1 (2018).
    
    

11. 11.van Zuydam, N. R. et al. A Genome-Wide Association Study of Diabetic Kidney Disease in Subjects With Type 2 Diabetes. Diabetes 67, 1414–1427 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2337/db18-1414-P&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29703844&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F29%2F2021.08.27.21262264.atom) 

12. 12.Pruim, R. J. et al. LocusZoom: regional visualization of genome-wide association scan results. Bioinformatics 26, 2336–2337 (2010).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btq419&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20634204&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F29%2F2021.08.27.21262264.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000281714100054&link_type=ISI) 

13. 13.Yang, J. et al. Conditional and joint multiple-SNP analysis of GWAS summary statistics identifies additional variants influencing complex traits. Nat Genet 44, 369–375, S1-3 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.2213&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22426310&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F29%2F2021.08.27.21262264.atom) 

14. 14.de Leeuw, C. A., Mooij, J. M., Heskes, T. & Posthuma, D. MAGMA: generalized gene-set analysis of GWAS data. PLoS Comput.Biol. 11, e1004219 (2015).
    
    

15. 15.Lamparter, D., Marbach, D., Rueedi, R., Kutalik, Z. & Bergmann, S. Fast and Rigorous Computation of Gene and Pathway Scores from SNP-Based Summary Statistics. PLoS Comput.Biol. 12, e1004714 (2016).
    
    

16. 16.Barbeira, A. N. et al. Exploring the phenotypic consequences of tissue specific gene expression variation inferred from GWAS summary statistics. Nat.Commun. 9, 1825–018-03621–1 (2018).
    
    

17. 17.Qiu, C. et al. Renal compartment-specific genetic variation analyses identify new pathways in chronic kidney disease. Nat.Med. (2018) doi:10.1038/s41591-018-0194-4.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-018-0194-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30275566&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F29%2F2021.08.27.21262264.atom) 

18. 18.Weeks, E. M. et al. Leveraging polygenic enrichments of gene features to predict genes underlying complex traits and diseases. medRxiv 2020.09.08.20190561 (2020) doi:10.1101/2020.09.08.20190561.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wOS4wOC4yMDE5MDU2MXYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDgvMjkvMjAyMS4wOC4yNy4yMTI2MjI2NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

19. 19.Fine, R. S., Pers, T. H., Amariuta, T., Raychaudhuri, S. & Hirschhorn, J. N. Benchmarker: An Unbiased, Association-Data-Driven Strategy to Evaluate Gene Prioritization Algorithms. Am J Hum Genet 104, 1025–1039 (2019).
    
    

20. 20.Võsa, U. et al. Unraveling the polygenic architecture of complex traits using blood eQTL metaanalysis. bioRxiv 447367 (2018) doi:10.1101/447367.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czo4OiI0NDczNjd2MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzA4LzI5LzIwMjEuMDguMjcuMjEyNjIyNjQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

21. 21.Ko, Y. A. et al. Genetic-Variation-Driven Gene-Expression Changes Highlight Genes with Important Functions for Kidney Disease. Am.J.Hum.Genet. 100, 940–953 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2017.05.004&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28575649&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F29%2F2021.08.27.21262264.atom) 

22. 22.GTEx Consortium. The GTEx Consortium atlas of genetic regulatory effects across human tissues. Science 369, 1318–1330 (2020).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzNjkvNjUwOS8xMzE4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDgvMjkvMjAyMS4wOC4yNy4yMTI2MjI2NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

23. 23.Gillies, C. E. et al. An eQTL Landscape of Kidney Tissue in Human Nephrotic Syndrome. Am J Hum Genet 103, 232–244 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2018.07.004&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30057032&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F29%2F2021.08.27.21262264.atom) 

24. 24.Willer, C. J., Li, Y. & Abecasis, G. R. METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics 26, 2190–2191 (2010).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btq340&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20616382&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F29%2F2021.08.27.21262264.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000281738900017&link_type=ISI) 

25. 25.Zhou, W., Triche, T. J., Laird, P. W. & Shen, H. SeSAMe: reducing artifactual detection of DNA methylation by Infinium BeadChips in genomic deletions. Nucleic Acids Res 46, e123 (2018).
    
    

26. 26.Shabalin, A. A. Matrix eQTL: ultra fast eQTL analysis via large matrix operations. Bioinformatics 28, 1353–1358 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/bts163&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22492648&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F29%2F2021.08.27.21262264.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000304053300009&link_type=ISI) 

27. 27.Sheng, X. et al. Systematic integrated analysis of genetic and epigenetic variation in diabetic kidney disease. Proc Natl Acad Sci U S A 117, 29013–29024 (2020).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE3LzQ2LzI5MDEzIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDgvMjkvMjAyMS4wOC4yNy4yMTI2MjI2NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

28. 28.Smyth, L. J. et al. Assessment of differentially methylated loci in individuals with end-stage kidney disease attributed to diabetic kidney disease: an exploratory study. Clin.Epigenetics 13, 99–021-01081–x (2021).
    
    

29. 29.Guan, Y. et al. A single genetic locus controls both expression of DPEP1/CHMP1A and kidney disease development via ferroptosis. Nat Commun 12, 5078 (2021).
    
    

30. 30.Nair, V. et al. A molecular morphometric approach to diabetic kidney disease can link structure to function and outcome. Kidney Int 93, 439–449 (2018).
    
    

31. 31.Looker, H. C. et al. Changes in Albuminuria But Not GFR are Associated with Early Changes in Kidney Structure in Type 2 Diabetes. J Am Soc Nephrol 30, 1049–1059 (2019).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiam5lcGhyb2wiO3M6NToicmVzaWQiO3M6OToiMzAvNi8xMDQ5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDgvMjkvMjAyMS4wOC4yNy4yMTI2MjI2NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

32. 32.Mifsud, B. et al. Mapping long-range promoter contacts in human cells with high-resolution capture Hi-C. Nat.Genet. 47, 598–606 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3286&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25938943&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F29%2F2021.08.27.21262264.atom) 

33. 33.Choy, M. K. et al. Promoter interactome of human embryonic stem cell-derived cardiomyocytes connects GWAS regions to cardiac gene networks. Nat.Commun. 9, 2526-018-04931–0 (2018).
    
    

34. 34.Javierre, B. M. et al. Lineage-Specific Genome Architecture Links Enhancers and Non-coding Disease Variants to Target Gene Promoters. Cell 167, 1369–1384.e19 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2016.09.037&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29428417&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F29%2F2021.08.27.21262264.atom) 

35. 35.Miguel-Escalada, I. et al. Human pancreatic islet 3D chromatin architecture provides insights into the genetics of type 2 diabetes. bioRxiv 400291 (2018).
    
    

36. 36.Sheng, X. et al. Mapping the genetic architecture of human traits to cell types in the kidney identifies mechanisms of disease and potential treatments. bioRxiv 2020.11.09.375592 (2020) doi:10.1101/2020.11.09.375592.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czoxOToiMjAyMC4xMS4wOS4zNzU1OTJ2MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzA4LzI5LzIwMjEuMDguMjcuMjEyNjIyNjQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

37. 37.Wilson, P. C. et al. The single-cell transcriptomic landscape of early human diabetic nephropathy. Proc Natl Acad Sci U S A 116, 19619–19625 (2019).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE2LzM5LzE5NjE5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDgvMjkvMjAyMS4wOC4yNy4yMTI2MjI2NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

38. 38.Boyle, A. P. et al. Annotation of functional variation in personal genomes using RegulomeDB. Genome Res. 22, 1790–1797 (2012).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ2Vub21lIjtzOjU6InJlc2lkIjtzOjk6IjIyLzkvMTc5MCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzA4LzI5LzIwMjEuMDguMjcuMjEyNjIyNjQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

39. 39.Ju, W. et al. Defining cell-type specificity at the transcriptional level in human disease. Genome Res. 23, 1862–1873 (2013).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ2Vub21lIjtzOjU6InJlc2lkIjtzOjEwOiIyMy8xMS8xODYyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDgvMjkvMjAyMS4wOC4yNy4yMTI2MjI2NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

40. 40.Woroniecka, K. I. et al. Transcriptome analysis of human diabetic kidney disease. Diabetes 60, 2354–2369 (2011).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZGlhYmV0ZXMiO3M6NToicmVzaWQiO3M6OToiNjAvOS8yMzU0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDgvMjkvMjAyMS4wOC4yNy4yMTI2MjI2NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

41. 41.Bulik-Sullivan, B. K. et al. LD Score regression distinguishes confounding from polygenicity in genome-wide association studies. Nat.Genet. 47, 291–295 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3211&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25642630&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F29%2F2021.08.27.21262264.atom) 

42. 42.Hemani, G. et al. The MR-Base platform supports systematic causal inference across the human phenome. Elife 7, e34408 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7554/eLife.34408&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29846171&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F29%2F2021.08.27.21262264.atom) 

43. 43.Bowden, J., Davey Smith, G., Haycock, P. C. & Burgess, S. Consistent Estimation in Mendelian Randomization with Some Invalid Instruments Using a Weighted Median Estimator. Genet Epidemiol 40, 304–314 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/gepi.21965&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27061298&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F29%2F2021.08.27.21262264.atom) 

44. 44.Tartey, S. et al. Akirin2 is critical for inducing inflammatory genes by bridging IκB-ζ and the SWI/SNF complex. EMBO J 33, 2332–2348 (2014).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiZW1ib2pubCI7czo1OiJyZXNpZCI7czoxMDoiMzMvMjAvMjMzMiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzA4LzI5LzIwMjEuMDguMjcuMjEyNjIyNjQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

45. 45.Marouli, E. et al. Rare and low-frequency coding variants alter human adult height. Nature 542, 186–190 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature21039&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28146470&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F29%2F2021.08.27.21262264.atom) 

46. 46.Ouni, M. et al. Epigenetic Changes in Islets of Langerhans Preceding the Onset of Diabetes. Diabetes 69, 2503–2517 (2020).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZGlhYmV0ZXMiO3M6NToicmVzaWQiO3M6MTA6IjY5LzExLzI1MDMiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wOC8yOS8yMDIxLjA4LjI3LjIxMjYyMjY0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

47. 47.Chandrakesan, P. et al. Regulatory Roles of Dclk1 in Epithelial Mesenchymal Transition and Cancer Stem Cells. J.Carcinog.Mutagen. 7, doi:10.4172/2157-2518.1000257. Epub 2016 Mar 7 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.4172/2157-2518.1000257&link_type=DOI) 

48. 48.Brennan, E. P. et al. Lipoxins Regulate the Early Growth Response-1 Network and Reverse Diabetic Kidney Disease. J.Am.Soc.Nephrol. 29, 1437–1448 (2018).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiam5lcGhyb2wiO3M6NToicmVzaWQiO3M6OToiMjkvNS8xNDM3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDgvMjkvMjAyMS4wOC4yNy4yMTI2MjI2NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

49. 49.Wuttke, M. et al. A catalog of genetic loci associated with kidney function from analyses of a million individuals. Nat Genet 51, 957–972 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-019-0407-x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31152163&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F29%2F2021.08.27.21262264.atom) 

50. 50.Goto, A. et al. Akirins are highly conserved nuclear proteins required for NF-kappaB-dependent gene expression in drosophila and mice. Nat.Immunol. 9, 97–104 (2008).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ni1543&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18066067&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F29%2F2021.08.27.21262264.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000251820800016&link_type=ISI) 

51. 51.Mezzano, S. et al. NF-κB activation and overexpression of regulated genes in human diabetic nephropathy. Nephrology Dialysis Transplantation 19, 2505–2512 (2004).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ndt/gfh207&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15280531&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F29%2F2021.08.27.21262264.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000224354300012&link_type=ISI) 

52. 52.Li, Y. et al. LSm14A is a processing body-associated sensor of viral nucleic acids that initiates cellular antiviral response in the early phase of viral infection. Proc Natl Acad Sci U S A 109, 11770–11775 (2012).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTA5LzI5LzExNzcwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDgvMjkvMjAyMS4wOC4yNy4yMTI2MjI2NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

53. 53.Lieberthal, W. & Levine, J. S. The role of the mammalian target of rapamycin (mTOR) in renal disease. Journal of the American Society of Nephrology 20, 2493–2502 (2009).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiam5lcGhyb2wiO3M6NToicmVzaWQiO3M6MTA6IjIwLzEyLzI0OTMiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wOC8yOS8yMDIxLjA4LjI3LjIxMjYyMjY0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

54. 54.Mariappan, M. M., Feliers, D., Mummidi, S., Choudhury, G. G. & Kasinath, B. S. High glucose, high insulin, and their combination rapidly induce laminin-beta1 synthesis by regulation of mRNA translation in renal epithelial cells. Diabetes 56, 476–485 (2007).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZGlhYmV0ZXMiO3M6NToicmVzaWQiO3M6ODoiNTYvMi80NzYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wOC8yOS8yMDIxLjA4LjI3LjIxMjYyMjY0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

55. 55.Saeki, K. et al. Targeted Disruption of the Protein Tyrosine Phosphatase-Like Molecule IA-2 Results in Alterations in Glucose Tolerance Tests and Insulin Secretion. Diabetes 51, 1842–1850 (2002).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZGlhYmV0ZXMiO3M6NToicmVzaWQiO3M6OToiNTEvNi8xODQyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDgvMjkvMjAyMS4wOC4yNy4yMTI2MjI2NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

56. 56.Kim, S. M. et al. Dense-core vesicle proteins IA-2 and IA-2{beta} affect renin synthesis and secretion through the {beta}-adrenergic pathway. Am J Physiol Renal Physiol 296, F382–389 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1152/ajprenal.90543.2008&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19019914&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F29%2F2021.08.27.21262264.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000262782700019&link_type=ISI) 

57. 57.Gao, P. et al. DsbA-L deficiency exacerbates mitochondrial dysfunction of tubular cells in diabetic kidney disease. Clin.Sci. 134, 677–694 (2020).
    
    

58. 58.Sheng, J. et al. NR4A1 Promotes Diabetic Nephropathy by Activating Mff-Mediated Mitochondrial Fission and Suppressing Parkin-Mediated Mitophagy. Cell.Physiol.Biochem. 48, 1675–1693 (2018).
    
    

59. 59.Voight, B. F. et al. Plasma HDL cholesterol and risk of myocardial infarction: a mendelian randomisation study. Lancet 380, 572–580 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(12)60312-2&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22607825&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F29%2F2021.08.27.21262264.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000307511400030&link_type=ISI) 

60. 60.Todd, J. N. et al. Genetic Evidence for a Causal Role of Obesity in Diabetic Kidney Disease. Diabetes 64, 4238–4246 (2015).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZGlhYmV0ZXMiO3M6NToicmVzaWQiO3M6MTA6IjY0LzEyLzQyMzgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wOC8yOS8yMDIxLjA4LjI3LjIxMjYyMjY0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

61. 61.Ahola, A. J. et al. The serum uric acid concentration is not causally linked to diabetic nephropathy in type 1 diabetes. Kidney Int. 91, 1178–1185 (2017).
    
    

62. 62.Jordan, D. M. et al. No causal effects of serum urate levels on the risk of chronic kidney disease: A Mendelian randomization study. PLoS medicine 16, e1002725 (2019).
    
    

63. 63.Fioretto, P. & Mauer, M. Histopathology of diabetic nephropathy. Semin.Nephrol. 27, 195–207 (2007).