Evolution of the CRISPR-Cas adaptive immunity systems in prokaryotes: models and observations on virus-host coevolution

Efficient coupling of within-and between-host infectious disease dynamics

Mathematical models of infectious disease transmission typically neglect within-host dynamics. Yet within-host dynamics - including pathogen replication, host immune responses, and interactions with microbiota - are crucial not only for determining the progression of disease at the individual level, but also for driving within-host evolution and onwards transmission, and therefore shape dynamics at the population level. Various approaches have been proposed to model both within- and between-host dynamics, but these typically require considerable simplifying assumptions to couple processes at contrasting scales (e.g., the within-host dynamics quickly reach a steady state) or are computationally intensive. Here we propose a novel, readily adaptable and broadly applicable method for modelling both within- and between-host processes which can fully couple dynamics across scales and is both realistic and computationally efficient. By individually tracking the deterministic within-host dynamics of infected individuals, and stochastically coupling these to continuous host state variables at the population-level, we take advantage of fast numerical methods at both scales while still capturing individual transient within-host dynamics and stochasticity in transmission between hosts. Our approach closely agrees with full stochastic individual-based simulations and is especially useful when the within-host dynamics do not rapidly reach a steady state or over longer timescales to track pathogen evolution. By applying our method to different pathogen growth scenarios we show how common simplifying assumptions fundamentally change epidemiological and evolutionary dynamics.

Two multidrug-resistant Proteus mirabilis clones carrying extended spectrum beta-lactamases revealed in a single hospital department by whole genome sequencing

Proteus mirabilis bacteria is a component of normal intestinal microflora of humans and animals, but can also be found in hospital settings causing urinary tract infections and sepsis. The problem of treating such infections is complicated by multidrug-resistant isolates producing extended spectrum beta-lactamases (ESBL), and the number of ESBL-carrying P. mirabilis strains has significantly increased recently. This study presents a detailed analysis of 12 multidrug-resistant P. mirabilis isolates obtained from the wounds of different patients in one surgical department of a multidisciplinary hospital in Moscow, Russia, using the short- and long-read whole genome sequencing. The isolates under investigation divided into two clusters (clones) C1 and C2 based on their genomic profiles and carried antimicrobial resistance (AMR) genes corresponding well with phenotypic profiles, which was the first case of reporting two different P. mirabilis clones obtained simultaneously from the same specimens at one hospital, to the best of our knowledge. Some genes, including ESBL encoding ones, were specific for either C1 or C2 (aac(6')-Ib10, ant(2″)-Ia, qnrA1, bla VEB-6 and fosA3, bla CTX -M-65 , correspondingly). Additionally, the Salmonella genomic islands 1 were found that differed in composition of multiple antibiotic resistance regions between C1 and C2 groups. CRISPR-Cas system type I-E was revealed only in C2 isolates, while the same set of virulence factors was determined for both P. mirabilis clones. Diversity of all genetic factors found in case of simultaneous existence of two clones collected from the same source at one department indicates high pathogenic potential of P. mirabilis and poses a requirement of proper spreading monitoring. The data obtained will facilitate the understanding of AMR transfer and dynamics within clinical P. mirabilis isolates and contribute to epidemiological surveillance of this pathogen.

SpacerPlacer: ancestral reconstruction of CRISPR arrays reveals the evolutionary dynamics of spacer deletions

Bacteria employ CRISPR-Cas systems for defense by integrating invader-derived sequences, termed spacers, into the CRISPR array, which constitutes an immunity memory. While spacer deletions occur randomly across the array, newly acquired spacers are predominantly integrated at the leader end. Consequently, spacer arrays can be used to derive the chronology of spacer insertions. Reconstruction of ancestral spacer acquisitions and deletions could help unravel the coevolution of phages and bacteria, the evolutionary dynamics in microbiomes, or track pathogens. However, standard reconstruction methods produce misleading results by overlooking insertion order and joint deletions of spacers. Here, we present SpacerPlacer, a maximum likelihood-based ancestral reconstruction approach for CRISPR array evolution. We used SpacerPlacer to reconstruct and investigate ancestral deletion events of 4565 CRISPR arrays, revealing that spacer deletions occur 374 times more frequently than mutations and are regularly deleted jointly, with an average of 2.7 spacers. Surprisingly, we observed a decrease in the spacer deletion frequency towards both ends of the reconstructed arrays. While the resulting trailer-end conservation is commonly observed, a reduced deletion frequency is now also detectable towards the variable leader end. Finally, our results point to the hypothesis that frequent loss of recently acquired spacers may provide a selective advantage.

Combination of nucleic acid amplification and CRISPR/Cas technology in pathogen detection

Major health events caused by pathogenic microorganisms are increasing, seriously jeopardizing human lives. Currently PCR and ITA are widely used for rapid testing in food, medicine, industry and agriculture. However, due to the non-specificity of the amplification process, researchers have proposed the combination of nucleic acid amplification technology with the novel technology CRISPR for detection, which improves the specificity and credibility of results. This paper summarizes the research progress of nucleic acid amplification technology in conjunction with CRISPR/Cas technology for the detection of pathogens, which provides a reference and theoretical basis for the subsequent application of nucleic acid amplification technology in the field of pathogen detection.

Dynamics of CRISPR-mediated virus-host interactions in the human gut microbiome

Arms races between mobile genetic elements and prokaryotic hosts are major drivers of ecological and evolutionary change in microbial communities. Prokaryotic defense systems such as CRISPR-Cas have the potential to regulate microbiome composition by modifying the interactions among bacteria, plasmids, and phages. Here, we used longitudinal metagenomic data from 130 healthy and diseased individuals to study how the interplay of genetic parasites and CRISPR-Cas immunity reflects on the dynamics and composition of the human gut microbiome. Based on the coordinated study of 80 000 CRISPR-Cas loci and their targets, we show that CRISPR-Cas immunity effectively modulates bacteriophage abundances in the gut. Acquisition of CRISPR-Cas immunity typically leads to a decrease in the abundance of lytic phages but does not necessarily cause their complete disappearance. Much smaller effects are observed for lysogenic phages and plasmids. Conversely, phage-CRISPR interactions shape bacterial microdiversity by producing weak selective sweeps that benefit immune host lineages. We also show that distal (and chronologically older) regions of CRISPR arrays are enriched in spacers that are potentially functional and target crass-like phages and local prophages. This suggests that exposure to reactivated prophages and other endemic viruses is a major selective pressure in the gut microbiome that drives the maintenance of long-lasting immune memory.

Virus diversity and activity is driven by snowmelt and host dynamics in a high-altitude watershed soil ecosystem

BACKGROUND

Viruses impact nearly all organisms on Earth, including microbial communities and their associated biogeochemical processes. In soils, highly diverse viral communities have been identified, with a global distribution seemingly driven by multiple biotic and abiotic factors, especially soil temperature and moisture. However, our current understanding of the stability of soil viral communities across time and their response to strong seasonal changes in environmental parameters remains limited. Here, we investigated the diversity and activity of environmental soil DNA and RNA viruses, focusing especially on bacteriophages, across dynamics' seasonal changes in a snow-dominated mountainous watershed by examining paired metagenomes and metatranscriptomes.

RESULTS

We identified a large number of DNA and RNA viruses taxonomically divergent from existing environmental viruses, including a significant proportion of fungal RNA viruses, and a large and unsuspected diversity of positive single-stranded RNA phages (Leviviricetes), highlighting the under-characterization of the global soil virosphere. Among these, we were able to distinguish subsets of active DNA and RNA phages that changed across seasons, consistent with a "seed-bank" viral community structure in which new phage activity, for example, replication and host lysis, is sequentially triggered by changes in environmental conditions. At the population level, we further identified virus-host dynamics matching two existing ecological models: "Kill-The-Winner" which proposes that lytic phages are actively infecting abundant bacteria, and "Piggyback-The-Persistent" which argues that when the host is growing slowly, it is more beneficial to remain in a dormant state. The former was associated with summer months of high and rapid microbial activity, and the latter with winter months of limited and slow host growth.

CONCLUSION

Taken together, these results suggest that the high diversity of viruses in soils is likely associated with a broad range of host interaction types each adapted to specific host ecological strategies and environmental conditions. As our understanding of how environmental and host factors drive viral activity in soil ecosystems progresses, integrating these viral impacts in complex natural microbiome models will be key to accurately predict ecosystem biogeochemistry. Video Abstract.

Life strategies for Aminicenantia in subseafloor oceanic crust

After decades studying the microbial "deep biosphere" in subseafloor oceanic crust, the growth and life strategies in this anoxic, low energy habitat remain poorly described. Using both single cell genomics and metagenomics, we reveal the life strategies of two distinct lineages of uncultivated Aminicenantia bacteria from the basaltic subseafloor oceanic crust of the eastern flank of the Juan de Fuca Ridge. Both lineages appear adapted to scavenge organic carbon, as each have genetic potential to catabolize amino acids and fatty acids, aligning with previous Aminicenantia reports. Given the organic carbon limitation in this habitat, seawater recharge and necromass may be important carbon sources for heterotrophic microorganisms inhabiting the ocean crust. Both lineages generate ATP via several mechanisms including substrate-level phosphorylation, anaerobic respiration, and electron bifurcation driving an Rnf ion translocation membrane complex. Genomic comparisons suggest these Aminicenantia transfer electrons extracellularly, perhaps to iron or sulfur oxides consistent with mineralogy of this site. One lineage, called JdFR-78, has small genomes that are basal to the Aminicenantia class and potentially use "primordial" siroheme biosynthetic intermediates for heme synthesis, suggesting this lineage retain characteristics of early evolved life. Lineage JdFR-78 contains CRISPR-Cas defenses to evade viruses, while other lineages contain prophage that may help prevent super-infection or no detectable viral defenses. Overall, genomic evidence points to Aminicenantia being well adapted to oceanic crust environments by taking advantage of simple organic molecules and extracellular electron transport.

Diversity and dynamics of the CRISPR-Cas systems associated with Bacteroides fragilis in human population

Background

CRISPR-Cas (clustered regularly interspaced short palindromic repeats—CRISPR-associated proteins) systems are adaptive immune systems commonly found in prokaryotes that provide sequence-specific defense against invading mobile genetic elements (MGEs). The memory of these immunological encounters are stored in CRISPR arrays, where spacer sequences record the identity and history of past invaders. Analyzing such CRISPR arrays provide insights into the dynamics of CRISPR-Cas systems and the adaptation of their host bacteria to rapidly changing environments such as the human gut.

Results

In this study, we utilized 601 publicly available Bacteroides fragilis genome isolates from 12 healthy individuals, 6 of which include longitudinal observations, and 222 available B. fragilis reference genomes to update the understanding of B. fragilis CRISPR-Cas dynamics and their differential activities. Analysis of longitudinal genomic data showed that some CRISPR array structures remained relatively stable over time whereas others involved radical spacer acquisition during some periods, and diverse CRISPR arrays (associated with multiple isolates) co-existed in the same individuals with some persisted over time. Furthermore, features of CRISPR adaptation, evolution, and microdynamics were highlighted through an analysis of host-MGE network, such as modules of multiple MGEs and hosts, reflecting complex interactions between B. fragilis and its invaders mediated through the CRISPR-Cas systems.

Conclusions

We made available of all annotated CRISPR-Cas systems and their target MGEs, and their interaction network as a web resource at https://omics.informatics.indiana.edu/CRISPRone/Bfragilis. We anticipate it will become an important resource for studying of B. fragilis, its CRISPR-Cas systems, and its interaction with mobile genetic elements providing insights into evolutionary dynamics that may shape the species virulence and lead to its pathogenicity.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-022-08770-8).

Predation strategies of the bacterium Bdellovibrio bacteriovorus result in overexploitation and bottlenecks

With increasing antimicrobial resistance, alternatives for treating infections or removing resistant bacteria are urgently needed, such as the bacterial predator Bdellovibrio bacteriovorus or bacteriophage. Therefore, we need to better understand microbial predator-prey dynamics. We developed mass-action mathematical models of predation for chemostats, which capture the low substrate concentration and slow growth typical for intended application areas of the predators such as wastewater treatment, aquaculture or the gut. Our model predicted that predator survival required a minimal prey cell size, explaining why Bdellovibrio is much smaller than its prey. A too good predator (attack rate too high, mortality too low) overexploited its prey leading to extinction (tragedy of the commons). Surprisingly, a predator taking longer to produce more offspring outcompeted a predator producing fewer offspring more rapidly (rate versus yield trade-off). Predation was only efficient in a narrow region around optimal parameters. Moreover, extreme oscillations under a wide range of conditions led to severe bottlenecks. These could be avoided when two prey species became available in alternating seasons. A bacteriophage outcompeted Bdellovibrio due to its higher burst size and faster life cycle. Together, results suggest that Bdellovibrio would struggle to survive on a single prey, explaining why it must be a generalist predator and suggesting it is better suited than phage to environments with multiple prey. Importance The discovery of antibiotics led to a dramatic drop in deaths due to infectious disease. Increasing levels of antimicrobial resistance, however, threaten to reverse this progress. There is thus a need for alternatives, such as therapies based on phage and predatory bacteria that kill bacteria regardless of whether they are pathogens or resistant to antibiotics. To best exploit them, we need to better understand what determines their effectiveness. By using a mathematical model to study bacterial predation in realistic slow growth conditions, we found that the generalist predator Bdellovibrio is most effective within a narrow range of conditions for each prey. For example, a minimum prey cell size is required, and the predator should not be too good as this would result in over-exploitation risking extinction. Together these findings give insights into the ecology of microbial predation and help explain why Bdellovibrio needs to be a generalist predator.

Searching for fat tails in CRISPR-Cas systems: Data analysis and mathematical modeling

Understanding CRISPR-Cas systems-the adaptive defence mechanism that about half of bacterial species and most of archaea use to neutralise viral attacks-is important for explaining the biodiversity observed in the microbial world as well as for editing animal and plant genomes effectively. The CRISPR-Cas system learns from previous viral infections and integrates small pieces from phage genomes called spacers into the microbial genome. The resulting library of spacers collected in CRISPR arrays is then compared with the DNA of potential invaders. One of the most intriguing and least well understood questions about CRISPR-Cas systems is the distribution of spacers across the microbial population. Here, using empirical data, we show that the global distribution of spacer numbers in CRISPR arrays across multiple biomes worldwide typically exhibits scale-invariant power law behaviour, and the standard deviation is greater than the sample mean. We develop a mathematical model of spacer loss and acquisition dynamics which fits observed data from almost four thousand metagenomes well. In analogy to the classical 'rich-get-richer' mechanism of power law emergence, the rate of spacer acquisition is proportional to the CRISPR array size, which allows a small proportion of CRISPRs within the population to possess a significant number of spacers. Our study provides an alternative explanation for the rarity of all-resistant super microbes in nature and why proliferation of phages can be highly successful despite the effectiveness of CRISPR-Cas systems.

Kinetic Fingerprinting Links Bacteria-Phage Interactions with Emergent Dynamics: Rapid Depletion of Klebsiella pneumoniae Indicates Phage Synergy

The specific temporal evolution of bacterial and phage population sizes, in particular bacterial depletion and the emergence of a resistant bacterial population, can be seen as a kinetic fingerprint that depends on the manifold interactions of the specific phage-host pair during the course of infection. We have elaborated such a kinetic fingerprint for a human urinary tract Klebsiella pneumoniae isolate and its phage vB_KpnP_Lessing by a modeling approach based on data from in vitro co-culture. We found a faster depletion of the initially sensitive bacterial population than expected from simple mass action kinetics. A possible explanation for the rapid decline of the bacterial population is a synergistic interaction of phages which can be a favorable feature for phage therapies. In addition to this interaction characteristic, analysis of the kinetic fingerprint of this bacteria and phage combination revealed several relevant aspects of their population dynamics: A reduction of the bacterial concentration can be achieved only at high multiplicity of infection whereas bacterial extinction is hardly accomplished. Furthermore the binding affinity of the phage to bacteria is identified as one of the most crucial parameters for the reduction of the bacterial population size. Thus, kinetic fingerprinting can be used to infer phage-host interactions and to explore emergent dynamics which facilitates a rational design of phage therapies.

The size of the immune repertoire of bacteria

Some bacteria possess an adaptive immune system that maintains a memory of past viral infections in the CRISPR loci of their genomes. This memory is used to mount targeted responses against later threats but is remarkably shallow: it remembers only a few dozen to a few hundred viruses. We present a statistical theory of CRISPR-based immunity that quantitatively predicts the depth of bacterial immune memory in terms of a tradeoff with fundamental constraints of the cellular biochemical machinery.

Some bacteria and archaea possess an immune system, based on the CRISPR-Cas mechanism, that confers adaptive immunity against viruses. In such species, individual prokaryotes maintain cassettes of viral DNA elements called spacers as a memory of past infections. Typically, the cassettes contain several dozen expressed spacers. Given that bacteria can have very large genomes and since having more spacers should confer a better memory, it is puzzling that so little genetic space would be devoted by prokaryotes to their adaptive immune systems. Here, assuming that CRISPR functions as a long-term memory-based defense against a diverse landscape of viral species, we identify a fundamental tradeoff between the amount of immune memory and effectiveness of response to a given threat. This tradeoff implies an optimal size for the prokaryotic immune repertoire in the observational range.

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.

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.

Genome-wide correlation analysis suggests different roles of CRISPR-Cas systems in the acquisition of antibiotic resistance genes in diverse species

CRISPR-Cas systems are widespread in bacterial and archaeal genomes, and in their canonical role in phage defence they confer a fitness advantage. However, CRISPR-Cas may also hinder the uptake of potentially beneficial genes. This is particularly true under antibiotic selection, where preventing the uptake of antibiotic resistance genes could be detrimental. Newly discovered features within these evolutionary dynamics are anti-CRISPR genes, which inhibit specific CRISPR-Cas systems. We hypothesized that selection for antibiotic resistance might have resulted in an accumulation of anti-CRISPR genes in genomes that harbour CRISPR-Cas systems and horizontally acquired antibiotic resistance genes. To assess that question, we analysed correlations between the CRISPR-Cas, anti-CRISPR and antibiotic resistance gene content of 104 947 reference genomes, including 5677 different species. In most species, the presence of CRISPR-Cas systems did not correlate with the presence of antibiotic resistance genes. However, in some clinically important species, we observed either a positive or negative correlation of CRISPR-Cas with antibiotic resistance genes. Anti-CRISPR genes were common enough in four species to be analysed. In Pseudomonas aeruginosa, the presence of anti-CRISPRs was associated with antibiotic resistance genes. This analysis indicates that the role of CRISPR-Cas and anti-CRISPRs in the spread of antibiotic resistance is likely to be very different in particular pathogenic species and clinical environments. This article is part of a discussion meeting issue 'The ecology and evolution of prokaryotic CRISPR-Cas adaptive immune systems'.

Programmable Molecular Scissors: Applications of a New Tool for Genome Editing in Biotech

Targeted genome editing is an advanced technique that enables precise modification of the nucleic acid sequences in a genome. Genome editing is typically performed using tools, such as molecular scissors, to cut a defined location in a specific gene. Genome editing has impacted various fields of biotechnology, such as agriculture; biopharmaceutical production; studies on the structure, regulation, and function of the genome; and the creation of transgenic organisms and cell lines. Although genome editing is used frequently, it has several limitations. Here, we provide an overview of well-studied genome-editing nucleases, including single-stranded oligodeoxynucleotides (ssODNs), transcription activator-like effector nucleases (TALENs), zinc-finger nucleases (ZFNs), and CRISPR-Cas9 RNA-guided nucleases (CRISPR-Cas9). To this end, we describe the progress toward editable nuclease-based therapies and discuss the minimization of off-target mutagenesis. Future prospects of this challenging scientific field are also discussed.

Environmental structure drives resistance to phages and antibiotics during phage therapy and to invading lysogens during colonisation

Microbial communities are shaped by bacteriophages through predation and lysogeny. A better understanding of the interactions between these processes across different types of environments is key to elucidate how phages mediate microbial competition and to design efficient phage therapies. We introduce an individual-based model (eVIVALDI) to investigate the role of environmental structure in the elimination of a population with a combined treatment of antibiotics and virulent phages, and in the invasion of a population of phage-sensitive bacteria by lysogens. We show that structured environments facilitate the emergence of double resistance, to antibiotics and phages, due to limited diffusion of phage particles and increased nutrient availability from dead cells. They also hinder phage amplification, thus decreasing the generation of phage genetic diversity and increasing the unpredictability of phage-bacteria arms-races. We used a machine learning approach to determine the variables most important for the invasion of sensitive populations by lysogens. They revealed that phage-associated traits and environmental structure are the key drivers of the process. Structured environments hinder invasions, and accounting for their existence improves the fit of the model to published in vivo experimental data. Our results underline environmental structure as key to understand in vivo phage-bacteria interactions.

CRISPRs for Strain Tracking and Their Application to Microbiota Transplantation Data Analysis

CRISPR-Cas systems are adaptive immune systems naturally found in bacteria and archaea. Prokaryotes use these immune systems to defend against invaders, which include phages, plasmids, and other mobile genetic elements. Relying on the integration of spacers derived from invader sequences (protospacers) into CRISPR loci (forming spacers flanked by repeats), CRISPR-Cas systems are able to store the memory of past immunological encounters. While CRISPR-Cas systems have evolved in response to invading mobile genetic elements, invaders have also developed mechanisms to avoid detection. As a result of an arms race between CRISPR-Cas systems and their targets, CRISPR arrays typically undergo rapid turnover of spacers through the acquisition and loss events. Additionally, microbiomes of different individuals rarely share spacers. Here, we present a computational pipeline, CRISPRtrack, for strain tracking based on CRISPR spacer content, and we applied it to fecal transplantation microbiome data to study the retention of donor strains in recipients. Our results demonstrate the potential use of CRISPRs as a simple yet effective tool for donor-strain tracking in fecal transplantation and as a general purpose tool for quantifying microbiome similarity.

CRISPR: a new principle of genome engineering linked to conceptual shifts in evolutionary biology

The CRISPR-Cas systems of bacterial and archaeal adaptive immunity have become a household name among biologists and even the general public thanks to the unprecedented success of the new generation of genome editing tools utilizing Cas proteins. However, the fundamental biological features of CRISPR-Cas are of no lesser interest and have major impacts on our understanding of the evolution of antivirus defense, host-parasite coevolution, self versus non-self discrimination and mechanisms of adaptation. CRISPR-Cas systems present the best known case in point for Lamarckian evolution, i.e. generation of heritable, adaptive genomic changes in response to encounters with external factors, in this case, foreign nucleic acids. CRISPR-Cas systems employ multiple mechanisms of self versus non-self discrimination but, as is the case with immune systems in general, are nevertheless costly because autoimmunity cannot be eliminated completely. In addition to the autoimmunity, the fitness cost of CRISPR-Cas systems appears to be determined by their inhibitory effect on horizontal gene transfer, curtailing evolutionary innovation. Hence the dynamic evolution of CRISPR-Cas loci that are frequently lost and (re)acquired by archaea and bacteria. Another fundamental biological feature of CRISPR-Cas is its intimate connection with programmed cell death and dormancy induction in microbes. In this and, possibly, other immune systems, active immune response appears to be coupled to a different form of defense, namely, "altruistic" shutdown of cellular functions resulting in protection of neighboring cells. Finally, analysis of the evolutionary connections of Cas proteins reveals multiple contributions of mobile genetic elements (MGE) to the origin of various components of CRISPR-Cas systems, furthermore, different biological systems that function by genome manipulation appear to have evolved convergently from unrelated MGE. The shared features of adaptive defense systems and MGE, namely the ability to recognize and cleave unique sites in genomes, make them ideal candidates for genome editing and engineering tools.

Non-conventional therapeutic technique to replace CRISPR bacteria from biofilm by inducible lysogen

Bacteriophage can be an effective means of regulating bacterial populations when conditions allow phage invasion of bacterial colonies. Phage can either infect and lyse a host cell, or insert their DNA into the host cell genome; the latter process is called lysogeny. The clustered regularly interspaced short palindromic repeat (CRISPR) system, linked with CRISPR-associated (Cas) genes, is a regulatory system present in a variety of bacteria which confers immunity against bacteriophage. Studies of the group behaviour of bacteria with CRISPR/Cas systems have provided evidence that CRISPR in lysogenized bacteria can cause an inability to form biofilm. This allows CRISPR-immune bacteria in biofilms to effectively resist phage therapy. Our recent work has described a potential therapeutic technique to eradicate CRISPR-immune bacteria from a biofilm by a continuous influx of lysogens carrying an identical phage sequence. However, this model predicted that the CRISPR-immune population could persist for long times before eradication. Our current focus is on the use of diverse lysogens against CRISPR-capable bacterial populations. The goal of this work is to find a suitable strategy which can eradicate bacteria with a CRISPR system through the influx of finite amounts of distinct lysogens over fixed intervals.

Study the Features of 57 Confirmed CRISPR Loci in 38 Strains of Staphylococcus aureus

Staphylococcus aureus is a foodborne pathogen that causes food contamination and food poisoning, which poses great harm to health, agriculture and other hosts. Clustered regularly interspaced short palindromic repeats (CRISPR) are a recently discovered bacterial immune system that resists foreign genes such as phage DNA. This system inhibits the transfer of specific movable genetic elements that match the CRISPR spacer sequences, thereby preventing the spread of drug-resistant genes between pathogens. In this study, 57 CRISPR loci were screened from 38 strains of S. aureus based on the CRISPR database, and bioinformatics tools were used to investigate the structural features and potential functions of S. aureus CRISPR loci. The results showed that most strains contained only one CRISPR locus, a few strains contained multiple loci with sparsely distributed sites. These loci mainly included highly conserved direct repeat sequences and highly variable spacer sequences, as well as polymorphic cas genes. In addition, the analysis of secondary structure of direct repeat RNA showed that all sites can form stable RNA secondary structure. The results of constructing phylogenetic tree based on spacer sequence showed that some strains contained a high degree of phylogenetic relationship, while the differences among other strains in evolutionary processes were quite obvious. Of the 57 CRISPR loci identified, only the cas gene was found near the 4 CRISPR loci.

The physicist's guide to one of biotechnology's hottest new topics: CRISPR-Cas

Clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated proteins (Cas) constitute a multi-functional, constantly evolving immune system in bacteria and archaea cells. A heritable, molecular memory is generated of phage, plasmids, or other mobile genetic elements that attempt to attack the cell. This memory is used to recognize and interfere with subsequent invasions from the same genetic elements. This versatile prokaryotic tool has also been used to advance applications in biotechnology. Here we review a large body of CRISPR-Cas research to explore themes of evolution and selection, population dynamics, horizontal gene transfer, specific and cross-reactive interactions, cost and regulation, non-immunological CRISPR functions that boost host cell robustness, as well as applicable mechanisms for efficient and specific genetic engineering. We offer future directions that can be addressed by the physics community. Physical understanding of the CRISPR-Cas system will advance uses in biotechnology, such as developing cell lines and animal models, cell labeling and information storage, combatting antibiotic resistance, and human therapeutics.

A Microfluidics and Agent-Based Modeling Framework for Investigating Spatial Organization in Bacterial Colonies: The Case of Pseudomonas Aeruginosa and H1-Type VI Secretion Interactions

The factors leading to changes in the organization of microbial assemblages at fine spatial scales are not well characterized or understood. However, they are expected to guide the succession of community development and function toward specific outcomes that could impact human health and the environment. In this study, we put forward a combined experimental and agent-based modeling framework and use it to interpret unique spatial organization patterns of H1-Type VI secretion system (T6SS) mutants of P. aeruginosa under spatial confinement. We find that key parameters, such as T6SS-mediated cell contact and lysis, spatial localization, relative species abundance, cell density and local concentrations of growth substrates and metabolites are influenced by spatial confinement. The model, written in the accessible programming language NetLogo, can be adapted to a variety of biological systems of interest and used to simulate experiments across a broad parameter space. It was implemented and run in a high-throughput mode by deploying it across multiple CPUs, with each simulation representing an individual well within a high-throughput microwell array experimental platform. The microfluidics and agent-based modeling framework we present in this paper provides an effective means by which to connect experimental studies in microbiology to model development. The work demonstrates progress in coupling experimental results to simulation while also highlighting potential sources of discrepancies between real-world experiments and idealized models.

Complete characterization of new isolates of Neptunomonas phycophila leads to emend its description and opens possibilities of biotechnological applications

Five strains were isolated from gonad of Great scallop (Pecten maximus) broodstock in a Norwegian hatchery. The study of 16S rRNA gene sequences showed that these isolates belong to Neptunomonas phycophila, a bacterium originally isolated from a symbiont of the anemone Aiptasia tagetes from Puerto Rico. The gyrB and rpoB genes sequences confirmed the affiliation of the scallop isolates to this species. Phenotypic characterization was performed and some differences between the Norwegian isolates and the type strain of N. phycophila were detected, such as ranges of temperature, pH, and tolerance to salinity or the use of several substrates as sole carbon source which lead to an emended description of the species. The strain 3CM2.5 showed phosphatidylethanolamine and phosphatidylglycerol as the major polar lipids. The whole genomes of the scallop strain 3CM2.5 and type strain of the species CECT 8716^T were obtained and the annotation of these genomes revealed the presence of genes involved in degradation of aromatic compounds in both strains. Results obtained not only widen the geographical and host ranges of N. phycophila, but also point out possible biotechnological applications for this bacterial species.

Structural Variation of Type I-F CRISPR RNA Guided DNA Surveillance

CRISPR-Cas systems are prokaryotic immune systems against invading nucleic acids. Type I CRISPR-Cas systems employ highly diverse, multi-subunit surveillance Cascade complexes that facilitate duplex formation between crRNA and complementary target DNA for R-loop formation, retention, and DNA degradation by the subsequently recruited nuclease Cas3. Typically, the large subunit recognizes bona fide targets through the PAM (protospacer adjacent motif), and the small subunit guides the non-target DNA strand. Here, we present the Apo- and target-DNA-bound structures of the I-Fv (type I-F variant) Cascade lacking the small and large subunits. Large and small subunits are functionally replaced by the 5' terminal crRNA cap Cas5fv and the backbone protein Cas7fv, respectively. Cas5fv facilitates PAM recognition from the DNA major groove site, in contrast to all other described type I systems. Comparison of the type I-Fv Cascade with an anti-CRISPR protein-bound I-F Cascade reveals that the type I-Fv structure differs substantially at known anti-CRISPR protein target sites and might therefore be resistant to viral Cascade interception.

Mathematical modelling of CRISPR-Cas system effects on biofilm formation

Clustered regularly interspaced short palindromic repeats (CRISPR), linked with CRISPR associated (Cas) genes, can confer adaptive immunity to bacteria, against bacteriophage infections. Thus from a therapeutic standpoint, CRISPR immunity increases biofilm resistance to phage therapy. Recently, however, CRISPR-Cas genes have been implicated in reducing biofilm formation in lysogenized cells. Thus CRISPR immunity can have complex effects on phage-host-lysogen interactions, particularly in a biofilm. In this contribution, we develop and analyse a series of dynamical systems to elucidate and disentangle these interactions. Two competition models are used to study the effects of lysogens (first model) and CRISPR-immune bacteria (second model) in the biofilm. In the third model, the effect of delivering lysogens to a CRISPR-immune biofilm is investigated. Using standard analyses of equilibria, stability and bifurcations, our models predict that lysogens may be able to displace CRISPR-immune bacteria in a biofilm, and thus suggest strategies to eliminate phage-resistant biofilms.

Evolutionary Genomics of Defense Systems in Archaea and Bacteria

Evolution of bacteria and archaea involves an incessant arms race against an enormous diversity of genetic parasites. Accordingly, a substantial fraction of the genes in most bacteria and archaea are dedicated to antiparasite defense. The functions of these defense systems follow several distinct strategies, including innate immunity; adaptive immunity; and dormancy induction, or programmed cell death. Recent comparative genomic studies taking advantage of the expanding database of microbial genomes and metagenomes, combined with direct experiments, resulted in the discovery of several previously unknown defense systems, including innate immunity centered on Argonaute proteins, bacteriophage exclusion, and new types of CRISPR-Cas systems of adaptive immunity. Some general principles of function and evolution of defense systems are starting to crystallize, in particular, extensive gain and loss of defense genes during the evolution of prokaryotes; formation of genomic defense islands; evolutionary connections between mobile genetic elements and defense, whereby genes of mobile elements are repeatedly recruited for defense functions; the partially selfish and addictive behavior of the defense systems; and coupling between immunity and dormancy induction/programmed cell death.

Embracing the enemy: the diversification of microbial gene repertoires by phage-mediated horizontal gene transfer

Bacteriophages and archaeal viruses contribute, through lysogenic conversion or transduction, to the horizontal transfer of genetic material between microbial genomes. Recent genomics, metagenomics, and single cell studies have shown that lysogenic conversion is widespread and provides hosts with adaptive traits often associated with biotic interactions. The quantification of the evolutionary impact of transduction has lagged behind and requires further theoretical and experimental work. Nevertheless, recent studies suggested that generalized transduction plays a role in the transfer of antibiotic resistance genes and in the acquisition of novel genes during intra-specific bacterial competition. The characteristics of transduction and lysogenic conversion complement those of other mechanisms of transfer, and could play a key role in the spread of adaptive genes between communities.

Reconstruction of the evolution of microbial defense systems

Background

Evolution of bacterial and archaeal genomes is a highly dynamic process that involves intensive loss of genes as well as gene gain via horizontal transfer, with a lesser contribution from gene duplication. The rates of these processes can be estimated by comparing genomes that are linked by an evolutionary tree. These estimated rates of genome dynamics events substantially differ for different functional classes of genes. The genes involved in defense against viruses and other invading DNA are among those that are gained and lost at the highest rates.

Results

We employed a stochastic birth-and-death model to obtain maximum likelihood estimates of the rates of gain and loss of defense genes in 35 groups of closely related bacterial genomes and one group of archaeal genomes. We find that on average, the defense genes experience 1.4 fold higher flux than the rest of microbial genes. This excessive flux of defense genes over the genomic mean is consistent across diverse microbial groups. The few exceptions include intracellular parasites with small, degraded genomes that possess few defense systems which are more stable than in other microbes. Generally, defense genes follow the previously established pattern of genome dynamics, with gene family loss being about 3 times more common than gain and an order of magnitude more common than expansion or contraction of gene families. Case by case analysis of the evolutionary dynamics of defense genes indicates frequent multiple events in the same locus and widespread involvement of mobile elements in the gain and loss of defense genes.

Conclusions

Evolution of microbial defense systems is highly dynamic but, notwithstanding the host-parasite arms race, generally follows the same trends that have been established for the rest of the genes. Apart from the paucity and the low flux of defense genes in parasitic bacteria with deteriorating genomes, there is no clear connection between the evolutionary regime of defense systems and microbial life style.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0942-y) contains supplementary material, which is available to authorized users.

A computational strategy for the search of regulatory small RNAs in Actinobacillus pleuropneumoniae

Bacterial regulatory small RNAs (sRNAs) play important roles in gene regulation and are frequently connected to the expression of virulence factors in diverse bacteria. Only a few sRNAs have been described for Pasteurellaceae pathogens and no in-depth analysis of sRNAs has been described for Actinobacillus pleuropneumoniae, the causative agent of porcine pleuropneumonia, responsible for considerable losses in the swine industry. To search for sRNAs in A. pleuropneumoniae, we developed a strategy for the computational analysis of the bacterial genome by using four algorithms with different approaches, followed by experimental validation. The coding strand and expression of 17 out of 23 RNA candidates were confirmed by Northern blotting, RT-PCR, and RNA sequencing. Among them, two are likely riboswitches, three are housekeeping regulatory RNAs, two are the widely studied GcvB and 6S sRNAs, and 10 are putative novel trans-acting sRNAs, never before described for any bacteria. The latter group has several potential mRNA targets, many of which are involved with virulence, stress resistance, or metabolism, and connect the sRNAs in a complex gene regulatory network. The sRNAs identified are well conserved among the Pasteurellaceae that are evolutionarily closer to A. pleuropneumoniae and/or share the same host. Our results show that the combination of newly developed computational programs can be successfully utilized for the discovery of novel sRNAs and indicate an intricate system of gene regulation through sRNAs in A. pleuropneumoniae and in other Pasteurellaceae, thus providing clues for novel aspects of virulence that will be explored in further studies.

Diverse evolutionary roots and mechanistic variations of the CRISPR-Cas systems

Adaptive immunity had been long thought of as an exclusive feature of animals. However, the discovery of the CRISPR-Cas defense system, present in almost half of prokaryotic genomes, proves otherwise. Because of the everlasting parasite-host arms race, CRISPR-Cas has rapidly evolved through horizontal transfer of complete loci or individual modules, resulting in extreme structural and functional diversity. CRISPR-Cas systems are divided into two distinct classes that each consist of three types and multiple subtypes. We discuss recent advances in CRISPR-Cas research that reveal elaborate molecular mechanisms and provide for a plausible scenario of CRISPR-Cas evolution. We also briefly describe the latest developments of a wide range of CRISPR-based applications.

Advancing microbial sciences by individual-based modelling

Remarkable technological advances have revealed ever more properties and behaviours of individual microorganisms, but the novel data generated by these techniques have not yet been fully exploited. In this Opinion article, we explain how individual-based models (IBMs) can be constructed based on the findings of such techniques and how they help to explore competitive and cooperative microbial interactions. Furthermore, we describe how IBMs have provided insights into self-organized spatial patterns from biofilms to the oceans of the world, phage-CRISPR dynamics and other emergent phenomena. Finally, we discuss how combining individual-based observations with IBMs can advance our understanding at both the individual and population levels, leading to the new approach of microbial individual-based ecology (μIBE).

The Contribution of Genetic Recombination to CRISPR Array Evolution

CRISPR (clustered regularly interspaced short palindromic repeats) is a microbial immune system against foreign DNA. Recognition sequences (spacers) encoded within the CRISPR array mediate the immune reaction in a sequence-specific manner. The known mechanisms for the evolution of CRISPR arrays include spacer acquisition from foreign DNA elements at the time of invasion and array erosion through spacer deletion. Here, we consider the contribution of genetic recombination between homologous CRISPR arrays to the evolution of spacer repertoire. Acquisition of spacers from exogenic arrays via recombination may confer the recipient with immunity against unencountered antagonists. For this purpose, we develop a novel method for the detection of recombination in CRISPR arrays by modeling the spacer order in arrays from multiple strains from the same species. Because the evolutionary signal of spacer recombination may be similar to that of pervasive spacer deletions or independent spacer acquisition, our method entails a robustness analysis of the recombination inference by a statistical comparison to resampled and perturbed data sets. We analyze CRISPR data sets from four bacterial species: two Gammaproteobacteria species harboring CRISPR type I and two Streptococcus species harboring CRISPR type II loci. We find that CRISPR array evolution in Escherichia coli and Streptococcus agalactiae can be explained solely by vertical inheritance and differential spacer deletion. In Pseudomonas aeruginosa, we find an excess of single spacers potentially incorporated into the CRISPR locus during independent acquisition events. In Streptococcus thermophilus, evidence for spacer acquisition by recombination is present in 5 out of 70 strains. Genetic recombination has been proposed to accelerate adaptation by combining beneficial mutations that arose in independent lineages. However, for most species under study, we find that CRISPR evolution is shaped mainly by spacer acquisition and loss rather than recombination. Since the evolution of spacer content is characterized by a rapid turnover, it is likely that recombination is not beneficial for improving phage resistance in the strains under study, or that it cannot be detected in the resolution of intraspecies comparisons.