Metagenomic Discovery of CRISPR-Associated Transposons

CRISPR-associated transposons (CASTs) co-opt Cas genes for RNA-guided transposition. CASTs are exceedingly rare in genomic databases; recent surveys have reported Tn7-like transposons that co-opt Type I-F, I-B, and V-K CRISPR effectors. Here, we expand the diversity of reported CAST systems via a bioinformatic search of metagenomic databases. We discover new architectures for all known CASTs, including novel arrangements of the Cascade effectors, new self-targeting modalities, and minimal V-K systems. We also describe new families of CASTs that have co-opted the Type I-C and Type IV CRISPR-Cas systems. Our search for non-Tn7 CASTs identifies putative candidates that co-opt Cas12a for horizontal gene transfer. These new systems shed light on how CRISPR systems have co-evolved with transposases and expand the programmable gene editing toolkit.

Draft Genome Sequence of Enterococcus faecalis AS003, a Strain Possessing All Three Type II-a CRISPR Loci

ABSTRACT Enterococcus faecalis is a clinically significant member of the human microbiome. Three CRISPR-Cas loci are located in conserved locations. Previous studies provide evidence that E. faecalis strains with functional CRISPR-Cas genes are negatively correlated with antibiotic resistance. Here, we report the genome sequence of an unusual strain possessing all three CRISPR-Cas loci.

The phylogenomics of CRISPR-Cas system and revelation of its features in Salmonella

Salmonellae display intricate evolutionary patterns comprising over 2500 serovars having diverse pathogenic profiles. The acquisition and/or exchange of various virulence factors influences the evolutionary framework. To gain insights into evolution of Salmonella in association with the CRISPR-Cas genes we performed phylogenetic surveillance across strains of 22 Salmonella serovars. The strains differed in their CRISPR1-leader and cas operon features assorting into two main clades, CRISPR1-STY/cas-STY and CRISPR1-STM/cas-STM, comprising majorly typhoidal and non-typhoidal Salmonella serovars respectively. Serovars of these two clades displayed better relatedness, concerning CRISPR1-leader and cas operon, across genera than between themselves. This signifies the acquisition of CRISPR1/Cas region could be through a horizontal gene transfer event owing to the presence of mobile genetic elements flanking CRISPR1 array. Comparison of CRISPR and cas phenograms with that of multilocus sequence typing (MLST) suggests differential evolution of CRISPR/Cas system. As opposed to broad-host-range, the host-specific serovars harbor fewer spacers. Mapping of protospacer sources suggested a partial correlation of spacer content with habitat diversity of the serovars. Some serovars like serovar Enteritidis and Typhimurium that inhabit similar environment/infect similar hosts hardly shared their protospacer sources.

Casboundary: automated definition of integral Cas cassettes

Motivation CRISPR-Cas are important systems found in most archaeal and many bacterial genomes, providing adaptive immunity against mobile genetic elements in prokaryotes. The CRISPR-Cas systems are encoded by a set of consecutive cas genes, here termed cassette. The identification of cassette boundaries is key for finding cassettes in CRISPR research field. This is often carried out by using Hidden Markov Models and manual annotation. In this article, we propose the first method able to automatically define the cassette boundaries. In addition, we present a Cas-type predictive model used by the method to assign each gene located in the region defined by a cassette’s boundaries a Cas label from a set of pre-defined Cas types. Furthermore, the proposed method can detect potentially new cas genes and decompose a cassette into its modules. Results We evaluate the predictive performance of our proposed method on data collected from the two most recent CRISPR classification studies. In our experiments, we obtain an average similarity of 0.86 between the predicted and expected cassettes. Besides, we achieve F-scores above 0.9 for the classification of cas genes of known types and 0.73 for the unknown ones. Finally, we conduct two additional study cases, where we investigate the occurrence of potentially new cas genes and the occurrence of module exchange between different genomes. Availability and implementation https://github.com/BackofenLab/Casboundary. Contact [email protected] or [email protected] Supplementary information Supplementary data are available at Bioinformatics online.

CasCollect: targeted assembly of CRISPR-associated operons from high-throughput sequencing data

CRISPR arrays and CRISPR-associated (Cas) proteins comprise a widespread adaptive immune system in bacteria and archaea. These systems function as a defense against exogenous parasitic mobile genetic elements that include bacteriophages, plasmids and foreign nucleic acids. With the continuous spread of antibiotic resistance, knowledge of pathogen susceptibility to bacteriophage therapy is becoming more critical. Additionally, gene-editing applications would benefit from the discovery of new cas genes with favorable properties. While next-generation sequencing has produced staggering quantities of data, transitioning from raw sequencing reads to the identification of CRISPR/Cas systems has remained challenging. This is especially true for metagenomic data, which has the highest potential for identifying novel cas genes. We report a comprehensive computational pipeline, CasCollect, for the targeted assembly and annotation of cas genes and CRISPR arrays—even isolated arrays—from raw sequencing reads. Benchmarking our targeted assembly pipeline demonstrates significantly improved timing by almost two orders of magnitude compared with conventional assembly and annotation, while retaining the ability to detect CRISPR arrays and cas genes. CasCollect is a highly versatile pipeline and can be used for targeted assembly of any specialty gene set, reconfigurable for user provided Hidden Markov Models and/or reference nucleotide sequences.

The Phylogenomics of CRISPR-Cas System in Salmonella: an evolutionary race with housekeeping genes

Salmonellae display intricate evolutionary patterns comprising over 2500 serovars having diverse pathogenic profiles. The acquisition/exchange of various virulence factors influence the evolutionary framework. To gain insights into evolution of Salmonella as a pathogen in association with the CRISPR-Cas genes we performed phylogenetic surveillance across strains of 22 Salmonella serovars. The strains assorted into two main clades, pertaining to the differences in their CRISPR1-leader and cas operon. Considering Salmonella enterica subsp. enterica serovar Typhimurium and serovar Typhi as signature serovars, we classified the clades as CRISPR1-STM/cas-STM and CRISPR1-STY/cas-STY, respectively. Serovars of the two clades displayed better relatedness, concerning CRISPR-1 leader and cas operon, across genera than between themselves. This signifies the acquisition of CRISPR1/Cas region a horizontal gene transfer event owing to the presence of mobile genetic elements flanking CRISPR1 array. The CRISPR2 tree does not show such relation. Spacer mapping of the two CRISPR arrays suggests the construct to be canonical, with only 8.8% spacer conservation among the serovars. As opposed to broad-host-range serovars, the host-specific serovars harbor fewer spacers. All typhoidal serovars have CRISPR1-STY/cas-STY system. Comparison of CRISPR and cas phenograms with that of multilocus sequence typing (MLST) suggests differential evolution of CRISPR/Cas system implying supplementary roles beyond immunity.

CRISPRcasIdentifier: Machine learning for accurate identification and classification of CRISPR-Cas systems

Background CRISPR-Cas genes are extraordinarily diverse and evolve rapidly when compared to other prokaryotic genes. With the rapid increase in newly sequenced archaeal and bacterial genomes, manual identification of CRISPR-Cas systems is no longer viable. Thus, an automated approach is required for advancing our understanding of the evolution and diversity of these systems and for finding new candidates for genome engineering in eukaryotic models. Results We introduce CRISPRcasIdentifier, a new machine learning–based tool that combines regression and classification models for the prediction of potentially missing proteins in instances of CRISPR-Cas systems and the prediction of their respective subtypes. In contrast to other available tools, CRISPRcasIdentifier can both detect cas genes and extract potential association rules that reveal functional modules for CRISPR-Cas systems. In our experimental benchmark on the most recently published and comprehensive CRISPR-Cas system dataset, CRISPRcasIdentifier was compared with recent and state-of-the-art tools. According to the experimental results, CRISPRcasIdentifier presented the best Cas protein identification and subtype classification performance. Conclusions Overall, our tool greatly extends the classification of CRISPR cassettes and, for the first time, predicts missing Cas proteins and association rules between Cas proteins. Additionally, we investigated the properties of CRISPR subtypes. The proposed tool relies not only on the knowledge of manual CRISPR annotation but also on models trained using machine learning.

Detection of CRISPR cassettes and cas genes in the Arabidopsis thaliana genome

The state of the art in the evolution of plant viruses allows the genetic foundations of antiviral immunity in higher (including the most important crops) plants to be categorized as one of the most pressing issues of genetics and selection. According to the endosymbiotic theory, mitochondria descended from alphaproteobacteria that had been absorbed but not degraded by the host cell. The discovery of CRISPR-Cas systems (clustered regularly interspaced short palindromic repeats (CRISPR)-associated proteins), which implement the adaptive immunity function in prokaryotes, raises the question whether such a mechanism of antiviral protection could be caught up by evolution and used by representatives of eukaryotes (in particular, plants). The purpose of this work was to analyze the complete sequences of nuclear, mitochondrial, and chloroplast genomes of Arabidopsis thaliana in order to search for genetic elements similar to those in CRISPR-Cas systems of bacteria and archaea. As a result, in silico methods helped us to detect a locus of regularly intermittent short direct repeats in the mitochondrial genome of A. thaliana ecotypes. The structure of this locus corresponds to the CRISPR locus of the prokaryotic adaptive antiviral immune system. The probable connection between the locus found in the mitochondrial genome of the higher plant and the function of adaptive immunity is indicated by a similarity between the spacer sequences in the CRISPR cassette found and the genome of Cauliflower mosaic virus affecting Arabidopsis plants. Sequences of repeats and spacers of CRISPR cassettes in Arabidopsis C24 and Ler lines are perfectly identical. However, the locations of the CRISPR locus in the mitochondrial genomes of these lines differ significantly. The CRISPR cassette in the Col-0 line was found to be completely broken as a result of four deletions and one insertion. Although cas genes were not detected in the mitochondrial genome of the studied Arabidopsis ecotypes, their presence was detected in the nuclear genome. Both cas genes and numerous CRISPR cassettes were found on all the five chromosomes in the nuclear genome of the Col-0 ecotype. The results suggest the existence of a system of adaptive immunity in plants, which is similar to the CRISPR immunity of bacteria and archaea.