Genome-wide association meta-analyses and fine-mapping elucidate pathways influencing albuminuria

Systematic Review of Genetic Modifiers Associated with the Development and/or Progression of Nephropathy in Patients with Sickle Cell Disease

Sickle cell nephropathy (SCN) is a common complication of sickle cell disease (SCD) that significantly contributes to morbidity and mortality. In addition to clinical and life-style factors, genetic variants influence this risk. We performed a systematic review, searching five databases. Studies evaluating the effect of genetic modifiers on SCN were eligible. Twenty-eight studies (fair-to-good quality) were included: one genome-wide association study, twenty-six case-control studies, and one article combining both approaches. APOL1 was significantly associated with albuminuria and hyperfiltration in children and with worse glomerular filtration in adults. On the other hand, alpha-thalassemia protected patients against albuminuria and hyperfiltration, while BCL11A variants were protective against albuminuria alone. The HMOX1 long GT-tandem repeat polymorphism led to a lower glomerular filtration rate. No modifiers for the risk of hyposthenuria were identified. A genome-wide association approach identified three new loci for proteinuria (CRYL1, VWF, and ADAMTS7) and nine loci were linked with eGFR (PKD1L2, TOR2A, CUBN, AGGF1, CYP4B1, CD163, LRP1B, linc02288, and FPGT-TNNI3K/TNNI3K). In conclusion, this systematic review supports the role of genetic modifiers in influencing the risk and progression of SCN. Incorporating and expanding this knowledge is crucial to improving the management and clinical outcomes of patients at risk.

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

SUMMARY

This article introduces the metaGWASmanager, which streamlines genome-wide association studies within large-scale meta-analysis consortia. It is a toolbox for both the central consortium analysis group and participating studies to generate homogeneous phenotypes, minimize unwanted variability from inconsistent methodologies, ensure high-quality association results, and implement time-efficient quality control workflows. The toolbox features a plug-in-based approach for customization of association testing.

RESULTS

The metaGWASmanager toolbox has been successfully deployed in both the CKDGen and MetalGWAS Initiative consortia across hundreds of participating studies, demonstrating its effectiveness in GWAS analysis optimization by automating routine tasks and ensuring the value and reliability of association results, thus, ultimately promoting scientific discovery. We provide a simulated data set with examples for script customization so that readers can reproduce the pipeline at their convenience.

AVAILABILITY AND IMPLEMENTATION

GitHub: https://github.com/genepi-freiburg/metaGWASmanager.

A clustering approach to improve our understanding of the genetic and phenotypic complexity of chronic kidney disease

Chronic kidney disease (CKD) is a complex disorder that causes a gradual loss of kidney function, affecting approximately 9.1% of the world's population. Here, we use a soft-clustering algorithm to deconstruct its genetic heterogeneity. First, we selected 322 CKD-associated independent genetic variants from published genome-wide association studies (GWAS) and added association results for 229 traits from the GWAS catalog. We then applied nonnegative matrix factorization (NMF) to discover overlapping clusters of related traits and variants. We computed cluster-specific polygenic scores and validated each cluster with a phenome-wide association study (PheWAS) on the BioMe biobank (n = 31,701). NMF identified nine clusters that reflect different aspects of CKD, with the top-weighted traits signifying areas such as kidney function, type 2 diabetes (T2D), and body weight. For most clusters, the top-weighted traits were confirmed in the PheWAS analysis. Results were found to be more significant in the cross-ancestry analysis, although significant ancestry-specific associations were also identified. While all alleles were associated with a decreased kidney function, associations with CKD-related diseases (e.g., T2D) were found only for a smaller subset of variants and differed across genetic ancestry groups. Our findings leverage genetics to gain insights into the underlying biology of CKD and investigate population-specific associations.

Genetic Association and Transferability for Urinary Albumin-Creatinine Ratio as a Marker of Kidney Disease in four Sub-Saharan African Populations and non-continental Individuals of African Ancestry

BACKGROUND

Genome-wide association studies (GWAS) have predominantly focused on populations of European and Asian ancestry, limiting our understanding of genetic factors influencing kidney disease in Sub-Saharan African (SSA) populations. This study presents the largest GWAS for urinary albumin-to-creatinine ratio (UACR) in SSA individuals, including 8,970 participants living in different African regions and an additional 9,705 non-resident individuals of African ancestry from the UK Biobank and African American cohorts.

METHODS

Urine biomarkers and genotype data were obtained from two SSA cohorts (AWI-Gen and ARK), and two non-resident African-ancestry studies (UK Biobank and CKD-Gen Consortium). Association testing and meta-analyses were conducted, with subsequent fine-mapping, conditional analyses, and replication studies. Polygenic scores (PGS) were assessed for transferability across populations.

RESULTS

Two genome-wide significant (P<5x10-8) UACR-associated loci were identified, one in the BMP6 region on chromosome 6, in the meta-analysis of resident African individuals, and another in the HBB region on chromosome 11 in the meta-analysis of non-resident SSA individuals, as well as the combined meta-analysis of all studies. Replication of previous significant results confirmed associations in known UACR-associated regions, including THB53, GATM, and ARL15. PGS estimated using previous studies from European ancestry, African ancestry, and multi-ancestry cohorts exhibited limited transferability of PGS across populations, with less than 1% of observed variance explained.

CONCLUSION

This study contributes novel insights into the genetic architecture of kidney disease in SSA populations, emphasizing the need for conducting genetic research in diverse cohorts. The identified loci provide a foundation for future investigations into the genetic susceptibility to chronic kidney disease in underrepresented African populations Additionally, there is a need to develop integrated scores using multi-omics data and risk factors specific to the African context to improve the accuracy of predicting disease outcomes.

METHODS

Urine biomarkers and genotype data were obtained from two SSA cohorts (AWI-Gen and ARK), and two non-resident African-ancestry studies (UK Biobank and CKD-Gen Consortium). Association testing and meta-analyses were conducted, with subsequent fine-mapping, conditional analyses, and replication studies. Polygenic scores (PGS) were assessed for transferability across populations.

RESULTS

Two genome-wide significant (P<5x10-8) UACR-associated loci were identified, one in the BMP6 region on chromosome 6, in the meta-analysis of resident African individuals, and another in the HBB region on chromosome 11 in the meta-analysis of non-resident SSA individuals, as well as the combined meta-analysis of all studies. Replication of previous significant results confirmed associations in known UACR-associated regions, including THB53, GATM, and ARL15. PGS estimated using previous studies from European ancestry, African ancestry, and multi-ancestry cohorts exhibited limited transferability of PGS across populations, with less than 1% of observed variance explained.

CONCLUSION

This study contributes novel insights into the genetic architecture of kidney function in SSA populations, emphasizing the need for conducting genetic research in diverse cohorts. The identified loci provide a foundation for future investigations into the genetic susceptibility to chronic kidney disease in underrepresented African populations.

A multi-ancestry GWAS of Fuchs corneal dystrophy highlights the contributions of laminins, collagen, and endothelial cell regulation

Fuchs endothelial corneal dystrophy (FECD) is a leading indication for corneal transplantation, but its molecular etiology remains poorly understood. We performed genome-wide association studies (GWAS) of FECD in the Million Veteran Program followed by multi-ancestry meta-analysis with the previous largest FECD GWAS, for a total of 3970 cases and 333,794 controls. We confirm the previous four loci, and identify eight novel loci: SSBP3, THSD7A, LAMB1, PIDD1, RORA, HS3ST3B1, LAMA5, and COL18A1. We further confirm the TCF4 locus in GWAS for admixed African and Hispanic/Latino ancestries and show an enrichment of European-ancestry haplotypes at TCF4 in FECD cases. Among the novel associations are low frequency missense variants in laminin genes LAMA5 and LAMB1 which, together with previously reported LAMC1, form laminin-511 (LM511). AlphaFold 2 protein modeling, validated through homology, suggests that mutations at LAMA5 and LAMB1 may destabilize LM511 by altering inter-domain interactions or extracellular matrix binding. Finally, phenome-wide association scans and colocalization analyses suggest that the TCF4 CTG18.1 trinucleotide repeat expansion leads to dysregulation of ion transport in the corneal endothelium and has pleiotropic effects on renal function.

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

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

SGLT2 inhibition, high-density lipoprotein, and kidney function: a mendelian randomization study

BACKGROUND

Sodium-glucose cotransporter 2 (SGLT2) inhibition is recognized for its evident renoprotective benefits in diabetic renal disease. Recent data suggest that SGLT2 inhibition also slows down kidney disease progression and reduces the risk of acute kidney injury, regardless of whether the patient has diabetes or not, but the mechanism behind these observed effects remains elusive. The objective of this study is to utilize a mendelian randomization (MR) methodology to comprehensively examine the influence of metabolites in circulation regarding the impact of SGLT2 inhibition on kidney function.

METHODS

We used a MR study to obtain associations between genetic proxies for SGLT2 inhibition and kidney function. We retrieved the most recent and comprehensive summary statistics from genome-wide association studies (GWAS) that have been previously published and involved kidney function parameters such as estimated glomerular filtration rate (eGFR), urine albumin-to-creatinine ratio (UACR), and albuminuria. Additionally, we included blood metabolite data from 249 biomarkers in the UK Biobank for a more comprehensive analysis. We performed MR analyses to explore the causal relationships between SGLT2 inhibition and kidney function and two-step MR to discover potential mediating metabolites.

RESULTS

The study found that a decrease in HbA1c levels by one standard deviation, which is genetically expected to result in SGLT2 inhibition, was linked to a decreased likelihood of developing type 2 diabetes mellitus (T2DM) (odds ratio [OR] = 0.55 [95% CI 0.35, 0.85], P = 0.007). Meanwhile, SGLT2 inhibition also protects eGFR (β = 0.05 [95% CI 0.03, 0.08], P = 2.45 × 10- 5) and decreased UACR (-0.18 [95% CI -0.33, -0.02], P = 0.025) and albuminuria (-1.07 [95% CI -1.58, -0.57], P = 3.60 × 10- 5). Furthermore, the study found that of the 249 metabolites present in the blood, only one metabolite, specifically the concentration of small high-density lipoprotein (HDL) particles, was significantly correlated with both SGLT2 inhibition and kidney function. This metabolite was found to play a crucial role in mediating the improvement of renal function through the use of SGLT2 inhibition (β = 0.01 [95% CI 0.005, 0.018], P = 0.001), with a mediated proportion of 13.33% (95% CI [5.71%, 26.67%], P = 0.020).

CONCLUSIONS

The findings of this investigation provide evidence in favor of a genetically anticipated biological linkage between the inhibition of SGLT2, the presence of circulating metabolites, and renal function. The findings demonstrate that the protective effect of SGLT2 inhibition on renal function is mostly mediated by HDL particle concentrations in circulating metabolites. These results offer significant theoretical support for both the preservation of renal function and a better comprehension of the mechanisms underlying SGLT2 inhibition.

Effects of Metabolic Syndrome and Its Components on Chronic Kidney Disease and Renal Function: A Two-Sample Mendelian Randomization Study

Objective: The association of metabolic syndrome (MetS) and its components with chronic kidney disease (CKD) and renal function remains controversial in observational studies. To comprehensively investigate the association between MetS and its components with CKD and renal function, a Mendelian randomization (MR) study was performed. Methods: The inverse variance weighting (IVW) of random effects was used as the main estimation method, while MR-Egger and weighted median analysis results were used for auxiliary judgments. Cochran's Q test, MR-Egger intercept test, leave-one-out analysis, and funnel plots were used to assess heterogeneity and pleiotropy. Results: The MR analyses of genetically predicted MetS and its components' association with CKD risk and renal function showed the following causal associations: hypertension with CKD risk; MetS and obesity with increased blood urea nitrogen and decreased estimated glomerular filtration rate based on cystatin C; hypertension and diabetes with increased urine albumin-creatinine ratio and increased risk of microalbuminuria; and CKD with increased triglyceride. Conclusion: Based on genetic data, this study demonstrated an association between hypertension and CKD risk and a causal association between other MetS components and renal function. The early diagnosis and prevention of MetS and its components might be essential for CKD management.

Novel genetic markers for chronic kidney disease in a geographically isolated population of Indigenous Australians: Individual and multiple phenotype genome-wide association study

BACKGROUND

Chronic kidney disease (CKD) is highly prevalent among Indigenous Australians, especially those in remote regions. The Tiwi population has been isolated from mainland Australia for millennia and exhibits unique genetic characteristics that distinguish them from other Indigenous and non-Indigenous populations. Notably, the rate of end-stage renal disease is up to 20 times greater in this population compared to non-Indigenous populations. Despite the identification of numerous genetic loci associated with kidney disease through GWAS, the Indigenous population such as Tiwi remains severely underrepresented and the increased prevalence of CKD in this population may be due to unique disease-causing alleles/genes.

METHODS

We used albumin-to-creatinine ratio (ACR) and estimated glomerular filtration rate (eGFR) to estimate the prevalence of kidney disease in the Tiwi population (N = 492) in comparison to the UK Biobank (UKBB) (N = 134,724) database. We then performed an exploratory factor analysis to identify correlations among 10 CKD-related phenotypes and identify new multi-phenotype factors. We subsequently conducted a genome-wide association study (GWAS) on all single and multiple phenotype factors using mixed linear regression models, adjusted for age, sex, population stratification, and genetic relatedness between individuals.

RESULTS

Based on ACR, 20.3% of the population was at severely increased risk of CKD progression and showed elevated levels of ACR compared to the UKBB population independent of HbA1c. A GWAS of ACR revealed novel association loci in the genes MEG3 (chr14:100812018:T:A), RAB36 (rs11704318), and TIAM2 (rs9689640). Additionally, multiple phenotypes GWAS of ACR, eGFR, urine albumin, and serum creatinine identified a novel variant that mapped to the gene MEIS2 (chr15:37218869:A:G). Most of the identified variants were found to be either absent or rare in the UKBB population.

CONCLUSIONS

Our study highlights the Tiwi population's predisposition towards elevated ACR, and the collection of novel genetic variants associated with kidney function. These associations may prove valuable in the early diagnosis and treatment of renal disease in this underrepresented population. Additionally, further research is needed to comprehensively validate the functions of the identified variants/genes.

Multi-omic analysis of human kidney tissue identified medulla-specific gene expression patterns

The kidney medulla is a specialized region with important homeostatic functions. It has been implicated in genetic and developmental disorders along with ischemic and drug-induced injuries. Despite its role in kidney function and disease, the medulla's baseline gene expression and epigenomic signatures have not been well described in the adult human kidney. Here we generated and analyzed gene expression (RNA-seq), chromatin accessibility (ATAC-seq), chromatin conformation (Hi-C) and spatial transcriptomic data from the adult human kidney cortex and medulla. Tissue samples were obtained from macroscopically dissected cortex and medulla of tumor-adjacent normal material in nephrectomy specimens from five male patients. We used these carefully annotated specimens to reassign incorrectly labeled samples in the larger public Genotype-Tissue Expression (GTEx) Project, and to extract meaningful medullary gene expression signatures. Using integrated analysis of gene expression, chromatin accessibility and conformation profiles, we found insights into medulla development and function and then validated this by spatial transcriptomics and immunohistochemistry. Thus, our datasets provide a valuable resource for functional annotation of variants from genome-wide association studies and are freely accessible through an epigenome browser portal.

Unraveling the epigenetic code: human kidney DNA methylation and chromatin dynamics in renal disease development

Epigenetic changes may fill a critical gap in our understanding of kidney disease development, as they not only reflect metabolic changes but are also preserved and transmitted during cell division. We conducted a genome-wide cytosine methylation analysis of 399 human kidney samples, along with single-nuclear open chromatin analysis on over 60,000 cells from 14 subjects, including controls, and diabetes and hypertension attributed chronic kidney disease (CKD) patients. We identified and validated differentially methylated positions associated with disease states, and discovered that nearly 30% of these alterations were influenced by underlying genetic variations, including variants known to be associated with kidney disease in genome-wide association studies. We also identified regions showing both methylation and open chromatin changes. These changes in methylation and open chromatin significantly associated gene expression changes, most notably those playing role in metabolism and expressed in proximal tubules. Our study further demonstrated that methylation risk scores (MRS) can improve disease state annotation and prediction of kidney disease development. Collectively, our results suggest a causal relationship between epigenetic changes and kidney disease pathogenesis, thereby providing potential pathways for the development of novel risk stratification methods.

Association between Psoriasis and Renal Functions: An Integration Study of Observational Study and Mendelian Randomization

Psoriasis is an autoimmune-mediated disease with several comorbidities in addition to typical skin lesions. Increasing evidence shows the relationships between psoriasis and renal functions, but the relationship and causality remain unclear. We aimed to investigate the associations and causality between psoriasis and four renal functions, including the estimated glomerular filtration rate (eGFR), blood urea nitrogen (BUN), urine albumin to creatinine ratio (UACR), and chronic kidney disease (CKD). For the population-based study, we analyzed the National Health and Nutrition Examination Survey (NHANES) data from five cycles (2003-2006 and 2009-2014) on psoriasis and renal functions. Subgroup analyses were conducted among different categories of participants. Meanwhile, a bidirectional two-sample Mendelian randomization (TSMR) study in European populations was also performed using summary-level genetic datasets. Causal effects were derived by conducting an inverse-variance weighted (MR-IVW) method. A series of pleiotropy-robust MR methods was employed to validate the robustness. Multivariable MR (MVMR) was conducted to complement the result when five competing risk factors were considered. A total of 20,244 participants were enrolled in the cross-sectional study, where 2.6% of them had psoriasis. In the fully adjusted model, participants with psoriasis had significantly lower eGFR (p = 0.025) compared with the healthy group. Individuals who are nonoverweight are more likely to be affected by psoriasis, leading to an elevation of BUN (Pint = 0.018). In the same line, TSMR showed a negative association between psoriasis and eGFR (p = 0.016), and sensitive analysis also consolidated the finding. No causality was identified between psoriasis and other renal functions, as well as the inverse causality (p > 0.05). The MVMR method further provided quite consistent results when adjusting five confounders (p = 0.042). We detected a significant negative effect of psoriasis on eGFR, with marginal association between BUN, UACR, and CKD. The adverse of psoriasis on the renal should merit further attention in clinical cares.

X-chromosome and kidney function: evidence from a multi-trait genetic analysis of 908,697 individuals reveals sex-specific and sex-differential findings in genes regulated by androgen response elements

X-chromosomal genetic variants are understudied but can yield valuable insights into sexually dimorphic human traits and diseases. We performed a sex-stratified cross-ancestry X-chromosome-wide association meta-analysis of seven kidney-related traits (n = 908,697), identifying 23 loci genome-wide significantly associated with two of the traits: 7 for uric acid and 16 for estimated glomerular filtration rate (eGFR), including four novel eGFR loci containing the functionally plausible prioritized genes ACSL4, CLDN2, TSPAN6 and the female-specific DRP2. Further, we identified five novel sex-interactions, comprising male-specific effects at FAM9B and AR/EDA2R, and three sex-differential findings with larger genetic effect sizes in males at DCAF12L1 and MST4 and larger effect sizes in females at HPRT1. All prioritized genes in loci showing significant sex-interactions were located next to androgen response elements (ARE). Five ARE genes showed sex-differential expressions. This study contributes new insights into sex-dimorphisms of kidney traits along with new prioritized gene targets for further molecular research.

Mendelian randomization indicates causal effects of estradiol levels on kidney function in males

Context

Chronic kidney disease (CKD) is a public health burden worldwide. Epidemiological studies observed an association between sex hormones, including estradiol, and kidney function.

Objective

We conducted a Mendelian randomization (MR) study to assess a possible causal effect of estradiol levels on kidney function in males and females.

Design

We performed a bidirectional two-sample MR using published genetic associations of serum levels of estradiol in men (n = 206,927) and women (n = 229,966), and of kidney traits represented by estimated glomerular filtration rate (eGFR, n = 567,460), urine albumin-to-creatinine ratio (UACR, n = 547,361), and CKD (n = 41,395 cases and n = 439,303 controls) using data obtained from the CKDGen Consortium. Additionally, we conducted a genome-wide association study using UK Biobank cohort study data (n = 11,798 men and n = 6,835 women) to identify novel genetic associations with levels of estradiol, and then used these variants as instruments in a one-sample MR.

Results

The two-sample MR indicated that genetically predicted estradiol levels are significantly associated with eGFR in men (beta = 0.077; p = 5.2E-05). We identified a single locus at chromosome 14 associated with estradiol levels in men being significant in the one-sample MR on eGFR (beta = 0.199; p = 0.017). We revealed significant results with eGFR in postmenopausal women and with UACR in premenopausal women, which did not reach statistical significance in the sensitivity MR analyses. No causal effect of eGFR or UACR on estradiol levels was found.

Conclusions

We conclude that serum estradiol levels may have a causal effect on kidney function. Our MR results provide starting points for studies to develop therapeutic strategies to reduce kidney disease.

The citrate transporter SLC13A5 as a therapeutic target for kidney disease: evidence from Mendelian randomization to inform drug development

BACKGROUND

Solute carrier family 13 member 5 (SLC13A5) is a Na+-coupled citrate co-transporter that mediates entry of extracellular citrate into the cytosol. SLC13A5 inhibition has been proposed as a target for reducing progression of kidney disease. The aim of this study was to leverage the Mendelian randomization paradigm to gain insight into the effects of SLC13A5 inhibition in humans, towards prioritizing and informing clinical development efforts.

METHODS

The primary Mendelian randomization analyses investigated the effect of SLC13A5 inhibition on measures of kidney function, including creatinine and cystatin C-based measures of estimated glomerular filtration rate (creatinine-eGFR and cystatin C-eGFR), blood urea nitrogen (BUN), urine albumin-creatinine ratio (uACR), and risk of chronic kidney disease and microalbuminuria. Secondary analyses included a paired plasma and urine metabolome-wide association study, investigation of secondary traits related to SLC13A5 biology, a phenome-wide association study (PheWAS), and a proteome-wide association study. All analyses were compared to the effect of genetically predicted plasma citrate levels using variants selected from across the genome, and statistical sensitivity analyses robust to the inclusion of pleiotropic variants were also performed. Data were obtained from large-scale genetic consortia and biobanks, with sample sizes ranging from 5023 to 1,320,016 individuals.

RESULTS

We found evidence of associations between genetically proxied SLC13A5 inhibition and higher creatinine-eGFR (p = 0.002), cystatin C-eGFR (p = 0.005), and lower BUN (p = 3 × 10-4). Statistical sensitivity analyses robust to the inclusion of pleiotropic variants suggested that these effects may be a consequence of higher plasma citrate levels. There was no strong evidence of associations of genetically proxied SLC13A5 inhibition with uACR or risk of CKD or microalbuminuria. Secondary analyses identified evidence of associations with higher plasma calcium levels (p = 6 × 10-13) and lower fasting glucose (p = 0.02). PheWAS did not identify any safety concerns.

CONCLUSIONS

This Mendelian randomization analysis provides human-centric insight to guide clinical development of an SLC13A5 inhibitor. We identify plasma calcium and citrate as biologically plausible biomarkers of target engagement, and plasma citrate as a potential biomarker of mechanism of action. Our human genetic evidence corroborates evidence from various animal models to support effects of SLC13A5 inhibition on improving kidney function.

Causal relationship from heart failure to kidney function and CKD: A bidirectional two-sample mendelian randomization study

BACKGROUND

Heart Failure (HF) is a widespread condition that affects millions of people, and it is caused by issues with the heart and blood vessels. Even though we know hypertension, coronary artery disease, obesity, diabetes, and genetics can increase the risk of HF and Chronic Kidney Disease (CKD), the exact cause of these conditions remains a mystery. To bridge this gap, we adopted Mendelian Randomization (MR), which relies on genetic variants as proxies.

METHODS

We used data from European populations for our Bidirectional Two-Sample MR Study, which included 930,014 controls and 47,309 cases of HF from the HERMES consortium, as well as 736,396 controls and 51,256 cases of CKD. We also employed several MR variations, including MR-Egger, Inverse Variance Weighted (IVW), and Weighted Median Estimator (WME), to guarantee the results were accurate and comprehensive.).

RESULTS

In this study, the MR analysis found that individuals with a genetic predisposition for HF have an elevated risk of CKD. Our study revealed a significant association between the genetic prediction of HF and the risk of CKD, as evidenced by the IVW method [with an odds ratio (OR) of 1.12 (95% CI, 1.03-1.21), p = 0.009] and the WME [with an OR of 1.14 (95% CI, 1.03-1.26), p = 0.008]. This causal relationship remained robust even after conducting MR analysis while adjusting for the effects of diabetes and hypertension, yielding ORs of 1.13 (IVW:95% CI, 1.03-1.23), 1.12 (MR-Egger: 95% CI, 0.85-1.48), and 1.15 (WME:95% CI, 1.04-1.27) (p = 0.008). However, in the reverse analysis aiming to explore CKD and renal function as exposures and HF as the outcome, we did not observe a statistically significant causal link between CKD and HF.

CONCLUSION

Our study demonstrates the significance of HF in CKD progression, thus having meaningful implications for treatment and the potential for discovering new therapies. To better understand the relationship between HF and CKD, we need to conduct research in a variety of populations.

GWAS for the composite traits of hematuria and albuminuria

Our GWAS of hematuria in the UK Biobank identified 6 loci, some of which overlap with loci for albuminuria suggesting pleiotropy. Since clinical syndromes are often defined by combinations of traits, generating a combined phenotype can improve power to detect loci influencing multiple characteristics. Thus the composite trait of hematuria and albuminuria was chosen to enrich for glomerular pathologies. Cases had both hematuria defined by ICD codes and albuminuria defined as uACR > 3 mg/mmol. Controls had neither an ICD code for hematuria nor an uACR > 3 mg/mmol. 2429 cases and 343,509 controls from the UK Biobank were included. eGFR was lower in cases compared to controls, with the exception of the comparison in females using CKD-EPI after age adjustment. Variants at 4 loci met genome-wide significance with the following nearest genes: COL4A4, TRIM27, ETV1 and CUBN. TRIM27 is part of the extended MHC locus. All loci with the exception of ETV1 were replicated in the Geisinger MyCode cohort. The previous GWAS of hematuria reported COL4A3-COL4A4 variants and HLA-B*0801 within MHC, which is in linkage disequilibrium with the TRIM27 variant (D' = 0.59). TRIM27 is highly expressed in the tubules. Additional loci included a coding sequence variant in CUBN (p.Ala2914Val, MAF = 0.014 (A), p = 3.29E-8, OR = 2.09, 95% CI = 1.61-2.72). Overall, GWAS for the composite trait of hematuria and albuminuria identified 4 loci, 2 of which were not previously identified in a GWAS of hematuria.

Caffeine Intake, Plasma Caffeine Level, and Kidney Function: A Mendelian Randomization Study

Caffeine is a psychoactive substance widely consumed worldwide, mainly via sources such as coffee and tea. The effects of caffeine on kidney function remain unclear. We leveraged the genetic variants in the CYP1A2 and AHR genes via the two-sample Mendelian randomization (MR) framework to estimate the association of genetically predicted plasma caffeine and caffeine intake on kidney traits. Genetic association summary statistics on plasma caffeine levels and caffeine intake were taken from genome-wide association study (GWAS) meta-analyses of 9876 and of >47,000 European ancestry individuals, respectively. Genetically predicted plasma caffeine levels were associated with a decrease in estimated glomerular filtration rate (eGFR) measured using either creatinine or cystatin C. In contrast, genetically predicted caffeine intake was associated with an increase in eGFR and a low risk of chronic kidney disease. The discrepancy is likely attributable to faster metabolizers of caffeine consuming more caffeine-containing beverages to achieve the same pharmacological effect. Further research is needed to distinguish whether the observed effects on kidney function are driven by the harmful effects of higher plasma caffeine levels or the protective effects of greater intake of caffeine-containing beverages, particularly given the widespread use of drinks containing caffeine and the increasing burden of kidney disease.

A clustering approach to improve our understanding of the genetic and phenotypic complexity of chronic kidney disease

Chronic kidney disease (CKD) is a complex disorder that causes a gradual loss of kidney function, affecting approximately 9.1% of the world's population. Here, we use a soft-clustering algorithm to deconstruct its genetic heterogeneity. First, we selected 322 CKD-associated independent genetic variants from published genome-wide association studies (GWAS) and added association results for 229 traits from the GWAS catalog. We then applied nonnegative matrix factorization (NMF) to discover overlapping clusters of related traits and variants. We computed cluster-specific polygenic scores and validated each cluster with a phenome-wide association study (PheWAS) on the BioMe biobank (n=31,701). NMF identified nine clusters that reflect different aspects of CKD, with the top-weighted traits signifying areas such as kidney function, type 2 diabetes (T2D), and body weight. For most clusters, the top-weighted traits were confirmed in the PheWAS analysis. Results were found to be more significant in the cross-ancestry analysis, although significant ancestry-specific associations were also identified. While all alleles were associated with a decreased kidney function, associations with CKD-related diseases (e.g., T2D) were found only for a smaller subset of variants and differed across genetic ancestry groups. Our findings leverage genetics to gain insights into the underlying biology of CKD and investigate population-specific associations.

A clustering approach to improve our understanding of the genetic and phenotypic complexity of chronic kidney disease

Chronic kidney disease (CKD) is a complex disorder that causes a gradual loss of kidney function, affecting approximately 9.1% of the world's population. Here, we use a soft-clustering algorithm to deconstruct its genetic heterogeneity. First, we selected 322 CKD-associated independent genetic variants from published genome-wide association studies (GWAS) and added association results for 229 traits from the GWAS catalog. We then applied nonnegative matrix factorization (NMF) to discover overlapping clusters of related traits and variants. We computed cluster-specific polygenic scores and validated each cluster with a phenome-wide association study (PheWAS) on the BioMe biobank (n=31,701). NMF identified nine clusters that reflect different aspects of CKD, with the top-weighted traits signifying areas such as kidney function, type 2 diabetes (T2D), and body weight. For most clusters, the top-weighted traits were confirmed in the PheWAS analysis. Results were found to be more significant in the cross-ancestry analysis, although significant ancestry-specific associations were also identified. While all alleles were associated with a decreased kidney function, associations with CKD-related diseases (e.g., T2D) were found only for a smaller subset of variants and differed across genetic ancestry groups. Our findings leverage genetics to gain insights into the underlying biology of CKD and investigate population-specific associations.

Genomics of Treatable Traits in Asthma

The astounding number of genetic variants revealed in the 15 years of genome-wide association studies of asthma has not kept pace with the goals of translational genomics. Moving asthma diagnosis from a nonspecific umbrella term to specific phenotypes/endotypes and related traits may provide insights into features that may be prevented or alleviated by therapeutical intervention. This review provides an overview of the different asthma endotypes and phenotypes and the genomic findings from asthma studies using patient stratification strategies and asthma-related traits. Asthma genomic research for treatable traits has uncovered novel and previously reported asthma loci, primarily through studies in Europeans. Novel genomic findings for asthma phenotypes and related traits may arise from multi-trait and specific phenotyping strategies in diverse populations.

Genome-wide meta-analysis identifies new candidate genes for sickle cell disease nephropathy

Sickle cell disease nephropathy (SCDN), a common SCD complication, is strongly associated with mortality. Polygenic risk scores calculated from recent transethnic meta-analyses of urinary albumin-to-creatinine ratio and estimated glomerular filtration rate (eGFR) trended toward association with proteinuria and eGFR in SCD but the model fit was poor (R2 < 0.01), suggesting that there are likely unique genetic risk factors for SCDN. Therefore, we performed genome-wide association studies (GWAS) for 2 critical manifestations of SCDN, proteinuria and decreased eGFR, in 2 well-characterized adult SCD cohorts, representing, to the best of our knowledge, the largest SCDN sample to date. Meta-analysis identified 6 genome-wide significant associations (false discovery rate, q ≤ 0.05): 3 for proteinuria (CRYL1, VWF, and ADAMTS7) and 3 for eGFR (LRP1B, linc02288, and FPGT-TNNI3K/TNNI3K). These associations are independent of APOL1 risk and represent novel SCDN loci, many with evidence for regulatory function. Moreover, GWAS SNPs in CRYL1, VWF, ADAMTS7, and linc02288 are associated with gene expression in kidney and pathways important to both renal function and SCD biology, supporting the hypothesis that SCDN pathophysiology is distinct from other forms of kidney disease. Together, these findings provide new targets for functional follow-up that could be tested prospectively and potentially used to identify patients with SCD who are at risk, before onset of kidney dysfunction.

Genome-Wide Association Studies of Coffee Intake in UK/US Participants of European Ancestry Uncover Gene-Cohort Influences

Coffee is one of the most widely consumed beverages. We performed a genome-wide association study (GWAS) of coffee intake in US-based 23andMe participants (N=130,153) and identified 7 significant loci, with many replicating in three multi-ancestral cohorts. We examined genetic correlations and performed a phenome-wide association study across thousands of biomarkers and health and lifestyle traits, then compared our results to the largest available GWAS of coffee intake from UK Biobank (UKB; N=334,659). The results of these two GWAS were highly discrepant. We observed positive genetic correlations between coffee intake and psychiatric illnesses, pain, and gastrointestinal traits in 23andMe that were absent or negative in UKB. Genetic correlations with cognition were negative in 23andMe but positive in UKB. The only consistent observations were positive genetic correlations with substance use and obesity. Our study shows that GWAS in different cohorts could capture cultural differences in the relationship between behavior and genetics.

GWAS of random glucose in 476,326 individuals provide insights into diabetes pathophysiology, complications and treatment stratification

Conventional measurements of fasting and postprandial blood glucose levels investigated in genome-wide association studies (GWAS) cannot capture the effects of DNA variability on 'around the clock' glucoregulatory processes. Here we show that GWAS meta-analysis of glucose measurements under nonstandardized conditions (random glucose (RG)) in 476,326 individuals of diverse ancestries and without diabetes enables locus discovery and innovative pathophysiological observations. We discovered 120 RG loci represented by 150 distinct signals, including 13 with sex-dimorphic effects, two cross-ancestry and seven rare frequency signals. Of these, 44 loci are new for glycemic traits. Regulatory, glycosylation and metagenomic annotations highlight ileum and colon tissues, indicating an underappreciated role of the gastrointestinal tract in controlling blood glucose. Functional follow-up and molecular dynamics simulations of lower frequency coding variants in glucagon-like peptide-1 receptor (GLP1R), a type 2 diabetes treatment target, reveal that optimal selection of GLP-1R agonist therapy will benefit from tailored genetic stratification. We also provide evidence from Mendelian randomization that lung function is modulated by blood glucose and that pulmonary dysfunction is a diabetes complication. Our investigation yields new insights into the biology of glucose regulation, diabetes complications and pathways for treatment stratification.

Hypoxia controls expression of kidney-pathogenic MUC1 variants

The interplay between genetic and environmental factors influences the course of chronic kidney disease (CKD). In this context, genetic alterations in the kidney disease gene MUC1 (Mucin1) predispose to the development of CKD. These variations comprise the polymorphism rs4072037, which alters splicing of MUC1 mRNA, the length of a region with variable number of tandem repeats (VNTR), and rare autosomal-dominant inherited dominant-negative mutations in or 5' to the VNTR that causes autosomal dominant tubulointerstitial kidney disease (ADTKD-MUC1). As hypoxia plays a pivotal role in states of acute and chronic kidney injury, we explored the effects of hypoxia-inducible transcription factors (HIF) on the expression of MUC1 and its pathogenic variants in isolated primary human renal tubular cells. We defined a HIF-binding DNA regulatory element in the promoter-proximal region of MUC1 from which hypoxia or treatment with HIF stabilizers, which were recently approved for an anti-anemic therapy in CKD patients, increased levels of wild-type MUC1 and the disease-associated variants. Thus, application of these compounds might exert unfavorable effects in patients carrying MUC1 risk variants.

Migraine, chronic kidney disease and kidney function: observational and genetic analyses

Epidemiological studies demonstrate an association between migraine and chronic kidney disease (CKD), while the genetic basis underlying the phenotypic association has not been investigated. We aimed to help avoid unnecessary interventions in individuals with migraine through the investigation of phenotypic and genetic relationships underlying migraine, CKD, and kidney function. We first evaluated phenotypic associations using observational data from UK Biobank (N = 255,896). We then investigated genetic relationships leveraging genomic data in European ancestry for migraine (Ncase/Ncontrol = 48,975/540,381), CKD (Ncase/Ncontrol = 41,395/439,303), and two traits of kidney function (estimated glomerular filtration rate [eGFR, N = 567,460] and urinary albumin-to-creatinine ratio [UACR, N = 547,361]). Observational analyses suggested no significant association of migraine with the risk of CKD (HR = 1.13, 95% CI = 0.85-1.50). While we did not find any global genetic correlation in general, we identified four specific genomic regions showing significant for migraine with eGFR. Cross-trait meta-analysis identified one candidate causal variant (rs1047891) underlying migraine, CKD, and kidney function. Transcriptome-wide association study detected 28 shared expression-trait associations between migraine and kidney function. Mendelian randomization analysis suggested no causal effect of migraine on CKD (OR = 1.03, 95% CI = 0.98-1.09; P = 0.28). Despite a putative causal effect of migraine on an increased level of UACR (log-scale-beta = 0.02, 95% CI = 0.01-0.04; P = 1.92 × 10-3), it attenuated to null when accounting for both correlated and uncorrelated pleiotropy. Our work does not find evidence supporting a causal association between migraine and CKD. However, our study highlights significant biological pleiotropy between migraine and kidney function. The value of a migraine prophylactic treatment for reducing future CKD in people with migraine is likely limited.

Kidney damage on fertility and pregnancy: A Mendelian randomization

BACKGROUND

Low fertility and adverse pregnancy outcomes are commonly observed in women with chronic kidney disease (CKD). However, a causal relationship between low fertility and adverse pregnancy outcomes with CKD remains unclear. Besides, whether mild kidney dysfunction can affect fertility and pregnancy still needs exploration. Hence, this study aimed to investigate the causal effect of kidney damage on fertility and pregnancy using Mendelian randomization (MR).

METHODS

We first used two-sample MR to examine the effects of kidney damage on fertility and pregnancy. Next, we introduced the Bayesian model averaging MR analysis to detect major causal relationships and render the results robust. The genetic instruments and outcome data were derived from various large genome-wide association studies.

RESULTS

Adverse pregnancy outcomes: Our analyses supported a suggestive causal effect of CKD and estimated glomerular filtration rate (eGFR) rapid on stillbirth, with CKD having an odds ratio (OR) of 1.020 [95% confidence interval (CI) 1.002 to 1.038] and eGFR rapid having an OR of 1.026 (95% CI 1.004-1.048). We also discovered a suggestive causal effect of eGFR on spontaneous abortion, with an OR of 2.63 (95% CI 1.269 to 5.450). Moreover, increased urinary albumin-to-creatinine ratio (UACR) was regarded as a potential risk factor for pre-eclampsia (OR = 1.936; 95% CI 1.065 to 3.517) and gestational hypertension (OR = 1.700; 95% CI 1.002 to 2.886). Fertility assessment: The results indicated that eGFR and UACR had a suggestive causal relationship with the anti-Müllerian hormone level (eGFR beta: 1.004; UACR beta: 0.405).

CONCLUSIONS

Our study used MR to demonstrate a suggestive causal relationship between kidney damage and fertility and pregnancy. We reported that mild kidney dysfunction might be a risk factor for reduced fertility and adverse pregnancy outcomes. Dynamic renal detection may help preserve fertility and reduce the risk of pregnancy loss.

The Effect of DNA Methylation in the Development and Progression of Chronic Kidney Disease in the General Population: An Epigenome-Wide Association Study Using the Korean Genome and Epidemiology Study Database

BACKGROUND

Although knowledge of the genetic factors influencing kidney disease is increasing, epigenetic profiles, which are associated with chronic kidney disease (CKD), have not been fully elucidated. We sought to identify the DNA methylation status of CpG sites associated with reduced kidney function and examine whether the identified CpG sites are associated with CKD development.

METHOD

We analyzed DNA methylation patterns of 440 participants in the Korean Genome and Epidemiology Study (KoGES) with estimated glomerular filtration rates (eGFRs) ≥ 60 mL/min/1.73 m² at baseline. CKD development was defined as a decrease in the eGFR of <60 at any time during an 8-year follow-up period ("CKD prediction" analysis). In addition, among the 440 participants, 49 participants who underwent a second methylation profiling were assessed for an association between a decline in kidney function and changes in the degree of methylation of CpG sites during the 8 years ("kidney function slope" analysis).

RESULTS

In the CKD prediction analysis, methylation profiles of a total of 403,129 CpG sites were evaluated at baseline in 440 participants, and increased and decreased methylation of 268 and 189 CpG sites, respectively, were significantly correlated with the development of CKD in multivariable logistic regression. During kidney function slope analysis using follow-up methylation profiles of 49 participants, the percent methylation changes in 913 CpG sites showed a linear relationship with the percent change in eGFR during 8 years. During functional enrichment analyses for significant CpG sites found in the CKD prediction and kidney function slope analyses, we found that those CpG sites represented MAPK, PI3K/Akt, and Rap1 pathways. In addition, three CpG sites from three genes, NPHS2, CHCHD4, and AHR, were found to be significant in the CKD prediction analysis and related to a decline in kidney function.

CONCLUSION

It is suggested that DNA methylation on specific genes is associated with the development of CKD and the deterioration of kidney function.

Genome-Wide Association Studies for Albuminuria of Nondiabetic Taiwanese Population

INTRODUCTION

Chronic kidney disease, which is defined by a reduced estimated glomerular filtration rate and albuminuria, imposes a large health burden worldwide. Ethnicity-specific associations are frequently observed in genome-wide association studies (GWAS). This study conducts a GWAS of albuminuria in the nondiabetic population of Taiwan.

METHODS

Nondiabetic individuals aged 30-70 years without a history of cancer were enrolled from the Taiwan Biobank. A total of 6,768 subjects were subjected to a spot urine examination. After quality control using PLINK and imputation using SHAPEIT and IMPUTE2, a total of 3,638,350 single-nucleotide polymorphisms (SNPs) remained for testing. SNPs with a minor allele frequency of less than 0.1% were excluded. Linear regression was used to determine the relationship between SNPs and log urine albumin-to-creatinine ratio.

RESULTS

Six suggestive loci are identified in or near the FCRL3 (p = 2.56 × 10-6), TMEM161 (p = 4.43 × 10-6), EFCAB1 (p = 2.03 × 10-6), ELMOD1 (p = 2.97 × 10-6), RYR3 (p = 1.34 × 10-6), and PIEZO2 (p = 2.19 × 10-7). Genetic variants in the FCRL3 gene that encode a secretory IgA receptor are found to be associated with IgA nephropathy, which can manifest as proteinuria. The PIEZO2 gene encodes a sensor for mechanical forces in mesangial cells and renin-producing cells. Five SNPs with a p-value between 5 × 10-6 and 5 × 10-5 are also identified in five genes that may have a biological role in the development of albuminuria.

CONCLUSION

Five new loci and one known suggestive locus for albuminuria are identified in the nondiabetic Taiwanese population.

Relationship of ultra-processed food consumption and new-onset chronic kidney diseases among participants with or without diabetes

BACKGROUND

Whether diabetes and genetic susceptibility of kidney diseases modifies the relationship between ultra-processed foods (UPF) consumption and incident chronic kidney disease (CKD) remains uncertain. We aimed to investigate the association between UPF consumption and new-onset CKD in participants with and without diabetes, and explore whether genetic risks of kidney diseases may modify the association.

METHODS

153,985 participants who were free of CKD at baseline and provided 24-h dietary recalls in the UK Biobank were included. UPF was defined according to the NOVA classification. The energy contribution of UPF was calculated by dividing the energy intake of UPF by the total energy intake. The study outcome was new-onset CKD, ascertained by self-report data and data linkage with primary care, hospital admissions, and death registry records.

RESULTS

During a median follow-up of 12.1 years, 4,058 participants developed new-onset CKD. There was a significant positive association between UPF consumption and new-onset CKD in total participants (per 10% increment, adjusted hazard ratio (HR) 1.04; 95% confidence interval (CI) [1.01;1.06]. The positive association between UPF consumption and risk of new-onset CKD was significantly stronger in participants with diabetes (per 10% increment, adjusted HR 1.11 [1.05;1.17]) than in those without diabetes (per 10% increment, adjusted HR 1.03 [1.00;1.05]; P-interaction = 0.005). Genetic risks of kidney diseases did not significantly modify the positive association in those with or without diabetes (all P-interactions > 0.05).

CONCLUSION

There was a significantly stronger positive association between UPF consumption and new-onset CKD in participants with diabetes compared with those without diabetes.

Multi-ancestry meta-analysis and fine-mapping in Alzheimer's disease

Genome-wide association studies (GWAS) of Alzheimer's disease are predominantly carried out in European ancestry individuals despite the known variation in genetic architecture and disease prevalence across global populations. We leveraged published GWAS summary statistics from European, East Asian, and African American populations, and an additional GWAS from a Caribbean Hispanic population using previously reported genotype data to perform the largest multi-ancestry GWAS meta-analysis of Alzheimer's disease and related dementias to date. This method allowed us to identify two independent novel disease-associated loci on chromosome 3. We also leveraged diverse haplotype structures to fine-map nine loci with a posterior probability >0.8 and globally assessed the heterogeneity of known risk factors across populations. Additionally, we compared the generalizability of multi-ancestry- and single-ancestry-derived polygenic risk scores in a three-way admixed Colombian population. Our findings highlight the importance of multi-ancestry representation in uncovering and understanding putative factors that contribute to risk of Alzheimer's disease and related dementias.

Transcriptome- and proteome-wide association studies nominate determinants of kidney function and damage

BACKGROUND

The pathophysiological causes of kidney disease are not fully understood. Here we show that the integration of genome-wide genetic, transcriptomic, and proteomic association studies can nominate causal determinants of kidney function and damage.

RESULTS

Through transcriptome-wide association studies (TWAS) in kidney cortex, kidney tubule, liver, and whole blood and proteome-wide association studies (PWAS) in plasma, we assess for effects of 12,893 genes and 1342 proteins on kidney filtration (glomerular filtration rate (GFR) estimated by creatinine; GFR estimated by cystatin C; and blood urea nitrogen) and kidney damage (albuminuria). We find 1561 associations distributed among 260 genomic regions that are supported as putatively causal. We then prioritize 153 of these genomic regions using additional colocalization analyses. Our genome-wide findings are supported by existing knowledge (animal models for MANBA, DACH1, SH3YL1, INHBB), exceed the underlying GWAS signals (28 region-trait combinations without significant GWAS hit), identify independent gene/protein-trait associations within the same genomic region (INHBC, SPRYD4), nominate tissues underlying the associations (tubule expression of NRBP1), and distinguish markers of kidney filtration from those with a role in creatinine and cystatin C metabolism. Furthermore, we follow up on members of the TGF-beta superfamily of proteins and find a prognostic value of INHBC for kidney disease progression even after adjustment for measured glomerular filtration rate (GFR).

CONCLUSION

In summary, this study combines multimodal, genome-wide association studies to generate a catalog of putatively causal target genes and proteins relevant to kidney function and damage which can guide follow-up studies in physiology, basic science, and clinical medicine.

Renoprotective effects of genetically proxied fibroblast growth factor 21: Mendelian randomization, proteome-wide and metabolome-wide association study

BACKGROUND

Fibroblast growth factor 21 (FGF21) has demonstrated efficacy for reducing liver fat and reversing non-alcoholic steatohepatitis in phase 2 clinical trials. It is also postulated to have anti-fibrotic effects and therefore may be amenable to repurposing for the prevention and treatment of chronic kidney disease (CKD).

METHODS

We leverage a missense genetic variant, rs739320 in the FGF21 gene, that associates with magnetic resonance imaging-derived liver fat as a clinically validated and biologically plausible instrumental variable for studying the effects of FGF21 analogs. Performing Mendelian randomization, we ascertain associations between instrumented FGF21 and kidney phenotypes, cardiometabolic disease risk factors, as well as the circulating proteome (Somalogic, 4907 aptamers) and metabolome (Nightingale platform, 249 metabolites).

RESULTS

We report consistent renoprotective associations of genetically proxied FGF21 effect, including higher glomerular filtration rates (p = 1.9 × 10^-4), higher urinary sodium excretion (p = 5.1 × 10^-11), and lower urine albumin-creatinine ratio (p = 3.6 × 10^-5). These favorable effects translated to lower CKD risk (odds ratio per rs739320 C-allele, 0.96; 95%CI, 0.94-0.98; p = 3.2 × 10^-4). Genetically proxied FGF21 effect was also associated with lower fasting insulin, waist-to-hip ratio, blood pressure (systolic and diastolic BP, p < 1.0 × 10^-07) and blood lipid (low-density lipoprotein cholesterol, triglycerides and apolipoprotein B, p < 6.5 × 10^-24) profiles. The latter associations are replicated in our metabolome-wide association study. Proteomic perturbations associated with genetically predicted FGF21 effect were consistent with fibrosis reduction.

CONCLUSION

This study highlights the pleiotropic effects of genetically proxied FGF21 and supports a re-purposing opportunity for the treatment and prevention of kidney disease specifically. Further work is required to triangulate these findings, towards possible clinical development of FGF21 towards the treatment and prevention of kidney disease.

Association Between Chinese Visceral Adipose Index and Albuminuria in Chinese Adults: A Cross-Sectional Study

Purpose

To explore the correlation between Chinese visceral adipose index (CVAI) and urinary microalbumin/creatinine ratio (UACR) and urinary albumin, and whether there is any difference in correlation between Han and Tujia ethnicity.

Methods

This cross-sectional study was conducted in Changde, Hunan, China from May 2021 to December 2021. Biochemical indicators including anthropometric parameters, blood pressure, blood glucose, blood lipids, and UACR of the participants were measured. Univariate analysis, multivariate analyses and multinomial logistic regression analysis were carried out to assess the association between CVAI and albuminuria. In addition, curve fitting and threshold effect analysis were used to explore the nonlinear association between CVAI and albuminuria, and to observe whether there were ethnic differences in this association.

Results

A total of 2026 adult residents were enrolled in this study, 500 of whom had albuminuria. Population-standardized prevalence of albuminuria is 19.06%. In the multivariable model adjusted for confounding factors, the odds ratio (OR) of albuminuria for pre-unit increase of CVAI and pre-SD increase of CVAI were 1.007 (1.003-1.010) and 1.298 (1.127-1.496), respectively. Multinomial logistic regression analysis confirmed the robustness and consistency of the results.The generalized additive model showed that CVAI and albuminuria had a nonlinear relationship with inflection point at 97.201 using the threshold effect. Compared with Han ethnic groups, the threshold between CVAI and albuminuria in Tujia people moved backward. The thresholds were 159.785 and 98.527, respectively.

Conclusion

There was a positive nonlinear dose-response relationship between increased CVAI and higher levels of albuminuria. Maintaining appropriate CVAI levels may be important for the prevention of albuminuria.

Alterations in the Circulating Proteome Associated with Albuminuria

SIGNIFICANCE STATEMENT

We describe circulating proteins associated with albuminuria in a population of African American Study of Kidney Disease and Hypertension with CKD (AASK) using the largest proteomic platform to date: nearly 7000 circulating proteins, representing approximately 2000 new targets. Findings were replicated in a subset of a general population cohort with kidney disease (ARIC) and a population with CKD Chronic Renal Insufficiency Cohort (CRIC). In cross-sectional analysis, 104 proteins were significantly associated with albuminuria in the Black group, of which 67 of 77 available proteins were replicated in ARIC and 68 of 71 available proteins in CRIC. LMAN2, TNFSFR1B, and members of the ephrin superfamily had the strongest associations. Pathway analysis also demonstrated enrichment of ephrin family proteins.

BACKGROUND

Proteomic techniques have facilitated understanding of pathways that mediate decline in GFR. Albuminuria is a key component of CKD diagnosis, staging, and prognosis but has been less studied than GFR. We sought to investigate circulating proteins associated with higher albuminuria.

METHODS

We evaluated the cross-sectional associations of the blood proteome with albuminuria and longitudinally with doubling of albuminuria in the African American Study of Kidney Disease and Hypertension (AASK; 38% female; mean GFR 46; median urine protein-to-creatinine ratio 81 mg/g; n =703) and replicated in two external cohorts: a subset of the Atherosclerosis Risk in Communities (ARIC) study with CKD and the Chronic Renal Insufficiency Cohort (CRIC).

RESULTS

In cross-sectional analysis, 104 proteins were significantly associated with albuminuria in AASK, of which 67 of 77 available proteins were replicated in ARIC and 68 of 71 available proteins in CRIC. Proteins with the strongest associations included LMAN2, TNFSFR1B, and members of the ephrin superfamily. Pathway analysis also demonstrated enrichment of ephrin family proteins. Five proteins were significantly associated with worsening albuminuria in AASK, including LMAN2 and EFNA4, which were replicated in ARIC and CRIC.

CONCLUSIONS

Among individuals with CKD, large-scale proteomic analysis identified known and novel proteins associated with albuminuria and suggested a role for ephrin signaling in albuminuria progression.

Genetic and epigenetic background of diabetic kidney disease

Diabetic kidney disease (DKD) is a severe diabetic complication that affects up to half of the individuals with diabetes. Elevated blood glucose levels are a key underlying cause of DKD, but DKD is a complex multifactorial disease, which takes years to develop. Family studies have shown that inherited factors also contribute to the risk of the disease. During the last decade, genome-wide association studies (GWASs) have emerged as a powerful tool to identify genetic risk factors for DKD. In recent years, the GWASs have acquired larger number of participants, leading to increased statistical power to detect more genetic risk factors. In addition, whole-exome and whole-genome sequencing studies are emerging, aiming to identify rare genetic risk factors for DKD, as well as epigenome-wide association studies, investigating DNA methylation in relation to DKD. This article aims to review the identified genetic and epigenetic risk factors for DKD.

Inflammatory bowel disease increases the levels of albuminuria and the risk of urolithiasis: a two-sample Mendelian randomization study

BACKGROUND

Alterations in kidney function and increased risk of kidney diseases in patients with inflammatory bowel disease (IBD) have been reported, but the causal relationship remains unclear. Herein, Mendelian randomization was employed to identify the causal effect of inflammatory bowel disease on kidney function and the risk of chronic kidney disease (CKD), urolithiasis, and IgA nephropathy.

METHODS

The International Inflammatory Bowel Disease Genetics Consortium provided the summary-level genome-wide association study (GWAS) data that correlates with Crohn's disease (CD) and ulcerative colitis (UC). GWAS data for estimated glomerular filtration rate from serum creatinine (eGFRcrea), urine albumin-creatinine ratio (uACR), and CKD were obtained from the CKDGen Consortium, and GWAS data for urolithiasis were obtained from the FinnGen consortium. The summary-level GWAS data for IgA nephropathy were obtained from the meta-analysis of UK-biobank, FinnGen, and Biobank Japan. Inverse-variance weighted was used as the primary estimate. Furthermore, the Steiger test was used to validate the direction of causality.

RESULTS

The inverse-variance weighted data revealed that genetically predicted UC significantly increased uACR levels, while genetically predicted CD significantly increased the risk of urolithiasis.

CONCLUSIONS

UC increases the levels of uACR, and CD increases the risk of urolithiasis.

Multi-ancestry GWAS of Fuchs corneal dystrophy highlights roles of laminins, collagen, and endothelial cell regulation

Fuchs endothelial corneal dystrophy (FECD) is a leading indication for corneal transplantation, but its molecular pathophysiology remains poorly understood. We performed genome-wide association studies (GWAS) of FECD in the Million Veteran Program (MVP) and meta-analyzed with the previous largest FECD GWAS, finding twelve significant loci (eight novel). We further confirmed the TCF4 locus in admixed African and Hispanic/Latino ancestries, and found an enrichment of European-ancestry haplotypes at TCF4 in FECD cases. Among the novel associations are low frequency missense variants in laminin genes LAMA5 and LAMB1 which, together with previously reported LAMC1, form laminin-511 (LM511). AlphaFold 2 protein modeling suggests that mutations at LAMA5 and LAMB1 may destabilize LM511 by altering inter-domain interactions or extracellular matrix binding. Finally, phenome-wide association scans and co-localization analyses suggest that the TCF4 CTG18.1 trinucleotide repeat expansion leads to dysregulation of ion transport in the corneal endothelium and has pleiotropic effects on renal function.

Mapping genomic regulation of kidney disease and traits through high-resolution and interpretable eQTLs

Expression quantitative trait locus (eQTL) studies illuminate genomic variants that regulate specific genes and contribute to fine-mapped loci discovered via genome-wide association studies (GWAS). Efforts to maximize their accuracy are ongoing. Using 240 glomerular (GLOM) and 311 tubulointerstitial (TUBE) micro-dissected samples from human kidney biopsies, we discovered 5371 GLOM and 9787 TUBE genes with at least one variant significantly associated with expression (eGene) by incorporating kidney single-nucleus open chromatin data and transcription start site distance as an "integrative prior" for Bayesian statistical fine-mapping. The use of an integrative prior resulted in higher resolution eQTLs illustrated by (1) smaller numbers of variants in credible sets with greater confidence, (2) increased enrichment of partitioned heritability for GWAS of two kidney traits, (3) an increased number of variants colocalized with the GWAS loci, and (4) enrichment of computationally predicted functional regulatory variants. A subset of variants and genes were validated experimentally in vitro and using a Drosophila nephrocyte model. More broadly, this study demonstrates that tissue-specific eQTL maps informed by single-nucleus open chromatin data have enhanced utility for diverse downstream analyses.

Association Between Obesity and Chronic Kidney Disease: Multivariable Mendelian Randomization Analysis and Observational Data From a Bariatric Surgery Cohort

Obesity is postulated to independently increase chronic kidney disease (CKD), even after adjusting for type 2 diabetes (T2D) and hypertension. Dysglycemia below T2D thresholds, frequently seen with obesity, also increases CKD risk. Whether obesity increases CKD independent of dysglycemia and hypertension is unknown and likely influences the optimal weight loss (WL) needed to reduce CKD. T2D remission rates plateau with 20-25% WL after bariatric surgery (BS), but further WL increases normoglycemia and normotension. We undertook bidirectional inverse variance weighted Mendelian randomization (IVWMR) to investigate potential independent causal associations between increased BMI and estimated glomerular filtration rate (eGFR) in CKD (CKDeGFR) (<60 mL/min/1.73 m2) and microalbuminuria (MA). In 5,337 BS patients, we assessed whether WL influences >50% decline in eGFR (primary outcome) or CKD hospitalization (secondary outcome), using <20% WL as a comparator. IVWMR results suggest that increased BMI increases CKDeGFR (b = 0.13, P = 1.64 × 10-4; odds ratio [OR] 1.14 [95% CI 1.07, 1.23]) and MA (b = 0.25; P = 2.14 × 10-4; OR 1.29 [1.13, 1.48]). After adjusting for hypertension and fasting glucose, increased BMI did not significantly increase CKDeGFR (b = -0.02; P = 0.72; OR 0.98 [0.87, 1.1]) or MA (b = 0.19; P = 0.08; OR 1.21 [0.98, 1.51]). Post-BS WL significantly reduced the primary outcome with 30 to <40% WL (hazard ratio [HR] 0.53 [95% CI 0.32, 0.87]) but not 20 to <30% WL (HR 0.72 [0.44, 1.2]) and ≥40% WL (HR 0.73 [0.41, 1.30]). For CKD hospitalization, progressive reduction was seen with increased WL, which was significant for 30 to <40% WL (HR 0.37 [0.17, 0.82]) and ≥40% WL (HR 0.24 [0.07, 0.89]) but not 20 to <30% WL (HR 0.60 [0.29, 1.23]). The data suggest that obesity is likely not an independent cause of CKD. WL thresholds previously associated with normotension and normoglycemia, likely causal mediators, may reduce CKD after BS.

A Genomic Risk Score Identifies Individuals at High Risk for Intracerebral Hemorrhage

BACKGROUND

Intracerebral hemorrhage (ICH) has an estimated heritability of 29%. We developed a genomic risk score for ICH and determined its predictive power in comparison to standard clinical risk factors.

METHODS

We combined genome-wide association data from individuals of European ancestry for ICH and related traits in a meta-genomic risk score ([metaGRS]; 2.6 million variants). We tested associations with ICH and its predictive performance in addition to clinical risk factors in a held-out validation dataset (842 cases and 796 controls). We tested associations with risk of incident ICH in the population-based UK Biobank cohort (486 784 individuals, 1526 events, median follow-up 11.3 years).

RESULTS

One SD increment in the metaGRS was significantly associated with 31% higher odds for ICH (95% CI, 1.16-1.48) in age-, sex- and clinical risk factor-adjusted models. The metaGRS identified individuals with almost 5-fold higher odds for ICH in the top score percentile (odds ratio, 4.83 [95% CI, 1.56-21.2]). Predictive models for ICH incorporating the metaGRS in addition to clinical predictors showed superior performance compared to the clinical risk factors alone (c-index, 0.695 versus 0.686). The metaGRS showed similar associations for lobar and nonlobar ICH, independent of the known APOE risk locus for lobar ICH. In the UK Biobank, the metaGRS was associated with higher risk of incident ICH (hazard ratio, 1.15 [95% CI, 1.09-1.21]). The associations were significant within both a relatively high-risk population of antithrombotic medications users, as well as among a relatively low-risk population with a good control of vascular risk factors and no use of anticoagulants.

CONCLUSIONS

We developed and validated a genomic risk score that predicts lifetime risk of ICH beyond established clinical risk factors among individuals of European ancestry. Whether implementation of the score in risk prognostication models for high-risk populations, such as patients under antithrombotic treatment, could improve clinical decision making should be explored in future studies.

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

Background

Congenital anomalies of the kidney and urinary tract (CAKUT) are the predominant cause for chronic kidney disease below 30 years of age. Many monogenic forms have been discovered mainly due to comprehensive genetic testing like exome sequencing (ES). However, disease-causing variants in known disease-associated genes still only explain a proportion of cases. Aim of this study was to unravel the underlying molecular mechanism of syndromic CAKUT in two multiplex families with presumed autosomal recessive inheritance.

Methods and Results

ES in the index individuals revealed two different rare homozygous variants in FOXD2, a transcription factor not previously implicated in CAKUT in humans: a frameshift in family 1 and a missense variant in family 2 with family segregation patterns consistent with autosomal-recessive inheritance. CRISPR/Cas9-derived Foxd2 knock-out (KO) mice presented with bilateral dilated renal pelvis accompanied by renal papilla atrophy while extrarenal features included mandibular, ophthalmologic, and behavioral anomalies, recapitulating the phenotype of humans with FOXD2 dysfunction. To study the pathomechanism of FOXD2-dysfunction-mediated developmental renal defects, in a complementary approach, we generated CRISPR/Cas9-mediated KO of Foxd2 in ureteric-bud-induced mouse metanephric mesenchyme cells. Transcriptomic analyses revealed enrichment of numerous differentially expressed genes important in renal/urogenital development, including Pax2 and Wnt4 as well as gene expression changes indicating a cell identity shift towards a stromal cell identity. Histology of Foxd2 KO mouse kidneys confirmed increased fibrosis. Further, GWAS data (genome-wide association studies) suggests that FOXD2 could play a role for maintenance of podocyte integrity during adulthood.

Conclusions

In summary, our data implicate that FOXD2 dysfunction is a very rare cause of autosomal recessive syndromic CAKUT and suggest disturbances of the PAX2-WNT4 cell signaling axis contribute to this phenotype.

Imputation-powered whole-exome analysis identifies genes associated with kidney function and disease in the UK Biobank

Genome-wide association studies have discovered hundreds of associations between common genotypes and kidney function but cannot comprehensively investigate rare coding variants. Here, we apply a genotype imputation approach to whole exome sequencing data from the UK Biobank to increase sample size from 166,891 to 408,511. We detect 158 rare variants and 105 genes significantly associated with one or more of five kidney function traits, including genes not previously linked to kidney disease in humans. The imputation-powered findings derive support from clinical record-based kidney disease information, such as for a previously unreported splice allele in PKD2, and from functional studies of a previously unreported frameshift allele in CLDN10. This cost-efficient approach boosts statistical power to detect and characterize both known and novel disease susceptibility variants and genes, can be generalized to larger future studies, and generates a comprehensive resource ( https://ckdgen-ukbb.gm.eurac.edu/ ) to direct experimental and clinical studies of kidney disease.

Whole-Exome Sequencing in Family Trios Reveals De Novo Mutations Associated with Type 1 Diabetes Mellitus

Type 1 diabetes mellitus (T1DM) is a chronic autoimmune disease characterized by insulin deficiency and loss of pancreatic islet β-cells. The objective of this study is to identify de novo mutations in 13 trios from singleton families that contribute to the genetic basis of T1DM through the application of whole-exome sequencing (WES). Of the 13 families sampled for this project, 12 had de novo variants, with Family 7 having the highest number (nine) of variants linked to T1DM/autoimmune pathways, whilst Family 4 did not have any variants past the filtering steps. There were 10 variants of 7 genes reportedly associated with T1DM (MST1; TDG; TYRO3; IFIHI; GLIS3; VEGFA; TYK2). There were 20 variants of 13 genes that were linked to endocrine, metabolic, or autoimmune diseases. Our findings demonstrate that trio-based WES is a powerful approach for identifying new candidate genes for the pathogenesis of T1D. Genotyping and functional annotation of the discovered de novo variants in a large cohort is recommended to ascertain their association with disease pathogenesis.

Association of polygenic scores with chronic kidney disease phenotypes in a longitudinal study of older adults

Risk of chronic kidney disease (CKD) is influenced by environmental and genetic factors, and increases sharply in individuals 70 years and older. Polygenic scores (PGS) for kidney disease-related traits have shown promise, but require validation in well-characterized cohorts. Here, we assessed the performance of recently-developed PGSs for CKD-related traits in a longitudinal cohort of healthy older individuals enrolled in the Australian ASPREE randomized controlled trial of daily low-dose aspirin with CKD risk at baseline and longitudinally. Among 11,813 genotyped participants aged 70 years or more with baseline eGFR measures, we tested associations between PGSs and measured eGFR at baseline, clinical phenotype of CKD, and longitudinal rate of eGFR decline spanning up to six years of follow-up per participant. A PGS for eGFR was associated with baseline eGFR, with a significant decrease of 3.9 mL/min/1.73m² (95% confidence interval -4.17 to -3.68) per standard deviation (SD) increase of the PGS. This PGS, as well as a PGS for CKD stage 3 were both associated with higher risk of baseline CKD stage 3 in cross-sectional analysis (Odds Ratio 1.75 per SD , 95% confidence interval 1.66-1.85, and Odds Ratio 1.51 per SD, 95% confidence interval 1.43-1.59, respectively). Longitudinally, two separate PGSs for eGFR slope were associated with significant kidney function decline during follow-up. Thus, our study demonstrates that kidney function has a considerable genetic component in older adults, and that new PGSs for kidney disease-related phenotypes may have potential utility for CKD risk prediction in advanced age.

Allometric body shape indices, type 2 diabetes and kidney function: A two-sample Mendelian randomization study

AIM

To examine the association between body mass index (BMI)-independent allometric body shape indices and kidney function.

MATERIALS AND METHODS

We performed a two-sample Mendelian randomization (MR) analysis, using summary statistics from UK Biobank, CKDGen and DIAGRAM. BMI-independent allometric body shape indices were: A Body Shape Index (ABSI), Waist-Hip Index (WHI) and Hip Index (HI). Kidney function outcomes were: urinary albumin-to-creatinine ratio (UACR), estimated glomerular filtration rate and blood urea nitrogen. Furthermore, we investigated type 2 diabetes (T2D) as a potential mediator on the pathway to albuminuria. The main analysis was inverse variance-weighted random-effects MR in participants of European ancestry. We also performed several sensitivity MR analyses.

RESULTS

A 1-standard deviation (SD) increase in genetically predicted ABSI and WHI levels was associated with higher UACR (β = 0.039 [95% confidence interval: 0.016, 0.063] log [UACR], P = 0.001 for ABSI, and β = 0.028 [0.012, 0.044] log [UACR], P = 6 x 10^-4 for WHI) in women, but not in men. Meanwhile, a 1-SD increase in genetically predicted HI was associated with lower UACR in women (β = -0.021 [-0.041, 0.000] log [UACR], P = 0.05) and in men (β = -0.026 [-0.058, 0.005] log [UACR], P = 0.10). Corresponding estimates in individuals with diabetes were substantially augmented. Risk of T2D increased for genetically high ABSI and WHI in women (P < 6 x 10^-19 ) only, but decreased for genetically high HI in both sexes (P < 9 x 10^-3 ). No other associations were observed.

CONCLUSIONS

Genetically high HI was associated with decreased risk of albuminuria, mediated through decreased T2D risk in both sexes. Opposite associations applied to genetically high ABSI and WHI in women only.

Harnessing Genomic Analysis to Explore the Role of Telomeres in the Pathogenesis and Progression of Diabetic Kidney Disease

The prevalence of diabetes is increasing globally, and this trend is predicted to continue for future decades. Research is needed to uncover new ways to manage diabetes and its co-morbidities. A significant secondary complication of diabetes is kidney disease, which can ultimately result in the need for renal replacement therapy, via dialysis or transplantation. Diabetic kidney disease presents a substantial burden to patients, their families and global healthcare services. This review highlights studies that have harnessed genomic, epigenomic and functional prediction tools to uncover novel genes and pathways associated with DKD that are useful for the identification of therapeutic targets or novel biomarkers for risk stratification. Telomere length regulation is a specific pathway gaining attention recently because of its association with DKD. Researchers are employing both observational and genetics-based studies to identify telomere-related genes associated with kidney function decline in diabetes. Studies have also uncovered novel functions for telomere-related genes beyond the immediate regulation of telomere length, such as transcriptional regulation and inflammation. This review summarises studies that have revealed the potential to harness therapeutics that modulate telomere length, or the associated epigenetic modifications, for the treatment of DKD, to potentially slow renal function decline and reduce the global burden of this disease.

Genome-wide association study of the risk of chronic kidney disease and kidney-related traits in the Japanese population: J-Kidney-Biobank

Chronic kidney disease (CKD) is a syndrome characterized by a gradual loss of kidney function with decreased estimated glomerular filtration rate (eGFR), which may be accompanied by an increase in the urine albumin-to-creatinine ratio (UACR). Although trans-ethnic genome-wide association studies (GWASs) have been conducted for kidney-related traits, there have been few analyses in the Japanese population, especially for the UACR trait. In this study, we conducted a GWAS to identify loci related to multiple kidney-related traits in Japanese individuals. First, to detect loci associated with CKD, eGFR, and UACR, we performed separate GWASs with the following two datasets: 475 cases of CKD diagnosed at seven university hospitals and 3471 healthy subjects (dataset 1) and 3664 cases of CKD-suspected individuals with eGFR <60 ml/min/1.73 m² or urinary protein ≥ 1+ and 5952 healthy subjects (dataset 2). Second, we performed a meta-analysis between these two datasets and detected the following associated loci: 10 loci for CKD, 9 loci for eGFR, and 22 loci for UACR. Among the loci detected, 22 have never been reported previously. Half of the significant loci for CKD were shared with those for eGFR, whereas most of the loci associated with UACR were different from those associated with CKD or eGFR. The GWAS of the Japanese population identified novel genetic components that were not previously detected. The results also suggest that the group primarily characterized by increased UACR possessed genetically different features from the group characterized by decreased eGFR.

Impairment of kidney function and kidney cancer: A bidirectional Mendelian randomization study

BACKGROUND

Many observational epidemiology studies discovered that kidney cancer and impaired kidney function have a bidirectional relationship. However, it remains unclear whether these two kinds of traits are causally linked. In this study, we aimed to investigate the bidirectional causal relation between kidney cancer and kidney function biomarkers (creatinine-based estimated glomerular filtration rate (eGFRcrea), cystatin C-based estimated glomerular filtration rate (eGFRcys), blood urea nitrogen (BUN), serum urate, and urinary albumin-to-creatinine ratio (UACR)).

METHODS

For both directions, single-nucleotide polymorphisms (SNPs), as genetic instruments, for the five kidney function traits were selected from up to 1,004,040 individuals, and SNPs for kidney cancer were from 408,786 participants(1338 cases). In the main analysis, we applied two state-of-the-art MR methods, namely, contamination mixture and Robust Adjusted Profile Score to downweight the effect of weak instrument bias, pleiotropy, and extreme outliers. We additionally conducted traditional MR analyses as sensitivity analyses. Summary-level data of European ancestry were extracted from UK Biobank, Chronic Kidney Disease Genetics Consortium, and Kaiser Permanente.

RESULTS

Based on 99 SNPs, we found that the eGFRcrea had a significant negative causal effect on the risk of kidney cancer (OR = 0.007, 95% CI:2.6 × 10-4 -0.569, p = 0.041). After adjusting for body composition or diabetes, urate had a significant negative causal effect on kidney cancer (OR <1, p < 0.05). For UACR, it showed a strong causal effect on kidney cancer, after adjusting for body composition (OR = 14.503, 95% CI: 2.546-96.001, p = 0.032). Due to lacking significant signals and effect power for the reverse MR, further investigations are warranted.

CONCLUSIONS

Our study suggested a potential causal effect of damaged kidney function on kidney cancer. EGFRcrea and UACR might be causally associated with kidney cancer, especially when patients were comorbid with obesity or diabetes. We called for larger sample-size studies to further unravel the underlying causal relationship and the exact mechanism.

Causal Effects of Specific Gut Microbiota on Chronic Kidney Diseases and Renal Function-A Two-Sample Mendelian Randomization Study

BACKGROUND

Targeting the gut microbiota may become a new therapeutic to prevent and delay the progression of chronic kidney disease (CKD). Nonetheless, the causal relationship between specific intestinal flora and CKD is still unclear.

MATERIALS AND METHOD

To identify genetically predicted microbiota, we used summary data from genome-wide association studies on gut microbiota in 18340 participants from 24 cohorts. Furthermore, we genetically predicted the causal relationship between 211 gut microbiotas and six phenotypes (outcomes) (CKD, estimated glomerular filtration rate (eGFR), urine albumin to creatinine ratio (UACR), dialysis, rapid progress to CKD, and rapid decline of eGFR). Four Mendelian randomization (MR) methods, including inverse variance weighted (IVW), MR-Egger, weighted median, and weighted mode were used to investigate the casual relationship between gut microbiotas and various outcomes. The result of IVW was deemed as the primary result. Then, Cochrane's Q test, MR-Egger, and MR-PRESSO Global test were used to detect heterogeneity and pleiotropy. The leave-one method was used for testing the stability of MR results and Bonferroni-corrected was used to test the strength of the causal relationship between exposure and outcome.

RESULTS

Through the MR analysis of 211 microbiotas and six clinical phenotypes, a total of 36 intestinal microflora were found to be associated with various outcomes. Among them, Class Bacteroidia (=-0.005, 95% CI: -0.001 to -0.008, p = 0.002) has a strong causality with lower eGFR after the Bonferroni-corrected test, whereas phylum Actinobacteria (OR = 1.0009, 95%CI: 1.0003-1.0015, p = 0.0024) has a strong causal relationship with dialysis. The Cochrane's Q test reveals that there is no significant heterogeneity between various single nucleotide polymorphisms. In addition, no significant level of pleiotropy was found according to MR-Egger and MR-PRESSO Global tests.

CONCLUSIONS

Through the two-sample MR analysis, we identified the specific intestinal flora that has a causal relationship with the incidence and progression of CKD at the level of gene prediction, which may provide helpful biomarkers for early disease diagnosis and potential therapeutic targets for CKD.