DNA Phosphorothioate Research Articles

Interrogating the Presence of RNA Phosphorothioate in Nature by a Highly Selective and Sensitive Fluorescence Method.

In nature, DNA phosphorothioate (PT) is found in the genomic materials of some prokaryotes. In contrast, whether there is natural RNA PT is still a question under debate. A groundbreaking study reported the discovery of RNA PT in cellular RNA samples from both prokaryotes and eukaryotes at contents of >100 PT per million nucleotides (PPM-nt) according to liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS) analysis. However, this finding was challenged by a later work showing that other RNA modifications, such as 2'-O-methylation, could give almost the same LC-MS/MS signal patterns as RNA PT. As the LC-MS/MS technique led to contradicting conclusions, another independent method is thus needed to interrogate the presence of RNA PT in nature. In this work, we have developed a highly selective and sensitive fluorescence method for RNA PT quantification based on a new RNA PT-specific conversion reaction. It can detect as low as 2.8 PPM-nt in RNA without interference from RNA thiobases or protein cysteines. We measured the total RNA samples from some bacteria and human cells using this method. None of these samples gave any RNA PT signal above the detection limit (2.8 PPM-nt), suggesting that the widespread presence of natural RNA PT at the 100 PPM-nt level or above is highly unlikely. Nevertheless, due to the limited number of cell species tested in this work, the possible existence of natural RNA PT cannot be excluded. The fluorescence method reported here is simple and low-cost; therefore, it should be an ideal assay for broadly screening various types of cells to search for the clue of RNA PT in nature.

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

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. 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. 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 invitro 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.

Genomic insights into symbiosis and host adaptation of sponge-associated novel bacterium, Rossellomorea orangium sp. nov.

Sponge-associated microorganisms play vital roles in marine sponge ecology. This study presents a genomic investigation of Rossellomorea sp. MCCB 382, isolated from Stelletta sp., reveals insights into its adaptations and symbiotic roles. Phylogenomic study and Overall Genomic Relatedness Index (OGRI) classify MCCB 382 as a novel species, Rossellomorea orangium sp. nov. The genome encodes numerous carbohydrate metabolism enzymes (CAZymes), likely aiding nutrient cycling in the sponge host. Unique eukaryotic-like protein domains hint at potential mechanisms of symbiosis. Defence mechanisms include CRISPR, restriction-modification systems, DNA phosphorothioation, toxin-antitoxin systems, and heavy metal and multidrug resistance genes, indicating adaptation to challenging marine environments. Unlike obligate mutualists, MCCB 382 shows no genome reduction. Furthermore, the presence of mobile genetic elements, horizontal gene transfer, and prophages suggest genetic versatility, implying flexible metabolic potential and capacity for rapid adaptation and symbiosis shifts. MCCB 382 possesses six biosynthetic gene clusters for secondary metabolites, including both type II and III polyketide synthases (PKS), terpenes, (NRPS), NRPS-independent-siderophore, and lassopeptide. Further genome mining using BiGScape revealed four distinct gene cluster families, T2PKS, NRPS-independent-siderophore, lasso peptide, and terpene, presenting opportunities for novel compound elucidation. Our study reveals a symbiotic lifestyle of MCCB 382 with the host sponge, highlighting symbiont factors that aid in establishing and sustaining this relationship. This is the pioneering genomic characterization of a novel Rossellomorea sp. within the sponge Stelletta sp. holobiont.

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.

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.

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.

Understanding base and backbone contributions of phosphorothioate DNA for molecular recognition with SBD proteins.

Bacterial DNA phosphorothioate (PT) modification provides a specific anchoring site for sulfur-binding proteins (SBDs). Besides, their recognition patterns include phosphate links and bases neighboring the PT-modified site, thereby bringing about genome sequence-dependent properties in PT-related epigenetics. Here, we analyze the contributions of the DNA backbone (phosphates and deoxyribose) and bases bound with two SBD proteins in Streptomyces pristinaespiralis and coelicolor (SBDSco and SBDSpr). The chalcogen-hydrophobic interactions remained constantly at the anchoring site while the adjacent bases formed conditional and distinctive non-covalent interactions. More importantly, SBD/PT-DNA interactions were not limited within the traditional "4-bp core" range from 5'-I to 3'-III but extended to upstream 5'-II and 5'-III bases and even 5''-I to 5''-III at the non-PT-modified complementary strand. From the epigenetic viewpoint, bases 3'-II, 5''-I, and 5''-III of SBDSpr and 3'-II, 5''-II, and 5''-III of SBDSco present remarkable differentiations in the molecular recognitions. From the protein viewpoint, H102 in SBDSpr and R191 in SBDSco contribute significantly while proline residues at the PT-bound site are strictly conserved for the PT-chalcogen bond. The mutual and make-up mutations are proposed to alter the SBD/PT-DNA recognition pattern, besides additional chiral phosphorothioate modifications on phosphates 5'-II, 5'-II, 3'-I, and 3'-II.

The functional coupling between restriction and DNA phosphorothioate modification systems underlying the DndFGH restriction complex

The functional coupling between restriction and DNA phosphorothioate modification systems underlying the DndFGH restriction complex

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′-CPSCA-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.

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.

Origin of iodine preferential attack at sulfur in phosphorothioate and subsequent P-O or P-S bond dissociation

Iodine-induced cleavage at phosphorothioate DNA (PT-DNA) is characterized by extremely high sensitivity (∼1 phosphorothioate link per 106 nucleotides), which has been used for detecting and sequencing PT-DNA in bacteria. Despite its foreseeable potential for wide applications, the cleavage mechanism at the PT-modified site has not been well established, and it remains unknown as to whether or not cleavage of the bridging P-O occurs at every PT-modified site. In this work, we conducted accurate ωB97X-D calculations and high-performance liquid chromatography-mass spectrometry to investigate the process of PT-DNA cleavage at the atomic and molecular levels. We have found that iodine chemoselectively binds to the sulfur atom of the phosphorothioate link via a strong halogen-chalcogen interaction (a type of halogen bond, with binding affinity as high as 14.9 kcal/mol) and thus triggers P-O bond cleavage via phosphotriester-like hydrolysis. Additionally, aside from cleavage of the bridging P-O bond, the downstream hydrolyses lead to unwanted P-S/P-O conversions and a loss of the phosphorothioate handle. The mechanism we outline helps to explain specific selectivity at the PT-modified site but also predicts the dynamic stoichiometry of P-S and P-O bond breaking. For instance, Tris is involved in the cascade derivation of S-iodo-phosphorothioate to S-amino-phosphorothioate, suppressing the S-iodo-phosphorothioate hydrolysis to a phosphate diester. However, hydrolysis of one-third of the Tris-O-grafting phosphotriester results in unwanted P-S/P-O conversions. Our study suggests that bacterial DNA phosphorothioation may more frequently occur than previous bioinformatic estimations have predicted from iodine-induced deep sequencing data.

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

ABSTRACTPhosphorothioate (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.

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

ABSTRACTAlthough 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.

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.

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.

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.

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.

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 modification facilitates the dissemination of mcr-1 and blaNDM-1 in drinking water supply systems.

The mechanism driving the dissemination of antibiotic resistance genes (ARGs) in drinking water supply systems (DWSSs) with multiple barriers remains poorly understood despite several recent efforts. Phosphorothioate (PT) modifications, governed by dndABCDE genes, occur naturally in various bacteria and involve the incorporation of sulfur into the DNA backbone. PT is regarded as a mild antioxidant invivo and is known to provide protection against bacterial genomes. We combined quantitative polymerase chain reaction, metagenomic, and network analyses for the water treatment process and laboratory-scale experiments for chlorine treatment using model strains to determine if DNA PT modification occurred in DWSS and facilitated the dissemination of mobilized colistin resistance-1 (mcr-1) and New Delhi metallo-β-lactamase-1 (blaNDM-1) in DWSS. Our results indicated that the relative abundance of dndB increased in the effluent, compared with the influent, in the water treatment plants. Presence of dndB copies had a positive correlation with the concentration of chloramine disinfectant. Network analysis revealed Bdellovibrio as a potential host for MCR genes, NDM genes, and dndB in the DWSS. E.coli DH10B (Wild-type with the dndABCDE gene cluster and ΔdndB) model strains were used to investigate resistance to chlorine treatment at the concentration range of 0.5-3mg/L. The resistance of the wild-type strain increased with increasing concentration of chlorine. DNA PT modification protected MCR- and NDM-carrying bacteria from chloramine disinfection during the water treatment process. The higher relative abundance of ARGs in the effluent of the water treatment plants may be due to the resistance of DNA PT modification to chloramine disinfection, thereby causing the enrichment of genera carrying MCR, NDM, and dndB. This study provides a new understanding on the mechanism of ARG dissemination in DWSS, which will help to improve the performance of drinking water treatment to control the risk associated with antibiotic-resistant bacteria.