Molecular characterization of the complete genome of a novel ormycovirus infecting the ectomycorrhizal fungus Hortiboletus rubellus.

Advancements in high-throughput sequencing and the development of new bioinformatics tools for large-scale data analysis play a crucial role in uncovering virus diversity and enhancing our understanding of virus evolution. The discovery of the ormycovirus clades, a group of RNA viruses that are phylogenetically distinct from all known Riboviria members and are found in fungi, highlights the value of these tools for the discovery of novel viruses. The aim of this study was to examine viral populations in fungal hosts to gain insights into the diversity, evolution, and classification of these viruses. Here, we report the molecular characterization of a newly discovered ormycovirus, which we have named "Hortiboletus rubellus ormycovirus 1" (HrOMV1), that was found in the ectomycorrhizal fungus Hortiboletus rubellus. The bipartite genome of HrOMV1, whose nucleotide sequence was determined by HTS and RLM-RACE, consists of two RNA segments (RNA1 and RNA2) that exhibit similarity to those of previously studied ormycoviruses in their organization and the proteins they encode. The presence of upstream, in-frame AUG triplets in the 5' termini of both RNA segments suggests that HrOMV1, like certain other ormycoviruses, employs a non-canonical translation initiation strategy. Phylogenetic analysis showed that HrOMV1 is positioned within the gammaormycovirus clade. Its putative RNA-dependent RNA polymerase (RdRp) exhibits sequence similarity to those of other gammaormycovirus members, the most similarity to that of Termitomyces ormycovirus 1, with 33.05% sequence identity. This protein was found to contain conserved motifs that are crucial for RNA replication, including the distinctive GDQ catalytic triad observed in gammaormycovirus RdRps. The results of this study underscore the significance of investigating the ecological role of mycoviruses in mycorrhizal fungi. This is the first report of an ormycovirus infecting a member of the ectomycorrhizal genus Hortiboletus.

Characterization of Mycoviruses in Armillaria ostoyae and A. cepistipes in the Czech Republic.

Members of the genus Armillaria are widespread forest pathogens against which effective protection has not yet been developed. Due to their longevity and the creation of large-scale cloning of Armillaria individuals, the use of mycoviruses as biocontrol agents (BCAs) against these pathogens could be an effective alternative. This work describes the detection and characterization of viruses in Armillaria spp. collected in the Czech Republic through the application of stranded total RNA sequencing. A total of five single-stranded RNA viruses were detected in Armillaria ostoyae and A. cepistipes, including viruses of the family Tymoviridae and four viruses belonging to the recently described "ambivirus" group with a circular ambisense genome arrangement. Both hammerhead (HHRz) and hairpin (HpRz) ribozymes were detected in all the ambiviricot sequences. Armillaria viruses were compared through phylogenetic analysis and confirmed their specific host by direct RT-PCR. One virus appears to infect both Armillaria species, suggesting the occurrence of interspecies transmission in nature.

U-CAN-seq: A Universal Competition Assay by Nanopore Sequencing.

RNA viruses quickly evolve subtle genotypic changes that can have major impacts on viral fitness and host range, with potential consequences for human health. It is therefore important to understand the evolutionary fitness of novel viral variants relative to well-studied genotypes of epidemic viruses. Competition assays are an effective and rigorous system with which to assess the relative fitness of viral genotypes. However, it is challenging to quickly and cheaply distinguish and quantify fitness differences between very similar viral genotypes. Here, we describe a protocol for using reverse transcription PCR in combination with commercial nanopore sequencing services to perform competition assays on untagged RNA viruses. Our assay, called the Universal Competition Assay by Nanopore Sequencing (U-CAN-seq), is relatively cheap and highly sensitive. We used a well-studied N24A mutation in the chikungunya virus (CHIKV) nsp3 gene to confirm that we could detect a competitive disadvantage using U-CAN-seq. We also used this approach to show that mutations to the CHIKV 5' conserved sequence element that disrupt sequence but not structure did not affect the fitness of CHIKV. However, similar mutations to an adjacent CHIKV stem loop (SL3) did cause a fitness disadvantage compared to wild-type CHIKV, suggesting that structure-independent, primary sequence determinants in this loop play an important role in CHIKV biology. Our novel findings illustrate the utility of the U-CAN-seq competition assay.

SPM4GAC: SPM based approach for genome analysis and classification of macromolecules.

Genome sequence analysis and classification play critical roles in properly understanding an organism's main characteristics, functionalities, and changing (evolving) nature. However, the rapid expansion of genomic data makes genome sequence analysis and classification a challenging task due to the high computational requirements, proper management, and understanding of genomic data. Recently proposed models yielded promising results for the task of genome sequence classification. Nevertheless, these models often ignore the sequential nature of nucleotides, which is crucial for revealing their underlying structure and function. To address this limitation, we present SPM4GAC, a sequential pattern mining (SPM)-based framework to analyze and classify the macromolecule genome sequences of viruses. First, a large dataset containing the genome sequences of various RNA viruses is developed and transformed into a suitable format. On the transformed dataset, algorithms for SPM are used to identify frequent sequential patterns of nucleotide bases. The obtained frequent sequential patterns of bases are then used as features to classify different viruses. Ten classifiers are employed, and their performance is assessed by using several evaluation measures. Finally, a performance comparison of SPM4GAC with state-of-the-art methods for genome sequence classification/detection reveals that SPM4GAC performs better than those methods.

Viromes of Haemaphysalis longicornis reveal different viral abundance and diversity in free and engorged ticks.

Haemaphysalis longicornis ticks, commonly found in East Asia, can transmit various pathogenic viruses, including the severe fever with thrombocytopenia syndrome virus (SFTSV) that has caused febrile diseases among humans in Hubei Province. However, understanding of the viromes of H. longicornis was limited, and the prevalence of viruses among H. longicornis ticks in Hubei was not well clarified. This study investigates the viromes of both engorged (fed) and free (unfed) H. longicornis ticks across three mountainous regions in Hubei Province from 2019 to 2020. RNA-sequencing analysis identified viral sequences that were related to 39 reference viruses belonging to unclassified viruses and seven RNA viral families, namely Chuviridae, Nairoviridae, Orthomyxoviridae, Parvoviridae, Phenuiviridae, Rhabdoviridae, and Totiviridae. Viral abundance and diversity in these ticks were analysed, and phylogenetic characteristics of the Henan tick virus (HNTV), Dabieshan tick virus (DBSTV), Okutama tick virus (OKTV), and Jingmen tick virus (JMTV) were elucidated based on their full genomic sequences. Prevalence analysis demonstrated that DBSTV was the most common virus found in individual H. longicornis ticks (12.59%), followed by HNTV (0.35%), whereas JMTV and OKTV were not detected. These results improve our understanding of H. longicornis tick viromes in central China and highlight the role of tick feeding status and geography in shaping the viral community. The findings of new viral strains and their potential impact on public health raise the need to strengthen surveillance efforts for comprehensively assessing their spillover potentials.

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.

Mycovirus Hunting Revealed the Presence of Diverse Viruses in a Single Isolate of the Phytopathogenic Fungus Diplodia seriata From Pakistan.

Diplodia seriata in the family Botryosphaeriaceae is a cosmopolitan phytopathogenic fungus and is responsible for causing cankers, fruit rot and leaf spots on economically important plants. In this study, we characterized the virome of a single Pakistani strain (L3) of D. seriata. Several viral-like contig sequences were obtained via a previously conducted next-generation sequencing analysis. Multiple infection of the L3 strain by eight RNA mycoviruses was confirmed through RT-PCR using total RNA samples extracted from this strain; the entire genomes were determined via Sanger sequencing of RT-PCR and RACE clones. A BLAST search and phylogenetic analyses indicated that these eight mycoviruses belong to seven different viral families. Four identified mycoviruses belong to double-stranded RNA viral families, including Polymycoviridae, Chrysoviridae, Totiviridae and Partitiviridae, and the remaining four identified mycoviruses belong to single-stranded RNA viral families, i.e., Botourmiaviridae, and two previously proposed families "Ambiguiviridae" and "Splipalmiviridae". Of the eight, five mycoviruses appear to represent new virus species. A morphological comparison of L3 and partially cured strain L3ht1 suggested that one or more of the three viruses belonging to Polymycoviridae, "Splipalmiviridae" and "Ambiguiviridae" are involved in the irregular colony phenotype of L3. To our knowledge, this is the first report of diverse virome characterization from D. seriata.

Efficiency, sensitivity and specificity of a quantitative real-time PCR assay for Tilapia Lake virus (TiLV).

Tilapia lake virus (TiLV) is an emerging viral pathogen of tilapiines worldwide in wild and farmed tilapia. TiLV is an orthomyxo-like, negative sense segmented RNA virus, belonging to genus Tilapinevirus, family Amnoonviridae. Here we developed a quantitative real-time PCR (qRT-PCR) assay testing primer sets targeting the 10 segments of TiLV. Sensitivity, specificity, efficiency and reproducibility of these assays were examined. Detection sensitivity was equivalent to 2 TCID50/ml when tested on supernatants from cell culture-grown TiLV. Specificity tests showed that all primer sets amplified their respective TiLV segments, and standard curves showed linear correlation of R2 > 0.998 and amplification efficiencies between 93 % and 98 %. Intra- and inter-assay coefficients of variation (CV %) were in the range of 0.0 %- 2.6 % and 0.0 %- 5.9 %, respectively. Sensitivity tests showed that primer sets targeting segments 1, 2, 3 and 4 had the highest detection sensitivities (100.301 TCID50/ml). The qRT-PCR used for detection of viral genome in TiLV infected organs gave virus titers equivalent to 3.80 log10, 3.94 log10 and 3.52 log10 TCID50/ml for brain, kidney and liver tissues, respectively as calculated on the basis of Ct values. These findings suggest that primer optimization for qPCR should not only focus on attaining high amplification efficiency but also sensitivity comparison of primer sets targeting different viral segments in order to develop a method with the highest sensitivity.

Diversity and ecological footprint of Global Ocean RNA viruses.

DNA viruses are increasingly recognized as influencing marine microbes and microbe-mediated biogeochemical cycling. However, little is known about global marine RNA virus diversity, ecology, and ecosystem roles. In this study, we uncover patterns and predictors of marine RNA virus community- and "species"-level diversity and contextualize their ecological impacts from pole to pole. Our analyses revealed four ecological zones, latitudinal and depth diversity patterns, and environmental correlates for RNA viruses. Our findings only partially parallel those of cosampled plankton and show unexpectedly high polar ecological interactions. The influence of RNA viruses on ecosystems appears to be large, as predicted hosts are ecologically important. Moreover, the occurrence of auxiliary metabolic genes indicates that RNA viruses cause reprogramming of diverse host metabolisms, including photosynthesis and carbon cycling, and that RNA virus abundances predict ocean carbon export.

A Capsidless Virus Is trans-Encapsidated by a Bisegmented Botybirnavirus.

RNA viruses usually have linear genomes and are encapsidated by their own capsids. Here, we newly identified four mycoviruses and two previously reported mycoviruses (a fungal reovirus and a botybirnavirus) in the hypovirulent strain SCH941 of Sclerotinia sclerotiorum. One of the newly discovered mycoviruses, Sclerotinia sclerotiorum yadokarivirus 1 (SsYkV1), with a nonsegmented positive-sense single-stranded RNA (+ssRNA) genome, was molecularly characterized. SsYkV1 is 5,256 nucleotides (nt) in length, excluding the poly(A) structure, and has a large open reading frame that putatively encodes a polyprotein with the RNA-dependent RNA polymerase (RdRp) domain and a 2A-like motif. SsYkV1 was phylogenetically positioned into the family Yadokariviridae and was most closely related to Rosellinia necatrix yadokarivirus 2 (RnYkV2), with 40.55% identity (78% coverage). Although SsYkV1 does not encode its own capsid protein, the RNA and RdRp of SsYkV1 are trans-encapsidated in virions of Sclerotinia sclerotiorum botybirnavirus 3 (SsBV3), a bisegmented double-stranded RNA (dsRNA) mycovirus within the genus Botybirnavirus. In this way, SsYkV1 likely replicates inside the heterocapsid comprised of the SsBV3 capsid protein, like a dsRNA virus. SsYkV1 has a limited impact on the biological features of S. sclerotiorum. This study represents an example of a yadokarivirus trans-encapsidated by an unrelated dsRNA virus, which greatly deepens our knowledge and understanding of the unique life cycles of RNA viruses. IMPORTANCE RNA viruses typically encase their linear genomes in their own capsids. However, a capsidless +ssRNA virus (RnYkV1) highjacks the capsid of a nonsegmented dsRNA virus for the trans-encapsidation of its own RNA and RdRp. RnYkV1 belongs to the family Yadokariviridae, which already contains more than a dozen mycoviruses. However, it is unknown whether other yadokariviruses except RnYkV1 are also hosted by a heterocapsid, although dsRNA viruses with capsid proteins were detected in fungi harboring yadokarivirus. It is noteworthy that almost all presumed partner dsRNA viruses of yadokariviruses belong to the order Ghabrivirales (most probably a totivirus or toti-like virus). Here, we found a capsidless +ssRNA mycovirus, SsYkV1, from hypovirulent strain SCH941 of S. sclerotiorum, and the RNA and RdRp of this mycovirus are trans-encapsidated in virions of a bisegmented dsRNA virus within the free-floating genus Botybirnavirus. Our results greatly expand our knowledge of the unique life cycles of RNA viruses.

A Novel Heptasegmented Positive-Sense Single-Stranded RNA Virus from the Phytopathogenic Fungus Colletotrichum fructicola.

In this study, a novel positive-sense single-stranded RNA (+ssRNA) mycovirus, tentatively named Colletotrichum fructicola RNA virus 1 (CfRV1), was identified in the phytopathogenic fungus Colletotrichum fructicola. CfRV1 has seven genomic components, encoding seven proteins from open reading frames (ORFs) flanked by highly conserved untranslated regions (UTRs). Proteins encoded by ORFs 1, 2, 3, 5, and 6 are more similar to the putative RNA-dependent RNA polymerase (RdRp), hypothetical protein (P2), methyltransferase, and two hypothetical proteins of Hadaka virus 1 (HadV1), a capsidless 10- or 11-segmented +ssRNA virus, while proteins encoded by ORFs 4 and 7 showed no detectable similarity to any known proteins. Notably, proteins encoded by ORFs 1 to 3 also share considerably high similarity with the corresponding proteins of polymycoviruses. Phylogenetic analysis conducted based on the amino acid sequence of CfRV1 RdRp and related viruses placed CfRV1 and HadV1 together in the same clade, close to polymycoviruses and astroviruses. CfRV1-infected C. fructicola strains demonstrate a moderately attenuated growth rate and virulence compared to uninfected isolates. CfRV1 is capsidless and potentially encapsulated in vesicles inside fungal cells, as revealed by transmission electron microscopy. CfRV1 and HadV1 are +ssRNA mycoviruses closely related to polymycoviruses and astroviruses, represent a new linkage between +ssRNA viruses and the intermediate double-stranded RNA (dsRNA) polymycoviruses, and expand our understanding of virus diversity, taxonomy, evolution, and biological traits. IMPORTANCE A scenario proposing that dsRNA viruses evolved from +ssRNA viruses is still considered controversial due to intergroup knowledge gaps in virus diversity. Recently, polymycoviruses and hadakaviruses were found as intermediate dsRNA and +ssRNA stages, respectively, between +ssRNA and dsRNA viruses. Here, we identified a novel +ssRNA mycovirus, Colletotrichum fructicola RNA virus 1 (CfRV1), isolated from Colletotrichum fructicola in China. CfRV1 is phylogenetically related to the 10- or 11-segmented Hadaka virus 1 (HadV1) but consists of only seven genomic segments encoding two novel proteins. CfRV1 is naked and may be encapsulated in vesicles inside fungal cells, representing a potential novel lifestyle for multisegmented RNA viruses. CfRV1 and HadV1 are intermediate +ssRNA mycoviruses in the linkage between +ssRNA viruses and the intermediate dsRNA polymycoviruses and expand our understanding of virus diversity, taxonomy, and evolution.

Inferring Protein Function in an Emerging Virus: Detection of the Nucleoprotein in Tilapia Lake Virus.

Emerging viruses impose global threats to animal and human populations and may bear novel genes with limited homology to known sequences, necessitating the development of novel approaches to infer and test protein functions. This challenge is dramatically evident in tilapia lake virus (TiLV), an emerging "orthomyxo-like" virus that threatens the global tilapia aquaculture and food security of millions of people. The majority of TiLV proteins have no homology to known sequences, impeding functionality assessments. Using a novel bioinformatics approach, we predicted that TiLV's Protein 4 encodes the nucleoprotein, a factor essential for viral RNA replication. Multiple methodologies revealed the expected properties of orthomyxoviral nucleoproteins. A modified yeast three-hybrid assay detected Protein 4-RNA interactions, which were independent of the RNA sequence, and identified specific positively charged residues involved. Protein 4-RNA interactions were uncovered by R-DeeP and XRNAX methodologies. Immunoelectron microscopy found that multiple Protein 4 copies localized along enriched ribonucleoproteins. TiLV RNA from cells and virions coimmunoprecipitated with Protein 4. Immunofluorescence microscopy detected Protein 4 in the cytoplasm and nuclei, and nuclear Protein 4 increased upon CRM1 inhibition, suggesting CRM1-dependent nuclear export of TiLV RNA. Together, these data reveal TiLV's nucleoprotein and highlight the ability to infer protein functionality, including novel RNA-binding proteins, in emerging pathogens. These are important in light of the expected discovery of many unknown viruses and the zoonotic potential of such pathogens. IMPORTANCE Tilapia is an important source of dietary protein, especially in developing countries. Massive losses of tilapia were identified worldwide, risking the food security of millions of people. Tilapia lake virus (TiLV) is an emerging pathogen responsible for these disease outbreaks. TiLV's genome encodes 10 major proteins, 9 of which show no homology to other known viral or cellular proteins, hindering functionality assessment of these proteins. Here, we describe a novel bioinformatics approach to infer the functionality of TiLV proteins, which predicted Protein 4 as the nucleoprotein, a factor essential for viral RNA replication. We provided experimental support for this prediction by applying multiple molecular, biochemical, and imaging approaches. Overall, we illustrate a strategy for functional analyses in viral discovery. The strategy is important in light of the expected discovery of many unknown viruses and the zoonotic potential of such pathogens.

One-Enzyme RTX-PCR for the Detection of RNA Viruses from Multiple Virus Genera and Crop Plants.

Reverse transcription PCR (RT-PCR) is a popular method for detecting RNA viruses in plants. RT-PCR is usually performed in a classical two-step procedure: in the first step, cDNA is synthesized by reverse transcriptase (RT), followed by PCR amplification by a thermostable polymerase in a separate tube in the second step. However, one-step kits containing multiple enzymes optimized for RT and PCR amplification in a single tube can also be used. Here, we describe an RT-PCR single-enzyme assay based on an RTX DNA polymerase that has both RT and polymerase activities. The expression plasmid pET_RTX_(exo-) was transferred to various E. coli genotypes that either compensated for codon bias (Rosetta-gami 2) or contained additional chaperones to promote solubility (BL21 (DE3) with plasmids pKJE8 or pTf2). The RTX enzyme was then purified and used for the RT-PCR assay. Several purified plant viruses (TMV, PVX, and PVY) were used to determine the efficiency of the assay compared to a commercial one-step RT-PCR kit. The RT-PCR assay with the RTX enzyme was validated for the detection of viruses from different genera using both total RNA and crude sap from infected plants. The detection endpoint of RTX-PCR for purified TMV was estimated to be approximately 0.01 pg of the whole virus per 25 µL reaction, corresponding to 6 virus particles/µL. Interestingly, the endpoint for detection of TMV from crude sap was also 0.01 pg per reaction in simulated crude plant extracts. The longest RNA fragment that could be amplified in a one-tube arrangement was 2379 bp long. The longest DNA fragment that could be amplified during a 10s extension was 6899 bp long. In total, we were able to detect 13 viruses from 11 genera using RTX-PCR. For each virus, two to three specific fragments were amplified. The RT-PCR assay using the RTX enzyme described here is a very robust, inexpensive, rapid, easy to perform, and sensitive single-enzyme assay for the detection of plant viruses.

Detailed Analyses of Molecular Interactions between Favipiravir and RNA Viruses In Silico.

There are currently no antiviral agents for human metapneumovirus (HMPV), respiratory syncytial virus (RSV), mumps virus (MuV), or measles virus (MeV). Favipiravir has been developed as an anti-influenza agent, and this agent may be effective against these viruses in vitro. However, the molecular mechanisms through which the agent affects virus replication remain to be fully elucidated. Thus, to clarify the detailed molecular interactions between favipiravir and the RNA-dependent RNA polymerase (RdRp) of HMPV, RSV, MuV, MeV, and influenza virus, we performed in silico studies using authentic bioinformatics technologies. As a result, we found that the active form of favipiravir (favipiravir ribofuranosyl-5'-triphosphate [F-RTP]) can bind to the RdRp active sites of HMPV, RSV, MuV, and MeV. The aspartic acid residue of RdRp active sites was involved in the interaction. Moreover, F-RTP was incorporated into the growing viral RNA chain in the presence of nucleotide triphosphate and magnesium ions. The results suggested that favipiravir shows two distinct mechanisms in various viruses: RdRp active site inhibition and/or genome replication inhibition.

Metagenome-Assembled Viral Genomes Analysis Reveals Diversity and Infectivity of the RNA Virome of Gerbillinae Species.

Rodents are a known reservoir for extensive zoonotic viruses, and also possess a propensity to roost in human habitation. Therefore, it is necessary to identify and catalogue the potentially emerging zoonotic viruses that are carried by rodents. Here, viral metagenomic sequencing was used for zoonotic virus detection and virome characterization on 32 Great gerbils of Rhombomys opimus, Meriones meridianus, and Meiiones Unguiculataus species in Xinjiang, Northwest China. In total, 1848 viral genomes that are potentially pathogenic to rodents and humans, as well as to other wildlife, were identified namely Retro-, Flavi-, Pneumo-, Picobirna-, Nairo-, Arena-, Hepe-, Phenui-, Rhabdo-, Calici-, Reo-, Corona-, Orthomyxo-, Peribunya-, and Picornaviridae families. In addition, a new genotype of rodent Hepacivirus was identified in heart and lung homogenates of seven viscera pools and phylogenetic analysis revealed the closest relationship to rodent Hepacivirus isolate RtMm-HCV/IM2014 that was previously reported to infect rodents from Inner Mongolia, China. Moreover, nine new genotype viral sequences that corresponded to Picobirnaviruses (PBVs), which have a bi-segmented genome and belong to the family Picobirnaviridae, comprising of three segment I and six segment II sequences, were identified in intestines and liver of seven viscera pools. In the two phylogenetic trees that were constructed using ORF1 and ORF2 of segment I, the three segment I sequences were clustered into distinct clades. Additionally, phylogenetic analysis showed that PBV sequences were distributed in the whole tree that was constructed using the RNA-dependent RNA polymerase (RdRp) gene of segment II with high diversity, sharing 68.42-82.67% nucleotide identities with other genogroup I and genogroup II PBV strains based on the partial RdRp gene. By RNA sequencing, we found a high degree of biodiversity of Retro-, Flavi-, Pneumo-, and Picobirnaridae families and other zoonotic viruses in gerbils, indicating that zoonotic viruses are a common presence in gerbils from Xinjiang, China. Therefore, further research is needed to determine the zoonotic potential of these viruses that are carried by other rodent species from different ecosystems and wildlife in general.

Viruses Infecting Greenhood Orchids (Pterostylidinae) in Eastern Australia.

The Australasian biogeographic realm is a major centre of diversity for orchids, with every subfamily of the Orchidaceae represented and high levels of endemism at the species rank. It is hypothesised that there is a commensurate diversity of viruses infecting this group of plants. In this study, we have utilised high-throughput sequencing to survey for viruses infecting greenhood orchids (Pterostylidinae) in New South Wales and the Australian Capital Territory. The main aim of this study was to characterise Pterostylis blotch virus (PtBV), a previously reported but uncharacterised virus that had been tentatively classified in the genus Orthotospovirus. This classification was confirmed by genome sequencing, and phylogenetic analyses suggested that PtBV is representative of a new species that is possibly indigenous to Australia as it does not belong to either the American or Eurasian clades of orthotospoviruses. Apart from PtBV, putative new viruses in the genera Alphaendornavirus, Amalgavirus, Polerovirus and Totivirus were discovered, and complete genome sequences were obtained for each virus. It is concluded that the polerovirus is likely an example of an introduced virus infecting a native plant species in its natural habitat, as this virus is probably vectored by an aphid, and Australia has a depauperate native aphid fauna that does not include any species that are host-adapted to orchids.

Four Novel Mycoviruses from the Hypovirulent Botrytis cinerea SZ-2-3y Isolate from Paris polyphylla: Molecular Characterisation and Mitoviral Sequence Transboundary Entry into Plants.

A hypovirulent SZ-2-3y strain isolated from diseased Paris polyphylla was identified as Botrytis cinerea. Interestingly, SZ-2-3y was coinfected with a mitovirus, two botouliviruses, and a 3074 nt fusarivirus, designated Botrytis cinerea fusarivirus 8 (BcFV8); it shares an 87.2% sequence identity with the previously identified Botrytis cinerea fusarivirus 6 (BcFV6). The full-length 2945 nt genome sequence of the mitovirus, termed Botrytis cinerea mitovirus 10 (BcMV10), shares a 54% sequence identity with Fusarium boothii mitovirus 1 (FbMV1), and clusters with fungus mitoviruses, plant mitoviruses and plant mitochondria; hence BcMV10 is a new Mitoviridae member. The full-length 2759 nt and 2812 nt genome sequences of the other two botouliviruses, named Botrytis cinerea botoulivirus 18 and 19 (BcBoV18 and 19), share a 40% amino acid sequence identity with RNA-dependent RNA polymerase protein (RdRp), and these are new members of the Botoulivirus genus of Botourmiaviridae. Horizontal transmission analysis showed that BcBoV18, BcBoV19 and BcFV8 are not related to hypovirulence, suggesting that BcMV10 may induce hypovirulence. Intriguingly, a partial BcMV10 sequence was detected in cucumber plants inoculated with SZ-2-3y mycelium or pXT1/BcMV10 agrobacterium. In conclusion, we identified a hypovirulent SZ-2-3y fungal strain from P. polyphylla, coinfected with four novel mycoviruses that could serve as potential biocontrol agents. Our findings provide evidence of cross-kingdom mycoviral sequence transmission.

Petabase-scale sequence alignment catalyses viral discovery.

Public databases contain a planetary collection of nucleic acid sequences, but their systematic exploration has been inhibited by a lack of efficient methods for searching this corpus, which (at the time of writing) exceeds 20 petabases and is growing exponentially1. Here we developed a cloud computing infrastructure, Serratus, to enable ultra-high-throughput sequence alignment at the petabase scale. We searched 5.7 million biologically diverse samples (10.2 petabases) for the hallmark gene RNA-dependent RNA polymerase and identified well over 105 novel RNA viruses, thereby expanding the number of known species by roughly an order of magnitude. We characterized novel viruses related to coronaviruses, hepatitis delta virus and huge phages, respectively, and analysed their environmental reservoirs. To catalyse the ongoing revolution of viral discovery, we established a free and comprehensive database of these data and tools. Expanding the known sequence diversity of viruses can reveal the evolutionary origins of emerging pathogens and improve pathogen surveillance for the anticipation and mitigation of future pandemics.

Virome of Bat-Infesting Arthropods: Highly Divergent Viruses in Different Vectors.

Bats are reservoirs of important zoonotic viruses like Nipah and SARS viruses. However, whether the blood-sucking arthropods on the body surface of bats also carry these viruses and the relationship between viruses carried by the blood-sucking arthropods and viruses carried by bats have not been reported. This study collected 686 blood-sucking arthropods on the body surface of bats from Yunnan Province, China, between 2012 and 2015, and they included wingless bat flies, bat flies, ticks, mites, and fleas. The viruses carried by these arthropods were analyzed using a meta-transcriptomic approach, and 144 highly diverse positive-sense single-stranded RNA, negative-sense single-stranded RNA, and double-stranded RNA viruses were found, of which 138 were potentially new viruses. These viruses were classified into 14 different virus families or orders, including Bunyavirales, Mononegavirales, Reoviridae, and Picornavirales. Further analyses found that Bunyavirales were the most abundant virus group (84% of total virus RNA) in ticks, whereas narnaviruses were the most abundant (52 to 92%) in the bat flies and wingless bat flies libraries, followed by solemoviruses (1 to 29%) and reoviruses (0 to 43%). These viruses were highly structured based on the arthropod types. It is worth noting that no bat-borne zoonotic viruses were found in the virome of bat-infesting arthropod, seemingly not supporting that bat surface arthropods are vectors of zoonotic viruses carried by bats. IMPORTANCE Bats are reservoirs of many important viral pathogens. To evaluate whether bat-parasitic blood-sucking arthropods participate in the circulation of these important viruses, it is necessary to conduct unbiased virome studies on these arthropods. We evaluated five types of blood-sucking parasitic arthropods on the surface of bats in Yunnan, China, and identified a variety of viruses, some of which had high prevalence and abundance levels, although there is limited overlap in virome between distant arthropods. While most of the virome discovered here is potentially arthropod-specific viruses, we identified three possible arboviruses, including one orthobunyavirus and two vesiculoviruses (family Rhabdoviridae), suggesting bat-parasitic arthropods carry viruses with risk of spillage, which warrants further study.