A catalog of genetic loci associated with kidney function from analyses of a million individuals

Polygenic association of glomerular filtration rate decline in world trade center responders

BACKGROUND

The factors associated with estimated glomerular filtrate rate (eGFR) decline in low risk adults remain relatively unknown. We hypothesized that a polygenic risk score (PRS) will be associated with eGFR decline.

METHODS

We analyzed genetic data from 1,601 adult participants with European ancestry in the World Trade Center Health Program (baseline age 49.68 ± 8.79 years, 93% male, 23% hypertensive, 7% diabetic and 1% with cardiovascular disease) with ≥ three serial measures of serum creatinine. PRSs were calculated from an aggregation of single nucleotide polymorphisms (SNPs) from a recent, large-scale genome-wide association study (GWAS) of rapid eGFR decline. Generalized linear models were used to evaluate the association of PRS with renal outcomes: baseline eGFR and CKD stage, rate of change in eGFR, stable versus declining eGFR over a 3-5-year observation period. eGFR decline was defined in separate analyses as "clinical" (> -1.0 ml/min/1.73 m2/year) or "empirical" (lower most quartile of eGFR slopes).

RESULTS

The mean baseline eGFR was ~ 86 ml/min/1.73 m2. Subjects with decline in eGFR were more likely to be diabetic. PRS was significantly associated with lower baseline eGFR (B = -0.96, p = 0.002), higher CKD stage (OR = 1.17, p = 0.010), decline in eGFR (OR = 1.14, p = 0.036) relative to stable eGFR, and the lower quartile of eGFR slopes (OR = 1.21, p = 0.008), after adjusting for established risk factors for CKD.

CONCLUSION

Common genetic variants are associated with eGFR decline in middle-aged adults with relatively low comorbidity burdens.

Effects of Selenium on Chronic Kidney Disease: A Mendelian Randomization Study

BACKGROUND

Previous observational studies have shown that there is a controversial association between selenium levels and chronic kidney disease (CKD). Our aim was to assess the causal relationship between selenium levels and CKD using Mendelian randomization (MR) analysis.

METHODS

We used the two-sample Mendelian randomization (MR) method to analyze the causal role of selenium levels on CKD risk. The variants associated with selenium levels were extracted from a large genome-wide association study (GWAS) meta-analysis of circulating selenium levels (n = 5477) and toenail selenium levels (n = 4162) in the European population. Outcome data were from the largest GWAS meta-analysis of European-ancestry participants for kidney function to date. Inverse variance weighted (IVW) method was used as the main analysis and a series of sensitivity analyses were carried out to detect potential violations of MR assumptions.

RESULTS

The MR analysis results indicate that the genetically predicted selenium levels were associated with decreased estimated glomerular filtration (eGFR) (effect = -0.0042, 95% confidence interval [CI]: -0.0053-0.0031, p = 2.186 × 10^-13) and increased blood urea nitrogen (BUN) (effect = 0.0029, 95% confidence interval [CI]: 0.0006-0.0052, p = 0.0136) with no pleiotropy detected.

CONCLUSIONS

The MR study indicated that an increased level of selenium is a causative factor for kidney function impairment.

Limited clinical utility for GWAS or polygenic risk score for postoperative acute kidney injury in non-cardiac surgery in European-ancestry patients

BACKGROUND

Prior studies support a genetic basis for postoperative acute kidney injury (AKI). We conducted a genome-wide association study (GWAS), assessed the clinical utility of a polygenic risk score (PRS), and estimated the heritable component of AKI in patients who underwent noncardiac surgery.

METHODS

We performed a retrospective large-scale genome-wide association study followed by a meta-analysis of patients who underwent noncardiac surgery at the Vanderbilt University Medical Center ("Vanderbilt" cohort) or Michigan Medicine, the academic medical center of the University of Michigan ("Michigan" cohort). In the Vanderbilt cohort, the relationship between polygenic risk score for estimated glomerular filtration rate and postoperative AKI was also tested to explore the predictive power of aggregating multiple common genetic variants associated with AKI risk. Similarly, in the Vanderbilt cohort genome-wide complex trait analysis was used to estimate the heritable component of AKI due to common genetic variants.

RESULTS

The study population included 8248 adults in the Vanderbilt cohort (mean [SD] 58.05 [15.23] years, 50.2% men) and 5998 adults in Michigan cohort (56.24 [14.76] years, 49% men). Incident postoperative AKI events occurred in 959 patients (11.6%) and in 277 patients (4.6%), respectively. No loci met genome-wide significance in the GWAS and meta-analysis. PRS for estimated glomerular filtration rate explained a very small percentage of variance in rates of postoperative AKI and was not significantly associated with AKI (odds ratio 1.050 per 1 SD increase in polygenic risk score [95% CI, 0.971-1.134]). The estimated heritability among common variants for AKI was 4.5% (SE = 4.5%) suggesting low heritability.

CONCLUSION

The findings of this study indicate that common genetic variation minimally contributes to postoperative AKI after noncardiac surgery, and likely has little clinical utility for identifying high-risk patients.

Differentiating Associations of Glycemic Traits With Atherosclerotic and Thrombotic Outcomes: Mendelian Randomization Investigation

We conducted a Mendelian randomization analysis to differentiate associations of four glycemic indicators with a broad range of atherosclerotic and thrombotic diseases. Independent genetic variants associated with fasting glucose (FG), 2 h glucose after an oral glucose challenge (2hGlu), fasting insulin (FI), and glycated hemoglobin (HbA1c) at the genome-wide significance threshold were used as instrumental variables. Summary-level data for 12 atherosclerotic and 4 thrombotic outcomes were obtained from large genetic consortia and the FinnGen and UK Biobank studies. Higher levels of genetically predicted glycemic traits were consistently associated with increased risk of coronary atherosclerosis-related diseases and symptoms. Genetically predicted glycemic traits except HbA1c showed positive associations with peripheral artery disease risk. Genetically predicted FI levels were positively associated with risk of ischemic stroke and chronic kidney disease. Genetically predicted FG and 2hGlu were positively associated with risk of large artery stroke. Genetically predicted 2hGlu levels showed positive associations with risk of small vessel stroke. Higher levels of genetically predicted glycemic traits were not associated with increased risk of thrombotic outcomes. Most associations for genetically predicted levels of 2hGlu and FI remained after adjustment for other glycemic traits. Increase in glycemic status appears to increase risks of coronary and peripheral artery atherosclerosis but not thrombosis.

Sodium-glucose cotransporter 1 inhibition and gout: Mendelian randomisation study

OBJECTIVE

Sodium-glucose cotransporter 2 inhibitors (SGLT2i) reduce serum urate, but their efficacy depends on renal function which is often impaired in people with gout. SGLT1 is primarily expressed in the small intestine and its inhibition may be a more suitable therapeutic target. We aimed to investigate the association of genetically proxied SGLT1i with gout risk, serum urate levels and cardiovascular safety using Mendelian randomisation (MR).

METHODS

Leveraging data from a genome-wide association study of 344,182 individuals in the UK Biobank, we identified a missense variant in the SLC5A1 gene that associated with glycated haemoglobin (HbA1c) to proxy SGLT1i. Outcome genetic data comprised 13,179 gout cases and 750,634 controls, 457,690 individuals for serum urate levels, and up to 977,323 individuals for cardiovascular safety outcomes. We applied the Wald ratio method and investigated potential genetic confounding using colocalization.

RESULTS

The rs17683430 missense variant was selected to instrument SGLT1i. Genetically proxied SGLT1i was associated with 75% reduction in gout risk (OR 0.25; 95%CI 0.06, 0.99; p = 0.048) and 32.0 μmol/L reduction in serum urate (95%CI -56.7, -7.3; p = 0.01), per 6.7 mmol/mol reduction in HbA1c. SGLT1i was associated with increased levels of low-density lipoprotein cholesterol (0.37 mmol/L; 95%CI 0.17, 0.56; p = 0.0002) but not risk of coronary heart disease, stroke, or chronic kidney disease. Colocalization did not suggest that results are attributable to genetic confounding.

CONCLUSION

SGLT1 inhibition may represent a novel therapeutic option for preventing gout in people with or without comorbid diabetes. Randomised trials are needed to formally investigate efficacy and safety.

Multi-context genetic modeling of transcriptional regulation resolves novel disease loci

A majority of the variants identified in genome-wide association studies fall in non-coding regions of the genome, indicating their mechanism of impact is mediated via gene expression. Leveraging this hypothesis, transcriptome-wide association studies (TWAS) have assisted in both the interpretation and discovery of additional genes associated with complex traits. However, existing methods for conducting TWAS do not take full advantage of the intra-individual correlation inherently present in multi-context expression studies and do not properly adjust for multiple testing across contexts. We introduce CONTENT-a computationally efficient method with proper cross-context false discovery correction that leverages correlation structure across contexts to improve power and generate context-specific and context-shared components of expression. We apply CONTENT to bulk multi-tissue and single-cell RNA-seq data sets and show that CONTENT leads to a 42% (bulk) and 110% (single cell) increase in the number of genetically predicted genes relative to previous approaches. We find the context-specific component of expression comprises 30% of heritability in tissue-level bulk data and 75% in single-cell data, consistent with cell-type heterogeneity in bulk tissue. In the context of TWAS, CONTENT increases the number of locus-phenotype associations discovered by over 51% relative to previous methods across 22 complex traits.

Whole-Exome Sequencing (WES) results of 50 patients with chronic kidney diseases: a perspective of Alport syndrome

OBJECTIVE:

Chronic kidney disease (CKD) remains one of the major common health problems, and the number of people affected by the disease is progressively increasing in Turkey and worldwide. This study aimed to investigate molecular defects in Alport syndrome (AS) and other genes in patients with clinically suspected CKD using whole-exome sequencing (WES).

METHODS:

Patients with clinical suspicion of CKD were included in the study. Molecular genetic analyses were performed on genomic DNA by using WES.

RESULTS:

A total of 15 with 5 different pathogenic or likely pathogenic variants were identified in CKD patients, with a diagnostic rate of 30%. Eight variants of uncertain significance were also detected. In this study, 10 variants were described for the first time. As a result, we detected variants associated with CKD in our study population and found AS as the most common CKD after other related kidney diseases.

CONCLUSIONS:

Our results suggest that in heterogeneous diseases such as CKD, WES analysis enables accurate identification of underlying molecular defects promptly. Although CKD accounts for 10–14% of all renal dysfunction, molecular genetic diagnosis is necessary for optimal long-term treatment, prognosis, and effective genetic counseling.

Multimodal single cell sequencing implicates chromatin accessibility and genetic background in diabetic kidney disease progression

The proximal tubule is a key regulator of kidney function and glucose metabolism. Diabetic kidney disease leads to proximal tubule injury and changes in chromatin accessibility that modify the activity of transcription factors involved in glucose metabolism and inflammation. Here we use single nucleus RNA and ATAC sequencing to show that diabetic kidney disease leads to reduced accessibility of glucocorticoid receptor binding sites and an injury-associated expression signature in the proximal tubule. We hypothesize that chromatin accessibility is regulated by genetic background and closely-intertwined with metabolic memory, which pre-programs the proximal tubule to respond differently to external stimuli. Glucocorticoid excess has long been known to increase risk for type 2 diabetes, which raises the possibility that glucocorticoid receptor inhibition may mitigate the adverse metabolic effects of diabetic kidney disease.

Serum bilirubin and kidney function: a Mendelian randomization study

Background

Further investigation is needed to determine the causal effects of serum bilirubin on the risk of chronic kidney disease (CKD).

Methods

This study is a Mendelian randomization (MR) analysis. Among the well-known single-nucleotide polymorphisms (SNPs) related to serum bilirubin levels, rs4149056 in the SLCO1B1 gene was selected as the genetic instrument for single-variant MR analysis, as it was found to be less related to possible confounders than other SNPs. The association between genetic predisposition for bilirubin levels and estimated glomerular filtration rate (eGFR) or CKD was assessed in 337 129 individuals of white British ancestry from the UK Biobank cohort. Two-sample MR based on summary-level data was also performed. SNPs related to total or direct bilirubin levels were collected from a previous genome-wide association study and confounder-associated SNPs were discarded. The independent CKDGen meta-analysis data for CKD were employed as the outcome summary statistics.

Results

The alleles of rs4149056 associated with higher bilirubin levels were associated with better kidney function in the UK Biobank data. In the summary-level MR, both of the genetically predicted total bilirubin {per 5 µmol/L increase; odds ratio [OR] 0.931 [95% confidence interval (CI) 0.871-0.995]} and direct bilirubin [per 1 µmol/L increase; OR 0.910 (95% CI 0.834-0.993)] levels were significantly associated with a lower risk of CKD, supported by the causal estimates from various MR sensitivity analyses.

Conclusion

Genetic predisposition for higher serum bilirubin levels is associated with better kidney function. This result suggests that higher serum bilirubin levels may have causal protective effects against kidney function impairment.

Visualizing novel connections and genetic similarities across diseases using a network-medicine based approach

Understanding the genetic relationships between human disorders could lead to better treatment and prevention strategies, especially for individuals with multiple comorbidities. A common resource for studying genetic-disease relationships is the GWAS Catalog, a large and well curated repository of SNP-trait associations from various studies and populations. Some of these populations are contained within mega-biobanks such as the Million Veteran Program (MVP), which has enabled the genetic classification of several diseases in a large well-characterized and heterogeneous population. Here we aim to provide a network of the genetic relationships among diseases and to demonstrate the utility of quantifying the extent to which a given resource such as MVP has contributed to the discovery of such relations. We use a network-based approach to evaluate shared variants among thousands of traits in the GWAS Catalog repository. Our results indicate many more novel disease relationships that did not exist in early studies and demonstrate that the network can reveal clusters of diseases mechanistically related. Finally, we show novel disease connections that emerge when MVP data is included, highlighting methodology that can be used to indicate the contributions of a given biobank.

Kidney disease progression and collider bias in GWAS

New genome-wide meta-analysis for longitudinal kidney function decline identified several genetic loci related to kidney disease progression. The study illustrated the complexity of modeling longitudinal traits in genome-wide association studies and highlighted the issue of a collider bias that can be introduced when a kidney disease progression phenotype is adjusted for baseline kidney function. Herein, we briefly outline the key findings of this study, their limitations, and implications for future studies.

Genetic insights into therapeutic targets for aortic aneurysms: A Mendelian randomization study

BACKGROUND

As aortic aneurysms (AAs) enlarge, they can become life-threatening if left undiagnosed or neglected. At present, there is a lack of radical treatments for preventing disease progression. Therefore, we aimed to identify effective drug targets that slow the progression of AAs.

METHODS

A Mendelian randomization (MR) analysis was conducted to identify therapeutic targets which are associated with AAs. Summary statistics for AAs were obtained from two datasets: the UK Biobank (2228 cases and 408,565 controls) and the FinnGen study (3658 cases and 244,907 controls). Cis-expression quantitative trait loci (cis-eQTL) for druggable genes were retrieved from the eQTLGen Consortium and used as genetic instrumental variables. Colocalization analysis was performed to determine the probability that single nucleotide polymorphisms (SNPs) associated with AAs and eQTL shared causal genetic variants.

FINDINGS

Four drug targets (BTN3A1, FASN, PLAU, and PSMA4) showed significant MR results in two independent datasets. Proteasome 20S subunit alpha 4 (PSMA4) and plasminogen activator, urokinase (PLAU) in particular, were found to have strong evidence for colocalization with AAs, and abdominal aortic aneurysm in particular. Additionally, except for the association between PSMA4 and intracranial aneurysms, no association between genetically proxied inhibition of PLAU and PSMA4 was detected in increasing the risk of other cardiometabolic risks and diseases.

INTERPRETATION

This study supports that drug-targeting PLAU and PSMA4 inhibition may reduce the risk of AAs.

FUNDING

This work was supported by National Key R&D Program of China (NO. 2017YFC0909400), Nature Science Foundation of China (No. 91839302, 81790624), Project supported by Shanghai Municipal Science and Technology Major Project (Grant No. 2017SHZDZX01), and Tongji Hospital Clinical Research Flagship Program (no. 2019CR207).

Genetic loci and prioritization of genes for kidney function decline derived from a meta-analysis of 62 longitudinal genome-wide association studies

Estimated glomerular filtration rate (eGFR) reflects kidney function. Progressive eGFR-decline can lead to kidney failure, necessitating dialysis or transplantation. Hundreds of loci from genome-wide association studies (GWAS) for eGFR help explain population cross section variability. Since the contribution of these or other loci to eGFR-decline remains largely unknown, we derived GWAS for annual eGFR-decline and meta-analyzed 62 longitudinal studies with eGFR assessed twice over time in all 343,339 individuals and in high-risk groups. We also explored different covariate adjustment. Twelve genome-wide significant independent variants for eGFR-decline unadjusted or adjusted for eGFR-baseline (11 novel, one known for this phenotype), including nine variants robustly associated across models were identified. All loci for eGFR-decline were known for cross-sectional eGFR and thus distinguished a subgroup of eGFR loci. Seven of the nine variants showed variant-by-age interaction on eGFR cross section (further about 350,000 individuals), which linked genetic associations for eGFR-decline with age-dependency of genetic cross-section associations. Clinically important were two to four-fold greater genetic effects on eGFR-decline in high-risk subgroups. Five variants associated also with chronic kidney disease progression mapped to genes with functional in-silico evidence (UMOD, SPATA7, GALNTL5, TPPP). An unfavorable versus favorable nine-variant genetic profile showed increased risk odds ratios of 1.35 for kidney failure (95% confidence intervals 1.03-1.77) and 1.27 for acute kidney injury (95% confidence intervals 1.08-1.50) in over 2000 cases each, with matched controls). Thus, we provide a large data resource, genetic loci, and prioritized genes for kidney function decline, which help inform drug development pipelines revealing important insights into the age-dependency of kidney function genetics.

Genetics of Chronic Kidney Disease in Low-Resource Settings

Advances in kidney genomics in the past 20 years has opened the door for more precise diagnosis of kidney disease and identification of new and specific therapeutic agents. Despite these advances, an imbalance exists between low-resource and affluent regions of the world. Individuals of European ancestry from the United States, United Kingdom, and Iceland account for 16% of the world's population, but represent more than 80% of all genome-wide association studies. South Asia, Southeast Asia, Latin America, and Africa together account for 57% of the world population but less than 5% of genome-wide association studies. Implications of this difference include limitations in new variant discovery, inaccurate interpretation of the effect of genetic variants in non-European populations, and unequal access to genomic testing and novel therapies in resource-poor regions. It also further introduces ethical, legal, and social pitfalls, and ultimately may propagate global health inequities. Ongoing efforts to reduce the imbalance in low-resource regions include funding and capacity building, population-based genome sequencing, population-based genome registries, and genetic research networks. More funding, training, and capacity building for infrastructure and expertise is needed in resource-poor regions. Focusing on this will ensure multiple-fold returns on investments in genomic research and technology.

Association of Chronic Kidney Disease With Risk of Intracerebral Hemorrhage

Importance

The evidence linking chronic kidney disease (CKD) to spontaneous intracerebral hemorrhage (ICH) is inconclusive owing to possible confounding by comorbidities that frequently coexist in patients with these 2 diseases.

Objective

To determine whether there is an association between CKD and ICH risk.

Design, Setting, and Participants

A 3-stage study that combined observational and genetic analyses was conducted. First, the association between CKD and ICH risk was tested in the Ethnic/Racial Variations of Intracerebral Hemorrhage (ERICH) study, a multicenter case-control study in the US. All participants with available data on CKD from ERICH were included. Second, this analysis was replicated in the UK Biobank (UKB), an ongoing population study in the UK. All participants in the UKB were included in this study. Third, mendelian randomization analyses were implemented in the UKB using 27 CKD-related genetic variants to test for genetic associations. ERICH was conducted from August 1, 2010, to August 1, 2017, and observed participants for 1 year. The UKB enrolled participants between 2006 and 2010 and will continue to observe them for 30 years. Data analysis was performed from November 11, 2019, to May 10, 2022.

Exposures

CKD stages 1 to 5.

Main Outcomes and Measures

The outcome of interest was ICH, ascertained in ERICH via expert review of neuroimages and in the UKB via a combination of self-reported data and International Statistical Classification of Diseases, Tenth Revision, codes.

Results

In the ERICH study, a total of 2914 participants with ICH and 2954 controls who had available data on CKD were evaluated (mean [SD] age, 61.6 [14.0] years; 2433 female participants [41.5%]; 3435 male participants [58.5%]); CKD was found to be independently associated with higher risk of ICH (odds ratio [OR], 1.95; 95% CI, 1.35-2.89; P < .001). This association was not modified by race and ethnicity. Replication in the UKB with 1341 participants with ICH and 501 195 controls (mean [SD] age, 56.5 [8.1] years; 273 402 female participants [54.4%]; 229 134 male participants [45.6%]) confirmed this association (OR, 1.28; 95% CI, 1.01-1.62; P = .04). Mendelian randomization analyses indicated that genetically determined CKD was associated with ICH risk (OR, 1.56; 95% CI, 1.13-2.16; P = .007).

Conclusions and Relevance

In this 3-stage study that combined observational and genetic analyses among study participants enrolled in 2 large observational studies with different characteristics and study designs, CKD was consistently associated with higher risk of ICH. Mendelian randomization analyses suggest that this association was causal. Further studies are needed to identify the specific biological pathways that mediate this association.

Genome-wide meta-analysis and omics integration identifies novel genes associated with diabetic kidney disease

AIMS/HYPOTHESIS

Diabetic kidney disease (DKD) is the leading cause of kidney failure and has a substantial genetic component. Our aim was to identify novel genetic factors and genes contributing to DKD by performing meta-analysis of previous genome-wide association studies (GWAS) on DKD and by integrating the results with renal transcriptomics datasets.

METHODS

We performed GWAS meta-analyses using ten phenotypic definitions of DKD, including nearly 27,000 individuals with diabetes. Meta-analysis results were integrated with estimated quantitative trait locus data from human glomerular (N=119) and tubular (N=121) samples to perform transcriptome-wide association study. We also performed gene aggregate tests to jointly test all available common genetic markers within a gene, and combined the results with various kidney omics datasets.

RESULTS

The meta-analysis identified a novel intronic variant (rs72831309) in the TENM2 gene associated with a lower risk of the combined chronic kidney disease (eGFR<60 ml/min per 1.73 m2) and DKD (microalbuminuria or worse) phenotype (p=9.8×10-9; although not withstanding correction for multiple testing, p>9.3×10-9). Gene-level analysis identified ten genes associated with DKD (COL20A1, DCLK1, EIF4E, PTPRN-RESP18, GPR158, INIP-SNX30, LSM14A and MFF; p<2.7×10-6). Integration of GWAS with human glomerular and tubular expression data demonstrated higher tubular AKIRIN2 gene expression in individuals with vs without DKD (p=1.1×10-6). The lead SNPs within six loci significantly altered DNA methylation of a nearby CpG site in kidneys (p<1.5×10-11). Expression of lead genes in kidney tubules or glomeruli correlated with relevant pathological phenotypes (e.g. TENM2 expression correlated positively with eGFR [p=1.6×10-8] and negatively with tubulointerstitial fibrosis [p=2.0×10-9], tubular DCLK1 expression correlated positively with fibrosis [p=7.4×10-16], and SNX30 expression correlated positively with eGFR [p=5.8×10-14] and negatively with fibrosis [p<2.0×10-16]).

CONCLUSIONS/INTERPRETATION

Altogether, the results point to novel genes contributing to the pathogenesis of DKD.

DATA AVAILABILITY

The GWAS meta-analysis results can be accessed via the type 1 and type 2 diabetes (T1D and T2D, respectively) and Common Metabolic Diseases (CMD) Knowledge Portals, and downloaded on their respective download pages ( https://t1d.hugeamp.org/downloads.html ; https://t2d.hugeamp.org/downloads.html ; https://hugeamp.org/downloads.html ).

Mechanisms of ferroptosis in chronic kidney disease

Ferroptosis is a newly identified form of regulated cell death characterized by iron accumulation and lipid peroxidation. Ferroptosis plays an essential role in the pathology of numerous diseases and has emerged as a key area of focus in studies of chronic kidney disease (CKD). CKD is a major public health problem with high incidence and mortality that is characterized by a gradual loss of kidney function over time. The severity and complexity of CKD combined with the limited knowledge of its underlying molecular mechanism(s) have led to increased interest in this disease area. Here, we summarize recent advances in our understanding of the regulatory mechanism(s) of ferroptosis and highlight recent studies describing its role in the pathogenesis and progression of CKD. We further discuss the potential therapeutic benefits of targeting ferroptosis for the treatment of CKD and the major hurdles to overcome for the translation of in vitro studies into the clinic.

An intermediate-effect size variant in UMOD confers risk for chronic kidney disease

The kidney-specific gene UMOD encodes for uromodulin, the most abundant protein excreted in normal urine. Rare large-effect variants in UMOD cause autosomal dominant tubulointerstitial kidney disease (ADTKD), while common low-impact variants strongly associate with kidney function and the risk of chronic kidney disease (CKD) in the general population. It is unknown whether intermediate-effect variants in UMOD contribute to CKD. Here, candidate intermediate-effect UMOD variants were identified using large-population and ADTKD cohorts. Biological and phenotypical effects were investigated using cell models, in silico simulations, patient samples, and international databases and biobanks. Eight UMOD missense variants reported in ADTKD are present in the Genome Aggregation Database (gnomAD), with minor allele frequency (MAF) ranging from 10-5 to 10-3. Among them, the missense variant p.Thr62Pro is detected in ∼1/1,000 individuals of European ancestry, shows incomplete penetrance but a high genetic load in familial clusters of CKD, and is associated with kidney failure in the 100,000 Genomes Project (odds ratio [OR] = 3.99 [1.84 to 8.98]) and the UK Biobank (OR = 4.12 [1.32 to 12.85). Compared with canonical ADTKD mutations, the p.Thr62Pro carriers displayed reduced disease severity, with slower progression of CKD and an intermediate reduction of urinary uromodulin levels, in line with an intermediate trafficking defect in vitro and modest induction of endoplasmic reticulum (ER) stress. Identification of an intermediate-effect UMOD variant completes the spectrum of UMOD-associated kidney diseases and provides insights into the mechanisms of ADTKD and the genetic architecture of CKD.

Multi-Trait Genetic Analysis Reveals Clinically Interpretable Hypertension Subtypes

Background:

Hypertension comprises a heterogeneous range of phenotypes. We asked whether underlying genetic structure could explain a part of this heterogeneity.

Methods:

Our study sample comprised N=198 148 FinnGen participants (56% women, mean age 58 years) and N=21 168 well-phenotyped FINRISK participants (53% women, mean age 50 years). First, we identified genetic hypertension components with an unsupervised Bayesian non-negative matrix factorization algorithm using public genome-wide association data for 144 genetic hypertension variants and 16 clinical traits. For these components, we computed their (1) cross-sectional associations with clinical traits in FINRISK using linear regression and (2) longitudinal associations with incident adverse outcomes in FinnGen using Cox regression.

Results:

We observed 4 genetic hypertension components corresponding to recognizable clinical phenotypes: obesity (high body mass index), dyslipidemia (low high-density lipoprotein cholesterol and high triglycerides), hypolipidemia (low low-density lipoprotein cholesterol and low total cholesterol), and short stature. In FINRISK, all hypertension components had robust associations with their respective clinical characteristics. In FinnGen, the Obesity component was associated with increased diabetes risk (hazard ratio per 1 SD increase 1.08 [Bonferroni corrected CI, 1.05–1.10]) and the Hypolipidemia component with increased autoimmune disease risk (hazard ratio per 1 SD increase 1.05 [Bonferroni corrected CI, 1.03–1.07]). In addition, all hypertension components were related to both hypertension and cardiovascular disease.

Conclusions:

Our unsupervised analysis demonstrates that the genetic basis of hypertension can be understood as a mixture of 4 broad, clinically interpretable components capturing disease heterogeneity. These components could be used to stratify individuals into specific genetic subtypes and, therefore, to benefit personalized health care and pharmaceutical research.

UMOD and the architecture of kidney disease

The identification of genetic factors associated with the risk, onset, and progression of kidney disease has the potential to provide mechanistic insights and therapeutic perspectives. In less than two decades, technological advances yielded a trove of information on the genetic architecture of chronic kidney disease. The spectrum of genetic influence ranges from (ultra)rare variants with large effect size, involved in Mendelian diseases, to common variants, often non-coding and with small effect size, which contribute to polygenic diseases. Here, we review the paradigm of UMOD, the gene coding for uromodulin, to illustrate how a kidney-specific protein of major physiological importance is involved in a spectrum of kidney disorders. This new field of investigation illustrates the importance of genetic variation in the pathogenesis and prognosis of disease, with therapeutic implications.

The Heritability of Kidney Function Using an Older Australian Twin Population

Introduction

Twin studies are unique population models which estimate observed rather than inferred genetic components of complex traits. Nonmonogenic chronic kidney disease (CKD) is a complex disease process with strong genetic and environmental influences, amenable to twin studies. We aimed to assess the heritability of CKD using twin analysis and modeling within Older Australian Twin Study (OATS) data.

Methods

OATS had 109 dizygotic (DZ) and 126 monozygotic (MZ) twin pairs with paired serum creatinine levels. Heritability of kidney function as estimated glomerular filtration rate (eGFR CKD Epidemiology Collaboration [CKD-EPI]) was modeled using the ACE model to estimate additive heritability (A), common (C), and unique (E) environmental factors. Intratwin pair analysis using mixed effects logistic regression allowed analysis of variation in eGFR from established CKD risk factors.

Results

The median age was 69.71 (interquartile range 78.4-83.0) years, with 65% female, and a mean CKD-EPI of 82.8 ml/min (SD 6.7). The unadjusted ACE model determined kidney function to be 33% genetically determined (A), 18% shared genetic-environmental (C), and 49% because of unique environment (E). This remained unchanged when adjusted for age, hypertension, and sex. Hypertension was associated with eGFR; however, intertwin variance in hypertension did not explain variance in eGFR. Two or more hypertension medications were associated with decreased eGFR (P = 0.009).

Conclusion

This study estimates observed heritability at 33%, notably higher than inferred heritability in genome-wide association study (GWAS) (7.1%-18%). Epigenetics and other genomic phenomena may explain this heritability gap. Difference in antihypertension medications explains part of unique environmental exposures, though discordance in hypertension and diabetes does not.

Phenome-Wide Association Study of UMOD Gene Variants and Differential Associations With Clinical Outcomes Across Populations in the Million Veteran Program a Multiethnic Biobank

Introduction

Common variants in the UMOD gene are considered an evolutionary adaptation against urinary tract infections (UTIs) and have been implicated in kidney stone formation, chronic kidney disease (CKD), and hypertension. However, differences in UMOD variant-phenotype associations across population groups are unclear.

Methods

We tested associations between UMOD/PDILT variants and up to 1528 clinical diagnosis codes mapped to phenotype groups in the Million Veteran Program (MVP), using published phenome-wide association study (PheWAS) methodology. Associations were tested using logistic regression adjusted for age, sex, and 10 principal components of ancestry. Bonferroni correction for multiple comparisons was applied.

Results

Among 648,593 veterans, mean (SD) age was 62 (14) years; 9% were female, 19% Black, and 8% Hispanic. In White patients, the rs4293393 UMOD risk variant associated with increased uromodulin was associated with increased odds of CKD (odds ratio [OR]: 1.22, 95% CI: 1.20-1.24, P = 5.90 × 10-111), end-stage kidney disease (OR: 1.17, 95% CI: 1.11-1.24, P = 2.40 × 10-09), and hypertension (OR: 1.03, 95% CI: 1.05-1.05, P = 2.11 × 10-06) and significantly lower odds of UTIs (OR: 0.94, 95% CI: 0.92-0.96, P = 1.21 × 10-10) and kidney calculus (OR: 0.85, 95% CI: 0.83-0.86, P = 4.27 × 10-69). Similar findings were observed across UMOD/PDILT variants. The rs77924615 PDILT variant had stronger associations with acute cystitis in White female (OR: 0.73, 95% CI: 0.59-0.91, P = 4.98 × 10-03) versus male (OR: 0.99, 95% CI: 0.89-1.11, P = 8.80 × 10-01) (P interaction = 0.01) patients. In Black patients, the rs77924615 PDILT variant was significantly associated with pyelonephritis (OR: 0.65, 95% CI: 0.54-0.79, P = 1.05 × 10-05), whereas associations with UMOD promoter variants were attenuated.

Conclusion

Robust associations were observed between UMOD/PDILT variants linked with increased uromodulin expression and lower odds of UTIs and calculus and increased odds of CKD and hypertension. However, these associations varied significantly across ancestry groups and sex.

Association of Genetically Predicted Fibroblast Growth Factor-23 with Heart Failure: A Mendelian Randomization Study

BACKGROUND AND OBJECTIVES

Elevated fibroblast growth factor-23 (FGF23) has been consistently associated with heart failure, particularly heart failure with preserved ejection fraction, among patients with CKD and in the general population. FGF23 may directly induce cardiac remodeling and heart failure. However, biases affecting observational studies impede robust causal inferences. Mendelian randomization leverages genetic determinants of a risk factor to examine causality. We performed a two-sample Mendelian randomization to assess causal associations between FGF23 and heart failure.

DESIGN, SETTING, PARTICIPANTS, & MEASUREMENTS

Genetic instruments were genome-wide significant genetic variants associated with FGF23, including variants near PIP5K1B, RGS14, LINC01229, and CYP24A1. We analyzed data from the Heart Failure Molecular Epidemiology for Therapeutic Targets and BioVU biobanks to examine associations of the four variants with overall heart failure, heart failure with preserved ejection fraction, and heart failure with reduced and mid-range ejection fraction. We developed an eGFR polygenic risk score using summary statistics from the Chronic Kidney Disease Genetics Consortium (CKDGen) genome-wide association study of eGFR in >1 million individuals and performed stratified analyses across eGFR polygenic risk score strata.

RESULTS

Genetically determined FGF23 was not associated with overall heart failure in the Heart Failure Molecular Epidemiology for Therapeutic Targets consortium (odds ratio, 1.13; 95% confidence interval, 0.89 to 1.42 per unit higher genetically predicted log FGF23) and the full BioVU sample (odds ratio, 1.32; 95% confidence interval, 0.95 to 1.84). In stratified analyses in BioVU, higher FGF23 was associated with overall heart failure (odds ratio, 3.09; 95% confidence interval, 1.38 to 6.91) among individuals with low eGFR-polygenic risk score (<1 SD below the mean), but not those with high eGFR-polygenic risk score (P interaction = 0.02). Higher FGF23 was also associated with heart failure with preserved ejection fraction among all BioVU participants (odds ratio, 1.47; 95% confidence interval, 1.01 to 2.14) and individuals with low eGFR-polygenic risk score (odds ratio, 7.20; 95% confidence interval, 2.80 to 18.49), but not those high eGFR-polygenic risk score (P interaction = 2.25 × 10^-4). No significant associations were observed with heart failure with reduced and midrange ejection fraction.

CONCLUSION

We found no association between genetically predicted FGF23 and heart failure in the Heart Failure Molecular Epidemiology for Therapeutic Targets consortium. In BioVU, genetically elevated FGF23 was associated with higher heart failure risk, specifically heart failure with preserved ejection fraction, particularly among individuals with low genetically predicted eGFR.

PODCAST

This article contains a podcast at https://www.asn-online.org/media/podcast/CJASN/2022_07_28_CJN00960122.mp3.

Evidence for correlations between BMI-associated SNPs and circRNAs

Circular RNAs (circRNAs) are regulators of processes like adipogenesis. Their expression can be modulated by SNPs. We analysed links between BMI-associated SNPs and circRNAs. First, we detected an enrichment of BMI-associated SNPs on circRNA genomic loci in comparison to non-significant variants. Analysis of sex-stratified GWAS data revealed that circRNA genomic loci encompassed more genome-wide significant BMI-SNPs in females than in males. To explore whether the enrichment is restricted to BMI, we investigated nine additional GWAS studies. We showed an enrichment of trait-associated SNPs in circRNAs for four analysed phenotypes (body height, chronic kidney disease, anorexia nervosa and autism spectrum disorder). To analyse the influence of BMI-affecting SNPs on circRNA levels in vitro, we examined rs4752856 located on hsa_circ_0022025. The analysis of heterozygous individuals revealed an increased level of circRNA derived from the BMI-increasing SNP allele. We conclude that genetic variation may affect the BMI partly through circRNAs.

Mendelian randomization reveals causal effects of kidney function on various biochemical parameters

The kidney is a vital organ with diverse biological effects and the burden of kidney function impairment is increasing in modern medicine. As the effects from kidney function on diverse biochemical parameters are yet fully understood, additional investigation to reveal the causal effects is warranted. Here we show the causal estimates from kidney function parameter, estimated glomerular filtration rate (eGFR), on 60 biochemical parameters by performing two-sample Mendelian randomization (MR) study in 337,138 white British UK Biobank participants. A higher genetically predicted eGFR was significantly associated with higher lymphocyte percentage, HDL cholesterol, and alanine aminotransferase. The causal estimates indicated that a higher genetically predicted eGFR was associated with lower urea, urate, insulin growth factor-1, and triglycerides levels. The parameters with significant but non-linear causal estimates were hemoglobin concentration, calcium, vitamin D, and urine creatinine values, identified by non-linear MR. Healthcare providers should understand that changes in eGFR may affect the identified biochemical parameters in diverse patterns. Future study is warranted to expand the knowledge of the mechanisms and clinical implications of the causal effects of eGFR on various biochemical parameters.

Benchmarking CKD: incidence of CKD in a European country with low prevalence of CKD and kidney replacement therapy

Iceland was one of six European countries with an adjusted incidence of kidney replacement therapy (KRT) in 2018 lower than 100 per million persons (pmp), along with Estonia, Montenegro, Russia, Serbia and Finland. It was also one of 10 countries with an adjusted KRT prevalence <900 pmp. Furthermore, the prevalence of chronic kidney disease (CKD) in Iceland is up to 2.44-fold lower and the death rate from CKD up to 3.44-fold lower than in other countries with a low incidence of KRT, suggesting that the low KRT incidence actually reflects a low need for KRT rather than low uptake or availability of KRT. This identifies Iceland as a benchmark for countries trying to reduce KRT incidence. Iceland also represents one of the best genetically characterized populations in the world, facilitating studies on the influence of the genetic background versus environment and lifestyle on CKD. This issue of CKJ reports the incidence and risk factors for CKD in Icelandic adults. Diabetes, acute kidney injury, hypertension, cardiovascular disease, chronic lung disease, malignancy and major psychiatric illness were associated with an increased risk of incident CKD, as were obesity and sleep apnea in women. However, in 75% of incident CKD cases, CKD was first detected in category G3 or higher, emphasizing the need for new tools that allow an earlier diagnosis of CKD that precedes the loss of >50% of the functioning kidney mass and/or wider use of albuminuria as a screening tool. The European Society of Cardiology just recommended assessing albuminuria for routine cardiovascular risk workups for all.

Epigenomic and transcriptomic analyses define core cell types, genes and targetable mechanisms for kidney disease

More than 800 million people suffer from kidney disease, yet the mechanism of kidney dysfunction is poorly understood. In the present study, we define the genetic association with kidney function in 1.5 million individuals and identify 878 (126 new) loci. We map the genotype effect on the methylome in 443 kidneys, transcriptome in 686 samples and single-cell open chromatin in 57,229 kidney cells. Heritability analysis reveals that methylation variation explains a larger fraction of heritability than gene expression. We present a multi-stage prioritization strategy and prioritize target genes for 87% of kidney function loci. We highlight key roles of proximal tubules and metabolism in kidney function regulation. Furthermore, the causal role of SLC47A1 in kidney disease is defined in mice with genetic loss of Slc47a1 and in human individuals carrying loss-of-function variants. Our findings emphasize the key role of bulk and single-cell epigenomic information in translating genome-wide association studies into identifying causal genes, cellular origins and mechanisms of complex traits.

ACLY and CKD: A Mendelian Randomization Analysis

Introduction

Adenosine triphosphate-citrate lyase (ACLY) inhibition is a therapeutic strategy under investigation for atherosclerotic cardiovascular disease, nonalcoholic steatohepatitis, and metabolic syndrome. Mouse models suggest that ACLY inhibition could reduce inflammation and kidney fibrosis. Genetic analysis of ACLY in chronic kidney disease (CKD) has not been performed.

Methods

We constructed a genetic instrument by selecting variants associated with ACLY expression in the expression quantitative trait loci genetics consortium (eQTLGen) from blood samples from 31,684 participants. In a 2-sample Mendelian randomization analysis, we evaluated the effect of genetically predicted ACLY expression on the risk of CKD, estimated glomerular filtration rate (eGFR), and albumin-to-creatinine ratio (ACR) using the CKD Genetics (CKDGen) consortium, UK Biobank, and the Finnish Genetics (FinnGen) consortium totaling 66,396 CKD cases and 958,517 controls.

Results

ACLY is constitutively expressed in all cell types including in whole blood. The genetic instrument included 13 variants and explained 1.5% of the variation in whole blood ACLY gene expression. A 34% reduction in ACLY expression score was associated with a 0.04 mmol/l reduced low-density lipoprotein (LDL) cholesterol (P = 3.4 × 10^-4) and a 9% reduced risk of CKD (stages 3, 4, 5, dialysis, or eGFR < 60 ml/min per 1.73 m²) (odds ratio [OR] = 0.91, 95% CI: 0.85-0.98, P = 0.008), but no association was observed with either eGFR or ACR.

Conclusion

Mendelian randomization analyses revealed that genetically reduced ACLY expression was associated with reduced risk of CKD but had no effect on either eGFR or ACR. Further evaluation of ACLY in kidney disease is warranted.

Genome-wide polygenic score to predict chronic kidney disease across ancestries

Chronic kidney disease (CKD) is a common complex condition associated with high morbidity and mortality. Polygenic prediction could enhance CKD screening and prevention; however, this approach has not been optimized for ancestrally diverse populations. By combining APOL1 risk genotypes with genome-wide association studies (GWAS) of kidney function, we designed, optimized and validated a genome-wide polygenic score (GPS) for CKD. The new GPS was tested in 15 independent cohorts, including 3 cohorts of European ancestry (n = 97,050), 6 cohorts of African ancestry (n = 14,544), 4 cohorts of Asian ancestry (n = 8,625) and 2 admixed Latinx cohorts (n = 3,625). We demonstrated score transferability with reproducible performance across all tested cohorts. The top 2% of the GPS was associated with nearly threefold increased risk of CKD across ancestries. In African ancestry cohorts, the APOL1 risk genotype and polygenic component of the GPS had additive effects on the risk of CKD.

Harnessing the Full Potential of Multi-Omic Analyses to Advance the Study and Treatment of Chronic Kidney Disease

Chronic kidney disease (CKD) was the 12th leading cause of death globally in 2017 with the prevalence of CKD estimated at ~9%. Early detection and intervention for CKD may improve patient outcomes, but standard testing approaches even in developed countries do not facilitate identification of patients at high risk of developing CKD, nor those progressing to end-stage kidney disease (ESKD). Recent advances in CKD research are moving towards a more personalised approach for CKD. Heritability for CKD ranges from 30% to 75%, yet identified genetic risk factors account for only a small proportion of the inherited contribution to CKD. More in depth analysis of genomic sequencing data in large cohorts is revealing new genetic risk factors for common diagnoses of CKD and providing novel diagnoses for rare forms of CKD. Multi-omic approaches are now being harnessed to improve our understanding of CKD and explain some of the so-called 'missing heritability'. The most common omic analyses employed for CKD are genomics, epigenomics, transcriptomics, metabolomics, proteomics and phenomics. While each of these omics have been reviewed individually, considering integrated multi-omic analysis offers considerable scope to improve our understanding and treatment of CKD. This narrative review summarises current understanding of multi-omic research alongside recent experimental and analytical approaches, discusses current challenges and future perspectives, and offers new insights for CKD.

Are Alterations in DNA Methylation Related to CKD Development?

The modifications in genomic DNA methylation are involved in the regulation of normal and pathological cellular processes. The epigenetic regulation stimulates biological plasticity as an adaptive response to variations in environmental factors. The role of epigenetic changes is vital for the development of some diseases, including atherogenesis, cancers, and chronic kidney disease (CKD). The results of studies presented in this review have suggested that altered DNA methylation can modulate the expression of pro-inflammatory and pro-fibrotic genes, as well those essential for kidney development and function, thus stimulating renal disease progression. Abnormally increased homocysteine, hypoxia, and inflammation have been suggested to alter epigenetic regulation of gene expression in CKD. Studies of renal samples have demonstrated the relationship between variations in DNA methylation and fibrosis and variations in estimated glomerular filtration rate (eGFR) in human CKD. The unravelling of the genetic-epigenetic profile would enhance our understanding of processes underlying the development of CKD. The understanding of multifaceted relationship between DNA methylation, genes expression, and disease development and progression could improve the ability to identify individuals at risk of CKD and enable the choice of appropriate disease management.

Sex-specific association between coffee consumption and incident chronic kidney disease: a population-based analysis of 359,906 participants from the UK Biobank

BACKGROUND

The risk for chronic kidney disease (CKD) is influenced by genetic predisposition, sex, and lifestyle. Previous research indicates that coffee is a potentially protective factor in CKD. The current study aims to investigate whether sex disparity exists in the coffee-CKD association, and whether genetic risk of CKD or genetic polymorphisms of caffeine metabolism affect this association.

METHODS

A total of 359,906 participants from the UK Biobank who were enrolled between 2006 and 2010 were included in this prospective cohort study, which aimed to estimate the hazard ratios for coffee intake and incident CKD using a Cox proportional hazard model. Allele scores of CKD and caffeine metabolism were additionally adjusted for in a subsample with qualified genetic data ( n = 255,343). Analyses stratified by genetic predisposition, comorbidities, and sex hormones were performed. Tests based on Bayesian model averaging were conducted to ascertain the robustness of the results.

RESULTS

Coffee was inversely associated with CKD in a dose-dependent manner. The effects of coffee did not differ across different strata of genetic risk for CKD, but were more evident among slower genetically predicted caffeine metabolizers. Significant sex disparity was observed ( P value for interaction = 0.013), in that coffee drinking was only associated with the risk reduction of CKD in females. Subgroup analysis revealed that testosterone and sex hormone-binding globulin (SHBG), but not estradiol, modified the coffee-CKD association.

CONCLUSIONS

In addition to the overall inverse coffee-CKD association that was observed in the general population, we could also establish that a sex disparity existed, in that females were more likely to experience the benefit of the association. Testosterone and SHBG may partly account for the sex disparity.

Differential and shared genetic effects on kidney function between diabetic and non-diabetic individuals

Reduced glomerular filtration rate (GFR) can progress to kidney failure. Risk factors include genetics and diabetes mellitus (DM), but little is known about their interaction. We conducted genome-wide association meta-analyses for estimated GFR based on serum creatinine (eGFR), separately for individuals with or without DM (nDM = 178,691, nnoDM = 1,296,113). Our genome-wide searches identified (i) seven eGFR loci with significant DM/noDM-difference, (ii) four additional novel loci with suggestive difference and (iii) 28 further novel loci (including CUBN) by allowing for potential difference. GWAS on eGFR among DM individuals identified 2 known and 27 potentially responsible loci for diabetic kidney disease. Gene prioritization highlighted 18 genes that may inform reno-protective drug development. We highlight the existence of DM-only and noDM-only effects, which can inform about the target group, if respective genes are advanced as drug targets. Largely shared effects suggest that most drug interventions to alter eGFR should be effective in DM and noDM.

Genetics in chronic kidney disease: conclusions from a Kidney Disease: Improving Global Outcomes (KDIGO) Controversies Conference

Numerous genes for monogenic kidney diseases with classical patterns of inheritance, as well as genes for complex kidney diseases that manifest in combination with environmental factors, have been discovered. Genetic findings are increasingly used to inform clinical management of nephropathies, and have led to improved diagnostics, disease surveillance, choice of therapy, and family counseling. All of these steps rely on accurate interpretation of genetic data, which can be outpaced by current rates of data collection. In March of 2021, Kidney Diseases: Improving Global Outcomes (KDIGO) held a Controversies Conference on "Genetics in Chronic Kidney Disease (CKD)" to review the current state of understanding of monogenic and complex (polygenic) kidney diseases, processes for applying genetic findings in clinical medicine, and use of genomics for defining and stratifying CKD. Given the important contribution of genetic variants to CKD, practitioners with CKD patients are advised to "think genetic," which specifically involves obtaining a family history, collecting detailed information on age of CKD onset, performing clinical examination for extrarenal symptoms, and considering genetic testing. To improve the use of genetics in nephrology, meeting participants advised developing an advanced training or subspecialty track for nephrologists, crafting guidelines for testing and treatment, and educating patients, students, and practitioners. Key areas of future research, including clinical interpretation of genome variation, electronic phenotyping, global representation, kidney-specific molecular data, polygenic scores, translational epidemiology, and open data resources, were also identified.

Causal Association between Chronic Kidney Disease and Risk of 19 Site-Specific Cancers: A Mendelian Randomization Study

BACKGROUND

Results of previous observational studies examining the risk of cancer among patients with chronic kidney disease (CKD) are conflicting. We here explore the causal relationship between estimated glomerular filtration rate (eGFR) and albuminuria, two principal measurements of CKD, and 19 site-specific cancers using Mendelian randomization (MR) analysis.

METHODS

Single-nucleotide polymorphisms reported to be strongly correlated with eGFR and albuminuria in recent large genome-wide association studies were used as instrumental variables to investigate the causal relationship with cancer using summary-level statistics from several cancer-specific consortia, as well as data of 347,408 participants in the UK Biobank and 260,405 participants in the FinnGen.

RESULTS

Our data showed that impaired kidney function was associated with higher odds of leukemia [OR = 1.23; 95% confidence interval (CI), 1.06-1.43; P = 0.007], cervical cancer (OR = 1.22; 95% CI, 1.04-1.43; P = 0.017), and female renal cell carcinoma (OR = 1.4; 95% CI, 1.12-1.77; P = 0.004), per 10% decrease in eGFR. The ORs were 1.21 (95% CI, 1.07-1.36; P = 0.002) for colorectal cancer and 0.76 (95% CI, 0.62-0.92; P = 0.006) for non-Hodgkin lymphoma, per doubling odds of albuminuria. In multivariable MR, effect sizes of eGFR-cervical cancer remained strong after adjusting for confounders.

CONCLUSIONS

The current study indicates that progression of CKD contributes to carcinogenesis of renal cell carcinoma, leukemia, cervical, and colorectal cancer.

IMPACT

The potential association of kidney function and albuminuria with certain cancers warrants further investigation in order to provide appropriate recommendations regarding cancer screening among patients with CKD.

A multiancestry genome-wide association study of unexplained chronic ALT elevation as a proxy for nonalcoholic fatty liver disease with histological and radiological validation

Nonalcoholic fatty liver disease (NAFLD) is a growing cause of chronic liver disease. Using a proxy NAFLD definition of chronic elevation of alanine aminotransferase (cALT) levels without other liver diseases, we performed a multiancestry genome-wide association study (GWAS) in the Million Veteran Program (MVP) including 90,408 cALT cases and 128,187 controls. Seventy-seven loci exceeded genome-wide significance, including 25 without prior NAFLD or alanine aminotransferase associations, with one additional locus identified in European American-only and two in African American-only analyses (P < 5 × 10-8). External replication in histology-defined NAFLD cohorts (7,397 cases and 56,785 controls) or radiologic imaging cohorts (n = 44,289) replicated 17 single-nucleotide polymorphisms (SNPs) (P < 6.5 × 10-4), of which 9 were new (TRIB1, PPARG, MTTP, SERPINA1, FTO, IL1RN, COBLL1, APOH and IFI30). Pleiotropy analysis showed that 61 of 77 multiancestry and all 17 replicated SNPs were jointly associated with metabolic and/or inflammatory traits, revealing a complex model of genetic architecture. Our approach integrating cALT, histology and imaging reveals new insights into genetic liability to NAFLD.

Improving lithium dose prediction using population pharmacokinetics and pharmacogenomics: a cohort genome-wide association study in Sweden

BACKGROUND

Lithium is the most effective treatment for bipolar disorder, resulting in strong suicide prevention effects. The therapeutic range of lithium, however, is narrow and treatment initiation requires individual titration to address inter-individual variability. We aimed to improve lithium dose prediction using clinical and genomic data.

METHODS

We performed a population pharmacokinetic study followed by a genome-wide association study (GWAS), including two clinical Swedish cohorts. Participants in cohort 1 were from specialised outpatient clinics at Huddinge Hospital, in Stockholm, Sweden, and participants in cohort 2 were identified using the Swedish National Quality Registry for Bipolar disorder (BipoläR). Patients who received a lithium dose corresponding to at least one tablet of lithium sulphate (6 mmol) per day and had clinically relevant plasma concentrations of lithium were included in the study. Data on age, sex, bodyweight, height, creatinine concentration, estimated glomerular filtration rate (eGFR), lithium preparation, number of tablets of lithium per day, serum lithium concentration, and medications affecting kidney function (C09 antihypertensives, C03 [except C03D] sodium-retaining diuretics, and non-steroidal anti-inflammatory drugs) were obtained retrospectively for several timepoints when possible from electronic health records, BipoläR, and the Swedish prescription registry. The median time between timepoints was 1·07 years for cohort 1 and 1·09 years for cohort 2. The primary outcome of interest was the natural logarithm of total body clearance for lithium (CL_Li) associated with the clinical variables. The residual effects after accounting for age and sex, representing the individual-level effects (CL_Li,age/sex), were used as the dependent variable in a GWAS.

FINDINGS

2357 patients who were administered lithium (1423 women [60·4%] and 934 men [39·6%]; mean age 53·6 years [range 17-89], mainly of European descent) were included and 5627 data points were obtained. Age (variance explained [R²]: R²_cohort1=0·41 and R²_cohort2=0·31; both p<0·0001), sex (R²_cohort1=0·0063 [p=0·045] and R²_cohort2=0·026 [p<0·0001]), eGFR (R²_cohort1=0·38 and R²_cohort2=0·20; both p<0·0001), comedication with diuretics (R²_cohort1=0·0058 [p=0·014] and R²_cohort2=0·0026 [p<0·0001]), and agents acting on the renin-aldosterone-angiotensin system (R²_cohort1=0·028 and R²_cohort2=0·015; both p<0·0001) were clinical predictors of CL_Li. Notably, an association between CL_Li and serum lithium was observed, with a lower CL_Li being associated with higher serum lithium (R²_cohort1=0·13 and R²_cohort2=0·15; both p<0·0001). In a GWAS of CL_Li,age/sex, one locus was associated with a change in CL_Li (rs583503; β=-0·053 [95% CI -0·071 to -0·034]; p<0·00000005). We also found enrichment of the associations with genes expressed in the medulla (p=0·0014, corrected FDR=0·04) and cortex of the kidney (p=0·0015, corrected FDR=0·04), as well as associations with polygenic risk scores for eGFR (p value threshold: 0·05, p=0·01), body-mass index (p value threshold: 0·05, p=0·00025), and blood urea nitrogen (p value threshold: 0·001, p=0·00043). The model based on six clinical predictors explained 61·4% of the variance in CL_Li in cohort 1 and 49·8% in cohort 2. Adding genetic markers did not lead to major improvement of the models: within the subsample of genotyped individuals, the variance explained only increased from 59·32% to 59·36% in cohort 1 and from 49·21% to 50·03% in cohort 2 when including rs583503 and the four first principal components.

INTERPRETATION

Our model predictors could be used clinically to better guide lithium dosage, shortening the time to reach therapeutic concentrations, thus improving care. Identification of the first genomic locus and PRS to be associated with CL_Li introduces the opportunity of individualised medicine in lithium treatment.

FUNDING

Stanley Medical Research Institute, Swedish Research Council, Swedish Foundation for Strategic Research, Swedish Brain Foundation, Swedish Research Council, Söderström-Königska Foundation, Bror Gadelius Minnesfond, Swedish Mental Health Fund, Karolinska Institutet and Hospital.

Genome-wide studies reveal factors associated with circulating uromodulin and its relationships to complex diseases

Uromodulin (UMOD) is a major risk gene for monogenic and complex forms of kidney disease. The encoded kidney-specific protein uromodulin is highly abundant in urine and related to chronic kidney disease, hypertension, and pathogen defense. To gain insights into potential systemic roles, we performed genome-wide screens of circulating uromodulin using complementary antibody-based and aptamer-based assays. We detected 3 and 10 distinct significant loci, respectively. Integration of antibody-based results at the UMOD locus with functional genomics data (RNA-Seq, ATAC-Seq, Hi-C) of primary human kidney tissue highlighted an upstream variant with differential accessibility and transcription in uromodulin-synthesizing kidney cells as underlying the observed cis effect. Shared association patterns with complex traits, including chronic kidney disease and blood pressure, placed the PRKAG2 locus in the same pathway as UMOD. Experimental validation of the third antibody-based locus, B4GALNT2, showed that the p.Cys466Arg variant of the encoded N-acetylgalactosaminyltransferase had a loss-of-function effect leading to higher serum uromodulin levels. Aptamer-based results pointed to enzymes writing glycan marks present on uromodulin and to their receptors in the circulation, suggesting that this assay permits investigating uromodulin's complex glycosylation rather than its quantitative levels. Overall, our study provides insights into circulating uromodulin and its emerging functions.

Association between Vitamin Intake and Chronic Kidney Disease According to a Variant Located Upstream of the PTGS1 Gene: A Cross-Sectional Analysis of Shika Study

Chronic kidney disease (CKD) patients have been advised to take vitamins; however, the effects have been controversial. The individual differences in developing CKD might involve genetic variants of inflammation, including variant rs883484 located upstream of the prostaglandin-endoperoxide synthase 1 (PTGS1) gene. We aimed to identify whether the 12 dietary vitamin intake interacts with genotypes of the rs883484 on developing CKD. The population-based, cross-sectional study had 684 Japanese participants (≥40 years old). The study used a validated, brief, self-administered diet history questionnaire to estimate the intake of the dietary vitamins. CKD was defined as estimated glomerular filtration < 60 mL/min/1.73 m2. The study participants had an average age of 62.1 ± 10.8 years with 15.4% minor homozygotes of rs883484, and 114 subjects had CKD. In the fully adjusted model, the higher intake of vitamins, namely niacin (odds ratio (OR) = 0.74, 95% confidence interval (CI): 0.57−0.96, p = 0.024), α-tocopherol (OR = 0.49, 95% CI: 0.26−0.95, p = 0.034), and vitamin C (OR = 0.97, 95% CI: 0.95−1.00, p = 0.037), was independently associated with lower CKD tendency in the minor homozygotes of rs883484. The results suggested the importance of dietary vitamin intake in the prevention of CKD in middle-aged to older-aged Japanese with minor homozygous of rs883484 gene variant.

Cardiovascular disease protein biomarkers are associated with kidney function: The Framingham Heart Study

Background

Biomarkers common to chronic kidney disease (CKD) and cardiovascular disease (CVD) may reflect early impairments underlying both diseases.

Methods

We evaluated associations of 71 CVD-related plasma proteins measured in 2,873 Framingham Heart Study (FHS) Offspring cohort participants with cross-sectional continuous eGFR and with longitudinal change in eGFR from baseline to follow-up (ΔeGFR). We also evaluated the associations of the 71 CVD proteins with the following dichotomous secondary outcomes: prevalent CKD stage ≥3 (cross-sectional), new-onset CKD stage ≥3 (longitudinal), and rapid decline in eGFR (longitudinal). Proteins significantly associated with eGFR and ΔeGFR were subsequently validated in 3,951 FHS Third Generation cohort participants and were tested using Mendelian randomization (MR) analysis to infer putatively causal relations between plasma protein biomarkers and kidney function.

Results

In cross-sectional analysis, 37 protein biomarkers were significantly associated with eGFR at FDR<0.05 in the FHS Offspring cohort and 20 of these validated in the FHS Third Generation cohort at p<0.05/37. In longitudinal analysis, 27 protein biomarkers were significantly associated with ΔeGFR at FDR<0.05 and 12 of these were validated in the FHS Third Generation cohort at p<0.05/27. Additionally, 35 protein biomarkers were significantly associated with prevalent CKD, five were significantly associated with new-onset CKD, and 17 were significantly associated with rapid decline in eGFR. MR suggested putatively causal relations of melanoma cell adhesion molecule (MCAM; -0.011±0.003 mL/min/1.73m², p = 5.11E-5) and epidermal growth factor-containing fibulin-like extracellular matrix protein 1 (EFEMP1; -0.006±0.002 mL/min/1.73m², p = 0.0001) concentration with eGFR.

Discussion/conclusions

Eight protein biomarkers were consistently associated with eGFR in cross-sectional and longitudinal analysis in both cohorts and may capture early kidney impairment; others were implicated in association and causal inference analyses. A subset of CVD protein biomarkers may contribute causally to the pathogenesis of kidney impairment and should be studied as targets for CKD treatment and early prevention.

Polygenic risk scores for CARDINAL study

The Cardiometabolic Disorders in African-Ancestry Populations (CARDINAL) study site is a well-powered, first-of-its-kind resource for developing, refining and validating methods for research into polygenic risk scores that accounts for local ancestry, to improve risk prediction in diverse populations.

Sex-specific effect of serum urate levels on coronary heart disease and myocardial infarction prevention: A Mendelian randomization study

BACKGROUND AND AIMS

Observational studies have examined serum urate levels in relation to coronary heart disease (CHD) and myocardial infarction (MI). Whether these associations are causal remains controversial, due to confounding factors and reverse causality. We aim to investigate the causality of these associations using Mendelian randomization method.

METHODS AND RESULTS

Instrumental variables were obtained from the largest genome-wide association studies of serum urate (457,690 individuals) to date. Summary statistics were from CARDIoGRAMplusC4D consortium (60,801 CHD cases; 43,676 MI cases), FinnGen (21,012 CHD cases; 12,801 MI cases), UK Biobank (10,157 CHD cases; 7018 MI cases), and Biobank Japan (29,319 CHD cases). Inverse-variance weighted method was applied as the main results. Other statistical methods and reverse MR analysis were conducted in the supplementary analyses. Elevated genetically determined serum urate levels were associated with increased risks of CHD and MI. The association pattern remained for the datasets in FinnGen, the combined results of three independent data sources (CHD: odds ratio (OR), 1.10; 95%CI, 1.06-1.15; p = 4.2 × 10^-6; MI: OR, 1.12; 95%CI, 1.07-1.18; p = 2.7 × 10^-6), and East Asian population. Interestingly, sex-specific subgroup analyses revealed that these associations kept in men only, but not among women in individuals of European ancestry. No consistent evidence was found for the causal effect of CHD or MI on serum urate levels.

CONCLUSION

We provide consistent evidence for the causal effect of genetically predicted serum urate levels on CHD and MI, but not the reverse effect. Urate-lowering therapy may be of cardiovascular benefit in the prevention of CHD and MI, especially for men.

Genetic Variations of Vitamin A-Absorption and Storage-Related Genes, and Their Potential Contribution to Vitamin A Deficiency Risks Among Different Ethnic Groups

Vitamin A, an essential fat-soluble micronutrient, plays a critical role in the body, by regulating vision, immune responses, and normal development, for instance. Vitamin A deficiency (VAD) is a major cause of xerophthalmia and increases the risk of death from infectious diseases. It is also emerging that prenatal exposure to VAD is associated with disease risks later in life. The overall prevalence of VAD has significantly declined over recent decades; however, the rate of VAD is still high in many low- and mid-income countries and even in high-income countries among specific ethnic/race groups. While VAD occurs when dietary intake is insufficient to meet demands, establishing a strong association between food insecurity and VAD, and vitamin A supplementation is the primary solution to treat VAD, genetic contributions have also been reported to effect serum vitamin A levels. In this review, we discuss genetic variations associated with vitamin A status and vitamin A bioactivity-associated genes, specifically those linked to uptake of the vitamin in the small intestine and its storage in the liver, as well as their potential contribution to vitamin A deficiency risks among different ethnic groups.

Therapeutic Targets for Heart Failure Identified Using Proteomics and Mendelian Randomization

BACKGROUND

Heart failure (HF) is a highly prevalent disorder for which disease mechanisms are incompletely understood. The discovery of disease-associated proteins with causal genetic evidence provides an opportunity to identify new therapeutic targets.

METHODS

We investigated the observational and causal associations of 90 cardiovascular proteins, which were measured using affinity-based proteomic assays. First, we estimated the associations of 90 cardiovascular proteins with incident heart failure by means of a fixed-effect meta-analysis of 4 population-based studies, composed of a total of 3019 participants with 732 HF events. The causal effects of HF-associated proteins were then investigated by Mendelian randomization, using cis-protein quantitative loci genetic instruments identified from genomewide association studies in more than 30 000 individuals. To improve the precision of causal estimates, we implemented an Mendelian randomization model that accounted for linkage disequilibrium between instruments and tested the robustness of causal estimates through a multiverse sensitivity analysis that included up to 120 combinations of instrument selection parameters and Mendelian randomization models per protein. The druggability of candidate proteins was surveyed, and mechanism of action and potential on-target side effects were explored with cross-trait Mendelian randomization analysis.

RESULTS

Forty-four of ninety proteins were positively associated with risk of incident HF (P<6.0×10-4). Among these, 8 proteins had evidence of a causal association with HF that was robust to multiverse sensitivity analysis: higher CSF-1 (macrophage colony-stimulating factor 1), Gal-3 (galectin-3) and KIM-1 (kidney injury molecule 1) were positively associated with risk of HF, whereas higher ADM (adrenomedullin), CHI3L1 (chitinase-3-like protein 1), CTSL1 (cathepsin L1), FGF-23 (fibroblast growth factor 23), and MMP-12 (matrix metalloproteinase-12) were protective. Therapeutics targeting ADM and Gal-3 are currently under evaluation in clinical trials, and all the remaining proteins were considered druggable, except KIM-1.

CONCLUSIONS

We identified 44 circulating proteins that were associated with incident HF, of which 8 showed evidence of a causal relationship and 7 were druggable, including adrenomedullin, which represents a particularly promising drug target. Our approach demonstrates a tractable roadmap for the triangulation of population genomic and proteomic data for the prioritization of therapeutic targets for complex human diseases.

Apolipoprotein L1 variability is associated with increased risk of renal failure in the Czech population

BACKGROUND

With stage 5 chronic kidney disease (CKD5) more prevalent in the Czech Republic than in most European countries, genetic susceptibility is potentially implicated.

METHODS

In a group of 1489 CKD5 kidney transplantation patients (93% with complete clinical characteristics; mean age 52.0 years, 37% females) and 2559 healthy controls (mean age 49.0 years, 51% females), we examined the prevalence of six APOL1 SNPs (rs73885319, rs71785313, rs13056427, rs136147, rs10854688 and rs9610473) and one newly detected 55-nucleotide insertion/deletion polymorphism.

RESULTS

The rs73885319 and rs71785313 variants were monomorphic in the Czech Caucasian population. Genotype frequencies of the three SNPs examined (rs13056427, rs136147 and rs9610473) were almost identical in patients and controls (all P values were between 0.39 and 0.91). Minor homozygotes of rs10854688 were more common between the patients (13.2%) than in controls (10.7%) (OR [95% CI]; 1.32 [1.08-1.64]; P < 0.01). Prevalence of the newly detected 55-bp APOL1 deletion was significantly higher in CKD5 patients (3.0% vs. 1.7%; OR [95% CI]; 1.80 [1.16-2.80]; P < 0.01) compared to controls. Frequencies of some individual APOL1 haplotypes were borderline different between patients and controls.

CONCLUSION

We found an association between rs10854688 SNP within the APOL1 gene and end-stage renal disease in the Czech Caucasian population. Further independent studies are required before a conclusive association between the newly detected APOL1 insertion/deletion polymorphism and CKD5 can be confirmed.

Genome-wide meta-analysis revealed several genetic loci associated with serum uric acid levels in Korean population: an analysis of Korea Biobank data

The serum uric acid (SUA) level is an important determinant of gout, hypertension, metabolic syndrome, and cardiovascular disease. Although previous genome-wide studies have identified multiple genetic variants associated with SUA, most genetic analyses have focused on individuals with European ancestry; thus, understanding of the genetic architecture of SUA is currently limited for Asian populations. We conducted a genome-wide meta-analysis based on Korea Biobank data consistent with three cohorts; namely, the Korean Genome and Epidemiology Study (KoGES) Ansan and Ansung, KoGES Health Examinee, and KoGES Cardiovascular Disease Association studies. In total, 60,585 participants aged ≥40 years were included in the analysis of the three cohorts. We used logistic regression analyses to perform genome-wide association study (GWAS) adjustments for confounding variables. Subsequently, a meta-analysis was conducted by combining the analyses of the three GWASs. We identified 8,105 variants at 22 genetic loci with a P value < 5 × 10^-8. Among these, six novel genetic loci associated with SUA in the Korean population were identified (rs4715517 in HCRTR2, rs145099458 in 3.2 kb 3' of MLXIPL, rs1137642 in B4GALT1, rs659107 in LOC105378410, rs7919329 in LOC107984274, and rs2240751 in MFSD12). Our meta-analysis provides insights into the genetic architecture of SUA in the Korean population. Further studies are warranted to replicate the study results and elucidate the specific role of these variants in SUA homeostasis.