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.

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

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

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.

Metadherin exon 11 skipping variant enhances metastatic spread of ovarian cancer.

Metastatic ovarian cancer has a dismal prognosis and current chemotherapeutic approaches have very limited success. Metadherin (MTDH) is expressed in human ovarian cancer tissue and its expression inversely correlates with patients overall survival. Consistent with these studies, we observed MTDH expression in tissue specimens of FIGO Stage III ovarian carcinomas (72/83 cases). However, we also observed this in normal human ovarian epithelial (OE) cells, which raised the question of whether MTDH-variants with functional differences exist. We identified a novel MTDH exon 11 skipping variant (MTDHdel) which was seen at higher levels in ovarian cancer compared to benign OE cells. We analyzed MTDH-binding partner interactions and found that 12 members of the small ribosomal subunit and several mRNA binding proteins bound stronger to MTDHdel than to wildtype MTDH which indicates differential effects on gene translation. Knockdown of MTDH in ovarian cancer cells reduced the amount of distant metastases and improved the survival of ovarian cancer-bearing mice. Selective overexpression of the MTDHdel enhanced murine and human ovarian cancer progression and caused a malignant phenotype in originally benign human OE cells. MTDHdel was detectable in microdissected ovarian cancer cells of some human tissue specimens of ovarian carcinomas. In summary, we have identified a novel MTDH exon 11 skipping variant that shows enhanced binding to small ribosomal subunit members and that caused reduced overall survival of ovarian cancer bearing mice. Based on the findings in the murine system and in human tissues, MTDHdel must be considered a major promalignant factor for ovarian cancer.

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.

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.

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

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

Therapeutic Effects of Inhibition of Sphingosine-1-Phosphate Signaling in HIF-2α Inhibitor-Resistant Clear Cell Renal Cell Carcinoma.

Specific inhibitors of HIF-2α have recently been approved for the treatment of ccRCC in VHL disease patients and have shown encouraging results in clinical trials for metastatic sporadic ccRCC. However, not all patients respond to therapy and pre-clinical and clinical studies indicate that intrinsic as well as acquired resistance mechanisms to HIF-2α inhibitors are likely to represent upcoming clinical challenges. It would be desirable to have additional therapeutic options for the treatment of HIF-2α inhibitor resistant ccRCCs. Here we investigated the effects on tumor growth and on the tumor microenvironment of three different direct and indirect HIF-α inhibitors, namely the HIF-2α-specific inhibitor PT2399, the dual HIF-1α/HIF-2α inhibitor Acriflavine, and the S1P signaling pathway inhibitor FTY720, in the autochthonous Vhl/Trp53/Rb1 mutant ccRCC mouse model and validated these findings in human ccRCC cell culture models. We show that FTY720 and Acriflavine exhibit therapeutic activity in several different settings of HIF-2α inhibitor resistance. We also identify that HIF-2α inhibition strongly suppresses T cell activation in ccRCC. These findings suggest prioritization of sphingosine pathway inhibitors for clinical testing in ccRCC patients and also suggest that HIF-2α inhibitors may inhibit anti-tumor immunity and might therefore be contraindicated for combination therapies with immune checkpoint inhibitors.