Ribonucleotide reduction - horizontal transfer of a required function spans all three domains

Genome sequences of the first Autographiviridae phages infecting marine Roseobacter

The ubiquitous and abundant marine phages play critical roles in shaping the composition and function of bacterial communities, impacting biogeochemical cycling in marine ecosystems. Autographiviridae is among the most abundant and ubiquitous phage families in the ocean. However, studies on the diversity and ecology of Autographiviridae phages in marine environments are restricted to isolates that infect SAR11 bacteria and cyanobacteria. In this study, ten new roseophages that infect marine Roseobacter strains were isolated from coastal waters. These new roseophages have a genome size ranging from 38 917 to 42 634 bp and G+C content of 44.6-50 %. Comparative genomics showed that they are similar to known Autographiviridae phages regarding gene content and architecture, thus representing the first Autographiviridae roseophages. Phylogenomic analysis based on concatenated conserved genes showed that the ten roseophages form three distinct subgroups within the Autographiviridae, and sequence analysis revealed that they belong to eight new genera. Finally, viromic read-mapping showed that these new Autographiviridae phages are widely distributed in global oceans, mostly inhabiting polar and estuarine locations. This study has expanded the current understanding of the genomic diversity, evolution and ecology of Autographiviridae phages and roseophages. We suggest that Autographiviridae phages play important roles in the mortality and community structure of roseobacters, and have broad ecological applications.

Coordinated transcriptional response to environmental stress by a Synechococcus virus

Viruses are a major control on populations of microbes. Often, their virulence is examined in controlled laboratory conditions. Yet, in nature, environmental conditions lead to changes in host physiology and fitness that may impart both costs and benefits on viral success. Phosphorus (P) is a major abiotic control on the marine cyanobacterium Synechococcus. Some viruses infecting Synechococcus have acquired, from their host, a gene encoding a P substrate binding protein (PstS), thought to improve virus replication under phosphate starvation. Yet, pstS is uncommon among cyanobacterial viruses. Thus, we asked how infections with viruses lacking PstS are affected by P scarcity. We show that the production of infectious virus particles of such viruses is reduced in low P conditions. However, this reduction in progeny is not caused by impaired phage genome replication, thought to be a major sink for cellular phosphate. Instead, transcriptomic analysis showed that under low P conditions, a PstS-lacking cyanophage increased the expression of a specific gene set that included mazG, hli2, and gp43 encoding a pyrophosphatase, a high-light inducible protein and DNA polymerase, respectively. Moreover, several of the upregulated genes were controlled by the host's phoBR two-component system. We hypothesize that recycling and polymerization of nucleotides liberates free phosphate and thus allows viral morphogenesis, albeit at lower rates than when phosphate is replete or when phages encode pstS. Altogether, our data show how phage genomes, lacking obvious P-stress-related genes, have evolved to exploit their host's environmental sensing mechanisms to coordinate their own gene expression in response to resource limitation.

Ubiquitous, B12-dependent virioplankton utilizing ribonucleotide-triphosphate reductase demonstrate interseasonal dynamics and associate with a diverse range of bacterial hosts in the pelagic ocean

Through infection and lysis of their coexisting bacterial hosts, viruses impact the biogeochemical cycles sustaining globally significant pelagic oceanic ecosystems. Currently, little is known of the ecological interactions between lytic viruses and their bacterial hosts underlying these biogeochemical impacts at ecosystem scales. This study focused on populations of lytic viruses carrying the B12-dependent Class II monomeric ribonucleotide reductase (RNR) gene, ribonucleotide-triphosphate reductase (Class II RTPR), documenting seasonal changes in pelagic virioplankton and bacterioplankton using amplicon sequences of Class II RTPR and the 16S rRNA gene, respectively. Amplicon sequence libraries were analyzed using compositional data analysis tools that account for the compositional nature of these data. Both virio- and bacterioplankton communities responded to environmental changes typically seen across seasonal cycles as well as shorter term upwelling-downwelling events. Defining Class II RTPR-carrying viral populations according to major phylogenetic clades proved a more robust means of exploring virioplankton ecology than operational taxonomic units defined by percent sequence homology. Virioplankton Class II RTPR populations showed positive associations with a broad phylogenetic diversity of bacterioplankton including dominant taxa within pelagic oceanic ecosystems such as Prochlorococcus and SAR11. Temporal changes in Class II RTPR virioplankton, occurring as both free viruses and within infected cells, indicated possible viral-host pairs undergoing sustained infection and lysis cycles throughout the seasonal study. Phylogenetic relationships inferred from Class II RTPR sequences mirrored ecological patterns in virio- and bacterioplankton populations demonstrating possible genome to phenome associations for an essential viral replication gene.

Pangenomic Analysis of Nucleo-Cytoplasmic Large DNA Viruses. I: The Phylogenetic Distribution of Conserved Oxygen-Dependent Enzymes Reveals a Capture-Gene Process

The Nucleo-Cytoplasmic Large DNA Viruses (NCLDVs) infect a wide range of eukaryotic species, including amoeba, algae, fish, amphibia, arthropods, birds, and mammals. This group of viruses has linear or circular double-stranded DNA genomes whose size spans approximately one order of magnitude, from 100 to 2500 kbp. The ultimate origin of this peculiar group of viruses remains an open issue. Some have argued that NCLDVs' origin may lie in a bacteriophage ancestor that increased its genome size by subsequent recruitment of eukaryotic and bacterial genes. Others have suggested that NCLDVs families originated from cells that underwent an irreversible process of genome reduction. However, the hypothesis that a number of NCLDVs sequences have been recruited from the host genomes has been largely ignored. In the present work, we have performed pangenomic analyses of each of the seven known NCLDVs families. We show that these families' core- and shell genes have cellular homologs, supporting possible escaping-gene events as part of its evolution. Furthermore, the detection of sequences that belong to two protein families (small chain ribonucleotide reductase and Erv1/Air) and to one superfamily [2OG-Fe(II) oxygenases] that are for distribution in all NCLDVs core and shell clusters encoding for oxygen-dependent enzymes suggests that the highly conserved core these viruses originated after the Proterozoic Great Oxidation Event that transformed the terrestrial atmosphere 2.4-2.3 Ga ago.

Viral Prevalence and Genomic Xenology in the Coevolution of HzNV-2 (Nudiviridae) with Host Helicoverpa zea (Lepidoptera: Noctuidae)

Insect viruses have been described from numerous lineages, yet patterns of genetic exchange and viral prevalence, which are essential to understanding host-virus coevolution, are rarely studied. In Helicoverpa zea, the virus HzNV-2 can cause deformity of male and female genitalia, resulting in sterility. Using ddPCR, we found that male H. zea with malformed genitalia (agonadal) contained high levels of HzNV-2 DNA, confirming previous work. HzNV-2 was found to be prevalent throughout the United States, at more than twice the rate of the baculovirus HaSNPV, and that it contained several host-acquired DNA sequences. HzNV-2 possesses four recently endogenized lepidopteran genes and several more distantly related genes, including one gene with a bacteria-like sequence found in both host and virus. Among the recently acquired genes is cytosolic serine hydroxymethyltransferase (cSHMT). In nearly all tested H. zea, cSHMT contained a 200 bp transposable element (TE) that was not found in cSHMT of the sister species H. armigera. No other virus has been found with host cSHMT, and the study of this shared copy, including possible interactions, may yield new insights into the function of this gene with possible applications to insect biological control, and gene editing.

Comprehensive phylogenetic analysis of the ribonucleotide reductase family reveals an ancestral clade

Ribonucleotide reductases (RNRs) are used by all free-living organisms and many viruses to catalyze an essential step in the de novo biosynthesis of DNA precursors. RNRs are remarkably diverse by primary sequence and cofactor requirement, while sharing a conserved fold and radical-based mechanism for nucleotide reduction. Here, we structurally aligned the diverse RNR family by the conserved catalytic barrel to reconstruct the first large-scale phylogeny consisting of 6779 sequences that unites all extant classes of the RNR family and performed evo-velocity analysis to independently validate our evolutionary model. With a robust phylogeny in-hand, we uncovered a novel, phylogenetically distinct clade that is placed as ancestral to the classes I and II RNRs, which we have termed clade Ø. We employed small-angle X-ray scattering (SAXS), cryogenic-electron microscopy (cryo-EM), and AlphaFold2 to investigate a member of this clade from Synechococcus phage S-CBP4 and report the most minimal RNR architecture to-date. Based on our analyses, we propose an evolutionary model of diversification in the RNR family and delineate how our phylogeny can be used as a roadmap for targeted future study.

Billions of years ago, the Earth’s atmosphere had very little oxygen. It was only after some bacteria and early plants evolved to harness energy from sunlight that oxygen began to fill the Earth’s environment. Oxygen is highly reactive and can interfere with enzymes and other molecules that are essential to life. Organisms living at this point in history therefore had to adapt to survive in this new oxygen-rich world.

An ancient family of enzymes known as ribonucleotide reductases are used by all free-living organisms and many viruses to repair and replicate their DNA. Because of their essential role in managing DNA, these enzymes have been around on Earth for billions of years. Understanding how they evolved could therefore shed light on how nature adapted to increasing oxygen levels and other environmental changes at the molecular level.

One approach to study how proteins evolved is to use computational analysis to construct a phylogenetic tree. This reveals how existing members of a family are related to one another based on the chain of molecules (known as amino acids) that make up each protein. Despite having similar structures and all having the same function, ribonucleotide reductases have remarkably diverse sequences of amino acids. This makes it computationally very demanding to build a phylogenetic tree.

To overcome this, Burnim, Spence, Xu et al. created a phylogenetic tree using structural information from a part of the enzyme that is relatively similar in many modern-day ribonucleotide reductases. The final result took seven continuous months on a supercomputer to generate, and includes over 6,000 members of the enzyme family.

The phylogenetic tree revealed a new distinct group of ribonucleotide reductases that may explain how one adaptation to increasing levels of oxygen emerged in some family members, while another adaptation emerged in others. The approach used in this work also opens up a new way to study how other highly diverse enzymes and other protein families evolved, potentially revealing new insights about our planet’s past.

Physicochemical classification of organisms

The hypervariable residues that compose the major part of proteins’ surfaces are generally considered outside evolutionary control. Yet, these “nonconserved” residues determine the outcome of stochastic encounters in crowded cells. It has recently become apparent that these encounters are not as random as one might imagine, but carefully orchestrated by the intracellular electrostatics to optimize protein diffusion, interactivity, and partner search. The most influential factor here is the protein surface-charge density, which takes different optimal values across organisms with different intracellular conditions. In this study, we examine how far the net-charge density and other physicochemical properties of proteomes will take us in terms of distinguishing organisms in general. The results show that these global proteome properties not only follow the established taxonomical hierarchy, but also provide clues to functional adaptation. In many cases, the proteome–property divergence is even resolved at species level. Accordingly, the variable parts of the genes are not as free to drift as they seem in sequence alignment, but present a complementary tool for functional, taxonomic, and evolutionary assignment.

Newly identified HMO-2011-type phages reveal genomic diversity and biogeographic distributions of this marine viral group

Viruses play critical roles in influencing biogeochemical cycles and adjusting host mortality, population structure, physiology, and evolution in the ocean. Marine viral communities are composed of numerous genetically distinct subfamily/genus-level viral groups. Among currently identified viral groups, the HMO-2011-type group is known to be dominant and broadly distributed. However, only four HMO-2011-type cultivated representatives that infect marine SAR116 and Roseobacter strains have been reported to date, and the genetic diversity, potential hosts, and ecology of this group remain poorly elucidated. Here, we present the genomes of seven HMO-2011-type phages that were isolated using four Roseobacter strains and one SAR11 strain, as well as additional 207 HMO-2011-type metagenomic viral genomes (MVGs) identified from various marine viromes. Phylogenomic and shared-gene analyses revealed that the HMO-2011-type group is a subfamily-level group comprising at least 10 discernible genus-level subgroups. Moreover, >2000 HMO-2011-type DNA polymerase sequences were identified, and the DNA polymerase phylogeny also revealed that the HMO-2011-type group contains diverse subgroups and is globally distributed. Metagenomic read-mapping results further showed that most HMO-2011-type phages are prevalent in global oceans and display distinct geographic distributions, with the distribution of most HMO-2011-type phages being associated with temperature. Lastly, we found that members in subgroup IX, represented by pelagiphage HTVC033P, were among the most abundant HMO-2011-type phages, which implies that SAR11 bacteria are crucial hosts for this viral group. In summary, our findings substantially expand current knowledge regarding the phylogenetic diversity, evolution, and distribution of HMO-2011-type phages, highlighting HMO-2011-type phages as major ecological agents that can infect certain key bacterial groups.

Class-II dihydroorotate dehydrogenases from three phylogenetically distant fungi support anaerobic pyrimidine biosynthesis

Background

In most fungi, quinone-dependent Class-II dihydroorotate dehydrogenases (DHODs) are essential for pyrimidine biosynthesis. Coupling of these Class-II DHODHs to mitochondrial respiration makes their in vivo activity dependent on oxygen availability. Saccharomyces cerevisiae and closely related yeast species harbor a cytosolic Class-I DHOD (Ura1) that uses fumarate as electron acceptor and thereby enables anaerobic pyrimidine synthesis. Here, we investigate DHODs from three fungi (the Neocallimastigomycete Anaeromyces robustus and the yeasts Schizosaccharomyces japonicus and Dekkera bruxellensis) that can grow anaerobically but, based on genome analysis, only harbor a Class-II DHOD.

Results

Heterologous expression of putative Class-II DHOD-encoding genes from fungi capable of anaerobic, pyrimidine-prototrophic growth (Arura9, SjURA9, DbURA9) in an S. cerevisiae ura1Δ strain supported aerobic as well as anaerobic pyrimidine prototrophy. A strain expressing DbURA9 showed delayed anaerobic growth without pyrimidine supplementation. Adapted faster growing DbURA9-expressing strains showed mutations in FUM1, which encodes fumarase. GFP-tagged SjUra9 and DbUra9 were localized to S. cerevisiae mitochondria, while ArUra9, whose sequence lacked a mitochondrial targeting sequence, was localized to the yeast cytosol. Experiments with cell extracts showed that ArUra9 used free FAD and FMN as electron acceptors. Expression of SjURA9 in S. cerevisiae reproducibly led to loss of respiratory competence and mitochondrial DNA. A cysteine residue (C265 in SjUra9) in the active sites of all three anaerobically active Ura9 orthologs was shown to be essential for anaerobic activity of SjUra9 but not of ArUra9.

Conclusions

Activity of fungal Class-II DHODs was long thought to be dependent on an active respiratory chain, which in most fungi requires the presence of oxygen. By heterologous expression experiments in S. cerevisiae, this study shows that phylogenetically distant fungi independently evolved Class-II dihydroorotate dehydrogenases that enable anaerobic pyrimidine biosynthesis. Further structure–function studies are required to understand the mechanistic basis for the anaerobic activity of Class-II DHODs and an observed loss of respiratory competence in S. cerevisiae strains expressing an anaerobically active DHOD from Sch. japonicus.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40694-021-00117-4.

Comparative genomic analysis reveals an 'open' pan-genome of African swine fever virus

The worldwide transmission of African swine fever virus (ASFV) drastically affects the pig industry and global trade. Development of vaccines is hindered by the lack of knowledge of the genomic characteristics of ASFV. In this study, we developed a pipeline for the de novo assembly of ASFV genome without virus isolation and purification. We then used a comparative genomics approach to systematically study 46 genomes of ASFVs to reveal the genomic characteristics. The analysis revealed that ASFV has an 'open' pan-genome based on both protein-coding genes and intergenic regions. Of the 151-174 genes found in the ASFV strains, only 86 were identified as core genes; the remainder were flexible accessory genes. Notably, 44 of the 86 core genes and 155 of the 324 accessory genes have been functionally annotated according to the known proteins. Interestingly, a dynamic number of taxis-related genes were identified in the accessory genes, and two potential virulence genes were identified in all ASFV isolates. The 'open' pan-genome of ASFV based on gene and intergenic regions reveals its pronounced natural diversity concerning genomic composition and regulation.

Non-host class II ribonucleotide reductase in Thermus viruses: sequence adaptation and host interaction

Ribonucleotide reductases (RNR) are essential enzymes for all known life forms. Their current taxonomic distribution suggests extensive horizontal gene transfer e.g., by processes involving viruses. To improve our understanding of the underlying processes, we characterized a monomeric class II RNR (NrdJm) enzyme from a Thermus virus, a subclass not present in any sequenced Thermus spp. genome. Phylogenetic analysis revealed a distant origin of the nrdJm gene with the most closely related sequences found in mesophiles or moderate thermophiles from the Firmicutes phylum. GC-content, codon usage and the ratio of coding to non-coding substitutions (dN/dS) suggest extensive adaptation of the gene in the virus in terms of nucleotide composition and amino acid sequence. The NrdJm enzyme is a monomeric B₁₂-dependent RNR with nucleoside triphosphate specificity. It exhibits a temperature optimum at 60–70 °C, which is in the range of the growth optimum of Thermus spp. Experiments in combination with the Thermus thermophilus thioredoxin system show that the enzyme is able to retrieve electrons from the host NADPH pool via host thioredoxin and thioredoxin reductases. This is different from other characterized viral RNRs such as T4 phage RNR, where a viral thioredoxin is present. We hence show that the monomeric class II RNR, present in Thermus viruses, was likely transferred from an organism phylogenetically distant from the one they were isolated from, and adapted to the new host in genetic signature and amino acids sequence.

Reannotation of the Ribonucleotide Reductase in a Cyanophage Reveals Life History Strategies Within the Virioplankton

Ribonucleotide reductases (RNRs) are ancient enzymes that catalyze the reduction of ribonucleotides to deoxyribonucleotides. They are required for virtually all cellular life and are prominent within viral genomes. RNRs share a common ancestor and must generate a protein radical for direct ribonucleotide reduction. The mechanisms by which RNRs produce radicals are diverse and divide RNRs into three major classes and several subclasses. The diversity of radical generation methods means that cellular organisms and viruses typically contain the RNR best-suited to the environmental conditions surrounding DNA replication. However, such diversity has also fostered high rates of RNR misannotation within subject sequence databases. These misannotations have resulted in incorrect translative presumptions of RNR biochemistry and have diminished the utility of this marker gene for ecological studies of viruses. We discovered a misannotation of the RNR gene within the Prochlorococcus phage P-SSP7 genome, which caused a chain of misannotations within commonly observed RNR genes from marine virioplankton communities. These RNRs are found in marine cyanopodo- and cyanosiphoviruses and are currently misannotated as Class II RNRs, which are O2-independent and require cofactor B12. In fact, these cyanoviral RNRs are Class I enzymes that are O2-dependent and may require a di-metal cofactor made of Fe, Mn, or a combination of the two metals. The discovery of an overlooked Class I β subunit in the P-SSP7 genome, together with phylogenetic analysis of the α and β subunits confirms that the RNR from P-SSP7 is a Class I RNR. Phylogenetic and conserved residue analyses also suggest that the P-SSP7 RNR may constitute a novel Class I subclass. The reannotation of the RNR clade represented by P-SSP7 means that most lytic cyanophage contain Class I RNRs, while their hosts, B12-producing Synechococcus and Prochlorococcus, contain Class II RNRs. By using a Class I RNR, cyanophage avoid a dependence on host-produced B12, a more effective strategy for a lytic virus. The discovery of a novel RNR β subunit within cyanopodoviruses also implies that some unknown viral genes may be familiar cellular genes that are too divergent for homology-based annotation methods to identify.

Hydroxylamine Derivatives as a New Paradigm in the Search of Antibacterial Agents

Serious infections caused by bacteria that are resistant to commonly used antibiotics have become a major global healthcare problem in the 21st century. Multidrug-resistant bacteria causing severe infections mainly grow in complex bacterial communities known as biofilms, in which bacterial resistance to antibacterial agents and to the host immune system is strengthened. As drug resistance is becoming a threatening problem, it is necessary to develop new antimicrobial agents with novel mechanisms of action. Here, we designed and synthesized a small library of N-substituted hydroxylamine (N-HA) compounds with antibacterial activity. These compounds, acting as radical scavengers, inhibit the bacterial ribonucleotide reductase (RNR) enzyme. RNR enzyme is essential for bacterial proliferation during infection, as it provides the building blocks for DNA synthesis and repair. We demonstrate the broad antimicrobial effect of several drug candidates against a variety of Gram-positive and Gram-negative bacteria, together with low toxicity toward eukaryotic cells. Furthermore, the most promising compounds can reduce the biomass of an established biofilm on Pseudomonas aeruginosa, Staphylococcus aureus, and Escherichia coli. This study settles the starting point to develop new N-hydroxylamine compounds as potential effective antibacterial agents to fight against drug-resistant pathogenic bacteria.

Family A DNA Polymerase Phylogeny Uncovers Diversity and Replication Gene Organization in the Virioplankton

Shotgun metagenomics, which allows for broad sampling of viral diversity, has uncovered genes that are widely distributed among virioplankton populations and show linkages to important biological features of unknown viruses. Over 25% of known dsDNA phage carry the DNA polymerase I (polA) gene, making it one of the most widely distributed phage genes. Because of its pivotal role in DNA replication, this enzyme is linked to phage lifecycle characteristics. Previous research has suggested that a single amino acid substitution might be predictive of viral lifestyle. In this study Chesapeake Bay virioplankton were sampled by shotgun metagenomic sequencing (using long and short read technologies). More polA sequences were predicted from this single viral metagenome (virome) than from 86 globally distributed virome libraries (ca. 2,100, and 1,200, respectively). The PolA peptides predicted from the Chesapeake Bay virome clustered with 69% of PolA peptides from global viromes; thus, remarkably the Chesapeake Bay virome captured the majority of known PolA peptide diversity in viruses. This deeply sequenced virome also expanded the diversity of PolA sequences, increasing the number of PolA clusters by 44%. Contigs containing polA sequences were also used to examine relationships between phylogenetic clades of PolA and other genes within unknown viral populations. Phylogenic analysis revealed five distinct groups of phages distinguished by the amino acids at their 762 (Escherichia coli IAI39 numbering) positions and replication genes. DNA polymerase I sequences from Tyr762 and Phe762 groups were most often neighbored by ring-shaped superfamily IV helicases and ribonucleotide reductases (RNRs). The Leu762 groups had non-ring shaped helicases from superfamily II and were further distinguished by an additional helicase gene from superfamily I and the lack of any identifiable RNR genes. Moreover, we found that the inclusion of ribonucleotide reductase associated with PolA helped to further differentiate phage diversity, chiefly within lytic podovirus populations. Altogether, these data show that DNA Polymerase I is a useful marker for observing the diversity and composition of the virioplankton and may be a driving factor in the divergence of phage replication components.

Identification and sequence analysis of two novel cryptic plasmids isolated from the vaginal mucosa of South African women

The vaginal mucosa is dominated by Gram positive, rod shaped lactobacilli which serve as a natural barrier against infection. In both healthy- and bacterial vaginosis (BV)-infected women Lactobacillus crispatus and Lactobacillus jensenii have been found to be the predominant Lactobacillus species. Many studies have been conducted to assess factors influencing lactobacilli dominance in the vaginal microbiome. In the present study two plasmids, pLc4 and pLc17, isolated from vaginal Lactobacillus strains of both healthy and BV-infected women were characterized. The smaller plasmid, pLc4 (4224 bp), was detected in both L. crispatus and L. jensenii strains, while pLc17 was only detected in L. crispatus. Based on its nucleotide sequence pLc4 appears highly novel, with its replication protein having 44% identity to the replication initiation protein of pSMQ173b_03. Phylogenetic analysis with other Rolling Circle Replication plasmids confirmed that pLc4 shows a low degree of similarity to these plasmids. Plasmid pLc17 (16,663 bp) appears to carry both a RCR replicon and a theta replicon, which is rare in naturally occurring plasmids. pLc4 was maintained at a high copy number of 29, while pLc17 appears to be a medium copy number plasmid maintained at 11 copies per chromosome. While sequence analysis is a valuable tool to study cryptic plasmids, further function-based analysis will be required in order to fully elucidate the role of these plasmids within the vaginal milieu.

A Critical Review on the Use of Support Values in Tree Viewers and Bioinformatics Toolkits

Phylogenetic trees are routinely visualized to present and interpret the evolutionary relationships of species. Most empirical evolutionary data studies contain a visualization of the inferred tree with branch support values. Ambiguous semantics in tree file formats can lead to erroneous tree visualizations and therefore to incorrect interpretations of phylogenetic analyses. Here, we discuss problems that arise when displaying branch values on trees after rerooting. Branch values are typically stored as node labels in the widely-used Newick tree format. However, such values are attributes of branches. Storing them as node labels can therefore yield errors when rerooting trees. This depends on the mostly implicit semantics that tools deploy to interpret node labels. We reviewed ten tree viewers and ten bioinformatics toolkits that can display and reroot trees. We found that 14 out of 20 of these tools do not permit users to select the semantics of node labels. Thus, unaware users might obtain incorrect results when rooting trees. We illustrate such incorrect mappings for several test cases and real examples taken from the literature. This review has already led to improvements in eight tools. We suggest tools should provide options that explicitly force users to define the semantics of node labels.

Ribonucleotide Reductases from Bifidobacteria Contain Multiple Conserved Indels Distinguishing Them from All Other Organisms: In Silico Analysis of the Possible Role of a 43 aa Bifidobacteria-Specific Insert in the Class III RNR Homolog

Bifidobacteria comprises an important group/order of bacteria whose members have widespread usage in the food and health industry due to their health-promoting activity in the human gastrointestinal tract. However, little is known about the underlying molecular properties that are responsible for the probiotic effects of these bacteria. The enzyme ribonucleotide reductase (RNR) plays a key role in all organisms by reducing nucleoside di- or tri- phosphates into corresponding deoxyribose derivatives required for DNA synthesis, and RNR homologs belonging to classes I and III are present in either most or all Bifidobacteriales. Comparative analyses of these RNR homologs have identified several novel sequence features in the forms of conserved signature indels (CSIs) that are exclusively found in bifidobacterial RNRs. Specifically, in the large subunit of the aerobic class Ib RNR, three CSIs have been identified that are uniquely found in the Bifidobacteriales homologs. Similarly, the large subunit of the anaerobic class III RNR contains five CSIs that are also distinctive characteristics of bifidobacteria. Phylogenetic analyses indicate that these CSIs were introduced in a common ancestor of the Bifidobacteriales and retained by all descendants, likely due to their conferring advantageous functional roles. The identified CSIs in the bifidobacterial RNR homologs provide useful tools for further exploration of the novel functional aspects of these important enzymes that are exclusive to these bacteria. We also report here the results of homology modeling studies, which indicate that most of the bifidobacteria-specific CSIs are located within the surface loops of the RNRs, and of these, a large 43 amino acid insert in the class III RNR homolog forms an extension of the allosteric regulatory site known to be essential for protein function. Preliminary docking studies suggest that this large CSI may be playing a role in enhancing the stability of the RNR dimer complex. The possible significances of the identified CSIs, as well as the distribution of RNR homologs in the Bifidobacteriales, are discussed.

Two distinct pools of B12 analogs reveal community interdependencies in the ocean

Organisms within all domains of life require the cofactor cobalamin (vitamin B₁₂), which is produced only by a subset of bacteria and archaea. On the basis of genomic analyses, cobalamin biosynthesis in marine systems has been inferred in three main groups: select heterotrophic Proteobacteria, chemoautotrophic Thaumarchaeota, and photoautotrophic Cyanobacteria. Culture work demonstrates that many Cyanobacteria do not synthesize cobalamin but rather produce pseudocobalamin, challenging the connection between the occurrence of cobalamin biosynthesis genes and production of the compound in marine ecosystems. Here we show that cobalamin and pseudocobalamin coexist in the surface ocean, have distinct microbial sources, and support different enzymatic demands. Even in the presence of cobalamin, Cyanobacteria synthesize pseudocobalamin-likely reflecting their retention of an oxygen-independent pathway to produce pseudocobalamin, which is used as a cofactor in their specialized methionine synthase (MetH). This contrasts a model diatom, Thalassiosira pseudonana, which transported pseudocobalamin into the cell but was unable to use pseudocobalamin in its homolog of MetH. Our genomic and culture analyses showed that marine Thaumarchaeota and select heterotrophic bacteria produce cobalamin. This indicates that cobalamin in the surface ocean is a result of de novo synthesis by heterotrophic bacteria or via modification of closely related compounds like cyanobacterially produced pseudocobalamin. Deeper in the water column, our study implicates Thaumarchaeota as major producers of cobalamin based on genomic potential, cobalamin cell quotas, and abundance. Together, these findings establish the distinctive roles played by abundant prokaryotes in cobalamin-based microbial interdependencies that sustain community structure and function in the ocean.

Stable isotope labeling methods for DNA

NMR is a powerful method for studying proteins and nucleic acids in solution. The study of nucleic acids by NMR is far more challenging than for proteins, which is mainly due to the limited number of building blocks and unfavorable spectral properties. For NMR studies of DNA molecules, (site specific) isotope enrichment is required to facilitate specific NMR experiments and applications. Here, we provide a comprehensive review of isotope-labeling strategies for obtaining stable isotope labeled DNA as well as specifically stable isotope labeled building blocks required for enzymatic DNA synthesis.

A phylogenomic data-driven exploration of viral origins and evolution

The origin of viruses remains mysterious because of their diverse and patchy molecular and functional makeup. Although numerous hypotheses have attempted to explain viral origins, none is backed by substantive data. We take full advantage of the wealth of available protein structural and functional data to explore the evolution of the proteomic makeup of thousands of cells and viruses. Despite the extremely reduced nature of viral proteomes, we established an ancient origin of the "viral supergroup" and the existence of widespread episodes of horizontal transfer of genetic information. Viruses harboring different replicon types and infecting distantly related hosts shared many metabolic and informational protein structural domains of ancient origin that were also widespread in cellular proteomes. Phylogenomic analysis uncovered a universal tree of life and revealed that modern viruses reduced from multiple ancient cells that harbored segmented RNA genomes and coexisted with the ancestors of modern cells. The model for the origin and evolution of viruses and cells is backed by strong genomic and structural evidence and can be reconciled with existing models of viral evolution if one considers viruses to have originated from ancient cells and not from modern counterparts.

Biochemical Characterization of the Split Class II Ribonucleotide Reductase from Pseudomonas aeruginosa

The opportunistic pathogen Pseudomonas aeruginosa can grow under both aerobic and anaerobic conditions. Its flexibility with respect to oxygen load is reflected by the fact that its genome encodes all three existing classes of ribonucleotides reductase (RNR): the oxygen-dependent class I RNR, the oxygen-indifferent class II RNR, and the oxygen-sensitive class III RNR. The P. aeruginosa class II RNR is expressed as two separate polypeptides (NrdJa and NrdJb), a unique example of a split RNR enzyme in a free-living organism. A split class II RNR is also found in a few closely related γ-Proteobacteria. We have characterized the P. aeruginosa class II RNR and show that both subunits are required for formation of a biologically functional enzyme that can sustain vitamin B12-dependent growth. Binding of the B12 coenzyme as well as substrate and allosteric effectors resides in the NrdJa subunit, whereas the NrdJb subunit mediates efficient reductive dithiol exchange during catalysis. A combination of activity assays and activity-independent methods like surface plasmon resonance and gas phase electrophoretic macromolecule analysis suggests that the enzymatically active form of the enzyme is a (NrdJa-NrdJb)2 homodimer of heterodimers, and a combination of hydrogen-deuterium exchange experiments and molecular modeling suggests a plausible region in NrdJa that interacts with NrdJb. Our detailed characterization of the split NrdJ from P. aeruginosa provides insight into the biochemical function of a unique enzyme known to have central roles in biofilm formation and anaerobic growth.

Welcome to pandoraviruses at the 'Fourth TRUC' club

Nucleocytoplasmic large DNA viruses, or representatives of the proposed order Megavirales, belong to families of giant viruses that infect a broad range of eukaryotic hosts. Megaviruses have been previously described to comprise a fourth monophylogenetic TRUC (things resisting uncompleted classification) together with cellular domains in the universal tree of life. Recently described pandoraviruses have large (1.9–2.5 MB) and highly divergent genomes. In the present study, we updated the classification of pandoraviruses and other reported giant viruses. Phylogenetic trees were constructed based on six informational genes. Hierarchical clustering was performed based on a set of informational genes from Megavirales members and cellular organisms. Homologous sequences were selected from cellular organisms using TimeTree software, comprising comprehensive, and representative sets of members from Bacteria, Archaea, and Eukarya. Phylogenetic analyses based on three conserved core genes clustered pandoraviruses with phycodnaviruses, exhibiting their close relatedness. Additionally, hierarchical clustering analyses based on informational genes grouped pandoraviruses with Megavirales members as a super group distinct from cellular organisms. Thus, the analyses based on core conserved genes revealed that pandoraviruses are new genuine members of the ‘Fourth TRUC’ club, encompassing distinct life forms compared with cellular organisms.

Diversity in Overall Activity Regulation of Ribonucleotide Reductase

Ribonucleotide reductase (RNR) catalyzes the reduction of ribonucleotides to the corresponding deoxyribonucleotides, which are used as building blocks for DNA replication and repair. This process is tightly regulated via two allosteric sites, the specificity site (s-site) and the overall activity site (a-site). The a-site resides in an N-terminal ATP cone domain that binds dATP or ATP and functions as an on/off switch, whereas the composite s-site binds ATP, dATP, dTTP, or dGTP and determines which substrate to reduce. There are three classes of RNRs, and class I RNRs consist of different combinations of α and β subunits. In eukaryotic and Escherichia coli class I RNRs, dATP inhibits enzyme activity through the formation of inactive α6 and α4β4 complexes, respectively. Here we show that the Pseudomonas aeruginosa class I RNR has a duplicated ATP cone domain and represents a third mechanism of overall activity regulation. Each α polypeptide binds three dATP molecules, and the N-terminal ATP cone is critical for binding two of the dATPs because a truncated protein lacking this cone could only bind dATP to its s-site. ATP activates the enzyme solely by preventing dATP from binding. The dATP-induced inactive form is an α4 complex, which can interact with β2 to form a non-productive α4β2 complex. Other allosteric effectors induce a mixture of α2 and α4 forms, with the former being able to interact with β2 to form active α2β2 complexes. The unique features of the P. aeruginosa RNR are interesting both from evolutionary and drug discovery perspectives.

Function of the Pseudomonas aeruginosa NrdR Transcription Factor: Global Transcriptomic Analysis and Its Role on Ribonucleotide Reductase Gene Expression

Ribonucleotide reductases (RNRs) are a family of sophisticated enzymes responsible for the synthesis of the deoxyribonucleotides (dNTPs), the building blocks for DNA synthesis and repair. Although any living cell must contain one RNR activity to continue living, bacteria have the capacity to encode different RNR classes in the same genome, allowing them to adapt to different environments and growing conditions. Pseudomonas aeruginosa is well known for its adaptability and surprisingly encodes all three known RNR classes (Ia, II and III). There must be a complex transcriptional regulation network behind this RNR activity, dictating which RNR class will be expressed according to specific growing conditions. In this work, we aim to uncover the role of the transcriptional regulator NrdR in P. aeruginosa. We demonstrate that NrdR regulates all three RNR classes, being involved in differential control depending on whether the growth conditions are aerobic or anaerobic. Moreover, we also identify for the first time that NrdR is not only involved in controlling RNR expression but also regulates topoisomerase I (topA) transcription. Finally, to obtain the entire picture of NrdR regulon, we performed a global transcriptomic analysis comparing the transcription profile of wild-type and nrdR mutant strains. The results provide many new data about the regulatory network that controls P. aeruginosa RNR transcription, bringing us a step closer to the understanding of this complex system.

Methyl-hydroxylamine as an efficacious antibacterial agent that targets the ribonucleotide reductase enzyme

The emergence of multidrug-resistant bacteria has encouraged vigorous efforts to develop antimicrobial agents with new mechanisms of action. Ribonucleotide reductase (RNR) is a key enzyme in DNA replication that acts by converting ribonucleotides into the corresponding deoxyribonucleotides, which are the building blocks of DNA replication and repair. RNR has been extensively studied as an ideal target for DNA inhibition, and several drugs that are already available on the market are used for anticancer and antiviral activity. However, the high toxicity of these current drugs to eukaryotic cells does not permit their use as antibacterial agents. Here, we present a radical scavenger compound that inhibited bacterial RNR, and the compound's activity as an antibacterial agent together with its toxicity in eukaryotic cells were evaluated. First, the efficacy of N-methyl-hydroxylamine (M-HA) in inhibiting the growth of different Gram-positive and Gram-negative bacteria was demonstrated, and no effect on eukaryotic cells was observed. M-HA showed remarkable efficacy against Mycobacterium bovis BCG and Pseudomonas aeruginosa. Thus, given the M-HA activity against these two bacteria, our results showed that M-HA has intracellular antimycobacterial activity against BCG-infected macrophages, and it is efficacious in partially disassembling and inhibiting the further formation of P. aeruginosa biofilms. Furthermore, M-HA and ciprofloxacin showed a synergistic effect that caused a massive reduction in a P. aeruginosa biofilm. Overall, our results suggest the vast potential of M-HA as an antibacterial agent, which acts by specifically targeting a bacterial RNR enzyme.

The origin and evolution of ribonucleotide reduction

Ribonucleotide reduction is the only pathway for de novo synthesis of deoxyribonucleotides in extant organisms. This chemically demanding reaction, which proceeds via a carbon-centered free radical, is catalyzed by ribonucleotide reductase (RNR). The mechanism has been deemed unlikely to be catalyzed by a ribozyme, creating an enigma regarding how the building blocks for DNA were synthesized at the transition from RNA- to DNA-encoded genomes. While it is entirely possible that a different pathway was later replaced with the modern mechanism, here we explore the evolutionary and biochemical limits for an origin of the mechanism in the RNA + protein world and suggest a model for a prototypical ribonucleotide reductase (protoRNR). From the protoRNR evolved the ancestor to modern RNRs, the urRNR, which diversified into the modern three classes. Since the initial radical generation differs between the three modern classes, it is difficult to establish how it was generated in the urRNR. Here we suggest a model that is similar to the B12-dependent mechanism in modern class II RNRs.

Ribonucleotide reductase NrdR as a novel regulator for motility and chemotaxis during adherent-invasive Escherichia coli infection

A critical step in the life cycle of all organisms is the duplication of the genetic material during cell division. Ribonucleotide reductases (RNRs) are essential enzymes for this step because they control the de novo production of the deoxyribonucleotides required for DNA synthesis and repair. Enterobacteriaceae have three functional classes of RNRs (Ia, Ib, and III), which are transcribed from separate operons and encoded by the genes nrdAB, nrdHIEF, and nrdDG, respectively. Here, we investigated the role of RNRs in the virulence of adherent-invasive Escherichia coli (AIEC) isolated from Crohn's disease (CD) patients. Interestingly, the LF82 strain of AIEC harbors four different RNRs (two class Ia, one class Ib, and one class III). Although the E. coli RNR enzymes have been extensively characterized both biochemically and enzymatically, little is known about their roles during bacterial infection. We found that RNR expression was modified in AIEC LF82 bacteria during cell infection, suggesting that RNRs play an important role in AIEC virulence. Knockout of the nrdR and nrdD genes, which encode a transcriptional regulator of RNRs and class III anaerobic RNR, respectively, decreased AIEC LF82's ability to colonize the gut mucosa of transgenic mice that express human CEACAM6 (carcinoembryonic antigen-related cell adhesion molecule 6). Microarray experiments demonstrated that NrdR plays an indirect role in AIEC virulence by interfering with bacterial motility and chemotaxis. Thus, the development of drugs targeting RNR classes, in particular NrdR and NrdD, could be a promising new strategy to control gut colonization by AIEC bacteria in CD patients.

The case for an early biological origin of DNA

All life generates deoxyribonucleotides, the building blocks of DNA, via ribonucleotide reductases (RNRs). The complexity of this reaction suggests it did not evolve until well after the advent of templated protein synthesis, which in turn suggests DNA evolved later than both RNA and templated protein synthesis. However, deoxyribonucleotides may have first been synthesised via an alternative, chemically simpler route—the reversal of the deoxyriboaldolase (DERA) step in deoxyribonucleotide salvage. In light of recent work demonstrating that this reaction can drive synthesis of deoxyribonucleosides, we consider what pressures early adoption of this pathway would have placed on cell metabolism. This in turn provides a rationale for the replacement of DERA-dependent DNA production by RNR-dependent production.

Antibacterial gene transfer across the tree of life

Though horizontal gene transfer (HGT) is widespread, genes and taxa experience biased rates of transferability. Curiously, independent transmission of homologous DNA to archaea, bacteria, eukaryotes, and viruses is extremely rare and often defies ecological and functional explanations. Here, we demonstrate that a bacterial lysozyme family integrated independently in all domains of life across diverse environments, generating the only glycosyl hydrolase 25 muramidases in plants and archaea. During coculture of a hydrothermal vent archaeon with a bacterial competitor, muramidase transcription is upregulated. Moreover, recombinant lysozyme exhibits broad-spectrum antibacterial action in a dose-dependent manner. Similar to bacterial transfer of antibiotic resistance genes, transfer of a potent antibacterial gene across the universal tree seemingly bestows a niche-transcending adaptation that trumps the barriers against parallel HGT to all domains. The discoveries also comprise the first characterization of an antibacterial gene in archaea and support the pursuit of antibiotics in this underexplored group.

DOI:http://dx.doi.org/10.7554/eLife.04266.001

Living things inherit most of their genetic material from their parents, so genes tend to be passed on from one generation to the next—from ancestors to descendants. Sometimes, however, DNA is transferred from one organism to another by other means. These events, collectively called horizontal gene transfer, are fairly common in nature; genes have been passed between different species as well as between different groups of organisms. For example, genes that confer resistance to antibacterial drugs have transferred from one species of bacteria to another, and other genes have also ‘jumped’ from bacteria to plants or animals.

Now Metcalf et al. have studied a gene that first arose in bacteria and that encodes an enzyme called a lysozyme. This enzyme breaks down the outer casing of a bacterial cell: a step that is required for a bacterium to reproduce and divide in two. When Metcalf et al. searched for relatives of the lysozyme gene, they found copies in many other species of bacteria and revealed that this gene has been repeatedly transferred between different bacteria. Members of the lysozyme gene family have also ‘jumped out’ of bacteria and into other organisms at least four times. Metcalf et al. found related lysozyme genes in a plant, an insect, many species of fungi, and a single-celled microbe (called an archaeon) that lives at hot, deep-sea vents.

A gene family being spread this widely across the tree of life has not been seen before. Nevertheless, as DNA is a common biological language to all living things, it is likely that all the different species that have received a lysozyme gene might use it for similar purposes.

Metcalf et al. reveal that the lysozyme could be being used as an antibacterial molecule. The archaeon lysozyme can kill a broad range of bacteria; and when the gene was transferred into Escherichia coli bacteria, only the bacteria that mutated the lysozyme gene to render it useless were able to survive. Metcalf et al. also revealed that the archaeon microbe produces more of the enzyme if bacteria are present, which allows it to outcompete these bacteria.

These findings suggest that there may be a number of horizontally transferred genes that have antibacterial activity against a wide range of bacteria. Searching for these genes—particularly in the largely underexplored group of archaea—might reveal new sources for antibiotic drugs to treat bacterial infections.

DOI:http://dx.doi.org/10.7554/eLife.04266.002

Antibiotic genes spread far and wide

The genes responsible for antibiotics can spread between the three domains of life—Archaea, Bacteria and Eukaryotes.

Ribonucleotide reductases reveal novel viral diversity and predict biological and ecological features of unknown marine viruses

Virioplankton play a crucial role in aquatic ecosystems as top-down regulators of bacterial populations and agents of horizontal gene transfer and nutrient cycling. However, the biology and ecology of virioplankton populations in the environment remain poorly understood. Ribonucleotide reductases (RNRs) are ancient enzymes that reduce ribonucleotides to deoxyribonucleotides and thus prime DNA synthesis. Composed of three classes according to O2 reactivity, RNRs can be predictive of the physiological conditions surrounding DNA synthesis. RNRs are universal among cellular life, common within viral genomes and virioplankton shotgun metagenomes (viromes), and estimated to occur within >90% of the dsDNA virioplankton sampled in this study. RNRs occur across diverse viral groups, including all three morphological families of tailed phages, making these genes attractive for studies of viral diversity. Differing patterns in virioplankton diversity were clear from RNRs sampled across a broad oceanic transect. The most abundant RNRs belonged to novel lineages of podoviruses infecting α-proteobacteria, a bacterial class critical to oceanic carbon cycling. RNR class was predictive of phage morphology among cyanophages and RNR distribution frequencies among cyanophages were largely consistent with the predictions of the "kill the winner-cost of resistance" model. RNRs were also identified for the first time to our knowledge within ssDNA viromes. These data indicate that RNR polymorphism provides a means of connecting the biological and ecological features of virioplankton populations.

The class III ribonucleotide reductase from Neisseria bacilliformis can utilize thioredoxin as a reductant

The class III anaerobic ribonucleotide reductases (RNRs) studied to date couple the reduction of ribonucleotides to deoxynucleotides with the oxidation of formate to CO2. Here we report the cloning and heterologous expression of the Neisseria bacilliformis class III RNR and show that it can catalyze nucleotide reduction using the ubiquitous thioredoxin/thioredoxin reductase/NADPH system. We present a structural model based on a crystal structure of the homologous Thermotoga maritima class III RNR, showing its architecture and the position of conserved residues in the active site. Phylogenetic studies suggest that this form of class III RNR is present in bacteria and archaea that carry out diverse types of anaerobic metabolism.

Ribonucleotide reductases: essential enzymes for bacterial life

Ribonucleotide reductase (RNR) is a key enzyme that mediates the synthesis of deoxyribonucleotides, the DNA precursors, for DNA synthesis in every living cell. This enzyme converts ribonucleotides to deoxyribonucleotides, the building blocks for DNA replication, and repair. Clearly, RNR enzymes have contributed to the appearance of genetic material that exists today, being essential for the evolution of all organisms on Earth. The strict control of RNR activity and dNTP pool sizes is important, as pool imbalances increase mutation rates, replication anomalies, and genome instability. Thus, RNR activity should be finely regulated allosterically and at the transcriptional level. In this review we examine the distribution, the evolution, and the genetic regulation of bacterial RNRs. Moreover, this enzyme can be considered an ideal target for anti-proliferative compounds designed to inhibit cell replication in eukaryotic cells (cancer cells), parasites, viruses, and bacteria.

Cellular regulation of ribonucleotide reductase in eukaryotes

Synthesis of deoxynucleoside triphosphates (dNTPs) is essential for both DNA replication and repair and a key step in this process is catalyzed by ribonucleotide reductases (RNRs), which reduce ribonucleotides (rNDPs) to their deoxy forms. Tight regulation of RNR is crucial for maintaining the correct levels of all four dNTPs, which is important for minimizing the mutation rate and avoiding genome instability. Although allosteric control of RNR was the first discovered mechanism involved in regulation of the enzyme, other controls have emerged in recent years. These include regulation of expression of RNR genes, proteolysis of RNR subunits, control of the cellular localization of the small RNR subunit, and regulation of RNR activity by small protein inhibitors. This review will focus on these additional mechanisms of control responsible for providing a balanced supply of dNTPs.

Unveiling of the diversity of Prasinoviruses (Phycodnaviridae) in marine samples by using high-throughput sequencing analyses of PCR-amplified DNA polymerase and major capsid protein genes

Viruses strongly influence the ecology and evolution of their eukaryotic hosts in the marine environment, but little is known about their diversity and distribution. Prasinoviruses infect an abundant and widespread class of phytoplankton, the Mamiellophyceae, and thereby exert a specific and important role in microbial ecosystems. However, molecular tools to specifically identify this viral genus in environmental samples are still lacking. We developed two primer sets, designed for use with polymerase chain reactions and 454 pyrosequencing technologies, to target two conserved genes, encoding the DNA polymerase (PolB gene) and the major capsid protein (MCP gene). While only one copy of the PolB gene is present in Prasinovirus genomes, there are at least seven paralogs for MCP, the copy we named number 6 being shared with other eukaryotic alga-infecting viruses. Primer sets for PolB and MCP6 were thus designed and tested on 6 samples from the Tara Oceans project. The results suggest that the MCP6 amplicons show greater richness but that PolB gave a wider coverage of Prasinovirus diversity. As a consequence, we recommend use of the PolB primer set, which will certainly reveal exciting new insights about the diversity and distribution of prasinoviruses at the community scale.

The genome, proteome and phylogenetic analysis of Sinorhizobium meliloti phage ΦM12, the founder of a new group of T4-superfamily phages

Phage ΦM12 is an important transducing phage of the nitrogen-fixing rhizobial bacterium Sinorhizobium meliloti. Here we report the genome, phylogenetic analysis, and proteome of ΦM12, the first report of the genome and proteome of a rhizobium-infecting T4-superfamily phage. The structural genes of ΦM12 are most similar to T4-superfamily phages of cyanobacteria. ΦM12 is the first reported T4-superfamily phage to lack genes encoding class I ribonucleotide reductase (RNR) and exonuclease dexA, and to possess a class II coenzyme B12-dependent RNR. ΦM12's novel collection of genes establishes it as the founder of a new group of T4-superfamily phages, fusing features of cyanophages and phages of enteric bacteria.

A novel clade of unique eukaryotic ribonucleotide reductase R2 subunits is exclusive to apicomplexan parasites

Apicomplexa are protist parasites of tremendous medical and economic importance, causing millions of deaths and billions of dollars in losses each year. Apicomplexan-related diseases may be controlled via inhibition of essential enzymes. Ribonucleotide reductase (RNR) provides the only de novo means of synthesizing deoxyribonucleotides, essential precursors for DNA replication and repair. RNR has long been the target of antibacterial and antiviral therapeutics. However, targeting this ubiquitous protein in eukaryotic pathogens may be problematic unless these proteins differ significantly from that of their respective host. The typical eukaryotic RNR enzymes belong to class Ia, and the holoenzyme consists minimally of two R1 and two R2 subunits (α₂β₂). We generated a comparative, annotated, structure-based, multiple-sequence alignment of R2 subunits, identified a clade of R2 subunits unique to Apicomplexa, and determined its phylogenetic position. Our analyses revealed that the apicomplexan-specific sequences share characteristics with both class I R2 and R2lox proteins. The putative radical-harboring residue, essential for the reduction reaction by class Ia R2-containing holoenzymes, was not conserved within this group. Phylogenetic analyses suggest that class Ia subunits are not monophyletic and consistently placed the apicomplexan-specific clade sister to the remaining class Ia eukaryote R2 subunits. Our research suggests that the novel apicomplexan R2 subunit may be a promising candidate for chemotherapeutic-induced inhibition as it differs greatly from known eukaryotic host RNRs and may be specifically targeted.

A bioinformatic analysis of ribonucleotide reductase genes in phage genomes and metagenomes

Background

Ribonucleotide reductase (RNR), the enzyme responsible for the formation of deoxyribonucleotides from ribonucleotides, is found in all domains of life and many viral genomes. RNRs are also amongst the most abundant genes identified in environmental metagenomes. This study focused on understanding the distribution, diversity, and evolution of RNRs in phages (viruses that infect bacteria). Hidden Markov Model profiles were used to analyze the proteins encoded by 685 completely sequenced double-stranded DNA phages and 22 environmental viral metagenomes to identify RNR homologs in cultured phages and uncultured viral communities, respectively.

Results

RNRs were identified in 128 phage genomes, nearly tripling the number of phages known to encode RNRs. Class I RNR was the most common RNR class observed in phages (70%), followed by class II (29%) and class III (28%). Twenty-eight percent of the phages contained genes belonging to multiple RNR classes. RNR class distribution varied according to phage type, isolation environment, and the host’s ability to utilize oxygen. The majority of the phages containing RNRs are Myoviridae (65%), followed by Siphoviridae (30%) and Podoviridae (3%). The phylogeny and genomic organization of phage and host RNRs reveal several distinct evolutionary scenarios involving horizontal gene transfer, co-evolution, and differential selection pressure. Several putative split RNR genes interrupted by self-splicing introns or inteins were identified, providing further evidence for the role of frequent genetic exchange. Finally, viral metagenomic data indicate that RNRs are prevalent and highly dynamic in uncultured viral communities, necessitating future research to determine the environmental conditions under which RNRs provide a selective advantage.

Conclusions

This comprehensive study describes the distribution, diversity, and evolution of RNRs in phage genomes and environmental viral metagenomes. The distinct distributions of specific RNR classes amongst phages, combined with the various evolutionary scenarios predicted from RNR phylogenies suggest multiple inheritance sources and different selective forces for RNRs in phages. This study significantly improves our understanding of phage RNRs, providing insight into the diversity and evolution of this important auxiliary metabolic gene as well as the evolution of phages in response to their bacterial hosts and environments.

A rare combination of ribonucleotide reductases in the social amoeba Dictyostelium discoideum

Ribonucleotide reductases (RNRs) catalyze the only pathway for de novo synthesis of deoxyribonucleotides needed for DNA replication and repair. The vast majority of eukaryotes encodes only a class I RNR, but interestingly some eukaryotes, including the social amoeba Dictyostelium discoideum, encode both a class I and a class II RNR. The amino acid sequence of the D. discoideum class I RNR is similar to other eukaryotic RNRs, whereas that of its class II RNR is most similar to the monomeric class II RNRs found in Lactobacillus spp. and a few other bacteria. Here we report the first study of RNRs in a eukaryotic organism that encodes class I and class II RNRs. Both classes of RNR genes were expressed in D. discoideum cells, although the class I transcripts were more abundant and strongly enriched during mid-development compared with the class II transcript. The quaternary structure, allosteric regulation, and properties of the diiron-oxo/radical cofactor of D. discoideum class I RNR are similar to those of the mammalian RNRs. Inhibition of D. discoideum class I RNR by hydroxyurea resulted in a 90% reduction in spore formation and decreased the germination viability of the surviving spores by 75%. Class II RNR could not compensate for class I inhibition during development, and an excess of vitamin B₁₂ coenzyme, which is essential for class II activity, did not improve spore formation. We suggest that class I is the principal RNR during D. discoideum development and growth and is important for spore formation, possibly by providing dNTPs for mitochondrial replication.

Use of structural phylogenetic networks for classification of the ferritin-like superfamily

In the postgenomic era, bioinformatic analysis of sequence similarity is an immensely powerful tool to gain insight into evolution and protein function. Over long evolutionary distances, however, sequence-based methods fail as the similarities become too low for phylogenetic analysis. Macromolecular structure generally appears better conserved than sequence, but clear models for how structure evolves over time are lacking. The exponential growth of three-dimensional structural information may allow novel structure-based methods to drastically extend the evolutionary time scales amenable to phylogenetics and functional classification of proteins. To this end, we analyzed 80 structures from the functionally diverse ferritin-like superfamily. Using evolutionary networks, we demonstrate that structural comparisons can delineate and discover groups of proteins beyond the "twilight zone" where sequence similarity does not allow evolutionary analysis, suggesting that considerable and useful evolutionary signal is preserved in three-dimensional structures.

DNA building blocks: keeping control of manufacture

Ribonucleotide reductase (RNR) is the only source for de novo production of the four deoxyribonucleoside triphosphate (dNTP) building blocks needed for DNA synthesis and repair. It is crucial that these dNTP pools are carefully balanced, since mutation rates increase when dNTP levels are either unbalanced or elevated. RNR is the major player in this homeostasis, and with its four different substrates, four different allosteric effectors and two different effector binding sites, it has one of the most sophisticated allosteric regulations known today. In the past few years, the structures of RNRs from several bacteria, yeast and man have been determined in the presence of allosteric effectors and substrates, revealing new information about the mechanisms behind the allosteric regulation. A common theme for all studied RNRs is a flexible loop that mediates modulatory effects from the allosteric specificity site (s-site) to the catalytic site for discrimination between the four substrates. Much less is known about the allosteric activity site (a-site), which functions as an on-off switch for the enzyme's overall activity by binding ATP (activator) or dATP (inhibitor). The two nucleotides induce formation of different enzyme oligomers, and a recent structure of a dATP-inhibited α(6)β(2) complex from yeast suggested how its subunits interacted non-productively. Interestingly, the oligomers formed and the details of their allosteric regulation differ between eukaryotes and Escherichia coli. Nevertheless, these differences serve a common purpose in an essential enzyme whose allosteric regulation might date back to the era when the molecular mechanisms behind the central dogma evolved.