The identification of microbial genes essential for survival as those with lethal knockout phenotype (LKP) is a common strategy for functional interrogation of genomes. However, interpretation of the LKP is complicated because a substantial fraction of the genes with this phenotype remains poorly functionally characterized. Furthermore, many genes can exhibit LKP not because their products perform essential cellular functions but because their knockout activates the toxicity of other genes (conditionally essential genes). We analyzed the sets of LKP genes for two archaea, Methanococcus maripaludis and Sulfolobus islandicus, using a variety of computational approaches aiming to differentiate between essential and conditionally essential genes and to predict at least a general function for as many of the proteins encoded by these genes as possible. This analysis allowed us to predict the functions of several LKP genes including previously uncharacterized subunit of the GINS protein complex with an essential function in genome replication and of the KEOPS complex that is responsible for an essential tRNA modification as well as GRP protease implicated in protein quality control. Additionally, several novel antitoxins (conditionally essential genes) were predicted, and this prediction was experimentally validated by showing that the deletion of these genes together with the adjacent genes apparently encoding the cognate toxins caused no growth defect. We applied principal component analysis based on sequence and comparative genomic features showing that this approach can separate essential genes from conditionally essential ones and used it to predict essential genes in other archaeal genomes.IMPORTANCEOnly a relatively small fraction of the genes in any bacterium or archaeon is essential for survival as demonstrated by the lethal effect of their disruption. The identification of essential genes and their functions is crucial for understanding fundamental cell biology. However, many of the genes with a lethal knockout phenotype remain poorly functionally characterized, and furthermore, many genes can exhibit this phenotype not because their products perform essential cellular functions but because their knockout activates the toxicity of other genes. We applied state-of-the-art computational methods to predict the functions of a number of uncharacterized genes with the lethal knockout phenotype in two archaeal species and developed a computational approach to predict genes involved in essential functions. These findings advance the current understanding of key functionalities of archaeal cells.

The survival of viruses depends on their ability to resist host defenses and, of all animal virus families, the poxviruses have the most antidefense genes. Orthopoxviruses (ORPV), a genus within the subfamily Chordopoxvirinae, infect diverse mammals and include one of the most devastating human pathogens, the now eradicated smallpox virus. ORPV encode ∼200 genes, of which roughly half are directly involved in virus genome replication and expression as well as virion morphogenesis. The remaining ∼100 "accessory" genes are responsible for virus-host interactions, particularly counter-defense of innate immunity. Complete sequences are currently available for several hundred ORPV genomes isolated from a variety of mammalian hosts, providing a rich resource for comparative genomics and reconstruction of ORPV evolution. To identify the provenance and evolutionary trends of the ORPV accessory genes, we constructed clusters including the orthologs of these genes from all chordopoxviruses. Most of the accessory genes were captured in three major waves early in chordopoxvirus evolution, prior to the divergence of ORPV and the sister genus Centapoxvirus from their common ancestor. The capture of these genes from the host was followed by extensive gene duplication, yielding several paralogous gene families. In addition, nine genes were gained during the evolution of ORPV themselves. In contrast, nearly every accessory gene was lost, some on multiple, independent occasions in numerous lineages of ORPV, so that no ORPV retains them all. A variety of functional interactions could be inferred from examination of pairs of ORPV accessory genes that were either often or rarely lost concurrently. IMPORTANCE Orthopoxviruses (ORPV) include smallpox (variola) virus, one of the most devastating human pathogens, and vaccinia virus, comprising the vaccine used for smallpox eradication. Among roughly 200 ORPV genes, about half are essential for genome replication and expression as well as virion morphogenesis, whereas the remaining half consists of accessory genes counteracting the host immune response. We reannotated the accessory genes of ORPV, predicting the functions of uncharacterized genes, and reconstructed the history of their gain and loss during the evolution of ORPV. Most of the accessory genes were acquired in three major waves antedating the origin of ORPV from chordopoxviruses. The evolution of ORPV themselves was dominated by gene loss, with numerous genes lost at the base of each major group of ORPV. Examination of pairs of ORPV accessory genes that were either often or rarely lost concurrently during ORPV evolution allows prediction of different types of functional interactions.

Gene transfer agents (GTAs) are virus-like elements integrated into bacterial genomes, particularly, those of Alphaproteobacteria The GTAs can be induced under conditions of nutritional stress, incorporate random fragments of bacterial DNA into miniphage particles, lyse the host cells, and infect neighboring bacteria, thus enhancing horizontal gene transfer. We show that GTA genes evolve under conditions of pronounced positive selection for the reduction of the energy cost of protein production as shown by comparison of the amino acid compositions with those of both homologous viral genes and host genes. The energy saving in GTA genes is comparable to or even more pronounced than that in the genes encoding the most abundant, essential bacterial proteins. In cases in which viruses acquire genes from GTAs, the bias in amino acid composition disappears in the course of evolution, showing that reduction of the energy cost of protein production is an important factor of evolution of GTAs but not bacterial viruses. These findings strongly suggest that GTAs represent bacterial adaptations rather than selfish, virus-like elements. Because GTA production kills the host cell and does not propagate the GTA genome, it appears likely that the GTAs are retained in the course of evolution via kin or group selection. Therefore, we hypothesize that GTAs facilitate the survival of bacterial populations under energy-limiting conditions through the spread of metabolic and transport capabilities via horizontal gene transfer and increases in nutrient availability resulting from the altruistic suicide of GTA-producing cells.IMPORTANCE Kin selection and group selection remain controversial topics in evolutionary biology. We argue that these types of selection are likely to operate in bacterial populations by showing that bacterial gene transfer agents (GTAs), but not related viruses, evolve under conditions of positive selection for the reduction of the energy cost of GTA particle production. We hypothesize that GTAs are dedicated devices mediating the survival of bacteria under conditions of nutrient limitation. The benefits conferred by GTAs under nutritional stress conditions appear to include horizontal dissemination of genes that could provide bacteria with enhanced capabilities for nutrient utilization and increases of nutrient availability occurring through the lysis of GTA-producing bacteria.

The emergence of the endomembrane system is a key step in the evolution of cellular complexity during eukaryogenesis. The endosomal sorting complex required for transport (ESCRT) machinery is essential and required for the endomembrane system functions in eukaryotic cells. Recently, genes encoding eukaryote-like ESCRT protein components have been identified in the genomes of Asgard archaea, a newly proposed archaeal superphylum that is thought to include the closest extant prokaryotic relatives of eukaryotes. However, structural and functional features of Asgard ESCRT remain uncharacterized. Here, we show that Vps4, Vps2/24/46, and Vps20/32/60, the core functional components of the Asgard ESCRT, coevolved eukaryote-like structural and functional features. Phylogenetic analysis shows that Asgard Vps4, Vps2/24/46, and Vps20/32/60 are closely related to their eukaryotic counterparts. Molecular dynamics simulation and biochemical assays indicate that Asgard Vps4 contains a eukaryote-like microtubule-interacting and transport (MIT) domain that binds the distinct type 1 MIT-interacting motif and type 2 MIT-interacting motif in Vps2/24/46 and Vps20/32/60, respectively. The Asgard Vps4 partly, but much more efficiently than homologs from other archaea, complements the vps4 null mutant of Saccharomyces cerevisiae, further supporting the functional similarity between the membrane remodeling machineries of Asgard archaea and eukaryotes. Thus, this work provides evidence that the ESCRT complexes from Asgard archaea and eukaryotes are evolutionarily related and functionally similar. Thus, despite the apparent absence of endomembranes in Asgard archaea, the eukaryotic ESCRT seems to have been directly inherited from an Asgard ancestor, to become a key component of the emerging endomembrane system.IMPORTANCE The discovery of Asgard archaea has changed the existing ideas on the origins of eukaryotes. Researchers propose that eukaryotic cells evolved from Asgard archaea. This hypothesis partly stems from the presence of multiple eukaryotic signature proteins in Asgard archaea, including homologs of ESCRT proteins that are essential components of the endomembrane system in eukaryotes. However, structural and functional features of Asgard ESCRT remain unknown. Our study provides evidence that Asgard ESCRT is functionally comparable to the eukaryotic counterparts, suggesting that despite the apparent absence of endomembranes in archaea, eukaryotic ESCRT was inherited from an Asgard archaeal ancestor, alongside the emergence of endomembrane system during eukaryogenesis.

The catalytic subunit of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA-dependent RNA polymerase (RdRp) Nsp12 has a unique nidovirus RdRp-associated nucleotidyltransferase (NiRAN) domain that transfers nucleoside monophosphates to the Nsp9 protein and the nascent RNA. The NiRAN and RdRp modules form a dynamic interface distant from their catalytic sites, and both activities are essential for viral replication. We report that codon-optimized (for the pause-free translation in bacterial cells) Nsp12 exists in an inactive state in which NiRAN-RdRp interactions are broken, whereas translation by slow ribosomes and incubation with accessory Nsp7/8 subunits or nucleoside triphosphates (NTPs) partially rescue RdRp activity. Our data show that adenosine and remdesivir triphosphates promote the synthesis of A-less RNAs, as does ppGpp, while amino acid substitutions at the NiRAN-RdRp interface augment activation, suggesting that ligand binding to the NiRAN catalytic site modulates RdRp activity. The existence of allosterically linked nucleotidyl transferase sites that utilize the same substrates has important implications for understanding the mechanism of SARS-CoV-2 replication and the design of its inhibitors. IMPORTANCEIn vitro interrogations of the central replicative complex of SARS-CoV-2, RNA-dependent RNA polymerase (RdRp), by structural, biochemical, and biophysical methods yielded an unprecedented windfall of information that, in turn, instructs drug development and administration, genomic surveillance, and other aspects of the evolving pandemic response. They also illuminated the vast disparity in the methods used to produce RdRp for experimental work and the hidden impact that this has on enzyme activity and research outcomes. In this report, we elucidate the positive and negative effects of codon optimization on the activity and folding of the recombinant RdRp and detail the design of a highly sensitive in vitro assay of RdRp-dependent RNA synthesis. Using this assay, we demonstrate that RdRp is allosterically activated by nontemplating phosphorylated nucleotides, including naturally occurring alarmone ppGpp and synthetic remdesivir triphosphate.

Dephospho-coenzyme A (dephospho-CoA) kinase (DPCK) catalyzes the ATP-dependent phosphorylation of dephospho-CoA, the final step in coenzyme A (CoA) biosynthesis. DPCK has been identified and characterized in bacteria and eukaryotes but not in archaea. The hyperthermophilic archaeon Thermococcus kodakarensis encodes two homologs of bacterial DPCK and the DPCK domain of eukaryotic CoA synthase, TK1334 and TK2192. We purified the recombinant TK1334 and TK2192 proteins and found that they lacked DPCK activity. Bioinformatic analyses showed that, in several archaea, the uncharacterized gene from arCOG04076 protein is fused with the gene for phosphopantetheine adenylyltransferase (PPAT), which catalyzes the reaction upstream of the DPCK reaction in CoA biosynthesis. This observation suggested that members of arCOG04076, both fused to PPAT and standalone, could be the missing archaeal DPCKs. We purified the recombinant TK1697 protein, a standalone member of arCOG04076 from T. kodakarensis, and demonstrated its GTP-dependent DPCK activity. Disruption of the TK1697 resulted in CoA auxotrophy, indicating that TK1697 encodes a DPCK that contributes to CoA biosynthesis in T. kodakarensis TK1697 homologs are widely distributed in archaea, suggesting that the arCOG04076 protein represents a novel family of DPCK that is not homologous to bacterial and eukaryotic DPCKs but is distantly related to bacterial and eukaryotic thiamine pyrophosphokinases. We also constructed and characterized gene disruption strains of TK0517 and TK2128, homologs of bifunctional phosphopantothenoylcysteine synthetase-phosphopantothenoylcysteine decarboxylase and PPAT, respectively. Both strains displayed CoA auxotrophy, indicating their contribution to CoA biosynthesis. Taken together with previous studies, the results experimentally validate the entire CoA biosynthesis pathway in T. kodakarensisIMPORTANCE CoA is utilized in a wide range of metabolic pathways, and its biosynthesis is essential for all life. Pathways for CoA biosynthesis in bacteria and eukaryotes have been established. In archaea, however, the enzyme that catalyzes the final step in CoA biosynthesis, dephospho-CoA kinase (DPCK), had not been identified. In the present study, bioinformatic analyses identified a candidate for the DPCK in archaea, which was biochemically and genetically confirmed in the hyperthermophilic archaeon Thermococcus kodakarensis Genetic analyses on genes presumed to encode bifunctional phosphopantothenoylcysteine synthetase-phosphopantothenoylcysteine decarboxylase and phosphopantetheine adenylyltransferase confirmed their involvement in CoA biosynthesis. Taken together with previous studies, the results reveal the entire pathway for CoA biosynthesis in a single archaeon and provide insight into the different mechanisms of CoA biosynthesis and their distribution in nature.