RFX2 is broadly required for ciliogenesis during vertebrate development

CilioGenics: an integrated method and database for predicting novel ciliary genes

Uncovering the full list of human ciliary genes holds enormous promise for the diagnosis of cilia-related human diseases, collectively known as ciliopathies. Currently, genetic diagnoses of many ciliopathies remain incomplete (1-3). While various independent approaches theoretically have the potential to reveal the entire list of ciliary genes, approximately 30% of the genes on the ciliary gene list still stand as ciliary candidates (4,5). These methods, however, have mainly relied on a single strategy to uncover ciliary candidate genes, making the categorization challenging due to variations in quality and distinct capabilities demonstrated by different methodologies. Here, we develop a method called CilioGenics that combines several methodologies (single-cell RNA sequencing, protein-protein interactions (PPIs), comparative genomics, transcription factor (TF) network analysis, and text mining) to predict the ciliary capacity of each human gene. Our combined approach provides a CilioGenics score for every human gene that represents the probability that it will become a ciliary gene. Compared to methods that rely on a single method, CilioGenics performs better in its capacity to predict ciliary genes. Our top 500 gene list includes 258 new ciliary candidates, with 31 validated experimentally by us and others. Users may explore the whole list of human genes and CilioGenics scores on the CilioGenics database (https://ciliogenics.com/).

Primary cilia promote the differentiation of human neurons through the WNT signaling pathway

BACKGROUND

Primary cilia emanate from most human cell types, including neurons. Cilia are important for communicating with the cell's immediate environment: signal reception and transduction to/from the ciliated cell. Deregulation of ciliary signaling can lead to ciliopathies and certain neurodevelopmental disorders. In the developing brain cilia play well-documented roles for the expansion of the neural progenitor cell pool, while information about the roles of cilia during post-mitotic neuron differentiation and maturation is scarce.

RESULTS

We employed ciliated Lund Human Mesencephalic (LUHMES) cells in time course experiments to assess the impact of ciliary signaling on neuron differentiation. By comparing ciliated and non-ciliated neuronal precursor cells and neurons in wild type and in RFX2 -/- mutant neurons with altered cilia, we discovered an early-differentiation "ciliary time window" during which transient cilia promote axon outgrowth, branching and arborization. Experiments in neurons with IFT88 and IFT172 ciliary gene knockdowns, leading to shorter cilia, confirm these results. Cilia promote neuron differentiation by tipping WNT signaling toward the non-canonical pathway, in turn activating WNT pathway output genes implicated in cyto-architectural changes.

CONCLUSIONS

We provide a mechanistic entry point into when and how ciliary signaling coordinates, promotes and translates into anatomical changes. We hypothesize that ciliary alterations causing neuron differentiation defects may result in "mild" impairments of brain development, possibly underpinning certain aspects of neurodevelopmental disorders.

Genetic mechanisms of multiciliated cell development: from fate choice to differentiation in zebrafish and other models

Multiciliated cells (MCCS) form bundles of cilia and their activities are essential for the proper development and physiology of many organ systems. Not surprisingly, defects in MCCs have profound consequences and are associated with numerous disease states. Here, we discuss the current understanding of MCC formation, with a special focus on the genetic and molecular mechanisms of MCC fate choice and differentiation. Furthermore, we cast a spotlight on the use of zebrafish to study MCC ontogeny and several recent advances made in understanding MCCs using this vertebrate model to delineate mechanisms of MCC emergence in the developing kidney.

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.

Bidirectional multiciliated cell extrusion is controlled by Notch-driven basal extrusion and Piezo1-driven apical extrusion

Xenopus embryos are covered with a complex epithelium containing numerous multiciliated cells (MCCs). During late-stage development, there is a dramatic remodeling of the epithelium that involves the complete loss of MCCs. Cell extrusion is a well-characterized process for driving cell loss while maintaining epithelial barrier function. Normal cell extrusion is typically unidirectional, whereas bidirectional extrusion is often associated with disease (e.g. cancer). We describe two distinct mechanisms for MCC extrusion, a basal extrusion driven by Notch signaling and an apical extrusion driven by Piezo1. Early in the process there is a strong bias towards basal extrusion, but as development continues there is a shift towards apical extrusion. Importantly, response to the Notch signal is age dependent and governed by the maintenance of the MCC transcriptional program such that extension of this program is protective against cell loss. In contrast, later apical extrusion is regulated by Piezo1, such that premature activation of Piezo1 leads to early extrusion while blocking Piezo1 leads to MCC maintenance. Distinct mechanisms for MCC loss underlie the importance of their removal during epithelial remodeling.

MOF-mediated histone H4 Lysine 16 acetylation governs mitochondrial and ciliary functions by controlling gene promoters

Histone H4 lysine 16 acetylation (H4K16ac), governed by the histone acetyltransferase MOF, orchestrates gene expression regulation and chromatin interaction. However, the roles of MOF and H4K16ac in controlling cellular function and regulating mammalian tissue development remain unclear. Here we show that conditional deletion of Mof in the skin, but not Kansl1, causes severe defects in the self-renewal of basal epithelial progenitors, epidermal differentiation, and hair follicle growth, resulting in barrier defects and perinatal lethality. MOF-regulated genes are highly enriched for essential functions in the mitochondria and cilia. Genetic deletion of Uqcrq, an essential subunit for the electron transport chain (ETC) Complex III, in the skin, recapitulates the defects in epidermal differentiation and hair follicle growth observed in MOF knockout mouse. Together, this study reveals the requirement of MOF-mediated epigenetic mechanism for regulating mitochondrial and ciliary gene expression and underscores the important function of the MOF/ETC axis for mammalian skin development.

Forkhead transcription factor FKH-8 cooperates with RFX in the direct regulation of sensory cilia in Caenorhabditis elegans

Cilia, either motile or non-motile (a.k.a primary or sensory), are complex evolutionarily conserved eukaryotic structures composed of hundreds of proteins required for their assembly, structure and function that are collectively known as the ciliome. Ciliome gene mutations underlie a group of pleiotropic genetic diseases known as ciliopathies. Proper cilium function requires the tight coregulation of ciliome gene transcription, which is only fragmentarily understood. RFX transcription factors (TF) have an evolutionarily conserved role in the direct activation of ciliome genes both in motile and non-motile cilia cell-types. In vertebrates, FoxJ1 and FoxN4 Forkhead (FKH) TFs work with RFX in the direct activation of ciliome genes, exclusively in motile cilia cell-types. No additional TFs have been described to act together with RFX in primary cilia cell-types in any organism. Here we describe FKH-8, a FKH TF, as a direct regulator of the sensory ciliome genes in Caenorhabditis elegans. FKH-8 is expressed in all ciliated neurons in C. elegans, binds the regulatory regions of ciliome genes, regulates ciliome gene expression, cilium morphology and a wide range of behaviors mediated by sensory ciliated neurons. FKH-8 and DAF-19 (C. elegans RFX) physically interact and synergistically regulate ciliome gene expression. C. elegans FKH-8 function can be replaced by mouse FOXJ1 and FOXN4 but not by other members of other mouse FKH subfamilies. In conclusion, RFX and FKH TF families act jointly as direct regulators of ciliome genes also in sensory ciliated cell types suggesting that this regulatory logic could be an ancient trait predating functional cilia sub-specialization.

Persistence of FoxJ1+ Pax6+ Sox2+ ependymal cells throughout life in the human spinal cord

Ependymal cells lining the central canal of the spinal cord play a crucial role in providing a physical barrier and in the circulation of cerebrospinal fluid. These cells express the FOXJ1 and SOX2 transcription factors in mice and are derived from various neural tube populations, including embryonic roof and floor plate cells. They exhibit a dorsal-ventral expression pattern of spinal cord developmental transcription factors (such as MSX1, PAX6, ARX, and FOXA2), resembling an embryonic-like organization. Although this ependymal region is present in young humans, it appears to be lost with age. To re-examine this issue, we collected 17 fresh spinal cords from organ donors aged 37-83 years and performed immunohistochemistry on lightly fixed tissues. We observed cells expressing FOXJ1 in the central region in all cases, which co-expressed SOX2 and PAX6 as well as RFX2 and ARL13B, two proteins involved in ciliogenesis and cilia-mediated sonic hedgehog signaling, respectively. Half of the cases exhibited a lumen and some presented portions of the spinal cord with closed and open central canals. Co-staining of FOXJ1 with other neurodevelopmental transcription factors (ARX, FOXA2, MSX1) and NESTIN revealed heterogeneity of the ependymal cells. Interestingly, three donors aged > 75 years exhibited a fetal-like regionalization of neurodevelopmental transcription factors, with dorsal and ventral ependymal cells expressing MSX1, ARX, and FOXA2. These results provide new evidence for the persistence of ependymal cells expressing neurodevelopmental genes throughout human life and highlight the importance of further investigation of these cells.

GEMC1 and MCIDAS interactions with SWI/SNF complexes regulate the multiciliated cell-specific transcriptional program

Multiciliated cells (MCCs) project dozens to hundreds of motile cilia from their apical surface to promote the movement of fluids or gametes in the mammalian brain, airway or reproductive organs. Differentiation of MCCs requires the sequential action of the Geminin family transcriptional activators, GEMC1 and MCIDAS, that both interact with E2F4/5-DP1. How these factors activate transcription and the extent to which they play redundant functions remains poorly understood. Here, we demonstrate that the transcriptional targets and proximal proteomes of GEMC1 and MCIDAS are highly similar. However, we identified distinct interactions with SWI/SNF subcomplexes; GEMC1 interacts primarily with the ARID1A containing BAF complex while MCIDAS interacts primarily with BRD9 containing ncBAF complexes. Treatment with a BRD9 inhibitor impaired MCIDAS-mediated activation of several target genes and compromised the MCC differentiation program in multiple cell based models. Our data suggest that the differential engagement of distinct SWI/SNF subcomplexes by GEMC1 and MCIDAS is required for MCC-specific transcriptional regulation and mediated by their distinct C-terminal domains.

Bidirectional multiciliated cell extrusion is controlled by Notch driven basal extrusion and Piezo 1 driven apical extrusion

Xenopus embryos are covered with a complex epithelium containing numerous multiciliated cells (MCCs). During late stage development there is a dramatic remodeling of the epithelium that involves the complete loss of MCCs. Cell extrusion is a well-characterized process for driving cell loss while maintaining epithelial barrier function. Normal cell extrusion is typically unidirectional whereas bidirectional extrusion is often associated with disease (e.g. cancer). We describe two distinct mechanisms for MCC extrusion, a basal extrusion driven by Notch signaling and an apical extrusion driven by Piezo1. Early in the process there is a strong bias towards basal extrusion, but as development continues there is a shift towards apical extrusion. Importantly, receptivity to the Notch signal is age-dependent and governed by the maintenance of the MCC transcriptional program such that extension of this program is protective against cell loss. In contrast, later apical extrusion is regulated by Piezo 1 such that premature activation of Piezo 1 leads to early extrusion while blocking Piezo 1 leads to MCC maintenance. Distinct mechansms for MCC loss underlie the importance of their removal during epithelial remodeling.

Summay Statement

Cell extrusion typically occurs unidirectionally. We have identified a single population of multiciliated cells that extrudes bidirectionally: Notch-driven basal extrusion and Piezo 1-mediated apical extrusion.

PCD Genes-From Patients to Model Organisms and Back to Humans

Primary ciliary dyskinesia (PCD) is a hereditary genetic disorder caused by the lack of motile cilia or the assembxly of dysfunctional ones. This rare human disease affects 1 out of 10,000-20,000 individuals and is caused by mutations in at least 50 genes. The past twenty years brought significant progress in the identification of PCD-causative genes and in our understanding of the connections between causative mutations and ciliary defects observed in affected individuals. These scientific advances have been achieved, among others, due to the extensive motile cilia-related research conducted using several model organisms, ranging from protists to mammals. These are unicellular organisms such as the green alga Chlamydomonas, the parasitic protist Trypanosoma, and free-living ciliates, Tetrahymena and Paramecium, the invertebrate Schmidtea, and vertebrates such as zebrafish, Xenopus, and mouse. Establishing such evolutionarily distant experimental models with different levels of cell or body complexity was possible because both basic motile cilia ultrastructure and protein composition are highly conserved throughout evolution. Here, we characterize model organisms commonly used to study PCD-related genes, highlight their pros and cons, and summarize experimental data collected using these models.

The First Report of a Missense Variant in RFX2 Causing Non-Syndromic Tooth Agenesis in a Consanguineous Pakistani Family

Background: The syndromic and non-syndromic congenital missing teeth phenotype is termed tooth agenesis. Since tooth agenesis is a heterogeneous disorder hence, the patients show diverse absent teeth phenotypes. Thus identifying novel genes involved in the morphogenesis of ectodermal appendages, including teeth, paves the way for establishing signaling pathways.

Methods and Results: We have recruited an autosomal recessive non-syndromic tooth agenesis family with two affected members. The exome sequencing technology identified a novel missense sequence variant c.1421T > C; p.(Ile474Thr) in a regulatory factor X (RFX) family member (RFX2, OMIM: 142,765). During the data analysis eight rare variants on various chromosomal locations were identified, but the co-segregation analysis using Sanger sequencing confirmed the segregation of only two variants RFX2: c.1421T > C; p.(Ile474Thr), DOHH: c.109C > G; p.(Pro37Ala) lying in a common 7.1 MB region of homozygosity on chromosome 19p13.3. Furthermore, the online protein prediction algorithms and protein modeling analysis verified the RFX2 variant as a damaging genetic alteration and ACMG pathogenicity criteria classified it as likely pathogenic. On the other hand, the DOHH variant showed benign outcomes.

Conclusion:RFX2 regulates the Hedgehog and fibroblast growth factor signaling pathways, which are involved in the epithelial and mesenchymal interactions during tooth development. Prior animal model studies have confirmed the expression of rfx2 at a developmental stage governing mouth formation. Moreover, its regulatory role and close association with ciliary and non-ciliary genes causing various dental malformations makes it a potential candidate gene for tooth agenesis phenotype. Further studies will contribute to exploring the direct role of RFX2 in human tooth development.

Cell-of-Origin and Genetic, Epigenetic, and Microenvironmental Factors Contribute to the Intra-Tumoral Heterogeneity of Pediatric Intracranial Ependymoma

Intra-tumoral heterogeneity (ITH) is a complex multifaceted phenomenon that posits major challenges for the clinical management of cancer patients. Genetic, epigenetic, and microenvironmental factors are concurrent drivers of diversity among the distinct populations of cancer cells. ITH may also be installed by cancer stem cells (CSCs), that foster unidirectional hierarchy of cellular phenotypes or, alternatively, shift dynamically between distinct cellular states. Ependymoma (EPN), a molecularly heterogeneous group of tumors, shows a specific spatiotemporal distribution that suggests a link between ependymomagenesis and alterations of the biological processes involved in embryonic brain development. In children, EPN most often arises intra-cranially and is associated with an adverse outcome. Emerging evidence shows that EPN displays large intra-patient heterogeneity. In this review, after touching on EPN inter-tumoral heterogeneity, we focus on the sources of ITH in pediatric intra-cranial EPN in the framework of the CSC paradigm. We also examine how single-cell technology has shed new light on the complexity and developmental origins of EPN and the potential impact that this understanding may have on the therapeutic strategies against this deadly pediatric malignancy.

Novel variations in spermatogenic transcription regulators RFX2 and TAF7 increase risk of azoospermia

PURPOSE

Genetic etiology of idiopathic male infertility is enigmatic owing to involvement of multiple gene regulatory networks in spermatogenesis process. Any change in optimal function of the transcription factors involved in this process owing to polymorphisms/mutations may increase the risk of infertility. We investigated polymorphisms/mutations of spermatogenic transcription regulators TAF7 and RFX2 and analysed their association with incidence of azoospermia among the men from West Bengal, India.

METHODS

Genotyping was carried by Sanger's dideoxy sequencing of 130 azoospermic men who were detected negative in Y chromosome microdeletion screening and 140 healthy controls. Association study was done by suitable statistical methods. In silico analysis was performed to infer the intuitive damaging effects of detected variants at transcripts and protein level.

RESULTS

We found significant association of TAF7 C16T (MW827584 G > A), RFX2 562delT (MZ560629delA), rs11547633 A > C, rs17606721 A > G, MW827583 C > T, and MZ379836 C > T variants with the incidence of azoospermia. In silico analysis predicted that the variants either alter the natural splice junctions of the transcript or cause probable damage in the structure of proteins of respective genes.

CONCLUSION

Polymorphisms/mutations of TAF7 and RFX2 genes increase risk of male infertility in Bengali population. The novel variants may be used as markers for male infertility screening in ART practise.

Identification and functional modelling of plausibly causative cis-regulatory variants in a highly-selected cohort with X-linked intellectual disability

Identifying causative variants in cis-regulatory elements (CRE) in neurodevelopmental disorders has proven challenging. We have used in vivo functional analyses to categorize rigorously filtered CRE variants in a clinical cohort that is plausibly enriched for causative CRE mutations: 48 unrelated males with a family history consistent with X-linked intellectual disability (XLID) in whom no detectable cause could be identified in the coding regions of the X chromosome (chrX). Targeted sequencing of all chrX CRE identified six rare variants in five affected individuals that altered conserved bases in CRE targeting known XLID genes and segregated appropriately in families. Two of these variants, FMR1CRE and TENM1CRE, showed consistent site- and stage-specific differences of enhancer function in the developing zebrafish brain using dual-color fluorescent reporter assay. Mouse models were created for both variants. In male mice Fmr1CRE induced alterations in neurodevelopmental Fmr1 expression, olfactory behavior and neurophysiological indicators of FMRP function. The absence of another likely causative variant on whole genome sequencing further supported FMR1CRE as the likely basis of the XLID in this family. Tenm1CRE mice showed no phenotypic anomalies. Following the release of gnomAD 2.1, reanalysis showed that TENM1CRE exceeded the maximum plausible population frequency of a XLID causative allele. Assigning causative status to any ultra-rare CRE variant remains problematic and requires disease-relevant in vivo functional data from multiple sources. The sequential and bespoke nature of such analyses renders them time-consuming and challenging to scale for routine clinical use.

Single-Cell RNA Sequencing of Childhood Ependymoma Reveals Neoplastic Cell Subpopulations That Impact Molecular Classification and Etiology

Ependymoma (EPN) is a brain tumor commonly presenting in childhood that remains fatal in most children. Intra-tumoral cellular heterogeneity in bulk-tumor samples significantly confounds our understanding of EPN biology, impeding development of effective therapy. We, therefore, use single-cell RNA sequencing, histology, and deconvolution to catalog cellular heterogeneity of the major childhood EPN subgroups. Analysis of PFA subgroup EPN reveals evidence of an undifferentiated progenitor subpopulation that either differentiates into subpopulations with ependymal cell characteristics or transitions into a mesenchymal subpopulation. Histological analysis reveals that progenitor and mesenchymal subpopulations co-localize in peri-necrotic zones. In conflict with current classification paradigms, relative PFA subpopulation proportions are shown to determine bulk-tumor-assigned subgroups. We provide an interactive online resource that facilitates exploration of the EPN single-cell dataset. This atlas of EPN cellular heterogeneity increases understanding of EPN biology.

Control of multiciliogenesis by miR-34/449 in the male reproductive tract through enforcing cell cycle exit

Multiciliated cells (MCCs) are terminally differentiated postmitotic cells that possess hundreds of motile cilia on their apical surface. Defects in cilia formation are associated with ciliopathies that affect many organs. In this study, we tested the role and mechanism of the miR-34/449 family in the regulation of multiciliogenesis in EDs using an miR-34b/c−/−; miR-449−/− double knockout (dKO) mouse model. MiR-34b/c and miR-449 depletion led to a reduced number of MCCs and abnormal cilia structure in the EDs starting from postnatal day (P)14. However, abnormal MCC differentiation in the dKO EDs could be observed as early as P7. RNA-seq analyses revealed that the aberrant development of MCCs in the EDs of dKO mice was associated with the upregulation of genes involved in cell cycle control. Using a cyclin-dependent kinase inhibitor to force cell cycle exit promoted MCC differentiation, and partially rescued the defective multiciliogenesis in the EDs of dKO mice. Taken together, our results suggest that miR-34b/c and miR-449 play an essential role in multiciliogenesis in EDs by regulating cell cycle exit.

Summary: Mutagenic, expression and histological analyses reveal an essential role for miR-34b/c and miR-449 in multiciliogenesis in efferent ductules of the male reproductive tract by regulating cell cycle exit.

RFX1: a promising therapeutic arsenal against cancer

Regulatory factor X1 (RFX1) is an evolutionary conserved transcriptional factor that influences a wide range of cellular processes such as cell cycle, cell proliferation, differentiation, and apoptosis, by regulating a number of target genes that are involved in such processes. On a closer look, these target genes also play a key role in tumorigenesis and associated events. Such observations paved the way for further studies evaluating the role of RFX1 in cancer. These studies were indispensable due to the failure of conventional chemotherapeutic drugs to target key cellular hallmarks such as cancer stemness, cellular plasticity, enhanced drug efflux, de-regulated DNA repair machinery, and altered pathways evading apoptosis. In this review, we compile significant evidence for the tumor-suppressive activities of RFX1 while also analyzing its oncogenic potential in some cancers. RFX1 induction decreased cellular proliferation, modulated the immune system, induced apoptosis, reduced chemoresistance, and sensitized cancer stem cells for chemotherapy. Thus, our review discusses the pleiotropic function of RFX1 in multitudinous gene regulations, decisive protein–protein interactions, and also its role in regulating key cell signaling events in cancer. Elucidation of these regulatory mechanisms can be further utilized for RFX1 targeted therapy.

Xenopus to the rescue: A model to validate and characterize candidate ciliopathy genes

Cilia are present on most vertebrate cells and play a central role in development, growth, and homeostasis. Thus, cilia dysfunction can manifest into an array of diseases, collectively termed ciliopathies, affecting millions of lives worldwide. Yet, our understanding of the gene regulatory networks that control cilia assembly and functions remain incomplete. With the advances in next-generation sequencing technologies, we can now rapidly predict pathogenic variants from hundreds of ciliopathy patients. While the pace of candidate gene discovery is exciting, most of these genes have never been previously implicated in cilia assembly or function. This makes assigning the disease causality difficult. This review discusses how Xenopus, a genetically tractable and high-throughput vertebrate model, has played a central role in identifying, validating, and characterizing candidate ciliopathy genes. The review is focused on multiciliated cells (MCCs) and diseases associated with MCC dysfunction. MCCs harbor multiple motile cilia on their apical surface to generate extracellular fluid flow inside the airway, the brain ventricles, and the oviduct. In Xenopus, these cells are external and present on the embryonic epidermal epithelia, facilitating candidate genes analysis in MCC development in vivo. The ability to introduce patient variants to study their effects on disease progression makes Xenopus a powerful model to improve our understanding of the underlying disease mechanisms and explain the patient phenotype.

Motile cilia and airway disease

A finely regulated system of airway epithelial development governs the differentiation of motile ciliated cells of the human respiratory tract, conferring the body's mucociliary clearance defence system. Human cilia dysfunction can arise through genetic mutations and this is a cause of debilitating disease morbidities that confer a greatly reduced quality of life. The inherited human motile ciliopathy disorder, primary ciliary dyskinesia (PCD), can arise from mutations in genes affecting various aspects of motile cilia structure and function through deficient production, transport and assembly of cilia motility components or through defective multiciliogenesis. Our understanding about the development of the respiratory epithelium, motile cilia biology and the implications for human pathology has expanded greatly over the past 20 years since isolation of the first PCD gene, rising to now nearly 50 genes. Systems level insights about cilia motility in health and disease have been made possible through intensive molecular and omics (genomics, transcriptomics, proteomics) research, applied in ciliate organisms and in animal and human disease modelling. Here, we review ciliated airway development and the genetic stratification that underlies PCD, for which the underlying genotype can increasingly be connected to biological mechanism and disease prognostics. Progress in this field can facilitate clinical translation of research advances, with potential for great medical impact, e.g. through improvements in ciliopathy disease diagnosis, management, family counselling and by enhancing the potential for future genetically tailored approaches to disease therapeutics.

Interplay of RFX transcription factors 1, 2 and 3 in motile ciliogenesis

Cilia assembly is under strict transcriptional control during animal development. In vertebrates, a hierarchy of transcription factors (TFs) are involved in controlling the specification, differentiation and function of multiciliated epithelia. RFX TFs play key functions in the control of ciliogenesis in animals. Whereas only one RFX factor regulates ciliogenesis in C. elegans, several distinct RFX factors have been implicated in this process in vertebrates. However, a clear understanding of the specific and redundant functions of different RFX factors in ciliated cells remains lacking. Using RNA-seq and ChIP-seq approaches we identified genes regulated directly and indirectly by RFX1, RFX2 and RFX3 in mouse ependymal cells. We show that these three TFs have both redundant and specific functions in ependymal cells. Whereas RFX1, RFX2 and RFX3 occupy many shared genomic loci, only RFX2 and RFX3 play a prominent and redundant function in the control of motile ciliogenesis in mice. Our results provide a valuable list of candidate ciliary genes. They also reveal stunning differences between compensatory processes operating in vivo and ex vivo.

Genome-Wide Analysis Reveals Mucociliary Remodeling of the Nasal Airway Epithelium Induced by Urban PM2.5

Air pollution particulate matter <2.5 μm (PM_2.5) exposure is associated with poor respiratory outcomes. Mechanisms underlying PM_2.5-induced lung pathobiology are poorly understood but likely involve cellular and molecular changes to the airway epithelium. We extracted and chemically characterized the organic and water-soluble components of air pollution PM_2.5 samples, then determined the whole transcriptome response of human nasal mucociliary airway epithelial cultures to a dose series of PM_2.5 extracts. We found that PM_2.5 organic extract (OE), but not water-soluble extract, elicited a potent, dose-dependent transcriptomic response from the mucociliary epithelium. Exposure to a moderate OE dose modified the expression of 424 genes, including activation of aryl hydrocarbon receptor signaling and an IL-1 inflammatory program. We generated an OE-response gene network defined by eight functional enrichment groups, which exhibited high connectivity through CYP1A1, IL1A, and IL1B. This OE exposure also robustly activated a mucus secretory expression program (>100 genes), which included transcriptional drivers of mucus metaplasia (SPDEF and FOXA3). Exposure to a higher OE dose modified the expression of 1,240 genes and further exacerbated expression responses observed at the moderate dose, including the mucus secretory program. Moreover, the higher OE dose significantly increased the MUC5AC/MUC5B gel-forming mucin expression ratio and strongly downregulated ciliated cell expression programs, including key ciliating cell transcription factors (e.g., FOXJ1 and MCIDAS). Chronic OE stimulation induced mucus metaplasia-like remodeling characterized by increases in MUC5AC⁺ secretory cells and MUC5AC mucus secretions. This epithelial remodeling may underlie poor respiratory outcomes associated with high PM_2.5 exposure.

UTX/KDM6A suppresses AP-1 and a gliogenesis program during neural differentiation of human pluripotent stem cells

BACKGROUND

UTX/KDM6A is known to interact and influence multiple different chromatin modifiers to promote an open chromatin environment to facilitate gene activation, but its molecular activities in developmental gene regulation remain unclear.

RESULTS

We report that in human neural stem cells, UTX binding correlates with both promotion and suppression of gene expression. These activities enable UTX to modulate neural stem cell self-renewal, promote neurogenesis, and suppress gliogenesis. In neural stem cells, UTX has a less influence over histone H3 lysine 27 and lysine 4 methylation but more predominantly affects histone H3 lysine 27 acetylation and chromatin accessibility. Furthermore, UTX suppresses components of AP-1 and, in turn, a gliogenesis program.

CONCLUSIONS

Our findings revealed that UTX coordinates dualistic gene regulation to govern neural stem cell properties and neurogenesis-gliogenesis switch.

Building a ciliated epithelium: Transcriptional regulation and radial intercalation of multiciliated cells

The epidermis of the Xenopus embryo has emerged as a powerful tool for studying the development of a ciliated epithelium. Interspersed throughout the epithelium are multiciliated cells (MCCs) with 100+ motile cilia that beat in a coordinated manner to generate fluid flow over the surface of the cell. MCCs are essential for various developmental processes and, furthermore, ciliary dysfunction is associated with numerous pathologies. Therefore, understanding the cellular mechanisms involved in establishing a ciliated epithelium are of particular interest. MCCs originate in the inner epithelial layer of Xenopus skin, where Notch signaling plays a critical role in determining which progenitors will adopt a ciliated cell fate. Then, activation of various transcriptional regulators, such as GemC1 and MCIDAS, initiate the MCC transcriptional program, resulting in centriole amplification and the formation of motile cilia. Following specification and differentiation, MCCs undergo the process of radial intercalation, where cells apically migrate from the inner layer to the outer epithelial layer. This process involves the cooperation of various cytoskeletal networks, activation of various signaling molecules, and changes in cell-ECM and cell-cell adhesion. Coordination of these cellular processes is required for complete incorporation into the outer epithelial layer and generation of a functional ciliated epithelium. Here, we highlight recent advances made in understanding the transcriptional cascades required for MCC specification and differentiation and the coordination of cellular processes that facilitate radial intercalation. Proper regulation of these signaling pathways and processes are the foundation for developing a ciliated epithelium.

Dissecting the cellular specificity of smoking effects and reconstructing lineages in the human airway epithelium

Cigarette smoke first interacts with the lung through the cellularly diverse airway epithelium and goes on to drive development of most chronic lung diseases. Here, through single cell RNA-sequencing analysis of the tracheal epithelium from smokers and non-smokers, we generate a comprehensive atlas of epithelial cell types and states, connect these into lineages, and define cell-specific responses to smoking. Our analysis infers multi-state lineages that develop into surface mucus secretory and ciliated cells and then contrasts these to the unique specification of submucosal gland (SMG) cells. Accompanying knockout studies reveal that tuft-like cells are the likely progenitor of both pulmonary neuroendocrine cells and CFTR-rich ionocytes. Our smoking analysis finds that all cell types, including protected stem and SMG populations, are affected by smoking through both pan-epithelial smoking response networks and hundreds of cell-specific response genes, redefining the penetrance and cellular specificity of smoking effects on the human airway epithelium.

Signaling in the primary cilium through the lens of the Hedgehog pathway

Cilia are microtubule-based, cell-surface projections whose machinery is evolutionarily conserved. In vertebrates, cilia are observed on almost every cell type and are either motile or immotile. Immotile sensory, or primary cilia, are responsive to extracellular ligands and signals. Cilia can be thought of as compartments, functionally distinct from the cell that provides an environment for signaling cascades. Hedgehog is a critical developmental signaling pathway which is functionally linked to primary cilia in vertebrates. The major components of the vertebrate Hedgehog signaling pathway dynamically localize to the ciliary compartment and ciliary membrane. Critically, G-protein coupled receptor (GPCR) Smoothened, the obligate transducer of the pathway, is enriched and activated in the cilium. While Smoothened is the most intensely studied ciliary receptor, many GPCRs localize within cilia. Understanding the link between Smoothened and cilia defines common features, and distinctions, of GPCR signaling within the primary cilium. This article is categorized under: Signaling Pathways > Global Signaling Mechanisms Signaling Pathways > Cell Fate Signaling.

JMJD3 and UTX determine fidelity and lineage specification of human neural progenitor cells

Neurogenesis, a highly orchestrated process, entails the transition from a pluripotent to neural state and involves neural progenitor cells (NPCs) and neuronal/glial subtypes. However, the precise epigenetic mechanisms underlying fate decision remain poorly understood. Here, we delete KDM6s (JMJD3 and/or UTX), the H3K27me3 demethylases, in human embryonic stem cells (hESCs) and show that their deletion does not impede NPC generation from hESCs. However, KDM6-deficient NPCs exhibit poor proliferation and a failure to differentiate into neurons and glia. Mechanistically, both JMJD3 and UTX are found to be enriched in gene loci essential for neural development in hNPCs, and KDM6 impairment leads to H3K27me3 accumulation and blockade of DNA accessibility at these genes. Interestingly, forced expression of neuron-specific chromatin remodelling BAF (nBAF) rescues the neuron/glia defect in KDM6-deficient NPCs despite H3K27me3 accumulation. Our findings uncover the differential requirement of KDM6s in specifying NPCs and neurons/glia and highlight the contribution of individual epigenetic regulators in fate decisions in a human development model.

Neurogenesis is an ordered transition from pluriptotent cells to neural precursor cells (NPCs) to neurons. Here the authors show that loss of the lysine demethylases JMJD3 and UTX leads reduced DNA accessibility at neurogenesis loci in human NPCs, and that the chromatin remodeller BAF can rescue differentiation defects.

Rare Human Diseases: Model Organisms in Deciphering the Molecular Basis of Primary Ciliary Dyskinesia

Primary ciliary dyskinesia (PCD) is a recessive heterogeneous disorder of motile cilia, affecting one per 15,000-30,000 individuals; however, the frequency of this disorder is likely underestimated. Even though more than 40 genes are currently associated with PCD, in the case of approximately 30% of patients, the genetic cause of the manifested PCD symptoms remains unknown. Because motile cilia are highly evolutionarily conserved organelles at both the proteomic and ultrastructural levels, analyses in the unicellular and multicellular model organisms can help not only to identify new proteins essential for cilia motility (and thus identify new putative PCD-causative genes), but also to elucidate the function of the proteins encoded by known PCD-causative genes. Consequently, studies involving model organisms can help us to understand the molecular mechanism(s) behind the phenotypic changes observed in the motile cilia of PCD affected patients. Here, we summarize the current state of the art in the genetics and biology of PCD and emphasize the impact of the studies conducted using model organisms on existing knowledge.

Characterizing smoking-induced transcriptional heterogeneity in the human bronchial epithelium at single-cell resolution

The human bronchial epithelium is composed of multiple distinct cell types that cooperate to defend against environmental insults. While studies have shown that smoking alters bronchial epithelial function and morphology, its precise effects on specific cell types and overall tissue composition are unclear. We used single-cell RNA sequencing to profile bronchial epithelial cells from six never and six current smokers. Unsupervised analyses led to the characterization of a set of toxin metabolism genes that localized to smoker ciliated cells, tissue remodeling associated with a loss of club cells and extensive goblet cell hyperplasia, and a previously unidentified peri-goblet epithelial subpopulation in smokers who expressed a marker of bronchial premalignant lesions. Our data demonstrate that smoke exposure drives a complex landscape of cellular alterations that may prime the human bronchial epithelium for disease.

Histone H2B monoubiquitination regulates heart development via epigenetic control of cilia motility

Genomic analyses of patients with congenital heart disease (CHD) have identified significant contribution from mutations affecting cilia genes and chromatin remodeling genes; however, the mechanism(s) connecting chromatin remodeling to CHD is unknown. Histone H2B monoubiquitination (H2Bub1) is catalyzed by the RNF20 complex consisting of RNF20, RNF40, and UBE2B. Here, we show significant enrichment of loss-of-function mutations affecting H2Bub1 in CHD patients (enrichment 6.01, P = 1.67 × 10^-03), some of whom had abnormal laterality associated with ciliary dysfunction. In Xenopus, knockdown of rnf20 and rnf40 results in abnormal heart looping, defective development of left-right (LR) asymmetry, and impaired cilia motility. Rnf20, Rnf40, and Ube2b affect LR patterning and cilia synergistically. Examination of global H2Bub1 level in Xenopus embryos shows that H2Bub1 is developmentally regulated and requires Rnf20. To examine gene-specific H2Bub1, we performed ChIP-seq of mouse ciliated and nonciliated tissues and showed tissue-specific H2Bub1 marks significantly enriched at cilia genes including the transcription factor Rfx3 Rnf20 knockdown results in decreased levels of rfx3 mRNA in Xenopus, and exogenous rfx3 can rescue the Rnf20 depletion phenotype. These data suggest that Rnf20 functions at the Rfx3 locus regulating cilia motility and cardiac situs and identify H2Bub1 as an upstream transcriptional regulator controlling tissue-specific expression of cilia genes. Our findings mechanistically link the two functional gene ontologies that have been implicated in human CHD: chromatin remodeling and cilia function.

Desmoplakin is required for epidermal integrity and morphogenesis in the Xenopus laevis embryo

Desmoplakin (Dsp) is a unique and critical desmosomal protein, that is integral to epidermal development. However, it is unclear whether this protein is required specifically for epidermal morphogenesis. Using morpholinos or Crispr/Cas9 mutagenesis we decreased the function of Dsp in frog embryos to better understand its role during epidermal development. Dsp morphant and mutant embryos had developmental defects such as epidermal fragility that mimicked what has been reported in mammals. Most importantly, we also uncovered a novel function for Dsp in the morphogenesis of the epidermis in X. laevis. In particular, Dsp is required during the process of radial intercalation where basally located cells move into the outer epidermal layer. Once inserted these newly intercalated cells expand their apical surface and then they differentiate into specific epidermal cell types. Decreased levels of Dsp resulted in the failure of the radially intercalating cells to expand their apical surface, thereby reducing the number of differentiated multiciliated and secretory cells. Such defects correlate with changes in E-cadherin levels and actin and microtubule localization which could explain the defects in apical expansion. A mutated form of Dsp that maintains cell-cell adhesion but eliminates the connections to the cytoskeleton results in the same epidermal morphogenesis defect. These results suggest a specific role for Dsp in the apical expansion of cells during radial intercalation. We have developed a novel system, in the frog, to demonstrate for the first time that desmosomes not only protect against mechanical stress but are also critical for epidermal morphogenesis.

Iroquois transcription factor irx2a is required for multiciliated and transporter cell fate decisions during zebrafish pronephros development

The genetic regulation of nephron patterning during kidney organogenesis remains poorly understood. Nephron tubules in zebrafish are composed of segment populations that have unique absorptive and secretory roles, as well as multiciliated cells (MCCs) that govern fluid flow. Here, we report that the transcription factor iroquois 2a (irx2a) is requisite for zebrafish nephrogenesis. irx2a transcripts localized to the developing pronephros and maturing MCCs, and loss of function altered formation of two segment populations and reduced MCC number. Interestingly, irx2a deficient embryos had reduced expression of an essential MCC gene ets variant 5a (etv5a), and were rescued by etv5a overexpression, supporting the conclusion that etv5a acts downstream of irx2a to control MCC ontogeny. Finally, we found that retinoic acid (RA) signaling affects the irx2a expression domain in renal progenitors, positioning irx2a downstream of RA. In sum, this work reveals new roles for irx2a during nephrogenesis, identifying irx2a as a crucial connection between RA signaling, segmentation, and the control of etv5a mediated MCC formation. Further investigation of the genetic players involved in these events will enhance our understanding of the molecular pathways that govern renal development, which can be used help create therapeutics to treat congenital and acquired kidney diseases.

Defects in efferent duct multiciliogenesis underlie male infertility in GEMC1-, MCIDAS- or CCNO-deficient mice

GEMC1 and MCIDAS are geminin family proteins that transcriptionally activate E2F4/5-target genes during multiciliogenesis, including Foxj 1 and Ccno Male mice that lacked Gemc1, Mcidas or Ccno were found to be infertile, but the origin of this defect has remained unclear. Here, we show that all three genes are necessary for the generation of functional multiciliated cells in the efferent ducts that are required for spermatozoa to enter the epididymis. In mice that are mutant for Gemc1, Mcidas or Ccno, we observed a similar spectrum of phenotypes, including thinning of the seminiferous tubule epithelia, dilation of the rete testes, sperm agglutinations in the efferent ducts and lack of spermatozoa in the epididymis (azoospermia). These data suggest that defective efferent duct development is the dominant cause of male infertility in these mouse models, and this likely extends to individuals with the ciliopathy reduced generation of multiple motile cilia with mutations in MCIDAS and CCNO.

WDR5 regulates left-right patterning via chromatin-dependent and -independent functions

Congenital heart disease (CHD) is a major cause of infant mortality and morbidity, yet the genetic causes and mechanisms remain opaque. In a patient with CHD and heterotaxy, a disorder of left-right (LR) patterning, a de novo mutation was identified in the chromatin modifier gene WDR5 WDR5 acts as a scaffolding protein in the H3K4 methyltransferase complex, but a role in LR patterning is unknown. Here, we show that Wdr5 depletion leads to LR patterning defects in Xenopus via its role in ciliogenesis. Unexpectedly, we find a dual role for WDR5 in LR patterning. First, WDR5 is expressed in the nuclei of monociliated cells of the LR organizer (LRO) and regulates foxj1 expression. LR defects in wdr5 morphants can be partially rescued with the addition of foxj1 Second, WDR5 localizes to the bases of cilia. Using a mutant form of WDR5, we demonstrate that WDR5 also has an H3K4-independent role in LR patterning. Guided by the patient phenotype, we identify multiple roles for WDR5 in LR patterning, providing plausible mechanisms for its role in ciliopathies like heterotaxy and CHD.

CDC20B is required for deuterosome-mediated centriole production in multiciliated cells

Multiciliated cells (MCCs) harbor dozens to hundreds of motile cilia, which generate hydrodynamic forces important in animal physiology. In vertebrates, MCC differentiation involves massive centriole production by poorly characterized structures called deuterosomes. Here, single-cell RNA sequencing reveals that human deuterosome stage MCCs are characterized by the expression of many cell cycle-related genes. We further investigated the uncharacterized vertebrate-specific cell division cycle 20B (CDC20B) gene, which hosts microRNA-449abc. We show that CDC20B protein associates to deuterosomes and is required for centriole release and subsequent cilia production in mouse and Xenopus MCCs. CDC20B interacts with PLK1, a kinase known to coordinate centriole disengagement with the protease Separase in mitotic cells. Strikingly, over-expression of Separase rescues centriole disengagement and cilia production in CDC20B-deficient MCCs. This work reveals the shaping of deuterosome-mediated centriole production in vertebrate MCCs, by adaptation of canonical and recently evolved cell cycle-related molecules.

Auto-fatty acylation of transcription factor RFX3 regulates ciliogenesis

Defects in cilia have been associated with an expanding human disease spectrum known as ciliopathies. Regulatory Factor X 3 (RFX3) is one of the major transcription factors required for ciliogenesis and cilia functions. In addition, RFX3 regulates pancreatic islet cell differentiation and mature β-cell functions. However, how RFX3 protein is regulated at the posttranslational level remains poorly understood. Using chemical reporters of protein fatty acylation and mass spectrometry analysis, here we show that RFX3 transcriptional activity is regulated by S-fatty acylation at a highly conserved cysteine residue in the dimerization domain. Surprisingly, RFX3 undergoes enzyme-independent, "self-catalyzed" auto-fatty acylation and displays preferences for 18-carbon stearic acid and oleic acid. The fatty acylation-deficient mutant of RFX3 shows decreased homodimerization; fails to promote ciliary gene expression, ciliogenesis, and elongation; and impairs Hedgehog signaling. Our findings reveal a regulation of RFX3 transcription factor and link fatty acid metabolism and protein lipidation to the regulation of ciliogenesis.

High doses of CRISPR/Cas9 ribonucleoprotein efficiently induce gene knockout with low mosaicism in the hydrozoan Clytia hemisphaerica through microhomology-mediated deletion

Targeted mutagenesis using CRISPR/Cas9 technology has been shown to be a powerful approach to examine gene function in diverse metazoan species. One common drawback is that mixed genotypes, and thus variable phenotypes, arise in the F0 generation because incorrect DNA repair produces different mutations amongst cells of the developing embryo. We report here an effective method for gene knockout (KO) in the hydrozoan Clytia hemisphaerica, by injection into the egg of Cas9/sgRNA ribonucleoprotein complex (RNP). Expected phenotypes were observed in the F0 generation when targeting endogenous GFP genes, which abolished fluorescence in embryos, or CheRfx123 (that codes for a conserved master transcriptional regulator for ciliogenesis) which caused sperm motility defects. When high concentrations of Cas9 RNP were used, the mutations in target genes at F0 polyp or jellyfish stages were not random but consisted predominantly of one or two specific deletions between pairs of short microhomologies flanking the cleavage site. Such microhomology-mediated (MM) deletion is most likely caused by microhomology-mediated end-joining (MMEJ), which may be favoured in early stage embryos. This finding makes it very easy to isolate uniform, largely non-mosaic mutants with predictable genotypes in the F0 generation in Clytia, allowing rapid and reliable phenotype assessment.

Characterization of the human RFX transcription factor family by regulatory and target gene analysis

BACKGROUND

Evolutionarily conserved RFX transcription factors (TFs) regulate their target genes through a DNA sequence motif called the X-box. Thereby they regulate cellular specialization and terminal differentiation. Here, we provide a comprehensive analysis of all the eight human RFX genes (RFX1-8), their spatial and temporal expression profiles, potential upstream regulators and target genes.

RESULTS

We extracted all known human RFX1-8 gene expression profiles from the FANTOM5 database derived from transcription start site (TSS) activity as captured by Cap Analysis of Gene Expression (CAGE) technology. RFX genes are broadly (RFX1-3, RFX5, RFX7) and specifically (RFX4, RFX6) expressed in different cell types, with high expression in four organ systems: immune system, gastrointestinal tract, reproductive system and nervous system. Tissue type specific expression profiles link defined RFX family members with the target gene batteries they regulate. We experimentally confirmed novel TSS locations and characterized the previously undescribed RFX8 to be lowly expressed. RFX tissue and cell type specificity arises mainly from differences in TSS architecture. RFX transcript isoforms lacking a DNA binding domain (DBD) open up new possibilities for combinatorial target gene regulation. Our results favor a new grouping of the RFX family based on protein domain composition. We uncovered and experimentally confirmed the TFs SP2 and ESR1 as upstream regulators of specific RFX genes. Using TF binding profiles from the JASPAR database, we determined relevant patterns of X-box motif positioning with respect to gene TSS locations of human RFX target genes.

CONCLUSIONS

The wealth of data we provide will serve as the basis for precisely determining the roles RFX TFs play in human development and disease.

Zebrafish Rfx4 controls dorsal and ventral midline formation in the neural tube

BACKGROUND

Rfx winged-helix transcription factors, best known as key regulators of core ciliogenesis, also play ciliogenesis-independent roles during neural development. Mammalian Rfx4 controls neural tube morphogenesis via both mechanisms.

RESULTS

We set out to identify conserved aspects of rfx4 gene function during vertebrate development and to establish a new genetic model in which to analyze these mechanisms further. To this end, we have generated frame-shift alleles in the zebrafish rfx4 locus using CRISPR/Cas9 mutagenesis. Using RNAseq-based transcriptome analysis, in situ hybridization and immunostaining we identified a requirement for zebrafish rfx4 in the forming midlines of the caudal neural tube. These functions are mediated, least in part, through transcriptional regulation of several zic genes in the dorsal hindbrain and of foxa2 in the ventral hindbrain and spinal cord (floor plate).

CONCLUSIONS

The midline patterning functions of rfx4 are conserved, because rfx4 regulates transcription of foxa2 and zic2 in zebrafish and in mouse. In contrast, zebrafish rfx4 function is dispensable for forebrain morphogenesis, while mouse rfx4 is required for normal formation of forebrain ventricles in a ciliogenesis-dependent manner. Collectively, this report identifies conserved aspects of rfx4 function and establishes a robust new genetic model for in-depth dissection of these mechanisms. Developmental Dynamics 247:650-659, 2018. © 2017 Wiley Periodicals, Inc.

Evolutionary Proteomics Uncovers Ancient Associations of Cilia with Signaling Pathways

Cilia are organelles specialized for movement and signaling. To infer when during evolution signaling pathways became associated with cilia, we characterized the proteomes of cilia from sea urchins, sea anemones, and choanoflagellates. We identified 437 high-confidence ciliary candidate proteins conserved in mammals and discovered that Hedgehog and G-protein-coupled receptor pathways were linked to cilia before the origin of bilateria and transient receptor potential (TRP) channels before the origin of animals. We demonstrated that candidates not previously implicated in ciliary biology localized to cilia and further investigated ENKUR, a TRP channel-interacting protein identified in the cilia of all three organisms. ENKUR localizes to motile cilia and is required for patterning the left-right axis in vertebrates. Moreover, mutation of ENKUR causes situs inversus in humans. Thus, proteomic profiling of cilia from diverse eukaryotes defines a conserved ciliary proteome, reveals ancient connections to signaling, and uncovers a ciliary protein that underlies development and human disease.

Systematic target function annotation of human transcription factors

Background

Transcription factors (TFs), the key players in transcriptional regulation, have attracted great experimental attention, yet the functions of most human TFs remain poorly understood. Recent capabilities in genome-wide protein binding profiling have stimulated systematic studies of the hierarchical organization of human gene regulatory network and DNA-binding specificity of TFs, shedding light on combinatorial gene regulation. We show here that these data also enable a systematic annotation of the biological functions and functional diversity of TFs.

Result

We compiled a human gene regulatory network for 384 TFs covering the 146,096 TF–target gene (TF–TG) relationships, extracted from over 850 ChIP-seq experiments as well as the literature. By integrating this network of TF–TF and TF–TG relationships with 3715 functional concepts from six sources of gene function annotations, we obtained over 9000 confident functional annotations for 279 TFs. We observe extensive connectivity between TFs and Mendelian diseases, GWAS phenotypes, and pharmacogenetic pathways. Further, we show that TFs link apparently unrelated functions, even when the two functions do not share common genes. Finally, we analyze the pleiotropic functions of TFs and suggest that the increased number of upstream regulators contributes to the functional pleiotropy of TFs.

Conclusion

Our computational approach is complementary to focused experimental studies on TF functions, and the resulting knowledge can guide experimental design for the discovery of unknown roles of TFs in human disease and drug response.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-017-0469-0) contains supplementary material, which is available to authorized users.

Conditional ablation of the RFX4 isoform 1 transcription factor: Allele dosage effects on brain phenotype

Regulatory factor X4 (RFX4) isoform 1 is a recently discovered isoform of the winged helix transcription factor RFX4, which can bind to X-box consensus sequences that are enriched in the promoters of cilia-related genes. Early insertional mutagenesis studies in mice first identified this isoform, and demonstrated that it was crucial for mouse brain development. RFX4 isoform 1 is the only RFX4 isoform significantly expressed in the mouse fetal and adult brain. In this study, we evaluated conditional knock-out (KO) mice in which one or two floxed alleles of Rfx4 were deleted early in development through the use of a Sox2-Cre transgene. Heterozygous deletion of Rfx4 resulted in severe, non-communicating congenital hydrocephalus associated with hypoplasia of the subcommissural organ. Homozygous deletion of Rfx4 resulted in formation of a single ventricle in the forebrain, and severe dorsoventral patterning defects in the telencephalon and midbrain at embryonic day 12.5, a collection of phenotypes that resembled human holoprosencephaly. No anatomical abnormalities were noted outside the brain in either case. At the molecular level, transcripts encoded by the cilia-related gene Foxj1 were significantly decreased, and Foxj1 was identified as a direct gene target of RFX4 isoform 1. The phenotypes were similar to those observed in the previous Rfx4 insertional mutagenesis studies. Thus, we provide a novel conditional KO animal model in which to investigate the downstream genes directly and/or indirectly regulated by RFX4 isoform 1. This model could provide new insights into the pathogenesis of obstructive hydrocephalus and holoprosencephaly in humans, both relatively common and disabling birth defects.

Vertebrate Left-Right Asymmetry: What Can Nodal Cascade Gene Expression Patterns Tell Us?

Laterality of inner organs is a wide-spread characteristic of vertebrates and beyond. It is ultimately controlled by the left-asymmetric activation of the Nodal signaling cascade in the lateral plate mesoderm of the neurula stage embryo, which results from a cilia-driven leftward flow of extracellular fluids at the left-right organizer. This scenario is widely accepted for laterality determination in wildtype specimens. Deviations from this norm come in different flavors. At the level of organ morphogenesis, laterality may be inverted (situs inversus) or non-concordant with respect to the main body axis (situs ambiguus or heterotaxia). At the level of Nodal cascade gene activation, expression may be inverted, bilaterally induced, or absent. In a given genetic situation, patterns may be randomized or predominantly lacking laterality (absence or bilateral activation). We propose that the distributions of patterns observed may be indicative of the underlying molecular defects, with randomizations being primarily caused by defects in the flow-generating ciliary set-up, and symmetrical patterns being the result of impaired flow sensing, on the left, the right, or both sides. This prediction, the reasoning of which is detailed in this review, pinpoints functions of genes whose role in laterality determination have remained obscure.

Controlling centriole numbers: Geminin family members as master regulators of centriole amplification and multiciliogenesis

To ensure that the genetic material is accurately passed down to daughter cells during mitosis, dividing cells must duplicate their chromosomes and centrosomes once and only once per cell cycle. The same key steps-licensing, duplication, and segregation-control both the chromosome and the centrosome cycle, which must occur in concert to safeguard genome integrity. Aberrations in genome content or centrosome numbers lead to genomic instability and are linked to tumorigenesis. Such aberrations, however, can also be part of the normal life cycle of specific cell types. Multiciliated cells best exemplify the deviation from a normal centrosome cycle. They are post-mitotic cells which massively amplify their centrioles, bypassing the rule for once-per-cell-cycle centriole duplication. Hundreds of centrioles dock to the apical cell surface and generate motile cilia, whose concerted movement ensures fluid flow across epithelia. The early steps that control the generation of multiciliated cells have lately started to be elucidated. Geminin and the vertebrate-specific GemC1 and McIdas are distantly related coiled-coil proteins, initially identified as cell cycle regulators associated with the chromosome cycle. Geminin is required to ensure once-per-cell-cycle genome replication, while McIdas and GemC1 bind to Geminin and are implicated in DNA replication control. Recent findings highlight Geminin family members as early regulators of multiciliogenesis. GemC1 and McIdas specify the multiciliate cell fate by forming complexes with the E2F4/5 transcription factors to switch on a gene expression program leading to centriole amplification and cilia formation. Positive and negative interactions among Geminin family members may link cell cycle control to centriole amplification and multiciliogenesis, acting close to the point of transition from proliferation to differentiation. We review key steps of centrosome duplication and amplification, present the role of Geminin family members in the centrosome and chromosome cycle, and discuss links with disease.

hmmr mediates anterior neural tube closure and morphogenesis in the frog Xenopus

Development of the central nervous system requires orchestration of morphogenetic processes which drive elevation and apposition of the neural folds and their fusion into a neural tube. The newly formed tube gives rise to the brain in anterior regions and continues to develop into the spinal cord posteriorly. Conspicuous differences between the anterior and posterior neural tube become visible already during neural tube closure (NTC). Planar cell polarity (PCP)-mediated convergent extension (CE) movements are restricted to the posterior neural plate, i.e. hindbrain and spinal cord, where they propagate neural fold apposition. The lack of CE in the anterior neural plate correlates with a much slower mode of neural fold apposition anteriorly. The morphogenetic processes driving anterior NTC have not been addressed in detail. Here, we report a novel role for the breast cancer susceptibility gene and microtubule (MT) binding protein Hmmr (Hyaluronan-mediated motility receptor, RHAMM) in anterior neurulation and forebrain development in Xenopus laevis. Loss of hmmr function resulted in a lack of telencephalic hemisphere separation, arising from defective roof plate formation, which in turn was caused by impaired neural tissue narrowing. hmmr regulated polarization of neural cells, a function which was dependent on the MT binding domains. hmmr cooperated with the core PCP component vangl2 in regulating cell polarity and neural morphogenesis. Disrupted cell polarization and elongation in hmmr and vangl2 morphants prevented radial intercalation (RI), a cell behavior essential for neural morphogenesis. Our results pinpoint a novel role of hmmr in anterior neural development and support the notion that RI is a major driving force for anterior neurulation and forebrain morphogenesis.

Genes and molecular pathways underpinning ciliopathies

Motile and non-motile (primary) cilia are nearly ubiquitous cellular organelles. The dysfunction of cilia causes diseases known as ciliopathies. The number of reported ciliopathies (currently 35) is increasing, as is the number of established (187) and candidate (241) ciliopathy-associated genes. The characterization of ciliopathy-associated proteins and phenotypes has improved our knowledge of ciliary functions. In particular, investigating ciliopathies has helped us to understand the molecular mechanisms by which the cilium-associated basal body functions in early ciliogenesis, as well as how the transition zone functions in ciliary gating, and how intraflagellar transport enables cargo trafficking and signalling. Both basic biological and clinical studies are uncovering novel ciliopathies and the ciliary proteins involved. The assignment of these proteins to different ciliary structures, processes and ciliopathy subclasses (first order and second order) provides insights into how this versatile organelle is built, compartmentalized and functions in diverse ways that are essential for human health.

The development and functions of multiciliated epithelia

Multiciliated cells are epithelial cells that are in contact with bodily fluids and are required for the proper function of major organs including the brain, the respiratory system and the reproductive tracts. Their multiple motile cilia beat unidirectionally to remove particles of external origin from their surface and/or drive cells or fluids into the lumen of the organs. Multiciliated cells in the brain are produced once, almost exclusively during embryonic development, whereas in respiratory tracts and oviducts they regenerate throughout life. In this Review, we provide a cell-to-organ overview of multiciliated cells and highlight recent studies that have greatly increased our understanding of the mechanisms driving the development and function of these cells in vertebrates. We discuss cell fate determination and differentiation of multiciliated cells, and provide a comprehensive account of their locations and functions in mammals.

La-related protein 6 controls ciliated cell differentiation

Background

La-related protein 6 (LARP6) is an evolutionally conserved RNA-binding protein. Vertebrate LARP6 binds the 5′ stem-loop found in mRNAs encoding type I collagen to regulate their translation, but other target mRNAs and additional functions for LARP6 are unknown. The aim of this study was to elucidate an additional function of LARP6 and to evaluate the importance of its function during development.

Methods

To uncover the role of LARP6 in development, we utilized Morpholino Oligos to deplete LARP6 protein in Xenopus embryos. Then, embryonic phenotypes and ciliary structures of LAPR6 morphants were examined. To identify the molecular mechanism underlying ciliogenesis regulated by LARP6, we tested the expression level of cilia-related genes, which play important roles in ciliogenesis, by RT-PCR or whole mount in situ hybridization (WISH).

Results

We knocked down LARP6 in Xenopus embryos and found neural tube closure defects. LARP6 mutant, which compromises the collagen synthesis, could rescue these defects. Neural tube closure defects are coincident with lack of cilia, antenna-like cellular organelles with motility- or sensory-related functions, in the neural tube. The absence of cilia at the epidermis was also observed in LARP6 morphants, and this defect was due to the absence of basal bodies which are formed from centrioles and required for ciliary assembly. In the process of multi-ciliated cell (MCC) differentiation, mcidas, which activates the transcription of genes required for centriole formation during ciliogenesis, could partially restore MCCs in LARP6 morphants. In addition, LARP6 likely controls the expression of mcidas in a Notch-independent manner.

Conclusions

La-related protein 6 is involved in ciliated cell differentiation during development by controlling the expression of cilia-related genes including mcidas. This LARP6 function involves a mechanism that is distinct from its established role in binding to collagen mRNAs and regulating their translation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13630-017-0047-7) contains supplementary material, which is available to authorized users.

Rfx2 Stabilizes Foxj1 Binding at Chromatin Loops to Enable Multiciliated Cell Gene Expression

Cooperative transcription factor binding at cis-regulatory sites in the genome drives robust eukaryotic gene expression, and many such sites must be coordinated to produce coherent transcriptional programs. The transcriptional program leading to motile cilia formation requires members of the DNA-binding forkhead (Fox) and Rfx transcription factor families and these factors co-localize to cilia gene promoters, but it is not clear how many cilia genes are regulated by these two factors, whether these factors act directly or indirectly, or how these factors act with specificity in the context of a 3-dimensional genome. Here, we use genome-wide approaches to show that cilia genes reside at the boundaries of topological domains and that these areas have low enhancer density. We show that the transcription factors Foxj1 and Rfx2 binding occurs in the promoters of more cilia genes than other known cilia transcription factors and that while Rfx2 binds directly to promoters and enhancers equally, Foxj1 prefers direct binding to enhancers and is stabilized at promoters by Rfx2. Finally, we show that Rfx2 and Foxj1 lie at the anchor endpoints of chromatin loops, suggesting that target genes are activated when Foxj1 bound at distal sites is recruited via a loop created by Rfx2 binding at both sites. We speculate that the primary function of Rfx2 is to stabilize distal enhancers with proximal promoters by operating as a scaffolding factor, bringing key regulatory domains bound by Foxj1 into close physical proximity and enabling coordinated cilia gene expression.

The multiciliated cell extends hundreds of motile cilia to produce fluid flow in the airways and other organ systems. The formation of this specialized cell type requires the coordinated expression of hundreds of genes in order to produce all the protein parts motile cilia require. While a relatively small number of transcription factors has been identified that promote gene expression during multiciliate cell differentiation, it is not clear how they work together to coordinate the expression of genes required for multiple motile ciliation. Here, we show that two transcription factors known to drive cilia formation, Foxj1 and Rfx2, play complementary roles wherein Foxj1 activates target genes but tends not to bind near them in the genome, whereas Rfx2 can’t activate target genes by itself but instead acts as a scaffold by localizing Foxj1 to the proper targets. These results suggest not only a mechanism by which complex gene expression is coordinated in multiciliated cells, but also how transcriptional programs in general could be modular and deployed across different cellular contexts with the same basic promoter configuration.