GenomeViz: visualizing microbial genomes

Genome sequence of Lactobacillus helveticus, an organism distinguished by selective gene loss and insertion sequence element expansion

Mobile genetic elements are major contributing factors to the generation of genetic diversity in prokaryotic organisms. For example, insertion sequence (IS) elements have been shown to specifically contribute to niche adaptation by promoting a variety of genetic rearrangements. The complete genome sequence of the cheese culture Lactobacillus helveticus DPC 4571 was determined and revealed significant conservation compared to three nondairy gut lactobacilli. Despite originating from significantly different environments, 65 to 75% of the genes were conserved between the commensal and dairy lactobacilli, which allowed key niche-specific gene sets to be described. However, the primary distinguishing feature was 213 IS elements in the DPC 4571 genome, 10 times more than for the other lactobacilli. Moreover, genome alignments revealed an unprecedented level of genome stability between these four Lactobacillus species, considering the number of IS elements in the L. helveticus genome. Comparative analysis also indicated that the IS elements were not the primary agents of niche adaptation for the L. helveticus genome. A clear bias toward the loss of genes reported to be important for gut colonization was observed for the cheese culture, but there was no clear evidence of IS-associated gene deletion and decay for the majority of genes lost. Furthermore, an extraordinary level of sequence diversity exists between copies of certain IS elements in the DPC 4571 genome, indicating they may represent an ancient component of the L. helveticus genome. These data suggest a special unobtrusive relationship between the DPC 4571 genome and its mobile DNA complement.

Pathogenomics of Listeria spp

This review provides an overview of recent progress in the exploration of genomic, transcriptomic, and proteomic data in Listeria spp. to understand genome evolution and diversity, as well as physiological aspects of metabolism utilized by the bacteria when growing in diverse and varied environments.

Contribution of exogenous genetic elements to the group A Streptococcus metagenome

Variation in gene content among strains of a bacterial species contributes to biomedically relevant differences in phenotypes such as virulence and antimicrobial resistance. Group A Streptococcus (GAS) causes a diverse array of human infections and sequelae, and exhibits a complex pathogenic behavior. To enhance our understanding of genotype-phenotype relationships in this important pathogen, we determined the complete genome sequences of four GAS strains expressing M protein serotypes (M2, M4, and 2 M12) that commonly cause noninvasive and invasive infections. These sequences were compared with eight previously determined GAS genomes and regions of variably present gene content were assessed. Consistent with the previously determined genomes, each of the new genomes is ∼1.9 Mb in size, with ∼10% of the gene content of each encoded on variably present exogenous genetic elements. Like the other GAS genomes, these four genomes are polylysogenic and prophage encode the majority of the variably present gene content of each. In contrast to most of the previously determined genomes, multiple exogenous integrated conjugative elements (ICEs) with characteristics of conjugative transposons and plasmids are present in these new genomes. Cumulatively, 242 new GAS metagenome genes were identified that were not present in the previously sequenced genomes. Importantly, ICEs accounted for 41% of the new GAS metagenome gene content identified in these four genomes. Two large ICEs, designated 2096-RD.2 (63 kb) and 10750-RD.2 (49 kb), have multiple genes encoding resistance to antimicrobial agents, including tetracycline and erythromycin, respectively. Also resident on these ICEs are three genes encoding inferred extracellular proteins of unknown function, including a predicted cell surface protein that is only present in the genome of the serotype M12 strain cultured from a patient with acute poststreptococcal glomerulonephritis. The data provide new information about the GAS metagenome and will assist studies of pathogenesis, antimicrobial resistance, and population genomics.

Plasmids pMOL28 and pMOL30 of Cupriavidus metallidurans are specialized in the maximal viable response to heavy metals

We fully annotated two large plasmids, pMOL28 (164 open reading frames [ORFs]; 171,459 bp) and pMOL30 (247 ORFs; 233,720 bp), in the genome of Cupriavidus metallidurans CH34. pMOL28 contains a backbone of maintenance and transfer genes resembling those found in plasmid pSym of C. taiwanensis and plasmid pHG1 of C. eutrophus, suggesting that they belong to a new class of plasmids. Genes involved in resistance to the heavy metals Co(II), Cr(VI), Hg(II), and Ni(II) are concentrated in a 34-kb region on pMOL28, and genes involved in resistance to Ag(I), Cd(II), Co(II), Cu(II), Hg(II), Pb(II), and Zn(II) occur in a 132-kb region on pMOL30. We identified three putative genomic islands containing metal resistance operons flanked by mobile genetic elements, one on pMOL28 and two on pMOL30. Transcriptomic analysis using quantitative PCR and microarrays revealed metal-mediated up-regulation of 83 genes on pMOL28 and 143 genes on pMOL30 that coded for all known heavy metal resistance proteins, some new heavy metal resistance proteins (czcJ, mmrQ, and pbrU), membrane proteins, truncated transposases, conjugative transfer proteins, and many unknown proteins. Five genes on each plasmid were down-regulated; for one of them, chrI localized on pMOL28, the down-regulation occurred in the presence of five cations. We observed multiple cross-responses (induction of specific metal resistance by other metals), suggesting that the cellular defense of C. metallidurans against heavy metal stress involves various regulons and probably has multiple stages, including a more general response and a more metal-specific response.

Genomorama: genome visualization and analysis

BACKGROUND

The ability to visualize genomic features and design experimental assays that can target specific regions of a genome is essential for modern biology. To assist in these tasks, we present Genomorama, a software program for interactively displaying multiple genomes and identifying potential DNA hybridization sites for assay design.

RESULTS

Useful features of Genomorama include genome search by DNA hybridization (probe binding and PCR amplification), efficient multi-scale display and manipulation of multiple genomes, support for many genome file types and the ability to search for and retrieve data from the National Center for Biotechnology Information (NCBI) Entrez server.

CONCLUSION

Genomorama provides an efficient computational platform for visualizing and analyzing multiple genomes.

Comparative genomic analysis for the presence of potential enterococcal virulence factors in the probiotic Enterococcus faecalis strain Symbioflor 1

Enterococci are members of the natural microbiota of animal and human intestinal tracts and are capable of causing opportunistic infections. They are also used as starter cultures in the food industry as well as in health supplements and probiotics by the pharmaceutical industry. This Janus-faced status requires a careful evaluation on the basis of pathogenic traits to ensure the safety of the strain used to produce food and pharmaceuticals. We performed gapped-genome sequencing of a probiotic strain Enterococcus faecalis Symbioflor 1 and present initial results deriving from comparative genome analysis with that of the previously sequenced pathogenic clinical isolate E. faecalis V583. There was strong overall conservation of synteny between both strains and a detailed analysis revealed the absence of large genomic regions from the chromosome of the probiotic strain, indicating gene loss. Genes absent from the Symbioflor 1 strain included those encoding the enterococcal cytolysin, enterococcal surface protein, and gelatinase (coccolysin) as well as hyaluronidase and the peptide antibiotic AS-48. This data was confirmed using PCR primers specific for the respective genes. However, other enterococcal determinants such as aggregation substance, collagen adhesion protein, the ability to resist oxygen anions as well as capsule formation were detected. The presence of these traits may be advantageous for the strain Symbioflor 1 since they potentially enable colonization and proliferation of the bacterium on mucosal surfaces thereby conferring on it probiotic traits.

The genome sequence of avian pathogenic Escherichia coli strain O1:K1:H7 shares strong similarities with human extraintestinal pathogenic E. coli genomes

Escherichia coli strains that cause disease outside the intestine are known as extraintestinal pathogenic E. coli (ExPEC) and include human uropathogenic E. coli (UPEC) and avian pathogenic E. coli (APEC). Regardless of host of origin, ExPEC strains share many traits. It has been suggested that these commonalities may enable APEC to cause disease in humans. Here, we begin to test the hypothesis that certain APEC strains possess potential to cause human urinary tract infection through virulence genotyping of 1,000 APEC and UPEC strains, generation of the first complete genomic sequence of an APEC (APEC O1:K1:H7) strain, and comparison of this genome to all available human ExPEC genomic sequences. The genomes of APEC O1 and three human UPEC strains were found to be remarkably similar, with only 4.5% of APEC O1's genome not found in other sequenced ExPEC genomes. Also, use of multilocus sequence typing showed that some of the sequenced human ExPEC strains were more like APEC O1 than other human ExPEC strains. This work provides evidence that at least some human and avian ExPEC strains are highly similar to one another, and it supports the possibility that a food-borne link between some APEC and UPEC strains exists. Future studies are necessary to assess the ability of APEC to overcome the hurdles necessary for such a food-borne transmission, and epidemiological studies are required to confirm that such a phenomenon actually occurs.

Burkholderia xenovorans LB400 harbors a multi-replicon, 9.73-Mbp genome shaped for versatility

Burkholderia xenovorans LB400 (LB400), a well studied, effective polychlorinated biphenyl-degrader, has one of the two largest known bacterial genomes and is the first nonpathogenic Burkholderia isolate sequenced. From an evolutionary perspective, we find significant differences in functional specialization between the three replicons of LB400, as well as a more relaxed selective pressure for genes located on the two smaller vs. the largest replicon. High genomic plasticity, diversity, and specialization within the Burkholderia genus are exemplified by the conservation of only 44% of the genes between LB400 and Burkholderia cepacia complex strain 383. Even among four B. xenovorans strains, genome size varies from 7.4 to 9.73 Mbp. The latter is largely explained by our findings that >20% of the LB400 sequence was recently acquired by means of lateral gene transfer. Although a range of genetic factors associated with in vivo survival and intercellular interactions are present, these genetic factors are likely related to niche breadth rather than determinants of pathogenicity. The presence of at least eleven "central aromatic" and twenty "peripheral aromatic" pathways in LB400, among the highest in any sequenced bacterial genome, supports this hypothesis. Finally, in addition to the experimentally observed redundancy in benzoate degradation and formaldehyde oxidation pathways, the fact that 17.6% of proteins have a better LB400 paralog than an ortholog in a different genome highlights the importance of gene duplication and repeated acquirement, which, coupled with their divergence, raises questions regarding the role of paralogs and potential functional redundancies in large-genome microbes.

Toward a more robust assessment of intraspecies diversity, using fewer genetic markers

Phylogenetic sequence analysis of single or multiple genes has dominated the study and census of the genetic diversity among closely related bacteria. It remains unclear, however, how the results based on a few genes in the genome correlate with whole-genome-based relatedness and what genes (if any) best reflect whole-genome-level relatedness and hence should be preferentially used to economize on cost and to improve accuracy. We show here that phylogenies of closely related organisms based on the average nucleotide identity (ANI) of their shared genes correspond accurately to phylogenies based on state-of-the-art analysis of their whole-genome sequences. We use ANI to evaluate the phylogenetic robustness of every gene in the genome and show that almost all core genes, regardless of their functions and positions in the genome, offer robust phylogenetic reconstruction among strains that show 80 to 95% ANI (16S rRNA identity, >98.5%). Lack of elapsed time and, to a lesser extent, horizontal transfer and recombination make the selection of genes more critical for applications that target the intraspecies level, i.e., strains that show >95% ANI according to current standards. A much more accurate phylogeny for the Escherichia coli group was obtained based on just three best-performing genes according to our analysis compared to the concatenated alignment of eight genes that are commonly employed for phylogenetic purposes in this group. Our results are reproducible within the Salmonella, Burkholderia, and Shewanella groups and therefore are expected to have general applicability for microevolution studies, including metagenomic surveys.

Genome sequence of the bioplastic-producing “Knallgas” bacterium Ralstonia eutropha H16

The H(2)-oxidizing lithoautotrophic bacterium Ralstonia eutropha H16 is a metabolically versatile organism capable of subsisting, in the absence of organic growth substrates, on H(2) and CO(2) as its sole sources of energy and carbon. R. eutropha H16 first attracted biotechnological interest nearly 50 years ago with the realization that the organism's ability to produce and store large amounts of poly[R-(-)-3-hydroxybutyrate] and other polyesters could be harnessed to make biodegradable plastics. Here we report the complete genome sequence of the two chromosomes of R. eutropha H16. Together, chromosome 1 (4,052,032 base pairs (bp)) and chromosome 2 (2,912,490 bp) encode 6,116 putative genes. Analysis of the genome sequence offers the genetic basis for exploiting the biotechnological potential of this organism and provides insights into its remarkable metabolic versatility.

Whole-genome sequence of Listeria welshimeri reveals common steps in genome reduction with Listeria innocua as compared to Listeria monocytogenes

We present the complete genome sequence of Listeria welshimeri, a nonpathogenic member of the genus Listeria. Listeria welshimeri harbors a circular chromosome of 2,814,130 bp with 2,780 open reading frames. Comparative genomic analysis of chromosomal regions between L. welshimeri, Listeria innocua, and Listeria monocytogenes shows strong overall conservation of synteny, with the exception of the translocation of an F(o)F(1) ATP synthase. The smaller size of the L. welshimeri genome is the result of deletions in all of the genes involved in virulence and of "fitness" genes required for intracellular survival, transcription factors, and LPXTG- and LRR-containing proteins as well as 55 genes involved in carbohydrate transport and metabolism. In total, 482 genes are absent from L. welshimeri relative to L. monocytogenes. Of these, 249 deletions are commonly absent in both L. welshimeri and L. innocua, suggesting similar genome evolutionary paths from an ancestor. We also identified 311 genes specific to L. welshimeri that are absent in the other two species, indicating gene expansion in L. welshimeri, including horizontal gene transfer. The species L. welshimeri appears to have been derived from early evolutionary events and an ancestor more compact than L. monocytogenes that led to the emergence of nonpathogenic Listeria spp.

Gain and loss of multiple genes during the evolution of Helicobacter pylori

Sequence diversity and gene content distinguish most isolates of Helicobacter pylori. Even greater sequence differences differentiate distinct populations of H. pylori from different continents, but it was not clear whether these populations also differ in gene content. To address this question, we tested 56 globally representative strains of H. pylori and four strains of Helicobacter acinonychis with whole genome microarrays. Of the weighted average of 1,531 genes present in the two sequenced genomes, 25% are absent in at least one strain of H. pylori and 21% were absent or variable in H. acinonychis. We extrapolate that the core genome present in all isolates of H. pylori contains 1,111 genes. Variable genes tend to be small and possess unusual GC content; many of them have probably been imported by horizontal gene transfer. Phylogenetic trees based on the microarray data differ from those based on sequences of seven genes from the core genome. These discrepancies are due to homoplasies resulting from independent gene loss by deletion or recombination in multiple strains, which distort phylogenetic patterns. The patterns of these discrepancies versus population structure allow a reconstruction of the timing of the acquisition of variable genes within this species. Variable genes that are located within the cag pathogenicity island were apparently first acquired en bloc after speciation. In contrast, most other variable genes are of unknown function or encode restriction/modification enzymes, transposases, or outer membrane proteins. These seem to have been acquired prior to speciation of H. pylori and were subsequently lost by convergent evolution within individual strains. Thus, the use of microarrays can reveal patterns of gene gain or loss when examined within a phylogenetic context that is based on sequences of core genes.

Amino Acid Cost and Codon-Usage Biases in 6 Prokaryotic Genomes: A Whole-Genome Analysis

For most prokaryotic organisms, amino acid biosynthesis represents a significant portion of their overall energy budget. The difference in the cost of synthesis between amino acids can be striking, differing by as much as 7-fold. Two prokaryotic organisms, Escherichia coli and Bacillus subtilis, have been shown to preferentially utilize less costly amino acids in highly expressed genes, indicating that parsimony in amino acid selection may confer a selective advantage for prokaryotes. This study confirms those findings and extends them to 4 additional prokaryotic organisms: Chlamydia trachomatis, Chlamydophila pneumoniae AR39, Synechocystis sp. PCC 6803, and Thermus thermophilus HB27. Adherence to codon-usage biases for each of these 6 organisms is inversely correlated with a coding region's average amino acid biosynthetic cost in a fashion that is independent of chemoheterotrophic, photoautotrophic, or thermophilic lifestyle. The obligate parasites C. trachomatis and C. pneumoniae AR39 are incapable of synthesizing many of the 20 common amino acids. Removing auxotrophic amino acids from consideration in these organisms does not alter the overall trend of preferential use of energetically inexpensive amino acids in highly expressed genes.

GRAST: a new way of genome reduction analysis using comparative genomics

MOTIVATION

Establishment of intra-cellular life involved a profound re-configuration of the genetic characteristics of bacteria, including genome reduction and rearrangements. Understanding the mechanisms underlying these phenomena will shed light on the genome rearrangements essential for the development of an intra-cellular lifestyle. Comparison of genomes with differences in their sizes poses statistical as well as computational problems. Little efforts have been made to develop flexible computational tools with which to analyse genome reduction and rearrangements.

RESULTS

Investigation of genome reduction and rearrangements in endosymbionts using a novel computational tool (GRAST) identified gathering of genes with similar functions. Conserved clusters of functionally related genes (CGSCs) were detected. Heterogeneous gene and gene cluster non-functionalization/loss are identified between genome regions, functional gene categories and during evolution. Results show that gene non-functionalisation has accelerated during the last 50 MY of Buchnera's evolution while CGSCs have been static.

GenomeDiagram: a python package for the visualization of large-scale genomic data

We present GenomeDiagram, a flexible, open-source Python module for the visualization of large-scale genomic, comparative genomic and other data with reference to a single chromosome or other biological sequence. GenomeDiagram may be used to generate publication-quality vector graphics, rastered images and in-line streamed graphics for webpages. The package integrates with datatypes from the BioPython project, and is available for Windows, Linux and Mac OS X systems.

AVAILABILITY

GenomeDiagram is freely available as source code (under GNU Public License) at http://bioinf.scri.ac.uk/lp/programs.html, and requires Python 2.3 or higher, and recent versions of the ReportLab and BioPython packages.

SUPPLEMENTARY INFORMATION

A user manual, example code and images are available at http://bioinf.scri.ac.uk/lp/programs.html.

Identification of genomic features using microsyntenies of domains: domain teams

The detection, across several genomes, of local conservation of gene content and proximity considerably helps the prediction of features of interest, such as gene fusions or physical and functional interactions. Here, we want to process realistic models of chromosomes, in which genes (or genomic segments of several genes) can be duplicated within a chromosome, or be absent from some other chromosome(s). Our approach adopts the technique of temporarily forgetting genes and working directly with protein "domains" such as those found in Pfam. This allows the detection of strings of domains that are conserved in their content, but not necessarily in their order, which we refer to as domain teams. The prominent feature of the method is that it relaxes the rigidity of the orthology criterion and avoids many of the pitfalls of gene-families identification methods, often hampered by multidomain proteins or low levels of sequence similarity. This approach, that allows both inter- and intrachromosomal comparisons, proves to be more sensitive than the classical methods based on pairwise sequence comparisons, particularly in the simultaneous treatment of many species. The automated and fast detection of domain teams, together with its increased sensitivity at identifying segments of identical (protein-coding) gene contents as well as gene fusions, should prove a useful complement to other existing methods.

ProGenExpress: visualization of quantitative data on prokaryotic genomes

Background

The integration of genomic information with quantitative experimental data is a key component of systems biology. An increasing number of microbial genomes are being sequenced, leading to an increasing amount of data from post-genomics technologies. The genomes of prokaryotes contain many structures of interest, such as operons, pathogenicity islands and prophage sequences, whose behaviour is of interest during infection and disease. There is a need for simple and novel tools to display and analyse data from these integrated datasets, and we have developed ProGenExpress as a tool for visualising arbitrarily complex numerical data in the context of prokaryotic genomes.

Results

Here we describe ProGenExpress, an R package that allows researchers to easily and quickly visualize quantitative measurements, such as those produced by microarray experiments, in the context of the genome organization of sequenced prokaryotes. Data from microarrays, proteomics or other whole-genome technologies can be accurately displayed on the genome. ProGenExpress can also search for novel regions of interest that consist of groups of adjacent genes that show similar patterns across the experimental data set. We demonstrate ProGenExpress with microarray data from a time-course experiment involving Salmonella typhimurium.

Conclusion

ProGenExpress can be used to visualize quantitative data from complex experiments in the context of the genome of sequenced prokaryotes, and to find novel regions of interest.