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/).

ConVarT: a search engine for matching human genetic variants with variants from non-human species.

The availability of genetic variants, together with phenotypic annotations from model organisms, facilitates comparing these variants with equivalent variants in humans. However, existing databases and search tools do not make it easy to scan for equivalent variants, namely 'matching variants' (MatchVars) between humans and other organisms. Therefore, we developed an integrated search engine called ConVarT (http://www.convart.org/) for matching variants between humans, mice, and Caenorhabditis elegans. ConVarT incorporates annotations (including phenotypic and pathogenic) into variants, and these previously unexploited phenotypic MatchVars from mice and C. elegans can give clues about the functional consequence of human genetic variants. Our analysis shows that many phenotypic variants in different genes from mice and C. elegans, so far, have no counterparts in humans, and thus, can be useful resources when evaluating a relationship between a new human mutation and a disease.

G-quadruplex prediction in E. coli genome reveals a conserved putative G-quadruplex-Hairpin-Duplex switch.

Many studies show that short non-coding sequences are widely conserved among regulatory elements. More and more conserved sequences are being discovered since the development of next generation sequencing technology. A common approach to identify conserved sequences with regulatory roles relies on topological changes such as hairpin formation at the DNA or RNA level. G-quadruplexes, non-canonical nucleic acid topologies with little established biological roles, are increasingly considered for conserved regulatory element discovery. Since the tertiary structure of G-quadruplexes is strongly dependent on the loop sequence which is disregarded by the generally accepted algorithm, we hypothesized that G-quadruplexes with similar topology and, indirectly, similar interaction patterns, can be determined using phylogenetic clustering based on differences in the loop sequences. Phylogenetic analysis of 52 G-quadruplex forming sequences in the Escherichia coli genome revealed two conserved G-quadruplex motifs with a potential regulatory role. Further analysis revealed that both motifs tend to form hairpins and G quadruplexes, as supported by circular dichroism studies. The phylogenetic analysis as described in this work can greatly improve the discovery of functional G-quadruplex structures and may explain unknown regulatory patterns.

Active transport and diffusion barriers restrict Joubert Syndrome-associated ARL13B/ARL-13 to an Inv-like ciliary membrane subdomain.

Cilia are microtubule-based cell appendages, serving motility, chemo-/mechano-/photo- sensation, and developmental signaling functions. Cilia are comprised of distinct structural and functional subregions including the basal body, transition zone (TZ) and inversin (Inv) compartments, and defects in this organelle are associated with an expanding spectrum of inherited disorders including Bardet-Biedl syndrome (BBS), Meckel-Gruber Syndrome (MKS), Joubert Syndrome (JS) and Nephronophthisis (NPHP). Despite major advances in understanding ciliary trafficking pathways such as intraflagellar transport (IFT), how proteins are transported to subciliary membranes remains poorly understood. Using Caenorhabditis elegans and mammalian cells, we investigated the transport mechanisms underlying compartmentalization of JS-associated ARL13B/ARL-13, which we previously found is restricted at proximal ciliary membranes. We now show evolutionary conservation of ARL13B/ARL-13 localisation to an Inv-like subciliary membrane compartment, excluding the TZ, in many C. elegans ciliated neurons and in a subset of mammalian ciliary subtypes. Compartmentalisation of C. elegans ARL-13 requires a C-terminal RVVP motif and membrane anchoring to prevent distal cilium and nuclear targeting, respectively. Quantitative imaging in more than 20 mutants revealed differential contributions for IFT and ciliopathy modules in defining the ARL-13 compartment; IFT-A/B, IFT-dynein and BBS genes prevent ARL-13 accumulation at periciliary membranes, whereas MKS/NPHP modules additionally inhibit ARL-13 association with TZ membranes. Furthermore, in vivo FRAP analyses revealed distinct roles for IFT and MKS/NPHP genes in regulating a TZ barrier to ARL-13 diffusion, and intraciliary ARL-13 diffusion. Finally, C. elegans ARL-13 undergoes IFT-like motility and quantitative protein complex analysis of human ARL13B identified functional associations with IFT-B complexes, mapped to IFT46 and IFT74 interactions. Together, these findings reveal distinct requirements for sequence motifs, IFT and ciliopathy modules in defining an ARL-13 subciliary membrane compartment. We conclude that MKS/NPHP modules comprise a TZ barrier to ARL-13 diffusion, whereas IFT genes predominantly facilitate ARL-13 ciliary entry and/or retention via active transport mechanisms.

CC2D1A causes ciliopathy, intellectual disability, heterotaxy, renal dysplasia, and abnormal CSF flow.

Intellectual and developmental disabilities result from abnormal nervous system development. Over a 1,000 genes have been associated with intellectual and developmental disabilities, driving continued efforts toward dissecting variant functionality to enhance our understanding of the disease mechanism. This report identified two novel variants in CC2D1A in a cohort of four patients from two unrelated families. We used multiple model systems for functional analysis, including Xenopus, Drosophila, and patient-derived fibroblasts. Our experiments revealed that cc2d1a is expressed explicitly in a spectrum of ciliated tissues, including the left-right organizer, epidermis, pronephric duct, nephrostomes, and ventricular zone of the brain. In line with this expression pattern, loss of cc2d1a led to cardiac heterotaxy, cystic kidneys, and abnormal CSF circulation via defective ciliogenesis. Interestingly, when we analyzed brain development, mutant tadpoles showed abnormal CSF circulation only in the midbrain region, suggesting abnormal local CSF flow. Furthermore, our analysis of the patient-derived fibroblasts confirmed defective ciliogenesis, further supporting our observations. In summary, we revealed novel insight into the role of CC2D1A by establishing its new critical role in ciliogenesis and CSF circulation.

RPI-1 (human DCDC2) displays functional redundancy with Nephronophthisis 4 in regulating cilia biogenesis in C. elegans.

Projecting from most cell surfaces, cilia serve as important hubs for sensory and signaling processes and have been linked to a variety of human disorders, including Bardet-Biedl Syndrome (BBS), Meckel-Gruber Syndrome (MKS), Nephronophthisis (NPHP), and Joubert Syndrome, and these diseases are collectively known as a ciliopathy. DCDC2 is a ciliopathy protein that localizes to cilia; nevertheless, our understanding of the role of DCDC2 in cilia is still limited. We employed C. elegans to investigate the function of C. elegans RPI-1, a Caenorhabditis elegans ortholog of human DCDC2, in cilia and found that C. elegans RPI-1 localizes to the entire ciliary axoneme, but is not present in the transition zone and basal body. We generated a null mutant of C. elegans rpi-1, and our analysis with a range of fluorescence-based ciliary markers revealed that DCDC2 and nephronophthisis 4 (NPHP-4/NPHP4) display functional redundant roles in regulating cilia length and cilia positions. Taken together, our analysis discovered a novel genetic interaction between two ciliopathy disease genes (RPI-1/DCDC2 and NPHP-4/NPHP4) in C. elegans.

CiliaMiner: an integrated database for ciliopathy genes and ciliopathies.

Cilia are found in eukaryotic species ranging from single-celled organisms, such as Chlamydomonas reinhardtii, to humans, but not in plants. The ability to respond to repellents and/or attractants, regulate cell proliferation and differentiation and provide cellular mobility are just a few examples of how crucial cilia are to cells and organisms. Over 30 distinct rare disorders generally known as ciliopathy are caused by abnormalities or functional impairments in cilia and cilia-related compartments. Because of the complexity of ciliopathies and the rising number of ciliopathies and ciliopathy genes, a ciliopathy-oriented and up-to-date database is required. Here, we present CiliaMiner, a manually curated ciliopathy database that includes ciliopathy lists collected from articles and databases. Analysis reveals that there are 55 distinct disorders likely related to ciliopathy, with over 4000 clinical manifestations. Based on comparative symptom analysis and subcellular localization data, diseases are classified as primary, secondary or atypical ciliopathies. CiliaMiner provides easy access to all of these diseases and disease genes, as well as clinical features and gene-specific clinical features, as well as subcellular localization of each protein. Additionally, the orthologs of disease genes are also provided for mice, zebrafish, Xenopus, Drosophila, Caenorhabditis elegans and Chlamydomonas reinhardtii. CiliaMiner (https://kaplanlab.shinyapps.io/ciliaminer) aims to serve the cilia community with its comprehensive content and highly enriched interactive heatmaps, and will be continually updated. Database URL: https://kaplanlab.shinyapps.io/ciliaminer/.

Matching variants for functional characterization of genetic variants.

Rapid and low-cost sequencing, as well as computer analysis, have facilitated the diagnosis of many genetic diseases, resulting in a substantial rise in the number of disease-associated genes. However, genetic diagnosis of many disorders remains problematic due to the lack of interpretation for many genetic variants, especially missenses, the infeasibility of high-throughput experiments on mammals, and the shortcomings of computational prediction technologies. Additionally, the available mutant databases are not well-utilized. Toward this end, we used Caenorhabditis elegans mutant resources to delineate the functions of eight missense variants (V444I, V517D, E610K, L732F, E817K, H873P, R1105K, and G1205E) and two stop codons (W937stop and Q1434stop), including several matching variants (MatchVar) with human in ciliopathy associated IFT-140 (also called CHE-11)//IFT140 (intraflagellar transport protein 140). Moreover, MatchVars carrying C. elegans mutants, including IFT-140(G680S) and IFT-140(P702A) for the human (G704S) (dbSNP: rs150745099) and P726A (dbSNP: rs1057518064 and a conflicting variation) were created using CRISPR/Cas9. IFT140 is a key component of IFT complex A (IFT-A), which is involved in the retrograde transport of IFT along cilia and the entrance of G protein-coupled receptors into cilia. Functional analysis of all 10 variants revealed that P702A and W937stop, but not others phenocopied the ciliary phenotypes (short cilia, IFT accumulations, mislocalization of membrane proteins, and cilia entry of nonciliary proteins) of the IFT-140 null mutant, indicating that both P702A and W937stop are phenotypic in C. elegans. Our functional data offered experimental support for interpreting human variants, by using ready-to-use mutants carrying MatchVars and generating MatchVars with CRISPR/Cas9.

WDR31 displays functional redundancy with GTPase-activating proteins (GAPs) ELMOD and RP2 in regulating IFT complex and recruiting the BBSome to cilium.

The correct intraflagellar transport (IFT) assembly at the ciliary base and the IFT turnaround at the ciliary tip are key for the IFT to perform its function, but we still have poor understanding about how these processes are regulated. Here, we identify WDR31 as a new ciliary protein, and analysis from zebrafish and Caenorhabditis elegans reveals the role of WDR31 in regulating the cilia morphology. We find that loss of WDR-31 together with RP-2 and ELMD-1 (the sole ortholog ELMOD1-3) results in ciliary accumulations of IFT Complex B components and KIF17 kinesin, with fewer IFT/BBSome particles traveling along cilia in both anterograde and retrograde directions, suggesting that the IFT/BBSome entry into the cilia and exit from the cilia are impacted. Furthermore, anterograde IFT in the middle segment travels at increased speed in wdr-31;rpi-2;elmd-1 Remarkably, a non-ciliary protein leaks into the cilia of wdr-31;rpi-2;elmd-1, possibly because of IFT defects. This work reveals WDR31-RP-2-ELMD-1 as IFT and BBSome trafficking regulators.

The AP-1 clathrin adaptor facilitates cilium formation and functions with RAB-8 in C. elegans ciliary membrane transport.

Clathrin adaptor (AP) complexes facilitate membrane trafficking between subcellular compartments. One such compartment is the cilium, whose dysfunction underlies disorders classified as ciliopathies. Although AP-1mu subunit (UNC-101) is linked to cilium formation and targeting of transmembrane proteins (ODR-10) to nematode sensory cilia at distal dendrite tips, these functions remain poorly understood. Here, using Caenorhabditis elegans sensory neurons and mammalian cell culture models, we find conservation of AP-1 function in facilitating cilium morphology, positioning and orientation, and microtubule stability and acetylation. These defects appear to be independent of IFT, because AP-1-depleted cells possess normal IFT protein localisation and motility. By contrast, disruption of chc-1 (clathrin) or rab-8 phenocopies unc-101 worms, preventing ODR-10 vesicle formation and causing misrouting of ODR-10 to all plasma membrane destinations. Finally, ODR-10 colocalises with RAB-8 in cell soma and they cotranslocate along dendrites, whereas ODR-10 and UNC-101 signals do not overlap. Together, these data implicate conserved roles for metazoan AP-1 in facilitating cilium structure and function, and suggest cooperation with RAB-8 to coordinate distinct early steps in neuronal ciliary membrane sorting and trafficking.

ConVarT: Search Engine for Missense Variants Between Humans and Other Organisms.

ConVarT (https://convart.org/) is a search engine for searching for conjugate variants between humans and other species. The search engine is based on matching conjugate variants called MatchVars between species. Matching equivalent variants requires correct alignment of orthologous proteins with the use of multiple sequence alignments (MSA). Indeed, the ConVarT pipeline has performed over a million MSAs and integrated variants and variant-specific annotations (pathogenicity, phenotypic variants; etc.) into the corresponding positions on MSAs. When a clinically relevant variant is discovered whose functional relevance is unknown, ConVarT offers clinician scientists the possibility to search for a MatchVar in other species and to look for functional data on that variant. Fortunately, ConVarT enables users to paste a protein sequence in FASTA format to search for human orthologous proteins. A pairwise sequence alignment (PSA) is then performed between the provided protein sequence and the human orthologous protein, allowing users to visualize human variants on the PSA. Here, we describe the step-by-step usage of ConVarT. © 2022 Wiley Periodicals LLC. Basic Protocol 1: Searching matching variants (MatchVar) with gene/protein identifiers. Basic Protocol 2: Searching with a FASTA sequence. Alternate Protocol: Search with gene name in multiple species. Basic Protocol 3: Search genes associated with a disease.

Protocol for determining the average speed and frequency of kinesin and dynein-driven intraflagellar transport (IFT) in C. elegans.

Here, we present a protocol to image a fluorescent-labeled intraflagellar transport (IFT) component in Caenorhabditis elegans with fluorescence microscopy, including steps of sample preparations, in vivo live-cell imaging, and post-microscopy analysis with kymographs. This protocol breaks down all processes into three categories: (1) pre-imaging preparations, (2) preparations for the time of image acquisition, and (3) post-imaging analyses. The protocol can be applied to determine the speed and frequency of moving particles. For complete details on the use and execution of this protocol, please refer to Cevik et al. (2021).

The Joubert syndrome protein CEP41 is excluded from the distal segment of cilia in C. elegans.

Rare diseases are a fundamental issue in today's world, affecting more than 300 million individuals worldwide. According to data from Orphanet and OMIM, about 50-60 new conditions are added to the list of over 6,000 clinically distinct diseases each year, rendering disease diagnosis and treatment even more challenging. Ciliopathies comprise a heterogeneous category of rare diseases made up of over 35 distinct diseases, including Joubert syndrome (JBTS; OMIM 213300), that are caused by functional and structural defects in cilia. JBTS is an autosomal recessive condition characterized by a range of symptoms, including cerebellar vermis hypoplasia and poor muscle tone. There are now a total of 38 genes that cause JBTS, almost all of which encode protein products that are found in cilia and cilia-associated compartments, such as the basal body and transition zone. CEP41 is a JBTS-associated protein that is found in cilia and the basal body of mammals, but its localization in other ciliary organisms remains elusive. C. elegans is an excellent model organism for studying the molecular mechanisms of rare diseases like JBTS. We, therefore, decided to use C. elegans to identify the localization of CEP41. Our microscopy analysis revealed that CEPH-41(CEntrosomal Protein Homolog 41) not only localizes to cilia but is excluded from the distal segment of the amphid and phasmid cilia in C. elegans. Furthermore, we discovered a putative X-box motif located in the promoter of ceph-41 and the expression of ceph-41 is regulated by DAF-19, a sole Regulatory Factor X (RFX) transcription factor.

Subcellular localization of the voltage-gated K+ channel EGL-36 , a member of the KV3 subfamily, in the ciliated sensory neurons in C. elegans.

Delineated as the first cellular organelle in 1675 by Antonie van Leeuwenhoek, cilia did not receive much attention until the 2000s, when it became apparent that cilia played a key role in the development of embryos, a variety of signaling pathways. Therefore, collective efforts by many scientists have led to the identification of many novel ciliopathy and cilia genes, while we are still far from disclosing the complete components of cilia.Here we used the ciliated sensory neurons in C. elegans as a model system that revealed the voltage-gated K+ channel EGL-36 (a member of the Shaw subfamily) as a new component associated with cilia. The confocal microscopy examination of fluorescence tagged EGL-36 together with ciliary (IFT-140) or transition zone (MKS-6) markers reveal that EGL-36 is only expressed in subsets of the ciliated sensory neurons, where it partially overlaps with the basal body signals and predominantly localizes to the periciliary membrane compartment. This expression pattern along with studies of egl-36 gain-of-function variants indicates that egl-36 is not essential for ciliogenesis in C. elegans. Our data identify the voltage-gated K+ channel EGL-36 as a new cilia-associated protein, and future studies should reveal the functional significance of EGL-36 in cilia biogenesis.

MSABrowser: dynamic and fast visualization of sequence alignments, variations and annotations.

Summary: Sequence alignment is an excellent way to visualize the similarities and differences between DNA, RNA or protein sequences, yet it is currently difficult to jointly view sequence alignment data with genetic variations, modifications such as post-translational modifications and annotations (i.e. protein domains). Here, we present the MSABrowser tool that makes it easy to co-visualize genetic variations, modifications and annotations on the respective positions of amino acids or nucleotides in pairwise or multiple sequence alignments. MSABrowser is developed entirely in JavaScript and works on any modern web browser at any platform, including Linux, Mac OS X and Windows systems without any installation. MSABrowser is also freely available for the benefit of the scientific community. Availability and implementation: MSABrowser is released as open-source and web-based software under MIT License. The visualizer, documentation, all source codes and examples are available at https://thekaplanlab.github.io/ and GitHub repository https://github.com/thekaplanlab/msabrowser. Supplementary information: Supplementary data are available at Bioinformatics Advances online.

Endocytosis genes facilitate protein and membrane transport in C. elegans sensory cilia.

BACKGROUND: Multiple intracellular transport pathways drive the formation, maintenance, and function of cilia, a compartmentalized organelle associated with motility, chemo-/mechano-/photosensation, and developmental signaling. These pathways include cilium-based intraflagellar transport (IFT) and poorly understood membrane trafficking events. Defects in ciliary transport contribute to the etiology of human ciliary disease such as Bardet-Biedl syndrome (BBS). In this study, we employ the genetically tractable nematode Caenorhabditis elegans to investigate whether endocytosis genes function in cilium formation and/or the transport of ciliary membrane or ciliary proteins. RESULTS: Here we show that localization of the clathrin light chain, AP-2 clathrin adaptor, dynamin, and RAB-5 endocytic proteins overlaps with a morphologically discrete periciliary membrane compartment associated with sensory cilia. In addition, ciliary transmembrane proteins such as G protein-coupled receptors concentrate at periciliary membranes. Disruption of endocytic gene function causes expansion of ciliary and/or periciliary membranes as well as defects in the ciliary targeting and/or transport dynamics of ciliary transmembrane and IFT proteins. Finally, genetic analyses reveal that the ciliary membrane expansions in dynamin and AP-2 mutants require bbs-8 and rab-8 function and that sensory signaling and endocytic genes may function in a common pathway to regulate ciliary membrane volume. CONCLUSIONS: These data implicate C. elegans endocytosis proteins localized at the ciliary base in regulating ciliary and periciliary membrane volume and suggest that membrane retrieval from these compartments is counterbalanced by BBS-8 and RAB-8-mediated membrane delivery.

Joubert syndrome Arl13b functions at ciliary membranes and stabilizes protein transport in Caenorhabditis elegans.

The small ciliary G protein Arl13b is required for cilium biogenesis and sonic hedgehog signaling and is mutated in patients with Joubert syndrome (JS). In this study, using Caenorhabditis elegans and mammalian cell culture systems, we investigated the poorly understood ciliary and molecular basis of Arl13b function. First, we show that Arl13b/ARL-13 localization is frequently restricted to a proximal ciliary compartment, where it associates with ciliary membranes via palmitoylation modification motifs. Next, we find that loss-of-function C. elegans arl-13 mutants possess defects in cilium morphology and ultrastructure, as well as defects in ciliary protein localization and transport; ciliary transmembrane proteins abnormally accumulate, PKD-2 ciliary abundance is elevated, and anterograde intraflagellar transport (IFT) is destabilized. Finally, we show that arl-13 interacts genetically with other ciliogenic and ciliary transport-associated genes in maintaining cilium structure/morphology and anterograde IFT stability. Together, these data implicate a role for JS-associated Arl13b at ciliary membranes, where it regulates ciliary transmembrane protein localizations and anterograde IFT assembly stability.