Skip to main content

Meta-analysis of African ancestry genome-wide association studies identified novel locus and validates multiple loci associated with kidney function

Abstract

Despite recent efforts to increase diversity in genome-wide association studies (GWASs), most loci currently associated with kidney function are still limited to European ancestry due to the underlying sample selection bias in available GWASs. We set out to identify susceptibility loci associated with estimated glomerular filtration rate (eGFRcrea) in 80027 individuals of African-ancestry from the UK Biobank (UKBB), Million Veteran Program (MVP), and Chronic Kidney Disease genetics (CKDGen) consortia.

We identified 8 lead SNPs, 7 of which were previously associated with eGFR in other populations. We identified one novel variant, rs77408001 which is an intronic variant mapped to the ELN gene. We validated three previously reported loci at GATM-SPATA5L1, SLC15A5 and AGPAT3. Fine-mapping analysis identified variants rs77121243 and rs201602445 as having a 99.9% posterior probability of being causal. Our results warrant designing bigger studies within individuals of African ancestry to gain new insights into the pathogenesis of Chronic Kidney Disease (CKD), and identify genomic variants unique to this ancestry that may influence renal function and disease.

Peer Review reports

Introduction

Chronic kidney disease (CKD) is defined as kidney damage or an estimated glomerular filtration rate (eGFR) < 60 ml/min/1.73m2 for three or more months [1]. The global prevalence of CKD is estimated to be ~ 10% in adults [2], and sub-Saharan Africa has an estimated prevalence of 13.1% [3]. Care and management of CKD impact a significant financial burden on healthcare systems worldwide, even more so in developing economies, majorly in Africa [4].

The genetic risk factors associated with CKD are ethnically driven including but not limited to telomere length, copy number variation (CNVs), and mitochondrial DNA [5]. Numerous loci related to estimated glomerular filtration rate (eGFR), which is a measure of kidney function for defining CKD, have mainly been discovered using Genome-wide association studies (GWAS) of European and Asian ancestry [6,7,8]. Owing to the systematic sample selection bias in global GWASs, many of the significant loci discovered in Eurocentric populations can hardly be replicated in African ancestry groups. This is attributed to differences in linkage disequilibrium (LD), differences in heritability, and differences in allele frequencies across ancestries [5, 9,10,11,12]. Polymorphisms in two critical genes, the plasma kallikrein gene (KLKB1), and the human homolog of the rodent renal failure gene, have been reported to be more associated with End-stage renal disease (ESRD) in the American black population than in whites [13].

In the past 15 years, efforts have been made to increase diversity in cohorts that perform GWASs of complex traits by increasing the number of previously under-represented ancestries, like the African-ancestry [14,15,16,17]. However, these efforts are still limited by various logistical challenges. Thus, low power due to the underlying sampling bias towards European ancestry individuals in most global consortia is still eminent [18, 19]. Most such studies have included the African-ancestry population as a small validation cohort in trans-ethnic analyses, mainly for comparison to a much larger European-ancestry cohort. This has not yielded much ancestry-specific discovery of loci to CKD. As a result, such study designs could have missed some of the loci unique to African-ancestry individuals within these global consortia.

Therefore, our study leveraged eGFRcrea GWAS summary data of only African ancestry individuals within three global consortia: Million Veteran Program (MVP), UK Biobank (UKBB), and Chronic Kidney Disease Genetics Consortium (CKDGen consortium), with a combined sample size of over 80 k African-ancestry participants. We aimed to uncover loci associated with CKD in these individuals, further characterize their functionality and phenotype association by performing a Phenome-wide association analysis, and determine the causality of variants using Bayesian Fine-mapping methods.

Results

Study samples

GWAS summary statistics of eGFR for 80,027 individuals of African ancestry from two global population-based consortia: UKBB [20] and CKDGen [21], and one Hospital-based cohort, MVP [22] were meta-analyzed. A summary of study sample characteristics is shown in Table 1. Further details of the cohorts contributing to this study are in the Methods section and the Supplementary material.

Table 1 Characteristics of the study cohorts

Genome-wide associated SNPs with eGFR

We analyzed ~ 22.6 M SNPs and 403 SNPs from those attained genome-wide significance (P < 5 X 10E-8). We identified lead SNPs by clumping at ± 500 kb window around the significant SNPs and identified 8 lead variants. Seven of them had previously been reported in association with eGFRcrea. The smallest and most genome-wide significant p-value from the meta-analysis was 1.326e-38 at marker 15:45592887. Figure 1 summarizes the meta-analysis results of African-ancestry eGFRcrea for the three cohorts. Our power calculation indicated an 80% power to determine genome-wide significance in our analyses.

Fig. 1
figure 1

Manhattan plot for Meta-analysis results showing -log10 (P value) versus genomic position plots for eGFRcrea. Known lead variants are plotted in red, and novel lead variants are in green

SNP rs201602445, located on chromosome 15, at base position 45592887, had the strongest association (p-value = 1.326e-38, Fig. 1, Supplementary Table s1). This SNP is located near the spermatogenesis-associated 5-like 1 (SPATA5L1) at the GATM-SPATA5L1 locus. We identified one intronic variant, rs77408001, mapping to the ELN gene (Table 2, Fig. 2), which had not previously been associated with eGFRcrea, and we considered this as a novel SNP henceforth. This SNP was replicated by proxy in the Uganda genome resource (UGR) dataset [25]. The proxied SNP was rs141339897 (chr7:73727705) (r = 1, D’ = 1 in MSL Haplotypes) as shown in Table 2. Other proxied replicated SNPs are shown in Supplementary Table s6, and the locus zoom plot for the most significant proxied SNP (rs3135654) is shown in Supplementary figure s3.

Table 2 Genome-wide significant locus (P < 5 × 10–8) associated with eGFR after meta-analysis
Fig. 2
figure 2

Regional association plot for eGFR in > 80 k African-ancestry individuals from UKBB and MVP consortium at ELN locus. The lead SNP rs77408001 is coloured purple. LD (r2) with other SNPs at the locus was calculated using an African Ancestry reference panel from 1000 Genomes. The arrows represent direction of transcription

Lead SNP rs10084572 is located near 1-acylglycerol-3-phosphate O-acyltransferase 3 (AGPAT3), found on chromosome 21, and lead SNP rs200950799 is located near Solute Carrier Family 15, Member 5 gene (SLC15A5) on chromosome 12 have been previously reported in association with eGFRcrea in non-Hispanic blacks. Further information about the genetic architecture of the lead SNPs is illustrated using regional association plots in Supplementary figure s1.

We applied genomic control correction to the meta-analysis results to account for inflation. We found no evidence of significant variation in expected allele frequencies across individuals in each cohort due to residual population structure as shown in Fig. 3. The inflation lambda was 1.07, indicating minimal bias in the test statistics.

Fig. 3
figure 3

Quantile–Quantile plot of genome-wide associations of eGFR based on serum creatinine

Bayesian fine-mapping

Bayesian fine-mapping was performed based on a ± 500 kb window on each of the lead SNPs from the meta-analysis. Two variants, rs77121243 and rs201602445, each had one credible set accounting for 99% of the posterior probability (PP) of driving the eGFR association in African-ancestry individuals. SNPs 12:17157119 and 21:45412872, each had a credible set-size of 2, accounting for 99% of the PP of being causal SNPs. Variants 7:73443012 (rs77408001), 12:17157119 (rs200950799), and 21:45412872 (rs10084572) had a credible set-size of 5, 2, and 2 respectively (Supplementary Table s2). We then proceeded to compare the 99% credible set sizes of these loci in EUR GWAS. They were smaller 99% credible set sizes in AFR compare to EUR GWAS except for 7:1286192 loci (Supplementary Table s3). This indicated that AFR GWAS enhanced the fine mapping resolution compared to the EUR GWAS.

eQTL analysis

We conducted kidney expression quantitative trait locus (eQTL) analysis in Genotype-Tissue Expression (GTEx), NephQTL [26], the Human Kidney eQTL Atlas [27], and RegulomeDB [28]. We found significant eQTLs in the GTEx database for 3/8 of the lead SNPs; rs9908131, rs10945657, and rs13230509 (Supplementary table s5) at (P < 5E-8). Variants: rs9908131 and rs10945657 had significant eQTLs in both GTEx and NephQTL resources. Variant rs13230509 was not significantly associated with any eQTLs in NephQTLs. For the three variants with significant eQTLs, we concluded that there is evidence of colocalization between eGFR association and eQTLs. All eight lead SNPs were not significant eQTLs in the Human Kidney Atlas. All the lead SNPs showing significant eQTLs were mapped to non-kidney tissues. Using RegulomeDB, we found no evidence of the regulatory impact of the lead SNPs and any of those in high LD.

SNP association with related phenotypes: PheWAS and MAGMA analysis

Of all the lead SNPs, 7 showed significant association with other phenotypes. SNP variants that showed the most significant association and their phenotypes are shown in supplementary Table s 4. Variants rs77121243 and rs13230509 were most significantly associated with metabolic and haematological traits, respectively. We further performed MAGMA tissue expression analysis, but none of the variants had significant expression in the kidney.

Discussion

We present the results from a meta-analysis of GWAS summary statistics of eGFR based on serum creatinine, in individuals of African ancestry. To the best of our knowledge, this is one of few studies to combine summary statistics of eGFRcrea from strictly African-ancestry individuals within these 3-major global GWAS consortia (UKBB, CKDGen, and MVP).

None of the lead SNPs identified in our analysis had previously been reported in European ancestry as associated with eGFRcrea. We validated seven variants as mapped to previously reported eGFR loci in individuals of African Ancestry; GATM/SPATA5L1, SLC15A5, and AGPAT3 [7, 8]. The lead SNP with the most significant association is located near the GATM/SPATA5L1 locus, previously mapped to creatinine synthesis in continental Africans elsewhere [25]. This lead SNP is different from those similarly mapped to the GATM locus in individuals of European ancestry [6, 24]. GATM is an enzyme involved in creatinine biosynthesis by encoding glycine amidinotransferase. SNPs at the GATM locus have been reported to have a role in serum creatinine levels without influencing susceptibility to kidney disease in individuals of European ancestry [7, 24]. A similar association has been recently reported in continental Africans [25]. Variant rs13230509 has previously been associated with eGFRcrea in both non-Hispanic blacks and non-Hispanic whites [24].

Locus SPATA5L1 is a known locus for eGFR [6]. This locus was previously reported to be associated with eGFR in individuals of European ancestry [8]. The mechanism of action of the SPATA5L1 encoded protein affects kidney function, and how abundant it is in kidney tissues is still not yet clear. Similarly, the SLC15A5 has previously been associated with eGFR in Europeans [24]. SLC15A5 encodes a protein that enables symporter activity. It is an integral component of the membrane and is predicted to be involved in peptide transport, protein transport, and transmembrane transport.

One lead SNP, rs77408001, has not yet been mapped to any eGFR loci and is located in the ELN gene. This gene encodes a protein, one of two elastic fibre components. The ELN gene, also known as the elastin gene, codes for a protein called tropoelastin from which Elastin is synthesized. Elastin is the major component of elastic fibres, which support the body's connective tissues [29]. A deletion mutation in this gene at 7q11.23 has been linked to Williams syndrome (WS), characterized by cardiovascular disease including elastin arteriopathy, peripheral pulmonary stenosis, supravalvar aortic stenosis, and hypertension [30].

The bayesian fine-mapping analysis identified two lead SNPs, rs77121243 and rs201602445, with a single credible set-size each showing a 99.9% posterior probability of driving the association at respective loci. The lead SNP; rs201602445, located near the more characterized locus; GATM, is significantly associated with serum creatinine in the African-ancestry population [25]. Overall, the aggregation of African Ancestry GWAS enhanced the fine mapping resolution compared to European GWAS.

Most of the lead variants were not significant expression quantitative trait loci (eQTLs) in publicly available kidney gene expression resources, and they did not co-localize to any kidney tissues. This is unsurprising considering the lack of generalizability across populations in these resources. For example, over 80% of the expression quantitative trait locus (eQTL) is from individuals of European ancestry [31]. The most significant PheWAS association was with eGFR, and two SNPs, rs10084572 and rs77408001, were associated with Height and Total cholesterol in large VLDL, respectively, but these were not significant associations.

Our study had some limitations, for example, diabetes and hypertension are common risk factors for kidney damage [32], however, we were not able to perform a stratified analysis based on underlying patient kidney disease comorbidities like diabetes and hypertension within these cohorts. This would enable us to identify which comorbidity groups contributed the most toward observed associations. Our largest contributing cohort, MVP, had previously reported diabetes-stratified GWAS of over 9 k individuals with 32 genome-wide significant variants [24]. They reported that one of the loci, SPATA5L1, identified in this study had the strongest association in the non-diabetics group. In this study, we only had access to GWAS summary statistics from these cohorts with no socio-demographics information, and thus couldn’t conduct this analysis.

Another potential limitation to our findings is the large effect size and standard error of the novel SNP rs77408001 in both the meta-analysis and the MVP cohort. We mainly attributed this to different phenotype transformations by individual studies. For example, the UKBB is Inverse-rank normal transformation, the MVP dataset has no transformation, and the CKDGen is log transformed. It is for this particular reason that we performed a p-value based meta-analysis.

However, another potential limitation is that the Imputation quality information for the SNP, rs77408001, associated with the novel loci, was not available in the summary statistics of either study. This, in addition to our inability to conduct a direct replication analysis for this SNP, makes it difficult to fully affirm the true genetic effect of this SNP on eGFR in African ancestry individuals. We recommend future similarly designed studies to critically characterize this variant in more detail.

Overall, our findings uncover a potential new locus for eGFR in African-ancestry individuals and further affirm those previously reported in other ancestries. Some of the variants reported here are unique to the African ancestry, and therefore this warrants larger ancestry-specific meta-analyses to decipher more loci associated with CKD in Africans. Furthermore, we justify here that even African-ancestry individuals outside Africa have genetic differences from continental Africans, and caution should be taken not to generalize individuals within this ancestry while conducting GWASs.

Methods

Study cohort characteristics

Million Veteran Program (MVP)

A detailed description of the MVP cohort has been described elsewhere [24]. Briefly, the MVP cohort recruited fully consented participants from 63 Veteran American (VA) medical facilities. Out of the several phenotypes recorded, 280722 participants have eGFR analysis available. Of these, 57336 are individuals of African ancestry. Most participants with eGFR analysis are non-diabetics, with just ~ 33% being diabetics. The participants are further stratified into hypertensives and non-hypertensives. About 93% of hypertensive diabetics are non-Hispanic blacks, but eGFR is generally higher in non-Hispanic blacks than non-Hispanic whites in this cohort. We included all individuals (n = 57336) of African ancestry with eGFR analyzed from the MVP cohort, and we didn’t stratify them based on any underlying clinical diagnoses. More information about this cohort can be found in the Supplementary methods.

Chronic Kidney Disease Genetics (CKDGen) consortium

This is a collaborative consortium of mainly population-based studies from different ethnicities. The consortium follows a centralized analysis plan where each participating study performs an independent analysis following a centralized analysis plan. Depending on the age of participants, the studies contributing to the CKDGen consortium generally estimate the eGFR based on serum creatinine using four-variable Modification of Diet in Renal Disease Study Equation [8] or the Schwartz equations [33]. The effect of outliers on eGFR values is harmonized by averaging at 15 and 200 ml/min/1.73 m2. Further information about the specific study within the consortium, from which participants included in our study were obtained can be found at https://doi.org/10.1038/ncomms10023 [8]. This study recruited 16474 participants of African ancestry, with information on eGFRcrea analysis. More details about the CKDGen consortium can be found in Supplementary methods.

UK BioBank (UKBB)

The UKBB is a large prospective cohort with > 50 k participants aged 40–73 years. The details of steps followed during enrolment have been previously described elsewhere [20, 23]. At baseline, each participant completed a detailed, computerized questionnaire. Biochemical tests and genotyping were conducted on archived blood specimen. For longitudinal follow-up, study data was linked to participant’s national health records. Phenotype expansion is still on-going. Estimation of glomerular filtration rate based on serum creatinine was calculated using the CKD-EPI equation [34]. From the UKBB cohort, our study included 6217 individuals of African ancestry with analyzed eGFRcrea. The supplementary methods have further details on this cohort. Ethnicity was coded as white, black, south Asian or other based on survey questionnaires. However, for purposes of eGFR calculation, ethnicity was coded as black or other.

Meta-analysis

We aggregated summary statistics by performing the Stouffer’s method implemented in METAL due to the different GWAS trait transformations that had been applied by the individual studies [35]. Two of the datasets (UKBB and MVP) had their genome coordinates in build GRCh37/hg19. The genome coordinates for the CKDGen dataset used in this analysis were reported in b36 but were converted to build GRCh37/hg19 before the meta-analysis. Briefly, this method converted the two-sided p-value to a z-statistic which reflected the direction of association based on the reference allele. Each z-score was then weighted; the squared weights were chosen to sum to 1, and each cohort-specific weight was proportional to the square root of the adequate number of individuals in that cohort. The weighted z-statistics were then summed across the three cohorts, and the summary z-score was converted to a two-sided p-value. We performed double GC correction once during the study-specific analyses and repeated the statistics from the meta-analysis. QQ-plots for the individual datasets and the meta-analysis are presented in the Supplementary data.

Fine-mapping

We used a Bayesian method to fine-map all SNPs within 500 kb of the eight lead SNPs [36]. This approach has been detailed elsewhere [36]. Briefly, we used the Z-score for the ith SNP (Zi) to calculate the Baye’s factor in favour of the association, denoted BFi, using equation (i) below:

$${BF}_{i}={e}^{\left[\frac{{Z}_{1}^{2}-\mathrm{log}(\mathrm{K})}{2}\right]},$$
(1)

where K is the number of studies.

The posterior probability (PP) for the ith SNP was calculated using equation (ii) below;

$${\pi }_{i}=\frac{{\mathrm{BF}}_{\mathrm{i}}}{{\Sigma }_{j}{\mathrm{BF}}_{\mathrm{j}}}$$
(2)

where \(\mathrm{\Sigma BFi}\) denotes the total sum of SNPs at that locus. We obtained the 99% credible sets using Bayes’ factors of the SNPs after rearranging them in descending order. We only counted those SNPs that had a posterior probability of ≥ 0.99.

Exploring functional mapping and annotation

FUMA analysis

To perform gene annotation and functional mapping for each lead SNP, we used a web-based platform called FUMA. This resource uses information from multiple biological resources to facilitate functional annotation of GWAS results, gene prioritization, and interactive visualization [37].

PheWAS analysis

We characterized the functionality of the lead SNPs using the global GWAS ATLAS resource [38]. This resource compiled a catalogue of > 4 k GWAS results across over 3000 unique traits in 28 domains. Leveraging this resource, for each lead SNP, we assessed the genetic nature of these variants in effect size distribution, minor allele frequency, and biological functionality of the traits associated with the lead variants. We also characterized the extent of pleiotropy at a trait-associated locus and associated genes and investigated trait polygenicity and SNP heritability.

MAGMA gene set analysis

Gene set analysis based on P-values was performed using MAGMA [39], a tool embroiled within FUMA. The method performs gene-based P-value analysis for protein-coding genes by mapping SNPs to genes if they are located within the same genes. The default MAGMA setting used an SNP-wise model for gene analysis and a competitive model for gene set analysis. Multiple testing was corrected using the Bonferroni correction for gene analysis and FDR for gene-set research.

eQTL analysis

We leveraged four publicly available kidney eQTL resources: GTEx [31], NephQTL [26], the Human Kidney eQTL Atlas, and RegulomeDB to perform expression quantitative trait loci analysis. For each lead SNP, we searched eQTLs associated with each lead SNP in the kidney cortex, glomerular and tubulointerstitial kidney tissue, whole kidney, glomerulus and tubules in GTEx, NephQTL, and Human kidney eQTL Atlas, respectively. We used RegulomeDB to determine the regulatory impact of each lead SNP.

Power analysis

The power to achieve genome-wide significance at p < 5 × 10–8 was estimated with the software MetaGAP [40]. The link to the tool can be found at http://www.devlaming.eu/metagap.html. For the genetic association analyses, a meta-analysis of ~ 80,000 individuals will achieve 80% power to detect association of a common genetic variant that explains 0.05% of the variance in eGFR at a stringent genome-wide significance level (p < 5 × 10–8).

Availability of data and materials

The MVP dataset was accessed upon request using the accession number dbgap proposal #30287. The UKBB eGFRcrea summary statistics dataset was accessed from the UK biobank using the phenocode 30700. The CKDGen dataset was accessed on the CKDGen consortium (https://ckdgen.imbi.uni-freiburg.de/).

All other materials and related data in this study are available upon request directed to segun.fatumo@lshtm.ac.uk.

References

  1. Eknoyan G, Lameire N, Eckardt K, Kasiske B, Wheeler D, Levin A, Stevens PE, Bilous RW, Lamb EJ, Coresh JJ. KDIGO 2012 clinical practice guideline for the evaluation and management of chronic kidney disease. Kidney int. 2013;3(1):5–14.

  2. Eckardt K-U, et al. Evolving importance of kidney disease: from subspecialty to global health burden. Lancet. 2013;382(9887):158–69.

    Article  PubMed  Google Scholar 

  3. Ajayi SO, et al. Prevalence of chronic kidney disease as a marker of hypertension target organ damage in Africa: a systematic review and meta-analysis. Int J Hypertens. 2021;2021:7243523.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Xie Y, et al. Analysis of the Global Burden of Disease study highlights the global, regional, and national trends of chronic kidney disease epidemiology from 1990 to 2016. Kidney Int. 2018;94(3):567–81.

    Article  PubMed  Google Scholar 

  5. Canadas-Garre M, et al. Genetic susceptibility to chronic kidney disease - some more pieces for the heritability puzzle. Front Genet. 2019;10:453.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Wuttke M, et al. A catalog of genetic loci associated with kidney function from analyses of a million individuals. Nat Genet. 2019;51(6):957–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Kottgen A, et al. Multiple loci associated with indices of renal function and chronic kidney disease. Nat Genet. 2009;41(6):712–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Pattaro C, et al. Genetic associations at 53 loci highlight cell types and biological pathways relevant for kidney function. Nat Commun. 2016;7:10023.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Fatumo S. The opportunity in African genome resource for precision medicine. EBioMedicine. 2020;54: 102721.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Deepti Gurdasani TC, Tekola-Ayele Fasil, Pagani Luca, Tachmazidou Ioanna, et al. The African Genome Variation Project shapes medical genetics in Africa. Nature. 2015;517:328.

    Google Scholar 

  11. Gurdasani D, et al. Uganda Genome Resource Enables Insights into Population History and Genomic Discovery in Africa. Cell. 2019;179(4):984-1002 e36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Udler MS, et al. Effect of Genetic African Ancestry on eGFR and Kidney Disease. J Am Soc Nephrol. 2015;26(7):1682–92.

    Article  CAS  PubMed  Google Scholar 

  13. Hsu CY, et al. Racial differences in the progression from chronic renal insufficiency to end-stage renal disease in the United States. J Am Soc Nephrol. 2003;14(11):2902–7.

    Article  PubMed  Google Scholar 

  14. Barroso I. The importance of increasing population diversity in genetic studies of type 2 diabetes and related glycaemic traits. Diabetologia. 2021;64(12):2653–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Kelleher, J., et al., 2018.

  16. Mulder N, et al. H3Africa: current perspectives. Pharmgenomics Pers Med. 2018;11:59–66.

    PubMed  PubMed Central  Google Scholar 

  17. Ramirez AH, et al. The All of Us Research Program: Data quality, utility, and diversity. Patterns (N Y). 2022;3(8):100570.

    Article  PubMed  Google Scholar 

  18. Improving equity in human genomics research. Commun Biol, 2022. 5(1): 281.

  19. Peterson RE, et al. Genome-wide association studies in ancestrally diverse populations: opportunities, methods, pitfalls, and recommendations. Cell. 2019;179(3):589–603.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Allen N, et al. UK Biobank: current status and what it means for epidemiology. Health Policy Technol. 2012;1(3):123–6.

    Article  Google Scholar 

  21. Kottgen A, Pattaro C. The CKDGen Consortium: ten years of insights into the genetic basis of kidney function. Kidney Int. 2020;97(2):236–42.

    Article  PubMed  Google Scholar 

  22. Gaziano JM, et al. Million Veteran Program: a mega-biobank to study genetic influences on health and disease. J Clin Epidemiol. 2016;70:214–23.

    Article  PubMed  Google Scholar 

  23. Cathie Sudlow, J.G., Naomi Allen, Valerie Beral, Paul Burton, JohnDanesh,PaulDowney,PaulElliott, Jane Green, Martin Landray, Bette Liu, Paul Matthews, Giok Ong,Jill Pell, Alan Silman, Alan Young, Tim Sprosen, TimPeakman,RoryCollins, UKBiobank: An Open Access Resource for Identifying the Causes of a Wide Range of Complex Diseases of Middle and Old Age. 2015.

  24. Hellwege JN, et al. Mapping eGFR loci to the renal transcriptome and phenome in the VA Million Veteran Program. Nat Commun. 2019;10(1):3842.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Fatumo S, et al. Discovery and fine-mapping of kidney function loci in first genome-wide association study in Africans. Hum Mol Genet. 2021;30(16):1559–68.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Gillies CE, et al. An eQTL Landscape of Kidney Tissue in Human Nephrotic Syndrome. Am J Hum Genet. 2018;103(2):232–44.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Sheng X, et al. Mapping the genetic architecture of human traits to cell types in the kidney identifies mechanisms of disease and potential treatments. Nat Genet. 2021;53(9):1322–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Dong, S., et al., 2022.

  29. Brooke B. New insights into elastin and vascular disease. Trends Cardiovasc Med. 2003;13(5):176–81.

    Article  CAS  PubMed  Google Scholar 

  30. Schubert C. The genomic basis of the Williams-Beuren syndrome. Cell Mol Life Sci. 2009;66(7):1178–97.

    Article  CAS  PubMed  Google Scholar 

  31. Consortium, G.T. The Genotype-Tissue Expression (GTEx) project. Nat Genet. 2013;45(6):580–5.

    Article  Google Scholar 

  32. Kazancioglu R. Risk factors for chronic kidney disease: an update. Kidney Int Suppl (2011). 2013;3(4):368–71.

    Article  PubMed  Google Scholar 

  33. Schwartz GJ, et al. Improved equations estimating GFR in children with chronic kidney disease using an immunonephelometric determination of cystatin C. Kidney Int. 2012;82(4):445–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Andrew S. Levey, L.A.S., Christopher H. Schmid, Yaping (Lucy) Zhang, Alejandro F. Castro, Harold I. Feldman, John W. Kusek, Paul Eggers, Frederick Van Lente, Tom Greene, and Josef Coresh, A New Equation to Estimate Glomerular Filtration Rate. 2009.

  35. Willer CJ, Li Y, Abecasis GR. METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics. 2010;26(17):2190–1.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Wellcome Trust Case Control, C, et al. Bayesian refinement of association signals for 14 loci in 3 common diseases. Nat Genet. 2012;44(12):1294–301.

    Article  Google Scholar 

  37. Watanabe K, et al. Functional mapping and annotation of genetic associations with FUMA. Nat Commun. 2017;8(1):1826.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Watanabe K, et al. A global overview of pleiotropy and genetic architecture in complex traits. Nat Genet. 2019;51(9):1339–48.

    Article  CAS  PubMed  Google Scholar 

  39. de Leeuw CA, et al. MAGMA: generalized gene-set analysis of GWAS data. PLoS Comput Biol. 2015;11(4):e1004219.

    Article  PubMed  PubMed Central  Google Scholar 

  40. de Vlaming R, et al. Meta-GWAS Accuracy and Power (MetaGAP) calculator shows that hiding heritability is partially due to imperfect genetic correlations across studies. PLoS Genet. 2017;13(1):e1006495.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors thank Million Veteran Program (MVP) staff, researchers, and volunteers, who have contributed to MVP, and especially participants who previously served their country in the military and now generously agreed to enroll in the study. (See https://www.research.va.gov/mvp/ for more details). The citation for MVP is Gaziano, J.M. et al. Million Veteran Program: A mega-biobank to study genetic influences on health and disease. J Clin Epidemiol 70, 214-23 (2016). This research is based on data from the Million Veteran Program, Office of Research and Development, Veterans Health Administration, and was supported by the Veterans Administration (VA) Cooperative Studies Program (CSP) award #G002. Data were accessed through approved dbGaP proposal #30287 entitled, “Genomic determinant of Complex Diseases in African ancestry individuals”.

Funding

This work was supported by the Wellcome Trust [grant number: 220740/Z/20/Z] awarded to Segun Fatumo. This work was supported by the UK Medical Research Council (MRC) and the UK Department for International Development (DFID) under the MRC/DFID Concordat agreement, through core funding to the MRC/UVRI and LSHTM Uganda Research Unit. TC is an international training fellow supported by the Wellcome Trust grant (214205/Z/18/Z).

Author information

Authors and Affiliations

Authors

Contributions

SF conceptualised the study. CK performed the main analyses with support of OS, TM, TC and SF. CK wrote the first draft of the manuscript. SF and TC supervised the project. All authors read and approved the final version of the manuscript.

Corresponding author

Correspondence to Segun Fatumo.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

 Supplementary methods. Table S1. Association results for meta-analysis lead SNPs. Table S2. Annotation-informed fine-mapping using the AFR meta-analyzed dataset. Table S3. Annotation-informed fine-mapping using the AFR and EUR datasets. Table S4. PheWAS association analysis. Supplementary Table s5. Colocalization with gene expression from GTEx and NephQTL. Table S6. Proxy SNPs with some evidence of association with lead SNP rs77408001 (+/- 500kb) within the CKDGen dataset. Figure S1. Regional association plots showing Genetic architecture of the genome-wide significant susceptibility variants for CKD (a)-(h). The most significant SNP in each region is plotted in blue. LD based on the 1000G sample is color-coded red (r2 to top SNP 0.8–1.0), orange (0.5–0.8), yellow (0.2–0.5) and blue (<0.2). Figure S2. The Phenome-wide association (PheWAS) plot shows the significant (p ≤ 0.05) associations of 7:73443012:C:A / rs77408001 for all available traits, generated by bottom-line integrative analysis across all datasets in the Portal. The triangle data points indicate direction of effect. Supplementary figure s3. Regional association plot for proxy SNP rs3135654 in 16474 African-ancestry individuals from the CKDGen consortium dataset at the ELN locus.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Kintu, C., Soremekun, O., Machipisa, T. et al. Meta-analysis of African ancestry genome-wide association studies identified novel locus and validates multiple loci associated with kidney function. BMC Genomics 24, 496 (2023). https://doi.org/10.1186/s12864-023-09601-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-023-09601-0

Keywords