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.

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

ABSTRACT 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. IMPORTANCE 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. To our surprise, ssDNA PT modification can also pair with a newly identified 3-gene sspFGH cassette to fend off phage infection with a different mode of action from that of SspE. Interestingly, both SspFGH and SspE can pair with the same SspABCD module for antiphage defense, and their combination provides Escherichia coli JM109 with additive phage resistance up to 105-fold compared to that for either barrier alone. This agrees with our observation that SspFGH and SspE coexist in 36 bacterial genomes, highlighting the diversity of the gene contents and molecular mechanisms of PT-based defense systems.

Development of Methods Derived from Iodine-Induced Specific Cleavage for Identification and Quantitation of DNA Phosphorothioate Modifications

DNA phosphorothioate (PT) modification is a novel modification that occurs on the DNA backbone, which refers to a non-bridging phosphate oxygen replaced by sulfur. This exclusive DNA modification widely distributes in bacteria but has not been found in eukaryotes to date. PT modification renders DNA nuclease tolerance and serves as a constitute element of bacterial restriction–modification (R–M) defensive system and more biological functions are awaiting exploration. Identification and quantification of the bacterial PT modifications are thus critical to better understanding their biological functions. This work describes three detailed methods derived from iodine-induced specific cleavage-an iodine-induced cleavage assay (ICA), a deep sequencing of iodine-induced cleavage at PT site (ICDS) and an iodine-induced cleavage PT sequencing (PT-IC-Seq)-for the investigation of PT modifications. Using these approaches, we have identified the presence of PT modifications and quantized the frequency of PT modifications in bacteria. These characterizations contributed to the high-resolution genomic mapping of PT modifications, in which the distribution of PT modification sites on the genome was marked accurately and the frequency of the specific modified sites was reliably obtained. Here, we provide time-saving and less labor-consuming methods for both of qualitative and quantitative analysis of genomic PT modifications. The application of these methodologies will offer great potential for better understanding the biology of the PT modifications and open the door to future further systematical study.

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 (CpsG, CpsT, ApsG, TpsG, GpsC, ApsT) 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.

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 GPSAAC/GPSTTC 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 G6mATC 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

ABSTRACT 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. IMPORTANCE 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.

RNA Phosphorothioate Modification in Prokaryotes and Eukaryotes

We report the first RNA phosphorothioate modification in both prokaryotes and eukaryotes. The GpsG modification exists in the Rp configuration and was quantified with liquid chromatography coupled with tandem mass spectrometry. Furthermore, we show the Dnd clusters that regulate DNA PT modification also play roles in RNA PT modification.

RNA Phosphorothioate Modification in Prokaryotes and Eukaryotes

We report the first RNA phosphorothioate modification in both prokaryotes and eukaryotes. The GpsG modification exists in the Rp configuration and was quantified with liquid chromatography coupled with tandem mass spectrometry. Furthermore, we show the Dnd clusters that regulate DNA PT modification also play roles in RNA PT modification.

RNA Phosphorothioate Modification in Prokaryotes and Eukaryotes

We report the first RNA phosphorothioate modification in both prokaryotes and eukaryotes. The GpsG modification exists in the Rp configuration and was quantified with liquid chromatography coupled with tandem mass spectrometry. Furthermore, we show the Dnd clusters that regulate DNA PT modification also play roles in RNA PT modification.