Toward functional genomics in bacteria: Analysis of gene expression in Escherichia coli from a bacterial artificial chromosome library of Bacillus cereus

As the study of microbes moves into the era of functional genomics, there is an increasing need for molecular tools for analysis of a wide diversity of microorganisms. Currently, biological study of many prokaryotes of agricultural, medical, and fundamental scientific interest is limited by the lack of adequate genetic tools. We report the application of the bacterial artificial chromosome (BAC) vector to prokaryotic biology as a powerful approach to address this need. We constructed a BAC library in Escherichia coli from genomic DNA of the Gram-positive bacterium Bacillus cereus. This library provides 5.75-fold coverage of the B. cereus genome, with an average insert size of 98 kb. To determine the extent of heterologous expression of B. cereus genes in the library, we screened it for expression of several B. cereus activities in the E. coli host. Clones expressing 6 of 10 activities tested were identified in the library, namely, ampicillin resistance, zwittermicin A resistance, esculin hydrolysis, hemolysis, orange pigment production, and lecithinase activity. We analyzed selected BAC clones genetically to identify rapidly specific B. cereus loci. These results suggest that BAC libraries will provide a powerful approach for studying gene expression from diverse prokaryotes.

Identification of eIF2Bγ and eIF2γ as cofactors of hepatitis C virus internal ribosome entry site-mediated translation using a functional genomics approach

The 5′-untranslated region of hepatitis C virus (HCV) is highly conserved, folds into a complex secondary structure, and functions as an internal ribosome entry site (IRES) to initiate translation of HCV proteins. We have developed a selection system based on a randomized hairpin ribozyme gene library to identify cellular factors involved in HCV IRES function. A retroviral vector ribozyme library with randomized target recognition sequences was introduced into HeLa cells, stably expressing a bicistronic construct encoding the hygromycin B phosphotransferase gene and the herpes simplex virus thymidine kinase gene (HSV-tk). Translation of the HSV-tk gene was mediated by the HCV IRES. Cells expressing ribozymes that inhibit HCV IRES-mediated translation of HSV-tk were selected via their resistance to both ganciclovir and hygromycin B. Two ribozymes reproducibly conferred the ganciclovir-resistant phenotype and were shown to inhibit IRES-mediated translation of HCV core protein but did not inhibit cap-dependent protein translation or cell growth. The functional targets of these ribozymes were identified as the gamma subunits of human eukaryotic initiation factors 2B (eIF2Bγ) and 2 (eIF2γ), respectively. The involvement of eIF2Bγ and eIF2γ in HCV IRES-mediated translation was further validated by ribozymes directed against additional sites within the mRNAs of these genes. In addition to leading to the identification of cellular IRES cofactors, ribozymes obtained from this cellular selection system could be directly used to specifically inhibit HCV viral translation, thereby facilitating the development of new antiviral strategies for HCV infection.

PlantsP: a functional genomics database for plant phosphorylation

The PlantsP database is a curated database that combines information derived from sequences with experimental functional genomics information. PlantsP focuses on plant protein kinases and protein phosphatases. The database will specifically provide a resource for information on a collection of T-DNA insertion mutants (knockouts) in each protein kinase and phosphatase in Arabidopsis thaliana. PlantsP also provides a curated view of each protein that includes a comprehensive annotation of functionally related sequence motifs, sequence family definitions, alignments and phylogenetic trees, and descriptive information drawn directly from the literature. PlantsP is available at http://PlantsP.sdsc.edu.

Functional genomics

Complete genome sequences are providing a framework to allow the investigation of biological processes by the use of comprehensive approaches. Genome analysis also is having a dramatic impact on medicine through its identification of genes and mutations involved in disease and the elucidation of entire microbial gene sets. Studies of the sequences of model organisms, such as that of the nematode worm Caenorhabditis elegans, are providing extraordinary insights into development and differentiation that aid the study of these processes in humans. The field of functional genomics seeks to devise and apply technologies that take advantage of the growing body of sequence information to analyze the full complement of genes and proteins encoded by an organism.

Functional genomics: Probing plant gene function and expression with transposons

Transposable elements provide a convenient and flexible means to disrupt plant genes, so allowing their function to be assessed. By engineering transposons to carry reporter genes and regulatory signals, the expression of target genes can be monitored and to some extent manipulated. Two strategies for using transposons to assess gene function are outlined here: First, the PCR can be used to identify plants that carry insertions into specific genes from among pools of heavily mutagenized individuals (site-selected transposon mutagenesis). This method requires that high copy transposons be used and that a relatively large number of reactions be performed to identify insertions into genes of interest. Second, a large library of plants, each carrying a unique insertion, can be generated. Each insertion site then can be amplified and sequenced systematically. These two methods have been demonstrated in maize, Arabidopsis, and other plant species, and the relative merits of each are discussed in the context of plant genome research.

Toward controlling gene expression at will: Specific regulation of the erbB-2/HER-2 promoter by using polydactyl zinc finger proteins constructed from modular building blocks

To create a universal system for the control of gene expression, we have studied methods for the construction of novel polydactyl zinc finger proteins that recognize extended DNA sequences. Elsewhere we have described the generation of zinc finger domains recognizing sequences of the 5′-GNN-3′ subset of a 64-member zinc finger alphabet. Here we report on the use of these domains as modular building blocks for the construction of polydactyl proteins specifically recognizing 9- or 18-bp sequences. A rapid PCR assembly method was developed that, together with this predefined set of zinc finger domains, provides ready access to 17 million novel proteins that bind the 5′-(GNN)6-3′ family of 18-bp DNA sites. To examine the efficacy of this strategy in gene control, the human erbB-2 gene was chosen as a model. A polydactyl protein specifically recognizing an 18-bp sequence in the 5′-untranslated region of this gene was converted into a transcriptional repressor by fusion with Krüppel-associated box (KRAB), ERD, or SID repressor domains. Transcriptional activators were generated by fusion with the herpes simplex VP16 activation domain or with a tetrameric repeat of VP16’s minimal activation domain, termed VP64. We demonstrate that both gene repression and activation can be achieved by targeting designed proteins to a single site within the transcribed region of a gene. We anticipate that gene-specific transcriptional regulators of the type described here will find diverse applications in gene therapy, functional genomics, and the generation of transgenic organisms.

Genetic footprinting with mariner-based transposition in Pseudomonas aeruginosa

The complete DNA sequence of Pseudomonas aeruginosa provides an opportunity to apply functional genomics to a major human pathogen. A comparative genomics approach combined with genetic footprinting was used as a strategy to identify genes required for viability in P. aeruginosa. Use of a highly efficient in vivo mariner transposition system in P. aeruginosa facilitated the analysis of candidate genes of this class. We have developed a rapid and efficient allelic exchange system by using the I-SceI homing endonuclease in conjunction with in vitro mariner mutagenesis to generate mutants within targeted regions of the P. aeruginosa chromosome for genetic footprinting analyses. This technique for generating transposon insertion mutants should be widely applicable to other organisms that are not naturally transformable or may lack well developed in vivo transposition systems. We tested this system with three genes in P. aeruginosa that have putative essential homologs in Haemophilus influenzae. We show that one of three H. influenzae essential gene homologs is needed for growth in P. aeruginosa, validating the practicality of this comparative genomics strategy to identify essential genes in P. aeruginosa.

Chromosomal transposition of a Tc1/mariner-like element in mouse embryonic stem cells

Mouse has become an increasingly important organism for modeling human diseases and for determining gene function in a mammalian context. Unfortunately, transposon-tagged mutagenesis, one of the most valuable tools for functional genomics, still is not available in this organism. On the other hand, it has long been speculated that members of the Tc1/mariner-like elements may be less dependent on host factors and, hence, can be introduced into heterologous organisms. However, this prediction has not been realized in mice. We report here the chromosomal transposition of the Sleeping Beauty (SB) element in mouse embryonic stem cells, providing evidence that it can be used as an in vivo mutagen in mice.

Statistical prediction of single-stranded regions in RNA secondary structure and application to predicting effective antisense target sites and beyond

Single-stranded regions in RNA secondary structure are important for RNA–RNA and RNA–protein interactions. We present a probability profile approach for the prediction of these regions based on a statistical algorithm for sampling RNA secondary structures. For the prediction of phylogenetically-determined single-stranded regions in secondary structures of representative RNA sequences, the probability profile offers substantial improvement over the minimum free energy structure. In designing antisense oligonucleotides, a practical problem is how to select a secondary structure for the target mRNA from the optimal structure(s) and many suboptimal structures with similar free energies. By summarizing the information from a statistical sample of probable secondary structures in a single plot, the probability profile not only presents a solution to this dilemma, but also reveals ‘well-determined’ single-stranded regions through the assignment of probabilities as measures of confidence in predictions. In antisense application to the rabbit β-globin mRNA, a significant correlation between hybridization potential predicted by the probability profile and the degree of inhibition of in vitro translation suggests that the probability profile approach is valuable for the identification of effective antisense target sites. Coupling computational design with DNA–RNA array technique provides a rational, efficient framework for antisense oligonucleotide screening. This framework has the potential for high-throughput applications to functional genomics and drug target validation.

The Medicago Genome Initiative: a model legume database

The Medicago Genome Initiative (MGI) is a database of EST sequences of the model legume Medicago truncatula. The database is available to the public and has resulted from a collaborative research effort between the Samuel Roberts Noble Foundation and the National Center for Genome Resources to investigate the genome of M.truncatula. MGI is part of the greater integrated Medicago functional genomics program at the Noble Foundation (http://www.noble .org), which is taking a global approach in studying the genetic and biochemical events associated with the growth, development and environmental interactions of this model legume. Our approach will include: large-scale EST sequencing, gene expression profiling, the generation of M.truncatula activation-tagged and promoter trap insertion mutants, high-throughput metabolic profiling, and proteome studies. These multidisciplinary information pools will be interfaced with one another to provide scientists with an integrated, holistic set of tools to address fundamental questions pertaining to legume biology. The public interface to the MGI database can be accessed at http://www.ncgr.org/research/mgi.

A rapid classification protocol for the CATH Domain Database to support structural genomics

In order to support the structural genomic initiatives, both by rapidly classifying newly determined structures and by suggesting suitable targets for structure determination, we have recently developed several new protocols for classifying structures in the CATH domain database (http://www.biochem.ucl.ac.uk/bsm/cath). These aim to increase the speed of classification of new structures using fast algorithms for structure comparison (GRATH) and to improve the sensitivity in recognising distant structural relatives by incorporating sequence information from relatives in the genomes (DomainFinder). In order to ensure the integrity of the database given the expected increase in data, the CATH Protein Family Database (CATH-PFDB), which currently includes 25 320 structural domains and a further 160 000 sequence relatives has now been installed in a relational ORACLE database. This was essential for developing more rigorous validation procedures and for allowing efficient querying of the database, particularly for genome analysis. The associated Dictionary of Homologous Superfamilies [Bray,J.E., Todd,A.E., Pearl,F.M.G., Thornton,J.M. and Orengo,C.A. (2000) Protein Eng., 13, 153–165], which provides multiple structural alignments and functional information to assist in assigning new relatives, has also been expanded recently and now includes information for 903 homo­logous superfamilies. In order to improve coverage of known structures, preliminary classification levels are now provided for new structures at interim stages in the classification protocol. Since a large proportion of new structures can be rapidly classified using profile-based sequence analysis [e.g. PSI-BLAST: Altschul,S.F., Madden,T.L., Schaffer,A.A., Zhang,J., Zhang,Z., Miller,W. and Lipman,D.J. (1997) Nucleic Acids Res., 25, 3389–3402], this provides preliminary classification for easily recognisable homologues, which in the latest release of CATH (version 1.7) represented nearly three-quarters of the non-identical structures.

Comparative genomics of the restriction-modification systems in Helicobacter pylori

Helicobacter pylori is a Gram-negative bacterial pathogen with a small genome of 1.64–1.67 Mb. More than 20 putative DNA restriction-modification (R-M) systems, comprising more than 4% of the total genome, have been identified in the two completely sequenced H. pylori strains, 26695 and J99, based on sequence similarities. In this study, we have investigated the biochemical activities of 14 Type II R-M systems in H. pylori 26695. Less than 30% of the Type II R-M systems in 26695 are fully functional, similar to the results obtained from strain J99. Although nearly 90% of the R-M genes are shared by the two H. pylori strains, different sets of these R-M genes are functionally active in each strain. Interestingly, all strain-specific R-M genes are active, whereas most shared genes are inactive. This agrees with the notion that strain-specific genes have been acquired more recently through horizontal transfer from other bacteria and selected for function. Thus, they are less likely to be impaired by random mutations. Our results also show that H. pylori has extremely diversified R-M systems in different strains, and that the diversity may be maintained by constantly acquiring new R-M systems and by inactivating and deleting the old ones.

Genes linked by fusion events are generally of the same functional category: A systematic analysis of 30 microbial genomes

Recent work in computational genomics has shown that a functional association between two genes can be derived from the existence of a fusion of the two as one continuous sequence in another genome. For each of 30 completely sequenced microbial genomes, we established all such fusion links among its genes and determined the distribution of links within and among 15 broad functional categories. We found that 72% of all fusion links related genes of the same functional category. A comparison of the distribution of links to simulations on the basis of a random model further confirmed the significance of intracategory fusion links. Where a gene of annotated function is linked to an unclassified gene, the fusion link suggests that the two genes belong to the same functional category. The predictions based on fusion links are shown here for Methanobacterium thermoautotrophicum, and another 661 predictions are available at http://fusion.bu.edu.