Rapid dissemination of host metabolism-manipulating genes via integrative and conjugative elements.

Integrative and conjugative elements (ICEs) are self-transmissible mobile elements that transfer functional genetic units across broad phylogenetic distances. Accessory genes shuttled by ICEs can make significant contributions to bacterial fitness. Most ICEs characterized to date encode readily observable phenotypes contributing to symbiosis, pathogenicity, and antimicrobial resistance, yet the majority of ICEs carry genes of unknown function. Recent observations of rapid acquisition of ICEs in a pandemic lineage of Pseudomonas syringae pv. actinidae led to investigation of the structural and functional diversity of these elements. Fifty-three unique ICE types were identified across the P. syringae species complex. Together they form a distinct family of ICEs (PsICEs) that share a distant relationship to ICEs found in Pseudomonas aeruginosa. PsICEs are defined by conserved backbone genes punctuated by an array of accessory cargo genes, are highly recombinogenic, and display distinct evolutionary histories compared to their bacterial hosts. The most common cargo is a recently disseminated 16-kb mobile genetic element designated Tn6212. Deletion of Tn6212 did not alter pathogen growth in planta, but mutants displayed fitness defects when grown on tricarboxylic acid (TCA) cycle intermediates. RNA-seq analysis of a set of nested deletion mutants showed that a Tn6212-encoded LysR regulator has global effects on chromosomal gene expression. We show that Tn6212 responds to preferred carbon sources and manipulates bacterial metabolism to maximize growth.

Ancient Darwinian replicators nested within eubacterial genomes.

Integrative mobile genetic elements (MGEs), such as transposons and insertion sequences, propagate within bacterial genomes, but persistence times in individual lineages are short. For long-term survival, MGEs must continuously invade new hosts by horizontal transfer. Theoretically, MGEs that persist for millions of years in single lineages, and are thus subject to vertical inheritance, should not exist. Here we draw attention to an exception - a class of MGE termed REPIN. REPINs are non-autonomous MGEs whose duplication depends on non-jumping RAYT transposases. Comparisons of REPINs and typical MGEs show that replication rates of REPINs are orders of magnitude lower, REPIN population size fluctuations correlate with changes in available genome space, REPIN conservation depends on RAYT function, and REPIN diversity accumulates within host lineages. These data lead to the hypothesis that REPINs form enduring, beneficial associations with eubacterial chromosomes. Given replicative nesting, our hypothesis predicts conflicts arising from the diverging effects of selection acting simultaneously on REPINs and host genomes. Evidence in support comes from patterns of REPIN abundance and diversity in two distantly related bacterial species. Together this bolsters the conclusion that REPINs are the genetic counterpart of mutualistic endosymbiotic bacteria.

Barcoding Populations of Pseudomonas fluorescens SBW25.

In recent years, evolutionary biologists have developed an increasing interest in the use of barcoding strategies to study eco-evolutionary dynamics of lineages within evolving populations and communities. Although barcoded populations can deliver unprecedented insight into evolutionary change, barcoding microbes presents specific technical challenges. Here, strategies are described for barcoding populations of the model bacterium Pseudomonas fluorescens SBW25, including the design and cloning of barcoded regions, preparation of libraries for amplicon sequencing, and quantification of resulting barcoded lineages. In so doing, we hope to aid the design and implementation of barcoding methodologies in a broad range of model and non-model organisms.

The effect of bottleneck size on evolution in nested Darwinian populations.

Previous work has shown how a minimal ecological structure consisting of patchily distributed resources and recurrent dispersal between patches can scaffold Darwinian properties onto collections of cells. When the timescale of dispersal is long compared with the time to consume resources, patch fitness increases but comes at a cost to cell growth rates. This creates conditions that initiate evolutionary transitions in individuality. A key feature of the scaffold is a bottleneck created during dispersal, causing patches to be founded by single cells. The bottleneck decreases competition within patches and, hence, creates a strong hereditary link at the level of patches. Here, we construct a fully stochastic model to investigate the effect of bottleneck size on the evolutionary dynamics of both cells and collectives. We show that larger bottlenecks simply slow the dynamics, but, at some point, which depends on the parameters of the within-patch model, the direction of evolution towards the equilibrium reverses. Introduction of random fluctuations in bottleneck sizes with some positive probability of smaller sizes counteracts this, even when the probability of smaller bottlenecks is minimal.

Meta-population structure and the evolutionary transition to multicellularity.

The evolutionary transition to multicellularity has occurred on numerous occasions, but transitions to complex life forms are rare. Here, using experimental bacterial populations as proxies for nascent multicellular organisms, we manipulate ecological factors shaping the evolution of groups. Groups were propagated under regimes requiring reproduction via a life cycle replete with developmental and dispersal (propagule) phases, but in one treatment lineages never mixed, whereas in a second treatment, cells from different lineages experienced intense competition during the dispersal phase. The latter treatment favoured traits promoting cell growth at the expense of traits underlying group fitness - a finding that is supported by results from a mathematical model. Our results show that the transition to multicellularity benefits from ecological conditions that maintain discreteness not just of the group (soma) phase, but also of the dispersal (germline) phase.

Repeated Phenotypic Evolution by Different Genetic Routes in Pseudomonas fluorescens SBW25.

Repeated evolution of functionally similar phenotypes is observed throughout the tree of life. The extent to which the underlying genetics are conserved remains an area of considerable interest. Previously, we reported the evolution of colony switching in two independent lineages of Pseudomonas fluorescens SBW25. The phenotypic and genotypic bases of colony switching in the first lineage (Line 1) have been described elsewhere. Here, we deconstruct the evolution of colony switching in the second lineage (Line 6). We show that, as for Line 1, Line 6 colony switching results from an increase in the expression of a colanic acid-like polymer (CAP). At the genetic level, nine mutations occur in Line 6. Only one of these-a nonsynonymous point mutation in the housekeeping sigma factor rpoD-is required for colony switching. In contrast, the genetic basis of colony switching in Line 1 is a mutation in the metabolic gene carB. A molecular model has recently been proposed whereby the carB mutation increases capsulation by redressing the intracellular balance of positive (ribosomes) and negative (RsmAE/CsrA) regulators of a positive feedback loop in capsule expression. We show that Line 6 colony switching is consistent with this model; the rpoD mutation generates an increase in ribosomal gene expression, and ultimately an increase in CAP expression.

Ribosome Provisioning Activates a Bistable Switch Coupled to Fast Exit from Stationary Phase.

Observations of bacteria at the single-cell level have revealed many instances of phenotypic heterogeneity within otherwise clonal populations, but the selective causes, molecular bases, and broader ecological relevance remain poorly understood. In an earlier experiment in which the bacterium Pseudomonas fluorescens SBW25 was propagated under a selective regime that mimicked the host immune response, a genotype evolved that stochastically switched between capsulation states. The genetic cause was a mutation in carB that decreased the pyrimidine pool (and growth rate), lowering the activation threshold of a preexisting but hitherto unrecognized phenotypic switch. Genetic components surrounding bifurcation of UTP flux toward DNA/RNA or UDP-glucose (a precursor of colanic acid forming the capsules) were implicated as key components. Extending these molecular analyses-and based on a combination of genetics, transcriptomics, biochemistry, and mathematical modeling-we show that pyrimidine limitation triggers an increase in ribosome biosynthesis and that switching is caused by competition between ribosomes and CsrA/RsmA proteins for the mRNA transcript of a positively autoregulated activator of colanic acid biosynthesis. We additionally show that in the ancestral bacterium the switch is part of a program that determines stochastic entry into a semiquiescent capsulated state, ensures that such cells are provisioned with excess ribosomes, and enables provisioned cells to exit rapidly from stationary phase under permissive conditions.

Life cycles, fitness decoupling and the evolution of multicellularity.

Cooperation is central to the emergence of multicellular life; however, the means by which the earliest collectives (groups of cells) maintained integrity in the face of destructive cheating types is unclear. One idea posits cheats as a primitive germ line in a life cycle that facilitates collective reproduction. Here we describe an experiment in which simple cooperating lineages of bacteria were propagated under a selective regime that rewarded collective-level persistence. Collectives reproduced via life cycles that either embraced, or purged, cheating types. When embraced, the life cycle alternated between phenotypic states. Selection fostered inception of a developmental switch that underpinned the emergence of collectives whose fitness, during the course of evolution, became decoupled from the fitness of constituent cells. Such development and decoupling did not occur when groups reproduced via a cheat-purging regime. Our findings capture key events in the evolution of Darwinian individuality during the transition from single cells to multicellularity.

Anaerobically Grown Escherichia coli Has an Enhanced Mutation Rate and Distinct Mutational Spectra.

Oxidative stress is a major cause of mutation but little is known about how growth in the absence of oxygen impacts the rate and spectrum of mutations. We employed long-term mutation accumulation experiments to directly measure the rates and spectra of spontaneous mutation events in Escherichia coli populations propagated under aerobic and anaerobic conditions. To detect mutations, whole genome sequencing was coupled with methods of analysis sufficient to identify a broad range of mutational classes, including structural variants (SVs) generated by movement of repetitive elements. The anaerobically grown populations displayed a mutation rate nearly twice that of the aerobic populations, showed distinct asymmetric mutational strand biases, and greater insertion element activity. Consistent with mutation rate and spectra observations, genes for transposition and recombination repair associated with SVs were up-regulated during anaerobic growth. Together, these results define differences in mutational spectra affecting the evolution of facultative anaerobes.

Causes and Biophysical Consequences of Cellulose Production by Pseudomonas fluorescens SBW25 at the Air-Liquid Interface.

Cellulose-overproducing wrinkly spreader mutants of Pseudomonas fluorescens SBW25 have been the focus of much investigation, but conditions promoting the production of cellulose in ancestral strain SBW25 and its effects and consequences have escaped in-depth investigation through lack of an in vitro phenotype. Here, using a custom-built device, we reveal that in static broth microcosms, ancestral SBW25 encounters environmental signals at the air-liquid interface that activate, via three diguanylate cyclase-encoding pathways (Wsp, Aws, and Mws), production of cellulose. Secretion of the polymer at the meniscus leads to modification of the environment and growth of numerous microcolonies that extend from the surface. Accumulation of cellulose and associated microbial growth leads to Rayleigh-Taylor instability resulting in bioconvection and rapid transport of water-soluble products over tens of millimeters. Drawing upon data, we built a mathematical model that recapitulates experimental results and captures the interactions between biological, chemical and physical processes.IMPORTANCE This work reveals a hitherto unrecognized behavior that manifests at the air-liquid interface that depends on production of cellulose and hints at undiscovered dimensions to bacterial life at surfaces. Additionally, the study links activation of known diguanylate cyclase-encoding pathways to cellulose expression and to signals encountered at the meniscus. Further significance stems from recognition of the consequences of fluid instabilities arising from surface production of cellulose for transport of water-soluble products over large distances.

Unravelling the complexity and redundancy of carbon catabolic repression in Pseudomonas fluorescens SBW25.

The two-component system CbrAB is the principal regulator for cellular metabolic balance in Pseudomonas fluorescens SBW25 and is necessary for growth on many substrates including xylose. To understand the regulatory linkage between CbrAB and genes for xylose utilization (xut), we performed transposon mutagenesis of ΔcbrB to select for Xut+ suppressors. This led to identification of crc and hfq. Subsequent genetic and biochemical analysis showed that Crc and Hfq are key mediators of succinate-provoked carbon catabolite repression (CCR). Specifically, Crc/Hfq sequentially bind to mRNAs of both the transcriptional activator and structural genes involved in xylose catabolism. However, in the absence of succinate, repression is relieved through competitive binding by two ncRNAs, CrcY and CrcZ, whose expression is activated by CbrAB. These findings provoke a model for CCR in which it is assumed that crc and hfq are functionally complementary, whereas crcY and crcZ are genetically redundant. Inactivation of either crcY or crcZ produced no effects on bacterial fitness in laboratory media, however, results of mathematical modelling predict that the co-existence of crcY and crcZ requires separate functional identity. Finally, we provide empirical evidence that CCR is advantageous in nutrient-complex environments where preferred carbon sources are present at high concentrations but fluctuate in their availability.

The ecological genetics of Pseudomonas syringae from kiwifruit leaves.

Interactions between commensal microbes and invading pathogens are understudied, despite their likely effects on pathogen population structure and infection processes. We describe the population structure and genetic diversity of a broad range of co-occurring Pseudomonas syringae isolated from infected and uninfected kiwifruit during an outbreak of bleeding canker disease caused by P. syringae pv. actinidiae (Psa) in New Zealand. Overall population structure was clonal and affected by ecological factors including infection status and cultivar. Most isolates are members of a new clade in phylogroup 3 (PG3a), also present on kiwifruit leaves in China and Japan. Stability of the polymorphism between pathogenic Psa and commensal P. syringae PG3a isolated from the same leaf was tested using reciprocal invasion from rare assays in vitro and in planta. P. syringae G33C (PG3a) inhibited Psa NZ54, while the presence of Psa NZ54 enhanced the growth of P. syringae G33C. This effect could not be attributed to virulence activity encoded by the Type 3 secretion system of Psa. Together our data contribute toward the development of an ecological perspective on the genetic structure of pathogen populations.

Bistability in a metabolic network underpins the de novo evolution of colony switching in Pseudomonas fluorescens.

Phenotype switching is commonly observed in nature. This prevalence has allowed the elucidation of a number of underlying molecular mechanisms. However, little is known about how phenotypic switches arise and function in their early evolutionary stages. The first opportunity to provide empirical insight was delivered by an experiment in which populations of the bacterium Pseudomonas fluorescens SBW25 evolved, de novo, the ability to switch between two colony phenotypes. Here we unravel the molecular mechanism behind colony switching, revealing how a single nucleotide change in a gene enmeshed in central metabolism (carB) generates such a striking phenotype. We show that colony switching is underpinned by ON/OFF expression of capsules consisting of a colanic acid-like polymer. We use molecular genetics, biochemical analyses, and experimental evolution to establish that capsule switching results from perturbation of the pyrimidine biosynthetic pathway. Of central importance is a bifurcation point at which uracil triphosphate is partitioned towards either nucleotide metabolism or polymer production. This bifurcation marks a cell-fate decision point whereby cells with relatively high pyrimidine levels favour nucleotide metabolism (capsule OFF), while cells with lower pyrimidine levels divert resources towards polymer biosynthesis (capsule ON). This decision point is present and functional in the wild-type strain. Finally, we present a simple mathematical model demonstrating that the molecular components of the decision point are capable of producing switching. Despite its simple mutational cause, the connection between genotype and phenotype is complex and multidimensional, offering a rare glimpse of how noise in regulatory networks can provide opportunity for evolution.

Fragmentation modes and the evolution of life cycles.

Reproduction is a defining feature of living systems. To reproduce, aggregates of biological units (e.g., multicellular organisms or colonial bacteria) must fragment into smaller parts. Fragmentation modes in nature range from binary fission in bacteria to collective-level fragmentation and the production of unicellular propagules in multicellular organisms. Despite this apparent ubiquity, the adaptive significance of fragmentation modes has received little attention. Here, we develop a model in which groups arise from the division of single cells that do not separate but stay together until the moment of group fragmentation. We allow for all possible fragmentation patterns and calculate the population growth rate of each associated life cycle. Fragmentation modes that maximise growth rate comprise a restrictive set of patterns that include production of unicellular propagules and division into two similar size groups. Life cycles marked by single-cell bottlenecks maximise population growth rate under a wide range of conditions. This surprising result offers a new evolutionary explanation for the widespread occurrence of this mode of reproduction. All in all, our model provides a framework for exploring the adaptive significance of fragmentation modes and their associated life cycles.

Evolution of copper resistance in the kiwifruit pathogen Pseudomonas syringae pv. actinidiae through acquisition of integrative conjugative elements and plasmids.

Horizontal gene transfer can precipitate rapid evolutionary change. In 2010 the global pandemic of kiwifruit canker disease caused by Pseudomonas syringae pv. actinidiae (Psa) reached New Zealand. At the time of introduction, the single clone responsible for the outbreak was sensitive to copper, however, analysis of a sample of isolates taken in 2015 and 2016 showed that a quarter were copper resistant. Genome sequences of seven strains showed that copper resistance - comprising czc/cusABC and copABCD systems - along with resistance to arsenic and cadmium, was acquired via uptake of integrative conjugative elements (ICEs), but also plasmids. Comparative analysis showed ICEs to have a mosaic structure, with one being a tripartite arrangement of two different ICEs and a plasmid that were isolated in 1921 (USA), 1968 (NZ) and 1988 (Japan), from P. syringae pathogens of millet, wheat and kiwifruit respectively. Two of the Psa ICEs were nearly identical to two ICEs isolated from kiwifruit leaf colonists prior to the introduction of Psa into NZ. Additionally, we show ICE transfer in vitro and in planta, analyze fitness consequences of ICE carriage, capture the de novo formation of novel recombinant ICEs, and explore ICE host-range.

Molecular mechanisms of xylose utilization by Pseudomonas fluorescens: overlapping genetic responses to xylose, xylulose, ribose and mannitol.

Bacterial degradation of xylose is sequentially mediated by two enzymes - an isomerase (XutA) and a xylulokinase (XutB) - with xylulose as an intermediate. Pseudomonas fluorescens SBW25, though capable of growth on xylose as a sole carbon source, encodes only one degradative enzyme XutA at the xylose utilization (xut) locus. Here, using site-directed mutagenesis and transcriptional assays, we have identified two functional xylulokinase-encoding genes (xutB1 and xutB2) and further show that expression of xutB1 is specifically induced by xylose. Surprisingly, xylose-induced xutB1 expression is mediated by the mannitol-responsive regulator MtlR, using xylulose rather than xylose as the direct inducer. In contrast, expression of the xutA operon is regulated by XutR - a transcriptional activator of the AraC family - in a xylose-, xylulose- and ribose-dependent manner. Detailed genetic and biochemical analyses of XutR, including DNase I footprinting assays, suggest an unconventional model of XutR regulation that does not involve DNA-looping, a mechanism typically found for AraC-type regulators from enteric bacteria. XutR functions as a dimer and recognizes two inverted repeat sequences, but binding to one half site is weak thus requiring an inducer molecule such as xylose for activation.

Founder niche constrains evolutionary adaptive radiation.

Adaptive radiation of a lineage into a range of organisms with different niches underpins the evolution of life's diversity. Although the role of the environment in shaping adaptive radiation is well established, theory predicts that the evolvability and niche of the founding ancestor are also of importance. Direct demonstration of a causal link requires resolving the independent effects of these additional factors. Here, we accomplish this using experimental bacterial populations and demonstrate how the dynamics of adaptive radiation are constrained by the niche of the founder. We manipulated the propensity of the founder to undergo adaptive radiation and resolved the underlying causal changes in both its evolvability and niche. Evolvability did not change, but the propensity for adaptive radiation was altered by changes in the position and breadth of the niche of the founder. These observations provide direct empirical evidence for a link between the niche of organisms and their propensity for adaptive radiation. This general mechanism may have rendered the evolutionary dynamics of extant adaptive radiations dependent on chance events that determined their founding ancestors.

Role of the Transporter-Like Sensor Kinase CbrA in Histidine Uptake and Signal Transduction.

UNLABELLED: CbrA is an atypical sensor kinase found in Pseudomonas. The autokinase domain is connected to a putative transporter of the sodium/solute symporter family (SSSF). CbrA functions together with its cognate response regulator, CbrB, and plays an important role in nutrient acquisition, including regulation of hut genes for the utilization of histidine and its derivative, urocanate. Here we report on the findings of a genetic and biochemical analysis of CbrA with a focus on the function of the putative transporter domain. The work was initiated with mutagenesis of histidine uptake-proficient strains to identify histidine-specific transport genes located outside the hut operon. Genes encoding transporters were not identified, but mutations were repeatedly found in cbrA. This, coupled with the findings of [(3)H]histidine transport assays and further mutagenesis, implicated CbrA in histidine uptake. In addition, mutations in different regions of the SSSF domain abolished signal transduction. Site-specific mutations were made at four conserved residues: W55 and G172 (SSSF domain), H766 (H box), and N876 (N box). The mutations W55G, G172H, and N876G compromised histidine transport but had minimal effects on signal transduction. The H766G mutation abolished both transport and signal transduction, but the capacity to transport histidine was restored upon complementation with a transport-defective allele of CbrA, most likely due to interdomain interactions. Our combined data implicate the SSSF domain of CbrA in histidine transport and suggest that transport is coupled to signal transduction. IMPORTANCE: Nutrient acquisition in bacteria typically involves membrane-bound sensors that, via cognate response regulators, determine the activity of specific transporters. However, nutrient perception and uptake are often coupled processes. Thus, from a physiological perspective, it would make sense for systems that couple the process of signaling and transport within a single protein and where transport is itself the stimulus that precipitates signal transduction to have evolved. The CbrA regulator in Pseudomonas represents a unique type of sensor kinase whose autokinase domain is connected to a transporter domain. We present genetic and biochemical evidence that suggests that CbrA plays a dual role in histidine uptake and sensing and that transport is dependent on signal transduction.

Automated reconstruction of whole-genome phylogenies from short-sequence reads.

Studies of microbial evolutionary dynamics are being transformed by the availability of affordable high-throughput sequencing technologies, which allow whole-genome sequencing of hundreds of related taxa in a single study. Reconstructing a phylogenetic tree of these taxa is generally a crucial step in any evolutionary analysis. Instead of constructing genome assemblies for all taxa, annotating these assemblies, and aligning orthologous genes, many recent studies 1) directly map raw sequencing reads to a single reference sequence, 2) extract single nucleotide polymorphisms (SNPs), and 3) infer the phylogenetic tree using maximum likelihood methods from the aligned SNP positions. However, here we show that, when using such methods to reconstruct phylogenies from sets of simulated sequences, both the exclusion of nonpolymorphic positions and the alignment to a single reference genome, introduce systematic biases and errors in phylogeny reconstruction. To address these problems, we developed a new method that combines alignments from mappings to multiple reference sequences and show that this successfully removes biases from the reconstructed phylogenies. We implemented this method as a web server named REALPHY (Reference sequence Alignment-based Phylogeny builder), which fully automates phylogenetic reconstruction from raw sequencing reads.

Genomic analysis of the Kiwifruit pathogen Pseudomonas syringae pv. actinidiae provides insight into the origins of an emergent plant disease.

The origins of crop diseases are linked to domestication of plants. Most crops were domesticated centuries--even millennia--ago, thus limiting opportunity to understand the concomitant emergence of disease. Kiwifruit (Actinidia spp.) is an exception: domestication began in the 1930s with outbreaks of canker disease caused by P. syringae pv. actinidiae (Psa) first recorded in the 1980s. Based on SNP analyses of two circularized and 34 draft genomes, we show that Psa is comprised of distinct clades exhibiting negligible within-clade diversity, consistent with disease arising by independent samplings from a source population. Three clades correspond to their geographical source of isolation; a fourth, encompassing the Psa-V lineage responsible for the 2008 outbreak, is now globally distributed. Psa has an overall clonal population structure, however, genomes carry a marked signature of within-pathovar recombination. SNP analysis of Psa-V reveals hundreds of polymorphisms; however, most reside within PPHGI-1-like conjugative elements whose evolution is unlinked to the core genome. Removal of SNPs due to recombination yields an uninformative (star-like) phylogeny consistent with diversification of Psa-V from a single clone within the last ten years. Growth assays provide evidence of cultivar specificity, with rapid systemic movement of Psa-V in Actinidia chinensis. Genomic comparisons show a dynamic genome with evidence of positive selection on type III effectors and other candidate virulence genes. Each clade has highly varied complements of accessory genes encoding effectors and toxins with evidence of gain and loss via multiple genetic routes. Genes with orthologs in vascular pathogens were found exclusively within Psa-V. Our analyses capture a pathogen in the early stages of emergence from a predicted source population associated with wild Actinidia species. In addition to candidate genes as targets for resistance breeding programs, our findings highlight the importance of the source population as a reservoir of new disease.