Mining metatranscriptomes reveals a vast world of viroid-like circular RNAs.

Viroids and viroid-like covalently closed circular (ccc) RNAs are minimal replicators that typically encode no proteins and hijack cellular enzymes for replication. The extent and diversity of viroid-like agents are poorly understood. We developed a computational pipeline to identify viroid-like cccRNAs and applied it to 5,131 metatranscriptomes and 1,344 plant transcriptomes. The search yielded 11,378 viroid-like cccRNAs spanning 4,409 species-level clusters, a 5-fold increase compared to the previously identified viroid-like elements. Within this diverse collection, we discovered numerous putative viroids, satellite RNAs, retrozymes, and ribozy-like viruses. Diverse ribozyme combinations and unusual ribozymes within the cccRNAs were identified. Self-cleaving ribozymes were identified in ambiviruses, some mito-like viruses and capsid-encoding satellite virus-like cccRNAs. The broad presence of viroid-like cccRNAs in diverse transcriptomes and ecosystems implies that their host range is far broader than currently known, and matches to CRISPR spacers suggest that some cccRNAs replicate in prokaryotes.

Expansion of the global RNA virome reveals diverse clades of bacteriophages.

High-throughput RNA sequencing offers broad opportunities to explore the Earth RNA virome. Mining 5,150 diverse metatranscriptomes uncovered >2.5 million RNA virus contigs. Analysis of >330,000 RNA-dependent RNA polymerases (RdRPs) shows that this expansion corresponds to a 5-fold increase of the known RNA virus diversity. Gene content analysis revealed multiple protein domains previously not found in RNA viruses and implicated in virus-host interactions. Extended RdRP phylogeny supports the monophyly of the five established phyla and reveals two putative additional bacteriophage phyla and numerous putative additional classes and orders. The dramatically expanded phylum Lenarviricota, consisting of bacterial and related eukaryotic viruses, now accounts for a third of the RNA virome. Identification of CRISPR spacer matches and bacteriolytic proteins suggests that subsets of picobirnaviruses and partitiviruses, previously associated with eukaryotes, infect prokaryotic hosts.

Modularity and diversity of target selectors in Tn7 transposons.

To spread, transposons must integrate into target sites without disruption of essential genes while avoiding host defense systems. Tn7-like transposons employ multiple mechanisms for target-site selection, including protein-guided targeting and, in CRISPR-associated transposons (CASTs), RNA-guided targeting. Combining phylogenomic and structural analyses, we conducted a broad survey of target selectors, revealing diverse mechanisms used by Tn7 to recognize target sites, including previously uncharacterized target-selector proteins found in newly discovered transposable elements (TEs). We experimentally characterized a CAST I-D system and a Tn6022-like transposon that uses TnsF, which contains an inactivated tyrosine recombinase domain, to target the comM gene. Additionally, we identified a non-Tn7 transposon, Tsy, encoding a homolog of TnsF with an active tyrosine recombinase domain, which we show also inserts into comM. Our findings show that Tn7 transposons employ modular architecture and co-opt target selectors from various sources to optimize target selection and drive transposon spread.

Small proteins can no longer be ignored.

Small proteins, here defined as proteins of 50 amino acids or fewer in the absence of processing, have traditionally been overlooked due to challenges in their annotation and biochemical detection. In the past several years, however, increasing numbers of small proteins have been identified either through the realization that mutations in intergenic regions are actually within unannotated small protein genes or through the discovery that some small, regulatory RNAs encode small proteins. These insights, together with comparative sequence analysis, indicate that tens if not hundreds of small proteins are synthesized in a given organism. This review summarizes what has been learned about the functions of several of these bacterial small proteins, most of which act at the membrane, illustrating the astonishing range of processes in which these small proteins act and suggesting several general conclusions. Important questions for future studies of these overlooked proteins are also discussed.

Mutation-selection balance and compensatory mechanisms in tumour evolution.

Intratumour heterogeneity and phenotypic plasticity, sustained by a range of somatic aberrations, as well as epigenetic and metabolic adaptations, are the principal mechanisms that enable cancers to resist treatment and survive under environmental stress. A comprehensive picture of the interplay between different somatic aberrations, from point mutations to whole-genome duplications, in tumour initiation and progression is lacking. We posit that different genomic aberrations generally exhibit a temporal order, shaped by a balance between the levels of mutations and selective pressures. Repeat instability emerges first, followed by larger aberrations, with compensatory effects leading to robust tumour fitness maintained throughout the tumour progression. A better understanding of the interplay between genetic aberrations, the microenvironment, and epigenetic and metabolic cellular states is essential for early detection and prevention of cancer as well as development of efficient therapeutic strategies.

Expanded diversity of Asgard archaea and their relationships with eukaryotes.

Asgard is a recently discovered superphylum of archaea that appears to include the closest archaeal relatives of eukaryotes1-5. Debate continues as to whether the archaeal ancestor of eukaryotes belongs within the Asgard superphylum or whether this ancestor is a sister group to all other archaea (that is, a two-domain versus a three-domain tree of life)6-8. Here we present a comparative analysis of 162 complete or nearly complete genomes of Asgard archaea, including 75 metagenome-assembled genomes that-to our knowledge-have not previously been reported. Our results substantially expand the phylogenetic diversity of Asgard and lead us to propose six additional phyla that include a deep branch that we have provisionally named Wukongarchaeota. Our phylogenomic analysis does not resolve unequivocally the evolutionary relationship between eukaryotes and Asgard archaea, but instead-depending on the choice of species and conserved genes used to build the phylogeny-supports either the origin of eukaryotes from within Asgard (as a sister group to the expanded Heimdallarchaeota-Wukongarchaeota branch) or a deeper branch for the eukaryote ancestor within archaea. Our comprehensive protein domain analysis using the 162 Asgard genomes results in a major expansion of the set of eukaryotic signature proteins. The Asgard eukaryotic signature proteins show variable phyletic distributions and domain architectures, which is suggestive of dynamic evolution through horizontal gene transfer, gene loss, gene duplication and domain shuffling. The phylogenomics of the Asgard archaea points to the accumulation of the components of the mobile archaeal 'eukaryome' in the archaeal ancestor of eukaryotes (within or outside Asgard) through extensive horizontal gene transfer.

Balance between asymmetric allocation and repair of somatic damage in unicellular life forms as an ancient form of r/K selection.

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.

Evolutionary entanglement of mobile genetic elements and host defence systems: guns for hire.

All cellular life forms are afflicted by diverse genetic parasites, including viruses and other types of mobile genetic elements (MGEs), and have evolved multiple, diverse defence systems that protect them from MGE assault via different mechanisms. Here, we provide our perspectives on how recent evidence points to tight evolutionary connections between MGEs and defence systems that reach far beyond the proverbial arms race. Defence systems incur a fitness cost for the hosts; therefore, at least in prokaryotes, horizontal mobility of defence systems, mediated primarily by MGEs, is essential for their persistence. Moreover, defence systems themselves possess certain features of selfish elements. Common components of MGEs, such as site-specific nucleases, are 'guns for hire' that can also function as parts of defence mechanisms and are often shuttled between MGEs and defence systems. Thus, evolutionary and molecular factors converge to mould the multifaceted, inextricable connection between MGEs and anti-MGE defence systems.

Compensatory relationship between low-complexity regions and gene paralogy in the evolution of prokaryotes.

The evolution of genomes in all life forms involves two distinct, dynamic types of genomic changes: gene duplication (and loss) that shape families of paralogous genes and extension (and contraction) of low-complexity regions (LCR), which occurs through dynamics of short repeats in protein-coding genes. Although the roles of each of these types of events in genome evolution have been studied, their co-evolutionary dynamics is not thoroughly understood. Here, by analyzing a wide range of genomes from diverse bacteria and archaea, we show that LCR and paralogy represent two distinct routes of evolution that are inversely correlated. The emergence of LCR is a prominent evolutionary mechanism in fast evolving, young protein families, whereas paralogy dominates the comparatively slow evolution of old protein families. The analysis of multiple prokaryotic genomes shows that the formation of LCR is likely a widespread, transient evolutionary mechanism that temporally and locally affects also ancestral functions, but apparently, fades away with time, under mutational and selective pressures, yielding to gene paralogy. We propose that compensatory relationships between short-term and longer-term evolutionary mechanisms are universal in the evolution of life.

Diversity, evolution, and classification of the RNA-guided nucleases TnpB and Cas12.

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.

Coevolution of reproducers and replicators at the origin of life and the conditions for the origin of genomes.

There are two fundamentally distinct but inextricably linked types of biological evolutionary units, reproducers and replicators. Reproducers are cells and organelles that reproduce via various forms of division and maintain the physical continuity of compartments and their content. Replicators are genetic elements (GE), including genomes of cellular organisms and various autonomous elements, that both cooperate with reproducers and rely on the latter for replication. All known cells and organisms comprise a union between replicators and reproducers. We explore a model in which cells emerged via symbiosis between primordial "metabolic" reproducers (protocells) which evolved, on short time scales, via a primitive form of selection and random drift, and mutualist replicators. Mathematical modeling identifies the conditions, under which GE-carrying protocells can outcompete GE-less ones, taking into account that, from the earliest stages of evolution, replicators split into mutualists and parasites. Analysis of the model shows that, for the GE-containing protocells to win the competition and to be fixed in evolution, it is essential that the birth-death process of the GE is coordinated with the rate of protocell division. At the early stages of evolution, random, high-variance cell division is advantageous compared with symmetrical division because the former provides for the emergence of protocells containing only mutualists, preventing takeover by parasites. These findings illuminate the likely order of key events on the evolutionary route from protocells to cells that involved the origin of genomes, symmetrical cell division, and antiparasite defense systems.

Widespread CRISPR-derived RNA regulatory elements in CRISPR-Cas systems.

CRISPR-cas loci typically contain CRISPR arrays with unique spacers separating direct repeats. Spacers along with portions of adjacent repeats are transcribed and processed into CRISPR(cr) RNAs that target complementary sequences (protospacers) in mobile genetic elements, resulting in cleavage of the target DNA or RNA. Additional, standalone repeats in some CRISPR-cas loci produce distinct cr-like RNAs implicated in regulatory or other functions. We developed a computational pipeline to systematically predict crRNA-like elements by scanning for standalone repeat sequences that are conserved in closely related CRISPR-cas loci. Numerous crRNA-like elements were detected in diverse CRISPR-Cas systems, mostly, of type I, but also subtype V-A. Standalone repeats often form mini-arrays containing two repeat-like sequence separated by a spacer that is partially complementary to promoter regions of cas genes, in particular cas8, or cargo genes located within CRISPR-Cas loci, such as toxins-antitoxins. We show experimentally that a mini-array from a type I-F1 CRISPR-Cas system functions as a regulatory guide. We also identified mini-arrays in bacteriophages that could abrogate CRISPR immunity by inhibiting effector expression. Thus, recruitment of CRISPR effectors for regulatory functions via spacers with partial complementarity to the target is a common feature of diverse CRISPR-Cas systems.

Thermodynamics of evolution and the origin of life.

We outline a phenomenological theory of evolution and origin of life by combining the formalism of classical thermodynamics with a statistical description of learning. The maximum entropy principle constrained by the requirement for minimization of the loss function is employed to derive a canonical ensemble of organisms (population), the corresponding partition function (macroscopic counterpart of fitness), and free energy (macroscopic counterpart of additive fitness). We further define the biological counterparts of temperature (evolutionary temperature) as the measure of stochasticity of the evolutionary process and of chemical potential (evolutionary potential) as the amount of evolutionary work required to add a new trainable variable (such as an additional gene) to the evolving system. We then develop a phenomenological approach to the description of evolution, which involves modeling the grand potential as a function of the evolutionary temperature and evolutionary potential. We demonstrate how this phenomenological approach can be used to study the "ideal mutation" model of evolution and its generalizations. Finally, we show that, within this thermodynamics framework, major transitions in evolution, such as the transition from an ensemble of molecules to an ensemble of organisms, that is, the origin of life, can be modeled as a special case of bona fide physical phase transitions that are associated with the emergence of a new type of grand canonical ensemble and the corresponding new level of description.

Toward a theory of evolution as multilevel learning.

We apply the theory of learning to physically renormalizable systems in an attempt to outline a theory of biological evolution, including the origin of life, as multilevel learning. We formulate seven fundamental principles of evolution that appear to be necessary and sufficient to render a universe observable and show that they entail the major features of biological evolution, including replication and natural selection. It is shown that these cornerstone phenomena of biology emerge from the fundamental features of learning dynamics such as the existence of a loss function, which is minimized during learning. We then sketch the theory of evolution using the mathematical framework of neural networks, which provides for detailed analysis of evolutionary phenomena. To demonstrate the potential of the proposed theoretical framework, we derive a generalized version of the Central Dogma of molecular biology by analyzing the flow of information during learning (back propagation) and predicting (forward propagation) the environment by evolving organisms. The more complex evolutionary phenomena, such as major transitions in evolution (in particular, the origin of life), have to be analyzed in the thermodynamic limit, which is described in detail in the paper by Vanchurin et al. [V. Vanchurin, Y. I. Wolf, E. V. Koonin, M. I. Katsnelson, Proc. Natl. Acad. Sci. U.S.A. 119, 10.1073/pnas.2120042119 (2022)].

Human pathogenic RNA viruses establish noncompeting lineages by occupying independent niches.

Many pathogenic viruses are endemic among human populations and can cause a broad variety of diseases, some potentially leading to devastating pandemics. How virus populations maintain diversity and what selective pressures drive population turnover is not thoroughly understood. We conducted a large-scale phylodynamic analysis of 27 human pathogenic RNA viruses spanning diverse life history traits, in search of unifying trends that shape virus evolution. For most virus species, we identify multiple, cocirculating lineages with low turnover rates. These lineages appear to be largely noncompeting and likely occupy semiindependent epidemiological niches that are not regionally or seasonally defined. Typically, intralineage mutational signatures are similar to interlineage signatures. The principal exception are members of the family Picornaviridae, for which mutations in capsid protein genes are primarily lineage defining. Interlineage turnover is slower than expected under a neutral model, whereas intralineage turnover is faster than the neutral expectation, further supporting the existence of independent niches. The persistence of virus lineages appears to stem from limited outbreaks within small communities, so that only a small fraction of the global susceptible population is infected at any time. As disparate communities become increasingly connected through globalization, interaction and competition between lineages might increase as well, which could result in changing selective pressures and increased diversification and/or pathogenicity. Thus, in addition to zoonotic events, ongoing surveillance of familiar, endemic viruses appears to merit global attention with respect to the prevention or mitigation of future pandemics.

Microbial diversity and ecological complexity emerging from environmental variation and horizontal gene transfer in a simple mathematical model.

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.

Defence systems and horizontal gene transfer in bacteria.

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.

A novel group of negative-sense RNA viruses associated with epizootics in managed and free-ranging freshwater turtles in Florida, USA.

Few aquatic animal negative-sense RNA viruses have been characterized, and their role in disease is poorly understood. Here, we describe a virus isolated from diseased freshwater turtles from a Florida farm in 2007 and from an ongoing epizootic among free-ranging populations of Florida softshell turtles (Apalone ferox), Florida red-bellied cooters (Pseudemys nelsoni), and peninsula cooters (Pseudemys peninsularis). Affected turtles presented with similar neurological signs, oral and genital ulceration, and secondary microbial infections. Microscopic lesions were most severe in the softshell turtles and included heterophilic/histiocytic meningoencephalitis, multi-organ vasculitis, and cytologic observation of leukocytic intracytoplasmic inclusions. The virus was isolated using Terrapene heart (TH-1) cells. Ultrastructurally, viral particles were round to pleomorphic and acquired an envelope with prominent surface projections by budding from the cell membrane. Viral genomes were sequenced from cDNA libraries of two nearly identical isolates and determined to be bi-segmented, with an ambisense coding arrangement. The larger segment encodes a predicted RNA-directed RNA polymerase (RdRP) and a putative zinc-binding matrix protein. The smaller segment encodes a putative nucleoprotein and an envelope glycoprotein precursor (GPC). Thus, the genome organization of this turtle virus resembles that of arenaviruses. Phylogenetic analysis shows that the RdRP of the turtle virus is highly diverged from the RdRPs of all known negative-sense RNA viruses and forms a deep branch within the phylum Negarnaviricota, that is not affiliated with any known group of viruses, even at the class level. In contrast, the GPC protein of the turtle virus is confidently affiliated with homologs from a distinct group of fish hantaviruses. Thus, the turtle virus is expected to become the founder of a new taxon of negative-sense RNA viruses, at least with a family rank, but likely, an order or even a class. These viruses probably evolved either by reassortment or by intrasegment recombination between a virus from a distinct branch of negarnaviruses distant from all known groups and a hanta-like aquatic virus. We suggest the provisional name Tosoviridae for the putative new family, with Turtle fraservirus 1 (TFV1) as the type species within the genus Fraservirus. A conventional RT-PCR assay, targeting the TFV1 RdRP, confirmed the presence of viral RNA in multiple tissues and exudates from diseased turtles. The systemic nature of the TFV1 infection was further supported by labeling of cells within lesions using in situ hybridization targeting the RNA of the TFV1 RdRP.

Molecular adaptations during viral epidemics.

In 1977, the world witnessed both the eradication of smallpox and the beginning of the modern age of genomics. Over the following half-century, 7 epidemic viruses of international concern galvanized virologists across the globe and led to increasingly extensive virus genome sequencing. These sequencing efforts exerted over periods of rapid adaptation of viruses to new hosts, in particular, humans provide insight into the molecular mechanisms underpinning virus evolution. Investment in virus genome sequencing was dramatically increased by the unprecedented support for phylogenomic analyses during the COVID-19 pandemic. In this review, we attempt to piece together comprehensive molecular histories of the adaptation of variola virus, HIV-1 M, SARS, H1N1-SIV, MERS, Ebola, Zika, and SARS-CoV-2 to the human host. Disruption of genes involved in virus-host interaction in animal hosts, recombination including genome segment reassortment, and adaptive mutations leading to amino acid replacements in virus proteins involved in host receptor binding and membrane fusion are identified as the key factors in the evolution of epidemic viruses.

The tRNA discriminator base defines the mutual orthogonality of two distinct pyrrolysyl-tRNA synthetase/tRNAPyl pairs in the same organism.

Site-specific incorporation of distinct non-canonical amino acids into proteins via genetic code expansion requires mutually orthogonal aminoacyl-tRNA synthetase/tRNA pairs. Pyrrolysyl-tRNA synthetase (PylRS)/tRNAPyl pairs are ideal for genetic code expansion and have been extensively engineered for developing mutually orthogonal pairs. Here, we identify two novel wild-type PylRS/tRNAPyl pairs simultaneously present in the deep-rooted extremely halophilic euryarchaeal methanogen Candidatus Methanohalarchaeum thermophilum HMET1, and show that both pairs are functional in the model halophilic archaeon Haloferax volcanii. These pairs consist of two different PylRS enzymes and two distinct tRNAs with dissimilar discriminator bases. Surprisingly, these two PylRS/tRNAPyl pairs display mutual orthogonality enabled by two unique features, the A73 discriminator base of tRNAPyl2 and a shorter motif 2 loop in PylRS2. In vivo translation experiments show that tRNAPyl2 charging by PylRS2 is defined by the enzyme's shortened motif 2 loop. Finally, we demonstrate that the two HMET1 PylRS/tRNAPyl pairs can simultaneously decode UAG and UAA codons for incorporation of two distinct noncanonical amino acids into protein. This example of a single base change in a tRNA leading to additional coding capacity suggests that the growth of the genetic code is not yet limited by the number of identity elements fitting into the tRNA structure.