Signatures of natural selection in tree topology shape of serially sampled viral phylogenies.

Tree shape metrics can be computed fast for trees of any size, which makes them promising alternatives to intensive statistical methods and parameter-rich evolutionary models in the era of massive data availability. Previous studies have demonstrated their effectiveness in unveiling important parameters in viral evolutionary dynamics, although the impact of natural selection on the shape of tree topologies has not been thoroughly investigated. We carried out a forward-time and individual-based simulation to investigate whether tree shape metrics of several kinds could predict the selection regime employed to generate the data. To examine the impact of the genetic diversity of the founder viral population, simulations were run under two opposing starting configurations of the genetic diversity of the infecting viral population. We found that four evolutionary regimes, namely, negative, positive, and frequency-dependent selection, as well as neutral evolution, were successfully distinguished by tree topology shape metrics. Two metrics from the Laplacian spectral density profile (principal eigenvalue and peakedness) and the number of cherries were the most informative for indicating selection type. The genetic diversity of the founder population had an impact on differentiating evolutionary scenarios. Tree imbalance, which has been frequently associated with the action of natural selection on intrahost viral diversity, was also characteristic of neutrally evolving serially sampled data. Metrics calculated from empirical analysis of HIV datasets indicated that most tree topologies exhibited shapes closer to the frequency-dependent selection or neutral evolution regimes.

Diversity and Pathobiology of an Ilarvirus Unexpectedly Detected in Diverse Plants and Global Sequencing Data.

High-throughput sequencing (HTS) and sequence mining tools revolutionized virus detection and discovery in recent years, and implementing them with classical plant virology techniques results in a powerful approach to characterize viruses. An example of a virus discovered through HTS is Solanum nigrum ilarvirus 1 (SnIV1) (Bromoviridae), which was recently reported in various solanaceous plants from France, Slovenia, Greece, and South Africa. It was likewise detected in grapevines (Vitaceae) and several Fabaceae and Rosaceae plant species. Such a diverse set of source organisms is atypical for ilarviruses, thus warranting further investigation. In this study, modern and classical virological tools were combined to accelerate the characterization of SnIV1. Through HTS-based virome surveys, mining of sequence read archive datasets, and a literature search, SnIV1 was further identified from diverse plant and non-plant sources globally. SnIV1 isolates showed relatively low variability compared with other phylogenetically related ilarviruses. Phylogenetic analyses showed a distinct basal clade of isolates from Europe, whereas the rest formed clades of mixed geographic origin. Furthermore, systemic infection of SnIV1 in Solanum villosum and its mechanical and graft transmissibility to solanaceous species were demonstrated. Near-identical SnIV1 genomes from the inoculum (S. villosum) and inoculated Nicotiana benthamiana were sequenced, thus partially fulfilling Koch's postulates. SnIV1 was shown to be seed-transmitted and potentially pollen-borne, has spherical virions, and possibly induces histopathological changes in infected N. benthamiana leaf tissues. Overall, this study provides information to better understand the diversity, global presence, and pathobiology of SnIV1; however, its possible emergence as a destructive pathogen remains uncertain. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Genetic Diversity, Evolution, and Diagnosis of Sugarcane Yellow Leaf Virus from 19 Sugarcane-Producing Locations Worldwide.

Sugarcane yellow leaf virus (SCYLV), the causal agent of yellow leaf, has been reported in an increasing number of sugarcane-growing locations since its first report in the 1990s in Brazil, Florida, and Hawaii. In this study, the genetic diversity of SCYLV was investigated using the genome coding sequence (5,561 to 5,612 nt) of 109 virus isolates from 19 geographical locations, including 65 new isolates from 16 geographical regions worldwide. These isolates were distributed in three major phylogenetic lineages (BRA, CUB, and REU), except for one isolate from Guatemala. Twenty-two recombination events were identified among the 109 isolates of SCYLV, thus confirming that recombination was a significant driving force in the genetic diversity and evolution of this virus. No temporal signal was found in the genomic sequence dataset, most likely because of the short temporal window of the 109 SCYLV isolates (1998 to 2020). Among 27 primers reported in the literature for the detection of the virus by RT-PCR, none matched 100% with all 109 SCYLV sequences, suggesting that the use of some primer pairs may not result in the detection of all virus isolates. Primers YLS111/YLS462, which were the first primer pair used by numerous research organizations to detect the virus by RT-PCR, failed to detect isolates belonging to the CUB lineage. In contrast, primer pair ScYLVf1/ScYLVr1 efficiently detected isolates of all three lineages. Continuous pursuit of knowledge of SCYLV genetic variability is therefore critical for effective diagnosis of yellow leaf, especially in virus-infected and mainly asymptomatic sugarcane plants.

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.

Ever-increasing viral diversity associated with the red imported fire ant Solenopsis invicta (Formicidae: Hymenoptera).

BACKGROUND: Advances in sequencing and analysis tools have facilitated discovery of many new viruses from invertebrates, including ants. Solenopsis invicta is an invasive ant that has quickly spread worldwide causing significant ecological and economic impacts. Its virome has begun to be characterized pertaining to potential use of viruses as natural enemies. Although the S. invicta virome is the best characterized among ants, most studies have been performed in its native range, with less information from invaded areas. METHODS: Using a metatranscriptome approach, we further identified and molecularly characterized virus sequences associated with S. invicta, in two introduced areas, U.S and Taiwan. The data set used here was obtained from different stages (larvae, pupa, and adults) of S. invicta life cycle. Publicly available RNA sequences from GenBank's Sequence Read Archive were downloaded and de novo assembled using CLC Genomics Workbench 20.0.1. Contigs were compared against the non-redundant protein sequences and those showing similarity to viral sequences were further analyzed. RESULTS: We characterized five putative new viruses associated with S. invicta transcriptomes. Sequence comparisons revealed extensive divergence across ORFs and genomic regions with most of them sharing less than 40% amino acid identity with those closest homologous sequences previously characterized. The first negative-sense single-stranded RNA virus genomic sequences included in the orders Bunyavirales and Mononegavirales are reported. In addition, two positive single-strand virus genome sequences and one single strand DNA virus genome sequence were also identified. While the presence of a putative tenuivirus associated with S. invicta was previously suggested to be a contamination, here we characterized and present strong evidence that Solenopsis invicta virus 14 (SINV-14) is a tenui-like virus that has a long-term association with the ant. Furthermore, based on virus sequence abundance compared to housekeeping genes, phylogenetic relationships, and completeness of viral coding sequences, our results suggest that four of five virus sequences reported, those being SINV-14, SINV-15, SINV-16 and SINV-17, may be associated to viruses actively replicating in the ant S. invicta. CONCLUSIONS: The present study expands our knowledge about viral diversity associated with S. invicta in introduced areas with potential to be used as biological control agents, which will require further biological characterization.

New Isolates of Pandoraviruses: Contribution to the Study of Replication Cycle Steps.

Giant viruses are complex members of the virosphere, exhibiting outstanding structural and genomic features. Among these viruses, the pandoraviruses are some of the most intriguing members, exhibiting giant particles and genomes presenting at up to 2.5 Mb, with many genes having no known function. In this work, we analyzed, by virological and microscopic methods, the replication cycle steps of three new pandoravirus isolates from samples collected in different regions of Brazil. Our data indicate that all analyzed pandoravirus isolates can deeply modify the Acanthamoeba cytoplasmic environment, recruiting mitochondria and membranes into and around the electron-lucent viral factories. We also observed that the viral factories start forming before the complete degradation of the cellular nucleus. Various patterns of pandoravirus particle morphogenesis were observed, and the assembly of the particles seemed to be started either by the apex or by the opposite side. On the basis of the counting of viral particles during the infection time course, we observed that pandoravirus particles could undergo exocytosis after their morphogenesis in a process that involved intense recruitment of membranes that wrapped the just-formed particles. The treatment of infected cells with brefeldin affected particle exocytosis in two of the three analyzed strains, indicating biological variability among isolates. Despite such particle exocytosis, the lysis of host cells also contributed to viral release. This work reinforces knowledge of and reveals important steps in the replication cycle of pandoraviruses.IMPORTANCE The emerging Pandoraviridae family is composed of some of the most complex viruses known to date. Only a few pandoravirus isolates have been described until now, and many aspects of their life cycle remain to be elucidated. A comprehensive description of the replication cycle is pivotal to a better understanding of the biology of the virus. For this report, we describe new pandoraviruses and used different methods to better characterize the steps of the replication cycle of this new group of viruses. Our results provide new information about the diversity and biology of these giant viruses.

Distantly Related Rotaviruses in Common Shrews, Germany, 2004-2014.

We screened samples from common shrews (Sorex araneus) collected in Germany during 2004-2014 and identified 3 genetically divergent rotaviruses. Virus protein 6 sequence similarities to prototype rotaviruses were low (64.5% rotavirus A, 50.1% rotavirus C [tentative species K], 48.2% rotavirus H [tentative species L]). Shrew-associated rotaviruses might have zoonotic potential.

Papaya Ringspot Virus Isolates From Papaya in Bangladesh: Detection, Characterization, and Distribution.

Papaya ringspot virus (PRSV) is the major constraint to papaya (Carica papaya) production in Bangladesh. Disease symptoms occurred in 90 to 100% of the plants surveyed. Full-length genomes of PRSV strains from severely infected papaya plants were determined using the Illumina NextSeq 500 platform, followed by Sanger DNA sequencing of viral genomes obtained by reverse-transcription PCR(RT-PCR). The genome sequences of two distinct PRSV strains, PRSV BD-1 (10,300 bp) and PRSV BD-2 (10,325 bp) were 74 and 83% identical to each other, respectively, at the nucleotide and amino acid levels. PRSV BD-1 and PRSV BD-2 were 74 to 75% and 79 to 88% identical, respectively, to other full-length PRSV sequences at the nucleotide level. Based on phylogenetic analysis, PRSV BD-2 was most closely related to PRSV-Meghalaya (MF356497) from papaya in India. PRSV BD-1 formed a branch distinct from the other PRSV sequences based on nucleotide and amino acid sequence comparisons. Comparisons of the genome sequences of these two strains with other sequenced PRSV genomes indicated two putative recombination events in PRSV BD-2. One recombinant event contained a 2,766-nucleotide fragment highly identical to PRSV-Meghalaya (MF356497). The other recombinant event contained a 5,105-nucleotide fragment highly identical to PRSV-China (KY933061). The occurrence rates of PRSV BD-1 and PRSV BD-2 in the sampled areas of Bangladesh were approximately 19 and 69%, respectively. Plants infected with both strains (11%) exhibited more severe symptoms than plants infected with either strain alone. The full-length genome sequences of these new PRSV strains and their distribution provide important information regarding the dynamics of papaya ringspot virus infections in papaya in Bangladesh.

Mechanisms and consequences of positive-strand RNA virus recombination.

Genetic recombination in positive-strand RNA viruses is a significant evolutionary mechanism that drives the creation of viral diversity by the formation of novel chimaeric genomes. The process and its consequences, for example the generation of viruses with novel phenotypes, has historically been studied by analysis of the end products. More recently, with an appreciation that there are both replicative and non-replicative mechanisms at work, and with new approaches and techniques to analyse intermediate products, the viral and cellular factors that influence the process are becoming understood. The major influence on replicative recombination is the fidelity of viral polymerase, although RNA structures and sequences may also have an impact. In replicative recombination the viral polymerase is necessary and sufficient, although roles for other viral or cellular proteins may exist. In contrast, non-replicative recombination appears to be mediated solely by cellular components. Despite these insights, the relative importance of replicative and non-replicative mechanisms is not clear. Using single-stranded positive-sense RNA viruses as exemplars, we review the current state of understanding of the processes and consequences of recombination.

Complete Genome Sequencing of Influenza A Viruses within Swine Farrow-to-Wean Farms Reveals the Emergence, Persistence, and Subsidence of Diverse Viral Genotypes.

Influenza A viruses (IAVs) are endemic in swine and represent a public health risk. However, there is limited information on the genetic diversity of swine IAVs within farrow-to-wean farms, which is where most pigs are born. In this longitudinal study, we sampled 5 farrow-to-wean farms for a year and collected 4,190 individual nasal swabs from three distinct pig subpopulations. Of these, 207 (4.9%) samples tested PCR positive for IAV, and 124 IAVs were isolated. We sequenced the complete genomes of 123 IAV isolates and found 31 H1N1, 26 H1N2, 63 H3N2, and 3 mixed IAVs. Based on the IAV hemagglutinin, seven different influenza A viral groups (VGs) were identified. Most of the remaining IAV gene segments allowed us to differentiate the same VGs, although an additional viral group was identified for gene segment 3 (PA). Moreover, the codetection of more than one IAV VG was documented at different levels (farm, subpopulation, and individual pigs), highlighting the environment for potential IAV reassortment. Additionally, 3 out of 5 farms contained IAV isolates (n = 5) with gene segments from more than one VG, and 79% of all the IAVs sequenced contained a signature mutation (S31N) in the matrix gene that has been associated with resistance to the antiviral amantadine. Within farms, some IAVs were detected only once, while others were detected for 283 days. Our results illustrate the maintenance and subsidence of different IAVs within swine farrow-to-wean farms over time, demonstrating that pig subpopulation dynamics are important to better understand the diversity and epidemiology of swine IAVs.IMPORTANCE On a global scale, swine are one of the main reservoir species for influenza A viruses (IAVs) and play a key role in the transmission of IAVs between species. Additionally, the 2009 IAV pandemics highlighted the role of pigs in the emergence of IAVs with pandemic potential. However, limited information is available regarding the diversity and distribution of swine IAVs on farrow-to-wean farms, where novel IAVs can emerge. In this study, we studied 5 swine farrow-to-wean farms for a year and characterized the genetic diversity of IAVs among three different pig subpopulations commonly housed on this type of farm. Using next-generation-sequencing technologies, we demonstrated the complex distribution and diversity of IAVs among the pig subpopulations studied. Our results demonstrated the dynamic evolution of IAVs within farrow-to-wean farms, which is crucial to improve health interventions to reduce the risk of transmission between pigs and from pigs to people.

Molecular epidemiology of Avian Rotaviruses Group A and D shed by different bird species in Nigeria.

BACKGROUND: Avian rotaviruses (RVs) cause gastrointestinal diseases of birds worldwide. However, prevalence, diversity, epidemiology and phylogeny of RVs remain largely under-investigated in Africa. METHODS: Fecal samples from 349 birds (158 symptomatic, 107 asymptomatic and 84 birds without recorded health status) were screened by reverse transcription PCR to detect RV groups A and D (RVA and RVD). Partial gene sequences of VP4, VP6, VP7 and NSP4 for RVA, and of VP6 and VP7 for RVD were obtained and analyzed to infer phylogenetic relationship. Fisher's exact test and logistic regression were applied to identify factors potentially influencing virus shedding in chickens. RESULTS: A high prevalence of RVA (36.1%; 126/349) and RVD (31.8%; 111/349) shedding was revealed in birds. In chickens, RV shedding was age-dependent and highest RVD shedding rates were found in commercial farms. No negative health effect could be shown, and RVA and RVD shedding was significantly more likely in asymptomatic chickens: RVA/RVD were detected in 51.9/48.1% of the asymptomatic chickens, compared to 18.9/29.7% of the symptomatic chickens (p < 0.001/p = 0.01). First RVA sequences were obtained from mallard ducks (Anas platyrhynchos) and guinea fowls (Numida meleagris). Phylogenetic analyses illustrated the high genetic diversity of RVA and RVD in Nigerian birds and suggested cross-species transmission of RVA, especially at live bird markets. Indeed, RVA strains highly similar to a recently published fox rotavirus (RVA/Fox-tc/ITA/288356/2011/G18P[17]) and distantly related to other avian RVs were detected in different bird species, including pigeons, ducks, guinea fowls, quails and chickens. CONCLUSION: This study provides new insights into epidemiology, diversity and classification of avian RVA and RVD in Nigeria. We show that cross-species transmission of host permissive RV strains occurs when different bird species are mixed.

Next generation sequencing technologies: tool to study avian virus diversity.

Increased globalisation, climatic changes and wildlife-livestock interface led to emergence of novel viral pathogens or zoonoses that have become serious concern to avian, animal and human health. High biodiversity and bird migration facilitate spread of the pathogen and provide reservoirs for emerging infectious diseases. Current classical diagnostic methods designed to be virus-specific or aim to be limited to group of viral agents, hinder identifying of novel viruses or viral variants. Recently developed approaches of next-generation sequencing (NGS) provide culture-independent methods that are useful for understanding viral diversity and discovery of novel virus, thereby enabling a better diagnosis and disease control. This review discusses the different possible steps of a NGS study utilizing sequence-independent amplification, high-throughput sequencing and bioinformatics approaches to identify novel avian viruses and their diversity. NGS lead to the identification of a wide range of new viruses such as picobirnavirus, picornavirus, orthoreovirus and avian gamma coronavirus associated with fulminating disease in guinea fowl and is also used in describing viral diversity among avian species. The review also briefly discusses areas of viral-host interaction and disease associated causalities with newly identified avian viruses.

Prediction of viral spillover risk based on the mass action principle.

Infectious zoonotic disease emergence, through spillover events, is of global concern and has the potential to cause significant harm to society, as recently demonstrated by COVID-19. More than 70% of the 400 infectious diseases that emerged in the past five decades have a zoonotic origin, including all recent pandemics. There have been several approaches used to predict the risk of spillover through some of the known or suspected infectious disease emergence drivers, largely using correlative approaches. Here, we predict the spatial distribution of spillover risk by approximating general transmission through animal and human interactions. These mass action interactions are approximated through the multiplication of the spatial distribution of zoonotic virus diversity and human population density. Although our results indicate higher risk in regions along the equator and in Southeast Asia where both virus diversity and human population density are high, it should be noted that this is primarily a conceptual exercise. We compared our spillover risk map to key factors, including the model inputs of zoonotic virus diversity estimate map, human population density map, and the spatial distribution of species richness. Despite the limitations of this approach, this viral spillover map is a step towards developing a more comprehensive spillover risk prediction system to inform global monitoring.

Integration of Rubus yellow net virus in the raspberry genome: A story centuries in the making.

Rubus yellow net virus (RYNV) belongs to genus Badnavirus. Badnaviruses are found in plants as endogenous, inactive or activatable sequences, and/or in episomal (infectious and active) forms. To assess the state of RYNV in Rubus germplasm, we sequenced the genomes of various cultivars and mined eight raspberry whole genome datasets. Bioinformatics analysis revealed the presence of a diverse array of endogenous RYNV (endoRYNV) sequences that differ significantly in their structure; some lineages have nearly complete, yet non-functional genomes whereas others have rudimentary, short sequence fragments. We developed assays to genotype the main lineages as well as the only known episomal lineage present in the United States. This study discloses the widespread presence of endoRYNVs in commercial raspberries, likely because breeding efforts have focused on a limited pool of germplasm that harbored endoRYNVs.

Diversity and connectedness of brine shrimp viruses in global hypersaline ecosystems.

Brine shrimp (Artemia) has existed on Earth for 400 million years and has major ecological importance in hypersaline ecosystems. As a crucial live food in aquaculture, brine shrimp cysts have become one of the most important aquatic products traded worldwide. However, our understanding of the biodiversity, prevalence and global connectedness of viruses in brine shrimp is still very limited. A total of 143 batches of brine shrimp (belonging to seven species) cysts were collected from six continents including 21 countries and more than 100 geographic locations worldwide during 1977-2019. In total, 55 novel RNA viruses were identified, which could be assigned to 18 different viral families and related clades. Eleven viruses were dsRNA viruses, 16 were +ssRNA viruses, and 28 were-ssRNA viruses. Phylogenetic analyses of the RNA-directed RNA polymerase (RdRp) showed that brine shrimp viruses were often grouped with viruses isolated from other invertebrates and fungi. Remarkably, most brine shrimp viruses were related to those from different hosts that might feed on brine shrimp or share the same ecological niche. A notable case was the novel brine shrimp noda-like virus 3, which shared 79.25% (RdRp) and 63.88% (capsid proteins) amino acid identity with covert mortality nodavirus (CMNV) that may cause losses in aquaculture. In addition, both virome composition and phylogenetic analyses revealed global connectedness in certain brine shrimp viruses, particularly among Asia and Northern America. This highlights the incredible species diversity of viruses in these ancient species and provides essential data for the prevalence of RNA viruses in the global aquaculture industry. More broadly, these findings provide novel insights into the previously unrecognized RNA virosphere in hypersaline ecosystems worldwide and demonstrate that human activity might have driven the global connectedness of brine shrimp viruses.

Synergism between remdesivir and ribavirin leads to SARS-CoV-2 extinction in cell culture.

BACKGROUND AND PURPOSE: There is a need for effective anti-COVID-19 treatments, mainly for individuals at risk of severe disease such as the elderly and the immunosuppressed. Drug repositioning has proved effective in identifying drugs that can find a new application for the control of coronavirus disease, in particular COVID-19. The purpose of the present study was to find synergistic antiviral combinations for COVID-19 based on lethal mutagenesis. EXPERIMENTAL APPROACH: The effect of combinations of remdesivir and ribavirin on the infectivity of SARS-CoV-2 in cell culture has been tested. Viral populations were monitored by ultra-deep sequencing, and the decrease of infectivity as a result of the treatment was measured. KEY RESULTS: Remdesivir and ribavirin exerted a synergistic inhibitory activity against SARS-CoV-2, quantified both by CompuSyn (Chou-Talalay method) and Synergy Finder (ZIP-score model). In serial passage experiments, virus extinction was readily achieved with remdesivir-ribavirin combinations at concentrations well below their cytotoxic 50 value, but not with the drugs used individually. Deep sequencing of treated viral populations showed that remdesivir, ribavirin, and their combinations evoked significant increases of the number of viral mutations and haplotypes, as well as modification of diversity indices that characterize viral quasi-species. CONCLUSION AND IMPLICATIONS: SARS-CoV-2 extinction can be achieved by synergistic combination treatments based on lethal mutagenesis. In addition, the results offer prospects of triple drug treatments for effective SARS-CoV-2 suppression.

Exploring Virus Diversity in the Potato leafhopper (Empoasca fabae), an Economically Important Agricultural Pest.

The potato leafhopper (Empoasca fabae, PLH) is a serious pest that feeds on a wide range of agricultural crops and is found throughout the United States but is not known to be a vector for plant-infecting viruses. We probed the diversity of virus sequences in field populations of PLH collected from four Midwestern states: Illinois, Indiana, Iowa, and Minnesota. High-throughput sequencing data from total RNAs extracted from PLH were used to assemble sequences of fifteen positive-stranded RNA viruses, two negative-stranded RNA viruses, and one DNA virus. These sequences included ten previously described plant viruses and eight putative insect-infecting viruses. All but one of the insect-specific viruses were novel and included three solemoviruses, one iflavirus, one phenuivirus, one lispivirus, and one ambidensovirus. Detailed analyses of the novel genome sequences and their evolutionary relationships with related family members were conducted. Our study revealed a diverse group of plant viruses circulating in the PLH population and discovered novel insect viruses, expanding knowledge on the untapped virus diversity in economically important crop pests. Our findings also highlight the importance of monitoring the emergence and circulation of plant-infecting viruses in agriculturally important arthropod pests.

Dothistroma septosporum and Dothistroma pini, the causal agents of Dothistroma needle blight, are infected by multiple viruses.

Dothistroma septosporum and Dothistroma pini are severe foliar pathogens of conifers. They infect a broad spectrum of hosts (mainly Pinus spp.), causing chlorosis, defoliation of needles, and eventually the death of pine trees in extreme cases. Mycoviruses represent a novel and innovative avenue for controlling pathogens. To search for possible viruses hosted by Dothistroma spp. we screened a subset of isolates (20 strains of D. septosporum and one D. pini) originating from the Czech Republic, Slovenia, Italy, Austria and Ireland for viral dsRNA segments. Only five of them showed the presence of dsRNA segments. A total of 21 fungal isolates were prepared for total RNA extractions. RNA samples were pooled, and two separate RNA libraries were constructed for stranded total RNA sequencing. RNA-Seq data processing, pairwise sequence comparisons (PASC) and phylogenetic analyses revealed the presence of thirteen novel putative viruses with varying genome types: seven negative-sense single-stranded RNA viruses, including six bunya-like viruses and one new member of the order Mononegavirales; three positive-sense single-stranded RNA viruses, two of which are similar to those of the family Narnaviridae, while the genome of the third correspond to those of the family Gammaflexiviridae; and three double-stranded RNA viruses, comprising two novel members of the family Chrysoviridae and a potentially new species of gammapartitivirus. The results were confirmed with RT-PCR screening that the fungal pathogens hosted all the viruses and showed that particular fungal strains harbour multiple virus infections and that they are transmitted vertically. In this study, we described the narnavirus infecting D. pini. To our knowledge, this is the first virus discovered in D. pini.

Virome analysis of an ectomycorrhizal fungus Suillus luteus revealing potential evolutionary implications.

Suillus luteus is a widespread edible ectomycorrhizal fungus that holds significant importance in both ecological and economic value. Mycoviruses are ubiquitous infectious agents hosted in different fungi, with some known to exert beneficial or detrimental effects on their hosts. However, mycoviruses hosted in ectomycorrhizal fungi remain poorly studied. To address this gap in knowledge, we employed next-generation sequencing (NGS) to investigate the virome of S. luteus. Using BLASTp analysis and phylogenetic tree construction, we identified 33 mycovirus species, with over half of them belonging to the phylum Lenarviricota, and 29 of these viruses were novel. These mycoviruses were further grouped into 11 lineages, with the discovery of a new negative-sense single-stranded RNA viral family in the order Bunyavirales. In addition, our findings suggest the occurrence of cross-species transmission (CST) between the fungus and ticks, shedding light on potential evolutionary events that have shaped the viral community in different hosts. This study is not only the first study to characterize mycoviruses in S. luteus but highlights the enormous diversity of mycoviruses and their implications for virus evolution.