Identification of the minimal bacterial 2'-deoxy-7-amido-7-deazaguanine synthesis machinery

Targeted genome mining with GATOR-GC maps the evolutionary landscape of biosynthetic diversity

Gene clusters, groups of physically adjacent genes that work collectively, are pivotal to bacterial fitness and valuable in biotechnology and medicine. While various genome mining tools can identify and characterize gene clusters, they often overlook their evolutionary diversity, a crucial factor in revealing novel cluster functions and applications. To address this gap, we developed GATOR-GC, a targeted genome mining tool that enables comprehensive and flexible exploration of gene clusters in a single execution. We show that GATOR-GC identified a diversity of over 4 million gene clusters similar to experimentally validated biosynthetic gene clusters (BGCs) that other tools fail to detect. To highlight the utility of GATOR-GC, we identified previously uncharacterized co-occurring conserved genes potentially involved in mycosporine-like amino acid biosynthesis and mapped the taxonomic and evolutionary patterns of genomic islands that modify DNA with 7-deazapurines. Additionally, with its proximity-weighted similarity scoring, GATOR-GC successfully differentiated BGCs of the FK-family of metabolites (e.g., rapamycin, FK506/520) according to their chemistries. We anticipate GATOR-GC will be a valuable tool to assess gene cluster diversity for targeted, exploratory, and flexible genome mining. GATOR-GC is available at https://github.com/chevrettelab/gator-gc.

Isolation, Characterization, and Genomic Analysis of Bacteriophages Against Pseudomonas aeruginosa Clinical Isolates from Early and Chronic Cystic Fibrosis Patients for Potential Phage Therapy

Pseudomonas aeruginosa is associated with both community and hospital-acquired infections. It colonizes the lungs of cystic fibrosis (CF) patients, establishing an ecological niche where it adapts and evolves from early to chronic stages, resulting in deteriorating lung function and frequent exacerbations. With antibiotics resistance on the rise, there is a pressing need for alternative personalized treatments (such as bacteriophage therapy) to combat P. aeruginosa infections. In this study, we aimed to isolate and characterize phages targeting both early and chronic P. aeruginosa isolates and evaluate their potential for phage therapy. Four highly virulent phages belonging to myoviral, podviral, and siphoviral morphotypes were isolated from sewage samples. These phages have a broad host range and effectively target 62.5% of the P. aeruginosa isolates with a positive correlation to the early isolates. All the phages have a virulence index of ≥0.90 (0.90-0.98), and one has a large burst size of 331 PFU/cell and a latency period of 30 min. All phages are stable under a wide range of temperature and pH conditions. Genomic analysis suggests the four phages are strictly lytic and devoid of identifiable temperate phage repressors and genes associated with antibiotic resistance and virulence. More significantly, two of the phages significantly delayed the onset of larval death when evaluated in a lethal Galleria mellonella infection model, suggesting their promise as phage therapy candidates for P. aeruginosa infections.

RNA-modification by Base Exchange: Structure, Function and Application of tRNA-guanine Transglycosylases

tRNA-guanine transglycosylases (TGT) occur in all domains of life. They are unique among RNA-modifying enzymes as they exchange a guanine base in the primary RNA transcript by various 7-substituted 7-deazaguanines leading to the modified nucleosides queuosine and archaeosine. Archaeosine is found in the D-loop of archaeal tRNAs, queuosine in the anticodon of bacterial and eukaryotic tRNAs specific for Asp, Asn, His and Tyr. Structural and functional studies revealed a common base-exchange mechanism for all TGTs. Nonetheless, there are also significant differences between TGTs, which will be discussed here. It concerns the specificity for different 7-deazaguanine substrates as well as the recognition of substrate tRNAs. For queuosine TGT an anticodon stem-loop containing the UGU recognition motif is a minimal substrate sufficient for binding to the active site, however, full-length tRNA is bound with higher affinity due to multiple interactions with the dimeric enzyme. Archaeal TGT also binds tRNAs as homodimer, even though the interaction pattern is very different and results in a large change of tRNA conformation. Interestingly, a closely related enzyme, DpdA, exchanges guanine by 7-cyano-7-deazguanine (preQ0) in double stranded DNA of several bacteria. Bacterial TGT is a target for structure-based drug design, as the virulence of Shigella depends on TGT activity, and mammalian TGT has been used for the treatment of murine experimental autoimmune encephalomyelitis, a model for chronic multiple sclerosis. Furthermore, TGT has become a valuable tool in nucleic acid chemistry, as it facilitates the incorporation of non-natural bases in tRNA molecules, e.g. for labelling or cross-linking purposes.

Biosynthesis and function of 7-deazaguanine derivatives in bacteria and phages

SUMMARYDeazaguanine modifications play multifaceted roles in the molecular biology of DNA and tRNA, shaping diverse yet essential biological processes, including the nuanced fine-tuning of translation efficiency and the intricate modulation of codon-anticodon interactions. Beyond their roles in translation, deazaguanine modifications contribute to cellular stress resistance, self-nonself discrimination mechanisms, and host evasion defenses, directly modulating the adaptability of living organisms. Deazaguanine moieties extend beyond nucleic acid modifications, manifesting in the structural diversity of biologically active natural products. Their roles in fundamental cellular processes and their presence in biologically active natural products underscore their versatility and pivotal contributions to the intricate web of molecular interactions within living organisms. Here, we discuss the current understanding of the biosynthesis and multifaceted functions of deazaguanines, shedding light on their diverse and dynamic roles in the molecular landscape of life.

Four additional natural 7-deazaguanine derivatives in phages and how to make them

Bacteriophages and bacteria are engaged in a constant arms race, continually evolving new molecular tools to survive one another. To protect their genomic DNA from restriction enzymes, the most common bacterial defence systems, double-stranded DNA phages have evolved complex modifications that affect all four bases. This study focuses on modifications at position 7 of guanines. Eight derivatives of 7-deazaguanines were identified, including four previously unknown ones: 2'-deoxy-7-(methylamino)methyl-7-deazaguanine (mdPreQ1), 2'-deoxy-7-(formylamino)methyl-7-deazaguanine (fdPreQ1), 2'-deoxy-7-deazaguanine (dDG) and 2'-deoxy-7-carboxy-7-deazaguanine (dCDG). These modifications are inserted in DNA by a guanine transglycosylase named DpdA. Three subfamilies of DpdA had been previously characterized: bDpdA, DpdA1, and DpdA2. Two additional subfamilies were identified in this work: DpdA3, which allows for complete replacement of the guanines, and DpdA4, which is specific to archaeal viruses. Transglycosylases have now been identified in all phages and viruses carrying 7-deazaguanine modifications, indicating that the insertion of these modifications is a post-replication event. Three enzymes were predicted to be involved in the biosynthesis of these newly identified DNA modifications: 7-carboxy-7-deazaguanine decarboxylase (DpdL), dPreQ1 formyltransferase (DpdN) and dPreQ1 methyltransferase (DpdM), which was experimentally validated and harbors a unique fold not previously observed for nucleic acid methylases.

Functional integration of a semi-synthetic azido-queuosine derivative into translation and a tRNA modification circuit

Substitution of the queuine nucleobase precursor preQ1 by an azide-containing derivative (azido-propyl-preQ1) led to incorporation of this clickable chemical entity into tRNA via transglycosylation in vitro as well as in vivo in Escherichia coli, Schizosaccharomyces pombe and human cells. The resulting semi-synthetic RNA modification, here termed Q-L1, was present in tRNAs on actively translating ribosomes, indicating functional integration into aminoacylation and recruitment to the ribosome. The azide moiety of Q-L1 facilitates analytics via click conjugation of a fluorescent dye, or of biotin for affinity purification. Combining the latter with RNAseq showed that TGT maintained its native tRNA substrate specificity in S. pombe cells. The semi-synthetic tRNA modification Q-L1 was also functional in tRNA maturation, in effectively replacing the natural queuosine in its stimulation of further modification of tRNAAsp with 5-methylcytosine at position 38 by the tRNA methyltransferase Dnmt2 in S. pombe. This is the first demonstrated in vivo integration of a synthetic moiety into an RNA modification circuit, where one RNA modification stimulates another. In summary, the scarcity of queuosinylation sites in cellular RNA, makes our synthetic q/Q system a 'minimally invasive' system for placement of a non-natural, clickable nucleobase within the total cellular RNA.

PADLOC: a web server for the identification of antiviral defence systems in microbial genomes

Most bacteria and archaea possess multiple antiviral defence systems that protect against infection by phages, archaeal viruses and mobile genetic elements. Our understanding of the diversity of defence systems has increased greatly in the last few years, and many more systems likely await discovery. To identify defence-related genes, we recently developed the Prokaryotic Antiviral Defence LOCator (PADLOC) bioinformatics tool. To increase the accessibility of PADLOC, we describe here the PADLOC web server (freely available at https://padloc.otago.ac.nz), allowing users to analyse whole genomes, metagenomic contigs, plasmids, phages and archaeal viruses. The web server includes a more than 5-fold increase in defence system types detected (since the first release) and expanded functionality enabling detection of CRISPR arrays and retron ncRNAs. Here, we provide user information such as input options, description of the multiple outputs, limitations and considerations for interpretation of the results, and guidance for subsequent analyses. The PADLOC web server also houses a precomputed database of the defence systems in > 230,000 RefSeq genomes. These data reveal two taxa, Campylobacterota and Spriochaetota, with unusual defence system diversity and abundance. Overall, the PADLOC web server provides a convenient and accessible resource for the detection of antiviral defence systems.

Hypermodified DNA in Viruses of E. coli and Salmonella

The DNA in bacterial viruses collectively contains a rich, yet relatively underexplored, chemical diversity of nucleobases beyond the canonical adenine, guanine, cytosine, and thymine. Herein, we review what is known about the genetic and biochemical basis for the biosynthesis of complex DNA modifications, also called DNA hypermodifications, in the DNA of tailed bacteriophages infecting Escherichia coli and Salmonella enterica. These modifications, and their diversification, likely arose out of the evolutionary arms race between bacteriophages and their cellular hosts. Despite their apparent diversity in chemical structure, the syntheses of various hypermodified bases share some common themes. Hypermodifications form through virus-directed synthesis of noncanonical deoxyribonucleotide triphosphates, direct modification DNA, or a combination of both. Hypermodification enzymes are often encoded in modular operons reminiscent of biosynthetic gene clusters observed in natural product biosynthesis. The study of phage-hypermodified DNA provides an exciting opportunity to expand what is known about the enzyme-catalyzed chemistry of nucleic acids and will yield new tools for the manipulation and interrogation of DNA.

Pantoea Bacteriophage vB_PagS_MED16-A Siphovirus Containing a 2'-Deoxy-7-amido-7-deazaguanosine-Modified DNA

A novel siphovirus, vB_PagS_MED16 (MED16) was isolated in Lithuania using Pantoea agglomerans strain BSL for the phage propagation. The double-stranded DNA genome of MED16 (46,103 bp) contains 73 predicted open reading frames (ORFs) encoding proteins, but no tRNA. Our comparative sequence analysis revealed that 26 of these ORFs code for unique proteins that have no reliable identity when compared to database entries. Based on phylogenetic analysis, MED16 represents a new genus with siphovirus morphology. In total, 35 MED16 ORFs were given a putative functional annotation, including those coding for the proteins responsible for virion morphogenesis, phage-host interactions, and DNA metabolism. In addition, a gene encoding a preQ0 DNA deoxyribosyltransferase (DpdA) is present in the genome of MED16 and the LC-MS/MS analysis indicates 2'-deoxy-7-amido-7-deazaguanosine (dADG)-modified phage DNA, which, to our knowledge, has never been experimentally validated in genomes of Pantoea phages. Thus, the data presented in this study provide new information on Pantoea-infecting viruses and offer novel insights into the diversity of DNA modifications in bacteriophages.

The DUF328 family member YaaA is a DNA-binding protein with a novel fold

DUF328 family proteins are present in many prokaryotes; however, their molecular activities are unknown. The Escherichia coli DUF328 protein YaaA is a member of the OxyR regulon and is protective against oxidative stress. Because uncharacterized proteins involved in prokaryotic oxidative stress response are rare, we sought to learn more about the DUF328 family. Using comparative genomics, we found a robust association between the DUF328 family and genes involved in DNA recombination and the oxidative stress response. In some proteins, DUF328 domains are fused to other domains involved in DNA binding, recombination, and repair. Cofitness analysis indicates that DUF328 family genes associate with recombination-mediated DNA repair pathways, particularly the RecFOR pathway. Purified recombinant YaaA binds to dsDNA, duplex DNA containing bubbles of unpaired nucleotides, and Holliday junction constructs in vitro with dissociation equilibrium constants of 200-300 nm YaaA binds DNA with positive cooperativity, forming multiple shifted species in electrophoretic mobility shift assays. The 1.65-Å resolution X-ray crystal structure of YaaA reveals that the protein possesses a new fold that we name the cantaloupe fold. YaaA has a positively charged cleft and a helix-hairpin-helix DNA-binding motif found in other DNA repair enzymes. Our results demonstrate that YaaA is a new type of DNA-binding protein associated with the oxidative stress response and that this molecular function is likely conserved in other DUF328 family members.

Deoxyinosine and 7-Deaza-2-Deoxyguanosine as Carriers of Genetic Information in the DNA of Campylobacter Viruses

Several reports have demonstrated that Campylobacter bacteriophage DNA is refractory to manipulation, suggesting that these phages encode modified DNA. The characterized Campylobacter jejuni phages fall into two phylogenetic groups within the Myoviridae: the genera Firehammervirus and Fletchervirus Analysis of genomic nucleosides from several of these phages by high-pressure liquid chromatography-mass spectrometry confirmed that 100% of the 2'-deoxyguanosine (dG) residues are replaced by modified bases. Fletcherviruses replace dG with 2'-deoxyinosine, while the firehammerviruses replace dG with 2'-deoxy-7-amido-7-deazaguanosine (dADG), noncanonical nucleotides previously described, but a 100% base substitution has never been observed to have been made in a virus. We analyzed the genome sequences of all available phages representing both groups to elucidate the biosynthetic pathway of these noncanonical bases. Putative ADG biosynthetic genes are encoded by the Firehammervirus phages and functionally complement mutants in the Escherichia coli queuosine pathway, of which ADG is an intermediate. To investigate the mechanism of DNA modification, we isolated nucleotide pools and identified dITP after phage infection, suggesting that this modification is made before nucleotides are incorporated into the phage genome. However, we were unable to observe any form of dADG phosphate, implying a novel mechanism of ADG incorporation into an existing DNA strand. Our results imply that Fletchervirus and Firehammervirus phages have evolved distinct mechanisms to express dG-free DNA.IMPORTANCE Bacteriophages are in a constant evolutionary struggle to overcome their microbial hosts' defenses and must adapt in unconventional ways to remain viable as infectious agents. One mode of adaptation is modifying the viral genome to contain noncanonical nucleotides. Genome modification in phages is becoming more commonly reported as analytical techniques improve, but guanosine modifications have been underreported. To date, two genomic guanosine modifications have been observed in phage genomes, and both are low in genomic abundance. The significance of our research is in the identification of two novel DNA modification systems in Campylobacter-infecting phages, which replace all guanosine bases in the genome in a genus-specific manner.

7-Deazaguanine modifications protect phage DNA from host restriction systems

Genome modifications are central components of the continuous arms race between viruses and their hosts. The archaeosine base (G+), which was thought to be found only in archaeal tRNAs, was recently detected in genomic DNA of Enterobacteria phage 9g and was proposed to protect phage DNA from a wide variety of restriction enzymes. In this study, we identify three additional 2'-deoxy-7-deazaguanine modifications, which are all intermediates of the same pathway, in viruses: 2'-deoxy-7-amido-7-deazaguanine (dADG), 2'-deoxy-7-cyano-7-deazaguanine (dPreQ0) and 2'-deoxy-7- aminomethyl-7-deazaguanine (dPreQ1). We identify 180 phages or archaeal viruses that encode at least one of the enzymes of this pathway with an overrepresentation (60%) of viruses potentially infecting pathogenic microbial hosts. Genetic studies with the Escherichia phage CAjan show that DpdA is essential to insert the 7-deazaguanine base in phage genomic DNA and that 2'-deoxy-7-deazaguanine modifications protect phage DNA from host restriction enzymes.

A Novel Benthic Phage Infecting Shewanella with Strong Replication Ability

The coastal sediments were considered to contain diverse phages playing important roles in driving biogeochemical cycles based on genetic analysis. However, till now, benthic phages in coastal sediments were very rarely isolated, which largely limits our understanding of their biological characteristics. Here, we describe a novel lytic phage (named Shewanella phage S0112) isolated from the coastal sediments of the Yellow Sea infecting a sediment bacterium of the genus Shewanella. The phage has a very high replication capability, with the burst size of ca. 1170 phage particles per infected cell, which is 5–10 times higher than that of most phages isolated before. Meanwhile, the latent period of this phage is relatively longer, which might ensure adequate time for phage replication. The phage has a double-stranded DNA genome comprising 62,286 bp with 102 ORFs, ca. 60% of which are functionally unknown. The expression products of 16 ORF genes, mainly structural proteins, were identified by LC-MS/MS analysis. Besides the general DNA metabolism and structure assembly genes in the phage genome, there is a cluster of auxiliary metabolic genes that may be involved in 7-cyano-7-deazaguanine (preQ₀) biosynthesis. Meanwhile, a pyrophosphohydrolase (MazG) gene being considered as a regulator of programmed cell death or involving in host stringer responses is inserted in this gene cluster. Comparative genomic and phylogenetic analysis both revealed a great novelty of phage S0112. This study represents the first report of a benthic phage infecting Shewanella, which also sheds light on the phage–host interactions in coastal sediments.

Discovery of novel bacterial queuine salvage enzymes and pathways in human pathogens

Queuosine (Q) is a complex tRNA modification widespread in eukaryotes and bacteria that contributes to the efficiency and accuracy of protein synthesis. Eukaryotes are not capable of Q synthesis and rely on salvage of the queuine base (q) as a Q precursor. While many bacteria are capable of Q de novo synthesis, salvage of the prokaryotic Q precursors preQ0 and preQ1 also occurs. With the exception of Escherichia coli YhhQ, shown to transport preQ0 and preQ1, the enzymes and transporters involved in Q salvage and recycling have not been well described. We discovered and characterized 2 Q salvage pathways present in many pathogenic and commensal bacteria. The first, found in the intracellular pathogen Chlamydia trachomatis, uses YhhQ and tRNA guanine transglycosylase (TGT) homologs that have changed substrate specificities to directly salvage q, mimicking the eukaryotic pathway. The second, found in bacteria from the gut flora such as Clostridioides difficile, salvages preQ1 from q through an unprecedented reaction catalyzed by a newly defined subgroup of the radical-SAM enzyme family. The source of q can be external through transport by members of the energy-coupling factor (ECF) family or internal through hydrolysis of Q by a dedicated nucleosidase. This work reinforces the concept that hosts and members of their associated microbiota compete for the salvage of Q precursors micronutrients.

Type II Restriction of Bacteriophage DNA With 5hmdU-Derived Base Modifications

To counteract bacterial defense systems, bacteriophages (phages) make extensive base modifications (substitutions) to block endonuclease restriction. Here we evaluated Type II restriction of three thymidine (T or 5-methyldeoxyuridine, 5mdU) modified phage genomes: Pseudomonas phage M6 with 5-(2-aminoethyl)deoxyuridine (5-NedU), Salmonella phage ViI (Vi1) with 5-(2-aminoethoxy)methyldeoxyuridine (5-NeOmdU) and Delftia phage phi W-14 (a.k.a. ΦW-14) with α-putrescinylthymidine (putT). Among >200 commercially available restriction endonucleases (REases) tested, phage M6, ViI, and phi W-14 genomic DNAs (gDNA) show resistance against 48.4, 71.0, and 68.8% of Type II restrictions, respectively. Inspection of the resistant sites indicates the presence of conserved dinucleotide TG or TC (TS, S=C, or G), implicating the specificity of TS sequence as the target that is converted to modified base in the genomes. We also tested a number of DNA methyltransferases (MTases) on these phage DNAs and found some MTases can fully or partially modify the DNA to confer more resistance to cleavage by REases. Phage M6 restriction fragments can be efficiently ligated by T4 DNA ligase. Phi W-14 restriction fragments show apparent reduced rate in E. coli exonuclease III degradation. This work extends previous studies that hypermodified T derived from 5hmdU provides additional resistance to host-encoded restrictions, in parallel to modified cytosines, guanine, and adenine in phage genomes. The results reported here provide a general guidance to use REases to map and clone phage DNA with hypermodified thymidine.