CRISPRTarget: bioinformatic prediction and analysis of crRNA targets

Characterization of CRISPR array in Salmonella enterica from asymptomatic people and patients

Salmonella enterica is a major foodborne pathogen causing symptomatic diseases or asymptomatic infections in humans. To reveal the genetic difference of Salmonella strains from patients to that from asymptomatic people, we used CRISPR typing to analyze the phylogenetic relationship of 180 clinical strains during 2017-2018 in Jiangsu, China. The CRISPR typing divided these isolates into 76 CRISPR types with a discriminatory power of 97.6%. S. Typhimurium and its monophasic variants of 6 CRISPR types are the significant serotypes causing both human diseases and asymptomatic infection, while S. Enteritidis mainly resulted in diseases and shared one CRISPR type. The spacer HadB20 displayed as a new molecular marker to differentiate ST34-S. Typhimurium monophasic variant from ST19-S. Typhimurium. S. Derby, S. London, and S. Senftenberg frequently caused asymptomatic infection with diverse CRISPR types, while S. Mbandaka and S. Meleagridis, occasionally isolated from patients, had conserved CRISPR types. Additionally, 30 of 516 newly identified spacers showed homology to sequences in both plasmids and bacteriophages. Interestingly, some spacers from one serotype showed homology to the correspondent prophage or plasmid sequences in another serotype; and more than two spacers identified in one strain showed homology to the sequences located in the identical plasmids or phages, revealing the constant evolution of Salmonella CRISPR arrays during the interactions between bacteria and phages or plasmids.

Lifestyle of sponge symbiont phages by host prediction and correlative microscopy

Bacteriophages (phages) are ubiquitous elements in nature, but their ecology and role in animals remains little understood. Sponges represent the oldest known extant animal-microbe symbiosis and are associated with dense and diverse microbial consortia. Here we investigate the tripartite interaction between phages, bacterial symbionts, and the sponge host. We combined imaging and bioinformatics to tackle important questions on who the phage hosts are and what the replication mode and spatial distribution within the animal is. This approach led to the discovery of distinct phage-microbe infection networks in sponge versus seawater microbiomes. A new correlative in situ imaging approach ('PhageFISH-CLEM') localised phages within bacterial symbiont cells, but also within phagocytotically active sponge cells. We postulate that the phagocytosis of free virions by sponge cells modulates phage-bacteria ratios and ultimately controls infection dynamics. Prediction of phage replication strategies indicated a distinct pattern, where lysogeny dominates the sponge microbiome, likely fostered by sponge host-mediated virion clearance, while lysis dominates in seawater. Collectively, this work provides new insights into phage ecology within sponges, highlighting the importance of tripartite animal-phage-bacterium interplay in holobiont functioning. We anticipate that our imaging approach will be instrumental to further understanding of viral distribution and cellular association in animal hosts.

Ecogenomics of Groundwater Phages Suggests Niche Differentiation Linked to Specific Environmental Tolerance

Viruses are ubiquitous microbiome components, shaping ecosystems via strain-specific predation, horizontal gene transfer and redistribution of nutrients through host lysis. Viral impacts are important in groundwater ecosystems, where microbes drive many nutrient fluxes and metabolic processes; however, little is known about the diversity of viruses in these environments. We analyzed four groundwater plasmidomes (the entire plasmid content of an environment) and identified 200 viral sequences, which clustered into 41 genus-level viral clusters (approximately equivalent to viral genera) including 9 known and 32 putative new genera. We used publicly available bacterial whole-genome sequences (WGS) and WGS from 261 bacterial isolates from this groundwater environment to identify potential viral hosts. We linked 76 of the 200 viral sequences to a range of bacterial phyla, the majority associated with Proteobacteria, followed by Firmicutes, Bacteroidetes, and Actinobacteria. The publicly available WGS enabled mapping bacterial hosts to several viral sequences. The WGS of groundwater isolates increased the depth of host prediction by allowing host identification at the strain level. The latter included 4 viruses that were almost entirely (>99% query coverage, >99% identity) identified as integrated in the genomes of Pseudomonas, Acidovorax, and Castellaniella strains, resulting in high-confidence host assignments. Lastly, 21 of these viruses carried putative auxiliary metabolite genes for metal and antibiotic resistance, which might drive their infection cycles and/or provide selective advantage to infected hosts. Exploring the groundwater virome provides a necessary foundation for integration of viruses into ecosystem models where they are key players in microbial adaption to environmental stress. IMPORTANCE To our knowledge, this is the first study to identify the bacteriophage distribution in a groundwater ecosystem shedding light on their prevalence and distribution across metal-contaminated and background sites. Our study is uniquely based on selective sequencing of solely the extrachromosomal elements of a microbiome followed by analysis for viral signatures, thus establishing a more focused approach for phage identifications. Using this method, we detected several novel phage genera along with those previously established. Our approach of using the whole-genome sequences of hundreds of bacterial isolates from the same site enabled us to make host assignments with high confidence, several at strain levels. Certain phage genes suggest that they provide an environment-specific selective advantage to their bacterial hosts. Our study lays the foundation for future research on directed phage isolations using specific bacterial host strains to further characterize groundwater phages, their life cycles, and their effects on groundwater microbiome and biogeochemistry.

Base editing-coupled survival screening enabled high-sensitive analysis of PAM compatibility and finding of the new possible off-target

Base editing (BE) is a promising genome engineering tool for modifying DNA or RNA and has been widely used in various microorganisms as well as eukaryotic cells. Despite the proximal protospacer adjacent motif (PAM) is critical to the targeting range and off-target effect of BE, there is still lack of a specific approach to analyze the PAM pattern in BE systems. Here, we developed a base editing-coupled survival screening method. Using dCas9 from Streptococcus pyogenes (SpdCas9) and its variants xdCas9 3.7 and dCas9 NG as example, their PAM patterns in BE systems were extensively characterized using the NNNN PAM library with high sensitivity. In addition to the typical PAM recognition features, we observed more unique PAMs exhibiting BE activity. These PAM patterns will boost the finding of potential off-target editing event arising from non-canonical PAMs and provide the guidelines for PAM usage in the BE system.

Pan-Genome Analysis Reveals Host-Specific Functional Divergences in Burkholderia gladioli

Burkholderia gladioli has high versatility and adaptability to various ecological niches. Here, we constructed a pan-genome using 14 genome sequences of B. gladioli, which originate from different niches, including gladiolus, rice, humans, and nature. Functional roles of core and niche-associated genomes were investigated by pathway enrichment analyses. Consequently, we inferred the uniquely important role of niche-associated genomes in (1) selenium availability during competition with gladiolus host; (2) aromatic compound degradation in seed-borne and crude oil-accumulated environments, and (3) stress-induced DNA repair system/recombination in the cystic fibrosis-niche. We also identified the conservation of the rhizomide biosynthetic gene cluster in all the B. gladioli strains and the concentrated distribution of this cluster in human isolates. It was confirmed the absence of complete CRISPR/Cas system in both plant and human pathogenic B. gladioli and the presence of the system in B. gladioli living in nature, possibly reflecting the inverse relationship between CRISPR/Cas system and virulence.

Programmable RNA editing with compact CRISPR-Cas13 systems from uncultivated microbes

Competitive coevolution between microbes and viruses has led to the diversification of CRISPR-Cas defense systems against infectious agents. By analyzing metagenomic terabase datasets, we identified two compact families (775 to 803 amino acids (aa)) of CRISPR-Cas ribonucleases from hypersaline samples, named Cas13X and Cas13Y. We engineered Cas13X.1 (775 aa) for RNA interference experiments in mammalian cell lines. We found Cas13X.1 could tolerate single-nucleotide mismatches in RNA recognition, facilitating prophylactic RNA virus inhibition. Moreover, a minimal RNA base editor, composed of engineered deaminase (385 aa) and truncated Cas13X.1 (445 aa), exhibited robust editing efficiency and high specificity to induce RNA base conversions. Our results suggest that there exist untapped bacterial defense systems in natural microbes that can function efficiently in mammalian cells, and thus potentially are useful for RNA-editing-based research.

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.

SpacePHARER: Sensitive identification of phages from CRISPR spacers in prokaryotic hosts

Summary

SpacePHARER (CRISPR Spacer Phage–Host Pair Finder) is a sensitive and fast tool for de novo prediction of phage–host relationships via identifying phage genomes that match CRISPR spacers in genomic or metagenomic data. SpacePHARER gains sensitivity by comparing spacers and phages at the protein level, optimizing its scores for matching very short sequences, and combining evidence from multiple matches, while controlling for false positives. We demonstrate SpacePHARER by searching a comprehensive spacer list against all complete phage genomes.

Availability and implementation

SpacePHARER is available as an open-source (GPLv3), user-friendly command-line software for Linux and macOS: https://github.com/soedinglab/spacepharer.

Supplementary information

Supplementary data are available at Bioinformatics online.

Genetic characteristics and phylogenetic analysis of Brazilian clinical strains of Pseudomonas aeruginosa harboring CRISPR/Cas systems

The CRISPR-Cas are adaptive immune systems found in archaea and bacteria, responsible for providing sequence-specific resistance against foreign DNA. Strains of Pseudomonas aeruginosa may carry CRISPR/Cas system types I-F, I-E and/or I-C; however, several aspects related to the epidemiology and functionality of these systems have not yet been revealed. Here, we report 13 genomes of clinical strains of P. aeruginosa from Brazil that were positive for CRISPR/Cas system types I-F and I-E, a rare feature in this species. The phylogenetic tree, which was constructed with 161 other publicly available genomes, suggested no direct relationship between positive strains, and the various types of CRISPR/Cas systems were spread throughout the tree. Comparative analysis of the genetic locations of type I-F and a specific orphan CRISPR array (without cas genes), named the LES locus, showed sequence similarities between this orphan locus and type I-F, but these LES loci were inserted in a different genomic location. We also report the presence of prophages, the presence of anti-CRISPR genes, and possibly the presence of self-targeting spacers. Here, we conclude that CRISPR/Cas is highly associated with certain lineages and is spread throughout the phylogenetic tree, showing no clear pattern of evolutionary distribution. Moreover, the LES locus might be a CRISPR1 locus related to type I-F that may have been misplaced and maintained over time. Furthermore, strains carrying I-F and I-E are rare, and not all of them are closely related. Further functional work is needed to better comprehend if aspects reported in this study are functional, including the LES locus, self-targeting spacers, anti-CRISPR protection, and I-F/I-E-carrying lineages.

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.

Scalable characterization of the PAM requirements of CRISPR-Cas enzymes using HT-PAMDA

The continued expansion of the genome-editing toolbox necessitates methods to characterize important properties of CRISPR-Cas enzymes. One such property is the requirement for Cas proteins to recognize a protospacer-adjacent motif (PAM) in DNA target sites. The high-throughput PAM determination assay (HT-PAMDA) is a method that enables scalable characterization of the PAM preferences of different Cas proteins. Here, we provide a step-by-step protocol for the method, discuss experimental design considerations, and highlight how the method can be used to profile naturally occurring CRISPR-Cas9 enzymes, engineered derivatives with improved properties, orthologs of different classes (e.g., Cas12a), and even different platforms (e.g., base editors). A distinguishing feature of HT-PAMDA is that the enzymes are expressed in a cell type or organism of interest (e.g., mammalian cells), permitting scalable characterization and comparison of hundreds of enzymes in a relevant setting. HT-PAMDA does not require specialized equipment or expertise and is cost effective for multiplexed characterization of many enzymes. The protocol enables comprehensive PAM characterization of dozens or hundreds of Cas enzymes in parallel in <2 weeks.

Mobile element warfare via CRISPR and anti-CRISPR in Pseudomonas aeruginosa

Bacteria deploy multiple defenses to prevent mobile genetic element (MGEs) invasion. CRISPR-Cas immune systems use RNA-guided nucleases to target MGEs, which counter with anti-CRISPR (Acr) proteins. Our understanding of the biology and co-evolutionary dynamics of the common Type I-C CRISPR-Cas subtype has lagged because it lacks an in vivo phage-host model system. Here, we show the anti-phage function of a Pseudomonas aeruginosa Type I-C CRISPR-Cas system encoded on a conjugative pKLC102 island, and its Acr-mediated inhibition by distinct MGEs. Seven genes with anti-Type I-C function (acrIC genes) were identified, many with highly acidic amino acid content, including previously described DNA mimic AcrIF2. Four of the acr genes were broad spectrum, also inhibiting I-E or I-F P. aeruginosa CRISPR-Cas subtypes. Dual inhibition comes at a cost, however, as simultaneous expression of Type I-C and I-F systems renders phages expressing the dual inhibitor AcrIF2 more sensitive to targeting. Mutagenesis of numerous acidic residues in AcrIF2 did not impair anti-I-C or anti-I-F function per se but did exacerbate inhibition defects during competition, suggesting that excess negative charge may buffer DNA mimics against competition. Like AcrIF2, five of the Acr proteins block Cascade from binding DNA, while two function downstream, likely preventing Cas3 recruitment or activity. One such inhibitor, AcrIC3, is found in an 'anti-Cas3' cluster within conjugative elements, encoded alongside bona fide Cas3 inhibitors AcrIF3 and AcrIE1. Our findings demonstrate an active battle between an MGE-encoded CRISPR-Cas system and its diverse MGE targets.

CRISPR-SE: a brute force search engine for CRISPR design

CRISPR is a revolutionary genome-editing tool that has been broadly used and integrated within novel biotechnologies. A major component of existing CRISPR design tools is the search engines that find the off-targets up to a predefined number of mismatches. Many CRISPR design tools adapted sequence alignment tools as the search engines to speed up the process. These commonly used alignment tools include BLAST, BLAT, Bowtie, Bowtie2 and BWA. Alignment tools use heuristic algorithm to align large amount of sequences with high performance. However, due to the seed-and-extend algorithms implemented in the sequence alignment tools, these methods are likely to provide incomplete off-targets information for ultra-short sequences, such as 20-bp guide RNAs (gRNA). An incomplete list of off-targets sites may lead to erroneous CRISPR design. To address this problem, we derived four sets of gRNAs to evaluate the accuracy of existing search engines; further, we introduce a search engine, namely CRISPR-SE. CRISPR-SE is an accurate and fast search engine using a brute force approach. In CRISPR-SE, all gRNAs are virtually compared with query gRNA, therefore, the accuracies are guaranteed. We performed the accuracy benchmark with multiple search engines. The results show that as expected, alignment tools reported an incomplete and varied list of off-target sites. CRISPR-SE performs well in both accuracy and speed. CRISPR-SE will improve the quality of CRISPR design as an accurate high-performance search engine.

Widespread distribution of prophages signaling the potential for adaptability and pathogenicity evolution of Ralstonia solanacearum species complex

Integrated bacteriophages (prophages) can impact host cells, affecting their lifestyle, genomic diversity, and fitness. However, many basic aspects of how these organisms affect the host cell remain poorly understood. Ralstonia solanacearum is a gram-negative plant pathogenic bacterium that encompasses a great diversity of ecotypes regarded as a species complex (R. solanacearum Species Complex - RSSC). RSSC genomes have a mosaic structure containing numerous elements, signaling the potential for its evolution through horizontal gene transfer. Here, we analyzed 120 Ralstonia spp. genomes from the public database to identify prophage sequences. In total, 379 prophage-like elements were found in the chromosome and megaplasmid of Ralstonia spp. These elements encode genes related to host fitness, virulence factors, antibiotic resistance, and niche adaptation, which might contribute to RSSC adaptability. Prophage-like elements are widespread into the complex in different species and geographic origins, suggesting that the RSSC phages are ancestrally acquired. Complete prophages belonging to the families Inoviridae, Myoviridae, and Siphoviridae were found, being the members of Inoviridae the most abundant. Analysis of CRISPR-Cas spacer sequences demonstrated the presence of prophages sequences that indicate successive infection events during bacterial evolution. Besides complete prophages, we also demonstrated 14 novel putative prophages integrated into Ralstonia spp. genomes. Altogether, our results provide insights into the diversity of prophages in RSSC genomes and suggest that these elements may deeply affect the virulence and host adaptation and shaping the genomes among the strains of this important pathogen.

Repurposing the Endogenous Type I-E CRISPR/Cas System for Gene Repression in Gluconobacter oxydans WSH-003

Gluconobacter oxydans is well-known for its incomplete oxidizing capacity and has been widely applied in industrial production. However, genetic tools in G. oxydans are still scarce compared with model microorganisms, limiting its metabolic engineering. This study aimed to develop a clustered regularly interspaced short palindromic repeats interference (CRISPRi) system based on the typical type I-E endogenous CRISPR/CRISPR-associated proteins (Cas) system in G. oxydans WSH-003. The nuclease Cas3 in this system was inactivated naturally and hence did not need to be knocked out. Subsequently, the CRISPRi effect was verified by repressing the expression of fluorescent proteins, revealing effective multiplex gene repression. Finally, the endogenous CRISPRi system was used to study the role of the central carbon metabolism pathway, including the pentose phosphate pathway (PPP) and Entner-Doudoroff pathway (EDP), in G. oxydans WSH-003. This was done to demonstrate a metabolic engineering application. The PPP was found to be important for cell growth and the substrate conversion rate. The development of the CRISPRi system enriched the gene regulation tools in G. oxydans and promoted the metabolic engineering modification of G. oxydans to improve its performance. In addition, it might have implications for metabolic engineering modification of other genetically recalcitrant strains.

Roles of bacteriophages, plasmids and CRISPR immunity in microbial community dynamics revealed using time-series integrated meta-omics

Viruses and plasmids (invasive mobile genetic elements (iMGEs)) have important roles in shaping microbial communities, but their dynamic interactions with CRISPR-based immunity remain unresolved. We analysed generation-resolved iMGE-host dynamics spanning one and a half years in a microbial consortium from a biological wastewater treatment plant using integrated meta-omics. We identified 31 bacterial metagenome-assembled genomes encoding complete CRISPR-Cas systems and their corresponding iMGEs. CRISPR-targeted plasmids outnumbered their bacteriophage counterparts by at least fivefold, highlighting the importance of CRISPR-mediated defence against plasmids. Linear modelling of our time-series data revealed that the variation in plasmid abundance over time explained more of the observed community dynamics than phages. Community-scale CRISPR-based plasmid-host and phage-host interaction networks revealed an increase in CRISPR-mediated interactions coinciding with a decrease in the dominant 'Candidatus Microthrix parvicella' population. Protospacers were enriched in sequences targeting genes involved in the transmission of iMGEs. Understanding the factors shaping the fitness of specific populations is necessary to devise control strategies for undesirable species and to predict or explain community-wide phenotypes.

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

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

Genomic Characterization of a Novel Tenericutes Bacterium from Deep-Sea Holothurian Intestine

Intestinal bacterial communities are highly relevant to the digestion, nutrition, growth, reproduction, and immunity of animals, but little is known about the composition and function of intestinal microbiota in deep-sea invertebrates. In this study, the intestinal microbiota of six holothurian Molpadia musculus were investigated, showing that their midguts were predominantly occupied by Izemoplasmatales bacteria. Using metagenomic sequencing, a draft genome of 1,822,181 bp was successfully recovered. After comparison with phylogenetically related bacteria, genes involved in saccharide usage and de novo nucleotide biosynthesis were reduced. However, a set of genes responsible for extracellular nucleoside utilization and 14 of 20 amino acid synthesis pathways were completely retained. Under oligotrophic condition, the gut-associated bacterium may make use of extracellular DNA for carbon and energy supplement, and may provide essential amino acids to the host. The clustered regularly interspaced short palindromic repeat (CRISPR) and restriction–modification (RM) systems presented in the genome may provide protection against invading viruses. A linear azol(in)e-containing peptide gene cluster for bacteriocin synthesize was also identified, which may inhibit the colonization and growth of harmful bacteria. Known virulence factors were not found by database searching. On the basis of its phylogenetic position and metabolic characteristics, we proposed that the bacterium represented a novel genus and a novel family within the Izemoplasmatales order and suggested it be named “Candidatus Bathyoplasma sp. NZ”. This was the first time describing host-associated Izemoplasmatales.

Multidrug-Resistant Proteus mirabilis Strain with Cointegrate Plasmid

Proteus mirabilis is a component of the normal intestinal microflora of humans and animals, but can cause urinary tract infections and even sepsis in hospital settings. In recent years, the number of multidrug-resistant P. mirabilis isolates, including the ones producing extended-spectrum β-lactamases (ESBLs), is increasing worldwide. However, the number of investigations dedicated to this species, especially, whole-genome sequencing, is much lower in comparison to the members of the ESKAPE pathogens group. This study presents a detailed analysis of clinical multidrug-resistant ESBL-producing P. mirabilis isolate using short- and long-read whole-genome sequencing, which allowed us to reveal possible horizontal gene transfer between Klebsiella pneumoniae and P. mirabilis plasmids and to locate the CRISPR-Cas system in the genome together with its probable phage targets, as well as multiple virulence genes. We believe that the data presented will contribute to the understanding of antibiotic resistance acquisition and virulence mechanisms for this important pathogen.

Discovery of multiple anti-CRISPRs highlights anti-defense gene clustering in mobile genetic elements

Many prokaryotes employ CRISPR-Cas systems to combat invading mobile genetic elements (MGEs). In response, some MGEs have developed strategies to bypass immunity, including anti-CRISPR (Acr) proteins; yet the diversity, distribution and spectrum of activity of this immune evasion strategy remain largely unknown. Here, we report the discovery of new Acrs by assaying candidate genes adjacent to a conserved Acr-associated (Aca) gene, aca5, against a panel of six type I systems: I-F (Pseudomonas, Pectobacterium, and Serratia), I-E (Pseudomonas and Serratia), and I-C (Pseudomonas). We uncover 11 type I-F and/or I-E anti-CRISPR genes encoded on chromosomal and extrachromosomal MGEs within Enterobacteriaceae and Pseudomonas, and an additional Aca (aca9). The acr genes not only associate with other acr genes, but also with genes encoding inhibitors of distinct bacterial defense systems. Thus, our findings highlight the potential exploitation of acr loci neighborhoods for the identification of previously undescribed anti-defense systems.

Phage gene expression and host responses lead to infection-dependent costs of CRISPR immunity

CRISPR-Cas immune systems are widespread in bacteria and archaea, but not ubiquitous. Previous work has demonstrated that CRISPR immunity is associated with an infection-induced fitness cost, which may help explain the patchy distribution observed. However, the mechanistic basis of this cost has remained unclear. Using Pseudomonas aeruginosa PA14 and its phage DMS3vir as a model, we perform a 30-day evolution experiment under phage mediated selection. We demonstrate that although CRISPR is initially selected for, bacteria carrying mutations in the phage receptor rapidly invade the population following subsequent reinfections. We then test three potential mechanisms for the observed cost of CRISPR: (1) autoimmunity from the acquisition of self-targeting spacers, (2) immunopathology or energetic costs from increased cas gene expression and (3) toxicity caused by phage gene expression prior to CRISPR-mediated cleavage. We find that phages can express genes before the immune system clears the infection and that expression of these genes can have a negative effect on host fitness. While infection does not lead to increased expression of cas genes, it does cause differential expression of multiple other host processes that may further contribute to the cost of CRISPR immunity. In contrast, we found little support for infection-induced autoimmunological and immunopathological effects. Phage gene expression prior to cleavage of the genome by the CRISPR-Cas immune system is therefore the most parsimonious explanation for the observed phage-induced fitness cost.

Rapid and assured genetic engineering methods applied to Acinetobacter baylyi ADP1 genome streamlining

One goal of synthetic biology is to improve the efficiency and predictability of living cells by removing extraneous genes from their genomes. We demonstrate improved methods for engineering the genome of the metabolically versatile and naturally transformable bacterium Acinetobacter baylyi ADP1 and apply them to a genome streamlining project. In Golden Transformation, linear DNA fragments constructed by Golden Gate Assembly are directly added to cells to create targeted deletions, edits, or additions to the chromosome. We tested the dispensability of 55 regions of the ADP1 chromosome using Golden Transformation. The 18 successful multiple-gene deletions ranged in size from 21 to 183 kb and collectively accounted for 23.4% of its genome. The success of each multiple-gene deletion attempt could only be partially predicted on the basis of an existing collection of viable ADP1 single-gene deletion strains and a new transposon insertion sequencing (Tn-Seq) dataset that we generated. We further show that ADP1’s native CRISPR/Cas locus is active and can be retargeted using Golden Transformation. We reprogrammed it to create a CRISPR-Lock, which validates that a gene has been successfully removed from the chromosome and prevents it from being reacquired. These methods can be used together to implement combinatorial routes to further genome streamlining and for more rapid and assured metabolic engineering of this versatile chassis organism.

Source Tracking Based on Core Genome SNV and CRISPR Typing of Salmonella enterica Serovar Heidelberg Isolates Involved in Foodborne Outbreaks in Québec, 2012

Whole-genome sequencing (WGS) is the method of choice for bacterial subtyping and it is rapidly replacing the more traditional methods such as pulsed-field gel electrophoresis (PFGE). Here we used the high-resolution core genome single nucleotide variant (cgSNV) typing method to characterize clinical and food from Salmonella enterica serovar Heidelberg isolates in the context of source attribution. Additionally, clustered regularly interspaced short palindromic repeats (CRISPR) analysis was included to further support this method. Our results revealed that cgSNV was highly discriminatory and separated the outbreak isolates into distinct clusters (0-4 SNVs). CRISPR analysis was also able to distinguish outbreak strains from epidemiologically unrelated isolates. Specifically, our data clearly demonstrated the strength of these two methods to determine the probable source(s) of a 2012 epidemiologically characterized outbreak of S. Heidelberg. Using molecular cut-off of 0-10 SNVs, the cgSNV analysis of 246 clinical and food isolates of S. Heidelberg collected in Québec, in the same year of the outbreak event, revealed that retail and abattoir chicken isolates likely represent an important source of human infection to S. Heidelberg. Interestingly, the isolates genetically related by cgSNV also harbored the same CRISPR as outbreak isolates and clusters. This indicates that CRISPR profiles can be useful as a complementary approach to determine source attribution in foodborne outbreaks. Use of the genomic analysis also allowed to identify a large number of cases that were missed by PFGE, indicating that most outbreaks are probably underestimated. Although epidemiological information must still support WGS-based results, cgSNV method is a highly discriminatory method for the resolution of outbreak events and the attribution of these events to their respective sources. CRISPR typing can serve as a complimentary tool to this analysis during source tracking.

Exploring potential of vaginal Lactobacillus isolates from South African women for enhancing treatment for bacterial vaginosis

Antibiotics continue to be the standard-of-care for bacterial vaginosis (BV), although recurrence rates are high. Vaginal probiotics may improve durability of BV treatment, although few probiotics for vaginal health contain Lactobacillus spp. that commonly colonize the lower female genital tract. Characteristics of vaginal Lactobacillus strains from South African women were evaluated for their probiotic potential in vitro compared to strains from commercial vaginal products, including growth at varying pHs, ability to lower pH, produce D-/L-lactate and H2O2, influence growth of BV-associated Gardnerella vaginalis and Prevotella bivia, adherence to cervical cells and susceptibility to antibiotics. Fifty-seven Lactobacillus strains were purified from cervico-vaginal fluid, including L. crispatus, L. jensenii, L. gasseri, L. mucosae, and L. vaginalis. L crispatus strains grew better at pHs below 4.5 and lowered pH more effectively than other strains. Production of D-/L-lactate and H2O2 varied between Lactobacillus species and strains. Lactobacillus strains generally inhibited P. bivia more uniformly than G. vaginalis isolates. All vaginal Lactobacillus isolates were resistant to metronidazole while susceptibility to clindamycin varied. Furthermore, vaginal Lactobacillus strains tended to be broadly susceptible to penicillin, amoxicillin, rifampicin and rifabutin. Whole-genome-sequencing of five of the best-performing vaginal Lactobacillus strains confirmed their likely safety, due to antimicrobial resistance elements being largely absent, while putative intact prophages were present in the genomes of two of the five strains. Overall, vaginal Lactobacillus strains largely performed better in these in vitro assays than probiotic strains currently used in probiotics for vaginal health. Including the best-performing vaginal Lactobacillus isolates in a region-specific probiotic for vaginal health may result in improved BV treatment options.

Repurposing the Native Type I-F CRISPR-Cas System in Pseudomonas aeruginosa for Genome Editing

Repurposing the broadly distributed native CRISPR-Cas systems in prokaryotes for genome editing is emerging as a new strategy for genetic manipulations. We recently reported the establishment of a single plasmid-mediated, one-step genome-editing technique in a multidrug-resistant genotype of the opportunistic pathogen Pseudomonas aeruginosa by harnessing its endogenous type I-F CRISPR-Cas system. The platform is readily applicable in additional type I-F CRISPR-containing clinical and environmental P. aeruginosa isolates. Herein, we provide the detailed protocol for the methodology.

For complete details on the establishment and exploitation of this protocol, please refer to Xu et al. (2019).

•

Workflow to repurpose the native I-F CRISPR-Cas system for genome editing in PA154197

•

Procedures to characterize native CRISPR-Cas systems in the strains of interest

•

Potential solutions for troubleshooting and expanding the method to other strains

Characterization and repurposing of the endogenous Type I-F CRISPR-Cas system of Zymomonas mobilis for genome engineering

Application of CRISPR-based technologies in non-model microorganisms is currently very limited. Here, we reported efficient genome engineering of an important industrial microorganism, Zymomonas mobilis, by repurposing the endogenous Type I-F CRISPR–Cas system upon its functional characterization. This toolkit included a series of genome engineering plasmids, each carrying an artificial self-targeting CRISPR and a donor DNA for the recovery of recombinants. Through this toolkit, various genome engineering purposes were efficiently achieved, including knockout of ZMO0038 (100% efficiency), cas2/3 (100%), and a genomic fragment of >10 kb (50%), replacement of cas2/3 with mCherry gene (100%), in situ nucleotide substitution (100%) and His-tagging of ZMO0038 (100%), and multiplex gene deletion (18.75%) upon optimal donor size determination. Additionally, the Type I-F system was further applied for CRISPRi upon Cas2/3 depletion, which has been demonstrated to successfully silence the chromosomally integrated mCherry gene with its fluorescence intensity reduced by up to 88%. Moreover, we demonstrated that genome engineering efficiency could be improved under a restriction–modification (R–M) deficient background, suggesting the perturbance of genome editing by other co-existing DNA targeting modules such as the R–M system. This study might shed light on exploiting and improving CRISPR–Cas systems in other microorganisms for genome editing and metabolic engineering practices.

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.

ICEberg 2.0: an updated database of bacterial integrative and conjugative elements

ICEberg 2.0 (http://db-mml.sjtu.edu.cn/ICEberg/) is an updated database that provides comprehensive information about bacterial integrative and conjugative elements (ICEs). Compared with the previous version, three major improvements were made. First, with the aid of text mining and manual curation, it now recorded the details of 1032 ICEs, including 270 with experimental supports and 762 from bioinformatics prediction. Second, as increasing evidence has shown that ICEs frequently mobilize the so-called ‘hitchhikers’, such as integrative and mobilizable elements (IMEs) and cis-mobilizable elements (CIMEs), 83 known transfer interactions between 49 IMEs and 7 CIMEs with 19 ICEs taken from the literature were included and illustrated with visually intuitive directed graphs. An expanded collection of 260 chromosome-borne IMEs and 235 CIMEs was also added. At last, ICEberg 2.0 provides an online tool ICEfinder to predict ICEs or IMEs in bacterial genome sequences. It combines a similarity search for the integrase, relaxase and/or type IV secretion system and the co-localization of these corresponding homologous genes. With the recent updates, ICEberg 2.0 might provide better support for understanding the biological traits of ICEs, especially as their interaction with cognate mobilizable elements may further promote horizontal gene flow.

Challenges and Advances in Genome Editing Technologies in Streptomyces

The genome of Streptomyces encodes a high number of natural product (NP) biosynthetic gene clusters (BGCs). Most of these BGCs are not expressed or are poorly expressed (commonly called silent BGCs) under traditional laboratory experimental conditions. These NP BGCs represent an unexplored rich reservoir of natural compounds, which can be used to discover novel chemical compounds. To activate silent BGCs for NP discovery, two main strategies, including the induction of BGCs expression in native hosts and heterologous expression of BGCs in surrogate Streptomyces hosts, have been adopted, which normally requires genetic manipulation. So far, various genome editing technologies have been developed, which has markedly facilitated the activation of BGCs and NP overproduction in their native hosts, as well as in heterologous Streptomyces hosts. In this review, we summarize the challenges and recent advances in genome editing tools for Streptomyces genetic manipulation with a focus on editing tools based on clustered regularly interspaced short palindrome repeat (CRISPR)/CRISPR-associated protein (Cas) systems. Additionally, we discuss the future research focus, especially the development of endogenous CRISPR/Cas-based genome editing technologies in Streptomyces.

Natural diversity of CRISPR spacers of Thermus: evidence of local spacer acquisition and global spacer exchange

We investigated the diversity of CRISPR spacers of Thermus communities from two locations in Italy, two in Chile and one location in Russia. Among the five sampling sites, a total of more than 7200 unique spacers belonging to different CRISPR-Cas systems types and subtypes were identified. Most of these spacers are not found in CRISPR arrays of sequenced Thermus strains. Comparison of spacer sets revealed that samples within the same area (separated by few to hundreds of metres) have similar spacer sets, which appear to be largely stable at least over the course of several years. While at further distances (hundreds of kilometres and more) the similarity of spacer sets is decreased, there are still multiple common spacers in Thermus communities from different continents. The common spacers can be reconstructed in identical or similar CRISPR arrays, excluding their independent appearance and suggesting an extensive migration of thermophilic bacteria over long distances. Several new Thermus phages were isolated in the sampling sites. Mapping of spacers to bacteriophage sequences revealed examples of local acquisition of spacers from some phages and distinct patterns of targeting of phage genomes by different CRISPR-Cas systems.

This article is part of a discussion meeting issue ‘The ecology and evolution of prokaryotic CRISPR-Cas adaptive immune systems’.

Repurposing the Streptococcus mutans CRISPR-Cas9 System to Understand Essential Gene Function

A recent genome-wide screen identified ~300 essential or growth-supporting genes in the dental caries pathogen Streptococcus mutans. To be able to study these genes, we built a CRISPR interference tool around the Cas9 nuclease (Cas9Smu) encoded in the S. mutans UA159 genome. Using a xylose-inducible dead Cas9Smu with a constitutively active single-guide RNA (sgRNA), we observed titratable repression of GFP fluorescence that compared favorably to that of Streptococcus pyogenes dCas9 (Cas9Spy). We then investigated sgRNA specificity and proto-spacer adjacent motif (PAM) requirements. Interference by sgRNAs did not occur with double or triple base-pair mutations, or if single base-pair mutations were in the 3' end of the sgRNA. Bioinformatic analysis of >450 S. mutans genomes allied with in vivo assays revealed a similar PAM recognition sequence as Cas9Spy. Next, we created a comprehensive library of sgRNA plasmids that were directed at essential and growth-supporting genes. We discovered growth defects for 77% of the CRISPRi strains expressing sgRNAs. Phenotypes of CRISPRi strains, across several biological pathways, were assessed using fluorescence microscopy. A variety of cell structure anomalies were observed, including segregational instability of the chromosome, enlarged cells, and ovococci-to-rod shape transitions. CRISPRi was also employed to observe how silencing of cell wall glycopolysaccharide biosynthesis (rhamnose-glucose polysaccharide, RGP) affected both cell division and pathogenesis in a wax worm model. The CRISPRi tool and sgRNA library are valuable resources for characterizing essential genes in S. mutans, some of which could prove to be promising therapeutic targets.

Comparative Genomics of Pediococcus pentosaceus Isolated From Different Niches Reveals Genetic Diversity in Carbohydrate Metabolism and Immune System

Pediococcus pentosaceus isolated from fermented food and the gastrointestinal tracts of humans and animals have been widely identified, and some strains have been reported to reduce inflammation, encephalopathy, obesity and fatty liver in animals. In this study, the genomes of 65 P. pentosaceus strains isolated from human and animal feces and different fermented food were sequenced and comparative genomics analysis was performed on all strains along with nine sequenced representative strains to preliminarily reveal the lifestyle of P. pentosaceus, and investigate the genomic diversity within this species. The results reveal that P. pentosaceus is not host-specific, and shares core genes encoding proteins related to translation, ribosomal structure and biogenesis and signal transduction mechanisms, while its genetic diversity relates mainly to carbohydrate metabolism, and horizontally transferred DNA, especially prophages and bacteriocins encoded on plasmids. Additionally, this is the first report of a type IIA CRISPR/Cas system in P. pentosaceus. This work provides expanded resources of P. pentosaceus genomes, and offers a framework for understanding the biotechnological potential of this species.

Genomic and ecological attributes of marine bacteriophages encoding bacterial virulence genes

Background

Bacteriophages encode genes that modify bacterial functions during infection. The acquisition of phage-encoded virulence genes is a major mechanism for the rise of bacterial pathogens. In coral reefs, high bacterial density and lysogeny has been proposed to exacerbate reef decline through the transfer of phage-encoded virulence genes. However, the functions and distribution of these genes in phage virions on the reef remain unknown.

Results

Here, over 28,000 assembled viral genomes from the free viral community in Atlantic and Pacific Ocean coral reefs were queried against a curated database of virulence genes. The diversity of virulence genes encoded in the viral genomes was tested for relationships with host taxonomy and bacterial density in the environment. These analyses showed that bacterial density predicted the profile of virulence genes encoded by phages. The Shannon diversity of virulence-encoding phages was negatively related with bacterial density, leading to dominance of fewer genes at high bacterial abundances. A statistical learning analysis showed that reefs with high microbial density were enriched in viruses encoding genes enabling bacterial recognition and invasion of metazoan epithelium. Over 60% of phages could not have their hosts identified due to limitations of host prediction tools; for those which hosts were identified, host taxonomy was not an indicator of the presence of virulence genes.

Conclusions

This study described bacterial virulence factors encoded in the genomes of bacteriophages at the community level. The results showed that the increase in microbial densities that occurs during coral reef degradation is associated with a change in the genomic repertoire of bacteriophages, specifically in the diversity and distribution of bacterial virulence genes. This suggests that phages are implicated in the rise of pathogens in disturbed marine ecosystems.

CRISPR-Cas bioinformatics

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

Efficient cleavage resolves PAM preferences of CRISPR-Cas in human cells

Clustered regularly interspaced short palindromic repeats and associated proteins (CRISPR-Cas) of bacterial adaptive immunity have been adopted as a powerful and versatile tool for manipulation of the genome. This paradigm has been widely applied in biological research and treatments of animal or cellular disease models. A critical feature of CRISPR-Cas is the protospacer adjacent motif (PAM), which dictates the DNA target recognition mechanism of Cas proteins. While, direct identifying functional PAM sequences in human cells remains a challenge. Here, we developed a positive screen system termed PAM-DOSE (PAM Definition by Observable Sequence Excision) to delineate the functional PAMs in human cells. Specifically, the PAM libraries for CRISPR-Cas (SpCas9, SpCas9-NG, FnCas12a, AsCas12a, LbCas12a and MbCas12a) were generated and the corresponding CRISPR-Cas mediated cleaved fragments with functional PAM in human cells were harvested for DNA sequencing, which could be tracked and visualized with either florescence microscopy or flow cytometry analysis. With this system, we identified the functional PAMs of CRISPR-Cas members. We also found that spacer sequence affects the PAM preference of Cas proteins. This method will facilitate identification of functional PAMs for Cas-mediated human genome editing applications.

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.

Diversity of CRISPR/Cas system in Clostridium perfringens

Clostridium perfringens is an important pathogen of human and livestock infections, posing a threat to health. The horizontal gene transfer (HGT) of plasmids that carry toxin-related genes is involved in C. perfringens pathogenicity. The CRISPR/Cas system, which has been identified in a wide range of prokaryotes, provides acquired immunity against HGT. However, information about the CRISPR/Cas system in Clostridium perfringens is still limited. In this study, 111 C. perfringens strains with publicly available genomes were used to analyze the occurrence and diversity of CRISPR/Cas system and evaluate the potential of CRISPR-based genotyping in this multi-host pathogen. A total of 59 out of the 111 genomes harbored at least one confirmed CRISPR array. Four CRISPR/Cas system subtypes, including subtypes IB, IIA, IIC, and IIID systems, were identified in 32 strains. Subtype IB system was the most prevalent in this species, which was subdivided into four subgroups displaying subgroup specificity in terms of cas gene content, repeat sequence content, and PAM. We showed that the CRISPR spacer polymorphism can be used for evolutionary studies, and that it can provide discriminatory power for typing strains. Nevertheless, the application of this approach was largely limited to strains that contain the CRISPR/Cas system. Spacer origin analysis revealed that approximately one-fifth of spacers showed significant matches to plasmids and phages, thereby suggesting the implication of CRISPR/Cas systems in controlling HGT. Collectively, our results provide new insights into the diversity and evolution of CRISPR/Cas system in C. perfringens.

Characterization of CRISPR-Cas Systems in Serratia marcescens Isolated from Rhynchophorus ferrugineus (Olivier, 1790) (Coleoptera: Curculionidae)

The CRISPR-Cas adaptive immune system has been attracting increasing scientific interest for biological functions and biotechnological applications. Data on the Serratia marcescens system are scarce. Here, we report a comprehensive characterisation of CRISPR-Cas systems identified in S. marcescens strains isolated as secondary symbionts of Rhynchophorus ferrugineus, also known as Red Palm Weevil (RPW), one of the most invasive pests of major cultivated palms. Whole genome sequencing was performed on four strains (S1, S5, S8, and S13), which were isolated from the reproductive apparatus of RPWs. Subtypes I-F and I-E were harboured by S5 and S8, respectively. No CRISPR-Cas system was detected in S1 or S13. Two CRISPR arrays (4 and 51 spacers) were detected in S5 and three arrays (11, 31, and 30 spacers) were detected in S8. The CRISPR-Cas systems were located in the genomic region spanning from ybhR to phnP, as if this were the only region where CRISPR-Cas loci were acquired. This was confirmed by analyzing the S. marcescens complete genomes available in the NCBI database. This region defines a genomic hotspot for horizontally acquired genes and/or CRISPR-Cas systems. This study also supplies the first identification of subtype I-E in S. marcescens.

CRISPR-Cas System of a Prevalent Human Gut Bacterium Reveals Hyper-targeting against Phages in a Human Virome Catalog

Bacteriophages are abundant within the human gastrointestinal tract, yet their interactions with gut bacteria remain poorly understood, particularly with respect to CRISPR-Cas immunity. Here, we show that the type I-C CRISPR-Cas system in the prevalent gut Actinobacterium Eggerthella lenta is transcribed and sufficient for specific targeting of foreign and chromosomal DNA. Comparative analyses of E. lenta CRISPR-Cas systems across (meta)genomes revealed 2 distinct clades according to cas sequence similarity and spacer content. We assembled a human virome database (HuVirDB), encompassing 1,831 samples enriched for viral DNA, to identify protospacers. This revealed matches for a majority of spacers, a marked increase over other databases, and uncovered "hyper-targeted" phage sequences containing multiple protospacers targeted by several E. lenta strains. Finally, we determined the positional mismatch tolerance of observed spacer-protospacer pairs. This work emphasizes the utility of merging computational and experimental approaches for determining the function and targets of CRISPR-Cas systems.

Genome editing using the endogenous type I CRISPR-Cas system in Lactobacillus crispatus

CRISPR-Cas systems are now widely used for genome editing and transcriptional regulation in diverse organisms. The compact and portable nature of class 2 single effector nucleases, such as Cas9 or Cas12, has facilitated directed genome modifications in plants, animals, and microbes. However, most CRISPR-Cas systems belong to the more prevalent class 1 category, which hinges on multiprotein effector complexes. In the present study, we detail how the native type I-E CRISPR-Cas system, with a 5'-AAA-3' protospacer adjacent motif (PAM) and a 61-nucleotide guide CRISPR RNA (crRNA) can be repurposed for efficient chromosomal targeting and genome editing in Lactobacillus crispatus, an important commensal and beneficial microbe in the vaginal and intestinal tracts. Specifically, we generated diverse mutations encompassing a 643-base pair (bp) deletion (100% efficiency), a stop codon insertion (36%), and a single nucleotide substitution (19%) in the exopolysaccharide priming-glycosyl transferase (p-gtf). Additional genetic targets included a 308-bp deletion (20%) in the prophage DNA packaging Nu1 and a 730-bp insertion of the green fluorescent protein gene downstream of enolase (23%). This approach enables flexible alteration of the formerly genetically recalcitrant species L. crispatus, with potential for probiotic enhancement, biotherapeutic engineering, and mucosal vaccine delivery. These results also provide a framework for repurposing endogenous CRISPR-Cas systems for flexible genome targeting and editing, while expanding the toolbox to include one of the most abundant and diverse systems found in nature.

Multifactorial control of the expression of a CRISPR-Cas system by an extracytoplasmic function σ/anti-σ pair and a global regulatory complex

Expression of CRISPR-Cas systems is a prerequisite for their defensive role against invading genetic elements. Yet, much remains unknown about how this crucial step is regulated. We describe a new mechanism controlling CRISPR-cas expression, which requires an extracytoplasmic function (ECF) σ factor (DdvS), its membrane-bound anti-σ (DdvA) and a global regulatory complex (CarD–CarG). Transcriptomic analyses revealed that the DdvS/CarD/CarG-dependent regulon comprises a type III-B CRISPR-Cas system in Myxococcus xanthus. We mapped four DdvS-driven CarD/CarG-dependent promoters, with one lying immediately upstream of the cas cluster. Consistent with direct action, DdvS and CarD–CarG localize at these promoters in vivo. The cas genes are transcribed as a polycistronic mRNA that reads through the leader into the CRISPR array, a putative σ^A-dependent promoter in the leader having negligible activity in vivo. Consequently, expression of the entire CRISPR-Cas system and mature CRISPR-RNA (crRNA) production is DdvS/CarD/CarG-dependent. DdvA likely uses its large C-terminal domain to sense and transduce the extracytoplasmic signal triggering CRISPR-cas expression, which we show is not starvation-induced multicellular development. An ECF-σ/anti-σ pair and a global regulatory complex provide an effective mechanism to coordinate signal-sensing with production of precursor crRNA, its processing Cas6 endoribonuclease and other Cas proteins for mature crRNA biogenesis and interference.

Long reads reveal the diversification and dynamics of CRISPR reservoir in microbiomes

BACKGROUND

Sequencing of microbiomes has accelerated the characterization of the diversity of CRISPR-Cas immune systems. However, the utilization of next generation short read sequences for the characterization of CRISPR-Cas dynamics remains limited due to the repetitive nature of CRISPR arrays. CRISPR arrays are comprised of short spacer segments (derived from invaders' genomes) interspaced between flanking repeat sequences. The repetitive structure of CRISPR arrays poses a computational challenge for the accurate assembly of CRISPR arrays from short reads. In this paper we evaluate the use of long read sequences for the analysis of CRISPR-Cas system dynamics in microbiomes.

RESULTS

We analyzed a dataset of Illumina's TruSeq Synthetic Long-Reads (SLR) derived from a gut microbiome. We showed that long reads captured CRISPR spacers at a high degree of redundancy, which highlights the spacer conservation of spacer sharing CRISPR variants, enabling the study of CRISPR array dynamics in ways difficult to achieve though short read sequences. We introduce compressed spacer graphs, a visual abstraction of spacer sharing CRISPR arrays, to provide a simplified view of complex organizational structures present within CRISPR array dynamics. Utilizing compressed spacer graphs, several key defining characteristics of CRISPR-Cas system dynamics were observed including spacer acquisition and loss events, conservation of the trailer end spacers, and CRISPR arrays' directionality (transcription orientation). Other result highlights include the observation of intense array contraction and expansion events, and reconstruction of a full-length genome for a potential invader (Faecalibacterium phage) based on identified spacers.

CONCLUSION

We demonstrate in an in silico system that long reads provide the necessary context for characterizing the organization of CRISPR arrays in a microbiome, and reveal dynamic and evolutionary features of CRISPR-Cas systems in a microbial population.

Transposon-encoded CRISPR-Cas systems direct RNA-guided DNA integration

Conventional CRISPR-Cas systems maintain genomic integrity by leveraging guide RNAs for the nuclease-dependent degradation of mobile genetic elements, including plasmids and viruses. Here we describe a notable inversion of this paradigm, in which bacterial Tn7-like transposons have co-opted nuclease-deficient CRISPR-Cas systems to catalyse RNA-guided integration of mobile genetic elements into the genome. Programmable transposition of Vibrio cholerae Tn6677 in Escherichia coli requires CRISPR- and transposon-associated molecular machineries, including a co-complex between the DNA-targeting complex Cascade and the transposition protein TniQ. Integration of donor DNA occurs in one of two possible orientations at a fixed distance downstream of target DNA sequences, and can accommodate variable length genetic payloads. Deep-sequencing experiments reveal highly specific, genome-wide DNA insertion across dozens of unique target sites. This discovery of a fully programmable, RNA-guided integrase lays the foundation for genomic manipulations that obviate the requirements for double-strand breaks and homology-directed repair.

Structural basis of Type IV CRISPR RNA biogenesis by a Cas6 endoribonuclease

Prokaryotic CRISPR-Cas adaptive immune systems rely on small non-coding RNAs derived from CRISPR loci to recognize and destroy complementary nucleic acids. However, the mechanism of Type IV CRISPR RNA (crRNA) biogenesis is poorly understood. To dissect the mechanism of Type IV CRISPR RNA biogenesis, we determined the x-ray crystal structure of the putative Type IV CRISPR associated endoribonuclease Cas6 from Mahella australiensis (Ma Cas6-IV) and characterized its enzymatic activity with RNA cleavage assays. We show that Ma Cas6-IV specifically cleaves Type IV crRNA repeats at the 3' side of a predicted stem loop, with a metal-independent, single-turnover mechanism that relies on a histidine and a tyrosine located within the putative endonuclease active site. Structure and sequence alignments with Cas6 orthologs reveal that although Ma Cas6-IV shares little sequence homology with other Cas6 proteins, all share common structural features that bind distinct crRNA repeat sequences. This analysis of Type IV crRNA biogenesis provides a structural and biochemical framework for understanding the similarities and differences of crRNA biogenesis across multi-subunit Class 1 CRISPR immune systems.

Genome sequencing and analysis of Salmonella enterica subsp. enterica serovar Stanley UPM 517: Insights on its virulence-associated elements and their potentials as vaccine candidates

Salmonella enterica subsp. enterica serovar Stanley (S. Stanley) is a pathogen that contaminates food, and is related to Salmonella outbreaks in a variety of hosts such as humans and farm animals through products like dairy items and vegetables. Despite the fact that several vaccines of Salmonella strains had been constructed, none of them were developed according to serovar Stanley up to this day. This study presents results of genome sequencing and analysis on our S. Stanley UPM 517 strain taken from fecal swabs of 21-day-old healthy commercial chickens in Perak, Malaysia and used Salmonella enterica subsp. enterica serovar Typhimurium LT2 (S. Typhimurium LT2) as a reference to be compared with. First, sequencing and assembling of the Salmonella Stanley UPM 517 genome into a contiguous form were done. The work was then continued with scaffolding and gap filling. Annotation and alignment of the draft genome was performed with S. Typhimurium LT2. The other elements of virulence estimated in this study included Salmonella pathogenicity islands, resistance genes, prophages, virulence factors, plasmid regions, restriction-modification sites and the CRISPR-Cas system. The S. Stanley UPM 517 draft genome had a length of 4,736,817 bp with 4,730 coding sequence and 58 RNAs. It was discovered via genomic analysis on this strain that there were antimicrobial resistance properties toward a wide variety of antibiotics. Tcf and ste, the two fimbrial virulence clusters related with human and broiler intestinal colonizations which were not found in S. Typhimurium LT2, were atypically discovered in the S. Stanley UPM 517 genome. These clusters are involved in the intestinal colonization of human and broilers, respectively. There were seven Salmonella pathogenicity islands (SPIs) within the draft genome, which contained the virulence factors associated with Salmonella infection (except SPI-14). Five intact prophage regions, mostly comprising of the protein encoding Gifsy-1, Fels-1, RE-2010 and SEN34 prophages, were also encoded in the draft genome. Also identified were Type I–III restriction-modification sites and the CRISPR-Cas system of the Type I–E subtype. As this strain exhibited resistance toward numerous antibiotics, we distinguished several genes that had the potential for removal in the construction of a possible vaccine candidate to restrain and lessen the pervasiveness of salmonellosis and to function as an alternative to antibiotics.

Long-read based de novo assembly of low-complexity metagenome samples results in finished genomes and reveals insights into strain diversity and an active phage system

BACKGROUND

Complete and contiguous genome assemblies greatly improve the quality of subsequent systems-wide functional profiling studies and the ability to gain novel biological insights. While a de novo genome assembly of an isolated bacterial strain is in most cases straightforward, more informative data about co-existing bacteria as well as synergistic and antagonistic effects can be obtained from a direct analysis of microbial communities. However, the complexity of metagenomic samples represents a major challenge. While third generation sequencing technologies have been suggested to enable finished metagenome-assembled genomes, to our knowledge, the complete genome assembly of all dominant strains in a microbiome sample has not been demonstrated. Natural whey starter cultures (NWCs) are used in cheese production and represent low-complexity microbiomes. Previous studies of Swiss Gruyère and selected Italian hard cheeses, mostly based on amplicon metagenomics, concurred that three species generally pre-dominate: Streptococcus thermophilus, Lactobacillus helveticus and Lactobacillus delbrueckii.

RESULTS

Two NWCs from Swiss Gruyère producers were subjected to whole metagenome shotgun sequencing using the Pacific Biosciences Sequel and Illumina MiSeq platforms. In addition, longer Oxford Nanopore Technologies MinION reads had to be generated for one to resolve repeat regions. Thereby, we achieved the complete assembly of all dominant bacterial genomes from these low-complexity NWCs, which was corroborated by a 16S rRNA amplicon survey. Moreover, two distinct L. helveticus strains were successfully co-assembled from the same sample. Besides bacterial chromosomes, we could also assemble several bacterial plasmids and phages and a corresponding prophage. Biologically relevant insights were uncovered by linking the plasmids and phages to their respective host genomes using DNA methylation motifs on the plasmids and by matching prokaryotic CRISPR spacers with the corresponding protospacers on the phages. These results could only be achieved by employing long-read sequencing data able to span intragenomic as well as intergenomic repeats.

CONCLUSIONS

Here, we demonstrate the feasibility of complete de novo genome assembly of all dominant strains from low-complexity NWCs based on whole metagenomics shotgun sequencing data. This allowed to gain novel biological insights and is a fundamental basis for subsequent systems-wide omics analyses, functional profiling and phenotype to genotype analysis of specific microbial communities.

Prediction of activity and specificity of CRISPR-Cpf1 using convolutional deep learning neural networks

BACKGROUND

CRISPR-Cpf1 has recently been reported as another RNA-guided endonuclease of class 2 CRISPR-Cas system, which expands the molecular biology toolkit for genome editing. However, most of the online tools and applications to date have been developed primarily for the Cas9. There are a limited number of tools available for the Cpf1.

RESULTS

We present DeepCpf1, a deep convolution neural networks (CNN) approach to predict Cpf1 guide RNAs on-target activity and off-target effects using their matched and mismatched DNA sequences. Trained on published data sets, DeepCpf1 is superior to other machine learning algorithms and reliably predicts the most efficient and less off-target effects guide RNAs for a given gene. Combined with a permutation importance analysis, the key features of guide RNA sequences are identified, which determine the activity and specificity of genome editing.

CONCLUSIONS

DeepCpf1 can significantly improve the accuracy of Cpf1-based genome editing and facilitates the generation of optimized guide RNAs libraries.

A phylogenetic test of the role of CRISPR-Cas in limiting plasmid acquisition and prophage integration in bacteria

CRISPR-Cas is a prokaryotic defense system capable of protecting the cell from damaging foreign genetic elements. However, some genetic elements can be beneficial, which suggests the hypothesis that bacteria with CRISPR-Cas incur a cost of reduced intake of mutualistic plasmids and prophage. Here we present the first robust test of this hypothesis that controls for phylogenic and ecological biases in the distribution of CRISPR-Cas. We filtered the available genomic data (~7000 strains from ~2100 species) by first selecting all pairs of conspecific strains, one with and one without CRISPR-Cas (controlling ecological bias), and second by retaining only one such pair per bacterial family (controlling phylogenetic bias), resulting in pairs representing 38 bacterial families. Analysis of these pairs of bacterial strains showed that on average the CRISPR-Cas strain of each pair contained significantly fewer plasmids than its CRISPR-Cas negative partner (0.86 vs. 1.86). It also showed that the CRISPR-Cas positive strains had 31% fewer intact prophage (1.17 vs. 1.75), but the effect was highly variable and not significant. These results support the hypothesis that CRISPR-Cas reduces the rate of plasmid-mediated HGT and, given the abundant evidence of beneficial genes carried by plasmids, provide a clear example of a cost associated with the CRISPR-Cas system.

Analysis of CRISPR/Cas system of Proteus and the factors affected the functional mechanism

BACKGROUND

The Proteus is one of the most common human and animal pathogens. Clustered regularly interspaced short palindromic repeats and CRISPR-associated proteins (CRISPR/Cas) are inheritable genetic elements found in a variety of archaea and bacteria in the evolution, providing immune function against foreign invasion.

OBJECTIVES

To analyze the characteristics and functions of the CRISPR/Cas system in Proteus genomes, as well as the internal and external factors affecting the system.

METHODS

CRISPR loci were identified and divided into groups based on the repeat sequence in 96 Proteus strains by identification. Compared the RNA secondary structure and minimum free energy of CRISPR loci through bioinformatics, the evolution of cas genes, and the effects of related elements were also discussed.

RESULTS

85 CRISPR loci were identified and divided into six groups based on the sequence of repeats, and the more stable the secondary structure of RNA, the smaller the minimum free energy, the fewer base mutations in the repeat, the more stable the CRISPR and the more complete the evolution of the system. In addition, Cas1 gene can be a symbol to distinguish species to some extent. Of all the influencing factors, CRISPR/Cas had the greatest impact on plasmids.

CONCLUSIONS

This study examined the diversity of CRISPR/Cas system in Proteus and found statistically significant positive/negative correlations between variety factors (the RNA stability, free energy, etc.) and the CRISPR locus, which played a vital role in regulating the CRISPR/Cas system.