Comparative analysis of CRISPR cassettes from the human gut metagenomic contigs

Identification of the EH CRISPR-Cas9 system on a metagenome and its application to genome engineering

Non-coding RNAs (crRNAs) produced from clustered regularly interspaced short palindromic repeats (CRISPR) loci and CRISPR-associated (Cas) proteins of the prokaryotic CRISPR-Cas systems form complexes that interfere with the spread of transmissible genetic elements through Cas-catalysed cleavage of foreign genetic material matching the guide crRNA sequences. The easily programmable targeting of nucleic acids enabled by these ribonucleoproteins has facilitated the implementation of CRISPR-based molecular biology tools for in vivo and in vitro modification of DNA and RNA targets. Despite the diversity of DNA-targeting Cas nucleases so far identified, native and engineered derivatives of the Streptococcus pyogenes SpCas9 are the most widely used for genome engineering, at least in part due to their catalytic robustness and the requirement of an exceptionally short motif (5'-NGG-3' PAM) flanking the target sequence. However, the large size of the SpCas9 variants impairs the delivery of the tool to eukaryotic cells and smaller alternatives are desirable. Here, we identify in a metagenome a new CRISPR-Cas9 system associated with a smaller Cas9 protein (EHCas9) that targets DNA sequences flanked by 5'-NGG-3' PAMs. We develop a simplified EHCas9 tool that specifically cleaves DNA targets and is functional for genome editing applications in prokaryotes and eukaryotic cells.

CRISPR-Cas Systems in Gut Microbiome of Children with Autism Spectrum Disorders

The human gut microbiome is associated with various diseases, including autism spectrum disorders (ASD). Variations of the taxonomical composition in the gut microbiome of children with ASD have been observed repeatedly. However, features and parameters of the microbiome CRISPR-Cas systems in ASD have not been investigated yet. Here, we demonstrate such an analysis in order to describe the overall changes in the microbiome CRISPR-Cas systems during ASD as well as to reveal their potential to be used in diagnostics and therapy. For the systems identification, we used a combination of the publicly available tools suited for completed genomes with subsequent filtrations. In the considered data, the microbiomes of children with ASD contained fewer arrays per Gb of assembly than the control group, but the arrays included more spacers on average. CRISPR arrays from the microbiomes of children with ASD differed from the control group neither in the fractions of spacers with protospacers from known genomes, nor in the sets of known bacteriophages providing protospacers. Almost all bacterial protospacers of the gut microbiome systems for both children with ASD and the healthy ones were located in prophage islands, leaving no room for the systems to participate in the interspecies competition.

SCRAMBLER: A Tool for De Novo CRISPR Array Reconstruction and Its Application for Analysis of the Structure of Prokaryotic Populations

CRISPR arrays are prokaryotic genomic loci consisting of repeat sequences alternating with unique spacers acquired from foreign nucleic acids. As one of the fastest-evolving parts of the genome, CRISPR arrays can be used to differentiate closely related prokaryotic lineages and track individual strains in prokaryotic communities. However, the assembly of full-length CRISPR arrays sequences remains a problem. Here, we developed SCRAMBLER, a tool that includes several pipelines for assembling CRISPR arrays from high-throughput short-read sequencing data. We assessed its performance with model data sets (Escherichia coli strains containing different CRISPR arrays and imitating prokaryotic communities of different complexities) and intestinal microbiomes of extant and extinct pachyderms. Evaluation of SCRAMBLER's performance using model data sets demonstrated its ability to assemble CRISPR arrays correctly from reads containing pairs of spacers, yielding a precision rate of >80% and a recall rate of 60-85% when checked against ground-truth data. Likewise, SCRAMBLER successfully assembled CRISPR arrays from the environmental samples, as attested by their matching with database entries. SCRAMBLER, an open-source software (github.com/biolab-tools/SCRAMBLER), can facilitate analysis of the composition and dynamics of CRISPR arrays in complex communities.

Investigation into the prevalent CRISPR-Cas systems among the Aeromonas genus

Clustered regularly interspaced short palindromic repeats (CRISPR) is a prokaryotic adaptive immune system that checks invasion by mobile genetic elements through nuclease targeting. In this study, we investigated the occurrence, diversity, and features of the CRISPR system in the genus Aeromonas using bioinformatics tools. Only 13 out of 122 complete genomes (10.66%) of the genus Aeromonas from the NCBI GenBank database harbored the CRISPR system. The Type I-F system was the most prevalent CRISPR system among the Aeromonads, followed by the Type I-E system. Only one strain harbored a Type I-C CRISPR system. Among the Aeromonads, Aeromonas caviae (22.7%) and Aeromonas veronii (20%) had a higher prevalence rate of the complete CRISPR system. The analysis of direct repeat (DR) sequences showed that all could form stable RNA secondary structures. A phylogenetic tree generated for the Cas1 protein classified CRISPR subtypes into three distinct clusters. Among the 748 spacers investigated, 41.98% and 17.25% showed perfect homology to phage and plasmid sequences, respectively. Some arrays had duplicated spacers. The CRISPR loci are closely linked to antibiotic resistance genes in most strains. Collectively, our results would contribute to research on antibiotic resistance in the Aeromonas group, and provide new insights into the diversity and evolution of the CRISPR-Cas system.

Evaluation of CRISPR Diversity in the Human Skin Microbiome for Personal Identification

Microbial community diversity analysis can be utilized to characterize the personal microbiome that varies between individuals. CRISPR sequences, which reflect virome structure, in the human skin environment may be highly personalized similar to the structures of individual viromes.

The highly personalized human skin microbiome may serve as a viable marker in personal identification. Amplicon sequencing resolution using 16S rRNA cannot identify bacterial communities sufficiently to discriminate between individuals. Thus, novel higher-resolution genetic markers are required for forensic purposes. The clustered regularly interspaced short palindromic repeats (CRISPRs) are prokaryotic genetic elements that can provide a history of infections encountered by the bacteria. The sequencing of CRISPR spacers may provide phylogenetic information with higher resolution than other markers. However, using spacer sequencing for discrimination of personal skin microbiome is difficult due to limited information on CRISPRs in human skin microbiomes. It remains unclear whether personal microbiome discrimination can be achieved using spacer diversity or which CRISPRs will be forensically relevant. We identified common CRISPRs in the human skin microbiome via metagenomic reconstruction and used amplicon sequencing for deep sequencing of spacers. We successfully reconstructed 24 putative CRISPR arrays using metagenomic data sets. A total of 1,223,462 reads from three CRISPR arrays revealed that spacers in the skin microbiome were highly personalized, and conserved repeats were commonly shared between individuals. These individual specificities observed using CRISPR typing were confirmed by comparing the CRISPR diversity to microbiome diversity assessed using 16S rRNA amplicon sequencing. CRISPR typing achieved 95.2% accuracy in personal classification, whereas 16S rRNA sequencing only achieved 52.6%. These results suggest that sequencing CRISPRs in the skin microbiome may be a more powerful approach for personal identification and ecological studies compared to conventional 16S rRNA sequencing.

IMPORTANCE Microbial community diversity analysis can be utilized to characterize the personal microbiome that varies between individuals. CRISPR sequences, which reflect virome structure, in the human skin environment may be highly personalized similar to the structures of individual viromes. In this study, we identified 24 putative CRISPR arrays using a shotgun metagenome data set of the human skin microbiome. The findings of this study expand our understanding of the nature of CRISPRs by identifying novel CRISPR candidates. We developed a method to efficiently determine the diversity of three CRISPR arrays. Our analysis revealed that the CRISPR spacer diversity in the human skin microbiome is highly personalized compared with the microbiome diversity assessed by 16S rRNA sequencing, providing a new perspective on the study of the skin microbiome.

Identification of Natural CRISPR Systems and Targets in the Human Microbiome

Many bacteria resist invasive DNA by incorporating sequences into CRISPR loci, which enable sequence-specific degradation. CRISPR systems have been well studied from isolate genomes, but culture-independent metagenomics provide a new window into their diversity. We profiled CRISPR loci and cas genes in the body-wide human microbiome using 2,355 metagenomes, yielding functional and taxonomic profiles for 2.9 million spacers by aligning the spacer content to each sample's metagenome and corresponding gene families. Spacer and repeat profiles agree qualitatively with those from isolate genomes but expand their diversity by approximately 13-fold, with the highest spacer load present in the oral microbiome. The taxonomy of spacer sequences parallels that of their source community, with functional targets enriched for viral elements. When coupled with cas gene systems, CRISPR-Cas subtypes are highly site and taxon specific. Our analysis provides a comprehensive collection of natural CRISPR-cas loci and targets in the human microbiome.

Extraordinary diversity of viruses in deep-sea sediments as revealed by metagenomics without prior virion separation

Our current knowledge of the virosphere in deep-sea sediments remains rudimentary. Here we investigated viral diversity at both gene and genomic levels in deep-sea sediments of Southwest Indian Ocean. Analysis of 19 676 106 non-redundant genes from the metagenomic DNA sequences revealed a large number of unclassified viral groups in these samples. A total of 1106 high-confidence viral contigs were obtained after two runs of assemblies, and 217 of these contigs with sizes up to ~120 kb were shown to represent complete viral genomes. These contigs are clustered with no known viral genomes, and over 2/3 of the ORFs on the viral contigs encode no known functions. Furthermore, most of the complete viral contigs show limited similarity to known viral genomes in genome organization. Most of the classified viral contigs are derived from dsDNA viruses belonging to the order Caudovirales, including primarily members of the families Myoviridae, Podoviridae and Siphoviridae. Most of these viruses infect Proteobacteria and, less frequently, Planctomycetes, Firmicutes, Chloroflexi, etc. Auxiliary metabolic genes (AMGs), present in abundance on the viral contigs, appear to function in modulating the host ability to sense environmental gradients and community changes, and to uptake and metabolize nutrients.

Probing the Mobilome: Discoveries in the Dynamic Microbiome

There has been an explosion of metagenomic data representing human, animal, and environmental microbiomes. This provides an unprecedented opportunity for comparative and longitudinal studies of many functional aspects of the microbiome that go beyond taxonomic classification, such as profiling genetic determinants of antimicrobial resistance, interactions with the host, potentially clinically relevant functions, and the role of mobile genetic elements (MGEs). One of the most important but least studied of these aspects are the MGEs, collectively referred to as the 'mobilome'. Here we elaborate on the benefits and limitations of using different metagenomic protocols, discuss the relative merits of various sequencing technologies, and highlight relevant bioinformatics tools and pipelines to predict the presence of MGEs and their microbial hosts.

Host and body site-specific adaptation of Lactobacillus crispatus genomes

Lactobacillus crispatus is a common inhabitant of both healthy poultry gut and human vaginal tract, and the absence of this species has been associated with a higher risk of developing infectious diseases. In this study, we analyzed 105 L. crispatus genomes isolated from a variety of ecological niches, including the human vaginal tract, human gut, chicken gut and turkey gut, to shed light on the genetic and functional features that drive evolution and adaptation of this important species. We performed in silico analyses to identify the pan and core genomes of L. crispatus, and to reveal the genomic differences and similarities associated with their origins of isolation. Our results demonstrated that, although a significant portion of the genomic content is conserved, human and poultry L. crispatus isolates evolved to encompass different genomic features (e.g. carbohydrate usage, CRISPR-Cas immune systems, prophage occurrence) in order to thrive in different environmental niches. We also observed that chicken and turkey L. crispatus isolates can be differentiated based on their genomic information, suggesting significant differences may exist between these two poultry gut niches. These results provide insights into host and niche-specific adaptation patterns in species of human and animal importance.

CRISPR-Cas bioinformatics

Clustered regularly interspaced short palindromic repeats (CRISPR) and their associated proteins (Cas) are essential genetic elements in many archaeal and bacterial genomes, playing a key role in a prokaryote adaptive immune system against invasive foreign elements. In recent years, the CRISPR-Cas system has also been engineered to facilitate target gene editing in eukaryotic genomes. Bioinformatics played an essential role in the detection and analysis of CRISPR systems and here we review the bioinformatics-based efforts that pushed the field of CRISPR-Cas research further. We discuss the bioinformatics tools that have been published over the last few years and, finally, present the most popular tools for the design of CRISPR-Cas9 guides.

Comparative metagenomic analysis of microbial taxonomic and functional variations in untreated surface and reclaimed waters used in irrigation applications

The use of irrigation water sourced from reclamation facilities and untreated surface water bodies may be a practical solution to attenuate the burden on diminishing groundwater aquifers. However, comprehensive microbial characterizations of these water sources are generally lacking, especially with regard to variations through time and across multiple water types. To address this knowledge gap we used a shotgun metagenomic approach to characterize the taxonomic and functional variations of microbial communities within two agricultural ponds, two freshwater creeks, two brackish rivers, and three water reclamation facilities located in the Mid-Atlantic, United States. Water samples (n = 24) were collected from all sites between October and November 2016, and filtered onto 0.2 μm membrane filters. Filters were then subjected to total DNA extraction and shotgun sequencing on the Illumina HiSeq platform. From these data, we found that Betaproteobacteria dominated the majority of freshwater sites, while Alphaproteobacteria were abundant at times in the brackish waters. One of these brackish sites was also host to a greater abundance of the bacterial genera Gimesia and Microcystis. Furthermore, predicted microbial features (e.g. antibiotic resistance genes (ARGs) and Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) arrays) varied based on specific site and sampling date. ARGs were found across samples, with the diversity and abundance highest in those from a reclamation facility and a wastewater-impacted freshwater creek. Additionally, we identified over 600 CRISPR arrays, containing ∼2600 unique spacers, suggestive of a diverse and often site-specific phage community. Overall, these results provide a better understanding of the complex microbial community in untreated surface and reclaimed waters, while highlighting possible environmental and human health impacts associated with their use in agriculture.

CRISPRs for Strain Tracking and Their Application to Microbiota Transplantation Data Analysis

CRISPR-Cas systems are adaptive immune systems naturally found in bacteria and archaea. Prokaryotes use these immune systems to defend against invaders, which include phages, plasmids, and other mobile genetic elements. Relying on the integration of spacers derived from invader sequences (protospacers) into CRISPR loci (forming spacers flanked by repeats), CRISPR-Cas systems are able to store the memory of past immunological encounters. While CRISPR-Cas systems have evolved in response to invading mobile genetic elements, invaders have also developed mechanisms to avoid detection. As a result of an arms race between CRISPR-Cas systems and their targets, CRISPR arrays typically undergo rapid turnover of spacers through the acquisition and loss events. Additionally, microbiomes of different individuals rarely share spacers. Here, we present a computational pipeline, CRISPRtrack, for strain tracking based on CRISPR spacer content, and we applied it to fecal transplantation microbiome data to study the retention of donor strains in recipients. Our results demonstrate the potential use of CRISPRs as a simple yet effective tool for donor-strain tracking in fecal transplantation and as a general purpose tool for quantifying microbiome similarity.

Insights into the Human Virome Using CRISPR Spacers from Microbiomes

Due to recent advances in next-generation sequencing over the past decade, our understanding of the human microbiome and its relationship to health and disease has increased dramatically. Yet, our insights into the human virome, and its interplay with important microbes that impact human health, is relatively limited. Prokaryotic and eukaryotic viruses are present throughout the human body, comprising a large and diverse population which influences several niches and impacts our health at various body sites. The presence of prokaryotic viruses like phages, has been documented at many different body sites, with the human gut being the richest ecological niche. Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) and associated proteins constitute the adaptive immune system of bacteria, which prevents attack by invasive nucleic acid. CRISPR-Cas systems function by uptake and integration of foreign genetic element sequences into the CRISPR array, which constitutes a genomic archive of iterative vaccination events. Consequently, CRISPR spacers can be investigated to reconstruct interplay between viruses and bacteria, and metagenomic sequencing data can be exploited to provide insights into host-phage interactions within a niche. Here, we show how the CRISPR spacer content of commensal and pathogenic bacteria can be used to determine the evidence of their phage exposure. This framework opens new opportunities for investigating host-virus dynamics in metagenomic data, and highlights the need to dedicate more efforts for virome sampling and sequencing.

Freiburg RNA tools: a central online resource for RNA-focused research and teaching

The Freiburg RNA tools webserver is a well established online resource for RNA-focused research. It provides a unified user interface and comprehensive result visualization for efficient command line tools. The webserver includes RNA-RNA interaction prediction (IntaRNA, CopraRNA, metaMIR), sRNA homology search (GLASSgo), sequence-structure alignments (LocARNA, MARNA, CARNA, ExpaRNA), CRISPR repeat classification (CRISPRmap), sequence design (antaRNA, INFO-RNA, SECISDesign), structure aberration evaluation of point mutations (RaSE), and RNA/protein-family models visualization (CMV), and other methods. Open education resources offer interactive visualizations of RNA structure and RNA-RNA interaction prediction as well as basic and advanced sequence alignment algorithms. The services are freely available at http://rna.informatik.uni-freiburg.de.

The Revolution Continues: Newly Discovered Systems Expand the CRISPR-Cas Toolkit

CRISPR-Cas systems defend prokaryotes against bacteriophages and mobile genetic elements and serve as the basis for revolutionary tools for genetic engineering. Class 2 CRISPR-Cas systems use single Cas endonucleases paired with guide RNAs to cleave complementary nucleic acid targets, enabling programmable sequence-specific targeting with minimal machinery. Recent discoveries of previously unidentified CRISPR-Cas systems have uncovered a deep reservoir of potential biotechnological tools beyond the well-characterized Type II Cas9 systems. Here we review the current mechanistic understanding of newly discovered single-protein Cas endonucleases. Comparison of these Cas effectors reveals substantial mechanistic diversity, underscoring the phylogenetic divergence of related CRISPR-Cas systems. This diversity has enabled further expansion of CRISPR-Cas biotechnological toolkits, with wide-ranging applications from genome editing to diagnostic tools based on various Cas endonuclease activities. These advances highlight the exciting prospects for future tools based on the continually expanding set of CRISPR-Cas systems.

Characterization and Exploitation of CRISPR Loci in Bifidobacterium longum

Diverse CRISPR-Cas systems provide adaptive immunity in many bacteria and most archaea, via a DNA-encoded, RNA-mediated, nucleic-acid targeting mechanism. Over time, CRISPR loci expand via iterative uptake of invasive DNA sequences into the CRISPR array during the adaptation process. These genetic vaccination cards thus provide insights into the exposure of strains to phages and plasmids in space and time, revealing the historical predatory exposure of a strain. These genetic loci thus constitute a unique basis for genotyping of strains, with potential of resolution at the strain-level. Here, we investigate the occurrence and diversity of CRISPR-Cas systems in the genomes of various Bifidobacterium longum strains across three sub-species. Specifically, we analyzed the genomic content of 66 genomes belonging to B. longum subsp. longum, B. longum subsp. infantis and B. longum subsp. suis, and identified 25 strains that carry 29 total CRISPR-Cas systems. We identify various Type I and Type II CRISPR-Cas systems that are widespread in this species, notably I-C, I-E, and II-C. Noteworthy, Type I-C systems showed extended CRISPR arrays, with extensive spacer diversity. We show how these hypervariable loci can be used to gain insights into strain origin, evolution and phylogeny, and can provide discriminatory sequences to distinguish even clonal isolates. By investigating CRISPR spacer sequences, we reveal their origin and implicate phages and prophages as drivers of CRISPR immunity expansion in this species, with redundant targeting of select prophages. Analysis of CRISPR spacer origin also revealed novel PAM sequences. Our results suggest that CRISPR-Cas immune systems are instrumental in mounting diversified viral resistance in B. longum, and show that these sequences are useful for typing across three subspecies.

An Enrichment of CRISPR and Other Defense-Related Features in Marine Sponge-Associated Microbial Metagenomes

Many marine sponges are populated by dense and taxonomically diverse microbial consortia. We employed a metagenomics approach to unravel the differences in the functional gene repertoire among three Mediterranean sponge species, Petrosia ficiformis, Sarcotragus foetidus, Aplysina aerophoba and seawater. Different signatures were observed between sponge and seawater metagenomes with regard to microbial community composition, GC content, and estimated bacterial genome size. Our analysis showed further a pronounced repertoire for defense systems in sponge metagenomes. Specifically, clustered regularly interspaced short palindromic repeats, restriction modification, DNA phosphorothioation and phage growth limitation systems were enriched in sponge metagenomes. These data suggest that defense is an important functional trait for an existence within sponges that requires mechanisms to defend against foreign DNA from microorganisms and viruses. This study contributes to an understanding of the evolutionary arms race between viruses/phages and bacterial genomes and it sheds light on the bacterial defenses that have evolved in the context of the sponge holobiont.

Characterization of CRISPR RNA transcription by exploiting stranded metatranscriptomic data

CRISPR-Cas systems are bacterial adaptive immune systems, each typically composed of a locus of cas genes and a CRISPR array of spacers flanked by repeats. Processed transcripts of CRISPR arrays (crRNAs) play important roles in the interference process mediated by these systems, guiding targeted immunity. Here we developed computational approaches that allow us to characterize the expression of many CRISPRs in their natural environments, using community RNA-seq (metatranscriptomic) data. By exploiting public human gut metatranscriptomic data sets, we studied the expression of 56 repeat-sequence types of CRISPRs, revealing that most CRISPRs are transcribed in one direction (producing crRNAs). In rarer cases, including a type II system associated with Bacteroides fragilis, CRISPRs are transcribed in both directions. Type III CRISPR-Cas systems were found in the microbiomes, but metatranscriptomic reads were barely found for their CRISPRs. We observed individual-level variation of the crRNA transcription, and an even greater transcription of a CRISPR from the antisense strand than the crRNA strand in one sample. The orientations of CRISPR expression implicated by metatranscriptomic data are largely in agreement with prior predictions for CRISPRs, with exceptions. Our study shows the promise of exploiting community RNA-seq data for investigating the transcription of CRISPR-Cas systems.

Deciphering and shaping bacterial diversity through CRISPR

Phage and bacteria have engaged in a sustainable arms race, a seemingly endless conflict, since the beginning of time. CRISPR-Cas systems shape and generate environmental diversity through evolution of both predator and prey genomes. Indeed, the gain or loss of CRISPR-mediated immunity and genome maintenance can spark speciation in bacteria. Alternatively, turning CRISPR-Cas on the host by targeting chromosomal DNA has led to the development of next-generation smart antimicrobials and genetic screening and engineering technologies. Although the ability to target and cleave DNA in a sequence-specific manner is a powerful mechanism utilized by bacteria to fend off phage, plasmids, and potentially harmful nucleic acids, it is also a promising technology for programmable targeting of undesirable bacteria in microbiome consortia.

Computational prediction of CRISPR cassettes in gut metagenome samples from Chinese type-2 diabetic patients and healthy controls

Background

CRISPR has been becoming a hot topic as a powerful technique for genome editing for human and other higher organisms. The original CRISPR-Cas (Clustered Regularly Interspaced Short Palindromic Repeats coupled with CRISPR-associated proteins) is an important adaptive defence system for prokaryotes that provides resistance against invading elements such as viruses and plasmids. A CRISPR cassette contains short nucleotide sequences called spacers. These unique regions retain a history of the interactions between prokaryotes and their invaders in individual strains and ecosystems. One important ecosystem in the human body is the human gut, a rich habitat populated by a great diversity of microorganisms. Gut microbiomes are important for human physiology and health. Metagenome sequencing has been widely applied for studying the gut microbiomes. Most efforts in metagenome study has been focused on profiling taxa compositions and gene catalogues and identifying their associations with human health. Less attention has been paid to the analysis of the ecosystems of microbiomes themselves especially their CRISPR composition.

Results

We conducted a preliminary analysis of CRISPR sequences in a human gut metagenomic data set of Chinese individuals of type-2 diabetes patients and healthy controls. Applying an available CRISPR-identification algorithm, PILER-CR, we identified 3169 CRISPR cassettes in the data, from which we constructed a set of 1302 unique repeat sequences and 36,709 spacers. A more extensive analysis was made for the CRISPR repeats: these repeats were submitted to a more comprehensive clustering and classification using the web server tool CRISPRmap. All repeats were compared with known CRISPRs in the database CRISPRdb. A total of 784 repeats had matches in the database, and the remaining 518 repeats from our set are potentially novel ones.

Conclusions

The computational analysis of CRISPR composition based contigs of metagenome sequencing data is feasible. It provides an efficient approach for finding potential novel CRISPR arrays and for analysing the ecosystem and history of human microbiomes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-015-0248-x) contains supplementary material, which is available to authorized users.

Global transcription of CRISPR loci in the human oral cavity

Background

Clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs) are active in acquired resistance against bacteriophage and plasmids in a number of environments. In the human mouth, CRISPR loci evolve to counteract oral phage, but the expression of these CRISPR loci has not previously been investigated. We sequenced cDNA from CRISPR loci found in numerous different oral bacteria and compared with oral phage communities to determine whether the transcription of CRISPR loci is specifically targeted towards highly abundant phage present in the oral environment.

Results

We found that of the 529,027 CRISPR spacer groups studied, 88 % could be identified in transcripts, indicating that the vast majority of CRISPR loci in the oral cavity were transcribed. There were no strong associations between CRISPR spacer repertoires and oral health status or nucleic acid type. We also compared CRISPR repertoires with oral bacteriophage communities, and found that there was no significant association between CRISPR transcripts and oral phage, regardless of the CRISPR type being evaluated. We characterized highly expressed CRISPR spacers and found that they were no more likely than other spacers to match oral phage. By reassembling the CRISPR-bearing reads into longer CRISPR loci, we found that the majority of the loci did not have spacers matching viruses found in the oral cavities of the subjects studied. For some CRISPR types, loci containing spacers matching oral phage were significantly more likely to have multiple spacers rather than a single spacer matching oral phage.

Conclusions

These data suggest that the transcription of oral CRISPR loci is relatively ubiquitous and that highly expressed CRISPR spacers do not necessarily target the most abundant oral phage.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1615-0) contains supplementary material, which is available to authorized users.