Type I restriction systems: sophisticated molecular machines (a legacy of Bertani and Weigle)

Evaluation of the enzymatic properties of DNA (cytosine-5)-methyltransferase M.ApeKI from archaea in the presence of metal ions

We previously identified M.ApeKI from Aeropyum pernix K1 as a highly thermostable DNA (cytosine-5)-methyltransferase. M.ApeKI uses the type II restriction-modification system (R-M system), among the best-studied R-M systems. Although endonucleases generally utilize Mg (II) as a cofactor, several reports have shown that MTases exhibit different reactions in the presence of metal ions. This study aim was to evaluate the enzymatic properties of DNA (cytosine-5)-methyltransferase M.ApeKI from archaea in the presence of metal ions. We evaluated the influence of metal ions on the catalytic activity and DNA binding of M.ApeKI. The catalytic activity was inhibited by Cu (II), Mg (II), Mn (II), and Zn (II), each at 5 m m. DNA binding was more strongly inhibited by 5 m m Cu (II) and 10 m m Zn (II). To our knowledge, this is the first report showing that DNA binding of type II MTase is inhibited by metal ions.

A Vibrio cholerae Type IV restriction system targets glucosylated 5-hydroxymethylcytosine to protect against phage infection

A major challenge faced by Vibrio cholerae is constant predation by bacteriophage (phage) in aquatic reservoirs and during infection of human hosts. To overcome phage predation, V. cholerae has acquired and/or evolved a myriad of phage defense systems. Although several novel defense systems have been discovered, we hypothesized that more were encoded in V. cholerae given the low diversity of phages that have been isolated, which infect this species. Using a V. cholerae genomic library, we identified a Type IV restriction system consisting of two genes within a 16-kB region of the Vibrio pathogenicity island-2, which we name TgvA and TgvB (Type I-embedded gmrSD-like system of VPI-2). We show that both TgvA and TgvB are required for defense against T2, T4, and T6 by targeting glucosylated 5-hydroxymethylcytosine (5hmC). T2 or T4 phages that lose the glucose modifications are resistant to TgvAB defense but exhibit a significant evolutionary tradeoff, becoming susceptible to other Type IV restriction systems that target unglucosylated 5hmC. We also show that the Type I restriction-modification system that embeds the tgvAB genes protects against phage T3, secΦ18, secΦ27, and λ, suggesting that this region is a phage defense island. Our study uncovers a novel Type IV restriction system in V. cholerae, increasing our understanding of the evolution and ecology of V. cholerae, while highlighting the evolutionary interplay between restriction systems and phage genome modification.IMPORTANCEBacteria are constantly being predated by bacteriophage (phage). To counteract this predation, bacteria have evolved a myriad of defense systems. Some of these systems specifically digest infecting phage by recognizing unique base modifications present on the phage DNA. In this study, we discover a Type IV restriction system encoded in V. cholerae, which we name TgvAB, and demonstrate it recognizes and restricts phage that have 5-hydroxymethylcytosine glucosylated DNA. Moreover, the evolution of resistance to TgvAB render phage susceptible to other Type IV restriction systems, demonstrating a significant evolutionary tradeoff. These results enhance our understanding of the evolution of V. cholerae and more broadly how bacteria evade phage predation.

Vibrio cholerae pathogenicity island 2 encodes two distinct types of restriction systems

In response to predation by bacteriophages and invasion by other mobile genetic elements such as plasmids, bacteria have evolved specialized defense systems that are often clustered together on genomic islands. The O1 El Tor strains of Vibrio cholerae responsible for the ongoing seventh cholera pandemic (7PET) contain a characteristic set of genomic islands involved in host colonization and disease, many of which contain defense systems. Notably, Vibrio pathogenicity island 2 contains several characterized defense systems as well as a putative type I restriction-modification (T1RM) system, which, interestingly, is interrupted by two genes of unknown function. Here, we demonstrate that the T1RM system is active, methylates the host genomes of a representative set of 7PET strains, and identify a specific recognition sequence that targets non-methylated plasmids for restriction. We go on to show that the two genes embedded within the T1RM system encode a novel two-protein modification-dependent restriction system related to the GmrSD family of type IV restriction enzymes. Indeed, we show that this system has potent anti-phage activity against diverse members of the Tevenvirinae, a subfamily of bacteriophages with hypermodified genomes. Taken together, these results expand our understanding of how this highly conserved genomic island contributes to the defense of pandemic V. cholerae against foreign DNA.

IMPORTANCE

Defense systems are immunity systems that allow bacteria to counter the threat posed by bacteriophages and other mobile genetic elements. Although these systems are numerous and highly diverse, the most common types are restriction enzymes that can specifically recognize and degrade non-self DNA. Here, we show that the Vibrio pathogenicity island 2, present in the pathogen Vibrio cholerae, encodes two types of restriction systems that use distinct mechanisms to sense non-self DNA. The first system is a classical Type I restriction-modification system, and the second is a novel modification-dependent type IV restriction system that recognizes hypermodified cytosines. Interestingly, these systems are embedded within each other, suggesting that they are complementary to each other by targeting both modified and non-modified phages.

Identifying functional multi-host shuttle plasmids to advance synthetic biology applications in Mesorhizobium and Bradyrhizobium

Ammonia availability has a crucial role in agriculture as it ensures healthy plant growth and increased crop yields. Since diazotrophs are the only organisms capable of reducing dinitrogen to ammonia, they have great ecological importance and potential to mitigate the environmental and economic costs of synthetic fertilizer use. Rhizobia are especially valuable being that they can engage in nitrogen-fixing symbiotic relationships with legumes, and they demonstrate great diversity and plasticity in genomic and phenotypic traits. However, few rhizobial species have sufficient genetic tractability for synthetic biology applications. This study established a basic genetic toolbox with antibiotic resistance markers, multi-host shuttle plasmids and a streamlined protocol for biparental conjugation with Mesorhizobium and Bradyrhizobium species. We identified two repABC origins of replication from Sinorhizobium meliloti (pSymB) and Rhizobium etli (p42d) that were stable across all three strains of interest. Furthermore, the NZP2235 genome was sequenced and phylogenetic analysis determined its reclassification to Mesorhizobium huakuii. These tools will enable the use of plasmid-based strategies for more advanced genetic engineering projects and ultimately contribute towards the development of more sustainable agriculture practices by means of novel nitrogen-fixing organelles, elite bioinoculants, or symbiotic association with nonlegumes.

CRISPR-Cas3 and type I restriction-modification team up against blaKPC-IncF plasmid transfer in Klebsiella pneumoniae

OBJECTIVE

We explored whether the Clustered regularly interspaced short palindromic repeat (CRISPR)-Cas and restriction-modification (R-M) systems are compatible and act together to resist plasmid attacks.

METHODS

932 global whole-genome sequences from GenBank, and 459 K. pneumoniae isolates from six provinces of China, were collected to investigate the co-distribution of CRISPR-Cas, R-M systems, and blaKPC plasmid. Conjugation and transformation assays were applied to explore the anti-plasmid function of CRISPR and R-M systems.

RESULTS

We found a significant inverse correlation between the presence of CRISPR and R-M systems and blaKPC plasmids in K. pneumoniae, especially when both systems cohabited in one host. The multiple matched recognition sequences of both systems in blaKPC-IncF plasmids (97%) revealed that they were good targets for both systems. Furthermore, the results of conjugation assay demonstrated that CRISPR-Cas and R-M systems in K. pneumoniae could effectively hinder blaKPC plasmid invasion. Notably, CRISPR-Cas and R-M worked together to confer a 4-log reduction in the acquisition of blaKPC plasmid in conjugative events, exhibiting robust synergistic anti-plasmid immunity.

CONCLUSIONS

Our results indicate the synergistic role of CRISPR and R-M in regulating horizontal gene transfer in K. pneumoniae and rationalize the development of antimicrobial strategies that capitalize on the immunocompromised status of KPC-KP.

RecA-dependent or independent recombination of plasmid DNA generates a conflict with the host EcoKI immunity by launching restriction alleviation

Bacterial defence systems are tightly regulated to avoid autoimmunity. In Type I restriction-modification (R-M) systems, a specific mechanism called restriction alleviation (RA) controls the activity of the restriction module. In the case of the Escherichia coli Type I R-M system EcoKI, RA proceeds through ClpXP-mediated proteolysis of restriction complexes bound to non-methylated sites that appear after replication or reparation of host DNA. Here, we show that RA is also induced in the presence of plasmids carrying EcoKI recognition sites, a phenomenon we refer to as plasmid-induced RA. Further, we show that the anti-restriction behavior of plasmid-borne non-conjugative transposons such as Tn5053, previously attributed to their ardD loci, is due to plasmid-induced RA. Plasmids carrying both EcoKI and Chi sites induce RA in RecA- and RecBCD-dependent manner. However, inactivation of both RecA and RecBCD restores RA, indicating that there exists an alternative, RecA-independent, homologous recombination pathway that is blocked in the presence of RecBCD. Indeed, plasmid-induced RA in a RecBCD-deficient background does not depend on the presence of Chi sites. We propose that processing of random dsDNA breaks in plasmid DNA via homologous recombination generates non-methylated EcoKI sites, which attract EcoKI restriction complexes channeling them for ClpXP-mediated proteolysis.

A Vibrio cholerae Type IV restriction system targets glucosylated 5-hydroxyl methyl cytosine to protect against phage infection

A major challenge faced by Vibrio cholerae is constant predation by bacteriophage (phage) in aquatic reservoirs and during infection of human hosts. To overcome phage predation, V. cholerae has evolved a myriad of phage defense systems. Although several novel defense systems have been discovered, we hypothesized more were encoded in V. cholerae given the relative paucity of phage that have been isolated which infect this species. Using a V. cholerae genomic library, we identified a Type IV restriction system consisting of two genes within a 16kB region of the Vibrio pathogenicity island-2 that we name TgvA and TgvB (Type I-embedded gmrSD-like system of VPI-2). We show that both TgvA and TgvB are required for defense against T2, T4, and T6 by targeting glucosylated 5-hydroxymethylcytosine (5hmC). T2 or T4 phages that lose the glucose modification are resistant to TgvAB defense but exhibit a significant evolutionary tradeoff becoming susceptible to other Type IV restriction systems that target unglucosylated 5hmC. We show that additional phage defense genes are encoded in VPI-2 that protect against other phage like T3, secΦ18, secΦ27 and λ. Our study uncovers a novel Type IV restriction system in V. cholerae, increasing our understanding of the evolution and ecology of V. cholerae while highlighting the evolutionary interplay between restriction systems and phage genome modification.

Recovering high-quality bacterial genomes from cross-contaminated cultures: a case study of marine Vibrio campbellii

BACKGROUND

Environmental monitoring of bacterial pathogens is critical for disease control in coastal marine ecosystems to maintain animal welfare and ecosystem function and to prevent significant economic losses. This requires accurate taxonomic identification of environmental bacterial pathogens, which often cannot be achieved by commonly used genetic markers (e.g., 16S rRNA gene), and an understanding of their pathogenic potential based on the information encoded in their genomes. The decreasing costs of whole genome sequencing (WGS), combined with newly developed bioinformatics tools, now make it possible to unravel the full potential of environmental pathogens, beyond traditional microbiological approaches. However, obtaining a high-quality bacterial genome, requires initial cultivation in an axenic culture, which is a bottleneck in environmental microbiology due to cross-contamination in the laboratory or isolation of non-axenic strains.

RESULTS

We applied WGS to determine the pathogenic potential of two Vibrio isolates from coastal seawater. During the analysis, we identified cross-contamination of one of the isolates and decided to use this dataset to evaluate the possibility of bioinformatic contaminant removal and recovery of bacterial genomes from a contaminated culture. Despite the contamination, using an appropriate bioinformatics workflow, we were able to obtain high quality and highly identical genomes (Average Nucleotide Identity value 99.98%) of one of the Vibrio isolates from both the axenic and the contaminated culture. Using the assembled genome, we were able to determine that this isolate belongs to a sub-lineage of Vibrio campbellii associated with several diseases in marine organisms. We also found that the genome of the isolate contains a novel Vibrio plasmid associated with bacterial defense mechanisms and horizontal gene transfer, which may offer a competitive advantage to this putative pathogen.

CONCLUSIONS

Our study shows that, using state-of-the-art bioinformatics tools and a sufficient sequencing effort, it is possible to obtain high quality genomes of the bacteria of interest and perform in-depth genomic analyses even in the case of a contaminated culture. With the new isolate and its complete genome, we are providing new insights into the genomic characteristics and functional potential of this sub-lineage of V. campbellii. The approach described here also highlights the possibility of recovering complete bacterial genomes in the case of non-axenic cultures or obligatory co-cultures.

Methicillin Resistance Elements in the Canine Pathogen Staphylococcus pseudintermedius and Their Association with the Peptide Toxin PSM-mec

Staphylococcus pseudintermedius is a frequent cause of infections in dogs. Infectious isolates of this coagulase-positive staphylococcal species are often methicillin- and multidrug-resistant, which complicates therapy. In staphylococci, methicillin resistance is encoded by determinants found on mobile genetic elements called Staphylococcal Chromosome Cassette mec (SCCmec), which, in addition to methicillin resistance factors, sometimes encode additional genes, such as further resistance factors and, rarely, virulence determinants. In this study, we analyzed SCCmec in a collection of infectious methicillin-resistant S. pseudintermedius (MRSP) isolates from predominant lineages in the United States. We found that several lineages characteristically have specific types of SCCmec elements and Agr types and harbor additional factors in their SCCmec elements that may promote virulence or affect DNA uptake. All isolates had SCCmec-encoded restriction-modification (R-M) systems of types I or II, and sequence types (STs) ST84 and ST64 had one type II and one type I R-M system, although the latter lacked a complete methylation enzyme gene. ST68 isolates also had an SCCmec-encoded CRISPR system. ST71 isolates had a psm-mec gene, which, in all but apparently Agr-dysfunctional isolates, produced a PSM-mec peptide toxin, albeit at relatively small amounts. This study gives detailed insight into the composition of SCCmec elements in infectious isolates of S. pseudintermedius and lays the genetic foundation for further efforts directed at elucidating the contribution of identified accessory SCCmec factors in impacting SCCmec-encoded and thus methicillin resistance-associated virulence and resistance to DNA uptake in this leading canine pathogen.

The Restriction-Modification Systems of Clostridium carboxidivorans P7

Clostridium carboxidivorans P7 (DSM 15243) is a bacterium that converts syngas (a mixture of CO, H2, and CO2) into hexanol. An optimized and scaled-up industrial process could therefore provide a renewable source of fuels and chemicals while consuming industry waste gases. However, the genetic engineering of this bacterium is hindered by its multiple restriction-modification (RM) systems: the genome of C. carboxidivorans encodes at least ten restriction enzymes and eight methyltransferases (MTases). To gain insight into the complex RM systems of C. carboxidivorans, we analyzed genomic methylation patterns using single-molecule real-time (SMRT) sequencing and bisulfite sequencing. We identified six methylated sequence motifs. To match the methylation sites to the predicted MTases of C. carboxidivorans, we expressed them individually in Escherichia coli for functional characterization. Recognition motifs were identified for all three Type I MTases (CAYNNNNNCTGC/GCAGNNNNNRTG, CCANNNNNNNNTCG/CGANNNNNNNNTGG and GCANNNNNNNTNNCG/CGNNANNNNNNNTGC), two Type II MTases (GATAAT and CRAAAAR), and a single Type III MTase (GAAAT). However, no methylated recognition motif was found for one of the three Type II enzymes. One recognition motif that was methylated in C. carboxidivorans but not in E. coli (AGAAGC) was matched to the remaining Type III MTase through a process of elimination. Understanding these enzymes and the corresponding recognition sites will facilitate the development of genetic tools for C. carboxidivorans that can accelerate the industrial exploitation of this strain.

T5-like phage BF23 evades host-mediated DNA restriction and methylation

Bacteriophage BF23 is a close relative of phage T5, a prototypical Tequintavirus that infects Escherichia coli. BF23 was isolated in the middle of the XXth century and was extensively studied as a model object. Like T5, BF23 carries long ∼9.7 kb terminal repeats, injects its genome into infected cell in a two-stage process, and carries multiple specific nicks in its double-stranded genomic DNA. The two phages rely on different host secondary receptors-FhuA (T5) and BtuB (BF23). Only short fragments of the BF23 genome, including the region encoding receptor interacting proteins, have been determined. Here, we report the full genomic sequence of BF23 and describe the protein content of its virion. T5-like phages represent a unique group that resist restriction by most nuclease-based host immunity systems. We show that BF23, like other Tequintavirus phages, resist Types I/II/III restriction-modification host immunity systems if their recognition sites are located outside the terminal repeats. We also demonstrate that the BF23 avoids host-mediated methylation. We propose that inhibition of methylation is a common feature of Tequintavirus and Epseptimavirus genera phages, that is not, however, associated with their antirestriction activity.

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.

Acquisition of Type I methyltransferase via horizontal gene transfer increases the drug resistance of Aeromonas veronii

Aeromonas veronii is an opportunistic pathogen that affects both fish and mammals, including humans, leading to bacteraemia, sepsis, meningitis and even death. The increasing virulence and drug resistance of A. veronii are of significant concern and pose a severe risk to public safety. The Type I restriction-modification (RM) system, which functions as a bacterial defence mechanism, can influence gene expression through DNA methylation. However, little research has been conducted to explore its origin, evolutionary path, and relationship to virulence and drug resistance in A. veronii. In this study, we analysed the pan-genome of 233 A. veronii strains, and the results indicated that it was 'open', meaning that A. veronii has acquired additional genes from other species. This suggested that A. veronii had the potential to adapt and evolve rapidly, which might have contributed to its drug resistance. One Type I methyltransferase (MTase) and two complete Type I RM systems were identified, namely AveC4I, AveC4II and AveC4III in A. veronii strain C4, respectively. Notably, AveC4I was exclusive to A. veronii C4. Phylogenetic analysis revealed that AveC4I was derived from horizontal gene transfer from Thiocystis violascens and exchanged genes with the human pathogen Comamonas kerstersii. Single molecule real-time sequencing was applied to identify the motif methylated by AveC4I, which was unique and not recognized by any reported MTases in the REBASE database. We also annotated the functions and pathways of the genes containing the motif, revealing that AveC4I may control drug resistance in A. veronii C4. Our findings provide new insight on the mechanisms underlying drug resistance in pathogenic bacteria. By identifying the specific genes and pathways affected by AveC4I, this study may aid in the development of new therapeutic approaches to combat A. veronii infections.

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.

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 Sequence and Analysis of a ST573 Multidrug-Resistant Methicillin-Resistant Staphylococcus aureus SauR3 Clinical Isolate from Terengganu, Malaysia

Methicillin-resistant Staphylococcus aureus (MRSA) is a World Health Organization-listed priority pathogen. Scarce genomic data are available for MRSA isolates from Malaysia. Here, we present the complete genome sequence of a multidrug-resistant MRSA strain SauR3, isolated from the blood of a 6-year-old patient hospitalized in Terengganu, Malaysia, in 2016. S. aureus SauR3 was resistant to five antimicrobial classes comprising nine antibiotics. The genome was sequenced on the Illumina and Oxford Nanopore platforms and hybrid assembly was performed to obtain its complete genome sequence. The SauR3 genome consists of a circular chromosome of 2,800,017 bp and three plasmids designated pSauR3-1 (42,928 bp), pSauR3-2 (3011 bp), and pSauR3-3 (2473 bp). SauR3 belongs to sequence type 573 (ST573), a rarely reported sequence type of the staphylococcal clonal complex 1 (CC1) lineage, and harbors a variant of the staphylococcal cassette chromosome mec (SCCmec) type V (5C2&5) element which also contains the aac(6')-aph(2″) aminoglycoside-resistance genes. pSauR3-1 harbors several antibiotic resistance genes in a 14,095 bp genomic island (GI), previously reported in the chromosome of other staphylococci. pSauR3-2 is cryptic, whereas pSauR3-3 encodes the ermC gene that mediates inducible resistance to macrolide-lincosamide-streptogramin B (iMLS_B). The SauR3 genome can potentially be used as a reference genome for other ST573 isolates.

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.

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.

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.

An Emerging Lineage of Uropathogenic Extended Spectrum β-Lactamase Escherichia coli ST127

Uropathogenic Escherichia coli (UPEC) is one of the most common causes of urinary tract infections. Here, we report for the first time the whole-genome sequencing (WGS) and analysis of four extended-spectrum β-lactamase (ESBL), UPEC sequence type (ST) 127 isolates that were recovered from patients in five hospitals in Armenia from January to August of 2019. A phylogenetic comparison revealed that our isolates were closely related to each other by their core and accessory genomes, despite having been isolated from different regions and hospitals in Armenia. We identified unique genes in our isolates and in a closely related isolate recovered in France. The unique genes (hemolysin E virulence gene, lactate utilization operon lutABC, and endonuclease restriction modification operon hsdMSR) were identified in three separate genomic regions that were adjacent to prophage genes, including one region containing the TonB-dependent iron siderophore receptor gene ireA, which was only found in 5 other ST127 isolates from the European Nucleotide Archive (ENA). We further identified that these isolates possessed unique virulence and metabolic genes and harbored antibiotic resistance genes, including the ESBL genes bla_CTX-M-3 (n = 3), bla_CTX-M-236 (n = 1), and bla_TEM-1 (n = 1), in addition to a quinolone resistance protein gene qnrD1 (n = 1), which was absent in the ST127 isolates obtained from the ENA. Moreover, a plasmid replicon gene IncI2 (n = 1) was unique to ARM88 of the Armenian isolates. Our findings demonstrate that at the time of this study, E. coli ST127 was a cause of urinary tract infections in patients in different regions of Armenia, with a possibility of cross-country transmission between Armenia and France. IMPORTANCE Whole-genome sequencing studies of pathogens causing infectious diseases are seriously lacking in Armenia, hampering global efforts to track, trace and contain infectious disease outbreaks. In this study, we report for the first-time the whole-genome sequencing and analysis of ESBL UPEC ST127 isolates recovered from hospitalized patients in Armenia and compare them with other E. coli ST127 retrieved from the ENA. We found close genetic similarities of the Armenian isolates, indicating that E. coli ST127 was potentially a dominant lineage causing urinary tract infections in Armenia. Furthermore, we identified unique genes that were horizontally acquired in the clusters of Armenian and French isolates that were absent in other ST127 isolates obtained from the ENA. Our findings highlight a possible cross-country transmission between Armenia and France and the idea that the implementation of WGS surveillance could contribute to global efforts in tackling antibiotic resistance, as bacteria carrying antimicrobial resistance (AMR) genes do not recognize borders.

The making of hypervirulent Klebsiella pneumoniae

Klebsiella pneumoniae is a notorious bacterium in clinical practice. Virulence, carbapenem-resistance and their convergence among K. pneumoniae are extensively discussed in this article. Hypervirulent K. pneumoniae (HvKP) has spread from the Asian Pacific Rim to the world, inducing various invasive infections, such as pyogenic liver abscess, endophthalmitis, and meningitis. Furthermore, HvKP has acquired more and more drug resistance. Among multidrug-resistant HvKP, hypervirulent carbapenem-resistant K. pneumoniae (Hv-CRKP), and carbapenem-resistant hypervirulent K. pneumoniae (CR-HvKP) are both devastating for their extreme drug resistance and virulence. The hypervirulence of HvKP is primarily attributed to hypercapsule, macromolecular exopolysaccharides, or excessive siderophores, although it has many other factors, for example, lipopolysaccharides, fimbriae, and porins. In contrast with classical determination of HvKP, that is, animal lethality test, molecular determination could be an optional and practical method after improvement. HvKP, including Hv-CRKP and CR-HvKP, has been progressing. R-M and CRISPR-Cas systems may play pivotal roles in such evolutions. Hv-CRKP and CR-HvKP, in particular the former, should be of severe concern due to their being more and more prevalent.

Molecular insights into DNA recognition and methylation by non-canonical type I restriction-modification systems

Type I restriction-modification systems help establish the prokaryotic DNA methylation landscape and provide protection against invasive DNA. In addition to classical m6A modifications, non-canonical type I enzymes catalyze both m6A and m4C using alternative DNA-modification subunits M1 and M2. Here, we report the crystal structures of the non-canonical PacII_M1M2S methyltransferase bound to target DNA and reaction product S-adenosylhomocysteine in a closed clamp-like conformation. Target DNA binds tightly within the central tunnel of the M1M2S complex and forms extensive contacts with all three protein subunits. Unexpectedly, while the target cytosine properly inserts into M2's pocket, the target adenine (either unmethylated or methylated) is anchored outside M1's pocket. A unique asymmetric catalysis is established where PacII_M1M2S has precisely coordinated the relative conformations of different subunits and evolved specific amino acids within M2/M1. This work provides insights into mechanisms of m6A/m4C catalysis and guidance for designing tools based on type I restriction-modification enzymes.

Phase-variable Type I methyltransferase M.NgoAV from Neisseria gonorrhoeae FA1090 regulates phasevarion expression and gonococcal phenotype

The restriction-modification (RM) systems are compared to a primitive, innate, prokaryotic immune system, controlling the invasion by foreign DNA, composed of methyltransferase (MTase) and restriction endonuclease. The biological significance of RM systems extends beyond their defensive function, but the data on the regulatory role of Type I MTases are limited. We have previously characterized molecularly a non-canonical Type I RM system, NgoAV, with phase-variable specificity, encoded by Neisseria gonorrhoeae FA1090. In the current work, we have investigated the impact of methyltransferase NgoAV (M.NgoAV) activity on gonococcal phenotype and on epigenetic control of gene expression. For this purpose, we have constructed and studied genetic variants (concerning activity and specificity) within M.NgoAV locus. Deletion of M.NgoAV or switch of its specificity had an impact on phenotype of N. gonorrhoeae. Biofilm formation and planktonic growth, the resistance to antibiotics, which target bacterial peptidoglycan or other antimicrobials, and invasion of human epithelial host cells were affected. The expression of genes was deregulated in gonococcal cells with knockout M.NgoAV gene and the variant with new specificity. For the first time, the existence of a phasevarion (phase-variable regulon), directed by phase-variable Type I MTase, is demonstrated.

Biological Features and In Planta Transcriptomic Analyses of a Microviridae Phage (CLasMV1) in "Candidatus Liberibacter asiaticus"

“Candidatus Liberibacter asiaticus” (CLas) is the causal agent of citrus Huanglongbing (HLB, also called citrus greening disease), a highly destructive disease threatening citrus production worldwide. A novel Microviridae phage (named CLasMV1) has been found to infect CLas, providing a potential therapeutic strategy for CLas/HLB control. However, little is known about the CLasMV1 biology. In this study, we analyzed the population dynamics of CLasMV1 between the insect vector of CLas, the Asian citrus psyllid (ACP, Diaphorina citri Kuwayama) and the holoparasitic dodder plant (Cuscuta campestris Yunck.); both acquired CLasMV1-infected CLas from an HLB citrus. All CLas-positive dodder samples were CLasMV1-positive, whereas only 32% of CLas-positive ACP samples were identified as CLasMV1-positive. Quantitative analyses showed a similar distribution pattern of CLasMV1 phage and CLas among eight citrus cultivars by presenting at highest abundance in the fruit pith and/or the center axis of the fruit. Transcriptome analyses revealed the possible lytic activity of CLasMV1 on CLas in fruit pith as evidenced by high-level expressions of CLasMV1 genes, and CLas genes related to cell wall biogenesis and remodeling to maintain the CLas cell envelope integrity. The up-regulation of CLas genes were involved in restriction–modification system that could involve possible phage resistance for CLas during CLasMV1 infection. In addition, the regulation of CLas genes involved in cell surface components and Sec pathway by CLasMV1 phage could be beneficial for phage infection. This study expanded our knowledge of CLasMV1 phage that will benefit further CLas phage research and HLB control.

Microbial Proteins in Stomach Biopsies Associated with Gastritis, Ulcer, and Gastric Cancer

(1) Background: Gastric cancer (GC) is the fourth leading cause of cancer-related deaths worldwide. Helicobacter pylori infection is a major risk factor, but other microbial species may also be involved. In the context of an earlier proteomics study of serum and biopsies of patients with gastroduodenal diseases, we explored here a simplified microbiome in these biopsies (H. pylori, Acinetobacter baumannii, Escherichia coli, Fusobacterium nucleatum, Bacteroides fragilis) on the protein level. (2) Methods: A cohort of 75 patients was divided into groups with respect to the findings of the normal gastric mucosa (NGM) and gastroduodenal disorders such as gastritis, ulcer, and gastric cancer (GC). The H. pylori infection status was determined. The protein expression analysis of the biopsy samples was carried out using high-definition mass spectrometry of the tryptic digest (label-free data-independent quantification and statistical analysis). (3) Results: The total of 304 bacterial protein matches were detected based on two or more peptide hits. Significantly regulated microbial proteins like virulence factor type IV secretion system protein CagE from H. pylori were found with more abundance in gastritis than in GC or NGM. This finding could reflect the increased microbial involvement in mucosa inflammation in line with current hypotheses. Abundant proteins across species were heat shock proteins and elongation factors. (4) Conclusions: Next to the bulk of human proteins, a number of species-specific bacterial proteins were detected in stomach biopsies of patients with gastroduodenal diseases, some of which, like those expressed by the cag pathogenicity island, may provide gateways to disease prevention without antibacterial intervention in order to reduce antibiotic resistance.

Metagenomic methylation patterns resolve bacterial genomes of unusual size and structural complexity

The plasticity of bacterial and archaeal genomes makes examining their ecological and evolutionary dynamics both exciting and challenging. The same mechanisms that enable rapid genomic change and adaptation confound current approaches for recovering complete genomes from metagenomes. Here, we use strain-specific patterns of DNA methylation to resolve complex bacterial genomes from long-read metagenomic data of a marine microbial consortium, the “pink berries” of the Sippewissett Marsh (USA). Unique combinations of restriction-modification (RM) systems encoded by the bacteria produced distinctive methylation profiles that were used to accurately bin and classify metagenomic sequences. Using this approach, we finished the largest and most complex circularized bacterial genome ever recovered from a metagenome (7.9 Mb with >600 transposons), the finished genome of Thiohalocapsa sp. PB-PSB1 the dominant bacteria in the consortia. From genomes binned by methylation patterns, we identified instances of horizontal gene transfer between sulfur-cycling symbionts (Thiohalocapsa sp. PB-PSB1 and Desulfofustis sp. PB-SRB1), phage infection, and strain-level structural variation. We also linked the methylation patterns of each metagenome-assembled genome with encoded DNA methyltransferases and discovered new RM defense systems, including novel associations of RM systems with RNase toxins.

Barriers to genetic manipulation of Enterococci: Current Approaches and Future Directions

Enterococcus faecalis and Enterococcus faecium are Gram-positive commensal gut bacteria that can also cause fatal infections. To study clinically relevant multi-drug resistant E. faecalis and E. faecium strains, methods are needed to overcome physical (thick cell wall) and enzymatic barriers that limit the transfer of foreign DNA and thus prevent facile genetic manipulation. Enzymatic barriers to DNA uptake identified in E. faecalis and E. faecium include type I, II and IV restriction modification systems and CRISPR-Cas. This review examines E. faecalis and E. faecium DNA defence systems and the methods with potential to overcome these barriers. DNA defence system bypass will allow the application of innovative genetic techniques to expedite molecular-level understanding of these important, but somewhat neglected, pathogens.

Morphological and Genetic Characterization of Eggerthella lenta Bacteriophage PMBT5

Eggerthella lenta is a common member of the human gut microbiome. We here describe the isolation and characterization of a putative virulent bacteriophage having E. lenta as host. The double-layer agar method for isolating phages was adapted to anaerobic conditions for isolating bacteriophage PMBT5 from sewage on a strictly anaerobic E. lenta strain of intestinal origin. For this, anaerobically grown E. lenta cells were concentrated by centrifugation and used for a 24 h phage enrichment step. Subsequently, this suspension was added to anaerobically prepared top (soft) agar in Hungate tubes and further used in the double-layer agar method. Based on morphological characteristics observed by transmission electron microscopy, phage PMBT5 could be assigned to the Siphoviridae phage family. It showed an isometric head with a flexible, noncontractile tail and a distinct single 45 nm tail fiber under the baseplate. Genome sequencing and assembly resulted in one contig of 30,930 bp and a mol% GC content of 51.3, consisting of 44 predicted protein-encoding genes. Phage-related proteins could be largely identified based on their amino acid sequence, and a comparison with metagenomes in the human virome database showed that the phage genome exhibits similarity to two distantly related phages.

Mining bacterial NGS data vastly expands the complete genomes of temperate phages

Temperate phages (active prophages induced from bacteria) help control pathogenicity, modulate community structure, and maintain gut homeostasis. Complete phage genome sequences are indispensable for understanding phage biology. Traditional plaque techniques are inapplicable to temperate phages due to their lysogenicity, curbing their identification and characterization. Existing bioinformatics tools for prophage prediction usually fail to detect accurate and complete temperate phage genomes. This study proposes a novel computational temperate phage detection method (TemPhD) mining both the integrated active prophages and their spontaneously induced forms (temperate phages) from next-generation sequencing raw data. Applying the method to the available dataset resulted in 192 326 complete temperate phage genomes with different host species, expanding the existing number of complete temperate phage genomes by more than 100-fold. The wet-lab experiments demonstrated that TemPhD can accurately determine the complete genome sequences of the temperate phages, with exact flanking sites, outperforming other state-of-the-art prophage prediction methods. Our analysis indicates that temperate phages are likely to function in the microbial evolution by (i) cross-infecting different bacterial host species; (ii) transferring antibiotic resistance and virulence genes and (iii) interacting with hosts through restriction-modification and CRISPR/anti-CRISPR systems. This work provides a comprehensively complete temperate phage genome database and relevant information, which can serve as a valuable resource for phage research.

Structural basis for broad anti-phage immunity by DISARM

In the evolutionary arms race against phage, bacteria have assembled a diverse arsenal of antiviral immune strategies. While the recently discovered DISARM (Defense Island System Associated with Restriction-Modification) systems can provide protection against a wide range of phage, the molecular mechanisms that underpin broad antiviral targeting but avoiding autoimmunity remain enigmatic. Here, we report cryo-EM structures of the core DISARM complex, DrmAB, both alone and in complex with an unmethylated phage DNA mimetic. These structures reveal that DrmAB core complex is autoinhibited by a trigger loop (TL) within DrmA and binding to DNA substrates containing a 5′ overhang dislodges the TL, initiating a long-range structural rearrangement for DrmAB activation. Together with structure-guided in vivo studies, our work provides insights into the mechanism of phage DNA recognition and specific activation of this widespread antiviral defense system.

DISARM (Defense Island System Associated with Restriction Modification) systems can provide bacteria with protection against a wide range of phage. Here, Bravo et al. determine cryo-EM structures of the core DISARM complex that shed light onto phage DNA recognition and activation of this widespread defense system.

Discovering Biological Conflict Systems Through Genome Analysis: Evolutionary Principles and Biochemical Novelty

Biological replicators, from genes within a genome to whole organisms, are locked in conflicts. Comparative genomics has revealed a staggering diversity of molecular armaments and mechanisms regulating their deployment, collectively termed biological conflict systems. These encompass toxins used in inter- and intraspecific interactions, self/nonself discrimination, antiviral immune mechanisms, and counter-host effectors deployed by viruses and intragenomic selfish elements. These systems possess shared syntactical features in their organizational logic and a set of effectors targeting genetic information flow through the Central Dogma, certain membranes, and key molecules like NAD⁺. These principles can be exploited to discover new conflict systems through sensitive computational analyses. This has led to significant advances in our understanding of the biology of these systems and furnished new biotechnological reagents for genome editing, sequencing, and beyond. We discuss these advances using specific examples of toxins, restriction-modification, apoptosis, CRISPR/second messenger-regulated systems, and other enigmatic nucleic acid-targeting systems. Expected final online publication date for the Annual Review of Biomedical Data Science, Volume 5 is August 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Expression and Purification of BsaXI Restriction Endonuclease and Engineering New Specificity From BsaXI Specificity Subunit

It is stated that BsaXI is a Type IIB restriction endonuclease (REase) that cleaves both sides of its recognition sequence 5′↓N9 AC N5 CTCC N10↓ 3′ (complement strand 5′ ↓N7 GGAG N5 GT N12↓ 3′), creating 3-base 3′ overhangs. Here we report the cloning and expression of bsaXIS and bsaXIRM genes in Escherichia coli. The BsaXI activity was successfully reconstituted by mixing the BsaXI RM fusion subunit with the BsaXI S subunit and the enzyme complex further purified by chromatography over 6 columns. As expected, the S subunit consisted of two subdomains encoding TRD1-CR1 [target recognition domain (TRD), conserved region (CR)] for 5′ AC 3′, and TRD2-CR2 presumably specifying 5′ CTCC 3′. TRD1-CR1 (TRD2-CR2 deletion) or duplication of TRD1 (TRD1-CR1-TRD1-CR2) both generated a new specificity 5′ AC N5 GT 3′ when the S variants were complexed with the RM subunits. The circular permutation of TRD1 and TRD2, i.e., the relocation of TRD2-CR2 to the N-terminus and TRD1-CR1 to the C-terminus generated the same specificity with the RM subunits, although some wobble cleavage was detected. The TRD2 domain in the BsaXI S subunit can be substituted by a close homolog (∼59% sequence identity) and generated the same specificity. However, TRD2-CR2 domain alone failed to express in E. coli, but CR1-TRD2-CR2 protein could be expressed and purified which showed partial nicking activity with the RM subunits. This work demonstrated that like Type I restriction systems, the S subunit of a Type IIB system could also be manipulated to create new specificities. The genome mining of BsaXI TRD2 homologs in GenBank found more than 36 orphan TRD2 homologs, implying that quite a few orphan TRD2s are present in microbial genomes that may be potentially paired with other TRDs to create new restriction specificities.

Impacts of Mycoplasma agalactiae restriction-modification systems on pan-epigenome dynamics and genome plasticity

DNA methylations play an important role in the biology of bacteria. Often associated with restriction modification (RM) systems, they are important drivers of bacterial evolution interfering in horizontal gene transfer events by providing a defence against foreign DNA invasion or by favouring genetic transfer through production of recombinogenic DNA ends. Little is known regarding the methylome of the Mycoplasma genus, which encompasses several pathogenic species with small genomes. Here, genome-wide detection of DNA methylations was conducted using single molecule real-time (SMRT) and bisulphite sequencing in several strains of Mycoplasma agalactiae, an important ruminant pathogen and a model organism. Combined with whole-genome analysis, this allowed the identification of 19 methylated motifs associated with three orphan methyltransferases (MTases) and eight RM systems. All systems had a homolog in at least one phylogenetically distinct Mycoplasma spp. Our study also revealed that several superimposed genetic events may participate in the M. agalactiae dynamic epigenomic landscape. These included (i) DNA shuffling and frameshift mutations that affect the MTase and restriction endonuclease content of a clonal population and (ii) gene duplication, erosion, and horizontal transfer that modulate MTase and RM repertoires of the species. Some of these systems were experimentally shown to play a major role in mycoplasma conjugative, horizontal DNA transfer. While the versatility of DNA methylation may contribute to regulating essential biological functions at cell and population levels, RM systems may be key in mycoplasma genome evolution and adaptation by controlling horizontal gene transfers.

Exploiting transposons in the study of Staphylococcus aureus pathogenesis and virulence

The opportunistic pathogen Staphylococcus aureus has an extremely complex relationship with humans. While the bacteria can exist as a commensal in many, it can cause a wide range of diseases and infections when turned pathogenic. Its presence is a determinant of chronicity and poor prognosis in numerous diseases, and its genomic plasticity causes S. aureus antimicrobial resistance to be one of the most dire contemporary medical problems to solve. Genetic manipulation of S. aureus has led to numerous findings that are vital in the fight against its pathogenesis. The utilisation of transposon mutant libraries for the systematic inspection of the S. aureus genome led to many landmark discoveries pertaining to the bacteria's pathogenicity, antimicrobial resistance acquisition, and virulence regulation. In this review, we describe mutant libraries, and their significant contributions, from various S. aureus strains created with commonly used transposons. The general workflow for the construction of libraries will be presented, along with a discussion of the challenges of undertaking the task of large-scale library construction. As the accessibility of transposon mutant library construction, screening, and analysis increases, this genetic tool could be further exploited in the study of the S. aureus genome.

Streptococcus mutans gene expression and functional profile in root caries: an RNA-seq study

The literature is still scarce on studies describing S. mutans global gene expression under clinical conditions such as those found on complex biofilms from sound root surfaces (SRS) and carious root surfaces (RC). This study aimed to investigate the S. mutans gene expression and functional profile within the metatranscriptome of biofilms from SRS and from RC in an attempt to identify enriched functional signatures potentially associated with healthy to disease transitioning process. Total RNA was extracted, and prepared libraries (SRS=10 and RC=9) were paired-end sequenced using the Illumina HiSeq2500. Read count assigned to each gene of the S. mutans UA159 strain were obtained. Differentially expressed genes (DEG) between SRS and RC were identified using the DESeq2 R package and weighted gene co-expression network analysis (WGCNA) was performed to explore and identify functional modules related to SRS and RC. We found seventeen DEG between SRS and RC samples, with three overexpressed in RC and related to membrane protein, alanyl-tRNA synthetase and GTP-binding protein with the remaining ones overexpressed in SRS samples and related to hypothetical protein, transposon integrase, histidine kinase, putative transporter, bacteriocin immunity protein, response regulator, 6-phospho-beta-galactosidase, purine metabolism and to transcriptional regulator. Key-functional modules were identified for SRS and RC conditions based on WCGNA, being 139 hub genes found on SRS key-module and 17 genes on RC key-module. Functional analysis of S. mutans within the metatranscriptome of biofilms from sound root and from carious root revealed a similar pattern of gene expression, and only a few genes have been differentially expressed between biofilms from sound root surfaces and from root carious lesions. However, S. mutans presented a greater functional abundance in the lesion samples. Some functional patterns related to sugar (starch, sucrose, fructose, mannose and lactose) and heterofermentative metabolisms, to cell-wall biosynthesis and to acid tolerance stress seem to be enriched on carious root surfaces conferring ecological advantages to S. mutans. Altogether, the present data suggest that a functional signature may be associated with carious root lesions.

Structural features of a minimal intact methyltransferase of a type I restriction-modification system

Type I restriction-modification enzymes are oligomeric proteins composed of methylation (M), DNA sequence-recognition (S), and restriction (R) subunits. The different bipartite DNA sequences of 2-4 consecutive bases are recognized by two discerned target recognition domains (TRDs) located at the two-helix bundle of the two conserved regions (CRs). Two M-subunits and a single S-subunit form an oligomeric protein that functions as a methyltransferase (M₂S₁ MTase). Here, we present the crystal structure of the intact MTase from Vibrio vulnificus YJ016 in complex with the DNA-mimicking Ocr protein and the S-adenosyl-L-homocysteine (SAH). This MTase includes the M-domain with a helix tail (M-tail helix) and the S_1/2-domain of a TRD and a CR α-helix. The Ocr binds to the cleft of the TRD surface and SAH is located in the pocket within the M-domain. The solution- and negative-staining electron microscopy-based reconstructed (M₁S_1/2)₂ structure reveals a symmetric (S_1/2)₂ assembly using two CR-helices and two M-tail helices as a pivot, which is plausible for recognizing two DNA regions of same sequence. The conformational flexibility of the minimal M₁S_1/2 MTase dimer indicates a particular state resembling the structure of M₂S₁ MTases.

Genomic Features and Pervasive Negative Selection in Rhodanobacter Strains Isolated from Nitrate and Heavy Metal Contaminated Aquifer

Rhodanobacter species dominate in the Oak Ridge Reservation (ORR) subsurface environments contaminated with acids, nitrate, metal radionuclides, and other heavy metals. To uncover the genomic features underlying adaptations to these mixed-waste environments and to guide genetic tool development, we sequenced the whole genomes of eight Rhodanobacter strains isolated from the ORR site. The genome sizes ranged from 3.9 to 4.2 Mb harboring 3,695 to 4,035 protein-coding genes and GC contents approximately 67%. Seven strains were classified as R. denitrificans and one strain, FW510-R12, as R. thiooxydans based on full length 16S rRNA sequences. According to gene annotation, the top two Cluster of Orthologous Groups (COGs) with high pan-genome expansion rates (Pan/Core gene ratio) were “replication, recombination and repair” and “defense mechanisms.” The denitrifying genes had high DNA homologies except the predicted protein structure variances in NosZ. In contrast, heavy metal resistance genes were diverse with between 7 to 34% of them were located in genomic islands, and these results suggested origins from horizontal gene transfer. Analysis of the methylation patterns in four strains revealed the unique 5mC methylation motifs. Most orthologs (78%) had ratios of nonsynonymous to synonymous substitutions (dN/dS) less than one when compared to the type strain 2APBS1, suggesting the prevalence of negative selection. Overall, the results provide evidence for the important roles of horizontal gene transfer and negative selection in genomic adaptation at the contaminated field site. The complex restriction-modification system genes and the unique methylation motifs in Rhodanobacter strains suggest the potential recalcitrance to genetic manipulation.

IMPORTANCE Despite the dominance of Rhodanobacter species in the subsurface of the contaminated Oak Ridge Reservation (ORR) site, very little is known about the mechanisms underlying their adaptions to the various stressors present at ORR. Recently, multiple Rhodanobacter strains have been isolated from the ORR groundwater samples from several wells with varying geochemical properties. Using Illumina, PacBio, and Oxford Nanopore sequencing platforms, we obtained the whole genome sequences of eight Rhodanobacter strains. Comparison of the whole genomes demonstrated the genetic diversity, and analysis of the long nanopore reads revealed the heterogeneity of methylation patterns in strains isolated from the same well. Although all strains contained a complete set of denitrifying genes, the predicted tertiary structures of NosZ differed. The sequence comparison results demonstrate the important roles of horizontal gene transfer and negative selection in adaptation. In addition, these strains may be recalcitrant to genetic manipulation due to the complex restriction-modification systems and methylations.

Transmission barrier of the blaKPC plasmid mediated by type I restriction-modification systems in Escherichia coli

BACKGROUND

Transportation of carbapenem-resistant plasmids contributes to carbapenem resistance in Gram-negative bacteria. KPC enzymes are the most clinically important enzymes among carbapenem-resistant Klebsiella pneumoniae, whereas the rate of blaKPC in Escherichia coli is low. The CRISPR-Cas system and restriction-modification system (R-M system) in bacteria defend against invading genomes. Currently, the role of the immune systems in the low rate of KPC-producing E. coli remains unclear.

OBJECTIVES

We investigated the relationship between immune systems and the low detection rate of blaKPC in E. coli.

METHODS

We searched for blaKPC among 1039 E. coli whole genomes available in GenBank using nucleotide BLAST. CRISPR-Cas systems and the R-M system were detected in all strains having the ST as blaKPC-positive strains. Nucleotide BLAST was used to search for protospacers on blaKPC plasmids. A conjugation assay was performed to determine whether the R-M system influences the acquisition of blaKPC plasmids by E. coli.

RESULTS

ST131 was the dominant ST of KPC-producing E. coli and IncN was the main plasmid type (12/32). CRISPR-Cas systems were frequently present in E. coli carrying blaKPC. Furthermore, CRISPR-Cas systems in E. coli didn't target plasmids with blaKPC. Type I R-M systems were rare in KPC-producing E. coli, but significantly over-represented in KPC-negative strains. E. coli DH5α with hsdR deletion accepted blaKPC-carrying plasmids, whereas those with hsdR complementation impeded blaKPC-carrying plasmid conjugation.

CONCLUSIONS

Horizontal transmission of blaKPC occurs among E. coli. The type I R-M system is associated with the defence against blaKPC plasmid transport into E. coli.

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

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

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

Recent Molecular Tools for the Genetic Manipulation of Highly Industrially Important Mucoromycota Fungi

Mucorales is the largest and most well-studied order of the phylum Mucormycota and is known for its rapid growth rate and various industrial applications. The Mucorales fungi are a fascinating group of filamentous organisms with many uses in research and the industrial and medical fields. They are widely used biotechnological producers of various secondary metabolites and other value-added products. Certain members of Mucorales are extensively used as model organisms for genetic and molecular investigation and have extended our understanding of the metabolisms of other members of this order as well. Compared with other fungal species, our understanding of Mucoralean fungi is still in its infancy, which could be linked to their lack of effective genetic tools. However, recent advancements in molecular tools and approaches, such as the construction of recyclable markers, silencing vectors, and the CRISPR-Cas9-based gene-editing system, have helped us to modify the genomes of these model organisms. Multiple genetic modifications have been shown to generate valuable products on a large scale and helped us to understand the morphogenesis, basic biology, pathogenesis, and host–pathogen interactions of Mucoralean fungi. In this review, we discuss various conventional and modern genetic tools and approaches used for efficient gene modification in industrially important members of Mucorales.

Advanced biotechnology using methyltransferase and its applications in bacteria: a mini review

Since prokaryotic restriction-modification (RM) systems protect the host by cleaving foreign DNA by restriction endonucleases, it is difficult to introduce engineered plasmid DNAs into newly isolated microorganisms whose RM system is not discovered. The prokaryotes also possess methyltransferases to protect their own DNA from the endonucleases. As those methyltransferases can be utilized to methylate engineered plasmid DNAs before transformation and to enhance the stability within the cells, the study on methyltransferases in newly isolated bacteria is essential for genetic engineering. Here, we introduce the mechanism of the RM system, specifically the methyltransferases and their biotechnological applications. These biotechnological strategies could facilitate plasmid DNA-based genetic engineering in bacteria strains that strongly defend against foreign DNA.

A Type I Restriction Modification System Influences Genomic Evolution Driven by Horizontal Gene Transfer in Paenibacillus polymyxa

Considered a “Generally Recognized As Safe” (GRAS) bacterium, the plant growth–promoting rhizobacterium Paenibacillus polymyxa has been widely applied in agriculture and animal husbandry. It also produces valuable compounds that are used in medicine and industry. Our previous work showed the presence of restriction modification (RM) system in P. polymyxa ATCC 842. Here, we further analyzed its genome and methylome by using SMRT sequencing, which revealed the presence of a larger number of genes, as well as a plasmid documented as a genomic region in a previous report. A number of mobile genetic elements (MGEs), including 78 insertion sequences, six genomic islands, and six prophages, were identified in the genome. A putative lysozyme-encoding gene from prophage P6 was shown to express lysin which caused cell lysis. Analysis of the methylome and genome uncovered a pair of reverse-complementary DNA methylation motifs which were widespread in the genome, as well as genes potentially encoding their cognate type I restriction-modification system PpoAI. Further genetic analysis confirmed the function of PpoAI as a RM system in modifying and restricting DNA. The average frequency of the DNA methylation motifs in MGEs was lower than that in the genome, implicating a role of PpoAI in restricting MGEs during genomic evolution of P. polymyxa. Finally, comparative analysis of R, M, and S subunits of PpoAI showed that homologs of the PpoAI system were widely distributed in species belonging to other classes of Firmicute, implicating a role of the ancestor of PpoAI in the genomic evolution of species beyond Paenibacillus.

DNA methylation by three Type I restriction modification systems of Escherichia coli does not influence gene regulation of the host bacterium

DNA methylation is a common epigenetic mark that influences transcriptional regulation, and therefore cellular phenotype, across all domains of life. In particular, both orphan methyltransferases and those from phasevariable restriction modification systems (RMSs) have been co-opted to regulate virulence epigenetically in many bacteria. We now show that three distinct non-phasevariable Type I RMSs in Escherichia coli have no measurable impact on gene expression, in vivo virulence, or any of 1190 in vitro growth phenotypes. We demonstrated this using both Type I RMS knockout mutants as well as heterologous installation of Type I RMSs into two E. coli strains. These data provide three clear and currently rare examples of restriction modification systems that have no impact on their host organism’s gene regulation. This leads to the possibility that other such nonregulatory methylation systems may exist, broadening our view of the potential role that RMSs may play in bacterial evolution.

SAMase of Bacteriophage T3 Inactivates Escherichia coli's Methionine S-Adenosyltransferase by Forming Heteropolymers

S-Adenosylmethionine lyase (SAMase) of bacteriophage T3 degrades the intracellular SAM pools of the host Escherichia coli cells, thereby inactivating a crucial metabolite involved in a plethora of cellular functions, including DNA methylation. SAMase is the first viral protein expressed upon infection, and its activity prevents methylation of the T3 genome. Maintenance of the phage genome in a fully unmethylated state has a profound effect on the infection strategy. It allows T3 to shift from a lytic infection under normal growth conditions to a transient lysogenic infection under glucose starvation. Using single-particle cryoelectron microscopy (cryo-EM) and biochemical assays, we demonstrate that SAMase performs its function by not only degrading SAM but also by interacting with and efficiently inhibiting the host's methionine S-adenosyltransferase (MAT), the enzyme that produces SAM. Specifically, SAMase triggers open-ended head-to-tail assembly of E. coli MAT into an unusual linear filamentous structure in which adjacent MAT tetramers are joined by two SAMase dimers. Molecular dynamics simulations together with normal mode analyses suggest that the entrapment of MAT tetramers within filaments leads to an allosteric inhibition of MAT activity due to a shift to low-frequency, high-amplitude active-site-deforming modes. The amplification of uncorrelated motions between active-site residues weakens MAT's substrate binding affinity, providing a possible explanation for the observed loss of function. We propose that the dual function of SAMase as an enzyme that degrades SAM and as an inhibitor of MAT activity has emerged to achieve an efficient depletion of the intracellular SAM pools. IMPORTANCE Self-assembly of enzymes into filamentous structures in response to specific metabolic cues has recently emerged as a widespread strategy of metabolic regulation. In many instances, filamentation of metabolic enzymes occurs in response to starvation and leads to functional inactivation. Here, we report that bacteriophage T3 modulates the metabolism of the host E. coli cells by recruiting a similar strategy: silencing a central metabolic enzyme by subjecting it to phage-mediated polymerization. This observation points to an intriguing possibility that virus-induced polymerization of the host metabolic enzymes is a common mechanism implemented by viruses to metabolically reprogram and subdue infected cells.

Comparative genomic analysis of obligately piezophilic Moritella yayanosii DB21MT-5 reveals bacterial adaptation to the Challenger Deep, Mariana Trench

Hadal trenches are the deepest but underexplored ecosystems on the Earth. Inhabiting the trench bottom is a group of micro-organisms termed obligate piezophiles that grow exclusively under high hydrostatic pressures (HHP). To reveal the genetic and physiological characteristics of their peculiar lifestyles and microbial adaptation to extreme high pressures, we sequenced the complete genome of the obligately piezophilic bacterium Moritella yayanosii DB21MT-5 isolated from the deepest oceanic sediment at the Challenger Deep, Mariana Trench. Through comparative analysis against pressure sensitive and deep-sea piezophilic Moritella strains, we identified over a hundred genes that present exclusively in hadal strain DB21MT-5. The hadal strain encodes fewer signal transduction proteins and secreted polysaccharases, but has more abundant metal ion transporters and the potential to utilize plant-derived saccharides. Instead of producing osmolyte betaine from choline as other Moritella strains, strain DB21MT-5 ferments on choline within a dedicated bacterial microcompartment organelle. Furthermore, the defence systems possessed by DB21MT-5 are distinct from other Moritella strains but resemble those in obligate piezophiles obtained from the same geographical setting. Collectively, the intensive comparative genomic analysis of an obligately piezophilic strain Moritella yayanosii DB21MT-5 demonstrates a depth-dependent distribution of energy metabolic pathways, compartmentalization of important metabolism and use of distinct defence systems, which likely contribute to microbial adaptation to the bottom of hadal trench.

Microbial Arsenal of Antiviral Defenses - Part I

Bacteriophages or phages are viruses that infect bacterial cells (for the scope of this review we will also consider viruses that infect Archaea). Constant threat of phage infection is a major force that shapes evolution of the microbial genomes. To withstand infection, bacteria had evolved numerous strategies to avoid recognition by phages or to directly interfere with phage propagation inside the cell. Classical molecular biology and genetic engineering have been deeply intertwined with the study of phages and host defenses. Nowadays, owing to the rise of phage therapy, broad application of CRISPR-Cas technologies, and development of bioinformatics approaches that facilitate discovery of new systems, phage biology experiences a revival. This review describes variety of strategies employed by microbes to counter phage infection, with a focus on novel systems discovered in recent years. First chapter covers defense associated with cell surface, role of small molecules, and innate immunity systems relying on DNA modification.

Improved transformation efficiency of group A Streptococcus by inactivation of a type I restriction modification system

Streptococcus pyogenes or group A Streptococcus (GAS) is a leading cause of bacterial pharyngitis, skin and soft tissue infections, life-threatening invasive infections, and the post-infectious autoimmune syndromes of acute rheumatic fever and post-streptococcal glomerulonephritis. Genetic manipulation of this important pathogen is complicated by resistance of the organism to genetic transformation. Very low transformation efficiency is attributed to recognition and degradation of introduced foreign DNA by a type I restriction-modification system encoded by the hsdRSM locus. DNA sequence analysis of this locus in ten GAS strains that had been previously transformed with an unrelated plasmid revealed that six of the ten harbored a spontaneous mutation in hsdR, S, or M. The mutations were all different, and at least five of the six were predicted to result in loss of function of the respective hsd gene product. The unexpected occurrence of such mutations in previously transformed isolates suggested that the process of transformation selects for spontaneous inactivating mutations in the Hsd system. We investigated the possibility of exploiting the increased transformability of hsd mutants by constructing a deletion mutation in hsdM in GAS strain 854, a clinical isolate representative of the globally dominant M1T1 clonal group. Mutant strain 854ΔhsdM exhibited a 5-fold increase in electrotransformation efficiency compared to the wild type parent strain and no obvious change in growth or off-target gene expression. We conclude that genetic transformation of GAS selects for spontaneous mutants in the hsdRSM restriction modification system. We propose that use of a defined hsdM mutant as a parent strain for genetic manipulation of GAS will enhance transformation efficiency and reduce the likelihood of selecting spontaneous hsd mutants with uncharacterized genotypes.

Genomes of Gut Bacteria from Nasonia Wasps Shed Light on Phylosymbiosis and Microbe-Assisted Hybrid Breakdown

Phylosymbiosis is a cross-system trend whereby microbial community relationships recapitulate the host phylogeny. In Nasonia parasitoid wasps, phylosymbiosis occurs throughout development, is distinguishable between sexes, and benefits host development and survival. Moreover, the microbiome shifts in hybrids as a rare Proteus bacterium in the microbiome becomes dominant. The larval hybrids then catastrophically succumb to bacterium-assisted lethality and reproductive isolation between the species. Two important questions for understanding phylosymbiosis and bacterium-assisted lethality in hybrids are (i) do the Nasonia bacterial genomes differ from other animal isolates and (ii) are the hybrid bacterial genomes the same as those in the parental species? Here, we report the cultivation, whole-genome sequencing, and comparative analyses of the most abundant gut bacteria in Nasonia larvae, Providencia rettgeri and Proteus mirabilis. Characterization of new isolates shows Proteus mirabilis forms a more robust biofilm than Providencia rettgeri and that, when grown in coculture, Proteus mirabilis significantly outcompetes Providencia rettgeri. Providencia rettgeri genomes from Nasonia are similar to each other and more divergent from pathogenic, human associates. Proteus mirabilis from Nasonia vitripennis, Nasonia giraulti, and their hybrid offspring are nearly identical and relatively distinct from human isolates. These results indicate that members of the larval gut microbiome within Nasonia are most similar to each other, and the strain of the dominant Proteus mirabilis in hybrids is resident in parental species. Holobiont interactions between shared, resident members of the wasp microbiome and the host underpin phylosymbiosis and hybrid breakdown.

IMPORTANCE Animal and plant hosts often establish intimate relationships with their microbiomes. In varied environments, closely related host species share more similar microbiomes, a pattern termed phylosymbiosis. When phylosymbiosis is functionally significant and beneficial, microbial transplants between host species and host hybridization can have detrimental consequences on host biology. In the Nasonia parasitoid wasp genus, which contains a phylosymbiotic gut community, both effects occur and provide evidence for selective pressures on the holobiont. Here, we show that bacterial genomes in Nasonia differ from other environments and harbor genes with unique functions that may regulate phylosymbiotic relationships. Furthermore, the bacteria in hybrids are identical to those in parental species, thus supporting a hologenomic tenet that the same members of the microbiome and the host genome impact phylosymbiosis, hybrid breakdown, and speciation.

Restriction Enzymes

Restriction enzymes provided the foundation on which molecular cloning was built, and they remain as essential tools in current recombinant DNA technology. The three classes of restriction enzymes and their features are introduced here.

Methylation analysis of Klebsiella pneumoniae from Portuguese hospitals

Klebsiella pneumoniae is an important nosocomial infectious agent with a high antimicrobial resistance (AMR) burden. The application of long read sequencing technologies is providing insights into bacterial chromosomal and putative extra-chromosomal genetic elements (PEGEs) associated with AMR, but also epigenetic DNA methylation, which is thought to play a role in cleavage of foreign DNA and expression regulation. Here, we apply the PacBio sequencing platform to eight Portuguese hospital isolates, including one carbapenemase producing isolate, to identify methylation motifs. The resulting assembled chromosomes were between 5.2 and 5.5Mbp in length, and twenty-six PEGEs were found. Four of our eight samples carry blaCTX-M-15, a dominant Extended Spectrum Beta Lactamase in Europe. We identified methylation motifs that control Restriction-Modification systems, including GATC of the DNA adenine methylase (Dam), which methylates N6-methyladenine (m6A) across all our K. pneumoniae assemblies. There was a consistent lack of methylation by Dam of the GATC motif downstream of two genes: fosA, a locus associated with low level fosfomycin resistance, and tnpB transposase on IncFIB(K) plasmids. Overall, we have constructed eight high quality reference genomes of K. pneumoniae, with insights into horizontal gene transfer and methylation m6A motifs.