Gut virome of mammals and birds reveals high genetic diversity of the family Microviridae

Metagenomic analysis reveals high diversity of gut viromes in yaks (Bos grunniens) from the Qinghai-Tibet Plateau

The Qinghai-Tibet Plateau (QTP), renowned for its exceptional biological diversity, is home to numerous endemic species. However, research on the virology of vulnerable vertebrates like yaks remains limited. In this study, our objective was to use metagenomics to provide a comprehensive understanding of the diversity and evolution of the gut virome in yak populations across different regions of the QTP. Our findings revealed a remarkably diverse array of viruses in the gut of yaks, including those associated with vertebrates and bacteriophages. Notably, some vertebrate-associated viruses, such as astrovirus and picornavirus, showed significant sequence identity across diverse yak populations. Additionally, we observed differences in the functional profiles of genes carried by the yak gut virome across different regions. Moreover, the virus-bacterium symbiotic network that we discovered holds potential significance in maintaining the health of yaks. Overall, this research expands our understanding of the viral communities in the gut of yaks and highlights the importance of further investigating the interactions between viruses and their hosts. These data will be beneficial for revealing the crucial role that viruses play in the yak gut ecology in future studies.

Diverse Circular DNA Viral Communities in Blood, Oral, and Fecal Samples of Captive Lemurs

Few studies have addressed viral diversity in lemurs despite their unique evolutionary history on the island of Madagascar and high risk of extinction. Further, while a large number of studies on animal viromes focus on fecal samples, understanding viral diversity across multiple sample types and seasons can reveal complex viral community structures within and across species. Groups of captive lemurs at the Duke Lemur Center (Durham, NC, USA), a conservation and research center, provide an opportunity to build foundational knowledge on lemur-associated viromes. We sampled individuals from seven lemur species, i.e., collared lemur (Eulemur collaris), crowned lemur (Eulemur coronatus), blue-eyed black lemur (Eulemur flavifrons), ring-tailed lemur (Lemur catta), Coquerel's sifaka (Propithecus coquereli), black-and-white ruffed lemur (Varecia variegata variegata), and red ruffed lemur (Varecia rubra), across two lemur families (Lemuridae, Indriidae). Fecal, blood, and saliva samples were collected from Coquerel's sifaka and black-and-white ruffed lemur individuals across two sampling seasons to diversify virome biogeography and temporal sampling. Using viral metagenomic workflows, the complete genomes of anelloviruses (n = 4), cressdnaviruses (n = 47), caudoviruses (n = 15), inoviruses (n = 34), and microviruses (n = 537) were determined from lemur blood, feces, and saliva. Many virus genomes, especially bacteriophages, identified in this study were present across multiple lemur species. Overall, the work presented here uses a viral metagenomics approach to investigate viral communities inhabiting the blood, oral cavity, and feces of healthy captive lemurs.

Highly diverse RNA viruses and phage sequences concealed within birds

The diversity of birds in most parts of the world is very high, and thus, they may carry different types of highly differentiated and unknown viruses. Thanks to advanced sequencing technologies, studies on the diversity of bird-associated viruses have increased over the past few years. In this study, a large-scale viral metagenomics survey was performed on cloacal swabs of 2,990 birds from nine provinces of the Chinese mainland. To detect undescribed RNA viruses in birds, more than 1,800 sequences sharing relatively low (<60%) amino acid sequence identity with the best match in the GenBank database were screened. Potentially novel viruses related to vertebrates have been identified, and several potential recombination signals were found. Additionally, hundreds of RNA viral sequences related to plants, fungi, and insects were detected, including previously unknown viruses. Furthermore, we investigated the novelty, functionality, and classification of the phages examined in this study. These viruses occupied topological positions on the evolutionary trees to a certain extent and might form novel putative families, genera, or species, thus providing information to fill the phylogenetic gaps of related viruses. These findings provided new insights into bird-associated viruses, but the interactions among these viruses remain unknown and require further investigation.IMPORTANCEStudying the diversity of RNA viruses in birds and mammals is crucial due to their potential impact on human health and the global ecosystem. Many RNA viruses, such as influenza and coronaviruses, have been shown to cross the species barrier and cause zoonotic diseases. In this metagenomics study involving 2,990 birds from at least 82 species, we identified over 1,800 RNA sequences with distant relationships to known viruses, some of which are rare in birds. The study highlights the scope and diversity of RNA viruses in birds, providing data to predict disease risks and monitor potential viral threats to wildlife, livestock, and human health. This information can aid in the development of strategies for disease prevention and control.

Characterizing the gut phageome and phage-borne antimicrobial resistance genes in pigs

BACKGROUND

Mammalian intestine harbors a mass of phages that play important roles in maintaining gut microbial ecosystem and host health. Pig has become a common model for biomedical research and provides a large amount of meat for human consumption. However, the knowledge of gut phages in pigs is still limited.

RESULTS

Here, we investigated the gut phageome in 112 pigs from seven pig breeds using PhaBOX strategy based on the metagenomic data. A total of 174,897 non-redundant gut phage genomes were assembled from 112 metagenomes. A total of 33,487 gut phage genomes were classified and these phages mainly belonged to phage families such as Ackermannviridae, Straboviridae, Peduoviridae, Zierdtviridae, Drexlerviridae, and Herelleviridae. The gut phages in seven pig breeds exhibited distinct communities and the gut phage communities changed with the age of pig. These gut phages were predicted to infect a broad range of 212 genera of prokaryotes, such as Candidatus Hamiltonella, Mycoplasma, Colwellia, and Lactobacillus. The data indicated that broad KEGG and CAZy functions were also enriched in gut phages of pigs. The gut phages also carried the antimicrobial resistance genes (ARGs) and the most abundant antimicrobial resistance genotype was diaminopyrimidine resistance.

CONCLUSIONS

Our research delineates a landscape for gut phages in seven pig breeds and reveals that gut phages serve as a key reservoir of ARGs in pigs. Video Abstract.

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

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

Microviruses: A World Beyond phiX174

Two decades of metagenomic analyses have revealed that in many environments, small (∼5 kb), single-stranded DNA phages of the family Microviridae dominate the virome. Although the emblematic microvirus phiX174 is ubiquitous in the laboratory, most other microviruses, particularly those of the gokushovirus and amoyvirus lineages, have proven to be much more elusive. This puzzling lack of representative isolates has hindered insights into microviral biology. Furthermore, the idiosyncratic size and nature of their genomes have resulted in considerable misjudgments of their actual abundance in nature. Fortunately, recent successes in microvirus isolation and improved metagenomic methodologies can now provide us with more accurate appraisals of their abundance, their hosts, and their interactions. The emerging picture is that phiX174 and its relatives are rather rare and atypical microviruses, and that a tremendous diversity of other microviruses is ready for exploration.

Using archived and biocollection samples towards deciphering the DNA virus diversity associated with rodent species in the families cricetidae and heteromyidae

Rodentia is the most speciose order of mammals, and they are known to harbor a wide range of viruses. Although there has been significant research on zoonotic viruses in rodents, research on the diversity of other viruses has been limited, especially for rodents in the families Cricetidae and Heteromyidae. In fecal and liver samples of nine species of rodents, we identify 346 distinct circular DNA viral genomes. Of these, a large portion are circular, single-stranded DNA viruses in the families Anelloviridae (n = 3), Circoviridae (n = 5), Genomoviridae (n = 7), Microviridae (n = 297), Naryaviridae (n = 4), Vilyaviridae (n = 15) and in the phylum Cressdnaviricota (n = 13) that cannot be assigned established families. We also identified two large bacteriophages of 36 and 50 kb that are part of the class Caudoviricetes. Some of these viruses are clearly those that infect rodents, however, most of these likely infect various organisms associated with rodents, their environment or their diet.

Exposing new taxonomic variation with inflammation - a murine model-specific genome database for gut microbiome researchers

BACKGROUND

The murine CBA/J mouse model widely supports immunology and enteric pathogen research. This model has illuminated Salmonella interactions with the gut microbiome since pathogen proliferation does not require disruptive pretreatment of the native microbiota, nor does it become systemic, thereby representing an analog to gastroenteritis disease progression in humans. Despite the value to broad research communities, microbiota in CBA/J mice are not represented in current murine microbiome genome catalogs.

RESULTS

Here we present the first microbial and viral genomic catalog of the CBA/J murine gut microbiome. Using fecal microbial communities from untreated and Salmonella-infected, highly inflamed mice, we performed genomic reconstruction to determine the impacts on gut microbiome membership and functional potential. From high depth whole community sequencing (~ 42.4 Gbps/sample), we reconstructed 2281 bacterial and 4516 viral draft genomes. Salmonella challenge significantly altered gut membership in CBA/J mice, revealing 30 genera and 98 species that were conditionally rare and unsampled in non-inflamed mice. Additionally, inflamed communities were depleted in microbial genes that modulate host anti-inflammatory pathways and enriched in genes for respiratory energy generation. Our findings suggest decreases in butyrate concentrations during Salmonella infection corresponded to reductions in the relative abundance in members of the Alistipes. Strain-level comparison of CBA/J microbial genomes to prominent murine gut microbiome databases identified newly sampled lineages in this resource, while comparisons to human gut microbiomes extended the host relevance of dominant CBA/J inflammation-resistant strains.

CONCLUSIONS

This CBA/J microbiome database provides the first genomic sampling of relevant, uncultivated microorganisms within the gut from this widely used laboratory model. Using this resource, we curated a functional, strain-resolved view on how Salmonella remodels intact murine gut communities, advancing pathobiome understanding beyond inferences from prior amplicon-based approaches. Salmonella-induced inflammation suppressed Alistipes and other dominant members, while rarer commensals like Lactobacillus and Enterococcus endure. The rare and novel species sampled across this inflammation gradient advance the utility of this microbiome resource to benefit the broad research needs of the CBA/J scientific community, and those using murine models for understanding the impact of inflammation on the gut microbiome more generally. Video Abstract.

Characterization of a novel picornavirus prevalent in experimental rabbits (Oryctolagus cuniculus)

Here, using viral metagenomic method a novel picornavirus (named UJS-2019picorna, GenBank accession number OP821762) was discovered in fecal and blood samples of experimental rabbits (Oryctolagus cuniculus). The complete genome size of UJS-2019picorna is 7832 bp excluding the poly(A)-tail, with GC content of 44.00% and a nucleotide composition of 28.0% A, 28.0% U, 21.5% G, and 22.5% C. The viral genome has a typical picornavirus organization pattern from the 5'-3' direction: VPg-5' UTR-(L)-P1, (VP4-VP2-VP3-VP1)-P2, (2 A-2B-2C)-P3, (3 A-3B-3C-3D)-3' UTR-poly(A). The P1 region of UJS-2019picorna is related to Erbovirus with amino acid identity of 37.31%, while the P2 and P3 regions are the closest to Bopivirus with amino acid identity of 35.66%-39.53%. According to the Picornaviridae Study Group guidelines, UJS-2019picorna should be presumed to be a new genus belonging to the Picornaviridae family. Epidemiologic study revealed that this novel picornavirus was prevalent in a cohort of experimental rabbits, with prevalence rate of 23.68% (9/38) in feces and 18.4% (7/38) in blood samples. Further work is required to elucidate whether this virus is pathogenic to rabbits and whether it has influence on studies using rabbits as experimental animal.

Virome diversity of ticks feeding on domestic mammals in China

Ticks are considered the second most common pathogen vectors transmitting a broad range of vital human and veterinary viruses. From 2017 to 2018, 640 ticks were collected in eight different provinces in central and western China. Six species were detected, including H.longicornis, De.everestianus, Rh.microplus, Rh.turanicus, Rh.sanguineous, and Hy.asiaticum. Sixty-four viral metagenomic libraries were constructed on the MiSeq Illumina platform, resulting in 13.44 ​G (5.88 ​× ​10⁷) of 250-bp-end reads, in which 2,437,941 are viral reads. We found 27 nearly complete genome sequences, including 16 genome sequences encoding entire protein-coding regions (lack of 3' or 5' end non-coding regions) and complete viral genomes, distributed in the arboviral family (Chuviridae, Rhabdoviridae, Nairoviridae, Phenuiviridae, Flaviviridae, Iflaviridae) as well as Parvoviridae and Polyomaviridae that cause disease in mammals and even humans. In addition, 13 virus sequences found in Chuviridae, Nairoviridae, Flaviviridae, Iflaviridae, Hepeviridae, Parvoviridae, and Polyomaviridae were identified as belonging to a new virus species in the identified viral genera. Besides, an epidemiological survey shows a high prevalence (9.38% and 15.63%) of two viruses (Ovine Copiparvovirus and Bovine parvovirus 2) in the tick cohort.

Viral metagenomics of the gut virome of diarrheal children with Rotavirus A infection

Diarrhea is a leading cause of morbidity and mortality in children worldwide and represents a major dysbiosis event. Rotavirus has been recognized as a global leading pathogen of diarrhea. This study is aimed at investigating differences in the gut virome between diarrheal children and healthy controls. In 2018, 76 diarrheal fecal samples and 27 healthy fecal samples in Shanghai and 40 diarrheal fecal samples and 19 healthy fecal samples in Taizhou were collected to investigate the composition of the gut virome. Viral metagenomic analyses revealed that the alpha diversity of the diarrheal virome was not significantly different from that of the healthy virome, and the beta diversity had a significant difference between diarrheal and healthy children. The diarrheal virome was mainly dominated by the families Adenoviridae, Astroviridae, Caliciviridae, and Picornaviridae. Meanwhile, the healthy virome also contains phages, including Microviridae and Caudovirales. The high prevalence of diverse enteric viruses in all samples and the little abundance of Microviridae and Caudovirales in diarrheal groups were identified. The study introduced a general overview of the gut virome in diarrheal children, revealed the compositional differences in the gut viral community compared to healthy controls, and provided a reference for efficient treatments and prevention of virus-infectious diarrhea in children.

Novel Virus Identification through Metagenomics: A Systematic Review

Metagenomic Next Generation Sequencing (mNGS) allows the evaluation of complex microbial communities, avoiding isolation and cultivation of each microbial species, and does not require prior knowledge of the microbial sequences present in the sample. Applications of mNGS include virome characterization, new virus discovery and full-length viral genome reconstruction, either from virus preparations enriched in culture or directly from clinical and environmental specimens. Here, we systematically reviewed studies that describe novel virus identification through mNGS from samples of different origin (plant, animal and environment). Without imposing time limits to the search, 379 publications were identified that met the search parameters. Sample types, geographical origin, enrichment and nucleic acid extraction methods, sequencing platforms, bioinformatic analytical steps and identified viral families were described. The review highlights mNGS as a feasible method for novel virus discovery from samples of different origins, describes which kind of heterogeneous experimental and analytical protocols are currently used and provides useful information such as the different commercial kits used for the purification of nucleic acids and bioinformatics analytical pipelines.

Expanding known viral diversity in plants: virome of 161 species alongside an ancient canal

BACKGROUND

Since viral metagenomic approach was applied to discover plant viruses for the first time in 2006, many plant viruses had been identified from cultivated and non-cultivated plants. These previous researches exposed that the viral communities (virome) of plants have still largely uncharacterized. Here, we investigated the virome in 161 species belonging to 38 plant orders found in a riverside ecosystem.

RESULTS

We identified 245 distinct plant-associated virus genomes (88 DNA and 157 RNA viruses) belonging to 27 known viral families, orders, or unclassified virus groups. Some viral genomes were sufficiently divergent to comprise new species, genera, families, or even orders. Some groups of viruses were detected that currently are only known to infect organisms other than plants. It indicates a wider host range for members of these clades than previously recognized theoretically. We cannot rule out that some viruses could be from plant contaminating organisms, although some methods were taken to get rid of them as much as possible. The same viral species could be found in different plants and co-infections were common.

CONCLUSIONS

Our data describe a complex viral community within a single plant ecosystem and expand our understanding of plant-associated viral diversity and their possible host ranges.

Organizing the Global Diversity of Microviruses

Microviruses encompass an astonishing array of small, single-stranded DNA phages that, due to the surge in metagenomic surveys, are now known to be prevalent in most environments. Current taxonomy concedes the considerable diversity within this lineage to a single family (the Microviridae), which has rendered it difficult to adequately and accurately assess the amount of variation that actually exists within this group. We amassed and curated the largest collection of microviral genomes to date and, through a combination of protein-sharing networks and phylogenetic analysis, discovered at least three meaningful taxonomic levels between the current ranks of family and genus. When considering more than 13,000 microviral genomes from recognized lineages and as-yet-unclassified microviruses in metagenomic samples, microviral diversity is better understood by elevating microviruses to the level of an order that consists of three suborders and at least 19 putative families, each with their respective subfamilies. These revisions enable fine-scale assessment of microviral dynamics: for example, in the human gut, there are considerable differences in the abundances of microviral families both between urban and rural populations and in individuals over time. In addition, our analysis of genome contents and gene exchange shows that microviral families carry no recognizable accessory metabolic genes and rarely, if ever, engage in horizontal gene transfer across microviral families or with their bacterial hosts. These insights bring microviral taxonomy in line with current developments in the taxonomy of other phages and increase the understanding of microvirus biology.

Virome in the cloaca of wild and breeding birds revealed a diversity of significant viruses

BACKGROUND

Wild birds may harbor and transmit viruses that are potentially pathogenic to humans, domestic animals, and other wildlife.

RESULTS

Using the viral metagenomic approach, we investigated the virome of cloacal swab specimens collected from 3182 birds (the majority of them wild species) consisting of > 87 different species in 10 different orders within the Aves classes. The virus diversity in wild birds was higher than that in breeding birds. We acquired 707 viral genomes from 18 defined families and 4 unclassified virus groups, with 265 virus genomes sharing < 60% protein sequence identities with their best matches in GenBank comprising new virus families, genera, or species. RNA viruses containing the conserved RdRp domain with no phylogenetic affinity to currently defined virus families existed in different bird species. Genomes of the astrovirus, picornavirus, coronavirus, calicivirus, parvovirus, circovirus, retrovirus, and adenovirus families which include known avian pathogens were fully characterized. Putative cross-species transmissions were observed with viruses in wild birds showing > 95% amino acid sequence identity to previously reported viruses in domestic poultry. Genomic recombination was observed for some genomes showing discordant phylogenies based on structural and non-structural regions. Mapping the next-generation sequencing (NGS) data respectively against the 707 genomes revealed that these viruses showed distribution pattern differences among birds with different habitats (breeding or wild), orders, and sampling sites but no significant differences between birds with different behavioral features (migratory and resident).

CONCLUSIONS

The existence of a highly diverse virome highlights the challenges in elucidating the evolution, etiology, and ecology of viruses in wild birds. Video Abstract.

Complete genome analysis of the newly isolated Vibrio phage vB_VpP_WS1 of the family Microviridae

Here, we describe the characterization and genome annotation of the newly isolated lytic Vibrio parahaemolyticus phage vB_VpP_WS1, isolated from sewage samples collected in Qingdao, China. Transmission electron microscopy revealed that vB_VpP_WS1 is about 22 nm in size and that the virions are isometric, likely icosahedral, particles similar to those of members of the Microviridae. The digestion patterns of phage nucleic acids and whole-genome sequencing analysis together revealed that phage vB_VpP_WS1 has a single-stranded DNA genome of 5564 nt. Eight open reading frames were identified, only four of which could be annotated. The proteins of vB_VpP_WS1 displayed low sequence similarity to their homologs encoded by other microviruses. Phylogenetic analysis based on the major capsid protein suggested that vB_VpP_WS1 is a tentative new member of the family Microviridae.

Gut Virome of the World's Highest-Elevation Lizard Species (Phrynocephalus erythrurus and Phrynocephalus theobaldi) Reveals Versatile Commensal Viruses

The gut virome is a reservoir of diverse symbiotic and pathogenic viruses coevolving with their hosts, and yet limited research has explored the gut viromes of highland-dwelling rare species. Using viral metagenomic analysis, the viral communities of the Phrynocephalus lizards living in the Qinghai-Tibet Plateau were investigated. Phage-encoded functional genes and antibiotic resistance genes (ARGs) were analyzed. The viral communities of different lizard species were all predominated by bacteriophages, especially the Caudovirales order. The virome of Phrynocephalus erythrurus living around the Namtso Lake possessed a unique structure, with the greatest abundance of the Parvoviridae family and the highest number of exclusive viral species. Several vertebrate-infecting viruses were discovered, including caliciviruses, astroviruses, and parvoviruses. Phylogenetic analyses demonstrated that the virus hallmark genes of bacteriophages possessed high genetic diversity. After functional annotation, the majority of phage-associated functional genes were classified in the energy metabolism category. In addition, plenty of ARGs belonging to the multidrug category were discovered, and five ARGs were exclusive to the virome from Phrynocephalus theobaldi. This study provided the first insight into the structure and function of the virome in highland lizards, contributing to the protection of threatened lizard species. Also, our research is of exemplary significance for the gut virome research of lizard species and other cold-blooded and highland-dwelling animals, prompting a better understanding of the interspecific differences and transmission of commensal viruses. IMPORTANCE The Phrynocephalus lizards inhabiting the Qinghai-Tibet Plateau (QTP) are considered to be the highest-altitude lizard species in the world, and they have been added to the IUCN list of threatened species. Living in the QTP with hypoxic, arid, natural conditions, the lizards presented a unique pattern of gut virome, which could provide both positive and negative effects, such as the enrichment of functional genes and the dissemination of antibiotic resistance genes (ARGs). This work provides the foundation for further research on the gut virome in these endangered lizard species and other cold-blooded and highland-dwelling animals, contributing to the maintenance of ecological balance on the plateau.

Comparison of Gut Viral Communities in Atopic Dermatitis and Healthy Children

The immune development and regulation of living individuals are affected by the gut microbiota. The imbalance of gut microbiota is considered to be a key factor that easily induces immune dysregulation and the development of atopic diseases. Atopic dermatitis (AD) is a chronic inflammatory skin disease that affects nearly 20% of children. To date, metagenomics research on AD has mainly focused on the skin and gut microbiome. However, here we assessed the composition of the virome in the gut of AD patients and healthy controls for the first time. This study has obtained possible dominant viruses at different viral classification levels. In terms of diversity, the alpha diversity of the patients group was significantly lower than that of the healthy controls group, and the beta diversity of the two groups was significantly different from phylum to family level. These findings provide a new perspective for us to better understand the effect of the gut microecological environment on AD.

Virome in Fecal Samples From Wild Giant Pandas (Ailuropoda Melanoleuca)

The giant panda (Ailuropoda melanoleuca) is one of the most endangered mammals in the world; anthropogenic habitat loss and poaching still threaten the survival of wild pandas. Viral infection has become one of the potential threats to the health of these animals, but the available information related to these infections is still limited. In order to detect possible vertebrate viruses, the virome in the fecal samples of seven wild giant pandas from Qinling Mountains was investigated by using the method of viral metagenomics. From the fecal virome of wild giant pandas, we determined six nearly complete genomes belonging to the order Picornavirales, two of which may be qualified as a novel virus family or genus. In addition, four complete genomes belonging to the Genomoviridae family were also fully characterized. This virological investigation has increased our understanding of the gut viral community in giant pandas. Whether these viruses detected in fecal samples can really infect giant panda needs further research.

Comparison of gut viral communities in diarrhoea and healthy dairy calves

Calf diarrhoea has been a major cause of economic losses in the global dairy industry. Many factors, including multiple pathogen infections, can directly or indirectly cause calf diarrhoea. This study compared the faecal virome between 15 healthy calves and 15 calves with diarrhoea. Significantly lower diversity of viruses was found in samples from animals with diarrhoea than those in the healthy ones, and this feature may also be related to the age of the calves. Viruses belonging to the families Astroviridae and Caliciviridae that may cause diarrhoea in dairy calves have been characterized, which revealed that reads of caliciviruses and astroviruses in diarrhoea calves were much higher than those in healthy calves. Five complete genomic sequences closely related to Smacoviridae have been identified, which may participate in the regulation of the gut virus community ecology of healthy hosts together with bacteriophages. This research provides a theoretical basis for further understanding of known or potential enteric pathogens related to calf diarrhoea.

Adaptation of the polony technique to quantify Gokushovirinae, a diverse group of single-stranded DNA phage

Advances in metagenomics have revealed the ubiquity of single-stranded DNA (ssDNA) phage belonging to the subfamily Gokushovirinae in the oceans; however, the abundance and ecological roles of this group are unknown. Here, we quantify gokushoviruses through adaptation of the polony method, in which viral template DNA is immobilized in a gel, amplified by PCR, and subsequently detected by hybridization. Primers and probes for this assay were designed based on PCR amplicon diversity of gokushovirus major capsid protein gene sequences from a depth profile in the Gulf of Aqaba, Red Sea sampled in September 2015. At ≥95% identity, these 87 gokushovirus sequences formed 14 discrete clusters with the largest clades showing distinct depth distributions. The application of the polony method enabled the first quantification of gokushoviruses in any environment. The gokushoviruses were most abundant in the upper 40 m of the stratified water column, with a subsurface peak in abundance of 1.26 × 10⁵ viruses ml^-1 . These findings suggest that discrete gokushovirus genotypes infect bacterial hosts that differentially partition in the water column. Since the designed primers and probe are conserved across marine ecosystems, this polony method can be applied broadly for the quantification of gokushoviruses throughout the global oceans.

Legionella pneumophila CRISPR-Cas Suggests Recurrent Encounters with One or More Phages in the Family Microviridae

Legionella pneumophila is a ubiquitous freshwater pathogen and the causative agent of Legionnaires' disease. L. pneumophila growth within protists provides a refuge from desiccation, disinfection, and other remediation strategies. One outstanding question has been whether this protection extends to phages. L. pneumophila isolates are remarkably devoid of prophages and to date no Legionella phages have been identified. Nevertheless, many L. pneumophila isolates maintain active CRISPR-Cas defenses. So far, the only known target of these systems is an episomal element that we previously named Legionella mobile element 1 (LME-1). The continued expansion of publicly available genomic data promises to further our understanding of the role of these systems. We now describe over 150 CRISPR-Cas systems across 600 isolates to establish the clearest picture yet of L. pneumophila's adaptive defenses. By searching for targets of 1,500 unique CRISPR-Cas spacers, LME-1 remains the only identified CRISPR-Cas targeted integrative element. We identified 3 additional LME-1 variants-all targeted by previously and newly identified CRISPR-Cas spacers-but no other similar elements. Notably, we also identified several spacers with significant sequence similarity to microviruses, specifically those within the subfamily Gokushovirinae. These spacers are found across several different CRISPR-Cas arrays isolated from geographically diverse isolates, indicating recurrent encounters with these phages. Our analysis of the extended Legionella CRISPR-Cas spacer catalog leads to two main conclusions: current data argue against CRISPR-Cas targeted integrative elements beyond LME-1, and the heretofore unknown L. pneumophila phages are most likely lytic gokushoviruses. IMPORTANCE Legionnaires' disease is an often-fatal pneumonia caused by Legionella pneumophila, which normally grows inside amoebae and other freshwater protists. L. pneumophila trades diminished access to nutrients for the protection and isolation provided by the host. One outstanding question is whether L. pneumophila is susceptible to phages, given the protection provided by its intracellular lifestyle. In this work, we use Legionella CRISPR spacer sequences as a record of phage infection to predict that the "missing" L. pneumophila phages belong to the microvirus subfamily Gokushovirinae. Gokushoviruses are known to infect another intracellular pathogen, Chlamydia. How do gokushoviruses access L. pneumophila (and Chlamydia) inside their "cozy niches"? Does exposure to phages happen during a transient extracellular period (during cell-to-cell spread) or is it indicative of a more complicated environmental lifestyle? One thing is clear, 100 years after their discovery, phages continue to hold important secrets about the bacteria upon which they prey.

Fishing for the Virome of Tropical Tuna

While planktonic viruses have received much attention in recent decades, knowledge of the virome of marine organisms, especially fish, still remains rudimentary. This is notably the case with tuna, which are among the most consumed fish worldwide and represent considerable economic, social and nutritional value. Yet the composition of the tuna virome and its biological and environmental determinants remain unknown. To begin to address this gap, we investigated the taxonomic diversity of viral communities inhabiting the skin mucus, gut and liver of two major tropical tuna species (skipjack and yellowfin) in individuals fished in the Atlantic and Indian Oceans. While we found significant differences in the virome composition between the organs, this was totally independent of the tuna species or sex. The tuna virome was mainly dominated by eukaryotic viruses in the digestive organs (gut and liver), while bacteriophages were predominant in the mucus. We observed the presence of specific viral families in each organ, some previously identified as fish or human pathogens (e.g., Iridoviridae, Parvoviridae, Alloherpesviridae, Papillomaviridae). Interestingly, we also detected a ‘core virome’ that was shared by all the organs and was mainly composed of Caudovirales, Microviridae and Circoviridae. These results show that tuna host a mosaic of viral niches, whose establishment, role and circulation remain to be elucidated.

Virome of respiratory secretion from children with unknown etiological acute respiratory disease revealed recombinant human parechovirus and other significant viruses

Using viral metagenomics, viral nucleic acid in 30 respiratory secretion samples collected from children with unknown etiological acute respiratory disease were investigated. Sequences showing similarity to human parainfluenza virus 1, anellovirus, bocavirus, coxsackievirus A4, human parechovirus (HPeV), and alphaflexivirus were recovered from these samples. Complete genomes of one anellovirus, one coxsackievirus A4, three parechoviruses were determined from these libraries. The anellovirus (MW267851) phylogenetically clustered with an unpublished anellovirus (MK212032) from respiratory sample of a Vietnamese patient, forming a separate branch neighboring to strains within the genus Betatorquevirus. The genome of coxsackievirus A4 (MW267852) shares the highest sequence identity of 96.4% to a coxsackievirus A4 (MN964079) which was identified in clinical samples from children with Hand, Foot, and Mouth Disease (HFMD). Two (MW267853 and MW267854) of the three parechoviruses belong to HPeV-1 and the other one (MW267855) belongs to HPeV-6. Recombination analysis indicated that an HPeV-1 (MW267854) identified in this study is a putative recombinant occurred between HPeV-1 and HPeV-3. Whether these viruses have association with specific respiratory disease calls for further investigation.

Genome Sequences of Microviruses Identified in a Sample from a Sewage Treatment Oxidation Pond

Oxidation ponds are often used in the treatment of sewage as an aeration step prior to discharge. We identified 99 microvirus genomes from a sample from a sewage oxidation pond. This diverse group of microviruses expands our knowledge of bacteriophages associated with sewage oxidation pond ecosystems.

Oxidation ponds are often used in the treatment of sewage as an aeration step prior to discharge. We identified 99 microvirus genomes from a sample from a sewage oxidation pond. This diverse group of microviruses expands our knowledge of bacteriophages associated with sewage oxidation pond ecosystems.

Thousands of previously unknown phages discovered in whole-community human gut metagenomes

BACKGROUND

Double-stranded DNA bacteriophages (dsDNA phages) play pivotal roles in structuring human gut microbiomes; yet, the gut virome is far from being fully characterized, and additional groups of phages, including highly abundant ones, continue to be discovered by metagenome mining. A multilevel framework for taxonomic classification of viruses was recently adopted, facilitating the classification of phages into evolutionary informative taxonomic units based on hallmark genes. Together with advanced approaches for sequence assembly and powerful methods of sequence analysis, this revised framework offers the opportunity to discover and classify unknown phage taxa in the human gut.

RESULTS

A search of human gut metagenomes for circular contigs encoding phage hallmark genes resulted in the identification of 3738 apparently complete phage genomes that represent 451 putative genera. Several of these phage genera are only distantly related to previously identified phages and are likely to found new families. Two of the candidate families, "Flandersviridae" and "Quimbyviridae", include some of the most common and abundant members of the human gut virome that infect Bacteroides, Parabacteroides, and Prevotella. The third proposed family, "Gratiaviridae," consists of less abundant phages that are distantly related to the families Autographiviridae, Drexlerviridae, and Chaseviridae. Analysis of CRISPR spacers indicates that phages of all three putative families infect bacteria of the phylum Bacteroidetes. Comparative genomic analysis of the three candidate phage families revealed features without precedent in phage genomes. Some "Quimbyviridae" phages possess Diversity-Generating Retroelements (DGRs) that generate hypervariable target genes nested within defense-related genes, whereas the previously known targets of phage-encoded DGRs are structural genes. Several "Flandersviridae" phages encode enzymes of the isoprenoid pathway, a lipid biosynthesis pathway that so far has not been known to be manipulated by phages. The "Gratiaviridae" phages encode a HipA-family protein kinase and glycosyltransferase, suggesting these phages modify the host cell wall, preventing superinfection by other phages. Hundreds of phages in these three and other families are shown to encode catalases and iron-sequestering enzymes that can be predicted to enhance cellular tolerance to reactive oxygen species.

CONCLUSIONS

Analysis of phage genomes identified in whole-community human gut metagenomes resulted in the delineation of at least three new candidate families of Caudovirales and revealed diverse putative mechanisms underlying phage-host interactions in the human gut. Addition of these phylogenetically classified, diverse, and distinct phages to public databases will facilitate taxonomic decomposition and functional characterization of human gut viromes. Video abstract.

The healthy human virome: from virus-host symbiosis to disease

Viruses are ubiquitous, essential components of any ecosystem, and of multicellular organism holobionts. Numerous viruses cause acute infection, killing the host or being cleared by immune system. In many other cases, viruses coexist with the host as symbionts, either temporarily or for the duration of the host's life. Apparently, virus-host relationships span the entire range from aggressive parasitism to mutualism. Here we attempt to delineate the healthy human virome, that is, the entirety of viruses that are present in a healthy human body. The bulk of the healthy virome consists of bacteriophages infecting bacteria in the intestine and other locations. However, a variety of viruses, such as anelloviruses and herpesviruses, and the numerous endogenous retroviruses, persist by replicating in human cells, and these are our primary focus. Crucially, the boundary between symbiotic and pathogenic viruses is fluid such that members of the healthy virome can become pathogens under changing conditions.

The ecogenomics of dsDNA bacteriophages in feces of stabled and feral horses

The viromes of the mammalian lower gut were shown to be heavily dominated by bacteriophages; however, only for humans were the composition and intervariability of the bacteriophage communities studied in depth. Here we present an ecogenomics survey of dsDNA bacteriophage diversity in the feces of horses (Equus caballus), comparing two groups of stabled horses, and a further group of feral horses that were isolated on an island. Our results indicate that the dsDNA viromes of the horse feces feature higher richness than in human viromes, with more even distribution of genotypes. No over-represented phage genotypes, such as CrAssphage-related viruses found in humans, were identified. Additionally, many bacteriophage genus-level clusters were found to be present in all three geographically isolated populations. The diversity of the horse intestinal bacteriophages is severely undersampled, and so consequently only a minor fraction of the phage contigs could be linked with the bacteriophage genomes. Our study indicates that bacteriophage ecological parameters in the intestinal ecosystems in horses and humans differ significantly, leading them to shape their corresponding viromes in different ways. Therefore, the diversity and structure of the intestinal virome in different animal species needs to be experimentally studied.

Single-gene lysis in the metagenomic era

The small lytic phages (Microviridae and Leviviridae), effect bacterial lysis with the product of a single gene. The three well-studied single-gene lysis (Sgl) proteins (E of φX174, A₂ of Qβ, and Lys^M of phage M) lack direct muralytic activity, and have been shown to function as 'protein antibiotics' by acting as noncompetitive inhibitors of conserved peptidoglycan (PG) biosynthesis enzymes, MurA, MraY, and MurJ respectively. The fourth, protein L of MS2, does not inhibit PG biosynthesis but instead is hypothesized to trigger host autolytic response through an unknown mechanism. Recent advances in meta-omics approaches have led to an explosion in the available genomes of small lytic phages. Of the thousands of new genomes, only one annotated Sgl shared some sequence similarity with a known Sgl (L of MS2), highlighting the diversity in Sgls. The newly available genomic space serves as an untapped resource for discovering novel Sgls.

Exploring the Remarkable Diversity of Culturable Escherichia coli Phages in the Danish Wastewater Environment

Phages drive bacterial diversity, profoundly influencing microbial communities, from microbiomes to the drivers of global biogeochemical cycling. Aiming to broaden our understanding of Escherichiacoli (MG1655, K-12) phages, we screened 188 Danish wastewater samples and isolated 136 phages. Ninety-two of these have genomic sequences with less than 95% similarity to known phages, while most map to existing genera several represent novel lineages. The isolated phages are highly diverse, estimated to represent roughly one-third of the true diversity of culturable virulent dsDNA Escherichia phages in Danish wastewater, yet almost half (40%) are not represented in metagenomic databases, emphasising the importance of isolating phages to uncover diversity. Seven viral families, Myoviridae, Siphoviridae, Podoviridae,Drexlerviridae,Chaseviridae,Autographviridae, and Microviridae, are represented in the dataset. Their genomes vary drastically in length from 5.3 kb to 170.8 kb, with a guanine and cytosine (GC) content ranging from 35.3% to 60.0%. Hence, even for a model host bacterium, substantial diversity remains to be uncovered. These results expand and underline the range of coliphage diversity and demonstrate how far we are from fully disclosing phage diversity and ecology.

Taxonomy of prokaryotic viruses: 2018-2019 update from the ICTV Bacterial and Archaeal Viruses Subcommittee

This article is a summary of the activities of the ICTV's Bacterial and Archaeal Viruses Subcommittee for the years 2018 and 2019. Highlights include the creation of a new order, 10 families, 22 subfamilies, 424 genera and 964 species. Some of our concerns about the ICTV's ability to adjust to and incorporate new DNA- and protein-based taxonomic tools are discussed.

Resurrection of a global, metagenomically defined gokushovirus

Gokushoviruses are single-stranded, circular DNA bacteriophages found in metagenomic datasets from diverse ecosystems worldwide, including human gut microbiomes. Despite their ubiquity and abundance, little is known about their biology or host range: Isolates are exceedingly rare, known only from three obligate intracellular bacterial genera. By synthesizing circularized phage genomes from prophages embedded in diverse enteric bacteria, we produced gokushoviruses in an experimentally tractable model system, allowing us to investigate their features and biology. We demonstrate that virions can reliably infect and lysogenize hosts by hijacking a conserved chromosome-dimer resolution system. Sequence motifs required for lysogeny are detectable in other metagenomically defined gokushoviruses; however, we show that even partial motifs enable phages to persist cytoplasmically without leading to collapse of their host culture. This ability to employ multiple, disparate survival strategies is likely key to the long-term persistence and global distribution of Gokushovirinae.

Single-Stranded DNA Viruses in Antarctic Cryoconite Holes

Antarctic cryoconite holes, or small melt-holes in the surfaces of glaciers, create habitable oases for isolated microbial communities with tightly linked microbial population structures. Viruses may influence the dynamics of polar microbial communities, but the viromes of the Antarctic cryoconite holes have yet to be characterized. We characterize single-stranded DNA (ssDNA) viruses from three cryoconite holes in the Taylor Valley, Antarctica, using metagenomics. Half of the assembled metagenomes cluster with those in the viral family Microviridae (n = 7), and the rest with unclassified circular replication associated protein (Rep)-encoding single-stranded (CRESS) DNA viruses (n = 7). An additional 18 virus-like circular molecules encoding either a Rep, a capsid protein gene, or other unidentified but viral-like open reading frames were identified. The samples from which the genomes were identified show a strong gradient in microbial diversity and abundances, and the number of viral genomes detected in each sample mirror that gradient. Additionally, one of the CRESS genomes assembled here shares ~90% genome-wide pairwise identity with a virus identified from a freshwater pond on the McMurdo Ice Shelf (Antarctica). Otherwise, the similarity of these viruses to those previously identified is relatively low. Together, these patterns are consistent with the presence of a unique regional virome present in fresh water host populations of the McMurdo Dry Valley region.

High Diversity and Novel Enteric Viruses in Fecal Viromes of Healthy Wild and Captive Thai Cynomolgus Macaques (Macaca fascicularis)

Cynomolgus macaques are common across South East Asian countries including Thailand. The National Primate Research Center of Thailand, Chulalongkorn University (NPRCT-CU) captures wild-borne cynomolgus macaque for research use. Limited information is available on the enteric viruses and possible zoonotic infections into or from cynomolgus macaques. We characterized and compare the fecal virome of two populations; healthy wild-originated captive cynomolgus macaques (n = 43) reared in NPRCT-CU and healthy wild cynomolgus macaques (n = 35). Over 90% of recognized viral sequence reads amplified from feces were from bacterial viruses. Viruses from seven families of mammalian viruses were also detected (Parvoviridae, Anelloviridae, Picornaviridae, Adenoviridae, Papillomaviridae, Herpesviridae, and Caliciviridae). The genomes of a member of a new picornavirus genus we named Mafapivirus, a primate chapparvovirus, and a circular Rep-encoding single-strand (CRESS) DNA virus were also characterized. Higher abundance of CRESS DNA viruses of unknown tropism and invertebrate-tropic ambidensovirus were detected in wild versus captive macaques likely reflecting dietary differences. Short term rearing in captivity did not have a pronounced effect on the diversity of mammalian viruses of wild cynomolgus macaques. This study is the first report of the fecal virome of cynomolgus macaques, non-human primates frequently used in biomedical research and vaccination studies.