Web-based visualization tools for bacterial genome alignments

Genome sequence of Cronobacter sakazakii BAA-894 and comparative genomic hybridization analysis with other Cronobacter species

BACKGROUND

The genus Cronobacter (formerly called Enterobacter sakazakii) is composed of five species; C. sakazakii, C. malonaticus, C. turicensis, C. muytjensii, and C. dublinensis. The genus includes opportunistic human pathogens, and the first three species have been associated with neonatal infections. The most severe diseases are caused in neonates and include fatal necrotizing enterocolitis and meningitis. The genetic basis of the diversity within the genus is unknown, and few virulence traits have been identified.

METHODOLOGY/PRINCIPAL FINDINGS

We report here the first sequence of a member of this genus, C. sakazakii strain BAA-894. The genome of Cronobacter sakazakii strain BAA-894 comprises a 4.4 Mb chromosome (57% GC content) and two plasmids; 31 kb (51% GC) and 131 kb (56% GC). The genome was used to construct a 387,000 probe oligonucleotide tiling DNA microarray covering the whole genome. Comparative genomic hybridization (CGH) was undertaken on five other C. sakazakii strains, and representatives of the four other Cronobacter species. Among 4,382 annotated genes inspected in this study, about 55% of genes were common to all C. sakazakii strains and 43% were common to all Cronobacter strains, with 10-17% absence of genes.

CONCLUSIONS/SIGNIFICANCE

CGH highlighted 15 clusters of genes in C. sakazakii BAA-894 that were divergent or absent in more than half of the tested strains; six of these are of probable prophage origin. Putative virulence factors were identified in these prophage and in other variable regions. A number of genes unique to Cronobacter species associated with neonatal infections (C. sakazakii, C. malonaticus and C. turicensis) were identified. These included a copper and silver resistance system known to be linked to invasion of the blood-brain barrier by neonatal meningitic strains of Escherichia coli. In addition, genes encoding for multidrug efflux pumps and adhesins were identified that were unique to C. sakazakii strains from outbreaks in neonatal intensive care units.

PipMaker: a World Wide Web server for genomic sequence alignments

PipMaker is a World-Wide Web site used to compare two long genomic sequences and identify conserved segments between them. This unit describes the use of the PipMaker server and explains the resulting output files. PipMaker provides an efficient method of aligning genomic sequences and returns a compact, but easy-to-interpret form of output, the percent identity plot (pip). For each aligning segment between two sequences the pip shows both the position relative to the first sequence and the degree of similarity. Optional annotations on the pip provide additional information to assist in the interpretation of the alignment. The default parameters of the underlying blastz alignment program are tuned for human-mouse alignments.

GenomeBlast: a web tool for small genome comparison

Background

Comparative genomics has become an essential approach for identifying homologous gene candidates and their functions, and for studying genome evolution. There are many tools available for genome comparisons. Unfortunately, most of them are not applicable for the identification of unique genes and the inference of phylogenetic relationships in a given set of genomes.

Results

GenomeBlast is a Web tool developed for comparative analysis of multiple small genomes. A new parameter called "coverage" was introduced and used along with sequence identity to evaluate global similarity between genes. With GenomeBlast, the following results can be obtained: (1) unique genes in each genome; (2) homologous gene candidates among compared genomes; (3) 2D plots of homologous gene candidates along the all pairwise genome comparisons; and (4) a table of gene presence/absence information and a genome phylogeny. We demonstrated the functions in GenomeBlast with an example of multiple herpesviral genome analysis and illustrated how GenomeBlast is useful for small genome comparison.

Conclusion

We developed a Web tool for comparative analysis of small genomes, which allows the user not only to identify unique genes and homologous gene candidates among multiple genomes, but also to view their graphical distributions on genomes, and to reconstruct genome phylogeny. GenomeBlast runs on a Linux server with 4 CPUs and 4 GB memory. The online version of GenomeBlast is available to public by using a Web browser with the URL .

Multistress regulation in Escherichia coli: expression of osmB involves two independent promoters responding either to sigmaS or to the RcsCDB His-Asp phosphorelay

Transcription of the Escherichia coli osmB gene is induced by several stress conditions. osmB is expressed from two promoters, osmBp1 and osmBp2. The downstream promoter, osmBp2, is induced after osmotic shock or upon entry into stationary phase in a sigma(S)-dependent manner. The upstream promoter, osmBp1, is independent of sigma(S) and is activated by RcsB, the response regulator of the His-Asp phosphorelay signal transduction system RcsCDB. RcsB is responsible for the induction of osmBp1 following treatment with chlorpromazine. Activation of osmBp1 by RcsB requires a sequence upstream of its -35 element similar to the RcsB binding site consensus, suggesting a direct regulatory role. osmB appears as another example of a multistress-responsive gene whose transcription involves both a sigma(S)-dependent promoter and a second one independent of sigma(S) but controlled by stress-specific transcription factors.

Multispecies comparative analysis of a mammalian-specific genomic domain encoding secretory proteins

The mammalian-specific casein gene cluster comprises 3 or 4 evolutionarily related genes and 1 physically linked gene with a functional association. To gain a better understanding of the mechanisms regulating the entire casein cluster at the genomic level we initiated a multispecies comparative sequence analysis. Despite the high level of divergence at the coding level, these studies have identified uncharacterized family members within two species and the presence at orthologous positions of previously uncharacterized genes. Also the previous suggestion that the histatin/statherin gene family, located in this region, was primate specific was ruled out. All 11 genes identified in this region appear to encode secretory proteins. Conservation of a number of noncoding regions was observed; one coincides with an element previously suggested to be important for beta-casein gene expression in human and cow. The conserved regions might have biological importance for the regulation of genes in this genomic "neighborhood."

Cassette-like variation of restriction enzyme genes in Escherichia coli C and relatives

A surprising result of comparative bacterial genomics has been the large amount of DNA found to be present in one strain but not in another of the same species. We examine in detail one location where gene content varies extensively, the restriction cluster in Escherichia coli. This region is designated the Immigration Control Region (ICR) for the density and variability of restriction functions found there. To better define the boundaries of this variable locus, we determined the sequence of the region from a restrictionless strain, E.coli C. Here we compare the 13.7 kb E.coli C sequence spanning the site of the ICR with corresponding sequences from five E.coli strains and Salmonella typhimurium LT2. To discuss this variation, we adopt the term 'framework' to refer to genes that are stable components of genomes within related lineages, while 'migratory' genes are transient inhabitants of the genome. Strikingly, seven different migratory DNA segments, encoding different sets of genes and gene fragments, alternatively occupy a single well-defined location in the seven strains examined. The flanking framework genes, yjiS and yjiA, display approximately normal patterns of conservation. The patterns observed are consistent with the action of a site-specific recombinase. Since no nearby gene codes for a likely recombinase of known families, such a recombinase must be of a new family or unlinked.

Comparative genomics

Comparing the genomes of two different species allow the exploration of a host of intriguing evolutionary and genetic questions

MGView: an alignment and visualization tool to enhance gap closure of microbial genomes

Gap closure is a challenging phase in microbial random shotgun genome sequencing projects, particularly since genome assemblies are often complicated by the presence of repeat elements, insertion sequences and other similar factors that contribute to sequence misassemblies. While it is well recognized that the conservation of genetic information between microbial genomes, combined with the exponential increase in available microbial sequences, can be exploited to increase the efficiency of gap closure, we lack the computational tools to aid in this process. We describe here a new tool, MGView, which was developed to create a graphical depiction of the alignment of a set of microbial contigs against a completed microbial genome. The results of our assembly of the Staphylococcus aureus RF122 genome show that MGView enables a considerable reduction in time and economic cost associated with closure. Together, the results also show that the application of MGView not only enables a reduction in fold-coverage requirements of the random shotgun sequence phase, but also provides interesting insights into differences in gene content and organization between finished and unfinished microbial genomes.

GenomeComp: a visualization tool for microbial genome comparison

We have developed a software tool, GenomeComp, for summarizing, parsing and visualizing the genome sequences comparison results derived from voluminous BLAST textual output. With GenomeComp, the variation between genomes can be easily highlighted, such as repeat regions, insertions, deletions and rearrangements of genomic segments. This software provides a new visualizing tool for microbe comparative genomics.

Lateral gene transfer in Salmonella

Comparative genomics and microarrays reveal that the genomes of different Salmonella enterica serovars are distinguished from each other by the presence or absence of hundreds of genes. The distribution of these variable genome regions is often not clonal. Therefore, lateral gene transfer (LGT) plays an important role in diversity among Salmonella. Overall, almost one quarter of the entire S. enterica sv Typhimurium genome may have been introduced by LGT.

NhaR and RcsB independently regulate the osmCp1 promoter of Escherichia coli at overlapping regulatory sites

Transcription of the Escherichia coli osmC gene is induced by several stress conditions. osmC is expressed from two overlapping promoters, osmCp1 and osmCp2. The proximal promoter, osmCp2, is transcribed at the entry into the stationary phase by the sigma(s) sigma factor. The distal promoter, osmCp1, is activated by NhaR and RcsB. NhaR is a positive regulator of the LysR family and is known to be an activator of the nhaA gene encoding an Na(+)/H(+) antiporter. RcsB is the response regulator of the RcsCDB His-Asp phosphorelay signal transduction system. Genetic data indicated that activation of osmCp1 by both NhaR and RcsB requires the same short sequences upstream of the -35 region of the promoter. Accordingly, DNase I footprint analysis indicated that both activators protect an overlapping region close to the -35 box of the promoter and suggested that the regulatory effect is direct. Despite the overlap of the binding sites, each activator acts independent of the other and is specific for a particular stress. NhaR can stimulate osmCp1 in response to an osmotic signal even in the absence of RcsB. RcsB is responsible for the induction of osmCp1 by alteration of the cell envelope, even in the absence of NhaR. osmCp1 as an example of multiple-stress-responsive promoter is discussed in light of a comparison of the NhaR and RcsB target regions in the Enterobacteriaceae.

EnteriX 2003: Visualization tools for genome alignments of Enterobacteriaceae

We describe EnteriX, a suite of three web-based visualization tools for graphically portraying alignment information from comparisons among several fixed and user-supplied sequences from related enterobacterial species, anchored on a reference genome (http://bio.cse.psu.edu/). The first visualization, Enteric, displays stacked pairwise alignments between a reference genome and each of the related bacteria, represented schematically as PIPs (Percent Identity Plots). Encoded in the views are large-scale genomic rearrangement events and functional landmarks. The second visualization, Menteric, computes and displays 1 Kb views of nucleotide-level multiple alignments of the sequences, together with annotations of genes, regulatory sites and conserved regions. The third, a Java-based tool named Maj, displays alignment information in two formats, corresponding roughly to the Enteric and Menteric views, and adds zoom-in capabilities. The uses of such tools are diverse, from examining the multiple sequence alignment to infer conserved sites with potential regulatory roles, to scrutinizing the commonalities and differences between the genomes for pathogenicity or phylogenetic studies. The EnteriX suite currently includes >15 enterobacterial genomes, generates views centered on four different anchor genomes and provides support for including user sequences in the alignments.

Bacterial genomics: potential for antimicrobial drug discovery

The sequencing of entire bacterial genomes is becoming increasingly routine, promising to revolutionise approaches to identifying putative antimicrobial drug targets. In silico methods can be used to identify putative gene products by comparing sequences of biochemically characterised enzymes and proteins with data produced by sequencing projects. Comparative genomics between a pathogenic bacterium versus nonpathogen as well as pathogen versus host can identify molecular targets that would be ideal for future investigation. The aim of these comparisons would be to identify genes that code for pathogenicity factors in the bacterium or genes essential for bacterial survival. The latter set of genes includes those that are nonfunctional or redundant in the host as well as genes absent from the host but essential in the pathogen. The products of these genes would be ideal targets for antimicrobial compounds. If compounds could be generated that disrupt the pathogen's ability to thrive but not affect the host, since there is a lack of the targeted protein, they could prove to be powerful therapeutics. An elegant example illustrating the power of comparative genomics involves comparison of the pathways of bacterial and eukaryotic aminoacyl-tRNA synthesis. Comparison of pathogenic bacterial genomes shows that many bacteria lack the genes encoding either one or two specific aminoacyl-tRNA synthetases, enzymes involved in ensuring correct aminoacylation of tRNA for subsequent translation of the genetic code. Bacteria have an alternative pathway by which amide aminoacyl-tRNAs are formed. Comparative genomics has demonstrated that this pathway is uniquely prokaryotic/archaeal and also relatively widely found in pathogenic bacteria, indicating the potential of the catalytic enzymes of the pathway as targets for novel antimicrobial drugs.

Strategies and tools for whole-genome alignments

The availability of the assembled mouse genome makes possible, for the first time, an alignment and comparison of two large vertebrate genomes. We investigated different strategies of alignment for the subsequent analysis of conservation of genomes that are effective for assemblies of different quality. These strategies were applied to the comparison of the working draft of the human genome with the Mouse Genome Sequencing Consortium assembly, as well as other intermediate mouse assemblies. Our methods are fast and the resulting alignments exhibit a high degree of sensitivity, covering more than 90% of known coding exons in the human genome. We obtained such coverage while preserving specificity. With a view towards the end user, we developed a suite of tools and Web sites for automatically aligning and subsequently browsing and working with whole-genome comparisons. We describe the use of these tools to identify conserved non-coding regions between the human and mouse genomes, some of which have not been identified by other methods.

F-MuLV acceleration of myelomonocytic tumorigenesis in SV40 large T antigen transgenic mice is accompanied by retroviral insertion at Fli1 and a novel locus, Fim4

We describe here the development of a murine system for the identification of genes involved in myelomonocytic neoplasms. Transgenic C57BL/6J mice expressing SV40 early region under a myelomonocytic promoter develop histiocytic sarcomas with a latency of 167 days. We used retroviral proviral tagging to accelerate tumorigenesis and to uncover genetic changes that contribute to tumor development. Infection of transgenic mice with Friend murine leukemia virus (F-MuLV) shortened the latency of morbidity to 103 days (P< 0.001); this was associated with clonal proviral integrations in tumor DNA. As expected for F-MuLV, proviral insertions occurred at Fli1 in both transgenic and nontransgenic tumors. Four insertions were found at a novel locus, termed Fim4, on chromosome 6. This region is syntenic to human 7q32, a region that is commonly deleted in human myelodysplastic syndrome and acute myeloid leukemia. A murine BAC containing Fim4 was sequenced and analyzed, and while there was significant human-mouse homology in the area of the insertions, no candidate gene has been identified. Thus we have established a system to identify genes involved in myelomonocytic tumors, and have used it to identify Fim4, a new common site of proviral insertion. Study of this locus may provide insight into genes involved in AML-associated 7q32 deletions in humans.

Genome compartimentation by a hybrid chromosome model (HXM). Application to Saccharomyces cerevisae subtelomeres

The aim of this paper is to present a new approach, called 'Hybrid Chromosome Model' (HXM), which allows both the extraction of regions of similarity between two sequences, and the compartimentation of a set of DNA sequences. The principle of the method consists in compacting a set of sequences (split into fragments of fixed length) into a 'hybrid chromosome', which results from the stacking of the whole sequence fragments. We have illustrated our approach on the 32 subtelomeres of Saccharomyces cerevisae. The compartimentation of these chromosome extremities into common regions of similarity has been carried out. The approach HXM is a fast and efficient tool for mapping entire genomes and for extracting ancient duplications within or between genomes.

Transcriptional regulation of pentose utilisation systems in the Bacillus/Clostridium group of bacteria

In Bacillus subtilis, utilisation of xylose, arabinose and ribose is controlled by the transcriptional factors XylR, AraR and RbsR, respectively. Here we apply the comparative approach to the analysis of these regulons in the Bacillus/Clostridium group. Evolutionary variability of operon structures is demonstrated and operator sites for the main transcription factors are predicted. The consensus sequences for the XylR and RbsR binding sites vary in different subgroups of genomes. The functional coupling of gene clusters and the conservation of regulatory sites allow for detection of non-orthologous gene displacement of ribulose kinase in Enterococcus faecium and Clostridium acetobutylicum. Moreover, candidate catabolite responsive elements found upstream of most pentose-utilising genes suggest CcpA-mediated catabolite repression.

Complete genome sequence of Salmonella enterica serovar Typhimurium LT2

Salmonella enterica subspecies I, serovar Typhimurium (S. typhimurium), is a leading cause of human gastroenteritis, and is used as a mouse model of human typhoid fever. The incidence of non-typhoid salmonellosis is increasing worldwide, causing millions of infections and many deaths in the human population each year. Here we sequenced the 4,857-kilobase (kb) chromosome and 94-kb virulence plasmid of S. typhimurium strain LT2. The distribution of close homologues of S. typhimurium LT2 genes in eight related enterobacteria was determined using previously completed genomes of three related bacteria, sample sequencing of both S. enterica serovar Paratyphi A (S. paratyphi A) and Klebsiella pneumoniae, and hybridization of three unsequenced genomes to a microarray of S. typhimurium LT2 genes. Lateral transfer of genes is frequent, with 11% of the S. typhimurium LT2 genes missing from S. enterica serovar Typhi (S. typhi), and 29% missing from Escherichia coli K12. The 352 gene homologues of S. typhimurium LT2 confined to subspecies I of S. enterica-containing most mammalian and bird pathogens-are useful for studies of epidemiology, host specificity and pathogenesis. Most of these homologues were previously unknown, and 50 may be exported to the periplasm or outer membrane, rendering them accessible as therapeutic or vaccine targets.

Comparisons of gene colinearity in genomes using GeneOrder2.0

Comparative genomics is enhanced by data mining the rapidly expanding DNA sequence databases. Because of the immense amount of data, computational tools and methods are needed to augment traditional manual visualizations and manipulations of these data. GeneOrder2.0, a Java-based interactive software programme, organizes genome sequence data into tabular and graphical visualizations of the extent of colinearity of genes between any two chromosome genomes of < or =250 kilobases. Both GenBank and proprietary data can be analyzed with this tool.

Comparative analysis of the gene-dense ACHE/TFR2 region on human chromosome 7q22 with the orthologous region on mouse chromosome 5

Chromosome 7q22 has been the focus of many cytogenetic and molecular studies aimed at delineating regions commonly deleted in myeloid leukemias and myelodysplastic syndromes. We have compared a gene-dense, GC-rich sub-region of 7q22 with the orthologous region on mouse chromosome 5. A physical map of 640 kb of genomic DNA from mouse chromosome 5 was derived from a series of overlapping bacterial artificial chromosomes. A 296 kb segment from the physical map, spanning ACHE: to Tfr2, was compared with 267 kb of human sequence. We identified a conserved linkage of 12 genes including an open reading frame flanked by ACHE: and Asr2, a novel cation-chloride cotransporter interacting protein Cip1, Ephb4, Zan and Perq1. While some of these genes have been previously described, in each case we present new data derived from our comparative sequence analysis. Adjacent unfinished sequence data from the mouse contains an orthologous block of 10 additional genes including three novel cDNA sequences that we subsequently mapped to human 7q22. Methods for displaying comparative genomic information, including unfinished sequence data, are becoming increasingly important. We supplement our printed comparative analysis with a new, Web-based program called Laj (local alignments with java). Laj provides interactive access to archived pairwise sequence alignments via the WWW. It displays synchronized views of a dot-plot, a percent identity plot, a nucleotide-level local alignment and a variety of relevant annotations. Our mouse-human comparison can be viewed at http://web.uvic.ca/~bioweb/laj.html. Laj is available at http://bio.cse.psu.edu/, along with online documentation and additional examples of annotated genomic regions.

Comparison of the Escherichia coli K-12 genome with sampled genomes of a Klebsiella pneumoniae and three salmonella enterica serovars, Typhimurium, Typhi and Paratyphi

The Escherichia coli K-12 genome (ECO) was compared with the sampled genomes of the sibling species Salmonella enterica serovars Typhimurium, Typhi and Paratyphi A (collectively referred to as SAL) and the genome of the close outgroup Klebsiella pneumoniae (KPN). There are at least 160 locations where sequences of >400 bp are absent from ECO but present in the genomes of all three SAL and 394 locations where sequences are present in ECO but close homologs are absent in all SAL genomes. The 394 sequences in ECO that do not occur in SAL contain 1350 (30.6%) of the 4405 ECO genes. Of these, 1165 are missing from both SAL and KPN. Most of the 1165 genes are concentrated within 28 regions of 10-40 kb, which consist almost exclusively of such genes. Among these regions were six that included previously identified cryptic phage. A hypothetical ancestral state of genomic regions that differ between ECO and SAL can be inferred in some cases by reference to the genome structure in KPN and the more distant relative Yersinia pestis. However, many changes between ECO and SAL are concentrated in regions where all four genera have a different structure. The rate of gene insertion and deletion is sufficiently high in these regions that the ancestral state of the ECO/SAL lineage cannot be inferred from the present data. The sequencing of other closely related genomes, such as S.bongori or Citrobacter, may help in this regard.