Aye Wollam

A catalog of reference genomes from the human microbiome.

News from the Inner Tube of Life A major initiative by the U.S. National Institutes of Health to sequence 900 genomes of microorganisms that live on the surfaces and orifices of the human body has established standardized protocols and methods for such large-scale reference sequencing. By combining previously accumulated data with new data, Nelson et al. (p. 994) present an initial analysis of 178 bacterial genomes. The sampling so far barely scratches the surface of the microbial diversity found on humans, but the work provides an important baseline for future analyses. Standardized protocols and methods are being established for large-scale sequencing of the microorganisms living on humans. The human microbiome refers to the community of microorganisms, including prokaryotes, viruses, and microbial eukaryotes, that populate the human body. The National Institutes of Health launched an initiative that focuses on describing the diversity of microbial species that are associated with health and disease. The first phase of this initiative includes the sequencing of hundreds of microbial reference genomes, coupled to metagenomic sequencing from multiple body sites. Here we present results from an initial reference genome sequencing of 178 microbial genomes. From 547,968 predicted polypeptides that correspond to the gene complement of these strains, previously unidentified (“novel”) polypeptides that had both unmasked sequence length greater than 100 amino acids and no BLASTP match to any nonreference entry in the nonredundant subset were defined. This analysis resulted in a set of 30,867 polypeptides, of which 29,987 (~97%) were unique. In addition, this set of microbial genomes allows for ~40% of random sequences from the microbiome of the gastrointestinal tract to be associated with organisms based on the match criteria used. Insights into pan-genome analysis suggest that we are still far from saturating microbial species genetic data sets. In addition, the associated metrics and standards used by our group for quality assurance are presented.

Characterizing a model human gut microbiota composed of members of its two dominant bacterial phyla

The adult human distal gut microbial community is typically dominated by 2 bacterial phyla (divisions), the Firmicutes and the Bacteroidetes. Little is known about the factors that govern the interactions between their members. Here, we examine the niches of representatives of both phyla in vivo. Finished genome sequences were generated from Eubacterium rectale and E. eligens, which belong to Clostridium Cluster XIVa, one of the most common gut Firmicute clades. Comparison of these and 25 other gut Firmicutes and Bacteroidetes indicated that the Firmicutes possess smaller genomes and a disproportionately smaller number of glycan-degrading enzymes. Germ-free mice were then colonized with E. rectale and/or a prominent human gut Bacteroidetes, Bacteroides thetaiotaomicron, followed by whole-genome transcriptional profiling, high-resolution proteomic analysis, and biochemical assays of microbial–microbial and microbial–host interactions. B. thetaiotaomicron adapts to E. rectale by up-regulating expression of a variety of polysaccharide utilization loci encoding numerous glycoside hydrolases, and by signaling the host to produce mucosal glycans that it, but not E. rectale, can access. E. rectale adapts to B. thetaiotaomicron by decreasing production of its glycan-degrading enzymes, increasing expression of selected amino acid and sugar transporters, and facilitating glycolysis by reducing levels of NADH, in part via generation of butyrate from acetate, which in turn is used by the gut epithelium. This simplified model of the human gut microbiota illustrates niche specialization and functional redundancy within members of its major bacterial phyla, and the importance of host glycans as a nutrient foundation that ensures ecosystem stability.

Genome Project Standards in a New Era of Sequencing

More detailed sequence standards that keep up with revolutionary sequencing technologies will aid the research community in evaluating data. For over a decade, genome sequences have adhered to only two standards that are relied on for purposes of sequence analysis by interested third parties (1, 2). However, ongoing developments in revolutionary sequencing technologies have resulted in a redefinition of traditional whole-genome sequencing that requires reevaluation of such standards. With commercially available 454 pyrosequencing (followed by Illumina, SOLiD, and now Helicos), there has been an explosion of genomes sequenced under the moniker “draft”; however, these can be very poor quality genomes (due to inherent errors in the sequencing technologies, and the inability of assembly programs to fully address these errors). Further, one can only infer that such draft genomes may be of poor quality by navigating through the databases to find the number and type of reads deposited in sequence trace repositories (and not all genomes have this available), or to identify the number of contigs or genome fragments deposited to the database. The difficulty in assessing the quality of such deposited genomes has created some havoc for genome analysis pipelines and has contributed to many wasted hours. Exponential leaps in raw sequencing capability and greatly reduced prices have further skewed the time- and cost-ratios of draft data generation versus the painstaking process of improving and finishing a genome. The result is an ever-widening gap between drafted and finished genomes that only promises to continue (see the figure, page 236); hence, there is an urgent need to distinguish good from poor data sets.

Modernizing Reference Genome Assemblies

I have read the journals policy and have the following conflicts: Paul Flicek is married to the deputy editor of PLoS Medicine, Melissa Norton. Evan Eichler is on the board of Pacific Biosciences. Support for this work came from the Intramural Research Program of the NIH, The National Library of Medicine, the European Molecular Biology Laboratory, the Wellcome Trust (grant number 077198), and the Howard Hughes Medical Institute (EEE). The funders had no role in study design, data collection and analysis, decision to publish or preparation of the manuscript.

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.

The genome of Cyanothece 51142, a unicellular diazotrophic cyanobacterium important in the marine nitrogen cycle

Unicellular cyanobacteria have recently been recognized for their contributions to nitrogen fixation in marine environments, a function previously thought to be filled mainly by filamentous cyanobacteria such as Trichodesmium. To begin a systems level analysis of the physiology of the unicellular N2-fixing microbes, we have sequenced to completion the genome of Cyanothece sp. ATCC 51142, the first such organism. Cyanothece 51142 performs oxygenic photosynthesis and nitrogen fixation, separating these two incompatible processes temporally within the same cell, while concomitantly accumulating metabolic products in inclusion bodies that are later mobilized as part of a robust diurnal cycle. The 5,460,377-bp Cyanothece 51142 genome has a unique arrangement of one large circular chromosome, four small plasmids, and one linear chromosome, the first report of a linear element in the genome of a photosynthetic bacterium. On the 429,701-bp linear chromosome is a cluster of genes for enzymes involved in pyruvate metabolism, suggesting an important role for the linear chromosome in fermentative processes. The annotation of the genome was significantly aided by simultaneous global proteomic studies of this organism. Compared with other nitrogen-fixing cyanobacteria, Cyanothece 51142 contains the largest intact contiguous cluster of nitrogen fixation-related genes. We discuss the implications of such an organization on the regulation of nitrogen fixation. The genome sequence provides important information regarding the ability of Cyanothece 51142 to accomplish metabolic compartmentalization and energy storage, as well as how a unicellular bacterium balances multiple, often incompatible, processes in a single cell.

Transferable vancomycin resistance in a community-associated MRSA lineage.

We report the case of a patient from Brazil with a bloodstream infection caused by a strain of methicillin-resistant Staphylococcus aureus (MRSA) that was susceptible to vancomycin (designated BR-VSSA) but that acquired the vanA gene cluster during antibiotic therapy and became resistant to vancomycin (designated BR-VRSA). Both strains belong to the sequence type (ST) 8 community-associated genetic lineage that carries the staphylococcal chromosomal cassette mec (SCCmec) type IVa and the S. aureus protein A gene (spa) type t292 and are phylogenetically related to MRSA lineage USA300. A conjugative plasmid of 55,706 bp (pBRZ01) carrying the vanA cluster was identified and readily transferred to other staphylococci. The pBRZ01 plasmid harbors DNA sequences that are typical of the plasmid-associated replication genes rep24 or rep21 described in community-associated MRSA strains from Australia (pWBG745). The presence and dissemination of community-associated MRSA containing vanA could become a serious public health concern.

Whole-Genome Analyses of Enterococcus faecium Isolates with Diverse Daptomycin MICs

ABSTRACT Daptomycin (DAP) is a lipopeptide antibiotic frequently used as a “last-resort” antibiotic against vancomycin-resistant Enterococcus faecium (VRE). However, an important limitation for DAP therapy against VRE is the emergence of resistance during therapy. Mutations in regulatory systems involved in cell envelope homeostasis are postulated to be important mediators of DAP resistance in E. faecium. Thus, in order to gain insights into the genetic bases of DAP resistance in E. faecium, we investigated the presence of changes in 43 predicted proteins previously associated with DAP resistance in enterococci and staphylococci using the genomes of 19 E. faecium with different DAP MICs (range, 3 to 48 μg/ml). Bodipy-DAP (BDP-DAP) binding to the cell membrane assays and time-kill curves (DAP alone and with ampicillin) were performed. Genetic changes involving two major pathways were identified: (i) LiaFSR, a regulatory system associated with the cell envelope stress response, and (ii) YycFGHIJ, a system involved in the regulation of cell wall homeostasis. Thr120→Ala and Trp73→Cys substitutions in LiaS and LiaR, respectively, were the most common changes identified. DAP bactericidal activity was abolished in the presence of liaFSR or yycFGHIJ mutations regardless of the DAP MIC and was restored in the presence of ampicillin, but only in representatives of the LiaFSR pathway. Reduced binding of BDP-DAP to the cell surface was the predominant finding correlating with resistance in isolates with DAP MICs above the susceptibility breakpoint. Our findings suggest that genotypic information may be crucial to predict response to DAP plus β-lactam combinations and continue to question the DAP breakpoint of 4 μg/ml.

Evolutionary Genomics of Salmonella enterica Subspecies

ABSTRACT Six subspecies are currently recognized in Salmonella enterica. Subspecies I (subspecies enterica) is responsible for nearly all infections in humans and warm-blooded animals, while five other subspecies are isolated principally from cold-blooded animals. We sequenced 21 phylogenetically diverse strains, including two representatives from each of the previously unsequenced five subspecies and 11 diverse new strains from S. enterica subspecies enterica, to put this species into an evolutionary perspective. The phylogeny of the subspecies was partly obscured by abundant recombination events between lineages and a relatively short period of time within which subspeciation took place. Nevertheless, a variety of different tree-building methods gave congruent evolutionary tree topologies for subspeciation. A total of 285 gene families were identified that were recruited into subspecies enterica, and most of these are of unknown function. At least 2,807 gene families were identified in one or more of the other subspecies that are not found in subspecies I or Salmonella bongori. Among these gene families were 13 new candidate effectors and 7 new candidate fimbrial clusters. A third complete type III secretion system not present in subspecies enterica (I) isolates was found in both strains of subspecies salamae (II). Some gene families had complex taxonomies, such as the type VI secretion systems, which were recruited from four different lineages in five of six subspecies. Analysis of nonsynonymous-to-synonymous substitution rates indicated that the more-recently acquired regions in S. enterica are undergoing faster fixation rates than the rest of the genome. Recently acquired AT-rich regions, which often encode virulence functions, are under ongoing selection to maintain their high AT content. IMPORTANCE We have sequenced 21 new genomes which encompass the phylogenetic diversity of Salmonella, including strains of the previously unsequenced subspecies arizonae, diarizonae, houtenae, salamae, and indica as well as new diverse strains of subspecies enterica. We have deduced possible evolutionary paths traversed by this very important zoonotic pathogen and identified novel putative virulence factors that are not found in subspecies I. Gene families gained at the time of the evolution of subspecies enterica are of particular interest because they include mechanisms by which this subspecies adapted to warm-blooded hosts. IMPORTANCE We have sequenced 21 new genomes which encompass the phylogenetic diversity of Salmonella, including strains of the previously unsequenced subspecies arizonae, diarizonae, houtenae, salamae, and indica as well as new diverse strains of subspecies enterica. We have deduced possible evolutionary paths traversed by this very important zoonotic pathogen and identified novel putative virulence factors that are not found in subspecies I. Gene families gained at the time of the evolution of subspecies enterica are of particular interest because they include mechanisms by which this subspecies adapted to warm-blooded hosts.

Phenotypic and genotypic analysis of Clostridium difficile isolates: A single-center study

ABSTRACT Clostridium difficile infections (CDI) are a growing concern in North America, because of their increasing incidence and severity. Using integrated approaches, we correlated pathogen genotypes and host clinical characteristics for 46 C. difficile infections in a tertiary care medical center during a 6-month interval from January to June 2010. Multilocus sequence typing (MLST) demonstrated 21 known and 2 novel sequence types (STs), suggesting that the institutions C. difficile strains are genetically diverse. ST-1 (which corresponds to pulsed-field gel electrophoresis strain type NAP1/ribotype 027) was the most prevalent (32.6%); 43.5% of the isolates were binary toxin gene positive, of which 75% were ST-1. All strains were ciprofloxacin resistant and metronidazole susceptible, and 8.3% and 13.0% of the isolates were resistant to clindamycin and tetracycline, respectively. The corresponding resistance loci, including potential novel mutations, were identified from the whole-genome sequencing (WGS) of the resistant strains. Core genome single nucleotide polymorphisms (SNPs) determining the phylogenetic relatedness of the 46 strains recapitulated MLST types and provided greater interstrain differentiation. The disease severity was greatest in patients infected with ST-1 and/or binary gene-positive strains, but genome-wide SNP analysis failed to provide additional associations with CDI severity within the same STs. We conclude that MLST and core genome SNP typing result in the same phylogenetic grouping of the 46 C. difficile strains collected in a single hospital. WGS also has the capacity to differentiate those strains within STs and allows the comparison of strains at the individual gene level and at the whole-genome level.