Genetic drivers of heterogeneity in type 2 diabetes pathophysiology.

Type 2 diabetes (T2D) is a heterogeneous disease that develops through diverse pathophysiological processes1,2 and molecular mechanisms that are often specific to cell type3,4. Here, to characterize the genetic contribution to these processes across ancestry groups, we aggregate genome-wide association study data from 2,535,601 individuals (39.7% not of European ancestry), including 428,452 cases of T2D. We identify 1,289 independent association signals at genome-wide significance (P < 5 × 10-8) that map to 611 loci, of which 145 loci are, to our knowledge, previously unreported. We define eight non-overlapping clusters of T2D signals that are characterized by distinct profiles of cardiometabolic trait associations. These clusters are differentially enriched for cell-type-specific regions of open chromatin, including pancreatic islets, adipocytes, endothelial cells and enteroendocrine cells. We build cluster-specific partitioned polygenic scores5 in a further 279,552 individuals of diverse ancestry, including 30,288 cases of T2D, and test their association with T2D-related vascular outcomes. Cluster-specific partitioned polygenic scores are associated with coronary artery disease, peripheral artery disease and end-stage diabetic nephropathy across ancestry groups, highlighting the importance of obesity-related processes in the development of vascular outcomes. Our findings show the value of integrating multi-ancestry genome-wide association study data with single-cell epigenomics to disentangle the aetiological heterogeneity that drives the development and progression of T2D. This might offer a route to optimize global access to genetically informed diabetes care.

metaGWASmanager: a toolbox for an automated workflow from phenotypes to meta-analysis in GWAS consortia.

SUMMARY: This article introduces the metaGWASmanager, which streamlines genome-wide association studies within large-scale meta-analysis consortia. It is a toolbox for both the central consortium analysis group and participating studies to generate homogeneous phenotypes, minimize unwanted variability from inconsistent methodologies, ensure high-quality association results, and implement time-efficient quality control workflows. The toolbox features a plug-in-based approach for customization of association testing. RESULTS: The metaGWASmanager toolbox has been successfully deployed in both the CKDGen and MetalGWAS Initiative consortia across hundreds of participating studies, demonstrating its effectiveness in GWAS analysis optimization by automating routine tasks and ensuring the value and reliability of association results, thus, ultimately promoting scientific discovery. We provide a simulated data set with examples for script customization so that readers can reproduce the pipeline at their convenience. AVAILABILITY AND IMPLEMENTATION: GitHub: https://github.com/genepi-freiburg/metaGWASmanager.

Implication of transcription factor FOXD2 dysfunction in syndromic congenital anomalies of the kidney and urinary tract (CAKUT).

Congenital anomalies of the kidney and urinary tract (CAKUT) are the predominant cause for chronic kidney disease below age 30 years. Many monogenic forms have been discovered due to comprehensive genetic testing like exome sequencing. However, disease-causing variants in known disease-associated genes only explain a proportion of cases. Here, we aim to unravel underlying molecular mechanisms of syndromic CAKUT in three unrelated multiplex families with presumed autosomal recessive inheritance. Exome sequencing in the index individuals revealed three different rare homozygous variants in FOXD2, encoding a transcription factor not previously implicated in CAKUT in humans: a frameshift in the Arabic and a missense variant each in the Turkish and the Israeli family with segregation patterns consistent with autosomal recessive inheritance. CRISPR/Cas9-derived Foxd2 knockout mice presented with a bilateral dilated kidney pelvis accompanied by atrophy of the kidney papilla and mandibular, ophthalmologic, and behavioral anomalies, recapitulating the human phenotype. In a complementary approach to study pathomechanisms of FOXD2-dysfunction-mediated developmental kidney defects, we generated CRISPR/Cas9-mediated knockout of Foxd2 in ureteric bud-induced mouse metanephric mesenchyme cells. Transcriptomic analyses revealed enrichment of numerous differentially expressed genes important for kidney/urogenital development, including Pax2 and Wnt4 as well as gene expression changes indicating a shift toward a stromal cell identity. Histology of Foxd2 knockout mouse kidneys confirmed increased fibrosis. Further, genome-wide association studies suggest that FOXD2 could play a role for maintenance of podocyte integrity during adulthood. Thus, our studies help in genetic diagnostics of monogenic CAKUT and in understanding of monogenic and multifactorial kidney diseases.

Meta-GWAS Reveals Novel Genetic Variants Associated with Urinary Excretion of Uromodulin.

BACKGROUND: Uromodulin, the most abundant protein excreted in normal urine, plays major roles in kidney physiology and disease. The mechanisms regulating the urinary excretion of uromodulin remain essentially unknown. METHODS: We conducted a meta-analysis of genome-wide association studies for raw (uUMOD) and indexed to creatinine (uUCR) urinary levels of uromodulin in 29,315 individuals of European ancestry from 13 cohorts. We tested the distribution of candidate genes in kidney segments and investigated the effects of keratin-40 (KRT40) on uromodulin processing. RESULTS: Two genome-wide significant signals were identified for uUMOD: a novel locus (P 1.24E-08) over the KRT40 gene coding for KRT40, a type 1 keratin expressed in the kidney, and the UMOD-PDILT locus (P 2.17E-88), with two independent sets of single nucleotide polymorphisms spread over UMOD and PDILT. Two genome-wide significant signals for uUCR were identified at the UMOD-PDILT locus and at the novel WDR72 locus previously associated with kidney function. The effect sizes for rs8067385, the index single nucleotide polymorphism in the KRT40 locus, were similar for both uUMOD and uUCR. KRT40 colocalized with uromodulin and modulating its expression in thick ascending limb (TAL) cells affected uromodulin processing and excretion. CONCLUSIONS: Common variants in KRT40, WDR72, UMOD, and PDILT associate with the levels of uromodulin in urine. The expression of KRT40 affects uromodulin processing in TAL cells. These results, although limited by lack of replication, provide insights into the biology of uromodulin, the role of keratins in the kidney, and the influence of the UMOD-PDILT locus on kidney function.

Genetic loci and prioritization of genes for kidney function decline derived from a meta-analysis of 62 longitudinal genome-wide association studies.

Estimated glomerular filtration rate (eGFR) reflects kidney function. Progressive eGFR-decline can lead to kidney failure, necessitating dialysis or transplantation. Hundreds of loci from genome-wide association studies (GWAS) for eGFR help explain population cross section variability. Since the contribution of these or other loci to eGFR-decline remains largely unknown, we derived GWAS for annual eGFR-decline and meta-analyzed 62 longitudinal studies with eGFR assessed twice over time in all 343,339 individuals and in high-risk groups. We also explored different covariate adjustment. Twelve genome-wide significant independent variants for eGFR-decline unadjusted or adjusted for eGFR-baseline (11 novel, one known for this phenotype), including nine variants robustly associated across models were identified. All loci for eGFR-decline were known for cross-sectional eGFR and thus distinguished a subgroup of eGFR loci. Seven of the nine variants showed variant-by-age interaction on eGFR cross section (further about 350,000 individuals), which linked genetic associations for eGFR-decline with age-dependency of genetic cross-section associations. Clinically important were two to four-fold greater genetic effects on eGFR-decline in high-risk subgroups. Five variants associated also with chronic kidney disease progression mapped to genes with functional in-silico evidence (UMOD, SPATA7, GALNTL5, TPPP). An unfavorable versus favorable nine-variant genetic profile showed increased risk odds ratios of 1.35 for kidney failure (95% confidence intervals 1.03-1.77) and 1.27 for acute kidney injury (95% confidence intervals 1.08-1.50) in over 2000 cases each, with matched controls). Thus, we provide a large data resource, genetic loci, and prioritized genes for kidney function decline, which help inform drug development pipelines revealing important insights into the age-dependency of kidney function genetics.

Assessment of significance of conditionally independent GWAS signals.

MOTIVATION: Multiple independently associated SNPs within a linkage disequilibrium region are a common phenomenon. Conditional analysis has been successful in identifying secondary signals. While conditional association tests are limited to specific genomic regions, they are benchmarked with genome-wide scale criterion, a conservative strategy. Within the weighted hypothesis testing framework, we developed a 'quasi-adaptive' method that uses the pairwise correlation (r2) and physical distance (d) from the index association to construct priority functions G =G(r2, d), which assign an SNP-specific α-threshold to each SNP. Family-wise error rate (FWER) and power of the approach were evaluated via simulations based on real GWAS data. We compared a series of different G-functions. RESULTS: Simulations under the null hypothesis on 1,100 primary SNPs confirmed appropriate empirical FWER for all G-functions. A G-function with optimal r2 = 0.3 between index and secondary SNP which down-weighted SNPs at higher distance step-wise-strong and gave more emphasis on d than on r2 had overall best power. It also gave the best results in application to the real datasets. As a proof of concept, 'quasi-adaptive' method was applied to GWAS on free thyroxine (FT4), inflammatory bowel disease (IBD) and human height. Application of the algorithm revealed 5 secondary signals in our example GWAS on FT4, 5 secondary signals in case of the IBD and 19 secondary signals on human height, that would have gone undetected with the established genome-wide threshold (α=5×10-8). AVAILABILITY AND IMPLEMENTATION: https://github.com/sghasemi64/Secondary-Signal. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Self-Reported Medication Use and Urinary Drug Metabolites in the German Chronic Kidney Disease (GCKD) Study.

BACKGROUND: Polypharmacy is common among patients with CKD, but little is known about the urinary excretion of many drugs and their metabolites among patients with CKD. METHODS: To evaluate self-reported medication use in relation to urine drug metabolite levels in a large cohort of patients with CKD, the German Chronic Kidney Disease study, we ascertained self-reported use of 158 substances and 41 medication groups, and coded active ingredients according to the Anatomical Therapeutic Chemical Classification System. We used a nontargeted mass spectrometry-based approach to quantify metabolites in urine; calculated specificity, sensitivity, and accuracy of medication use and corresponding metabolite measurements; and used multivariable regression models to evaluate associations and prescription patterns. RESULTS: Among 4885 participants, there were 108 medication-drug metabolite pairs on the basis of reported medication use and 78 drug metabolites. Accuracy was excellent for measurements of 36 individual substances in which the unchanged drug was measured in urine (median, 98.5%; range, 61.1%-100%). For 66 pairs of substances and their related drug metabolites, median measurement-based specificity and sensitivity were 99.2% (range, 84.0%-100%) and 71.7% (range, 1.2%-100%), respectively. Commonly prescribed medications for hypertension and cardiovascular risk reduction-including angiotensin II receptor blockers, calcium channel blockers, and metoprolol-showed high sensitivity and specificity. Although self-reported use of prescribed analgesics (acetaminophen, ibuprofen) was <3% each, drug metabolite levels indicated higher usage (acetaminophen, 10%-26%; ibuprofen, 10%-18%). CONCLUSIONS: This comprehensive screen of associations between urine drug metabolite levels and self-reported medication use supports the use of pharmacometabolomics to assess medication adherence and prescription patterns in persons with CKD, and indicates under-reported use of medications available over the counter, such as analgesics.

Meta-analysis uncovers genome-wide significant variants for rapid kidney function decline.

Rapid decline of glomerular filtration rate estimated from creatinine (eGFRcrea) is associated with severe clinical endpoints. In contrast to cross-sectionally assessed eGFRcrea, the genetic basis for rapid eGFRcrea decline is largely unknown. To help define this, we meta-analyzed 42 genome-wide association studies from the Chronic Kidney Diseases Genetics Consortium and United Kingdom Biobank to identify genetic loci for rapid eGFRcrea decline. Two definitions of eGFRcrea decline were used: 3 mL/min/1.73m2/year or more ("Rapid3"; encompassing 34,874 cases, 107,090 controls) and eGFRcrea decline 25% or more and eGFRcrea under 60 mL/min/1.73m2 at follow-up among those with eGFRcrea 60 mL/min/1.73m2 or more at baseline ("CKDi25"; encompassing 19,901 cases, 175,244 controls). Seven independent variants were identified across six loci for Rapid3 and/or CKDi25: consisting of five variants at four loci with genome-wide significance (near UMOD-PDILT (2), PRKAG2, WDR72, OR2S2) and two variants among 265 known eGFRcrea variants (near GATM, LARP4B). All these loci were novel for Rapid3 and/or CKDi25 and our bioinformatic follow-up prioritized variants and genes underneath these loci. The OR2S2 locus is novel for any eGFRcrea trait including interesting candidates. For the five genome-wide significant lead variants, we found supporting effects for annual change in blood urea nitrogen or cystatin-based eGFR, but not for GATM or LARP4B. Individuals at high compared to those at low genetic risk (8-14 vs. 0-5 adverse alleles) had a 1.20-fold increased risk of acute kidney injury (95% confidence interval 1.08-1.33). Thus, our identified loci for rapid kidney function decline may help prioritize therapeutic targets and identify mechanisms and individuals at risk for sustained deterioration of kidney function.

Integration of GWAS Summary Statistics and Gene Expression Reveals Target Cell Types Underlying Kidney Function Traits.

BACKGROUND: Genetic variants identified in genome-wide association studies (GWAS) are often not specific enough to reveal complex underlying physiology. By integrating RNA-seq data and GWAS summary statistics, novel computational methods allow unbiased identification of trait-relevant tissues and cell types. METHODS: The CKDGen consortium provided GWAS summary data for eGFR, urinary albumin-creatinine ratio (UACR), BUN, and serum urate. Genotype-Tissue Expression Project (GTEx) RNA-seq data were used to construct the top 10% specifically expressed genes for each of 53 tissues followed by linkage disequilibrium (LD) score-based enrichment testing for each trait. Similar procedures were performed for five kidney single-cell RNA-seq datasets from humans and mice and for a microdissected tubule RNA-seq dataset from rat. Gene set enrichment analyses were also conducted for genes implicated in Mendelian kidney diseases. RESULTS: Across 53 tissues, genes in kidney function-associated GWAS loci were enriched in kidney (P=9.1E-8 for eGFR; P=1.2E-5 for urate) and liver (P=6.8·10-5 for eGFR). In the kidney, proximal tubule was enriched in humans (P=8.5E-5 for eGFR; P=7.8E-6 for urate) and mice (P=0.0003 for eGFR; P=0.0002 for urate) and confirmed as the primary cell type in microdissected tubules and organoids. Gene set enrichment analysis supported this and showed enrichment of genes implicated in monogenic glomerular diseases in podocytes. A systematic approach generated a comprehensive list of GWAS genes prioritized by cell type-specific expression. CONCLUSIONS: Integration of GWAS statistics of kidney function traits and gene expression data identified relevant tissues and cell types, as a basis for further mechanistic studies to understand GWAS loci.

Genetic associations of hemoglobin in children with chronic kidney disease in the PediGFR Consortium.

BACKGROUND: Genome-wide association studies (GWAS) in healthy populations have identified variants associated with erythrocyte traits, but genetic causes of hemoglobin variation in children with CKD are incompletely understood. METHODS: The Pediatric Investigation of Genetic Factors Linked with Renal Progression (PediGFR) Consortium comprises three pediatric CKD cohorts: Chronic Kidney Disease in Children (CKiD), Effect of Strict Blood Pressure Control and ACE Inhibition on the Progression of CRF in Pediatric Patients (ESCAPE), and Cardiovascular Comorbidity in Children with CKD (4C). We performed cross-sectional and longitudinal association studies of single-nucleotide polymorphisms (SNPs) in 1125 patients. RESULTS: Children of European (n = 725) or Turkish (n = 400) ancestry (EA or TA) were included. In cross-sectional analysis, two SNPs (rs10758658 and rs12718597) previously associated with RBC traits were significantly associated with hemoglobin levels in children of EA and TA. In longitudinal analysis, SNP rs2540917 was nominally associated with hemoglobin in EA and TA children. CONCLUSIONS: SNPs associated with erythrocyte traits in healthy populations were marginally significant for an association with hemoglobin. Further analyses/replication studies are needed in larger CKD cohorts to investigate SNPs of unknown significance associated with hemoglobin. Functional studies will be required to confirm that the observed associations between SNPs and clinical phenotype are causal.

Genome-Wide Association Studies of Metabolites in Patients with CKD Identify Multiple Loci and Illuminate Tubular Transport Mechanisms.

Background The kidneys have a central role in the generation, turnover, transport, and excretion of metabolites, and these functions can be altered in CKD. Genetic studies of metabolite concentrations can identify proteins performing these functions.Methods We conducted genome-wide association studies and aggregate rare variant tests of the concentrations of 139 serum metabolites and 41 urine metabolites, as well as their pairwise ratios and fractional excretions in up to 1168 patients with CKD.Results After correction for multiple testing, genome-wide significant associations were detected for 25 serum metabolites, two urine metabolites, and 259 serum and 14 urinary metabolite ratios. These included associations already known from population-based studies. Additional findings included an association for the uremic toxin putrescine and variants upstream of an enzyme catalyzing the oxidative deamination of polyamines (AOC1, P-min=2.4×10-12), a relatively high carrier frequency (2%) for rare deleterious missense variants in ACADM that are collectively associated with serum ratios of medium-chain acylcarnitines (P-burden=6.6×10-16), and associations of a common variant in SLC7A9 with several ratios of lysine to neutral amino acids in urine, including the lysine/glutamine ratio (P=2.2×10-23). The associations of this SLC7A9 variant with ratios of lysine to specific neutral amino acids were much stronger than the association with lysine concentration alone. This finding is consistent with SLC7A9 functioning as an exchanger of urinary cationic amino acids against specific intracellular neutral amino acids at the apical membrane of proximal tubular cells.Conclusions Metabolomic indices of specific kidney functions in genetic studies may provide insight into human renal physiology.

X-Linked Glomerulopathy Due to COL4A5 Founder Variant.

Alport syndrome is a rare hereditary disorder caused by rare variants in 1 of 3 genes encoding for type IV collagen. Rare variants in COL4A5 on chromosome Xq22 cause X-linked Alport syndrome, which accounts for ∼80% of the cases. Alport syndrome has a variable clinical presentation, including progressive kidney failure, hearing loss, and ocular defects. Exome sequencing performed in 2 affected related males with an undefined X-linked glomerulopathy characterized by global and segmental glomerulosclerosis, mesangial hypercellularity, and vague basement membrane immune complex deposition revealed a COL4A5 sequence variant, a substitution of a thymine by a guanine at nucleotide 665 (c.T665G; rs281874761) of the coding DNA predicted to lead to a cysteine to phenylalanine substitution at amino acid 222, which was not seen in databases cataloguing natural human genetic variation, including dbSNP138, 1000 Genomes Project release version 01-11-2004, Exome Sequencing Project 21-06-2014, or ExAC 01-11-2014. Review of the literature identified 2 additional families with the same COL4A5 variant leading to similar atypical histopathologic features, suggesting a unique pathologic mechanism initiated by this specific rare variant. Homology modeling suggests that the substitution alters the structural and dynamic properties of the type IV collagen trimer. Genetic analysis comparing members of the 3 families indicated a distant relationship with a shared haplotype, implying a founder effect.

Genetic risk variants for membranous nephropathy: extension of and association with other chronic kidney disease aetiologies.

Background: Membranous nephropathy (MN) is a common cause of nephrotic syndrome in adults. Previous genome-wide association studies (GWAS) of 300 000 genotyped variants identified MN-associated loci at HLA-DQA1 and PLA2R1. Methods: We used a combined approach of genotype imputation, GWAS, human leucocyte antigen (HLA) imputation and extension to other aetiologies of chronic kidney disease (CKD) to investigate genetic MN risk variants more comprehensively. GWAS using 9 million high-quality imputed genotypes and classical HLA alleles were conducted for 323 MN European-ancestry cases and 345 controls. Additionally, 4960 patients with different CKD aetiologies in the German Chronic Kidney Disease (GCKD) study were genotyped for risk variants at HLA-DQA1 and PLA2R1. Results: In GWAS, lead variants in known loci [rs9272729, HLA-DQA1, odds ratio (OR) = 7.3 per risk allele, P = 5.9 × 10-27 and rs17830558, PLA2R1, OR = 2.2, P = 1.9 × 10-8] were significantly associated with MN. No novel signals emerged in GWAS of X-chromosomal variants or in sex-specific analyses. Classical HLA alleles (DRB1*0301-DQA1*0501-DQB1*0201 haplotype) were associated with MN but provided little additional information beyond rs9272729. Associations were replicated in 137 GCKD patients with MN (HLA-DQA1: P = 6.4 × 10-24; PLA2R1: P = 5.0 × 10-4). MN risk increased steeply for patients with high-risk genotype combinations (OR > 79). While genetic variation in PLA2R1 exclusively associated with MN across 19 CKD aetiologies, the HLA-DQA1 risk allele was also associated with lupus nephritis (P = 2.8 × 10-6), type 1 diabetic nephropathy (P = 6.9 × 10-5) and focal segmental glomerulosclerosis (P = 5.1 × 10-5), but not with immunoglobulin A nephropathy. Conclusions: PLA2R1 and HLA-DQA1 are the predominant risk loci for MN detected by GWAS. While HLA-DQA1 risk variants show an association with other CKD aetiologies, PLA2R1 variants are specific to MN.

Genetic loci associated with renal function measures and chronic kidney disease in children: the Pediatric Investigation for Genetic Factors Linked with Renal Progression Consortium.

BACKGROUND: Chronic kidney disease (CKD) in children is characterized by rapid progression and a high incidence of end-stage renal disease and therefore constitutes an important health problem. While unbiased genetic screens have identified common risk variants influencing renal function and CKD in adults, the presence and identity of such variants in pediatric CKD are unknown. METHODS: The international Pediatric Investigation for Genetic Factors Linked with Renal Progression (PediGFR) Consortium comprises three pediatric CKD cohorts: Chronic Kidney Disease in Children (CKiD), Effect of Strict Blood Pressure Control and ACE Inhibition on the Progression of CRF in Pediatric Patients (ESCAPE) and Cardiovascular Comorbidity in Children with CKD (4C). Clean genotype data from > 10 million genotyped or imputed single-nucleotide polymorphisms (SNPs) were available for 1136 patients with measurements of serum creatinine at study enrollment. Genome-wide association studies were conducted to relate the SNPs to creatinine-based estimated glomerular filtration rate (eGFR crea) and proteinuria (urinary albumin- or protein-to-creatinine ratio ≥ 300 and ≥ 500 mg/g, respectively). In addition, European-ancestry PediGFR patients (cases) were compared with 1347 European-ancestry children without kidney disease (controls) to identify genetic variants associated with the presence of CKD. RESULTS: SNPs with suggestive association P-values < 1 × 10(-5) were identified in 10 regions for eGFR crea, four regions for proteinuria and six regions for CKD including some plausible biological candidates. No SNP was associated at genome-wide significance (P < 5 × 10(-8)). Investigation of the candidate genes for proteinuria in adults from the general population provided support for a region on chromosome 15 near RSL24D1/UNC13C/RAB27A. Conversely, targeted investigation of genes harboring GFR-associated variants in adults from the general population did not reveal significantly associated SNPs in children with CKD. CONCLUSIONS: Our findings suggest that larger collaborative efforts will be needed to draw reliable conclusions about the presence and identity of common variants associated with eGFR, proteinuria and CKD in pediatric populations.

Genome-wide association studies in pediatric chronic kidney disease.

The genome-wide association study (GWAS) has become an established scientific method that provides an unbiased screen for genetic loci potentially associated with phenotypes of clinical interest, such as chronic kidney disease (CKD). Thus, GWAS provides opportunities to gain new perspectives regarding the genetic architecture of CKD progression by identifying new candidate genes and targets for intervention. As such, it has become an important arm of translational science providing a complementary line of investigation to identify novel therapeutics to treat CKD. In this review, we describe the method and the challenges of performing GWAS in the pediatric CKD population. We also provide an overview of successful GWAS for kidney disease, and we discuss the established pediatric CKD cohorts in North America and Europe that are poised to identify genetic risk variants associated with CKD progression.

Genome-wide association studies in nephrology: using known associations for data checks.

Prior to conducting genome-wide association studies (GWAS) of renal traits and diseases, systematic checks to ensure data integrity and analytical work flow should be conducted. Using positive controls (ie, known associations between a single-nucleotide polymorphism [SNP] and a corresponding trait) allows for identifying errors that are not apparent solely from global evaluation of summary statistics. Strong genetic control associations of chronic kidney disease (CKD), as derived from GWAS, are lacking in the non-African ancestry CKD population; thus, in this perspective, we provide examples of and considerations for using positive controls among patients with CKD. Using data from individuals with CKD who participated in the CRIC (Chronic Renal Insufficiency Cohort) Study or PediGFR (Pediatric Investigation for Genetic Factors Linked to Renal Progression) Consortium, we evaluated 2 kinds of positive control traits: traits unrelated to kidney function (bilirubin level and body height) and those related to kidney function (cystatin C and urate levels). For the former, the proportion of variance in the control trait that is explained by the control SNP is the main determinant of the strength of the observable association, irrespective of adjustment for kidney function. For the latter, adjustment for kidney function can be effective in uncovering known associations among patients with CKD. For instance, in 1,092 participants in the PediGFR Consortium, the P value for the association of cystatin C concentrations and rs911119 in the CST3 gene decreased from 2.7×10(-3) to 2.4×10(-8) upon adjustment for serum creatinine-based estimated glomerular filtration rate. In this perspective, we give recommendations for the appropriate selection of control traits and SNPs that can be used for data checks prior to conducting GWAS among patients with CKD.