CRISPR Visualizer: rapid identification and visualization of CRISPR loci via an automated high-throughput processing pipeline

Unleashing the potential: type I CRISPR-Cas systems in actinomycetes for genome editing

Covering: up to the end of 2023Type I CRISPR-Cas systems are widely distributed, found in over 40% of bacteria and 80% of archaea. Among genome-sequenced actinomycetes (particularly Streptomyces spp.), 45.54% possess type I CRISPR-Cas systems. In comparison to widely used CRISPR systems like Cas9 or Cas12a, these endogenous CRISPR-Cas systems have significant advantages, including better compatibility, wide distribution, and ease of operation (since no exogenous Cas gene delivery is needed). Furthermore, type I CRISPR-Cas systems can simultaneously edit and regulate genes by adjusting the crRNA spacer length. Meanwhile, most actinomycetes are recalcitrant to genetic manipulation, hindering the discovery and engineering of natural products (NPs). The endogenous type I CRISPR-Cas systems in actinomycetes may offer a promising alternative to overcome these barriers. This review summarizes the challenges and recent advances in CRISPR-based genome engineering technologies for actinomycetes. It also presents and discusses how to establish and develop genome editing tools based on type I CRISPR-Cas systems in actinomycetes, with the aim of their future application in gene editing and the discovery of NPs in actinomycetes.

CRISPR genotyping methods: Tracing the evolution from spoligotyping to machine learning

CRISPR (clustered regularly interspaced short palindromic repeats)-Cas (CRISPR-associated) systems provide prokaryotes with adaptive immunity defenses against foreign genetic invaders. The identification of CRISPR-Cas function is among the most impactful discoveries of recent decades that have shaped the development of genome editing in various organisms paving the way for a plethora of promising applications in biotechnology and health. Even before the discovery of CRISPR-Cas biological role, the particular structure of CRISPR loci has been explored for epidemiological genotyping of bacterial pathogens. CRISPR-Cas loci are arranged in CRISPR arrays of mostly identical direct repeats intercalated with invader-derived spacers and an operon of cas genes encoding the Cas protein components. Each small CRISPR RNA (crRNA) encoded within the CRISPR array constitutes a key functional unit of this RNA-based CRISPR-Cas defense system guiding the Cas effector proteins toward the foreign nucleic acids for their destruction. The information acquired from prior invader encounters and stored within CRISPR arrays turns out to be extremely valuable in tracing the microevolution and epidemiology of major bacterial pathogens. We review here the history of CRISPR-based typing strategies highlighting the first PCR-based methods that have set the stage for recent developments of high-throughput sequencing and machine learning-based approaches. A great amount of whole genome sequencing and metagenomic data accumulated in recent years opens up new avenues for combining experimental and computational approaches of high-resolution CRISPR-based typing.

Comparative genomics of Loigolactobacillus coryniformis with an emphasis on L. coryniformis strain FOL-19 isolated from cheese

Loigolactobacillus coryniformis is a member of lactic acid bacteria isolated from various ecological niches. We isolated a novel L. coryniformis strain FOL-19 from artisanal Tulum cheese and performed the whole-genome sequencing for FOL-19. Then, genomic characterization of FOL-19 against ten available whole genome sequences of the same species isolated from kimchi, silage, fermented meat, air of cowshed, dairy, and pheasant chyme was performed to uncover the genetic diversity and biotechnological potential of overall species. The average genome size of 2.93 ± 0.1 Mb, GC content of 42.96% ± 0.002, number of CDS of 2905 ± 165, number of tRNA of 56 ± 10, and number of CRISPR elements of 6.55 ± 1.83 was found. Both Type I and II Cas clusters were observed in L. coryniformis. No bacteriocin biosynthesis gene clusters were found. All strains harbored at least one plasmid except KCTC 3167. All strains were predicted to carry multiple IS elements. The most common origin of the IS elements was belong to Lactiplantibacillus plantarum. Comparative genomic analysis of L. coryniformis revealed hypervariability at the strain level and the presence of CRISPR/Cas suggests that L. coryniformis holds a promising potential for being a reservoir for new CRISPR-based tools. All L. coryniformis strains except PH-1 were predicted to harbor pdu and cbi-cob-hem gene clusters encoding industrially relevant traits of reuterin and cobalamin biosynthesis, respectively. These findings put a step forward for the genomic characterization of L. coryniformis strains for biotechnological applications via genome-guided strain selection to identify industrially relevant traits.

Genomic and phenotypic characterization of Salmonella enterica serovar Kentucky

Non-typhoidal Salmonella are extremely diverse and different serovars can exhibit varied phenotypes, including host adaptation and the ability to cause clinical illness in animals and humans. In the USA, Salmonella enterica serovar Kentucky is infrequently found to cause human illness, despite being the top serovar isolated from broiler chickens. Conversely, in Europe, this serovar falls in the top 10 serovars linked to human salmonellosis. Serovar Kentucky is polyphyletic and has two lineages, Kentucky-I and Kentucky-II; isolates belonging to Kentucky-I are frequently isolated from poultry in the USA, while Kentucky-II isolates tend to be associated with human illness. In this study, we analysed whole-genome sequences and associated metadata deposited in public databases between 2017 and 2021 by federal agencies to determine serovar Kentucky incidence across different animal and human sources. Of 5151 genomes, 90.3 % were from isolates that came from broilers, while 5.9 % were from humans and 3.0 % were from cattle. Kentucky-I isolates were associated with broilers, while isolates belonging to Kentucky-II and a new lineage, Kentucky-III, were more commonly associated with cattle and humans. Very few serovar Kentucky isolates were associated with turkey and swine sources. Phylogenetic analyses showed that Kentucky-III genomes were more closely related to Kentucky-I, and this was confirmed by CRISPR-typing and multilocus sequence typing (MLST). In a macrophage assay, serovar Kentucky-II isolates were able to replicate over eight times better than Kentucky-I isolates. Analysis of virulence factors showed unique patterns across these three groups, and these differences may be linked to their association with different hosts.

CRISPR Comparison Toolkit: Rapid Identification, Visualization, and Analysis of CRISPR Array Diversity

CRISPR-Cas systems provide immunity against mobile genetic elements (MGEs) through sequence-specific targeting by spacer sequences encoded in CRISPR arrays. Spacers are highly variable between microbial strains and can be acquired rapidly, making them well suited for use in strain typing of closely related organisms. However, no tools are currently available to automate the process of reconstructing strain histories using CRISPR spacers. We therefore developed the CRISPR Comparison Toolkit (CCTK) to enable analyses of array relationships. The CCTK includes tools to identify arrays, analyze relationships between arrays using CRISPRdiff and CRISPRtree, and predict targets of spacers. CRISPRdiff visualizes arrays and highlights the similarities between them. CRISPRtree infers a phylogenetic tree from array relationships and presents a hypothesis of the evolutionary history of the arrays. The CCTK unifies several CRISPR analysis tools into a single command line application, including the first tool to infer phylogenies from array relationships.

Despite Shared Geography, Campylobacter Isolated from Surface Water Are Genetically Distinct from Campylobacter Isolated from Chickens

We tested the hypothesis that Campylobacter isolated from chicken ceca and river water in an overlapping geographic area would share genetic information. Isolates of C. jejuni from chicken ceca were collected from a commercial slaughter plant and isolates of C. jejuni were also collected from rivers and creeks in the same watershed. Isolates were subjected to whole-genome sequencing and the data were used for core genome multilocus sequence typing (cgMLST). Cluster analysis showed that there were four distinct subpopulations, two from chickens and two from water. Calculation of fixation statistic (Fst) showed that all four subpopulations were significantly distinct. Greater than 90% of the loci were differentiated by subpopulation. Only two genes showed clear differentiation of both chicken subpopulations from both water subpopulations. Sequence fragments of the CJIE4 bacteriophage family were found frequently in the main chicken subpopulation and the water outgroup subpopulation but were sparsely found in the main water population and not at all in the chicken outgroup. CRISPR spacers that targeted the phage sequences were common in the main water subpopulation, only once in the main chicken subpopulation, and not at all in the chicken or water outgroups. Restriction enzyme genes also showed a biased distribution. These data suggest that there is little transfer of C. jejuni genetic material between chickens and nearby river water. Campylobacter differentiation according to these two sources does not show clear evidence of evolutionary selection; the differentiation is probably due to geospatial isolation, genetic drift, and the action of CRISPRs and restriction enzymes. IMPORTANCE Campylobacter jejuni causes gastroenteritis in humans, and chickens and environmental water are leading sources of infection. We tested the hypothesis that Campylobacter isolated from chicken ceca and river water in an overlapping geographic area would share genetic information. Isolates of Campylobacter were collected from water and chicken sources in the same watershed and their genomes were sequenced and analyzed. Four distinct subpopulations were found. There was no evidence of sharing genetic material between the subpopulations. Phage profiles, CRISPR profiles and restriction systems differed by subpopulation.

Buffered Peptone Water Formulation Does Not Influence Growth of pESI-positive Salmonella enterica Serovar Infantis

Salmonella enterica is a major cause of human foodborne illness and is often attributed to poultry food sources. S. enterica serovar Infantis, specifically those carrying the pESI plasmid, has become a frequently isolated serotype from poultry meat samples at processing and has caused numerous recent human infections. In 2016, the USDA-Food Safety and Inspection Service changed the official sampling method for raw poultry products from BPW to using neutralizing BPW (nBPW) as the rinsing agent in order to prevent residual antimicrobial effects from acidifying and oxidizing processing aids. This change was contemporaneous to the emergence of pESI-positive ser. Infantis as a prevalent serovar in poultry, prompting some to question if nBPW could be selecting for this prevalent serovar. We performed two experiments: a comparison of ser. Infantis growth in BPW versus nBPW, and a simulation of regulatory sampling methods. We found that when inoculated into both broths, ser. Infantis initially grows slightly slower in nBPW than in BPW but little difference was seen in abundance after 6 h of growth. Additionally, the use of nBPW to simulate poultry rinse sample and overnight cold shipping to a regulatory lab did not affect the survival or subsequent growth of ser. Infantis in BPW. We concluded that the change in USDA-FSIS methodology to include nBPW in sampling procedures has likely not affected the emergence of S. ser. Infantis as a prevalent serovar in chicken and turkey meat product samples.

Analyzing the genetic diversity and biotechnological potential of Leuconostoc pseudomesenteroides by comparative genomics

Leuconostoc pseudomesenteroides is a lactic acid bacteria species widely exist in fermented dairy foods, cane juice, sourdough, kimchi, apple dumpster, caecum, and human adenoid. In the dairy industry, Ln. pseudomesenteroides strains are usually found in mesophilic starter cultures with lactococci. This species plays a crucial role in the production of aroma compounds such as acetoin, acetaldehyde, and diacetyl, thus beneficially affecting dairy technology. We performed genomic characterization of 38 Ln. pseudomesenteroides from diverse ecological niches to evaluate this species' genetic diversity and biotechnological potential. A mere ~12% of genes conserved across 38 Ln. pseudomesenteroides genomes indicate that accessory genes are the driving force for genotypic distinction in this species. Seven main clades were formed with variable content surrounding mobile genetic elements, namely plasmids, transposable elements, IS elements, prophages, and CRISPR-Cas. All but three genomes carried CRISPR-Cas system. Furthermore, a type IIA CRISPR-Cas system was found in 80% of the CRISPR-Cas positive strains. AMBR10, CBA3630, and MGBC116435 were predicted to encode bacteriocins. Genes responsible for citrate metabolism were found in all but five strains belonging to cane juice, sourdough, and unknown origin. On the contrary, arabinose metabolism genes were only available in nine strains isolated from plant-related systems. We found that Ln. pseudomesenteroides genomes show evolutionary adaptation to their ecological environment due to niche-specific carbon metabolism and forming closely related phylogenetic clades based on their isolation source. This species was found to be a reservoir of type IIA CRISPR-Cas system. The outcomes of this study provide a framework for uncovering the biotechnological potential of Ln. pseudomesenteroides and its future development as starter or adjunct culture for dairy industry.

Comparative Genomics of Lentilactobacillus parabuchneri isolated from dairy, KEM complex, Makgeolli, and Saliva Microbiomes

BACKGROUND

Lentilactobacillus parabuchneri is of particular concern in fermented food bioprocessing due to causing unwanted gas formation, cracks, and off-flavor in fermented dairy foods. This species is also a known culprit of histamine poisonings because of decarboxylating histidine to histamine in ripening cheese. Twenty-eight genomes in NCBI GenBank were evaluated via comparative analysis to determine genomic diversity within this species and identify potential avenues for reducing health associated risks and economic losses in the food industry caused by these organisms.

RESULT

Core genome-based phylogenetic analysis revealed four distinct major clades. Eight dairy isolates, two strains from an unknown source, and a saliva isolate formed the first clade. Three out of five strains clustered on clade 2 belonged to dairy, and the remaining two strains were isolated from the makgeolli and Korean effective microorganisms (KEM) complex. The third and fourth clade members were isolated from Tete de Moine and dairy-associated niches, respectively. Whole genome analysis on twenty-eight genomes showed ~ 40% of all CDS were conserved across entire strains proposing a considerable diversity among L. parabuchneri strains analyzed. After assigning CDS to their corresponding function, ~ 79% of all strains were predicted to carry putative intact prophages, and ~ 43% of the strains harbored at least one plasmid; however, all the strains were predicted to encode genomic island, insertion sequence, and CRISPR-Cas system. A type I-E CRISPR-Cas subgroup was identified in all the strains, with the exception of DSM15352, which carried a type II-A CRISPR-Cas system. Twenty strains were predicted to encode histidine decarboxylase gene cluster that belongs to not only dairy but also saliva, KEM complex, and unknown source. No bacteriocin-encoding gene(s) or antibiotic resistome was found in any of the L. parabuchneri strains screened.

CONCLUSION

The findings of the present work provide in-depth knowledge of the genomics of L. parabuchneri by comparing twenty-eight genomes available to date. For example, the hdc gene cluster was generally reported in cheese isolates; however, our findings in the current work indicated that it could also be encoded in those strains isolated from saliva, KEM complex, and unknown source. We think prophages are critical mobile elements of L. parabuchneri genomes that could pave the way for developing novel tools to reduce the occurrence of this unwanted species in the food industry.

Comparative genomics of Leuconostoc lactis strains isolated from human gastrointestinal system and fermented foods microbiomes

Background

Leuconostoc lactis forms a crucial member of the genus Leuconostoc and has been widely used in the fermentation industry to convert raw material into acidified and flavored products in dairy and plant-based food systems. Since the ecological niches that strains of Ln. lactis being isolated from were truly diverse such as the human gut, dairy, and plant environments, comparative genome analysis studies are needed to better understand the strain differences from a metabolic adaptation point of view across diverse sources of origin. We compared eight Ln. lactis strains of 1.2.28, aa_0143, BIOML-A1, CBA3625, LN19, LN24, WIKIM21, and WiKim40 using bioinformatics to elucidate genomic level characteristics of each strain for better utilization of this species in a broad range of applications in food industry.

Results

Phylogenomic analysis of twenty-nine Ln. lactis strains resulted in nine clades. Whole-genome sequence analysis was performed on eight Ln. lactis strains representing human gastrointestinal tract and fermented foods microbiomes. The findings of the present study are based on comparative genome analysis against the reference Ln. lactis CBA3625 genome. Overall, a ~ 41% of all CDS were conserved between all strains. When the coding sequences were assigned to a function, mobile genetic elements, mainly insertion sequences were carried by all eight strains. All strains except LN24 and WiKim40 harbor at least one intact putative prophage region, and two of the strains contained CRISPR-Cas system. All strains encoded Lactococcin 972 bacteriocin biosynthesis gene clusters except for CBA3625.

Conclusions

The findings in the present study put forth new perspectives on genomics of Ln. lactis via complete genome sequence based comparative analysis and further determination of genomic characteristics. The outcomes of this work could potentially pave the way for developing elements for future strain engineering applications.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-022-01074-6.

Clustered Regularly Interspaced Short Palindromic Repeats Genotyping of Multidrug-Resistant Salmonella Heidelberg Strains Isolated From the Poultry Production Chain Across Brazil

Salmonella enterica subsp. enterica serovar Heidelberg has been associated with a broad host range, such as poultry, dairy calves, swine, wild birds, environment, and humans. The continuous evolution of S. Heidelberg raises a public health concern since there is a global dispersal of lineages harboring a wide resistome and virulome on a global scale. Here, we characterized the resistome, phylogenetic structure and clustered regularly interspaced short palindromic repeats (CRISPR) array composition of 81 S. Heidelberg strains isolated from broiler farms (n = 16), transport and lairage (n = 5), slaughterhouse (n = 22), and retail market (n = 38) of the poultry production chain in Brazil, between 2015 and 2016 using high-resolution approaches including whole-genome sequencing (WGS) and WGS-derived CRISPR genotyping. More than 91% of the S. Heidelberg strains were multidrug-resistant. The total antimicrobial resistance (AMR) gene abundances did not vary significantly across regions and sources suggesting the widespread distribution of antibiotic-resistant strains from farm to market. The highest AMR gene abundance was observed for fosA7, aac(6′)-Iaa, sul2, tet(A), gyrA, and parC for 100% of the isolates, followed by 88.8% for bla_CMY–2. The β-lactam resistance was essentially driven by the presence of the plasmid-mediated AmpC (pAmpC) bla_CMY–2 gene, given the isolates which did not carry this gene were susceptible to cefoxitin (FOX). Most S. Heidelberg strains were classified within international lineages, which were phylogenetically nested with Salmonella strains from European countries; while CRISPR genotyping analysis revealed that the spacer content was overall highly conserved, but distributed into 13 distinct groups. In summary, our findings underscore the potential role of S. Heidelberg as a key pathogen disseminated from farm to fork in Brazil and reinforce the importance of CRISPR-based genotyping for salmonellae. Hence, we emphasized the need for continuous mitigation programs to monitor the dissemination of this high-priority pathogen.

CrisprVi: a software for visualizing and analyzing CRISPR sequences of prokaryotes

Background

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.

Results

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.

Conclusions

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04716-9.

Evidence for common ancestry and microevolution of passerine-adapted Salmonella enterica serovar Typhimurium in the UK and USA

The evolution of Salmonella enterica serovar Typhimurium (S. Typhimurium) within passerines has resulted in pathoadaptation of this serovar to the avian host in Europe. Recently, we identified an S. Typhimurium lineage from passerines in North America. The emergence of passerine-adapted S. Typhimurium in Europe and North America raises questions regarding its evolutionary origin. Here, we demonstrated that the UK and US passerine-adapted S. Typhimurium shared a common ancestor from ca. 1838, and larids played a key role in the clonal expansion by disseminating the common ancestor between North America and Europe. Further, we identified virulence gene signatures common in the passerine- and larid-adapted S. Typhimurium, including conserved pseudogenes in fimbrial gene lpfD and Type 3 Secretion System (T3SS) effector gene steC. However, the UK and US passerine-adapted S. Typhimurium also possessed unique virulence gene signatures (i.e. pseudogenes in fimbrial gene fimC and T3SS effector genes sspH2, gogB, sseJ and sseK2), and the majority of them (38/47) lost a virulence plasmid pSLT that was present in the larid-adapted S. Typhimurium. These results provide evidence that passerine-adapted S. Typhimurium share a common ancestor with those from larids, and the divergence of passerine- and larid-adapted S. Typhimurium might be due to pseudogenization or loss of specific virulence genes.

Salmonella enterica serovar Typhimurium from Wild Birds in the United States Represent Distinct Lineages Defined by Bird Type

Salmonella enterica serovar Typhimurium is typically considered a host generalist; however, certain isolates are associated with specific hosts and show genetic features of host adaptation. Here, we sequenced 131 S. Typhimurium isolates from wild birds collected in 30 U.S. states during 1978–2019. We found that isolates from broad taxonomic host groups including passerine birds, water birds (Aequornithes), and larids (gulls and terns) represented three distinct lineages and certain S. Typhimurium CRISPR types presented in individual lineages. We also showed that lineages formed by wild bird isolates differed from most isolates originating from domestic animal sources, and that genomes from these lineages substantially improved source attribution of Typhimurium genomes to wild birds by a machine learning classifier. Furthermore, virulence gene signatures that differentiated S. Typhimurium from passerines, water birds, and larids were detected. Passerine isolates tended to lack S. Typhimurium-specific virulence plasmids. Isolates from the passerine, water bird, and larid lineages had close genetic relatedness with human clinical isolates, including those from a 2021 U.S. outbreak linked to passerine birds. These observations indicate that S. Typhimurium from wild birds in the United States are likely host-adapted, and the representative genomic data set examined in this study can improve source prediction and facilitate outbreak investigation.

IMPORTANCE Within-host evolution of S. Typhimurium may lead to pathovars adapted to specific hosts. Here, we report the emergence of disparate avian S. Typhimurium lineages with distinct virulence gene signatures. The findings highlight the importance of wild birds as a reservoir for S. Typhimurium and contribute to our understanding of the genetic diversity of S. Typhimurium from wild birds. Our study indicates that S. Typhimurium may have undergone adaptive evolution within wild birds in the United States. The representative S. Typhimurium genomes from wild birds, together with the virulence gene signatures identified in these bird isolates, are valuable for S. Typhimurium source attribution and epidemiological surveillance.

Comparative Genomic Analyses and CRISPR-Cas Characterization of Cutibacterium acnes Provide Insights Into Genetic Diversity and Typing Applications

Cutibacterium acnes is an important member of the human skin microbiome and plays a critical role in skin health and disease. C. acnes encompasses different phylotypes that have been found to be associated with different skin phenotypes, suggesting a genetic basis for their impact on skin health. Here, we present a comprehensive comparative analysis of 255 C. acnes genomes to provide insights into the species genetic diversity and identify unique features that define various phylotypes. Results revealed a relatively small and open pan genome (6,240 genes) with a large core genome (1,194 genes), and three distinct phylogenetic clades, with multiple robust sub-clades. Furthermore, we identified several unique gene families driving differences between distinct C. acnes clades. Carbohydrate transporters, stress response mechanisms and potential virulence factors, potentially involved in competitive growth and host colonization, were detected in type I strains, which are presumably responsible for acne. Diverse type I-E CRISPR-Cas systems and prophage sequences were detected in select clades, providing insights into strain divergence and adaptive differentiation. Collectively, these results enable to elucidate the fundamental differences among C. acnes phylotypes, characterize genetic elements that potentially contribute to type I-associated dominance and disease, and other key factors that drive the differentiation among clades and sub-clades. These results enable the use of comparative genomics analyses as a robust method to differentiate among the C. acnes genotypes present in the skin microbiome, opening new avenues for the development of biotherapeutics to manipulate the skin microbiota.

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.

Antimicrobial Resistance and CRISPR Typing Among Salmonella Isolates From Poultry Farms in China

Although knowledge of the clustered regularly interspaced short palindromic repeat (CRISPR)-Cas system has been applied in many research areas, comprehensive studies of this system in Salmonella, particularly in analysis of antibiotic resistance, have not been reported. In this work, 75 Salmonella isolates obtained from broilers or broilers products were characterized to determine their antimicrobial susceptibilities, antibiotic resistance gene profiles, and CRISPR array diversities, and genotyping was explored. In total, 80.00% (60/75) of the strains were multidrug resistant, and the main pattern observed in the isolates was CN-AZM-AMP-AMC-CAZ-CIP-ATM-TE-SXT-FOS-C. The resistance genes of streptomycin (aadA), phenicol (floR-like and catB3-like), β-lactams (bla _TEM, bla _OXA, and bla _CTX), tetracycline [tet(A)-like], and sulfonamides (sul1 and sul2) appeared at higher frequencies among the corresponding resistant isolates. Subsequently, we analyzed the CRISPR arrays and found 517 unique spacer sequences and 31 unique direct repeat sequences. Based on the CRISPR spacer sequences, we developed a novel typing method, CRISPR locus three spacer sequences typing (CLTSST), to help identify sources of Salmonella outbreaks especially correlated with epidemiological data. Compared with multi-locus sequence typing (MLST), conventional CRISPR typing (CCT), and CRISPR locus spacer pair typing (CLSPT), discrimination using CLTSST was weaker than that using CCT but stronger than that using MLST and CLSPT. In addition, we also found that there were no close correlations between CRISPR loci and antibiotics but had close correlations between CRISPR loci and antibiotic resistance genes in Salmonella isolates.

CRISPRclassify: Repeat-Based Classification of CRISPR Loci

Detection and classification of CRISPR-Cas systems in metagenomic data have become increasingly prevalent in recent years due to their potential for diverse applications in genome editing. Traditionally, CRISPR-Cas systems are classified through reference-based identification of proximate cas genes. Here, we present a machine learning approach for the detection and classification of CRISPR loci using repeat sequences in a cas-independent context, enabling identification of unclassified loci missed by traditional cas-based approaches. Using biological attributes of the CRISPR repeat, the core element in CRISPR arrays, and leveraging methods from natural language processing, we developed a machine learning model capable of accurate classification of CRISPR loci in an extensive set of metagenomes, resulting in an F1 measure of 0.82 across all predictions and an F1 measure of 0.97 when limiting to classifications with probabilities >0.85. Furthermore, assessing performance on novel repeats yielded an F1 measure of 0.96. Although the performance of cas-based identification will exceed that of a repeat-based approach in many cases, CRISPRclassify provides an efficient approach to classification of CRISPR loci for cases in which cas gene information is unavailable, such as metagenomes and fragmented genome assemblies.

Comparative genomic analysis of Mycoplasma anatis strains

Background

The Gram-negative intracellular bacterium Mycoplasma anatis is a pathogen of respiratory infectious diseases in ducks and has caused significant economic losses in the poultry industry.

Objective

This study, as the first report of the structure and function of the pan-genome of Mycoplasma anatis, may provide a valuable genetic basis for many aspects of future research on the pathogens of waterfowl.

Methods

We sequenced the whole genomes of 15 Mycoplasma anatis isolated from ducks in China. Draft genome sequencing was carried out and whole-genome sequencing was performed by the sequencers of the PacBio Sequel and an IonTorrent Personal Genome Machine (PGM). Then the common genic elements of protein-coding genes, tRNAs, and rRNAs of Mycoplasma anatis genomes were predicted by using the pipeline Prokka v1.13.7. To investigate homologous protein clusters across Mycoplasma anatis genomes, we adopted Roary v3.13.0 to cluster orthologous genes (OGs) based on the following criteria.

Results

We obtained one complete genome and 14 genome sketches. Microbial mobile genetic element analysis revealed the distribution of insertion sequences (IS30, IS3, and IS1634), prophage regions, and CRISPR arrays in the genome of Mycoplasma anatis. Comparative genomic analysis decoded the genetic components and functional classification of the pan-genome of Mycoplasma anatis that comprised 646 core genes, 231 dispensable genes and among them 110 was strain-specific. Virulence-related gene profiles of Mycoplasma anatis were systematically identified, and the products of these genes included bacterial ABC transporter systems, iron transport proteins, toxins, and secretion systems.

Conclusion

A complete virulence-related gene profile of Mycoplasma anatis has been identified, most of the genes are highly conserved in all strains. Sequencing results are relevant to the molecular mechanisms of drug resistance, adaptive evolution of pathogens, population structure, and vaccine development.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13258-021-01129-5.

Diversity of the type I-U CRISPR-Cas system in Bifidobacterium

The CRISPR-Cas system is widely distributed in prokaryotes and plays an important role in the adaptive immunity of bacteria and archaea. Bifidobacterium is an important component of the intestinal flora of humans and animals, and some species of this bacterium can be employed as food additives. However, the Bifidobacterium CRISPR-Cas system has not been fully elucidated to date. In this study, the genomes of 110 strains of Bifidobacterium were employed to research the diversity of the type I-U system. The 110 strains were divided into five groups according to the genes adjacent to the CRISPR locus, including group A, B, C, D and E. Strains in the intergroup had unique species classifications and MLST types. An evolutionary tree was constructed based on the conserved cas4/cas1 fusion gene. The results showed that group A had a different evolutionary branch compared with the other groups and had a relatively low spacer number. Notably, group B, C and E had exhibited ABC transporter regulators in the genes adjacent to the CRISPR locus. ABC transporters play important roles in the exocytosis of many antibiotics and are involved in horizontal gene transfer. This mechanism may have promoted the evolution of Bifidobacterium and the horizontal gene transfer of the type I-U system, which may have promoted the generation of system diversity. In summary, our results help to elucidate the role of the type I-U system in the evolution of Bifidobacterium.

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.

Whole-genome sequencing analysis and CRISPR genotyping of rare antibiotic-resistant Salmonella enterica serovars isolated from food and related sources

For decades, Salmonella Typhimurium and Salmonella Enteritidis have prevailed in several countries as agents of salmonellosis outbreaks. In Brazil, the largest exporter of poultry meat, relatively little attention has been paid to infrequent serovars. Here, we report the emergence and characterization of rare serovars isolated from food and related sources collected between 2014 and 2016 in Brazil. Twenty-two Salmonella enterica isolates were analyzed through the use of whole-genome sequencing (WGS) and clustered regularly interspaced short palindromic repeats (CRISPR) genotyping. These isolates were classified into 10 infrequent serovars, including S. Abony, S. Isangi, S. Rochdale, S. Saphra, S. Orion, S. Ouakam, S. Grumpensis, S. Carrau, S. Abaetetuba, and S. Idikan. The presence of six antimicrobial resistance (AMR) genes, qnrB19, bla_CMY-2, tetA, aac(6')-Iaa, sul2 and fosA7, which encode resistance to quinolones, third-generation cephalosporin, tetracycline, aminoglycoside, sulfonamide and fosfomycin, respectively, were confirmed by WGS. All S. Isangi harbored qnrB19 with conserved genomic context across strains, while S. Abony harbored bla_CMY-2. Twelve (54.5%) strains displayed chromosomal mutations in parC (Thr57→Ser). Most serovars were classified as independent lineages, except S. Abony and S. Abaetetuba, which phylogenetically nested with Salmonella strains from different countries. CRISPR analysis revealed that the spacer content was strongly correlated with serovar and multi-locus sequence type for all strains, independently confirming the observed phylogenetic patterns, and highlighting the value of CRISPR-based genotyping for Salmonella. These findings add valuable information to the epidemiology of S. enterica in Brazil, where the emergency of antibiotic-resistant Salmonella continues to evolve.

Harnessing the type I CRISPR-Cas systems for genome editing in prokaryotes

Genetic analysis is crucial to the understanding, exploitation, and control of microorganisms. The advent of CRISPR-Cas-based genome-editing techniques, particularly those mediated by the single-effector (Cas9 and Cas12a) class 2 CRISPR-Cas systems, has revolutionized the genetics in model eukaryotic organisms. However, their applications in prokaryotes are rather limited, largely owing to the exceptional diversity of DNA homeostasis in microorganisms and severe cytotoxicity of overexpressing these nuclease proteins in certain genotypes. Remarkably, CRISPR-Cas systems belonging to different classes and types are continuously identified in prokaryotic genomes and serve as a deep reservoir for expansion of the CRISPR-based genetic toolkits. ~90% of the CRISPR-Cas systems identified so far belong to the class 1 system which hinges on multi-protein effector complexes for DNA interference. Harnessing these widespread native CRISPR-Cas systems for 'built-in' genome editing represents an emerging and powerful genetic tool in prokaryotes, especially in the genetically recalcitrant non-model species and strains. In this progress review, we introduce the general workflow of this emerging editing platform and summarize its establishment in a growing number of prokaryotes by harnessing the most widespread, diverse type I CRISPR-Cas systems present in their genomes. We also discuss the various factors affecting the success and efficiency of this editing platform and the corresponding solutions.

Applications of CRISPR-Cas systems in lactic acid bacteria

As a phenotypically and phylogenetically diverse group, lactic acid bacteria are found in a variety of natural environments and occupy important roles in medicine, biotechnology, food and agriculture. The widespread use of lactic acid bacteria across these industries fuels the need for new and functionally diverse strains that may be utilized as starter cultures or probiotics. Originally characterized in lactic acid bacteria, CRISPR-Cas systems and derived molecular machines can be used natively or exogenously to engineer new strains with enhanced functional attributes. Research on CRISPR-Cas biology and its applications has exploded over the past decade with studies spanning from the initial characterization of CRISPR-Cas immunity in Streptococcus thermophilus to the use of CRISPR-Cas for clinical gene therapies. Here, we discuss CRISPR-Cas classification, overview CRISPR biology and mechanism of action, and discuss current and future applications in lactic acid bacteria, opening new avenues for their industrial exploitation and manipulation of microbiomes.

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

In Archaea and Bacteria, the arrays called CRISPRs for 'clustered regularly interspaced short palindromic repeats' and the CRISPR associated genes or cas provide adaptive immunity against viruses, plasmids and transposable elements. Short sequences called spacers, corresponding to fragments of invading DNA, are stored in-between repeated sequences. The CRISPR-Cas systems target sequences homologous to spacers leading to their degradation. To facilitate investigations of CRISPRs, we developed 12 years ago a website holding the CRISPRdb. We now propose CRISPRCasdb, a completely new version giving access to both CRISPRs and cas genes. We used CRISPRCasFinder, a program that identifies CRISPR arrays and cas genes and determine the system's type and subtype, to process public whole genome assemblies. Strains are displayed either in an alphabetic list or in taxonomic order. The database is part of the CRISPR-Cas++ website which also offers the possibility to analyse submitted sequences and to download programs. A BLAST search against lists of repeats and spacers extracted from the database is proposed. To date, 16 990 complete prokaryote genomes (16 650 bacteria from 2973 species and 340 archaea from 300 species) are included. CRISPR-Cas systems were found in 36% of Bacteria and 75% of Archaea strains. CRISPRCasdb is freely accessible at https://crisprcas.i2bc.paris-saclay.fr/.

Comprehensive Mining and Characterization of CRISPR-Cas Systems in Bifidobacterium

The clustered regularly interspaced short palindromic repeats (CRISPR)-Cas (CRISPR-associated cas) systems constitute the adaptive immune system in prokaryotes, which provides resistance against bacteriophages and invasive genetic elements. The landscape of applications in bacteria and eukaryotes relies on a few Cas effector proteins that have been characterized in detail. However, there is a lack of comprehensive studies on naturally occurring CRISPR-Cas systems in beneficial bacteria, such as human gut commensal Bifidobacterium species. In this study, we mined 954 publicly available Bifidobacterium genomes and identified CRIPSR-Cas systems in 57% of these strains. A total of five CRISPR-Cas subtypes were identified as follows: Type I-E, I-C, I-G, II-A, and II-C. Among the subtypes, Type I-C was the most abundant (23%). We further characterized the CRISPR RNA (crRNA), tracrRNA, and PAM sequences to provide a molecular basis for the development of new genome editing tools for a variety of applications. Moreover, we investigated the evolutionary history of certain Bifidobacterium strains through visualization of acquired spacer sequences and demonstrated how these hypervariable CRISPR regions can be used as genotyping markers. This extensive characterization will enable the repurposing of endogenous CRISPR-Cas systems in Bifidobacteria for genome engineering, transcriptional regulation, genotyping, and screening of rare variants.

Genomic characterization of Lactobacillus fermentum DSM 20052

Background

Lactobacillus fermentum, a member of the lactic acid bacteria complex, has recently garnered increased attention due to documented antagonistic properties and interest in assessing the probiotic potential of select strains that may provide human health benefits. Here, we genomically characterize L. fermentum using the type strain DSM 20052 as a canonical representative of this species.

Results

We determined the polished whole genome sequence of this type strain and compared it to 37 available genome sequences within this species. Results reveal genetic diversity across nine clades, with variable content encompassing mobile genetic elements, CRISPR-Cas immune systems and genomic islands, as well as numerous genome rearrangements. Interestingly, we determined a high frequency of occurrence of diverse Type I, II, and III CRISPR-Cas systems in 72% of the genomes, with a high level of strain hypervariability.

Conclusions

These findings provide a basis for the genetic characterization of L. fermentum strains of scientific and commercial interest. Furthermore, our study enables genomic-informed selection of strains with specific traits for commercial product formulation, and establishes a framework for the functional characterization of features of interest.

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.

Prevalence, genetic analysis and CRISPR typing of Cronobacter spp. isolated from meat and meat products in China

Cronobacter spp. are important foodborne pathogens that infections occur in all age groups, especially cause serious life-threatening diseases in infants. This study aimed to acquire data on Cronobacter spp. contamination of meat and meat products (n = 588) in China during 2011 to 2016, and investigated the use of CRISPR typing technology as an approach for characterizing the genetics of Cronobacter spp. The overall contamination rate for Cronobacter spp. was determined to be 9.18% (54/588). Of the positive samples, 90.74% (49/54) had <10 MPN/g, with duck samples had a relatively high contamination rate (15.69%, 8/51) and highest contamination level (28.90 MPN/g). Four species and nine serotypes were identified among 69 isolates, of which C. sakazakii was the major species (n = 50) and C. sakazakii serogroup O1 and O2 (n = 17) were the primary serotypes. The majority of Cronobacter spp. strains were found to be susceptible to most antibiotics except exhibited high resistance to cephalothin (76.81%, 53/69), and total two multi-drug resistant C. sakazakii strains were isolated from duck. The genetic diversity of Cronobacter spp. was remarkably high, as evidenced by the identification of 40 sequence types (STs) and 60 CRISPR types (CTs). C. sakazakii ST64 (n = 7) was the predominant genotype and was further divided into two sub-lineages based on CRISPR diversity, showing different antibiotic resistance profile. These results demonstrate that CRISPR typing results have a good correspondence with bacterial phenotypes, and it will be a tremendously useful approach for elucidating inter-subtyping during molecular epidemiological investigations while interpreting the divergent evolution of Cronobacter. The presence of Cronobacter spp. in meat and meat product is a potential threat to human public health.

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

In Archaea and Bacteria, the arrays called CRISPRs for 'clustered regularly interspaced short palindromic repeats' and the CRISPR associated genes or cas provide adaptive immunity against viruses, plasmids and transposable elements. Short sequences called spacers, corresponding to fragments of invading DNA, are stored in-between repeated sequences. The CRISPR-Cas systems target sequences homologous to spacers leading to their degradation. To facilitate investigations of CRISPRs, we developed 12 years ago a website holding the CRISPRdb. We now propose CRISPRCasdb, a completely new version giving access to both CRISPRs and cas genes. We used CRISPRCasFinder, a program that identifies CRISPR arrays and cas genes and determine the system's type and subtype, to process public whole genome assemblies. Strains are displayed either in an alphabetic list or in taxonomic order. The database is part of the CRISPR-Cas++ website which also offers the possibility to analyse submitted sequences and to download programs. A BLAST search against lists of repeats and spacers extracted from the database is proposed. To date, 16 990 complete prokaryote genomes (16 650 bacteria from 2973 species and 340 archaea from 300 species) are included. CRISPR-Cas systems were found in 36% of Bacteria and 75% of Archaea strains. CRISPRCasdb is freely accessible at https://crisprcas.i2bc.paris-saclay.fr/.

Comparative genomics of eight Lactobacillus buchneri strains isolated from food spoilage

Abstract

Background

Lactobacillus buchneri is a lactic acid bacterium frequently associated with food bioprocessing and fermentation and has been found to be either beneficial or detrimental to industrial food processes depending on the application. The ability to metabolize lactic acid into acetic acid and 1,2-propandiol makes L. buchneri invaluable to the ensiling process, however, this metabolic activity leads to spoilage in other applications, and is especially damaging to the cucumber fermentation industry. This study aims to augment our genomic understanding of L. buchneri in order to make better use of the species in a wide range of applicable industrial settings.

Results

Whole-genome sequencing (WGS) was performed on seven phenotypically diverse strains isolated from spoiled, fermented cucumber and the ATCC type strain for L. buchneri, ATCC 4005. Here, we present our findings from the comparison of eight newly-sequenced and assembled genomes against two publicly available closed reference genomes, L. buchneri CD034 and NRRL B-30929. Overall, we see ~ 50% of all coding sequences are conserved across these ten strains. When these coding sequences are clustered by functional description, the strains appear to be enriched in mobile genetic elements, namely transposons. All isolates harbor at least one CRISPR-Cas system, and many contain putative prophage regions, some of which are targeted by the host’s own DNA-encoded spacer sequences.

Conclusions

Our findings provide new insights into the genomics of L. buchneri through whole genome sequencing and subsequent characterization of genomic features, building a platform for future studies and identifying elements for potential strain manipulation or engineering.

Cronobacter sakazakii, Cronobacter malonaticus, and Cronobacter dublinensis Genotyping Based on CRISPR Locus Diversity

Cronobacter strains harboring CRISPR-Cas systems are important foodborne pathogens that cause serious neonatal infections. CRISPR typing is a new molecular subtyping method to track the sources of pathogenic bacterial outbreaks and shows a promise in typing Cronobacter, however, this molecular typing procedure using routine PCR method has not been established. Therefore, the purpose of this study was to establish such methodology, 257 isolates of Cronobacter sakazakii, C. malonaticus, and C. dublinensis were used to verify the feasibility of the method. Results showed that 161 C. sakazakii strains could be divided into 129 CRISPR types (CTs), among which CT15 (n = 7) was the most prevalent CT followed by CT6 (n = 4). Further, 65 C. malonaticus strains were divided into 42 CTs and CT23 (n = 8) was the most prevalent followed by CT2, CT3, and CT13 (n = 4). Finally, 31 C. dublinensis strains belonged to 31 CTs. There was also a relationship among CT, sequence type (ST), food types, and serotype. Compared to multi-locus sequence typing (MLST), this new molecular method has greater power to distinguish similar strains and had better accordance with whole genome sequence typing (WGST). More importantly, some lineages were found to harbor conserved ancestral spacers ahead of their divergent specific spacer sequences; this can be exploited to infer the divergent evolution of Cronobacter and provide phylogenetic information reflecting common origins. Compared to WGST, CRISPR typing method is simpler and more affordable, it could be used to identify sources of Cronobacter food-borne outbreaks, from clinical cases to food sources and the production sites.

Predicting and visualizing features of CRISPR-Cas systems

Pervasive application of CRISPR-Cas systems in genome editing has prompted an increase in both interest and necessity to further elucidate existing systems as well as discover putative novel systems. The ubiquity and power of current computational platforms have made in silico approaches to CRISPR-Cas identification and characterization accessible to a wider audience and increasingly amenable for processing extensive data sets. Here, we describe in silico methods for predicting and visualizing notable features of CRISPR-Cas systems, including Cas domain determination, CRISPR array visualization, and inference of the protospacer-adjacent motif. The efficiency of these tools enables rapid exploration of CRISPR-Cas diversity across prokaryotic genomes and supports scalable analysis of large genomic data sets.

CRISPRStudio: A User-Friendly Software for Rapid CRISPR Array Visualization

The CRISPR-Cas system biologically serves as an adaptive defense mechanism against phages. However, there is growing interest in exploiting the hypervariable nature of the CRISPR locus, often of viral origin, for microbial typing and tracking. Moreover, the spacer content of any given strain provides a phage resistance profile. Large-scale CRISPR typing studies require an efficient method for showcasing CRISPR array similarities across multiple isolates. Historically, CRISPR arrays found in microbes have been represented by colored shapes based on nucleotide sequence identity and, while this approach is now routinely used, only scarce computational resources are available to automate the process, making it very time-consuming for large datasets. To alleviate this tedious task, we introduce CRISPRStudio, a command-line tool developed to accelerate CRISPR analysis and standardize the preparation of CRISPR array figures. It first compares nucleotide spacer sequences present in a dataset and then clusters them based on sequence similarity to assign a meaningful representative color. CRISPRStudio offers versatility to suit different biological contexts by including options such as automatic sorting of CRISPR loci and highlighting of shared spacers, while remaining fast and user-friendly.