Diverse CRISPRs evolving in human microbiomes

Temporal turnover of viral biodiversity and functional potential in intertidal wetlands

As the central members of the microbiome networks, viruses regulate the composition of microbial communities and drive the nutrient cycles of ecosystems by lysing host cells. Therefore, uncovering the dynamic patterns and the underlying ecological mechanisms mediating the tiniest viral communities across space and through time in natural ecosystems is of crucial importance for better understanding the complex microbial world. Here, the temporal dynamics of intertidal viral communities were investigated via a time-series sampling effort. A total of 1911 viral operational taxonomic units were recovered from 36 bimonthly collected shotgun metagenomes. Functionally important auxiliary metabolic genes involved in carbohydrate, sulfur, and phosphorus metabolism were detected, some of which (e.g., cysH gene) were stably present within viral genomes over time. Over the sampling period, strong and comparable temporal turnovers were observed for intertidal viromes and their host microbes. Winter was determined as the pivotal point for the shifts in viral diversity patterns. Notably, the viral micro-diversity covaried with the macro-diversity, following similar temporal patterns. The relative abundances of viral taxa also covaried with their host prokaryotes. Meanwhile, the virus-host relationships at the whole community level were relatively stable. Further statistical analyses demonstrated that the dynamic patterns of viral communities were highly deterministic, for which temperature was the major driver. This study provided valuable mechanistic insights into the temporal turnover of viral communities in complex ecosystems such as intertidal wetlands.

Virus-prokaryote infection pairs associated with prokaryotic production in a freshwater lake

Viruses infect and kill prokaryotic populations in a density- or frequency-dependent manner and affect carbon cycling. However, the effects of the stratification transition, including the stratified and de-stratified periods, on the changes in prokaryotic and viral communities and their interactions remain unclear. We conducted a monthly survey of the surface and deep layers of a large and deep freshwater lake (Lake Biwa, Japan) for a year and analyzed the prokaryotic production and prokaryotic and viral community composition. Our analysis revealed that, in the surface layer, 19 prokaryotic species, accounting for approximately 40% of the total prokaryotic abundance, could potentially contribute to the majority of prokaryotic production, which is the highest during the summer and is suppressed by viruses. This suggests that a small fraction of prokaryotes and phages were the key infection pairs during the peak period of prokaryotic activity in the freshwater lake. We also found that approximately 50% of the dominant prokaryotic and viral species in the deep layer were present throughout the study period. This suggests that the "kill the winner" model could explain the viral impact on prokaryotes in the surface layer, but other dynamics may be at play in the deep layer. Furthermore, we found that annual vertical mixing could result in a similar rate of community change between the surface and deep layers. These findings may be valuable in understanding how communities and the interaction among them change when freshwater lake stratification is affected by global warming in the future.IMPORTANCEViral infection associated with prokaryotic production occurs in a density- or frequency-dependent manner and regulates the prokaryotic community. Stratification transition and annual vertical mixing in freshwater lakes are known to affect the prokaryotic community and the interaction between prokaryotes and viruses. By pairing measurements of virome analysis and prokaryotic production of a 1-year survey of the depths of surface and deep layers, we revealed (i) the prokaryotic infection pairs associated with prokaryotic production and (ii) the reset in prokaryotic and viral communities through annual vertical mixing in a freshwater lake. Our results provide a basis for future work into changes in stratification that may impact the biogeochemical cycling in freshwater lakes.

Eco-evolutionary dynamics of gut phageome in wild gibbons (Hoolock tianxing) with seasonal diet variations

It has been extensively studied that the gut microbiome provides animals flexibility to adapt to food variability. Yet, how gut phageome responds to diet variation of wild animals remains unexplored. Here, we analyze the eco-evolutionary dynamics of gut phageome in six wild gibbons (Hoolock tianxing) by collecting individually-resolved fresh fecal samples and parallel feeding behavior data for 15 consecutive months. Application of complementary viral and microbial metagenomics recovers 39,198 virulent and temperate phage genomes from the feces. Hierarchical cluster analyses show remarkable seasonal diet variations in gibbons. From high-fruit to high-leaf feeding period, the abundances of phage populations are seasonally fluctuated, especially driven by the increased abundance of virulent phages that kill the Lachnospiraceae hosts, and a decreased abundance of temperate phages that piggyback the Bacteroidaceae hosts. Functional profiling reveals an enrichment through horizontal gene transfers of toxin-antitoxin genes on temperate phage genomes in high-leaf season, potentially conferring benefits to their prokaryotic hosts. The phage-host ecological dynamics are driven by the coevolutionary processes which select for tail fiber and DNA primase genes on virulent and temperate phage genomes, respectively. Our results highlight complex phageome-microbiome interactions as a key feature of the gibbon gut microbial ecosystem responding to the seasonal diet.

Macroalgal virosphere assists with host-microbiome equilibrium regulation and affects prokaryotes in surrounding marine environments

The microbiomes in macroalgal holobionts play vital roles in regulating macroalgal growth and ocean carbon cycling. However, the virospheres in macroalgal holobionts remain largely underexplored, representing a critical knowledge gap. Here we unveil that the holobiont of kelp (Saccharina japonica) harbors highly specific and unique epiphytic/endophytic viral species, with novelty (99.7% unknown) surpassing even extreme marine habitats (e.g. deep-sea and hadal zones), indicating that macroalgal virospheres, despite being closest to us, are among the least understood. These viruses potentially maintain microbiome equilibrium critical for kelp health via lytic-lysogenic infections and the expression of folate biosynthesis genes. In-situ kelp mesocosm cultivation and metagenomic mining revealed that kelp holobiont profoundly reshaped surrounding seawater and sediment virus-prokaryote pairings through changing surrounding environmental conditions and virus-host migrations. Some kelp epiphytic viruses could even infect sediment autochthonous bacteria after deposition. Moreover, the presence of ample viral auxiliary metabolic genes for kelp polysaccharide (e.g. laminarin) degradation underscores the underappreciated viral metabolic influence on macroalgal carbon cycling. This study provides key insights into understanding the previously overlooked ecological significance of viruses within macroalgal holobionts and the macroalgae-prokaryotes-virus tripartite relationship.

Genomic and transcriptomic insights into complex virus-prokaryote interactions in marine biofilms

Marine biofilms are complex communities of microorganisms that play a crucial ecological role in oceans. Although prokaryotes are the dominant members of these biofilms, little is known about their interactions with viruses. By analysing publicly available and newly sequenced metagenomic data, we identified 2446 virus-prokaryote connections in 84 marine biofilms. Most of these connections were between the bacteriophages in the Uroviricota phylum and the bacteria of Proteobacteria, Cyanobacteria and Bacteroidota. The network of virus-host pairs is complex; a single virus can infect multiple prokaryotic populations or a single prokaryote is susceptible to several viral populations. Analysis of genomes of paired prokaryotes and viruses revealed the presence of 425 putative auxiliary metabolic genes (AMGs), 239 viral genes related to restriction-modification (RM) systems and 38,538 prokaryotic anti-viral defence-related genes involved in 15 defence systems. Transcriptomic evidence from newly established biofilms revealed the expression of viral genes, including AMGs and RM, and prokaryotic defence systems, indicating the active interplay between viruses and prokaryotes. A comparison between biofilms and seawater showed that biofilm prokaryotes have more abundant defence genes than seawater prokaryotes, and the defence gene composition differs between biofilms and the surrounding seawater. Overall, our study unveiled active viruses in natural biofilms and their complex interplay with prokaryotes, which may result in the blooming of defence strategists in biofilms. The detachment of bloomed defence strategists may reduce the infectivity of viruses in seawater and result in the emergence of a novel role of marine biofilms.

Combined virome analysis and metagenomic sequencing to reveal the viral communities and risk of virus-associated antibiotic resistance genes during composting

The issue of antibiotic resistance genes (ARGs) pollution in manure has garnered significant attention, with viruses now being recognized as crucial carriers and disseminators of ARGs. However, the virus-associated ARG profiles and potential health risks in composts are still unclear. In this study, the viral communities and associated ARGs in biogas residue and pig faeces composts were profiled by virome analysis. The viral communities were dominated by Caudovirales, and non-thermophilic viruses were inactivated during composting. The diversity and abundance of ARGs were lower in virome than in metagenome, while ARGs' risk was greater in virome than in metagenome. There were six bacterial genera identified as viral hosts at the genomic level, Pseudomonas and Clostridium carried high-risk ARGs. Virus-associated ARGs in viral hosts had a higher risk rank than non-virus-associated ARGs. Composting reduced the diversity, abundance and risk of viral ARGs. The risk of ARGs in biogas residues was significantly lower than that of pig faeces in the initial period of composting, and the two different substracts equally less harmful after composting. These results revealed that viruses play a non-negligible role in spreading ARGs, posing high risk to environmental and human health.

Spatiotemporal Dynamics of Coastal Viral Community Structure and Potential Biogeochemical Roles Affected by an Ulva prolifera Green Tide

The world's largest macroalgal green tide, caused by Ulva prolifera, has resulted in serious consequences for coastal waters of the Yellow Sea, China. Although viruses are considered to be one of the key factors in controlling microalgal bloom demise, understanding of the relationship between viral communities and the macroalgal green tide is still poor. Here, a Qingdao coastal virome (QDCV) time-series data set was constructed based on the metagenomic analysis of 17 DNA viromes along three coastal stations of the Yellow Sea, covering different stages of the green tide from Julian days 165 to 271. A total of 40,076 viral contigs were detected and clustered into 28,058 viral operational taxonomic units (vOTUs). About 84% of the vOTUs could not be classified, and 62% separated from vOTUs in other ecosystems. Green tides significantly influenced the spatiotemporal dynamics of the viral community structure, diversity, and potential functions. For the classified vOTUs, the relative abundance of Pelagibacter phages declined with the arrival of the bloom and rebounded after the bloom, while Synechococcus and Roseobacter phages increased, although with a time lag from the peak of their hosts. More than 80% of the vOTUs reached peaks in abundance at different specific stages, and the viral peaks were correlated with specific hosts at different stages of the green tide. Most of the viral auxiliary metabolic genes (AMGs) were associated with carbon and sulfur metabolism and showed spatiotemporal dynamics relating to the degradation of the large amount of organic matter released by the green tide. IMPORTANCE To the best of our knowledge, this study is the first to investigate the responses of viruses to the world's largest macroalgal green tide. It revealed the spatiotemporal dynamics of the unique viral assemblages and auxiliary metabolic genes (AMGs) following the variation and degradation of Ulva prolifera. These findings demonstrate a tight coupling between viral assemblages, and prokaryotic and eukaryotic abundances were influenced by the green tide.

Monthly dynamics of microbial communities and variation of nitrogen-cycling genes in an industrial-scale expanded granular sludge bed reactor

Introduction

The expanded granular sludge bed (EGSB) is a major form of anaerobic digestion system during wastewater treatment. Yet, the dynamics of microbial and viral communities and members functioning in nitrogen cycling along with monthly changing physicochemical properties have not been well elucidated.

Methods

Here, by collecting the anaerobic activated sludge samples from a continuously operating industrial-scale EGSB reactor, we conducted 16S rRNA gene amplicon sequencing and metagenome sequencing to reveal the microbial community structure and variation with the ever-changing physicochemical properties along within a year.

Results

We observed a clear monthly variation of microbial community structures, while COD, the ratio of volatile suspended solids (VSS) to total suspended solids (TSS) (VSS/TSS ratio), and temperature were predominant factors in shaping community dissimilarities examined by generalized boosted regression modeling (GBM) analysis. Meanwhile, a significant correlation was found between the changing physicochemical properties and microbial communities (p <0.05). The alpha diversity (Chao1 and Shannon) was significantly higher (p <0.05) in both winter (December, January, and February) and autumn (September, October, and November) with higher organic loading rate (OLR), higher VSS/TSS ratio, and lower temperature, resulting higher biogas production and nutrition removal efficiency. Further, 18 key genes covering nitrate reduction, denitrification, nitrification, and nitrogen fixation pathways were discovered, the total abundance of which was significantly associated with the changing environmental factors (p <0.05). Among these pathways, the dissimilatory nitrate reduction to ammonia (DNRA) and denitrification had the higher abundance contributed by the top highly abundant genes narGH, nrfABCDH, and hcp. The COD, OLR, and temperature were primary factors in affecting DNRA and denitrification by GBM evaluation. Moreover, by metagenome binning, we found the DNRA populations mainly belonged to Proteobacteria, Planctomycetota, and Nitrospirae, while the denitrifying bacteria with complete denitrification performance were all Proteobacteria. Besides, we detected 3,360 non-redundant viral sequences with great novelty, in which Siphoviridae, Podoviridae, and Myoviridae were dominant viral families. Interestingly, viral communities likewise depicted clear monthly variation and had significant associations with the recovered populations (p <0.05).

Discussion

Our work highlights the monthly variation of microbial and viral communities during the continuous operation of EGSB affected by the predominant changing COD, OLR, and temperature, while DNRA and denitrification pathways dominated in this anaerobic system. The results also provide a theoretical basis for the optimization of the engineered system.

Significant Differences in Planktonic Virus Communities Between "Cellular Fraction" (0.22 ~ 3.0 µm) and "Viral Fraction" (< 0.22 μm) in the Ocean

Compared to free-living viruses (< 0.22 m) in the ocean, planktonic viruses in the "cellular fraction" (0.22 ~ 3.0 μm) are now far less well understood, and the differences between them remain largely unexplored. Here, we revealed that even in the same seawater samples, the "cellular fraction" comprised significantly distinct virus communities from the free virioplankton, with only 13.87% overlap in viral contigs at the species level. Compared to the viral genomes deposited in NCBI RefSeq database, 99% of the assembled viral genomes in the "cellular fraction" represented novel genera. Notably, the assembled (near-) complete viral genomes within the "cellular fraction" were significantly larger than that in the "viral fraction," and the "cellular fraction" contained three times more species of giant viruses or jumbo phages with genomes > 200 kb than the "viral fraction." The longest complete genomes of jumbo phage (~ 252 kb) and giant virus (~ 716 kb) were both detected only in the "cellular fraction." Moreover, a relatively higher proportion of proviruses were predicted within the "cellular fraction" than "viral fraction." Besides the substantial divergence in viral community structure, the different fractions also contained their unique viral auxiliary metabolic genes; e.g., those potentially participating in inorganic carbon fixation in deep sea were detected only in the "cellular-fraction" viromes. In addition, there was a considerable divergence in the community structure of both "cellular fraction" and "viral fraction" viromes between the surface and deep-sea habitats, suggesting that they might have similar environmental adaptation properties. The findings deepen our understanding of the complexity of viral community structure and function in the ocean.

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).

Patterns and ecological drivers of viral communities in acid mine drainage sediments across Southern China

Recent advances in environmental genomics have provided unprecedented opportunities for the investigation of viruses in natural settings. Yet, our knowledge of viral biogeographic patterns and the corresponding drivers is still limited. Here, we perform metagenomic deep sequencing on 90 acid mine drainage (AMD) sediments sampled across Southern China and examine the biogeography of viruses in this extreme environment. The results demonstrate that prokaryotic communities dictate viral taxonomic and functional diversity, abundance and structure, whereas other factors especially latitude and mean annual temperature also impact viral populations and functions. In silico predictions highlight lineage-specific virus-host abundance ratios and richness-dependent virus-host interaction structure. Further functional analyses reveal important roles of environmental conditions and horizontal gene transfers in shaping viral auxiliary metabolic genes potentially involved in phosphorus assimilation. Our findings underscore the importance of both abiotic and biotic factors in predicting the taxonomic and functional biogeographic dynamics of viruses in the AMD sediments.

The biogeography of viral communities in extreme environments remains understudied. Here, the authors use metagenomic sequencing on 90 acid mine drainage sediments sampled across Southern China, showing the predominant effects of prokaryotic communities and the influence of environmental variables on viral taxonomy and function.

Function and Application of the CRISPR-Cas System in the Plant Pathogen Erwinia amylovora

Phage-based biocontrol is an emerging method for managing the plant pathogen Erwinia amylovora. Control of E. amylovora in North America is achieved chiefly through the application of streptomycin and has led to the development of streptomycin resistance. Resistant E. amylovora can be tracked through the analysis of CRISPR spacer sequences. An alternative to antibiotics are bacterial viruses, known as phages, which lyse their hosts during replication to control the bacterial population. Endogenous CRISPR-Cas systems act as phage resistance mechanisms however, preliminary genomic analysis suggests this activity is limited in E. amylovora. This leaves the functionality of the CRISPR-Cas system, any clade-based differences, and the impact which this system may have on phage-based biocontrol in question. In this study, the CRISPR arrays from 127 newly available genomic sequences of E. amylovora were analyzed through a novel bioinformatic pipeline. Through this, the Eastern and Western North American clades were shown to be incompatible with the current PCR-based approaches for tracking E. amylovora given the size and composition of their CRISPR arrays. Two artificial CRISPR arrays were designed to investigate the functionality of the CRISPR-Cas system in E. amylovora. This system was capable of curing a targeted plasmid and providing phage resistance but was not the source of phage resistance observed within the controls. This suggests that while the CRISPR-Cas system is an important defense mechanism for invasive plasmids, an as yet unidentified mechanism is the primary source of phage resistance in E. amylovora. IMPORTANCE Erwinia amylovora is an economically significant agricultural pathogen found throughout the world. In North America, E. amylovora has developed streptomycin resistance and therefore alternative treatments using phages have received increased attention. In this study, we analyzed recently published genomes to determine that two significant groups of E. amylovora are poorly identified using the current, CRISPR-based tracking methods. We also showed that the CRISPR-Cas system and an unidentified mechanism work together to provide a significant degree of resistance against one of the phages proposed for phage-based biocontrol.

Transkingdom Analysis of the Female Reproductive Tract Reveals Bacteriophages form Communities

The female reproductive tract (FRT) microbiome plays a vital role in maintaining vaginal health. Viruses are key regulators of other microbial ecosystems, but little is known about how the FRT viruses (virome), particularly bacteriophages that comprise the phageome, impact FRT health and dysbiosis. We hypothesize that bacterial vaginosis (BV) is associated with altered FRT phageome diversity, transkingdom interplay, and bacteriophage discriminate taxa. Here, we conducted a retrospective, longitudinal analysis of vaginal swabs collected from 54 BV-positive and 46 BV-negative South African women. Bacteriome analysis revealed samples clustered into five distinct bacterial community groups (CGs), and further, bacterial alpha diversity was significantly associated with BV. Virome analysis on a subset of baseline samples showed FRT bacteriophages clustering into novel viral state types (VSTs), a viral community clustering system based on virome composition and abundance. Distinct BV bacteriophage signatures included increased alpha diversity along with discriminant Bacillus, Burkholderia, and Escherichia bacteriophages. Bacteriophage-bacteria transkingdom associations were also identified between Bacillus and Burkholderia viruses and BV-associated bacteria, providing key insights for future studies elucidating the transkingdom interactions driving BV-associated microbiome perturbations. In this cohort, bacteriophage-bacterial associations suggest complex interactions, which may play a role in the establishment and maintenance of BV.

Genomic Analysis of Molecular Bacterial Mechanisms of Resistance to Phage Infection

To optimize phage therapy, we need to understand how bacteria evolve against phage attacks. One of the main problems of phage therapy is the appearance of bacterial resistance variants. The use of genomics to track antimicrobial resistance is increasingly developed and used in clinical laboratories. For that reason, it is important to consider, in an emerging future with phage therapy, to detect and avoid phage-resistant strains that can be overcome by the analysis of metadata provided by whole-genome sequencing. Here, we identified genes associated with phage resistance in 18 Acinetobacter baumannii clinical strains belonging to the ST-2 clonal complex during a decade (Ab2000 vs. 2010): 9 from 2000 to 9 from 2010. The presence of genes putatively associated with phage resistance was detected. Genes detected were associated with an abortive infection system, restriction–modification system, genes predicted to be associated with defense systems but with unknown function, and CRISPR-Cas system. Between 118 and 171 genes were found in the 18 clinical strains. On average, 26% of these genes were detected inside genomic islands in the 2000 strains and 32% in the 2010 strains. Furthermore, 38 potential CRISPR arrays in 17 of 18 of the strains were found, as well as 705 proteins associated with CRISPR-Cas systems. A moderately higher presence of these genes in the strains of 2010 in comparison with those of 2000 was found, especially those related to the restriction–modification system and CRISPR-Cas system. The presence of these genes in genomic islands at a higher rate in the strains of 2010 compared with those of 2000 was also detected. Whole-genome sequencing and bioinformatics could be powerful tools to avoid drawbacks when a personalized therapy is applied. In this study, it allows us to take care of the phage resistance in A. baumannii clinical strains to prevent a failure in possible phage therapy.

The female reproductive tract virome: understanding the dynamic role of viruses in gynecological health and disease

The human body is inhabited by a large and complex network of commensal and predatory eukaryotic viruses and bacteriophages collectively termed the virome. Despite being the most abundant and genetically diverse biological entities on the planet, the impact of viruses on human health especially within the female reproductive tract (FRT) remains understudied. To better appreciate current knowledge regarding the dynamic role of viruses in FRT health and disease, in this review we highlight the known constituents of the FRT virome, transkingdom interactions within the FRT and their influence on gynecological disease. A better understanding of the FRT virome may pave the way toward improved outcomes in gynecological, reproductive, and neonatal health.

Using metagenomic data to boost protein structure prediction and discovery

Over the past decade, metagenomic sequencing approaches have been providing an ever-increasing amount of protein sequence data at an astonishing rate. These constitute an invaluable source of information which has been exploited in various research fields such as the study of the role of the gut microbiota in human diseases and aging. However, only a small fraction of all metagenomic sequences collected have been functionally or structurally characterized, leaving much of them completely unexplored. Here, we review how this information has been used in protein structure prediction and protein discovery. We begin by presenting some widely used metagenomic databases and analyze in detail how metagenomic data has contributed to the impressive improvement in the accuracy of structure prediction methods in recent years. We then examine how metagenomic information can be exploited to annotate protein sequences. More specifically, we focus on the role of metagenomes in the discovery of enzymes and new CRISPR-Cas systems, and in the identification of antibiotic resistance genes. With this review, we provide an overview of how metagenomic data is currently revolutionizing our understanding of protein science.

Saline lakes on the Qinghai-Tibet Plateau harbor unique viral assemblages mediating microbial environmental adaption

The highest plateau on Earth, Qinghai-Tibet Plateau, contains thousands of lakes with broad salinity and diverse and unique microbial communities. However, little is known about their co-occurring viruses. Herein, we identify 4,560 viral Operational Taxonomic Units (vOTUs) from six viromes of three saline lakes on Qinghai-Tibet Plateau, with less than 1% that could be classified. Most of the predicted vOTUs were associated with the dominant bacterial and archaeal phyla. Virus-encoded auxiliary metabolic genes suggest that viruses influence microbial metabolisms of carbon, nitrogen, sulfur, and lipid; the antibiotic resistance mediation; and their salinity adaption. The six viromes clustered together with the ice core viromes and bathypelagic ocean viromes and might represent a new viral habitat. This study has revealed the unique characteristics and potential ecological roles of DNA viromes in the lakes of the highest plateau and established a foundation for the recognition of the viral roles in plateau lake ecosystems.

Rapid screening and identification of viral pathogens in metagenomic data

BACKGROUND

Virus screening and viral genome reconstruction are urgent and crucial for the rapid identification of viral pathogens, i.e., tracing the source and understanding the pathogenesis when a viral outbreak occurs. Next-generation sequencing (NGS) provides an efficient and unbiased way to identify viral pathogens in host-associated and environmental samples without prior knowledge. Despite the availability of software, data analysis still requires human operations. A mature pipeline is urgently needed when thousands of viral pathogen and viral genome reconstruction samples need to be rapidly identified.

RESULTS

In this paper, we present a rapid and accurate workflow to screen metagenomics sequencing data for viral pathogens and other compositions, as well as enable a reference-based assembler to reconstruct viral genomes. Moreover, we tested our workflow on several metagenomics datasets, including a SARS-CoV-2 patient sample with NGS data, pangolins tissues with NGS data, Middle East Respiratory Syndrome (MERS)-infected cells with NGS data, etc. Our workflow demonstrated high accuracy and efficiency when identifying target viruses from large scale NGS metagenomics data. Our workflow was flexible when working with a broad range of NGS datasets from small (kb) to large (100 Gb). This took from a few minutes to a few hours to complete each task. At the same time, our workflow automatically generates reports that incorporate visualized feedback (e.g., metagenomics data quality statistics, host and viral sequence compositions, details about each of the identified viral pathogens and their coverages, and reassembled viral pathogen sequences based on their closest references).

CONCLUSIONS

Overall, our system enabled the rapid screening and identification of viral pathogens from metagenomics data, providing an important piece to support viral pathogen research during a pandemic. The visualized report contains information from raw sequence quality to a reconstructed viral sequence, which allows non-professional people to screen their samples for viruses by themselves (Additional file 1).

Engineered CRISPR-Cas systems for the detection and control of antibiotic-resistant infections

Antibiotic resistance is spreading rapidly around the world and seriously impeding efforts to control microbial infections. Although nucleic acid testing is widely deployed for the detection of antibiotic resistant bacteria, the current techniques-mainly based on polymerase chain reaction (PCR)-are time-consuming and laborious. There is an urgent need to develop new strategies to control bacterial infections and the spread of antimicrobial resistance (AMR). The CRISPR-Cas system is an adaptive immune system found in many prokaryotes that presents attractive opportunities to target and edit nucleic acids with high precision and reliability. Engineered CRISPR-Cas systems are reported to effectively kill bacteria or even revert bacterial resistance to antibiotics (resensitizing bacterial cells to antibiotics). Strategies for combating antimicrobial resistance using CRISPR (i.e., Cas9, Cas12, Cas13, and Cas14) can be of great significance in detecting bacteria and their resistance to antibiotics. This review discusses the structures, mechanisms, and detection methods of CRISPR-Cas systems and how these systems can be engineered for the rapid and reliable detection of bacteria using various approaches, with a particular focus on nanoparticles. In addition, we summarize the most recent advances in applying the CRISPR-Cas system for virulence modulation of bacterial infections and combating antimicrobial resistance.

Comparison of CRISPR-Cas Immune Systems in Healthcare-Related Pathogens

The ESKAPE pathogens (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, and Enterobacter species) and Clostridium difficile have been identified as the leading global cause of multidrug-resistant bacterial infections in hospitals. CRISPR-Cas systems are bacterial immune systems, empowering the bacteria with defense against invasive mobile genetic elements that may carry the antimicrobial resistance (AMR) genes, among others. On the other hand, the CRISPR-Cas systems are themselves mobile. In this study, we annotated and compared the CRISPR-Cas systems in these pathogens, utilizing their publicly available large numbers of sequenced genomes (e.g., there are more than 12 thousands of S. aureus genomes). The presence of CRISPR-Cas systems showed a very broad spectrum in these pathogens: S. aureus has the least tendency of obtaining the CRISPR-Cas systems with only 0.55% of its isolates containing CRISPR-Cas systems, whereas isolates of C. difficile we analyzed have CRISPR-Cas systems each having multiple CRISPRs. Statistical tests show that CRISPR-Cas containing isolates tend to have more AMRs for four of the pathogens (A. baumannii, E. faecium, P. aeruginosa, and S. aureus). We made available all the annotated CRISPR-Cas systems in these pathogens with visualization at a website (https://omics.informatics.indiana.edu/CRISPRone/pathogen), which we believe will be an important resource for studying the pathogens and their arms-race with invaders mediated through the CRISPR-Cas systems, and for developing potential clinical applications of the CRISPR-Cas systems for battles against the antibiotic resistant pathogens.

Arms race in a cell: genomic, transcriptomic, and proteomic insights into intracellular phage-bacteria interplay in deep-sea snail holobionts

BACKGROUND

Deep-sea animals in hydrothermal vents often form endosymbioses with chemosynthetic bacteria. Endosymbionts serve essential biochemical and ecological functions, but the prokaryotic viruses (phages) that determine their fate are unknown.

RESULTS

We conducted metagenomic analysis of a deep-sea vent snail. We assembled four genome bins for Caudovirales phages that had developed dual endosymbiosis with sulphur-oxidising bacteria (SOB) and methane-oxidising bacteria (MOB). Clustered regularly interspaced short palindromic repeat (CRISPR) spacer mapping, genome comparison, and transcriptomic profiling revealed that phages Bin1, Bin2, and Bin4 infected SOB and MOB. The observation of prophages in the snail endosymbionts and expression of the phage integrase gene suggested the presence of lysogenic infection, and the expression of phage structural protein and lysozyme genes indicated active lytic infection. Furthermore, SOB and MOB appear to employ adaptive CRISPR-Cas systems to target phage DNA. Additional expressed defence systems, such as innate restriction-modification systems and dormancy-inducing toxin-antitoxin systems, may co-function and form multiple lines for anti-viral defence. To counter host defence, phages Bin1, Bin2, and Bin3 appear to have evolved anti-restriction mechanisms and expressed methyltransferase genes that potentially counterbalance host restriction activity. In addition, the high-level expression of the auxiliary metabolic genes narGH, which encode nitrate reductase subunits, may promote ATP production, thereby benefiting phage DNA packaging for replication.

CONCLUSIONS

This study provides new insights into phage-bacteria interplay in intracellular environments of a deep-sea vent snail. Video Abstract.

CRISPRclassify: Repeat-Based Classification of CRISPR Loci

Detection and classification of CRISPR-Cas systems in metagenomic data have become increasingly prevalent in recent years due to their potential for diverse applications in genome editing. Traditionally, CRISPR-Cas systems are classified through reference-based identification of proximate cas genes. Here, we present a machine learning approach for the detection and classification of CRISPR loci using repeat sequences in a cas-independent context, enabling identification of unclassified loci missed by traditional cas-based approaches. Using biological attributes of the CRISPR repeat, the core element in CRISPR arrays, and leveraging methods from natural language processing, we developed a machine learning model capable of accurate classification of CRISPR loci in an extensive set of metagenomes, resulting in an F1 measure of 0.82 across all predictions and an F1 measure of 0.97 when limiting to classifications with probabilities >0.85. Furthermore, assessing performance on novel repeats yielded an F1 measure of 0.96. Although the performance of cas-based identification will exceed that of a repeat-based approach in many cases, CRISPRclassify provides an efficient approach to classification of CRISPR loci for cases in which cas gene information is unavailable, such as metagenomes and fragmented genome assemblies.

Immunity, virus evolution, and effectiveness of SARS-CoV-2 vaccines

Phylogenetic and pathogenesis studies of the severe acute respiratory syndrome-related coronaviruses (SARS-CoVs) strains have highlighted some specific mutations that could confer the RNA genome fitness advantages and immunological resistance for their rapid spread in the human population. The analyses of 30 kb RNA SARS-CoVs genome sequences, protein structures, and functions have provided us a perspective of how host-virus protein-protein complexes act to mediate virus infection. The open reading frame (ORF)1a and ORF1b translation yields 16 non-structural (nsp1-16) and 6 accessory proteins (p6, p7a, p8ab, p9b) with multiple functional domains. Viral proteins recruit over 300 host partners forming hetero-oligomeric complexes enabling the viral RNA synthesis, packing, and virion release. Many cellular host factors and the innate immune cells through pattern-recognition receptors and intracellular RNA sensor molecules act to inhibit virus entry and intracellular replication. However, non-structural ORF proteins hijack them and suppress interferon synthesis and its antiviral effects. Pro-inflammatory chemokines and cytokines storm leads to dysfunctional inflammation, lung injury, and several clinical symptoms in patients. During the global pandemic, COVID-19 patients were identified with non-synonymous substitution of G614D in the spike protein, indicating virus co-evolution in host cells. We review findings that suggest that host RNA editing and DNA repair systems, while carrying on recombination, mutation, and repair of viral RNA intermediates, may facilitate virus evolution. Understanding how the host cell RNA replication process may be driven by SARS-CoV-2 RNA genome fitness will help the testing of vaccines effectiveness to multiple independent mutated coronavirus strains that will emerge.

CRISPRbuilder-TB: "CRISPR-builder for tuberculosis". Exhaustive reconstruction of the CRISPR locus in mycobacterium tuberculosis complex using SRA

Mycobacterium tuberculosis complex (MTC) CRISPR locus diversity has long been studied solely investigating the presence/absence of a known set of spacers. Unveiling the genetic mechanisms of its evolution requires a more exhaustive reconstruction in a large amount of representative strains.

In this article, we point out and resolve, with a new pipeline, the problem of CRISPR reconstruction based directly on short read sequences in M. tuberculosis. We first show that the process we set up, that we coin as “CRISPRbuilder-TB” (https://github.com/cguyeux/CRISPRbuilder-TB), allows an efficient reconstruction of simulated or real CRISPRs, even when including complex evolutionary steps like the insertions of mobile elements. Compared to more generalist tools, the whole process is much more precise and robust, and requires only minimal manual investigation. Second, we show that more than 1/3 of the currently complete genomes available for this complex in the public databases contain largely erroneous CRISPR loci. Third, we highlight how both the classical experimental in vitro approach and the basic in silico spoligotyping provided by existing analytic tools miss a whole diversity of this locus in MTC, by not capturing duplications, spacer and direct repeats variants, and IS6110 insertion locations. This description is extended in a second article that describes MTC-CRISPR diversity and suggests general rules for its evolution. This work opens perspectives for an in-depth exploration of M. tuberculosis CRISPR loci diversity and of mechanisms involved in its evolution and its functionality, as well as its adaptation to other CRISPR locus-harboring bacterial species.

In this article, we tackle the bioinformatical issue of the reconstruction of the Mycobacterium tuberculosis complex CRISPR locus using short read sequences without requiring genome assembly. We first show that many complete genomes, as found in public databases and often reconstructed by de novo assemblies, often contain errors on this locus as well as on other repeated sequences. We provide an in-depth description of our new method, designated as ‘CRISPRbuilder-TB’, and we show that our method provides much more exhaustive and reliable information (on DR variants, spacer diversity, global structure) than Crass and CRISPR_detector. The new and unsuspected genomic diversity we detected is described in a companion paper. Scripts are available to adapt the tool to other species.

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.

Evaluation of CRISPR Diversity in the Human Skin Microbiome for Personal Identification

Microbial community diversity analysis can be utilized to characterize the personal microbiome that varies between individuals. CRISPR sequences, which reflect virome structure, in the human skin environment may be highly personalized similar to the structures of individual viromes.

The highly personalized human skin microbiome may serve as a viable marker in personal identification. Amplicon sequencing resolution using 16S rRNA cannot identify bacterial communities sufficiently to discriminate between individuals. Thus, novel higher-resolution genetic markers are required for forensic purposes. The clustered regularly interspaced short palindromic repeats (CRISPRs) are prokaryotic genetic elements that can provide a history of infections encountered by the bacteria. The sequencing of CRISPR spacers may provide phylogenetic information with higher resolution than other markers. However, using spacer sequencing for discrimination of personal skin microbiome is difficult due to limited information on CRISPRs in human skin microbiomes. It remains unclear whether personal microbiome discrimination can be achieved using spacer diversity or which CRISPRs will be forensically relevant. We identified common CRISPRs in the human skin microbiome via metagenomic reconstruction and used amplicon sequencing for deep sequencing of spacers. We successfully reconstructed 24 putative CRISPR arrays using metagenomic data sets. A total of 1,223,462 reads from three CRISPR arrays revealed that spacers in the skin microbiome were highly personalized, and conserved repeats were commonly shared between individuals. These individual specificities observed using CRISPR typing were confirmed by comparing the CRISPR diversity to microbiome diversity assessed using 16S rRNA amplicon sequencing. CRISPR typing achieved 95.2% accuracy in personal classification, whereas 16S rRNA sequencing only achieved 52.6%. These results suggest that sequencing CRISPRs in the skin microbiome may be a more powerful approach for personal identification and ecological studies compared to conventional 16S rRNA sequencing.

IMPORTANCE Microbial community diversity analysis can be utilized to characterize the personal microbiome that varies between individuals. CRISPR sequences, which reflect virome structure, in the human skin environment may be highly personalized similar to the structures of individual viromes. In this study, we identified 24 putative CRISPR arrays using a shotgun metagenome data set of the human skin microbiome. The findings of this study expand our understanding of the nature of CRISPRs by identifying novel CRISPR candidates. We developed a method to efficiently determine the diversity of three CRISPR arrays. Our analysis revealed that the CRISPR spacer diversity in the human skin microbiome is highly personalized compared with the microbiome diversity assessed by 16S rRNA sequencing, providing a new perspective on the study of the skin microbiome.

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

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

Identification of Natural CRISPR Systems and Targets in the Human Microbiome

Many bacteria resist invasive DNA by incorporating sequences into CRISPR loci, which enable sequence-specific degradation. CRISPR systems have been well studied from isolate genomes, but culture-independent metagenomics provide a new window into their diversity. We profiled CRISPR loci and cas genes in the body-wide human microbiome using 2,355 metagenomes, yielding functional and taxonomic profiles for 2.9 million spacers by aligning the spacer content to each sample's metagenome and corresponding gene families. Spacer and repeat profiles agree qualitatively with those from isolate genomes but expand their diversity by approximately 13-fold, with the highest spacer load present in the oral microbiome. The taxonomy of spacer sequences parallels that of their source community, with functional targets enriched for viral elements. When coupled with cas gene systems, CRISPR-Cas subtypes are highly site and taxon specific. Our analysis provides a comprehensive collection of natural CRISPR-cas loci and targets in the human microbiome.

In silico Method in CRISPR/Cas System: An Expedite and Powerful Booster

The CRISPR/Cas system has stood in the center of attention in the last few years as a revolutionary gene editing tool with a wide application to investigate gene functions. However, the labor-intensive workflow requires a sophisticated pre-experimental and post-experimental analysis, thus becoming one of the hindrances for the further popularization of practical applications. Recently, the increasing emergence and advancement of the in silico methods play a formidable role to support and boost experimental work. However, various tools based on distinctive design principles and frameworks harbor unique characteristics that are likely to confuse users about how to choose the most appropriate one for their purpose. In this review, we will present a comprehensive overview and comparisons on the in silico methods from the aspects of CRISPR/Cas system identification, guide RNA design, and post-experimental assistance. Furthermore, we establish the hypotheses in light of the new trends around the technical optimization and hope to provide significant clues for future tools development.

Depth-related variability in viral communities in highly stratified sulfidic mine tailings

BACKGROUND

Recent studies have significantly expanded our knowledge of viral diversity and functions in the environment. Exploring the ecological relationships between viruses, hosts, and the environment is a crucial first step towards a deeper understanding of the complex and dynamic interplays among them.

RESULTS

Here, we obtained extensive 16S rRNA gene amplicon, metagenomics sequencing, and geochemical datasets from different depths of two highly stratified sulfidic mine tailings cores with steep geochemical gradients especially pH, and explored how variations in viral community composition and functions were coupled to the co-existing prokaryotic assemblages and the varying environmental conditions. Our data showed that many viruses in the mine tailings represented novel genera, based on gene-sharing networks. Siphoviridae, Podoviridae, and Myoviridae dominated the classified viruses in the surface tailings and deeper layers. Both viral richness and normalized coverage increased with depth in the tailings cores and were significantly correlated with geochemical properties, for example, pH. Viral richness was also coupled to prokaryotic richness (Pearson's r = 0.65, P = 0.032). The enrichment of prophages in the surface mine tailings suggested a preference of lysogenic viral lifestyle in more acidic conditions. Community-wide comparative analyses clearly showed that viruses in the surface tailings encoded genes mostly with unknown functions while viruses in the deeper layers contained genes mainly annotated as conventional functions related to metabolism and structure. Notably, significantly abundant assimilatory sulfate reduction genes were identified from the deeper tailings layers and they were widespread in viruses predicted to infect diverse bacterial phyla.

CONCLUSIONS

Overall, our results revealed a depth-related distribution of viral populations in the extreme and heterogeneous tailings system. The viruses may interact with diverse hosts and dynamic environmental conditions and likely play a role in the functioning of microbial community and modulate sulfur cycles in situ. Video Abstract.

CRISPR-Cas bioinformatics

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

Comparative metagenomic analysis of microbial taxonomic and functional variations in untreated surface and reclaimed waters used in irrigation applications

The use of irrigation water sourced from reclamation facilities and untreated surface water bodies may be a practical solution to attenuate the burden on diminishing groundwater aquifers. However, comprehensive microbial characterizations of these water sources are generally lacking, especially with regard to variations through time and across multiple water types. To address this knowledge gap we used a shotgun metagenomic approach to characterize the taxonomic and functional variations of microbial communities within two agricultural ponds, two freshwater creeks, two brackish rivers, and three water reclamation facilities located in the Mid-Atlantic, United States. Water samples (n = 24) were collected from all sites between October and November 2016, and filtered onto 0.2 μm membrane filters. Filters were then subjected to total DNA extraction and shotgun sequencing on the Illumina HiSeq platform. From these data, we found that Betaproteobacteria dominated the majority of freshwater sites, while Alphaproteobacteria were abundant at times in the brackish waters. One of these brackish sites was also host to a greater abundance of the bacterial genera Gimesia and Microcystis. Furthermore, predicted microbial features (e.g. antibiotic resistance genes (ARGs) and Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) arrays) varied based on specific site and sampling date. ARGs were found across samples, with the diversity and abundance highest in those from a reclamation facility and a wastewater-impacted freshwater creek. Additionally, we identified over 600 CRISPR arrays, containing ∼2600 unique spacers, suggestive of a diverse and often site-specific phage community. Overall, these results provide a better understanding of the complex microbial community in untreated surface and reclaimed waters, while highlighting possible environmental and human health impacts associated with their use in agriculture.

Prevalence of viral photosynthesis genes along a freshwater to saltwater transect in Southeast USA

Bacteriophages encode host-acquired functional genes known as auxiliary metabolic genes (AMGs). Photosynthesis AMGs are commonly found in marine cyanobacteria-infecting Myoviridae and Podoviridae cyanophages, but their ecology remains understudied in freshwater environments. To advance knowledge of this issue, we analysed viral metagenomes collected in the summertime for four years from five lakes and two estuarine locations interconnected by the Chattahoochee River, Southeast USA. Sequences representing ten different AMGs were recovered and found to be prevalent in all sites. Most freshwater AMGs were 10-fold less abundant than estuarine and marine AMGs and were encoded by novel Myoviridae and Podoviridae cyanophage genera. Notably, several of the corresponding viral genomes showed endemism to a specific province along the river. This translated into psbA gene phylogenetic clustering patterns that matched a marine vs. freshwater origin indicating that psbA may serve as a robust classification and source-tracking biomarker. Genomes classified in a novel viral lineage represented by isolate S-EIVl contained psbA, which is unprecedented for this lineage. Collectively, our findings indicated that the acquisition of photosynthesis AMGs is a widespread strategy used by cyanophages in aquatic ecosystems, and further indicated the existence of viral provinces in which certain viral species and/or genotypes are locally abundant.

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.

Phages and Human Health: More Than Idle Hitchhikers

Bacteriophages, or phages, are viruses that infect bacteria and archaea. Phages have diverse morphologies and can be coded in DNA or RNA and as single or double strands with a large range of genome sizes. With the increasing use of metagenomic sequencing approaches to analyze complex samples, many studies generate massive amounts of "viral dark matter", or sequences of viral origin unable to be classified either functionally or taxonomically. Metagenomic analysis of phages is still in its infancy, and uncovering novel phages continues to be a challenge. Work over the past two decades has begun to uncover key roles for phages in different environments, including the human gut. Recent studies in humans have identified expanded phage populations in both healthy infants and in inflammatory bowel disease patients, suggesting distinct phage activity during development and in specific disease states. In this review, we examine our current knowledge of phage biology and discuss recent efforts to improve the analysis and discovery of novel phages. We explore the roles phages may play in human health and disease and discuss the future of phage research.

The Integrative Human Microbiome Project

The NIH Human Microbiome Project (HMP) has been carried out over ten years and two phases to provide resources, methods, and discoveries that link interactions between humans and their microbiomes to health-related outcomes. The recently completed second phase, the Integrative Human Microbiome Project, comprised studies of dynamic changes in the microbiome and host under three conditions: pregnancy and preterm birth; inflammatory bowel diseases; and stressors that affect individuals with prediabetes. The associated research begins to elucidate mechanisms of host-microbiome interactions under these conditions, provides unique data resources (at the HMP Data Coordination Center), and represents a paradigm for future multi-omic studies of the human microbiome.

The Integrative Human Microbiome Project

The NIH Human Microbiome Project (HMP) has been carried out over ten years and two phases to provide resources, methods, and discoveries that link interactions between humans and their microbiomes to health-related outcomes. The recently completed second phase, the Integrative Human Microbiome Project, comprised studies of dynamic changes in the microbiome and host under three conditions: pregnancy and preterm birth; inflammatory bowel diseases; and stressors that affect individuals with prediabetes. The associated research begins to elucidate mechanisms of host-microbiome interactions under these conditions, provides unique data resources (at the HMP Data Coordination Center), and represents a paradigm for future multi-omic studies of the human microbiome.

CRISPR Spacers Indicate Preferential Matching of Specific Virioplankton Genes

Viral infection exerts selection pressure on marine microbes, as virus-induced cell lysis causes 20 to 50% of cell mortality, resulting in fluxes of biomass into oceanic dissolved organic matter. Archaeal and bacterial populations can defend against viral infection using the clustered regularly interspaced short palindromic repeat (CRISPR)-associated (Cas) system, which relies on specific matching between a spacer sequence and a viral gene. If a CRISPR spacer match to any gene within a viral genome is equally effective in preventing lysis, no viral genes should be preferentially matched by CRISPR spacers. However, if there are differences in effectiveness, certain viral genes may demonstrate a greater frequency of CRISPR spacer matches. Indeed, homology search analyses of bacterioplankton CRISPR spacer sequences against virioplankton sequences revealed preferential matching of replication proteins, nucleic acid binding proteins, and viral structural proteins. Positive selection pressure for effective viral defense is one parsimonious explanation for these observations. CRISPR spacers from virioplankton metagenomes preferentially matched methyltransferase and phage integrase genes within virioplankton sequences. These virioplankton CRISPR spacers may assist infected host cells in defending against competing phage. Analyses also revealed that half of the spacer-matched viral genes were unknown, some genes matched several spacers, and some spacers matched multiple genes, a many-to-many relationship. Thus, CRISPR spacer matching may be an evolutionary algorithm, agnostically identifying those genes under stringent selection pressure for sustaining viral infection and lysis. Investigating this subset of viral genes could reveal those genetic mechanisms essential to virus-host interactions and provide new technologies for optimizing CRISPR defense in beneficial microbes.IMPORTANCE The CRISPR-Cas system is one means by which bacterial and archaeal populations defend against viral infection which causes 20 to 50% of cell mortality in the ocean. We tested the hypothesis that certain viral genes are preferentially targeted for the initial attack of the CRISPR-Cas system on a viral genome. Using CASC, a pipeline for CRISPR spacer discovery, and metagenome data from oceanic microbes and viruses, we found a clear subset of viral genes with high match frequencies to CRISPR spacers. Moreover, we observed a many-to-many relationship of spacers and viral genes. These high-match viral genes were involved in nucleotide metabolism, DNA methylation, and viral structure. It is possible that CRISPR spacer matching is an evolutionary algorithm pointing to those viral genes most important to sustaining infection and lysis. Studying these genes may advance the understanding of virus-host interactions in nature and provide new technologies for leveraging CRISPR-Cas systems in beneficial microbes.

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.

Predicting and visualizing features of CRISPR-Cas systems

Pervasive application of CRISPR-Cas systems in genome editing has prompted an increase in both interest and necessity to further elucidate existing systems as well as discover putative novel systems. The ubiquity and power of current computational platforms have made in silico approaches to CRISPR-Cas identification and characterization accessible to a wider audience and increasingly amenable for processing extensive data sets. Here, we describe in silico methods for predicting and visualizing notable features of CRISPR-Cas systems, including Cas domain determination, CRISPR array visualization, and inference of the protospacer-adjacent motif. The efficiency of these tools enables rapid exploration of CRISPR-Cas diversity across prokaryotic genomes and supports scalable analysis of large genomic data sets.

Le système CRISPR-Cas

CRISPR-Cas est un système immunitaire adaptatif utilisé par de nombreux microbes pour se défendre contre l’invasion d’acides nucléiques tels que les génomes viraux et autres éléments génétiques mobiles. Le système microbien utilise son locus CRISPR pour stocker de l’information génétique afin de produire des ARN guides. Ces derniers, de concert avec des endonucléases (Cas), empêchent des invasions futures. Des parties de ce système microbien ont été exploitées pour développer un puissant outil d’édition des génomes dans une panoplie d’organismes. La capacité de CRISPR-Cas9 à couper efficacement et à des endroits très précis de l’ADN pourrait peut-être permettre un jour de guérir certaines maladies génétiques humaines. La malléabilité de cet outil d’édition rend possible une variété d’applications allant de la modulation de l’expression de gènes à des modifications épigénétiques. Les locus CRISPR représentent également une mine d’informations pouvant servir de méthode de typage de souches microbiennes ou encore une façon d’étudier les interactions entre les bactéries et leurs habitats.

Genome Editing: A New Horizon for Oral and Craniofacial Research

Precise and efficient genetic manipulations have enabled researchers to understand gene functions in disease and development, providing a platform to search for molecular cures. Over the past decade, the unprecedented advancement of genome editing techniques has revolutionized the biological research fields. Early genome editing strategies involved many naturally occurring nucleases, including meganucleases, zinc finger nucleases, and transcription activator-like effector-based nucleases. More recently, the clustered regularly interspaced short palindromic repeats (CRISPR) / CRISPR-associated nucleases (CRISPR/Cas) system has greatly enriched genetic manipulation methods in conducting research. Those nucleases generate double-strand breaks in the target gene sequences and then utilize DNA repair mechanisms to permit precise yet versatile genetic manipulations. The oral and craniofacial field harbors a plethora of diseases and developmental defects that require genetic models that can exploit these genome editing techniques. This review provides an overview of the genome editing techniques, particularly the CRISPR/Cas9 technique, for the oral and craniofacial research community. We also discuss the details about the emerging applications of genome editing in oral and craniofacial biology.

Metagenomic analysis with strain-level resolution reveals fine-scale variation in the human pregnancy microbiome

Recent studies suggest that the microbiome has an impact on gestational health and outcome. However, characterization of the pregnancy-associated microbiome has largely relied on 16S rRNA gene amplicon-based surveys. Here, we describe an assembly-driven, metagenomics-based, longitudinal study of the vaginal, gut, and oral microbiomes in 292 samples from 10 subjects sampled every three weeks throughout pregnancy. Nonhuman sequences in the amount of 1.53 Gb were assembled into scaffolds, and functional genes were predicted for gene- and pathway-based analyses. Vaginal assemblies were binned into 97 draft quality genomes. Redundancy analysis (RDA) of microbial community composition at all three body sites revealed gestational age to be a significant source of variation in patterns of gene abundance. In addition, health complications were associated with variation in community functional gene composition in the mouth and gut. The diversity of Lactobacillus iners-dominated communities in the vagina, unlike most other vaginal community types, significantly increased with gestational age. The genomes of co-occurring Gardnerella vaginalis strains with predicted distinct functions were recovered in samples from two subjects. In seven subjects, gut samples contained strains of the same Lactobacillus species that dominated the vaginal community of that same subject and not other Lactobacillus species; however, these within-host strains were divergent. CRISPR spacer analysis suggested shared phage and plasmid populations across body sites and individuals. This work underscores the dynamic behavior of the microbiome during pregnancy and suggests the potential importance of understanding the sources of this behavior for fetal development and gestational outcome.

Insights into the Human Virome Using CRISPR Spacers from Microbiomes

Due to recent advances in next-generation sequencing over the past decade, our understanding of the human microbiome and its relationship to health and disease has increased dramatically. Yet, our insights into the human virome, and its interplay with important microbes that impact human health, is relatively limited. Prokaryotic and eukaryotic viruses are present throughout the human body, comprising a large and diverse population which influences several niches and impacts our health at various body sites. The presence of prokaryotic viruses like phages, has been documented at many different body sites, with the human gut being the richest ecological niche. Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) and associated proteins constitute the adaptive immune system of bacteria, which prevents attack by invasive nucleic acid. CRISPR-Cas systems function by uptake and integration of foreign genetic element sequences into the CRISPR array, which constitutes a genomic archive of iterative vaccination events. Consequently, CRISPR spacers can be investigated to reconstruct interplay between viruses and bacteria, and metagenomic sequencing data can be exploited to provide insights into host-phage interactions within a niche. Here, we show how the CRISPR spacer content of commensal and pathogenic bacteria can be used to determine the evidence of their phage exposure. This framework opens new opportunities for investigating host-virus dynamics in metagenomic data, and highlights the need to dedicate more efforts for virome sampling and sequencing.

Murine colitis reveals a disease-associated bacteriophage community

The dysregulation of intestinal microbial communities is associated with inflammatory bowel diseases (IBD). Studies aimed at understanding the contribution of the microbiota to inflammatory diseases have primarily focused on bacteria, yet the intestine harbors a viral component dominated by prokaryotic viruses known as bacteriophages (phages). Phage numbers are elevated at the intestinal mucosal surface and phages increase in abundance during IBD, suggesting that phages play an unidentified role in IBD. We used a sequence independent approach for the selection of viral contigs and then applied quantitative metagenomics to study intestinal phages in a mouse model of colitis. We discovered that during colitis the intestinal phage population is altered and transitions from an ordered state to a stochastic dysbiosis. We identified phages specific to pathobiotic hosts associated with intestinal disease, whose abundances are significantly altered during colitis. Additionally, phage populations in healthy and diseased mice overlapped with phages from healthy humans and humans with IBD. Our findings indicate that intestinal phage communities are altered during inflammatory disease establishing a platform for investigating phage involvement in IBD.

Structured Populations of Sulfolobus acidocaldarius with Susceptibility to Mobile Genetic Elements

The impact of a structured environment on genome evolution can be determined through comparative population genomics of species that live in the same habitat. Recent work comparing three genome sequences of Sulfolobus acidocaldarius suggested that highly structured, extreme, hot spring environments do not limit dispersal of this thermoacidophile, in contrast to other co-occurring Sulfolobus species. Instead, a high level of conservation among these three S. acidocaldarius genomes was hypothesized to result from rapid, global-scale dispersal promoted by low susceptibility to viruses that sets S. acidocaldarius apart from its sister Sulfolobus species. To test this hypothesis, we conducted a comparative analysis of 47 genomes of S. acidocaldarius from spatial and temporal sampling of two hot springs in Yellowstone National Park. While we confirm the low diversity in the core genome, we observe differentiation among S. acidocaldarius populations, likely resulting from low migration among hot spring “islands” in Yellowstone National Park. Patterns of genomic variation indicate that differing geological contexts result in the elimination or preservation of diversity among differentiated populations. We observe multiple deletions associated with a large genomic island rich in glycosyltransferases, differential integrations of the Sulfolobus turreted icosahedral virus, as well as two different plasmid elements. These data demonstrate that neither rapid dispersal nor lack of mobile genetic elements result in low diversity in the S. acidocaldarius genomes. We suggest instead that significant differences in the recent evolutionary history, or the intrinsic evolutionary rates, of sister Sulfolobus species result in the relatively low diversity of the S. acidocaldarius genome.

Temporal dynamics of uncultured viruses: a new dimension in viral diversity

Recent work has vastly expanded the known viral genomic sequence space, but the seasonal dynamics of viral populations at the genome level remain unexplored. Here we followed the viral community in a freshwater lake for 1 year using genome-resolved viral metagenomics, combined with detailed analyses of the viral community structure, associated bacterial populations and environmental variables. We reconstructed 8950 complete and partial viral genomes, the majority of which were not persistent in the lake throughout the year, but instead continuously succeeded each other. Temporal analysis of 732 viral genus-level clusters demonstrated that one-fifth were undetectable at specific periods of the year. Based on host predictions for a subset of reconstructed viral genomes, we for the first time reveal three distinct patterns of host-pathogen dynamics, where the viruses may peak before, during or after the peak in their host's abundance, providing new possibilities for modelling of their interactions. Time series metagenomics opens up a new dimension in viral profiling, which is essential to understand the full scale of viral diversity and evolution, and the ecological roles of these important factors in the global ecosystem.

Population-Genomic Insights into Variation in Prevotella intermedia and Prevotella nigrescens Isolates and Its Association with Periodontal Disease

High-throughput sequencing has helped to reveal the close relationship between Prevotella and periodontal disease, but the roles of subspecies diversity and genomic variation within this genus in periodontal diseases still need to be investigated. We performed a comparative genome analysis of 48 Prevotella intermedia and Prevotella nigrescens isolates that from the same cohort of subjects to identify the main drivers of their pathogenicity and adaptation to different environments. The comparisons were done between two species and between disease and health based on pooled sequences. The results showed that both P. intermedia and P. nigrescens have highly dynamic genomes and can take up various exogenous factors through horizontal gene transfer. The major differences between disease-derived and health-derived samples of P. intermedia and P. nigrescens were factors related to genome modification and recombination, indicating that the Prevotella isolates from disease sites may be more capable of genomic reconstruction. We also identified genetic elements specific to each sample, and found that disease groups had more unique virulence factors related to capsule and lipopolysaccharide synthesis, secretion systems, proteinases, and toxins, suggesting that strains from disease sites may have more specific virulence, particularly for P. intermedia. The differentially represented pathways between samples from disease and health were related to energy metabolism, carbohydrate and lipid metabolism, and amino acid metabolism, consistent with data from the whole subgingival microbiome in periodontal disease and health. Disease-derived samples had gained or lost several metabolic genes compared to healthy-derived samples, which could be linked with the difference in virulence performance between diseased and healthy sample groups. Our findings suggest that P. intermedia and P. nigrescens may serve as “crucial substances” in subgingival plaque, which may reflect changes in microbial and environmental dynamics in subgingival microbial ecosystems. This provides insight into the potential of P. intermedia and P. nigrescens as new predictive biomarkers and targets for effective interventions in periodontal disease.

MetaCRAST: reference-guided extraction of CRISPR spacers from unassembled metagenomes

Clustered regularly interspaced short palindromic repeat (CRISPR) systems are the adaptive immune systems of bacteria and archaea against viral infection. While CRISPRs have been exploited as a tool for genetic engineering, their spacer sequences can also provide valuable insights into microbial ecology by linking environmental viruses to their microbial hosts. Despite this importance, metagenomic CRISPR detection remains a major challenge. Here we present a reference-guided CRISPR spacer detection tool (Metagenomic CRISPR Reference-Aided Search Tool-MetaCRAST) that constrains searches based on user-specified direct repeats (DRs). These DRs could be expected from assembly or taxonomic profiles of metagenomes. We compared the performance of MetaCRAST to those of two existing metagenomic CRISPR detection tools-Crass and MinCED-using both real and simulated acid mine drainage (AMD) and enhanced biological phosphorus removal (EBPR) metagenomes. Our evaluation shows MetaCRAST improves CRISPR spacer detection in real metagenomes compared to the de novo CRISPR detection methods Crass and MinCED. Evaluation on simulated metagenomes show it performs better than de novo tools for Illumina metagenomes and comparably for 454 metagenomes. It also has comparable performance dependence on read length and community composition, run time, and accuracy to these tools. MetaCRAST is implemented in Perl, parallelizable through the Many Core Engine (MCE), and takes metagenomic sequence reads and direct repeat queries (FASTA or FASTQ) as input. It is freely available for download at https://github.com/molleraj/MetaCRAST.