Beyond Restriction Modification: Epigenomic Roles of DNA Methylation in Prokaryotes.

The amount of bacterial and archaeal genome sequence and methylome data has greatly increased over the last decade, enabling new insights into the functional roles of DNA methylation in these organisms. Methyltransferases (MTases), the enzymes responsible for DNA methylation, are exchanged between prokaryotes through horizontal gene transfer and can function either as part of restriction-modification systems or in apparent isolation as single (orphan) genes. The patterns of DNA methylation they confer on the host chromosome can have significant effects on gene expression, DNA replication, and other cellular processes. Some processes require very stable patterns of methylation, resulting in conservation of persistent MTases in a particular lineage. Other processes require patterns that are more dynamic yet more predictable than what is afforded by horizontal gene transfer and gene loss, resulting in phase-variable or recombination-driven MTase alleles. In this review, we discuss what is currently known about the functions of DNA methylation in prokaryotes in light of these evolutionary patterns.

A novel family of sugar-specific phosphodiesterases that remove zwitterionic modifications of GlcNAc.

The zwitterions phosphorylcholine (PC) and phosphoethanolamine (PE) are often found esterified to certain sugars in polysaccharides and glycoconjugates in a wide range of biological species. One such modification involves PC attachment to the 6-carbon of N-acetylglucosamine (GlcNAc-6-PC) in N-glycans and glycosphingolipids (GSLs) of parasitic nematodes, a modification that helps the parasite evade host immunity. Knowledge of enzymes involved in the synthesis and degradation of PC and PE modifications is limited. More detailed studies on such enzymes would contribute to a better understanding of the function of PC modifications and have potential application in the structural analysis of zwitterion-modified glycans. In this study, we used functional metagenomic screening to identify phosphodiesterases encoded in a human fecal DNA fosmid library that remove PC from GlcNAc-6-PC. A novel bacterial phosphodiesterase was identified and biochemically characterized. This enzyme (termed GlcNAc-PDase) shows remarkable substrate preference for GlcNAc-6-PC and GlcNAc-6-PE, with little or no activity on other zwitterion-modified hexoses. The identified GlcNAc-PDase protein sequence is a member of the large endonuclease/exonuclease/phosphatase superfamily where it defines a distinct subfamily of related sequences of previously unknown function, mostly from Clostridium bacteria species. Finally, we demonstrate use of GlcNAc-PDase to confirm the presence of GlcNAc-6-PC in N-glycans and GSLs of the parasitic nematode Brugia malayi in a glycoanalytical workflow.

Rapid identification of methylase specificity (RIMS-seq) jointly identifies methylated motifs and generates shotgun sequencing of bacterial genomes.

DNA methylation is widespread amongst eukaryotes and prokaryotes to modulate gene expression and confer viral resistance. 5-Methylcytosine (m5C) methylation has been described in genomes of a large fraction of bacterial species as part of restriction-modification systems, each composed of a methyltransferase and cognate restriction enzyme. Methylases are site-specific and target sequences vary across organisms. High-throughput methods, such as bisulfite-sequencing can identify m5C at base resolution but require specialized library preparations and single molecule, real-time (SMRT) sequencing usually misses m5C. Here, we present a new method called RIMS-seq (rapid identification of methylase specificity) to simultaneously sequence bacterial genomes and determine m5C methylase specificities using a simple experimental protocol that closely resembles the DNA-seq protocol for Illumina. Importantly, the resulting sequencing quality is identical to DNA-seq, enabling RIMS-seq to substitute standard sequencing of bacterial genomes. Applied to bacteria and synthetic mixed communities, RIMS-seq reveals new methylase specificities, supporting routine study of m5C methylation while sequencing new genomes.

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

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

COMBREX-DB: an experiment centered database of protein function: knowledge, predictions and knowledge gaps.

The COMBREX database (COMBREX-DB; combrex.bu.edu) is an online repository of information related to (i) experimentally determined protein function, (ii) predicted protein function, (iii) relationships among proteins of unknown function and various types of experimental data, including molecular function, protein structure, and associated phenotypes. The database was created as part of the novel COMBREX (COMputational BRidges to EXperiments) effort aimed at accelerating the rate of gene function validation. It currently holds information on ∼ 3.3 million known and predicted proteins from over 1000 completely sequenced bacterial and archaeal genomes. The database also contains a prototype recommendation system for helping users identify those proteins whose experimental determination of function would be most informative for predicting function for other proteins within protein families. The emphasis on documenting experimental evidence for function predictions, and the prioritization of uncharacterized proteins for experimental testing distinguish COMBREX from other publicly available microbial genomics resources. This article describes updates to COMBREX-DB since an initial description in the 2011 NAR Database Issue.

The complex methylome of the human gastric pathogen Helicobacter pylori.

The genome of Helicobacter pylori is remarkable for its large number of restriction-modification (R-M) systems, and strain-specific diversity in R-M systems has been suggested to limit natural transformation, the major driving force of genetic diversification in H. pylori. We have determined the comprehensive methylomes of two H. pylori strains at single base resolution, using Single Molecule Real-Time (SMRT®) sequencing. For strains 26695 and J99-R3, 17 and 22 methylated sequence motifs were identified, respectively. For most motifs, almost all sites occurring in the genome were detected as methylated. Twelve novel methylation patterns corresponding to nine recognition sequences were detected (26695, 3; J99-R3, 6). Functional inactivation, correction of frameshifts as well as cloning and expression of candidate methyltransferases (MTases) permitted not only the functional characterization of multiple, yet undescribed, MTases, but also revealed novel features of both Type I and Type II R-M systems, including frameshift-mediated changes of sequence specificity and the interaction of one MTase with two alternative specificity subunits resulting in different methylation patterns. The methylomes of these well-characterized H. pylori strains will provide a valuable resource for future studies investigating the role of H. pylori R-M systems in limiting transformation as well as in gene regulation and host interaction.

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.

The methylomes of six bacteria.

Six bacterial genomes, Geobacter metallireducens GS-15, Chromohalobacter salexigens, Vibrio breoganii 1C-10, Bacillus cereus ATCC 10987, Campylobacter jejuni subsp. jejuni 81-176 and C. jejuni NCTC 11168, all of which had previously been sequenced using other platforms were re-sequenced using single-molecule, real-time (SMRT) sequencing specifically to analyze their methylomes. In every case a number of new N(6)-methyladenine ((m6)A) and N(4)-methylcytosine ((m4)C) methylation patterns were discovered and the DNA methyltransferases (MTases) responsible for those methylation patterns were assigned. In 15 cases, it was possible to match MTase genes with MTase recognition sequences without further sub-cloning. Two Type I restriction systems required sub-cloning to differentiate their recognition sequences, while four MTase genes that were not expressed in the native organism were sub-cloned to test for viability and recognition sequences. Two of these proved active. No attempt was made to detect 5-methylcytosine ((m5)C) recognition motifs from the SMRT® sequencing data because this modification produces weaker signals using current methods. However, all predicted (m6)A and (m4)C MTases were detected unambiguously. This study shows that the addition of SMRT sequencing to traditional sequencing approaches gives a wealth of useful functional information about a genome showing not only which MTase genes are active but also revealing their recognition sequences.

COMBREX: a project to accelerate the functional annotation of prokaryotic genomes.

COMBREX (http://combrex.bu.edu) is a project to increase the speed of the functional annotation of new bacterial and archaeal genomes. It consists of a database of functional predictions produced by computational biologists and a mechanism for experimental biochemists to bid for the validation of those predictions. Small grants are available to support successful bids.

A Survey of Archaeal Restriction-Modification Systems.

When compared with bacteria, relatively little is known about the restriction-modification (RM) systems of archaea, particularly those in taxa outside of the haloarchaea. To improve our understanding of archaeal RM systems, we surveyed REBASE, the restriction enzyme database, to catalog what is known about the genes and activities present in the 519 completely sequenced archaeal genomes currently deposited there. For 49 (9.4%) of these genomes, we also have methylome data from Single-Molecule Real-Time (SMRT) sequencing that reveal the target recognition sites of the active m6A and m4C DNA methyltransferases (MTases). The gene-finding pipeline employed by REBASE is trained primarily on bacterial examples and so will look for similar genes in archaea. Nonetheless, the organizational structure and protein sequence of RM systems from archaea are highly similar to those of bacteria, with both groups acquiring systems from a shared genetic pool through horizontal gene transfer. As in bacteria, we observe numerous examples of "persistent" DNA MTases conserved within archaeal taxa at different levels. We experimentally validated two homologous members of one of the largest "persistent" MTase groups, revealing that methylation of C(m5C)WGG sites may play a key epigenetic role in Crenarchaea. Throughout the archaea, genes encoding m6A, m4C, and m5C DNA MTases, respectively, occur in approximately the ratio 4:2:1.

Genomic Analysis of Haloarchaea from Diverse Environments, including Permian Halite, Reveals Diversity of Ultraviolet Radiation Survival and DNA Photolyase Gene Variants.

Ultraviolet (UV) radiation responses of extremophilic and archaeal microorganisms are of interest from evolutionary, physiological, and astrobiological perspectives. Previous studies determined that the halophilic archaeon, Halobacterium sp. NRC-1, which survives in multiple extremes, is highly tolerant of UV radiation. Here, Halobacterium sp. NRC-1 UV tolerance was compared to taxonomically diverse Haloarchaea isolated from high-elevation salt flats, surface warm and cold hypersaline lakes, and subsurface Permian halite deposits. Haloterrigena/Natrinema spp. from subsurface halite deposits were the least tolerant after exposure to photoreactivating light. This finding was attributed to deviation of amino acid residues in key positions in the DNA photolyase enzyme or to the complete absence of the photolyase gene. Several Halobacterium, Halorubrum and Salarchaeum species from surface environments exposed to high solar irradiance were found to be the most UV tolerant, and Halorubrum lacusprofundi from lake sediment was of intermediate character. These results indicate that high UV tolerance is not a uniform character trait of Haloarchaea and is likely reflective of UV exposure experienced in their environment. This is the first report correlating natural UV tolerance to photolyase gene functionality among Haloarchaea and provides insights into their survival in ancient halite deposits and potentially on the surface of Mars.

Functional characterization of the YmcB and YqeV tRNA methylthiotransferases of Bacillus subtilis.

Methylthiotransferases (MTTases) are a closely related family of proteins that perform both radical-S-adenosylmethionine (SAM) mediated sulfur insertion and SAM-dependent methylation to modify nucleic acid or protein targets with a methyl thioether group (-SCH(3)). Members of two of the four known subgroups of MTTases have been characterized, typified by MiaB, which modifies N(6)-isopentenyladenosine (i(6)A) to 2-methylthio-N(6)-isopentenyladenosine (ms(2)i(6)A) in tRNA, and RimO, which modifies a specific aspartate residue in ribosomal protein S12. In this work, we have characterized the two MTTases encoded by Bacillus subtilis 168 and find that, consistent with bioinformatic predictions, ymcB is required for ms(2)i(6)A formation (MiaB activity), and yqeV is required for modification of N(6)-threonylcarbamoyladenosine (t(6)A) to 2-methylthio-N(6)-threonylcarbamoyladenosine (ms(2)t(6)A) in tRNA. The enzyme responsible for the latter activity belongs to a third MTTase subgroup, no member of which has previously been characterized. We performed domain-swapping experiments between YmcB and YqeV to narrow down the protein domain(s) responsible for distinguishing i(6)A from t(6)A and found that the C-terminal TRAM domain, putatively involved with RNA binding, is likely not involved with this discrimination. Finally, we performed a computational analysis to identify candidate residues outside the TRAM domain that may be involved with substrate recognition. These residues represent interesting targets for further analysis.

Genome Sequence of the Early 20th-Century Extreme Halophile Halobacterium sp. Strain NRC-34001.

Halobacterium sp. strain NRC-34001 is a red, extremely halophilic archaeon isolated in Canada in 1934. Single-molecule real-time sequencing revealed a 2.3-Mbp genome with a 2-Mbp chromosome and two plasmids (235 kb and 43 kb). The genome encodes all conserved core haloarchaeal groups (cHOGs) and a highly acidic proteome.

RimO, a MiaB-like enzyme, methylthiolates the universally conserved Asp88 residue of ribosomal protein S12 in Escherichia coli.

Ribosomal protein S12 undergoes a unique posttranslational modification, methylthiolation of residue D88, in Escherichia coli and several other bacteria. Using mass spectrometry, we have identified the enzyme responsible for this modification in E. coli, the yliG gene product. This enzyme, which we propose be called RimO, is a radical-S-adenosylmethionine protein that bears strong sequence similarity to MiaB, which methylthiolates tRNA. We show that RimO and MiaB represent two of four subgroups of a larger, ancient family of likely methylthiotransferases, the other two of which are typified by Bacillus subtilis YqeV and Methanococcus jannaschii Mj0867, and we predict that RimO is unique among these subgroups in its modification of protein as opposed to tRNA. Despite this, RimO has not significantly diverged from the other three subgroups at the sequence level even within the C-terminal TRAM domain, which in the methyltransferase RumA is known to bind the RNA substrate and which we presume to be responsible for substrate binding and recognition in all four subgroups of methylthiotransferases. To our knowledge, RimO and MiaB represent the most extreme known case of resemblance between enzymes modifying protein and nucleic acid. The initial results presented here constitute a bioinformatics-driven prediction with preliminary experimental validation that should serve as the starting point for several interesting lines of further inquiry.

Genome-wide identification of 5-methylcytosine sites in bacterial genomes by high-throughput sequencing of MspJI restriction fragments.

Single-molecule Real-Time (SMRT) sequencing can easily identify sites of N6-methyladenine and N4-methylcytosine within DNA sequences, but similar identification of 5-methylcytosine sites is not as straightforward. In prokaryotic DNA, methylation typically occurs within specific sequence contexts, or motifs, that are a property of the methyltransferases that "write" these epigenetic marks. We present here a straightforward, cost-effective alternative to both SMRT and bisulfite sequencing for the determination of prokaryotic 5-methylcytosine methylation motifs. The method, called MFRE-Seq, relies on excision and isolation of fully methylated fragments of predictable size using MspJI-Family Restriction Enzymes (MFREs), which depend on the presence of 5-methylcytosine for cleavage. We demonstrate that MFRE-Seq is compatible with both Illumina and Ion Torrent sequencing platforms and requires only a digestion step and simple column purification of size-selected digest fragments prior to standard library preparation procedures. We applied MFRE-Seq to numerous bacterial and archaeal genomic DNA preparations and successfully confirmed known motifs and identified novel ones. This method should be a useful complement to existing methodologies for studying prokaryotic methylomes and characterizing the contributing methyltransferases.

Complete Genome Sequence of an Extremely Halophilic Archaeon from Great Salt Lake, Halobacterium sp. GSL-19.

An extremely halophilic archaeon, Halobacterium sp. GSL-19, was isolated from the north arm of Great Salt Lake in Utah. Single-molecule real-time (SMRT) sequencing was used to establish a GC-rich 2.3-Mbp genome composed of a circular chromosome and 2 plasmids, with 2,367 predicted genes, including 1 encoding a CTAG-methylase widely distributed among Haloarchaea.

Genome Sequence of Halobacterium sp. Strain BOL4-2, Isolated and Cultured from Salar de Uyuni, Bolivia.

Halobacterium sp. strain BOL4-2 was isolated from an Andean salt flat, Salar de Uyuni, in Bolivia. Single-molecule real-time (SMRT) sequencing revealed a 2.4-Mbp genome with a 2.0-Mbp chromosome and four plasmids (2 to 299 kb). Its isolation from an environment experiencing multiple extremes makes the strain interesting for astrobiology.

Genome-wide methylome analysis of two strains belonging to the hypervirulent Neisseria meningitidis serogroup W ST-11 clonal complex.

A rising incidence of meningococcal serogroup W disease has been evident in many countries worldwide. Serogroup W isolates belonging to the sequence type (ST)-11 clonal complex have been associated with atypical symptoms and increased case fatality rates. The continued expansion of this clonal complex in the later part of the 2010s has been largely due to a shift from the so-called original UK strain to the 2013 strain. Here we used single-molecule real-time (SMRT) sequencing to determine the methylomes of the two major serogroup W strains belonging to ST-11 clonal complex. Five methylated motifs were identified in this study, and three of the motifs, namely 5'-GATC-3', 5'-GAAGG-3', 5'-GCGCGC-3', were found in all 13 isolates investigated. The results showed no strain-specific motifs or difference in active restriction modification systems between the two strains. Two phase variable methylases were identified and the enrichment or depletion of the methylation motifs generated by these methylases varied between the two strains. Results from this work give further insight into the low diversity of methylomes in highly related strains and encourage further research to decipher the role of regions with under- or overrepresented methylation motifs.

Genome Sequence and Methylation Pattern of Haloterrigena salifodinae BOL5-1, an Extremely Halophilic Archaeon from a Bolivian Salt Mine.

The halophilic archaeon Haloterrigena salifodinae BOL5-1 was isolated from a Bolivian salt mine and sequenced using single-molecule real-time sequencing. The GC-rich genome was 5.1 Mbp, with a 4.2-Mbp chromosome and 5 plasmids ranging from 96 to 281 kbp. The genome annotation was incorporated into HaloWeb (https://halo.umbc.edu), and the methylation patterns were incorporated into REBASE (http://tools.neb.com/genomes/view.php?seq_id=99167&list=1).