Long-run bacteria-phage coexistence dynamics under natural habitat conditions in an environmental biotechnology system

Bacterial viruses are widespread and abundant across natural and engineered habitats. They influence ecosystem functioning through interactions with their hosts. Laboratory studies of phage–host pairs have advanced our understanding of phenotypic and genetic diversification in bacteria and phages. However, the dynamics of phage–host interactions have been seldom recorded in complex natural environments. We conducted an observational metagenomic study of the dynamics of interaction between Gordonia and their phages using a three-year data series of samples collected from a full-scale wastewater treatment plant. The aim was to obtain a comprehensive picture of the coevolution dynamics in naturally evolving populations at relatively high time resolution. Coevolution was followed by monitoring changes over time in the CRISPR loci of Gordonia metagenome-assembled genome, and reciprocal changes in the viral genome. Genome-wide analysis indicated low strain variability of Gordonia, and almost clonal conservation of the trailer end of the CRISPR loci. Incorporation of newer spacers gave rise to multiple coexisting bacterial populations. The host population carrying a shorter CRISPR locus that contain only ancestral spacers, which has not acquired newer spacers against the coexisting phages, accounted for more than half of the total host abundance in the majority of samples. Phages genome co-evolved by introducing directional changes, with no preference for mutations within the protospacer and PAM regions. Metagenomic reconstruction of time-resolved variants of host and viral genomes revealed how the complexity at the population level has important consequences for bacteria-phage coexistence.Subject terms: Microbial ecology, Metagenomics, Bacteriophages     

Highly diverse flavobacterial phages isolated from North Sea spring blooms

It is generally recognized that phages are a mortality factor for their bacterial hosts. This could be particularly true in spring phytoplankton blooms, which are known to be closely followed by a highly specialized bacterial community. We hypothesized that phages modulate these dense heterotrophic bacteria successions following phytoplankton blooms. In this study, we focused on Flavobacteriia, because they are main responders during these blooms and have an important role in the degradation of polysaccharides. A cultivation-based approach was used, obtaining 44 lytic flavobacterial phages (flavophages), representing twelve new species from two viral realms. Taxonomic analysis allowed us to delineate ten new phage genera and ten new families, from which nine and four, respectively, had no previously cultivated representatives. Genomic analysis predicted various life styles and genomic replication strategies. A likely eukaryote-associated host habitat was reflected in the gene content of some of the flavophages. Detection in cellular metagenomes and by direct-plating showed that part of these phages were actively replicating in the environment during the 2018 spring bloom. Furthermore, CRISPR/Cas spacers and re-isolation during two consecutive years suggested that, at least part of the new flavophages are stable components of the microbial community in the North Sea. Together, our results indicate that these diverse flavophages have the potential to modulate their respective host populations.Subject terms: Bacteriophages, Water microbiology     

Efficient dilution-to-extinction isolation of novel virus–host model systems for fastidious heterotrophic bacteria

Microbes and their associated viruses are key drivers of biogeochemical processes in marine and soil biomes. While viruses of phototrophic cyanobacteria are well-represented in model systems, challenges of isolating marine microbial heterotrophs and their viruses have hampered experimental approaches to quantify the importance of viruses in nutrient recycling. A resurgence in cultivation efforts has improved the availability of fastidious bacteria for hypothesis testing, but this has not been matched by similar efforts to cultivate their associated bacteriophages. Here, we describe a high-throughput method for isolating important virus–host systems for fastidious heterotrophic bacteria that couples advances in culturing of hosts with sequential enrichment and isolation of associated phages. Applied to six monthly samples from the Western English Channel, we first isolated one new member of the globally dominant bacterial SAR11 clade and three new members of the methylotrophic bacterial clade OM43. We used these as bait to isolate 117 new phages, including the first known siphophage-infecting SAR11, and the first isolated phage for OM43. Genomic analyses of 13 novel viruses revealed representatives of three new viral genera, and infection assays showed that the viruses infecting SAR11 have ecotype-specific host ranges. Similar to the abundant human-associated phage ɸCrAss001, infection dynamics within the majority of isolates suggested either prevalent lysogeny or chronic infection, despite a lack of associated genes, or host phenotypic bistability with lysis putatively maintained within a susceptible subpopulation. Broader representation of important virus–host systems in culture collections and genomic databases will improve both our understanding of virus–host interactions, and accuracy of computational approaches to evaluate ecological patterns from metagenomic data.Subject terms: Bacteriophages, Microbial ecology     

Phage response to CRISPR-encoded resistance in Streptococcus thermophilus

Clustered regularly interspaced short palindromic repeats (CRISPR) and their associated genes are linked to a mechanism of acquired resistance against bacteriophages. Bacteria can integrate short stretches of phage-derived sequences (spacers) within CRISPR loci to become phage resistant. In this study, we further characterized the efficiency of CRISPR1 as a phage resistance mechanism in Streptococcus thermophilus. First, we show that CRISPR1 is distinct from previously known phage defense systems and is effective against the two main groups of S. thermophilus phages. Analyses of 30 bacteriophage-insensitive mutants of S. thermophilus indicate that the addition of one new spacer in CRISPR1 is the most frequent outcome of a phage challenge and that the iterative addition of spacers increases the overall phage resistance of the host. The added new spacers have a size of between 29 to 31 nucleotides, with 30 being by far the most frequent. Comparative analysis of 39 newly acquired spacers with the complete genomic sequences of the wild-type phages 2972, 858, and DT1 demonstrated that the newly added spacer must be identical to a region (named proto-spacer) in the phage genome to confer a phage resistance phenotype. Moreover, we found a CRISPR1-specific sequence (NNAGAAW) located downstream of the proto-spacer region that is important for the phage resistance phenotype. Finally, we show through the analyses of 20 mutant phages that virulent phages are rapidly evolving through single nucleotide mutations as well as deletions, in response to CRISPR1.     

CRISPR-Cas systems target a diverse collection of invasive mobile genetic elements in human microbiomes

Background

Bacteria and archaea develop immunity against invading genomes by incorporating pieces of the invaders'' sequences, called spacers, into a clustered regularly interspaced short palindromic repeats (CRISPR) locus between repeats, forming arrays of repeat-spacer units. When spacers are expressed, they direct CRISPR-associated (Cas) proteins to silence complementary invading DNA. In order to characterize the invaders of human microbiomes, we use spacers from CRISPR arrays that we had previously assembled from shotgun metagenomic datasets, and identify contigs that contain these spacers'' targets.

Results

We discover 95,000 contigs that are putative invasive mobile genetic elements, some targeted by hundreds of CRISPR spacers. We find that oral sites in healthy human populations have a much greater variety of mobile genetic elements than stool samples. Mobile genetic elements carry genes encoding diverse functions: only 7% of the mobile genetic elements are similar to known phages or plasmids, although a much greater proportion contain phage- or plasmid-related genes. A small number of contigs share similarity with known integrative and conjugative elements, providing the first examples of CRISPR defenses against this class of element. We provide detailed analyses of a few large mobile genetic elements of various types, and a relative abundance analysis of mobile genetic elements and putative hosts, exploring the dynamic activities of mobile genetic elements in human microbiomes. A joint analysis of mobile genetic elements and CRISPRs shows that protospacer-adjacent motifs drive their interaction network; however, some CRISPR-Cas systems target mobile genetic elements lacking motifs.

Conclusions

We identify a large collection of invasive mobile genetic elements in human microbiomes, an important resource for further study of the interaction between the CRISPR-Cas immune system and invaders.     

Adaptation to abiotic conditions drives local adaptation in bacteria and viruses coevolving in heterogeneous environments

Parasite local adaptation, the greater performance of parasites on their local compared with foreign hosts, has important consequences for the maintenance of diversity and epidemiology. While the abiotic environment may significantly affect local adaptation, most studies to date have failed either to incorporate the effects of the abiotic environment, or to separate them from those of the biotic environment. Here, we tease apart biotic and abiotic components of local adaptation using the bacterium Pseudomonas fluorescens and its viral parasite bacteriophage Φ2. We coevolved replicate populations of bacteria and phages at three different temperatures, and determined their performance against coevolutionary partners from the same and different temperatures. Crucially, we measured performance at different assay temperatures, which allowed us to disentangle adaptation to biotic and abiotic habitat components. Our results show that bacteria and phages are more resistant and infectious, respectively, at the temperature at which they previously coevolved, confirming that local adaptation to abiotic conditions can play a crucial role in determining parasite infectivity and host resistance. Our work underlines the need to assess host–parasite interactions across multiple relevant abiotic environments, and suggests that microbial adaption to local temperatures can create ecological barriers to dispersal across temperature gradients.     

Comparisons of clustered regularly interspaced short palindromic repeats and viromes in human saliva reveal bacterial adaptations to salivary viruses

Explorations of human microbiota have provided substantial insight into microbial community composition; however, little is known about interactions between various microbial components in human ecosystems. In response to the powerful impact of viral predation, bacteria have acquired potent defences, including an adaptive immune response based on the clustered regularly interspaced short palindromic repeats (CRISPRs)/Cas system. To improve our understanding of the interactions between bacteria and their viruses in humans, we analysed 13?977 streptococcal CRISPR sequences and compared them with 2?588?172 virome reads in the saliva of four human subjects over 17 months. We found a diverse array of viruses and CRISPR spacers, many of which were specific to each subject and time point. There were numerous viral sequences matching CRISPR spacers; these matches were highly specific for salivary viruses. We determined that spacers and viruses coexist at the same time, which suggests that streptococcal CRISPR/Cas systems are under constant pressure from salivary viruses. CRISPRs in some subjects were just as likely to match viral sequences from other subjects as they were to match viruses from the same subject. Because interactions between bacteria and viruses help to determine the structure of bacterial communities, CRISPR-virus analyses are likely to provide insight into the forces shaping the human microbiome.     

Cytotoxic Chromosomal Targeting by CRISPR/Cas Systems Can Reshape Bacterial Genomes and Expel or Remodel Pathogenicity Islands

In prokaryotes, clustered regularly interspaced short palindromic repeats (CRISPRs) and their associated (Cas) proteins constitute a defence system against bacteriophages and plasmids. CRISPR/Cas systems acquire short spacer sequences from foreign genetic elements and incorporate these into their CRISPR arrays, generating a memory of past invaders. Defence is provided by short non-coding RNAs that guide Cas proteins to cleave complementary nucleic acids. While most spacers are acquired from phages and plasmids, there are examples of spacers that match genes elsewhere in the host bacterial chromosome. In Pectobacterium atrosepticum the type I-F CRISPR/Cas system has acquired a self-complementary spacer that perfectly matches a protospacer target in a horizontally acquired island (HAI2) involved in plant pathogenicity. Given the paucity of experimental data about CRISPR/Cas–mediated chromosomal targeting, we examined this process by developing a tightly controlled system. Chromosomal targeting was highly toxic via targeting of DNA and resulted in growth inhibition and cellular filamentation. The toxic phenotype was avoided by mutations in the cas operon, the CRISPR repeats, the protospacer target, and protospacer-adjacent motif (PAM) beside the target. Indeed, the natural self-targeting spacer was non-toxic due to a single nucleotide mutation adjacent to the target in the PAM sequence. Furthermore, we show that chromosomal targeting can result in large-scale genomic alterations, including the remodelling or deletion of entire pre-existing pathogenicity islands. These features can be engineered for the targeted deletion of large regions of bacterial chromosomes. In conclusion, in DNA–targeting CRISPR/Cas systems, chromosomal interference is deleterious by causing DNA damage and providing a strong selective pressure for genome alterations, which may have consequences for bacterial evolution and pathogenicity.     

Association between living environment and human oral viral ecology

The human oral cavity has an indigenous microbiota known to include a robust community of viruses. Very little is known about how oral viruses are spread throughout the environment or to which viruses individuals are exposed. We sought to determine whether shared living environment is associated with the composition of human oral viral communities by examining the saliva of 21 human subjects; 11 subjects from different households and 10 unrelated subjects comprising 4 separate households. Although there were many viral homologues shared among all subjects studied, there were significant patterns of shared homologues in three of the four households that suggest shared living environment affects viral community composition. We also examined CRISPR (clustered regularly interspaced short palindromic repeat) loci, which are involved in acquired bacterial and archaeal resistance against invading viruses by acquiring short viral sequences. We analyzed 2 065 246 CRISPR spacers from 5 separate repeat motifs found in oral bacterial species of Gemella, Veillonella, Leptotrichia and Streptococcus to determine whether individuals from shared living environments may have been exposed to similar viruses. A significant proportion of CRISPR spacers were shared within subjects from the same households, suggesting either shared ancestry of their oral microbiota or similar viral exposures. Many CRISPR spacers matched virome sequences from different subjects, but no pattern specific to any household was found. Our data on viromes and CRISPR content indicate that shared living environment may have a significant role in determining the ecology of human oral viruses.     

Foreign DNA acquisition by the I-F?CRISPR–Cas system requires all components of the interference machinery

CRISPR immunity depends on acquisition of fragments of foreign DNA into CRISPR arrays. For type I-E CRISPR–Cas systems two modes of spacer acquisition, naïve and primed adaptation, were described. Naïve adaptation requires just two most conserved Cas1 and Cas2 proteins; it leads to spacer acquisition from both foreign and bacterial DNA and results in multiple spacers incapable of immune response. Primed adaptation requires all Cas proteins and a CRISPR RNA recognizing a partially matching target. It leads to selective acquisition of spacers from DNA molecules recognized by priming CRISPR RNA, with most spacers capable of protecting the host. Here, we studied spacer acquisition by a type I-F CRISPR–Cas system. We observe both naïve and primed adaptation. Both processes require not just Cas1 and Cas2, but also intact Csy complex and CRISPR RNA. Primed adaptation shows a gradient of acquisition efficiency as a function of distance from the priming site and a strand bias that is consistent with existence of single-stranded adaption intermediates. The results provide new insights into the mechanism of spacer acquisition and illustrate surprising mechanistic diversity of related CRISPR–Cas systems.     

Potential use of host-derived genome signatures to root virus phylogenies

Two methods are presented that provide independent evidence with which to test virus-phylogeny roots. The methods may be applied to phages and other viruses in which at least one horizontal transfer between hosts is inferred across the virus phylogeny. The methods are based upon the inference that viral DNA sequences acquire similar genome signatures (or “genomic signatures”) to those of their hosts. One of the two methods may be applied to horizontal virus transfers between three or more hosts. Both methods may potentially be extended to rooting plasmid and transposable-element phylogenies. The effect of using different word lengths (2-bp, 3-bp, or 4-bp) to calculate genome signatures, which host genomes (of bacteria) or genome regions (of eukaryotes) were used, which virus sequences were used, and whether distance, counts, or Z-scores are applied were examined using empirical datasets.     

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.     

Top-down controls on bacterial community structure: microbial network analysis of bacteria,T4-like viruses and protists

Characterizing ecological relationships between viruses, bacteria and protists in the ocean are critical to understanding ecosystem function, yet these relationships are infrequently investigated together. We evaluated these relationships through microbial association network analysis of samples collected approximately monthly from March 2008 to January 2011 in the surface ocean (0–5 m) at the San Pedro Ocean Time series station. Bacterial, T4-like myoviral and protistan communities were described by Automated Ribosomal Intergenic Spacer Analysis and terminal restriction fragment length polymorphism of the gene encoding the major capsid protein (g23) and 18S ribosomal DNA, respectively. Concurrent shifts in community structure suggested similar timing of responses to environmental and biological parameters. We linked T4-like myoviral, bacterial and protistan operational taxonomic units by local similarity correlations, which were then visualized as association networks. Network links (correlations) potentially represent synergistic and antagonistic relationships such as viral lysis, grazing, competition or other interactions. We found that virus–bacteria relationships were more cross-linked than protist–bacteria relationships, suggestive of increased taxonomic specificity in virus–bacteria relationships. We also found that 80% of bacterial–protist and 74% of bacterial–viral correlations were positive, with the latter suggesting that at monthly and seasonal timescales, viruses may be following their hosts more often than controlling host abundance.     

Diverse CRISPRs evolving in human microbiomes

CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats) loci, together with cas (CRISPR-associated) genes, form the CRISPR/Cas adaptive immune system, a primary defense strategy that eubacteria and archaea mobilize against foreign nucleic acids, including phages and conjugative plasmids. Short spacer sequences separated by the repeats are derived from foreign DNA and direct interference to future infections. The availability of hundreds of shotgun metagenomic datasets from the Human Microbiome Project (HMP) enables us to explore the distribution and diversity of known CRISPRs in human-associated microbial communities and to discover new CRISPRs. We propose a targeted assembly strategy to reconstruct CRISPR arrays, which whole-metagenome assemblies fail to identify. For each known CRISPR type (identified from reference genomes), we use its direct repeat consensus sequence to recruit reads from each HMP dataset and then assemble the recruited reads into CRISPR loci; the unique spacer sequences can then be extracted for analysis. We also identified novel CRISPRs or new CRISPR variants in contigs from whole-metagenome assemblies and used targeted assembly to more comprehensively identify these CRISPRs across samples. We observed that the distributions of CRISPRs (including 64 known and 86 novel ones) are largely body-site specific. We provide detailed analysis of several CRISPR loci, including novel CRISPRs. For example, known streptococcal CRISPRs were identified in most oral microbiomes, totaling ～8,000 unique spacers: samples resampled from the same individual and oral site shared the most spacers; different oral sites from the same individual shared significantly fewer, while different individuals had almost no common spacers, indicating the impact of subtle niche differences on the evolution of CRISPR defenses. We further demonstrate potential applications of CRISPRs to the tracing of rare species and the virus exposure of individuals. This work indicates the importance of effective identification and characterization of CRISPR loci to the study of the dynamic ecology of microbiomes.     

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.     

前噬菌体研究进展

噬菌体是地球上最多的生物实体,一直被认为是细菌的天敌.然而随着基因组学和分子生物学等技术的快速发展,人们发现噬菌体与宿主之间存在微妙而复杂的关系.前噬菌体是指溶原性细菌内存在的整套噬菌体DNA基因组,广泛分布在细菌基因组中,对调节细菌宿主生理具有重要作用,如参与调节宿主的毒力、影响生物膜形成、赋予宿主免疫力等.有趣的是...     

CRISPR–Cas9 gRNA efficiency prediction: an overview of predictive tools and the role of deep learning

The clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated protein 9 (Cas9) system has become a successful and promising technology for gene-editing. To facilitate its effective application, various computational tools have been developed. These tools can assist researchers in the guide RNA (gRNA) design process by predicting cleavage efficiency and specificity and excluding undesirable targets. However, while many tools are available, assessment of their application scenarios and performance benchmarks are limited. Moreover, new deep learning tools have been explored lately for gRNA efficiency prediction, but have not been systematically evaluated. Here, we discuss the approaches that pertain to the on-target activity problem, focusing mainly on the features and computational methods they utilize. Furthermore, we evaluate these tools on independent datasets and give some suggestions for their usage. We conclude with some challenges and perspectives about future directions for CRISPR–Cas9 guide design.     

Host population diversity as a driver of viral infection cycle in wild populations of green sulfur bacteria with long standing virus-host interactions

Temperate phages are viruses of bacteria that can establish two types of infection: a lysogenic infection in which the virus replicates with the host cell without producing virions, and a lytic infection where the host cell is eventually destroyed, and new virions are released. While both lytic and lysogenic infections are routinely observed in the environment, the ecological and evolutionary processes regulating these viral dynamics are still not well understood, especially for uncultivated virus-host pairs. Here, we characterized the long-term dynamics of uncultivated viruses infecting green sulfur bacteria (GSB) in a model freshwater lake (Trout Bog Lake, TBL). As no GSB virus has been formally described yet, we first used two complementary approaches to identify new GSB viruses from TBL; one in vitro based on flow cytometry cell sorting, the other in silico based on CRISPR spacer sequences. We then took advantage of existing TBL metagenomes covering the 2005–2018 period to examine the interactions between GSB and their viruses across years and seasons. From our data, GSB populations in TBL were constantly associated with at least 2-8 viruses each, including both lytic and temperate phages. The dominant GSB population in particular was consistently associated with two prophages with a nearly 100% infection rate for >10 years. We illustrate with a theoretical model that such an interaction can be stable given a low, but persistent, level of prophage induction in low-diversity host populations. Overall, our data suggest that lytic and lysogenic viruses can readily co-infect the same host population, and that host strain-level diversity might be an important factor controlling virus-host dynamics including lytic/lysogeny switch.Subject terms: Bacteriophages, Metagenomics     

Prophages are associated with extensive CRISPR–Cas auto-immunity

CRISPR–Cas systems require discriminating self from non-self DNA during adaptation and interference. Yet, multiple cases have been reported of bacteria containing self-targeting spacers (STS), i.e. CRISPR spacers targeting protospacers on the same genome. STS has been suggested to reflect potential auto-immunity as an unwanted side effect of CRISPR–Cas defense, or a regulatory mechanism for gene expression. Here we investigated the incidence, distribution, and evasion of STS in over 100 000 bacterial genomes. We found STS in all CRISPR–Cas types and in one fifth of all CRISPR-carrying bacteria. Notably, up to 40% of I-B and I-F CRISPR–Cas systems contained STS. We observed that STS-containing genomes almost always carry a prophage and that STS map to prophage regions in more than half of the cases. Despite carrying STS, genetic deterioration of CRISPR–Cas systems appears to be rare, suggesting a level of escape from the potentially deleterious effects of STS by other mechanisms such as anti-CRISPR proteins and CRISPR target mutations. We propose a scenario where it is common to acquire an STS against a prophage, and this may trigger more extensive STS buildup by primed spacer acquisition in type I systems, without detrimental autoimmunity effects as mechanisms of auto-immunity evasion create tolerance to STS-targeted prophages.