- Open Access
CrisprVi: a software for visualizing and analyzing CRISPR sequences of prokaryotes
BMC Bioinformatics volume 23, Article number: 172 (2022)
Clustered regularly interspaced short palindromic repeats (CRISPR) and their spacers are important components of prokaryotic CRISPR-Cas systems. In order to analyze the CRISPR loci of multiple genomes more intuitively and comparatively, here we propose a visualization analysis tool named CrisprVi.
CrisprVi is a Python package consisting of a graphic user interface (GUI) for visualization, a module for commands parsing and data transmission, local SQLite and BLAST databases for data storage and a functions layer for data processing. CrisprVi can not only visually present information of CRISPR direct repeats (DRs) and spacers, such as their orders on the genome, IDs, start and end coordinates, but also provide interactive operation for users to display, label and align the CRISPR sequences, which help researchers investigate the locations, orders and components of the CRISPR sequences in a global view. In comparison to other CRISPR visualization tools such as CRISPRviz and CRISPRStudio, CrisprVi not only improves the interactivity and effects of the visualization, but also provides basic statistics of the CRISPR sequences, and the consensus sequences of DRs/spacers across the input strains can be inspected from a clustering heatmap based on the BLAST results of the CRISPR sequences hitting against the genomes.
CrisprVi is a convenient tool for visualizing and analyzing the CRISPR sequences and it would be helpful for users to inspect novel CRISPR-Cas systems of prokaryotes.
Prokaryotic (including half of bacterial and most archaeal) genomes contain CRISPR-Cas systems mainly composed of clustered regularly interspaced short palindromic repeats (CRISPR) array, the promoter for its transcription (the leader) and a set of CRISPR-associated (Cas) genes [1,2,3]. The CRISPR-Cas systems are formed by the fight against invasive genetic elements (e.g. phages and plasmids) during the evolution of prokaryotes, and they are sort of adaptive immunity that can block the infection process [4,5,6]. Due to the functional properties of the CRISPR-Cas systems, current research around them has expanded to a broad area, covering their origins and structures , prokaryotic population , CRISPR typing , gene editing [10, 11], etc.
As two key elements of the CRISPR-Cas system, direct repeats (DRs) of 24–50 bp and similar sized spacers between the DRs records what kind of RNA molecules can be derived from the array for activating the adaptive immunity . Thus, it is important to investigate the loci and sequences of the DRs and spacers comprehensively. For a prokaryotic genome, it is straight forward to predict its CRISPR array/sequences using computational methods, such as CRT , PILER-CR , CRISPRFinder , CRISPRCasFinder , MinCED , CRISPRDetect , MetaCRAST , and CRISPRdisco . And several CRISPR related databases, such as CRISPRdb , CRISPRI , CRISPRCasdb , also integrate programs for CRISPR identification, some of which are those mentioned above. However, there are few tools specifically designed to visualize, manipulate, and analyze the CRISPR arrays. In the past, Excel macro programs were commonly used to visualize the CRISPR sequences, but the programming is complicated and not interactive. In recent years, two computational tools, namely CRISPRviz  and CRISPRStudio , were developed for analyzing the CRISPR sequences interactively. For example, CRISPRviz can be used to predict, visualize and manipulate the CRISPR sequences by a web interface. Specifically, CRISPRviz can extract CRISPR DRs and spacers using MinCED , compare the sequence graphics visually and conduct alignment of spacer arrays. However, when the input includes a lot of strains with complex composition of CRISPR sequences, CRISPRviz will produce a variety of combinations of colors and symbols for the CRISPR sequences, which may be quite confusing. Unfortunately, CRISPRviz does not provide methods to help users analyze such complex scenario. Another problem of CRISPRviz is that it strongly depends on MinCED for CRISPR detection, which affects the accuracy of CRISPR visualization. Another tool, namely CRISPRStudio, can only present the spacers in a graphical way. Neither CRISPRviz nor CRISPRStudio provides comprehensive functions for the users to manipulate and analyze the DRs/spacers, such as changing colors of the visulization. For investigating and comparing the CRISPR loci of prokaryotic strains more intuitively and interactively, here we propose a novel visualization tool named CrisprVi.
Framework of CrisprVi
CrisprVi mainly consists of a graphic user interface (GUI) for visualization, a module for commands parsing and data transmission, local SQLite and BLAST [25, 26] databases for data storage and a functions layer supporting data processing (see Fig. 1). The GUI including input and output layers based on PyQt5  is used to load, access and manipulate the CRISPR sequences interactively. The annotation data of CRISPR sequences in general feature format (GFF, see documentation of CrisprVi)  is extracted and stored in the SQLite database. And optionally raw genomes sequences can be converted to BLAST databases and be further used to find consensus sequences of DRs/spacers. The functions layer of CrisprVi is composed of CrisprVi core algorithms, BLAST and several Python packages, such as PyQt5, pandas , Numpy , matplotlib , seaborn , Biopython , etc. The pandas and Numpy packages are used for data processing and computation while the matplotlib and seborn packages are used to visualize results of the statistical analysis and consensus DR/spacer sequence finding.
The GUI provides users with a convenient way to load the original data of CRISPR sequences, to present the data in various graphical ways, to manipulate the graphics, to conduct statistical analysis of the CRISPRs, and to find consensus DR/spacer sequences across strains.
The input layer consists of four modules, namely DataLoad.ui, Manipulate.ui, StatisticOp.ui and ConsSeqFind.ui, for the users to input data and operation commands through the graphic interface. To facilitate data loading, CrisprVi provides a wizard guiding the users to load GFF files and to set input/output directories step by step. First, the users need to load CRISPR annotation files in GFF, which contain all information of DRs and spacers of comparative strains, such as ID of the source genome, type (DR/spacer), coordinates on the genome, strand/direction, sequence and IDs. By default, the GFF file is the standard output by several CRISPR prediction programs, e.g. CRISPRCasFinder and CRISPRDetect. It is noting that the users should prepare the GFF files before starting CrisperVi. The CRISPR annotation information imported by the users will be stored in the SQLite database for subsequent query, delete, add and other operation.
After loading the GFF files of CRISPR annotations, the graphics of DRs and spacers will be shown in the main window by an output module Genomic.display, which includes three sub-modules, namely ‘DRs and Spacers’, ‘Spacers’ and ‘DRs’, for presenting CRISPR information in different ways. Then the users can directly manipulate the DRs and spacers by clicking on the graphics or specific buttons on the board provided by Manipulate.ui. The general manipulation includes showing CRISPR information (ID, source, start, end, strand, sequence, etc.), zooming in/out, deleting strains, sorting strains, aligning spacers, adding/deleting gaps.
Moreover, the users can conduct several statistical analysis on the CRISPR annotation data, such as counting DRs and spacers of strains, and calculating GC contents of DRs and spacers, etc., by commands provided by the module StatisticOp.ui. The results of the statistical analysis will be shown in figures by Statistics.display.
When the users are interested in the occurrences of consensus (identical or similar) sequences of DRs/spacers across strains, the module ConsSeqFind.ui can first find out the hits of DRs/spacers on the input genomes using BLAST, and then conduct clustering based on the BLAST results (see “Consensus sequences finding” section for details), which will generate a clustering heatmap through ConsPatt.display for the users to analyze the consensus patterns of CRISPR sequences.
Commands parsing and data transmission
CrisprVi provides several functions to parse different operation commands sent from the users via GUI before triggering different functions, and the results will be returned and displayed by the output layer. Meanwhile, the data channels are used for data transmission between the caller and specific executive functions.
The CRISPR information, such as id, source, start, end, strand, and sequence, in the input GFF files are extracted and stored in the SQLite database, where the users can add, delete, and search for the query sequences. Table 1 shows the format of the CRISPR database. On the other hand, the genome sequences of comparative strains will be converted to BLAST databases using commands provided by the BLAST package.
CrisprVi core algorithms
To support several functions of interactive operation of the CRISPR sequences, we integrate several core algorithms in CrisprVi as follows.
Alignment of spacer arrays
It is a modified pairwise alignment method named SpacerAlign for comparing the similarity of spacer arrays of input strains. Specifically, it first constructs a diagonal score matrix by conducting pairwise alignment (the Needleman-Wunsch algorithm [33, 34]) on each pair of spacer arrays. Then the matrix is used to construct an UPGMA tree , which guides the alignment of multiple spacer arrays progressively.
Consensus sequences finding
It is an algorithm designed for detecting consensus sequences of DRs/spacers across input strains using local BLAST and clustering. In the module of consensus sequences finding, each DR/spacer of the input strains is searched for in the genome database using BLAST. Then a DR/spacer by genome matrix (DG or SG matrix) is constructed by a calculation method based on the BLAST results. Specifically, for each BLAST output file, the highest product of hit score and identity percentage (%) fills a cell of the DG/SG matrix corresponding to each pair of DR and genome. After that, the constructed DG/SG matrix is fed to the seaborn package for showing a clustering heatmap, from which the users can find some patterns of the CRISPR consensus sequences across strains.
Several other algorithms for such as custom sequence extraction, SQLite data access and conversion and so on, are also included in the CrisprVi core algorithms.
To evaluate CrisprVi and other CRISPR visualization tools, we first prepared dataset-I containing core genomes of 12 Campylobacter coli (C. coli) and 12 Campylobacter jejuni (C. jejuni) strains (see Additional file 1 for summary of dataset-I), which were downloaded from NCBI (https://www.ncbi.nlm.nih.gov/)  in March 2020. To evaluate the tools in visualizing large datasets, we prepared dataset-II containing 100 prokaryotic DNA sequences (see Additional file 2 for summary of dataset-II), which were downloaded from NCBI in Jun 2021. To conduct the comparison more fairly, we used two strategies as follows to obtain the CRISPR annotations. The strain counts of dataset-I and dataset-II are summarized in Table 2.
CRISPR annotation based on CRISPRCasFinder
The online tool CRISPRCasFinder  was used to generate GFF files containing CRISPR annotations for 12 C. coli and 12 C. jejuni of dataset-I respectively. Then the GFF files were input to CrisprVi for visualization and analysis. It should be noticed that the GFF files cannot be import to CRISPRviz, which only requires the FASTA files of genomes as input.
CRISPR annotation based on MinCED
On the other hand, we imported the genomes of dataset-I to CRISPRviz by running its crisprviz.sh script (with parameters -pxc), which calls MinCED  to annotate CRISPRs. It is noticed that MinCED just generated CRISPR annotation files for 9 C. coli and 8 C. jejuni strains on dataset-I, which were all converted to GFF files using the minced2gff.py script packaged in CrisprVi for further loading to CrisprVi (Table 2). On dataset-II, we used the same strategy based on MinCED to obtain 80 CRISPR annotation files, which were all converted to GFF files for inputting to CrisprVi.
CrisprVi provides GUI for analyzing CRISPR sequences
CrisprVi is a graphically interactive software including several function modules, such as showing DRs and/or spacers, alignment of spacers, statistics of DRs and/or spacers, and consensus patterns of DRs/spacers, as shown in Fig. 2. In comparison with other visualization tools such as CRISPRviz that can predict the CRISPR sequences using specific build-in tool, our CrisprVi just focuses on visualizing the CRISPR sequences of comparative strains extracted from the CRISPR annotations in GFF. Thus, CrisprVi can be more flexible as the users can choose any suitable tool to predict the CRISPR sequences before further visualization.
Interactive operation of DRs and/or spacers
By clicking on different buttons/menus on the panel of CrisprVi, the users can change and manipulate the graphics denoting the DRs/spacers in several ways, such as showing DRs and spacers together or separately, and aligning spacer arrays (see Fig. 2). In CrisprVi, each spacer is denoted by a colored rectangular graphic with an inner number plus strand inside while the DR is represented by a colored diamond graphic with an inner number followed by an apostrophe plus strand inside. It is noting that the inner number is assigned to the DRs/spacers having the same sequence. To discriminate between different DRs/spacers visually, CrisprVi automatically assigns different combinations of colors and numbers for the DRs/spacers with different nucleotide composition. If the users do not like the color automatically assigned, they can right-click on the graphic and change the color by the color palette (see Fig. 3A). In the display area, the DRs and/or spacers of each strain or CRISPR is displayed in a track, which can be deleted selectively. And the users can sort the tracks according to their lengths (see Fig. 3B), namely the counts of DRs and/or spacers of the tracks. Once moving the mouse arrow onto specific DR/spacer, the information of the DR/spacer will be displayed over the graphic. Likewise, if the users click on the graphic, the information of the DR/spacer will be displayed on the panel under the scrollbar area, and all graphics having the same sequence as the selected one will be highlighted by red boarders (Fig. 3C), which can help the users view the global distribution of specific DR/spacer across strains more clearly.
Statistical analysis of DRs and spacers
To obtain basic statistics of the DRs/spacers quickly, the users can perform statistical analysis of the DRs/spacers on the selected strains of interest (Figs. 4 and 5). First, the DRs occurred in the strains can be counted and displayed in histograms (see Fig. 4). As seen from the histogram in Fig. 4A, the repeat DR2 (type = ‘dr’, inner number = ‘2’) occurs the most frequently in all of the 12 C. coli strains while DR6 (type = ‘dr’, inner number = ‘6’) of C. jejuni is the most frequent repeat (Fig. 4B). The original source IDs of the DR with an inner number can be traced back by clicking on the graphics with the same inner number in other visualization modules mentioned previously. Second, the spacer counts across strains can be calculated and visualized in histograms, as shown in Fig. 5. Third, the users can compare the distribution of GC contents between DRs and spacers globally for the selected strains using boxplots (Fig. 6). In addition, the users can inspect details of the plots using commands in the menu under the plot. As seen from the statistical figures of the CRISPRs of comparative C. coli and C. jejuni strains from dataset-I, the C. coli strains represented more diverse properties than C. jejuni in terms of DRs/spacers categories and frequencies. Previous study also showed that the C. coli isolates presented more diverse allelic distribution than C. jejuni . All these discoveries reveal the evolutionary differences between C. coli and C. jejuni strains.
Analysis of consensus DR/spacer sequences
CrisprVi provides a method for analyzing consensus sequences of DRs/spacers in all input genomes based on BLAST and clustering (see “Consensus sequences finding” section). Take dataset-I as an example, we first established local BLAST libraries for the genomes of C. coli and C. jejuni strains respectively, and then CrisprVi was used to align the DR/spacer sequences against the genomes in the libraries by calling the local BLAST. Based on the BLAST results, a DG/SG matrix was generated and visualized in a clustering heatmap, from which we can see how consensus each DR/spacer sequence occurs in the input strains, as shown in Fig. 7. First, the resulting heatmaps can be used to compare consensus sequence patterns across strains between DRs/spacers. For example, the DR sequences of the 12 C. coli strains represent a more diverse pattern than that of the C. jejuni strains. Second, the heatmap can help detect consensus pattern of some DR/spacer across the input strains. For example, Fig. 7C illustrates that two DR sequences (inner number = ‘1’ and ‘6’) represented consensus patterns across the C. jejuni strains. Then we traced back to their IDs and sequences using CrisprVi, which resulted in ‘1’ for ‘DR_1433807’ etc.
and ‘6’ for ‘DR_1533204’ etc.
Furthermore, we used the two sequences to query the local or online BLAST databases to validate the discovery. The blastn results showed that both of the sequences are C. jejuni specific without being found in C. coli or other species, which would be valuable for further study.
Comparison of CRISPR visualization methods
Several CRISPR visualization methods, e.g. CRISPRviz and CRISPRStudio, integrate CRISPR finding algorithms, which makes it difficult to compare these tools with our CrisprVi directly since current CrisprVi version does not contain the module for CRISPR finding. Since CRISPRStudio cannot visualize DR arrays, here we conducted comparison experiments between CrisprVi and CRISPRviz on dataset-I and dataset-II. As seen from Figs. 8 and 9, CrisprVi and CRISPRviz can both represent DRs and spacers in several ways, but they differ in some aspects.
CrisprVi and CRISPRviz show different color-symbol systems and operation
CrisprVi and CRISPRviz use different color-symbol systems and operation. Specifically, CRISPRviz displays the DRs and spacers in random combinations of colors and symbols, whereas in CrisprVi the colors of DR/spacer graphics are automatically assigned and can be changed manually (Fig. 3A). In addition, CrisprVi also uses inner number plus strand to denote one type of DRs/spacers uniquely (Fig. 3B). The color-symbol system of CRISPRviz or CrisprVi would become more and more complex as the number of DRs/spacers increases while CrisprVi can highlight the identical DRs/spacers by red boarders (Fig. 3C), which can help the users view the DRs/spacers of interest more clearly. Furthermore, the inner number used by CrisprVi are more useful for tracing back the meta data of the CRISPRs in such complex scenario.
CrisprVi accepts output from different CRISPR finding methods
Different CRISPR finding methods, e.g. CRISPRCasFinder  and MinCED , may give different prediction results. Once the users choose CRISPRviz for CRISPR visualization, they must acquiesce in the prediction accuracy of MinCED called by CRISPRviz. In comparison, CrisprVi does not connect to any specific CRISPR finding method, but the input files should follow the GFF format defined by CrisprVi or CRISPRCasFinder (see documentation of CrisprVi). Thus, the output GFF files of CRISPRCasFinder can be directly loaded to CrisprVi, while the CRISPR annotation files output by CRISPRviz/MinCED should be converted to our GFF format using the script minced2gff.py packaged in CrisprVi. Such conversion is important for the visualization comparison between CRISPRviz and CrisprVi as we can see the differences of CRISPR annotation in a visual way. As seen from Fig. 9, CrisprVi can arrange the two CRISPR arrays of a C. coli strain (NZ_CP007183_1) in the same track (Fig. 9A and B) while CRISPRviz only shows the two CRISPRs in two different tracks (Fig. 9C). The benefit of showing all CRISPRs of a strain in a line in CrisprVi is that the users can see the relative positions of the CRISPRs, and the boundary of two CRISPRs is just at the point of two connected DRs (Fig. 9A and B). Optionally, the users can arrange different CRISPR arrays of a strain in different tracks in CrisprVi as in CRISPRviz by selecting the CRISPR button on the board. Obviously, CRISPRCasFinder and MinCED gave different CRISPR prediction results on the same strain (e.g. NZ_CP007183_1). Moreover, CRISPRviz cannot show the array orientation, which is ignored by MinCED. Thus, the CrisprVi visualization based on MinCED prediction uses dot (‘.’) to indicate such missing information. But CRISPRviz does not present such comprehensive information of the CRISPR arrays as CrisprVi.
Differences in visualizing DRs/spacers
In terms of visualizing the DRs/spacers, CRISPRviz and CrisprVi have different operation. For example, Fig. 10 illustrates a DR array composed of five repeats in the strain with ID ‘NC_022132.1’ (see Table 3 for details). As seen from Table 3, the DR array just has three types of sequences with slight differences, two of which are represented by the same color-symbol in CRISPRviz (Fig. 10A). Although the DRs with the same color-symbol (the green graphics in Fig. 10A) implies that they are similar to each other, their differences has been hidden (the sequences of No. 2 to 5 in Table 3), especially when the DR array have great diversity of DRs. Unfortunately, in such scenario the users cannot find any functions in CRISPRviz to inspect the differences between the similar DRs, which can only be dug out from the original annotation file. In comparison, CrisprVi assigns different color-number combinations to the DRs (Fig. 10B), which makes it easy to distinguish between the different DRs. In CrisprVi, the nucleotide composition of specific DR can be showed on the panel.
Visualizing large datasets
Both CrisprVi and CRISPRviz can be used to visualize large datasets, which was validated on dataset-II containing 80 available strains with CRISPRs annotated by MinCED (see Additional files 3 and 4). Meanwhile, the visualization differences between the two tools still exist in such scenario.
Running time comparison
In terms of running time, CrisprVi and CRISPRviz include different components of running time. Specifically, CRISPRviz mainly includes spending time for DR/Spacer detection and data loading via web service while our CrisprVi include the time for data preprocessing, data loading locally, statistical analysis, and consensus sequencing finding. Thus, we only compared CrisprVi and CRISPRviz in terms of data loading via local network (see Table 4) on a laptop (Win10 64bit, Intel Core i5 CPU, 8 GB RAM). Meanwhile, CRISPRviz ran CRISPR detection on the same laptop with virtual Linux OS (Ubuntu 18.04). Five sub-datasets with varying size (20, 40, 60, 80 and 100) were generated by sampling from the original 100 strains of dataset-II. The CRISPRviz was run on each sub-dataset for visualization. Meanwhile, the CRISPR annotation files generated by MinCED/CRISPRviz were converted to GFF files, which were input to CrisprVi for visualization. The results manifest that CRISPRviz spent a lot of time for CRISPR detection and its data loading was much faster than our CrisprVi. Since CrisprVi integrates several types of functional objects for modelling, manipulating and analyzing the CRISPRs, and the number of the objects as well as the time for data loading increase as the sample size. The results indicate that the underlying efficiency of our visualization framework based on PyQt5 should be improved in the future.
Summary of CrisprVi, CRISPRviz and CRISPRStudio visualization
In compared with CrisprVi and CRISPRviz, CRISPRStudio just implements visualization of the spacers. Table 5 summarizes the similarities and differences of the main functions between them in terms of input files, DRs and spacers extraction methods, color-symbol system, alignment of spacer arrays, guide tree, statistics, etc.
Since the role of CRISPR-Cas system was revealed a decade ago, a diversity of CRISPR-Cas systems with various physiological functions have been discovered in prokaryotes. Currently the CRISPR-Cas systems can be categorized into two classes with six types and additional subtypes, and new types might remain to be discovered . From an evolutionary view, it is of importance to investigate the CRISPR distribution across strains of interest. Fortunately, computational tools such as CRISPRviz and CRISPRStudio have provided researchers with a large scale analysis to find and visualize the CRISPR arrays . However, the functions of current tools are not complete enough and few tools are available for CRISPR visualization. Therefore, novel tools for CRISPR analysis are still in demand.
To provide an efficient CRISPR visualization tool to the community, here we proposed a Python package named CrisprVi. By the GUI of CrisprVi, the information of CRISPR sequences including DRs and spacers can be visually displayed, and the users can also manipulate the CRISPR arrays on the screen. In comparison to other tools, CrisprVi not only improves functions and effects for CRISPR visualization, but also helps conduct several statistical analysis of the CRISPR sequences, and show the results in different ways. In addition, CrisprVi provides a method for analysing consensus sequences of DRs/spacers based on BLAST and clustering. Overall, CrisprVi can be used to compare and analyze the CRISPR information of multiple strains intuitively. In the experiment, CrisprVi was successfully tested on different magnitude of datasets and scenarios, which also indicates that our tool can be used for analysing other prokaryotic CRISPR strains. Although CrisprVi does not provide the function of CRISPR detection, it could be more flexible than those one-step processing tools that hides CRISPR detection, as the users can choose any available CRISPR finding tool and conduct curation on the automatically predicted results before CrisprVi visualization. In the future, CrisprVi is planned to be integrated with methods for complete CRISPR-Cas system identification and strain typing, including CRISPR finding, cas genes detection, etc.
CrisprVi is a convenient tool for visualizing and analyzing the CRISPR sequences and it would be helpful for researchers to inspect novel CRISPR-Cas systems of prokaryotes.
Availability and requirements
Project name: CrisprVi
Project home page: https://sourceforge.net/projects/crisprvi
Operating system(s): Platform independent
Programming language: Python
Other requirements: Python 3.7, PyQt5, Numpy, pandas, matplotlib, seaborn, Biopython, BLAST 2.9.0
License: GNU GPL (version 3)
Any restrictions to use by non-academics: No
Availability of data and materials
The datasets analyzed during the current study are available at: https://sourceforge.net/projects/crisprvi/files/data/dataset-1.zip and https://sourceforge.net/projects/crisprvi/files/data/dataset-2.zip.
Clustered regularly interspaced short palindromic repeats
Graphic user interface
Basic Local Alignment Search Tool
General feature format
- C. coli :
- C. jejuni :
Lander ES. The heroes of CRISPR. Cell. 2016;164(1):18–28.
Westra ER, Buckling A, Fineran PC. CRISPR-Cas systems: beyond adaptive immunity. Nat Rev Microbiol. 2014;12(5):317–26.
Barrangou R, Marraffini L. CRISPR-Cas systems: prokaryotes upgrade to adaptive immunity. Mol Cell. 2014;54(2):234–44.
Barrangou R, Fremaux C, Deveau H, Richards M, Boyaval P, Moineau S, et al. CRISPR provides acquired resistance against viruses in prokaryotes. Science. 2007;315(5819):1709–12.
Barrangou R, Horvath P. CRISPR: new horizons in phage resistance and strain identification. Annu Rev Food Sci Technol. 2012;3(1):143.
Deveau H, Barrangou R, Garneau JE, Labonte J, Fremaux C, Boyaval P, et al. Phage response to CRISPR-encoded resistance in Streptococcus thermophilus. J Bacteriol. 2008;190(4):1390–400.
Ogrodzki P, Forsythe SJ. CRISPR–cas loci profiling of Cronobacter sakazakii pathovars. Future Microbiol. 2016;00(12):1507–19.
Rodriguez-Valera F, Martin-Cuadrado AB, Rodriguez-Brito B, Pasić L, Thingstad TF, Rohwer F, et al. Explaining microbial population genomics through phage predation. Nat Rev Microbiol. 2009;7(11):828–36.
Mokrousov I, Rastogi N. Spacer-based macroarrays for CRISPR genotyping. In: Lundgren M, Charpentier E, Fineran PC, editors. CRISPR: methods and protocols. New York: Springer; 2015. p. 111–31.
Doudna JA, Charpentier E. The new frontier of genome engineering with CRISPR-Cas9. Science. 2014;346(6213):1258096.
Ran FA, Hsu PD, Wright J, Agarwala V, Scott DA, Zhang F. Genome engineering using the CRISPR-Cas9 system. Nat Protoc. 2013;8(11):2281–308.
Couvin D, Bernheim A, Toffano-Nioche C, Touchon M, Michalik J, Néron B, et al. CRISPRCasFinder, an update of CRISRFinder, includes a portable version, enhanced performance and integrates search for Cas proteins. Nucleic Acids Res. 2018;46(W1):W246–51.
Bland C, Ramsey TL, Sabree F, Lowe M, Brown K, Kyrpides NC, et al. CRISPR Recognition Tool (CRT): a tool for automatic detection of clustered regularly interspaced palindromic repeats. BMC Bioinformatics. 2007;8(1):209.
Edgar RC. PILER-CR: fast and accurate identification of CRISPR repeats. BMC Bioinform. 2007;8(1):18.
Grissa I, Vergnaud G, Pourcel C. CRISPRFinder: a web tool to identify clustered regularly interspaced short palindromic repeats. Nucleic Acids Res. 2007;35(Web Server issue):W52–7.
Skennerton CT, Soranzo N, Angly F. MinCED. https://github.com/ctSkennerton/minced. Accessed 31 Mar 2021.
Biswas A, Staals RHJ, Morales SE, Fineran PC, Brown CM. CRISPRDetect: a flexible algorithm to define CRISPR arrays. BMC Genomics. 2016;17(1):356.
Moller AG, Liang C. MetaCRAST: reference-guided extraction of CRISPR spacers from unassembled metagenomes. PeerJ. 2017;5: e3788.
Crawley AB, Henriksen JR, Barrangou R. CRISPRdisco: an automated pipeline for the discovery and analysis of CRISPR-Cas systems. The CRISPR J. 2018:crispr.2017.0022.
Grissa I, Vergnaud G, Pourcel C. The CRISPRdb database and tools to display CRISPRs and to generate dictionaries of spacers and repeats. BMC Bioinform. 2007;8(1):172.
Rousseau C, Gonnet M, Le Romancer M, Nicolas J. CRISPI: a CRISPR interactive database. Bioinformatics. 2009;25(24):3317–8.
Pourcel C, Touchon M, Villeriot N, Vernadet J-P, Couvin D, Toffano-Nioche C, et al. CRISPRCasdb a successor of CRISPRdb containing CRISPR arrays and cas genes from complete genome sequences, and tools to download and query lists of repeats and spacers. Nucleic Acids Res. 2019;48(D1):D535–44.
Nethery MA, Barrangou R. CRISPR Visualizer: rapid identification and visualization of CRISPR loci via an automated high-throughput processing pipeline. RNA Biol. 2019;16(4):577–84.
Dion MB, Labrie SJ, Shah SA, Moineau S. CRISPRStudio: a user-friendly software for rapid CRISPR array visualization. Viruses. 2018;10(11):602.
Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, et al. BLAST+: architecture and applications. BMC Bioinform. 2009;10:421.
Zhang Z, Schwartz S, Wagner L, Miller W. A greedy algorithm for aligning DNA sequences. J Comput Biol. 2000;7(1–2):203–14.
Willman JM. Getting started with PyQt: beginning PyQt; 2020.
McKinney W, et al., editors. Data structures for statistical computing in python. In: Proceedings of the 9th python in science conference; 2010.
Harris CR, Millman KJ, van der Walt SJ, Gommers R, Virtanen P, Cournapeau D, et al. Array programming with NumPy. Nature. 2020;585(7825):357–62.
Hunter JD. Matplotlib: a 2D graphics environment. Comput Sci Eng. 2007;9(3):90–5.
Waskom M, The seaborn development t. mwaskom/seaborn. Zenodo; 2020.
Cock PJA, Antao T, Chang JT, Chapman BA, Cox CJ, Dalke A, et al. Biopython: freely available Python tools for computational molecular biology and bioinformatics. Bioinformatics. 2009;25(11):1422–3.
Koyutürk M, Kim Y, Topkara U, Subramaniam S, Szpankowski W, Grama A. Pairwise alignment of protein interaction networks. J Comput Biol. 2006;13(2):182–99.
Needleman SB. A general method applicable to the search for similarities in the amino acid sequence of two proteins. J Mol Biol. 1970;48(3):443–53.
Sokal RR, Michener CD. A statistical method for evaluating systematic relationships. Univ Kansas Sci Bull. 1958;38:1409–38.
Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 2018;46(D1):D8–D13.
Sheppard SK, Maiden MCJ. The evolution of Campylobacter jejuni and Campylobacter coli. Cold Spring Harb Perspect Biol. 2015;7(8): a018119.
Makarova KS, Wolf YI, Iranzo J, Shmakov SA, Alkhnbashi OS, Brouns SJJ, et al. Evolutionary classification of CRISPR-Cas systems: a burst of class 2 and derived variants. Nat Rev Microbiol. 2020;18(2):67–83.
Charpentier E, Elsholz A, Marchfelder A. CRISPR-Cas: more than ten years and still full of mysteries. RNA Biol. 2019;16(4):377–9.
Pearson W, Crusoe MR. fasta36. https://github.com/wrpearson/fasta36. Accessed 31 Mar 2021.
The authors thank Dr. Jie Ding and Dr. Yunliang Yue for their suggestions and help on this study.
About this supplement
This article has been published as part of BMC Bioinformatics Volume 23 Supplement 3, 2022: Selected articles from the International Conference on Intelligent Biology and Medicine (ICIBM 2021): bioinformatics. The full contents of the supplement are available online at https://bmcbioinformatics.biomedcentral.com/articles/supplements/volume-23-supplement-3.
The design of the study is funded by National Key Research and Development Program of China (2018YFD0500500) and Open Project Program of Jiangsu Key Laboratory of Zoonosis (R1805). Collection, analysis, and interpretation of data and manuscript writing are supported by Natural Science Foundation of Jiangsu Province (BK20190878) and Cross-Disciplinary Project of the Animal Science Special Discipline of Yangzhou University (yzuxk202015). Publication costs are funded by Open Project Program of Jiangsu Key Laboratory of Zoonosis (R1805).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1. Table S1:
Summary of dataset-I. Dataset-I includes core genomes of 12 Campylobacter coli (C. coli) and 12 Campylobacter jejuni (C. jejuni) strains.
Additional file 2. Table S2:
Summary of dataset-II. Dataset-II includes DNA sequences of 100 prokaryotic strains.
Additional file 3. Fig. S1:
Snapshot of visualizing CRISPRs of 80 strains on CrisprVi.
Additional file 4. Fig. S2:
Snapshot of visualizing CRISPRs of 80 strains on CRISPRviz.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Sun, L., Wang, J., Yan, F. et al. CrisprVi: a software for visualizing and analyzing CRISPR sequences of prokaryotes. BMC Bioinformatics 23 (Suppl 3), 172 (2022). https://doi.org/10.1186/s12859-022-04716-9
- Direct repeat
- Consensus sequence