Ocean plankton. Patterns and ecological drivers of ocean viral communities

A Novel and Ubiquitous Marine Methylophage Provides Insights into Viral-Host Coevolution and Possible Host-Range Expansion in Streamlined Marine Heterotrophic Bacteria

The methylotrophic OM43 clade are Gammaproteobacteria that comprise some of the smallest free-living cells known and have highly streamlined genomes. OM43 represents an important microbial link between marine primary production and remineralization of carbon back to the atmosphere. Bacteriophages shape microbial communities and are major drivers of mortality and global marine biogeochemistry. Recent cultivation efforts have brought the first viruses infecting members of the OM43 clade into culture. Here, we characterize a novel myophage infecting OM43 called Melnitz. Melnitz was isolated independently from water samples from a subtropical ocean gyre (Sargasso Sea) and temperate coastal (Western English Channel) systems. Metagenomic recruitment from global ocean viromes confirmed that Melnitz is globally ubiquitous, congruent with patterns of host abundance. Bacteria with streamlined genomes such as OM43 and the globally dominant SAR11 clade use riboswitches as an efficient method to regulate metabolism. Melnitz encodes a two-piece tmRNA (ssrA), controlled by a glutamine riboswitch, providing evidence that riboswitch use also occurs for regulation during phage infection of streamlined heterotrophs. Virally encoded tRNAs and ssrA found in Melnitz were phylogenetically more closely related to those found within the alphaproteobacterial SAR11 clade and their associated myophages than those within their gammaproteobacterial hosts. This suggests the possibility of an ancestral host transition event between SAR11 and OM43. Melnitz and a related myophage that infects SAR11 were unable to infect hosts of the SAR11 and OM43, respectively, suggesting host transition rather than a broadening of host range. IMPORTANCE Isolation and cultivation of viruses are the foundations on which the mechanistic understanding of virus-host interactions and parameterization of bioinformatic tools for viral ecology are based. This study isolated and characterized the first myophage known to infect the OM43 clade, expanding our knowledge of this understudied group of microbes. The nearly identical genomes of four strains of Melnitz isolated from different marine provinces and the global abundance estimations from metagenomic data suggest that this viral population is globally ubiquitous. Genome analysis revealed several unusual features in Melnitz and related genomes recovered from viromes, such as a curli operon and virally encoded tmRNA controlled by a glutamine riboswitch, neither of which are found in the host. Further phylogenetic analysis of shared genes indicates that this group of viruses infecting the gammaproteobacterial OM43 shares a recent common ancestor with viruses infecting the abundant alphaproteobacterial SAR11 clade. Host ranges are affected by compatible cell surface receptors, successful circumvention of superinfection exclusion systems, and the presence of required accessory proteins, which typically limits phages to singular narrow groups of closely related bacterial hosts. This study provides intriguing evidence that for streamlined heterotrophic bacteria, virus-host transitioning may not be necessarily restricted to phylogenetically related hosts but is a function of shared physical and biochemical properties of the cell.

There's more to RNA viruses than diseases

A simple, pervasive biological entity in the ocean sheds light on evolution.

Comparative Analysis of the Gut Microbiota of Mongolian Gazelle (Procapra gutturosa) Under Fragmented Habitats

The existence of man-made facilities such as pasture fences makes the grassland ecosystem fragmented and endangers the survival of local wild animals. The Mongolian gazelle is highly sensitive to hunting and habitat destruction, and is one of the most threatened artiodactyls in Eurasia. It provides a critical model to studying gut microbiota under fragmented habitats. Therefore, we applied metagenomics sequencing to analyze the gut microbiota communities and functions of Mongolian gazelle under fragmented habitats. The results demonstrated that there were no significant differences in gut microbial communities between the different groups at both the phylum and genus level. The functional analyses showed that the Mongolian gazelle in fragmented habitat had a stronger ability to degrade naphthalene, but their ability to absorb carbohydrates was weaker. This study provided fundamental information about the gut microbiota of Mongolian gazelle, and we recommend reducing habitat fragmentation to better protect the Mongolian gazelle.

Viral community analysis in a marine oxygen minimum zone indicates increased potential for viral manipulation of microbial physiological state

Microbial communities in oxygen minimum zones (OMZs) are known to have significant impacts on global biogeochemical cycles, but viral influence on microbial processes in these regions are much less studied. Here we provide baseline ecological patterns using microscopy and viral metagenomics from the Eastern Tropical North Pacific (ETNP) OMZ region that enhance our understanding of viruses in these climate-critical systems. While extracellular viral abundance decreased below the oxycline, viral diversity and lytic infection frequency remained high within the OMZ, demonstrating that viral influences on microbial communities were still substantial without the detectable presence of oxygen. Viral community composition was strongly related to oxygen concentration, with viral populations in low-oxygen portions of the water column being distinct from their surface layer counterparts. However, this divergence was not accompanied by the expected differences in viral-encoded auxiliary metabolic genes (AMGs) relating to nitrogen and sulfur metabolisms that are known to be performed by microbial communities in these low-oxygen and anoxic regions. Instead, several abundant AMGs were identified in the oxycline and OMZ that may modulate host responses to low-oxygen stress. We hypothesize that this is due to selection for viral-encoded genes that influence host survivability rather than modulating host metabolic reactions within the ETNP OMZ. Together, this study shows that viruses are not only diverse throughout the water column in the ETNP, including the OMZ, but their infection of microorganisms has the potential to alter host physiological state within these biogeochemically important regions of the ocean.

An Inducible Microbacterium Prophage vB_MoxS-R1 Represents a Novel Lineage of Siphovirus

Lytic and lysogenic infections are the main strategies used by viruses to interact with microbial hosts. The genetic information of prophages provides insights into the nature of phages and their potential influences on hosts. Here, the siphovirus vB_MoxS-R1 was induced from a Microbacterium strain isolated from an estuarine Synechococcus culture. vB_MoxS-R1 has a high replication capability, with an estimated burst size of 2000 virions per cell. vB_MoxS-R1 represents a novel phage genus-based genomic analysis. Six transcriptional regulator (TR) genes were predicted in the vB_MoxS-R1 genome. Four of these TR genes are involved in stress responses, virulence and amino acid transportation in bacteria, suggesting that they may play roles in regulating the host cell metabolism in response to external environmental changes. A glycerophosphodiester phosphodiesterase gene related to phosphorus acquisition was also identified in the vB_MoxS-R1 genome. The presence of six TR genes and the phosphorus-acquisition gene suggests that prophage vB_MoxS-R1 has the potential to influence survival and adaptation of its host during lysogeny. Possession of four endonuclease genes in the prophage genome suggests that vB_MoxS-R1 is likely involved in DNA recombination or gene conversion and further influences host evolution.

Microbial Diversity in the Indian Ocean Sediments: An Insight into the Distribution and Associated Factors

Indian Ocean is the third largest oceanic division of the world and shelter to a huge microbial diversity. These microbes play an important role in the metabolism of carbon, sulfur, nitrogen, and phosphorus in the ocean water. They are also major contributors of carbon fixing and sequestration, as much as terrestrial plants to achieve CO₂ emissions reduction. The prokaryotic community in the East Indian Ocean primarily comprises of heterotrophic bacteria like Alphaproteobacteria and Gammaproteobacteria, followed by Firmicutes and Actinobacteria. The Arabian Sea and the Bay of Bengal are typically characterized by presence of vast areas of oxygen minimum zones (OMZs) and have been witnessing a shift in the microbial diversity due to the changing conditions in the ocean water. Several canonical correspondence analyses reveal temperature, salinity, and phosphate levels as crucial environmental factors in propelling the distribution of diazotrophs. The viral consortia are dominated by the Caudovirales, an order of tailed bacteriophages. Due to the rapid change in the environmental factors such as topography, temperature, and sunlight contributing toward climate change, their role in sustaining the chemical composition of the ocean can be drastically affected especially with the evidence of several bacterial and fungal communities responding to latitudinal and temperature change. Therefore, we aim to critically review the status of microbial diversity in Indian Ocean to predict their response toward climate change as they are the sentinels of change in marine life and to understand the dynamics of microbial communities in the various locations of Indian Ocean.

Metagenomics-based insights into the microbial community profiling and flavor development potentiality of baijiu Daqu and huangjiu wheat Qu

Daqu and wheat Qu are saccharification and fermenting agents in Chinese huangjiu and baijiu production. This study aimed to investigate the difference between Daqu and wheat Qu in physicochemical indices, microbial communities, functional genes, and the metabolic network of key microbes responsible for flavor synthesis by whole-metagenome sequencing and metabolite analysis. Herein, physicochemical indices indicated that compared with wheat Qu, Daqu exhibited higher protease and cellulase activity and acidity, and lower glucoamylase and amylase enzyme activity. Metagenomic sequencing reveals that although Daqu and wheat Qu community composition have significant differences at species level, they have similar functional genes. Daqu were enriched in Pediococcus pentosaceus, Weissella paramesenteroides, Rasamsonia emersonii and Byssochlamys spectabilis (22.48% of the total abundance), while wheat Qu harbored greater abundances of Saccharopolyspora (54.78%, Saccharopolyspora rectivirgula, Saccharopolyspora shandongensis, Saccharopolyspora hirsuta, Saccharopolyspora spinose, and Saccharopolyspora erythraea). From a functional perspective, the important functions of Daqu and wheat Qu are both amino acid metabolism and carbohydrate metabolism. Meanwhile, a combined analysis among microbiota, functional genes, and dominant flavors indicated S. shandongensis, S. rectivirgula, and S. spinose might be the main contributor to the synthesis of flavor compounds in wheat Qu, while R. emersonii, W. paramesenteroides, Leuconostoc citreum, Leuconostoc mesenteroides, Weissella cibaria and P. pentosaceus may make the greatest contribution to flavor compounds synthesis in Daqu. This study reveals the microbial and functional dissimilarities of Daqu and wheat Qu, and helps elucidating different metabolic roles of microbes during flavor formation.

Independent host- and bacterium-based determinants protect a model symbiosis from phage predation

Bacteriophages (phages) are diverse and abundant constituents of microbial communities worldwide, capable of modulating bacterial populations in diverse ways. Here, we describe the phage HNL01, which infects the marine bacterium Vibrio fischeri. We use culture-based approaches to demonstrate that mutations in the exopolysaccharide locus of V. fischeri render this bacterium resistant to infection by HNL01, highlighting the extracellular matrix as a key determinant of HNL01 infection. Additionally, using the natural symbiosis between V. fischeri and the squid Euprymna scolopes, we show that, during colonization, V. fischeri is protected from phages present in the ambient seawater. Taken together, these findings shed light on independent yet synergistic host- and bacterium-based strategies for resisting symbiosis-disrupting phage predation, and we present important implications for understanding these strategies in the context of diverse host-associated microbial ecosystems.

A holistic genome dataset of bacteria, archaea and viruses of the Pearl River estuary

Estuaries are one of the most important coastal ecosystems. While microbiomes and viromes have been separately investigated in some estuaries, few studies holistically deciphered the genomes and connections of viruses and their microbial hosts along an estuarine salinity gradient. Here we applied deep metagenomic sequencing on microbial and viral communities in surface waters of the Pearl River estuary, one of China’s largest estuaries with strong anthropogenic impacts. Overall, 1,205 non-redundant prokaryotic genomes with ≥50% completeness and ≤10% contamination, and 78,502 non-redundant viral-like genomes were generated from samples of three size fractions and five salinity levels. Phylogenomic analysis and taxonomy classification show that majority of these estuarine prokaryotic and viral genomes are novel at species level according to public databases. Potential connections between the microbial and viral populations were further investigated by host-virus matching. These combined microbial and viral genomes provide an important complement of global marine genome datasets and should greatly facilitate our understanding of microbe-virus interactions, evolution and their implications in estuarine ecosystems.

Measurement(s)	bacteria • Archaea • viruses
Technology Type(s)	Shotgun Sequencing
Sample Characteristic - Organism	estuary metagenome
Sample Characteristic - Environment	subtropical estuarine
Sample Characteristic - Location	Pearl river estuary

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.17139998

Metagenomics of the midgut microbiome of Rhipicephalus microplus from China

Background

Ticks, which are ectoparasites of animals, may carry multiple pathogens. The cattle tick Rhipicephalus microplus is an important bovine parasite in China. However, the midgut microbiome of R. microplus from China has not been characterized via metagenomic methods.

Methods

Rhipicephalus microplus were collected from cattle in the city of Changsha in Hunan province, China. The DNA of the midgut contents was extracted from fully engorged adult female R. microplus. A DNA library was constructed and sequenced using an Illumina HiSeq sequencing platform. SOAPdenovo software was used to assemble and analyze the clean data. The latent class analysis algorithm applied to system classification by MEGAN software was used to annotate the information on the species’ sequences. DIAMOND software was used to compare unigenes with the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, and functional annotation was carried out based on the results of the comparison.

Results

The dominant phyla in the five samples were Firmicutes, Proteobacteria, and Actinobacteria. Streptococcus, Mycobacterium, Anaplasma, Enterococcus, Shigella, Lactobacillus, Brachyspira, Pseudomonas, Enterobacter, Bacillus, and Lactococcus were the dominant genera in the five samples. The endosymbiotic bacterium Wolbachia was also detected in all of the samples. Mycobacterium malmesburyense, Streptococcus pneumoniae, Anaplasma phagocytophilum, Enterococcus faecium, Shigella sonnei, Enterococcus faecalis, Lactobacillus casei, Brachyspira hampsonii, Pseudomonas syringae, Enterobacter cloacae, and Lactococcus garvieae were the dominant species in the five samples. In addition to these bacterial species, we also detected some eukaryotes, such as Rhizophagus irregularis, Enterospora canceri, Smittium culicis, Zancudomyces culisetae, Trachipleistophora hominis, and viruses such as orf virus, human endogenous retrovirus type W, enzootic nasal tumor virus of goats, bovine retrovirus CH15, and galidia endogenous retrovirus in all of the samples at the species level. The results of the annotated KEGG pathway predictions for the gene functions of the midgut microflora of R. microplus indicated genes involved in lipid and amino acid metabolism, infectious diseases (e.g., Streptococcuspneumonia infection, human granulocytic anaplasmosis, Shigellasonnei infection, Salmonella enterica infection, and pathogenic Escherichia coli infection), and cancer.

Conclusions

Our study revealed that the midgut microbiome of R. microplus is not only composed of a large number of bacteria, but that a portion also comprises eukaryotes and viruses. The data presented here enhance our understanding of this tick’s midgut microbiome and provide fundamental information for the control of ticks and tick-borne diseases.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05161-6.

Genomic Characterization of a Novel Freshwater Cyanophage Reveals a New Lineage of Cyanopodovirus

Cyanobacteria are one of the dominant autotrophs in tropical freshwater communities, yet phages infecting them remain poorly characterized. Here we present the characterization of cyanophage S-SRP02, isolated from a tropical freshwater lake in Singapore, which infects Synechococcus sp. Strain SR-C1 isolated from the same lake. S-SRP02 represents a new evolutionary lineage of cyanophage. Out of 47 open reading frames (ORFs), only 20 ORFs share homology with genes encoding proteins of known function. There is lack of auxiliary metabolic genes which was commonly found as core genes in marine cyanopodoviruses. S-SRP02 also harbors unique structural genes highly divergent from other cultured phages. Phylogenetic analysis and viral proteomic tree further demonstrate the divergence of S-SRP02 from other sequenced phage isolates. Nonetheless, S-SRP02 shares synteny with phage genes of uncultured phages obtained from the Mediterranean Sea deep chlorophyll maximum fosmids, indicating the ecological importance of S-SRP02 and its related viruses. This is further supported by metagenomic mapping of environmental viral metagenomic reads onto the S-SRP02 genome.

Temporal and Spatial Signaling Mediating the Balance of the Plankton Microbiome

The annual patterns of plankton succession in the ocean determine ecological and biogeochemical cycles. The temporally fluctuating interplay between photosynthetic eukaryotes and the associated microbiota balances the composition of aquatic planktonic ecosystems. In addition to nutrients and abiotic factors, chemical signaling determines the outcome of interactions between phytoplankton and their associated microbiomes. Chemical mediators control essential processes, such as the development of key morphological, physiological, behavioral, and life-history traits during algal growth. These molecules thus impact species succession and community composition across time and space in processes that are highlighted in this review. We focus on spatial, seasonal, and physiological dynamics that occur during the early association of algae with bacteria, the exponential growth of a bloom, and its decline and recycling. We also discuss how patterns from field data and global surveys might be linked to the actions of metabolic markers in natural phytoplankton assemblages.

Putative Invertebrate, Plant, and Wastewater Derived ssRNA Viruses in Plankton of the Anthropogenically Impacted Anacostia River, District of Columbia, USA

The Anacostia River is a highly impacted watershed in the Northeastern United States which experiences combined sewage outfall in downstream waters. We examined the composition of RNA viruses at three sites in the river using viral metagenomics. Viromes had well represented Picornaviruses, Tombusviruses, Wolframviruses, Nodaviruses, with fewer Tobamoviruses, Sobemoviruses, and Densoviruses (ssDNA). Phylogenetic ana-lyses of detected viruses provide evidence for putatively autochthonous and allochthonous invertebrate, plant, and vertebrate host origin. The number of viral genomes matching Ribovaria increased downstream, and assemblages were most disparate between distant sites, suggesting impacts of the combined sewage overflows at these sites. Additionally, we recovered a densovirus genome fragment which was highly similar to the Clinch ambidensovirus 1, which has been attributed to mass mortality of freshwater mussels in Northeastern America. Taken together, these data suggest that RNA viromes of the Anacostia River reflect autochthonous production of virus particles by benthic metazoan and plants, and inputs from terrestrial habitats including sewage.

A catalogue of 1,167 genomes from the human gut archaeome

The human gut microbiome plays an important role in health, but its archaeal diversity remains largely unexplored. In the present study, we report the analysis of 1,167 nonredundant archaeal genomes (608 high-quality genomes) recovered from human gastrointestinal tract, sampled across 24 countries and rural and urban populations. We identified previously undescribed taxa including 3 genera, 15 species and 52 strains. Based on distinct genomic features, we justify the split of the Methanobrevibacter smithii clade into two separate species, with one represented by the previously undescribed 'Candidatus Methanobrevibacter intestini'. Patterns derived from 28,581 protein clusters showed significant associations with sociodemographic characteristics such as age groups and lifestyle. We additionally show that archaea are characterized by specific genomic and functional adaptations to the host and carry a complex virome. Our work expands our current understanding of the human archaeome and provides a large genome catalogue for future analyses to decipher its impact on human physiology.

Tianwang Buxin Granules Influence the Intestinal Flora in Perimenopausal Insomnia

Methods

The subjects included 13 PI patients from the Hubei Provincial Hospital of TCM, Hubei University of TCM, and Wuhan Traditional Chinese Medicine Hospital, and the corresponding noninsomniac spouses of the patients were selected as controls. TWBXG was continuously administered for 4 weeks. The feces of PI patients and their noninsomniac spouses before and after treatment with TWBXG were collected. The intestinal flora composition of each group was detected by metagenomic sequencing, and the efficacy of TWBXG was evaluated by the PSQI scale.

Results

Compared with the control group, the model group showed an increase in the abundance of Roseburia faecis, Ruminococcus, Prevotella copri, Fusicatenibacter saccharivorans, and Blautia obeum, while those of Bacteroides, fecal Bacteroidetes, and Faecalibacterium prausnitzii were decreased. Compared with pretreatment, the PSQI score was significantly reduced (P < 0.05), the abundance of Bacteroides, fecal Bacteroidetes, and Faecalibacterium prausnitzii increased, and that of Roseburia faecis, Ruminococcus, Prevotella copri, Fusicatenibacter saccharivorans, and Blautia obeum decreased after treatment. However, there was still a certain gap in the abundance of related flora in the treatment group compared with the control.

Conclusion

PI is associated with disturbances in the intestinal flora and is mainly related to the disorders of Roseburia faecis, Ruminococcus, Prevotella copri, Fusicatenibacter saccharivorans, Blautia obeum, Bacteroides, fecal Bacteroidetes, and Faecalibacterium prausnitzii. TWBXG can effectively treat PI, and its effect may be achieved by regulating the disordered intestinal flora. Clinical Trials. The study was registered in the Chinese clinical trial registry and approved by the World Health Organization clinical trial registration platform (Effects of the modified Tianwang Buxin granule and modified Tianwang Buxin decoction pieces on insomnia: a randomized, controlled trial, ChiCTR-IPR-17011549).

Unveiling Ecological and Genetic Novelty within Lytic and Lysogenic Viral Communities of Hot Spring Phototrophic Microbial Mats

Viruses exert diverse ecosystem impacts by controlling their host community through lytic predator-prey dynamics. However, the mechanisms by which lysogenic viruses influence their host-microbial community are less clear. In hot springs, lysogeny is considered an active lifestyle, yet it has not been systematically studied in all habitats, with phototrophic microbial mats (PMMs) being particularly not studied. We carried out viral metagenomics following in situ mitomycin C induction experiments in PMMs from Porcelana hot spring (Northern Patagonia, Chile). The compositional changes of viral communities at two different sites were analyzed at the genomic and gene levels. Furthermore, the presence of integrated prophage sequences in environmental metagenome-assembled genomes from published Porcelana PMM metagenomes was analyzed. Our results suggest that virus-specific replicative cycles (lytic and lysogenic) were associated with specific host taxa with different metabolic capacities. One of the most abundant lytic viral groups corresponded to cyanophages, which would infect the cyanobacteria Fischerella, the most active and dominant primary producer in thermophilic PMMs. Likewise, lysogenic viruses were related exclusively to chemoheterotrophic bacteria from the phyla Proteobacteria, Firmicutes, and Actinobacteria. These temperate viruses possess accessory genes to sense or control stress-related processes in their hosts, such as sporulation and biofilm formation. Taken together, these observations suggest a nexus between the ecological role of the host (metabolism) and the type of viral lifestyle in thermophilic PMMs. This has direct implications in viral ecology, where the lysogenic-lytic switch is determined by nutrient abundance and microbial density but also by the metabolism type that prevails in the host community.

IMPORTANCE Hot springs harbor microbial communities dominated by a limited variety of microorganisms and, as such, have become a model for studying community ecology and understanding how biotic and abiotic interactions shape their structure. Viruses in hot springs are shown to be ubiquitous, numerous, and active components of these communities. However, lytic and lysogenic viral communities of thermophilic phototrophic microbial mats (PMMs) remain largely unexplored. In this work, we use the power of viral metagenomics to reveal changes in the viral community following a mitomycin C induction experiment in PMMs. The importance of our research is that it will improve our understanding of viral lifestyles in PMMs via exploring the differences in the composition of natural and induced viral communities at the genome and gene levels. This novel information will contribute to deciphering which biotic and abiotic factors may control the transitions between lytic and lysogenic cycles in these extreme environments.

Antimicrobial resistance bacteria and genes detected in hospital sewage provide valuable information in predicting clinical antimicrobial resistance

Extensive use of antibiotics is significantly associated with development of antibiotic-resistant (AR) bacteria. However, their causal relationships have not been adequately investigated, especially in human population and hospitals. Our aims were to understand clinical AR through revealing co-occurrence patterns between antibiotic-resistant bacteria and genes (ARB and ARGs), and their association with antibiotic use, and to consider impact of ARB and ARGs on environmental and human health. Antibiotic usage was calculated based on the actual consumption in our target hospital. ARB was identified by culture. In isolates collected from hospital sewage, bacterial-specific DNA sequences and ARGs were determined using metagenomics. Our data revealed that the use of culture-based single-indicator-strain approaches only captured ARB in 16.17% of the infectious samples. On the other hand, 1573 bacterial species and 885 types of ARGs were detected in the sewage. Furthermore, hospital use of antibiotics influenced the resistance profiles, but the strength varied among bacteria. From our metagenomics analyses, ARGs for aminoglycosides were the most common, followed by sulfonamide, tetracycline, phenicol, macrolides, and quinolones, comprising 82.6% of all ARGs. Association analyses indicated that 519 pairs of ARGs were significantly correlated with ARB species (r > 0.8). The co-occurrence patterns of bacteria-ARGs mirrored the AR in the clinic. In conclusion, our systematic investigation further emphasized that antibiotic usage in hospital significantly influenced the abundance and types of ARB and ARGs in dose- and time-dependent manners which, in turn, mirrored clinical AR. In addition, our data provide novel information on development of certain ARB with multiple antibiotic resistance. These ARB and ARGs from sewage can also be disseminated into the environment and communities to create health problems. Therefore, it would be helpful to use such data to develop improved predictive risk model of AR, to enhance effective use of antibiotics, and to reduce environmental pollution.

Viruses infecting a warm water picoeukaryote shed light on spatial co-occurrence dynamics of marine viruses and their hosts

The marine picoeukaryote Bathycoccus prasinos has been considered a cosmopolitan alga, although recent studies indicate two ecotypes exist, Clade BI (B. prasinos) and Clade BII. Viruses that infect Bathycoccus Clade BI are known (BpVs), but not that infect BII. We isolated three dsDNA prasinoviruses from the Sargasso Sea against Clade BII isolate RCC716. The BII-Vs do not infect BI, and two (BII-V2 and BII-V3) have larger genomes (~210 kb) than BI-Viruses and BII-V1. BII-Vs share ~90% of their proteins, and between 65% to 83% of their proteins with sequenced BpVs. Phylogenomic reconstructions and PolB analyses establish close-relatedness of BII-V2 and BII-V3, yet BII-V2 has 10-fold higher infectivity and induces greater mortality on host isolate RCC716. BII-V1 is more distant, has a shorter latent period, and infects both available BII isolates, RCC716 and RCC715, while BII-V2 and BII-V3 do not exhibit productive infection of the latter in our experiments. Global metagenome analyses show Clade BI and BII algal relative abundances correlate positively with their respective viruses. The distributions delineate BI/BpVs as occupying lower temperature mesotrophic and coastal systems, whereas BII/BII-Vs occupy warmer temperature, higher salinity ecosystems. Accordingly, with molecular diagnostic support, we name Clade BII Bathycoccus calidus sp. nov. and propose that molecular diversity within this new species likely connects to the differentiated host-virus dynamics observed in our time course experiments. Overall, the tightly linked biogeography of Bathycoccus host and virus clades observed herein supports species-level host specificity, with strain-level variations in infection parameters.

Genomic Characterization of Two Novel RCA Phages Reveals New Insights into the Diversity and Evolution of Marine Viruses

Viruses are the most abundant living entities in marine ecosystems, playing critical roles in altering the structure and function of microbial communities and driving ocean biogeochemistry. Phages that infect Roseobacter clade-affiliated (RCA) cluster strains are an important component of marine viral communities. Here, we characterize the genome sequences of two new RCA phages, CRP-9 and CRP-13, which infect RCA strain FZCC0023. Genomic analysis reveals that CRP-9 and CRP-13 represent a novel evolutionary lineage of marine phages. They both have a DNA replication module most similar to those in Cobavirus group phages. In contrast, their morphogenesis and packaging modules are distinct from those in cobaviruses but homologous to those in HMO-2011-type phages. The genomic architecture of CRP-9 and CRP-13 suggests a genomic recombination event between distinct phage groups. Metagenomic data sets were examined for metagenome-assembled viral genomes (MAVGs) with similar recombinant genome architectures. Fifteen CRP-9-type MAVGs were identified from marine viromes. Additionally, 158 MAVGs were identified containing HMO-2011-type morphogenesis and packaging modules with other types of DNA replication genes, providing more evidence that recombination between different phage groups is a major driver of phage evolution. Altogether, this study significantly expands the understanding of diversity and evolution of marine roseophages. Meanwhile, the analysis of these novel RCA phages and MAVGs highlights the critical role of recombination in shaping phage diversity. These phage sequences are valuable resources for inferring the evolutionary connection of distinct phage groups. IMPORTANCE Diversity and evolution of phages that infect the relatively slow-growing but dominant Roseobacter lineages are largely unknown. In this study, RCA phages CRP-9 and CRP-13 have been isolated on a Roseobacter RCA strain and shown to have a unique genomic architecture, which appears to be the result of a recombination event. CRP-9 and CRP-13 have a DNA replication module most similar to those in Cobavirus group phages and morphogenesis and packaging modules most similar to those in HMO-2011-type phages. HMO-2011-type morphogenesis and packaging modules are found in combination with distinct types of DNA replication genes, suggesting compatibility with various DNA replication modules. Altogether, this study contributes toward a better understanding of marine viral diversity and evolution.

Viral Characteristics of the Warm Atlantic and Cold Arctic Water Masses in the Nordic Seas

Nordic Seas are the subarctic seas connecting the Arctic Ocean and North Atlantic Ocean with complex water masses, experiencing an abrupt climate change. Though knowledge of the marine virosphere has expanded rapidly, the diversity of viruses and their relationships with host cells and water masses in the Nordic Seas remain to be fully revealed. Here, we establish the Nordic Sea DNA virome (NSV) data set of 55,315 viral contigs including 1,478 unique viral populations from seven stations influenced by both the warm Atlantic and cold Arctic water masses. Caudovirales dominated in the seven NSVs, especially in the warm Atlantic waters. The major giant nucleocytoplasmic large DNA viruses (NCLDVs) contributed a significant proportion of the classified viral contigs in the NSVs (32.2%), especially in the cold Arctic waters (44.9%). The distribution patterns of Caudovirales and NCLDVs were a reflection of the community structure of their hosts in the corresponding water masses and currents. Latitude, pH, and flow speed were found to be key factors influencing the microbial communities and coinfluencing the variation of viral communities. Network analysis illustrated the tight coupling between the variation of viral communities and microbial communities in the Nordic Seas. This study suggests a probable linkage between viromes, host cells, and surface water masses from both the cool Arctic and warm Atlantic Oceans. IMPORTANCE This is a systematic study of Nordic Sea viromes using metagenomic analysis. The viral diversity, community structure, and their relationships with host cells and the complex water masses from both the cool Arctic and the warm Atlantic oceans were illustrated. The NCLDVs and Caudovirales are proposed as the viral characteristics of the cold Arctic and warm Atlantic waters, respectively. This study provides an important background for the viromes in the subarctic seas connecting the Arctic Ocean and North Atlantic Ocean and sheds light on their responses to abrupt climate change in the future.

Immune lag is a major cost of prokaryotic adaptive immunity during viral outbreaks

Clustered regularly interspaced short palindromic repeat (CRISPR)-Cas adaptive immune systems enable bacteria and archaea to efficiently respond to viral pathogens by creating a genomic record of previous encounters. These systems are broadly distributed across prokaryotic taxa, yet are surprisingly absent in a majority of organisms, suggesting that the benefits of adaptive immunity frequently do not outweigh the costs. Here, combining experiments and models, we show that a delayed immune response which allows viruses to transiently redirect cellular resources to reproduction, which we call ‘immune lag’, is extremely costly during viral outbreaks, even to completely immune hosts. Critically, the costs of lag are only revealed by examining the early, transient dynamics of a host–virus system occurring immediately after viral challenge. Lag is a basic parameter of microbial defence, relevant to all intracellular, post-infection antiviral defence systems, that has to-date been largely ignored by theoretical and experimental treatments of host-phage systems.

DNase Treatment Improves Viral Enrichment in Agricultural Soil Viromes

The small genomes of most viruses make it difficult to fully capture viral diversity in metagenomes dominated by DNA from cellular organisms. Viral size fraction metagenomics (viromics) protocols facilitate the enrichment of viral DNA from environmental samples, and these protocols typically include DNase treatment of the post-0.2-μm-filtered viromic fraction to remove contaminating free DNA prior to virion lysis. However, DNase may also remove desirable viral genomic DNA (e.g., contained in virions compromised due to frozen storage or laboratory processing), suggesting that DNase-untreated viromes might be useful in some cases. In order to understand how virome preparation with and without DNase treatment influences the resultant data, here, we compared 15 soil viromes (7 DNase treated and 8 untreated) from 8 samples collected from agricultural fields prior to tomato planting. DNase-treated viromes yielded significantly more assembled viral contigs, contained significantly less nonviral microbial DNA, and recovered more viral populations (viral operational taxonomic units [vOTUs]) through read mapping. However, DNase-treated and untreated viromes were statistically indistinguishable in terms of ecological patterns across viral communities. Although the results suggest that DNase treatment is preferable where possible, in comparison to previously reported total metagenomes from the same samples, both DNase-treated and untreated viromes were significantly enriched in viral signatures by all metrics compared, including a 225-times-higher proportion of viral reads in untreated viromes compared to total metagenomes. Thus, even without DNase treatment, viromics was preferable to total metagenomics for capturing viral diversity in these soils, suggesting that preparation of DNase-untreated viromes can be worthwhile when DNase treatment is not possible. IMPORTANCE Viromics is becoming an increasingly popular method for characterizing soil viral communities. DNase treatment of the viral size fraction prior to DNA extraction is meant to reduce contaminating free DNA and is a common step within viromics protocols to ensure that sequences are of viral origin. However, some samples may not be amenable to DNase treatment due to viral particles being compromised either in storage (i.e., frozen) or during other sample processing steps. To date, the effect of DNase treatment on the recovery of viruses and downstream ecological interpretations of soil viral communities is not thoroughly understood. This work sheds light on these questions and indicates that while DNase treatment of soil viromes improves the recovery of viral populations, this improvement is modest in comparison to the gains made by viromics over total soil metagenomics. Furthermore, DNase treatment may not be necessary to observe the ecological patterns structuring soil viral communities.

Imaging Techniques for Detecting Prokaryotic Viruses in Environmental Samples

Viruses are the most abundant biological entities on Earth with an estimate of 1031 viral particles across all ecosystems. Prokaryotic viruses-bacteriophages and archaeal viruses-influence global biogeochemical cycles by shaping microbial communities through predation, through the effect of horizontal gene transfer on the host genome evolution, and through manipulating the host cellular metabolism. Imaging techniques have played an important role in understanding the biology and lifestyle of prokaryotic viruses. Specifically, structure-resolving microscopy methods, for example, transmission electron microscopy, are commonly used for understanding viral morphology, ultrastructure, and host interaction. These methods have been applied mostly to cultivated phage-host pairs. However, recent advances in environmental genomics have demonstrated that the majority of viruses remain uncultivated, and thus microscopically uncharacterized. Although light- and structure-resolving microscopy of viruses from environmental samples is possible, quite often the link between the visualization and the genomic information of uncultivated prokaryotic viruses is missing. In this minireview, we summarize the current state of the art of imaging techniques available for characterizing viruses in environmental samples and discuss potential links between viral imaging and environmental genomics for shedding light on the morphology of uncultivated viruses and their lifestyles in Earth's ecosystems.

A new family of "megaphages" abundant in the marine environment

Megaphages, bacteriophages harbouring extremely large genomes, have recently been found to be ubiquitous, being described from a variety of microbiomes ranging from the animal gut to soil and freshwater systems. However, no complete marine megaphage has been identified to date. Here, using both short and long read sequencing, we assembled >900 high-quality draft viral genomes from water in the English Channel. One of these genomes included a novel megaphage, Mar_Mega_1 at >650 Kb, making it one of the largest phage genomes assembled to date. Utilising phylogenetic and network approaches, we found this phage represents a new family of megaphages. Genomic analysis showed Mar_Mega_1 shares relatively few homologues with its closest relatives, but, as with other megaphages Mar_Mega_1 contained a variety of auxiliary metabolic genes responsible for carbon metabolism and nucleotide biosynthesis, including a NADP-dependent isocitrate dehydrogenase [Idh] and nicotinamide-nucleotide amidohydrolase [PncC], which have not previously been identified in megaphages. Mar_Mega_1 was abundant in a marine virome sample and related phages are widely prevalent in the oceans.

Metagenomic Analysis Identifies Sex-Related Cecal Microbial Gene Functions and Bacterial Taxa in the Quail

Background: Japanese quail (Coturnix japonica) are important and widely distributed poultry in China. Researchers continue to pursue genetic selection for heavier quail. The intestinal microbiota plays a substantial role in growth promotion; however, the mechanisms involved in growth promotion remain unclear.

Results: We generated 107.3 Gb of cecal microbiome data from ten Japanese quail, providing a series of quail gut microbial gene catalogs (1.25 million genes). We identified a total of 606 main microbial species from 1,033,311 annotated genes distributed among the ten quail. Seventeen microbial species from the genera Anaerobiospirillum, Alistipes, Barnesiella, and Butyricimonas differed significantly in their abundances between the female and male gut microbiotas. Most of the functional gut microbial genes were involved in metabolism, primarily in carbohydrate transport and metabolism, as well as some active carbohydrate-degrading enzymes. We also identified 308 antibiotic-resistance genes (ARGs) from the phyla Bacteroidetes, Firmicutes and Euryarchaeota. Studies of the differential gene functions between sexes indicated that abundances of the gut microbes that produce carbohydrate-active enzymes varied between female and male quail. Bacteroidetes was the predominant ARG-containing phylum in female quail; Euryarchaeota was the predominant ARG-containing phylum in male quail.

Conclusion: This article provides the first description of the gene catalog of the cecal bacteria in Japanese quail as well as insights into the bacterial taxa and predictive metagenomic functions between male and female quail to provide a better understanding of the microbial genes in the quail ceca.

Cyanolichen microbiome contains novel viruses that encode genes to promote microbial metabolism

Lichen thalli are formed through the symbiotic association of a filamentous fungus and photosynthetic green alga and/or cyanobacterium. Recent studies have revealed lichens also host highly diverse communities of secondary fungal and bacterial symbionts, yet few studies have examined the viral component within these complex symbioses. Here, we describe viral biodiversity and functions in cyanolichens collected from across North America and Europe. As current machine-learning viral-detection tools are not trained on complex eukaryotic metagenomes, we first developed efficient methods to remove eukaryotic reads prior to viral detection and a custom pipeline to validate viral contigs predicted with three machine-learning methods. Our resulting high-quality viral data illustrate that every cyanolichen thallus contains diverse viruses that are distinct from viruses in other terrestrial ecosystems. In addition to cyanobacteria, predicted viral hosts include other lichen-associated bacterial lineages and algae, although a large fraction of viral contigs had no host prediction. Functional annotation of cyanolichen viral sequences predicts numerous viral-encoded auxiliary metabolic genes (AMGs) involved in amino acid, nucleotide, and carbohydrate metabolism, including AMGs for secondary metabolism (antibiotics and antimicrobials) and fatty acid biosynthesis. Overall, the diversity of cyanolichen AMGs suggests that viruses may alter microbial interactions within these complex symbiotic assemblages.

A systematic review from basics to omics on bacteriophage applications in poultry production and processing

The growing human population is currently facing an unprecedented challenge on global food production and sustainability. Despite recognizing poultry as one of the most successful and rapidly growing food industries to address this challenge; poultry health and safety remain major issues that entail immediate attention. Bacterial diseases including colibacillosis, salmonellosis, and necrotic enteritis have become increasingly prevalent during poultry production. Likewise, outbreaks caused by consumption of undercooked poultry products contaminated with zoonotic bacterial pathogens such as Salmonella, Campylobacter and Listeria, are a serious public health concern. With antimicrobial resistance problem and restricted use of antibiotics in food producing animals, bacteriophages are increasingly recognized as an attractive natural antibacterial alternative. Bacteriophages have recently shown promising results to treat diseases in poultry, reduce contamination of carcasses, and enhance the safety of poultry products. Omics technologies have been successfully employed to accurately characterize bacteriophages and their genes/proteins important for interaction with bacterial hosts. In this review, the potential of using lytic bacteriophages to mitigate the risk of major poultry-associated bacterial pathogens are explored. This study also explores challenges associated with the adoption of this technology by industries. Furthermore, the impact of omics approaches on studying bacteriophages, their host interaction and applications is discussed.

Diversity and distribution of viruses inhabiting the deepest ocean on Earth

As the most abundant biological entities on the planet, viruses significantly influence the overall functioning of marine ecosystems. The abundance, distribution, and biodiversity of viral communities in the upper ocean have been relatively well studied, but our understanding of viruses in the hadal biosphere remains poor. Here, we established the oceanic trench viral genome dataset (OTVGD) by analysing 19 microbial metagenomes derived from seawater and sediment samples of the Mariana, Yap, and Kermadec Trenches. The trench viral communities harbored remarkably high novelty, and they were predicted to infect ecologically important microbial clades, including Thaumarchaeota and Oleibacter. Significant inter-trench and intra-trench exchange of viral communities was proposed. Moreover, viral communities in different habitats (seawater/sediment and depth-stratified ocean zones) exhibited distinct niche-dependent distribution patterns and genomic properties. Notably, microbes and viruses in the hadopelagic seawater seemed to preferably adopt lysogenic lifestyles compared to those in the upper ocean. Furthermore, niche-specific auxiliary metabolic genes were identified in the hadal viral genomes, and a novel viral D-amino acid oxidase was functionally and phylogenetically characterized, suggesting the contribution of these genes in the utilization of refractory organic matter. Together, these findings highlight the genomic novelty, dynamic movement, and environment-driven diversification of viral communities in oceanic trenches, and suggest that viruses may influence the hadal ecosystem by reprogramming the metabolism of their hosts and modulating the community of keystone microbes.

Analysis of Microorganism Diversity in Haemaphysalis longicornis From Shaanxi, China, Based on Metagenomic Sequencing

Ticks are dangerous ectoparasites of humans and animals, as they are important disease vectors and serve as hosts for various microorganisms (including a variety of pathogenic microorganisms). Diverse microbial populations coexist within the tick body. Metagenomic next-generation sequencing (mNGS) has been suggested to be useful for rapidly and accurately obtaining microorganism abundance and diversity data. In this study, we performed mNGS to analyze the microbial diversity of Haemaphysalis longicornis from Baoji, Shaanxi, China, with the Illumina HiSeq platform. We identified 189 microbial genera (and 284 species) from ticks in the region; the identified taxa included Anaplasma spp., Rickettsia spp., Ehrlichia spp., and other important tick-borne pathogens at the genus level as well as symbiotic microorganisms such as Wolbachia spp., and Candidatus Entotheonella. The results of this study provide insights into possible tick-borne diseases and reveal new tick-borne pathogens in this region. Additionally, valuable information for the biological control of ticks is provided. In conclusion, this study provides reference data for guiding the development of prevention and control strategies targeting ticks and tick-borne diseases in the region, which can improve the effectiveness of tick and tick-borne disease control.

High-throughput identification of viral termini and packaging mechanisms in virome datasets using PhageTermVirome

Viruses that infect bacteria (phages) are increasingly recognized for their importance in diverse ecosystems but identifying and annotating them in large-scale sequence datasets is still challenging. Although efficient scalable virus identification tools are emerging, defining the exact ends (termini) of phage genomes is still particularly difficult. The proper identification of termini is crucial, as it helps in characterizing the packaging mechanism of bacteriophages and provides information on various aspects of phage biology. Here, we introduce PhageTermVirome (PTV) as a tool for the easy and rapid high-throughput determination of phage termini and packaging mechanisms using modern large-scale metagenomics datasets. We successfully tested the PTV algorithm on a mock virome dataset and then used it on two real virome datasets to achieve the rapid identification of more than 100 phage termini and packaging mechanisms, with just a few hours of computing time. Because PTV allows the identification of free fully formed viral particles (by recognition of termini present only in encapsidated DNA), it can also complement other virus identification softwares to predict the true viral origin of contigs in viral metagenomics datasets. PTV is a novel and unique tool for high-throughput characterization of phage genomes, including phage termini identification and characterization of genome packaging mechanisms. This software should help researchers better visualize, map and study the virosphere. PTV is freely available for downloading and installation at https://gitlab.pasteur.fr/vlegrand/ptv.

Identification of tick-borne pathogens using metagenomic analyses in H. longicornis feeding on humans in downtown Beijing

On August 14th, 2018, a Beijing resident living in Xicheng District found a female H. longicornis tick attached to the skin at the front of his upper shin. On examination, the patient was afebrile and appeared well. The species of the tick was identified through morphological characteristics and phylogenetic analysis based on cytochrome C oxidase subunit I. This H. longicornis tick was screened for tick-borne pathogens such as viruses, bacteria and parasites. RNA pathogens were screened by PCR and sequencing, while DNA pathogens were screened by metagenomic analyses. It was found that the tick was positive for the DNA sequences of zoonotic and animal pathogens such as A. phagocytophilum, Ehrlichia minasensis and C. burnetii. Considering the good health condition of the patient, we hypothesized that the pathogens originated from the tick specimen itself rather than host blood meal. For the first time, our study reveals the possible risk of transmission of tick-borne pathogens to human beings through tick bit in downtown Beijing. Further research is needed to screen for tick-borne pathogens among unfed ticks collected from central Beijing.

Visualizing active viral infection reveals diverse cell fates in synchronized algal bloom demise

Marine viruses are the most abundant biological entity in the ocean and are considered as major evolutionary drivers of microbial life [C. A. Suttle, Nat. Rev. Microbiol. 5, 801-812 (2007)]. Yet, we lack quantitative approaches to assess their impact on the marine ecosystem. Here, we provide quantification of active viral infection in the bloom forming single-celled phytoplankton Emiliania huxleyi infected by the large virus EhV, using high-throughput single-molecule messenger RNA in situ hybridization (smFISH) of both virus and host transcripts. In natural samples, viral infection reached only 25% of the population despite synchronized bloom demise exposing the coexistence of infected and noninfected subpopulations. We prove that photosynthetically active cells chronically release viral particles through nonlytic infection and that viral-induced cell lysis can occur without viral release, thus challenging major assumptions regarding the life cycle of giant viruses. We could also assess active infection in cell aggregates linking viral infection and carbon export to the deep ocean [C. P. Laber et al., Nat. Microbiol. 3, 537-547 (2018)] and suggest a potential host defense strategy by enrichment of infected cells in sinking aggregates. Our approach can be applied to diverse marine microbial systems, opening a mechanistic dimension to the study of biotic interactions in the ocean.

Ecogenomics and Adaptation Strategies of Southern Ocean Viral Communities

The Southern Ocean (SO) represents up to one-fifth of the total carbon drawdown worldwide. Intense selective pressures (low temperature, high UV radiation, and strong seasonality) and physical isolation characterize the SO, serving as a "natural" laboratory for the study of ecogenomics and unique adaptations of endemic viral populations. Here, we report 2,416 novel viral genomes from the SO, obtained from newly sequenced viral metagenomes in combination with mining of publicly available data sets, which represents a 25% increase in the SO viral genomes reported to date. They comprised 567 viral clusters (defined as approximately genus-level groups), with 186 genera endemic to the SO, demonstrating that the SO viral community is predominantly constituted by a large pool of genetically divergent viral species from widespread viral families. The predicted proteome from SO viruses revealed that several protein clusters related to cold-shock-event responses and quorum-sensing mechanisms involved in the lysogenic-lytic cycle shift decision were under positive selection, which is ultimately important for fine adaptation of viral populations in response to the strong selective pressures of the SO. Finally, changes in the hydrophobicity patterns and amino acid frequencies suggested marked temperature-driven genetic selection of the SO viral proteome. Our data provide valuable insights into how viruses adapt and remain successful in this extreme polar marine environment. IMPORTANCE Viruses are the most abundant biologic entities in marine systems and strongly influence the microbial community composition and diversity. However, little is known about viral communities' adaptation and diversification in the ocean. In this work, we take advantage of the geographical isolation and the intense selective pressures of the SO, to which viruses are exposed, to identify potential viral adaptations due to positive environmental selection and dispersal limitation. To that end, we recovered more than two thousand novel viral genomes, revealing a high degree of divergence in these SO endemic communities. Furthermore, we describe remarkable viral adaptations in amino acid frequencies and accessory proteins related to cold shock response and quorum sensing that allow them to thrive at lower temperatures. Consequently, our work greatly expands the understanding of the diversification of the viral communities of the SO and their particular adaptations to low temperatures.

Spatial and Temporal Dynamics of Prokaryotic and Viral Community Assemblages in a Lotic System (Manatee Springs, Florida)

Flow from high-magnitude springs fed by the Floridan aquifer system contributes hundreds of liters of water per second to rivers, creating unique lotic systems. Despite their importance as freshwater sources and their contributions to the state's major rivers, little is known about the composition and spatiotemporal variability of prokaryotic and viral communities of these spring systems or their influence on downstream river sites. At four time points throughout a year, we determined the abundance and diversity of prokaryotic and viral communities at three sites within the first-magnitude Manatee Springs system (the spring head where water emerges from the aquifer, a mixed region where the spring run ends, and a downstream site in the Suwannee River). The abundance of prokaryotes and virus-like particles increased 100-fold from the spring head to the river and few members from the head communities persisted in the river at low abundance, suggesting the springs play a minor role in seeding downstream communities. Prokaryotic and viral communities within Manatee Springs clustered by site, with seasonal variability likely driven by flow. As water flowed through the system, microbial community composition was affected by changes in physiochemical parameters and community coalescence. Evidence of species sorting and mass effects could be seen in the assemblages. Greater temporal fluctuations were observed in prokaryotic and viral community composition with increasing distance from the spring outflow, reflecting the relative stability of the groundwater environment, and comparisons to springs from prior work reaffirmed that distinct first-magnitude springs support unique communities. IMPORTANCE Prokaryotic and viral communities are central to food webs and biogeochemical processes in aquatic environments, where they help maintain ecosystem health. The Floridan aquifer system (FAS), which is the primary drinking water source for millions of people in the southeastern United States, contributes large amounts of freshwater to major river systems in Florida through its springs. However, there is a paucity of information regarding the spatiotemporal dynamics of microbial communities in these essential flowing freshwater systems. This work explored the prokaryotic and viral communities in a first-magnitude spring system fed by the FAS that discharges millions of liters of water per day into the Suwannee River. This study examined microbial community composition through space and time as well as the environmental parameters and metacommunity assembly mechanisms that shape these communities, providing a foundational understanding for monitoring future changes.

Metagenomics-Based Analysis of the Age-Related Cumulative Effect of Antibiotic Resistance Genes in Gut Microbiota

Antibiotic resistance in bacteria has become a major global health problem. One of the main reservoirs of antibiotic resistance genes is the human gut microbiota. To characterise these genes, a metagenomic approach was used. In this study, a comprehensive antibiotic resistome catalog was established using fecal samples from 246 healthy individuals from world's longevity township in Jiaoling, China. In total, 606 antibiotic resistance genes were detected. Our results indicated that antibiotic resistance genes in the human gut microbiota accumulate and become more complex with age as older groups harbour the highest abundance of these genes. Tetracycline resistance gene type tetQ was the most abundant group of antibiotic resistance genes in gut microbiota, and the main carrier of antibiotic resistance genes was Bacteroides. Antibiotic efflux, inactivation, and target alteration were found to be the dominant antimicrobial resistance mechanisms. This research may help to establish a comprehensive antibiotic resistance catalog that includes extremely long-lived healthy people such as centenarians, and may provide potential recommendations for controlling the use of antibiotics.

Revealing RNA virus diversity and evolution in unicellular algae transcriptomes

Remarkably little is known about the diversity and evolution of RNA viruses in unicellular eukaryotes. We screened a total of 570 transcriptomes from the Marine Microbial Eukaryote Transcriptome Sequencing Project that encompasses a wide diversity of microbial eukaryotes, including most major photosynthetic lineages (i.e. the microalgae). From this, we identified thirty new and divergent RNA virus species, occupying a range of phylogenetic positions within the overall diversity of RNA viruses. Approximately one-third of the newly described viruses comprised single-stranded positive-sense RNA viruses from the order Lenarviricota associated with fungi, plants, and protists, while another third were related to the order Ghabrivirales, including members of the protist and fungi-associated Totiviridae. Other viral species showed sequence similarity to positive-sense RNA viruses from the algae-associated Marnaviridae, the double-stranded RNA (ds-RNA) Partitiviridae, as well as tentative evidence for one negative-sense RNA virus related to the Qinviridae. Importantly, we were able to identify divergent RNA viruses from distant host taxa, revealing the ancestry of these viral families and greatly extending our knowledge of the RNA viromes of microalgal cultures. Both the limited number of viruses detected per sample and the low sequence identity to known RNA viruses imply that additional microalgal viruses exist that could not be detected at the current sequencing depth or were too divergent to be identified using sequence similarity. Together, these results highlight the need for further investigation of algal-associated RNA viruses as well as the development of new tools to identify RNA viruses that exhibit very high levels of sequence divergence.

Identifying Potential Hosts of Short-Branch Microsporidia

Microsporidia are obligate parasites that are closely related to Fungi. While the widely known "long-branch" Microsporidia infect mostly metazoans, the hosts of "short-branch" Microsporidia are only partially characterized or not known at all. Here, we used network analyses from Neotropical rainforest soil metabarcoding data, to infer co-occurrences between environmental lineages of short-branch microsporidians and their potential hosts. We found significant co-occurrences with several taxa, especially with Apicomplexa, Cercozoa, and Fungi, as well as some Metazoa. Our results are the first step to identify potential hosts of the environmental lineages of short-branch microsporidians, which can be targeted in future molecular and microscopic studies.

Expanding our understanding of marine viral diversity through metagenomic analyses of biofilms

Recent metagenomics surveys have provided insights into the marine virosphere. However, these surveys have focused solely on viruses in seawater, neglecting those associated with biofilms. By analyzing 1.75 terabases of biofilm metagenomic data, 3974 viral sequences were identified from eight locations around the world. Over 90% of these viral sequences were not found in previously reported datasets. Comparisons between biofilm and seawater metagenomes identified viruses that are endemic to the biofilm niche. Analysis of viral sequences integrated within biofilm-derived microbial genomes revealed potential functional genes for trimeric autotransporter adhesin and polysaccharide metabolism, which may contribute to biofilm formation by the bacterial hosts. However, more than 70% of the genes could not be annotated. These findings show marine biofilms to be a reservoir of novel viruses and have enhanced our understanding of natural virus-bacteria ecosystems.

Glacier ice archives nearly 15,000-year-old microbes and phages

BACKGROUND

Glacier ice archives information, including microbiology, that helps reveal paleoclimate histories and predict future climate change. Though glacier-ice microbes are studied using culture or amplicon approaches, more challenging metagenomic approaches, which provide access to functional, genome-resolved information and viruses, are under-utilized, partly due to low biomass and potential contamination.

RESULTS

We expand existing clean sampling procedures using controlled artificial ice-core experiments and adapted previously established low-biomass metagenomic approaches to study glacier-ice viruses. Controlled sampling experiments drastically reduced mock contaminants including bacteria, viruses, and free DNA to background levels. Amplicon sequencing from eight depths of two Tibetan Plateau ice cores revealed common glacier-ice lineages including Janthinobacterium, Polaromonas, Herminiimonas, Flavobacterium, Sphingomonas, and Methylobacterium as the dominant genera, while microbial communities were significantly different between two ice cores, associating with different climate conditions during deposition. Separately, ~355- and ~14,400-year-old ice were subject to viral enrichment and low-input quantitative sequencing, yielding genomic sequences for 33 vOTUs. These were virtually all unique to this study, representing 28 novel genera and not a single species shared with 225 environmentally diverse viromes. Further, 42.4% of the vOTUs were identifiable temperate, which is significantly higher than that in gut, soil, and marine viromes, and indicates that temperate phages are possibly favored in glacier-ice environments before being frozen. In silico host predictions linked 18 vOTUs to co-occurring abundant bacteria (Methylobacterium, Sphingomonas, and Janthinobacterium), indicating that these phages infected ice-abundant bacterial groups before being archived. Functional genome annotation revealed four virus-encoded auxiliary metabolic genes, particularly two motility genes suggest viruses potentially facilitate nutrient acquisition for their hosts. Finally, given their possible importance to methane cycling in ice, we focused on Methylobacterium viruses by contextualizing our ice-observed viruses against 123 viromes and prophages extracted from 131 Methylobacterium genomes, revealing that the archived viruses might originate from soil or plants.

CONCLUSIONS

Together, these efforts further microbial and viral sampling procedures for glacier ice and provide a first window into viral communities and functions in ancient glacier environments. Such methods and datasets can potentially enable researchers to contextualize new discoveries and begin to incorporate glacier-ice microbes and their viruses relative to past and present climate change in geographically diverse regions globally. Video Abstract.

Responses of bacterial and bacteriophage communities to long-term exposure to antimicrobial agents in wastewater treatment systems

The occurrence of antibacterial agents has received increasing concern due to their possible threats to human health. However, the effects of antibacterial residues on the evolution and dynamics between bacteria and bacteriophages in wastewater treatment systems have seldom been researched. Especially for phages, little is known about their response to antimicrobial exposure. In this study, two identical anoxic-aerobic wastewater treatment systems were established to evaluate the responses of bacterial and phage communities to long-term exposure to antimicrobial agents. The results indicated simultaneous exposure to combined antimicrobials significantly inhibited (p < 0.05) the abundance of phages and bacteria. Metagenomic sequencing analysis indicated the community of bacteria and phages changed greatly at the genus level due to combined antibacterial exposure. Additionally, long-term exposure to antimicrobial agents promoted the attachment of receptor-binding protein genes to Klebsiella, Escherichia and Salmonella (which were all members of Enterobacteriaceae). Compared to that in the control system, the numbers of receptor-binding protein genes on their possible phages (such as Lambdalikevirus and P2likevirus) were also obviously higher when the microorganisms were exposed to antimicrobials. The results are helpful to understanding the microbial communities and tracking the relationship of phage-bacterial host systems, especially under the pressure of antimicrobial exposure.

Monsavirus in monkey rectal swab and throat swab specimens in China: Proposal for Posaliviridae as a new family in Picornavirales

Posa-like viruses have been detected in the fecal samples of several host species and are considered unclassified members of Picornavirales. Here, we identified genomic fragments of novel posa-like viruses (monsaviruses) in monkey specimens through next generation sequencing and obtained 11 full-length genomes. This monsavirus shared 88.5-89.2% nucleotide similarity with the Tottori-HG1 strain (GenBank accession LC123275). In total, 713 nucleotide polymorphism sites were identified, indicating their persistent evolution during circulation. The genomic organization and phylogenetic relationship of monsavirus were determined. Subsequent phylogenetic analysis of the conserved replication block of Hel-Pro-RdRp and core RNA-dependent RNA polymerase domain-based analysis of posa-like viruses showed significant separation compared with other known families. Further, posa-like virus genomes possessed the classical replication block of picornavirus in the 5' part of genome and picorna-like capsid domains at the structural coding region of 3' part of genome. Based on these results, we proposed the new family Posaliviridae, within Picornavirales. Four genera, which showed 68.6-75.5% amino acid distances but similar genomic organization including the conserved replication block of Hel-Pro-RdRp, the same order of the genomic coding region, and picorna-like capsid domains, were identified. The flexible genomic organization strategy and a large evolutionary scale of Posaliviridae was explicit. This study provides novel information on monsaviruses and important taxonomic data for the family Posaliviridae.

Spatial patterns in phage-Rhizobium coevolutionary interactions across regions of common bean domestication

Bacteriophages play significant roles in the composition, diversity, and evolution of bacterial communities. Despite their importance, it remains unclear how phage diversity and phage-host interactions are spatially structured. Local adaptation may play a key role. Nitrogen-fixing symbiotic bacteria, known as rhizobia, have been shown to locally adapt to domesticated common bean at its Mesoamerican and Andean sites of origin. This may affect phage-rhizobium interactions. However, knowledge about the diversity and coevolution of phages with their respective Rhizobium populations is lacking. Here, through the study of four phage-Rhizobium communities in Mexico and Argentina, we show that both phage and host diversity is spatially structured. Cross-infection experiments demonstrated that phage infection rates were higher overall in sympatric rhizobia than in allopatric rhizobia except for one Argentinean community, indicating phage local adaptation and host maladaptation. Phage-host interactions were shaped by the genetic identity and geographic origin of both the phage and the host. The phages ranged from specialists to generalists, revealing a nested network of interactions. Our results suggest a key role of local adaptation to resident host bacterial communities in shaping the phage genetic and phenotypic composition, following a similar spatial pattern of diversity and coevolution to that in the host.

Dynamics of Baltic Sea phages driven by environmental changes

Phage predation constitutes a major mortality factor for bacteria in aquatic ecosystems, and thus, directly impacts nutrient cycling and microbial community dynamics. Yet, the population dynamics of specific phages across time scales from days to months remain largely unexplored, which limits our understanding of their influence on microbial succession. To investigate temporal changes in diversity and abundance of phages infecting particular host strains, we isolated 121 phage strains that infected three bacterial hosts during a Baltic Sea mesocosm experiment. Genome analysis revealed a novel Flavobacterium phage genus harboring gene sets putatively coding for synthesis of modified nucleotides and glycosylation of bacterial cell surface components. Another novel phage genus revealed a microdiversity of phage species that was largely maintained during the experiment and across mesocosms amended with different nutrients. In contrast to the newly described Flavobacterium phages, phages isolated from a Rheinheimera strain were highly similar to previously isolated genotypes, pointing to genomic consistency in this population. In the mesocosm experiment, the investigated phages were mainly detected after a phytoplankton bloom peak. This concurred with recurrent detection of the phages in the Baltic Proper during summer months, suggesting an influence on the succession of heterotrophic bacteria associated with phytoplankton blooms.

Diverse Viruses in Deep-Sea Hydrothermal Vent Fluids Have Restricted Dispersal across Ocean Basins

In the ocean, viruses impact microbial mortality, regulate biogeochemical cycling, and alter the metabolic potential of microbial lineages. At deep-sea hydrothermal vents, abundant viruses infect a wide range of hosts among the archaea and bacteria that inhabit these dynamic habitats. However, little is known about viral diversity, host range, and biogeography across different vent ecosystems, which has important implications for how viruses manipulate microbial function and evolution. Here, we examined viral diversity, viral and host distribution, and virus-host interactions in microbial metagenomes generated from venting fluids from several vent sites within three different geochemically and geographically distinct hydrothermal systems: Piccard and Von Damm vent fields at the Mid-Cayman Rise in the Caribbean Sea, and at several vent sites within Axial Seamount in the Pacific Ocean. Analysis of viral sequences and clustered regularly interspaced short palindromic repeat (CRISPR) spacers revealed highly diverse viral assemblages and evidence of active infection. Network analysis revealed that viral host range was relatively narrow, with very few viruses infecting multiple microbial lineages. Viruses were largely endemic to individual vent sites, indicating restricted dispersal, and in some cases, viral assemblages persisted over time. Thus, we show that hydrothermal vent fluids are home to novel, diverse viral assemblages that are highly localized to specific regions and taxa. IMPORTANCE Viruses play important roles in manipulating microbial communities and their evolution in the ocean, yet not much is known about viruses in deep-sea hydrothermal vents. However, viral ecology and evolution are of particular interest in hydrothermal vent habitats because of their unique nature: previous studies have indicated that most viruses in hydrothermal vents are temperate rather than lytic, and it has been established that rates of horizontal gene transfer (HGT) are particularly high among thermophilic vent microbes, and viruses are common vectors for HGT. If viruses have broad host range or are widespread across vent sites, they have increased potential to act as gene-sharing "highways" between vent sites. By examining viral diversity, distribution, and infection networks across disparate vent sites, this study provides the opportunity to better characterize and constrain the viral impact on hydrothermal vent microbial communities. We show that viruses in hydrothermal vents are diverse and apparently active, but most have restricted host range and are not widely distributed among vent sites. Thus, the impacts of viral infection are likely to be highly localized and constrained to specific taxa in these habitats.

Viral footprints across Gulfs of Kathiawar Peninsula and Arabian Sea: Unraveled from pelagic sediment metagenomic data

Marine biosphere is one of the largest, diverse and dynamic system hosting numerous of microorganisms. Viruses being the most abundant under explored lifeforms in ocean, represent a reservoir of great genetic diversity. We report the metagenomic insights on the viral communities in the deep sediments of the two Gulfs of Gujarat i.e. Gulf of Khambhat and Gulf of Kutch, with one sample from Arabian Sea, treated as open sea control. The viral reads were filtered from the whole dataset, assembled and studied for viral diversity, which was visualized by Pavian. The sequences were checked for the viral abundance, diversity and functionality. The resulting viral taxonomic classification contained 6 orders, 8 families and 47 genera. The results revealed that the phages infecting Cyanobacterium, Bacillus and Vibrio dominated the sediments. Further, it was observed that majority of viral sequences belonged to double-stranded DNA phages. The present study attempts to provide a primary insight of the viral signals and potential genetic content in the Gulfs of Kathiawar.

RaFAH: Host prediction for viruses of Bacteria and Archaea based on protein content

Culture-independent approaches have recently shed light on the genomic diversity of viruses of prokaryotes. One fundamental question when trying to understand their ecological roles is: which host do they infect? To tackle this issue we developed a machine-learning approach named Random Forest Assignment of Hosts (RaFAH), that uses scores to 43,644 protein clusters to assign hosts to complete or fragmented genomes of viruses of Archaea and Bacteria. RaFAH displayed performance comparable with that of other methods for virus-host prediction in three different benchmarks encompassing viruses from RefSeq, single amplified genomes, and metagenomes. RaFAH was applied to assembled metagenomic datasets of uncultured viruses from eight different biomes of medical, biotechnological, and environmental relevance. Our analyses led to the identification of 537 sequences of archaeal viruses representing unknown lineages, whose genomes encode novel auxiliary metabolic genes, shedding light on how these viruses interfere with the host molecular machinery. RaFAH is available at https://sourceforge.net/projects/rafah/.

•

RaFAH was developed to predict the hosts of viruses of Bacteria and Archaea

•

RaFAH displayed comparable or superior performance to other host-prediction tools

•

RaFAH performed well across viromes from eight different ecosystems

•

RaFAH identified hundreds of genomic sequences as derived from viruses of Archaea

Viruses that infect Bacteria and Archaea are ubiquitous and extremely abundant. Recent advances have led to the discovery of many thousands of complete and partial genomes of these biological entities. Understanding the biology of these viruses and how they influence their ecosystems depends on knowing which hosts they infect. We developed a tool that uses data from complete or fragmented genomes to predict the hosts of viruses using a machine-learning approach. Our tool, RaFAH, displayed performance comparable with or superior to that of other host-prediction tools. In addition, it identified hundreds of sequences as derived from the genomes of viruses of Archaea, which are one of the least characterized fractions of the global virosphere.

Expanding standards in viromics: in silico evaluation of dsDNA viral genome identification, classification, and auxiliary metabolic gene curation

BACKGROUND

Viruses influence global patterns of microbial diversity and nutrient cycles. Though viral metagenomics (viromics), specifically targeting dsDNA viruses, has been critical for revealing viral roles across diverse ecosystems, its analyses differ in many ways from those used for microbes. To date, viromics benchmarking has covered read pre-processing, assembly, relative abundance, read mapping thresholds and diversity estimation, but other steps would benefit from benchmarking and standardization. Here we use in silico-generated datasets and an extensive literature survey to evaluate and highlight how dataset composition (i.e., viromes vs bulk metagenomes) and assembly fragmentation impact (i) viral contig identification tool, (ii) virus taxonomic classification, and (iii) identification and curation of auxiliary metabolic genes (AMGs).

RESULTS

The in silico benchmarking of five commonly used virus identification tools show that gene-content-based tools consistently performed well for long (≥3 kbp) contigs, while k-mer- and blast-based tools were uniquely able to detect viruses from short (≤3 kbp) contigs. Notably, however, the performance increase of k-mer- and blast-based tools for short contigs was obtained at the cost of increased false positives (sometimes up to ∼5% for virome and ∼75% bulk samples), particularly when eukaryotic or mobile genetic element sequences were included in the test datasets. For viral classification, variously sized genome fragments were assessed using gene-sharing network analytics to quantify drop-offs in taxonomic assignments, which revealed correct assignations ranging from ∼95% (whole genomes) down to ∼80% (3 kbp sized genome fragments). A similar trend was also observed for other viral classification tools such as VPF-class, ViPTree and VIRIDIC, suggesting that caution is warranted when classifying short genome fragments and not full genomes. Finally, we highlight how fragmented assemblies can lead to erroneous identification of AMGs and outline a best-practices workflow to curate candidate AMGs in viral genomes assembled from metagenomes.

CONCLUSION

Together, these benchmarking experiments and annotation guidelines should aid researchers seeking to best detect, classify, and characterize the myriad viruses 'hidden' in diverse sequence datasets.

Going to extremes - a metagenomic journey into the dark matter of life

The Virus-X-Viral Metagenomics for Innovation Value-project was a scientific expedition to explore and exploit uncharted territory of genetic diversity in extreme natural environments such as geothermal hot springs and deep-sea ocean ecosystems. Specifically, the project was set to analyse and exploit viral metagenomes with the ultimate goal of developing new gene products with high innovation value for applications in biotechnology, pharmaceutical, medical, and the life science sectors. Viral gene pool analysis is also essential to obtain fundamental insight into ecosystem dynamics and to investigate how viruses influence the evolution of microbes and multicellular organisms. The Virus-X Consortium, established in 2016, included experts from eight European countries. The unique approach based on high throughput bioinformatics technologies combined with structural and functional studies resulted in the development of a biodiscovery pipeline of significant capacity and scale. The activities within the Virus-X consortium cover the entire range from bioprospecting and methods development in bioinformatics to protein production and characterisation, with the final goal of translating our results into new products for the bioeconomy. The significant impact the consortium made in all of these areas was possible due to the successful cooperation between expert teams that worked together to solve a complex scientific problem using state-of-the-art technologies as well as developing novel tools to explore the virosphere, widely considered as the last great frontier of life.

Seasonal dynamics of mycoplankton in the Yellow Sea reflect the combined effect of riverine inputs and hydrographic conditions

Little is known about how multiple factors including land-based inputs and ocean currents affect the spatiotemporal distribution of the mycoplankton in coastal regions. To explore the seasonal changes of mycoplanktonic communities and potential environmental drivers, we collected water samples from the Yellow Sea, used here as a model for subtropical sea habitats, in different seasons over two years. Compared with winter and spring, summer exhibited higher levels of fungal richness and community heterogeneity in the water column. The seasonal shifts in mycoplankton diversity and community composition were mainly ascribed to freshwater inputs, the Cold Water Mass and invasion of the Yellow Sea Warm Current. Among the physicochemical variables tested, temperature was the primary determinant of fungal diversity and showed contrasting influences on fungal richness in the surface and bottom waters during summer. In addition, we provide evidence for the community similarity and dissolved nutrients of different water bodies to highlight the potential origin of the Cold Water Mass. Our findings bring new understanding on the factors determining the dynamics of mycoplankton communities by modelling the influence of physicochemical variables and tracking the geographical distribution of certain fungal taxa.