<title>Background</title> <p>The growing availability of the whole-genome sequences of eukaryotes has accelerated large-scale functional studies to understand the mechanisms of animal development and evolution [##REF##16000021##1##, ####REF##12777501##2##, ##REF##11336667##3##, ##REF##12459721##4####12459721##4##]. Many of these studies have highlighted the importance of regulatory evolution and the fundamental role that transcription factors (TFs) play in this process. Alterations in TF function and regulation are linked to phenotypic variation [##REF##11117736##5##, ####REF##17632547##6##, ##REF##14758362##7####14758362##7##] as well as numerous pathologies, including cancers [##REF##17957143##8##,##REF##17661355##9##]. Therefore, a detailed analysis of sequence and function of TFs across animal phyla will provide important information about their evolutionary patterns, thereby increasing our ability to understand the molecular basis of diseases and organismal complexity. The nematode <italic>Caenorhabditis elegans </italic>serves as a powerful model organism to unravel TF function due to the wealth of available resources and the ease with which it can be reared, maintained, and manipulated in the laboratory [##REF##17549065##10##]. The completion of its genome sequence has aided in the design of large-scale experiments that are beginning to elucidate the complexity of transcriptional regulation and gene interaction networks in multicelllular eukaryotes [##REF##12142361##11##,##REF##14704431##12##]. The recent releases of the genome sequence of two other <italic>Caenorhabditid </italic>species, <italic>C. briggsae </italic>[##REF##14624247##13##] and <italic>C. remanei </italic>[##UREF##0##14##], provide an excellent opportunity for genome-wide study of the conservation and evolution of transcription factors across nematodes. These three species are estimated to have shared a common ancestor between 20–120 million years ago [##REF##14624247##13##, ####UREF##0##14##, ##REF##14635896##15####14635896##15##] and while they are morphologically similar, studies have shown differences in development and behavior [##REF##18050493##16##].</p> <p>As a first step in facilitating the comparative study of TFs in nematodes, we have compiled an updated list of putative TF genes in <italic>C. elegans </italic>and used it to identify orthologs in <italic>C. briggsae </italic>and <italic>C. remanei</italic>. Our results show that two-thirds of all <italic>C. elegans </italic>TF genes have 3-way one-to-one best reciprocal orthologs in the other two species, whereas the remaining third are either species-specific paralogs or too divergent to assign proper orthologous relationships. We observed that among <italic>Caenorhabditid </italic>species, although TF genes have a greater sequence divergence than the non-TF genes, they exhibit significantly more detectable interspecific orthologs than non-TF genes. We also identified 150 best reciprocal orthologs of the TF genes conserved among nematodes in fruit fly (<italic>Drosophila melanogaster</italic>), mouse (<italic>Mus musculus</italic>), and human (<italic>Homo sapiens</italic>) many of which are associated with known disorders. We also examined the relationship between gene function and interactions, the results of which demonstrate that conserved TF genes exhibit a significantly greater number of interactions and are more likely to be associated with mutant phenotypes when compared to those that lack detectable orthologs. Our findings provide a framework for future studies of nematode TFs and facilitate the development of resources allowing us to study morphological and developmental diversity in metazoans.</p>
<title>Methods</title> <title><italic>C. elegans</italic>, <italic>C. briggsae </italic>and <italic>C. remanei </italic>TF gene sets</title> <p>The <italic>C. elegans </italic>TF-encoding genes were searched using 8 GO terms (Table ##TAB##0##1##) within WS173 release of Wormbase. The <italic>C. briggsae </italic>and <italic>C. remanei </italic>TFs were identified using the HMMER [##UREF##3##22##,##REF##9918945##65##] and InParanoid programs [##REF##11743721##21##]. The complete genome sequences of each of the three <italic>Caenorhabditid </italic>species were downloaded from WormBase (<italic>C. elegans </italic>release WS173, <italic>C. briggsae </italic>release WS173 and <italic>C. remanei </italic>release 11/29/2005) [##UREF##1##17##]. As the <italic>C. remanei </italic>predicted peptide dataset is known to contain redundant copies of genes due to heterozygosity in the sequenced genome, (E. Schwartz, personal communication) we used the CD-HIT program (version 2007-0131) [##REF##16731699##66##] in order to cluster and remove all additional transcripts that had greater than or equal to 98% sequence similarity to other transcripts at the protein level. The original dataset of 25,948 transcripts was truncated down to 24,267 non-redundant transcripts that were used in further analysis [##REF##18281268##27##].</p> <p>InParanoid was run with default values, using blastall version 2.2.14 with –VT emulation, on all three complete genome predicted peptide datasets in pairwise comparisons. The results were collected and placed into species-specific paralogs, 2- and 3-way best-hit reciprocal ortholog categories using custom PERL scripts. Each category was searched for genes from the <italic>C. elegans </italic>TF set and the number of TFs in each category was identified (Additional files ##SUPPL##5##6## and ##SUPPL##6##7##). HMM alignment-based searches were carried out on the <italic>C. briggsae </italic>and <italic>C. remanei </italic>predicted peptides using previously established techniques [##UREF##3##22##,##REF##11697912##67##]. The HMMER signature files (profiles) of known DNA binding domains were retrieved from Pfam [##UREF##6##68##]. In most cases, a cut-off score of 0.1 was used. If a HMMER predicted TF gene in non-<italic>elegans </italic>species lacked a homolog in <italic>C. elegans</italic>, it was considered false positive and therefore removed creating the final, conservative datasets that were used in the study.</p> <p>The <italic>C. elegans </italic>orthologs of <italic>D. melanogaster</italic>, <italic>M. musculus </italic>and <italic>H. sapiens </italic>TFs were retrieved using the data available on the InParanoid database [##REF##15608241##23##,##UREF##7##69##].</p> <title>Identification of the TF gene families</title> <p>Genes were grouped into different families based on the presence of known DNA-binding domains according to the WormBase [##UREF##0##14##], Pfam [##UREF##6##68##], and InterPro [##UREF##8##70##] databases. Only well defined and unambiguous domains that are known to be involved in transcriptional regulation were considered. Families with fewer than 5 members were placed together in a miscellaneous category. The TF families shown in Figure ##FIG##2##3## are as follows. AP2: Activator protein-2 family; AT hook: AT hook DNA binding motif (preference to A/T rich region) family; bHLH: basic helix-loop-helix family; bZIP: basic leucine zipper family; CBFB/NF-YA: CCAAT binding factor family; CSD: Cold shock DNA binding domain family; HMG box: High mobility group box family; HOX: Homeobox family; MADF: Myb DNA binding domain family; SAND: DNA binding domain family named after Sp100, AIRE-1, NucP41/75, DEAF-1; SANT: Myb-like DNA binding domain; SMAD: SMAD (Mothers against decapentaplegic (MAD) homolog) domain family; T-box: T-box family; WH: Winged-helix family; WH-FH: Winged-helix and Forkhead domain family; WH-ETS: Winged-helix and ETS domain family; ZF-C2H2: C2H2-type zinc finger protein family; ZF-C2H2-BED: C2H2 and BED-type zinc finger protein family; ZF-BED: BED-type zinc finger family; ZF-C2H2-RING: C2H2 and RING-type zinc finger protein family; ZF-C4/NHR: C4-type zinc finger/Nuclear hormone receptor family; ZF-CCCH: C-x8-C-x5-C-x3-H class of zinc finger family; ZF-DHHC: DHHC-type zinc finger family; ZF-FLYWCH: FLYWCH-type of zinc finger family; ZF-GATA: GATA class of zinc finger family; ZF-PHD: C4HC3 zinc-finger-like motif family; ZF-others: zinc finger family members not listed above; ZF-DM: DM (dsx and mab-3) zinc finger family; ZF, AT hook: AT hook and zinc finger domain family; ZF, SANT: SANT and zinc finger domain family; Misc: Miscellaneous TF family not listed above.</p> <title>Generation of the chromosomal map</title> <p>The physical locations of <italic>C. elegans </italic>and <italic>C. briggsae </italic>TF and non-TF genes were retrieved from Wormbase (WS173 release) and grouped into non-overlapping windows of 200 kb (similar to the 250 kb used by [##REF##15268860##33##]). A 400 kb window analysis was also performed and the conclusions remain the same (data not shown). Since many genes are alternatively spliced, we eliminated transcript-specific bias by focusing on single open reading frame for each transcription factor. In the case of <italic>C. briggsae</italic>, a total of 1329 genes were not assigned to any of the chromosomes and hence were excluded from the analysis. For simplicity, we only used the average between the start and end positions as a proxy for the gene position. The significance of TF clustering on chromosomes was determined by comparing their frequency with the overall frequency of genes in a given window using a χ<sup>2 </sup>test [##REF##15268860##33##]. Clusters with <italic>p </italic>value less than 0.05 were considered significant.</p> <title>Phylogenetic analysis of the nematode NHR genes</title> <p>The predicted <italic>C. elegans </italic>NHR gene dataset (283 genes) was used to identify orthologs and paralogs in <italic>C. briggsae </italic>and <italic>C. remanei </italic>using the complete genome INPARANOID datasets (see above). 204 and 152 potential homologs were identified in <italic>C. briggsae </italic>and <italic>C. remanei</italic>, respectively. The peptide dataset was aligned using Dialign 2.2 [##REF##10222408##71##] and then manually inspected. We identified two large conserved blocks within most predicted peptides and removed all sequences that did not align within these blocks. The remaining sequences were then realigned with Dialign 2.2 and truncated only to retain the two conserved domains. As per Robinson-Rechavi et al. [##REF##15983867##29##] we chose to use only ungapped sites and removed first sequences missing significant portions of the conserved domains and finally excluded all gapped sites. In the end, we retained 437 sequences (213 <italic>C. elegans</italic>, 106 <italic>C. briggsae </italic>and 118 <italic>C. remanei</italic>) for phylogenetic analysis.</p> <p>The phylogeny was constructed using a maximum likelihood based method as implemented in PhyML [##REF##14530136##72##] using the JTT substitution model [##REF##1633570##73##] with the default proportion of invariable sites (0.0) and rate heterogeneity between sites corrected by a gamma law (using the default gamma parameter of 1.0 and eight rate categories). The phylogeny was then bootstrapped by generating 1000 randomized datasets using SEQBOOT and assessing the percentage of consensus trees using CONSENSE, both in the PHYLIP package [##UREF##9##74##].</p> <title>Calculation of TF divergence</title> <p>DNA sequences from <italic>C. elegans</italic>, <italic>C. briggsae </italic>and <italic>C. remanei </italic>were aligned according to their protein alignment using Dialign 2.2 [##REF##15215344##75##] and RevTrans 1.4 [##REF##12824361##76##]. Rates of synonymous substitutions per synonymous site (<italic>d</italic><sub><italic>S</italic></sub>) and non-synonymous substitutions per non-synonymous site (<italic>d</italic><sub><italic>N</italic></sub>) were estimated using codeml from PAML [##REF##12032247##77##]. Evolutionary rates between TF and non-TF data sets were compared using a permuted Kruskal-Wallis rank sum test using 10,000 permutations.</p> <title>Curation of the mutant phenotypes of TFs</title> <p>The RNAi phenotypes of all known <italic>C. elegans </italic>genes were retrieved from Wormbase (WS170 release). A total of 13,648 phenotypes associated with 4,351 genes were analyzed and sorted into 82 different categories (Unc, Dpy, Vul etc.) (Additional files ##SUPPL##12##13## and ##SUPPL##13##14##).</p> <p>For phenotypes associated with <italic>C. elegans </italic>TF orthologs in fly, mouse, and human, we searched Flybase [##UREF##10##78##], NCBI OMIM [##UREF##11##79##], PubMed [##UREF##12##80##], and other public databases (<ext-link ext-link-type="uri" xlink:href="http://www.informatics.jax.org"/>, <ext-link ext-link-type="uri" xlink:href="http://www.bioscience.org/knockout/alphabet.htm"/>, <ext-link ext-link-type="uri" xlink:href="http://www.dsi.univ-paris5.fr/genatlas"/>, <ext-link ext-link-type="uri" xlink:href="http://www.genetests.org"/>). Only those phenotypes that were unambiguous and did not show discrepancy between different published sources were included. In order to reduced any effect linked to a differential amount of genes annotated as involved in particular mutant phenotypes, all the analyses were performed within each phenotypic class by comparing the distribution of genes with mutant phenotypes among the different sets (non-TF genes, TF genes, <italic>C. elegans </italic>TF genes, TF genes conserved among the three nematode species, and TF genes with orthologs in nematodes, fly, mouse, and human).</p> <title>Construction of TF interaction network</title> <p>The <italic>C. elegans </italic>gene network was built using the genetic and protein-protein interaction data for transcription factors curated by BioGRID (version 2.0.27 release) [##UREF##4##37##,##REF##16381927##38##]. The network was visualized by using Cytoscape [##REF##14597658##81##].</p>
<title>Results</title> <title>The <italic>C. elegans </italic>TF gene set</title> <p>As a first step in the identification of TFs in <italic>Caenorhabditid </italic>species, we generated an updated list of putative <italic>C. elegans </italic>TF genes by searching its annotated genome sequence (Wormbase WS173 release) [##UREF##1##17##] for gene ontology (GO) terms associated with transcription factors. This led to the identification of 1271 putative TF genes (Table ##TAB##0##1##). Since our criteria for selecting a TF was the presence of a well-defined DNA binding domain that selectively modulates gene transcription (for example, bHLH or homeobox), we manually inspected the above list of putative TFs. This allowed us to reject 564 genes as false positives since these encode factors that are associated with the basal transcriptional apparatus (for example, DNA polymerases), chromatin alterations, DNA packaging (histones), as well as entries that were incorrectly curated in Wormbase (Additional files ##SUPPL##0##1##, ##SUPPL##1##2##, ##SUPPL##2##3##). To the remaining genes (707), we added 281 TF encoding genes found in published literature and other public database entries that were not identified in our initial search (See Materials and Methods). The final <italic>C. elegans </italic>TF set included a total of 988 genes (Table ##TAB##1##2## and additional file ##SUPPL##3##4##), of which 917 are shared with the previously annotated <italic>C. elegans </italic>TF set (wTF2.0, 934 genes) [##REF##16420670##18##]. The 17 genes in the wTF2.0 set that are not shared in our updated set either lack a known DNA binding domain or are annotated as pseudogenes (Additional file ##SUPPL##4##5##). The increased number of genes in the present TF set likely results from the availability of annotations published since the compilation of wTF2.0.</p> <title>Identification of transcription factors in nematodes and other phyla</title> <p>We used the newly defined <italic>C. elegans </italic>TF set to search for homologs in the fully sequenced genomes of <italic>C. briggsae </italic>(CB3 release) and <italic>C. remanei </italic>(11/29/2005 release) [##UREF##1##17##,##REF##17608563##19##,##UREF##2##20##]. We used InParanoid [##REF##11743721##21##] to identify 713 and 703 best reciprocal hit orthologs in <italic>C. briggsae </italic>and <italic>C. remanei</italic>, respectively (Table ##TAB##1##2##, see Material and Methods). To these lists, we added 282 <italic>C. briggsae </italic>and 390 <italic>C. remanei </italic>putative TF genes that were identified through Hidden Markov Model (HMM)-based searches [##UREF##3##22##]. Altogether, a total of 995 and 1093 potential TF encoding genes were identified in <italic>C. briggsae </italic>and <italic>C. remanei</italic>, respectively (Table ##TAB##1##2## and additional files ##SUPPL##5##6## and ##SUPPL##6##7##). Among the TF orthologs in the three nematode species, we identified 652 genes that exhibit a 3-way best reciprocal BLAST orthologous relationship (hereafter referred to as the nematode TF set) (Figure ##FIG##0##1##). The proportion of <italic>C. elegans </italic>TF genes with detectable orthologs in <italic>C. briggsae </italic>(713/995, 71.7%) is significantly higher compared to the proportion of all conserved genes between the two species (12858/20621, 62.4%) [##REF##14624247##13##] (χ<sup>2 </sup>= 7.56, df = 1, <italic>p </italic>= 6.0 × 10<sup>-3</sup>), which may indicate strong selective pressure to maintain these genes.</p> <p>To examine the evolutionary conservation of the nematode TF set of genes in other phyla, we searched for their orthologs in the genomes of fruit fly(<italic>D. melanogaster</italic>), mouse (<italic>M. musculus</italic>), and human (<italic>H. sapiens</italic>). Using the InParanoid database [##REF##15608241##23##] we identified a total of 150 TFs that exhibit reciprocal orthologous relationships between three nematode species and are conserved in fly, mouse, and human (Additional files ##SUPPL##3##4## and ##SUPPL##7##8##).</p> <title>Coding sequence divergence in nematode TF genes</title> <p>Best-hit reciprocal orthologs could not be identified for 215 TF genes in <italic>C. elegans</italic>, 211 in <italic>C. briggsae</italic>, and 310 in <italic>C. remanei </italic>(Figure ##FIG##0##1## and additional files ##SUPPL##3##4##, ##SUPPL##5##6##, and ##SUPPL##6##7##). It should be pointed out that <italic>C. briggsae </italic>and <italic>C. remanei </italic>TF genes are based on computational predictions and that the <italic>C. remanei </italic>genome has yet to be assembled; hence while many of the TF genes without detectable orthologs may have arisen by lineage-specific gene duplication, others could result from incomplete annotation of the <italic>C. briggsae </italic>and <italic>C. remanei </italic>genomes. Therefore, the actual number of divergent TF genes in these species is likely to be smaller than the numbers we have estimated. To further study this set of genes in <italic>C. briggsae </italic>(211), we searched for their closest homologs in <italic>C. elegans</italic>. This revealed 30 genes with weak sequence similarity (BLASTP E-value &gt; 10<sup>-10</sup>) suggesting that these most likely represent candidate <italic>C. briggsae</italic>-specific TF genes (Additional file ##SUPPL##8##9##). The remaining 181 appear to be species-specific paralogs, of which 69 are zinc finger-C4/nuclear hormone receptor (ZF-C4/NHR) family members (see below).</p> <p>Previous studies in humans and other organisms have shown that TF genes tend to evolve more rapidly than non-transcription factor (non-TF) genes [##REF##16525476##24##, ####REF##16237444##25##, ##REF##17994087##26####17994087##26##], therefore we performed a similar analysis in nematodes by analyzing their coding sequence divergence and comparing it to non-TF genes. Due to the large divergence times between the three species [##REF##14624247##13##,##REF##14635896##15##], the rate of synonymous substitution per synonymous site (<italic>d</italic><sub>S</sub>) for many genes is likely to be saturated (<italic>d</italic><sub>S </sub>&gt; 3, Figure ##FIG##1##2A##). Therefore we restricted our analysis to the rate of non-synonymous substitution per non-synonymous site (<italic>d</italic><sub>N</sub>), which does not show such saturation [##REF##18281268##27##]. We found a significantly higher <italic>d</italic><sub>N </sub>for TF genes conserved among nematodes (652) when compared to 3-way conserved non-TF gene orthologs (10,827 genes; Kruskal-Wallis rank sum test, <italic>p </italic>&lt; 2.2 × 10<sup>-16</sup>; Figure ##FIG##1##2B##), whereas no difference was detected between TF genes with orthologs in nematodes, fly, mouse, and human (150) and non-TF genes (Kruskal-Wallis test, <italic>p </italic>= 0.3498).</p> <title>Distribution of TF families in <italic>C. elegans</italic></title> <p>We studied the distribution of protein families among <italic>C. elegans </italic>TFs based on known DNA binding domains. This analysis revealed more than 50 distinct families of which 30 were found to contain 5 or more members (Figure ##FIG##2##3##). No significant difference was observed in the representation of various families between the <italic>C. elegans </italic>set and the set conserved among three nematodes (χ<sup>2 </sup>= 35.05, df = 30, <italic>p </italic>= 0.2408), indicating that the distributions of TF families in these species may be similar. The majority of genes in <italic>C. elegans </italic>and nematode TF sets (28.6% and 20.5%, respectively) were found to encode the nuclear hormone receptors (NHRs), a C4-type sub-family of zinc finger proteins that play key roles in development and homeostasis [##REF##11427696##28##]. The NHR genes was previously shown to have undergone extensive lineage-specific expansion in <italic>C. elegans </italic>[##REF##15983867##29##]. Besides NHR, HOX genes that regulate cell fate specification and embryogenesis [##REF##16341070##30##] are also among the largest TF families in nematodes (11% of <italic>C. elegans </italic>TF genes, and 11.6% of nematode TF genes) (Figure ##FIG##2##3A, C##). In contrast, the distribution of TF families among the divergent <italic>C. elegans </italic>gene set (215 genes) differs significantly from that observed among the entire <italic>C. elegans </italic>TF set (χ<sup>2 </sup>= 83.91, df = 30, <italic>p </italic>= 5.33 × 10<sup>-7</sup>) due to its high proportion of NHR genes (52.6%, Figure ##FIG##2##3B##). Likewise, the representation of different families among TF genes with orthologs in nematodes, fly, mouse, and human also differs from that of the <italic>C. elegans </italic>TF set (χ<sup>2 </sup>= 152.27, df = 30, <italic>p </italic>= 0, Figure ##FIG##2##3D##). Interestingly, the single largest conserved family represented among the orthologs in three different phyla is HOX (17.3%), supporting multiple previous studies indicating the importance of this family among all metazoans [##REF##16420670##18##,##REF##16341070##30##,##REF##17244357##31##].</p> <title>Chromosomal distribution of TF genes in <italic>C. elegans </italic>and <italic>C. briggsae</italic></title> <p>Studies in <italic>C. elegans </italic>as well as other organisms have shown that genes that are co-expressed and/or functionally related are frequently clustered together on chromosomes [##REF##12214599##32##, ####REF##15268860##33##, ##REF##12478293##34##, ##REF##12144710##35##, ##REF##11992122##36####11992122##36##]. To investigate whether TF genes in nematodes exhibit a similar pattern, we plotted the physical locations of <italic>C. elegans </italic>and <italic>C. briggsae </italic>TF genes using non-overlapping windows of 200 kb (the genome of <italic>C. remanei </italic>has not yet been assembled and therefore was not used in this analysis). Figure ##FIG##3##4## shows that TF genes in <italic>C. elegans </italic>and <italic>C. briggsae </italic>as well as those that are conserved among nematodes are non-randomly distributed on chromosomes. A total of 183 <italic>C. elegans </italic>TF genes were found to be located in 25 distinct clusters (marked with stars in Figure ##FIG##3##4A##, Table ##TAB##2##3##). A similar pattern was observed in <italic>C. briggsae </italic>(184 genes in 27 clusters) (Figure ##FIG##3##4B## and Table ##TAB##2##3##). Chromosome V carries highest number of clusters (and genes) in both species (<italic>C. elegans</italic>: 97 genes in 10 clusters; <italic>C. briggsae</italic>: 64 genes in 7 clusters) (Table ##TAB##2##3##) that are primarily composed of NHR family members (92% in <italic>C. elegans </italic>and 84% in <italic>C. briggsae</italic>, red bars in Figure ##FIG##3##4##).</p> <p>The analysis of the chromosomal distribution of TF genes also revealed that many members of the large TF families, such as ZF-C4/NHR, ZF-C2H2, T-box and HOX are arranged in perfect tandem arrays (defined as having a contiguous repetition of TF genes) (267 <italic>C. elegans </italic>genes in 103 arrays, 235 <italic>C. briggsae </italic>genes in 107 arrays), the largest of which consists of 8 NHR genes in <italic>C. elegans </italic>and 6 NHR genes in <italic>C. briggsae </italic>(Figure ##FIG##4##5## and additional files ##SUPPL##9##10## and ##SUPPL##10##11##). Although the majority of such arrays consists of genes of the same TF family (76.7% in <italic>C. elegans </italic>and 67.3% in <italic>C. briggsae</italic>), less than half of all genes found in such arrays have best-reciprocal hit orthologs between the two species (31.6% in <italic>C. elegans</italic>, 47.2% in <italic>C. briggsae</italic>) (Additional files ##SUPPL##9##10## and ##SUPPL##10##11##) suggesting significant lineage-specific duplication and expansion of the tandem arrays.</p> <title>Evolution of the Nuclear Hormone Receptor family in nematodes</title> <p>Our findings extend Robinson-Rechavi et al.'s analysis of the extensive lineage-specific expansion of NHR genes in <italic>C. elegans </italic>[##REF##15983867##29##] to the other two <italic>Caenorhabditid </italic>species (283, 232, and 256 NHR genes in <italic>C. elegans</italic>, <italic>C. briggsae</italic>, and <italic>C. remanei</italic>, respectively). The sequence analyses revealed a total of 134 NHR genes having 3-way best-reciprocal orthologs among the nematode species (Additional files ##SUPPL##3##4##, ##SUPPL##5##6##, and ##SUPPL##6##7##). The remaining NHRs are composed of what appear to be lineage-specific paralogs and those that have diverged sufficiently in sequence such that orthologous relationships could no longer be assigned.</p> <p>We constructed a phylogenetic tree of the nematode NHR family members (437 genes, see Materials and Methods) to study their inter– as well as intra-specific relationships. The most striking feature of the phylogeny is the frequent presence of several closely related NHRs located tandemly on the same chromosome (Additional file ##SUPPL##11##12##). Such groupings suggest the presence of extensive tandem duplications, which could explain the mechanism behind the expansion of the NHR gene family, and perhaps the independent occurrence of some NHR genes in the lineages of each of these species. In the case of <italic>C. elegans </italic>NHRs, we found at least 15 distinct groups on chromosome V including 7 that are located in one large cluster of the phylogeny (Additional file ##SUPPL##11##12##).</p> <p>The presence of NHRs in chromosomal clusters prompted us to study their distribution in further detail. We identified a total of 47 tandem arrays composed of contiguous repetitions of NHR genes in <italic>C. elegans</italic>, which are found on all chromosomes with the exception of chromosome III (Additional file ##SUPPL##9##10##). These include 10 arrays that are comprised of 5 or more genes, all of which are located on chromosome V. A similar analysis in <italic>C. briggsae </italic>identified 30 NHR arrays having 6 or fewer genes (Additional file ##SUPPL##10##11##). In total, 9 NHR arrays were partially or completely conserved between <italic>C. elegans </italic>and <italic>C. briggsae</italic>. One of these arrays, for instance, consists of 7 genes in <italic>C. elegans </italic>(<italic>nhr-136, nhr-153, nhr-154, nhr-206, nhr-207, nhr-208</italic>, and <italic>nhr-209</italic>) and the corresponding 4 in <italic>C. briggsae </italic>(<italic>CBG23383/Cbr-nhr-136, CBG23380/Cbr-nhr-153, CBG23380/Cbr-nhr-154 </italic>and <italic>CBG23379/Cbr-nhr-209</italic>). This suggests that either the array has expanded in <italic>C. elegans </italic>or perhaps lost 3 of the genes in <italic>C. briggsae</italic>. Examination of the <italic>C. remanei </italic>TFs revealed the presence of best reciprocal hit orthologs for all array members found in <italic>C. elegans </italic>with the exception of <italic>nhr-206 </italic>leading us to propose that <italic>nhr-207 </italic>and <italic>nhr-208 </italic>were most likely lost in the <italic>C. briggsae </italic>lineage. This analysis, however, carries a caveat in that the annotations of the <italic>C. briggsae </italic>and <italic>C. remanei </italic>genomes are based on computational predictions and lack experimental validation.</p> <p>Finally, we found that 7 tandem arrays in <italic>C. briggsae </italic>are composed of NHR genes that lack best reciprocal hit orthologs in <italic>C. elegans </italic>and <italic>C. remanei </italic>(Additional file ##SUPPL##10##11##). The largest of these is comprised of 6 NHR genes (<italic>CBG01243, CBG01244, CBG01245, CBG01246, CBG01247, CBG01248</italic>) (Figure ##FIG##4##5##). These <italic>C. briggsae</italic>-specific arrays may be caused by lineage-specific expansion although the possibility of a selective loss of their orthologs in other species cannot be ruled out.</p> <title>Comparison of TF gene sequence conservation and function in <italic>C. elegans</italic></title> <p>We investigated the relationship between sequence conservation and function of TF genes in <italic>C. elegans</italic>. From a comprehensive list of 13,647 RNAi phenotypes associated with 4,351 genes [##UREF##0##14##], we identified 281 TFs that exhibit one or more mutant phenotypes (Additional file ##SUPPL##12##13##). These consist of more than half of all TF genes conserved among nematodes, fly, mouse, and human (52.7%, 79 of 150), over one-third of genes conserved among the three nematode species (36.5%, 238 of 652), and one-fifth of the TF genes in <italic>C. elegans </italic>that did not have identifiable orthologs in the other nematode species (20%, 43 of 215). We also determined the number of distinct mutant phenotypes associated with TF genes in each of the above three groups as well as with non-TF genes. This analysis revealed that TF genes conserved among nematodes, fly, mouse, and human are linked to a significantly greater number of mutant phenotypes in <italic>C. elegans </italic>when compared to the other sets (4.38 ± 2.31, 3.36 ± 2.09, 2.91 ± 1.82 and 3.19 ± 1.84 phenotypes per gene for TF genes conserved across phyla, conserved in nematodes, <italic>C. elegans </italic>TF genes without detectable orthologs in the other nematode species and non-TF genes, respectively; Kruskal-Wallis rank sum test, <italic>p </italic>= 8.58 × 10<sup>-3</sup>, 1.8 × 10<sup>-3</sup>, 1.32 × 10<sup>-15</sup>, respectively, after Bonferroni correction). No difference was found in pairwise comparisons between the other gene sets (Kruskal-Wallis rank sum test <italic>p </italic>= 1, in all comparisons after Bonferroni correction).</p> <p>To further analyze the roles of <italic>C. elegans </italic>TF genes in specific tissues and developmental processes, RNAi phenotypes were sorted into six broad categories: viability (embryonic and post-embryonic growth and survival), fertility (germline and germ cells), sex (sex determination and reproductive system), vulva (vulval cell proliferation and morphogenesis), body (cuticle, size, and morphology), and behavior (movement and feeding) (Additional files ##SUPPL##12##13## and ##SUPPL##13##14##). Among the six categories, "viability" ranks highest in terms of the proportion of TF and non-TF genes (Figure ##FIG##5##6##). However, it is important to keep in mind that this may be linked to a greater interest in identifying transcription factors that are involved in growth and survival of <italic>C. elegans</italic>. A further sub-classification of this category into "embryonic viability" and "post-embryonic viability" (based on the phenotype when lethality occurs in RNAi-treated animals) revealed that among the "embryonic viability" class TF genes conserved among nematode species are significantly under-represented when compared to the non-TFs (χ<sup>2 </sup>= 7.39, df = 1, <italic>p </italic>= 6.56 × 10<sup>-3</sup>) (Figure ##FIG##5##6##), while no difference was observed among the datasets for genes affecting post-embryonic viability (χ<sup>2 </sup>= 0.77, df = 1, <italic>p </italic>= 0.38). By contrast, the TF genes conserved among nematodes, fly, mouse, and human showed no enrichment for any of these two categories (χ<sup>2 </sup>= 3.59 and 0.39, df = 1, <italic>p </italic>= 0.059 and 0.53, respectively). Among other categories, we observed an over-representation of mutant phenotypes for nematode-conserved as well as nematode-fly-mouse-human-conserved TF gene sets associated with vulval development (χ<sup>2 </sup>= 8.24 and 11.75, df = 1, <italic>p </italic>= 4.1 × 10<sup>-3 </sup>and 6.08 × 10<sup>-4</sup>, respectively) and sex determination and reproductive system-related processes (χ<sup>2 </sup>= 9.51 and 8.78, df = 1, <italic>p </italic>= 2.04 × 10<sup>-3 </sup>and 3.59 × 10<sup>-3</sup>, respectively) when compared to the non-TF gene set.</p> <title>Phenotypes associated with nematode TF orthologs in fly, mouse, and human</title> <p>The above findings that more than half of <italic>C. elegans </italic>TF genes conserved across phyla are associated with RNAi phenotypes prompted us to examine their mutant phenotypes in other organisms. We found that 69 (46%) of TF genes conserved among nematodes, fly, mouse, and human are associated with lethal phenotype in fly (Additional file ##SUPPL##14##15##). In the case of mouse, out of a total of 81 orthologs for which knock out and mutant phenotypes are described (see Materials and Methods), 75 (92.6%) exhibit defects ranging from mild to gross abnormalities, including lethality (Additional file ##SUPPL##14##15##). A similar analysis in human revealed 35 TF genes linked to various diseases and genetic disorders (Table ##TAB##3##4##). In total, 44 (29.3%) TF genes regulating <italic>C. elegans </italic>development and behavior are also essential either in mouse or human or both. These include 30 genes that control viability in the fruit fly. Overall, 121 (80.7%) TF genes conserved among nematodes, fly, mouse, and human play important roles in at least one of these organisms. This is likely an underestimate due to technical limitations of RNAi experiments (e.g., strains, redundancy of factors or pathways) and that comparisons between organisms involve different experimental approaches (e.g., RNAi in <italic>C elegans </italic>and chromosomal mutations in <italic>D. melanogaster</italic>). Thus, functional studies of conserved TFs in <italic>C. elegans </italic>promise to elucidate mechanisms involved in biological processes conserved across phyla.</p> <title>Analysis of TF interaction networks in <italic>C. elegans</italic></title> <p>To further explore the mechanism of transcription factor function in metazoans, we generated an interaction map of <italic>C. elegans </italic>TF genes based on known physical and genetic interactions [##UREF##4##37##,##REF##16381927##38##]. The map consists of 1594 interactions involving 277 TF genes and their direct non-TF interactors (Figure ##FIG##6##7A## and additional file ##SUPPL##15##16##). The network appears to be scale free as seen by the presence of several nodes with high degree of connectivity (such as <italic>lin-35</italic>, which shows the highest number of interactions and is connected to more than one-third of all existing nodes; 521 of 1340) (Figure ##FIG##6##7B##). <italic>lin-35 </italic>is an ortholog of the human <italic>Retinoblastoma </italic>(<italic>Rb</italic>) gene which plays an important role in cell proliferation [##REF##9875852##39##,##REF##10502774##40##]. Among the <italic>lin-35 </italic>interacting genes, 43 (8%) encode TFs, of which 18 have best reciprocal hit orthologs in mouse and human. Other prominent hubs include <italic>pal-1 </italic>(conserved among nematode species), as well as other TF genes with orthologs in nematodes, fly, mouse, and human: <italic>tag-331, eya-1</italic>, and <italic>sma-4 </italic>(Figure ##FIG##6##7B##). Each of these genes plays important role in <italic>C. elegans </italic>development, and RNAi-mediated knock-downs cause defects such as slow growth (<italic>pal-1</italic>), lethality (<italic>pal-1, tag-331</italic>), larval arrest (<italic>eya-1, tag-331</italic>), uncoordinated movement (<italic>eya-1</italic>), and small size (<italic>sma-4</italic>) [##REF##16154558##41##, ####UREF##5##42##, ##REF##15791247##43##, ##REF##12529635##44####12529635##44##]. Interestingly, the subnetwork comprising of the hub gene <italic>tag-331 </italic>(human ortholog RNF113A) and its 32 direct interactors appears to be largely isolated. A closer examination revealed that two-thirds of these 22 genes are conserved in nematodes yet lack best reciprocal hit orthologs in fly, mouse, or human genomes. The remaining third includes four genes conserved in nematodes, fly, mouse, and human (<italic>zfp-1, R11F4.1, apl-1 </italic>and <italic>fcd-2</italic>) and whose human homologs are linked to disorders (AF10/MLLT10: leukemia, Glycerol kinase: hyperglycerolemia, APP: Alzheimer's, and FANCD2: Fanconi anemia). It remains to be determined if the human genes interact with RNF113A as well as whether RNF113A is involved in any of these diseases. Among the remaining hub genes, the <italic>eya-1 </italic>mammalian orthologs promote development of tissues and organs [##REF##16154558##41##,##REF##10655545##45##,##REF##17098221##46##] whereas the <italic>sma-4 </italic>ortholog SMAD4/DCP4 acts as a tumor suppressor [##REF##8553070##47##, ####REF##10340381##48##, ##REF##11791187##49####11791187##49##].</p> <p>In addition to analyzing the prominent hubs in the interaction network, we also examined the relationship between connectivity of TFs, sequence conservation, and known function. The results revealed a significantly greater number of interactions among <italic>C. elegans </italic>TFs that are conserved in nematodes, fly, mouse, and human, as compared to those that are not (Kruskal-Wallis rank sum test, <italic>p </italic>= 0.0207, after Bonferroni correction). We also found that TF genes associated with mutant phenotypes in RNAi assays exhibit significantly more interactions when compared to those that lack a detectable phenotype (Kruskal-Wallis rank sum test, <italic>p </italic>= 0.0069). These results are consistent with previous studies showing that highly connected hubs tend to be enriched in essential genes [##REF##16751849##50##,##REF##15616139##51##].</p>
<title>Discussion</title> <p>This paper presents the first genome-wide comparative study of TF genes in nematodes and their orthologs in fly (<italic>D. melanogaster</italic>), mouse (<italic>M. musculus</italic>), and human (<italic>H. sapiens</italic>). We took both computational and manual curation approaches to compile sets of TF genes in three <italic>Caenorhabditid </italic>species, leading to the identification of 988 genes in <italic>C. elegans</italic>, 995 in <italic>C. briggsae </italic>and 1093 in <italic>C. remanei</italic>. A comparison of these data sets has revealed 652 3-way best reciprocal orthologs among these species. Furthermore, using currently available genome annotations, we identified 150 TF gene orthologs shared among nematodes, fly, mouse, and human and shown that according to mutant phenotypes or associated disorders, many of these genes are functionally important. It should be noted that many of the TF genes identified in <italic>C. elegans </italic>as well as most of those identified as orthologs, paralogs, and divergent in the other two nematode species are based entirely on computational predictions, and thus await experimental validation. However, the results of our study suggest the most likely group of candidate genes from which further experimental tests of TF activity can be designed. In contrast, the majority of the orthologs identified in the two other phyla are annotated as TF genes themselves, owing to the extensive experimental validation performed in these organisms.</p> <p>The sequence comparison of orthologs among nematodes has revealed that TF genes conserved among the three nematodes species (652 genes) are evolving more rapidly than non-TF genes, which is in agreement with earlier reports from other species in which TF genes have been shown to be evolving more rapidly than the coding genome average, and that significantly more TF genes have been found to be evolving under positive selection when compared to the rest of the genome [##REF##16525476##24##,##REF##16237444##25##,##REF##17665086##52##,##REF##10022975##53##]. While our observation of a greater number of conserved orthologs among all three nematode species, coupled to an accelerated rate of divergence may seem paradoxical, it may be suggestive of widespread positive selection, and thus divergence, acting on genes that are otherwise functionally important. Given the wide estimates of the divergence time between the three nematode species considered in this study, it is unsurprising that the rate of synonymous substitution (<italic>d</italic><sub>S</sub>) is saturated, and is therefore not amenable for use in analyses that could test the hypothesis of widespread positive selection among TF genes. Additional data, such as a large-scale polymorphism analysis among multiple <italic>Caenorabditid </italic>nematodes could provide the sensitivity to test for evidence of differential selective pressure affecting specific gene groups.</p> <p>The analysis of TF families in nematodes has revealed several interesting features, such as the high proportion of C2H2 and C4/NHR class of zinc-finger family members relative to the other TF families in all three species (see Figure ##FIG##2##3##). It was previously shown that the NHR family has undergone significant lineage-specific expansion in <italic>C. elegans </italic>and <italic>C. briggsae </italic>[##REF##10022975##53##]. Considering, for example, that <italic>Drosophila </italic>and humans carry less than 50 identified NHR genes (21 and 48, respectively) [##REF##18050471##54##], the presence of more than 200 genes in <italic>Caenorhabditid </italic>species is striking. Although it remains to be seen whether all of these have important roles to play, studies in <italic>C. elegans </italic>have shown that roughly 10% of NHRs mediate diverse processes including molting (<italic>nhr-23, nhr-25, nhr-67</italic>), neuronal differentiation (<italic>unc-55, fax-1</italic>), sex determination (<italic>sex-1</italic>), and dauer formation (<italic>daf-12</italic>) [##REF##18050471##54##]. We found that roughly half of all NHRs in each of the <italic>Caenorhabditid </italic>species are conserved as 3-way best reciprocal orthologs and another 10% exhibit 2-way orthologous relationships with at least one of the other nematode species. The remaining NHRs are likely to have arisen from lineage-specific gene duplications, suggesting that this class of TF may have a significant role in many of those differences that make individual nematode species unique. While the expansion of the NHR family in nematodes is certainly unusual, other TF families show interesting lineage-specific features as well. Previous studies as well as results presented here indicate that TF families such as ZF-C2H2, HOX and T-box have also diverged between the <italic>C. elegans </italic>and <italic>C. briggsae </italic>lineages (see Figure ##FIG##2##3B## and additional file ##SUPPL##8##9##) [##REF##17244357##31##].</p> <p>Our work demonstrates that TF genes are non-randomly distributed in the genomes of both <italic>C. elegans </italic>and <italic>C. briggsae</italic>. We found that members of gene families such as NHR, HOX, and T-box are frequently clustered and present in tandem arrays. A subset of the rapidly evolving NHR family of TF genes in <italic>C. elegans </italic>was previously shown to be located on chromosome V [##REF##10022975##53##,##REF##11275326##55##,##REF##16291650##56##]. We have shown not only that <italic>C. briggsae </italic>exhibits a similar pattern, but also that the majority of the chromosome V NHRs in both species is tandemly arrayed. Our finding that many NHRs appear to be lineage-specific paralogs suggests that gene duplication has played a significant role in the expansion of this gene family in nematodes. The phenomenon of gene clustering has been observed not only in <italic>C. elegans</italic>, but also in other species such as <italic>D. melanogaster </italic>and mouse [##REF##12214599##32##, ####REF##15268860##33##, ##REF##12478293##34####12478293##34##,##REF##11275326##55##,##REF##11279525##57##], and in some cases these clusters are composed of genes that are co-expressed [##REF##12214599##32##,##REF##12478293##34##]. While the precise mechanism of the origin of such clusters remains to be determined, these may be caused by small-scale regional translocations and illegitimate recombination events leading to tandem gene duplications [##REF##10198121##58##,##REF##14704166##59##].</p> <p>Our study has revealed that <italic>C. elegans </italic>TF genes conserved across multiple phyla are more likely to be associated with mutant phenotypes when compared to the remaining TF and non-TF genes. Likewise, the fly, mouse, and human orthologs of <italic>C. elegans </italic>TF genes are enriched in essential genes when compared to <italic>C. elegans </italic>TF genes without detectable orthologs (46%, 50% and 23.3%, respectively). The analysis of the relationship between gene function and interactions revealed that TF genes conserved across phyla exhibit greater number of interactions and mutant phenotypes when compared to those that are divergent. Among the TFs with described interactions, <italic>lin-35 </italic>(human <italic>Rb </italic>ortholog) appears to have an exceptionally large number of interactions. <italic>lin-35 </italic>is known to interact with cell cycle-related and chromatin remodeling factors to regulate tissue growth and morphology [##REF##15574590##60##,##REF##17075059##61##]. We found that among the <italic>lin-35 </italic>interacting genes, 43 (8%) encode TFs, of which 18 have best reciprocal hit orthologs in mouse and human. It is important to keep in mind that conservation in sequence does not indicate the roles of orthologous genes in regulating similar biological processes. Instead, it simply means that genes that are evolutionarily conserved are very likely to play important roles in the development and functioning of the organism. Our results are also consistent with studies in other organisms that have found a significant correlation between connectivity, rate of evolution and gene dispensability (according to lethal or sterile phenotype), even across multiple metazoan phyla. In general, hubs with high degree of connectivity tend to be enriched in essential genes and appear to evolve relatively slower than genes with lower connectivity [##REF##18281268##27##,##REF##16751849##50##,##REF##12769820##62##, ####REF##11976460##63##, ##REF##15746013##64####15746013##64##].</p>
<title>Conclusion</title> <p>This study describes a genome-wide analysis of TF genes in three <italic>Caenorhabditid </italic>nematode species (<italic>C. elegans</italic>, <italic>C. briggsae </italic>and <italic>C. remanei</italic>) as well as their orthologs in fruit fly (<italic>D. melanogaster</italic>), mouse (<italic>M. musculus</italic>) and human (<italic>H. sapiens</italic>). We observed a significantly higher conservation of orthology for the TF genes among <italic>Caenorhabditid </italic>species, while also noting that the coding sequence of TF genes diverges more rapidly than the coding genome average. Finally, the analyses of sequence conservation, gene interactions, and function revealed that TF set conserved in nematodes, fly, mouse, and human is significantly more enriched in essential genes compared to those that lack orthologs in other phyla. Our findings will serve as a resource in aiding us to understand transcriptional networks and their conservation and divergence among metazoa. The compilation of the TF sets also serves as a stepping-stone in generating various resources such as knock-out mutants, cDNA and promoter clones, and reporter gene expressing lines, with the intent of systematically mapping and studying TF function in nematodes. In parallel with many of ongoing initiatives in <italic>C. elegans </italic>these resources will provide foundation for future studies of the conservation of TF function and interaction across the breadth of biodiversity.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>Much of the morphological diversity in eukaryotes results from differential regulation of gene expression in which transcription factors (TFs) play a central role. The nematode <italic>Caenorhabditis elegans </italic>is an established model organism for the study of the roles of TFs in controlling the spatiotemporal pattern of gene expression. Using the fully sequenced genomes of three <italic>Caenorhabditid </italic>nematode species as well as genome information from additional more distantly related organisms (fruit fly, mouse, and human) we sought to identify orthologous TFs and characterized their patterns of evolution.</p> <title>Results</title> <p>We identified 988 TF genes in <italic>C. elegans</italic>, and inferred corresponding sets in <italic>C. briggsae </italic>and <italic>C. remanei</italic>, containing 995 and 1093 TF genes, respectively. Analysis of the three gene sets revealed 652 3-way reciprocal 'best hit' orthologs (nematode TF set), approximately half of which are zinc finger (ZF-C2H2 and ZF-C4/NHR types) and HOX family members. Examination of the TF genes in <italic>C. elegans </italic>and <italic>C. briggsae </italic>identified the presence of significant tandem clustering on chromosome V, the majority of which belong to ZF-C4/NHR family. We also found evidence for lineage-specific duplications and rapid evolution of many of the TF genes in the two species. A search of the TFs conserved among nematodes in <italic>Drosophila melanogaster</italic>, <italic>Mus musculus </italic>and <italic>Homo sapiens </italic>revealed 150 reciprocal orthologs, many of which are associated with important biological processes and human diseases. Finally, a comparison of the sequence, gene interactions and function indicates that nematode TFs conserved across phyla exhibit significantly more interactions and are enriched in genes with annotated mutant phenotypes compared to those that lack orthologs in other species.</p> <title>Conclusion</title> <p>Our study represents the first comprehensive genome-wide analysis of TFs across three nematode species and other organisms. The findings indicate substantial conservation of transcription factors even across distant evolutionary lineages and form the basis for future experiments to examine TF gene function in nematodes and other divergent phyla.</p>
<title>Abbreviations</title> <p><italic>d</italic><sub>N:</sub> non-synonymous substitutions per non-synonymous site; <italic>d</italic><sub>S:</sub> synonymous substitutions per synonymous site; NHR: Nuclear hormone receptor; TF: Transcription factor; ZF: Zinc finger.</p> <title>Authors' contributions</title> <p>The laboratories of BPG and RSS contributed to this publication. BPG and NK identified the <italic>C. elegans </italic>TF set, protein families, chromosomal maps, and mutant phenotypes. WH identified nematode TF orthologs in fly, mouse, and human and carried out most of the sequence alignments, and interaction network analysis. CA performed the InParanoid orthology searches creating the <italic>C. briggsae </italic>and <italic>C. remanei </italic>TF gene sets and constructed NHR phylogenetic tree. BPG, WH and CA drafted the manuscript. BPG conceived and coordinated the study. All authors read and approved the final manuscript.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>We thank Phil Cumbo, Eric Schwartz, Jack Chen, and anonymous reviewers for constructive comments and advice. This work was supported by Natural Sciences and Engineering Research Council of Canada (NSERC) funds to BPG and RSS. CA is supported by an NSERC Post-Graduate Doctoral Scholarship and NK was an NSERC undergraduate summer trainee.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>TF-encoding genes in <italic>C. elegans</italic>, <italic>C. briggsae </italic>and <italic>C. remanei</italic>.</bold> The total number of TF genes in each of the species is given inside the brackets. The numbers of divergent TF genes and those conserved among the three nematode species are shown along the vertices and inside of the triangle, respectively.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p><bold>Sequence divergence of transcription factors in <italic>C. elegans</italic>.</bold> Rates of synonymous substitutions per synonymous site (<italic>d</italic><sub>S</sub>) (A) and non-synonymous substitutions per non-synonymous site (<italic>d</italic><sub>N</sub>) (B) as calculated under model 0 in PAML (Yang) for non-TF (10,827), TF genes conserved in nematodes, fly, mouse, and human (150) and TF genes conserved among the three nematode species (652) are shown. The boxplot indicates the first and third quartiles and the dotted lines the 5th and 95th percentiles. The notches indicate the level of uncertainty associated with the median.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p><bold>Distribution of TF gene families in <italic>C. elegans</italic>.</bold> The pie charts show distributions of all (A), divergent (B), nematode-conserved (C), and nematode-fly-mouse-human-conserved (D) TF genes in <italic>C. elegans</italic>. For details on various gene families please refer to Materials and Methods.</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p><bold>Chromosomal distribution of TF genes in <italic>C. elegans </italic>(A) and <italic>C. briggsae </italic>(B).</bold> The maps have been plotted by taking all TF genes in non-overlapping 200 kb windows. The color codes are as follows. Red: NHR genes, blue: non-NHR TF genes, green: TF genes conserved among the three nematode species. Gene clustering was analyzed by comparing the numbers of TF and non-TF genes located in each window using a χ<sup>2 </sup>test. Gene clusters that are significantly enriched have been marked with stars (*: P &lt; 0.05; ** P &lt; 0.01; ***: P &lt; 0.0001).</p></caption></fig> <fig position="float" id="F5"><label>Figure 5</label><caption><p><bold>Tandem arrays of NHR genes in <italic>C. elegans </italic>and <italic>C. briggsae</italic>.</bold> The snapshots of the genomic regions, visualized by Wormbase genome browser, show 8 genes in <italic>C. elegans </italic>and 6 in <italic>C. briggsae</italic>. The colors of the open reading frames indicate their orientation (blue: leftward, pink: rightward).</p></caption></fig> <fig position="float" id="F6"><label>Figure 6</label><caption><p><bold>Functional classification of <italic>C. elegans </italic>TF genes.</bold> The six broad categories are based on the mutant phenotypes in RNAi studies. Non-TF genes have been plotted for comparison. Viability-E and viability-PE are based on the embryonic and post-embryonic stage lethality phenotypes in RNAi assays, respectively. Refer to text for the description of other categories.</p></caption></fig> <fig position="float" id="F7"><label>Figure 7</label><caption><p><bold>The interaction network of <italic>C. elegans </italic>TF genes.</bold> (A) The breakdowns of TF nodes, non-TF nodes and gene interactions in the network for each of the TF categories. The <italic>C. elegans</italic>-divergent category refers to TF genes that lack unique reciprocal orthologs in other nematode species. (B) The network exhibits several high degree nodes, five of which – <italic>tag-331, eya-1, lin-35, sma-4</italic>, and <italic>pal-1 </italic>– are boxed and shown at high magnification on the right (marked by arrows). The node colors mark different TF genes (red: conserved among nematodes, fly, mouse, and human; yellow: conserved among the three nematode species; green: <italic>C. elegans</italic>-divergent). The network was visualized by using Cytoscape [##REF##14597658##81##].</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>GO term-based searches of TF genes in <italic>C. elegans</italic>.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center">GO ID</td><td align="left">Term</td><td align="center">TF genes</td><td align="center">Unique</td></tr></thead><tbody><tr><td align="center">0003700</td><td align="left">Transcription Factor activity</td><td align="center">614</td><td align="center">6</td></tr><tr><td align="center">0043565</td><td align="left">DNA binding, sequence specific</td><td align="center">515</td><td align="center">6</td></tr><tr><td align="center">0003677</td><td align="left">DNA binding</td><td align="center">858</td><td align="center">352</td></tr><tr><td align="center">0030528</td><td align="left">Transcription regulator activity</td><td align="center">75</td><td align="center">9</td></tr><tr><td align="center">0006355</td><td align="left">Regulation of transcription, DNA dependent</td><td align="center">768</td><td align="center">78</td></tr><tr><td align="center">0045449</td><td align="left">Regulation of transcription</td><td align="center">199</td><td align="center">19</td></tr><tr><td align="center">0000122</td><td align="left">Negative regulation of transcription from RNA pol II promoter</td><td align="center">8</td><td align="center">4</td></tr><tr><td align="center">004544</td><td align="left">Positive regulation of transcription from RNA pol II promoter</td><td align="center">24</td><td align="center">10</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>The breakdowns of TF genes in each of the nematode species genomes based on various search categories.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>Search method</bold></td><td align="center" colspan="3"><bold>Number of TF genes</bold></td></tr><tr><td/><td colspan="3"><hr/></td></tr><tr><td/><td align="center"><bold><italic>C. elegans</italic></bold></td><td align="center"><bold><italic>C. briggsae</italic></bold></td><td align="center"><bold><italic>C. remanei</italic></bold></td></tr></thead><tbody><tr><td align="left">GO term-based</td><td align="center">707</td><td align="center">ND</td><td align="center">NA</td></tr><tr><td align="left">Orthologs (InParanoid and reciprocal BLAST)</td><td align="center">ND</td><td align="center">713</td><td align="center">703</td></tr><tr><td align="left">Manual curation</td><td align="center">281</td><td align="center">NA</td><td align="center">NA</td></tr><tr><td align="left">HMM alignments</td><td align="center">ND</td><td align="center">282</td><td align="center">390</td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left"><bold>TOTAL:</bold></td><td align="center"><bold>988</bold></td><td align="center"><bold>995</bold></td><td align="center"><bold>1093</bold></td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p>Chromosome-wise breakdown of TF gene clusters in <italic>C. elegans </italic>and <italic>C. briggsae</italic>.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center"><bold>Chromosome</bold></td><td align="center" colspan="2"><bold><italic>C. elegans</italic></bold></td><td align="center" colspan="2"><bold><italic>C. briggsae</italic></bold></td></tr><tr><td/><td colspan="2"><hr/></td><td colspan="2"><hr/></td></tr><tr><td/><td align="center"><bold>Number of clusters</bold></td><td align="center"><bold>Number of genes</bold></td><td align="center"><bold>Number of clusters</bold></td><td align="center"><bold>Number of genes</bold></td></tr></thead><tbody><tr><td align="center">1</td><td align="center">2</td><td align="center">12</td><td align="center">3</td><td align="center">19</td></tr><tr><td align="center">2</td><td align="center">4</td><td align="center">24</td><td align="center">5</td><td align="center">34</td></tr><tr><td align="center">3</td><td align="center">1</td><td align="center">5</td><td align="center">3</td><td align="center">17</td></tr><tr><td align="center">4</td><td align="center">5</td><td align="center">29</td><td align="center">4</td><td align="center">21</td></tr><tr><td align="center">5</td><td align="center">10</td><td align="center">97</td><td align="center">7</td><td align="center">64</td></tr><tr><td align="center">X</td><td align="center">3</td><td align="center">16</td><td align="center">5</td><td align="center">29</td></tr><tr><td colspan="5"><hr/></td></tr><tr><td align="center"><bold>TOTAL</bold></td><td align="center"><bold>25</bold></td><td align="center"><bold>183</bold></td><td align="center"><bold>27</bold></td><td align="center"><bold>184</bold></td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T4"><label>Table 4</label><caption><p>Genetic disorders linked to human TF genes conserved among nematodes, fly, mouse, and human.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold><italic>C. elegans </italic>gene</bold></td><td align="left"><bold>Human ortholog</bold></td><td align="left"><bold>Human disorder</bold></td></tr></thead><tbody><tr><td align="left"><italic>vab-3</italic></td><td align="left"><italic>Pax6</italic></td><td align="left">Aniridia type II, Peters anomaly with cataract, foveal hypoplasia</td></tr><tr><td align="left"><italic>Y38H8A.5</italic></td><td align="left"><italic>FEZF1</italic></td><td align="left">Beckwith-Wiedemann syndrome</td></tr><tr><td align="left"><italic>ceh-33</italic></td><td align="left"><italic>SIX1</italic></td><td align="left">Branchiootic syndrome 3</td></tr><tr><td align="left"><italic>mab-9</italic></td><td align="left"><italic>TBX20</italic></td><td align="left">Cardiomyopathy, atrial septal defect 1</td></tr><tr><td align="left"><italic>tag-192</italic></td><td align="left"><italic>CHD7</italic></td><td align="left">CHARGE syndrome</td></tr><tr><td align="left"><italic>ceh-24</italic></td><td align="left"><italic>TITF1</italic></td><td align="left">Congenital hypothyroidism, neonatal respiratory insufficiency</td></tr><tr><td align="left"><italic>dve-1</italic></td><td align="left"><italic>SATB2</italic></td><td align="left">Cleft palate isolated</td></tr><tr><td align="left"><italic>ceh-14</italic></td><td align="left"><italic>LHX3</italic></td><td align="left">Combined pituitary hormone deficiency 3</td></tr><tr><td align="left"><italic>unc-86</italic></td><td align="left"><italic>Pou4f3</italic></td><td align="left">DFNA15 syndrome</td></tr><tr><td align="left"><italic>elt-1</italic></td><td align="left"><italic>GATA1</italic></td><td align="left">Dyserythropoietic anemia with thrombocytopenia</td></tr><tr><td align="left"><italic>K02H8.1</italic></td><td align="left"><italic>MBNL2</italic></td><td align="left">Dystrophia myotonica 1</td></tr><tr><td align="left"><italic>fax-1</italic></td><td align="left"><italic>Nr2e3</italic></td><td align="left">Enhanced s-cone syndrome</td></tr><tr><td align="left"><italic>ceh-17</italic></td><td align="left"><italic>PHOX2A</italic></td><td align="left">Congenital fibrosis of the extraocular muscles 2</td></tr><tr><td align="left"><italic>ceh-32</italic></td><td align="left"><italic>SIX3</italic></td><td align="left">Holoprosencephaly 2</td></tr><tr><td align="left"><italic>sbp-1</italic></td><td align="left"><italic>Srebf1</italic></td><td align="left">Hypercholesterolemia, familial</td></tr><tr><td align="left"><italic>lin-28</italic></td><td align="left"><italic>LIN28B</italic></td><td align="left">Hypomyelination and cataract</td></tr><tr><td align="left"><italic>alr-1</italic></td><td align="left"><italic>ARX</italic></td><td align="left">Lissencephaly, X-linked, with ambiguous genitalia</td></tr><tr><td align="left"><italic>hmg-5</italic></td><td align="left"><italic>Tfam</italic></td><td align="left">Kearns-Sayre syndrome</td></tr><tr><td align="left"><italic>cnd-1</italic></td><td align="left"><italic>NEUROD1</italic></td><td align="left">Maturity-onset diabetes of the young</td></tr><tr><td align="left"><italic>lim-6</italic></td><td align="left"><italic>LMX1B</italic></td><td align="left">Nail patella syndrome NPS1</td></tr><tr><td align="left"><italic>grh-1</italic></td><td align="left"><italic>GRHL2</italic></td><td align="left">Neurosensory deafness 28</td></tr><tr><td align="left"><italic>sma-4</italic></td><td align="left"><italic>Smad4</italic></td><td align="left">Pancreatic cancer, Hemorrhagic Telangiectasia Syndrome (HTT)</td></tr><tr><td align="left"><italic>nhr-6</italic></td><td align="left"><italic>NR4A2</italic></td><td align="left">PARK14</td></tr><tr><td align="left"><italic>ceh-6</italic></td><td align="left"><italic>POU3F3</italic></td><td align="left">Perilymphatic gusher-deafness syndrome</td></tr><tr><td align="left"><italic>zag-1</italic></td><td align="left"><italic>ZEB1</italic></td><td align="left">Posterior polymorphous corneal dystrophy 3</td></tr><tr><td align="left"><italic>eor-1</italic></td><td align="left"><italic>MYNN</italic></td><td align="left">Promyelocytic leukemia</td></tr><tr><td align="left"><italic>R07E5.3</italic></td><td align="left"><italic>Smarcb1</italic></td><td align="left">Rhabdoid tumor</td></tr><tr><td align="left"><italic>cbp-1</italic></td><td align="left"><italic>CREBBP</italic></td><td align="left">Rubinstein-taybi syndrome, acute myeloid leukemia</td></tr><tr><td align="left"><italic>ceh-43</italic></td><td align="left"><italic>DLX5</italic></td><td align="left">Split-hand/foot malformation</td></tr><tr><td align="left"><italic>ing-3</italic></td><td align="left"><italic>ING3</italic></td><td align="left">Squamous cell carcinoma</td></tr><tr><td align="left"><italic>ast-1</italic></td><td align="left"><italic>FLI1</italic></td><td align="left">Thrombocytopenia, Paris-Trousseau type</td></tr><tr><td align="left"><italic>nhr-64</italic></td><td align="left"><italic>HNF4A</italic></td><td align="left">Maturity-onset diabetes of the young</td></tr><tr><td align="left"><italic>tbx-2</italic></td><td align="left"><italic>Tbx2</italic></td><td align="left">Ulnar mammary syndrome</td></tr><tr><td align="left"><italic>K02D7.2</italic></td><td align="left"><italic>SNAI2</italic></td><td align="left">Waardenburg syndrome, piebaldism</td></tr><tr><td align="left"><italic>F53F8.1</italic></td><td align="left"><italic>KLF3</italic></td><td align="left">Wilms tumor</td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p>List of 314 incorrect entries (non-TFs) in <italic>C. elegans</italic>.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional file 2</title><p>List of 167 genes that encode chromatin remodeling, general transcription and DNA/RNA binding factors in <italic>C. elegans</italic>.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S3"><caption><title>Additional file 3</title><p>List of 83 histone-encoding genes in <italic>C. elegans</italic>.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S4"><caption><title>Additional file 4</title><p>List of 988 TF genes in <italic>C. elegans</italic>.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S5"><caption><title>Additional file 5</title><p>List of 17 false positives in wTF2.0.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S6"><caption><title>Additional file 6</title><p>List of 995 TF genes in <italic>C. briggsae</italic>.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S7"><caption><title>Additional file 7</title><p>List of 1093 TF genes in <italic>C. remanei</italic>.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S8"><caption><title>Additional file 8</title><p>List of 150 TF orthologs in <italic>Drosophila melanogaster</italic>, <italic>M. musculus</italic>, and <italic>H. sapiens</italic>.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S9"><caption><title>Additional file 9</title><p>BLASTP hits of <italic>C. briggsae</italic>-divergent TFs in <italic>C. elegans </italic>genome.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S10"><caption><title>Additional file 10</title><p>Tandem arrays of TF genes in <italic>C. elegans</italic>.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S11"><caption><title>Additional file 11</title><p>Tandem arrays of TF genes in <italic>C. briggsae</italic>.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S12"><caption><title>Additional file 12</title><p>Phylogenetic tree of NHR genes in <italic>Caenorhabditid </italic>nematode species. Colors mark NHR genes in different species (blue: <italic>C. elegans</italic>, red: <italic>C. briggsae </italic>and light green: <italic>C. remanei</italic>). Tandemly along chromosomes and phylogenetically clustered genes are indicated by vertical bars. Chromosomes carrying NHR clusters are indicated by roman numerals. The sub-branch comprised of seven groups of NHR genes on chromosome V has been marked by a star (*). Scale bar represents 0.5 substitutions per site.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S13"><caption><title>Additional file 13</title><p><italic>C. elegans </italic>genes and their mutant phenotypes in RNAi assays.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S14"><caption><title>Additional file 14</title><p>RNAi phenotypes sorted into six broad categories.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S15"><caption><title>Additional file 15</title><p>List of the worm (<italic>C. elegans</italic>) fly (<italic>D. melanogaster</italic>) and mouse (<italic>M. musculus</italic>) mutant phenotypes associated with conserved TF genes.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S16"><caption><title>Additional file 16</title><p>List of TF genes and their interactors. The columns A and B merely list gene pairs that interact with each other and do not indicate the direction of regulation.</p></caption></supplementary-material>
<table-wrap-foot><p>The table lists numbers of all TF genes as well as those uniquely identified by each of the terms.</p></table-wrap-foot> <table-wrap-foot><p>NA: not applicable, ND: not done.</p></table-wrap-foot>
<graphic xlink:href="1471-2164-9-399-1"/> <graphic xlink:href="1471-2164-9-399-2"/> <graphic xlink:href="1471-2164-9-399-3"/> <graphic xlink:href="1471-2164-9-399-4"/> <graphic xlink:href="1471-2164-9-399-5"/> <graphic xlink:href="1471-2164-9-399-6"/> <graphic xlink:href="1471-2164-9-399-7"/>
<media xlink:href="1471-2164-9-399-S1.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-399-S2.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-399-S3.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-399-S4.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-399-S5.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-399-S6.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-399-S7.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-399-S8.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-399-S9.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-399-S10.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-399-S11.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-399-S12.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-399-S13.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-399-S14.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-399-S15.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-399-S16.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media>
[{"article-title": ["Wormbase"]}, {"article-title": ["Wormbase FTP site"]}, {"italic": ["C. remanei"]}, {"article-title": ["The HMMER software package"]}, {"article-title": ["The BioGRID"]}, {"surname": ["Simmer", "Moorman", "van der Linden", "Kuijk", "van den Berghe", "Kamath", "Fraser", "Ahringer", "Plasterk"], "given-names": ["F", "C", "AM", "E", "PV", "RS", "AG", "J", "RH"], "article-title": ["Genome-wide RNAi of "], "italic": ["C. elegans"], "source": ["PLoS Biol"], "year": ["2003"], "volume": ["1"], "fpage": ["77"], "lpage": ["84"], "pub-id": ["10.1371/journal.pbio.0000012"]}, {"article-title": ["Pfam"]}, {"article-title": ["The InParanoid database"], "fpage": ["Eukaryotic ortholog groups"]}, {"article-title": ["Interpro"]}, {"surname": ["Felsenstein"], "given-names": ["J"], "article-title": ["The PHYLIP package"]}, {"article-title": ["Flybase"]}, {"article-title": ["OMIM - Online Mendelian Inheritance in Man"]}, {"article-title": ["The NCBI PubMed"]}]
BMC Genomics. 2008 Aug 27; 9:399