Abstract
We investigated the incidence of RNA viruses in a collection of Halophytophthora spp. from estuarine ecosystems in southern Portugal. The first approach to detect the presence of viruses was based on the occurrence of dsRNA, typically considered as a viral molecule in plants and fungi. Two dsRNA-banding patterns (∼7 and 9 kb) were observed in seven of 73 Halophytophthora isolates tested (9.6%). Consequently, two dsRNA-hosting isolates were chosen to perform stranded RNA sequencing for de novo virus sequence assembly. A total of eight putative novel virus species with genomic affinities to members of the order Bunyavirales were detected and their full-length RdRp gene characterized by RACE. Based on the direct partial amplification of their RdRp gene by RT-PCR multiple viral infections occur in both isolates selected. Likewise, the screening of those viruses in the whole collection of Halophytophthora isolates showed that their occurrence is limited to one single Halophytophthora species. To our knowledge, this is the first report demonstrating the presence of negative (−) ssRNA viruses in marine oomycetes.
Introduction
Halophytophthora species are fungal-like oomycetes with similar morphology and life cycles as members from their well-known plant pathogenic sister genus Phytophthora (Sullivan et al., 2018). As oomycetes they are distantly related to brown algae and belong to the Kingdom Stramenopila (Heterokonta). In pre-molecular times most marine oomycetes were assigned to the genus Halophytophthora but recent phylogenetic studies demonstrated that Halophytophthora is polyphyletic. Consequently, numerous Halophytophthora species were transferred to several new genera including Calycofera, Phytopythium, and Salisapilia (; ; ; ). Currently Halophytophthora sensu stricto contains nine described species. Most Halophytophthora species live in brackish and salt water habitats and have been traditionally described as saprophytes playing a key role as decomposers mainly in mangrove ecosystems (, ). However, they may also be pathogenic. Several studies have already illustrated their pathogenicity on the marine eelgrass Zostera marina (), which has a key ecological role along shores of North America and Eurasia (). Halophytophthora zostera was shown to restrict the viability of Z. marina seeds and seedling development (). Knowledge on Halophytophthora in temperate ecosystems has been scarce for a long time, hence, the conditions explaining why these decomposers may turn pathogenic are still unknown (Sullivan et al., 2018).
Virus diversity in marine ecosystems appears predominant (Suttle, 2005). With a total estimation of ∼1030 viruses, they have been implicated as the most abundant pathogens in the oceans (Suttle, 2005). Viruses infect all organisms from bacteria to whales but the majority infect bacteria (phages), which control microbial abundance and release dissolved organic matter, influencing global biogeochemical cycles (Wommack and Colwell, 2000; Weinbauer, 2004; ). Regarding lower eukaryotes, metagenomic and metatranscriptomic studies are rapidly expanding the knowledge about virus diversity and can detect the presence of a range of DNA or RNA genomes of various architectural type and size (). Viruses with large double-stranded (ds) DNA genomes are typically found in phytoplankton (Weynberg et al., 2017). Single-stranded (ss) DNA viruses have been described in diatom algae (; ), corals (Vega Thurber et al., 2008) and stromatolites (). DsRNA viruses infect photosynthetic flagellates, heterotrophic protists () and also fungi. Recently, dsRNA viruses hosted by marine fungi were isolated from the seagrass Posidonia oceanica () and sea cucumber Holothuria polii (), and identified a negative (−) ssRNA virus in Penicillium roseopurpureum.
The order Bunyavirales is comprised of viruses with segmented, linear (−) ssRNA genomes (). The morphology of Bunyaviruses varies from symmetric icosahedral particles, similar to phleboviruses (), to spherical and pleomorphic particles, similar to orthobunyaviruses () and hantaviruses (). A model bunyavirus genome consists of three RNA segments, the smallest of which (S) may encode nucleocapside proteins (NP) and silencing suppressors (NS) while the medium sized (M) segment usually encodes glycoproteins (Gn and Gc) and movement proteins (). Finally, the largest segment (L) appears to possess a single open reading frame (ORF), which encodes the RNA dependent RNA polymerase (RdRp). According to the last update from the International Committee on Taxonomy of Viruses (ICTV), there are 12 bunyaviral families, four subfamilies, 46 genera and 287 species registered (). However, due to more intensive samplings and, in particular, because of the use of de novo virus detection with metagenomics sequencing the diversity of bunyaviruses is constantly increasing (; Shi et al., 2016). Currently, bunyaviruses have been described in plants, invertebrate and vertebrate hosts and are transmitted by arthropod and mammalian vectors (). Viruses with genomic similarities to bunyaviruses have also been described from diverse fungal hosts, including the phytopathogenic fungi Botrytis cinerea (Botrytis cinerea negative-stranded RNA virus 1, BcNSRV-1) () and Macrophomina phaseolina (Macrophomina phaseolina negative-stranded RNA virus 1, MpNSRV1) (). Entoleuca phenui-like virus 1 (EnPLV1) has been identified from an avirulent isolate of Entoleuca sp. collected from avocado rhizosphere (Velasco et al., 2019). Coniothyrium diplodiella negative-stranded RNA virus 1 (CdNsRV1), Alternaria tenuissima negative-stranded RNA virus 2 (AtNsRV2) and Cladosporium cladosporioides negative-stranded RNA virus 2 (CcNsRV2), were isolated from the grapevine wood-inhabiting endophytes Coniothyrium diplodiella, Alternaria tenuissima, and Cladosporium cladosporioides, respectively (). Likewise, the shiitake mushroom (Lentinula edodes) hosts Lentinula edodes negative-stranded virus 2 (LeNSRV2), which is related to phenuiviruses (); a marine fungus, P. roseopurpureum hosts Penicillium roseopurpureum negative ssRNA virus 1 (PrNSRV1) () and the oomycete Pythium polare, Pythium polare bunya-like RNA virus 1 (PpBRV1) (). However, no effects of bunyavirus infection on growth or virulence of their fungal or oomycete hosts have been recorded yet.
Knowledge on viruses infecting oomycetes is relatively limited (Sutela et al., 2019). Several viruses have been described in downy mildews (biotrophic plant parasites) including unclassified (+) ssRNA viruses on Sclerophthora macrospora with similarities to a Noda-Tombus-like virus (Yokoi et al., 2003) and in Plasmopara parasitica, which causes hypovirulence on its host (; ). The genus Pythium is comprised of water and soilborne oomycetes causing moderate to significant damages in plant roots (Sutela et al., 2019). Virus-like particles and/or dsRNA have been described in Pythium irregular (). Recently, an unclassified gammapartitivirus was reported in Pythium nuun (Shiba et al., 2018), a toti-like virus was also characterized from two strains of Globisporangium splendens (formerly Pythium splendes) and three virus-like sequences, Pythium polare RNA virus 1 (PpRV1), Pythium polare RNA virus 2 (PpRV2) and PpBRV1 were detected in Pythium polare infecting mosses in the Arctic (). In the genus Phytophthora the first virus to be reported in the USA () was classified as an alphaendornavirus. It was found in an isolate of the undescribed Phytophthora taxon “douglasfir.” Later, similar virus strains were detected in P. ramorum isolates from several hosts in Europe () but no further investigations were performed. Phytophthora infestans, the causal agent of potato late blight, also harbors four RNA viruses. Phytophthora infestans RNA virus 1 and 2, (PiRV-1 and PiRV-2, respectively) apparently represent novel virus families (, ), Phytophthora infestans RNA virus 3 (PiRV-3) is clustered with the newly proposed family “Fusagraviridae” () and Phytophthora infestans RNA virus 4 (PiRV-4) is an unclassified member of Narnaviridae (). Recently, PiRV2, known to be 100% transmittable through asexual spores, was shown to stimulate sporangia production and enhance the virulence of P. infestans (). Phytophthora cactorum RNA virus 1 (PcRV1), a toti-like virus, has been newly described in an isolate of Phytophthora cactorum collected from a trunk lesion on silver birch in Denmark ().
Because the genus Halophytophthora represents a unique group of marine microorganisms with uncertain roles as phytopathogens on coastal and estuarine grasses, we wished to expand current research to better understand their behavior by investigating their potential virome. More precisely, our study had two goals: (i) to confirm the presence of viruses in these marine oomycetes by combining traditional and state-of-art technologies, and (ii) to assess their abundance and genetic variability in a collection of Halophytophthora isolates from Portugal.
Materials and Methods
Oomycete Isolates
All Halophytophthora isolates studied were collected from seven localities in southern Portugal (Table 1) using an in situ baiting technique (). At each site, 15–20 non-wounded young leaves of three tree species, Ceratonia siliqua, Citrus sinensis and Quercus suber, were placed as baits in a 25 × 30 cm raft, prepared using fly mesh and styrofoam, and the raft put to float in the tidal zone. The rafts were collected after 3 days. Baiting leaves were washed in distilled water and blotted dry on filter paper. Five to ten pieces (approximately 2 × 2 mm) were cut from the margins of each watersoaked or necrotic lesion of each leaf, blotted on filter paper and plated onto selective PARPNH agar (V8-juice agar (V8A) amended with 10 μg/ml pimaricin, 200 μg/ml ampicillin, 10 μg/ml rifampicin, 25 μg/ml pentachloronitrobenzene (PCNB), 50 μg/ml nystatin and 50 μg/ml hymexazol). The Halophytophthora isolates are part of the collection of the Phytophthora Research Centre (Mendel University), located in Brno, Czech Republic. The Supplementary Table S1 provides detailed information and species identification of all Halophytophthora isolates, which is currently under investigation. For the different purposes of the present study all isolates were grown for 7–21 days in darkness on V8A media, covered with cellophane (EJA08-100; Gel Company, Inc., CA, United States).
TABLE 1
| Locality, municipality | N& | dsRNA screening$ | dsRNA patterns | RT-PCR screening+ | Virus found by RT-PCR* |
| Ribeira de Odelouca, Silves | 14 | 14 | ∼9 kb | 14 | None |
| Rio Séqua, Tavira | 9 | 9 | None | 9 | None |
| Parque Natural da Ria Formosa, Santa Luzia, Tavira | 24 | 18 | ∼7 and 9 kb | 21 | HRV1-8 |
| Parque Natural da Ria Formosa, Quelfes, Olhão | 6 | 5 | ∼9 kb | 6 | None |
| Ria de Alvor, Alvor, Portimão | 8 | 6 | None | 8 | None |
| Parque Natural da Ria Formosa, Almancil, Loulé | 9 | 8 | ∼9 kb | 8 | None |
| Sapal de Castro Marim/Rio Guadiana, Castro Marim | 32 | 13 | None | 29 | None |
Relation of isolates per sampling locality and virus detection.
&Total number of isolates screened for viruses by both dsRNA isolation and partial amplification of their RdRp by RT-PCR using specific primers for each of them; $number of isolates used in the dsRNA screening; +number of isolates screened by RT-PCR; *viruses detected by RT-PCR. See Supplementary Table S1 for the complete information of the isolates, DsRNA and RT-PCR screening.
DsRNA Isolation
A total of 73 isolates from Portugal were screened for dsRNA (Figure 1), which was purified using a modified version of the protocol of . Approximately 2 g of fresh mycelium were transferred to a 50 ml Falcon tube and disrupted by vortexing for 3 min with two stainless steel balls (diameter 10 mm). The rest of the protocol was performed as described by Tuomivirta et al. (2002).
FIGURE 1
Stranded RNA Sequencing of Samples BD641 and BD647
Total RNA was purified from approximately 50 mg of fresh mycelium using SPLIT RNA Extraction Kit (Lexogen, Austria) and treated with DNase I (ThermoFisher Scientific). RNA quantity and quality were checked using respectively a Qubit® 2.0 Fluorometer (Invitrogen) and Tape Station 4200 (Agilent) resulting in a RNA integrity number (RIN) of 10.
Approximately 1 μg of total RNA eluted in RNase-free water was sent to Fasteris SA (Plan-les-Ouates, Switzerland) for RNA library construction and deep sequencing. The library preparation was performed using the Illumina TruSeq® stranded RNA Sample Preparation Kit (Illumina, San Diego, CA, United States). The library was sequenced in pair-end (2x 75 nt) runs on an Illumina NextSeq 500 machine. Prior to the library preparation, Ribo-Zero rRNA Removal Kit (Human/Mouse/Rat) was used successfully. An “in-lane” PhiX control spike was included in each lane of the flow-cell. 93.94% of the reads had a quality value (Q30) ≥ 30, i.e., less than 1 error in 1000 bases. The raw data was deposited in the Sequence Read Archive (SRA) with the BioProject ID PRJNA619952 and BioSample ID SAMN14856044.
Bioinformatics Analysis
Adapter trimming was not needed as the adapter sequence was only present in less than 1% of the reads. All reads were compared to the genome sequence of Phytophthora cinnamomi1 with BWA 0.7.5a2. Unmapped reads were selected from the mapping files and saved in a fastq format. BAM post-processing was performed with toolbox for manipulation of SAM/BAM files V. 1.13 and BEDTOOLS V. 2.21.04. De novo assembly was performed with VELVET V1.2.105. As the repeat resolution module of VELVET assumes linearity and uniform coverage distribution, it produces fragmented transcriptome assemblies. To consider the unequal expression levels and alternative splicing breakpoints, the preliminary assemblies produced by VELVET were inputted to OASES 0.2.086, which exploits read sequence and pairing information (if available) to produce transcript isoforms. When possible, OASES also detects and reports standard alternative splicing events. The first 1M paired-reads of each library were mapped on each OASES assembly. The alignment was done using the mapping software BWA 0.7.5a and SAMTOOLS 1.17. Reads mapping to several positions on the reference sequence with the same mapping quality were attributed at random to one of the positions with a mapping quality of 0. When an input read had N’s in their nucleotide sequence, BWA replaced the Ns with a random nucleotide. De novo assembly and validation mapping were done using BWA 0.7.5a. The result of the mapping of the reads on the reference sequences is summarized in Supplementary Table S2. The assemblies with a hash of 51 gave the highest representability. The numbers of reads aligning and coverage depth for the final viral sequence were calculated using Geneious Prime® 2020.0.4.
The sequence alignment and blast output post-processing with BLAST ncbi-blast-2.2.267. The contig files were aligned to the NCBI database using a local installation of the BLAST software. The results from the BLASTn, BLASTX and the BLAST search on the viral reference dataset (RefSeq) were compared. The output file was then parsed using in-house (Fasteris) scripts.
Rapid Amplification of cDNA Ends (RACE) and Full-Length Viral Sequences
In order to determine the 5′- and 3′- terminal sequences and lengths of the putative viral genomes the SMARTer® RACE 5′/3′ KIT (Takara Bio USA, Inc.) was used with total RNA extracted as described above using SPLIT RNA Extraction Kit (Lexogen, Austria). The 3′ Poly(A) Tailing of RNA was performed at 37°C for 20 min using yeast Poly(A) Polymerase (MCLAB, San Francisco, CA). Thereafter, approximately 1 μg of total RNA was used for RACE First-strand cDNA synthesis of both 5′- and 3′- termini was performed using specific primers designed for the eight putative viruses in a 5′- and 3′ – orientation (Supplementary Table S3) as described by the manufacturer. Then, 5′-RACE and 3′-RACE PCR amplification was performed to generate the corresponding cDNA fragments using SeqAmp DNA Polymerase (Takara Bio USA, Inc.) as described by the manufacturer. Amplicons were extracted and purified from the gel using NucleoSpin® Gel and PCR Clean-Up Kit (Macherey-Nagel GmbH & Co. KG). All amplicons were cloned using In-Fusion HD Cloning Kit and Stellar Competent Cells (Takara Bio USA, Inc.). Recombinant plasmids were extracted with Thermo Scientific GeneJET Plasmid Miniprep Kit and sequenced by GATC Biotech, Germany. Nested PCR was often necessary in order to precisely amplify the virus 3′ termini. Here, the screening forward primers (Supplementary Table S3) in combination with the Universal Primer mix (UPM) from the SMARTer® RACE 5′/3′ KIT (Takara Bio USA, Inc.) were used to generate amplicons to act as template for nested PCR. The nested PCR RACE primers and the Short Universal primer provided by the kit were used. To obtain the 3′ end terminal sequence of virus 8 a different primer was designed.
Screening for Virus Incidence by RT-PCR Amplification
Mycelium from individual isolates was collected and homogenized using a bead tube holder devide (740469; Macherey-Nagel; Germany). RNA was isolated and purified using the Monarch Total RNA Miniprep Kit (T2010S; New England Biolabs, MA, United States). Efficient RNA extraction was achieved using the manufacturer’s recommended protocol for tough-to-lyse samples. Contaminating host DNA was removed from the extracts using a combination of gDNA removal columns and DNase I treatment. RNA was eluted in 30 μl volumes and stored at −80°C.
cDNA was synthesized using the ProtoScript II First Strand cDNA Synthesis Kit (E6560; New England Biolabs, MA, United States). Here, oligo d(T)23VN was incubated with the RNA at 65°C for 5 min. Next, random primer mix was added together with ProtoScript II Enzyme and Reaction Mix, followed by incubation at 25°C for 5 min, 42°C for 60 min and denatured at 80°C for 5 min.
Virus incidence in the Halophytophthora ssp. isolates was tested using PCR amplification with virus-specific primers (Supplementary Table S4). All primers were used to amplify a fragment of the RdRp gene of each virus and were designed by Primer 3 2.3.7 under Geneious Prime® 2020.0.4. Primer sequences, amplicon sizes and annealing temperatures are shown in Supplementary Table S4. PCR amplification was performed with 12.5 μl OneTaq Quick-Load 2X Master Mix with Standard Buffer (M0486; New England Biolabs, MA, United States), 0.5 μl of each 10 mM primer, 3 μl of cDNA in a total volume of 25 μl. Cycling conditions were used according to manufacturer’s recommendations and the annealing temperature of each primer set was calculated using the on-line tool8 (v1.9.9 May 30).
Amplicons were visualized and separated by electrophoresis (300 V; 10 min) through 1.5% agarose gels in TBE buffer (106177; Merck KGaA, Germany) and with DNA Stain G (39803; SERVA; Germany). Amplicons of the expected lengths were purified and sequenced in both directions by GATC BioTech (Eurofins; Konstanz, Germany) with the primers used in the initial PCR amplification. All the amplicon sequences were deposited in the GenBank under the accession number MT277331-349. As an internal control for successful PCR amplification from viral RNA templates routine amplification of an actin housekeeping gene (Weiland and Sundsbak, 2000), was performed simultaneously with primers MIDFWACT and MIDREVACT (Supplementary Table S4) in all experiments under identical conditions (Supplementary Figure S1).
Genetic Variability and Phylogenetic Analysis
Pairwise identities of the nucleotide and amino acid sequences (Supplementary Tables S5, S6) were obtained after aligning the eight viral nucleotide and amino acid sequences by MUSCLE () and calculated using Geneious Prime® 2020.0.4. DnaSP v5 () was used to estimate genetic diversity parameters for strains of the different Halophytophthora viruses.
In order to search for conserved domains within the putative viral proteins the NCBI CD-search tool was used9 (). Viral protein sequences were aligned by MUSCLE () using Geneious Prime® 2020.0.4 (Figure 2).
FIGURE 2
A maximum likelihood phylogenetic tree was constructed using a rapid bootstrapping algorithm (Stamatakis et al., 2008) in RAxML-HPC v.8 on XSEDE conducted in CIPRES Science Gateway () (Figure 3). 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 the substitution model for proteins. Bootstrapping was configurated with the recommended parameters provided by CIPRES Science Gateway. The resulting data were visualized using the software FIGTREE software version 1.4.410.
FIGURE 3
Results
Identification of the Mycoviruses Infecting Selected Halophytophthora Isolates
A total of 7 out of 73 Halophytophthora isolates (9.6 %) from three different Portuguese localities were found to contain one or two dsRNA segments as illustrated by gel electrophoresis (Table 1 and Supplementary Table S1). Those isolates which contained dsRNA apparently belong to three different previously undescribed Halophytophthora species (Supplementary Table S1). Isolates BD093 and BD094 (Halophytophthora sp. 01) from Ribeira de Odelouca (Silves) contained a ∼9 kb-dsRNA segment. Isolates BD641, BD647 and BD654 (H. sp. 04) from Parque Natural da Ria Formosa (Tavira) contained two different banding patterns with ∼7 and/or ∼9 kb dsRNA bands (Figure 1). BD665 (H. sp. 03) from Parque Natural da Ria Formosa in Olhão and BD685 from Parque Natural da Ria Formosa in Almancil, Loulé, both contained a ∼7 kb segment. Two of these positive isolates, BD641 and BD647, were chosen for stranded RNA sequencing.
One sequencing library was prepared from rRNA-depleted total RNA and generated 3.6 × 1011 75-nucleotide (nt) paired end (PE) sequence reads (Supplementary Table S2). When they were assembled against the genome of P. cinamomi as host reference sequence and 1.15 × 108 PEs mapped. The unmapped 1.12 × 108 PE reads were selected for the further investigation. The final contig file had 35,621 contigs (including 548 undetermined bases) with an average length of 894 nt and a maximum contig length of 21,602 nt. Assembly of the sequence reads from the final contig file revealed eight potential virus sequences as evidenced by significant E-values and identity percentage of their predicted amino acid sequences (Table 2). The contigs representing the eight potential viruses had different read numbers and coverage (Table 3), HRV8 had the deepest coverage and HRV2 the least. The eight sequences identified were similar to members of the order Bunyavirales hosted by fungal pathogens, including Botrytis cinerea (Botrytis cinerea negative-stranded RNA virus 1, BcNSRV-1, YP_009182153), Macrophomina phaseolina (Macrophomina phaseolina negative-stranded RNA virus 1, MpNSRV1, ALD89106.2), oomycetes, such as Pythium polare (Pythium polare bunya-like RNA virus 1, PpRV1, YP_09551341.1) and arthropod-like insects including Asellus sp. (Wuhan insect virus 3, AJG39263.1) and crabs (Beihai sesarmid crab virus 5, APG79283.1) (Table 2). According to the convention of the International Committee on Taxonomy of Viruses (ICTV) the eight putative viruses were designated as Halophytophthora RNA virus (HRV) 1–8 (Table 2).
TABLE 2
| Virus name | Acronym | GenBank accession numbers | Most similar virus in GenBank | E value | Query cover (%) | Identity (%) |
| Halophytophthora RNA virus 1 | HRV1 | MT277350 | Beihai sesarmid crab virus 5 | 0.0 | 72 | 32.73 |
| Halophytophthora RNA virus 2 | HRV2 | MT277351 | Botrytis cinerea negative-stranded RNA virus 1 | 0.0 | 52 | 32.40 |
| Halophytophthora RNA virus 3 | HRV3 | MT277352 | Botrytis cinerea negative-stranded RNA virus 1 | 0.0 | 54 | 32.13 |
| Halophytophthora RNA virus 4 | HRV4 | MT277353 | Pythium polare bunya-like RNA virus 1 | 0.0 | 82 | 38.35 |
| Halophytophthora RNA virus 5 | HRV5 | MT277354 | Pythium polare bunya-like RNA virus 1 | 0.0 | 84 | 37.84 |
| Halophytophthora RNA virus 6 | HRV6 | MT277355 | Pythium polare bunya-like RNA virus 1 | 0.0 | 86 | 38.30 |
| Halophytophthora RNA virus 7 | HRV7 | MT277356 | Macrophomina phaseolina negative-stranded RNA virus 1 | 7e–126 | 68 | 28.76 |
| Halophytophthora RNA virus 8 | HRV8 | MT277357 | Wuhan Insect virus 3 | 0.0 | 60 | 58.07 |
Identification of Halophytophthora viruses’ most similar RdRp sequences in the GenBank based on BLASTX search.
TABLE 3
| Initial contig length (nt) | Genome size (nt) | Largest ORF (nt) | Largest ORF (aa) | Function | Read Number | Coverage | |
| Halophytophthora RNA virus 1 | 9,341 | 9,340 | 8,184 | 2,728 | RdRp | 630,934 | 5,066 |
| Halophytophthora RNA virus 2 | 9,134 | 9,152 | 8,361 | 2,787 | RdRp | 15,298 | 126 |
| Halophytophthora RNA virus 3 | 9,113 | 9,184 | 8,985 | 2,995 | RdRp | 716,385 | 5,850 |
| Halophytophthora RNA virus 4 | 7,794 | 7,822 | 7,716 | 2,572 | RdRp | 341,627 | 3,275 |
| Halophytophthora RNA virus 5 | 7,103 | 7,735 | 6,681 | 2,227 | RdRp | 376,150 | 3,647 |
| Halophytophthora RNA virus 6 | 6,326 | 7,742 | 5,796 | 1,932 | RdRp | 843,893 | 8,175 |
| Halophytophthora RNA virus 7 | 5,644 | 6,789 | 6,705 | 2,235 | RdRp | 264,248 | 2,919 |
| Halophytophthora RNA virus 8* | 3,299 | 7,874 | 7,788 | 2,596 | RdRp | 720,641 | 6,864 |
Parameters of the genome organization of Halophytophthora viruses and their initial contigs.
*The length of the 3′ end remains unsure.
Virus Genome Organization
Based on our analyses, only ORFs encoding an RdRp gene were found.
Based on the sequence of the original contigs oligonucleotide primers were designed to amplify and confirm the terminal RdRp sequences of the eight viruses (Supplementary Table S4) some of which, including HRV 1, 2, 3, 4, and 5, were complete or nearly complete. The remainder of the viral terminal sequences of HRV6, 7, and 8 were completed by RACE. The eight virus genomes ranged in size from 7.8 up to 9.3 kb (Table 3) including a single large ORF encoding an RdRp flanked by 5′ and 3′-untranslated regions (UTRs) (Figure 4). The longest genome corresponds to HRV1 (9,340 bp) and the shortest to HRV6 (6,816 bp) (Figure 3 and Table 3).
FIGURE 4
Based on the amino acid (aa) sequence analysis of the eight viral ORFs all of them contained conserved regions belonging to pfam04196, Bunyavirus RdRp, which is the solitary member of the superfamily cl20265. Taken in turn, HRV1 has a conserved region ranging from aa 1454–1759 (expect value 1.01 e– 11); HRV2 from aa 1036–1394 (E-value: 9.37 e–13); HRV3 from aa 1239–1588 (E-value: 5.11 e–10); HRV4 from aa 735–1208 (E-value: 1.51 e–07); HRV5 from aa 716–1301 (E-value: 1.53 e–07); HRV6 from aa 1132- 1421 (E-value: 4.89 e–06); HRV7 from aa 1204-1771 (E-value: 6.24e–10); HRV8 from aa 894 to 1515 (E-value: 1.27 e–30) and from aa 615 to 806 (E-value: 7.55 e–06). The alignment of the RdRps of HRV1-8 with other viral sequences indicated the presence of the typical conserved motifs in RdRps of bunyaviruses () with certain variability (Figure 2): motif A (DxxxWx), where HRV1-3 had an arginine (R) instead of lysine (K) before the tryptophan (W), motif B (XGxxNxxSS), motif C (SDD), motif D (KK) and motif E (ExxSx). Premotif A with the three basic residues inside (K, R and R/K) and, downstream, the glutamic acid (E), were also identified. In addition, the conserved aa triplet TPD, typical of bunyaviruses (), was identified in HRV1-7 (but not in HRV8) in positions 160, 107, 321, 83, 87, 282, and 160, respectively. The doublet RY was strictly conserved in HRV1, 4-8 in positions 894, 754, 757, 970, 1212 and 894, respectively, and partly conserved in HRV2 and 3 (positions 205 and 1019, respectively), which had R (arginine) and N (asparagine) instead of Y (tyrosine).
Phylogenetic Relationships Between the Halophytophthora Viruses and Other (−) ssRNA Viruses
An examination of the phylogenetic relationships between the Halophytophthora viruses and other (−) ssRNA viruses retrieved from the GenBank shows that they cluster with unclassified viruses with genomic affinities of the order Bunyavirales (Figure 3). HRV1 to 7 are very closely related to one another and form a cluster with MpNSRV1, BcNSRV-1 and PpRV1. This cluster includes other bunyaviruses described in invertebrate species representing different metazoan phyla including sesarmid crabs (Beihai sesarmid crab virus 5) and nematodes (Soybean cyst nematode bunya-like virus 1, SCN-BLV1). HRV8 differs from HRV1-7 and is grouped in a different cluster with a virus apparently hosted by Aselus sp. insects (Wuhan insect virus 3), which are phylogenetically closer to reported members of bunyavirus families including Tospoviridae, Fimoviridae, Peribunyaviridae, Hantaviridae, and Phenuiviridae. All HRV1-8 appeared distanced from members of the family Mymonaviridae in the order Mononegavirales.
HRV Occurrence in the Halophytophthora Isolate Collection
A total of four out of 95 Halophytophthora isolates (4.2%) from one sampling locality (Parque Natural da Ria Formosa, Santa Luzia, Tavira) resulted to host one or more of the eight Halophytophthora viruses described in this study. The most abundant was HRV6, which was present in 4 isolates (BD641, BD647, BD650, and BD654); HRV1 and HRV8 in 3 isolates (BD641, BD647, and BD654). HRV7 in isolates BD641 and BD647; HRV3 in isolates BD641 and BD654; HRV5 is in isolates BD647 and BD654, and HRV2 and HRV4 seem to be only present in isolate BD647.
While isolates BD093, BD094, BD665, and BD685 contained dsRNA elements (results not resulted) they were unrelated to the bunyaviral dsRNAs described as determined by RT-PCR (Supplementary Table S1). In some cases, e.g., HRV6 in BD650, dsRNA elements were not visualized by gel staining but RT-PCR confirmed the presence of a virus (Supplementary Table S1). Isolates BD641 and BD647 contained two and one dsRNA elements, respectively (Table 1 and Figure 1), which have not been correlated with HRV1-8 by RT-PCR using dsRNA as a template.
Genetic Diversity of HRV
At the nucleotide level (Supplementary Table S5), the highest pairwise identity between viruses was HRV2 and 3 (57.05%) followed by HRV2 and 3 (39.00%). Conversely, the lowest pairwise identity was between HRV1 and 7 (23.43%). At the amino acid (aa) level (Supplementary Table S6), HRV2 and 3 also have the highest pairwise identity (47.48%), followed by HRV5 and 4 (39.91%). And, HRV7 and 8, the lowest (9.00%).
The genetic variability of the partially sequenced RdRp genes of the eight Halophytophthora virus strains was also assessed for those viruses that were hosted in more than one isolate (Table 4). All of the isolates appear to have very high haplotype diversity with HRV5 possessing the highest nucleotide diversity (0.12) and segregating sites (80). Conversely, HRV7 seems to be the least genetically diverse virus with a low nucleotide diversity (0.05) and the lowest number of segregating sites (34).
TABLE 4
| Species | N | Sites | Net sites | S | h | Hd | π | AvNumDif |
| HRV1 | 3 | 708 | 704 | 51 | 3 | 1.00 | 0.05 | 34 |
| HRV2 | 1 | – | – | – | – | – | – | – |
| HRV3 | 2 | 701 | 700 | 52 | 2 | 1.00 | 0.07 | 52 |
| HRV4 | 1 | – | – | – | – | – | – | – |
| HRV5 | 2 | 655 | 654 | 80 | 2 | 1.00 | 0.12 | 80 |
| HRV6 | 4 | 523 | 523 | 50 | 4 | 1.00 | 0.05 | 28 |
| HRV7 | 3 | 473 | 473 | 34 | 3 | 1.00 | 0.05 | 23 |
| HRV8 | 3 | 659 | 652 | 65 | 3 | 1.00 | 0.07 | 44 |
Genetic parameters of the partial RdRp nucleotide sequences of HRV strains.
(–) not assessed; (N) Number of isolates; (Sites) length of the sequence, (Net sites) length of the sequence analyzed; (S) number of segregating sites; (h) number of haplotypes; (π), nucleotide diversity estimated by the average number of differences per site between two sequences; (AvNumDif) average number of differences.
Discussion
Traditional dsRNA extraction procedures were used to identify the potential presence of viruses in a collection of Halophytophthora spp. isolates from estuarine ecosystems in southern Portugal. Subsequently, stranded RNA sequencing followed by de novo contig assembling and comparative BLAST searches were used to identify eight putative novel viral sequences in two dsRNA-hosting isolates (Figure 1) belonging to the same Halophytophthora species (H. sp. 04, Supplementary Table S1). Primer-specific RT-PCR and 5′ and 3′ RACE were performed to confirm the presence and complete the sequences of each virus. All eight viral sequences contained elements of RdRp genes and bunyavirus motifs (Figure 2). Their eight virus genomes ranged in size from 6.8 to 9.3 kb (Table 3). BLAST analysis revealed that the eight viruses were significantly similar in sequence (ca. 30% identity; Table 2) to a number of unclassified bunya-like viruses isolated from different fungi, oomycetes, insects and crabs (Table 2). HRV8 differs from HRV1-7 and appears in a phylogenetic outgroup and is 58.07% similar to Wuhan insect virus 3 which is apparently hosted by Asellus sp. in China.
The overall pairwise nucleotide and aa identity of the eight viruses was 30.7 and 16.6%, respectively. HRV2 and 3 appear to one another apparently share identities of 57.05 and 47.80% at the nucleotide and aa levels, respectively. The remainder of the pairwise sequence comparisons (PASC) between the HRV isolates reveal lower identities of 40 and 20% at the nucleotide and aa levels (Supplementary Tables S5, S6). Analysis of the RdRp of HRV1-8 illustrate that they differ > 10% between one another and with their most similar matches in the GenBank. Since the primary classification criteria for genus and species used currently are based on pairwise sequence comparisons (PASC) and phylogenetic analyses, the viruses discovered in this investigation may constitute eight novel virus species, designated as Halophytophthora RNA Virus 1-8 (HRV1-8).
Bunyaviruses are enveloped viruses with a genome consisting of three ssRNA segments (called L, M, and S), the S RNA encodes the nucleocapsid protein, the M glycoproteins and the L segments encode the RNA polymerase. Each genome segment is coated by the viral nucleoproteins (NPs) and the polymerase (L protein) to form a functional ribonucleoprotein (RNP) complex, which is necessary for the RNA replication and gene transcription (reviewed in ). However, in our study we have only discovered the L segment. Our result does not categorically rule out the existence of M and S segments but their copy number might be very low compared to the polymerase fragment. The NGS performed in this investigation may not have analyzed sufficient reads to identify smaller HRV genomic components but this is unlikely because it has been established that 100 M reads is sufficient to identify all RNAs of interest and the rRNA depletion worked successfully. A quality check for the presence of rRNAs was performed by mapping the reads on rRNA animal and human databases and <5% of rRNA reads were detected. Similar to HRV1-8 a number of other bunya-like mycoviruses, including BcNSRV-1 (), MpNSRV1 (), and PpRV1 () apparently only possess an RdRp. However, it is more plausible that the genome description of these viruses (including HRV1-8) is incomplete. The putative NP and other non-structural (Ns) associated proteins are likely not conserved enough to be detected by homology, in contrast to what has been observed in other viruses including PrNSRV1 () and LeNSRV2 ().
The phylogenetic tree shown in Figure 2, which includes HRV1-8 and 42 (−) ssRNA viruses, illustrates that HRV1-7 cluster with several unclassified viruses found in a variety of fungal, oomycete, nematode, crab and insect hosts, and a virus detected by NGS approach from the stool of a rhesus monkey (Macaca mulatta) (Zhao et al., 2017). More specifically, HRV4-6 are significantly similar in sequence to PpRV1 from the Arctic and Antarctic moss pathogen Pythium polare (). HRV2, 3, and 7 are closer to BcNSRV-1 from the air-borne fungal pathogen Botrytis cinerea () and MpNSRV1 from the soil-borne pathogen Macrophomina phaseolina (), and HRV1 is grouped with Combu negative-strand RNA mycovirus isolate C3B from the soil fungus Mucor irregularis in Belém, a port municipality in the Brasilian Amazon (accession number MH990635, unpublished), and Beihai sesarmid crab virus 5 detected by NGS from a sesarmid crab mix from Beihai, China. Sesarmid crabs have an important ecological role in mangrove ecosystems because they consume large amounts of leaf litter (). HRV8 has some sequence similarity with classified bunyaviruses and, particularly, with Wuhan insect virus 3, hosted by an individual of Asellus sp. This genus of isopod crustaceans is known to feed primarily on decaying vegetation, microscopic algae and small invertebrates (reviewed in ). Since Halophytophthora spp. are the first colonizers of fallen mangrove leaves (; ) and estuarine grasses due to their ability to produce large amounts of chemotactic zoospores (reviewed in ) they likely share their habitat with insects, nematodes and crustaceans, where they undoubtedly interact. As an example, the class Oomycota includes marine holocarpic pathogens of nematodes, algae, crustaceans and molluscs (Thines and Choi, 2016), thus some of the viruses discovered by NGS of invertebrates may have originated from their gut mycoflora and/or parasites. Even though arthropods are thought to be the ancestral hosts of bunyaviruses () if viruses with similar genomes continue to be discovered in other invertebrates, protists and fungi (Shi et al., 2018) this theory may have to be redefined.
The phylogenetic tree shown in Figure 2 has similarities with previous phylogenetic studies that group novel bunya-like mycoviruses separated from most of the currently classified bunyavirus families (; ; ; ). And, it also supports the proposal of different mycobunyaviral families (): family Deltamycobunyaviridae with MpNSRV1, BcNSRV1 also includes PpRV1 and HRV1-7; two other families (proposed alfamycobunyaviridae and betamycobunyaviridae) including myco-phlebo-like viruses (; Velasco et al., 2019; ) and the family Gammamycobunyaviridae with CcNSRV2, RsNSV4 and Plasmopara associated mycobunyavirales-like RNA Virus 2 (; ). In addition, we propose a fifth family (Epsilonmycobunyaviridae) for HRV8 and Wuhan insect virus 3. As indicated previously (), this classification is not formal as it has not been accepted by ICTV but it sheds light on the current knowledge of novel bunya-like mycoviruses.
The genetic variability of the partial RdRp sequences of HRV1-8 showed that they have a relatively low nucleotide but high haplotype diversity since all the strains were different (Table 4). In this regard, HRV5 possessed the highest nucleotide variability (0.12). It is well established that bunyavirus replication is error-prone and results in genome modification. Such modifications can accumulate over the time either due either to random genetic drift or as genetic adaptations of the virus to a new environment and/or a new host (Schneider and Roossinck, 2001; ; ). Novel viral genotypes can be generated through mutation, recombination and reassortment. Viral reassortment seems to be a powerful mechanism underlying the evolution of the Bunyavirales order (). pointed out that most bunyaviruses described so far are actually reassortants of existing or extinguished viruses. This process leads to the generation of progeny viruses with novel genomic organizations as a result of gene shuffling between coinfecting closely related bunyaviruses (). The same might exist for HRV1-8, which are only hosted by four Halophytophthora isolates, BD641, BD647, BD650 and BD654. These isolates belong to the same species (H. sp. 04) and were sampled from the same site, Parque Natural da Ria Formosa, Santa Luzia, Tavira, Portugal. As this yet undescribed Halophytophthora species has a functional homothallic sexual system it is feasible that the viruses were transferred between isolates during the mating process or by simple contact between vegetative hyphae, resulting in the coinfection of their viruses. Interestingly, HRV1-8 only occur in H. sp. 04 despite the presence of at least seven Halophytophthora species at the Santa Luzia site. The viruses were not found in isolates of H. sp. 04 at the Ria de Alvor, Alvor, Portimão site, suggesting that HRV1-8 might only be transmitted intraspecifically between co-occurring Halophytophthora isolates. Little is known about virus transmission in oomycetes. The only studies concerning virus transmission in oomycetes were performed in P. infestans, where several viruses are stably maintained (). For instance, individual zoospores show 100% inheritance of PiRV-3 (). In addition, PiRV-2 is readily horizontally transmitted by hyphal anastomosis, and vertically transmitted by asexual reproduction through sporangia (). However, attempts to transfer PiRV-2 into apparently vegetative incompatible P. infestans isolates failed. Although several studies have demonstrated that interspecies transmission does occur naturally (Van Diepeningen et al., 1998; ; Vainio et al., 2011), mycovirus transmission typically occurs between strains of the same fungal species. Viruses spread readily through fungal hyphal networks crossing pores contained in compartmenting septa. However, as potential detrimental cytoplasmic elements, their transmission between strains may be restrained by a genetic self/non-self recognition system termed vegetative incompatibility (vic), mating type incompatibility or intersterility (). Oomycetes lack septa and in some cases, as in the genus Phytophthora, they are prompted to interspecific hybridizations which play a major role in speciation and species radiations in diverse natural ecosystems (; ). However, viruses might also be recognized as invasive compounds of the cytoplasmatic entities limiting their transfer to another species. While vic and mating systems serve as antiviral defense mechanisms at the population level, RNA silencing or RNA interference (RNAi) provides a fungal antiviral defense response at the cellular level (). Antiviral RNA silencing has been demonstrated in different types of fungi including the chestnut pathogen Cryphonectria parasitica (Segers et al., 2007) or the arbuscular mycorrhizal fungus Gigaspora margarita (Silvestri et al., 2020) but also in oomycetes (). Since the RNAi machinery targets possible detrimental non-self-nucleic acids, virus-infected host organisms are normally enriched with viral small interfering (si) RNA. For instance, in an analysis of the virus-derived small RNAs following high-throughput sequencing of the Halophytophthora isolate BD647 small RNA reads were mapped to HRV6 (unpublished data).
Conclusion
A combination of traditional and new technologies has been used to identify and sequence eight bunya-like mycoviruses that coinfect isolates belonging to the same species of Halophytophthora from southern Portuguese estuaries. However, any relationships between the different viruses and any effects of the viruses on the phenotype, virulence and host range of their oomycete host remain unknown.
Statements
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 at: https://www.ncbi.nlm.nih.gov/bioproject/PRJNA619952.
Author contributions
LB and TJ: original idea. CM, TJ, and MJ: sampling and phylogenetic characterization of Halophytophthora isolates. LB, JJ, and MR: methodology. LB: formal analysis, writing—original draft preparation. LB, TJ, and JJ: writing—review and editing.
Funding
This research was funded by the Project Phytophthora Research Centre Reg. No. CZ.02.1.01/0.0/0.0/15_003/0000453, co-financed by the European Regional Development Fund.
Acknowledgments
LB would like to acknowledge Dr. María Ángeles Ayllón from Centro de Biotecnología y Genómica de Plantas (CBGP-UPM-INIA), Universidad Politécnica de Madrid (UPM) (Spain) for her useful suggestions.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.01467/full#supplementary-material
TABLE S1Complete list of isolates used in the present study, indicating species identification, location, dsRNA pattern and results of the RT-PCR screening title.
TABLE S2Mapping results for P. cinnamomi as reference sequence.
TABLE S3Gene-Specific Primers (GSPs) used in the SMARTer RACE reactions.
TABLE S4Specific primers used for the partial amplification of the viral RdRp in order to screen the presence of each virus in all the isolates of the collection.
TABLE S5Percentages of pairwise identities between the eight viral full-length nucleotide genome sequences.
Footnotes
1.^ftp://ftp.ncbi.nlm.nih.gov/genomes/all/GCA/001/314/365/GCA_001314365.1_MP94-48v2/GCA_001314365.1_MP94-48v2_genomic.fna.gz
2.^http://bio-bwa.sourceforge.net
4.^http://code.google.com/p/bedtools/
5.^https://www.ebi.ac.uk/~zerbino/velvet/
6.^https://www.ebi.ac.uk/~zerbino/oases/
7.^http://www.htslib.org/doc/samtools-11.html
8.^http://tmcalculator.neb.com/
9.^https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi
References
1
AbudurexitiA.AdkinsS.AliotoD.AlkhovskyS. V.Avšič-ŽupancT.BallingerM. J.et al (2019). Taxonomy of the order bunyavirales: update 2019.Arch. Virol.1641949–1965.
2
BarrJ. N.FearnsR. (2010). How RNA viruses maintain their genome integrity.J. Gen. Virol.911373–1387. 10.1099/vir.0.020818-20810
3
BennettR. M.de CockA. W. A. M.LévesqueC. A.ThinesM. (2017). Calycofera gen. nov., an estuarine sister taxon to Phytopythium.Peronosporaceae. Mycol. Prog.16947–954. 10.1007/s11557-017-1326-1329
4
BennettR. M.ThinesM. (2019). Revisiting salisapiliaceae.FUSE3171–184. 10.3114/fuse.2019.03.10
5
BreitbartM. (2012). Marine viruses: truth or dare.Ann. Rev. Mar. Sci.4425–448. 10.1146/annurev-marine-120709-142805
6
BrieseT.CalisherC. H.HiggsS. (2013). Viruses of the family Bunyaviridae : are all available isolates reassortants?Virology446207–216. 10.1016/j.virol.2013.07.030
7
CaiG.FryW. E.HillmanB. I. (2019). PiRV-2 stimulates sporulation in Phytophthora infestans.Virus Res271:197674. 10.1016/J.VIRUSRES.2019.197674
8
CaiG.KrychiwJ. F.MyersK.FryW. E.HillmanB. I. (2013). A new virus from the plant pathogenic oomycete Phytophthora infestans with an 8 kb dsRNA genome: the sixth member of a proposed new virus genus.Virology435341–349.
9
CaiG.MyersK.FryW. E.HillmanB. I. (2012). A member of the virus family Narnaviridae from the plant pathogenic oomycete Phytophthora infestans.Arch. Virol.157165–169.
10
CaiG.MyersK.HillmanB. I.FryW. E. (2009). A novel virus of the late blight pathogen, Phytophthora infestans, with two RNA segments and a supergroup 1 RNA-dependent RNA polymerase.Virology39252–61. 10.1016/j.virol.2009.06.040
11
ChiapelloM.Rodríguez-RomeroJ.NervaL.ForgiaM.ChitarraW.AyllónM. A.et al (2020). Putative new plant viruses associated with Plasmopara viticola -infected grapevine samples.Ann. Appl. Biol.176180–191. 10.1111/aab.12563
12
CoupeauD.BayrouC.BaillieuxP.MarichalA.LenaertsA. C.CatyC.et al (2019). Host-dependence of in vitro reassortment dynamics among the sathuperi and shamonda simbuviruses.Emerg. Microbes Infect.8381–395. 10.1080/22221751.2019.1586410
13
CoyS.GannE.PoundH.ShortS.WilhelmS. (2018). Viruses of eukaryotic algae: diversity, methods for detection, and future directions.Viruses10:487. 10.3390/v10090487
14
DesnuesC.Rodriguez-BritoB.RayhawkS.KelleyS.TranT.HaynesM.et al (2008). Biodiversity and biogeography of phages in modern stromatolites and thrombolites.Nature452340–343. 10.1038/nature06735
15
DonaireL.PagánI.AyllónM. A. (2016). Characterization of Botrytis cinerea negative-stranded RNA virus 1, a new mycovirus related to plant viruses, and a reconstruction of host pattern evolution in negative-sense ssRNA viruses.Virology499212–218. 10.1016/J.VIROL.2016.09.017
16
EdgarR. C. (2004). MUSCLE: multiple sequence alignment with high accuracy and high throughput.Nucleic Acids Res.321792–1797. 10.1093/nar/gkh340
17
FahlgrenN.BollmannS. R.KasschauK. D.CuperusJ. T.PressC. M.SullivanC. M.et al (2013). Phytophthora have distinct endogenous small RNA populations that include short interfering and microRNAs.PLoS One8:e77181. 10.1371/journal.pone.0077181
18
FerronF.WeberF.de la TorreJ. C.RegueraJ. (2017). Transcription and replication mechanisms of Bunyaviridae and Arenaviridae L proteins.Virus Res.234118–134. 10.1016/j.virusres.2017.01.018
19
GillingsM. R.TesorieroL. A.GunnL. V. (1993). Detection of double-stranded RNA and virus-like particles in Australian isolates of Pythium irregulare.Plant Pathol.426–15. 10.1111/j.1365-3059.1993.tb01466.x
20
GoversL. L.Man In ‘T VeldW. A.MeffertJ. P.BoumaT. J.van RijswickP. C. J.HeusinkveldJ. H. T.et al (2016). Marine Phytophthora species can hamper conservation and restoration of vegetated coastal ecosystems.Proc. R. Soc. B Biol. Sci.283:20160812. 10.1098/rspb.2016.0812
21
GrasseW.SpringO. (2017). ssRNA viruses from biotrophic Oomycetes form a new phylogenetic group between Nodaviridae and Tombusviridae. Arch. Virol.1621319–1324. 10.1007/s00705-017-3243-3242
22
GrasseW.ZipperR.TotskaM.SpringO. (2013). Plasmopara halstedii virus causes hypovirulence in Plasmopara halst edii, the downy mildew pathogen of the sunflower.Fungal Genet. Biol.5742–47. 10.1016/j.fgb.2013.05.009
23
HackerC. V.BrasierC. M.BuckK. W. (2005). A double-stranded RNA from a Phytophthora species is related to the plant endornaviruses and contains a putative UDP glycosyltransferase gene.J. Gen. Virol.861561–1570. 10.1099/vir.0.80808-80800
24
HalldorssonS.LiS.LiM.HarlosK.BowdenT. A.HuiskonenJ. T. (2018). Shielding and activation of a viral membrane fusion protein.Nat. Commun.91–9. 10.1038/s41467-017-02789-2782
25
HuiskonenJ. T.HepojokiJ.LaurinmakiP.VaheriA.LankinenH.ButcherS. J.et al (2010). Electron cryotomography of tula hantavirus suggests a unique assembly paradigm for enveloped viruses.J. Virol.844889–4897. 10.1128/jvi.00057-10
26
JungT.ChangT. T.BakonyiJ.SeressD.Pérez-SierraA.YangX.et al (2017a). Diversity of Phytophthora species in natural ecosystems of Taiwan and association with disease symptoms.Plant Pathol.66194–211. 10.1111/ppa.12564
27
JungT.ScanuB.BakonyiJ.SeressD.KovácsG. M.DuránA.et al (2017b). Nothophytophthora gen. nov., a new sister genus of Phytophthora from natural and semi-natural ecosystems.Persoonia39143–174. 10.3767/persoonia.2017.39.07
28
KormelinkR.GarciaM. L.GoodinM.SasayaT.HaenniA.-L. (2011). Negative-strand RNA viruses: the plant-infecting counterparts.Virus Res.162184–202. 10.1016/j.virusres.2011.09.028
29
KozlakidisZ.BrownN. A.JamalA.PhoonX.CouttsR. H. A. (2010). Incidence of endornaviruses in Phytophthora taxon douglasfir and Phytophthora ramorum. Virus Genes40130–134. 10.1007/s11262-009-0421-427
30
LeeS. Y. (1998). Ecological role of grapsid crabs in mangrove ecosystems: a review.Mar. Freshw. Res.49335–343. 10.1071/mf97179
31
LeslieJ. F.ZellerK. A. (1996). Heterokaryon incompatibility in fungi—more than just another way to die.J. Genet.75415–424. 10.1007/BF02966319
32
LiC.-X.ShiM.TianJ.-H.LinX.-D.KangY.-J.ChenL.-J.et al (2015). Unprecedented genomic diversity of RNA viruses in arthropods reveals the ancestry of negative-sense RNA viruses.Elife4:e05378. 10.7554/eLife.05378
33
LiH.RoossinckM. J. (2004). Genetic bottlenecks reduce population variation in an experimental rna virus population.Society7810582–10587. 10.1128/JVI.78.19.10582
34
LinY.-H.FujitaM.ChibaS.HyodoK.AndikaI. B.SuzukiN.et al (2019). Two novel fungal negative-strand RNA viruses related to mymonaviruses and phenuiviruses in the shiitake mushroom (Lentinula edodes).Virology533125–136. 10.1016/j.virol.2019.05.008
35
LuS.WangJ.ChitsazF.DerbyshireM. K.GeerR. C.GonzalesN. R.et al (2020). CDD/SPARCLE: the conserved domain database in 2020.Nucleic Acids Res.48D265–D268. 10.1093/nar/gkz991
36
Man in ’t VeldW. A.RosendahlK. C. H. M.van RijswickP. C. J.MeffertJ. P.BoerE.WestenbergM.et al (2019). Multiple Halophytophthora spp. and Phytophthora spp. including P. gemini, P. inundata and P. chesapeakensis sp. nov. isolated from the seagrass Zostera marina in the Northern hemisphere.Eur. J. Plant Pathol.153341–357. 10.1007/s10658-018-1561-1561
37
MaranoA. V.JesusA. L.de SouzaJ. I.LeañoE. M.JamesT. Y.JerônimoG. H.et al (2014). A new combination in Phytopythium: P. kandeliae (Oomycetes. Straminipila).Mycosphere5510–522. 10.5943/mycosphere/5/4/3
38
MargariaP.BoscoL.VallinoM.CiuffoM.MautinoG. C.TavellaL.et al (2014). The NSs Protein of tomato spotted wilt virus is required for persistent infection and transmission by Frankliniella occidentalis.J. Virol.885788–5802. 10.1128/jvi.00079-14
39
MarklewitzM.ZirkelF.KurthA.DrostenC.JunglenaS. (2015). Evolutionary and phenotypic analysis of live virus isolates suggests arthropod origin of a pathogenic RNA virus family.Proc. Natl. Acad. Sci. U.S.A.1127536–7541. 10.1073/pnas.1502036112
40
MarzanoS.-Y. L.NelsonB. D.Ajayi-OyetundeO.BradleyC. A.HughesT. J.HartmanG. L.et al (2016). Identification of diverse mycoviruses through metatranscriptomics characterization of the viromes of five major fungal plant pathogens.J. Virol.906846–6863. 10.1128/jvi.00357-316
41
MelzerM. S.IkedaS. S.BolandG. J. (2002). Interspecific Transmission of Double-Stranded RNA and Hypovirulence from Sclerotinia sclerotiorum to S. minor.Phytopathology92780–784. 10.1094/PHYTO.2002.92.7.780
42
MillerM. A.PfeifferW.SchwartzT. (2010). “Creating the CIPRES Science Gateway for inference of large phylogenetic trees,” in 2010 Gateway Computing Environments Workshop, GCE 2010.10.1109/GCE.2010.5676129
43
MorrisT. J.DoddsJ. A. (1979). Isolation and analysis of double-stranded RNA from virus-infected plant and fungal tissue.Phytopathology69854. 10.1094/Phyto-69-854
44
MullerR.PochO.DelarueM.BishopD. H. L.BouloyM. (1994). Rift valley fever virus L segment: correction of the sequence and possible functional role of newly identified regions conserved in RNA-dependent polymerases.J. Gen. Virol75(Pt 6), 1345–1352. 10.1099/0022-1317-75-6-1345
45
NagasakiK.TomaruY.TakaoY.NishidaK.ShiraiY.SuzukiH.et al (2005). Previously unknown virus infects marine diatom.Appl. Environ. Microbiol713528–3535. 10.1128/AEM.71.7.3528-3535.2005
46
NervaL.CiuffoM.VallinoM.MargariaP.VareseG. C.GnaviG.et al (2015). Multiple approaches for the detection and characterization of viral and plasmid symbionts from a collection of marine fungi.Virus Res21922–38. 10.1016/j.virusres.2015.10.028
47
NervaL.ForgiaM.CiuffoM.ChitarraW.ChiapelloM.VallinoM.et al (2019a). The mycovirome of a fungal collection from the sea cucumber Holothuria polii.Virus Res.273:197737. 10.1016/j.virusres.2019.197737
48
NervaL.TurinaM.ZanzottoA.GardimanM.GaiottiF.GambinoG.et al (2019b). Isolation, molecular characterization and virome analysis of culturable wood fungal endophytes in esca symptomatic and asymptomatic grapevine plants.Environ. Microbiol.212886–2904. 10.1111/1462-2920.14651
49
NewellS. Y.FellJ. W. (1992). Distribution and experimental responses to substrate of marine oomycetes (Halophytophthora spp.) in mangrove ecosystems.Mycol. Res.96851–856. 10.1016/S0953-7562(09)81030-81037
50
NewellS. Y.FellJ. W. (1997). Competition among mangrove oomycotes, and between oomycotes and other microbes.Aquat. Microb. Ecol.1221–28. 10.3354/ame012021
51
NussD. L. (2011). Mycoviruses. RNA silencing, and viral RNA recombination.Adv. Virus Res.8025–48. 10.1016/B978-0-12-385987-7.00002-6
52
ObijeskiJ. F.BishopD. H.MurphyF. A.PalmerE. L. (1976). Structural proteins of La Crosse virus.J. Virol.19985–997. 10.1128/jvi.19.3.985-997.1976
53
O’CallaghanI.HarrisonS.FitzpatrickD.SullivanT. (2019). The freshwater isopod Asellus aquaticus as a model biomonitor of environmental pollution: a review.Chemosphere235498–509. 10.1016/j.chemosphere.2019.06.217
54
PoimalaA.VainioE. J. (2020). Complete genome sequence of a novel toti - like virus from the plant - pathogenic oomycete Phytophthora cactorum.Arch. Virol.165, 1679–1682. 10.1007/s00705-020-04642-4642
55
RozasJ.Ferrer-MataA.Sanchez-Del BarrioJ. C.Guirao-RicoS.LibradoP.Ramos-OnsinsS. E.et al (2017). DnaSP 6: DNA sequence polymorphism analysis of large data sets.Mol. Biol. Evol.343299–3302. 10.1093/molbev/msx248
56
SasaiS.TamuraK.TojoM.HerreroM.-L.HoshinoT.OhkiS. T.et al (2018). A novel non-segmented double-stranded RNA virus from an Arctic isolate of Pythium polare. Virology522234–243. 10.1016/j.virol.2018.07.012
57
SchardlC. L.CravenK. D. (2003). Interspecific hybridization in plant-associated fungi and oomycetes: a review.Mol. Ecol.122861–2873. 10.1046/j.1365-294X.2003.01965.x
58
SchneiderW. L.RoossinckM. J. (2001). Genetic diversity in RNA virus quasispecies is controlled by host-virus interactions.Society756566–6571. 10.1128/JVI.75.14.6566
59
SegersG. C.ZhangX.DengF.SunQ.NussD. L. (2007). Evidence that RNA silencing functions as an antiviral defense mechanism in fungi.Proc. Natl. Acad. Sci. U.S.A.10412902–12906. 10.1073/pnas.0702500104
60
ShiM.LinX. D.TianJ. H.ChenL. J.ChenX.LiC. X.et al (2016). Redefining the invertebrate RNA virosphere.Nature540539–543. 10.1038/nature20167
61
ShiM.ZhangY. Z.HolmesE. C. (2018). Meta-transcriptomics and the evolutionary biology of RNA viruses.Virus Res.24383–90. 10.1016/j.virusres.2017.10.016
62
ShibaK.HattaC.SasaiS.TojoM.OhkiS. T.MochizukiT. (2018). Genome sequence of a novel partitivirus identified from the oomycete Pythium nunn.Arch. Virol.1632561–2563.
63
SilvestriA.TurinaM.FiorilliV.MiozziL.VeniceF.BonfanteP.et al (2020). Different genetic sources contribute to the small RNA Population in the arbuscular mycorrhizal fungus Gigaspora margarita.Front. Microbiol.11:395. 10.3389/fmicb.2020.00395
64
StamatakisA.HooverP.RougemontJ. (2008). A rapid bootstrap algorithm for the RAxML Web Servers.Syst. Biol.57758–771. 10.1080/10635150802429642
65
SullivanB. K.Trevathan-TackettS. M.NeuhauserS.GoversL. L. (2018). Review: host-pathogen dynamics of seagrass diseases under future global change.Mar. Pollut. Bull.13475–88. 10.1016/j.marpolbul.2017.09.030
66
SutelaS.PoimalaA.VainioE. J. (2019). Viruses of fungi and oomycetes in the soil environment.FEMS Microbiol. Ecol.951–18. 10.1093/femsec/fiz119
67
SuttleC. A. (2005). Viruses in the sea.Nature437356–361. 10.1038/nature04160
68
ThinesM.ChoiY. J. (2016). Evolution, diversity, and taxonomy of the Peronosporaceae, with focus on the genus Peronospora.Phytopathology1066–18. 10.1094/PHYTO-05-15-0127-RVW
69
TuomivirtaT. T.UotilaA.HantulaJ. (2002). Two independent double-stranded RNA patterns occur in the Finnish Gremmeniella abietina var. abietina type A.For. Pathol.32197–205. 10.1046/j.1439-0329.2002.00285.x
70
VainioE. J.HakanpääJ.DaiY.-C.HansenE.KorhonenK.HantulaJ. (2011). Species of Heterobasidion host a diverse pool of partitiviruses with global distribution and interspecies transmission.Fungal Biol.1151234–1243. 10.1016/j.funbio.2011.08.008
71
Van DiepeningenA. D.DebetsA. J. M.HoekstraR. F. (1998). Intra- and interspecies virus transfer in Aspergilli via protoplast fusion.Fungal Genet. Biol.25171–180. 10.1006/fgbi.1998.1096
72
Vega ThurberR. L.BarottK. L.HallD.LiuH.Rodriguez-MuellerB.DesnuesC.et al (2008). Metagenomic analysis indicates that stressors induce production of herpes-like viruses in the coral Porites compressa.Proc. Natl. Acad. Sci. U.S.A.10518413–18418. 10.1073/pnas.0808985105
73
VelascoL.Arjona-GironaI.CretazzoE.López-HerreraC. (2019). Viromes in Xylariaceae fungi infecting avocado in Spain.Virology53211–21. 10.1016/j.virol.2019.03.021
74
WeilandJ. J.SundsbakJ. L. (2000). Differentiation and detection of sugar beet fungal pathogens using pcr amplification of actin coding sequences and its region of the rRNA gene.Plant Dis.84475–482. 10.1094/PDIS.2000.84.4.475
75
WeinbauerM. G. (2004). Ecology of prokaryotic viruses.FEMS Microbiol. Rev.28127–181. 10.1016/j.femsre.2003.08.001
76
WeynbergK.AllenM.WilsonW. (2017). Marine prasinoviruses and their tiny plankton hosts: a review.Viruses9:43. 10.3390/v9030043
77
WommackK. E.ColwellR. R. (2000). Virioplankton: viruses in aquatic ecosystems.Microbiol. Mol. Biol. Rev.6469–114. 10.1128/MMBR.64.1.69-114.2000
78
YokoiT.YamashitaS.HibiT. (2003). The nucleotide sequence and genome organization of Sclerophthora macrospora virus A.Virology311394–399.
79
ZhaoG.WuG.LimE. S.DroitL.KrishnamurthyS.BarouchD. H.et al (2017). VirusSeeker, a computational pipeline for virus discovery and virome composition analysis.Virology50321–30. 10.1016/j.virol.2017.01.005
Summary
Keywords
dsRNA, (-) ssRNA, RdRp, mycovirus, marine virology, estuaries
Citation
Botella L, Janoušek J, Maia C, Jung MH, Raco M and Jung T (2020) Marine Oomycetes of the Genus Halophytophthora Harbor Viruses Related to Bunyaviruses. Front. Microbiol. 11:1467. doi: 10.3389/fmicb.2020.01467
Received
04 April 2020
Accepted
04 June 2020
Published
15 July 2020
Volume
11 - 2020
Edited by
Ioly Kotta-Loizou, Imperial College London, United Kingdom
Reviewed by
Robert Henry Arnold Coutts, University of Hertfordshire, United Kingdom; Massimo Turina, National Research Council (CNR), Italy
Updates
Copyright
© 2020 Botella, Janoušek, Maia, Jung, Raco and Jung.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Leticia Botella, qqbotell@mendelu.cz
This article was submitted to Virology, a section of the journal Frontiers in Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.