Restriction-modification systems as genomic parasites in competition for specific sequences

Dependence of post-segregational killing mediated by Type II restriction-modification systems on the lifetime of restriction endonuclease effective activity

Plasmid-borne Type II restriction-modification (RM) systems mediate post-segregational killing (PSK). PSK is thought to be caused by the dilution of restriction and modification enzymes during cell division, resulting in accumulation of unmethylated DNA recognition sites and their cleavage by restriction endonucleases. PSK is the likely reason for stabilization of plasmids carrying RM systems in the absence of selection for plasmid maintenance. In this study, we developed a CRISPR interference-based method to eliminate RM-carrying plasmids and study PSK-related phenomena with minimal perturbation to the Escherichia coli host. Plasmids carrying the EcoRV, Eco29kI, and EcoRI RM systems were highly stable, and their loss resulted in SOS response and PSK. In contrast, plasmids carrying the Esp1396I system were poorly stabilized; their loss led to a temporary cessation of growth, followed by full recovery. We demonstrate that this unusual behavior is due to a limited lifetime of the Esp1396I restriction endonuclease activity, which, upon Esp1396I plasmid loss, disappears approximately after two cycles of cell division, i.e., before unmethylated sites appear in significant numbers. Our results indicate that whenever PSK induced by a loss of RM systems, and, possibly, other toxin-antitoxin systems, is considered, the lifetimes of individual system components and the growth rate of host cells shall be taken in account. Mathematical modeling shows, that unlike the situation with classical toxin-antitoxin systems, RM system-mediated PSK is possible when the lifetimes of restriction endonuclease and methyltransferase activities are similar, as long as the toxic restriction endonuclease activity persists for more than two chromosome replication cycles.IMPORTANCEIt is widely accepted that many Type II restriction-modification (RM) systems mediate post-segregational killing (PSK) if plasmids that encode them are lost. In this study, we harnessed an inducible CRISPR-Cas system to remove RM plasmids from Escherichia coli cells to study PSK while minimally perturbing cell physiology. We demonstrate that PSK depends on restriction endonuclease activity lifetime and is not observed when it is less than two replication cycles. We present a mathematical model that explains experimental data and shows that unlike the case of toxin-antitoxin-mediated PSK, the loss of an RM system induced PSK even when the RM enzymes have identical lifetimes.

Restriction modification systems in archaea: A panoramic outlook

Restriction modification (RM) systems are one of the ubiquitous yet primitive defense responses employed by bacteria and archaea with the primary role of safeguarding themselves against invading bacteriophages. Protection of the host occurs by the cleavage of the invading foreign DNA via restriction endonucleases with concomitant methylation of host DNA with the aid of a methyltransferase counterpart. RM systems have been extensively studied in bacteria, however, in the case of archaea there are limited reports of RM enzymes that are investigated to date owing to their inhospitable growth demands. This review aims to broaden the knowledge about what is known about the diversity of RM systems in archaea and encapsulate the current knowledge on restriction and modification enzymes characterized in archaea so far and the role of RM systems in the milieu of archaeal biology.

A mathematician's guide to plasmids: an introduction to plasmid biology for modellers

Plasmids, extrachromosomal DNA molecules commonly found in bacterial and archaeal cells, play an important role in bacterial genetics and evolution. Our understanding of plasmid biology has been furthered greatly by the development of mathematical models, and there are many questions about plasmids that models would be useful in answering. In this review, we present an introductory, yet comprehensive, overview of the biology of plasmids suitable for modellers unfamiliar with plasmids who want to get up to speed and to begin working on plasmid-related models. In addition to reviewing the diversity of plasmids and the genes they carry, their key physiological functions, and interactions between plasmid and host, we also highlight selected plasmid topics that may be of particular interest to modellers and areas where there is a particular need for theoretical development. The world of plasmids holds a great variety of subjects that will interest mathematical biologists, and introducing new modellers to the subject will help to expand the existing body of plasmid theory.

OUP accepted manuscript

Epigenetic DNA methylation plays an important role in bacteria by influencing gene expression and allowing discrimination between self-DNA and intruders such as phages and plasmids. Restriction–modification (RM) systems use a methyltransferase (MTase) to modify a specific sequence motif, thus protecting host DNA from cleavage by a cognate restriction endonuclease (REase) while leaving invading DNA vulnerable. Other REases occur solitarily and cleave methylated DNA. REases and RM systems are frequently mobile, influencing horizontal gene transfer by altering the compatibility of the host for foreign DNA uptake. However, whether mobile defence systems affect pre-existing host defences remains obscure. Here, we reveal an epigenetic conflict between an RM system (PcaRCI) and a methylation-dependent REase (PcaRCII) in the plant pathogen Pectobacterium carotovorum RC5297. The PcaRCI RM system provides potent protection against unmethylated plasmids and phages, but its methylation motif is targeted by the methylation-dependent PcaRCII. This potentially lethal co-existence is enabled through epigenetic silencing of the PcaRCII-encoding gene via promoter methylation by the PcaRCI MTase. Comparative genome analyses suggest that the PcaRCII-encoding gene was already present and was silenced upon establishment of the PcaRCI system. These findings provide a striking example for selfishness of RM systems and intracellular competition between different defences.

Evolutionary Ecology and Interplay of Prokaryotic Innate and Adaptive Immune Systems

Like many organisms, bacteria and archaea have both innate and adaptive immune systems to defend against infection by viruses and other parasites. Innate immunity most commonly relies on the endonuclease-mediated cleavage of any incoming DNA that lacks a specific epigenetic modification, through a system known as restriction-modification. CRISPR-Cas-mediated adaptive immunity relies on the insertion of short DNA sequences from parasite genomes into CRISPR arrays on the host genome to provide sequence-specific protection. The discovery of each of these systems has revolutionised our ability to carry out genetic manipulations, and, as a consequence, the enzymes involved have been characterised in exquisite detail. In comparison, much less is known about the importance of these two arms of the defence for the ecology and evolution of prokaryotes and their parasites. Here, we review our current ecological and evolutionary understanding of these systems in isolation, and discuss the need to study how innate and adaptive immune responses are integrated when they coexist in the same cell.

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.

Toxin-Antidote Elements Across the Tree of Life

In life's constant battle for survival, it takes one to kill but two to conquer. Toxin-antitoxin or toxin-antidote (TA) elements are genetic dyads that cheat the laws of inheritance to guarantee their transmission to the next generation. This seemingly simple genetic arrangement—a toxin linked to its antidote—is capable of quickly spreading and persisting in natural populations. TA elements were first discovered in bacterial plasmids in the 1980s and have recently been characterized in fungi, plants, and animals, where they underlie genetic incompatibilities and sterility in crosses between wild isolates. In this review, we provide a unified view of TA elements in both prokaryotic and eukaryotic organisms and highlight their similarities and differences at the evolutionary, genetic, and molecular levels. Finally, we propose several scenarios that could explain the paradox of the evolutionary origin of TA elements and argue that these elements may be key evolutionary players and that the full scope of their roles is only beginning to be uncovered.

Intracellular Competitions Reveal Determinants of Plasmid Evolutionary Success

Plasmids are autonomously replicating genetic elements that are ubiquitous in all taxa and habitats where they constitute an integral part of microbial genomes. The stable inheritance of plasmids depends on their segregation during cell division and their long-term persistence in a host population is thought to largely depend on their impact on the host fitness. Nonetheless, many plasmids found in nature are lacking a clear trait that is advantageous to their host; the determinants of plasmid evolutionary success in the absence of plasmid benefit to the host remain understudied. Here we show that stable plasmid inheritance is an important determinant of plasmid evolutionary success. Borrowing terminology from evolutionary biology of cellular living forms, we hypothesize that Darwinian fitness is key for the plasmid evolutionary success. Performing intracellular plasmid competitions between non-mobile plasmids enables us to compare the evolutionary success of plasmid genotypes within the host, i.e., the plasmid fitness. Intracellular head-to-head competitions between stable and unstable variants of the same model plasmid revealed that the stable plasmid variant has a higher fitness in comparison to the unstable plasmid. Preemptive plasmid competitions reveal that plasmid fitness may depend on the order of plasmid arrival in the host. Competitions between plasmids characterized by similar stability of inheritance reveal plasmid fitness differences depending on the plasmid-encoded trait. Our results further reveal that competing plasmids can be maintained in coexistence following plasmid fusions that maintain unstable plasmid variants over time. Plasmids are not only useful accessory genetic elements to their host but they are also evolving and replicating entities, similarly to cellular living forms. There is a clear link between plasmid genetics and plasmid evolutionary success – hence plasmids are evolving entities whose fitness is quantifiable.

'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.

Escherichia coli ItaT is a type II toxin that inhibits translation by acetylating isoleucyl-tRNAIle

Prokaryotic toxin-antitoxin (TA) modules are highly abundant and are involved in stress response and drug tolerance. The most common type II TA modules consist of two interacting proteins. The type II toxins are diverse enzymes targeting various essential intracellular targets. The antitoxin binds to cognate toxin and inhibits its function. Recently, TA modules whose toxins are GNAT-family acetyltransferases were described. For two such systems, the target of acetylation was shown to be aminoacyl-tRNA: the TacT toxin targets aminoacylated elongator tRNAs, while AtaT targets the amino acid moiety of initiating tRNAMet. We show that the itaRT gene pair from Escherichia coli encodes a TA module with acetyltransferase toxin ItaT that specifically and exclusively acetylates Ile-tRNAIle thereby blocking translation and inhibiting cell growth. ItaT forms a tight complex with the ItaR antitoxin, which represses the transcription of itaRT operon. A comprehensive bioinformatics survey of GNAT acetyltransferases reveals that enzymes encoded by validated or putative TA modules are common and form a distinct branch of the GNAT family tree. We speculate that further functional analysis of such TA modules will result in identification of enzymes capable of specifically targeting many, perhaps all, aminoacyl tRNAs.

CRISPR-assisted multi-dimensional regulation for fine-tuning gene expression in Bacillus subtilis

Fine-tuning of gene expression is crucial for protein expression and pathway construction, but it still faces formidable challenges due to the hierarchical gene regulation at multiple levels in a context-dependent manner. In this study, we defined the optimal targeting windows for CRISPRa and CRISPRi of the dCas9-α/ω system, and demonstrated that this system could act as a single master regulator to simultaneously activate and repress the expression of different genes by designing position-specific gRNAs. The application scope of dCas9-ω was further expanded by a newly developed CRISPR-assisted Oligonucleotide Annealing based Promoter Shuffling (OAPS) strategy, which could generate a high proportion of functional promoter mutants and facilitate the construction of effective promoter libraries in microorganisms with low transformation efficiency. Combing OAPS and dCas9-ω, the influences of promoter-based transcription, molecular chaperone-assisted protein folding and protease-mediated degradation on the expression of amylase BLA in Bacillus subtilis were systematically evaluated, and a 260-fold enhancement of BLA production was obtained. The success of the OAPS strategy and dCas9-ω for BLA production in this study thus demonstrated that it could serve as a powerful tool kit to regulate the expression of multiple genes multi-directionally and multi-dimensionally in bacteria.

The Toxin-Antidote Model of Cytoplasmic Incompatibility: Genetics and Evolutionary Implications

Wolbachia bacteria inhabit the cells of about half of all arthropod species, an unparalleled success stemming in large part from selfish invasive strategies. Cytoplasmic incompatibility (CI), whereby the symbiont makes itself essential to embryo viability, is the most common of these and constitutes a promising weapon against vector-borne diseases. After decades of theoretical and experimental struggle, major recent advances have been made toward a molecular understanding of this phenomenon. As pieces of the puzzle come together, from yeast and Drosophila fly transgenesis to CI diversity patterns in natural mosquito populations, it becomes clearer than ever that the CI induction and rescue stem from a toxin-antidote (TA) system. Further, the tight association of the CI genes with prophages provides clues to the possible evolutionary origin of this phenomenon and the levels of selection at play.

Restriction glycosylases: involvement of endonuclease activities in the restriction process

All restriction enzymes examined are phosphodiesterases generating 3΄-OH and 5΄-P ends, but one restriction enzyme (restriction glycosylase) excises unmethylated bases from its recognition sequence. Whether its restriction activity involves endonucleolytic cleavage remains unclear. One report on this enzyme, R.PabI from a hyperthermophile, ascribed the breakage to high temperature while another showed its weak AP lyase activity generates atypical ends. Here, we addressed this issue in mesophiles. We purified R.PabI homologs from Campylobacter coli (R.CcoLI) and Helicobacter pylori (R.HpyAXII) and demonstrated their DNA cleavage, DNA glycosylase and AP lyase activities in vitro at 37°C. The AP lyase activity is more coupled with glycosylase activity in R.CcoLI than in R.PabI. R.CcoLI/R.PabI expression caused restriction of incoming bacteriophage/plasmid DNA and endogenous chromosomal DNA within Escherichia coli at 37°C. The R.PabI-mediated restriction was promoted by AP endonuclease action in vivo or in vitro. These results reveal the role of endonucleolytic DNA cleavage in restriction and yet point to diversity among the endonucleases. The cleaved ends are difficult to repair in vivo, which may indicate their biological significance. These results support generalization of the concept of restriction–modification system to the concept of self-recognizing epigenetic system, which combines any epigenetic labeling and any DNA damaging.

Plasmid stability analysis based on a new theoretical model employing stochastic simulations

Here, we present a simple theoretical model to study plasmid stability, based on one input parameter which is the copy number of plasmids present in a host cell. The Monte Carlo approach was used to analyze random fluctuations affecting plasmid replication and segregation leading to gradual reduction in the plasmid population within the host cell. This model was employed to investigate maintenance of pEC156 derivatives, a high-copy number ColE1-type Escherichia coli plasmid that carries an EcoVIII restriction-modification system. Plasmid stability was examined in selected Escherichia coli strains (MG1655, wild-type; MG1655 pcnB, and hyper-recombinogenic JC8679 sbcA). We have compared the experimental data concerning plasmid maintenance with the simulations and found that the theoretical stability patterns exhibited an excellent agreement with those empirically tested. In our simulations, we have investigated the influence of replication fails (α parameter) and uneven partition as a consequence of multimer resolution fails (δ parameter), and the post-segregation killing factor (β parameter). All of these factors act at the same time and affect plasmid inheritance at different levels. In case of pEC156-derivatives we concluded that multimerization is a major determinant of plasmid stability. Our data indicate that even small changes in the fidelity of segregation can have serious effects on plasmid stability. Use of the proposed mathematical model can provide a valuable description of plasmid maintenance, as well as enable prediction of the probability of the plasmid loss.

Exploring the parameters of post-segregational killing using heterologous expression of secreted toxin barnase and antitoxin barstar in an Escherichia coli case study

Post-segregational killing (PSK) is a phenotype determined by plasmids using a toxin and an antitoxin gene pair. Loss of the genes depletes the cell's reserve of antitoxin and allows the toxin to act upon the cell. PSK benefits mobile elements when it increases reproductive success relative to other mobile competitors. A side effect of PSK is that plasmids become refractory to displacement from the cell during growth as a monoculture. Most PSK systems use a cytoplasmic toxin, but the external toxins of bacteriocins also have a PSK-like effect. It may be that any toxin and antitoxin gene pair can demonstrate PSK when it is on a plasmid. The secreted ribonuclease barnase and its protein inhibitor barstar have features in common with PSK modules, though their native context is chromosomal. We hypothesized that their recruitment to a plasmid could produce an emergent PSK phenotype. Others had shown that secreted barnase could exert a lethal effect on susceptible bacteria similarly to bacteriocins. However, barnase toxicity did not occur under the conditions tested, suggesting that barnase is toxic to neighbouring cells only under very specific conditions. Bacteriocins are only produced under some conditions, and some conditionality on toxin function or release may be advantageous in general to PSKs with external toxins because it would prevent killing of potential plasmid-naive hosts. Too much conditionality, however, would limit how advantageous the gene pair was to mobile elements, making the genes unlikely to be recruited as a PSK system.

Regulation of genetic flux between bacteria by restriction-modification systems

Restriction-modification (R-M) systems are often regarded as bacteria's innate immune systems, protecting cells from infection by mobile genetic elements (MGEs). Their diversification has been recently associated with the emergence of particularly virulent lineages. However, we have previously found more R-M systems in genomes carrying more MGEs. Furthermore, it has been suggested that R-M systems might favor genetic transfer by producing recombinogenic double-stranded DNA ends. To test whether R-M systems favor or disfavor genetic exchanges, we analyzed their frequency with respect to the inferred events of homologous recombination and horizontal gene transfer within 79 bacterial species. Genetic exchanges were more frequent in bacteria with larger genomes and in those encoding more R-M systems. We created a recognition target motif predictor for Type II R-M systems that identifies genomes encoding systems with similar restriction sites. We found more genetic exchanges between these genomes, independently of their evolutionary distance. Our results reconcile previous studies by showing that R-M systems are more abundant in promiscuous species, wherein they establish preferential paths of genetic exchange within and between lineages with cognate R-M systems. Because the repertoire and/or specificity of R-M systems in bacterial lineages vary quickly, the preferential fluxes of genetic transfer within species are expected to constantly change, producing time-dependent networks of gene transfer.

High-resolution genetic analysis of the requirements for horizontal transmission of the ESBL plasmid from Escherichia coli O104:H4

Horizontal dissemination of the genes encoding extended spectrum beta-lactamases (ESBLs) via conjugative plasmids is facilitating the increasingly widespread resistance of pathogens to beta-lactam antibiotics. However, there is relatively little known about the regulatory factors and mechanisms that govern the spread of these plasmids. Here, we carried out a high-throughput, transposon insertion site sequencing analysis (TnSeq) to identify genes that enable the maintenance and transmission of pESBL, an R64 (IncI1)-related resistance plasmid that was isolated from Escherichia coli O104:H4 linked to a recent large outbreak of gastroenteritis. With a few exceptions, the majority of the genes identified as required for maintenance and transmission of pESBL matched those of their previously defined R64 counterparts. However, our analyses of the high-density transposon insertion library in pESBL also revealed two very short and linked regions that constitute a previously unrecognized regulatory system controlling spread of IncI1 plasmids. In addition, we investigated the function of the pESBL-encoded M.EcoGIX methyltransferase, which is also encoded by many other IncI1 and IncF plasmids. This enzyme proved to protect pESBL from restriction in new hosts, suggesting it aids in expanding the plasmid's host range. Collectively, our work illustrates the power of the TnSeq approach to enable rapid and comprehensive analyses of plasmid genes and sequences that facilitate the dissemination of determinants of antibiotic resistance.

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.

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.

Proteome variability among Helicobacter pylori isolates clustered according to genomic methylation

AIMS

To understand whether the variability found in the proteome of Helicobacter pylori relates to the genomic methylation, virulence and associated gastric disease.

METHODS AND RESULTS

We applied the Minimum-Common-Restriction-Modification (MCRM) algorithm to genomic methylation data of 30 Portuguese H. pylori strains, obtained by genome sensitivity to Type II restriction enzymes' digestion. All the generated dendrograms presented three clusters with no association with gastric disease. Comparative analysis of two-dimensional gel electrophoresis (2DE) maps obtained for total protein extracts of 10 of these strains, representative of the three main clusters, revealed that among 70 matched protein spots (in a universe of 300), 16 were differently abundant (P < 0·05) among clusters. Of these, 13 proteins appear to be related to the cagA genotype or gastric disease. The abundance of three protein species, DnaK, GlnA and HylB, appeared to be dictated by the methylation status of their gene promoter.

CONCLUSIONS

Variations in the proteome profile of strains with common geographic origin appear to be related to differences in cagA genotype or gastric disease, rather than to clusters organized according to strain genomic methylation.

SIGNIFICANCE AND IMPACT OF THE STUDY

The simultaneous study of the genomic methylation and proteome is important to correlate epigenetic modifications with gene expression and pathogen virulence.

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.

Screening for catalytically active Type II restriction endonucleases using segregation-induced methylation deficiency

Type II restriction endonucleases (REases) are one of the basic tools of recombinant DNA technology. They also serve as models for elucidation of mechanisms for both site-specific DNA recognition and cleavage by proteins. However, isolation of catalytically active mutants from their libraries is challenging due to the toxicity of REases in the absence of protecting methylation, and techniques explored so far had limited success. Here, we present an improved SOS induction-based approach for in vivo screening of active REases, which we used to isolate a set of active variants of the catalytic mutant, Cfr10I(E204Q). Detailed characterization of plasmids from 64 colonies screened from the library of ∼200,000 transformants revealed 29 variants of cfr10IR gene at the level of nucleotide sequence and 15 variants at the level of amino acid sequence, all of which were able to induce SOS response. Specific activity measurements of affinity-purified mutants revealed >200-fold variance among them, ranging from 100% (wild-type isolates) to 0.5% (S188C mutant), suggesting that the technique is equally suited for screening of mutants possessing high or low activity and confirming that it may be applied for identification of residues playing a role in catalysis.

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.

Cleavage of a model DNA replication fork by a methyl-specific endonuclease

Epigenetic DNA methylation is involved in many biological processes. An epigenetic status can be altered by gain or loss of a DNA methyltransferase gene or its activity. Repair of DNA damage can also remove DNA methylation. In response to such alterations, DNA endonucleases that sense DNA methylation can act and may cause cell death. Here, we explored the possibility that McrBC, a methylation-dependent DNase of Escherichia coli, cleaves DNA at a replication fork. First, we found that in vivo restriction by McrBC of bacteriophage carrying a foreign DNA methyltransferase gene is increased in the absence of homologous recombination. This suggests that some cleavage events are repaired by recombination and must take place during or after replication. Next, we demonstrated that the enzyme can cleave a model DNA replication fork in vitro. Cleavage of a fork required methylation on both arms and removed one, the other or both of the arms. Most cleavage events removed the methylated sites from the fork. This result suggests that acquisition of even rarely occurring modification patterns will be recognized and rejected efficiently by modification-dependent restriction systems that recognize two sites. This process might serve to maintain an epigenetic status along the genome through programmed cell death.

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.

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.

Plasmid addiction systems: perspectives and applications in biotechnology

Biotechnical production processes often operate with plasmid-based expression systems in well-established prokaryotic and eukaryotic hosts such as Escherichia coli or Saccharomyces cerevisiae, respectively. Genetically engineered organisms produce important chemicals, biopolymers, biofuels and high-value proteins like insulin. In those bioprocesses plasmids in recombinant hosts have an essential impact on productivity. Plasmid-free cells lead to losses in the entire product recovery and decrease the profitability of the whole process. Use of antibiotics in industrial fermentations is not an applicable option to maintain plasmid stability. Especially in pharmaceutical or GMP-based fermentation processes, deployed antibiotics must be inactivated and removed. Several plasmid addiction systems (PAS) were described in the literature. However, not every system has reached a full applicable state. This review compares most known addiction systems and is focusing on biotechnical applications.

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.

CRISPR-Cas: an adaptive immunity system in prokaryotes

Most of the archaea and numerous bacteria possess an elaborate system of adaptive immunity to mobile genetic elements known as the CRISPR (clustered regularly interspaced short palindromic repeats)-associated system (CRISPR-Cas), which consists of arrays of short repeats interspersed with unique DNA spacers and adjacent operons encompassing CRISPR-associated (cas) genes with predicted and, in some cases, experimentally validated nuclease, helicase, and polymerase activities. The system functions by integrating fragments of alien DNA between the repeats and employing their transcripts to degrade the DNA of the respective invading elements via an RNA interference-like mechanism. The CRISPR-Cas system is a case of apparent Lamarckian inheritance.

Geographic distribution of methyltransferases of Helicobacter pylori: evidence of human host population isolation and migration

Background

Helicobacter pylori colonizes the human stomach and is associated with gastritis, peptic ulcer, and gastric cancer. This ubiquitous association between H. pylori and humans is thought to be present since the origin of modern humans. The H. pylori genome encodes for an exceptional number of restriction and modifications (R-M) systems. To evaluate if R-M systems are an adequate tool to determine the geographic distribution of H. pylori strains, we typed 221 strains from Africa, America, Asia, and Europe, and evaluated the expression of different 29 methyltransferases.

Results

Independence tests and logistic regression models revealed that ten R-M systems correlate with geographical localization. The distribution pattern of these methyltransferases may have been originated by co-divergence of regional H. pylori after its human host migrated out of Africa. The expression of specific methyltransferases in the H. pylori population may also reflect the genetic and cultural background of its human host. Methyltransferases common to all strains, M. HhaI and M. NaeI, are likely conserved in H. pylori, and may have been present in the bacteria genome since the human diaspora out of Africa.

Conclusion

This study indicates that some methyltransferases are useful geomarkers, which allow discrimination of bacterial populations, and that can be added to our tools to investigate human migrations.

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.

Cleavage of a model DNA replication fork by a Type I restriction endonuclease

Cleavage of a DNA replication fork leads to fork restoration by recombination repair. In prokaryote cells carrying restriction-modification systems, fork passage reduces genome methylation by the modification enzyme and exposes the chromosome to attack by the restriction enzyme. Various observations have suggested a relationship between the fork and Type I restriction enzymes, which cleave DNA at a distance from a recognition sequence. Here, we demonstrate that a Type I restriction enzyme preparation cleaves a model replication fork at its branch. The enzyme probably tracks along the DNA from an unmethylated recognition site on the daughter DNA and cuts the fork upon encountering the branch point. Our finding suggests that these restriction-modification systems contribute to genome maintenance through cell death and indicates that DNA replication fork cleavage represents a critical point in genome maintenance to choose between the restoration pathway and the destruction pathway.

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.

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.

Bioinformatic and partial functional analysis of pEspA and pEspB, two plasmids from Exiguobacterium arabatum sp. nov. RFL1109

The complete nucleotide sequences of two plasmids from Exiguobacterium arabatum sp. nov. RFL1109, pEspA (4563bp) and pEspB (38,945bp), have been determined. Five ORFs were identified in the pEspA plasmid, and putative functions were assigned to two of them. Using deletion mapping approach, the Rep-independent replication region of pEspA, which functions in Bacillus subtilis, was localized within a 0.6kb DNA region. Analysis of the pEspB sequence revealed 42 ORFs. From these, function of two genes encoding enzymes of the Lsp1109I restriction-modification system was confirmed experimentally, while putative functions of another 18 ORFs were suggested based on comparative analysis. Three functional regions have been proposed for the pEspB plasmid: the putative conjugative transfer region, the region involved in plasmid replication and maintenance, and the region responsible for transposition of the IS21 family-like transposable elements.

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.

A new algorithm for cluster analysis of genomic methylation: the Helicobacter pylori case

MOTIVATION

The genomic methylation analysis is useful to type bacteria that have a high number of expressed type II methyltransferases. Methyltransferases are usually committed to Restriction and Modification (R-M) systems, in which the restriction endonuclease imposes high pressure on the expression of the cognate methyltransferase that hinder R-M system loss. Conventional cluster methods do not reflect this tendency. An algorithm was developed for dendrogram construction reflecting the propensity for conservation of R-M Type II systems.

RESULTS

The new algorithm was applied to 52 Helicobacter pylori strains from different geographical regions and compared with conventional clustering methods. The algorithm works by first grouping strains that share a common minimum set of R-M systems and gradually adds strains according to the number of the R-M systems acquired. Dendrograms revealed a cluster of African strains, which suggest that R-M systems are present in H.pylori genome since its human host migrates from Africa.

AVAILABILITY

The software files are available at http://www.ff.ul.pt/paginas/jvitor/Bioinformatics/MCRM_algorithm.zip.

A putative DNA adenine methyltransferase is involved in Yersinia pseudotuberculosis pathogenicity

Some adenine methyltransferases have been shown not only to protect specific DNA restriction sites from cleavage by a restriction endonuclease, but also to play a role in various bacterial processes and sometimes in bacterial virulence. This study focused on a type I restriction-modification system (designated yrmI) of Y. pseudotuberculosis. This system is composed of three adjacent genes which could potentially encode an N6-adenine DNA methylase (YamA), an enzyme involved in site-specific recognition (YrsA) and a restriction endonuclease (YreA). Screening of 85 isolates of Y. pestis and Y. pseudotuberculosis indicated that the yrmI system has been lost by Y. pestis and that yamA (but not yrsA or yreA) is present in all Y. pseudotuberculosis strains tested, suggesting that it may be important at some stages of the epidemiological cycle of this species. To further investigate the role of yamA in Y. pseudotuberculosis survival, multiplication or virulence, a DeltayamA mutant of Y. pseudotuberculosis IP32953 was constructed by allelic exchange with a kanamycin cassette. The fact that DeltayamA mutants were obtained indicated that this gene is not essential for Y. pseudotuberculosis viability. The IP32953DeltayamA mutant strain grew as well as the wild-type in a rich medium at both 28 degrees C and 37 degrees C. It also grew normally in a chemically defined medium at 28 degrees C, but exhibited a growth defect at 37 degrees C. In contrast to the Dam adenine methyltransferase, a mutation in yamA did not impair the functions of DNA repair or resistance to detergents. However, the DeltayamA mutant exhibited a virulence defect in a mouse model of intragastric infection. The in silico analysis indicated that the chromosomal region carrying the Y. pseudotuberculosis yrmI locus has been replaced in Y. pestis by a horizontally acquired region which potentially encodes another methyltransferase. YamA might thus be dispensable for Y. pestis growth and virulence because this species has acquired another gene fulfilling the same functions.

Negative regulation of the EcoRI restriction enzyme gene is associated with intragenic reverse promoters

Type II restriction-modification systems are expected to possess mechanisms for tight regulation of their expression to suppress the potential of lethal attack on their host bacteria when they establish and maintain themselves within them. Although the EcoRI restriction enzyme has been well characterized, regulation of its expression is still poorly understood. In this study, mutational analysis with lacZ gene fusion and primer extension assay identified a promoter for the transcription of the ecoRIR gene. Further analyses revealed that an intragenic region containing two overlapping reverse promoter-like elements acted as a negative regulator for ecoRIR gene expression. The activity of these putative reverse promoters was verified by transcriptional gene fusion, primer extension and in vitro transcription. Mutations in these reverse promoters resulted in increased gene expression in both translational and transcriptional gene fusions. An RNase protection assay revealed that the transcript level of the wild type relative to that of the reverse promoter mutant at the downstream regions was much lower than the level at the upstream regions. This suggests that these reverse promoter-like elements affect their downstream transcript level. The possible mechanisms of this kind of negative regulation, in addition to their possible biological roles, are discussed.

Genomic methylation: a tool for typing Helicobacter pylori isolates

The genome sequences of three Helicobacter pylori strains revealed an abundant number of putative restriction and modification (R-M) systems within a small genome (1.60 to 1.67 Mb). Each R-M system includes an endonuclease that cleaves a specific DNA sequence and a DNA methyltransferase that methylates either adenosine or cytosine within the same DNA sequence. These are believed to be a defense mechanism, protecting bacteria from foreign DNA. They have been classified as selfish genetic elements; in some instances it has been shown that they are not easily lost from their host cell. Possibly because of this phenomenon, the H. pylori genome is very rich in R-M systems, with considerable variation in potential recognition sequences. For this reason the protective aspect of the methyltransferase gene has been proposed as a tool for typing H. pylori isolates. We studied the expression of H. pylori methyltransferases by digesting the genomic DNAs of 50 strains with 31 restriction endonucleases. We conclude that methyltransferase diversity is sufficiently high to enable the use of the genomic methylation status as a typing tool. The stability of methyltransferase expression was assessed by comparing the methylation status of genomic DNAs from strains that were isolated either from the same patient at different times or from different stomach locations (antrum and corpus). We found a group of five methyltransferases common to all tested strains. These five may be characteristic of the genetic pool analyzed, and their biological role may be important in the host/bacterium interaction.

Chromosomal toxin-antitoxin loci can diminish large-scale genome reductions in the absence of selection

Superintegrons (SIs) are chromosomal genetic elements containing assemblies of genes, each flanked by a recombination sequence (attC site) targeted by the integron integrase. SIs may contain hundreds of attC sites and intrinsic instability is anticipated; yet SIs are remarkably stable. This implies that either selective pressure maintains the genes or mechanisms exist which favour their persistence in the absence of selection. Toxin/antitoxin (TA) systems encode a stable toxin and a specific, unstable antitoxin. Once activated, the continued synthesis of the unstable antitoxin is necessary for cell survival. A bioinformatic search of accessible microbial genomes for SIs and TA systems revealed that large SIs harboured TA gene cassettes while smaller SIs did not. We demonstrated the function of TA loci in different genomic contexts where large-scale deletions can occur; in SIs and in a 165 kb dispensable region of the Escherichia coli genome. When devoid of TA loci, large-scale genome loss was evident in both environments. The inclusion of two TA loci, relBE1 and parDE1, which we identified in the Vibrio vulnificus SI rendered these environments refractory to gene loss. Thus, chromosomal TA loci can stabilize massive SI arrays and limit the extensive gene loss that is a hallmark of reductive evolution.

Theoretical model of restriction endonuclease HpaI in complex with DNA, predicted by fold recognition and validated by site-directed mutagenesis

Type II restriction enzymes are commercially important deoxyribonucleases and very attractive targets for protein engineering of new specificities. At the same time they are a very challenging test bed for protein structure prediction methods. Typically, enzymes that recognize different sequences show little or no amino acid sequence similarity to each other and to other proteins. Based on crystallographic analyses that revealed the same PD-(D/E)XK fold for more than a dozen case studies, they were nevertheless considered to be related until the combination of bioinformatics and mutational analyses has demonstrated that some of these proteins belong to other, unrelated folds PLD, HNH, and GIY-YIG. As a part of a large-scale project aiming at identification of a three-dimensional fold for all type II REases with known sequences (currently approximately 1000 proteins), we carried out preliminary structure prediction and selected candidates for experimental validation. Here, we present the analysis of HpaI REase, an ORFan with no detectable homologs, for which we detected a structural template by protein fold recognition, constructed a model using the FRankenstein monster approach and identified a number of residues important for the DNA binding and catalysis. These predictions were confirmed by site-directed mutagenesis and in vitro analysis of the mutant proteins. The experimentally validated model of HpaI will serve as a low-resolution structural platform for evolutionary considerations in the subgroup of blunt-cutting REases with different specificities. The research protocol developed in the course of this work represents a streamlined version of the previously used techniques and can be used in a high-throughput fashion to build and validate models for other enzymes, especially ORFans that exhibit no sequence similarity to any other protein in the database.

Evolution of new protein topologies through multistep gene rearrangements

New protein folds have emerged throughout evolution, but it remains unclear how a protein fold can evolve while maintaining its function, particularly when fold changes require several sequential gene rearrangements. Here, we explored hypothetical evolutionary pathways linking different topological families of the DNA-methyltransferase superfamily. These pathways entail successive gene rearrangements through a series of intermediates, all of which should be sufficiently active to maintain the organism's fitness. By means of directed evolution, and starting from HaeIII methyltransferase (M.HaeIII), we selected all the required intermediates along these paths (a duplicated fused gene and duplicates partially truncated at their 5' or 3' coding regions) that maintained function in vivo. These intermediates led to new functional genes that resembled natural methyltransferases from three known classes or that belonged to a new class first seen in our evolution experiments and subsequently identified in natural genomes. Our findings show that new protein topologies can evolve gradually through multistep gene rearrangements and provide new insights regarding these processes.

Type III restriction is alleviated by bacteriophage (RecE) homologous recombination function but enhanced by bacterial (RecBCD) function

Previous works have demonstrated that DNA breaks generated by restriction enzymes stimulate, and are repaired by, homologous recombination with an intact, homologous DNA region through the function of lambdoid bacteriophages lambda and Rac. In the present work, we examined the effect of bacteriophage functions, expressed in bacterial cells, on restriction of an infecting tester phage in a simple plaque formation assay. The efficiency of plaque formation on an Escherichia coli host carrying EcoRI, a type II restriction system, is not increased by the presence of Rac prophage-presumably because, under the single-infection conditions of the plaque assay, a broken phage DNA cannot find a homologue with which to recombine. To our surprise, however, we found that the efficiency of plaque formation in the presence of a type III restriction system, EcoP1 or EcoP15, is increased by the bacteriophage-mediated homologous recombination functions recE and recT of Rac prophage. This type III restriction alleviation does not depend on lar on Rac, unlike type I restriction alleviation. On the other hand, bacterial RecBCD-homologous recombination function enhances type III restriction. These results led us to hypothesize that the action of type III restriction enzymes takes place on replicated or replicating DNA in vivo and leaves daughter DNAs with breaks at nonallelic sites, that bacteriophage-mediated homologous recombination reconstitutes an intact DNA from them, and that RecBCD exonuclease blocks this repair by degradation from the restriction breaks.

Stability of EcoRI restriction-modification enzymes in vivo differentiates the EcoRI restriction-modification system from other postsegregational cell killing systems

Certain type II restriction modification gene systems can kill host cells when these gene systems are eliminated from the host cells. Such ability to cause postsegregational killing of host cells is the feature of bacterial addiction modules, each of which consists of toxin and antitoxin genes. With these addiction modules, the differential stability of toxin and antitoxin molecules in cells plays an essential role in the execution of postsegregational killing. We here examined in vivo stability of the EcoRI restriction enzyme (toxin) and modification enzyme (antitoxin), the gene system of which has previously been shown to cause postsegregational host killing in Escherichia coli. Using two different methods, namely, quantitative Western blot analysis and pulse-chase immunoprecipitation analysis, we demonstrated that both the EcoRI restriction enzyme and modification enzyme are as stable as bulk cellular proteins and that there is no marked difference in their stability. The numbers of EcoRI restriction and modification enzyme molecules present in a host cell during the steady-state growth were estimated. We monitored changes in cellular levels of the EcoRI restriction and modification enzymes during the postsegregational killing. Results from these analyses together suggest that the EcoRI gene system does not rely on differential stability between the toxin and the antitoxin molecules for execution of postsegregational cell killing. Our results provide insights into the mechanism of postsegregational killing by restriction-modification systems, which seems to be distinct from mechanisms of postsegregational killing by other bacterial addiction modules.

Selection for plasmid post-segregational killing depends on multiple infection: evidence for the selection of more virulent parasites through parasite-level competition

Is the virulence of parasites an outcome of optimized infection? Virulence has often been considered an inevitable consequence of parasite reproduction when the cost incurred by the parasite in reducing the fitness of its current host is offset by increased infection of new hosts. More recent models have focused on how competition occurring between parasites during co-infection might effect selection of virulence. For example, if co-infection was common, parasites with higher intrinsic growth rates might be selected, even at the expense of being optimally adapted to infect new hosts. If growth rate is positively correlated with virulence, then competition would select increased virulence. We tested these models using a plasmid-encoded virulence determinant. The virulence determinant did not contribute to the plasmid's reproduction within or between hosts. Despite this, virulent plasmids were more successful than avirulent derivatives during selection in an environment allowing within-host competition. To explain these findings we propose and test a model in which virulent parasites are selected by reducing the reproduction of competitors.

Characterization of the Phd repressor-antitoxin boundary

The P1 plasmid addiction operon (a classic toxin-antitoxin system) encodes Phd, an unstable 73-amino-acid repressor-antitoxin protein, and Doc, a stable toxin. It was previously shown by deletion analysis that the N terminus of Phd was required for repressor activity and that the C terminus was required for antitoxin activity. Since only a quarter of the protein or less was required for both activities, it was hypothesized that Phd might have a modular organization. To further test the modular hypothesis, we constructed and characterized a set of 30 point mutations in the third and fourth quarters of Phd. Four mutations (PhdA36H, V37A, I38A, and F44A) had major defects in repressor activity. Five mutations (PhdD53A, D53R, E55A, F56A, and F60A) had major defects in antitoxin activity. As predicted by the modular hypothesis, point mutations affecting each activity belonged to disjoint, rather than overlapping, sets and were separated rather than interspersed within the linear sequence. A final deletion experiment demonstrated that the C-terminal 24 amino acid residues of Phd (preceded by a methionine) retained full antitoxin activity.

Toxin-antitoxin loci are highly abundant in free-living but lost from host-associated prokaryotes

Prokaryotic chromosomes code for toxin–antitoxin (TA) loci, often in multiple copies. In E.coli, experimental evidence indicates that TA loci are stress-response elements that help cells survive unfavorable growth conditions. The first gene in a TA operon codes for an antitoxin that combines with and neutralizes a regulatory ‘toxin’, encoded by the second gene. RelE and MazF toxins are regulators of translation that cleave mRNA and function, in interplay with tmRNA, in quality control of gene expression. Here, we present the results from an exhaustive search for TA loci in 126 completely sequenced prokaryotic genomes (16 archaea and 110 bacteria). We identified 671 TA loci belonging to the seven known TA gene families. Surprisingly, obligate intracellular organisms were devoid of TA loci, whereas free-living slowly growing prokaryotes had particularly many (38 in Mycobacterium tuberculosis and 43 in Nitrosomonas europaea). In many cases, TA loci were clustered and closely linked to mobile genetic elements. In the most extreme of these cases, all 13 TA loci of Vibrio cholerae were bona fide integron elements located in the V.cholerae mega-integron. These observations strongly suggest that TA loci are mobile cassettes that move frequently within and between chromosomes and also lend support to the hypothesis that TA loci function as stress-response elements beneficial to free-living prokaryotes.