New Insights into How Yersinia pestis Adapts to Its Mammalian Host during Bubonic Plague

Bubonic plague (a fatal, flea-transmitted disease) remains an international public health concern. Although our understanding of the pathogenesis of bubonic plague has improved significantly over the last few decades, researchers have still not been able to define the complete set of Y. pestis genes needed for disease or to characterize the mechanisms that enable infection. Here, we generated a library of Y. pestis mutants, each lacking one or more of the genes previously identified as being up-regulated in vivo. We then screened the library for attenuated virulence in rodent models of bubonic plague. Importantly, we tested mutants both individually and using a novel, “per-pool” screening method that we have developed. Our data showed that in addition to genes involved in physiological adaption and resistance to the stress generated by the host, several previously uncharacterized genes are required for virulence. One of these genes (ympt1.66c, which encodes a putative helicase) has been acquired by horizontal gene transfer. Deletion of ympt1.66c reduced Y. pestis'' ability to spread to the lymph nodes draining the dermal inoculation site – probably because loss of this gene decreased the bacteria''s ability to survive inside macrophages. Our results suggest that (i) intracellular survival during the early stage of infection is important for plague and (ii) horizontal gene transfer was crucial in the acquisition of this ability.     

Draft Genome Sequence of Yersinia pestis Strain 2501, an Isolate from the Great Gerbil Plague Focus in Xinjiang, China

We deciphered the genome of Yersinia pestis strain 2501, isolated from the Junggar Basin, a newly discovered great gerbil plague focus in Xinjiang, China. The total length of assembly was 4,597,322 bp, and 4,265 coding sequences were predicted within the genome. It is the first Y. pestis genome from this plague focus.     

Genotyping and Phylogenetic Analysis of Yersinia pestis by MLVA: Insights into the Worldwide Expansion of Central Asia Plague Foci

Background

The species Yersinia pestis is commonly divided into three classical biovars, Antiqua, Medievalis, and Orientalis, belonging to subspecies pestis pathogenic for human and the (atypical) non-human pathogenic biovar Microtus (alias Pestoides) including several non-pestis subspecies. Recent progress in molecular typing methods enables large-scale investigations in the population structure of this species. It is now possible to test hypotheses about its evolution which were proposed decades ago. For instance the three classical biovars of different geographical distributions were suggested to originate from Central Asia. Most investigations so far have focused on the typical pestis subspecies representatives found outside of China, whereas the understanding of the emergence of this human pathogen requires the investigation of strains belonging to subspecies pestis from China and to the Microtus biovar.

Methodology/Principal Findings

Multi-locus VNTR analysis (MLVA) with 25 loci was performed on a collection of Y. pestis isolates originating from the majority of the known foci worldwide and including typical rhamnose-negative subspecies pestis as well as rhamnose-positive subspecies pestis and biovar Microtus. More than 500 isolates from China, the Former Soviet Union (FSU), Mongolia and a number of other foci around the world were characterized and resolved into 350 different genotypes. The data revealed very close relationships existing between some isolates from widely separated foci as well as very high diversity which can conversely be observed between nearby foci.

Conclusions/Significance

The results obtained are in full agreement with the view that the Y. pestis subsp. pestis pathogenic for humans emerged in the Central Asia region between China, Kazakhstan, Russia and Mongolia, only three clones of which spread out of Central Asia. The relationships among the strains in China, Central Asia and the rest of the world based on the MLVA25 assay provide an unprecedented view on the expansion and microevolution of Y. pestis.     

Ecological Aspects of Evolution of the Plague Microbe Yersinia pestis and the Genesis of Natural Foci

A new hypothesis of the origin of the plague microbe in the Mongolian bobak (Marmota sibirica Radde, 1862) populations in Central Asia during the Pleistocene is based on the ideas of its relative phylogenetic recency. The Late Pleistocene cooling, which induced a deep freezing of the grounds in southern Siberia, Mongolia, and Manchuria, is considered as an inducer of speciation. The main ecological factors of the plague microbe evolution include the species specific behavior of the Mongolian bobak as it prepared to hibernate related to its occurrence in arid petrophytic landscapes and the larval parasitism of the flea Oropsylla silantiewi Wagn., 1898 in winter. Genesis of the plague foci is divided into two periods: natural-historical and biosocial. During the first period, the primary natural foci in Eurasia were formed and, during the second period, synanthropic (rat) and secondary natural foci appeared with the participation of humans in Africa, The New World, and on some tropical islands.     

Draft Genome Sequences of Yersinia pestis Isolates from Natural Foci of Endemic Plague in China

To gain insights into the evolutionary origin, emergence, and pathogenicity of the etiologic agent of plague, we have sequenced the genomes of four Yersinia pestis strains isolated from the zoonotic rodent reservoir in foci of endemic plague in China. These resources enable in-depth studies of Y. pestis sequence variations and detailed whole-genome comparisons of very closely related genomes from the supposed site of the origin and the emergence of global pandemics of plague.Here we report on the genomes of Yersinia pestis strains B42003004, K1973002, E1979001, and F1991016, which represent a sample of the genetic diversity found in four foci of endemic plague in China (24). Y. pestis bv. orientalis strain F1991016 was isolated in 1991 from Cangyuan County, China, from a rat (Rattus flavipectus), and Y. pestis bv. antiqua strain E1979001 was isolated in 1979 from Jianchuan, China, from a vole (Eothenomys miletus). Both Y. pestis strains K1973002 and B42003004 of biovars medievalis and antiqua, respectively, originate from marmota species (Marmota himalayana Hetian 1973; Marmota baibacina Wenquan 2003) (24). Genome analyses of these key isolates outline the details of microevolution of the plague bacterium, as these isolates represent important evolutionary milestones of the species, which is thought to have originated in Central Asia as a clonal descendant of Yersinia pseudotuberculosis (1). Genomic DNA was subjected to whole-genome shotgun sequencing and closure strategies as previously described (15). Plasmid (pHOS2) and fosmid (pCC1fos) libraries were constructed, with insert sizes of 4 to 6 kb and 30 to 40 kb, respectively. An average of 67,000 high-quality Sanger reads (total, 268,160) was obtained with an 860-bp average read length. The genomes with an average 12-fold read coverage depth were assembled using a Celera Assembler (11) and manually annotated using Manatee (http://manatee.sourceforge.net/). Genomic architectures were compared using Mauve (5, 18), and proteomes were analyzed with the BLAST score ratio tool (17).The young evolutionary history of the species and resulting homogenous population structure is reflected in a high degree of proteome conservation between the sequenced isolates and the modern strain CO92 (16). Y. pestis pathogenicity is anchored in its mobile inventory, and typically, isolates harbor three virulence plasmids, the species-specific plasminogen activator and murine toxin plasmids and the low-calcium-response plasmid pCD (23). Their pCD-borne lcrV antigen shows a genetic makeup identical to that of CO92 (2, 16). The insertion sequence element expansion clearly distinguishes these Central Asian isolates from the progenitor Y. pseudotuberculosis (3, 8). Comprehensive analyses reveal a lack of genome-wide synteny and suggest massive intrachromosomal rearrangements, a characteristic feature of Y. pestis genome evolution (6, 8). Besides insertion sequence element abundance, we observed isolate-specific propagation patterns that not only shaped the reorganization of the genomic architecture but also are known to drive microevolutionary adaptation in Y. pestis (4, 9, 14, 21, 24). Based upon the phenotypic and genotypic features that differentiate these isolates (13, 20, 24), B42003004 belongs to the most ancient Y. pestis lineage known to exist in China; hence, it is phylogenetically thought to be closest to the species progenitor Y. pseudotuberculosis (22). We studied metabolic genes that determine their biovar classification and investigated the underlying genetic determinants (24). Isolate K1973002 is defective in the nitrate reductase napA gene, similar to strain KIM (7), and represents the results of the evolutionary processes implicated in the biovar conversion from antiqua to medievalis. Isolate F1991016 carries an in-frame deletion in the glycerol-3-phosphate dehydrogenase glpD gene (19), similar to strain CO92 (16), and characteristic of the antiqua-to-orientalis conversion. The observed genetic traits strengthen the hypothesis that biovars medievalis and orientalis arose through parallel evolution from a glycerol- and nitrate-positive antiqua progenitor due to the acquisition of independent mutations (1, 10, 14). Variable-number tandem-nucleotide-repeat alleles (12) (allele K, K1973002; allele K, B42003004; allele P, E1979001; allele G, F1991016) are not biovar specific and are not discriminative enough to differentiate these isolates, which clearly supports a population-based phylogeny, as introduced by Achtman et al. (1).The whole-genome draft sequences of these evolutionary key isolates of Y. pestis will facilitate additional bioinformatic and phylogenetic analyses. The availability of high-quality Sanger sequences is crucial to resolve the genetically homogenous population structure and to shed light on Y. pestis speciation. Understanding the plasticity and genome dynamics further aids in forensic and epidemiological analyses by setting up the basis for an accurate and robust typing system for plague surveillance and promotes diagnostics development and control measures.     

Genome Sequence of the Solvent-Producing Bacterium Clostridium carboxidivorans Strain P7T

Clostridium carboxidivorans strain P7^T is a strictly anaerobic acetogenic bacterium that produces acetate, ethanol, butanol, and butyrate. The C. carboxidivorans genome contains all the genes for the carbonyl branch of the Wood-Ljungdahl pathway for CO₂ fixation, and it encodes enzymes for conversion of acetyl coenzyme A into butanol and butyrate.Clostridium carboxidivorans strain P7^T (equivalent to ATCC BAA-624^T and DSM 15243^T) is an obligate anaerobe that can grow autotrophically with H₂ and CO₂ or CO (fixing carbon via the Wood-Ljungdahl pathway), or it can grow chemoorganotrophically with simple sugars (1). Acetate, ethanol, butanol, and butyrate are end products of metabolism.For slow-growing strict anaerobes such as Clostridium carboxidivorans, genome sequencing provides a rapid theoretical characterization of its metabolism compared to traditional methods. We isolated and amplified genomic C. carboxidivorans DNA using the Wizard genomic DNA purification kit (Promega, Madison, WI) and the REPLI-g kit (Qiagen). A single shotgun pyrosequencing run using a Genome Sequencer FLX system (454 Life Sciences, Branford, CT) resulted in 429,680 high-quality reads (mean read length, 231.6 bp) that were assembled using Newbler software (454 Life Sciences) into 225 contigs >500 bp long. Paired-end sequencing produced 111,154 reads (mean read length, 256.3 bp). Assembly of the paired-end and shotgun reads produced 73 scaffolds containing 216 large contigs with a mean sequence depth of 16.33 reads. PCR amplification and Sanger sequencing were conducted, followed by scaffold assembly using Sequencher (Gene Codes, Ann Arbor, MI). The 4.4-Mb final assembly has 33 scaffolds containing 69 contigs with a Phred-equivalent quality score of 40 or above (accuracy, >99.99%) (GenBank accession no. {"type":"entrez-nucleotide","attrs":{"text":"ADEK00000000","term_id":"295883218","term_text":"ADEK00000000"}}ADEK00000000).The sequence was annotated using Annotation Engine (J. Craig Venter Institute) and manually curated using Manatee (http://manatee.sourceforge.net/). The genome has 29.7% G+C content and contains 4,174 protein-coding sequences, 3 rRNA operons, 1 tmRNA (dual tRNA-like and mRNA-like nature), 6 noncoding RNAs (ncRNAs), and 48 tRNA genes. (6). Comparison of 16S rRNA genes showed that C. carboxidivorans is closely related to Clostridium scatologenes ATCC 25775^T (97% sequence identity) and Clostridium drakei type strain SL1^T (99% sequence identity). C. carboxidivorans shares 94% 16S rRNA sequence identity with Clostridium ljungdahlii (4.6 Mb), another solventogenic species.Pathway analyses indicated that C. carboxidivorans is similar to other anaerobic acetogens, such as Moorella thermoacetica (8), in having an incomplete reductive tricarboxylic acid (TCA) cycle where fumarate reductase is absent. Like other acetogenic clostridia, C. carboxidivorans uses the Wood-Ljungdahl pathway for fixing carbon dioxide to organic carbon via acetyl coenzyme A (acetyl-CoA) (5). Two of these genes encode carbon monoxide dehydrogenase (CODH) and acetyl-CoA synthase (ACS), which form a complex to catalyze the carbonyl branch of the pathway for carbon fixation and acetyl-CoA production. C. carboxidivorans has genes that encode phosphotransacetylase and acetate kinase for converting acetyl-CoA into acetate, yielding ATP (2).C. carboxidivorans is unique among other known acetogenic clostridia because it can fix carbon via the Wood-Ljungdahl pathway and convert acetyl-CoA into butanol, which is more energy dense than ethanol. Both C. carboxidivorans and Clostridium acetobutylicum encode NADPH-dependent butanol dehydrogenase (74% identity) to convert acetyl-CoA into butanol (3, 4), but C. acetobutylicum cannot fix CO₂ or CO into acetyl-CoA. Conversely, C. ljungdahlii can fix CO and CO₂, but it lacks butanol dehydrogenase and cannot convert acetyl-CoA into butanol. Therefore, P7 includes beneficial properties of both these industrially important strains. The genome sequence of C. carboxidivorans P7 could potentially accelerate research allowing its industrial application for biofuel production or to enable some of its pathways to be used directly in synthetic biology for biofuel production.     

Draft Genome Sequence of the Cyanide-Utilizing Bacterium Pseudomonas fluorescens Strain NCIMB 11764

We report here the 6.97-Mb draft genome sequence of Pseudomonas fluorescens strain NCIMB 11764, which is capable of growth on cyanide as the sole nitrogen source. The draft genome sequence allowed the discovery of several genes implicated in enzymatic cyanide turnover and provided additional information contributing to a better understanding of this organism''s unique cyanotrophic ability. This is the first sequenced genome of a cyanide-assimilating bacterium.     

Genome Sequence of the Aerobic Bacterium Bacillus sp. Strain FJAT-13831

Bacillus sp. strain FJAT-13831 was isolated from the no. 1 pit soil of Emperor Qin''s Terracotta Warriors in Xi''an City, People''s Republic of China. The isolate showed a close relationship to the Bacillus cereus group. The draft genome sequence of Bacillus sp. FJAT-13831 was 4,425,198 bp in size and consisted of 5,567 genes (protein-coding sequences [CDS]) with an average length of 782 bp and a G+C value of 36.36%.     

Completed Genome Sequence of the Anaerobic Iron-Oxidizing Bacterium Acidovorax ebreus Strain TPSY

Acidovorax ebreus strain TPSY is the first anaerobic nitrate-dependent Fe(II) oxidizer for which there is a completed genome sequence. Preliminary protein annotation revealed an organism optimized for survival in a complex environmental system. Here, we briefly report the completed and annotated genome sequence of strain TPSY.Microorganisms from diverse anoxic environments are capable of nitrate-dependent Fe(II) oxidation at circumneutral pH (4, 11, 17, 18, 20, 21). Despite their geochemical importance (22), little is known of the underlying biochemical and genetic mechanisms. Genome sequencing of several nitrate-dependent Fe(II) oxidizers will provide insight into this process. By comparing Fe(II) oxidation mechanisms in various organisms, we hope to identify both the conserved and disparate aspects of the metabolism. The genome of Acidovorax ebreus strain TPSY is the first of these to be sequenced.Strain TPSY is a motile, Gram-negative facultative anaerobe isolated from groundwater collected from the U.S. Department of Energy Integrated Field Research Challenge site at Oak Ridge, TN. Growth experiments performed as previously described (21) revealed TPSY''s incapacity for lithoautotrophic growth, which was supported by a lack of genes in the genome encoding any known CO₂ fixation pathways. TPSY did grow mixotrophically with Fe(II) as the electron donor and a 0.1 mM acetate carbon source. 16S rRNA gene sequence analysis placed TPSY in the class Betaproteobacteria with 99.8% similarity to Acidovorax sp. strain JS42 in the family Comamonadaceae.The completed genome consisted of a single circular chromosome of 3,796,573 bp with an average 66.8% G+C content. A total of 3,479 protein-encoding genes were predicted, and 34 (0.98%) had no similarity to public database sequences. Sequencing performed at the Department of Energy Joint Genome Institute (JGI) used Sanger sequencing and 454 pyrosequencing to a depth of 20× coverage. All JGI library construction and sequencing techniques can be found at http://www.jgi.doe.gov/. Sequence assembly, quality assessment, and annotation were performed using the software Phred/Phrap/Consed (www.phrap.com) (6-8), Dupfinisher (10), CRITICA (2), GLIMMER, and GENERATION (5) and the JGI Integrated Microbial Genomes site (12). The completed genome sequence contained 33,341 reads and had an average of ninefold coverage per base and an error rate of <1 in 100,000.TPSY was named in part for its meandering motility, and its genome confirmed the twitching phenotype with the presence of pilT, pilU, and a complete set of flagellar and chemotaxis genes. The ability of TPSY to oxidize simple alcohols and acids with oxygen or nitrate respiration was confirmed by the genome. In addition, biosynthetic pathways for all amino acids except tyrosine and phenylalanine were present. No homologues of chorismate mutase (EC 5.4.99.5), an enzyme required for tyrosine and phenylalanine anabolism, were identified. The genome contained both intact Embden-Meyerhof-Parnas and Entner-Doudoroff pathways, in addition to a pentose phosphate pathway and a trichloroacetic acid cycle.In support of its facultative anaerobicity, a complete set of genes for denitrification and three different terminal oxidases (cytochrome aa₃, cbb₃, and cytochrome d quinol oxidase) were present. The cbb₃ and cytochrome d oxidases, with their high oxygen affinity, putatively enable survival in microaerobic environments (14).TPSY had sequences encoding 30 transposases, 11 integrases, and 11 phage/prophage-related genes. A region of particular interest putatively conferred resistance to lead, arsenate, and mercury: pbrRATARTBC, arsRDAB, and merRPCADE. Evidence suggests horizontal transfer and insertion of this region, as it was flanked on the 5′ end by λ prophage-related genes and the 3′ end encoded a putative Tn21 transposase. Phenotypic studies by the method of Wang et al. (19) revealed MICs of 16 μM phenylmercuric acetate and 250 μM MgCl₂. TPSY was also capable of growth in the presence of arsenate (10 mM) but did not use it as an electron acceptor.Related to phage infection, one CRISPR (clustered, regularly interspaced, short palindromic repeats) region (3, 16) was predicted. The core proteins, the cas1 and cas2 genes, and a csn1 gene formed the CRISPR subtype Nmeni, which is associated with vertebrate pathogens and commensals (9). However, the lack of typical pathogenic type I or III secretion systems such as the hec cluster of Dickeya chrysanthemi (15) or the inv/spa system of Salmonella enterica serovar Typhimurium (13) indicated that TPSY would probably not exhibit a pathogenic lifestyle.     

Complete Genome Sequence of the Probiotic Bacterium Bifidobacterium bifidum Strain BGN4

Bifidobacterium bifidum, a common endosymbiotic inhabitant of the human gut, is considered a prominent probiotic microorganism that may promote health. We completely decrypted the 2.2-Mb genome sequence of B. bifidum BGN4, a strain that had been isolated from the fecal sample of a healthy breast-fed infant, and annotated 1,835 coding sequences.     

Draft Genome Sequence of the Biocontrol Bacterium Bacillus amyloliquefaciens Strain M27

Bacillus amyloliquefaciens strain M27 is a biocontrol agent with antagonistic activities against a wide range of fungal pathogens. Here we present the 3.86-Mb draft genome sequence of the bacterium with the aims of providing insights into the genomic basis of its antifungal mechanism and facilitating its application in the biocontrol of plant diseases.     

The Solanum commersonii Genome Sequence Provides Insights into Adaptation to Stress Conditions and Genome Evolution of Wild Potato Relatives

Here, we report the draft genome sequence of Solanum commersonii, which consists of ∼830 megabases with an N50 of 44,303 bp anchored to 12 chromosomes, using the potato (Solanum tuberosum) genome sequence as a reference. Compared with potato, S. commersonii shows a striking reduction in heterozygosity (1.5% versus 53 to 59%), and differences in genome sizes were mainly due to variations in intergenic sequence length. Gene annotation by ab initio prediction supported by RNA-seq data produced a catalog of 1703 predicted microRNAs, 18,882 long noncoding RNAs of which 20% are shown to target cold-responsive genes, and 39,290 protein-coding genes with a significant repertoire of nonredundant nucleotide binding site-encoding genes and 126 cold-related genes that are lacking in S. tuberosum. Phylogenetic analyses indicate that domesticated potato and S. commersonii lineages diverged ∼2.3 million years ago. Three duplication periods corresponding to genome enrichment for particular gene families related to response to salt stress, water transport, growth, and defense response were discovered. The draft genome sequence of S. commersonii substantially increases our understanding of the domesticated germplasm, facilitating translation of acquired knowledge into advances in crop stability in light of global climate and environmental changes.     

Genome Sequence of the Agar-Degrading Marine Bacterium Alteromonadaceae sp. Strain G7

Here, we present the high-quality draft genome sequence of the agar-degrading marine gammaproteobacterium Alteromonadaceae sp. strain G7, which was isolated from coastal seawater to be utilized as a bioresource for production of agar-derived biofuels. The 3.91-Mb genome contains a number of genes encoding algal polysaccharide-degrading enzymes such as agarases and sulfatases.     

Complete Genome Sequence of the Endophytic Bacterium Burkholderia sp. Strain KJ006

Endophytes live inside plant tissues without causing any harm and may even benefit plants. Here, we provide the high-quality genome sequence of Burkholderia sp. strain KJ006, an endophytic bacterium of rice with antifungal activity. The 6.6-Mb genome, consisting of three chromosomes and a single plasmid, contains genes related to plant growth promotion or degradation of aromatic compounds.     

Complete Genome Sequence of the Pyrene-Degrading Bacterium Cycloclasticus sp. Strain P1

Cycloclasticus sp. strain P1 was isolated from deep-sea sediments of the Pacific Ocean and characterized as a unique bacterium in the degradation of pyrene, a four-ring polycyclic aromatic hydrocarbon (PAH). Here we report the complete genome of P1 and genes associated with PAH degradation.     

The Genome of Polaromonas sp. Strain JS666: Insights into the Evolution of a Hydrocarbon- and Xenobiotic-Degrading Bacterium, and Features of Relevance to Biotechnology

Polaromonas sp. strain JS666 can grow on cis-1,2-dichloroethene (cDCE) as a sole carbon and energy source and may be useful for bioremediation of chlorinated solvent-contaminated sites. Analysis of the genome sequence of JS666 (5.9 Mb) shows a bacterium well adapted to pollution that carries many genes likely to be involved in hydrocarbon and xenobiotic catabolism and metal resistance. Clusters of genes coding for haloalkane, haloalkanoate, n-alkane, alicyclic acid, cyclic alcohol, and aromatic catabolism were analyzed in detail, and growth on acetate, catechol, chloroacetate, cyclohexane carboxylate, cyclohexanol, ferulate, heptane, 3-hydroxybenzoate, hydroxyquinol, gentisate, octane, protocatechuate, and salicylate was confirmed experimentally. Strain JS666 also harbors diverse putative mobile genetic elements, including retrons, inteins, a miniature inverted-repeat transposable element, insertion sequence transposases from 14 families, eight genomic islands, a Mu family bacteriophage, and two large (338- and 360-kb) plasmids. Both plasmids are likely to be self-transferable and carry genes for alkane, alcohol, aromatic, and haloacid metabolism. Overall, the JS666 genome sequence provides insights into the evolution of pollutant-degrading bacteria and provides a toolbox of catabolic genes with utility for biotechnology.     

Genome Sequence of the Alkane-Degrading Bacterium Alcanivorax hongdengensis Type Strain A-11-3

Alcanivorax hongdengensis A-11-3^T was isolated from an oil-enriched consortium enriched from the surface seawater of Hong-Deng dock in the Straits of Malacca and Singapore. Strain A-11-3^T can degrade n-alkane and produce a lipopeptide biosurfactant. Here we report the genome of A-11-3^T and the genes associated with alkane degradation.