Recombination and population structure in Salmonella enterica

Chlamydiaceae Genomics Reveals Interspecies Admixture and the Recent Evolution of Chlamydia abortus Infecting Lower Mammalian Species and Humans

Chlamydiaceae are obligate intracellular bacteria that cause a diversity of severe infections among humans and livestock on a global scale. Identification of new species since 1989 and emergence of zoonotic infections, including abortion in women, underscore the need for genome sequencing of multiple strains of each species to advance our knowledge of evolutionary dynamics across Chlamydiaceae. Here, we genome sequenced isolates from avian, lower mammalian and human hosts. Based on core gene phylogeny, five isolates previously classified as Chlamydia abortus were identified as members of Chlamydia psittaci and Chlamydia pecorum. Chlamydia abortus is the most recently emerged species and is a highly monomorphic group that lacks the conserved virulence-associated plasmid. Low-level recombination and evidence for adaptation to the placenta echo evolutionary processes seen in recently emerged, highly virulent niche-restricted pathogens, such as Bacillus anthracis. In contrast, gene flow occurred within C. psittaci and other Chlamydiaceae species. The C. psittaci strain RTH, isolated from a red-tailed hawk (Buteo jamaicensis), is an outlying strain with admixture of C. abortus, C. psittaci, and its own population markers. An average nucleotide identity of less than 94% compared with other Chlamydiaceae species suggests that RTH belongs to a new species intermediary between C. psittaci and C. abortus. Hawks, as scavengers and predators, have extensive opportunities to acquire multiple species in their intestinal tract. This could facilitate transformation and homologous recombination with the potential for new species emergence. Our findings indicate that incubator hosts such as birds-of-prey likely promote Chlamydiaceae evolution resulting in novel pathogenic lineages.

Formal Comment to Pettengill: The Time to Most Recent Common Ancestor Does Not (Usually) Approximate the Date of Divergence

In 2013 Zhou et al. concluded that Salmonella enterica serovar Agona represents a genetically monomorphic lineage of recent ancestry, whose most recent common ancestor existed in 1932, or earlier. The Abstract stated 'Agona consists of three lineages with minimal mutational diversity: only 846 single nucleotide polymorphisms (SNPs) have accumulated in the non-repetitive, core genome since Agona evolved in 1932 and subsequently underwent a major population expansion in the 1960s.' These conclusions have now been criticized by Pettengill, who claims that the evolutionary models used to date Agona may not have been appropriate, the dating estimates were inaccurate, and the age of emergence of Agona should have been qualified by an upper limit reflecting the date of its divergence from an outgroup, serovar Soerenga. We dispute these claims. Firstly, Pettengill's analysis of Agona is not justifiable on technical grounds. Secondly, an upper limit for divergence from an outgroup would only be meaningful if the outgroup were closely related to Agona, but close relatives of Agona are yet to be identified. Thirdly, it is not possible to reliably date the time of divergence between Agona and Soerenga. We conclude that Pettengill's criticism is comparable to a tempest in a teapot.

Characterization of the cytolethal distending toxin (typhoid toxin) in non-typhoidal Salmonella serovars

BACKGROUND

For many putative Salmonella enterica subsp. enterica virulence genes, functional characterization across serovars has been limited. Cytolethal distending toxin B (CdtB) is an incompletely characterized virulence factor that is found not only in Salmonella enterica subsp. enterica serovar Typhi (Salmonella Typhi) and dozens of Gram negative bacterial pathogens, but also in non-typhoidal Salmonella (NTS) serovars.

METHODS

A comparative genomics approach was performed to characterize sequence conservation of the typhoid toxin (TT), encoded in the CdtB-islet, between Salmonella Typhi and NTS serovars. The cytotoxic activity of representative Salmonella enterica subsp. enterica serovars Javiana, Montevideo and Schwarzengrund strains and their respective isogenic cdtB mutants was determined in human intestinal epithelial Henle-407 cells by assessment of cell cycle progression of infected cells using fluorescence-activated cell sorting (FACS). Two-way analysis of variance (ANOVA) was used to determine whether cdtB deletion had a significant (p < 0.05) effect on the percentage of Henle-407 cells at each stage of the cell cycle.

RESULTS

Here we show that a CdtB-islet encoding the cytolethal distending toxin B (CdtB), pertussis-like toxin A (PltA), and pertussis-like toxin B (PltB) is present in a dozen NTS serovars and that these proteins have a high level of sequence conservation and each form monophyletic clades with corresponding Salmonella Typhi genes. Human epithelial Henle-407 cells infected with three representative CdtB-encoding NTS serovars displayed G2/M phase cell cycle arrest that was absent in cells infected with corresponding isogenic cdtB null mutants (p < 0.0001 for the factor ∆cdtB deletion).

CONCLUSION

Our results show that CdtB encoded by NTS serovars has a genomic organization, amino acid sequence conservation and biological activity similar to the TT, and thus, may contribute to disease pathogenesis.

The Harvest suite for rapid core-genome alignment and visualization of thousands of intraspecific microbial genomes

Whole-genome sequences are now available for many microbial species and clades, however existing whole-genome alignment methods are limited in their ability to perform sequence comparisons of multiple sequences simultaneously. Here we present the Harvest suite of core-genome alignment and visualization tools for the rapid and simultaneous analysis of thousands of intraspecific microbial strains. Harvest includes Parsnp, a fast core-genome multi-aligner, and Gingr, a dynamic visual platform. Together they provide interactive core-genome alignments, variant calls, recombination detection, and phylogenetic trees. Using simulated and real data we demonstrate that our approach exhibits unrivaled speed while maintaining the accuracy of existing methods. The Harvest suite is open-source and freely available from: http://github.com/marbl/harvest.

Inter-individual differences in the gene content of human gut bacterial species

BACKGROUND

Gene content differences in human gut microbes can lead to inter-individual phenotypic variations such as digestive capacity. It is unclear whether gene content variation is caused by differences in microbial species composition or by the presence of different strains of the same species; the extent of gene content variation in the latter is unknown. Unlike pan-genome studies of cultivable strains, the use of metagenomic data can provide an unbiased view of structural variation of gut bacterial strains by measuring them in their natural habitats, the gut of each individual in this case, representing native boundaries between gut bacterial populations. We analyzed publicly available metagenomic data from fecal samples to characterize inter-individual variation in gut bacterial species.

RESULTS

A comparison of 11 abundant gut bacterial species showed that the gene content of strains from the same species differed, on average, by 13% between individuals. This number is based on gene deletions only and represents a lower limit, yet the variation is already in a similar range as observed between completely sequenced strains of cultivable species. We show that accessory genes that differ considerably between individuals can encode important functions, such as polysaccharide utilization and capsular polysaccharide synthesis loci.

CONCLUSION

Metagenomics can yield insights into gene content variation of strains in complex communities, which cannot be predicted by phylogenetic marker genes alone. The large degree of inter-individual variability in gene content implies that strain resolution must be considered in order to fully assess the functional potential of an individual's human gut microbiome.

Progressive genomic convergence of two Helicobacter pylori strains during mixed infection of a patient with chronic gastritis

OBJECTIVE

To study the detailed nature of genomic microevolution during mixed infection with multiple Helicobacter pylori strains in an individual.

DESIGN

We sampled 18 isolates from a single biopsy from a patient with chronic gastritis and nephritis. Whole-genome sequencing was applied to these isolates, and statistical genetic tools were used to investigate their evolutionary history.

RESULTS

The genomes fall into two clades, reflecting colonisation of the stomach by two distinct strains, and these lineages have accumulated diversity during an estimated 2.8 and 4.2 years of evolution. We detected about 150 clear recombination events between the two clades. Recombination between the lineages is a continuous ongoing process and was detected on both clades, but the effect of recombination in one clade was nearly an order of magnitude higher than in the other. Imputed ancestral sequences also showed evidence of recombination between the two strains prior to their diversification, and we estimate that they have both been infecting the same host for at least 12 years. Recombination tracts between the lineages were, on average, 895 bp in length, and showed evidence for the interspersion of recipient sequences that has been observed in in vitro experiments. The complex evolutionary history of a phage-related protein provided evidence for frequent reinfection of both clades by a single phage lineage during the past 4 years.

CONCLUSIONS

Whole genome sequencing can be used to make detailed conclusions about the mechanisms of genetic change of H. pylori based on sampling bacteria from a single gastric biopsy.

Epidemic Clones, Oceanic Gene Pools, and Eco-LD in the Free Living Marine Pathogen Vibrio parahaemolyticus

We investigated global patterns of variation in 157 whole-genome sequences of Vibrio parahaemolyticus, a free-living and seafood associated marine bacterium. Pandemic clones, responsible for recent outbreaks of gastroenteritis in humans, have spread globally. However, there are oceanic gene pools, one located in the oceans surrounding Asia and another in the Mexican Gulf. Frequent recombination means that most isolates have acquired the genetic profile of their current location. We investigated the genetic structure in the Asian gene pool by calculating the effective population size in two different ways. Under standard neutral models, the two estimates should give similar answers but we found a 27-fold difference. We propose that this discrepancy is caused by the subdivision of the species into a hundred or more ecotypes which are maintained stably in the population. To investigate the genetic factors involved, we used 51 unrelated isolates to conduct a genome-wide scan for epistatically interacting loci. We found a single example of strong epistasis between distant genome regions. A majority of strains had a type VI secretion system associated with bacterial killing. The remaining strains had genes associated with biofilm formation and regulated by cyclic dimeric GMP signaling. All strains had one or other of the two systems and none of isolate had complete complements of both systems, although several strains had remnants. Further "top down" analysis of patterns of linkage disequilibrium within frequently recombining species will allow a detailed understanding of how selection acts to structure the pattern of variation within natural bacterial populations.

Molecular characteristics of Salmonella enterica Paratyphi A in Yunnan Province, southwest China

Previously, the prevalence of Salmonella enterica Paratyphi A in Yunnan was high; and recently Yunnan was the predominant endemic province in China. To identify the molecular epidemiology, antibiotic resistance profile and genotypic diversity of the S. Paratyphi A isolates from 1995 to 2013 in Yunnan, we performed the study. Antibiotic susceptibility tests, pulse-field gel electrophoresis (PFGE) and multilocus sequence typing (MLST) were used to identify the characteristics of the bacterial isolates. The results showed from 1995 to 2013, 366 S. Paratyphi A were isolated: 295 isolates (80.6%) from Yuxi and 68 isolates (18.58%) from Honghe. All of the strains were resistant to nalidixic acid, and some were resistant to ampicillin and trimethoprim/sulfamethoxazole in different years. All the isolates were sensitive to cefotaxime and ciprofloxacin. Identical PFGE with two enzyme digestion patterns were found for 339 isolates. Some environmental isolates in different years were homologous with the strains isolated from food and patients. MLST showed 349 strains were ST85, only 17 isolates were ST129. S. Paratyphi A isolates from Yunnan showed a high similarity, and we found the pathogen isolated from patients, the environment and food had the close epidemiological relationship, forming a transmission circulation. These findings have important implications for paratyphoid-control strategies.

Genome-wide patterns of recombination in the opportunistic human pathogen Pseudomonas aeruginosa

The bacterium Pseudomonas aeruginosa is a significant cause of acute nosocomial infections as well as chronic respiratory infections in patients with cystic fibrosis (CF). Recent reports of the intercontinental spread of a CF-specific epidemic strain, combined with high intrinsic levels of antibiotic resistance, have made this opportunistic pathogen an important public health concern. Strain-specific differences correlate with variation in clinical outcomes of infected CF patients, increasing the urgency to understand the evolutionary origin of genetic factors conferring important phenotypes that enable infection, virulence, or resistance. Here, we describe the genome-wide patterns of homologous and nonhomologous recombination in P. aeruginosa, and the extent to which the genomes are affected by these diversity-generating processes. Based on whole-genome sequence data from 32 clinical isolates of P. aeruginosa, we examined the rate and distribution of recombination along the genome, and its effect on the reconstruction of phylogenetic relationships. Multiple lines of evidence suggested that recombination was common and usually involves short stretches of DNA (200-300 bp). Although mutation was the main source of nucleotide diversity, the import of polymorphisms by homologous recombination contributed nearly as much. We also identified the genomic regions with frequent recombination, and the specific sequences of recombinant origin within epidemic strains. The functional characteristics of the genes contained therein were examined for potential associations with a pathogenic lifestyle or adaptation to the CF lung environment. A common link between many of the high-recombination genes was their functional affiliation with the cell wall, suggesting that the products of recombination may be maintained by selection for variation in cell-surface molecules that allows for evasion of the host immune system.

Population genetic analysis of Streptomyces albidoflavus reveals habitat barriers to homologous recombination in the diversification of streptomycetes

Examining the population structure and the influence of recombination and ecology on microbial populations makes great sense for understanding microbial evolution and speciation. Streptomycetes are a diverse group of bacteria that are widely distributed in nature and a rich source of useful bioactive compounds; however, they are rarely subjected to population genetic investigations. In this study, we applied a five-gene-based multilocus sequence analysis (MLSA) scheme to 41 strains of Streptomyces albidoflavus derived from diverse sources, mainly insects, sea, and soil. Frequent recombination was detected in S. albidoflavus, supported by multiple lines of evidence from the pairwise homoplasy index (Φw) test, phylogenetic discordance, the Shimodaira-Hasegawa (SH) test, and network analysis, underpinning the predominance of homologous recombination within Streptomyces species. A strong habitat signal was also observed in both phylogenetic and Structure 2.3.3 analyses, indicating the importance of ecological difference in shaping the population structure. Moreover, all three habitat-associated groups, particularly the entomic group, demonstrated significantly reduced levels of gene flow with one another, generally revealing habitat barriers to recombination. Therefore, a combined effect of homologous recombination and ecology is inferred for S. albidoflavus, where dynamic evolution is at least partly balanced by the extent that differential distributions of strains among habitats limit genetic exchange. Our study stresses the significance of ecology in microbial speciation and reveals the coexistence of homologous recombination and ecological divergence in the evolution of streptomycetes.

The Harvest suite for rapid core-genome alignment and visualization of thousands of intraspecific microbial genomes

Whole-genome sequences are now available for many microbial species and clades, however existing whole-genome alignment methods are limited in their ability to perform sequence comparisons of multiple sequences simultaneously. Here we present the Harvest suite of core-genome alignment and visualization tools for the rapid and simultaneous analysis of thousands of intraspecific microbial strains. Harvest includes Parsnp, a fast core-genome multi-aligner, and Gingr, a dynamic visual platform. Together they provide interactive core-genome alignments, variant calls, recombination detection, and phylogenetic trees. Using simulated and real data we demonstrate that our approach exhibits unrivaled speed while maintaining the accuracy of existing methods. The Harvest suite is open-source and freely available from: http://github.com/marbl/harvest.

Rapid phylogenetic analysis of large samples of recombinant bacterial whole genome sequences using Gubbins

The emergence of new sequencing technologies has facilitated the use of bacterial whole genome alignments for evolutionary studies and outbreak analyses. These datasets, of increasing size, often include examples of multiple different mechanisms of horizontal sequence transfer resulting in substantial alterations to prokaryotic chromosomes. The impact of these processes demands rapid and flexible approaches able to account for recombination when reconstructing isolates' recent diversification. Gubbins is an iterative algorithm that uses spatial scanning statistics to identify loci containing elevated densities of base substitutions suggestive of horizontal sequence transfer while concurrently constructing a maximum likelihood phylogeny based on the putative point mutations outside these regions of high sequence diversity. Simulations demonstrate the algorithm generates highly accurate reconstructions under realistically parameterized models of bacterial evolution, and achieves convergence in only a few hours on alignments of hundreds of bacterial genome sequences. Gubbins is appropriate for reconstructing the recent evolutionary history of a variety of haploid genotype alignments, as it makes no assumptions about the underlying mechanism of recombination. The software is freely available for download at github.com/sanger-pathogens/Gubbins, implemented in Python and C and supported on Linux and Mac OS X.

Phylogenetic diversity of the enteric pathogen Salmonella enterica subsp. enterica inferred from genome-wide reference-free SNP characters

The enteric pathogen Salmonella enterica is one of the leading causes of foodborne illness in the world. The species is extremely diverse, containing more than 2,500 named serovars that are designated for their unique antigen characters and pathogenicity profiles—some are known to be virulent pathogens, while others are not. Questions regarding the evolution of pathogenicity, significance of antigen characters, diversity of clustered regularly interspaced short palindromic repeat (CRISPR) loci, among others, will remain elusive until a strong evolutionary framework is established. We present the first large-scale S. enterica subsp. enterica phylogeny inferred from a new reference-free k-mer approach of gathering single nucleotide polymorphisms (SNPs) from whole genomes. The phylogeny of 156 isolates representing 78 serovars (102 were newly sequenced) reveals two major lineages, each with many strongly supported sublineages. One of these lineages is the S. Typhi group; well nested within the phylogeny. Lineage-through-time analyses suggest there have been two instances of accelerated rates of diversification within the subspecies. We also found that antigen characters and CRISPR loci reveal different evolutionary patterns than that of the phylogeny, suggesting that a horizontal gene transfer or possibly a shared environmental acquisition might have influenced the present character distribution. Our study also shows the ability to extract reference-free SNPs from a large set of genomes and then to use these SNPs for phylogenetic reconstruction. This automated, annotation-free approach is an important step forward for bacterial disease tracking and in efficiently elucidating the evolutionary history of highly clonal organisms.

Burkholderia pseudomallei sequencing identifies genomic clades with distinct recombination, accessory, and epigenetic profiles

Burkholderia pseudomallei (Bp) is the causative agent of the infectious disease melioidosis. To investigate population diversity, recombination, and horizontal gene transfer in closely related Bp isolates, we performed whole-genome sequencing (WGS) on 106 clinical, animal, and environmental strains from a restricted Asian locale. Whole-genome phylogenies resolved multiple genomic clades of Bp, largely congruent with multilocus sequence typing (MLST). We discovered widespread recombination in the Bp core genome, involving hundreds of regions associated with multiple haplotypes. Highly recombinant regions exhibited functional enrichments that may contribute to virulence. We observed clade-specific patterns of recombination and accessory gene exchange, and provide evidence that this is likely due to ongoing recombination between clade members. Reciprocally, interclade exchanges were rarely observed, suggesting mechanisms restricting gene flow between clades. Interrogation of accessory elements revealed that each clade harbored a distinct complement of restriction-modification (RM) systems, predicted to cause clade-specific patterns of DNA methylation. Using methylome sequencing, we confirmed that representative strains from separate clades indeed exhibit distinct methylation profiles. Finally, using an E. coli system, we demonstrate that Bp RM systems can inhibit uptake of non-self DNA. Our data suggest that RM systems borne on mobile elements, besides preventing foreign DNA invasion, may also contribute to limiting exchanges of genetic material between individuals of the same species. Genomic clades may thus represent functional units of genetic isolation in Bp, modulating intraspecies genetic diversity.

Hidden diversity in honey bee gut symbionts detected by single-cell genomics

Microbial communities in animal guts are composed of diverse, specialized bacterial species, but little is known about how gut bacteria diversify to produce genetically and ecologically distinct entities. The gut microbiota of the honey bee, Apis mellifera, presents a useful model, because it consists of a small number of characteristic bacterial species, each showing signs of diversification. Here, we used single-cell genomics to study the variation within two species of the bee gut microbiota: Gilliamella apicola and Snodgrassella alvi. For both species, our analyses revealed extensive variation in intraspecific divergence of protein-coding genes but uniformly high levels of 16S rRNA similarity. In both species, the divergence of 16S rRNA loci appears to have been curtailed by frequent recombination within populations, while other genomic regions have continuously diverged. Furthermore, gene repertoires differ markedly among strains in both species, implying distinct metabolic capabilities. Our results show that, despite minimal divergence at 16S rRNA genes, in situ diversification occurs within gut communities and generates bacterial lineages with distinct ecological niches. Therefore, important dimensions of microbial diversity are not evident from analyses of 16S rRNA, and single cell genomics has potential to elucidate processes of bacterial diversification.

Gut microbial communities are often complex, consisting of bacteria from divergent phyla as well as multiple strains of each of the constituent species. But because the composition of these communities is typically assessed using 16S rRNA analyses, little is known about genomic changes associated with in situ diversification of bacterial lineages in animal guts. We undertook a single-cell genomic approach to investigate the diversification within two species of the gut microbiota of honey bees. Each species exhibited a surprisingly high level of genomic variation, despite uniformity in the 16S rRNA sequences. Our data indicate that genetically and ecologically distinct lineages can evolve in the gut of the same host species in the presence of frequent recombination at 16S rRNA genes. These findings parallel observations from mammals, suggesting that in situ diversification of a few bacterial lineages is a common pattern in the evolution of gut communities.

To tree or not to tree? Genome-wide quantification of recombination and reticulate evolution during the diversification of strict intracellular bacteria

It is well known that horizontal gene transfer (HGT) is a major force in the evolution of prokaryotes. During the adaptation of a bacterial population to a new ecological niche, and particularly for intracellular bacteria, selective pressures are shifted and ecological niches reduced, resulting in a lower rate of genetic connectivity. HGT and positive selection are therefore two important evolutionary forces in microbial pathogens that drive adaptation to new hosts. In this study, we use genomic distance analyses, phylogenomic networks, tree topology comparisons, and Bayesian inference methods to investigate to what extent HGT has occurred during the evolution of the genus Rickettsia, the effect of the use of different genomic regions in estimating reticulate evolution and HGT events, and the link of these to host range. We show that ecological specialization restricts recombination occurrence in Rickettsia, but other evolutionary processes and genome architecture are also important for the occurrence of HGT. We found that recombination, genomic rearrangements, and genome conservation all show evidence of network-like evolution at whole-genome scale. We show that reticulation occurred mainly, but not only, during the early Rickettsia radiation, and that core proteome genes of every major functional category have experienced reticulated evolution and possibly HGT. Overall, the evolution of Rickettsia bacteria has been tree-like, with evidence of HGT and reticulated evolution for around 10–25% of the core Rickettsia genome. We present evidence of extensive recombination/incomplete lineage sorting (ILS) during the radiation of the genus, probably linked with the emergence of intracellularity in a wide range of hosts.

Phylogenetic analyses on the diversity of Aspergillus fumigatus sensu lato based on five orthologous loci

One hundred isolates of Aspergillus fumigatus sensu lato mainly from China, as well as from Australia, France, India, Indonesia, Ireland, UK, and USA were analyzed to infer their sequence types (STs) and population diversity based on partial calmodulin, calcineurin regulatory subunit B, beta-tubulin, cytochrome C and calcineurin catalytic subunit A genes as well as their mating types, using ClonalFrame, Structure and MEGA software. Our results inferred 48 STs and showed that most of the STs or lineages evolved independently and without clear population structure among them. Whereas one lineage was recognized that could be a true population and in which one clade might diverge into another distinct lineage, namely, a cryptic species, A. neoellipticus. In addition, we found that mutation, parasexual, and sexual recombination could, respectively, play specific roles in the evolution of these fungi. Our results also showed that MAT1-1/MAT1-2 mating type ratios of A. fumigatus sensu lato was biased to nearly 1:1.4 (20/28) when clone-corrected, but when not clone-corrected, the ratio of MAT1-1/MAT1-2 was so biased as near 1:2 (35/65), which might mean that isolates with MAT1-2 are in the process of losing sexual ability preceding those with MAT1-1.

The Dynamic Interactions between Salmonella and the Microbiota, within the Challenging Niche of the Gastrointestinal Tract

Understanding how Salmonella species establish successful infections remains a foremost research priority. This gastrointestinal pathogen not only faces the hostile defenses of the host's immune system, but also faces fierce competition from the large and diverse community of microbiota for space and nutrients. Salmonella have solved these challenges ingeniously. To jump-start growth, Salmonella steal hydrogen produced by the gastrointestinal microbiota. Type 3 effector proteins are subsequently secreted by Salmonella to trigger potent inflammatory responses, which generate the alternative terminal electron acceptors tetrathionate and nitrate. Salmonella exclusively utilize these electron acceptors for anaerobic respiration, permitting metabolic access to abundant substrates such as ethanolamine to power growth blooms. Chemotaxis and flagella-mediated motility enable the identification of nutritionally beneficial niches. The resulting growth blooms also promote horizontal gene transfer amongst the resident microbes. Within the gastrointestinal tract there are opportunities for chemical signaling between host cells, the microbiota, and Salmonella. Host produced catecholamines and bacterial autoinducers form components of this chemical dialogue leading to dynamic interactions. Thus, Salmonella have developed remarkable strategies to initially shield against host defenses and to transiently compete against the intestinal microbiota leading to successful infections. However, the immunocompetent host is subsequently able to reestablish control and clear the infection.

Paratyphoid fever: splicing the global analyses

The incidence of enteric fever caused by Salmonella enterica serovar Paratyphi A (S. Paratyphi A) is increasing in many parts of the world. Although there is no major outbreak of paratyphoid fever in recent years, S. Paratyphi A infection still remains a public health problem in many tropical countries. Therefore, surveillance studies play an important role in monitoring infections and the emergence of multidrug resistance, especially in endemic countries such as India, Nepal, Pakistan and China. In China, enteric fever was caused predominantly by S. Paratyphi A rather than by Salmonella enterica serovar Typhi (S. Typhi). Sometimes, S. Paratyphi A infection can evolve into a carrier state which increases the risk of transmission for travellers. Hence, paratyphoid fever is usually classified as a "travel-associated" disease. To date, diagnosis of paratyphoid fever based on the clinical presentation is not satisfactory as it resembles other febrile illnesses, and could not be distinguished from S. Typhi infection. With the availability of Whole Genome Sequencing technology, the genomes of S. Paratyphi A could be studied in-depth and more specific targets for detection will be revealed. Hence, detection of S. Paratyphi A with Polymerase Chain Reaction (PCR) method appears to be a more reliable approach compared to the Widal test. On the other hand, due to increasing incidence of S. Paratyphi A infections worldwide, the need to produce a paratyphoid vaccine is essential and urgent. Hence various vaccine projects that involve clinical trials have been carried out. Overall, this review provides the insights of S. Paratyphi A, including the bacteriology, epidemiology, management and antibiotic susceptibility, diagnoses and vaccine development.

Efficient inference of recombination hot regions in bacterial genomes

In eukaryotes, detailed surveys of recombination rates have shown variation at multiple genomic scales and the presence of “hotspots” of highly elevated recombination. In bacteria, studies of recombination rate variation are less developed, in part because there are few analysis methods that take into account the clonal context within which bacterial evolution occurs. Here, we focus in particular on identifying “hot regions” of the genome where DNA is transferred frequently between isolates. We present a computationally efficient algorithm based on the recently developed “chromosome painting” algorithm, which characterizes patterns of haplotype sharing across a genome. We compare the average genome wide painting, which principally reflects clonal descent, with the painting for each site which additionally reflects the specific deviations at the site due to recombination. Using simulated data, we show that hot regions have consistently higher deviations from the genome wide average than normal regions. We applied our approach to previously analyzed Escherichia coli genomes and revealed that the new method is highly correlated with the number of recombination events affecting each site inferred by ClonalOrigin, a method that is only applicable to small numbers of genomes. Furthermore, we analyzed recombination hot regions in Campylobacter jejuni by using 200 genomes. We identified three recombination hot regions, which are enriched for genes related to membrane proteins. Our approach and its implementation, which is downloadable from https://github.com/bioprojects/orderedPainting, will help to develop a new phase of population genomic studies of recombination in prokaryotes.

Frequent recombination shapes the epidemic population structure of Planktothrix (Cyanoprokaryota) in Italian subalpine lakes

The planktonic genus Planktothrix, as other cyanobacteria, shows signals of both homologous and nonhomologous recombination. However, the frequency of recombination and its effect on Planktothrix population structuring is unknown. We isolated 290 Planktothrix strains from seven neighboring lakes in the subalpine Italian region and analyzed these using multilocus sequence typing. Four of six loci analyzed were polymorphic, resulting in 20 distinct multilocus genotypes. Association indices among alleles at different loci were suggestive of an "epidemic population structure," resulting from an explosive (and temporary) dominance of one genotype against a panmictic background. ClonalFrame analyses supported this view by detecting: (i) three major clades affected by three distinct recombination events, (ii) a recombination rate about equal to the mutation rate, and (iii) the fact that recombination had an impact on introducing molecular diversity more than double the mutation rate. Furthermore, analysis of molecular variance over an annual cycle in three of seven lakes revealed that both local clonal expansion and recombination processes affected among-lake diversity. Our observations suggest that recombination affects microevolution of Planktothrix and that an epidemic structure can emerge in populations of this genus.

Salmonella pathogenicity and host adaptation in chicken-associated serovars

Enteric pathogens such as Salmonella enterica cause significant morbidity and mortality. S. enterica serovars are a diverse group of pathogens that have evolved to survive in a wide range of environments and across multiple hosts. S. enterica serovars such as S. Typhi, S. Dublin, and S. Gallinarum have a restricted host range, in which they are typically associated with one or a few host species, while S. Enteritidis and S. Typhimurium have broad host ranges. This review examines how S. enterica has evolved through adaptation to different host environments, especially as related to the chicken host, and continues to be an important human pathogen. Several factors impact host range, and these include the acquisition of genes via horizontal gene transfer with plasmids, transposons, and phages, which can potentially expand host range, and the loss of genes or their function, which would reduce the range of hosts that the organism can infect. S. Gallinarum, with a limited host range, has a large number of pseudogenes in its genome compared to broader-host-range serovars. S. enterica serovars such as S. Kentucky and S. Heidelberg also often have plasmids that may help them colonize poultry more efficiently. The ability to colonize different hosts also involves interactions with the host's immune system and commensal organisms that are present. Thus, the factors that impact the ability of Salmonella to colonize a particular host species, such as chickens, are complex and multifactorial, involving the host, the pathogen, and extrinsic pressures. It is the interplay of these factors which leads to the differences in host ranges that we observe today.

Molecular methods for serovar determination of Salmonella

Salmonella is a diverse foodborne pathogen, which has more than 2600 recognized serovars. Classification of Salmonella isolates into serovars is essential for surveillance and epidemiological investigations; however, determination of Salmonella serovars, by traditional serotyping, has some important limitations (e.g. labor intensive, time consuming). To overcome these limitations, multiple methods have been investigated to develop molecular serotyping schemes. Currently, molecular methods to predict Salmonella serovars include (i) molecular subtyping methods (e.g. PFGE, MLST), (ii) classification using serovar-specific genomic markers and (iii) direct methods, which identify genes encoding antigens or biosynthesis of antigens used for serotyping. Here, we reviewed reported methodologies for Salmonella molecular serotyping and determined the "serovar-prediction accuracy", as the percentage of isolates for which the serovar was correctly classified by a given method. Serovar-prediction accuracy ranged from 0 to 100%, 51 to 100% and 33 to 100% for molecular subtyping, serovar-specific genomic markers and direct methods, respectively. Major limitations of available schemes are errors in predicting closely related serovars (e.g. Typhimurium and 4,5,12:i:-), and polyphyletic serovars (e.g. Newport, Saintpaul). The high diversity of Salmonella serovars represents a considerable challenge for molecular serotyping approaches. With the recent improvement in sequencing technologies, full genome sequencing could be developed into a promising molecular approach to serotype Salmonella.

Inference of the properties of the recombination process from whole bacterial genomes

Patterns of linkage disequilibrium, homoplasy, and incompatibility are difficult to interpret because they depend on several factors, including the recombination process and the population structure. Here we introduce a novel model-based framework to infer recombination properties from such summary statistics in bacterial genomes. The underlying model is sequentially Markovian so that data can be simulated very efficiently, and we use approximate Bayesian computation techniques to infer parameters. As this does not require us to calculate the likelihood function, the model can be easily extended to investigate less probed aspects of recombination. In particular, we extend our model to account for the bias in the recombination process whereby closely related bacteria recombine more often with one another. We show that this model provides a good fit to a data set of Bacillus cereus genomes and estimate several recombination properties, including the rate of bias in recombination. All the methods described in this article are implemented in a software package that is freely available for download at http://code.google.com/p/clonalorigin/.

Nonagricultural reservoirs contribute to emergence and evolution of Pseudomonas syringae crop pathogens

While the existence of environmental reservoirs of human pathogens is well established, less is known about the role of nonagricultural environments in emergence, evolution, and spread of crop pathogens. Here, we analyzed phylogeny, virulence genes, host range, and aggressiveness of Pseudomonas syringae strains closely related to the tomato pathogen P. syringae pv. tomato (Pto), including strains isolated from snowpack and streams. The population of Pto relatives in nonagricultural environments was estimated to be large and its diversity to be higher than that of the population of Pto and its relatives on crops. Ancestors of environmental strains, Pto, and other genetically monomorphic crop pathogens were inferred to have frequently recombined, suggesting an epidemic population structure for P. syringae. Some environmental strains have repertoires of type III-secreted effectors very similar to Pto, are almost as aggressive on tomato as Pto, but have a wider host range than typical Pto strains. We conclude that crop pathogens may have evolved through a small number of evolutionary events from a population of less aggressive ancestors with a wider host range present in nonagricultural environments.

Genomic evolution and transmission of Helicobacter pylori in two South African families

Helicobacter pylori infects the stomachs of one in two humans and can cause sequelae that include ulcers and cancer. Here we sequenced the genomes of 97 H. pylori isolates from 52 members of two families living in rural conditions in South Africa. From each of 45 individuals, two H. pylori strains were isolated from the antrum and corpus parts of the stomach, and comparisons of their genomes enabled us to study within-host evolution. In 5 of these 45 hosts, the two genomes were too distantly related to be derived from each other and therefore represented evidence of multiple infections. From the remaining 40 genome pairs, we estimated that the synonymous mutation rate was 1.38 × 10(-5) per site per year, with a low effective population size within host probably reflecting population bottlenecks and immune selection. Some individuals showed very little evidence for recombination, whereas in others, recombination introduced up to 100-times more substitutions than mutation. These differences may reflect unequal opportunities for recombination depending on the presence or absence of multiple infections. Comparing the genomes carried by distinct individuals enabled us to establish probable transmission links. Transmission events were found significantly more frequently between close relatives, and between individuals living in the same house. We found, however, that a majority of individuals (27/52) were not linked by transmission to other individuals. Our results suggest that transmission does not always occur within families, and that coinfection with multiple strains is frequent and evolutionarily important despite a fast turnover of the infecting strains within-host.

Multilocus sequence typing reveals evidence of homologous recombination linked to antibiotic resistance in the genus Salinispora

The three closely related species that currently comprise the genus Salinispora were analyzed using a multilocus sequence typing approach targeting 48 strains derived from four geographic locations. Phylogenetic congruence and a well-supported concatenated tree provide strong support for the delineation of the three species as currently described and the basal relationship of Salinispora arenicola to the more recently diverged sister taxa S. tropica and S. pacifica. The phylogeny of the initial region of the rpoB gene sequenced was atypical, placing the related genera Micromonospora and Verrucosispora within the Salinispora clade. This phylogenetic incongruence was subsequently ascribed to a homologous-recombination event in a portion of the gene associated with resistance to compounds in the rifamycin class, which target RpoB. All S. arenicola strains produced compounds in this class and possessed resistance-conferring amino acid changes in RpoB. The phylogeny of a region of the rpoB gene that is not associated with rifamycin resistance was congruent with the other housekeeping genes. The link between antibiotic resistance and homologous recombination suggests that incongruent phylogenies provide opportunities to identify the molecular targets of secondary metabolites, an observation with potential relevance for drug discovery efforts. Low ratios of interspecies recombination to mutation, even among cooccurring strains, coupled with high levels of within-species recombination suggest that the three species have been described in accordance with natural barriers to recombination.

Differing lifestyles of Staphylococcus epidermidis as revealed through Bayesian clustering of multilocus sequence types

Staphylococcus epidermidis is part of the normal bacterial flora of human skin and a leading cause of infections associated with indwelling medical devices. Previous phylogenetic analyses of subgenomic data have been unable to distinguish between S. epidermidis strains with nosocomial or commensal lifestyles, despite the identification of specific phenotypes and accessory genes that may contribute to such lifestyles. To attempt to better define the population structure of this species, the international S. epidermidis multilocus sequence typing database was analyzed with the Bayesian clustering programs STRUCTURE and BAPS. A total of six genetic clusters (GCs) were identified. A local population of S. epidermidis from clinical specimens was classified according to these six GCs, and further characterized for antibiotic susceptibilities, biofilm, and various genetic markers. GC5 was abundant and significantly enriched for isolates that were resistant to four classes of antibiotics, high biofilm production, and positive for the virulence markers icaA, IS256, and sesD/bhp, indicating its potential clinical relevance. In contrast, GC2 was rare and contained the only isolates positive for the putative commensal marker, fdh. GC1 and GC6 were abundant but not significantly associated with any of the examined characteristics, except for sesF/aap and GC6. GC3 was rare and identified as a potential genetic sink that received, but did not donate, core genetic material from other GCs. In conclusion, population genetics analyses were essential for identifying clusters of strains that may differ in their adaptation to nosocomial or commensal lifestyles. These results provide a new, population genetics framework for studying S. epidermidis.

Population structure and evolution of non-O1/non-O139 Vibrio cholerae by multilocus sequence typing

Pathogenic non-O1/non-O139 Vibrio cholerae strains can cause sporadic outbreaks of cholera worldwide. In this study, multilocus sequence typing (MLST) of seven housekeeping genes was applied to 55 non-O1/non-O139 isolates from clinical and environmental sources. Data from five published O1 isolates and 17 genomes were also included, giving a total of 77 isolates available for analysis. There were 66 sequence types (STs), with the majority being unique, and only three clonal complexes. The V. cholerae strains can be divided into four subpopulations with evidence of recombination among the subpopulations. Subpopulations I and III contained predominantly clinical strains. PCR screening for virulence factors including Vibrio pathogenicity island (VPI), cholera toxin prophage (CTXΦ), type III secretion system (T3SS), and enterotoxin genes (rtxA and sto/stn) showed that combinations of these factors were present in the clinical isolates with 85.7% having rtxA, 51.4% T3SS, 31.4% VPI, 31.4% sto/stn (NAG-ST) and 11.4% CTXΦ. These factors were also present in environmental isolates but at a lower frequency. Five strains previously mis-identified as V. cholerae serogroups O114 to O117 were also analysed and formed a separate population with V. mimicus. The MLST scheme developed in this study provides a framework to identify sporadic cholera isolates by genetic identity.

Chromosome painting in silico in a bacterial species reveals fine population structure

Identifying population structure forms an important basis for genetic and evolutionary studies. Most current methods to identify population structure have limitations in analyzing haplotypes and recombination across the genome. Recently, a method of chromosome painting in silico has been developed to overcome these shortcomings and has been applied to multiple human genome sequences. This method detects the genome-wide transfer of DNA sequence chunks through homologous recombination. Here, we apply it to the frequently recombining bacterial species Helicobacter pylori that has infected Homo sapiens since their birth in Africa and shows wide phylogeographic divergence. Multiple complete genome sequences were analyzed including sequences from Okinawa, Japan, that we recently sequenced. The newer method revealed a finer population structure than revealed by a previous method that examines only MLST housekeeping genes or a phylogenetic network analysis of the core genome. Novel subgroups were found in Europe, Amerind, and East Asia groups. Examination of genetic flux showed some singleton strains to be hybrids of subgroups and revealed evident signs of population admixture in Africa, Europe, and parts of Asia. We expect this approach to further our understanding of intraspecific bacterial evolution by revealing population structure at a finer scale.

Minimum core genome sequence typing of bacterial pathogens: a unified approach for clinical and public health microbiology

Bacterial pathogens impose a heavy health burden worldwide. In the new era of high-throughput sequencing and online bioinformatics, real-time genome typing of infecting agents, and in particular those with potential severe clinical outcomes, holds promise for guiding clinical care to limit the detrimental effects of infections and to prevent potential local or global outbreaks. Here, we sequenced and compared 85 isolates of Streptococcus suis, a zoonotic human and swine pathogen, wherein we analyzed 32 recognized serotypes and 75 sequence types representing the diversity of the species and the human clinical isolates with high public health significance. We found that 1,077 of the 2,469 genes are shared by all isolates. Excluding 201 common but mobile genes, 876 genes were defined as the minimum core genome (MCG) of the species. Of 190,894 single-nucleotide polymorphisms (SNPs) identified, 58,501 were located in the MCG genes and were referred to as MCG SNPs. A population structure analysis of these MCG SNPs classified the 85 isolates into seven MCG groups, of which MCG group 1 includes all isolates from human infections and outbreaks. Our MCG typing system for S. suis provided a clear separation of groups containing human-associated isolates from those containing animal-associated isolates. It also separated the group containing outbreak isolates, including those causing life-threatening streptococcal toxic shock-like syndrome, from sporadic or less severe meningitis or bacteremia-only isolates. The typing system facilitates the application of genome data to the fields of clinical medicine and epidemiology and to the surveillance of S. suis. The MCG groups may also be used as the taxonomical units of S. suis to define bacterial subpopulations with the potential to cause severe clinical infections and large-scale outbreaks.

Genetic Characterization of Legionella pneumophila Isolated from a Common Watershed in Comunidad Valenciana, Spain

Legionella pneumophila infects humans to produce legionellosis and Pontiac fever only from environmental sources. In order to establish control measures and study the sources of outbreaks it is essential to know extent and distribution of strain variants of this bacterium in the environment. Sporadic and outbreak-related cases of legionellosis have been historically frequent in the Comunidad Valenciana region (CV, Spain), with a high prevalence in its Southeastern-most part (BV). Environmental investigations for the detection of Legionella pneumophila are performed in this area routinely. We present a population genetics study of 87 L. pneumophila strains isolated in 13 different localities of the BV area irrigated from the same watershed and compare them to a dataset of 46 strains isolated in different points of the whole CV. Our goal was to compare environmental genetic variation at two different geographic scales, at county and regional levels. Genetic diversity, recombination and population structure were analyzed with Sequence-Based Typing data and three intergenic regions. The results obtained reveal a low, but detectable, level of genetic differentiation between both datasets, mainly, but not only, attributed to the occurrence of unusual variants of the neuA locus present in the BV populations. This differentiation is still detectable when the 10 loci considered are analyzed independently, despite the relatively high incidence of the most common genetic variant in this species, sequence type 1 (ST-1). However, when the genetic data are considered without their associated geographic information, four major groups could be inferred at the genetic level which did not show any correlation with sampling locations. The overall results indicate that the population structure of these environmental samples results from the joint action of a global, widespread ST-1 along with genetic differentiation at shorter geographic distances, which in this case are related to the common watershed for the BV localities.

Meta-analysis of pulsed-field gel electrophoresis fingerprints based on a constructed Salmonella database

A database was constructed consisting of 45,923 Salmonella pulsed-field gel electrophoresis (PFGE) patterns. The patterns, randomly selected from all submissions to CDC PulseNet during 2005 to 2010, included the 20 most frequent serotypes and 12 less frequent serotypes. Meta-analysis was applied to all of the PFGE patterns in the database. In the range of 20 to 1100 kb, serotype Enteritidis averaged the fewest bands at 12 bands and Paratyphi A the most with 19, with most serotypes in the 13−15 range among the 32 serptypes. The 10 most frequent bands for each of the 32 serotypes were sorted and distinguished, and the results were in concordance with those from distance matrix and two-way hierarchical cluster analyses of the patterns in the database. The hierarchical cluster analysis divided the 32 serotypes into three major groups according to dissimilarity measures, and revealed for the first time the similarities among the PFGE patterns of serotype Saintpaul to serotypes Typhimurium, Typhimurium var. 5-, and I 4,[5],12:i:-; of serotype Hadar to serotype Infantis; and of serotype Muenchen to serotype Newport. The results of the meta-analysis indicated that the pattern similarities/dissimilarities determined the serotype discrimination of PFGE method, and that the possible PFGE markers may have utility for serotype identification. The presence of distinct, serotype specific patterns may provide useful information to aid in the distribution of serotypes in the population and potentially reduce the need for laborious analyses, such as traditional serotyping.

Evolutionary Genomics of Salmonella enterica Subspecies

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

We have sequenced 21 new genomes which encompass the phylogenetic diversity of Salmonella, including strains of the previously unsequenced subspecies arizonae, diarizonae, houtenae, salamae, and indica as well as new diverse strains of subspecies enterica. We have deduced possible evolutionary paths traversed by this very important zoonotic pathogen and identified novel putative virulence factors that are not found in subspecies I. Gene families gained at the time of the evolution of subspecies enterica are of particular interest because they include mechanisms by which this subspecies adapted to warm-blooded hosts.

Phylogenetics and differentiation of Salmonella Newport lineages by whole genome sequencing

Salmonella Newport has ranked in the top three Salmonella serotypes associated with foodborne outbreaks from 1995 to 2011 in the United States. In the current study, we selected 26 S. Newport strains isolated from diverse sources and geographic locations and then conducted 454 shotgun pyrosequencing procedures to obtain 16–24 × coverage of high quality draft genomes for each strain. Comparative genomic analysis of 28 S. Newport strains (including 2 reference genomes) and 15 outgroup genomes identified more than 140,000 informative SNPs. A resulting phylogenetic tree consisted of four sublineages and indicated that S. Newport had a clear geographic structure. Strains from Asia were divergent from those from the Americas. Our findings demonstrated that analysis using whole genome sequencing data resulted in a more accurate picture of phylogeny compared to that using single genes or small sets of genes. We selected loci around the mutS gene of S. Newport to differentiate distinct lineages, including those between invH and mutS genes at the 3′ end of Salmonella Pathogenicity Island 1 (SPI-1), ste fimbrial operon, and Clustered, Regularly Interspaced, Short Palindromic Repeats (CRISPR) associated-proteins (cas). These genes in the outgroup genomes held high similarity with either S. Newport Lineage II or III at the same loci. S. Newport Lineages II and III have different evolutionary histories in this region and our data demonstrated genetic flow and homologous recombination events around mutS. The findings suggested that S. Newport Lineages II and III diverged early in the serotype evolution and have evolved largely independently. Moreover, we identified genes that could delineate sublineages within the phylogenetic tree and that could be used as potential biomarkers for trace-back investigations during outbreaks. Thus, whole genome sequencing data enabled us to better understand the genetic background of pathogenicity and evolutionary history of S. Newport and also provided additional markers for epidemiological response.

Reproductive clonality of pathogens: a perspective on pathogenic viruses, bacteria, fungi, and parasitic protozoa

We propose that clonal evolution in micropathogens be defined as restrained recombination on an evolutionary scale, with genetic exchange scarce enough to not break the prevalent pattern of clonal population structure, a definition already widely used for all kinds of pathogens, although not clearly formulated by many scientists and rejected by others. The two main manifestations of clonal evolution are strong linkage disequilibrium (LD) and widespread genetic clustering ("near-clading"). We hypothesize that this pattern is not mainly due to natural selection, but originates chiefly from in-built genetic properties of pathogens, which could be ancestral and could function as alternative allelic systems to recombination genes ("clonality/sexuality machinery") to escape recombinational load. The clonal framework of species of pathogens should be ascertained before any analysis of biomedical phenotypes (phylogenetic character mapping). In our opinion, this model provides a conceptual framework for the population genetics of any micropathogen.

Role of homologous recombination in adaptive diversification of extraintestinal Escherichia coli

The contribution of homologous exchange (recombination) of core genes in the adaptive evolution of bacterial pathogens is not well understood. To investigate this, we analyzed fully assembled genomes of two Escherichia coli strains from sequence type 131 (ST131), a clonal group that is both the leading cause of extraintestinal E. coli infections and the main source of fluoroquinolone-resistant E. coli. Although the sequences of each of the seven multilocus sequence typing genes were identical in the two ST131 isolates, the strains diverged from one another by homologous recombination that affected at least 9% of core genes. This was on a par with the contribution to genomic diversity of horizontal gene transfer and point gene mutation. The genomic positions of recombinant and mobile genetic regions were partially linked, suggesting their concurrent occurrence. One of the genes affected by homologous recombination was fimH, which encodes mannose-specific type 1 fimbrial adhesin, resulting in functionally distinct copies of the gene in ST131 strains. One strain, a uropathogenic isolate, had a pathoadaptive variant of fimH that was acquired by homologous replacement into the commensal strain background. Close examination of FimH structure and function in additional ST131 isolates revealed that recombination led to acquisition of several functionally distinct variants that, upon homologous exchange, were targeted by a variety of pathoadaptive mutations under strong positive selection. Different recombinant fimH strains also showed a strong clonal association with ST131 isolates that had distinct fluoroquinolone resistance profiles. Thus, homologous recombination of core genes plays a significant role in adaptive diversification of bacterial pathogens, especially at the level of clonally related groups of isolates.

Developing insights into the mechanisms of evolution of bacterial pathogens from whole-genome sequences

Evolution of bacterial pathogen populations has been detected in a variety of ways including phenotypic tests, such as metabolic activity, reaction to antisera and drug resistance and genotypic tests that measure variation in chromosome structure, repetitive loci and individual gene sequences. While informative, these methods only capture a small subset of the total variation and, therefore, have limited resolution. Advances in sequencing technologies have made it feasible to capture whole-genome sequence variation for each sample under study, providing the potential to detect all changes at all positions in the genome from single nucleotide changes to large-scale insertions and deletions. In this review, we focus on recent work that has applied this powerful new approach and summarize some of the advances that this has brought in our understanding of the details of how bacterial pathogens evolve.

From Exit to Entry: Long-term Survival and Transmission of Salmonella

Salmonella spp. are a leading cause of human infectious disease worldwide and pose a serious health concern. While we have an improving understanding of pathogenesis and the host-pathogen interactions underlying the infection process, comparatively little is known about the survival of pathogenic Salmonella outside their hosts. This review focuses on three areas: (1) in vitro evidence that Salmonella spp. can survive for long periods of time under harsh conditions; (2) observations and conclusions about Salmonella persistence obtained from human outbreaks; and (3) new information revealed by genomic- and population-based studies of Salmonella and related enteric pathogens. We highlight the mechanisms of Salmonella persistence and transmission as an essential part of their lifecycle and a prerequisite for their evolutionary success as human pathogens.

Human infections attributable to the D-tartrate-fermenting variant of Salmonella enterica serovar Paratyphi B in Germany originate in reptiles and, on rare occasions, poultry

In this study, the population structure, incidence, and potential sources of human infection caused by the d-tartrate-fermenting variant of Salmonella enterica serovar Paratyphi B [S. Paratyphi B (dT+)] was investigated. In Germany, the serovar is frequently isolated from broilers. Therefore, a selection of 108 epidemiologically unrelated S. enterica serovar Paratyphi B (dT+) strains isolated in Germany between 2002 and 2010 especially from humans, poultry/poultry meat, and reptiles was investigated by phenotypic and genotypic methods. Strains isolated from poultry and products thereof were strongly associated with multilocus sequence type ST28 and showed antimicrobial multiresistance profiles. Pulsed-field gel electrophoresis XbaI profiles were highly homogeneous, with only a few minor XbaI profile variants. All strains isolated from reptiles, except one, were strongly associated with ST88, another distantly related type. Most of the strains were susceptible to antimicrobial agents, and XbaI profiles were heterogeneous. Strains isolated from humans yielded seven sequence types (STs) clustering in three distantly related lineages. The first lineage, comprising five STs, represented mainly strains belonging to ST43 and ST149. The other two lineages were represented only by one ST each, ST28 and ST88. The relatedness of strains based on the pathogenicity gene repertoire (102 markers tested) was mostly in agreement with the multilocus sequence type. Because ST28 was frequently isolated from poultry but rarely in humans over the 9-year period investigated, overall, this study indicates that in Germany S. enterica serovar Paratyphi B (dT+) poses a health risk preferentially by contact with reptiles and, to a less extent, by exposure to poultry or poultry meat.

Pan-genomic analysis provides insights into the genomic variation and evolution of Salmonella Paratyphi A

Salmonella Paratyphi A (S. Paratyphi A) is a highly adapted, human-specific pathogen that causes paratyphoid fever. Cases of paratyphoid fever have recently been increasing, and the disease is becoming a major public health concern, especially in Eastern and Southern Asia. To investigate the genomic variation and evolution of S. Paratyphi A, a pan-genomic analysis was performed on five newly sequenced S. Paratyphi A strains and two other reference strains. A whole genome comparison revealed that the seven genomes are collinear and that their organization is highly conserved. The high rate of substitutions in part of the core genome indicates that there are frequent homologous recombination events. Based on the changes in the pan-genome size and cluster number (both in the core functional genes and core pseudogenes), it can be inferred that the sharply increasing number of pseudogene clusters may have strong correlation with the inactivation of functional genes, and indicates that the S. Paratyphi A genome is being degraded.

Large plasmids of Escherichia coli and Salmonella encode highly diverse arrays of accessory genes on common replicon families

Plasmids are important in evolution and adaptation of host bacteria, yet we lack a comprehensive picture of their own natural variation. We used replicon typing and RFLP analysis to assess diversity and distribution of plasmids in the ECOR, SARA, SARB and SARC reference collections of Escherichia coli and Salmonella. Plasmids, especially large (≥30 kb) plasmids, are abundant in these collections. Host species and genotype clearly impact plasmid prevalence; plasmids are more abundant in ECOR than SAR, but, within ECOR, subgroup B2 strains have the fewest large plasmids. The majority of large plasmids have unique RFLP patterns, suggesting high variation, even within dominant replicon families IncF and IncI1. We found only four conserved plasmid types within ECOR, none of which are widely distributed. Within SAR, conserved plasmid types are primarily serovar-specific, including a pSLT-like plasmid in 13 Typhimurium strains. Conservation of pSLT contrasts with variability of other plasmids, suggesting evolution of serovar-specific virulence plasmids is distinct from that of most enterobacterial plasmids. We sequenced a conserved serovar Heidelberg plasmid but did not detect virulence or antibiotic resistance genes. Our data illustrate the high degree of natural variation in large plasmids of E. coli and Salmonella, even among plasmids sharing backbone genes.

Population genomics of Chlamydia trachomatis: insights on drift, selection, recombination, and population structure

The large number of sexually transmitted diseases and ocular trachoma cases that are caused globally each year by Chlamydia trachomatis has made this organism a World Health Organization priority for vaccine development. However, there is no gene transfer system for Chlamydia to help identify potential vaccine targets. To accelerate discoveries toward this goal, here we analyzed the broadest diversity of C. trachomatis genomes to date, including 25 geographically dispersed clinical and seven reference strains representing 14 of the 19 known serotypes. Strikingly, all 32 genomes were found to have evidence of DNA acquisition by homologous recombination in their history. Four distinct clades were identified, which correspond to all C. trachomatis disease phenotypes: lymphogranuloma venereum (LGV; Clade 1); noninvasive urogenital infections (Clade 2); ocular trachoma (Clade 3); and protocolitis (Clade 4; also includes some noninvasive urogenital infections). Although the ancestral relationship between clades varied, most strains acted as donor and recipient of recombination with no evidence for barriers to genetic exchange. The niche-specific LGV and trachoma clades have undergone less recombination, although the opportunity for mixing with strains from other clades that infect the rectal and ocular mucosa, respectively, is evident. Furthermore, there are numerous occasions for gene conversion events through sequential infections at the same anatomic sites. The size of recombinant segments is relatively small (~357 bp) compared with in vitro experiments of various C. trachomatis strains but is consistent with in vitro estimates for other bacterial species including Escherichia coli and Helicobacter pylori. Selection has also played a crucial role during the diversification of the organism. Clade 2 had the lowest nonsynonymous to synonymous ratio (dN/dS) but the highest effect of recombination, which is consistent with the widespread occurrence of synonymous substitutions in recombined genomic segments. The trachoma Clade 3 had the highest dN/dS estimates, which may be caused by an increased effect of genetic drift from niche specialization and a reduced effective population size. The degree of drift, selection, and recombination in C. trachomatis suggests that the challenge will remain to identify genomic regions that are stable and cross protective for the development of an efficacious vaccine.

Identification and characterization of novel Salmonella mobile elements involved in the dissemination of genes linked to virulence and transmission

The genetic diversity represented by >2,500 different Salmonella serovars provides a yet largely uncharacterized reservoir of mobile elements that can contribute to the frequent emergence of new pathogenic strains of this important zoonotic pathogen. Currently, our understanding of Salmonella mobile elements is skewed by the fact that most studies have focused on highly virulent or common serovars. To gain a more global picture of mobile elements in Salmonella, we used prediction algorithms to screen for mobile elements in 16 sequenced Salmonella genomes representing serovars for which no prior genome scale mobile element data were available. From these results, selected mobile elements underwent further analyses in the form of validation studies, comparative analyses, and PCR-based population screens. Through this analysis we identified a novel plasmid that has two cointegrated replicons (IncI1-IncFIB); this plasmid type was found in four genomes representing different Salmonella serovars and contained a virulence gene array that had not been previously identified. A Salmonella Montevideo isolate contained an IncHI and an IncN2 plasmid, which both encoded antimicrobial resistance genes. We also identified two novel genomic islands (SGI2 and SGI3), and 42 prophages with mosaic architecture, seven of them harboring known virulence genes. Finally, we identified a novel integrative conjugative element (ICE) encoding a type IVb pilus operon in three non-typhoidal Salmonella serovars. Our analyses not only identified a considerable number of mobile elements that have not been previously reported in Salmonella, but also found evidence that these elements facilitate transfer of genes that were previously thought to be limited in their distribution among Salmonella serovars. The abundance of mobile elements encoding pathogenic properties may facilitate the emergence of strains with novel combinations of pathogenic traits.

Convergent molecular evolution of genomic cores in Salmonella enterica and Escherichia coli

One of the strongest signals of adaptive molecular evolution of proteins is the occurrence of convergent hot spot mutations: repeated changes in the same amino acid positions. We performed a comparative genome-wide analysis of mutation-driven evolution of core (omnipresent) genes in 17 strains of Salmonella enterica subspecies I and 22 strains of Escherichia coli. More than 20% of core genes in both Salmonella and E. coli accumulated hot spot mutations, with a predominance of identical changes having recent evolutionary origin. There is a significant overlap in the functional categories of the adaptively evolving genes in both species, although mostly via separate molecular mechanisms. As a strong evidence of the link between adaptive mutations and virulence in Salmonella, two human-restricted serovars, Typhi and Paratyphi A, shared the highest number of genes with serovar-specific hot spot mutations. Many of the core genes affected by Typhi/Paratyphi A-specific mutations have known virulence functions. For each species, a list of nonrecombinant core genes (and the hot spot mutations therein) under positive selection is provided.

Impact of loci nature on estimating recombination and mutation rates in Chlamydia trachomatis

The knowledge of the frequency and relative weight of mutation and recombination events in evolution is essential for understanding how microorganisms reach fitted phenotypes. Traditionally, these evolutionary parameters have been inferred by using data from multilocus sequence typing (MLST), which is known to have yielded conflicting results. In the near future, these estimations will certainly be performed by computational analyses of full-genome sequences. However, it is not known whether this approach will yield accurate results as bacterial genomes exhibit heterogeneous representation of loci categories, and it is not clear how loci nature impacts such estimations. Therefore, we assessed how mutation and recombination inferences are shaped by loci with different genetic features, using the bacterium Chlamydia trachomatis as the study model. We found that loci assigning a high number of alleles and positively selected genes yielded nonconvergent estimates and incongruent phylogenies and thus are more prone to confound algorithms. Unexpectedly, for the model under evaluation, housekeeping genes and noncoding regions shaped estimations in a similar manner, which points to a nonrandom role of the latter in C. trachomatis evolution. Although the present results relate to a specific bacterium, we speculate that microbe-specific genomic architectures (such as coding capacity, polymorphism dispersion, and fraction of positively selected loci) may differentially buffer the effect of the confounding factors when estimating recombination and mutation rates and, thus, influence the accuracy of using full-genome sequences for such purpose. This putative bias associated with in silico inferences should be taken into account when discussing the results obtained by the analyses of full-genome sequences, in which the “one size fits all” approach may not be applicable.

Multilocus sequence typing as a replacement for serotyping in Salmonella enterica

Salmonella enterica subspecies enterica is traditionally subdivided into serovars by serological and nutritional characteristics. We used Multilocus Sequence Typing (MLST) to assign 4,257 isolates from 554 serovars to 1092 sequence types (STs). The majority of the isolates and many STs were grouped into 138 genetically closely related clusters called eBurstGroups (eBGs). Many eBGs correspond to a serovar, for example most Typhimurium are in eBG1 and most Enteritidis are in eBG4, but many eBGs contained more than one serovar. Furthermore, most serovars were polyphyletic and are distributed across multiple unrelated eBGs. Thus, serovar designations confounded genetically unrelated isolates and failed to recognize natural evolutionary groupings. An inability of serotyping to correctly group isolates was most apparent for Paratyphi B and its variant Java. Most Paratyphi B were included within a sub-cluster of STs belonging to eBG5, which also encompasses a separate sub-cluster of Java STs. However, diphasic Java variants were also found in two other eBGs and monophasic Java variants were in four other eBGs or STs, one of which is in subspecies salamae and a second of which includes isolates assigned to Enteritidis, Dublin and monophasic Paratyphi B. Similarly, Choleraesuis was found in eBG6 and is closely related to Paratyphi C, which is in eBG20. However, Choleraesuis var. Decatur consists of isolates from seven other, unrelated eBGs or STs. The serological assignment of these Decatur isolates to Choleraesuis likely reflects lateral gene transfer of flagellar genes between unrelated bacteria plus purifying selection. By confounding multiple evolutionary groups, serotyping can be misleading about the disease potential of S. enterica. Unlike serotyping, MLST recognizes evolutionary groupings and we recommend that Salmonella classification by serotyping should be replaced by MLST or its equivalents.

Microbiologists have used serological and nutritional characteristics to subdivide pathogenic bacteria for nearly 100 years. These subdivisions in Salmonella enterica are called serovars, some of which are thought to be associated with particular diseases and epidemiology. We used MultiLocus Sequence-based Typing (MLST) to identify clusters of S. enterica isolates that are related by evolutionary descent. Some clusters correspond to serovars on a one to one basis. But many clusters include multiple serovars, which is of public health significance, and most serovars span multiple, unrelated clusters. Despite its broad usage, serological typing of S. enterica has resulted in confusing systematics, with a few exceptions. We recommend that serotyping for strain discrimination of S. enterica be replaced by a DNA-based method, such as MLST. Serotyping and other non-sequence based typing methods are routinely used for detecting outbreaks and to support public health responses. Moving away from these methods will require a major shift in thinking by public health microbiology laboratories as well as national and international agencies. However, a transition to the routine use of MLST, supplemented where appropriate by even more discriminatory sequence-based typing methods based on entire genomes, will provide a clearer picture of long-term transmission routes of Salmonella, facilitate data transfer and support global control measures.

Impact of homologous and non-homologous recombination in the genomic evolution of Escherichia coli

Background

Escherichia coli is an important species of bacteria that can live as a harmless inhabitant of the guts of many animals, as a pathogen causing life-threatening conditions or freely in the non-host environment. This diversity of lifestyles has made it a particular focus of interest for studies of genetic variation, mainly with the aim to understand how a commensal can become a deadly pathogen. Many whole genomes of E. coli have been fully sequenced in the past few years, which offer helpful data to help understand how this important species evolved.

Results

We compared 27 whole genomes encompassing four phylogroups of Escherichia coli (A, B1, B2 and E). From the core-genome we established the clonal relationships between the isolates as well as the role played by homologous recombination during their evolution from a common ancestor. We found strong evidence for sexual isolation between three lineages (A+B1, B2, E), which could be explained by the ecological structuring of E. coli and may represent on-going speciation. We identified three hotspots of homologous recombination, one of which had not been previously described and contains the aroC gene, involved in the essential shikimate metabolic pathway. We also described the role played by non-homologous recombination in the pan-genome, and showed that this process was highly heterogeneous. Our analyses revealed in particular that the genomes of three enterohaemorrhagic (EHEC) strains within phylogroup B1 have converged from originally separate backgrounds as a result of both homologous and non-homologous recombination.

Conclusions

Recombination is an important force shaping the genomic evolution and diversification of E. coli, both by replacing fragments of genes with an homologous sequence and also by introducing new genes. In this study, several non-random patterns of these events were identified which correlated with important changes in the lifestyle of the bacteria, and therefore provide additional evidence to explain the relationship between genomic variation and ecological adaptation.