Strain and lineage-level methylome heterogeneity in the multi-drug resistant pathogenic Escherichia coli ST101 clone

Escherichia coli Sequence Type (ST)101 is an emerging, multi-drug resistant lineage associated with carbapenem resistance. We recently completed a comprehensive genomics study on mobile genetic elements (MGEs) and their role in blaNDM-1 dissemination within the ST101 lineage. DNA methyltransferases (MTases) are also frequently associated with MGEs, with DNA methylation guiding numerous biological processes including genomic defence against foreign DNA and regulation of gene expression. The availability of Pacific Biosciences Single Molecule Real Time Sequencing data for seven ST101 strains enabled us to investigate the role of DNA methylation on a genome-wide scale (methylome). We defined the methylome of two complete (MS6192 and MS6193) and five draft (MS6194, MS6201, MS6203, MS6204, MS6207) ST101 genomes. Our analysis identified 14 putative MTases and eight N6-methyladenine DNA recognition sites, with one site that has not been described previously. Furthermore, we identified a Type I MTase encoded within a Transposon 7-like Transposon and show its acquisition leads to differences in the methylome between two almost identical isolates. Genomic comparisons with 13 previously published ST101 draft genomes identified variations in MTase distribution, consistent with MGE differences between genomes, highlighting the diversity of active MTases within strains of a single E. coli lineage. It is well established that MGEs can contribute to the evolution of E. coli due to their virulence and resistance gene repertoires. This study emphasises the potential for mobile genetic elements to also enable highly similar bacterial strains to rapidly acquire genome-wide functional differences via changes to the methylome.