Evolution of sequence recognition by restriction-modification enzymes: selective pressure for specificity decrease

Genetic conflicts and the origin of self/nonself-discrimination in the vertebrate immune system

Genetic conflicts shape the genomes of prokaryotic and eukaryotic organisms. Here, we argue that some of the key evolutionary novelties of adaptive immune systems of vertebrates are descendants of prokaryotic toxin-antitoxin (TA) systems. Cytidine deaminases and RAG recombinase have evolved from genotoxic enzymes to programmable editors of host genomes, supporting the astounding discriminatory capability of variable lymphocyte receptors of jawless vertebrates, as well as immunoglobulins and T cell receptors of jawed vertebrates. The evolutionarily recent lymphoid lineage is uniquely sensitive to mutations of the DNA maintenance methylase, which is an orphaned distant relative of prokaryotic restriction-modification systems. We discuss how the emergence of adaptive immunity gave rise to higher order genetic conflicts between genetic parasites and their vertebrate host.

The Influence of 5'R and 5'S cdA and cdG on the Activity of BsmAI and SspI Restriction Enzymes

Restriction endonucleases (REs) are intra-bacterial scissors that are considered tools in the fight against foreign genetic material. SspI and BsmAI, examined in this study, cleave dsDNA at their site of recognition or within a short distance of it. Both enzymes are representatives of type II REs, which have played an extremely important role in research on the genetics of organisms and molecular biology. Therefore, the study of agents affecting their activity has become highly important. Ionizing radiation may damage basic cellular mechanisms by inducing lesions in the genome, with 5',8-cyclo-2'-deoxypurines (cdPus) as a model example. Since cdPus may become components of clustered DNA lesions (CDLs), which are unfavorable for DNA repair pathways, their impact on other cellular mechanisms is worthy of attention. This study investigated the influence of cdPus on the elements of the bacterial restriction-modification system. In this study, it was shown that cdPus present in DNA affect the activity of REs. SspI was blocked by any cdPu lesion present at the enzyme's recognition site. When lesions were placed near the recognition sequence, the SspI was inhibited up to 46%. Moreover, (5'S)-5',8-cyclo-2'-deoxyadenosine (ScdA) present in the oligonucleotide sequence lowered BsmAI activity more than (5'R)-5',8-cyclo-2'-deoxyadenosine (RcdA). Interestingly, in the case of 5',8-cyclo-2'-deoxyguanosine (cdG), both 5'S and 5'R diastereomers inhibited BsmAI activity (up to 55% more than cdA). The inhibition was weaker when cdG was present at the recognition site rather than the cleavage site.

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.

Restriction enzymes use a 24 dimensional coding space to recognize 6 base long DNA sequences

Restriction enzymes recognize and bind to specific sequences on invading bacteriophage DNA. Like a key in a lock, these proteins require many contacts to specify the correct DNA sequence. Using information theory we develop an equation that defines the number of independent contacts, which is the dimensionality of the binding. We show that EcoRI, which binds to the sequence GAATTC, functions in 24 dimensions. Information theory represents messages as spheres in high dimensional spaces. Better sphere packing leads to better communications systems. The densest known packing of hyperspheres occurs on the Leech lattice in 24 dimensions. We suggest that the single protein EcoRI molecule employs a Leech lattice in its operation. Optimizing density of sphere packing explains why 6 base restriction enzymes are so common.

The thermostability and specificity of ancient proteins

Were ancient proteins systematically different than modern proteins? The answer to this question is profoundly important, shaping how we understand the origins of protein biochemical, biophysical, and functional properties. Ancestral sequence reconstruction (ASR), a phylogenetic approach to infer the sequences of ancestral proteins, may reveal such trends. We discuss two proposed trends: a transition from higher to lower thermostability and a tendency for proteins to acquire higher specificity over time. We review the evidence for elevated ancestral thermostability and discuss its possible origins in a changing environmental temperature and/or reconstruction bias. We also conclude that there is, as yet, insufficient data to support a trend from promiscuity to specificity. Finally, we propose future work to understand these proposed evolutionary trends.

Lifespan of restriction-modification systems critically affects avoidance of their recognition sites in host genomes

Background

Avoidance of palindromic recognition sites of Type II restriction-modification (R-M) systems was shown for many R-M systems in dozens of prokaryotic genomes. However the phenomenon has not been investigated systematically for all presently available genomes and annotated R-M systems. We have studied all known recognition sites in thousands of prokaryotic genomes and found factors that influence their avoidance.

Results

Only Type II R-M systems consisting of independently acting endonuclease and methyltransferase (called ‘orthodox’ here) cause avoidance of their sites, both palindromic and asymmetric, in corresponding prokaryotic genomes; the avoidance takes place for ~ 50 % of 1774 studied cases. It is known that prokaryotes can acquire and lose R-M systems. Thus it is possible to talk about the lifespan of an R-M system in a genome. We have shown that the recognition site avoidance correlates with the lifespan of R-M systems. The sites of orthodox R-M systems that are encoded in host genomes for a long time are avoided more often (up to 100 % in certain cohorts) than the sites of recently acquired ones. We also found cases of site avoidance in absence of the corresponding R-M systems in the genome. An analysis of closely related bacteria shows that such avoidance can be a trace of lost R-M systems. Sites of Type I, IIС/G, IIM, III, and IV R-M systems are not avoided in vast majority of cases.

Conclusions

The avoidance of orthodox Type II R-M system recognition sites in prokaryotic genomes is a widespread phenomenon. Presence of an R-M system without an underrepresentation of its site may indicate that the R-M system was acquired recently. At the same time, a significant underrepresentation of a site may be a sign of presence of the corresponding R-M system in this organism or in its ancestors for a long time. The drastic difference between site avoidance for orthodox Type II R-M systems and R-M systems of other types can be explained by a higher rate of specificity changes or a less self-toxicity of the latter.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2288-4) contains supplementary material, which is available to authorized users.

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.

Promiscuous restriction is a cellular defense strategy that confers fitness advantage to bacteria

Most bacterial genomes harbor restriction-modification systems, encoding a REase and its cognate MTase. On attack by a foreign DNA, the REase recognizes it as nonself and subjects it to restriction. Should REases be highly specific for targeting the invading foreign DNA? It is often considered to be the case. However, when bacteria harboring a promiscuous or high-fidelity variant of the REase were challenged with bacteriophages, fitness was maximal under conditions of catalytic promiscuity. We also delineate possible mechanisms by which the REase recognizes the chromosome as self at the noncanonical sites, thereby preventing lethal dsDNA breaks. This study provides a fundamental understanding of how bacteria exploit an existing defense system to gain fitness advantage during a host-parasite coevolutionary "arms race."

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.

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.

Activity and specificity of the bacterial PD-(D/E)XK homing endonuclease I-Ssp6803I

The restriction endonuclease fold [a three-layer alpha-beta sandwich containing variations of the PD-(D/E)XK nuclease motif] has been greatly diversified during evolution, facilitating its use for many biological functions. Here we characterize DNA binding and cleavage by the PD-(D/E)XK homing endonuclease I-Ssp6803I. Unlike most restriction endonucleases harboring the same core fold, the specificity profile of this enzyme extends over a long (17 bp) target site. The DNA binding and cleavage specificity profiles of this enzyme were independently determined and found to be highly correlated. However, the DNA target sequence contains several positions where binding and cleavage activities are not tightly coupled: individual DNA base-pair substitutions at those positions that significantly decrease cleavage activity have minor effects on binding affinity. These changes in the DNA target sequence appear to correspond to substitutions that uniquely increase the free energy change between the ground state and the transition state, rather than simply decreasing the overall DNA binding affinity. The specificity of the enzyme reflects constraints on its host gene and limitations imposed by the enzyme's quaternary structure and illustrate the highly diverse repertoire of DNA recognition specificities that can be adopted by the related folds surrounding the PD-(D/E)XK nuclease motif.

Structures of the rare-cutting restriction endonuclease NotI reveal a unique metal binding fold involved in DNA binding

The structure of the rare-cutting restriction endonuclease NotI, which recognizes the 8 bp target 5'-GCGGCCGC-3', has been solved with and without bound DNA. Because of its specificity (recognizing a site that occurs once per 65 kb), NotI is used to generate large genomic fragments and to map DNA methylation status. NotI contains a unique metal binding fold, found in a variety of putative endonucleases, occupied by an iron atom coordinated within a tetrahedral Cys4 motif. This domain positions nearby protein elements for DNA recognition, and serves a structural role. While recognition of the central six base pairs of the target is accomplished via a saturated hydrogen bond network typical of restriction enzymes, the most peripheral base pairs are engaged in a single direct contact in the major groove, reflecting reduced pressure to recognize those positions. NotI may represent an evolutionary intermediate between mobile endonucleases (which recognize longer target sites) and canonical restriction endonucleases.

Direct monitoring of allosteric recognition of type IIE restriction endonuclease EcoRII

EcoRII is a homodimer with two domains consisting of a DNA-binding N terminus and a catalytic C terminus and recognizes two specific sequences on DNA. It shows a relatively complicated cleavage reaction in bulk solution. After binding to either recognition site, EcoRII cleaves the other recognition site of the same DNA (cis-binding) strand and/or the recognition site of the other DNA (trans-binding) strand. Although it is difficult to separate these two reactions in bulk solution, we could simply obtain the binding and cleavage kinetics of only the cis-binding by following the frequency (mass) changes of a DNA-immobilized quartz-crystal microbalance (QCM) responding to the addition of EcoRII in aqueous solution. We obtained the maximum binding amounts (Deltam(max)), the dissociation constants (K(d)), the binding and dissociation rate constants (k(on) and k(off)), and the catalytic cleavage reaction rate constants (k(cat)) for wild-type EcoRII, the N-terminal-truncated form (EcoRII N-domain), and the mutant derivatives in its C-terminal domain (K263A and R330A). It was determined from the kinetic analyses that the N-domain, which covers the catalytic C-domain in the absence of DNA, preferentially binds to the one DNA recognition site while transforming EcoRII into an active form allosterically, and then the secondary C-domain binds to and cleaves the other recognition site of the DNA strand.

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.

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.

Common patterns in type II restriction enzyme binding sites

Restriction enzymes are among the best studied examples of DNA binding proteins. In order to find general patterns in DNA recognition sites, which may reflect important properties of protein–DNA interaction, we analyse the binding sites of all known type II restriction endonucleases. We find a significantly enhanced GC content and discuss three explanations for this phenomenon. Moreover, we study patterns of nucleotide order in recognition sites. Our analysis reveals a striking accumulation of adjacent purines (R) or pyrimidines (Y). We discuss three possible reasons: RR/YY dinucleotides are characterized by (i) stronger H-bond donor and acceptor clusters, (ii) specific geometrical properties and (iii) a low stacking energy. These features make RR/YY steps particularly accessible for specific protein–DNA interactions. Finally, we show that the recognition sites of type II restriction enzymes are underrepresented in host genomes and in phage genomes.

Multiplication of a restriction-modification gene complex

Previous works have suggested that some gene complexes encoding a restriction (R) enzyme and a cognate modification (M) enzyme may behave as selfish mobile genetic elements. RM gene complexes, which destroy 'non-self' elements marked by the absence of proper methylation, are often associated with mobile genetic elements and are involved in various genome rearrangements. Here, we found amplification of a restriction-modification gene complex. BamHI gene complex inserted into the Bacillus chromosome showed resistance to replacement by a homologous stretch of DNA. Some cells became transformed with the donor without losing BamHI. In most of these transformants, multiple copies of BamHI and the donor allele were arranged as tandem repeats. When a clone carrying one copy of each allele was propagated, extensive amplification of BamHI and the donor unit was observed in a manner dependent on restriction enzyme gene. This suggests that restriction cutting of the genome participates in the amplification. Visualization by fluorescent in situ hybridization revealed that the amplification occurred in single cells in a burst-like fashion that is reminiscent of induction of provirus replication. The multiplication ability in a bacterium with natural capacity for DNA release, uptake and transformation will be discussed in relation to spreading of RM gene -complexes.

A DNA methyltransferase can protect the genome from postdisturbance attack by a restriction-modification gene complex

In prokaryotic genomes, some DNA methyltransferases form a restriction-modification gene complex, but some others are present by themselves. Dcm gene product, one of these orphan methyltransferases found in Escherichia coli and related bacteria, methylates DNA to generate 5'-C(m)CWGG just as some of its eukaryotic homologues do. Vsr mismatch repair function of an adjacent gene prevents C-to-T mutagenesis enhanced by this methylation but promotes other types of mutation and likely has affected genome evolution. The reason for the existence of the dcm-vsr gene pair has been unclear. Earlier we found that several restriction-modification gene complexes behave selfishly in that their loss from a cell leads to cell killing through restriction attack on the genome. There is also increasing evidence for their potential mobility. EcoRII restriction-modification gene complex recognizes the same sequence as Dcm, and its methyltransferase is phylogenetically related to Dcm. In the present work, we found that stabilization of maintenance of a plasmid by linkage of EcoRII gene complex, likely through postsegregational cell killing, is diminished by dcm function. Disturbance of EcoRII restriction-modification gene complex led to extensive chromosome degradation and severe loss of cell viability. This cell killing was partially suppressed by chromosomal dcm and completely abolished by dcm expressed from a plasmid. Dcm, therefore, can play the role of a "molecular vaccine" by defending the genome against parasitism by a restriction-modification gene complex.

Sequence permutations in the molecular evolution of DNA methyltransferases

BACKGROUND

DNA methyltransferases (MTases), unlike MTases acting on other substrates, exhibit sequence permutation. Based on the sequential order of the cofactor-binding subdomain, the catalytic subdomain, and the target recognition domain (TRD), several classes of permutants have been proposed. The majority of known DNA MTases fall into the alpha, beta, and gamma classes. There is only one member of the zeta class known and no members of the delta and epsilon classes have been identified to date. Two mechanisms of permutation have been proposed: one involving gene duplication and in-frame fusion, and the other involving inter- and intragenic shuffling of gene segments.

RESULTS

Two novel cases of sequence permutation in DNA MTases implicated in restriction-modification systems have been identified, which suggest that members of the delta and zeta classes (M.MwoI and M.TvoORF1413P, respectively) evolved from beta-class MTases. This is the first identification of the delta-class MTase and the second known zeta-class MTase (the first zeta-class member among DNA:m4C and m6A-MTases).

CONCLUSIONS

Fragmentation of a DNA MTase gene may result from attack of nucleases, for instance when the RM system invades a new cell. Its reassembly into a functional form, the order of motifs notwithstanding, may be strongly selected for, if the cognate ENase gene remains active and poses a threat to the host's chromosome. The "cut-and-paste" mechanism is proposed for beta-delta permutation, which is non-circular and involves relocation of one segment of a gene. The circular beta-zeta permutation may be explained both by gene duplication or shuffling of gene fragments. These two mechanisms are not mutually exclusive and probably both played a role in the evolution of permuted DNA MTases.

Behavior of restriction-modification systems as selfish mobile elements and their impact on genome evolution

Restriction-modification (RM) systems are composed of genes that encode a restriction enzyme and a modification methylase. RM systems sometimes behave as discrete units of life, like viruses and transposons. RM complexes attack invading DNA that has not been properly modified and thus may serve as a tool of defense for bacterial cells. However, any threat to their maintenance, such as a challenge by a competing genetic element (an incompatible plasmid or an allelic homologous stretch of DNA, for example) can lead to cell death through restriction breakage in the genome. This post-segregational or post-disturbance cell killing may provide the RM complexes (and any DNA linked with them) with a competitive advantage. There is evidence that they have undergone extensive horizontal transfer between genomes, as inferred from their sequence homology, codon usage bias and GC content difference. They are often linked with mobile genetic elements such as plasmids, viruses, transposons and integrons. The comparison of closely related bacterial genomes also suggests that, at times, RM genes themselves behave as mobile elements and cause genome rearrangements. Indeed some bacterial genomes that survived post-disturbance attack by an RM gene complex in the laboratory have experienced genome rearrangements. The avoidance of some restriction sites by bacterial genomes may result from selection by past restriction attacks. Both bacteriophages and bacteria also appear to use homologous recombination to cope with the selfish behavior of RM systems. RM systems compete with each other in several ways. One is competition for recognition sequences in post-segregational killing. Another is super-infection exclusion, that is, the killing of the cell carrying an RM system when it is infected with another RM system of the same regulatory specificity but of a different sequence specificity. The capacity of RM systems to act as selfish, mobile genetic elements may underlie the structure and function of RM enzymes.

Insertion with long target duplication: a mechanism for gene mobility suggested from comparison of two related bacterial genomes

The complete genome sequences of two closely related organisms--two Helicobacter pylori strains--have recently become available. Comparison of these genomes at single base pair level has suggested the presence of a mechanism for bacterial gene mobility--insertion with long target duplications. This mechanism is formally similar to classical transposon insertion, but the duplication is much longer, often in the range of 100bp. Restriction and/or modification enzyme genes are often within or adjacent to the insertion. A similar process may have mediated insertion of the cag(+) pathogenicity island in H. pylori. A similar structure was identified in comparisons between Neisseria meningitidis and Neisseria gonorrhoeae genomes. We hypothesize that this mechanism, as well as two other types of polymorphism linked with restriction-modification genes (insertion accompanied by target deletion and a tripartite structure composed of substitution/inversion/deletion), have resulted from attack by restriction enzymes on the chromosome.