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.

The untapped potential of phage model systems as therapeutic agents

With the emergence of widespread antibiotic resistance, phages are an appealing alternative to antibiotics in the fight against multidrug-resistant bacteria. Over the past few years, many phages have been isolated from various environments to treat bacterial pathogens. While isolating novel phages for treatment has had some success for compassionate use, developing novel phages into a general therapeutic will require considerable time and financial resource investments. These investments may be less significant for well-established phage model systems. The knowledge acquired from decades of research on their structure, life cycle, and evolution ensures safe application and efficient handling. However, one major downside of the established phage model systems is their inability to infect pathogenic bacteria. This problem is not insurmountable; phage host range can be extended through genetic engineering or evolution experiments. In the future, breeding model phages to infect pathogens could provide a new avenue to develop phage therapeutic agents.

Stepwise evolution of E. coli C and ΦX174 reveals unexpected lipopolysaccharide (LPS) diversity

Phage therapy is a promising method for the treatment of multi-drug-resistant bacterial infections. However, its long-term efficacy depends on understanding the evolutionary effects of the treatment. Current knowledge of such evolutionary effects is lacking, even in well-studied systems. We used the bacterium Escherichia coli C and its bacteriophage ΦX174, which infects cells using host lipopolysaccharide (LPS) molecules. We first generated 31 bacterial mutants resistant to ΦX174 infection. Based on the genes disrupted by these mutations, we predicted that these E. coli C mutants collectively produce eight unique LPS structures. We then developed a series of evolution experiments to select for ΦX174 mutants capable of infecting the resistant strains. During phage adaptation, we distinguished two types of phage resistance: one that was easily overcome by ΦX174 with few mutational steps ("easy" resistance), and one that was more difficult to overcome ("hard" resistance). We found that increasing the diversity of the host and phage populations could accelerate the adaptation of phage ΦX174 to overcome the hard resistance phenotype. From these experiments, we isolated 16 ΦX174 mutants that, together, can infect all 31 initially resistant E. coli C mutants. Upon determining the infectivity profiles of these 16 evolved phages, we uncovered 14 distinct profiles. Given that only eight profiles are anticipated if the LPS predictions are correct, our findings highlight that the current understanding of LPS biology is insufficient to accurately forecast the evolutionary outcomes of bacterial populations infected by phage.

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.

Transposable elements promote the evolution of genome streamlining

Eukaryotes and prokaryotes have distinct genome architectures, with marked differences in genome size, the ratio of coding/non-coding DNA, and the abundance of transposable elements (TEs). As TEs replicate independently of their hosts, the proliferation of TEs is thought to have driven genome expansion in eukaryotes. However, prokaryotes also have TEs in intergenic spaces, so why do prokaryotes have small, streamlined genomes? Using an in silico model describing the genomes of single-celled asexual organisms that coevolve with TEs, we show that TEs acquired from the environment by horizontal gene transfer can promote the evolution of genome streamlining. The process depends on local interactions and is underpinned by rock-paper-scissors dynamics in which populations of cells with streamlined genomes beat TEs, which beat non-streamlined genomes, which beat streamlined genomes, in continuous and repeating cycles. Streamlining is maladaptive to individual cells, but improves lineage viability by hindering the proliferation of TEs. Streamlining does not evolve in sexually reproducing populations because recombination partially frees TEs from the deleterious effects they cause. This article is part of the theme issue 'The secret lives of microbial mobile genetic elements'.

How sequence populations persist inside bacterial genomes

Compared to their eukaryotic counterparts, bacterial genomes are small and contain extremely tightly packed genes. Repetitive sequences are rare but not completely absent. One of the most common repeat families is REPINs. REPINs can replicate in the host genome and form populations that persist for millions of years. Here, we model the interactions of these intragenomic sequence populations with the bacterial host. We first confirm well-established results, in the presence and absence of horizontal gene transfer (hgt) sequence populations either expand until they drive the host to extinction or the sequence population gets purged from the genome. We then show that a sequence population can be stably maintained, when each individual sequence provides a benefit that decreases with increasing sequence population size. Maintaining a sequence population of stable size also requires the replication of the sequence population to be costly to the host, otherwise the sequence population size will increase indefinitely. Surprisingly, in regimes with high hgt rates, the benefit conferred by the sequence population does not have to exceed the damage it causes to its host. Our analyses provide a plausible scenario for the persistence of sequence populations in bacterial genomes. We also hypothesize a limited biologically relevant parameter range for the provided benefit, which can be tested in future experiments.

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.

Parallel evolution of HIV-1 in a long-term experiment

One of the most intriguing puzzles in biology is the degree to which evolution is repeatable. The repeatability of evolution, or parallel evolution, has been studied in a variety of model systems, but has rarely been investigated with clinically relevant viruses. To investigate parallel evolution of HIV-1, we passaged two replicate HIV-1 populations for almost 1 year in each of two human T-cell lines. For each of the four evolution lines, we determined the genetic composition of the viral population at nine time points by deep sequencing the entire genome. Mutations that were carried by the majority of the viral population accumulated continuously over 1 year in each evolution line. Many majority mutations appeared in more than one evolution line, that is, our experiments showed an extreme degree of parallel evolution. In one of the evolution lines, 62% of the majority mutations also occur in another line. The parallelism impairs our ability to reconstruct the evolutionary history by phylogenetic methods. We show that one can infer the correct phylogenetic topology by including minority mutations in our analysis. We also find that mutation diversity at the beginning of the experiment is predictive of the frequency of majority mutations at the end of the experiment.

Dissecting HIV Virulence: Heritability of Setpoint Viral Load, CD4+ T-Cell Decline, and Per-Parasite Pathogenicity

Pathogen strains may differ in virulence because they attain different loads in their hosts, or because they induce different disease-causing mechanisms independent of their load. In evolutionary ecology, the latter is referred to as “per-parasite pathogenicity”. Using viral load and CD4+ T-cell measures from 2014 HIV-1 subtype B-infected individuals enrolled in the Swiss HIV Cohort Study, we investigated if virulence—measured as the rate of decline of CD4+ T cells—and per-parasite pathogenicity are heritable from donor to recipient. We estimated heritability by donor–recipient regressions applied to 196 previously identified transmission pairs, and by phylogenetic mixed models applied to a phylogenetic tree inferred from HIV pol sequences. Regressing the CD4+ T-cell declines and per-parasite pathogenicities of the transmission pairs did not yield heritability estimates significantly different from zero. With the phylogenetic mixed model, however, our best estimate for the heritability of the CD4+ T-cell decline is 17% (5–30%), and that of the per-parasite pathogenicity is 17% (4–29%). Further, we confirm that the set-point viral load is heritable, and estimate a heritability of 29% (12–46%). Interestingly, the pattern of evolution of all these traits differs significantly from neutrality, and is most consistent with stabilizing selection for the set-point viral load, and with directional selection for the CD4+ T-cell decline and the per-parasite pathogenicity. Our analysis shows that the viral genotype affects virulence mainly by modulating the per-parasite pathogenicity, while the indirect effect via the set-point viral load is minor.

Estimating the mutational fitness effects distribution during early HIV infection

The evolution of HIV during acute infection is often considered a neutral process. Recent analysis of sequencing data from this stage of infection, however, showed high levels of shared mutations between independent viral populations. This suggests that selection might play a role in the early stages of HIV infection. We adapted an existing model for random evolution during acute HIV-infection to include selection. Simulations of this model were used to fit a global mutational fitness effects distribution to previously published sequencing data of the env gene of individuals with acute HIV infection. Measures of sharing between viral populations were used as summary statistics to compare the data to the simulations. We confirm that evolution during acute infection is significantly different from neutral. The distribution of mutational fitness effects is best fit by a distribution with a low, but significant fraction of beneficial mutations and a high fraction of deleterious mutations. While most mutations are neutral or deleterious in this model, about 5% of mutations are beneficial. These beneficial mutations will, on average, result in a small but significant increase in fitness. When assuming no epistasis, this indicates that, at the moment of transmission, HIV is near, but not on the fitness peak for early infection.

Identification and Characterization of Domesticated Bacterial Transposases

Selfish genetic elements, such as insertion sequences and transposons are found in most genomes. Transposons are usually identifiable by their high copy number within genomes. In contrast, REP-associated tyrosine transposases (RAYTs), a recently described class of bacterial transposase, are typically present at just one copy per genome. This suggests that RAYTs no longer copy themselves and thus they no longer function as a typical transposase. Motivated by this possibility we interrogated thousands of fully sequenced bacterial genomes in order to determine patterns of RAYT diversity, their distribution across chromosomes and accessory elements, and rate of duplication. RAYTs encompass exceptional diversity and are divisible into at least five distinct groups. They possess features more similar to housekeeping genes than insertion sequences, are predominantly vertically transmitted and have persisted through evolutionary time to the point where they are now found in 24% of all species for which at least one fully sequenced genome is available. Overall, the genomic distribution of RAYTs suggests that they have been coopted by host genomes to perform a function that benefits the host cell.

Recombination of Virulence Genes in Divergent Acidovorax avenae Strains That Infect a Common Host

Bacterial etiolation and decline (BED), caused by Acidovorax avenae, is an emerging disease of creeping bentgrass on golf courses in the United States. We performed the first comprehensive analysis of A. avenae on a nationwide collection of turfgrass- and maize-pathogenic A. avenae. Surprisingly, our results reveal that the turfgrass-pathogenic A. avenae in North America are not only highly divergent but also belong to two distinct phylogroups. Both phylogroups specifically infect turfgrass but are more closely related to maize pathogens than to each other. This suggests that, although the disease is only recently reported, it has likely been infecting turfgrass for a long time. To identify a genetic basis for the host specificity, we searched for genes closely related among turfgrass strains but distantly related to their homologs from maize strains. We found a cluster of 11 such genes generated by three ancient recombination events within the type III secretion system (T3SS) pathogenicity island. Ever since the recombination, the cluster has been conserved by strong purifying selection, hinting at its selective importance. Together our analyses suggest that BED is an ancient disease that may owe its host specificity to a highly conserved cluster of 11 T3SS genes.

Discovering Complete Quasispecies in Bacterial Genomes

Mobile genetic elements can be found in almost all genomes. Possibly the most common nonautonomous mobile genetic elements in bacteria are repetitive extragenic palindromic doublets forming hairpins (REPINs) that can occur hundreds of times within a genome. The sum of all REPINs in a genome can be viewed as an evolving population because REPINs replicate and mutate. In contrast to most other biological populations, we know the exact composition of the REPIN population and the sequence of each member of the population. Here, we model the evolution of REPINs as quasispecies. We fit our quasispecies model to 10 different REPIN populations from 10 different bacterial strains and estimate effective duplication rates. Our estimated duplication rates range from ∼5 × 10-9 to 15 × 10-9 duplications per bacterial generation per REPIN. The small range and the low level of the REPIN duplication rates suggest a universal trade-off between the survival of the REPIN population and the reduction of the mutational load for the host genome. The REPIN populations we investigated also possess features typical of other natural populations. One population shows hallmarks of a population that is going extinct, another population seems to be growing in size, and we also see an example of competition between two REPIN populations.

Curiosities of REPINs and RAYTs

Repetitive extragenic palindromic (REP) sequences are a ubiquitous feature of bacterial genomes. Recent work shows that REPs are remnants of a larger mobile genetic element termed a REPIN. REPINs consists of two REP sequences in inverted orientation separated by a spacer region and are thought to be non-autonomous mobile genetic elements that exploit the transposase encoded by REP-Associated tYrosine Transposases (RAYTs). Complimentarity between the two ends of the REPIN suggests that the element forms hairpin structures in single stranded DNA or RNA. In addition to REPINs, other more complex arrangements of REPs have been identified in bacterial genomes, including the genome of the model organism Pseudomonas fluorescens SBW25. Here, we summarize existing knowledge and present new data concerning REPIN diversity. We also consider factors affecting the evolution of REPIN diversity, the ease with which REPINs might be co-opted by host genomes and the consequences of REPIN activity for the structure of bacterial genomes.

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.

Within-genome evolution of REPINs: a new family of miniature mobile DNA in bacteria

Repetitive sequences are a conserved feature of many bacterial genomes. While first reported almost thirty years ago, and frequently exploited for genotyping purposes, little is known about their origin, maintenance, or processes affecting the dynamics of within-genome evolution. Here, beginning with analysis of the diversity and abundance of short oligonucleotide sequences in the genome of Pseudomonas fluorescens SBW25, we show that over-represented short sequences define three distinct groups (GI, GII, and GIII) of repetitive extragenic palindromic (REP) sequences. Patterns of REP distribution suggest that closely linked REP sequences form a functional replicative unit: REP doublets are over-represented, randomly distributed in extragenic space, and more highly conserved than singlets. In addition, doublets are organized as inverted repeats, which together with intervening spacer sequences are predicted to form hairpin structures in ssDNA or mRNA. We refer to these newly defined entities as REPINs (REP doublets forming hairpins) and identify short reads from population sequencing that reveal putative transposition intermediates. The proximal relationship between GI, GII, and GIII REPINs and specific REP-associated tyrosine transposases (RAYTs), combined with features of the putative transposition intermediate, suggests a mechanism for within-genome dissemination. Analysis of the distribution of REPs in a range of RAYT–containing bacterial genomes, including Escherichia coli K-12 and Nostoc punctiforme, show that REPINs are a widely distributed, but hitherto unrecognized, family of miniature non-autonomous mobile DNA.

DNA sequences that copy themselves throughout genomes, and make no specific contribution to reproductive success, are by definition “selfish.” Such DNA is a feature of the genomes of all organisms and evident by virtue of its repetitive nature. In bacteria the predominant repetitive sequences are short (∼20 bp), extragenic, and palindromic. These so-called REP sequences may occur many hundreds of times per genome, but their origins and means of dissemination have been a longstanding mystery. We show that REPs are components of higher-order replicative entities termed REPINs, which are themselves thought to be derived from REP sequences that flanked an ancestral autonomous selfish element. In this ancestral state the REP sequences were likely to have been critical for the movement of the selfish element, but were devoid of any capacity to replicate independently. REPINs, on the other hand, have evolved to have a life of their own, albeit one that exploits—even enslaves—a genetic element upon which their existence depends. REPINs are the ultimate non-autonomous, super-streamlined, selfish element and are widespread among bacteria.

Global identification of yeast chromosome interactions using Genome conformation capture

The association of chromosomes with each other and other nuclear components plays a critical role in nuclear organization and Genome function. Here, using a novel and generally applicable methodology (Genome conformation capture [GCC]), we reveal the network of chromosome interactions for the yeast Saccharomyces cerevisiae. Inter- and intra-chromosomal interactions are non-random and the number of interactions per open reading frame depends upon the dispensability of the gene product. Chromosomal interfaces are organized and provide evidence of folding within chromosomes. Interestingly, the genomic connections also involve the 2 microm plasmid and the mitochondrial genome. Mitochondrial interaction partners include genes of alpha-proteobacterial origin and the ribosomal DNA. Organization of the 2 microm plasmid aligns two inverted repeats (IR1 and IR2) and displays the stability locus on a prominent loop thus making it available for plasmid clustering. Our results form the first global map of chromosomal interactions in a eukaryotic nucleus and demonstrate the highly connected nature of the yeast genome. These results have significant implications for understanding eukaryotic genome organization.

The perils of plenty: what are we going to do with all these genes?

This new century's biology promises more of everything--more genes, more organisms, more species and, in short, more data. The flood of data challenges us to find better and quicker ways to summarize and analyse. Here, we present preliminary results and proofs of concept from three of our research projects that are motivated by our search for solutions to the perils of plenty. First, we discuss how models of evolution can accommodate change to better reflect the dynamics of sequence diversity, particularly when it is becoming a lot easier to obtain sequences at different times and across intervals where the probability of new mutations contributing to this diversity is high. Second, we describe our work on the use of a single locus for species delimitation; this research targets the new DNA-barcoding approach that aims to catalogue the entirety of life. We have developed a single-locus test based on the coalescent that tests the null hypothesis of panmixis. Finally, we discuss new sequencing technologies, the types of data available and the efficacy of alignment-free methods to estimate pairwise distances for phylogenetic analyses.