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

Genetic background and transplantation outcomes: insights from genome-wide association studies

PURPOSE OF REVIEW

The current review summarizes recent advances in the genetic studies of transplantation outcomes, including new genome-wide association studies for acute rejection, allograft survival, pharmacogenomics, and common transplant comorbidities.

RECENT FINDINGS

Genetic studies of kidney transplantation outcomes have begun to address the question of genetic compatibility beyond human leukocyte antigens, including the role of genome-wide mismatches in missense variants, and the 'genomic collision' hypothesis under which the risk of rejection may be increased in recipients homozygous for loss-of-function variants with grafts from nonhomozygous donors. In recent pilot studies, missense mismatch scores for transmembrane and secreted proteins were associated with antibodies against the mismatched peptides and reduced allograft survival. A 'genomic collision' at the LIMS1 locus involving a common deletion near LIMS1 gene was associated with anti-LIMS1 antibody response and increased risk of rejection. Additional genetic factors under active investigation include genome-wide polygenic risk scores for renal function and apolipoprotein L1 risk genotypes in African-American kidney donors. Due to the heterogeneity and complexity of clinical outcomes, new genome-wide association studies for rejection, allograft survival, and specific transplant comorbidities will require larger multicenter meta-analyses.

SUMMARY

Genetic compatibilities between donor and recipient represent an important determinant of rejection and long-term allograft survival. Genetic background of transplant donors may be additionally predictive of allograft function, while recipient's genomes are likely determinant of a wide range of transplantation outcomes, from rejection susceptibility to pharmacogenetics and various comorbidities related to prolonged immunosuppression.

Causal associations between urinary sodium with body mass, shape and composition: a Mendelian randomization study

Observational studies have found associations between urinary sodium (UNa) with obesity, body shape and composition; but the findings may be biased by residual confounding. The objective of this two-sample Mendelian randomization (MR) study was to analyze their causal associations in both sex-combined and sex-specific models. Genome-wide association studies of UNa, body mass index (BMI), BMI-adjusted waist-to-hip ratio (WHR), body fat (BF) percentage and estimated glomerular filtration rate (eGFR) were identified. We initially extracted fifty SNPs associated with UNa at significance level of 5 × 10^–8, but further removed those SNPs with potential horizontal pleiotropy. Univariable and multivariable MR with adjustment for eGFR were performed. Inverse-variance weighted MR was performed as the primary analysis, with MR-Egger methods as sensitivity analysis. The potential bidirectional association between BMI and UNa was investigated. All exposure and outcomes were continuous, and the effect measure was regression coefficients (beta) and their 95% confidence intervals (95% CI). The total sample size was up to 322 154. UNa was causally associated with increased BMI in both men [eGFR-adjusted beta 0.443 (0.163–0.724)] and women [0.594 (0.333–0.855)]. UNa caused BF percentage increase in men [0.622 (0.268–0.976)] and women [0.334 (0.007–0.662)]. UNa significantly elevated BMI-adjusted WHR in men [0.321 (0.094–0.548)], but not in women [0.170 (− 0.052 to 0.391)]. Additionally, we found that BMI causally increased UNa [0.043 (0.023–0.063)]. UNa increased BMI and BF percentage. Salt intake affects male body shape by increasing BMI-adjusted WHR, but showed no effects on female body shape. The bidirectional association between BMI and UNa suggested that salt reduction measures and weight reduction measures should be implemented simultaneously to break the vicious cycle and gain more health benefits.

Short or Long Sleep Duration and CKD: A Mendelian Randomization Study

BACKGROUND

Studies have found sleeping behaviors, such as sleep duration, to be associated with kidney function and cardiovascular disease risk. However, whether short or long sleep duration is a causative factor for kidney function impairment has been rarely studied.

METHODS

We studied data from participants aged 40-69 years in the UK Biobank prospective cohort, including 25,605 self-reporting short-duration sleep (<6 hours per 24 hours), 404,550 reporting intermediate-duration sleep (6-8 hours), and 35,659 reporting long-duration sleep (≥9 hours) in the clinical analysis. Using logistic regression analysis, we investigated the observational association between the sleep duration group and prevalent CKD stages 3-5, analyzed by logistic regression analysis. We performed Mendelian randomization (MR) analysis involving 321,260 White British individuals using genetic instruments (genetic variants linked with short- or long-duration sleep behavior as instrumental variables). We performed genetic risk score analysis as a one-sample MR and extended the finding with a two-sample MR analysis with CKD outcome information from the independent CKDGen Consortium genome-wide association study meta-analysis.

RESULTS

Short or long sleep duration clinically associated with higher prevalence of CKD compared with intermediate duration. The genetic risk score for short (but not long) sleep was significantly related to CKD (per unit reflecting a two-fold increase in the odds of the phenotype; adjusted odds ratio, 1.80; 95% confidence interval, 1.25 to 2.60). Two-sample MR analysis demonstrated causal effects of short sleep duration on CKD by the inverse variance weighted method, supported by causal estimates from MR-Egger regression.

CONCLUSIONS

These findings support an adverse effect of a short sleep duration on kidney function. Clinicians may encourage patients to avoid short-duration sleeping behavior to reduce CKD risk.

Genetic Pleiotropy of Bone-Related Phenotypes: Insights from Osteoporosis

PURPOSE OF REVIEW

We summarize recent evidence on the shared genetics within and outside the musculoskeletal system (mostly related to bone density and osteoporosis).

RECENT FINDINGS

Osteoporosis is determined by an interplay between multiple genetic and environmental factors. Significant progress has been made regarding its genetic background revealing a number of robustly validated loci and respective pathways. However, pleiotropic factors affecting bone and other tissues are not well understood. The analytical methods proposed to test for potential associations between genetic variants and multiple phenotypes can be applied to bone-related data. A number of recent genetic studies have shown evidence of pleiotropy between bone density and other different phenotypes (traits, conditions, or diseases), within and outside the musculoskeletal system. Power benefits of combining correlated phenotypes, as well as unbiased discovery, make these studies promising. Studies in humans are supported by evidence from animal models. Drug development and repurposing should benefit from the pleiotropic approach. We believe that future studies should take into account shared genetics between the bone and related traits.

A multi-ethnic meta-analysis identifies novel genes, including ACSL5, associated with amyotrophic lateral sclerosis

Amyotrophic lateral sclerosis (ALS) is a devastating progressive motor neuron disease that affects people of all ethnicities. Approximately 90% of ALS cases are sporadic and thought to have multifactorial pathogenesis. To understand the genetics of sporadic ALS, we conducted a genome-wide association study using 1,173 sporadic ALS cases and 8,925 controls in a Japanese population. A combined meta-analysis of our Japanese cohort with individuals of European ancestry revealed a significant association at the ACSL5 locus (top SNP p = 2.97 × 10⁻⁸). We validated the association with ACSL5 in a replication study with a Chinese population and an independent Japanese population (1941 ALS cases, 3821 controls; top SNP p = 1.82 × 10⁻⁴). In the combined meta-analysis, the intronic ACSL5 SNP rs3736947 showed the strongest association (p = 7.81 × 10⁻¹¹). Using a gene-based analysis of the full multi-ethnic dataset, we uncovered additional genes significantly associated with ALS: ERGIC1, RAPGEF5, FNBP1, and ATXN3. These results advance our understanding of the genetic basis of sporadic ALS.

Gen Sobue, Masao Nagasaki and colleagues report a genome-wide association study for amyotrophic lateral sclerosis (ALS) in a large, multi-ethnic cohort comprising Japanese, Chinese, and European ancestry populations. They find a significant association to variants within the ACSL5 gene and identify novel associations with 4 additional genes using a gene-based approach.

Multi-omics study for interpretation of genome-wide association study

Genome-wide association studies (GWASs) have identified thousands of genetic loci associated with complex traits, including a wide variety of diseases. Despite the successful identification of associated loci, interpreting GWAS findings remains challenging and requires other biological resources. Omics, including genomics, transcriptomics, proteomics, metabolomics, and epigenomics, are increasingly used in a broad range of research fields. Integrative analyses applying GWAS with these omics data are expected to expand our knowledge of complex traits and provide insight into the pathogenesis of complex diseases and their causative factors. Recently, associations between genetic variants and omics data have been comprehensively evaluated, providing new information on the influence of genetic variants on omics. Furthermore, recent advances in analytic methods, including single-cell technologies, have revealed previously unknown disease mechanisms. To advance our understanding of complex traits, integrative analysis using GWAS with multi-omics data is needed. In this review, I describe successful examples of integrative analyses based on omics and GWAS, discuss the limitations of current multi-omics analyses, and provide a perspective on future integrative studies.

An investigation of causal relationships between prediabetes and vascular complications

Prediabetes is a state of glycaemic dysregulation below the diagnostic threshold of type 2 diabetes (T2D). Globally, ~352 million people have prediabetes, of which 35-50% develop full-blown diabetes within five years. T2D and its complications are costly to treat, causing considerable morbidity and early mortality. Whether prediabetes is causally related to diabetes complications is unclear. Here we report a causal inference analysis investigating the effects of prediabetes in coronary artery disease, stroke and chronic kidney disease, complemented by a systematic review of relevant observational studies. Although the observational studies suggest that prediabetes is broadly associated with diabetes complications, the causal inference analysis revealed that prediabetes is only causally related with coronary artery disease, with no evidence of causal effects on other diabetes complications. In conclusion, prediabetes likely causes coronary artery disease and its prevention is likely to be most effective if initiated prior to the onset of diabetes.

Associations of maternal early-pregnancy blood glucose and insulin concentrations with DNA methylation in newborns

BACKGROUND

Intrauterine exposure to a disturbed maternal glucose metabolism is associated with adverse offspring outcomes. DNA methylation is a potential mechanism underlying these associations. We examined whether maternal early-pregnancy glucose and insulin concentrations are associated with newborn DNA methylation. In a population-based prospective cohort study among 935 pregnant women, maternal plasma concentrations of non-fasting glucose and insulin were measured at a median of 13.1 weeks of gestation (95% range 9.4-17.4). DNA methylation was measured using the Infinium HumanMethylation450 BeadChip (Ilumina). We analyzed associations of maternal early-pregnancy glucose and insulin concentrations with single-CpG DNA methylation using robust linear regression models. Differentially methylated regions were analyzed using the dmrff package in R. We stratified the analyses on normal weight versus overweight or obese women. We also performed a look-up of CpGs and differently methylated regions from previous studies to be associated with maternal gestational diabetes, hyperglycemia or hyperinsulinemia, or with type 2 diabetes in adults.

RESULTS

Maternal early-pregnancy glucose and insulin concentrations were not associated with DNA methylation at single CpGs nor with differentially methylated regions in the total group. In analyses stratified on maternal BMI, maternal early-pregnancy glucose concentrations were associated with DNA methylation at one CpG (cg03617420, XKR6) among normal weight women and at another (cg12081946, IL17D) among overweight or obese women. No stratum-specific associations were found for maternal early-pregnancy insulin concentrations. The two CpGs were not associated with birth weight or childhood glycemic measures (p values > 0.1). Maternal early-pregnancy insulin concentrations were associated with one CpG known to be related to adult type 2 diabetes. Enrichment among nominally significant findings in our maternal early-pregnancy glucose concentrations was found for CpGs identified in a previous study on adult type 2 diabetes.

CONCLUSIONS

Maternal early-pregnancy glucose concentrations, but not insulin concentrations, were associated with DNA methylation at one CpG each in the subgroups of normal weight and of overweight or obese women. No associations were present in the full group. The role of these CpGs in mechanisms underlying offspring health outcomes needs further study. Future studies should replicate our results in larger samples with early-pregnancy information on maternal fasting glucose metabolism.

The Use of Genomics to Drive Kidney Disease Drug Discovery and Development

As opposed to diseases such as cancer, autoimmune disease, and diabetes, identifying drugs to treat CKD has proven significantly more challenging. Over the past 2 decades, new potential therapeutic targets have been identified as genetically altered proteins involved in rare monogenetic kidney diseases. Other possible target genes have been implicated through common genetic polymorphisms associated with CKD in the general population. Significant challenges remain before translating these genetic insights into clinical therapies for CKD. This paper will discuss how genetic variants may be leveraged to develop drugs and will especially focus on those genes associated with CKD to exemplify the value and challenges in including genetic information in the drug development pipeline.

Cardiometabolic Risk Factors for COVID-19 Susceptibility and Severity: A Mendelian Randomization Analysis

IMPORTANCE

Early epidemiological studies report associations of diverse cardiometabolic conditions especially body mass index (BMI), with COVID-19 susceptibility and severity, but causality has not been established. Identifying causal risk factors is critical to inform preventive strategies aimed at modifying disease risk.

OBJECTIVE

We sought to evaluate the causal associations of cardiometabolic conditions with COVID-19 susceptibility and severity.

DESIGN

Two-sample Mendelian Randomization (MR) Study.

SETTING

Population-based cohorts that contributed to the genome-wide association study (GWAS) meta-analysis by the COVID-19 Host Genetics Initiative.

PARTICIPANTS

Patients hospitalized with COVID-19 diagnosed by RNA PCR, serologic testing, or clinician diagnosis. Population controls defined as anyone who was not a case in the cohorts. Exposures: Selected genetic variants associated with 17 cardiometabolic diseases, including diabetes, coronary artery disease, stroke, chronic kidney disease, and BMI, at p<5 x 10-8 from published largescale GWAS.

MAIN OUTCOMES

We performed an inverse-variance weighted averages of variant-specific causal estimates for susceptibility, defined as people who tested positive for COVID-19 vs. population controls, and severity, defined as patients hospitalized with COVID-19 vs. population controls, and repeated the analysis for BMI using effect estimates from UKBB. To estimate direct and indirect causal effects of BMI through obesity-related cardiometabolic diseases, we performed pairwise multivariable MR. We used p<0.05/17 exposure/2 outcomes=0.0015 to declare statistical significance.

RESULTS

Genetically increased BMI was causally associated with testing positive for COVID-19 [6,696 cases / 1,073,072 controls; p=6.7 x 10-4, odds ratio and 95% confidence interval 1.08 (1.03, 1.13) per kg/m2] and a higher risk of COVID-19 hospitalization [3,199 cases/897,488 controls; p=8.7 x 10-4, 1.12 (1.04, 1.21) per kg/m2]. In the multivariable MR, the direct effect of BMI was abolished upon conditioning on the effect on type 2 diabetes but persisted when conditioning on the effects on coronary artery disease, stroke, chronic kidney disease, and c-reactive protein. No other cardiometabolic exposures tested were associated with a higher risk of poorer COVID-19 outcomes.

CONCLUSIONS AND RELEVANCE

Genetic evidence supports BMI as a causal risk factor for COVID-19 susceptibility and severity. This relationship may be mediated via type 2 diabetes. Obesity may have amplified the disease burden of the COVID-19 pandemic either single-handedly or through its metabolic consequences.

Impaired renal hemodynamics and glomerular hyperfiltration contribute to hypertension-induced renal injury

Recently, we reported a mutation in γ-adducin (ADD3) was associated with an impaired myogenic response of the afferent arteriole and hypertension-induced chronic kidney disease (CKD) in fawn hooded hypertensive (FHH) rats. However, the mechanisms by which altered renal blood flow (RBF) autoregulation promotes hypertension-induced renal injury remain to be determined. The present study compared the time course of changes in renal hemodynamics and the progression of CKD during the development of DOCA-salt hypertension in FHH 1^BN congenic rats [wild-type (WT)] with an intact myogenic response versus FHH 1^BN Add3KO (Add3KO) rats, which have impaired myogenic response. RBF was well autoregulated in WT rats but not in Add3KO rats. Glomerular capillary pressure rose by 6 versus 14 mmHg in WT versus Add3KO rats when blood pressure increased from 100 to 150 mmHg. After 1 wk of hypertension, glomerular filtration rate increased by 38% and glomerular nephrin expression decreased by 20% in Add3KO rats. Neither were altered in WT rats. Proteinuria doubled in WT rats versus a sixfold increase in Add3KO rats. The degree of renal injury was greater in Add3KO than WT rats after 3 wk of hypertension. RBF, glomerular filtration rate, and glomerular capillary pressure were lower by 20%, 28%, and 19% in Add3KO rats than in WT rats, which was associated with glomerular matrix expansion and loss of capillary filtration area. The results indicated that impaired RBF autoregulation and eutrophic remodeling of preglomerular arterioles increase the transmission of pressure to glomeruli, which induces podocyte loss and accelerates the progression of CKD in hypertensive Add3KO rats.

Integration of GWAS Summary Statistics and Gene Expression Reveals Target Cell Types Underlying Kidney Function Traits

BACKGROUND

Genetic variants identified in genome-wide association studies (GWAS) are often not specific enough to reveal complex underlying physiology. By integrating RNA-seq data and GWAS summary statistics, novel computational methods allow unbiased identification of trait-relevant tissues and cell types.

METHODS

The CKDGen consortium provided GWAS summary data for eGFR, urinary albumin-creatinine ratio (UACR), BUN, and serum urate. Genotype-Tissue Expression Project (GTEx) RNA-seq data were used to construct the top 10% specifically expressed genes for each of 53 tissues followed by linkage disequilibrium (LD) score-based enrichment testing for each trait. Similar procedures were performed for five kidney single-cell RNA-seq datasets from humans and mice and for a microdissected tubule RNA-seq dataset from rat. Gene set enrichment analyses were also conducted for genes implicated in Mendelian kidney diseases.

RESULTS

Across 53 tissues, genes in kidney function-associated GWAS loci were enriched in kidney (P=9.1E-8 for eGFR; P=1.2E-5 for urate) and liver (P=6.8·10^-5 for eGFR). In the kidney, proximal tubule was enriched in humans (P=8.5E-5 for eGFR; P=7.8E-6 for urate) and mice (P=0.0003 for eGFR; P=0.0002 for urate) and confirmed as the primary cell type in microdissected tubules and organoids. Gene set enrichment analysis supported this and showed enrichment of genes implicated in monogenic glomerular diseases in podocytes. A systematic approach generated a comprehensive list of GWAS genes prioritized by cell type-specific expression.

CONCLUSIONS

Integration of GWAS statistics of kidney function traits and gene expression data identified relevant tissues and cell types, as a basis for further mechanistic studies to understand GWAS loci.

Sex-specific associations of insulin resistance with chronic kidney disease and kidney function: a bi-directional Mendelian randomisation study

AIMS/HYPOTHESIS

Reasons for the sexual disparity in chronic kidney disease (CKD) are unclear. To provide insight we contextualised these differences within evolutionary biology, and explored sex-specific effects of insulin resistance because it may have sex-specific effects on the reproductive axis. Impaired kidney function may also cause insulin resistance. We assessed these possibilities using bi-directional, sex-specific, two-sample Mendelian randomisation (MR).

METHODS

Given that fasting insulin, fasting glucose and HbA_1c are related, we used MR-Bayesian model averaging (MR-BMA) to identify the best-fitting model and most influential exposure. Genetic associations with glycaemic traits were obtained from genome-wide association studies (GWAS) in Europeans without diabetes (n = 108,557 for fasting insulin, as a proxy for insulin resistance, and for fasting glucose, n = 123,665 for HbA_1c in the Meta-Analyses of Glucose and Insulin-related traits Consortium [MAGIC]), and applied to GWAS of 480,698 Europeans for overall associations with CKD (cases n = 41,395) and eGFR. We also used sex-specific individual information in white British (179,917 men, 6016 CKD cases; 212,079 women, 5958 CKD cases) from the UK Biobank. Univariable or multivariable MR was used to assess the role of glycaemic trait(s) selected by MR-BMA in CKD and kidney function. Genetic variants predicting eGFR were used to assess the role of kidney function in the most influential exposure(s).

RESULTS

Fasting insulin was selected as the most likely exposure by both overall and sex-specific MR-BMA. It increased CKD in men (OR 7.23 per pmol/l higher fasting insulin [95% CI 2.46, 21.2]) but not in women (OR 1.05 [95% CI 0.21, 5.21]), and reduced eGFR in men (-0.04 [95% CI -0.07, -0.01]) but not in women (0.01 [95% CI -0.02, 0.03]). Genetically predicted eGFR was unrelated to fasting insulin.

CONCLUSIONS/INTERPRETATION

Genetically predicted fasting insulin was sex-specifically associated with CKD and unhealthier kidney function but was not affected by kidney function. Graphical abstract.

The relationship between blood metabolites of the tryptophan pathway and kidney function: a bidirectional Mendelian randomization analysis

Blood metabolites of the tryptophan pathway were found to be associated with kidney function and disease in observational studies. In order to evaluate causal relationship and direction, we designed a study using a bidirectional Mendelian randomization approach. The analyses were based on published summary statistics with study sizes ranging from 1,960 to 133,413. After correction for multiple testing, results provided no evidence of an effect of metabolites of the tryptophan pathway on estimated glomerular filtration rate (eGFR). Conversely, lower eGFR was related to higher levels of four metabolites: C-glycosyltryptophan (effect estimate = − 0.16, 95% confidence interval [CI] (− 0.22; − 0.1); p = 9.2e−08), kynurenine (effect estimate = − 0.18, 95% CI (− 0.25; − 0.11); p = 1.1e−06), 3-indoxyl sulfate (effect estimate = − 0.25, 95% CI (− 0.4; − 0.11); p = 6.3e−04) and indole-3-lactate (effect estimate = − 0.26, 95% CI (− 0.38; − 0.13); p = 5.4e−05). Our study supports that lower eGFR causes higher blood metabolite levels of the tryptophan pathway including kynurenine, C-glycosyltryptophan, 3-indoxyl sulfate, and indole-3-lactate. These findings aid the notion that metabolites of the tryptophan pathway are a consequence rather than a cause of reduced eGFR. Further research is needed to specifically examine relationships with respect to chronic kidney disease (CKD) progression among patients with existing CKD.

Genomewide Association Study of Platelet Reactivity and Cardiovascular Response in Patients Treated With Clopidogrel: A Study by the International Clopidogrel Pharmacogenomics Consortium

Antiplatelet response to clopidogrel shows wide variation, and poor response is correlated with adverse clinical outcomes. CYP2C19 loss‐of‐function alleles play an important role in this response, but account for only a small proportion of variability in response to clopidogrel. An aim of the International Clopidogrel Pharmacogenomics Consortium (ICPC) is to identify other genetic determinants of clopidogrel pharmacodynamics and clinical response. A genomewide association study (GWAS) was performed using DNA from 2,750 European ancestry individuals, using adenosine diphosphate‐induced platelet reactivity and major cardiovascular and cerebrovascular events as outcome parameters. GWAS for platelet reactivity revealed a strong signal for CYP2C19*2 (P value = 1.67e−33). After correction for CYP2C19*2 no other single‐nucleotide polymorphism reached genomewide significance. GWAS for a combined clinical end point of cardiovascular death, myocardial infarction, or stroke (5.0% event rate), or a combined end point of cardiovascular death or myocardial infarction (4.7% event rate) showed no significant results, although in coronary artery disease, percutaneous coronary intervention, and acute coronary syndrome subgroups, mutations in SCOS5P1, CDC42BPA, and CTRAC1 showed genomewide significance (lowest P values: 1.07e−09, 4.53e−08, and 2.60e−10, respectively). CYP2C19*2 is the strongest genetic determinant of on‐clopidogrel platelet reactivity. We identified three novel associations in clinical outcome subgroups, suggestive for each of these outcomes.

Mendelian Randomization Analysis Reveals a Causal Effect of Urinary Sodium/Urinary Creatinine Ratio on Kidney Function in Europeans

Salt restriction was recommended in clinical practice guideline for chronic kidney disease (CKD) treatment, but its effect on kidney outcomes remains conflicting. We aimed to test the causal effect of salt intake, using estimated 24-h sodium excretion from spot urinary sodium/urinary creatinine (UNa/UCr) ratio as a surrogate, on renal function using two-sample Mendelian randomization (MR). Genetic instruments for UNa/UCr were derived from a recent genome-wide association study of 218,450 European-descent individuals in the UK Biobank. Kidney outcomes were creatinine-based estimated glomerular filtration rate (eGFRcrea) (N = 567,460) and CKD (eGFRcrea < 60 ml/min/1.73 m², N cases = 41,395, N controls = 439,303) from the CKDGen consortium. Cystatin C-based eGFR (eGFRcys) and eGFRcrea single-nucleotide polymorphisms associated with blood urea nitrogen (BUN) were used for sensitivity analyses. MR revealed a causal effect of UNa/UCr on higher eGFRcrea [β = 0.14, unit change in log ml/min/1.73 m² per UNa/UCr ratio; 95% confidence interval (CI) = 0.07 - 0.20, P = 2.15 × 10^-5] and a protective effect against CKD risk (odds ratio = 0.24, 95% CI = 0.14 to 0.41, P = 1.20 × 10^-7). The MR findings were confirmed by MR-Egger regression, weighted median MR, and mode estimate MR, with less evidence of existence of horizontal pleiotropy. Consistent positive causal effect of UNa/UCr on eGFRcys was also detected. On the other hand, bidirectional MR suggested inconclusive results of CKD, eGFRcrea, eGFRcrea (BUN associated), and eGFRcys on UNa/UCr. The average 24-h sodium excretion was estimated to be approximately 2.6 g per day for women and 3.7 g per day for men. This study provides evidence that sodium excretion, well above the recommendation of <2 g per day of sodium intake, might not have a harmful effect on kidney function. Clinical trials are warranted to evaluate the sodium restriction target on kidney function.

Complexities of Understanding Function from CKD-Associated DNA Variants

Genome-wide association studies (GWASs) have facilitated the unbiased discovery of hundreds of genomic loci associated with CKD and kidney function. The vast majority of disease-associated DNA variants are noncoding. Those that are causal in CKD pathogenesis likely modulate transcription of target genes in a cell type-specific manner. To gain novel biological insights into mechanisms driving the development of CKD, the causal variants (which are usually not the most significant variant reported in a GWAS), their target genes, and causal cell types need to be identified. This functional validation requires a large number of new data sets, complex bioinformatics analyses, and experimental cellular and in vivo studies. Here, we review the basic principles and some of the current approaches being leveraged to assign functional significance to a genotype-phenotype association.

Genetics of diabetes mellitus and diabetes complications

Diabetes is one of the fastest growing diseases worldwide, projected to affect 693 million adults by 2045. Devastating macrovascular complications (cardiovascular disease) and microvascular complications (such as diabetic kidney disease, diabetic retinopathy and neuropathy) lead to increased mortality, blindness, kidney failure and an overall decreased quality of life in individuals with diabetes. Clinical risk factors and glycaemic control alone cannot predict the development of vascular complications; numerous genetic studies have demonstrated a clear genetic component to both diabetes and its complications. Early research aimed at identifying genetic determinants of diabetes complications relied on familial linkage analysis suited to strong-effect loci, candidate gene studies prone to false positives, and underpowered genome-wide association studies limited by sample size. The explosion of new genomic datasets, both in terms of biobanks and aggregation of worldwide cohorts, has more than doubled the number of genetic discoveries for both diabetes and diabetes complications. We focus herein on genetic discoveries for diabetes and diabetes complications, empowered primarily through genome-wide association studies, and emphasize the gaps in research for taking genomic discovery to the next level.

Plasma Metabolomics Identifies Markers of Impaired Renal Function: A Meta-analysis of 3089 Persons with Type 2 Diabetes

CONTEXT

There is a need for novel biomarkers and better understanding of the pathophysiology of diabetic kidney disease.

OBJECTIVE

To investigate associations between plasma metabolites and kidney function in people with type 2 diabetes (T2D).

DESIGN

3089 samples from individuals with T2D, collected between 1999 and 2015, from 5 independent Dutch cohort studies were included. Up to 7 years follow-up was available in 1100 individuals from 2 of the cohorts.

MAIN OUTCOME MEASURES

Plasma metabolites (n = 149) were measured by nuclear magnetic resonance spectroscopy. Associations between metabolites and estimated glomerular filtration rate (eGFR), urinary albumin-to-creatinine ratio (UACR), and eGFR slopes were investigated in each study followed by random effect meta-analysis. Adjustments included traditional cardiovascular risk factors and correction for multiple testing.

RESULTS

In total, 125 metabolites were significantly associated (PFDR = 1.5×10-32 - 0.046; β = -11.98-2.17) with eGFR. Inverse associations with eGFR were demonstrated for branched-chain and aromatic amino acids (AAAs), glycoprotein acetyls, triglycerides (TGs), lipids in very low-density lipoproteins (VLDL) subclasses, and fatty acids (PFDR < 0.03). We observed positive associations with cholesterol and phospholipids in high-density lipoproteins (HDL) and apolipoprotein A1 (PFDR < 0.05). Albeit some metabolites were associated with UACR levels (P < 0.05), significance was lost after correction for multiple testing. Tyrosine and HDL-related metabolites were positively associated with eGFR slopes before adjustment for multiple testing (PTyr = 0.003; PHDLrelated < 0.05), but not after.

CONCLUSIONS

This study identified metabolites associated with impaired kidney function in T2D, implying involvement of lipid and amino acid metabolism in the pathogenesis. Whether these processes precede or are consequences of renal impairment needs further investigation.

GADL1 is a multifunctional decarboxylase with tissue-specific roles in β-alanine and carnosine production

Carnosine and related β-alanine-containing peptides are believed to be important antioxidants, pH buffers, and neuromodulators. However, their biosynthetic routes and therapeutic potential are still being debated. This study describes the first animal model lacking the enzyme glutamic acid decarboxylase-like 1 (GADL1). We show that Gadl1^-/- mice are deficient in β-alanine, carnosine, and anserine, particularly in the olfactory bulb, cerebral cortex, and skeletal muscle. Gadl1^-/- mice also exhibited decreased anxiety, increased levels of oxidative stress markers, alterations in energy and lipid metabolism, and age-related changes. Examination of the GADL1 active site indicated that the enzyme may have multiple physiological substrates, including aspartate and cysteine sulfinic acid. Human genetic studies show strong associations of the GADL1 locus with plasma levels of carnosine, subjective well-being, and muscle strength. Together, this shows the multifaceted and organ-specific roles of carnosine peptides and establishes Gadl1 knockout mice as a versatile model to explore carnosine biology and its therapeutic potential.

Genetic architecture of host proteins interacting with SARS-CoV-2

Strategies to develop therapeutics for SARS-CoV-2 infection may be informed by experimental identification of viral-host protein interactions in cellular assays and measurement of host response proteins in COVID-19 patients. Identification of genetic variants that influence the level or activity of these proteins in the host could enable rapid 'in silico' assessment in human genetic studies of their causal relevance as molecular targets for new or repurposed drugs to treat COVID-19. We integrated large-scale genomic and aptamer-based plasma proteomic data from 10,708 individuals to characterize the genetic architecture of 179 host proteins reported to interact with SARS-CoV-2 proteins or to participate in the host response to COVID-19. We identified 220 host DNA sequence variants acting in cis (MAF 0.01-49.9%) and explaining 0.3-70.9% of the variance of 97 of these proteins, including 45 with no previously known protein quantitative trait loci (pQTL) and 38 encoding current drug targets. Systematic characterization of pQTLs across the phenome identified protein-drug-disease links, evidence that putative viral interaction partners such as MARK3 affect immune response, and establish the first link between a recently reported variant for respiratory failure of COVID-19 patients at the ABO locus and hypercoagulation, i.e. maladaptive host response. Our results accelerate the evaluation and prioritization of new drug development programmes and repurposing of trials to prevent, treat or reduce adverse outcomes. Rapid sharing and dynamic and detailed interrogation of results is facilitated through an interactive webserver ( https://omicscience.org/apps/covidpgwas/ ).

Genetic drug target validation using Mendelian randomisation

Mendelian randomisation (MR) analysis is an important tool to elucidate the causal relevance of environmental and biological risk factors for disease. However, causal inference is undermined if genetic variants used to instrument a risk factor also influence alternative disease-pathways (horizontal pleiotropy). Here we report how the 'no horizontal pleiotropy assumption' is strengthened when proteins are the risk factors of interest. Proteins are typically the proximal effectors of biological processes encoded in the genome. Moreover, proteins are the targets of most medicines, so MR studies of drug targets are becoming a fundamental tool in drug development. To enable such studies, we introduce a mathematical framework that contrasts MR analysis of proteins with that of risk factors located more distally in the causal chain from gene to disease. We illustrate key model decisions and introduce an analytical framework for maximising power and evaluating the robustness of analyses.

Elevated body mass index is associated with an increased risk of infectious disease admissions and mortality: a mendelian randomization study

OBJECTIVE

The effect of body mass index (BMI) on the risk of infectious diseases admissions and mortality is unclear and is difficult to study given the risks of confounding variables.

METHODS

We used genome-wide association studies (GWASs) with mendelian randomization (MR) to obtain causal inference of BMI on the following infectious diseases outcomes: hospital admissions for pneumonia, sepsis, urinary tract infections, skin and soft tissue infections (SSTIs) or all-cause infections. For patients with pneumonia and sepsis, we also analysed their 28-day and 90-day mortalities. The UK Biobank (UKB) cohort (n > 500 000) provided data for GWASs on infectious diseases. The GIANT consortium (n = 681 265) GWAS was used to identify single-nucleotide polymorphisms (SNPs) associated with BMI.

RESULTS

Genetically increased BMI, by one standard deviation, was associated with higher rates of admission due to all infectious disease. The effect was most important for SSTIs (OR: 1.11, 95%CI: 1.09, 1.12). Increasing BMI by one standard deviation was associated with higher pneumonia mortality, especially at 28 days (OR: 1.03, 95%CI: 1.01, 1.05). BMI was not clearly associated with sepsis mortality, although interpretation of the results was limited by a small sample size. There were consistent findings in sensitivity analysis performed by removing highly pleiotropic SNPs and multivariate MR including type-2 diabetes mellitus, estimated glomerular filtration rate, high-density lipoprotein, educational attainment, and a history of smoking.

CONCLUSIONS

Increased BMI was associated with increased risk of admission for infectious disease and mortality. While the pathophysiology behind this phenomenon remains unknown, increasing BMI may influence immune dysregulation.

A genome-wide association study implicates multiple mechanisms influencing raised urinary albumin-creatinine ratio

Raised albumin–creatinine ratio (ACR) is an indicator of microvascular damage and renal disease. We aimed to identify genetic variants associated with raised ACR and study the implications of carrying multiple ACR-raising alleles with metabolic and vascular-related disease. We performed a genome-wide association study of ACR using 437 027 individuals from the UK Biobank in the discovery phase, 54 527 more than previous studies, and followed up our findings in independent studies. We identified 62 independent associations with ACR across 56 loci (P < 5 × 10^–8), of which 20 were not previously reported. Pathway analyses and the identification of 20 of the 62 variants (at r² > 0.8) coinciding with signals for at least 16 related metabolic and vascular traits, suggested multiple pathways leading to raised ACR levels. After excluding variants at the CUBN locus, known to alter ACR via effects on renal absorption, an ACR genetic risk score was associated with a higher risk of hypertension, and less strongly, type 2 diabetes and stroke. For some rare genotype combinations at the CUBN locus, most individuals had ACR levels above the microalbuminuria clinical threshold. Contrary to our hypothesis, individuals carrying more CUBN ACR-raising alleles, and above the clinical threshold, had a higher frequency of vascular disease. The CUBN allele effects on ACR were twice as strong in people with diabetes—a result robust to an optimization-algorithm approach to simulating interactions, validating previously reported gene–diabetes interactions (P ≤ 4 × 10^–5). In conclusion, a variety of genetic mechanisms and traits contribute to variation in ACR.

Genetic Risk Scores for Complex Disease Traits in Youth

BACKGROUND

For most disease-related traits the magnitude of the contribution of genetic factors in adolescents remains unclear.

METHODS

Twenty continuous traits related to anthropometry, cardiovascular and renal function, metabolism, and inflammation were selected from the ongoing prospective Tracking Adolescents' Individual Lives Survey cohort in the Netherlands with measurements of up to 5 waves from age 11 to 22 years (n=1354, 47.6% males) and all traits available at the third wave (mean age [SD]=16.22 [0.66]). For each trait, unweighted and weighted genetic risk scores (GRSs) were generated based on significantly associated single nucleotide polymorphisms identified from literature. The variance explained by the GRSs in adolescents were estimated by linear regression after adjustment for covariates.

RESULTS

Except for ALT (alanine transaminase), all GRSs were significantly associated with their traits. The trait variance explained by the GRSs was highest for lipoprotein[a] (39.59%) and varied between 0.09% (ALT) and 18.49% (LDL [low-density lipoprotein]) for the other traits. For most traits, the variances explained in adolescents were comparable with or slightly smaller than those in adults. Significant increases of trait levels (except ALT) and increased risks for overweight/obesity (odds ratio, 6.41 [95% CI, 2.95-15.56]) and hypertension (odds ratio, 2.86 [95% CI, 1.39-6.17]) were found in individuals in the top GRS decile compared with those at the bottom decile.

CONCLUSIONS

Variances explained by adult-based GRSs for disease-related traits in adolescents, although still relatively modest, were comparable with or slightly smaller than in adults offering promise for improved risk prediction at early ages.

The role of testosterone in chronic kidney disease and kidney function in men and women: a bi-directional Mendelian randomization study in the UK Biobank

BACKGROUND

Chronic kidney disease (CKD) has an apparent sex disparity, with a more rapid progress in men than in women. Whether the well-established sex-specific evolutionary biology trade-off between reproduction and longevity might inform CKD has not previously been considered. Relevant evidence from randomized controlled trials (RCTs) is not available.

METHODS

We used a bi-directional Mendelian randomization study to obtain unconfounded estimates using the UK Biobank. Single nucleotide polymorphisms (SNPs) that strongly (p value < 5 × 10^-8) predicted testosterone in a sex-specific manner were applied to 179,916 white British men (6016 CKD cases) and 212,079 white British women (5958 CKD cases) to obtain sex-specific associations with CKD, albuminuria, and estimated glomerular filtration rate (eGFR). We also used multivariable MR to control for sex hormone binding globulin (SHBG). For validation, we similarly examined their role in hemoglobin and high-density lipoprotein cholesterol (HDL-c). We also assessed the role of kidney function in serum testosterone, by applying eGFR-related SNPs to testosterone in the UK Biobank.

RESULTS

Genetically predicted testosterone was associated with CKD in men (odds ratio (OR) for bioavailable testosterone 1.17 per standard deviation, 95% confidence interval (CI) 1.03 to 1.33) based on 125 SNPs but not in women (OR 1.02, 95% CI 0.92 to 1.14 for total testosterone) based on 254 SNPs. Multivariable MR allowing for SHBG showed consistent patterns. Genetically predicted bioavailable testosterone in men and women and genetically predicted total testosterone in women increased hemoglobin and lowered HDL-c, as seen in RCTs. Genetically predicted eGFR was not related to serum testosterone in men or in women.

CONCLUSIONS

Genetically predicted testosterone was associated with CKD and worse kidney function in men, whilst not affected by kidney function. Identifying drivers of testosterone and the underlying pathways could provide new insights into CKD prevention and treatment.

Genetics of Biochemical Phenotypes

Biomarkers diagnose, predict or assess the risk of disease, and studies of the effects of genetic variation on biomarker phenotypes in the general population complement studies on patients diagnosed with disease. This paper traces the evolution of studies on biomarker genetics over the past 40 years through examples drawn from the work of Professor Martin and his colleagues.

A bidirectional Mendelian randomization study supports causal effects of kidney function on blood pressure

Blood pressure and kidney function have a bidirectional relation. Hypertension has long been considered as a risk factor for kidney function decline. However, whether intensive blood pressure control could promote kidney health has been uncertain. The kidney is known to have a major role in affecting blood pressure through sodium extraction and regulating electrolyte balance. This bidirectional relation makes causal inference between these two traits difficult. Therefore, to examine the causal relations between these two traits, we performed two-sample Mendelian randomization analyses using summary statistics of large-scale genome-wide association studies. We selected genetic instruments more likely to be specific for kidney function using meta-analyses of complementary kidney function biomarkers (glomerular filtration rate estimated from serum creatinine [eGFRcr], and blood urea nitrogen from the CKDGen Consortium). Systolic and diastolic blood pressure summary statistics were from the International Consortium for Blood Pressure and UK Biobank. Significant evidence supported the causal effects of higher kidney function on lower blood pressure. Based on the mode-based Mendelian randomization method, the effect estimates for one standard deviation (SD) higher in log-transformed eGFRcr was -0.17 SD unit (95 % confidence interval: -0.09 to -0.24) in systolic blood pressure and -0.15 SD unit (95% confidence interval: -0.07 to -0.22) in diastolic blood pressure. In contrast, the causal effects of blood pressure on kidney function were not statistically significant. Thus, our results support causal effects of higher kidney function on lower blood pressure and suggest preventing kidney function decline can reduce the public health burden of hypertension.

Kidney Disease, Race, and GFR Estimation

Assessment of GFR is central to clinical practice, research, and public health. Current Kidney Disease Improving Global Outcomes guidelines recommend measurement of serum creatinine to estimate GFR as the initial step in GFR evaluation. Serum creatinine is influenced by creatinine metabolism as well as GFR; hence, all equations to estimate GFR from serum creatinine include surrogates for muscle mass, such as age, sex, race, height, or weight. The guideline-recommended equation in adults (the 2009 Chronic Kidney Disease Epidemiology Collaboration creatinine equation) includes a term for race (specified as black versus nonblack), which improves the accuracy of GFR estimation by accounting for differences in non-GFR determinants of serum creatinine by race in the study populations used to develop the equation. In that study, blacks had a 16% higher average measured GFR compared with nonblacks with the same age, sex, and serum creatinine. The reasons for this difference are only partly understood, and the use of race in GFR estimation has limitations. Some have proposed eliminating the race coefficient, but this would induce a systematic underestimation of measured GFR in blacks, with potential unintended consequences at the individual and population levels. We propose a more cautious approach that maintains and improves accuracy of GFR estimates and avoids disadvantaging any racial group. We suggest full disclosure of use of race in GFR estimation, accommodation of those who decline to identify their race, and shared decision making between health care providers and patients. We also suggest mindful use of cystatin C as a confirmatory test as well as clearance measurements. It would be preferable to avoid specification of race in GFR estimation if there was a superior, evidence-based substitute. The goal of future research should be to develop more accurate methods for GFR estimation that do not require use of race or other demographic characteristics.

A Mendelian Randomization Study Provides Evidence That Adiposity and Dyslipidemia Lead to Lower Urinary Albumin-to-Creatinine Ratio, a Marker of Microvascular Function

Urinary albumin-to-creatinine ratio (ACR) is a marker of diabetic nephropathy and microvascular damage. Metabolic-related traits are observationally associated with ACR, but their causal role is uncertain. Here, we confirmed ACR as a marker of microvascular damage and tested whether metabolic-related traits have causal relationships with ACR. The association between ACR and microvascular function (responses to acetylcholine [ACH] and sodium nitroprusside) was tested in the SUMMIT study. Two-sample Mendelian randomization (MR) was used to infer the causal effects of 11 metabolic risk factors, including glycemic, lipid, and adiposity traits, on ACR. MR was performed in up to 440,000 UK Biobank and 54,451 CKDGen participants. ACR was robustly associated with microvascular function measures in SUMMIT. Using MR, we inferred that higher triglyceride (TG) and LDL cholesterol (LDL-C) levels caused elevated ACR. A 1 SD higher TG and LDL-C level caused a 0.062 (95% CI 0.040, 0.083) and a 0.026 (95% CI 0.008, 0.044) SD higher ACR, respectively. There was evidence that higher body fat and visceral body fat distribution caused elevated ACR, while a metabolically "favorable adiposity" phenotype lowered ACR. ACR is a valid marker for microvascular function. MR suggested that seven traits have causal effects on ACR, highlighting the role of adiposity-related traits in causing lower microvascular function.

Using Genetic and Species Diversity to Tackle Kidney Disease

Progress in the identification of causal genes and understanding of the mechanism underlying kidney disease is hindered by the almost exclusive use of a few animal models with restrictive monogenic backgrounds that may be more resistant to kidney disease compared with humans and, therefore, poor models. Exploring the large genetic diversity in classical animal models, such as mice and rats, and leveraging species diversity will allow us to use the genetic advantages of zebrafish, Drosophila, and other species, to develop both new animal models that are more relevant to the study of human kidney disease and potential therapies.

Results from the German Chronic Kidney Disease (GCKD) study support association of relative telomere length with mortality in a large cohort of patients with moderate chronic kidney disease

Telomere length is known to be inversely associated with aging and has been proposed as a marker for aging-related diseases. Telomere attrition can be accelerated by oxidative stress and inflammation, both commonly present in patients with chronic kidney disease. Here, we investigated whether relative telomere length is associated with mortality in a large cohort of patients with chronic kidney disease stage G3 and A1-3 or G1-2 with overt proteinuria (A3) at enrollment. Relative telomere length was quantified in peripheral blood by a quantitative PCR method in 4,955 patients from the GCKD study, an ongoing prospective observational cohort. Complete four-year follow-up was available from 4,926 patients in whom we recorded 354 deaths. Relative telomere length was a strong and independent predictor of all-cause mortality. Each decrease of 0.1 relative telomere length unit was highly associated with a 14% increased risk of death (hazard ratio1.14 [95% confidence interval 1.06-1.22]) in a model adjusted for age, sex, baseline eGFR, urine albumin/creatinine ratio, diabetes mellitus, prevalent cardiovascular disease, LDL-cholesterol, HDL-cholesterol, smoking, body mass index, systolic and diastolic blood pressure, C-reactive protein and serum albumin. This translated to a 75% higher risk for those in the lowest compared to the highest quartile of relative telomere length. The association was mainly driven by 117 cardiovascular deaths (1.20 [1.05-1.35]) as well as 67 deaths due to infections (1.27 [1.07-1.50]). Thus, our findings support an association of shorter telomere length with all-cause mortality, cardiovascular mortality and death due to infections in patients with moderate chronic kidney disease.

Applying single-cell technologies to clinical pathology: progress in nephropathology

Cells represent the basic building blocks of living organisms. Accurate characterisation of cellular phenotype, intercellular signalling networks, and the spatial organisation of cells within organs is crucial to deliver a better understanding of the processes underpinning physiology, and the perturbations that lead to disease. Single-cell methodologies have increased rapidly in scale and scope in recent years and are set to generate important insights into human disease. Here, we review current practices in nephropathology, which are dominated by relatively simple morphological descriptions of tissue biopsies based on their appearance using light microscopy. Bulk transcriptomics have more recently been used to explore glomerular and tubulointerstitial kidney disease, renal cancer, and the responses to injury and alloimmunity in kidney transplantation, generating novel disease insights and prognostic biomarkers. These studies set the stage for single-cell transcriptomic approaches that reveal cell-type-specific gene expression patterns in health and disease. These technologies allow genome-wide disease susceptibility genes to be interpreted with the knowledge of the specific cell populations within organs that express them, identifying candidate cell types for further study. Single-cell technologies are also moving beyond assaying individual cellular transcriptomes, to measuring the epigenetic landscape of single cells. Single-cell antigen-receptor gene sequencing also enables specific T- and B-cell clones to be tracked in different tissues and disease states. In the coming years these rich 'multi-omic' descriptions of kidney disease will enable histopathological descriptions to be comprehensively integrated with molecular phenotypes, enabling better disease classification and prognostication and the application of personalised treatment strategies. © 2020 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

Dnmt3a and Dnmt3b-Decommissioned Fetal Enhancers are Linked to Kidney Disease

BACKGROUND

Cytosine methylation is an epigenetic mark that dictates cell fate and response to stimuli. The timing and establishment of methylation logic during kidney development remains unknown. DNA methyltransferase 3a and 3b are the enzymes capable of establishing de novo methylation.

METHODS

We generated mice with genetic deletion of Dnmt3a and Dnmt3b in nephron progenitor cells (Six2CreDnmt3a/3b) and kidney tubule cells (KspCreDnmt3a/3b). We characterized KspCreDnmt3a/3b mice at baseline and after injury. Unbiased omics profiling, such as whole genome bisulfite sequencing, reduced representation bisulfite sequencing and RNA sequencing were performed on whole-kidney samples and isolated renal tubule cells.

RESULTS

KspCreDnmt3a/3b mice showed no obvious morphologic and functional alterations at baseline. Knockout animals exhibited increased resistance to cisplatin-induced kidney injury, but not to folic acid-induced fibrosis. Whole-genome bisulfite sequencing indicated that Dnmt3a and Dnmt3b play an important role in methylation of gene regulatory regions that act as fetal-specific enhancers in the developing kidney but are decommissioned in the mature kidney. Loss of Dnmt3a and Dnmt3b resulted in failure to silence developmental genes. We also found that fetal-enhancer regions methylated by Dnmt3a and Dnmt3b were enriched for kidney disease genetic risk loci. Methylation patterns of kidneys from patients with CKD showed defects similar to those in mice with Dnmt3a and Dnmt3b deletion.

CONCLUSIONS

Our results indicate a potential locus-specific convergence of genetic, epigenetic, and developmental elements in kidney disease development.

Mechanisms of cognitive dysfunction in CKD

Cognitive impairment is an increasingly recognized major cause of chronic disability and is commonly found in patients with chronic kidney disease (CKD). Knowledge of the relationship between kidney dysfunction and impaired cognition may improve our understanding of other forms of cognitive dysfunction. Patients with CKD are at an increased risk (compared with the general population) of both dementia and its prodrome, mild cognitive impairment (MCI), which are characterized by deficits in executive functions, memory and attention. Brain imaging in patients with CKD has revealed damage to white matter in the prefrontal cortex and, in animal models, in the subcortical monoaminergic and cholinergic systems, accompanied by widespread macrovascular and microvascular damage. Unfortunately, current interventions that target cardiovascular risk factors (such as anti-hypertensive drugs, anti-platelet agents and statins) seem to have little or no effect on CKD-associated MCI, suggesting that the accumulation of uraemic neurotoxins may be more important than disturbed haemodynamic factors or lipid metabolism in MCI pathogenesis. Experimental models show that the brain monoaminergic system is susceptible to uraemic neurotoxins and that this system is responsible for the altered sleep pattern commonly observed in patients with CKD. Neural progenitor cells and the glymphatic system, which are important in Alzheimer disease pathogenesis, may also be involved in CKD-associated MCI. More detailed study of CKD-associated MCI is needed to fully understand its clinical relevance, underlying pathophysiology, possible means of early diagnosis and prevention, and whether there may be novel approaches and potential therapies with wider application to this and other forms of cognitive decline.

The genetic architecture of membranous nephropathy and its potential to improve non-invasive diagnosis

Membranous Nephropathy (MN) is a rare autoimmune cause of kidney failure. Here we report a genome-wide association study (GWAS) for primary MN in 3,782 cases and 9,038 controls of East Asian and European ancestries. We discover two previously unreported loci, NFKB1 (rs230540, OR = 1.25, P = 3.4 × 10-12) and IRF4 (rs9405192, OR = 1.29, P = 1.4 × 10-14), fine-map the PLA2R1 locus (rs17831251, OR = 2.25, P = 4.7 × 10-103) and report ancestry-specific effects of three classical HLA alleles: DRB1*1501 in East Asians (OR = 3.81, P = 2.0 × 10-49), DQA1*0501 in Europeans (OR = 2.88, P = 5.7 × 10-93), and DRB1*0301 in both ethnicities (OR = 3.50, P = 9.2 × 10-23 and OR = 3.39, P = 5.2 × 10-82, respectively). GWAS loci explain 32% of disease risk in East Asians and 25% in Europeans, and correctly re-classify 20-37% of the cases in validation cohorts that are antibody-negative by the serum anti-PLA2R ELISA diagnostic test. Our findings highlight an unusual genetic architecture of MN, with four loci and their interactions accounting for nearly one-third of the disease risk.

Protecting the kidney in systemic lupus erythematosus: from diagnosis to therapy

Lupus nephritis (LN) is a common manifestation of systemic lupus erythematosus that can lead to irreversible renal impairment. Although the prognosis of LN has improved substantially over the past 50 years, outcomes have plateaued in the USA in the past 20 years as immunosuppressive therapies have failed to reverse disease in more than half of treated patients. This failure might reflect disease complexity and heterogeneity, as well as social and economic barriers to health-care access that can delay intervention until after damage has already occurred. LN progression is still poorly understood and involves multiple cell types and both immune and non-immune mechanisms. Single-cell analysis of intrinsic renal cells and infiltrating cells from patients with LN is a new approach that will help to define the pathways of renal injury at a cellular level. Although many new immune-modulating therapies are being tested in the clinic, the development of therapies to improve regeneration of the injured kidney and to prevent fibrosis requires a better understanding of the mechanisms of LN progression. This mechanistic understanding, together with the development of clinical measures to evaluate risk and detect early disease and better access to expert health-care providers, should improve outcomes for patients with LN.

White Blood Cells and Blood Pressure: A Mendelian Randomization Study

BACKGROUND

High blood pressure (BP) is a risk factor for cardiovascular morbidity and mortality. While BP is regulated by the function of kidney, vasculature, and sympathetic nervous system, recent experimental data suggest that immune cells may play a role in hypertension.

METHODS

We studied the relationship between major white blood cell types and blood pressure in the UK Biobank population and used Mendelian randomization (MR) analyses using the ≈750 000 UK-Biobank/International Consortium of Blood Pressure-Genome-Wide Association Studies to examine which leukocyte populations may be causally linked to BP.

RESULTS

A positive association between quintiles of lymphocyte, monocyte, and neutrophil counts, and increased systolic BP, diastolic BP, and pulse pressure was observed (eg, adjusted systolic BP mean±SE for 1st versus 5th quintile respectively: 140.13±0.08 versus 141.62±0.07 mm Hg for lymphocyte, 139.51±0.08 versus 141.84±0.07 mm Hg for monocyte, and 137.96±0.08 versus 142.71±0.07 mm Hg for neutrophil counts; all P<10^-50). Using 121 single nucleotide polymorphisms in MR, implemented through the inverse-variance weighted approach, we identified a potential causal relationship of lymphocyte count with systolic BP and diastolic BP (causal estimates: 0.69 [95% CI, 0.19-1.20] and 0.56 [95% CI, 0.23-0.90] of mm Hg per 1 SD genetically elevated lymphocyte count, respectively), which was directionally concordant to the observational findings. These inverse-variance weighted estimates were consistent with other robust MR methods. The exclusion of rs3184504 SNP in the SH2B3 locus attenuated the magnitude of the signal in some of the MR analyses. MR in the reverse direction found evidence of positive effects of BP indices on counts of monocytes, neutrophils, and eosinophils but not lymphocytes or basophils. Subsequent MR testing of lymphocyte count in the context of genetic correlation with renal function or resting and postexercise heart rate demonstrated a positive association of lymphocyte count with urine albumin-to-creatinine ratio.

CONCLUSIONS

Observational and genetic analyses demonstrate a concordant, positive and potentially causal relationship of lymphocyte count with systolic BP and diastolic BP.

Hypothyroidism and Kidney Function: A Mendelian Randomization Study

Background: Uncertainty in the mechanism and directionality of observational associations between thyroid function and kidney function may be addressed by genetic analysis with an instrumental variable method termed bidirectional Mendelian randomization (MR). Methods: In the Women's Genome Health Study (WGHS), observational associations between thyroid measures and kidney function were evaluated. Genetic instruments for MR were from recent genome-wide association studies (GWAS) of hypothyroidism, thyrotropin (TSH), and free thyroxine (fT4) concentrations within the reference range, thyroid peroxidase antibodies (TPOAb), estimated glomerular filtration rate from creatinine (eGFR_crea), eGFR from cystatin C (eGFR_cys), and chronic kidney disease (CKD). In WGHS individual-level data, these instruments were used for bidirectional MR between thyroid (N = 3336) and kidney (N = 23,186) functions. To increase power, MR was also performed using GWAS summary statistics from the Chronic Kidney Disease Genetics Consortium (CKDGen) for eGFR_crea (N = 567,460), eGFR_cys (N = 24,063), CKD [N(total) = 480,698, N(cases) = 41,395], and urinary albumin/creatinine ratio (UACR/N = 54,450). Results: In the WGHS, hypothyroidism was observationally associated with decreased eGFR_crea [beta (standard error, SE): -0.024 (0.009) ln(mL/min/1.73 m²), p = 0.01]. By MR, hypothyroidism was associated with decreased eGFR_crea in the WGHS [beta (SE): -0.007 (0.002) per doubled odds hypothyroidism, p = 1.7 × 10^-3] and in CKDGen [beta (SE): -0.004 (0.0005), p = 2.0 × 10^-22], and robust to sensitivity analysis. Hypothyroidism was also associated by MR with increased CKD in CKDGen (odds ratio, OR [confidence interval, CI]: 1.05 [1.03-1.08], p = 3.3 × 10^-5), but not in the WGHS (OR [CI]: 1.02 [0.95-1.10], p = 0.57). Increased TSH within the reference range had an MR association with increased eGFR_crea in the WGHS [beta (SE): -0.018 (0.007) ln(mL/min/1.73 m²)/standard deviation, SD, p = 6.5 × 10^-3] and CKDGen [beta (SE): -0.008 (0.001) ln(mL/min/1.73 m²)/SD, p = 6.8 × 10^-17], and with CKD in CKDGen (OR [CI]: 1.10 [1.04-1.15], p = 3.1 × 10^-4). There were no MR associations of hypothyroidism or TSH with eGFR_cys or UACR, and MR associations of fT4 in the reference range with kidney function were inconsistent in both the WGHS and CKDGen. However, by MR in CKDGen, TPOAb were robustly associated with decreased eGFR_crea [beta (SE): -0.041 (0.009), p = 6.2 × 10^-6] and decreased eGFR_cys [beta (SE): -0.294 (0.065), p = 6.2 × 10^-6]. TPOAb were less robustly associated with CKD but not associated with UACR. In reverse MR in the WGHS, kidney function was not consistently associated with thyroid function. Conclusions: Bidirectional MR supports a directional association from hypothyroidism, increased TSH, and TPOAb, but not fT4, to decreased eGFR_crea and increased CKD.

A Mutation in γ-Adducin Impairs Autoregulation of Renal Blood Flow and Promotes the Development of Kidney Disease

BACKGROUND

The genes and mechanisms involved in the association between diabetes or hypertension and CKD risk are unclear. Previous studies have implicated a role for γ-adducin (ADD3), a cytoskeletal protein encoded by Add3.

METHODS

We investigated renal vascular function in vitro and in vivo and the susceptibility to CKD in rats with wild-type or mutated Add3 and in genetically modified rats with overexpression or knockout of ADD3. We also studied glomeruli and primary renal vascular smooth muscle cells isolated from these rats.

RESULTS

This study identified a K572Q mutation in ADD3 in fawn-hooded hypertensive (FHH) rats-a mutation previously reported in Milan normotensive (MNS) rats that also develop kidney disease. Using molecular dynamic simulations, we found that this mutation destabilizes a critical ADD3-ACTIN binding site. A reduction of ADD3 expression in membrane fractions prepared from the kidney and renal vascular smooth muscle cells of FHH rats was associated with the disruption of the F-actin cytoskeleton. Compared with renal vascular smooth muscle cells from Add3 transgenic rats, those from FHH rats had elevated membrane expression of BKα and BK channel current. FHH and Add3 knockout rats exhibited impairments in the myogenic response of afferent arterioles and in renal blood flow autoregulation, which were rescued in Add3 transgenic rats. We confirmed these findings in a genetic complementation study that involved crossing FHH and MNS rats that share the ADD3 mutation. Add3 transgenic rats showed attenuation of proteinuria, glomerular injury, and kidney fibrosis with aging and mineralocorticoid-induced hypertension.

CONCLUSIONS

This is the first report that a mutation in ADD3 that alters ACTIN binding causes renal vascular dysfunction and promotes the susceptibility to kidney disease.

The Genetic Architecture of Kidney Disease

The kidney is subject to a wide range of abnormalities, many of which have a significant hereditable component. Next generation sequencing is increasingly bringing the genetic drivers of Mendelian disease into focus at the base pair level, whereas inexpensive genotyping arrays have surveyed hundreds of thousands of individuals to identify common variants that predispose to kidney dysfunction. In this first article in a CJASN series on kidney genomics, we review how both rare and common variants contribute to kidney disease, explore how evolution may influence the genetic variants that affect kidney function, consider how genetic information is and will be used in the clinic, and identify some of the most important future directions for kidney disease research. Forthcoming articles in the series will elaborate on many of these themes.

Human C-terminal CUBN variants associate with chronic proteinuria and normal renal function

BACKGROUNDProteinuria is considered an unfavorable clinical condition that accelerates renal and cardiovascular disease. However, it is not clear whether all forms of proteinuria are damaging. Mutations in CUBN cause Imerslund-Gräsbeck syndrome (IGS), which is characterized by intestinal malabsorption of vitamin B12 and in some cases proteinuria. CUBN encodes for cubilin, an intestinal and proximal tubular uptake receptor containing 27 CUB domains for ligand binding.METHODSWe used next-generation sequencing for renal disease genes to genotype cohorts of patients with suspected hereditary renal disease and chronic proteinuria. CUBN variants were analyzed using bioinformatics, structural modeling, and epidemiological methods.RESULTSWe identified 39 patients, in whom biallelic pathogenic variants in the CUBN gene were associated with chronic isolated proteinuria and early childhood onset. Since the proteinuria in these patients had a high proportion of albuminuria, glomerular diseases such as steroid-resistant nephrotic syndrome or Alport syndrome were often the primary clinical diagnosis, motivating renal biopsies and the use of proteinuria-lowering treatments. However, renal function was normal in all cases. By contrast, we did not found any biallelic CUBN variants in proteinuric patients with reduced renal function or focal segmental glomerulosclerosis. Unlike the more N-terminal IGS mutations, 37 of the 41 proteinuria-associated CUBN variants led to modifications or truncations after the vitamin B12-binding domain. Finally, we show that 4 C-terminal CUBN variants are associated with albuminuria and slightly increased GFR in meta-analyses of large population-based cohorts.CONCLUSIONCollectively, our data suggest an important role for the C-terminal half of cubilin in renal albumin reabsorption. Albuminuria due to reduced cubilin function could be an unexpectedly common benign condition in humans that may not require any proteinuria-lowering treatment or renal biopsy.FUNDINGATIP-Avenir program, Fondation Bettencourt-Schueller (Liliane Bettencourt Chair of Developmental Biology), Agence Nationale de la Recherche (ANR) Investissements d'avenir program (ANR-10-IAHU-01) and NEPHROFLY (ANR-14-ACHN-0013, to MS), Steno Collaborative Grant 2018 (NNF18OC0052457, to TSA and MS), Heisenberg Professorship of the German Research Foundation (KO 3598/5-1, to AK), Deutsche Forschungsgemeinschaft (DFG) Collaborative Research Centre (SFB) KIDGEM 1140 (project 246781735, to CB), and Federal Ministry of Education and Research (BMB) (01GM1515C, to CB).

Human C-terminal CUBN variants associate with chronic proteinuria and normal renal function

BACKGROUNDProteinuria is considered an unfavorable clinical condition that accelerates renal and cardiovascular disease. However, it is not clear whether all forms of proteinuria are damaging. Mutations in CUBN cause Imerslund-Gräsbeck syndrome (IGS), which is characterized by intestinal malabsorption of vitamin B12 and in some cases proteinuria. CUBN encodes for cubilin, an intestinal and proximal tubular uptake receptor containing 27 CUB domains for ligand binding.METHODSWe used next-generation sequencing for renal disease genes to genotype cohorts of patients with suspected hereditary renal disease and chronic proteinuria. CUBN variants were analyzed using bioinformatics, structural modeling, and epidemiological methods.RESULTSWe identified 39 patients, in whom biallelic pathogenic variants in the CUBN gene were associated with chronic isolated proteinuria and early childhood onset. Since the proteinuria in these patients had a high proportion of albuminuria, glomerular diseases such as steroid-resistant nephrotic syndrome or Alport syndrome were often the primary clinical diagnosis, motivating renal biopsies and the use of proteinuria-lowering treatments. However, renal function was normal in all cases. By contrast, we did not found any biallelic CUBN variants in proteinuric patients with reduced renal function or focal segmental glomerulosclerosis. Unlike the more N-terminal IGS mutations, 37 of the 41 proteinuria-associated CUBN variants led to modifications or truncations after the vitamin B12-binding domain. Finally, we show that 4 C-terminal CUBN variants are associated with albuminuria and slightly increased GFR in meta-analyses of large population-based cohorts.CONCLUSIONCollectively, our data suggest an important role for the C-terminal half of cubilin in renal albumin reabsorption. Albuminuria due to reduced cubilin function could be an unexpectedly common benign condition in humans that may not require any proteinuria-lowering treatment or renal biopsy.FUNDINGATIP-Avenir program, Fondation Bettencourt-Schueller (Liliane Bettencourt Chair of Developmental Biology), Agence Nationale de la Recherche (ANR) Investissements d'avenir program (ANR-10-IAHU-01) and NEPHROFLY (ANR-14-ACHN-0013, to MS), Steno Collaborative Grant 2018 (NNF18OC0052457, to TSA and MS), Heisenberg Professorship of the German Research Foundation (KO 3598/5-1, to AK), Deutsche Forschungsgemeinschaft (DFG) Collaborative Research Centre (SFB) KIDGEM 1140 (project 246781735, to CB), and Federal Ministry of Education and Research (BMB) (01GM1515C, to CB).

The Urinary Excretion of Uromodulin is Regulated by the Potassium Channel ROMK

Uromodulin, the most abundant protein in normal urine, is produced by cells lining the thick ascending limb (TAL) of the loop of Henle. Uromodulin regulates the activity of the potassium channel ROMK in TAL cells. Common variants in KCNJ1, the gene encoding ROMK, are associated with urinary levels of uromodulin in population studies. Here, we investigated the functional link between ROMK and uromodulin in Kcnj1 knock-out mouse models, in primary cultures of mouse TAL (mTAL) cells, and in patients with Bartter syndrome due to KCNJ1 mutations. Both global and kidney-specific Kcnj1 knock-out mice showed reduced urinary levels of uromodulin paralleled by increased levels in the kidney, compared to wild-type controls. Pharmacological inhibition and genetic deletion of ROMK in mTAL cells caused a reduction in apical uromodulin excretion, reflected by cellular accumulation. In contrast, NKCC2 inhibition showed no effect on uromodulin processing. Patients with Bartter syndrome type 2 showed reduced urinary uromodulin levels compared to age and gender matched controls. These results demonstrate that ROMK directly regulates processing and release of uromodulin by TAL cells, independently from NKCC2. They support the functional link between transport activity and uromodulin in the TAL, relevant for blood pressure control and urinary concentrating ability.