Four Novel Mycoviruses from the Hypovirulent Botrytis cinerea SZ-2-3y Isolate from Paris polyphylla: Molecular Characterisation and Mitoviral Sequence Transboundary Entry into Plants

A hypovirulent SZ-2-3y strain isolated from diseased Paris polyphylla was identified as Botrytis cinerea. Interestingly, SZ-2-3y was coinfected with a mitovirus, two botouliviruses, and a 3074 nt fusarivirus, designated Botrytis cinerea fusarivirus 8 (BcFV8); it shares an 87.2% sequence identity with the previously identified Botrytis cinerea fusarivirus 6 (BcFV6). The full-length 2945 nt genome sequence of the mitovirus, termed Botrytis cinerea mitovirus 10 (BcMV10), shares a 54% sequence identity with Fusarium boothii mitovirus 1 (FbMV1), and clusters with fungus mitoviruses, plant mitoviruses and plant mitochondria; hence BcMV10 is a new Mitoviridae member. The full-length 2759 nt and 2812 nt genome sequences of the other two botouliviruses, named Botrytis cinerea botoulivirus 18 and 19 (BcBoV18 and 19), share a 40% amino acid sequence identity with RNA-dependent RNA polymerase protein (RdRp), and these are new members of the Botoulivirus genus of Botourmiaviridae. Horizontal transmission analysis showed that BcBoV18, BcBoV19 and BcFV8 are not related to hypovirulence, suggesting that BcMV10 may induce hypovirulence. Intriguingly, a partial BcMV10 sequence was detected in cucumber plants inoculated with SZ-2-3y mycelium or pXT1/BcMV10 agrobacterium. In conclusion, we identified a hypovirulent SZ-2-3y fungal strain from P. polyphylla, coinfected with four novel mycoviruses that could serve as potential biocontrol agents. Our findings provide evidence of cross-kingdom mycoviral sequence transmission.


Introduction
High-throughput sequencing and bioinformatics have revealed the presence of numerous mycoviruses in filamentous fungi, especially mitovirus and ourmia-like viral sequences [1][2][3], and the observed genetic variability has implications for understanding the evolution of viruses in fungi [4][5][6][7]. Mitoviruses and ourmia-like mycoviruses belong to Mitoviridae and Botourmiaviridae, respectively. Members of these two families are positive-sense single-stranded RNA (+ssRNA) viruses, with genome sequences ranging from 2.2 to 3.6 kb that only encode a single and essential RNA-dependent RNA polymerase protein (RdRp) [8][9][10]. Botourmiaviridae includes six genera, in which Ourmiavirus infects plants, Scleroulivirus infects both fungi and plants, and the other four genera (Botoulivirus, Magoulivirus, Penoulivirus and Rhizoulivirus) mainly infect fungal hosts such Paris polyphylla is a horticultural crop with high value in traditional Chinese medicines. In June 2018, various P. polyphylla plants were killed by grey mold disease in Sui County, Hubei Province, China. The SZ-2-3y strain was isolated from the diseased leaves of P. polyphylla and further identified as a B. cinerea strain. The SZ-2-3y isolate showed relatively normal morphology, but exhibited hypovirulent characteristics in in vitro pathogenicity tests, including the leaves of P. polyphylla and tobacco (Nicotiana benthamiana), and the fruits of tomato and apple. By screening viral contigs and dsRNA patterns, and performing RT-PCR, two ourmia-like viruses designated BcBoV18 and BcBoV19, one mitovirus designated BcMV10, and a fusarivirus designated BcFV8, were identified in SZ-2-3y.
To further explore the potential of these mycoviruses for biological control, the molecular characterisation of mycoviruses and their biological functions in the B. cinerea strain SZ-2-3y were conducted in the present work. Interestingly, BcMV10 sequences were also detected in the new leaves of cucumber plants inoculated with SZ-2-3y mycelium, as well as agrobacterium-mediated inoculation of pXT1/BcMV10. The findings provide new insight into the introduction of mitovirus sequences from fungi into plants, mediated by HGT.

Isolation of the Pathogenic SZ-2-3y Strain and Detection in P. polyphylla
From April to June 2018, various severely diseased P. polyphylla plants with a grey mold layer on their flowers and leaves were identified in Sui County, Hubei, China. The pathogen was isolated from the diseased leaves and flowers. Genomic DNA was extracted from the above isolated mycelium and designated the SZ-2-3y strain, based on the Cetylramethylammonium bromide (CTAB) method, with slight modifications. Internal transcribed space (ITS)-rDNA region sequence analysis, in combination with PCR detection using B. cinerea-specific PCR primers, was used to identify Botrytis species [45]. Meanwhile, the B. cinerea B05.10 isolate was provided by Prof. Li [15].

Characterisation of SZ-2-3y Strain Growth Morphology and Pathogenicity
The isolated SZ-2-3y strain mycelial plugs, collected from actively growing colonies on 2-day-old agar plates, and the B05.10 reference isolate were incubated on potato dextrose agar (PDA) at 20 • C and 25 • C in darkness, to determine radial growth rates, sclerotia production, and quantities over 12 days. The colony diameter was measured by the cross method to calculate the growth rate. For determination of SZ-2-3y pathogenicity, fresh mycelial agar plugs with 5 mm diameter colonies were placed on fully expanded leaves detached from potted P. polyphylla and N. benthamiana, on tomato fruits from a supermarket by intact inoculation, and on apple fruits (Malus mali) from a supermarket by wound inoculation. Meanwhile, the B05.10 reference strain and blank PDA agar inoculation served as controls. The samples were placed in a plastic tray with moist white cotton gauze and a transparent plastic covering to retain humidity, and the diseased lesion diameter was measured at 24 h intervals. The results were analysed by one-way analysis of variance (ANOVA) at a significance level of 0.05 using SPSS software (version 22.0) to determine differences between treatments, and the number of sclerotia and lesion diameters were analysed by Student's t test at p < 0.05 and p < 0.01, respectively.

DsRNA Extraction and Detection
Strain SZ-2-3y hyphae were cultured on PDA covered with cellophane for 5 to 7 days. The mycelium mass was harvested and quickly ground to a fine powder in liquid nitrogen. Two methods were employed to extract dsRNA: a cross-column method developed in our lab, and a cellulose (CF-11) chromatography method reported previously [46,47]. The obtained dsRNAs were separated by 1.2% agarose gel electrophoresis to visualise patterns, and stored for subsequent experiments.

High-Throughput Sequencing and Analysis
High-throughput sequencing was used to identify novel mycoviruses in fungi. Four mycelium samples from three Alternaria spp. strains and SZ-2-3y from the present study were used for Illumina sequencing to screen for mycoviruses. Total RNAs from mixtures of cultured mycelial mass were extracted and further assessed by Nanodrop and Agilent 2100 instruments, and by gel electrophoresis. The rRNA-depleted RNA library was constructed using the Epicentre Ribo-Zero rRNA removal kit (Epicentre, Madison, WI, USA) and a TruSeq RNA Sample Prep Kit v2 (Illumina, San Diego, CA, USA). Sequencing was performed by the Illumina HiSeq X-ten sequencing platform with PE150 bp (Biomarker Technologies, Beijing, China). Clean reads of a high quality were obtained by filtering low-quality raw data, including adapter-contaminated reads, reads with unknown bases (N), and short reads. After quality trimming, de novo assembly was performed with IDBA_ud software. Subsequently, assembled contigs were further analysed to search for homology with mycoviral sequences using BLASTn and BLASTx in the NCBI database. The obtained mycoviral sequences from the library presence were detected and confirmed by RT-PCR for SZ-2-3y in this study (Supplementary Table S1).

Full-Length Cloning of Mycovirus cDNA
Intermediate fragments of mycoviruses were obtained by reverse transcription PCR (RT-PCR) using virus-specific primers, designed based on the assembled mycovirus contig sequences. The 5 -terminal and 3 -terminal sequences of dsRNA were determined as previously described [48]. The adaptor primer RACE-OLIGO was ligated to the 3 -terminus of each dsRNA strand using T4 RNA ligase (New England BioLabs, Beijing, China) at 25 • C, for 2 h. The adaptor-ligated dsRNAs were purified, reverse transcribed using O5RACE-1 primer (the RT primer complementary to the adaptor primer RACE-OLIGO sequence). Synthetic cDNAs were amplified using primers O5RACE-2 or O5RACE-3, complementary to the RNA ligation oligonucleotide and mycovirus sequence-specific primers corresponding to the 5 -and 3 -terminal sequences of dsRNAs, to perform amplification of viral cDNA ends. The experiment was repeated at least three times independently. Primers used in this study are listed in Supplementary Table S1. All PCR products were purified and ligated into the pMD18-T vector (Takara, Dalian, China), then transformed into Escherichia coli Top 10 competent cells for DNA sequencing.

Analysis of Mycoviral Genome Sequences
The sequences obtained from clones were assembled using DNAMAN. The viral open reading frames (ORFs) were predicted using ORF finder (https://www.ncbi.nlm. nih.gov/orffinder, accessed on 28 July 2021). The mycoviral sequences were searched for using the NCBI database. Multiple mycoviral sequences were aligned using MAFFT (https://mafft.cbrc.jp/alignment/server/index.html, accessed on 28 July 2021). The phylogenetic analysis was conducted via neighbour-joining (NJ) and maximum likelihood (ML) methods with 1000 replicates using MEGA 5.2, and the RdRp motif sequences were analysed using GeneDoc software [49,50]. The bootstrap values (>50%) are labelled on branches. The information for the selected mycoviral isolates used in the phylogenetic tree and multiple alignment comparisons of RdRp amino acid (aa) sequences is listed in Supplementary Tables S2-S5. Potential secondary structures of terminus sequences from the three mycovirus genomes were predicated by RNAfold (http: //mfold.rna.albany.edu/?q=DINAMelt/Quickfold, accessed on 28 July 2021) [51]. The resulting RNA reads were mapped to mycoviral genomes using HISAT version 2.22 [52]. The number of mapped reads for each mycovirus was calculated by Samtools [53], and duplicates were identified using Picard version 2.2.4 (https://broadinstitute.github.io/picard/, accessed on 28 July 2021). The sequencing depth of viral sites was calculated and visualised using Circos software.

Horizontal Transmission of Mycoviruses in SZ-2-3y
The dsRNA-containing hypovirulent strain SZ-2-3y, as the donor, was horizontally transmitted to the virulent isolate B05.10, as the recipient, on PDA using the pair-culturing technique. From colonies of each recipient in 7-day-old cultures, six derivatives (Yb-1, Yb-2, Yb-3, Yb-7, Yb-8 and Yb-11) were tested for the presence of the four mycoviruses by RT-PCR with mycovirus-specific primer pairs (Supplementary Table S1). Additionally, these derivatives were tested for biological characteristics, namely the growth rate and phenotype on PDA, and for pathogenicity by intact inoculation on tobacco leaves using the above-described methods. Parental SZ-2-3y and B05.10 strains were included as controls.

SZ-2-3y Mycelium Inoculation on Plants
For inoculation of the SZ-2-3y strain to cucumber plants, mycelia-containing PDA plugs, picked up from the edge of 2-or 3-day-old culture colonies with active infection ability, were placed on the cotyledons of cucumber plants. The inoculated parts of the leaves were wrapped with plastic wrap for 48 h. The healthy cucumber plants and plants inoculated with blank PDA or B05.10 served as controls. The inoculated plants were grown in pots under a 16 h light/8 h dark cycle at 24 ± 2 • C, in a growth chamber. Two or three independent experiments were conducted. Cucumber seeds from the 'Suyo' cultivar (cv) were provided by Prof. Shifang Li from the Chinese Academy of Agricultural Science (CAAS), and seeds of the 'Cuiyu' cultivar were bought from online supermarkets.

Construction of the pXT1/BcMV10 Vector
The modified binary plant expression vector pXT1 (designated as pCB301-2 × 35S-MCS-HDVRZ-NOS), derived from the original vector pCB301 and used as the backbone for the construction of infectious clones, was kindly supplied by Prof. Xiaorong Tao (Nanjing Agricultural University) [54]. The homologous recombination method was employed for construction [55,56]. The full-length BcMV10 cDNA was amplified, including the 30 nt vector sequence; the 5 -end of the first viral genome fragment was joined through a 15 nt overlap to the end of the 35 S promoter of pXT1, and the 3 -end viral fragment was joined through a 15 nt overlap to the 5 end of the HDV ribozyme of pXT1 (Supplementary Table S6). The pXT1 plasmid was linearised with restriction enzymes StuI and SmaI (Takara, Dalian, China). The resulting PCR fragments and linear pXT1 were purified to generate the full-length BcMV10 clone using a ClonExpress II One Step Cloning Kit (Vazyme Biotech Co. Ltd., Nanjing, China) according to the manufacturer's procedure, with slight modifications. The obtained assembled DNA products were transformed into E. coli strain Top10, and positive colonies were confirmed by DNA sequencing. The resulting recombinant plasmid (named pXT1/BcMV10) was used to transform Agrobacterium tumefaciens GV3101 by heat-shock transformation for infiltration experiments [57].

Ago-Infiltration of Cucumber Plants with pXT1/BcMV10
The obtained positive plasmids designated pXT1/BcMV10, and the pXT1 control, were separately transformed into GV3101, and the Agrobacteria cells were grown at 28 • C for 2 days on LB solid medium with 50 µg mL −1 kanamycin and 25 µg mL −1 rifampicin antibiotics. The positive monoclonal culture was shaken at 28 • C overnight, and cells were washed with sterile water three times after centrifugation. The culture was suspended in inoculation buffer including 0.1 M MgCl 2 , 0.1 M 2-Morpholineethanesulfonic acid (MES), 0.015 M Acetosyringone (AS) and the cell density was adjusted to an OD600 value of 1.0. The culture was incubated for 4 h at room temperature then used for inoculation. The cultures (300 µL) were agroinfiltrated into the abaxial surface of two cotyledons from 'Cuiyu' and 'Suyo' seedlings, using 1 mL syringes without a needle. The inoculated plants were grown under 16 h/8 h light/dark conditions at 24 • C. Two or three independent experiments were conducted.

Total RNA Extraction and RT-PCR Analysis of Inoculated Plants
At about 2 weeks post-inoculation (dpi), the leaves and roots of cucumber shoots were evaluated to detect mycoviral sequences by RT-PCR using viral-specific primer pairs (Supplementary Table S6). Total RNAs were extracted using TRIzol (Aidlab Biotechnologies Co., Ltd., Beijing, China) from: the local and new leaves or roots of cucumber plants, inoculated with pXT1/BcMV10 or SZ-2-3y mycelium; new leaves of control plants, inoculated with PDA or pXT1; and healthy cucumber plants as negative controls. cDNAs were synthesised by M-MLV reverse transcriptase (Promega, Madison, WI, USA). BcMV10 sequence-specific primers and 2× DNA polymerase mixture (ComWin Biotech Co., Ltd., Beijing, China) were used for RT-PCR, and the resulting PCR products were verified by gel electrophoresis analysis and DNA sequencing.

Identification and Biological Characteristics of the Hypovirulent SZ-2-3y Strain
In June 2018, some severely diseased P. polyphylla plants were identified in Sui County, Hubei, China, with obvious grey mold symptoms on their leaves and flowers (Supplementary Figure S1). The SZ-2-3y strain was isolated from P. polyphylla leaves with typical grey mold disease, and identified based on cultural morphology and PCR detection. The rDNA-ITS nucleotide sequence from the mycelium of the SZ-2-3y strain was 99% identical to that of the B. cinerea strains. Meanwhile, a 327 bp specific DNA fragment was amplified using reported B. cinerea-specific primers Bc-f/Bc-r and sequenced, revealing that the SZ-2-3y isolate was indeed B. cinerea (Supplementary Table S1). Compared with B05.10 controls-with 1.9 cm/day after incubation at 20 • C, for 5 days, in darkness-the colony morphology of the SZ-2-3y strain was normal, with a slightly slower growth rate of 1.6 cm/day ( Figure 1A,B); furthermore, more sclerotia were produced (143/plate compared with 50/plate for B05.10) after 12 days ( Figure 1C). Upon incubation at 25 • C for 4 days, in darkness, the colony morphology of the SZ-2-3y strain was abnormal, mycelia were sparse, and the diameter of colonies was smaller, with an average growth rate of 0.8 cm/day compared with 1.9 cm/day for B05.10. Thus, temperature had a significant effect on SZ-2-3y morphology.
Pathogenicity tests were also conducted to investigate the virulence of SZ-2-3y. No lesions were observed on plants that were inoculated with the SZ-2-3y strain and grown at room temperature after 4 dpi, regardless of the inoculation method (intact inoculation for P. polyphylla, tobacco leaves and tomato fruits, or wound inoculation for apple fruits). By contrast, B05.10 caused large necrotic lesions on the leaves of P. polyphylla and tobacco, and on the fruits of tomato and apple, with average lesion diameters of 0.86, 2.56, 2.54 and 3.28 cm, respectively ( Figure 1D,E). These results reveal that the SZ-2-3y strain was hypovirulent.

Identification of Four Mycoviruses in the SZ-2-3y Strain
Since the SZ-2-3y strain exhibited hypovirulence and apparent phenotypic changes under the influence of temperature, we questioned whether it could infect mycoviruses, and if mycoviruses were responsible for the hypovirulence of the SZ-2-3y strain. We first extracted dsRNAs from SZ-2-3y mycelia, and obtained bands corresponding to fragments~8 kb and 3 kb in length following SI nuclease (Thermo Fisher Scientific, Waltham, MA, USA) digestion or long-time storage, to trigger rRNA degradation (Supplementary Figure S2). To identify the sequences amplified from the SZ-2-3y strain, high-throughput sequencing was carried out. The obtained high-quality reads were further assembled, and obtained contigs were aligned using BLASTx in NCBI. To clarify whether SZ-2-3y carries mycoviruses, RT-PCR was performed using total RNAs or dsRNA extracted from the SZ-2-3y strain as template, in combination with specific primers designed from the obtained contigs. RT-PCR results and sequencing analysis resulted in 587 nt, 578 nt, 523 nt and 260 nt fragments, in accordance with the four mycovirus contigs obtained by high-throughput sequencing, which were confirmed to be derived from the SZ-2-3y strain (Supplementary Table S1). We obtained six mycoviral contigs (Contig111, Contig334, Contig420, Contig25026, Con-tig12923 and Contig46148; Supplementary File S1). A homology search using BLASTx showed that RdRps of the four novel mycoviruses were most closely related to the homolog in Fusarium boothii mitovirus 1 (FbMV1; 54%), Sclerotinia sclerotiorum ourmia-like virus 13 (79%), Botrytis cinerea ourmia-like virus 12 (79%) and Botrytis cinerea fusarivirus 5 (85%), as shown in Table 1. Since nine mitoviruses in BcMV1-9 have been identified in B. cinerea and grape, the mitovirus found in the SZ-2-3y strain was designated as Botrytis cinerea mitovirus 10 (BcMV10). Meanwhile, 17 Botrytis cinerea ourmia-like viruses have been identified in B. cinerea and reported in the NCBI database and the literature. The botouliviruses identified in SZ-2-3y were designated Botrytis cinerea botoulivirus 18 (BcBoV18) and Botrytis cinerea botoulivirus 19 (BcBoV19), based on ICTV classification (https://talk.ictvonline.org/taxonomy/, accessed on 28 July 2021). Pathogenicity tests were also conducted to investigate the virulence of SZ-2-3y. No lesions were observed on plants that were inoculated with the SZ-2-3y strain and grown at room temperature after 4 dpi, regardless of the inoculation method (intact inoculation for P. polyphylla, tobacco leaves and tomato fruits, or wound inoculation for apple fruits). By contrast, B05.10 caused large necrotic lesions on the leaves of P. polyphylla and tobacco, and on the fruits of tomato and apple, with average lesion diameters of 0.86, 2.56, 2.54 and 3.28 cm, respectively ( Figure 1D,E). These results reveal that the SZ-2-3y strain was hypovirulent. Results are means ± standard deviation (SD; n = 7); (C) histogram showing the number of sclerotia for the two isolates (** p < 0.01); (D) pathogenicity of B. cinerea isolates on leaves of Paris polyphylla and tobacco, fruits of tomato following intact inoculation with mycelia, and apple fruits following wound inoculation with the mycelia using isolates of SZ-2-3y and B05.10, respectively (20 • C, 4 dpi); (E) histogram showing average lesion diameters of leaves of P. polyphylla and tobacco, and fruits of tomato and apple, caused by the two isolates. B05.10 served as a reference strain in the above experiments.  Figure 2).

Analysis of Full-Length cDNA Genome Sequences: A Novel Mitovirus and Two New Botouliviruses
The sequences of full-length cDNAs were determined by assembling partial cDNAs amplified from the purified dsRNAs or tRNAs as templates, using RT-PCR with random primers and cDNA end amplification with adaptor primers. RT-PCR, combined with rapid amplification of cDNA ends (RACE) cloning and sequencing, yielded intermediate sequences of three viruses (2452 nt for BcMV10, 1899 and 935 nt for BcBoV18, and 1731 and 1342 nt for BcBoV19) from specific designed primers. The obtained 5′-and 3′-RACE sequences were 452 and 301 nt for BcMV10, 202 and 287 nt for BcBoV18, and 387 and 192 nt for BcBoV19 (Supplementary Table S1). The obtained full-length genome sequences have been deposited in GenBank under accession number OK634394 for BcMV10, OK634395 for BcBoV18 and OK634396 for BcBoV19 (Supplementary File S3).

A Novel Mitovirus
The genomic structure of the mitovirus is shown in Figure 3A. The full-length genome sequence of BcMV10 is 2945 bp with an AU-rich element, and it contains a single ORF, spanning positions 361 to 2817 nt, which is predicted to encode an RdRp protein of 818 aa. The 5′-untranslated region (UTR) and 3′-UTR are 360 nt and 128 nt, respectively. The highest RdRp aa sequence identity of BcMV10 is 54%, shared with FbMV1; this is

Analysis of Full-Length cDNA Genome Sequences: A Novel Mitovirus and Two New Botouliviruses
The sequences of full-length cDNAs were determined by assembling partial cDNAs amplified from the purified dsRNAs or tRNAs as templates, using RT-PCR with random primers and cDNA end amplification with adaptor primers. RT-PCR, combined with rapid amplification of cDNA ends (RACE) cloning and sequencing, yielded intermediate sequences of three viruses (2452 nt for BcMV10, 1899 and 935 nt for BcBoV18, and 1731 and 1342 nt for BcBoV19) from specific designed primers. The obtained 5 -and 3 -RACE sequences were 452 and 301 nt for BcMV10, 202 and 287 nt for BcBoV18, and 387 and 192 nt for BcBoV19 (Supplementary Table S1). The obtained full-length genome sequences have been deposited in GenBank under accession number OK634394 for BcMV10, OK634395 for BcBoV18 and OK634396 for BcBoV19 (Supplementary File S3).

A Novel Mitovirus
The genomic structure of the mitovirus is shown in Figure 3A. The full-length genome sequence of BcMV10 is 2945 bp with an AU-rich element, and it contains a single ORF, span- ning positions 361 to 2817 nt, which is predicted to encode an RdRp protein of 818 aa. The 5 -untranslated region (UTR) and 3 -UTR are 360 nt and 128 nt, respectively. The highest RdRp aa sequence identity of BcMV10 is 54%, shared with FbMV1; this is followed by other mitoviruses in fungi hosts, including Nigrosp oraoryzae mitovirus 1 (NoMV1; 50%); other plant mitoviruses, including Solanum chacoense mitovirus 1 (19%, 115/591 RdRp aa); and plant mitochondria, including Solanum tuberosum (24%, 74/310 RdRp aa) (Supplementary Table S3). Both the 5 -and 3 -UTRs of the BcMV10 genome are clearly folded into two potential stem-loop structures. In addition, the long panhandle structure is attached to the 3 -and 5 -termini sequences through base pairing (Supplementary Figure S3A).

Two Novel Botouliviruses
The full-length genome sequences of BcBoV18 and BcBoV19 are 2759 nt and 2812, respectively, and they possess a single ORF encoding RdRp of 1986 nt and 2085 nt, respectively ( Figure 3B,C). The genomes of BcBoV18 and BcBoV19 contain AUs with 51.3% and 56% identity, respectively. The RdRp encoded by the ORF of BcBoV18 is most closely related to that of BcOlV17, while that of BcBoV19 is most closely related to that of Botrytis cinerea ourmia-like virus 12 (BcOlV12), with RdRp aa identities of 78%, followed by members of the Botoulivirus genus (Supplementary Tables S4 and S5). The 5′-UTR and 3′-UTR are 43 nt and 730 nt for BcBoV18, and 54 nt and 673 nt for BcBoV19. The 5′-and 3′-UTRs of BcBoV18 and BcBoV19 genomes could form potential terminal stable stem-loop structures (Supplementary Figure S3B,C).

Molecular Phylogenetic Analysis of Viral RdRps from BcMV10, BcBoV18, BcBoV19 and Their Relatives
A phylogenetic tree and multiple alignment analyses of the RdRp sequences were performed on the three newly discovered mycoviruses, 25 representative members belonging to Mitoviridae, and 24 representative members in Botourmiaviridae from the NCBI database.

A Novel Mitovirus
We searched for the presence of BcMV10-like sequences in the NCBI database using the putative BcMV10-encoded protein as the query, revealing several significant matches,

Two Novel Botouliviruses
The full-length genome sequences of BcBoV18 and BcBoV19 are 2759 nt and 2812, respectively, and they possess a single ORF encoding RdRp of 1986 nt and 2085 nt, respectively ( Figure 3B,C). The genomes of BcBoV18 and BcBoV19 contain AUs with 51.3% and 56% identity, respectively. The RdRp encoded by the ORF of BcBoV18 is most closely related to that of BcOlV17, while that of BcBoV19 is most closely related to that of Botrytis cinerea ourmia-like virus 12 (BcOlV12), with RdRp aa identities of 78%, followed by members of the Botoulivirus genus (Supplementary Tables S4 and S5). The 5 -UTR and 3 -UTR are 43 nt and 730 nt for BcBoV18, and 54 nt and 673 nt for BcBoV19. The 5 -and 3 -UTRs of BcBoV18 and BcBoV19 genomes could form potential terminal stable stem-loop structures (Supplementary Figure S3B,C).

Molecular Phylogenetic Analysis of Viral RdRps from BcMV10, BcBoV18, BcBoV19 and Their Relatives
A phylogenetic tree and multiple alignment analyses of the RdRp sequences were performed on the three newly discovered mycoviruses, 25 representative members belonging to Mitoviridae, and 24 representative members in Botourmiaviridae from the NCBI database.

A Novel Mitovirus
We searched for the presence of BcMV10-like sequences in the NCBI database using the putative BcMV10-encoded protein as the query, revealing several significant matches, mainly in fungi associated with plants (Supplementary Table S3). The phylogenetic tree indicated two subgroups supported by high bootstrap values. One subgroup includes fungi mitoviruses consisting of three clades (Clades I, II and III); another subgroup contains plant mitoviruses such as Solanum chacoense mitovirus 1, Chenopodium quinoa mitovirus 1 and Cannabis sativa mitovirus 1, as well as plant mitochondria such as Anthurium amnicola, Solanum chacoense, Solanum tuberosum and Eucalyptus grandis, which clustered into one large clade belonging to the proposed Mitoviridae family ( Figure 4A). BcMV10 is part of a separate branch clustered with fungi mitoviruses forming clade I. Mitoviruses from B. cinerea strains are widely distributed in Clades I, II and III. Based on Conserved Domain (CD) search results and multiple-sequence alignment analysis, BcMV10 RdRps contain typical conserved motifs (I-VI) which are characteristic of the mitoviral RdRp domain superfamily (pfam05919), similar to fungi mitoviruses, plant mitoviruses, and plant mitochondria to some extent ( Figure 4B).

Two Novel Botouliviruses
Phylogenetic tree analysis separated BcBoV18, BcBoV19 and related members into six groups in accordance with the genera Botoulivirus, Penoulivirus, Scleroulivirus, Magoulivirus, Rhizoulivirus and Ourmiavirus in the Botourmiaviridae family. The sequence identities of BcBoV18 and BcBoV19 with members in the Botourmiaviridae family are in accordance with the phylogenetic tree results (Supplementary Tables S4 and S5). BcBoV18 and BcBoV19 are divided into different branches in the Botoulivirus genus, which share 40% of their identity with RdRp aa sequences. Taken separately, BcBoV18 with BcOLV17 share high similarity, while BcBoV19, BcOLV12 and Sclerotinia minor botoulivirus 1 (SmBV1) are more closely related. However, mycoviruses from B. cinerea strains were divided into the genera Penoulivirus and Magoulivirus, except for Botoulivirus, revealing the presence of diverse botourmiaviruses in the B. cinerea natural population ( Figure 5A). Based on a CD search of the NCBI database, putative conserved domains were detected in BcBoV18, while no putative conserved domains were detected in BcBoV19. Nevertheless, we tried to predict putative conserved domains in BcBoV19 using the MOTIF search tool (https://www. genome.jp/tools/motif, accessed on 28 July 2021). Multiple alignments of putative RdRp regions showed that BcBoV18, BcBoV19 and members of the family Botourmiaviridae contain eight conserved motifs (I to VIII), including a highly conserved GDD signature ( Figure 5B). Therefore, BcBoV18 and BcBoV19 are two novel members of the Botoulivirus genus of the Botourmiaviridae family. Interestingly, nuclear localisation signals (NLSs) in RdRps of BcBoV18 and BcBoV19 were predicated using the online cNLS Mapper (Supplementary  Table S7).
which clustered into one large clade belonging to the proposed Mitoviridae family ( Figure  4A). BcMV10 is part of a separate branch clustered with fungi mitoviruses forming clade I. Mitoviruses from B. cinerea strains are widely distributed in Clades I, II and III. Based on Conserved Domain (CD) search results and multiple-sequence alignment analysis, BcMV10 RdRps contain typical conserved motifs (I-VI) which are characteristic of the mitoviral RdRp domain superfamily (pfam05919), similar to fungi mitoviruses, plant mitoviruses, and plant mitochondria to some extent ( Figure 4B).

Distribution of Mycovirus-Derived RNA Reads along the Viral Genome
The RNA sequencing (RNA-seq) library was constructed and sequenced from the mixture samples. Mapping of RNA reads obtained from the RNA-seq library to the fulllength genomes showed that 1,411,244 (3.42% overall alignment rate), 143,869 (0.035%) and 131,690 (0.32%) reads aligned with the obtained BcMV10, BcBoV18 and BcBoV19 genomic sequences, respectively. The RNA reads almost covered the full-length BcMV10, BcBoV18 and BcBoV19 RNAs. Much greater RNA read depth was achieved for the 3′-UTRs in the full-length genomes of BcBoV18 and BcBoV19. By contrast, no apparent read peaks were observed throughout the whole BcMV10 genome (Figure 6).

Distribution of Mycovirus-Derived RNA Reads along the Viral Genome
The RNA sequencing (RNA-seq) library was constructed and sequenced from the mixture samples. Mapping of RNA reads obtained from the RNA-seq library to the fulllength genomes showed that 1,411,244 (3.42% overall alignment rate), 143,869 (0.035%) and 131,690 (0.32%) reads aligned with the obtained BcMV10, BcBoV18 and BcBoV19 genomic sequences, respectively. The RNA reads almost covered the full-length BcMV10, BcBoV18 and BcBoV19 RNAs. Much greater RNA read depth was achieved for the 3 -UTRs in the full-length genomes of BcBoV18 and BcBoV19. By contrast, no apparent read peaks were observed throughout the whole BcMV10 genome ( Figure 6).

BcMV10 Sequences Are Present in Cucumber Plants Inoculated with Strain SZ-2-3y Mycelia
We inoculated SZ-2-3y mycelium agar onto cotyledons of cucumber plants. No SZ-2-3y mycelia were isolated from either the leaves or roots of cucumber plants (the data is

BcMV10 Sequences Are Present in Cucumber Plants Inoculated with Strain SZ-2-3y Mycelia
We inoculated SZ-2-3y mycelium agar onto cotyledons of cucumber plants. No SZ-2-3y mycelia were isolated from either the leaves or roots of cucumber plants (the data is not shown). Following RNA extraction from leaves or roots, nested RT-PCR, using different sets of BcMV10-based primers, yielded amplicons of expected sizes, suggesting that the BcMV10 sequence was present in cucumber plants. By contrast, no BcMV10-specific amplicons were obtained from healthy leaf tRNAs ( Figure 8A, Supplementary Figure S4A). Associated experiments were conducted and (nested) RT-PCR analysis consistently confirmed the presence of the BcMV10 sequence in systemic leaves of 'Cuiyu' at 2 weeks post-inoculation. The average infection efficiency rate was~50% (1/2) in the local leaves; 50% (2/4) and 33% (3/9) in the new leaves; and 2/4 (50%) in the roots. Meanwhile, it was positive for 10% (1/9) of the new leaves from inoculated 'Suyo' plants ( Table 3). The presence of the BcMV10 sequence was also verified by Sanger sequencing of the above, randomly selected RT-PCR amplicons.   The full-length BcMV10 cDNA clone and pXT1 were incorporated by homologous recombination [54]. Briefly, total RNA from the SZ-2-3y isolate was extracted and reversetranscribed into cDNA to obtain the template. The full-length BcMV10 genome of 2975 nt was obtained and purified, and inserted into the linearised pXT1 expression vector. They were assembled to generate BcMV10 by transformation of E. coli and Agrobacterium cells. PCR detected 560 bp, 825 bp and 3522 bp target fragments that were used to screen positive clones from the above plasmids, which were designated pXT1/BcMV10 (Supplementary  Table S6, Figure 8B). All positive clones were verified by DNA sequencing.
Cucumis sativus 'Suyo' and 'Cuiyu' plants were agroinoculated into the abaxial intercellular leaf space of two cotyledons with infectious clones of pXT1/BcMV10. An empty pXT1 vector served as a negative control ( Figure 8B, Supplementary Figure S4B). At nearly  (Table 3). These results suggest that the BcMV10 sequence was present in the new leaves of inoculated cucumber shoots.

Discussion
Currently, grey mold disease caused by Botrytis species threatens the horticultural crop, the P. polyphylla plant, by diminishing both its economic and medical value. Botrytis species can infect a large variety of hosts, including field-crop plants, fruit trees, vegetables and flowers [27,30,39,40]. Grey mold disease is a serious problem for agriculture worldwide, and current control relies on the excessive use of fungicides [27][28][29][30]40]. Unfortunately, cultivars/varieties that possess a high resistance to Botrytis species infection are not yet available. Nevertheless, mycovirus resources might be potential biocontrol agents for the control of fungal disease, because various mycoviruses have been identified and successfully used for chestnut blight disease and sclerotium disease [25,[58][59][60][61][62].
In the present study, four mycoviruses (BcMV10, BcBoV18, BcBoV19 and BcFV8) were identified in the B. cinerea SZ-2-3y strain, the causal agent of P. polyphylla grey mold disease ( Figure 1, Table 1). Based on genome structure, sequence homology, and phylogenetic analysis, all four are novel mycoviruses belonging to Mitovirus, Botoulivirus and Fusarivirus genera, respectively (Figures 2-5). However, RdRp aa reveals pairwise identities of BcMV10, with the six closest mitoviruses clustering in the same branch ranging from 54% to 32% (Figure 4, Supplementary Table S3). These results revealed that mitoviruses are widely distributed and ancestrally present in fungal hosts.
Mycoviruses infecting a plant pathogenic Botrytis species are classified in different genera of Botourmiaviridae. Herein, BcBoV18 and BcBoV19 are clustered in two separate branches in the genus Botoulivirus ( Figure 5, Supplementary Tables S4 and S5). The BcMV10 reads had no apparent mapping depth along the whole genome, but BcBoV18 and BcBoV19 had high mapping depth to RNA reads in the 3 -termini of the full-length genome ( Figure 6). The distinguishing distribution profiles could also be implicated in the different viral expression strategies for Botouliviruses and Mitoviruses. In conclusion, we characterised two novel botouliviruses and a novel mitovirus from B. cinerea strain SZ-2-3y, and provide insights into their molecular diversity and evolutionary origins.
It is well known that co-infections can result in interactions between viruses and hosts, displayed as antagonism or synergism [63][64][65][66]. It will be interesting to see whether interactions exist between the co-infected mitovirus, fusarivirus and the two botouliviruses. Confrontation cultures of SZ-2-3y and B05.10 showed an absence of BcMV10 but not of BcBoV19, and the mixture of BcBoV18, BcBoV19 and BcFV8 did not affect the phenotype and virulence of B05.10 on tobacco (Figure 7, Table 2). This result indicates that BcBoV18, BcBoV19 and BcFV8 might transmit horizontally, and have no apparent influence on the growth and pathogenicity of B. cinerea strains. Thus, we speculate that BcMV10 might be related to the hypovirulent SZ-2-3y strain. It is reported that BcMV1, SsMV1/HC025, SsMV3/NZ1 and SsMV4/NZ1 are hypovirulent, but their effects on fungal phenotype and biological characterisation are poorly understood due to the difficulty in isolating and transmitting to virus-free strains, as well as restricting the host antiviral RNA silencing machinery [6,15,17,22]. To our knowledge, BcMV10 is the tenth reported mitovirus infecting B. cinerea, and it is predicted to be related to the hypovirulence of SZ-2-3y; however, this needs to be further explored. In the future, we will attempt to obtain virus-free and BcMV10 strains using protoplast regeneration, which could shed light on the direct impact of BcMV10 on host biology and pathogenicity.
Bioinformatic analyses showed that related ourmia-like virus and mitovirus sequences are present in different organisms, suggesting a wide prevalence of these types of viruses across kingdoms [23,67,68]. Mitovirus sequences may transfer between fungi and plants, likely via fungus-mediated HGT events [19,20]. However, there is no direct evidence that mitovirus sequences can be transferred between fungi and plants. In the present study, we used the BcMV10-encoded protein in BLASTx searches for the presence of BcMV10-like sequences in the TSA-nr database. We identified several significant matches, with identity ranging from 26% to 37% in Pseudo-nitzschia australis, Solanum chacoense and Anthurium Amnicola (Supplementary Table S3). Viral genome organisation and phylogenetic analyses of BcMV10 indicate that this virus is closely related to fungal mitoviruses. Together with plant mitoviruses and plant mitochondria, they form a monophyletic cluster and belong to the proposed Mitoviridae family (Figure 4), consistent with their descent from a common ancestor [17,19]. We wondered whether the BcMV10 genome sequence is transmitted to host plants via HGT. Therefore, we tried to artificially inoculate cucumber plants with BcMV10 in different ways. Interestingly, partial BcMV10 genome sequences were detected in local and new leaves, and even the roots of cucumber plants inoculated with strain SZ-2-3y mycelia or Agrobacteria harbouring pXT1/BcMV10, respectively ( Figure 8, Table 3, Supplementary Figure S4). In contrast, BcBoV18, BcBoV19 and BcFV8 sequences were not detected in plants inoculated with strain SZ-2-3y mycelia. This reveals that the B. cinerea mitoviral sequence may exhibit transboundary entry into plants. Nevertheless, how mitoviral sequences are transmitted into plants, and whether they may be related to genetic material exchange, needs to be explored in future research.
Recently, cross-kingdom existence in different organisms has been explored. Some experimental evidence indicates that some plant viruses can replicate in fungus cells. Cucumber mosaic virus (CMV) may infect the phytopathogenic fungus Rhizoctonia solani in nature [69]. Additionally, plant viroid transcripts may infect phytopathogenic fungi based on artificial plant viroid infectious clones [70]. It is known that virus and virus-like organisms transfer between plants and fungi. An insect can transmit mycoviruses, which could provide opportunities for horizontal virus transfer between different organisms [71]. Thus, infectious mycovirus CHV1 may infect plants with the help of tobacco mosaic virus (TMV) and movement protein (MP), suggesting that mycovirus may horizontally transfer between fungi and plants [72]. Our results demonstrate the entry of a partial B. cinerea mitovirus genome sequence into cucumber plants, and suggest that BcMV10 could be closely linked to the horizontal transmission of mitoviruses between plants and fungi. We tried to inoculate B05.10 mycelia onto the new leaves of cucumber positive for BcMV10. The new leaves exhibited symptoms at 2 days after inoculation. We further isolated B05.10 mycelia from diseased leaves, which were negative for BcMV10 (Supplementary Figure S5). However, it still not known how mitoviral sequences transmit between fungi and plants, and if mitoviruses can infect plants and exert specific biological functions. We believe that the discovery of BcMV10 significantly contributes to our understanding of the complex relationship between fungi, mitoviruses and plants.

Conclusions
Herein, we describe Botrytis cinerea strain SZ-2-3y, the causal agent of P. polyphylla grey mold disease in China. A novel mitovirus, two novel botouliviruses and a novel fusarivirus, isolated from the hypovirulent B. cinerea SZ-2-3y strain, were identified and molecularly characterised. Based on phylogenetic tree analysis, mitoviruses infecting B. cinerea species cluster into different clades. Conversely, some mitoviruses isolated from different fungal species share relatively close relationships. These results imply that mitoviruses may have evolved before fungal hosts diverged. In addition, alignment of RdRp aa sequences showed that BcMV10 is closely related to members of fungal mitoviruses. Together with plant mitoviruses and plant mitochondria, they form a monophyletic cluster and belong to the proposed Mitoviridae family, revealing their descent from a common ancestor. The BcMV10 genome sequence is present in tissues of cucumber plants experimentally inoculated with strain SZ-2-3y mycelia and agrobacterium-mediated pXT1/BcMV10. This provides evidence and new insight into the occurrence of cross-kingdom transfer of mitovirus sequences from fungi to plants, and illuminates aspects of mycoviral diversity and evolution. Whether or not BcMV10 plays a significant role in the relationship between fungi, mitoviruses and plants remain to be explored.

Supplementary Materials:
The following are available online at https://www.mdpi.com/article/ 10.3390/v14010151/s1: Table S1-list of the oligonucleotide primers used in this study; Table S2sequence identities of BcFV8 to other viruses used for multiple sequence alignment and phylogenetic analysis; Table S3-sequence identities of BcMV10 to other viruses used for multiple sequence alignment and phylogenetic analysis; Table S4-sequence identities of BcBoV18 to other viruses used for multiple sequence alignment and phylogenetic analysis; Table S5-sequence identities of BcBoV19 to other viruses used for multiple sequence alignment and phylogenetic analysis; Table S6-primers used for construction of pXT1/BcMV10 clone and detection of BcMV10 sequences in cucumber plants; Table S7-predicated NLSs in mycoviruses RdRp using the online website (http://nls-mapper. iab.keio.ac.jp, accessed on 20 September 2021); Figure S1-the symptoms of Paris polyphylla plants infected with grey mold diseases as a whole (A) and the grey mold symptoms on single P. polyphylla plant (B); Figure S2-1.2% agarose electrophoreses analysis of dsRNAs; Figure S3-the stem-loop second structures of terminal, non-coding RNA regions of BcMV10 (A), BcBoV18 (B) and BcBoV19 (C), respectively; Figure S4-the BcMV10 sequences from the new leaves of 'Cuiyu' cv inoculated with SZ-2-3y mycelia (A), and agro-infiltration with pXT1/BcMV10 clone (B) were detected by RT-PCR, respectively. Line 1-9: total RNAs of the leaves from 9 'Cuiyu' cv plants, respectively, were extracted as templates to detect BdMV10; Ck2-: total RNA from 'Cuiyu' cv plant inoculated by empty PDA as negative control; Ck1-: pXT1-inoculated plants as negative control; Ck+: Total RNA from SZ-2-3y mycelia on PDA; M: Marker III; Figure S5-schematic diagram of intact inoculation of B05.10 strain on the new leaves of cucumber plants and RT-PCR detection of the BcMV10 sequence from the isolated strains. The symptoms were exhibited on new leaves of 'Cuiyu' cucumbers at 2 days (A) and 7 days (B) post-inoculation of B05.10 strain; the colony morphology of strains isolated from the leaves at two weeks post-inoculation of B05.10 strain (1-5, the isolated strains) (C); BcMV10 sequences from 5 isolated strains were detected by RT-PCR (D). Supplementary File S1-the sequences information of six viral contigs from Botrytis cinerea SZ-2-3y strain used in this study; Supplementary File S2-the sequences of partial genome and encoded RdRp aa from BcFV8; Supplementary File S3the sequences of full-length genome and encoded complete RdRp aa from BcMV10, BcBoV18 and BcBoV19, respectively.