Search

The Clusters of Orthologous Genes (COG) database, originally created in 1997, has been updated to reflect the constantly growing collection of completely sequenced prokaryotic genomes. This update increased the genome coverage from 1309 to 2296 species, including 2103 bacteria and 193 archaea, in most cases, with a single representative genome per genus. This set covers all genera of bacteria and archaea that included organisms with 'complete genomes' as per NCBI databases in November 2023. The number of COGs has been expanded from 4877 to 4981, primarily by including protein families involved in bacterial protein secretion. Accordingly, COG pathways and functional groups now include secretion systems of types II through X, as well as Flp/Tad and type IV pili. These groupings allow straightforward identification and examination of the prokaryotic lineages that encompass-or lack-a particular secretion system. Other developments include improved annotations for the rRNA and tRNA modification proteins, multi-domain signal transduction proteins, and some previously uncharacterized protein families. The new version of COGs is available at https://www.ncbi.nlm.nih.gov/research/COG, as well as on the NCBI FTP site https://ftp.ncbi.nlm.nih.gov/pub/COG/, which also provides archived data from previous COG releases., (Published by Oxford University Press on behalf of Nucleic Acids Research 2024.)

Bacterial and archaeal genomes encompass numerous operons that typically consist of two to five genes. On larger scales, however, gene order is poorly conserved through the evolution of prokaryotes. Nevertheless, non-random localization of different classes of genes on prokaryotic chromosomes could reflect important functional and evolutionary constraints. We explored the patterns of genomic localization of evolutionarily conserved (ancient) and variable (young) genes across the diversity of bacteria and archaea. Nearly all bacterial and archaeal chromosomes were found to encompass large segments of 100-300 kb that were significantly enriched in either ancient or young genes. Similar clustering of genes with lethal knockout phenotype (essential genes) was observed as well. Mathematical modeling of genome evolution suggests that this long-range gene clustering in prokaryotic chromosomes reflects perpetual genome rearrangement driven by a combination of selective and neutral processes rather than evolutionary conservation., (Published by Oxford University Press on behalf of Nucleic Acids Research 2024.)

The most abundant clustered regularly interspaced short palindromic repeats (CRISPR) type I systems employ a multisubunit RNA-protein effector complex (Cascade), with varying protein composition and activity. The Escherichia coli Cascade complex consists of 11 protein subunits and functions as an effector through CRISPR RNA (crRNA) binding, protospacer adjacent motif (PAM)-specific double-stranded DNA targeting, R-loop formation, and Cas3 helicase-nuclease recruitment for target DNA cleavage. Here, we present a biochemical reconstruction of the E. coli Cascade from purified Cas proteins and analyze its activities including crRNA binding, dsDNA targeting, R-loop formation, and Cas3 recruitment. Affinity purification of 6His-tagged Cas7 coexpressed with untagged Cas5 revealed the physical association of these proteins, thus producing the Cas5-Cas7 subcomplex that was able to bind specifically to type I-E crRNA with an efficiency comparable to that of the complete Cascade. The crRNA-loaded Cas5-7 was found to bind specifically to the target dsDNA in a PAM-independent manner, albeit with a lower affinity than the complete Cascade, with both spacer sequence complementarity and repeat handles contributing to the DNA targeting specificity. The crRNA-loaded Cas5-7 targeted the complementary dsDNA with detectable formation of R-loops, which was stimulated by the addition of Cas8 and/or Cas11 acting synergistically. Cascade activity reconstitution using purified Cas5-7 and other Cas proteins showed that Cas8 was essential for specific PAM recognition, whereas the addition of Cas11 was required for Cas3 recruitment and target DNA nicking. Thus, although the core Cas5-7 subcomplex is sufficient for specific crRNA binding and basal DNA targeting, both Cas8 and Cas11 make unique contributions to efficient target recognition and cleavage., (© 2024 Federation of European Biochemical Societies.)

SUMMARYPolintons are 15-20 kb-long self-synthesizing transposons that are widespread in eukaryotic, and in particular protist, genomes. Apart from a transposase and a protein-primed DNA polymerase, polintons encode homologs of major and minor jelly-roll capsid proteins, DNA-packaging ATPases, and proteases involved in capsid maturation of diverse eukaryotic viruses of kingdom Bamfordvirae . Given the conservation of these structural and morphogenetic proteins among polintons, these elements are predicted to alternate between transposon and viral lifestyles and, although virions have thus far not been detected, are classified as viruses (class Polintoviricetes ) in the phylum Preplasmiviricota . Related to polintoviricetes are vertebrate adenovirids; unclassified polinton-like viruses (PLVs) identified in various environments or integrated into diverse protist genomes; virophages ( Maveriviricetes ), which are part of tripartite hyperparasitic systems including protist hosts and giant viruses; and capsid-less derivatives, such as cytoplasmic linear DNA plasmids of fungi and transpovirons. Phylogenomic analysis indicates that the polinton-like supergroup of viruses bridges bacterial tectivirids (preplasmiviricot class Tectiliviricetes ) to the phylum Nucleocytoviricota that includes large and giant eukaryotic DNA viruses. Comparative structural analysis of proteins encoded by polinton-like viruses led to the discovery of previously undetected functional domains, such as terminal proteins and distinct proteases implicated in DNA polymerase processing, and clarified the evolutionary relationships within Polintoviricetes . Here, we leverage these insights into the evolution of the polinton-like supergroup to develop an amended megataxonomy that groups Polintoviricetes , PLVs (new class ' Aquintoviricetes '), and virophages (renamed class ' Virophaviricetes ') together with Adenoviridae (new class ' Pharingeaviricetes ') in a preplasmiviricot subphylum ' Polisuviricotina ' sister to a subphylum including Tectiliviricetes (' Prepoliviricotina ')., Competing Interests: The authors declare no conflict of interest.

Background: Viruses with double-stranded (ds) DNA genomes in the realm Duplodnaviria share a conserved structural gene module but show a broad range of variation in their repertoires of DNA replication proteins. Some of the duplodnaviruses encode (nearly) complete replication systems whereas others lack (almost) all genes required for replication, relying on the host replication machinery. DNA polymerases (DNAPs) comprise the centerpiece of the DNA replication apparatus. The replicative DNAPs are classified into 4 unrelated or distantly related families (A-D), with the protein structures and sequences within each family being, generally, highly conserved. More than half of the duplodnaviruses encode a DNAP of family A, B or C. We showed previously that multiple pairs of closely related viruses in the order Crassvirales encode DNAPs of different families., Methods: Groups of phages in which DNAP swapping likely occurred were identified as subtrees of a defined depth in a comprehensive evolutionary tree of tailed bacteriophages that included phages with DNAPs of different families. The DNAP swaps were validated by constrained tree analysis that was performed on phylogenetic tree of large terminase subunits, and the phage genomes encoding swapped DNAPs were aligned using Mauve. The structures of the discovered unusual DNAPs were predicted using AlphaFold2., Results: We identified four additional groups of tailed phages in the class Caudoviricetes in which the DNAPs apparently were swapped on multiple occasions, with replacements occurring both between families A and B, or A and C, or between distinct subfamilies within the same family. The DNAP swapping always occurs "in situ", without changes in the organization of the surrounding genes. In several cases, the DNAP gene is the only region of substantial divergence between closely related phage genomes, whereas in others, the swap apparently involved neighboring genes encoding other proteins involved in phage genome replication. In addition, we identified two previously undetected, highly divergent groups of family A DNAPs that are encoded in some phage genomes along with the main DNAP implicated in genome replication., Conclusions: Replacement of the DNAP gene by one encoding a DNAP of a different family occurred on many independent occasions during the evolution of different families of tailed phages, in some cases, resulting in very closely related phages encoding unrelated DNAPs. DNAP swapping was likely driven by selection for avoidance of host antiphage mechanisms targeting the phage DNAP that remain to be identified, and/or by selection against replicon incompatibility., (© 2024. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

The phylum Nucleocytoviricota consists of large and giant viruses that range in genome size from about 100 kilobases (kb) to more than 2.5 megabases. Here, using metagenome mining followed by extensive phylogenomic analysis and protein structure comparison, we delineate a distinct group of viruses with double-stranded (ds) DNA genomes in the range of 35-45 kb that appear to be related to the Nucleocytoviricota . In phylogenetic trees of the conserved double jelly-roll major capsid proteins (MCPs) and DNA packaging ATPases, these viruses do not show affinity to any particular branch of the Nucleocytoviricota and accordingly would comprise a class which we propose to name " Mriyaviricetes " (after Ukrainian "mriya," dream). Structural comparison of the MCP suggests that, among the extant virus lineages, mriyaviruses are the closest one to the ancestor of the Nucleocytoviricota . In the phylogenetic trees, mriyaviruses split into two well-separated branches, the family Yaraviridae and proposed new family " Gamadviridae. " The previously characterized members of these families, yaravirus and Pleurochrysis sp. endemic viruses, infect amoeba and haptophytes, respectively. The genomes of the rest of the mriyaviruses were assembled from metagenomes from diverse environments, suggesting that mriyaviruses infect various unicellular eukaryotes. Mriyaviruses lack DNA polymerase, which is encoded by all other members of the Nucleocytoviricota , and RNA polymerase subunits encoded by all cytoplasmic viruses among the Nucleocytoviricota , suggesting that they replicate in the host cell nuclei. All mriyaviruses encode a HUH superfamily endonuclease that is likely to be essential for the initiation of virus DNA replication via the rolling circle mechanism., Importance: The origin of giant viruses of eukaryotes that belong to the phylum Nucleocytoviricota is not thoroughly understood and remains a matter of major interest and debate. Here, we combine metagenome database searches with extensive protein sequence and structure analysis to describe a distinct group of viruses with comparatively small genomes of 35-45 kilobases that appear to comprise a distinct class within the phylum Nucleocytoviricota that we provisionally named " Mriyaviricetes ." Mriyaviruses appear to be the closest identified relatives of the ancestors of the Nucleocytoviricota . Analysis of proteins encoded in mriyavirus genomes suggests that they replicate their genome via the rolling circle mechanism that is unusual among viruses with double-stranded DNA genomes and so far not described for members of Nucleocytoviricota ., Competing Interests: The authors declare no conflict of interest.

Background: Microbiomes are generally characterized by high diversity of coexisting microbial species and strains, and microbiome composition typically remains stable across a broad range of conditions. However, under fixed conditions, microbial ecology conforms with the exclusion principle under which two populations competing for the same resource within the same niche cannot coexist because the less fit population inevitably goes extinct. Therefore, the long-term persistence of microbiome diversity calls for an explanation., Results: To explore the conditions for stabilization of microbial diversity, we developed a simple mathematical model consisting of two competing populations that could exchange a single gene allele via horizontal gene transfer (HGT). We found that, although in a fixed environment, with unbiased HGT, the system obeyed the exclusion principle, in an oscillating environment, within large regions of the phase space bounded by the rates of reproduction and HGT, the two populations coexist. Moreover, depending on the parameter combination, all three major types of symbiosis were obtained, namely, pure competition, host-parasite relationship, and mutualism. In each of these regimes, certain parameter combinations provided for synergy, that is, a greater total abundance of both populations compared to the abundance of the winning population in the fixed environment., Conclusions: The results of this modeling study show that basic phenomena that are universal in microbial communities, namely, environmental variation and HGT, provide for stabilization and persistence of microbial diversity, and emergence of ecological complexity., (© 2024. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

The phylum Preplasmiviricota (kingdom Bamfordvirae , realm Varidnaviria ) is a broad assemblage of diverse viruses with comparatively short double-stranded DNA genomes (<50 kbp) that produce icosahedral capsids built from double jelly-roll major capsid proteins. Preplasmiviricots infect hosts from all cellular domains, testifying to their ancient origin, and, in particular, are associated with six of the seven supergroups of eukaryotes. Preplasmiviricots comprise four major groups of viruses, namely, polintons, polinton-like viruses (PLVs), virophages, and adenovirids. We used protein structure modeling and analysis to show that protein-primed DNA polymerases (pPolBs) of polintons, virophages, and cytoplasmic linear plasmids encompass an N-terminal domain homologous to the terminal proteins (TPs) of prokaryotic PRD1-like tectivirids and eukaryotic adenovirids that are involved in protein-primed replication initiation, followed by a viral ovarian tumor-like cysteine deubiquitinylase (vOTU) domain. The vOTU domain is likely responsible for the cleavage of the TP from the large pPolB polypeptide and is inactivated in adenovirids, in which TP is a separate protein. Many PLVs and transpovirons encode a distinct derivative of polinton-like pPolB that retains the TP, vOTU, and pPolB polymerization palm domains but lacks the exonuclease domain and instead contains a superfamily 1 helicase domain. Analysis of the presence/absence and inactivation of the vOTU domains and replacement of pPolB with other DNA polymerases in eukaryotic preplasmiviricots enabled us to outline a complete scenario for their origin and evolution., Competing Interests: Competing interests statement:The authors declare no competing interest.

Over the course of multiple divisions, cells accumulate diverse nongenetic, somatic damage including misfolded and aggregated proteins and cell wall defects. If the rate of damage accumulation exceeds the rate of dilution through cell growth, a dedicated mitigation strategy is required to prevent eventual population collapse. Strategies for somatic damage control can be divided into two categories, asymmetric allocation and repair, which are not, in principle, mutually exclusive. We explore a mathematical model to identify the optimal strategy, maximizing the total cell number, over a wide range of environmental and physiological conditions. The optimal strategy is primarily determined by extrinsic, damage-independent mortality and the physiological model for damage accumulation that can be either independent (linear) or increasing (exponential) with respect to the prior accumulated damage. Under the linear regime, the optimal strategy is either exclusively repair or asymmetric allocation, whereas under the exponential regime, the optimal strategy is a combination of asymmetry and repair. Repair is preferred when extrinsic mortality is low, whereas at high extrinsic mortality, asymmetric damage allocation becomes the strategy of choice. We hypothesize that at an early stage of life evolution, optimization over repair and asymmetric allocation of somatic damage gave rise to r and K selection strategists., Competing Interests: Competing interests statement:The authors declare no competing interest.

A comprehensive census of McrBC systems, among the most common forms of prokaryotic Type IV restriction systems, followed by phylogenetic analysis, reveals their enormous abundance in diverse prokaryotes and a plethora of genomic associations. We focus on a previously uncharacterized branch, which we denote co iled- co il nu clease t andems (CoCoNuTs) for their salient features: the presence of extensive coiled-coil structures and tandem nucleases. The CoCoNuTs alone show extraordinary variety, with three distinct types and multiple subtypes. All CoCoNuTs contain domains predicted to interact with translation system components, such as OB-folds resembling the SmpB protein that binds bacterial transfer-messenger RNA (tmRNA), YTH-like domains that might recognize methylated tmRNA, tRNA, or rRNA, and RNA-binding Hsp70 chaperone homologs, along with RNases, such as HEPN domains, all suggesting that the CoCoNuTs target RNA. Many CoCoNuTs might additionally target DNA, via McrC nuclease homologs. Additional restriction systems, such as Type I RM, BREX, and Druantia Type III, are frequently encoded in the same predicted superoperons. In many of these superoperons, CoCoNuTs are likely regulated by cyclic nucleotides, possibly, RNA fragments with cyclic termini, that bind associated CARF ( C RISPR- A ssociated R ossmann F old) domains. We hypothesize that the CoCoNuTs, together with the ancillary restriction factors, employ an echeloned defense strategy analogous to that of Type III CRISPR-Cas systems, in which an immune response eliminating virus DNA and/or RNA is launched first, but then, if it fails, an abortive infection response leading to PCD/dormancy via host RNA cleavage takes over., Competing Interests: RB, HS, KM, YW, EK No competing interests declared

Trypanosomatids (Euglenozoa) are a diverse group of unicellular flagellates predominately infecting insects (monoxenous species) or circulating between insects and vertebrates or plants (dixenous species). Monoxenous trypanosomatids harbor a wide range of RNA viruses belonging to the families Narnaviridae, Totiviridae , Qinviridae , Leishbuviridae, and a putative group of tombus-like viruses. Here, we focus on the subfamily Blastocrithidiinae, a previously unexplored divergent group of monoxenous trypanosomatids comprising two related genera: Obscuromonas and Blastocrithidia . Members of the genus Blastocrithidia employ a unique genetic code, in which all three stop codons are repurposed to encode amino acids, with TAA also used to terminate translation. Obscuromonas isolates studied here bear viruses of three families: Narnaviridae, Qinviridae , and Mitoviridae . The latter viral group is documented in trypanosomatid flagellates for the first time. While other known mitoviruses replicate in the mitochondria, those of trypanosomatids appear to reside in the cytoplasm. Although no RNA viruses were detected in Blastocrithidia spp., we identified an endogenous viral element in the genome of B. triatomae indicating its past encounter(s) with tombus-like viruses., Competing Interests: None declared., (© The Author(s) 2024. Published by Oxford University Press.)

In silico identification of viral anti-CRISPR proteins (Acrs) has relied largely on the guilt-by-association method using known Acrs or anti-CRISPR associated proteins (Acas) as the bait. However, the low number and limited spread of the characterized archaeal Acrs and Aca hinders our ability to identify Acrs using guilt-by-association. Here, based on the observation that the few characterized archaeal Acrs and Aca are transcribed immediately post viral infection, we hypothesize that these genes, and many other unidentified anti-defense genes (ADG), are under the control of conserved regulatory sequences including a strong promoter, which can be used to predict anti-defense genes in archaeal viruses. Using this consensus sequence based method, we identify 354 potential ADGs in 57 archaeal viruses and 6 metagenome-assembled genomes. Experimental validation identified a CRISPR subtype I-A inhibitor and the first virally encoded inhibitor of an archaeal toxin-antitoxin based immune system. We also identify regulatory proteins potentially akin to Acas that can facilitate further identification of ADGs combined with the guilt-by-association approach. These results demonstrate the potential of regulatory sequence analysis for extensive identification of ADGs in viruses of archaea and bacteria., (© 2024. The Author(s).)

Viruses with double-stranded (ds) DNA genomes in the realm Duplodnaviria share a conserved structural gene module but show a broad range of variation in their repertoires of DNA replication proteins. Some of the duplodnaviruses encode (nearly) complete replication systems whereas others lack (almost) all genes required for replication, relying on the host replication machinery. DNA polymerases (DNAPs) comprise the centerpiece of the DNA replication apparatus. The replicative DNAPs are classified into 4 unrelated or distantly related families (A-D), with the protein structures and sequences within each family being, generally, highly conserved. More than half of the duplodnaviruses encode a DNAP of family A, B or C. We showed previously that multiple pairs of closely related viruses in the order Crassvirales encode DNAPs of different families. Here we identify four additional groups of tailed phages in the class Caudoviricetes in which the DNAPs apparently were swapped on multiple occasions, with replacements occurring both between families A and B, or A and C, or between distinct subfamilies within the same family. The DNAP swapping always occurs "in situ", without changes in the organization of the surrounding genes. In several cases, the DNAP gene is the only region of substantial divergence between closely related phage genomes, whereas in others, the swap apparently involved neighboring genes encoding other proteins involved in phage replication. We hypothesize that DNAP swapping is driven by selection for avoidance of host antiphage mechanisms targeting the phage DNAP that remain to be identified, and/or by selection against replicon incompatibility. In addition, we identified two previously undetected, highly divergent groups of family A DNAPs that are encoded in some phage genomes along with the main DNAP implicated in genome replication.

Type VI CRISPR-Cas systems are among the few CRISPR varieties that target exclusively RNA. The CRISPR RNA-guided, sequence-specific binding of target RNAs, such as phage transcripts, activates the type VI effector, Cas13. Once activated, Cas13 causes collateral RNA cleavage, which induces bacterial cell dormancy, thus protecting the host population from the phage spread. We show here that the principal form of collateral RNA degradation elicited by Leptotrichia shahii Cas13a expressed in Escherichia coli cells is the cleavage of anticodons in a subset of transfer RNAs (tRNAs) with uridine-rich anticodons. This tRNA cleavage is accompanied by inhibition of protein synthesis, thus providing defense from the phages. In addition, Cas13a-mediated tRNA cleavage indirectly activates the RNases of bacterial toxin-antitoxin modules cleaving messenger RNA, which could provide a backup defense. The mechanism of Cas13a-induced antiphage defense resembles that of bacterial anticodon nucleases, which is compatible with the hypothesis that type VI effectors evolved from an abortive infection module encompassing an anticodon nuclease.

Bacterial defense against phage predation involves diverse defense systems acting individually and concurrently, yet their interactions remain poorly understood. We investigated >100 defense systems in 42,925 bacterial genomes and identified numerous instances of their non-random co-occurrence and negative association. For several pairs of defense systems significantly co-occurring in Escherichia coli strains, we demonstrate synergistic anti-phage activity. Notably, Zorya II synergizes with Druantia III and ietAS defense systems, while tmn exhibits synergy with co-occurring systems Gabija, Septu I, and PrrC. For Gabija, tmn co-opts the sensory switch ATPase domain, enhancing anti-phage activity. Some defense system pairs that are negatively associated in E. coli show synergy and significantly co-occur in other taxa, demonstrating that bacterial immune repertoires are largely shaped by selection for resistance against host-specific phages rather than negative epistasis. Collectively, these findings demonstrate compatibility and synergy between defense systems, allowing bacteria to adopt flexible strategies for phage defense., Competing Interests: Declaration of interests Y.E.G. is a full-time employee of SNIPR Biome., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Horizontal gene transfer (HGT) is a fundamental process in prokaryotic evolution, contributing significantly to diversification and adaptation. HGT is typically facilitated by mobile genetic elements (MGEs), such as conjugative plasmids and phages, which often impose fitness costs on their hosts. However, a considerable number of bacterial genes are involved in defence mechanisms that limit the propagation of MGEs, suggesting they may actively restrict HGT. In our study, we investigated whether defence systems limit HGT by examining the relationship between the HGT rate and the presence of 73 defence systems across 12 bacterial species. We discovered that only six defence systems, three of which were different CRISPR-Cas subtypes, were associated with a reduced gene gain rate at the species evolution scale. Hosts of these defence systems tend to have a smaller pangenome size and fewer phage-related genes compared to genomes without these systems. This suggests that these defence mechanisms inhibit HGT by limiting prophage integration. We hypothesize that the restriction of HGT by defence systems is species-specific and depends on various ecological and genetic factors, including the burden of MGEs and the fitness effect of HGT in bacterial populations., (Published 2024. This article is a U.S. Government work and is in the public domain in the USA. Environmental Microbiology published by John Wiley & Sons Ltd.)

Endosomal sorting complexes required for transport (ESCRT) play key roles in protein sorting between membrane-bounded compartments of eukaryotic cells. Homologs of many ESCRT components are identifiable in various groups of archaea, especially in Asgardarchaeota, the archaeal phylum that is currently considered to include the closest relatives of eukaryotes, but not in bacteria. We performed a comprehensive search for ESCRT protein homologs in archaea and reconstructed ESCRT evolution using the phylogenetic tree of Vps4 ATPase (ESCRT IV) as a scaffold and using sensitive protein sequence analysis and comparison of structural models to identify previously unknown ESCRT proteins. Several distinct groups of ESCRT systems in archaea outside of Asgard were identified, including proteins structurally similar to ESCRT-I and ESCRT-II, and several other domains involved in protein sorting in eukaryotes, suggesting an early origin of these components. Additionally, distant homologs of CdvA proteins were identified in Thermoproteales which are likely components of the uncharacterized cell division system in these archaea. We propose an evolutionary scenario for the origin of eukaryotic and Asgard ESCRT complexes from ancestral building blocks, namely, the Vps4 ATPase, ESCRT-III components, wH (winged helix-turn-helix fold) and possibly also coiled-coil, and Vps28-like domains. The last archaeal common ancestor likely encompassed a complex ESCRT system that was involved in protein sorting. Subsequent evolution involved either simplification, as in the TACK superphylum, where ESCRT was co-opted for cell division, or complexification as in Asgardarchaeota. In Asgardarchaeota, the connection between ESCRT and the ubiquitin system that was previously considered a eukaryotic signature was already established.IMPORTANCEAll eukaryotic cells possess complex intracellular membrane organization. Endosomal sorting complexes required for transport (ESCRT) play a central role in membrane remodeling which is essential for cellular functionality in eukaryotes. Recently, it has been shown that Asgard archaea, the archaeal phylum that includes the closest known relatives of eukaryotes, encode homologs of many components of the ESCRT systems. We employed protein sequence and structure comparisons to reconstruct the evolution of ESCRT systems in archaea and identified several previously unknown homologs of ESCRT subunits, some of which can be predicted to participate in cell division. The results of this reconstruction indicate that the last archaeal common ancestor already encoded a complex ESCRT system that was involved in protein sorting. In Asgard archaea, ESCRT systems evolved toward greater complexity, and in particular, the connection between ESCRT and the ubiquitin system that was previously considered a eukaryotic signature was established., Competing Interests: The authors declare no conflict of interest.

Cell division is fundamental to all cellular life. Most archaea depend on either the prokaryotic tubulin homologue FtsZ or the endosomal sorting complex required for transport for division but neither system has been robustly characterized. Here, we show that three of the four photosynthesis reaction centre barrel domain proteins of Haloferax volcanii (renamed cell division proteins B1/2/3 (CdpB1/2/3)) play important roles in cell division. CdpB1 interacts directly with the FtsZ membrane anchor SepF and is essential for cell division, whereas deletion of cdpB2 and cdpB3 causes a major and a minor division defect, respectively. Orthologues of CdpB proteins are also involved in cell division in other haloarchaea, indicating a conserved function of these proteins. Phylogenetic analysis shows that photosynthetic reaction centre barrel proteins are widely distributed among archaea and appear to be central to cell division in most if not all archaea., (© 2024. The Author(s), under exclusive licence to Springer Nature Limited.)

Competing Interests: Competing interests statement:The author declares no competing interest.

Microproteins encoded by small open reading frames (smORFs) comprise the "dark matter" of proteomes. Although functional microproteins were identified in diverse organisms from all three domains of life, bacterial smORFs remain poorly characterized. In this comprehensive study of intergenic smORFs (ismORFs, 15-70 codons) in 5,668 bacterial genomes of the family Enterobacteriaceae, we identified 67,297 clusters of ismORFs subject to purifying selection. The ismORFs mainly code for hydrophobic, potentially transmembrane, unstructured, or minimally structured microproteins. Using AlphaFold Multimer, we predicted interactions of some of the predicted microproteins encoded by transcribed ismORFs with proteins encoded by neighboring genes, revealing the potential of microproteins to regulate the activity of various proteins, particularly, under stress. We compiled a catalog of predicted microprotein families with different levels of evidence from synteny analysis, structure prediction, and transcription and translation data. This study offers a resource for investigation of biological functions of microproteins., Competing Interests: Declaration of interests The authors declare no competing interests.

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., Competing Interests: The authors declare no conflict of interest.

Metatranscriptome sequencing expanded the known diversity of the bacterial RNA virome, suggesting that additional riboviruses infecting bacterial hosts remain to be discovered. Here we employed double-stranded RNA sequencing to recover complete genome sequences of two ribovirus groups from acidic hot springs in Japan. One group, denoted hot spring riboviruses (HsRV), consists of viruses with distinct RNA-directed RNA polymerases (RdRPs) that seem to be intermediates between typical ribovirus RdRPs and viral reverse transcriptases. This group forms a distinct phylum, Artimaviricota, or even kingdom within the realm Riboviria. We identified viruses encoding HsRV-like RdRPs in marine water, river sediments and salt marshes, indicating that this group is widespread beyond extreme ecosystems. The second group, denoted hot spring partiti-like viruses (HsPV), forms a distinct branch within the family Partitiviridae. The genome architectures of HsRV and HsPV and their identification in bacteria-dominated habitats suggest that these viruses infect thermoacidophilic bacteria., (© 2024. The Author(s).)

Kolmioviridae is a family for negative-sense RNA viruses with circular, viroid-like genomes of about 1.5-1.7 kb that are maintained in mammals, amphibians, birds, fish, insects and reptiles. Deltaviruses, for instance, can cause severe hepatitis in humans. Kolmiovirids encode delta antigen (DAg) and replicate using host-cell DNA-directed RNA polymerase II and ribozymes encoded in their genome and antigenome. They require evolutionary unrelated helper viruses to provide envelopes and incorporate helper virus proteins for infectious particle formation. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the family Kolmioviridae , which is available at ictv.global/report/kolmioviridae.

The archaeal ancestor of eukaryotes apparently belonged to the phylum Asgardarchaeota, but the ecology and evolution of Asgard archaea are poorly understood. The optimal GDP-binding temperature of a translation elongation factor (EF-1A or EF-Tu) has been previously shown to correlate with the optimal growth temperature of diverse prokaryotes. Here, we reconstruct ancestral EF-1A sequences and experimentally measure the optimal GDP-binding temperature of EF-1A from ancient and extant Asgard archaea, to infer the evolution of optimal growth temperatures in Asgardarchaeota. Our results suggest that the Asgard ancestor of eukaryotes was a moderate thermophile, with an optimal growth temperature around 53 °C. The origin of eukaryotes appears to coincide with a transition from thermophilic to mesophilic lifestyle during the evolution of Asgard archaea., (© 2024. The Author(s).)

Bathyarchaeia represent a class of archaea common and abundant in sedimentary ecosystems. Here we report 56 metagenome-assembled genomes of Bathyarchaeia viruses identified in metagenomes from different environments. Gene sharing network and phylogenomic analyses led to the proposal of four virus families, including viruses of the realms Duplodnaviria and Adnaviria , and archaea-specific spindle-shaped viruses. Genomic analyses uncovered diverse CRISPR elements in these viruses. Viruses of the proposed family " Fuxiviridae " harbor an atypical Type IV-B CRISPR-Cas system and a Cas4 protein that might interfere with host immunity. Viruses of the family " Chiyouviridae " encode a Cas2-like endonuclease and two mini-CRISPR arrays, one with a repeat identical to that in the host CRISPR array, potentially allowing the virus to recruit the host CRISPR adaptation machinery to acquire spacers that could contribute to competition with other mobile genetic elements or to inhibit host defenses. These findings present an outline of the Bathyarchaeia virome and offer a glimpse into their counter-defense mechanisms., Competing Interests: The authors declare no competing interests., (© The Author(s) 2024. Published by Oxford University Press on behalf of the International Society for Microbial Ecology.)

Nearly all organisms are hosts to multiple viruses that collectively appear to be the most abundant biological entities in the biosphere. With recent advances in metagenomics and metatranscriptomics, the known diversity of viruses substantially expanded. Comparative analysis of these viruses using advanced computational methods culminated in the reconstruction of the evolution of major groups of viruses and enabled the construction of a virus megataxonomy, which has been formally adopted by the International Committee on Taxonomy of Viruses. This comprehensive taxonomy consists of six virus realms, which are aspired to be monophyletic and assembled based on the conservation of hallmark proteins involved in capsid structure formation or genome replication. The viruses in different major taxa substantially differ in host range and accordingly in ecological niches. In this review article, we outline the latest developments in virus megataxonomy and the recent discoveries that will likely lead to reassessment of some major taxa, in particular, split of three of the current six realms into two or more independent realms. We then discuss the correspondence between virus taxonomy and the distribution of viruses among hosts and ecological niches, as well as the abundance of viruses versus cells in different habitats. The distribution of viruses across environments appears to be primarily determined by the host ranges, i.e. the virome is shaped by the composition of the biome in a given habitat, which itself is affected by abiotic factors., (© The Author(s) 2024. Published by Oxford University Press on behalf of the International Society for Microbial Ecology.)

Thermus thermophilus bacteriophage P23-45 encodes a giant 5,002-residue tail tape measure protein (TMP) that defines the length of its extraordinarily long tail. Here, we show that the N-terminal portion of P23-45 TMP is an unusual RNA polymerase (RNAP) homologous to cellular RNAPs. The TMP-fused virion RNAP transcribes pre-early phage genes, including a gene that encodes another, non-virion RNAP, that transcribes early and some middle phage genes. We report the crystal structures of both P23-45 RNAPs. The non-virion RNAP has a crab-claw-like architecture. By contrast, the virion RNAP adopts a unique flat structure without a clamp. Structure and sequence comparisons of the P23-45 RNAPs with other RNAPs suggest that, despite the extensive functional differences, the two P23-45 RNAPs originate from an ancient gene duplication in an ancestral phage. Our findings demonstrate striking adaptability of RNAPs that can be attained within a single virus species., (© 2024. The Author(s).)

Although social interactions are known to drive pathogen transmission, the contributions of socially transmissible host-associated mutualists and commensals to host health and disease remain poorly explored. We use the concept of the social microbiome-the microbial metacommunity of a social network of hosts-to analyze the implications of social microbial transmission for host health and disease. We investigate the contributions of socially transmissible microbes to both eco-evolutionary microbiome community processes (colonization resistance, the evolution of virulence, and reactions to ecological disturbance) and microbial transmission-based processes (transmission of microbes with metabolic and immune effects, inter-specific transmission, transmission of antibiotic-resistant microbes, and transmission of viruses). We consider the implications of social microbial transmission for communicable and non-communicable diseases and evaluate the importance of a socially transmissible component underlying canonically non-communicable diseases. The social transmission of mutualists and commensals may play a significant, under-appreciated role in the social determinants of health and may act as a hidden force in social evolution., Competing Interests: Declaration of interests A.S., C.J.A.M., A.R., N.G.O.I., K.V.-A.J., M.V.-C., I.L.B., J.H., E.A.A., L.B.B., B.B.F., R.N.C., and A.H.M. declare no competing interests. S.H. is a current employee of SilverCloud Health and SilverCloud Health is a subsidiary of Amwell. S.H. holds shares in Amwell. F.S.G. is listed as an inventor on patent applications pertaining to the treatment of cancer or the study of the human microbiome. E.V.K. is listed as an inventor on patents and patent applications pertaining to CRISPR enzymes and systems or protein regulation., (Copyright © 2023 Elsevier Inc. All rights reserved.)

CRISPR systems mediate adaptive immunity in bacteria and archaea through diverse effector mechanisms and have been repurposed for versatile applications in therapeutics and diagnostics thanks to their facile reprogramming with RNA guides. RNA-guided CRISPR-Cas targeting and interference are mediated by effectors that are either components of multisubunit complexes in class 1 systems or multidomain single-effector proteins in class 2. The compact class 2 CRISPR systems have been broadly adopted for multiple applications, especially genome editing, leading to a transformation of the molecular biology and biotechnology toolkit. The diversity of class 2 effector enzymes, initially limited to the Cas9 nuclease, was substantially expanded via computational genome and metagenome mining to include numerous variants of Cas12 and Cas13, providing substrates for the development of versatile, orthogonal molecular tools. Characterization of these diverse CRISPR effectors uncovered many new features, including distinct protospacer adjacent motifs (PAMs) that expand the targeting space, improved editing specificity, RNA rather than DNA targeting, smaller crRNAs, staggered and blunt end cuts, miniature enzymes, promiscuous RNA and DNA cleavage, etc. These unique properties enabled multiple applications, such as harnessing the promiscuous RNase activity of the type VI effector, Cas13, for supersensitive nucleic acid detection. class 1 CRISPR systems have been adopted for genome editing, as well, despite the challenge of expressing and delivering the multiprotein class 1 effectors. The rich diversity of CRISPR enzymes led to rapid maturation of the genome editing toolbox, with capabilities such as gene knockout, base editing, prime editing, gene insertion, DNA imaging, epigenetic modulation, transcriptional modulation, and RNA editing. Combined with rational design and engineering of the effector proteins and associated RNAs, the natural diversity of CRISPR and related bacterial RNA-guided systems provides a vast resource for expanding the repertoire of tools for molecular biology and biotechnology.

Polintons are double-stranded DNA, virus-like self-synthesizing transposons widely found in eukaryotic genomes. Recent metagenomic discoveries of Polinton-like viruses are consistent with the hypothesis that Polintons invade eukaryotic host genomes through infectious viral particles. Nematode genomes contain multiple copies of Polintons and provide an opportunity to explore the natural distribution and evolution of Polintons during this process. We performed an extensive search of Polintons across nematode genomes, identifying multiple full-length Polinton copies in several species. We provide evidence of both ancient Polinton integrations and recent mobility in strains of the same nematode species. In addition to the major nematode Polinton family, we identified a group of Polintons that are overall closely related to the major family but encode a distinct protein-primed DNA polymerase B (pPolB) that is related to homologs from a different group of Polintons present outside of the Nematoda. Phylogenetic analyses on the pPolBs support the evolutionary scenarios in which these extrinsic pPolBs that seem to derive from Polinton families present in oomycetes and molluscs replaced the canonical pPolB in subsets of Polintons found in terrestrial and marine nematodes, respectively, suggesting interphylum horizontal gene transfers. The pPolBs of the terrestrial nematode and oomycete Polintons share a unique feature, an insertion of an HNH nuclease domain, whereas the pPolBs in the marine nematode Polintons share an insertion of a VSR nuclease domain with marine mollusc pPolBs. We hypothesize that horizontal gene transfer occurs among Polintons from widely different but cohabiting hosts., (© The Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Most life forms harbor multiple, diverse mobile genetic elements (MGE) that widely differ in their rates and mechanisms of mobility. Recent findings on two classes of MGE in prokaryotes revealed a novel mechanism, RNA-guided transposition, where a transposon-encoded guide RNA directs the transposase to a unique site in the host genome. Tn7-like transposons, on multiple occasions, recruited CRISPR systems that lost the capacity to cleave target DNA and instead mediate RNA-guided transposition via CRISPR RNA. Conversely, the abundant transposon-associated, RNA-guided nucleases IscB and TnpB that appear to promote proliferation of IS 200 /IS 605 and IS 607 transposons were the likely evolutionary ancestors of type II and type V CRISPR systems, respectively. Thus, RNA-guided target recognition is a major biological phenomenon that connects MGE with host defense mechanisms. More RNA-guided defensive and MGE-associated functionalities are likely to be discovered., Competing Interests: Conflict of interest statement The authors declare they have no conflict of interest.

The TnpB proteins are transposon-associated RNA-guided nucleases that are among the most abundant proteins encoded in bacterial and archaeal genomes, but whose functions in the transposon life cycle remain unknown. TnpB appears to be the evolutionary ancestor of Cas12, the effector nuclease of type V CRISPR-Cas systems. We performed a comprehensive census of TnpBs in archaeal and bacterial genomes and constructed a phylogenetic tree on which we mapped various features of these proteins. In multiple branches of the tree, the catalytic site of the TnpB nuclease is rearranged, demonstrating structural and probably biochemical malleability of this enzyme. We identified numerous cases of apparent recruitment of TnpB for other functions of which the most common is the evolution of type V CRISPR-Cas effectors on about 50 independent occasions. In many other cases of more radical exaptation, the catalytic site of the TnpB nuclease is apparently inactivated, suggesting a regulatory function, whereas in others, the activity appears to be retained, indicating that the recruited TnpB functions as a nuclease, for example, as a toxin. These findings demonstrate remarkable evolutionary malleability of the TnpB scaffold and provide extensive opportunities for further exploration of RNA-guided biological systems as well as multiple applications., Competing Interests: Competing interests statement:F.Z. is a scientific advisor and cofounder of Editas Medicine, Beam Therapeutics, Pairwise Plants, Arbor Biotechnologies, and Aera Therapeutics. F.Z. is also a scientific advisor for Octant.

Microbial systems underpin many biotechnologies, including CRISPR, but the exponential growth of sequence databases makes it difficult to find previously unidentified systems. In this work, we develop the fast locality-sensitive hashing-based clustering (FLSHclust) algorithm, which performs deep clustering on massive datasets in linearithmic time. We incorporated FLSHclust into a CRISPR discovery pipeline and identified 188 previously unreported CRISPR-linked gene modules, revealing many additional biochemical functions coupled to adaptive immunity. We experimentally characterized three HNH nuclease-containing CRISPR systems, including the first type IV system with a specified interference mechanism, and engineered them for genome editing. We also identified and characterized a candidate type VII system, which we show acts on RNA. This work opens new avenues for harnessing CRISPR and for the broader exploration of the vast functional diversity of microbial proteins.

Genomes of bacteria and archaea contain a much larger fraction of unidirectional (serial) gene pairs than convergent or divergent gene pairs. Many of the unidirectional gene pairs have short overlaps of -4 nt and -1 nt. As shown previously, translation of the genes in overlapping unidirectional gene pairs is tightly coupled. Two alternative models for the fate of the post-termination ribosome predict either that overlaps or very short intergenic distances are essential for translational coupling or that the undissociated post-termination ribosome can scan through long intergenic regions, up to hundreds of nucleotides. We aimed to experimentally resolve the contradiction between the two models by analyzing three native gene pairs from the model archaeon Haloferax volcanii and three native pairs from Escherichia coli . A two reporter gene system was used to quantify the reinitiation frequency, and several stop codons in the upstream gene were introduced to increase the intergenic distances. For all six gene pairs from two species, an extremely strong dependence of the reinitiation efficiency on the intergenic distance was unequivocally demonstrated, such that even short intergenic distances of about 20 nt almost completely abolished translational coupling. Bioinformatic analysis of the intergenic distances in all unidirectional gene pairs in the genomes of H. volcanii and E. coli and in 1,695 prokaryotic species representative of 49 phyla showed that intergenic distances of -4 nt or -1 nt (= short gene overlaps of 4 nt or 1 nt) were by far most common in all these groups of archaea and bacteria. A small set of genes in E. coli , but not in H. volcanii , had intergenic distances of around +10 nt. Our experimental and bioinformatic analyses clearly show that translational coupling requires short gene overlaps, whereas scanning of intergenic regions by the post-termination ribosome occurs rarely, if at all. Short overlaps are enriched among genes that encode subunits of heteromeric complexes, and co-translational complex formation requiring precise subunit stoichiometry likely confers an evolutionary advantage that drove the formation and conservation of overlapping gene pairs during evolution., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision., (Copyright © 2023 Huber, Vogel, Borst, Pfeiffer, Karamycheva, Wolf, Koonin and Soppa.)

Yueviridae is a family of negative-sense RNA viruses with genomes of 7.8-8.2 kb that have been associated with crustaceans, insects, stramenopiles and plants. The yuevirid genome consist of two segments, each with at least one ORF. The large (L) segment ORF encodes a large protein containing an RNA-directed RNA polymerase domain. The small (S) segment ORF encodes a nucleocapsid protein. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the family Yueviridae , which is available at http://www.ictv.global/report/yueviridae.

Amnoonviridae is a family of negative-sense RNA viruses with genomes totalling about 10.3 kb. These viruses have been found in fish. The amnoonvirid genome consists of 10 segments, each with at least 1 open reading frame (ORF). The RNA1-3 ORFs encode the three subunits of the viral polymerase. The RNA4 ORF encodes a nucleoprotein. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the family Amnoonviridae , which is available at ictv.global/report/amnoonviridae.

Tosoviridae is a family of negative-sense RNA viruses with genomes totaling about 12.3 kb that have been found in turtles. The tosovirid genome consists of two segments, each with two open reading frames (ORFs) in ambisense orientation. The small (S) segment encodes a nucleoprotein (NP) and a glycoprotein precursor (GPC); the large (L) segment encodes an L protein containing an RNA-directed RNA polymerase (RdRP) domain and a zinc-binding (Z) protein. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the family Tosoviridae , which is available at ictv.global/report/tosoviridae.