DksA is a conserved master regulator of stress response in Acinetobacter baumannii.

Coordination of bacterial stress response mechanisms is critical for long-term survival in harsh environments for successful host infection. The general and specific stress responses of well-studied Gram-negative pathogens like Escherichia coli are controlled by alternative sigma factors, archetypically RpoS. The deadly hospital pathogen Acinetobacter baumannii is notoriously resistant to environmental stresses, yet it lacks RpoS, and the molecular mechanisms driving this incredible stress tolerance remain poorly defined. Here, using functional genomics, we identified the transcriptional regulator DksA as a master regulator for broad stress protection and virulence in A. baumannii. Transcriptomics, phenomics and in vivo animal studies revealed that DksA controls ribosomal protein expression, metabolism, mutation rates, desiccation, antibiotic resistance, and host colonization in a niche-specific manner. Phylogenetically, DksA was highly conserved and well-distributed across Gammaproteobacteria, with 96.6% containing DksA, spanning 88 families. This study lays the groundwork for understanding DksA as a major regulator of general stress response and virulence in this important pathogen.

Genomic and phenotypic analyses of diverse non-clinical Acinetobacter baumannii strains reveals strain-specific virulence and resistance capacity.

Acinetobacter baumannii is a critically important pathogen known for its widespread antibiotic resistance and ability to persist in hospital-associated environments. Whilst the majority of A. baumannii infections are hospital-acquired, infections from outside the hospital have been reported with high mortality. Despite this, little is known about the natural environmental reservoir(s) of A. baumannii and the virulence potential underlying non-clinical strains. Here, we report the complete genome sequences of six diverse strains isolated from environments such as river, soil, and industrial sites around the world. Phylogenetic analyses showed that four of these strains were unrelated to representative nosocomial strains and do not share a monophyletic origin, whereas two had sequence types belonging to the global clone lineages GC1 and GC2. Further, the majority of these strains harboured genes linked to virulence and stress protection in nosocomial strains. These genotypic properties correlated well with in vitro virulence phenotypic assays testing resistance to abiotic stresses, serum survival, and capsule formation. Virulence potential was confirmed in vivo, with most environmental strains able to effectively kill Galleria mellonella greater wax moth larvae. Using phenomic arrays and antibiotic resistance profiling, environmental and nosocomial strains were shown to have similar substrate utilisation patterns although environmental strains were distinctly more sensitive to antibiotics. Taken together, these features of environmental A. baumannii strains suggest the existence of a strain-specific distinct gene pools for niche specific adaptation. Furthermore, environmental strains appear to be equally virulent as contemporary nosocomial strains but remain largely antibiotic sensitive.

Escherichia coli mutation rates and spectra with combinations of environmental limitations.

Micro-organisms often face multiple stresses in natural habitats. Individual stresses are well known to influence mutation rates and the spectra of mutational types, but the extent to which multiple stresses affect the genetic variation in populations is unknown. Here we investigate pair-wise combinations of nutritional stresses in Escherichia coli to determine their effect on mutation rates and mutational types. Environmental interactions modified both the rate and spectrum of mutations in double-limited environments, but the effects were not additive or synergistic relative to single stresses. Generally, bacteria in the mixed environments behaved as if one of the two single-stress stimuli was more dominant and the genetic variation seen with every dual limitation was intermediate between known patterns with individual stresses. The composition of mutational types with double stresses was also intermediate between individual stress patterns. At least with mutations, the single stressor results available are reasonable indicators of stress-induced genetic variation in multifaceted natural habitats. With the influence of 11 conditions available on mutational patterns, we can now also see the clustering of mutational types as a function of these environments.

The impact of growth rate and environmental factors on mutation rates and spectra in Escherichia coli.

Genetic variation in bacterial populations is remarkably sensitive to environmental influences, including simple, nutritional differences. Not only the rate but also the kind of mutational changes is biased by the nutritional state of bacteria. Here we investigate the mutational consequences of a universal variable for free-living bacteria, namely the growth rate. By controlling growth in chemostats, the rate and mix of mutations was investigated for populations of Escherichia coli subject to different specific growth rates. Both aerobic and anaerobic cultures were compared with see if growth rate is a factor in the commonest respiratory conditions for E. coli. We find mutation rates are raised markedly with decreasing growth rate. Base pair substitutions and 1-bp insertions and deletions increase with reduced growth rate, but less so in anaerobic cultures. Insertion sequence movements are particularly sensitive to growth rate, with IS2 being optimal at intermediate growth rates whereas IS1 and IS150 movements are highest at the slowest tested growth rate. A comprehensive comparison of growth rate effects, as well as six other environmental factors, provides the most complete picture yet of the range of mutational signatures in bacterial genetic variation.

A shifting mutational landscape in 6 nutritional states: Stress-induced mutagenesis as a series of distinct stress input-mutation output relationships.

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.

Natural Escherichia coli isolates rapidly acquire genetic changes upon laboratory domestication.

The adaptation of environmental bacteria to laboratory conditions was analysed through the exploration of genomic changes in four strains of Escherichia coli freshly isolated from their natural habitats and belonging to different taxonomic clusters. Up to 25 mutations were present in all cultures of natural isolates within 10 days of transfer in rich media or with a single growth cycle involving an extended stationary phase. Among numerous individual mutations, two genes were affected in parallel in distinct backgrounds. Mutations in rpoS (encoding sigma factor RpoS), altering a multiplication-survival trade-off in E. coli, were present in isolates derived from all four different ancestors. More surprisingly, two different natural isolates acquired mutations in mutL, affecting DNA mismatch repair, and a third also involved higher mutation rates. The elevated mutation rates in these isolates indicate the danger of increased genetic instability arising from laboratory domestication. Neither rpoS nor mutator mutations were detected in the already-acclimatized MG1655 laboratory strain; only one or no new mutations were present in the laboratory strain under the same culture conditions. Our results indicate rapid adaptation to the laboratory environment. Ancestor-specific responses also arise in the laboratory and mutational events are also sensitive to culture conditions such as extended stationary phase. To maintain natural isolates in a stable state, our data suggest that the transition of strains to the laboratory should minimize culture cycles and extended stationary phase.

Simple phenotypic sweeps hide complex genetic changes in populations.

Changes in allele frequencies and the fixation of beneficial mutations are central to evolution. The precise relationship between mutational and phenotypic sweeps is poorly described however, especially when multiple alleles are involved. Here, we investigate these relationships in a bacterial population over 60 days in a glucose-limited chemostat in a large population. High coverage metagenomic analysis revealed a disconnection between smooth phenotypic sweeps and the complexity of genetic changes in the population. Phenotypic adaptation was due to convergent evolution and involved soft sweeps by 7-26 highly represented alleles of several genes in different combinations. Allele combinations spread from undetectably low baselines, indicating that minor subpopulations provide the basis of most innovations. A hard sweep was also observed, involving a single combination of rpoS, mglD, malE, sdhC, and malT mutations sweeping to greater than 95% of the population. Other mutant genes persisted but at lower abundance, including hfq, consistent with its demonstrated frequency-dependent fitness under glucose limitation. Other persistent, newly identified low-frequency mutations were in the aceF, galF, ribD and asm genes, in noncoding regulatory regions, three large indels and a tandem duplication; these were less affected by fluctuations involving more dominant mutations indicating separate evolutionary paths. Our results indicate a dynamic subpopulation structure with a minimum of 42 detectable mutations maintained over 60 days. We also conclude that the massive population-level mutation supply in combination with clonal interference leads to the soft sweeps observed, but not to the exclusion of an occasional hard sweep.

A metabolic trade-off between phosphate and glucose utilization in Escherichia coli.

Getting the most out of available nutrients is a key challenge that all organisms face. Little is known about how they optimize and balance the simultaneous utilization of multiple elemental resources. We investigated the effects of long-term phosphate limitation on carbon metabolism of the model organism Escherichia coli using chemostat cultures. We profiled metabolic changes in the growth medium over time and found evidence for an increase in fermentative metabolism despite the aerobic conditions. Using full-genome sequencing and competition experiments, we found that fitness under phosphate-limiting conditions was reproducibly increased by a mutation preventing flux through succinate in the tricarboxylic acid cycle. In contrast, these mutations reduced competitive ability under carbon limitation, and thus reveal a conflicting metabolic benefit in the role of the TCA cycle in environments limited by inorganic phosphate and glucose.

Global population structure and evolution of Bordetella pertussis and their relationship with vaccination.

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. IMPORTANCE Whooping cough is mainly caused by Bordetella pertussis, and current vaccines are targeted against this organism. Recently, there have been increasing outbreaks of whooping cough, even where vaccine coverage is high. Analysis of the genomes of 343 B. pertussis isolates from around the world over the last 100 years suggests that the organism has emerged within the last 500 years, consistent with historical records. We show that global transmission of new strains is very rapid and that the worldwide population of B. pertussis is evolving in response to vaccine introduction, potentially enabling vaccine escape.

A case of adaptation through a mutation in a tandem duplication during experimental evolution in Escherichia coli.

BACKGROUND: DNA duplications constitute important precursors for genome variation. Here we analyzed an unequal duplication harboring a beneficial mutation that may provide alternative evolutionary outcomes. RESULTS: We characterized this evolutionary event during experimental evolution for only 100 generations of an Escherichia coli strain under glucose limitation within chemostats. By combining Insertion Sequence based Restriction Length Polymorphism experiments, pulsed field gel electrophoresis and two independent genome re-sequencing experiments, we identified an evolved lineage carrying a 180 kb duplication of the 46' region of the E. coli chromosome. This evolved duplication revealed a heterozygous state, with one copy harboring a 2668 bp deletion that included part of the ogrK gene and both the yegR and yegS genes. By genetically manipulating ancestral and evolved strains, we showed that the single yegS inactivation was sufficient to confer a frequency dependent fitness increase under the chemostat selective conditions in both the ancestor and evolved genetic contexts, implying that the duplication itself was not a direct fitness contributor. Nonetheless, the heterozygous duplicated state was relatively stable in the conditions prevailing during evolution in chemostats, in striking contrast to non selective conditions in which the duplication resolved at high frequency into either its ancestral or deleted copy. CONCLUSIONS: Our results suggest that the duplication state may constitute a second order selection process providing higher evolutionary potential. Moreover, its heterozygous nature may provide differential evolutionary opportunities in alternating environments. Our results also highlighted how careful analyses of whole genome data are needed to identify such complex rearrangements.

Mutation accumulation and fitness in mutator subpopulations of Escherichia coli.

Bacterial populations in clinical and laboratory settings contain a significant proportion of mutants with elevated mutation rates (mutators). Mutators have a particular advantage when multiple beneficial mutations are needed for fitness, as in antibiotic resistance. Nevertheless, high mutation rates potentially lead to increasing numbers of deleterious mutations and subsequently to the decreased fitness of mutators. To test how fitness changed with mutation accumulation, genome sequencing and fitness assays of nine Escherichia coli mutY mutators were undertaken in an evolving chemostat population at three time points. Unexpectedly, the fitness in members of the mutator subpopulation became constant despite a growing number of mutations over time. To test if the accumulated mutations affected fitness, we replaced each of the known beneficial mutations with wild-type alleles in a mutator isolate. We found that the other 25 accumulated mutations were not deleterious. Our results suggest that isolates with deleterious mutations are eliminated by competition in a continuous culture, leaving mutators with mostly neutral mutations. Interestingly, the mutator-non-mutator balance in the population reversed after the fitness plateau of mutators was reached, suggesting that the mutator-non-mutator ratio in populations has more to do with competition between members of the population than the accumulation of deleterious mutations.

The multiplicity of divergence mechanisms in a single evolving population.

BACKGROUND: Evolutionary divergence is common within bacterial species and populations, even during a single bacterial infection. We use large-scale genomic and phenotypic analysis to identify the extent of diversification in controlled experimental populations and apply these data to differentiate between several potential mechanisms of evolutionary divergence. RESULTS: We defined testable differences between five proposed mechanisms and used experimental evolution studies to follow eight glucose-limited Escherichia coli chemostat populations at two growth rates. Simple phenotypic tests identified 11 phenotype combinations evolving under glucose limitation. Each evolved population exhibited 3 to 5 different combinations of the 11 phenotypic clusters. Genome sequencing of a representative of each phenotypic cluster from each population identified 193 mutations in 48 isolates. Only two of the 48 strains had evolved identically. Convergent paths to the same phenotype occurred, but two pleiotropic mutations were unique to slow-growing bacteria, permitting them greater phenotypic variance. Indeed, greater diversity arose in slower-growing, more stressed cultures. Mutation accumulation, hypermutator presence and fitness mechanisms varied between and within populations, with the evolved fitness considerably more uniform with fast growth cultures. Negative frequency-dependent fitness was shown by a subset of isolates. CONCLUSIONS: Evolutionary diversity is unlikely to be explained by any one of the available mechanisms. For a large population as used in this study, our results suggest that multiple mechanisms contribute to the mix of phenotypes and evolved fitness types in a diversifying population. Another major conclusion is that the capacity of a population to diversify is a function of growth rate.

The uncertain consequences of transferring bacterial strains between laboratories - rpoS instability as an example.

BACKGROUND: Microbiological studies frequently involve exchanges of strains between laboratories and/or stock centers. The integrity of exchanged strains is vital for archival reasons and to ensure reproducible experimental results. For at least 50 years, one of the most common means of shipping bacteria was by inoculating bacterial samples in agar stabs. Long-term cultures in stabs exhibit genetic instabilities and one common instability is in rpoS. The sigma factor RpoS accumulates in response to several stresses and in the stationary phase. One consequence of RpoS accumulation is the competition with the vegetative sigma factor σ70. Under nutrient limiting conditions mutations in rpoS or in genes that regulate its expression tend to accumulate. Here, we investigate whether short-term storage and mailing of cultures in stabs results in genetic heterogeneity. RESULTS: We found that samples of the E. coli K-12 strain MC4100TF exchanged on three separate occasions by mail between our laboratories became heterogeneous. Reconstruction studies indicated that LB-stabs exhibited mutations previously found in GASP studies in stationary phase LB broth. At least 40% of reconstructed stocks and an equivalent proportion of actually mailed stock contained these mutations. Mutants with low RpoS levels emerged within 7 days of incubation in the stabs. Sequence analysis of ten of these segregants revealed that they harboured each of three different rpoS mutations. These mutants displayed the classical phenotypes of bacteria lacking rpoS. The genetic stability of MC4100TF was also tested in filter disks embedded in glycerol. Under these conditions, GASP mutants emerge only after a 3-week period. We also confirm that the intrinsic high RpoS level in MC4100TF is mainly due to the presence of an IS1 insertion in rssB. CONCLUSIONS: Given that many E. coli strains contain high RpoS levels similar to MC4100TF, the integrity of such strains during transfers and storage is questionable. Variations in important collections may be due to storage-transfer related issues. These results raise important questions on the integrity of bacterial archives and transferred strains, explain variation like in the ECOR collection between laboratories and indicate a need for the development of better methods of strain transfer.

Insertion sequence-driven evolution of Escherichia coli in chemostats.

Insertion sequence (IS) elements are present in almost all bacterial genomes and are mobile enough to provide genomic tools to differentiate closely related isolates. They can be used to estimate genetic diversity and identify fitness-enhancing mutations during evolution experiments. Here, we determined the genomic distribution of eight IS elements in 120 genomes sampled from Escherichia coli populations that evolved in glucose- and phosphate-limited chemostats by comparison to the ancestral pattern. No significant differential transposition of the various IS types was detected across the environments. The phylogenies revealed substantial diversity amongst clones sampled from each chemostat, consistent with the phenotypic diversity within populations. Two IS-related changes were common to independent chemostats, suggesting parallel evolution. One of them corresponded to insertions of IS1 elements within rpoS encoding the master regulator of stress conditions. The other parallel event was an IS5-dependent deletion including mutY involved in DNA repair, thereby providing the molecular mechanism of generation of mutator clones in these evolving populations. These deletions occurred in different co-existing genotypes within single populations and were of various sizes. Moreover, differential locations of IS elements combined with their transpositional activity provided evolved clones with different phenotypic landscapes. Therefore, IS elements strongly influenced the evolutionary processes in continuous E. coli cultures by providing a way to modify both the global regulatory network and the mutation rates of evolving cells.

Insight into evolution of Bordetella pertussis from comparative genomic analysis: evidence of vaccine-driven 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.

Divergence involving global regulatory gene mutations in an Escherichia coli population evolving under phosphate limitation.

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.

Bordetella pertussis clones identified by multilocus variable-number tandem-repeat analysis.

Multilocus variable-number tandem-repeat analysis (MLVA) of 316 Bordetella pertussis isolates collected over 40 years from Australia and 3 other continents identified 66 MLVA types (MTs), including 6 predominant MTs. Typing of genes encoding acellular vaccine antigens showed changes that may be vaccine driven in 2 MTs prevalent in Australia.

Rapid and accurate typing of Bordetella pertussis targeting genes encoding acellular vaccine antigens using real time PCR and High Resolution Melt analysis.

Real Time-PCR (RT-PCR) and high resolution melt (HRM) analyses were used for rapid typing of genes encoding components of the pertussis acellular vaccine, namely prn, ptxA, fhaB, fim2 and fim3. The length polymorphisms in prn were detected by RT-PCR followed by HRM; single nucleotide polymorphisms in prn and other genes were detected by hairpin primer RT-PCR. These rapid methods are suitable for large-scale studies of vaccine-driven evolution of Bordetella pertussis.

Genome-wide analysis of single nucleotide polymorphisms in Bordetella pertussis using comparative genomic sequencing.

Bordetella pertussis is known to be a genotypically homogeneous pathogen but the extent of homogeneity at the genomic level is unknown. A currently circulating B. pertussis isolate from Australia was compared with the genome-sequenced Tohama I strain isolated in Japan in the 1950s from a distantly related lineage. Microarray-based comparative genome sequencing (CGS) was used to detect single nucleotide polymorphisms (SNPs) in a total of 1.4 Mb of the 4.09 Mb genome, including 1012 coding-regions, 217 pseudogenes and 268 intergenic regions. The CGS analysis, followed by validation using real-time PCR and DNA sequencing, identified 70 SNPs and five 1-3 bp indels, giving an overall frequency of base changes of 1 per 20 kb. Thirty-two of the 56 SNPs in coding regions were non-synonymous, including five located in virulence-associated genes. The data also allowed us to compare genomic diversity with other "clonal" human pathogens such as Mycobacterium tuberculosis and Yersinia pestis, showing that B. pertussis may be one of the least variable pathogenic bacterial species.