Skip to main content
Download PDF

Genetic variants affecting mitochondrial function provide further insights for kidney disease

BMC Genomics volume 25, Article number: 576 (2024) Cite this article

Abstract

Background

Chronic kidney disease (CKD) is a complex disorder that has become a high prevalence global health problem, with diabetes being its predominant pathophysiologic driver. Autosomal genetic variation only explains some of the predisposition to kidney disease. Variations in the mitochondrial genome (mtDNA) and nuclear-encoded mitochondrial genes (NEMG) are implicated in susceptibility to kidney disease and CKD progression, but they have not been thoroughly explored. Our aim was to investigate the association of variation in both mtDNA and NEMG with CKD (and related traits), with a particular focus on diabetes.

Methods

We used the UK Biobank (UKB) and UK-ROI, an independent collection of individuals with type 1 diabetes mellitus (T1DM) patients.

Results

Fourteen mitochondrial variants were associated with estimated glomerular filtration rate (eGFR) in UKB. Mitochondrial variants and haplogroups U, H and J were associated with eGFR and serum variables. Mitochondrial haplogroup H was associated with all the serum variables regardless of the presence of diabetes. Mitochondrial haplogroup X was associated with end-stage kidney disease (ESKD) in UKB. We confirmed the influence of several known NEMG on kidney disease and function and found novel associations for SLC39A13, CFL1, ACP2 or ATP5G1 with serum variables and kidney damage, and for SLC4A1, NUP210 and MYH14 with ESKD. The G allele of TBC1D32-rs113987180 was associated with higher risk of ESKD in patients with diabetes (OR:9.879; CI95%:4.440–21.980; P = 2.0E-08). In UK-ROI, AGXT2-rs71615838 and SURF1-rs183853102 were associated with diabetic nephropathies, and TFB1M-rs869120 with eGFR.

Conclusions

We identified novel variants both in mtDNA and NEMG which may explain some of the missing heritability for CKD and kidney phenotypes. We confirmed the role of MT-ND5 and mitochondrial haplogroup H on renal disease (serum variables), and identified the MT-ND5-rs41535848G variant, along with mitochondrial haplogroup X, associated with higher risk of ESKD. Despite most of the associations were independent of diabetes, we also showed potential roles for NEMG in T1DM.

Peer Review reports

Introduction

Chronic kidney disease (CKD) is increasing in prevalence, especially in older populations, and is a major global health problem [1,2,3,4]. CKD is predicted to become the fifth leading cause of death by 2040 [5, 6]. The most severe form of CKD is known as CKD stage 5 or end-stage kidney disease (ESKD), which can be managed by kidney replacement therapy (KRT), such as chronic dialysis and/or kidney transplantation. KRT is expensive, complex, and is often unavailable for persons with ESKD in low- and middle-income countries [7].

CKD is a complex heterogeneous disease whose causality is driven both by genetic and environmental factors. Diabetes and hypertension are major aetiologies contributing to CKD burden [3, 8]. Although CKD heritability can be up to 75% [9,10,11], molecular markers, identified mainly by meta-analyses of genome-wide association studies (GWAS), do not account for all the inherited susceptibility to CKD [12, 13]. It is therefore plausible that other genetic factors, beyond single nucleotide changes identified from nuclear GWAS, may contribute to CKD [14, 15].

Mitochondria are double-membrane-bound organelles responsible for generating the necessary energy for cellular metabolism [16,17,18]. Mitochondria contain several copies of their own genome, a circular DNA molecule of ≈16.6 kb in humans including a total of 37 genes: 13 code for the subunits of respiratory complexes I, III, IV and V [19]; 22 code for transfer RNAs (tRNAs) for the 20 standard amino acids, plus an extra gene for leucine and serine [16, 20, 21], and two for ribosomal RNAs (rRNAs) [22]. The replication origin(s) and promoters for mitochondrial DNA (mtDNA) are contained in an additional displacement loop (D-loop). The functions of mitochondria are also regulated by nuclear genes encoding proteins related to mtDNA transcription, replication, cell apoptosis and mitophagy, nucleotide biosynthesis, metabolism, and iron and calcium homeostasis [23, 24].

Mitochondrial dysfunction in kidney tissue can severely impact kidney health and has previously been implicated in CKD development [25,26,27,28,29,30,31,32,33,34,35]. Maintenance of mitochondrial integrity has been highlighted in limiting the progression of acute kidney injury to CKD [36,37,38]. Lower mtDNA copy number in peripheral blood has been associated with higher risk of diabetes and microalbuminuria, two important risk factors for CKD progression, and with a higher incidence of CKD [39, 40]. Genetic variants in the mitochondrial D-loop have been proposed as predictors of kidney survival in CKD patients, helping to identify CKD patient subgroups at higher risk of poor outcomes [41,42,43,44]. Variants in MT-ND5 have been implicated in adult-onset kidney disease [45]. Genetic variation in mtDNA may be inherited or acquired; human mtDNA is particularly susceptible to acquiring somatic mutations due to its close location to the generation of mutagenic reactive oxygen species through oxidative phosphorylation [46, 47], and the limited nucleic acid repair mechanisms. Higher mutational rates in mtDNA have been reported in tumours, which may correspond to the increased level of reactive oxidative species in renal parenchymal cells in ESKD [48]. Despite the mitochondrial genome being minimally investigated in relation to CKD, some mitochondrial proteins, encoded by nuclear-encoded mitochondrial genes (NEMG), and specific mtDNA variations in MT-HV2, MT-HV3, MT-ND5 and MT-RNR2 have been associated with kidney disease and/or well-established serum clinical biomarkers of CKD, such as serum creatinine (SCr) levels, and estimated glomerular filtration rate (eGFR) [14, 25, 40, 49,50,51]. Among NEMG, NAT8, CPS1, GATM, SLC22A2, WDR72 and AGXT2 are known susceptibility loci for CKD, progression and/or kidney function [52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71].

Aim

The aim of this study was to investigate the role of genetic variants influencing mitochondrial function on CKD (and related traits). The investigations explore genetic variants in both mtDNA and NEMG. To gain insights into the impact of diabetes on the association of mtDNA and NEMG variants with CKD, we investigated renal phenotypes in a large population cohort (UK Biobank) stratified by diabetes and a collection of individuals with type 1 diabetes mellitus (T1DM) and known kidney function who were explicitly recruited to explore molecular associations with diabetic kidney disease (UK-ROI).

Methodology

Study design and populations

This cross-sectional study included only participants of European ethnicity with body-mass index (BMI) values within 18.5 and 40 kg/m2, corresponding to nutritional status between normal weight and obesity class II, according to the World Health Organisation [72]. Two populations were included in the study (Fig. 1): the UK Biobank (UKB), a population cohort, and the UK-ROI collection of individuals diagnosed with T1DM and known kidney status, as detailed in Supplementary Methods, within the Supplementary Material (Supplementary Methods, Populations). To ensure there was no overlap between these populations, only UK-ROI participants from the Republic of Ireland and Northern Ireland were included.

Fig. 1
figure 1

Design of the study, showing final number of participants after quality control, mitochondrial and NEMG variants included, along with the kidney phenotypes studied in each cohort. Abbreviations: eGFR: estimated globular filtration rate. NEMG: nuclear-encoded mitochondrial genes. T1DM: type 1 diabetes mellitus

Full size image

Phenotypic variables

Outcome variables included CKD, ESKD, kidney damage (any pathology indicating kidney injury), serum variables, diabetic kidney disease (DKD), and kidney replacement therapy (KRT). Serum variables included serum creatinine (SCr), serum cystatin C (SCysC), eGFR based on SCr (eGFRCrea), SCysC (eGFRCysC) and SCr/SCysC combined equation (eGFRCreaCysC)). Detailed definitions are provided within the ‘Supplementary Methods’ section of the Supplementary Material. Not every variable could be considered for each population (Fig. 1 and Supplementary Table 1).

Analyses were adjusted for age, sex, genotyping batch, smoking habit (defined as ‘yes’ if the patient had ever smoked), hypertension and diabetes (in the analysis of the Overall Cohort for CKD). Data Fields from the UKB used to calculate variables can be found in Supplementary Table 2, along with detailed information on calculations in Supplementary Table 3 and Supplementary Table 4.

Genotyping and quality control

UK Biobank

The Applied Biosystems™ UK Biobank Axiom™ and UK BiLEVE Axiom™ Affymetrix Arrays were used for genotyping by the UK Biobank. Genotypes were imputed by the UKB using a combination of the Haplotype Reference Consortium and merged UK10K and 1000 Genomes phase 3 reference panels [73]. PLINK 1.90 beta and PLINK 2.00 alpha were used to perform quality control (QC) and association analysis [74, 75]. Before QC, the study was comprised of 488,377 participants, 711,188 SNPs in NEMG and 265 mitochondrial directly genotyped variants. Individuals with high missingness rate or call rate lower than 95% were removed. Related individuals (identity by kinship coefficient > 0.0884) and principal component analysis (PCA) outliers, as calculated by the UK Biobank, were also removed [73]. Variants with minor allele frequency (MAF) < 1%, minor allele count (MAC) < 20 or SNP call rate < 95% were removed from the analysis. Autosomal SNPs not fulfilling Hardy–Weinberg equilibrium (HWE) (p < 1e-20) or imputation score under 0.3 were also excluded. After QC, 374,760 participants (Overall Cohort), 371,542 variants in 2,527 NEMG, along with 192 variants in 32 mitochondrial genes for the combined arrays and 93 variants in 28 genes for the BiLEVE array remained.

UK-ROI

Genotyping of the UK-ROI data collection was performed as described by Sandholm and colleagues [76]. Briefly, DNA samples collected from 1,804 white individuals with T1DM were genotyped with the HumanOmni1-Quad array and imputed to 1000 Genomes Phase 3 reference panel. Before QC, the study was comprised of 1,804 participants, 2,484,108 autosomal variants and 212 mitochondrial DNA directly genotyped variants. Variants with MAF < 1%, MAC < 20 or SNP call rate < 90% were removed from the analysis. Autosomal SNPs not fulfilling HWE (p < 1e-20) or imputation score under 0.3 were also excluded. After QC, 1,468 individuals, 423,906 variants in 2,527 NEMG and 60 variants in 23 mitochondrial genes passed filters.

Mitochondrial haplogroups

Mitochondrial haplogroups were estimated using HaploGrep2 [77], based on PhyloTree17 [78]. Only the major European haplogroups H, V, HV, J, T, U, K, Z, W, X, I and N were considered, grouping the remaining options in the “Other” category.

Selection of nuclear-encoded mitochondrial genes

A total of 2,527 unique autosomal genes coding for 22,713 transcripts were investigated. The selection process produced 2,448 unique genes returned from database searches with a further 180 genes identified from literature searches for genes influencing mitochondrial function [50]. Briefly, several online databases and literature resources were searched for NEMGs: Mitoproteome [79,80,81,82], MitoMiner [83], MitoMap [84], Ensembl [85] and UniProt [86]. Genes extracted from individual sources were reviewed and duplicates were excluded. Gene names were then screened to ensure there was no duplication between the database searches and literature searches. Genes were annotated with their official HUGO Gene Nomenclature Committee (HGNC) gene symbol [87] using Ensembl BioMart release 67 (May 2012) based on the February 2009 Homo sapiens high coverage assembly GRCh37 from the Genome Reference Consortium [85]. Any genes not found in the BioMart [85] search were manually annotated according to their official HGNC gene symbol [87]. The list of genes was then checked again for duplicates based on HGNC symbols, known pseudonyms and gene positions. Only genes found in autosomes were included in the analysis. Any genes on sex chromosomes, non-human genes, or bacterial artificial chromosomes were excluded from the final list of genes encoding proteins required for mitochondrial function.

eQTL

Tissue data from kidney-cortex used for the analyses described in this manuscript were downloaded from the Portal of The Genotype-Tissue Expression (GTEx) project (GTEx Analysis V8 release; dbGaP accession number phs000424.v8.p2) on 31/01/2022.

Statistical analysis

Descriptive analyses were performed using R [88]. Qualitative variables were expressed as percentage (%) of their total. Non-normally distributed variables were expressed as the median and the interquartile range (Q3-Q1). Normality was assessed with Kolmogorov–Smirnov test.

Association analysis

Association analysis for individual variants including sex, age, genotyping batch, smoking habit, hypertension, diabetes (Overall Cohort) and 10 PCAs as covariates was performed with PLINK 2.00 alpha, using the ‘–glm’ flag [74]. For binary phenotypes (CKD, ESKD and kidney damage) –glm fits a logistic or Firth regression model [74]. For quantitative phenotypes, –glm fits the linear model [74]. Quantitative outcome variables were natural logarithmic transformed and analysed using the additional ‘–pheno-quantile-normalize’ flag, to force quantitative phenotypes to a N(0, 1) distribution, preserving only the original rank orders [74]. The –glm flag performs a multicollinearity check before each regression, which skips and reports 'NA' results when it fails. Fixed effects meta-analysis was performed for associated variants in the UK-ROI collection showing consistent direction of effect for the same phenotype with UKB using METAL [89]. Between-study heterogeneity was assessed with the I2 statistic [90] and random effects meta-analysis was performed if showing a significant heterogeneity (p > 0.05), calculating tau-square and the random effects parameters.

Multiple comparisons correction

To correct for multiple testing, a Bonferroni correction for the number of independent variants (estimated using a pruning procedure of our data; r2 < 0.2, window size 50 bp, offset 5 bp) after QC was used [91]. The pruning estimated 47 independent variants for the mitochondrial chromosome for the combined arrays of the UK Biobank (35 when only the BiLEVE array was considered), yielding a threshold of 1E-03, and 57,457 variants for NEMG, yielding a threshold of 9E-07. Our estimation for the mitochondrial chromosome is similar to that calculated by Kraja et al., obtained by permutation analysis, which concluded that 49 variants represented the number of independent genetic effects in the mitochondrial chromosome [92]. For the UK-ROI collection, 60 variants in 23 genes passed QC for the mitochondrial chromosome and 27 were estimated after pruning (0.05/27 = 1.9E-03), with 40,236 variants for NEMG (0.05/40,236 = 1.2E-06).

Clumping & annotation

Independent loci were identified using PLINK 1.90 beta clumping procedure (–clump-p1 5e-05 –clump-r2 0.1 –clump-kb 500) [74]. A physical distance threshold for clumping of 1 kb was used for the mitochondrial chromosome. The independent loci were annotated using the variant annotation tool TabAnno [93] complemented with NCBI dbSNP database [94] and SNPnexus [95,96,97,98,99].

Mitochondrial haplogroups

Association analysis for mitochondrial haplogroups was performed using logistic regression in R version 3.6.0 (2019–04-26) [88], including as covariates duration of T1DM, sex, hypertension (UK-ROI); age, sex, genotyping batch, diabetes (Overall Cohort), smoking habit and hypertension (UK Biobank). Each haplogroup was analysed separately using all the other haplogroups as reference, after constructing dummy variables taking the values of 0 and 1, with the R package "fastDummies" [100]. Principal components were not used as covariates to account for ancestry because of their potential correlation with haplogroups. The Bonferroni correction was applied to account for multiple comparisons, adjusting the p-value threshold, dividing by the number of haplogroups in each dataset (0.05/number of haplogroups).

Results

A descriptive analysis of the populations included in the study can be found in Table 1, showing the different strata for the UKB and the diabetic kidney disease (DKD) for the UK-ROI. UKB participants in the subgroup with diabetes were older than those without diabetes and there was a higher percentage of females. In the UK-ROI collection, those with DKD were older and more likely to be males.

Table 1 Descriptive analysis of the populations included in the study. For the UK Biobank, a description of the different variables for the overall cohort and the two subgroups after stratification are shown (with and without diabetes)
Full size table

Mitochondrial Variants

UK Biobank

A higher risk of ESKD was associated with the G allele of the MT-ND5-rs41535848 variant in the Non-Diabetic Cohort of the UKB (OR: 4.971; CI95%: 1.990–12.415; P = 6.0E-04). No other mitochondrial variants were associated with CKD, ESKD or kidney damage in any cohort of the UKB after multiple comparisons correction (p < 1E-03).

Fourteen and fifteen mitochondrial variants were associated with eGFRCrea (Supplementary Tables 5 and 10) and eGFRCreaCysC (Supplementary Tables 7 and 12) in the Overall and Non-Diabetic Cohorts, respectively. Ten mitochondrial variants were associated with eGFRCysC in both cohorts (Supplementary Tables 6 and 11).

Fourteen mitochondrial variants in ten genes were consistently associated with eGFR (any equation) both in the Overall and Non-Diabetic Cohorts; additionally, MT-CO1-G6734A and MT-CO3-rs41482146 were associated with eGFRCrea and MT-ND4-rs2853493 with eGFRCreaCysC, in the Non-Diabetic Cohort (Supplementary Tables 5–7 and 10–12). None of the mitochondrial variants showed associations in the Diabetic Cohort. Fifteen mitochondrial variants were associated with SCr (Supplementary Tables 8 and 13) and ten with SCysC (Supplementary Tables 9 and 14) in the Overall and Non-Diabetic Cohorts.

UK-ROI

No mitochondrial variants were associated with CKD or serum variables after multiple comparisons correction (p < 1.9E-03); significant variants from the UKB array rs2854131, rs28359178, rs2853506 and rs2857290 showed consistent directions of effect for eGFRCrea and SCr.

Mitochondrial haplogroups

UK Biobank

Six mitochondrial haplogroups were present in the UKB cohort with a frequency over 3% (H, U, J, T, K and I) (Supplementary Table 15). The mitochondrial haplogroups X, H, U, J and I were associated with a number of the phenotypes investigated at the corrected p-value of 0.05/13 = 0.003846 (Supplementary Table 16). Figure 2 shows the main results for mitochondrial variants and haplogroups associated with some phenotype in the UKB.

Fig. 2
figure 2

Main results for mitochondrial variants and haplogroups associated with some phenotype in the UK Biobank (Expressed as odds ratio (95% confidence interval) for qualitative variables and betas ± standard error for quantitative variables; Colour code: green increases and red decreases the beta estimate; NA: Not Associated). eGFRCrea: eGFR based on serum creatinine. eGFRCreaCysC: eGFR based on both serum creatinine and cystatin C. eGFRCysC: eGFR based on serum cystatin C. SCr: serum creatinine. SCysC: serum cystatin C. Mitochondrial variants are described as Gene | Base Position (GRCh37) | rs identifier | Effect Allele | Alternative Allele

Full size image

The mitochondrial haplogroup X, defined by the mutations T6221C, C6371T, A13966G, T14470C, T16189C and C16278T, was associated with higher risk of ESKD in both the Overall and Non-Diabetic cohorts of the UKB, but not in the Diabetic Cohort (Supplementary Table 16), as consistently shown by the association of the A13966G mutation (rs41535848) with ESKD in the Non-Diabetic Cohort (OR: 4.971;CI95%: 1.990–12.415; P = 6.0E-04).

The mitochondrial haplogroup H, defined by the mutations G2706A (rs2854128) and T7028C, was associated with all the serum variables across all cohorts (Supplementary Table 16), as consistently shown by those two mutations in the Overall and Non-Diabetic cohorts (Supplementary Tables 6–14).

The mitochondrial haplogroup U, defined by the A11467G, A12308G and G12372A mutations, was associated with all the serum variables in both the Overall and Non-Diabetic cohorts (Supplementary Table 16), as seen with the A11467G (rs2853493) for eGFRCreaCysC in the Non-Diabetic cohort (Supplementary Table 12) and G12372A (rs2853499) for all serum variables in the Overall and Non-Diabetic cohorts (Supplementary Tables 6–11;13;14).

The mitochondrial haplogroup J, defined by the C295T, T489C, A10398G, A12612G, G13708A and C16069T mutations, was associated with SCr, SCysC and eGFRCreaCysC in the Non-Diabetic Cohort, but only with the latter in the Overall Cohort (Supplementary Table 16). Consistent associations were found for the individual mutations C295T (rs41528348) and G13708A (rs28359178) (Supplementary Tables 6–14).

The mitochondrial haplogroup I, defined by the mutations T10034C and G16129A, was associated with SCr in both the Overall and Non-Diabetic cohorts (Supplementary Table 16).

The category “Other”, encompassing all the mitochondrial haplogroups not common in Europeans, was associated with SCysC and its derived eGFR equation in both the Overall and Non-Diabetic cohorts (Supplementary Table 16).

UK-ROI

The Z mitochondrial haplogroup was not present in the UK-ROI collection, leaving 12 categories to compare (corrected p-value of 0.05/12 = 0.004). The mitochondrial haplogroup J, defined by the C295T, T489C, A10398G, A12612G, G13708A and C16069T mutations, showed a trend to higher risk of kidney replacement therapy (OR: 1.957; CI95%: 1.188–3.178; P = 0.00727), lower eGFR (BETA: -0.115; CI95%: -0.211-(-0.019); P = 0.0186) and higher SCr levels (BETA: 0.098; CI95%: 0.015–0.181; P = 0.0201) although not significant after multiple comparisons correction.

NEMG

UK Biobank

Up to 360 NEMG variants were associated with any serum variable in the UKB dataset (Supplementary Table 17). Among them, only the intronic variant in the NADH:ubiquinone oxidoreductase subunit B10, NDUFB10-rs338788C, associated with higher eGFRCrea in the Non-Diabetic Cohort, is a known expression quantitative trait locus (eQTL) in kidney cortex, in particular, increasing the expression of the ribosomal protein L26 (RPL26) pseudogene, according to GTEx (Supplementary Table 18). Three SNPs (SLC39A13-rs2293576, LMNA-rs4641 and TRMT1-rs35601737), associated with different kidney phenotypes, are known splicing quantitative trait loci (sQTL) in kidney cortex (Supplementary Table 18). The SLC39A13-rs10742802 variant, in strong LD with SLC39A13-rs2293576 (R2 = 0.7632) was also associated with Kidney Damage and the serum variables eGFRCrea, eGFRCreaCysC, SCr and SCysC (Supplementary Table 18).

Overall cohort

Five index SNPs were associated with higher risk of CKD with a p-value < 9E-07 in 373,164 participants of the Overall Cohort of the UKB (NAT8-rs13538A, CPS1-rs1047891A, SLC22A2-rs3127573G, GATM-rs58764877C, WDR72-rs72747347C) (Supplementary Table 19). NAT8-rs13538 and CPS1-rs1047891 were also associated with all the serum variables in the Overall and Non-Diabetic Cohort, along with GATM-rs58764877 for some of the variables. They were also associated with at least one of the serum variables in the Non-Diabetic Cohort. RAB24-rs80237806, associated with all the serum variables in the Overall and Non-Diabetic Cohorts, had a p-value for CKD of 2.3E-06.

Two clusters were associated with higher risk of ESKD with a p-value < 9E-07 (NUP210-rs144856263T and SLC4A1-rs116844389A) for the index SNP in the Overall Cohort of the UKB (Supplementary Table 19).

Twenty-one index SNPs were associated with kidney damage with a p-value < 9E-07 in 352,722 participants of the Overall Cohort of the UKB (Supplementary Table 19), three in common with serum variables.

CFL1-rs117624356, NOS3-rs3918226 and SLC39A13-rs10742802 were associated with all the serum variables in the Overall Cohort; ACP2-rs75393320, ATP5G1-rs1800632 with SCysC-related variables and SQOR-rs629024 with SCr-related variables.

Non-diabetic cohort

Four index SNPs were associated with CKD with a p-value < 9E-07 in 340,185 participants of the Non-Diabetic Cohort of the UKB (Supplementary Table 20), three in common with the Overall Cohort (NAT8-rs13538, GATM-rs58764877 and WDR72-rs72747347). NAT8-rs13538 was also associated with all the serum variables in the Non-Diabetic Cohort. GATM-rs58764877 was associated with some of the serum variables in the Non-Diabetic Cohort.

Two index SNPs (NUP210-rs144856263T and MYH14-rs148695576T) were associated with higher risk of ESKD with a p-value < 9E-07 in 340,185 participants of the Non-Diabetic Cohort of the UKB (Supplementary Table 20), the first one in common with the Overall Cohort, none in common with serum variables.

Eleven index SNPs were associated with kidney damage with a p-value < 9E-07 in 320,718 participants of the Non-Diabetic Cohort of the UKB (Supplementary Table 20), all in common with the Overall Cohort. Among them, SLC39A13-rs10742802 was associated with all serum variables and NOS3-rs3918226 and SQOR-rs629024 with the SCr-related variables in the Non-Diabetic Cohort.

Diabetic cohort

No variants were found to be associated with CKD or kidney damage with a p-value < 9E-07 in 32,979 participants of the Diabetic Cohort of the UK Biobank. Patients with diabetes carrying the G allele of TBC1D32-rs113987180 variant (MAF 3%) were at a much higher risk of developing ESKD than patients without diabetes (OR: 9.879; SE: 0.408; P = 2.0E-08), whereas the effect in the other UKB cohorts was weaker and non-statistically significant, indicating a potential interaction with the presence of diabetes (Non-Diabetic: OR: 1.226; SE:0.456; P = 6.6E-01 and Overall: OR: 2.653; SE: 0.288; P = 7.16E-04).

Eleven variants were in common among the different serum variables within the Diabetic Cohort (Supplementary Table 21). In particular, RBKS-rs13023003C, ERBB4-rs10168303A and GATM-15:45672447_GAA_GG (variant included the GATM-rs58764877 LD block, associated with CKD in the Overall and Non-Diabetic Cohorts) were associated with most of the variables derived from SCr. These three genotyped variants were associated with higher SCr levels and consequently lower values of eGFRCrea. Directions of effects across cohorts were consistent.

UK-ROI

Three NEMG variants were associated with at least one phenotype in the UK-ROI collection with a p-value < 1.2e-6 (Supplementary Table 22). The G-allele of AGXT2-rs71615838 was associated with lower risk of DKD, SURF1-rs183853102A increased the risk of ESKD and TFB1M-rs71575026C increased eGFRCrea (Supplementary Table 22).

The AGXT2-rs71615838G variant, associated with lower risk of DKD in the UK-ROI collection, showed a different effect in the Diabetic Cohort of the UKB (increased ESKD risk and SCysC levels, decreased eGFRCysC and eGFRCreaCysC), but this was not significant after multiple comparisons correction (Table 2).

Table 2 Comparison of the values for variants in nuclear-encoded mitochondrial genes associated with any variable in the UK-ROI collection with those in the UKB Cohort (only for associations with p < 0.05)
Full size table

The SURF1-rs183853102A variant, associated with higher risk of ESKD in the UK-ROI collection, showed consistent effects in the Diabetic Cohort of the UKB (increased CKD risk and SCr levels, decreased eGFRCrea and eGFRCreaCysC), but not significant after multiple comparisons correction (Table 2). Combined OR for random effects meta-analysis for ESKD was 4.656 (SE:3.139; P = 0.1378; Supplementary Table 23).

The TFB1M-rs71575026C variant, associated with higher eGFRCrea in the UK-ROI collection, showed the opposite effect in the UKB cohorts (non-statistically significant).

Discussion

CKD is a complex heterogeneous disease with a strong genetic component [9,10,11]. However, CKD heritability is not fully accounted for by current GWAS data [12, 13], indicating that additional genetic factors may be responsible for CKD susceptibility [14, 15]. Mitochondria are crucial to kidney health [25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40], but genetic variation in mtDNA [41,42,43,44,45] and in NEMG [52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70] has not been explored fully [14]. This paper assessed NEMG and mtDNA variants for their association with CKD and kidney phenotypes in individuals with and without diabetes, the most important global risk factor for CKD. Main results are depicted in Fig. 3.

Fig. 3
figure 3

Genes with the associations most consistently shown across cohorts and phenotypes

Full size image

Mitochondrial variants

MT-ND5 codes for the mitochondrially encoded NADH:Ubiquinone Oxidoreductase Core Subunit 5, essential for the catalytic activity and assembly of complex I, which catalyses electron transfer from NADH through the respiratory chain, using ubiquinone as an electron acceptor [101]. This gene has been associated with eye diseases as Leber Hereditary Optic Neuropathy (LHON) [102]. Another variant in this gene, rs267606894A (12770A > G), has also been associated with age-related macular degeneration (AMD) in 17,832 controls and 16,144 advanced AMD cases of European ancestry from the International AMD Genomics Consortium (IAMDGC) dataset [103]. In kidney, two MT-ND5 pathogenic variants, the m.13513G > A and the m.13514A > G, have been involved in adult-onset kidney disease in three unrelated patients [45]. In our study, four polymorphisms in MT-ND5 (rs2853499, rs28359172, rs2853503 and rs28359178) were associated with eGFR both in the Overall and Non-Diabetic Cohorts. A recent analysis in the UKB by Yonova-Doing et al. showed 16 mitochondrial variants associated with parameters of renal function (SCr, SCysC and eGFR) [51], 10 of them confirmed by our results (MT-HV2-rs869183622, MT-RNR2-rs2854128, MT-RNR2-rs3928306, MT-RNR2-rs2854131, MT-CO1-7028C, MT-ND5-rs2853499, MT-ND5-rs2853503, MT-CYB-rs2853506, MT-5-rs2857290 and MT-ND4-rs2853493). Interestingly, our analysis, not only confirms the association of rs2853499 and rs2853503 variants in MT-ND5 with serum variables, but also newly identifies the MT-ND5-rs41535848G variant, along with the mitochondrial haplogroup X, associated with higher risk of ESKD in the Overall/Non-Diabetic Cohorts, emphasising the role of this gene in kidney disease. Two variants in MT-ND5, G12372A (rs2853499) and G13708A (rs28359178) were associated with practically all serum variables in the Overall and Non-Diabetic cohorts, along with mitochondrial haplogroup U, in line with the study by Yonova-Doing et al. [51]. Another defining mutation for haplogroup U, MT-ND4-G11467A (rs2853493), was associated with eGFRCreaCysC in the Non-Diabetic Cohort. MT-ND4 codes the core subunit 4 of the mitochondrial membrane respiratory chain NADH dehydrogenase (Complex I), which catalyses electron transfer from NADH through the respiratory chain, using ubiquinone as an electron acceptor, essential for the catalytic activity and assembly of complex I [101, 104, 105].

In our study, mitochondrial haplogroup H, along with its defining mutations MT-RNR2-G2706A (rs2854128) and MT-CO1-T7028C, were associated with all the serum variables across all cohorts, confirming the results by Yonova-Doing et al. for eGFR, SCr and SCysC [51]. The study by Yonova-Doing et al. provided a comprehensive atlas of mitochondrial associations, for a set of 877 complex traits, including kidney phenotypes, but mostly related to kidney function, such as SCr, SCysC and related eGFR [51]. Our study, with a focus on kidney diseases, adds new and specific evidence on CKD, DKD and ESKD and those kidney conditions that may cause injury in the kidneys through the ‘kidney damage’ composite variable. Therefore, our analyses provide more precise evidence on the susceptibility of individuals to kidney disease in different stages and also in a more comprehensive context, as we stratified by diabetes to account for the influence of this condition and identify markers of susceptibility in this specific group of patients.

SNPs in the D-loop were previously identified as potential predictors of kidney survival and poor disease outcomes in 119 CKD patients and 159 controls [41,42,43]. Among the D-loop mitochondrial variants proposed, the MT-HV2-rs869183622A-allele increased CKD susceptibility (99.1% vs 0%; p < 0.001) [41], while the analysis in the UKB by Yonova-Doing et al. showed association with SCr and eGFRCreaCysC (p < 5·10–5) [51]. Our analysis showed a consistent effect of MT-HV2-rs869183622A increasing SCr and SCysC levels and decreasing their corresponding eGFR equations.

NEMG variants

Among the NEMG investigated in our study, NAT8-rs13538 and CPS1-rs1047891 stood out not only as predictors of CKD, but also showing associations with all the serum variables in the Overall and Non-Diabetic cohorts. The influence of CPS1-rs1047891A variant on lower SCr levels and consequently higher scores of eGFR equations based on SCr across was consistent among all cohorts. NAT8 codes for the N-acetyltransferase 8 (NAT8) enzyme, which catalyses the last step of mercapturic acid formation by acetylating cysteine S-conjugates to mercapturic acids [106] and plays an important role in the development and maintenance of normal kidney and liver structure and function [107], where is abundantly and specifically expressed [107], in particular by tubular cells of the renal cortex [108]. NAT8 (rs13538) is a known susceptibility locus for CKD and kidney function [13, 52, 54,55,56, 108,109,110]. The variant rs10206899 (close to NAT8 and in LD with rs13538, R2 = 0.988) was associated with SCr, eGFR, SCysC and CKD in a GWAS meta-analysis of nine studies encompassing 23,812 European white participants, by Chambers et al. [108]. NAT8-rs13538 results in a non-conservative amino acid change (F143S) within the acetyl-coenzyme A binding site, an effect predicted to influence acetylation by NAT8 [108], a key metabolic pathway for the detoxification of nephrotoxic substances [111, 112]. Association with serum metabolites was reported for NAT8-rs13538 (N-acetylornithine) and GATM-rs2433610 (located in 15:45,686,091, relatively close to GATM-15:45672447_GAA_GG) and CPS1-rs1047891 (glycine) in 1,260 African Americans from the Atherosclerosis Risk in Communities (ARIC) study [113]. These SNP-metabolite associations had also been seen in Europeans, in participants of the German KORA F4 study (n = 1,768) and the British TwinsUK study (n = 1,052) [114]. The fact that the risk allele associated here with higher serum levels of N-acetylornithine that is also associated with higher risk of CKD has pinpointed a role for ornithine acetylation in the aetiology of CKD [114]. In support of this, the circulating and urinary levels of 14 N-acetylated amino acids were associated with NAT8-rs13538 variant in 962 participants of the African American Study of Kidney Disease and Hypertension, 1050 from the Atherosclerosis Risk in Communities study and 680 from the electronic health record-linked biorepository BioMe [115, 116]. Higher circulating levels of five of these N-acetylated amino acids predicted kidney failure in the combined meta-analysis [115]. However, none of the urinary levels of these N-acetylated amino acids were associated with kidney failure in 1624 participants from the German CKD study [115].

CPS1 codes for the carbamoyl phosphate synthetase I, the rate-limiting enzyme catalysing the first committed step of the hepatic urea cycle by synthesizing carbamoyl phosphate from ammonia, bicarbonate, and 2 molecules of ATP [117]. CPS1-rs1047891, previously associated with eGFRCrea [71, 109], was associated with glycine levels among African Americans from the Atherosclerosis Risk in Communities (ARIC) Study (P = 4E-12) [113]. Glycine is metabolically related to carbamoyl phosphate, the product of CPS1 and the entry point of ammonia into the urea cycle. Other SNPs in this gene were also associated with glycine concentrations in 2,820 individuals from two large population-based European cohorts, as KORA F4 (rs2371015; P = 3E-09; R2 = 0.276) and TwinsUK (rs4673553: P = 2E-23; R2 = 0.461; rs4673558: P = 4.3E-11; R2 = 0.256) [114]. In the meta-analysis, CPS1-rs2216405 variant (LD with rs1047891: R2 = 0.385) was also associated with serum levels of glycine, along with serum levels of creatine, produced from glycine, probably indicating an altered ammonia metabolism [114]. Another SNP, located in the 3′ untranslated region of the CPS1 gene, rs715 (in LD with rs1047891: R2 = 0.907), previously associated with CKD [109], also showed an association with glycine levels in 1,004 nondiabetic individuals from the RISC study (P = 3.3E-50) [118]. The influence of this SNP on glycine can be seen not only in serum. In urine samples from 3,861 participants of the SHIP-0 cohort and 1,691 subjects of the KORA F4 cohort, CPS1-rs715C was associated with higher glycine/threonine ratio (Beta: 0.141; P = 8.5E-31; SHIP-0 dataset) [119]. In the same study, AGXT2-rs37369T predicted higher urinary levels of 3-aminoisobutyrate in 3,828 individuals of the SHIP-0 dataset (Beta: 1.277; P = 7.5E-26) [119], confirming previous results in urine [120,121,122] and serum [123]. Other studies had also related this variant with increased serum levels of symmetric/asymmetric dimethylarginine [124] and decreased homoarginine [125]. In our study, a different variant in the same gene, AGXT2-rs71615838, predicted higher risk of DKD in T1DM patients from the UK-ROI dataset; however, the effect was not seen in the Diabetic Cohort of the UK Biobank.

In our study, GATM-rs58764877 was associated with CKD in the Overall and Non-Diabetic Cohorts, and GATM-15:45672447_GAA_GG (in LD with rs58764877) with most of the variables derived from SCr in the Diabetic Cohort. In addition to the GATM variant discussed above (rs2433610), another variant in SPATA5L1, close to GATM, was associated with eGFRCrea in 2,388 CKD cases included among the participants of four different population-based cohorts of European-ancestry (Beta: -0.013; SE: 0.002; P = 6.2E-14) [54]. GATM encodes the enzyme L-arginine:glycine amidinotransferase, involved in creatine biosynthesis. Although SNPs at this locus have been proposed to affect levels of SCr without influencing susceptibility to kidney disease [54], our results show GATM variants can predict both SCr levels and CKD.

In our participants, two variants in SLC22A2 showed an influence on kidney phenotypes. In particular, SLC22A2-rs3127573G increased CKD risk in the Overall Cohort, and SLC22A2-rs1554261092(delT) lowered SCr levels, with consequently higher eGFRCrea. SLC22A2 has been previously associated with kidney traits. The SLC22A2-rs3127573 variant (located in 6:160,681,393) was associated with SCr and eGFR in the GWAS meta-analysis by Chambers et al. [108]. Other variants as have been associated with eGFRCrea and/or CKD [57, 109]. SLC22A2 codes for the solute carrier family 22 member 2, that are polyspecific organic cation transporters expressed in the liver, kidney, intestine, and other organs and are critical for elimination of many endogenous small organic cations as well as a wide array of drugs and environmental toxins [126].

Genes like SLC22A2 [57, 108, 109, 127, 128], CPS1 [109], WDR72 and WDR37 [57, 109, 127,128,129] have shown specific association with eGFRCrea, therefore have been proposed to be involved in creatinine secretion, rather than kidney function, and consequently may not be representative of CKD. Functional analysis and other measures of eGFR, such as SCysC, may throw some light on this, to elucidate whether variants associated with creatinine biosynthesis or secretion are really relevant to predict CKD [12]. In our patients, polymorphisms in these three genes were associated with kidney traits not only derived from SCr, but also from SCysC, such as eGFRCreaCysC, eGFRCysC and SCysC, along with CKD, showing that their relationship may be indicative of kidney disease and function, not only creatinine production or secretion.

In our study, one of the loci predictive of CKD in both the Overall and Non-Diabetic Cohorts was WDR72-rs72747347. Variants in WDR72 (rs491567; R2 = 0.202 with rs72747347) have been previously associated with eGFRCrea (P = 2.7E-13) although not with eGFRCysC [109]. In 14,700 Japanese individuals, WDR72-rs10518733 (R2 = 0.193 with rs72747347) was associated with SCr levels (Beta: -0.068; SE: 0.012; P = 1.7E-08) [58]. The intronic WDR72-rs491567A variant (R2 = 0.202 with rs72747347) decreased eGFRCrea levels in up to 3,282 American Indian individuals from the Strong Heart Family Study (Beta: -0.09; SE: 0.03; P = 4.5E-04) [59]. WDR72-rs17730436 (R2 = 0.193 with rs72747347) was associated with both SCr (Beta: 0.0005; SE: 0.0007; P = 1.2E-13) and eGFRCrea (Beta: -0.0057; SE: 0.0009; P = 6E-13) and WDR72-rs17730281 (R2 = 0.201 with rs72747347 but R2 = 0.915 with rs17730436) with blood urea nitrogen (BUN) (Beta: 0.0051; SE: 0.0008; P = 3E-11) in a GWAS meta-analysis including up to 71,149 East Asian subjects to investigate kidney function-related traits [60]. In 490 unrelated Emirati nationals with type 2 diabetes mellitus (T2DM), variants in WDR72 showed a trend to association with SCr (rs1031755) or eGFRCrea (rs4776168 and rs10518733), but it was not significant after correction for multiple comparisons (R2 < 0.2 with rs72747347) [61]. The fact that multiple associations within WDR72 have been shown across different ethnicities suggest that this gene may have relevance for kidney function across multiple ancestries. A recent transethnic GWAS meta-analysis, including participants from the UK and Japan Biobanks, to understand the common genetic factors contributing to nephrolithiasis, identified a potential role of the intronic rs578595 variant in WDR72 in calcium-sensing receptor (CaSR) signalling [130]. WDR72, WD repeat-containing protein 72, encodes a protein with eight WD-40 repeats, highly expressed in kidney [58], which plays an important role in enamel mineralization, possibly due to endocytic vesicle trafficking [131, 132] and is thought to play a role in clathrin-mediated endocytosis, a central process to sustained intracellular CaSR signaling [132]. WDR72 is also causative of hereditary distal renal tubular acidosis (dRTA), a rare genetic disease, in an autosomal recessive manner [133,134,135]. One of the genes strongly implicated in the pathogenesis of dRTA is SLC4A1 [133, 136], encoding for the solute carrier family 4 member 1, part of the anion exchanger family and expressed in the erythrocyte plasma membrane, where it functions as a chloride/bicarbonate exchanger involved in carbon dioxide transport from tissues to lungs. In participants from the Overall Cohort, SLC4A1-rs116844389A was associated with ESKD, but not with serum variables. Renal complications of dRTA include nephrocalcinosis, nephrolithiasis, medullary cysts, and impaired kidney function, and it is common to develop moderate to severe CKD over time [137, 138]. WDR72-rs77593734T (OR:1.102; P = 1.4E-11; R2 = 0.251 with rs72747347) and rs690428A (OR:1.078; P = 1.5E-5; R2 = 0.279 with rs72747347) variants were recently identified to increase rapid eGFRCrea decline in a GWAS meta-analysis including 42 studies [139]. In the same study, WDR72-rs77593734 was also associated with eGFRSCysC (P = 1.9E-16) and BUN in UKB participants [139], confirming previous results in BUN for 416,076 participants from the Chronic Kidney Disease Genetics (CKDGen) Consortium (P = 8.5E-17) [71]. In Africans, rs12906891 (R2 = 0.209 with rs72747347) and rs11070992 (R2 = 0.210 with rs72747347) variants in WDR72 have been associated with proliferative diabetic retinopathy (PDR), a sight-threatening complication of diabetes that is associated with longer duration of diabetes and poor glycemic control (P = 9.7E-10; OR:1.46 and P = 4.2E-08; OR:1.28, respectively) [140]. WDR72-rs551225A (R2 = 0.747 with rs72747347) was associated with low urine pH (Beta: -0.03; CI95%: -0.03-(-0.02); P = 2.6E-15) and increased risk of kidney stones in a set of 150,274 Icelanders (OR:1.09; CI95%: 1.06–1.12; P = 4.8E-08) [141].

In our study, ERBB4-rs10168303A, an intronic variant located between introns 1–2 of human ERBB4, was associated with most of the variables derived from SCr in the Diabetic Cohort of the UKB, showing higher SCr levels and consequently lower values of eGFR. The ERBB4 gene is a member of the type I receptor tyrosine kinase subfamily, encoding a receptor for NDF/heregulin (NRG1) [142]. Recently, Erbb4-IR, a lncRNA located within the intron region between the first and second exons of ERBB4 on chromosome 1 of the mouse genome, has been found to induce renal fibrosis [143] and to promote diabetic kidney Injury in db/db Mice by targeting miR-29b [144]. This is not the first time that a variant in this gene has been associated with kidney disease in patients with diabetes. The top SNP associated with the primary DKD phenotype identified in a combined GWAS meta-analysis of discovery and second stages in DKD individuals (6, 691) was rs7588550 (R2 < 0.1 with rs10168303), an intronic SNP in the ERBB4 gene, which demonstrated consistent protective effects in the replication samples (OR = 0.66; CI95%: 0.56–0.77; P = 2.1E-07) [76].

We also investigated potential eQTLs among our findings. The SLC39A13-rs10742802 variant was associated with kidney damage and the serum variables eGFRCrea, eGFRCreaCysC, SCr and SCysC in both the Overall and Non-Diabetic Cohorts. The SLC39A13-rs2293576 variant, in strong LD with SLC39A13-rs10742802 (R2 = 0.763) is a known sQTL in kidney cortex. SLC39A13, a member of the SLC39A family of zinc transporters, responsible for zinc influx [145], is over-expressed in kidney during dietary zinc deficiency in Wistar rats, playing a role in transporting zinc into cells to avoid zinc depletion [146]. Kidneys contribute to zinc homeostasis in the body by reabsorbing the portion of zinc entering the glomerular filtrate in the nephron [147]. In nonpathological conditions, zinc urinary loss is minimised; however, abnormal kidney function and diabetes, among other conditions, result in reduced serum zinc concentrations and increased urinary excretion [148].

Two other SNPs (LMNA-rs4641 and TRMT1-rs35601737) associated with different kidney phenotypes in our study, were also found to be sQTL in kidney cortex. The LMNA gene encodes the structural protein components of the nuclear lamina, lamin A and lamin C a protein network underlying the inner nuclear membrane that determines nuclear shape and size. Pathogenic variants in LMNA are highly pleiotropic and are responsible for many laminopathies [149] but have also been linked to other unrelated phenotypes, such as CKD (P = 1.13E-06), even independently of primary cardiomyopathy (P = 1.33E-03) in 11,451 unselected individuals from the Penn Medicine Biobank (Pennsylvania, US), suggesting an independent pathophysiological mechanism for renal failure in the context of loss of function in LMNA [150]. Germline LMNA mutations carried across generations in focal segmental glomerulosclerosis (FSGS) patients also points to a physiological role of LMNA in the maintenance of glomerular structure and function [151]. TRMT1 has been recently validated as part of a 5-gene prognostic signature for kidney renal clear cell carcinoma (KIRC) in The Cancer Genome Atlas (TCGA) and E-MTAB-1980 cohorts [152, 153].

Among the novel potential associations of NEMG with kidney traits found in our study, the PTPN11-rs11614544 variant was consistently associated with SCysC levels and consequently eGFRSCysC. Interestingly, MYH14-rs148695576T and NUP210-rs144856263T were markers associated with a higher risk of ESKD. We also had a particular interest in variants that could identify kidney damage. Among the NEMG predictive of kidney damage along with serum phenotypes in our patients (NOS3, ACP2 and SLC39A13), some of them have previously been linked to kidney disease. The transancestry meta-analysis of GWAS for eGFR performed Wuttke et al. [71], using nearly a million individuals, identified 147 loci relevant for kidney function based on associations with the alternative kidney function marker BUN, which subsequently were tested as a GRS in clinically diagnosed CKD and CKD-related outcomes in the UKB. None of the novel signals in our study made it to the GRS in the UKB, but interestingly, SQOR-rs629024 and SLC239A13-rs10742802 were associated with eGFRSCr (P = 7.5E-31 and 1. 9E-10, respectively). Furthermore, when these results were later meta-analysed joining data from the UKB (n = 1,201,909) [154], among their results for eGFRSCys, ATP5G1-rs1800632 and PTPN11-rs11614544 showed p-values of 5.6E-11 and 2.3E-23 respectively, whereas SQOR-rs629024 and SLC239A13-rs10742802 were confirmed for eGFRSCr (P = 7.4E-34 and 1.4E-12, respectively), although these associations were not discussed in the article. The meta-analysis also revealed NOS3-rs3918226 associated with eGFRSCr (P = 5.4E-09) and SLC239A13-rs10742802 with eGFRSCys (P = 3.4E-15) [154]. For other variants in our study, there was not information available from previous studies, despite some of the GWAS analysis having been performed on UKB data. There are several reasons why the novel associations identified in our study may not have been seen in previous GWAS of kidney traits. Different study designs, outcome measures and/or different ancestries may have generated different results to our study. The candidate-gene design of our study, including a lower number of variants than a GWAS, has allowed us to use a threshold of 9E-07 to identify associations, lower than the widely accepted GWAS-significant threshold of 5E-08. This may have caused these associations to have been deemed non-significant in such studies. For instance, Kintu et al. recently performed a GWAS meta-analysis to identify susceptibility loci associated with eGFRcrea in 80,027 individuals of African-ancestry from the UKB, Million Veteran Program, and CKDGen consortia [155]. Despite the different ancestry, this study identified 8 lead SNPs, 7 previously associated with eGFRCrea in other populations. However, the authors use a threshold for significance of 5E-08, that may have missed potential variants in common with our study. The summary statistics for these analyses are not available, therefore we cannot confirm whether our associations may have p-values close to the threshold in this case.

It is remarkable that CFL1-rs117624356 was associated with all the serum variables in our Overall Cohort, especially given that cofilin-1, encoded by CFL1, has been proposed to be integral to the development of proteinuria, being necessary for modulating actin dynamics in podocytes [156, 157]. Podocyte alterations in actin architecture may initiate or aid the progression of glomerular diseases [158]. Cofilin1 has also been involved in hypertensive nephropathy by modulating the nuclear translocation of NF-κB and the expression of its downstream inflammatory factors in renal tubular epithelial cells [159]. Up-regulation of cofilin-1 in HK-2 cells treated with calcium oxalate monohydrate, the major crystalline composition of most kidney stones, suggests it may play a key role in response to the COM crystals adhesion [160]. Animals with mutations in CFL1 usually display abnormalities including ureter duplication, renal hypoplasia, and abnormal kidney shape [161].

Two variants in NOS3 (rs3918226 and rs891511) were associated with kidney damage in the UKB cohort. NOS3-rs3918226 was also associated with some serum traits. NOS3 encodes the nitric oxide synthase 3, involved in endothelin pathways and the EGF/EGFR signalling pathway, whose inhibition protects against the development of DKD [162]. There is a clear link between NOS3 gene variation and ESKD [62, 63], CKD [64] and CKD progression [65]. A recent meta-analysis of 13 studies found that two NOS3 gene polymorphisms (rs1799983, R2 = 0.189 with rs3918226 and the intron 4 VNTR a/b polymorphism) significantly increased ESKD risk in autosomal dominant polycystic kidney disease (ADPKD) patients [62]. The rs1799983 variant has also showed a role in development of CKD in a recent meta-analysis [64]. The NOS3-rs2070744 gene polymorphism (R2 = 0.140 with rs3918226) increased risk of ESKD among 100 Egyptian patients compared to 100 healthy controls (P < 0.001) [63]. NOS3-rs7830 (R2 = 0.494 with rs891511) has been associated with the risk of DKD in T2DM patients of Greek Caucasian origin (OR: 1.598; CI95%: 1.152–2.217; 121 DKD/220 T2DM) [66]. SNPs in NOS3 have been proposed as potential molecular markers to predict the risk of T2DM and DKD in Chinese Han population [67]. In 490 T2DM patients and 485 healthy controls, NOS3-rs3918188 was associated to susceptibility to T2DM; the rs1800783 polymorphism (R2 = 0.1404 with rs3918226) predicted DKD, and family history of diabetes was closely associated with rs11771443 (R2 < 0.1 with rs3918226) polymorphism in DKD [67]. The ACP2 protein was included in a prognostic model of DKD including 35 DKD patients with good and 19 with poor prognosis [68]. On the other hand, associations in ATP5G1, the ATP Synthase Membrane Subunit C Locus 1 represents one of the novel findings in our study [62].

One of the aims of our study was to explore the impact of diabetes on the association of mtDNA and NEMG variants with kidney phenotypes. Remarkably, our findings have made evident that the risk of kidney outcomes may be exacerbated by the presence of diabetes. That is the case of TBC1D32, not previously linked to kidney traits, accentuating a high risk of ESKD and emerging as a potential predictor of worse outcome in patients with diabetes.

Our study also identified three NEMG predictive of DKD and ESKD among T1DM patients, two of them for the first time. In the UK-ROI collection, individuals carrying the AGXT2-rs71615838G variant had decreased DKD risk, whereas carrying the SURF1-rs183853102A allele was associated with an increased risk of ESKD. Also, individuals with the TFB1M-rs869120C allele had lower serum levels of eGFR. AGXT2 encodes the enzyme alanine:glyoxylate aminotransferase 2, expressed primarily in the kidney, which converts L-homoarginine (hArg) into 6-guanidino-2-oxocaproic acid (GOCA) [163]. Decreased plasma concentrations of hArg have recently become an emerging marker for clinical status and prognosis in kidney disease [164,165,166]. In 527 patients with different stages of CKD from the CARE FOR HOMe study confirmed this and found that a decreased ratio between hArg and GOCA predicted even more pronounced risks for kidney events [69]. The mitochondrial aminotransferase encoded by AGXT2 also catalyses the reaction of β-aminoisobutyrate with pyruvate to form 2-methyl-3-oxopropanoate and alanine.

Overall, our study has found many associations, both in mitochondrial and nuclear genes, with potential implications for the identification of patients with higher susceptibility to kidney diseases. In particular for ESKD, we have identified three genes that might help identify people with a higher risk of developing this more aggressive form of kidney disease, both in the overall population (NUP210 and SLC4A1) and in people with diabetes (SURF1). Their actual impact on risk stratification in kidney disease needs further investigation. However, given that most of the effect sizes, especially for the mitochondrial genes, were small, their use in isolation may be limited. To improve kidney disease risk prediction may require carefully validated combinations of clinical risk factors and a variety of ‘-omics’ based predictors.

Limitations

This is a well-powered cross-sectional study exploring an extensive selection of kidney phenotypes in one of the largest populations cohorts available, the UKB cohort. However, it is not exempt of limitations. The UK-ROI collection and the Diabetic Cohort within the UKB are small compared to the overall UKB cohort and may have reduced the power to detect significant associations. Furthermore, most individuals in the UKB Diabetic cohort had type 2 diabetes, whereas the UK-ROI collection is composed exclusively of T1DM individuals, which may partially explain why findings in T1DM from the UK-ROI collection were not replicated in the UKB.

The eGFR was calculated using the 2009 CKD-EPI creatinine equation, which uses coefficients for age, sex, and race in addition to SCr or SCysC, a premise that has recently come under scrutiny and criticism [167, 168], since race adjustment may have overestimated eGFR (≈16% higher) in Black patients, systematically resulting in an unfavourable bias that potentially may have reduced their access to pre-emptive kidney transplantation. As our study only included participants of European ancestry, no correction for race was applied in the eGFR equations or as covariate in the association analyses, therefore the ethnicity issue should not be a concern.

Conclusions

In summary, our study has brought to light the influence of variants both in mtDNA and NEMG which may explain some of the missing heritability in CKD. The consistent effects on eGFR and serum variables observed for mitochondrial variants and haplogroups, in particular haplogroup H, associated with all serum variables across all cohorts, support a role of variation in mitochondrial genes for kidney disease (Fig. 3). Our results confirm previously identified associations of NEMG (NAT8, GATM, SLC22A2, WDR72, NOS3, AGXT2 and CPS1) with kidney disease and kidney function, providing new information that expands these associations to other kidney disease biomarkers. We also provide new evidence for association of variants in SLC39A13, CFL1, ACP2, or ATP5G1, with serum phenotypes and kidney injury. We identified several associations with the most severe stage of CKD—ESKD, both in nuclear (SLC4A1, NUP210 MYH14) and mitochondrial genes (MT-ND5, defining variant for haplogroup X).

Our findings also highlight variants in TBC1D32 and SURF1 that are associated with a higher risk of ESKD in individuals with diabetes.

Nevertheless, there are still many ‘gaps’ in our knowledge, beyond SNP variations in nuclear and mitochondrial DNA, to fully understand the genetic predisposition to CKD.

Availability of data and materials

The primary data used and/or analysed during the current study are available from the UK Biobank. This research was conducted using the UK Biobank Resource under Application Number 14259. Detailed phenotyping information and summative data are presented in this paper.

Abbreviations

ADPKD:

Autosomal dominant polycystic kidney disease

BMI:

Body mass index

CKD:

Chronic kidney disease

CKD-EPI:

Chronic Kidney Disease Epidemiology Collaboration

DBP:

Diastolic blood pressure

DKD:

Diabetic kidney disease

D-loop:

Displacement loop

DM:

Diabetes mellitus

DNCRI:

Diabetic Nephropathy Collaborative Research Initiative

eGFR:

Estimated glomerular filtration rate

eGFRCrea :

eGFR based on serum creatinine

eGFRCreaCysC :

eGFR based on the combined serum creatinine and cystatin C equation

eGFRCysC :

eGFR based on serum cystatin C

eGFREPI :

Estimated glomerular filtration rate based on the CKD-EPI equation

eQTL:

Expression quantitative trait loci

ESKD:

End-stage kidney disease

GTEx:

Gene tissue expression

GWAS:

Genome-wide association studies

HGNC:

HUGO gene nomenclature committee

HWE:

Hardy–Weinberg equilibrium

KRT:

Kidney replacement therapy

MAC:

Minor allele count

MAF:

Minor allele frequency

mtDNA:

Mitochondrial DNA

NEMG:

Nuclear-encoded mitochondrial genes

PCA:

Principal component analysis

QC:

Quality control

rRNAs:

Ribosomal rRNAs

tRNAs:

Transfer RNAs

SBP:

Systolic blood pressure

SCr:

Serum creatinine

SCysC:

Serum cystatin C

SNP:

Single nucleotide polymorphism

sQTL:

Splicing quantitative trait loci

T1DM:

Type 1 diabetes mellitus

T2DM:

Type 2 diabetes mellitus

UKB:

United Kingdom Biobank

UK-ROI:

United Kingdom-Republic of Ireland

References

  1. Levey AS, Atkins R, Coresh J, Cohen EP, Collins AJ, Eckardt K-U, et al. Chronic kidney disease as a global public health problem: Approaches and initiatives – a position statement from Kidney Disease Improving Global Outcomes. Kidney Int. 2007;72:247–59.

    Article  CAS  PubMed  Google Scholar 

  2. Centers for Disease Control and Prevention. Chronic Kidney Disease Initiative | CDC. CKD Topics. https://www.cdc.gov/kidneydisease/index.html. Accessed 3 Dec 2018.

  3. Bash LD, Coresh J, Köttgen A, Parekh RS, Fulop T, Wang Y, et al. Defining incident chronic kidney disease in the research setting: The ARIC Study. Am J Epidemiol. 2009;170:414–24.

    Article  PubMed  PubMed Central  Google Scholar 

  4. GBD 2017 Disease and Injury Incidence and Prevalence Collaborators. G 2017 D and II and P. Global, regional, and national incidence, prevalence, and years lived with disability for 354 diseases and injuries for 195 countries and territories, 1990–2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet. 2017;2018(392):1789–858.

    Google Scholar 

  5. Foreman KJ, Marquez N, Dolgert A, Fukutaki K, Fullman N, McGaughey M, et al. Forecasting life expectancy, years of life lost, and all-cause and cause-specific mortality for 250 causes of death: reference and alternative scenarios for 2016–40 for 195 countries and territories. Lancet (London, England). 2018;392:2052–90.

    Article  PubMed  Google Scholar 

  6. GBD 2017 Causes of Death Collaborators. Global, regional, and national age-sex-specific mortality for 282 causes of death in 195 countries and territories, 1980–2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet. 2017;2018(392):1736–88.

    Google Scholar 

  7. Kerr M. Chronic Kidney Disease in England: The Human and Financial Cost. 2012. PMID: 22815543.

  8. Kazancioğlu R. Risk factors for chronic kidney disease: an update. Kidney Int Suppl. 2013;3:368–71.

    Article  Google Scholar 

  9. Regele F, Jelencsics K, Shiffman D, Paré G, McQueen MJ, Mann JFE, et al. Genome-wide studies to identify risk factors for kidney disease with a focus on patients with diabetes. Nephrol Dial Transplant. 2015;30(Suppl 4):iv26-3426–34.

    Article  Google Scholar 

  10. O’Seaghdha CM, Fox CS. Genome-wide association studies of chronic kidney disease: what have we learned? Nat Rev Nephrol. 2011;8:89–99.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Satko SG, Freedman BI. The Familial Clustering of Renal Disease and Related Phenotypes. Med Clin North Am. 2005;89:447–56.

    Article  CAS  PubMed  Google Scholar 

  12. Cañadas-Garre M, Anderson K, McGoldrick J, Maxwell AP, McKnight AJ. Genomic approaches in the search for molecular biomarkers in chronic kidney disease. J Transl Med. 2018;16:292.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Cañadas-Garre M, Anderson K, McGoldrick J, Maxwell AP, McKnight AJ. Proteomic and metabolomic approaches in the search for biomarkers in chronic kidney disease. J Proteomics. 2019;193:93–122.

    Article  PubMed  Google Scholar 

  14. Canadas-Garre M, Anderson K, Cappa R, Skelly R, Smyth LJ, McKnight AJ, 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 

  15. Anderson K, Canadas-Garre M, Chambers R, Maxwell AP, McKnight AJ. The Challenges of Chromosome Y Analysis and the Implications for Chronic Kidney Disease. Front Genet. 2019;10:781.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Cooper G. The Cell: A Molecular Approach. 2nd ed. Sunderland (MA): Sinauer Associates; 2000.

    Google Scholar 

  17. Chaban Y, Boekema EJ, Dudkina NV. Structures of mitochondrial oxidative phosphorylation supercomplexes and mechanisms for their stabilisation. Biochim Biophys Acta. 2014;1837:418–26.

    Article  CAS  PubMed  Google Scholar 

  18. Lodish H, Berk A, Zipursky S, Lawrence MP, Baltimore D, Darnell J. Electron Transport and Oxidative Phosphorylation  In: Molecular Cell Biology. 4th ed. New York: Freeman; 2012. p. 474.

    Google Scholar 

  19. Meiklejohn CD, Holmbeck MA, Siddiq MA, Abt DN, Rand DM, Montooth KL. An incompatibility between a mitochondrial tRNA and its nuclear-encoded tRNA synthetase compromises development and fitness in Drosophila. PLoS Genet. 2013;9: e1003238.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Taanman J-W. The mitochondrial genome: structure, transcription, translation and replication. Biochim Biophys Acta - Bioenerg. 1999;1410:103–23.

    Article  CAS  Google Scholar 

  21. Gray MW, Gray MW, Burger G, Lang BF. Mitochondrial Evolution. Science. 2008;1476:1–16.

    Google Scholar 

  22. Chan DC. Mitochondria: Dynamic Organelles in Disease, Aging, and Development. Cell. 2006;125:1241–52.

    Article  CAS  PubMed  Google Scholar 

  23. Dolezal P, Likic V, Tachezy J, Lithgow T. Evolution of the molecular machines for protein import into mitochondria. Science. 2006;313:314–8.

    Article  CAS  PubMed  Google Scholar 

  24. Timmis JN, Ayliffe MA, Huang CY, Martin W. Endosymbiotic gene transfer: organelle genomes forge eukaryotic chromosomes. Nat Rev Genet. 2004;5:123.

    Article  CAS  PubMed  Google Scholar 

  25. Skelly R, Maxwell A, McKnight A. Mitochondria and chronic kidney disease: a molecular update. SPG BioMed. 2019. https://doi.org/10.32392/biomed.36.

  26. Cappa R, de Campos C, Maxwell A, McKnight A. “Mitochondrial toolbox” - A review of online resources to explore mitochondrial genomics. Front Genet. 2020. PMID: 32457801.

  27. Douglas AP, Vance DR, Kenny EM, Morris DW, Maxwell AP, McKnight AJ. Next-generation sequencing of the mitochondrial genome and association with IgA nephropathy in a renal transplant population. Sci Rep. 2014;4:7379.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Rahman S, Hall AM. Mitochondrial disease–an important cause of end-stage renal failure. Pediatric nephrology (Berlin, Germany). 2013;28:357–61.

    Article  CAS  PubMed  Google Scholar 

  29. Wallace DC. A mitochondrial bioenergetic etiology of disease. J Clin Invest. 2013;123:1405–12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Swan EJ, Salem RM, Sandholm N, Tarnow L, Rossing P, Lajer M, et al. Genetic risk factors affecting mitochondrial function are associated with kidney disease in people with Type 1 diabetes. Diabet Med. 2015;32:1104–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Zhan M, Brooks C, Liu F, Sun L, Dong Z. Mitochondrial dynamics: regulatory mechanisms and emerging role in renal pathophysiology. Kidney Int. 2013;83:568–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Che R, Yuan Y, Huang S, Zhang A. Mitochondrial dysfunction in the pathophysiology of renal diseases. Am J Physiol Ren Physiol. 2014;306. PMID: 24305473.

  33. Galvan DL, Green NH, Danesh FR. The hallmarks of mitochondrial dysfunction in chronic kidney disease. Kidney Int. 2017;92:1051–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Ma L, Chou JW, Snipes JA, Bharadwaj MS, Craddock AL, Cheng D, et al. APOL1 renal-risk variants induce mitochondrial dysfunction. J Am Soc Nephrol. 2017;28:1093–105.

    Article  CAS  PubMed  Google Scholar 

  35. Braga PC, Alves MG, Rodrigues AS, Oliveira PF. Mitochondrial Pathophysiology on Chronic Kidney Disease. Int J Mol Sci. 2022;23:1776.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. De Seigneux S, Martin PY. Preventing the progression of AKI to CKD: The role of mitochondria. J Am Soc Nephrol. 2017;28:1327–9.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Szeto HH, Liu S, Soong Y, Seshan SV, Cohen-Gould L, Manichev V, et al. Mitochondria protection after acute ischemia prevents prolonged upregulation of IL-1b and IL-18 and arrests CKD. J Am Soc Nephrol. 2017;28:1437–49.

    Article  CAS  PubMed  Google Scholar 

  38. Tang C, Dong Z. Mitochondria in Kidney Injury: When the Power Plant Fails. J Am Soc Nephrol. 2016;27:1869–72.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Tin A, Grams ME, Ashar FN, Lane JA, Rosenberg AZ, Grove ML, et al. Association between mitochondrial DNA copy number in peripheral blood and incident CKD in the atherosclerosis risk in communities study. J Am Soc Nephrol. 2016;27:2467–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Govers LP, Toka HR, Hariri A, Walsh SB, Bockenhauer D. Mitochondrial DNA mutations in renal disease: an overview. Pediatric Nephrology. 2020. PMID: 31925537.

  41. Xu J, Guo Z, Bai Y, Zhang J, Cui L, Zhang H, et al. Single nucleotide polymorphisms in the D-loop region of mitochondrial DNA is associated with the kidney survival time in chronic kidney disease patients. Ren Fail. 2015;37:108–12.

    Article  CAS  PubMed  Google Scholar 

  42. Bai Y, Guo Z, Xu J, Zhang J, Cui L, Zhang H, et al. Association of sequence polymorphism in the mitochondrial D-loop with chronic kidney disease. Ren Fail. 2014;36:781–4.

    Article  CAS  PubMed  Google Scholar 

  43. Bai Y, Guo Z, Xu J, Zhang J, Cui L, Zhang H, et al. Single nucleotide polymorphisms in the D-loop region of mitochondrial DNA and age-at-onset of patients with chronic kidney disease. Chin Med J (Engl). 2014;127:3088–91.

    Article  PubMed  Google Scholar 

  44. Chen JB, Yang YH, Lee WC, Liou CW, Lin TK, Chung YH, et al. Sequence-based polymorphisms in the mitochondrial D-loop and potential SNP predictors for chronic dialysis. PLoS One. 2012;7. PMID: 22815937.

  45. Bakis H, Trimouille A, Vermorel A, Redonnet I, Goizet C, Boulestreau R, et al. Adult onset tubulo-interstitial nephropathy in MT-ND5-related phenotypes. Clin Genet. 2020;97:628–33.

    Article  CAS  PubMed  Google Scholar 

  46. Nass MMK. Mitochondrial DNA. I. Intramitochondrial distribution and structural relations of single- and double-length circular DNA. J Mol Biol. 1969;42. PMID: 5816965.

  47. Albring M, Griffith J, Attardi G. Association of a protein structure of probable membrane derivation with HeLa cell mitochondrial DNA near its origin of replication. Proc Natl Acad Sci U S A. 1977;74:1348–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Nagy A, Wilhelm M, Kovacs G. Mutations of mtDNA in renal cell tumours arising in end-stage renal disease. J Pathol. 2003;199:237–42.

    Article  CAS  PubMed  Google Scholar 

  49. Cappa R, Smyth L, Cañadas Garre M, de Campos CP, Skelly R, Cruise S, Young I, Kee F, McGuiness B, Godson C, Maxwell AP MA. Genetic and epigenetic analysis in genes affecting mitochondrial function reveal loci which are associated with chronic kidney disease in an older population. ASN Kidney Week 2018. 2018. https://www.asn-online.org/education/kidneyweek/2018/program-abstract.aspx?controlId=3014450.

  50. Skelly R. Next Generation Sequencing and Genome-Wide Association Studies to Identify Mitochondrial Genomic Features Associated with Diabetic Kidney Disease. 2020. https://pureadmin.qub.ac.uk/ws/portalfiles/portal/211677741/Next_Generation_Sequencing_and_Genone_Wide_Association_Studies_to_Identify_Mitochondrial_Genomic_Features_Associated.

  51. Yonova-Doing E, Calabrese C, Gomez-Duran A, Schon K, Wei W, Karthikeyan S, et al. An atlas of mitochondrial DNA genotype–phenotype associations in the UK Biobank. Nat Genet. 2021;53:982–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Cañadas-Garre M, Anderson K, McGoldrick J, Maxwell AP, McKnight AJ. Genomic approaches in the search for molecular biomarkers in chronic kidney disease 11 Medical and Health Sciences 1103 Clinical Sciences. Journal of Translational Medicine. 2018;16. PMID: 30359254.

  53. Böger CA, Heid IM. Chronic kidney disease: novel insights from genome-wide association studies. Kidney Blood Press Res. 2011;34:225–34.

    Article  PubMed  Google Scholar 

  54. Köttgen A, Glazer NL, Dehghan A, Hwang S-J, Katz R, Li M, et al. Multiple loci associated with indices of renal function and chronic kidney disease. Nat Genet. 2009;41:712–7.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Köttgen A. Genome-wide association studies in nephrology research. Am J Kidney Dis. 2010;56:743–58.

    Article  PubMed  Google Scholar 

  56. Tin A, Colantuoni E, Boerwinkle E, Kottgen A, Franceschini N, Astor BC, et al. Using Multiple Measures for Quantitative Trait Association Analyses: Application toEstimated Glomerular Filtration Rate (eGFR). J Hum Genet. 2013;58:461.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Pattaro C, Teumer A, Gorski M, Chu AY, Li M, Mijatovic V, 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 

  58. Kamatani Y, Matsuda K, Okada Y, Kubo M, Hosono N, Daigo Y, et al. Genome-wide association study of hematological and biochemical traits in a Japanese population. Nat Genet. 2010;42:210–5.

    Article  CAS  PubMed  Google Scholar 

  59. Okada Y, Sim X, Go MJ, Wu J-Y, Gu D, Takeuchi F, et al. Meta-analysis identifies multiple loci associated with kidney function-related traits in east Asian populations. Nat Genet. 2012;44:904–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Paterson AD, Waggott D, Boright AP, Hosseini SM, Shen E, Sylvestre MP, et al. A Genome-Wide Association Study Identifies a Novel Major Locus for Glycemic Control in Type 1 Diabetes, as Measured by Both A1C and Glucose. Diabetes. 2010;59:539.

    Article  CAS  PubMed  Google Scholar 

  61. Osman WM, Jelinek HF, Tay GK, Khandoker AH, Khalaf K, Almahmeed W, et al. Clinical and genetic associations of renal function and diabetic kidney disease in the United Arab Emirates: a cross-sectional study. BMJ Open. 2018;8. PMID: 30552240.

  62. Padhi UN, Mulkalwar M, Saikrishna L, Verma HK, Bhaskar L. NOS3 gene intron 4 a/b polymorphism is associated with ESRD in autosomal dominant polycystic kidney disease patients. J Bras Nefrol. 2022. https://doi.org/10.1590/2175-8239-JBN-2021-0089.

    Article  PubMed  PubMed Central  Google Scholar 

  63. Elsaid A, Samir eid O, Said SB, Zahran RF. Association of NOS3 (rs 2070744) and SOD2Val16Ala (rs4880) gene polymorphisms with increased risk of ESRD among Egyptian patients. J Genet Eng Biotechnol. 2021;19. PMID: 34661767.

  64. Gunawan A, Fajar JK, Tamara F, Mahendra AI, Ilmawan M, Purnamasari Y, et al. Nitride oxide synthase 3 and klotho gene polymorphisms in the pathogenesis of chronic kidney disease and age-related cognitive impairment: a systematic review and meta-analysis. F1000Research. 2020;9. PMID: 34035901.

  65. Medina AM, Zubero EE, Jiménez MAA, Barragan SAA, García CAL, Ramos JJG, et al. NOS3 Polymorphisms and Chronic Kidney Disease. J Bras Nefrol. 2018;40:273–7.

  66. Roumeliotis A, Roumeliotis S, Tsetsos F, Georgitsi M, Georgianos PI, Stamou A, et al. Oxidative Stress Genes in Diabetes Mellitus Type 2: Association with Diabetic Kidney Disease. Oxid Med Cell Longev. 2021;2021. PMID: 34545296.

  67. Chen F, Li Y-M, Yang L-Q, Zhong C-G, Zhuang Z-X. Association of NOS2 and NOS3 gene polymorphisms with susceptibility to type 2 diabetes mellitus and diabetic nephropathy in the Chinese Han population. IUBMB Life. 2016;68:516–25.

    Article  CAS  PubMed  Google Scholar 

  68. Ahn HS, Kim JH, Jeong H, Yu J, Yeom J, Song SH, et al. Differential Urinary Proteome Analysis for Predicting Prognosis in Type 2 Diabetes Patients with and without Renal Dysfunction. Int J Mol Sci. 2020;21:1–19.

    Article  Google Scholar 

  69. Martens-Lobenhoffer J, Emrich IE, Zawada AM, Fliser D, Wagenpfeil S, Heine GH, et al. L-Homoarginine and its AGXT2-metabolite GOCA in chronic kidney disease as markers for clinical status and prognosis. Amino Acids. 2018;50:1347–56.

    Article  CAS  PubMed  Google Scholar 

  70. Hu X-L, Zeng W-J, Li M-P, Yang Y-L, Kuang D-B, Li H, et al. AGXT2 rs37369 polymorphism predicts the renal function in patients with chronic heart failure. Gene. 2017;637:145–51.

    Article  CAS  PubMed  Google Scholar 

  71. Wuttke M, Li Y, Li M, Sieber KB, Feitosa MF, Gorski M, et al. A catalog of genetic loci associated with kidney function from analyses of a million individuals. Nat Genet. 2019;51:957–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. WHO :: Global Database on Body Mass Index. 2000. http://apps.who.int/bmi/index.jsp?introPage=intro_3.html. Accessed 12 Apr 2018.

  73. Bycroft C, Freeman C, Petkova D, Band G, Elliott LT, Sharp K, et al. The UK Biobank resource with deep phenotyping and genomic data. Nature. 2018;562:203–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Chang CC, Chow CC, Tellier LCAM, Vattikuti S, Purcell SM, Lee JJ. Second-generation PLINK: Rising to the challenge of larger and richer datasets. Gigascience. 2015;4:1–16.

    Article  Google Scholar 

  75. Chang C, GRAIL I, Human Longevity I, Department of Biomedical Data Science S. PLINK 2.00 alpha. 2020. https://www.cog-genomics.org/plink/2.0/.

  76. Sandholm N, Salem RM, McKnight AJ, Brennan EP, Forsblom C, Isakova T, et al. New susceptibility loci associated with kidney disease in type 1 diabetes. PLoS Genet. 2012;8: e1002921.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  77. Weissensteiner H, Pacher D, Kloss-Brandstätter A, Forer L, Specht G, Bandelt H-J, et al. HaploGrep 2: mitochondrial haplogroup classification in the era of high-throughput sequencing. Nucleic Acids Res. 2016;44:W58-63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  78. van Oven M. PhyloTree Build 17: Growing the human mitochondrial DNA tree. Forensic Sci Int Genet Suppl Ser. 2015;5:e392–4.

    Article  Google Scholar 

  79. Calvo SE, Clauser KR, Mootha VK. MitoCarta2.0: an updated inventory of mammalian mitochondrial proteins. Nucleic Acids Res. 2016;44:D1251-7.

    Article  CAS  PubMed  Google Scholar 

  80. Pagliarini DJ, Calvo SE, Chang B, Sheth SA, Vafai SB, Ong S-E, et al. A mitochondrial protein compendium elucidates complex I disease biology. Cell. 2008;134:112–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  81. Cotter D, Guda P, Fahy E, Subramaniam S. MitoProteome: mitochondrial protein sequence database and annotation system. Nucleic Acids Res. 2004;32(Database issue):D463-7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  82. Taylor SW, Fahy E, Zhang B, Glenn GM, Warnock DE, Wiley S, et al. Characterization of the human heart mitochondrial proteome. Nat Biotechnol. 2003;21:281–6.

    Article  CAS  PubMed  Google Scholar 

  83. Smith AC, Robinson AJ. MitoMiner v3.1, an update on the mitochondrial proteomics database. Nucleic Acids Res. 2016;44:1258–61.

    Article  Google Scholar 

  84. Brandon MC, Lott MT, Nguyen KC, Spolim S, Navathe SB, Baldi P, et al. MITOMAP: a human mitochondrial genome database—2004 update. Nucleic Acids Res. 2005;33(Database Issue):D611-3.

    Article  CAS  PubMed  Google Scholar 

  85. Zerbino DR, Achuthan P, Akanni W, Amode MR, Barrell D, Bhai J, et al. Ensembl 2018. Nucleic Acids Res. 2018;46:D754–61.

    Article  CAS  PubMed  Google Scholar 

  86. The UniProt Consortium. UniProt: the universal protein knowledgebase. Nucleic Acids Res. 2017;45:D158–69.

    Article  Google Scholar 

  87. Wain HM, Bruford EA, Lovering RC, Lush MJ, Wright MW, Povey S. Guidelines for human gene nomenclature. Genomics. 2002;79:464–70.

    Article  CAS  PubMed  Google Scholar 

  88. R Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing. 2022. https://www.R-project.org.

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  90. Higgins JPT, Thompson SG, Deeks JJ, Altman DG. Measuring inconsistency in meta-analyses. BMJ. 2003;327:557–60.

    Article  PubMed  PubMed Central  Google Scholar 

  91. Wuttke M, Wong CS, Wuhl E, Epting D, Luo L, Hoppmann A, et al. Genetic loci associated with renal function measures and chronic kidney disease in children: the Pediatric Investigation for Genetic Factors Linked with Renal Progression Consortium. Nephrol Dial Transplant. 2016;31:262–9.

    CAS  PubMed  Google Scholar 

  92. Kraja AT, Liu C, Fetterman JL, Graff M, Have CT, Gu C, et al. Associations of Mitochondrial and Nuclear Mitochondrial Variants and Genes with Seven Metabolic Traits. Am J Hum Genet. 2019;104:112–38.

    Article  CAS  PubMed  Google Scholar 

  93. Zhan X, Abecasis G. TabAnno. 2012. https://github.com/zhanxw/anno.

  94. Sherry ST. dbSNP: the NCBI database of genetic variation. Nucleic Acids Res. 2001;29:308–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  95. Dayem Ullah AZ, Oscanoa J, Wang J, Nagano A, Lemoine NR, Chelala C. SNPnexus: assessing the functional relevance of genetic variation to facilitate the promise of precision medicine. Nucleic Acids Res. 2018;46:W109–13.

    Article  PubMed  PubMed Central  Google Scholar 

  96. Dayem Ullah AZ, Lemoine NR, Chelala C. A practical guide for the functional annotation of genetic variations using SNPnexus. Brief Bioinform. 2013;14:437–47.

    Article  CAS  PubMed  Google Scholar 

  97. Dayem Ullah AZ, Lemoine NR, Chelala C. SNPnexus: a web server for functional annotation of novel and publicly known genetic variants (2012 update). Nucleic Acids Res. 2012;40:W65-70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  98. Chelala C, Khan A, Lemoine NR. SNPnexus: a web database for functional annotation of newly discovered and public domain single nucleotide polymorphisms. Bioinformatics. 2009;25:655–61.

    Article  CAS  PubMed  Google Scholar 

  99. Oscanoa J, Sivapalan L, Gadaleta E, Dayem Ullah AZ, Lemoine NR, Chelala C. SNPnexus: a web server for functional annotation of human genome sequence variation (2020 update). Nucleic Acids Res. 2020;48:W185–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  100. Kaplan J. fastDummies: Fast Creation of Dummy (Binary) Columns and Rows from Categorical Variables. 2020. https://cran.r-project.org/package=fastDummies.

  101. Bourges I, Ramus C, Mousson De Camaret B, Beugnot R, Remacle C, Cardol P, et al. tructural organization of mitochondrial human complex I: role of the ND4 and ND5 mitochondria-encoded subunits and interaction with prohibitin. Biochem J. 2004;383((Pt.3)):491–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  102. Crimi M, Galbiati S, Moroni I, Bordoni A, Perini MP, Lamantea E, et al. A missense mutation in the mitochondrial ND5 gene associated with a Leigh-MELAS overlap syndrome. Neurology. 2003;60:1857–61.

    Article  PubMed  Google Scholar 

  103. Persad PJ, Heid IM, Weeks DE, Baird PN, de Jong EK, Haines JL, et al. Joint Analysis of Nuclear and Mitochondrial Variants in Age-Related Macular Degeneration Identifies Novel Loci TRPM1 and ABHD2/RLBP1. Invest Ophthalmol Vis Sci. 2017;58:4027–38.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  104. Hofhaus G, Attardi G. Lack of assembly of mitochondrial DNA-encoded subunits of respiratory NADH dehydrogenase and loss of enzyme activity in a human cell mutant lacking the mitochondrial ND4 gene product. EMBO J. 1993;12:3043–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  105. De Vries DD, Went LN, Bruyn GW, Scholte HR, Hofstra RM, Bolhuis PA, et al. Genetic and biochemical impairment of mitochondrial complex I activity in a family with Leber hereditary optic neuropathy and hereditary spastic dystonia. Am J Hum Genet. 1996;58:703–11.

    PubMed  PubMed Central  Google Scholar 

  106. Veiga-da-Cunha M, Tyteca D, Stroobant V, Courtoy PJ, Opperdoes FR, Van Schaftingen E. Molecular identification of NAT8 as the enzyme that acetylates cysteine S-conjugates to mercapturic acids. J Biol Chem. 2010;285:18888–98.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  107. Ozaki K, Fujiwara T, Nakamura Y, Takahashi EI. Isolation and mapping of a novel human kidney- and liver-specific gene homologous to the bacterial acetyltransferases. J Hum Genet. 1998;43:255–8.

    Article  CAS  PubMed  Google Scholar 

  108. Chambers JC, Zhang W, Lord GM, van der Harst P, Lawlor DA, Sehmi JS, et al. Genetic loci influencing kidney function and chronic kidney disease. Nat Genet. 2010;42:373–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  109. Köttgen A, Pattaro C, Böger CA, Fuchsberger C, Olden M, Glazer NL, et al. New loci associated with kidney function and chronic kidney disease. Nat Genet. 2010;42:376–84.

    Article  PubMed  PubMed Central  Google Scholar 

  110. Böger CA, Gorski M, Li M, Hoffmann MM, Huang C, Yang Q, et al. Association of eGFR-Related Loci Identified by GWAS with Incident CKD and ESRD. PLoS Genet. 2011;7.

  111. Kharasch ED. Adverse drug reactions with halogenated anesthetics. Clin Pharmacol Ther. 2008;84:158–62.

  112. Lash LH, Fisher JW, Lipscomb JC, Parker JC. Metabolism of trichloroethylene. Environ Health Perspect. 2000;108(Suppl):177–200.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  113. Yu B, Zheng Y, Alexander D, Morrison AC, Coresh J, Boerwinkle E. Genetic determinants influencing human serum metabolome among African Americans. PLoS Genet. 2014;10: e1004212.

    Article  PubMed  PubMed Central  Google Scholar 

  114. Suhre K, Shin S-Y, Petersen A-K, Mohney RP, Meredith D, Wägele B, et al. Human metabolic individuality in biomedical and pharmaceutical research. Nature. 2011;477:54–60.

    Article  CAS  PubMed  Google Scholar 

  115. Luo S, Surapaneni A, Zheng Z, Rhee EP, Coresh J, Hung AM, et al. NAT8 Variants, N-Acetylated Amino Acids, and Progression of CKD. Clin J Am Soc Nephrol. 2020;16:37–47.

    Article  PubMed  PubMed Central  Google Scholar 

  116. Luo S, Feofanova EV, Tin A, Tung S, Rhee EP, Coresh J, et al. Genome-wide association study of serum metabolites in the African American Study of Kidney Disease and Hypertension. Kidney Int. 2021;100:430–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  117. Häberle J, Shchelochkov OA, Wang J, Katsonis P, Hall L, Reiss S, et al. Molecular defects in human carbamoy phosphate synthetase I: mutational spectrum, diagnostic and protein structure considerations. Hum Mutat. 2011;32:579–89.

    Article  PubMed  PubMed Central  Google Scholar 

  118. Xie W, Wood AR, Lyssenko V, Weedon MN, Knowles JW, Alkayyali S, et al. Genetic variants associated with glycine metabolism and their role in insulin sensitivity and type 2 diabetes. Diabetes. 2013;62:2141–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  119. Raffler J, Friedrich N, Arnold M, Kacprowski T, Rueedi R, Altmaier E, et al. Genome-Wide Association Study with Targeted and Non-targeted NMR Metabolomics Identifies 15 Novel Loci of Urinary Human Metabolic Individuality. PLoS Genet. 2015;11.

  120. Tanaka T, Shen J, Abecasis GR, Kisialiou A, Ordovas JM, Guralnik JM, et al. Genome-wide association study of plasma polyunsaturated fatty acids in the InCHIANTI Study. PLoS Genet. 2009;5: e1000338.

    Article  PubMed  PubMed Central  Google Scholar 

  121. Suhre K, Wallaschofski H, Raffler J, Friedrich N, Haring R, Michael K, et al. A genome-wide association study of metabolic traits in human urine. Nat Genet. 2011;43:565–9.

    Article  CAS  PubMed  Google Scholar 

  122. Rueedi R, Ledda M, Nicholls AW, Salek RM, Marques-Vidal P, Morya E, et al. Genome-wide association study of metabolic traits reveals novel gene-metabolite-disease links. PLoS Genet. 2014;10: e1004132.

    Article  PubMed  PubMed Central  Google Scholar 

  123. Rhee EP, Ho JE, Chen M-H, Shen D, Cheng S, Larson MG, et al. A genome-wide association study of the human metabolome in a community-based cohort. Cell Metab. 2013;18:130–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  124. Seppälä I, Kleber ME, Lyytikäinen LP, Hernesniemi JA, Mäkelä KM, Oksala N, et al. Genome-wide association study on dimethylarginines reveals novel AGXT2 variants associated with heart rate variability but not with overall mortality. Eur Heart J. 2014;35:524–30.

    Article  PubMed  Google Scholar 

  125. Kleber ME, Seppälä I, Pilz S, Hoffmann MM, Tomaschitz A, Oksala N, et al. Genome-wide association study identifies 3 genomic loci significantly associated with serum levels of homoarginine: the AtheroRemo Consortium. Circ Cardiovasc Genet. 2013;6:505–13.

    Article  CAS  PubMed  Google Scholar 

  126. Gründemann D, Schömig E. Gene structures of the human non-neuronal monoamine transporters EMT and OCT2. Hum Genet. 2000;106:627–35.

    Article  PubMed  Google Scholar 

  127. Sveinbjornsson G, Mikaelsdottir E, Palsson R, Indridason OS, Holm H, Jonasdottir A, et al. Rare mutations associating with serum creatinine and chronic kidney disease. Hum Mol Genet. 2014;23:6935–43.

    Article  CAS  PubMed  Google Scholar 

  128. Pattaro C, Köttgen A, Teumer A, Garnaas M, Böger CA, Fuchsberger C, et al. Genome-wide association and functional follow-up reveals new loci for kidney function. PLoS Genet. 2012;8: e1002584.

    Article  PubMed  PubMed Central  Google Scholar 

  129. Parsa A, Kanetsky PA, Xiao R, Gupta J, Mitra N, Limou S, et al. Genome-Wide Association of CKD Progression: The Chronic Renal Insufficiency Cohort Study. https://doi.org/10.1681/ASN.2015101152.

  130. Howles SA, Wiberg A, Goldsworthy M, Bayliss AL, Gluck AK, Ng M, et al. Genetic variants of calcium and vitamin D metabolism in kidney stone disease. Nat Commun. 2019;10:1–10.

    Article  CAS  Google Scholar 

  131. Katsura KA, Horst JA, Chandra D, Le TQ, Nakano Y, Zhang Y, et al. WDR72 models of structure and function: a stage-specific regulator of enamel mineralization. Matrix Biol. 2014;38:48–58.

    Article  CAS  PubMed  Google Scholar 

  132. Wang S, Hu Y, Yang J, Smith CE, Nunez SM, Richardson AS, et al. Critical roles for WDR72 in calcium transport and matrix protein removal during enamel maturation. Mol Genet genomic Med. 2015;3:302–19.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  133. Alexander RT, Law L, Gil-Peña H, Greenbaum LA, Santos F. Hereditary Distal Renal Tubular Acidosis. GeneReviews®. 1993.

  134. Rungroj N, Nettuwakul C, Sawasdee N, Sangnual S, Deejai N, Misgar RA, et al. Distal renal tubular acidosis caused by tryptophan-aspartate repeat domain 72 (WDR72) mutations. Clin Genet. 2018;94:409–18. https://www.ncbi.nlm.nih.gov/books/NBK547595/#:~:text=Hereditary%20dRTA%20includes%20both%20%22complete,a%20mild%20renal%20acidification%20defect.

  135. Jobst-Schwan T, Klämbt V, Tarsio M, Heneghan JF, Majmundar AJ, Shril S, et al. Whole exome sequencing identified ATP6V1C2 as anovel candidate gene for recessive distal renal tubular acidosis. Kidney Int. 2020;97:567.

    Article  CAS  PubMed  Google Scholar 

  136. D’Ambrosio V, Azzarà A, Sangiorgi E, Gurrieri F, Hess B, Gambaro G, et al. Results of a Gene Panel Approach in a Cohort of Patients with Incomplete Distal Renal Tubular Acidosis and Nephrolithiasis. Kidney Blood Press Res. 2021;46:469–74.

    Article  PubMed  Google Scholar 

  137. Watanabe T. Improving outcomes for patients with distal renal tubular acidosis: recent advances and challenges ahead. Pediatr Heal Med Ther. 2018;9:181–90.

    CAS  Google Scholar 

  138. Gómez-Conde S, García-Castaño A, Aguirre M, Herrero M, Gondra L, García-Pérez N, et al. Molecular aspects and long-term outcome of patients with primary distal renal tubular acidosis. Pediatr Nephrol. 2021;36:3133–42.

    Article  PubMed  Google Scholar 

  139. Gorski M, Jung B, Li Y, Matias-Garcia PR, Wuttke M, Coassin S, et al. Meta-analysis uncovers genome-wide significant variants for rapid kidney function decline. Kidney Int. 2021;99:926–39.

    Article  CAS  PubMed  Google Scholar 

  140. Liu C, Chen G, Bentley AR, Doumatey A, Zhou J, Adeyemo A, et al. Genome-wide association study for proliferative diabetic retinopathy in Africans. NPJ Genomic Med. 2019;4.

  141. Benonisdottir S, Kristjansson RP, Oddsson A, Steinthorsdottir V, Mikaelsdottir E, Kehr B, et al. Sequence variants associating with urinary biomarkers. Hum Mol Genet. 2019;28:1199.

    Article  CAS  PubMed  Google Scholar 

  142. Takahashi Y, Fukuda Y, Yoshimura J, Toyoda A, Kurppa K, Moritoyo H, et al. ERBB4 mutations that disrupt the neuregulin-ErbB4 pathway cause amyotrophic lateral sclerosis type 19. Am J Hum Genet. 2013;93:900–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  143. Feng M, Tang PMK, Huang XR, Sun SF, You YK, Xiao J, et al. TGF-β Mediates Renal Fibrosis via the Smad3-Erbb4-IR Long Noncoding RNA Axis. Mol Ther. 2018;26:148–61.

    Article  CAS  PubMed  Google Scholar 

  144. Sun SF, Tang PMK, Feng M, Xiao J, Huang XR, Li P, et al. Novel lncRNA Erbb4-IR Promotes Diabetic Kidney Injury in db/db Mice by Targeting miR-29b. Diabetes. 2018;67:731–44.

    Article  CAS  PubMed  Google Scholar 

  145. Eide DJ. The SLC39 family of metal ion transporters. Pflugers Arch. 2004;447:796–800.

    Article  CAS  PubMed  Google Scholar 

  146. Guo L, Hu X, Xu T, Qi X, Wan Y, Liu X, et al. Over-expression of Zip-13 mRNA in kidney and lung during dietary zinc deficiency in Wistar rats. Mol Biol Rep. 2011;38:1869–74.

    Article  CAS  PubMed  Google Scholar 

  147. Ranaldi G, Perozzi G, Truong-Tran A, Zalewski P, Murgia C. Intracellular distribution of labile Zn(II) and zinc transporter expression in kidney and MDCK cells. Am J Physiol - Ren Physiol. 2002;283(6):1365–75.

    Article  Google Scholar 

  148. Brandão-Neto J, Silva CAB, Shuhama T, Silva JA, Oba L. Renal handling of zinc in insulin-dependent diabetes mellitus patients. Biometals. 2001;14:75–80.

    Article  PubMed  Google Scholar 

  149. Worman HJ, Bonne G. “Laminopathies:” a wide spectrum of human diseases. Exp Cell Res. 2007;313:2121.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  150. Park J, Levin MG, Haggerty CM, Hartzel DN, Judy R, Kember RL, et al. A genome-first approach to aggregating rare genetic variants in LMNA for association with electronic health record phenotypes. Genet Med. 2020;22:102.

    Article  CAS  PubMed  Google Scholar 

  151. Thong KM, Xu Y, Cook J, Takou A, Wagner B, Kawar B, et al. Cosegregation of focal segmental glomerulosclerosis in a family with familial partial lipodystrophy due to a mutation in LMNA. Nephron Clin Pract. 2013;124:31–7.

    Article  CAS  PubMed  Google Scholar 

  152. Huang Y, Chen S, Xiao L, Qin W, Li L, Wang Y, et al. A Novel Prognostic Signature for Survival Prediction and Immune Implication Based on SARS-CoV-2-Related Genes in Kidney Renal Clear Cell Carcinoma. Front Bioeng Biotechnol. 2022;9.

  153. Wu Y, Wei X, Feng H, Hu B, Liu B, Luan Y, et al. Transcriptome Analyses Identify an RNA Binding Protein Related Prognostic Model for Clear Cell Renal Cell Carcinoma. Front Genet. 2021;11. PMID: 35141213.

  154. Stanzick KJ, Li Y, Schlosser P, Gorski M, Wuttke M, Thomas LF, et al. Discovery and prioritization of variants and genes for kidney function in >1.2 million individuals. Nat Commun. 2021;12. PMID: 33488680.

  155. Kintu C, Soremekun O, Machipisa T, Mayanja R, Kalyesubula R, Bagaya BS, et al. Meta-analysis of African ancestry genome-wide association studies identified novel locus and validates multiple loci associated with kidney function. BMC Genomics. 2023;24:1–9.

    Article  Google Scholar 

  156. Teng B, Lukasz A, Schiffer M. The ADF/Cofilin-Pathway and Actin Dynamics in Podocyte Injury. Int J Cell Biol. 2012;2012.

  157. Garg P, Verma R, Cook L, Soofi A, Venkatareddy M, George B, et al. Actin-depolymerizing factor cofilin-1 is necessary in maintaining mature podocyte architecture. J Biol Chem. 2010;285:22676–88.

    Article  PubMed  PubMed Central  Google Scholar 

  158. Kriz W, Lemley KV. The role of the podocyte in glomerulosclerosis. Curr Opin Nephrol Hypertens. 1999;8:489–97.

    Article  CAS  PubMed  Google Scholar 

  159. Wang QZ, Gao HQ, Liang Y, Zhang J, Wang J, Qiu J. Cofilin1 is involved in hypertension-induced renal damage via the regulation of NF-κB in renal tubular epithelial cells. J Transl Med. 2015;13:323.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  160. Wang Z, Li MX, Xu CZ, Zhang Y, Deng Q, Sun R, et al. Comprehensive study of altered proteomic landscape in proximal renal tubular epithelial cells in response to calcium oxalate monohydrate crystals. BMC Urol. 2020;20.

  161. Kuure S, Cebrian C, Machingo Q, Lu BC, Chi X, Hyink D, et al. Actin Depolymerizing Factors Cofilin1 and Destrin Are Required for Ureteric Bud Branching Morphogenesis. PLoS Genet. 2010;6:1–11.

    Article  Google Scholar 

  162. Li Y, Pan Y, Cao S, Sasaki K, Wang Y, Niu A, et al. Podocyte EGFR Inhibits Autophagy Through Upregulation of Rubicon in Type 2 Diabetic Nephropathy. Diabetes. 2021;70:562–76.

    Article  CAS  PubMed  Google Scholar 

  163. Rodionov RN, Oppici E, Martens-Lobenhoffer J, Jarzebska N, Brilloff S, Burdin D, et al. A Novel Pathway for Metabolism of the Cardiovascular Risk Factor Homoarginine by alanine:glyoxylate aminotransferase 2. Sci Rep. 2016;6:35277.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  164. Kayacelebi AA, Minović I, Hanff E, Frenay ARS, de Borst MH, Feelisch M, et al. Low plasma homoarginine concentration is associated with high rates of all-cause mortality in renal transplant recipients. Amino Acids. 2017;49:1193–202.

    Article  CAS  PubMed  Google Scholar 

  165. Ravani P, Maas R, Malberti F, Pecchini P, Mieth M, Quinn R, et al. Homoarginine and mortality in pre-dialysis chronic kidney disease (CKD) patients. PLoS One. 2013;8.

  166. Drechsler C, Kollerits B, Meinitzer A, März W, Ritz E, König P, et al. Homoarginine and progression of chronic kidney disease: results from the Mild to Moderate Kidney Disease Study. PLoS One. 2013;8. PMID: 24023762.

  167. Diao JA, Inker LA, Levey AS, Tighiouart H, Powe NR, Manrai AK. In Search of a Better Equation - Performance and Equity in Estimates of Kidney Function. N Engl J Med. 2021;384:396–9.

    Article  PubMed  PubMed Central  Google Scholar 

  168. Doshi MD, Singh N, Hippen BE, Woodside KJ, Mohan P, Byford HL, et al. Transplant Clinician Opinions on Use of Race in the Estimation of Glomerular Filtration Rate. Clin J Am Soc Nephrol. 2021;16:1552–9.

    Article  PubMed  PubMed Central  Google Scholar 

  169. Bailie C, Kilner J, Maxwell AP, McKnight AJ. Development of next generation sequencing panel for UMOD and association with kidney disease. PLoS ONE. 2017;12: e0178321.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

This research has been conducted using the UK Biobank Resourceunder Application Number14259. Some elements of this research were included in PhD theses submitted by RC & RS at Queen´s University Belfast, UK, and Master’s theses by BBJ and JJM at University of Granada, Spain.

The authors would like to thank Prof. Chris Patterson for his statistical advice on the analysis of this work.

The Genotype-Tissue Expression (GTEx) Project was supported by the Common Fund of the Office of the Director of the National Institutes of Health, and by NCI, NHGRI, NHLBI, NIDA, NIMH, and NINDS. The data used for the analyses described in this manuscript were obtained from the GTEx Portal on 31/01/2022 (GTEx Analysis V8 release; dbGaP accession number phs000424.v8.p2)

Funding

Funding for the development of this work was provided by Science Foundation Ireland and the Department for the Economy, Northern Ireland Investigator Program Partnership Award (15/IA/3152). MCG and CH are supported by 15/IA/3152. LJS is supported by funding from the National Institutes of Health (R01AG068937), Northern Ireland Health and Social Care Research and Development Office (STL/5569/19), and the Medical Research Council (MC_PC_20026). RS & RC were in receipt of PhD studentships from the NI Department for the Economy. BBJ and JJM were supported by the Erasmus+ EU programme for education, training, youth and sport for master students.

Author information

Author notes

  1. Blanca Baños-Jaime and Joaquín J. Maqueda contributed equally to this work.

Authors and Affiliations

  1. Molecular Epidemiology and Public Health Research Group, Centre for Public Health,, Queen’s University Belfast, Institute for Clinical Sciences A, Royal Victoria Hospital, Belfast, BT12 6BA, UK

    Marisa Cañadas-Garre, Blanca Baños-Jaime, Joaquín J. Maqueda, Laura J. Smyth, Ruaidhri Cappa, Ryan Skelly, Claire Hill, Alexander P. Maxwell & Amy Jayne McKnight

  2. Genomic Oncology Area, Centre for Genomics and Oncological Research: Pfizer, GENYO, University of Granada-Andalusian Regional Government, PTS Granada. Avenida de La Ilustración 114, 18016, Granada, Spain

    Marisa Cañadas-Garre

  3. Hematology Department, Hospital Universitario Virgen de Las Nieves, Avenida de Las Fuerzas Armadas 2, 18014, Granada, Spain

    Marisa Cañadas-Garre

  4. Instituto de Investigación Biosanitaria de Granada (Ibs.GRANADA), Avda. de Madrid, 15, 18012, Granada, Spain

    Marisa Cañadas-Garre

  5. Instituto de Investigaciones Químicas (IIQ), Centro de Investigaciones Científicas Isla de La Cartuja (cicCartuja), Consejo Superior de Investigaciones Científicas (CSIC), Universidad de Sevilla, Avda. Américo Vespucio 49, 41092, Seville, Spain

    Blanca Baños-Jaime

  6. Experimental Oncology Laboratory, IRCCS Rizzoli Orthopaedic Institute, 40136, Bologna, Italy

    Joaquín J. Maqueda

  7. Department of Experimental, Diagnostic and Specialty Medicine (DIMES), University of Bologna, 40126, Bologna, Italy

    Joaquín J. Maqueda

  8. UCD Diabetes Complications Research Centre, Conway Institute of Biomolecular and Biomedical Research, University College Dublin, Dublin, D04 V1W8, Ireland

    Eoin P. Brennan, Ross Doyle & Catherine Godson

  9. School of Medicine, University College Dublin, Dublin, D04 V1W8, Ireland

    Eoin P. Brennan, Ross Doyle & Catherine Godson

  10. Mater Misericordiae University Hospital, Eccles St, Dublin, D07 R2WY, Ireland

    Ross Doyle

  11. Regional Nephrology Unit, Belfast City Hospital, Level 11Lisburn Road, Belfast, BT9 7AB, UK

    Alexander P. Maxwell

Authors
  1. Marisa Cañadas-Garre
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Blanca Baños-Jaime
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Joaquín J. Maqueda
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Laura J. Smyth
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Ruaidhri Cappa
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Ryan Skelly
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Claire Hill
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Eoin P. Brennan
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Ross Doyle
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Catherine Godson
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Alexander P. Maxwell
    View author publications

    You can also search for this author in PubMed Google Scholar

  12. Amy Jayne McKnight
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

MCG, AJM and APM conceived and designed the study. MCG performed the statistical analyses, interpreted the data and drafted the manuscript. MCG, BBJ, JJM, LJS, RC, RS, EPB, RD, CG, APM and AJM provided important intellectual content, contributed to data interpretation, critically revised, and approved the final version of the manuscript.

Corresponding author

Correspondence to Marisa Cañadas-Garre.

Ethics declarations

Ethics approval and consent to participate

This research was developed in accordance with the Declaration of Helsinki. The work was partly conducted using the UK Biobank Resource under Application Number 14259, involving secondary use of publicly available data, so did not require separate ethical approval. Similarly, data analysed for the UK-ROI collection was also novel secondary use of publicly available data. Both collections employed informed consent for participants with relevant ethical approvals for the primary studies [73, 169]. Ethical approval from Queens University Belfast and the Office of Research Ethics Northern Ireland was obtained for the UK-ROI study [169]. Ethics approval for the primary UK Biobank study was obtained from the North West Centre for Research Ethics Committee (11/NW/0382) [73].

Consent for publication

Not applicable. This manuscript does not contain any individual person’s data in any form.

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

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

Cañadas-Garre, M., Baños-Jaime, B., Maqueda, J.J. et al. Genetic variants affecting mitochondrial function provide further insights for kidney disease. BMC Genomics 25, 576 (2024). https://doi.org/10.1186/s12864-024-10449-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-024-10449-1

Keywords

BMC Genomics

ISSN: 1471-2164

Contact us