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

Genomic reconstruction of unclassified microorganisms: Analysis of CRISPR arrays and genes involved in defense mechanisms

The constant battle between bacteria and viruses has led to the development of sophisticated antiviral defense strategies by bacteria to defend themselves against phages. This study analyzed a marshland metagenome to identify and characterize bacterial antiviral defense systems and phage interactions. We assembled 210 metagenome-assembled genomes (MAGs) from environmental DNA extracted from Pallikaranai marshland soil and 37 unclassified MAGs were filtered. MIMAG standards were followed, 2 high-quality and 15 medium-quality unclassified MAGs were picked. MINCED was used to identify 137 CRISPR arrays in the quality MAGs, and ViroBLAST was used to identify the phages that interact with the bacteria. About 242 spacer sequences were extracted from the CRISPR arrays, of which 54 had significant matches in the ViroBLAST database. 7 unverified bacteriophage species were also detected in the MAGs. The viral group of Caudoviricetes phage elements were identified as a frequent genome terminal repeat. The PADLOC identified 11 genes involved as a defense system in the MAGs. The PD-T4-6 defense system was found to be prevalent in 15 different unclassified MAGs. This study presents valuable insights intothe adaptations of unclassified bacteria to bacteriophages, as well as the genes used by these bacteria as a defense mechanism.

The Cas10 nuclease activity relieves host dormancy to facilitate spacer acquisition and retention during type III-A CRISPR immunity

A hallmark of CRISPR immunity is the acquisition of short viral DNA sequences, known as spacers, that are transcribed into guide RNAs to recognize complementary sequences. The staphylococcal type III-A CRISPR-Cas system uses guide RNAs to locate viral transcripts and start a response that displays two mechanisms of immunity. When immunity is triggered by an early-expressed phage RNA, degradation of viral ssDNA can cure the host from infection. In contrast, when the RNA guide targets a late-expressed transcript, defense requires the activity of Csm6, a non-specific RNase. Here we show that Csm6 triggers a growth arrest of the host that provides immunity at the population level which hinders viral propagation to allow the replication of non-infected cells. We demonstrate that this mechanism leads to defense against not only the target phage but also other viruses present in the population that fail to replicate in the arrested cells. On the other hand, dormancy limits the acquisition and retention of spacers that trigger it. We found that the ssDNase activity of type III-A systems is required for the re-growth of a subset of the arrested cells, presumably through the degradation of the phage DNA, ending target transcription and inactivating the immune response. Altogether, our work reveals a built-in mechanism within type III-A CRISPR-Cas systems that allows the exit from dormancy needed for the subsistence of spacers that provide broad-spectrum immunity.

Isolation and genomic analysis of a type IV pili-independent Thermus thermophilus phage, φMN1 from a Japanese hot spring

A Thermus thermophilus lytic phage was isolated from a Japanese hot spring using a type IV pili-deficient strain as an indicator host, and designated as φMN1. Electron microscopic (EM) examination revealed that φMN1 had an icosahedral head and a contractile tail, suggesting that φMN1 belonged to Myoviridae. An EM analysis focused on φMN1 adsorption to the Thermus host cell showed that the receptor molecules for the phage were uniformly distributed on the outer surface of the cells. The circular double-stranded DNA of φMN1 was 76,659 base pairs in length, and the guanine and cytosine content was 61.8%. It was predicted to contain 99 open reading frames, and its putative distal tail fiber protein, which is essential for non-piliated host cell surface receptor recognition, was dissimilar in terms of sequence and length with its counterpart in the type IV pili-dependent φYS40. A phage proteomic tree revealed that φMN1 and φYS40 are in the same cluster, but many genes had low sequence similarities and some seemed to be derived from both mesophilic and thermophilic organisms. The gene organization suggested that φMN1 evolved from a non-Thermus phage through large-scale recombination events of the genes determining the host specificity, followed by gradual evolution by recombination of both the thermophilic and mesophilic DNAs assimilated by the host Thermus cells. This newly isolated phage will provide evolutionary insights into thermophilic phages.

Functional biology and biotechnology of thermophilic viruses

Viruses have developed sophisticated biochemical and genetic mechanisms to manipulate and exploit their hosts. Enzymes derived from viruses have been essential research tools since the first days of molecular biology. However, most viral enzymes that have been commercialized are derived from a small number of cultivated viruses, which is remarkable considering the extraordinary diversity and abundance of viruses revealed by metagenomic analysis. Given the explosion of new enzymatic reagents derived from thermophilic prokaryotes over the past 40 years, those obtained from thermophilic viruses should be equally potent tools. This review discusses the still-limited state of the art regarding the functional biology and biotechnology of thermophilic viruses with a focus on DNA polymerases, ligases, endolysins, and coat proteins. Functional analysis of DNA polymerases and primase-polymerases from phages infecting Thermus, Aquificaceae, and Nitratiruptor has revealed new clades of enzymes with strong proofreading and reverse transcriptase capabilities. Thermophilic RNA ligase 1 homologs have been characterized from Rhodothermus and Thermus phages, with both commercialized for circularization of single-stranded templates. Endolysins from phages infecting Thermus, Meiothermus, and Geobacillus have shown high stability and unusually broad lytic activity against Gram-negative and Gram-positive bacteria, making them targets for commercialization as antimicrobials. Coat proteins from thermophilic viruses infecting Sulfolobales and Thermus strains have been characterized, with diverse potential applications as molecular shuttles. To gauge the scale of untapped resources for these proteins, we also document over 20,000 genes encoded by uncultivated viral genomes from high-temperature environments that encode DNA polymerase, ligase, endolysin, or coat protein domains.

High-efficiency genome editing of an extreme thermophile Thermus thermophilus using endogenous type I and type III CRISPR-Cas systems

Thermus thermophilus is an attractive species in the bioindustry due to its valuable natural products, abundant thermophilic enzymes, and promising fermentation capacities. However, efficient and versatile genome editing tools are not available for this species. In this study, we developed an efficient genome editing tool for T. thermophilus HB27 based on its endogenous type I-B, I-C, and III-A/B CRISPR-Cas systems. First, we systematically characterized the DNA interference capabilities of the different types of the native CRISPR-Cas systems in T. thermophilus HB27. We found that genomic manipulations such as gene deletion, mutation, and in situ tagging could be easily implemented by a series of genome-editing plasmids carrying an artificial self-targeting mini-CRISPR and a donor DNA responsible for the recombinant recovery. We also compared the genome editing efficiency of different CRISPR-Cas systems and the editing plasmids with donor DNAs of different lengths. Additionally, we developed a reporter gene system for T. thermophilus based on a heat-stable β-galactosidase gene TTP0042, and constructed an engineered strain with a high production capacity of superoxide dismutases by genome modification.

Crystal structure of DNA polymerase I from Thermus phage G20c

This study describes the structure of DNA polymerase I from Thermus phage G20c, termed PolI_G20c. This is the first structure of a DNA polymerase originating from a group of related thermophilic bacteriophages infecting Thermus thermophilus, including phages G20c, TSP4, P74-26, P23-45 and phiFA and the novel phage Tth15-6. Sequence and structural analysis of PolI_G20c revealed a 3′–5′ exonuclease domain and a DNA polymerase domain, and activity screening confirmed that both domains were functional. No functional 5′–3′ exonuclease domain was present. Structural analysis also revealed a novel specific structure motif, here termed SβαR, that was not previously identified in any polymerase belonging to the DNA polymerases I (or the DNA polymerase A family). The SβαR motif did not show any homology to the sequences or structures of known DNA polymerases. The exception was the sequence conservation of the residues in this motif in putative DNA polymerases encoded in the genomes of a group of thermophilic phages related to Thermus phage G20c. The structure of PolI_G20c was determined with the aid of another structure that was determined in parallel and was used as a model for molecular replacement. This other structure was of a 3′–5′ exonuclease termed ExnV1. The cloned and expressed gene encoding ExnV1 was isolated from a thermophilic virus metagenome that was collected from several hot springs in Iceland. The structure of ExnV1, which contains the novel SβαR motif, was first determined to 2.19 Å resolution. With these data at hand, the structure of PolI_G20c was determined to 2.97 Å resolution. The structures of PolI_G20c and ExnV1 are most similar to those of the Klenow fragment of DNA polymerase I (PDB entry 2kzz) from Escherichia coli, DNA polymerase I from Geobacillus stearothermophilus (PDB entry 1knc) and Taq polymerase (PDB entry 1bgx) from Thermus aquaticus.

Varidnaviruses in the Human Gut: A Major Expansion of the Order Vinavirales

Bacteriophages play key roles in the dynamics of the human microbiome. By far the most abundant components of the human gut virome are tailed bacteriophages of the realm Duplodnaviria, in particular, crAss-like phages. However, apart from duplodnaviruses, the gut virome has not been dissected in detail. Here we report a comprehensive census of a minor component of the gut virome, the tailless bacteriophages of the realm Varidnaviria. Tailless phages are primarily represented in the gut by prophages, that are mostly integrated in genomes of Alphaproteobacteria and Verrucomicrobia and belong to the order Vinavirales, which currently consists of the families Corticoviridae and Autolykiviridae. Phylogenetic analysis of the major capsid proteins (MCP) suggests that at least three new families should be established within Vinavirales to accommodate the diversity of prophages from the human gut virome. Previously, only the MCP and packaging ATPase genes were reported as conserved core genes of Vinavirales. Here we report an extended core set of 12 proteins, including MCP, packaging ATPase, and previously undetected lysis enzymes, that are shared by most of these viruses. We further demonstrate that replication system components are frequently replaced in the genomes of Vinavirales, suggestive of selective pressure for escape from yet unknown host defenses or avoidance of incompatibility with coinfecting related viruses. The results of this analysis show that, in a sharp contrast to marine viromes, varidnaviruses are a minor component of the human gut virome. Moreover, they are primarily represented by prophages, as indicated by the analysis of the flanking genes, suggesting that there are few, if any, lytic varidnavirus infections in the gut at any given time. These findings complement the existing knowledge of the human gut virome by exploring a group of viruses that has been virtually overlooked in previous work.

Analysis of Probiotic Bacteria Genomes: Comparison of CRISPR/Cas Systems and Spacer Acquisition Diversity

Clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated (Cas) genes constitute an adaptive (acquired) defense system of bacteria and archaea. Here 72 probiotic bacteria genomes were investigated in terms of the presence of CRISPR/Cas systems and phage/plasmid invaders through spacer analysis. 49 CRISPR/Cas systems were detected within probiotic strains, namely,17 type II-A, 10 type I-C, 8 type I-E, 5 Type I-U (I-G), 4 type III-A, 2 type I-B, 1 type I-A, 1 type IV-B, and 1 type II-C. The predicted target of spacers was determined in 25 strains and consequently, three different spacer and target patterns were revealed. The diversity of CRISPR spacers provides insight and understanding to determine strain-specific invaders of probiotic bacteria as well as their relationships between strains. CRISPR systems were clarified in many studies for genomic characterization. However, recently, endogenous genome editing with CRISPR has provided an approach for various genome editing projects. Thus, in the future, producing strain-specific phage-resistant starter cultures or probiotics by endogenous genome editing methods according to phage/plasmid survey can be utilized for industrial and pharmaceutical applications. Therefore, this study intended a comprehensive investigation of CRISPR systems of probiotic bacteria.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12088-021-00971-1.

SCRAMBLER: A Tool for De Novo CRISPR Array Reconstruction and Its Application for Analysis of the Structure of Prokaryotic Populations

CRISPR arrays are prokaryotic genomic loci consisting of repeat sequences alternating with unique spacers acquired from foreign nucleic acids. As one of the fastest-evolving parts of the genome, CRISPR arrays can be used to differentiate closely related prokaryotic lineages and track individual strains in prokaryotic communities. However, the assembly of full-length CRISPR arrays sequences remains a problem. Here, we developed SCRAMBLER, a tool that includes several pipelines for assembling CRISPR arrays from high-throughput short-read sequencing data. We assessed its performance with model data sets (Escherichia coli strains containing different CRISPR arrays and imitating prokaryotic communities of different complexities) and intestinal microbiomes of extant and extinct pachyderms. Evaluation of SCRAMBLER's performance using model data sets demonstrated its ability to assemble CRISPR arrays correctly from reads containing pairs of spacers, yielding a precision rate of >80% and a recall rate of 60-85% when checked against ground-truth data. Likewise, SCRAMBLER successfully assembled CRISPR arrays from the environmental samples, as attested by their matching with database entries. SCRAMBLER, an open-source software (github.com/biolab-tools/SCRAMBLER), can facilitate analysis of the composition and dynamics of CRISPR arrays in complex communities.

Evolutionary stasis of a deep subsurface microbial lineage

Sulfate-reducing bacteria Candidatus Desulforudis audaxviator (CDA) were originally discovered in deep fracture fluids accessed via South African gold mines and have since been found in geographically widespread deep subsurface locations. In order to constrain models for subsurface microbial evolution, we compared CDA genomes from Africa, North America and Eurasia using single cell genomics. Unexpectedly, 126 partial single amplified genomes from the three continents, a complete genome from of an isolate from Eurasia, and metagenome-assembled genomes from Africa and Eurasia shared >99.2% average nucleotide identity, low frequency of SNP's, and near-perfectly conserved prophages and CRISPRs. Our analyses reject sample cross-contamination, recent natural dispersal, and unusually strong purifying selection as likely explanations for these unexpected results. We therefore conclude that the analyzed CDA populations underwent only minimal evolution since their physical separation, potentially as far back as the breakup of Pangea between 165 and 55 Ma ago. High-fidelity DNA replication and repair mechanisms are the most plausible explanation for the highly conserved genome of CDA. CDA presents a stark contrast to the current model organisms in microbial evolutionary studies, which often develop adaptive traits over far shorter periods of time.

PAM-repeat associations and spacer selection preferences in single and co-occurring CRISPR-Cas systems

BACKGROUND

The adaptive CRISPR-Cas immune system stores sequences from past invaders as spacers in CRISPR arrays and thereby provides direct evidence that links invaders to hosts. Mapping CRISPR spacers has revealed many aspects of CRISPR-Cas biology, including target requirements such as the protospacer adjacent motif (PAM). However, studies have so far been limited by a low number of mapped spacers in the database.

RESULTS

By using vast metagenomic sequence databases, we map approximately one-third of more than 200,000 unique CRISPR spacers from a variety of microbes and derive a catalog of more than two hundred unique PAM sequences associated with specific CRISPR-Cas subtypes. These PAMs are further used to correctly assign the orientation of CRISPR arrays, revealing conserved patterns between the last nucleotides of the CRISPR repeat and PAM. We could also deduce CRISPR-Cas subtype-specific preferences for targeting either template or coding strand of open reading frames. While some DNA-targeting systems (type I-E and type II systems) prefer the template strand and avoid mRNA, other DNA- and RNA-targeting systems (types I-A and I-B and type III systems) prefer the coding strand and mRNA. In addition, we find large-scale evidence that both CRISPR-Cas adaptation machinery and CRISPR arrays are shared between different CRISPR-Cas systems. This could lead to simultaneous DNA and RNA targeting of invaders, which may be effective at combating mobile genetic invaders.

CONCLUSIONS

This study has broad implications for our understanding of how CRISPR-Cas systems work in a wide range of organisms for which only the genome sequence is known.

Investigation into the prevalent CRISPR-Cas systems among the Aeromonas genus

Clustered regularly interspaced short palindromic repeats (CRISPR) is a prokaryotic adaptive immune system that checks invasion by mobile genetic elements through nuclease targeting. In this study, we investigated the occurrence, diversity, and features of the CRISPR system in the genus Aeromonas using bioinformatics tools. Only 13 out of 122 complete genomes (10.66%) of the genus Aeromonas from the NCBI GenBank database harbored the CRISPR system. The Type I-F system was the most prevalent CRISPR system among the Aeromonads, followed by the Type I-E system. Only one strain harbored a Type I-C CRISPR system. Among the Aeromonads, Aeromonas caviae (22.7%) and Aeromonas veronii (20%) had a higher prevalence rate of the complete CRISPR system. The analysis of direct repeat (DR) sequences showed that all could form stable RNA secondary structures. A phylogenetic tree generated for the Cas1 protein classified CRISPR subtypes into three distinct clusters. Among the 748 spacers investigated, 41.98% and 17.25% showed perfect homology to phage and plasmid sequences, respectively. Some arrays had duplicated spacers. The CRISPR loci are closely linked to antibiotic resistance genes in most strains. Collectively, our results would contribute to research on antibiotic resistance in the Aeromonas group, and provide new insights into the diversity and evolution of the CRISPR-Cas system.

Characterization of CRISPR Spacer and Protospacer Sequences in Paenibacillus larvae and Its Bacteriophages

The bacterium Paenibacillus larvae is the causative agent of American foulbrood, the most devastating bacterial disease of honeybees. Because P. larvae is antibiotic resistant, phages that infect it are currently used as alternative treatments. However, the acquisition by P. larvae of CRISPR spacer sequences from the phages could be an obstacle to treatment efforts. We searched nine complete genomes of P. larvae strains and identified 714 CRISPR spacer sequences, of which 384 are unique. Of the four epidemiologically important P. larvae strains, three of these have fewer than 20 spacers, while one strain has over 150 spacers. Of the 384 unique spacers, 18 are found as protospacers in the genomes of 49 currently sequenced P. larvae phages. One P. larvae strain does not have any protospacers found in phages, while another has eight. Protospacer distribution in the phages is uneven, with two phages having up to four protospacers, while a third of phages have none. Some phages lack protospacers found in closely related phages due to point mutations, indicating a possible escape mechanism. This study serve a point of reference for future studies on the CRISPR-Cas system in P. larvae as well as for comparative studies of other phage–host systems.

Spacer acquisition by Type III CRISPR-Cas system during bacteriophage infection of Thermus thermophilus

Type III CRISPR–Cas systems provide immunity to foreign DNA by targeting its transcripts. Target recognition activates RNases and DNases that may either destroy foreign DNA directly or elicit collateral damage inducing death of infected cells. While some Type III systems encode a reverse transcriptase to acquire spacers from foreign transcripts, most contain conventional spacer acquisition machinery found in DNA-targeting systems. We studied Type III spacer acquisition in phage-infected Thermus thermophilus, a bacterium that lacks either a standalone reverse transcriptase or its fusion to spacer integrase Cas1. Cells with spacers targeting a subset of phage transcripts survived the infection, indicating that Type III immunity does not operate through altruistic suicide. In the absence of selection spacers were acquired from both strands of phage DNA, indicating that no mechanism ensuring acquisition of RNA-targeting spacers exists. Spacers that protect the host from the phage demonstrate a very strong strand bias due to positive selection during infection. Phages that escaped Type III interference accumulated deletions of integral number of codons in an essential gene and much longer deletions in a non-essential gene. This and the fact that Type III immunity can be provided by plasmid-borne mini-arrays open ways for genomic manipulation of Thermus phages.

Thermus thermophilus Argonaute Functions in the Completion of DNA Replication

In many eukaryotes, Argonaute proteins, guided by short RNA sequences, defend cells against transposons and viruses. In the eubacterium Thermus thermophilus, the DNA-guided Argonaute TtAgo defends against transformation by DNA plasmids. Here, we report that TtAgo also participates in DNA replication. In vivo, TtAgo binds 15- to 18-nt DNA guides derived from the chromosomal region where replication terminates and associates with proteins known to act in DNA replication. When gyrase, the sole T. thermophilus type II topoisomerase, is inhibited, TtAgo allows the bacterium to finish replicating its circular genome. In contrast, loss of gyrase and TtAgo activity slows growth and produces long sausage-like filaments in which the individual bacteria are linked by DNA. Finally, wild-type T. thermophilus outcompetes an otherwise isogenic strain lacking TtAgo. We propose that the primary role of TtAgo is to help T. thermophilus disentangle the catenated circular chromosomes generated by DNA replication.