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.

Electronic supplementary material

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

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.     

A Mimicking-of-DNA-Methylation-Patterns Pipeline for Overcoming the Restriction Barrier of Bacteria

Genetic transformation of bacteria harboring multiple Restriction-Modification (R-M) systems is often difficult using conventional methods. Here, we describe a mimicking-of-DNA-methylation-patterns (MoDMP) pipeline to address this problem in three difficult-to-transform bacterial strains. Twenty-four putative DNA methyltransferases (MTases) from these difficult-to-transform strains were cloned and expressed in an Escherichia coli strain lacking all of the known R-M systems and orphan MTases. Thirteen of these MTases exhibited DNA modification activity in Southwestern dot blot or Liquid Chromatography–Mass Spectrometry (LC–MS) assays. The active MTase genes were assembled into three operons using the Saccharomyces cerevisiae DNA assembler and were co-expressed in the E. coli strain lacking known R-M systems and orphan MTases. Thereafter, results from the dot blot and restriction enzyme digestion assays indicated that the DNA methylation patterns of the difficult-to-transform strains are mimicked in these E. coli hosts. The transformation of the Gram-positive Bacillus amyloliquefaciens TA208 and B. cereus ATCC 10987 strains with the shuttle plasmids prepared from MoDMP hosts showed increased efficiencies (up to four orders of magnitude) compared to those using the plasmids prepared from the E. coli strain lacking known R-M systems and orphan MTases or its parental strain. Additionally, the gene coding for uracil phosphoribosyltransferase (upp) was directly inactivated using non-replicative plasmids prepared from the MoDMP host in B. amyloliquefaciens TA208. Moreover, the Gram-negative chemoautotrophic Nitrobacter hamburgensis strain X14 was transformed and expressed Green Fluorescent Protein (GFP). Finally, the sequence specificities of active MTases were identified by restriction enzyme digestion, making the MoDMP system potentially useful for other strains. The effectiveness of the MoDMP pipeline in different bacterial groups suggests a universal potential. This pipeline could facilitate the functional genomics of the strains that are difficult to transform.     

Helicobacter pylori interstrain restriction-modification diversity prevents genome subversion by chromosomal DNA from competing strains

Helicobacter pylori, bacteria that colonize the human gastric mucosa, possess a large number of genes for restriction-modification (R-M) systems, and essentially, every strain possesses a unique complement of functional and partial R-M systems. Nearly half of the H.pylori strains studied possess an active type IIs R-M system, HpyII, with the recognition sequence GAAGA. Recombination between direct repeats that flank the R-M cassette allows for its deletion whereas strains lacking hpyIIRM can acquire this cassette through natural transformation. We asked whether strains lacking HpyII R-M activity can acquire an active hpyIIRM cassette [containing a 1.4 kb kanamycin resistance (aphA) marker], whether such acquisition is DNase sensitive or resistant and whether restriction barriers limit acquisition of chromosomal DNA. Our results indicate that natural transformation and conjugation-like mechanisms may contribute to the transfer of large (4.8 kb) insertions of chromosomal DNA between H.pylori strains, that inactive or partial R-M systems can be reactivated upon recombination with a functional allele, consistent with their being contingency genes, and that H.pylori R-M diversity limits acquisition of chromosomal DNA fragments of ≥1 kb.     

The Epigenomic Landscape of Prokaryotes

DNA methylation acts in concert with restriction enzymes to protect the integrity of prokaryotic genomes. Studies in a limited number of organisms suggest that methylation also contributes to prokaryotic genome regulation, but the prevalence and properties of such non-restriction-associated methylation systems remain poorly understood. Here, we used single molecule, real-time sequencing to map DNA modifications including m6A, m4C, and m5C across the genomes of 230 diverse bacterial and archaeal species. We observed DNA methylation in nearly all (93%) organisms examined, and identified a total of 834 distinct reproducibly methylated motifs. This data enabled annotation of the DNA binding specificities of 620 DNA Methyltransferases (MTases), doubling known specificities for previously hard to study Type I, IIG and III MTases, and revealing their extraordinary diversity. Strikingly, 48% of organisms harbor active Type II MTases with no apparent cognate restriction enzyme. These active ‘orphan’ MTases are present in diverse bacterial and archaeal phyla and show motif specificities and methylation patterns consistent with functions in gene regulation and DNA replication. Our results reveal the pervasive presence of DNA methylation throughout the prokaryotic kingdoms, as well as the diversity of sequence specificities and potential functions of DNA methylation systems.     

Geographic distribution of methyltransferases of Helicobacter pylori: evidence of human host population isolation and migration

Background  

Helicobacter pylori colonizes the human stomach and is associated with gastritis, peptic ulcer, and gastric cancer. This ubiquitous association between H. pylori and humans is thought to be present since the origin of modern humans. The H. pylori genome encodes for an exceptional number of restriction and modifications (R-M) systems. To evaluate if R-M systems are an adequate tool to determine the geographic distribution of H. pylori strains, we typed 221 strains from Africa, America, Asia, and Europe, and evaluated the expression of different 29 methyltransferases.     

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⁶-methyladenine (^m6A) and N⁴-methylcytosine (^m4C) 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 (^m5C) recognition motifs from the SMRT® sequencing data because this modification produces weaker signals using current methods. However, all predicted ^m6A and ^m4C 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.     

A new algorithm for cluster analysis of genomic methylation: the Helicobacter pylori case

Motivation: The genomic methylation analysis is useful to typebacteria that have a high number of expressed type II methyltransferases.Methyltransferases are usually committed to Restriction andModification (R-M) systems, in which the restriction endonucleaseimposes high pressure on the expression of the cognate methyltransferasethat hinder R-M system loss. Conventional cluster methods donot reflect this tendency. An algorithm was developed for dendrogramconstruction reflecting the propensity for conservation of R-MType II systems. Results: The new algorithm was applied to 52 Helicobacter pyloristrains from different geographical regions and compared withconventional clustering methods. The algorithm works by firstgrouping strains that share a common minimum set of R-M systemsand gradually adds strains according to the number of the R-Msystems acquired. Dendrograms revealed a cluster of Africanstrains, which suggest that R-M systems are present in H.pylorigenome since its human host migrates from Africa. Availability: The software files are available at http://www.ff.ul.pt/paginas/jvitor/Bioinformatics/MCRM_algorithm.zip Contact: filipavale{at}fe.ucp.pt Supplementary information: Supplementary data are availableat Bioinformatics online. Associate Editor: Martin Bishop     

Complementation system for Helicobacter pylori

Previously Langford et al. (2006) developed the pIR203C04 complementation system for Helicobacter pylori, which can be used to complement and restore phenotypic effects in H. pylori mutant, and furthermore they used the complementation system in vivo experiments to animals without altering the ability of strain SSI to colonize mice. In their previous study, the pIR203C04 was able to transform 26695, SSI, J99, and 43504 H. pylori strains by an electroporation method. However, in the present study using a natural transformation the pIR203C04 transformed only 26695 H. pylori but not SSI, J99, 7.13, and G27 H. pylori strains. Since the useful complementation system has a limitation of narrow selection among H. pylori strains, we redesigned the complementation system for the improvement. The same intergenic chromosomal site between hp0203 and hp0204 was utilized for the new complementation system because the insertion at the intergenic site didn’t show any polar effects and disruption of other H. pylori genes. The genome sequence analysis showed that the intergenic regions among H. pylori strains may have too low homology to each others to do a homologous recombination. Thus, in addition to the short intergenic region, the fragments of the new complementation system included 3′ conserved parts of hp0203 and hp0204 coding regions. Between the fragments there are a chloramphenicol acetyltransferase cassette and multicloning sites, resulting in pKJMSH. DNA fragment of the interest can be cloned into the multicloning sites of pKJMSH and the fragment can be integrated at the intergenic region of H. pylori chromosome by the homologous recombination. Indeed, by the natural transformation, pKJMSH was able to transform all five H. pylori strains of 26695, SSI, J99, 7.13, and G27, which are common for the investigation of molecular pathogenesis. Thus, the new pKJMSH complementation system is applicable to most H. pylori wild-type stains.     

Engineering selectivity of Cutibacterium acnes phages by epigenetic imprinting

Cutibacterium acnes (C. acnes) is a gram-positive bacterium and a member of the human skin microbiome. Despite being the most abundant skin commensal, certain members have been associated with common inflammatory disorders such as acne vulgaris. The availability of the complete genome sequences from various C. acnes clades have enabled the identification of putative methyltransferases, some of them potentially belonging to restriction-modification (R-M) systems which protect the host of invading DNA. However, little is known on whether these systems are functional in the different C. acnes strains. To investigate the activity of these putative R-M and their relevance in host protective mechanisms, we analyzed the methylome of six representative C. acnes strains by Oxford Nanopore Technologies (ONT) sequencing. We detected the presence of a 6-methyladenine modification at a defined DNA consensus sequence in strain KPA171202 and recombinant expression of this R-M system confirmed its methylation activity. Additionally, a R-M knockout mutant verified the loss of methylation properties of the strain. We studied the potential of one C. acnes bacteriophage (PAD20) in killing various C. acnes strains and linked an increase in its specificity to phage DNA methylation acquired upon infection of a methylation competent strain. We demonstrate a therapeutic application of this mechanism where phages propagated in R-M deficient strains selectively kill R-M deficient acne-prone clades while probiotic ones remain resistant to phage infection.     

Strain-specific genes of Helicobacter pylori: distribution, function and dynamics

Whole-genome clustering of the two available genome sequences of Helicobacter pylori strains 26695 and J99 allows the detection of 110 and 52 strain-specific genes, respectively. This set of strain-specific genes was compared with the sets obtained with other computational approaches of direct genome comparison as well as experimental data from microarray analysis. A considerable number of novel function assignments is possible using database-driven sequence annotation, although the function of the majority of the identified genes remains unknown. Using whole-genome clustering, it is also possible to detect species-specific genes by comparing the two H.pylori strains against the genome sequence of Campylobacter jejuni. It is interesting that the majority of strain-specific genes appear to be species specific. Finally, we introduce a novel approach to gene position analysis by employing measures from directional statistics. We show that although the two strains exhibit differences with respect to strain-specific gene distributions, this is due to the extensive genome rearrangements. If these are taken into account, a common pattern for the genome dynamics of the two Helicobacter strains emerges, suggestive of certain spatial constraints that may act as control mechanisms of gene flux.     

Cofactor Requirement of HpyAV Restriction Endonuclease

Background

Helicobacter pylori is the etiologic agent of common gastritis and a risk factor for gastric cancer. It is also one of the richest sources of Type II restriction-modification (R-M) systems in microorganisms.

Principal Findings

We have cloned, expressed and purified a new restriction endonuclease HpyAV from H. pylori strain 26695. We determined the HpyAV DNA recognition sequence and cleavage site as CCTTC 6/5. In addition, we found that HpyAV has a unique metal ion requirement: its cleavage activity is higher with transition metal ions than in Mg⁺⁺. The special metal ion requirement of HpyAV can be attributed to the presence of a HNH catalytic site similar to ColE9 nuclease instead of the canonical PD-X-D/EXK catalytic site found in many other REases. Site-directed mutagenesis was carried out to verify the catalytic residues of HpyAV. Mutation of the conserved metal-binding Asn311 and His320 to alanine eliminated cleavage activity. HpyAV variant H295A displayed approximately 1% of wt activity.

Conclusions/Significance

Some HNH-type endonucleases have unique metal ion cofactor requirement for optimal activities. Homology modeling and site-directed mutagenesis confirmed that HpyAV is a member of the HNH nuclease family. The identification of catalytic residues in HpyAV paved the way for further engineering of the metal binding site. A survey of sequenced microbial genomes uncovered 10 putative R-M systems that show high sequence similarity to the HpyAV system, suggesting lateral transfer of a prototypic HpyAV-like R-M system among these microorganisms.     

Choosing a Benchtop Sequencing Machine to Characterise Helicobacter pylori Genomes

The fully annotated genome sequence of the European strain, 26695 was first published in 1997 and, in 1999, it was directly compared to the USA isolate J99, promoting two standard laboratory isolates for Helicobacter pylori (H. pylori) research. With the genomic scaffolds available from these important genomes and the advent of benchtop high-throughput sequencing technology, a bacterial genome can now be sequenced within a few days. We sequenced and analysed strains J99 and 26695 using the benchtop-sequencing machines Ion Torrent PGM and the Illumina MiSeq Nextera and Nextera XT methodologies. Using publically available algorithms, we analysed the raw data and interrogated both genomes by mapping the data and by de novo assembly. We compared the accuracy of the coding sequence assemblies to the originally published sequences. With the Ion Torrent PGM, we found an inherently high-error rate in the raw sequence data. Using the Illumina MiSeq, we found significantly more non-covered nucleotides when using the less expensive Illumina Nextera XT compared with the Illumina Nextera library creation method. We found the most accurate de novo assemblies using the Nextera technology, however, extracting an accurate multi-locus sequence type was inconsistent compared to the Ion Torrent PGM. We found the cagPAI failed to assemble onto a single contig in all technologies but was more accurate using the Nextera. Our results indicate the Illumina MiSeq Nextera method is the most accurate for de novo whole genome sequencing of H. pylori.     

DNA Methylation Assessed by SMRT Sequencing Is Linked to Mutations in Neisseria meningitidis Isolates

The Gram-negative bacterium Neisseria meningitidis features extensive genetic variability. To present, proposed virulence genotypes are also detected in isolates from asymptomatic carriers, indicating more complex mechanisms underlying variable colonization modes of N. meningitidis.We applied the Single Molecule, Real-Time (SMRT) sequencing method from Pacific Biosciences to assess the genome-wide DNA modification profiles of two genetically related N. meningitidis strains, both of serogroup A. The resulting DNA methylomes revealed clear divergences, represented by the detection of shared and of strain-specific DNA methylation target motifs. The positional distribution of these methylated target sites within the genomic sequences displayed clear biases, which suggest a functional role of DNA methylation related to the regulation of genes.DNA methylation in N. meningitidis has a likely underestimated potential for variability, as evidenced by a careful analysis of the ORF status of a panel of confirmed and predicted DNA methyltransferase genes in an extended collection of N. meningitidis strains of serogroup A. Based on high coverage short sequence reads, we find phase variability as a major contributor to the variability in DNA methylation. Taking into account the phase variable loci, the inferred functional status of DNA methyltransferase genes matched the observed methylation profiles.Towards an elucidation of presently incompletely characterized functional consequences of DNA methylation in N. meningitidis, we reveal a prominent colocalization of methylated bases with Single Nucleotide Polymorphisms (SNPs) detected within our genomic sequence collection. As a novel observation we report increased mutability also at 6mA methylated nucleotides, complementing mutational hotspots previously described at 5mC methylated nucleotides.These findings suggest a more diverse role of DNA methylation and Restriction-Modification (RM) systems in the evolution of prokaryotic genomes.     

Haemophilus influenzae phasevarions have evolved from type III DNA restriction systems into epigenetic regulators of gene expression

Phase variably expressed (randomly switching) methyltransferases associated with type III restriction-modification (R-M) systems have been identified in a variety of pathogenic bacteria. We have previously shown that a phase variable methyltransferase (Mod) associated with a type III R-M system in Haemophilus influenzae strain Rd coordinates the random switching of expression of multiple genes, and constitutes a phase variable regulon—‘phasevarion’. We have now identified the recognition site for the Mod methyltransferase in H. influenzae strain Rd as 5′-CGAAT-3′. This is the same recognition site as the previously described HinfIII system. A survey of 59 H. influenzae strains indicated significant sequence heterogeneity in the central, variable region of the mod gene associated with target site recognition. Intra- and inter-strain transformation experiments using Mod methylated or non-methylated plasmids, and a methylation site assay demonstrated that the sequence heterogeneity seen in the region encoding target site specificity does correlate to distinct target sites. Mutations were identified within the res gene in several strains surveyed indicating that Res is not functional. These data suggest that evolution of this type III R-M system into an epigenetic mechanism for controlling gene expression has, in some strains, resulted in loss of the DNA restriction function.     

A Single Nucleotide Change Affects Fur-Dependent Regulation of sodB in H. pylori

Helicobacter pylori is a significant human pathogen that has adapted to survive the many stresses found within the gastric environment. Superoxide Dismutase (SodB) is an important factor that helps H. pylori combat oxidative stress. sodB was previously shown to be repressed by the Ferric Uptake Regulator (Fur) in the absence of iron (apo-Fur regulation) [1]. Herein, we show that apo regulation is not fully conserved among all strains of H. pylori. apo-Fur dependent changes in sodB expression are not observed under iron deplete conditions in H. pylori strains G27, HPAG1, or J99. However, Fur regulation of pfr and amiE occurs as expected. Comparative analysis of the Fur coding sequence between G27 and 26695 revealed a single amino acid difference, which was not responsible for the altered sodB regulation. Comparison of the sodB promoters from G27 and 26695 also revealed a single nucleotide difference within the predicted Fur binding site. Alteration of this nucleotide in G27 to that of 26695 restored apo-Fur dependent sodB regulation, indicating that a single base difference is at least partially responsible for the difference in sodB regulation observed among these H. pylori strains. Fur binding studies revealed that alteration of this single nucleotide in G27 increased the affinity of Fur for the sodB promoter. Additionally, the single base change in G27 enabled the sodB promoter to bind to apo-Fur with affinities similar to the 26695 sodB promoter. Taken together these data indicate that this nucleotide residue is important for direct apo-Fur binding to the sodB promoter.