In the constant evolutionary battle against mobile genetic elements (MGEs), bacteria have developed several defense mechanisms, some of which target the incoming, foreign nucleic acids e.g. restriction-modification (R-M) or CRISPR-Cas systems. Some of these MGEs, including bacteriophages, have in turn evolved different strategies to evade these hurdles. It was recently shown that the siphophage CAjan and 180 other viruses use 7-deazaguanine modifications in their DNA to evade bacterial R-M systems. Among others, phage CAjan genome contains a gene coding for a DNA-modifying homolog of a tRNA-deazapurine modification enzyme, together with four 7-cyano-7-deazaguanine synthesis genes. Using the CRISPR-Cas9 genome editing tool combined with the Nanopore Sequencing (ONT) we showed that the 7-deazaguanine modification in the CAjan genome is dependent on phage-encoded genes. The modification is also site-specific and is found mainly in two separate DNA sequence contexts: GA and GGC. Homology modeling of the modifying enzyme DpdA provides insight into its probable DNA binding surface and general mode of DNA recognition.

The modified nucleosides 2'-deoxy-7-cyano- and 2'-deoxy-7-amido-7-deazaguanosine (dPreQ0 and dADG, respectively) recently discovered in DNA are the products of the bacterial queuosine tRNA modification pathway and the dpd gene cluster, the latter of which encodes proteins that comprise the elaborate Dpd restriction-modification system present in diverse bacteria. Recent genetic studies implicated the dpdA, dpdB and dpdC genes as encoding proteins necessary for DNA modification, with dpdD-dpdK contributing to the restriction phenotype. Here we report the in vitro reconstitution of the Dpd modification machinery from Salmonella enterica serovar Montevideo, the elucidation of the roles of each protein and the X-ray crystal structure of DpdA supported by small-angle X-ray scattering analysis of DpdA and DpdB, the former bound to DNA. While the homology of DpdA with the tRNA-dependent tRNA-guanine transglycosylase enzymes (TGT) in the queuosine pathway suggested a similar transglycosylase activity responsible for the exchange of a guanine base in the DNA for 7-cyano-7-deazaguanine (preQ0), we demonstrate an unexpected ATPase activity in DpdB necessary for insertion of preQ0 into DNA, and identify several catalytically essential active site residues in DpdA involved in the transglycosylation reaction. Further, we identify a modification site for DpdA activity and demonstrate that DpdC functions independently of DpdA/B in converting preQ0-modified DNA to ADG-modified DNA.

Queuosine (Q) stands out as the sole tRNA modification that can be synthesized via salvage pathways. Comparative genomic analyses identified specific bacteria that showed a discrepancy between the projected Q salvage route and the predicted substrate specificities of the two identified salvage proteins: 1) the distinctive enzyme tRNA guanine-34 transglycosylase (TGT), responsible for inserting precursor bases into target tRNAs; and 2) Queuosine Precursor Transporter (QPTR) , a transporter protein that imports Q precursors. Organisms like the facultative intracellular pathogen Bartonella henselae, which possess only TGT and QPTR but lack predicted enzymes for converting preQ1 to Q, would be expected to salvage the queuine (q) base, mirroring the scenario for the obligate intracellular pathogen Chlamydia trachomatis. However, sequence analyses indicate that the substrate-specificity residues of their TGTs resemble those of enzymes inserting preQ1 rather than q. Intriguingly, mass spectrometry analyses of tRNA modification profiles in B. henselae reveal trace amounts of preQ1, previously not observed in a natural context. Complementation analysis demonstrates that B. henselae TGT and QPTR not only utilize preQ1, akin to their E. coli counterparts, but can also process q when provided at elevated concentrations. The experimental and phylogenomic analyses suggest that the Q pathway in B. henselae could represent an evolutionary transition among intracellular pathogens-from ancestors that synthesized Q de novo to a state prioritizing the salvage of q. Another possibility that will require further investigations is that the insertion of preQ1 has fitness advantages when B. henselae is growing outside a mammalian host.

Gene duplication followed by mutation is a classic mechanism of neofunctionalization, producing gene families with functional diversity. In some cases, a single point mutation is sufficient to change the substrate specificity and/or the chemistry performed by an enzyme, making it difficult to accurately separate enzymes with identical functions from homologs with different functions. Because sequence similarity is often used as a basis for assigning functional annotations to genes, non-isofunctional gene families pose a great challenge for genome annotation pipelines. Here we describe how integrating evolutionary and functional information such as genome context, phylogeny, metabolic reconstruction and signature motifs may be required to correctly annotate multifunctional families. These integrative analyses can also lead to the discovery of novel gene functions, as hints from specific subgroups can guide the functional characterization of other members of the family. We demonstrate how careful manual curation processes using comparative genomics can disambiguate subgroups within large multifunctional families and discover their functions. We present the COG0720 protein family as a case study. We also discuss strategies to automate this process to improve the accuracy of genome functional annotation pipelines.

The universal tRNA modification t6A is found at position 37 of nearly all tRNAs decoding ANN codons. The absence of t6A37 leads to severe growth defects in baker's yeast, phenotypes similar to those caused by defects in mcm5s2U34 synthesis. Mutants in mcm5s2U34 can be suppressed by overexpression of tRNALysUUU, but we show t6A phenotypes could not be suppressed by expressing any individual ANN decoding tRNA, and t6A and mcm5s2U are not determinants for each other's formation. Our results suggest that t6A deficiency, like mcm5s2U deficiency, leads to protein folding defects, and show that the absence of t6A led to stress sensitivities (heat, ethanol, salt) and sensitivity to TOR pathway inhibitors. Additionally, L-homoserine suppressed the slow growth phenotype seen in t6A-deficient strains, and proteins aggregates and Advanced Glycation End-products (AGEs) were increased in the mutants. The global consequences on translation caused by t6A absence were examined by ribosome profiling. Interestingly, the absence of t6A did not lead to global translation defects, but did increase translation initiation at upstream non-AUG codons and increased frame-shifting in specific genes. Analysis of codon occupancy rates suggests that one of the major roles of t6A is to homogenize the process of elongation by slowing the elongation rate at codons decoded by high abundance tRNAs and I34:C3 pairs while increasing the elongation rate of rare tRNAs and G34:U3 pairs. This work reveals that the consequences of t6A absence are complex and multilayered and has set the stage to elucidate the molecular basis of the observed phenotypes.

Gene fusions are the most powerful type of in silico-derived functional associations. However, many fusion compilations were made when <100 genomes were available, and algorithms for identifying fusions need updating to handle the current avalanche of sequenced genomes. The availability of a large fusion dataset would help probe functional associations and enable systematic analysis of where and why fusion events occur.

Tetrahydrofolate (THF) and its one-carbon derivatives, collectively termed folates, are essential cofactors, but are inherently unstable. While it is clear that chemical oxidation can cleave folates or damage their pterin precursors, very little is known about enzymatic damage to these molecules or about whether the folate biosynthesis pathway responds adaptively to damage to its end-products. The presence of a duplication of the gene encoding the folate biosynthesis enzyme 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase (FolK) in many sequenced bacterial genomes combined with a strong chromosomal clustering of the folK gene with panB, encoding the 5,10-methylene-THF-dependent enzyme ketopantoate hydroxymethyltransferase, led us to infer that PanB has a side activity that cleaves 5,10-methylene-THF, yielding a pterin product that is recycled by FolK. Genetic and metabolic analyses of Escherichia coli strains showed that overexpression of PanB leads to accumulation of the likely folate cleavage product 6-hydroxymethylpterin and other pterins in cells and medium, and-unexpectedly-to a 46% increase in total folate content. In silico modeling of the folate biosynthesis pathway showed that these observations are consistent with the in vivo cleavage of 5,10-methylene-THF by a side-activity of PanB, with FolK-mediated recycling of the pterin cleavage product, and with regulation of folate biosynthesis by folates or their damage products.

The emergence of drug-resistant bacteria poses a serious threat to human health. In the case of several antibiotics, including those of the quinolone and rifamycin classes, bacteria rapidly acquire resistance through mutation of chromosomal genes during therapy. In this work, we show that preventing induction of the SOS response by interfering with the activity of the protease LexA renders pathogenic Escherichia coli unable to evolve resistance in vivo to ciprofloxacin or rifampicin, important quinolone and rifamycin antibiotics. We show in vitro that LexA cleavage is induced during RecBC-mediated repair of ciprofloxacin-mediated DNA damage and that this results in the derepression of the SOS-regulated polymerases Pol II, Pol IV and Pol V, which collaborate to induce resistance-conferring mutations. Our findings indicate that the inhibition of mutation could serve as a novel therapeutic strategy to combat the evolution of antibiotic resistance.

The release of the 1000th complete microbial genome will occur in the next two to three years. In anticipation of this milestone, the Fellowship for Interpretation of Genomes (FIG) launched the Project to Annotate 1000 Genomes. The project is built around the principle that the key to improved accuracy in high-throughput annotation technology is to have experts annotate single subsystems over the complete collection of genomes, rather than having an annotation expert attempt to annotate all of the genes in a single genome. Using the subsystems approach, all of the genes implementing the subsystem are analyzed by an expert in that subsystem. An annotation environment was created where populated subsystems are curated and projected to new genomes. A portable notion of a populated subsystem was defined, and tools developed for exchanging and curating these objects. Tools were also developed to resolve conflicts between populated subsystems. The SEED is the first annotation environment that supports this model of annotation. Here, we describe the subsystem approach, and offer the first release of our growing library of populated subsystems. The initial release of data includes 180 177 distinct proteins with 2133 distinct functional roles. This data comes from 173 subsystems and 383 different organisms.

C-1 carriers are essential cofactors in all domains of life, and in Archaea, these can be derivatives of tetrahydromethanopterin (H(4)-MPT) or tetrahydrofolate (H(4)-folate). Their synthesis requires 6-hydroxymethyl-7,8-dihydropterin diphosphate (6-HMDP) as the precursor, but the nature of pathways that lead to its formation were unknown until the recent discovery of the GTP cyclohydrolase IB/MptA family that catalyzes the first step, the conversion of GTP to dihydroneopterin 2',3'-cyclic phosphate or 7,8-dihydroneopterin triphosphate [El Yacoubi, B.; et al. (2006) J. Biol. Chem., 281, 37586-37593 and Grochowski, L. L.; et al. (2007) Biochemistry46, 6658-6667]. Using a combination of comparative genomics analyses, heterologous complementation tests, and in vitro assays, we show that the archaeal protein families COG2098 and COG1634 specify two of the missing 6-HMDP synthesis enzymes. Members of the COG2098 family catalyze the formation of 6-hydroxymethyl-7,8-dihydropterin from 7,8-dihydroneopterin, while members of the COG1634 family catalyze the formation of 6-HMDP from 6-hydroxymethyl-7,8-dihydropterin. The discovery of these missing genes solves a long-standing mystery and provides novel examples of convergent evolutions where proteins of dissimilar architectures perform the same biochemical function.

The reduction of molecular oxygen to water is catalyzed by complicated membrane-bound metallo-enzymes containing variable numbers of subunits, called cytochrome c oxidases or quinol oxidases. We previously described the cytochrome c oxidase II from the hyperthermophilic bacterium Aquifex aeolicus as a ba(3)-type two-subunit (subunits I and II) enzyme and showed that it is included in a supercomplex involved in the sulfide-oxygen respiration pathway. It belongs to the B-family of the heme-copper oxidases, enzymes that are far less studied than the ones from family A. Here, we describe the presence in this enzyme of an additional transmembrane helix "subunit IIa", which is composed of 41 amino acid residues with a measured molecular mass of 5105 Da. Moreover, we show that subunit II, as expected, is in fact longer than the originally annotated protein (from the genome) and contains a transmembrane domain. Using Aquifex aeolicus genomic sequence analyses, N-terminal sequencing, peptide mass fingerprinting and mass spectrometry analysis on entire subunits, we conclude that the B-type enzyme from this bacterium is a three-subunit complex. It is composed of subunit I (encoded by coxA(2)) of 59000 Da, subunit II (encoded by coxB(2)) of 16700 Da and subunit IIa which contain 12, 1 and 1 transmembrane helices respectively. A structural model indicates that the structural organization of the complex strongly resembles that of the ba(3) cytochrome c oxidase from the bacterium Thermus thermophilus, the IIa helical subunit being structurally the lacking N-terminal transmembrane helix of subunit II present in the A-type oxidases. Analysis of the genomic context of genes encoding oxidases indicates that this third subunit is present in many of the bacterial oxidases from B-family, enzymes that have been described as two-subunit complexes.

The EKC/KEOPS complex is universally conserved in Archaea and Eukarya and has been implicated in several cellular processes, including transcription, telomere homeostasis and genomic instability. However, the molecular function of the complex has remained elusive so far. We analyzed the transcriptome of EKC/KEOPS mutants and observed a specific profile that is highly enriched in targets of the Gcn4p transcriptional activator. GCN4 expression was found to be activated at the translational level in mutants via the defective recognition of the inhibitory upstream ORFs (uORFs) present in its leader. We show that EKC/KEOPS mutants are defective for the N6-threonylcarbamoyl adenosine modification at position 37 (t(6)A(37)) of tRNAs decoding ANN codons, which affects initiation at the inhibitory uORFs and provokes Gcn4 de-repression. Structural modeling reveals similarities between Kae1 and bacterial enzymes involved in carbamoylation reactions analogous to t(6)A(37) formation, supporting a direct role for the EKC in tRNA modification. These findings are further supported by strong genetic interactions of EKC mutants with a translation initiation factor and with threonine biosynthesis genes. Overall, our data provide a novel twist to understanding the primary function of the EKC/KEOPS and its impact on several essential cellular functions like transcription and telomere homeostasis.

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.

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.

Members of the DUF34 (domain of unknown function 34) family, also known as the NIF3 protein superfamily, are ubiquitous across superkingdoms. Proteins of this family have been widely annotated as "GTP cyclohydrolase I type 2" through electronic propagation based on one study. Here, the annotation status of this protein family was examined through a comprehensive literature review and integrative bioinformatic analyses that revealed varied pleiotropic associations and phenotypes. This analysis combined with functional complementation studies strongly challenges the current annotation and suggests that DUF34 family members may serve as metal ion insertases, chaperones, or metallocofactor maturases. This general molecular function could explain how DUF34 subgroups participate in highly diversified pathways such as cell differentiation, metal ion homeostasis, pathogen virulence, redox, and universal stress responses.

Over the last 25 years, biology has entered the genomic era and is becoming a science of 'big data'. Most interpretations of genomic analyses rely on accurate functional annotations of the proteins encoded by more than 500 000 genomes sequenced to date. By different estimates, only half the predicted sequenced proteins carry an accurate functional annotation, and this percentage varies drastically between different organismal lineages. Such a large gap in knowledge hampers all aspects of biological enterprise and, thereby, is standing in the way of genomic biology reaching its full potential. A brainstorming meeting to address this issue funded by the National Science Foundation was held during 3-4 February 2022. Bringing together data scientists, biocurators, computational biologists and experimentalists within the same venue allowed for a comprehensive assessment of the current state of functional annotations of protein families. Further, major issues that were obstructing the field were identified and discussed, which ultimately allowed for the proposal of solutions on how to move forward.

Protein synthesis is a complex and highly coordinated process requiring many different protein factors as well as various types of nucleic acids. All translation machinery components require multiple maturation events to be functional. These include post-transcriptional and post-translational modification steps and methylations are the most frequent among these events. In eukaryotes, Trm112, a small protein (COG2835) conserved in all three domains of life, interacts and activates four methyltransferases (Bud23, Trm9, Trm11 and Mtq2) that target different components of the translation machinery (rRNA, tRNAs, release factors). To clarify the function of Trm112 in archaea, we have characterized functionally and structurally its interaction network using Haloferax volcanii as model system. This led us to unravel that methyltransferases are also privileged Trm112 partners in archaea and that this Trm112 network is much more complex than anticipated from eukaryotic studies. Interestingly, among the identified enzymes, some are functionally orthologous to eukaryotic Trm112 partners, emphasizing again the similarity between eukaryotic and archaeal translation machineries. Other partners display some similarities with bacterial methyltransferases, suggesting that Trm112 is a general partner for methyltransferases in all living organisms.

Folates are tripartite molecules comprising pterin, para-aminobenzoate (PABA), and glutamate moieties, which are essential cofactors involved in DNA and amino acid synthesis. The obligately intracellular Chlamydia species have lost several biosynthetic pathways for essential nutrients which they can obtain from their host but have retained the capacity to synthesize folate. In most bacteria, synthesis of the pterin moiety of folate requires the FolEQBK enzymes, while synthesis of the PABA moiety is carried out by the PabABC enzymes. Bioinformatic analyses reveal that while members of Chlamydia are missing the genes for FolE (GTP cyclohydrolase) and FolQ, which catalyze the initial steps in de novo synthesis of the pterin moiety, they have genes for the rest of the pterin pathway. We screened a chlamydial genomic library in deletion mutants of Escherichia coli to identify the "missing genes" and identified a novel enzyme, TrpFCtL2, which has broad substrate specificity. TrpFCtL2, in combination with GTP cyclohydrolase II (RibA), the first enzyme of riboflavin synthesis, provides a bypass of the first two canonical steps in folate synthesis catalyzed by FolE and FolQ. Notably, TrpFCtL2 retains the phosphoribosyl anthranilate isomerase activity of the original annotation. Additionally, we independently confirmed the recent discovery of a novel enzyme, CT610, which uses an unknown precursor to synthesize PABA and complements E. coli mutants with deletions of pabA, pabB, or pabC. Thus, Chlamydia species have evolved a variant folate synthesis pathway that employs a patchwork of promiscuous and adaptable enzymes recruited from other biosynthetic pathways. Importance: Collectively, the involvement of TrpFCtL2 and CT610 in the tetrahydrofolate pathway completes our understanding of folate biosynthesis in Chlamydia. Moreover, the novel roles for TrpFCtL2 and CT610 in the tetrahydrofolate pathway are sophisticated examples of how enzyme evolution plays a vital role in the adaptation of obligately intracellular organisms to host-specific niches. Enzymes like TrpFCtL2 which possess an enzyme fold common to many other enzymes are highly versatile and possess the capacity to evolve to catalyze related reactions in two different metabolic pathways. The continued identification of unique enzymes such as these in bacterial pathogens is important for development of antimicrobial compounds, as drugs that inhibit such enzymes would likely not have any targets in the host or the host's normal microbial flora.

tRNA molecules contain numerous chemically altered nucleosides, which are formed by enzymatic modification of the primary transcripts during the complex tRNA maturation process. Some of the modifications are introduced by single reactions, while other require complex series of reactions carried out by several different enzymes. The location and distribution of various types of modifications vary greatly between different tRNA molecules, organisms and organelles. We have developed a computational method tRNAmodpred, for predicting modifications in tRNA sequences. Briefly, our method takes as an input one or more unmodified tRNA sequences and a set of protein sequences corresponding to a proteome of a cell. Subsequently it identifies homologs of known tRNA modification enzymes in the proteome, predicts tRNA modification activities and maps them onto known pathways of RNA modification from the MODOMICS database. Thereby, theoretically possible modification pathways are identified, and products of these modification reactions are proposed for query tRNAs. This method allows for predicting modification patterns for newly sequenced genomes as well as for checking tentative modification status of tRNAs from one species treated with enzymes from another source, e.g. to predict the possible modifications of eukaryotic tRNAs expressed in bacteria. tRNAmodpred is freely available as a web server at http://genesilico.pl/trnamodpred/.

The obligate-heritable endosymbionts of insects possess some of the smallest known bacterial genomes. This is likely due to loss of genomic material during symbiosis. The mode and rate of this erosion may change over evolutionary time: faster in newly formed associations and slower in long-established ones. The endosymbionts of human and anthropoid primate lice present a unique opportunity to study genome erosion in newly established (or young) symbionts. This is because we have a detailed phylogenetic history of these endosymbionts with divergence dates for closely related species. This allows for genome evolution to be studied in detail and rates of change to be estimated in a phylogenetic framework. Here, we sequenced the genome of the chimpanzee louse endosymbiont (Candidatus Riesia pediculischaeffi) and compared it with the closely related genome of the human body louse endosymbiont. From this comparison, we found evidence for recent genome erosion leading to gene loss in these endosymbionts. Although gene loss was detected, it was not significantly greater than in older endosymbionts from aphids and ants. Additionally, we searched for genes associated with B-vitamin synthesis in the two louse endosymbiont genomes because these endosymbionts are believed to synthesize essential B vitamins absent in the louse's diet. All of the expected genes were present, except those involved in thiamin synthesis. We failed to find genes encoding for proteins involved in the biosynthesis of thiamin or any complete exogenous means of salvaging thiamin, suggesting there is an undescribed mechanism for the salvage of thiamin. Finally, genes encoding for the pantothenate de novo biosynthesis pathway were located on a plasmid in both taxa along with a heat shock protein. Movement of these genes onto a plasmid may be functionally and evolutionarily significant, potentially increasing production and guarding against the deleterious effects of mutation. These data add to a growing resource of obligate endosymbiont genomes and to our understanding of the rate and mode of genome erosion in obligate animal-associated bacteria. Ultimately sequencing additional louse p-endosymbiont genomes will provide a model system for studying genome evolution in obligate host associated bacteria.