The Epigenomic Landscape of Prokaryotes

Mapping-based genome size estimation

While the size of chromosomes can be measured under a microscope, obtaining the exact size of a genome remains a challenge. Biochemical methods and k-mer distribution-based approaches allow only estimations. An alternative approach to estimate the genome size based on high contiguity assemblies and read mappings is presented here. Analyses of Arabidopsis thaliana and Beta vulgaris data sets are presented to show the impact of different parameters. Oryza sativa, Brachypodium distachyon, Solanum lycopersicum, Vitis vinifera, and Zea mays were also analyzed to demonstrate the broad applicability of this approach. Further, MGSE was also used to analyze Escherichia coli, Saccharomyces cerevisiae, and Caenorhabditis elegans datasets to show its utility beyond plants. Mapping-based Genome Size Estimation (MGSE) and additional scripts are available on GitHub: https://github.com/bpucker/MGSE . MGSE predicts genome sizes based on short reads or long reads requiring a minimal coverage of 5-fold.

Comprehensive comparison of the third-generation sequencing tools for bacterial 6mA profiling

DNA N6-methyladenine (6mA) serves as an intrinsic and principal epigenetic marker in prokaryotes, impacting various biological processes. To date, limited advanced sequencing technologies and analyzing tools are available for bacterial DNA 6mA. Here, we evaluate eight tools designed for the 6mA identification or de novo methylation detection. This assessment includes Nanopore (R9 and R10), Single-Molecule Real-Time (SMRT) Sequencing, and cross-reference with 6mA-IP-seq and DR-6mA-seq. Our multi-dimensional evaluation report encompasses motif discovery, site-level accuracy, single-molecule accuracy, and outlier detection across six bacteria strains. While most tools correctly identify motifs, their performance varies at single-base resolution, with SMRT and Dorado consistently delivering strong performance. Our study indicates that existing tools cannot accurately detect low-abundance methylation sites. Additionally, we introduce an optimized method for advancing 6mA prediction, which substantially improves the detection performance of Dorado. Overall, our study provides a robust and detailed examination of computational tools for bacterial 6mA profiling, highlighting insights for further tool enhancement and epigenetic research.

Decoding bacterial methylomes in four public health-relevant microbial species: nanopore sequencing enables reproducible analysis of DNA modifications

Investigating bacterial methylation profiles provides essential complementary information to the native DNA sequence, significantly extending our understanding of how DNA modifications influence virulence, antibiotic resistance, and the ability of bacteria to evade the immune system. Recent advancements in real-time Nanopore sequencing and basecalling algorithms have enabled the direct detection of modified bases from raw signal data, eliminating the need for bisulfite treatment of DNA. However, decoding methylation signals remains challenging due to rapid technological and methodological progress. In this study, we focus on public health-relevant bacterial strains to analyze their methylation profiles and identify methylation motifs. Our dataset includes samples from Staphylococcus aureus, Listeria monocytogenes, Enterococcus faecium, and Klebsiella pneumoniae, sequenced on the Nanopore GridION platform using the latest flow cell chemistry (R10.4.1) and modification basecalling models (Dorado basecalling SUP model v5). We investigated distinct methylation patterns within and between species, focusing on heavily modified genes or genomic regions. Our results reveal distinct species-specific methylation profiles, with each strain exhibiting unique modification patterns. We developed a modular pipeline using Nextflow and the Nanopore Modkit tool to streamline the detection of methylated motifs. We compared the results with outputs from MicrobeMod, a recent toolkit for exploring prokaryotic methylation and base modifications in nanopore sequencing. Our pipeline is publicly available for further use (github.com/rki-mf1/ont-methylation). We identified known methylation motifs already described in the literature and novel de novo motifs, providing deeper insights into the diversity of bacterial DNA modifications. Furthermore, we identified genomic regions that are extensively methylated, which could have implications for bacterial behavior and pathogenicity. We also assess improvements in basecalling accuracy, specifically how methylated bases can influence neighboring basecalls. Recent advances in basecalling models, particularly v5 models as part of Dorado, have reduced these issues, improving the reliability of methylation detection in bacterial genomes. In conclusion, our study highlights the potential of current nanopore sequencing tools for detecting DNA modifications in prokaryotes. By making our pipeline and results publicly available, we facilitate further research into bacterial DNA modifications and their role in microbial pathogenesis.

Evolution of a Plasmid Regulatory Circuit Ameliorates Plasmid Fitness Cost

Plasmids promote adaptation of bacteria by facilitating horizontal transfer of diverse genes, notably those conferring antibiotic resistance. Some plasmids, like those of the incompatibility group IncP-1, are known to replicate and persist in a broad range of bacteria. We investigated a poorly understood exception, the IncP-1β plasmid pBP136 from a clinical Bordetella pertussis isolate, which quickly became extinct in laboratory Escherichia coli populations. Through experimental evolution, we found that the inactivation of a previously uncharacterized plasmid gene, upf31, drastically improved plasmid persistence in E. coli. The gene inactivation caused alterations in the plasmid regulatory system, including decreased transcription of the global plasmid regulators (korA, korB, and korC) and numerous genes in their regulons. This is consistent with our findings that Upf31 represses its own transcription. It also caused secondary transcriptional changes in many chromosomal genes. In silico analyses predicted that Upf31 interacts with the plasmid regulator KorB at its C-terminal dimerization domain (CTD). We showed experimentally that adding the CTD of upf31/pBP136 to the naturally truncated upf31 allele of the stable IncP-1β archetype R751 results in plasmid destabilization in E. coli. Moreover, mutagenesis showed that upf31 alleles encoded on nearly half of the sequenced IncP-1β plasmids also possess this destabilization phenotype. While Upf31 might be beneficial in many hosts, we show that in E. coli some alleles have harmful effects that can be rapidly alleviated with a single mutation. Thus, broad-host-range plasmid adaptation to new hosts can involve fine-tuning their transcriptional circuitry through evolutionary changes in a single gene.

Decoding Adenine DNA Methylation Effects in Streptococcus Mutans: Insights Into Self-DNA Protection and Autoaggregation

Streptococcus mutans, a key player in dental caries, faces multiple environmental challenges within the oral cavity, including oxidative stress, nutrient scarcity, and acidic pH. To survive and thrive, S. mutans has evolved intricate mechanisms, including the CSP-ComDE quorum sensing system, which coordinates responses to environmental cues. The CSP-ComDE system enables S. mutans to communicate with neighboring cells via its CSP pheromone. Under stress conditions, the CSP pheromone production increases, triggering a cascade of events. Notably, our research demonstrated that the CSP pheromone activates the expression of a Type II restriction-modification (R-M) system. Type II R-M systems are well-known tools in molecular biology and genetic engineering and consist of two distinct enzymes: a restriction enzyme and a methyltransferase. An increasing number of studies have revealed that bacterial adenine methylation (Dam methylation) has a broader role beyond mere DNA protection. In fact, the marks introduced into the DNA provide signals for a variety of physiological processes. Our results highlight a conserved chromosomal locus in S. mutans encoding the DpnII R-M system. DpnII R-M methylates DNA at 5'-GATC target sites within the S. mutans genome and cleaves unmarked DNA. Furthermore, our findings suggest that Dam methylation significantly impacts foreign DNA acquisition via natural transformation and modulates mutanobactin expression-a secondary metabolite linked to oxidative stress tolerance. Collectively, our findings suggest that Dam methylation bridges epigenetics and bacterial fitness, potentially opening new avenues in bacterial epigenetics. As we explore this intricate biological process, we may uncover novel therapeutic strategies to combat bacterial infections.

Epigenetics in evolution and adaptation to environmental challenges: pathways for disease prevention and treatment

Adaptation to challenging environmental conditions is crucial for the survival/fitness of all organisms. Alongside genetic mutations that provide adaptive potential during environmental challenges, epigenetic modifications offer dynamic, reversible, and rapid mechanisms for regulating gene expression in response to environmental changes in both evolution and daily life, without altering DNA sequences or relying on accidental favorable mutations. The widespread conservation of diverse epigenetic mechanisms - like DNA methylation, histone modifications, and RNA interference across diverse species, including plants - underscores their significance in evolutionary biology. Remarkably, environmentally induced epigenetic alterations are passed to daughter cells and inherited transgenerationally through germline cells, shaping offspring phenotypes while preserving adaptive epigenetic memory. Throughout anthropoid evolution, epigenetic modifications have played crucial roles in: i) suppressing transposable elements and viral genomes intruding into the host genome; ii) inactivating one of the X chromosomes in female cells to balance gene dosage; iii) genetic imprinting to ensure expression from one parental allele; iv) regulating functional alleles to compensate for dysfunctional ones; and v) modulating the epigenome and transcriptome in response to influence from the gut microbiome among other functions. Understanding the interplay between environmental factors and epigenetic processes may provide valuable insights into developmental plasticity, evolutionary dynamics, and disease susceptibility.

Decoding the genome and epigenome of avian Escherichia coli strains by R10.4.1 nanopore sequencing

Avian pathogenic Escherichia coli (APEC) causes colibacillosis in poultry, which is a very important disease worldwide. Despite well-documented genomic traits and diversity of APEC, its epigenomic characteristics are less understood. This study utilized the high throughput and long-read capabilities of Oxford Nanopore Technology (ONT) to elucidate the genome structures and methylation modifications of three E. coli isolates of avian origin: one intestinal isolate from a healthy wild bird and two systemic isolates from clinically affected chickens. Three complete genomes, each comprising a single chromosome and multiple plasmids were assembled. Diverse virulence-associated genes, antimicrobial resistance genes, mobile genetic elements plasmids and integrons were characterized from the genomes. Despite a limited sample size, our whole genome sequencing (WGS) data highlighted significant genomic diversity among the E. coli strains and enriched repertoire of gene clusters related to APEC pathogenicity. From the epigenetic analysis, multiple methylation modifications, including three N5-methylcytosine (5mC), eight N6-methyladenine (6mA) and two N4-methylcytosine (4mC) modification motifs were identified within all three isolates. Furthermore, common GATC and CCWGG methylation motifs were predominantly distributed within regulatory regions, suggesting a role in epigenetic transcription regulation. This study opens the avenue for future research into pathogenesis, diagnostic and therapeutic strategies of APEC considering epigenetic analysis.

Plasmid-borne mcr-1 and replicative transposition of episomal and chromosomal blaNDM-1, blaOXA-69, and blaOXA-23 carbapenemases in a clinical Acinetobacter baumannii isolate

A multidrug-resistant clinical Acinetobacter baumannii isolate with resistance to most antibiotics was isolated from a patient at an intensive care unit. The genetic environment, transcriptome, mobile, and resistome were characterized. The MicroScan system, disc diffusion, and broth microdilution were used to determine the resistance profile of the isolate. A multiplex PCR assay was also used to screen for carbapenemases and mcr-1 to -5 resistance genes. Efflux-pump inhibitors were used to evaluate efflux activity. The resistome, mobilome, epigenome, and transcriptome were characterized. There was phenotypic resistance to 22 of the 25 antibiotics tested, intermediate resistance to levofloxacin and nalidixic acid, and susceptibility to tigecycline, which corresponded to the 27 resistance genes found in the genome, most of which occurred in multiple copies through replicative transposition. A plasmid-borne (pR-B2.MM_C3) mcr-1 and chromosomal blaPER-7, blaOXA-69, blaOXA-23 (three copies), blaADC-25, blaTEM-1B, and blaNDM-1 were found within composite transposons, ISs, and/or class 1 and 2 integrons on genomic islands. Types I and II methylases and restriction endonucleases were in close synteny to these resistance genes within the genomic islands; chromosomal genomic islands aligned with known plasmids. There was a closer evolutionary relationship between the strain and global strains but not local or regional strains; the resistomes also differed. Significantly expressed/repressed genes (6.2%) included resistance genes, hypothetical proteins, mobile elements, methyltransferases, transcription factors, and membrane and efflux proteins. The genomic evolution observed in this strain explains its adaptability and pandrug resistance and shows its genomic plasticity on exposure to antibiotics.

IMPORTANCE

A pandrug-resistant pathogen that was susceptible to only a single antibiotic, tigecycline, was isolated from a middle-aged patient in an ICU. This pathogen contained two plasmids and had a chromosome that contained portions that were integrated externally from plasmids. These genomic islands were rich with resistance genes, mobile genetic elements, and restriction-modification systems that protected the pathogen and facilitated gene regulation. The strain contained 35 resistance genes and 12 virulence genes. The strain was of closer evolutionary distance to several international strains suggesting that it was imported into South Africa. However, its resistome was unique, suggesting an independent evolution on exposure to antibiotic therapy mediated by epigenomic factors and MGE transposition events. The varied mechanisms available to this strain to overcome antibiotic resistance and spread to other areas and/or transfer its resistance determinants are worrying. This is ultimately a risk to public health, evincing the need for antibiotic stewardship.

DNA methylome regulates virulence and metabolism in Pseudomonas syringae

Bacterial pathogens employ epigenetic mechanisms, including DNA methylation, to adapt to environmental changes, and these mechanisms play important roles in various biological processes. Pseudomonas syringae is a model phytopathogenic bacterium, but its methylome is less well known than that of other species. In this study, we conducted single-molecule real-time sequencing to profile the DNA methylation landscape in three model pathovars of P. syringae. We identified one Type I restriction-modification system (HsdMSR), including the conserved sequence motif associated with N6-methyladenine (6mA). About 25-40% of the genes involved in DNA methylation were conserved in two or more of the strains, revealing the functional conservation of methylation in P. syringae. Subsequent transcriptomic analysis highlighted the involvement of HsdMSR in virulent and metabolic pathways, including the Type III secretion system, biofilm formation, and translational efficiency. The regulatory effect of HsdMSR on transcription was dependent on both strands being fully 6mA methylated. Overall, this work illustrated the methylation profile in P. syringae and the critical involvement of DNA methylation in regulating virulence and metabolism. Thus, this work contributes to a deeper understanding of epigenetic transcriptional control in P. syringae and related bacteria.

Streamlining whole genome sequencing for clinical diagnostics with ONT technology

Recent advances in whole-genome sequencing (WGS) have increased the accessibility of this tool, offering substantial potential for pathogen surveillance, outbreak response, and diagnostics. However, the routine clinical adoption of WGS is hindered by factors such as high costs, technical complexity, and the requirement for bioinformatics expertise for data analysis. To address these challenges, we propose RapidONT, a workflow designed for cost-effective and accessible WGS-based pathogen analysis. RapidONT employs a mechanical shearing-based DNA extraction protocol, followed by library construction by using a multiplexing Oxford nanopore technologies (ONT) rapid barcoding kit. Flye software is used for de novo assembly without manual intervention, followed by basic assembly polishing using Medaka and Homopolish. The polished assemblies are then analyzed using the user-friendly web-based platform Pathogenwatch, which facilitates species identification, molecular typing, and antimicrobial resistance (AMR) prediction, all while requiring minimal bioinformatics expertise. The efficacy of RapidONT was evaluated using nine clinically relevant pathogens, encompassing a total of 90 gram-positive and gram-negative bacterial strains. The workflow demonstrated high accuracy in critical tasks such as multilocus sequence typing (MLST) and AMR identification, using only ONT R9.4.1 flowcell data. Notably, limitations were observed with Salmonella spp. and Neisseria gonorrhoeae. Furthermore, RapidONT enabled the generation of genomic information for 48 bacterial isolates by using a single flow cell, significantly reducing sequencing costs. This approach eliminates the need for extensive experimentation in obtaining crucial genomic information. This workflow facilitates broader WGS implementation in clinical pathogen analysis and diagnostics.

Methylation of foreign DNA overcomes the restriction barrier of Flavobacterium psychrophilum and allows efficient genetic manipulation

Flavobacterium psychrophilum causes bacterial cold-water disease (BCWD) in salmonids and other fish, resulting in substantial economic losses in aquaculture worldwide. The mechanisms F. psychrophilum uses to cause disease are poorly understood. Despite considerable effort, most strains of F. psychrophilum have resisted attempts at genetic manipulation. F. psychrophilum restriction-modification (R-M) systems may contribute to this resistance. Restriction endonucleases (REases) rapidly degrade nonself DNA if it is not properly methylated by their cognate DNA methyltransferases (MTases). We used comparative genomics to show that R-M systems are abundant in F. psychrophilum and that strain-specific variations partially align with phylogeny. We identified two critical type II R-M systems, HpaII-like (FpsJI) and ScrFI-like (FpsJVI), that are conserved in most of the sequenced strains. Protection of foreign DNA against HpaII and ScrFI was achieved by expression of the MTases M.FpsJI and M.FpsJVI in Escherichia coli. Furthermore, deleting the two REase genes from F. psychrophilum resulted in efficient conjugative DNA transfer from E. coli into the otherwise genetically intractable F. psychrophilum strain CSF259-93. This allowed us to construct a CSF259-93 mutant lacking gldN, a core component of the type IX protein secretion system. The pre-methylation system developed in this study functions in all tested F. psychrophilum strains harboring HpaII-like and ScrFI-like REases. These newly developed genetic tools may allow the identification of key virulence factors and facilitate the development of live attenuated vaccines or other measures to control BCWD.

IMPORTANCE

Bacterial cold-water disease (BCWD) caused by Flavobacterium psychrophilum is a problem for salmonid aquaculture worldwide, and current control measures are inadequate. An obstacle in understanding and controlling BCWD is that most F. psychrophilum strains resist DNA transfer, thus limiting genetic studies of their virulence mechanisms. F. psychrophilum restriction enzymes that destroy foreign DNA were suspected to contribute to this problem. Here, we used F. psychrophilum DNA methyltransferases to modify and protect foreign DNA from digestion. This allowed efficient conjugative DNA transfer into nine diverse F. psychrophilum strains that had previously resisted DNA transfer. Using this approach, we constructed a gene deletion mutant that failed to cause disease in rainbow trout. Further genetic studies could help determine the molecular factors involved in pathogenesis and may aid development of innovative BCWD control strategies.

Metagenomic global survey and in-depth genomic analyses of Ruminococcus gnavus reveal differences across host lifestyle and health status

Ruminococcus gnavus is a gut bacterium found in > 90% of healthy individuals, but its increased abundance is also associated with chronic inflammatory diseases, particularly Crohn's disease. Nevertheless, its global distribution and intraspecies genomic variation remain understudied. By surveying 12,791 gut metagenomes, we recapitulated known associations with metabolic diseases and inflammatory bowel disease. We uncovered a higher prevalence and abundance of R. gnavus in Westernized populations and observed bacterial relative abundances up to 83% in newborns. Next, we built a resource of R. gnavus isolates (N = 45) from healthy individuals and Crohn's disease patients and generated complete R. gnavus genomes using PacBio circular consensus sequencing. Analysis of these genomes and publicly available high-quality draft genomes (N = 333 genomes) revealed multiple clades which separated Crohn's-derived isolates from healthy-derived isolates. Presumed R. gnavus virulence factors could not explain this separation. Bacterial genome-wide association study revealed that Crohn's-derived isolates were enriched in genes related to mobile elements and mucin foraging. Together, we present a large R. gnavus resource that will be available to the scientific community and provide novel biological insights into the global distribution and genomic variation of R. gnavus.

The Complex Epigenetic Panorama in the Multipartite Genome of the Nitrogen-Fixing Bacterium Sinorhizobium meliloti

In prokaryotes, DNA methylation plays roles in DNA repair, gene expression, cell cycle progression, and immune recognition of foreign DNA. Genome-wide methylation patterns can vary between strains, influencing phenotype, and gene transfer. However, broader evolutionary studies on bacterial epigenomic variation remain limited. In this study, we conducted an epigenomic analysis using single-molecule real-time sequencing on 21 strains of Sinorhizobium meliloti, a facultative plant nitrogen-fixing alphaproteobacterium. This species is notable for its multipartite genome structure, consisting of a chromosome, chromid, and megaplasmid, leading to significant genomic and phenotypic diversity. We identified 16 palindromic and nonpalindromic methylated DNA motifs, including N4-methylcytosine and N6-methyladenine modifications, and analyzed their associated methyltransferases. Some motifs were methylated across all strains, forming a core set of epigenomic signatures, while others exhibited variable methylation frequencies, indicating a dispensable (shell) epigenome. Additionally, we observed differences in methylation frequency between replicons and within coding sequences versus regulatory regions, suggesting that methylation patterns may reflect multipartite genome evolution and influence gene regulation. Overall, our findings reveal extensive epigenomic diversity in S. meliloti, with complex epigenomic signatures varying across replicons and genomic regions. These results enhance our understanding of multipartite genome evolution and highlight the potential role of epigenomic diversity in phenotypic variation.

Methionine-driven methylation modification overcomes plasmid-mediated high-level tigecycline resistance

Tigecycline is a last-resort antibiotic to treat complicated infections caused by multidrug-resistant pathogens, while the emergence of plasmid-mediated tet(X) family severely compromises its clinical efficacy. Novel antimicrobial strategies not limited to new antibiotics in pharmaceutical pipeline are urgently needed. Herein, we reveal the metabolic disparities between tet(X)-negative and -positive E. coli, including distinct energy demand patterns under tigecycline exposure. In particular, the cysteine and methionine metabolism pathway is remarkably downregulated in tet(X)-positive bacteria. More importantly, we find that the addition of exogenous L-methionine (Met) effectively resensitizes tet(X)-positive pathogens to tigecycline. Our mechanistic analysis demonstrates that exogenous Met promotes intracellular tigecycline accumulation by upregulating bacterial proton motive force. Moreover, Met accelerates the conversion to S-adenosyl-L-methionine, an essential methyl donor, thereby enhancing 5mC methylation modification in the promoter region of tet(X4) gene and reducing its expression. Consistently, the potentiation of Met to tigecycline is abolished in tet(X4)-carrying E. coli Δdcm but restored in dcm-complementary bacteria, which encodes DNA-cytosine methyltransferase. In multiple animal models of infection, Met markedly potentiates the effectiveness of tigecycline against pathogenic E. coli and K. pneumoniae. Overall, this work highlights the therapeutic potential of Met in overcoming plasmid-mediated high-level tigecycline resistance, and provides a new paradigm to enhance antibiotic efficacy by harnessing cellular metabolic networks as well as epigenetic modifications.

Host RNA N6-methyladenosine and incoming DNA N6-methyldeoxyadenosine modifications cooperatively elevate the condensation potential of DNA to activate immune surveillance

Self-non-self discrimination is fundamental to life, thereby even microbes can apply DNA modifications to recognize non-self DNA. However, mammalian cytosolic DNA sensors indiscriminately bind DNA, necessitating specific mechanism(s) for self-non-self discrimination. Here, we show that mammalian RNA N6-methyladenosine (m6A) and incoming DNA N6-methyldeoxyadenosine (6mdA) cooperatively elevate the condensation potential of DNA to activate immunosurveillance. RNA m6A modification was found to enhance the activation of cyclic guanosine monophosphate-AMP synthase (cGAS) via increasing DNA phase separation. And 6mdA further increased the phase separation potential of DNA. Consistently, host RNA m6A and incoming DNA 6mdA modifications cooperatively elevated the incoming DNA condensation and cGAS activation. Moreover, we developed a prodrug, QKY-613. QKY-613 promoted a discriminative incorporation of 6mdA into viral DNAs to elevate host immune surveillance, and decreased mortality in virus-infected aged mice. Our results link nucleic acid modification diversity with immune surveillance via phase separation, which might be targeted for therapeutic intervention.

Sequencing-based analysis of microbiomes

Microbiomes occupy a range of niches and, in addition to having diverse compositions, they have varied functional roles that have an impact on agriculture, environmental sciences, and human health and disease. The study of microbiomes has been facilitated by recent technological and analytical advances, such as cheaper and higher-throughput DNA and RNA sequencing, improved long-read sequencing and innovative computational analysis methods. These advances are providing a deeper understanding of microbiomes at the genomic, transcriptional and translational level, generating insights into their function and composition at resolutions beyond the species level.

Epigenetic control of tetrapyrrole biosynthesis by m4C DNA methylation in a cyanobacterium

Epigenetic DNA modifications are pivotal in eukaryotic gene expression, but their regulatory significance in bacteria is less understood. In Synechocystis 6803, the DNA methyltransferase M.Ssp6803II modifies the first cytosine in the GGCC motif, forming N4-methylcytosine (GGm4CC). Deletion of the sll0729 gene encoding M.Ssp6803II (∆sll0729) caused a bluish phenotype due to reduced chlorophyll levels, which was reversed by suppressor mutations. Re-sequencing of 7 suppressor clones revealed a common GGCC to GGTC mutation in the slr1790 promoter's discriminator sequence, encoding protoporphyrinogen IX oxidase, HemJ, crucial for tetrapyrrole biosynthesis. Transcriptomic and qPCR analyses indicated aberrant slr1790 expression in ∆sll0729 mutants. This aberration led to the accumulation of coproporphyrin III and protoporphyrin IX, indicative of impaired HemJ activity. To confirm the importance of DNA methylation in hemJ expression, hemJ promoter variants with varying discriminator sequences were introduced into the wild type, followed by sll0729 deletion. The sll0729 deletion segregated in strains with the GGTC discriminator motif, resulting in wild-type-like pigmentation, whereas freshly prepared ∆sll0729 mutants with the native hemJ promoter exhibited the bluish phenotype. These findings demonstrate that hemJ is tightly regulated in Synechocystis and that N4-methylcytosine is essential for proper hemJ expression. Thus, cytosine N4-methylation is a relevant epigenetic marker in Synechocystis and likely other cyanobacteria.

Complete genomes of Asgard archaea reveal diverse integrated and mobile genetic elements

Asgard archaea are of great interest as the progenitors of Eukaryotes, but little is known about the mobile genetic elements (MGEs) that may shape their ongoing evolution. Here, we describe MGEs that replicate in Atabeyarchaeia, a wetland Asgard archaea lineage represented by two complete genomes. We used soil depth-resolved population metagenomic data sets to track 18 MGEs for which genome structures were defined and precise chromosome integration sites could be identified for confident host linkage. Additionally, we identified a complete 20.67 kbp circular plasmid and two family-level groups of viruses linked to Atabeyarchaeia, via CRISPR spacer targeting. Closely related 40 kbp viruses possess a hypervariable genomic region encoding combinations of specific genes for small cysteine-rich proteins structurally similar to restriction-homing endonucleases. One 10.9 kbp integrative conjugative element (ICE) integrates genomically into the Atabeyarchaeum deiterrae-1 chromosome and has a 2.5 kbp circularizable element integrated within it. The 10.9 kbp ICE encodes an expressed Type IIG restriction-modification system with a sequence specificity matching an active methylation motif identified by Pacific Biosciences (PacBio) high-accuracy long-read (HiFi) metagenomic sequencing. Restriction-modification of Atabeyarchaeia differs from that of another coexisting Asgard archaea, Freyarchaeia, which has few identified MGEs but possesses diverse defense mechanisms, including DISARM and Hachiman, not found in Atabeyarchaeia. Overall, defense systems and methylation mechanisms of Asgard archaea likely modulate their interactions with MGEs, and integration/excision and copy number variation of MGEs in turn enable host genetic versatility.

Bacterial epigenetics and its implication for agriculture, probiotics development, and biotechnology design

In their natural habitats, organisms encounter numerous external stimuli and must be able to sense and adapt to those stimuli to survive. Unlike mutations, epigenetic changes do not alter the underlying DNA sequence. Instead, they create modifications that promote or silence gene expression. Bacillus subtilis has long been a model organism in studying genetics and development. It is beneficial for numerous biotechnological applications where it is included as a probiotic, in fermentation, or in bio-concrete design. This bacterium has also emerged recently as a model organism for studying bacterial epigenetic adaptation. In this review, we examine the evolving knowledge of epigenetic regulation (restriction-modification systems (RM), orphan methyltransferases, and chromosome condensation) in B. subtilis and related bacteria, and utilize it as a case study to test their potential roles and future applications in genetic engineering and microbial biotechnology. Finally, we suggest how the implementation of these fundamental findings promotes the design of synthetic epigenetic memory circuits and their future applications in agriculture, medicine, and biotechnology.

Precision DNA methylation typing via hierarchical clustering of Nanopore current signals and attention-based neural network

Decoding DNA methylation sites through nanopore sequencing has emerged as a cutting-edge technology in the field of DNA methylation research, as it enables direct sequencing of native DNA molecules without the need for prior enzymatic or chemical treatments. During nanopore sequencing, methylation modifications on DNA bases cause changes in electrical current intensity. Therefore, constructing deep neural network models to decode the electrical signals of nanopore sequencing has become a crucial step in methylation site identification. In this study, we utilized nanopore sequencing data containing diverse DNA methylation types and motif sequence diversity. We proposed a feature encoding method based on current signal clustering and leveraged the powerful attention mechanism in the Transformer framework to construct the PoreFormer model for identifying DNA methylation sites in nanopore sequencing. The model demonstrated excellent performance under conditions of multi-class methylation and motif sequence diversity, offering new insights into related research fields.

Changes in DNA methylation contribute to rapid adaptation in bacterial plant pathogen evolution

Adaptation is usually explained by beneficial genetic mutations that are transmitted from parents to offspring and become fixed in the adapted population. However, genetic mutation analysis alone is not sufficient to fully explain the adaptive processes, and several studies report the existence of nongenetic (or epigenetic) inheritance that can enable adaptation to new environments. In the present work, we tested the hypothesis of the role of DNA methylation, a form of epigenetic modification, in adaptation of the plant pathogen Ralstonia pseudosolanacearum to the host during experimental evolution. Using SMRT-seq technology, we analyzed the methylomes of 31 experimentally evolved clones obtained after serial passages on 5 different plant species during 300 generations. Comparison with the methylome of the ancestral clone revealed a list of 50 differential methylated sites (DMSs) at the GTWWAC motif. Gene expression analysis of the 39 genes targeted by these DMSs revealed limited correlation between differential methylation and differential expression of the corresponding genes. Only 1 gene showed a correlation, the RSp0338 gene encoding the EpsR regulator protein. The MSRE-qPCR technology, used as an alternative approach for DNA methylation analysis, also found the 2 DMSs upstream RSp0338. Using site-directed mutagenesis, we demonstrated the contribution of these 2 DMSs in host adaptation. As these DMSs appeared very early in the experimental evolution, we hypothesize that such fast epigenetic changes can allow rapid adaptation to the plant stem environment. In addition, we found that the change in DNA methylation upstream RSp0338 remains stable at least for 100 generations outside the host and thus can contribute to long-term adaptation to the host plant. To our knowledge, this is the first study showing a direct link between bacterial epigenetic variation and adaptation to a new environment.

Depletion of CpG dinucleotides in bacterial genomes may represent an adaptation to high temperatures

Dinucleotide biases have been widely investigated in the genomes of eukaryotes and viruses, but not in bacteria. We assembled a dataset of bacterial genomes (>15 000), which are representative of the genetic diversity in the kingdom Eubacteria, and we analyzed dinucleotide biases in relation to different traits. We found that TpA dinucleotides are the most depleted and that CpG dinucleotides show the widest dispersion. The abundances of both dinucleotides vary with genomic G + C content and show a very strong phylogenetic signal. After accounting for G + C content and phylogenetic inertia, we analyzed different bacterial lifestyle traits. We found that temperature preferences associate with the abundance of CpG dinucleotides, with thermophiles/hyperthemophiles being particularly depleted. Conversely, the TpA dinucleotide displays a bias that only depends on genomic G + C composition. Using predictions of intrinsic cyclizability we also show that CpG depletion may associate with higher DNA bendability in both thermophiles/hyperthermophiles and mesophiles, and that the former are predicted to have significantly more flexible genomes than the latter. We suggest that higher bendability is advantageous at high temperatures because it facilitates DNA positive supercoiling and that, through modulation of DNA mechanical properties, local or global CpG depletion controls genome organization, most likely not only in bacteria.

Bacmethy: A novel and convenient tool for investigating bacterial DNA methylation pattern and their transcriptional regulation effects

DNA methylation serves as the primary mode of epigenetic regulation in prokaryotes, particularly through transcriptional regulation. With the rapid implementation of third-generation sequencing technology, we are currently experiencing a golden age of bacterial epigenomics. However, there has been a lack of comprehensive research exploring the versatility and consequential impact of bacterial DNA methylome on cellular and physiological functions. There is a critical need for a user-friendly bioinformatics tool that can effectively characterize DNA methylation modification features and predict the regulation patterns. To address this gap, the current study introduces Bacmethy, an innovative tool that utilizes SMRT-seq data and offers a range of analytical modules. First, the tool classifies methylation sites in the genome, highlighting the distinct regulations present under varying modification fractions and location enrichment. Furthermore, this tool enables us to identify regulatory region methylation and potential cis and trans interactions between methylation sites and regulatory effectors. Using benchmark data sets and our data, we show that our tool facilitates the understanding of the distinctive traits of DNA methylation modifications and predicts transcriptional regulation effects on important physiological and pathological functions. Bacmethy code is freely available, and the Docker image is downloadable. Bacmethy has been made available as a user-friendly web server interface at https://bacmethy.med.sustech.edu.cn.

Restriction modification systems in archaea: A panoramic outlook

Restriction modification (RM) systems are one of the ubiquitous yet primitive defense responses employed by bacteria and archaea with the primary role of safeguarding themselves against invading bacteriophages. Protection of the host occurs by the cleavage of the invading foreign DNA via restriction endonucleases with concomitant methylation of host DNA with the aid of a methyltransferase counterpart. RM systems have been extensively studied in bacteria, however, in the case of archaea there are limited reports of RM enzymes that are investigated to date owing to their inhospitable growth demands. This review aims to broaden the knowledge about what is known about the diversity of RM systems in archaea and encapsulate the current knowledge on restriction and modification enzymes characterized in archaea so far and the role of RM systems in the milieu of archaeal biology.

Prophage-encoded methyltransferase drives adaptation of community-acquired methicillin-resistant Staphylococcus aureus

We recently described the evolution of a community-acquired methicillin-resistant Staphylococcus aureus (CA-MRSA) USA300 variant responsible for an outbreak of skin and soft tissue infections. Acquisition of a mosaic version of the Φ11 prophage (mΦ11) that increases skin abscess size was an early step in CA-MRSA adaptation that primed the successful spread of the clone. The present report shows how prophage mΦ11 exerts its effect on virulence for skin infection without encoding a known toxin or fitness genes. Abscess size and skin inflammation were associated with DNA methylase activity of an mΦ11-encoded adenine methyltransferase (designated pamA). pamA increased expression of fibronectin-binding protein A (fnbA; FnBPA), and inactivation of fnbA eliminated the effect of pamA on abscess virulence without affecting strains lacking pamA. Thus, fnbA is a pamA-specific virulence factor. Mechanistically, pamA was shown to promote biofilm formation in vivo in skin abscesses, a phenotype linked to FnBPA's role in biofilm formation. Collectively, these data reveal a novel mechanism-epigenetic regulation of staphylococcal gene expression-by which phage can regulate virulence to drive adaptive leaps by S. aureus.

Important role of DNA methylation hints at significant potential in tuberculosis

Tuberculosis (TB), an infectious disease caused by Mycobacterium tuberculosis (Mtb) infection, has persisted as a major global public health threat for millennia. Until now, TB continues to challenge efforts aimed at controlling it, with drug resistance and latent infections being the two main factors hindering treatment efficacy. The scientific community is still striving to understand the underlying mechanisms behind Mtb's drug resistance and latent infection. DNA methylation, a critical epigenetic modification occurring throughout an individual's growth and development, has gained attention following advances in high-throughput sequencing technologies. Researchers have observed abnormal DNA methylation patterns in the host genome during Mtb infection. Given the escalating issue of drug-resistant Mtb, delving into the role of DNA methylation in TB's development is crucial. This review article explores DNA methylation's significance in human growth, development and disease, and its role in regulating Mtb's evolution and infection processes. Additionally, it discusses potential applications of DNA methylation research in tuberculosis.

Dinucleotide biases in the genomes of prokaryotic and eukaryotic dsDNA viruses and their hosts

The genomes of cellular organisms display CpG and TpA dinucleotide composition biases. Such biases have been poorly investigated in dsDNA viruses. Here, we show that in dsDNA virus, bacterial, and eukaryotic genomes, the representation of TpA and CpG dinucleotides is strongly dependent on genomic G + C content. Thus, the classical observed/expected ratios do not fully capture dinucleotide biases across genomes. Because a larger portion of the variance in TpA frequency was explained by G + C content, we explored which additional factors drive the distribution of CpG dinucleotides. Using the residuals of the linear regressions as a measure of dinucleotide abundance and ancestral state reconstruction across eukaryotic and prokaryotic virus trees, we identified an important role for phylogeny in driving CpG representation. Nonetheless, phylogenetic ANOVA analyses showed that few host associations also account for significant variations. Among eukaryotic viruses, most significant differences were observed between arthropod-infecting viruses and viruses that infect vertebrates or unicellular organisms. However, an effect of viral DNA methylation status (either driven by the host or by viral-encoded methyltransferases) is also likely. Among prokaryotic viruses, cyanobacteria-infecting phages resulted to be significantly CpG-depleted, whereas phages that infect bacteria in the genera Burkolderia and Staphylococcus were CpG-rich. Comparison with bacterial genomes indicated that this effect is largely driven by the general tendency for phages to resemble the host's genomic CpG content. Notably, such tendency is stronger for temperate than for lytic phages. Our data shed light into the processes that shape virus genome composition and inform manipulation strategies for biotechnological applications.

Genome-wide DNA N6-methyladenosine in Aeromonas veronii and Helicobacter pylori

BACKGROUND

DNA N6-methyladenosine (6mA), as an important epigenetic modification, widely exists in bacterial genomes and participates in the regulation of toxicity, antibiotic resistance, and antioxidant. With the continuous development of sequencing technology, more 6mA sites have been identified in bacterial genomes, but few studies have focused on the distribution characteristics of 6mA at the whole-genome level and its association with gene expression and function.

RESULTS

This study conducted an in-depth analysis of the 6mA in the genomes of two pathogenic bacteria, Aeromonas veronii and Helicobacter pylori. The results showed that the 6mA was widely distributed in both strains. In A. veronii, 6mA sites were enriched at 3' end of protein-coding genes, exhibiting a certain inhibitory effect on gene expression. Genes with low 6mA density were associated with cell motility. While in H. pylori, 6mA sites were enriched at 5' end of protein-coding genes, potentially enhancing gene expression. Genes with low 6mA density were closely related to defense mechanism.

CONCLUSIONS

This study elucidated the distribution characteristics of 6mA in A. veronii and H. pylori, highlighting the effects of 6mA on gene expression and function. These findings provide valuable insights into the epigenetic regulation and functional characteristics of A. veronii and H. pylori.

DNA Sequence Control of Enzyme Filamentation and Activation of the SgrAI Endonuclease

Enzyme polymerization (also known as filamentation) has emerged as a new layer of enzyme regulation. SgrAI is a sequence-dependent DNA endonuclease that forms polymeric filaments with enhanced DNA cleavage activity as well as altered DNA sequence specificity. To better understand this unusual regulatory mechanism, full global kinetic modeling of the reaction pathway, including the enzyme filamentation steps, has been undertaken. Prior work with the primary DNA recognition sequence cleaved by SgrAI has shown how the kinetic rate constants of each reaction step are tuned to maximize activation and DNA cleavage while minimizing the extent of DNA cleavage to the host genome. In the current work, we expand on our prior study by now including DNA cleavage of a secondary recognition sequence, to understand how the sequence of the bound DNA modulates filamentation and activation of SgrAI. The work shows that an allosteric equilibrium between low and high activity states is modulated by the sequence of bound DNA, with primary sequences more prone to activation and filament formation, while SgrAI bound to secondary recognition sequences favor the low (and nonfilamenting) state by up to 40-fold. In addition, the degree of methylation of secondary sequences in the host organism, Streptomyces griseus, is now reported for the first time and shows that as predicted, these sequences are left unprotected from the SgrAI endonuclease making sequence specificity critical in this unusual filament-forming enzyme.

Characterization of Limnospira platensis PCC 9108 R-M and CRISPR-Cas systems

The filamentous cyanobacterium Limnospira platensis, formerly known as Arthrospira platensis or spirulina, is one of the most commercially important species of microalgae. Due to its high nutritional value, pharmacological and industrial applications it is extensively cultivated on a large commercial scale. Despite its widespread use, its precise manipulation is still under development due to the lack of effective genetic protocols. Genetic transformation of Limnospira has been attempted but the methods reported have not been generally reproducible in other laboratories. Knowledge of the transformation defense mechanisms is essential for understanding its physiology and for broadening their applications. With the aim to understand more about the genetic defenses of L. platensis, in this work we have identified the restriction-modification and CRISPR-Cas systems and we have cloned and characterized thirteen methylases. In parallel, we have also characterized the methylome and orphan methyltransferases using genome-wide analysis of DNA methylation patterns and RNA-seq. The identification and characterization of these enzymes will be a valuable resource to know how this strain avoids being genetically manipulated and for further genomics studies.

mEnrich-seq: methylation-guided enrichment sequencing of bacterial taxa of interest from microbiome

Metagenomics has enabled the comprehensive study of microbiomes. However, many applications would benefit from a method that sequences specific bacterial taxa of interest, but not most background taxa. We developed mEnrich-seq (in which 'm' stands for methylation and seq for sequencing) for enriching taxa of interest from metagenomic DNA before sequencing. The core idea is to exploit the self versus nonself differentiation by natural bacterial DNA methylation and rationally choose methylation-sensitive restriction enzymes, individually or in combination, to deplete host and background taxa while enriching targeted taxa. This idea is integrated with library preparation procedures and applied in several applications to enrich (up to 117-fold) pathogenic or beneficial bacteria from human urine and fecal samples, including species that are hard to culture or of low abundance. We assessed 4,601 bacterial strains with mapped methylomes so far and showed broad applicability of mEnrich-seq. mEnrich-seq provides microbiome researchers with a versatile and cost-effective approach for selective sequencing of diverse taxa of interest.

The Helicobacter pylori methylome is acid-responsive due to regulation by the two-component system ArsRS and the type I DNA methyltransferase HsdM1 (HP0463)

In addition to its role in genome protection, DNA methylation can regulate gene expression. In this study, we characterized the impact of acidity, phase variation, and the ArsRS TCS on the expression of the Type I m6A DNA methyltransferase HsdM1 (HP0463) of Helicobacter pylori 26695 and their subsequent effects on the methylome. Transcription of hsdM1 increases at least fourfold in the absence of the sensory histidine kinase ArsS, the major acid-sensing protein of H. pylori. hsdM1 exists in the phase-variable operon hsdR1-hsdM1. Phase-locking hsdR1 (HP0464), the restriction endonuclease gene, has significant impacts on the transcription of hsdM1. To determine the impacts of methyltransferase transcription patterns on the methylome, we conducted methylome sequencing on samples cultured at pH 7 or pH 5. We found differentially methylated motifs between these growth conditions and that deletions of arsS and/or hsdM1 interfere with the epigenetic acid response. Deletion of arsS leads to altered activity of HsdM1 and multiple other methyltransferases under both pH conditions indicating that the ArsRS TCS, in addition to direct effects on regulon transcription during acid acclimation, may also indirectly impact gene expression via regulation of the methylome. We determined the target motif of HsdM1 (HP0463) to be the complementary bipartite sequence pair 5'-TCAm6AVN6TGY-3' and 3'-AGTN6GAm6ACA-5'. This complex regulation of DNA methyltransferases, and thus differential methylation patterns, may have implications for the decades-long persistent infection by H. pylori. IMPORTANCE This study expands the possibilities for complex, epigenomic regulation in Helicobacter pylori. We demonstrate that the H. pylori methylome is plastic and acid sensitive via the two-component system ArsRS and the DNA methyltransferase HsdM1. The control of a methyltransferase by ArsRS may allow for a layered response to changing acidity. Likely, an early response whereby ArsR~P affects regulon expression, including the methyltransferase hsdM1. Then, a somewhat later effect as the altered methylome, due to altered HsdM1 expression, subsequently alters the expression of other genes involved in acclimation. The intermediate methylation of certain motifs supports the hypothesis that methyltransferases play a regulatory role. Untangling this additional web of regulation could play a key role in understanding H. pylori colonization and persistence.

Isolation of a widespread giant virus implicated in cryptophyte bloom collapse

Photosynthetic cryptophytes are ubiquitous protists that are major participants in the freshwater phytoplankton bloom at the onset of spring. Mortality due to change in environmental conditions and grazing have been recognized as key factors contributing to bloom collapse. In contrast, the role of viral outbreaks as factors terminating phytoplankton blooms remains unknown from freshwaters. Here, we isolated and characterized a cryptophyte virus contributing to the annual collapse of a natural cryptophyte spring bloom population. This viral isolate is also representative for a clade of abundant giant viruses (phylum Nucleocytoviricota) found in freshwaters all over the world.

A sequential one-pot approach for rapid and convenient characterization of putative restriction-modification systems

IMPORTANCE

The elucidation of the molecular basis of virus-host coevolutionary interactions is boosted with state-of-the-art sequencing technologies. However, the sequence-only information is often insufficient to output a conclusive argument without biochemical characterizations. We proposed a 1-day and one-pot approach to confirm the exact function of putative restriction-modification (R-M) genes that presumably mediate microbial coevolution. The experiments mainly focused on a series of putative R-M enzymes from a deep-sea virus and its host bacterium. The results quickly unveiled unambiguous substrate specificities, superior catalytic performance, and unique sequence preferences for two new restriction enzymes (capable of cleaving DNA) and two new methyltransferases (capable of modifying DNA with methyl groups). The reality of the functional R-M system reinforced a model of mutually beneficial interactions with the virus in the deep-sea microbial ecosystem. The cell culture-independent approach also holds great potential for exploring novel and biotechnologically significant R-M enzymes from microbial dark matter.

PacBio sequencing of human fecal samples uncovers the DNA methylation landscape of 22 673 gut phages

Gut phages have an important impact on human health. Methylation plays key roles in DNA recognition, gene expression regulation and replication for phages. However, the DNA methylation landscape of gut phages is largely unknown. Here, with PacBio sequencing (2120×, 4785 Gb), we detected gut phage methylation landscape based on 22 673 gut phage genomes, and presented diverse methylation motifs and methylation differences in genomic elements. Moreover, the methylation rate of phages was associated with taxonomy and host, and N6-methyladenine methylation rate was higher in temperate phages than in virulent phages, suggesting an important role for methylation in phage-host interaction. In particular, 3543 (15.63%) phage genomes contained restriction-modification system, which could aid in evading clearance by the host. This study revealed the DNA methylation landscape of gut phage and its potential roles, which will advance the understanding of gut phage survival and human health.

Enhanced transformation efficiency in Treponema denticola enabled by SyngenicDNA-based plasmids lacking restriction-modification target motifs

Oral spirochetes are among a small group of keystone pathogens contributing to dysregulation of tissue homeostatic processes that leads to breakdown of the tissue and bone supporting the teeth in periodontal disease. Additionally, our group has recently demonstrated that Treponema are among the dominant microbial genera detected intracellularly in tumor specimens from patients with oral squamous cell carcinoma. While over 60 species and phylotypes of oral Treponema have been detected, T. denticola is one of the few that can be grown in culture and the only one in which genetic manipulation is regularly performed. Thus, T. denticola is a key model organism for studying spirochete metabolic processes, interactions with other microbes, and host cell and tissue responses relevant to oral diseases, as well as venereal and nonvenereal treponematoses whose agents lack workable genetic systems. We previously demonstrated improved transformation efficiency using an Escherichia coli-T. denticola shuttle plasmid and its utility for expression in T. denticola of an exogenous fluorescent protein that is active under anaerobic conditions. Here, we expand on this work by characterizing T. denticola Type I and Type II restriction-modification (R-M) systems and designing a high-efficiency R-M-silent "SyngenicDNA" shuttle plasmid resistant to all T. denticola ATCC 35405 R-M systems. Resequencing of the ATCC 33520 genome revealed an additional Type I R-M system consistent with the relatively low transformation efficiency of the shuttle plasmid in this strain. Using SyngenicDNA approaches, we optimized shuttle plasmid transformation efficiency in T. denticola and used it to complement a defined T. denticola ΔfhbB mutant strain. We further report the first high-efficiency transposon mutagenesis of T. denticola using an R-M-silent, codon-optimized, himarC9 transposase-based plasmid. Thus, use of SyngenicDNA-based strategies and tools can enable further mechanistic examinations of T. denticola physiology and behavior.

4mC-CGRU: Identification of N4-Methylcytosine (4mC) sites using convolution gated recurrent unit in Rosaceae genome

An epigenetic modification is DNA N4-methylcytosine (4mC) that affects several biological functions without altering the DNA nucleotides, including DNA conformation, cell development, replication, stability, and DNA structural changes. To prevent restriction enzyme from damaging self-DNA, 4mC performs a critical role in restriction-modification functions. Existing studies mainly focused on finding hand-crafted features to identify 4mC locations, but these methods are inefficient due to high time consuming and high costs. In our research work, we propose a 4mC-CGRU which is a deep learning-based computational model with a standard encoding method to identify the 4mC sites from DNA sequences that learned autonomous feature selection in the Rosaceae genome, particularly in Rosa chinensis (R. chinensis) and Fragaria vesca (F. vesca). The proposed model consists of a convolutional neural network (CNN) and a gated recurrent unit network (GRU)-based model for identifying 4mC sites from Fragaria vesca and Rosa chinensis in the genomes. The CNN model extracts useful features from the datasets and the GRU classifies the DNA sequences. Thus, our approach can automatically extract important features to detect relative sites from DNA sequence. The performance analysis shows that the proposed model consistently outperforms over the state-of-the-art works in detecting 4mC sites.

Differential adenine methylation analysis reveals increased variability in 6mA in the absence of methyl-directed mismatch repair

IMPORTANCE

Methylation greatly influences the bacterial genome by guiding DNA repair and regulating pathogenic and stress-response phenotypes. But, the rate of epigenetic changes and their consequences on molecular phenotypes are underexplored. Through a detailed characterization of genome-wide adenine methylation in a commonly used laboratory strain of Escherichia coli, we reveal that mismatch repair deficient populations experience an increase in epimutations resulting in a genome-wide reduction of 6mA methylation in a manner consistent with genetic drift. Our findings highlight how methylation patterns evolve and the constraints on epigenetic evolution due to post-replicative DNA repair, contributing to a deeper understanding of bacterial genome evolution and how epimutations may introduce semi-permanent variation that can influence adaptation.

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.

Long-Read Sequencing Reveals Extensive DNA Methylations in Human Gut Phagenome Contributed by Prevalently Phage-Encoded Methyltransferases

DNA methylation plays a crucial role in the survival of bacteriophages (phages), yet the understanding of their genome methylation remains limited. In this study, DNA methylation patterns are analyzed in 8848 metagenome-assembled high-quality phages from 104 fecal samples using single-molecule real-time sequencing. The results demonstrate that 97.60% of gut phages exhibit methylation, with certain factors correlating with methylation densities. Phages with higher methylation densities appear to have potential viability advantages. Strikingly, more than one-third of the phages possess their own DNA methyltransferases (MTases). Increased MTase copies are associated with higher genome methylation densities, specific methylation motifs, and elevated prevalence of certain phage groups. Notably, the majority of these MTases share close homology with those encoded by gut bacteria, suggesting their exchange during phage-bacterium interactions. Furthermore, these MTases can be employed to accurately predict phage-host relationships. Overall, the findings indicate the widespread utilization of DNA methylation by gut DNA phages as an evasion mechanism against host defense systems, with a substantial contribution from phage-encoded MTases.

Methylome evolution suggests lineage-dependent selection in the gastric pathogen Helicobacter pylori

The bacterial pathogen Helicobacter pylori, the leading cause of gastric cancer, is genetically highly diverse and harbours a large and variable portfolio of restriction-modification systems. Our understanding of the evolution and function of DNA methylation in bacteria is limited. Here, we performed a comprehensive analysis of the methylome diversity in H. pylori, using a dataset of 541 genomes that included all known phylogeographic populations. The frequency of 96 methyltransferases and the abundance of their cognate recognition sequences were strongly influenced by phylogeographic structure and were inter-correlated, positively or negatively, for 20% of type II methyltransferases. Low density motifs were more likely to be affected by natural selection, as reflected by higher genomic instability and compositional bias. Importantly, direct correlation implied that methylation patterns can be actively enriched by positive selection and suggests that specific sites have important functions in methylation-dependent phenotypes. Finally, we identified lineage-specific selective pressures modulating the contraction and expansion of the motif ACGT, revealing that the genetic load of methylation could be dependent on local ecological factors. Taken together, natural selection may shape both the abundance and distribution of methyltransferases and their specific recognition sequences, likely permitting a fine-tuning of genome-encoded functions not achievable by genetic variation alone.

Roving methyltransferases generate a mosaic epigenetic landscape and influence evolution in Bacteroides fragilis group

Three types of DNA methyl modifications have been detected in bacterial genomes, and mechanistic studies have demonstrated roles for DNA methylation in physiological functions ranging from phage defense to transcriptional control of virulence and host-pathogen interactions. Despite the ubiquity of methyltransferases and the immense variety of possible methylation patterns, epigenomic diversity remains unexplored for most bacterial species. Members of the Bacteroides fragilis group (BFG) reside in the human gastrointestinal tract as key players in symbiotic communities but also can establish anaerobic infections that are increasingly multi-drug resistant. In this work, we utilize long-read sequencing technologies to perform pangenomic (n = 383) and panepigenomic (n = 268) analysis of clinical BFG isolates cultured from infections seen at the NIH Clinical Center over four decades. Our analysis reveals that single BFG species harbor hundreds of DNA methylation motifs, with most individual motif combinations occurring uniquely in single isolates, implying immense unsampled methylation diversity within BFG epigenomes. Mining of BFG genomes identified more than 6000 methyltransferase genes, approximately 1000 of which were associated with intact prophages. Network analysis revealed substantial gene flow among disparate phage genomes, implying a role for genetic exchange between BFG phages as one of the ultimate sources driving BFG epigenome diversity.

Shaping eukaryotic epigenetic systems by horizontal gene transfer

DNA methylation constitutes one of the pillars of epigenetics, relying on covalent bonds for addition and/or removal of chemically distinct marks within the major groove of the double helix. DNA methyltransferases, enzymes which introduce methyl marks, initially evolved in prokaryotes as components of restriction-modification systems protecting host genomes from bacteriophages and other invading foreign DNA. In early eukaryotic evolution, DNA methyltransferases were horizontally transferred from bacteria into eukaryotes several times and independently co-opted into epigenetic regulatory systems, primarily via establishing connections with the chromatin environment. While C5-methylcytosine is the cornerstone of plant and animal epigenetics and has been investigated in much detail, the epigenetic role of other methylated bases is less clear. The recent addition of N4-methylcytosine of bacterial origin as a metazoan DNA modification highlights the prerequisites for foreign gene co-option into the host regulatory networks, and challenges the existing paradigms concerning the origin and evolution of eukaryotic regulatory systems.

Navigating the pitfalls of mapping DNA and RNA modifications

Chemical modifications to nucleic acids occur across the kingdoms of life and carry important regulatory information. Reliable high-resolution mapping of these modifications is the foundation of functional and mechanistic studies, and recent methodological advances based on next-generation sequencing and long-read sequencing platforms are critical to achieving this aim. However, mapping technologies may have limitations that sometimes lead to inconsistent results. Some of these limitations are technical in nature and specific to certain types of technology. Here, however, we focus on common (yet not always widely recognized) pitfalls that are shared among frequently used mapping technologies and discuss strategies to help technology developers and users mitigate their effects. Although the emphasis is primarily on DNA modifications, RNA modifications are also discussed.

The extremophile Picrophilus torridus carries a DNA adenine methylase M.PtoI that is part of a Type I restriction-modification system

DNA methylation events mediated by orphan methyltransferases modulate various cellular processes like replication, repair and transcription. Bacteria and archaea also harbor DNA methyltransferases that are part of restriction-modification systems, which serve to protect the host genome from being cleaved by the cognate restriction enzyme. While DNA methylation has been exhaustively investigated in bacteria it remains poorly understood in archaea. Picrophilus torridus is a euryarchaeon that can thrive under conditions of extremely low pH (0.7), and thus far no reports have been published regarding DNA methylation in this extremophile. This study reports the first experimentation examining DNA methylation in P. torridus. We find the genome to carry methylated adenine (m6A) but not methylated cytosine (m5C) residues. The m6A modification is absent at GATC sites, indicating the absence of an active Dam methylase even though the dam gene has been annotated in the genome sequence. Two other methylases have also been annotated in the P. torridus genome sequence. One of these is a part of a Type I restriction-modification system. Considering that all Type I modification methylases characterized to date target adenine residues, the modification methylase of this Type I system has been examined. The genes encoding the S subunit (that is responsible for DNA recognition) and M subunit (that is responsible for DNA methylation) have been cloned and the recombinant protein purified from E.coli, and regions involved in M-S interactions have been identified. The M.PtoI enzyme harbors all the motifs that typify Type I modification methylases, and displays robust adenine methylation in in vitro assays under a variety of conditions. Interestingly, magnesium is essential for enzyme activity. The enzyme displays substrate inhibition at higher concentrations of AdoMet. Mutational analyses reveal that Motif I plays a role in AdoMet binding, and Motif IV is critical for methylation activity. The data presented here lays the foundation for further research in the area of DNA methylation and restriction-modification research in this most unusual microorganism.

DNA Methylation and RNA-DNA Hybrids Regulate the Single-Molecule Localization of a DNA Methyltransferase on the Bacterial Nucleoid

Bacterial DNA methyltransferases (MTases) function in restriction modification systems, cell cycle control, and the regulation of gene expression. DnmA is a recently described DNA MTase that forms N6-methyladenosine at nonpalindromic 5'-GACGAG-3' sites in Bacillus subtilis, yet how DnmA activity is regulated is unknown. To address DnmA regulation, we tested substrate binding in vitro and found that DnmA binds poorly to methylated DNA and to an RNA-DNA hybrid with the DNA recognition sequence. Further, DnmA variants with amino acid substitutions that disrupt cognate sequence recognition or catalysis also bind poorly to DNA. Using superresolution fluorescence microscopy and single-molecule tracking of DnmA-PAmCherry, we characterized the subcellular DnmA diffusion and detected its preferential localization to the replisome region and the nucleoid. Under conditions where the chromosome is highly methylated, upon RNA-DNA hybrid accumulation, or with a DnmA variant with severely limited DNA binding activity, DnmA is excluded from the nucleoid, demonstrating that prior methylation or accumulation of RNA-DNA hybrids regulates the association of DnmA with the chromosome in vivo. Furthermore, despite the high percentage of methylated recognition sites and the proximity to putative endonuclease genes conserved across bacterial species, we find that DnmA fails to protect B. subtilis against phage predation, suggesting that DnmA is functionally an orphan MTase involved in regulating gene expression. Our work explores the regulation of a bacterial DNA MTase and identifies prior methylation and RNA-DNA hybrids as regulators of MTase localization. These MTase regulatory features could be common across biology. IMPORTANCE DNA methyltransferases (MTases) influence gene expression, cell cycle control, and host defense through DNA modification. Predicted MTases are pervasive across bacterial genomes, but the vast majority remain uncharacterized. Here, we show that in the soil microorganism Bacillus subtilis, the DNA MTase dnmA and neighboring genes are remnants of a phage defense system that no longer protects against phage predation. This result suggests that portions of the bacterial methylome may originate from inactive restriction modification systems that have maintained methylation activity. Analysis of DnmA movement in vivo shows that active DnmA localizes in the nucleoid, suggesting that DnmA can search for recognition sequences throughout the nucleoid region with some preference for the replisome. Our results further show that prior DNA methylation and RNA-DNA hybrids regulate DnmA dynamics and nucleoid localization, providing new insight into how DNA methylation is coordinated within the cellular environment.

Bacterial DNA methyltransferase: A key to the epigenetic world with lessons learned from proteobacteria

Epigenetics modulates expression levels of various important genes in both prokaryotes and eukaryotes. These epigenetic traits are heritable without any change in genetic DNA sequences. DNA methylation is a universal mechanism of epigenetic regulation in all kingdoms of life. In bacteria, DNA methylation is the main form of epigenetic regulation and plays important roles in affecting clinically relevant phenotypes, such as virulence, host colonization, sporulation, biofilm formation et al. In this review, we survey bacterial epigenomic studies and focus on the recent developments in the structure, function, and mechanism of several highly conserved bacterial DNA methylases. These methyltransferases are relatively common in bacteria and participate in the regulation of gene expression and chromosomal DNA replication and repair control. Recent advances in sequencing techniques capable of detecting methylation signals have enabled the characterization of genome-wide epigenetic regulation. With their involvement in critical cellular processes, these highly conserved DNA methyltransferases may emerge as promising targets for developing novel epigenetic inhibitors for biomedical applications.

Identification and Mechanism of Action of the Global Secondary Metabolism Regulator SaraC in Stereum hirsutum

DNA methylation is an important factor in the regulation of gene expression. In analyzing genomic data of Stereum hirsutum FP-91666, we found a hypothetical bifunctional transcription regulator/O6Meguanine-DNA methyltransferase (named SaraC), which is widely present in both bacteria and fungi, and confirmed that its function in bacteria is mainly for DNA reparation. In this paper, we confirmed that SaraC has the function of DNA binding and demethylation through surface plasma resonance and reaction experiments in vitro. Then, we achieved the overexpression of SaraC (OES) in S. hirsutum, sequenced the methylation and transcription levels of the whole-genome, and further conducted untargeted metabolomics analyses of the OES transformants and the wild type (WT). The results confirmed that the overall-methylation levels of the transformants were significantly downregulated, and various genes related to secondary metabolism were upregulated. Through comparative untargeted metabolomic analyses, it showed that OES SA6 transformant produced a greater number of hybrid polyketides, and we identified 2 novel hybrid polyketides from the fermentation products of SA6. Our results show that overexpression SaraC can effectively stimulate the expression of secondary-metabolism-related genes, which could be a broad-spectrum tool for discovery of metabolites due to its cross-species conservation. IMPORTANCE Fungi are one of the important sources of active compounds. However, in fungi, most of the secondary metabolic biosynthetic gene clusters are weakly expressed or silenced under conventional culture conditions. How to efficiently excavate potential new compounds contained in fungi is becoming a research hot spot in the world. In this study, we found a DNA demethylation protein (SaraC) and confirmed that it is a global secondary metabolism regulator in Stereum hirsutum FP-91666. In the past, SaraC-like proteins were mainly regarded as DNA repair proteins, but our findings proved that it will be a powerful tool for mining secondary metabolites for overexpression of SaraC, which can effectively stimulate the expression of genes related to secondary metabolism.

Growth condition-dependent differences in methylation imply transiently differentiated DNA methylation states in Escherichia coli

DNA methylation in bacteria frequently serves as a simple immune system, allowing recognition of DNA from foreign sources, such as phages or selfish genetic elements. However, DNA methylation also affects other cell phenotypes in a heritable manner (i.e. epigenetically). While there are several examples of methylation affecting transcription in an epigenetic manner in highly localized contexts, it is not well-established how frequently methylation serves a more general epigenetic function over larger genomic scales. To address this question, here we use Oxford Nanopore sequencing to profile DNA modification marks in three natural isolates of Escherichia coli. We first identify the DNA sequence motifs targeted by the methyltransferases in each strain. We then quantify the frequency of methylation at each of these motifs across the entire genome in different growth conditions. We find that motifs in specific regions of the genome consistently exhibit high or low levels of methylation. Furthermore, we show that there are replicable and consistent differences in methylated regions across different growth conditions. This suggests that during growth, E. coli transiently differentiate into distinct methylation states that depend on the growth state, raising the possibility that measuring DNA methylation alone can be used to infer bacterial growth states without additional information such as transcriptome or proteome data. These results show the utility of using Oxford Nanopore sequencing as an economic means to infer DNA methylation status. They also provide new insights into the dynamics of methylation during bacterial growth and provide evidence of differentiated cell states, a transient analog to what is observed in the differentiation of cell types in multicellular organisms.

Role of DNA modifications in Mycoplasma gallisepticum

The epigenetics of bacteria, and bacteria with a reduced genome in particular, is of great interest, but is still poorly understood. Mycoplasma gallisepticum, a representative of the class Mollicutes, is an excellent model of a minimal cell because of its reduced genome size, lack of a cell wall, and primitive cell organization. In this study we investigated DNA modifications of the model object Mycoplasma gallisepticum and their roles. We identified DNA modifications and methylation motifs in M. gallisepticum S6 at the genome level using single molecule real time (SMRT) sequencing. Only the ANCNNNNCCT methylation motif was found in the M. gallisepticum S6 genome. The studied bacteria have one functional system for DNA modifications, the Type I restriction-modification (RM) system, MgaS6I. We characterized its activity, affinity, protection and epigenetic functions. We demonstrated the protective effects of this RM system. A common epigenetic signal for bacteria is the m6A modification we found, which can cause changes in DNA-protein interactions and affect the cell phenotype. Native methylation sites are underrepresented in promoter regions and located only near the -35 box of the promoter, which does not have a significant effect on gene expression in mycoplasmas. To study the epigenetics effect of m6A for genome-reduced bacteria, we constructed a series of M. gallisepticum strains expressing EGFP under promoters with the methylation motifs in their different elements. We demonstrated that m6A modifications of the promoter located only in the -10-box affected gene expression and downregulated the expression of the corresponding gene.