Variants in tubule epithelial regulatory elements mediate most heritable differences in human kidney function.

Kidney failure, the decrease of kidney function below a threshold necessary to support life, is a major cause of morbidity and mortality. We performed a genome-wide association study (GWAS) of 406,504 individuals in the UK Biobank, identifying 430 loci affecting kidney function in middle-aged adults. To investigate the cell types affected by these loci, we integrated the GWAS with human kidney candidate cis-regulatory elements (cCREs) identified using single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq). Overall, 56% of kidney function heritability localized to kidney tubule epithelial cCREs and an additional 7% to kidney podocyte cCREs. Thus, most heritable differences in adult kidney function are a result of altered gene expression in these two cell types. Using enhancer assays, allele-specific scATAC-seq and machine learning, we found that many kidney function variants alter tubule epithelial cCRE chromatin accessibility and function. Using CRISPRi, we determined which genes some of these cCREs regulate, implicating NDRG1, CCNB1 and STC1 in human kidney function.

Variants in tubule epithelial regulatory elements mediate most heritable differences in human kidney function.

Kidney disease is highly heritable; however, the causal genetic variants, the cell types in which these variants function, and the molecular mechanisms underlying kidney disease remain largely unknown. To identify genetic loci affecting kidney function, we performed a GWAS using multiple kidney function biomarkers and identified 462 loci. To begin to investigate how these loci affect kidney function, we generated single-cell chromatin accessibility (scATAC-seq) maps of the human kidney and identified candidate cis-regulatory elements (cCREs) for kidney podocytes, tubule epithelial cells, and kidney endothelial, stromal, and immune cells. Kidney tubule epithelial cCREs explained 58% of kidney function SNP-heritability and kidney podocyte cCREs explained an additional 6.5% of SNP-heritability. In contrast, little kidney function heritability was explained by kidney endothelial, stromal, or immune cell-specific cCREs. Through functionally informed fine-mapping, we identified putative causal kidney function variants and their corresponding cCREs. Using kidney scATAC-seq data, we created a deep learning model (which we named ChromKid) to predict kidney cell type-specific chromatin accessibility from sequence. ChromKid and allele specific kidney scATAC-seq revealed that many fine-mapped kidney function variants locally change chromatin accessibility in tubule epithelial cells. Enhancer assays confirmed that fine-mapped kidney function variants alter tubule epithelial regulatory element function. To map the genes which these regulatory elements control, we used CRISPR interference (CRISPRi) to target these regulatory elements in tubule epithelial cells and assessed changes in gene expression. CRISPRi of enhancers harboring kidney function variants regulated NDRG1 and RBPMS expression. Thus, inherited differences in tubule epithelial NDRG1 and RBPMS expression may predispose to kidney disease in humans. We conclude that genetic variants affecting tubule epithelial regulatory element function account for most SNP-heritability of human kidney function. This work provides an experimental approach to identify the variants, regulatory elements, and genes involved in polygenic disease.

Cilia-enriched oxysterol 7ß,27-DHC is required for polycystin ion channel activation.

Polycystin-1 (PC-1) and PC-2 form a heteromeric ion channel complex that is abundantly expressed in primary cilia of renal epithelial cells. This complex functions as a non-selective cation channel, and mutations within the polycystin complex cause autosomal dominant polycystic kidney disease (ADPKD). The spatial and temporal regulation of the polycystin complex within the ciliary membrane remains poorly understood. Using both whole-cell and ciliary patch-clamp recordings, we identify a cilia-enriched oxysterol, 7ß,27-dihydroxycholesterol (DHC), that serves as a necessary activator of the polycystin complex. We further identify an oxysterol-binding pocket within PC-2 and showed that mutations within this binding pocket disrupt 7ß,27-DHC-dependent polycystin activation. Pharmacologic and genetic inhibition of oxysterol synthesis reduces channel activity in primary cilia. In summary, our findings reveal a regulator of the polycystin complex. This oxysterol-binding pocket in PC-2 may provide a specific target for potential ADPKD therapeutics.

An alternative cell cycle coordinates multiciliated cell differentiation.

The canonical mitotic cell cycle coordinates DNA replication, centriole duplication and cytokinesis to generate two cells from one1. Some cells, such as mammalian trophoblast giant cells, use cell cycle variants like the endocycle to bypass mitosis2. Differentiating multiciliated cells, found in the mammalian airway, brain ventricles and reproductive tract, are post-mitotic but generate hundreds of centrioles, each of which matures into a basal body and nucleates a motile cilium3,4. Several cell cycle regulators have previously been implicated in specific steps of multiciliated cell differentiation5,6. Here we show that differentiating multiciliated cells integrate cell cycle regulators into a new alternative cell cycle, which we refer to as the multiciliation cycle. The multiciliation cycle redeploys many canonical cell cycle regulators, including cyclin-dependent kinases (CDKs) and their cognate cyclins. For example, cyclin D1, CDK4 and CDK6, which are regulators of mitotic G1-to-S progression, are required to initiate multiciliated cell differentiation. The multiciliation cycle amplifies some aspects of the canonical cell cycle, such as centriole synthesis, and blocks others, such as DNA replication. E2F7, a transcriptional regulator of canonical S-to-G2 progression, is expressed at high levels during the multiciliation cycle. In the multiciliation cycle, E2F7 directly dampens the expression of genes encoding DNA replication machinery and terminates the S phase-like gene expression program. Loss of E2F7 causes aberrant acquisition of DNA synthesis in multiciliated cells and dysregulation of multiciliation cycle progression, which disrupts centriole maturation and ciliogenesis. We conclude that multiciliated cells use an alternative cell cycle that orchestrates differentiation instead of controlling proliferation.

Emerging mechanistic understanding of cilia function in cellular signalling.

Primary cilia are solitary, immotile sensory organelles present on most cells in the body that participate broadly in human health, physiology and disease. Cilia generate a unique environment for signal transduction with tight control of protein, lipid and second messenger concentrations within a relatively small compartment, enabling reception, transmission and integration of biological information. In this Review, we discuss how cilia function as signalling hubs in cell-cell communication using three signalling pathways as examples: ciliary G-protein-coupled receptors (GPCRs), the Hedgehog (Hh) pathway and polycystin ion channels. We review how defects in these ciliary signalling pathways lead to a heterogeneous group of conditions known as 'ciliopathies', including metabolic syndromes, birth defects and polycystic kidney disease. Emerging understanding of these pathways' transduction mechanisms reveals common themes between these cilia-based signalling pathways that may apply to other pathways as well. These mechanistic insights reveal how cilia orchestrate normal and pathophysiological signalling outputs broadly throughout human biology.

GPR161 structure uncovers the redundant role of sterol-regulated ciliary cAMP signaling in the Hedgehog pathway.

The orphan G protein-coupled receptor (GPCR) GPR161 plays a central role in development by suppressing Hedgehog signaling. The fundamental basis of how GPR161 is activated remains unclear. Here, we determined a cryogenic-electron microscopy structure of active human GPR161 bound to heterotrimeric Gs. This structure revealed an extracellular loop 2 that occupies the canonical GPCR orthosteric ligand pocket. Furthermore, a sterol that binds adjacent to transmembrane helices 6 and 7 stabilizes a GPR161 conformation required for Gs coupling. Mutations that prevent sterol binding to GPR161 suppress Gs-mediated signaling. These mutants retain the ability to suppress GLI2 transcription factor accumulation in primary cilia, a key function of ciliary GPR161. By contrast, a protein kinase A-binding site in the GPR161 C terminus is critical in suppressing GLI2 ciliary accumulation. Our work highlights how structural features of GPR161 interface with the Hedgehog pathway and sets a foundation to understand the role of GPR161 function in other signaling pathways.

ADPKD-Causing Missense Variants in Polycystin-1 Disrupt Cell Surface Localization or Polycystin Channel Function.

Autosomal dominant polycystic kidney disease (ADPKD) is the leading monogenic cause of kidney failure and affects millions of people worldwide. Despite the prevalence of this monogenic disorder, our limited mechanistic understanding of ADPKD has hindered therapeutic development. Here, we successfully developed bioassays that functionally classify missense variants in polycystin-1 (PC1). Strikingly, ADPKD pathogenic missense variants cluster into two major categories: 1) those that disrupt polycystin cell surface localization or 2) those that attenuate polycystin ion channel activity. We found that polycystin channels with defective surface localization could be rescued with a small molecule. We propose that small-molecule-based strategies to improve polycystin cell surface localization and channel function will be effective therapies for ADPKD patients.

Primary cilia control oligodendrocyte precursor cell proliferation in white matter injury via Hedgehog-independent CREB signaling.

Remyelination after white matter injury (WMI) often fails in diseases such as multiple sclerosis because of improper recruitment and repopulation of oligodendrocyte precursor cells (OPCs) in lesions. How OPCs elicit specific intracellular programs in response to a chemically and mechanically diverse environment to properly regenerate myelin remains unclear. OPCs construct primary cilia, specialized signaling compartments that transduce Hedgehog (Hh) and G-protein-coupled receptor (GPCR) signals. We investigated the role of primary cilia in the OPC response to WMI. Removing cilia from OPCs genetically via deletion of Ift88 results in OPCs failing to repopulate WMI lesions because of reduced proliferation. Interestingly, loss of cilia does not affect Hh signaling in OPCs or their responsiveness to Hh signals but instead leads to dysfunctional cyclic AMP (cAMP)-dependent cAMP response element-binding protein (CREB)-mediated transcription. Because inhibition of CREB activity in OPCs reduces proliferation, we propose that a GPCR/cAMP/CREB signaling axis initiated at OPC cilia orchestrates OPC proliferation during development and in response to WMI.

Rab35 is required for embryonic development and kidney and ureter homeostasis through regulation of epithelial cell junctions.

Background: Rab35 is a member of a GTPase family of endocytic trafficking proteins. Studies in cell lines have indicated that Rab35 participates in cell adhesion, polarity, cytokinesis, and primary cilia length and composition. Additionally, sea urchin Rab35 regulates actin organization and is required for gastrulation. In mice, loss of Rab35 in the CNS disrupts hippocampal development and neuronal organization. Outside of the CNS, the functions of mammalian Rab35 in vivo are unknown. Methods: We generated and analyzed the consequences of both congenital and conditional null Rab35 mutations in mice. Using a LacZ reporter allele, we assessed Rab35 expression during development and postnatally. We assessed Rab35 loss in the kidney and ureter using histology, immunofluorescence microscopy, and western blotting. Results: Congenital Rab35 loss of function caused embryonic lethality: homozygous mutants arrested at E7.5 with cardiac edema. Conditional loss of Rab35, either during gestation or postnatally, caused hydronephrosis. The kidney and ureter phenotype were associated with disrupted actin cytoskeletal architecture, altered Arf6 epithelial polarity, reduced adherens junctions, loss of tight junction formation, defects in EGFR expression and localization, disrupted cell differentiation, and shortened primary cilia. Conclusion: Rab35 is essential for mammalian development and the maintenance of kidney and ureter architecture. Loss of Rab35 leads to non-obstructive hydronephrosis, making the Rab35 mutant mouse a novel mammalian model to study mechanisms underlying this disease.

Polycystin-2-dependent transcriptome reveals early response of autosomal dominant polycystic kidney disease.

Autosomal dominant polycystic kidney disease (ADPKD) is caused by mutations in polycystin genes, Pkd1 and Pkd2, but the underlying pathogenic mechanisms are poorly understood. To identify genes and pathways that operate downstream of polycystin-2 (PC2), a comprehensive gene expression database was created, cataloging changes in the transcriptome immediately following PC2 protein depletion. To explore cyst initiation processes, an immortalized mouse inner medullary collecting duct line was developed with the ability to knock out the Pkd2 gene conditionally. Genome-wide transcriptome profiling was performed using RNA sequencing in the cells immediately after PC2 was depleted and compared with isogenic control cells. Differentially expressed genes were identified, and a bioinformatic analysis pipeline was implemented. Altered expression of candidate cystogenic genes was validated in Pkd2 knockout mice. The expression of nearly 900 genes changed upon PC2 depletion. Differentially expressed genes were enriched for genes encoding components of the primary cilia, the canonical Wnt pathway, and MAPK signaling. Among the PC2-dependent ciliary genes, the transcription factor Glis3 was significantly downregulated. MAPK signaling formed a key node at the epicenter of PC2-dependent signaling networks. Activation of Wnt and MAPK signaling, concomitant with the downregulation of Glis3, was corroborated in Pkd2 knockout mice. The data identify a PC2 cilia-to-nucleus signaling axis and dysregulation of the Gli-similar subfamily of transcription factors as a potential initiator of cyst formation in ADPKD. The catalog of PC2-regulated genes should provide a valuable resource for future ADPKD research and new opportunities for drug development.NEW & NOTEWORTHY Autosomal dominant polycystic kidney disease (ADPKD) is the most common inherited kidney disease. Mutations in polycystin genes cause the disease, but the underlying mechanisms of cystogenesis are unknown. To help fill this knowledge gap, we created an inducible cell model of ADPKD and assembled a catalog of genes that respond in immediate proximity to polycystin-2 depletion using transcriptomic profiling. The catalog unveils a ciliary signaling-to-nucleus axis proximal to polycystin-2 dysfunction, highlighting Glis, Wnt, and MAPK signaling.

An RFX transcription factor regulates ciliogenesis in the closest living relatives of animals.

Cilia allowed our protistan ancestors to sense and explore their environment, avoid predation, and capture bacterial prey.1,2,3 Regulated ciliogenesis was likely critical for early animal evolution,2,4,5,6 and in modern animals, deploying cilia in the right cells at the right time is crucial for development and physiology. Two transcription factors, RFX and FoxJ1, coordinate ciliogenesis in animals7,8,9 but are absent from the genomes of many other ciliated eukaryotes, raising the question of how the regulation of ciliogenesis in animals evolved.10,11 By comparing the genomes of animals with those of their closest living relatives, the choanoflagellates, we found that the genome of their last common ancestor encoded at least three RFX paralogs and a FoxJ1 homolog. Disruption of the RFX homolog cRFXa in the model choanoflagellate Salpingoeca rosetta resulted in delayed cell proliferation and aberrant ciliogenesis, marked by the collapse and resorption of nascent cilia. In cRFXa mutants, ciliogenesis genes and foxJ1 were significantly downregulated. Moreover, the promoters of S. rosetta ciliary genes are enriched for DNA motifs matching those bound by the cRFXa protein in vitro. These findings suggest that an ancestral cRFXa homolog coordinated ciliogenesis in the progenitors of animals and choanoflagellates and that the selective deployment of the RFX regulatory module may have been necessary to differentiate ciliated from non-ciliated cell types during early animal evolution.

Seriously cilia: A tiny organelle illuminates evolution, disease, and intercellular communication.

The borders between cell and developmental biology, which have always been permeable, have largely dissolved. One manifestation is the blossoming of cilia biology, with cell and developmental approaches (increasingly complemented by human genetics, structural insights, and computational analysis) fruitfully advancing understanding of this fascinating, multifunctional organelle. The last eukaryotic common ancestor probably possessed a motile cilium, providing evolution with ample opportunity to adapt cilia to many jobs. Over the last decades, we have learned how non-motile, primary cilia play important roles in intercellular communication. Reflecting their diverse motility and signaling functions, compromised cilia cause a diverse range of diseases collectively called "ciliopathies." In this review, we highlight how cilia signal, focusing on how second messengers generated in cilia convey distinct information; how cilia are a potential source of signals to other cells; how evolution may have shaped ciliary function; and how cilia research may address thorny outstanding questions.

Physiological Condition-Dependent Changes in Ciliary GPCR Localization in the Brain.

Primary cilia are cellular appendages critical for diverse types of Signaling. They are found on most cell types, including cells throughout the CNS. Cilia preferentially localize certain G-protein-coupled receptors (GPCRs) and are critical for mediating the signaling of these receptors. Several of these neuronal GPCRs have recognized roles in feeding behavior and energy homeostasis. Cell and model systems, such as Caenorhabditis elegans and Chlamydomonas, have implicated both dynamic GPCR cilia localization and cilia length and shape changes as key for signaling. It is unclear whether mammalian ciliary GPCRs use similar mechanisms in vivo and under what conditions these processes may occur. Here, we assess two neuronal cilia GPCRs, melanin-concentrating hormone receptor 1 (MCHR1) and neuropeptide-Y receptor 2 (NPY2R), as mammalian model ciliary receptors in the mouse brain. We test the hypothesis that dynamic localization to cilia occurs under physiological conditions associated with these GPCR functions. Both receptors are involved in feeding behaviors, and MCHR1 is also associated with sleep and reward. Cilia were analyzed with a computer-assisted approach allowing for unbiased and high-throughput analysis. We measured cilia frequency, length, and receptor occupancy. We observed changes in ciliary length, receptor occupancy, and cilia frequency under different conditions for one receptor but not another and in specific brain regions. These data suggest that dynamic cilia localization of GPCRs depends on properties of individual receptors and cells where they are expressed. A better understanding of subcellular localization dynamics of ciliary GPCRs could reveal unknown molecular mechanisms regulating behaviors like feeding.

Centriolar satellites expedite mother centriole remodeling to promote ciliogenesis.

Centrosomes are orbited by centriolar satellites, dynamic multiprotein assemblies nucleated by Pericentriolar material 1 (PCM1). To study the requirement for centriolar satellites, we generated mice lacking PCM1, a crucial component of satellites. Pcm1-/- mice display partially penetrant perinatal lethality with survivors exhibiting hydrocephalus, oligospermia, and cerebellar hypoplasia, and variably expressive phenotypes such as hydronephrosis. As many of these phenotypes have been observed in human ciliopathies and satellites are implicated in cilia biology, we investigated whether cilia were affected. PCM1 was dispensable for ciliogenesis in many cell types, whereas Pcm1-/- multiciliated ependymal cells and human PCM1-/- retinal pigmented epithelial 1 (RPE1) cells showed reduced ciliogenesis. PCM1-/- RPE1 cells displayed reduced docking of the mother centriole to the ciliary vesicle and removal of CP110 and CEP97 from the distal mother centriole, indicating compromised early ciliogenesis. Similarly, Pcm1-/- ependymal cells exhibited reduced removal of CP110 from basal bodies in vivo. We propose that PCM1 and centriolar satellites facilitate efficient trafficking of proteins to and from centrioles, including the departure of CP110 and CEP97 to initiate ciliogenesis, and that the threshold to trigger ciliogenesis differs between cell types.

MRAP2 regulates energy homeostasis by promoting primary cilia localization of MC4R.

The G protein-coupled receptor melanocortin-4 receptor (MC4R) and its associated protein melanocortin receptor-associated protein 2 (MRAP2) are essential for the regulation of food intake and body weight in humans. MC4R localizes and functions at the neuronal primary cilium, a microtubule-based organelle that senses and relays extracellular signals. Here, we demonstrate that MRAP2 is critical for the weight-regulating function of MC4R neurons and the ciliary localization of MC4R. More generally, our study also reveals that GPCR localization to primary cilia can require specific accessory proteins that may not be present in heterologous cell culture systems. Our findings further demonstrate that targeting of MC4R to neuronal primary cilia is essential for the control of long-term energy homeostasis and suggest that genetic disruption of MC4R ciliary localization may frequently underlie inherited forms of obesity.

Zika virus alters centrosome organization to suppress the innate immune response.

Zika virus (ZIKV) is a flavivirus transmitted via mosquitoes and sex to cause congenital neurodevelopmental defects, including microcephaly. Inherited forms of microcephaly (MCPH) are associated with disrupted centrosome organization. Similarly, we found that ZIKV infection disrupted centrosome organization. ZIKV infection disrupted the organization of centrosomal proteins including CEP63, a MCPH-associated protein. The ZIKV nonstructural protein NS3 bound CEP63, and expression of NS3 was sufficient to alter centrosome architecture and CEP63 localization. Loss of CEP63 suppressed ZIKV-induced centrosome disorganization, indicating that ZIKV requires CEP63 to disrupt centrosome organization. ZIKV infection or CEP63 loss decreased the centrosomal localization and stability of TANK-binding kinase 1 (TBK1), a regulator of the innate immune response. ZIKV infection also increased the centrosomal accumulation of the CEP63 interactor DTX4, a ubiquitin ligase that degrades TBK1. Therefore, we propose that ZIKV disrupts CEP63 function to increase centrosomal DTX4 localization and destabilization of TBK1, thereby tempering the innate immune response.

The Intimate Connection Between Lipids and Hedgehog Signaling.

Hedgehog (HH) signaling is an intercellular communication pathway involved in directing the development and homeostasis of metazoans. HH signaling depends on lipids that covalently modify HH proteins and participate in signal transduction downstream. In many animals, the HH pathway requires the primary cilium, an organelle with a specialized protein and lipid composition. Here, we review the intimate connection between HH signaling and lipids. We highlight how lipids in the primary cilium can create a specialized microenvironment to facilitate signaling, and how HH and components of the HH signal transduction pathway use lipids to communicate between cells.

Endoderm development requires centrioles to restrain p53-mediated apoptosis in the absence of ERK activity.

Centrioles comprise the heart of centrosomes, microtubule-organizing centers. To study the function of centrioles in lung and gut development, we genetically disrupted centrioles throughout the mouse endoderm. Surprisingly, removing centrioles from the endoderm did not disrupt intestinal growth or development but blocked lung branching. In the lung, acentriolar SOX2-expressing airway epithelial cells apoptosed. Loss of centrioles activated p53, and removing p53 restored survival of SOX2-expressing cells, lung branching, and mouse viability. To investigate how endodermal p53 activation specifically killed acentriolar SOX2-expressing cells, we assessed ERK, a prosurvival cue. ERK was active throughout the intestine and in the distal lung buds, correlating with tolerance to centriole loss. Pharmacologically inhibiting ERK activated apoptosis in acentriolar cells, revealing that ERK activity protects acentriolar cells from apoptosis. Therefore, centrioles are largely dispensable for endodermal growth and the spatial distribution of ERK activity in the endoderm shapes the developmental consequences of centriolar defects and p53 activation.