C. sacchari 125 was originally isolated from diseased sugarcane that was collected from Khuzestan province, Iran (Fotouhifar, 2007). Isolate 125 was identified as a strain of C. sacchari by amplifying and sequencing of the internal transcribed spacer (ITS) region, ITS1-5.8S-ITS2, and comparing the sequence (GenBank accession no. EF447388) with those of other strains in Genbank (Fotouhifar et al., 2010). C. sacchari 125 was cultured on potato dextrose agar (PDA) at 25-27°C for 4-5 days. Then five mycelial agar discs from the PDA plate were inoculated to 150 ml potato dextrose broth (PDB) and cultured for 3-4 days at 25-27°C with shaking. The biomass was extracted with Buchner funnel, lyophilized and then kept in −20°C for extracting dsRNA.

Viral dsRNA was extracted from isolate 125 by CF11 (Whatman, UK) cellulose chromatography as described by Dodds et al. (1984). In brief, mycelial tissue was ground by mortar and pestle in liquid nitrogen to a fine powder. Total nucleic acids were extracted from 40 mg of mycelial powder using 800 μl 2×STE buffer (20 ml of 1 M Tris (pH 8.0), 40 ml of 5 M NaCl, 4 ml of 0.5 M EDTA (pH 8.0) and 936 ml of sterile distilled water) and proteins were precipitated with 80 μl 10% sodium dodecyl sulfate (SDS) and 800 μl phenol-chloroform-isoamylalcohol (12:12:1). The resulting mixture was centrifuged at 5,000×g for 10 min, and exactly 700 μl of the aqueous phase was mixed with 700 μl chloroform-isoamylalcohol (24:1), followed by centrifugation to separate phases for 5 min, 13,000×g. Exactly 600 μl of the aqueous phase was transferred to a 1.5 ml microfuge tube, and the extracted nucleic acids adjusted with absolute ethanol to a final concentration of 15%. The dsRNA was purified by chromatography on CF11 cellulose; the cellulose CF-11 was washed with 1×STE/15% ethanol and the dsRNA was eluted with 1×STE. The dsRNA was precipitated in 0.3 M sodium acetate (pH 4) and 70% ethanol and was collected by centrifugation. The pellet was rinsed with 70% ethanol. The final dsRNA pellet was dissolved in 20 μl RNase-free water and digested with DNase I and RNase A at high and low salt concentrations as described by Howitt et al. (1995). The final dsRNA was electrophoresed on a 0.8% agarose gel in TBE buffer stained with ethidium bromide, and observed under UV illumination. The dsRNA sample was stored at −20°C for further steps.

Cycloheximide treatments were used to eliminate dsRNA in isolate 125. These treatments were carried out on PDA plates containing 25 and 50 μg cycloheximide (Sigma, USA) per milliliter (Hammar et al., 1989). Mycelial plugs of isolate 125 were transferred onto treated PDA plates. Two replicates of each concentration were incubated at 25-27°C. After 3 days, hyphal tips from the margins of the colonies were transferred to fresh PDA. This was repeated two times. Then mycelial plugs were cultured in PDB for 3-4 days. The dsRNA was extracted to investigate the effect of cycloheximide and hyphal tip culture on mycelial infection as described above.

DsRNA from isolate 125 was used for cDNA synthesis with tagged random hexamer primers using a protocol based on that of Roossinck et al. (2010). A 9 μl aliquot of dsRNA was mixed with 1 μl 10 mM TE buffer (10 mM Tris, 10 mM EDTA, pH 7.5) and 2 μl (20 μM) random primer 5′-CCT GAA TTC GGA TCC TCC N6-3′. The dsRNA was placed 2 min in boiling water for denaturing, and then chilled on ice. Eight μl of reverse transcription mixture containing 4 μl 5× superscript buffer, 1 μl dNTPs (10 mM each), 2 μl dithiothreitol (supplied by manufacturer) and 1 μl Superscript III (Invitrogen, USA) was added to the sample, and the mixture was incubated on ice for 15 min, and then placed at 50°C for 2 h. The samples were removed from 50°C and treated with 1 μl (10 mg/ml) of boiled ribonuclease A (Sigma, USA), and incubated at room temperature for 15 min. The sample then was heated to 85°C for 2 min, and the primers removed using a cycle pure kit (Omega, Bio-Tek, USA) according to the manufacturer’s instructions. The sample was eluted in 30 μl water. A 1.5 μl aliquot of the cDNA was amplified in a total volume of 50 μl containing a final concentration of 10× PCR buffer (30 mM MgCl₂, Idaho Technologies), 0.2 mM dNTPs, 2 μM Tag primer 5′-CCT GAA TTC GGA TCC TCC-3′, and 1 unit of Taq DNA polymerase (Promega, USA). The amplification was done in an Idaho Technologies Rapid Cycler II, using capillary tubes and the following program: 94°C for 1 min; 65°C 0 s; 72°C 5 s, with a slope of 9, followed by 40 cycles of 94°C 0 s; 45°C 0 s; 72°C 30 s, with a slope of 5, and a final 5 min at 72°C and 5 min at 37°C. The amplified cDNAs were removed from capillary tubes, and were viewed after separation on 1.2% agarose gel and staining with ethidium bromide.

For cloning, the double stranded cDNAs were purified using a cycle pure kit (Omega, Bio-tek), and were ligated with pGEM-T vector using T4 DNA ligase according to manufacturer’s instructions (Promega, USA). The recombinant vectors were transformed into competent cells of Escherichia coli strain JM101. The inserted clones were screened and amplified with M13 forward and reverse universal primers, using the program: 94°C for 15 s, followed by 30 cycles of 94°C 15 s; 54°C 6 s; 72°C 3 min and a final 5 min at 72°C and 5 min at 37°C. The amplified clones were sequenced using Applied Biosystems 3730XL. All clones were sequenced in both directions.

PCR was conducted to determine missing sequence gaps using specific primers, designed to amplifying a 300 bp fragment (Table 1), and Taq DNA polymerase. The PCR conditions included an initial denaturation step of 15 s at 94°C, followed by 35 cycles of 15 s at 94°C, 6 s at 53°C and 45 s at 72°C, with a final of 5 min at 72°C and 5 min at 37°C. PCR products, with a predicted size of 300 bp, were fractionated by gel electrophoresis on 1.2% agarose gels and stained with ethidium bromide. PCR products were cloned and sequenced as described above.

Sequence data were assembled using Geneious software version 6.1.2 (Drummond et al., 2011). The results were compared to databases with the both BlastN and BlastX against the GenBank nucleotide (nt) and non-redundant protein (nr) databases, respectively. The BLAST program is available in The National Center for Biotechnology Information (NCBI). The putative amino acid sequence of the isolate 125 RdRp gene was compared to databases with BlastP. This putative RdRp and those of closely related viruses were obtained from the NCBI database, and aligned using the ClustalW program. After hand editing of the alignments these were used for phylogenetic analyses and tree construction using MrBayes (Huelsenbeck and Ronquist, 2001) via the Geneious plugin. The rate matrix was set to a Poisson distribution with a gamma rate variation. The total chain length was 10,000 and branch lengths were unconstrained.

Viruses of the family Partitiviridae have bi-segmented genomes containing two linear dsRNA molecules, which encode two proteins: the smaller segment encodes the CP, and the larger one encodes the RdRp (Nibert et al., 2013). Extracting dsRNA from mycelial tissue of C. sacchari 125, recovered from sugarcane in Khuzestan, Iran, indicated a potential mycovirus infection. The electrophoresis analysis showed the extracted dsRNAs consist of more than the two linear dsRNAs of approximately 1.85, and 1.65 kbp that are presumed to be the viral genomic elements; they also contain a smaller dsRNA segment (about 1.27 kbp) that could be a satellite dsRNA or defective-interfering dsRNA (King et al., 2012) (Fig. 1). In some partitiviruses such as those infecting Rosellinia necatrix, Discula destructive, Gremminiella abietina and Atkinsonella hypoxylon, more than two linear dsRNAs, have been detected, thought to be satellites or defective RNAs (Chiba et al., 2013; Ghabrial et al., 2008). Analysis of the effects of subviral RNAs in partitiviruses are difficult because no reverse genetics system has been established (Chiba et al., 2013).

For confirming the identity of these fragments as dsRNAs, they were treated with DNase I and RNase A in high and low ionic strength buffer. These fragments were resistant to treatment by DNase I and RNase A in high salt concentrations, but they were sensitive to RNase A treatment in low salt concentrations (not shown). These results confirmed the identity of these elements as dsRNAs.

Many methods for eliminating dsRNA viruses from fungi have been reported (Márquez et al., 2007; Hammar et al., 1989). These methods met with varying degrees of success. Hammar et al. (1989) used cyclohexamide to eliminate dsRNAs in a strain of Leucostoma persoonii, some segments of dsRNAs were eliminated and partially cured strain was obtained. Márquez et al. (2007) successfully eliminated dsRNAs from Curvularia protuberata by desiccation and freeze-thawing. Elias and Coty (1996) used a single conidial transfer method for removing dsRNAs from six isolates of Aspergillus section flavi; except one isolate, these were cured completely or partially.

We attempted to eliminate CsPV from isolate 125 using treatment with cycloheximide in two concentrations (25 and 50 μg) and hyphal tip culture. The results showed that neither 25 μg nor 50 μg of cycloheximide affected the potential mycovirus infection in isolate 125, and also hyphal tip culture had no effect on infected mycelia, and the dsRNAs were retained in C. sacchari isolate 125 (data not shown). Details of mycovirus biology are still scant, but it is clear that there is a great deal of variation, with some viruses being extremely unstable in culture, and others being recalcitrant to curing (M. Roossinck, unpublished). This makes detailed analysis of the effects of these viruses on their hosts very challenging.

The partial nucleotide sequence of isolate 125 was obtained by sequencing recombinant plasmid clones derived from a cDNA library. Twenty-four cDNA plasmid clones were randomly selected for sequencing in two directions. Sequence assembly of these cDNA clones and comparison with known mycoviruses using the BLAST program indicated that there was a gap, which was filled by PCR with specific primers, and we obtained a 1,284 base pair (bp) sequence. The analysis of the sequence with BlastN and BlastX programs revealed significant sequence similarities with the RNA-dependent RNA polymerase (RdRp) of viruses of the family Partitiviridae. Sequence similarity among partitiviruses and the obtained sequence of isolate 125 indicate that it is a novel virus in the genus Partitivirus in the family Partitiviridae; we have tentatively designated it as Cytospora sacchari partitivirus (CsPV). The partial nucleotide sequence of CsPV has been deposited in GeneBank (accession no. KF697164).

The obtained sequence of CsPV can potentially encode a 405 amino acid (aa) protein. The search in BLASTP showed that this protein has similarity with RdRps of a number of fungal partitiviruses. Pairwise aa identities between CsPV and previously identified partitiviruses are presented in Table 2. The highest aa identity (76%) is between CsPV and Penicillium stoloniferum virus S (PsV-S). Blast search with deduced amino acid sequence reveals that this sequence has a single ORF. Search of the conserved domain database (CDD) revealed a significant match with conserved domains related to the Partitivirus RdRp (RNA-dep-RNAP, cd01699). To define the relationship among the putative RdRp of this virus and the previously characterized partitiviruses, multiple sequence alignment was conducted. The alignment shows four amino acid conserved domains (I-IV) characteristic of partitivirus RdRps including the highly conserved GDD motif typical of RdRps in motif IV (Fig. 2).

A Bayesian phylogenetic analysis was conducted based on alignments of the putative amino acid sequence of RdRps of CsPV and selected members of the genus Partitivirus (Fig. 3). As expected, CsPV was included in a clade within the partitiviruses. Based on this analysis CsPV was clustered with Aspergillus ochraceous virus 1 (ABC86749.1) in a clade with 0.9631 posterior probability.

Based on a recent phylogenetic analysis it seems that the taxonomy of the family Partitiviridae is due for some revisions (Crawford et al., 2006; Nibert et al., 2013). In the current classification of Partitiviridae, host range has an important role for genus divisions; fungal partitiviruses have been classified in the genus Partitivirus while plant partitiviruses have been placed in Alphacryptovirus or Betacryptovirus. However the discovery of plant-infecting partitiviruses in Primula malacoides (Li et al., 2009), and Cannabis sativa (Ziegler et al., 2012), and the recent phylogenetic results suggest that this criterion may not be appropriate for genus divisions in the family Partitiviridae (Nibert et al., 2013; Roossinck, 2010). The existence of plant-infecting viruses in the genus Partitivirus supports the possibility of horizontal transfer of these viruses between plants and fungi (Crawford et al., 2006; Ghabrial et al., 2008; Roossinck, 1997; Roossinck, 2010). In addition, evidence has been recently found for integration of partiti-virus sequences into plant genomes (Chiba et al., 2011; Liu et al., 2010). Phylogenetic analysis based on aa sequences of RdRps of members of family Partitiviridae led to identification of two clusters in the genus Partitivirus, comprised of a mixture of fungal viruses and plant viruses (Ghabrial et al., 2008; Nibert et al., 2013), these two clusters may be marked as two subgroups (subgroup 1 and 2) in this genus (Crawford et al., 2006; Ghabrial et al., 2008). One of the taxonomic revisions based on these results is to divide the genus Partitivirus into two new genera (Nibert et al., 2013), although other revisions of the Partitiviridae family are probably also necessary (Roossinck, 2010).

Phylogenetic analysis of aa sequences of RdRps of the established members of genus Partitivirus and putative RdRp of CsPV showed that CsPV would be placed in one of these subgroups that contains fungal viruses, whereas the other one consists of a mixture of fungal viruses and plant viruses. Based on suggested revision, each of these subgroups could be considered as a new genus in the family Partitiviridae. CsPV with AoV1 form a cluster in one of these new genera.

In summary, we have identified a putative partitivirus associated with C. sacchari that we tentatively named CsPV. The attempts for removing the dsRNAs with treatment by cyclohexamide and hyphal tip culture were unsuccessful. This study is the first report of partitivirus infection in the plant pathogenic fungus C. sacchari.