Characterization of the core and accessory genomes of Pseudomonas aeruginosa using bioinformatic tools Spine and AGEnt

Background

Pseudomonas aeruginosa is an important opportunistic pathogen responsible for many infections in hospitalized and immunocompromised patients. Previous reports estimated that approximately 10% of its 6.6 Mbp genome varies from strain to strain and is therefore referred to as “accessory genome”. Elements within the accessory genome of P. aeruginosa have been associated with differences in virulence and antibiotic resistance. As whole genome sequencing of bacterial strains becomes more widespread and cost-effective, methods to quickly and reliably identify accessory genomic elements in newly sequenced P. aeruginosa genomes will be needed.

Results

We developed a bioinformatic method for identifying the accessory genome of P. aeruginosa. First, the core genome was determined based on sequence conserved among the completed genomes of twelve reference strains using Spine, a software program developed for this purpose. The core genome was 5.84 Mbp in size and contained 5,316 coding sequences. We then developed an in silico genome subtraction program named AGEnt to filter out core genomic sequences from P. aeruginosa whole genomes to identify accessory genomic sequences of these reference strains. This analysis determined that the accessory genome of P. aeruginosa ranged from 6.9-18.0% of the total genome, was enriched for genes associated with mobile elements, and was comprised of a majority of genes with unknown or unclear function. Using these genomes, we showed that AGEnt performed well compared to other publically available programs designed to detect accessory genomic elements. We then demonstrated the utility of the AGEnt program by applying it to the draft genomes of two previously unsequenced P. aeruginosa strains, PA99 and PA103.

Conclusions

The P. aeruginosa genome is rich in accessory genetic material. The AGEnt program accurately identified the accessory genomes of newly sequenced P. aeruginosa strains, even when draft genomes were used. As P. aeruginosa genomes become available at an increasingly rapid pace, this program will be useful in cataloging the expanding accessory genome of this bacterium and in discerning correlations between phenotype and accessory genome makeup. The combination of Spine and AGEnt should be useful in defining the accessory genomes of other bacterial species as well.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-737) contains supplementary material, which is available to authorized users.     

Comparative genomic analyses of seventeen Streptococcus pneumoniae strains: insights into the pneumococcal supragenome

The distributed-genome hypothesis (DGH) states that pathogenic bacteria possess a supragenome that is much larger than the genome of any single bacterium and that these pathogens utilize genetic recombination and a large, noncore set of genes as a means of diversity generation. We sequenced the genomes of eight nasopharyngeal strains of Streptococcus pneumoniae isolated from pediatric patients with upper respiratory symptoms and performed quantitative genomic analyses among these and nine publicly available pneumococcal strains. Coding sequences from all strains were grouped into 3,170 orthologous gene clusters, of which 1,454 (46%) were conserved among all 17 strains. The majority of the gene clusters, 1,716 (54%), were not found in all strains. Genic differences per strain pair ranged from 35 to 629 orthologous clusters, with each strain's genome containing between 21 and 32% noncore genes. The distribution of the orthologous clusters per genome for the 17 strains was entered into the finite-supragenome model, which predicted that (i) the S. pneumoniae supragenome contains more than 5,000 orthologous clusters and (ii) 99% of the orthologous clusters (~3,000) that are represented in the S. pneumoniae population at frequencies of ≥0.1 can be identified if 33 representative genomes are sequenced. These extensive genic diversity data support the DGH and provide a basis for understanding the great differences in clinical phenotype associated with various pneumococcal strains. When these findings are taken together with previous studies that demonstrated the presence of a supragenome for Streptococcus agalactiae and Haemophilus influenzae, it appears that the possession of a distributed genome is a common host interaction strategy.     

Intraspecies comparison of Streptomyces pratensis genomes reveals high levels of recombination and gene conservation between strains of disparate geographic origin

Background

Streptomyces are widespread bacteria that contribute to the terrestrial carbon cycle and produce the majority of clinically useful antibiotics. While interspecific genomic diversity has been investigated among Streptomyces, information is lacking on intraspecific genomic diversity. Streptomyces pratensis has high rates of homologous recombination but the impact of such gene exchange on genome evolution and the evolution of natural product gene clusters remains uncharacterized.

Results

We report draft genome sequences of four S. pratensis strains and compare to the complete genome of Streptomyces flavogriseus IAF-45-CD (=ATCC 33331), a strain recently reclassified to S. pratensis. Despite disparate geographic origins, the genomes are highly similar with 85.9% of genes present in the core genome and conservation of all natural product gene clusters. Natural products include a novel combination of carbapenem and beta-lactamase inhibitor gene clusters. While high intraspecies recombination rates abolish the phylogenetic signal across the genome, intraspecies recombination is suppressed in two genomic regions. The first region is centered on an insertion/deletion polymorphism and the second on a hybrid NRPS-PKS gene. Finally, two gene families accounted for over 25% of the divergent genes in the core genome. The first includes homologs of bldB (required for spore development and antibiotic production) while the second includes homologs of an uncharacterized protein with a helix-turn-helix motif (hpb). Genes from these families co-occur with fifteen pairs spread across the genome. These genes have evidence for co-evolution of co-localized pairs, supporting previous assertions that these genes may function akin to a toxin-antitoxin system.

Conclusions

S. pratensis genomes are highly similar with exceptional levels of recombination which erase phylogenetic signal among strains of the species. This species has a large core genome and variable terminal regions that are smaller than those found in interspecies comparisons. There is no geographic differentiation between these strains, but there is evidence for local linkage disequilibrium affecting two genomic regions. We have also shown further observational evidence that the DUF397-HTH (bldB and hpb) are a novel toxin-antitoxin pair.     

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.     

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.     

The Pseudomonas aeruginosa Pathogenicity Island PAPI-1 Is Transferred via a Novel Type IV Pilus

Pseudomonas aeruginosa is a major cause of nosocomial infections, particularly in immunocompromised patients or in individuals with cystic fibrosis. The notable ability of P. aeruginosa to inhabit a broad range of environments, including humans, is in part due to its large and diverse genomic repertoire. The genomes of most strains contain a significant number of large and small genomic islands, including those carrying virulence determinants (pathogenicity islands). The pathogenicity island PAPI-1 of strain PA14 is a cluster of 115 genes, and some have been shown to be responsible for virulence phenotypes in a number of infection models. We have previously demonstrated that PAPI-1 can be transferred to other P. aeruginosa strains following excision from the chromosome of the donor. Here we show that PAPI-1 is transferred into recipient P. aeruginosa by a conjugative mechanism, via a type IV pilus, encoded in PAPI-1 by a 10-gene cluster which is closely related to the genes in the enterobacterial plasmid R64. We also demonstrate that the precursor of the major pilus subunit, PilS2, is processed by the chromosomally encoded prepillin peptidase PilD but not its paralog FppA. Our results suggest that the pathogenicity island PAPI-1 may have evolved by acquisition of a conjugation system but that because of its dependence on an essential chromosomal determinant, its transfer is restricted to P. aeruginosa or other species capable of providing a functional prepilin peptidase.The genomes of a number of microorganisms, primarily those that have a capability of changing and adapting to a wide range of environments, evolve by acquisition of novel genetic information in blocks of genes via a process referred to as horizontal gene transfer (HGT). Other bacterial species change their genetic repertoire minimally, principally those that have adapted to a particular environment and, in the case of pathogenic bacteria, to a specific host. For HGT-mediated acquisition of genes to occur, a recipient has to be in an environment where donor genetic material is available, such as different strains of the same species cohabitating a shared niche or growing in a large and diverse community of several hundred different microorganisms. Moreover, for bacteria to become successful recipients of foreign genetic material, they have to posses one of three mechanisms of HGT: natural competence for uptake of foreign DNA (transformation), the ability to be infected by transducing bacteriophages (transduction), or serving as recipients during conjugation of plasmids or mobilized chromosomal DNA (conjugation). Acquired genetic material can consist of individual genes, where they recombine into homologous sequences in the recipient genome and thus increase the genetic diversity. However, large blocks of hundreds of contiguous genes in elements called genomic islands can be also transferred between bacteria, allowing the recipient microorganisms to acquire a number of new traits by a single HGT event.Previous studies comparing genomes of the opportunistic pathogen Pseudomonas aeruginosa pointed toward HGT as an important factor in its evolution (23). The genomes of all strains sequenced to date contain a significant fraction of horizontally acquired genes, in genomic islands and prophages, consisting of a few to several hundred. These islands can be recognized by the presence of certain signature features, such as an atypical nucleotide composition relative to the rest of the genome, location within predicted sites of chromosomal integration (att sites), and the presence of genes encoding bacteriophages and conjugation machineries. We have recently demonstrated that PAPI-1, a large P. aeruginosa genomic (pathogenicity) island, can be excised from its tRNA att site and that a copy can be transferred into a recipient, where it integrates into the same tRNA gene (27). Inspection of the genes in PAPI-1 and features of the transfer process, namely, an integrase-dependent excision and formation of a circular intermediate, suggested that PAPI-1 is an integrative and conjugative element and that it is likely transferred by a conjugative mechanism.Here we extended our analysis of PAPI-1 by testing its transfer from a preselected group of P. aeruginosa PA14 mutants with insertions in each of the genes on the island. Among those mutants that were defective in PAPI-1 transfer, one group of genes encode homologs of type IV pilus proteins. While type IV pili have been found to be involved primarily in bacterial adhesion and twitching motility (24), the PAPI-1-encoded pilus is closely related to the conjugative apparatus of plasmid R64 (14). Moreover, we show that an essential posttranslational modification reaction, converting the precursor of the major pilin subunit encoded in PAPI-1 into a mature protein, is carried out by an enzyme encoded in the chromosome of the donor cells. The acquisition and adaptation of groups of genes and subsequent loss of an essential function may represent a novel evolutionary strategy, limiting horizontal transfer to a specific bacterial species.     

Identification and analysis of whole microcystin synthetase genes from two Korean strains of the cyanobacterium Microcystis aeruginosa

Microcystins are cyanobacterial hepatotoxins, and are produced by nonribosomal enzyme complexes, mcy gene cluster. In this study, we report on whole mcy gene clusters from two Korean strains of M. aeruginosa that were blooming in Lake Paldang (FCY-26) and Geum river (FCY-28). Their specific gene locus, amino acid information, and sub-cluster orientation were also characterized in both strains. Both gene clusters are of 55 kb, and also each length, number and the arrangement are identical. Their sequence analysis revealed a cluster of 10 genes (mcyA, B, C, D, E, F, G, H, I, and J) involved in the biosynthesis of microcystin, and mcyABC and mcyDEFGHIJ formed two polycistronic operon structures that are transcribed bidirectionally from a central promoter region between mcyA and mcyD. The analysis of SNPs provided different nucleotide composition and amino acid variations in two Korean strains of M. aeruginosa. This approach is useful to develop genetic indicators identifying toxic cyanobacteria and their cyanotoxins, and helpful for a better understanding of the diversities of mcy gene clusters, the biosynthesis of microcystin, and the mediation of environmental parameters causing algal blooming and HABs.     

Core and accessory genome architecture in a group of Pseudomonas aeruginosa Mu-like phages

Background

Bacteriophages that infect the opportunistic pathogen Pseudomonas aeruginosa have been classified into several groups. One of them, which includes temperate phage particles with icosahedral heads and long flexible tails, bears genomes whose architecture and replication mechanism, but not their nucleotide sequences, are like those of coliphage Mu. By comparing the genomic sequences of this group of P. aeruginosa phages one could draw conclusions about their ontogeny and evolution.

Results

Two newly isolated Mu-like phages of P. aeruginosa are described and their genomes sequenced and compared with those available in the public data banks. The genome sequences of the two phages are similar to each other and to those of a group of P. aeruginosa transposable phages. Comparing twelve of these genomes revealed a common genomic architecture in the group. Each phage genome had numerous genes with homologues in all the other genomes and a set of variable genes specific for each genome. The first group, which comprised most of the genes with assigned functions, was named “core genome”, and the second group, containing mostly short ORFs without assigned functions was called “accessory genome”. Like in other phage groups, variable genes are confined to specific regions in the genome.

Conclusion

Based on the known and inferred functions for some of the variable genes of the phages analyzed here, they appear to confer selective advantages for the phage survival under particular host conditions. We speculate that phages have developed a mechanism for horizontally acquiring genes to incorporate them at specific loci in the genome that help phage adaptation to the selective pressures imposed by the host.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1146) contains supplementary material, which is available to authorized users.     

Multilocus Sequence Typing and Phylogenetic Analyses of Pseudomonas aeruginosa Isolates from the Ocean

Recent isolation of Pseudomonas aeruginosa strains from the open ocean and subsequent pulsed-field gel electrophoresis analyses indicate that these strains have a unique genotype (N. H. Khan, Y. Ishii, N. Kimata-Kino, H. Esaki, T. Nishino, M. Nishimura, and K. Kogure, Microb. Ecol. 53:173-186, 2007). We hypothesized that ocean P. aeruginosa strains have a unique phylogenetic position relative to other strains. The objective of this study was to clarify the intraspecies phylogenetic relationship between marine strains and other strains from various geographical locations. Considering the advantages of using databases, multilocus sequence typing (MLST) was chosen for the typing and discrimination of ocean P. aeruginosa strains. Seven housekeeping genes (acsA, aroE, guaA, mutL, nuoD, ppsA, and trpE) were analyzed, and the results were compared with data on the MLST website. These genes were also used for phylogenetic analysis of P. aeruginosa. Rooted and unrooted phylogenetic trees were generated for each gene locus and the concatenated gene fragments. MLST data showed that all the ocean strains were new. Trees constructed for individual and concatenated genes revealed that ocean P. aeruginosa strains have clusters distinct from those of other P. aeruginosa strains. These clusters roughly reflected the geographical locations of the isolates. These data support our previous findings that P. aeruginosa strains are present in the ocean. It can be concluded that the ocean P. aeruginosa strains have diverged from other isolates and form a distinct cluster based on MLST and phylogenetic analyses of seven housekeeping genes.     

Genomic sequencing of different sequevars of Ralstonia solanacearum belonging to the Moko ecotype

Banana vascular wilt or Moko is a disease caused by Ralstonia solanacearum. This study aimed to sequence, assemble, annotate, and compare the genomes of R. solanacearum Moko ecotypes of different sequevar strains from Brazil. Average nucleotide identity analyses demonstrated a high correlation (> 96%) between the genome sequences of strains CCRMRs277 (sequevar IIA-24), CCRMRs287 (IIB-4), CCRMRs304 (IIA-24), and CCRMRsB7 (IIB-25), which were grouped into phylotypes IIA and IIB. The number of coding sequences present in chromosomes and megaplasmids varied from 3,070 to 3,521 and 1,669 to 1,750, respectively. Pangenome analysis identified 3,378 clusters in the chromosomes, of which 2,604 were shared by all four analyzed genomes and 2,580 were single copies. In megaplasmids, 1,834 clusters were identified, of which 1,005 were shared by all four genomes and 992 were identified as single copies. Strains CCRMRsB7 and CCRMRs287 differed from the others by having unique clusters in both their chromosomes and megaplasmids, and CCRMRsB7 possessed the largest genome among all Moko ecotype strains sequenced to date. Therefore, the genomic information obtained in this study provides a theoretical basis for the identification, characterization, and phylogenetic analysis of R. solanacearum Moko ecotypes.     

Intraspecific variability in response to phosphorus depleted conditions in the cyanobacteria Microcystis aeruginosa and Raphidiopsis raciborskii

Phosphorus loading plays an important role in the occurrence of cyanobacterial blooms and understanding how this nutrient affects the physiology of cyanobacteria is imperative to manage these phenomena. Microcystis aeruginosa and Raphidiopsis raciborskii are cyanobacterial species that form potentially toxic blooms in freshwater ecosystems worldwide. Blooms comprise numerous strains with high trait variability, which can contribute to the widespread distribution of these species. Here, we explored the intraspecific variability in response to phosphorus depleted conditions (P-) testing five strains of each species. Strains could be differentiated by cell volume or genetic profiles except for those of the same species, sampling location and date, though these presented differences in their response to (P-). Although differently affected by (P-) over 10 days, all strains were able to grow and maintain photosynthetic activity. For most M. aeruginosa and R. raciborskii strains growth rates were not significantly different comparing (P+) and (P-) conditions. After ten days in (P-), only one M. aeruginosa strain and two R. raciborskii strains showed reduction in biovolume yield as compared to (P+) but in most strains chlorophyll-a concentrations were lower in (P-) than in (P+). Reduced photosystem II efficiency was found for only one R. raciborskii strain while all M. aeruginosa strains were affected. Only two M. aeruginosa and one R. raciborskii strain increased alkaline phosphatase activity under (P-) as compared to (P+). Variation in P-uptake was also observed but comparison among strains yielded homogeneous groups comprised of representatives of both species. Comparing the response of each species as a whole, the (P-) condition affected growth rate, biovolume yield and chlorophyll yield. However, these parameters revealed variation among strains of the same species to the extent that differences between M. aeruginosa and R. raciborskii were not significant. Taken together, these results do not support the idea that R. raciborskii, as a species, can withstand phosphorus limitation better than M. aeruginosa and also point that the level of intraspecific variation may preclude generalizations based on studies that use only one or few strains.     

What It Takes to Be a Pseudomonas aeruginosa? The Core Genome of the Opportunistic Pathogen Updated

Pseudomonas aeruginosa is an opportunistic bacterial pathogen able to thrive in highly diverse ecological niches and to infect compromised patients. Its genome exhibits a mosaic structure composed of a core genome into which accessory genes are inserted en bloc at specific sites. The size and the content of the core genome are open for debate as their estimation depends on the set of genomes considered and the pipeline of gene detection and clustering. Here, we redefined the size and the content of the core genome of P. aeruginosa from fully re-analyzed genomes of 17 reference strains. After the optimization of gene detection and clustering parameters, the core genome was defined at 5,233 orthologs, which represented ~ 88% of the average genome. Extrapolation indicated that our panel was suitable to estimate the core genome that will remain constant even if new genomes are added. The core genome contained resistance determinants to the major antibiotic families as well as most metabolic, respiratory, and virulence genes. Although some virulence genes were accessory, they often related to conserved biological functions. Long-standing prophage elements were subjected to a genetic drift to eventually display a G+C content as higher as that of the core genome. This contrasts with the low G+C content of highly conserved ribosomal genes. The conservation of metabolic and respiratory genes could guarantee the ability of the species to thrive on a variety of carbon sources for energy in aerobiosis and anaerobiosis. Virtually all the strains, of environmental or clinical origin, have the complete toolkit to become resistant to the major antipseudomonal compounds and possess basic pathogenic mechanisms to infect humans. The knowledge of the genes shared by the majority of the P. aeruginosa isolates is a prerequisite for designing effective therapeutics to combat the wide variety of human infections.     

Comparative Genomic Analysis of Sulfurospirillum cavolei MES Reconstructed from the Metagenome of an Electrosynthetic Microbiome

Sulfurospirillum spp. play an important role in sulfur and nitrogen cycling, and contain metabolic versatility that enables reduction of a wide range of electron acceptors, including thiosulfate, tetrathionate, polysulfide, nitrate, and nitrite. Here we describe the assembly of a Sulfurospirillum genome obtained from the metagenome of an electrosynthetic microbiome. The ubiquity and persistence of this organism in microbial electrosynthesis systems suggest it plays an important role in reactor stability and performance. Understanding why this organism is present and elucidating its genetic repertoire provide a genomic and ecological foundation for future studies where Sulfurospirillum are found, especially in electrode-associated communities. Metabolic comparisons and in-depth analysis of unique genes revealed potential ecological niche-specific capabilities within the Sulfurospirillum genus. The functional similarities common to all genomes, i.e., core genome, and unique gene clusters found only in a single genome were identified. Based upon 16S rRNA gene phylogenetic analysis and average nucleotide identity, the Sulfurospirillum draft genome was found to be most closely related to Sulfurospirillum cavolei. Characterization of the draft genome described herein provides pathway-specific details of the metabolic significance of the newly described Sulfurospirillum cavolei MES and, importantly, yields insight to the ecology of the genus as a whole. Comparison of eleven sequenced Sulfurospirillum genomes revealed a total of 6246 gene clusters in the pan-genome. Of the total gene clusters, 18.5% were shared among all eleven genomes and 50% were unique to a single genome. While most Sulfurospirillum spp. reduce nitrate to ammonium, five of the eleven Sulfurospirillum strains encode for a nitrous oxide reductase (nos) cluster with an atypical nitrous-oxide reductase, suggesting a utility for this genus in reduction of the nitrous oxide, and as a potential sink for this potent greenhouse gas.     

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.     

Diversity of Pseudomonas Genomes,Including Populus-Associated Isolates,as Revealed by Comparative Genome Analysis

The Pseudomonas genus contains a metabolically versatile group of organisms that are known to occupy numerous ecological niches, including the rhizosphere and endosphere of many plants. Their diversity influences the phylogenetic diversity and heterogeneity of these communities. On the basis of average amino acid identity, comparative genome analysis of >1,000 Pseudomonas genomes, including 21 Pseudomonas strains isolated from the roots of native Populus deltoides (eastern cottonwood) trees resulted in consistent and robust genomic clusters with phylogenetic homogeneity. All Pseudomonas aeruginosa genomes clustered together, and these were clearly distinct from other Pseudomonas species groups on the basis of pangenome and core genome analyses. In contrast, the genomes of Pseudomonas fluorescens were organized into 20 distinct genomic clusters, representing enormous diversity and heterogeneity. Most of our 21 Populus-associated isolates formed three distinct subgroups within the major P. fluorescens group, supported by pathway profile analysis, while two isolates were more closely related to Pseudomonas chlororaphis and Pseudomonas putida. Genes specific to Populus-associated subgroups were identified. Genes specific to subgroup 1 include several sensory systems that act in two-component signal transduction, a TonB-dependent receptor, and a phosphorelay sensor. Genes specific to subgroup 2 contain hypothetical genes, and genes specific to subgroup 3 were annotated with hydrolase activity. This study justifies the need to sequence multiple isolates, especially from P. fluorescens, which displays the most genetic variation, in order to study functional capabilities from a pangenomic perspective. This information will prove useful when choosing Pseudomonas strains for use to promote growth and increase disease resistance in plants.     

Nonclinical and Clinical Enterococcus faecium Strains,but Not Enterococcus faecalis Strains,Have Distinct Structural and Functional Genomic Features

Certain strains of Enterococcus faecium and Enterococcus faecalis contribute beneficially to animal health and food production, while others are associated with nosocomial infections. To determine whether there are structural and functional genomic features that are distinct between nonclinical (NC) and clinical (CL) strains of those species, we analyzed the genomes of 31 E. faecium and 38 E. faecalis strains. Hierarchical clustering of 7,017 orthologs found in the E. faecium pangenome revealed that NC strains clustered into two clades and are distinct from CL strains. NC E. faecium genomes are significantly smaller than CL genomes, and this difference was partly explained by significantly fewer mobile genetic elements (ME), virulence factors (VF), and antibiotic resistance (AR) genes. E. faecium ortholog comparisons identified 68 and 153 genes that are enriched for NC and CL strains, respectively. Proximity analysis showed that CL-enriched loci, and not NC-enriched loci, are more frequently colocalized on the genome with ME. In CL genomes, AR genes are also colocalized with ME, and VF are more frequently associated with CL-enriched loci. Genes in 23 functional groups are also differentially enriched between NC and CL E. faecium genomes. In contrast, differences were not observed between NC and CL E. faecalis genomes despite their having larger genomes than E. faecium. Our findings show that unlike E. faecalis, NC and CL E. faecium strains are equipped with distinct structural and functional genomic features indicative of adaptation to different environments.     

Transcriptional profiling of Actinobacillus pleuropneumoniae during the acute phase of a natural infection in pigs

Background

The marine cyanobacterium Prochlorococcus marinus, having multiple ecotypes of distinct genotypic/phenotypic traits and being the first documented example of genome shrinkage in free-living organisms, offers an ideal system for studying niche-driven molecular micro-diversity in closely related microbes. The present study, through an extensive comparative analysis of various genomic/proteomic features of 6 high light (HL) and 6 low light (LL) adapted strains, makes an attempt to identify molecular determinants associated with their vertical niche partitioning.

Results

Pronounced strand-specific asymmetry in synonymous codon usage is observed exclusively in LL strains. Distinct dinucleotide abundance profiles are exhibited by 2 LL strains with larger genomes and G+C-content ≈ 50% (group LLa), 4 LL strains having reduced genomes and G+C-content ≈ 35-37% (group LLb), and 6 HL strains. Taking into account the emergence of LLa, LLb and HL strains (based on 16S rRNA phylogeny), a gradual increase in average aromaticity, pI values and beta- & coil-forming propensities and a decrease in mean hydrophobicity, instability indices and helix-forming propensities of core proteins are observed. Greater variations in orthologous gene repertoire are found between LLa and LLb strains, while higher number of positively selected genes exist between LL and HL strains.

Conclusion

Strains of different Prochlorococcus groups are characterized by distinct compositional, physicochemical and structural traits that are not mere remnants of a continuous genetic drift, but are potential outcomes of a grand scheme of niche-oriented stepwise diversification, that might have driven them chronologically towards greater stability/fidelity and invoked upon them a special ability to inhabit diverse oceanic environments.     

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.     

Comparison of 61 Sequenced Escherichia coli Genomes

Escherichia coli is an important component of the biosphere and is an ideal model for studies of processes involved in bacterial genome evolution. Sixty-one publically available E. coli and Shigella spp. sequenced genomes are compared, using basic methods to produce phylogenetic and proteomics trees, and to identify the pan- and core genomes of this set of sequenced strains. A hierarchical clustering of variable genes allowed clear separation of the strains into clusters, including known pathotypes; clinically relevant serotypes can also be resolved in this way. In contrast, when in silico MLST was performed, many of the various strains appear jumbled and less well resolved. The predicted pan-genome comprises 15,741 gene families, and only 993 (6%) of the families are represented in every genome, comprising the core genome. The variable or ‘accessory’ genes thus make up more than 90% of the pan-genome and about 80% of a typical genome; some of these variable genes tend to be co-localized on genomic islands. The diversity within the species E. coli, and the overlap in gene content between this and related species, suggests a continuum rather than sharp species borders in this group of Enterobacteriaceae.