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.

The majority of mycoviruses lead to asymptomatic infections in their hosts . However, some mycoviruses induce hypovirulence in phytopathogenic fungi, as exemplified by Chryphonectria hypovirus 1 (CHV1, a hypovirus), which serves as a biological control agent against the destructive chestnut blight (Rigling and Prospero, 2018). An increasing number of mycoviruses are now known to induce phenotypic alterations such as decrease in vitro fungal growth and/or virulence, when their hosts are pathogenic to higher organisms Okada et al., 2018;Garcia-Pedrajas et al., 2019). A unique example is an ssDNA virus that has been shown to alter the host fungus Sclerotinia sclerotiorum from a parasite to an endophyte (Zhang et al., 2020). Moreover, the enhancement of in vitro growth, sporulation and/or virulence by mycoviruses has been reported (Kotta-Loizou and Coutts, 2017;Okada et al., 2018;Shah et al., 2020). However, investigation of phenotypic effects by mycoviruses on their host fungi is often hampered by coinfections that are common in various fungi (Hillman et al., 2018). Several types of virus/virus interactions, e.g., synergistic, antagonistic, and mutualistic interactions, have been reported in fungal hosts (Eusebio-Cope and Suzuki, 2015;Hillman et al., 2018). Coinfections of single fungal strains by over ten mycoviruses have been reported for several phytopathogenic fungi. Examples include a strain of Fusarium poae coinfected with 16 RNA mycoviruses belonging to 11 viral families (Osaki et al., 2016), and a strain of Kickxella alabastrina (the subdivision of Kickxellomycotina) co-infected with 11 RNA mycoviruses (Myers et al., 2020). How these co-infecting mycoviruses interplay remains unknown in these cases. We also have screened many different Pakistani fungal strains from different sources for mycovirus hunting and have reported molecular characterization of some of the discovered novel mycoviruses (Khan et al., 2021a;Khan et al., 2022). The tested fungal strains included the SP1 strain of Fusarium mangiferae that was found to be co-infected with 11 mycoviruses (Khan et al., 2021b).
The ascomycetous family Botryosphaeriaceae includes many important phyotpathogenic fungi such as members of the genera Botryosphaeria, Neofusicoccum, and Diploidia. Whereas some members of the first two genera have been explored as viral hosts (Wang et al., 2014;Marais et al., 2021), the third one has not. Diplodia seriata (anamorph of Botryosphaeria obtusa) is a phytopathogenic ascomycete that causes fruits rot, leaf spot, canker, leaf chlorosis and trunk dieback in many woody and herbaceous plants (Phillips et al., 2007;Farr and Rossman, 2022). In the current study, as an extension of our previous studies, we characterized the virome of a Pakistani strain (L3) of D. seriata. The fungal strain was shown to be co-infected by a total of eight RNA mycoviruses spanning four dsRNA virus families and two (+)ssRNA virus families with encapsidated and putative capsidless nature, some of which accommodate members with new virus lifestyles. Hyphal tipping resulted in the elimination of three viruses from the original fungal strain and the restoration of the colony growth.

Isolation of Fungal Strains and Extraction of dsRNA/Total RNA Fractions
The fungal strain L3 (D. seriata) was collected and identified during a mycoviral screening survey conducted in 2018 from a diseased leaf of the Lokath plant (Khan et al., 2021b). Some mycoviruses identified as a result of this screening survey have already been published (Khan et al., 2021a;Khan et al., 2021b;Khan et al., 2022). For isolating purified cultures of L3, surface sterilization of leaves was performed with 1% sodium hypochlorite solution and further washing was done using autoclaved sterile distilled water. Samples were air-dried in a safety cabinet and infected parts were inoculated on a fresh potato dextrose agar (PDA, Becton, Dickinson and Co.) plate to obtain fungal cultures. Incubation of fungal cultures was performed at 25°C for 4 to 8 days for growth purposes. Glycerol stocks (40% v/v) of the pure colonies were prepared and stored at -80°C for future use. Fungal identification was performed based on the morphology ( Figure 1) and internal transcribed spacer (ITS) region sequencing (Supplementary Figure 1) (White et al., 1990). Total RNA fractions were obtained as described earlier (Eusebio-Cope and Suzuki, 2015). The classical dsRNA extraction method using cellulose chromatography column was employed to extract dsRNA fractions from the fungal samples (Bhatti et al., 2012;Eusebio-Cope and Suzuki, 2015). The obtained dsRNA fractions were then treated with RQ1 DNase (Promega Corp.) and S1 Nuclease (Thermo Fisher Scientific Inc.) to eliminate genomic DNA and ribosomal RNA contaminations. under a dissecting microscope in a clean bench, single hyphal tips from single hyphae were excised with a scalpel and were allowed to grow on new water agar plates. Single sub-isolates were obtained after four rounds of hyphal tipping. All regenerated isolates were cultured on cellophane overlaid PDA to examine for the presence of mycoviruses via dsRNA extraction or RT-PCR assay. A total of 30 independent hyphal tip isolates were obtained. Phenotype of the parental strain (L3) and a partially cured hyphal tip subculture (L3ht1) were compared. A C B FIGURE 1 | Diverse mycoviruses were detected from Diplodia seriata (the strain L3). (A) Colony morphology of the strain L3 and L3ht1 grown of PDA media. L3 was partially cured by hyphal tipping to obtain L3ht1. (B) Agarose gel electrophoresis of dsRNA and total nucleic acid fractions. DsRNA was purified from an equal amount (5 mg) of total RNA fractions of the two strains and resuspended in an equal amount of distilled water. The genomic dsRNA of MyRV1 (Suzuki et al., 2004) was used as the size standard. (C) RT-PCR analysis of the two fungal strains L3 and L3ht1. Total RNA was isolated from the strains and subjected to RT-PCR using primer sets shown in Supplementary Table S1 Next-Generation Sequencing of Viral RNA Total RNA fractions were separately obtained from three dsRNApositive fungal isolates, L3 (D. seriata, a phytopathogen), SP1 (F. mangiferae, a phytopathogen) and HK (Geotrichum candidum, a soil-inhabitant) and were pooled together as already described by Khan et al. (Khan et al., 2021a;Khan et al., 2021b). Ribosomal RNA depletion, cDNA library construction and RNA sequencing were performed by Macrogen Inc. (Tokyo, Japan) as already described by Khan et al. (Khan et al., 2021b). The assembly of RNA sequencing reads was performed using the Read Mapping algorithm with a default parameter in the CLC Genomics Workbench (version 11, CLC Bio-Qiagen) assembled 42,770 contig sequences were subjected to local BLAST searches to check sequence similarities with deposited virus reference sequences (RefSeq) in the National Center for Biotechnology Information (NCBI) database (data not shown). The presence of mycoviruses corresponding to virus-like contigs was confirmed via RT-PCR amplification from total RNA fractions of the L3 strain with virus specific primers (Supplementary Table S1).
The terminal sequences of viral genomes or segments were achieved through the 3′ RNA ligase mediated rapid amplification of complementary cDNA ends (RLM-RACE) procedure as described by Suzuki et al. (Suzuki et al., 2004). RT-PCR and RACE amplicons were subjected to Sanger sequencing.  -Bauer et al., 2017). Conservation at 5′-and 3′-terminal sequences of the genomes or segments and the motif prediction (I-VIII) for RNA-dependent RNA polymerase (RdRP) core domains were achieved by conducting nucleotide and amino acid sequence alignments, respectively, t h r o u g h t h e U G E N E p r o g r a m ( h t t p : / / u g e n e . n e t / ) (Okonechnikov et al., 2012). The viral sequences were deposited in the GenBank/ENA/DDBJ database. GC content (%) was calculated through the online server of ENDMEMO program (http://www.endmemo.com/bio/gc.php). The predicted molecular mass of the putative protein encoded by viral genomes or their segments was calculated by an online protein molecular weight estimation program (https://www.bioinformatics.org/ sms/prot_mw.html). Transmembrane domains were predicted by TMHMM -2.0 (https://services.healthtech.dtu.dk/service. php?TMHMM-2.0).

Bioinformatic Analysis
Phylogenetic analyses were performed using RdRP deduced amino acid sequences. Multiple sequence alignments were constructed using the online server of MAFFT (version 7) (https://mafft.cbrc.jp/alignment/server/) (Katoh and Toh, 2008). Unreliable regions of alignments were trimmed using Gblocks ver. 0.91b (Talavera and Castresana, 2007). Maximum likelihood (ML) trees were constructed using the online version (3.0) of PhyML with automatic model selection by SMS (Smart Model Selection) (Lefort et al., 2017) (http://www.atgcmontpellier.fr/phyml-sms/). Bootstrap values with 1000 replicates were used for the analysis and values less than 500 were masked. Trees were rooted at the mid-point for defining the out-group. Final ML trees were visualized using the iTOLL; Interactive Tree Of Life online server (https://itol.embl.de/ upload.cgi).

RESULTS AND DISCUSSION
The current study describes the molecular characterization of diverse mycoviruses from a single isolate of D. seriata (strain L3) ( Figure 1 and Supplementary Figure S1), which was identified as dsRNA positive in the previous screening of Pakistani fungal collections (Khan et al., 2021b). The agarose gel profile of dsRNA extracted from L3 showed the presence of multiple segments with sizes varying from 1.2 to 6.0 kbp ( Figure 1). We first used the RNA sequencing approach to obtain genomic sequences of mycoviruses harbored by strain L3. NGS sequencing provided many virus-like sequence contigs. Because the NGS data were from three fungal strains (see Materials and Methods), sequence contigs associated with previously reported mycoviruses derived from two other fungi (Khan et al., 2021a;Khan et al., 2021b) were eliminated from the current analyses. RT-PCR analyses assigned a total of 17 virus-like sequence contigs to strain L3. Their RT-PCR pattern is shown in Figure 1. Nucleotide sequence and molecular phylogenetic analyses confirmed the presence of eight putative novel mycoviruses belonging to seven viral families (four dsRNA and four (+) ssRNA viruses, respectively) (see below). Identified mycoviruses belonged to the families Chrysoviridae, Polymycoviridae, Totiviridae, Partitiviridae, Botourmiaviridae, as well as proposed families "Ambiguiviridae" and "Splipalmiviridae" ( Table 1). The complete genome sequences were determined for almost all these viruses, while a few genomic segments have yet to be completed and are in the coding-complete status. The accession numbers of the newly determined sequences are presented in Table 1.
Phylogenetic analysis, based on the complete aa sequence of the RdRP of DsCV1 and selected members of the Chrysoviridae family indicated that DsCV1 was phylogenetically close to MpCV1 and clustered with alphachrysoviruses ( Figure 2). These results are in accordance with the BLASTp search result ( Table 1). According to ICTV species demarcation criteria for chrysoviruses (≤70 % and ≤53 % identity in RdRP and CP) (Kotta-Loizou et al., 2020), DsCV1 appears to represent a new species of the genus Alphachrysovirus, the first chrysovirus infecting the phytopathogenic fungus D. seriata. Since CnCV1 is likely to have a very narrow host range (Shahi et al., 2021), it is worth interesting whether DsCV1 and other chrysoviruses have narrow or broad host ranges.
Although the experimental data are not available, the presence of this domain suggests that P3 is involved in capping genomic dsRNA or messenger RNAs as proposed earlier (Kanhayuwa et al., 2015). DsPmV1 dsRNA4 encodes a protein of 275 aa (designated P4, 28.6 kDa) (Figure 3). This segment was predicted to contain a proline-alanine-serine rich protein (PASrp) [P(21), A(41), and S(21) residues] and showed 77.0% similarity with P4 of BdRV1 (Table 1). Phylogenetic analysis, based on the complete aa sequence of the RdRPs, suggested that DsPmV1 shows close phylogenetic affinity to BdRV1 in the family Polymycoviridae (Figure 3). According to polymycovirus species demarcation criteria (≤70 % aa sequence identity in the RdRP) (https://talk.ictvonline.org/ files/ictv_official_taxonomy_updates_since_the_8th_report/m/ fungal-official/9407), DsPmV1 should be classified as the same species as BdRV1 based on the high level of aa sequence identity. It will be an interesting question whether these two viruses might have co-evolved with their host fungi or whether they might have been horizontally transferred between the two phylogenetically closely related fungi.
An RdRP-based phylogenetic tree showed that DsPV1 belongs to the newly proposed genus "Epsilonpartitivirus" in the family Partitiviridae together with MbPV1, CePV1, PaPV1 and many others (Nerva et al., 2017) (Figure 4). The ICTV-set species threshold for partitiviruses is ≤ 90% aa sequence identity in RdRP and ≤ 80% aa sequence identity in CP (Vainio et al., 2018). Based on this criterion, DsPV1 is considered to represent a new species in this proposed genus. CP and RdRP shared 85.2% and 81.6% aa sequence identity with those of Sphaeropsis sapinea RNA virus 2 (SsRV2) from a pine pathogen Diplodia pinea (Preisig et al., 1998) (Table 1). A multiple nucleotide sequence alignment showed that DsVV1 shared conservation at both termini with reported vivtoriviruses, Magnaporthe oryzae virus 2 (MoV2; accession no. LC573906) and Sphaeropsis sapinea RNA virus 2 (SsRV2; accession no. AF039080) (Supplementary Figure S5A). NCBI CDD search showed that ORF1 (CP) and ORF2 (RdRP) were predicted to have a conserved domain of totivirus CP (pfam05518, spanning from 15 aa to 707 aa) and RdRP domains (pfam02123, spanning from 101 to 574 aa), respectively (Supplementary Figure S5B). At the junction of ORF1 and ORF2, a conserved tetranucleotide sequence 5′-AUGA-3′ (a potential re-initiation start codon is underlined) was observed ( Figure 5). This tetranucleotide could be involved in translational termination and re-initiation of the two ORFs, which can express two proteins CP and RdRP separately, as for those of Helminthosporium victoriae 190S virus (Hv190SV) and some other victroiviruses (Li et al., 2015). A pseudoknot structure important for stop/re-initiation translation is predicted to form directly upstream of the tetranucleotide sequence 5′-AUGA-3′ at a position similar to that in Hv190SV (Li et al., 2015). As a putative stop/restart facilitator, some victoriviruses and hypoviruses have a pentanucleotide sequence at their junction (5'-UAAUG-3') (Guo et al., 2009;Chiba et al., 2013;Jamal et al., 2019).

A Novel Totivirid
Phylogenetic analysis, based on the complete aa sequence of the RdRp (ORF2) of DsVV1 and selected members of the Totiviridae family indicated that DsVV1 made clade with reported victoriviruses, though with a low supporting value 50.2% for victoriviruses ( Figure 5). The ICTV species threshold for victoriviruses is set low: <60%. Thus, it is reasonable to classify DsVV1 in the same species as SsRV2, both of which were isolated from Diplodia fungi.

A Novel Splipalmivirus
Two contig sequences (contig numbers 218 and 310) showed resemblance with narna-like viruses [(+)ssRNA genome], particularly members of the newly proposed family "Splipalmiviridae" (Sutela et al., 2020;Sato et al., 2022) (Table 1). Splipalmiviruses have been coined for viruses with the split palm (italicized letters are used in the term "splipalmivirus") domains of RdRP that are encoded separately by different genomic segments (Sutela et al., 2020). The obtained viral sequences likely corresponded to the genomic segments of a splipalmivirus tentatively named Diplodia seriata splipalmivirus 1 (DsSpV1). The complete sequences of DsSpV1 RNA1 and RNA2 shared 12 conserved nucleotides (5′-UUUUGCUUGCGA-3′) and four nucleotides (5′-GUUU-3′) at the 5′ and 3′ termini, respectively (Supplementary Figure S6A), which suggest that both contigs represent the genome of a single virus. During data analyses, a third contig sequence (contig number 175) has been detected as a candidate of the third RNA segment of DsSpV1 (Table 1), however their terminal sequences were still not determined.
To estimate the phylogenetic relationships, the split RdRP sequences encoded by RNA1 and RNA2 of DsSpV1 and other splipalmiviruses, and the representative members of the family Narnaviridae were used to construct phylogenetic trees. The results showed that among the three proposed splipalmiviral genera , DsSpV1 was placed in a clade with high branch support values of 75.9% (for P1) and 99.9% (for P2) with duasplipalmiviruses ( Figure 6C). Currently, there is no established species demarcation criteria set by ICTV for splipalmiviruses. DsSpV1 is divergent enough from previously reported duasplipalmiviruses (~58.3% or~45.3% for the split RdRP aa sequence identity) to classify in a new species in the proposed genus "Duasplipalmivirus". DsSpV1 will be among rare splipalmivirids that are characterized biologically (see below); its horizontal transmission was confirmed.
Phylogenetically, DsBOV1 made clade with scleroliviruses (genus Sclerolivirus) and DsBOV2 made clade with penoliviruses (genus Penolivirus) within family Botourmiaviridae (Figure 8). Based on the ICTV species demarcation criteria (<90% RdRP aa sequence identity) , DsBOV1 should represent a novel species within the genus Sclerolivirus in the family Botourmiaviridae, while DsBOV2 should be classified as a strain belonging to the same species as NpOLV1 within the genus Penolivirus.
Phylogenetic analysis with RdRP sequences showed that DsAV1 is grouped together with other reported fungal RNA viruses (ambiguiviruses) and well-separated from other related viruses from plants (family Tombusviridae) and insects (unassigned viruses) ( Figure 8B). Classification of possible members of the recently proposed family "Ambugiviridae" has not been established. Based on the high degree of sequence similarity (~60.0% RdRP aa sequence identity), both DsAV1 and EnAAV1 should belong to a same novel species in the proposed family. Biological information is unavailable for all reported ambiguiviruses. In this regard, it is of importance that DsAV1 A B C FIGURE 6 | Genome architecture and phylogenetic position of Diplodia seriata splipalmivirus 1 (DsSpV1). (A) Schematic representation of the DsSpV1genome. The DsSpV1 genome has three segments encoding for split RdRP (P1 and P2) and a hypothetical protein (P3). While the RNA1 and RNA2 sequences were completely determined, the third segment, RNA3, has yet to be fully determined. (B, C) RdRP-based phylogenetic analysis of DsSpV1 and other splipalmiviruses (proposed family "Splipalmiviridae") and narna, and narna-like viruses. The best fit model (LG+G+I+F for P1 and VT+G+F for P2) was selected as a substitution model for evaluating the phylogenetic relationships. The tree was rooted at mid-point, and classical narnaviruses formed an out-group (B, Figure 7). could be eliminated from and re-introduced into the host fungus (see below).

Elimination of Some Viruses From D. seriata Strain L3
As shown above, D. seriata strain L3 was found to harbor at least eight mycoviruses. To check the possible effect of these viruses on the host fungal strain, we attempted to cure it of the viruses by a few methods. While single spore isolation is often used for this purpose, L3 turned out not to produce spores on PDA media. We then performed hyphal tipping and over 30 subcultures were tested for virus infection. One subculture termed L3ht1, showed a colony morphology different from L3 (Figure 1), and was predicted to have lost a few mycoviruses originally carried in the L3 strain. This could be seen from a dsRNA profile, showing the absence of expected dsRNAs of DsCV1 (3594 bp, 3265 bp, 3058 bp and 3202 bp) and DsPmV1(2439DsPmV1( bp, 2189DsPmV1( bp, 2014DsPmV1( bp, and 1133 in L3ht1 that were present in L3 (Figure 1). Note that some (+)ssRNA viruses, such as splipalmiviruses do not accumulate replicative dsRNA forms in infected fungal host cells . To confirm which mycoviruses were eliminated in L3ht1, RT-PCR was carried out using the nine primer sets for eight mycoviruses, including split RdRP encoded segments, as listed in Supplementary Table S1. Consequently, the L3ht1 was shown to harbor five mycoviruses, DsVV1, DsCV1, DsPV1, DsBOV1 and DsBOV2, and to have lost three other mycoviruses, i.e., DsPmV1, DsSpV1 and DsAV1 ( Figure 1). For multisegmented RNA viruses such as DsCV1, DsPV1, DsPmV1 and DsSpV1, consistent RT-PCR results were obtained in which different primer sets targeting different genomic segments of a virus of interest were used (Figure 1). Gel electrophoretic analyses of dsRNA and total nucleic acids suggested that whereas L3ht1 accumulated DsVV1 dsRNA more than L3; an opposite trend was observed for DsCV1 ( Figure 1). This interesting phenomenon warrants further investigation but suggests antagonistic and synergistic interactions between the eliminated viruses and the viruses whose accumulation was enhanced and suppressed in L3ht1, respectively.
Although we could not isolate virus-free subcultures, phenotypic comparison of the partially cured strain L3ht1 showed milder symptoms, increased colony diameter and regular colony margin, as compared to the original strain L3 (data not shown). Thus, this suggests that one or more of the three mycoviruses is involved in the symptom induction. However, it remains unknown how the respective eliminated mycoviruses contribute to the phenotypic alterations.

CONCLUSION
Thus far, only one mycovirus (Diplodia seriata betaendornavirus 1, DsEV1) has been reported from D. seriata (Nerva et al., 2019b). Here we report the characterization of co-infection of a single strain of D. seriata from Pakistan by eight mycoviruses. Of the eight identified mycoviruses, DsCV1, DsPV1, DsSpV1, DsBOV1 and DsAV1 appear to represent new mycovirus species based on the species demarcation criteria set by the ICTV or the high degree of their sequence diversity. The remaining three mycoviruses are putative new strains of previously recognized mycoviral species. The identified mycoviruses belong to the families Chrysoviridae, Polymycoviridae, Totiviridae, Partitiviridae and Botourmiaviridae, and the proposed families "Ambiguiviridae" and "Splipalmiviridae". As stated above, fungal virus hunting studies revealed new genome organizations and new virus lifestyles. Particularly, polymycoviruses and splipalmiviruses are considerably peculiar among the viruses discovered in this study. Polymycoviruses are phylogenetically closely related to hadakaviruses and relatively distantly related to caliciviruses with (+)ssRNA genomes (Kanhayuwa et al., 2015;Sato et al., 2020b;Khan et al., 2021a). However, polymycoviruses are classified as dsRNA viruses, because they are likely infectious in the form of dsRNA deproteinized and associated with PASrp as well (Kanhayuwa et al., 2015). It is noteworthy that DsPmV1 belongs to the family Polymycoviridae, whose infectivity as dsRNA needs to be verified. The closest relative of DsPmV1 is Botryosphaeria dothidea virus 1 (BdRV1) (Zhai et al., 2016), which has five genomic dsRNA segments. We failed to detect a fifth dsRNA of DsPmV1 by a similarity-based search of our NGS data.
"Splipalmivirus" is a recently proposed group of (+)ssRNA viruses with bifurcate RdRP domains. Phylogenetically, splipalmiviruses belong to the phylum Lenarviricota, which has been expanded rapidly, and most closely related to the established family Narnaviridae. A number of eukaryotic Lenarviricota members including splipalmiviruses, have been exponentially reported from various filamentous fungi, insects and plants, but most of them have yet to be characterized biologically (Sutela et al., 2020;Chiba et al., 2021a;Chiba et al., 2021b;Jia et al., 2021;Ruiz-Padilla et al., 2021;Sato et al., 2022). One of the best-studied fungal viruses is Saccharomyces 23S RNA narnavirus (ScNV23S), which has an undivided (+)ssRNA genome and exists in infected cells as a capsidless RNA/RdRP complex at a 1:1 molecular ratio (Solorzano et al., 2000;Esteban and Fujimura, 2003;Wickner et al., 2013). A variety of genome types are now detectable that include narnaviruses with multi-segmented or ambisense nature (Sutela et al., 2020;Kondo et al., 2022). The most peculiar narnavirus type is splipalmiviruses because no viruses with divided RdRP domains are found in other members of the kingdom Orthornavirae. Their RNA segment numbers vary from two to seven. Although three segments of DsSpV1 have been identified in this study, we cannot rule out that the virus has additional unidentified segments. There are several other open questions about splipalmiviruses. For example, where do splipalmiviruses replicate within host cells? Is association of the two splipalmivirus-encoded proteins with the separate RdRP domains necessary for RNA synthesis?
Interesting virus/virus interactions are expected in strain L3 coinfected with eight viruses, because in partially cured strain L3ht1 an enhanced accumulation of DsVV1 was observed (Figure 1). An antagonistic interaction appeared to exist between DsVV1 and one or more of the three removed viruses, i.e., DsPmV1, DsSpV1 and DsAV1. One of the best-studied antagonistic virus/virus interactions in fungi is that between Rosellinia necatrix victorivirus 1 (RnVV1, a victorivirus) and Cryphonectria hypovirus 1 Dp69 (CHV1-Dp69, a hypovirus mutant lacking an RNA silencing suppressor) or mycoreovirus 1 (MyRV1, a mycoreovirus) (Chiba and Suzuki, 2015). RnVV1 is eliminated or severely suppressed in replication by the RNA silencing status highly activated by coinfecting CHV1-Dp69 or MyRV1. A synergistic interaction was anticipated between DsCV1 and one or more of the three eliminated viruses, because in L3ht1 DsCV1accumulation was decreased ( Figure 1). This may be a reminiscent of the synergism between a hypovirus and a mycoreovirus (Sun et al., 2006). One or more of the three eliminated viruses are implicated in the abnormal colony morphology (Figure 1). As DsPmV1 and DsSpV1 belong to newly established or proposed virus families, for most members, biological information is unavailable. Therefore, this study provides a platform for further investigations into molecular and biologic virus/ virus interactions among the detected viruses.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

AUTHOR CONTRIBUTIONS
HAK: Sampling, investigation, writing -original draft. PT: Investigation and results interpretation. HK: Investigation, results interpretation, reviewing and editing original draft. MFB and NS: Conceptualization, supervision, results interpretation, writingreviewing and editing-original draft. All authors contributed to the article and approved the submitted version.

FUNDING
This investigation was partly supported by HEC NRPU grant (No-20-4109/R&D/HEC/14/847 awarded to MFB) and NUST students research, Grants-in-Aid for Scientific Research (A) and Grants-in-Aid for Scientific Research on Innovative Areas from the Japanese Ministry of Education, Culture, Sports, Science and Technology (KAKENHI 21H05035, 17H01463, 16H06436, 16H06429 and 16K21723 to NS and HK).