Botryosphaeria dothidea is the important and widespread trunk pathogen, which causes pear ring rot disease in the pear production regions in China. Botryosphaeria dothidea with an extensive host range not only harms pear trees by causing severe symptoms such as stem wart and canker, branch dieback and fruit rot, but also other fruit trees including apple, Jujube and grapevine [1–6]. The ring rot disease is becoming increasingly serious and distributes throughout all pear producing areas at present, which seriously affects the yield and quality of pear to further hind the healthy development of pear industry [6–7]. It is urgent to find safe and effective measure to biological control of pear ring rot diseases. As mycoviruses can infect and proliferate in fungi and some mycoviruses present hypovirulence-associated traits, they are potential biological control agents for the control of these fungal disease. Indeed, many mycoviruses with potential biocontrol value have been found and identified from pathogenic fungi infecting fruit trees [8]. Therefore, it is particularly important to find and identify mycoviruses in Botryosphaeria dothidea strains casing pear ring rot disease. To identify and characterize mycovirus, high through-put sequencing techniques and subsequent bioinformatics have been shown as a very effective approach [9–11]. So far, eleven mycoviruses have been identified in B. dothidea strains and characterized. In our lab, Botryosphaeria dothidea chrysovirus 1 (BdCV1), Bipolaris maydis botybirnavirus 1 strain BdEW220 (BmRV1-BdEW220), Botryosphaeria dothidea RNA virus 1 (BdRV1) involved in hypovirulence [12–14], Botryosphaeria dothidea partitivirus 1 (BdPV1), Botryosphaeria dothidea victorivirus 1 (BdVV1), Botryosphaeria dothidea victorivirus 2 (BdVV2) and Botryosphaeria dothidea botourmiavirus 1 (BdBOV-1) with no obvious effect on host, have been identified in B. dothidea strains causing pear ring rot disease [12, 15–17]. Botryosphaeria dothidea bipartite mycovirus 1 (BdBMV1), Botryosphaeria dothidea mitovirus 1 (BdMV1), Botryosphaeria dothidea fusarivirus 1 (BdFV1) and BdCV1-G1 have been also identified in B. dothidea strains isolated from Kerria japonica, pepper and apple hosts, respectively [18–21]. The Narnaviridae genome has positive single strand RNA (+ssRNA) with approximately 2.2–3.6 kb in length and encodes only one RNA-dependent RNA polymerase (RdRp). The Narnaviridae family is divided into two genera, Narnavirus and Mitovirus (Mitochondrial viruses), which replicate in the cytoplasm and mitochondria, respectively. The mycoviruses in Narnaviridae exist widely in fungi, but it is not clear if they have any effect on all kinds of fungi [22–23]. Here, we report the identification of a novel mitovirus isolated from B. dothidea strain, the causal pathogen of pear ring rot disease. We further determined its +ssRNA full-length genome sequence and molecular characteristics. This mycovirus was tentatively designated as Botryosphaeria dothidea mitovirus 2 (BdMV2).
Provenance and sequencing of BdMV2
B. dothidea strain FJ was isolated from pear stem bark exhibiting wart symptoms in Fujian Province of China, in 2010, and the hypha blocks of FJ strain were stored in 25% glycerin solution at -80℃. The mycelium was grown on potato dextrose agar (PDA) plates overlaid with a layer of cellophane at 25 ℃ in the dark. The mycelium of B. dothidea strain FJ and other Valsa spp strains were together collected to be used for Illumina sequencing (Beckman company, Beijing, China). The obtained high-quality clean reads were assembled contigs, which were analyzed to search mycoviruses using Blastx. Among these contigs, Contig1603 with 2.5 kb in length shares the highest sequence similarity with 67.14% in comparison to that of Macrophomina phaseolina mitovirus 4. Total RNA was extracted by Trizol method (Aidlab Biotechnologies Co., Ltd, Beijing, China) as template for synthesize cDNA. Finally, the putative mycoviral genome intermediate sequence was determined by RT-PCR and obtained 2042 nt by sequencing in B. dothidea strain FJ using the designed primers based on Contig1603 sequence (Supplementary Table 1). The dsRNA was extracted from the cultured hypha by the column separation method [24]. The adaptor RACE-OLIGO was ligated to the 3’ terminus of each strand of dsRNA, which were used for template to reversely transcribed cDNA using Oligo REV primer as described previously [25]. The terminal sequence amplification was performed by nested PCR using the primers designed based on obtained viral sequence and sequence-specific adaptor primers of O5RACE-2 and O5RACE-3, respectively (Supplementary Table 1). The obtained PCR product was purified, linked to pMD18-T vector (Takara, Dalian, China), and transformed into Escherichia coli Top10. For each PCR product, at least three independent identified positive clones were sequenced in both directions (Sangon Biotech Co., Ltd, Shanghai, China). Sequences obtained from clones were assembled using DNAMAN. The viral ORF was predicted using the ORF finder online website of https://www.ncbi.nlm.nih.gov/orffinder. Mycoviral sequences search was performed using NCBI database. Multiple mycoviral sequences were aligned using CLUSTAL X [26]. Phylogenetic analysis was conducted by using NJ method with 1000 replicates of MEGA 7.0, and RdRp domain motif sequences were analyzed using the GeneDoc software [27–28]. The secondary structures of BdMV2 genome terminus were predicated by online RNAfold software.
Sequence properties
The colony morphology of B. dothidea strain FJ was shown, which exhibited normal mycelia and the mycelia could grow to cover the whole plate in 5 days (Fig. 1A). RT-PCR and RNA-ligase-mediated rapid amplification of cDNA ends (RLM-RACE) in combination with cloning and sequence assembly were used to obtain the BdMV2 full-length genome sequence from B. dothidea strain FJ. Its complete genome sequence was deposited in the NCBI GenBank database (Acc. No. MW553911), which was 2538 nt in length with a CG content of 47.24% for the genome. BdMV2 contains a single big ORF with AUG as the initiation condon at positions 176–178 nt and UAA as the termination condon at positions 2243–2245 nt of the genome by fungal mitochondrial codon usage. The ORF (2,070 nt in length) potentially encoded 689 aa RdRp protein with predicated molecular weight of 77.35 kDa. The 5’ and 3’ untranslated regions (UTRs) are 175 and 293 nt long, respectively (Fig. 1B). Using the RNAfold program, two potential stem-loop structures in the 5’ and 3’ UTRs and a potential panhandle structure of the BdMV2 positive-sense strand UTRs were predicted, respectively (Fig. 1C), which showed typical characteristic of reported members of the genus Mitovirus [22]. The BdMV2 contains a conserved RdRp domain at positions 164–495 aa of the genome (Mitovirus RNA-dependent RNA polymeras; pfam05919) by the Conserved Domain Database (CDD) search of NCBI. Multiple sequence alignment analysis of RdRp domain from BdMV2 with those members in Mitovirus showed that BdMV2 RdRp domain contains all six conserved motifs (I to VI), which are characteristics of members in the Mitovirus (Fig. 2A). A BLASTp comparison of the RdRp sequence of BdMV2 showed the highest similarity of 69.36% with that of RsMV10, and had identity 28.91–68.79% (query coverage more than 90%) to those of members in the genus Mitovirus (Supplementary Table 2). To determine the taxonomic status of BdMV2, phylogenetic analysis was performed using the RdRp aa sequences encoded by BdMV2 and other selected mycoviruses in Mitovirus and Narnavirus of the family Narnaviridae. The results showed that all mitoviruses were classified into two clades of clade I and clade II. It revealed that BdMV2 existed in a separate branch, which was closely clustered in the same subclade with those of RsMV10 and MpMV4 belonging to clade I in Mitovirus genus, distinct from members of the genus Narnavirus (Fig. 2B). In conclusion, based on its genome organization, sequence alignment and phylogenetic analysis, BdMV2 is a novel mitovirus in the B. dothidea pathogen causing pear ring rot disease.