Characterization and genomic analysis of the Lyme disease spirochete bacteriophage ϕBB-1.

Lyme disease is a tick-borne infection caused by the spirochete Borrelia (Borreliella) burgdorferi. Borrelia species have highly fragmented genomes composed of a linear chromosome and a constellation of linear and circular plasmids some of which are required throughout the enzootic cycle. Included in this plasmid repertoire by almost all Lyme disease spirochetes are the 32-kb circular plasmid cp32 prophages that are capable of lytic replication to produce infectious virions called ϕBB-1. While the B. burgdorferi genome contains evidence of horizontal transfer, the mechanisms of gene transfer between strains remain unclear. While we know that ϕBB-1 transduces cp32 and shuttle vector DNA during in vitro cultivation, the extent of ϕBB-1 DNA transfer is not clear. Herein, we use proteomics and long-read sequencing to further characterize ϕBB-1 virions. Our studies identified the cp32 pac region and revealed that ϕBB-1 packages linear cp32s via a headful mechanism with preferential packaging of plasmids containing the cp32 pac region. Additionally, we find ϕBB-1 packages fragments of the linear chromosome and full-length plasmids including lp54, cp26, and others. Furthermore, sequencing of ϕBB-1 packaged DNA allowed us to resolve the covalently closed hairpin telomeres for the linear B. burgdorferi chromosome and most linear plasmids in strain CA-11.2A. Collectively, our results shed light on the biology of the ubiquitous ϕBB-1 phage and further implicates ϕBB-1 in the generalized transduction of diverse genes and the maintenance of genetic diversity in Lyme disease spirochetes.

Polyamines and linear DNA mediate bacterial threat assessment of bacteriophage infection.

Monitoring the extracellular environment for danger signals is a critical aspect of cellular survival. However, the danger signals released by dying bacteria and the mechanisms bacteria use for threat assessment remain largely unexplored. Here, we show that lysis of Pseudomonas aeruginosa cells releases polyamines that are subsequently taken up by surviving cells via a mechanism that relies on Gac/Rsm signaling. While intracellular polyamines spike in surviving cells, the duration of this spike varies according to the infection status of the cell. In bacteriophage-infected cells, intracellular polyamines are maintained at high levels, which inhibits replication of the bacteriophage genome. Many bacteriophages package linear DNA genomes and linear DNA is sufficient to trigger intracellular polyamine accumulation, suggesting that linear DNA is sensed as a second danger signal. Collectively, these results demonstrate how polyamines released by dying cells together with linear DNA allow P. aeruginosa to make threat assessments of cellular injury.

Evolutionary dynamics and structural consequences of de novo beneficial mutations and mutant lineages arising in a constant environment.

BACKGROUND: Microbial evolution experiments can be used to study the tempo and dynamics of evolutionary change in asexual populations, founded from single clones and growing into large populations with multiple clonal lineages. High-throughput sequencing can be used to catalog de novo mutations as potential targets of selection, determine in which lineages they arise, and track the fates of those lineages. Here, we describe a long-term experimental evolution study to identify targets of selection and to determine when, where, and how often those targets are hit. RESULTS: We experimentally evolved replicate Escherichia coli populations that originated from a mutator/nonsense suppressor ancestor under glucose limitation for between 300 and 500 generations. Whole-genome, whole-population sequencing enabled us to catalog 3346 de novo mutations that reached > 1% frequency. We sequenced the genomes of 96 clones from each population when allelic diversity was greatest in order to establish whether mutations were in the same or different lineages and to depict lineage dynamics. Operon-specific mutations that enhance glucose uptake were the first to rise to high frequency, followed by global regulatory mutations. Mutations related to energy conservation, membrane biogenesis, and mitigating the impact of nonsense mutations, both ancestral and derived, arose later. New alleles were confined to relatively few loci, with many instances of identical mutations arising independently in multiple lineages, among and within replicate populations. However, most never exceeded 10% in frequency and were at a lower frequency at the end of the experiment than at their maxima, indicating clonal interference. Many alleles mapped to key structures within the proteins that they mutated, providing insight into their functional consequences. CONCLUSIONS: Overall, we find that when mutational input is increased by an ancestral defect in DNA repair, the spectrum of high-frequency beneficial mutations in a simple, constant resource-limited environment is narrow, resulting in extreme parallelism where many adaptive mutations arise but few ever go to fixation.

Fitness and Productivity Increase with Ecotypic Diversity among Escherichia coli Strains That Coevolved in a Simple, Constant Environment.

The productivity of a biological community often correlates with its diversity. In the microbial world this phenomenon can sometimes be explained by positive, density-dependent interactions such as cross-feeding and syntrophy. These metabolic interactions help account for the astonishing variety of microbial life and drive many of the biogeochemical cycles without which life as we know it could not exist. While it is difficult to recapitulate experimentally how these interactions evolved among multiple taxa, we can explore in the laboratory how they arise within one. These experiments provide insight into how different bacterial ecotypes evolve and from these, possibly new "species." We have previously shown that in a simple, constant environment a single clone of Escherichia coli can give rise to a consortium of genetically and phenotypically differentiated strains, in effect, a set of ecotypes, that coexist by cross-feeding. We marked these different ecotypes and their shared ancestor by integrating fluorescent protein into their genomes and then used flow cytometry to show that each evolved strain is more fit than the shared ancestor, that pairs of evolved strains are fitter still, and that the entire consortium is the fittest of all. We further demonstrate that the rank order of fitness values agrees with estimates of yield, indicating that an experimentally evolved consortium more efficiently converts primary and secondary resources to offspring than its ancestor or any member acting in isolation.IMPORTANCE Polymicrobial consortia occur in both environmental and clinical settings. In many cases, diversity and productivity correlate in these consortia, especially when sustained by positive, density-dependent interactions. However, the evolutionary history of such entities is typically obscure, making it difficult to establish the relative fitness of consortium partners and to use those data to illuminate the diversity-productivity relationship. Here, we dissect an Escherichia coli consortium that evolved under continuous glucose limitation in the laboratory from a single common ancestor. We show that a partnership consisting of cross-feeding ecotypes is better able to secure primary and secondary resources and to convert those resources to offspring than the ancestral clone. Such interactions may be a prelude to a special form of syntrophy and are likely determinants of microbial community structure in nature, including those having clinical significance such as chronic infections.

Stability of Cross-Feeding Polymorphisms in Microbial Communities.

Cross-feeding, a relationship wherein one organism consumes metabolites excreted by another, is a ubiquitous feature of natural and clinically-relevant microbial communities and could be a key factor promoting diversity in extreme and/or nutrient-poor environments. However, it remains unclear how readily cross-feeding interactions form, and therefore our ability to predict their emergence is limited. In this paper we developed a mathematical model parameterized using data from the biochemistry and ecology of an E. coli cross-feeding laboratory system. The model accurately captures short-term dynamics of the two competitors that have been observed empirically and we use it to systematically explore the stability of cross-feeding interactions for a range of environmental conditions. We find that our simple system can display complex dynamics including multi-stable behavior separated by a critical point. Therefore whether cross-feeding interactions form depends on the complex interplay between density and frequency of the competitors as well as on the concentration of resources in the environment. Moreover, we find that subtly different environmental conditions can lead to dramatically different results regarding the establishment of cross-feeding, which could explain the apparently unpredictable between-population differences in experimental outcomes. We argue that mathematical models are essential tools for disentangling the complexities of cross-feeding interactions.

Ex uno plures: clonal reinforcement drives evolution of a simple microbial community.

A major goal of genetics is to define the relationship between phenotype and genotype, while a major goal of ecology is to identify the rules that govern community assembly. Achieving these goals by analyzing natural systems can be difficult, as selective pressures create dynamic fitness landscapes that vary in both space and time. Laboratory experimental evolution offers the benefit of controlling variables that shape fitness landscapes, helping to achieve both goals. We previously showed that a clonal population of E. coli experimentally evolved under continuous glucose limitation gives rise to a genetically diverse community consisting of one clone, CV103, that best scavenges but incompletely utilizes the limiting resource, and others, CV101 and CV116, that consume its overflow metabolites. Because this community can be disassembled and reassembled, and involves cooperative interactions that are stable over time, its genetic diversity is sustained by clonal reinforcement rather than by clonal interference. To understand the genetic factors that produce this outcome, and to illuminate the community's underlying physiology, we sequenced the genomes of ancestral and evolved clones. We identified ancestral mutations in intermediary metabolism that may have predisposed the evolution of metabolic interdependence. Phylogenetic reconstruction indicates that the lineages that gave rise to this community diverged early, as CV103 shares only one Single Nucleotide Polymorphism with the other evolved clones. Underlying CV103's phenotype we identified a set of mutations that likely enhance glucose scavenging and maintain redox balance, but may do so at the expense of carbon excreted in overflow metabolites. Because these overflow metabolites serve as growth substrates that are differentially accessible to the other community members, and because the scavenging lineage shares only one SNP with these other clones, we conclude that this lineage likely served as an "engine" generating diversity by creating new metabolic niches, but not the occupants themselves.

Characterization and genomic analysis of the Lyme disease spirochete bacteriophage ϕBB-1.

Lyme disease is a tick-borne infection caused by the spirochete Borrelia (Borreliella) burgdorferi. Borrelia species have highly fragmented genomes composed of a linear chromosome and a constellation of linear and circular plasmids some of which are required throughout the enzootic cycle. Included in this plasmid repertoire by almost all Lyme disease spirochetes are the 32-kb circular plasmid cp32 prophages that are capable of lytic replication to produce infectious virions called ϕBB-1. While the B. burgdorferi genome contains evidence of horizontal transfer, the mechanisms of gene transfer between strains remain unclear. While we know that ϕBB-1 transduces cp32 and shuttle vector DNA during in vitro cultivation, the extent of ϕBB-1 DNA transfer is not clear. Herein, we use proteomics and long-read sequencing to further characterize ϕBB-1 virions. Our studies identified the cp32 pac region and revealed that ϕBB-1 packages linear cp32s via a headful mechanism with preferentially packaging of plasmids containing the cp32 pac region. Additionally, we find ϕBB-1 packages fragments of the linear chromosome and full-length plasmids including lp54, cp26, and others. Furthermore, sequencing of ϕBB-1 packaged DNA allowed us to resolve the covalently closed hairpin telomeres for the linear B. burgdorferi chromosome and most linear plasmids in strain CA-11.2A. Collectively, our results shed light on the biology of the ubiquitous ϕBB-1 phage and further implicates ϕBB-1 in the generalized transduction of diverse genes and the maintenance of genetic diversity in Lyme disease spirochetes.

Longitudinal map of transcriptome changes in the Lyme pathogen Borrelia burgdorferi during tick-borne transmission.

Borrelia burgdorferi (Bb), the causative agent of Lyme disease, adapts to vastly different environments as it cycles between tick vector and vertebrate host. During a tick bloodmeal, Bb alters its gene expression to prepare for vertebrate infection; however, the full range of transcriptional changes that occur over several days inside of the tick are technically challenging to capture. We developed an experimental approach to enrich Bb cells to longitudinally define their global transcriptomic landscape inside nymphal Ixodes scapularis ticks during a transmitting bloodmeal. We identified 192 Bb genes that substantially change expression over the course of the bloodmeal from 1 to 4 days after host attachment. The majority of upregulated genes encode proteins found at the cell envelope or proteins of unknown function, including 45 outer surface lipoproteins embedded in the unusual protein-rich coat of Bb. As these proteins may facilitate Bb interactions with the host, we utilized mass spectrometry to identify candidate tick proteins that physically associate with Bb. The Bb enrichment methodology along with the ex vivo Bb transcriptomes and candidate tick interacting proteins presented here provide a resource to facilitate investigations into key determinants of Bb priming and transmission during the tick stage of its unique transmission cycle.

Coupled induction of prophage and virulence factors during tick transmission of the Lyme disease spirochete.

The alternative sigma factor RpoS plays a central role in the critical host-adaptive response of the Lyme disease spirochete, Borrelia burgdorferi. We previously identified bbd18 as a negative regulator of RpoS but could not inactivate bbd18 in wild-type spirochetes. In the current study we employed an inducible bbd18 gene to demonstrate the essential nature of BBD18 for viability of wild-type spirochetes in vitro and at a unique point in vivo. Transcriptomic analyses of BBD18-depleted cells demonstrated global induction of RpoS-dependent genes prior to lysis, with the absolute requirement for BBD18, both in vitro and in vivo, circumvented by deletion of rpoS. The increased expression of plasmid prophage genes and the presence of phage particles in the supernatants of lysing cultures indicate that RpoS regulates phage lysis-lysogeny decisions. Through this work we identify a mechanistic link between endogenous prophages and the RpoS-dependent adaptive response of the Lyme disease spirochete.

E Unibus Plurum: genomic analysis of an experimentally evolved polymorphism in Escherichia coli.

Microbial populations founded by a single clone and propagated under resource limitation can become polymorphic. We sought to elucidate genetic mechanisms whereby a polymorphism evolved in Escherichia coli under glucose limitation and persisted because of cross-feeding among multiple adaptive clones. Apart from a 29 kb deletion in the dominant clone, no large-scale genomic changes distinguished evolved clones from their common ancestor. Using transcriptional profiling on co-evolved clones cultured separately under glucose-limitation we identified 180 genes significantly altered in expression relative to the common ancestor grown under similar conditions. Ninety of these were similarly expressed in all clones, and many of the genes affected (e.g., mglBAC, mglD, and lamB) are in operons coordinately regulated by CRP and/or rpoS. While the remaining significant expression differences were clone-specific, 93% were exhibited by the majority clone, many of which are controlled by global regulators, CRP and CpxR. When transcriptional profiling was performed on adaptive clones cultured together, many expression differences that distinguished the majority clone cultured in isolation were absent, suggesting that CpxR may be activated by overflow metabolites removed by cross-feeding strains in co-culture. Relative to their common ancestor, shared expression differences among adaptive clones were partly attributable to early-arising shared mutations in the trans-acting global regulator, rpoS, and the cis-acting regulator, mglO. Gene expression differences that distinguished clones may in part be explained by mutations in trans-acting regulators malT and glpK, and in cis-acting sequences of acs. In the founder, a cis-regulatory mutation in acs (acetyl CoA synthetase) and a structural mutation in glpR (glycerol-3-phosphate repressor) likely favored evolution of specialists that thrive on overflow metabolites. Later-arising mutations that led to specialization emphasize the importance of compensatory rather than gain-of-function mutations in this system. Taken together, these findings underscore the importance of regulatory change, founder genotype, and the biotic environment in the adaptive evolution of microbes.

Complete Genome Sequence of the N4-like Pseudomonas aeruginosa Bacteriophage vB_PaeP_CMS1.

vB_PaeP_CMS1 is a lytic bacteriophage that infects Pseudomonas aeruginosa. Whole-genome sequencing revealed that vB_PaeP_CMS1 has a linear double-stranded DNA (dsDNA) genome composed of 72,673 bp, with a GC content of 54.9%, that harbors 92 predicted open reading frames, with the greatest genome homology to members of the family Podoviridae, genus Litunavirus.