Synergistic epistasis among cancer drivers can rescue early tumors from the accumulation of deleterious passengers.

Epistasis among driver mutations is pervasive and explains relevant features of cancer, such as differential therapy response and convergence towards well-characterized molecular subtypes. Furthermore, a growing body of evidence suggests that tumor development could be hampered by the accumulation of slightly deleterious passenger mutations. In this work, we combined empirical epistasis networks, computer simulations, and mathematical models to explore how synergistic interactions among driver mutations affect cancer progression under the burden of slightly deleterious passengers. We found that epistasis plays a crucial role in tumor development by promoting the transformation of precancerous clones into rapidly growing tumors through a process that is analogous to evolutionary rescue. The triggering of epistasis-driven rescue is strongly dependent on the intensity of epistasis and could be a key rate-limiting step in many tumors, contributing to their unpredictability. As a result, central genes in cancer epistasis networks appear as key intervention targets for cancer therapy.

UG/Abi: a highly diverse family of prokaryotic reverse transcriptases associated with defense functions.

Reverse transcriptases (RTs) are enzymes capable of synthesizing DNA using RNA as a template. Within the last few years, a burst of research has led to the discovery of novel prokaryotic RTs with diverse antiviral properties, such as DRTs (Defense-associated RTs), which belong to the so-called group of unknown RTs (UG) and are closely related to the Abortive Infection system (Abi) RTs. In this work, we performed a systematic analysis of UG and Abi RTs, increasing the number of UG/Abi members up to 42 highly diverse groups, most of which are predicted to be functionally associated with other gene(s) or domain(s). Based on this information, we classified these systems into three major classes. In addition, we reveal that most of these groups are associated with defense functions and/or mobile genetic elements, and demonstrate the antiphage role of four novel groups. Besides, we highlight the presence of one of these systems in novel families of human gut viruses infecting members of the Bacteroidetes and Firmicutes phyla. This work lays the foundation for a comprehensive and unified understanding of these highly diverse RTs with enormous biotechnological potential.

High density of unrepaired genomic ribonucleotides leads to Topoisomerase 1-mediated severe growth defects in absence of ribonucleotide reductase.

Cellular levels of ribonucleoside triphosphates (rNTPs) are much higher than those of deoxyribonucleoside triphosphates (dNTPs), thereby influencing the frequency of incorporation of ribonucleoside monophosphates (rNMPs) by DNA polymerases (Pol) into DNA. RNase H2-initiated ribonucleotide excision repair (RER) efficiently removes single rNMPs in genomic DNA. However, processing of rNMPs by Topoisomerase 1 (Top1) in absence of RER induces mutations and genome instability. Here, we greatly increased the abundance of genomic rNMPs in Saccharomyces cerevisiae by depleting Rnr1, the major subunit of ribonucleotide reductase, which converts ribonucleotides to deoxyribonucleotides. We found that in strains that are depleted of Rnr1, RER-deficient, and harbor an rNTP-permissive replicative Pol mutant, excessive accumulation of single genomic rNMPs severely compromised growth, but this was reversed in absence of Top1. Thus, under Rnr1 depletion, limited dNTP pools slow DNA synthesis by replicative Pols and provoke the incorporation of high levels of rNMPs in genomic DNA. If a threshold of single genomic rNMPs is exceeded in absence of RER and presence of limited dNTP pools, Top1-mediated genome instability leads to severe growth defects. Finally, we provide evidence showing that accumulation of RNA/DNA hybrids in absence of RNase H1 and RNase H2 leads to cell lethality under Rnr1 depletion.

Cancer-mutation network and the number and specificity of driver mutations.

Cancer genomics has produced extensive information on cancer-associated genes, but the number and specificity of cancer-driver mutations remains a matter of debate. We constructed a bipartite network in which 7,665 tumors from 30 cancer types are connected via shared mutations in 198 previously identified cancer genes. We show that about 27% of the tumors can be assigned to statistically supported modules, most of which encompass one or two cancer types. The rest of the tumors belong to a diffuse network component suggesting lower gene specificity of driver mutations. Linear regression of the mutational loads in cancer genes was used to estimate the number of drivers required for the onset of different cancers. The mean number of drivers in known cancer genes is approximately two, with a range of one to five. Cancers that are associated with modules had more drivers than those from the diffuse network component, suggesting that unidentified and/or interchangeable drivers exist in the latter.

Disentangling the effects of selection and loss bias on gene dynamics.

We combine mathematical modeling of genome evolution with comparative analysis of prokaryotic genomes to estimate the relative contributions of selection and intrinsic loss bias to the evolution of different functional classes of genes and mobile genetic elements (MGE). An exact solution for the dynamics of gene family size was obtained under a linear duplication-transfer-loss model with selection. With the exception of genes involved in information processing, particularly translation, which are maintained by strong selection, the average selection coefficient for most nonparasitic genes is low albeit positive, compatible with observed positive correlation between genome size and effective population size. Free-living microbes evolve under stronger selection for gene retention than parasites. Different classes of MGE show a broad range of fitness effects, from the nearly neutral transposons to prophages, which are actively eliminated by selection. Genes involved in antiparasite defense, on average, incur a fitness cost to the host that is at least as high as the cost of plasmids. This cost is probably due to the adverse effects of autoimmunity and curtailment of horizontal gene transfer caused by the defense systems and selfish behavior of some of these systems, such as toxin-antitoxin and restriction modification modules. Transposons follow a biphasic dynamics, with bursts of gene proliferation followed by decay in the copy number that is quantitatively captured by the model. The horizontal gene transfer to loss ratio, but not duplication to loss ratio, correlates with genome size, potentially explaining increased abundance of neutral and costly elements in larger genomes.

Myosin-driven transport network in plants.

We investigate the myosin XI-driven transport network in Arabidopsis using protein-protein interaction, subcellular localization, gene knockout, and bioinformatics analyses. The two major groups of nodes in this network are myosins XI and their membrane-anchored receptors (MyoB) that, together, drive endomembrane trafficking and cytoplasmic streaming in the plant cells. The network shows high node connectivity and is dominated by generalists, with a smaller fraction of more specialized myosins and receptors. We show that interaction with myosins and association with motile vesicles are common properties of the MyoB family receptors. We identify previously uncharacterized myosin-binding proteins, putative myosin adaptors that belong to two unrelated families, with four members each (MadA and MadB). Surprisingly, MadA1 localizes to the nucleus and is rapidly transported to the cytoplasm, suggesting the existence of myosin XI-driven nucleocytoplasmic trafficking. In contrast, MadA2 and MadA3, as well as MadB1, partition between the cytosolic pools of motile endomembrane vesicles that colocalize with myosin XI-K and diffuse material that does not. Gene knockout analysis shows that MadB1-4 contribute to polarized root hair growth, phenocopying myosins, whereas MadA1-4 are redundant for this process. Phylogenetic analysis reveals congruent evolutionary histories of the myosin XI, MyoB, MadA, and MadB families. All these gene families emerged in green algae and show concurrent expansions via serial duplication in flowering plants. Thus, the myosin XI transport network increased in complexity and robustness concomitantly with the land colonization by flowering plants and, by inference, could have been a major contributor to this process.

Bipartite Network Analysis of the Archaeal Virosphere: Evolutionary Connections between Viruses and Capsidless Mobile Elements.

Archaea and particularly hyperthermophilic crenarchaea are hosts to many unusual viruses with diverse virion shapes and distinct gene compositions. As is typical of viruses in general, there are no universal genes in the archaeal virosphere. Therefore, to obtain a comprehensive picture of the evolutionary relationships between viruses, network analysis methods are more productive than traditional phylogenetic approaches. Here we present a comprehensive comparative analysis of genomes and proteomes from all currently known taxonomically classified and unclassified, cultivated and uncultivated archaeal viruses. We constructed a bipartite network of archaeal viruses that includes two classes of nodes, the genomes and gene families that connect them. Dissection of this network using formal community detection methods reveals strong modularity, with 10 distinct modules and 3 putative supermodules. However, compared to similar previously analyzed networks of eukaryotic and bacterial viruses, the archaeal virus network is sparsely connected. With the exception of the tailed viruses related to bacteriophages of the order Caudovirales and the families Turriviridae and Sphaerolipoviridae that are linked to a distinct supermodule of eukaryotic and bacterial viruses, there are few connector genes shared by different archaeal virus modules. In contrast, most of these modules include, in addition to viruses, capsidless mobile elements, emphasizing tight evolutionary connections between the two types of entities in archaea. The relative contributions of distinct evolutionary origins, in particular from nonviral elements, and insufficient sampling to the sparsity of the archaeal virus network remain to be determined by further exploration of the archaeal virosphere. IMPORTANCE: Viruses infecting archaea are among the most mysterious denizens of the virosphere. Many of these viruses display no genetic or even morphological relationship to viruses of bacteria and eukaryotes, raising questions regarding their origins and position in the global virosphere. Analysis of 5,740 protein sequences from 116 genomes allowed dissection of the archaeal virus network and showed that most groups of archaeal viruses are evolutionarily connected to capsidless mobile genetic elements, including various plasmids and transposons. This finding could reflect actual independent origins of the distinct groups of archaeal viruses from different nonviral elements, providing important insights into the emergence and evolution of the archaeal virome.

Immunity, suicide or both? Ecological determinants for the combined evolution of anti-pathogen defense systems.

BACKGROUND: Parasite-host arms race is one of the key factors in the evolution of life. Most cellular life forms, in particular prokaryotes, possess diverse forms of defense against pathogens including innate immunity, adaptive immunity and programmed cell death (altruistic suicide). Coevolution of these different but interacting defense strategies yields complex evolutionary regimes. RESULTS: We develop and extensively analyze a computational model of coevolution of different defense strategies to show that suicide as a defense mechanism can evolve only in structured populations and when the attainable degree of immunity against pathogens is limited. The general principle of defense evolution seems to be that hosts do not evolve two costly defense mechanisms when one is sufficient. Thus, the evolutionary interplay of innate immunity, adaptive immunity and suicide, leads to an equilibrium state where the combination of all three defense strategies is limited to a distinct, small region of the parameter space. The three strategies can stably coexist only if none of them are highly effective. Coupled adaptive immunity-suicide systems, the existence of which is implied by the colocalization of genes for the two types of defense in prokaryotic genomes, can evolve either when immunity-associated suicide is more efficacious than other suicide systems or when adaptive immunity functionally depends on the associated suicide system. CONCLUSIONS: Computational modeling reveals a broad range of outcomes of coevolution of anti-pathogen defense strategies depending on the relative efficacy of different mechanisms and population structure. Some of the predictions of the model appear compatible with recent experimental evolution results and call for additional experiments.

Large-scale genomic analysis suggests a neutral punctuated dynamics of transposable elements in bacterial genomes.

Insertion sequences (IS) are the simplest and most abundant form of transposable DNA found in bacterial genomes. When present in multiple copies, it is thought that they can promote genomic plasticity and genetic exchange, thus being a major force of evolutionary change. The main processes that determine IS content in genomes are, though, a matter of debate. In this work, we take advantage of the large amount of genomic data currently available and study the abundance distributions of 33 IS families in 1811 bacterial chromosomes. This allows us to test simple models of IS dynamics and estimate their key parameters by means of a maximum likelihood approach. We evaluate the roles played by duplication, lateral gene transfer, deletion and purifying selection. We find that the observed IS abundances are compatible with a neutral scenario where IS proliferation is controlled by deletions instead of purifying selection. Even if there may be some cases driven by selection, neutral behavior dominates over large evolutionary scales. According to this view, IS and hosts tend to coexist in a dynamic equilibrium state for most of the time. Our approach also allows for a detection of recent IS expansions, and supports the hypothesis that rapid expansions constitute transient events-punctuations-during which the state of coexistence of IS and host becomes perturbated.

Tempo and mode of inhibitor-mutagen antiviral therapies: a multidisciplinary approach.

The continuous emergence of drug-resistant viruses is a major obstacle for the successful treatment of viral infections, thus representing a persistent spur to the search for new therapeutic strategies. Among them, multidrug treatments are currently at the forefront of pharmaceutical, clinical, and computational investigation. Still, there are many unknowns in the way that different drugs interact among themselves and with the pathogen that they aim to control. Inspired by experimental studies with picornavirus, here, we discuss the performance of sequential vs. combination therapies involving two dissimilar drugs: the mutagen ribavirin and an inhibitor of viral replication, guanidine. Because a systematic analysis of viral response to drug doses demands a precious amount of time and resources, we present and analyze an in silico model describing the dynamics of the viral population under the action of the two drugs. The model predicts the response of the viral population to any dose combination, the optimal therapy to be used in each case, and the way to minimize the probability of appearance of resistant mutants. In agreement with the theoretical predictions, in vitro experiments with foot-and-mouth disease virus confirm that the suitability of simultaneous or sequential administration depends on the drug doses. In addition, intrinsic replicative characteristics of the virus (e.g., replication through RNA only or a DNA intermediate) play a key role to determine the appropriateness of a sequential or combination therapy. Knowledge of several model parameters can be derived by means of few, simple experiments, such that the model and its predictions can be extended to other viral systems.

Dynamics of CRISPR-mediated virus-host interactions in the human gut microbiome.

Arms races between mobile genetic elements and prokaryotic hosts are major drivers of ecological and evolutionary change in microbial communities. Prokaryotic defense systems such as CRISPR-Cas have the potential to regulate microbiome composition by modifying the interactions among bacteria, plasmids, and phages. Here, we used longitudinal metagenomic data from 130 healthy and diseased individuals to study how the interplay of genetic parasites and CRISPR-Cas immunity reflects on the dynamics and composition of the human gut microbiome. Based on the coordinated study of 80 000 CRISPR-Cas loci and their targets, we show that CRISPR-Cas immunity effectively modulates bacteriophage abundances in the gut. Acquisition of CRISPR-Cas immunity typically leads to a decrease in the abundance of lytic phages but does not necessarily cause their complete disappearance. Much smaller effects are observed for lysogenic phages and plasmids. Conversely, phage-CRISPR interactions shape bacterial microdiversity by producing weak selective sweeps that benefit immune host lineages. We also show that distal (and chronologically older) regions of CRISPR arrays are enriched in spacers that are potentially functional and target crass-like phages and local prophages. This suggests that exposure to reactivated prophages and other endemic viruses is a major selective pressure in the gut microbiome that drives the maintenance of long-lasting immune memory.

Evolutionary dynamics of the prokaryotic adaptive immunity system CRISPR-Cas in an explicit ecological context.

A stochastic, agent-based mathematical model of the coevolution of the archaeal and bacterial adaptive immunity system, CRISPR-Cas, and lytic viruses shows that CRISPR-Cas immunity can stabilize the virus-host coexistence rather than leading to the extinction of the virus. In the model, CRISPR-Cas immunity does not specifically promote viral diversity, presumably because the selection pressure on each single proto-spacer is too weak. However, the overall virus diversity in the presence of CRISPR-Cas grows due to the increase of the host and, accordingly, the virus population size. Above a threshold value of total viral diversity, which is proportional to the viral mutation rate and population size, the CRISPR-Cas system becomes ineffective and is lost due to the associated fitness cost. Our previous modeling study has suggested that the ubiquity of CRISPR-Cas in hyperthermophiles, which contrasts its comparative low prevalence in mesophiles, is due to lower rates of mutation fixation in thermal habitats. The present findings offer a complementary, simpler perspective on this contrast through the larger population sizes of mesophiles compared to hyperthermophiles, because of which CRISPR-Cas can become ineffective in mesophiles. The efficacy of CRISPR-Cas sharply increases with the number of proto-spacers per viral genome, potentially explaining the low information content of the proto-spacer-associated motif (PAM) that is required for spacer acquisition by CRISPR-Cas because a higher specificity would restrict the number of spacers available to CRISPR-Cas, thus hampering immunity. The very existence of the PAM might reflect the tradeoff between the requirement of diverse spacers for efficient immunity and avoidance of autoimmunity.

Protocol for comparing gene-level selection on coding mutations between two groups of samples with Coselens.

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 Coselens package to compare gene-level selection between two groups of samples. After installing Coselens and preparing the datasets, a typical run on a laptop takes less than 10 min. Coselens is best suited to analyze somatic mutations and data from experimental evolution, for which independently evolved samples are available. For complete details on the use and execution of this protocol, please refer to Iranzo et al. (2022).1.

Comparison of gene clustering criteria reveals intrinsic uncertainty in pangenome analyses.

BACKGROUND: A key step for comparative genomics is to group open reading frames into functionally and evolutionarily meaningful gene clusters. Gene clustering is complicated by intraspecific duplications and horizontal gene transfers that are frequent in prokaryotes. In consequence, gene clustering methods must deal with a trade-off between identifying vertically transmitted representatives of multicopy gene families, which are recognizable by synteny conservation, and retrieving complete sets of species-level orthologs. We studied the implications of adopting homology, orthology, or synteny conservation as formal criteria for gene clustering by performing comparative analyses of 125 prokaryotic pangenomes. RESULTS: Clustering criteria affect pangenome functional characterization, core genome inference, and reconstruction of ancestral gene content to different extents. Species-wise estimates of pangenome and core genome sizes change by the same factor when using different clustering criteria, allowing robust cross-species comparisons regardless of the clustering criterion. However, cross-species comparisons of genome plasticity and functional profiles are substantially affected by inconsistencies among clustering criteria. Such inconsistencies are driven not only by mobile genetic elements, but also by genes involved in defense, secondary metabolism, and other accessory functions. In some pangenome features, the variability attributed to methodological inconsistencies can even exceed the effect sizes of ecological and phylogenetic variables. CONCLUSIONS: Choosing an appropriate criterion for gene clustering is critical to conduct unbiased pangenome analyses. We provide practical guidelines to choose the right method depending on the research goals and the quality of genome assemblies, and a benchmarking dataset to assess the robustness and reproducibility of future comparative studies.

Evolutionary dynamics of genome segmentation in multipartite viruses.

Multipartite viruses are formed by a variable number of genomic fragments packed in independent viral capsids. This fact poses stringent conditions on their transmission mode, demanding, in particular, a high multiplicity of infection (MOI) for successful propagation. The actual advantages of the multipartite viral strategy are as yet unclear. The origin of multipartite viruses represents an evolutionary puzzle. While classical theories suggested that a faster replication rate or higher replication fidelity would favour shorter segments, recent experimental results seem to point to an increased stability of virions with incomplete genomes as a factor able to compensate for the disadvantage of mandatory complementation. Using as main parameters differential stability as a function of genome length and MOI, we calculate the conditions under which a set of complementary segments of a viral genome would outcompete the non-segmented variant. Further, we examine the likeliness that multipartite viral forms could be the evolutionary outcome of the competition among the defective genomes of different lengths that spontaneously arise under replication of a complete, wild-type genome. We conclude that only multipartite viruses with a small number of segments could be produced in our scenario, and discuss alternative hypotheses for the origin of multipartite viruses with more than four segments.

Pervasive conditional selection of driver mutations and modular epistasis networks in cancer.

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. Coselens infers 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.

The spatial Ultimatum game revisited.

We revisit the issue of the emergence of fair behavior in the framework of the spatial Ultimatum game, adding many important results and insights to the pioneering work by Page et al. [2000. The spatial Ultimatum game. Proc. R. Soc. London B 267, 2177], who showed in a specific example that on a two-dimensional setup evolution may lead to strategies with some degree of fairness. Within this spatial framework, we carry out a thorough simulation study and show that the emergence of altruism is a very generic phenomenon whose details depend on the dynamics considered. A very frequent feature is the spontaneous emergence and fixation of quasiempathetic individuals, whose offers are very close to their acceptance thresholds. We present analytical arguments that allow an understanding of our results and give insights on the manner in which local effects in evolution may lead to such non-rational or apparently maladaptive behaviors.

Game-Theoretical Modeling of Interviral Conflicts Mediated by Mini-CRISPR Arrays.

All cellular organisms coevolve with multiple viruses, so that both virus-host and intervirus conflicts are major factors of evolution. Accordingly, hosts evolve multiple, elaborate defense systems and viruses respond by evolving means of antidefense. Although less thoroughly characterized, several dedicated mechanisms of intervirus competition have been described as well. Recently, the genomes of some bacterial and archaeal viruses have been shown to harbor CRISPR mini-arrays that typically contain a single spacer targeting a closely related virus. The involvement of mini-arrays in an intervirus conflict has been experimentally demonstrated for a pair of archaeal viruses. We model the evolution of virus-encoded CRISPR mini-arrays using a game theoretical approach. Analysis of the model reveals multiple equilibria that include mutual targeting, unidirectional targeting, no targeting, cyclic polymorphism, and bistability. The choice between these evolutionary regimes depends on the model parameters including the coinfection frequency, differential productivity of the conflicting viruses, and the fitness cost of mini-arrays. At high coinfection frequencies, the model becomes a version of the Prisoner's dilemma in which defection, i.e., mutual targeting between the competing viruses, is the winning strategy.

Genomes of the "Candidatus Actinomarinales" Order: Highly Streamlined Marine Epipelagic Actinobacteria.

"Candidatus Actinomarinales" was defined as a subclass of exclusively marine Actinobacteria with small cells and genomes. We have collected all the available genomes in databases to assess the diversity included in this group and analyzed it by comparative genomics. We have found the equivalent of five genera and 18 genomospecies. They have genome reduction parameters equal to those of freshwater actinobacterial "Candidatus Nanopelagicales" or marine alphaproteobacterial Pelagibacterales Genome recruitment shows that they are found only in the photic zone and mainly in surface waters, with only one genus that is found preferentially at or below the deep chlorophyll maximum. "Ca Actinomarinales" show a highly conserved core genome (80% of the gene families conserved for the whole order) with a saturation of genomic diversity of the flexible genome at the genomospecies level. We found only a flexible genomic island preserved throughout the order; it is related to the sugar decoration of the envelope and uses several tRNAs as hot spots to increase its genomic diversity. Populations had a discrete level of sequence diversity similar to other marine microbes but drastically different from the much higher levels found for Pelagibacterales Genomic analysis suggests that they are all aerobic photoheterotrophs with one type 1 rhodopsin and a heliorhodopsin. Like other actinobacteria, they possess the F420 coenzyme biosynthesis pathway, and its lower reduction potential could provide access to an increased range of redox chemical transformations. Last, sequence analysis revealed the first "Ca Actinomarinales" phages, including a prophage, with metaviromic islands related to sialic acid cleavage.IMPORTANCE Microbiology is in a new age in which sequence databases are primary sources of information about many microbes. However, in-depth analysis of environmental genomes thus retrieved is essential to substantiate the new knowledge. Here, we study 182 genomes belonging to the only known exclusively marine pelagic group of the phylum Actinobacteria The aquatic branch of this phylum is largely known from environmental sequencing studies (single-amplified genomes [SAGs] and metagenome-assembled genomes [MAGs]), and we have collected and analyzed the available information present in databases about the "Ca. Actinomarinales." They are among the most streamlined microbes to live in the epipelagic zone of the ocean, and their study is critical to obtain a proper view of the diversity of Actinobacteria and their role in aquatic ecosystems.

Evolutionary classification of CRISPR-Cas systems: a burst of class 2 and derived variants.

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 cas genes, 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.