Genome-wide association study of estimated glomerular filtration rate using repeated measurements in the Taiwan Biobank

Background: Chronic kidney disease (CKD) is a major global public health issue, with genetic factors playing a significant role in kidney function. Although genome-wide association studies (GWAS) have identified numerous loci associated with estimated glomerular filtration rate (eGFR), most studies relied on a single time-point measurement, which limits the capacity to account for within-individual measurement variability. Methods: We performed a repeated-measurement GWAS in the prospective Taiwan Biobank (Taiwanese ancestry; n = 25,004) using two repeated creatinine-based eGFR measurements. Repeated eGFR values were analyzed using a linear mixed-effects model with a subject-specific random intercept and time-varying covariates, providing a more precise estimate of eGFR level. Identified loci underwent functional annotation (expression quantitative trait locus, deleteriousness prediction, and epigenetic markers) and were compared with results from a single-measurement GWAS. Results: Six loci associated with eGFR were identified, including four previously reported regions (1q22, 4q21.1, 11p14.1, and 17q21.2) and two additional loci (6p21.32 and 15q24.2). Functional annotation implicated several candidate genes-such as MUC1/EFNA1, SHROOM3, HLA-DQB1, MPPED2, NRG4, and PGAP3/FBXL20-in the regulation of kidney function. Conclusion: Incorporating repeated eGFR measurements into GWAS may improve phenotypic precision for identifying genetic associations with kidney function. This study identified eGFR-associated loci and biologically plausible candidate genes in a Taiwanese population, which require further replication and functional validation.