A gene cluster for the synthesis of serotype g-specific polysaccharide antigen in Aggregatibacter actinomycetemcomitans

Aggregatibacter actinomycetemcomitans is an important pathogen related to aggressively progressive periodontal breakdown in adolescents and adults. The species can be divided into six serotypes (a–f) according to their surface carbohydrate antigens. Recently, a new serotype g of A. actinomycetemcomitans was proposed. The aim of the present study was to sequence the gene cluster associated with the biosynthesis of the serotype g-specific polysaccharide antigen and develop serotype-specific primers for PCR assay to identify serotype g strains of A. actinomycetemcomitans. The serotype-specific polysaccharide (SSPS) gene cluster of the NUM-Aa 4039 strain contained 21 genes in 21,842-bp nucleotides. The similarity of the SSPS gene cluster sequence was 96.7 % compared with that of the serotype e strain. Seventeen serotype g genes showed more than 90 % homology both in nucleotide and amino acids to the serotype e strain. Three additional genes with 1,579 bp in NUM-Aa 4039 were inserted into the corresponding ORF13 of the serotype e strain. The serotype g-specific primers were designed from the insertion region of NUM-Aa 4039. Serotypes of the a–f strains were not amplified by serotype-specific g primers; only NUM-Aa 4039 showed an amplicon band. The NUM-Aa 4039 strain was three genes in the SSPS gene cluster different from those of serotype e strain. The specific primers derived from these different regions are useful for identification and distribution of serotype g strain among A. actinomycetemcomitans from clinical samples.     

Complete Genome Analysis of Three Acinetobacter baumannii Clinical Isolates in China for Insight into the Diversification of Drug Resistance Elements

Background

The emergence and rapid spreading of multidrug-resistant Acinetobacter baumannii strains has become a major health threat worldwide. To better understand the genetic recombination related with the acquisition of drug-resistant elements during bacterial infection, we performed complete genome analysis on three newly isolated multidrug-resistant A. baumannii strains from Beijing using next-generation sequencing technology.

Methodologies/Principal Findings

Whole genome comparison revealed that all 3 strains share some common drug resistant elements including carbapenem-resistant bla _OXA-23 and tetracycline (tet) resistance islands, but the genome structures are diversified among strains. Various genomic islands intersperse on the genome with transposons and insertions, reflecting the recombination flexibility during the acquisition of the resistant elements. The blood-isolated BJAB07104 and ascites-isolated BJAB0868 exhibit high similarity on their genome structure with most of the global clone II strains, suggesting these two strains belong to the dominant outbreak strains prevalent worldwide. A large resistance island (RI) of about 121-kb, carrying a cluster of resistance-related genes, was inserted into the ATPase gene on BJAB07104 and BJAB0868 genomes. A 78-kb insertion element carrying tra-locus and bla _OXA-23 island, can be either inserted into one of the tniB gene in the 121-kb RI on the chromosome, or transformed to conjugative plasmid in the two BJAB strains. The third strains of this study, BJAB0715, which was isolated from spinal fluid, exhibit much more divergence compared with above two strains. It harbors multiple drug-resistance elements including a truncated AbaR-22-like RI on its genome. One of the unique features of this strain is that it carries both bla _OXA-23 and bla _OXA-58 genes on its genome. Besides, an Acinetobacter lwoffii adeABC efflux element was found inserted into the ATPase position in BJAB0715.

Conclusions

Our comparative analysis on currently completed Acinetobacter baumannii genomes revealed extensive and dynamic genome organizations, which may facilitate the bacteria to acquire drug-resistance elements into their genomes.     

Comparative genomic analysis of pyrene-degrading <Emphasis Type="Italic">Mycobacterium</Emphasis> species: Genomic islands and ring-hydroxylating dioxygenases involved in pyrene degradation

The genome sequences of two pyrene-degrading bacterial strains of Mycobacterium spp. PYR10 and PYR15, isolated from the estuarine wetland of the Han river, South Korea, were determined using the PacBio RS II sequencing platform. The complete genome of strain PYR15 was 6,037,017 bp in length with a GC content of 66.5%, and contained 5,933 protein-coding genes. The genome of strain PYR10 was 5,999,427 bp in length with a GC content of 67.7%, and contained 5,767 protein-coding genes. Based on the average nucleotide identity values, these strains were designated as M. gilvum PYR10 and M. pallens PYR15. A genomic comparison of these pyrene-degrading Mycobacterium strains with pyrene-non-degrading strains revealed that the genomes of pyrene-degrading strains possessed similar repertoires of ringhydroxylating dioxygenases (RHDs), including the pyrenehydroxylating dioxygenases encoded by nidA and nidA3, which could be readily distinguished from those of pyrenenon-degraders. Furthermore, genomic islands, containing catabolic gene clusters, were shared only among the pyrenedegrading Mycobacterium strains and these gene clusters contained RHD genes, including nidAB and nidA3B3. Our genome data should facilitate further studies on the evolution of the polycyclic aromatic hydrocarbon-degradation pathways in the genus Mycobacterium.     

Pathogenic and antimicrobial resistance genes in Streptococcus oralis strains revealed by comparative genome analysis

Streptococcus oralis is an early colonizer bacterium in dental plaques and is considered a potential pathogen of infective endocarditis (IE) disease. In this study, we built a complete genome map of Streptococcus oralis strain SOT, Streptococcus oralis strain SOD and Streptococcus infantis strain SO and performed comparative genomic analysis among these three strains. The results showed that there are five genomic islands (GIs) in strain SOT and one CRISPR in strain SOD. Each genome harbors various pathogenic genes related to diseases and drug resistance, while the antibiotic resistance genes in strains SOT and SOD were quite similar but different from those in strain SO. In addition, we identified 17 main virulence factors and capsule-related genes in three strains. These results suggest the pathogenic potential of Streptococcus strains, which lay a foundation for the prevention and treatment of a Streptococcus oralis infection.     

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.     

Identification of Genes and Genomic Islands Correlated with High Pathogenicity in Streptococcus suis Using Whole Genome Tilling Microarrays

Streptococcus suis is an important zoonotic pathogen that can cause meningitis and sepsis in both pigs and humans. Infections in humans have been sporadic worldwide but two severe outbreaks occurred in China in recent years, while infections in pigs are a major problem in the swine industry. Some S. suis strains are more pathogenic than others with 2 sequence types (ST), ST1 and ST7, being well recognized as highly pathogenic. We analyzed 31 isolates from 23 serotypes and 25 STs by NimbleGen tiling microarray using the genome of a high pathogenicity (HP) ST1 strain, GZ1, as reference and a new algorithm to detect gene content difference. The number of genes absent in a strain ranged from 49 to 225 with a total of 632 genes absent in at least one strain, while 1346 genes were found to be invariably present in all strains as the core genome of S. suis, accounting for 68% of the GZ1 genome. The majority of genes are located in chromosomal blocks with two or more contiguous genes. Sixty two blocks are absent in two or more strains and defined as regions of difference (RDs), among which 26 are putative genomic islands (GIs). Clustering and statistical analyses revealed that 8 RDs including 6 putative GIs and 21 genes within these RDs are significantly associated with HP. Three RDs encode known virulence related factors including the extracellular factor, the capsular polysaccharide and a SrtF pilus. The strains were divided into 5 groups based on population genetic analysis of multilocus sequence typing data and the distribution of the RDs among the groups revealed gain and loss of RDs in different groups. Our study elucidated the gene content diversity of S. suis and identified genes that potentially promote HP.     

A horizontal gene transfer event defines two distinct groups within Burkholderia pseudomallei that have dissimilar geographic distributions

Burkholderia pseudomallei is the etiologic agent of melioidosis. Many disease manifestations are associated with melioidosis, and the mechanisms causing this variation are unknown; genomic differences among strains offer one explanation. We compared the genome sequences of two strains of B. pseudomallei: the original reference strain K96243 from Thailand and strain MSHR305 from Australia. We identified a variable homologous region between the two strains. This region was previously identified in comparisons of the genome of B. pseudomallei strain K96243 with the genome of strain E264 from the closely related B. thailandensis. In that comparison, K96243 was shown to possess a horizontally acquired Yersinia-like fimbrial (YLF) gene cluster. Here, we show that the homologous genomic region in B. pseudomallei strain 305 is similar to that previously identified in B. thailandensis strain E264. We have named this region in B. pseudomallei strain 305 the B. thailandensis-like flagellum and chemotaxis (BTFC) gene cluster. We screened for these different genomic components across additional genome sequences and 571 B. pseudomallei DNA extracts obtained from regions of endemicity. These alternate genomic states define two distinct groups within B. pseudomallei: all strains contained either the BTFC gene cluster (group BTFC) or the YLF gene cluster (group YLF). These two groups have distinct geographic distributions: group BTFC is dominant in Australia, and group YLF is dominant in Thailand and elsewhere. In addition, clinical isolates are more likely to belong to group YLF, whereas environmental isolates are more likely to belong to group BTFC. These groups should be further characterized in an animal model.     

Genomic Adaptation of the Lactobacillus casei Group

Lactobacillus casei, L. paracasei, and L. rhamnosus form a closely related taxonomic group (Lactobacillus casei group) within the facultatively heterofermentative lactobacilli. Here, we report the complete genome sequences of L. paracasei JCM 8130 and L. casei ATCC 393, and the draft genome sequence of L. paracasei COM0101, all of which were isolated from daily products. Furthermore, we re-annotated the genome of L. rhamnosus ATCC 53103 (also known as L. rhamnosus GG), which we have previously reported. We confirmed that ATCC 393 is distinct from other strains previously described as L. paracasei. The core genome of 10 completely sequenced strains of the L. casei group comprised 1,682 protein-coding genes. Although extensive genome-wide synteny was found among the L. casei group, the genomes of ATCC 53103, JCM 8130, and ATCC 393 contained genomic islands compared with L. paracasei ATCC 334. Several genomic islands, including carbohydrate utilization gene clusters, were found at the same loci in the chromosomes of the L. casei group. The spaCBA pilus gene cluster, which was first identified in GG, was also found in other strains of the L. casei group, but several L. paracasei strains including COM0101 contained truncated spaC gene. ATCC 53103 encoded a higher number of proteins involved in carbohydrate utilization compared with intestinal lactobacilli, and extracellular adhesion proteins, several of which are absent in other strains of the L. casei group. In addition to previously fully sequenced L. rhamnosus and L. paracasei strains, the complete genome sequences of L. casei will provide valuable insights into the evolution of the L. casei group.     

Methylobacterium Genome Sequences: A Reference Blueprint to Investigate Microbial Metabolism of C1 Compounds from Natural and Industrial Sources

Background

Methylotrophy describes the ability of organisms to grow on reduced organic compounds without carbon-carbon bonds. The genomes of two pink-pigmented facultative methylotrophic bacteria of the Alpha-proteobacterial genus Methylobacterium, the reference species Methylobacterium extorquens strain AM1 and the dichloromethane-degrading strain DM4, were compared.

Methodology/Principal Findings

The 6.88 Mb genome of strain AM1 comprises a 5.51 Mb chromosome, a 1.26 Mb megaplasmid and three plasmids, while the 6.12 Mb genome of strain DM4 features a 5.94 Mb chromosome and two plasmids. The chromosomes are highly syntenic and share a large majority of genes, while plasmids are mostly strain-specific, with the exception of a 130 kb region of the strain AM1 megaplasmid which is syntenic to a chromosomal region of strain DM4. Both genomes contain large sets of insertion elements, many of them strain-specific, suggesting an important potential for genomic plasticity. Most of the genomic determinants associated with methylotrophy are nearly identical, with two exceptions that illustrate the metabolic and genomic versatility of Methylobacterium. A 126 kb dichloromethane utilization (dcm) gene cluster is essential for the ability of strain DM4 to use DCM as the sole carbon and energy source for growth and is unique to strain DM4. The methylamine utilization (mau) gene cluster is only found in strain AM1, indicating that strain DM4 employs an alternative system for growth with methylamine. The dcm and mau clusters represent two of the chromosomal genomic islands (AM1: 28; DM4: 17) that were defined. The mau cluster is flanked by mobile elements, but the dcm cluster disrupts a gene annotated as chelatase and for which we propose the name “island integration determinant” (iid).

Conclusion/Significance

These two genome sequences provide a platform for intra- and interspecies genomic comparisons in the genus Methylobacterium, and for investigations of the adaptive mechanisms which allow bacterial lineages to acquire methylotrophic lifestyles.     

Comparative genomic analysis of two Burkholderia glumae strains from different geographic origins reveals a high degree of plasticity in genome structure associated with genomic islands

Burkholderia glumae is the major causal agent of bacterial panicle blight of rice, a growing disease problem in global rice production. To better understand its genome-scale characteristics, the genome of the highly virulent B. glumae strain 336gr-1 isolated from Louisiana, USA was sequenced using the Illumina Genome Analyser II system. De novo assembled 336gr-1 contigs were aligned and compared with the previously sequenced genome of B. glumae strain BGR1, which was isolated from an infected rice plant in South Korea. Comparative analysis of the whole genomes of B. glumae 336gr-1 and B. glumae BGR1 revealed numerous unique genomic regions present only in one of the two strains. These unique regions contained accessory genes including mobile elements and phage-related genes, and some of the unique regions in B. glumae BGR1 corresponded to predicted genomic islands. In contrast, little variation was observed in known and potential virulence genes between the two genomes. The considerable amount of plasticity largely based on accessory genes and genome islands observed from the comparison of the genomes of these two strains of B. glumae may explain the versatility of this bacterial species in various environmental conditions and geographic locations.     

Genome Sequence of Aggregatibacter actinomycetemcomitans Serotype c Strain D11S-1

Aggregatibacter actinomycetemcomitans is a major etiological agent of periodontitis. Here we report the complete genome sequence of serotype c strain D11S-1, which was recovered from the subgingival plaque of a patient diagnosed with generalized aggressive periodontitis.Aggregatibacter actinomycetemcomitans is a major etiologic agent of human periodontal disease, in particular aggressive periodontitis (12). The natural population of A. actinomycetemcomitans is clonal (7). Six A. actinomycetemcomitans serotypes are distinguished based on the structural and serological characteristics of the O antigen of LPS (6, 7). Three of the serotypes (a, b, and c) comprise >80% of all strains, and each serotype represents a distinct clonal lineage (1, 6, 7). Serotype c strain D11S-1 was cultured from a subgingival plaque sample of a patient diagnosed with generalized aggressive periodontitis. The complete genome sequencing of the strain was determined by 454 pyrosequencing (10), which achieved 25× coverage. Assembly was performed using the Newbler assembler (454, Branford, CT) and generated 199 large contigs, with 99.3% of the bases having a quality score of 40 and above. The contigs were aligned with the genome of the sequenced serotype b strain HK1651 (http://www.genome.ou.edu/act.html) using software written in house. The putative contig gaps were then closed by primer walking and sequencing of PCR products over the gaps. The final genome assembly was further confirmed by comparison of an in silico NcoI restriction map to the experimental map generated by optical mapping (8). The genome structure of the D11S-1 strain was compared to that of the sequenced strain HK1651 using the program MAUVE (2, 3). The automated annotation was done using a protocol similar to the annotation engine service at The Institute for Genomic Research/J. Craig Venter Institute with some local modifications. Briefly, protein-coding genes were identified using Glimmer3 (4). Each protein sequence was then annotated by comparing to the GenBank nonredundant protein database. BLAST-Extend-Repraze was applied to the predicted genes to identify genes that might have been truncated due to a frameshift mutation or premature stop codon. tRNA and rRNA genes were identified by using tRNAScan-SE (9) and a similarity search to our in-house RNA database, respectively.The D11S-1 circular genome contains 2,105,764 nucleotides, a GC content of 44.55%, 2,134 predicted coding sequences, and 54 tRNA and 19 rRNA genes (see additional data at http://expression.washington.edu/bumgarnerlab/publications.php). The distribution of predicted genes based on functional categories was similar between D11S-1 and HK1651 (http://expression.washington.edu/bumgarnerlab/publications.php). One hundred six and 86 coding sequences were unique to strain D11S-1 and HK1651, respectively (http://expression.washington.edu/bumgarnerlab/publications.php). Genomic islands were identified based on annotations for strain HK1651 and based on manual inspection of contiguous D11S-1 specific DNA regions with G+C bias (http://expression.washington.edu/bumgarnerlab/publications.php). Among 12 identified genomics islands, 5 (B, C, D, E and G; cytolethal distending toxin gene cluster, tight adherence gene cluster, O-antigen biosynthesis and transport gene cluster, leukotoxin gene cluster, and lipoligosaccharide biosynthesis enzyme gene, respectively) correspond to islands 2 to 5 and 8 of strain HK1651 (http://www.oralgen.lanl.gov/) (5). Island F (∼5 kb) is homologous to a portion of the 12.5-kb island 7 in HK1651. Five genomic islands (H to L) were unique to strain D11S-1. The remaining island (A) is a fusion of genomic islands 1 and 6, in strain HK1651. The genome of D11S-1 is largely in synteny with the genome of the sequenced serotype b strain HK1651 but contained several large-scale genomic rearrangements.Strain D11S-1 harbors a 43-kb bacteriophage and two plasmids of 31 and 23 kb (http://expression.washington.edu/bumgarnerlab/publications.php). Excluding an ∼9-kb region of low homology, the phage showed >90% nucleotide sequence identity with AaΦ23 (11). A 49-bp attB site (11) was identified at coordinates 2,024,825 to 2,024,873. The location of the inserted phage was identified in the optical map of strain D11S-1 and further confirmed by PCR amplification and sequencing of the regions flanking the insertion site. A closed circular form of the phage was also detected in strain D11S-1 by PCR analysis of the phage ends. The 23-kb plasmid is homologous to pVT745 (92% nucleotide identities). The 31-kb plasmid is a novel plasmid. It has significant homologies in short regions (<2 kb) to Haemophilus influenzae biotype aegyptius plasmid pF1947 and other plasmids.     

Characterization of Bacteriophages Cp1 and Cp2, the Strain-Typing Agents for Xanthomonas axonopodis pv. citri

The strains of Xanthomonas axonopodis pv. citri, the causative agent of citrus canker, are historically classified based on bacteriophage (phage) sensitivity. Nearly all X. axonopodis pv. citri strains isolated from different regions in Japan are lysed by either phage Cp1 or Cp2; Cp1-sensitive (Cp1^s) strains have been observed to be resistant to Cp2 (Cp2^r) and vice versa. In this study, genomic and molecular characterization was performed for the typing agents Cp1 and Cp2. Morphologically, Cp1 belongs to the Siphoviridae. Genomic analysis revealed that its genome comprises 43,870-bp double-stranded DNA (dsDNA), with 10-bp 3′-extruding cohesive ends, and contains 48 open reading frames. The genomic organization was similar to that of Xanthomonas phage phiL7, but it lacked a group I intron in the DNA polymerase gene. Cp2 resembles morphologically Escherichia coli T7-like phages of Podoviridae. The 42,963-bp linear dsDNA genome of Cp2 contained terminal repeats. The Cp2 genomic sequence has 40 open reading frames, many of which did not show detectable homologs in the current databases. By proteomic analysis, a gene cluster encoding structural proteins corresponding to the class III module of T7-like phages was identified on the Cp2 genome. Therefore, Cp1 and Cp2 were found to belong to completely different virus groups. In addition, we found that Cp1 and Cp2 use different molecules on the host cell surface as phage receptors and that host selection of X. axonopodis pv. citri strains by Cp1 and Cp2 is not determined at the initial stage by binding to receptors.     

The genome of Borrelia recurrentis, the agent of deadly louse-borne relapsing fever, is a degraded subset of tick-borne Borrelia duttonii

In an effort to understand how a tick-borne pathogen adapts to the body louse, we sequenced and compared the genomes of the recurrent fever agents Borrelia recurrentis and B. duttonii. The 1,242,163–1,574,910-bp fragmented genomes of B. recurrentis and B. duttonii contain a unique 23-kb linear plasmid. This linear plasmid exhibits a large polyT track within the promoter region of an intact variable large protein gene and a telomere resolvase that is unique to Borrelia. The genome content is characterized by several repeat families, including antigenic lipoproteins. B. recurrentis exhibited a 20.4% genome size reduction and appeared to be a strain of B. duttonii, with a decaying genome, possibly due to the accumulation of genomic errors induced by the loss of recA and mutS. Accompanying this were increases in the number of impaired genes and a reduction in coding capacity, including surface-exposed lipoproteins and putative virulence factors. Analysis of the reconstructed ancestral sequence compared to B. duttonii and B. recurrentis was consistent with the accelerated evolution observed in B. recurrentis. Vector specialization of louse-borne pathogens responsible for major epidemics was associated with rapid genome reduction. The correlation between gene loss and increased virulence of B. recurrentis parallels that of Rickettsia prowazekii, with both species being genomic subsets of less-virulent strains.     

Genome sequence analysis of the Indian strain Mannheimia haemolytica serotype A2 from ovine pneumonic pasteurellosis

Mannheimia haemolytica is a leading causative agent of pasteurellosis in ruminants. Genome of M. haemolytica strains from different hosts has been sequenced worldwide to understand its pathogenesis. There are only few reports on the isolation of M. haemolytica in India with limited information on its molecular characteristics. The present study focuses on genome sequence analysis of a M. haemolytica strain isolated from pneumonic sheep. Mannheimia haemolytica A2 strain NIVEDI/MH/1 was isolated and identified by species and serotype-specific PCRs. Whole genome sequencing was performed using the Ion Torrent Personal Genome Machine. A comparative genomic analysis was performed to understand the virulence determinants of the Indian strain and its phylogenetic relationship with other global strains. Sequence data revealed a draft genome of 2,211,426 bp size with 41.3% GC content, assembled into 17 contigs, and contained 2379 genes. Five genomic islands identified in the genome showed high sequence identity with other respiratory pathogens of the Pasteurellaceae family. Phylogenetic analysis showed M. haemolytica A2 NIVEDI/MH/1 is very close to a M. haemolytica A2 strain from pneumonic calf. Further, the analysis revealed the presence of virulence, metal-, and multidrug resistance genes needed for pathogenesis and survival of the bacteria during infection. Also, we identified the presence of type I-C and type II-C of CRISPR-Cas arrays in the present sequenced genome. The study emphasizes the role of M. haemolytica in respiratory infections of ruminants in the Indian subcontinent and indicates the role of vertical and horizontal gene pools in pathogenicity and survivability of the bacteria.     

A Large Gene Cluster Encoding Several Magnetosome Proteins Is Conserved in Different Species of Magnetotactic Bacteria

In magnetotactic bacteria, a number of specific proteins are associated with the magnetosome membrane (MM) and may have a crucial role in magnetite biomineralization. We have cloned and sequenced the genes of several of these polypeptides in the magnetotactic bacterium Magnetospirillum gryphiswaldense that could be assigned to two different genomic regions. Except for mamA, none of these genes have been previously reported to be related to magnetosome formation. Homologous genes were found in the genome sequences of M. magnetotacticum and magnetic coccus strain MC-1. The MM proteins identified display homology to tetratricopeptide repeat proteins (MamA), cation diffusion facilitators (MamB), and HtrA-like serine proteases (MamE) or bear no similarity to known proteins (MamC and MamD). A major gene cluster containing several magnetosome genes (including mamA and mamB) was found to be conserved in all three of the strains investigated. The mamAB cluster also contains additional genes that have no known homologs in any nonmagnetic organism, suggesting a specific role in magnetosome formation.     

Complete Genome Sequence of the Biocontrol Strain Pseudomonas protegens Cab57 Discovered in Japan Reveals Strain-Specific Diversity of This Species

The biocontrol strain Pseudomonas sp. Cab57 was isolated from the rhizosphere of shepherd’s purse growing in a field in Hokkaido by screening the antibiotic producers. The whole genome sequence of this strain was obtained by paired-end and whole-genome shotgun sequencing, and the gaps between the contigs were closed using gap-spanning PCR products. The P. sp. Cab57 genome is organized into a single circular chromosome with 6,827,892 bp, 63.3% G+C content, and 6,186 predicted protein-coding sequences. Based on 16S rRNA gene analysis and whole genome analysis, strain Cab57 was identified as P. protegens. As reported in P. protegens CHA0 and Pf-5, four gene clusters (phl, prn, plt, and hcn) encoding the typical antibiotic metabolites and the reported genes associated with Gac/Rsm signal transduction pathway of these strains are fully conserved in the Cab57 genome. Actually strain Cab57 exhibited typical Gac/Rsm activities and antibiotic production, and these activities were enhanced by knocking out the retS gene (for a sensor kinase acting as an antagonist of GacS). Two large segments (79 and 115 kb) lacking in the Cab57 genome, as compared with the Pf-5 genome, accounted for the majority of the difference (247 kb) between these genomes. One of these segments was the complete rhizoxin analog biosynthesis gene cluster (ca. 79 kb) and another one was the 115-kb mobile genomic island. A whole genome comparison of those relative strains revealed that each strain has unique gene clusters involved in metabolism such as nitrite/nitrate assimilation, which was identified in the Cab57 genome. These findings suggest that P. protegens is a ubiquitous bacterium that controls its biocontrol traits while building up strain-specific genomic repertoires for the biosynthesis of secondary metabolites and niche adaptation.     

Whole-Genome Single-Nucleotide-Polymorphism Analysis for Discrimination of Clostridium botulinum Group I Strains

Clostridium botulinum is a genetically diverse Gram-positive bacterium producing extremely potent neurotoxins (botulinum neurotoxins A through G [BoNT/A-G]). The complete genome sequences of three strains harboring only the BoNT/A1 nucleotide sequence are publicly available. Although these strains contain a toxin cluster (HA⁺ OrfX⁻) associated with hemagglutinin genes, little is known about the genomes of subtype A1 strains (termed HA⁻ OrfX⁺) that lack hemagglutinin genes in the toxin gene cluster. We sequenced the genomes of three BoNT/A1-producing C. botulinum strains: two strains with the HA⁺ OrfX⁻ cluster (69A and 32A) and one strain with the HA⁻ OrfX⁺ cluster (CDC297). Whole-genome phylogenic single-nucleotide-polymorphism (SNP) analysis of these strains along with other publicly available C. botulinum group I strains revealed five distinct lineages. Strains 69A and 32A clustered with the C. botulinum type A1 Hall group, and strain CDC297 clustered with the C. botulinum type Ba4 strain 657. This study reports the use of whole-genome SNP sequence analysis for discrimination of C. botulinum group I strains and demonstrates the utility of this analysis in quickly differentiating C. botulinum strains harboring identical toxin gene subtypes. This analysis further supports previous work showing that strains CDC297 and 657 likely evolved from a common ancestor and independently acquired separate BoNT/A1 toxin gene clusters at distinct genomic locations.     

Lactococcus garvieae: Where Is It From? A First Approach to Explore the Evolutionary History of This Emerging Pathogen

The population structure and diversity of Lactococcus garvieae, an emerging pathogen of increasing clinical significance, was determined at both gene and genome level. Selected lactococcal isolates of various origins were analyzed by a multi locus sequence typing (MLST). This gene-based analysis was compared to genomic characteristics, estimated through the complete genome sequences available in database. The MLST identified two branches containing the majority of the strains and two branches bearing one strain each. One strain was particularly differentiated from the other L. garvieae strains, showing a significant genetic distance. The genomic characteristics, correlated to the MLST-based phylogeny, indicated that this “separated strain” appeared first and could be considered the evolutionary intermediate between Lactococcus lactis and L. garvieae main clusters. A preliminary genome analysis of L. garvieae indicated a pan-genome constituted of about 4100 genes, which included 1341 core genes and 2760 genes belonging to the dispensable genome. A total of 1491 Clusters of Orthologous Genes (COGs) were found to be specific to the 11 L. garvieae genomes, with the genome of the “separated strain” showing the highest presence of unique genes.     

Strain-specific genes of Helicobacter pylori: genome evolution driven by a novel type IV secretion system and genomic island transfer

The availability of multiple bacterial genome sequences has revealed a surprising extent of variability among strains of the same species. The human gastric pathogen Helicobacter pylori is known as one of the most genetically diverse species. We have compared the genome sequence of the duodenal ulcer strain P12 and six other H. pylori genomes to elucidate the genetic repertoire and genome evolution mechanisms of this species. In agreement with previous findings, we estimate that the core genome comprises about 1200 genes and that H. pylori possesses an open pan-genome. Strain-specific genes are preferentially located at potential genome rearrangement sites or in distinct plasticity zones, suggesting two different mechanisms of genome evolution. The P12 genome contains three plasticity zones, two of which encode type IV secretion systems and have typical features of genomic islands. We demonstrate for the first time that one of these islands is capable of self-excision and horizontal transfer by a conjugative process. We also show that excision is mediated by a protein of the XerD family of tyrosine recombinases. Thus, in addition to its natural transformation competence, conjugative transfer of genomic islands has to be considered as an important source of genetic diversity in H. pylori.