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

Exploring the comparative genome of rice pathogen Burkholderia plantarii: unveiling virulence, fitness traits, and a potential type III secretion system effector

This study presents a comprehensive genomic analysis of Burkholderia plantarii, a rice pathogen that causes blight and grain rot in seedlings. The entire genome of B. plantarii KACC 18964 was sequenced, followed by a comparative genomic analysis with other available genomes to gain insights into its virulence, fitness, and interactions with rice. Multiple secondary metabolite gene clusters were identified. Among these, 12 demonstrated varying similarity levels to known clusters linked to bioactive compounds, whereas eight exhibited no similarity, indicating B. plantarii as a source of potentially novel secondary metabolites. Notably, the genes responsible for tropolone and quorum sensing were conserved across the examined genomes. Additionally, B. plantarii was observed to possess three complete CRISPR systems and a range of secretion systems, exhibiting minor variations among the analyzed genomes. Genomic islands were analyzed across the four genomes, and a detailed study of the B. plantarii KACC 18964 genome revealed 59 unique islands. These islands were thoroughly investigated for their gene contents and potential roles in virulence. Particular attention has been devoted to the Type III secretion system (T3SS), a crucial virulence factor. An in silico analysis of potential T3SS effectors identified a conserved gene, aroA. Further mutational studies, in planta and in vitro analyses validated the association between aroA and virulence in rice. Overall, this study enriches our understanding of the genomic basis of B. plantarii pathogenicity and emphasizes the potential role of aroA in virulence. This understanding may guide the development of effective disease management strategies.

Escherichia coli CRISPR arrays from early life fecal samples preferentially target prophages

CRISPR-Cas systems are defense mechanisms against phages and other nucleic acids that invade bacteria and archaea. In Escherichia coli, it is generally accepted that CRISPR-Cas systems are inactive in laboratory conditions due to a transcriptional repressor. In natural isolates, it has been shown that CRISPR arrays remain stable over the years and that most spacer targets (protospacers) remain unknown. Here, we re-examine CRISPR arrays in natural E. coli isolates and investigate viral and bacterial genomes for spacer targets using a bioinformatics approach coupled to a unique biological dataset. We first sequenced the CRISPR1 array of 1769 E. coli isolates from the fecal samples of 639 children obtained during their first year of life. We built a network with edges between isolates that reflect the number of shared spacers. The isolates grouped into 34 modules. A search for matching spacers in bacterial genomes showed that E. coli spacers almost exclusively target prophages. While we found instances of self-targeting spacers, those involving a prophage and a spacer within the same bacterial genome were rare. The extensive search for matching spacers also expanded the library of known E. coli protospacers to 60%. Altogether, these results favor the concept that E. coli's CRISPR-Cas is an antiprophage system and highlight the importance of reconsidering the criteria use to deem CRISPR-Cas systems active.

The CRISPR-Cas system in clinical strains of Acinetobacter baumannii: an in-silico analysis

Acinetobacter baumannii is a relevant bacterium due to its high-resistance profile. It is well known that antimicrobial resistance is primarily linked to mutations and the acquisition of external genomic material, such as plasmids or phages, to which the Clustered Regularly Interspaced Short Palindromic Repeats associated with Cas proteins, or CRISPR-Cas, system is related. It is known that the system can influence the acquisition of foreign genetic material and play a role in various physiological pathways. In this study, we conducted an in-silico analysis using 91 fully assembled genomes of clinical strains obtained from the NCBI database. Among the analyzed genomes, the I-F1 subtype of the CRISPR-Cas system was detected showcasing variations in architecture and phylogeny. Using bioinformatic tools, we determined the presence, distribution, and specific characteristics of the CRISPR-Cas system. We found a possible association of the system with resistance genes but not with virulence determinants. Analysis of the system's components, including spacer sequences, suggests its potential role in protecting against phage infections, highlighting its protective function.

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.

Pangenomic Study of Fusobacterium nucleatum Reveals the Distribution of Pathogenic Genes and Functional Clusters at the Subspecies and Strain Levels

Fusobacterium nucleatum is a prevalent periodontal pathogen and is associated with many systemic diseases. Our knowledge of the genomic characteristics and pathogenic effectors of different F. nucleatum strains is limited. In this study, we completed the whole genome assembly of the 4 F. nucleatum strains and carried out a comprehensive pangenomic study of 30 strains with their complete genome sequences. Phylogenetic analysis revealed that the F. nucleatum strains are mainly divided into 4 subspecies, while 1 of the sequenced strains was classified into a new subspecies. Gene composition analysis revealed that a total of 517 "core/soft-core genes" with housekeeping functions widely distributed in almost all the strains. Each subspecies had a unique gene cluster shared by strains within the subspecies. Analysis of the virulence factors revealed that many virulence factors were widely distributed across all the strains, with some present in multiple copies. Some virulence genes showed no consistent occurrence rule at the subspecies level and were specifically distributed in certain strains. The genomic islands mainly revealed strain-specific characteristics instead of subspecies level consistency, while CRISPR types and secondary metabolite biosynthetic gene clusters were identically distributed in F. nucleatum strains from the same subspecies. The variation in amino acid sites in the adhesion protein FadA did not affect the monomer and dimer 3D structures, but it may affect the binding surface and the stability of binding to host receptors. This study provides a basis for the pathogenic study of F. nucleatum at the subspecies and strain levels. IMPORTANCE We used F. nucleatum as an example to analyze the genomic characteristics of oral pathogens at the species, subspecies, and strain levels and elucidate the similarities and differences in functional genes and virulence factors among different subspecies/strains of the same oral pathogen. We believe that the unique biological characteristics of each subspecies/strain can be attributed to the differences in functional gene clusters or the presence/absence of certain virulence genes. This study showed that F. nucleatum strains from the same subspecies had similar functional gene compositions, CRISPR types, and secondary metabolite biosynthetic gene clusters, while pathogenic genes, such as virulence genes, antibiotic resistance genes, and GIs, had more strain level specificity. The findings of this study suggest that, for microbial pathogenicity studies, we should carefully consider the subspecies/strains being used, as different strains may vary greatly.

CRISPR-Cas phage defense systems and prophages in Candidatus Accumulibacter

Candidatus Accumulibacter plays a major role in enhanced biological phosphorus removal (EBPR) from wastewater. Although bacteriophages have been shown to represent fatal threats to Ca. Accumulibacter organisms and thus interfere with the stability of the EBPR process, little is known about the ability of different Ca. Accumulibacter strains to resist phage infections. We conducted a systematic analysis of the occurrence and characteristics of clustered regularly interspaced short palindromic repeats and associated proteins (CRISPR-Cas) systems and prophages in Ca. Accumulibacter lineage members (43 in total, including 10 newly recovered genomes). Results indicate that 28 Ca. Accumulibacter genomes encode CRISPR-Cas systems. They were likely acquired via horizontal gene transfer, conveying a distinct adaptivity to phage predation to different Ca. Accumulibacter members. Major differences in the number of spacers show the unique phage resistance of these members. A comparison of the spacers in closely related Ca. Accumulibacter members from distinct geographical locations indicates that habitat isolation may have resulted in the acquisition of resistance to different phages by different Ca. Accumulibacter. Long-term operation of three laboratory-scale EBPR bioreactors revealed high relative abundances of Ca. Accumulibacter with CRISPSR-Cas systems. Their specific resistance to phages in these reactors was indicated by spacer analysis. Metatranscriptomic analyses showed the activation of the CRISPR-Cas system under both anaerobic and aerobic conditions. Additionally, 133 prophage regions were identified in 43 Ca. Accumulibacter genomes. Twenty-seven of them (in 19 genomes) were potentially active. Major differences in the occurrence of CRISPR-Cas systems and prophages in Ca. Accumulibacter will lead to distinct responses to phage predation. This study represents the first systematic analysis of CRISPR-Cas systems and prophages in the Ca. Accumulibacter lineage, providing new perspectives on the potential impacts of phages on Ca. Accumulibacter and EBPR systems.

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.

Dairy streptococcal cell wall and exopolysaccharide genome diversity

The large-scale and high-intensity application of Streptococcus thermophilus species in milk fermentation processes is associated with a persistent threat of (bacterio)phage infection. Phage infection of starter cultures may cause inconsistent, slow or even failed fermentations with consequent diminished product quality and/or output. The phage life cycle commences with the recognition of, and binding to, a specific host-encoded and surface-exposed receptor, which in the case of S. thermophilus can be the rhamnose-glucose polysaccharide (RGP; specified by the rgp gene cluster) or exopolysaccharide (EPS; specified by the eps gene cluster). The genomic diversity of 23 S . thermophilus strains isolated from unpasteurized dairy products was evaluated, including a detailed analysis of the rgp and eps loci. In the present study, five novel eps genotypes were identified while variations of currently recognized rgp gene cluster types were also observed. Furthermore, the diversity of rgp genotypes amongst retrieved isolates positively correlated with phage diversity based on phageome analysis of eight representative dairy products. Our findings therefore substantially expand our knowledge on S. thermophilus’ strain and phage diversity in (artisanal) dairy products and highlight the merit of phageome analysis of artisanal and traditional fermented foods as a sensitive marker of dominant microbiota involved in the fermentation.

Large-Scale Phylogenetic Analysis Reveals a New Genetic Clade among Escherichia coli O26 Strains

Shiga toxin-producing Escherichia coli (STEC) O26 is the predominant non-O157 serogroup causing hemolytic uremic syndrome worldwide. Moreover, the serogroup is highly dynamic and harbors several pathogenic clones. Here, we investigated the phylogenetic relationship of STEC O26 at a global level based on 1,367 strains from 20 countries deposited in NCBI and Enterobase databases. The whole-genome-based analysis identified a new genetic clade, called ST29C4. The new clade was unique in terms of multilocus sequence type (ST29), CRISPR (group Ia), and dominant plasmid gene profile (ehxA+/katP-/espP-/etpD-). Moreover, the combination of multiple typing methods (core genome single nucleotide polymorphism [SNP] typing, CRISPR typing, and virulence genes analysis) demonstrated that this new lineage ST29C4 was in the intermediate phylogenetic position between ST29C3 and other non-ST29C3 strains. Besides, we observed that ST29C4 harbored extraintestinal pathogenic E. coli (ExPEC)-related virulence gene (VG), tsh, and STEC-associated VG, stx2a, suggesting the emergence of a hybrid pathogen. The ST29C4 strains also exhibited high similarity in stx2a-prophage and integrase with the O104:H4 strain, further demonstrating its potential risk to human health. Collectively, the large-scale phylogenetic analysis extends the understanding of the clonal structure of O26 strains and provides new insights for O26 strain microevolution.

IMPORTANCE Shiga toxin-producing Escherichia coli (STEC) O26 is the second prevalent STEC serogroup only to O157, which can cause a series of diseases ranging from mild diarrhea to life-threatening hemolytic uremic syndrome (HUS). The serogroup is highly diverse and multiple clones are characterized, including ST29C1-C3 and ST21C1-C2. However, the phylogenetic relationship of these clones remains fully unclear. In this study, we revealed a new genetic clade among O26 strains, ST29C4, which was unique in terms of CRISPR, multilocus sequence type (MLST), and plasmid gene profile (PGP). Moreover, the combination of multiple typing methods demonstrated that this new clone was located in the intermediate phylogenetic position between ST29C3 and other non-ST29C3 strains (i.e., ST29C1-C2 and ST21C1-C2). Overall, the large-scale phylogenetic analysis extends our current understanding of O26 microevolution.

First Comparative Analysis of Clostridium septicum Genomes Provides Insights Into the Taxonomy, Species Genetic Diversity, and Virulence Related to Gas Gangrene

Clostridium septicum is a Gram-positive, toxin-producing, and spore-forming bacterium that is recognized, together with C. perfringens, as the most important etiologic agent of progressive gas gangrene. Clostridium septicum infections are almost always fatal in humans and animals. Despite its clinical and agricultural relevance, there is currently limited knowledge of the diversity and genome structure of C. septicum. This study presents the complete genome sequence of C. septicum DSM 7534T type strain as well as the first comparative analysis of five C. septicum genomes. The taxonomy of C. septicum, as revealed by 16S rRNA analysis as well as by genomic wide indices such as protein-based phylogeny, average nucleotide identity, and digital DNA–DNA hybridization indicates a stable clade. The composition and presence of prophages, CRISPR elements and accessory genetic material was variable in the investigated genomes. This is in contrast to the limited genetic variability described for the phylogenetically and phenotypically related species Clostridium chauvoei. The restriction-modification (RM) systems between two C. septicum genomes were heterogeneous for the RM types they encoded. C. septicum has an open pangenome with 2,311 genes representing the core genes and 1,429 accessory genes. The core genome SNP divergence between genome pairs varied up to 4,886 pairwise SNPs. A vast arsenal of potential virulence genes was detected in the genomes studied. Sequence analysis of these genes revealed that sialidase, hemolysin, and collagenase genes are conserved compared to the α-toxin and hyaluronidase genes. In addition, a conserved gene found in all C. septicum genomes was predicted to encode a leucocidin homolog (beta-channel forming cytolysin) similar (71.10% protein identity) to Clostridium chauvoei toxin A (CctA), which is a potent toxin. In conclusion, our results provide first, valuable insights into strain relatedness and genomic plasticity of C. septicum and contribute to our understanding of the virulence mechanisms of this important human and animal pathogen.

Genome Sequence Analysis of Clostridium chauvoei Strains of European Origin and Evaluation of Typing Options for Outbreak Investigations

Black quarter caused by Clostridium (C.) chauvoei is an important bacterial disease that affects cattle and sheep with high mortality. A comparative genomics analysis of 64 C. chauvoei strains, most of European origin and a few of non-European and unknown origin, was performed. The pangenome analysis showed limited new gene acquisition for the species. The accessory genome involved prophages and genomic islands, with variations in gene composition observed in a few strains. This limited accessory genome may indicate that the species replicates only in the host or that an active CRISPR/Cas system provides immunity to foreign genetic elements. All strains contained a CRISPR type I-B system and it was confirmed that the unique spacer sequences therein can be used to differentiate strains. Homologous recombination events, which may have contributed to the evolution of this pathogen, were less frequent compared to other related species from the genus. Pangenome single nucleotide polymorphism (SNP) based phylogeny and clustering indicate diverse clusters related to geographical origin. Interestingly the identified SNPs were mostly non-synonymous. The study demonstrates the possibility of the existence of polymorphic populations in one host, based on strain variability observed for strains from the same animal and strains from different animals of one outbreak. The study also demonstrates that new outbreak strains are mostly related to earlier outbreak strains from the same farm/region. This indicates the last common ancestor strain from one farm can be crucial to understand the genetic changes and epidemiology occurring at farm level. Known virulence factors for the species were highly conserved among the strains. Genetic elements involved in Nicotinamide adenine dinucleotide (NAD) precursor synthesis (via nadA, nadB, and nadC metabolic pathway) which are known as potential anti-virulence loci are completely absent in C. chauvoei compared to the partial inactivation in C. septicum. A novel core-genome MLST based typing method was compared to sequence typing based on CRISPR spacers to evaluate the usefulness of the methods for outbreak investigations.

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.

Legionella pneumophila CRISPR-Cas Suggests Recurrent Encounters with One or More Phages in the Family Microviridae

Legionella pneumophila is a ubiquitous freshwater pathogen and the causative agent of Legionnaires' disease. L. pneumophila growth within protists provides a refuge from desiccation, disinfection, and other remediation strategies. One outstanding question has been whether this protection extends to phages. L. pneumophila isolates are remarkably devoid of prophages and to date no Legionella phages have been identified. Nevertheless, many L. pneumophila isolates maintain active CRISPR-Cas defenses. So far, the only known target of these systems is an episomal element that we previously named Legionella mobile element 1 (LME-1). The continued expansion of publicly available genomic data promises to further our understanding of the role of these systems. We now describe over 150 CRISPR-Cas systems across 600 isolates to establish the clearest picture yet of L. pneumophila's adaptive defenses. By searching for targets of 1,500 unique CRISPR-Cas spacers, LME-1 remains the only identified CRISPR-Cas targeted integrative element. We identified 3 additional LME-1 variants-all targeted by previously and newly identified CRISPR-Cas spacers-but no other similar elements. Notably, we also identified several spacers with significant sequence similarity to microviruses, specifically those within the subfamily Gokushovirinae. These spacers are found across several different CRISPR-Cas arrays isolated from geographically diverse isolates, indicating recurrent encounters with these phages. Our analysis of the extended Legionella CRISPR-Cas spacer catalog leads to two main conclusions: current data argue against CRISPR-Cas targeted integrative elements beyond LME-1, and the heretofore unknown L. pneumophila phages are most likely lytic gokushoviruses. IMPORTANCE Legionnaires' disease is an often-fatal pneumonia caused by Legionella pneumophila, which normally grows inside amoebae and other freshwater protists. L. pneumophila trades diminished access to nutrients for the protection and isolation provided by the host. One outstanding question is whether L. pneumophila is susceptible to phages, given the protection provided by its intracellular lifestyle. In this work, we use Legionella CRISPR spacer sequences as a record of phage infection to predict that the "missing" L. pneumophila phages belong to the microvirus subfamily Gokushovirinae. Gokushoviruses are known to infect another intracellular pathogen, Chlamydia. How do gokushoviruses access L. pneumophila (and Chlamydia) inside their "cozy niches"? Does exposure to phages happen during a transient extracellular period (during cell-to-cell spread) or is it indicative of a more complicated environmental lifestyle? One thing is clear, 100 years after their discovery, phages continue to hold important secrets about the bacteria upon which they prey.

Characterization and Analysis of Clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs) in Pandemic and Non-Pandemic Vibrio parahaemolyticus Isolates from Seafood Sources

Vibrio parahaemolyticus is one of the significant seafood-borne pathogens causing gastroenteritis in humans. Clustered regularly interspaced short palindromic repeats (CRISPR) are commonly detected in the genomes of V. parahaemolyticus and the polymorphism of CRISPR patterns has been applied as a genetic marker for tracking its evolution. In this work, a total of 15 pandemic and 36 non-pandemic V. parahaemolyticus isolates obtained from seafood between 2000 and 2012 were characterized based on hemolytic activity, antimicrobial susceptibility, and CRISPR elements. The results showed that 15/17 of the V. parahaemolyticus seafood isolates carrying the thermostable direct hemolysin gene (tdh+) were Kanagawa phenomenon (KP) positive. The Multiple Antibiotic Resistance (MAR) index ranged between 0.1 and 0.4, and 45% of the isolates have an MAR index ≥ 0.2. A total of 19 isolates were positive for CRISPR detection, including all tdh+ trh− isolates, two of tdh− trh+, and each of tdh+ trh+ and tdh− trh−. Four spacer types (Sp1 to Sp4) were identified, and CRISPR-positive isolates had at least one type of spacer homolog to the region of Vibrio alginolyticus megaplasmid. It is of interest that a specific CRISPR profile and spacer sequence type was observed with correlations to the hemolysin genotype (tdh/trh). Thus, these provide essential data on the exposure of foreign genetic elements and indicate shared ancestry within different genotypes of V. parahaemolyticus isolates.

CRISPR-Cas systems in Proteus mirabilis

The Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) is a bacterial defense mechanism against bacteriophages composed of two different parts: the CRISPR array and the Cas genes. The spacer acquisition is done by the adaptation module consisting of the hallmark Cas1 Cas2 proteins, which inserts new spacers into the CRISPR array. Here we aimed to describe the CRISPR-Cas system in Proteus mirabilis (P. mirabilis) isolates. CRISPR loci was observed in 30 genomic contents of 109 P. mirabilis isolates that each locus was consisted of two CRISPR arrays and each array had a different preserved leader sequences. Only the type I-E CRISPR-Cas system was common in these isolates. The source of the spacers was identified, including phages and prophages. CRISPR spacer origin analysis also identified a conserved PAM sequence of 5'-AAG-3' nucleotide stretch. Through collecting spacers, CRISPR arrays of P. mirabilis isolates were expanded mostly by integration of bacteriophageal source of spacers. This study shows novel findings in the area of the P-mirabilis CRISPR-Cas system. In this regard, among analyzed genome of P. mirabilis isolates, Class I CRISR-Cas systems were dominant, and all belonged to type I-E. In the flanks of the CRISPR, some other elements with regulatory role were also found. A motif of 11 nt size was found to be preserved among the analyzed genome. We believe that it might has a CRISPR-Cas system transcription facilitator by targeting the Rho element.

Genomic diversity and CRISPR-Cas systems in the cyanobacterium Nostoc in the High Arctic

Nostoc (Nostocales, Cyanobacteria) has a global distribution in the Polar Regions. However, the genomic diversity of Nostoc is little known and there are no genomes available for polar Nostoc. Here we carried out the first genomic analysis of the Nostoc commune morphotype with a recent sample from the High Arctic and a herbarium specimen collected during the British Arctic Expedition (1875-76). Comparisons of the polar genomes with 26 present-day non-polar members of the Nostocales family highlighted that there are pronounced genetic variations among Nostoc strains and species. Osmoprotection and other stress genes were found in all Nostoc strains, but the two Arctic strains had markedly higher numbers of biosynthetic gene clusters for uncharacterised non-ribosomal peptide synthetases, suggesting a high diversity of secondary metabolites. Since viral-host interactions contribute to microbial diversity, we analysed the CRISPR-Cas systems in the Arctic and two temperate Nostoc species. There were a large number of unique repeat-spacer arrays in each genome, indicating diverse histories of viral attack. All Nostoc strains had a subtype I-D system, but the polar specimens also showed evidence of a subtype I-B system that has not been previously reported in cyanobacteria, suggesting diverse cyanobacteria-virus interactions in the Arctic.

Characterization of type I-F CRISPR-Cas system in Laribacter hongkongensis isolates from animals, the environment and diarrhea patients

Laribacter hongkongensis is a foodborne organism that is associated with gastroenteritis and diarrhea in humans. Here we describe the structural characteristics and potential function of CRISPR systems to obtain insight into the genotypic diversity of L. hongkongensis. Specifically, we analyzed the genomic content of six L. hongkongensis genomes and identified two CRISPR loci (CRISPR1 and CRISPR2) belonging to the I-F subtype of CRISPR systems. CRISPR1 was flanked on one side by cas genes and a 170 bp-long putative leader sequence, while CRISPR2 arrays located further and processed by the same cas genes. Then a combination of PCR and sequencing was used to determine the prevalence and distribution of the two CRISPR arrays in 112 L. hongkongensis strains isolated from patients, animals, and water reservoirs. In total, the CRISPR1-Cas system of complete subtype I-F was detected in 91.5% (108/118) of the isolates, whereas CRISPR2 locus existed in 72.0% (85/118). Ten strains only possessed part of the cas genes of subtype I-F and four of them with CRISPR2 array. The two loci contained highly conserved and identical direct repeat sequences which were stable in their RNA secondary structure. Additionally, 2564 total spacers including 980 unique spacers arranged in 59 alleles were identified. Homology analysis showed only 1.8% (18/980) of the spacers matched with plasmid or phage. CRISPR polymorphism present in human isolates and frog isolates was more closely related and more extensive than that of fish isolates based on spacer polymorphism. The elucidation of the structural characteristics of the CRISPR-Cas system may be helpful for further studying the specific mechanism of adaptive immunity and other biological functions mediated by CRISPR in L. hongkongensis. The conservation of CRISPR loci and hypervariable repeat-spacer arrays imply the potential for molecular typing of L. hongkongensis.

Ectopic Spacer Acquisition in Streptococcus thermophilus CRISPR3 Array

Streptococcus thermophilus relies heavily on two type II-A CRISPR-Cas systems, CRISPR1 and CRISPR3, to resist siphophage infections. One hallmark of these systems is the integration of a new spacer at the 5' end of the CRISPR arrays following phage infection. However, we have previously shown that ectopic acquisition of spacers can occur within the CRISPR1 array. Here, we present evidence of the acquisition of new spacers within the array of CRISPR3 of S. thermophilus. The analysis of randomly selected bacteriophage-insensitive mutants of the strain Uy01 obtained after phage infection, as well as the comparison with other S. thermophilus strains with similar CRISPR3 content, showed that a specific spacer within the array could be responsible for misguiding the adaptation complex. These results also indicate that while the vast majority of new spacers are added at the 5' end of the CRISPR array, ectopic spacer acquisition is a common feature of both CRISPR1 and CRISPR3 systems in S. thermophilus, and it can still provide phage resistance. Ectopic spacer acquisition also appears to have occurred naturally in some strains of Streptococcus pyogenes, suggesting that it is a general phenomenon, at least in type II-A systems.

The endless battle between phages and CRISPR-Cas systems in Streptococcus thermophilus

This review describes the contribution of basic research on phage-bacteria interactions to the understanding of CRISPR-Cas systems and their various applications. It focuses on the natural function of CRISPR-Cas systems as adaptive defense mechanisms against mobile genetic elements such as bacteriophage genomes and plasmids. Some of the advances in the characterization of the type II-A CRISPR-Cas system of Streptococcus thermophilus and Streptococcus pyogenes led to the development of the CRISPR-Cas9 genome-editing technology. We mostly discuss the 3 stages of the CRISPR-Cas system in S. thermophilus, namely the adaptation stage, which is unique to this resistance mechanism; the CRISPR RNA biogenesis; and the DNA-cutting activity in the interference stage to protect bacteria against phages. Finally, we look into applications of CRISPR-Cas in microbiology, including overcoming limitations in genome editing.

Comparative Genomic Analysis of Mycobacteriaceae Reveals Horizontal Gene Transfer-Mediated Evolution of the CRISPR-Cas System in the Mycobacterium tuberculosis Complex

Clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated (Cas) genes are conserved genetic elements in many prokaryotes, including Mycobacterium tuberculosis, the causative agent of tuberculosis. Although knowledge of CRISPR locus variability has been utilized in M. tuberculosis strain genotyping, its evolutionary path in Mycobacteriaceae is not well understood. In this study, we have performed a comparative analysis of 141 mycobacterial genomes and identified the exclusive presence of the CRISPR-Cas type III-A system in M. tuberculosis complex (MTBC). Our global phylogenetic analysis of CRISPR repeats and Cas10 proteins offers evidence of horizontal gene transfer (HGT) of the CRISPR-Cas module in the last common ancestor of MTBC and Mycobacterium canettii from a Streptococcus-like environmental bacterium. Additionally, our results show that the variation of CRISPR-Cas organization in M. tuberculosis lineages, especially in the Beijing sublineage of lineage 2, is due to the transposition of insertion sequence IS6110 The direct repeat (DR) region of the CRISPR-Cas locus acts as a hot spot for IS6110 insertion. We show in M. tuberculosis H37Rv that the repeat at the 5' end of CRISPR1 of the forward strand is an atypical repeat made up partly of IS-terminal inverted repeat and partly CRISPR DR. By tracing an undetectable spacer sequence in the DR region, the two CRISPR loci could theoretically be joined to reconstruct the ancestral single CRISPR-Cas locus organization, as seen in M. canettii This study retracing the evolutionary events of HGT and IS6110-driven genomic deletions helps us to better understand the strain-specific variations in M. tuberculosis lineages.IMPORTANCE Comparative genomic analysis of prokaryotes has led to a better understanding of the biology of several pathogenic microorganisms. One such clinically important pathogen is M. tuberculosis, the leading cause of bacterial infection worldwide. Recent evidence on the functionality of the CRISPR-Cas system in M. tuberculosis has brought back focus on these conserved genetic elements, present in many prokaryotes. Our study advances understanding of mycobacterial CRISPR-Cas origin and its diversity among the different species. We provide phylogenetic evidence of acquisition of CRISPR-Cas type III-A in the last common ancestor shared between MTBC and M. canettii, by HGT-mediated events. The most likely source of HGT was an environmental Firmicutes bacterium. Genomic mapping of the CRISPR loci showed the IS6110 transposition-driven variations in M. tuberculosis strains. Thus, this study offers insights into events related to the evolution of CRISPR-Cas in M. tuberculosis lineages.

Analysis of CRISPR-Cas systems in Gardnerella suggests its potential role in the mechanisms of bacterial vaginosis

Bacterial vaginosis (BV) is the principal cause of vaginal discharge among women, and it can lead to many comorbidities with a negative impact in women's daily activities. Despite the fact that the pathophysiological process of BV remains unclear, great advances had been achieved in determining consequences of the shift in the vaginal community, and it was defined that Gardnerella spp., plays a key role in the pathogenesis of BV. Interactions of vaginal phage communities and bacterial hosts may be relevant in eubiosis/dysbiosis states, so defense mechanisms in Gardnerella spp., against phage infections could be relevant in BV development. In this study, we analyzed CRISPR-Cas systems among the 13 Gardnerella species recently classified, considering that these systems act as prokaryotic immune systems against phages, plasmids, and other mobile genetic elements. In silico analyses for CRISPR-Cas systems mining over the 81 Gardnerella spp., strains genomes analyzed led to the identification of subtypes I-E and II-C. Spacers analyses showed a hypervariable region across species, providing a high resolution level in order to distinguish clonality in strains, which was supported with phylogenomic analyses based on Virtual Genomic Fingerprinting. Moreover, most of the spacers revealed interactions between Gardnerella spp., strains and prophages over the genus. Furthermore, virulence traits of the 13 species showed insights of potential niche specificity in the vaginal microbiome. Overall, our results suggest that the CRISPR-Cas systems in the genus Gardnerella may play an important role in the mechanisms of the development and maintenance of BV, considering that the Gardnerella species occupies different niches in the vaginal community; in addition, spacer sequences can be used for genotyping studies.

Utilization of Clustered Regularly Interspaced Short Palindromic Repeats to Genotype Escherichia coli Serogroup O80

The hypervariable nature of clustered regularly interspaced short palindromic repeats (CRISPRs) makes them valuable biomarkers for subtyping and epidemiological investigation of Escherichia coli. Shiga toxin-producing E. coli (STEC) serogroup O80 is one hybrid pathotype that is emerging recently in Europe and is involved in hemolytic uremic syndrome with bacteremia. However, whether STEC O80 strains can be genotyped using CRISPR has not been evaluated. In this study, we aimed to characterize the genetic diversity of 81 E. coli serogroup O80 isolates deposited in the National Center for Biotechnology Information databases using CRISPR typing and to explore the association between virulence potential and CRISPR types (CTs). A total of 21 CTs were identified in 80 O80 strains. CRISRP typing provided discrimination with variants of a single serotype, which suggested a stronger discriminatory power. Based on CRISPR spacer profiles, 70 O80:H2 isolates were further divided into four lineages (lineage LI, LII, LIII, and LIV), which correlated well with whole-genome single nucleotide polymorphisms typing and virulence gene profiles. Moreover, the association between CRISPR lineages and virulence gene profiles hinted that STEC O80:H2 strains may originate from O80:H19 or O80:H26 and that lineage LI may have been evolved from lineage LII. CT2 and CT13 were shared by human and cattle isolates, suggesting that there might be the potential transmission between cattle and human. Collectively, CRISPR typing is one technology that can be used to monitor the transmission of STEC O80 strains and provide new insights into microevolution of serogroup O80.

Novel Genus of Phages Infecting Streptococcus thermophilus: Genomic and Morphological Characterization

Streptococcus thermophilus is a lactic acid bacterium commonly used for the manufacture of yogurt and specialty cheeses. Virulent phages represent a major risk for milk fermentation processes worldwide, as they can inactivate the added starter bacterial cells, leading to low-quality fermented dairy products. To date, four genetically distinct groups of phages infecting S. thermophilus have been described. Here, we describe a fifth group. Phages P738 and D4446 are virulent siphophages that infect a few industrial strains of S. thermophilus The genomes of phages P738 and D4446 were sequenced and found to contain 34,037 and 33,656 bp as well as 48 and 46 open reading frames, respectively. Comparative genomic analyses revealed that the two phages are closely related to each other but display very limited similarities to other S. thermophilus phages. In fact, these two novel S. thermophilus phages share similarities with streptococcal phages of nondairy origin, suggesting that they emerged recently in the dairy environment.IMPORTANCE Despite decades of research and adapted antiphage strategies such as CRISPR-Cas systems, virulent phages are still a persistent risk for the milk fermentation industry worldwide, as they can cause manufacturing failures and alter product quality. Phages P738 and D4446 are novel virulent phages that infect the food-grade Gram-positive bacterial species Streptococcus thermophilus These two related viruses represent a fifth group of S. thermophilus phages, as they are significantly distinct from other known S. thermophilus phages. Both phages share similarities with phages infecting nondairy streptococci, suggesting their recent emergence and probable coexistence in dairy environments. These findings highlight the necessity of phage surveillance programs as the phage population evolves in response to the application of antiphage strategies.

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/.

Milk microbial composition of Brazilian dairy cows entering the dry period and genomic comparison between Staphylococcus aureus strains susceptible to the bacteriophage vB_SauM-UFV_DC4

Brazil has the second-largest dairy cattle herd in the world, and bovine mastitis still can cause significant losses for dairy farmers. Despite this fact, little information is available about milk microbial composition of Brazilian dairy cows, as well as the potential use of bacteriophages in the control of S. aureus. Here, we investigated milk bacterial composition of 28 Holstein Fresian cows (109 teats), selected in the dry-off period, using 16S rRNA analysis. Furthermore, a representative S. aureus strain (UFV2030RH1) was obtained at drying-off for isolation of a bacteriophage (vB_SauM-UFV_DC4, UFV_DC4) and bacterial genomic comparison purposes. Our outcomes revealed that Staphylococcus was the third most prevalent genus and positively correlated with subclinical mastitis events. As a major finding, genomic analyses showed the presence of adhesive matrix molecules that recognize microbial surface components (MSCRAMM) in UFV2030RH1 and might indicate great biofilm formation capability. A minimum inhibitory concentration (MIC) assay showed that resistance to ampicillin was the highest among the antibiotic tested in S. aureus 3059 and UFV2030RH1, displaying values four and sixteen times greater than MIC resistance breakpoint, respectively. Together, our results suggest that Staphylococcus is highly prevalent in dairy cows at drying-off and the use of the phage UFV_DC4 as a biocontrol agent must be investigated in future studies.

Type I-F CRISPR-Cas Distribution and Array Dynamics in Legionella pneumophila

In bacteria and archaea, several distinct types of CRISPR-Cas systems provide adaptive immunity through broadly similar mechanisms: short nucleic acid sequences derived from foreign DNA, known as spacers, engage in complementary base pairing with invasive genetic elements setting the stage for nucleases to degrade the target DNA. A hallmark of type I CRISPR-Cas systems is their ability to acquire spacers in response to both new and previously encountered invaders (naïve and primed acquisition, respectively). Our phylogenetic analyses of 43 L. pneumophila type I-F CRISPR-Cas systems and their resident genomes suggest that many of these systems have been horizontally acquired. These systems are frequently encoded on plasmids and can co-occur with nearly identical chromosomal loci. We show that two such co-occurring systems are highly protective and undergo efficient primed acquisition in the lab. Furthermore, we observe that targeting by one system’s array can prime spacer acquisition in the other. Lastly, we provide experimental and genomic evidence for a model in which primed acquisition can efficiently replenish a depleted type I CRISPR array following a mass spacer deletion event.

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/.

Monitoring of Lactobacillus sanfranciscensis strains during wheat and rye sourdough fermentations by CRISPR locus length polymorphism PCR

Lactobacillus (L.) sanfranciscensis is a competitive key species in sourdough fermentations. However, the principles involved in establishing the commonly observed phenomenon of strain dominance are unresolved. This has been studied little because the methods for fast and reliable differentiation of strains and their monitoring during fermentation are tedious and cannot be done with large numbers of isolates. In this contribution, we present a strain-specific, PCR-based typing method that uses length heterogeneities of the clustered regularly interspaced short palindromic repeats (CRISPR) loci as they occur in the genomes of different strains. In silico analysis of 21 genomes revealed 14 different CRISPR genotypes. We then designed a primer set to simultaneously detect different strains in a multiplex PCR assay designated CRISPR locus length polymorphism PCR (CLLP-PCR). The usefulness of this method was evaluated in lab-scale sourdough fermentations conducted with rye and wheat flours. First, the flour was mixed with water to a dough yield of 200. Then each dough was inoculated with four different L. sanfranciscensis strains (TMW 1.1150, TMW 1.392, TMW 1.2142, and TMW 1.2138) at levels of 10⁹ cfu/g each. Sourdoughs were propagated at 28 °C for 5 days by back slopping 5% to the flour mass every 24 h. Samples were collected each day; DNA was isolated, and the presence of strains was detected qualitatively in the sourdoughs with PCR. L. sanfranciscensis TMW 1.392 became dominant as early as 2 days into the fermentation and remained the only detectable strain for the rest of the sampling period. CLLP-PCR proved to be useful in investigating the assertiveness of different strains of L. sanfranciscensis in sourdoughs. Therefore, CLLP-PCR may be used as a tool to investigate assertiveness of microorganisms in food fermentations at the strain level.

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.

Polymorphism of Type I-F CRISPR/Cas system in Escherichia coli of phylogenetic group B2 and its application in genotyping

E. coli of phylogenetic group B2 is responsible for many extraintestinal infections, posing a great threat to health. The relatively polymorphic nature of CRISPR in phylogenetically related E. coli strains makes them potential markers for bacterial typing and evolutionary studies. In the current work, we investigated the occurrence and diversity of CRISPR/Cas system and explored its potential for genotyping. Type I-F CRISPR/Cas systems were found in 413 of 1190 strains of E. coli and exhibited the clustering within certain CCs and STs. And CRISPR spacer contents correlated well with MLST types. The divergence analysis of CRISPR showed stronger discriminatory power than MLST, and CRISPR polymorphism was instrumental for differentiating highly closely related strains. The timeline of spacer acquisition and deletion provided important information for inferring the evolution model between distinct serotypes. Identical spacer sequences were shared by strains with the same H-antigen type but not strains with the same O-antigen type. The homology between spacers and antibiotic-resistant plasmids demonstrated the role of Type I-F system in limiting the acquisition of antimicrobial resistance. Collectively, our data presents the dynamic nature of Type I-F CRISPR in E. coli of phylogenetic group B2 and provides new insights into the application of CRISPR-based typing in the species.