Characterization of a Novel Mitovirus Infecting Melanconiella theae Isolated From Tea Plants

Abstract

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.

Introduction

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¹, 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).

Materials and Methods

Fungus Isolation and Identification

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.

Double-Stranded RNA (dsRNA) Extraction

Viral dsRNA was extracted from the fungal mycelial mass using the silica spin column-based method as described previously (Yang et al., 2017). 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).

Complementary DNA (cDNA) Cloning and Sequencing, and RT-PCR Amplication

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, 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². 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).

Phylogenetic and Sequence Analyses

The viral ORF was identified using Expasy³ website, and the amino acid (aa) sequences encoded by the ORF were deduced. 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⁴ provided by NCBI, and multiple amino acid alignments.

Association of Strain WJB-5 With Viral-Like Particles

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 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 phenol-chloroform-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 carbon-formvar 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.

Attempts to Cure Melanconiella Strains

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.

Analysis of Mycovirus Transmission Through Conidia

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⁷ conidia/mL. The resulting conidial suspension was poured on PDA plates and incubated at 20°C in darkness, and used for culturing of single-conidium 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.

Growth Rate and Morphological Assessments

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 virus-eliminated 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).

Results

Extraction of dsRNA Segments From Melanconiella theae Strains

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 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).

FIGURE 1

Sequence and Molecular Characterization of WJB-5 dsRNA

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 tentatively named Melanconiella theae mitovirus 1 (MtMV1). The complete nucleotide sequence of MtMV1 has been deposited in GenBank with the accession number MW802251.

Genomic Organization and Phylogenetic Analysis

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).

FIGURE 2

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).

FIGURE 3

Predicted Secondary Structures of the Terminal Sequence

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).

FIGURE 4

Attempts to Check the Viral-Like Particles From WJB-5 Strain

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.

Attempts to Cure WJB-5 Strain From Melanconiella Theae Mitovirus 1

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).

FIGURE 5

Vertical Transmission of Melanconiella Theae Mitovirus 1 Through Conidia

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).

Melanconiella theae Mitovirus 1 Confers No Obvious Effects on the Growth Rates and Morphologies of the Host Fungus

To access the effect of MtMV1 on the host fungus, six MtMV1-infected 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 MtMV1-infected 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).

Discussion

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 sub-isolates 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.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

• 1

  AbidM.KhanM. A. U.MushtaqS.AfzaalS.HaiderM. S. (2018). A Comprehensive Review on Mycoviruses as Biological Control Agent.World J. Biol. Biotechnol.3187–192. 10.33865/wjb.003.02.0146

  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  ChenY.ShangH. H.YangH. Q.da GaoB.ZhongJ. (2017). A mitovirus isolated from the phytopathogenic fungus Alternaria brassicicola.Arch. Virol.1622869–2874. 10.1007/s00705-017-3371-8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 3

  ColeT. E.HongY.BrasierC. M.BuckK. W. (2000). Detection of an RNA-dependent RNA polymerase in mitochondria from a mitovirus-infected isolate of the dutch elm disease fungus. Ophiostoma novo- ulmi.Virology268239–243. 10.1006/viro.1999.0097

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 4

  DuZ.FanX. L.YangQ.TianC. M. (2017). Host and geographic range extensions of Melanconiella, with a new species M. Cornuta in China.Phytotaxa327252–260. 10.11646/phytotaxa.327.3.4

  • CrossRef
  • Google Scholar
  • View reference in article

• 5

  FonsecaP.FerreiraF.da SilvaF.OliveiraL. S.MarquesJ. T.Goes-NetoA.et al (2021). Characterization of a novel mitovirus of the sand fly Lutzomyia longipalpis using genomic and virus–host interaction signatures.Viruses13:9. 10.3390/v13010009

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  GhabrialS. A. (1998). Origin, adaptation and evolutionary pathways of fungal viruses.Virus Genes16119–131. 10.1023/a:1007966229595

  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  GhabrialS. A.SuzukiN. (2009). Viruses of plant pathogenic fungi.Annu. Rev. Phytopathol.47353–384. 10.1146/annurev-phyto-080508-081932

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 8

  HeH.ChenX.LiP.QiuD.GuoL. (2018). Complete genome sequence of a Fusarium graminearum double-stranded RNA virus in a newly proposed family, Alternaviridae.Genome Announc.6e00064–18. 10.1128/genomeA.00064-18

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 9

  HillmanB. I.CaiG. (2013). The family Narnaviridae: simplest of RNA viruses.Adv. Virus Res.86149–176. 10.1016/B978-0-12-394315-6.00006-4

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  HongY.ColeT. E.BrasierC. M.BuckK. W. (1998). 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.Virology246158–169. 10.1006/viro.1998.9178

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 11

  HongY.DoverS. L.ColeT. E.BrasierC. M.BuckK. W. (1999). Multiple mitochondrial viruses in an isolate of the dutch elm disease fungus Ophiostoma novo-ulmi.Virology258118–127. 10.1006/viro.1999.9691

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 12

  JacquatA. G.TheumerM. G.CañizaresM. C.DebatH. J.IglesiasJ.PedrajasM. D. G.et al (2020). A survey of mycoviral infection in Fusarium spp. isolated from maize and sorghum in argentina identifies the first mycovirus from Fusarium verticillioides.Viruses12:1161. 10.3390/v12101161

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  JiaH.DongK.ZhouL.WangG.HongN.JiangD.et al (2017). A dsRNA virus with filamentous viral particles.Nat. Commun.8:168. 10.1038/s41467-017-00237-9

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  KamaruzzamanM.HeG.WuM.ZhangJ.YangL.ChenW.et al (2019). A novel partitivirus in the hypovirulent isolate qt5-19 of the plant pathogenic fungus Botrytis cinerea.Viruses11:24. 10.3390/v11010024

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  KhalifaM. E.PearsonM. N. (2013). Molecular characterization of three mitoviruses co-infecting a hypovirulent isolate of Sclerotinia sclerotiorum fungus.Virology44122–30. 10.1016/j.virol.2013.03.002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  KocanováM.EichmeierA.BotellaL. (2020). A novel mitovirus detected in Diaporthe rudis, a fungus associated with Phomopsis dieback on grapevines.Arch. Virol.1652405–2408. 10.1007/s00705-020-04755-8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  LimW. S.JeongJ. H.JeongR. D.YooY. B.YieS. W.KimK. H. (2005). Complete nucleotide sequence and genome organization of a dsRNA partitivirus infecting Pleurotus ostreatus.Virus Res.108111–119. 10.1016/j.virusres.2004.08.017

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 18

  LiuC. H.DiY. P. (2020). “Analysis of RNA sequencing data using CLC genomics workbench,” in Methods in Molecular Biology, edsKeohavongP.SinghK. P.GaoW. (New York: Springer), 61–113. 10.1007/978-1-0716-0223-2_4

  • CrossRef
  • Google Scholar
  • View reference in article

• 19

  LiuH.LiuM.ZhuH.ZhongJ.LiaoX.ZhouQ. (2021). Molecular characterization of a novel mitovirus from the plant-pathogenic fungus Botryosphaeria dothidea.Arch. Virol.166633–637. 10.1007/s00705-020-04886-y

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  LiuH.LiuR.LiC. X.WangH.ZhuH. J.da GaoB.et al (2019). A victorivirus and two novel mitoviruses co-infected the plant pathogen Nigrospora oryzae.Viruses11:83. 10.3390/v11010083

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 21

  LiuL.XieJ.ChengJ.FuY.LiG.YiX.et al (2014). Fungal negative-stranded RNA virus that is related to bornaviruses and nyaviruses.Proc. Natl. Acad. Sci. U. S. A.11112205–12210. 10.1073/pnas.1401786111

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  MagaeY. (2012). Molecular characterization of a novel mycovirus in the cultivated mushroom, Lentinula edodes.Virol. J.9:60. 10.1186/1743-422X-9-60

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  MhatreS. S. T.NaikS.PatravaleV. (2021). Antiviral activity of green tea and black tea polyphenols in prophylaxis and treatment of COVID-19: a review.Phytomedicine85:153286. 10.1016/j.phymed.2020.153286

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 24

  MizutaniY.AbrahamA.UesakaK.KondoH.SugaH.SuzukiN.et al (2018). Novel mitoviruses and a unique tymo-like virus in hypovirulent and virulent strains of the Fusarium head blight fungus, Fusarium boothii.Viruses10:584. 10.3390/v10110584

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  NervaL.ViganiG.di SilvestreD.CiuffoM.ForgiaM.ChitarraW.et al (2019). Biological and molecular characterization of Chenopodium quinoa mitovirus 1 reveals a distinct small RNA response compared to those of cytoplasmic RNA viruses.J. Virol.93e01998–18. 10.1128/jvi.01998-18

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 26

  NibertM. L. (2017). Mitovirus UGA(Trp) codon usage parallels that of host mitochondria.Virology50796–100. 10.1016/j.virol.2017.04.010

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 27

  NibertM. L.DebatH. J.MannyA. R.GrigorievI. V.de Fine LichtH. H. (2019). Mitovirus and mitochondrial coding sequences from basal fungus Entomophthora muscae.Viruses11:351. 10.3390/v11040351

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 28

  NibertM. L.VongM.FugateK. K.DebatH. J. (2018). Evidence for contemporary plant mitoviruses.Virology51814–24. 10.1016/j.virol.2018.02.005

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 29

  NussD. L. (2005). Hypovirulence: mycoviruses at the fungal-plant interface.Nat. Rev. Microbiol.3632–642. 10.1038/nrmicro1206

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 30

  OkadaF. (1978). Antiviral effects of tea catechins and black tea theaflavins on plant viruses.JARQ1227–32.

  • Google Scholar
  • View reference in article

• 31

  OsakiH.NakamuraH.NomuraK.MatsumotoN.YoshidaK. (2005). Nucleotide sequence of a mitochondrial RNA virus from the plant pathogenic fungus, Helicobasidium mompa Tanaka.Virus Res.10739–46. 10.1016/j.virusres.2004.06.010

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 32

  OsakiH.SasakiA.NomiyamaK.SekiguchiH.TomiokaK.TakeharaT. (2015). Isolation and characterization of two mitoviruses and a putative alphapartitivirus from Fusarium spp.Virus Genes50466–473. 10.1007/s11262-015-1182-0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 33

  PearsonM. N.BeeverR. E.BoineB.ArthurK. (2009). Mycoviruses of filamentous fungi and their relevance to plant pathology.Mol. Plant Pathol.10115–128. 10.1111/j.1364-3703.2008.00503.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 34

  PolashockJ. J.HillmanB. I. (1994). 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.Proc. Natl. Acad. Sci. U. S. A.918680–8684. 10.1073/pnas.91.18.8680

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 35

  RanH.LiuL.LiB.JiasenC.FuY.DaohongJ.et al (2016). Co-infection of a hypovirulent isolate of Sclerotinia sclerotiorum with a new botybirnavirus and a strain of a mitovirus.Virol. J.13:92. 10.1186/s12985-016-0550-2

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 36

  RoH. S.KangE. J.YuJ. S.LeeT. S.LeeC. W.LeeH. S. (2007). Isolation and characterization of a novel mycovirus, PeSV, in Pleurotus eryngii and the development of a diagnostic system for it.Biotechnol. Lett.29129–135. 10.1007/s10529-006-9206-4

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 37

  RomeraloC.BezosD.Martínez-álvarezP.DiezJ. J. (2018). Vertical transmission of Fusarium circinatum mitoviruses FcMV1 and FcMV2-2 via microconidia.Forests9:356. 10.3390/f9060356

  • CrossRef
  • Google Scholar

• 38

  SahinE.AkataI. (2019). Complete genome sequence of a novel mitovirus from the ectomycorrhizal fungus Geopora sumneriana.Arch. Virol.1642853–2857. 10.1007/s00705-019-04367-x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 39

  ShahiS.Eusebio-CopeA.KondoHidekiB.HillmanI.SuzukiN. (2019). Investigation of host range of and host defense against a mitochondrially replicating mitovirus.J. Virol.93e1503–e1518.

  • Google Scholar

• 40

  StielowB.KlenkH. P.WinterS.MenzelW. (2011). A novel Tuber aestivum (Vittad.) mitovirus.Arch. Virol.1561107–1110. 10.1007/s00705-011-0998-8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 41

  StraussE. E.LakshmanD. K.TavantzisS. M. (2000). Molecular characterization of the genome of a partitivirus from the basidiomycete Rhizoctonia solani.J. Gen. Virol.81549–555.

  • Google Scholar

• 42

  ToniL. S.GarciaA. M.JeffreyD. A.JiangX.StaufferB. L.MiyamotoS. D.et al (2018). Optimization of phenol-chloroform RNA extraction.MethodsX5599–608. 10.1016/j.mex.2018.05.011

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 43

  Torres-TrenasA.Pérez-ArtésE. (2020). Characterization and incidence of the first member of the genus Mitovirus identified in the phytopathogenic species Fusarium oxysporum.Viruses12:279. 10.3390/v12030279

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 44

  VoglmayrH.RossmanA. Y.CastleburyL. A.JaklitschW. M. (2012). Multigene phylogeny and taxonomy of the genus Melanconiella (Diaporthales).Fungal Divers.571–44. 10.1007/s13225-012-0175-8

  • CrossRef
  • Google Scholar

• 45

  VongM.MannyA. R.SmithK. L.GaoW.NibertM. L. (2019). Beta vulgaris mitovirus 1 in diverse cultivars of beet and chard.Virus Res.26580–87. 10.1016/j.virusres.2019.02.008

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 46

  WangH.LiuH.ZhouQ. (2021). The complete genome sequence of a new mitovirus from the phytopathogenic fungus Colletotrichum higginsianum.Arch. Virol.1661481–1484. 10.1007/s00705-021-04996-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 47

  WangY.XuZ.ChengJ.FuY.JiangD.XieJ. (2019). Characterization of a novel RNA virus from the phytopathogenic fungus Leptosphaeria biglobosa related to members of the genus Mitovirus.Arch. Virol.164913–916. 10.1007/s00705-018-04143-3

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 48

  WuM. D.ZhangL.LiG. Q.JiangD. H.HouM. S.HuangH. C. (2007). Hypovirulence and double-stranded RNA in Botrytis cinerea.Phytopathology971590–1599. 10.1094/PHYTO-97-12-1590

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 49

  XieJ.GhabrialS. A. (2012). Molecular characterizations of two mitoviruses co-infecting a hyovirulent isolate of the plant pathogenic fungus Sclerotinia sclerotiorum.Virology42877–85. 10.1016/j.virol.2012.03.015

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 50

  XieJ.XiaoX.FuY.LiuH.ChengJ.GhabrialS. A.et al (2011). A novel mycovirus closely related to hypoviruses that infects the plant pathogenic fungus Sclerotinia sclerotiorum.Virology41849–56. 10.1016/j.virol.2011.07.008

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 51

  XuZ.WuS.LiuL.ChengJ.FuY.JiangD.et al (2015). A mitovirus related to plant mitochondrial gene confers hypovirulence on the phytopathogenic fungus Sclerotinia sclerotiorum.Virus Res.197127–136. 10.1016/j.virusres.2014.12.023

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 52

  YangD.WuM.ZhangJ.ChenW.LiG.YangL. (2018). Sclerotinia minor endornavirus 1, a novel pathogenicity debilitation-associated mycovirus with a wide spectrum of horizontal transmissibility.Viruses10:589. 10.3390/v10110589

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 53

  YangF.WangG.XuW.HongN. (2017). A rapid silica spin column-based method of RNA extraction from fruit trees for RT-PCR detection of viruses.J. Virol. Methods24761–67. 10.1016/j.jviromet.2017.05.020

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 54

  YangY.ZhangT. (2019). Antimicrobial activities of tea polyphenol on phytopathogens: a review.Molecules24:816. 10.3390/molecules24040816

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 55

  YuX.LiB.FuY.JiangD.GhabrialS. A.LiG.et al (2010). A geminivirus-related DNA mycovirus that confers hypovirulence to a plant pathogenic fungus.Proc. Natl. Acad. Sci. U. S. A.1078387–8392. 10.1073/pnas.0913535107

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 56

  ZhangX. N.LiaoY. W.WangX. R.ZhangL.AhammedG. J.LiQ. Y.et al (2020). Epigallocatechin-3-gallate enhances tomato resistance to tobacco mosaic virus by modulating RBOH1-dependent H2O2 signaling.Plant Physiol. Biochem.150263–269. 10.1016/j.plaphy.2020.03.008

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 57

  ZhouL.LiX.Kotta-LoizouI.DongK.LiS.NiD.et al (2021). A mycovirus modulates the endophytic and pathogenic traits of a plant associated fungus.ISME J.151893–1906. 10.1038/s41396-021-00892-3

  • Pubmed Abstract
  • CrossRef
  • Google Scholar