Natural Populations from the Phytophthora palustris Complex Show a High Diversity and Abundance of ssRNA and dsRNA Viruses

We explored the virome of the “Phytophthora palustris complex”, a group of aquatic specialists geographically limited to Southeast and East Asia, the native origin of many destructive invasive forest Phytophthora spp. Based on high-throughput sequencing (RNAseq) of 112 isolates of “P. palustris” collected from rivers, mangroves, and ponds, and natural forests in subtropical and tropical areas in Indonesia, Taiwan, and Japan, 52 putative viruses were identified, which, to varying degrees, were phylogenetically related to the families Botybirnaviridae, Narnaviridae, Tombusviridae, and Totiviridae, and the order Bunyavirales. The prevalence of all viruses in their hosts was investigated and confirmed by RT-PCR. The rich virus composition, high abundance, and distribution discovered in our study indicate that viruses are naturally infecting taxa from the “P. palustris complex” in their natural niche, and that they are predominant members of the host cellular environment. Certain Indonesian localities are the viruses’ hotspots and particular “P. palustris” isolates show complex multiviral infections. This study defines the first bi-segmented bunya-like virus together with the first tombus-like and botybirna-like viruses in the genus Phytophthora and provides insights into the spread and evolution of RNA viruses in the natural populations of an oomycete species.


Introduction
Global trade has broken down the natural distribution ranges of species by enabling the long-distance movement of living organisms (including plant pests and pathogens) around the world and their establishment in new territories [1][2][3][4]. Alien invasive tree pathogens are recognized world-wide as a rising hazard to biodiversity and ecosystem functioning [2,5,6]. Therefore, the development of successful management strategies requires a profound understanding of the biology, epidemiology, and pathways of spread of the target organism [7]. Although challenging, identifying the original geographical and ecological niche of a given pathogen can not only help to predict the places where it may survive, but also identify natural antagonists, such as hyperparasites, which should be more frequent in the native range of their natural hosts than elsewhere [8,9].
A considerable proportion of known destructive alien forest and crop pathogens are oomycetes belonging to the genus Phytophthora, a group of eukaryotes phylogenetically related to brown algae, diatoms, and other Stramenopiles [10], which share much of their ecological niches, lifestyle, basic structural features, and virulence strategies with fungi [11]. Prominent examples of devastating emerging diseases include "chestnut ink disease" (caused by Phytophthora cinnamomi and P. × cambivora), "sudden oak death" in KAPA Library Quantification Kit for Illumina platform was used for absolute, qPCR-based quantification of the Illumina libraries flanked by the P5 and P7 flow cell oligo sequences. Libraries underwent paired-end (PE) (2 × 150 nt) sequencing on a NovaSeq6000 (DS-150) (Illumina, San Diego, CA, USA) using a NovaSeq S4 v1.5 reagent kit. An "in-lane" PhiX control spike was included in each lane of the flow cell.
Consequent analyses, including the determination of the final virus sequences, primer design, detection of open reading frames (ORF), protein translation, pairwise (pw) sequence comparison (PASC), and read number calculation, were performed using the platform Geneious Prime ® 2021.0.4. and 2022.0.1. For the calculation of the coverage depth, we used the following formula: (Total reads mapped to the final identified virus * average read length)/virus genome or contig length).

Genetic Variability Analyses and Conserved Domains
Pairwise identities of the nucleotide and amino acid sequences were obtained after aligning the viral nucleotide and amino acid sequences by MAFFT V1. 4.0 and calculated in Geneious Prime ® 2021.0.4. In order to search for conserved domains within the putative viral proteins, the NCBI CDD-search tool was used (https://www.ncbi.nlm.nih.gov/ Structure/cdd/wrpsb.cgi, accessed on 28 August 2022).

Phylogenetic Trees
Maximum likelihood (ML) phylogenetic trees were constructed using a rapid bootstrapping algorithm [50] in RAxML-HPC v.8 on XSEDE, conducted in CIPRES Science Gateway [51]. Tree search was enabled under the GAMMA model to avoid thorough optimization of the best scoring ML tree at the end of the run. The Jones-Taylor-Thornton (JTT) model was chosen as a substitution model for proteins. Bootstrapping was configured with the recommended parameters given by CIPRES Science Gateway. The resulting data were visualized using the software FIGTREE version 1.4.4.

Rapid Amplification of cDNA Ends (RACE) and Confirmation of Viruses' Occurrence by Direct Reverse-Transcription Polymerase Chain Reaction (RT-PCR)
To confirm the length of the RNA 1 of botybirna-like virus 2 and RNA 1 and 2 of bunyalike virus 11 (see Results), we used the SMARTer RACE 5 /3 KIT (TAKARABIO USA, Inc., Mountain View, CA, USA), as described in [46], and following the producer's instructions. The occurrence of each identified virus was confirmed by direct reverse-transcription polymerase chain reaction (RT-PCR) with specific primers using total RNA as a template. The screening did not differentiate among variants, and the variant-level prevalence was not studied. A High-Capacity cDNA Reverse Transcription Kit (Applied Biosciences, Park Ave, NY, USA) was used for the cDNA synthesis. PCRs were performed with a Hot Start Taq 2× Master Mix (New England BioLabs, Ipswich, MA, USA) including 25 µL Master Mix, 1 µL of each primer (10 mM), and 4 µL of cDNA in a total volume of 50 µL. RT-PCR products were visualized using gel electrophoresis (120 V; 60 min). Analyzed fragments were separated on 1.5% agarose gel prepared with a TBE 1X buffer (Merck KGaA, Gernsheim, Germany) and stained by Ethidium bromide (SIGMA-Aldrich, Steinheim, Germany). PCR products showing the amplicons of expected length were purified and sequenced by GATC BioTech (Eurofins; Konstanz, Germany) by both directions with the primers used for PCR amplification. All the primers used for the partial amplification of the RNA-dependent RNA polymerase (RdRP) of each virus (Table S2) were designed by Primer 3 2.3.7 under Geneious Prime ® 2020.0.4.

Graphs and Visualizations
Maps and graphs were created based on physically measured GPS coordinates using the R version 4.0.5 [52] programming environment using the special libraries ggplot2, sf, mapplots, ggspatial, scatterpie, and ggrepel.

Virus Identification
A total of 10 RNA libraries containing 10-15 RNAs pooled from "P. palustris" isolates from Japan, Taiwan, and Indonesia (Kalimantan and Sumatra) were constructed and sequenced (Tables 1 and S1). After trimming and quality checking,~5200 million reads were obtained (Table S3). Further removal of reads mapped to the host genome yielded a total of 259,086,400 reads. The number of host-mapped reads was very variable among the different pools (Table S4) but generally very high, always higher than 84%. The highest percentage was obtained in pools JP2, SU2, and SU4 (>98%) and the lowest in pool KA1 (84.41%). A total of 730,972 possible viral contigs (length > 0-25,000 bp) were obtained by de novo assembly of viral contigs (Table S5). After being aligned using BLAST+ to all the databases mentioned in the Materials and Methods section, virus contigs were finally identified to represent 52 viruses related to five virus families and orders: Totiviridae and Botybirnaviridae with dsRNA genomes; Bunyavirales with (−)ssRNA genomes; and Narnaviridae and Tombusviridae as (+)ssRNA genomes (Table 1). No DNA viruses were identified applying our pipeline.

(−)ssRNA Viruses
Based on the BLASTX comparison, a total of 18 viral contigs with affinities with the L (large) segment of putative members of the order Bunyavirales were assembled and characterized in Sumatran and Kalimantan pools (SU1-4 and KA1-3) ( Table 1). All of them enclosed a unique large open reading frame (ORF) encoding the RNA-dependent RNA polymerase (RdRP) (Figure 1a). In addition, one contig showed a high identity percentage with nucleocapsid proteins (NC) of the family Phenuiviridae. The 17 putative viruses were designated as Phytophthora palustris bunya-like viruses (PpaBLV) ( Table 1). The PASC of the 17 bunya-like RdRP nucleotide (nt) and amino (aa) sequences showed an overall pw identity of 34 and 39.4%, respectively. When all the RdRP nt and aa sequences were compared, several contigs were seen to have a pw identity >90% (Tables S8 and S9); therefore, those sequences were designated as variants of the same putative virus. Thus, both PpaBLV 9 and 14 had three variants. In ICTV, there are no primary classification and delimitation criteria for genus and species in the order Bunyavirales, and PASC and phylogenetic analyses seem to be the main point of reference to name new bunyaviruses.  A comprehensive reconstruction of the phylogenetic relationships of "P. palustris" bunya-like viruses' RdRP with other bunyaviruses deposited in the GenBank is shown in Figure 2a. The majority of "P. palustris" bunya-like viruses cluster together and seem to be closer to Phytophthora condilina negative-stranded RNA virus 4 (PcoNSRV4), Halophytophthora RNA virus 1 (HRV1), both viruses described in marine oomycetes collected Conserved domains (CDD) of the Bunya_RdRP superfamily cl20265 were found in bunya-like virus contigs 1, 2, 3, 4, 5, 6, 8, 9 (variant 1), 10, 11, and 13 ( Table 1). The amino acid alignment of the CDD regions ( Figure S1a) showed high similarities of the premotif A and motifs A (DxxxWx), B (XGxxNxxSS), C (SDD), D (KK), and E (ExxSx) with the rest of the bunyaviruses included in the alignment. Premotif A with the three basic residues inside (K, R, and R/K) and, downstream, the glutamic acid (E), were also identified. The conserved aa triplet TPD (threonine), typical of bunyaviruses, was also detected [53]. The aa sequence of the NC of the Phytophthora palustris bunya-like virus 11 encloses a conserved motif of pfam05733, the only member of the superfamily cl05345 with a significant e-value (2.08 × 10 −7 ) (Table 1, Figure 1a). This family consists of several Tenuivirus and Phlebovirus nucleocapsid proteins ( Figure S1, Table 1). In addition to the conserved domain of Bunya_RdRP, an L-protein N-terminus (and endonuclease domain) present in the N-terminus of many bunyavirus L proteins was found at nts 7912-8166.
On the basis of the RACE analysis, we determined the possible full viral genome sequence of PpaBLV11, which comprised RNA1 (6448 nt) and RNA2 (864 nt) (Figure 1a). The RNA2 sequence is shorter than the NGS contig (ca. 1 kb) because only 2 relatively short sequences of 3 were obtained by the RACE analysis. Moreover, the coverage of the ends of the NGS-RNA2 contig was very low and had very low quality. The 5 and 3 termini of RNA1 and RNA2 had 9 and 8 conserved nts, respectively ( Figure 1b). A total of 9 sequential nts at the 5 and 3 termini of RNA1 and 8 nts in RNA 2 were complementary to each other ( Figure 1c). The comparison of the 5 and 3 terminal nucleotides of PpaBLV11 to those of other bunyaviruses showed a high identity. Higher variation appears to occur in the 3 terminus of the RNA2, having a C, which is lacking in the rest of the 3 termini of the compared bunyaviruses; similarly, no complementary G is found in the 5 terminus of PpaBLV1 RNA1.
ORF1, the unique large ORF found in RNA1, encodes a putative protein of 2084 aa (codons) and a molecular weight of 241.396 kDa (p241), and 5 and 3 untranslated regions (UTRs) comprising 72 and 121 nt, respectively ( Figure 1a). ORF2, the only ORF found in RNA2, encodes a putative protein of 230 aa and a molecular weight of 26.459 kDa (p26).
A comprehensive reconstruction of the phylogenetic relationships of "P. palustris" bunya-like viruses' RdRP with other bunyaviruses deposited in the GenBank is shown in Figure 2a. The majority of "P. palustris" bunya-like viruses cluster together and seem to be closer to Phytophthora condilina negative-stranded RNA virus 4 (PcoNSRV4), Halophytophthora RNA virus 1 (HRV1), both viruses described in marine oomycetes collected in estuarine ecosystems in southern Portugal [29,42], and to Phytophthora cactorum virus 1 (PcBV2), described in several isolates of P. cactorum causing crown rot in strawberries in Finland [33]. PpaBLV5 appears jointly with Phytophthora cactorum bunya virus 1 [33] in a closely related cluster, which just hosts bunya-like viruses found in oomycetes. PpaBLV6 is the most distinct taxon, clustering with Phytophthora condilina negative-stranded RNA virus 6 (PcoNSRV6) [29]. These results indicate a close evolutionary relationship among bunyalike viruses, described in different species of oomycetes collected from different environments.
The NC of PpaBLV11 appears to be phylogenetically related to NCs of an insect virus (genus Phasivirus) and a plant virus, and not related to those of fungal phenuiviruses (Figure 2b).

(+)ssRNA Viruses
Tombus-like virus. One viral contig of 4180 bp found in the Sumatran pool SU2 appeared to have phylogenetic affinities with unclassified members of the family Tombusviridae. Although the best hit subjects in the BLASTX comparison were sequences obtained in a metatranscriptomic study from soil dominated by stands of Avena species (Table 1), other related results were Leuven tombus-like virus 3, described from larvae of predatory wasps in Belgium (e-value 8 × 10 −49 ; identity ~32%), and Serdyukov virus (e-value 5 × 10 −45 ; identity ~36%), reported in the ticks of Antarctic penguins [54].
When using the standard translation code, the contig encodes two positive-sense ORFs (Figure 3a), and the ORF1 or 5′ end-proximal ORF encodes an unknown protein with no conserved domains detected. The ORF2 or 3′ end-proximal ORF encodes the RdRP. The conserved domain of the RT-like superfamily (cl2808) is located at nt positions 1797-3191 (e-value 2.11 × 10 −57 ), particularly, the superfamily member pfam00998 (Viral RdRP_3) found in Hepatitis C virus and other plant viruses ( Figure S2). A third antisense ORF (nts 635-63) is detected proximal to 5′ end and just separated by one proline (Pro) from the next ORF1. The resultant protein (191 aa) does not have any similar records when

(+)ssRNA Viruses
Tombus-like virus. One viral contig of 4180 bp found in the Sumatran pool SU2 appeared to have phylogenetic affinities with unclassified members of the family Tombusviridae. Although the best hit subjects in the BLASTX comparison were sequences obtained in a metatranscriptomic study from soil dominated by stands of Avena species (Table 1), other related results were Leuven tombus-like virus 3, described from larvae of predatory wasps in Belgium (e-value 8 × 10 −49 ; identity ~32%), and Serdyukov virus (e-value 5 × 10 −45 ; identity ~36%), reported in the ticks of Antarctic penguins [54].
When using the standard translation code, the contig encodes two positive-sense ORFs (Figure 3a), and the ORF1 or 5′ end-proximal ORF encodes an unknown protein with no conserved domains detected. The ORF2 or 3′ end-proximal ORF encodes the RdRP. The conserved domain of the RT-like superfamily (cl2808) is located at nt positions 1797-3191 (e-value 2.11 × 10 −57 ), particularly, the superfamily member pfam00998 (Viral RdRP_3) found in Hepatitis C virus and other plant viruses ( Figure S2). A third antisense ORF (nts 635-63) is detected proximal to 5′ end and just separated by one proline (Pro) from the next ORF1. The resultant protein (191 aa) does not have any similar records when

(+)ssRNA Viruses
Tombus-like virus. One viral contig of 4180 bp found in the Sumatran pool SU2 appeared to have phylogenetic affinities with unclassified members of the family Tombusviridae. Although the best hit subjects in the BLASTX comparison were sequences obtained in a metatranscriptomic study from soil dominated by stands of Avena species (Table 1), other related results were Leuven tombus-like virus 3, described from larvae of predatory wasps in Belgium (e-value 8 × 10 −49 ; identity ~32%), and Serdyukov virus (e-value 5 × 10 −45 ; identity ~36%), reported in the ticks of Antarctic penguins [54].
When using the standard translation code, the contig encodes two positive-sense ORFs (Figure 3a), and the ORF1 or 5′ end-proximal ORF encodes an unknown protein with no conserved domains detected. The ORF2 or 3′ end-proximal ORF encodes the RdRP. The conserved domain of the RT-like superfamily (cl2808) is located at nt positions 1797-3191 (e-value 2.11 × 10 −57 ), particularly, the superfamily member pfam00998 (Viral RdRP_3) found in Hepatitis C virus and other plant viruses ( Figure S2

(+)ssRNA Viruses
Tombus-like virus. One viral contig of 4180 bp found in the Sumatran pool SU2 appeared to have phylogenetic affinities with unclassified members of the family Tombusviridae. Although the best hit subjects in the BLASTX comparison were sequences obtained in a metatranscriptomic study from soil dominated by stands of Avena species (Table 1), other related results were Leuven tombus-like virus 3, described from larvae of predatory wasps in Belgium (e-value 8 × 10 −49 ; identity ~32%), and Serdyukov virus (e-value 5 × 10 −45 ; identity ~36%), reported in the ticks of Antarctic penguins [54].
When using the standard translation code, the contig encodes two positive-sense ORFs (Figure 3a), and the ORF1 or 5′ end-proximal ORF encodes an unknown protein with no conserved domains detected. The ORF2 or 3′ end-proximal ORF encodes the RdRP. The conserved domain of the RT-like superfamily (cl2808) is located at nt positions 1797-3191 (e-value 2.11 × 10 −57 ), particularly, the superfamily member pfam00998 (Viral RdRP_3) found in Hepatitis C virus and other plant viruses ( Figure S2

(+)ssRNA Viruses
Tombus-like virus. One viral contig of 4180 bp found in the Sumatran pool SU2 appeared to have phylogenetic affinities with unclassified members of the family Tombusviridae. Although the best hit subjects in the BLASTX comparison were sequences obtained in a metatranscriptomic study from soil dominated by stands of Avena species (Table 1), other related results were Leuven tombus-like virus 3, described from larvae of predatory wasps in Belgium (e-value 8 × 10 −49 ; identity ~32%), and Serdyukov virus (e-value 5 × 10 −45 ; identity ~36%), reported in the ticks of Antarctic penguins [54].
When using the standard translation code, the contig encodes two positive-sense ORFs (Figure 3a), and the ORF1 or 5′ end-proximal ORF encodes an unknown protein with no conserved domains detected. The ORF2 or 3′ end-proximal ORF encodes the RdRP. The conserved domain of the RT-like superfamily (cl2808) is located at nt positions 1797-3191 (e-value 2.11 × 10 −57 ), particularly, the superfamily member pfam00998 (Viral RdRP_3) found in Hepatitis C virus and other plant viruses ( Figure S2

(+)ssRNA Viruses
Tombus-like virus. One viral contig of 4180 bp found in the Sumatran pool SU2 appeared to have phylogenetic affinities with unclassified members of the family Tombusviridae. Although the best hit subjects in the BLASTX comparison were sequences obtained in a metatranscriptomic study from soil dominated by stands of Avena species (Table 1), other related results were Leuven tombus-like virus 3, described from larvae of predatory wasps in Belgium (e-value 8 × 10 −49 ; identity ~32%), and Serdyukov virus (e-value 5 × 10 −45 ; identity ~36%), reported in the ticks of Antarctic penguins [54].
When using the standard translation code, the contig encodes two positive-sense ORFs (Figure 3a), and the ORF1 or 5′ end-proximal ORF encodes an unknown protein with no conserved domains detected. The ORF2 or 3′ end-proximal ORF encodes the RdRP. The conserved domain of the RT-like superfamily (cl2808) is located at nt positions 1797-3191 (e-value 2.11 × 10 −57 ), particularly, the superfamily member pfam00998 (Viral RdRP_3) found in Hepatitis C virus and other plant viruses ( Figure S2). A third antisense ORF (nts 635-63) is detected proximal to 5′ end and just separated by one proline (Pro) from the next ORF1. The resultant protein (191 aa) does not have any similar records when

(+)ssRNA Viruses
Tombus-like virus. One viral contig of 4180 bp found in the Sumatran pool SU2 appeared to have phylogenetic affinities with unclassified members of the family Tombusviridae. Although the best hit subjects in the BLASTX comparison were sequences obtained in a metatranscriptomic study from soil dominated by stands of Avena species (Table 1), other related results were Leuven tombus-like virus 3, described from larvae of predatory wasps in Belgium (e-value 8 × 10 −49 ; identity ~32%), and Serdyukov virus (e-value 5 × 10 −45 ; identity ~36%), reported in the ticks of Antarctic penguins [54].
When using the standard translation code, the contig encodes two positive-sense ORFs (Figure 3a), and the ORF1 or 5′ end-proximal ORF encodes an unknown protein with no conserved domains detected. The ORF2 or 3′ end-proximal ORF encodes the RdRP. The conserved domain of the RT-like superfamily (cl2808) is located at nt positions 1797-3191 (e-value 2.11 × 10 −57 ), particularly, the superfamily member pfam00998 (Viral RdRP_3) found in Hepatitis C virus and other plant viruses ( Figure S2

(+)ssRNA Viruses
Tombus-like virus. One viral contig of 4180 bp found in the Sumatran pool SU2 appeared to have phylogenetic affinities with unclassified members of the family Tombusviridae. Although the best hit subjects in the BLASTX comparison were sequences obtained in a metatranscriptomic study from soil dominated by stands of Avena species (Table 1), other related results were Leuven tombus-like virus 3, described from larvae of predatory wasps in Belgium (e-value 8 × 10 −49 ; identity ~32%), and Serdyukov virus (e-value 5 × 10 −45 ; identity ~36%), reported in the ticks of Antarctic penguins [54].
When using the standard translation code, the contig encodes two positive-sense ORFs (Figure 3a), and the ORF1 or 5′ end-proximal ORF encodes an unknown protein with no conserved domains detected. The ORF2 or 3′ end-proximal ORF encodes the RdRP. The conserved domain of the RT-like superfamily (cl2808) is located at nt positions 1797-3191 (e-value 2.11 × 10 −57 ), particularly, the superfamily member pfam00998 (Viral RdRP_3) found in Hepatitis C virus and other plant viruses ( Figure S2). A third antisense ORF (nts 635-63) is detected proximal to 5′ end and just separated by one proline (Pro) from the next ORF1. The resultant protein (191 aa) does not have any similar records when

(+)ssRNA Viruses
Tombus-like virus. One viral contig of 4180 bp found in the Sumatran pool SU2 appeared to have phylogenetic affinities with unclassified members of the family Tombusviridae. Although the best hit subjects in the BLASTX comparison were sequences obtained in a metatranscriptomic study from soil dominated by stands of Avena species (Table 1), other related results were Leuven tombus-like virus 3, described from larvae of predatory wasps in Belgium (e-value 8 × 10 −49 ; identity ~32%), and Serdyukov virus (e-value 5 × 10 −45 ; identity ~36%), reported in the ticks of Antarctic penguins [54].
When using the standard translation code, the contig encodes two positive-sense ORFs (Figure 3a), and the ORF1 or 5′ end-proximal ORF encodes an unknown protein with no conserved domains detected. The ORF2 or 3′ end-proximal ORF encodes the RdRP. The conserved domain of the RT-like superfamily (cl2808) is located at nt positions 1797-3191 (e-value 2.11 × 10 −57 ), particularly, the superfamily member pfam00998 (Viral RdRP_3) found in Hepatitis C virus and other plant viruses ( Figure S2

(+)ssRNA Viruses
Tombus-like virus. One viral contig of 4180 bp found in the Sumatran pool SU2 appeared to have phylogenetic affinities with unclassified members of the family Tombusviridae. Although the best hit subjects in the BLASTX comparison were sequences obtained in a metatranscriptomic study from soil dominated by stands of Avena species (Table 1), other related results were Leuven tombus-like virus 3, described from larvae of predatory wasps in Belgium (e-value 8 × 10 −49 ; identity ~32%), and Serdyukov virus (e-value 5 × 10 −45 ; identity ~36%), reported in the ticks of Antarctic penguins [54].
When using the standard translation code, the contig encodes two positive-sense ORFs (Figure 3a), and the ORF1 or 5′ end-proximal ORF encodes an unknown protein with no conserved domains detected. The ORF2 or 3′ end-proximal ORF encodes the RdRP. The conserved domain of the RT-like superfamily (cl2808) is located at nt positions 1797-3191 (e-value 2.11 × 10 −57 ), particularly, the superfamily member pfam00998 (Viral RdRP_3) found in Hepatitis C virus and other plant viruses ( Figure S2). A third antisense ORF (nts 635-63) is detected proximal to 5′ end and just separated by one proline (Pro) from the next ORF1. The resultant protein (191 aa) does not have any similar records when Excavata. Scale bars represent expected changes per site per branch.

(+)ssRNA Viruses
Tombus-like virus. One viral contig of 4180 bp found in the Sumatran pool SU2 appeared to have phylogenetic affinities with unclassified members of the family Tombusviridae. Although the best hit subjects in the BLASTX comparison were sequences obtained in a metatranscriptomic study from soil dominated by stands of Avena species (Table 1), other related results were Leuven tombus-like virus 3, described from larvae of predatory wasps in Belgium (e-value 8 × 10 −49 ; identity~32%), and Serdyukov virus (e-value 5 × 10 −45 ; identity~36%), reported in the ticks of Antarctic penguins [54].
When using the standard translation code, the contig encodes two positive-sense ORFs (Figure 3a), and the ORF1 or 5 end-proximal ORF encodes an unknown protein with no conserved domains detected. The ORF2 or 3 end-proximal ORF encodes the RdRP. The conserved domain of the RT-like superfamily (cl2808) is located at nt positions 1797-3191 (e-value 2.11 × 10 −57 ), particularly, the superfamily member pfam00998 (Viral RdRP_3) found in Hepatitis C virus and other plant viruses ( Figure S2). A third antisense ORF (nts 635-63) is detected proximal to 5 end and just separated by one proline (Pro) from the next ORF1. The resultant protein (191 aa) does not have any similar records when blasted to GenBank. This contig was named Phytophthora palustris tombus-like virus 1 (PpaTbLV1). The phylogenetic tree (Figure 3b) shows that PpaTbLV1 is more related to tombusviruses infecting insects rather than those infecting the obligate biotrophic oomycetes Sclerophthora macrospora [55] and Plasmopara halstedii [42]. blasted to GenBank. This contig was named Phytophthora palustris tombus-like virus 1 (PpaTbLV1). The phylogenetic tree (Figure 3b) shows that PpaTbLV1 is more related to tombusviruses infecting insects rather than those infecting the obligate biotrophic oomycetes Sclerophthora macrospora [55] and Plasmopara halstedii [42]. Narna-like viruses. Virus contigs resembling members of the family Narnaviridae were abundant in all pools. After removing redundant contigs and selecting the longest ones, 10 viruses were determined. Eight of them had lengths longer than 2.5 kb and

(+)ssRNA Viruses
Tombus-like virus. One viral contig of 4180 bp found in the Sumatran pool SU2 appeared to have phylogenetic affinities with unclassified members of the family Tombusviridae. Although the best hit subjects in the BLASTX comparison were sequences obtained in a metatranscriptomic study from soil dominated by stands of Avena species (Table 1), other related results were Leuven tombus-like virus 3, described from larvae of predatory wasps in Belgium (e-value 8 × 10 −49 ; identity ~32%), and Serdyukov virus (e-value 5 × 10 −45 ; identity ~36%), reported in the ticks of Antarctic penguins [54].
When using the standard translation code, the contig encodes two positive-sense ORFs (Figure 3a), and the ORF1 or 5′ end-proximal ORF encodes an unknown protein with no conserved domains detected. The ORF2 or 3′ end-proximal ORF encodes the RdRP. The conserved domain of the RT-like superfamily (cl2808) is located at nt positions (+)ssRNA Viruses ombus-like virus. One viral contig of 4180 bp found in the Sumatran pool SU2 apto have phylogenetic affinities with unclassified members of the family Tombusvirilthough the best hit subjects in the BLASTX comparison were sequences obtained in atranscriptomic study from soil dominated by stands of Avena species (Table 1), related results were Leuven tombus-like virus 3, described from larvae of predatory in Belgium (e-value 8 × 10 −49 ; identity ~32%), and Serdyukov virus (e-value 5 × 10 −45 ; ty ~36%), reported in the ticks of Antarctic penguins [54]. hen using the standard translation code, the contig encodes two positive-sense (Figure 3a), and the ORF1 or 5′ end-proximal ORF encodes an unknown protein o conserved domains detected. The ORF2 or 3′ end-proximal ORF encodes the . The conserved domain of the RT-like superfamily (cl2808) is located at nt positions  (Table 1), elated results were Leuven tombus-like virus 3, described from larvae of predatory in Belgium (e-value 8 × 10 −49 ; identity ~32%), and Serdyukov virus (e-value 5 × 10 −45 ; y ~36%), reported in the ticks of Antarctic penguins [54]. hen using the standard translation code, the contig encodes two positive-sense Figure 3a), and the ORF1 or 5′ end-proximal ORF encodes an unknown protein o conserved domains detected. The ORF2 or 3′ end-proximal ORF encodes the The conserved domain of the RT-like superfamily (cl2808) is located at nt positions −57 Arthropoda, blasted to GenBank. This contig was named Phytophthora palustris tombus-like virus 1 (PpaTbLV1). The phylogenetic tree (Figure 3b) shows that PpaTbLV1 is more related to tombusviruses infecting insects rather than those infecting the obligate biotrophic oomycetes Sclerophthora macrospora [55] and Plasmopara halstedii [42].

(+)ssRNA Viruses
Tombus-like virus. One viral contig of 4180 bp found in the Sumatran pool SU2 appeared to have phylogenetic affinities with unclassified members of the family Tombusviridae. Although the best hit subjects in the BLASTX comparison were sequences obtained in a metatranscriptomic study from soil dominated by stands of Avena species (Table 1), other related results were Leuven tombus-like virus 3, described from larvae of predatory wasps in Belgium (e-value 8 × 10 −49 ; identity ~32%), and Serdyukov virus (e-value 5 × 10 −45 ; identity ~36%), reported in the ticks of Antarctic penguins [54].
When using the standard translation code, the contig encodes two positive-sense ORFs (Figure 3a), and the ORF1 or 5′ end-proximal ORF encodes an unknown protein with no conserved domains detected. The ORF2 or 3′ end-proximal ORF encodes the RdRP. The conserved domain of the RT-like superfamily (cl2808) is located at nt positions Narna-like viruses. Virus contigs resembling members of the family Narnaviridae were abundant in all pools. After removing redundant contigs and selecting the longest ones, 10 viruses were determined. Eight of them had lengths longer than 2.5 kb and enclosed a complete single large ORF, likely representing nearly full-length virus genomes (Table 1; Figure 4a). enclosed a complete single large ORF, likely representing nearly full-length virus genomes (Table 1; Figure 4a). Two contigs contained incomplete ORFs and represented partial genomes of two potential virus isolates. No narnavirus conserved domains were detected in any of the

(+)ssRNA Viruses
Tombus-like virus. One viral contig of 4 peared to have phylogenetic affinities with un dae. Although the best hit subjects in the BLAS a metatranscriptomic study from soil domin other related results were Leuven tombus-like wasps in Belgium (e-value 8 × 10 −49 ; identity ~3 identity ~36%), reported in the ticks of Antarc When using the standard translation co Viruses ike virus. One viral contig of 4180 bp found in the Sumatran pool SU2 apphylogenetic affinities with unclassified members of the family Tombusvirithe best hit subjects in the BLASTX comparison were sequences obtained in ptomic study from soil dominated by stands of Avena species (Table 1), sults were Leuven tombus-like virus 3, described from larvae of predatory um (e-value 8 × 10 −49 ; identity ~32%), and Serdyukov virus (e-value 5 × 10 −45 ; , reported in the ticks of Antarctic penguins [54]. Viruses ke virus. One viral contig of 4180 bp found in the Sumatran pool SU2 apphylogenetic affinities with unclassified members of the family Tombusvirihe best hit subjects in the BLASTX comparison were sequences obtained in tomic study from soil dominated by stands of Avena species (Table 1), sults were Leuven tombus-like virus 3, described from larvae of predatory m (e-value 8 × 10 −49 ; identity ~32%), and Serdyukov virus (e-value 5 × 10 −45 ; reported in the ticks of Antarctic penguins [54]. blasted to GenBank. This contig was named Phytophthora palustris tombus-like virus 1 (PpaTbLV1). The phylogenetic tree (Figure 3b) shows that PpaTbLV1 is more related to tombusviruses infecting insects rather than those infecting the obligate biotrophic oomycetes Sclerophthora macrospora [55] and Plasmopara halstedii [42].

. (+)ssRNA Viruses
Tombus-like virus. One viral contig of 4180 bp found in the Sumatran pool SU2 aped to have phylogenetic affinities with unclassified members of the family Tombusvirilthough the best hit subjects in the BLASTX comparison were sequences obtained in tatranscriptomic study from soil dominated by stands of Avena species (Table 1), r related results were Leuven tombus-like virus 3, described from larvae of predatory s in Belgium (e-value 8 × 10 −49 ; identity ~32%), and Serdyukov virus (e-value 5 × 10 −45 ; tity ~36%), reported in the ticks of Antarctic penguins [54].

(+)ssRNA Viruses
Tombus-like virus. One viral contig of 4180 bp found in the Sumatran pool SU peared to have phylogenetic affinities with unclassified members of the family Tombu dae. Although the best hit subjects in the BLASTX comparison were sequences obtain a metatranscriptomic study from soil dominated by stands of Avena species (Tab other related results were Leuven tombus-like virus 3, described from larvae of preda wasps in Belgium (e-value 8 × 10 −49 ; identity ~32%), and Serdyukov virus (e-value 5 × identity ~36%), reported in the ticks of Antarctic penguins [54].
When using the standard translation code, the contig encodes two positive-s Ochrophyta (Heterokonta),  Figure 4a).  Two contigs contained incomplete ORFs and represented partial genomes of two potential virus isolates. No narnavirus conserved domains were detected in any of the putative virus contigs in the Conserved Domain Database (CDD) but the BLASTX search confirmed their similarity with unclassified narna-like viruses found in the fungal pathogen Eryshiphe necator and in invertebrates, and with Phytophthora infestans RNA virus 4 ( Table 1). The alignment of the conserved motifs of RdRP of the narna-like viruses is shown in Figure S4. The overall PASC of the eight narna-like contig nt and RdRP aa sequences were 34.2 and 27.1%, respectively. Based on PASC in the nt and aa sequences (Tables S10 and S11) and the species demarcation criterion proposed in the ICTV for the family Narnaviridae, eight possible viruses designated as Phytophthora palustris narna-like viruses 1-8 (PpaNLV1-8) were represented in these 10 contigs. PpaNLV3 was represented by three variants as the pw identity of their aa sequences was higher than 90%.
The phylogenetic analyses, based on the aa sequences of their RdRP (Figure 4b), showed that "P. palustris" narna-like viruses are grouped in two different clusters. PpaNLV2, 3 (including all 3 variants), 4, 5, and 8 appeared in a detached group, maybe representing a different taxon within the "alphanarnavirus" clade, and closer to Saccharomycces 23S RNA narnavirus and Saccharomycces 20S RNA narnavirus. Within the clade, "betanarnaviruses" PpaNLV1 and 6 appear together with narnaviruses described in Phytophthora spp. and PpaNLV7 groups with two unclassified viruses described in an insect and in symptomatic grapevine tissue associated with the downy mildew Plasmopara viticola.

dsRNA Viruses
Botybirna-like viruses. Four contigs >6 kb in length enclosing one large ORF resembling the RdRP segment of unclassified botybirnaviruses present in fungi and oomycetes ( Figure 5a; Table 1) were found in two pools from Sumatra (SU1, SU3) and one from Kalimantan (KA3). No second segment was detected in our data. These viral contigs were named as Phytophthora palustris botybirna-like viruses (PpaBbLV). The PASC of the nt and aa sequences of the four contigs showed high identity between them (Tables S6 and S7), with particularly high identity levels (>92%) between the Sumatran contigs. Because the ICTV does not specify virus species demarcation criteria for this family, we classified the four sequences as belonging to two viruses: three contigs represent three Sumatran variants belonging to PpaBbL1 while one represents a Kalimantan botybirna-like virus. CDD belonging to the RT-like superfamily (cl2808) and, particularly, the superfamily member pfam02123 (RdRP_4, which includes RdRPs from Luteovirus, Totivirus and Rotavirus) were also found (Figure 5a and Figure S3; Table 1). The reconstruction of the phylogenetic relationships of PpaBbLV1 and 2 with other dsRNA viruses (Figure 5b) shows that PpaB-bLV1 and 2 form a novel cluster separated from other botybirnaviruses found in fungi and oomycetes, suggesting an evolutionary diversification among the mycoviruses of this family.  Table 1). The alignment of the conserved motifs of RdRP of the narna-like viruses is shown in Figure S4. The overall PASC of the eight narna-like contig nt and RdRP aa sequences were 34.2 and 27.1%, respectively. Based on PASC in the nt and aa sequences (Tables S10 and S11) and the species demarcation criterion proposed in the ICTV for the family Narnaviridae, eight possible viruses designated as Phytophthora palustris narnalike viruses 1-8 (PpaNLV1-8) were represented in these 10 contigs. PpaNLV3 was represented by three variants as the pw identity of their aa sequences was higher than 90%. The phylogenetic analyses, based on the aa sequences of their RdRP (Figure 4b), showed that "P. palustris" narna-like viruses are grouped in two different clusters. PpaNLV2, 3 (including all 3 variants), 4, 5, and 8 appeared in a detached group, maybe representing a different taxon within the "alphanarnavirus" clade, and closer to Saccharomycces 23S RNA narnavirus and Saccharomycces 20S RNA narnavirus. Within the clade, "betanarnaviruses" PpaNLV1 and 6 appear together with narnaviruses described in Phytophthora spp. and PpaNLV7 groups with two unclassified viruses described in an insect and in symptomatic grapevine tissue associated with the downy mildew Plasmopara viticola.

dsRNA Viruses
Botybirna-like viruses. Four contigs >6 kb in length enclosing one large ORF resembling the RdRP segment of unclassified botybirnaviruses present in fungi and oomycetes ( Figure 5a; Table 1) were found in two pools from Sumatra (SU1, SU3) and one from Kalimantan (KA3). No second segment was detected in our data. These viral contigs were named as Phytophthora palustris botybirna-like viruses (PpaBbLV). The PASC of the nt and aa sequences of the four contigs showed high identity between them (Tables S6 and  S7), with particularly high identity levels (>92%) between the Sumatran contigs. Because the ICTV does not specify virus species demarcation criteria for this family, we classified the four sequences as belonging to two viruses: three contigs represent three Sumatran variants belonging to PpaBbL1 while one represents a Kalimantan botybirna-like virus. CDD belonging to the RT-like superfamily (cl2808) and, particularly, the superfamily member pfam02123 (RdRP_4, which includes RdRPs from Luteovirus, Totivirus and Rotavirus) were also found (Figures 5a and S3; Table 1). The reconstruction of the phylogenetic relationships of PpaBbLV1 and 2 with other dsRNA viruses (Figure 5b) shows that PpaBbLV1 and 2 form a novel cluster separated from other botybirnaviruses found in fungi and oomycetes, suggesting an evolutionary diversification among the mycoviruses of this family.   Toti-like viruses. Virus contigs resembling members the family Totiviridae were also abundant in all pools. After removing redundant and shorter ones, 20 contigs were characterized and named Phytophthora palustris toti-like viruses (PpaTLV). Based on the PASC (Tables S12-S14) and ICTV criterion for the species demarcation within the family Totiviridae, 14 totiviruses were described (Table 1), and those viral contigs with pw identity higher than 90% in both nt and aa sequences were considered variants of the same virus: PpaTLV1 (variants 1 and 2), PpaTLV2 (variants 1, 2, and 3), PpaTLV3 (variants 1, 2, and 3), PpaTLV4, PpaTLV5 (variants 1 and 2), PpaTLV6,7,8,9,10,11, and PpaTLV12 (variants 1, 2, and 3). Most of the contigs consisted of~5-6 kb sequences enclosing two ORFs. Only PpaTVL4 and 8 were shorter and contained one ORF. Conserved domains belonging to the superfamily member pfam02123 (RdRP_4) were found in the 3 -terminus proximal ORF in all the contigs (Figures 5a and S3, Table 1). In addition, conserved domains of the totivirus coat superfamily, cl25797 and pfam05518, were found in the 5 -terminus proximal ORF of PpaTLV1-1, 1-2, 2-1, 2-2, 2-3, 5-1, 5-2, 6, 9, 10, and 11 ( Figure S3, Table 1).
The phylogenetic relationships of PpaTLVs, reconstructed based on their aa RdRP sequences, show a high degree of evolutionary differentiation among them. Whilst PpaTLV1, 2, 5, 6, 10, and 11 appear to cluster with viruses belonging to the genus Victorivirus described in different fungi (Figure 5b), PpaTLV3 falls into the cluster with mycoviruses classified in the genus Totivirus. In a separated cluster, PpaTLV12 is phylogenetically closely related to viruses from Plasmopara viticola, Pythium, and Phytophthora spp. from diverse ecosystems, which form an oomycete-specific cluster related to Giardia canis virus (GCL). PpaTLV4, 7, and 8 cluster with unclassified toti-like viruses reported from fungi, oomycetes (including P. condilina), and diatoms. Comparisons of the geographical distribution of "P. palustris" viruses in Sumatra and Kalimantan show that some viruses have a wider distribution than others (Figures 8 and  9). Thus, five viruses, including three toti-like viruses (PpaTLV3, PpaTLV5, and PpaTLV12) and one narna-like virus (PpaNLV3), even occurred in "P. palustris" isolates from both Sumatra and Kalimantan. PpaTLV3 was found in seven isolates collected at five sites (three in Kalimantan and two in Sumatra). PpaTLV5 was the most widely distributed virus, infecting eight isolates from six different localities (two in Kalimantan and four in Sumatran). Finally, PpaTLV12 infected six isolates from five localities (one from Kaliman- In Sumatra, 9 (SU_F11, F42, R02, R04, R05, R09, R13, R16, R29) out of 30 sampling sites (30%) did not appear to host any of the screened viruses. However, in 21 of the locations (70%), at least one virus was detected. A total of 30 "P. palustris" isolates (58%) hosted at least one virus, and 15 of them showed viral coinfections or multiviral infections. In particular, isolate SU0376 appeared to be infected by seven distinct viruses while SU1474 was infected by six viruses (Figures 7a, 8 and 9). In Kalimantan, isolates from four out the five sites screened (KA_R03, R05, R06, R08) contained viruses. Out of 30 "P. palustris" isolates, 15 isolates hosted at least one virus (50%), and seven of them (47%) showed coinfections and multiviral infections. "P. palustris" isolates KA0119, KA0139, KA0146, and KA0156 hosted at least 11 viruses (Figures 7b, 8 and 9).   The virus richness appears to be greater in Sumatra (23 viruses) than in Kalimantan (17 viruses) (Figures 8 and 9). Moreover, in our samples from Sumatra, the virus family composition was more diverse (Figures 7-9), as a tombus-like virus (PpaTbV1) was found in two isolates collected from two Sumatran sites (SU_R06 and R07). Conversely, with 62 and 66 confirmed viruses, respectively, total virus abundance in "P. palustris" seems to be similar in Kalimantan and Sumatra. However, virus abundance in relation to the number of isolates tested from each of the two islands showed considerable differences: Kalimantan 63/30 = 2.1, Sumatra 66/52 = 1.27. Notably, with 14 different viruses, the Kalimantan site R06 hosted the highest virus diversity of the whole study.
Comparisons of the geographical distribution of "P. palustris" viruses in Sumatra and Kalimantan show that some viruses have a wider distribution than others (Figures 8 and 9). Thus, five viruses, including three toti-like viruses (PpaTLV3, PpaTLV5, and PpaTLV12) and one narna-like virus (PpaNLV3), even occurred in "P. palustris" isolates from both Sumatra and Kalimantan. PpaTLV3 was found in seven isolates collected at five sites (three in Kalimantan and two in Sumatra). PpaTLV5 was the most widely distributed virus, infecting eight isolates from six different localities (two in Kalimantan and four in Sumatran). Finally, PpaTLV12 infected six isolates from five localities (one from Kalimantan and four in Sumatran). Within the same island, PpaTLV11 was present in 12 isolates from two Kalimantan localities, and PpaBLV9 infected 7 isolates obtained from seven Sumatran localities. Finally, the results show that, in general, toti-like viruses are more abundant and diverse in Kalimantan whereas bunya-like viruses are more prevalent and diverse in Sumatra.
According to the average pw identity percentages of the alignment of the amplicon sequences (variants) with the virus contig sequences obtained by NGS (Table S15), some viruses also have more genetic stability than others, independently of their occurrence in two or more "P. palustris" isolates, i.e., the three variants of botybirna-like virus 1 (PpaBbLV1) are more diverse than the four variants of PpaBbLV2.

Virus Read Abundance per Pool and Island
There is a large variation in the total viral read numbers between Japan, Taiwan, and Indonesia (Kalimantan and Sumatra) as well as between specific virus contigs. As shown in Figure S5, Sumatran libraries contained the highest percentage of viral reads, while Japanese and Taiwanese libraries contained the lowest. Likewise, the narna-like viruses constituted the highest read numbers (Table 1, Figure S5), with the PpaNLV 3-1, 3-2, 3-3, and 4 being clearly the most abundant viruses, occurring in libraries KA1, KA2, and SU1-SU3. Bunya-like viruses' reads also appeared to be very abundant, in particular, the three variants of PpaBLV9 (found in libraries SU1-SU4). Interestingly, despite not having been confirmed by RT-PCR, a total of 608 reads were mapped to bunyavirus 8 (PpaBLV8) in all three libraries JP-TW, JP1, and JP2, while 68 reads were mapped to narna-like virus 8 (PpaNLV8). These reads scarcely covered these two virus sequences, repeatedly mapping short regions of those viruses. The abundance of toti-like and tombus-like viruses' reads did not appear to be correlated with their distribution and richness. Compared to the high diversity and distribution of toti-like viruses, the number of reads for each contig was low, while the only tombus-like virus contig, which occurred mostly in SU2, had nearly 200,000 reads.

Virus Read Abundance per Pool and Island
There is a large variation in the total viral read numbers between Japan, Taiwan, and Indonesia (Kalimantan and Sumatra) as well as between specific virus contigs. As shown in Figure S5, Sumatran libraries contained the highest percentage of viral reads, while Japanese and Taiwanese libraries contained the lowest. Likewise, the narna-like viruses constituted the highest read numbers (Table 1, Figure S5), with the PpaNLV 3-1, 3-2, 3-3, and 4 being clearly the most abundant viruses, occurring in libraries KA1, KA2, and SU1-SU3. Bunya-like viruses' reads also appeared to be very abundant, in particular, the three

Discussion
The virus composition, abundance, and distribution discovered in our study indicate that viruses related to the families Totiviridae, Tombusviridae, Narnaviridae, and Botybirnaviridae and the order Bunyavirales are naturally infecting "P. palustris" isolates in their niches and must be readily carried by "P. palustris" propagules, revealing that they are primarily members of the host cellular environment.

Evolutionary Insights into Novel Taxa of (+)ssRNA Viruses
PpaTbLV1 is the first tombus-like virus described in a Phytophthora species, but not the first one in oomycetes. Four tombus-like viruses have been identified in Plasmoparaviticola-associated lesions on grapevines in Italy [43]. Plasmopara halstedii virus (PhV) and Sclerophthora macrospora virus A (SmV-A) have been detected in the downy mildew of sunflower Plasmopara halstedii [55] and the grass pathogen Sclerophthora macrospora [42], respectively, although PhV and SmV-A have been assigned to a new putative viral group between Nodaviridae and Tombusviridae. Based on the BLASTX and the phylogenetic analyses of this study (Figure 3b), PpaTbLV1 does not appear to be closely related to these viruses but to insect viruses and a virus found by a metagenomic approach in soils of wild oat (Avena fatua), an annual grass common in Mediterranean climates [56]. Interestingly, PpaTbLV1 is not grouped either with recognized plant tombusviruses, such as tomato bush stunt virus (TBSV) or with fungal tombus-like viruses (mycotombusviruses), such as Diaporthe RNA virus (DRV) [57]. However, PpaTbVLV1 genomic properties are closer to the latter, as shown by similar genome size (~4.2 kb), apparent absence of movement proteins typical of plant viruses, and the presence of two main ORFs encoding a hypothetical protein of unknown function and RdRP (Figure 3a). The RdRP is presumed to be expressed as a fusion product with the 5 -proximal protein via readthrough of the amber termination codon, UAG, of the upstream ORF [58]. However, PpaTbLV1 possesses GDD as the catalytic triplet in the motif-C of its RdRP ( Figure S2), in contrast with a common feature in mycotombusviruses, which have GDN as the RdRP catalytic triplet [59]. Intriguingly, a smaller antisense ORF encoding an unknown protein is detected proximal to the 5 -terminus, which seems to be a unique property, not yet found in other tombus-like viruses. Thanks to high-throughput sequencing (HTS) studies, ambisense coding genomes emerge more and more commonly within RNA viruses and, particularly, in fungal viruses, such as ambiviruses [60] and "ambinarnaviruses" [61].
PpaTbLV1 was confirmed to be carried by two isolates of "P. palustris" collected in two Sumatran localities, within a distance of a few kilometers (Figures 8 and 9, Table S1). Hence, it seems possible that PpaTbV1 was transported by swimming zoospores and then exchanged via anastomosis between compatible individuals of "P. palustris". Members of the "P. palustris" complex are aquatic oomycetes, and, interestingly, several reports have identified tombusviruses in association with rivers and lakes throughout the world [62,63]. TBSV and Carnation Italian ringspot virus (CIRV), have been isolated from waterbodies draining forest areas in Northrhine-Westfalia, Germany [64]. In addition, tombus-like viruses are also hosted by unicellular protists [65], they have been associated with the holobiont of freshwater shrimp [66], and have been discovered on marine and freshwater RNA viromes using metagenomics [67].
The family Narnaviridae is undergoing an intense reconstruction, as more and more new narna-like viruses with very diverse genomic properties have been discovered lately [58]. Until recently, it was generally understood that narnaviruses had monopartite (+) ssRNA genomes enclosing solely an essential RdRp gene encoded by a single ORF. However, it has been discovered that some narnaviruses have divided RdRPs ("splipalmiviruses"), i.e., viruses hosted by deep-sea fungi, such as Aspergillus tennesseensis. Others have biand multi-segmented genomes (reviewed in [58]) or possess long-reverse-frame ORFs, i.e., Plasmopara-viticola-associated narnaviruses [43]. Although we cannot rule out the possibility of segmented genomes, "P. palustris" narna-like viruses seem to be phylogenetically related to traditional monopartite positive-coding narnaviruses included in both "alpha-narnaviruses" and "betanarnaviruses" clades ( Figure 4b). Likewise, no reverse large ORFs were detected in any of the final contigs. The phylogenetic tree and the pairwise comparisons (Tables S10 and S11) show moderate levels of genetic variability. Whilst the narna-like viruses' richness in "P. palustris" is remarkable (eight putative viruses are described), five of them are closely related and cluster together, suggesting they might have originated from the same "P. palustris" host and gradually incorporated genetic modifications over time [68]), probably during the replication process due to random mutations and genetic adaptations to a new geographical area [69][70][71]. A common ancestor might also explain the evolutionary relationship of PpaNLV1 with a narnavirus found in Phytophthora castaneae in Vietnam [40], and PpaNLV6 with a narnavirus found in P. infestans [34]. Narnaviruses seem to be effectively transmitted throughout "P. palustris" populations and are widespread in both Sumatra and Kalimantan. A certain relation between the high narnavirus abundance and read number obtained, in particular, for PpaNLV2, 3, and 4 (Table 1, Figure S5), agrees with the frequent occurrence of these viruses on both islands.

Discovery of New dsRNA Viruses and High Abundance of Toti-Like Viruses
Kalimantan and Sumatran isolates of "P. palustris" have been found to host two putative botybirna-like viruses (PpaBbLV1 and 2). While genetically highly similar (Tables S6 and S7), each PpaBbLV occurs in several isolates and more than one site but on separate islands, PpaBbLV1 in Sumatra and PpaBbLV2 in Kalimantan (Figures 8 and 9). This constitutes the first discovery of a putative botybirnavirus in a Phytophthora species. Previously, two botybirnaviruses were identified in grapevine lesions associated with Plasmopara viticola in Spain and Italy [43], and in several fungal plant pathogens, including Botrytis porri [72] and Sclerotinia sclerotiorum [73]. Botybirnaviruses have a bipartite genome, with two types of structural proteins (p85/80, p70) and an RdRP. However, only the segment enclosing the RdRP was identified in our study. Both PpaBbLV 1 and 2 appear to be phylogenetically distinct from other mycobotybirnaviruses, representing a potential new taxon. Moreover, aside from motifs I and II, the rest of the RdRP domains closely resemble those of totiviruses ( Figure S3). "P. palustris" isolates are mainly infected by totivirus-related viruses. They are widespread in both Kalimantan and Sumatra (Figures 8 and 9) and novel Bona-fide totiviruses, victoriviruses, and toti-like viruses show paraphyletic relationships across the phylogenetic tree and considerable genetic diversification (Figure 5b). Some viruses reside in an oomycete-specific cluster while others cluster with unclassified toti-like viruses reported from fungi, oomycetes and diatoms dwelling in a wide range of ecosystems. This result may suggest long-term coevolution between "P. palustris" and its totiviruslike viruses.

Persistence of Bunya-Like Viruses and the First Oomycete Bunyavirus with NC Protein
Similar to recent studies on marine oomycetes [33,46], high variability and abundance of the L segment (RNA 1) of putative bunyaviruses were also found in "P. palustris", most of them phylogenetically closely related but with an overall high pairwise genetic variability that allows the description of 13 viruses with different variants. Such virus richness can be the consequence of mutation, recombination and reassortment during the course of evolution. Viral reassortment seems to be a powerful mechanism underlying the evolution of the Bunyavirales order, and it has been pointed out that most bunyaviruses described so far are actually reassortants of extant or extinct viruses [74]. Alternatively, it cannot be ruled out that the P. palustris bunya-like contig sequences would not represent true biological entities. Since bunyaviruses might be reassortants, which may sometimes be seen in their genomes as having highly conserved sequence stretches alternating with highly variable ones, putatively resulting in chimeric contigs in the de novo assembly when they are found abundantly in NGS-pools (e.g., [37]). The phylogenetic analysis (Figure 2a) shows the increasingly complex taxonomy of mycobunyaviruses with strong cluster support.
PpaBLV11 is the first bunya-like virus with bipartite genome reported in Phytophthora. It is related to phleboviruses, which seem to have a tripartite genome (three RNA segments). In this study, we only found two RNAs, RNA1 and RNA2, respectively, corresponding to the L and S segments of phleboviruses. Neither the M segment nor M segment-like contigs were detected during the NGS analysis. If PpaBLV11 has a tripartite genome, the M segment may not be considered essential. This is in accordance with the most similar virus, tulip streak virus (TuSV) [75] and fungal phleboviruses, as Lentinula edodes negative-strand RNA virus 2-HG3 (LeNSRV2-HG3) [76]. However, similar to the rest of the Phytophthora bunyaviruses, which seem to have monopartite genomes, it is more plausible that the genome description of these viruses is incomplete. The putative NC and other non-structural (Ns) associated proteins are likely not conserved enough to be detected by homology.

Virus Population Structure and Transmission Efficiency
Our study indicates that "P. palustris" viral populations adhere to a certain structure. Although the virus family composition is similar on both Kalimantan and Sumatra islands, the virus structure differs and appears to be specific to each island. In Kalimantan, totilike viruses predominate, whereas bunya-like viruses are more prevalent in Sumatra. Furthermore, the absence of viruses in Japanese and Taiwanese "P. palustris" isolates suggests a degree of divergence regarding Indonesian members of the "P. palustris" complex and their viruses. However, it cannot be ruled out entirely that Japanese and Taiwanese isolates also host viruses. In fact, various viral reads have been mapped to PpaNV6 and PpaBLV8 in the RNA pools JP1, JP2, and TW-JP ( Figure S5). If there are any viruses, they are not those found in Kalimantan and Sumatra. Since mycoviruses do not have an extracellular stage, they strongly depend on their host for dispersal, and an analogous evolutionary scenario is expected for both the host ("P. palustris") and the hyperparasite (mycoviruses) [26]. However, due to geographical isolation, local adaptation, and drift, "P. palustris" seems to have evolved into new lineages, and the species P. sp. palustris in Japan and Taiwan is very different from all Indonesian taxa, indicating they have been separated from each other millions of years ago (T. Jung and M. Horta Jung, unpublished data), similarly to evolutionary divergent lineages of P. ramorum [16] and P. lateralis [13]. As a consequence of speciation processes in "P. palustris" hosts, the possibilities of the associated mycoviruses to coevolve and be always co-introduced with their host would significantly decrease. Multiple introduction routes do not necessarily imply that hyperparasites are always introduced along all of them [77]. Alternatively, the absence of viruses in these pools could be also due to the limited availability of similar viral sequences in GenBank; hence, a data re-analysis is highly recommended in a few years.
Our results suggest reasonable efficient viral transmission across Indonesian "P. palustris" populations. However, some "P. palustris" viruses are more widespread than others, merging in distinct Indonesian localities of the same island and even on different islands (PpaNLV3, and PpaTLV3, 4, and 12). Most likely these viruses evolved towards optimal virus-host interactions that enabled them to replicate inside their "P. palustris" hosts and move more efficiently through their host populations [77]. Members of the "P. palustris" complex are aquatic Phytophthora species and, hence, able to produce large amounts of motile chemotactic zoospores [78], which in waterbodies can be transported passively over long distances. However, "P. palustris" shares its aquatic habitat with nematodes, algae, crustaceans, and mollusks [79], with which they undoubtedly interact, and which might eventually serve as transport means of potential virus-hosting oomycete individuals. On the contrary, the lower occurrence of other viruses might be related to the lower efficiency of interspecies virus transmission, because five species of the "P. palustris" complex occur in Sumatra and Kalimantan (T. Jung and M. Horta Jung, unpublished data). This study did not observe a virus/"P. palustris" species correlation but it should be further investigated. Previous studies concerning virus transmission in P. infestans demonstrated 100% inheritance of Phytophthora infestans RNA virus 3 (PiRV3) in individual zoospores [80]. Phytophthora infestans RNA virus 2 (PiRV2) was readily transmitted by hyphal anastomosis, and by asexual reproduction through sporangia [80]. However, attempts to transfer PiRV-2 into apparently vegetatively incompatible P. infestans isolates failed.

Multiple Viral Infections at Isolate and Site Level
Particular Indonesian localities and "P. palustris" isolates appear to be big reservoirs of viruses (Figures 7-9). From an ecological perspective, in those localities with higher virus diversity and abundance, the number of contact events between and among host individuals must be higher [26]. Thus, the probability of a virus-hosting individual having contact with another compatible virus-hosting "P. palustris" isolate increases with host population density. From an individual perspective, multiple viral infections are, indeed, known to occur in oomycetes in nature [23,37,40,43,46]. However, the interactions between different viruses and potential limits for the accumulation of multiple viruses in a single host remain unclear. In a native environment, long-term co-existence between virus and host should lead to a non-lethal equilibrium. Thus, the cytoplasmic exchange that takes place during fungal (and oomycete) growth promotes the accumulation of multiple viruses within a hypha, colony, or hyphal network [81]. Unlike fungi, oomycete hyphae lack septa, potentially facilitating the accumulation and exchange of viruses within and between compatible individuals. In addition, many Phytophthora species, in particular those thriving in aquatic environments, are prompted to interspecific hybridizations, which play a major role in speciation and species radiations in diverse natural ecosystems [82][83][84]. Interspecific hybridizations might enable virus transmission between different Phytophthora species.
Supplementary Materials: The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/jof8111118/s1, Table S1. List of "Phytophthora palustris" isolates used in the present study, Table S2. Primers used for the confirmation of each virus presence in the collection of "Phytophthora palustris" isolates by RT-PCR and RACE, Table S3. Trimmed low quality and adaptor data, Table S4. Mapped reads to Phytophthora parasitica genome (used as a reference host genome), Table S5. Data of the resulting contigs obtained by the de novo assembly in each sequenced RNA library. Pairwise identities of the RdRp nucleotide (Table S6) and amino acid (Table S7) sequences of "Phytophthora palustris" botybirnaviruses, Table S8. Pairwise identities of the RdRp amino acid sequences of "Phytophthora palustris" bunyaviruses, Table S9. Pairwise identities of the RdRp nucleotide sequences of "Phytophthora palustris" bunyaviruses, Table S10. Pairwise identities of the RdRp nucleotide sequences of "Phytophthora palustris" putative narnaviruses, Table S11. Pairwise identities of the RdRp amino acid sequences of "Phytophthora palustris" narnaviruses, Table S12. Pairwise identities of the RdRp nucleotide sequences of "Phytophthora palustris" toti-like viruses, Table S13. Pairwise identities of the RdRp amino acid sequences of "Phytophthora palustris" toti-like viruses, Table S14. Pairwise identities of the CP amino acid sequences of "Phytophthora palustris" toti-like viruses, Table S15. Average pairwise (PW) identity of the alignment of the partial viral sequences (variants) with the virus contig sequences obtained by NGS, Figure S1. A. Amino acid alignment showing Bunya_RdRp conserved motifs A to E and premotif A within the RdRp of PpaBLV1-14 and selected bunyaviruses, B. Amino acid alignment showing the conserved motifs of Tenuivirus/Phlebovirus nucleocapsid protein within PpaBLV11 and selected bunyaviruses, Figure S2. Amino acid alignment showing the conserved RdRp conserved motifs of PpaTbVL1 and selected tombusviruses, Figure S3. Conserved aa sequence motifs (I-VIII) of RdRp of totiviruses and botybirnaviruses in P. palustris, Figure S4. Conserved aa sequence motifs (III-VII) of RdRp of narna-like viruses, Figure S5