A nomenclature for restriction enzymes, DNA methyltransferases, homing endonucleases and their genes

Dam and Dcm methylations prevent gene transfer into Clostridium pasteurianum NRRL B-598: development of methods for electrotransformation, conjugation, and sonoporation

BACKGROUND

Butanol is currently one of the most discussed biofuels. Its use provides many benefits in comparison to bio-ethanol, but the price of its fermentative production is still high. Genetic improvements could help solve many problems associated with butanol production during ABE fermentation, such as its toxicity, low concentration achievable in the cultivation medium, the need for a relatively expensive substrate, and many more. Clostridium pasteurianum NRRL B-598 is non-type strain producing butanol, acetone, and a negligible amount of ethanol. Its main benefits are high oxygen tolerance, utilization of a wide range of carbon and nitrogen sources, and the availability of its whole genome sequence. However, there is no established method for the transfer of foreign DNA into this strain; this is the next step necessary for progress in its use for butanol production.

RESULTS

We have described functional protocols for conjugation and transformation of the bio-butanol producer C. pasteurianum NRRL B-598 by foreign plasmid DNA. We show that the use of unmethylated plasmid DNA is necessary for efficient transformation or successful conjugation. Genes encoding DNA methylation and those for restriction-modification systems and antibiotic resistance were searched for in the whole genome sequence and their homologies with other clostridial bacteria were determined. Furthermore, activity of described novel type I restriction system was proved experimentally. The described electrotransformation protocol achieved an efficiency 1.2 × 10(2) cfu/μg DNA after step-by-step optimization and an efficiency of 1.6 × 10(2) cfu/μg DNA was achieved by the sonoporation technique using a standard laboratory ultrasound bath. The highest transformation efficiency was achieved using a combination of these approaches; sono/electroporation led to an increase in transformation efficiency, to 5.3 × 10(2) cfu/μg DNA.

CONCLUSIONS

Both Dam and Dcm methylations are detrimental for transformation of C. pasteurianum NRRL B-598. Methods for conjugation, electroporation, sonoporation, and a combined method for sono/electroporation were established for this strain. The methods described could be used for genetic improvement of this strain, which is suitable for bio-butanol production.

SMRT sequencing of the Campylobacter coli BfR-CA-9557 genome sequence reveals unique methylation motifs

BACKGROUND

Campylobacter species are the most prevalent bacterial pathogen causing acute enteritis worldwide. In contrast to Campylobacter jejuni, about 5 % of Campylobacter coli strains exhibit susceptibility to restriction endonuclease digestion by DpnI cutting specifically 5'-G(m)ATC-3' motifs. This indicates significant differences in DNA methylation between both microbial species. The goal of the study was to analyze the methylome of a C. coli strain susceptible to DpnI digestion, to identify its methylation motifs and restriction modification systems (RM-systems), and compare them to related organisms like C. jejuni and Helicobacter pylori.

RESULTS

Using one SMRT cell and the PacBio RS sequencing technology followed by PacBio Modification and Motif Analysis the complete genome of the DpnI susceptible strain C. coli BfR-CA-9557 was sequenced to 500-fold coverage and assembled into a single contig of 1.7 Mbp. The genome contains a CJIE1-like element prophage and is phylogenetically closer to C. coli clade 1 isolates than clade 3. 45,881 6-methylated adenines (ca. 2.7 % of genome positions) that are predominantly arranged in eight different methylation motifs and 1,788 4-methylated cytosines (ca. 0.1 %) have been detected. Only two of these motifs correspond to known restriction modification motifs. Characteristic for this methylome was the very high fraction of methylation of motifs with mostly above 99 %.

CONCLUSIONS

Only five dominant methylation motifs have been identified in C. jejuni, which have been associated with known RM-systems. C. coli BFR-CA-9557 shares one (RAATTY) of these, but four ORFs could be assigned to putative Type I RM-systems, seven ORFs to Type II RM-systems and three ORFs to Type IV RM-systems. In accordance with DpnI prescreening RM-system IIP, methylation of GATC motifs was detected in C. coli BfR-CA-9557. A homologous IIP RM-system has been described for H. pylori. The remaining methylation motifs are specific for C. coli BfR-CA-9557 and have been neither detected in C. jejuni nor in H. pylori. The results of this study give us new insights into epigenetics of Campylobacteraceae and provide the groundwork to resolve the function of RM-systems in C. coli.

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.

Adenine methylation in eukaryotes: Apprehending the complex evolutionary history and functional potential of an epigenetic modification

While N⁶‐methyladenosine (m⁶A) is a well‐known epigenetic modification in bacterial DNA, it remained largely unstudied in eukaryotes. Recent studies have brought to fore its potential epigenetic role across diverse eukaryotes with biological consequences, which are distinct and possibly even opposite to the well‐studied 5‐methylcytosine mark. Adenine methyltransferases appear to have been independently acquired by eukaryotes on at least 13 occasions from prokaryotic restriction‐modification and counter‐restriction systems. On at least four to five instances, these methyltransferases were recruited as RNA methylases. Thus, m⁶A marks in eukaryotic DNA and RNA might be more widespread and diversified than previously believed. Several m⁶A‐binding protein domains from prokaryotes were also acquired by eukaryotes, facilitating prediction of potential readers for these marks. Further, multiple lineages of the AlkB family of dioxygenases have been recruited as m⁶A demethylases. Although members of the TET/JBP family of dioxygenases have also been suggested to be m⁶A demethylases, this proposal needs more careful evaluation.

Also watch the Video Abstract.

Whole genome sequencing provides possible explanations for the difference in phage susceptibility among two Salmonella Typhimurium phage types (DT8 and DT30) associated with a single foodborne outbreak

BACKGROUND

Phage typing has been used for decades as a rapid, low cost approach for the epidemiological surveillance of Salmonella enterica subsp. enterica serovar Typhimurium. Although molecular methods are replacing phage typing the system is still in use and provides a valuable model for study of phage-host interaction. Phage typing depends on the pattern of bacterial resistance or sensitivity to a panel of specific bacteriophages. In the phage typing scheme, S. Typhimurium definitive phage types (DT) 8 and 30 differ greatly in their susceptibility to the 30 typing phages of S. Typhimurium; DT8 is susceptible to 11 phages whereas DT30 is resistant to all typing phages except one phage although both DT8 and DT30 were reported to be associated with a single foodborne salmonellosis outbreak in Ireland between 2009 and 2011. We wished to study the genomic correlates of the DT8 and DT30 difference in phage susceptibility using the whole genome sequence (WGS) of S. Typhimurium DT8 and DT30 representatives.

RESULTS

Comparative genome analysis revealed that both S. Typhimurium DT8 and DT30 are lysogenic for three prophages including two S. Typhimurium associated prophages (Gifsy-2 and ST64B) and one S. Enteritidis associated prophage (Enteritidis lysogenic phage S) which has not been detected previously in S. Typhimurium. Furthermore, DT8 and DT30 contain identical clustered regularly interspaced short palindromic repeats (CRISPRs). Interestingly, S. Typhimurium DT8 harbours an accessory genome represented by a virulence plasmid that is highly related to the pSLT plasmid of S. Typhimurium strain LT2 (phage typed as DT4) and codes a unique methyltransferase (MTase); M.EcoGIX related MTase. This plasmid is not detected in DT30. On the other hand, DT30 carries a unique genomic island similar to the integrative and conjugative element (ICE) of Enterotoxigenic Escherichia coli (ETEC) and encodes type IV secretion pathway system (T4SS) and several hypothetical proteins. This genomic island is not detected in DT8.

CONCLUSIONS

We suggest that differences in phage susceptibility between DT8 and DT30 may be related to acquisition of ICE in DT30 and loss of pSLT like plasmid that might be associated with DT30 resistance to almost all phages used in the typing scheme. Additional studies are required to determine the significance of the differences among DT8 and DT30 in relation to the difference in phage susceptibility. This study represents an initial step toward understanding the molecular basis of this host-phage relationship.

Lineage-Specific Methyltransferases Define the Methylome of the Globally Disseminated Escherichia coli ST131 Clone

Escherichia coli sequence type 131 (ST131) is a clone of uropathogenic E. coli that has emerged rapidly and disseminated globally in both clinical and community settings. Members of the ST131 lineage from across the globe have been comprehensively characterized in terms of antibiotic resistance, virulence potential, and pathogenicity, but to date nothing is known about the methylome of these important human pathogens. Here we used single-molecule real-time (SMRT) PacBio sequencing to determine the methylome of E. coli EC958, the most-well-characterized completely sequenced ST131 strain. Our analysis of 52,081 methylated adenines in the genome of EC958 discovered three (m6)A methylation motifs that have not been described previously. Subsequent SMRT sequencing of isogenic knockout mutants identified the two type I methyltransferases (MTases) and one type IIG MTase responsible for (m6)A methylation of novel recognition sites. Although both type I sites were rare, the type IIG sites accounted for more than 12% of all methylated adenines in EC958. Analysis of the distribution of MTase genes across 95 ST131 genomes revealed their prevalence is highly conserved within the ST131 lineage, with most variation due to the presence or absence of mobile genetic elements on which individual MTase genes are located.

IMPORTANCE

DNA modification plays a crucial role in bacterial regulation. Despite several examples demonstrating the role of methyltransferase (MTase) enzymes in bacterial virulence, investigation of this phenomenon on a whole-genome scale has remained elusive until now. Here we used single-molecule real-time (SMRT) sequencing to determine the first complete methylome of a strain from the multidrug-resistant E. coli sequence type 131 (ST131) lineage. By interrogating the methylome computationally and with further SMRT sequencing of isogenic mutants representing previously uncharacterized MTase genes, we defined the target sequences of three novel ST131-specific MTases and determined the genomic distribution of all MTase target sequences. Using a large collection of 95 previously sequenced ST131 genomes, we identified mobile genetic elements as a major factor driving diversity in DNA methylation patterns. Overall, our analysis highlights the potential for DNA methylation to dramatically influence gene regulation at the transcriptional level within a well-defined E. coli clone.

Phylogenomics and sequence-structure-function relationships in the GmrSD family of Type IV restriction enzymes

BACKGROUND

GmrSD is a modification-dependent restriction endonuclease that specifically targets and cleaves glucosylated hydroxymethylcytosine (glc-HMC) modified DNA. It is encoded either as two separate single-domain GmrS and GmrD proteins or as a single protein carrying both domains. Previous studies suggested that GmrS acts as endonuclease and NTPase whereas GmrD binds DNA.

METHODS

In this work we applied homology detection, sequence conservation analysis, fold recognition and homology modeling methods to study sequence-structure-function relationships in the GmrSD restriction endonucleases family. We also analyzed the phylogeny and genomic context of the family members.

RESULTS

Results of our comparative genomics study show that GmrS exhibits similarity to proteins from the ParB/Srx fold which can have both NTPase and nuclease activity. In contrast to the previous studies though, we attribute the nuclease activity also to GmrD as we found it to contain the HNH endonuclease motif. We revealed residues potentially important for structure and function in both domains. Moreover, we found that GmrSD systems exist predominantly as a fused, double-domain form rather than as a heterodimer and that their homologs are often encoded in regions enriched in defense and gene mobility-related elements. Finally, phylogenetic reconstructions of GmrS and GmrD domains revealed that they coevolved and only few GmrSD systems appear to be assembled from distantly related GmrS and GmrD components.

CONCLUSIONS

Our study provides insight into sequence-structure-function relationships in the yet poorly characterized family of Type IV restriction enzymes. Comparative genomics allowed to propose possible role of GmrD domain in the function of the GmrSD enzyme and possible active sites of both GmrS and GmrD domains. Presented results can guide further experimental characterization of these enzymes.

MmoSTI restriction endonuclease, isolated from Morganella morganii infecting a tropical moth, Actias selene, cleaving 5'-|CCNGG-3' sequences

A type II restriction endonuclease, MmoSTI, from the pathogenic bacterium Morganella morganii infecting a tropical moth, Actias selene, has been detected and biochemically characterized, as a potential etiological differentiation factor. The described REase recognizes interrupted palindromes, i.e., 5′-CCNGG-3′ sequences and cleaves DNA leaving 5-nucleotide (nt) long, single-stranded (ss), 5′-cohesive ends, which was determined by three complementary methods: (i) cleavage of custom and standard DNA substrates, (ii) run-off sequencing of cleavage products, and (iii) shotgun cloning and sequencing of bacteriophage lambda (λ) DNA digested with MmoSTI. MmoSTI, the first 5′-CCNGG-3′ REase characterized from M. morganii, is a neoschizomer of ScrFI, which cleaves DNA leaving 1-nt long, ss, 5′-cohesive ends. It is a high-frequency cutter and can be isolated from easily cultured bacteria, thus it can potentially serve as a tool for DNA manipulations.

Functional significance of protein assemblies predicted by the crystal structure of the restriction endonuclease BsaWI

Type II restriction endonuclease BsaWI recognizes a degenerated sequence 5′-W/CCGGW-3′ (W stands for A or T, ‘/’ denotes the cleavage site). It belongs to a large family of restriction enzymes that contain a conserved CCGG tetranucleotide in their target sites. These enzymes are arranged as dimers or tetramers, and require binding of one, two or three DNA targets for their optimal catalytic activity. Here, we present a crystal structure and biochemical characterization of the restriction endonuclease BsaWI. BsaWI is arranged as an ‘open’ configuration dimer and binds a single DNA copy through a minor groove contacts. In the crystal primary BsaWI dimers form an indefinite linear chain via the C-terminal domain contacts implying possible higher order aggregates. We show that in solution BsaWI protein exists in a dimer-tetramer-oligomer equilibrium, but in the presence of specific DNA forms a tetramer bound to two target sites. Site-directed mutagenesis and kinetic experiments show that BsaWI is active as a tetramer and requires two target sites for optimal activity. We propose BsaWI mechanism that shares common features both with dimeric Ecl18kI/SgrAI and bona fide tetrameric NgoMIV/SfiI enzymes.

Genetic Stabilization of the Drug-Resistant PMEN1 Pneumococcus Lineage by Its Distinctive DpnIII Restriction-Modification System

The human pathogen Streptococcus pneumoniae (pneumococcus) exhibits a high degree of genomic diversity and plasticity. Isolates with high genomic similarity are grouped into lineages that undergo homologous recombination at variable rates. PMEN1 is a pandemic, multidrug-resistant lineage. Heterologous gene exchange between PMEN1 and non-PMEN1 isolates is directional, with extensive gene transfer from PMEN1 strains and only modest transfer into PMEN1 strains. Restriction-modification (R-M) systems can restrict horizontal gene transfer, yet most pneumococcal strains code for either the DpnI or DpnII R-M system and neither limits homologous recombination. Our comparative genomic analysis revealed that PMEN1 isolates code for DpnIII, a third R-M system syntenic to the other Dpn systems. Characterization of DpnIII demonstrated that the endonuclease cleaves unmethylated double-stranded DNA at the tetramer sequence 5′ GATC 3′, and the cognate methylase is a C5 cytosine-specific DNA methylase. We show that DpnIII decreases the frequency of recombination under in vitro conditions, such that the number of transformants is lower for strains transformed with unmethylated DNA than in those transformed with cognately methylated DNA. Furthermore, we have identified two PMEN1 isolates where the DpnIII endonuclease is disrupted, and phylogenetic work by Croucher and colleagues suggests that these strains have accumulated genomic differences at a higher rate than other PMEN1 strains. We propose that the R-M locus is a major determinant of genetic acquisition; the resident R-M system governs the extent of genome plasticity.

Pneumococcus is one of the most important community-acquired bacterial pathogens. Pneumococcal strains can develop resistance to antibiotics and to serotype vaccines by acquiring genes from other strains or species. Thus, genomic plasticity is associated with strain adaptability and pneumococcal success. PMEN1 is a widespread and multidrug-resistant highly pathogenic pneumococcal lineage, which has evolved over the past century and displays a relatively stable genome. In this study, we characterize DpnIII, a restriction-modification (R-M) system that limits recombination. DpnIII is encountered in the PMEN1 lineage, where it replaces other R-M systems that do not decrease plasticity. Our hypothesis is that this genomic region, where different pneumococcal lineages code for variable R-M systems, plays a role in the fine-tuning of the extent of genomic plasticity. It is possible that well-adapted lineages such as PMEN1 have a mechanism to increase genomic stability, rather than foster genomic plasticity.

The complete methylome of Helicobacter pylori UM032

BACKGROUND

The genome of the human gastric pathogen Helicobacter pylori encodes a large number of DNA methyltransferases (MTases), some of which are shared among many strains, and others of which are unique to a given strain. The MTases have potential roles in the survival of the bacterium. In this study, we sequenced a Malaysian H. pylori clinical strain, designated UM032, by using a combination of PacBio Single Molecule, Real-Time (SMRT) and Illumina MiSeq next generation sequencing platforms, and used the SMRT data to characterize the set of methylated bases (the methylome).

RESULTS

The N4-methylcytosine and N6-methyladenine modifications detected at single-base resolution using SMRT technology revealed 17 methylated sequence motifs corresponding to one Type I and 16 Type II restriction-modification (R-M) systems. Previously unassigned methylation motifs were now assigned to their respective MTases-coding genes. Furthermore, one gene that appears to be inactive in the H. pylori UM032 genome during normal growth was characterized by cloning.

CONCLUSION

Consistent with previously-studied H. pylori strains, we show that strain UM032 contains a relatively large number of R-M systems, including some MTase activities with novel specificities. Additional studies are underway to further elucidating the biological significance of the R-M systems in the physiology and pathogenesis of H. pylori.

Comparative genomics reveals diversified CRISPR-Cas systems of globally distributed Microcystis aeruginosa, a freshwater bloom-forming cyanobacterium

Microcystis aeruginosa is one of the most common and dominant bloom-forming cyanobacteria in freshwater lakes around the world. Microcystis cells can produce toxic secondary metabolites, such as microcystins, which are harmful to human health. Two M. aeruginosa strains were isolated from two highly eutrophic lakes in China and their genomes were sequenced. Comparative genomic analysis was performed with the 12 other available M. aeruginosa genomes and closely related unicellular cyanobacterium. Each genome of M. aeruginosa containing at least one clustered regularly interspaced short palindromic repeat (CRISPR) locus and total 71 loci were identified, suggesting it is ubiquitous in M. aeruginosa genomes. In addition to the previously reported subtype I-D cas gene sets, three CAS subtypes I-A, III-A and III-B were identified and characterized in this study. Seven types of CRISPR direct repeat have close association with CAS subtype, confirming that different and specific secondary structures of CRISPR repeats are important for the recognition, binding and process of corresponding cas gene sets. Homology search of the CRISPR spacer sequences provides a history of not only resistance to bacteriophages and plasmids known to be associated with M. aeruginosa, but also the ability to target much more exogenous genetic material in the natural environment. These adaptive and heritable defense mechanisms play a vital role in keeping genomic stability and self-maintenance by restriction of horizontal gene transfer. Maintaining genomic stability and modulating genomic plasticity are both important evolutionary strategies for M. aeruginosa in adaptation and survival in various habitats.

Complete Genome Sequence Analysis of Bacillus subtilis T30

The complete genome sequence of Bacillus subtilis T30 was determined by SMRT sequencing. The entire genome contains 4,138 predicted genes. The genome carries one intact prophage sequence (37.4 kb) similar to Bacillus phage SPBc2 and one incomplete prophage genome of 39.9 kb similar to Bacillus phage phi105.

The TET/JBP Family of Nucleic Acid Base-Modifying 2-Oxoglutarate and Iron-Dependent Dioxygenases

The TET/JBP family of enzymes includes 2-oxoglutarate- and Fe(ii)-dependent dioxygenases that oxidize 5-methylpyrimidines in nucleic acids. They include euglenozoan JBP enzymes that catalyse the first step in the biosynthesis of the hypermodified thymine, base J, and metazoan TET enzymes that generate oxidized 5-methylcytosines (hydroxy-, formyl- and carboxymethylcytosine) in DNA. Recent studies suggest that these modified bases function as epigenetic marks and/or as potential intermediates for DNA demethylation during resetting of epigenetic 5mC marks upon zygote formation and in primordial germ cell development. Studies in mammalian models also point to an important role for these enzymes in haematopoiesis, tumour suppression, cell differentiation and neural behavioural adaptation. The TET/JBP family has undergone extensive gene expansion in fungi, such as mushrooms, in conjunction with a novel class of transposons and might play a role in genomic plasticity and speciation. Certain versions from stramenopiles and chlorophytes are likely to modify RNA and often show fusions to other RNA-modifying enzymatic domains. The ultimate origin of the TET/JBP family lies in bacteriophages where the enzymes are likely to catalyse formation of modified bases with key roles in DNA packaging and evasion of host restriction.

Specificity of the ModA11, ModA12 and ModD1 epigenetic regulator N(6)-adenine DNA methyltransferases of Neisseria meningitidis

Phase variation (random ON/OFF switching) of gene expression is a common feature of host-adapted pathogenic bacteria. Phase variably expressed N(6)-adenine DNA methyltransferases (Mod) alter global methylation patterns resulting in changes in gene expression. These systems constitute phase variable regulons called phasevarions. Neisseria meningitidis phasevarions regulate genes including virulence factors and vaccine candidates, and alter phenotypes including antibiotic resistance. The target site recognized by these Type III N(6)-adenine DNA methyltransferases is not known. Single molecule, real-time (SMRT) methylome analysis was used to identify the recognition site for three key N. meningitidis methyltransferases: ModA11 (exemplified by M.NmeMC58I) (5'-CGY M6A: G-3'), ModA12 (exemplified by M.Nme77I, M.Nme18I and M.Nme579II) (5'-AC M6A: CC-3') and ModD1 (exemplified by M.Nme579I) (5'-CC M6A: GC-3'). Restriction inhibition assays and mutagenesis confirmed the SMRT methylome analysis. The ModA11 site is complex and atypical and is dependent on the type of pyrimidine at the central position, in combination with the bases flanking the core recognition sequence 5'-CGY M6A: G-3'. The observed efficiency of methylation in the modA11 strain (MC58) genome ranged from 4.6% at 5'-GCGC M6A: GG-3' sites, to 100% at 5'-ACGT M6A: GG-3' sites. Analysis of the distribution of modified sites in the respective genomes shows many cases of association with intergenic regions of genes with altered expression due to phasevarion switching.

Next-generation sequencing as an approach to dairy starter selection

Lactococcal and streptococcal starter strains are crucial ingredients to manufacture fermented dairy products. As commercial starter culture suppliers and dairy producers attempt to overcome issues of phage sensitivity and develop new product ranges, there is an ever increasing need to improve technologies for the rational selection of novel starter culture blends. Whole genome sequencing, spurred on by recent advances in next-generation sequencing platforms, is a promising approach to facilitate rapid identification and selection of such strains based on gene-trait matching. This review provides a comprehensive overview of the available methodologies to analyse the technological potential of candidate starter strains and highlights recent advances in the area of dairy starter genomics.

Genome-wide methylation patterns in Salmonella enterica Subsp. enterica Serovars

The methylation of DNA bases plays an important role in numerous biological processes including development, gene expression, and DNA replication. Salmonella is an important foodborne pathogen, and methylation in Salmonella is implicated in virulence. Using single molecule real-time (SMRT) DNA-sequencing, we sequenced and assembled the complete genomes of eleven Salmonella enterica isolates from nine different serovars, and analysed the whole-genome methylation patterns of each genome. We describe 16 distinct N⁶-methyladenine (m6A) methylated motifs, one N⁴-methylcytosine (m4C) motif, and one combined m6A-m4C motif. Eight of these motifs are novel, i.e., they have not been previously described. We also identified the methyltransferases (MTases) associated with 13 of the motifs. Some motifs are conserved across all Salmonella serovars tested, while others were found only in a subset of serovars. Eight of the nine serovars contained a unique methylated motif that was not found in any other serovar (most of these motifs were part of Type I restriction modification systems), indicating the high diversity of methylation patterns present in Salmonella.

Genome-wide DNA methylation analysis of Haloferax volcanii H26 and identification of DNA methyltransferase related PD-(D/E)XK nuclease family protein HVO_A0006

Restriction-modification (RM) systems have evolved to protect the cell from invading DNAs and are composed of two enzymes: a DNA methyltransferase and a restriction endonuclease. Although RM systems are present in both archaeal and bacterial genomes, DNA methylation in archaea has not been well defined. In order to characterize the function of RM systems in archaeal species, we have made use of the model haloarchaeon Haloferax volcanii. A genomic DNA methylation analysis of H. volcanii strain H26 was performed using PacBio single molecule real-time (SMRT) sequencing. This analysis was also performed on a strain of H. volcanii in which an annotated DNA methyltransferase gene HVO_A0006 was deleted from the genome. Sequence analysis of H26 revealed two motifs which are modified in the genome: C(m4)TAG and GCA(m6)BN6VTGC. Analysis of the ΔHVO_A0006 strain indicated that it exhibited reduced adenine methylation compared to the parental strain and altered the detected adenine motif. However, protein domain architecture analysis and amino acid alignments revealed that HVO_A0006 is homologous only to the N-terminal endonuclease region of Type IIG RM proteins and contains a PD-(D/E)XK nuclease motif, suggesting that HVO_A0006 is a PD-(D/E)XK nuclease family protein. Further bioinformatic analysis of the HVO_A0006 gene demonstrated that the gene is rare among the Halobacteria. It is surrounded by two transposition genes suggesting that HVO_A0006 is a fragment of a Type IIG RM gene, which has likely been acquired through gene transfer, and affects restriction-modification activity by interacting with another RM system component(s). Here, we present the first genome-wide characterization of DNA methylation in an archaeal species and examine the function of a DNA methyltransferase related gene HVO_A0006.

Genome Modification in Enterococcus faecalis OG1RF Assessed by Bisulfite Sequencing and Single-Molecule Real-Time Sequencing

Enterococcus faecalis is a Gram-positive bacterium that natively colonizes the human gastrointestinal tract and opportunistically causes life-threatening infections. Multidrug-resistant (MDR) E. faecalis strains have emerged, reducing treatment options for these infections. MDR E. faecalis strains have large genomes containing mobile genetic elements (MGEs) that harbor genes for antibiotic resistance and virulence determinants. Bacteria commonly possess genome defense mechanisms to block MGE acquisition, and we hypothesize that these mechanisms have been compromised in MDR E. faecalis. In restriction-modification (R-M) defense, the bacterial genome is methylated at cytosine (C) or adenine (A) residues by a methyltransferase (MTase), such that nonself DNA can be distinguished from self DNA. A cognate restriction endonuclease digests improperly modified nonself DNA. Little is known about R-M in E. faecalis. Here, we use genome resequencing to identify DNA modifications occurring in the oral isolate OG1RF. OG1RF has one of the smallest E. faecalis genomes sequenced to date and possesses few MGEs. Single-molecule real-time (SMRT) and bisulfite sequencing revealed that OG1RF has global 5-methylcytosine (m5C) methylation at 5'-GCWGC-3' motifs. A type II R-M system confers the m5C modification, and disruption of this system impacts OG1RF electrotransformability and conjugative transfer of an antibiotic resistance plasmid. A second DNA MTase was poorly expressed under laboratory conditions but conferred global N(4)-methylcytosine (m4C) methylation at 5'-CCGG-3' motifs when expressed in Escherichia coli. Based on our results, we conclude that R-M can act as a barrier to MGE acquisition and likely influences antibiotic resistance gene dissemination in the E. faecalis species.

IMPORTANCE

The horizontal transfer of antibiotic resistance genes among bacteria is a critical public health concern. Enterococcus faecalis is an opportunistic pathogen that causes life-threatening infections in humans. Multidrug resistance acquired by horizontal gene transfer limits treatment options for these infections. In this study, we used innovative DNA sequencing methodologies to investigate how a model strain of E. faecalis discriminates its own DNA from foreign DNA, i.e., self versus nonself discrimination. We also assess the role of an E. faecalis genome modification system in modulating conjugative transfer of an antibiotic resistance plasmid. These results are significant because they demonstrate that differential genome modification impacts horizontal gene transfer frequencies in E. faecalis.

Improving transformation of Staphylococcus aureus belonging to the CC1, CC5 and CC8 clonal complexes

Methicillin resistant Staphylococcus aureus (MRSA) is an opportunistic pathogen found in hospital and community environments that can cause serious infections. A major barrier to genetic manipulations of clinical isolates has been the considerable difficulty in transforming these strains with foreign plasmids, such as those from E. coli, in part due to the type I and IV Restriction Modification (R-M) barriers. Here we combine a Plasmid Artificial Modification (PAM) system with DC10B E. coli cells (dcm mutants) to bypass the barriers of both type I and IV R-M of S. aureus, thus allowing E. coli plasmid DNA to be transformed directly into clinical MRSA strains MW2, N315 and LAC, representing three of the most common clonal complexes. Successful transformation of clinical S. aureus isolates with E. coli-derived plasmids should greatly increase the ability to genetically modify relevant S. aureus strains and advance our understanding of S. aureus pathogenesis.

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.

Restriction-modification system with methyl-inhibited base excision and abasic-site cleavage activities

The restriction-modification systems use epigenetic modification to distinguish between self and nonself DNA. A modification enzyme transfers a methyl group to a base in a specific DNA sequence while its cognate restriction enzyme introduces breaks in DNA lacking this methyl group. So far, all the restriction enzymes hydrolyze phosphodiester bonds linking the monomer units of DNA. We recently reported that a restriction enzyme (R.PabI) of the PabI superfamily with half-pipe fold has DNA glycosylase activity that excises an adenine base in the recognition sequence (5′-GTAC). We now found a second activity in this enzyme: at the resulting apurinic/apyrimidinic (AP) (abasic) site (5′-GT#C, # = AP), its AP lyase activity generates an atypical strand break. Although the lyase activity is weak and lacks sequence specificity, its covalent DNA–R.PabI reaction intermediates can be trapped by NaBH₄ reduction. The base excision is not coupled with the strand breakage and yet causes restriction because the restriction enzyme action can impair transformation ability of unmethylated DNA even in the absence of strand breaks in vitro. The base excision of R.PabI is inhibited by methylation of the target adenine base. These findings expand our understanding of genetic and epigenetic processes linking those in prokaryotes and eukaryotes.

A comparative analysis of methylome profiles of Campylobacter jejuni sheep abortion isolate and gastroenteric strains using PacBio data

Campylobacter jejuni is a leading cause of human gastrointestinal disease and small ruminant abortions in the United States. The recent emergence of a highly virulent, tetracycline-resistant C. jejuni subsp. jejuni sheep abortion clone (clone SA) in the United States, and that strain's association with human disease, has resulted in a heightened awareness of the zoonotic potential of this organism. Pacific Biosciences' Single Molecule, Real-Time sequencing technology was used to explore the variation in the genome-wide methylation patterns of the abortifacient clone SA (IA3902) and phenotypically distinct gastrointestinal-specific C. jejuni strains (NCTC 11168 and 81-176). Several notable differences were discovered that distinguished the methylome of IA3902 from that of 11168 and 81-176: identification of motifs novel to IA3902, genome-specific hypo- and hypermethylated regions, strain level variability in genes methylated, and differences in the types of methylation motifs present in each strain. These observations suggest a possible role of methylation in the contrasting disease presentations of these three C. jejuni strains. In addition, the methylation profiles between IA3902 and a luxS mutant were explored to determine if variations in methylation patterns could be identified that might explain the role of LuxS-dependent methyl recycling in IA3902 abortifacient potential.

Genomic insights into the rapid emergence and evolution of MDR in Staphylococcus pseudintermedius

OBJECTIVES

MDR methicillin-resistant Staphylococcus pseudintermedius (MRSP) strains have emerged rapidly as major canine pathogens and present serious treatment issues and concerns to public health due to their, albeit low, zoonotic potential. A further understanding of the genetics of resistance arising from a broadly susceptible background of S. pseudintermedius is needed.

METHODS

We sequenced the genomes of 12 S. pseudintermedius isolates of varied STs and resistance phenotypes.

RESULTS

Nine distinct clonal lineages had acquired either staphylococcal cassette chromosome (SCC) mec elements and/or Tn5405-like elements carrying up to five resistance genes [aphA3, sat, aadE, erm(B), dfrG] to generate MRSP, MDR methicillin-susceptible S. pseudintermedius and MDR MRSP populations. The most successful and clinically problematic MDR MRSP clones, ST68 SCCmecV(T) and ST71 SCCmecII-III, have further accumulated mutations in gyrA and grlA conferring resistance to fluoroquinolones. The carriage of additional mobile genetic elements (MGEs) was highly variable, suggesting that horizontal gene transfer is frequent in S. pseudintermedius populations.

CONCLUSIONS

Importantly, the data suggest that MDR MRSP evolved rapidly by the acquisition of a very limited number of MGEs and mutations, and that the use of many classes of antimicrobials may co-select for the spread and emergence of MDR and XDR strains. Antimicrobial stewardship will need to be comprehensive, encompassing human medicine and veterinary disciplines to successfully preserve antimicrobial efficacy.

DNA cleavage by CgII and NgoAVII requires interaction between N- and R-proteins and extensive nucleotide hydrolysis

The stress-sensitive restriction-modification (RM) system CglI from Corynebacterium glutamicum and the homologous NgoAVII RM system from Neisseria gonorrhoeae FA1090 are composed of three genes: a DNA methyltransferase (M.CglI and M.NgoAVII), a putative restriction endonuclease (R.CglI and R.NgoAVII, or R-proteins) and a predicted DEAD-family helicase/ATPase (N.CglI and N.NgoAVII or N-proteins). Here we report a biochemical characterization of the R- and N-proteins. Size-exclusion chromatography and SAXS experiments reveal that the isolated R.CglI, R.NgoAVII and N.CglI proteins form homodimers, while N.NgoAVII is a monomer in solution. Moreover, the R.CglI and N.CglI proteins assemble in a complex with R2N2 stoichiometry. Next, we show that N-proteins have ATPase activity that is dependent on double-stranded DNA and is stimulated by the R-proteins. Functional ATPase activity and extensive ATP hydrolysis (∼170 ATP/s/monomer) are required for site-specific DNA cleavage by R-proteins. We show that ATP-dependent DNA cleavage by R-proteins occurs at fixed positions (6-7 nucleotides) downstream of the asymmetric recognition sequence 5'-GCCGC-3'. Despite similarities to both Type I and II restriction endonucleases, the CglI and NgoAVII enzymes may employ a unique catalytic mechanism for DNA cleavage.

Two-stage gene assembly/cloning of a member of the TspDTI subfamily of bifunctional restriction endonucleases, TthHB27I

The Thermus sp. family of bifunctional type IIS/IIG/IIC restriction endonucleases (REase)-methyltransferases (MTase) comprises thermo-stable TaqII, TspGWI, TspDTI, TsoI, Tth111II/TthHB27I enzymes as well as a number of putative enzymes/open reading frames (ORFs). All of the family members share properties including a large protein size (ca. 120kDa), amino acid (aa) sequence homologies, enzymatic activity modulation by S-adenosylmethionine (SAM), recognition of similar asymmetric cognate DNA sites and cleavage at a distance of 11/9 nt. Analysis of the enzyme aa sequences and domain/motif organisation led to further Thermus sp. family division into the TspDTI and TspGWI subfamilies. The latter exhibits an unprecedented phenomenon of DNA recognition change upon substitution of SAM by its analogue, sinefungin (SIN), towards a very frequent DNA cleavage. We report cloning in Escherichia coli (E. coli), using a two-stage procedure and a putative tthHB27IRM gene, detected by bioinformatics analysis of the Thermus thermophilus HB27 (T. thermophilus) genome. The functionality of a 3366 base pair (bp)-/1121 aa-long, high GC content ORF was validated experimentally through the expression in E. coli. Protein features corroborated with the reclassification of TthHB27I into the TspDTI subfamily, which manifested in terms of aa-sequence/motif homologies and insensitivity to SIN-induced specificity shift. However, both SAM and SIN stimulated the REase DNA cleavage activity by at least 16-32 times; the highest was observed for the Thermus sp. family. The availability of TthHB27I and the need to include SAM or SIN in the reaction in order to convert the enzyme from "hibernation" status to efficient DNA cleavage is of practical significance in molecular biotechnology, extending the palette of available REase specificities.

A new prototype IIS/IIC/IIG endonuclease-methyltransferase TsoI from the thermophile Thermus scotoductus, recognising 5'-TARCCA(N11/9)-3' sequences

The Thermus sp. family of IIS/IIG/IIC enzymes includes the thermostable, bifunctional, fused restriction endonuclease (REase)-methyltransferases (MTase): TaqII, Tth111II/TthHB27I, TspGWI, TspDTI and TsoI. The enzymes are large proteins (approximately 120kDa), their enzymatic activities are affected by S-adenosylmethionine (SAM), they recognise similar asymmetric cognate sites and cleave at a distance of 11/9 nucleotides (nt). The enzymes exhibit similarities of their amino acid (aa) sequences and DNA catalytic motifs. Thermus sp. enzymes are an example of functional aa sequence homologies among REases recognising different, yet related DNA sequences. The family consists of TspGWI- and TspDTI-subfamilies. TsoI appears to be a non-identical 'triplet', related to TspDTI and Tth111II/TthHB27I. The discovery of TsoI, purified from Thermus scotoductus, is described. This prototype, displaying a novel specificity, which was determined by: (i) cleavage of a reference plasmid and bacteriophage DNA, (ii) cleavage of custom PCR DNA substrates, (iii) run-off sequencing of cleavage products and (iv) shotgun cloning and sequencing of bacteriophage lambda (λ) DNA digested with TsoI. The enzyme recognises a degenerated 5'-TARCCA-3' sequence, whereas DNA strands are cut 11/9 nt downstream. The discovery of the TsoI prototype is of practical importance in biotechnology, as it extends the palette of cleavage specificities for gene cloning.

Crystal structure of the R-protein of the multisubunit ATP-dependent restriction endonuclease NgoAVII

The restriction endonuclease (REase) NgoAVII is composed of two proteins, R.NgoAVII and N.NgoAVII, and shares features of both Type II restriction enzymes and Type I/III ATP-dependent restriction enzymes (see accompanying paper Zaremba et al., 2014). Here we present crystal structures of the R.NgoAVII apo-protein and the R.NgoAVII C-terminal domain bound to a specific DNA. R.NgoAVII is composed of two domains: an N-terminal nucleolytic PLD domain; and a C-terminal B3-like DNA-binding domain identified previously in BfiI and EcoRII REases, and in plant transcription factors. Structural comparison of the B3-like domains of R.NgoAVII, EcoRII, BfiI and the plant transcription factors revealed a conserved DNA-binding surface comprised of N- and C-arms that together grip the DNA. The C-arms of R.NgoAVII, EcoRII, BfiI and plant B3 domains are similar in size, but the R.NgoAVII N-arm which makes the majority of the contacts to the target site is much longer. The overall structures of R.NgoAVII and BfiI are similar; however, whilst BfiI has stand-alone catalytic activity, R.NgoAVII requires an auxiliary cognate N.NgoAVII protein and ATP hydrolysis in order to cleave DNA at the target site. The structures we present will help formulate future experiments to explore the molecular mechanisms of intersubunit crosstalk that control DNA cleavage by R.NgoAVII and related endonucleases.

The EcoKI type I restriction-modification system in Escherichia coli affects but is not an absolute barrier for conjugation

The rapid evolution of bacteria is crucial to their survival and is caused by exchange, transfer, and uptake of DNA, among other things. Conjugation is one of the main mechanisms by which bacteria share their DNA, and it is thought to be controlled by varied bacterial immune systems. Contradictory results about restriction-modification systems based on phenotypic studies have been presented as reasons for a barrier to conjugation with and other means of uptake of exogenous DNA. In this study, we show that inactivation of the R.EcoKI restriction enzyme in strain Escherichia coli K-12 strain MG1655 increases the conjugational transfer of plasmid pOLA52, which carriers two EcoKI recognition sites. Interestingly, the results were not absolute, and uptake of unmethylated pOLA52 was still observed in the wild-type strain (with an intact hsdR gene) but at a reduction of 85% compared to the uptake of the mutant recipient with a disrupted hsdR gene. This leads to the conclusion that EcoKI restriction-modification affects the uptake of DNA by conjugation but is not a major barrier to plasmid transfer.

REBASE--a database for DNA restriction and modification: enzymes, genes and genomes

REBASE is a comprehensive and fully curated database of information about the components of restriction-modification (RM) systems. It contains fully referenced information about recognition and cleavage sites for both restriction enzymes and methyltransferases as well as commercial availability, methylation sensitivity, crystal and sequence data. All genomes that are completely sequenced are analyzed for RM system components, and with the advent of PacBio sequencing, the recognition sequences of DNA methyltransferases (MTases) are appearing rapidly. Thus, Type I and Type III systems can now be characterized in terms of recognition specificity merely by DNA sequencing. The contents of REBASE may be browsed from the web http://rebase.neb.com and selected compilations can be downloaded by FTP (ftp.neb.com). Monthly updates are also available via email.

Mobility of DNA sequence recognition domains in DNA methyltransferases suggests epigenetics-driven adaptive evolution

DNA methylation is one of the best studied epigenetic modifications observed in prokaryotes as well as eukaryotes. It affects nearby gene expression. Most DNA methylation reactions in prokaryotes are catalyzed by a DNA methyltransferase, the modification enzyme of a restriction-modification (RM) system. Its target recognition domain (TRD) recognizes a specific DNA sequence for methylation. In this commentary, we review recent evidence for movement of TRDs between non-orthologous genes and movement within a gene. These movements are likely mediated by DNA recombination machinery, and are expected to alter the methylation status of a genome. Such alterations potentially lead to changes in global gene expression pattern and various phenotypes. The targets of natural selection in adaptive evolution might be these diverse methylomes rather than diverse genome sequences, the target according to the current paradigm in biology. This “epigenetics-driven adaptive evolution” hypothesis can explain several observations in the evolution of prokaryotes and eukaryotes.

Methyltransferases acquired by lactococcal 936-type phage provide protection against restriction endonuclease activity

Background

So-called 936-type phages are among the most frequently isolated phages in dairy facilities utilising Lactococcus lactis starter cultures. Despite extensive efforts to control phage proliferation and decades of research, these phages continue to negatively impact cheese production in terms of the final product quality and consequently, monetary return.

Results

Whole genome sequencing and in silico analysis of three 936-type phage genomes identified several putative (orphan) methyltransferase (MTase)-encoding genes located within the packaging and replication regions of the genome. Utilising SMRT sequencing, methylome analysis was performed on all three phages, allowing the identification of adenine modifications consistent with N-6 methyladenine sequence methylation, which in some cases could be attributed to these phage-encoded MTases. Heterologous gene expression revealed that M.Phi145I/M.Phi93I and M.Phi93DAM, encoded by genes located within the packaging module, provide protection against the restriction enzymes HphI and DpnII, respectively, representing the first functional MTases identified in members of 936-type phages.

Conclusions

SMRT sequencing technology enabled the identification of the target motifs of MTases encoded by the genomes of three lytic 936-type phages and these MTases represent the first functional MTases identified in this species of phage. The presence of these MTase-encoding genes on 936-type phage genomes is assumed to represent an adaptive response to circumvent host encoded restriction-modification systems thereby increasing the fitness of the phages in a dynamic dairy environment.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-831) contains supplementary material, which is available to authorized users.

Burkholderia pseudomallei sequencing identifies genomic clades with distinct recombination, accessory, and epigenetic profiles

Burkholderia pseudomallei (Bp) is the causative agent of the infectious disease melioidosis. To investigate population diversity, recombination, and horizontal gene transfer in closely related Bp isolates, we performed whole-genome sequencing (WGS) on 106 clinical, animal, and environmental strains from a restricted Asian locale. Whole-genome phylogenies resolved multiple genomic clades of Bp, largely congruent with multilocus sequence typing (MLST). We discovered widespread recombination in the Bp core genome, involving hundreds of regions associated with multiple haplotypes. Highly recombinant regions exhibited functional enrichments that may contribute to virulence. We observed clade-specific patterns of recombination and accessory gene exchange, and provide evidence that this is likely due to ongoing recombination between clade members. Reciprocally, interclade exchanges were rarely observed, suggesting mechanisms restricting gene flow between clades. Interrogation of accessory elements revealed that each clade harbored a distinct complement of restriction-modification (RM) systems, predicted to cause clade-specific patterns of DNA methylation. Using methylome sequencing, we confirmed that representative strains from separate clades indeed exhibit distinct methylation profiles. Finally, using an E. coli system, we demonstrate that Bp RM systems can inhibit uptake of non-self DNA. Our data suggest that RM systems borne on mobile elements, besides preventing foreign DNA invasion, may also contribute to limiting exchanges of genetic material between individuals of the same species. Genomic clades may thus represent functional units of genetic isolation in Bp, modulating intraspecies genetic diversity.

Diversity, ecology and intestinal function of bifidobacteria

The human gastrointestinal tract represents an environment which is a densely populated home for a microbiota that has evolved to positively contribute to host health. At birth the essentially sterile gastrointestinal tract (GIT) is rapidly colonized by microorganisms that originate from the mother and the surrounding environment. Within a short timeframe a microbiota establishes within the (breastfed) infant's GIT where bifidobacteria are among the dominant members, although their numerical dominance disappears following weaning. The numerous health benefits associated with bifidobacteria, and the consequent commercial relevance resulting from their incorporation into functional foods, has led to intensified research aimed at the molecular understanding of claimed probiotic attributes of this genus. In this review we provide the current status on the diversity and ecology of bifidobacteria. In addition, we will discuss the molecular mechanisms that allow this intriguing group of bacteria to colonize and persist in the GIT, so as to facilitate interaction with its host.

The basic building blocks and evolution of CRISPR-CAS systems

CRISPR (clustered regularly interspaced short palindromic repeats)-Cas (CRISPR-associated) is an adaptive immunity system in bacteria and archaea that functions via a distinct self/non-self recognition mechanism that involves unique spacers homologous with viral or plasmid DNA and integrated into the CRISPR loci. Most of the Cas proteins evolve under relaxed purifying selection and some underwent dramatic structural rearrangements during evolution. In many cases, CRISPR-Cas system components are replaced either by homologous or by analogous proteins or domains in some bacterial and archaeal lineages. However, recent advances in comparative sequence analysis, structural studies and experimental data suggest that, despite this remarkable evolutionary plasticity, all CRISPR-Cas systems employ the same architectural and functional principles, and given the conservation of the principal building blocks, share a common ancestry. We review recent advances in the understanding of the evolution and organization of CRISPR-Cas systems. Among other developments, we describe for the first time a group of archaeal cas1 gene homologues that are not associated with CRISPR-Cas loci and are predicted to be involved in functions other than adaptive immunity.

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

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

Type II restriction endonucleases--a historical perspective and more

This article continues the series of Surveys and Summaries on restriction endonucleases (REases) begun this year in Nucleic Acids Research. Here we discuss 'Type II' REases, the kind used for DNA analysis and cloning. We focus on their biochemistry: what they are, what they do, and how they do it. Type II REases are produced by prokaryotes to combat bacteriophages. With extreme accuracy, each recognizes a particular sequence in double-stranded DNA and cleaves at a fixed position within or nearby. The discoveries of these enzymes in the 1970s, and of the uses to which they could be put, have since impacted every corner of the life sciences. They became the enabling tools of molecular biology, genetics and biotechnology, and made analysis at the most fundamental levels routine. Hundreds of different REases have been discovered and are available commercially. Their genes have been cloned, sequenced and overexpressed. Most have been characterized to some extent, but few have been studied in depth. Here, we describe the original discoveries in this field, and the properties of the first Type II REases investigated. We discuss the mechanisms of sequence recognition and catalysis, and the varied oligomeric modes in which Type II REases act. We describe the surprising heterogeneity revealed by comparisons of their sequences and structures.

Structural basis of the methylation specificity of R.DpnI

R.DpnI consists of N-terminal catalytic and C-terminal winged helix domains that are separately specific for the Gm6ATC sequences in Dam-methylated DNA. Here we present a crystal structure of R.DpnI with oligoduplexes bound to the catalytic and winged helix domains and identify the catalytic domain residues that are involved in interactions with the substrate methyl groups. We show that these methyl groups in the Gm6ATC target sequence are positioned very close to each other. We further show that the presence of the two methyl groups requires a deviation from B-DNA conformation to avoid steric conflict. The methylation compatible DNA conformation is complementary with binding sites of both R.DpnI domains. This indirect readout of methylation adds to the specificity mediated by direct favorable interactions with the methyl groups and solvation/desolvation effects. We also present hydrogen/deuterium exchange data that support ‘crosstalk’ between the two domains in the identification of methylated DNA, which should further enhance R.DpnI methylation specificity.

Structure of 5-hydroxymethylcytosine-specific restriction enzyme, AbaSI, in complex with DNA

AbaSI, a member of the PvuRts1I-family of modification-dependent restriction endonucleases, cleaves deoxyribonucleic acid (DNA) containing 5-hydroxymethylctosine (5hmC) and glucosylated 5hmC (g5hmC), but not DNA containing unmodified cytosine. AbaSI has been used as a tool for mapping the genomic locations of 5hmC, an important epigenetic modification in the DNA of higher organisms. Here we report the crystal structures of AbaSI in the presence and absence of DNA. These structures provide considerable, although incomplete, insight into how this enzyme acts. AbaSI appears to be mainly a homodimer in solution, but interacts with DNA in our structures as a homotetramer. Each AbaSI subunit comprises an N-terminal, Vsr-like, cleavage domain containing a single catalytic site, and a C-terminal, SRA-like, 5hmC-binding domain. Two N-terminal helices mediate most of the homodimer interface. Dimerization brings together the two catalytic sites required for double-strand cleavage, and separates the 5hmC binding-domains by ∼70 Å, consistent with the known activity of AbaSI which cleaves DNA optimally between symmetrically modified cytosines ∼22 bp apart. The eukaryotic SET and RING-associated (SRA) domains bind to DNA containing 5-methylcytosine (5mC) in the hemi-methylated CpG sequence. They make contacts in both the major and minor DNA grooves, and flip the modified cytosine out of the helix into a conserved binding pocket. In contrast, the SRA-like domain of AbaSI, which has no sequence specificity, contacts only the minor DNA groove, and in our current structures the 5hmC remains intra-helical. A conserved, binding pocket is nevertheless present in this domain, suitable for accommodating 5hmC and g5hmC. We consider it likely, therefore, that base-flipping is part of the recognition and cleavage mechanism of AbaSI, but that our structures represent an earlier, pre-flipped stage, prior to actual recognition.

Elimination of inter-domain interactions increases the cleavage fidelity of the restriction endonuclease DraIII

DraIII is a type IIP restriction endonucleases (REases) that recognizes and creates a double strand break within the gapped palindromic sequence CAC↑NNN↓GTG of double-stranded DNA (↑ indicates nicking on the bottom strand; ↓ indicates nicking on the top strand). However, wild type DraIII shows significant star activity. In this study, it was found that the prominent star site is CAT↑GTT↓GTG, consisting of a star 5′ half (CAT) and a canonical 3′ half (GTG). DraIII nicks the 3′ canonical half site at a faster rate than the 5′ star half site, in contrast to the similar rate with the canonical full site. The crystal structure of the DraIII protein was solved. It indicated, as supported by mutagenesis, that DraIII possesses a ββα-metal HNH active site. The structure revealed extensive intra-molecular interactions between the N-terminal domain and the C-terminal domain containing the HNH active site. Disruptions of these interactions through site-directed mutagenesis drastically increased cleavage fidelity. The understanding of fidelity mechanisms will enable generation of high fidelity REases.

Methylome diversification through changes in DNA methyltransferase sequence specificity

Epigenetic modifications such as DNA methylation have large effects on gene expression and genome maintenance. Helicobacter pylori, a human gastric pathogen, has a large number of DNA methyltransferase genes, with different strains having unique repertoires. Previous genome comparisons suggested that these methyltransferases often change DNA sequence specificity through domain movement--the movement between and within genes of coding sequences of target recognition domains. Using single-molecule real-time sequencing technology, which detects N6-methyladenines and N4-methylcytosines with single-base resolution, we studied methylated DNA sites throughout the H. pylori genome for several closely related strains. Overall, the methylome was highly variable among closely related strains. Hypermethylated regions were found, for example, in rpoB gene for RNA polymerase. We identified DNA sequence motifs for methylation and then assigned each of them to a specific homology group of the target recognition domains in the specificity-determining genes for Type I and other restriction-modification systems. These results supported proposed mechanisms for sequence-specificity changes in DNA methyltransferases. Knocking out one of the Type I specificity genes led to transcriptome changes, which suggested its role in gene expression. These results are consistent with the concept of evolution driven by DNA methylation, in which changes in the methylome lead to changes in the transcriptome and potentially to changes in phenotype, providing targets for natural or artificial selection.

Exploring bacterial epigenomics in the next-generation sequencing era: a new approach for an emerging frontier

Epigenetics has an important role for the success of foodborne pathogen persistence in diverse host niches. Substantial challenges exist in determining DNA methylation to situation-specific phenotypic traits. DNA modification, mediated by restriction-modification systems, functions as an immune response against antagonistic external DNA, and bacteriophage-acquired methyltransferases (MTase) and orphan MTases - those lacking the cognate restriction endonuclease - facilitate evolution of new phenotypes via gene expression modulation via DNA and RNA modifications, including methylation and phosphorothioation. Recent establishment of large-scale genome sequencing projects will result in a significant increase in genome availability that will lead to new demands for data analysis including new predictive bioinformatics approaches that can be verified with traditional scientific rigor. Sequencing technologies that detect modification coupled with mass spectrometry to discover new adducts is a powerful tactic to study bacterial epigenetics, which is poised to make novel and far-reaching discoveries that link biological significance and the bacterial epigenome.

TstI, a Type II restriction-modification protein with DNA recognition, cleavage and methylation functions in a single polypeptide

Type II restriction–modification systems cleave and methylate DNA at specific sequences. However, the Type IIB systems look more like Type I than conventional Type II schemes as they employ the same protein for both restriction and modification and for DNA recognition. Several Type IIB proteins, including the archetype BcgI, are assemblies of two polypeptides: one with endonuclease and methyltransferase roles, another for DNA recognition. Conversely, some IIB proteins express all three functions from separate segments of a single polypeptide. This study analysed one such single-chain protein, TstI. Comparison with BcgI showed that the one- and the two-polypeptide systems differ markedly. Unlike the heterologous assembly of BcgI, TstI forms a homotetramer. The tetramer bridges two recognition sites before eventually cutting the DNA in both strands on both sides of the sites, but at each site the first double-strand break is made long before the second. In contrast, BcgI cuts all eight target bonds at two sites in a single step. TstI also differs from BcgI in either methylating or cleaving unmodified sites at similar rates. The site may thus be modified before TstI can make the second double-strand break. TstI MTase acts best at hemi-methylated sites.

Restriction enzyme body doubles and PCR cloning: on the general use of type IIs restriction enzymes for cloning

The procedure described here allows the cloning of PCR fragments containing a recognition site of the restriction endonuclease (Type IIP) used for cloning in the sequence of the insert. A Type IIS endonuclease - a Body Double of the Type IIP enzyme - is used to generate the same protruding palindrome. Thus, the insert can be cloned to the Type IIP site of the vector without digesting the PCR product with the same Type IIP enzyme. We achieve this by incorporating the recognition site of a Type IIS restriction enzyme that cleaves the DNA outside of its recognition site in the PCR primer in such a way that the cutting positions straddle the desired overhang sequence. Digestion of the PCR product by the Body Double generates the required overhang. Hitherto the use of Type IIS restriction enzymes in cloning reactions has only been used for special applications, the approach presented here makes Type IIS enzymes as useful as Type IIP enzymes for general cloning purposes. To assist in finding Body Double enzymes, we summarised the available Type IIS enzymes which are potentially useful for Body Double cloning and created an online program (http://group.szbk.u-szeged.hu/welkergr/body_double/index.html) for the selection of suitable Body Double enzymes and the design of the appropriate primers.

Type III restriction endonucleases are heterotrimeric: comprising one helicase-nuclease subunit and a dimeric methyltransferase that binds only one specific DNA

Fundamental aspects of the biochemistry of Type III restriction endonucleases remain unresolved despite being characterized by numerous research groups in the past decades. One such feature is the subunit stoichiometry of these hetero-oligomeric enzyme complexes, which has important implications for the reaction mechanism. In this study, we present a series of results obtained by native mass spectrometry and size exclusion chromatography with multi-angle light scattering consistent with a 1:2 ratio of Res to Mod subunits in the EcoP15I, EcoPI and PstII complexes as the main holoenzyme species and a 1:1 stoichiometry of specific DNA (sDNA) binding by EcoP15I and EcoPI. Our data are also consistent with a model where ATP hydrolysis activated by recognition site binding leads to release of the enzyme from the site, dissociation from the substrate via a free DNA end and cleavage of the DNA. These results are discussed critically in the light of the published literature, aiming to resolve controversies and discuss consequences in terms of the reaction mechanism.

Characterization of cleavage intermediate and star sites of RM.Tth111II

Tth111II is a thermostable Type IIGS restriction enzyme that recognizes DNA sites CAARCA (R = A or G) and cleaves downstream at N11/N9. Here, the tth111IIRM gene was cloned and expressed in E. coli, and Tth111II was purified. The purified enzyme contains internally-bound S-adenosylmethionine (SAM). When the internal SAM was removed, the endonuclease activity was stimulated by adding SAM or its analog sinefungin. The cleavage intermediate is mostly top-strand nicked DNA on a single-site plasmid. Addition of duplex oligos with a cognate site stimulates cleavage activity of the one-site substrate. Tth111II cleaves a two-site plasmid DNA with equal efficiency regardless of site orientation. We propose the top-strand nicking is carried out by a Tth111II monomer and bottom-strand cleavage is carried out by a transient dimer. Tth111II methylates cleavage product-like duplex oligos CAAACAN9, but the modification rate is estimated to be much slower than the top-strand nicking rate. We cloned and sequenced a number of Tth111II star sites which are 1-bp different from the cognate sites. A biochemical pathway is proposed for the restriction and methylation activities of Tth111II.

Cellular localization and dynamics of the Mrr type IV restriction endonuclease of Escherichia coli

In this study, we examined the intracellular whereabouts of Mrr, a cryptic type IV restriction endonuclease of Escherichia coli K12, in response to different conditions. In absence of stimuli triggering its activity, Mrr was found to be strongly associated with the nucleoid as a number of discrete foci, suggesting the presence of Mrr hotspots on the chromosome. Previously established elicitors of Mrr activity, such as exposure to high (hydrostatic) pressure (HP) or expression of the HhaII methyltransferase, both caused nucleoid condensation and an unexpected coalescence of Mrr foci. However, although the resulting Mrr/nucleoid complex was stable when triggered with HhaII, it tended to be only short-lived when elicited with HP. Moreover, HP-mediated activation of Mrr typically led to cellular blebbing, suggesting a link between chromosome and cellular integrity. Interestingly, Mrr variants could be isolated that were specifically compromised in either HhaII- or HP-dependent activation, underscoring a mechanistic difference in the way both triggers activate Mrr. In general, our results reveal that Mrr can take part in complex spatial distributions on the nucleoid and can be engaged in distinct modes of activity.