Genomic mapping of phosphorothioates reveals partial modification of short consensus sequences

Bacterial phosphorothioate (PT) DNA modifications are incorporated by Dnd proteins A-E and often function with DndF-H as a restriction-modification (R-M) system, as in Escherichia coli B7A. However, bacteria such as Vibrio cyclitrophicus FF75 lack dndF-H, which points to other PT functions. Here we report two novel, orthogonal technologies to map PTs across the genomes of B7A and FF75 with >90% agreement: single molecule, real-time sequencing and deep sequencing of iodine-induced cleavage at PT (ICDS). In B7A, we detect PT on both strands of G[subscript ps]AAC/G[subscript ps]TTC motifs, but with only 12% of 40,701 possible sites modified. In contrast, PT in FF75 occurs as a single-strand modification at C[subscript ps]CA, again with only 14% of 160,541 sites modified. Single-molecule analysis indicates that modification could be partial at any particular genomic site even with active restriction by DndF-H, with direct interaction of modification proteins with GAAC/GTTC sites demonstrated with oligonucleotides. These results point to highly unusual target selection by PT-modification proteins and rule out known R-M mechanisms.National Natural Science Foundation (China)Ministry of Science and Technology of the People's Republic of China (973 and 863 Programs)Shanghai Municipal Council of Science and Technology. Shanghai Pujiang ProgramNational Science Foundation (U.S.) (CHE-1019990)National Institute of Environmental Health Sciences (ES002109)Singapore. National Research Foundation (Singapore-MIT Alliance for Research and Technology