In vitro analysis of phosphorothioate modification of DNA reveals substrate recognition by a multiprotein complex

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

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.

Novel Iodine-induced Cleavage Real-time PCR Assay for Accurate Quantification of Phosphorothioate Modified Sites in Bacterial DNA

DNA Phosphorothioate (PT), replacing a non-bridging phosphate oxygen atom with a sulfur atom, is one kind of common DNA modification in bacteria. Whole genome scale description of the location and frequency of PT modification is the key to understand its biological function. Herein we developed a novel method, named with iodine-induced cleavage quantitative real-time PCR (IC-qPCR), to evaluate the frequency of PT modification at a given site in bacterial DNA. The efficiency, dynamic range, sensitivity, reproducibility and accuracy of IC-qPCR were well tested and verified employing an E. coli B7A strain as example. The amplification efficiency of IC-qPCR assay ranged from 91% to 99% with a high correlation coefficient ≥0.99. The limit of quantification was determined as low as 10 copies per reaction for the 607710 and 1818096 sites, and 5 copies for the 302695 and 4120753 sites. Based on the developed IC-qPCR method, the modification frequency of four PTs in E. coli B7A was determined with high accuracy, and the results showed that the PT modification was partial and that the modification frequency varied among investigated PT sites. All these results showed that IC-qPCR was suitable for evaluating the PT modification, which would be helpful to further understand the biological function of PT modification.

Quantitative mapping of DNA phosphorothioatome reveals phosphorothioate heterogeneity of low modification frequency

Phosphorothioate (PT) modifications of the DNA backbone, widespread in prokaryotes, are first identified in bacterial enteropathogens Escherichia coli B7A more than a decade ago. However, methods for high resolution mapping of PT modification level are still lacking. Here, we developed the PT-IC-seq technique, based on iodine-induced selective cleavage at PT sites and high-throughput next generation sequencing, as a mean to quantitatively characterizing the genomic landscape of PT modifications. Using PT-IC-seq we foud that most PT sites are partially modified at a lower PT frequency (< 5%) in E. coli B7A and Salmonella enterica serovar Cerro 87, and both show a heterogeneity pattern of PT modification similar to those of the typical methylation modification. Combining the iodine-induced cleavage and absolute quantification by droplet digital PCR, we developed the PT-IC-ddPCR technique to further measure the PT modification level. Consistent with the PT-IC-seq measurements, PT-IC-ddPCR analysis confirmed the lower PT frequency in E. coli B7A. Our study has demonstrated the heterogeneity of PT modification in the bacterial population and we also established general tools for rigorous mapping and characterization of PT modification events at whole genome level. We describe to our knowledge the first genome-wide quantitative characterization of PT landscape and provides appropriate strategies for further functional studies of PT modification.

Structural basis for the recognition of sulfur in phosphorothioated DNA

There have been very few reports on protein domains that specifically recognize sulfur. Here we present the crystal structure of the sulfur-binding domain (SBD) from the DNA phosphorothioation (PT)-dependent restriction endonuclease ScoMcrA. SBD contains a hydrophobic surface cavity that is formed by the aromatic ring of Y164, the pyrolidine ring of P165, and the non-polar side chains of four other residues that serve as lid, base, and wall of the cavity. The SBD and PT-DNA undergo conformational changes upon binding. The S¹⁸⁷RGRR¹⁹¹ loop inserts into the DNA major groove to make contacts with the bases of the G_PSGCC core sequence. Mutating key residues of SBD impairs PT-DNA association. More than 1000 sequenced microbial species from fourteen phyla contain SBD homologs. We show that three of these homologs bind PT-DNA in vitro and restrict PT-DNA gene transfer in vivo. These results show that SBD-like PT-DNA readers exist widely in prokaryotes.

Solid-Phase Synthesis of Phosphorothioate Oligonucleotides Using Sulfurization Byproducts for in Situ Capping

Oligonucleotides containing phosphorothioate (PS) linkages have recently demonstrated significant clinical utility. PS oligonucleotides are manufactured via a solid-phase chain elongation process in which a four-reaction cycle consisting of detritylation, coupling, sulfurization, and failure sequence capping with Ac₂O is repeated. In the capping step, uncoupled sequences are acetylated at the 5'-OH to stop the chain growth and control the levels of deletion, or ( n-1), impurities. Herein, we report that the byproducts of commonly used sulfurization reagents react with the 5'-OH and cap the failure sequences. The standard Ac₂O capping step can therefore be eliminated, and this 3-reaction cycle process affords a higher yield and higher or comparable overall purity compared to the conventional 4-reaction synthesis. This improvement results in reducing the number of reactions from ∼80 to ∼60 for the synthesis of a typical length 20-mer oligonucleotide. For every kilogram of an oligonucleotide intermediate synthesized, > 500 L of reagents and organic solvents is saved, and the E-factor is decreased to <1500 from ∼2000.

Genome and epigenome of a novel marine Thaumarchaeota strain suggest viral infection, phosphorothioation DNA modification and multiple restriction systems

Marine Thaumarchaeota are abundant ammonia-oxidizers but have few representative laboratory-cultured strains. We report the cultivation of Candidatus Nitrosomarinus catalina SPOT01, a novel strain that is less warm-temperature tolerant than other cultivated Thaumarchaeota. Using metagenomic recruitment, strain SPOT01 comprises a major portion of Thaumarchaeota (4-54%) in temperate Pacific waters. Its complete 1.36 Mbp genome possesses several distinguishing features: putative phosphorothioation (PT) DNA modification genes; a region containing probable viral genes; and putative urea utilization genes. The PT modification genes and an adjacent putative restriction enzyme (RE) operon likely form a restriction modification (RM) system for defence from foreign DNA. PacBio sequencing showed >98% methylation at two motifs, and inferred PT guanine modification of 19% of possible TGCA sites. Metagenomic recruitment also reveals the putative virus region and PT modification and RE genes are present in 18-26%, 9-14% and <1.5% of natural populations at 150 m with ≥85% identity to strain SPOT01. The presence of multiple probable RM systems in a highly streamlined genome suggests a surprising importance for defence from foreign DNA for dilute populations that infrequently encounter viruses or other cells. This new strain provides new insights into the ecology, including viral interactions, of this important group of marine microbes.

Mechanistic Investigation on ROS Resistance of Phosphorothioated DNA

Phosphorothioated DNA (PT-DNA) exhibits a mild anti-oxidant property both in vivo and in vitro. It was found that 8-OHdG and ROS levels were significantly lower in dnd+ (i.e. S⁺) E. coli., compared to a dnd− (i.e. S⁻) strain. Furthermore, different from traditional antioxidants, phosphorothioate compound presents an unexpectedly high capacity to quench hydroxyl radical. Oxidative product analysis by liquid chromatography-mass spectrometry and quantum mechanistic computation supported its unique anti-oxidant characteristic of the hydroxyl selectivity: phosphorothioate donates an electron to either hydroxyl radical or guanine radical derived from hydroxyl radical, leading to a PS^• radical; a complex of PS^• radical and OH⁻ (i.e. the reductive product of hydroxyl radical) releases a highly reductive HS^• radical, which scavenges more equivalents of oxidants in the way to high-covalent sulphur compounds such as sulphur, sulphite and sulphate. The PS-PO conversion (PS and PO denote phosphorus-sulphur and phosphorus-oxygen compounds, respectively) made a switch of extremely oxidative OH^• to highly reductive HS^• species, endowing PT-DNA with the observed high capacity in hydroxyl-radical neutralization. This plausible mechanism provides partial rationale as to why bacteria develop the resource-demanding PT modification on guanine-neighboring phosphates in genome.

DndEi Exhibits Helicase Activity Essential for DNA Phosphorothioate Modification and ATPase Activity Strongly Stimulated by DNA Substrate with a GAAC/GTTC Motif

Phosphorothioate (PT) modification of DNA, in which the non-bridging oxygen of the backbone phosphate group is replaced by sulfur, is governed by the DndA-E proteins in prokaryotes. To better understand the biochemical mechanism of PT modification, functional analysis of the recently found PT-modifying enzyme DndEi, which has an additional domain compared with canonical DndE, from Riemerella anatipestifer is performed in this study. The additional domain is identified as a DNA helicase, and functional deletion of this domain in vivo leads to PT modification deficiency, indicating an essential role of helicase activity in PT modification. Subsequent analysis reveals that the additional domain has an ATPase activity. Intriguingly, the ATPase activity is strongly stimulated by DNA substrate containing a GAAC/GTTC motif (i.e. the motif at which PT modifications occur in R. anatipestifer) when the additional domain and the other domain (homologous to canonical DndE) are co-expressed as a full-length DndEi. These results reveal that PT modification is a biochemical process with DNA strand separation and intense ATP hydrolysis.