Bordetella pertussis causes pertussis, a respiratory disease that is most severe for infants. Vaccination was introduced in the 1950s, and in recent years, a resurgence of disease was observed worldwide, with significant mortality in infants. Possible causes for this include the switch from whole-cell vaccines (WCVs) to less effective acellular vaccines (ACVs), waning immunity, and pathogen adaptation. Pathogen adaptation is suggested by antigenic divergence between vaccine strains and circulating strains and by the emergence of strains with increased pertussis toxin production. We applied comparative genomics to a worldwide collection of 343 B. pertussis strains isolated between 1920 and 2010. The global phylogeny showed two deep branches; the largest of these contained 98% of all strains, and its expansion correlated temporally with the first descriptions of pertussis outbreaks in Europe in the 16th century. We found little evidence of recent geographical clustering of the strains within this lineage, suggesting rapid strain flow between countries. We observed that changes in genes encoding proteins implicated in protective immunity that are included in ACVs occurred after the introduction of WCVs but before the switch to ACVs. Furthermore, our analyses consistently suggested that virulence-associated genes and genes coding for surface-exposed proteins were involved in adaptation. However, many of the putative adaptive loci identified have a physiological role, and further studies of these loci may reveal less obvious ways in which B. pertussis and the host interact. This work provides insight into ways in which pathogens may adapt to vaccination and suggests ways to improve pertussis vaccines.

The evolution of new combinations of bacterial properties contributes to biodiversity and the emergence of new diseases. We investigated the capacity for bacterial divergence with a chemostat culture of Escherichia coli. A clonal population radiated into more than five phenotypic clusters within 26 days, with multiple variations in global regulation, metabolic strategies, surface properties, and nutrient permeability pathways. Most isolates belonged to a single ecotype, and neither periodic selection events nor ecological competition for a single niche prevented an adaptive radiation with a single resource. The multidirectional exploration of fitness space is an underestimated ingredient to bacterial success even in unstructured environments.

Environmental stresses increase genetic variation in bacteria, plants, and human cancer cells. The linkage between various environments and mutational outcomes has not been systematically investigated, however. Here, we established the influence of nutritional stresses commonly found in the biosphere (carbon, phosphate, nitrogen, oxygen, or iron limitation) on both the rate and spectrum of mutations in Escherichia coli. We found that each limitation was associated with a remarkably distinct mutational profile. Overall mutation rates were not always elevated, and nitrogen, iron, and oxygen limitation resulted in major spectral changes but no net increase in rate. Our results thus suggest that stress-induced mutagenesis is a diverse series of stress input–mutation output linkages that is distinct in every condition. Environment-specific spectra resulted in the differential emergence of traits needing particular mutations in these settings. Mutations requiring transpositions were highest under iron and oxygen limitation, whereas base-pair substitutions and indels were highest under phosphate limitation. The unexpected diversity of input–output effects explains some important phenomena in the mutational biases of evolving genomes. The prevalence of bacterial insertion sequence transpositions in the mammalian gut or in anaerobically stored cultures is due to environmentally determined mutation availability. Likewise, the much-discussed genomic bias towards transition base substitutions in evolving genomes can now be explained as an environment-specific output. Altogether, our conclusion is that environments influence genetic variation as well as selection.

Despite high vaccine coverage, pertussis incidence has increased substantially in recent years in many countries. A significant factor that may be contributing to this increase is adaptation to the vaccine by Bordetella pertussis, the causative agent of pertussis. In this study, we first assessed the genetic diversity of B. pertussis by microarray-based comparative genome sequencing of 10 isolates representing diverse genotypes and different years of isolation. We discovered 171 single nucleotide polymorphisms (SNPs) in a total of 1.4 Mb genome analyzed. The frequency of base changes was estimated as one per 32 kb per isolate, confirming that B. pertussis is one of the least variable bacterial pathogens. We then analyzed an international collection of 316 B. pertussis isolates using a subset of 65 of the SNPs and identified 42 distinct SNP profiles (SPs). Phylogenetic analysis grouped the SPs into six clusters. The majority of recent isolates belonged to clusters I-IV and were descendants of a single prevaccine lineage. Cluster I appeared to be a major clone with a worldwide distribution. Typing of genes encoding acellular vaccine (ACV) antigens, ptxA, prn, fhaB, fim2, and fim3 revealed the emergence and increasing incidence of non-ACV alleles occurring in clusters I and IV, which may have been driven by ACV immune selection. Our findings suggest that B. pertussis, despite its high population homogeneity, is evolving in response to vaccination pressure with recent expansion of clones carrying variants of genes encoding ACV antigens.

Many of the important changes in evolution are regulatory in nature. Sequenced bacterial genomes point to flexibility in regulatory circuits but we do not know how regulation is remodeled in evolving bacteria. Here, we study the regulatory changes that emerge in populations evolving under controlled conditions during experimental evolution of Escherichia coli in a phosphate-limited chemostat culture. Genomes were sequenced from five clones with different combinations of phenotypic properties that coexisted in a population after 37 days. Each of the distinct isolates contained a different mutation in 1 of 3 highly pleiotropic regulatory genes (hfq, spoT, or rpoS). The mutations resulted in dissimilar proteomic changes, consistent with the documented effects of hfq, spoT, and rpoS mutations. The different mutations do share a common benefit, however, in that the mutations each redirect cellular resources away from stress responses that are redundant in a constant selection environment. The hfq mutation lowers several individual stress responses as well the small RNA–dependent activation of rpoS translation and hence general stress resistance. The spoT mutation reduces ppGpp levels, decreasing the stringent response as well as rpoS expression. The mutations in and upstream of rpoS resulted in partial or complete loss of general stress resistance. Our observations suggest that the degeneracy at the core of bacterial stress regulation provides alternative solutions to a common evolutionary challenge. These results can explain phenotypic divergence in a constant environment and also how evolutionary jumps and adaptive radiations involve altered gene regulation.