Mycoviruses as a part of the global virome: Diversity, evolutionary links and lifestyle

Knowledge of mycovirus diversity, evolution, horizontal gene transfer and shared ancestry with viruses infecting distantly related hosts, such as plants and arthropods, has increased vastly during the last few years due to advances in the high throughput sequencing methodologies. This also has enabled the discovery of novel mycoviruses with previously unknown genome types, mainly new positive and negative single-stranded RNA mycoviruses ((+) ssRNA and (-) ssRNA) and single-stranded DNA mycoviruses (ssDNA), and has increased our knowledge of double-stranded RNA mycoviruses (dsRNA), which in the past were thought to be the most common viruses infecting fungi. Fungi and oomycetes (Stramenopila) share similar lifestyles and also have similar viromes. Hypothesis about the origin and cross-kingdom transmission events of viruses have been raised and are supported by phylogenetic analysis and by the discovery of natural exchange of viruses between different hosts during virus-fungus coinfection in planta. In this review we make a compilation of the current information on the genome organization, diversity and taxonomy of mycoviruses, discussing their possible origins. Our focus is in recent findings suggesting the expansion of the host range of many viral taxa previously considered to be exclusively fungal, but we also address factors affecting virus transmissibility and coexistence in single fungal or oomycete isolates, as well as the development of synthetic mycoviruses and their use in investigating mycovirus replication cycles and pathogenicity.

Bunyaviruses Affect Growth, Sporulation, and Elicitin Production in Phytophthora cactorum

Phytophthora cactorum is an important oomycetous plant pathogen with numerous host plant species, including garden strawberry (Fragaria × ananassa) and silver birch (Betula pendula). P. cactorum also hosts mycoviruses, but their phenotypic effects on the host oomycete have not been studied earlier. In the present study, we tested polyethylene glycol (PEG)-induced water stress for virus curing and created an isogenic virus-free isolate for testing viral effects in pair with the original isolate. Phytophthora cactorum bunya-like viruses 1 and 2 (PcBV1 & 2) significantly reduced hyphal growth of the P. cactorum host isolate, as well as sporangia production and size. Transcriptomic and proteomic analyses revealed an increase in the production of elicitins due to bunyavirus infection. However, the presence of bunyaviruses did not seem to alter the pathogenicity of P. cactorum. Virus transmission through anastomosis was unsuccessful in vitro.

High Diversity of Novel Viruses in the Tree Pathogen Phytophthora castaneae Revealed by High-Throughput Sequencing of Total and Small RNA

Phytophthora castaneae, an oomycete pathogen causing root and trunk rot of different tree species in Asia, was shown to harbor a rich diversity of novel viruses from different families. Four P. castaneae isolates collected from Chamaecyparis hodginsii in a semi-natural montane forest site in Vietnam were investigated for viral presence by traditional and next-generation sequencing (NGS) techniques, i.e., double-stranded RNA (dsRNA) extraction and high-throughput sequencing (HTS) of small RNAs (sRNAs) and total RNA. Genome organization, sequence similarity, and phylogenetic analyses indicated that the viruses were related to members of the order Bunyavirales and families Endornaviridae, Megabirnaviridae, Narnaviridae, Totiviridae, and the proposed family “Fusagraviridae.” The study describes six novel viruses: Phytophthora castaneae RNA virus 1–5 (PcaRV1-5) and Phytophthora castaneae negative-stranded RNA virus 1 (PcaNSRV1). All six viruses were detected by sRNA sequencing, which demonstrates an active RNA interference (RNAi) system targeting viruses in P. castaneae. To our knowledge, this is the first report of viruses in P. castaneae and the whole Phytophthora major Clade 5, as well as of the activity of an RNAi mechanism targeting viral genomes among Clade 5 species. PcaRV1 is the first megabirnavirus described in oomycetes and the genus Phytophthora.

Phenotypic Recovery of a Heterobasidion Isolate Infected by a Debilitation-Associated Virus Is Related to Altered Host Gene Expression and Reduced Virus Titer

The fungal genus Heterobasidion includes forest pathogenic species hosting a diverse group of partitiviruses. They include the host debilitating Heterobasidion partitivirus 13 strain an1 (HetPV13-an1), which was originally observed in a slowly growing H. annosum strain 94233. In this study, a relatively fast-growing sector strain 94233-RC3 was isolated from a highly debilitated mycelial culture of 94233, and its gene expression and virus transcript quantities as well as the genomic sequence of HetPV13-an1 were examined. The sequence of HetPV13-an1 genome in 94233-RC3 was identical to that in the original 94233, and thus not the reason for the partial phenotypic recovery. According to RNA-seq analysis, the HetPV13-an1 infected 94233-RC3 transcribed eight genes differently from the partitivirus-free 94233-32D. Three of these genes were downregulated and five upregulated. The number of differentially expressed genes was considerably lower and the changes in their expression were small compared to those of the highly debilitated original strain 94233 with the exception of the most highly upregulated ones, and therefore viral effects on the host transcriptome correlated with the degree of the virus-caused debilitation. The amounts of RdRp and CP transcripts of HetPV13-an1 were considerably lower in 94233-RC3 and also in 94233 strain infected by a closely related mildly debilitating virus HetPV13-an2, suggesting that the virus titer would have a role in determining the effect of HetPV13 viruses on their hosts.

Discovery and Community Dynamics of Novel ssRNA Mycoviruses in the Conifer Pathogen Heterobasidion parviporum

Heterobasidion species are highly destructive basidiomycetous conifer pathogens of the Boreal forest region. Earlier studies have revealed dsRNA virus infections of families Curvulaviridae and Partitiviridae in Heterobasidion strains, and small RNA deep sequencing has also identified infections of Mitoviridae members in these fungi. In this study, the virome of Heterobasidion parviporum was examined for the first time by RNA-Seq using total RNA depleted of rRNA. This method successfully revealed new viruses representing two established (+)ssRNA virus families not found earlier in Heterobasidion: Narnaviridae and Botourmiaviridae. In addition, we identified the presence of a recently described virus group tentatively named “ambiviruses” in H. parviporum. The H. parviporum isolates included in the study originated from experimental forest sites located within 0.7 km range from each other, and a population analysis including 43 isolates was conducted at one of the experimental plots to establish the prevalence of the newly identified viruses in clonally spreading H. parviporum individuals. Our results indicate that viral infections are considerably more diverse and common among Heterobasidion isolates than known earlier and include ssRNA viruses with high prevalence and interspecies variation.

RNA-Dependent RNA Polymerase from Heterobasidion RNA Virus 6 Is an Active Replicase In Vitro

Heterobasidion RNA virus 6 (HetRV6) is a double-stranded (ds)RNA mycovirus and a member of the recently established genus Orthocurvulavirus within the family Orthocurvulaviridae. The purpose of the study was to determine the biochemical requirements for RNA synthesis catalyzed by HetRV6 RNA-dependent RNA polymerase (RdRp). HetRV6 RdRp was expressed in Escherichia coli and isolated to near homogeneity using liquid chromatography. The enzyme activities were studied in vitro using radiolabeled UTP. The HetRV6 RdRp was able to initiate RNA synthesis in a primer-independent manner using both virus-related and heterologous single-stranded (ss)RNA templates, with a polymerization rate of about 46 nt/min under optimal NTP concentration and temperature. NTPs with 2'-fluoro modifications were also accepted as substrates in the HetRV6 RdRp-catalyzed RNA polymerization reaction. HetRV6 RdRp transcribed viral RNA genome via semi-conservative mechanism. Furthermore, the enzyme demonstrated terminal nucleotidyl transferase (TNTase) activity. Presence of Mn²⁺ was required for the HetRV6 RdRp catalyzed enzymatic activities. In summary, our study shows that HetRV6 RdRp is an active replicase in vitro that can be potentially used in biotechnological applications, molecular biology, and biomedicine.

Towards the Forest Virome: High-Throughput Sequencing Drastically Expands Our Understanding on Virosphere in Temperate Forest Ecosystems

Thanks to the development of HTS technologies, a vast amount of genetic information on the virosphere of temperate forests has been gained in the last seven years. To estimate the qualitative/quantitative impact of HTS on forest virology, we have summarized viruses affecting major tree/shrub species and their fungal associates, including fungal plant pathogens, mutualists and saprotrophs. The contribution of HTS methods is extremely significant for forest virology. Reviewed data on viral presence in holobionts allowed us a first attempt to address the role of virome in holobionts. Forest health is dependent on the variability of microorganisms interacting with the host tree/holobiont; symbiotic microbiota and pathogens engage in a permanent interplay, which influences the host. Through virus–virus interplays synergistic or antagonistic relations may evolve, which may drastically affect the health of the holobiont. Novel insights of these interplays may allow practical applications for forest plant protection based on endophytes and mycovirus biocontrol agents. The current analysis is conceived in light of the prospect that novel viruses may initiate an emergent infectious disease and that measures for the avoidance of future outbreaks in forests should be considered.

Viral diversity in Phytophthora cactorum population infecting strawberry

Eighty-eight Phytophthora cactorum strains isolated from crown or leather rot of strawberry in 1971-2019 were screened for viruses using RNA-seq and RT-PCR. Remarkably, all but one isolate were virus-infected, most of them harbouring more than one virus of different genera or species. The most common virus occurring in 94% of the isolates was the Phytophthora cactorum RNA virus 1 (PcRV1) resembling members of Totiviridae. Novel viruses related to members of Endornaviridae, named Phytophthora cactorum alphaendornaviruses 1-3 (PcAEV1-3), were found in 57% of the isolates. Four isolates hosted viruses with affinities to Bunyaviridae, named Phytophthora cactorum bunyaviruses 1-3 (PcBV1-3), and a virus resembling members of the proposed genus 'Ustivirus', named Phytophthora cactorum usti-like virus (PcUV1), was found in a single isolate. Most of the virus species were represented by several distinct strains sharing ≥81.4% aa sequence identity. We found no evidence of spatial differentiation but some temporal changes in the P. cactorum virus community were observed. Some isolates harboured two or more closely related strains of the same virus (PcAEV1 or PcRV1) sharing 86.6%-96.4% nt identity in their polymerase sequence. This was surprising as viruses with such a high similarity are typically mutually exclusive.

The virome from a collection of endomycorrhizal fungi reveals new viral taxa with unprecedented genome organization

Mutualistic plant-associated fungi are recognized as important drivers in plant evolution, diversity, and health. The discovery that mycoviruses can take part and play important roles in symbiotic tripartite interactions has prompted us to study the viromes associated with a collection of ericoid and orchid mycorrhizal (ERM and ORM, respectively) fungi. Our study, based on high-throughput sequencing of transcriptomes (RNAseq) from fungal isolates grown in axenic cultures, revealed in both ERM and ORM fungi the presence of new mycoviruses closely related to already classified virus taxa, but also new viruses that expand the boundaries of characterized RNA virus diversity to previously undescribed evolutionary trajectories. In ERM fungi, we provide first evidence of a bipartite virus, distantly related to narnaviruses, that splits the RNA-dependent RNA polymerase (RdRP) palm domain into two distinct proteins, encoded by each of the two segments. Furthermore, in one isolate of the ORM fungus Tulasnella spp. we detected a 12 kb genomic fragment coding for an RdRP with features of bunyavirus-like RdRPs. However, this 12 kb genomic RNA has the unique features, for Bunyavirales members, of being tri-cistronic and carrying ORFs for the putative RdRP and putative nucleocapsid in ambisense orientation on the same genomic RNA. Finally, a number of ORM fungal isolates harbored a group of ambisense bicistronic viruses with a genomic size of around 5 kb, where we could identify a putative RdRP palm domain that has some features of plus strand RNA viruses; these new viruses may represent a new lineage in the Riboviria, as they could not be reliably assigned to any of the branches in the recently derived monophyletic tree that includes most viruses with an RNA genome.

Complete genome sequence of a novel toti-like virus from the plant-pathogenic oomycete Phytophthora cactorum

This report describes the complete genome sequence of a double-stranded RNA (dsRNA) virus infecting the oomycetous plant pathogen Phytophthora cactorum. The virus genome consists of a single dsRNA segment of 5699 bp with two open reading frames predicted to overlap with each other and encoding a putative capsid protein of 705 aa and an RNA-dependent RNA polymerase of 779 aa. Sequence comparisons indicated that this virus, designated as "Phytophthora cactorum RNA virus 1" (PcRV1), shares the highest sequence similarity with the unclassified Pythium splendens RNA virus 1 (58% RdRp aa sequence identity). Phylogenetic analysis revealed that these two oomycete viruses group together with Giardia lamblia virus (GVL; family Totiviridae) and several unclassified toti-like viruses from arthropods, fish and fungi. This is the first report of a toti-like virus in a member of the genus Phytophthora and the first virus characterized in P. cactorum.

Mixed infection by a partitivirus and a negative-sense RNA virus related to mymonaviruses in the polypore fungus Bondarzewia berkeleyi

Virus communities of forest fungi remain poorly characterized. In this study, we detected two new viruses co-infecting an isolate of the polypore fungus Bondarzewia berkeleyi using high-throughput sequencing. One of them was a putative new partitivirus designated as Bondarzewia berkeleyi partitivirus 1 (BbPV1), with two linear dsRNA genome segments of 1928 and 1863 bp encoding a putative RNA-dependent RNA polymerase (RdRP) of 591 aa and a putative capsid protein of 538 aa. The other virus, designated as Bondarzewia berkeleyi negative-strand RNA virus 1 (BbNSRV1), had a non-segmented negative-sense RNA genome of 10,983 nt and was related to members of family Mymonaviridae. The BbNSRV1 genome includes six predicted open reading frames (ORFs) of 279, 425, 230, 174, 200 and 1970 aa. The longest ORF contained conserved regions corresponding to Mononegavirales RdRP and mRNA-capping enzyme region V constituting the mononegavirus Large protein. In addition, a low level of sequence identity was detected between the putative nucleocapsid protein-coding ORF2 of Lentinula edodes negative-strand RNA virus 1 and BbNSRV1. The viruses characterized in this study are the first ones described in Bondarzewia spp., and BbNSRV1 is the second mymona-like virus described in a basidiomycete host.

Virus population structure in the ectomycorrhizal fungi Lactarius rufus and L. tabidus at two forest sites in Southern Finland

Lactarius fungi belong to the Russulaceae family and have an important ecological role as ectomycorrhizal symbionts of coniferous and deciduous trees. Two Lactarius species, L. tabidus and L. rufus have been shown to harbor bisegmented dsRNA viruses belonging to an unclassified virus group including the mutualistic Curvularia thermal tolerance virus (CThTV). In this study, we characterized the first complete genome sequences of these viruses designated as Lactarius tabidus RNA virus 1 (LtRV1) and Lactarius rufus RNA virus 1 (LrRV1), both of which included two genome segments of 2241 and 2049 bp. We also analyzed spatial distribution and sequence diversity of the viruses in sixty host strains at two forest sites, and showed that the viruses are species-specific at sites where both host species co-occur. We also found that single virus isolates inhabited several different conspecific host strains, and were involved in persistent infections during up to eight years.

Mitoviruses in the conifer root rot pathogens Heterobasidion annosum and H. parviporum

Mitoviral infections are highly common among fungi, but so far only one mitovirus has been described in Heterobasidion spp. conifer pathogens. Here, the occurrence of further mitoviruses was investigated using a previously published RNA-Seq dataset for de novo contig assembly. This allowed the identification of two additional mitovirus strains designated as Heterobasidion mitovirus 2 (HetMV2) and HetMV3 with genome lengths of ca. 2.9 and 5.0 kb. Furthermore, the occurrence of similar viruses was screened among a collection of Heterobasidion isolates using RT-PCR. Mitoviruses were detected in six more fungal isolates and two different host species, H. annosum and H. parviporum.

Viruses of fungi and oomycetes in the soil environment

Soils support a myriad of organisms hosting highly diverse viromes. In this minireview, we focus on viruses hosted by true fungi and oomycetes (members of Stamenopila, Chromalveolata) inhabiting bulk soil, rhizosphere and litter layer, and representing different ecological guilds, including fungal saprotrophs, mycorrhizal fungi, mutualistic endophytes and pathogens. Viruses infecting fungi and oomycetes are characterized by persistent intracellular nonlytic lifestyles and transmission via spores and/or hyphal contacts. Almost all fungal and oomycete viruses have genomes composed of single-stranded or double-stranded RNA, and recent studies have revealed numerous novel viruses representing yet unclassified family-level groups. Depending on the virus–host combination, infections can be asymptomatic, beneficial or detrimental to the host. Thus, mycovirus infections may contribute to the multiplex interactions of hosts, therefore likely affecting the dynamics of fungal communities required for the functioning of soil ecosystems. However, the effects of fungal and oomycete viruses on soil ecological processes are still mostly unknown. Interestingly, new metagenomics data suggest an extensive level of horizontal virus transfer between plants, fungi and insects.

Alphapartitiviruses of Heterobasidion Wood Decay Fungi Affect Each Other's Transmission and Host Growth

Heterobasidion spp. root rot fungi are highly destructive forest pathogens of the northern boreal forests, and are known to host a diverse community of partitiviruses. The transmission of these mycoviruses occurs horizontally among host strains via mycelial anastomoses. We revealed using dual cultures that virus transmission rates are affected by pre-existing virus infections among two strains of H. annosum. The transmission efficacy of mycovirus HetPV15-pa1 to a pre-infected host was elevated from zero to 50% by the presence of HetPV13-an1, and a double infection of these viruses in the donor resulted in an overall transmission rate of 90% to a partitivirus-free recipient. On contrary, pre-existing virus infections of two closely related strains of HetPV11 hindered each other's transmission, but had unexpectedly dissimilar effects on the transmission of more distantly related viruses. The co-infection of HetPV13-an1 and HetPV15-pa1 significantly reduced host growth, whereas double infections including HetPV11 strains had variable effects. Moreover, the results showed that RdRp transcripts are generally more abundant than capsid protein (CP) transcripts and the four different virus strains express unique transcripts ratios of RdRp and CP. Taken together, the results show that the interplay between co-infecting viruses and their host is extremely complex and highly unpredictable.

Virus Detection by High-Throughput Sequencing of Small RNAs: Large-Scale Performance Testing of Sequence Analysis Strategies

Recent developments in high-throughput sequencing (HTS), also called next-generation sequencing (NGS), technologies and bioinformatics have drastically changed research on viral pathogens and spurred growing interest in the field of virus diagnostics. However, the reliability of HTS-based virus detection protocols must be evaluated before adopting them for diagnostics. Many different bioinformatics algorithms aimed at detecting viruses in HTS data have been reported but little attention has been paid thus far to their sensitivity and reliability for diagnostic purposes. Therefore, we compared the ability of 21 plant virology laboratories, each employing a different bioinformatics pipeline, to detect 12 plant viruses through a double-blind large-scale performance test using 10 datasets of 21- to 24-nucleotide small RNA (sRNA) sequences from three different infected plants. The sensitivity of virus detection ranged between 35 and 100% among participants, with a marked negative effect when sequence depth decreased. The false-positive detection rate was very low and mainly related to the identification of host genome-integrated viral sequences or misinterpretation of the results. Reproducibility was high (91.6%). This work revealed the key influence of bioinformatics strategies for the sensitive detection of viruses in HTS sRNA datasets and, more specifically (i) the difficulty in detecting viral agents when they are novel or their sRNA abundance is low, (ii) the influence of key parameters at both assembly and annotation steps, (iii) the importance of completeness of reference sequence databases, and (iv) the significant level of scientific expertise needed when interpreting pipeline results. Overall, this work underlines key parameters and proposes recommendations for reliable sRNA-based detection of known and unknown viruses.

Distribution of Viruses Inhabiting Heterobasidion annosum in a Pine-Dominated Forest Plot in Southern Finland

We investigated the diversity and spatial distribution of viruses infecting strains of the root rot fungus Heterobasidion annosum collected from pine stumps at a heavily infected forest site. Four different partitiviruses were detected in 14 H. annosum isolates at the study site, constituting approximately 29% of all Heterobasidion isolates investigated (N = 48). Two of the viruses detected were new partitiviruses designated here as Heterobasidion partitivirus 16 (HetPV16) and HetPV20, and two were previously known partitiviruses: HetPV7 and HetPV13. The two new partitiviruses found, HetPV16-an1 and HetPV20-an1, shared ~70% RdRp nucleotide sequence identity with the alphapartitivirus Rosellinia necatrix partitivirus 2, and less than 40% identity with known viruses of Heterobasidion spp. HetPV7-an1 was closely similar to HetPV7-pa1 isolated earlier from Heterobasidion parviporum, supporting the view of conspecific virus pools in different Heterobasidion species. Three fungal isolates were found to be co-infected with two different partitivirus strains (HetPV7-an1 and HetPV13-an2 or HetPV16-an1 and HetPV20-an1). Different isolates representing each host clone had variable virus compositions, and virus strains occurring in more than one host clone showed minor sequence variations between clones.

Intraspecific comparative genomics of isolates of the Norway spruce pathogen (Heterobasidion parviporum) and identification of its potential virulence factors

BACKGROUND

Heterobasidion parviporum is an economically most important fungal forest pathogen in northern Europe, causing root and butt rot disease of Norway spruce (Picea abies (L.) Karst.). The mechanisms underlying the pathogenesis and virulence of this species remain elusive. No reference genome to facilitate functional analysis is available for this species.

RESULTS

To better understand the virulence factor at both phenotypic and genomic level, we characterized 15 H. parviporum isolates originating from different locations across Finland for virulence, vegetative growth, sporulation and saprotrophic wood decay. Wood decay capability and latitude of fungal origins exerted interactive effects on their virulence and appeared important for H. parviporum virulence. We sequenced the most virulent isolate, the first full genome sequences of H. parviporum as a reference genome, and re-sequenced the remaining 14 H. parviporum isolates. Genome-wide alignments and intrinsic polymorphism analysis showed that these isolates exhibited overall high genomic similarity with an average of at least 96% nucleotide identity when compared to the reference, yet had remarkable intra-specific level of polymorphism with a bias for CpG to TpG mutations. Reads mapping coverage analysis enabled the classification of all predicted genes into five groups and uncovered two genomic regions exclusively present in the reference with putative contribution to its higher virulence. Genes enriched for copy number variations (deletions and duplications) and nucleotide polymorphism were involved in oxidation-reduction processes and encoding domains relevant to transcription factors. Some secreted protein coding genes based on the genome-wide selection pressure, or the presence of variants were proposed as potential virulence candidates.

CONCLUSION

Our study reported on the first reference genome sequence for this Norway spruce pathogen (H. parviporum). Comparative genomics analysis gave insight into the overall genomic variation among this fungal species and also facilitated the identification of several secreted protein coding genes as putative virulence factors for the further functional analysis. We also analyzed and identified phenotypic traits potentially linked to its virulence.

Characterization of small RNAs originating from mitoviruses infecting the conifer pathogen Fusarium circinatum

Deep sequencing of small RNAs has proved effective in the diagnosis of mycovirus infections. In this study, the presence of mycoviruses in ten isolates of the phytopathogenic fungus Fusarium circinatum was investigated by high-throughput sequencing (HTS) of small RNAs. The contigs resulting from de novo assembly of the reads were aligned to viral genome sequences. The presence of each mycovirus detected in the isolates was confirmed by RT-PCR analysis with four previously described primer pairs and seven new pairs designed on the basis of sequencing data. The findings demonstrate the potential use of HTS for reconstructing previously identified mitoviruses infecting F. circinatum.

ICTV Virus Taxonomy Profile: Partitiviridae

The Partitiviridae is a family of small, isometric, non-enveloped viruses with bisegmented double-stranded (ds) RNA genomes of 3–4.8 kbp. The two genome segments are individually encapsidated. The family has five genera, with characteristic hosts for members of each genus: either plants or fungi for genera Alphapartitivirus and Betapartitivirus, fungi for genus Gammapartitivirus, plants for genus Deltapartitivirus and protozoa for genus Cryspovirus. Partitiviruses are transmitted intracellularly via seeds (plants), oocysts (protozoa) or hyphal anastomosis, cell division and sporogenesis (fungi); there are no known natural vectors. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the taxonomy of the Partitiviridae, which is available at www.ictv.global/report/partitiviridae.

A Framework for the Evaluation of Biosecurity, Commercial, Regulatory, and Scientific Impacts of Plant Viruses and Viroids Identified by NGS Technologies

Recent advances in high-throughput sequencing technologies and bioinformatics have generated huge new opportunities for discovering and diagnosing plant viruses and viroids. Plant virology has undoubtedly benefited from these new methodologies, but at the same time, faces now substantial bottlenecks, namely the biological characterization of the newly discovered viruses and the analysis of their impact at the biosecurity, commercial, regulatory, and scientific levels. This paper proposes a scaled and progressive scientific framework for efficient biological characterization and risk assessment when a previously known or a new plant virus is detected by next generation sequencing (NGS) technologies. Four case studies are also presented to illustrate the need for such a framework, and to discuss the scenarios.

Taxonomy, biogeography and importance of Heterobasidion viruses

The genus Heterobasidion consists of several species of necrotrophic and saprotrophic fungi, and includes some of the most detrimental organisms in boreal conifer forests. These fungi host a widespread and diverse mycovirus community composed of more than 16 species of Partitiviridae, a species of Narnaviridae and one taxonomically unassigned virus related to the Curvularia thermal tolerance virus. These viruses are able to cross species borders, co-infect single host strains and cause phenotypic changes in their hosts. The abundance of viruses increases over time in Heterobasidion infection centers, and they are targeted by fungal RNA interference. Long-term field studies are essential for obtaining a comprehensive view of virus effects in the nature.

Description and prevalence of a putative novel mycovirus within the conifer pathogen Gremmeniella abietina

The European race of Gremmeniella abietina (Lagerberg) Morelet is the causal agent of stem canker and shoot blight on numerous conifers in Europe and North America. It comprises different species and biotypes in which the presence of mycoviruses has been determined. In this report, we describe the full-length sequence of the RNA-dependent RNA polymerase (RdRp) of a putative novel virus, Gremmeniella abietina RNA virus 6 (GaRV6), with 2165 nt and a GC content of 54.7 %. A BLASTp search using the deduced RdRp amino acid sequence confirmed GaRV6 to be related to members of a still unassigned virus taxon, which includes, e.g., Fusarium graminearum dsRNA mycovirus 4 (FgV-4) and the mutualistic Curvularia thermal tolerance virus (CThTV). The prevalence and genetic diversity of GaRV6 was also studied within the European race of G. abietina. We examined 162 isolates originating from Canada, the Czech Republic, Finland, Italy, Montenegro, Serbia, Spain, Switzerland, Turkey and the United States. According to direct specific reverse transcription (RT) PCR screening based on the RdRp sequence, the virus appears to be present only in Spain, where it is relatively abundant but genetically highly uniform.

Viruses accumulate in aging infection centers of a fungal forest pathogen

Fungal viruses (mycoviruses) with RNA genomes are believed to lack extracellular infective particles. These viruses are transmitted laterally among fungal strains through mycelial anastomoses or vertically via their infected spores, but little is known regarding their prevalence and patterns of dispersal under natural conditions. Here, we examined, in detail, the spatial and temporal changes in a mycovirus community and its host fungus Heterobasidion parviporum, the most devastating fungal pathogen of conifers in the Boreal forest region. During the 7-year sampling period, viruses accumulated in clonal host individuals as a result of indigenous viruses spreading within and between clones as well as novel strains arriving via airborne spores. Viral community changes produced pockets of heterogeneity within large H. parviporum clones. The appearance of novel viral infections in aging clones indicated that transient cell-to-cell contacts between Heterobasidion strains are likely to occur more frequently than what was inferred from genotypic analyses. Intraspecific variation was low among the three partitivirus species at the study site, whereas the unassigned viral species HetRV6 was highly polymorphic. The accumulation of point mutations during persistent infections resulted in viral diversification, that is, the presence of nearly identical viral sequence variants within single clones. Our results also suggest that co-infections by distantly related viral species are more stable than those between conspecific strains, and mutual exclusion may play a role in determining mycoviral communities.

Diagnosis and discovery of fungal viruses using deep sequencing of small RNAs

Analysis of virus-derived small RNAs with high-throughput sequencing has been successful for detecting novel viruses in plants and invertebrates. However, the applicability of this method has not been demonstrated in fungi, although fungi were among the first organisms reported to utilize RNA silencing. Here, we used virus-infected isolates of the fungal species complex Heterobasidion annosum sensu lato as a model system to test whether mycovirus genome segments can be detected with small RNA deep sequencing. Species of the genus Heterobasidion are some of the most devastating forest pathogens in boreal forests. These fungi cause wood decay and are commonly infected with species of the family Partitiviridae and the unassigned virus species Heterobasidion RNA virus 6. Small RNA deep sequencing allowed the simultaneous detection of all eight double-stranded RNA virus strains known to be present in the tested samples and one putative mitovirus species (family Narnaviridae) with a single-stranded RNA genome, designated here as Heterobasidion mitovirus 1. Prior to this study, no members of the family Narnaviridae had been described as infecting species of Heterobasidion. Quantification of viral double- and single-stranded RNA with quantitative PCR indicated that co-infecting viral species and viruses with segmented genomes can be detected with small RNA deep sequencing despite vast differences in the amount of RNA. This is the first study demonstrating the usefulness of this method for detecting fungal viruses. Moreover, the results suggest that viral genomes are processed into small RNAs by different species of Heterobasidion.

Molecular characterization of HetRV8-ir1, a partitivirus of the invasive conifer pathogenic fungus Heterobasidion irregulare

The Heterobasidion annosum (Fr.) Bref. complex includes some of the most destructive conifer pathogenic fungi in the Boreal forest region. H. irregulare, formerly known as the North American pine type of H. annosum, was introduced from North America into Italy during the Second World War and occurs as an invasive pathogen in Pinus pinea stands together with the native European species H. annosum sensu stricto. We describe the complete nucleotide sequence of a new putative partitivirus from an Italian strain of H. irregulare. The bisegmented genome of HetRV8-ir1 encodes an RNA-dependent RNA polymerase of 704 aa and a capsid protein of 638 aa. The polymerase and capsid aa sequences are relatively similar (59-78 %) to those of Fusarium poae virus 1, Pleurotus ostreatus virus 1, and grapevine-associated partitivirus 1. HetRV8-ir1 is the first virus described from H. irregulare, and it is distantly related to previously known partitiviruses of Heterobasidion species.

Virus community dynamics in the conifer pathogenic fungus Heterobasidion parviporum following an artificial introduction of a partitivirus

Viruses infecting the conifer pathogenic fungus Heterobasidion annosum sensu lato are intracellular and spread via anastomosis contacts. In the laboratory, these viruses transmit readily even between somatically incompatible isolates, but their dispersal capacity in natural conditions has not been previously studied. We introduced a mycovirus to a heavily diseased forest site by inoculating Norway spruce stumps with heartrot decay using a mycelial suspension of Heterobasidion parviporum strain RT3.49C hosting the partitivirus strain HetRV4-pa1. The Heterobasidion population at the sample plot was screened for mycoviruses prior to and after the inoculation. Based on sequence analysis, the resident H. parviporum strains harbored six different strains of the virus species Heterobasidion RNA virus 6 (HetRV6) and one strain of HetRV4 prior to the inoculation. After three growth seasons, the inoculated H. parviporum host strain was not detected, but the introduced virus had infected two resident H. parviporum genets. The presence of a preexisting HetRV6 infection did not hinder spread of the introduced partitivirus but resulted in coinfections instead. The resident HetRV6 virus population seemed to be highly stable during the incubation period, while the single indigenous HetRV4 infection was not detected after the inoculation. In laboratory infection experiments, the introduced virus could be transmitted successfully into all of the resident H. parviporum genets. This study shows for the first time transmission of a Heterobasidion virus between somatically incompatible hosts in natural conditions.

A novel putative partitivirus of the saprotrophic fungus Heterobasidion ecrustosum infects pathogenic species of the Heterobasidion annosum complex

We characterized the bisegmented genome of a putative double-stranded (ds) RNA virus from a Chinese isolate of the fungus Heterobasidion ecrustosum, a member of the Heterobasidion insulare species complex. The larger genomic segment of 1885bp encoded a putative RNA dependent RNA polymerase (RdRp, 585aa), and the smaller one for a putative coat protein of 521aa (1826bp). Phylogenetic analyses suggest that this novel virus species, named as 'Heterobasidion RNA virus 3 from H. ecrustosum, strain 1' (HetRV3-ec1), can be assigned to the family Partitiviridae, being most similar to the Helicobasidium mompa dsRNA mycovirus with RdRp amino acid similarity of 54%. The similarity to known viruses of other Heterobasidion species was notably low (25-39%). The virus could be experimentally transmitted to members of the Heterobasidion annosum complex: the European Heterobasidion abietinum and North American Heterobasidion occidentale, and the original host strain could be cured from the virus by thermal treatment. Microscopical observations showed that hyphae of H. ecrustosum anastomosed occasionally with H. abietinum and H. occidentale, and suggested a possible route for horizontal transmission between these sexually incompatible species. The virus infection seemed to cause variable effects on the growth rate of its fungal hosts, but the results were strongly dependent on fungal strain, growth medium and incubation temperature.

Direct analysis of ribosomal DNA in denaturing gradients: application on the effects of Phlebiopsis gigantea treatment on fungal communities of conifer stumps

The aim of this study was to test the usefulness of direct PCR-amplification in analysing fungal diversity in stumps. The analysis was conducted on stumps treated against Heterobasidion spp. using a commercial formulation of Phlebiopsis gigantea (Rotstop), and carried out using denaturing gradient gel electrophoresis (DGGE) of small subunit (SSU) ribosomal DNA (rDNA) fragments PCR-amplified directly from wood DNA samples using two separate fungus-specific primer pairs. On average, two (range 0-9) different amplification products were observed by DGGE in single wood samples of approximately 500 mm3. The PCR products were classified into operational taxonomical unit (OTU) groups based on their DGGE mobility. Six OTUs could be affiliated with a known species based on a reference fungal collection of 37 species: Heterobasidion annosum, H. parviporum. Hypholoma capnoides, P. gigantea, Resinicium bicolor and Stereum sanguinolentum. Sequence analyses did not give further identifications. OTU profiles from old (6 yr-old) and fresh (1-year-old) Scots pine and Norway spruce stumps from treated and untreated forest plots were compared statistically, and some significant differences were observed in the species composition between the treated and untreated plots. However, the frequency of most of the OTUs did not seem to be affected, and the treatment did not seem to have reduced the overall level of fungal diversity. Based on these results, direct PCR-amplification seems to be useful in analyses of fungal communities in decaying conifer stumps.

Comparison of partial 16S rRNA gene sequences obtained from activated sludge bacteria

The cultivated and uncultivated bacterial communities of an activated sludge plant were studied. Two samples were taken and a total of 516 bacterial isolates were classified into groups using their whole-cell protein patterns. The distribution of bacteria into protein-pattern groups differed significantly between the two samples, suggesting variation in culturable bacterial flora. Partial 16S rRNA gene sequences were determined for representatives of the commonest protein-pattern groups. Most of the sequences obtained were previously unknown, but relatively closely related to known sequences of organisms belonging to the alpha, beta or gamma subclasses of the proteobacteria, the first two subclasses being predominant. This classification of bacteria isolated on a diluted nutrient-rich medium differed from recent culture-dependent studies using nutrient-rich media. The uncultivated bacterial community was studied by analyzing ten partial 16S rRNA gene sequences cloned directly from activated sludge. None of the cloned sequences was identical to those determined for culturable organisms; or to those in the GenBank database. They were, however, related to the alpha or beta subclasses of the proteobacteria, or to the gram-positive bacteria with a high G + C DNA content.