A compendium of Caenorhabditis elegans regulatory transcription factors: a resource for mapping transcription regulatory networks

Functional modularity of nuclear hormone receptors in a Caenorhabditis elegans metabolic gene regulatory network

We present the first gene regulatory network (GRN) that pertains to post-developmental gene expression. Specifically, we mapped a transcription regulatory network of Caenorhabditis elegans metabolic gene promoters using gene-centered yeast one-hybrid assays. We found that the metabolic GRN is enriched for nuclear hormone receptors (NHRs) compared with other gene-centered regulatory networks, and that these NHRs organize into functional network modules.

The NHR family has greatly expanded in nematodes; C. elegans has 284 NHRs, whereas humans have only 48. We show that the NHRs in the metabolic GRN have metabolic phenotypes, suggesting that they do not simply function redundantly.

The mediator subunit MDT-15 preferentially interacts with NHRs that occur in the metabolic GRN.

We describe an NHR circuit that responds to nutrient availability and propose a model for the evolution and organization of NHRs in C. elegans metabolic regulatory networks.

Physical and/or regulatory interactions between transcription factors (TFs) and their target genes are essential to establish body plans of multicellular organisms during development, and these interactions have been studied extensively in the context of GRNs. The precise control of differential gene expression is also of critical importance to maintain physiological homeostasis, and many metabolic disorders such as obesity and diabetes coincide with substantial changes in gene expression. Much work has focused on the GRNs that control metazoan development; however, the design principles and organization of the GRNs that control systems physiology remain largely unexplored.

In this study, we present the first gene-centered GRN that includes ∼70 genes involved in C. elegans metabolism and physiology, 100 TFs and more than 500 protein–DNA interactions between them. The resulting metabolic GRN is enriched for NHRs, compared with other gene-centered regulatory networks. NHRs are well-known regulators of lipid meta-qj;bolism in mammals. The transcriptional activity of NHRs can be modified by diffusible ligands, which allows these TFs to function as molecular sensors and rapidly alter the expression of their target genes. Interestingly, NHRs comprise the largest family of TFs in nematodes; the C. elegans genome encodes 284 NHRs, most of which are uncharacterized. Furthermore, their organization in GRNs has not yet been investigated. In our study, we show that the C. elegans NHRs that we retrieved in the metabolic GRN organize into network modules, and that most of these NHRs function to maintain lipid homeostasis in the nematode. Interestingly, network modularity has been proposed to facilitate rapid and robust changes in gene expression. Our results suggest that the C. elegans metabolic GRN may have evolved by combining NHR family expansion with the specific modular wiring of NHRs to enable the rapid adaptation of the animal to different environmental cues.

NHRs can interact with transcriptional cofactors such as chromatin remodeling complexes and Mediator components. For instance, the C. elegans Mediator subunit, MDT-15, can interact with NHR-49 to regulate the expression of its target genes. To find all the TFs that MDT-15 can interact with, we performed systematic yeast two-hybrid assays with MDT-15 versus 755 full-length TFs. We found that MDT-15 preferentially associates with NHRs, and specifically with those NHRs that confer a metabolic phenotype and that occur in the metabolic GRN. This illustrates the central role of MDT-15 in the regulation of metabolic gene expression.

Using a variety of genetic and biochemical approaches, we characterized NHR-86 in more detail. NHR-86 participates in one of the two NHR modules, and has a high-flux capacity; that is it has both a high incoming and a high outgoing degree. We obtained an nhr-86 mutant and generated an NHR-86 antibody, and showed that NHR-86 functions as an auto-repressor in vivo and that nhr-86 mutant animals store abnormally high levels of body fat.

Finally, we discovered a novel NHR circuit that responds to nutrient availability. In this circuit NHR-45 regulates the activity of nhr-178 promoter in two distinct physiologically important tissues: the intestine and the hypodermis. Both of these NHRs are required to maintain lipid homeostasis in C. elegans. The expression of nhr-178 is responsive to the nutritional status of the animal, which switches between ON and OFF states in the hypodermis. We found that NHR-45 activity is necessary to control this switch in the hypodermis. Interestingly, NHR-45 has opposite effects on the activity of the nhr-178 promoter in these tissues: NHR-45 activates this promoter in the intestine, but represses it in the hypodermis.

Altogether our study leads to a model in which the expansion of the NHR family, TFs that have the capacity to act as fast molecular sensors, is combined with a modular network organization to enable rapid and robust responses to various environmental cues.

Gene regulatory networks (GRNs) provide insights into the mechanisms of differential gene expression at a systems level. GRNs that relate to metazoan development have been studied extensively. However, little is still known about the design principles, organization and functionality of GRNs that control physiological processes such as metabolism, homeostasis and responses to environmental cues. In this study, we report the first experimentally mapped metazoan GRN of Caenorhabditis elegans metabolic genes. This network is enriched for nuclear hormone receptors (NHRs). The NHR family has greatly expanded in nematodes: humans have 48 NHRs, but C. elegans has 284, most of which are uncharacterized. We find that the C. elegans metabolic GRN is highly modular and that two GRN modules predominantly consist of NHRs. Network modularity has been proposed to facilitate a rapid response to different cues. As NHRs are metabolic sensors that are poised to respond to ligands, this suggests that C. elegans GRNs evolved to enable rapid and adaptive responses to different cues by a concurrence of NHR family expansion and modular GRN wiring.

Integrating one-dimensional and three-dimensional maps of genomes

Genomes exist in vivo as complex physical structures, and their functional output (i.e. the gene expression profile of a cell) is related to their spatial organization inside the nucleus as well as to local chromatin status. Chromatin modifications and chromosome conformation are distinct in different tissues and cell types, which corresponds closely with the diversity in gene-expression patterns found in different tissues of the body. The biological processes and mechanisms driving these general correlations are currently the topic of intense study. An emerging theme is that genome compartmentalization - both along the linear length of chromosomes, and in three dimensions by the spatial colocalization of chromatin domains and genomic loci from across the genome - is a crucial parameter in regulating genome expression. In this Commentary, we propose that a full understanding of genome regulation requires integrating three different types of data: first, one-dimensional data regarding the state of local chromatin - such as patterns of protein binding along chromosomes; second, three-dimensional data that describe the population-averaged folding of chromatin inside cells and; third, single-cell observations of three-dimensional spatial colocalization of genetic loci and trans factors that reveal information about their dynamics and frequency of colocalization.

Cell fate specification in the C. elegans embryo

Cell specification requires that particular subsets of cells adopt unique expression patterns that ultimately define the fates of their descendants. In C. elegans, cell fate specification involves the combinatorial action of multiple signals that produce activation of a small number of "blastomere specification" factors. These initiate expression of gene regulatory networks that drive development forward, leading to activation of "tissue specification" factors. In this review, the C. elegans embryo is considered as a model system for studies of cell specification. The techniques used to study cell fate in this species, and the themes that have emerged, are described.

Recent advances in understanding the molecular mechanisms regulating C. elegans transcription

We review recent studies that have advanced our understanding of the molecular mechanisms regulating transcription in the nematode C. elegans. Topics covered include: (i) general properties of C. elegans promoters; (ii) transcription factors and transcription factor combinations involved in cell fate specification and cell differentiation; (iii) new roles for general transcription factors; (iv) nucleosome positioning in C. elegans "chromatin"; and (v) some characteristics of histone variants and histone modifications and their possible roles in controlling C. elegans transcription.

The matrix peptide exporter HAF-1 signals a mitochondrial UPR by activating the transcription factor ZC376.7 in C. elegans

Genetic analyses previously implicated the matrix-localized protease ClpP in signaling the stress of protein misfolding in the mitochondrial matrix to activate nuclear-encoded mitochondrial chaperone genes in C. elegans (UPR(mt)). Here, we report that haf-1, a gene encoding a mitochondria-localized ATP-binding cassette protein, is required for signaling within the UPR(mt) and for coping with misfolded protein stress. Peptide efflux from isolated mitochondria was ATP dependent and required HAF-1 and the protease ClpP. Defective UPR(mt) signaling in the haf-1-deleted worms was associated with failure of the bZIP protein, ZC376.7, to localize to nuclei in worms with perturbed mitochondrial protein folding, whereas zc376.7(RNAi) strongly inhibited the UPR(mt). These observations suggest a simple model whereby perturbation of the protein-folding environment in the mitochondrial matrix promotes ClpP-mediated generation of peptides whose haf-1-dependent export from the matrix contributes to UPR(mt) signaling across the mitochondrial inner membrane.

bZIP transcription factor zip-2 mediates an early response to Pseudomonas aeruginosa infection in Caenorhabditis elegans

Very little is known about how animals discriminate pathogens from innocuous microbes. To address this question, we examined infection-response gene induction in the nematode Caenorhabditis elegans. We focused on genes that are induced in C. elegans by infection with the bacterial pathogen Pseudomonas aeruginosa, but are not induced by an isogenic attenuated gacA mutant. Most of these genes are induced independently of known immunity pathways. We generated a GFP reporter for one of these genes, infection response gene 1 (irg-1), which is induced strongly by wild-type P. aeruginosa strain PA14, but not by other C. elegans pathogens or by other wild-type P. aeruginosa strains that are weakly pathogenic to C. elegans. To identify components of the pathway that induces irg-1 in response to infection, we performed an RNA interference screen of C. elegans transcription factors. This screen identified zip-2, a bZIP transcription factor that is required for inducing irg-1, as well as several other genes, and is important for defense against infection by P. aeruginosa. These data indicate that zip-2 is part of a specialized pathogen response pathway that is induced by virulent strains of P. aeruginosa and provides defense against this pathogen.

Gene-centered regulatory networks

Differential gene expression plays a critical role in the development and physiology of multicellular organisms. At a 'systems level' (e.g. at the level of a tissue, organ or whole organism), this process can be studied using gene regulatory network (GRN) models that capture physical and regulatory interactions between genes and their regulators. In the past years, significant progress has been made toward the mapping of GRNs using a variety of experimental and computational approaches. Here, we will discuss gene-centered approaches that we employed to characterize GRNs and describe insights that we have obtained into the global design principles of gene regulation in complex metazoan systems.

DNA-centered approaches to characterize regulatory protein-DNA interaction complexes

Gene regulation is mediated by site-specific DNA-binding proteins or transcription factors (TFs), which form protein complexes at regulatory loci either to activate or repress the expression of a target gene. The study of the dynamic properties of these regulatory DNA-binding complexes has so far been dominated by protein-centered methodologies, aiming to characterize the DNA-binding behavior of one specific protein at a time. With the emerging evidence for a role of DNA in allosterically influencing DNA-binding protein complex formation, there is renewed interest in DNA-centered approaches to capture protein complexes on defined regulatory loci and to correlate changes in their composition with alterations in target gene expression. In this review, we present the current state-of-the-art in such DNA-centered approaches and evaluate recent technological improvements in the purification as well as in the identification of regulatory DNA-binding protein complexes within or outside their biological context. Finally, we suggest possible areas of improvement and assess the putative impact of DNA-centered methodologies on the gene regulation field for the forthcoming years.

In silico analysis of transcription factor repertoire and prediction of stress responsive transcription factors in soybean

Sequence-specific DNA-binding transcription factors (TFs) are often termed as 'master regulators' which bind to DNA and either activate or repress gene transcription. We have computationally analysed the soybean genome sequence data and constructed a proper set of TFs based on the Hidden Markov Model profiles of DNA-binding domain families. Within the soybean genome, we identified 4342 loci encoding 5035 TF models which grouped into 61 families. We constructed a database named SoybeanTFDB (http://soybeantfdb.psc.riken.jp) containing the full compilation of soybean TFs and significant information such as: functional motifs, full-length cDNAs, domain alignments, promoter regions, genomic organization and putative regulatory functions based on annotations of gene ontology (GO) inferred by comparative analysis with Arabidopsis. With particular interest in abiotic stress signalling, we analysed the promoter regions for all of the TF encoding genes as a means to identify abiotic stress responsive cis-elements as well as all types of cis-motifs provided by the PLACE database. SoybeanTFDB enables scientists to easily access cis-element and GO annotations to aid in the prediction of TF function and selection of TFs with functions of interest. This study provides a basic framework and an important user-friendly public information resource which enables analyses of transcriptional regulation in soybean.

Comparison of diverse developmental transcriptomes reveals that coexpression of gene neighbors is not evolutionarily conserved

Genomic analyses have shown that adjacent genes are often coexpressed. However, it remains unclear whether the observed coexpression is a result of functional organization or a consequence of adjacent active chromatin or transcriptional read-through, which may be free of selective biases. Here, we compare temporal expression profiles of one-to-one orthologs in conserved or divergent genomic positions in two genetically distant nematode species-Caenorhabditis elegans and C. briggsae-that share a near-identical developmental program. We find, for all major patterns of temporal expression, a substantive amount of gene expression divergence. However, this divergence is not random: Genes that function in essential developmental processes show less divergence than genes whose functions are not required for viability. Coexpression of gene neighbors in either species is highly divergent in the other, in particular when the neighborhood is not conserved. Interestingly, essential genes appear to maintain their expression profiles despite changes in neighborhoods suggesting exposure to stronger selection. Our results suggest that a significant fraction of the coexpression observed among gene neighbors may be accounted for by neutral processes, and further that these may be distinguished by comparative gene expression analyses.

A 3D digital atlas of C. elegans and its application to single-cell analyses

We built a digital nuclear atlas of the newly hatched, first larval stage (L1) of the wild-type hermaphrodite of Caenorhabditis elegans at single-cell resolution from confocal image stacks of 15 individual worms. The atlas quantifies the stereotypy of nuclear locations and provides other statistics on the spatial patterns of the 357 nuclei that could be faithfully segmented and annotated out of the 558 present at this developmental stage. We then developed an automated approach to assign cell names to each nucleus in a three-dimensional image of an L1 worm. We achieved 86% accuracy in identifying the 357 nuclei automatically. This computational method will allow high-throughput single-cell analyses of the post-embryonic worm, such as gene expression analysis, or ablation or stimulation of cells under computer control in a high-throughput functional screen.

Unlocking the secrets of the genome

Despite the successes of genomics, little is known about how genetic information produces complex organisms. A look at the crucial functional elements of fly and worm genomes could change that.

The National Human Genome Research Institute's modENCODE project (the model organism ENCyclopedia Of DNA Elements) was set up in 2007 with the goal of identifying all the sequence-based functional elements in the genomes of two important experimental organisms, Caenorhabditis elegans and Drosophila melanogaster. Armed with modENCODE data, geneticists will be able to undertake the comprehensive molecular studies of regulatory networks that hold the key to how complex multicellular organisms arise from the list of instructions coded in the genome. In this issue, modENCODE team members outline their plan of campaign. Data from the project are to be made available on http://www.modencode.org and elsewhere as the work progresses.

A multiparameter network reveals extensive divergence between C. elegans bHLH transcription factors

Differences in expression, protein interactions, and DNA binding of paralogous transcription factors ("TF parameters") are thought to be important determinants of regulatory and biological specificity. However, both the extent of TF divergence and the relative contribution of individual TF parameters remain undetermined. We comprehensively identify dimerization partners, spatiotemporal expression patterns, and DNA-binding specificities for the C. elegans bHLH family of TFs, and model these data into an integrated network. This network displays both specificity and promiscuity, as some bHLH proteins, DNA sequences, and tissues are highly connected, whereas others are not. By comparing all bHLH TFs, we find extensive divergence and that all three parameters contribute equally to bHLH divergence. Our approach provides a framework for examining divergence for other protein families in C. elegans and in other complex multicellular organisms, including humans. Cross-species comparisons of integrated networks may provide further insights into molecular features underlying protein family evolution. For a video summary of this article, see the PaperFlick file available with the online Supplemental Data.

Experimental advances in the characterization of metazoan gene regulatory networks

Gene regulatory networks (GRNs) play a vital role in metazoan development and function, and deregulation of these networks is often implicated in disease. GRNs depict the dynamic interactions between genomic and regulatory state components. The genomic components comprise genes and their associated cis-regulatory elements. The regulatory state components consist primarily of transcriptional complexes that bind the latter elements. With the availability of complete genome sequences, several approaches have recently been developed which promise to significantly enhance our ability to identify either the genomic or regulatory state components, or the interactions between these two. In this review, I provide an in-depth overview of these approaches and detail how each contributes to a more comprehensive understanding of GRN composition and function.

Rapid sequence evolution of transcription factors controlling neuron differentiation in Caenorhabditis

Whether phenotypic evolution proceeds predominantly through changes in regulatory sequences is a controversial issue in evolutionary genetics. Ample evidence indicates that the evolution of gene regulatory networks via changes in cis-regulatory sequences is an important determinant of phenotypic diversity. However, recent experimental work suggests that the role of transcription factor (TF) divergence in developmental evolution may be underestimated. In order to help understand what levels of constraints are acting on the coding sequence of developmental regulatory genes, evolutionary rates were investigated among 48 TFs required for neuronal development in Caenorhabditis elegans. Allelic variation was then sampled for 28 of these genes within a population of the related species Caenorhabditis remanei. Neuronal TFs are more divergent, both within and between species, than structural genes. TFs affecting different neuronal classes are under different levels of selective constraints. The regulatory genes controlling the differentiation of chemosensory neurons evolve particularly fast and exhibit higher levels of within- and between-species nucleotide variation than TFs required for the development of several neuronal classes and TFs required for motorneuron differentiation. The TFs affecting chemosensory neuron development are also more divergent than chemosensory genes expressed in the neurons they differentiate. These results illustrate that TFs are not as highly constrained as commonly thought and suggest that the role of divergence in developmental regulatory genes during the evolution of gene regulatory networks requires further attention.

Universal protein-binding microarrays for the comprehensive characterization of the DNA-binding specificities of transcription factors

Protein-binding microarray (PBM) technology provides a rapid, high-throughput means of characterizing the in vitro DNA-binding specificities of transcription factors (TFs). Using high-density, custom-designed microarrays containing all 10-mer sequence variants, one can obtain comprehensive binding-site measurements for any TF, regardless of its structural class or species of origin. Here, we present a protocol for the examination and analysis of TF-binding specificities at high resolution using such 'all 10-mer' universal PBMs. This procedure involves double-stranding a commercially synthesized DNA oligonucleotide array, binding a TF directly to the double-stranded DNA microarray and labeling the protein-bound microarray with a fluorophore-conjugated antibody. We describe how to computationally extract the relative binding preferences of the examined TF for all possible contiguous and gapped 8-mers over the full range of affinities, from highest affinity sites to nonspecific sites. Multiple proteins can be tested in parallel in separate chambers on a single microarray, enabling the processing of a dozen or more TFs in a single day.

A census of human transcription factors: function, expression and evolution

Transcription factors are key cellular components that control gene expression: their activities determine how cells function and respond to the environment. Currently, there is great interest in research into human transcriptional regulation. However, surprisingly little is known about these regulators themselves. For example, how many transcription factors does the human genome contain? How are they expressed in different tissues? Are they evolutionarily conserved? Here, we present an analysis of 1,391 manually curated sequence-specific DNA-binding transcription factors, their functions, genomic organization and evolutionary conservation. Much remains to be explored, but this study provides a solid foundation for future investigations to elucidate regulatory mechanisms underlying diverse mammalian biological processes.

The interplay between transcription factors and microRNAs in genome-scale regulatory networks

Metazoan genomes contain thousands of protein-coding and non-coding RNA genes, most of which are differentially expressed, i.e., at different locations, at different times during development, or in response to environmental signals. Differential gene expression is achieved through complex regulatory networks that are controlled in part by two types of trans-regulators: transcription factors (TFs) and microRNAs (miRNAs). TFs bind to cis-regulatory DNA elements that are often located in or near their target genes, while miRNAs hybridize to cis-regulatory RNA elements mostly located in the 3' untranslated region of their target mRNAs. Here, we describe how these trans-regulators interact with each other in the context of gene regulatory networks to coordinate gene expression at the genome-scale level, and discuss future challenges of integrating these networks with other types of functional networks.

TFCat: the curated catalog of mouse and human transcription factors

TFCat is a catalog of mouse and human transcription factors based on a reliable core collection of annotations obtained by expert review of the scientific literature

Unravelling regulatory programs governed by transcription factors (TFs) is fundamental to understanding biological systems. TFCat is a catalog of mouse and human TFs based on a reliable core collection of annotations obtained by expert review of the scientific literature. The collection, including proven and homology-based candidate TFs, is annotated within a function-based taxonomy and DNA-binding proteins are organized within a classification system. All data and user-feedback mechanisms are available at the TFCat portal .

A C. elegans genome-scale microRNA network contains composite feedback motifs with high flux capacity

MicroRNAs (miRNAs) and transcription factors (TFs) are primary metazoan gene regulators. Whereas much attention has focused on finding the targets of both miRNAs and TFs, the transcriptional networks that regulate miRNA expression remain largely unexplored. Here, we present the first genome-scale Caenorhabditis elegans miRNA regulatory network that contains experimentally mapped transcriptional TF --> miRNA interactions, as well as computationally predicted post-transcriptional miRNA --> TF interactions. We find that this integrated miRNA network contains 23 miRNA <--> TF composite feedback loops in which a TF that controls a miRNA is itself regulated by that same miRNA. By rigorous network randomizations, we show that such loops occur more frequently than expected by chance and, hence, constitute a genuine network motif. Interestingly, miRNAs and TFs in such loops are heavily regulated and regulate many targets. This "high flux capacity" suggests that loops provide a mechanism of high information flow for the coordinate and adaptable control of miRNA and TF target regulons.

The FLYWCH transcription factors FLH-1, FLH-2, and FLH-3 repress embryonic expression of microRNA genes in C. elegans

MicroRNAs (miRNAs) are small noncoding RNAs that regulate gene expression post-transcriptionally via antisense base-pairing. Although miRNAs are involved in a variety of important biological functions, little is known about their transcriptional regulation. Using yeast one-hybrid assays, we identified transcription factors with a FLYWCH Zn-finger DNA-binding domain that bind to the promoters of several Caenorhabditis elegans miRNA genes. The products of the flh-1 and flh-2 genes function redundantly to repress embryonic expression of lin-4, mir-48, and mir-241, miRNA genes that are normally expressed only post-embryonically. Although single mutations in either flh-1 or flh-2 genes result in a viable phenotype, double mutation of flh-1 and flh-2 results in early larval lethality and an enhanced derepression of their target miRNAs in embryos. Double mutations in flh-2 and a third FLYWCH Zn-finger-containing transcription factor, flh-3, also result in enhanced precocious expression of target miRNAs. Mutations of lin-4 or mir-48&mir-241 do not rescue the lethal flh-1; flh-2 double-mutant phenotype, suggesting that the inviability is not solely the result of precocious expression of these miRNAs. Therefore, the FLH-1 and FLH-2 proteins likely play a more general role in regulating gene expression in embryos.

Comparative analysis of function and interaction of transcription factors in nematodes: extensive conservation of orthology coupled to rapid sequence evolution

Background

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 Caenorhabditis elegans 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 Caenorhabditid 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.

Results

We identified 988 TF genes in C. elegans, and inferred corresponding sets in C. briggsae and C. remanei, 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 C. elegans and C. briggsae 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 Drosophila melanogaster, Mus musculus and Homo sapiens 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.

Conclusion

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.

Anti-fungal innate immunity in C. elegans is enhanced by evolutionary diversification of antimicrobial peptides

Encounters with pathogens provoke changes in gene transcription that are an integral part of host innate immune responses. In recent years, studies with invertebrate model organisms have given insights into the origin, function, and evolution of innate immunity. Here, we use genome-wide transcriptome analysis to characterize the consequence of natural fungal infection in Caenorhabditis elegans. We identify several families of genes encoding putative antimicrobial peptides (AMPs) and proteins that are transcriptionally up-regulated upon infection. Many are located in small genomic clusters. We focus on the nlp-29 cluster of six AMP genes and show that it enhances pathogen resistance in vivo. The same cluster has a different structure in two other Caenorhabditis species. A phylogenetic analysis indicates that the evolutionary diversification of this cluster, especially in cases of intra-genomic gene duplications, is driven by natural selection. We further show that upon osmotic stress, two genes of the nlp-29 cluster are strongly induced. In contrast to fungus-induced nlp expression, this response is independent of the p38 MAP kinase cascade. At the same time, both involve the epidermal GATA factor ELT-3. Our results suggest that selective pressure from pathogens influences intra-genomic diversification of AMPs and reveal an unexpected complexity in AMP regulation as part of the invertebrate innate immune response.

We are interested in how exactly the nematode Caenorhabditi elegans, widely used in biological research, defends itself against fungal infection. Like most animals, this worm responds to infection by switching on defense genes. We used DNA chips to measure the levels of all the worm's 20,000 genes and discovered new inducible defense genes. Many of them encode small proteins or peptides that can probably kill microbes. By looking in other nematode species, we saw that these antimicrobial peptide genes are evolving rapidly. This means that they could be important for the worms' survival in their natural environment. We also looked at how some of these genes are regulated and uncovered a sophisticated control mechanism involving a series of proteins called kinases that relay signals within cells. The genes we looked at are active in the worm's skin. Some of the antimicrobial peptide genes that we looked at are also switched on in the skin by high salt, but in this case, the regulation doesn't involve the same cascade of kinases. The responses to both infection and high salt do, however, require the same transcription factor (the protein that actually switches genes on), in this case called a GATA factor.

Genome evolution in Caenorhabditis

Since the completion of the Caenorhabditis elegans genome sequence 10 years ago, efforts of the large community of C. elegans geneticists have resulted in a high-quality annotation of the structures and sequence relatedness of nearly all the protein encoding and RNA genes. Based on increasingly accurate gene counts in other species, it now appears that C. elegans has more functional genes than most insects and approximately the same number as most mammals. In the last few years, draft genome sequences for several other nematodes have been published (C. briggsae and Brugia malayi) or publicly released (C. remanei, C. brenneri, C. japonica, Pristionchus pacificus, Trichinella spiralis and Haemonchus contortus). Comparisons of gene content within the phylum and to other phyla reveal complex patterns of genome evolution. These patterns include substantial numbers of genes conserved across all the major metazoan phyla (core metazoan genes) and many nematode-specific genes and gene families. Nematode-specific genes are located predominantly on autosomal arms, which also have higher recombination rates. It appears that evolutionary innovations occur mostly in these regions, probably facilitated by higher recombination. Few of these genes have gross phenotypes when knocked down by RNAi, suggesting that many of them function in specific aspects of nematode biology that were not tested, including chemosensation, pathogen response and xenobiotic detoxification.

Towards a mutation in every gene in Caenorhabditis elegans

The combined efforts of the Caenorhabditis elegans Knockout Consortium and individuals within the worm community are moving us closer to the goal of identifying mutations in every gene in the nematode C. elegans. At present, we count about 7000 deletion alleles that fall within 5500 genes. The principal method used to detect deletion mutations in the nematode utilizes polymerase chain reaction (PCR). More recently, the Moerman group has incorporated array comparative genome hybridization (aCGH) to detect deletions across the entire coding genome. Other methods used to detect mutant alleles in C. elegans include targeting induced local lesion in genomes (TILLING), transposon tagging, using either Tc1 or Mos1 and resequencing. These combined strategies have improved the overall throughput of the gene-knockout labs, and have broadened the types of mutations that we, and others, can identify. In this review, we will discuss these different approaches.

Transcription factor functionality and transcription regulatory networks

Now that numerous high-quality complete genome sequences are available, many efforts are focusing on the "second genomic code", namely the code that determines how the precise temporal and spatial expression of each gene in the genome is achieved. In this regard, the elucidation of transcription regulatory networks that describe combined transcriptional circuits for an organism of interest has become valuable to our understanding of gene expression at a systems level. Such networks describe physical and regulatory interactions between transcription factors (TFs) and the target genes they regulate under different developmental, physiological, or pathological conditions. The mapping of high-quality transcription regulatory networks depends not only on the accuracy of the experimental or computational method chosen, but also relies on the quality of TF predictions. Moreover, the total repertoire of TFs is not only determined by the protein-coding capacity of the genome, but also by different protein properties, including dimerization, co-factor interactions and post-translational modifications. Here, we discuss the factors that influence TF functionality and, hence, the functionality of the networks in which they operate.

Text-mining assisted regulatory annotation

Text-mining technologies can be integrated with genome annotation systems, increasing the availability of annotated cis-regulatory data.

Background

Decoding transcriptional regulatory networks and the genomic cis-regulatory logic implemented in their control nodes is a fundamental challenge in genome biology. High-throughput computational and experimental analyses of regulatory networks and sequences rely heavily on positive control data from prior small-scale experiments, but the vast majority of previously discovered regulatory data remains locked in the biomedical literature.

Results

We develop text-mining strategies to identify relevant publications and extract sequence information to assist the regulatory annotation process. Using a vector space model to identify Medline abstracts from papers likely to have high cis-regulatory content, we demonstrate that document relevance ranking can assist the curation of transcriptional regulatory networks and estimate that, minimally, 30,000 papers harbor unannotated cis-regulatory data. In addition, we show that DNA sequences can be extracted from primary text with high cis-regulatory content and mapped to genome sequences as a means of identifying the location, organism and target gene information that is critical to the cis-regulatory annotation process.

Conclusion

Our results demonstrate that text-mining technologies can be successfully integrated with genome annotation systems, thereby increasing the availability of annotated cis-regulatory data needed to catalyze advances in the field of gene regulation.

Matrix and Steiner-triple-system smart pooling assays for high-performance transcription regulatory network mapping

Yeast one-hybrid (Y1H) assays provide a gene-centered method for the identification of interactions between gene promoters and regulatory transcription factors (TFs). To date, Y1H assays have involved library screens that are relatively expensive and laborious. We present two Y1H strategies that allow immediate prey identification: matrix assays that use an array of 755 individual Caenorhabditis elegans TFs, and smart-pool assays that use TF multiplexing. Both strategies simplify the Y1H pipeline and reduce the cost of protein-DNA interaction identification. We used a Steiner triple system (STS) to create smart pools of 4-25 TFs. Notably, we uniplexed a small number of highly connected TFs to allow efficient assay deconvolution. Both strategies outperform library screens in terms of coverage, confidence and throughput. These versatile strategies can be adapted both to TFs in other systems and, likely, to other biomolecules and assays as well.

Transcription factor modularity in a gene-centered C. elegans core neuronal protein-DNA interaction network

Transcription regulatory networks play a pivotal role in the development, function, and pathology of metazoan organisms. Such networks are comprised of protein-DNA interactions between transcription factors (TFs) and their target genes. An important question pertains to how the architecture of such networks relates to network functionality. Here, we show that a Caenorhabditis elegans core neuronal protein-DNA interaction network is organized into two TF modules. These modules contain TFs that bind to a relatively small number of target genes and are more systems specific than the TF hubs that connect the modules. Each module relates to different functional aspects of the network. One module contains TFs involved in reproduction and target genes that are expressed in neurons as well as in other tissues. The second module is enriched for paired homeodomain TFs and connects to target genes that are often exclusively neuronal. We find that paired homeodomain TFs are specifically expressed in C. elegans and mouse neurons, indicating that the neuronal function of paired homeodomains is evolutionarily conserved. Taken together, we show that a core neuronal C. elegans protein-DNA interaction network possesses TF modules that relate to different functional aspects of the complete network.

PlanTAPDB, a phylogeny-based resource of plant transcription-associated proteins

Diversification of transcription-associated protein (TAP) families during land plant evolution is a key process yielding increased complexity of plant life. Understanding the evolutionary relationships between these genes is crucial to gain insight into plant evolution. We have determined a substantial set of TAPs that are focused on, but not limited to, land plants using PSI-BLAST searches and subsequent filtering and clustering steps. Phylogenies were created in an automated way using a combination of distance and maximum likelihood methods. Comparison of the data to previously published work confirmed their accuracy and usefulness for the majority of gene families. Evidence is presented that the flowering plant apical stem cell regulator WUSCHEL evolved from an ancestral homeobox gene that was already present after the water-to-land transition. The presence of distinct expanded gene families, such as COP1 and HIT in moss, is discussed within the evolutionary backdrop. Comparative analyses revealed that almost all angiosperm transcription factor families were already present in the earliest land plants, whereas many are missing among unicellular algae. A global analysis not only of transcription factors but also of transcriptional regulators and novel putative families is presented. A wealth of data about plant TAP families and all data accrued throughout their automated detection and analysis are made available via the PlanTAPDB Web interface. Evolutionary relationships of these genes are readily accessible to the nonexpert at a mouse-click. Initial analyses of selected gene families revealed that PlanTAPDB can easily be exerted for knowledge discovery.

Identification of muscle-specific regulatory modules in Caenorhabditis elegans

Transcriptional regulation is the major regulatory mechanism that controls the spatial and temporal expression of genes during development. This is carried out by transcription factors (TFs), which recognize and bind to their cognate binding sites. Recent studies suggest a modular organization of TF-binding sites, in which clusters of transcription-factor binding sites cooperate in the regulation of downstream gene expression. In this study, we report our computational identification and experimental verification of muscle-specific cis-regulatory modules in Caenorhabditis elegans. We first identified a set of motifs that are correlated with muscle-specific gene expression. We then predicted muscle-specific regulatory modules based on clusters of those motifs with characteristics similar to a collection of well-studied modules in other species. The method correctly identifies 88% of the experimentally characterized modules with a positive predictive value of at least 65%. The prediction accuracy of muscle-specific expression on an independent test set is highly significant (P<0.0001). We performed in vivo experimental tests of 12 predicted modules, and 10 of those drive muscle-specific gene expression. These results suggest that our method is highly accurate in identifying functional sequences important for muscle-specific gene expression and is a valuable tool for guiding experimental designs.

Insight into transcription factor gene duplication from Caenorhabditis elegans Promoterome-driven expression patterns

BACKGROUND

The C. elegans Promoterome is a powerful resource for revealing the regulatory mechanisms by which transcription is controlled pan-genomically. Transcription factors will form the core of any systems biology model of genome control and therefore the promoter activity of Promoterome inserts for C. elegans transcription factor genes was examined, in vivo, with a reporter gene approach.

RESULTS

Transgenic C. elegans strains were generated for 366 transcription factor promoter/gfp reporter gene fusions. GFP distributions were determined, and then summarized with reference to developmental stage and cell type. Reliability of these data was demonstrated by comparison to previously described gene product distributions. A detailed consideration of the results for one C. elegans transcription factor gene family, the Six family, comprising ceh-32, ceh-33, ceh-34 and unc-39 illustrates the value of these analyses. The high proportion of Promoterome reporter fusions that drove GFP expression, compared to previous studies, led to the hypothesis that transcription factor genes might be involved in local gene duplication events less frequently than other genes. Comparison of transcription factor genes of C. elegans and Caenorhabditis briggsae was therefore carried out and revealed very few examples of functional gene duplication since the divergence of these species for most, but not all, transcription factor gene families.

CONCLUSION

Examining reporter expression patterns for hundreds of promoters informs, and thereby improves, interpretation of this data type. Genes encoding transcription factors involved in intrinsic developmental control processes appear acutely sensitive to changes in gene dosage through local gene duplication, on an evolutionary time scale.

EDGEdb: a transcription factor-DNA interaction database for the analysis of C. elegans differential gene expression

Background

Transcription regulatory networks are composed of protein-DNA interactions between transcription factors and their target genes. A long-term goal in genome biology is to map protein-DNA interaction networks of all regulatory regions in a genome of interest. Both transcription factor -and gene-centered methods can be used to systematically identify such interactions. We use high-throughput yeast one-hybrid assays as a gene-centered method to identify protein-DNA interactions between regulatory sequences (e.g. gene promoters) and transcription factors in the nematode Caenorhabditis elegans. We have already mapped several hundred protein-DNA interactions and analyzed the transcriptional consequences of some by examining differential gene expression of targets in the presence or absence of an upstream regulator. The rapidly increasing amount of protein-DNA interaction data at a genome scale requires a database that facilitates efficient data storage, retrieval and integration.

Description

Here, we report the implementation of a C. elegans differential gene expression database (EDGEdb). This database enables the storage and retrieval of protein-DNA interactions and other data that relate to differential gene expression. Specifically, EDGEdb contains: i) sequence information of regulatory elements, including gene promoters, ii) sequence information of all 934 predicted transcription factors, their DNA binding domains, and, where available, their dimerization partners and consensus DNA binding sites, iii) protein-DNA interactions between regulatory elements and transcription factors, and iv) expression patterns conferred by regulatory elements, and how such patterns are affected by interacting transcription factors.

Conclusion

EDGEdb provides a protein-DNA -and protein-protein interaction resource for C. elegans transcription factors and a framework for similar databases for other organisms. The database is available at .

Parallel evolution of conserved non-coding elements that target a common set of developmental regulatory genes from worms to humans

BACKGROUND

The human genome contains thousands of non-coding sequences that are often more conserved between vertebrate species than protein-coding exons. These highly conserved non-coding elements (CNEs) are associated with genes that coordinate development, and have been proposed to act as transcriptional enhancers. Despite their extreme sequence conservation in vertebrates, sequences homologous to CNEs have not been identified in invertebrates.

RESULTS

Here we report that nematode genomes contain an alternative set of CNEs that share sequence characteristics, but not identity, with their vertebrate counterparts. CNEs thus represent a very unusual class of sequences that are extremely conserved within specific animal lineages yet are highly divergent between lineages. Nematode CNEs are also associated with developmental regulatory genes, and include well-characterized enhancers and transcription factor binding sites, supporting the proposed function of CNEs as cis-regulatory elements. Most remarkably, 40 of 156 human CNE-associated genes with invertebrate orthologs are also associated with CNEs in both worms and flies.

CONCLUSION

A core set of genes that regulate development is associated with CNEs across three animal groups (worms, flies and vertebrates). We propose that these CNEs reflect the parallel evolution of alternative enhancers for a common set of developmental regulatory genes in different animal groups. This 're-wiring' of gene regulatory networks containing key developmental coordinators was probably a driving force during the evolution of animal body plans. CNEs may, therefore, represent the genomic traces of these 'hard-wired' core gene regulatory networks that specify the development of each alternative animal body plan.

GATA transcription factor required for immunity to bacterial and fungal pathogens

In the past decade, Caenorhabditis elegans has been used to dissect several genetic pathways involved in immunity; however, little is known about transcription factors that regulate the expression of immune effectors. C. elegans does not appear to have a functional homolog of the key immune transcription factor NF-κB. Here we show that that the intestinal GATA transcription factor ELT-2 is required for both immunity to Salmonella enterica and expression of a C-type lectin gene, clec-67, which is expressed in the intestinal cells and is a good marker of S. enterica infection. We also found that ELT-2 is required for immunity to Pseudomonas aeruginosa, Enterococcus faecalis, and Cryptococcus neoformans. Lack of immune inhibition by DAF-2, which negatively regulates the FOXO transcription factor DAF-16, rescues the hypersusceptibility to pathogens phenotype of elt-2(RNAi) animals. Our results indicate that ELT-2 is part of a multi-pathogen defense pathway that regulates innate immunity independently of the DAF-2/DAF-16 signaling pathway.

Defining the transcriptional redundancy of early bodywall muscle development in C. elegans: evidence for a unified theory of animal muscle development

Myogenic regulatory factors (MRFs) are required for mammalian skeletal myogenesis. In contrast, bodywall muscle is readily detectable in Caenorhabditis elegans embryos lacking activity of the lone MRF ortholog HLH-1, indicating that additional myogenic factors must function in the nematode. We find that two additional C. elegans proteins, UNC-120/SRF and HND-1/HAND, can convert naïve blastomeres to muscle when overproduced ectopically in the embryo. In addition, we have used genetic null mutants to demonstrate that both of these factors act in concert with HLH-1 to regulate myogenesis. Loss of all three factors results in embryos that lack detectable bodywall muscle differentiation, identifying this trio as a set that is both necessary and sufficient for bodywall myogenesis in C. elegans. In mammals, SRF and HAND play prominent roles in regulating smooth and cardiac muscle development. That C. elegans bodywall muscle development is dependent on transcription factors that are associated with all three types of mammalian muscle supports a theory that all animal muscle types are derived from a common ancestral contractile cell type.

Unraveling transcription regulatory networks by protein-DNA and protein-protein interaction mapping

Metazoan genomes contain thousands of protein-coding and noncoding RNA genes, most of which are differentially expressed, i.e., at different locations or at different times during development, function, or pathology of the organism. Differential gene expression is achieved in part by the action of regulatory transcription factors (TFs) that bind to cis-regulatory elements that are often located in or near their target genes. Each TF likely regulates many targets in the context of intricate transcription regulatory networks. Up to 10% of a genome may encode TFs, but only a handful of these have been studied in detail. Here, I will discuss the different steps involved in the mapping and analysis of transcription regulatory networks, including the identification of network nodes (TFs and their target sequences) and edges (TF-TF dimers and TF-DNA target interactions), integration with other data types, and network properties and emerging principles that provide insights into differential gene expression.

Gateway-compatible yeast one-hybrid screens

INTRODUCTIONProtein-DNA interactions (PDIs) between transcription factors (TFs) and their target genes form the backbone of transcription regulatory networks. Such PDIs can be identified with either a TF or a gene as a starting point. The Gateway-compatible yeast one-hybrid (Y1H) system provides a high-throughput, gene-centered method for the identification of interactions between a "DNA bait" (e.g., cis-regulatory DNA elements or gene promoters) and "protein preys" (e.g., TFs). The Y1H system is a genetic system to detect PDIs based on selection of reporter gene expression in yeast. DNA baits are fused by Gateway cloning to two reporter genes, HIS3 and lacZ, and the resulting DNA bait::reporter constructs are subsequently integrated into the genome of the host yeast strain. After integration, baits are examined for self-activation (i.e., their ability to drive reporter gene expression in the absence of an exogenous prey protein). Subsequently, each DNA bait is screened for interacting proteins by transforming a library of preys into the corresponding Y1H DNA bait yeast strain. Preys are hybrid proteins composed of a protein from the organism of interest and a heterologous transcription activation domain. When a prey protein binds to the DNA bait, the heterologous activation domain activates reporter gene expression. Thus, physical interactions between both repressors and activators and their DNA targets can be identified.

The C2H2 zinc finger genes of Strongylocentrotus purpuratus and their expression in embryonic development

The C2H2 zinc finger is one of the most abundant protein domains and is thought to have been extensively replicated in diverse animal clades. Some well-studied proteins that contain this domain are transcriptional regulators. As part of an attempt to delineate all transcription factors encoded in the Strongylocentrotus purpuratus genome, we identified the C2H2 zinc finger genes indicated in the sequence, and examined their involvement in embryonic development. We found 377 zinc finger genes in the sea urchin genome, about half the number found in mice or humans. Their expression was measured by quantitative PCR. Up to the end of gastrulation less than a third of these genes is expressed, and about 75% of the expressed genes are maternal; both parameters distinguish these from all other classes of regulatory genes as measured in other studies. Spatial expression pattern was determined by whole mount in situ hybridization for 43 genes transcribed at a sufficient level, and localized expression was observed in diverse embryonic tissues. These genes may execute important regulatory functions in development. However, the functional meaning of the majority of this large gene family remains undefined.

A gene-centered C. elegans protein-DNA interaction network

Transcription regulatory networks consist of physical and functional interactions between transcription factors (TFs) and their target genes. The systematic mapping of TF-target gene interactions has been pioneered in unicellular systems, using "TF-centered" methods (e.g., chromatin immunoprecipitation). However, metazoan systems are less amenable to such methods. Here, we used "gene-centered" high-throughput yeast one-hybrid (Y1H) assays to identify 283 interactions between 72 C. elegans digestive tract gene promoters and 117 proteins. The resulting protein-DNA interaction (PDI) network is highly connected and enriched for TFs that are expressed in the digestive tract. We provide functional annotations for approximately 10% of all worm TFs, many of which were previously uncharacterized, and find ten novel putative TFs, illustrating the power of a gene-centered approach. We provide additional in vivo evidence for multiple PDIs and illustrate how the PDI network provides insights into metazoan differential gene expression at a systems level.

Specification of the C. elegans MS blastomere by the T-box factor TBX-35

In C. elegans, many mesodermal cell types are made by descendants of the progenitor MS, born at the seven-cell stage of embryonic development. Descendants of MS contribute to body wall muscle and to the posterior half of the pharynx. We have previously shown that MS is specified by the activity of the divergent MED-1,2 GATA factors. We report that the MED-1,2 target gene tbx-35, which encodes a T-box transcription factor, specifies the MS fate. Embryos homozygous for a putative tbx-35-null mutation fail to generate MS-derived pharynx and body muscle, and instead generate ectopic PAL-1-dependent muscle and hypodermis, tissues normally made by the C blastomere. Conversely, overexpression of tbx-35 results in the generation of ectopic pharynx and muscle tissue. The MS and E sister cells are made different by transduction of a Wnt/MAPK/Src pathway signal through the nuclear effector TCF/POP-1. We show that in E, tbx-35 is repressed in a Wnt-dependent manner that does not require activity of TCF/POP-1, suggesting that an additional nuclear Wnt effector functions in E to repress MS development. Genes of the T-box family are known to function in protostomes and deuterostomes in the specification of mesodermal fates. Our results show that this role has been evolutionarily conserved in the early C. elegans embryo, and that a progenitor of multiple tissue types can be specified by a surprisingly simple gene cascade.

Genomic clusters, putative pathogen recognition molecules, and antimicrobial genes are induced by infection of C. elegans with M. nematophilum

The interaction between the nematode Caenorhabditis elegans and a Gram-positive bacterial pathogen, Microbacterium nematophilum, provides a model for an innate immune response in nematodes. This pathogen adheres to the rectal and post-anal cuticle of the worm, causing slowed growth, constipation, and a defensive swelling response of rectal hypodermal cells. To explore the genomic responses that the worm activates after pathogenic attack we used microarray analysis of transcriptional changes induced after 6-h infection, comparing virulent with avirulent infection. We defined 89 genes with statistically significant expression changes of at least twofold, of which 68 were up-regulated and 21 were down-regulated. Among the former, those encoding C-type lectin domains were the most abundant class. Many of the 89 genes exhibit genomic clustering, and we identified one large cluster of 62 genes, of which most were induced in response to infection. We tested 41 of the induced genes for involvement in immunity using mutants or RNAi, finding that six of these are required for the swelling response and five are required more generally for defense. Our results indicate that C-type lectins and other putative pathogen-recognition molecules are important for innate immune defense in C. elegans. We also found significant induction of genes encoding lysozymes, proteases, and defense-related proteins, as well as various domains of unknown function. The genes induced during infection by M. nematophilum appear largely distinct from genes induced by other pathogens, suggesting that C. elegans mounts pathogen-specific responses to infection.