Next Article in Journal
Establishment and Comparison of Pathogenicity and Related Neurotropism in Two Age Groups of Immune Competent Mice, C57BL/6J Using an Indian Isolate of Chikungunya Virus (CHIKV)
Previous Article in Journal
A DNA Prime Immuno-Potentiates a Modified Live Vaccine against the Porcine Reproductive and Respiratory Syndrome Virus but Does Not Improve Heterologous Protection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Novel Victorivirus from a Pakistani Isolate of Alternaria alternata Lacking a Typical Translational Stop/Restart Sequence Signature

1
Institute of Plant Science and Resources, Okayama University, Chuou 2-20-1, Kurashiki, Okayama 710-0046, Japan
2
Crop Diseases Research Institute, National Agricultural Research Centre, Park Road, Islamabad 45500, Pakistan
*
Authors to whom correspondence should be addressed.
Viruses 2019, 11(6), 577; https://doi.org/10.3390/v11060577
Submission received: 4 June 2019 / Revised: 19 June 2019 / Accepted: 21 June 2019 / Published: 25 June 2019
(This article belongs to the Section Viruses of Plants, Fungi and Protozoa)

Abstract

:
The family Totiviridae currently contains five genera Totivirus, Victorivirus, Leishmavirus, Trichomonasvirus, and Giardiavirus. Members in this family generally have a set of two-open reading frame (ORF) elements in their genome with the 5′-proximal ORF (ORF1) encoding a capsid protein (CP) and the 3′-proximal one (ORF2) for RNA-dependent RNA polymerase (RdRp). How the downstream open reading frames (ORFs) are expressed is genus-specific. All victoriviruses characterized thus far appear to use the stop/restart translation mechanism, allowing for the expression of two separate protein products from bicitronic genome-sized viral mRNA, while the totiviruses use a −1 ribosomal frame-shifting that leads to a fusion product of CP and RdRp. We report the biological and molecular characterization of a novel victorivirus termed Alternaria alternata victorivirus 1 (AalVV1) isolated from Alternaria alternata in Pakistan. The phylogenetic and molecular analyses showed AalVV1 to be distinct from previously reported victoriviruses. AalVV1 appears to have a sequence signature required for the −1 frame-shifting at the ORF1/2 junction region, rather than a stop/restart key mediator. By contrast, SDS–polyacrylamide gel electrophoresis and peptide mass fingerprinting analyses of purified virion preparations suggested the expression of two protein products, not a CP-RdRp fusion product. How these proteins are expressed is discussed in this study. Possible effects of infection by this virus were tested in two fungal species: A. alternata and RNA silencing proficient and deficient strains of Cryphonectria parasitica, a model filamentous fungus. AalVV1 showed symptomless infection in all of these fungal strains, even in the RNA silencing deficient C. parasitica strain.

1. Introduction

There are many double-stranded RNA (dsRNA) viruses with a two-open reading frame (ORF), nonsegmented genome organization, one of the simplest genome types, reported from plants, fungi, protozoa, and insects, as exemplified by totiviruses (members of the genus Totivirus) and members of other genera in the Totiviridae [1], amalgaviruses [2,3], phlegiviruses [4], fusagraviruses [5], megatotiviruses [6], and many unclassified viruses. The genome sizes of such viruses range from 3.5 (amalgaviruses) to ~12 kbp (megatotiviruses). With some exceptions such as amalgaviruses and phlegiviruses, these viruses typically encode capsid protein (CP) on the 5′-proximal ORF and RNA-dependent RNA polymerase (RdRp) on the 3′-proximal ORF. The 3′ proximal ORFs of these viruses are expressed by noncanonical translation mechanisms such as stop/restart translation [7,8] and ribosomal frame-shifting [9,10,11], though not substantiated, for many of these viruses. These viruses are expected have icosahedral T = 1 capsids as shown for totiviruses [12,13].
Members of the family Totiviridae, comprising five genera (Totivirus, Victorivirus, Leishmaniavirus, Trichomonasvirus, and Giardiavirus), are characterized by undivided dsRNA genomes with two ORFs encoding a CP and an RdRp or RdRp domain of a fusion protein. Among them, most members of which (belonging to genera Totivirus and Victorivirus) are isolated from fungi. Most totiviruses (members of the genus Totivirus) characterized thus far have been reported from yeast as exemplified by Saccharomyces cerevisiae virus L-A (ScV-L-A). However, recently totiviruses were reported from obligate filamentous phytopathogenic fungi [14] and plants [15] (see [16] for review). The characterized victoriviruses were from filamentous fungi but not from yeast. Members of the two genera are different from each other in several features such as the translation strategy and phylogenetic placement. Yeast totiviruses snatch the 5′ cap structure from host mRNAs and the snatching activity resides in the CP and is then mapped to His-154 that is located in the trench on the outer-surface or the cytoplasmic surface of the virion [17,18]. Note that SsV-L-A is the first dsRNA fungal virus whose 3-D structure model of T = 1 capsids with 120 homodimers was proposed [12,19,20]. The 5′ proximal ORF is translated according to the scanning model while the 3′-proximal ORF is expressed as a CP-RdRp fusion product by the −1 ribosomal frame-shifting mediated by slippery sequences and downstream pseudoknot or stem-loop structures. The rate of the frame-shifting ribosomes was estimated as 1.9% of ribosomes that have translated CP [11].
Like totiviruses, victoriviruses form icosahedral particles of ~40 nm in diameter with a T = 1 lattice composed of 60 CP dimers (120 molecules) [21]. Victoriviruses appear to have internal ribosomal entry sites (IRESs) in the 5′ untranslated region for the CP translation [22], while the RdRp is translated from bicistronic viral mRNA by the stop/reinitiation or coupled termination/reinitiation mechanism as a separate product from CP, not a fusion product, unlike totiviruses [7,8]. The reinitiation in victoriviruses is most likely to be mediated largely by the “AUGA” or “UAAUG” [7,8], in which the underlined triplet serves as the stop codon for the CP ORF while the italicized triplet is the initiation codon for the RdRp ORF. The UAAUG is also used by the well-studied prototypic hypovirus Cryphonectria hypovirus 1 (CHV1, a single-strand RNA virus) [23]. Li et al. identified two RNA sequence elements, the AUGA motif combined with an upstream pseudoknot RNA structure, as the necessary sequence elements for the stop/restart translation of victoriviruses [8]. Therefore, members of the two genera are different in the translation strategy for both the first and second ORFs.
Almost all victorivirus show symptomless infections. Exceptions include Helminthosporium victoriae virus 190S (HvV190S) that induces debilitation in the natural fungal host (Helminthosporium victoriae) and an RNA silencing mutant of Cryphonectria parasitica [24], a model filamentous fungus for mycovirus research [25]. Another victorivirus from Rosellinia necatrix (Rosellinia necatrix victorivirus 1, RnVV1) causes no overt phenotypic alterations in the natural host and the standard strain of C. parasitica [26]. However, RnVV1 induces a growth defect in an antiviral RNA silencing deficient strain in which the antiviral host defense is compromised [26]. Interestingly, RnVV1 is very susceptible to RNA silencing and eliminated by coinfection with another RNA mycoviruses or transgenic expression of the hairpin dsRNA, both of which can induce host antiviral RNA silencing [27].
Here we report the molecular and biological characterization of a victorivirus strain, which was isolated from a field strain, A-16, of Alternaria alternata, a phytopathogenic ascomycete, isolated from a wheat grain in Pakistan. The virus was tentatively termed Alternaria alternata victorivirus 1 (AalVV1). Surprisingly, the virus lacks a typical stop/restart signature sequence for the translation of RdRp, thought to be the rule for victoriviruses, but instead likely has a sequence feature often found for the −1 frame-shifting at the ORF1 and ORF2 junction. This paper represents the second report on a fungal virus from the Indian subcontinent.

2. Materials and Methods

2.1. Fungal Culture and Isolation Condition

The fungal strain designated as A-16 was isolated from a wheat grain at the National Agricultural Research Centre (NARC), Islamabad, Pakistan by a standard blotter method as described by the International Seed Testing Association) ISTA [28].
The strain was identified as A. alternata based on morphological characteristics of spores and on the sequence of the internal transcribed spacer (ITS) regions of the fungal ribosomal DNA (rDNA) amplified by using a primer pair ITS1 (TCCGTAGGTGAACCTGCGG) and ITS4 (TCCTCCGCTTATTGATATGC) [29]. The culture was maintained on potato dextrose agar (PDA, BD Difco Laboratories, Detroit, MI, USA) at 25 °C and for RNA extractions (small and large scale), mycelial plugs were inoculated in potato dextrose broth (PDB, BD Difco Laboratories) at 25 °C for two weeks. To extract dsRNA, mycelial plugs were grown in PDB under shaking (150 rpm) at 25 °C for seven days for mycelial mass collection [30]. A virus-free Japanese strain, Ally-12 of A. alternata f. sp. lycopersici [31], was provided by Dr. Masatoki Taga, Okayama University, while the standard strain EP155 of C. parasitica and its mutant ∆dcl2 [32] were a generous gift from Dr. Donald L. Nuss, University of Maryland.

2.2. RNA Preparation

Total nucleic acid was extracted from the mycelia of the fungal strain A-16 cultured in 20 mL PDB while the dsRNA fraction was isolated from mycelia grown on PDA overlaid with cellophane (PDA-cellophane) until it covered the plate [30]. Harvested mycelia were ground in mortar with liquid nitrogen followed by the addition of buffer (100 mM Tris-HCl pH 8.0, 4 mM EDTA, 200 mM NaCl and 2% SDS) and clarification with one round each of phenol-chloroform and chloroform in 2ml tubes. Total nucleic acid fractions were obtained by ethanol precipitation. For dsRNA isolation, the resulting extract was mixed with 16% ethanol, STE buffer (10 mM Tris-HCl pH 8.0, 1 mM EDTA and 150 mM NaCl), and 0.05 g of cellulose powder with 200–300 mesh (Advantech, Tokyo, Japan), followed by one hour of incubation with continuous rotation at room temperature. The sample was washed thrice with STE and ethanol (16%), vortexed and centrifuged between washes and the resulting dsRNA was eluted once with STE buffer from dried cellulose powder and precipitated with ethanol and sodium acetate. The sample was analyzed using 1% agarose gel electrophoresis. A fragment of approximately 5 kbp was observed and the dsRNA nature was confirmed after treatment with RNase free DNase I and S1 nuclease (Takara, Shiga, Japan).

2.3. Viral Genome Sequencing using RNA-Seq

An equal amount of dsRNA samples prepared from three virus infected strains of different fungi including A. alternata strain A-16 and unidentified fungi were pooled (namely “pool A1”) and 5.9 µg of dsRNA (91 ng/µL) was subjected to next generation sequencing. Viruses detected from the other two strains of the fungi will be reported elsewhere. The cDNA library construction and subsequent deep sequencing on the Illumina platform (HiSeq 2500, 50 bp single-end reads) were carried out by Macrogen Inc. (Tokyo, Japan). After sequencing, the adaptor sequences are trimmed and then the adaptor-trimmed reads (29,875,302 reads) were assembled de novo into 14,209 contigs (~488–6785 nt in length) using CLC Genomics Workbench (version 11, CLC Bio-Qiagen, Aarhus, Demark). These contigs were subjected to local BLAST searches against the viral reference sequence (RefSeq) dataset of the National Center for Biotechnology Information (NCBI). The complete nucleotide sequence of the virus genome in this article has been deposited with the EMBL/GenBank/DDBJ Data Library under Accession No. LC477336.

2.4. Database Search and Sequence Analysis

Potential ORFs were identified using the EnzymeX version 3.3.3 [33] Sequence similarities were calculated using the protein BLAST program (blastp) available from NCBI. The RNA secondary structure (pseudoknots) was predicted using the DotKnot program [34] and drawn with PseudoViewer3 [35].

2.5. Phylogenetic Analysis

Phylogenetic analysis was based on a maximum likelihood (ML) method as described previously [36] with minor modifications. Virus sequences were aligned with the MAFFT version 7 [37] and subsequently ambiguously aligned regions were eliminated using Gblocks 0.91b [38]. We then used PhyML 3.0 [39] to construct the ML phylogenetic trees with the automatic model selection by Smart Model Selection (SMS) [40]. The phylogenetic tree was visualized and refined with the FigTree graphical viewer interface version 1.3.1 [41].

2.6. RT-PCR for Confirmation of Contigs and RLM-RACE for Terminal Sequence Determination

Reverse transcription (RT)-polymerase chain reaction (PCR) was carried out to confirm that contig (A1-19) was derived from the genome of a victorivirus. Approximately 10 µg of purified dsRNA along with 10 pmol random hexamer (N6) was denatured by dimethyl sulfoxide (DMSO) and reverse transcribed at 37 °C for an hour in a 20-µL reaction mixture containing 50 mM Tris-HCl (pH 8.3), 50 mM KCl, 4 mM MgCl2, 10 mM dithiothreitol, 40 U of moloney murine leukemia virus reverse transcriptase (Thermo Fisher Scientific, Invitrogen, Waltham, MA, USA), and 20 U of RNase inhibitor (Toyobo, Osaka, Japan). One µL of the reaction mixture was subjected to PCR with sequence-specific primer sets, AJ-19F (1)/AJ-19R (1) and AJ-19F (2) and AJ-19R (2), respectively (Table S1).
For determination of the terminal sequences 5′ and 3′, 3′ RNA ligase-mediated amplification of cDNA ends (RLM-RACE) was done according to the procedure described by Suzuki et al. [42] with slight modifications. A 5′-phosphorylated, 3′-amino-linked oligodeoxy-nucleotide (LIG-Rev, 5′-PO4-GATCCAACTAGTTCTAGAGCGG-NH2-3′) was ligated to the 3′-end of the dsRNA using T4 RNA ligase (Takara) at 16 °C for 16 h. The ligated dsRNA was then denatured and used as the template in the first strand cDNA synthesis at 42 °C for one h in which an oligo-deoxynucleotide (LIG-FOR, 5′-CCGCTCTAGAACTAGTTGGATC-3′), complementary to the adapter, and Superscript II RNase H-Reverse Transcriptase (Thermo Fisher Scientific) were used. Resulting cDNA was used for PCR with sequence-specific primer (Table S1) along with primer LIG-FOR. PCR products were cloned and then sequenced using the Sanger sequencing method.

2.7. Purification of Virus Particles

Viral particles were purified as described by Shamsi et al. [43]. This preparation was used for molecular analysis, transfection, and electron microscopy. For electron microscopy (EM) observations, purified virus particles were negatively strained with an EM stain (EM stainer, an alternative for uranyl acetate, Nissin EM Co., Tokyo, Japan). Prepared specimens were observed in a Hitachi model H-7650 transmission electron microscope (Hitachi, Tokyo, Japan) and photographed as digital images. SDS–polyacrylamide gel electrophoresis (PAGE) analysis followed with Coomassie Brilliant Blue (CBB) and silver staining was also performed as previously described [43].
The protein bands were silver-stained and subjected to in-gel tryptic digestion followed by liquid chromatography–tandem mass spectrometry (LC-MS/MS) analysis as described by Shamsi et al. [43].

2.8. Transfection of C. parasitica Spheroplasts

Transfection of C. parasitica was performed as described previously [44,45]. Spheroplasts of the C. parasitica Δdcl2 strain were prepared using the method of Eusebio-Cope et al. [30]. A-16 (a victorivirus) virions were purified according to the method described above by differential and cesium chloride (CsCl) gradient centrifugation. Purified particles were transfected into C. parasitica Δdcl2 spheroplasts using polyethylene glycol (PEG). After regeneration of the spheroplasts, mycelial plugs from multiple positions were transferred onto new PDA plates and propagated. Transfectants were analyzed by agarose gel electrophoretic analysis of their dsRNA-enriched fractions.

2.9. Protoplast Fusion

Protoplasts of the Japanese A. alternata strain Ally-12 were prepared with the method as described by Shahi et al. [46] and transformed by pCPXHY3, carrying a hygromycin resistance gene (hygromycin B phosphotransferase), to obtain a hygromycin-resistant (HygR) strain. The resulting transformant was used as a recipient, while the Pakistani strain A-16 (carrying a victorivirus) was used as a donor. Equal amounts of protoplast (2 × 107) from strains Ally-12 and A-16 were mixed and fused in the presence of PEG. Fungal colonies showing a HygR phenotype were tested for virus infection.

3. Results

3.1. Genome Organization of a Victorivirus Alternaria alternata victorivirus 1

A Pakistani fungal strain A-16 (Figure 1A) which carries a mycovirus-like dsRNA element (Figure 1B) was identified as Alternaria alternata morphologically and molecularly. The nucleotide sequences of the internal transcribed spacer (ITS) regions of A-16 were compared using [47], and were found to be matched best with the already deposited sequences of A. alternata (KM233267, identity: 94.6%).
The raw next generation sequencing (NGS) reads obtained by the HiSeq sequencing were assembled into different contigs. A local BLAST analysis showed the presence of a virus-like contig (A1-19), with an average coverage of 2680 (total 138,061 reads), and a size of 5159 bp in length (Figure 1C), which corresponded well to the length of the dsRNA purified from A-16. A NCBI BLASTn search of the 5159-bp contig revealed significant similarity to the complete genome sequence of Nigrospora oryzae victorivirus 1 (NoVV1, accession number KT428155) (E-value, 0; nt identity 68.3%) and Phomopsis vexans RNA virus (PvRV, accession number KP090346) (E-value, 0; nt identity 68.3%), which are unclassified dsRNA mycoviruses that belong to the genus Victorivirus. The nucleotide sequence of this putative mycovirus, which we called “Alternaria alternata victorivirus 1” (AalVV1) has been deposited in the GenBank/EMBL/DDBJ databases under accession number LC477336.
RNA ligase-mediated amplification of cDNA ends (RLM-RACE) and RT-PCR analyses revealed that the AalVV1 genomic dsRNA was 5120 bp long and contained two overlapping open reading frames (ORFs) on its plus-strand (Figure 1D). ORF1 was 2298 nt long (positions 267 to 2564 nt) and was predicted to encode a protein with a calculated molecular mass of 80.2 kDa. ORF2 was 2157 nt long (positions 2887 to 5043 nt) and could encode a 78.5-kDa protein (Figure 1D). The 5′- and 3′- untranslated regions (UTR) were found to be 266 nt and 77 nt long, respectively. A tandem stem loops was predicted upstream of the AalVV1 ORF2 (Figure 1E). Although there is no clear evidence for conserved sequences at the genome termini of different victoriviruses [48], the 5′ UTR of AalVV1 showed conservation in a few nucleotides at the beginning of the 5′ UTR with NoVV1 and PvRV respectively, whereas in the 3′-UTR, a small stretch of the conservation with NoVV1 was also observed. In the 5′-UTR, the octanucleotide sequence 5′-AGGGUUCC-3′, which is conserved in several victoriviruses, including Coniothyrium minitans RNA virus (CmRV, accession number AF527633) and Epichloe festucae virus 1 (EfV1, accession number AM261427), was found at nt 204–211 in the AlaVV1, 55 nt upstream of the AUG initiation codon of ORF1. Interestingly, instead of the octanucleotide, an extra seven conserved nucleotides (5′-AGGGUUCCGUUGAUC-3′, bold letters) were noted in AalVV1 and several other victoriviruses (Figure S1). One of the victoriviruses (CmRV) has this fifteen nucleotide sequence at the start of the 5′-UTR. These sequences might have an important role to play as they are found in the majority of the victoriviruses (Figure S1A, boxed). The victorivirus CPs are unique in the family Totiviridae in having an Ala/Gly/Pro-rich region predicted near their C-termini [48]. The C-terminal sequence of AalVV1 showed the similar Ala/Gly/Pro-rich region as shown in the prototype victorivirus Helminthosporium victoriae virus 190S (HvV190S) (Figure S1C).
Victoriviruses generally utilize stop/restart strategies to translate the RdRp-encoding ORF2 whose in which AUGA or UAAUG facilitate it together with pseudoknot structures located at immediately upstream of the facilitators [7,8]. Visual examination failed to detect cotranslational reinitiation facilitators in AalVV1, while a H-type pseudoknot structure was predicted upstream of the CP stop codon by the DotKnot program (Figure 1E). Instead, a possible -1 frameshift slippery site might present at the positions 2545 to 2562 followed by a second pseudoknot structure (Figure 1E). However, no typical heptanucleotide slippery sequence “X XXY YYZ” (spaced triplets represent preframeshift codons), which is a key signal element for -1 ribosomal frame-shifting [49], was found in the region (Figure 1E). However, a slippery-like sequence, “C CaA AAU,” was detected at the nt positions 2555–2561. RT-PCR fragments spanning the ORF1 and ORF2 junction, where the NGS read coverage was significantly low (Figure 1C, arrow), were cloned and sequenced. All six of the RT-PCR clones sequenced conformed to the genetic organization discussed above.
BLASTp searches of the deduced amino acid (aa) sequences of ORF1 of the contig A1-19 dsRNA revealed the highest identity, 69.8% (E-value, 0.0), to the putative CP of PvRV (Table 1). ORF2 had the highest identity, 61.4% (E-value, 0.0), to the putative RdRp of NoVV1 (Table 1). The aa sequence analysis of the putative RdRp of the AalVV1 revealed the presence of eight motifs that were conserved in the RdRps of the members of the family Totiviridae (Figure S1B). Accordingly, phylogenetic analysis based on the aligned motif sequences in RdRp showed that AalVV1 clustered together with victoriviruses (genus Victorivirus) and was separated from the members of the genera Leishmaniavirus, Trichomonasvirus and Totivirus in the family Totiviridae, and it showed the closet homology with PvRV and NoVV1 amongst all the known victoriviruses (Figure 2A). Similarly, aligned CP sequences also showed AalVV1 in a separate cluster with all other victoriviruses (Figure 2B). A comprehensive up-to-date table has been arranged with all the information of the victoriviruses reported so far indicating the name of the virus, size of dsRNA, encoded protein regions, host, GenBank accession numbers, and the presence of conserved sequences in the 5′-UTR (Table S2).
Based on the genomic organization, global amino acid sequence similarities of CP and RdRp, and the phylogeny, AalVV1 appears to represent the genome of a novel victorivirus, based on the demarcation criteria for mycoviruses in this genus. This study represents the first report of the full-length genome sequence of a victorivirus infecting A. alternata.

3.2. Morphology and Protein Components of Purified AalVV1 Particles

To obtain purified AalVV1 particles, crude extracts of strain A-16 were subjected to CsCl density-gradient centrifugation. The virus-free strain Ally-12 was also subjected to the experiment as a negative control. In the virus-particle (VP) fraction where the AalVV1 dsRNA was enriched, nonenveloped spherical particles of ~40 nm diameter were observed with transmission electron microscope (Figure 3A). Protein components of the VP fraction were detected by SDS-PAGE followed with CBB and silver staining. In the VP fraction of the A-16 strain, a major band corresponding to approximately 80 kDa in size was detected by both staining methods (Figure 3B). In addition, some minor bands sized larger than 250 kDa and smaller than 80 kDa were detected by silver staining of the VP fraction of A-16 (Figure 3B). These bands were undetectable for the samples of a virus-free strain (Ally-12) (Figure 3B). No protein bands in an estimated size of a potential CP-RdRp fusion protein (171 kDa) were observed (Figure 3B).
Sequences of peptide fragments obtained after in-gel tryptic digestion of the major 80 kDa protein band of A-16 was predicted by liquid chromatography/tandem mass spectrometry (LC-MS/MS). As a result, the band contained the peptides derived from both ORF1 and ORF2 of the AalVV1 genome (Figure 4 and Figure S2). Furthermore, the band also included the peptides translated from frame 1, the same frame as ORF2, of the region between two ORFs (Figure 4A and Figure S2). Note that such peptide sequences are underlined in Figure 4B. A large portion of peptide fragments detected by peptide mass fingerprinting (PMF) were derived from CP-encoding ORF1, while a small but significant portion, approximately one fourth of the former, was derived from RdRp-encoding ORF2 (Figure S2). Considering the estimated protein size of each ORF (Figure 1D), the major 80 kDa band likely contained at least two proteins translated from each ORF. It remains to be elucidated how RdRp and the interval peptides downstream of CP are translated. In any case, purified AalVV1 particles were composed of a major CP and RdRp that were encoded in ORF1 and ORF2, respectively.

3.3. Asymptomatic Infection of A. alternata by AalVV1

In order to determine the effect of the virus on the host, attempts were made to create virus-free isogenic line and for that purpose, treatment of a translation inhibition drug cycloheximide, which has been previously used to cure fungus of a mycovirus infection [50], along with the single spore isolation technique was done. A. alternata Pakistani A-16 strain was grown on PDA plates amended with different concentrations of cycloheximide (10, 50 and 100 mM) and then single spore isolation was done from the mycelial plugs. However, our attempts to cure A-16 were unsuccessful. Thus, we took a protoplast fusion approach which was used for different viruses.
First, a virus-free Japanese strain (Ally-12) of A. alternata was transformed with a HygR gene cassette as described by Shahi et al. [46] and used as a primary recipient. After protoplast fusion and hygromycin selection, we conducted several rounds of anastomosis using the untransformed Ally-12 strain (hygromycin B-susceptible, HygS) as a secondary recipient. The virus AalVV1 was finally moved to the original Ally-12 strain (HygS) and compared with the isogenic virus-free Ally-12 strain. Consequently, the two strains showed indistinguishable colony phenotype on PDA, indicating no symptoms induced by AalVV1 (Figure 5). Note that the Japanese virus-free strain Ally-12 showed indistinguishable colony morphology to strain A-16 (see Figure 1A).

3.4. AalVV1 is Infectious to the Fungus C. parasitica as Particles and Shows Symptomless Infection

Effects of AalVV1 infection on the fungal phenotype was also tested in a model fungus, C. parasitica, as an experimental host, which proved useful for studying virus/host and virus/virus interactions [25]. Spheroplasts of the standard virus-free C. parasitica strain EP155 and an RNA silencing-deficient strain, ∆dcl2, derived from EP155 were transfected by semipurified AalVV1 particles. Infected colonies were readily obtained in both the EP155 and ∆dcl2 background and showed an unaltered phenotype which was indistinguishable from virus-free EP155 or ∆dcl2 (Figure S3).
Agarose gel electrophoresis of total RNA fractions from transfectants showed AalVV1 accumulation in ∆dcl2, comparable to that in A. alternata (Figure 6A). A considerable low level of AalVV1 accumulation was observed in wild type EP155. AalVV1 was detectable only by RT-PCR (Figure 6B, see the 1 kbp-specific band).

4. Discussion

A number of RNA viruses have been reported from Alternaria spp., including dsRNA viruses: an alternavirus [51], a partitivirus [52], botybirnaviruses [43,53,54], chrysoviruses [55], and victoriviruses (a different virus than AalVV1 from Alternaria arborescens) [56]; positive-sense (+) single-stranded RNA (ssRNA) viruses: mitoviruses [57,58], a fusarivirus [59], and an endornavirus [60]; and other unidentified viruses [61]. This study describes the characterization of a victorivirus, AalVV1, from A. alternata, and contributes to the growing database of viruses of Alternaria spp. To our knowledge, this study represents the second report on a fungal virus from the Indian subcontinent next to the paper by Shamsi et al. [43] that reports a thorough characterization of a botybirnavirus (dsRNA virus) of A. alternata.
AalVV1 has interesting features similar and dissimilar to previously characterized victoriviruses. For the RdRp translation from bicistronic genome-sized viral mRNA, members of the family Totiviridae use −1 frame-shifting (totiviruses, trichomonasviruses, giardiaviruses), −2 frame-shifting (trichomonasvirses), presumable +1 frame-shifting or ribosomal hopping (leishmaniaviruses) and stop/restart strategies (victoriviruses) [10,62] (Figure 2A). The stop/restart translational regulation of victoriviruses was thoroughly explored for HvV190S, the prototype victorivirus [7,8] in which the 38-nt region including the AUGA reinitiation facilitator and an H-type pseudoknot structure immediately upstream of the tetra-nucleotide is necessary and sufficient for the stop/restart translation. Similar translational strategies are used in other fungal and animal viruses with (+)ssRNA genomes, such as hypoviruses [23] and caliciviruses [63,64], where 4–20% of ribosomes terminating translation at the stop codon of the upstream ORF reinitiate translation of RdRp. This translational efficiency is similar to that observed for −1 frame-shifting of the prototype totivirus, ScV-L-A [11]. The slippery sites and pseudoknot structures up- and down-stream of the termination codon facilitate −1 frame-shifting.
This study showed AalVV1 to have the sequence features frequently found in −1 frame-shifting and stop/restart mechanisms. AalVV1 has an H-type pseudoknot structure upstream of the ORF1 termination codon (Figure 1) as in the case for HvV190S, and yet lacks the stop/restart facilitator sequence such as AUGA or UAAUG. Simultaneously, the virus has a signature-like sequence “C CaA AAU” for possible −1 frame-shifting, slippery sequence conforming to the consensus sequence (X XXY YYZ) at map positions 2555–2561 immediately upstream of the stop codon, UAG, of ORF1, and a predicted stem-loop structure downstream of the stop codon (Figure 1D). Interestingly, the aforementioned H-type pseudoknot structure is retained in AalVV1 immediately upstream of the −1 frame-shifting slippery site (Figure 1D). This study provided some clues, but failed to determine which mechanism is used for the translation of the AalVV1 RdRp ORF. The SDS-PAGE analysis failed to detect the CP-RdRp fusion product expected to be expressed via −1 ribosomal frame-shifting at an expected migration position of 171 kDa. Note that the migration position of a minor band corresponding to >250 kDa is much slower than that expected for the fusion product (171 kDa) (Figure 3B), indicating that the >250 kDa band did not represent the CP-RdRp fusion protein. The PMF indicated that the RdRp peptide fragments were detected in polypeptides of approximately 80 kDa, and that the peptides included amino acid sequences encoded by the ORF1/2 junction region (map positions 2565 to 2886). Namely, it was shown that translation of the −1 frame of ORF2 started upstream of the AUG codon at 2887. These SDS-PAGE and PMF data favor the stop/restart mechanism, but AalVV1 lacks a typical stop/restart facilitator such as AUGA and UAAUG around the ORF1/2 junction. If the virus utilizes the −1 frameshift, proteolytic processing or degradation must occur to account for the size (80 kDa) of the detected proteins. If the virus uses termination-coupled reinitiation, there must be an unidentified noncanonical start codon downstream the H-type pseudoknot structure and upstream of the nucleotide sequences (map positions 2593 to 2622) encoding amino acids TASDLDLYLK detected by peptide mass fingerprinting (Figure 4B and Figure S2).
Both the phylogenetic trees with the CP and RdRp sequences showed two unclassified victoriviruses, Nigrospora oryzae victorivirus 1 (NoVV1) and Phomopsis vexans RNA virus (PvRV), to make up a sister clade to the one containing AalVV1 (Figure 2A,B). NoVV1 and PvRV were isolated from fungi phylogenetically distant from A. alternata, the original host of AalVV1. This clade does not include another Alternaria victorivirus (Alternaria arborescens victorivirus 1) from A. arborescens, the tomato pathotype of A. alternata [56]. This is not surprising given that victoriviruses have broad host ranges as in the case of partitviruses [24,26,27]. There might have been horizontal transmission as suggested for other fungal viruses [65,66,67,68]. An additional interesting finding is a possible difference in the translational strategy for ORF2. NoVV1 appears to use the stop/restart mechanism like many other victoriviruses. However, PvRV interestingly lacks a reinitiation facilitator like AalVV1, and ORF2 is situated in the −2 or +1 frame relative to ORF1 (Figure S4).
Victoriviruses are ubiquitous in filamentous fungi [48]. However, only a few victoriviruses have been explored at the molecular level. As aforementioned, the prototypic HvV190S has been utilized for the investigation of viral gene expression [7,8,69]. Biological characterization has been limited to a few victoriviruses largely because inoculation methods had been unavailable until recently. Chiba et al. showed RnVV1 by virion transfection to induce asymptomatic infection in the original fungal host, R. necatrix and a model filamentous fungus, C. parasitica [26]. Subsequently, Xie et al. succeeded in the transfection of the same fungus with HvV190S as a new experimental host [24]. To our knowledge, HvV190S is the only victorivirus that was shown to induce debilitation and hypovirulence on their original fungal host (H. victoriae) [24]. In this study, AalVV1 from a Pakistani A. alternata strain (A-16) was shown to cause symptomless infection in wild type strains of a Japanese A. alternata (Ally-12) and the RNA silencing-proficient standard strain of C. parasitica (EP155). Moreover, unlike RnVV1 and HvV190S [24,26], AalVV1 induced no phenotypic alterations even in an antiviral RNA silencing deficient strain of C. parastica, ∆dcl2, in spite of the enhanced accumulation of AalVV1 in ∆dcl2 (Figure 6 and Figure S3), indicating that antiviral RNA silencing operating in C. parasitica targets this virus like other viruses [27,32,70,71]. This is a rare phenomenon: while many fungal viruses transfected cause minor symptoms or little in the wild type RNA silencing proficient strain EP155, ∆dcl2 infected by those viruses show enhanced symptom infection as exemplified by three dsRNA mycoviruses, RnVV1, Rosellinia necatrix partitivirus 2 (RnPV2), and RnPV6 [26,45,72]. A case similar to AalVV1 is a hypovirus (Cryphonectria hypovirus 4) that causes no phenotypic changes in ∆dcl2 or EP155 [73]. The occasional loss (data not shown) and considerably low accumulation of AalVV1 in EP155 (Figure 6) suggest its greater susceptibility to antiviral RNA silencing. Despite these attributes, this study clearly showed the infectivity of AalVV1 to C. parasitica, belonging to a class, Sordariomycetes, different from another one, Dothideomycetes, which include the original host, A. alternata. The previous and current studies suggest broad host ranges of victoriviruses in contrast to the relatively narrow host ranges of two other C. parasitica viruses, a chrysovirus, Cryphonectria nitschkei chrysovirus 1 [74] and a mitovirus, Cryphonectria parasitica mitovirus 1 [46].

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4915/11/6/577/s1, Figure S1: Multiple alignments of victorivirus nucleotide and amino acid sequences, Figure S2: The Mascot search result of the LC-MS/MS analysis of the major 80 kDa protein band of AalVV1, Figure S3: Colony morphology of Cryphonectria parasitica infected AalVV1, Figure S4: Predicted secondary structure of the ORF1/ORF2 junction sequences of the two AalVV1 related victoriviruses, Table S1: List of primers used in the confirmation and completion of the AalVV1 sequence, Table S2: Some properties of reported victoriviruses.

Author Contributions

A.J. and N.S. designed the experiments; A.J., Y.S., S.S., and W.S. performed the experiments; Y.S. and H.K. analyzed the data; N.S. wrote the paper; A.J., Y.S., and H.K. revised the manuscript.

Funding

This study was supported in part by Yomogi Inc. and Grants-in-Aid for Scientific Research on Innovative Areas from the Japanese Ministry of Education, Culture, Sports, Science and Technology (KAKENHI 25252011 and 16H06436, 16H06429 and 16K21723 to N.S. and H.K.).

Acknowledgments

The authors are grateful to Donald L. Nuss and Masatoki Taga for generous gifts of C. parasitica strains EP155 and ∆dcl2 and A. alternata strain Ally-12 and Ida Bagus Andika and Sakae Hisano for fruitful discussion. The authors are also grateful to Tsugumi Shiokawa and Hiroko Tada at Division of Instrumental Analysis, Okayama University for the protein analyses.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Wickner, R.B.; Ghabrial, S.A.; Nibert, M.L.; Patterson, J.L.; Wang, C.C. Family Totiviridae. In Virus Taxonomy: Ninth Report of the International Committee for the Taxonomy of Viruses; King, A.M.Q., Adams, M.J., Carstens, E.B., Lefkowits, E.J., Eds.; Elsevier, Academic Press: New York, NY, USA, 2011; pp. 639–650. [Google Scholar]
  2. Sabanadzovic, S.; Valverde, R.A.; Brown, J.K.; Martin, R.R.; Tzanetakis, I.E. Southern tomato virus: The link between the families Totiviridae and Partitiviridae. Virus Res. 2009, 140, 130–137. [Google Scholar] [CrossRef] [PubMed]
  3. Martin, R.R.; Zhou, J.; Tzanetakis, I.E. Blueberry latent virus: An amalgam of the Partitiviridae and Totiviridae. Virus Res. 2011, 155, 175–180. [Google Scholar] [CrossRef] [PubMed]
  4. Kozlakidis, Z.; Hacker, C.V.; Bradley, D.; Jamal, A.; Phoon, X.; Webber, J.; Brasier, C.M.; Buck, K.W.; Coutts, R.H. Molecular characterisation of two novel double-stranded RNA elements from Phlebiopsis gigantea. Virus Genes 2009, 39, 132–136. [Google Scholar] [CrossRef] [PubMed]
  5. Liu, H.; Fu, Y.; Xie, J.; Cheng, J.; Ghabrial, S.A.; Li, G.; Peng, Y.; Yi, X.; Jiang, D. Evolutionary genomics of mycovirus-related dsRNA viruses reveals cross-family horizontal gene transfer and evolution of diverse viral lineages. BMC Evol. Biol. 2012, 12, 91. [Google Scholar] [CrossRef] [PubMed]
  6. Arjona-Lopez, J.M.; Telengech, P.; Jamal, A.; Hisano, S.; Kondo, H.; Yelin, M.D.; Arjona-Girona, M.I.; Kanematsu, S.; Lopez-Herrera, C.; Suzuki, N. Novel, diverse RNA viruses from Mediterranean isolates of the phytopathogenic fungus, Rosellinia necatrix: Insights into evolutionary biology of fungal viruses. Env. Microbiol 2018, 20, 1464–1483. [Google Scholar] [CrossRef] [PubMed]
  7. Li, H.; Havens, W.M.; Nibert, M.L.; Ghabrial, S.A. RNA sequence determinants of a coupled termination-reinitiation strategy for downstream open reading frame translation in Helminthosporium victoriae virus 190S and other victoriviruses (family Totiviridae). J. Virol. 2011, 85, 7343–7352. [Google Scholar] [CrossRef] [PubMed]
  8. Li, H.; Havens, W.M.; Nibert, M.L.; Ghabrial, S.A. An RNA cassette from Helminthosporium victoriae virus 190S necessary and sufficient for stop/restart translation. Virology 2015, 474, 131–143. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Nibert, M.L.; Pyle, J.D.; Firth, A.E. A+1 ribosomal frameshifting motif prevalent among plant amalgaviruses. Virology 2016, 498, 201–208. [Google Scholar] [CrossRef]
  10. Atkins, J.F.; Loughran, G.; Bhatt, P.R.; Firth, A.E.; Baranov, P.V. Ribosomal frameshifting and transcriptional slippage: From genetic steganography and cryptography to adventitious use. Nucleic Acids Res. 2016, 44, 7007–7078. [Google Scholar] [CrossRef] [Green Version]
  11. Dinman, J.D.; Icho, T.; Wickner, R.B. A-1 ribosomal frameshift in a double-stranded RNA virus of yeast forms a gag-pol fusion protein. Proc. Natl. Acad. Sci. USA 1991, 88, 174–178. [Google Scholar] [CrossRef]
  12. Caston, J.R.; Trus, B.L.; Booy, F.P.; Wickner, R.B.; Wall, J.S.; Steven, A.C. Structure of L-A virus: A specialized compartment for the transcription and replication of double-stranded RNA. J. Cell Biol. 1997, 138, 975–985. [Google Scholar] [CrossRef] [PubMed]
  13. Luque, D.; Mata, C.P.; Suzuki, N.; Ghabrial, S.A.; Caston, J.R. Capsid structure of dsRNA fungal viruses. Viruses 2018, 10, 481. [Google Scholar] [CrossRef] [PubMed]
  14. Kondo, H.; Hisano, S.; Chiba, S.; Maruyama, K.; Andika, I.B.; Toyoda, K.; Fujimori, F.; Suzuki, N. Sequence and phylogenetic analyses of novel totivirus-like double-stranded RNAs from field-collected powdery mildew fungi. Virus Res. 2016, 213, 353–364. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Chen, S.; Cao, L.; Huang, Q.; Qian, Y.; Zhou, X. The complete genome sequence of a novel maize-associated totivirus. Arch. Virol. 2016, 161, 487–490. [Google Scholar] [CrossRef] [PubMed]
  16. Roossinck, M.J. Evolutionary and ecological links between plant and fungal viruses. New Phytol. 2019, 221, 86–92. [Google Scholar] [CrossRef] [PubMed]
  17. Fujimura, T.; Esteban, R. Cap-snatching mechanism in yeast L-A double-stranded RNA virus. Proc. Natl. Acad. Sci. USA 2011, 108, 17667–17671. [Google Scholar] [CrossRef] [Green Version]
  18. Wickner, R.B.; Fujimura, T.; Esteban, R. Viruses and prions of Saccharomyces cerevisiae. Adv. Virus Res. 2013, 86, 1–36. [Google Scholar]
  19. Cheng, R.H.; Caston, J.R.; Wang, G.J.; Gu, F.; Smith, T.J.; Baker, T.S.; Bozarth, R.F.; Trus, B.L.; Cheng, N.; Wickner, R.B.; et al. Fungal virus capsids, cytoplasmic compartments for the replication of double-stranded RNA, formed as icosahedral shells of asymmetric Gag dimers. J. Mol. Biol. 1994, 244, 255–258. [Google Scholar] [CrossRef]
  20. Naitow, H.; Tang, J.; Canady, M.; Wickner, R.B.; Johnson, J.E. L-A virus at 3.4 A resolution reveals particle architecture and mRNA decapping mechanism. Nat. Struct. Biol. 2002, 9, 725–728. [Google Scholar] [CrossRef]
  21. Caston, J.R.; Luque, D.; Trus, B.L.; Rivas, G.; Alfonso, C.; Gonzalez, J.M.; Carrascosa, J.L.; Annamalai, P.; Ghabrial, S.A. Three-dimensional structure and stoichiometry of Helmintosporium victoriae190S totivirus. Virology 2006, 347, 323–332. [Google Scholar] [CrossRef]
  22. Chiba, S.; Jamal, A.; Suzuki, N. First evidence for internal ribosomal entry sites in diverse fungal virus genomes. MBio 2018, 9. [Google Scholar] [CrossRef] [PubMed]
  23. Guo, L.H.; Sun, L.; Chiba, S.; Araki, H.; Suzuki, N. Coupled termination/reinitiation for translation of the downstream open reading frame B of the prototypic hypovirus CHV1-EP713. Nucleic Acids Res. 2009, 37, 3645–3659. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Xie, J.; Havens, W.M.; Lin, Y.H.; Suzuki, N.; Ghabrial, S.A. The victorivirus Helminthosporium victoriae virus 190S is the primary cause of disease/hypovirulence in its natural host and a heterologous host. Virus Res. 2016, 213, 238–245. [Google Scholar] [CrossRef] [PubMed]
  25. Eusebio-Cope, A.; Sun, L.; Tanaka, T.; Chiba, S.; Kasahara, S.; Suzuki, N. The chestnut blight fungus for studies on virus/host and virus/virus interactions: From a natural to a model host. Virology 2015, 477, 164–175. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Chiba, S.; Lin, Y.H.; Kondo, H.; Kanematsu, S.; Suzuki, N. A novel victorivirus from a phytopathogenic fungus, Rosellinia necatrix is infectious as particles and targeted by RNA silencing. J. Virol. 2013, 87, 6727–6738. [Google Scholar] [CrossRef] [PubMed]
  27. Chiba, S.; Suzuki, N. Highly activated RNA silencing via strong induction of dicer by one virus can interfere with the replication of an unrelated virus. Proc. Natl. Acad. Sci. USA 2015, 112, E4911–E4918. [Google Scholar] [CrossRef] [PubMed]
  28. International Rules for Seed Testing 2018. Available online: https://www.seedtest.org/upload/cms/user/ISTARules2018SHmethod7-014_updated20171109.pdf (accessed on 22 June 2019).
  29. White, T.J.; Bruns, T.; Lee, S.; Taylor, J. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In PCR Protocols: A Guide to Methods and Applications; Innis, M.A., Gelfand, D.H., Sninsky, J.J., White, T.J., Eds.; Academic Press: New York, NY, USA, 1990; pp. 315–322. [Google Scholar]
  30. Eusebio-Cope, A.; Suzuki, N. Mycoreovirus genome rearrangements associated with RNA silencing deficiency. Nucleic Acids Res. 2015, 43, 3802–3813. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Taga, M.; Murata, M. Visualization of mitotic chromosomes in filamentous fungi by fluorescence staining and fluorescence in situ hybridization. Chromosoma 1994, 103, 408–413. [Google Scholar] [CrossRef] [PubMed]
  32. Segers, G.C.; Zhang, X.; Deng, F.; Sun, Q.; Nuss, D.L. Evidence that RNA silencing functions as an antiviral defense mechanism in fungi. Proc. Natl. Acad. Sci. USA 2007, 104, 12902–12906. [Google Scholar] [CrossRef] [Green Version]
  33. EnzymeX version 3.3.3. Available online: http://nucleobytes.com/enzymex/index.html (accessed on 22 June 2019).
  34. Sperschneider, J.; Datta, A. DotKnot: Pseudoknot prediction using the probability dot plot under a refined energy model. Nucleic Acids Res. 2010, 38, e103. [Google Scholar] [CrossRef]
  35. Byun, Y.; Han, K. PseudoViewer3: Generating planar drawings of large-scale RNA structures with pseudoknots. Bioinformatics 2009, 25, 1435–1437. [Google Scholar] [CrossRef] [PubMed]
  36. Kondo, H.; Chiba, S.; Suzuki, N. Detection and analysis of non-retroviral RNA virus-like elements in plant, fungal, and insect genomes. Methods Mol. Biol. 2015, 1236, 73–88. [Google Scholar] [PubMed]
  37. Katoh, K.; Standley, D.M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef] [PubMed]
  38. Talavera, G.; Castresana, J. Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Syst. Biol. 2007, 56, 564–577. [Google Scholar] [CrossRef]
  39. Guindon, S.; Dufayard, J.F.; Lefort, V.; Anisimova, M.; Hordijk, W.; Gascuel, O. New algorithms and methods to estimate maximum-likelihood phylogenies: Assessing the performance of PhyML 3.0. Syst. Biol. 2010, 59, 307–321. [Google Scholar] [CrossRef]
  40. Lefort, V.; Longueville, J.E.; Gascuel, O. SMS: Smart Model Selection in PhyML. Mol. Biol. Evol. 2017, 34, 2422–2424. [Google Scholar] [CrossRef] [Green Version]
  41. FigTree version 1.3.1. Available online: http://tree.bio.ed.ac.uk/software/ (accessed on 22 June 2019).
  42. Suzuki, N.; Supyani, S.; Maruyama, K.; Hillman, B.I. Complete genome sequence of Mycoreovirus-1/Cp9B21, a member of a novel genus within the family Reoviridae, isolated from the chestnut blight fungus Cryphonectria parasitica. J. Gen. Virol. 2004, 85, 3437–3448. [Google Scholar] [CrossRef]
  43. Shamsi, W.; Sato, Y.; Jamal, A.; Shahi, S.; Kondo, H.; Suzuki, N.; Bhatti, M.F. Molecular and biological characterization of a novel botybirnavirus identified from a Pakistani isolate of Alternaria alternata. Virus Res. 2019, 263, 119–128. [Google Scholar] [CrossRef]
  44. Hillman, B.I.; Supyani, S.; Kondo, H.; Suzuki, N. A reovirus of the fungus Cryphonectria parasitica that is infectious as particles and related to the Coltivirus genus of animal pathogens. J. Virol. 2004, 78, 892–898. [Google Scholar] [CrossRef]
  45. Chiba, S.; Lin, Y.H.; Kondo, H.; Kanematsu, S.; Suzuki, N. A novel betapartitivirus RnPV6 from Rosellinia necatrix tolerates host RNA silencing but is interfered by its defective RNAs. Virus Res. 2016, 219, 62–72. [Google Scholar] [CrossRef]
  46. Shahi, S.; Eusebio-Cope, A.; Kondo, H.; Hillman, B.I.; Suzuki, N. Investigation of host range of and host defense against a mitochondrially replicating mitovirus. J. Virol. 2019, 93, e01503-18. [Google Scholar] [CrossRef] [PubMed]
  47. CBS-KNAW Collections in short. Available online: http://www.westerdijkinstitute.nl/collections/ (accessed on 3 April 2018).
  48. Ghabrial, S.A.; Nibert, M.L. Victorivirus, a new genus of fungal viruses in the family Totiviridae. Arch. Virol. 2009, 154, 373–379. [Google Scholar] [CrossRef] [PubMed]
  49. Staple, D.W.; Butcher, S.E. Pseudoknots: RNA structures with diverse functions. PLoS Biol. 2005, 3, e213. [Google Scholar] [CrossRef] [PubMed]
  50. Fulbright, D.W. Effect of eliminating dsRNA in hypovirulent Endothia parasitica. Phytopathology 1984, 74, 722–724. [Google Scholar] [CrossRef]
  51. Aoki, N.; Moriyama, H.; Kodama, M.; Arie, T.; Teraoka, T.; Fukuhara, T. A novel mycovirus associated with four double-stranded RNAs affects host fungal growth in Alternaria alternata. Virus Res. 2009, 140, 179–187. [Google Scholar] [CrossRef]
  52. Xavier, A.D.S.; Barros, A.P.O.; Godinho, M.T.; Zerbini, F.M.; Souza, F.O.; Bruckner, F.P.; Alfenas-Zerbini, P. A novel mycovirus associated to Alternaria alternata comprises a distinct lineage in Partitiviridae. Virus Res. 2018, 244, 21–26. [Google Scholar] [CrossRef] [PubMed]
  53. Xiang, J.; Fu, M.; Hong, N.; Zhai, L.; Xiao, F.; Wang, G. Characterization of a novel botybirnavirus isolated from a phytopathogenic Alternaria fungus. Arch. Virol. 2017, 162, 3907–3911. [Google Scholar] [CrossRef]
  54. Ma, G.; Liang, Z.; Hua, H.; Zhou, T.; Wu, X. Complete genome sequence of a new botybirnavirus isolated from a phytopathogenic Alternaria alternata in China. Arch. Virol. 2019, 164, 1225–1228. [Google Scholar] [CrossRef]
  55. Okada, R.; Ichinose, S.; Takeshita, K.; Urayama, S.I.; Fukuhara, T.; Komatsu, K.; Arie, T.; Ishihara, A.; Egusa, M.; Kodama, M.; et al. Molecular characterization of a novel mycovirus in Alternaria alternata manifesting two-sided effects: Down-regulation of host growth and up-regulation of host plant pathogenicity. Virology 2018, 519, 23–32. [Google Scholar] [CrossRef]
  56. Komatsu, K.; Katayama, Y.; Omatsu, T.; Mizutani, T.; Fukuhara, T.; Kodama, M.; Arie, T.; Teraoka, T.; Moriyama, H. Genome sequence of a novel victorivirus identified in the phytopathogenic fungus Alternaria arborescens. Arch. Virol. 2016, 161, 1701–1704. [Google Scholar] [CrossRef]
  57. Komatsu, K.; Katayama, Y.; Omatsu, T.; Mizutani, T.; Fukuhara, T.; Kodama, M.; Arie, T.; Teraoka, T.; Moriyama, H. Genome sequence of a novel mitovirus identified in the phytopathogenic fungus Alternaria arborescens. Arch. Virol. 2016, 161, 2627–2631. [Google Scholar] [CrossRef] [PubMed]
  58. Chen, Y.; Shang, H.H.; Yang, H.Q.; Da Gao, B.; Zhong, J. A mitovirus isolated from the phytopathogenic fungus Alternaria brassicicola. Arch. Virol. 2017, 162, 2869–2874. [Google Scholar] [CrossRef] [PubMed]
  59. Zhong, J.; Shang, H.H.; Zhu, C.X.; Zhu, J.Z.; Zhu, H.J.; Hu, Y.; Gao, B.D. Characterization of a novel single-stranded RNA virus, closely related to fusariviruses, infecting the plant pathogenic fungus Alternaria brassicicola. Virus Res. 2016, 217, 1–7. [Google Scholar] [CrossRef]
  60. Shang, H.H.; Zhong, J.; Zhang, R.J.; Chen, C.Y.; Gao, B.D.; Zhu, H.J. Genome sequence of a novel endornavirus from the phytopathogenic fungus Alternaria brassicicola. Arch. Virol. 2015, 160, 1827–1830. [Google Scholar] [CrossRef] [PubMed]
  61. Hayashi, N.; Tsuge, T.; Kobayashi, H.; Nishimura, S. The presence of double-stranded RNAs in Alternaria alternata Japanese pear pathotype and their participation in AK-toxin productivity. Ann. Phytopath. Soc. Jpn. 1988, 54, 250–252. [Google Scholar] [CrossRef]
  62. Ghabrial, S. Totiviruses. In Encyclopedia of Virology, 3rd ed.; Mahy, B.W.J., Van Regenmortel, M.H.V., Eds.; Elsevier: Oxford, UK, 2008; Volume 5, pp. 163–174. [Google Scholar]
  63. Meyers, G. Translation of the minor capsid protein of a calicivirus is initiated by a novel termination-dependent reinitiation mechanism. J. Biol. Chem. 2003, 278, 34051–34060. [Google Scholar] [CrossRef] [PubMed]
  64. Meyers, G. Characterization of the sequence element directing translation reinitiation in RNA of the calicivirus rabbit hemorrhagic disease virus. J. Virol. 2007, 81, 9623–9632. [Google Scholar] [CrossRef] [PubMed]
  65. Liu, Y.C.; Linder-Basso, D.; Hillman, B.I.; Kaneko, S.; Milgroom, M.G. Evidence for interspecies transmission of viruses in natural populations of filamentous fungi in the genus Cryphonectria. Mol. Ecol. 2003, 12, 1619–1628. [Google Scholar] [CrossRef] [PubMed]
  66. Deng, F.; Xu, R.; Boland, G.J. Hypovirulence-associated double-stranded RNA from Sclerotinia homoeocarpa is conspecific with Ophiostoma novo-ulmi mitovirus 3a-Ld. Phytopathology 2003, 93, 1407–1414. [Google Scholar] [CrossRef]
  67. Yaegashi, H.; Nakamura, H.; Sawahata, T.; Sasaki, A.; Iwanami, Y.; Ito, T.; Kanematsu, S. Appearance of mycovirus-like double-stranded RNAs in the white root rot fungus, Rosellinia necatrix, in an apple orchard. Fems. Microbiol. Ecol. 2013, 83, 49–62. [Google Scholar] [CrossRef]
  68. Kashif, M.; Hyder, R.; de Vega Perez, D.; Hantula, J.; Vainio, E.J. Heterobasidion wood decay fungi host diverse and globally distributed viruses related to Helicobasidium mompa partitivirus V70. Virus Res. 2015, 195, 119–123. [Google Scholar] [CrossRef] [PubMed]
  69. Huang, S.; Ghabrial, S.A. Organization and expression of the double-stranded RNA genome of Helminthosporium victoriae 190S virus, a totivirus infecting a plant pathogenic filamentous fungus. Proc. Natl. Acad. Sci. USA 1996, 93, 12541–12546. [Google Scholar] [CrossRef] [PubMed]
  70. Sun, Q.; Choi, G.H.; Nuss, D.L. A single Argonaute gene is required for induction of RNA silencing antiviral defense and promotes viral RNA recombination. Proc. Natl. Acad. Sci. USA 2009, 106, 17927–17932. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Andika, I.B.; Jamal, A.; Kondo, H.; Suzuki, N. SAGA complex mediates the transcriptional up-regulation of antiviral RNA silencing. Proc. Natl. Acad. Sci. USA 2017, 114, E3499–E3506. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Chiba, S.; Lin, Y.H.; Kondo, H.; Kanematsu, S.; Suzuki, N. Effects of defective-interfering RNA on symptom induction by, and replication of a novel partitivirus from a phytopathogenic fungus Rosellinia. necatrix. J. Virol. 2013, 87, 2330–2341. [Google Scholar] [CrossRef] [PubMed]
  73. Aulia, A.; Andika, I.B.; Kondo, H.; Hillman, B.I.; Suzuki, N. A symptomless hypovirus, CHV4, facilitates stable infection of the chestnut blight fungus by a coinfecting reovirus likely through suppression of antiviral RNA silencing. Virology 2019, 533, 99–107. [Google Scholar] [CrossRef] [PubMed]
  74. Shahi, S.; Chiba, C.; Kondo, H.; Suzuki, N. Cryphonectria nitschkei chrysovirus 1 with unique molecular features and a very narrow host range. 2019. (in preparation) [Google Scholar]
Figure 1. Molecular characterization of Alternaria alternata victorivirus 1 (AalVV1). (A) Colony morphology of a Pakistani Alternaria alternata strain A-16. The fungal isolate was grown on potato dextrose agar for one week and photographed. (B) Agarose gel electrophoresis of AalVV1 (A-16) genomic dsRNA. AalVV1 was separated on 1.0% agarose gel and stained with ethidium bromide (EtBr). The AalVV1 genomic dsRNA of 5.1 kbp is shown by a red arrow. (C) Read mapping of the contig A1-19 (AalVV1). The reads are differentiated by strands with red and green, respectively. A region with significantly low read coverage is shown by a black arrow (D) Genomic organization of AalVV1. The two nonoverlapping open reading frames (ORFs) and the untranslated regions (UTRs) are shown by red and blue boxes and a black line, respectively. ORF 1 (765aa) encodes capsid protein (CP) of 80.2 kDa, whereas ORF2 (718aa) encodes an RNA dependent RNA polymerase (RdRp). The putative slippery site for the −1 frameshift (boxed) is shown at the junction of the two ORF’s with a black arrow along with the start codon of the second ORF. (E) Schematic representation of the translation strategy of the downstream ORF (ORF2) employed by AalVV1. Predicted −1 frameshift site, a putative slippery region (blue bracket) at the junction of the two ORF’s at map positions 2556 to 2560 with two possible H-type pseudoknots structures (2526–2551 and 2575–2608) up- and downstream of the termination codon, respectively, along with the estimated free energy for both structures.
Figure 1. Molecular characterization of Alternaria alternata victorivirus 1 (AalVV1). (A) Colony morphology of a Pakistani Alternaria alternata strain A-16. The fungal isolate was grown on potato dextrose agar for one week and photographed. (B) Agarose gel electrophoresis of AalVV1 (A-16) genomic dsRNA. AalVV1 was separated on 1.0% agarose gel and stained with ethidium bromide (EtBr). The AalVV1 genomic dsRNA of 5.1 kbp is shown by a red arrow. (C) Read mapping of the contig A1-19 (AalVV1). The reads are differentiated by strands with red and green, respectively. A region with significantly low read coverage is shown by a black arrow (D) Genomic organization of AalVV1. The two nonoverlapping open reading frames (ORFs) and the untranslated regions (UTRs) are shown by red and blue boxes and a black line, respectively. ORF 1 (765aa) encodes capsid protein (CP) of 80.2 kDa, whereas ORF2 (718aa) encodes an RNA dependent RNA polymerase (RdRp). The putative slippery site for the −1 frameshift (boxed) is shown at the junction of the two ORF’s with a black arrow along with the start codon of the second ORF. (E) Schematic representation of the translation strategy of the downstream ORF (ORF2) employed by AalVV1. Predicted −1 frameshift site, a putative slippery region (blue bracket) at the junction of the two ORF’s at map positions 2556 to 2560 with two possible H-type pseudoknots structures (2526–2551 and 2575–2608) up- and downstream of the termination codon, respectively, along with the estimated free energy for both structures.
Viruses 11 00577 g001
Figure 2. Maximum likelihood (ML) phylogenetic trees based on the RdRp (A) and CP sequences (B) from members of the genus Totivirus and related viruses. The unrooted ML trees were generated by using PhyML 3.0 [39]. The virus names (blue, approved virus species) and their GenBank accession numbers, and the potential translational stop/restart sequences or the ORF switching modes (frame-shifting) are shown in the trees. The selected members of the genera Leishmaniavirus and Trichomonasvirus, infecting protozoan parasites, were used as the outgroups. The numbers on the branches are the bootstrap values (only values greater than 50% are shown).
Figure 2. Maximum likelihood (ML) phylogenetic trees based on the RdRp (A) and CP sequences (B) from members of the genus Totivirus and related viruses. The unrooted ML trees were generated by using PhyML 3.0 [39]. The virus names (blue, approved virus species) and their GenBank accession numbers, and the potential translational stop/restart sequences or the ORF switching modes (frame-shifting) are shown in the trees. The selected members of the genera Leishmaniavirus and Trichomonasvirus, infecting protozoan parasites, were used as the outgroups. The numbers on the branches are the bootstrap values (only values greater than 50% are shown).
Viruses 11 00577 g002
Figure 3. Components of purified virus particles of AalVV1. (A) Transmission electron micrograph of negatively stained particles of AalVV1. The white scale bar denotes 100 nm. (B) Structural proteins included in the VP fraction of A. alternata strain Ally-12 and A-16. Proteins were denatured in modified Laemmli’s sample buffer containing 0.6% 2-mercaptoethanol and electrophoresed in 10% polyacrylamide gels. The gels were stained with Coomassie Brilliant Blue R-250 (CBB) or silver nitrate. The lane M shows approximate protein molecular size with the Precision Plus Protein Dual Color Standards (Bio-Rad Laboratories, Inc., Hercules, California, USA). The black arrows are pointed to the major 80 kDa protein band subjected to LC-MS/MS analysis.
Figure 3. Components of purified virus particles of AalVV1. (A) Transmission electron micrograph of negatively stained particles of AalVV1. The white scale bar denotes 100 nm. (B) Structural proteins included in the VP fraction of A. alternata strain Ally-12 and A-16. Proteins were denatured in modified Laemmli’s sample buffer containing 0.6% 2-mercaptoethanol and electrophoresed in 10% polyacrylamide gels. The gels were stained with Coomassie Brilliant Blue R-250 (CBB) or silver nitrate. The lane M shows approximate protein molecular size with the Precision Plus Protein Dual Color Standards (Bio-Rad Laboratories, Inc., Hercules, California, USA). The black arrows are pointed to the major 80 kDa protein band subjected to LC-MS/MS analysis.
Viruses 11 00577 g003
Figure 4. Liquid chromatography–tandem mass spectrometry (LC-MS/MS) analysis of the major 80 kDa protein band of AalVV1. (A) The map of peptides obtained by LC-MS/MS analysis of the 80 kDa protein band of AalVV1. The peptides mapped are shown as the black bars on the putative amino acid sequence of frame 3 (267–2564 nt, balck bars in the cyan background) and frame 1 (2560–5043 nt, balck bars in the magenta background) or bridged region between ORF1 and ORF2 (frame 1) on the positive strand of the AalVV1 genome. The vertical axis shows the number of the peptides hit on the indicated positions on the amino acid sequence. The lower diagram represents the AalVV1 genome organization (see Figure 1 legend). The raw data of the LC-MS/MS analysis is shown in Figure S2. (B) The nucleic acid sequence (+ sense RNA) and the amino acid sequence (frame +1 and frame +3) of the bridged region between ORF1 and ORF2 of AalVV1 genome. The blue colored letters indicate the predicted amino acid sequence of ORF1, the red colored letters indicate the predicted amino acid sequence of ORF2, and the beige colored letters indicate the predicted amino acid sequence encoded by the bridged sequence of frame 1 between ORF1 and ORF2. The peptides mapped on the bridged region were indicated in bold letters with underline. The stop codon of ORF1 and the start codon of ORF2 are pointed to with the black arrow heads.
Figure 4. Liquid chromatography–tandem mass spectrometry (LC-MS/MS) analysis of the major 80 kDa protein band of AalVV1. (A) The map of peptides obtained by LC-MS/MS analysis of the 80 kDa protein band of AalVV1. The peptides mapped are shown as the black bars on the putative amino acid sequence of frame 3 (267–2564 nt, balck bars in the cyan background) and frame 1 (2560–5043 nt, balck bars in the magenta background) or bridged region between ORF1 and ORF2 (frame 1) on the positive strand of the AalVV1 genome. The vertical axis shows the number of the peptides hit on the indicated positions on the amino acid sequence. The lower diagram represents the AalVV1 genome organization (see Figure 1 legend). The raw data of the LC-MS/MS analysis is shown in Figure S2. (B) The nucleic acid sequence (+ sense RNA) and the amino acid sequence (frame +1 and frame +3) of the bridged region between ORF1 and ORF2 of AalVV1 genome. The blue colored letters indicate the predicted amino acid sequence of ORF1, the red colored letters indicate the predicted amino acid sequence of ORF2, and the beige colored letters indicate the predicted amino acid sequence encoded by the bridged sequence of frame 1 between ORF1 and ORF2. The peptides mapped on the bridged region were indicated in bold letters with underline. The stop codon of ORF1 and the start codon of ORF2 are pointed to with the black arrow heads.
Viruses 11 00577 g004
Figure 5. Asymptomatic infection of a Japanese strain of Alternaria alternata by AalVV1. (A) Colony morphology of a virus free and AlaVV1-infected substrains derived from a Japanese A. alternata strain Ally-12. Protoplasts of the virus-free Ally-12 were fused with those of the AalVV1-infected a Pakistani strain of A. alternata (A-16). Virus-infected recipient strains with the Ally-12 genetic background were obtained as described in the Materials and Methods section. Virus-free (VF) and virus-infected Ally-12 (V1 and V2) were cultured for five days and photographed. (B) Confirmation of AalVV1 infection of Ally-12. DsRNA was isolated from Ally12-V1 (lane 1), Ally-12-V2 (lane 2), A-16 (lane 3), and Ally-12-VF (lane 4), and electrophoresed in agarose gel and stained by EtBr. The 1-kb DNA ladder was used as size standards.
Figure 5. Asymptomatic infection of a Japanese strain of Alternaria alternata by AalVV1. (A) Colony morphology of a virus free and AlaVV1-infected substrains derived from a Japanese A. alternata strain Ally-12. Protoplasts of the virus-free Ally-12 were fused with those of the AalVV1-infected a Pakistani strain of A. alternata (A-16). Virus-infected recipient strains with the Ally-12 genetic background were obtained as described in the Materials and Methods section. Virus-free (VF) and virus-infected Ally-12 (V1 and V2) were cultured for five days and photographed. (B) Confirmation of AalVV1 infection of Ally-12. DsRNA was isolated from Ally12-V1 (lane 1), Ally-12-V2 (lane 2), A-16 (lane 3), and Ally-12-VF (lane 4), and electrophoresed in agarose gel and stained by EtBr. The 1-kb DNA ladder was used as size standards.
Viruses 11 00577 g005
Figure 6. Infection of AalVV1 to a model filamentous fungal host, Cryphonectria parasitica. (A) Agarose gel electrophoresis of dsRNA from AalVV1-infected C. parasitica. Semipurified AalVV1 was transfected into the RNA silencing-competent C. parasitica EP155 and -deficient Δdcl2 strains. Total RNA extracted from transfectants with either of the genetic backgrounds was electrophoresed in 1.2% agarose gel, and stained by EtBr. AalVV1-infected A. alternata strain A-16 and virus-free A. alternata strain Ally-12 were analyzed in parallel. The AalVV1 genomic dsRNA of 5.1 kbp is shown by a red arrow. M refers to 1-kb DNA-ladder markers. (B) RT-PCR analysis of the EP155 infected by AalVV1. Total RNA fractions were obtained from virus-free Δdcl2 and AalVV1-trasnfected EP155 and Δdcl2, and subjected to RT-PCR. Amplified fragments were analyzed in agarose gel electrophoresis. The AalVV1-specific band corresponding to 1 kbp is shown by a red arrow. The migration positions of 1-kb DNA-ladder markers are shown on the left.
Figure 6. Infection of AalVV1 to a model filamentous fungal host, Cryphonectria parasitica. (A) Agarose gel electrophoresis of dsRNA from AalVV1-infected C. parasitica. Semipurified AalVV1 was transfected into the RNA silencing-competent C. parasitica EP155 and -deficient Δdcl2 strains. Total RNA extracted from transfectants with either of the genetic backgrounds was electrophoresed in 1.2% agarose gel, and stained by EtBr. AalVV1-infected A. alternata strain A-16 and virus-free A. alternata strain Ally-12 were analyzed in parallel. The AalVV1 genomic dsRNA of 5.1 kbp is shown by a red arrow. M refers to 1-kb DNA-ladder markers. (B) RT-PCR analysis of the EP155 infected by AalVV1. Total RNA fractions were obtained from virus-free Δdcl2 and AalVV1-trasnfected EP155 and Δdcl2, and subjected to RT-PCR. Amplified fragments were analyzed in agarose gel electrophoresis. The AalVV1-specific band corresponding to 1 kbp is shown by a red arrow. The migration positions of 1-kb DNA-ladder markers are shown on the left.
Viruses 11 00577 g006
Table 1. BlastP results for Alternaria alternata victorivirus 1 (AalVV1) proteins.
Table 1. BlastP results for Alternaria alternata victorivirus 1 (AalVV1) proteins.
Query/Virus nameQC *E-ValueIdentityAccession
Query: AalVV1 ORF1 (putative capsid protein, CP) **
Phomopsis vexans RNA virus192%0.069.82%YP_009115491.1
Nigrospora oryzae victorivirus 192%0.068.86%YP_009254735.1
Fusarium poae victorivirus 193%0.066.34%YP_009272904.1
Ustilaginoidea virens RNA virus 399%0.064.01%YP_009004155.1
Colletotrichum caudatum totivirus 198%0.062.97%AZT88630.1
Alternaria arborescens victorivirus 192%0.067.37%YP_009553477.1
Coniothyrium minitans RNA virus 92%0.066.95%ALM62230.1
Ustilaginoidea virens RNA virus 694%0.060.00%ART91348.1
Magnaporthe oryzae virus 394%0.060.22%YP_009143306.1
Gremmeniella abietina RNA virus L293%0.061.53%YP_044806.1
Query: AalVV1 ORF2 (putative RNA-dependent RNA-polymerase, RdRp) **
Nigrospora oryzae victorivirus 199%0.061.36%YP_009254736.1
Phomopsis vexans RNA virus99%0.059.53%YP_009115492.1
Ustilaginoidea virens RNA virus 399%0.056.39%YP_009004156.1
Ustilaginoidea virens RNA virus 699%0.054.86%ART91349.1
Fusarium poae victorivirus 199%0.056.39%YP_009272905.1
Colletotrichum caudatum totivirus 199%0.053.69%AZT88631.1
Coniothyrium minitans RNA virus99%0.050.21%YP_392467.1
Phomopsis longicolla totivirus 199%0.050.14%ALD89108.1
Aspergillus mycovirus 17899%0.048.68%ABX79995.1
Magnaporthe oryzae virus 399%0.049.58%YP_009143307.1
*: query coverage; **: selected top ten hits.

Share and Cite

MDPI and ACS Style

Jamal, A.; Sato, Y.; Shahi, S.; Shamsi, W.; Kondo, H.; Suzuki, N. Novel Victorivirus from a Pakistani Isolate of Alternaria alternata Lacking a Typical Translational Stop/Restart Sequence Signature. Viruses 2019, 11, 577. https://doi.org/10.3390/v11060577

AMA Style

Jamal A, Sato Y, Shahi S, Shamsi W, Kondo H, Suzuki N. Novel Victorivirus from a Pakistani Isolate of Alternaria alternata Lacking a Typical Translational Stop/Restart Sequence Signature. Viruses. 2019; 11(6):577. https://doi.org/10.3390/v11060577

Chicago/Turabian Style

Jamal, Atif, Yukiyo Sato, Sabitree Shahi, Wajeeha Shamsi, Hideki Kondo, and Nobuhiro Suzuki. 2019. "Novel Victorivirus from a Pakistani Isolate of Alternaria alternata Lacking a Typical Translational Stop/Restart Sequence Signature" Viruses 11, no. 6: 577. https://doi.org/10.3390/v11060577

APA Style

Jamal, A., Sato, Y., Shahi, S., Shamsi, W., Kondo, H., & Suzuki, N. (2019). Novel Victorivirus from a Pakistani Isolate of Alternaria alternata Lacking a Typical Translational Stop/Restart Sequence Signature. Viruses, 11(6), 577. https://doi.org/10.3390/v11060577

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop