DNA Phosphorothioate Modifications Are Widely Distributed in the Human Microbiome

Overview of Phage Defense Systems in Bacteria and Their Applications

As natural parasites of bacteria, phages have greatly contributed to bacterial evolution owing to their persistent threat. Diverse phage resistance systems have been developed in bacteria during the coevolutionary process with phages. Conversely, phage contamination has a devastating effect on microbial fermentation, resulting in fermentation failure and substantial economic loss. Accordingly, natural defense systems derived from bacteria can be employed to obtain robust phage-resistant host cells that can overcome the threats posed by bacteriophages during industrial bacterial processes. In this review, diverse phage resistance mechanisms, including the remarkable research progress and potential applications, are systematically summarized. In addition, the development prospects and challenges of phage-resistant bacteria are discussed. This review provides a useful reference for developing phage-resistant bacteria.

Phosphorothioate DNA modification by BREX Type 4 systems in the human gut microbiome

Among dozens of microbial DNA modifications regulating gene expression and host defense, phosphorothioation (PT) is the only known backbone modification, with sulfur inserted at a non-bridging oxygen by dnd and ssp gene families. Here we explored the distribution of PT genes in 13,663 human gut microbiome genomes, finding that 6.3% possessed dnd or ssp genes predominantly in Bacillota, Bacteroidota, and Pseudomonadota. This analysis uncovered several putative new PT synthesis systems, including Type 4 Bacteriophage Exclusion (BREX) brx genes, which were genetically validated in Bacteroides salyersiae. Mass spectrometric analysis of DNA from 226 gut microbiome isolates possessing dnd, ssp, and brx genes revealed 8 PT dinucleotide settings confirmed in 6 consensus sequences by PT-specific DNA sequencing. Genomic analysis showed PT enrichment in rRNA genes and depletion at gene boundaries. These results illustrate the power of the microbiome for discovering prokaryotic epigenetics and the widespread distribution of oxidation-sensitive PTs in gut microbes.

PT-seq: A method for metagenomic analysis of phosphorothioate epigenetics in complex microbial communities

Among dozens of known epigenetic marks, naturally occurring phosphorothioate (PT) DNA modifications are unique in replacing a non-bridging phosphate oxygen with redox-active sulfur and function in prokaryotic restriction-modification and transcriptional regulation. Interest in PTs has grown due to the widespread distribution of the dnd, ssp, and brx genes among bacteria and archaea, as well as the discovery of PTs in 5-10% of gut microbes. Efforts to map PTs in complex microbiomes using existing next-generation and direct sequencing technologies have failed due to poor sensitivity. Here we developed PT-seq as a high-sensitivity method to quantitatively map PTs across genomes and metagenomically identify PT-containing microbes in complex genomic mixtures. Like other methods for mapping PTs in individual genomes, PT-seq exploits targeted DNA strand cleavage at PTs by iodine, followed by sequencing library construction using ligation or template switching approaches. However, PT-specific sequencing reads are dramatically increased by adding steps to heat denature the DNA, block pre-existing 3'-ends, fragment DNA after T-tailing, and enrich iodine-induced breaks using biotin-labeling and streptavidin beads capture. Iterative optimization of the sensitivity and specificity of PT-seq is demonstrated with individual bacteria and human fecal DNA.

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

Background

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.

Results

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.

Conclusions

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.

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.

Qiang Huang et al. fed C. elegans with E. coli containing phosphorothioate (PT) DNA or a control strain and evaluated the impact on animal physiology. They observed that worms fed PT( + ) diets exhibited low reactive oxygen species, more active movement, and a longer lifespan compared to controls, suggesting that PT-DNA may have a positive effect on animal health.

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.

Strategies to Avoid Artifacts in Mass Spectrometry-Based Epitranscriptome Analyses

In this report, we perform structure validation of recently reported RNA phosphorothioate (PT) modifications, a new set of epitranscriptome marks found in bacteria and eukaryotes including humans. By comparing synthetic PT-containing diribonucleotides with native species in RNA hydrolysates by high-resolution mass spectrometry (MS), metabolic stable isotope labeling, and PT-specific iodine-desulfurization, we disprove the existence of PTs in RNA from E. coli, S. cerevisiae, human cell lines, and mouse brain. Furthermore, we discuss how an MS artifact led to the initial misidentification of 2'-O-methylated diribonucleotides as RNA phosphorothioates. To aid structure validation of new nucleic acid modifications, we present a detailed guideline for MS analysis of RNA hydrolysates, emphasizing how the chosen RNA hydrolysis protocol can be a decisive factor in discovering and quantifying RNA modifications in biological samples.

Acinetobacter defence mechanisms against biological aggressors and their use as alternative therapeutic applications

Several Acinetobacter strains are important nosocomial pathogens, with Acinetobacter baumannii being the species of greatest worldwide concern due to its multi-drug resistance and the recent appearance of hyper-virulent strains in the clinical setting. Colonisation of this environment is associated with a multitude of bacterial factors, and the molecular features that promote environmental persistence in abiotic surfaces, including intrinsic desiccation resistance, biofilm formation and motility, have been previously addressed. On the contrary, mechanisms enabling Acinetobacter spp. survival when faced against other biological competitors are starting to be characterised. Among them, secretion systems (SS) of different types, such as the T5bSS (Contact-dependent inhibition systems) and the T6SS, confer adaptive advantages against bacterial aggressors. Regarding mechanisms of defence against bacteriophages, such as toxin-antitoxin, restriction-modification, Crispr-Cas and CBASS, among others, have been identified but remain poorly characterised. In view of this, we aimed to summarise the present knowledge on defence mechanisms that enable niche establishment in members of the Acinetobacter genus. Different proposals are also described for the use of some components of these systems as molecular tools to treat Acinetobacter infections.

Nucleic acid binding by SAMHD1 contributes to the antiretroviral activity and is enhanced by the GpsN modification

SAMHD1 impedes infection of myeloid cells and resting T lymphocytes by retroviruses, and the enzymatic activity of the protein-dephosphorylation of deoxynucleotide triphosphates (dNTPs)-implicates enzymatic dNTP depletion in innate antiviral immunity. Here we show that the allosteric binding sites of the enzyme are plastic and can accommodate oligonucleotides in place of the allosteric activators, GTP and dNTP. SAMHD1 displays a preference for oligonucleotides containing phosphorothioate bonds in the Rp configuration located 3' to G nucleotides (GpsN), the modification pattern that occurs in a mechanism of antiviral defense in prokaryotes. In the presence of GTP and dNTPs, binding of GpsN-containing oligonucleotides promotes formation of a distinct tetramer with mixed occupancy of the allosteric sites. Mutations that impair formation of the mixed-occupancy complex abolish the antiretroviral activity of SAMHD1, but not its ability to deplete dNTPs. The findings link nucleic acid binding to the antiretroviral activity of SAMHD1, shed light on the immunomodulatory effects of synthetic phosphorothioated oligonucleotides and raise questions about the role of nucleic acid phosphorothioation in human innate immunity.

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.