Search

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.

CrAssphage and related viruses of the order Crassvirales(hereafter referred to as crassviruses) were originally discovered by cross-assembly of metagenomic sequences. They are the most abundant viruses in the human gut, are found in the majority of individual gut viromes, and account for up to 95% of the viral sequences in some individuals1–4. Crassviruses are likely to have major roles in shaping the composition and functionality of the human microbiome, but the structures and roles of most of the virally encoded proteins are unknown, with only generic predictions resulting from bioinformatic analyses4,5. Here we present a cryo-electron microscopy reconstruction of Bacteroides intestinalisvirus ΦcrAss0016, providing the structural basis for the functional assignment of most of its virion proteins. The muzzle protein forms an assembly about 1 MDa in size at the end of the tail and exhibits a previously unknown fold that we designate the ‘crass fold’, that is likely to serve as a gatekeeper that controls the ejection of cargos. In addition to packing the approximately 103 kb of virus DNA, the ΦcrAss001 virion has extensive storage space for virally encoded cargo proteins in the capsid and, unusually, within the tail. One of the cargo proteins is present in both the capsid and the tail, suggesting a general mechanism for protein ejection, which involves partial unfolding of proteins during their extrusion through the tail. These findings provide a structural basis for understanding the mechanisms of assembly and infection of these highly abundant crassviruses.

The COVID-19 pandemic highlighted the need to understand the emergence of viral variants, given that these can have implications for vaccination success. A bioinformatics tool offers a way to predict viral evolution.

RNA-guided systems, such as CRISPR–Cas, combine programmable substrate recognition with enzymatic function, a combination that has been used advantageously to develop powerful molecular technologies1,2. Structural studies of these systems have illuminated how the RNA and protein jointly recognize and cleave their substrates, guiding rational engineering for further technology development3. Recent work identified a new class of RNA-guided systems, termed OMEGA, which include IscB, the likely ancestor of Cas9, and the nickase IsrB, a homologue of IscB lacking the HNH nuclease domain4. IsrB consists of only around 350 amino acids, but its small size is counterbalanced by a relatively large RNA guide (roughly 300-nt ωRNA). Here, we report the cryogenic-electron microscopy structure of Desulfovirgula thermocuniculiIsrB (DtIsrB) in complex with its cognate ωRNA and a target DNA. We find the overall structure of the IsrB protein shares a common scaffold with Cas9. In contrast to Cas9, however, which uses a recognition (REC) lobe to facilitate target selection, IsrB relies on its ωRNA, part of which forms an intricate ternary structure positioned analogously to REC. Structural analyses of IsrB and its ωRNA as well as comparisons to other RNA-guided systems highlight the functional interplay between protein and RNA, advancing our understanding of the biology and evolution of these diverse systems.

CRISPR–Cas interference is mediated by Cas effector nucleases that are either components of multisubunit complexes—in class 1 CRISPR–Cas systems—or domains of a single protein—in class 2 systems1–3. Here we show that the subtype III-E effector Cas7-11 is a single-protein effector in the class 1 CRISPR–Cas systems originating from the fusion of a putative Cas11 domain and multiple Cas7 subunits that are derived from subtype III-D. Cas7-11 from Desulfonema ishimotonii(DiCas7-11), when expressed in Escherichia coli, has substantial RNA interference effectivity against mRNAs and bacteriophages. Similar to many class 2 effectors—and unique among class 1 systems—DiCas7-11 processes pre-CRISPR RNA into mature CRISPR RNA (crRNA) and cleaves RNA at positions defined by the target:spacer duplex, without detectable non-specific activity. We engineered Cas7-11 for RNA knockdown and editing in mammalian cells. We show that Cas7-11 has no effects on cell viability, whereas other RNA-targeting tools (such as short hairpin RNAs and Cas13) show substantial cell toxicity4,5. This study illustrates the evolution of a single-protein effector from multisubunit class 1 effector complexes, expanding our understanding of the diversity of CRISPR systems. Cas7-11 provides the basis for new programmable RNA-targeting tools that are free of collateral activity and cell toxicity.

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.

CrAss-like phages are a recently described expansive group of viruses that includes the most abundant virus in the human gut1–3. The genomes of all crAss-like phages encode a large virion-packaged protein2,4that contains a DFDxD sequence motif, which forms the catalytic site in cellular multisubunit RNA polymerases (RNAPs)5. Here, using Cellulophaga balticacrAss-like phage phi14:2 as a model system, we show that this protein is a DNA-dependent RNAP that is translocated into the host cell along with the phage DNA and transcribes early phage genes. We determined the crystal structure of this 2,180-residue enzyme in a self-inhibited state, which probably occurs before virion packaging. This conformation is attained with the help of a cleft-blocking domain that interacts with the active site and occupies the cavity in which the RNA–DNA hybrid binds. Structurally, phi14:2 RNAP is most similar to eukaryotic RNAPs that are involved in RNA interference6,7, although most of the phi14:2 RNAP structure (nearly 1,600 residues) maps to a new region of the protein fold space. Considering this structural similarity, we propose that eukaryal RNA interference polymerases have their origins in phage, which parallels the emergence of the mitochondrial transcription apparatus8.

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.

The last universal cellular ancestor (LUCA) is the most recent population of organisms from which all cellular life on Earth descends. The reconstruction of the genome and phenotype of the LUCA is a major challenge in evolutionary biology. Given that all life forms are associated with viruses and/or other mobile genetic elements, there is no doubt that the LUCA was a host to viruses. Here, by projecting back in time using the extant distribution of viruses across the two primary domains of life, bacteria and archaea, and tracing the evolutionary histories of some key virus genes, we attempt a reconstruction of the LUCA virome. Even a conservative version of this reconstruction suggests a remarkably complex virome that already included the main groups of extant viruses of bacteria and archaea. We further present evidence of extensive virus evolution antedating the LUCA. The presence of a highly complex virome implies the substantial genomic and pan-genomic complexity of the LUCA itself.

The number and diversity of known CRISPR–Cas systems have substantially increased in recent years. Here, we provide an updated evolutionary classification of CRISPR–Cas systems and casgenes, with an emphasis on the major developments that have occurred since the publication of the latest classification, in 2015. The new classification includes 2 classes, 6 types and 33 subtypes, compared with 5 types and 16 subtypes in 2015. A key development is the ongoing discovery of multiple, novel class 2 CRISPR–Cas systems, which now include 3 types and 17 subtypes. A second major novelty is the discovery of numerous derived CRISPR–Cas variants, often associated with mobile genetic elements that lack the nucleases required for interference. Some of these variants are involved in RNA-guided transposition, whereas others are predicted to perform functions distinct from adaptive immunity that remain to be characterized experimentally. The third highlight is the discovery of numerous families of ancillary CRISPR-linked genes, often implicated in signal transduction. Together, these findings substantially clarify the functional diversity and evolutionary history of CRISPR–Cas.

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.

The signature component of type III CRISPR-Cas systems is the Cas10 protein that consists of two Palm domains homologous to those of DNA and RNA polymerases and nucleotide cyclases and an HD nuclease domain. However, until very recently, the activity of the Palm domains and their role in CRISPR function have not been experimentally established. Most of the type III CRISPR-Cas systems and some type I systems also encompass proteins containing the CARF (CRISPR-associated Rossmann fold) domain that has been predicted to regulate CRISPR functions vianucleotide binding, but its function in CRISPR-Cas remained obscure. Two independent recent studies show that the Palm domain of Cas10 catalyzes synthesis of oligoadenylates, which bind the CARF domain of the Csm6 protein and activate its RNase domain that cleaves foreign transcripts enabling interference by type III CRISPR-Cas. In one coup, these findings resolved two long-standing puzzles of CRISPR biology and reveal a new regulatory pathway that governs the CRISPR response. However, the full extent of this pathway, and especially the driving forces behind the evolution of this complex mechanism of CRISPR-Cas activation, remains to be uncovered.

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.

Functionally linked genes in bacterial and archaeal genomes are often organized into operons. However, the composition and architecture of operons are highly variable and frequently differ even among closely related genomes. Therefore, to efficiently extract reliable functional predictions for uncharacterized genes from comparative analyses of the rapidly growing genomic databases, dedicated computational approaches are required. We developed a protocol to systematically and automatically identify genes that are likely to be functionally associated with a ‘bait’ gene or locus by using relevance metrics. Given a set of bait loci and a genomic database defined by the user, this protocol compares the genomic neighborhoods of the baits to identify genes that are likely to be functionally linked to the baits by calculating the abundance of a given gene within and outside the bait neighborhoods and the distance to the bait. We exemplify the performance of the protocol with three test cases, namely, genes linked to CRISPR–Cas systems using the ‘CRISPRicity’ metric, genes associated with archaeal proviruses and genes linked to Argonaute genes in halobacteria. The protocol can be run by users with basic computational skills. The computational cost depends on the sizes of the genomic dataset and the list of reference loci and can vary from one CPU-hour to hundreds of hours on a supercomputer.

ABSTRACTTrans-activating CRISPR (tracr) RNA is a distinct RNA species that interacts with the CRISPR (cr) RNA to form the dual guide (g) RNA in type II and subtype V-B CRISPR-Cas systems. The tracrRNA-crRNA interaction is essential for pre-crRNA processing as well as target recognition and cleavage. The tracrRNA consists of an antirepeat, which forms an imperfect hybrid with the repeat in the crRNA, and a distal region containing a Rho-independent terminator. Exhaustive comparative analysis of the sequences and predicted structures of the Class 2 CRISPR guide RNAs shows that all these guide RNAs share distinct structural features, in particular, the nexus stem-loop that separates the repeat-antirepeat hybrid from the distal portion of the tracrRNA and the conserved GU pair at that end of the hybrid. These structural constraints might ensure full exposure of the spacer for target recognition. Reconstruction of tracrRNA evolution for 4 tight bacterial groups demonstrates random drift of repeat-antirepeat complementarity within a window of hybrid stability that is, apparently, maintained by selection. An evolutionary scenario is proposed whereby tracrRNAs evolved on multiple occasions, viarearrangement of a CRISPR array to form the antirepeat in different locations with respect to the array. A functional tracrRNA would form if, in the new location, the antirepeat is flanked by sequences that meet the minimal requirements for a promoter and a Rho-independent terminator. Alternatively, or additionally, the antirepeat sequence could be occasionally ‘reset’ by recombination with a repeat, restoring the functionality of tracrRNAs that drift beyond the required minimal hybrid stability.

One of the most prominent features of archaea is the extraordinary diversity of their DNA viruses. Many archaeal viruses differ substantially in morphology from bacterial and eukaryotic viruses and represent unique virus families. The distinct nature of archaeal viruses also extends to the gene composition and architectures of their genomes and the properties of the proteins that they encode. Environmental research has revealed prominent roles of archaeal viruses in influencing microbial communities in ocean ecosystems, and recent metagenomic studies have uncovered new groups of archaeal viruses that infect extremophiles and mesophiles in diverse habitats. In this Review, we summarize recent advances in our understanding of the genomic and morphological diversity of archaeal viruses and the molecular biology of their life cycles and virus–host interactions, including interactions with archaeal CRISPR–Cas systems. We also examine the potential origins and evolution of archaeal viruses and discuss their place in the global virosphere.

Class 2 CRISPR–Cas systems are characterized by effector modules that consist of a single multidomain protein, such as Cas9 or Cpf1. We designed a computational pipeline for the discovery of novel class 2 variants and used it to identify six new CRISPR–Cas subtypes. The diverse properties of these new systems provide potential for the development of versatile tools for genome editing and regulation. In this Analysis article, we present a comprehensive census of class 2 types and class 2 subtypes in complete and draft bacterial and archaeal genomes, outline evolutionary scenarios for the independent origin of different class 2 CRISPR–Cas systems from mobile genetic elements, and propose an amended classification and nomenclature of CRISPR–Cas.

The number and diversity of viral sequences that are identified in metagenomic data far exceeds that of experimentally characterized virus isolates. In a recent workshop, a panel of experts discussed the proposal that, with appropriate quality control, viruses that are known only from metagenomic data can, and should be, incorporated into the official classification scheme of the International Committee on Taxonomy of Viruses (ICTV). Although a taxonomy that is based on metagenomic sequence data alone represents a substantial departure from the traditional reliance on phenotypic properties, the development of a robust framework for sequence-based virus taxonomy is indispensable for the comprehensive characterization of the global virome. In this Consensus Statement article, we consider the rationale for why metagenomic sequence data should, and how it can, be incorporated into the ICTV taxonomy, and present proposals that have been endorsed by the Executive Committee of the ICTV.

The study of genes that evolve under conditional selection can shed light on the genomic underpinnings of adaptation, revealing epistasis and phenotypic plasticity. This protocol describes how to use the Coselenspackage to compare gene-level selection between two groups of samples. After installing Coselensand preparing the datasets, a typical run on a laptop takes less than 10 min. Coselensis best suited to analyze somatic mutations and data from experimental evolution, for which independently evolved samples are available.

Viruses hold solutions to a lot of problems, so let’s fund and reward cataloguing, urge Jens H. Kuhn and colleagues.

Cancer driver mutations often display mutual exclusion or co-occurrence, underscoring the key role of epistasis in carcinogenesis. However, estimating the magnitude of epistasis and quantifying its effect on tumor evolution remains a challenge. We develop a method (Coselens) to quantify conditional selection on the excess of nonsynonymous substitutions in cancer genes. Coselensinfers the number of drivers per gene in different partitions of a cancer genomics dataset using covariance-based mutation models and determines whether coding mutations in a gene affect selection for drivers in any other gene. Using Coselens, we identify 296 conditionally selected gene pairs across 16 cancer types in the TCGA dataset. Conditional selection affects 25%–50% of driver substitutions in tumors with >2 drivers. Conditionally co-selected genes form modular networks, whose structures challenge the traditional interpretation of within-pathway mutual exclusivity and across-pathway synergy, suggesting a more complex scenario where gene-specific across-pathway epistasis shapes differentiated cancer subtypes.