Interactive exploration of adverse events and multimorbidity in CKD.

BACKGROUND AND HYPOTHESIS: Persons with chronic kidney disease (CKD) are at increased risk of adverse events, early mortality, and multimorbidity. A detailed overview of adverse event types and rates from a large CKD cohort under regular nephrological care is missing. We generated an interactive tool to enable exploration of adverse events and their combinations in the prospective, observational German CKD (GCKD) study. METHODS: The GCKD study enrolled 5217 participants under regular nephrological care with an estimated glomerular filtration rate of 30-60 or >60 mL/min/1.73m2 and an overt proteinuria. Cardio-, cerebro- and peripheral vascular, kidney, infection, and cancer events, as well as deaths were adjudicated following a standard operation procedure. We summarized these time-to-event data points for exploration in interactive graphs within an R shiny app. Multivariable adjusted Cox models for time to first event were fitted. Cumulative incidence functions, Kaplan-Meier curves and intersection plots were used to display main adverse events and their combinations by sex and CKD etiology. RESULTS: Over a median of 6.5 years, 10 271 events occurred in total and 680 participants (13.0%) died while 2947 participants (56.5%) experienced any event. The new publicly available interactive platform enables readers to scrutinize adverse events and their combinations as well as mortality trends as a gateway to better understand multimorbidity in CKD: incident rates per 1000 patient-years varied by event type, CKD etiology, and baseline characteristics. Incidence rates for the most frequent events and their recurrence were 113.6 (cardiovascular), 75.0 (kidney), and 66.0 (infection). Participants with diabetic kidney disease and men were more prone to experiencing events. CONCLUSION: This comprehensive explorative tool to visualize adverse events (https://gckd.diz.uk-erlangen.de/), their combination, mortality, and multimorbidity among persons with CKD may manifest as a valuable resource for patient care, identification of high-risk groups, health services, and public health policy planning.

SLC17 transporters mediate renal excretion of Lac-Phe in mice and humans.

N-lactoyl-phenylalanine (Lac-Phe) is a lactate-derived metabolite that suppresses food intake and body weight. Little is known about the mechanisms that mediate Lac-Phe transport across cell membranes. Here we identify SLC17A1 and SLC17A3, two kidney-restricted plasma membrane-localized solute carriers, as physiologic urine Lac-Phe transporters. In cell culture, SLC17A1/3 exhibit high Lac-Phe efflux activity. In humans, levels of Lac-Phe in urine exhibit a strong genetic association with the SLC17A1-4 locus. Urine Lac-Phe levels are also increased following a Wingate sprint test. In mice, genetic ablation of either SLC17A1 or SLC17A3 reduces urine Lac-Phe levels. Despite these differences, both knockout strains have normal blood Lac-Phe and body weights, demonstrating that urine and plasma Lac-Phe pools are functionally and biochemically de-coupled. Together, these data establish SLC17 family members as the physiologic urine transporters for Lac-Phe and uncover a biochemical pathway for the renal excretion of this signaling metabolite.

Integrated Proteomic and Metabolomic Modules associated with Risk of Kidney Disease Progression.

BACKGROUND: Proteins and metabolites play crucial roles in various biological functions and are frequently interconnected through enzymatic or transport processes. METHODS: We present an integrated analysis of 4,091 proteins and 630 metabolites in the Chronic Renal Insufficiency Cohort Study (N=1,708; average follow-up for kidney failure [KF], 9.5 years, with 537 events). Proteins and metabolites were integrated using an unsupervised clustering method and we assessed associations between clusters and CKD progression and kidney failure using Cox proportional hazards models. Analyses were adjusted for demographics and risk factors including the estimated glomerular filtration rate (eGFR) and urine protein-creatinine ratio. Associations were identified in a discovery sample (random two-thirds, N=1139) and then evaluated in a replication sample (one-third, N=569). RESULTS: We identified 139 modules of correlated proteins and metabolites, which were represented by their principal components (PC). Modules and PC loadings were projected onto the replication sample which demonstrated a consistent network structure. Two modules, representing a total of 236 proteins and 82 metabolites, were robustly associated with both CKD progression and kidney failure in both discovery and validation samples. Using gene set enrichment, several transmembrane related terms were identified as over-represented in these modules. Transmembrane-ephrin receptor activity displayed the largest odds (OR = 13.2, P-value = 5.5×10 -5 ). A module containing CRIM1 and NPNT expressed in podocytes demonstrated particularly strong associations with kidney failure (P-value = 2.6×10 -5 ). CONCLUSIONS: This study demonstrates that integration of the proteome and metabolome can identify functions of pathophysiologic importance in kidney disease.

Clonal hematopoiesis of indeterminate potential is associated with acute kidney injury.

Age is a predominant risk factor for acute kidney injury (AKI), yet the biological mechanisms underlying this risk are largely unknown. Clonal hematopoiesis of indeterminate potential (CHIP) confers increased risk for several chronic diseases associated with aging. Here we sought to test whether CHIP increases the risk of AKI. In three population-based epidemiology cohorts, we found that CHIP was associated with a greater risk of incident AKI, which was more pronounced in patients with AKI requiring dialysis and in individuals with somatic mutations in genes other than DNMT3A, including mutations in TET2 and JAK2. Mendelian randomization analyses supported a causal role for CHIP in promoting AKI. Non-DNMT3A-CHIP was also associated with a nonresolving pattern of injury in patients with AKI. To gain mechanistic insight, we evaluated the role of Tet2-CHIP and Jak2V617F-CHIP in two mouse models of AKI. In both models, CHIP was associated with more severe AKI, greater renal proinflammatory macrophage infiltration and greater post-AKI kidney fibrosis. In summary, this work establishes CHIP as a genetic mechanism conferring impaired kidney function recovery after AKI via an aberrant inflammatory response mediated by renal macrophages.

Imaging Markers Derived From MRI-Based Automated Kidney Segmentation­an Analysis of Data From the German National Cohort (NAKO Gesundheitsstudie)

BACKGROUND: Population-wide research on potential new imaging biomarkers of the kidney depends on accurate automated segmentation of the kidney and its compartments (cortex, medulla, and sinus). METHODS: We developed a robust deep-learning framework for kidney (sub-)segmentation based on a hierarchical, three-dimensional convolutional neural network (CNN) that was optimized for multi-scale problems of combined localization and segmentation. We applied the CNN to abdominal magnetic resonance images from the population-based German National Cohort (NAKO) study. RESULTS: There was good to excellent agreement between the model predictions and manual segmentations. The median values for the body-surface normalized total kidney, cortex, medulla, and sinus volumes of 9934 persons were 158, 115, 43, and 24 mL/m2. Distributions of these markers are provided both for the overall study population and for a subgroup of persons without kidney disease or any associated conditions. Multivariable adjusted regression analyses revealed that diabetes, male sex, and a higher estimated glomerular filtration rate (eGFR) are important predictors of higher total and cortical volumes. Each increase of eGFR by one unit (i.e., 1 mL/min per 1.73 m2 body surface area) was associated with a 0.98 mL/m2 increase in total kidney volume, and this association was significant. Volumes were lower in persons with eGFR-defined chronic kidney disease. CONCLUSION: The extraction of image-based biomarkers through CNN-based renal sub-segmentation using data from a population-based study yields reliable results, forming a solid foundation for future investigations.

Sex differences in kidney metabolism may reflect sex-dependent outcomes in human diabetic kidney disease.

Diabetic kidney disease (DKD) is the main cause of chronic kidney disease (CKD) and progresses faster in males than in females. We identify sex-based differences in kidney metabolism and in the blood metabolome of male and female individuals with diabetes. Primary human proximal tubular epithelial cells (PTECs) from healthy males displayed increased mitochondrial respiration, oxidative stress, apoptosis, and greater injury when exposed to high glucose compared with PTECs from healthy females. Male human PTECs showed increased glucose and glutamine fluxes to the TCA cycle, whereas female human PTECs showed increased pyruvate content. The male human PTEC phenotype was enhanced by dihydrotestosterone and mediated by the transcription factor HNF4A and histone demethylase KDM6A. In mice where sex chromosomes either matched or did not match gonadal sex, male gonadal sex contributed to the kidney metabolism differences between males and females. A blood metabolomics analysis in a cohort of adolescents with or without diabetes showed increased TCA cycle metabolites in males. In a second cohort of adults with diabetes, females without DKD had higher serum pyruvate concentrations than did males with or without DKD. Serum pyruvate concentrations positively correlated with the estimated glomerular filtration rate, a measure of kidney function, and negatively correlated with all-cause mortality in this cohort. In a third cohort of adults with CKD, male sex and diabetes were associated with increased plasma TCA cycle metabolites, which correlated with all-cause mortality. These findings suggest that differences in male and female kidney metabolism may contribute to sex-dependent outcomes in DKD.

Associations of Baseline and Longitudinal Serum Uromodulin With Kidney Failure and Mortality: Results From the African American Study of Kidney Disease and Hypertension (AASK) Trial.

RATIONALE & OBJECTIVE: Uromodulin (UMOD) is the most abundant protein found in urine and has emerged as a promising biomarker of tubule health. Circulating UMOD is also detectable, but at lower levels. We evaluated whether serum UMOD levels were associated with the risks of incident kidney failure with replacement therapy (KFRT) and mortality. STUDY DESIGN: Prospective cohort. SETTING & PARTICIPANTS: Participants in AASK (the African American Study of Kidney Disease and Hypertension) with available stored serum samples from the 0-, 12-, and 24-month visits for biomarker measurement. PREDICTORS: Baseline log-transformed UMOD and change in UMOD over 2 years. OUTCOMES: KFRT and mortality. ANALYTICAL APPROACH: Cox proportional hazards and mixed-effects models. RESULTS: Among 500 participants with baseline serum UMOD levels (mean age, 54y; 37% female), 161 KFRT events occurred during a median of 8.5 years. After adjusting for baseline demographic factors, clinical factors, glomerular filtration rate, log-transformed urine protein-creatinine ratio, and randomized treatment groups, a 50% lower baseline UMOD level was independently associated with a 35% higher risk of KFRT (adjusted HR, 1.35; 95% CI, 1.07-1.70). For annual UMOD change, each 1-standard deviation lower change was associated with a 67% higher risk of KFRT (adjusted HR, 1.67; 95% CI, 1.41-1.99). Baseline UMOD and UMOD change were not associated with mortality. UMOD levels declined more steeply for metoprolol versus ramipril (P<0.001) as well as for intensive versus standard blood pressure goals (P = 0.002). LIMITATIONS: Small sample size and limited generalizability. CONCLUSIONS: Lower UMOD levels at baseline and steeper declines in UMOD over time were associated with a higher risk of subsequent KFRT in a cohort of African American adults with chronic kidney disease and hypertension. PLAIN-LANGUAGE SUMMARY: Prior studies of uromodulin (UMOD), the most abundant protein in urine, and kidney disease have focused primarily on urinary UMOD levels. The present study evaluated associations of serum UMOD levels with the risks of kidney failure with replacement therapy (KFRT) and mortality in a cohort of African American adults with hypertension and chronic kidney disease. It found that participants with lower levels of UMOD at baseline were more likely to experience KFRT even after accounting for baseline kidney measures. Similarly, participants who experienced steeper annual declines in UMOD also had a heightened risk of kidney failure. Neither baseline nor annual change in UMOD was associated with mortality. Serum UMOD is a promising biomarker of kidney health.

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

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

A systematic review of metabolomic findings in adult and pediatric renal disease.

Chronic kidney disease (CKD) affects over 0.5 billion people worldwide across their lifetimes. Despite a growingly ageing world population, an increase in all-age prevalence of kidney disease persists. Adult-onset forms of kidney disease often result from lifestyle-modifiable metabolic illnesses such as type 2 diabetes. Pediatric and adolescent forms of renal disease are primarily caused by morphological abnormalities of the kidney, as well as immunological, infectious and inherited metabolic disorders. Alterations in energy metabolism are observed in CKD of varying causes, albeit the molecular mechanisms underlying pathology are unclear. A systematic indexing of metabolites identified in plasma and urine of patients with kidney disease alongside disease enrichment analysis uncovered inborn errors of metabolism as a framework that links features of adult and pediatric kidney disease. The relationship of genetics and metabolism in kidney disease could be classified into three distinct landscapes: (i) Normal genotypes that develop renal damage because of lifestyle and / or comorbidities; (ii) Heterozygous genetic variants and polymorphisms that result in unique metabotypes that may predispose to the development of kidney disease via synergistic heterozygosity, and (iii) Homozygous genetic variants that cause renal impairment by perturbing metabolism, as found in children with monogenic inborn errors of metabolism. Interest in the identification of early biomarkers of onset and progression of CKD has grown steadily in the last years, though it has not translated into clinical routine yet. This systematic review indexes findings of differential concentration of metabolites and energy pathway dysregulation in kidney disease and appraises their potential use as biomarkers.

Mapping of the gene network that regulates glycan clock of ageing.

Glycans are an essential structural component of immunoglobulin G (IgG) that modulate its structure and function. However, regulatory mechanisms behind this complex posttranslational modification are not well known. Previous genome-wide association studies (GWAS) identified 29 genomic regions involved in regulation of IgG glycosylation, but only a few were functionally validated. One of the key functional features of IgG glycosylation is the addition of galactose (galactosylation), a trait which was shown to be associated with ageing. We performed GWAS of IgG galactosylation (N=13,705) and identified 16 significantly associated loci, indicating that IgG galactosylation is regulated by a complex network of genes that extends beyond the galactosyltransferase enzyme that adds galactose to IgG glycans. Gene prioritization identified 37 candidate genes. Using a recently developed CRISPR/dCas9 system we manipulated gene expression of candidate genes in the in vitro IgG expression system. Upregulation of three genes, EEF1A1, MANBA and TNFRSF13B, changed the IgG glycome composition, which confirmed that these three genes are involved in IgG galactosylation in this in vitro expression system.

Meta-GWAS on PCSK9 concentrations reveals associations of novel loci outside the PCSK9 locus in White populations.

BACKGROUND AND AIMS: Proprotein convertase subtilisin/kexin type 9 (PCSK9) is a key regulator of lipid homeostasis. A few earlier genome-wide association studies (GWAS) investigated genetic variants associated with circulating PCSK9 concentrations. However, uncertainty remains about some of the genetic loci discovered beyond the PCSK9 locus. By conducting the largest PCSK9 meta-analysis of GWAS (meta-GWAS) so far, we aimed to identify novel loci and validate the previously reported loci that regulate PCSK9 concentrations. METHODS: We performed GWAS for PCSK9 concentrations in two large cohorts (GCKD (n = 4,963) and KORA F3 (n = 2,895)). These were meta-analyzed with previously published data encompassing together 20,579 individuals. We further conducted a second meta-analysis in statin-naïve individuals (n = 15,390). A genetic risk score (GRS) was constructed on PCSK9-increasing SNPs and assessed its impact on the risk for coronary artery disease (CAD) in 394,943 statin-naïve participants (17,077 with events) of the UK Biobank by performing CAD-free survival analysis. RESULTS: Nine loci were genome-wide significantly associated with PCSK9 concentrations. These included the previously described PCSK9, APOB, KCNA1/KCNA5, and TM6SF2/SUGP1 loci. All imputed SNPs in the PCSK9 locus account for ∼15% of variance of PCSK9 concentrations. We further identified FADS2 as a novel locus that was also found in statin-naïve participants. All imputed SNPs within the FADS2 locus explain ∼1.2% of variance of PCSK9 concentrations. Additionally, four further loci (a region on chromosome 5, SDK1, SPATA16 and HPR) were genome-wide significant in either the main model or the statin-naïve subset. The linear increase in a PCSK9 genetic risk score was associated with 1.41-fold (95%CI 1.16-1.72, p < 0.001) higher risk for incident CAD. CONCLUSIONS: We identified five novel loci (FADS2, SPATA16, SDK1, HPR and a region on chromosome 5) for PCSK9 concentrations that would require further research. Additionally, we confirm the genome-wide significant loci that were previously detected.

Estimated Glomerular Filtration Rate, Albuminuria, and Adverse Outcomes: An Individual-Participant Data Meta-Analysis.

Importance: Chronic kidney disease (low estimated glomerular filtration rate [eGFR] or albuminuria) affects approximately 14% of adults in the US. Objective: To evaluate associations of lower eGFR based on creatinine alone, lower eGFR based on creatinine combined with cystatin C, and more severe albuminuria with adverse kidney outcomes, cardiovascular outcomes, and other health outcomes. Design, Setting, and Participants: Individual-participant data meta-analysis of 27 503 140 individuals from 114 global cohorts (eGFR based on creatinine alone) and 720 736 individuals from 20 cohorts (eGFR based on creatinine and cystatin C) and 9 067 753 individuals from 114 cohorts (albuminuria) from 1980 to 2021. Exposures: The Chronic Kidney Disease Epidemiology Collaboration 2021 equations for eGFR based on creatinine alone and eGFR based on creatinine and cystatin C; and albuminuria estimated as urine albumin to creatinine ratio (UACR). Main Outcomes and Measures: The risk of kidney failure requiring replacement therapy, all-cause mortality, cardiovascular mortality, acute kidney injury, any hospitalization, coronary heart disease, stroke, heart failure, atrial fibrillation, and peripheral artery disease. The analyses were performed within each cohort and summarized with random-effects meta-analyses. Results: Within the population using eGFR based on creatinine alone (mean age, 54 years [SD, 17 years]; 51% were women; mean follow-up time, 4.8 years [SD, 3.3 years]), the mean eGFR was 90 mL/min/1.73 m2 (SD, 22 mL/min/1.73 m2) and the median UACR was 11 mg/g (IQR, 8-16 mg/g). Within the population using eGFR based on creatinine and cystatin C (mean age, 59 years [SD, 12 years]; 53% were women; mean follow-up time, 10.8 years [SD, 4.1 years]), the mean eGFR was 88 mL/min/1.73 m2 (SD, 22 mL/min/1.73 m2) and the median UACR was 9 mg/g (IQR, 6-18 mg/g). Lower eGFR (whether based on creatinine alone or based on creatinine and cystatin C) and higher UACR were each significantly associated with higher risk for each of the 10 adverse outcomes, including those in the mildest categories of chronic kidney disease. For example, among people with a UACR less than 10 mg/g, an eGFR of 45 to 59 mL/min/1.73 m2 based on creatinine alone was associated with significantly higher hospitalization rates compared with an eGFR of 90 to 104 mL/min/1.73 m2 (adjusted hazard ratio, 1.3 [95% CI, 1.2-1.3]; 161 vs 79 events per 1000 person-years; excess absolute risk, 22 events per 1000 person-years [95% CI, 19-25 events per 1000 person-years]). Conclusions and Relevance: In this retrospective analysis of 114 cohorts, lower eGFR based on creatinine alone, lower eGFR based on creatinine and cystatin C, and more severe UACR were each associated with increased rates of 10 adverse outcomes, including adverse kidney outcomes, cardiovascular diseases, and hospitalizations.

Clonal Hematopoiesis of Indeterminate Potential is Associated with Acute Kidney Injury.

Age is a predominant risk factor for acute kidney injury (AKI), yet the biological mechanisms underlying this risk are largely unknown and to date no genetic mechanisms for AKI have been established. Clonal hematopoiesis of indeterminate potential (CHIP) is a recently recognized biological mechanism conferring risk of several chronic aging diseases including cardiovascular disease, pulmonary disease and liver disease. In CHIP, blood stem cells acquire mutations in myeloid cancer driver genes such as DNMT3A, TET2, ASXL1 and JAK2 and the myeloid progeny of these mutated cells contribute to end-organ damage through inflammatory dysregulation. We sought to establish whether CHIP causes acute kidney injury (AKI). To address this question, we first evaluated associations with incident AKI events in three population-based epidemiology cohorts (N = 442,153). We found that CHIP was associated with a greater risk of AKI (adjusted HR 1.26, 95% CI: 1.19-1.34, p<0.0001), which was more pronounced in patients with AKI requiring dialysis (adjusted HR 1.65, 95% CI: 1.24-2.20, p=0.001). The risk was particularly high in the subset of individuals where CHIP was driven by mutations in genes other than DNMT3A (HR: 1.49, 95% CI: 1.37-1.61, p<0.0001). We then examined the association between CHIP and recovery from AKI in the ASSESS-AKI cohort and identified that non-DNMT3A CHIP was more common among those with a non-resolving pattern of injury (HR 2.3, 95% CI: 1.14-4.64, p = 0.03). To gain mechanistic insight, we evaluated the role of Tet2-CHIP to AKI in ischemia-reperfusion injury (IRI) and unilateral ureteral obstruction (UUO) mouse models. In both models, we observed more severe AKI and greater post-AKI kidney fibrosis in Tet2-CHIP mice. Kidney macrophage infiltration was markedly increased in Tet2-CHIP mice and Tet2-CHIP mutant renal macrophages displayed greater proinflammatory responses. In summary, this work establishes CHIP as a genetic mechanism conferring risk of AKI and impaired kidney function recovery following AKI via an aberrant inflammatory response in CHIP derived renal macrophages.

Genetic studies of paired metabolomes reveal enzymatic and transport processes at the interface of plasma and urine.

The kidneys operate at the interface of plasma and urine by clearing molecular waste products while retaining valuable solutes. Genetic studies of paired plasma and urine metabolomes may identify underlying processes. We conducted genome-wide studies of 1,916 plasma and urine metabolites and detected 1,299 significant associations. Associations with 40% of implicated metabolites would have been missed by studying plasma alone. We detected urine-specific findings that provide information about metabolite reabsorption in the kidney, such as aquaporin (AQP)-7-mediated glycerol transport, and different metabolomic footprints of kidney-expressed proteins in plasma and urine that are consistent with their localization and function, including the transporters NaDC3 (SLC13A3) and ASBT (SLC10A2). Shared genetic determinants of 7,073 metabolite-disease combinations represent a resource to better understand metabolic diseases and revealed connections of dipeptidase 1 with circulating digestive enzymes and with hypertension. Extending genetic studies of the metabolome beyond plasma yields unique insights into processes at the interface of body compartments.

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

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

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

Background: Congenital anomalies of the kidney and urinary tract (CAKUT) are the predominant cause for chronic kidney disease below 30 years of age. Many monogenic forms have been discovered mainly due to comprehensive genetic testing like exome sequencing (ES). However, disease-causing variants in known disease-associated genes still only explain a proportion of cases. Aim of this study was to unravel the underlying molecular mechanism of syndromic CAKUT in two multiplex families with presumed autosomal recessive inheritance. Methods and Results: ES in the index individuals revealed two different rare homozygous variants in FOXD2, a transcription factor not previously implicated in CAKUT in humans: a frameshift in family 1 and a missense variant in family 2 with family segregation patterns consistent with autosomal-recessive inheritance. CRISPR/Cas9-derived Foxd2 knock-out (KO) mice presented with bilateral dilated renal pelvis accompanied by renal papilla atrophy while extrarenal features included mandibular, ophthalmologic, and behavioral anomalies, recapitulating the phenotype of humans with FOXD2 dysfunction. To study the pathomechanism of FOXD2-dysfunction-mediated developmental renal defects, in a complementary approach, we generated CRISPR/Cas9-mediated KO of Foxd2 in ureteric-bud-induced mouse metanephric mesenchyme cells. Transcriptomic analyses revealed enrichment of numerous differentially expressed genes important in renal/urogenital development, including Pax2 and Wnt4 as well as gene expression changes indicating a cell identity shift towards a stromal cell identity. Histology of Foxd2 KO mouse kidneys confirmed increased fibrosis. Further, GWAS data (genome-wide association studies) suggests that FOXD2 could play a role for maintenance of podocyte integrity during adulthood. Conclusions: In summary, our data implicate that FOXD2 dysfunction is a very rare cause of autosomal recessive syndromic CAKUT and suggest disturbances of the PAX2-WNT4 cell signaling axis contribute to this phenotype.

KLKB1 and CLSTN2 are associated with HDL-mediated cholesterol efflux capacity in a genome-wide association study.

BACKGROUND AND AIMS: HDL-mediated cholesterol efflux capacity (CEC) may protect from cardiovascular disease. Thus, we aimed to identify its genetic and non-genetic determinants. METHODS: We measured CEC to 2% apolipoprotein B-depleted serum using BODIPY-cholesterol and cAMP-stimulated J774A.1 macrophages using serum samples from 4,981 participants in the German Chronic Kidney Disease (GCKD) study. Variance of CEC explained by clinical and biochemical parameters in a multivariable linear regression model was calculated by proportional marginal variance decomposition. A genome-wide association study with 7,746,917 variants was performed based on an additive genetic model. The main model was adjusted for age, sex and principal components 1-10. Further models were selected for sensitivity analysis and to reduce residual variance by known CEC pathways. RESULTS: Variables that explained 1% and more of the variance of CEC were concentrations of triglycerides (12.9%), HDL-cholesterol (11.8%), LDL-cholesterol (3.0%), apolipoprotein A-IV (2.8%), PCSK9 (1.0%), and eGFR (1.0%). The KLKB1 (chr4) and APOE/C1 (chr19) loci were genome-wide significantly (p < 5x10-8) associated with CEC in our main model (p = 8.8x10-10 and p = 3.3x10-10, respectively). KLKB1 remained significantly associated after additional adjustment for either kidney parameters, HDL-cholesterol, triglycerides or apolipoprotein A-IV concentrations, while the APOE/C1 locus was not significantly associated anymore after adjustment for triglycerides. Adjustment for triglycerides also revealed an association with the CLSTN2 locus (chr3; p = 6.0x10-9). CONCLUSIONS: We identified HDL-cholesterol and triglycerides as the main determinants of CEC. Furthermore, we newly found a significant association of CEC with the KLKB1 and the CLSTN2 locus and confirmed the association with the APOE/C1 locus, likely mediated by triglycerides.

Epigenome-Wide Association Study Reveals CpG Sites Associated with Thyroid Function and Regulatory Effects on KLF9.

Background: Thyroid hormones play a key role in differentiation and metabolism and are known regulators of gene expression through both genomic and epigenetic processes including DNA methylation. The aim of this study was to examine associations between thyroid hormones and DNA methylation. Methods: We carried out a fixed-effect meta-analysis of epigenome-wide association study (EWAS) of blood DNA methylation sites from 8 cohorts from the ThyroidOmics Consortium, incorporating up to 7073 participants of both European and African ancestry, implementing a discovery and replication stage. Statistical analyses were conducted using normalized beta CpG values as dependent and log-transformed thyrotropin (TSH), free thyroxine, and free triiodothyronine levels, respectively, as independent variable in a linear model. The replicated findings were correlated with gene expression levels in whole blood and tested for causal influence of TSH and free thyroxine by two-sample Mendelian randomization (MR). Results: Epigenome-wide significant associations (p-value <1.1E-7) of three CpGs for free thyroxine, five for free triiodothyronine, and two for TSH concentrations were discovered and replicated (combined p-values = 1.5E-9 to 4.3E-28). The associations included CpG sites annotated to KLF9 (cg00049440) and DOT1L (cg04173586) that overlap with all three traits, consistent with hypothalamic-pituitary-thyroid axis physiology. Significant associations were also found for CpGs in FKBP5 for free thyroxine, and at CSNK1D/LINCO1970 and LRRC8D for free triiodothyronine. MR analyses supported a causal effect of thyroid status on DNA methylation of KLF9. DNA methylation of cg00049440 in KLF9 was inversely correlated with KLF9 gene expression in blood. The CpG at CSNK1D/LINC01970 overlapped with thyroid hormone receptor alpha binding peaks in liver cells. The total additive heritability of the methylation levels of the six significant CpG sites was between 25% and 57%. Significant methylation QTLs were identified for CpGs at KLF9, FKBP5, LRRC8D, and CSNK1D/LINC01970. Conclusions: We report novel associations between TSH, thyroid hormones, and blood-based DNA methylation. This study advances our understanding of thyroid hormone action particularly related to KLF9 and serves as a proof-of-concept that integrations of EWAS with other -omics data can provide a valuable tool for unraveling thyroid hormone signaling in humans by complementing and feeding classical in vitro and animal studies.

SLC26A1 is a major determinant of sulfate homeostasis in humans.

Sulfate plays a pivotal role in numerous physiological processes in the human body, including bone and cartilage health. A role of the anion transporter SLC26A1 (Sat1) for sulfate reabsorption in the kidney is supported by the observation of hyposulfatemia and hypersulfaturia in Slc26a1-knockout mice. The impact of SLC26A1 on sulfate homeostasis in humans remains to be defined. By combining clinical genetics, functional expression assays, and population exome analysis, we identify SLC26A1 as a sulfate transporter in humans and experimentally validate several loss-of-function alleles. Whole-exome sequencing from a patient presenting with painful perichondritis, hyposulfatemia, and renal sulfate wasting revealed a homozygous mutation in SLC26A1, which has not been previously described to the best of our knowledge. Whole-exome data analysis of more than 5,000 individuals confirmed that rare, putatively damaging SCL26A1 variants were significantly associated with lower plasma sulfate at the population level. Functional expression assays confirmed a substantial reduction in sulfate transport for the SLC26A1 mutation of our patient, which we consider to be novel, as well as for the additional variants detected in the population study. In conclusion, combined evidence from 3 complementary approaches supports SLC26A1 activity as a major determinant of sulfate homeostasis in humans. In view of recent evidence linking sulfate homeostasis with back pain and intervertebral disc disorder, our study identifies SLC26A1 as a potential target for modulation of musculoskeletal health.