Characterization of pEC156, a ColE1-type plasmid from Escherichia coli E1585-68 that carries genes of the EcoVIII restriction-modification system

Regulator-dependent temporal dynamics of a restriction-modification system's gene expression upon entering new host cells: single-cell and population studies

Restriction-modification (R-M) systems represent a first line of defense against invasive DNAs, such as bacteriophage DNAs, and are widespread among bacteria and archaea. By acquiring a Type II R-M system via horizontal gene transfer, the new hosts generally become more resistant to phage infection, through the action of a restriction endonuclease (REase), which cleaves DNA at or near specific sequences. A modification methyltransferase (MTase) serves to protect the host genome against its cognate REase activity. The production of R-M system components upon entering a new host cell must be finely tuned to confer protective methylation before the REase acts, to avoid host genome damage. Some type II R-M systems rely on a third component, the controller (C) protein, which is a transcription factor that regulates the production of REase and/or MTase. Previous studies have suggested C protein effects on the dynamics of expression of an R-M system during its establishment in a new host cell. Here, we directly examine these effects. By fluorescently labelling REase and MTase, we demonstrate that lack of a C protein reduces the delay of REase production, to the point of being simultaneous with, or even preceding, production of the MTase. Single molecule tracking suggests that a REase and a MTase employ different strategies for their target search within host cells, with the MTase spending much more time diffusing in proximity to the nucleoid than does the REase. This difference may partially ameliorate the toxic effects of premature REase expression.

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.

Natural tuning of restriction endonuclease synthesis by cluster of rare arginine codons

Restriction–modification (R-M) systems are highly widespread among bacteria and archaea, and they appear to play a pivotal role in modulating horizontal gene transfer, as well as in protecting the host organism against viruses and other invasive DNA particles. Type II R-M systems specify two independent enzymes: a restriction endonuclease (REase) and protective DNA methyltransferase (MTase). If the cell is to survive, the counteracting activities as toxin and antitoxin, must be finely balanced in vivo. The molecular basis of this regulatory process remains unclear and current searches for regulatory elements in R-M modules are focused mainly at the transcription step. In this report, we show new aspects of REase control that are linked to translation. We used the EcoVIII R-M system as a model. Both, the REase and MTase genes for this R-M system contain an unusually high number of rare arginine codons (AGA and AGG) when compared to the rest of the E. coli K-12 genome. Clusters of these codons near the N-terminus of the REase greatly affect the translational efficiency. Changing these to higher frequency codons for E. coli (CGC) improves the REase synthesis, making the R-M system more potent to defend its host against bacteriophages. However, this improved efficiency in synthesis reduces host fitness due to increased autorestriction. We hypothesize that expression of the endonuclease gene can be modulated depending on the host genetic context and we propose a novel post-transcriptional mode of R–M system regulation that alleviates the potential lethal action of the restriction enzyme.

Isospecific adenine DNA methyltransferases show distinct preferences towards DNA substrates

Here, we report results on systematic analysis of DNA substrate preferences of three N6-adenine β-class DNA methyltransferases that are part of the type II restriction-modification systems. The studied enzymes were: M.EcoVIII, M.HindIII and M.LlaCI, which although found in phylogenetically distant bacteria (γ-proteobacteria and low-GC Gram-positive bacteria), recognize the same palindromic specific sequence 5′-AAGCTT-3′ and catalyze formation of N6-methyladenine at the first A-residue. As expected overall the enzymes share the most analyzed features, but they show also some distinct differences in substrate recognition. Therefore DNA methylation reactions were carried out not only under standard, but also under relaxed conditions using DMSO or glycerol. We found that all of these enzymes preferred DNA containing a hemimethylated target site, but differ in modification of ssDNA, especially more pronounced for M.EcoVIII under relaxed conditions. In these conditions they also have shown varied preferences toward secondary sites, which differ by one nucleotide from specific sequence. They preferred sequences with substitutions at the 1^st (A¹ → G/C) and at the 2^nd position (A² → C), while sites with substitutions at the 3^rd position (G³ → A/C) were modified less efficiently. Kinetic parameters of the methylation reaction carried out by M.EcoVIII were determined. Methylation efficiency (k_cat/K_m) of secondary sites was 4.5–10 times lower when compared to the unmethylated specific sequences, whilst efficiency observed for the hemimethylated substrate was almost 4.5 times greater. We also observed a distinct effect of analyzed enzymes on unspecific interaction with DNA phosphate backbone. We concluded that for all three enzymes the most critical is the phosphodiester bond between G³-C⁴ nucleotides at the center of the target site.

Isolation and sequence analysis of pCS36-4CPA, a small plasmid from Citrobacter sp. 36-4CPA

A small plasmid designated pCS36-4CPA with a size of 5217 base pairs and G-C content of 50.74% was isolated from Citrobacter sp. 36-4CPA. The origin of replication (ori) of the plasmid was identified as a region of about 800 bp in length with an identity of 67.1% to the ColE1 plasmid at the nucleotide level. The replication region contained typical elements of ColE1-like plasmids: RNA I and RNA II with their corresponding -10 and -35 boxes, a single-strand initiation site (ssi), and a lagging-strand termination site (terH). As seen in other ColE1-like plasmids, pCS36-4CPA carried mobilisation machinery that include mobABCD genes but it did not possess the rom gene. Analysis of the multimer resolution site (mrs) was performed and XerC and XerD binding sites were identified. Also, the 70-nt transcript Rcd of pCS36-4CPA was predicted and similarity of the transcript's secondary structure with those of the ColE1-family was shown. The cargo module of pCS36-4CPA contained three open reading frames (ORFs). Two of them (ORF5 and ORF6) showed no significant homology to any known gene sequences but contained putative THAP DNA-binding (DBD) and type II restriction endonuclease EcoO109I domains. The seventh open reading frame (ORF7) encodes YhdJ-like DNA modification methylase. The region highly homologous to pCS36-4CPA was found in the Salmonella phage SE2 genome.

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.

Plasmid pEC156, a Naturally Occurring Escherichia coli Genetic Element That Carries Genes of the EcoVIII Restriction-Modification System, Is Mobilizable among Enterobacteria

Type II restriction-modification systems are ubiquitous in prokaryotes. Some of them are present in naturally occurring plasmids, which may facilitate the spread of these systems in bacterial populations by horizontal gene transfer. However, little is known about the routes of their dissemination. As a model to study this, we have chosen an Escherichia coli natural plasmid pEC156 that carries the EcoVIII restriction modification system. The presence of this system as well as the cis-acting cer site involved in resolution of plasmid multimers determines the stable maintenance of pEC156 not only in Escherichia coli but also in other enterobacteria. We have shown that due to the presence of oriT-type F and oriT-type R64 loci it is possible to mobilize pEC156 by conjugative plasmids (F and R64, respectively). The highest mobilization frequency was observed when pEC156-derivatives were transferred between Escherichia coli strains, Enterobacter cloacae and Citrobacter freundii representing coliform bacteria. We found that a pEC156-derivative with a functional EcoVIII restriction-modification system was mobilized in enterobacteria at a frequency lower than a plasmid lacking this system. In addition, we found that bacteria that possess the EcoVIII restriction-modification system can efficiently release plasmid content to the environment. We have shown that E. coli cells can be naturally transformed with pEC156-derivatives, however, with low efficiency. The transformation protocol employed neither involved chemical agents (e.g. CaCl₂) nor temperature shift which could induce plasmid DNA uptake.

Natural C-independent expression of restriction endonuclease in a C protein-associated restriction-modification system

Restriction-modification (R-M) systems are highly prevalent among bacteria and archaea, and appear to play crucial roles in modulating horizontal gene transfer and protection against phage. There is much to learn about these diverse enzymes systems, especially their regulation. Type II R-M systems specify two independent enzymes: a restriction endonuclease (REase) and protective DNA methyltransferase (MTase). Their activities need to be finely balanced in vivo Some R-M systems rely on specialized transcription factors called C (controller) proteins. These proteins play a vital role in the temporal regulation of R-M gene expression, and function to indirectly modulate the horizontal transfer of their genes across the species. We report novel regulation of a C-responsive R-M system that involves a C protein of a poorly-studied structural class - C.Csp231I. Here, the C and REase genes share a bicistronic transcript, and some of the transcriptional auto-control features seen in other C-regulated R-M systems are conserved. However, separate tandem promoters drive most transcription of the REase gene, a distinctive property not seen in other tested C-linked R-M systems. Further, C protein only partially controls REase expression, yet plays a role in system stability and propagation. Consequently, high REase activity was observed after deletion of the entire C gene, and cells bearing the ΔC R-M system were outcompeted in mixed culture assays by those with the WT R-M system. Overall, our data reveal unexpected regulatory variation among R-M systems.

Relaxed specificity of prokaryotic DNA methyltransferases results in DNA site-specific modification of RNA/DNA heteroduplexes

RNA/DNA hybrid duplexes regularly occur in nature, for example in transcriptional R loops. Their susceptibility to modification by DNA-specific or RNA-specific enzymes is, thus, a biologically relevant question, which, in addition, has possible biotechnological implications. In this study, we investigated the activity of four isospecific DNA methyltransferases (M.EcoVIII, M.LlaCI, M.HindIII, M.BstZ1II) toward an RNA/DNA duplex carrying one 5'-AAGCUU-3'/3'-TTCGAA-5' target sequence. The analyzed enzymes belong to the β-group of adenine N6-methyltransferases and recognize the palindromic DNA sequence 5'-AAGCTT-3'/3'-TTCGAA-5'. Under standard conditions, none of these isospecific enzymes could detectibly methylate the RNA/DNA duplex. However, the addition of agents that generally relax specificity, such as dimethyl sulfoxide (DMSO) and glycerol, resulted in substantial methylation of the RNA/DNA duplex by M.EcoVIII and M.LlaCI. Only the DNA strand of the RNA/DNA duplex was methylated. The same was not observed for M.HindIII or M.BstZ1II. This is, to our knowledge, the first report that demonstrates such activity by prokaryotic DNA methyltransferases. Possible applications of these findings in a laboratory practice are also discussed.

Genetic analysis of maintenance of pEC156, a naturally occurring Escherichia coli plasmid that carries genes of the EcoVIII restriction-modification system

In the present study the role of the mechanisms responsible for maintenance of a natural plasmid pEC156, that carries genes of the EcoVIII restriction-modification system was investigated. Analysis of this plasmid's genetic content revealed the presence of genetic determinants suggesting two such mechanisms. The first of them relies on site specific recombination utilizing the Xer/cer molecular machinery, while the second involves a restriction-modification system as an addiction module. Our analysis indicated that three factors affect the maintenance of pEC156: (i) a cis-acting cer site involved in resolution of plasmid multimers, (ii) a gene coding for EcoVIII endonuclease, and (iii) plasmid copy number control. The lowest stability was observed with pEC156 derivatives deprived of the cer site. Decreased stability of pEC156 derivatives was also observed in E.coli strains deficient in genes coding for proteins involved in plasmid multimer resolution (XerC, XerD, ArgR and PepA). A similar effect, but to a much lesser extent was observed for the pEC156 derivative without a functional gene coding for EcoVIII endonuclease. Our results indicate that the presence of the cer site is more important for pEC156 stable maintenance than the presence of a functional gene coding for EcoVIII endonuclease. In our work we also tested maintenance of pEC156 possessing a ColE1-type replicon in bacteria belonging to Enterobacteriaceae family. We have found that pEC156 was most stably maintained in Enterobacter cloacae and Klebsiella oxytoca representing coli-type enterobacteria. We have found that in all enterobacteria tested pEC156 derivatives deficient in the cer site were significantly less stably maintained than cer(+) variants.

Plasmid Diversity and Adaptation Analyzed by Massive Sequencing of Escherichia coli Plasmids

Whole-genome sequencing is revolutionizing the analysis of bacterial genomes. It leads to a massive increase in the amount of available data to be analyzed. Bacterial genomes are usually composed of one main chromosome and a number of accessory chromosomes, called plasmids. A recently developed methodology called PLACNET (for pla smid c onstellation net works) allows the reconstruction of the plasmids of a given genome. Thus, it opens an avenue for plasmidome analysis on a global scale. This work reviews our knowledge of the genetic determinants for plasmid propagation (conjugation and related functions), their diversity, and their prevalence in the variety of plasmids found by whole-genome sequencing. It focuses on the results obtained from a collection of 255 Escherichia coli plasmids reconstructed by PLACNET. The plasmids found in E. coli represent a nonaleatory subset of the plasmids found in proteobacteria. Potential reasons for the prevalence of some specific plasmid groups will be discussed and, more importantly, additional questions will be posed.

Molecular characterization of plasmid pMbo4.6 of Moraxella bovis ATCC 10900

We report the characterization of a small cryptic plasmid unlike any previously described from Moraxella bovis ATCC 10900, a Gram-negative bacterium belonging to the family Moraxellaceae. The complete nucleotide sequence of the plasmid pMbo4.6 was determined. The plasmid was analyzed and found to be 4658 in size with a G+C content of 38.6 mol %. Computer analysis of the sequence data revealed four major open reading frames encoding putative proteins of 10.1 (ORF1), 64.2 (ORF2), 45.7 (ORF3), and 12.1 kDa (ORF4). ORF1 and ORF2 encode proteins that show a high level of amino acid sequence similarity (44 %) with some mobilization proteins. ORF3 encodes a protein showing a relatively high amino acid sequence similarity (about 40 %) with several plasmid replication initiator proteins. Upstream of ORF3, a 320-bp intergenic region, constituting the putative origin of replication that contained an AT-rich region followed by four direct repeats, was identified. This set of repeated sequences resembles iteron structures and plays an important role in the control of plasmid replication by providing a target site for the initiation of transcription and replication factors (IHF and RepA). Several palindromic sequences, inverted repeats, and hairpin-loop structures, which might confer regulatory effects on the replication of the plasmid, were also noted. ORF4 encodes an uncharacterized protein, conserved in bacteria, belonging to the DUF497 family. Sequence analysis and structural features indicate that pMbo4.6 replicates by a theta mechanism.

A rapid and efficient method for cloning genes of type II restriction-modification systems by use of a killer plasmid

We present a method for cloning restriction-modification (R-M) systems that is based on the use of a lethal plasmid (pKILLER). The plasmid carries a functional gene for a restriction endonuclease having the same DNA specificity as the R-M system of interest. The first step is the standard preparation of a representative, plasmid-borne genomic library. Then this library is transformed with the killer plasmid. The only surviving bacteria are those which carry the gene specifying a protective DNA methyltransferase. Conceptually, this in vivo selection approach resembles earlier methods in which a plasmid library was selected in vitro by digestion with a suitable restriction endonuclease, but it is much more efficient than those methods. The new method was successfully used to clone two R-M systems, BstZ1II from Bacillus stearothermophilus 14P and Csp231I from Citrobacter sp. strain RFL231, both isospecific to the prototype HindIII R-M system.

Sequence analysis of the plasmid pColG from the Escherichia coli strain CA46

The complete 4715 nucleotide sequence of the plasmid pColG from the Escherichia coli strain CA46, which was originally assumed to code for colicin G activity, has been determined. Based on the nucleotide sequence homology of the 1828bp replication region, with an average G+C content of 48%, pColG was classified as a ColE1-like plasmid. Computer assisted analysis of the remaining 2887bp nucleotide sequence with an average G+C content of 34% revealed three putative OFRs. To find out whether one or all of the three ORFs code for a possible bacteriocin, a DNA fragment encompassing these ORFs has been cloned and the recombinant colonies tested for bacteriocin production. None of the colonies had an inhibitory activity against E. coli strains DH5, HB101 and MC4100. The assumption that the plasmid pColG from the E. coli strain CA46 codes for a bacteriocin thus could not be confirmed.

Salmonella enterica subsp. enterica serovar Enteritidis harbours ColE1, ColE2, and rolling-circle-like replicating plasmids

Using DNA hybridization, at least three distinct groups of low molecular mass plasmids were identified in Salmonella enterica subsp. enterica serovar Enteritidis. After sequencing representative plasmids from each group, we concluded that they belonged to ColE1, ColE2, and rolling-circle-like replicating plasmids. Plasmid pK (4245 bp) is a representative of widely distributed ColE1 plasmids. Plasmid pP (4301 bp) is homologous to ColE2 plasmids and was present predominantly in single-stranded DNA form. The smallest plasmids pJ (2096 bp) and pB (1983 bp) were classified as rolling-circle-like replicating plasmids. Both encoded only a single protein essential for their own replication, and they must have existed in an unusual molecular structure, as (i) they were capable of hybridization without denaturation, (ii) their DNA could be linearized with S1 nuclease, and (iii) even after such treatment, the ability to hybridize without denaturation did not disappear.

Characterization of the LlaCI methyltransferase from Lactococcus lactis subsp. cremoris W15 provides new insights into the biology of type II restriction-modification systems

The gene encoding the LlaCI methyltransferase (M.LlaCI) from Lactococcus lactis subsp. cremoris W15 was overexpressed in Escherichia coli. The enzyme was purified to apparent homogeneity using three consecutive steps of chromatography on phosphocellulose, blue-agarose and Superose 12HR, yielding a protein of M(r) 31 300+/-1000 under denaturing conditions. The exact position of the start codon AUG was determined by protein microsequencing. This enzyme recognizes the specific palindromic sequence 5'-AAGCTT-3'. Purified M.LlaCI was characterized. Unlike many other methyltransferases, M.LlaCI exists in solution predominantly as a dimer. It modifies the first adenine residue at the 5' end of the specific sequence to N(6)-methyladenine and thus is functionally identical to the corresponding methyltransferases of the HindIII (Haemophilus influenzae Rd) and EcoVIII (Escherichia coli E1585-68) restriction-modification systems. This is reflected in the identity of M.LlaCI with M.HindIII and M.EcoVIII noted at the amino acid sequence level (50 % and 62 %, respectively) and in the presence of nine sequence motifs conserved among N(6)-adenine beta-class methyltransferases. However, polyclonal antibodies raised against M.EcoVIII cross-reacted with M.LlaCI but not with M.HindIII. Restriction endonucleases require Mg(2+) for phosphodiester bond cleavage. Mg(2+) was shown to be a strong inhibitor of the M.LlaCI enzyme and its isospecific homologues. This observation suggests that sensitivity of the M.LlaCI to Mg(2+) may strengthen the restriction activity of the cognate endonuclease in the bacterial cell. Other biological implications of this finding are also discussed.

Genetic organization and molecular analysis of the EcoVIII restriction-modification system of Escherichia coli E1585-68 and its comparison with isospecific homologs

The EcoVIII restriction-modification (R-M) system is carried by the Escherichia coli E1585-68 natural plasmid pEC156 (4,312 bp). The two genes were cloned and characterized. The G+C content of the EcoVIII R-M system is 36.1%, which is significantly lower than the average G+C content of either plasmid pEC156 (43.6%) or E. coli genomic DNA (50.8%). The difference suggests that there is a possibility that the EcoVIII R-M system was recently acquired by the genome. The 921-bp EcoVIII endonuclease (R. EcoVIII) gene (ecoVIIIR) encodes a 307-amino-acid protein with an M(r) of 35,554. The convergently oriented EcoVIII methyltransferase (M. EcoVIII) gene (ecoVIIIM) consists of 912 bp that code for a 304-amino-acid protein with an M(r) of 33,930. The exact positions of the start codon AUG were determined by protein microsequencing. Both enzymes recognize the specific palindromic sequence 5'-AAGCTT-3'. Preparations of EcoVIII R-M enzymes purified to homogeneity were characterized. R. EcoVIII acts as a dimer and cleaves a specific sequence between two adenine residues, leaving 4-nucleotide 5' protruding ends. M. EcoVIII functions as a monomer and modifies the first adenine residue at the 5' end of the specific sequence to N(6)-methyladenine. These enzymes are thus functionally identical to the corresponding enzymes of the HindIII (Haemophilus influenzae Rd) and LlaCI (Lactococcus lactis subsp. cremoris W15) R-M systems. This finding is reflected by the levels of homology of M. EcoVIII with M. HindIII and M. LlaCI at the amino acid sequence level (50 and 62%, respectively) and by the presence of nine sequence motifs conserved among m(6) N-adenine beta-class methyltransferases. The deduced amino acid sequence of R. EcoVIII shows weak homology with its two isoschizomers, R. HindIII (26%) and R. LlaCI (17%). A catalytic sequence motif characteristic of restriction endonucleases was found in the primary structure of R. EcoVIII (D(108)X(12)DXK(123)), as well as in the primary structures of R. LlaCI and R. HindIII. Polyclonal antibodies raised against R. EcoVIII did not react with R. HindIII, while anti-M. EcoVIII antibodies cross-reacted with M. LlaCI but not with M. HindIII. R. EcoVIII requires Mg(II) ions for phosphodiester bond cleavage. We found that the same ions are strong inhibitors of the M. EcoVIII enzyme. The biological implications of this finding are discussed.

Mobility of a restriction-modification system revealed by its genetic contexts in three hosts

The flow of genes among prokaryotes plays a fundamental role in shaping bacterial evolution, and restriction-modification systems can modulate this flow. However, relatively little is known about the distribution and movement of restriction-modification systems themselves. We have isolated and characterized the genes for restriction-modification systems from two species of Salmonella, S. enterica serovar Paratyphi A and S. enterica serovar Bareilly. Both systems are closely related to the PvuII restriction-modification system and share its target specificity. In the case of S. enterica serovar Paratyphi A, the restriction endonuclease is inactive, apparently due to a mutation in the subunit interface region. Unlike the chromosomally located Salmonella systems, the PvuII system is plasmid borne. We have completed the sequence characterization of the PvuII plasmid pPvu1, originally from Proteus vulgaris, making this the first completely sequenced plasmid from the genus Proteus. Despite the pronounced similarity of the three restriction-modification systems, the flanking sequences in Proteus and Salmonella are completely different. The SptAI and SbaI genes lie between an equivalent pair of bacteriophage P4-related open reading frames, one of which is a putative integrase gene, while the PvuII genes are adjacent to a mob operon and a XerCD recombination (cer) site.