The Genetic Landscape of Renal Complications in Type 1 Diabetes
Total Page:16
File Type:pdf, Size:1020Kb
Supplementary information: Genome-wide dissection of diabetic kidney disease The Genetic Landscape of Renal Complications in Type 1 Diabetes SUPPLEMENTAL INFORMATION Niina Sandholm, Natalie Van Zuydam, Emma Ahlqvist, Thorhildur Juliusdottir, Harshal A. Deshmukh, N. William Rayner, Barbara Di Camillo, Carol Forsblom, Joao Fadista, Daniel Ziemek, Rany M. Salem, Linda T. Hiraki, Marcus Pezzolesi, David Trégouët, Emma Dahlström, Erkka Valo, Nikolay Oskolkov, Claes Ladenvall, M. Loredana Marcovecchio, Jason Cooper, Francesco Sambo, Alberto Malovini, Marco Manfrini, Amy Jayne McKnight, Maria Lajer, Valma Harjutsalo, Daniel Gordin, Maija Parkkonen, Jaakko Tuomilehto, Valeriya Lyssenko, Paul M. McKeigue, Stephen S. Rich, Mary Julia Brosnan, Eric Fauman, Riccardo Bellazzi, Peter Rossing, Samy Hadjadj, Andrzej Krolewski, Andrew D. Paterson, Jose C. Florez, Joel N. Hirschhorn, Alexander P. Maxwell, David Dunger, The DCCT/EDIC Study Group, GENIE Consortium, The FinnDiane Study Group, Claudio Cobelli, Helen M. Colhoun, Leif Groop, Mark I. McCarthy, Per- Henrik Groop, on behalf of The SUMMIT Consortium 1 Supplementary information: Genome-wide dissection of diabetic kidney disease COMPLETE METHODS ........................................................................................................................................... 4 Subjects with Type 1 diabetes (T1D) .............................................................................................................................. 4 Genotypes ....................................................................................................................................................................... 7 Statistical methods ......................................................................................................................................................... 9 SUPPLEMENTAL TABLES ...................................................................................................................................... 14 Supplemental Table 1: Proportion of phenotypic variance explained by the GWAS genotypes in FinnDiane, estimated with GCTA .................................................................................................................................................... 14 Supplemental Table 2: Information on genotyping, and clinical characteristics of the discovery and replication patients divided based on the seven case – control definitions .................................................................................. 15 Supplemental Table 3: Number of subjects included in the analysis in the discovery and in silico replication cohorts ...................................................................................................................................................................................... 16 Supplemental Table 4: Statistical power to detect association with ‘combined DKD’ with genome-wide significance (p<5×10-8) at the discovery stage. ............................................................................................................................... 17 Supplemental Table 5: Association analysis results for the 101 GWAS SNPs selected for in silico replication. .......... 18 Supplemental Table 6: Statistical power to detect association with ‘Combined DKD’ with genome-wide significance (p<5x10-8) with the two-stage study design. ................................................................................................................ 19 Supplemental Table 7: Association at the AFF3 locus with ‘ESRD vs. non-ESRD’ phenotypic comparison, conditional on the previously reported lead SNP rs7583877 .......................................................................................................... 20 Supplemental Table 8: Statistical power to detect association with the ‘Late DKD’ phenotype for varying odds ratio and risk allele frequency. .............................................................................................................................................. 21 Supplemental Table 9: Evaluation of previously reported candidate genes or GWAS loci on kidney complications in type 1 and type 2 diabetes, or GWAS on CKD in the general population. ................................................................... 22 Supplemental Table 10: Association between diabetic kidney complications and genetic risk scores of related phenotypes ................................................................................................................................................................... 24 Supplemental Table 11: MAGENTA Gene set enrichment results with FDR<0.05 ....................................................... 28 Supplemental Table 12: Characteristics of the patients selected for the whole exome sequencing .......................... 30 Supplemental Table 13. Top 20 results of single variant analysis for WES ‘ESRD vs. no DKD’ using the score test .... 31 Supplemental Table 14. Top 20 associations for WES ‘Late DKD’ using the single variant score test ......................... 32 Supplemental Table 15: Top 10 associations to WES ‘Late DKD’ with VT with 4 masks .............................................. 33 Supplemental Table 16: Top 10 associations to WES ‘Late DKD’ using SKAT-O with 4 different masks. ..................... 34 Supplemental Table 17: Top 10 associations to WES ‘Late DKD’ using SKAT with 4 masks ......................................... 35 Supplemental table 18: Top 10 associations to WES ‘ESRD vs. no DKD’ with VT and 4 masks .................................... 36 Supplemental Table 19: Top 10 associations to WES ‘ESRD vs. no DKD’ with SKAT-O with 4 masks .......................... 37 Supplemental Table 20: Top 10 associations to WES ‘ESRD vs. no DKD’ with SKAT with 4 masks .............................. 38 Supplemental Table 21: Gene-sets which showed enrichment in WES association data on the ‘Late DKD’ phenotype, with permutation (N=100) results ................................................................................................................................ 39 Supplemental Table 22: ABACUS association analysis results for the top SNPs from the GWAS discovery ................ 41 2 Supplementary information: Genome-wide dissection of diabetic kidney disease Supplemental Table 23: ABACUS association analysis results for the top SNPs from the WES discovery ................... 43 Supplemental Table 24: DAVID functional clustering of the genes mapped by the SNPs selected by the ABACUS software on GWAS data ............................................................................................................................................... 44 Supplemental Table 25: DAVID functional clustering of the genes mapped by the SNPs selected by the ABACUS software on WES data .................................................................................................................................................. 46 Supplemental Table 26: Phenotype definitions. Table A: albuminuria- and eGFR based definitions. Table B: Case – control phenotypes. ..................................................................................................................................................... 48 Supplemental Table 27: Membership of the GENIE Consortium ................................................................................. 49 Supplemental Table 28: List of the FinnDiane centers and participating physicians and nurses. ............................... 50 Supplemental Table 29: Membership of the SUMMIT Consortium ............................................................................. 53 SUPPLEMENTAL FIGURES .................................................................................................................................... 61 Supplemental Figure 1: Schematic picture of the study plan. In the GWAS setting, the stage 1 included T1D patients from the FinnDiane, EURODIAB, SDR, and Cambridge studies. Stage 1 GWAS meta-analysis results were used for evaluation of the previously reported loci, analysis of genetic risk scores, LD score regression, and for the pathway analyses. Stage 2 included patients from the UK-ROI, GoKinD US, French-Danish effort, DCCT/EDIC, and Joslin studies. Stage 3 replication consisted of additional T1D FinnDiane patients not part of the FinnDiane GWAS. Whole exome sequencing (WES) included patients from the FinnDiane, SDR, and Steno studies. Finally, the bivariate association analyses were performed in all GWAS stage 1 studies and in WES studies. ............................................. 61 Supplemental Figure 2: Manhattan and QQ-plots for the seven studied phenotype definitions. .............................. 62 Supplemental Figure 3: LocusZoom and Forest plots of the top loci. .......................................................................... 65 Supplemental Figure 4: P-value distribution of association at the previously reported loci for DKD or CKD in the general population. ...................................................................................................................................................... 66 Supplemental Figure 5: Association at previously reported loci plotted by the previously reported A) p-values and by B) Z-scores. ............................................................................................................................................................... 67 Supplemental Figure 6: Genome-wide comparison of the