Characteristics and genomic epidemiology of colistin-resistant Enterobacterales from farmers, swine, and hospitalized patients in Thailand, 2014-2017

BACKGROUND

Colistin is one of the last resort therapeutic options for treating carbapenemase-producing Enterobacterales, which are resistant to a broad range of beta-lactam antibiotics. However, the increased use of colistin in clinical and livestock farming settings in Thailand and China, has led to the inevitable emergence of colistin resistance. To better understand the rise of colistin-resistant strains in each of these settings, we characterized colistin-resistant Enterobacterales isolated from farmers, swine, and hospitalized patients in Thailand.

METHODS

Enterobacterales were isolated from 149 stool samples or rectal swabs collected from farmers, pigs, and hospitalized patients in Thailand between November 2014-December 2017. Confirmed colistin-resistant isolates were sequenced. Genomic analyses included species identification, multilocus sequence typing, and detection of antimicrobial resistance determinants and plasmids.

RESULTS

The overall colistin-resistant Enterobacterales colonization rate was 26.2% (n = 39/149). The plasmid-mediated colistin-resistance gene (mcr) was detected in all 25 Escherichia coli isolates and 9 of 14 (64.3%) Klebsiella spp. isolates. Five novel mcr allelic variants were also identified: mcr-2.3, mcr-3.21, mcr-3.22, mcr-3.23, and mcr-3.24, that were only detected in E. coli and Klebsiella spp. isolates from farmed pigs.

CONCLUSION

Our data confirmed the presence of colistin-resistance genes in combination with extended spectrum beta-lactamase genes in bacterial isolates from farmers, swine, and patients in Thailand. Differences between the colistin-resistance mechanisms of Escherichia coli and Klebsiella pneumoniae in hospitalized patients were observed, as expected. Additionally, we identified mobile colistin-resistance mcr-1.1 genes from swine and patient isolates belonging to plasmids of the same incompatibility group. This supported the possibility that horizontal transmission of bacterial strains or plasmid-mediated colistin-resistance genes occurs between humans and swine.

DNAmix 2021: Laboratory policies, procedures, and casework scenarios summary and dataset

DNAmix 2021 is a large-scale study conducted to evaluate the extent of consistency and variation among forensic laboratories in the interpretation of DNA mixtures, and to assess the effects of various potential sources of variability. This study utilized a multi-phasic approach designed to collect information about participating laboratories, laboratory policies, and their standard operating procedures (SOPs). It also characterizes the degree of variation in assessments of suitability and number of contributors as well as in comparisons and statistical analyses of DNA mixture profiles. This paper specifically details the study design and the data collected in the first two phases of the study: the Policies & Procedures (P&P) Questionnaire and the Casework Scenarios Questionnaire (CSQ). We report on the variation in policies and SOPs for 86 forensic laboratories-including information about their DNA workflows, systems, and type of statistics reported. We also provide details regarding various case-scenario specific decisions and the nature of mixture casework for 83 forensic laboratories. The data discussed in this article provide insight into the state of the field for forensic DNA mixture interpretation policies and SOPs at the time of the study (2021-2022).

AbGRI4, a novel antibiotic resistance island in multiply antibiotic-resistant Acinetobacter baumannii clinical isolates

OBJECTIVES

To investigate the genomic context of a novel resistance island (RI) in multiply antibiotic-resistant Acinetobacter baumannii clinical isolates and global isolates.

METHODS

Using a combination of long and short reads generated from the Oxford Nanopore and Illumina platforms, contiguous chromosomes and plasmid sequences were determined. BLAST-based analysis was used to identify the RI insertion target.

RESULTS

Genomes of four multiply antibiotic-resistant A. baumannii clinical strains, from a US hospital system, belonging to prevalent MLST ST2 (Pasteur scheme) and ST281 (Oxford scheme) clade F isolates were sequenced to completion. A class 1 integron carrying aadB (tobramycin resistance) and aadA2 (streptomycin/spectinomycin resistance) was identified. The class 1 integron was 6.8 kb, bounded by IS26 at both ends, and embedded in a new target location between an α/β-hydrolase and a reductase. Due to its novel insertion site and unique RI composition, we suggest naming this novel RI AbGRI4. Molecular analysis of global A. baumannii isolates identified multiple AbGRI4 RI variants in non-ST2 clonal lineages, including variations in the resistance gene cassettes, integron backbone and insertion breakpoints at the hydrolase gene.

CONCLUSIONS

A novel RI insertion target harbouring a class 1 integron was identified in a subgroup of ST2/ST281 clinical isolates. Variants of the RI suggested evolution and horizontal transfer of the RI across clonal lineages. Long- and short-read hybrid assembly technology completely resolved the genomic context of IS-bounded RIs, which was not possible using short reads alone.

Genomic and clinical characterisation of multidrug-resistant carbapenemase-producing ST231 and ST16 Klebsiella pneumoniae isolates colonising patients at Siriraj hospital, Bangkok, Thailand from 2015 to 2017

BACKGROUND

Infections caused by carbapenemase-producing Enterobacteriaceae (CPE) have continually grown as a global public health threat, with significant mortality rates observed across the world. We examined the clinical data from patients with CPE infections and their outcomes, concentrating on Klebsiella pneumoniae isolates. We analysed the clinical information, performed antimicrobial susceptibility testing, and conducted molecular epidemiological and genomic analyses on the isolates to identify patterns in the data.

METHODS

The clinical characteristics of 33 hospitalised patients with confirmed CPE, including patient-related factors associated with the development of CPE infections, were examined. Patients were divided according to whether they were "colonised" or "infected" with CPE and by the timing and frequency of their rectal swab collections, from which 45 swabs were randomly selected for analysis. CPE isolates were purified, and antimicrobial susceptibility tests performed. Whole genome sequences of these isolates were determined and analysed to compute bacterial multilocus sequence types and plasmid replicon types, infer phylogenetic relationships, and identify antimicrobial resistance and virulence genes.

RESULTS

Altogether, 88.9% (40/45) of the CPE isolates were K. pneumoniae. The most abundant carbapenemase gene family in the K. pneumoniae isolates (33/39) was bla_OXA-232, with bla_NDM-1 additionally identified in 19 of them. All CPE isolates carrying either bla_OXA-232 or bla_NDM-1 were resistant to meropenem, but only 40 from 45 were susceptible to colistin. Among the CPE-infected patients (n = 18) and CPE-colonised patients who developed CPE infections during the study (n = 3), all but one received standard colistin-based combination therapy. Phylogenetic analysis revealed the polyclonal spread of carbapenemase-producing K. pneumoniae (CPKP) within the patient population, with the following two major subclades identified: ST16 (n = 15) and ST231 (n = 14). CPKP-ST231 had the highest virulence score of 4 and was associated with primary bacteraemia. The siderophores yersiniabactin and aerobactin, considered to be important virulence factors, were only identified in the CPKP-ST231 genomes.

CONCLUSIONS

This study has revealed the genomic features of colonising CPE isolates, focusing on antimicrobial resistance and virulence determinants. This type of multi-layered analysis can be further exploited in Thailand and elsewhere to modify the regimes used for empirical antibiotic treatment and improve the management strategies for CPE infections in hospitalised patients.

Large-scale comparative analysis of microbial pan-genomes using PanOCT

Summary

The JCVI pan-genome pipeline is a collection of programs to run PanOCT and tools that support and extend the capabilities of PanOCT. PanOCT (pan-genome ortholog clustering tool) is a tool for pan-genome analysis of closely related prokaryotic species or strains. The JCVI Pan-Genome Pipeline wrapper invokes command-line utilities that prepare input genomes, invoke third-party tools such as NCBI Blast+, run PanOCT, generate a consensus pan-genome, annotate features of the pan-genome, detect sets of genes of interest such as antimicrobial resistance (AMR) genes and generate figures, tables and html pages to visualize the results. The pipeline can run in a hierarchical mode, lowering the RAM and compute resources used.

Availability and implementation

Source code, demo data, and detailed documentation are freely available at https://github.com/JCVenterInstitute/PanGenomePipeline.

GGRaSP: a R-package for selecting representative genomes using Gaussian mixture models

Motivation

The vast number of available sequenced bacterial genomes occasionally exceeds the facilities of comparative genomic methods or is dominated by a single outbreak strain, and thus a diverse and representative subset is required. Generation of the reduced subset currently requires a priori supervised clustering and sequence-only selection of medoid genomic sequences, independent of any additional genome metrics or strain attributes.

Results

The Gaussian Genome Representative Selector with Prioritization (GGRaSP) R-package described below generates a reduced subset of genomes that prioritizes maintaining genomes of interest to the user as well as minimizing the loss of genetic variation. The package also allows for unsupervised clustering by modeling the genomic relationships using a Gaussian mixture model to select an appropriate cluster threshold. We demonstrate the capabilities of GGRaSP by generating a reduced list of 315 genomes from a genomic dataset of 4600 Escherichia coli genomes, prioritizing selection by type strain and by genome completeness.

Availability and implementaion

GGRaSP is available at https://github.com/JCVenterInstitute/ggrasp/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Emergence of New Delhi Metallo-β-Lactamase (NDM-5) in Klebsiella quasipneumoniae from Neonates in a Nigerian Hospital

Outbreaks of infection occur more often than they are reported in most developing countries, largely due to poor diagnostic services. A Klebsiella species bacteremia outbreak in a newborn unit with high mortality was recently encountered at a location being surveilled for childhood bacteremia. These surveillance efforts offered the opportunity to determine the cause of this neonatal outbreak. In this report, we present the whole-genome sequences of New Delhi metallo-β-lactamase (NDM-5)-containing Klebsiella quasipneumoniae subsp. similipneumoniae bloodstream isolates from a neonatal bacteremia outbreak at a tertiary hospital in Nigeria and as part of the largest collection of K. pneumoniae bloodstream isolates from children in Africa. Comparative analysis of the genetic environment surrounding the NDM-5 genes revealed nearly perfect sequence identity to blaNDM-5-bearing IncX3-type plasmids from other members of the EnterobacteriaceaeIMPORTANCE Carbapenem-resistant Klebsiella pneumoniae is of global health importance, yet there is a paucity of genome-based studies in Africa. Here we report fatal blood-borne NDM-5-producing K. quasipneumoniae subsp. similipneumoniae infections from Nigeria, Africa. New Delhi metallo-β-lactamase (NDM)-producing Klebsiella spp. are responsible for high mortality and morbidity, with the NDM-5 variant showing elevated carbapenem resistance. The prevalence of NDM-5 in Klebsiella has been limited primarily to K. pneumoniae, with only one isolate being collected from Africa. During an outbreak of sepsis in a teaching hospital in Nigeria, five NDM-5-producing K. quasipneumoniae subsp. similipneumoniae sequence type 476 isolates were identified. Given the increased resistance profile of these strains, this study highlights the emerging threat of blaNDM-5 dissemination in hospital environments. The observation of these NDM-5-producing isolates in Africa stresses the urgency to improve monitoring and clinical practices to reduce or prevent the further spread of resistance.

PanACEA: a bioinformatics tool for the exploration and visualization of bacterial pan-chromosomes

BACKGROUND

Bacterial pan-genomes, comprised of conserved and variable genes across multiple sequenced bacterial genomes, allow for identification of genomic regions that are phylogenetically discriminating or functionally important. Pan-genomes consist of large amounts of data, which can restrict researchers ability to locate and analyze these regions. Multiple software packages are available to visualize pan-genomes, but currently their ability to address these concerns are limited by using only pre-computed data sets, prioritizing core over variable gene clusters, or by not accounting for pan-chromosome positioning in the viewer.

RESULTS

We introduce PanACEA (Pan-genome Atlas with Chromosome Explorer and Analyzer), which utilizes locally-computed interactive web-pages to view ordered pan-genome data. It consists of multi-tiered, hierarchical display pages that extend from pan-chromosomes to both core and variable regions to single genes. Regions and genes are functionally annotated to allow for rapid searching and visual identification of regions of interest with the option that user-supplied genomic phylogenies and metadata can be incorporated. PanACEA's memory and time requirements are within the capacities of standard laptops. The capability of PanACEA as a research tool is demonstrated by highlighting a variable region important in differentiating strains of Enterobacter hormaechei.

CONCLUSIONS

PanACEA can rapidly translate the results of pan-chromosome programs into an intuitive and interactive visual representation. It will empower researchers to visually explore and identify regions of the pan-chromosome that are most biologically interesting, and to obtain publication quality images of these regions.

Multidrug resistant pathogens respond differently to the presence of co-pathogen, commensal, probiotic and host cells

In light of the ongoing antimicrobial resistance crisis, there is a need to understand the role of co-pathogens, commensals, and the local microbiome in modulating virulence and antibiotic resistance. To identify possible interactions that influence the expression of virulence or survival mechanisms in both the multidrug-resistant organisms (MDROs) and human host cells, unique cohorts of clinical isolates were selected for whole genome sequencing with enhanced assembly and full annotation, pairwise co-culturing, and transcriptome profiling. The MDROs were co-cultured in pairwise combinations either with: (1) another MDRO, (2) skin commensals (Staphylococcus epidermidis and Corynebacterium jeikeium), (3) the common probiotic Lactobacillus reuteri, and (4) human fibroblasts. RNA-Seq analysis showed distinct regulation of virulence and antimicrobial resistance gene responses across different combinations of MDROs, commensals, and human cells. Co-culture assays demonstrated that microbial interactions can modulate gene responses of both the target and pathogen/commensal species, and that the responses are specific to the identity of the pathogen/commensal species. In summary, bacteria have mechanisms to distinguish between friends, foe and host cells. These results provide foundational data and insight into the possibility of manipulating the local microbiome when treating complicated polymicrobial wound, intra-abdominal, or respiratory infections.

Enterobacterhormaechei subsp. hoffmannii subsp. nov., Enterobacter hormaechei subsp. xiangfangensis comb. nov., Enterobacter roggenkampii sp. nov., and Enterobacter muelleri is a later heterotypic synonym of Enterobacter asburiae based on computational analysis of sequenced Enterobacter genomes

Background: The predominant species in clinical Enterobacter isolates is E. hormaechei. Many articles, clinicians, and GenBank submissions misname these strains as E. cloacae. The lack of sequenced type strains or named species/subspecies for some clades in the E. cloacae complex complicate the issue. Methods: The genomes of the type strains for Enterobacter hormaechei subsp.  oharae, E.  hormaechei subsp.  steigerwaltii, and E. xiangfangensis, and two strains from Hoffmann clusters III and IV of the E. cloacae complex were sequenced. These genomes, the E.  hormaechei subsp.  hormaechei type strain, and other available Enterobacter type strains were analysed in conjunction with all extant Enterobacter genomes in NCBI's RefSeq using Average Nucleotide Identity (ANI). Results: There were five recognizable subspecies of E. hormaechei: E. hormaechei subsp. hoffmannii subsp. nov., E. hormaechei subsp. xiangfangensis comb. nov., and the three previously known subspecies. One of the strains sequenced from the E. cloacae complex was not a novel E. hormaechei subspecies but rather a member of a clade of a novel species: E. roggenkampii sp. nov.. E. muelleri was determined to be a later heterotypic synonym of E. asburiae which should take precedence. Conclusion: The phylogeny of the Enterobacter genus, particularly the cloacae complex, was re-evaluated based on the type strain genome sequences and all other available Enterobacter genomes in RefSeq.

Enterobacterhormaechei subsp. hoffmannii subsp. nov., Enterobacter hormaechei subsp. xiangfangensis comb. nov., Enterobacter roggenkampii sp. nov., and Enterobacter muelleri is a later heterotypic synonym of Enterobacter asburiae based on computational analysis of sequenced Enterobacter genomes

Background: The predominant species in clinical Enterobacter isolates is E. hormaechei. Many articles, clinicians, and GenBank submissions misname these strains as E. cloacae. The lack of sequenced type strains or named species/subspecies for some clades in the E. cloacae complex complicate the issue. Methods: The genomes of the type strains for Enterobacter hormaechei subsp.  oharae, E.  hormaechei subsp.  steigerwaltii, and E. xiangfangensis, and two strains from Hoffmann clusters III and IV of the E. cloacae complex were sequenced. These genomes, the E.  hormaechei subsp.  hormaechei type strain, and other available Enterobacter type strains were analysed in conjunction with all extant Enterobacter genomes in NCBI's RefSeq using Average Nucleotide Identity (ANI). Results: There were five recognizable subspecies of E. hormaechei: E. hormaechei subsp. hoffmannii subsp. nov., E. hormaechei subsp. xiangfangensis comb. nov., and the three previously known subspecies. One of the strains sequenced from the E. cloacae complex was not a novel E. hormaechei subspecies but rather a member of a clade of a novel species: E. roggenkampii sp. nov.. E. muelleri was determined to be a later heterotypic synonym of E. asburiae which should take precedence. Conclusion: The phylogeny of the Enterobacter genus, particularly the cloacae complex, was re-evaluated based on the type strain genome sequences and all other available Enterobacter genomes in RefSeq.

LOCUST: a custom sequence locus typer for classifying microbial isolates

Summary

LOCUST is a custom sequence locus typer tool for classifying microbial genomes. It provides a fully automated opportunity to customize the classification of genome-wide nucleotide variant data most relevant to biological research.

Availability and Implementation

Source code, demo data, and detailed documentation are freely available at http://sourceforge.net/projects/locustyper .

Contact

lbrinkac@jcvi.org.

Supplementary information

Supplementary data are available at Bioinformatics online.

Identification and taxonomic characterization of Bordetella pseudohinzii sp. nov. isolated from laboratory-raised mice

Bordetella hinzii is known to cause respiratory disease in poultry and has been associated with a variety of infections in immunocompromised humans. In addition, there are several reports of B. hinzii infections in laboratory-raised mice. Here we sequenced and analysed the complete genome sequences of multiple B. hinzii-like isolates, obtained from vendor-supplied C57BL/6 mice in animal research facilities on different continents, and we determined their taxonomic relationship to other Bordetella species. The whole-genome based and 16S rRNA gene based phylogenies each identified two separate clades in B. hinzii, one was composed of strains isolated from poultry, humans and a rabbit whereas the other clade was restricted to isolates from mice. Distinctly different estimated DNA–DNA hybridization values, average nucleotide identity scores, gene content, metabolic profiles and host specificity all provide compelling evidence for delineation of the two species, B. hinzii – from poultry, humans and rabbit – and Bordetella pseudohinzii sp. nov. type strain 8-296-03^T (=NRRL B-59942^T=NCTC 13808^T) that infect mice.

Draft genome sequences for the obligate bacterial predators Bacteriovorax spp. of four phylogenetic clusters

Bacteriovorax is the halophilic genus of the obligate bacterial predators, Bdellovibrio and like organisms. The predators are known for their unique biphasic life style in which they search for and attack their prey in the free living phase; penetrate, grow, multiply and lyse the prey in the intraperiplasmic phase. Bacteriovorax isolates representing four phylogenetic clusters were selected for genomic sequencing. Only one type strain genome has been published so far from the genus Bacteriovorax. We report the genomes from non-type strains isolated from aquatic environments. Here we describe and compare the genomic features of the four strains, together with the classification and annotation.

Comparative genomics of plant-associated Pseudomonas spp.: insights into diversity and inheritance of traits involved in multitrophic interactions

We provide here a comparative genome analysis of ten strains within the Pseudomonas fluorescens group including seven new genomic sequences. These strains exhibit a diverse spectrum of traits involved in biological control and other multitrophic interactions with plants, microbes, and insects. Multilocus sequence analysis placed the strains in three sub-clades, which was reinforced by high levels of synteny, size of core genomes, and relatedness of orthologous genes between strains within a sub-clade. The heterogeneity of the P. fluorescens group was reflected in the large size of its pan-genome, which makes up approximately 54% of the pan-genome of the genus as a whole, and a core genome representing only 45–52% of the genome of any individual strain. We discovered genes for traits that were not known previously in the strains, including genes for the biosynthesis of the siderophores achromobactin and pseudomonine and the antibiotic 2-hexyl-5-propyl-alkylresorcinol; novel bacteriocins; type II, III, and VI secretion systems; and insect toxins. Certain gene clusters, such as those for two type III secretion systems, are present only in specific sub-clades, suggesting vertical inheritance. Almost all of the genes associated with multitrophic interactions map to genomic regions present in only a subset of the strains or unique to a specific strain. To explore the evolutionary origin of these genes, we mapped their distributions relative to the locations of mobile genetic elements and repetitive extragenic palindromic (REP) elements in each genome. The mobile genetic elements and many strain-specific genes fall into regions devoid of REP elements (i.e., REP deserts) and regions displaying atypical tri-nucleotide composition, possibly indicating relatively recent acquisition of these loci. Collectively, the results of this study highlight the enormous heterogeneity of the P. fluorescens group and the importance of the variable genome in tailoring individual strains to their specific lifestyles and functional repertoire.

We sequenced the genomes of seven strains of the Pseudomonas fluorescens group that colonize plant surfaces and function as biological control agents, protecting plants from disease. In this study, we demonstrated the genomic diversity of the group by comparing these strains to each other and to three other strains that were sequenced previously. Only about half of the genes in each strain are present in all of the other strains, and each strain has hundreds of unique genes that are not present in the other genomes. We mapped the genes that contribute to biological control in each genome and found that most of the biological control genes are in the variable regions of the genome, which are not shared by all of the other strains. This finding is consistent with our knowledge of the distinctive biology of each strain. Finally, we looked for new genes that are likely to confer antimicrobial traits needed to suppress plant pathogens, but have not been identified previously. In each genome, we discovered many of these new genes, which provide avenues for future discovery of new traits with the potential to manage plant diseases in agriculture or natural ecosystems.

New assembly, reannotation and analysis of the Entamoeba histolytica genome reveal new genomic features and protein content information

Background

In order to maintain genome information accurately and relevantly, original genome annotations need to be updated and evaluated regularly. Manual reannotation of genomes is important as it can significantly reduce the propagation of errors and consequently diminishes the time spent on mistaken research. For this reason, after five years from the initial submission of the Entamoeba histolytica draft genome publication, we have re-examined the original 23 Mb assembly and the annotation of the predicted genes.

Principal Findings

The evaluation of the genomic sequence led to the identification of more than one hundred artifactual tandem duplications that were eliminated by re-assembling the genome. The reannotation was done using a combination of manual and automated genome analysis. The new 20 Mb assembly contains 1,496 scaffolds and 8,201 predicted genes, of which 60% are identical to the initial annotation and the remaining 40% underwent structural changes. Functional classification of 60% of the genes was modified based on recent sequence comparisons and new experimental data. We have assigned putative function to 3,788 proteins (46% of the predicted proteome) based on the annotation of predicted gene families, and have identified 58 protein families of five or more members that share no homology with known proteins and thus could be entamoeba specific. Genome analysis also revealed new features such as the presence of segmental duplications of up to 16 kb flanked by inverted repeats, and the tight association of some gene families with transposable elements.

Significance

This new genome annotation and analysis represents a more refined and accurate blueprint of the pathogen genome, and provides an upgraded tool as reference for the study of many important aspects of E. histolytica biology, such as genome evolution and pathogenesis.

Entamoeba histolytica is an anaerobic parasitic protozoan that causes amoebic dysentery. The parasites colonize the large intestine, but under some circumstances may invade the intestinal mucosa, enter the bloodstream and lead to the formation of abscesses such amoebic liver abscesses. The draft genome of E. histolytica, published in 2005, provided the scientific community with the first comprehensive view of the gene set for this parasite and important tools for elucidating the genetic basis of Entamoeba pathogenicity. Because complete genetic knowledge is critical for drug discovery and potential vaccine development for amoebiases, we have re-examined the original draft genome for E. histolytica. We have corrected the sequence assembly, improved the gene predictions and refreshed the functional gene assignments. As a result, this effort has led to a more accurate gene annotation, and the discovery of novel features, such as the presence of genome segmental duplications and the close association of some gene families with transposable elements. We believe that continuing efforts to improve genomic data will undoubtedly help to identify and characterize potential targets for amoebiasis control, as well as to contribute to a better understanding of genome evolution and pathogenesis for this parasite.

Meeting report: a workshop on Best Practices in Genome Annotation

Efforts to annotate the genomes of a wide variety of model organisms are currently carried out by sequencing centers, model organism databases and academic/institutional laboratories around the world. Different annotation methods and tools have been developed over time to meet the needs of biologists faced with the task of annotating biological data. While standardized methods are essential for consistent curation within each annotation group, methods and tools can differ between groups, especially when the groups are curating different organisms. Biocurators from several institutes met at the Third International Biocuration Conference in Berlin, Germany, April 2009 and hosted the ‘Best Practices in Genome Annotation: Inference from Evidence’ workshop to share their strategies, pipelines, standards and tools. This article documents the material presented in the workshop.

The Protein Naming Utility: a rules database for protein nomenclature

Generation of syntactically correct and unambiguous names for proteins is a challenging, yet vital task for functional annotation processes. Proteins are often named based on homology to known proteins, many of which have problematic names. To address the need to generate high-quality protein names, and capture our significant experience correcting protein names manually, we have developed the Protein Naming Utility (PNU, http://www.jcvi.org/pn-utility). The PNU is a web-based database for storing and applying naming rules to identify and correct syntactically incorrect protein names, or to replace synonyms with their preferred name. The PNU allows users to generate and manage collections of naming rules, optionally building upon the growing body of rules generated at the J. Craig Venter Institute (JCVI). Since communities often enforce disparate conventions for naming proteins, the PNU supports grouping rules into user-managed collections. Users can check their protein names against a selected PNU rule collection, generating both statistics and corrected names. The PNU can also be used to correct GenBank table files prior to submission to GenBank. Currently, the database features 3080 manual rules that have been entered by JCVI Bioinformatics Analysts as well as 7458 automatically imported names.

Pathema: a clade-specific bioinformatics resource center for pathogen research

Pathema (http://pathema.jcvi.org) is one of the eight Bioinformatics Resource Centers (BRCs) funded by the National Institute of Allergy and Infectious Disease (NIAID) designed to serve as a core resource for the bio-defense and infectious disease research community. Pathema strives to support basic research and accelerate scientific progress for understanding, detecting, diagnosing and treating an established set of six target NIAID Category A-C pathogens: Category A priority pathogens; Bacillus anthracis and Clostridium botulinum, and Category B priority pathogens; Burkholderia mallei, Burkholderia pseudomallei, Clostridium perfringens and Entamoeba histolytica. Each target pathogen is represented in one of four distinct clade-specific Pathema web resources and underlying databases developed to target the specific data and analysis needs of each scientific community. All publicly available complete genome projects of phylogenetically related organisms are also represented, providing a comprehensive collection of organisms for comparative analyses. Pathema facilitates the scientific exploration of genomic and related data through its integration with web-based analysis tools, customized to obtain, display, and compute results relevant to ongoing pathogen research. Pathema serves the bio-defense and infectious disease research community by disseminating data resulting from pathogen genome sequencing projects and providing access to the results of inter-genomic comparisons for these organisms.

Genome degradation in Brucella ovis corresponds with narrowing of its host range and tissue tropism

Brucella ovis is a veterinary pathogen associated with epididymitis in sheep. Despite its genetic similarity to the zoonotic pathogens B. abortus, B. melitensis and B. suis, B. ovis does not cause zoonotic disease. Genomic analysis of the type strain ATCC25840 revealed a high percentage of pseudogenes and increased numbers of transposable elements compared to the zoonotic Brucella species, suggesting that genome degradation has occurred concomitant with narrowing of the host range of B. ovis. The absence of genomic island 2, encoding functions required for lipopolysaccharide biosynthesis, as well as inactivation of genes encoding urease, nutrient uptake and utilization, and outer membrane proteins may be factors contributing to the avirulence of B. ovis for humans. A 26.5 kb region of B. ovis ATCC25840 Chromosome II was absent from all the sequenced human pathogenic Brucella genomes, but was present in all of 17 B. ovis isolates tested and in three B. ceti isolates, suggesting that this DNA region may be of use for differentiating B. ovis from other Brucella spp. This is the first genomic analysis of a non-zoonotic Brucella species. The results suggest that inactivation of genes involved in nutrient acquisition and utilization, cell envelope structure and urease may have played a role in narrowing of the tissue tropism and host range of B. ovis.

Comparative genomic evidence for a close relationship between the dimorphic prosthecate bacteria Hyphomonas neptunium and Caulobacter crescentus

The dimorphic prosthecate bacteria (DPB) are alpha-proteobacteria that reproduce in an asymmetric manner rather than by binary fission and are of interest as simple models of development. Prior to this work, the only member of this group for which genome sequence was available was the model freshwater organism Caulobacter crescentus. Here we describe the genome sequence of Hyphomonas neptunium, a marine member of the DPB that differs from C. crescentus in that H. neptunium uses its stalk as a reproductive structure. Genome analysis indicates that this organism shares more genes with C. crescentus than it does with Silicibacter pomeroyi (a closer relative according to 16S rRNA phylogeny), that it relies upon a heterotrophic strategy utilizing a wide range of substrates, that its cell cycle is likely to be regulated in a similar manner to that of C. crescentus, and that the outer membrane complements of H. neptunium and C. crescentus are remarkably similar. H. neptunium swarmer cells are highly motile via a single polar flagellum. With the exception of cheY and cheR, genes required for chemotaxis were absent in the H. neptunium genome. Consistent with this observation, H. neptunium swarmer cells did not respond to any chemotactic stimuli that were tested, which suggests that H. neptunium motility is a random dispersal mechanism for swarmer cells rather than a stimulus-controlled navigation system for locating specific environments. In addition to providing insights into bacterial development, the H. neptunium genome will provide an important resource for the study of other interesting biological processes including chromosome segregation, polar growth, and cell aging.

Genome sequence of Synechococcus CC9311: Insights into adaptation to a coastal environment

Coastal aquatic environments are typically more highly productive and dynamic than open ocean ones. Despite these differences, cyanobacteria from the genus Synechococcus are important primary producers in both types of ecosystems. We have found that the genome of a coastal cyanobacterium, Synechococcus sp. strain CC9311, has significant differences from an open ocean strain, Synechococcus sp. strain WH8102, and these are consistent with the differences between their respective environments. CC9311 has a greater capacity to sense and respond to changes in its (coastal) environment. It has a much larger capacity to transport, store, use, or export metals, especially iron and copper. In contrast, phosphate acquisition seems less important, consistent with the higher concentration of phosphate in coastal environments. CC9311 is predicted to have differences in its outer membrane lipopolysaccharide, and this may be characteristic of the speciation of some cyanobacterial groups. In addition, the types of potentially horizontally transferred genes are markedly different between the coastal and open ocean genomes and suggest a more prominent role for phages in horizontal gene transfer in oligotrophic environments.

Skewed genomic variability in strains of the toxigenic bacterial pathogen, Clostridium perfringens

Clostridium perfringens is a Gram-positive, anaerobic spore-forming bacterium commonly found in soil, sediments, and the human gastrointestinal tract. C. perfringens is responsible for a wide spectrum of disease, including food poisoning, gas gangrene (clostridial myonecrosis), enteritis necroticans, and non-foodborne gastrointestinal infections. The complete genome sequences of Clostridium perfringens strain ATCC 13124, a gas gangrene isolate and the species type strain, and the enterotoxin-producing food poisoning strain SM101, were determined and compared with the published C. perfringens strain 13 genome. Comparison of the three genomes revealed considerable genomic diversity with >300 unique "genomic islands" identified, with the majority of these islands unusually clustered on one replichore. PCR-based analysis indicated that the large genomic islands are widely variable across a large collection of C. perfringens strains. These islands encode genes that correlate to differences in virulence and phenotypic characteristics of these strains. Significant differences between the strains include numerous novel mobile elements and genes encoding metabolic capabilities, strain-specific extracellular polysaccharide capsule, sporulation factors, toxins, and other secreted enzymes, providing substantial insight into this medically important bacterial pathogen.

Comparative genomics of emerging human ehrlichiosis agents

Anaplasma (formerly Ehrlichia) phagocytophilum, Ehrlichia chaffeensis, and Neorickettsia (formerly Ehrlichia) sennetsu are intracellular vector-borne pathogens that cause human ehrlichiosis, an emerging infectious disease. We present the complete genome sequences of these organisms along with comparisons to other organisms in the Rickettsiales order. Ehrlichia spp. and Anaplasma spp. display a unique large expansion of immunodominant outer membrane proteins facilitating antigenic variation. All Rickettsiales have a diminished ability to synthesize amino acids compared to their closest free-living relatives. Unlike members of the Rickettsiaceae family, these pathogenic Anaplasmataceae are capable of making all major vitamins, cofactors, and nucleotides, which could confer a beneficial role in the invertebrate vector or the vertebrate host. Further analysis identified proteins potentially involved in vacuole confinement of the Anaplasmataceae, a life cycle involving a hematophagous vector, vertebrate pathogenesis, human pathogenesis, and lack of transovarial transmission. These discoveries provide significant insights into the biology of these obligate intracellular pathogens.

Ehrlichiosis is an acute disease that triggers flu-like symptoms in both humans and animals. It is caused by a range of bacteria transmitted by ticks or flukes. Because these bacteria are difficult to culture, however, the organisms are poorly understood. The genomes of three emerging human pathogens causing ehrlichiosis were sequenced. A database was designed to allow the comparison of these three genomes to sixteen other bacteria with similar lifestyles. Analysis from this database reveals new species-specific and disease-specific genes indicating niche adaptations, pathogenic traits, and other features. In particular, one of the organisms contains more than 100 copies of a single gene involved in interactions with the host(s). These comparisons also enabled a reconstruction of the metabolic potential of five representative genomes from these bacteria and their close relatives. With this work, scientists can study these emerging pathogens in earnest.

Life in hot carbon monoxide: the complete genome sequence of Carboxydothermus hydrogenoformans Z-2901

We report here the sequencing and analysis of the genome of the thermophilic bacterium Carboxydothermus hydrogenoformans Z-2901. This species is a model for studies of hydrogenogens, which are diverse bacteria and archaea that grow anaerobically utilizing carbon monoxide (CO) as their sole carbon source and water as an electron acceptor, producing carbon dioxide and hydrogen as waste products. Organisms that make use of CO do so through carbon monoxide dehydrogenase complexes. Remarkably, analysis of the genome of C. hydrogenoformans reveals the presence of at least five highly differentiated anaerobic carbon monoxide dehydrogenase complexes, which may in part explain how this species is able to grow so much more rapidly on CO than many other species. Analysis of the genome also has provided many general insights into the metabolism of this organism which should make it easier to use it as a source of biologically produced hydrogen gas. One surprising finding is the presence of many genes previously found only in sporulating species in the Firmicutes Phylum. Although this species is also a Firmicutes, it was not known to sporulate previously. Here we show that it does sporulate and because it is missing many of the genes involved in sporulation in other species, this organism may serve as a “minimal” model for sporulation studies. In addition, using phylogenetic profile analysis, we have identified many uncharacterized gene families found in all known sporulating Firmicutes, but not in any non-sporulating bacteria, including a sigma factor not known to be involved in sporulation previously.

Carboxydothermus hydrogenoformans, a bacterium isolated from a Russian hotspring, is studied for three major reasons: it grows at very high temperature, it lives almost entirely on a diet of carbon monoxide (CO), and it converts water to hydrogen gas as part of its metabolism. Understanding this organism's unique biology gets a boost from the decoding of its genome, reported in this issue of PLoS Genetics. For example, genome analysis reveals that it encodes five different forms of the protein machine carbon monoxide dehydrogenase (CODH). Most species have no CODH and even species that utilize CO usually have only one or two. The five CODH in C. hydrogenoformans likely allow it to both use CO for diverse cellular processes and out-compete for it when it is limiting. The genome sequence also led the researchers to experimentally document new aspects of this species' biology including the ability to form spores. The researchers then used comparative genomic analysis to identify conserved genes found in all spore-forming species, including Bacillus anthracis, and not in any other species. Finally, the genome sequence and analysis reported here will aid in those trying to develop this and other species into systems to biologically produce hydrogen gas from water.

Whole-genome sequence analysis of Pseudomonas syringae pv. phaseolicola 1448A reveals divergence among pathovars in genes involved in virulence and transposition

Pseudomonas syringae pv. phaseolicola, a gram-negative bacterial plant pathogen, is the causal agent of halo blight of bean. In this study, we report on the genome sequence of P. syringae pv. phaseolicola isolate 1448A, which encodes 5,353 open reading frames (ORFs) on one circular chromosome (5,928,787 bp) and two plasmids (131,950 bp and 51,711 bp). Comparative analyses with a phylogenetically divergent pathovar, P. syringae pv. tomato DC3000, revealed a strong degree of conservation at the gene and genome levels. In total, 4,133 ORFs were identified as putative orthologs in these two pathovars using a reciprocal best-hit method, with 3,941 ORFs present in conserved, syntenic blocks. Although these two pathovars are highly similar at the physiological level, they have distinct host ranges; 1448A causes disease in beans, and DC3000 is pathogenic on tomato and Arabidopsis. Examination of the complement of ORFs encoding virulence, fitness, and survival factors revealed a substantial, but not complete, overlap between these two pathovars. Another distinguishing feature between the two pathovars is their distinctive sets of transposable elements. With access to a fifth complete pseudomonad genome sequence, we were able to identify 3,567 ORFs that likely comprise the core Pseudomonas genome and 365 ORFs that are P. syringae specific.

Genome analysis of multiple pathogenic isolates of Streptococcus agalactiae: implications for the microbial "pan-genome"

The development of efficient and inexpensive genome sequencing methods has revolutionized the study of human bacterial pathogens and improved vaccine design. Unfortunately, the sequence of a single genome does not reflect how genetic variability drives pathogenesis within a bacterial species and also limits genome-wide screens for vaccine candidates or for antimicrobial targets. We have generated the genomic sequence of six strains representing the five major disease-causing serotypes of Streptococcus agalactiae, the main cause of neonatal infection in humans. Analysis of these genomes and those available in databases showed that the S. agalactiae species can be described by a pan-genome consisting of a core genome shared by all isolates, accounting for approximately 80% of any single genome, plus a dispensable genome consisting of partially shared and strain-specific genes. Mathematical extrapolation of the data suggests that the gene reservoir available for inclusion in the S. agalactiae pan-genome is vast and that unique genes will continue to be identified even after sequencing hundreds of genomes.

The psychrophilic lifestyle as revealed by the genome sequence of Colwellia psychrerythraea 34H through genomic and proteomic analyses

The completion of the 5,373,180-bp genome sequence of the marine psychrophilic bacterium Colwellia psychrerythraea 34H, a model for the study of life in permanently cold environments, reveals capabilities important to carbon and nutrient cycling, bioremediation, production of secondary metabolites, and cold-adapted enzymes. From a genomic perspective, cold adaptation is suggested in several broad categories involving changes to the cell membrane fluidity, uptake and synthesis of compounds conferring cryotolerance, and strategies to overcome temperature-dependent barriers to carbon uptake. Modeling of three-dimensional protein homology from bacteria representing a range of optimal growth temperatures suggests changes to proteome composition that may enhance enzyme effectiveness at low temperatures. Comparative genome analyses suggest that the psychrophilic lifestyle is most likely conferred not by a unique set of genes but by a collection of synergistic changes in overall genome content and amino acid composition.

Complete genome sequence of the plant commensal Pseudomonas fluorescens Pf-5

Pseudomonas fluorescens Pf-5 is a plant commensal bacterium that inhabits the rhizosphere and produces secondary metabolites that suppress soilborne plant pathogens. The complete sequence of the 7.1-Mb Pf-5 genome was determined. We analyzed repeat sequences to identify genomic islands that, together with other approaches, suggested P. fluorescens Pf-5's recent lateral acquisitions include six secondary metabolite gene clusters, seven phage regions and a mobile genomic island. We identified various features that contribute to its commensal lifestyle on plants, including broad catabolic and transport capabilities for utilizing plant-derived compounds, the apparent ability to use a diversity of iron siderophores, detoxification systems to protect from oxidative stress, and the lack of a type III secretion system and toxins found in related pathogens. In addition to six known secondary metabolites produced by P. fluorescens Pf-5, three novel secondary metabolite biosynthesis gene clusters were also identified that may contribute to the biocontrol properties of P. fluorescens Pf-5.

Genome sequence of the PCE-dechlorinating bacterium Dehalococcoides ethenogenes

Dehalococcoides ethenogenes is the only bacterium known to reductively dechlorinate the groundwater pollutants, tetrachloroethene (PCE) and trichloroethene, to ethene. Its 1,469,720-base pair chromosome contains large dynamic duplicated regions and integrated elements. Genes encoding 17 putative reductive dehalogenases, nearly all of which were adjacent to genes for transcription regulators, and five hydrogenase complexes were identified. These findings, plus a limited repertoire of other metabolic modes, indicate that D. ethenogenes is highly evolved to utilize halogenated organic compounds and H2. Diversification of reductive dehalogenase functions appears to have been mediated by recent genetic exchange and amplification. Genome analysis provides insights into the organism's complex nutrient requirements and suggests that an ancestor was a nitrogen-fixing autotroph.

Major structural differences and novel potential virulence mechanisms from the genomes of multiple campylobacter species

Sequencing and comparative genome analysis of four strains of Campylobacter including C. lari RM2100, C. upsaliensis RM3195, and C. coli RM2228 has revealed major structural differences that are associated with the insertion of phage- and plasmid-like genomic islands, as well as major variations in the lipooligosaccharide complex. Poly G tracts are longer, are greater in number, and show greater variability in C. upsaliensis than in the other species. Many genes involved in host colonization, including racR/S, cadF, cdt, ciaB, and flagellin genes, are conserved across the species, but variations that appear to be species specific are evident for a lipooligosaccharide locus, a capsular (extracellular) polysaccharide locus, and a novel Campylobacter putative licABCD virulence locus. The strains also vary in their metabolic profiles, as well as their resistance profiles to a range of antibiotics. It is evident that the newly identified hypothetical and conserved hypothetical proteins, as well as uncharacterized two-component regulatory systems and membrane proteins, may hold additional significant information on the major differences in virulence among the species, as well as the specificity of the strains for particular hosts.

Structural flexibility in the Burkholderia mallei genome

The complete genome sequence of Burkholderia mallei ATCC 23344 provides insight into this highly infectious bacterium's pathogenicity and evolutionary history. B. mallei, the etiologic agent of glanders, has come under renewed scientific investigation as a result of recent concerns about its past and potential future use as a biological weapon. Genome analysis identified a number of putative virulence factors whose function was supported by comparative genome hybridization and expression profiling of the bacterium in hamster liver in vivo. The genome contains numerous insertion sequence elements that have mediated extensive deletions and rearrangements of the genome relative to Burkholderia pseudomallei. The genome also contains a vast number (>12,000) of simple sequence repeats. Variation in simple sequence repeats in key genes can provide a mechanism for generating antigenic variation that may account for the mammalian host's inability to mount a durable adaptive immune response to a B. mallei infection.

Genome Properties: a system for the investigation of prokaryotic genetic content for microbiology, genome annotation and comparative genomics

MOTIVATION

The presence or absence of metabolic pathways and structures provide a context that makes protein annotation far more reliable. Compiling such information across microbial genomes improves the functional classification of proteins and provides a valuable resource for comparative genomics.

RESULTS

We have created a Genome Properties system to present key aspects of prokaryotic biology using standardized computational methods and controlled vocabularies. Properties reflect gene content, phenotype, phylogeny and computational analyses. The results of searches using hidden Markov models allow many properties to be deduced automatically, especially for families of proteins (equivalogs) conserved in function since their last common ancestor. Additional properties are derived from curation, published reports and other forms of evidence. Genome Properties system was applied to 156 complete prokaryotic genomes, and is easily mined to find differences between species, correlations between metabolic features and families of uncharacterized proteins, or relationships among properties.

AVAILABILITY

Genome Properties can be found at http://www.tigr.org/Genome_Properties

SUPPLEMENTARY INFORMATION

http://www.tigr.org/tigr-scripts/CMR2/genome_properties_references.spl.

Whole genome comparisons of serotype 4b and 1/2a strains of the food-borne pathogen Listeria monocytogenes reveal new insights into the core genome components of this species

The genomes of three strains of Listeria monocytogenes that have been associated with food-borne illness in the USA were subjected to whole genome comparative analysis. A total of 51, 97 and 69 strain-specific genes were identified in L.monocytogenes strains F2365 (serotype 4b, cheese isolate), F6854 (serotype 1/2a, frankfurter isolate) and H7858 (serotype 4b, meat isolate), respectively. Eighty-three genes were restricted to serotype 1/2a and 51 to serotype 4b strains. These strain- and serotype-specific genes probably contribute to observed differences in pathogenicity, and the ability of the organisms to survive and grow in their respective environmental niches. The serotype 1/2a-specific genes include an operon that encodes the rhamnose biosynthetic pathway that is associated with teichoic acid biosynthesis, as well as operons for five glycosyl transferases and an adenine-specific DNA methyltransferase. A total of 8603 and 105 050 high quality single nucleotide polymorphisms (SNPs) were found on the draft genome sequences of strain H7858 and strain F6854, respectively, when compared with strain F2365. Whole genome comparative analyses revealed that the L.monocytogenes genomes are essentially syntenic, with the majority of genomic differences consisting of phage insertions, transposable elements and SNPs.

The genome sequence of the anaerobic, sulfate-reducing bacterium Desulfovibrio vulgaris Hildenborough

Desulfovibrio vulgaris Hildenborough is a model organism for studying the energy metabolism of sulfate-reducing bacteria (SRB) and for understanding the economic impacts of SRB, including biocorrosion of metal infrastructure and bioremediation of toxic metal ions. The 3,570,858 base pair (bp) genome sequence reveals a network of novel c-type cytochromes, connecting multiple periplasmic hydrogenases and formate dehydrogenases, as a key feature of its energy metabolism. The relative arrangement of genes encoding enzymes for energy transduction, together with inferred cellular location of the enzymes, provides a basis for proposing an expansion to the 'hydrogen-cycling' model for increasing energy efficiency in this bacterium. Plasmid-encoded functions include modification of cell surface components, nitrogen fixation and a type-III protein secretion system. This genome sequence represents a substantial step toward the elucidation of pathways for reduction (and bioremediation) of pollutants such as uranium and chromium and offers a new starting point for defining this organism's complex anaerobic respiration.

Comparison of the genome of the oral pathogen Treponema denticola with other spirochete genomes

We present the complete 2,843,201-bp genome sequence of Treponema denticola (ATCC 35405) an oral spirochete associated with periodontal disease. Analysis of the T. denticola genome reveals factors mediating coaggregation, cell signaling, stress protection, and other competitive and cooperative measures, consistent with its pathogenic nature and lifestyle within the mixed-species environment of subgingival dental plaque. Comparisons with previously sequenced spirochete genomes revealed specific factors contributing to differences and similarities in spirochete physiology as well as pathogenic potential. The T. denticola genome is considerably larger in size than the genome of the related syphilis-causing spirochete Treponema pallidum. The differences in gene content appear to be attributable to a combination of three phenomena: genome reduction, lineage-specific expansions, and horizontal gene transfer. Genes lost due to reductive evolution appear to be largely involved in metabolism and transport, whereas some of the genes that have arisen due to lineage-specific expansions are implicated in various pathogenic interactions, and genes acquired via horizontal gene transfer are largely phage-related or of unknown function.

Phylogenomics of the reproductive parasite Wolbachia pipientis wMel: a streamlined genome overrun by mobile genetic elements

The complete sequence of the 1,267,782 bp genome of Wolbachia pipientis wMel, an obligate intracellular bacteria of Drosophila melanogaster, has been determined. Wolbachia, which are found in a variety of invertebrate species, are of great interest due to their diverse interactions with different hosts, which range from many forms of reproductive parasitism to mutualistic symbioses. Analysis of the wMel genome, in particular phylogenomic comparisons with other intracellular bacteria, has revealed many insights into the biology and evolution of wMel and Wolbachia in general. For example, the wMel genome is unique among sequenced obligate intracellular species in both being highly streamlined and containing very high levels of repetitive DNA and mobile DNA elements. This observation, coupled with multiple evolutionary reconstructions, suggests that natural selection is somewhat inefficient in wMel, most likely owing to the occurrence of repeated population bottlenecks. Genome analysis predicts many metabolic differences with the closely related Rickettsia species, including the presence of intact glycolysis and purine synthesis, which may compensate for an inability to obtain ATP directly from its host, as Rickettsia can. Other discoveries include the apparent inability of wMel to synthesize lipopolysaccharide and the presence of the most genes encoding proteins with ankyrin repeat domains of any prokaryotic genome yet sequenced. Despite the ability of wMel to infect the germline of its host, we find no evidence for either recent lateral gene transfer between wMel and D. melanogaster or older transfers between Wolbachia and any host. Evolutionary analysis further supports the hypothesis that mitochondria share a common ancestor with the α-Proteobacteria, but shows little support for the grouping of mitochondria with species in the order Rickettsiales. With the availability of the complete genomes of both species and excellent genetic tools for the host, the wMel–D. melanogaster symbiosis is now an ideal system for studying the biology and evolution of Wolbachia infections.

The genome sequence of Wolbachia provides insights into the origins of mitochondria, as well as the ecology and evolution of endosymbiosis

Genome of Geobacter sulfurreducens: metal reduction in subsurface environments

The complete genome sequence of Geobacter sulfurreducens, a delta-proteobacterium, reveals unsuspected capabilities, including evidence of aerobic metabolism, one-carbon and complex carbon metabolism, motility, and chemotactic behavior. These characteristics, coupled with the possession of many two-component sensors and many c-type cytochromes, reveal an ability to create alternative, redundant, electron transport networks and offer insights into the process of metal ion reduction in subsurface environments. As well as playing roles in the global cycling of metals and carbon, this organism clearly has the potential for use in bioremediation of radioactive metals and in the generation of electricity.

The genome sequence of Bacillus anthracis Ames and comparison to closely related bacteria

Bacillus anthracis is an endospore-forming bacterium that causes inhalational anthrax. Key virulence genes are found on plasmids (extra-chromosomal, circular, double-stranded DNA molecules) pXO1 (ref. 2) and pXO2 (ref. 3). To identify additional genes that might contribute to virulence, we analysed the complete sequence of the chromosome of B. anthracis Ames (about 5.23 megabases). We found several chromosomally encoded proteins that may contribute to pathogenicity--including haemolysins, phospholipases and iron acquisition functions--and identified numerous surface proteins that might be important targets for vaccines and drugs. Almost all these putative chromosomal virulence and surface proteins have homologues in Bacillus cereus, highlighting the similarity of B. anthracis to near-neighbours that are not associated with anthrax. By performing a comparative genome hybridization of 19 B. cereus and Bacillus thuringiensis strains against a B. anthracis DNA microarray, we confirmed the general similarity of chromosomal genes among this group of close relatives. However, we found that the gene sequences of pXO1 and pXO2 were more variable between strains, suggesting plasmid mobility in the group. The complete sequence of B. anthracis is a step towards a better understanding of anthrax pathogenesis.

Complete genome sequence of the Q-fever pathogen Coxiella burnetii

The 1,995,275-bp genome of Coxiella burnetii, Nine Mile phase I RSA493, a highly virulent zoonotic pathogen and category B bioterrorism agent, was sequenced by the random shotgun method. This bacterium is an obligate intracellular acidophile that is highly adapted for life within the eukaryotic phagolysosome. Genome analysis revealed many genes with potential roles in adhesion, invasion, intracellular trafficking, host-cell modulation, and detoxification. A previously uncharacterized 13-member family of ankyrin repeat-containing proteins is implicated in the pathogenesis of this organism. Although the lifestyle and parasitic strategies of C. burnetii resemble that of Rickettsiae and Chlamydiae, their genome architectures differ considerably in terms of presence of mobile elements, extent of genome reduction, metabolic capabilities, and transporter profiles. The presence of 83 pseudogenes displays an ongoing process of gene degradation. Unlike other obligate intracellular bacteria, 32 insertion sequences are found dispersed in the chromosome, indicating some plasticity in the C. burnetii genome. These analyses suggest that the obligate intracellular lifestyle of C. burnetii may be a relatively recent innovation.

The Brucella suis genome reveals fundamental similarities between animal and plant pathogens and symbionts

The 3.31-Mb genome sequence of the intracellular pathogen and potential bioterrorism agent, Brucella suis, was determined. Comparison of B. suis with Brucella melitensis has defined a finite set of differences that could be responsible for the differences in virulence and host preference between these organisms, and indicates that phage have played a significant role in their divergence. Analysis of the B. suis genome reveals transport and metabolic capabilities akin to soil/plant-associated bacteria. Extensive gene synteny between B. suis chromosome 1 and the genome of the plant symbiont Mesorhizobium loti emphasizes the similarity between this animal pathogen and plant pathogens and symbionts. A limited repertoire of genes homologous to known bacterial virulence factors were identified.

Complete genome sequence and comparative genomic analysis of an emerging human pathogen, serotype V Streptococcus agalactiae

The 2,160,267 bp genome sequence of Streptococcus agalactiae, the leading cause of bacterial sepsis, pneumonia, and meningitis in neonates in the U.S. and Europe, is predicted to encode 2,175 genes. Genome comparisons among S. agalactiae, Streptococcus pneumoniae, Streptococcus pyogenes, and the other completely sequenced genomes identified genes specific to the streptococci and to S. agalactiae. These in silico analyses, combined with comparative genome hybridization experiments between the sequenced serotype V strain 2603 V/R and 19 S. agalactiae strains from several serotypes using whole-genome microarrays, revealed the genetic heterogeneity among S. agalactiae strains, even of the same serotype, and provided insights into the evolution of virulence mechanisms.