SOS2 and ACP1 Loci Identified through Large-Scale Exome Chip Analysis Regulate Kidney Development and Function

Multivariate canonical correlation analysis identifies additional genetic variants for chronic kidney disease

Chronic kidney diseases (CKD) have genetic associations with kidney function. Univariate genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) associated with estimated glomerular filtration rate (eGFR) and blood urea nitrogen (BUN), two complementary kidney function markers. However, it is unknown whether additional SNPs for kidney function can be identified by multivariate statistical analysis. To address this, we applied canonical correlation analysis (CCA), a multivariate method, to two individual-level CKD genotype datasets, and metaCCA to two published GWAS summary statistics datasets. We identified SNPs previously associated with kidney function by published univariate GWASs with high replication rates, validating the metaCCA method. We then extended discovery and identified previously unreported lead SNPs for both kidney function markers, jointly. These showed expression quantitative trait loci (eQTL) colocalisation with genes having significant differential expression between CKD and healthy individuals. Several of these identified lead missense SNPs were predicted to have a functional impact, including in SLC14A2. We also identified previously unreported lead SNPs that showed significant correlation with both kidney function markers, jointly, in the European ancestry CKDGen, National Unified Renal Translational Research Enterprise (NURTuRE)-CKD and Salford Kidney Study (SKS) datasets. Of these, rs3094060 colocalised with FLOT1 gene expression and was significantly more common in CKD cases in both NURTURE-CKD and SKS, than in the general population. Overall, by using multivariate analysis by CCA, we identified additional SNPs and genes for both kidney function and CKD, that can be prioritised for further CKD analyses.

Genome-wide Association Analysis of Schizophrenia and Vitamin D Levels Shows Shared Genetic Architecture and Identifies Novel Risk Loci

Low vitamin D (vitD) levels have been consistently reported in schizophrenia (SCZ) suggesting a role in the etiopathology. However, little is known about the role of underlying shared genetic mechanisms. We applied a conditional/conjunctional false discovery rate approach (FDR) on large, nonoverlapping genome-wide association studies for SCZ (N cases = 53 386, N controls = 77 258) and vitD serum concentration (N = 417 580) to evaluate shared common genetic variants. The identified genomic loci were characterized using functional analyses and biological repositories. We observed cross-trait SNP enrichment in SCZ conditioned on vitD and vice versa, demonstrating shared genetic architecture. Applying the conjunctional FDR approach, we identified 72 loci jointly associated with SCZ and vitD at conjunctional FDR < 0.05. Among the 72 shared loci, 40 loci have not previously been reported for vitD, and 9 were novel for SCZ. Further, 64% had discordant effects on SCZ-risk and vitD levels. A mixture of shared variants with concordant and discordant effects with a predominance of discordant effects was in line with weak negative genetic correlation (rg = -0.085). Our results displayed shared genetic architecture between SCZ and vitD with mixed effect directions, suggesting overlapping biological pathways. Shared genetic variants with complex overlapping mechanisms may contribute to the coexistence of SCZ and vitD deficiency and influence the clinical picture.

Effectiveness of zebrafish models in understanding human diseases-A review of models

Understanding the detailed mechanism behind every human disease, disorder, defect, and deficiency is a daunting task concerning the clinical diagnostic tools for patients. Hence, a closely resembling living or simulated model is of paramount interest for the development and testing of a probable novel drug for rectifying the conditions pertaining to the various ailments. The animal model that can be easily genetically manipulated to suit the study of the therapeutic motive is an indispensable asset and within the last few decades, the zebrafish models have proven their effectiveness by becoming such potent human disease models with their use being extended to various avenues of research to understand the underlying mechanisms of the diseases. As zebrafish are explored as model animals in understanding the molecular basis and genetics of many diseases owing to the 70% genetic homology between the human and zebrafish genes; new and fascinating facts about the diseases are being surfaced, establishing it as a very powerful tool for upcoming research. These prospective research areas can be explored in the near future using zebrafish as a model. In this review, the effectiveness of the zebrafish as an animal model against several human diseases such as osteoporosis, atrial fibrillation, Noonan syndrome, leukemia, autism spectrum disorders, etc. has been discussed.

Biomarkers of kidney function and cognitive ability: A Mendelian randomization study

BACKGROUND

Estimated glomerular filtration rate (eGFR), albuminuria and serum uric acid (SUA) are markers of kidney function that have been associated with cognitive ability. However, whether these associations are causal is unclear.

METHODS

We performed one-sample Mendelian randomization (MR) to estimate the effects of kidney function markers on cognitive performance using data from the UK Biobank. Polygenic scores for SUA, urine albumin to creatinine ratio (ACR), estimated glomerular filtration rate based on serum creatinine (eGFRcre) and serum cystatin C (eGFRcys) were used as instrumental variables, and cognitive function outcomes included tests of verbal-numeric reasoning, reaction time, visual memory, and numeric memory.

RESULTS

We found no evidence of a causal effect of genetically determined SUA, eGFRcre or eGFRcys on cognitive function outcomes. There was no association between a polygenic score for ACR and verbal-numeric reasoning or numeric memory. However, there was suggestive evidence of a relationship between genetically increased ACR and slower reaction time and worse visual memory. ACR was no longer significantly associated with visual memory in analyses using an unweighted polygenic score and in analyses stratified by sex and age category. Pleiotropy adjusted estimates were directionally consistent with those of the principal analysis but overlapped with the null.

CONCLUSIONS

This MR study does not support causal effects of SUA, eGFRcre or eGFRcys on cognitive performance. Genetically increased ACR was associated with slower processing speed and visual memory, but results need confirmation in independent samples.

First prenatal case of Noonan syndrome with SOS2 mutation: Implications of early diagnosis for genetic counseling

RASopathies are a group of syndromes with partially overlapping clinical features caused by germline mutations of the RAS/MAPK signaling pathway genes. The most common disorder is Noonan syndrome (NS; MIM 163950). We report the first prenatal case of NS with SOS2 (NM_006939.4) mutation in a euploid fetus with a severe increase in nuchal translucency (NT > 12 mm). Trio-based custom next-generation sequencing detected a de novo heterozygous missense mutation in the SOS2 gene: c.800 T > A (p.Met267Lys). Owing to the marked variable expressivity of NS and the scarcity of SOS2 mutation-related NS cases reported in the literature, it is difficult to provide appropriate genetic counseling. Several issues such as the best management technique and optimal NT cutoff have been discussed. In addition, in general, the fine balance between the advantages of an early prenatal diagnosis and the challenge of determining if the detected gene variant is pathogenic and, primarily, the stress of the counselees when providing a genetic counseling with limited information on the prenatal phenotype have been discussed. A prenatal path comprising examinations and multidisciplinary counseling is essential to support couples in a shared decision-making process.

Whole genome sequence analyses of eGFR in 23,732 people representing multiple ancestries in the NHLBI trans-omics for precision medicine (TOPMed) consortium

Background

Genetic factors that influence kidney traits have been understudied for low frequency and ancestry-specific variants.

Methods

We combined whole genome sequencing (WGS) data from 23,732 participants from 10 NHLBI Trans-Omics for Precision Medicine (TOPMed) Program multi-ethnic studies to identify novel loci for estimated glomerular filtration rate (eGFR). Participants included European, African, East Asian, and Hispanic ancestries. We applied linear mixed models using a genetic relationship matrix estimated from the WGS data and adjusted for age, sex, study, and ethnicity.

Findings

When testing single variants, we identified three novel loci driven by low frequency variants more commonly observed in non-European ancestry (PRKAA2, rs180996919, minor allele frequency [MAF] 0.04%, P = 6.1 × 10⁻¹¹; METTL8, rs116951054, MAF 0.09%, P = 4.5 × 10⁻⁹; and MATK, rs539182790, MAF 0.05%, P = 3.4 × 10⁻⁹). We also replicated two known loci for common variants (rs2461702, MAF=0.49, P = 1.2 × 10⁻⁹, nearest gene GATM, and rs71147340, MAF=0.34, P = 3.3 × 10⁻⁹, CDK12). Testing aggregated variants within a gene identified the MAF gene. A statistical approach based on local ancestry helped to identify replication samples for ancestry-specific variants.

Interpretation

This study highlights challenges in studying variants influencing kidney traits that are low frequency in populations and more common in non-European ancestry.

AHR canonical pathway: in vivo findings to support novel antihypertensive strategies

Essential hypertension (HTN) is a disease where genetic and environmental factors interact to produce a high prevalent set of almost indistinguishable phenotypes. The weak definition of what is under the umbrella of HTN is a consequence of the lack of knowledge on the players involved in environment-gene interaction and their impact on blood pressure (BP) and mechanisms. The disclosure of these mechanisms that sense and (mal)adapt to toxic-environmental stimuli might at least determine some phenotypes of essential HTN and will have important therapeutic implications. In the present manuscript, we looked closer to the environmental sensor aryl hydrocarbon receptor (AHR), a ligand-activated transcription factor involved in cardiovascular physiology, but better known by its involvement in biotransformation of xenobiotics through its canonical pathway. This review aims to disclose the contribution of the AHR-canonical pathway to HTN. For better mirror the complexity of the mechanisms involved in BP regulation, we privileged evidence from in vivo studies. Here we ascertained the level of available evidence and a comprehensive characterization of the AHR-related phenotype of HTN. We reviewed clinical and rodent studies on AHR-HTN genetic association and on AHR ligands and their impact on BP. We concluded that AHR is a druggable mechanistic linker of environmental exposure to HTN. We conclude that is worth to investigate the canonical pathway of AHR and the expression/polymorphisms of its related genes and/or other biomarkers (e.g. tryptophan-related ligands), in order to identify patients that may benefit from an AHR-centered antihypertensive treatment.

Discovery of rare variants associated with blood pressure regulation through meta-analysis of 1.3 million individuals

Genetic studies of blood pressure (BP) to date have mainly analyzed common variants (minor allele frequency > 0.05). In a meta-analysis of up to ~1.3 million participants, we discovered 106 new BP-associated genomic regions and 87 rare (minor allele frequency ≤ 0.01) variant BP associations (P < 5 × 10-8), of which 32 were in new BP-associated loci and 55 were independent BP-associated single-nucleotide variants within known BP-associated regions. Average effects of rare variants (44% coding) were ~8 times larger than common variant effects and indicate potential candidate causal genes at new and known loci (for example, GATA5 and PLCB3). BP-associated variants (including rare and common) were enriched in regions of active chromatin in fetal tissues, potentially linking fetal development with BP regulation in later life. Multivariable Mendelian randomization suggested possible inverse effects of elevated systolic and diastolic BP on large artery stroke. Our study demonstrates the utility of rare-variant analyses for identifying candidate genes and the results highlight potential therapeutic targets.

Genome-Wide Association Studies of CKD and Related Traits

The past few years have seen major advances in genome-wide association studies (GWAS) of CKD and kidney function-related traits in several areas: increases in sample size from >100,000 to >1 million, enabling the discovery of >250 associated genetic loci that are highly reproducible; the inclusion of participants not only of European but also of non-European ancestries; and the use of advanced computational methods to integrate additional genomic and other unbiased, high-dimensional data to characterize the underlying genetic architecture and prioritize potentially causal genes and variants. Together with other large-scale biobank and genetic association studies of complex traits, these GWAS of kidney function-related traits have also provided novel insight into the relationship of kidney function to other diseases with respect to their genetic associations, genetic correlation, and directional relationships. A number of studies also included functional experiments using model organisms or cell lines to validate prioritized potentially causal genes and/or variants. In this review article, we will summarize these recent GWAS of CKD and kidney function-related traits, explain approaches for downstream characterization of associated genetic loci and the value of such computational follow-up analyses, and discuss related challenges along with potential solutions to ultimately enable improved treatment and prevention of kidney diseases through genetics.

SOS GEFs in health and disease

SOS1 and SOS2 are the most universal and widely expressed family of guanine exchange factors (GEFs) capable or activating RAS or RAC1 proteins in metazoan cells. SOS proteins contain a sequence of modular domains that are responsible for different intramolecular and intermolecular interactions modulating mechanisms of self-inhibition, allosteric activation and intracellular homeostasis. Despite their homology, analyses of SOS1/2-KO mice demonstrate functional prevalence of SOS1 over SOS2 in cellular processes including proliferation, migration, inflammation or maintenance of intracellular redox homeostasis, although some functional redundancy cannot be excluded, particularly at the organismal level. Specific SOS1 gain-of-function mutations have been identified in inherited RASopathies and various sporadic human cancers. SOS1 depletion reduces tumorigenesis mediated by RAS or RAC1 in mouse models and is associated with increased intracellular oxidative stress and mitochondrial dysfunction. Since WT RAS is essential for development of RAS-mutant tumors, the SOS GEFs may be considered as relevant biomarkers or therapy targets in RAS-dependent cancers. Inhibitors blocking SOS expression, intrinsic GEF activity, or productive SOS protein-protein interactions with cellular regulators and/or RAS/RAC targets have been recently developed and shown preclinical and clinical effectiveness blocking aberrant RAS signaling in RAS-driven and RTK-driven tumors.

Combined Utility of 25 Disease and Risk Factor Polygenic Risk Scores for Stratifying Risk of All-Cause Mortality

While genome-wide association studies have identified susceptibility variants for numerous traits, their combined utility for predicting broad measures of health, such as mortality, remains poorly understood. We used data from the UK Biobank to combine polygenic risk scores (PRS) for 13 diseases and 12 mortality risk factors into sex-specific composite PRS (cPRS). These cPRS were moderately associated with all-cause mortality in independent data within the UK Biobank: the estimated hazard ratios per standard deviation were 1.10 (95% confidence interval: 1.05, 1.16) and 1.15 (1.10, 1.19) for women and men, respectively. Differences in life expectancy between the top and bottom 5% of the cPRS were estimated to be 4.79 (1.76, 7.81) years and 6.75 (4.16, 9.35) years for women and men, respectively. These associations were substantially attenuated after adjusting for non-genetic mortality risk factors measured at study entry (i.e., middle age for most participants). The cPRS may be useful in counseling younger individuals at higher genetic risk of mortality on modification of non-genetic factors.

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.

Genetic overlap and causal inferences between kidney function and cerebrovascular disease

OBJECTIVE

Leveraging large-scale genetic data, we aimed to identify shared pathogenic mechanisms and causal relationships between impaired kidney function and cerebrovascular disease phenotypes.

METHODS

We used summary statistics from genome-wide association studies (GWAS) of kidney function traits (chronic kidney disease diagnosis, estimated glomerular filtration rate [eGFR], and urinary albumin-to-creatinine ratio [UACR]) and cerebrovascular disease phenotypes (ischemic stroke and its subtypes, intracerebral hemorrhage [ICH], and white matter hyperintensities [WMH] on brain MRI). We (1) tested the genetic overlap between them with polygenic risk scores (PRS), (2) searched for common pleiotropic loci with pairwise GWAS analyses, and (3) explored causal associations by employing 2-sample Mendelian randomization.

RESULTS

A PRS for lower eGFR was associated with higher large artery stroke (LAS) risk (p = 1 × 10-4). Multiple pleiotropic loci were identified between kidney function traits and cerebrovascular disease phenotypes, with 12q24 associated with eGFR and both LAS and small vessel stroke (SVS), and 2q33 associated with UACR and both SVS and WMH. Mendelian randomization revealed associations of both lower eGFR (odds ratio [OR] per 1-log decrement, 2.10; 95% confidence interval [CI], 1.38-3.21) and higher UACR (OR per 1-log increment, 2.35; 95% CI, 1.12-4.94) with a higher risk of LAS, as well as between higher UACR and higher risk of ICH.

CONCLUSIONS

Impaired kidney function, as assessed by decreased eGFR and increased UACR, may be causally involved in the pathogenesis of LAS. Increased UACR, previously proposed as a marker of systemic small vessel disease, is involved in ICH risk and shares a genetic risk factor at 2q33 with manifestations of cerebral small vessel disease.

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.

Family-Based Quantitative Trait Meta-Analysis Implicates Rare Noncoding Variants in DENND1A in Polycystic Ovary Syndrome

CONTEXT

Polycystic ovary syndrome (PCOS) is among the most common endocrine disorders of premenopausal women, affecting 5% to15% of this population depending on the diagnostic criteria applied. It is characterized by hyperandrogenism, ovulatory dysfunction, and polycystic ovarian morphology. PCOS is highly heritable, but only a small proportion of this heritability can be accounted for by the common genetic susceptibility variants identified to date.

OBJECTIVE

The objective of this study was to test whether rare genetic variants contribute to PCOS pathogenesis.

DESIGN, PATIENTS, AND METHODS

We performed whole-genome sequencing on DNA from 261 individuals from 62 families with one or more daughters with PCOS. We tested for associations of rare variants with PCOS and its concomitant hormonal traits using a quantitative trait meta-analysis.

RESULTS

We found rare variants in DENND1A (P = 5.31 × 10-5, adjusted P = 0.039) that were significantly associated with reproductive and metabolic traits in PCOS families.

CONCLUSIONS

Common variants in DENND1A have previously been associated with PCOS diagnosis in genome-wide association studies. Subsequent studies indicated that DENND1A is an important regulator of human ovarian androgen biosynthesis. Our findings provide additional evidence that DENND1A plays a central role in PCOS and suggest that rare noncoding variants contribute to disease pathogenesis.

Mapping eGFR loci to the renal transcriptome and phenome in the VA Million Veteran Program

Chronic kidney disease (CKD), defined by low estimated glomerular filtration rate (eGFR), contributes to global morbidity and mortality. Here we conduct a transethnic Genome-Wide Association Study of eGFR in 280,722 participants of the Million Veteran Program (MVP), with replication in 765,289 participants from the Chronic Kidney Disease Genetics (CKDGen) Consortium. We identify 82 previously unreported variants, confirm 54 loci, and report interesting findings including association of the sickle cell allele of betaglobin among non-Hispanic blacks. Our transcriptome-wide association study of kidney function in healthy kidney tissue identifies 36 previously unreported and nine known genes, and maps gene expression to renal cell types. In a Phenome-Wide Association Study in 192,868 MVP participants using a weighted genetic score we detect associations with CKD stages and complications and kidney stones. This investigation reinterprets the genetic architecture of kidney function to identify the gene, tissue, and anatomical context of renal homeostasis and the clinical consequences of dysregulation.

The impact of donor and recipient common clinical and genetic variation on estimated glomerular filtration rate in a European renal transplant population

Genetic variation across the human leukocyte antigen loci is known to influence renal-transplant outcome. However, the impact of genetic variation beyond the human leukocyte antigen loci is less clear. We tested the association of common genetic variation and clinical characteristics, from both the donor and recipient, with posttransplant eGFR at different time-points, out to 5 years posttransplantation. We conducted GWAS meta-analyses across 10 844 donors and recipients from five European ancestry cohorts. We also analyzed the impact of polygenic risk scores (PRS), calculated using genetic variants associated with nontransplant eGFR, on posttransplant eGFR. PRS calculated using the recipient genotype alone, as well as combined donor and recipient genotypes were significantly associated with eGFR at 1-year posttransplant. Thirty-two percent of the variability in eGFR at 1-year posttransplant was explained by our model containing clinical covariates (including weights for death/graft-failure), principal components and combined donor-recipient PRS, with 0.3% contributed by the PRS. No individual genetic variant was significantly associated with eGFR posttransplant in the GWAS. This is the first study to examine PRS, composed of variants that impact kidney function in the general population, in a posttransplant context. Despite PRS being a significant predictor of eGFR posttransplant, the effect size of common genetic factors is limited compared to clinical variables.

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

Chronic kidney disease (CKD) is responsible for a public health burden with multi-systemic complications. Through trans-ancestry meta-analysis of genome-wide association studies of estimated glomerular filtration rate (eGFR) and independent replication (n = 1,046,070), we identified 264 associated loci (166 new). Of these, 147 were likely to be relevant for kidney function on the basis of associations with the alternative kidney function marker blood urea nitrogen (n = 416,178). Pathway and enrichment analyses, including mouse models with renal phenotypes, support the kidney as the main target organ. A genetic risk score for lower eGFR was associated with clinically diagnosed CKD in 452,264 independent individuals. Colocalization analyses of associations with eGFR among 783,978 European-ancestry individuals and gene expression across 46 human tissues, including tubulo-interstitial and glomerular kidney compartments, identified 17 genes differentially expressed in kidney. Fine-mapping highlighted missense driver variants in 11 genes and kidney-specific regulatory variants. These results provide a comprehensive priority list of molecular targets for translational research.

Zebrafish as a model for kidney function and disease

Kidney disease is a global problem with around three million people diagnosed in the UK alone and the incidence is rising. Research is critical to develop better treatments. Animal models can help to better understand the pathophysiology behind the various kidney diseases and to screen for therapeutic compounds, but the use especially of mammalian models should be minimised in the interest of animal welfare. Zebrafish are increasingly used, as they are genetically tractable and have a basic renal anatomy comparable to mammalian kidneys with glomerular filtration and tubular filtration processing. Here, we discuss how zebrafish have advanced the study of nephrology and the mechanisms underlying kidney disease.

Sex-specific and pleiotropic effects underlying kidney function identified from GWAS meta-analysis

Chronic kidney disease (CKD) is a growing health burden currently affecting 10–15% of adults worldwide. Estimated glomerular filtration rate (eGFR) as a marker of kidney function is commonly used to diagnose CKD. We analyze eGFR data from the Nord-Trøndelag Health Study and Michigan Genomics Initiative and perform a GWAS meta-analysis with public summary statistics, more than doubling the sample size of previous meta-analyses. We identify 147 loci (53 novel) associated with eGFR, including genes involved in transcriptional regulation, kidney development, cellular signaling, metabolism, and solute transport. Additionally, sex-stratified analysis identifies one locus with more significant effects in women than men. Using genetic risk scores constructed from these eGFR meta-analysis results, we show that associated variants are generally predictive of CKD with only modest improvements in detection compared with other known clinical risk factors. Collectively, these results yield additional insight into the genetic factors underlying kidney function and progression to CKD.

Estimated glomerular filtration rate (eGFR) is a measure of kidney function and used to characterize chronic kidney disease. Here, Graham et al. identify 53 novel loci for eGFR in a GWAS meta-analysis, a subset of which are associated with other common diseases, such as diabetes and hypertension, based on PheWAS.

A novel rare CUBN variant and three additional genes identified in Europeans with and without diabetes: results from an exome-wide association study of albuminuria

AIMS/HYPOTHESIS

Identifying rare coding variants associated with albuminuria may open new avenues for preventing chronic kidney disease and end-stage renal disease, which are highly prevalent in individuals with diabetes. Efforts to identify genetic susceptibility variants for albuminuria have so far been limited, with the majority of studies focusing on common variants.

METHODS

We performed an exome-wide association study to identify coding variants in a two-stage (discovery and replication) approach. Data from 33,985 individuals of European ancestry (15,872 with and 18,113 without diabetes) and 2605 Greenlanders were included.

RESULTS

We identified a rare (minor allele frequency [MAF]: 0.8%) missense (A1690V) variant in CUBN (rs141640975, β = 0.27, p = 1.3 × 10^-11) associated with albuminuria as a continuous measure in the combined European meta-analysis. The presence of each rare allele of the variant was associated with a 6.4% increase in albuminuria. The rare CUBN variant had an effect that was three times stronger in individuals with type 2 diabetes compared with those without (p_interaction = 7.0 × 10^-4, β with diabetes = 0.69, β without diabetes = 0.20) in the discovery meta-analysis. Gene-aggregate tests based on rare and common variants identified three additional genes associated with albuminuria (HES1, CDC73 and GRM5) after multiple testing correction (p_Bonferroni < 2.7 × 10^-6).

CONCLUSIONS/INTERPRETATION

The current study identifies a rare coding variant in the CUBN locus and other potential genes associated with albuminuria in individuals with and without diabetes. These genes have been implicated in renal and cardiovascular dysfunction. The findings provide new insights into the genetic architecture of albuminuria and highlight target genes and pathways for the prevention of diabetes-related kidney disease.

Genome-Wide Association Studies of Metabolite Concentrations (mGWAS): Relevance for Nephrology

Metabolites are small molecules that are intermediates or products of metabolism, many of which are freely filtered by the kidneys. In addition, the kidneys have a central role in metabolite anabolism and catabolism, as well as in active metabolite reabsorption and/or secretion during tubular passage. This review article illustrates how the coupling of genomics and metabolomics in genome-wide association analyses of metabolites can be used to illuminate mechanisms underlying human metabolism, with a special focus on insights relevant to nephrology. First, genetic susceptibility loci for reduced kidney function and chronic kidney disease (CKD) were reviewed systematically for their associations with metabolite concentrations in metabolomics studies of blood and urine. Second, kidney function and CKD-associated metabolites reported from observational studies were interrogated for metabolite-associated genetic variants to generate and discuss complementary insights. Finally, insights originating from the simultaneous study of both blood and urine or by modeling intermetabolite relationships are summarized. We also discuss methodologic questions related to the study of metabolite concentrations in urine as well as among CKD patients. In summary, genome-wide association analyses of metabolites using metabolite concentrations quantified from blood and/or urine are a promising avenue of research to illuminate physiological and pathophysiological functions of the kidney.

Comparative Analysis and Refinement of Human PSC-Derived Kidney Organoid Differentiation with Single-Cell Transcriptomics

Kidney organoids derived from human pluripotent stem cells have great utility for investigating organogenesis and disease mechanisms and, potentially, as a replacement tissue source, but how closely organoids derived from current protocols replicate adult human kidney is undefined. We compared two directed differentiation protocols by single-cell transcriptomics of 83,130 cells from 65 organoids with single-cell transcriptomes of fetal and adult kidney cells. Both protocols generate a diverse range of kidney cells with differing ratios, but organoid-derived cell types are immature, and 10%-20% of cells are non-renal. Reconstructing lineage relationships by pseudotemporal ordering identified ligands, receptors, and transcription factor networks associated with fate decisions. Brain-derived neurotrophic factor (BDNF) and its cognate receptor NTRK2 were expressed in the neuronal lineage during organoid differentiation. Inhibiting this pathway improved organoid formation by reducing neurons by 90% without affecting kidney differentiation, highlighting the power of single-cell technologies to characterize and improve organoid differentiation.

Large-scale whole-exome sequencing association studies identify rare functional variants influencing serum urate levels

Elevated serum urate levels can cause gout, an excruciating disease with suboptimal treatment. Previous GWAS identified common variants with modest effects on serum urate. Here we report large-scale whole-exome sequencing association studies of serum urate and kidney function among ≤19,517 European ancestry and African-American individuals. We identify aggregate associations of low-frequency damaging variants in the urate transporters SLC22A12 (URAT1; p = 1.3 × 10⁻⁵⁶) and SLC2A9 (p = 4.5 × 10⁻⁷). Gout risk in rare SLC22A12 variant carriers is halved (OR = 0.5, p = 4.9 × 10⁻³). Selected rare variants in SLC22A12 are validated in transport studies, confirming three as loss-of-function (R325W, R405C, and T467M) and illustrating the therapeutic potential of the new URAT1-blocker lesinurad. In SLC2A9, mapping of rare variants of large effects onto the predicted protein structure reveals new residues that may affect urate binding. These findings provide new insights into the genetic architecture of serum urate, and highlight molecular targets in SLC22A12 and SLC2A9 for lowering serum urate and preventing gout.

Elevated serum urate levels are a risk factor for gout. Here, Tin et al. perform whole-exome sequencing in 19,517 individuals and detect low-frequency genetic variants in urate transporter genes, SLC22A12 and SLC2A9, associated with serum urate levels and confirm their damaging nature in vitro and in silico.

Lessons from CKD-Related Genetic Association Studies-Moving Forward

Over the past decade, genetic association studies have uncovered numerous determinants of kidney function in the general, diabetic, hypertensive, CKD, ESRD, and GN-based study populations (e.g., IgA nephropathy, membranous nephropathy, FSGS). These studies have led to numerous novel and unanticipated findings, which are helping improve our understanding of factors and pathways affecting both normal and pathologic kidney function. In this review, we report on major discoveries and advances resulting from this rapidly progressing research domain. We also predict some of the next steps the nephrology community should embrace to accelerate the identification of genetic and molecular processes leading to kidney dysfunction, pathophysiologically based disease subgroups, and specific therapeutic targets, as we attempt to transition toward a more precision-based medicine approach.

Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

Genome-wide association studies (GWAS) have identified >250 loci for body mass index (BMI), implicating pathways related to neuronal biology. Most GWAS loci represent clusters of common, noncoding variants from which pinpointing causal genes remains challenging. Here we combined data from 718,734 individuals to discover rare and low-frequency (minor allele frequency (MAF) < 5%) coding variants associated with BMI. We identified 14 coding variants in 13 genes, of which 8 variants were in genes (ZBTB7B, ACHE, RAPGEF3, RAB21, ZFHX3, ENTPD6, ZFR2 and ZNF169) newly implicated in human obesity, 2 variants were in genes (MC4R and KSR2) previously observed to be mutated in extreme obesity and 2 variants were in GIPR. The effect sizes of rare variants are ~10 times larger than those of common variants, with the largest effect observed in carriers of an MC4R mutation introducing a stop codon (p.Tyr35Ter, MAF = 0.01%), who weighed ~7 kg more than non-carriers. Pathway analyses based on the variants associated with BMI confirm enrichment of neuronal genes and provide new evidence for adipocyte and energy expenditure biology, widening the potential of genetically supported therapeutic targets in obesity.

The UMOD Locus: Insights into the Pathogenesis and Prognosis of Kidney Disease

The identification of genetic factors associated with kidney disease has the potential to provide critical insights into disease mechanisms. Genome-wide association studies have uncovered genomic regions associated with renal function metrics and risk of CKD. UMOD is among the most outstanding loci associated with CKD in the general population, because it has a large effect on eGFR and CKD risk that is consistent across different ethnic groups. The relevance of UMOD for CKD is clear, because the encoded protein, uromodulin (Tamm-Horsfall protein), is exclusively produced by the kidney tubule and has specific biochemical properties that mediate important functions in the kidney and urine. Rare mutations in UMOD are the major cause of autosomal dominant tubulointerstitial kidney disease, a condition that leads to CKD and ESRD. In this brief review, we use the UMOD paradigm to describe how population genetic studies can yield insight into the pathogenesis and prognosis of kidney diseases.

Identification of CDC42BPG as a novel susceptibility locus for hyperuricemia in a Japanese population

Chronic kidney disease and hyperuricemia are serious global health problems. Recent genome-wide association studies have identified various genetic variants related to these disorders. However, most studies have been conducted in a cross-sectional manner. To identify novel susceptibility loci for chronic kidney disease or hyperuricemia, we performed longitudinal exome-wide association studies (EWASs), using ~ 244,000 genetic variants and clinical data of Japanese individuals who had undergone annual health checkups for several years. After establishing quality controls, the association of renal function-related traits in 5648 subjects (excluding patients with dialysis and population outliers) with 24,579 single nucleotide variants (SNVs) for three genetic models (P < 3.39 × 10^− 7) was tested using generalized estimating equation models. The longitudinal EWASs revealed novel relations of five SNVs to renal function-related traits. Cross-sectional data for renal function-related traits in 7699 Japanese subjects were examined in a replication study. Among the five SNVs, rs55975541 in CDC42BPG was significantly (P < 4.90 × 10^− 4) related to the serum concentration of uric acid in the replication cohort. We also examined the SNVs detected in our longitudinal EWASs with the information on P values in GKDGEN meta-analysis data. Four SNVs in SLC15A2 were significantly associated with the estimated glomerular filtration rate in European ancestry populations, although these SNVs were related to the serum concentration of uric acid with borderline significance in our longitudinal EWASs. Our findings indicate that CDC42BPG may be a novel susceptibility locus for hyperuricemia.

The Loss of GSTM1 Associates with Kidney Failure and Heart Failure

Glutathione S-transferase mu 1 (GSTM1) encodes an enzyme that catalyzes the conjugation of electrophilic compounds with glutathione to facilitate their degradation or excretion. The loss of one or both copies of GSTM1 is common in many populations and has been associated with CKD progression. With the hypothesis that the loss of GSTM1 is also associated with incident kidney failure and heart failure, we estimated GSTM1 copy number using exome sequencing reads in the Atherosclerosis Risk in Communities (ARIC) Study, a community-based prospective cohort of white and black participants. Overall, 51.2% and 39.8% of white participants and 25.6% and 48.5% of black participants had zero or one copy of GSTM1, respectively. Over a median follow-up of 24.6 years, 256 kidney failure events occurred in 5715 participants without prevalent kidney failure, and 1028 heart failure events occurred in 5368 participants without prevalent heart failure. In analysis adjusted for demographics, diabetes, and hypertension, having zero or one copy of GSTM1 associated with higher risk of kidney failure and heart failure (adjusted hazard ratio [95% confidence interval] for zero or one versus two copies of GSTM1: kidney failure, 1.66 [1.27 to 2.17]; heart failure, 1.16 [1.04 to 1.29]). Risk did not differ significantly between participants with zero and one copy of GSTM1 (P>0.10). In summary, the loss of GSTM1 was significantly associated with incident kidney and heart failure, independent of traditional risk factors. These results suggest GSTM1 function is a potential treatment target for the prevention of kidney and heart failure.