Polymorphic toxin systems: Comprehensive characterization of trafficking modes, processing, mechanisms of action, immunity and ecology using comparative genomics

Background: Proteinaceous toxins are observed across all levels of inter-organismal and intra-genomic conflicts. These include recently discovered prokaryotic polymorphic toxin systems implicated in intra-specific conflicts. They are characterized by a remarkable diversity of C-terminal toxin domains generated by recombination with standalone toxin-coding cassettes. Prior analysis revealed a striking diversity of nuclease and deaminase domains among the toxin modules. We systematically investigated polymorphic toxin systems using comparative genomics, sequence and structure analysis. Results: Polymorphic toxin systems are distributed across all major bacterial lineages and are delivered by at least eight distinct secretory systems. In addition to type-II, these include type-V, VI, VII (ESX), and the poorly characterized "Photorhabdus virulence cassettes (PVC)", PrsW-dependent and MuF phage-capsid-like systems. We present evidence that trafficking of these toxins is often accompanied by autoproteolytic processing catalyzed by HINT, ZU5, PrsW, caspase-like, papain-like, and a novel metallopeptidase associated with the PVC system. We identified over 150 distinct toxin domains in these systems. These span an extraordinary catalytic spectrum to include 23 distinct clades of peptidases, numerous previously unrecognized versions of nucleases and deaminases, ADP-ribosyltransferases, ADP ribosyl cyclases, RelA/SpoT-like nucleotidyltransferases, glycosyltranferases and other enzymes predicted to modify lipids and carbohydrates, and a pore-forming toxin domain. Several of these toxin domains are shared with host-directed effectors of pathogenic bacteria. Over 90 families of immunity proteins might neutralize anywhere between a single to at least 27 distinct types of toxin domains. In some organisms multiple tandem immunity genes or immunity protein domains are organized into polyimmunity loci or polyimmunity proteins. Gene-neighborhood-analysis of polymorphic toxin systems predicts the presence of novel trafficking-related components, and also the organizational logic that allows toxin diversification through recombination. Domain architecture and protein-length analysis revealed that these toxins might be deployed as secreted factors, through directed injection, or via inter-cellular contact facilitated by filamentous structures formed by RHS/YD, filamentous hemagglutinin and other repeats. Phyletic pattern and life-style analysis indicate that polymorphic toxins and polyimmunity loci participate in cooperative behavior and facultative 'cheating' in several ecosystems such as the human oral cavity and soil. Multiple domains from these systems have also been repeatedly transferred to eukaryotes and their viruses, such as the nucleo-cytoplasmic large DNA viruses. Conclusions: Along with a comprehensive inventory of toxins and immunity proteins, we present several testable predictions regarding active sites and catalytic mechanisms of toxins, their processing and trafficking and their role in intra-specific and inter-specific interactions between bacteria. These systems provide insights regarding the emergence of key systems at different points in eukaryotic evolution, such as ADP ribosylation, interaction of myosin VI with cargo proteins, mediation of apoptosis, hyphal heteroincompatibility, hedgehog signaling, arthropod toxins, cell-cell interaction molecules like teneurins and different signaling messengers.

4-chlorophenol biodegradation by Arthrobacter chlorophenolicus A6

A microorganism was isolated which could grow on unusually high concentrations of the toxic pollutant 4-chlorophenol. Taxonomic studies showed that the microorganism constituted a novel species within the genus Arthrobacter and it was named Arthrobacter chlorophenolicus A6. A. chlorophenolicus A6 was chromosomally tagged with either the gfp gene, encoding the green fluorescent protein (GFP), or the luc gene, encoding firefly luciferase. When the tagged cells were inoculated into 4-chlorophenol contaminated soil they could completely remove 175 µg/g 4-chlorophenol within 10 days, whereas no loss of 4-chlorophenol was observed in the uninoculated control microcosms. During these experiments the gfp and luc marker genes allowed monitoring of cell number and metabolic status. When A. chlorophenolicus A6 was grown on mixtures of phenolic compounds, the strain exhibited a preference for 4-nitrophenol over 4-chlorophenol, which in turn was preferred over phenol. Analysis of growth and degradation data indicated that the same enzyme system was used for removal of 4-chlorophenol and 4-nitrophenol. However, degradation of unbstituted phenol appeared to be mediated by another or an additional enzyme system. The luc-tagged A. chlorophenolicus A6 gave valuable information about growth, substrate depletion and toxicity of the phenolic compounds in substrate mixtures. The 4-chlorophenol degradation pathway in A. chlorophenolicus A6 was elucidated. The metabolic intermediate subject to ring cleavage was found to be hydroxyquinol and two different pathway branches led from 4-chlorophenol to hydroxyquinol. A gene cluster involved in 4-chlorophenol degradation was cloned from A. chlorophenolicus A6. The cluster contained two functional hydroxyquinol 1,2-dioxygenase genes and a number of other open reading frames presumed to encode enzymes involved in 4-chlorophenol catabolism. Analysis of the DNA sequence suggested that the gene cluster had partly been assembled by horizontal gene transfer. In summary, 4-chlorophenol degradation by A. chlorophenolicus A6 was studied from a number of angles. This organism has several interesting and useful traits such as the ability to degrade high concentrations of 4-chlorophenol and other phenols alone and in mixtures, an unusual and effective 4-chlorophenol degradation pathway and demonstrated ability to remove 4-chlorophenol from contaminated soil.

Meccanismi evolutivi dei parapoxvirus: caratterizzazione genomica di pseudocowpoxvirus e messa a punto di sistemi per lo studio delle ricombinazioni

The Poxviruses are a family of double stranded DNA (dsDNA) viruses that cause disease in many species, both vertebrate and invertebrate. Their genomes range in size from 135 to 365 kbp and show conservation in both organization and content. In particular, the central genomic regions of the chordopoxvirus subfamily (those capable of infecting vertebrates) contain 88 genes which are present in all the virus species characterised to date and which mostly occur in the same order and orientation. In contrast, however, the terminal regions of the genomes frequently contain genes that are species or genera-specific and that are not essential for the growth of the virus in vitro but instead often encode factors with important roles in vivo including modulation of the host immune response to infection and determination of the host range of the virus. The Parapoxviruses (PPV), of which Orf virus is the prototypic species, represent a genus within the chordopoxvirus subfamily of Poxviridae and are characterised by their ability to infect ruminants and humans. The genus currently contains four recognised species of virus, bovine papular stomatitis virus (BPSV) and pseudocowpox virus (PCPV) both of which infect cattle, orf virus (OV) that infects sheep and goats, and parapoxvirus of red deer in New Zealand (PVNZ). The ORFV genome has been fully sequenced, as has that of BPSV, and is ~138 kb in length encoding ~132 genes. The vast majority of these genes allow the virus to replicate in the cytoplasm of the infected host cell and therefore encode proteins involved in replication, transcription and metabolism of nucleic acids. These genes are well conserved between all known genera of poxviruses. There is however another class of genes, located at either end of the linear dsDNA genome, that encode proteins which are non-essential for replication and generally dictate host range and virulence of the virus. The non-essential genes are often the most variable within and between species of virus and therefore are potentially useful for diagnostic purposes. Given their role in subverting the host-immune response to infection they are also targets for novel therapeutics. The function of only a relatively small number of these proteins has been elucidated and there are several genes whose function still remains obscure principally because there is little similarity between them and proteins of known function in current sequence databases. It is thought that by selectively removing some of the virulence genes, or at least neutralising the proteins in some way, current vaccines could be improved. The evolution of poxviruses has been proposed to be an adaptive process involving frequent events of gene gain and loss, such that the virus co-evolves with its specific host. Gene capture or horizontal gene transfer from the host to the virus is considered an important source of new viral genes including those likely to be involved in host range and those enabling the virus to interfere with the host immune response to infection. Given the low rate of nucleotide substitution, recombination can be seen as an essential evolutionary driving force although it is likely underestimated. Recombination in poxviruses is intimately linked to DNA replication with both viral and cellular proteins participate in this recombination-dependent replication. It has been shown, in other poxvirus genera, that recombination between isolates and perhaps even between species does occur, thereby providing another mechanism for the acquisition of new genes and for the rapid evolution of viruses. Such events may result in viruses that have a selective advantage over others, for example in re-infections (a characteristic of the PPV), or in viruses that are able to jump the species barrier and infect new hosts. Sequence data related to viral strains isolated from goats suggest that possible recombination events may have occurred between OV and PCPV (Ueda et al. 2003). The recombination events are frequent during poxvirus replication and comparative genomic analysis of several poxvirus species has revealed that recombinations occur frequently on the right terminal region. Intraspecific recombination can occur between strains of the same PPV species, but also interspecific recombination can happen depending on enough sequence similarity to enable recombination between distinct PPV species. The most important pre-requisite for a successful recombination is the coinfection of the individual host by different virus strains or species. Consequently, the following factors affecting the distribution of different viruses to shared target cells need to be considered: dose of inoculated virus, time interval between inoculation of the first and the second virus, distance between the marker mutations, genetic homology. At present there are no available data on the replication dynamics of PPV in permissive and non permissive hosts and reguarding co-infetions there are no information on the interference mechanisms occurring during the simultaneous replication of viruses of different species. This work has been carried out to set up permissive substrates allowing the replication of different PPV species, in particular keratinocytes monolayers and organotypic skin cultures. Furthermore a method to isolate and expand ovine skin stem cells was has been set up to indeep further aspects of viral cellular tropism during natural infection. The study produced important data to elucidate the replication dynamics of OV and PCPV virus in vitro as well as the mechanisms of interference that can arise during co-infection with different viral species. Moreover, the analysis carried on the genomic right terminal region of PCPV 1303/05 contributed to a better knowledge of the viral genes involved in host interaction and pathogenesis as well as to locate recombination breakpoints and genetic homologies between PPV species. Taken together these data filled several crucial gaps for the study of interspecific recombinations of PPVs which are thought to be important for a better understanding of the viral evolution and to improve the biosafety of antiviral therapy and PPV-based vectors.

Influence of pollination and seed dispersal on the genetic population structure of two Commiphora species in Madagascar and South Africa

Pollination and seed dispersal are important ecological processes for the regeneration of plant populations and both vectors for gene exchange between plant populations. For my thesis, I studied the pollination ecology of the South African tree Commiphora harveyi (Burseraceae) and compared it with C. guillauminii from Madagascar. Both species have low visitation rates and a low number of pollinating insect species, resulting in a low fruit set. While their pollination ecology is very similar, they differ in their seed dispersal with a low seed dispersal rate in the Malagasy and a high seed dispersal rate in the South African species. This should be reflected in a stronger genetic differentiation among populations in the Malagasy than in the South African species. My results, based on AFLP markers, contradict these expectations, the overall differentiation was lower in the Malagasy (FST = 0.05) than in the South African species (FST = 0.16). However, at a smaller spatial scale (below 3 km), the Malagasy species was genetically more strongly differentiated than the South African species, which was reflected by the high inter-population variance within the sample site (C. guillauminii: 72.2 - 85.5 %; C. harveyi: 8.4 - 14.5 %). This strong differentiation could arise from limited gene flow, which was confirmed by spatial autocorrelation analyses. The shape of the autocorrelogram suggested that gene exchange between individuals occurred only up to 3 km in the Malagasy species, whereas up to 30 km in the South African species. These results on the genetic structure correspond to the expectations based on seed dispersal data. Thus, seed dispersal seems to be a key factor for the genetic structure in plant populations on a local scale.

Functional and Genetic characterization of new genomic islands from an E. coli strain associated with neonatal meningitis.

Enterobacteriaceae genomes evolve through mutations, rearrangements and horizontal gene transfer (HGT). The latter evolutionary pathway works through the acquisition DNA (GEI) modules of foreign origin that enhances fitness of the host to a given environment. The genome of E. coli IHE3034, a strain isolated from a case of neonatal meningitis, has recently been sequenced and its subsequent sequence analysis has predicted 18 possible GEIs, of which: 8 have not been previously described, 5 fully meet the pathogenic island definition and at least 10 that seem to be of prophagic origin. In order to study the GEI distribution of our reference strain, we screened for the presence 18 GEIs a panel of 132 strains, representative of E. coli diversity. Also, using an inverse nested PCR approach we identified 9 GEI that can form an extrachromosomal circular intermediate (CI) and their respective attachment sites (att). Further, we set up a qPCR approach that allowed us to determine the excision rates of 5 genomic islands in different growth conditions. Four islands, specific for strains appertaining to the sequence type complex 95 (STC95), have been deleted in order to assess their function in a Dictyostelium discoideum grazing assays. Overall, the distribution data presented here indicate that 16 IHE3034 GEIs are more associated to the STC95 strains. Also the functional and genetic characterization has uncovered that GEI 13, 17 and 19 are involved in the resistance to phagocitation by Dictyostelium d thus suggesting a possible role in the adaptation of the pathogen during certain stages of infection.

Multiple colonization with S. pneumoniae before and after introduction of the seven-valent conjugated pneumococcal polysaccharide vaccine

Simultaneous carriage of more than one strain of Streptococcus pneumoniae promotes horizontal gene transfer events and may lead to capsule switch and acquisition of antibiotic resistance. We studied the epidemiology of cocolonization with S. pneumoniae before and after introduction of the seven-valent conjugated pneumococcal vaccine (PCV7).

Detection of Streptococcus pneumoniae strain cocolonization in the nasopharynx

Colonization with more than one distinct strain of the same species, also termed cocolonization, is a prerequisite for horizontal gene transfer between pneumococcal strains that may lead to change of the capsular serotype. Capsule switch has become an important issue since the introduction of conjugated pneumococcal polysaccharide vaccines. There is, however, a lack of techniques to detect multiple colonization by S. pneumoniae strains directly in nasopharyngeal samples. Two hundred eighty-seven nasopharyngeal swabs collected during the prevaccine era within a nationwide surveillance program were analyzed by a novel technique for the detection of cocolonization, based on PCR amplification of a noncoding region adjacent to the pneumolysin gene (plyNCR) and restriction fragment length polymorphism (RFLP) analysis. The numbers of strains and their relative abundance in cocolonized samples were determined by terminal RFLP. The pneumococcal carriage rate found by PCR was 51.6%, compared to 40.0% found by culture. Cocolonization was present in 9.5% (10/105) of samples, most (9/10) of which contained two strains in a ratio of between 1:1 and 17:1. Five of the 10 cocolonized samples showed combinations of vaccine types only (n = 2) or combinations of nonvaccine types only (n = 3). Carriers of multiple pneumococcal strains had received recent antibiotic treatment more often than those colonized with a single strain (33% versus 9%, P = 0.025). This new technique allows for the rapid and economical study of pneumococcal cocolonization in nasopharyngeal swabs. It will be valuable for the surveillance of S. pneumoniae epidemiology under vaccine selection pressure.

Evolution of sensory complexity recorded in a myxobacterial genome.

Myxobacteria are single-celled, but social, eubacterial predators. Upon starvation they build multicellular fruiting bodies using a developmental program that progressively changes the pattern of cell movement and the repertoire of genes expressed. Development terminates with spore differentiation and is coordinated by both diffusible and cell-bound signals. The growth and development of Myxococcus xanthus is regulated by the integration of multiple signals from outside the cells with physiological signals from within. A collection of M. xanthus cells behaves, in many respects, like a multicellular organism. For these reasons M. xanthus offers unparalleled access to a regulatory network that controls development and that organizes cell movement on surfaces. The genome of M. xanthus is large (9.14 Mb), considerably larger than the other sequenced delta-proteobacteria. We suggest that gene duplication and divergence were major contributors to genomic expansion from its progenitor. More than 1,500 duplications specific to the myxobacterial lineage were identified, representing >15% of the total genes. Genes were not duplicated at random; rather, genes for cell-cell signaling, small molecule sensing, and integrative transcription control were amplified selectively. Families of genes encoding the production of secondary metabolites are overrepresented in the genome but may have been received by horizontal gene transfer and are likely to be important for predation.

Decrease in pneumococcal co-colonization following vaccination with the seven-valent pneumococcal conjugate vaccine.

Understanding the epidemiology of pneumococcal co-colonization is important for monitoring vaccine effectiveness and the occurrence of horizontal gene transfer between pneumococcal strains. In this study we aimed to evaluate the impact of the seven-valent pneumococcal conjugate vaccine (PCV7) on pneumococcal co-colonization among Portuguese children. Nasopharyngeal samples from children up to 6 years old yielding a pneumococcal culture were clustered into three groups: pre-vaccine era (n = 173), unvaccinated children of the vaccine era (n = 169), and fully vaccinated children (4 doses; n = 150). Co-colonization, serotype identification, and relative serotype abundance were detected by analysis of DNA of the total bacterial growth of the primary culture plate using the plyNCR-RFLP method and a molecular serotyping microarray-based strategy. The plyNCR-RFLP method detected an overall co-colonization rate of 20.1%. Microarray analysis confirmed the plyNCR-RFLP results. Vaccination status was the only factor found to be significantly associated with co-colonization: co-colonization rates were significantly lower (p = 0.004; Fisher's exact test) among fully vaccinated children (8.0%) than among children from the pre-PCV7 era (17.3%) or unvaccinated children of the PCV7 era (18.3%). In the PCV7 era there were significantly less non-vaccine type (NVT) co-colonization events than would be expected based on the NVT distribution observed in the pre-PCV7 era (p = 0.024). In conclusion, vaccination with PCV7 resulted in a lower co-colonization rate due to an asymmetric distribution between NVTs found in single and co-colonized samples. We propose that some NVTs prevalent in the PCV7 era are more competitive than others, hampering their co-existence in the same niche. This result may have important implications since a decrease in co-colonization events is expected to translate in decreased opportunities for horizontal gene transfer, hindering pneumococcal evolution events such as acquisition of antibiotic resistance determinants or capsular switch. This might represent a novel potential benefit of conjugate vaccines.

Copying, emergence, and digital reproduction

This talk combines concepts of the copying of manuscripts in the age before mechanical reproduction (as described in W. Benjamin’s famous Artwork-essay), of the emergence of non-vertical evolution (as to be found in textual ‘contamination’ and in horizontal gene transfer alike), and of the electronic reproduction of such phenomena (as presented in scholarly digital editions and phylogenetic trees). The guiding idea is that ‘copying’, ‘emergence’, and ‘(digital) reproduction’ enable variation depending on particular physical or biological forms, on contextual or environmental conditions, as well as on user habits or receptive fields. – Is it possible to develop a theory of copying and reproduction on this base? The material of the talk will be drawn from the Parzival-Project, a critical electronic edition of an Arthurian romance, composed by the German poet Wolfram von Eschenbach shortly after 1200 and transmitted in over eighty witnesses.

The sequence comparison index: A novel method for comparing proteins and proteomes

Historically morphological features were used as the primary means to classify organisms. However, the age of molecular genetics has allowed us to approach this field from the perspective of the organism's genetic code. Early work used highly conserved sequences, such as ribosomal RNA. The increasing number of complete genomes in the public data repositories provides the opportunity to look not only at a single gene, but at organisms' entire parts list. ^ Here the Sequence Comparison Index (SCI) and the Organism Comparison Index (OCI), algorithms and methods to compare proteins and proteomes, are presented. The complete proteomes of 104 sequenced organisms were compared. Over 280 million full Smith-Waterman alignments were performed on sequence pairs which had a reasonable expectation of being related. From these alignments a whole proteome phylogenetic tree was constructed. This method was also used to compare the small subunit (SSU) rRNA from each organism and a tree constructed from these results. The SSU rRNA tree by the SCI/OCI method looks very much like accepted SSU rRNA trees from sources such as the Ribosomal Database Project, thus validating the method. The SCI/OCI proteome tree showed a number of small but significant differences when compared to the SSU rRNA tree and proteome trees constructed by other methods. Horizontal gene transfer does not appear to affect the SCI/OCI trees until the transferred genes make up a large portion of the proteome. ^ As part of this work, the Database of Related Local Alignments (DaRLA) was created and contains over 81 million rows of sequence alignment information. DaRLA, while primarily used to build the whole proteome trees, can also be applied shared gene content analysis, gene order analysis, and creating individual protein trees. ^ Finally, the standard BLAST method for analyzing shared gene content was compared to the SCI method using 4 spirochetes. The SCI system performed flawlessly, finding all proteins from one organism against itself and finding all the ribosomal proteins between organisms. The BLAST system missed some proteins from its respective organism and failed to detect small ribosomal proteins between organisms. ^

Whole genome analysis of Fusobacterium nucleatum, Rickettsia typhi and Francisella tularensis

The genomes of Fusobacterium nucleatum subspecies polymorphum strain ATCC 10953, Rickettsia typhi strain Wilmington, and Francisella tularensis subspecies holarctica strain OSU18 were sequenced, annotated, and analyzed. Each genome was then compared to the sequenced genomes of closely related bacteria. The genome of F. nucleatum ATCC 10953 was compared to two additional F. nucleatum subspecies, subspecies nucleatum and subspecies vincentii. This analysis revealed substantial evidence of horizontal gene transfer along with considerable genetic diversity within the species of F. nucleatum. R. typhi was compared to R. prowazekii and R. conorii. This analysis uncovered a hotspot for chromosomal rearrangements in the Spotted Fever Group but not the Typhus Group Rickettsia and revealed the close genetic relationship between the Typhus Group rickettsial species. F. tularensis OSU18 was compared to two additional F. tularensis strains. These comparisons uncovered significant chromosomal rearrangements between F. tularensis subspecies due to recombination between insertion sequence elements. ^

Modeling, simulation and application of bacterial transduction in genetic algorithms

At present, all methods in Evolutionary Computation are bioinspired by the fundamental principles of neo-Darwinism, as well as by a vertical gene transfer. Virus transduction is one of the key mechanisms of horizontal gene propagation in microorganisms (e.g. bacteria). In the present paper, we model and simulate a transduction operator, exploring the possible role and usefulness of transduction in a genetic algorithm. The genetic algorithm including transduction has been named PETRI (abbreviation of Promoting Evolution Through Reiterated Infection). Our results showed how PETRI approaches higher fitness values as transduction probability comes close to 100%. The conclusion is that transduction improves the performance of a genetic algorithm, assuming a population divided among several sub-populations or ?bacterial colonies?.

Interpreting the universal phylogenetic tree

The universal phylogenetic tree not only spans all extant life, but its root and earliest branchings represent stages in the evolutionary process before modern cell types had come into being. The evolution of the cell is an interplay between vertically derived and horizontally acquired variation. Primitive cellular entities were necessarily simpler and more modular in design than are modern cells. Consequently, horizontal gene transfer early on was pervasive, dominating the evolutionary dynamic. The root of the universal phylogenetic tree represents the first stage in cellular evolution when the evolving cell became sufficiently integrated and stable to the erosive effects of horizontal gene transfer that true organismal lineages could exist.

Postsymbiotic plasmid acquisition and evolution of the repA1-replicon in Buchnera aphidicola

Buchnera aphidicola is an obligate, strictly vertically transmitted, bacterial symbiont of aphids. It supplies its host with essential amino acids, nutrients required by aphids but deficient in their diet of plant phloem sap. Several lineages of Buchnera show adaptation to their nutritional role in the form of plasmid-mediated amplification of key-genes involved in the biosynthesis of tryptophan (trpEG) and leucine (leuABCD). Phylogenetic analyses of these plasmid-encoded functions have thus far suggested the absence of horizontal plasmid exchange among lineages of Buchnera. Here, we describe three new Buchnera plasmids, obtained from species of the aphid host families Lachnidae and Pemphigidae. All three plasmids belong to the repA1 family of Buchnera plasmids, which is characterized by the presence of a repA1-replicon responsible for replication initiation. A comprehensive analysis of this family of plasmids unexpectedly revealed significantly incongruent phylogenies for different plasmid and chromosomally encoded loci. We infer from these incongruencies a case of horizontal plasmid transfer in Buchnera. This process may have been mediated by secondary endosymbionts, which occasionally undergo horizontal transmission in aphids.