An atlas of healthy and injured cell states and niches in the human kidney.

Understanding kidney disease relies on defining the complexity of cell types and states, their associated molecular profiles and interactions within tissue neighbourhoods1. Here we applied multiple single-cell and single-nucleus assays (>400,000 nuclei or cells) and spatial imaging technologies to a broad spectrum of healthy reference kidneys (45 donors) and diseased kidneys (48 patients). This has provided a high-resolution cellular atlas of 51 main cell types, which include rare and previously undescribed cell populations. The multi-omic approach provides detailed transcriptomic profiles, regulatory factors and spatial localizations spanning the entire kidney. We also define 28 cellular states across nephron segments and interstitium that were altered in kidney injury, encompassing cycling, adaptive (successful or maladaptive repair), transitioning and degenerative states. Molecular signatures permitted the localization of these states within injury neighbourhoods using spatial transcriptomics, while large-scale 3D imaging analysis (around 1.2 million neighbourhoods) provided corresponding linkages to active immune responses. These analyses defined biological pathways that are relevant to injury time-course and niches, including signatures underlying epithelial repair that predicted maladaptive states associated with a decline in kidney function. This integrated multimodal spatial cell atlas of healthy and diseased human kidneys represents a comprehensive benchmark of cellular states, neighbourhoods, outcome-associated signatures and publicly available interactive visualizations.

Mapping the NPHP-JBTS-MKS protein network reveals ciliopathy disease genes and pathways.

Nephronophthisis (NPHP), Joubert (JBTS), and Meckel-Gruber (MKS) syndromes are autosomal-recessive ciliopathies presenting with cystic kidneys, retinal degeneration, and cerebellar/neural tube malformation. Whether defects in kidney, retinal, or neural disease primarily involve ciliary, Hedgehog, or cell polarity pathways remains unclear. Using high-confidence proteomics, we identified 850 interactors copurifying with nine NPHP/JBTS/MKS proteins and discovered three connected modules: "NPHP1-4-8" functioning at the apical surface, "NPHP5-6" at centrosomes, and "MKS" linked to Hedgehog signaling. Assays for ciliogenesis and epithelial morphogenesis in 3D renal cultures link renal cystic disease to apical organization defects, whereas ciliary and Hedgehog pathway defects lead to retinal or neural deficits. Using 38 interactors as candidates, linkage and sequencing analysis of 250 patients identified ATXN10 and TCTN2 as new NPHP-JBTS genes, and our Tctn2 mouse knockout shows neural tube and Hedgehog signaling defects. Our study further illustrates the power of linking proteomic networks and human genetics to uncover critical disease pathways.

Pax protein depletion in proximal tubules triggers conserved mechanisms of resistance to acute ischemic kidney injury preventing transition to chronic kidney disease.

Acute kidney injury (AKI) is a common condition that lacks effective treatments. In part, this shortcoming is due to an incomplete understanding of the genetic mechanisms that control pathogenesis and recovery. Identifying the molecular and genetic regulators unique to nephron segments that dictate vulnerability to injury and regenerative potential could lead to new therapeutic targets to treat ischemic kidney injury. Pax2 and Pax8 are homologous transcription factors with overlapping functions that are critical for kidney development and are re-activated in AKI. Here, we examined the role of Pax2 and Pax8 in recovery from ischemic AKI and found them upregulated after severe AKI and correlated with chronic injury. Surprisingly, proximal-tubule-selective deletion of Pax2 and Pax8 resulted in a less severe chronic injury phenotype. This effect was mediated by protection against the acute insult, similar to pre-conditioning. Prior to injury, Pax2 and Pax8 mutant mice develop a unique subpopulation of proximal tubule cells in the S3 segment that displayed features usually seen only in acute or chronic injury. The expression signature of these cells was strongly enriched with genes associated with other mechanisms of protection against ischemic AKI including caloric restriction, hypoxic pre-conditioning, and female sex. Thus, our results identified a novel role for Pax2 and Pax8 in mature proximal tubules that regulates critical genes and pathways involved in both the injury response and protection from ischemic AKI.

Sodium glucose co-transporter 2 inhibition increases epidermal growth factor expression and improves outcomes in patients with type 2 diabetes.

Underlying molecular mechanisms of the kidney protective effects of sodium glucose co-transporter 2 (SGLT2) inhibitors are not fully elucidated. Therefore, we studied the association between urinary epidermal growth factor (uEGF), a mitogenic factor involved in kidney repair, and kidney outcomes in patients with type 2 diabetes (T2D). The underlying molecular mechanisms of the SGLT2 inhibitor canagliflozin on EGF using single-cell RNA sequencing from kidney tissue were examined. Urinary EGF-to-creatinine ratio (uEGF/Cr) was measured in 3521 CANagliflozin cardioVascular Assessment Study (CANVAS) participants at baseline and week 52. Associations of uEGF/Cr with kidney outcome were assessed using multivariable-adjusted Cox regression models. Single-cell RNA sequencing was performed using protocol kidney biopsy tissue from ten young patients with T2D on SGLT2i, six patients with T2D on standard care only, and six healthy controls (HCs). In CANVAS, each doubling in baseline uEGF/Cr was associated with a 12% (95% confidence interval 1-22) decreased risk of kidney outcome. uEGF/Cr decreased after 52 weeks with placebo and remained stable with canagliflozin (between-group difference +7.3% (2.0-12.8). In young persons with T2D, EGF mRNA was primarily expressed in the thick ascending loop of Henle. Expression in biopsies from T2D without SGLT2i was significantly lower compared to HCs, whereas treatment with SGLT2i increased EGF levels closer to the healthy state. In young persons with T2D without SGLT2i, endothelin-1 emerged as a key regulator of the EGF co-expression network. SGLT2i treatment was associated with a shift towards normal EGF expression. Thus, decreased uEGF represents increased risk of kidney disease progression in patients with T2D. Canagliflozin increased kidney tissue expression of EGF and was associated with a downstream signaling cascade linked to tubular repair and reversal of tubular injury.

Glomerular endothelial cell-podocyte stresses and crosstalk in structurally normal kidney transplants.

Increased podocyte detachment begins immediately after kidney transplantation and is associated with long-term allograft failure. We hypothesized that cell-specific transcriptional changes in podocytes and glomerular endothelial cells after transplantation would offer mechanistic insights into the podocyte detachment process. To test this, we evaluated cell-specific transcriptional profiles of glomerular endothelial cells and podocytes from 14 patients of their first-year surveillance biopsies with normal histology from low immune risk recipients with no post-transplant complications and compared these to biopsies of 20 healthy living donor controls. Glomerular endothelial cells from these surveillance biopsies were enriched for genes related to fluid shear stress, angiogenesis, and interferon signaling. In podocytes, pathways were enriched for genes in response to growth factor signaling and actin cytoskeletal reorganization but also showed evidence of podocyte stress as indicated by reduced nephrin (adhesion protein) gene expression. In parallel, transcripts coding for proteins required to maintain podocyte adherence to the underlying glomerular basement membrane were downregulated, including the major glomerular podocyte integrin α3 and the actin cytoskeleton-related gene synaptopodin. The reduction in integrin α3 protein expression in surveillance biopsies was confirmed by immunoperoxidase staining. The combined growth and stress response of patient allografts post-transplantation paralleled similar changes in a rodent model of nephrectomy-induced glomerular hypertrophic stress that progress to develop proteinuria and glomerulosclerosis with shortened kidney life span. Thus, even among patients with apparently healthy allografts with no detectable histologic abnormality including alloimmune injury, transcriptomic changes reflecting cell stresses are already set in motion that could drive hypertrophy-associated glomerular disease progression.

A multimodal and integrated approach to interrogate human kidney biopsies with rigor and reproducibility: guidelines from the Kidney Precision Medicine Project.

Comprehensive and spatially mapped molecular atlases of organs at a cellular level are a critical resource to gain insights into pathogenic mechanisms and personalized therapies for diseases. The Kidney Precision Medicine Project (KPMP) is an endeavor to generate three-dimensional (3-D) molecular atlases of healthy and diseased kidney biopsies by using multiple state-of-the-art omics and imaging technologies across several institutions. Obtaining rigorous and reproducible results from disparate methods and at different sites to interrogate biomolecules at a single-cell level or in 3-D space is a significant challenge that can be a futile exercise if not well controlled. We describe a "follow the tissue" pipeline for generating a reliable and authentic single-cell/region 3-D molecular atlas of human adult kidney. Our approach emphasizes quality assurance, quality control, validation, and harmonization across different omics and imaging technologies from sample procurement, processing, storage, shipping to data generation, analysis, and sharing. We established benchmarks for quality control, rigor, reproducibility, and feasibility across multiple technologies through a pilot experiment using common source tissue that was processed and analyzed at different institutions and different technologies. A peer review system was established to critically review quality control measures and the reproducibility of data generated by each technology before their being approved to interrogate clinical biopsy specimens. The process established economizes the use of valuable biopsy tissue for multiomics and imaging analysis with stringent quality control to ensure rigor and reproducibility of results and serves as a model for precision medicine projects across laboratories, institutions and consortia.

An eQTL Landscape of Kidney Tissue in Human Nephrotic Syndrome.

Expression quantitative trait loci (eQTL) studies illuminate the genetics of gene expression and, in disease research, can be particularly illuminating when using the tissues directly impacted by the condition. In nephrology, there is a paucity of eQTL studies of human kidney. Here, we used whole-genome sequencing (WGS) and microdissected glomerular (GLOM) and tubulointerstitial (TI) transcriptomes from 187 individuals with nephrotic syndrome (NS) to describe the eQTL landscape in these functionally distinct kidney structures. Using MatrixEQTL, we performed cis-eQTL analysis on GLOM (n = 136) and TI (n = 166). We used the Bayesian "Deterministic Approximation of Posteriors" (DAP) to fine-map these signals, eQTLBMA to discover GLOM- or TI-specific eQTLs, and single-cell RNA-seq data of control kidney tissue to identify the cell type specificity of significant eQTLs. We integrated eQTL data with an IgA Nephropathy (IgAN) GWAS to perform a transcriptome-wide association study (TWAS). We discovered 894 GLOM eQTLs and 1,767 TI eQTLs at FDR < 0.05. 14% and 19% of GLOM and TI eQTLs, respectively, had >1 independent signal associated with its expression. 12% and 26% of eQTLs were GLOM specific and TI specific, respectively. GLOM eQTLs were most significantly enriched in podocyte transcripts and TI eQTLs in proximal tubules. The IgAN TWAS identified significant GLOM and TI genes, primarily at the HLA region. In this study, we discovered GLOM and TI eQTLs, identified those that were tissue specific, deconvoluted them into cell-specific signals, and used them to characterize known GWAS alleles. These data are available for browsing and download via our eQTL browser, "nephQTL."

Single-cell analysis of progenitor cell dynamics and lineage specification in the human fetal kidney.

The mammalian kidney develops through reciprocal interactions between the ureteric bud and the metanephric mesenchyme to give rise to the entire collecting system and the nephrons. Most of our knowledge of the developmental regulators driving this process arises from the study of gene expression and functional genetics in mice and other animal models. In order to shed light on human kidney development, we have used single-cell transcriptomics to characterize gene expression in different cell populations, and to study individual cell dynamics and lineage trajectories during development. Single-cell transcriptome analyses of 6414 cells from five individual specimens identified 11 initial clusters of specific renal cell types as defined by their gene expression profile. Further subclustering identifies progenitors, and mature and intermediate stages of differentiation for several renal lineages. Other lineages identified include mesangium, stroma, endothelial and immune cells. Novel markers for these cell types were revealed in the analysis, as were components of key signaling pathways driving renal development in animal models. Altogether, we provide a comprehensive and dynamic gene expression profile of the developing human kidney at the single-cell level.

Hypertension induces glomerulosclerosis in phospholipase C-ε1 deficiency.

Loss-of-function mutations in phospholipase C-ε1 (PLCE1) have been detected in patients with nephrotic syndrome, but other family members with the same mutation were asymptomatic, suggesting additional stressor are required to cause the full phenotype. Consistent with these observations, we determined that global Plce1-deficient mice have histologically normal glomeruli and no albuminuria at baseline. Angiotensin II (ANG II) is known to induce glomerular damage in genetically susceptible individuals. Therefore, we tested whether ANG II enhances glomerular damage in Plce1-deficient mice. ANG II increased blood pressure equally in Plce1-deficient and wild-type littermates. Additionally, it led to 20-fold increased albuminuria and significantly more sclerotic glomeruli in Plce1-deficient mice compared with wild-type littermates. Furthermore, Plce1-deficient mice demonstrated diffuse mesangial expansion, podocyte loss, and focal podocyte foot process effacement. To determine whether these effects are mediated by hypertension and hyperfiltration, rather than directly through ANG II, we raised blood pressure to a similar level using DOCA + salt + uninephrectomy and norepinephrine. This caused a fivefold increase in albuminuria in Plce1-deficient mice and a significant increase in the number of sclerotic glomeruli. Consistent with previous findings in mice, we detected strong PLCE1 transcript expression in podocytes using single cell sequencing of human kidney tissue. In hemagglutinin-tagged Plce1 transgenic mice, Plce1 was detected in podocytes and also in glomerular arterioles using immunohistochemistry. Our data demonstrate that Plce1 deficiency in mice predisposes to glomerular damage secondary to hypertensive insults.

SARS-CoV-2 receptor networks in diabetic and COVID-19-associated kidney disease.

COVID-19 morbidity and mortality are increased via unknown mechanisms in patients with diabetes and kidney disease. SARS-CoV-2 uses angiotensin-converting enzyme 2 (ACE2) for entry into host cells. Because ACE2 is a susceptibility factor for infection, we investigated how diabetic kidney disease and medications alter ACE2 receptor expression in kidneys. Single cell RNA profiling of kidney biopsies from healthy living donors and patients with diabetic kidney disease revealed ACE2 expression primarily in proximal tubular epithelial cells. This cell-specific localization was confirmed by in situ hybridization. ACE2 expression levels were unaltered by exposures to renin-angiotensin-aldosterone system inhibitors in diabetic kidney disease. Bayesian integrative analysis of a large compendium of public -omics datasets identified molecular network modules induced in ACE2-expressing proximal tubular epithelial cells in diabetic kidney disease (searchable at hb.flatironinstitute.org/covid-kidney) that were linked to viral entry, immune activation, endomembrane reorganization, and RNA processing. The diabetic kidney disease ACE2-positive proximal tubular epithelial cell module overlapped with expression patterns seen in SARS-CoV-2-infected cells. Similar cellular programs were seen in ACE2-positive proximal tubular epithelial cells obtained from urine samples of 13 hospitalized patients with COVID-19, suggesting a consistent ACE2-coregulated proximal tubular epithelial cell expression program that may interact with the SARS-CoV-2 infection processes. Thus SARS-CoV-2 receptor networks can seed further research into risk stratification and therapeutic strategies for COVID-19-related kidney damage.

Exome-wide Association Study Identifies GREB1L Mutations in Congenital Kidney Malformations.

Renal agenesis and hypodysplasia (RHD) are major causes of pediatric chronic kidney disease and are highly genetically heterogeneous. We conducted whole-exome sequencing in 202 case subjects with RHD and identified diagnostic mutations in genes known to be associated with RHD in 7/202 case subjects. In an additional affected individual with RHD and a congenital heart defect, we found a homozygous loss-of-function (LOF) variant in SLIT3, recapitulating phenotypes reported with Slit3 inactivation in the mouse. To identify genes associated with RHD, we performed an exome-wide association study with 195 unresolved case subjects and 6,905 control subjects. The top signal resided in GREB1L, a gene implicated previously in Hoxb1 and Shha signaling in zebrafish. The significance of the association, which was p = 2.0 × 10-5 for novel LOF, increased to p = 4.1 × 10-6 for LOF and deleterious missense variants combined, and augmented further after accounting for segregation and de novo inheritance of rare variants (joint p = 2.3 × 10-7). Finally, CRISPR/Cas9 disruption or knockdown of greb1l in zebrafish caused specific pronephric defects, which were rescued by wild-type human GREB1L mRNA, but not mRNA containing alleles identified in case subjects. Together, our study provides insight into the genetic landscape of kidney malformations in humans, presents multiple candidates, and identifies SLIT3 and GREB1L as genes implicated in the pathogenesis of RHD.

Genetic Drivers of Kidney Defects in the DiGeorge Syndrome.

BACKGROUND: The DiGeorge syndrome, the most common of the microdeletion syndromes, affects multiple organs, including the heart, the nervous system, and the kidney. It is caused by deletions on chromosome 22q11.2; the genetic driver of the kidney defects is unknown. METHODS: We conducted a genomewide search for structural variants in two cohorts: 2080 patients with congenital kidney and urinary tract anomalies and 22,094 controls. We performed exome and targeted resequencing in samples obtained from 586 additional patients with congenital kidney anomalies. We also carried out functional studies using zebrafish and mice. RESULTS: We identified heterozygous deletions of 22q11.2 in 1.1% of the patients with congenital kidney anomalies and in 0.01% of population controls (odds ratio, 81.5; P=4.5×10-14). We localized the main drivers of renal disease in the DiGeorge syndrome to a 370-kb region containing nine genes. In zebrafish embryos, an induced loss of function in snap29, aifm3, and crkl resulted in renal defects; the loss of crkl alone was sufficient to induce defects. Five of 586 patients with congenital urinary anomalies had newly identified, heterozygous protein-altering variants, including a premature termination codon, in CRKL. The inactivation of Crkl in the mouse model induced developmental defects similar to those observed in patients with congenital urinary anomalies. CONCLUSIONS: We identified a recurrent 370-kb deletion at the 22q11.2 locus as a driver of kidney defects in the DiGeorge syndrome and in sporadic congenital kidney and urinary tract anomalies. Of the nine genes at this locus, SNAP29, AIFM3, and CRKL appear to be critical to the phenotype, with haploinsufficiency of CRKL emerging as the main genetic driver. (Funded by the National Institutes of Health and others.).

Glycine Amidinotransferase (GATM), Renal Fanconi Syndrome, and Kidney Failure.

Background For many patients with kidney failure, the cause and underlying defect remain unknown. Here, we describe a novel mechanism of a genetic order characterized by renal Fanconi syndrome and kidney failure.Methods We clinically and genetically characterized members of five families with autosomal dominant renal Fanconi syndrome and kidney failure. We performed genome-wide linkage analysis, sequencing, and expression studies in kidney biopsy specimens and renal cells along with knockout mouse studies and evaluations of mitochondrial morphology and function. Structural studies examined the effects of recognized mutations.Results The renal disease in these patients resulted from monoallelic mutations in the gene encoding glycine amidinotransferase (GATM), a renal proximal tubular enzyme in the creatine biosynthetic pathway that is otherwise associated with a recessive disorder of creatine deficiency. In silico analysis showed that the particular GATM mutations, identified in 28 members of the five families, create an additional interaction interface within the GATM protein and likely cause the linear aggregation of GATM observed in patient biopsy specimens and cultured proximal tubule cells. GATM aggregates-containing mitochondria were elongated and associated with increased ROS production, activation of the NLRP3 inflammasome, enhanced expression of the profibrotic cytokine IL-18, and increased cell death.Conclusions In this novel genetic disorder, fully penetrant heterozygous missense mutations in GATM trigger intramitochondrial fibrillary deposition of GATM and lead to elongated and abnormal mitochondria. We speculate that this renal proximal tubular mitochondrial pathology initiates a response from the inflammasome, with subsequent development of kidney fibrosis.

DCDC2 mutations cause a renal-hepatic ciliopathy by disrupting Wnt signaling.

Nephronophthisis-related ciliopathies (NPHP-RC) are recessive diseases characterized by renal dysplasia or degeneration. We here identify mutations of DCDC2 as causing a renal-hepatic ciliopathy. DCDC2 localizes to the ciliary axoneme and to mitotic spindle fibers in a cell-cycle-dependent manner. Knockdown of Dcdc2 in IMCD3 cells disrupts ciliogenesis, which is rescued by wild-type (WT) human DCDC2, but not by constructs that reflect human mutations. We show that DCDC2 interacts with DVL and DCDC2 overexpression inhibits ß-catenin-dependent Wnt signaling in an effect additive to Wnt inhibitors. Mutations detected in human NPHP-RC lack these effects. A Wnt inhibitor likewise restores ciliogenesis in 3D IMCD3 cultures, emphasizing the importance of Wnt signaling for renal tubulogenesis. Knockdown of dcdc2 in zebrafish recapitulates NPHP-RC phenotypes, including renal cysts and hydrocephalus, which is rescued by a Wnt inhibitor and by WT, but not by mutant, DCDC2. We thus demonstrate a central role of Wnt signaling in the pathogenesis of NPHP-RC, suggesting an avenue for potential treatment of NPHP-RC.

Mutations in EMP2 cause childhood-onset nephrotic syndrome.

Nephrotic syndrome (NS) is a genetically heterogeneous group of diseases that are divided into steroid-sensitive NS (SSNS) and steroid-resistant NS (SRNS). SRNS inevitably leads to end-stage kidney disease, and no curative treatment is available. To date, mutations in more than 24 genes have been described in Mendelian forms of SRNS; however, no Mendelian form of SSNS has been described. To identify a genetic form of SSNS, we performed homozygosity mapping, whole-exome sequencing, and multiplex PCR followed by next-generation sequencing. We thereby detected biallelic mutations in EMP2 (epithelial membrane protein 2) in four individuals from three unrelated families affected by SRNS or SSNS. We showed that EMP2 exclusively localized to glomeruli in the kidney. Knockdown of emp2 in zebrafish resulted in pericardial effusion, supporting the pathogenic role of mutated EMP2 in human NS. At the cellular level, we showed that knockdown of EMP2 in podocytes and endothelial cells resulted in an increased amount of CAVEOLIN-1 and decreased cell proliferation. Our data therefore identify EMP2 mutations as causing a recessive Mendelian form of SSNS.

Mutations of CEP83 cause infantile nephronophthisis and intellectual disability.

Ciliopathies are a group of hereditary disorders associated with defects in cilia structure and function. The distal appendages (DAPs) of centrioles are involved in the docking and anchoring of the mother centriole to the cellular membrane during ciliogenesis. The molecular composition of DAPs was recently elucidated and mutations in two genes encoding DAPs components (CEP164/NPHP15, SCLT1) have been associated with human ciliopathies, namely nephronophthisis and orofaciodigital syndrome. To identify additional DAP components defective in ciliopathies, we independently performed targeted exon sequencing of 1,221 genes associated with cilia and 5 known DAP protein-encoding genes in 1,255 individuals with a nephronophthisis-related ciliopathy. We thereby detected biallelic mutations in a key component of DAP-encoding gene, CEP83, in seven families. All affected individuals had early-onset nephronophthisis and four out of eight displayed learning disability and/or hydrocephalus. Fibroblasts and tubular renal cells from affected individuals showed an altered DAP composition and ciliary defects. In summary, we have identified mutations in CEP83, another DAP-component-encoding gene, as a cause of infantile nephronophthisis associated with central nervous system abnormalities in half of the individuals.

Evaluating Mendelian nephrotic syndrome genes for evidence for risk alleles or oligogenicity that explain heritability.

BACKGROUND: More than 30 genes can harbor rare exonic variants sufficient to cause nephrotic syndrome (NS), and the number of genes implicated in monogenic NS continues to grow. However, outside the first year of life, the majority of affected patients, particularly in ancestrally mixed populations, do not have a known monogenic form of NS. Even in those children classified with a monogenic form of NS, there is phenotypic heterogeneity. Thus, we have only discovered a fraction of the heritability of NS-the underlying genetic factors contributing to phenotypic variation. Part of the "missing heritability" for NS has been posited to be explained by patients harboring coding variants across one or more previously implicated NS genes, insufficient to cause NS in a classical Mendelian manner, but that nonetheless have a sufficient impact on protein function to cause disease. However, systematic evaluation in patients with NS for rare or low-frequency risk alleles within single genes, or in combination across genes ("oligogenicity"), has not been reported. To determine whether, compared with a reference population, patients with NS have either a significantly increased burden of protein-altering variants ("risk-alleles"), or a unique combination of them ("oligogenicity"), in a set of 21 genes implicated in Mendelian forms of NS. METHODS: In 303 patients with NS enrolled in the Nephrotic Syndrome Study Network (NEPTUNE), we performed targeted amplification paired with next-generation sequencing of 21 genes implicated in monogenic NS. We created a high-quality variant call set and compared it with a variant call set of the same genes in a reference population composed of 2,535 individuals from phase 3 of the 1000 Genomes Project. We created both a "stringent" and a "relaxed" pathogenicity-filtering pipeline, applied them to both cohorts, and computed the burden of variants in the entire gene set per cohort, the burden of variants in the entire gene set per individual, the burden of variants within a single gene per cohort, and unique combinations of variants across two or more genes per cohort. RESULTS: With few exceptions when using the relaxed filter, and which are likely the result of confounding by population stratification, NS patients did not have a significantly increased burden of variants in Mendelian NS genes in comparison to a reference cohort, nor was there any evidence for oligogenicity. This was true when using both the relaxed and the stringent variant pathogenicity filter. CONCLUSION: In our study, there were no significant differences in the burden or particular combinations of low-frequency or rare protein-altering variants in a previously implicated Mendelian NS genes cohort between North American patients with NS and a reference population. Studies in larger independent cohorts or meta-analyses are needed to assess the generalizability of our discoveries and also address whether there is in fact small but significant enrichment of risk alleles or oligogenicity in NS cases that was undetectable with this current sample size. It is still possible that rare protein-altering variants in these genes, insufficient to cause Mendelian disease, still contribute to NS as risk alleles and/or via oligogenicity. However, we suggest that more accurate bioinformatic analyses and the incorporation of functional assays would be necessary to identify bona fide instances of this form of genetic architecture as a contributor to the heritability of NS.

Large-scale targeted sequencing comparison highlights extreme genetic heterogeneity in nephronophthisis-related ciliopathies.

BACKGROUND: The term nephronophthisis-related ciliopathies (NPHP-RC) describes a group of rare autosomal-recessive cystic kidney diseases, characterised by broad genetic and clinical heterogeneity. NPHP-RC is frequently associated with extrarenal manifestations and accounts for the majority of genetically caused chronic kidney disease (CKD) during childhood and adolescence. Generation of a molecular diagnosis has been impaired by this broad genetic heterogeneity. However, recently developed high-throughput exon sequencing techniques represent powerful and efficient tools to screen large cohorts for dozens of causative genes. METHODS: Therefore, we performed massively multiplexed targeted sequencing using the modified molecular inversion probe strategy (MIPs) in an international cohort of 384 patients diagnosed with NPHP-RC. RESULTS: As a result, we established the molecular diagnoses in 81/384 unrelated individuals (21.1%). We detected 127 likely disease-causing mutations in 18 of 34 evaluated NPHP-RC genes, 22 of which were novel. We further compared a subgroup of current findings to the results of a previous study in which we used an array-based microfluidic PCR technology in the same cohort. While 78 likely disease-causing mutations were previously detected by the array-based microfluidic PCR, the MIPs approach identified 94 likely pathogenic mutations. Compared with the previous approach, MIPs redetected 66 out of 78 variants and 28 previously unidentified variants, for a total of 94 variants. CONCLUSIONS: In summary, we demonstrate that the modified MIPs technology is a useful approach to screen large cohorts for a multitude of established NPHP genes in order to identify the underlying molecular cause. Combined application of two independent library preparation and sequencing techniques, however, may still be indicated for Mendelian diseases with extensive genetic heterogeneity in order to further increase diagnostic sensitivity.

MKS1 regulates ciliary INPP5E levels in Joubert syndrome.

BACKGROUND: Joubert syndrome (JS) is a recessive ciliopathy characterised by a distinctive brain malformation 'the molar tooth sign'. Mutations in >27 genes cause JS, and mutations in 12 of these genes also cause Meckel-Gruber syndrome (MKS). The goals of this work are to describe the clinical features of MKS1-related JS and determine whether disease causing MKS1 mutations affect cellular phenotypes such as cilium number, length and protein content as potential mechanisms underlying JS. METHODS: We measured cilium number, length and protein content (ARL13B and INPP5E) by immunofluorescence in fibroblasts from individuals with MKS1-related JS and in a three-dimensional (3D) spheroid rescue assay to test the effects of disease-related MKS1 mutations. RESULTS: We report MKS1 mutations (eight of them previously unreported) in nine individuals with JS. A minority of the individuals with MKS1-related JS have MKS features. In contrast to the truncating mutations associated with MKS, all of the individuals with MKS1-related JS carry ≥ 1 non-truncating mutation. Fibroblasts from individuals with MKS1-related JS make normal or fewer cilia than control fibroblasts, their cilia are more variable in length than controls, and show decreased ciliary ARL13B and INPP5E. Additionally, MKS1 mutant alleles have similar effects in 3D spheroids. CONCLUSIONS: MKS1 functions in the transition zone at the base of the cilium to regulate ciliary INPP5E content, through an ARL13B-dependent mechanism. Mutations in INPP5E also cause JS, so our findings in patient fibroblasts support the notion that loss of INPP5E function, due to either mutation or mislocalisation, is a key mechanism underlying JS, downstream of MKS1 and ARL13B.