Expansion of the global RNA virome reveals diverse clades of bacteriophages.

High-throughput RNA sequencing offers broad opportunities to explore the Earth RNA virome. Mining 5,150 diverse metatranscriptomes uncovered >2.5 million RNA virus contigs. Analysis of >330,000 RNA-dependent RNA polymerases (RdRPs) shows that this expansion corresponds to a 5-fold increase of the known RNA virus diversity. Gene content analysis revealed multiple protein domains previously not found in RNA viruses and implicated in virus-host interactions. Extended RdRP phylogeny supports the monophyly of the five established phyla and reveals two putative additional bacteriophage phyla and numerous putative additional classes and orders. The dramatically expanded phylum Lenarviricota, consisting of bacterial and related eukaryotic viruses, now accounts for a third of the RNA virome. Identification of CRISPR spacer matches and bacteriolytic proteins suggests that subsets of picobirnaviruses and partitiviruses, previously associated with eukaryotes, infect prokaryotic hosts.

Viruses direct carbon cycling in lake sediments under global change.

Global change is altering the vast amount of carbon cycled by microbes between land and freshwater, but how viruses mediate this process is poorly understood. Here, we show that viruses direct carbon cycling in lake sediments, and these impacts intensify with future changes in water clarity and terrestrial organic matter (tOM) inputs. Using experimental tOM gradients within sediments of a clear and a dark boreal lake, we identified 156 viral operational taxonomic units (vOTUs), of which 21% strongly increased with abundances of key bacteria and archaea, identified via metagenome-assembled genomes (MAGs). MAGs included the most abundant prokaryotes, which were themselves associated with dissolved organic matter (DOM) composition and greenhouse gas (GHG) concentrations. Increased abundances of virus-like particles were separately associated with reduced bacterial metabolism and with shifts in DOM toward amino sugars, likely released by cell lysis rather than higher molecular mass compounds accumulating from reduced tOM degradation. An additional 9.6% of vOTUs harbored auxiliary metabolic genes associated with DOM and GHGs. Taken together, these different effects on host dynamics and metabolism can explain why abundances of vOTUs rather than MAGs were better overall predictors of carbon cycling. Future increases in tOM quantity, but not quality, will change viral composition and function with consequences for DOM pools. Given their importance, viruses must now be explicitly considered in efforts to understand and predict the freshwater carbon cycle and its future under global environmental change.

Why are there so few (or so many) circulating coronaviruses?

Despite vast diversity in non-human hosts and conspicuous recent spillover events, only a small number of coronaviruses have been observed to persist in human populations. This puzzling mismatch suggests substantial barriers to establishment. We detail hypotheses that might contribute to explain the low numbers of endemic coronaviruses, despite their considerable evolutionary and emergence potential. We assess possible explanations ranging from issues of ascertainment, historically lower opportunities for spillover, aspects of human demographic changes, and features of pathogen biology and pre-existing adaptive immunity to related viruses. We describe how successful emergent viral species must triangulate transmission, virulence, and host immunity to maintain circulation. Characterizing the factors that might shape the limits of viral persistence can delineate promising research directions to better understand the combinations of pathogens and contexts that are most likely to lead to spillover.

Viral Infection Leads to a Unique Suite of Allelopathic Chemical Signals in Three Diatom Host-Virus Pairs.

Ecophysiological stress and the grazing of diatoms are known to elicit the production of chemical defense compounds called oxylipins, which are toxic to a wide range of marine organisms. Here we show that (1) the viral infection and lysis of diatoms resulted in oxylipin production; (2) the suite of compounds produced depended on the diatom host and the infecting virus; and (3) the virus-mediated oxylipidome was distinct, in both magnitude and diversity, from oxylipins produced due to stress associated with the growth phase. We used high-resolution accurate-mass mass spectrometry to observe changes in the dissolved lipidome of diatom cells infected with viruses over 3 to 4 days, compared to diatom cells in exponential, stationary, and decline phases of growth. Three host virus pairs were used as model systems: Chaetoceros tenuissimus infected with CtenDNAV; C. tenuissimus infected with CtenRNAV; and Chaetoceros socialis infected with CsfrRNAV. Several of the compounds that were significantly overproduced during viral infection are known to decrease the reproductive success of copepods and interfere with microzooplankton grazing. Specifically, oxylipins associated with allelopathy towards zooplankton from the 6-, 9-, 11-, and 15-lipogenase (LOX) pathways were significantly more abundant during viral lysis. 9-hydroperoxy hexadecatetraenoic acid was identified as the strongest biomarker for the infection of Chaetoceros diatoms. C. tenuissimus produced longer, more oxidized oxylipins when lysed by CtenRNAV compared to CtenDNAV. However, CtenDNAV caused a more statistically significant response in the lipidome, producing more oxylipins from known diatom LOX pathways than CtenRNAV. A smaller set of compounds was significantly more abundant in stationary and declining C. tenuissimus and C. socialis controls. Two allelopathic oxylipins in the 15-LOX pathway and essential fatty acids, arachidonic acid (ARA), eicosapentaenoic acid (EPA), and docosahexaenoic acid (DHA) were more abundant in the stationary phase than during the lysis of C. socialis. The host-virus pair comparisons underscore the species-level differences in oxylipin production and the value of screening more host-virus systems. We propose that the viral infection of diatoms elicits chemical defense via oxylipins which deters grazing with downstream trophic and biogeochemical effects.

Using Discarded Facial Tissues to Monitor and Diagnose Viral Respiratory Infections.

Molecular biology amplification enables sensitive detection of most respiratory viruses through nasopharyngeal swabbing. We developed an innovative approach to detect viral genomes on used facial tissues. In 2 communities of children, used tissues were collected once weekly for 1 year. Pooled analysis of tissues enabled detection of successive virus circulation in 4 age groups over time and forecasted by several weeks the circulation of influenza in the general population. At the individual level, in a proof-of-concept study of 30 volunteers with influenza-like signs/symptoms, we identified common respiratory viruses. The signals for SARS-CoV-2 obtained in parallel from 15 facial tissues and swab samples were similar and often higher for the tissues (11/15). Individual analysis of tissues offers a noninvasive, sensitive, and affordable alternative to self-sampling without a medical care requirement. Pooled analyses may be used to detect virus spread in specific communities, predict seasonal epidemics, and alert the population to viral infections.

Glycosylation and the global virome.

The sugars that coat the outsides of viruses and host cells are key to successful disease transmission, but they remain understudied compared to other molecular features. Understanding the comparative zoology of glycosylation - and harnessing it for predictive science - could help close the molecular gap in zoonotic risk assessment.

Analysis of Different Size Fractions Provides a More Complete Perspective of Viral Diversity in a Freshwater Embayment.

Inspired by recent discoveries of the prevalence of large viruses in the environment, we reassessed the longstanding approach of filtering water through small-pore-size filters to separate viruses from cells before metagenomic analysis. We collected samples from three sites in Hamilton Harbour, an embayment of Lake Ontario, and studied 6 data sets derived from <0.45-µm- and >0.45-µm-size fractions to compare the diversity of viruses in these fractions. At the level of virus order/family, we observed highly diverse and distinct virus communities in the >0.45-µm-size fractions, whereas the <0.45-µm-size fractions were composed primarily of Caudovirales The relative abundances of Caudovirales for which hosts could be inferred varied widely between size fractions, with higher relative abundances of cyanophages in the >0.45-µm-size fractions, potentially indicating replication within cells during ongoing infections. Many viruses of eukaryotes, such as Mimiviridae, Phycodnaviridae, Iridoviridae, and Poxviridae, were detected exclusively in the often-disregarded >0.45-µm-size fractions. In addition to observing unique virus communities associated with each size fraction from every site we examined, we detected viruses common to both fractions, suggesting that these are candidates for further exploration because they could be the product of ongoing or recent lytic events. Most importantly, our observations indicate that analysis of either fraction alone provides only a partial perspective of double-stranded DNA (dsDNA) viruses in the environment, highlighting the need for more comprehensive approaches for analyzing virus communities inferred from metagenomic sequencing.IMPORTANCE Most studies of aquatic virus communities analyze DNA sequences derived from the smaller-size "free-virus" fraction. Our study demonstrates that analysis of virus communities using only the smaller-size fraction can lead to erroneously low diversity estimates for many of the larger viruses such as Mimiviridae, Phycodnaviridae, Iridoviridae, and Poxviridae, whereas analyzing only the larger->0.45-µm-size fraction can lead to underestimates of Caudovirales diversity and relative abundance. Similarly, our data show that examining only the smaller-size fraction can lead to underestimations of virophage and cyanophage relative abundances that could, in turn, cause researchers to assume their limited ecological importance. Given the considerable differences we observed in this study, we recommend cautious interpretations of environmental virus community assemblages and dynamics when based on metagenomic data derived from different size fractions.

[A century of research on bacteriophages]. / Un siècle de recherche sur les bactériophages.

Bacteriophages have a prominent place in the living world. They participate to our understanding of the living world through three main aspects : (i) the dissection of the most intimist aspects of viral infection molecular mechanisms (molecular biology), (ii) the description and functioning mechanisms of ecosystems (ecology), and (iii) the adaptive dynamics of integrated viral and host-cell populations (evolution). This review looks back at the genesis of these fundamental findings and draws a picture of the most active fields of current research.

A Survey on Plant Viruses in Natural Brassicaceae Communities Using RNA-Seq.

Studies on plant viruses are biased towards crop diseases and little is known about viruses in natural vegetation. We conducted extensive surveys of plant viruses in wild Brassicaceae plants occurring in three local plant communities in central Japan. We applied RNA-Seq with selective depletion of rRNA, which allowed us to detect infections of all genome-reported viruses simultaneously. Infections of Turnip mosaic virus (TuMV), Cucumber mosaic virus (CMV), Brassica yellows virus, Pelargonium zonate spot virus, and Arabidopsis halleri partitivirus 1 were detected from the two perennial species, Arabidopsis halleri subsp. gemmifera and Rorippa indica. De novo assembly further detected partial sequences of a putative novel virus in Arabis fragellosa. Virus species composition and infection rate differed depending on site and plant species. Viruses were most frequently detected from the perennial clonal plant, A. halleri, in which a high clonal transmission rate of viruses across multiple years was confirmed. Phylogenetic analysis of TuMV and CMV showed that virus strains from wild Brassicaceae were included as a major clade of these viruses with other reported strains from crop plants, suggesting that viruses were shared among wild plants and crops. Our studies indicated that distribution of viruses in natural plant populations are determined by the combinations of life histories of viruses and hosts. Revealing viral distribution in the natural plant communities improves our knowledge on the ecology of plant viruses.

Assessing circovirus gene flow in multiple spill-over events.

The establishment of viral pathogens in new host environments following spillover events probably requires adaptive changes within both the new host and pathogen. After many generations, signals for ancient cross-species transmission may become lost and a strictly host-adapted phylogeny may mimic true co-divergence while the virus may retain an inherent ability to jump host species. The mechanistic basis for such processes remains poorly understood. To study the dynamics of virus-host co-divergence and the arbitrary chances of spillover in various reservoir hosts with equal ecological opportunity, we examined structural constraints of capsid protein in extant populations of Beak and feather disease virus (BFDV) during known spillover events. By assessing reservoir-based genotype stratification, we identified co-divergence defying signatures in the evolution BFDV which highlighted primordial processes of cryptic host adaptation and competing forces of host co-divergence and cross-species transmission. We demonstrate that, despite extensive surface plasticity gathered over a longer span of evolution, structural constraints of the capsid protein allow opportunistic host switching in host-adapted populations. This study provides new insights into how small populations of endangered psittacine species may face multidirectional forces of infection from reservoirs with apparently co-diverging genotypes.

Serologic and Molecular Evidence of Vaccinia Virus Circulation among Small Mammals from Different Biomes, Brazil.

Vaccinia virus (VACV) is a zoonotic agent that causes a disease called bovine vaccinia, which is detected mainly in milking cattle and humans in close contact with these animals. Even though many aspects of VACV infection have been described, much is still unknown about its circulation in the environment and its natural hosts/reservoirs. To investigate the presence of Orthopoxvirus antibodies or VACV DNA, we captured small rodents and marsupials in 3 areas of Minas Gerais state, Brazil, and tested their samples in a laboratory. A total of 336 animals were tested; positivity ranged from 18.1% to 25.5% in the 3 studied regions located in different biomes, including the Atlantic Forest and the Cerrado. Analysis of nucleotide sequences indicated co-circulation of VACV groups I and II. Our findings reinforce the possible role played by rodents and marsupials in VACV maintenance and its transmission chain.

Ribonucleotide reductases reveal novel viral diversity and predict biological and ecological features of unknown marine viruses.

Virioplankton play a crucial role in aquatic ecosystems as top-down regulators of bacterial populations and agents of horizontal gene transfer and nutrient cycling. However, the biology and ecology of virioplankton populations in the environment remain poorly understood. Ribonucleotide reductases (RNRs) are ancient enzymes that reduce ribonucleotides to deoxyribonucleotides and thus prime DNA synthesis. Composed of three classes according to O2 reactivity, RNRs can be predictive of the physiological conditions surrounding DNA synthesis. RNRs are universal among cellular life, common within viral genomes and virioplankton shotgun metagenomes (viromes), and estimated to occur within >90% of the dsDNA virioplankton sampled in this study. RNRs occur across diverse viral groups, including all three morphological families of tailed phages, making these genes attractive for studies of viral diversity. Differing patterns in virioplankton diversity were clear from RNRs sampled across a broad oceanic transect. The most abundant RNRs belonged to novel lineages of podoviruses infecting α-proteobacteria, a bacterial class critical to oceanic carbon cycling. RNR class was predictive of phage morphology among cyanophages and RNR distribution frequencies among cyanophages were largely consistent with the predictions of the "kill the winner-cost of resistance" model. RNRs were also identified for the first time to our knowledge within ssDNA viromes. These data indicate that RNR polymorphism provides a means of connecting the biological and ecological features of virioplankton populations.

Temperature effects on virion volume and genome length in dsDNA viruses.

Heterogeneity in rates of survival, growth and reproduction among viruses is related to virus particle (i.e. virion) size, but we have little understanding of the factors that govern the four to five orders of magnitude in virus size variation. Here, we analyse variation in virion size in 67 double-stranded DNA viruses (i.e. dsDNA) that span all major biomes, and infect organisms ranging from single-celled prokaryotes to multicellular eukaryotes. We find that two metrics of virion size (i.e. virion volume and genome length) decrease by about 55-fold as the temperature of occurrence increases from 0 to 40°C. We also find that gene overlap increases exponentially with temperature, such that smaller viruses have proportionally greater gene overlap at higher temperatures. These results indicate dsDNA virus size increases with environmental temperature in much the same way as the cell or genome size of many host species.

The virus of my virus is my friend: ecological effects of virophage with alternative modes of coinfection.

Virophages are viruses that rely on the replication machinery of other viruses to reproduce within eukaryotic hosts. Two different modes of coinfection have been posited based on experimental observation. In one mode, the virophage and the virus enter the host independently. In the other mode, the virophage adheres to the virus so both virophage and virus enter the host together. Here we ask: what are the ecological effects of these different modes of coinfection? In particular, what ecological effects are common to both infection modes, and what are the differences particular to each mode? We develop a pair of biophysically motivated ODE models of viral-host population dynamics, corresponding to dynamics arising from each mode of infection. We find that both modes of coinfection allow for the coexistence of the virophage, virus, and host either at a stable fixed point or through cyclical dynamics. In both models, virophage tends to be the most abundant population and their presence always reduces the viral abundance and increases the host abundance. However, we do find qualitative differences between models. For example, via extensive sampling of biologically relevant parameter space, we only observe bistability when the virophage and the virus enter the host together. We discuss how such differences may be leveraged to help identify modes of infection in natural environments from population level data.

Host-virus coevolution drives soil microbial function succession along a millennium land reclamation chronosequence.

INTRODUCTION: Gene exchange between viruses and hosts plays an important role in driving virus-host coevolution, enabling adaptation of both viruses and hosts to environmental changes. However, the mechanisms and functional significance of virus-host gene exchanges over long-term scales remain largely unexplored. OBJECTIVE: The present study aimed to gain insights into the role of viruses in virus-host interactions and coevolution by monitoring virome dynamics along a millennium-long land reclamation chronosequence. METHODS: We collected 24 soil samples from 8 stages of a millennium-long land reclamation chronosequence, including non-reclamation, and reclamation periods of 10, 50, 100, 300, 500, 700, and 1000 years. We characterized their metagenomes, and identified DNA viruses within these metagenomes. RESULTS: Our findings reveal a significant shift in viral community composition after 50 years of land reclamation, but soil viral diversity reached a stable phase approximately 300 years after the initial reclamation. Analysis of the virus-host network showed a scale-free degree distribution and a reduction in complexity over time, with generalist viruses emerging as key facilitators of horizontal gene transfer. CONCLUSION: These findings highlight the integral role of viruses, especially generalist types, in mediating gene exchanges between viruses and hosts, thereby influencing the coevolutionary dynamics in soil ecosystems over significant timescales. This study offers novel insights into long-term virus-host interactions, showing how the virome responds to environmental changes, driving shifts in various microbial functions in reclaimed land.

The diverse liver viromes of Australian geckos and skinks are dominated by hepaciviruses and picornaviruses and reflect host taxonomy and habitat.

Lizards have diverse ecologies and evolutionary histories, and represent a promising group to explore how hosts shape virome structure and virus evolution. Yet, little is known about the viromes of these animals. In Australia, squamates (lizards and snakes) comprise the most diverse order of vertebrates, and Australia hosts the highest diversity of lizards globally, with the greatest breadth of habitat use. We used meta-transcriptomic sequencing to determine the virome of nine co-distributed, tropical lizard species from three taxonomic families in Australia and analyzed these data to identify host traits associated with viral abundance and diversity. We show that lizards carry a large diversity of viruses, identifying more than thirty novel, highly divergent vertebrate-associated viruses. These viruses were from nine viral families, including several that contain well known pathogens, such as the Flaviviridae, Picornaviridae, Bornaviridae, Iridoviridae, and Rhabdoviridae. Members of the Flaviviridae were particularly abundant across species sampled here, largely belonging to the genus Hepacivirus: fourteen novel hepaciviruses were identified, broadening the known diversity of this group and better defining its evolution by uncovering new reptilian clades. The evolutionary histories of the viruses studied here frequently aligned with the biogeographic and phylogenetic histories of the hosts, indicating that exogenous viruses may help infer host evolutionary history if sampling is strategic and sampling density high enough. Notably, analysis of alpha and beta diversity revealed that virome composition and richness in the animals sampled here was shaped by host taxonomy and habitat. In sum, we identified a diverse range of reptile viruses that broadly contributes to our understanding of virus-host ecology and evolution.

Rapid taxonomic categorization of short, abundant virus sequences for ecological analyses.

Public health concerns about recent viral epidemics have motivated researchers to seek novel ways to understand pathogen infection in native, wildlife hosts. With its deep history of tools and perspectives for understanding the abundance and distribution of organisms, ecology can shed new light on viral infection dynamics. However, datasets allowing deep explorations of viral communities from an ecological perspective are lacking. We sampled 1086 bats from two, adjacent Puerto Rican caves and tested them for infection by herpesviruses, resulting in 3131 short, viral sequences. Using percent identity of nucleotides and a machine learning algorithm (affinity propagation), we categorized herpesviruses into 43 operational taxonomic units (OTUs) to be used in place of species in subsequent ecological analyses. Herpesvirus metacommunities demonstrated long-tailed rank frequency distributions at all analyzed levels of host organization (i.e., individual, population, and community). Although 13 herpesvirus OTUs were detected in more than one host species, OTUs generally exhibited host specificity by infecting a single core host species at a significantly higher prevalence than in all satellite species combined. We describe the natural history of herpesvirus metacommunities in Puerto Rican bats and suggest that viruses follow the general law that communities comprise few common and many rare species. To guide future efforts in the field of viral ecology, hypotheses are presented regarding mechanisms that contribute to these patterns.

Targeted viromes and total metagenomes capture distinct components of bee gut phage communities.

BACKGROUND: Despite being among the most abundant biological entities on earth, bacteriophage (phage) remain an understudied component of host-associated systems. One limitation to studying host-associated phage is the lack of consensus on methods for sampling phage communities. Here, we compare paired total metagenomes and viral size fraction metagenomes (viromes) as methods for investigating the dsDNA viral communities associated with the GI tract of two bee species: the European honey bee Apis mellifera and the eastern bumble bee Bombus impatiens. RESULTS: We find that viromes successfully enriched for phage, thereby increasing phage recovery, but only in honey bees. In contrast, for bumble bees, total metagenomes recovered greater phage diversity. Across both bee species, viromes better sampled low occupancy phage, while total metagenomes were biased towards sampling temperate phage. Additionally, many of the phage captured by total metagenomes were absent altogether from viromes. Comparing between bees, we show that phage communities in commercially reared bumble bees are significantly reduced in diversity compared to honey bees, likely reflecting differences in bacterial titer and diversity. In a broader context, these results highlight the complementary nature of total metagenomes and targeted viromes, especially when applied to host-associated environments. CONCLUSIONS: Overall, we suggest that studies interested in assessing total communities of host-associated phage should consider using both approaches. However, given the constraints of virome sampling, total metagenomes may serve to sample phage communities with the understanding that they will preferentially sample dominant and temperate phage. Video Abstract.

Host species and geography impact bee-associated RNA virus communities with evidence for isolation by distance in viral populations.

Virus symbionts are important mediators of ecosystem function, yet we know little of their diversity and ecology in natural populations. The alarming decline of pollinating insects in many regions of the globe, especially the European honey bee, Apis mellifera, has been driven in part by worldwide transmission of virus pathogens. Previous work has examined the transmission of known honey bee virus pathogens to wild bee populations, but only a handful of studies have investigated the native viromes associated with wild bees, limiting epidemiological predictors associated with viral pathogenesis. Further, variation among different bee species might have important consequences in the acquisition and maintenance of bee-associated virome diversity. We utilized comparative metatranscriptomics to develop a baseline description of the RNA viromes associated with wild bee pollinators and to document viral diversity, community composition, and structure. Our sampling includes five wild-caught, native bee species that vary in social behavior as well as managed honey bees. We describe 26 putatively new RNA virus species based on RNA-dependent RNA polymerase phylogeny and show that each sampled bee species was associated with a specific virus community composition, even among sympatric populations of distinct host species. From 17 samples of a single host species, we recovered a single virus species despite over 600 km of distance between host populations and found strong evidence for isolation by distance in associated viral populations. Our work adds to the small number of studies examining viral prevalence and community composition in wild bees.

Virome release of an invasive exotic plant species in southern France.

The increase in human-mediated introduction of plant species to new regions has resulted in a rise of invasive exotic plant species (IEPS) that has had significant effects on biodiversity and ecosystem processes. One commonly accepted mechanism of invasions is that proposed by the enemy release hypothesis (ERH), which states that IEPS free from their native herbivores and natural enemies in new environments can outcompete indigenous species and become invasive. We here propose the virome release hypothesis (VRH) as a virus-centered variant of the conventional ERH that is only focused on enemies. The VRH predicts that vertically transmitted plant-associated viruses (PAV, encompassing phytoviruses and mycoviruses) should be co-introduced during the dissemination of the IEPS, while horizontally transmitted PAV of IEPS should be left behind or should not be locally transmitted in the introduced area due to a maladaptation of local vectors. To document the VRH, virome richness and composition as well as PAV prevalence, co-infection, host range, and transmission modes were compared between indigenous plant species and an invasive grass, cane bluestem (Bothriochloa barbinodis), in both its introduced range (southern France) and one area of its native range (Sonoran Desert, Arizona, USA). Contrary to the VRH, we show that invasive populations of B. barbinodis in France were not associated with a lower PAV prevalence or richness than native populations of B. barbinodis from the USA. However, comparison of virome compositions and network analyses further revealed more diverse and complex plant-virus interactions in the French ecosystem, with a significant richness of mycoviruses. Setting mycoviruses apart, only one putatively vertically transmitted phytovirus (belonging to the Amalgaviridae family) and one putatively horizontally transmitted phytovirus (belonging to the Geminiviridae family) were identified from B. barbinodis plants in the introduced area. Collectively, these characteristics of the B. barbinodis-associated PAV community in southern France suggest that a virome release phase may have immediately followed the introduction of B. barbinodis to France in the 1960s or 1970s, and that, since then, the invasive populations of this IEPS have already transitioned out of this virome release phase, and have started interacting with several local mycoviruses and a few local plant viruses.