Chronically increased echogenicity on renal ultrasound is a sensitive early finding of chronic kidney disease that can be detected before manifestation of other symptoms. Increased echogenicity, however, is not specific for a certain etiology of chronic kidney disease. Here, we performed whole exome sequencing in 79 consanguineous or familial cases of suspected nephronophthisis in order to determine the underlying molecular disease cause. In 50 cases, there was a causative mutation in a known monogenic disease gene. In 32 of these cases whole exome sequencing confirmed the diagnosis of a nephronophthisis-related ciliopathy. In 8 cases it revealed the diagnosis of a renal tubulopathy. The remaining 10 cases were identified as Alport syndrome (4), autosomal-recessive polycystic kidney disease (2), congenital anomalies of the kidney and urinary tract (3), and APECED syndrome (1). In 5 families, in whom mutations in known monogenic genes were excluded, we applied homozygosity mapping for variant filtering and identified 5 novel candidate genes (RBM48, FAM186B, PIAS1, INCENP, and RCOR1) for renal ciliopathies. Thus, whole exome sequencing allows the detection of the causative mutation in 2/3 of affected individuals, thereby presenting the etiologic diagnosis, and allows identification of novel candidate genes.

Primary ciliary dyskinesia (PCD) is caused when defects of motile cilia lead to chronic airway infections, male infertility, and situs abnormalities. Multiple causative PCD mutations account for only 65% of cases, suggesting that many genes essential for cilia function remain to be discovered. By using zebrafish morpholino knockdown of PCD candidate genes as an in vivo screening platform, we identified c21orf59, ccdc65, and c15orf26 as critical for cilia motility. c21orf59 and c15orf26 knockdown in zebrafish and planaria blocked outer dynein arm assembly, and ccdc65 knockdown altered cilia beat pattern. Biochemical analysis in Chlamydomonas revealed that the C21orf59 ortholog FBB18 is a flagellar matrix protein that accumulates specifically when cilia motility is impaired. The Chlamydomonas ida6 mutant identifies CCDC65/FAP250 as an essential component of the nexin-dynein regulatory complex. Analysis of 295 individuals with PCD identified recessive truncating mutations of C21orf59 in four families and CCDC65 in two families. Similar to findings in zebrafish and planaria, mutations in C21orf59 caused loss of both outer and inner dynein arm components. Our results characterize two genes associated with PCD-causing mutations and elucidate two distinct mechanisms critical for motile cilia function: dynein arm assembly for C21orf59 and assembly of the nexin-dynein regulatory complex for CCDC65.

Nephronophthisis is an autosomal recessive cystic kidney disease that leads to renal failure in childhood or adolescence. Most NPHP gene products form molecular networks. Here we identify ANKS6 as a new NPHP family member that connects NEK8 (NPHP9) to INVS (NPHP2) and NPHP3. We show that ANKS6 localizes to the proximal cilium and confirm its role in renal development through knockdown experiments in zebrafish and Xenopus laevis. We also identify six families with ANKS6 mutations affected by nephronophthisis, including severe cardiovascular abnormalities, liver fibrosis and situs inversus. The oxygen sensor HIF1AN hydroxylates ANKS6 and INVS and alters the composition of the ANKS6-INVS-NPHP3 module. Knockdown of Hif1an in Xenopus results in a phenotype that resembles loss of other NPHP proteins. Network analyses uncovered additional putative NPHP proteins and placed ANKS6 at the center of this NPHP module, explaining the overlapping disease manifestation caused by mutation in ANKS6, NEK8, INVS or NPHP3.

Joubert syndrome and related disorders (JSRD) are primarily autosomal-recessive conditions characterized by hypotonia, ataxia, abnormal eye movements, and intellectual disability with a distinctive mid-hindbrain malformation. Variable features include retinal dystrophy, cystic kidney disease, and liver fibrosis. JSRD are included in the rapidly expanding group of disorders called ciliopathies, because all six gene products implicated in JSRD (NPHP1, AHI1, CEP290, RPGRIP1L, TMEM67, and ARL13B) function in the primary cilium/basal body organelle. By using homozygosity mapping in consanguineous families, we identify loss-of-function mutations in CC2D2A in JSRD patients with and without retinal, kidney, and liver disease. CC2D2A is expressed in all fetal and adult tissues tested. In ciliated cells, we observe localization of recombinant CC2D2A at the basal body and colocalization with CEP290, whose cognate gene is mutated in multiple hereditary ciliopathies. In addition, the proteins can physically interact in vitro, as shown by yeast two-hybrid and GST pull-down experiments. A nonsense mutation in the zebrafish CC2D2A ortholog (sentinel) results in pronephric cysts, a hallmark of ciliary dysfunction analogous to human cystic kidney disease. Knockdown of cep290 function in sentinel fish results in a synergistic pronephric cyst phenotype, revealing a genetic interaction between CC2D2A and CEP290 and implicating CC2D2A in cilium/basal body function. These observations extend the genetic spectrum of JSRD and provide a model system for studying extragenic modifiers in JSRD and other ciliopathies.

Proteins related to the ubiquitin-proteasome system play an important role during the differentiation and ciliogenesis of photoreceptor cells. Mutations in several genes involved in ubiquitination and proteostasis have been identified as causative of inherited retinal dystrophies (IRDs) and ciliopathies. USP48 is a deubiquitinating enzyme whose role in the retina is still unexplored although previous studies indicate its relevance for neurosensory organs. In this work, we describe that a pool of endogenous USP48 localises to the basal body in retinal cells and provide data that supports the function of USP48 in the photoreceptor cilium. We also demonstrate that USP48 interacts with the IRD-associated proteins ARL3 and UNC119a, and stabilise their protein levels using different mechanisms. Our results suggest that USP48 may act in the regulation/stabilisation of key ciliary proteins for photoreceptor function, in the modulation of intracellular protein transport, and in ciliary trafficking to the photoreceptor outer segment.

Nephronophthisis-related ciliopathies (NPHP-RC) are degenerative recessive diseases that affect kidney, retina, and brain. Genetic defects in NPHP gene products that localize to cilia and centrosomes defined them as "ciliopathies." However, disease mechanisms remain poorly understood. Here, we identify by whole-exome resequencing, mutations of MRE11, ZNF423, and CEP164 as causing NPHP-RC. All three genes function within the DNA damage response (DDR) pathway. We demonstrate that, upon induced DNA damage, the NPHP-RC proteins ZNF423, CEP164, and NPHP10 colocalize to nuclear foci positive for TIP60, known to activate ATM at sites of DNA damage. We show that knockdown of CEP164 or ZNF423 causes sensitivity to DNA damaging agents and that cep164 knockdown in zebrafish results in dysregulated DDR and an NPHP-RC phenotype. Our findings link degenerative diseases of the kidney and retina, disorders of increasing prevalence, to mechanisms of DDR.

Chronic kidney disease (CKD) represents a major health burden. Its central feature of renal fibrosis is not well understood. By exome sequencing, we identified mutations in FAN1 as a cause of karyomegalic interstitial nephritis (KIN), a disorder that serves as a model for renal fibrosis. Renal histology in KIN is indistinguishable from that of nephronophthisis, except for the presence of karyomegaly. The FAN1 protein has nuclease activity and acts in DNA interstrand cross-link (ICL) repair within the Fanconi anemia DNA damage response (DDR) pathway. We show that cells from individuals with FAN1 mutations have sensitivity to the ICL-inducing agent mitomycin C but do not exhibit chromosome breakage or cell cycle arrest after diepoxybutane treatment, unlike cells from individuals with Fanconi anemia. We complemented ICL sensitivity with wild-type FAN1 but not with cDNA having mutations found in individuals with KIN. Depletion of fan1 in zebrafish caused increased DDR, apoptosis and kidney cysts. Our findings implicate susceptibility to environmental genotoxins and inadequate DNA repair as novel mechanisms contributing to renal fibrosis and CKD.

Primary ciliary dyskinesia (PCD) is a genetically heterogeneous, autosomal-recessive disorder, characterized by oto-sino-pulmonary disease and situs abnormalities. PCD-causing mutations have been identified in 20 genes, but collectively they account for only ∼65% of all PCDs. To identify mutations in additional genes that cause PCD, we performed exome sequencing on three unrelated probands with ciliary outer and inner dynein arm (ODA+IDA) defects. Mutations in SPAG1 were identified in one family with three affected siblings. Further screening of SPAG1 in 98 unrelated affected individuals (62 with ODA+IDA defects, 35 with ODA defects, 1 without available ciliary ultrastructure) revealed biallelic loss-of-function mutations in 11 additional individuals (including one sib-pair). All 14 affected individuals with SPAG1 mutations had a characteristic PCD phenotype, including 8 with situs abnormalities. Additionally, all individuals with mutations who had defined ciliary ultrastructure had ODA+IDA defects. SPAG1 was present in human airway epithelial cell lysates but was not present in isolated axonemes, and immunofluorescence staining showed an absence of ODA and IDA proteins in cilia from an affected individual, thus indicating that SPAG1 probably plays a role in the cytoplasmic assembly and/or trafficking of the axonemal dynein arms. Zebrafish morpholino studies of spag1 produced cilia-related phenotypes previously reported for PCD-causing mutations in genes encoding cytoplasmic proteins. Together, these results demonstrate that mutations in SPAG1 cause PCD with ciliary ODA+IDA defects and that exome sequencing is useful to identify genetic causes of heterogeneous recessive disorders.

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.

Primary cilia are required for Smoothened to transduce vertebrate Hedgehog signals, but how Smoothened accumulates in cilia and is activated is incompletely understood. Here, we identify cilia-associated oxysterols that promote Smoothened accumulation in cilia and activate the Hedgehog pathway. Our data reveal that cilia-associated oxysterols bind to two distinct Smoothened domains to modulate Smoothened accumulation in cilia and tune the intensity of Hedgehog pathway activation. We find that the oxysterol synthase HSD11β2 participates in the production of Smoothened-activating oxysterols and promotes Hedgehog pathway activity. Inhibiting oxysterol biosynthesis impedes oncogenic Hedgehog pathway activation and attenuates the growth of Hedgehog pathway-associated medulloblastoma, suggesting that targeted inhibition of Smoothened-activating oxysterol production may be therapeutically useful for patients with Hedgehog-associated cancers.

Cilia have a unique diffusion barrier ("gate") within their proximal region, termed transition zone (TZ), that compartmentalises signalling proteins within the organelle. The TZ is known to harbour two functional modules/complexes (Meckel syndrome [MKS] and Nephronophthisis [NPHP]) defined by genetic interaction, interdependent protein localisation (hierarchy), and proteomic studies. However, the composition and molecular organisation of these modules and their links to human ciliary disease are not completely understood. Here, we reveal Caenorhabditis elegans CEP-290 (mammalian Cep290/Mks4/Nphp6 orthologue) as a central assembly factor that is specific for established MKS module components and depends on the coiled coil region of MKS-5 (Rpgrip1L/Rpgrip1) for TZ localisation. Consistent with a critical role in ciliary gate function, CEP-290 prevents inappropriate entry of membrane-associated proteins into cilia and keeps ARL-13 (Arl13b) from leaking out of cilia via the TZ. We identify a novel MKS module component, TMEM-218 (Tmem218), that requires CEP-290 and other MKS module components for TZ localisation and functions together with the NPHP module to facilitate ciliogenesis. We show that TZ localisation of TMEM-138 (Tmem138) and CDKL-1 (Cdkl1/Cdkl2/Cdkl3/Cdlk4 related), not previously linked to a specific TZ module, similarly depends on CEP-290; surprisingly, neither TMEM-138 or CDKL-1 exhibit interdependent localisation or genetic interactions with core MKS or NPHP module components, suggesting they are part of a distinct, CEP-290-associated module. Lastly, we show that families presenting with Oral-Facial-Digital syndrome type 6 (OFD6) have likely pathogenic mutations in CEP-290-dependent TZ proteins, namely Tmem17, Tmem138, and Tmem231. Notably, patient fibroblasts harbouring mutated Tmem17, a protein not yet ciliopathy-associated, display ciliogenesis defects. Together, our findings expand the repertoire of MKS module-associated proteins--including the previously uncharacterised mammalian Tmem80--and suggest an MKS-5 and CEP-290-dependent assembly pathway for building a functional TZ.

The Meckel syndrome (MKS) complex functions at the transition zone, located between the basal body and axoneme, to regulate the localization of ciliary membrane proteins. We investigated the role of Tmem231, a two-pass transmembrane protein, in MKS complex formation and function. Consistent with a role in transition zone function, mutation of mouse Tmem231 disrupts the localization of proteins including Arl13b and Inpp5e to cilia, resulting in phenotypes characteristic of MKS such as polydactyly and kidney cysts. Tmem231 and B9d1 are essential for each other and other complex components such as Mks1 to localize to the transition zone. As in mouse, the Caenorhabditis elegans orthologue of Tmem231 localizes to and controls transition zone formation and function, suggesting an evolutionarily conserved role for Tmem231. We identified TMEM231 mutations in orofaciodigital syndrome type 3 (OFD3) and MKS patients that compromise transition zone function. Thus, Tmem231 is critical for organizing the MKS complex and controlling ciliary composition, defects in which cause OFD3 and MKS.

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.

Organoids derived from human pluripotent stem cells are a potentially powerful tool for high-throughput screening (HTS), but the complexity of organoid cultures poses a significant challenge for miniaturization and automation. Here, we present a fully automated, HTS-compatible platform for enhanced differentiation and phenotyping of human kidney organoids. The entire 21-day protocol, from plating to differentiation to analysis, can be performed automatically by liquid-handling robots, or alternatively by manual pipetting. High-content imaging analysis reveals both dose-dependent and threshold effects during organoid differentiation. Immunofluorescence and single-cell RNA sequencing identify previously undetected parietal, interstitial, and partially differentiated compartments within organoids and define conditions that greatly expand the vascular endothelium. Chemical modulation of toxicity and disease phenotypes can be quantified for safety and efficacy prediction. Screening in gene-edited organoids in this system reveals an unexpected role for myosin in polycystic kidney disease. Organoids in HTS formats thus establish an attractive platform for multidimensional phenotypic screening.

The identification of recessive disease-causing genes by homozygosity mapping is often restricted by lack of suitable consanguineous families. To overcome these limitations, we apply homozygosity mapping to single affected individuals from outbred populations. In 72 individuals of 54 kindred ascertained worldwide with known homozygous mutations in 13 different recessive disease genes, we performed total genome homozygosity mapping using 250,000 SNP arrays. Likelihood ratio Z-scores (ZLR) were plotted across the genome to detect ZLR peaks that reflect segments of homozygosity by descent, which may harbor the mutated gene. In 93% of cases, the causative gene was positioned within a consistent ZLR peak of homozygosity. The number of peaks reflected the degree of inbreeding. We demonstrate that disease-causing homozygous mutations can be detected in single cases from outbred populations within a single ZLR peak of homozygosity as short as 2 Mb, containing an average of only 16 candidate genes. As many specialty clinics have access to cohorts of individuals from outbred populations, and as our approach will result in smaller genetic candidate regions, the new strategy of homozygosity mapping in single outbred individuals will strongly accelerate the discovery of novel recessive disease genes.

The transcription factor NRF2 is a master regulator of cellular antioxidant and detoxification responses, but it also regulates other processes such as autophagy and pluripotency. In human embryonic stem cells (hESCs), NRF2 antagonizes neuroectoderm differentiation, which only occurs after NRF2 is repressed via a Primary Cilia-Autophagy-NRF2 (PAN) axis. However, the functional connections between NRF2 and primary cilia, microtubule-based plasma membrane protrusions that function as cellular antennae, remain poorly understood. For instance, nothing is known about whether NRF2 affects cilia, or whether cilia regulation of NRF2 extends beyond hESCs. Here, we show that NRF2 and primary cilia reciprocally regulate each other. First, we demonstrate that fibroblasts lacking primary cilia have higher NRF2 activity, which is rescued by autophagy-activating mTOR inhibitors, indicating that the PAN axis also operates in differentiated cells. Furthermore, NRF2 controls cilia formation and function. NRF2-null cells grow fewer and shorter cilia and display impaired Hedgehog signaling, a cilia-dependent pathway. These defects are not due to increased oxidative stress or ciliophagy, but rather to NRF2 promoting expression of multiple ciliogenic and Hedgehog pathway genes. Among these, we focused on GLI2 and GLI3, the transcription factors controlling Hh pathway output. Both their mRNA and protein levels are reduced in NRF2-null cells, consistent with their gene promoters containing consensus ARE sequences predicted to bind NRF2. Moreover, GLI2 and GLI3 fail to accumulate at the ciliary tip of NRF2-null cells upon Hh pathway activation. Given the importance of NRF2 and ciliary signaling in human disease, our data may have important biomedical implications.

Imerslund-Gräsbeck Syndrome (IGS) is a rare autosomal recessive disease characterized by intestinal vitamin B12 malabsorption. Clinical features include megaloblastic anemia, recurrent infections, failure to thrive, and proteinuria. Recessive mutations in cubilin (CUBN) and in amnionless (AMN) have been shown to cause IGS. To date, there are only about 300 cases described worldwide with only 37 different mutations found in CUBN and 30 different in the AMN gene.

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.

Hypophosphatemic rickets (HR) is a heterogeneous genetic phosphate wasting disorder. The disease is most commonly caused by mutations in the PHEX gene located on the X-chromosome or by mutations in CLCN5, DMP1, ENPP1, FGF23, and SLC34A3. The aims of this study were to perform molecular diagnostics for four patients with HR of Indian origin (two independent families) and to describe their clinical features.

Nephronophthisis-related ciliopathies (NPHP-RC) are recessive disorders that feature dysplasia or degeneration occurring preferentially in the kidney, retina and cerebellum. Here we combined homozygosity mapping with candidate gene analysis by performing 'ciliopathy candidate exome capture' followed by massively parallel sequencing. We identified 12 different truncating mutations of SDCCAG8 (serologically defined colon cancer antigen 8, also known as CCCAP) in 10 families affected by NPHP-RC. We show that SDCCAG8 is localized at both centrioles and interacts directly with OFD1 (oral-facial-digital syndrome 1), which is associated with NPHP-RC. Depletion of sdccag8 causes kidney cysts and a body axis defect in zebrafish and induces cell polarity defects in three-dimensional renal cell cultures. This work identifies loss of SDCCAG8 function as a cause of a retinal-renal ciliopathy and validates exome capture analysis for broadly heterogeneous single-gene disorders.