1型CRISPR位点间区序列分析在嗜热链球菌菌株分型中的应用

为了建立适用于嗜热链球菌菌株资源多样性调查的菌株分型方法,尝试将1型CRISPR位点间区序列分析用于嗜热链球菌的菌株分型,并与常用ERIC-PCR指纹图谱方法进行了比较。结果表明,1型CRISPR位点间区序列分析可以把30株从三个不同样品中分离的嗜热链球菌分成三种差异明显的类型:不同类型菌株之间没有相同的间区序列;而同一类型菌株之间,间区序列的组成和排列则基本一致,并且上述分型的结果与用ERIC-PCR指纹图谱技术获得的结果完全一致。因此,1型CRISPR位点间区序列分析是嗜热链球菌分型鉴定的可靠方法,并适用于大量菌株的分型鉴定和多样性调查。     

Diversity, activity, and evolution of CRISPR loci in Streptococcus thermophilus

Clustered regularly interspaced short palindromic repeats (CRISPR) are hypervariable loci widely distributed in prokaryotes that provide acquired immunity against foreign genetic elements. Here, we characterize a novel Streptococcus thermophilus locus, CRISPR3, and experimentally demonstrate its ability to integrate novel spacers in response to bacteriophage. Also, we analyze CRISPR diversity and activity across three distinct CRISPR loci in several S. thermophilus strains. We show that both CRISPR repeats and cas genes are locus specific and functionally coupled. A total of 124 strains were studied, and 109 unique spacer arrangements were observed across the three CRISPR loci. Overall, 3,626 spacers were analyzed, including 2,829 for CRISPR1 (782 unique), 173 for CRISPR2 (16 unique), and 624 for CRISPR3 (154 unique). Sequence analysis of the spacers revealed homology and identity to phage sequences (77%), plasmid sequences (16%), and S. thermophilus chromosomal sequences (7%). Polymorphisms were observed for the CRISPR repeats, CRISPR spacers, cas genes, CRISPR motif, locus architecture, and specific sequence content. Interestingly, CRISPR loci evolved both via polarized addition of novel spacers after exposure to foreign genetic elements and via internal deletion of spacers. We hypothesize that the level of diversity is correlated with relative CRISPR activity and propose that the activity is highest for CRISPR1, followed by CRISPR3, while CRISPR2 may be degenerate. Globally, the dynamic nature of CRISPR loci might prove valuable for typing and comparative analyses of strains and microbial populations. Also, CRISPRs provide critical insights into the relationships between prokaryotes and their environments, notably the coevolution of host and viral genomes.     

Intricate interactions between the bloom-forming cyanobacterium Microcystis aeruginosa and foreign genetic elements, revealed by diversified clustered regularly interspaced short palindromic repeat (CRISPR) signatures

Clustered regularly interspaced short palindromic repeats (CRISPR) confer sequence-dependent, adaptive resistance in prokaryotes against viruses and plasmids via incorporation of short sequences, called spacers, derived from foreign genetic elements. CRISPR loci are thus considered to provide records of past infections. To describe the host-parasite (i.e., cyanophages and plasmids) interactions involving the bloom-forming freshwater cyanobacterium Microcystis aeruginosa, we investigated CRISPR in four M. aeruginosa strains and in two previously sequenced genomes. The number of spacers in each locus was larger than the average among prokaryotes. All spacers were strain specific, except for a string of 11 spacers shared in two closely related strains, suggesting diversification of the loci. Using CRISPR repeat-based PCR, 24 CRISPR genotypes were identified in a natural cyanobacterial community. Among 995 unique spacers obtained, only 10 sequences showed similarity to M. aeruginosa phage Ma-LMM01. Of these, six spacers showed only silent or conservative nucleotide mutations compared to Ma-LMM01 sequences, suggesting a strategy by the cyanophage to avert CRISPR immunity dependent on nucleotide identity. These results imply that host-phage interactions can be divided into M. aeruginosa-cyanophage combinations rather than pandemics of population-wide infectious cyanophages. Spacer similarity also showed frequent exposure of M. aeruginosa to small cryptic plasmids that were observed only in a few strains. Thus, the diversification of CRISPR implies that M. aeruginosa has been challenged by diverse communities (almost entirely uncharacterized) of cyanophages and plasmids.     

Comparative Analysis of the Orphan CRISPR2 Locus in 242 Enterococcus faecalis Strains

Clustered, Regularly Interspaced Short Palindromic Repeats and their associated Cas proteins (CRISPR-Cas) provide prokaryotes with a mechanism for defense against mobile genetic elements (MGEs). A CRISPR locus is a molecular memory of MGE encounters. It contains an array of short sequences, called spacers, that generally have sequence identity to MGEs. Three different CRISPR loci have been identified among strains of the opportunistic pathogen Enterococcus faecalis. CRISPR1 and CRISPR3 are associated with the cas genes necessary for blocking MGEs, but these loci are present in only a subset of E. faecalis strains. The orphan CRISPR2 lacks cas genes and is ubiquitous in E. faecalis, although its spacer content varies from strain to strain. Because CRISPR2 is a variable locus occurring in all E. faecalis, comparative analysis of CRISPR2 sequences may provide information about the clonality of E. faecalis strains. We examined CRISPR2 sequences from 228 E. faecalis genomes in relationship to subspecies phylogenetic lineages (sequence types; STs) determined by multilocus sequence typing (MLST), and to a genome phylogeny generated for a representative 71 genomes. We found that specific CRISPR2 sequences are associated with specific STs and with specific branches on the genome tree. To explore possible applications of CRISPR2 analysis, we evaluated 14 E. faecalis bloodstream isolates using CRISPR2 analysis and MLST. CRISPR2 analysis identified two groups of clonal strains among the 14 isolates, an assessment that was confirmed by MLST. CRISPR2 analysis was also used to accurately predict the ST of a subset of isolates. We conclude that CRISPR2 analysis, while not a replacement for MLST, is an inexpensive method to assess clonality among E. faecalis isolates, and can be used in conjunction with MLST to identify recombination events occurring between STs.     

Polymorphism at the ovine major histocompatibility complex class II loci

Southern hybridization analysis of the ovine major histocompatibility complex (MHC) ( MhcOvar ) class II region, using sheep-specific probes for the DQA1, DQA2, DQB and DRA loci, has revealed extensive polymorphism. DQA1 and DQAP had eight and 16 alleles respectively, DQB had six and DRA had three alleles. Little information was derived from the DRB locus owing to extensive cross-hybridization between the DRB probe and the DQB locus. Differences in allele frequency between breeds were revealed. At the DQA1 locus a null allele (DQA1-N) was observed with a frequency of between 27% and 45%, making this the most common DQA1 allele in all breeds examined. The frequency of DQA1-N homozygotes was between 11% and 18%, raising questions as to the functional significance of the DQA1 gene. Linkage analysis between the DQA1, DQA2, DQB and DRA loci did not reveal any recombination.     

A data-mining approach to spacer oligonucleotide typing of Mycobacterium tuberculosis

MOTIVATION: The Direct Repeat (DR) locus of Mycobacterium tuberculosis is a suitable model to study (i) molecular epidemiology and (ii) the evolutionary genetics of tuberculosis. This is achieved by a DNA analysis technique (genotyping), called sp acer oligo nucleotide typing (spoligotyping ). In this paper, we investigated data analysis methods to discover intelligible knowledge rules from spoligotyping, that has not yet been applied on such representation. This processing was achieved by applying the C4.5 induction algorithm and knowledge rules were produced. Finally, a Prototype Selection (PS) procedure was applied to eliminate noisy data. This both simplified decision rules, as well as the number of spacers to be tested to solve classification tasks. In the second part of this paper, the contribution of 25 new additional spacers and the knowledge rules inferred were studied from a machine learning point of view. From a statistical point of view, the correlations between spacers were analyzed and suggested that both negative and positive ones may be related to potential structural constraints within the DR locus that may shape its evolution directly or indirectly. RESULTS: By generating knowledge rules induced from decision trees, it was shown that not only the expert knowledge may be modeled but also improved and simplified to solve automatic classification tasks on unknown patterns. A practical consequence of this study may be a simplification of the spoligotyping technique, resulting in a reduction of the experimental constraints and an increase in the number of samples processed.     

Ultra-high resolution HLA genotyping and allele discovery by highly multiplexed cDNA amplicon pyrosequencing

ABSTRACT: BACKGROUND: High-resolution HLA genotyping is a critical diagnostic and research assay. Current methods rarely achieve unambiguous high-resolution typing without making population-specific frequency inferences due to a lack of locus coverage and difficulty in exon-phase matching. Achieving high-resolution typing is also becoming more challenging with traditional methods as the database of known HLA alleles increases. RESULTS: We designed a cDNA amplicon-based pyrosequencing method to capture 94% of the HLA class I open-reading-frame with only two amplicons per sample, and an analogous method for class II HLA genes, with a primary focus on sequencing the DRB loci. We present a novel Galaxy server-based analysis workflow for determining genotype. During assay validation, we performed two GS Junior sequencing runs to determine the accuracy of the HLA class I amplicons and DRB amplicon at different levels of multiplexing. When 116 amplicons were multiplexed, we unambiguously resolved 99%of class I alleles to four- or six-digit resolution, as well as 100% unambiguous DRB calls. The second experiment, with 271 multiplexed amplicons, missed some alleles, but generated high-resolution, concordant typing for 93% of class I alleles, and 96% for DRB1 alleles. In a third, preliminary experiment we attempted to sequence novel amplicons for other class II loci with mixed success. CONCLUSIONS: The presented assay is higher-throughput and higher-resolution than existing HLA genotyping methods, and suitable for allele discovery or large cohort sampling. The validated class I and DRB primers successfully generated unambiguously high-resolution genotypes, while further work is needed to validate additional class II genotyping amplicons.     

Heterozygosity excess at the cattle DRB locus revealed by large scale genotyping of two closely linked microsatellites

A method for MHC DRB typing in cattle based on two closely linked and highly polymorphic microsatellites is described. The two microsatellites DRBP1ms and DRB3ms are located in intron 2 of the corresponding DRB gene. The very strong linkage disequilibrium between the two loci made it possible to establish DRB microsatellite haplotypes. The typing results with this method on reference samples followed closely that obtained with RFLP and direct sequence analysis of DRB3 exon 2. The method is well suited for large scale genotyping and was successfully applied for typing more than 600 unrelated animals representing 23 breeds. The data were used to test whether the observed DRB allele frequency distributions were consistent with that expected for selectively neutral alleles in populations at mutation-drift equilibrium. A significant heterozygosity excess was detected and there was an obvious trend across breeds towards a more even allele frequency distribution than expected. The deviation may be due to balancing selection acting on the DRB locus or by recent population bottlenecks.     

CRISPR analysis of bacteriophage‐insensitive mutants (BIMs) of industrial Streptococcus thermophilus– implications for starter design

Aims: An efficient approach for generation of bacteriophage‐insensitive mutants (BIMs) of Streptococcus thermophilus starters was described in our laboratory [Mills et al. (2007) J Microbiol Methods 70 , 159–164]. The aim of this study was to analyse the phage resistance mechanism responsible for BIM formation. Methods and Results: Three clustered regularly interspaced short palindromic repeat (CRISPR) regions have been identified in Strep. thermophilus, and Strep. thermophilus can integrate novel spacers into these loci in response to phage attack. Characterization of three sets of BIMs indicated that two sets had altered CRISPR1 and/or CRISPR3 loci. A range of BIMs of yoghurt starter CSK938 were generated with the same phage in different phage challenge experiments, and each acquired unique spacer regions ranging between one and four new spacers in CRISPR1. In addition, the BIM that acquired only one new spacer in CRISPR1 also acquired an additional spacer in CRISPR3. A fourth BIM, generated with a different phage, had two spacers deleted from CRISPR1 but acquired two spacers in CRISPR3. Analysis of the Mozzarella starter CSK939 and its associated BIMs indicated that formation of second generation BIMs does not lead to increases in spacer number but to alterations in spacer regions. BIMs of an exopolysaccharide (EPS)‐producing strain that lost the ability to produce EPS did not harbour an altered CRISPR, suggesting that phage sensitivity may be related to the EPS‐producing phenotype. Conclusions: Acquisition/deletion of new spacers in CRISPR loci in response to phage attack generates distinctly individual variants. It also demonstrates that other modifications may be responsible for the phage resistance of Strep. thermophilus BIMs. Significance and Impact of the Study: Isolation of individual BIMs that have unique spacers towards the leader region of the CRISPR locus may be a very useful approach for rotation strategies with the same starter backbone. Upon phage infection, BIMs ‘in reserve’ can be slotted into the rotation scheme.     

Genotyping of Genetically Monomorphic Bacteria: DNA Sequencing in Mycobacterium tuberculosis Highlights the Limitations of Current Methodologies

Because genetically monomorphic bacterial pathogens harbour little DNA sequence diversity, most current genotyping techniques used to study the epidemiology of these organisms are based on mobile or repetitive genetic elements. Molecular markers commonly used in these bacteria include Clustered Regulatory Short Palindromic Repeats (CRISPR) and Variable Number Tandem Repeats (VNTR). These methods are also increasingly being applied to phylogenetic and population genetic studies. Using the Mycobacterium tuberculosis complex (MTBC) as a model, we evaluated the phylogenetic accuracy of CRISPR- and VNTR-based genotyping, which in MTBC are known as spoligotyping and Mycobacterial Interspersed Repetitive Units (MIRU)-VNTR-typing, respectively. We used as a gold standard the complete DNA sequences of 89 coding genes from a global strain collection. Our results showed that phylogenetic trees derived from these multilocus sequence data were highly congruent and statistically robust, irrespective of the phylogenetic methods used. By contrast, corresponding phylogenies inferred from spoligotyping or 15-loci-MIRU-VNTR were incongruent with respect to the sequence-based trees. Although 24-loci-MIRU-VNTR performed better, it was still unable to detect all strain lineages. The DNA sequence data showed virtually no homoplasy, but the opposite was true for spoligotyping and MIRU-VNTR, which was consistent with high rates of convergent evolution and the low statistical support obtained for phylogenetic groupings defined by these markers. Our results also revealed that the discriminatory power of the standard 24 MIRU-VNTR loci varied by strain lineage. Taken together, our findings suggest strain lineages in MTBC should be defined based on phylogenetically robust markers such as single nucleotide polymorphisms or large sequence polymorphisms, and that for epidemiological purposes, MIRU-VNTR loci should be used in a lineage-dependent manner. Our findings have implications for strain typing in other genetically monomorphic bacteria.     

Molecular inversion probes for sensitive detection of Mycobacterium tuberculosis

Nucleic acid-based detection of Mycobacterium tuberculosis infections has the potential to improve the analysis of the tuberculosis epidemiology and patient care by increasing the specificity and sensitivity of diagnosis. One potential diagnostic sequence, the DR locus, is present in all isolates of M. tuberculosis complex bacteria. It encodes no known gene product but is useful for molecular typing of M. tuberculosis because of its fortuitous absence in non-tuberculosis strains of mycobacteria. The DR locus contains a variable number of short direct repeats interspersed with non-repetitive spacers and is commonly used as a target for the spoligotyping method, a technique based on the detection of the presence or absence of distinct spacers between the repeats. In this study, we attempted to combine the specificity of molecular inversion probe (MIP) technology with the sensitivity of modified pyrosequencing readout in order to detect a short conserved 18 bp sequence included in DR locus in 25 isolates of M. tuberculosis. Additional sensitivity was obtained by introducing modifications in pyrosequencing methodology, by these means we achieved to detect 500 fg of M. tuberculosis DNA.     

CRISPR inhibition of prophage acquisition in Streptococcus pyogenes

Streptococcus pyogenes, one of the major human pathogens, is a unique species since it has acquired diverse strain-specific virulence properties mainly through the acquisition of streptococcal prophages. In addition, S. pyogenes possesses clustered regularly interspaced short palindromic repeats (CRISPR)/Cas systems that can restrict horizontal gene transfer (HGT) including phage insertion. Therefore, it was of interest to examine the relationship between CRISPR and acquisition of prophages in S. pyogenes. Although two distinct CRISPR loci were found in S. pyogenes, some strains lacked CRISPR and these strains possess significantly more prophages than CRISPR harboring strains. We also found that the number of spacers of S. pyogenes CRISPR was less than for other streptococci. The demonstrated spacer contents, however, suggested that the CRISPR appear to limit phage insertions. In addition, we found a significant inverse correlation between the number of spacers and prophages in S. pyogenes. It was therefore suggested that S. pyogenes CRISPR have permitted phage insertion by lacking its own spacers. Interestingly, in two closely related S. pyogenes strains (SSI-1 and MGAS315), CRISPR activity appeared to be impaired following the insertion of phage genomes into the repeat sequences. Detailed analysis of this prophage insertion site suggested that MGAS315 is the ancestral strain of SSI-1. As a result of analysis of 35 additional streptococcal genomes, it was suggested that the influences of the CRISPR on the phage insertion vary among species even within the same genus. Our results suggested that limitations in CRISPR content could explain the characteristic acquisition of prophages and might contribute to strain-specific pathogenesis in S. pyogenes.     

Genetic variation and evolutionary origin of the direct repeat locus of Mycobacterium tuberculosis complex bacteria

The direct repeat region in Mycobacterium tuberculosis complex strains is composed of multiple direct variant repeats (DVRs), each of which is composed of a 36-bp direct repeat (DR) plus a nonrepetitive spacer sequence of similar size. It has been shown previously that clinical isolates show extensive polymorphism in the DR region by the variable presence of DVRs, and this polymorphism has been used in the epidemiology of tuberculosis. In an attempt to better understand the evolutionary scenario leading to polymorphic DR loci and to improve strain differentiation by spoligotyping, we characterized and compared the DNA sequences of the complete DR region and its flanking DNA of M. tuberculosis complex strains. We identified 94 different spacer sequences among 26 M. tuberculosis complex strains. No sequence homology was found between any of these spacers and M. tuberculosis DNA outside of the DR region or with any other known bacterial sequence. Although strains differed extensively in the presence or absence of DVRs, the order of the spacers in the DR locus was found to be well conserved. The data strongly suggest that the polymorphism in clinical isolates is the result of successive deletions of single discrete DVRs or of multiple contiguous DVRs from a primordial DR region containing many more DVRs than seen in present day isolates and that virtually no scrambling of DVRs took place during evolution. Because the majority of the novel spacer sequences identified in this study were confined to isolates of the rare Mycobacterium canettii taxon, the use of the novel spacers in spoligotyping led only to a slight improvement of strain differentiation by spoligotyping.     

Selective and hyperactive uptake of foreign DNA by adaptive immune systems of an archaeon via two distinct mechanisms

Central to the disparate adaptive immune systems of archaea and bacteria are clustered regularly interspaced short palindromic repeats (CRISPR). The spacer regions derive from invading genetic elements and, via RNA intermediates and associated proteins, target and cleave nucleic acids of the invader. Here we demonstrate the hyperactive uptake of hundreds of unique spacers within CRISPR loci associated with type I and IIIB immune systems of a hyperthermophilic archaeon. Infection with an environmental virus mixture resulted in the exclusive uptake of protospacers from a co-infecting putative conjugative plasmid. Spacer uptake occurred by two distinct mechanisms in only one of two CRISPR loci subfamilies present. In two loci, insertions, often multiple, occurred adjacent to the leader while in a third locus single spacers were incorporated throughout the array. Protospacer DNAs were excised from the invading genetic element immediately after CCN motifs, on either strand, with the secondary cut apparently produced by a ruler mechanism. Over a 10-week period, there was a gradual decrease in the number of wild-type cells present in the culture but the virus and putative conjugative plasmid were still propagating. The results underline the complex dynamics of CRISPR-based immune systems within a population infected with genetic elements.     

The forest behind the tree: phylogenetic exploration of a dominant Mycobacterium tuberculosis strain lineage from a high tuberculosis burden country

Background

Genotyping of Mycobacterium tuberculosis isolates is a powerful tool for epidemiological control of tuberculosis (TB) and phylogenetic exploration of the pathogen. Standardized PCR-based typing, based on 15 to 24 mycobacterial interspersed repetitive unit-variable number of tandem repeat (MIRU-VNTR) loci combined with spoligotyping, has been shown to have adequate resolution power for tracing TB transmission and to be useful for predicting diverse strain lineages in European settings. Its informative value needs to be tested in high TB-burden countries, where the use of genotyping is often complicated by dominance of geographically specific, genetically homogeneous strain lineages.

Methodology/Principal Findings

We tested this genotyping system for molecular epidemiological analysis of 369 M. tuberculosis isolates from 3 regions of Brazil, a high TB-burden country. Deligotyping, targeting 43 large sequence polymorphisms (LSPs), and the MIRU-VNTRplus identification database were used to assess phylogenetic predictions. High congruence between the different typing results consistently revealed the countrywide supremacy of the Latin-American-Mediterranean (LAM) lineage, comprised of three main branches. In addition to an already known RDRio branch, at least one other branch characterized by a phylogenetically informative LAM3 spoligo-signature seems to be globally distributed beyond Brazil. Nevertheless, by distinguishing 321 genotypes in this strain population, combined MIRU-VNTR typing and spoligotyping demonstrated the presence of multiple distinct clones. The use of 15 to 24 loci discriminated 21 to 25% more strains within the LAM lineage, compared to a restricted lineage-specific locus set suggested to be used after SNP analysis. Noteworthy, 23 of the 28 molecular clusters identified were exclusively composed of patient isolates from a same region, consistent with expected patterns of mostly local TB transmission.

Conclusions/Significance

Standard MIRU-VNTR typing combined with spoligotyping can reveal epidemiologically meaningful clonal diversity behind a dominant M. tuberculosis strain lineage in a high TB-burden country and is useful to explore international phylogenetical ramifications.     

High resolution genetic typing of the class II HLA-DRB1 locus using group-specific amplification and SSO-hybridisation in microplates

The HLA-DRB1 locus is one of the most polymorphic HLA class II loci and rapid and accurate typing of this polymorphism is important both in bone-marrow transplantation, analysis of disease association and in forensic medicine. The allelic variation at DRB1 is characterized by combinations of a limited number of amino-acid motifs, reducing the resolution of a typing strategy based on a single PCR and subsequent analysis of polymorphic motifs. In the present paper we describe a strategy for typing of DRB1 based on eight allele-specific PCRs followed by sandwich hybridization to immobilized probes in a microplate format. The combined approach results in a rapid typing system with very high resolution. Using a rapid DNA extraction protocol, a complete HLA-DRB1 typing can be performed in less than a day.     

Polymorphisms in MHC-DRA and -DRB alleles of water buffalo (Bubalus bubalis) reveal different features from cattle DR alleles

Seventy-five individuals of Bubalus bubalis belonging to four different breeds, three of river buffalo and one of swamp buffalo, were studied for polymorphism in MHC DRB (Bubu-DRB) and DRA (Bubu-DRA) loci. Eight alleles of Bubu-DRB were found, and all alleles in the swamp type were shared with the three river breeds. All alleles sampled from the breed of European origin (Mediterranean) were present in breeds sampled in Brazil, thus variability of this locus may have been preserved to a great extent in the more recently founded Brazilian population. Bubu-DRB alleles contained higher proportions of synonymous vs. non-synonymous substitutions in the non-peptide-binding sites (PBS) region, in contrast to the pattern of variation found in BoLA-DRB3, the orthologous locus in cattle. This indicated that either the first domain exon (exon 2) of Bubu-DRB has not undergone as much recombination and/or gene conversion as in cattle alleles, or Bubu-DRB may be more ancient than BoLA-DRB3 alleles. Phylogenetic analysis of DRB alleles from Bubalus, Syncerus c. caffer, the Cape buffalo, and domestic cattle demonstrated transspecies polymorphism. Water buffalo contained two alleles of DRA that differed from each other in two amino acid positions, including one in the PBS (alpha22) that was also shared with Anoa depressicornis, the anoa. Discovery of variation in DRA was surprising as the first domain of DRA is a highly conserved polypeptide in mammals in general and especially in ruminants, where no other substitution in PBS was seen.     

Structure and Genetic Content of the Megaplasmids of Neurotoxigenic Clostridium butyricum Type E Strains from Italy

We determined the genetic maps of the megaplasmids of six neutoroxigenic Clostridium butyricum type E strains from Italy using molecular and bioinformatics techniques. The megaplasmids are circular, not linear as we had previously proposed. The differently-sized megaplasmids share a genetic region that includes structural, metabolic and regulatory genes. In addition, we found that a 168 kb genetic region is present only in the larger megaplasmids of two tested strains, whereas it is absent from the smaller megaplasmids of the four remaining strains. The genetic region unique to the larger megaplasmids contains, among other features, a locus for clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR associated (cas) genes, i.e. a bacterial adaptive immune system providing sequence-specific protection from invading genetic elements. Some CRISPR spacer sequences of the neurotoxigenic C. butyricum type E strains showed homology to prophage, phage and plasmid sequences from closely related clostridia species or from distant species, all sharing the intestinal habitat, suggesting that the CRISPR locus might be involved in the microorganism adaptation to the human or animal intestinal environment. Besides, we report here that each of four distinct CRISPR spacers partially matched DNA sequences of different prophages and phages, at identical nucleotide locations. This suggests that, at least in neurotoxigenic C. butyricum type E, the CRISPR locus is potentially able to recognize the same conserved DNA sequence of different invading genetic elements, besides targeting sequences unique to previously encountered invading DNA, as currently predicted for a CRISPR locus. Thus, the results of this study introduce the possibility that CRISPR loci can provide resistance to a wider range of invading DNA elements than previously appreciated. Whether it is more advantageous for the peculiar neurotoxigenic C. butyricum type E strains to maintain or to lose the CRISPR-cas system remains an open question.