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

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.

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.

Sequence variation in the human transcription factor gene POU5F1.

BACKGROUND: POU5F1 expression is required to maintain stem cell pluripotency and for primordial germ cells to retain proliferative capability in embryonic development. Recent evidence suggests that POU5F1 may also be a testicular germ cell carcinoma (TGCC) oncogene, and POU5F1 variation may influence TGCC risk. As an important first step to a genetic association study, we sought to identify all common sequence variants in an 11.3 kb region containing POU5F1, and to describe the linkage disequilibrium patterns, using DNA from individuals of African-descent (AD) and European-descent (ED). RESULTS: A higher number of polymorphisms was observed in the AD (n = 102) versus ED (n = 82) population. Among the 41 observed haplotypes, 21 (51%) and 12 (29%) were unique to the AD and ED populations, respectively, while 8 (20%) were observed in both. The number of tagging polymorphisms necessary to explain at least 80% of common variation (minor allele frequency > or = 0.10) due to the remaining untyped polymorphisms was 17 for an AD and 10 for an ED population, providing a 4.0- and 7.0-fold gain in genotyping efficiency for characterizing nucleotide variation, respectively. CONCLUSION: POU5F1 is highly polymorphic, however a smaller subset of polymorphisms can tag the observed genetic variation with little loss of information.

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.

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.

Towards a proteome-scale map of the human protein-protein interaction network.

Systematic mapping of protein-protein interactions, or 'interactome' mapping, was initiated in model organisms, starting with defined biological processes and then expanding to the scale of the proteome. Although far from complete, such maps have revealed global topological and dynamic features of interactome networks that relate to known biological properties, suggesting that a human interactome map will provide insight into development and disease mechanisms at a systems level. Here we describe an initial version of a proteome-scale map of human binary protein-protein interactions. Using a stringent, high-throughput yeast two-hybrid system, we tested pairwise interactions among the products of approximately 8,100 currently available Gateway-cloned open reading frames and detected approximately 2,800 interactions. This data set, called CCSB-HI1, has a verification rate of approximately 78% as revealed by an independent co-affinity purification assay, and correlates significantly with other biological attributes. The CCSB-HI1 data set increases by approximately 70% the set of available binary interactions within the tested space and reveals more than 300 new connections to over 100 disease-associated proteins. This work represents an important step towards a systematic and comprehensive human interactome project.

C. elegans ORFeome version 3.1: increasing the coverage of ORFeome resources with improved gene predictions.

The first version of the Caenorhabditis elegans ORFeome cloning project, based on release WS9 of Wormbase (August 1999), provided experimental verifications for approximately 55% of predicted protein-encoding open reading frames (ORFs). The remaining 45% of predicted ORFs could not be cloned, possibly as a result of mispredicted gene boundaries. Since the release of WS9, gene predictions have improved continuously. To test the accuracy of evolving predictions, we attempted to PCR-amplify from a highly representative worm cDNA library and Gateway-clone approximately 4200 ORFs missed earlier and for which new predictions are available in WS100 (May 2003). In this set we successfully cloned 63% of ORFs with supporting experimental data ("touched" ORFs), and 42% of ORFs with no supporting experimental evidence ("untouched" ORFs). Approximately 2000 full-length ORFs were cloned in-frame, 13% of which were corrected in their exon/intron structure relative to WS100 predictions. In total, approximately 12,500 C. elegans ORFs are now available as Gateway Entry clones for various reverse proteomics (ORFeome v3.1). This work illustrates why the cloning of a complete C. elegans ORFeome, and likely the ORFeomes of other multicellular organisms, needs to be an iterative process that requires multiple rounds of experimental validation together with gradually improving gene predictions.

Human ORFeome version 1.1: a platform for reverse proteomics.

The advent of systems biology necessitates the cloning of nearly entire sets of protein-encoding open reading frames (ORFs), or ORFeomes, to allow functional studies of the corresponding proteomes. Here, we describe the generation of a first version of the human ORFeome using a newly improved Gateway recombinational cloning approach. Using the Mammalian Gene Collection (MGC) resource as a starting point, we report the successful cloning of 8076 human ORFs, representing at least 7263 human genes, as mini-pools of PCR-amplified products. These were assembled into the human ORFeome version 1.1 (hORFeome v1.1) collection. After assessing the overall quality of this version, we describe the use of hORFeome v1.1 for heterologous protein expression in two different expression systems at proteome scale. The hORFeome v1.1 represents a central resource for the cloning of large sets of human ORFs in various settings for functional proteomics of many types, and will serve as the foundation for subsequent improved versions of the human ORFeome.

A first version of the Caenorhabditis elegans Promoterome.

An important aspect of the development of systems biology approaches in metazoans is the characterization of expression patterns of nearly all genes predicted from genome sequences. Such "localizome" maps should provide information on where (in what cells or tissues) and when (at what stage of development or under what conditions) genes are expressed. They should also indicate in what cellular compartments the corresponding proteins are localized. Caenorhabditis elegans is particularly suited for the development of a localizome map since all its 959 adult somatic cells can be visualized by microscopy, and its cell lineage has been completely described. Here we address one of the challenges of C. elegans localizome mapping projects: that of obtaining a genome-wide resource of C. elegans promoters needed to generate transgenic animals expressing localization markers such as the green fluorescent protein (GFP). To ensure high flexibility for future uses, we utilized the newly developed MultiSite Gateway system. We generated and validated "version 1.1" of the Promoterome: a resource of approximately 6000 C. elegans promoters. These promoters can be transferred easily into various Gateway Destination vectors to drive expression of markers such as GFP, alone (promoter::GFP constructs), or in fusion with protein-encoding open reading frames available in ORFeome resources (promoter::ORF::GFP).

Generation of the Brucella melitensis ORFeome version 1.1.

The bacteria of the Brucella genus are responsible for a worldwide zoonosis called brucellosis. They belong to the alpha-proteobacteria group, as many other bacteria that live in close association with a eukaryotic host. Importantly, the Brucellae are mainly intracellular pathogens, and the molecular mechanisms of their virulence are still poorly understood. Using the complete genome sequence of Brucella melitensis, we generated a database of protein-coding open reading frames (ORFs) and constructed an ORFeome library of 3091 Gateway Entry clones, each containing a defined ORF. This first version of the Brucella ORFeome (v1.1) provides the coding sequences in a user-friendly format amenable to high-throughput functional genomic and proteomic experiments, as the ORFs are conveniently transferable from the Entry clones to various Expression vectors by recombinational cloning. The cloning of the Brucella ORFeome v1.1 should help to provide a better understanding of the molecular mechanisms of virulence, including the identification of bacterial protein-protein interactions, but also interactions between bacterial effectors and their host's targets.

A map of the interactome network of the metazoan C. elegans.

To initiate studies on how protein-protein interaction (or "interactome") networks relate to multicellular functions, we have mapped a large fraction of the Caenorhabditis elegans interactome network. Starting with a subset of metazoan-specific proteins, more than 4000 interactions were identified from high-throughput, yeast two-hybrid (HT=Y2H) screens. Independent coaffinity purification assays experimentally validated the overall quality of this Y2H data set. Together with already described Y2H interactions and interologs predicted in silico, the current version of the Worm Interactome (WI5) map contains approximately 5500 interactions. Topological and biological features of this interactome network, as well as its integration with phenome and transcriptome data sets, lead to numerous biological hypotheses.