key: cord-0070694-xsaj9047
authors: Shafik, Karim; Umer, Muhammad; You, Huafeng; Aboushedida, Hamdy; Wang, Zhenhua; Ni, Dejiang; Xu, Wenxing
title: Characterization of a Novel Mitovirus Infecting Melanconiella theae Isolated From Tea Plants
date: 2021-11-17
journal: Front Microbiol
DOI: 10.3389/fmicb.2021.757556
sha: 26667b77a23b7e5a8adbd60ee29274385363a482
doc_id: 70694
cord_uid: xsaj9047

A dsRNA segment was identified in the fungus Melanconiella theae isolated from tea plants. The complete dsRNA sequence, determined by random cloning together with RACE protocol, is 2,461 bp in length with an AU-rich content (62.37%) and comprises a single ORF of 2,265-nucleotides encoding an RNA-dependent RNA-polymerase (RdRp, 754 amino acids in size). The terminus sequences can fold into predicted stable stem-loop structures. A BLASTX and phylogenetic analysis revealed the dsRNA genome shows similarities with the RdRp sequences of mitoviruses, with the highest identity of 48% with those of grapevine-associated mitovirus 20 and Colletotrichum fructicola mitovirus 1. Our results reveal a novel member, tentatively named Melanconiella theae mitovirus 1 (MtMV1), belongs to the family Mitoviridae. MtMV1 is capsidless as examined by transmission electron microscope, efficiently transmitted through conidia as 100 conidium-generated colonies were analyzed, and easily eliminated by hyphal tipping method combined with green-leaf tea powder. MtMV1 has a genomic sequence obviously divergent from those of most members in the family Mitoviridae and some unique characteristics unreported in known members. This is the first report of a mycovirus infecting Melanconiella fungi to date.

Mycoviruses (also known as fungal viruses) are widespread in all major taxa of fungi (Strauss et al., 2000; Lim et al., 2005; Ro et al., 2007; Stielow et al., 2011; Magae, 2012) . Mycoviruses are usually associated with symptomless (cryptic or latent symptoms) infections (Kamaruzzaman et al., 2019) . However, some can attenuate the virulence of their fungal hosts and contribute as biocontrol agents for preventing the diseases induced by phytopathogenic fungi (Wu et al., 2007; Ghabrial and Suzuki, 2009; Pearson et al., 2009; Abid et al., 2018; Mizutani et al., 2018) . According to the International Committee on Taxonomy of Virus 1 , mycoviruses have been taxonomically grouped into 19 families, including six with dsRNA genomes (Chrysoviridae, Partitiviridae, Quadriviridae, Reoviridae, Totiviridae, and Megabirnaviridae), nine with positive-sense single-stranded RNA (+ ssRNA) (Alphaflexiviridae, Barnaviridae, Botourmiaviridae, Deltafleiviridae, Endornaviridae Gammaflexiviridae, Hypoviridae, Mitoviridae, and Narnaviridae) (Kamaruzzaman et al., 2019) , one with negative-sense ssRNA (-ssRNA; Mymonaviridae) (Liu et al., 2014) , one with ssDNA (Genomoviridae) (Yu et al., 2010) , two with reverse-transcribing RNA (Metaviridae and Pseudoviridae), and unclassified taxa (Liu et al., 2014; He et al., 2018) . The number of mycovirus primarily increased in recent years through deep sequencing analysis, and most of the mycoviruses are derived from major phytopathogenic fungi isolated from crops, fruits, and forest trees. In contrast, to date, there are much less reports on mycoviruses associated with fungi isolated from tea plants [Camellia sinensis (L.) O. Kuntze], which originated from China and has been cultivated for over 3,000 years.

Mitoviruses (family Mitoviridae) represent the simplest group of all RNA viruses (Hillman and Cai, 2013) , with a unique genome of + ssRNA which ranges in size between 2.0 and 4.5 kb, encompassing a single long open reading frame (ORF) which encodes a putative RNA-dependent RNA polymerase (RdRp) with six conserved amino acid motifs (I-IV) (Osaki et al., 2015; Ran et al., 2016; Chen et al., 2017; Wang et al., 2019; Kocanová et al., 2020; Liu et al., 2021) . Their Genomes are AU rich, especially in the third position of codons (Hong et al., 1998; Xie and Ghabrial, 2012; Liu et al., 2019; Torres-Trenas and Pérez-Artés, 2020; Wang et al., 2021) . Mitoviral genomes are characterized by multiple UGA codons encode tryptophan (Trp) instead of stop translation (Nibert, 2017; Nibert et al., 2019) ; however, the termination codon is usually UAA and sometimes UAG (Nibert, 2017) . The terminus sequences could be folded into stable genomic secondary structures that have crucial roles in replicating and protecting their genome (Khalifa and Pearson, 2013; Sahin and Akata, 2019; Torres-Trenas and Pérez-Artés, 2020; Liu et al., 2021) . However, mitoviruses are initially found infecting fungi, they have recently been reported from plants as expected infectious entities (Nibert et al., 2018; Nerva et al., 2019; Vong et al., 2019) . They are transmitted through horizontal transmission (hyphal anastomosis and heterokaryosis) or vertical transmission (spores) (Romeralo et al., 2018; Kamaruzzaman et al., 2019) . Mitoviruses are the only viruses associated with the mitochondria in their hosts (Cole et al., 2000; Osaki et al., 2005; Nibert, 2017; Fonseca et al., 2021) and often led to hypovirulence in their fungal hosts (Nuss, 2005; Xu et al., 2015; Yang et al., 2018) .

The genus Melanconiella (order Diaporthales, ascomycetes) involves saprophytes and endophytes as well as phytopathogenic fungi causing cankers, diebacks, and rots on crops, ornamentals, and forest trees (Du et al., 2017) . Members of the fungus Melanconiella were mostly observed to be highly host-specific pathogenic and endophytic fungi of the host family Betulaceae especially in the north temperate zone (Voglmayr et al., 2012) . To our knowledge, no mycoviruses have been reported from Melanconiella theae. Here, we characterized a novel mitovirus from M. theae that represent the first mycovirus from Melanconiella fungi. It will expand our understanding of the origin, ecology, and evolutionary pathways of mycoviruses (Ghabrial, 1998; Ghabrial and Suzuki, 2009 ).

Melanconiella theae strains WJB-5 (Supplementary Figure 1) , WJB-1-1, WJB-1-2, WJB-1-9, WJB-1-12, WJB-1-20, WJT-1-1, WJT-1-3, WJT-1-7, WJT-1-9, WJT-1-12, and WJT-1-13 were isolated from naturally infected tea plants showing necrotic symptoms (Supplementary Figure 1A) collected from Xuanen county, Hubei province, China, and identified based on the internal transcribed spacer (ITS) region. Of these strains, WJB-1-9 and WJB-1-12 have the same genetic background with WJT-1-9 and WJT-1-12, respectively, since they were derived from the same strain by culturing single mycelium. Fungal mycelia of these strains were grown on PDA for a week at 25 • C and then stored in a 25% (v/v) glycerol solution at −70 • C until use.

Viral dsRNA was extracted from the fungal mycelial mass using the silica spin column-based method as described previously . The dsRNAs extracted from M. theae strains WJB-1-1, WJB-1-2, WJB-1-9, WJB-1-12, WJB-1-20, and WJB-5, and those from strain WJB-5 were used for viral genome characterization and treated with RNase-free DNase I and S1 Nuclease to remove any DNA and ssRNA remains, fractionated by agarose gel (1%, w/v) electrophoresis, and detected by UV transillumination after staining with ethidium bromide (0.1 mg/mL).

Purified dsRNA from WJB-5 was used as a template for cDNA synthesis following the previously described method (Xie et al., 2011) using a cDNA synthesis kit (Promega, Madison, WI, United States) and tagged random dN6 primer, 05RACE-3RT (Supplementary Table 1 ). The resulting cDNA was PCR-amplified using 2 × Gimmico GimiTaq PCR Mix (Gimmico Biotech Co., Wuhan, China) and 05RACE-3 primer (Supplementary Table 1 ). RT-PCR products were purified and then ligated to the pMD18-T vector (TaKaRa, Beijing, China), then transformed into competent cells of Escherichia coli DH5α. Positive clones with cDNA inserts were selected on Luria-Bertani (LB) agar medium containing ampicillin (50 µg/mL), checked with PCR using M13F-47/M13R-48 primers (Supplementary Table 1) , and then submitted for sequencing. Specific primers, MtMV1-F and MtMV1-R (Supplementary Table 1) were designed for rapid amplification of cDNA ends (RACE) to amplify terminal sequences of the genome that were not cloned by the initial random cDNA synthesis. For RACE, the PC3-T7 loop adapter (Supplementary Table 1 ) was ligated to the 3 end of each strand in the purified dsRNA using T4 RNA ligase (TaKaRa, FIGURE 1 | Agarose gel electrophoresis of nucleic acids extracted from Melanconiella theae strains, enzyme digest and RT-PCR identification of MtMV1. (A) Nucleic acids extracted from strains WJB-1-1, WJB-1-2, WJB-1-9, WJB-1-12, WJB-1-20, and WJB-5. (B) The nucleic acids extracted from strain WJB-5 were untreated (Lane 1) or treated with S1 nuclease (Lane 2). (C) RT-PCR identification of MtMV1 in strain WJB-5. M: DL5000 DNA marker (100; 250; 500; 750; 1,000; 1,500; 2,000; 3,000; and 5,000 bp). The arrows indicate the target bands.

Beijing, China) at 16 • C for 18 h, as described previously (Ran et al., 2016) . The oligonucleotide-ligated dsRNA was purified, denatured in DMSO, and subsequently reverse transcribed. RACE and nested PCR amplification were performed using sets of MtMV1-F/PC2 and PC2/MtMV1-R primers (Supplementary Table 1 ) to amplify the 3 -and 5 -termini respectively. PCR amplicons were purified, cloned into a pMD18-T vector, and sequenced. All amplicons were sequenced by Tsingke Biological Technology Co., LTD., Wuhan, China. Sequence homology searches were performed with the BLASTX search tool on the National Center for Biotechnology Information, NCBI 2 . Sequence assemblies from the cDNA and the RACE clones were manipulated using DNAMAN software (version 7.0.2; Lynnon Corp., Vaudreuil-Dorion, QC, Canada) to obtain the complete genome sequence of the dsRNA extracted from the strain WJB-5.

RT-PCR amplification was performed using a specific primer pair (MtMV1-F1: 5 -GTCTTATCCGTGTATTCTGGG-3 ; MtMV1-R1: 5 -GGCTTGAGGAACATTGAGA-3 ), generating a 1981-bp fragment. An annealing temperature of 56 • C was used with a PCR Thermal Cycler (Model A200, LongGene, China).

The viral ORF was identified using Expasy 3 website, and the amino acid (aa) sequences encoded by the ORF were deduced. 2 http://www.ncbi.nlm.nih.gov/ 3 https://web.expasy.org/translate/ The yeast mitochondrial genetic code was used for the translation of the ORF. The mycovirus sequence was searched using the NCBI BLASTX program to identify similar mycoviruses for phylogenetic analysis. Maximum-likelihood (ML) phylogenetic tree with 1,000 bootstrap replicates was constructed using DNAMAN software (version 7.0.2; Lynnon Corp., Vaudreuil-Dorion, QC, Canada) for the amino acid sequences of RdRp from WJB-5 dsRNA along with 26 mitoviruses and 4 narnaviruses retrieved from the GenBank database. The evolutionary distances were computed using the Poisson correction method (uniform rates between sites). The potential RNA secondary structures of the 5 -and 3 -terminal sequences were predicted, and the free energy ( G) was estimated using CLC genomics workbench software (version 21.0.4; QIAGEN Aarhus, Aarhus, Denmark) as described previously (Liu and Di, 2020) . Mitoviral RdRp characteristic motifs were searched using the conserved domain database (CDD) search 4 provided by NCBI, and multiple amino acid alignments.

Mycelium plugs of the strain WJB-5 were grown on an autoclaved cellophane membrane placed on media for 8 days (d) at 25 • C. Minor modifications were adopted to isolate and purify viral Table 1 . The bootstrap values were deduced based on 1,000 replicates and indicated beside the branches. The evolutionary distances were computed using the Poisson correction method and the branch lengths correspond to the genetic distance; the scale bar refers to the genetic distance of 0.05. A tringle indicates MtMV1. particles in the previously described method (Ran et al., 2016) . In brief, 35 g mycelia of WJB-5 strain were harvested and crushed in the presence of liquid nitrogen. After that, the crushed mycelia were transferred into new sterile tubes containing 3 volumes of cold phosphate buffer (0.1 M sodium phosphate, pH 7.0, containing 0.2 M KCl and 0.5% β-mercaptoethanol) and gently mixed. The homogenate was centrifuged at 10,000 rpm and 4 • C for 15 min to remove the hyphal debris. After three ultracentrifugation cycles, virus particles were purified using a stepwise sucrose gradient of 20-50% (W/V). Following sucrose gradient ultracentrifugation, sucrose fractions were collected, diluted, and subjected to ultracentrifugation to pellet virus particles, and the pellet was resuspended again in 150 µL of phosphate buffer (0.1 mol/L sodium phosphate, pH 7.0). Viral dsRNA was recovered from purified viral particles by a phenolchloroform-based technique as described previously (Toni et al., 2018) , fractionated by agarose gel (1%, w/v) electrophoresis and detected by UV transillumination. For negative staining, 10 µL purified suspension of the virus was loaded on a hydrophobic parafilm surface, and subsequently copper grids with a carbonformvar coating (200-mesh) were floated onto them for 3 min. After drying with a filter paper, the grids plated with virus particles were immediately refloated for 3 min in a drop of 2% (W/V) phosphotungstic acid solution. The excess suspension was removed using filter paper, and the grids were then air-dried for a few min. Finally, transmission electron microscopy (TEM) was used to visualize the samples.

To eliminate the mycovirus from WJB-5 and five other (WJB-1-1, WJB-1-2, WJB-1-9, WJB-1-12, and WJB-1-20) strains of Melanconiella theae, hyphal tipping method together with thermal/dark treatment were conducted. All strains were separately incubated on PDA-T [potato dextrose agar medium containing 0.4% (w/v) green-leaf tea powder] (Yang and Zhang, 2019) media in the dark at 20 • C for 3 days. The hyphal tips in the colony margin were individually cut using a fine needle, transferred again to PDA-T in new Petri dishes (1 hyphal tip per dish), and incubated simultaneously. The colonies were subcultured on PDA-T at 20 • C in dark conditions for three consecutive subculturing. Finally, the subculturing was achieved on PDA through incubation at 20 • C and then used to detect MtMV1 dsRNA. To emphasize that the M. theae strains got cured from MtMV1, total RNA was extracted and RT-PCR amplification was performed using specific primers (MtMV1-F1: 5 -GTCTTATCCGTGTATTCTGGG-3 ; MtMV1-R1: 5 -GGCTTGAGGAACATTGAGA-3 ), as previously described.

Conidia of WJB-5 strain (harvested from a 30-d-old culture grown on PDA at 20 • C in the dark) were suspended in sterile distilled water to the final concentration of 3 × 10 7 conidia/mL. The resulting conidial suspension was poured on PDA plates and incubated at 20 • C in darkness, and used for culturing of singleconidium isolates. A total of 100 single-conidium isolates derived from WJB-5 were obtained and separately incubated on PDA plates at 20 • C for mycelial growth. After the colonies had grown, they were subjected to dsRNA extraction to detect the described virus infection.

Morphological traits of virus-infected strains (WJT-1-1, WJT-1-3, WJT-1-7, WJT-1-9, WJT-1-12, and WJT-1-13) and viruseliminated strains (WJB-5, WJB-1-1, WJB-1-2, WJB-1-9, WJB-1-12, and WJB-1-20) of M. theae were assessed as previously described (Jia et al., 2017) . Briefly, freshly grown mycelial discs (5 mm in diameter) were transferred onto PDA and incubated in darkness at 25 • C. Colony diameter was measured daily for 9 days to calculate growth rates (mm/d).

To check whether Melanconiella fungi harbor mycoviruses, six strains (WJB-1-1, WJB-1-2, WJB-1-9, WJB-1-12, WJB-1-20, and WJB-5) from M. theae isolated from Hubei province were randomly chosen from the fungal collection and subjected to dsRNA extraction. It revealed that all strains of the FIGURE 3 | Multiple alignments of the conserved motifs of RdRp in MtMV1 and other mitoviruses. Lines above the aligned sequences indicate to the positions of the motifs (I-VI). Numbers between every two motifs correspond to the number of aa residues separating the motifs. Identical, conserved, and semi-conserved aa residues are color-highlighted with homology level as 100% (black), 75% (magenta), 50% (turquoise), and 33% (yellow). phytopathogenic fungus M. theae contain mycovirus-like dsRNA segments (Figure 1 ). After digestion with DNase I and S1 nuclease enzymes, for the strain WJB-5, a clear band of dsRNA approximately 2.40 kb was detected (Figure 1 ).

The entire cDNA sequence of the dsRNA fragment of the strain WJB-5 was obtained through random priming cDNA synthesis, RT-PCR, and RACE cloning. Sequence assembly of thirty random cDNA clones revealed that the genome was composed of a single dsRNA segment with a molecular size of 2,461 bp and with an AU-rich content of 62.37%. Searching of the full-length genome and the deduced amino acid sequences using BLASTX and BLASTP, the coding protein shared 37.77-48.65% identities with the RdRp sequences of known mitoviruses belonging to the family Mitoviridae, and the highest identities are 48.65 and 48.35% with those of grapevine-associated mitovirus 20 (GaMV20) and Colletotrichum fructicola mitovirus 1 (CfMV1), respectively (Supplementary Table 2) . Thus, the dsRNA component is proposed as a novel mitovirus and 

Based on the yeast mitochondrial code, the genomic organization analysis revealed that the MtMV1 genome contains a single large ORF on the positive sense strand of the dsRNA (Figure 2A) .

The ORF starts and terminates at nucleotide positions 46 and 2,310, respectively, of 2,265 bp in size, and potentially codes for an RdRp protein with 754 aa in size. The ORF contains multiple codons (10 UGA and 2 UGG) encoding Trp, and the UAG codon stops the translation. Moreover, the 5 -and 3 -untranslated regions (UTRs) are 45 and 151 bp long, respectively ( Table 1) . The phylogenetic analysis separated the known mitoviruses into two distinct clusters and placed MtMV1 together with CfMV1 and GaMV20 in cluster I, which is distantly separated from other members (Figure 2B) . Based on the conserved domain database (CDD) search and multiple amino acid alignments, the six typical amino acid motifs (I-IV) characterized for the mitoviral RdRp domain superfamily were found to be conserved in MtMV1 (Figure 3) . However, MtMV1 contains many unique amino acids as compared with the other members in the conserved motifs (e.g., positions 12 and 20 in motif I, 4, and 19 in motif II, 18 in motif III, 3 and 4 in motif IV), which were composed by threonines (Thr) for MtMV1 while leucines (Leu), isoleucines (Ile) or valines (Val) for the others (Figure 3) . Moreover, an unique serine (Ser) in position 16 in motif IV and a unique histidine (His) in position 7 in motif V were observed for MtMV1, while alanine in motif IV and aspartic acids, asparagines or threonines in motif V in these positions for the others (Figure 3) .

The 5 -and 3 -terminus sequences of the MtMV1 genome were subjected to the prediction of their secondary structures using the CLC genomic workbench software. The results showed that both 5 -and 3 -UTRs could be folded into potential stem-loop structures that stabilize the RNA genome, with G values of −14.3 and −56.5 kcal/mol, respectively (Figures 4A,B) . Additionally, the partial regions of both terminal sequences are reverse complementary, with positions 1 to 27 of 5 -UTR paired with positions 100 to 126 of 3 -UTR, and potentially fold into a panhandle structure with G values of −105 kcal/mol ( Figure 4C ).

To determine whether virions are associated with MtMV1, possible viral proteins were extracted from the mycelia of strain WJB-5 and subjected to ultracentrifugation in stepwise sucrose gradients. The sucrose fraction containing the MtMV1 dsRNAs was further purified and subjected to observation with a transmission electron microscope (TEM). No virus-like particles were observed under TEM, confirming that MtMV1 dsRNAs are naked like other mitoviruses. 

To cure the MtMV1 from the WJB-5 strain and the other five strains of M. theae (WJB-1-1, WJB-1-2, WJB-1-9, WJB-1-12, and WJB-1-20), the hyphal tipping method combined with antiviral extracts (incorporating green-leaf tea powder in the culture medium) was used. After three sequential rounds of subculturing, MtMV1 was detected by dsRNA extraction. The results showed that MtMV1 was cured from all detected strains (Figure 5B) , which was further confirmed by RT-PCR identification resulting in no target bands (1.98 kb in size) (Supplementary Figure 2) .

To investigate whether MtMV1 can be vertically transmitted to the next generation through fungal conidia, conidial suspension prepared from the strain WJB-5 were smeared on PDA, and a total of 100 single-conidium generated colonies were subjected to dsRNA extraction. The results showed that MtMV1 was retained in all these subisolates, suggesting that MtMV1 has efficiently transmitted through the conidia of M. theae ( Figure 5A ).

To access the effect of MtMV1 on the host fungus, six MtMV1infected and six MtMV1-free strains of Melanconiella spp. were cultured on PDA in darkness to access their growth rates and morphologies. Most of these strains had similar morphologies and growth rates, ranging from 4.4 to 6.1 mm/d for the MtMV1infected strains and 4.1 to 5.4 mm/d for the MtMV1-free ones, suggesting that no obvious effects were observed for MtMV1 related to the morphologies and growth rates of their host strains (Supplementary Figure 3) .

A novel mycovirus, named MtMV1, was identified infecting the phytopathogenic fungus M. theae. MtMV1 shares most characteristics of members in the family Mitoviridae considering that mitoviruses have the simplest positive ssRNA genomes that harbor a single ORF encoding a putative viral RdRp for their replication (Cole et al., 2000; Hillman and Cai, 2013; Khalifa and Pearson, 2013) , capsidless (Hillman and Cai, 2013) , and their RdRps contain six conserved amino acid motifs (I-IV) (Osaki et al., 2015; Ran et al., 2016; Chen et al., 2017; Wang et al., 2019; Kocanová et al., 2020; Liu et al., 2021) . Consequently, MtMV1 together with GaMV20 and CFMV1 were placed in cluster I as the third member in this group.

Since mitoviruses essentially replicate within mitochondria of their hosts (Polashock and Hillman, 1994; Cole et al., 2000; Osaki et al., 2005; Nibert, 2017; Fonseca et al., 2021) , MtMV1 genome shares the most features of mitochondrial codons described before, e.g., the MtMV1 genome contains ten UGA codons, which encode Trp rather than function as a stop codon (Nibert, 2017; Vong et al., 2019) , while only two UGG, and has rich of A-U content (62.37%) like other members (59.55 to 73.25%). Moreover, the initiation codon of MtMV1 ORF has a preference to be either A or U with a frequency of 68.3% in the third nucleotide position, similar to those of other members (Liu et al., 2019; Torres-Trenas and Pérez-Artés, 2020; Wang et al., 2021) and has a UAG stop codon, similar to the observation in other members with either UAA or UAG (Nibert, 2017) . Additionally, the 5 -and 3 -terminus sequences of MtMV1 could be folded into stem-loop structures, and their reverse complementary regions allow to form a stable panhandle structure, which is also observed in other members (Hong et al., 1999; Khalifa and Pearson, 2013; Wang et al., 2021) . Whereas, the size of MtMV1 3 -UTR (151 nt) is relatively longer than its 5 -UTR (45 nt), it is different as compared with other mitoviruses since their 3 -UTR sizes are all shorter than their 5 -UTR (ranging from 1 to 413 nucleotide longer) (Stielow et al., 2011; Xie and Ghabrial, 2012; Nibert et al., 2018) . It is worth noting that MtMV1 has the smallest 5 -UTR in size compared with other members (ranging from 49 to 472 nt) in the family Mitoviridae whether some biological functions related to the unique traits of MtMV1 terminal sequences require further studies. Moreover, MtMV1 harbors a genomic sequence obviously divergent from those of most members in the family Mitoviridae as exemplified by having Thr-rich RdRp motifs as compared with other members (Figure 3) .

MtMV1 is efficiently transmitted through all conidia of strain WJB-5 (100 conidium-generated colonies were checked) in agreement with those mitoviruses detected in a high percentage of the individual spore isolates, as exemplified by Fusarium circinatum mitovirus 1 (FcMV1) and Fusarium circinatum mitovirus 2-2 (FcMV2-2) that showed vertical transmission rates between 60 and 100% depending on the fungal isolate (Romeralo et al., 2018) , however, both Fusarium verticillioides mitovirus 1 (FvMV1) and Fusarium andiyazi mitovirus 1 strain 162 (FaMV1-162) showed 100% vertical transmission rate through fungal conidia (Jacquat et al., 2020) . In contrast, some mycoviruses are less efficiently transmitted through conidia, as exemplified by Colletotrichum camelliae filamentous virus 1 (CcFV-1) and Pestalotiopsis theae chrysovirus-1 (PtCV1) (Jia et al., 2017; Zhou et al., 2021) . Generally, the efficient transmission of mitoviruses through asexual spores and sclerotia has hindered to understand mitovirus biology, partly due to lacking mitovirus-free subisolates used for infection, although a protoplast fusion-based protocol for horizontal transmission of a mitovirus has been developed (Shahi et al., 2019) . In light of several reports of tea polyphenols exhibiting antiviral activities against various viruses, especially positive-sense single-stranded RNA viruses (Okada, 1978; Zhang et al., 2020; Mhatre et al., 2021) , we tried the elimination of MtMV1 from the fungal host, it is worth noting that MtMV1 was easily cured by the hyphal tipping method combined with green-leaf tea powder in the culture medium as examined by dsRNA extraction and RT-PCR identification, which is expected to provide an efficient approach to obtain mitovirus-free host strains. Moreover, MtMV1 confers no obvious effects on the growth rates and morphologies as demonstrated by some homologous and isogenous strains. Thus, we did not proceed to viral transmission for more strict assessment of its biological effects.

Collectively, a novel mitovirus (MtMV1) from M. theae was characterized based on its taxonomic analysis, molecular characterization, genome organization, electron microscope observation, and some biological traits. MtMV1 harbors a genomic sequence obviously divergent with those of most members in the family Mitoviridae and some unique characteristics unreported in known members, contributing helpful information to understand the evolution, molecular and biological traits of mycoviruses.

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/ Supplementary Material. 

A Comprehensive Review on Mycoviruses as

A mitovirus isolated from the phytopathogenic fungus Alternaria brassicicola

Detection of an RNA-dependent RNA polymerase in mitochondria from a mitovirus-infected isolate of the dutch elm disease fungus. Ophiostoma novo-ulmi

Host and geographic range extensions of Melanconiella, with a new species M. Cornuta in China

Characterization of a novel mitovirus of the sand fly Lutzomyia longipalpis using genomic and virus-host interaction signatures

Origin, adaptation and evolutionary pathways of fungal viruses

Viruses of plant pathogenic fungi

Complete genome sequence of a Fusarium graminearum double-stranded RNA virus in a newly proposed family

The family Narnaviridae: simplest of RNA viruses

Evolutionary relationships among putative RNA-dependent RNA polymerases encoded by a mitochondrial virus-like RNA in the dutch elm disease fungus, Ophiostoma novo-ulmi, by other viruses and virus-like RNAs and by the Arabidopsis mitochondrial genome

Multiple mitochondrial viruses in an isolate of the dutch elm disease fungus Ophiostoma novo-ulmi

A survey of mycoviral infection in Fusarium spp. isolated from maize and sorghum in argentina identifies the first mycovirus from Fusarium verticillioides

A dsRNA virus with filamentous viral particles

A novel partitivirus in the hypovirulent isolate qt5-19 of the plant pathogenic fungus Botrytis cinerea

Molecular characterization of three mitoviruses co-infecting a hypovirulent isolate of Sclerotinia sclerotiorum fungus

A novel mitovirus detected in Diaporthe rudis, a fungus associated with Phomopsis dieback on grapevines

Complete nucleotide sequence and genome organization of a dsRNA partitivirus infecting Pleurotus ostreatus

Analysis of RNA sequencing data using CLC genomics workbench

Molecular characterization of a novel mitovirus from the plant-pathogenic fungus Botryosphaeria dothidea

A victorivirus and two novel mitoviruses co-infected the plant pathogen Nigrospora oryzae

Fungal negative-stranded RNA virus that is related to bornaviruses and nyaviruses

Molecular characterization of a novel mycovirus in the cultivated mushroom, Lentinula edodes

Antiviral activity of green tea and black tea polyphenols in prophylaxis and treatment of COVID-19: a review

Novel mitoviruses and a unique tymo-like virus in hypovirulent and virulent strains of the Fusarium head blight fungus, Fusarium boothii

Biological and molecular characterization of Chenopodium quinoa mitovirus 1 reveals a distinct small RNA response compared to those of cytoplasmic RNA viruses

Mitovirus UGA(Trp) codon usage parallels that of host mitochondria

Mitovirus and mitochondrial coding sequences from basal fungus Entomophthora muscae

Evidence for contemporary plant mitoviruses

Hypovirulence: mycoviruses at the fungal-plant interface

Antiviral effects of tea catechins and black tea theaflavins on plant viruses

Nucleotide sequence of a mitochondrial RNA virus from the plant pathogenic fungus, Helicobasidium mompa Tanaka

Isolation and characterization of two mitoviruses and a putative alphapartitivirus from Fusarium spp

Mycoviruses of filamentous fungi and their relevance to plant pathology

A small mitochondrial double-stranded (ds) RNA element associated with a hypovirulent strain of the chestnut blight fungus and ancestrally related to yeast cytoplasmic T and W dsRNAs

Co-infection of a hypovirulent isolate of Sclerotinia sclerotiorum with a new botybirnavirus and a strain of a mitovirus

Isolation and characterization of a novel mycovirus, PeSV, in Pleurotus eryngii and the development of a diagnostic system for it

Vertical transmission of Fusarium circinatum mitoviruses FcMV1 and FcMV2-2 via microconidia

Complete genome sequence of a novel mitovirus from the ectomycorrhizal fungus Geopora sumneriana

Investigation of host range of and host defense against a mitochondrially replicating mitovirus

A novel Tuber aestivum (Vittad.) mitovirus. Arch. Virol

Molecular characterization of the genome of a partitivirus from the basidiomycete Rhizoctonia solani

Optimization of phenol-chloroform RNA extraction

Characterization and incidence of the first member of the genus Mitovirus identified in the phytopathogenic species Fusarium oxysporum

Multigene phylogeny and taxonomy of the genus Melanconiella (Diaporthales)

Beta vulgaris mitovirus 1 in diverse cultivars of beet and chard

The complete genome sequence of a new mitovirus from the phytopathogenic fungus Colletotrichum higginsianum

Characterization of a novel RNA virus from the phytopathogenic fungus Leptosphaeria biglobosa related to members of the genus Mitovirus

Hypovirulence and double-stranded RNA in Botrytis cinerea

Molecular characterizations of two mitoviruses co-infecting a hyovirulent isolate of the plant pathogenic fungus Sclerotinia sclerotiorum

A novel mycovirus closely related to hypoviruses that infects the plant pathogenic fungus Sclerotinia sclerotiorum

A mitovirus related to plant mitochondrial gene confers hypovirulence on the phytopathogenic fungus Sclerotinia sclerotiorum

Sclerotinia minor endornavirus 1, a novel pathogenicity debilitation-associated mycovirus with a wide spectrum of horizontal transmissibility

A rapid silica spin columnbased method of RNA extraction from fruit trees for RT-PCR detection of viruses

Antimicrobial activities of tea polyphenol on phytopathogens: a review

A geminivirusrelated DNA mycovirus that confers hypovirulence to a plant pathogenic fungus

Epigallocatechin-3-gallate enhances tomato resistance to tobacco mosaic virus by modulating RBOH1-dependent H2O2 signaling

A mycovirus modulates the endophytic and pathogenic traits of a plant associated fungus

The authors thank the financial support by the National Natural Science Foundation of China and Hubei Hongshan Laboratory.

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.

Conflict of Interest: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.