Genomic mapping of phosphorothioates reveals partial modification of short consensus sequences

DNA phosphorothioate modification-a new multi-functional epigenetic system in bacteria

Synthetic phosphorothioate (PT) internucleotide linkages, in which a nonbridging oxygen is replaced by a sulphur atom, share similar physical and chemical properties with phosphodiesters but confer enhanced nuclease tolerance on DNA/RNA, making PTs a valuable biochemical and pharmacological tool. Interestingly, PT modification was recently found to occur naturally in bacteria in a sequence-selective and RP configuration-specific manner. This oxygen-sulphur swap is catalysed by the gene products of dndABCDE, which constitute a defence barrier with DndFGH in some bacterial strains that can distinguish and attack non-PT-modified foreign DNA, resembling DNA methylation-based restriction-modification (R-M) systems. Despite their similar defensive mechanisms, PT- and methylation-based R-M systems have evolved to target different consensus contexts in the host cell because when they share the same recognition sequences, the protective function of each can be impeded. The redox and nucleophilic properties of PT sulphur render PT modification a versatile player in the maintenance of cellular redox homeostasis, epigenetic regulation and environmental fitness. The widespread presence of dnd systems is considered a consequence of extensive horizontal gene transfer, whereas the lability of PT during oxidative stress and the susceptibility of PT to PT-dependent endonucleases provide possible explanations for the ubiquitous but sporadic distribution of PT modification in the bacterial world.

Gut microbiome interventions in human health and diseases

Ongoing studies have determined that the gut microbiota is a major factor influencing both health and disease. Host genetic factors and environmental factors contribute to differences in gut microbiota composition and function. Intestinal dysbiosis is a cause or a contributory cause for diseases in multiple body systems, ranging from the digestive system to the immune, cardiovascular, respiratory, and even nervous system. Investigation of pathogenesis has identified specific species or strains, bacterial genes, and metabolites that play roles in certain diseases and represent potential drug targets. As research progresses, gut microbiome-based diagnosis and therapy are proposed and applied, which might lead to considerable progress in precision medicine. We further discuss the limitations of current studies and potential solutions.

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

Archaea and Bacteria have evolved different defence strategies that target virtually all steps of the viral life cycle. The diversified virion morphotypes and genome contents of archaeal viruses result in a highly complex array of archaea-virus interactions. However, our understanding of archaeal antiviral activities lags far behind our knowledges of those in bacteria. Here we report a new archaeal defence system that involves DndCDEA-specific DNA phosphorothioate (PT) modification and the PbeABCD-mediated halt of virus propagation via inhibition of DNA replication. In contrast to the breakage of invasive DNA by DndFGH in bacteria, DndCDEA-PbeABCD does not degrade or cleave viral DNA. The PbeABCD-mediated PT defence system is widespread and exhibits extensive interdomain and intradomain gene transfer events. Our results suggest that DndCDEA-PbeABCD is a new type of PT-based virus resistance system, expanding the known arsenal of defence systems as well as our understanding of host-virus interactions.

Phosphorothioated DNA Is Shielded from Oxidative Damage

DNA is the carrier of genetic information. DNA modifications play a central role in essential physiological processes. Phosphorothioation (PT) modification involves the replacement of an oxygen atom on the DNA backbone with a sulfur atom. PT modification can cause genomic instability in Salmonella enterica under hypochlorous acid stress. This modification restores hydrogen peroxide (H₂O₂) resistance in the catalase-deficient Escherichia coli Hpx^- strain. Here, we report biochemical characterization results for a purified PT modification protein complex (DndCDE) from S. enterica We observed multiplex oligomeric states of DndCDE by using native PAGE. This protein complex bound avidly to PT-modified DNA. DndCDE with an intact iron-sulfur cluster (DndCDE-FeS) possessed H₂O₂ decomposition activity, with a V _max of 10.58 ± 0.90 mM min^-1 and a half-saturation constant, K _0.5S, of 31.03 mM. The Hill coefficient was 2.419 ± 0.59 for this activity. The protein's activity toward H₂O₂ was observed to be dependent on the intact DndCDE and on the formation of an iron-sulfur (Fe-S) cluster on the DndC subunit. In addition to cysteine residues that mediate the formation of this Fe-S cluster, other cysteine residues play a catalytic role. Finally, catalase activity was also detected in DndCDE from Pseudomonas fluorescens Pf0-1. The data and conclusions presented suggest that DndCDE-FeS is a short-lived catalase. Our experiments also indicate that the complex binds to PT sites, shielding PT DNA from H₂O₂ damage. This catalase shield might be able to extend from PT sites to the entire bacterial genome.IMPORTANCE DNA phosphorothioation has been reported in many bacteria. These PT-hosting bacteria live in very different environments, such as the human body, soil, or hot springs. The physiological function of DNA PT modification is still elusive. A remarkable property of PT modification is that purified genomic PT DNA is susceptible to oxidative cleavage. Among the oxidants, hypochlorous acid and H₂O₂ are of physiological relevance for human pathogens since they are generated during the human inflammation response to bacterial infection. However, expression of PT genes in the catalase-deficient E. coli Hpx^- strain restores H₂O₂ resistance. Here, we seek to solve this obvious paradox. We demonstrate that DndCDE-FeS is a short-lived catalase that binds tightly to PT DNA. It is thus possible that by docking to PT sites the catalase activity protects the bacterial genome against H₂O₂ damage.

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.

Advances in CRISPR-Cas systems for RNA targeting, tracking and editing

Clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated protein (Cas) systems, especially type II (Cas9) systems, have been widely used in gene/genome targeting. Modifications of Cas9 enable these systems to become platforms for precise DNA manipulations. However, the utilization of CRISPR-Cas systems in RNA targeting remains preliminary. The discovery of type VI CRISPR-Cas systems (Cas13) shed light on RNA-guided RNA targeting. Cas13d, the smallest Cas13 protein, with a length of only ~930 amino acids, is a promising platform for RNA targeting compatible with viral delivery systems. Much effort has also been made to develop Cas9, Cas13a and Cas13b applications for RNA-guided RNA targeting. The discovery of new RNA-targeting CRISPR-Cas systems as well as the development of RNA-targeting platforms with Cas9 and Cas13 will promote RNA-targeting technology substantially. Here, we review new advances in RNA-targeting CRISPR-Cas systems as well as advances in applications of these systems in RNA targeting, tracking and editing. We also compare these Cas protein-based technologies with traditional technologies for RNA targeting, tracking and editing. Finally, we discuss remaining questions and prospects for the future.

Tight control of genomic phosphorothioate modification by the ATP-modulated autoregulation and reusability of DndB

DNA phosphorothioate (PT) modification was recently identified to occur naturally in diverse bacteria and to be governed by DndABCDE proteins. The nuclease resistance as well as the redox and nucleophilic properties of PT sulfur make PT modification a versatile player in restriction-modification (R-M) defense, epigenetic regulation, environmental fitness and the maintenance of cellular redox homeostasis. In this study, we discovered that tight control of PT levels is mediated by the ATPase activity of DndB. The ATP-binding activity of DndB stimulates the dissociation of the DndB-DNA complex, allowing transcriptional initiation, whereas its ATP hydrolysis activity promotes the conversion of DndB-ATP to free DndB that is capable of rebinding to promoter DNA for transcriptional inhibition. Since sulfur incorporation is an ATP-consuming process, these activities provide an economical way to fine-tune PT modification in an ATP-sensing manner. To our knowledge, this ATP-mediated regulation is a rare example among DNA epigenetic modification systems; the features of autoregulation and the repeated usage of DndB allow the dedicated regulation of PT levels in response to cellular ATP concentrations, providing insight into PT function and its role in physiology.

Engineering and modification of microbial chassis for systems and synthetic biology

Engineering and modifying synthetic microbial chassis is one of the best ways not only to unravel the fundamental principles of life but also to enhance applications in the health, medicine, agricultural, veterinary, and food industries. The two primary strategies for constructing a microbial chassis are the top-down approach (genome reduction) and the bottom-up approach (genome synthesis). Research programs on this topic have been funded in several countries. The 'Minimum genome factory' (MGF) project was launched in 2001 in Japan with the goal of constructing microorganisms with smaller genomes for industrial use. One of the best examples of the results of this project is E. coli MGF-01, which has a reduced-genome size and exhibits better growth and higher threonine production characteristics than the parental strain [1]. The 'cell factory' project was carried out from 1998 to 2002 in the Fifth Framework Program of the EU (European Union), which tried to comprehensively understand microorganisms used in the application field. One of the outstanding results of this project was the elucidation of proteins secreted by Bacillus subtilis, which was summarized as the 'secretome' [2]. The GTL (Genomes to Life) program began in 2002 in the United States. In this program, researchers aimed to create artificial cells both in silico and in vitro, such as the successful design and synthesis of a minimal bacterial genome by John Craig Venter's group [3]. This review provides an update on recent advances in engineering, modification and application of synthetic microbial chassis, with particular emphasis on the value of learning about chassis as a way to better understand life and improve applications.

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.

DiscMLA: An Efficient Discriminative Motif Learning Algorithm over High-Throughput Datasets

The transcription factors (TFs) can activate or suppress gene expression by binding to specific sites, hence are crucial regulatory elements for transcription. Recently, series of discriminative motif finders have been tailored to offering promising strategy for harnessing the power of large quantities of accumulated high-throughput experimental data. However, in order to achieve high speed, these algorithms have to sacrifice accuracy by employing simplified statistical models during the searching process. In this paper, we propose a novel approach named Discriminative Motif Learning via AUC (DiscMLA) to discover motifs on high-throughput datasets. Unlike previous approaches, DiscMLA tries to optimize with a more comprehensive criterion (AUC) during motifs searching. In addition, based on an experimental observation of motif identification on large-scale datasets, some novel procedures are designed to accelerate DiscMLA. The experimental results on 52 real-world datasets demonstrate that our approach substantially outperforms previous methods on discriminative motif learning problems. DiscMLA' stability, discriminability, and validity will help to exploit high-throughput datasets and answer many fundamental biological questions.

Regulation of dndB Gene Expression in Streptomyces lividans

DNA sulfur modification is a unique modification occurring in the sugar-phosphate backbone of DNA, with a nonbridging oxygen atom substituted with sulfur in a sequence-specific and Rp stereo-specific manner. Bioinformatics, RNA-seq, and in vitro transcriptional analyses have shown that DNA sulfur modification may be involved in epigenetic regulation. However, the in vivo evidence supporting this assertion is not convincing. Here, we aimed to characterize two sulfur-modified sites near the dndB promoter region in Streptomyces lividans. Single mutation of either site had no effect on dndB transcription, whereas double mutation of both sites significantly elevated dndB expression. These findings suggested that DNA sulfur modification affected gene expression, and the role of DNA sulfur modification in epigenetic regulation depended on the number of sulfur-modified sites. We also identified an inverted repeat, the R repeat sequence, and showed that this sequence participated in the positive regulation of dndB gene expression.

Genetic mechanisms of arsenic detoxification and metabolism in bacteria

Arsenic, distributed pervasively in the natural environment, is an extremely toxic substance which can severely impair the normal functions of living cells. Research on the genetic mechanisms of arsenic metabolism is of great importance for remediating arsenic-contaminated environments. Many organisms, including bacteria, have developed various strategies to tolerate arsenic, by either detoxifying this harmful element or utilizing it for energy generation. This review summarizes arsenic detoxification as well as arsenic respiratory metabolic pathways in bacteria and discusses novel arsenic resistance pathways in various bacterial strains. This knowledge provides insights into the mechanisms of arsenic biotransformation in bacteria. Multiple detoxification strategies among bacteria imply possible functional relationships among different arsenic detoxification/metabolism pathways. In addition, this review sheds light on the bioremediation of arsenic-contaminated environments and prevention of antibiotic resistance.

Instantaneous Iodine-Assisted DNAzyme Cleavage of Phosphorothioate RNA

Metal ions play a critical role in the RNA-cleavage reaction by interacting with the scissile phosphate and stabilizing the highly negatively charged transition state. Many metal-dependent DNAzymes have been selected for RNA cleavage. Herein, we report that the Ce13d DNAzyme can use nonmetallic iodine (I₂) to cleave a phosphorothioate (PS)-modified substrate. The cleavage yield exceeded 60% for both the R_p and S_p stereoisomers in 10 s, while the yield without the enzyme strand was only ∼10%. The Ce13d cleavage with I₂ also required Na⁺, consistent with the property of Ce13d and confirming the similar role of I₂ as a metal ion. Ce13d had the highest yield among eight tested DNAzymes, with the second highest DNAzyme showing only 20% cleavage. The incomplete cleavage was due to competition from desulfurization and isomerization reactions. This DNAzyme was engineered for fluorescence-based I₂ detection. With EDTA for masking metal ions, I₂ was selectively detected down to 4.7 nM. Oxidation of I^- with Fe³⁺ produced I₂ in situ, allowing detection of Fe³⁺ down to 78 nM. By harnessing nonelectrostatic interactions, such as the I₂/sulfur interaction observed here, more nonmetal species might be discovered to assist DNAzyme-based RNA cleavage.

Identification of a conserved DNA sulfur recognition domain by characterizing the phosphorothioate-specific endonuclease SprMcrA from Streptomyces pristinaespiralis

Streptomyces species have been valuable models for understanding the phenomenon of DNA phosphorothioation in which sulfur replaces a non-bridging oxygen in the phosphate backbone of DNA. We previously reported that the restriction endonuclease ScoMcrA from Streptomyces coelicolor cleaves phosphorothioate DNA and Dcm-methylated DNA at sites 16-28 nucleotides away from the modification sites. However, cleavage of modified DNA by ScoMcrA is always incomplete and accompanied by severe promiscuous activity on unmodified DNA. These features complicate the studies of recognition and cleavage of phosphorothioate DNA. For these reasons, we here characterized SprMcrA from Streptomyces pristinaespiralis, a much smaller homolog of ScoMcrA with a rare HRH motif, a variant of the HNH motif that forms the catalytic center of these endonucleases. The sulfur-binding domain of SprMcrA and its phosphorothioation recognition site were determined. Compared to ScoMcrA, SprMcrA has higher specificity in discerning phosphorothioate DNA from unmodified DNA, and this enzyme generally cuts both strands at a distance of 11-14 nucleotides from the 5' side of the recognition site. The HRH/HNH motif has its own sequence specificity in DNA hydrolysis, leading to failure of cleavage at some phosphorothioated sites. An R248N mutation of the central residue in HRH resulted in 30-fold enhancement in cleavage activity of phosphorothioate DNA and altered the cleavage efficiency at some sites, whereas mutation of both His residues abolished restriction activity. This is the first report of a recognition domain for phosphorothioate DNA and phosphorothioate-dependent and sequence-specific restriction activity.

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.

Signature Arsenic Detoxification Pathways in Halomonas sp. Strain GFAJ-1

Since the original report that Halomonas sp. strain GFAJ-1 was capable of using arsenic instead of phosphorus to sustain growth, additional studies have been conducted, and GFAJ-1 is now considered a highly arsenic-resistant but phosphorus-dependent bacterium. However, the mechanisms supporting the extreme arsenic resistance of the GFAJ-1 strain remain unknown. In this study, we show that GFAJ-1 has multiple distinct arsenic resistance mechanisms. It lacks the genes to reduce arsenate, which is the essential step in the well-characterized resistance mechanism of arsenate reduction coupled to arsenite extrusion. Instead, GFAJ-1 has two arsenic resistance operons, arsH1-acr3-2-arsH2 and mfs1-mfs2-gapdh, enabling tolerance to high levels of arsenate. mfs2 and gapdh encode proteins homologous to Pseudomonas aeruginosa ArsJ and glyceraldehyde-3-phosphate dehydrogenase (GAPDH), respectively, which constitute the equivalent of an As(V) efflux system to catalyze the transformation of inorganic arsenate to pentavalent organoarsenical 1-arseno-3-phosphoglycerate and its subsequent extrusion. Surprisingly, the arsH1-acr3-2-arsH2 operon seems to consist of typical arsenite resistance genes, but this operon is sufficient to confer both arsenite and arsenate resistance on Escherichia coli AW3110 even in the absence of arsenate reductase, suggesting a novel pathway of arsenic detoxification. The simultaneous occurrence of these two unusual detoxification mechanisms enables the adaptation of strain GFAJ-1 to the particularly arsenic-rich environment of Mono Lake.

Halomonas sp. strain GFAJ-1 was previously reported to use arsenic as a substitute for phosphorus to sustain life under phosphate-limited conditions. Although this claim was later undermined by several groups, how GFAJ-1 can thrive in environments with high arsenic concentrations remains unclear. Here, we determined that this ability can be attributed to the possession of two arsenic detoxification operons, arsH1-acr3-2-arsH2 and mfs1-mfs2-gapdh. mfs2 and gapdh encode proteins homologous to ArsJ and GAPDH in Pseudomonas aeruginosa; these proteins create an arsenate efflux pathway to reduce cellular arsenate accumulation. Interestingly, the combination of acr3-2 with either arsH gene was sufficient to confer resistance to both arsenite and arsenate in E. coli AW3110, even in the absence of arsenate reductase, suggesting a new strategy for bacterial arsenic detoxification. This study concludes that the survival of GFAJ-1 in high arsenic concentrations is attributable to the cooccurrence of these two unusual arsenic detoxification mechanisms.

Occurrence, evolution, and functions of DNA phosphorothioate epigenetics in bacteria

The chemical diversity of physiological DNA modifications has expanded with the identification of phosphorothioate (PT) modification in which the nonbridging oxygen in the sugar-phosphate backbone of DNA is replaced by sulfur. Together with DndFGH as cognate restriction enzymes, DNA PT modification, which is catalyzed by the DndABCDE proteins, functions as a bacterial restriction-modification (R-M) system that protects cells against invading foreign DNA. However, the occurrence of dnd systems across a large number of bacterial genomes and their functions other than R-M are poorly understood. Here, a genomic survey revealed the prevalence of bacterial dnd systems: 1,349 bacterial dnd systems were observed to occur sporadically across diverse phylogenetic groups, and nearly half of these occur in the form of a solitary dndBCDE gene cluster that lacks the dndFGH restriction counterparts. A phylogenetic analysis of 734 complete PT R-M pairs revealed the coevolution of M and R components, despite the observation that several PT R-M pairs appeared to be assembled from M and R parts acquired from distantly related organisms. Concurrent epigenomic analysis, transcriptome analysis, and metabolome characterization showed that a solitary PT modification contributed to the overall cellular redox state, the loss of which perturbed the cellular redox balance and induced Pseudomonas fluorescens to reconfigure its metabolism to fend off oxidative stress. An in vitro transcriptional assay revealed altered transcriptional efficiency in the presence of PT DNA modification, implicating its function in epigenetic regulation. These data suggest the versatility of PT in addition to its involvement in R-M protection.

DNA Backbone Sulfur-Modification Expands Microbial Growth Range under Multiple Stresses by its anti-oxidation function

DNA phosphorothioate (PT) modification is a sulfur modification on the backbone of DNA introduced by the proteins DndA-E. It has been detected within many bacteria isolates and metagenomic datasets, including human pathogens, and is considered to be widely distributed in nature. However, little is known about the physiological function of this modification, and thus its evolutionary significance and application potential remains largely a mystery. In this study, we focused on the advantages of DNA PT modification to bacterial cells coping with environmental stresses. We show that the mesophile Escherichia coli and the extremophile Shewanella piezotolerans both expanded their growth ranges following exposure to extreme temperature, salinity, pH, pressure, UV, X-ray and heavy metals as a result of DNA phophorothioation. The phophorothioated DNA reacted to both H₂O₂ and hydroxyl radicals in vivo, and protected genomic DNA as well as sensitive enzymes from intracellular oxidative damage. We further demonstrate that this process has evolved separate from its associated role in DNA restriction and modification. These findings provide a physiological role for a covalent modification widespread in nature and suggest possible applications in biotechnology and biomedicine.

Oxidation of phosphorothioate DNA modifications leads to lethal genomic instability

Genomic modification with sulfur as phosphorothioate (PT) is widespread among prokaryotes, including human pathogens. Apart from its physiological functions, the redox and nucleophilic properties of PT sulfur suggest effects on bacterial fitness in stressful environments. Here we show that PTs are dynamic and labile DNA modifications that cause genomic instability during oxidative stress. Using coupled isotopic labeling-mass spectrometry, we observed sulfur replacement in PTs at a rate of ~2%/h in unstressed Escherichia coli and Salmonella enterica. While PT levels were unaffected by exposure to hydrogen peroxide (H₂O₂) or hypochlorous acid (HOCl), PT turnover increased to 3.8–10%/h for HOCl and was unchanged for H₂O₂, consistent with repair of HOCl-induced sulfur damage. PT-dependent HOCl sensitivity extended to cytotoxicity and DNA strand-breaks, which occurred at orders-of-magnitude lower doses of HOCl than H₂O₂. The genotoxicity of HOCl in PT-containing bacteria suggests reduced fitness in competition with HOCl-producing organisms and during human infections.

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.

Convergence of DNA methylation and phosphorothioation epigenetics in bacterial genomes

Explosive growth in the study of microbial epigenetics has revealed a diversity of chemical structures and biological functions of DNA modifications in restriction-modification (R-M) and basic genetic processes. Here, we describe the discovery of shared consensus sequences for two seemingly unrelated DNA modification systems, ^6mA methylation and phosphorothioation (PT), in which sulfur replaces a nonbridging oxygen in the DNA backbone. Mass spectrometric analysis of DNA from Escherichia coli B7A and Salmonella enterica serovar Cerro 87, strains possessing PT-based R-M genes, revealed d(G_PS^6mA) dinucleotides in the G_PS^6mAAC consensus representing ∼5% of the 1,100 to 1,300 PT-modified d(G_PSA) motifs per genome, with ^6mA arising from a yet-to-be-identified methyltransferase. To further explore PT and ^6mA in another consensus sequence, G_PS^6mATC, we engineered a strain of E. coli HST04 to express Dnd genes from Hahella chejuensis KCTC2396 (PT in G_PSATC) and Dam methyltransferase from E. coli DH10B (^6mA in G^6mATC). Based on this model, in vitro studies revealed reduced Dam activity in G_PSATC-containing oligonucleotides whereas single-molecule real-time sequencing of HST04 DNA revealed ^6mA in all 2,058 G_PSATC sites (5% of 37,698 total GATC sites). This model system also revealed temperature-sensitive restriction by DndFGH in KCTC2396 and B7A, which was exploited to discover that ^6mA can substitute for PT to confer resistance to restriction by the DndFGH system. These results point to complex but unappreciated interactions between DNA modification systems and raise the possibility of coevolution of interacting systems to facilitate the function of each.

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.

DNA Phosphorothioate Modification Plays a Role in Peroxides Resistance in Streptomyces lividans

DNA phosphorothioation, conferred by dnd genes, was originally discovered in the soil-dwelling bacterium Streptomyces lividans, and thereafter found to exist in various bacterial genera. However, the physiological significance of this sulfur modification of the DNA backbone remains unknown in S. lividans. Our studies indicate that DNA phosphorothioation has a major role in resistance to oxidative stress in the strain. Although Streptomyces species express multiple catalase/peroxidase and organic hydroperoxide resistance genes to protect them against peroxide damage, a wild type strain of S. lividans exhibited two-fold to 10-fold higher survival, compared to a dnd⁻ mutant, following treatment with peroxides. RNA-seq experiments revealed that, catalase and organic hydroperoxide resistance gene expression were not up-regulated in the wild type strain, suggesting that the resistance to oxidative stress was not due to the up-regulation of these genes by DNA phosphorothioation. Quantitative RT-PCR analysis was conducted to trace the expression of the catalase and the organic hydroperoxide resistance genes after peroxides treatments. A bunch of these genes were activated in the dnd⁻ mutant rather than the wild type strain in response to peroxides. Moreover, the organic hydroperoxide peracetic acid was scavenged more rapidly in the presence than in the absence of phosphorothioate modification, both in vivo and in vitro. The dnd gene cluster can be up-regulated by the disulfide stressor diamide. Overall, our observations suggest that DNA phosphorothioate modification functions as a peroxide resistance system in S. lividans.

Quantification of total phosphorothioate in bacterial DNA by a bromoimane-based fluorescent method

The discovery of phosphorothioate (PT) modifications in bacterial DNA has challenged our understanding of conserved phosphodiester backbone structure of cellular DNA. This exclusive DNA modification in bacteria is not found in animal cells yet, and its biological function in bacteria is still poorly understood. Quantitative information about the bacterial PT modifications is thus important for the investigation of their possible biological functions. In this study, we have developed a simple fluorescence method for selective quantification of total PTs in bacterial DNA, based on fluorescent labeling of PTs and subsequent release of the labeled fluorophores for absolute quantification. The method was highly selective to PTs and not interfered by the presence of reactive small molecules or proteins. The quantification of PTs in an E. coli DNA sample was successfully achieved using our method and gave a result of about 455 PTs per million DNA nucleotides, while almost no detectable PTs were found in a mammalian calf thymus DNA. With this new method, the content of phosphorothioate in bacterial DNA could be successfully quantified, serving as a simple method suitable for routine use in biological phosphorothioate related studies.

DNA Methylation Assessed by SMRT Sequencing Is Linked to Mutations in Neisseria meningitidis Isolates

The Gram-negative bacterium Neisseria meningitidis features extensive genetic variability. To present, proposed virulence genotypes are also detected in isolates from asymptomatic carriers, indicating more complex mechanisms underlying variable colonization modes of N. meningitidis. We applied the Single Molecule, Real-Time (SMRT) sequencing method from Pacific Biosciences to assess the genome-wide DNA modification profiles of two genetically related N. meningitidis strains, both of serogroup A. The resulting DNA methylomes revealed clear divergences, represented by the detection of shared and of strain-specific DNA methylation target motifs. The positional distribution of these methylated target sites within the genomic sequences displayed clear biases, which suggest a functional role of DNA methylation related to the regulation of genes. DNA methylation in N. meningitidis has a likely underestimated potential for variability, as evidenced by a careful analysis of the ORF status of a panel of confirmed and predicted DNA methyltransferase genes in an extended collection of N. meningitidis strains of serogroup A. Based on high coverage short sequence reads, we find phase variability as a major contributor to the variability in DNA methylation. Taking into account the phase variable loci, the inferred functional status of DNA methyltransferase genes matched the observed methylation profiles. Towards an elucidation of presently incompletely characterized functional consequences of DNA methylation in N. meningitidis, we reveal a prominent colocalization of methylated bases with Single Nucleotide Polymorphisms (SNPs) detected within our genomic sequence collection. As a novel observation we report increased mutability also at 6mA methylated nucleotides, complementing mutational hotspots previously described at 5mC methylated nucleotides. These findings suggest a more diverse role of DNA methylation and Restriction-Modification (RM) systems in the evolution of prokaryotic genomes.

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.

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

A wide variety of prokaryotes possess DNA modifications consisting of sequence-specific phosphorothioates (PT) inserted by members of a five-gene cluster. Recent genome mapping studies revealed two unusual features of PT modifications: short consensus sequences and partial modification of a specific genomic site in a population of bacteria. To better understand the mechanism of target selection of PT modifications that underlies these features, we characterized the substrate recognition of the PT-modifying enzymes termed DptC, D and E in a cell extract system from Salmonella. The results revealed that double-stranded oligodeoxynucleotides underwent de novo PT modification in vitro, with the same modification pattern as in vivo, i. e., G_psAAC/G_psTTC motif. Unexpectedly, in these in vitro analyses we observed no significant effect on PT modification by sequences flanking GAAC/GTTC motif, while PT also occurred in the GAAC/GTTC motif that could not be modified in vivo. Hemi-PT DNA also served as substrate of the PT-modifying enzymes, but not single-stranded DNA. The PT-modifying enzymes were then found to function as a large protein complex, with all of three subunits in tetrameric conformations. This study provided the first demonstration of in vitro DNA PT modification by PT-modifying enzymes that function as a large protein complex.

Regulation of DNA phosphorothioate modification in Salmonella enterica by DndB

DNA phosphorothioate (PT) modification, in which the non-bridging oxygen of the sugar-phosphate backbone is substituted by sulfur, occurs naturally in diverse bacteria and archaea and is regulated by the DndABCDE proteins. DndABCDE and the restriction cognate DndFGHI constitute a new type of defense system that prevents the invasion of foreign DNA in Salmonella enterica serovar Cerro 87. GAAC/GTTC consensus contexts across genomes were found to possess partial PT modifications even in the presence of restriction activity, indicating the regulation of PT. The abundance of PT in cells must be controlled to suit cellular activities. However, the regulatory mechanism of PT modification has not been characterized. The result here indicated that genomic PT modification in S. enterica is controlled by the transcriptional regulator DndB, which binds to two regions in the promoter, each possessing a 5'-TACGN(10)CGTA-3' palindromic motif, to regulate the transcription of dndCDE and its own gene. Site-directed mutagenesis showed that the Cys29 residue of DndB plays a key role in its DNA-binding activity or conformation. Proteomic analysis identified changes to a number of cellular proteins upon up-regulation and loss of PT. Considering the genetic conservation of dnd operons, regulation of PT by DndB might be widespread in diverse organisms.

Regulation of DNA phosphorothioate modifications by the transcriptional regulator DptB in Salmonella

DNA phosphorothioate (PT) modifications, with one non-bridging phosphate oxygen replaced with sulfur, are widely but sporadically distributed in prokaryotic genomes. Short consensus sequences surround the modified linkage in each strain, although each site is only partially modified. The mechanism that maintains this low-frequency modification status is still unknown. In Salmonella enterica serovar Cerro 87, PT modification is mediated by a four-gene cluster called dptBCDE. Here, we found that deletion of dptB led to a significant increase in intracellular PT modification level. In this deletion, transcription of downstream genes was elevated during rapid cell growth. Restoration of dptB on a plasmid restored wild-type levels of expression of downstream genes and PT modification. In vitro, DptB directly protected two separate sequences within the dpt promoter region from DNase I cleavage. Each protected sequence contained a direct repeat (DR). Mutagenesis assays of the DRs demonstrated that each DR was essential for DptB binding. The observation of two shifted species by gel-shift analysis suggests dimer conformation of DptB protein. These DRs are conserved among the promoter regions of dptB homologs, suggesting that this regulatory mechanism is widespread. These findings demonstrate that PT modification is regulated at least in part at the transcriptional level.

Single molecule-level detection and long read-based phasing of epigenetic variations in bacterial methylomes

Beyond its role in host defense, bacterial DNA methylation also plays important roles in the regulation of gene expression, virulence and antibiotic resistance. Bacterial cells in a clonal population can generate epigenetic heterogeneity to increase population-level phenotypic plasticity. Single molecule, real-time (SMRT) sequencing enables the detection of N6-methyladenine and N4-methylcytosine, two major types of DNA modifications comprising the bacterial methylome. However, existing SMRT sequencing-based methods for studying bacterial methylomes rely on a population-level consensus that lacks the single-cell resolution required to observe epigenetic heterogeneity. Here, we present SMALR (single-molecule modification analysis of long reads), a novel framework for single molecule-level detection and phasing of DNA methylation. Using seven bacterial strains, we show that SMALR yields significantly improved resolution and reveals distinct types of epigenetic heterogeneity. SMALR is a powerful new tool that enables de novo detection of epigenetic heterogeneity and empowers investigation of its functions in bacterial populations.

Bacterial DNA methylation is involved in many processes, from host defense to antibiotic resistance, however current methods for examining methylated genomes lack single-cell resolution. Here Beaulaurier et al. present Single Molecule Modification Analysis of Long Reads, a new tool for de novo detection of epigenetic heterogeneity.

The dnd operon for DNA phosphorothioation modification system in Escherichia coli is located in diverse genomic islands

Background

Strains of Escherichia coli that are non-typeable by pulsed-field gel electrophoresis (PFGE) due to in-gel degradation can influence their molecular epidemiological data. The DNA degradation phenotype (Dnd⁺) is mediated by the dnd operon that encode enzymes catalyzing the phosphorothioation of DNA, rendering the modified DNA susceptible to oxidative cleavage during a PFGE run. In this study, a PCR assay was developed to detect the presence of the dnd operon in Dnd⁺E. coli strains and to improve their typeability. Investigations into the genetic environments of the dnd operon in various E. coli strains led to the discovery that the dnd operon is harboured in various diverse genomic islands.

Results

The dndBCDE genes (dnd operon) were detected in all Dnd⁺E. coli strains by PCR. The addition of thiourea improved the typeability of Dnd⁺E. coli strains to 100% using PFGE and the Dnd⁺ phenotype can be observed in both clonal and genetically diverse E. coli strains.

Genomic analysis of 101 dnd operons from genome sequences of Enterobacteriaceae revealed that the dnd operons of the same bacterial species were generally clustered together in the phylogenetic tree. Further analysis of dnd operons of 52 E. coli genomes together with their respective immediate genetic environments revealed a total of 7 types of genetic organizations, all of which were found to be associated with genomic islands designated dnd-encoding GIs. The dnd-encoding GIs displayed mosaic structure and the genomic context of the 7 islands (with 1 representative genome from each type of genetic organization) were also highly variable, suggesting multiple recombination events. This is also the first report where two dnd operons were found within a strain although the biological implication is unknown. Surprisingly, dnd operons were frequently found in pathogenic E. coli although their link with virulence has not been explored.

Conclusion

Genomic islands likely play an important role in facilitating the horizontal gene transfer of the dnd operons in E. coli with 7 different types of islands discovered so far.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1421-8) contains supplementary material, which is available to authorized users.

Crystallization and preliminary X-ray analysis of the type IV restriction endonuclease ScoMcrA from Streptomyces coelicolor, which cleaves both Dcm-methylated DNA and phosphorothioated DNA

ScoMcrA is a type IV modification-dependent restriction endonuclease found in the model strain Streptomyces coelicolor. Unlike type I, II and III restriction endonucleases, which cleave unmodified DNA, type IV restriction endonucleases cleave modified DNA, including methylated, hydroxymethylated, glucosyl-hydroxymethylated and phosphorothioated DNA. ScoMcrA targets both Dcm-methylated DNA and phosphorothioated DNA, and makes double-strand breaks 16-28 nt away from the modified nucleotides or the phosphorothioate links. However, the mechanism by which ScoMcrA recognizes these two entirely different types of modification remains unclear. In this study, the ScoMcrA protein was overexpressed, purified and crystallized. The crystals diffracted to 3.35 Å resolution and belonged to space group P2(1)2(1)2(1). The unit-cell parameters were determined to be a=130.19, b=139.36, c=281.01 Å, α=β=γ=90°. These results will facilitate the detailed structural analysis of ScoMcrA and further elucidation of its biochemical mechanism.

Theoretical study on the relationship between Rp-phosphorothioation and base-step in S-DNA: based on energetic and structural analysis

Phosphorothioation (PT), previously used in synthetic antisense drugs to arrest the transcription or translation process, is also a novel physiological modification in bacteria DNAs. In the previous study, we reported that Rp-phosphorothioation (Rp-PT) destabilizes B-type helix significantly, using a quantum-mechanics-based energy scoring function developed with a dinucleotide model ( Zhang et al. J. Phys. Chem. B , 2012 , 116 , 10639 - 10648 ). A consequent question surfaces in the field of the phosphorothioated DNA (S-DNA) research: does the endogenous chemical modification interact with the base sequence in the bacterial genomes, e.g., in terms of the most common structure of the B-type helix? In this work, we carried out further energetic analysis on the backbone relative energies calculated with the scoring function according to 16 groups of base-step classifications. Moreover, we conducted molecular dynamics simulations of the B-helical structure with the different base-pair steps, to investigate the detailed structural changes upon the O-/S-substitution. As a result, the Rp-PT modification definitively enhances the stiffness of the backbone and differentiates backbone stability as an interaction with base-steps. Furthermore, certain exceptional sequences such as GT and CC were highlighted in the structural analysis of the sulfur local contacts and relative orientation of double strands, indicating that Rp-PT can cross-talk with particular base-steps. The special effects between the phosphorothioation and base-step may be related to the conservative consensus observed highly frequently in bacterial genomes.

DNA phosphorothioate modifications influence the global transcriptional response and protect DNA from double-stranded breaks

The modification of DNA by phosphorothioate (PT) occurs when the non-bridging oxygen in the sugar-phosphate backbone of DNA is replaced with sulfur. This DNA backbone modification was recently discovered and is governed by the dndABCDE genes in a diverse group of bacteria and archaea. However, the biological function of DNA PT modifications is poorly understood. In this study, we employed the RNA-seq analysis to characterize the global transcriptional changes in response to PT modifications. Our results show that DNA without PT protection is susceptible to DNA damage caused by the dndFGHI gene products. The DNA double-stranded breaks then trigger the SOS response, cell filamentation and prophage induction. Heterologous expression of dndBCDE conferring DNA PT modifications at G_PSA and G_PST prevented the damage in Salmonella enterica. Our data provide insights into the physiological role of the DNA PT system.

In vivo mutational characterization of DndE involved in DNA phosphorothioate modification

DNA phosphorothioate (PT) modification is a recently identified epigenetic modification that occurs in the sugar-phosphate backbone of prokaryotic DNA. Previous studies have demonstrated that DNA PT modification is governed by the five DndABCDE proteins in a sequence-selective and R_P stereo-specific manner. Bacteria may have acquired this physiological modification along with dndFGH as a restriction-modification system. However, little is known about the biological function of Dnd proteins, especially the smallest protein, DndE, in the PT modification pathway. DndE was reported to be a DNA-binding protein with a preference for nicked dsDNA in vitro; the binding of DndE to DNA occurs via six positively charged lysine residues on its surface. The substitution of these key lysine residues significantly decreased the DNA binding affinities of DndE proteins to undetectable levels. In this study, we conducted site-directed mutagenesis of dndE on a plasmid and measured DNA PT modifications under physiological conditions by mass spectrometry. We observed distinctive differences from the in vitro binding assays. Several mutants with lysine residues mutated to alanine decreased the total frequency of PT modifications, but none of the mutants completely eliminated PT modification. Our results suggest that the nicked dsDNA-binding capacity of DndE may not be crucial for PT modification and/or that DndE may have other biological functions in addition to binding to dsDNA.

Pathological phenotypes and in vivo DNA cleavage by unrestrained activity of a phosphorothioate-based restriction system in Salmonella

Prokaryotes protect their genomes from foreign DNA with a diversity of defence mechanisms, including a widespread restriction-modification (R-M) system involving phosphorothioate (PT) modification of the DNA backbone. Unlike classical R-M systems, highly partial PT modification of consensus motifs in bacterial genomes suggests an unusual mechanism of PT-dependent restriction. In Salmonella enterica, PT modification is mediated by four genes dptB-E, while restriction involves additional three genes dptF-H. Here, we performed a series of studies to characterize the PT-dependent restriction, and found that it presented several features distinct with traditional R-M systems. The presence of restriction genes in a PT-deficient mutant was not lethal, but instead resulted in several pathological phenotypes. Subsequent transcriptional profiling revealed the expression of > 600 genes was affected by restriction enzymes in cells lacking PT, including induction of bacteriophage, SOS response and DNA repair-related genes. These transcriptional responses are consistent with the observation that restriction enzymes caused extensive DNA cleavage in the absence of PT modifications in vivo. However, overexpression of restriction genes was lethal to the host in spite of the presence PT modifications. These results point to an unusual mechanism of PT-dependent DNA cleavage by restriction enzymes in the face of partial PT modification.