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

Single-base resolution quantitative genome methylation analysis in the model bacterium Helicobacter pylori by enzymatic methyl sequencing (EM-Seq) reveals influence of strain, growth phase, and methyl homeostasis

BACKGROUND

Bacterial epigenetics is a rapidly expanding research field. DNA methylation by diverse bacterial methyltransferases (MTases) contributes to genomic integrity and replication, and many recent studies extended MTase function also to global transcript regulation and phenotypic variation. Helicobacter pylori is currently one of those bacterial species which possess the highest number and the most variably expressed set of DNA MTases. Next-generation sequencing technologies can directly detect DNA base methylation. However, they still have limitations in their quantitative and qualitative performance, in particular for cytosine methylation.

RESULTS

As a complementing approach, we used enzymatic methyl sequencing (EM-Seq), a technology recently established that has not yet been fully evaluated for bacteria. Thereby, we assessed quantitatively, at single-base resolution, whole genome cytosine methylation for all methylated cytosine motifs in two different H. pylori strains and isogenic MTase mutants. EM-Seq reliably detected both m5C and m4C methylation. We demonstrated that three different active cytosine MTases in H. pylori provide considerably different levels of average genome-wide single-base methylation, in contrast to isogenic mutants which completely lost specific motif methylation. We found that strain identity and changed environmental conditions, such as growth phase and interference with methyl donor homeostasis, significantly influenced quantitative global and local genome-wide methylation in H. pylori at specific motifs. We also identified significantly hyper- or hypo-methylated cytosines, partially linked to overlapping MTase target motifs. Notably, we revealed differentially methylated cytosines in genome-wide coding regions under conditions of methionine depletion, which can be linked to transcript regulation.

CONCLUSIONS

This study offers new knowledge on H. pylori global and local genome-wide methylation and establishes EM-Seq for quantitative single-site resolution analyses of bacterial cytosine methylation.

Cost-effective DNA methylation profiling by FML-seq

Current methods for profiling DNA methylation require costly reagents, sequencing, and labor time. We introduce fragmentation at methylated loci and sequencing (FML-seq), a sequencing library protocol that greatly reduces all these costs. Relative to other techniques tested on the same human cell lines, FML-seq produces similar measurements of absolute and differential cytosine methylation at a fraction of the price. FML-seq enables inexpensive, high-throughput experimental designs for large-scale epigenetics research projects.

The DNA cytosine methylome revealed two methylation motifs in the upstream regions of genes related to morphological and physiological differentiation in Streptomyces coelicolor A(3)2 M145

DNA methylation is an epigenetic modification detected in both prokaryotic and eukaryotic genomic DNAs. In bacteria, the importance of 5-methylcytosine (m5C) in gene expression has been less investigated than in eukaryotic systems. Through dot-blot analysis employing m5C antibodies against chromosomal DNA, we have previously demonstrated that m5C influences the differentiation of Streptomyces coelicolor A(3)2 M145 in solid sporulating and liquid non-sporulating complex media. Here, we mapped the methylated cytosines of the M145 strain growing in the defined Maltose Glutamate (MG) liquid medium. Sequencing of the M145 genome after bisulfite treatment (BS-sequencing) evidenced 3360 methylated cytosines and the two methylation motifs, GGCmCGG and GCCmCG, in the upstream regions of 321 genes. Besides, the role of cytosine methylation was investigated using the hypo-methylating agent 5'-aza-2'-deoxycytidine (5-aza-dC) in S. coelicolor cultures, demonstrating that m5C affects both growth and antibiotic biosynthesis. Finally, quantitative reverse-transcription polymerase-chain-reaction (RT-qPCR) analysis of genes containing the methylation motifs in the upstream regions showed that 5-aza-dC treatment influenced their transcriptional levels and those of the regulatory genes for two antibiotics. To the best of our knowledge, this is the first study that reports the cytosine methylome of S. coelicolor M145, supporting the crucial role ascribed to cytosine methylation in controlling bacterial gene expression.

Extension of bacterial rDNA sequencing for simultaneous methylation detection and its application in microflora analysis

Although polymerase chain reaction (PCR) amplification and sequencing of the bacterial 16S rDNA region has numerous scientific applications, it does not provide DNA methylation information. Herein, we propose a simple extension for bisulfite sequencing to investigate 5-methylcytosine residues in the bacterial 16S rDNA region from clinical isolates or flora. Multiple displacement amplification without DNA denaturation was used to preferentially pre-amplify single-stranded bacterial DNA after bisulfite conversion. Following the pre-amplification, the 16S rDNA region was analyzed using nested bisulfite PCR and sequencing, enabling the simultaneous identification of DNA methylation status and sequence data. We used this approach (termed sm16S rDNA PCR/sequencing) to identify novel methylation sites and a methyltransferase (M. MmnI) in Morganella morganii and different methylation motifs among Enterococcus faecalis strains from small volumes of clinical specimens. Further, our analysis suggested that M. MmnI may be correlated to erythromycin resistance. Thus, sm16S rDNA PCR/sequencing is a useful extension method for analyzing the DNA methylation of 16S rDNA regions in a microflora, providing additional information not provided by conventional PCR. Given the relationship between DNA methylation status and drug resistance in bacteria, we believe this technique can be effectively applied in clinical sample testing.

Anti-restriction protein ArdA promotes clinical Klebsiella pneumoniae carbapenemase (KPC)-producing K. pneumoniae spread and its molecular mechanism

BACKGROUND

Klebsiella pneumoniae carbapenemase (KPC)-producing K. pneumoniae (KPC-KP) has spread worldwide and has become a major threat to public health. The restriction modification system provides an innate defence of bacteria against plasmids or transposons, while many different types of plasmid encoding the anti-restriction protein ArdA can specifically affect the restriction activity in bacteria.

OBJECTIVES

To detect the codistribution of ArdA and blaKPC-2 plasmids in KPC-KP and explore the molecular mechanism of ArdA promoting KPC-KP spread.

METHODS

We collected 65 clinical CRKP isolates from Ningbo, China, and 68 cases of plasmid complete sequences in GenBank to determine the prevalence of ArdA gene on the K. pneumoniae blaKPC-2 plasmid. The anti-restriction function of ArdA in promoting horizontal gene transfer (HGT) was verified by transformation, conjugation and transduction methods, and the pull-down experiment was used to investigate the molecular mechanism of ArdA protein in vitro.

RESULTS

We found that ArdA was widely distributed in KPC-KP in 100% of cases, which was detected in 0% of drug susceptible K. pneumoniae, and the plasmids containing the ArdA gene in 90% of the 30 cases randomly retrieved from the database. We also verified that ArdA has a good anti-restriction function (P < 0.05) through two aspects of HGT (transformation, transduction), and explored the non-occurrence interaction of ArdA and the hsdM subunit protein of EcoKI enzyme from the perspective of protein molecules.

CONCLUSIONS

These findings suggest that the coexistence advantage of ArdA with the blaKPC-2 plasmids may provide KPC-producing K. pneumoniae with a very efficient evasion of the restriction of type I systems, which not only favours ArdA-containing mobile genetic elements in the same species HGT between bacteria also facilitates HGT between other bacterial species.

Characterization of the Type I Restriction Modification System Broadly Conserved among Group A Streptococci

Although prokaryotic DNA methylation investigations have long focused on immunity against exogenous DNA, it has been recently recognized that DNA methylation impacts gene expression and phase variation in Streptococcus pneumoniae and Streptococcus suis. A comprehensive analysis of DNA methylation is lacking for beta-hemolytic streptococci, and thus we sought to examine DNA methylation in the major human pathogen group A Streptococcus (GAS). Using a database of 224 GAS genomes encompassing 80 emm types, we found that nearly all GAS strains encode a type I restriction modification (RM) system that lacks the hsdS′ alleles responsible for impacting gene expression in S. pneumoniae and S. suis. The GAS type I system is located on the core chromosome, while sporadically present type II orphan methyltransferases were identified on prophages. By combining single-molecule real-time (SMRT) analyses of 10 distinct emm types along with phylogenomics of 224 strains, we were able to assign 13 methylation patterns to the GAS population. Inactivation of the type I RM system, occurring either naturally through phage insertion or through laboratory-induced gene deletion, abrogated DNA methylation detectable via either SMRT or MinION sequencing. Contrary to a previous report, inactivation of the type I system did not impact transcript levels of the gene (mga) encoding the key multigene activator protein (Mga) or Mga-regulated genes. Inactivation of the type I system significantly increased plasmid transformation rates. These data delineate the breadth of the core chromosomal type I RM system in the GAS population and clarify its role in immunity rather than impacting Mga regulon expression.

IMPORTANCE The advent of whole-genome approaches capable of detecting DNA methylation has markedly expanded appreciation of the diverse roles of epigenetic modification in prokaryotic physiology. For example, recent studies have suggested that DNA methylation impacts gene expression in some streptococci. The data described herein are from the first systematic analysis of DNA methylation in a beta-hemolytic streptococcus and one of the few analyses to comprehensively characterize DNA methylation across hundreds of strains of the same bacterial species. We clarify that DNA methylation in group A Streptococcus (GAS) is primarily due to a type I restriction modification (RM) system present in the core genome and does not impact mga-regulated virulence gene expression, but does impact immunity against exogenous DNA. The identification of the DNA motifs recognized by each type I RM system may assist with optimizing methods for GAS genetic manipulation and help us understand how bacterial pathogens acquire exogenous DNA elements.

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.

Genomic Stability and Genetic Defense Systems in Dolosigranulum pigrum, a Candidate Beneficial Bacterium from the Human Microbiome

Dolosigranulum pigrum is positively associated with indicators of health in multiple epidemiological studies of human nasal microbiota. Knowledge of the basic biology of D. pigrum is a prerequisite for evaluating its potential for future therapeutic use; however, such data are very limited. To gain insight into D. pigrum's chromosomal structure, pangenome, and genomic stability, we compared the genomes of 28 D. pigrum strains that were collected across 20 years. Phylogenomic analysis showed closely related strains circulating over this period and closure of 19 genomes revealed highly conserved chromosomal synteny. Gene clusters involved in the mobilome and in defense against mobile genetic elements (MGEs) were enriched in the accessory genome versus the core genome. A systematic analysis for MGEs identified the first candidate D. pigrum prophage and insertion sequence. A systematic analysis for genetic elements that limit the spread of MGEs, including restriction modification (RM), CRISPR-Cas, and deity-named defense systems, revealed strain-level diversity in host defense systems that localized to specific genomic sites, including one RM system hot spot. Analysis of CRISPR spacers pointed to a wealth of MGEs against which D. pigrum defends itself. These results reveal a role for horizontal gene transfer and mobile genetic elements in strain diversification while highlighting that in D. pigrum this occurs within the context of a highly stable chromosomal organization protected by a variety of defense mechanisms. IMPORTANCE Dolosigranulum pigrum is a candidate beneficial bacterium with potential for future therapeutic use. This is based on its positive associations with characteristics of health in multiple studies of human nasal microbiota across the span of human life. For example, high levels of D. pigrum nasal colonization in adults predicts the absence of Staphylococcus aureus nasal colonization. Also, D. pigrum nasal colonization in young children is associated with healthy control groups in studies of middle ear infections. Our analysis of 28 genomes revealed a remarkable stability of D. pigrum strains colonizing people in the United States across a 20-year span. We subsequently identified factors that can influence this stability, including genomic stability, phage predators, the role of MGEs in strain-level variation, and defenses against MGEs. Finally, these D. pigrum strains also lacked predicted virulence factors. Overall, these findings add additional support to the potential for D. pigrum as a therapeutic bacterium.