Kidney Res Clin Pract > Epub ahead of print
Lee, Moon, Song, Lee, Kim, Cho, Gil, Lee, and Park: Genome-wide association study and fine-mapping on Korean biobank to discover renal trait-associated variants



Chronic kidney disease is a significant health burden worldwide, with increasing incidence. Although several genome-wide association studies (GWAS) have investigated single nucleotide polymorphisms (SNP) associated with kidney trait, most studies were focused on European ancestry.


We utilized clinical and genetic information collected from the Korean Genome and Epidemiology Study (KoGES).


More than five million SNPs from 58,406 participants were analyzed. After meta-GWAS, 1,360 loci associated with estimated glomerular filtration rate (eGFR) at a genome-wide significant level (p = 5 × 10–8) were identified. Among them, 399 loci were validated with at least one other biomarker (blood urea nitrogen [BUN] or eGFRcysC) and 149 loci were validated using both markers. Among them, 18 SNPs (nine known ones and nine novel ones) with 20 putative genes were found. The aggregated effect of genes estimated by MAGMA gene analysis showed that these significant genes were enriched in kidney-associated pathways, with the kidney and liver being the most enriched tissues.


In this study, we conducted GWAS for more than 50,000 Korean individuals and identified several variants associated with kidney traits, including eGFR, BUN, and eGFRcysC. We also investigated functions of relevant genes using computational methods to define putative causal variants.


Chronic kidney disease (CKD) is a significant health issue with a globally increasing incidence, affecting over 850 million individuals with kidney diseases worldwide [1]. To uncover the pathogenesis of CKD, global consensus initiatives have conducted large-scale genome-wide association studies (GWAS) and meta-analyses of GWAS (meta-GWAS). During the earlier stage of the GWAS era, several loci associated with renal function and kidney disease have been identified [2,3]. These studies have demonstrated the vast potential of GWAS in CKD research. Recently, several fine-mapping studies have been conducted in conjunction with GWAS to perform functional annotations of identified variants, shedding light on the pathogenesis of CKD [48]. However, most GWAS were focused on European ancestry.
The Korean Genome and Epidemiology Study (KoGES) is a nationwide cohort that has collected clinical and genetic information since 2001 [9]. Using this cohort, GWAS have been conducted on various traits, including alcoholic liver disease [10], serum uric acid [11], and muscle mass [12]. Estimated glomerular filtration rate (eGFR) is a widely accepted kidney trait for GWAS. A few studies have been conducted to identify single nuclear polymorphisms (SNP) associated with eGFR in the Korean population [13,14]. However, these studies have reported associated SNPs without functional fine-mapping or functional annotations. Given that SNPs discovered by GWAS are regulatory variants associated with complex traits and diseases, uncovering functional annotation and fine-mapping are important [15]. Therefore, we conducted a GWAS to identify loci associated with eGFR and a fine-mapping study to reveal putative causal SNPs in the Korean populations.


Participants and genotyping

This study protocol was reviewed and approved by the Institutional Review Board (IRB) of Soonchunhyang University Cheonan Hospital (Cheonan, Korea) and the need for informed consent was waived by IRB (No. SCHCA 2021-11-035). This study complied with the principles of the Declaration of Helsinki.
All genotype data were obtained from the Korea Biobank Array Project managed by the Korea National Institute of Health. All samples and clinical data were collected as part of the KoGES [9]. The cohort comprised three sub-cohorts, namely the Korean Association Resource (KARE), Health Examinee (HEXA), and Cardiovascular Disease Association Study (CAVAS). The methodology for genotyping across all cohorts has been detailed elsewhere [16]. Shortly, these cohorts comprised community-dwellers aged ≥40 years at baseline. The KARE cohort, also known as the Ansan and Anseong cohort, included participants from the Ansan and Anseong regions representing urban and rural areas, respectively. The HEXA study recruited participants who attended regular health check-ups in urban areas. The CAVAS cohort had participants from rural areas to investigate cardiovascular diseases.
The three sub-cohorts consisted of 10,030, 177,357, and 28,338 participants, respectively. Among them, 8,840, 58,694, and 8,105 subjects were genotyped using Affymetrix Genome-Wide Human SNP array version 5.0 (Affymetrix) [17] for KARE and Korea Biobank Array [16] for HEXA and CAVAS. Although the quality control process of genotyping was initially conducted before distributing cohort data to researchers, as previously reported [1618], we additionally performed rigorous quality control processes [19]. These processes involved quality controls based on predetermined criteria, such as SNPs including low genotype calls (<0.01), individuals with high rates of genotype missingness (<0.05), sex discrepancy-based X-chromosome homozygosity, low minor allele frequency (MAF, <0.05), Hardy-Weinberg equilibrium (p < 1 × 10–6), individuals with heterozygosity rate deviated ±3 standard deviation from the mean, cryptic relatedness based on pi-hat threshold of 0.2, and population stratification. Population stratification was analyzed using a multidimensional scaling (MDS) approach for KARE and CAVAS cohorts. However, principal component analysis was used to calculate population stratification in the HEXA cohort due to the large sample size that made MDS analysis unfeasible.

Genotype imputation

After completing quality control, the imputation of all cohorts was carried out using the 1000 Genome Project Phase 3 reference panel which comprised an Asian population [20]. The imputation was performed using Beagle 5.4 and involved haplotype phasing [21] and imputation [22]. Following the imputation process, additional quality control was executed based on imputation information quality score of >0.8 and MAF of >0.05.

Phenotype and covariates

In the HEXA cohort, hemoglobin A1c (HbA1c) was available for a subset (only 54%) of participants at baseline. Due to the limited availability of laboratories in the HEXA cohort at baseline, we utilized data from the first follow-up where almost all individuals had HbA1c results. Hence, the first follow-up data from the HEXA cohort were utilized for this analysis. For the KARE and CAVAS cohorts, baseline data were used for analysis.
The primary phenotype in all cohorts was a quantitative trait, which was defined by eGFR using the Chronic Kidney Disease Epidemiology Collaboration equation [23]. Creatinine levels were measured using a Hitachi Analyzer 7600 (Hitachi) in KARE and an ADVIA 1650 (Siemens Healthcare) in the CAVAS cohort. Serum cystatin C was measured using a Cobas c702 (Roche). All KARE and CAVAS cohort samples were measured at Seoul Clinical Laboratories. It should be noted that the HEXA cohort was based on data from a nationwide health examination where each laboratory result was measured at the institution where the health examination was performed. Thus, specific methods used by each institution were unavailable.
To estimate an unbiased effect of genotype on renal function (i.e., eGFR), subjects with diabetes mellitus (DM) or albuminuria were excluded. As covariates, age, sex, body mass index, systolic blood pressure (SBP), and past medical history of hypertension were used. Past medical history of hypertension was defined as SBP of ≥140 mmHg, diastolic blood pressure of ≥90 mmHg, and self-report by participants. History of DM was defined by HbA1c of ≥6.5%, fasting glucose of ≥126 mg/dL, and self-report by participants. Albuminuria was defined by ≥1+ in the urine dipstick test. Individuals with missing values in covariates were excluded from further analysis. Population structure was visually inspected and found to be homogenous, with all participants having Korean ancestry (Supplementary Fig. 1, available online). Nonetheless, to adjust for population stratification, 10 principal components calculated based on genotypes were included as covariates.
In addition to eGFR, blood urea nitrogen (BUN) and eGFR calculated using serum cystatin C (eGFRcysC) [24] were used for validation, although serum cystatin C was available for the CAVAS cohort. All phenotypes were quantitative traits since the prevalence of CKD based on eGFR criteria was too low to perform a case-control analysis. When the histogram of each variable was normally distributed, original values were used as phenotype. Otherwise, naturally transformed values were used.

Genome-wide association study and meta-analysis

Most GWAS were performed using plink version 1.9 [25]. Quantitative trait was analyzed using linear regression with the plink command “--linear” and prespecified covariates with the assumption of an additive genetic model. An association study was performed for the three cohorts individually, then a meta-GWAS analysis was performed using METAL (Meta-Analysis Tool for Genome-wide Association Scans) with genomic control correction [26]. Genome-wide significant (GWS) level was defined as p of <5 × 10–8. Distribution of observed p-values and estimated p-values of given SNPs were depicted using quantile-quantile (QQ) plots. QQ plots and Manhattan plots were drawn using qqman R packages of R software (R Foundation for Statistical Computing).

Functional mapping and annotation

Following meta-GWAS, variants associated with eGFR with GWS level were selected for validation analysis. To confirm associations of these variants with other biomarkers, results from meta-GWAS of log-transformed BUN (log-BUN) and eGFRcysC were utilized. Variants that had both a directionally opposite, nominally significant association (p < 0.05) with log-BUN and a directionally concordant, nominally significant association with eGFRcysC were defined as validated loci.
We used FUMA (Functional Mapping and Annotation of Genome-Wide Association Studies) to perform functional mapping and annotation with specific purposes: (1) to specify genomic risk loci that were independently significant with other SNPs (i.e., clumping); (2) to annotate genes; and (3) to validate enrichment in tissue expression data [27]. At first, FUMA selected independent significant SNPs that exhibited a GWS association (p < 5 × 10–8) and were independent (r2 < 0.6). Thus, independent significant SNPs were equivalent to SNPs that remained after clumping GWAS-tagged SNPs with the same p-value and r2 threshold. Based on these SNPs, lead SNPs were identified when they were independent of other SNPs with a level of r2 < 0.1. If the linkage disequilibrium (LD) blocks of those independent significant SNPs were closely positioned (within ±250 kb from the first and last of the LD block), they were integrated into a genomic risk region. FUMA defined the lead SNP of the genomic risk region as the genomic risk loci. FUMA used ANNOVAR to annotate and map identified variants [28]. Additionally, we performed MAGMA (Multi-marker Analysis of Genomic Annotation) gene analysis [29] and conditional and interaction gene-set analysis (i.e., gene-property analysis) [30] using FUMA. The gene-property analysis was based on tissue expression data from Genotype-Tissue Expression (GTEx) v8 [31]. Of all 54 GTEx tissues, 49 tissues with a sample size of ≥70 were used.

Gene-set and tissue enrichment analysis

Additional gene-set enrichment analysis was performed using the DOSE R package [32]. To investigate tissue type-specific enrichment based on GWAS summary statistics, stratified LD score regression applying to specifically expressed genes (LDSC-SEG) was performed [33]. Two gene expression datasets from the GTEx project [31] and Franke lab [34], of which annotation data had previously been curated specifically for East Asians by LDSC-SEG builders, were used for LDSC-SEG [35].
When more than two independent significant SNPs were identified and remained after validation with other biomarkers (i.e., log-BUN and eGFRcysC) within a genomic risk region, conditional analysis was performed using a Genome-wide Complex Analysis (GCTA) tool [36]. Loci with p of <5 × 10–8 were retained after conditional analysis for the most significant SNPs (lowest p-value) within a genomic risk region. Genotype data from the HEXA cohort, which had the majority of data, were used as the LD reference for conditional analysis.


Baseline characteristics and phenotypes

After genotype quality control, 8,384, 58,079, and 7,966 subjects remained in the KARE, HEXA, and CAVAS cohorts, respectively. Among the 58,079 subjects in the HEXA cohort, 51,902 visited at the first follow-up schedule. Individuals with missing values and those with DM or albuminuria were then removed. Finally, 6,848, 44,787, and 6,771 individuals remained in the KARE, HEXA, and CAVAS cohorts, respectively (Supplementary Fig. 2, available online).
Table 1 shows the baseline characteristics of the KoGES cohort. Median serum creatinine levels in KARE, HEXA, and CAVAS cohorts were 0.80 (interquartile range [IQR], 0.70–1.00), 0.75 (IQR, 0.66–0.89), and 0.91 (IQR, 0.83–1.03), respectively. The median age of the KARE cohort was lower than that of the HEXA cohort or the CAVAS cohort. As all cohorts targeted the general population, the proportion of CKD patients was low. The highest proportion of CKD, defined by an eGFR of <60 mL/min/1.73 m2, was observed in the CAVAS cohort.
Histograms indicated that eGFR and eGFRcysC of each cohort were normally distributed, while BUN was more likely to be normally distributed when it was naturally log-transformed (Supplementary Fig. 3, available online). Hence, eGFR and natural log-BUN were used as quantitative traits.


After quality control and genotype imputation processes, 4,864,729, 5,741,581, and 5,746,961 variants remained in the KARE, HEXA, and CAVAS cohorts, respectively. Following meta-GWAS analysis, 1,360 loci were found to be associated with eGFR at the GWS level. Manhattan plots of meta-GWAS for eGFR revealed several regions with significant association (outer circle of Fig. 1). Results from the meta-GWAS on log-BUN demonstrated that 2,454 loci were associated at the GWS level (Supplementary Fig. 4, available online). In addition, 143 SNPs were associated with eGFRcysC at the GWS level. Of 1,360 variants associated with eGFR, 399 were validated by at least one biomarker (i.e., log-BUN or eGFRcysC) and 149 were validated by both biomarkers (Supplementary Table 1, available online).
Validated loci (yellow points in Fig. 1) showed a similar pattern to that of a previously reported large GWAS study (inner circle of Fig. 1) [37]. The inner plot in Fig. 1 illustrates 399 validated loci for at least one biomarker. It was worth noting that the effect (beta) of validated variants was more prominent in cases where the mean allele frequency was lower.

Functional mapping and annotation

FUMA identified 87 independent significant SNPs, consisting of lead SNPs and 33 genomic risk regions. Of these, 31 loci were validated by at least one biomarker (log-BUN or eGFRcysC) (Fig. 2A). Fourteen loci were validated by both eGFRcysC and log-BUN. Ten loci were validated by log-BUN and seven loci were validated by eGFRcysC. Scatterplots of effect (beta) between eGFR and other biomarkers are also depicted in Fig. 2 (red dots validated by both log-BUN and eGFRcysC, blue dots validated by log-BUN, and green dots validated by eGFRcysC).
Among the 31 validated loci, 11 lead SNPs were identified within 15 genomic risk regions. When loci within a single genomic risk region had ≥2, a conditional analysis was performed (see Methods section). As a result, 18 loci remained, and three loci (rs4665985, rs62141288, and rs35578578) remained significant after conditional analysis (Table 2).
The median eGFR showed significant differences according to the dosage of loci, i.e., the number of effect alleles (Supplementary Fig. 5, available online). In the HEXA cohorts, all variants in Table 2 showed significant differences between genotypes after Bonferroni correction. Only seven variants (rs1260326, rs33921462, rs744103, rs35578578, rs35449439, rs2240736, and rs549752) were validated as significant in the CAVAS cohort. None was validated in the KARE cohort, although this was likely due to the smaller sample size of the KARE cohort.
We also investigated whether the effect (beta) of variants discovered in our study was concordant with that of a previous large study [37]. Among the 399 loci validated by log-BUN or eGFRcysC, 333 variants were also observed. The degree and direction of the effect of loci seemed to be concordant with the previous report (Fig. 3). Of 18 variants shown in Table 2, 12 loci were noted.

Gene-set and tissue enrichment study

MAGMA gene analysis was also performed using FUMA, resulting in 65 genes that remained statistically significant after Bonferroni correction (Pbon < 0.05). Among these genes, 33 were located within the genomic risk region (Fig. 4A). Gene-property analysis based on GTEx v8 revealed that the kidney cortex was the most enriched tissue (Fig. 4B), out of the 49 GTEx tissues (nominal p = 0.002).
To further investigate the significance of these genes, overrepresentation analysis was performed using the disease gene network (DieGeNet) [38] for genes that remained significant after Bonferroni correction from the MAGMA gene analysis (Fig. 4C, D). Overrepresented genes were enriched in kidney-associated terms, such as glomerular filtration rate, creatinine measurement, uric acid measurement, and BUN measurements (Fig. 4C). A network plot of enriched terms revealed several genes shown in Table 2 as well as known genes associated with renal traits, such as TBX2.
To investigate tissue enrichment, heritability enrichment using LDSC-SEG was performed (see Methods section). The liver was the most significantly expressed tissue in GTEx (false discovery rate [FDR], <0.05) (Fig. 5A), while the kidney cortex showed nominal significance (nominal p = 0.037, FDR, 0.098). Tissue enrichment in kidneys was significant when using the Franke lab dataset (Fig. 5B).


In this study, we conducted a GWAS for the trait of eGFR using a Korea Biobank Array on the Korean population. There were some previous GWAS studies for Koreans. However, one study has only involved the KARE cohort of KoGES [13]. Another study has reported GWAS results only for selected patients [14]. In addition, previous Korean GWAS studies did not perform fine-mapping analysis. We attempted to define functional annotation and fine-mapping and validate discovered loci in this study by comparing them with those found in other previous studies. We identified 18 SNPs (nine novel ones and nine previously reported ones) across 15 genomic risk regions spanning 20 genes (Table 2). Furthermore, our GWAS results demonstrated statistically significant enrichment in kidney-related diseases, as confirmed in the DieGeNet database. We also observed enrichment in kidney and liver tissues, in agreement with a previous study [37].
Among genes found in this study, some were concordant with previous reports and their functional roles were also demonstrated. For example, LRP2, also known as megalin, has been reported to be a target molecule associated with anti-brush border antibodies and renal failure (ABBA disease) [39]. SHROOM3 is well known as a GFR-associated gene [2,37]. It is associated with the development of kidneys in an animal model [40]. A recent study has shown the role of lysosomal beta-mannosidase (MANBA) expression in kidney disease. Manba and Ubed2d were expressed in kidney tubule cells and fibroblasts, respectively [41]. The association of SLC34A1 with kidney disease has already been reported [3]. Given that mutations in SLC34A1 are associated with nephrolithiasis, the association between genotype and phenotype might be due to renal stone-related damage to the kidney [42]. PDILT is known to be located near UMOD and to regulate uromodulin expression [6]. NFATC1 is associated with tumor necrosis factor-associated podocyte injury by NFATC1/ABCA1-dependent mechanism [43]. The role of WDR72 remains unclear, although it has been reported in previous a meta-GWAS study [4]. Recent work showed WDR72 might have a role associated with uromodulin along with the UMOD-PDILT locus [44].
We searched whether 18 SNPs discovered in this study were available in the GWAS catalog [45]. Only nine SNPs (rs1260326, rs4665985, rs223471, rs33921462, rs744103, rs58063923, rs11066132, rs2240736, and rs549752) could be found in the GWAS catalog (searched in February 2023). This was quite different from what was expected, while previous GWAS discovered many loci associated with kidney traits. This might be attributed to the fact that previous studies were almost based on European ancestry. Indeed, the population structure in the KoGES cohort was quite different from those of other populations (Supplementary Fig. 1, available online). Hence, we examined the possibility of novel SNPs being in LD with previously reported SNPs, excluding rs62141288 and rs35578578, which were determined to be conditional on rs1260326 and rs744103, respectively (both rs1260326 and rs744103 were previously identified as kidney-associated SNPs). Of the remaining SNPs (n = 7), all were found to be in nearly complete LD with SNPs that were previously known to be associated with kidneys (Supplementary Fig. 6, available online). In addition, we conducted an investigation on the associations between nine SNPs that were not previously identified in the GWAS catalog and kidney-related traits using the BioBank Japan database [46]. Among them, six SNPs (rs3790638, rs62141288, rs77366165, rs7677847, rs572528, and rs35449439) exhibited a significant association (p < 5 × 10–8) with kidney-related traits (Supplementary Table 2, available online). Given the similarity in ethnicity between Koreans and Japanese, the novel SNPs discovered in this study may be extrapolated as specific to the northeastern Asian population.
Our study has a limitation. It only involved Koreans. Thus, significant loci were different from previous studies and the GWAS catalog. Ethnic differences in population structure and genetic architecture might have influenced our results as described above. Most loci discovered by GWAS were known to be non-coding variants that might exert regulatory functions [47]. Gene expression varies according to their ancestry [48]. Given that regulatory variants might be different between ancestries, lead SNPs associated with the trait (i.e., eGFR) were different from other studies based on most of the European population. However, enrichment analysis showed that tissue-specific expression was not significantly different from the previous report. Target genes were expressed in a kidney-specific manner (Fig. 4, 5), although the concrete pattern of expression might differ slightly from results from European ancestry. Therefore, Korean eQTL data for other diseases should be curated in the future. In addition, the sample size in our study was not sufficient to draw robust conclusions. Particularly, the SNPs we discovered did not demonstrate significance in the KARE or CAVAS cohorts, which increases the risk of false positives. The number of participants needs to be increased in future studies.
In conclusion, we discovered several SNPs associated with kidney traits in the Korean population based on KoGES, the largest cohort in Korea. We also discovered that variants were validated in other enrichment analyses. Although a detailed causality and associated mechanisms should be elucidated in the future, we found not only concordant results with previous GWAS but also novel loci that might be specific to the Korean population.


Conflicts of interest

All authors have no conflicts of interest to declare.


This work was supported by the Basic Science Research Program (NRF-2022R1F1A1071128 to Samel Park and 2021R1C1C1007810 to Jong-Seok Moon) from the National Research Foundation of Korea.

Data sharing statement

The data presented in this study was obtained from the Korean Genome-wide Epidemiology Study (KoGES). It could be requested on the website (

Authors’ contributions

Conceptualization, Formal analysis, Visualization: SP

Data curation: DSK, NJC

Funding acquisition, Project administration: JSM, SP

Investigation: DJL, SP

Methodology: DKS, YSL

Resources: DJL, NJC

Software: DJL, DSK

Supervision: JSM, YSL, HWG, EYL

Writing–original draft: DJL, SP

Writing–review & editing: DJL, SP

All authors read and approved the final manuscript.


This study was conducted with bioresources from the National Biobank of Korea, the Korea Disease Control and Prevention Agency, Republic of Korea (NBK-21120702-01-01).

Figure 1.

Circos plot for GWAS meta-analysis.

The outer circle depicts the loci associated with eGFR in KoGES, while the inner circle represents those in a previously published GWAS by Stanzick et al. [37]. Yellow points indicate validated loci by at least one biomarker (log-BUN or eGFRcysC). The Y-axis denotes –log10(p) for association with eGFR. In the inner circle, the Y-axis was truncated at 40. Red dotted lines indicate a genome-wide significant level (p = 5 × 10–8). The inner plot illustrates a correlation between beta and mean allele frequency of validated loci (red points, lead SNPs; blue points, independent significant SNPs; gray points, validated by both log-BUN and eGFRcysC).
BUN, blood urine nitrogen; eGFR, estimated glomerular filtration rate calculated using serum creatinine by CKD-EPI equation; eGFRcysC, eGFR calculated using serum cystatin C by CKD-EPI equation; GWAS, genome-wide association study; KoGES, the Korean Genome and Epidemiology Study; SNP, single nucleotide polymorphism.
Figure 2.

Independent significant SNPs validated by at least one biomarker, log-BUN or eGFRcysC.

FUMA identified 87 independent significant SNPs. (A) Among them, 14 SNPs were validated by both log-BUN and eGFRcysC, whereas 56 were validated by none of them. The concordance between eGFR and log-BUN (B) and that between eGFR and eGFRcysC (C) are represented. Red dots correspond to independent significant SNPs validated by both log-BUN and eGFRcysC, while blue dots represent loci validated by log-BUN (B) and green dots represent those validated by eGFRcysC (C).
BUN, blood urine nitrogen; eGFR, estimated glomerular filtration rate calculated using serum creatinine by CKD-EPI equation; eGFRcysC, eGFR calculated using serum cystatin C by CKD-EPI equation; SNP, single nucleotide polymorphism.
Figure 3.

Scatter plot for concordance between studies.

Of 399 loci validated by log-BUN or eGFRcysC, 333 variants were also identified in results from a meta-GWAS by Stanzick et al. [37]. The x- and y-axis shows beta and 95% confidence interval from our study (KoGES) and a previous large meta-GWAS (by Stanzick et al.), respectively. Gray dots represent validated only loci. Blue dots indicate validated independent significant single nucleotide polymorphisms. Red dots represent validated genomic significant loci.
BUN, blood urine nitrogen; eGFR, estimated glomerular filtration rate calculated using serum creatinine by CKD-EPI equation; eGFRcysC, eGFR calculated using serum cystatin C by CKD-EPI equation; GWAS, genome-wide association study; KoGES, the Korean Genome and Epidemiology Study.
Figure 4.

MAGMA gene analysis and pathway analysis based on significant genes.

MAGMA gene analysis identified significant genes based on GWAS summary statistics. (A) Manhattan plots showing 33 genes located within the genomic risk region that were significant after Bonferroni correction (pbon < 0.05). (B) To investigate enrichment in tissues, gene-property analysis was performed. The top 20 enrichment tissues (out of 54 GTEx tissues, with 49 tissues having sample sizes above 70) are depicted to be ordered according to their significance. The dashed line indicates pbon < 0.05 and the solid line indicates nominal p < 0.05. (C, D) Overrepresentation analysis based on disease gene network (DieGeNet).
GTEx, Genotype-Tissue Expression; GWAS, genome-wide association study.
Figure 5.

Heritability enrichment in tissues from GTEx and Franke using LDSC-SEG.

Tissue-specific enrichment was performed using gene expression dataset from GTEx (A) and from Franke lab (B). The red dashed line indicates significance after Bonferroni correction (pbon < 0.05). The blue line represents nominal significance (p < 0.05).
LDSC-SEG, linkage disequilibrium score regression applied to specifically expressed genes; GTEx, Genotype-Tissue Expression.
Table 1.
Baseline characteristics of HEXA, KARE, and CAVAS cohorts
Characteristic KARE HEXA CAVAS
No. of patients 6,848 44,787 6,771
Age (yr) 49 (44–59) 58 (52–64) 59 (51–65)
 40–65 6,281 (91.7) 35,414 (79.1) 5,137 (75.9)
 >65 567 (8.3) 9,373 (20.9) 1,634 (24.1)
Male sex 3,163 (46.2) 14,796 (33.0) 2,518 (37.2)
Hypertension 1,949 (28.5) 8,252 (18.4) 2,039 (30.1)
SBP (mmHg) 118 (107–130) 121 (112–132) 122 (111–134)
DBP (mmHg) 79 (71–87) 74 (68–80) 78 (70–85)
BMI (kg/m2) 24.3 (22.4–26.3) 23.5 (21.8–25.4) 24.3 (22.4–26.3)
Waist circumference (cm) 82.0 (75.7–87.9) 80.5 (74.8–86.2) 84.0 (78.0–89.8)
Hemoglobin (g/dL) 13.5 (12.5–14.7) 13.8 (13.0–14.8) 13.7 (12.8–14.7)
Albumin (g/dL) 4.5 (4.3–4.7) 4.6 (4.5–4.8) 4.5 (4.3–4.6)
Glucose (mg/dL) 87 (82–92) 95 (89–101) 92 (87–99)
HbA1c (%) 5.5 (5.3–5.8) 5.5 (5.2–5.7) 5.5 (5.2–5.7)
AST (IU/L) 22.0 (19.0–28.0) 23.0 (20.0–27.0) 24.0 (21.0–28.0)
ALT (IU/L) 19.0 (14.0–26.0) 19.0 (15.0–24.0) 20.0 (16.0–27.0)
BUN (mg/dL) 13.4 (11.2–16.0) 14.5 (12.2–17.1) 15.2 (12.7–18.0)
Creatinine (mg/dL) 0.80 (0.70–1.00) 0.75 (0.66–0.89) 0.91 (0.83–1.03)
eGFR (mL/min/1.73m2) 95.5 (81.8–105.4) 92.7 (83.4–99.6) 76.1 (68.9–83.4)
eGFRcysC (mL/min/1.73m2) 91.0 (78.0–102.2)
eGFR < 60 mL/min 125 (1.8) 865 (1.9) 491 (7.3)

Data are expressed as number only, median (interquartile range ), or number (%).

There were some missing values: 1) 6 for waist circumference, 24 for glucose level, 1 for HbA1c, and 3 for ALT in KARE cohort; 2) 2 for DBP, 7 for waist circumference, 5 for hemoglobin, 5 for HbA1c, 1 for AST, 29 for ALT in HEXA cohort; 3) 14 for waist circumference, 475 for HbA1c, 5 for ALT, and 80 for cystatin C in CAVAS cohort.

ALT, alanine transferase; AST, aspartate transferase; BUN, blood urine nitrogen; BMI, body mass index; CAVAS, Cardiovascular Disease Association Study; DBP, diastolic blood pressure; eGFR, estimated glomerular filtration rate calculated using serum creatinine by CKD-EPI equation; eGFRcysC, estimated glomerular filtration rate calculated using serum cystatin C by CKD-EPI equation; HbA1c, glycated hemoglobin; HEXA, Health Examinee; KARE, Korean Association Resource; SBP, systolic blood pressure.

Table 2.
Validated SNPs within each genomic risk region and mapped gene list
CHR Genomic risk region rsID Position (b37) EA/NEA (EAF) BETA p-value Cond P Mapped gene lista
1 1:10,702,266–10,718,377 rs3790638 10,707,812 A/G (0.05) 0.84 3.6 × 10–8 RP4-734G22.3
2 2:27,598,097–27,844,601 rs1260326b 27,730,940 C/T (0.45) 0.56 9.6 × 10–17 GCKR
rs4665985c 27,753,878 A/C (0.47) –0.47 2.2 × 10–12 2.4 × 10–37 GCKR, AC109829.1
rs62141288 27,783,198 A/G (0.33) –0.41 7.6 × 10–9 7.6 × 10–25 AC109829.1
2 2:170,165,283–170,206,062 rs77366165 170,170,804 A/G (0.10) 0.90 3.9 × 10–15 LRP2
4 4:77,363,639–77,414,988 rs7677847 77,364,126 A/G (0.35) –0.39 2.6 × 10–8 SHROOM3
4 4:103,675,108–103,954,851 rs223471b 103,698,786 C/G (0.36) 0.42 2.9 × 10–9 MANBA, UBE2D3
5 5:176,757,841–176,842,474 rs33921462b 176,814,656 A/G (0.31) –0.70 8.0 × 10–21 SLC34A1
6 6:43,804,571–43,829,941 rs744103b 43,805,362 T/A (0.14) –0.63 2.3 × 10–10 VEGFA, RP11-344J7.2
rs35578578 43,810,526 G/GC (0.10) 0.65 4.5 × 10–9 6.6 × 10–25 VEGFA, RP11-344J7.2
7 7:1,243,525–1,299,800 rs58063923c 1,270,738 T/C (0.33) 0.41 3.7 × 10–8 AC091729.9, UNCX
10 10:847,688–1,081,293 rs17159964 913,064 T/G (0.09) –0.88 3.2 × 10–14 LARP4B
11 11:30,749,090–30,777,790 rs56870952 30,750,092 T/TACAAAACAAA (0.33) –0.58 2.9 × 10–16 RP5-1024C24.1, DCDC1
12 12:111,301,027–113,117,897 rs11066132b 112,468,206 T/C (0.16) –0.84 9.0 × 10–20 NAA25
15 15:53,882,419–54,006,275 rs572528 53,972,484 A/G (0.38) 0.50 2.9 × 10–13 WDR72
16 16:20,383,049–20,407,196 rs35449439 20,385,182 C/G (0.18) 0.51 3.0 × 10–8 PDILT
17 17:59,239,149–59,498,250 rs2240736b 59,485,393 T/C (0.45) 0.57 2.8 × 10–17 RP11-332H18.4
18 18:77,156,103–77,160,235 rs549752b 77,158,225 A/G (0.32) –0.52 3.8 × 10–12 NFATC1

CHR, chromosome; EA, effect allele; EAF, effect allele frequency; NEA, non-effect allele; rsID, reference SNP cluster ID; SNP, single nucleotide polymorphism.

a The mapped gene was annotated by ANNOVAR when performing FUMA.

b These data denote an SNP previously linked to kidney trait (e.g., estimated glomerular filtration rate, uric acid, blood pressure, or serum creatinine).

c These data indicate an SNP that has been searched for in the genome-wide association study catalog without exhibiting an association with kidney traits (rs4665985 for alcohol consumption, triglyceride, and liver fat; rs58063923 for height).


1. Jager KJ, Kovesdy C, Langham R, Rosenberg M, Jha V, Zoccali C. A single number for advocacy and communication-worldwide more than 850 million individuals have kidney diseases. Kidney Int 2019;96:1048–1050.
crossref pmid
2. Köttgen A, Glazer NL, Dehghan A, et al. Multiple loci associated with indices of renal function and chronic kidney disease. Nat Genet 2009;41:712–717.
pmid pmc
3. Köttgen A, Pattaro C, Böger CA, et al. New loci associated with kidney function and chronic kidney disease. Nat Genet 2010;42:376–384.
pmid pmc
4. Pattaro C, Teumer A, Gorski M, et al. Genetic associations at 53 loci highlight cell types and biological pathways relevant for kidney function. Nat Commun 2016;7:10023.
pmid pmc
5. Broekema RV, Bakker OB, Jonkers IH. A practical view of fine-mapping and gene prioritization in the post-genome-wide association era. Open Biol 2020;10:190221.
crossref pmid pmc pdf
6. Wuttke M, Li Y, Li M, et al. A catalog of genetic loci associated with kidney function from analyses of a million individuals. Nat Genet 2019;51:957–972.
pmid pmc
7. Hellwege JN, Velez Edwards DR, Giri A, et al. Mapping eGFR loci to the renal transcriptome and phenome in the VA Million Veteran Program. Nat Commun 2019;10:3842.
crossref pmid pmc pdf
8. Teumer A, Li Y, Ghasemi S, et al. Genome-wide association meta-analyses and fine-mapping elucidate pathways influencing albuminuria. Nat Commun 2019;10:4130.
pmid pmc
9. Kim Y, Han BG; KoGES group. Cohort profile: The Korean Genome and Epidemiology Study (KoGES) Consortium. Int J Epidemiol 2017;46:e20.
crossref pmid
10. Kim KY, Kim JO, Kim YS, et al. Genome-wide association of individual vulnerability with alcohol-associated liver disease: a Korean genome and epidemiology study. Hepatology 2022;75:391–402.
crossref pmid pdf
11. Park JS, Kim Y, Kang J. Genome-wide meta-analysis revealed several genetic loci associated with serum uric acid levels in Korean population: an analysis of Korea Biobank data. J Hum Genet 2022;67:231–237.
crossref pmid pdf
12. Gim JA, Lee S, Kim SC, Baek KW, Yoo JI. Demographic and genome wide association analyses according to muscle mass using data of the Korean Genome and Epidemiology Study. J Korean Med Sci 2022;37:e346.
crossref pmid pmc pdf
13. Lee J, Lee Y, Park B, Won S, Han JS, Heo NJ. Genome-wide association analysis identifies multiple loci associated with kidney disease-related traits in Korean populations. PLoS One 2018;13:e0194044.
crossref pmid pmc
14. Kim HR, Jin HS, Eom YB. A genome-wide association study for hypertensive kidney disease in Korean men. Genes (Basel) 2021;12:751.
crossref pmid pmc
15. Albert FW, Kruglyak L. The role of regulatory variation in complex traits and disease. Nat Rev Genet 2015;16:197–212.
crossref pmid pdf
16. Moon S, Kim YJ, Han S, et al. The Korea Biobank Array: design and identification of coding variants associated with blood biochemical traits. Sci Rep 2019;9:1382.
crossref pmid pmc pdf
17. Cho YS, Go MJ, Kim YJ, et al. A large-scale genome-wide association study of Asian populations uncovers genetic factors influencing eight quantitative traits. Nat Genet 2009;41:527–534.
crossref pmid pdf
18. Nam K, Kim J, Lee S. Genome-wide study on 72,298 individuals in Korean biobank data for 76 traits. Cell Genom 2022;2:100189.
crossref pmid pmc
19. Marees AT, de Kluiver H, Stringer S, et al. A tutorial on conducting genome-wide association studies: quality control and statistical analysis. Int J Methods Psychiatr Res 2018;27:e1608.
crossref pmid pmc pdf
20. 1000 Genomes Project Consortium, Auton A, Brooks LD, et al. A global reference for human genetic variation. Nature 2015;526:68–74.
pmid pmc
21. Browning BL, Tian X, Zhou Y, Browning SR. Fast two-stage phasing of large-scale sequence data. Am J Hum Genet 2021;108:1880–1890.
crossref pmid pmc
22. Browning BL, Zhou Y, Browning SR. A one-penny imputed genome from next-generation reference panels. Am J Hum Genet 2018;103:338–348.
crossref pmid pmc
23. Levey AS, Stevens LA, Schmid CH, et al. A new equation to estimate glomerular filtration rate. Ann Intern Med 2009;150:604–612.
crossref pmid pmc
24. Inker LA, Schmid CH, Tighiouart H, et al. Estimating glomerular filtration rate from serum creatinine and cystatin C. N Engl J Med 2012;367:20–29.
crossref pmid pmc
25. Chang CC, Chow CC, Tellier LC, Vattikuti S, Purcell SM, Lee JJ. Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience 2015;4:7.
crossref pmid pmc pdf
26. Willer CJ, Li Y, Abecasis GR. METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics 2010;26:2190–2191.
crossref pmid pmc pdf
27. Watanabe K, Taskesen E, van Bochoven A, Posthuma D. Functional mapping and annotation of genetic associations with FUMA. Nat Commun 2017;8:1826.
crossref pmid pmc pdf
28. Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res 2010;38:e164.
crossref pmid pmc
29. de Leeuw CA, Mooij JM, Heskes T, Posthuma D. MAGMA: generalized gene-set analysis of GWAS data. PLoS Comput Biol 2015;11:e1004219.
crossref pmid pmc
30. de Leeuw CA, Stringer S, Dekkers IA, Heskes T, Posthuma D. Conditional and interaction gene-set analysis reveals novel functional pathways for blood pressure. Nat Commun 2018;9:3768.
crossref pmid pmc pdf
31. GTEx Consortium. The GTEx Consortium atlas of genetic regulatory effects across human tissues. Science 2020;369:1318–1330.
pmid pmc
32. Yu G, Wang LG, Yan GR, He QY. DOSE: an R/Bioconductor package for disease ontology semantic and enrichment analysis. Bioinformatics 2015;31:608–609.
crossref pmid pdf
33. Finucane HK, Reshef YA, Anttila V, et al. Heritability enrichment of specifically expressed genes identifies disease-relevant tissues and cell types. Nat Genet 2018;50:621–629.
crossref pmid pmc pdf
34. Pers TH, Karjalainen JM, Chan Y, et al. Biological interpretation of genome-wide association studies using predicted gene functions. Nat Commun 2015;6:5890.
crossref pmid pdf
35. bulik/ldsc. ldsc wiki [Internet]. rkwalters; c2019 [cited 2023 Mar 26]. Available from:
36. Yang J, Ferreira T, Morris AP, et al. Conditional and joint multiple-SNP analysis of GWAS summary statistics identifies additional variants influencing complex traits. Nat Genet 2012;44:369–375.
crossref pmid pmc pdf
37. Stanzick KJ, Li Y, Schlosser P, et al. Discovery and prioritization of variants and genes for kidney function in >1.2 million individuals. Nat Commun 2021;12:4350.
pmid pmc
38. Piñero J, Queralt-Rosinach N, Bravo À, et al. DisGeNET: a discovery platform for the dynamical exploration of human diseases and their genes. Database (Oxford) 2015;2015:bav028.
crossref pmid pmc
39. Larsen CP, Trivin-Avillach C, Coles P, et al. LDL Receptor-related protein 2 (Megalin) as a target antigen in human kidney anti-brush border antibody disease. J Am Soc Nephrol 2018;29:644–653.
crossref pmid
40. Khalili H, Sull A, Sarin S, et al. Developmental origins for kidney disease due to Shroom3 deficiency. J Am Soc Nephrol 2016;27:2965–2973.
crossref pmid pmc
41. Gu X, Yang H, Sheng X, et al. Kidney disease genetic risk variants alter lysosomal beta-mannosidase (MANBA) expression and disease severity. Sci Transl Med 2021;13:eaaz1458.
crossref pmid pmc
42. Daga A, Majmundar AJ, Braun DA, et al. Whole exome sequencing frequently detects a monogenic cause in early onset nephrolithiasis and nephrocalcinosis. Kidney Int 2018;93:204–213.
crossref pmid
43. Pedigo CE, Ducasa GM, Leclercq F, et al. Local TNF causes NFATc1-dependent cholesterol-mediated podocyte injury. J Clin Invest 2016;126:3336–3350.
crossref pmid pmc
44. Joseph CB, Mariniello M, Yoshifuji A, et al. Meta-GWAS reveals novel genetic variants associated with urinary excretion of uromodulin. J Am Soc Nephrol 2022;33:511–529.
crossref pmid pmc
45. Sollis E, Mosaku A, Abid A, et al. The NHGRI-EBI GWAS Catalog: knowledgebase and deposition resource. Nucleic Acids Res 2023;51:D977–D985.
crossref pmid pdf
46. Sakaue S, Kanai M, Tanigawa Y, et al. A cross-population atlas of genetic associations for 220 human phenotypes. Nat Genet 2021;53:1415–1424.
crossref pmid pdf
47. Maurano MT, Humbert R, Rynes E, et al. Systematic localization of common disease-associated variation in regulatory DNA. Science 2012;337:1190–1195.
crossref pmid pmc
48. Mogil LS, Andaleon A, Badalamenti A, et al. Genetic architecture of gene expression traits across diverse populations. PLoS Genet 2018;14:e1007586.
crossref pmid pmc

Editorial Office
#301, (Miseung Bldg.) 23, Apgujenog-ro 30-gil, Gangnam-gu, Seoul 06022, Korea
Tel: +82-2-3486-8736    Fax: +82-2-3486-8737    E-mail:                

Copyright © 2024 by The Korean Society of Nephrology.

Developed in M2PI

Close layer