Bacterial evolution by genomic island transfer occurs via DNA transformation in planta

Our understanding of the evolution of microbial pathogens has been advanced by the discovery of "islands" of DNA that differ from core genomes and contain determinants of virulence [1, 2]. The acquisition of genomic islands (GIs) by horizontal gene transfer (HGT) is thought to have played a major role in microbial evolution. There are, however, few practical demonstrations of the acquisition of genes that control virulence, and, significantly, all have been achieved outside the animal or plant host. Loss of a GI from the bean pathogen Pseudomonas syringae pv. phaseolicola (Pph) is driven by exposure to the stress imposed by the plant's resistance response [3]. Here, we show that the complete episomal island, which carries pathogenicity genes including the effector avrPphB, transfers between strains of Pph by transformation in planta and inserts at a specific att site in the genome of the recipient. Our results show that the evolution of bacterial pathogens by HGT may be achieved via transformation, the simplest mechanism of DNA exchange. This process is activated by exposure to plant defenses, when the pathogen is in greatest need of acquiring new genetic traits to alleviate the antimicrobial stress imposed by plant innate immunity [4].

Calibrating bacterial evolution

Attempts to calibrate bacterial evolution have relied on the assumption that rates of molecular sequence divergence in bacteria are similar to those of higher eukaryotes, or to those of the few bacterial taxa for which ancestors can be reliably dated from ecological or geological evidence. Despite similarities in the substitution rates estimated for some lineages, comparisons of the relative rates of evolution at different classes of nucleotide sites indicate no basis for their universal application to all bacteria. However, there is evidence that bacteria have a constant genome-wide mutation rate on an evolutionary time scale but that this rate differs dramatically from the rate estimated by experimental methods.

Bayesian statistical analysis of bacterial diversity

Bacteria play an important role in many ecological systems. The molecular characterization of bacteria using either cultivation-dependent or cultivation-independent methods reveals the large scale of bacterial diversity in natural communities, and the vastness of subpopulations within a species or genus. Understanding how bacterial diversity varies across different environments and also within populations should provide insights into many important questions of bacterial evolution and population dynamics. This thesis presents novel statistical methods for analyzing bacterial diversity using widely employed molecular fingerprinting techniques. The first objective of this thesis was to develop Bayesian clustering models to identify bacterial population structures. Bacterial isolates were identified using multilous sequence typing (MLST), and Bayesian clustering models were used to explore the evolutionary relationships among isolates. Our method involves the inference of genetic population structures via an unsupervised clustering framework where the dependence between loci is represented using graphical models. The population dynamics that generate such a population stratification were investigated using a stochastic model, in which homologous recombination between subpopulations can be quantified within a gene flow network. The second part of the thesis focuses on cluster analysis of community compositional data produced by two different cultivation-independent analyses: terminal restriction fragment length polymorphism (T-RFLP) analysis, and fatty acid methyl ester (FAME) analysis. The cluster analysis aims to group bacterial communities that are similar in composition, which is an important step for understanding the overall influences of environmental and ecological perturbations on bacterial diversity. A common feature of T-RFLP and FAME data is zero-inflation, which indicates that the observation of a zero value is much more frequent than would be expected, for example, from a Poisson distribution in the discrete case, or a Gaussian distribution in the continuous case. We provided two strategies for modeling zero-inflation in the clustering framework, which were validated by both synthetic and empirical complex data sets. We show in the thesis that our model that takes into account dependencies between loci in MLST data can produce better clustering results than those methods which assume independent loci. Furthermore, computer algorithms that are efficient in analyzing large scale data were adopted for meeting the increasing computational need. Our method that detects homologous recombination in subpopulations may provide a theoretical criterion for defining bacterial species. The clustering of bacterial community data include T-RFLP and FAME provides an initial effort for discovering the evolutionary dynamics that structure and maintain bacterial diversity in the natural environment.

Within-host evolution of Burkholderia pseudomallei over a twelve-year chronic carriage infection

UNLABELLED: Burkholderia pseudomallei causes the potentially fatal disease melioidosis. It is generally accepted that B. pseudomallei is a noncommensal bacterium and that any culture-positive clinical specimen denotes disease requiring treatment. Over a 23-year study of melioidosis cases in Darwin, Australia, just one patient from 707 survivors has developed persistent asymptomatic B. pseudomallei carriage. To better understand the mechanisms behind this unique scenario, we performed whole-genome analysis of two strains isolated 139 months apart. During this period, B. pseudomallei underwent several adaptive changes. Of 23 point mutations, 78% were nonsynonymous and 43% were predicted to be deleterious to gene function, demonstrating a strong propensity for positive selection. Notably, a nonsense mutation inactivated the universal stress response sigma factor RpoS, with pleiotropic implications. The genome underwent substantial reduction, with four deletions in chromosome 2 resulting in the loss of 221 genes. The deleted loci included genes involved in secondary metabolism, environmental survival, and pathogenesis. Of 14 indels, 11 occurred in coding regions and 9 resulted in frameshift mutations that dramatically affected predicted gene products. Disproportionately, four indels affected lipopolysaccharide biosynthesis and modification. Finally, we identified a frameshift mutation in both P314 isolates within wcbR, an important component of the capsular polysaccharide I locus, suggesting virulence attenuation early in infection. Our study illustrates a unique clinical case that contrasts a high-consequence infectious agent with a long-term commensal infection and provides further insights into bacterial evolution within the human host.

IMPORTANCE: Some bacterial pathogens establish long-term infections that are difficult or impossible to eradicate with current treatments. Rapid advances in genome sequencing technologies provide a powerful tool for understanding bacterial persistence within the human host. Burkholderia pseudomallei is considered a highly pathogenic bacterium because infection is commonly fatal. Here, we document within-host evolution of B. pseudomallei in a unique case of human infection with ongoing chronic carriage. Genomic comparison of isolates obtained 139 months (11.5 years) apart showed a strong signal of adaptation within the human host, including inactivation of virulence and immunogenic factors, and deletion of pathways involved in environmental survival. Two global regulatory genes were mutated in the 139-month isolate, indicating extensive regulatory changes favoring bacterial persistence. Our study provides insights into B. pseudomallei pathogenesis and, more broadly, identifies parallel evolutionary mechanisms that underlie chronic persistence of all bacterial pathogens.

Use of bacterial whole-genome sequencing to investigate local persistence and spread in bovine tuberculosis

Mycobacterium bovis is the causal agent of bovine tuberculosis, one of the most important diseases currently facing the UK cattle industry. Here, we use high-density whole genome sequencing (WGS) in a defined sub-population of M. bovis in 145 cattle across 66 herd breakdowns to gain insights into local spread and persistence. We show that despite low divergence among isolates, WGS can in principle expose contributions of under-sampled host populations to M. bovis transmission. However, we demonstrate that in our data such a signal is due to molecular type switching, which had been previously undocumented for M. bovis. Isolates from farms with a known history of direct cattle movement between them did not show a statistical signal of higher genetic similarity. Despite an overall signal of genetic isolation by distance, genetic distances also showed no apparent relationship with spatial distance among affected farms over distances <5 km. Using simulations, we find that even over the brief evolutionary timescale covered by our data, Bayesian phylogeographic approaches are feasible. Applying such approaches showed that M. bovis dispersal in this system is heterogeneous but slow overall, averaging 2 km/year. These results confirm that widespread application of WGS to M. bovis will bring novel and important insights into the dynamics of M. bovis spread and persistence, but that the current questions most pertinent to control will be best addressed using approaches that more directly integrate WGS with additional epidemiological data.

Convergent evolution of PICIs-mediated phage interference

Staphylococcal pathogenicity islands (SaPIs), the prototype members of the family of phage inducible chromosomal islands (PICIs), are extremely mobile phage satellites, which are transferred between bacterial hosts after their induction by a helper phage. The intimate relationship between SaPIs and their helper phages is one of the most studied examples of virus satellite interactions in prokaryotic cells. SaPIs encode and disseminate virulence and fitness factors, representing a driving force for bacterial adaptation and pathogenesis. Many SaPIs encode a conserved morphogenetic operon, including a core set of genes whose function allows them to parasitize and exploit the phage life cycle. One of the central mechanisms of this molecular piracy is the specific packaging of the SaPI genomes into reduced sized capsid structures derived from phage proteins. Pac phages were classically thought to be the only phages involved in the mobilisation of phage-mediated virulence genes, including the transfer of SaPIs within related and non-related bacteria. This study presents the involvement of S. aureus cos phages in the intra- and intergeneric transfer of cos SaPIs for the first time. A novel example of molecular parasitism is shown, by which this newly characterised group of cos SaPIs uses two distinct and complementary mechanisms to take over the helper phage packaging machinery for their own reproduction. SaPIbov5, the prototype of the cos SaPIs, does not encode the characteristic morphogenetic operon found in pac SaPIs. However, cos SaPIs features both pac and cos phage cleavage sequences in their genome, ensuring SaPI packaging in small- and full-sized phage particles, depending on the helper phage. Moreover, cos-site packaging in S. aureus was shown to require the activity of a phage HNH nuclease. The HNH protein functions together with the large terminase subunit, triggering cleavage and melting of the cos-site sequence. In addition, a novel piracy strategy, severely interfering with the helper phage reproduction, was identified in cos SaPIs and characterised. This mechanism of piracy depends on the cos SaPI-encoded ccm gene, which encodes a capsid protein involved in the formation of small phage particles, modifying the assembling process via a scaffolding mechanism. This strategy resembles the ones described for pac SaPIs and represents a remarkable example of convergent evolution. A further convergent mechanism of capsid size-reduction was identified and characterised for the Enterococcus faecalis EfCIV583 pathogenicity island, another member of the PICI family. In this case, the self-encoded CpmE conducts this molecular piracy through a putative scaffolding function. Similar to cos SaPIs, EfCIV583 carries the helper phage cleavage sequence in its genome enabling its mobilisation by the phage terminase complex. The results presented in this thesis show how two examples of non-related members of the PICI family follow the same evolutionary convergent strategy to interfere with their helper phage. These findings could indicate that the described strategies might be widespread among PICIs and implicate a significant impact of PICIs mediated-virulence gene transfer in bacterial evolution and the emergence of pathogenic bacteria.

Genotypic and phenotypic diversity of the noncapsulated Haemophilus influenzae:Adaptation and pathogenesis in the human airways

The human respiratory tract contains a highly adapted microbiota including commensal and opportunistic pathogens. Noncapsulated or nontypable Haemophilus influenzae (NTHi) is a human-restricted member of the normal airway microbiota in healthy carriers and an opportunistic pathogen in immunocompromised individuals. The duality of NTHi as a colonizer and as a symptomatic infectious agent is closely related to its adaptation to the host, which in turn greatly relies on the genetic plasticity of the bacterium and is facilitated by its condition as a natural competent. The variable genotype of NTHi accounts for its heterogeneous gene expression and variable phenotype, leading to differential host-pathogen interplay among isolates. Here we review our current knowledge of NTHi diversity in terms of genotype, gene expression, antigenic variation, and the phenotypes associated with colonization and pathogenesis. The potential benefits of NTHi diversity studies discussed herein include the unraveling of pathogenicity clues, the generation of tools to predict virulence from genomic data, and the exploitation of a unique natural system for the continuous monitoring of long-term bacterial evolution in human airways exposed to noxious agents. Finally, we highlight the challenge of monitoring both the pathogen and the host in longitudinal studies, and of applying comparative genomics to clarify the meaning of the vast NTHi genetic diversity and its translation to virulence phenotypes.

Mobile elements drive recombination hotspots in the core genome of Staphylococcus aureus

Horizontal gene transfer is an important driver of bacterial evolution, but genetic exchange in the core genome of clonal species, including the major pathogen Staphylococcus aureus, is incompletely understood. Here we reveal widespread homologous recombination in S. aureus at the species level, in contrast to its near-complete absence between closely related strains. We discover a patchwork of hotspots and coldspots at fine scales falling against a backdrop of broad-scale trends in rate variation. Over megabases, homoplasy rates fluctuate 1.9-fold, peaking towards the origin-of-replication. Over kilobases, we find core recombination hotspots of up to 2.5-fold enrichment situated near fault lines in the genome associated with mobile elements. The strongest hotspots include regions flanking conjugative transposon ICE6013, the staphylococcal cassette chromosome (SCC) and genomic island νSaα. Mobile element-driven core genome transfer represents an opportunity for adaptation and challenges our understanding of the recombination landscape in predominantly clonal pathogens, with important implications for genotype–phenotype mapping.

Genome sequence of Fusobacterium nucleatum subspecies polymorphum - a genetically tractable fusobacterium.

Fusobacterium nucleatum is a prominent member of the oral microbiota and is a common cause of human infection. F. nucleatum includes five subspecies: polymorphum, nucleatum, vincentii, fusiforme, and animalis. F. nucleatum subsp. polymorphum ATCC 10953 has been well characterized phenotypically and, in contrast to previously sequenced strains, is amenable to gene transfer. We sequenced and annotated the 2,429,698 bp genome of F. nucleatum subsp. polymorphum ATCC 10953. Plasmid pFN3 from the strain was also sequenced and analyzed. When compared to the other two available fusobacterial genomes (F. nucleatum subsp. nucleatum, and F. nucleatum subsp. vincentii) 627 open reading frames unique to F. nucleatum subsp. polymorphum ATCC 10953 were identified. A large percentage of these mapped within one of 28 regions or islands containing five or more genes. Seventeen percent of the clustered proteins that demonstrated similarity were most similar to proteins from the clostridia, with others being most similar to proteins from other gram-positive organisms such as Bacillus and Streptococcus. A ten kilobase region homologous to the Salmonella typhimurium propanediol utilization locus was identified, as was a prophage and integrated conjugal plasmid. The genome contains five composite ribozyme/transposons, similar to the CdISt IStrons described in Clostridium difficile. IStrons are not present in the other fusobacterial genomes. These findings indicate that F. nucleatum subsp. polymorphum is proficient at horizontal gene transfer and that exchange with the Firmicutes, particularly the Clostridia, is common.

Experimental surgery to create subgenomes of Bacillus subtilis 168

The 4,188-kb circular genome of Bacillus subtilis 168 was artificially dissected into two stable circular chromosomes in vivo, one being the 3,878-kb main genome and the other the 310-kb subgenome that was recovered as covalently closed circular DNA in CsCl-ethidium bromide ultracentrifugation. The minimal requirements to physically separate the 310-kb DNA segment out of the genome were two interrepeat homologous sequences and an origin of DNA replication between them. The subgenome originated from the 1,255–1,551-kb region of the B. subtilis genome was essential for the cell to survive because the subgenome was not lost from the cell. The finding that the B. subtilis genome has a potential to be divided and the resulting two replicons stably maintained may shed light on origins and formation mechanisms of giant plasmids or second chromosomes present in many bacteria. Similar excision or its reversal process, i.e., integration of large sized covalently closed circular DNA pieces into the main genome, implies significant roles of subgenomes in the exchange of genetic information and size variation of bacterial genomes in bacterial evolution.

Intracellular compartmentation in planctomycetes

The phylum Planctomycetes of the domain Bacteria consists of budding, peptidoglycan-less organisms important for understanding the origins of complex cell organization. Their significance for cell biology lies in their possession of intracellular membrane compartmentation. All planctomycetes share a unique cell plan, in which the cell cytoplasm is divided into compartments by one or more membranes, including a major cell compartment containing the nucleoid. Of special significance is Gemmata obscuriglobus, in which the nucleoid is enveloped in two membranes to form a nuclear body that is analogous to the structure of a eukaryotic nucleus. Planctomycete compartmentation may have functional physiological roles, as in the case of anaerobic ammonium-oxidizing anammox planctomycetes, in which the anammoxosome harbors specialized enzymes and is wrapped in an envelope possessing unique ladderane lipids. Organisms in phyla other than the phylum Planctomycetes may possess compartmentation similar to that of some planctomycetes, as in the case of members of the phylum Poribacteria from marine sponges.

An updated evolutionary classification of CRISPR–Cas systems

The evolution of CRISPR–cas loci, which encode adaptive immune systems in archaea and bacteria, involves rapid changes, in particular numerous rearrangements of the locus architecture and horizontal transfer of complete loci or individual modules. These dynamics complicate straightforward phylogenetic classification, but here we present an approach combining the analysis of signature protein families and features of the architecture of cas loci that unambiguously partitions most CRISPR–cas loci into distinct classes, types and subtypes. The new classification retains the overall structure of the previous version but is expanded to now encompass two classes, five types and 16 subtypes. The relative stability of the classification suggests that the most prevalent variants of CRISPR–Cas systems are already known. However, the existence of rare, currently unclassifiable variants implies that additional types and subtypes remain to be characterized.

Evolution of bacterial transcription factors: how proteins take on new tasks, but do not always stop doing the old ones

Many bacterial transcription factors do not behave as per the textbook operon model. We draw on whole genome work, as well as reported diversity across different bacteria, to argue that transcription factors may have evolved from nucleoid-associated proteins. This view would explain a large amount of recent data gleaned from high-throughput sequencing and bioinformatic analyses.