A new type of DNA phosphorothioation-based antiviral system in archaea

Parabacteroides absconsus sp. nov., isolated from feces of a child with atopic dermatitis

An obligate anaerobic, Gram-negative rod-shaped bacterial strain designated AD58T was isolated from the feces of a 3-year-old boy with atopic dermatitis. The closest species is Parabacteroides fecalis with 96.98% 16S rRNA gene identity. The average nucleotide identity value between AD58T and P. fecalis AGMB00274T is 85.0-85.4%. The circular genome sequence is 3.77 Mbp long with 43.5 mol% G + C content. The strain AD58T grows at 32-42 °C, its pH range for growth is 6.0-7.5. No growth is observed in the presence of 1% or higher NaCl concentrations. The major fatty acids are anteiso-C15: 0, iso-C15: 0, and summed feature 11 (iso-C17: 0 3-OH and/or C18: 2 DMA), and the predominant respiratory quinone is MK-9. Conditioned media from AD58T increased expression of IL-8 but decreased expression of TLR-4 in HT29 cells. Based on the described properties, we propose AD58T as the type strain of Parabacteroides absconsus sp. nov. (= VKM B-3630T = JCM 35468T).

Molecular basis of the phosphorothioation-sensing antiphage defense system IscS-DndBCDE-DndI

Phosphorothioation serves as a DNA backbone modification mechanism, wherein a sulfur atom substitutes the nonbridging oxygen atom within the phosphodiester, facilitated by the gene products of dndABCDE or sspABCD. The combination of dndABCDE with dndFGH forms a bona fide defense system, where the DndFGH protein complex exhibits DNA nickase and DNA translocase activities to prevent phage invasion. In this study, we identified that dndI, co-transcribed with dndFGH, can independently couple with iscS-dndBCDE as an anti-phage defense system. Moreover, we resolved the crystal structure of DndI from Salmonella at a resolution of 3.10 Å. We discovered that its residue Y25, residing within a hydrophobic region of DndI, is involved in phosphorothioate (PT) sensing. Upon sensing PT modifications at 5'-GPSAAC-3'/5'-GPSTTC-3', the ATPase activity of DndI is stimulated, which subsequently triggers a conformational transition, facilitating the dissociation of DndI from self-DNA, thereby allowing DndI to avoid cleaving self-DNA while restricting PT-deficient phage DNA. This research broadens the knowledge of the mechanistic diversity underlying PT-based defense systems and highlights their complexity in the course of evolution.

A widespread phage-encoded kinase enables evasion of multiple host antiphage defence systems

DNA degradation (Dnd) is a widespread bacterial antiphage defence system that relies on DNA phosphorothioate (PT) modification for self/non-self discrimination and subsequent degradation of unmodified DNA. Phages employ counterstrategies to evade host immunity, but anti-Dnd immunity has not been characterized. Here we report an immune evasion protein encoded by the Salmonella phage JSS1 that contributes to subverting Dnd and other defence systems. Using quantitative proteomic and phosphoproteomic analyses, we show that the protein JSS1_004 employs N-terminal Ser/Thr/Tyr protein kinase activity to catalyse the multisite phosphorylation of host DndFGH. Notably, JSS1_004 also phosphorylates other bacterial immune systems to varying degrees, including CRISPR‒Cas, QatABCD, SIR2+HerA and DUF4297+HerA. Given that JSS1_004 and its homologues are widespread in phylogenetically diverse phages, we suggest that this strategy constitutes a family of immune evasion proteins that increases the chances of phage proliferation even when a host deploys multiple defence systems.

DNA Phosphorothioate Modification Systems and Associated Phage Defense Systems

In contrast to the well-known DNA methylation of nucleobases, DNA phosphorothioate (PT) modification occurs in the DNA sugar-phosphate backbone. The non-bridging oxygen is replaced by a sulfur atom, which increases the nuclease tolerance of the DNA. In recent years, we have witnessed advances in understanding of PT modification enzymes, the features of PT modification across prokaryotic genomes, and PT-related physiological functions. Although only a small fraction of modifiable recognition sites across bacterial genomes undergo PT modification, enzymes such as DndFGH and SspE can use this modification as a recognition marker to differentiate between self- and non-self-DNA, thus destroying PT-lacking invasive DNA and preventing autoimmunity. We highlight the molecular mechanisms of PT modification-associated defense systems. We also describe notable applications of PT systems in the engineering of phage-resistant bacterial strains, RNA editing, and nucleic acid detection.

Overcoming Bacteriophage Contamination in Bioprocessing: Strategies and Applications

Bacteriophage contamination has a devastating impact on the viability of bacterial hosts and can significantly reduce the productivity of bioprocesses in biotechnological industries. The consequences range from widespread fermentation failure to substantial economic losses, highlighting the urgent need for effective countermeasures. Conventional prevention methods, which focus primarily on the physical removal of bacteriophages from equipment, bioprocess units, and the environment, have proven ineffective in preventing phage entry and contamination. The coevolutionary dynamics between phages and their bacterial hosts have spurred the development of a diverse repertoire of antiviral defense mechanisms within microbial communities. These naturally occurring defense strategies can be harnessed through genetic engineering to convert phage-sensitive hosts into robust, phage-resistant cell factories, providing a strategic approach to mitigate the threats posed by bacteriophages to industrial bacterial processes. In this review, an overview of the various defense strategies and immune systems that curb the propagation of bacteriophages and highlight their applications in fermentation bioprocesses to combat phage contamination is provided. Additionally, the tactics employed by phages to circumvent these defense strategies are also discussed, as preventing the emergence of phage escape mutants is a key component of effective contamination management.

Catalytically active prokaryotic Argonautes employ phospholipase D family proteins to strengthen immunity against different genetic invaders

Prokaryotic Argonautes (pAgos) provide bacteria and archaea with immunity against plasmids and viruses. Catalytically active pAgos utilize short oligonucleotides as guides to directly cleave foreign nucleic acids, while inactive pAgos lacking catalytic residues employ auxiliary effectors, such as nonspecific nucleases, to trigger abortive infection upon detection of foreign nucleic acids. Here, we report a unique group of catalytically active pAgo proteins that frequently associate with a phospholipase D (PLD) family protein. We demonstrate that this particular system employs the catalytic center of the associated PLD protein rather than that of pAgo to restrict plasmid DNA, while interestingly, its immunity against a single-stranded DNA virus relies on the pAgo catalytic center and is enhanced by the PLD protein. We also find that this system selectively suppresses viral DNA propagation without inducing noticeable abortive infection outcomes. Moreover, the pAgo protein alone enhances gene editing, which is unexpectedly inhibited by the PLD protein. Our data highlight the ability of catalytically active pAgo proteins to employ auxiliary proteins to strengthen the targeted eradication of different genetic invaders and underline the trend of PLD nucleases to participate in host immunity.

Genomic analysis of Oceanotoga teriensis strain UFV_LIMV02, a multidrug-resistant thermophilic bacterium isolated from an offshore oil reservoir

Bacteria of the species Oceanotoga teriensis belong to the family Petrotogaceae, are Gram-negative bacilli, are moderately thermophilic and are included in the group of thiosulfate-reducing bacteria, being capable of significantly accelerating corrosion in metallic structures. However, no in-depth study on the genome, antibiotic resistance and mobile elements has been carried out so far. In this work, the isolation, phenotypic and genotypic characterization of the multi-resistant O. teriensis UFV_LIMV02 strain was carried out, from water samples from an offshore oil extraction platform in Rio de Janeiro (Brazil). We determined that the isolate has a genome of 2 812 778 bp in size, with 26 % GC content, organized into 34 contigs. Genomic annotation using Rapid Annotation using Subsystem Technology revealed the presence of genes related to resistance to antibiotics and heavy metals. By evaluating the antimicrobial resistance of the isolate using the disc diffusion technique, resistance was verified for the classes of antibiotics, beta-lactams, fluoroquinolones, aminoglycosides, sulfonamides, lincosamides and rifamycins, a total of 14 antibiotics. The search for genomic islands, prophages and defence systems against phage infection revealed the presence of five genomic islands in its genome, containing genes related to resistance to heavy metals and antibiotics, most of which are efflux pumps and several transposases. No prophage was found in its genome; however, nine different defence systems against phage infection were detected. When analysing the clustered regularly interspaced short palindromic repeat (CRISPR) systems, four CRISPR arrays, classified as types I-B and III-B, with 272 spacers, can provide the strain with immunity to different mobile genetic elements and bacteriophage infection. The results found in this study show that the isolate UFV_LIVM02 is an environmental bacterium, resistant to different classes of antibiotics, and that the proteins encoded by the predicted genomic islands may be associated with the development of greater resistance to antibiotics and heavy metals. They provide evidence that environmental bacteria found in offshore oil exploration residues may pose a risk for the spread of antibiotic resistance genes. More comprehensive studies on the microbial community present in oil waste are needed to assess the risks of horizontal gene transfer.

Phosphorothioate DNA modification by BREX Type 4 systems in the human gut microbiome

Among dozens of microbial DNA modifications regulating gene expression and host defense, phosphorothioation (PT) is the only known backbone modification, with sulfur inserted at a non-bridging oxygen by dnd and ssp gene families. Here we explored the distribution of PT genes in 13,663 human gut microbiome genomes, finding that 6.3% possessed dnd or ssp genes predominantly in Bacillota, Bacteroidota, and Pseudomonadota. This analysis uncovered several putative new PT synthesis systems, including Type 4 Bacteriophage Exclusion (BREX) brx genes, which were genetically validated in Bacteroides salyersiae. Mass spectrometric analysis of DNA from 226 gut microbiome isolates possessing dnd, ssp, and brx genes revealed 8 PT dinucleotide settings confirmed in 6 consensus sequences by PT-specific DNA sequencing. Genomic analysis showed PT enrichment in rRNA genes and depletion at gene boundaries. These results illustrate the power of the microbiome for discovering prokaryotic epigenetics and the widespread distribution of oxidation-sensitive PTs in gut microbes.

PT-seq: A method for metagenomic analysis of phosphorothioate epigenetics in complex microbial communities

Among dozens of known epigenetic marks, naturally occurring phosphorothioate (PT) DNA modifications are unique in replacing a non-bridging phosphate oxygen with redox-active sulfur and function in prokaryotic restriction-modification and transcriptional regulation. Interest in PTs has grown due to the widespread distribution of the dnd, ssp, and brx genes among bacteria and archaea, as well as the discovery of PTs in 5-10% of gut microbes. Efforts to map PTs in complex microbiomes using existing next-generation and direct sequencing technologies have failed due to poor sensitivity. Here we developed PT-seq as a high-sensitivity method to quantitatively map PTs across genomes and metagenomically identify PT-containing microbes in complex genomic mixtures. Like other methods for mapping PTs in individual genomes, PT-seq exploits targeted DNA strand cleavage at PTs by iodine, followed by sequencing library construction using ligation or template switching approaches. However, PT-specific sequencing reads are dramatically increased by adding steps to heat denature the DNA, block pre-existing 3'-ends, fragment DNA after T-tailing, and enrich iodine-induced breaks using biotin-labeling and streptavidin beads capture. Iterative optimization of the sensitivity and specificity of PT-seq is demonstrated with individual bacteria and human fecal DNA.

Temporal dynamics and metagenomics of phosphorothioate epigenomes in the human gut microbiome

Background

Epigenetic regulation of gene expression and host defense is well established in microbial communities, with dozens of DNA modifications comprising the epigenomes of prokaryotes and bacteriophage. Phosphorothioation (PT) of DNA, in which a chemically-reactive sulfur atom replaces a non-bridging oxygen in the sugar-phosphate backbone, is catalyzed by dnd and ssp gene families widespread in bacteria and archaea. However, little is known about the role of PTs or other microbial epigenetic modifications in the human microbiome. Here we optimized and applied fecal DNA extraction, mass spectrometric, and metagenomics technologies to characterize the landscape and temporal dynamics of gut microbes possessing PT modifications.

Results

Exploiting the nuclease-resistance of PTs, mass spectrometric analysis of limit digests of PT-containing DNA reveals PT dinucleotides as part of genomic consensus sequences, with 16 possible dinucleotide combinations. Analysis of mouse fecal DNA revealed a highly uniform spectrum of 11 PT dinucleotides in all littermates, with PTs estimated to occur in 5-10% of gut microbes. Though at similar levels, PT dinucleotides in fecal DNA from 11 healthy humans possessed signature combinations and levels of individual PTs. Comparison with a widely distributed microbial epigenetic mark, m6dA, suggested temporal dynamics consistent with expectations for gut microbial communities based on Taylor's Power Law. Application of PT-seq for site-specific metagenomic analysis of PT-containing bacteria in one fecal donor revealed the larger consensus sequences for the PT dinucleotides in Bacteroidota, Firmicutes, Actinobacteria, and Proteobacteria, which differed from unbiased metagenomics and suggested that the abundance of PT-containing bacteria did not simply mirror the spectrum of gut bacteria. PT-seq further revealed low abundance PT sites not detected as dinucleotides by mass spectrometry, attesting to the complementarity of the technologies.

Conclusions

The results of our studies provide a benchmark for understanding the behavior of an abundant and chemically-reactive epigenetic mark in the human gut microbiome, with implications for inflammatory conditions of the gut.

The binding affinity-dependent inhibition of cell growth and viability by DNA sulfur-binding domains

Increasing evidence suggests that DNA phosphorothioate (PT) modification serves several purposes in the bacterial host, and some restriction enzymes specifically target PT-DNA. PT-dependent restriction enzymes (PDREs) bind PT-DNA through their DNA sulfur binding domain (SBD) with dissociation constants (KD) of 5 nM~1 μM. Here, we report that SprMcrA, a PDRE, failed to dissociate from PT-DNA after cleavage due to high binding affinity, resulting in low DNA cleavage efficiency. Expression of SBDs in Escherichia coli cells with PT modification induced a drastic loss of cell viability at 25°C when both DNA strands of a PT site were bound, with one SBD on each DNA strand. However, at this temperature, SBD binding to only one PT DNA strand elicited a severe growth lag rather than lethality. This cell growth inhibition phenotype was alleviated by raising the growth temperature. An in vitro assay mimicking DNA replication and RNA transcription demonstrated that the bound SBD hindered the synthesis of new DNA and RNA when using PT-DNA as the template. Our findings suggest that DNA modification-targeting proteins might regulate cellular processes involved in DNA metabolism in addition to being components of restriction-modification systems and epigenetic readers.

Phage treatment of Pseudomonas aeruginosa yields a phage-resistant population with different susceptibility to innate immune responses and mild effects on metabolic profiles

In this study, we have investigated innate immune activation capacity and metabolic features of a population of P. aeruginosa PAO1 phage-resistant mutants with diverse genetic modification (large genomic deletions and point mutations) arising after exposure to phages targetting lipopolysaccharide (LPS) or Type-4 pili (T4P). Deletions led to the loss of genes involved in LPS synthesis, cell envelope permeability, efflux systems, biofilm production, oxidative stress tolerance, and DNA repair. Loss of LPS O antigen resulted in bacterial sensitivity to serum complement and stimulation of inflammatory cascades but did not cause increased phagocytosis, while T4P phage-resistant mutants were more effectively phagocytized than LPS-defective mutants. Changes in the utilization of different carbon, nitrogen, sulphur, and phosphorus sources were identified, especially in mutants where the two phage DNA persisted in the bacterial population (pseudolysogeny). However, the metabolic changes did not directly correlate with single-gene mutations or the large gene deletions, suggesting they reflect adaptive changes to the gene modifications that arise during the selection of resistant mutants. In contrast, phage-resistant mutants were susceptible to humoral innate immune responses, suggesting that phage resistance may be a beneficial outcome of phage therapy.

Ru-Catalyzed and Selenium-Directed Selective Formation of ortho- and Dialkenylated Selanes, Mixed Organoselenoethers, and Isoselenochromenes

Herein, the Ru-catalyzed chemo- and regioselective formation of four novel organoselenium compounds is established. Mono- and dialkenylated selanes were formed by the Se-directed o-C-H activation of benzyl(phenyl)selanes with alkynes. Unprecedented debenzylative/dearylative hydroselenations of alkynes were performed by slightly varying the amount of catalyst and temperature. Catalyst-driven direct formation of novel isoselenochromenes is also recorded. Altogether, 45 new organoseleno compounds were synthesized in good amounts with varieties of alkynes and selanes. This is the first report of its kind to deal with the synthesis of novel, challenging, and unusual organoseleno compounds in one reaction.

Reappraisal of the DNA phosphorothioate modification machinery: uncovering neglected functional modalities and identification of new counter-invader defense systems

The DndABCDE systems catalysing the unusual phosphorothioate (PT) DNA backbone modification, and the DndFGH systems, which restrict invasive DNA, have enigmatic and paradoxical features. Using comparative genomics and sequence-structure analyses, we show that the DndABCDE module is commonly functionally decoupled from the DndFGH module. However, the modification gene-neighborhoods encode other nucleases, potentially acting as the actual restriction components or suicide effectors limiting propagation of the selfish elements. The modification module's core consists of a coevolving gene-pair encoding the DNA-scanning apparatus - a DndD/CxC-clade ABC ATPase and DndE with two ribbon-helix-helix (MetJ/Arc) DNA-binding domains. Diversification of DndE's DNA-binding interface suggests a multiplicity of target specificities. Additionally, many systems feature DNA cytosine methylase genes instead of PT modification, indicating the DndDE core can recruit other nucleobase modifications. We show that DndFGH is a distinct counter-invader system with several previously uncharacterized domains, including a nucleotide kinase. These likely trigger its restriction endonuclease domain in response to multiple stimuli, like nucleotides, while blocking protective modifications by invader methylases. Remarkably, different DndH variants contain a HerA/FtsK ATPase domain acquired from multiple sources, including cellular genome-segregation systems and mobile elements. Thus, we uncovered novel HerA/FtsK-dependent defense systems that might intercept invasive DNA during replication, conjugation, or packaging.

The highly diverse antiphage defence systems of bacteria

Bacteria and their viruses have coevolved for billions of years. This ancient and still ongoing arms race has led bacteria to develop a vast antiphage arsenal. The development of high-throughput screening methods expanded our knowledge of defence systems from a handful to more than a hundred systems, unveiling many different molecular mechanisms. These findings reveal that bacterial immunity is much more complex than previously thought. In this Review, we explore recently discovered bacterial antiphage defence systems, with a particular focus on their molecular diversity, and discuss the ecological and evolutionary drivers and implications of the existing diversity of antiphage defence mechanisms.

Characterization of an archaeal virus-host system reveals massive genomic rearrangements in a laboratory strain

Halophilic archaea (haloarchaea) are known to exhibit multiple chromosomes, with one main chromosome and one or several smaller secondary chromosomes or megaplasmids. Halorubrum lacusprofundi, a model organism for studying cold adaptation, exhibits one secondary chromosome and one megaplasmid that include a large arsenal of virus defense mechanisms. We isolated a virus (Halorubrum tailed virus DL1, HRTV-DL1) infecting Hrr. lacusprofundi, and present an in-depth characterization of the virus and its interactions with Hrr. lacusprofundi. While studying virus-host interactions between Hrr. lacusprofundi and HRTV-DL1, we uncover that the strain in use (ACAM34_UNSW) lost the entire megaplasmid and about 38% of the secondary chromosome. The loss included the majority of virus defense mechanisms, making the strain sensitive to HRTV-DL1 infection, while the type strain (ACAM34_DSMZ) appears to prevent virus replication. Comparing infection of the type strain ACAM34_DSMZ with infection of the laboratory derived strain ACAM34_UNSW allowed us to identify host responses to virus infection that were only activated in ACAM34_UNSW upon the loss of virus defense mechanisms. We identify one of two S-layer proteins as primary receptor for HRTV-DL1 and conclude that the presence of two different S-layer proteins in one strain provides a strong advantage in the arms race with viruses. Additionally, we identify archaeal homologs to eukaryotic proteins potentially being involved in the defense against virus infection.

A DNA phosphorothioation-based Dnd defense system provides resistance against various phages and is compatible with the Ssp defense system

DndABCDE-catalyzed DNA phosphorothioation (PT), in which the nonbridging oxygen is swapped with a sulfur atom, was first identified in the bacterial genome. Usually, this modification gene cluster is paired with a restriction module consisting of DndF, DndG, and DndH. Although the mechanisms for the antiphage activity conferred by this Dnd-related restriction and modification (R-M) system have been well characterized, several features remain unclear, including the antiphage spectrum and potential interference with DNA methylation. Recently, a novel PT-related R-M system, composed of the modification module SspABCD paired with a single restriction enzyme, SspE, was revealed to be widespread in the bacterial kingdom, which aroused our interest in the interaction between Dnd- and Ssp-based R-M systems. In this study, we discussed the action of Dnd-related R-M systems against phages and demonstrated that the host could benefit from the protection provided by Dnd-related R-M systems against infection by various lytic phages as well as temperate phages. However, this defense barrier would fail against lysogenic phages. Interestingly, DNA methylation, even in the consensus sequence recognized by the Dnd system, could not weaken the restriction efficiency. Finally, we explored the interaction between Dnd- and Ssp-based R-M systems and found that these two systems were compatible. This study not only expands our knowledge of Dnd-associated R-M systems but also reveals a complex interaction between different defense barriers that coexist in the cell. IMPORTANCE Recently, we decoded the mechanism of Dnd-related R-M systems against genetic parasites. In the presence of exogenous DNA that lacks PT, the macromolecular machine consisting of DndF, DndG, and DndH undergoes conformational changes to perform DNA binding, translocation, and DNA nicking activities and scavenge the foreign DNA. However, several questions remain unanswered, including questions regarding the antiphage spectrum, potential interference by DNA methylation, and interplay with other PT-dependent R-M systems. Here, we revealed that the host could benefit from Dnd-related R-M systems for a broad range of antiphage activities, regardless of the presence of DNA methylation. Furthermore, we demonstrated that the convergence of Dnd- and Ssp-related R-M systems could confer to the host a stronger antiphage ability through the additive suppression of phage replication. This study not only deepens our understanding of PT-related defense barriers but also expands our knowledge of the arms race between bacteria and their predators.

Recurring and emerging themes in prokaryotic innate immunity

A resurgence of interest in the pathways that bacteria use to protect against their viruses (i.e. phages) has led to the discovery of dozens of new antiphage defenses. Given the sheer abundance and diversity of phages - the ever-evolving targets of immunity - it is not surprising that these newly described defenses are also remarkably diverse. However, as their mechanisms slowly come into focus, some common strategies and themes are also beginning to emerge. This review highlights recurring and emerging themes in the mechanisms of innate immunity in bacteria and archaea, with an emphasis on recently described systems that have undergone more thorough mechanistic characterization.

The DNA Phosphorothioation Restriction-Modification System Influences the Antimicrobial Resistance of Pathogenic Bacteria

Bacterial defense barriers, such as DNA methylation-associated restriction-modification (R-M) and the CRISPR-Cas system, play an important role in bacterial antimicrobial resistance (AMR). Recently, a novel R-M system based on DNA phosphorothioate (PT) modification has been shown to be widespread in the kingdom of Bacteria as well as Archaea. However, the potential role of the PT R-M system in bacterial AMR remains unclear. In this study, we explored the role of PT R-Ms in AMR with a series of common clinical pathogenic bacteria. By analyzing the distribution of AMR genes related to mobile genetic elements (MGEs), it was shown that the presence of PT R-M effectively reduced the distribution of horizontal gene transfer (HGT)-derived AMR genes in the genome, even in the bacteria that did not tend to acquire AMR genes by HGT. In addition, unique gene variation analysis based on pangenome analysis and MGE prediction revealed that the presence of PT R-M could suppress HGT frequency. Thus, this is the first report showing that the PT R-M system has the potential to repress HGT-derived AMR gene acquisition by reducing the HGT frequency. IMPORTANCE In this study, we demonstrated the effect of DNA PT modification-based R-M systems on horizontal gene transfer of AMR genes in pathogenic bacteria. We show that there is no apparent association between the genetic background of the strains harboring PT R-Ms and the number of AMR genes or the kinds of gene families. The strains equipped with PT R-M harbor fewer plasmid-derived, prophage-derived, or integrating mobile genetic element (iMGE)-related AMR genes and have a lower HGT frequency, but the degree of inhibition varies among different bacteria. In addition, compared with Salmonella enterica and Escherichia coli, Klebsiella pneumoniae prefers to acquire MGE-derived AMR genes, and there is no coevolution between PT R-M clusters and bacterial core genes.

Nicking mechanism underlying the DNA phosphorothioate-sensing antiphage defense by SspE

DNA phosphorothioate (PT) modification, with a nonbridging phosphate oxygen substituted by sulfur, represents a widespread epigenetic marker in prokaryotes and provides protection against genetic parasites. In the PT-based defense system Ssp, SspABCD confers a single-stranded PT modification of host DNA in the 5'-C_PSCA-3' motif and SspE impedes phage propagation. SspE relies on PT modification in host DNA to exert antiphage activity. Here, structural and biochemical analyses reveal that SspE is preferentially recruited to PT sites mediated by the joint action of its N-terminal domain (NTD) hydrophobic cavity and C-terminal domain (CTD) DNA binding region. PT recognition enlarges the GTP-binding pocket, thereby increasing GTP hydrolysis activity, which subsequently triggers a conformational switch of SspE from a closed to an open state. The closed-to-open transition promotes the dissociation of SspE from self PT-DNA and turns on the DNA nicking nuclease activity of CTD, enabling SspE to accomplish self-nonself discrimination and limit phage predation, even when only a small fraction of modifiable consensus sequences is PT-protected in a bacterial genome.

Systematic strategies for developing phage resistant Escherichia coli strains

Phages are regarded as powerful antagonists of bacteria, especially in industrial fermentation processes involving bacteria. While bacteria have developed various defense mechanisms, most of which are effective against a narrow range of phages and consequently exert limited protection from phage infection. Here, we report a strategy for developing phage-resistant Escherichia coli strains through the simultaneous genomic integration of a DNA phosphorothioation-based Ssp defense module and mutations of components essential for the phage life cycle. The engineered E. coli strains show strong resistance against diverse phages tested without affecting cell growth. Additionally, the resultant engineered phage-resistant strains maintain the capabilities of producing example recombinant proteins, D-amino acid oxidase and coronavirus-encoded nonstructural protein nsp8, even under high levels of phage cocktail challenge. The strategy reported here will be useful for developing engineered E. coli strains with improved phage resistance for various industrial fermentation processes for producing recombinant proteins and chemicals of interest.

Phage contamination is a persistent problem in industrial biotechnology processes employing bacterial strains. Here, the authors report the construction of E. coli host strains with broad antiphase activities via the genomic integration of the Ssp defense system and mutations of components essential for phage infection cycles.

Microfluidics-Assisted Fluorescence Mapping of DNA Phosphorothioation

As the key player of a new restriction modification system, DNA phosphorothioate (PT) modification, which swaps oxygen for sulfur on the DNA backbone, protects the bacterial host from foreign DNA invasion. The identification of PT sites helps us understand its physiological defense mechanisms, but accurately quantifying this dynamic modification remains a challenge. Herein, we report a simple quantitative analysis method for optical mapping of PT sites in the single bacterial genome. DNA molecules are fully stretched and immobilized in a microfluidic chip by capillary flow and electrostatic interactions, improving the labeling efficiency by maximizing exposure of PT sites on DNA while avoiding DNA loss and damage. After screening 116 candidates, we identified a bifunctional chemical compound, iodoacetyl-polyethylene glycol-biotin, that can noninvasively and selectively biotinylate PT sites, enabling further labeling with streptavidin fluorescent nanoprobes. With this method, PT sites in PT⁺ DNA can be easily detected by fluorescence, while almost no detectable ones were found in PT^- DNA, achieving real-time visualization of PT sites on a single DNA molecule. Collectively, this facile genome-wide PT site detection method directly characterizes the distribution and frequency of DNA modification, facilitating a better understanding of its modification mechanism that can be potentially extended to label DNAs in different species.

PADLOC: a web server for the identification of antiviral defence systems in microbial genomes

Most bacteria and archaea possess multiple antiviral defence systems that protect against infection by phages, archaeal viruses and mobile genetic elements. Our understanding of the diversity of defence systems has increased greatly in the last few years, and many more systems likely await discovery. To identify defence-related genes, we recently developed the Prokaryotic Antiviral Defence LOCator (PADLOC) bioinformatics tool. To increase the accessibility of PADLOC, we describe here the PADLOC web server (freely available at https://padloc.otago.ac.nz), allowing users to analyse whole genomes, metagenomic contigs, plasmids, phages and archaeal viruses. The web server includes a more than 5-fold increase in defence system types detected (since the first release) and expanded functionality enabling detection of CRISPR arrays and retron ncRNAs. Here, we provide user information such as input options, description of the multiple outputs, limitations and considerations for interpretation of the results, and guidance for subsequent analyses. The PADLOC web server also houses a precomputed database of the defence systems in > 230,000 RefSeq genomes. These data reveal two taxa, Campylobacterota and Spriochaetota, with unusual defence system diversity and abundance. Overall, the PADLOC web server provides a convenient and accessible resource for the detection of antiviral defence systems.

Horizontal Gene Transfer in Archaea-From Mechanisms to Genome Evolution

Archaea remains the least-studied and least-characterized domain of life despite its significance not just to the ecology of our planet but also to the evolution of eukaryotes. It is therefore unsurprising that research into horizontal gene transfer (HGT) in archaea has lagged behind that of bacteria. Indeed, several archaeal lineages may owe their very existence to large-scale HGT events, and thus understanding both the molecular mechanisms and the evolutionary impact of HGT in archaea is highly important. Furthermore, some mechanisms of gene exchange, such as plasmids that transmit themselves via membrane vesicles and the formation of cytoplasmic bridges that allows transfer of both chromosomal and plasmid DNA, may be archaea specific. This review summarizes what we know about HGT in archaea, and the barriers that restrict it, highlighting exciting recent discoveries and pointing out opportunities for future research. Expected final online publication date for the Annual Review of Microbiology, Volume 76 is September 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Molecular recognition between bacterial phosphorothioate DNA and sulfur-binding domain (SBD): competition between the water cage and chalcogen-hydrophobic packet

Bacterial DNA phosphorothioation (PT) physiologically and stereo-specifically replaces a non-bridging oxygen in a phosphate link with a sulfur atom, which can be recognized by a highly conserved sulfur-binding domain (SBD). Here we conducted thermodynamic integration (TI), molecular dynamics simulation, and quantum chemical calculations to decipher the specific molecular interactions between PT-DNA and SBD in Streptomyces coelicolor type IV restriction enzyme ScoMcrA. The TI-calculated binding affinity of (5'-CCG_Rp-PSGCCGG-3')₂ is larger than that of (5'-CCGGCCGG-3')₂ by about 7.4-7.7 kcal mol^-1. The binding difference dominantly stems from hydration energy of non-phosphorothioate DNA (9.8-10.6 kcal mol^-1) in aqueous solution, despite the persistent preference of 2.6-3.2 kcal mol^-1 in the DNA-SBD MD simulations. Furthermore, the quantum chemical calculations reveal an unusual non-covalent interaction in the phosphorothioate-binding scenario, where the ^PS⋯N^P165 chalcogen bond prevails the ^PS⋯HC_β vdW interactions from the adjacent residues H116-R117-Y164-P165-A168. Thus, the chalcogen-hydrophobic interaction pulls PT-DNA into the SBD binding pocket while the water cage pulls a normal DNA molecule out. The synergetic mechanism suggests the special roles of the proline pyrrolidine group in the SBD proteins, consistent with the experimental observations in the X-ray crystallography and structural bioinformatics analysis.

Involvement of the DNA Phosphorothioation System in TorR Binding and Anaerobic TMAO Respiration in Salmonella enterica

Although the phosphorothioate (PT) modification, in which the nonbridging oxygen in the DNA sugar-phosphate backbone is replaced by sulfur, has been reported to play versatile roles in multiple cellular processes, very little data have been obtained to define the role of PT in epigenetic regulation. In this study, we report that the PT system in Salmonella enterica serovar Cerro 87 is involved in the transcriptional regulation of the torCAD operon encoding the trimethylamine N-oxide (TMAO) respiration machinery that enables the use of TMAO as a terminal electron acceptor for respiration when oxygen is not available. In vitro, PT enhanced the binding of the transcriptional activator of the torCAD operon, namely, TorR, to its DNA substrate (tor boxes). However, in vivo, the PT modification protein complex DndCDE downregulated torCAD transcription through competing with the binding of TorR to the tor boxes. The altered expression of torCAD caused by PT modification proteins affected cell growth that relied on TMAO respiration. To our knowledge, this is the first report supporting that PT proteins participate in transcriptional regulation, showing a new function of PT systems.

Structural and Functional Analysis of DndE Involved in DNA Phosphorothioation in the Haloalkaliphilic Archaea Natronorubrum bangense JCM10635

Phosphorothioate (PT) modification, a sequence-specific modification that replaces the nonbridging oxygen atom with sulfur in a DNA phosphodiester through the gene products of dndABCDE or sspABCD, is widely distributed in prokaryotes. DNA PT modification functions together with gene products encoded by dndFGH, pbeABCD, or sspE to form defense systems that can protect against invasion by exogenous DNA particles. While the functions of the multiple enzymes in the PT system have been elucidated, the exact role of DndE in the PT process is still obscure. Here, we solved the crystal structure of DndE from the haloalkaliphilic archaeal strain Natronorubrum bangense JCM10635 at a resolution of 2.31 Å. Unlike the tetrameric conformation of DndE in Escherichia coli B7A, DndE from N. bangense JCM10635 exists in a monomeric conformation and can catalyze the conversion of supercoiled DNA to nicked or linearized products. Moreover, DndE exhibits preferential binding affinity to nicked DNA by virtue of the R19- and K23-containing positively charged surface. This work provides insight into how DndE functions in PT modification and the potential sulfur incorporation mechanism of DNA PT modification.

Diverse DNA modification in marine prokaryotic and viral communities

DNA chemical modifications, including methylation, are widespread and play important roles in prokaryotes and viruses. However, current knowledge of these modification systems is severely biased towards a limited number of culturable prokaryotes, despite the fact that a vast majority of microorganisms have not yet been cultured. Here, using single-molecule real-time sequencing, we conducted culture-independent 'metaepigenomic' analyses (an integrated analysis of metagenomics and epigenomics) of marine microbial communities. A total of 233 and 163 metagenomic-assembled genomes (MAGs) were constructed from diverse prokaryotes and viruses, respectively, and 220 modified motifs and 276 DNA methyltransferases (MTases) were identified. Most of the MTase genes were not genetically linked with the endonuclease genes predicted to be involved in defense mechanisms against extracellular DNA. The MTase-motif correspondence found in the MAGs revealed 10 novel pairs, 5 of which showed novel specificities and experimentally confirmed the catalytic specificities of the MTases. We revealed novel alternative specificities in MTases that are highly conserved in Alphaproteobacteria, which may enhance our understanding of the co-evolutionary history of the methylation systems and the genomes. Our findings highlight diverse unexplored DNA modifications that potentially affect the ecology and evolution of prokaryotes and viruses in nature.

Phosphorothioate-DNA bacterial diet reduces the ROS levels in C. elegans while improving locomotion and longevity

DNA phosphorothioation (PT) is widely distributed in the human gut microbiome. In this work, PT-diet effect on nematodes was studied with PT-bioengineering bacteria. We found that the ROS level decreased by about 20–50% and the age-related lipofuscin accumulation was reduced by 15–25%. Moreover, the PT-feeding worms were more active at all life periods, and more resistant to acute stressors. Intriguingly, their lifespans were prolonged by ~21.7%. Comparative RNA-seq analysis indicated that many gene expressions were dramatically regulated by PT-diet, such as cysteine-rich protein (scl-11/12/13), sulfur-related enzyme (cpr-2), longevity gene (jnk-1) and stress response (sod-3/5, gps-5/6, gst-18/20, hsp-12.6). Both the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis suggested that neuroactivity pathways were upregulated, while phosphoryl transfer and DNA-repair pathways were down-regulated in good-appetite young worms. The findings pave the way for pro-longevity of multicellular organisms by PT-bacterial interference.

Qiang Huang et al. fed C. elegans with E. coli containing phosphorothioate (PT) DNA or a control strain and evaluated the impact on animal physiology. They observed that worms fed PT( + ) diets exhibited low reactive oxygen species, more active movement, and a longer lifespan compared to controls, suggesting that PT-DNA may have a positive effect on animal health.

The origin and impeded dissemination of the DNA phosphorothioation system in prokaryotes

Phosphorothioate (PT) modification by the dnd gene cluster is the first identified DNA backbone modification and constitute an epigenetic system with multiple functions, including antioxidant ability, restriction modification, and virus resistance. Despite these advantages for hosting dnd systems, they are surprisingly distributed sporadically among contemporary prokaryotic genomes. To address this ecological paradox, we systematically investigate the occurrence and phylogeny of dnd systems, and they are suggested to have originated in ancient Cyanobacteria after the Great Oxygenation Event. Interestingly, the occurrence of dnd systems and prophages is significantly negatively correlated. Further, we experimentally confirm that PT modification activates the filamentous phage SW1 by altering the binding affinity of repressor and the transcription level of its encoding gene. Competition assays, concurrent epigenomic and transcriptomic sequencing subsequently show that PT modification affects the expression of a variety of metabolic genes, which reduces the competitive fitness of the marine bacterium Shewanella piezotolerans WP3. Our findings strongly suggest that a series of negative effects on microorganisms caused by dnd systems limit horizontal gene transfer, thus leading to their sporadic distribution. Overall, our study reveals putative evolutionary scenario of the dnd system and provides novel insights into the physiological and ecological influences of PT modification.

Evolution of Type IV CRISPR-Cas Systems: Insights from CRISPR Loci in Integrative Conjugative Elements of Acidithiobacillia

Type IV CRISPR-Cas are a distinct variety of highly derived CRISPR-Cas systems that appear to have evolved from type III systems through the loss of the target-cleaving nuclease and partial deterioration of the large subunit of the effector complex. All known type IV CRISPR-Cas systems are encoded on plasmids, integrative and conjugative elements (ICEs), or prophages, and are thought to contribute to competition between these elements, although the mechanistic details of their function remain unknown. There is a clear parallel between the compositions and likely origin of type IV and type I systems recruited by Tn7-like transposons and mediating RNA-guided transposition. We investigated the diversity and evolutionary relationships of type IV systems, with a focus on those in Acidithiobacillia, where this variety of CRISPR is particularly abundant and always found on ICEs. Our analysis revealed remarkable evolutionary plasticity of type IV CRISPR-Cas systems, with adaptation and ancillary genes originating from different ancestral CRISPR-Cas varieties, and extensive gene shuffling within the type IV loci. The adaptation module and the CRISPR array apparently were lost in the type IV ancestor but were subsequently recaptured by type IV systems on several independent occasions. We demonstrate a high level of heterogeneity among the repeats with type IV CRISPR arrays, which far exceed the heterogeneity of any other known CRISPR repeats and suggest a unique adaptation mechanism. The spacers in the type IV arrays, for which protospacers could be identified, match plasmid genes, in particular those encoding the conjugation apparatus components. Both the biochemical mechanism of type IV CRISPR-Cas function and their role in the competition among mobile genetic elements remain to be investigated.

Sulfur modification in natural RNA and therapeutic oligonucleotides

Sulfur modifications have been discovered on both DNA and RNA. Sulfur substitution of oxygen atoms at nucleobase or backbone locations in the nucleic acid framework led to a wide variety of sulfur-modified nucleosides and nucleotides. While the discovery, regulation and functions of DNA phosphorothioate (PS) modification, where one of the non-bridging oxygen atoms is replaced by sulfur on the DNA backbone, are important topics, this review focuses on the sulfur modification in natural cellular RNAs and therapeutic nucleic acids. The sulfur modifications on RNAs exhibit diversity in terms of modification location and cellular function, but the various sulfur modifications share common biosynthetic strategies across RNA species, cell types and domains of life. The first section reviews the post-transcriptional sulfur modifications on nucleobases with an emphasis on thiouridine on tRNA and phosphorothioate modification on RNA backbones, as well as the functions of the sulfur modifications on different species of cellular RNAs. The second section reviews the biosynthesis of different types of sulfur modifications and summarizes the general strategy for the biosynthesis of sulfur-containing RNA residues. One of the main goals of investigating sulfur modifications is to aid the genomic drug development pipeline and enhance our understandings of the rapidly growing nucleic acid-based gene therapies. The last section of the review focuses on the current drug development strategies employing sulfur substitution of oxygen atoms in therapeutic RNAs.

Single-molecule optical mapping of the distribution of DNA phosphorothioate epigenetics

DNA phosphorothioate (PT) modifications, with the nonbridging phosphate oxygen replaced by sulfur, governed by DndABCDE or SspABCD, are widely distributed in prokaryotes and have a highly unusual feature of occupying only a small portion of available consensus sequences in a genome. Despite the presence of plentiful non-PT-protected consensuses, DNA PT modification is still employed as a recognition tag by the restriction cognate, for example, DndFGH or SspE, to discriminate and destroy PT-lacking foreign DNA. This raises a fundamental question about how PT modifications are distributed along DNA molecules to keep the restriction components in check. Here, we present two single-molecule strategies that take advantage of the nucleophilicity of PT in combination with fluorescent markers for optical mapping of both single- and double-stranded PT modifications across individual DNA molecules. Surprisingly, PT profiles vary markedly from molecule to molecule, with different PT locations and spacing distances between PT pairs, even in the presence of DndFGH or SspE. The results revealed unprecedented PT modification features previously obscured by ensemble averaging, providing novel insights into the riddles regarding unusual target selection by PT modification and restriction components.

SspABCD-SspFGH Constitutes a New Type of DNA Phosphorothioate-Based Bacterial Defense System

We recently found that SspABCD, catalyzing single-stranded (ss) DNA phosphorothioate (PT) modification, coupled with SspE provides protection against phage infection. SspE performs both PT-simulated NTPase and DNA-nicking nuclease activities to damage phage DNA, rendering SspA-E a PT-sensing defense system.

Unlike nucleobase modifications in canonical restriction-modification systems, DNA phosphorothioate (PT) epigenetic modification occurs in the DNA sugar-phosphate backbone when the nonbridging oxygen is replaced by sulfur in a double-stranded (ds) or single-stranded (ss) manner governed by DndABCDE or SspABCD, respectively. SspABCD coupled with SspE constitutes a defense barrier in which SspE depends on sequence-specific PT modifications to exert its antiphage activity. Here, we identified a new type of ssDNA PT-based SspABCD-SspFGH defense system capable of providing protection against phages through a mode of action different from that of SspABCD-SspE. We provide further evidence that SspFGH damages non-PT-modified DNA and exerts antiphage activity by suppressing phage DNA replication. Despite their different defense mechanisms, SspFGH and SspE are compatible and pair simultaneously with one SspABCD module, greatly enhancing the protection against phages. Together with the observation that the sspBCD-sspFGH cassette is widely distributed in bacterial genomes, this study highlights the diversity of PT-based defense barriers and expands our knowledge of the arsenal of phage defense mechanisms.

Resistance and Adaptation of Bacteria to Non-Antibiotic Antibacterial Agents: Physical Stressors, Nanoparticles, and Bacteriophages

Antimicrobial resistance is a significant threat to human health worldwide, forcing scientists to explore non-traditional antibacterial agents to support rapid interventions and combat the emergence and spread of drug resistant bacteria. Many new antibiotic-free approaches are being developed while the old ones are being revised, resulting in creating unique solutions that arise at the interface of physics, nanotechnology, and microbiology. Specifically, physical factors (e.g., pressure, temperature, UV light) are increasingly used for industrial sterilization. Nanoparticles (unmodified or in combination with toxic compounds) are also applied to circumvent in vivo drug resistance mechanisms in bacteria. Recently, bacteriophage-based treatments are also gaining momentum due to their high bactericidal activity and specificity. Although the number of novel approaches for tackling the antimicrobial resistance crisis is snowballing, it is still unclear if any proposed solutions would provide a long-term remedy. This review aims to provide a detailed overview of how bacteria acquire resistance against these non-antibiotic factors. We also discuss innate bacterial defense systems and how bacteriophages have evolved to tackle them.

Genomic Insights Into the Lifestyles of Thaumarchaeota Inside Sponges

Sponges are among the oldest metazoans and their success is partly due to their abundant and diverse microbial symbionts. They are one of the few animals that have Thaumarchaeota symbionts. Here we compare genomes of 11 Thaumarchaeota sponge symbionts, including three new genomes, to free-living ones. Like their free-living counterparts, sponge-associated Thaumarchaeota can oxidize ammonia, fix carbon, and produce several vitamins. Adaptions to life inside the sponge host include enrichment in transposases, toxin-antitoxin systems and restriction modifications systems, enrichments previously reported also from bacterial sponge symbionts. Most thaumarchaeal sponge symbionts lost the ability to synthesize rhamnose, which likely alters their cell surface and allows them to evade digestion by the host. All but one archaeal sponge symbiont encoded a high-affinity, branched-chain amino acid transporter system that was absent from the analyzed free-living thaumarchaeota suggesting a mixotrophic lifestyle for the sponge symbionts. Most of the other unique features found in sponge-associated Thaumarchaeota, were limited to only a few specific symbionts. These features included the presence of exopolyphosphatases and a glycine cleavage system found in the novel genomes. Thaumarchaeota have thus likely highly specific interactions with their sponge host, which is supported by the limited number of host sponge species to which each of these symbionts is restricted.

Comprehensive classification of ABC ATPases and their functional radiation in nucleoprotein dynamics and biological conflict systems

ABC ATPases form one of the largest clades of P-loop NTPase fold enzymes that catalyze ATP-hydrolysis and utilize its free energy for a staggering range of functions from transport to nucleoprotein dynamics. Using sensitive sequence and structure analysis with comparative genomics, for the first time we provide a comprehensive classification of the ABC ATPase superfamily. ABC ATPases developed structural hallmarks that unambiguously distinguish them from other P-loop NTPases such as an alternative to arginine-finger-based catalysis. At least five and up to eight distinct clades of ABC ATPases are reconstructed as being present in the last universal common ancestor. They underwent distinct phases of structural innovation with the emergence of inserts constituting conserved binding interfaces for proteins or nucleic acids and the adoption of a unique dimeric toroidal configuration for DNA-threading. Specifically, several clades have also extensively radiated in counter-invader conflict systems where they serve as nodal nucleotide-dependent sensory and energetic components regulating a diversity of effectors (including some previously unrecognized) acting independently or together with restriction-modification systems. We present a unified mechanism for ABC ATPase function across disparate systems like RNA editing, translation, metabolism, DNA repair, and biological conflicts, and some unexpected recruitments, such as MutS ATPases in secondary metabolism.

Defense Mechanism of Phosphorothioated DNA under Peroxynitrite-Mediated Oxidative Stress

DNA phosphorothioation (PT) exists in many pathogenic bacteria; however, the mechanism of PT-DNA resistance to the immune response is unclear. In this work, we meticulously investigated the peroxynitrite (PN) tolerance using PT-bioengineered E. coli strains. The in vivo experiment confirms that the S⁺ strain survives better than the S^- strain under moderately oxidative stress. The LCMS, IC, and GCMS experiments demonstrated that phosphorothioate partially converted to phosphate, and the byproduct included sulfate and elemental sulfur. When O,O-diethyl thiophosphate ester (DETP) was used, the reaction rate k₁ was determined to be 4.3 ± 0.5 M^-1 s^-1 in the first-order for both phosphorothioate and peroxynitrite at 35 °C and pH of 8.0. The IC₅₀ values of phosphorothioate dinucleotides are dramatically increased by 400-700-fold compared to DETP. The SH/OH Yin-Yang mechanism rationalizes the in situ DNA self-defense against PN-mediated oxidative stress at the extra bioenergetic cost of DNA modification.

Detection and Quantification of RNA Phosphorothioate Modifications Using Mass Spectrometry

This article describes a protocol for detecting and quantifying RNA phosphorothioate modifications in cellular RNA samples. Starting from solid-phase synthesis of phosphorothioate RNA dinucleotides, followed by purification with reversed-phase HPLC, phosphorothioate RNA dinucleotide standards are prepared for UPLC-MS and LC-MS/MS methods. RNA samples are extracted from cells using TRIzol reagent, then digested with a nuclease mixture and analyzed by mass spectrometry. UPLC-MS is employed first to identify RNA phosphorothioate modifications. An optimized LC-MS/MS method is then employed to quantify the frequency of RNA phosphorothioate modifications in a series of model cells. © 2020 Wiley Periodicals LLC. Basic Protocol 1: Synthesis, purification, and characterization of RNA phosphorothioate dinucleotides Basic Protocol 2: Digestion of RNA samples extracted from cells Basic Protocol 3: Detection and quantification of RNA phosphorothioate modifications by mass spectrometry.

Protein Domain Guided Screen for Sequence Specific and Phosphorothioate-Dependent Restriction Endonucleases

Modification dependent restriction endonucleases (MDREs) restrict modified DNA, typically with limited sequence specificity (∼2-4 bp). Here, we focus on MDREs that have an SRA and/or SBD (sulfur binding domain) fused to an HNH endonuclease domain, cleaving cytosine modified or phosphorothioated (PT) DNA. We independently characterized the SBD-SRA-HNH endonuclease ScoMcrA, which preferentially cleaves 5hmC modified DNA. We report five SBD-HNH endonucleases, all recognizing GpsAAC/GpsTTC sequence and cleaving outside with a single nucleotide 3' stagger: EcoWI (N7/N6), Ksp11411I (N5/N4), Bsp305I (N6/N4-5), Mae9806I [N(8-10)/N(8-9)], and Sau43800I [N(8-9)/N(7-8)]. EcoWI and Bsp305I are more specific for PT modified DNA in Mg2+ buffer, and promiscuous with Mn2+. Ksp11411I is more PT specific with Ni2+. EcoWI and Ksp11411I cleave fully- and hemi-PT modified oligos, while Bsp305I cleaves only fully modified ones. EcoWI forms a dimer in solution and cleaves more efficiently in the presence of two modified sites. In addition, we demonstrate that EcoWI PT-dependent activity has biological function: EcoWI expressing cells restrict dnd+ GpsAAC modified plasmid strongly, and GpsGCC DNA weakly. This work establishes a framework for biotechnology applications of PT-dependent restriction endonucleases (PTDRs).

DNA Phosphorothioate Modifications Are Widely Distributed in the Human Microbiome

The DNA phosphorothioate (PT) modification existing in many prokaryotes, including bacterial pathogens and commensals, confers multiple characteristics, including restricting gene transfer, influencing the global transcriptional response, and reducing fitness during exposure to chemical mediators of inflammation. While PT-containing bacteria have been investigated in a variety of environments, they have not been studied in the human microbiome. Here, we investigated the distribution of PT-harboring strains and verified their existence in the human microbiome. We found over 2000 PT gene-containing strains distributed in different body sites, especially in the gastrointestinal tract. PT-modifying genes are preferentially distributed within several genera, including Pseudomonas, Clostridioides, and Escherichia, with phylogenic diversities. We also assessed the PT modification patterns and found six new PT-linked dinucleotides (C_psG, C_psT, A_psG, T_psG, G_psC, A_psT) in human fecal DNA. To further investigate the PT in the human gut microbiome, we analyzed the abundance of PT-modifying genes and quantified the PT-linked dinucleotides in the fecal DNA. These results confirmed that human microbiome is a rich reservoir for PT-containing microbes and contains a wide variety of PT modification patterns.

ORF4 of the Temperate Archaeal Virus SNJ1 Governs the Lysis-Lysogeny Switch and Superinfection Immunity

Recent environmental and metagenomic studies have considerably increased the repertoire of archaeal viruses and suggested that they play important roles in nutrient cycling in the biosphere. However, very little is known about how they regulate their life cycles and interact with their hosts. Here, we report that the life cycle of the temperate haloarchaeal virus SNJ1 is controlled by the product ORF4, a small protein belonging to the antitoxin MazE superfamily. We show that ORF4 controls the lysis-lysogeny switch of SNJ1 and mediates superinfection immunity by repression of genomic DNA replication of the superinfecting viruses. Bioinformatic analysis shows that ORF4 is highly conserved in two SNJ1-like proviruses, suggesting that the mechanisms for lysis-lysogeny switch and superinfection immunity are conserved in this group of viruses. As the lysis-lysogeny switch and superinfection immunity of archaeal viruses have been poorly studied, we suggest that SNJ1 could serve as a model system to study these processes.IMPORTANCE Archaeal viruses are important parts of the virosphere. Understanding how they regulate their life cycles and interact with host cells provide crucial insights into their biological functions and the evolutionary histories of viruses. However, mechanistic studies of the life cycle of archaeal viruses are scarce due to a lack of genetic tools and demanding cultivation conditions. Here, we discover that the temperate haloarchaeal virus SNJ1, which infects Natrinema sp. strain J7, employs a lysis-lysogeny switch and establishes superinfection immunity like bacteriophages. We show that its ORF4 is critical for both processes and acts as a repressor of the replication of SNJ1. These results establish ORF4 as a master regulator of SNJ1 life cycle and provides novel insights on the regulation of life cycles by temperate archaeal viruses and on their interactions with host cells.

RNA Phosphorothioate Modification in Prokaryotes and Eukaryotes

RNA modifications play important roles in RNA structures and regulation of gene expression and translation. We report the first RNA modification on the phosphate, the RNA phosphorothioate (PS) modification, discovered in both prokaryotes and eukaryotes. The PS modification is also first reported on nucleic acids of eukaryotes. The GpsG modification exists in the Rp configuration and was quantified with liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS). By knocking out the DndA gene in E. coli, we show the Dnd clusters that regulate DNA PS modification may also play roles in RNA PS modification. We also show that the GpsG modification locates on rRNA in E. coli, L. lactis, and HeLa cells, and it is not detected in rRNA-depleted total RNAs from these cells.

Epigenetic competition reveals density-dependent regulation and target site plasticity of phosphorothioate epigenetics in bacteria

Phosphorothioate (PT) DNA modifications-in which a nonbonding phosphate oxygen is replaced with sulfur-represent a widespread, horizontally transferred epigenetic system in prokaryotes and have a highly unusual property of occupying only a small fraction of available consensus sequences in a genome. Using Salmonella enterica as a model, we asked a question of fundamental importance: How do the PT-modifying DndA-E proteins select their G_PSAAC/G_PSTTC targets? Here, we applied innovative analytical, sequencing, and computational tools to discover a novel behavior for DNA-binding proteins: The Dnd proteins are "parked" at the G^6mATC Dam methyltransferase consensus sequence instead of the expected GAAC/GTTC motif, with removal of the ^6mA permitting extensive PT modification of GATC sites. This shift in modification sites further revealed a surprising constancy in the density of PT modifications across the genome. Computational analysis showed that GAAC, GTTC, and GATC share common features of DNA shape, which suggests that PT epigenetics are regulated in a density-dependent manner partly by DNA shape-driven target selection in the genome.

Structural Analysis of an l-Cysteine Desulfurase from an Ssp DNA Phosphorothioation System

Apart from its roles in Fe-S cluster assembly, tRNA thiolation, and sulfur-containing cofactor biosynthesis, cysteine desulfurase serves as a sulfur donor in the DNA PT modification, in which a sulfur atom substitutes a nonbridging oxygen in the DNA phosphodiester backbone. The initial sulfur mobilization from l-cysteine is catalyzed by the SspA cysteine desulfurase in the SspABCD-mediated DNA PT modification system. By determining the crystal structure of SspA, the study presents the molecular mechanism that SspA employs to recognize its cysteine substrate and PLP cofactor. To overcome the long distance (8.9 Å) between the catalytic Cys314 and the cysteine substrate, a conformational change occurs to bring Cys314 to the vicinity of the substrate, allowing for nucleophilic attack.

DNA phosphorothioate (PT) modification, in which the nonbridging oxygen in the sugar-phosphate backbone is substituted by sulfur, is catalyzed by DndABCDE or SspABCD in a double-stranded or single-stranded manner, respectively. In Dnd and Ssp systems, mobilization of sulfur in PT formation starts with the activation of the sulfur atom of cysteine catalyzed by the DndA and SspA cysteine desulfurases, respectively. Despite playing the same biochemical role, SspA cannot be functionally replaced by DndA, indicating its unique physiological properties. In this study, we solved the crystal structure of Vibrio cyclitrophicus SspA in complex with its natural substrate, cysteine, and cofactor, pyridoxal phosphate (PLP), at a resolution of 1.80 Å. Our solved structure revealed the molecular mechanism that SspA employs to recognize its cysteine substrate and PLP cofactor, suggesting a common binding mode shared by cysteine desulfurases. In addition, although the distance between the catalytic Cys314 and the substrate cysteine is 8.9 Å, which is too far for direct interaction, our structural modeling and biochemical analysis revealed a conformational change in the active site region toward the cysteine substrate to move them close to each other to facilitate the nucleophilic attack. Finally, the pulldown analysis showed that SspA could form a complex with SspD, an ATP pyrophosphatase, suggesting that SspD might potentially accept the activated sulfur atom directly from SspA, providing further insights into the biochemical pathway of Ssp-mediated PT modification.

Comparative Genomics of Two New HF1-like Haloviruses

Few genomes of the HF1-group of viruses are currently available, and further examples would enhance the understanding of their evolution, improve their gene annotation, and assist in understanding gene function and regulation. Two novel HF1-group haloviruses, Serpecor1 and Hardycor2, were recovered from widely separated hypersaline lakes in Australia. Both are myoviruses with linear dsDNA genomes and infect the haloarchaeon Halorubrum coriense. Both genomes possess long, terminal direct repeat (TDR) sequences (320 bp for Serpecor1 and 306 bp for Hardycor2). The Serpecor1 genome is 74,196 bp in length, 57.0% G+C, and has 126 annotated coding sequences (CDS). Hardycor2 has a genome of 77,342 bp, 55.6% G+C, and 125 annotated CDS. They show high nucleotide sequence similarity to each other (78%) and with HF1 (>75%), and carry similar intergenic repeat (IR) sequences to those originally described in HF1 and HF2. Hardycor2 carries a DNA methyltransferase gene in the same genomic neighborhood as the methyltransferase genes of HF1, HF2 and HRTV-5, but is in the opposite orientation, and the inferred proteins are only distantly related. Comparative genomics allowed us to identify the candidate genes mediating cell attachment. The genomes of Serpecor1 and Hardycor2 encode numerous small proteins carrying one or more CxxC motifs, a signature feature of zinc-finger domain proteins that are known to participate in diverse biomolecular interactions.

An in vitro DNA phosphorothioate modification reaction

Phosphorothioation (PT) involves the replacement of a nonbridging phosphate oxygen on the DNA backbone with sulfur. In bacteria, the procedure is both sequence- and stereo-specific. We reconstituted the PT reaction using purified DndCDE from Salmonella enterica and IscS from Escherichia coli. We determined that the in vitro process of PT was oxygen sensitive. Only one strand on a double-stranded (ds) DNA substrate was modified in the reaction. The modification was dominant between G and A in the GAAC/GTTC conserved sequence. The modification between G and T required the presence of PT between G and A on the opposite strand. Cysteine, S-adenosyl methionine (SAM) and the formation of an iron-sulfur cluster in DndCDE (DndCDE-FeS) were essential for the process. Results from SAM cleavage reactions support the supposition that PT is a radical SAM reaction. Adenosine triphosphate (ATP) promoted the reaction but was not essential. The data and conclusions presented suggest that the PT reaction in bacteria involves three steps. The first step is the binding of DndCDE-FeS to DNA and searching for the modification sequence, possibly with the help of ATP. Cysteine locks DndCDE-FeS to the modification site with an appropriate protein conformation. SAM triggers the radical SAM reaction to complete the oxygen-sulfur swapping.