Restriction-modification System Research Articles

Structure and functions of a multireplicon genome of Antarctic Psychrobacter sp. ANT_H3: characterization of the genetic modules suitable for the construction of the plasmid-vectors for cold-active bacteria

Among Psychrobacter spp., there are several multireplicon strains, carrying more than two plasmids. Psychrobacter sp. ANT_H3 carries as many as 11 extrachromosomal replicons, which is the highest number in Psychrobacter spp. Plasmids of this strain were subjected to detailed genomic analysis, which enables an insight into the structure and functioning of this multireplicon genome. The replication and conjugal transfer modules of ANT_H3 plasmids were analyzed functionally to discover their potential for being used as building blocks for the construction of novel plasmid-vectors for cold-active bacteria. It was shown that two plasmids have a narrow host range as they were not able to replicate in species other than Psychrobacter, while remaining plasmids had a wider host range and were functional in various Alpha- and Gammaproteobacteria. Moreover, it was confirmed that mobilization modules of seven plasmids were functional, i.e., could be mobilized for conjugal transfer by the RK2 conjugation system. Auxiliary genes were also distinguished in ANT_H3 plasmids, including these encoding putative DNA-protecting protein DprA, multidrug efflux SMR transporter of EmrE family, glycine cleavage system T protein, MscS small-conductance mechanosensitive channel protein, and two type II restriction-modification systems. Finally, all genome-retrieved plasmids of Psychrobacter spp. were subjected to complex genome- and proteome-based comparative analyses showing that Antarctic replicons are significantly different from plasmids from other locations.

The association of prokaryotic antiviral systems and symbiotic phage communities in drinking water microbiomes

Prokaryotic antiviral systems are important mediators for prokaryote-phage interactions, which have significant implications for the survival of prokaryotic community. However, the prokaryotic antiviral systems under environmental stress are poorly understood, limiting the understanding of microbial adaptability. Here, we systematically investigated the profile of the prokaryotic antiviral systems at the community level and prokaryote-phage interactions in the drinking water microbiome. Chlorine disinfectant was revealed as the main ecological driver for the difference in prokaryotic antiviral systems and prokaryote-phage interactions. Specifically, the prokaryotic antiviral systems in the microbiome exhibited a higher abundance, broader antiviral spectrum, and lower metabolic burden under disinfectant stress. Moreover, significant positive correlations were observed between phage lysogenicity and enrichment of antiviral systems (e.g., Type IIG and IV restriction-modification (RM) systems, and Type II CRISPR-Cas system) in the presence of disinfection, indicating these antiviral systems might be more compatible with lysogenic phages and prophages. Accordingly, there was a stronger prokaryote-phage symbiosis in disinfected microbiome, and the symbiotic phages carried more auxiliary metabolic genes (AMGs) related to prokaryotic adaptability as well as antiviral systems, which might further enhance prokaryote survival in drinking water distribution systems. Overall, this study demonstrates that the prokaryotic antiviral systems had a close association with their symbiotic phages, which provides novel insights into prokaryote-phage interactions and microbial environmental adaptation.

Uncovering the link between the SpnIII restriction modification system and LuxS in Streptococcus pneumoniae meningitis isolates.

Streptococcus pneumoniae is capable of randomly switching their genomic DNA methylation pattern between six distinct bacterial subpopulations (A-F) via recombination of a type 1 restriction-modification locus, spnIII. These pneumococcal subpopulations exhibit phenotypic changes which favor carriage or invasive disease. In particular, the spnIIIB allele has been associated with increased nasopharyngeal carriage and the downregulation of the luxS gene. The LuxS/AI-2 QS system represent a universal language for bacteria and has been linked to virulence and biofilm formation in S. pneumoniae. In this work, we have explored the link between spnIII alleles, the luxS gene and virulence in two clinical pneumococcal isolates from the blood and cerebrospinal fluid (CSF) of one pediatric meningitis patient. The blood and CSF strains showed different virulence profiles in mice. Analysis of the spnIII system of these strains recovered from the murine nasopharynx showed that the system switched to different alleles commensurate with the initial source of the isolate. Of note, the blood strain showed high expression of spnIIIB allele, previously linked with less LuxS protein production. Importantly, strains with deleted luxS displayed different phenotypic profiles compared to the wildtype, but similar to the strains recovered from the nasopharynx of infected mice. This study used clinically relevant S. pneumoniae strains to demonstrate that the regulatory network between luxS and the type 1 restriction-modification system play a key role in infections and may support different adaptation to specific host niches.

Broadness and specificity: ArdB, ArdA, and Ocr against various restriction-modification systems.

ArdB, ArdA, and Ocr proteins inhibit the endonuclease activity of the type I restriction-modification enzymes (RMI). In this study, we evaluated the ability of ArdB, ArdA, and Ocr to inhibit different subtypes of Escherichia coli RMI systems (IA, IB, and IC) as well as two Bacillus licheniformis RMI systems. Furthermore we explored, the antirestriction activity of ArdA, ArdB, and Ocr against a type III restriction-modification system (RMIII) EcoPI and BREX. We found that DNA-mimic proteins, ArdA and Ocr exhibit different inhibition activity, depending on which RM system tested. This effect might be linked to the DNA mimicry nature of these proteins. In theory, DNA-mimic might competitively inhibit any DNA-binding proteins; however, the efficiency of inhibition depend on the ability to imitate the recognition site in DNA or its preferred conformation. In contrast, ArdB protein with an undescribed mechanism of action, demonstrated greater versatility against various RMI systems and provided similar antirestriction efficiency regardless of the recognition site. However, ArdB protein could not affect restriction systems that are radically different from the RMI such as BREX or RMIII. Thus, we assume that the structure of DNA-mimic proteins allows for selective inhibition of any DNA-binding proteins depending on the recognition site. In contrast, ArdB-like proteins inhibit RMI systems independently of the DNA recognition site.

Complete genome assembly and methylome dissection of Methanococcus aeolicus PL15/H p.

Although restriction-modification systems are found in both Eubacterial and Archaeal kingdoms, comparatively less is known about patterns of DNA methylation and genome defense systems in archaea. Here we report the complete closed genome sequence and methylome analysis of Methanococcus aeolicus PL15/H p , a strain of the CO2-reducing methanogenic archaeon and a commercial source for MaeI, MaeII, and MaeIII restriction endonucleases. The M. aeolicus PL15/H p genome consists of a 1.68 megabase circular chromosome predicted to contain 1,615 protein coding genes and 38 tRNAs. A combination of methylome sequencing, homology-based genome annotation, and recombinant gene expression identified five restriction-modification systems encoded by this organism, including the methyltransferase and site-specific endonuclease of MaeIII. The MaeIII restriction endonuclease was recombinantly expressed, purified and shown to have site-specific DNA cleavage activity in vitro.

7-Deazaguanines in DNA: functional and structural elucidation of a DNA modification system.

The modified nucleosides 2'-deoxy-7-cyano- and 2'-deoxy-7-amido-7-deazaguanosine (dPreQ0 and dADG, respectively) recently discovered in DNA are the products of the bacterial queuosine tRNA modification pathway and the dpd gene cluster, the latter of which encodes proteins that comprise the elaborate Dpd restriction-modification system present in diverse bacteria. Recent genetic studies implicated the dpdA, dpdB and dpdC genes as encoding proteins necessary for DNA modification, with dpdD-dpdK contributing to the restriction phenotype. Here we report the in vitro reconstitution of the Dpd modification machinery from Salmonella enterica serovar Montevideo, the elucidation of the roles of each protein and the X-ray crystal structure of DpdA supported by small-angle X-ray scattering analysis of DpdA and DpdB, the former bound to DNA. While the homology of DpdA with the tRNA-dependent tRNA-guanine transglycosylase enzymes (TGT) in the queuosine pathway suggested a similar transglycosylase activity responsible for the exchange of a guanine base in the DNA for 7-cyano-7-deazaguanine (preQ0), we demonstrate an unexpected ATPase activity in DpdB necessary for insertion of preQ0 into DNA, and identify several catalytically essential active site residues in DpdA involved in the transglycosylation reaction. Further, we identify a modification site for DpdA activity and demonstrate that DpdC functions independently of DpdA/B in converting preQ0-modified DNA to ADG-modified DNA.

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.

Investigation of a Listeria monocytogenes Chromosomal Immigration Control Region Reveals Diverse Restriction Modification Systems with Complete Sequence Type Conservation.

Listeria monocytogenes is a Gram-positive pathogen responsible for the severe foodborne disease listeriosis. A chromosomal hotspot between lmo0301 and lmo0305 has been noted to harbor diverse restriction modification (RM) systems. Here, we analyzed 872 L. monocytogenes genomes to better understand the prevalence and types of RM systems in this region, designated the immigration control region (ICR). Type I, II, III and IV RM systems were found in 86.1% of strains inside the ICR and in 22.5% of strains flanking the ICR. ICR content was completely conserved within the same multilocus sequence typing-based sequence type (ST), but the same RM system could be identified in diverse STs. The intra-ST conservation of ICR content suggests that this region may drive the emergence of new STs and promote clone stability. Sau3AI-like, LmoJ2 and LmoJ3 type II RM systems as well as type I EcoKI-like, and type IV AspBHI-like and mcrB-like systems accounted for all RM systems in the ICR. A Sau3AI-like type II RM system with specificity for GATC was harbored in the ICR of many STs, including all strains of the ancient, ubiquitous ST1. The extreme paucity of GATC recognition sites in lytic phages may reflect ancient adaptation of these phages to preempt resistance associated with the widely distributed Sau3AI-like systems. These findings indicate that the ICR has a high propensity for RM systems which are intraclonaly conserved and may impact bacteriophage susceptibility as well as ST emergence and stability.

Characterization of the Staphylococcus xylosus methylome reveals a new variant of type I restriction modification system in staphylococci.

Restriction modification (RM) systems are known to provide a strong barrier to the exchange of DNA between and within bacterial species. Likewise, DNA methylation is known to have an important function in bacterial epigenetics regulating essential pathways such as DNA replication and the phase variable expression of prokaryotic phenotypes. To date, research on staphylococcal DNA methylation focused mainly on the two species Staphylococcus aureus and S. epidermidis. Less is known about other members of the genus such as S. xylosus, a coagulase-negative commensal of mammalian skin. The species is commonly used as starter organism in food fermentations but is also increasingly considered to have an as yet elusive function in bovine mastitis infections. We analyzed the methylomes of 14 S. xylosus strains using single-molecular, real-time (SMRT) sequencing. Subsequent in silico sequence analysis allowed identification of the RM systems and assignment of the respective enzymes to the discovered modification patterns. Hereby the presence of type I, II, III and IV RM systems in varying numbers and combinations among the different strains was revealed, clearly distinguishing the species from what is known for other members of the genus so far. In addition, the study characterizes a newly discovered type I RM system, encoded by S. xylosus but also by a variety of other staphylococcal species, with a hitherto unknown gene arrangement that involves two specificity units instead of one (hsdRSMS). Expression of different versions of the operon in E. coli showed proper base modification only when genes encoding both hsdS subunits were present. This study provides new insights into the general understanding of the versatility and function of RM systems as well as the distribution and variations in the genus Staphylococcus.

Proteomic Study of the Interactions between Phages and the Bacterial Host Klebsiella pneumoniae.

Phages and bacteria have acquired resistance mechanisms for protection. In this context, the aims of the present study were to analyze the proteins isolated from 21 novel lytic phages of Klebsiella pneumoniae in search of defense mechanisms against bacteria and also to determine the infective capacity of the phages. A proteomic study was also conducted to investigate the defense mechanisms of two clinical isolates of K. pneumoniae infected by phages. For this purpose, the 21 lytic phages were sequenced and de novo assembled. The host range was determined in a collection of 47 clinical isolates of K. pneumoniae, revealing the variable infective capacity of the phages. Genome sequencing showed that all of the phages were lytic phages belonging to the order Caudovirales. Phage sequence analysis revealed that the proteins were organized in functional modules within the genome. Although most of the proteins have unknown functions, multiple proteins were associated with defense mechanisms against bacteria, including the restriction-modification system, the toxin-antitoxin system, evasion of DNA degradation, blocking of host restriction and modification, the orphan CRISPR-Cas system, and the anti-CRISPR system. Proteomic study of the phage-host interactions (i.e., between isolates K3574 and K3320, which have intact CRISPR-Cas systems, and phages vB_KpnS-VAC35 and vB_KpnM-VAC36, respectively) revealed the presence of several defense mechanisms against phage infection (prophage, defense/virulence/resistance, oxidative stress and plasmid proteins) in the bacteria, and of the Acr candidate (anti-CRISPR protein) in the phages. IMPORTANCE Researchers, including microbiologists and infectious disease specialists, require more knowledge about the interactions between phages and their bacterial hosts and about their defense mechanisms. In this study, we analyzed the molecular mechanisms of viral and bacterial defense in phages infecting clinical isolates of K. pneumoniae. Viral defense mechanisms included restriction-modification system evasion, the toxin-antitoxin (TA) system, DNA degradation evasion, blocking of host restriction and modification, and resistance to the abortive infection system, anti-CRISPR and CRISPR-Cas systems. Regarding bacterial defense mechanisms, proteomic analysis revealed expression of proteins involved in the prophage (FtsH protease modulator), plasmid (cupin phosphomannose isomerase protein), defense/virulence/resistance (porins, efflux pumps, lipopolysaccharide, pilus elements, quorum network proteins, TA systems, and methyltransferases), oxidative stress mechanisms, and Acr candidates (anti-CRISPR protein). The findings reveal some important molecular mechanisms involved in the phage-host bacterial interactions; however, further study in this field is required to improve the efficacy of phage therapy.

Endogenous Linear Plasmids lp28-4 and lp25 Are Required for Infectivity and Restriction Protection in the Lyme Disease Spirochete Borrelia mayonii.

Borrelia mayonii is a newly recognized causative agent of Lyme disease in the Upper Midwestern United States, with distinct clinical presentations compared to classical Lyme disease caused by other Lyme Borrelia species. However, little is known about the B. mayonii genetic determinants required for establishing infection or perpetuating disease in mammals. Extrachromosomal plasmids in Borrelia species often encode proteins necessary for infection and pathogenesis, and spontaneous loss of these plasmids can lead to the identification of virulence determinant genes. Here, we describe infection of Lyme disease-susceptible C3H mice with B. mayonii, and show bacterial dissemination and persistence in peripheral tissues. Loss of endogenous plasmids, including lp28-4, lp25, and lp36 correlated with reduced infectivity in mice. The apparent requirement for lp28-4 during murine infection suggests the presence of a novel virulence determinant, as this plasmid does not encode homologs of any known virulence determinant. We also describe transformation and stable maintenance of a self-replicating shuttle vector in B. mayonii, and show that loss of either lp25 or lp28-4 correlated with increased transformation competency. Finally, we demonstrate that linear plasmids lp25 and lp28-4 each encode functional restriction modification systems with distinct but partially overlapping target modification sequences, which likely accounts for the observed decrease in transformation efficiency when those plasmids are present. Taken together, this study describes a role for endogenous plasmids in mammalian infection and restriction protection in the Lyme disease spirochete Borrelia mayonii.

A Grid Search-Based Multilayer Dynamic Ensemble System to Identify DNA N4-Methylcytosine Using Deep Learning Approach.

DNA (Deoxyribonucleic Acid) N4-methylcytosine (4mC), a kind of epigenetic modification of DNA, is important for modifying gene functions, such as protein interactions, conformation, and stability in DNA, as well as for the control of gene expression throughout cell development and genomic imprinting. This simply plays a crucial role in the restriction-modification system. To further understand the function and regulation mechanism of 4mC, it is essential to precisely locate the 4mC site and detect its chromosomal distribution. This research aims to design an efficient and high-throughput discriminative intelligent computational system using the natural language processing method "word2vec" and a multi-configured 1D convolution neural network (1D CNN) to predict 4mC sites. In this article, we propose a grid search-based multi-layer dynamic ensemble system (GS-MLDS) that can enhance existing knowledge of each level. Each layer uses a grid search-based weight searching approach to find the optimal accuracy while minimizing computation time and additional layers. We have used eight publicly available benchmark datasets collected from different sources to test the proposed model's efficiency. Accuracy results in test operations were obtained as follows: 0.978, 0.954, 0.944, 0.961, 0.950, 0.973, 0.948, 0.952, 0.961, and 0.980. The proposed model has also been compared to 16 distinct models, indicating that it can accurately predict 4mC.

Evolution of restriction-modification systems with one restriction endonuclease and two DNA methyltransferases

Some restriction-modification systems contain two DNA methyltransferases. In the present work, we have classified such systems according to the families of catalytic domains present in restriction endonucleases and both DNA methyltransferases. The evolution of restriction-modification systems of one class was studied in detail. Systems in this class include an endonuclease with a NOV_C family domain and two DNA methyltransferases, both with DNA_methylase family domains. The phylogenetic tree of DNA methyltransferases from systems of this class consists of two clades of the same size. Two DNA methyltransferases of each restriction-modification system of the class belong to different clades. This indicates independent evolution of the two methyltransferases. We detected multiple cross-species horizontal transfers of systems as a whole, as well as cases of gene transfer between systems.

Restriction endonuclease cleavage of phage DNA enables resuscitation from Cas13-induced bacterial dormancy.

Type VI CRISPR systems protect against phage infection using the RNA-guided nuclease Cas13 to recognize viral messenger RNA. Upon target recognition, Cas13 cleaves phage and host transcripts non-specifically, leading to cell dormancy that is incompatible with phage propagation. However, whether and how infected cells recover from dormancy is unclear. Here we show that type VI CRISPR and DNA-cleaving restriction-modification (RM) systems frequently co-occur and synergize to clear phage infections and resuscitate cells. In the natural type VI CRISPR host Listeria seeligeri, we show that RM cleaves the phage genome, thus removing the source of phage transcripts and enabling cells to recover from Cas13-induced cellular dormancy. We find that phage infections are neutralized more effectively when Cas13 and RM systems operate together. Our work reveals that type VI CRISPR immunity is cell-autonomous and non-abortive when paired with RM, and hints at other synergistic roles for the diverse host-directed immune systems in bacteria.

Evolution of Restriction–Modification Systems Consisting of One Restriction Endonuclease and Two DNA Methyltransferases

Some restriction-modification systems contain two DNA methyltransferases. In the present work, we have classified such systems according to the families of catalytic domains present in the restriction endonucleases and both DNA methyltransferases. Evolution of the restriction-modification systems containing an endonuclease with a NOV_C family domain and two DNA methyltransferases, both with DNA_methylase family domains, was investigated in detail. Phylogenetic tree of DNA methyltransferases from the systems of this class consists of two clades of the same size. Two DNA methyltransferases of each restriction-modification system of this class belong to the different clades. This indicates independent evolution of the two methyltransferases. We detected multiple cross-species horizontal transfers of the systems as a whole, as well as the cases of gene transfer between the systems.

Exploring a High-Efficiency Genetic Transformation System for Engineering Saccharopolyspora pogona ASAGF58 To Improve Butenyl-Spinosyn Production

Butenyl-spinosyn is a potent insecticide potentially useful as a broad-spectrum pesticide. Because it is relatively nontoxic to mammals and does not damage the environment, there has been considerable interest in the use of butenyl-spinosyn for increasing agricultural production. However, genetically engineering Saccharopolyspora pogona ASAGF58 to increase its relatively low butenyl-spinosyn content remains challenging because it cannot be transformed efficiently. In this study, genes encoding novel methyltransferases (04455 and 28970) were identified in the Sa. pogona ASAGF58 genome through a bioinformatic-based analysis. Additionally, the transformation efficiency increased by 5.8- and 16.4-fold when foreign DNA was pre-methylated in ET-28970 and ET-04455, respectively, through bypassing of the restriction–modification system. A comparative proteomic analysis of Sa. pogona and Saccharopolyspora spinosa revealed that acetyl-CoA synthetase may be useful for improving butenyl-spinosyn production. The fermentation results indicated that compared with the wild-type butenyl-spinosyn content, overexpressing the acetyl-CoA synthetase gene increased butenyl-spinosyn production by 2-fold. The findings presented herein suggest that the strategy employed in this study may be applicable for the genetic engineering of other nonmodel Saccharopolyspora strains and increase target product yields.

The Evolutionary History of a DNA Methylase Reveals Frequent Horizontal Transfer and Within-Gene Recombination.

Inteins, often referred to as protein introns, are highly mobile genetic elements that invade conserved genes throughout the tree of life. Inteins have been found to invade a wide variety of key genes within actinophages. While in the process of conducting a survey of these inteins in actinophages, we discovered that one protein family of methylases contained a putative intein, and two other unique insertion elements. These methylases are known to occur commonly in phages as orphan methylases (possibly as a form of resistance to restriction-modification systems). We found that the methylase family is not conserved within phage clusters and has a disparate distribution across divergent phage groups. We determined that two of the three insertion elements have a patchy distribution within the methylase protein family. Additionally, we found that the third insertion element is likely a second homing endonuclease, and that all three elements (the intein, the homing endonuclease, and what we refer to as the ShiLan domain) have different insertion sites that are conserved in the methylase gene family. Furthermore, we find strong evidence that both the intein and ShiLan domain are partaking in long-distance horizontal gene transfer events between divergent methylases in disparate phage hosts within the already dispersed methylase distribution. The reticulate evolutionary history of methylases and their insertion elements reveals high rates of gene transfer and within-gene recombination in actinophages.

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.

Systematic identification of gene-altering programmed inversions across the bacterial domain.

Programmed chromosomal inversions allow bacteria to generate intra-population genotypic and functional heterogeneity, a bet-hedging strategy important in changing environments. Some programmed inversions modify coding sequences, producing different alleles in several gene families, most notably in specificity-determining genes such as Type I restriction-modification systems, where systematic searches revealed cross phylum abundance. Yet, a broad, gene-independent, systematic search for gene-altering programmed inversions has been absent, and little is known about their genomic sequence attributes and prevalence across gene families. Here, identifying intra-species variation in genomes of over 35 000 species, we develop a predictive model of gene-altering inversions, revealing key attributes of their genomic sequence attributes, including gene-pseudogene size asymmetry and orientation bias. The model predicted over 11,000 gene-altering loci covering known targeted gene families, as well as novel targeted families including Type II restriction-modification systems, a protein of unknown function, and a fusion-protein containing conjugative-pilus and phage tail domains. Publicly available long-read sequencing datasets validated representatives of these newly predicted inversion-targeted gene families, confirming intra-population genetic heterogeneity. Together, these results reveal gene-altering programmed inversions as a key strategy adopted across the bacterial domain, and highlight programmed inversions that modify Type II restriction-modification systems as a possible new mechanism for maintaining intra-population heterogeneity.

Natural Recombination among Type I Restriction-Modification Systems Creates Diverse Genomic Methylation Patterns among Xylella fastidiosa Strains.

Xylella fastidiosa is an important bacterial plant pathogen causing high-consequence diseases in agricultural crops around the world. Although as a species X. fastidiosa can infect many host plants, there is significant variability between strains regarding virulence on specific host plant species and other traits. Natural competence and horizontal gene transfer are believed to occur frequently in X. fastidiosa and likely influence the evolution of this pathogen. However, some X. fastidiosa strains are difficult to manipulate genetically using standard transformation techniques. Several type I restriction-modification (R-M) systems are encoded in the X. fastidiosa genome, which may influence horizontal gene transfer and recombination. Type I R-M systems themselves may undergo recombination, exchanging target recognition domains (TRDs) between specificity subunits (hsdS) to generate novel alleles with new target specificities. In this study, several conserved type I R-M systems were compared across 129 X. fastidiosa genome assemblies representing all known subspecies and 32 sequence types. Forty-four unique TRDs were identified among 50 hsdS alleles, which are arrayed in 31 allele profiles that are generally conserved within a monophyletic cluster of strains. Inactivating mutations were identified in type I R-M systems of specific strains, showing heterogeneity in the complements of functional type I R-M systems across X. fastidiosa. Genomic DNA methylation patterns were characterized in 20 X. fastidiosa strains and associated with type I R-M system allele profiles. Overall, these data suggest hsdS genes recombine among Xylella strains and/or unknown donors, and the resulting TRD reassortment establishes differential epigenetic modifications across Xylella lineages. IMPORTANCE Economic impacts on agricultural production due to X. fastidiosa have been severe in the Americas, Europe, and parts of Asia. Despite a long history of research on this pathogen, certain fundamental questions regarding the biology, pathogenicity, and evolution of X. fastidiosa have still not been answered. Wide-scale whole-genome sequencing has begun to provide more insight into X. fastidiosa genetic diversity and horizontal gene transfer, but the mechanics of genomic recombination in natural settings and the extent to which this directly influences bacterial phenotypes such as plant host range are not well understood. Genome methylation is an important factor in horizontal gene transfer and bacterial recombination that has not been comprehensively studied in X. fastidiosa. This study characterizes methylation associated with type I restriction-modification systems across a wide range of X. fastidiosa strains and lays the groundwork for a better understanding of X. fastidiosa biology and evolution through epigenetics.