Genetic addiction: selfish gene's strategy for symbiosis in the genome

Spatial structure and benefits to hosts allow plasmids with and without post-segregational killing systems to coexist

To persist, a plasmid relies on being passed on to a daughter cell, but this does not always occur. Plasmids with post-segregational killing (PSK) systems kill a daughter cell if the plasmid has not been passed on. By killing the host, it also kills competing plasmids in the same host, something competing plasmids without a similar system cannot do. Accordingly, plasmids with PSK systems can displace other plasmids. In nature, plasmids with and without PSK systems coexist and prior theory has suggested this is expected to be very rare or unstable, such that one or the other type of plasmid eventually takes over. Here, we show that if there is spatial structure and plasmids confer benefits to hosts, coexistence of plasmids occurs broadly. Often plasmids confer benefits (even ones with a PSK system) and bacteria are often spatially structured. So, our results may be generally applicable.

Segregational instability of multicopy plasmids: A population genetics approach

Plasmids are extra-chromosomal genetic elements that encode a wide variety of phenotypes and can be maintained in bacterial populations through vertical and horizontal transmission, thus increasing bacterial adaptation to hostile environmental conditions like those imposed by antimicrobial substances. To circumvent the segregational instability resulting from randomly distributing plasmids between daughter cells upon division, nontransmissible plasmids tend to be carried in multiple copies per cell, with the added benefit of exhibiting increased gene dosage and resistance levels. But carrying multiple copies also results in a high metabolic burden to the bacterial host, therefore reducing the overall fitness of the population. This trade-off poses an existential question for plasmids: What is the optimal plasmid copy number? In this manuscript, we address this question by postulating and analyzing a population genetics model to evaluate the interaction between selective pressure, the number of plasmid copies carried by each cell, and the metabolic burden associated with plasmid bearing in the absence of selection for plasmid-encoded traits. Parameter values of the model were estimated experimentally using Escherichia coli K12 carrying a multicopy plasmid encoding for a fluorescent protein and bla _TEM-1, a gene conferring resistance to β-lactam antibiotics. By numerically determining the optimal plasmid copy number for constant and fluctuating selection regimes, we show that plasmid copy number is a highly optimized evolutionary trait that depends on the rate of environmental fluctuation and balances the benefit between increased stability in the absence of selection with the burden associated with carrying multiple copies of the plasmid.

Landscape of mobile genetic elements and their antibiotic resistance cargo in prokaryotic genomes

Prokaryotic Mobile Genetic Elements (MGEs) such as transposons, integrons, phages and plasmids, play important roles in prokaryotic evolution and in the dispersal of cargo functions like antibiotic resistance. However, each of these MGE types is usually annotated and analysed individually, hampering a global understanding of phylogenetic and environmental patterns of MGE dispersal. We thus developed a computational framework that captures diverse MGE types, their cargos and MGE-mediated horizontal transfer events, using recombinases as ubiquitous MGE marker genes and pangenome information for MGE boundary estimation. Applied to ∼84k genomes with habitat annotation, we mapped 2.8 million MGE-specific recombinases to six operational MGE types, which together contain on average 13% of all the genes in a genome. Transposable elements (TEs) dominated across all taxa (∼1.7 million occurrences), outnumbering phages and phage-like elements (<0.4 million). We recorded numerous MGE-mediated horizontal transfer events across diverse phyla and habitats involving all MGE types, disentangled and quantified the extent of hitchhiking of TEs (17%) and integrons (63%) with other MGE categories, and established TEs as dominant carriers of antibiotic resistance genes. We integrated all these findings into a resource (proMGE.embl.de), which should facilitate future studies on the large mobile part of genomes and its horizontal dispersal.

Mathematical Models of Plasmid Population Dynamics

With plasmid-mediated antibiotic resistance thriving and threatening to become a serious public health problem, it is paramount to increase our understanding of the forces that enable the spread and maintenance of drug resistance genes encoded in mobile genetic elements. The relevance of plasmids as vehicles for the dissemination of antibiotic resistance genes, in addition to the extensive use of plasmid-derived vectors for biotechnological and industrial purposes, has promoted the in-depth study of the molecular mechanisms controlling multiple aspects of a plasmids' life cycle. This body of experimental work has been paralleled by the development of a wealth of mathematical models aimed at understanding the interplay between transmission, replication, and segregation, as well as their consequences in the ecological and evolutionary dynamics of plasmid-bearing bacterial populations. In this review, we discuss theoretical models of plasmid dynamics that span from the molecular mechanisms of plasmid partition and copy-number control occurring at a cellular level, to their consequences in the population dynamics of complex microbial communities. We conclude by discussing future directions for this exciting research topic.

Low-level expression of the Type II restriction-modification system confers potent bacteriophage resistance in Escherichia coli

Restriction–modification systems (R–M) are one of the antiviral defense tools used by bacteria, and those of the Type II family are composed of a restriction endonuclease (REase) and a DNA methyltransferase (MTase). Most entering DNA molecules are usually cleaved by the REase before they can be methylated by MTase, although the observed level of fragmented DNA may vary significantly. Using a model EcoRI R–M system, we report that the balance between DNA methylation and cleavage may be severely affected by transcriptional signals coming from outside the R–M operon. By modulating the activity of the promoter, we obtained a broad range of restriction phenotypes for the EcoRI R–M system that differed by up to 4 orders of magnitude in our biological assays. Surprisingly, we found that high expression levels of the R–M proteins were associated with reduced restriction of invading bacteriophage DNA. Our results suggested that the regulatory balance of cleavage and methylation was highly sensitive to fluctuations in transcriptional signals both up- and downstream of the R–M operon. Our data provided further insights into Type II R–M system maintenance and the potential conflict within the host bacterium.

'Bacterial Programmed Cell Death': cellular altruism or genetic selfism?

Cell-dependent propagation of the 'self' is the driver of all species, organisms and even genes. Conceivably, elimination of these entities is caused by cellular death. Then, how can genes that cause the death of the same cell evolve? Programmed cell death (PCD) is the gene-dependent self-inflicted death. In multicellular organisms, PCD of a cell confers fitness to the surviving rest of the organism, which thereby allows the selection of genes responsible for PCD. However, PCD in free-living bacteria is intriguing; the death of the cell is the death of the organism. How can such PCD genes be selected in unicellular organisms? The bacterial PCD in a population is proposed to confer fitness to the surviving kin in the form of sporulation, nutrition, infection-containment and matrix materials. While the cell-centred view leading to propositions of 'altruism' is enticing, the gene-centred view of 'selfism' is neglected. In this opinion piece, we reconceptualize the PCD propositions as genetic selfism (death due to loss/mutation of selfish genes) rather than cellular altruism (death for the conferment of fitness to kin). Within the scope and the available evidence, we opine that some of the PCD-like observations in bacteria seem to be the manifestation of genetic selfism by Restriction-Modification systems and Toxin-Antitoxin systems.

Quantifying plasmid dynamics using single-cell microfluidics and image bioinformatics

Multicopy plasmids play an important role in bacterial ecology and evolution by accelerating the rate of adaptation and providing a platform for rapid gene amplification and evolutionary rescue. Despite the relevance of plasmids in bacterial evolutionary dynamics, evaluating the population-level consequences of randomly segregating and replicating plasmids in individual cells remains a challenging problem, both in theory and experimentally. In recent years, technological advances in fluorescence microscopy and microfluidics have allowed studying temporal changes in gene expression by quantifying the fluorescent intensity of individual cells under controlled environmental conditions. In this paper, we will describe the manufacture, experimental setup, and data analysis pipeline of different microfluidic systems that can be used to study plasmid dynamics, both in single-cells and in populations. To illustrate the benefits and limitations of microfluidics to study multicopy plasmid dynamics, we will use an experimental model system consisting on Escherichia coli K12 carrying non-conjugative, multicopy plasmids (19 copies per cell, in average) encoding different fluorescent markers and β-lactam resistance genes. First, we will use an image-based flow cytometer to estimate changes in the allele distribution of a heterogeneous population under different selection regimes. Then we will use a mothermachine microfluidic device to obtain time-series of fluorescent intensity of individual cells to argue that plasmid segregation and replication dynamics are inherently stochastic processes. Finally, using a microchemostat, we track thousands of cells in time to reconstruct bacterial lineages and evaluate the allele frequency distributions that emerge in response to a range of selective pressures.

Genome-inferred spatio-temporal resolution of an uncultivated Roizmanbacterium reveals its ecological preferences in groundwater

Subsurface ecosystems like groundwater harbour diverse microbial communities, including small-sized, putatively symbiotic organisms of the Candidate Phyla Radiation, yet little is known about their ecological preferences and potential microbial partners. Here, we investigated a member of the superphylum Microgenomates (Cand. Roizmanbacterium ADI133) from oligotrophic groundwater using mini-metagenomics and monitored its spatio-temporal distribution using 16S rRNA gene analyses. A Roizmanbacteria-specific quantitative PCR assay allowed us to track its abundance over the course of 1 year within eight groundwater wells along a 5.4 km hillslope transect, where Roizmanbacteria reached maximum relative abundances of 2.3%. In-depth genomic analyses suggested that Cand. Roizmanbacterium ADI133 is a lactic acid fermenter, potentially able to utilize a range of complex carbon substrates, including cellulose. We hypothesize that it attaches to host cells using a trimeric autotransporter adhesin and inhibits their cell wall biosynthesis using a toxin-antitoxin system. Network analyses based on correlating Cand. Roizmanbacterium ADI133 abundances with amplicon sequencing-derived microbial community profiles suggested one potential host organism, classified as a member of the class Thermodesulfovibrionia (Nitrospirae). By providing lactate as an electron donor Cand. Roizmanbacterium ADI133 potentially mediates the transfer of carbon to other microorganisms and thereby is an important connector in the microbial community.

Evolutionary entanglement of mobile genetic elements and host defence systems: guns for hire

All cellular life forms are afflicted by diverse genetic parasites, including viruses and other types of mobile genetic elements (MGEs), and have evolved multiple, diverse defence systems that protect them from MGE assault via different mechanisms. Here, we provide our perspectives on how recent evidence points to tight evolutionary connections between MGEs and defence systems that reach far beyond the proverbial arms race. Defence systems incur a fitness cost for the hosts; therefore, at least in prokaryotes, horizontal mobility of defence systems, mediated primarily by MGEs, is essential for their persistence. Moreover, defence systems themselves possess certain features of selfish elements. Common components of MGEs, such as site-specific nucleases, are 'guns for hire' that can also function as parts of defence mechanisms and are often shuttled between MGEs and defence systems. Thus, evolutionary and molecular factors converge to mould the multifaceted, inextricable connection between MGEs and anti-MGE defence systems.

In vitro exploration of the Xanthomonas hortorum pv. vitians genome using transposon insertion sequencing and comparative genomics to discriminate between core and contextual essential genes

The essential genome of a bacterium encompasses core genes associated with basic cellular processes and conditionally essential genes dependent upon environmental conditions or the genetic context. Comprehensive knowledge of those gene sets allows for a better understanding of fundamental bacterial biology and offers new perspectives for antimicrobial drug research against detrimental bacteria such as pathogens. We investigated the essential genome of Xanthomonas hortorum pv. vitians, a gammaproteobacterial plant pathogen of lettuce (Lactuca sativa L.) which belongs to the plant-pathogen reservoir genus Xanthomonas and is affiliated to the family Xanthomonadaceae. No practical means of disease control or prevention against this pathogen is currently available, and its molecular biology is virtually unknown. To reach a comprehensive overview of the essential genome of X. hortorum pv. vitians LM16734, we developed a mixed approach combining high-quality full genome sequencing, saturated transposon insertion sequencing (Tn-Seq) in optimal growth conditions, and coupled computational analyses such as comparative genomics, synteny assessment and phylogenomics. Among the 370 essential loci identified by Tn-Seq, a majority was bound to critical cell processes conserved across bacteria. The remaining genes were either related to specific ecological features of Xanthomonas or Xanthomonadaceae species, or acquired through horizontal gene transfer of mobile genetic elements and associated with ancestral parasitic gene behaviour and bacterial defence systems. Our study sheds new light on our usual concepts about gene essentiality and is pioneering in the molecular and genomic study of X. hortorum pv. vitians.

Evolution of the Quorum network and the mobilome (plasmids and bacteriophages) in clinical strains of Acinetobacter baumannii during a decade

In this study, we compared eighteen clinical strains of A. baumannii belonging to the ST-2 clone and isolated from patients in the same intensive care unit (ICU) in 2000 (9 strains referred to collectively as Ab_GEIH-2000) and 2010 (9 strains referred to collectively as Ab_GEIH-2010), during the GEIH-REIPI project (Umbrella BioProject PRJNA422585). We observed two main molecular differences between the Ab_GEIH-2010 and the Ab_GEIH-2000 collections, acquired over the course of the decade long sampling interval and involving the mobilome: i) a plasmid harbouring genes for bla_{OXA 24/40} ß-lactamase and abKA/abkB proteins of a toxin-antitoxin system; and ii) two temperate bacteriophages, Ab105-1ϕ (63 proteins) and Ab105-2ϕ (93 proteins), containing important viral defence proteins. Moreover, all Ab_GEIH-2010 strains contained a Quorum functional network of Quorum Sensing (QS) and Quorum Quenching (QQ) mechanisms, including a new QQ enzyme, AidA, which acts as a bacterial defence mechanism against the exogenous 3-oxo-C12-HSL. Interestingly, the infective capacity of the bacteriophages isolated in this study (Ab105-1ϕ and Ab105-2ϕ) was higher in the Ab_GEIH-2010 strains (carrying a functional Quorum network) than in the Ab_GEIH-2000 strains (carrying a deficient Quorum network), in which the bacteriophages showed little or no infectivity. This is the first study about the evolution of the Quorum network and the mobilome in clinical strains of Acinetobacter baumannii during a decade.

Programmed cell death can increase the efficacy of microbial bet -hedging

Programmed cell death (PCD) occurs in both unicellular and multicellular organisms. While PCD plays a key role in the development and maintenance of multicellular organisms, explaining why single-celled organisms would evolve to actively commit suicide has been far more challenging. Here, we explore the potential for PCD to act as an accessory to microbial bet-hedging strategies that utilize stochastic phenotype switching. We consider organisms that face unpredictable and recurring disasters, in which fitness depends on effective phenotypic diversification. We show that when reproductive opportunities are limited by carrying capacity, PCD drives population turnover, providing increased opportunities for phenotypic diversification through stochastic phenotype switching. The main cost of PCD, providing resources for growth to a PCD(−) competitor, is ameliorated by genetic assortment in spatially structured populations. Using agent -based simulations, we explore how basic demographic factors, namely bottlenecks and local dispersal, can generate sufficient spatial structure to favor the evolution of high PCD rates.

Viruses and mobile elements as drivers of evolutionary transitions

The history of life is punctuated by evolutionary transitions which engender emergence of new levels of biological organization that involves selection acting at increasingly complex ensembles of biological entities. Major evolutionary transitions include the origin of prokaryotic and then eukaryotic cells, multicellular organisms and eusocial animals. All or nearly all cellular life forms are hosts to diverse selfish genetic elements with various levels of autonomy including plasmids, transposons and viruses. I present evidence that, at least up to and including the origin of multicellularity, evolutionary transitions are driven by the coevolution of hosts with these genetic parasites along with sharing of ‘public goods’. Selfish elements drive evolutionary transitions at two distinct levels. First, mathematical modelling of evolutionary processes, such as evolution of primitive replicator populations or unicellular organisms, indicates that only increasing organizational complexity, e.g. emergence of multicellular aggregates, can prevent the collapse of the host–parasite system under the pressure of parasites. Second, comparative genomic analysis reveals numerous cases of recruitment of genes with essential functions in cellular life forms, including those that enable evolutionary transitions.

This article is part of the themed issue ‘The major synthetic evolutionary transitions’.

The third restriction-modification system from Thermus aquaticus YT-1: solving the riddle of two TaqII specificities

Two restriction-modification systems have been previously discovered in Thermus aquaticus YT-1. TaqI is a 263-amino acid (aa) Type IIP restriction enzyme that recognizes and cleaves within the symmetric sequence 5'-TCGA-3'. TaqII, in contrast, is a 1105-aa Type IIC restriction-and-modification enzyme, one of a family of Thermus homologs. TaqII was originally reported to recognize two different asymmetric sequences: 5'-GACCGA-3' and 5'-CACCCA-3'. We previously cloned the taqIIRM gene, purified the recombinant protein from Escherichia coli, and showed that TaqII recognizes the 5'-GACCGA-3' sequence only. Here, we report the discovery, isolation, and characterization of TaqIII, the third R-M system from T. aquaticus YT-1. TaqIII is a 1101-aa Type IIC/IIL enzyme and recognizes the 5'-CACCCA-3' sequence previously attributed to TaqII. The cleavage site is 11/9 nucleotides downstream of the A residue. The enzyme exhibits striking biochemical similarity to TaqII. The 93% identity between their aa sequences suggests that they have a common evolutionary origin. The genes are located on two separate plasmids, and are probably paralogs or pseudoparalogs. Putative positions and aa that specify DNA recognition were identified and recognition motifs for 6 uncharacterized Thermus-family enzymes were predicted.

Novel toxins from type II toxin-antitoxin systems with acetyltransferase activity

Type II toxin-antitoxin (TA) systems are widespread in bacterial and archeal genomes. These modules are very dynamic and participate in bacterial genome evolution through horizontal gene transfer. TA systems are commonly composed of a labile antitoxin and a stable toxin. Toxins appear to preferentially inhibit the protein synthesis process. Toxins use a variety of molecular mechanisms and target nearly every step of translation to achieve their inhibitory function. This review focuses on a recently identified TA family that includes acetyltransferase toxins. The AtaT and TacT toxins are the best-characterized to date in this family. AtaT and TacT both inhibit translation by acetylating the amino acid charged on tRNAs. However, the specificities of these 2 toxins are different as AtaT inhibits translation initiation by acetylation of the initiator tRNA whereas TacT acetylates elongator tRNAs. The molecular mechanisms of these toxins are discussed, as well as the functions and possible evolutionary origins of this diverse toxin family.

The Evolution of Antiviral Defense Systems

Self-replicating genetic material presumably provided the architecture necessary for generating the last universal ancestor of all nucleic-acid-based life. As biological complexity increased in the billions of years that followed, the same genetic material also morphed into a wide spectrum of viruses and other parasitic genetic elements. The resulting struggle for existence drove the evolution of host defenses, giving rise to a perpetual arms race. This Perspective summarizes the antiviral mechanisms evident across the tree of life, discussing each in their evolutionary context to postulate how the coevolution of host and pathogen shaped the cellular antiviral defenses we know today.

Evolutionary Ecology of Prokaryotic Immune Mechanisms

Bacteria have a range of distinct immune strategies that provide protection against bacteriophage (phage) infections. While much has been learned about the mechanism of action of these defense strategies, it is less clear why such diversity in defense strategies has evolved. In this review, we discuss the short- and long-term costs and benefits of the different resistance strategies and, hence, the ecological conditions that are likely to favor the different strategies alone and in combination. Finally, we discuss some of the broader consequences, beyond resistance to phage and other genetic elements, resulting from the operation of different immune strategies.

Toxin-Antitoxin Systems in Clinical Pathogens

Toxin-antitoxin (TA) systems are prevalent in bacteria and archaea. Although not essential for normal cell growth, TA systems are implicated in multiple cellular functions associated with survival under stress conditions. Clinical strains of bacteria are currently causing major human health problems as a result of their multidrug resistance, persistence and strong pathogenicity. Here, we present a review of the TA systems described to date and their biological role in human pathogens belonging to the ESKAPE group (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa and Enterobacter spp.) and others of clinical relevance (Escherichia coli, Burkholderia spp., Streptococcus spp. and Mycobacterium tuberculosis). Better understanding of the mechanisms of action of TA systems will enable the development of new lines of treatment for infections caused by the above-mentioned pathogens.

PsasM2I, a type II restriction-modification system in Pseudomonas savastanoi pv. savastanoi: differential distribution of carrier strains in the environment and the evolutionary history of homologous RM systems in the Pseudomonas syringae complex

A type II restriction-modification system was found in a native plasmid of Pseudomonas savastanoi pv. savastanoi MLLI2. Functional analysis of the methyltransferase showed that the enzyme acts by protecting the DNA sequence CTGCAG from cleavage. Restriction endonuclease expression in recombinant Escherichia coli cells resulted in mutations in the REase sequence or transposition of insertion sequence 1A in the coding sequence, preventing lethal gene expression. Population screening detected homologous RM systems in other P. savastanoi strains and in the Pseudomonas syringae complex. An epidemiological survey carried out by sampling olive and oleander knots in two Italian regions showed an uneven diffusion of carrier strains, whose presence could be related to a selective advantage in maintaining the RM system in particular environments or subpopulations. Moreover, carrier strains can coexist in the same orchards, plants, and knot tissues with non-carriers, revealing unexpected genetic variability on a very small spatial scale. Phylogenetic analysis of the RM system and housekeeping gene sequences in the P. syringae complex demonstrated the ancient acquisition of the RM systems. However, the evolutionary history of the gene complex also showed the involvement of horizontal gene transfer between related strains and recombination events.

The interplay of restriction-modification systems with mobile genetic elements and their prokaryotic hosts

The roles of restriction-modification (R-M) systems in providing immunity against horizontal gene transfer (HGT) and in stabilizing mobile genetic elements (MGEs) have been much debated. However, few studies have precisely addressed the distribution of these systems in light of HGT, its mechanisms and its vectors. We analyzed the distribution of R-M systems in 2261 prokaryote genomes and found their frequency to be strongly dependent on the presence of MGEs, CRISPR-Cas systems, integrons and natural transformation. Yet R-M systems are rare in plasmids, in prophages and nearly absent from other phages. Their abundance depends on genome size for small genomes where it relates with HGT but saturates at two occurrences per genome. Chromosomal R-M systems might evolve under cycles of purifying and relaxed selection, where sequence conservation depends on the biochemical activity and complexity of the system and total gene loss is frequent. Surprisingly, analysis of 43 pan-genomes suggests that solitary R-M genes rarely arise from the degradation of R-M systems. Solitary genes are transferred by large MGEs, whereas complete systems are more frequently transferred autonomously or in small MGEs. Our results suggest means of testing the roles for R-M systems and their associations with MGEs.

Characterization of plasmids carrying the blaOXA-24/40 carbapenemase gene and the genes encoding the AbkA/AbkB proteins of a toxin/antitoxin system

BACKGROUND

Carbapenem-resistant Acinetobacter baumannii (CRAb) is a major source of nosocomial infections in Spain associated with the production of OXA-58-like or OXA-24/40-like β-lactamase enzymes. We analysed the plasmids carrying the bla(OXA-24/40)-like gene in CRAb isolates obtained a decade apart.

METHODS

The presence of β-lactamases was screened for by PCR (metallo-β-lactamases, carbapenem-hydrolysing class D β-lactamases, GES and KPC) in 101 CRAb isolates obtained in two multicentre studies (GEIH/REIPI-Ab-2000 and GEIH/REIPI-Ab-2010; n = 493 Acinetobacter spp). We analysed the distribution and characterization of the plasmids carrying the bla(OXA-24/40)-like gene and sequenced two plasmids, AbATCC223p (2000) and AbATCC329p (2010) from A. baumannii ATCC 17978 transformants.

RESULTS

Acquisition of the bla(OXA-24/40)-like gene was the main mechanism underlying resistance to carbapenems (48.7% in 2000 compared with 51.6% in 2010). This gene was mainly isolated in ST2 A. baumannii strains in both studies, although some novel STs (ST79 and ST80) appeared in 2010. The gene was located in plasmids (8-12 kbp) associated with the repAci2 or repAci2/repGR12 types. The sequences of AbATCC223p (8840 bp) and AbATCC329p (8842 bp) plasmids were similar, particularly regarding the presence of the genes encoding the AbkA/AbkB proteins associated with the toxin/antitoxin system. Moreover, the abkA/abkB gene sequences (>96% identity) were also located in plasmids harbouring the bla(OXA-58)-like gene.

CONCLUSIONS

The action of OXA-24/40 and OXA-58 β-lactamase-like enzymes represents the main mechanism underlying resistance to carbapenems in Spain in the last decade. AbkA/AbkB proteins in the toxin/antitoxin system may be involved in the successful dissemination of plasmids carrying the bla(OXA-24/40)-like gene, and probably also the bla(OXA-58)-like gene, thus contributing to the plasmid stability.

Highlights of the DNA cutters: a short history of the restriction enzymes

In the early 1950's, 'host-controlled variation in bacterial viruses' was reported as a non-hereditary phenomenon: one cycle of viral growth on certain bacterial hosts affected the ability of progeny virus to grow on other hosts by either restricting or enlarging their host range. Unlike mutation, this change was reversible, and one cycle of growth in the previous host returned the virus to its original form. These simple observations heralded the discovery of the endonuclease and methyltransferase activities of what are now termed Type I, II, III and IV DNA restriction-modification systems. The Type II restriction enzymes (e.g. EcoRI) gave rise to recombinant DNA technology that has transformed molecular biology and medicine. This review traces the discovery of restriction enzymes and their continuing impact on molecular biology and medicine.

To be or not to be: regulation of restriction-modification systems and other toxin-antitoxin systems

One of the simplest classes of genes involved in programmed death is that containing the toxin–antitoxin (TA) systems of prokaryotes. These systems are composed of an intracellular toxin and an antitoxin that neutralizes its effect. These systems, now classified into five types, were initially discovered because some of them allow the stable maintenance of mobile genetic elements in a microbial population through postsegregational killing or the death of cells that have lost these systems. Here, we demonstrate parallels between some TA systems and restriction–modification systems (RM systems). RM systems are composed of a restriction enzyme (toxin) and a modification enzyme (antitoxin) and limit the genetic flux between lineages with different epigenetic identities, as defined by sequence-specific DNA methylation. The similarities between these systems include their postsegregational killing and their effects on global gene expression. Both require the finely regulated expression of a toxin and antitoxin. The antitoxin (modification enzyme) or linked protein may act as a transcriptional regulator. A regulatory antisense RNA recently identified in an RM system can be compared with those RNAs in TA systems. This review is intended to generalize the concept of TA systems in studies of stress responses, programmed death, genetic conflict and epigenetics.

Diverse functions of restriction-modification systems in addition to cellular defense

Restriction-modification (R-M) systems are ubiquitous and are often considered primitive immune systems in bacteria. Their diversity and prevalence across the prokaryotic kingdom are an indication of their success as a defense mechanism against invading genomes. However, their cellular defense function does not adequately explain the basis for their immaculate specificity in sequence recognition and nonuniform distribution, ranging from none to too many, in diverse species. The present review deals with new developments which provide insights into the roles of these enzymes in other aspects of cellular function. In this review, emphasis is placed on novel hypotheses and various findings that have not yet been dealt with in a critical review. Emerging studies indicate their role in various cellular processes other than host defense, virulence, and even controlling the rate of evolution of the organism. We also discuss how R-M systems could have successfully evolved and be involved in additional cellular portfolios, thereby increasing the relative fitness of their hosts in the population.

Cryptic Streptococcus mutans 5.6-kb plasmids encode a toxin-antitoxin system for plasmid stabilization

BACKGROUND

In all Streptococcus mutans strains, 5-13% carry a 5.6-kb plasmid. Despite its frequency, little is known about its mediated functions with most of the information coming from a single study focussing on plasmid pUA140.

OBJECTIVE

Here, we describe the sequence and genetic organization of two S. mutans 5.6-kb plasmids, pDC09 and pNC101.

RESULTS

Based on PicoGreen dsDNA quantification and Real-Time quantitative PCR (RTQ-PCR), the plasmid copy number was found to range between 10 and 74, depending on the strain tested. In contrast to literature, we identified six instead of five open reading frames (ORFs). While the putative gene products of ORF1 (as a Rep-protein) and ORF2 (as a Mob-protein) could be confirmed as being identical to those from pUA140, the functions of ORF3 (unknown) and ORF 4 (possibly AtpE homologue) could not be further revealed. However, the product of ORF5 showed a fairly high identity (38-50%) and structural similarity (58-74%) to RelE of Streptococcus pneumoniae, Streptococcus equi, and Streptococcus downei. In addition, we identified a functionally corresponding ORF6 encoding a protein with 61-68% identity (81-86% similarity) to the S. equi and S. downei antitoxin of the RelB family. RelE and RelB together form a plasmid-encoded toxin-antitoxin (TA) system, RelBE(plas). Despite its rather limited sequence similarity with chromosomal TA systems in S. mutans (RelBE(chro), MazEF, HicBA), we found similar tertiary structures applying I-Tasser protein prediction analysis.

CONCLUSION

Type II-toxins, as the plasmid-encoded RelE, are RNA endonucleases. Depending on their mRNA cleavage activity, they might 1) kill every plasmid-free progeny, thereby stabilizing plasmid transfer at the expense of the host and/or 2) help S. mutans enter a dormant state and survive unfavourable environmental conditions. Whilst a function in plasmid stabilization has been confirmed, a function in persistence under nutritional stress, tested here by inducing amino acid starvation, could not be demonstrated so far.

The coevolution of toxin and antitoxin genes drives the dynamics of bacterial addiction complexes and intragenomic conflict

Bacterial genomes commonly contain 'addiction' gene complexes that code for both a toxin and a corresponding antitoxin. As long as both genes are expressed, cells carrying the complex can remain healthy. However, loss of the complex (including segregational loss in daughter cells) can entail death of the cell. We develop a theoretical model to explore a number of evolutionary puzzles posed by toxin-antitoxin (TA) population biology. We first extend earlier results demonstrating that TA complexes can spread on plasmids, as an adaptation to plasmid competition in spatially structured environments, and highlight the role of kin selection. We then considered the emergence of TA complexes on plasmids from previously unlinked toxin and antitoxin genes. We find that one of these traits must offer at least initially a direct advantage in some but not all environments encountered by the evolving plasmid population. Finally, our study predicts non-transitive 'rock-paper-scissors' dynamics to be a feature of intragenomic conflict mediated by TA complexes. Intragenomic conflict could be sufficient to select deleterious genes on chromosomes and helps to explain the previously perplexing observation that many TA genes are found on bacterial chromosomes.

Movement of DNA sequence recognition domains between non-orthologous proteins

Comparisons of proteins show that they evolve through the movement of domains. However, in many cases, the underlying mechanisms remain unclear. Here, we observed the movements of DNA recognition domains between non-orthologous proteins within a prokaryote genome. Restriction-modification (RM) systems, consisting of a sequence-specific DNA methyltransferase and a restriction enzyme, contribute to maintenance/evolution of genomes/epigenomes. RM systems limit horizontal gene transfer but are themselves mobile. We compared Type III RM systems in Helicobacter pylori genomes and found that target recognition domain (TRD) sequences are mobile, moving between different orthologous groups that occupy unique chromosomal locations. Sequence comparisons suggested that a likely underlying mechanism is movement through homologous recombination of similar DNA sequences that encode amino acid sequence motifs that are conserved among Type III DNA methyltransferases. Consistent with this movement, incongruence was observed between the phylogenetic trees of TRD regions and other regions in proteins. Horizontal acquisition of diverse TRD sequences was suggested by detection of homologs in other Helicobacter species and distantly related bacterial species. One of these RM systems in H. pylori was inactivated by insertion of another RM system that likely transferred from an oral bacterium. TRD movement represents a novel route for diversification of DNA-interacting proteins.

Success of a suicidal defense strategy against infection in a structured habitat

Pathogen infection often leads to the expression of virulence and host death when the host-pathogen symbiosis seems more beneficial for the pathogen. Previously proposed explanations have focused on the pathogen's side. In this work, we tested a hypothesis focused on the host strategy. If a member of a host population dies immediately upon infection aborting pathogen reproduction, it can protect the host population from secondary infections. We tested this "Suicidal Defense Against Infection" (SDAI) hypothesis by developing an experimental infection system that involves a huge number of bacteria as hosts and their virus as pathogen, which is linked to modeling and simulation. Our experiments and simulations demonstrate that a population with SDAI strategy is successful in the presence of spatial structure but fails in its absence. The infection results in emergence of pathogen mutants not inducing the host suicide in addition to host mutants resistant to the pathogen.

Domain movement within a gene: a novel evolutionary mechanism for protein diversification

A protein function is carried out by a specific domain localized at a specific position. In the present study, we report that, within a gene, a specific amino acid sequence can move between a certain position and another position. This was discovered when the sequences of restriction-modification systems within the bacterial species Helicobacter pylori were compared. In the specificity subunit of Type I restriction-modification systems, DNA sequence recognition is mediated by target recognition domain 1 (TRD1) and TRD2. To our surprise, several sequences are shared by TRD1 and TRD2 of genes (alleles) at the same locus (chromosomal location); these domains appear to have moved between the two positions. The gene/protein organization can be represented as x-(TRD1)-y-x-(TRD2)-y, where x and y represent repeat sequences. Movement probably occurs by recombination at these flanking DNA repeats. In accordance with this hypothesis, recombination at these repeats also appears to decrease two TRDs into one TRD or increase these two TRDs to three TRDs (TRD1-TRD2-TRD2) and to allow TRD movement between genes even at different loci. Similar movement of domains between TRD1 and TRD2 was observed for the specificity subunit of a Type IIG restriction enzyme. Similar movement of domain between TRD1 and TRD2 was observed for Type I restriction-modification enzyme specificity genes in two more eubacterial species, Streptococcus pyogenes and Mycoplasma agalactiae. Lateral domain movements within a protein, which we have designated DOMO (domain movement), represent novel routes for the diversification of proteins.

Antisense RNA associated with biological regulation of a restriction-modification system

Restriction–modification systems consist of a modification enzyme that methylates a specific DNA sequence and a restriction endonuclease that cleaves DNA lacking this epigenetic signature. Their gene expression should be finely regulated because their potential to attack the host bacterial genome needs to be controlled. In the EcoRI system, where the restriction gene is located upstream of the modification gene in the same orientation, we previously identified intragenic reverse promoters affecting gene expression. In the present work, we identified a small (88 nt) antisense RNA (Rna0) transcribed from a reverse promoter (P_REV0) at the 3′ end of the restriction gene. Its antisense transcription, as measured by transcriptional gene fusion, appeared to be terminated by the P_M1,M2 promoter. P_M1,M2 promoter-initiated transcription, in turn, appeared to be inhibited by P_REV0. Mutational inactivation of P_REV0 increased expression of the restriction gene. The biological significance of this antisense transcription is 2-fold. First, a mutation in P_REV0 increased restriction of incoming DNA. Second, the presence of the antisense RNA gene (ecoRIA) in trans alleviated cell killing after loss of the EcoRI plasmid (post-segregational killing). Taken together, these results strongly suggested the involvement of an antisense RNA in the biological regulation of this restriction–modification system.

IS-linked movement of a restriction-modification system

Potential mobility of restriction-modification systems has been suggested by evolutionary/bioinformatic analysis of prokaryotic genomes. Here we demonstrate in vivo movement of a restriction-modification system within a genome under a laboratory condition. After blocking replication of a temperature-sensitive plasmid carrying a PaeR7I restriction-modification system in Escherichia coli cells, the plasmid was found integrated into the chromosome of the surviving cells. Sequence analysis revealed that, in the majority of products, the restriction-modification system was linked to chromosomal insertion sequences (ISs). Three types of products were: (I) apparent co-integration of the plasmid and the chromosome at a chromosomal IS1 or IS5 copy (24/28 analyzed); (II) de novo insertion of IS1 with the entire plasmid except for a 1–3 bp terminal deletion (2/28); and (III) reciprocal crossing-over between the plasmid and the chromosome involving 1–3 bp of sequence identity (2/28). An R-negative mutation apparently decreased the efficiency of successful integration by two orders of magnitude. Reconstruction experiments demonstrated that the restriction-dependence was mainly due to selection against cells without proper integration: their growth was inhibited by the restriction enzyme action. These results demonstrate collaboration of a mobile element and a restriction-modification system for successful joint migration. This collaboration may have promoted the spread and, therefore, the long-term persistence of these complexes and restriction-modification systems in a wide range of prokaryotes.

What traits are carried on mobile genetic elements, and why?

Although similar to any other organism, prokaryotes can transfer genes vertically from mother cell to daughter cell, they can also exchange certain genes horizontally. Genes can move within and between genomes at fast rates because of mobile genetic elements (MGEs). Although mobile elements are fundamentally self-interested entities, and thus replicate for their own gain, they frequently carry genes beneficial for their hosts and/or the neighbours of their hosts. Many genes that are carried by mobile elements code for traits that are expressed outside of the cell. Such traits are involved in bacterial sociality, such as the production of public goods, which benefit a cell's neighbours, or the production of bacteriocins, which harm a cell's neighbours. In this study we review the patterns that are emerging in the types of genes carried by mobile elements, and discuss the evolutionary and ecological conditions under which mobile elements evolve to carry their peculiar mix of parasitic, beneficial and cooperative genes.

Conflicts targeting epigenetic systems and their resolution by cell death: novel concepts for methyl-specific and other restriction systems

Epigenetic modification of genomic DNA by methylation is important for defining the epigenome and the transcriptome in eukaryotes as well as in prokaryotes. In prokaryotes, the DNA methyltransferase genes often vary, are mobile, and are paired with the gene for a restriction enzyme. Decrease in a certain epigenetic methylation may lead to chromosome cleavage by the partner restriction enzyme, leading to eventual cell death. Thus, the pairing of a DNA methyltransferase and a restriction enzyme forces an epigenetic state to be maintained within the genome. Although restriction enzymes were originally discovered for their ability to attack invading DNAs, it may be understood because such DNAs show deviation from this epigenetic status. DNAs with epigenetic methylation, by a methyltransferase linked or unlinked with a restriction enzyme, can also be the target of DNases, such as McrBC of Escherichia coli, which was discovered because of its methyl-specific restriction. McrBC responds to specific genome methylation systems by killing the host bacterial cell through chromosome cleavage. Evolutionary and genomic analysis of McrBC homologues revealed their mobility and wide distribution in prokaryotes similar to restriction–modification systems. These findings support the hypothesis that this family of methyl-specific DNases evolved as mobile elements competing with specific genome methylation systems through host killing. These restriction systems clearly demonstrate the presence of conflicts between epigenetic systems.

Within-host competition selects for plasmid-encoded toxin-antitoxin systems

Toxin-antitoxin (TA) systems are commonly found on bacterial plasmids. The antitoxin inhibits toxin activity unless the system is lost from the cell. Then the shorter lived antitoxin degrades and the cell becomes susceptible to the toxin. Selection for plasmid-encoded TA systems was initially thought to result from their reducing the number of plasmid-free cells arising during growth in monoculture. However, modelling and experiments have shown that this mechanism can only explain the success of plasmid TA systems under a restricted set of conditions. Previously, we have proposed and tested an alternative model explaining the success of plasmid TA systems as a consequence of competition occurring between plasmids during co-infection of bacterial hosts. Here, we test a further prediction of this model, that competition between plasmids will lead to the biased accumulation of TA systems on plasmids relative to chromosomes. Transposon-encoded TA systems were added to populations of plasmid-containing cells, such that TA systems could insert into either plasmids or chromosomes. These populations were enriched for transposon-containing cells and then incubated in environments that did, or did not, allow effective within-host plasmid competition to occur. Changes in the ratio of plasmid- to chromosome-encoded TA systems were monitored. In agreement with our model, we found that plasmid-encoded TA systems had a competitive advantage, but only when host cells were sensitive to the effect of TA systems. This result demonstrates that within-host competition between plasmids can select for TA systems.

Rickettsia felis infection in a common household insect pest, Liposcelis bostrychophila (Psocoptera: Liposcelidae)

Many species of Rickettsia are well-known mammalian pathogens transmitted by blood-feeding arthropods. However, molecular surveys are continually uncovering novel Rickettsia species, often in unexpected hosts, including many arthropods that do not feed on blood. This study reports a systematic molecular characterization of a Rickettsia infecting the psocid Liposcelis bostrychophila (Psocoptera: Liposcelidae), a common and cosmopolitan household pest. Surprisingly, the psocid Rickettsia is shown to be Rickettsia felis, a human pathogen transmitted by fleas that causes serious morbidity and occasional mortality. The plasmid from the psocid R. felis was sequenced and was found to be virtually identical to the one in R. felis from fleas. As Liposcelis insects are often intimately associated with humans and other vertebrates, it is speculated that they acquired R. felis from fleas. Whether the R. felis in psocids causes disease in vertebrates is not known and warrants further study.

Genome comparison and context analysis reveals putative mobile forms of restriction-modification systems and related rearrangements

The mobility of restriction–modification (RM) gene complexes and their association with genome rearrangements is a subject of active investigation. Here we conducted systematic genome comparisons and genome context analysis on fully sequenced prokaryotic genomes to detect RM-linked genome rearrangements. RM genes were frequently found to be linked to mobility-related genes such as integrase and transposase homologs. They were flanked by direct and inverted repeats at a significantly high frequency. Insertion by long target duplication was observed for I, II, III and IV restriction types. We found several RM genes flanked by long inverted repeats, some of which had apparently inserted into a genome with a short target duplication. In some cases, only a portion of an apparently complete RM system was flanked by inverted repeats. We also found a unit composed of RM genes and an integrase homolog that integrated into a tRNA gene. An allelic substitution of a Type III system with a linked Type I and IV system pair, and allelic diversity in the putative target recognition domain of Type IIG systems were observed. This study revealed the possible mobility of all types of RM systems, and the diversity in their mobility-related organization.

Evolutionary maintenance of selfish homing endonuclease genes in the absence of horizontal transfer

Homing endonuclease genes are "selfish" mobile genetic elements whose endonuclease promotes the spread of its own gene by creating a break at a specific target site and using the host machinery to repair the break by copying and inserting the gene at this site. Horizontal transfer across the boundary of a species or population within which mating takes place has been thought to be necessary for their evolutionary persistence. This is based on the assumption that they will become fixed in a host population, where opportunities of homing will disappear, and become susceptible to degeneration. To test this hypothesis, we modeled behavior of a homing endonuclease gene that moves during meiosis through double-strand break repair. We mathematically explored conditions for persistence of the homing endonuclease gene and elucidated their parameter dependence as phase diagrams. We found that, if the cost of the pseudogene is lower than that of the homing endonuclease gene, the 2 forms can persist in a population through autonomous periodic oscillation. If the cost of the pseudogene is higher, 2 types of dynamics appear that enable evolutionary persistence: bistability dependent on initial frequency or fixation irrespective of initial frequency. The prediction of long persistence in the absence of horizontal transfer was confirmed by stochastic simulations in finite populations. The average time to extinction of the endonuclease gene was found to be thousands of meiotic generations or more based on realistic parameter values. These results provide a solid theoretical basis for an understanding of these and other extremely selfish elements.

Contribution of RecFOR machinery of homologous recombination to cell survival after loss of a restriction-modification gene complex

Loss of a type II restriction-modification (RM) gene complex, such as EcoRI, from a bacterial cell leads to death of its descendent cells through attack by residual restriction enzymes on undermethylated target sites of newly synthesized chromosomes. Through such post-segregational host killing, these gene complexes impose their maintenance on their host cells. This finding led to the rediscovery of type II RM systems as selfish mobile elements. The host prokaryote cells were found to cope with such attacks through a variety of means. The RecBCD pathway of homologous recombination in Escherichia coli repairs the lethal lesions on the chromosome, whilst it destroys restricted non-self DNA. recBCD homologues, however, appear very limited in distribution among bacterial genomes, whereas homologues of the RecFOR proteins, responsible for another pathway, are widespread in eubacteria, just like the RM systems. In the present work, therefore, we examined the possible contribution of the RecFOR pathway to cell survival after loss of an RM gene complex. A recF mutation reduced survival in an otherwise rec-positive background and, more severely, in a recBC sbcBC background. We also found that its effect is prominent in the presence of specific non-null mutant forms of the RecBCD enzyme: the resistance to killing seen with recC1002, recC1004, recC2145 and recB2154 is severely reduced to the level of a null recBC allele when combined with a recF, recO or recR mutant allele. Such resistance was also dependent on RecJ and RecQ functions. UV resistance of these non-null recBCD mutants is also reduced by recF, recJ or recQ mutation. These results demonstrate that the RecFOR pathway of recombination can contribute greatly to resistance to RM-mediated host killing, depending on the genetic background.

Bacterial toxin-antitoxin systems: more than selfish entities?

Bacterial toxin–antitoxin (TA) systems are diverse and widespread in the prokaryotic kingdom. They are composed of closely linked genes encoding a stable toxin that can harm the host cell and its cognate labile antitoxin, which protects the host from the toxin's deleterious effect. TA systems are thought to invade bacterial genomes through horizontal gene transfer. Some TA systems might behave as selfish elements and favour their own maintenance at the expense of their host. As a consequence, they may contribute to the maintenance of plasmids or genomic islands, such as super-integrons, by post-segregational killing of the cell that loses these genes and so suffers the stable toxin's destructive effect. The function of the chromosomally encoded TA systems is less clear and still open to debate. This Review discusses current hypotheses regarding the biological roles of these evolutionarily successful small operons. We consider the various selective forces that could drive the maintenance of TA systems in bacterial genomes.

From damaged genome to cell surface: transcriptome changes during bacterial cell death triggered by loss of a restriction-modification gene complex

Genetically programmed cell deaths play important roles in unicellular prokaryotes. In postsegregational killing, loss of a gene complex from a cell leads to its descendants' deaths. With type II restriction-modification gene complexes, such death is triggered by restriction endonuclease's attacks on under-methylated chromosomes. Here, we examined how the Escherichia coli transcriptome changes after loss of PaeR7I gene complex. At earlier time points, activation of SOS genes and sigma(E)-regulon was noticeable. With time, more SOS genes, stress-response genes (including sigma(S)-regulon, osmotic-, oxidative- and periplasmic-stress genes), biofilm-related genes, and many hitherto uncharacterized genes were induced, and genes for energy metabolism, motility and outer membrane biogenesis were repressed. As expected from the activation of sigma(E)-regulon, the death was accompanied by cell lysis and release of cellular proteins. Expression of several sigma(E)-regulon genes indeed led to cell lysis. We hypothesize that some signal was transduced, among multiple genes involved, from the damaged genome to the cell surface and led to its disintegration. These results are discussed in comparison with other forms of programmed deaths in bacteria and eukaryotes.

The Helicobacter pylori HpyAXII restriction-modification system limits exogenous DNA uptake by targeting GTAC sites but shows asymmetric conservation of the DNA methyltransferase and restriction endonuclease components

The naturally competent organism Helicobacter pylori encodes a large number of restriction-modification (R-M) systems that consist of a restriction endonuclease and a DNA methyltransferase. R-M systems are not only believed to limit DNA exchange among bacteria but may also have other cellular functions. We report a previously uncharacterized H. pylori type II R-M system, M.HpyAXII/R.HpyAXII. We show that this system targets GTAC sites, which are rare in the H. pylori chromosome but numerous in ribosomal RNA genes. As predicted, this type II R-M system showed attributes of a selfish element. Deletion of the methyltransferase M.HpyAXII is lethal when associated with an active endonuclease R.HpyAXII unless compensated by adaptive mutation or gene amplification. R.HpyAXII effectively restricted both unmethylated plasmid and chromosomal DNA during natural transformation and was predicted to belong to the novel 'half pipe' structural family of endonucleases. Analysis of a panel of clinical isolates revealed that R.HpyAXII was functional in a small number of H. pylori strains (18.9%, n = 37), whereas the activity of M.HpyAXII was highly conserved (92%, n = 50), suggesting that GTAC methylation confers a selective advantage to H. pylori. However, M.HpyAXII activity did not enhance H. pylori fitness during stomach colonization of a mouse infection model.

Cell death upon epigenetic genome methylation: a novel function of methyl-specific deoxyribonucleases

BACKGROUND

Alteration in epigenetic methylation can affect gene expression and other processes. In Prokaryota, DNA methyltransferase genes frequently move between genomes and present a potential threat. A methyl-specific deoxyribonuclease, McrBC, of Escherichia coli cuts invading methylated DNAs. Here we examined whether McrBC competes with genome methylation systems through host killing by chromosome cleavage.

RESULTS

McrBC inhibited the establishment of a plasmid carrying a PvuII methyltransferase gene but lacking its recognition sites, likely through the lethal cleavage of chromosomes that became methylated. Indeed, its phage-mediated transfer caused McrBC-dependent chromosome cleavage. Its induction led to cell death accompanied by chromosome methylation, cleavage and degradation. RecA/RecBCD functions affect chromosome processing and, together with the SOS response, reduce lethality. Our evolutionary/genomic analyses of McrBC homologs revealed: a wide distribution in Prokaryota; frequent distant horizontal transfer and linkage with mobility-related genes; and diversification in the DNA binding domain. In these features, McrBCs resemble type II restriction-modification systems, which behave as selfish mobile elements, maintaining their frequency by host killing. McrBCs are frequently found linked with a methyltransferase homolog, which suggests a functional association.

CONCLUSIONS

Our experiments indicate McrBC can respond to genome methylation systems by host killing. Combined with our evolutionary/genomic analyses, they support our hypothesis that McrBCs have evolved as mobile elements competing with specific genome methylation systems through host killing. To our knowledge, this represents the first report of a defense system against epigenetic systems through cell death.

Type II restriction endonuclease R.Hpy188I belongs to the GIY-YIG nuclease superfamily, but exhibits an unusual active site

BACKGROUND

Catalytic domains of Type II restriction endonucleases (REases) belong to a few unrelated three-dimensional folds. While the PD-(D/E)XK fold is most common among these enzymes, crystal structures have been also determined for single representatives of two other folds: PLD (R.BfiI) and half-pipe (R.PabI). Bioinformatics analyses supported by mutagenesis experiments suggested that some REases belong to the HNH fold (e.g. R.KpnI), and that a small group represented by R.Eco29kI belongs to the GIY-YIG fold. However, for a large fraction of REases with known sequences, the three-dimensional fold and the architecture of the active site remain unknown, mostly due to extreme sequence divergence that hampers detection of homology to enzymes with known folds.

RESULTS

R.Hpy188I is a Type II REase with unknown structure. PSI-BLAST searches of the non-redundant protein sequence database reveal only 1 homolog (R.HpyF17I, with nearly identical amino acid sequence and the same DNA sequence specificity). Standard application of state-of-the-art protein fold-recognition methods failed to predict the relationship of R.Hpy188I to proteins with known structure or to other protein families. In order to increase the amount of evolutionary information in the multiple sequence alignment, we have expanded our sequence database searches to include sequences from metagenomics projects. This search resulted in identification of 23 further members of R.Hpy188I family, both from metagenomics and the non-redundant database. Moreover, fold-recognition analysis of the extended R.Hpy188I family revealed its relationship to the GIY-YIG domain and allowed for computational modeling of the R.Hpy188I structure. Analysis of the R.Hpy188I model in the light of sequence conservation among its homologs revealed an unusual variant of the active site, in which the typical Tyr residue of the YIG half-motif had been substituted by a Lys residue. Moreover, some of its homologs have the otherwise invariant Arg residue in a non-homologous position in sequence that nonetheless allows for spatial conservation of the guanidino group potentially involved in phosphate binding.

CONCLUSION

The present study eliminates a significant "white spot" on the structural map of REases. It also provides important insight into sequence-structure-function relationships in the GIY-YIG nuclease superfamily. Our results reveal that in the case of proteins with no or few detectable homologs in the standard "non-redundant" database, it is useful to expand this database by adding the metagenomic sequences, which may provide evolutionary linkage to detect more remote homologs.

Genomics of bacteria and archaea: the emerging dynamic view of the prokaryotic world

The first bacterial genome was sequenced in 1995, and the first archaeal genome in 1996. Soon after these breakthroughs, an exponential rate of genome sequencing was established, with a doubling time of approximately 20 months for bacteria and approximately 34 months for archaea. Comparative analysis of the hundreds of sequenced bacterial and dozens of archaeal genomes leads to several generalizations on the principles of genome organization and evolution. A crucial finding that enables functional characterization of the sequenced genomes and evolutionary reconstruction is that the majority of archaeal and bacterial genes have conserved orthologs in other, often, distant organisms. However, comparative genomics also shows that horizontal gene transfer (HGT) is a dominant force of prokaryotic evolution, along with the loss of genetic material resulting in genome contraction. A crucial component of the prokaryotic world is the mobilome, the enormous collection of viruses, plasmids and other selfish elements, which are in constant exchange with more stable chromosomes and serve as HGT vehicles. Thus, the prokaryotic genome space is a tightly connected, although compartmentalized, network, a novel notion that undermines the ‘Tree of Life’ model of evolution and requires a new conceptual framework and tools for the study of prokaryotic evolution.

Maintenance forced by a restriction-modification system can be modulated by a region in its modification enzyme not essential for methyltransferase activity

Several type II restriction-modification gene complexes can force their maintenance on their host bacteria by killing cells that have lost them in a process called postsegregational killing or genetic addiction. It is likely to proceed by dilution of the modification enzyme molecule during rounds of cell division following the gene loss, which exposes unmethylated recognition sites on the newly replicated chromosomes to lethal attack by the remaining restriction enzyme molecules. This process is in apparent contrast to the process of the classical types of postsegregational killing systems, in which built-in metabolic instability of the antitoxin allows release of the toxin for lethal action after the gene loss. In the present study, we characterize a mutant form of the EcoRII gene complex that shows stronger capacity in such maintenance. This phenotype is conferred by an L80P amino acid substitution (T239C nucleotide substitution) mutation in the modification enzyme. This mutant enzyme showed decreased DNA methyltransferase activity at a higher temperature in vivo and in vitro than the nonmutated enzyme, although a deletion mutant lacking the N-terminal 83 amino acids did not lose activity at either of the temperatures tested. Under a condition of inhibited protein synthesis, the activity of the L80P mutant was completely lost at a high temperature. In parallel, the L80P mutant protein disappeared more rapidly than the wild-type protein. These results demonstrate that the capability of a restriction-modification system in forcing maintenance on its host can be modulated by a region of its antitoxin, the modification enzyme, as in the classical postsegregational killing systems.

Modelling the spatial dynamics of plasmid transfer and persistence

Bacterial plasmids are extra-chromosomal genetic elements that code for a wide variety of phenotypes in their bacterial hosts and are maintained in bacterial communities through both vertical and horizontal transfer. Current mathematical models of plasmid-bacteria dynamics, based almost exclusively on mass-action differential equations that describe these interactions in completely mixed environments, fail to adequately explain phenomena such as the long-term persistence of plasmids in natural and clinical bacterial communities. This failure is, at least in part, due to the absence of any spatial structure in these models, whereas most bacterial populations are spatially structured in microcolonies and biofilms. To help bridge the gap between theoretical predictions and observed patterns of plasmid spread and persistence, an individual-based lattice model (interacting particle system) that provides a predictive framework for understanding the dynamics of plasmid-bacteria interactions in spatially structured populations is presented here. To assess the accuracy and flexibility of the model, a series of experiments that monitored plasmid loss and horizontal transfer of the IncP-1beta plasmid pB10 : : rfp in Escherichia coli K12 and other bacterial populations grown on agar surfaces were performed. The model-based visual patterns of plasmid loss and spread, as well as quantitative predictions of the effects of different initial parental strain densities and incubation time on densities of transconjugants formed on a 2D grid, were in agreement with this and previously published empirical data. These results include features of spatially structured populations that are not predicted by mass-action differential equation models.