Search

Ticks are blood-feeding ectoparasites that also transmit various pathogens, posing severe risks to human and animal health. DNA viruses play a crucial role in the microbial ecology of ticks, but their distribution and ecological significance remain largely undetermined. Here, we assembled an extensive catalog encompassing 4320 viral operational taxonomic units (vOTUs) from six main dominant tick species in China, of which 94.8% have not been found in any other environment. To bridge the knowledge gap in tick DNA virus research and provide a crucial resource platform, we developed the Tick DNA Virus Database. This database includes the vOTUs that are known to cause diseases. Most of the predicted vOTUs are associated with dominant bacterial and archaeal phyla. We identified 105 virus-encoded putative auxiliary metabolic genes (AMGs) that are involved in host metabolism and environmental adaptation, potentially influencing ticks through both top-down and bottom-up mechanisms. The identification of microbial communities and antibiotic resistance in wild tick species suggests that wild ticks are reservoirs of antibiotic resistance and potential spreaders of antibiotic resistance. These findings reveal the potential role of tick viruses in ecosystems, highlighting the importance of monitoring tick microbiomes to address global public health challenges. [ABSTRACT FROM AUTHOR]

Purpose: Viruses are the largest genetic repository and most diverse host-associated replicating organisms in a wide range of ecosystems in the entire biosphere. Although playing crucial roles in numerous processes such as nutrient cycles, viral shunt, and population control within marine ecosystems, viral communities continue to be the least explored biological entities, particularly in extreme marine habitats. Therefore, details about the viral diversity and their metabolic potential from Arctic sediments has great importance. Methods: Sediment samples were collected from two Arctic fjords, Kongsfjorden and Krossfjorden, in Svalbard. Viral nucleic acids were extracted from both fjord sediments after a pre-processing step. Illumina Novaseq based shotgun sequencing was performed and bioinformatic analysis was done to explore viral communities and their functional potential. Results: An array of DNA and RNA viruses especially bacteriophages along with viruses infecting algae, plants, molluscs, shrimps, fishes, amphibians, reptiles, birds, humans, and other mammals were detected from both fjords. The double-stranded DNA viruses were the most abundant (Siphoviridae, Podoviridae, Myoviridae, and Phycodnaviridae), followed by single-stranded DNA (Inoviridae) and RNA (Retroviridae and Betaflexiviridae) viruses from both metavirome. Functional analysis explored genes encoding virus structures, enzymes for phage replication, integration and excision and proteins related to phage regulation of gene expression. Conclusions: Fjord metavirome analysis revealed the occurrence of virus groups endemic to Arctic freshwater and marine habitats along with an extra pool of unclassified or unassigned virus reads. Viruses infecting a variety of bacterial groups, and other higher trophic levels in fjord environments were explored. Functional annotation revealed the abundance of phage-related structural genes and metabolic genes. Comparative analysis revealed the abundance of dsDNA viruses from Caudovirales (Myoviridae, Podoviridae and Siphoviridae) in all the analysed Arctic samples in different proportions. The present study enhances our understanding of viral diversity and their metabolic potential in Arctic fjord sediments. [ABSTRACT FROM AUTHOR]

ABSTRACT The surface of smear-ripened cheeses constitutes a dynamic microbial ecosystem resulting from the successive development of different microbial groups such as lactic acid bacteria, fungi, and ripening bacteria. Recent studies indicate that a viral community, mainly composed of bacteriophages, also represents a common and substantial part of the cheese microbiome. However, the composition of this community, its temporal variations, and associations between bacteriophages and their hosts remain poorly characterized. Here, we studied a French smear-ripened cheese by both viral metagenomics and 16S metabarcoding approaches to assess both the succession of phages and bacterial communities on the cheese surface during cheese ripening and their temporal variations in ready-to-eat cheeses over the years of production. We observed a clear transition of the phage community structure during ripening with a decreased relative abundance of viral species (vOTUs) associated with Lactococcus phages, which were replaced by vOTUs associated with phages infecting ripening bacteria such as Brevibacterium, Glutamicibacter, Pseudoalteromonas, and Vibrio. The dynamics of the phage community was strongly associated with bacterial successions observed on the cheese surface. Finally, while some variations in the distribution of phages were observed in ready-to-eat cheeses produced at different dates spanning more than 4 years of production, the most abundant phages were detected throughout. This result revealed the long-term persistence of the dominant phages in the cheese production environment. Together, these findings offer novel perspectives on the ecology of bacteriophages in smear-ripened cheese and emphasize the significance of incorporating bacteriophages in the microbial ecology studies of fermented foods.IMPORTANCEThe succession of diverse microbial populations is critical for ensuring the production of high-quality cheese. We observed a temporal succession of phages on the surface of a smear-ripened cheese, with new phage communities showing up when ripening bacteria start covering this surface. Interestingly, the final phage community of this cheese is also consistent over large periods of time, as the same bacteriophages were found in cheese products from the same manufacturer made over 4 years. This research highlights the importance of considering these bacteriophages when studying the microbial life of fermented foods like cheese.

Ticks are blood-feeding ectoparasites that also transmit various pathogens, posing severe risks to human and animal health. DNA viruses play a crucial role in the microbial ecology of ticks, but their distribution and ecological significance remain largely undetermined. Here, we assembled an extensive catalog encompassing 4320 viral operational taxonomic units (vOTUs) from six main dominant tick species in China, of which 94.8% have not been found in any other environment. To bridge the knowledge gap in tick DNA virus research and provide a crucial resource platform, we developed the Tick DNA Virus Database. This database includes the vOTUs that are known to cause diseases. Most of the predicted vOTUs are associated with dominant bacterial and archaeal phyla. We identified 105 virus-encoded putative auxiliary metabolic genes (AMGs) that are involved in host metabolism and environmental adaptation, potentially influencing ticks through both top-down and bottom-up mechanisms. The identification of microbial communities and antibiotic resistance in wild tick species suggests that wild ticks are reservoirs of antibiotic resistance and potential spreaders of antibiotic resistance. These findings reveal the potential role of tick viruses in ecosystems, highlighting the importance of monitoring tick microbiomes to address global public health challenges.

Abstract Background There is an increasing interest in investigating the human gut virome for its influence on the gut bacterial community and its putative influence on the trajectory towards health or disease. Most gut virome studies are based on sequencing of stored fecal samples. However, relatively little is known about how conventional storage buffers and storage conditions affect the infectivity of bacteriophages and influence the downstream metavirome sequencing. Results We demonstrate that the infectivity and genome recovery rate of different spiked bacteriophages (T4, c2 and Phi X174) are variable and highly dependent on storage buffers. Regardless of the storage temperature and timespan, all tested phages immediately lost 100% (DNA/RNA Shield) or more than 90% (StayRNA and RNAlater) of their infectivity. Generally, in SM buffer at 4 °C phage infectivity was preserved for up to 30 days and phage DNA integrity was maintained for up to 100 days. While in CANVAX, the most effective buffer, all spiked phage genomes were preserved for at least 100 days. Prolonged storage time (500 days) at – 80 °C impacted viral diversity differently in the different buffers. Samples stored in CANVAX or DNA/RNA Shield buffer had the least shifts in metavirome composition, after prolonged storage, but they yielded more contigs classified as “uncharacterised”. Moreover, in contrast to the SM buffer, these storage buffers yielded a higher fraction of bacterial DNA in metavirome-sequencing libraries. We demonstrated that the latter was due to inactivation of the DNases employed to remove extra-cellular DNA during virome extraction. The latter could be partly avoided by employing additional washing steps prior to virome extraction. Conclusion Fecal sample storage buffers and storage conditions (time and temperature) strongly influence bacteriophage infectivity and viral composition as determined by plaque assay and metavirome sequencing. The choice of buffer had a larger effect than storage temperature and storage time on the quality of the viral sequences and analyses. Based on these results, we recommend storage of fecal virome samples at in SM buffer at 4 °C for the isolation of viruses and at – 80 °C for metagenomic applications if practically feasible (i.e., access to cold storage). For fecal samples stored in other buffers, samples should be cleared of these buffers before viral extraction and sequencing. Video Abstract

Background: There is an increasing interest in investigating the human gut virome for its influence on the gut bacterial community and its putative influence on the trajectory towards health or disease. Most gut virome studies are based on sequencing of stored fecal samples. However, relatively little is known about how conventional storage buffers and storage conditions affect the infectivity of bacteriophages and influence the downstream metavirome sequencing. Results: We demonstrate that the infectivity and genome recovery rate of different spiked bacteriophages (T4, c2 and Phi X174) are variable and highly dependent on storage buffers. Regardless of the storage temperature and timespan, all tested phages immediately lost 100% (DNA/RNA Shield) or more than 90% (StayRNA and RNAlater) of their infectivity. Generally, in SM buffer at 4 °C phage infectivity was preserved for up to 30 days and phage DNA integrity was maintained for up to 100 days. While in CANVAX, the most effective buffer, all spiked phage genomes were preserved for at least 100 days. Prolonged storage time (500 days) at – 80 °C impacted viral diversity differently in the different buffers. Samples stored in CANVAX or DNA/RNA Shield buffer had the least shifts in metavirome composition, after prolonged storage, but they yielded more contigs classified as "uncharacterised". Moreover, in contrast to the SM buffer, these storage buffers yielded a higher fraction of bacterial DNA in metavirome-sequencing libraries. We demonstrated that the latter was due to inactivation of the DNases employed to remove extra-cellular DNA during virome extraction. The latter could be partly avoided by employing additional washing steps prior to virome extraction. Conclusion: Fecal sample storage buffers and storage conditions (time and temperature) strongly influence bacteriophage infectivity and viral composition as determined by plaque assay and metavirome sequencing. The choice of buffer had a larger effect than storage temperature and storage time on the quality of the viral sequences and analyses. Based on these results, we recommend storage of fecal virome samples at in SM buffer at 4 °C for the isolation of viruses and at – 80 °C for metagenomic applications if practically feasible (i.e., access to cold storage). For fecal samples stored in other buffers, samples should be cleared of these buffers before viral extraction and sequencing. 2F_FbCzVMdkP5vWvM-fCta Video Abstract [ABSTRACT FROM AUTHOR]

A wide diversity of pathogenic mosquito-borne viruses circulate in the Brazilian Amazon, and the intense deforestation can contribute to the spread of these viruses. In this context, this study aimed to investigate the viral diversity in mosquitoes of the genera Aedes, Culex, Haemagogus, and Sabethes from a transition area between the Amazon, Cerrado, and Caatinga biomes in Brazil. Metagenomic high-throughput sequencing was used to characterize the virome of 20 mosquito pools. A total of 15 virus-like genomes were identified, comprising species genomically close to insect-specific viruses of the families Iflaviridae, Metaviridae, Lispiviridae, Rhabdoviridae, Xinmoviridae, and Parvoviridae and species of plant viruses of the families Solemoviridae, Virgaviridae, and Partitiviridae. However, sequences of viruses associated with human and animal diseases were not detected. Most of the recovered genomes were divergent from those previously described. These findings reveal that there are a large number of unknown viruses to be explored in the middle-north of Brazil. [ABSTRACT FROM AUTHOR]

Despite their abundance and ecological importance, little is known about the diversity of marine viruses, in part because most cannot be cultured in the laboratory. Here, we used high-throughput viral metagenomics of uncultivated viruses to investigate the dynamics of DNA viruses in tropical seawater sampled from Chuuk State, Federated States of Micronesia, in March, June, and December 2014. Among the identified viruses, 71–79% were bacteriophages belonging to the families Myoviridae, Siphoviridae, and Podoviridae (Caudoviriales), listed in order of abundance at all sampling times. Although the measured environmental factors (temperature, salinity, and pH) remained unchanged in the seawater over time, viral dynamics changed. The proportion of cyanophages (34.7%) was highest in June, whereas the proportion of mimiviruses, phycodnaviruses, and other nucleo-cytoplasmic large DNA viruses (NCLDVs) was higher in March and December. Although host species were not analysed, the dramatic viral community change observed in June was likely due to changes in the abundance of cyanophage-infected cyanobacteria, whereas that in NCLDVs was likely due to the abundance of potential eukaryote-infected hosts. These results serve as a basis for comparative analyses of other marine viral communities, and guide policy-making when considering marine life care in Chuuk State. [ABSTRACT FROM AUTHOR]

Ticks are the second most harmful infectious agent in the world after mosquitoes and can transmit a variety of pathogens. The surveillance of pathogen in ticks is one of the important means for the prevention and control of tick-borne diseases and tick-borne pathogens. This paper reviewed the literature published in the past two decades on the principles， advantages， disadvantages and applications of current tick pathogen detection technologies. The paper introduced the development and application of various PCR technologies， and looks forward to the future direction. Finally， it points out that metagenomics represented by high-throughput sequencing will become the mainstream development direction of pathogen detection technology in the future.

Background: Metagenomics is a promising approach to discovering novel sequences of microorganisms in environmental samples. A recently published RNA-Seq data of Solanum tuberosum from China was used for a metavirome study. Methods: RNA-seq data of a BioSample project of S. tuberosum containing sequence read archive (SRA) of six plant samples were imported into the Galaxy server. Transcriptome data were de novo assembled for viral sequences with rnaviralSPAdes. The contig files were further organized by VirHunter and Kraken tools. Raw SRA data were trimmed and assembled for severe acute respiratory syndrome coronavirus 2 (SARS‑CoV‑2) genome by coronaSPAdes. The scaffolds were arranged by pairwise alignment against the SARS-CoV-2 reference genome (NC_045512.2). Coronavirus Typing Tool, Nextstrain, and Pangolin platforms were used to further investigate the SARS-CoV-2 genotype, phylogenetic analysis, and mutation estimations. Results: Several environmentally related non-intact virus sequence reads from forest animals, moths, bacteria, and amoeba were detected. Further investigation resulted in non-indigenous sequences of SARS-CoV-2 genomes of lineage B with novel substitutions. Three polymorphisms, including A22D and A36V in the envelope protein, and Q498H in the spike (S) glycoprotein that were recently reported from a mice-adopted strain of SARS-CoV-2 with enhanced virulence were detected in all samples. Further novel substitutions at ORF1ab were also uncovered. These were L1457V, D4553N, W6538S, I1525T, D1585Y, D6928G, N3414K, and T3432S. Two unexpected frameshifts, ORF1a:2338-4401 and ORF1a:3681-4401, were also detected. Conclusion: The findings of the presented study highlight the threats of the emerged potentially severe genotypes of SARS-CoV-2 bearing substitutions that are not yet clinically reported.

Research suggests a potential role of the oral-neuro and gut-brain axes in schizophrenia, involving non-brain microbiomes such as salivary and gut microbiomes. However, the blood-brain barrier effectively prevents microorganism entry. Additionally, despite approximately 8% of the human genome consisting of retroviruses and the established link between viral infections and schizophrenia, the presence of a resident virome (a viral component of the microbiome) in the brain and its association with mental disorders remain unexplored. Methods: Whole-genome sequencing raw data from postmortem Brodmann Area 46 (BA46) tissue from 49 individuals (20 healthy controls [HCs], 29 with schizophrenia [SCZs]) obtained from the NCBI SRA database from BioProject: PRJNA422380.Virome profiles were retrieved using Metaphlan3, and viral signatures were identified using linear discriminant analysis effect size (LEfSe). Mann-Whitney tests and receiver operating characteristic curve validated the viral signatures. Results: In BA46, 30 distinct species representing 9 phyla, 10 classes, 10 orders, 13 families, and 19 genera were identified. HCs exhibited greater alpha diversity, and there were significant differences in beta diversity between the groups. LEfSe analysis highlighted distinct viral levels, including Escherichia virus Lambda , Escherichia virus phiV10 , Human endogenous retrovirus K , Taterapox virus , Alcelaphine gammaherpesvirus 1 , and Bovine gammaherpesvirus 4 in HCs, while Glypta fumiferanae ichnovirus and unknown virus showed higher levels in schizophrenia. Conclusion: This is the first study to identify a human brain virome associated with schizophrenia in BA46. Brain virome dysbiosis may be associated with mental illness, and viral signatures may serve as biomarkers for the early detection of schizophrenia. [ABSTRACT FROM AUTHOR]

Background: Metagenomics is a promising approach to discovering novel sequences of microorganisms in environmental samples. A recently published RNA-Seq data of Solanum tuberosum from China was used for a metavirome study. Methods: RNA-seq data of a BioSample project of S. tuberosum containing sequence read archive (SRA) of six plant samples were imported into the Galaxy server. Transcriptome data were de novo assembled for viral sequences with rnaviralSPAdes. The contig files were further organized by VirHunter and Kraken tools. Raw SRA data were trimmed and assembled for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome by coronaSPAdes. The scaffolds were arranged by pairwise alignment against the SARS-CoV-2 reference genome (NC_045512.2). Coronavirus Typing Tool, Nextstrain, and Pangolin platforms were used to further investigate the SARS-CoV-2 genotype, phylogenetic analysis, and mutation estimations. Results: Several environmentally related non-intact virus sequences reads from forest animals, moths, bacteria, and amoeba were detected. Further investigation resulted in non-indigenous sequences of SARS-CoV-2 genomes of lineage B with novel substitutions. Three polymorphisms, including A22D and A36V in the envelope protein, and Q498H in the spike (S) glycoprotein that were recently reported from a mice-adopted strain of SARS-CoV-2 with enhanced virulence were detected in all samples. Further novel substitutions at ORF1ab were also uncovered. These were L1457V, D4553N, W6538S, I1525T, D1585Y, D6928G, N3414K, and T3432S. Two unexpected frameshifts, ORF1a:2338-4401 and ORF1a:3681-4401, were also detected. Conclusion: The findings of the presented study highlight the threats of the emerged potentially severe genotypes of SARS-CoV-2 bearing substitutions that are not yet clinically reported. [ABSTRACT FROM AUTHOR]

ABSTRACT Microbial starter cultures are used in the production of many cheeses around the world, such as Parmigiano-Reggiano, in Italy, Époisses, in France, and Canastra, in Brazil, providing many of the unique features of these cheeses. Bacteriophages (phages) are ubiquitous and well known to modulate the structure of bacterial communities, and recent data indicate that cheeses contain a high abundance of naturally occurring phages. Here, we analyze the viral and bacterial metagenomes of Canastra cheese: a traditional artisanal Brazilian cheese produced using an endogenous starter culture and raw milk. Over 1,200 viral operational taxonomic units were recovered using both isolated viral-like particles and complete metagenomic DNA. Common viral families identified included Siphoviridae and Myoviridae, with 40% of putative phage genomes unidentified at the family level of classification. We observed very high phage diversity, which varied greatly across different cheese producers, with 28% of phage genomes detected in only one producer. Several metagenome-assembled genomes were recovered for lactic acid-producing bacteria, as well as nonstarter bacterial species, and we identified several phage-bacterium interactions, at the strain level of resolution, varying across distinct cheese producers. We postulate that at least one bacterial strain detected could be endogenous and unique to the Canastra cheese-producing region in Brazil and that its growth seems to be modulated by autochthonous phages present in this artisanal production system. This phage-host relationship is likely to influence the fermentation dynamics and ultimately the sensorial profile of these cheeses, with implications for other similar cheese production systems around the world. IMPORTANCE Our work demonstrated a dynamic yet stable microbial ecosystem during cheese production using an endogenous starter culture. This was observed across several distinct producers and was marked by genomic evidence of continued phage-bacterium interactions, such as the presence of bacterial defense mechanisms. Furthermore, we provide evidence of unique microbial signatures for each individual cheese producer studied in the region, a fact that may have profound consequences on product traceability. This was the first effort to describe and understand the bacteriophage composition and ecological dynamics within the Brazilian Canastra cheese production system. The study of this prototypical backslopping production system provides a solid background for further mechanistic studies of the production of many cheeses around the world.

Viruses are key players in the environment, and recent metagenomic studies have revealed their diversity and genetic complexity. Despite progress in understanding the ecology of viruses in extreme environments, viruses' dynamics and functional roles in dryland ecosystems, which cover about 45% of the Earth's land surfaces, remain largely unexplored. This study characterizes virus sequences in the metagenomes of endolithic (within rock) microbial communities ubiquitously found in hyper-arid deserts. Taxonomic classification and network construction revealed the presence of novel and diverse viruses in communities inhabiting calcite, gypsum, and ignimbrite rocks. Viral genome maps show a high level of protein diversity within and across endolithic communities and the presence of virus-encoded auxiliary metabolic genes. Phage-host relationships were predicted by matching tRNA, CRISPR spacer, and protein sequences in the viral and microbial metagenomes. Primary producers and heterotrophic bacteria were found to be putative hosts to some viruses. Intriguingly, viral diversity was not correlated with microbial diversity across rock substrates. [ABSTRACT FROM AUTHOR]

Abstract Background The full biosphere structure and functional exploration of the microbial communities of the Challenger Deep of the Mariana Trench, the deepest known hadal zone on Earth, lag far behind that of other marine realms. Results We adopt a deep metagenomics approach to investigate the microbiome in the sediment of Challenger Deep, Mariana Trench. We construct 178 metagenome-assembled genomes (MAGs) representing 26 phyla, 16 of which are reported from hadal sediment for the first time. Based on the MAGs, we find the microbial community functions are marked by enrichment and prevalence of mixotrophy and facultative anaerobic metabolism. The microeukaryotic community is found to be dominated by six fungal groups that are characterized for the first time in hadal sediment to possess the assimilatory and dissimilatory nitrate/sulfate reduction, and hydrogen sulfide oxidation pathways. By metaviromic analysis, we reveal novel hadal Caudovirales clades, distinctive virus-host interactions, and specialized auxiliary metabolic genes for modulating hosts’ nitrogen/sulfur metabolism. The hadal microbiome is further investigated by large-scale cultivation that cataloged 1070 bacterial and 19 fungal isolates from the Challenger Deep sediment, many of which are found to be new species specialized in the hadal habitat. Conclusion Our hadal MAGs and isolates increase the diversity of the Challenger Deep sediment microbial genomes and isolates present in the public. The deep metagenomics approach fills the knowledge gaps in structure and diversity of the hadal microbiome, and provides novel insight into the ecology and metabolism of eukaryotic and viral components in the deepest biosphere on earth.

Despite their abundance and ecological importance, little is known about the diversity of marine viruses, in part because most cannot be cultured in the laboratory. Here, we used high-throughput viral metagenomics of uncultivated viruses to investigate the dynamics of DNA viruses in tropical seawater sampled from Chuuk State, Federated States of Micronesia, in March, June, and December 2014. Among the identified viruses, 71–79% were bacteriophages belonging to the families Myoviridae, Siphoviridae, and Podoviridae (Caudoviriales), listed in order of abundance at all sampling times. Although the measured environmental factors (temperature, salinity, and pH) remained unchanged in the seawater over time, viral dynamics changed. The proportion of cyanophages (34.7%) was highest in June, whereas the proportion of mimiviruses, phycodnaviruses, and other nucleo-cytoplasmic large DNA viruses (NCLDVs) was higher in March and December. Although host species were not analysed, the dramatic viral community change observed in June was likely due to changes in the abundance of cyanophage-infected cyanobacteria, whereas that in NCLDVs was likely due to the abundance of potential eukaryote-infected hosts. These results serve as a basis for comparative analyses of other marine viral communities, and guide policy-making when considering marine life care in Chuuk State.

Viruses are the most prodigious repertory of the genetic material on the earth. They are elusive, breakneck, evolutionary life particles that constitute a riveting concealed world. Environmental viruses have been obscurely explored, and hence, such an intriguing world of viruses was studied in the Himalayan Geothermal Belt of Indian peninsula at Sikkim corridor through hot springs. The hot springs located at the North Sikkim district were selected for the current study. The solfataric mud sediment samples were pooled from both the hot springs. The virus community showed significant diversity among the two hot springs of Yume Samdung. Reads for viruses among the mud sediments at Old Yume Samdung hot springs (OYS) was observed to be 11% and in the case of New Yume Samdung hot springs (NYS) it was 6%. Both the hot springs were abundant in dsDNA viromes. The metavirome reads in both the OYS and NYS hot spring mud sediments showed the predominance of Caudovirales; Herpesvirales; Ortervirales among which viral reads from Siphoviridae, Myoviridae, Phycodnaviridae and Podoviridae were abundantly present. Other viral communities belonged to families like Baculoviridae, Mimiviridae, Parvoviridae, Marseilleviridae etc. Interestingly, in the case of NYS, the unassigned group reads belonged to some unclassified giant DNA viruses like genera Pandoravirus and Pithovirus. Other interesting findings were – reads for Badnavirus having ds (RT-DNA) was exclusively found in NYS whereas Rubulavirus having ss(-)RNA was exclusively found in OYS sample. This is the first ever report on viruses from any hot springs of Sikkim till date.

The surface of smear-ripened cheeses constitutes a dynamic microbial ecosystem resulting from the successive development of different microbial groups such as lactic acid bacteria, fungi, and ripening bacteria. Recent studies indicate that a viral community, mainly composed of bacteriophages, also represents a common and substantial part of the cheese microbiome. However, the composition of this community, its temporal variations, and associations between bacteriophages and their hosts remain poorly characterized. Here, we studied a French smear-ripened cheese by both viral metagenomics and 16S metabarcoding approaches to assess both the succession of phages and bacterial communities on the cheese surface during cheese ripening and their temporal variations in ready-to-eat cheeses over the years of production. We observed a clear transition of the phage community structure during ripening with a decreased relative abundance of viral species (vOTUs) associated with Lactococcus phages, which were replaced by vOTUs associated with phages infecting ripening bacteria such as Brevibacterium, Glutamicibacter, Pseudoalteromonas, and Vibrio . The dynamics of the phage community was strongly associated with bacterial successions observed on the cheese surface. Finally, while some variations in the distribution of phages were observed in ready-to-eat cheeses produced at different dates spanning more than 4 years of production, the most abundant phages were detected throughout. This result revealed the long-term persistence of the dominant phages in the cheese production environment. Together, these findings offer novel perspectives on the ecology of bacteriophages in smear-ripened cheese and emphasize the significance of incorporating bacteriophages in the microbial ecology studies of fermented foods.IMPORTANCEThe succession of diverse microbial populations is critical for ensuring the production of high-quality cheese. We observed a temporal succession of phages on the surface of a smear-ripened cheese, with new phage communities showing up when ripening bacteria start covering this surface. Interestingly, the final phage community of this cheese is also consistent over large periods of time, as the same bacteriophages were found in cheese products from the same manufacturer made over 4 years. This research highlights the importance of considering these bacteriophages when studying the microbial life of fermented foods like cheese., Competing Interests: The authors declare no conflict of interest.

Viruses are key players in the environment, and recent metagenomic studies have revealed their diversity and genetic complexity. Despite progress in understanding the ecology of viruses in extreme environments, viruses’ dynamics and functional roles in dryland ecosystems, which cover about 45% of the Earth’s land surfaces, remain largely unexplored. This study characterizes virus sequences in the metagenomes of endolithic (within rock) microbial communities ubiquitously found in hyper-arid deserts. Taxonomic classification and network construction revealed the presence of novel and diverse viruses in communities inhabiting calcite, gypsum, and ignimbrite rocks. Viral genome maps show a high level of protein diversity within and across endolithic communities and the presence of virus-encoded auxiliary metabolic genes. Phage-host relationships were predicted by matching tRNA, CRISPR spacer, and protein sequences in the viral and microbial metagenomes. Primary producers and heterotrophic bacteria were found to be putative hosts to some viruses. Intriguingly, viral diversity was not correlated with microbial diversity across rock substrates.

Twenty complete genomes (29-63 kb) and 29 genomes with an estimated completeness of over 90 % (30-90 kb) were identified for novel dsDNA viruses in the Yangshan Harbor metavirome. These newly discovered viruses contribute to the expansion of viral taxonomy by introducing 46 potential new families. Except for one virus, all others belong to the class Caudoviricetes . The exception is a novel member of the recently characterized viral group known as Gossevirus. Fifteen viruses were predicted to be temperate. The predicted hosts for the viruses appear to be involved in various aspects of the nitrogen cycle, including nitrogen fixation, oxidation and denitrification. Two viruses were identified to have a host of Flavobacterium and Tepidimonas fonticaldi , respectively, by matching CRISPR spacers with viral protospacers. Our findings provide an overview for characterizing and identifying specific viruses from Yangshan Harbor. The Gossevirus-like virus uncovered emphasizes the need for further comprehensive isolation and investigation of polinton-like viruses.

Sponges are remarkable holobionts harboring extremely diverse microbial and viral communities. However, the interactions between the components within holobionts and between a holobiont and environment are largely unknown, especially for polar organisms. To investigate possible interactions within and between sponge-associated communities, we probed the microbiomes and viromes of cold-water sympatric sponges Isodictya palmata (n = 2), Halichondria panicea (n = 3), and Halichondria sitiens (n = 3) by 16S and shotgun metagenomics. We showed that the bacterial and viral communities associated with these White Sea sponges are species-specific and different from the surrounding water. Extensive mining of bacterial antiphage defense systems in the metagenomes revealed a variety of defense mechanisms. The abundance of defense systems was comparable in the metagenomes of the sponges and the surrounding water, thus distinguishing the White Sea sponges from those inhabiting the tropical seas. We developed a network-based approach for the combined analysis of CRISPR-spacers and protospacers. Using this approach, we showed that the virus–host interactions within the sponge-associated community are typically more abundant (three out of four interactions studied) than the inter-community interactions. Additionally, we detected the occurrence of viral exchanges between the communities. Our work provides the first insight into the metagenomics of the three cold-water sponge species from the White Sea and paves the way for a comprehensive analysis of the interactions between microbial communities and associated viruses.

Metagenomic sequencing is a promising method to determine the virus diversity in environmental samples such as sewage or shellfish. However, to identify the short RNA genomes of human enteric viruses among the large diversity of nucleic acids present in such complex matrices, method optimization is still needed. This work presents methodological developments focused on norovirus, a small ssRNA non-enveloped virus known as the major cause of human gastroenteritis worldwide and frequently present in human excreta and sewage. Different elution protocols were applied and Illumina MiSeq technology were used to study norovirus diversity. A double approach, agnostic deep sequencing and a capture-based approach (VirCapSeq-VERT) was used to identify norovirus in environmental samples. Family-specific viral contigs were classified and sorted by SLIM and final norovirus contigs were genotyped using the online Norovirus genotyping tool v2.0. From sewage samples, 14 norovirus genogroup I sequences were identified of which six were complete genomes. For norovirus genogroup II, nine sequences were identified and three of them comprised more than half of the genome. In oyster samples bioaccumulated with these sewage samples, only the use of an enrichment step during library preparation allowed successful identification of nine different sequences of norovirus genogroup I and four for genogroup II (>500 bp). This study demonstrates the importance of method development to increase virus recovery, and the interest of a capture-based approach to be able to identify viruses present at low concentrations.

Metagenomic sequencing is a promising method to determine the virus diversity in environmental samples such as sewage or shellfish. However, to identify the short RNA genomes of human enteric viruses among the large diversity of nucleic acids present in such complex matrices, method optimization is still needed. This work presents methodological developments focused on norovirus, a small ssRNA non-enveloped virus known as the major cause of human gastroenteritis worldwide and frequently present in human excreta and sewage. Different elution protocols were applied and Illumina MiSeq technology were used to study norovirus diversity. A double approach, agnostic deep sequencing and a capture-based approach (VirCapSeq-VERT) was used to identify norovirus in environmental samples. Family-specific viral contigs were classified and sorted by SLIM and final norovirus contigs were genotyped using the online Norovirus genotyping tool v2.0. From sewage samples, 14 norovirus genogroup I sequences were identified of which six were complete genomes. For norovirus genogroup II, nine sequences were identified and three of them comprised more than half of the genome. In oyster samples bioaccumulated with these sewage samples, only the use of an enrichment step during library preparation allowed successful identification of nine different sequences of norovirus genogroup I and four for genogroup II (>500 bp). This study demonstrates the importance of method development to increase virus recovery, and the interest of a capture-based approach to be able to identify viruses present at low concentrations. [ABSTRACT FROM AUTHOR]

With the emergence of Next Generation Sequencing, major advances were made with regard to identifying viruses in natural environments. However, bioinformatical research on viruses is still limited because of the low amounts of viral DNA that can be obtained for analysis. To overcome this limitation, DNA is often amplified with multiple displacement amplification (MDA), which may cause an unavoidable bias. Here, we describe a case study in which the virome of a bioreactor is sequenced using Ion Torrent technology. DNA-spiking of samples is compared with MDA-amplified samples. DNA for spiking was obtained by amplifying a bacterial 16S rRNA gene. After sequencing, the 16S rRNA gene reads were removed by mapping to the Silva database. Three samples were tested, a whole genome from Enterobacteria P1 Phage and two viral metagenomes from an infected bioreactor. For one sample, the new DNA-spiking protocol was compared with the MDA technique. When MDA was applied, the overall GC content of the reads showed a bias towards lower GC%, indicating a change in composition of the DNA sample. Assemblies using all available reads from both MDA and the DNA-spiked samples resulted in six viral genomes. All six genomes could be almost completely retrieved (97.9%–100%) when mapping the reads from the DNA-spiked sample to those six genomes. In contrast, 6.3%–77.7% of three viral genomes was covered by reads obtained using the MDA amplification method and only three were nearly fully covered (97.4%–100%). This case study shows that DNA-spiking could be a simple and inexpensive alternative with very low bias for sequencing of metagenomes for which low amounts of DNA are available.

The Formosan subterranean termite; Coptotermes formosanus is nutritionally dependent on the complex and diverse community of bacteria and protozoa in their gut. Although, there have been many studies to decipher the taxonomic and functional diversity of bacterial communities in the guts of termites, their bacteriophages remain unstudied. We sequenced the metavirome of the guts of Formosan subterranean termite workers to study the diversity of bacteriophages and other associated viruses. Results showed that the termites harbor a virome in their gut comprised of varied and previously unknown bacteriophages. Between 87–90% of the predicted dsDNA virus genes by Metavir showed similarity to the tailed bacteriophages (Caudovirales). Many predicted genes from the virome matched to bacterial prophage regions. These data are suggestive of a virome dominated by temperate bacteriophages. We predicted the genomes of seven novel Caudovirales bacteriophages from the termite gut. Three of these predicted bacteriophage genomes were found in high proportions in all the three termite colonies tested. Two bacteriophages are predicted to infect endosymbiotic bacteria of the gut protozoa. The presence of these putative bacteriophages infecting endosymbionts of the gut protozoa, suggests a quadripartite relationship between the termites their symbiotic protozoa, endosymbiotic bacteria of the protozoa and their bacteriophages. Other than Caudovirales, ss-DNA virus related genes were also present in the termite gut. We predicted the genomes of 12 novel Microviridae phages from the termite gut and seven of those possibly represent a new proposed subfamily. Circovirus like genomes were also assembled from the termite gut at lower relative abundance. We predicted 10 novel circovirus genomes in this study. Whether these circoviruses infect the termites remains elusive at the moment. The functional and taxonomical annotations suggest that the termites may harbor a core virome comprised of the bacteriophages infecting endosymbionts of the gut protozoa.

With the emergence of Next Generation Sequencing, major advances were made with regard to identifying viruses in natural environments. However, bioinformatical research on viruses is still limited because of the low amounts of viral DNA that can be obtained for analysis. To overcome this limitation, DNA is often amplified with multiple displacement amplification (MDA), which may cause an unavoidable bias. Here, we describe a case study in which the virome of a bioreactor is sequenced using Ion Torrent technology. DNA-spiking of samples is compared with MDA-amplified samples. DNA for spiking was obtained by amplifying a bacterial 16S rRNA gene. After sequencing, the 16S rRNA gene reads were removed by mapping to the Silva database. Three samples were tested, a whole genome from Enterobacteria P1 Phage and two viral metagenomes from an infected bioreactor. For one sample, the new DNA-spiking protocol was compared with the MDA technique. When MDA was applied, the overall GC content of the reads showed a bias towards lower GC%, indicating a change in composition of the DNA sample. Assemblies using all available reads from both MDA and the DNA-spiked samples resulted in six viral genomes. All six genomes could be almost completely retrieved (97.9%-100%) when mapping the reads from the DNA- spiked sample to those six genomes. In contrast, 6.3%-77.7% of three viral genomes was covered by reads obtained using the MDA amplification method and only three were nearly fully covered (97.4%-100%). This case study shows that DNA-spiking could be a simple and inexpensive alternative with very low bias for sequencing of metagenomes for which low amounts of DNA are available. [ABSTRACT FROM AUTHOR]

The Formosan subterranean termite; Coptotermes formosanus is nutritionally dependent on the complex and diverse community of bacteria and protozoa in their gut. Although, there have been many studies to decipher the taxonomic and functional diversity of bacterial communities in the guts of termites, their bacteriophages remain unstudied. We sequenced the metavirome of the guts of Formosan subterranean termite workers to study the diversity of bacteriophages and other associated viruses. Results showed that the termites harbor a virome in their gut comprised of varied and previously unknown bacteriophages. Between 87-90% of the predicted dsDNA virus genes by Metavir showed similarity to the tailed bacteriophages (Caudovirales). Many predicted genes from the virome matched to bacterial prophage regions. These data are suggestive of a virome dominated by temperate bacteriophages. We predicted the genomes of seven novel Caudovirales bacteriophages from the termite gut. Three of these predicted bacteriophage genomes were found in high proportions in all the three termite colonies tested. Two bacteriophages are predicted to infect endosymbiotic bacteria of the gut protozoa. The presence of these putative bacteriophages infecting endosymbionts of the gut protozoa, suggests a quadripartite relationship between the termites their symbiotic protozoa, endosymbiotic bacteria of the protozoa and their bacteriophages. Other than Caudovirales, ss-DNA virus related genes were also present in the termite gut. We predicted the genomes of 12 novel Microviridae phages from the termite gut and seven of those possibly represent a new proposed subfamily. Circovirus like genomes were also assembled from the termite gut at lower relative abundance. We predicted 10 novel circovirus genomes in this study. Whether these circoviruses infect the termites remains elusive at the moment. The functional and taxonomical annotations suggest that the termites may harbor a core virome comprised of the bacteriophages infecting endosymbionts of the gut protozoa. [ABSTRACT FROM AUTHOR]

For over a century, viruses have been known as the most abundant and diverse group of organisms on Earth, forming a virosphere. Based on extensive meta-analyses, we present, for the first time, a wide and complete overview of virus–host network, covering all known viral species. Our data indicate that most of known viral species, regardless of their genomic category, have an intriguingly narrow host range, infecting only 1 or 2 host species. Our data also show that the known virosphere has expanded based on viruses of human interest, related to economical, medical or biotechnological activities. In addition, we provide an overview of the distribution of viruses on different environments on Earth, based on meta-analyses of available metaviromic data, showing the contrasting ubiquity of head-tailed phages against the specificity of some viral groups in certain environments. Finally, we uncovered all human viral species, exploring their diversity and the most affected organic systems. The virus–host network presented here shows an anthropocentric view of the virology. It is therefore clear that a huge effort and change in perspective is necessary to see more than the tip of the iceberg when it comes to virology.

The current view of virus diversity in terrestrial hot springs is limited to a few sampling sites. To expand our current understanding of hot spring viral community diversity, this study aimed to investigate the first African hot spring (Brandvlei hot spring; 60 °C, pH 5.7) by means of electron microscopy and sequencing of the virus fraction. Microscopy analysis revealed a mixture of regular- and 'jumbo'-sized tailed morphotypes (Caudovirales), lemon-shaped virions (Fuselloviridae-like; salterprovirus-like) and pleiomorphic virus-like particles. Metavirome analysis corroborated the presence of His1-like viruses and has expanded the current clade of salterproviruses using a polymerase B gene phylogeny. The most represented viral contig was to a cyanophage genome fragment, which may underline basic ecosystem functioning provided by these viruses. Furthermore, a putative Gemmata-related phage was assembled with high coverage, a previously undocumented phage-host association. This study demonstrated that a moderately thermophilic spring environment contained a highly novel pool of viruses and should encourage future characterization of a wider temperature range of hot springs throughout the world. [ABSTRACT FROM AUTHOR]

Despite being potentially highly useful for characterizing the biodiversity of phages, metagenomic studies are currently not available for dairy bacteriophages, partly due to the lack of a standard procedure for phage extraction. We optimized an extraction method that allows the removal of the bulk protein from whey and milk samples with losses of less than 50% of spiked phages. The protocol was applied to extract phages from whey in order to test the notion that members of Lactococcus lactis 936 (now Sk1virus), P335, c2 (now C2virus) and Leuconostoc phage groups are the most frequently encountered in the dairy environment. The relative abundance and diversity of phages in eight and four whey mixtures from dairies using undefined mesophilic mixed-strain cultures containing Lactococcus lactis subsp. lactis biovar diacetylactis and Leuconostoc species (i.e., DL starter cultures) and defined cultures, respectively, were assessed. Results obtained from transmission electron microscopy and high-throughput sequence analyses revealed the dominance of Lc. lactis 936 phages (order Caudovirales, family Siphoviridae) in dairies using undefined DL starter cultures and Lc. lactis c2 phages (order Caudovirales, family Siphoviridae) in dairies using defined cultures. The 936 and Leuconostoc phages demonstrated limited diversity. Possible coinduction of temperate P335 prophages and satellite phages in one of the whey mixtures was also observed. [ABSTRACT FROM AUTHOR]

For over a century, viruses have been known as the most abundant and diverse group of organisms on Earth, forming a virosphere. Based on extensive meta-analyses, we present, for the first time, a wide and complete overview of virus-host network, covering all known viral species. Our data indicate that most of known viral species, regardless of their genomic category, have an intriguingly narrow host range, infecting only 1 or 2 host species. Our data also show that the known virosphere has expanded based on viruses of human interest, related to economical, medical or biotechnological activities. In addition, we provide an overview of the distribution of viruses on different environments on Earth, based on meta-analyses of available metaviromic data, showing the contrasting ubiquity of head-tailed phages against the specificity of some viral groups in certain environments. Finally, we uncovered all human viral species, exploring their diversity and the most affected organic systems. The virus-host network presented here shows an anthropocentric view of the virology. It is therefore clear that a huge effort and change in perspective is necessary to see more than the tip of the iceberg when it comes to virology. [ABSTRACT FROM AUTHOR]

The current view of virus diversity in terrestrial hot springs is limited to a few sampling sites. To expand our current understanding of hot spring viral community diversity, this study aimed to investigate the first African hot spring (Brandvlei hot spring; 60 °C, pH 5.7) by means of electron microscopy and sequencing of the virus fraction. Microscopy analysis revealed a mixture of regular- and ‘jumbo’-sized tailed morphotypes (Caudovirales), lemon-shaped virions (Fuselloviridae-like; salterprovirus-like) and pleiomorphic virus-like particles. Metavirome analysis corroborated the presence of His1-like viruses and has expanded the current clade of salterproviruses using a polymerase B gene phylogeny. The most represented viral contig was to a cyanophage genome fragment, which may underline basic ecosystem functioning provided by these viruses. Furthermore, a putative Gemmata-related phage was assembled with high coverage, a previously undocumented phage-host association. This study demonstrated that a moderately thermophilic spring environment contained a highly novel pool of viruses and should encourage future characterization of a wider temperature range of hot springs throughout the world.

Microbial starter cultures are used in the production of many cheeses around the world, such as Parmigiano-Reggiano, in Italy, Époisses, in France, and Canastra, in Brazil, providing many of the unique features of these cheeses. Bacteriophages (phages) are ubiquitous and well known to modulate the structure of bacterial communities, and recent data indicate that cheeses contain a high abundance of naturally occurring phages. Here, we analyze the viral and bacterial metagenomes of Canastra cheese: a traditional artisanal Brazilian cheese produced using an endogenous starter culture and raw milk. Over 1,200 viral operational taxonomic units were recovered using both isolated viral-like particles and complete metagenomic DNA. Common viral families identified included Siphoviridae and Myoviridae , with 40% of putative phage genomes unidentified at the family level of classification. We observed very high phage diversity, which varied greatly across different cheese producers, with 28% of phage genomes detected in only one producer. Several metagenome-assembled genomes were recovered for lactic acid-producing bacteria, as well as nonstarter bacterial species, and we identified several phage-bacterium interactions, at the strain level of resolution, varying across distinct cheese producers. We postulate that at least one bacterial strain detected could be endogenous and unique to the Canastra cheese-producing region in Brazil and that its growth seems to be modulated by autochthonous phages present in this artisanal production system. This phage-host relationship is likely to influence the fermentation dynamics and ultimately the sensorial profile of these cheeses, with implications for other similar cheese production systems around the world. IMPORTANCE Our work demonstrated a dynamic yet stable microbial ecosystem during cheese production using an endogenous starter culture. This was observed across several distinct producers and was marked by genomic evidence of continued phage-bacterium interactions, such as the presence of bacterial defense mechanisms. Furthermore, we provide evidence of unique microbial signatures for each individual cheese producer studied in the region, a fact that may have profound consequences on product traceability. This was the first effort to describe and understand the bacteriophage composition and ecological dynamics within the Brazilian Canastra cheese production system. The study of this prototypical backslopping production system provides a solid background for further mechanistic studies of the production of many cheeses around the world.

The article discusses a study which investigated the sequencing of marine metagenomic fosmids that resulted in the discovery of new complete phage genomes. Topics discussed include the sequencing of a much larger set of fosmids, the existence of a biogeography of microbes, and the complexity of phage behavior. Also cited is the need to tap the largest genetic reservoir of the Earth using the advances in sequencing.

The diversity and function of viruses in coral holobionts has only recently received attention. The non-reef building gorgonian octocoral, Gorgonia ventalina, is a major constituent of Caribbean reefs. We investigated viral communities associated with G. ventalina tissues to understand their role in gorgonian ecology. Pyrosequencing was used to prepare a total of 514,632 sequence reads of DNA- and RNA-based mixed-community viral genomes (metaviromes). RNA viral assemblages were comprised of primarily unidentifiable reads, with most matching host transcripts and other RNA metaviromes. DNA metaviromes were similar between healthy and diseased tissues and comprised of contiguous sequences (contigs) that matched primarily metazoan and bacterial proteins. Only ~5% of contigs matched viral proteins that were primarily cyanophage and viruses of Chlorella and Ostreococcus. Our results confirm that DNA and RNA viruses comprise a component of the gorgonian holobiont, suggesting that they may play a role in the ecology of G. ventalina. [ABSTRACT FROM AUTHOR]

Sponges are remarkable holobionts harboring extremely diverse microbial and viral communities. However, the interactions between the components within holobionts and between a holobiont and environment are largely unknown, especially for polar organisms. To investigate possible interactions within and between sponge-associated communities, we probed the microbiomes and viromes of cold-water sympatric sponges Isodictya palmata (n = 2), Halichondria panicea (n = 3), and Halichondria sitiens (n = 3) by 16S and shotgun metagenomics. We showed that the bacterial and viral communities associated with these White Sea sponges are species-specific and different from the surrounding water. Extensive mining of bacterial antiphage defense systems in the metagenomes revealed a variety of defense mechanisms. The abundance of defense systems was comparable in the metagenomes of the sponges and the surrounding water, thus distinguishing the White Sea sponges from those inhabiting the tropical seas. We developed a network-based approach for the combined analysis of CRISPR-spacers and protospacers. Using this approach, we showed that the virus–host interactions within the sponge-associated community are typically more abundant (three out of four interactions studied) than the inter-community interactions. Additionally, we detected the occurrence of viral exchanges between the communities. Our work provides the first insight into the metagenomics of the three cold-water sponge species from the White Sea and paves the way for a comprehensive analysis of the interactions between microbial communities and associated viruses. [ABSTRACT FROM AUTHOR]

The COVID-19 pandemic caused by SARS-CoV-2 has become a major threat across the globe. Here, we developed machine learning approaches to identify key pathogenic regions in coronavirus genomes. We trained and evaluated 7,562,625 models on 3,665 genomes including SARS-CoV-2, MERS-CoV, SARS-CoV, and other coronaviruses of human and animal origins to return quantitative and biologically interpretable signatures at nucleotide and amino acid resolutions. We identified hotspots across the SARS-CoV-2 genome, including previously unappreciated features in spike, RdRp, and other proteins. Finally, we integrated pathogenicity genomic profiles with B cell and T cell epitope predictions for enrichment of sequence targets to help guide vaccine development. These results provide a systematic map of predicted pathogenicity in SARS-CoV-2 that incorporates sequence, structural, and immunologic features, providing an unbiased collection of genetic elements for functional studies. This metavirome-based framework can also be applied for rapid characterization of new coronavirus strains or emerging pathogenic viruses., Competing Interests: The authors declare no competing interests., (© 2021 The Author(s).)

Sponges are remarkable holobionts harboring extremely diverse microbial and viral communities. However, the interactions between the components within holobionts and between a holobiont and environment are largely unknown, especially for polar organisms. To investigate possible interactions within and between sponge-associated communities, we probed the microbiomes and viromes of cold-water sympatric sponges Isodictya palmata (n = 2), Halichondria panicea (n = 3), and Halichondria sitiens (n = 3) by 16S and shotgun metagenomics. We showed that the bacterial and viral communities associated with these White Sea sponges are species-specific and different from the surrounding water. Extensive mining of bacterial antiphage defense systems in the metagenomes revealed a variety of defense mechanisms. The abundance of defense systems was comparable in the metagenomes of the sponges and the surrounding water, thus distinguishing the White Sea sponges from those inhabiting the tropical seas. We developed a network-based approach for the combined analysis of CRISPR-spacers and protospacers. Using this approach, we showed that the virus-host interactions within the sponge-associated community are typically more abundant (three out of four interactions studied) than the inter-community interactions. Additionally, we detected the occurrence of viral exchanges between the communities. Our work provides the first insight into the metagenomics of the three cold-water sponge species from the White Sea and paves the way for a comprehensive analysis of the interactions between microbial communities and associated viruses.