Genomic identification and functional analysis of essential genes in Caenorhabditis elegans.

BACKGROUND: Essential genes are required for an organism's viability and their functions can vary greatly, spreading across many pathways. Due to the importance of essential genes, large scale efforts have been undertaken to identify the complete set of essential genes and to understand their function. Studies of genome architecture and organization have found that genes are not randomly disturbed in the genome. RESULTS: Using combined genetic mapping, Illumina sequencing, and bioinformatics analyses, we successfully identified 44 essential genes with 130 lethal mutations in genomic regions of C. elegans of around 7.3 Mb from Chromosome I (left). Of the 44 essential genes, six of which were genes not characterized previously by mutant alleles, let-633/let-638 (B0261.1), let-128 (C53H9.2), let-511 (W09C3.4), let-162 (Y47G6A.18), let-510 (Y47G6A.19), and let-131 (Y71G12B.6). Examine essential genes with Hi-C data shows that essential genes tend to cluster within TAD units rather near TAD boundaries. We have also shown that essential genes in the left half of chromosome I in C. elegans function in enzyme and nucleic acid binding activities during fundamental processes, such as DNA replication, transcription, and translation. From protein-protein interaction networks, essential genes exhibit more protein connectivity than non-essential genes in the genome. Also, many of the essential genes show strong expression in embryos or early larvae stages, indicating that they are important to early development. CONCLUSIONS: Our results confirmed that this work provided a more comprehensive picture of the essential gene and their functional characterization. These genetic resources will offer important tools for further heath and disease research.

Systematic functional analysis of the Ras GTPase family unveils a conserved network required for anterograde protein trafficking.

Phylogeny is often used to compare entire families of genes/proteins. We previously showed that classification of Caenorhabditis elegans Rho GTPases on the basis of their enzymatic properties was significantly different from sequence alignments. To further develop this concept, we have developed an integrated approach to classify C. elegans small GTPases based on functional data comprising affinity for GTP, sub-cellular localization, tissue distribution and silencing impact. This analysis led to establish a novel functional classification for small GTPases. To test the relevance of this classification in mammals, we focused our attention on the human orthologs of small GTPases from a specific group comprising arf-1.2, evl-20, arl-1, Y54E10BR.2, unc-108 and rab-7. We then tested their involvement in protein secretion and membrane traffic in mammalian systems. Using this approach we identify a novel network containing 18 GTPases, and 23 functionally interacting proteins, conserved between C. elegans and mammals, which is involved in membrane traffic and protein secretion.

Cytoplasmic LSM-1 protein regulates stress responses through the insulin/IGF-1 signaling pathway in Caenorhabditis elegans.

Genes coding for members of the Sm-like (LSm) protein family are conserved through evolution from prokaryotes to humans. These proteins have been described as forming homo- or heterocomplexes implicated in a broad range of RNA-related functions. To date, the nuclear LSm2-8 and the cytoplasmic LSm1-7 heteroheptamers are the best characterized complexes in eukaryotes. Through a comprehensive functional study of the LSm family members, we found that lsm-1 and lsm-3 are not essential for C. elegans viability, but their perturbation, by RNAi or mutations, produces defects in development, reproduction, and motility. We further investigated the function of lsm-1, which encodes the distinctive protein of the cytoplasmic complex. RNA-seq analysis of lsm-1 mutants suggests that they have impaired Insulin/IGF-1 signaling (IIS), which is conserved in metazoans and involved in the response to various types of stress through the action of the FOXO transcription factor DAF-16. Further analysis using a DAF-16::GFP reporter indicated that heat stress-induced translocation of DAF-16 to the nuclei is dependent on lsm-1. Consistent with this, we observed that lsm-1 mutants display heightened sensitivity to thermal stress and starvation, while overexpression of lsm-1 has the opposite effect. We also observed that under stress, cytoplasmic LSm proteins aggregate into granules in an LSM-1-dependent manner. Moreover, we found that lsm-1 and lsm-3 are required for other processes regulated by the IIS pathway, such as aging and pathogen resistance.

xnd-1 regulates the global recombination landscape in Caenorhabditis elegans.

Meiotic crossover (CO) recombination establishes physical linkages between homologous chromosomes that are required for their proper segregation into developing gametes, and promotes genetic diversity by shuffling genetic material between parental chromosomes. COs require the formation of double strand breaks (DSBs) to create the substrate for strand exchange. DSBs occur in small intervals called hotspots and significant variation in hotspot usage exists between and among individuals. This variation is thought to reflect differences in sequence identity and chromatin structure, DNA topology and/ or chromosome domain organization. Chromosomes show different frequencies of nondisjunction (NDJ), reflecting inherent differences in meiotic crossover control, yet the underlying basis of these differences remains elusive. Here we show that a novel chromatin factor, X non-disjunction factor 1 (xnd-1), is responsible for the global distribution of COs in C. elegans. xnd-1 is also required for formation of double-strand breaks (DSBs) on the X, but surprisingly XND-1 protein is autosomally enriched. We show that xnd-1 functions independently of genes required for X chromosome-specific gene silencing, revealing a novel pathway that distinguishes the X from autosomes in the germ line, and further show that xnd-1 exerts its effects on COs, at least in part, by modulating levels of H2A lysine 5 acetylation.

Identification of 526 conserved metazoan genetic innovations exposes a new role for cofactor E-like in neuronal microtubule homeostasis.

The evolution of metazoans from their choanoflagellate-like unicellular ancestor coincided with the acquisition of novel biological functions to support a multicellular lifestyle, and eventually, the unique cellular and physiological demands of differentiated cell types such as those forming the nervous, muscle and immune systems. In an effort to understand the molecular underpinnings of such metazoan innovations, we carried out a comparative genomics analysis for genes found exclusively in, and widely conserved across, metazoans. Using this approach, we identified a set of 526 core metazoan-specific genes (the 'metazoanome'), approximately 10% of which are largely uncharacterized, 16% of which are associated with known human disease, and 66% of which are conserved in Trichoplax adhaerens, a basal metazoan lacking neurons and other specialized cell types. Global analyses of previously-characterized core metazoan genes suggest a prevalent property, namely that they act as partially redundant modifiers of ancient eukaryotic pathways. Our data also highlights the importance of exaptation of pre-existing genetic tools during metazoan evolution. Expression studies in C. elegans revealed that many metazoan-specific genes, including tubulin folding cofactor E-like (TBCEL/coel-1), are expressed in neurons. We used C. elegans COEL-1 as a representative to experimentally validate the metazoan-specific character of our dataset. We show that coel-1 disruption results in developmental hypersensitivity to the microtubule drug paclitaxel/taxol, and that overexpression of coel-1 has broad effects during embryonic development and perturbs specialized microtubules in the touch receptor neurons (TRNs). In addition, coel-1 influences the migration, neurite outgrowth and mechanosensory function of the TRNs, and functionally interacts with components of the tubulin acetylation/deacetylation pathway. Together, our findings unveil a conserved molecular toolbox fundamental to metazoan biology that contains a number of neuronally expressed and disease-related genes, and reveal a key role for TBCEL/coel-1 in regulating microtubule function during metazoan development and neuronal differentiation.

RNA-seq analysis of the C. briggsae transcriptome.

Curation of a high-quality gene set is the critical first step in genome research, enabling subsequent analyses such as ortholog assignment, cis-regulatory element finding, and synteny detection. In this project, we have reannotated the genome of Caenorhabditis briggsae, the best studied sister species of the model organism Caenorhabditis elegans. First, we applied a homology-based gene predictor genBlastG to annotate the C. briggsae genome. We then validated and further improved the C. briggsae gene annotation through RNA-seq analysis of the C. briggsae transcriptome, which resulted in the first validated C. briggsae gene set (23,159 genes), among which 7347 genes (33.9% of all genes with introns) have all of their introns confirmed. Most genes (14,812, or 68.3%) have at least one intron validated, compared with only 3.9% in the most recent WormBase release (WS228). Of all introns in the revised gene set (103,083), 61,503 (60.1%) have been confirmed. Additionally, we have identified numerous trans-splicing leaders (SL1 and SL2 variants) in C. briggsae, leading to the first genome-wide annotation of operons in C. briggsae (1105 operons). The majority of the annotated operons (564, or 51.0%) are perfectly conserved in C. elegans, with an additional 345 operons (or 31.2%) somewhat divergent. Additionally, RNA-seq analysis revealed over 10 thousand small-size assembly errors in the current C. briggsae reference genome that can be readily corrected. The revised C. briggsae genome annotation represents a solid platform for comparative genomics analysis and evolutionary studies of Caenorhabditis species.

High-throughput capturing and characterization of mutations in essential genes of Caenorhabditis elegans.

BACKGROUND: Essential genes are critical for the development of all organisms and are associated with many human diseases. These genes have been a difficult category to study prior to the availability of balanced lethal strains. Despite the power of targeted mutagenesis, there are limitations in identifying mutations in essential genes. In this paper, we describe the identification of coding regions for essential genes mutated using forward genetic screens in Caenorhabditis elegans. The lethal mutations described here were isolated and maintained by a wild-type allele on a rescuing duplication. RESULTS: We applied whole genome sequencing to identify the causative molecular lesion resulting in lethality in existing C. elegans mutant strains. These strains are balanced and can be easily maintained for subsequent characterization. Our method can be effectively used to analyze mutations in a large number of essential genes. We describe here the identification of 64 essential genes in a region of chromosome I covered by the duplication sDp2. Of these, 42 are nonsense mutations, six are splice signal mutations, one deletion, and 15 are non-synonymous mutations. Many of the essential genes in this region function in cell cycle, transcriptional regulation, and RNA processing. CONCLUSIONS: The essential genes identified here are represented by mutant strains, many of which have more than one mutant allele. The genetic resource can be utilized to further our understanding of essential gene function and will be applicable to the study of C. elegans development, conserved cellular function, and ultimately lead to improved human health.

Fine tuning of RFX/DAF-19-regulated target gene expression through binding to multiple sites in Caenorhabditis elegans.

In humans, mutations of a growing list of regulatory factor X (RFX) target genes have been associated with devastating genetics disease conditions including ciliopathies. However, mechanisms underlying RFX transcription factors (TFs)-mediated gene expression regulation, especially differential gene expression regulation, are largely unknown. In this study, we explore the functional significance of the co-existence of multiple X-box motifs in regulating differential gene expression in Caenorhabditis elegans. We hypothesize that the effect of multiple X-box motifs is not a simple summation of binding effect to individual X-box motifs located within a same gene. To test this hypothesis, we identified eight C. elegans genes that contain two or more X-box motifs using comparative genomics. We examined one of these genes, F25B4.2, which contains two 15-bp X-box motifs. F25B4.2 expression in ciliated neurons is driven by the proximal motif and its expression is repressed by the distal motif. Our data suggest that two X-box motifs cooperate together to regulate the expression of F25B4.2 in location and intensity. We propose that multiple X-box motifs might be required to tune specific expression level. Our identification of genes with multiple X-box motifs will also improve our understanding of RFX/DAF-19-mediated regulation in C. elegans and in other organisms including humans.

The ortholog of the human proto-oncogene ROS1 is required for epithelial development in C. elegans.

The orphan receptor ROS1 is a human proto-oncogene, mutations of which are found in an increasing number of cancers. Little is known about the role of ROS1, however in vertebrates it has been implicated in promoting differentiation programs in specialized epithelial tissues. In this study we show that the C. elegans ortholog of ROS1, the receptor tyrosine kinase ROL-3, has an essential role in orchestrating the morphogenesis and development of specialized epidermal tissues, highlighting a potentially conserved function in coordinating crosstalk between developing epithelial cells. We also provide evidence of a direct relationship between ROL-3, the mucin SRAP-1, and BCC-1, the homolog of mRNA regulating protein Bicaudal-C. This study answers a longstanding question as to the developmental function of ROL-3, identifies three new genes that are expressed and function in the developing epithelium of C. elegans, and introduces the nematode as a potentially powerful model system for investigating the increasingly important, yet poorly understood, human oncogene ROS1.

Protein phosphatase 2A cooperates with the autophagy-related kinase UNC-51 to regulate axon guidance in Caenorhabditis elegans.

UNC-51 is a serine/threonine protein kinase conserved from yeast to humans. The yeast homolog Atg1 regulates autophagy (catabolic membrane trafficking) required for surviving starvation. In C. elegans, UNC-51 regulates the axon guidance of many neurons by a different mechanism than it and its homologs use for autophagy. UNC-51 regulates the subcellular localization (trafficking) of UNC-5, a receptor for the axon guidance molecule UNC-6/Netrin; however, the molecular details of the role for UNC-51 are largely unknown. Here, we report that UNC-51 physically interacts with LET-92, the catalytic subunit of serine/threonine protein phosphatase 2A (PP2A-C), which plays important roles in many cellular functions. A low allelic dose of LET-92 partially suppressed axon guidance defects of weak, but not severe, unc-51 mutants, and a low allelic dose of PP2A regulatory subunits A (PAA-1/PP2A-A) and B (SUR-6/PP2A-B) partially enhanced the weak unc-51 mutants. We also found that LET-92 can work cell-non-autonomously on axon guidance in neurons, and that LET-92 colocalized with UNC-51 in neurons. In addition, PP2A dephosphorylated phosphoproteins that had been phosphorylated by UNC-51. These results suggest that, by forming a complex, PP2A cooperates with UNC-51 to regulate axon guidance by regulating phosphorylation. This is the first report of a serine/threonine protein phosphatase functioning in axon guidance in vivo.

An antibiotic selection marker for nematode transgenesis.

We have developed a nematode transformation vector carrying the bacterial neomycin resistance gene (NeoR) and shown that it could confer resistance to G-418 on both wild-type Caenorhabditis elegans and C. briggsae. This selection system allows hands-off maintenance and enrichment of transgenic worms carrying non-integrated transgenes on selective plates. We also show that this marker can be used for Mos1-mediated single-copy insertion in wild-type genetic backgrounds (MosSCI-biotic).

Mutations in a member of the Ras superfamily of small GTP-binding proteins causes Bardet-Biedl syndrome.

RAB, ADP-ribosylation factors (ARFs) and ARF-like (ARL) proteins belong to the Ras superfamily of small GTP-binding proteins and are essential for various membrane-associated intracellular trafficking processes. None of the approximately 50 known members of this family are linked to human disease. Using a bioinformatic screen for ciliary genes in combination with mutational analyses, we identified ARL6 as the gene underlying Bardet-Biedl syndrome type 3, a multisystemic disorder characterized by obesity, blindness, polydactyly, renal abnormalities and cognitive impairment. We uncovered four different homozygous substitutions in ARL6 in four unrelated families affected with Bardet-Biedl syndrome, two of which disrupt a threonine residue important for GTP binding and function of several related small GTP-binding proteins. Analysis of the Caenorhabditis elegans ARL6 homolog indicates that it is specifically expressed in ciliated cells, and that, in addition to the postulated cytoplasmic functions of ARL proteins, it undergoes intraflagellar transport. These findings implicate a small GTP-binding protein in ciliary transport and the pathogenesis of a pleiotropic disorder.

The α-arrestin SUP-13/ARRD-15 promotes isoform turnover of actin-interacting protein 1 in Caenorhabditis elegans striated muscle.

Precise arrangement of actin, myosin, and other regulatory components in a sarcomeric pattern is critical for producing contractile forces in striated muscles. Actin-interacting protein 1 (AIP1), also known as WD-repeat protein 1 (WDR1), is one of essential factors that regulate sarcomeric assembly of actin filaments. In the nematode Caenorhabditis elegans, mutation in unc-78, encoding one of the two AIP1 isoforms, causes severe disorganization of sarcomeric actin filaments and near paralysis, but mutation in sup-13 suppresses the unc-78-mutant phenotypes to restore nearly normal sarcomeric actin organization and worm motility. Here, we identified that sup-13 is a nonsense allele of arrd-15 encoding an α-arrestin. The sup-13/arrd-15 mutation suppressed the phenotypes of unc-78 null mutant but required aipl-1 that encodes a second AIP1 isoform. aipl-1 was normally expressed highly in embryos and downregulated in mature muscle. However, in the sup-13/arrd-15 mutant, the AIPL-1 protein was maintained at high levels in adult muscle to compensate for the absence of the UNC-78 protein. The sup-13/arrd-15 mutation caused accumulation of ubiquitinated AIPL-1 protein, suggesting that a normal function of sup-13/arrd-15 is to enhance degradation of ubiquitinated AIPL-1, thereby promoting transition of AIP1 isoforms from AIPL-1 to UNC-78 in developing muscle. These results suggest that α-arrestin is a novel factor to promote isoform turnover by enhancing protein degradation.

The atm-1 gene is required for genome stability in Caenorhabditis elegans.

The Ataxia-telangiectasia-mutated (ATM) gene in humans was identified as the basis of a rare autosomal disorder leading to cancer susceptibility and is now well known as an important signal transducer in response to DNA damage. An approach to understanding the conserved functions of this gene is provided by the model system, Caenorhabditis elegans. In this paper we describe the structure and loss of function phenotype of the ortholog atm-1. Using bioinformatic and molecular analysis we show that the atm-1 gene was previously misannotated. We find that the transcript is in fact a product of three gene predictions, Y48G1BL.2 (atm-1), K10E9.1, and F56C11.4 that together make up the complete coding region of ATM-1. We also characterize animals that are mutant for two available knockout alleles, gk186 and tm5027. As expected, atm-1 mutant animals are sensitive to ionizing radiation. In addition, however, atm-1 mutants also display phenotypes associated with genomic instability, including low brood size, reduced viability and sterility. We document several chromosomal fusions arising from atm-1 mutant animals. This is the first time a mutator phenotype has been described for atm-1 in C. elegans. Finally we demonstrate the use of a balancer system to screen for and capture atm-1-derived mutational events. Our study establishes C. elegans as a model for the study of ATM as a mutator potentially leading to the development of screens to identify therapeutic targets in humans.

Conserved elements associated with ribosomal genes and their trans-splice acceptor sites in Caenorhabditis elegans.

The recent publication of the Caenorhabditis elegans cisRED database has provided an extensive catalog of upstream elements that are conserved between nematode genomes. We have performed a secondary analysis to determine which subsequences of the cisRED motifs are found in multiple locations throughout the C. elegans genome. We used the word-counting motif discovery algorithm DME to form the motifs into groups based on sequence similarity. We then examined the genes associated with each motif group using DAVID and Ontologizer to determine which groups are associated with genes that also have significant functional associations in the Gene Ontology and other gene annotation sources. Of the 3265 motif groups formed, 612 (19%) had significant functional associations with respect to GO terms. Eight of the first 20 motif groups based on frequent dodecamers among the cisRED motif sequences were specifically associated with ribosomal protein genes; two of these were similar to mouse EBP-45, rat HNF3-family and Drosophila Zeste transcription factor binding sites. Additionally, seven motif groups were extensions of the canonical C. elegans trans-splice acceptor site. One motif group was tested for regulatory function in a series of green fluorescent protein expression experiments and was shown to be involved in pharyngeal expression.

An integrated strategy to study muscle development and myofilament structure in Caenorhabditis elegans.

A crucial step in the development of muscle cells in all metazoan animals is the assembly and anchorage of the sarcomere, the essential repeat unit responsible for muscle contraction. In Caenorhabditis elegans, many of the critical proteins involved in this process have been uncovered through mutational screens focusing on uncoordinated movement and embryonic arrest phenotypes. We propose that additional sarcomeric proteins exist for which there is a less severe, or entirely different, mutant phenotype produced in their absence. We have used Serial Analysis of Gene Expression (SAGE) to generate a comprehensive profile of late embryonic muscle gene expression. We generated two replicate long SAGE libraries for sorted embryonic muscle cells, identifying 7,974 protein-coding genes. A refined list of 3,577 genes expressed in muscle cells was compiled from the overlap between our SAGE data and available microarray data. Using the genes in our refined list, we have performed two separate RNA interference (RNAi) screens to identify novel genes that play a role in sarcomere assembly and/or maintenance in either embryonic or adult muscle. To identify muscle defects in embryos, we screened specifically for the Pat embryonic arrest phenotype. To visualize muscle defects in adult animals, we fed dsRNA to worms producing a GFP-tagged myosin protein, thus allowing us to analyze their myofilament organization under gene knockdown conditions using fluorescence microscopy. By eliminating or severely reducing the expression of 3,300 genes using RNAi, we identified 122 genes necessary for proper myofilament organization, 108 of which are genes without a previously characterized role in muscle. Many of the genes affecting sarcomere integrity have human homologs for which little or nothing is known.

Transgenic nematodes as biosensors for metal stress in soil pore water samples.

Caenorhabditis elegans strains carrying stress-reporter green fluorescent protein transgenes were used to explore patterns of response to metals. Multiple stress pathways were induced at high doses by most metals tested, including members of the heat shock, oxidative stress, metallothionein (mtl) and xenobiotic response gene families. A mathematical model (to be published separately) of the gene regulatory circuit controlling mtl production predicted that chemically similar divalent metals (classic inducers) should show additive effects on mtl gene induction, whereas chemically dissimilar metals should show interference. These predictions were verified experimentally; thus cadmium and mercury showed additive effects, whereas ferric iron (a weak inducer) significantly reduced the effect of mercury. We applied a similar battery of tests to diluted samples of soil pore water extracted centrifugally after mixing 20% w/w ultrapure water with air-dried soil from an abandoned lead/zinc mine in the Murcia region of Spain. In addition, metal contents of both soil and soil pore water were determined by ICP-MS, and simplified mixtures of soluble metal salts were tested at equivalent final concentrations. The effects of extracted soil pore water (after tenfold dilution) were closely mimicked by mixtures of its principal component ions, and even by the single most prevalent contaminant (zinc) alone, though other metals modulated its effects both positively and negatively. In general, mixtures containing similar (divalent) metal ions exhibited mainly additive effects, whereas admixture of dissimilar (e.g. trivalent) ions often resulted in interference, reducing overall levels of stress-gene induction. These findings were also consistent with model predictions.

A strawberry notch homolog, let-765/nsh-1, positively regulates lin-3/egf expression to promote RAS-dependent vulval induction in C. elegans.

The specification and patterning of vulval precursor cells (VPCs) in Caenorhabditiselegans is achieved using a conserved EGFR/RAS signaling pathway that is activated by the ligand lin-3/EGF, which is secreted by the neighboring somatic gonad. Previous work has demonstrated that the expression of lin-3 must be tightly regulated to ensure that only three of six equivalent VPCs are induced to differentiate into the mature vulva. Here, we have identified a novel regulator of EGFR/RAS signaling, let-765/nsh-1, that functions upstream of the pathway to promote vulval induction. let-765 encodes a conserved DExD/H box helicase protein and is the C. elegans ortholog of Drosophila strawberry notch. By investigating genetic interactions between let-765 and RAS pathway genes as well as with synthetic multivulva (synMuv) genes, we have demonstrated that let-765 positively regulates the RAS pathway and antagonizes synMuv activity at the level of lin-3/EGF. In support of these proposals, we found that LET-765 is required for producing wild-type levels of lin-3 mRNA. Mutations in let-765 result in pleiotropic phenotypes that imply its function must be required in multiple developmental processes and, together with data presented here, suggest that LET-765 promotes the expression of diverse targets, potentially through interactions with transcriptional activator or repressor complexes.

Evaluation of the therapeutic potential of carbonic anhydrase inhibitors in two animal models of dystrophin deficient muscular dystrophy.

Duchenne Muscular Dystrophy is an inherited muscle degeneration disease for which there is still no efficient treatment. However, compounds active on the disease may already exist among approved drugs but are difficult to identify in the absence of cellular models. We used the Caenorhabditis elegans animal model to screen a collection of 1000 already approved compounds. Two of the most active hits obtained were methazolamide and dichlorphenamide, carbonic anhydrase inhibitors widely used in human therapy. In C. elegans, these drugs were shown to interact with CAH-4, a putative carbonic anhydrase. The therapeutic efficacy of these compounds was further validated in long-term experiments on mdx mice, the mouse model of Duchenne Muscular Dystrophy. Mice were treated for 120 days with food containing methazolamide or dichlorphenamide at two doses each. Musculus tibialis anterior and diaphragm muscles were histologically analyzed and isometric muscle force was measured in M. extensor digitorum longus. Both substances increased the tetanic muscle force in the treated M. extensor digitorum longus muscle group, dichlorphenamide increased the force significantly by 30%, but both drugs failed to increase resistance of muscle fibres to eccentric contractions. Histological analysis revealed a reduction of centrally nucleated fibers in M. tibialis anterior and diaphragm in the treated groups. These studies further demonstrated that a C. elegans-based screen coupled with a mouse model validation strategy can lead to the identification of potential pharmacological agents for rare diseases.

Convergent evolution of RFX transcription factors and ciliary genes predated the origin of metazoans.

BACKGROUND: Intraflagellar transport (IFT) genes, which are critical for the development and function of cilia and flagella in metazoans, are tightly regulated by the Regulatory Factor X (RFX) transcription factors (TFs). However, how and when their evolutionary relationship was established remains unknown. RESULTS: We have identified evidence suggesting that RFX TFs and IFT genes evolved independently and their evolution converged before the first appearance of metazoans. Both ciliary genes and RFX TFs exist in all metazoans as well as some unicellular eukaryotes. However, while RFX TFs and IFT genes are found simultaneously in all sequenced metazoan genomes, RFX TFs do not co-exist with IFT genes in most pre-metazoans and thus do not regulate them in these organisms. For example, neither the budding yeast nor the fission yeast possesses cilia although both have well-defined RFX TFs. Conversely, most unicellular eukaryotes, including the green alga Chlamydomonas reinhardtii, have typical cilia and well conserved IFT genes but lack RFX TFs. Outside of metazoans, RFX TFs and IFT genes co-exist only in choanoflagellates including M. brevicollis, and only one fungus Allomyces macrogynus of the 51 sequenced fungus genomes. M. brevicollis has two putative RFX genes and a full complement of ciliary genes. CONCLUSIONS: The evolution of RFX TFs and IFT genes were independent in pre-metazoans. We propose that their convergence in evolution, or the acquired transcriptional regulation of IFT genes by RFX TFs, played a pivotal role in the establishment of metazoan.