Skip to main content

ORIGINAL RESEARCH article

Front. Microbiol., 25 November 2024
Sec. Microbial Symbioses

Unveiling mycoviral diversity in Ophiocordyceps sinensis through transcriptome analyses

Qin Kang,Qin Kang1,2Jihong Zhang,Jihong Zhang1,3Fangzhou ChenFangzhou Chen4Caihong DongCaihong Dong5Qilian Qin,Qilian Qin1,3Xuan Li,Xuan Li1,3Hongtuo Wang,Hongtuo Wang1,3Huan Zhang,
Huan Zhang1,3*Qian Meng,
Qian Meng1,3*
  • 1State Key Laboratory of Integrated Management of Pest Insects and Rodents, Institute of Zoology, Chinese Academy of Sciences, Beijing, China
  • 2University of Chinese Academy of Sciences, Beijing, China
  • 3Beijing Yun’an Bio-tech Co. Ltd., Beijing, China
  • 4China Pharmaceutical University, School of Pharmacy, Nanjing, China
  • 5State Key Laboratory of Mycology, Institute of Microbiology, Chinese Academy of Sciences, Beijing, China

Ophiocordyceps sinensis, an entomopathogenic fungus, infects larvae from the Lepidoptera: Hepialidae family, forming the valuable Chinese cordyceps. Mycoviruses are widespread across major lineages of filamentous fungi, oomycetes, and yeasts and have the potential to influence fungal biology and ecology. This study aimed to detect mycovirus within O. sinensis by isolating double-stranded RNA from six stains for transcriptomic sequencing and analyzing publicly available transcriptome data from 13 O. sinensis representative samples. Our analysis revealed 13 mycoviruses, with nine reported for the first time in O. sinensis. These mycoviruses are distributed across five families—Partitiviridae, Mitoviridae, Narnaviridae, Botourmiaviridae, Deltaflexiviridae—and two unclassified lineages, Ormycovirus and Vivivirus. This study also revealed frequent coinfections within individual O. sinensis strains and dynamic shifts in viral composition during fungal development. These findings enhance our knowledge of mycovirus diversity within O. sinensis and provide new insights into their taxonomy.

1 Introduction

Chinese cordyceps is an insect-fungus complex generated by Ophiocordyceps sinensis (Ophiocordyceps, Ophiocordycipitaceae, Hypocreales, Sordariomycetes, Ascomycota) infecting insect larvae (Lepidoptera, Hepialidae). Beyond its natural soil habitat, O. sinensis can endogenously colonize plants (Baral, 2017). Chinese cordyceps is widely distributed across the alpine grasslands of the Himalayas and the mountainous regions of the Tibetan Plateau (Zhang et al., 2012). As a traditional Chinese medicine for centuries, Chinese cordyceps has become one of the most expensive natural products, sometimes surpassing gold in price, and is key to the local economy (Fan and He, 2024). Artificial cultivation of Chinese cordyceps is an effective way to alleviate supply and demand pressures and ecological impacts. However, technological bottlenecks still exist, such as low natural infection rates, long production cycles, and disease susceptibility caused by other pathogens (Qin et al., 2018).

Mycoviruses infect and replicate within fungi. They have been increasingly identified since the first discovery in Agaricus bisporus in 1962 (Hollings, 1962). Advances in high throughput sequencing (HTS) have facilitated the identification of mycoviruses in various fungal species, leading to a substantial increase in the number of described viruses infecting fungi in recent years (Villan Larios et al., 2023). Most mycoviruses have genomes consisting of positive-sense (+) or negative-sense (−) single-stranded RNA (ssRNA), double-stranded RNA (dsRNA) or even DNA (Sato and Suzuki, 2023). Currently, mycoviruses can be identified via dsRNA analysis and HTS. With the accessibility of fungal transcriptome data on platforms like NCBI Sequence Read Archive (SRA), the presence of mycoviruses in less-studied fungal populations can be explored (Hough et al., 2023). Although most mycovirus infections are asymptomatic, certain mycoviruses can influence various aspects of fungal physiology, including growth rate, reproduction modes, adaptability, and pathogenicity (Kotta-Loizou, 2021; Sato and Suzuki, 2023). In entomopathogenic fungi, particularly in the Metarhizium and Beauveria bassiana, mycoviruses have been relatively well-documented (Filippou et al., 2018; Hwang et al., 2023; Guo et al., 2024; Shi et al., 2024). In recent years, three different viruses—Cordyceps chanhua alternavirus 1 (CcAV1), Cordyceps chanhua victorivirus 1 (CcV1), and Cordyceps chanhua partitivirus 1 (CchPV1)—have been identified within Cordyceps chanhua. CchPV1 has been shown to affect development and multi-stress tolerance in the host fungus (Zhang et al., 2022; Zhu et al., 2022a; Zhu et al., 2022b). However, research on mycoviruses in O. sinensis remains limited.

While residing in extreme ecosystems, O. sinensis interacts with insects, plants, and other microorganisms, where mycoviruses likely influence its development and prevalence. Early studies mining public databases suggest the existence of mycoviruses in O. sinensis, but systematic investigation is lacking (Forgia et al., 2022; Gilbert et al., 2019). To address this gap, we selected six rapidly growing O. sinensis strains from our laboratory for dsRNA screening and transcriptomic sequencing. Additionally, we analyzed publicly available O. sinensis transcriptome data in NCBI SRA to assess the diversity of mycoviruses. Our findings provide valuable insights into the diversity, origins and evolutionary relationships of mycoviruses in O. sinensis, as well as their genomic structure and potential functions.

2 Materials and methods

2.1 Isolates and growth conditions

Strains IOZ14 and IOZ23 were isolated from sclerotium samples of O. sinensis, collected in Guoluo Tibetan Autonomous Prefecture, Qinghai Province, China. Strains IOZ25, IOZ26, and IOZ31 were derived from fresh Chinese cordyceps harvested in Xiaojin County, Sichuan Province, China, while strain IOZ29 was derived from fresh Chinese cordyceps collected in Yushu Tibetan Autonomous Prefecture, Qinghai Province, China. O. sinensis was cultured in peptone potato dextrose agar (PPDA) medium at 18°C to harvest conidia and mycelium, following previously described methods (Wu et al., 2022a).

2.2 Extraction and screening of dsRNA

After 1 month of cultivation, fresh mycelium samples were collected for nucleic acid extraction. The dsRNA was extracted using the CF11 cellulose powder adsorption method (Morris and Dodds, 1979), and then digested by S1 nuclease and DNase I (TaKaRa, Dalian, China) to remove DNA and ssRNA. The purified dsRNA was analyzed on a 1% agarose gel to detect mycovirus presence.

2.3 Extraction of total RNA and sequencing

Total RNA was extracted using TRIzol reagent (Thermo Fisher Scientific, Waltham, USA) following the manufacturer’s instructions, and then treated with DNase I (TaKaRa, Dalian, China) to eliminate genomic DNA. Two μg of RNA from each IOZ strain was pooled and sent to Novogene Company (Beijing, China) for library construction (using ribosomal RNA deletion method) and sequencing on the Illumina MiSeq 2000/2500.

2.4 Assembly and annotation of mycovirus sequences

The raw sequencing data were processed using fastp (version 0.23.2) with default parameters to filter out low-quality reads and adapter sequences (Chen et al., 2018). The resulting clean data were aligned to the O. sinensis genome (GCA_012934285.1) using HISAT2 (Kim et al., 2019). Subsequently, the unmapped reads were extracted and assembled de novo using Megahit (v1.2.9) with default parameters (Li D. et al., 2016). After removing redundancy, the assembled sequences were annotated by the Non-Redundant Protein Sequence Database (nr) database of NCBI to identify partial mycovirus sequences. To calculate the number of reads for the individually identified viruses in each sample, the clean reads were mapped to the final fragments of each virus using samtools (v1.20) (Danecek et al., 2021). The coverage of each mycovirus identified in the public transcriptome was calculated based on the length of the longest contig assembled for each mycovirus in the sample as a percentage of the length of the final viral genome length.

2.5 Reverse-transcription PCR (RT-PCR) and mycovirus full-length determination

To further confirm the presence of mycoviruses in the IOZ strains, reverse transcriptase M-MLV (TaKaRa, Dalian, China) and random primers were used to synthesize first-strand cDNA from 500 ng of RNA. Virus detection was carried out through PCR amplification using TransStart FastPfu DNA Polymerase (TransGen, Beijing, China). Primers specific for individual mycoviruses were designed based on the assembled sequences (Supplementary Table S1). The rapid amplification of cDNA ends (RACE) technique was employed to determine the complete sequences of all mycoviruses infecting the IOZ strain (Lambden et al., 1992). The primers used in the RACE analysis are listed in Supplementary Table S2.

2.6 Mining publicly available Ophiocordyceps sinensis transcriptomic datasets

The search term “(Ophiocordyceps sinensis) AND bioproject_sra[filter] NOT bioproject_gap[filter]” and “(Chinese [All Fields] AND (“Cordyceps”[Organism] OR cordyceps [All Fields])) AND “bioproject sra”[Filter]” were used to retrieve all available O. sinensis raw read files from the NCBI SRA database accessed on March 2, 20241. Initially, a total of 64 bioprojects were returned. Through manual curation, samples were refined based on specific criteria, including the combination of samples from the same strain, exclusion of mixed RNA samples, and removal of samples with short or unavailable sequences. Data for further analysis were randomly drawn from one biological replicate of each control sample or wild-type strain. Additionally, different developmental stages of the same strain of Chinese cordyceps were also analyzed for virus screening.

2.7 Mycoviral genome and phylogenetic analysis

Open reading frames (ORFs) within the viral sequences were determined using NCBI ORF Finder program2. Conserved structural domains were predicted using NCBI Conserved Domain Search (CD-Search) tool3. RNA dependent RNA polymerase (RdRp) protein sequences were retrieved from NCBI database, followed by multiple sequence comparisons using MUSCLE (Edgar, 2022). Subsequently, the sequences were trimmed by trimAl (v1.4.rev15) with the parameters of gt 0.6 and cons 60 (Capella-Gutiérrez et al., 2009). Maximum likelihood phylogenetic analysis was performed using IQ-TREE 2.1.2, with 1,000 bootstrap replicates (Minh et al., 2020). The phylogenetic tree was visualized using MEGA 11 (Tamura et al., 2021).

3 Results

3.1 DsRNA screening

Six IOZ strains were screened for the presence of dsRNA molecules. At least one dsRNA fragment with a size between 2.0–3.0 kilobase (kb) was detected in IOZ14 and IOZ29 treated with DNase I and S1 nuclease (Figure 1A), indicating the existence of mycoviruses in these two strains. The fragment may be the genome of dsRNA viruses or the replicative form of single-stranded (ss) RNA viruses. Considering that the titers of certain RNA viruses are too low to detect (Kondo et al., 2022), a pooled sample of the six strains was subjected to RNA sequencing.

Figure 1
www.frontiersin.org

Figure 1. Detection of viruses in IOZ strains. (A) Detection of total viral dsRNAs in the IOZ strains. The dsRNA samples with (+) or without (−) DNase I and S1 nuclease treatment were analyzed on a 1% agarose gel. (B) RT-PCR amplification detection of viruses in IOZ strains. O. sinensis actin (Osactin) was used as the reference gene. OsOMV1RNA, RNA of Ophiocordyceps sinensis ormycovirus 1; OsOVA, Ophiocordyceps ourmiavirus A; OsMV, Ophiocordyceps sinensis mitovirus; OsNV, Ophiocordyceps sinensis narnavirus.

3.2 Viral sequences identified in the transcriptome of Ophiocordyceps sinensis

Transcriptomic sequencing of the IOZ stains yielded 36,302,134 raw reads. After quality filtering, a total of 32,038,562 clean reads remained, which were assembled into 70 contigs. The contigs that showed greater than 90% similarity with the reference virus at the amino acid level were considered to be the homology. The RT-PCR analyses showed that strains IOZ14 and IOZ29 contained eight viruses, while the other four strains each had seven viruses (Figure 1B). All of the detected viruses were (+) ssRNA viruses. Three of them were reported before in O. sinensis and the others were newly identified. They belong to four distinct taxa: Botourmiaviridae, Mitoviridae, Narnaviridae and Ormycovirus.

3.3 Viral sequences identified from publicly available Ophiocordyceps sinensis transcriptomes

Querying the NCBI SRA database for O. sinensis transcriptomic datasets initially returned 64 BioProjects. After manual screening, there are 13 SRA samples left for further analysis (Table 1). A total of 101 contigs were assembled from the 13 SRA samples and 13 viruses were identified (Supplementary Table S3). The number of reads and coverage validated the presence of each identified virus (Supplementary Table S4). These viruses included the eight identified in our transcriptome analysis, plus five others: Ophiocordyceps sinensis mitovirus 1 (OsMV1) and four novel viruses. Ophiocordyceps sinensis ormycovirus 1 (OsOMV1) was relatively the most prevalent, found in 8 SRA samples. Following closely were Ophiocordyceps ourmiavirus A (OsOVA), Ophiocordyceps sinensis mitovirus 3 (OsMV3) and Ophiocordyceps sinensis narnavirus 2 (OsNV2), each of which appeared in 5 samples. Ophiocordyceps sinensis narnavirus 3 (OsNV3), Ophiocordyceps sinensis partitivirus 1 (OsPV1) and Ophiocordyceps sinensis vivivirus 1 (OsVV1) were each identified in only a single, distinct sample (Table 2).

Table 1
www.frontiersin.org

Table 1. Summary of O. sinensis datasets retrieved from the NCBI SRA database.

Table 2
www.frontiersin.org

Table 2. Summary of mycoviral genomes identified from O. sinensis transcriptomic datasets.

3.4 One virus characterized in the Ormycovirus

Ormycovirus is a newly defined group of viruses that contain two RNA fragments. The RNA1 encodes an RdRP, while the RNA2 encodes a putative protein (Pagnoni et al., 2023). By utilizing tBLASTn analysis against the NCBI available fungal Transcriptome Shotgun Assembly (TSA) databases, ormycovirus was identified in O. sinensis and classified within the group betormycovirus, named Ophiocordyceps sinensis ormycovirus 1 (OsOMV1) (Forgia et al., 2022). The OsOMV1 sequences, previously incomplete in the public databases, were fully characterized in this study using RACE. RNA1 and RNA2 of OsOMV1 were determined to be 2,506nt and 2,156nt long, respectively (Figure 2A). The RNA1 was predicted to encode an RdRp of 820 aa, sharing 43.7% sequence identity with the RdRp encoded by Erysiphe lesion-associated ormycovirus 2 (query coverage: 88%; E-value: 0; and accession: USW07207.1). The RNA2 appeared to encode a 663-aa protein, displaying 36.0% identity with that of Verticillium dahliae ormycovirus (query coverage: 71%; E-value: 3e−93; and accession: WPV08071.1). Both of the 5′ and 3′ termini in the two RNA fragments are highly conserved (Figure 2B). In addition, the terminal regions exhibited a high similarity to those of Starmerella bacillaris ormycovirus 1 (SbOMV1) (Forgia et al., 2022).

Figure 2
www.frontiersin.org

Figure 2. Genome organization and terminal conservation of OsOMV1. (A) Genomic structure of OsOMV1. ORFs are shown as boxes. (B) Terminal conservation between OsOMV1 RNA1 and RNA2 is shown by sequence alignment.

3.5 One virus characterized in the family Botourmiaviridae

Viruses of the family Botourmiaviridae are capable of infecting plants and filamentous fungi. Botourmiaviridae encompass 12 genera with (+) ssRNA genomes that can be either single or multiple segments (Ayllón et al., 2020). A virus from this family has been identified in O. sinensis. The viral partial sequence has been deposited in the NCBI database and designated as Ophiocordyceps ourmiavirus A (OsOVA). Its complete genome sequence was obtained through RACE. Notably, the OsOVA possesses a single segment with a full length of 2,967 nt, encoding a 682-aa RdRp (Figure 3A). The nucleotide sequence showed 51.0% identity with the RdRp sequence from Hulunbuir botou tick virus 5 (query coverage: 67%; E-value: 0; and accession: UYL95443.1). Phylogenetic analysis indicated that OsOVA exhibits a close relationship with Plasmopara viticola lesion associated ourmia-like virus 49, and 50, both belonging to the genus Gammascleroulivirus (Figure 3B).

Figure 3
www.frontiersin.org

Figure 3. Genome organization and phylogenetic analysis of OsOVA. (A) Genomic structure of OsOVA. ORF is shown as a box. (B) Evolutionary relationships of the viruses in the Botourmiaviridae family. RdRp amino acid sequences from different genera of Botourmiaviridae were aligned to derive a phylogenetic tree using the maximum likelihood method with the LG + I + G substitution model. Narnaviridae is set as the outgroup. Bootstrap values (%) of 1,000 replicates that are greater than 50% are labeled on the branches. The virus OsOVA is marked with a red triangle.

3.6 Three viruses characterized in the family Mitoviridae

Mitoviridae, a family of capsid-less (+) ssRNA viruses, associates with the mitochondria of fungi, plants or invertebrates. Viruses in the family are predicted to harbor a single continuous ORF adhering to the mitochondrial translation system, where UGA encodes tryptophan (Begeman et al., 2023). Gilbert et al. previously identified Ophiocordyceps sinensis mitovirus 1 and 2 (OsMV1, OsMV2) through analysis of data available in NCBI SRA database (Gilbert et al., 2019). We detected two mitoviral genomes in the IOZ strains. One corresponded to OsMV2, while the other represented a novel mitovirus, designated as Ophiocordyceps sinensis mitovirus 3 (OsMV3). The complete genome of OsMV3 spans 2,600 nt and encodes a 738-aa RdRp protein (Figure 4A). The BLASTx analysis revealed 50.5% identity between the RdRp of OsMV3 and that of OsMV1 (query coverage: 83%; E-value: 0; and accession: AZT88623.1).

Figure 4
www.frontiersin.org

Figure 4. Genome organization and phylogenetic analysis of OsMV2–4. (A) Genomic structures of OsMV2–4. ORFs are shown as boxes. (B) Evolutionary relationships of the viruses in the Mitoviridae family. RdRp amino acid sequences from different genera of Mitoviridae were aligned to derive a phylogenetic tree using the maximum likelihood method with the LG + I + G substitution model. Narnaviridae is set as the outgroup. Bootstrap values (%) of 1,000 replicates that are greater than 50% are labeled on the branches. The viruses OsMV2–4 detected in this study are marked with a red triangle.

Additionally, we identified another novel mitovirus in samples SRR5282569, SRR5428527, SRR8258357 and SRR9290661 and named it Ophiocordyceps sinensis mitovirus 4 (OsMV4). The viral length is 2,501 nt which encodes a 732-aa RdRp protein. The completeness of the terminal sequence of OsMV4 is currently uncertain. BLASTx analysis showed an identity of 60.9% between the RdRp of OsMV4 and that of Erysiphe necator-associated mitovirus 8 (query coverage: 85%; E-value: 0; and accession: QHD64819.1). According to the International Committee on Taxonomy of Viruses (ICTV) classification, there are four genera in the family Mitoviridae. Phylogenetic analysis suggested that OsMV1-4 belong to the genus Unuamitovirus (Figure 4B).

3.7 Four novel viruses in the family Narnaviridae

According to the ICTV classification, the Narnaviridae has only one genus, the narnavirus. A typical narnavirus lacks a capsid and encodes a single RdRp protein. We have identified four new narnaviruses in O. sinensis. All the viruses exhibit a single large ORF encoding an RdRp. They were designated as Ophiocordyceps sinensis narnavirus 1, 2, 3 and 4 (OsNV1, OsNV2, OsNV3, OsNV4). OsNV1 is comprised of 3,391 nt and encodes a 1,068-aa RdRp (Figure 5A). BLASTx analysis revealed that the RdRp of OsNV1 shares 54.1% identity with that of Guiyang Paspalum thunbergii narna-like virus 1 (query coverage: 93%; E-value: 0; and accession: UUW20993.1). OsNV2 has a full length of 2,975 nt and encodes a 964-aa RdRp protein showing the highest sequence identity of 39.0% with the RdRp of Downy mildew lesion associated orfanplasmovirus 1 (query coverage: 71%; E-value: 3e−133; and accession: QNQ74063.1). Interestingly, ambigrammatic narnaviruses can harbor an additional, reverse frame ORF (rORF) overlapping the RdRp ORF. The rORF spans almost the entire length of the viral genome (Retallack et al., 2021; Zhang et al., 2023). Instead of a large rORF, OsNV2 presents two smaller reverse-frame ORFs with sizes of 78-aa and 89-aa. OsNV3 and OsNV4, with full lengths of 2,174 nt and 2029 nt, respectively. OsNV3 encodes a 649-aa RdRp and shares 45.6% identity with the RdRp of Downy mildew lesion associated splipalmivirus 46 (query coverage: 88%; E-value: 0; and accession: WLJ60693.1), while OsNV4’s 597-aa RdRp shows 63.5% identity with Suillus luteus narnavirus 6 (query coverage: 86%; E-value: 0; and accession: WLK77414.1). The phylogenetic analysis supported the classification of OsNV1–4 as new species within the family Narnaviridae (Figure 5B).

Figure 5
www.frontiersin.org

Figure 5. Genome organization and phylogenetic analysis of OsNV1–4. (A) Genomic structures of OsNV1–4. ORFs are shown as boxes. (B) Evolutionary relationships of the viruses in the Narnaviridae family. RdRp amino acid sequences from different genera of Narnaviridae were aligned to derive a phylogenetic tree using the maximum likelihood method with the LG + I + G substitution model. Botourmiaviridae is set as the outgroup. Bootstrap values (%) of 1,000 replicates that are greater than 50% are labeled on the branches. The viruses OsNV1–4 are marked with a red triangle.

3.8 Partitiviridae-related sequences

A new dsRNA partitivirus, named Ophiocordyceps sinensis partitivirus 1 (OsPV1), was discovered in sample SRR5428527. This marks the first report of a partitivirus in O. sinensis. Typically, partitiviruses have two dsRNA genomic segments, one of which encodes an RdRp and the other encodes a CP (Zhu et al., 2024). However, only one contig of 2,152 nt was found in OsPV1 (Figure 6A). It encodes an incomplete RdRp sequence with 76.9% identity to the RdRp in Rosellinia necatrix partitivirus 15 (query coverage: 99%; E-value: 0; and accession: BBU59838.1). The segment encoding the CP was not detected in OsPV1. Based on phylogenetic analysis, OsPV1 showed a close evolutionary relationship with members of the genus Betapartitivirus (Figure 6B).

Figure 6
www.frontiersin.org

Figure 6. Genome organization and phylogenetic analysis of OsPV1. (A) Genomic structure of OsPV1. ORF is shown as a box. The sequence of the genome is incomplete. (B) Evolutionary relationships of the viruses in the Partitiviridae family. RdRp nucleotide sequences from different genera of Partitiviridae were aligned to derive a phylogenetic tree using the maximum likelihood method with the GTR + I + G substitution model. Polymycoviridae is set as the outgroup. Bootstrap values (%) of 1,000 replicates that are greater than 50% are labeled on the branches. The virus OsPV1 is marked with a red triangle.

3.9 Vivivirus-related sequences

A new vivivirus, named Ophiocordyceps sinensis vivivirus 1 (OsVV1), was identified in this study. It is the first time for vivivirus identified among entomopathogenic fungi. Viviviruses, related phylogenetically to plant virgaviruses, have been discovered in phytopathogenic fungi. Vivivirus typically possesses two to three (+) ssRNA genomic segments, one of which expresses methyltransferase (Mtr) and helicase (Hel), the other expresses RdRp and sometimes expresses another Mtr (Chiapello et al., 2020; Degola et al., 2021; Kondo et al., 2022). Sample SRR8258357 contained two contigs, with lengths of 1,276 nt and 2,482 nt, respectively (Figure 7A). BLASTx analysis revealed that the shorter contig had 52.4% identity to RdRp in Sisal-associated virgavirus A (query coverage: 93%; E-value: 1e−122; and accession: QYJ09848.1). Phylogenetic analysis showed a close relationship between OsVV1 and other viviviruses, supporting its identification (Figure 7B). The other contig showed 27.9% identity to a hypothetical protein of Aspergillus flavus vivivirus 1 (query coverage: 56%; E-value: 6e−20; and accession: BED98313.1). In addition, the structural domain encoding Mtr was not found.

Figure 7
www.frontiersin.org

Figure 7. Genome organization and phylogenetic analysis of OsVV1. (A) Genome structures of OsVV1. ORFs are shown as boxes. The sequences of the genomes are all incomplete. (B) Evolutionary relationships of the viruses in the viviviruses. RdRp nucleotide sequences from different viviviruses were aligned to derive a phylogenetic tree using the maximum likelihood method with the GTR + I + G substitution model. Closteroviridae is set as the outgroup. Bootstrap values (%) of 1,000 replicates that are greater than 50% are labeled on the branches. The virus OsVV1 is marked with a red triangle.

3.10 Deltaflexiviridae-related sequences

A contig of 7,877 nt from samples SRR12952889 and SRR8258357 showed 56.9% identity to the RdRp of Erysiphe necator-associated flexivirus 1 (EnaFV1) (query coverage: 69%; E-value: 0; and accession: QKN22686.1). This contig contains a large ORF encoding 1926 aa, which was predicted to present three structural domains: Mtr, Hel1 and RdRP (Figure 8A). Flexivirues can be classified into five families: Alphaflexiviridae, Betaflexiviridae, Deltaflexiviridae, Gammaflexiviridae and Tymoviridae (Ye et al., 2023). Based on phylogenetic analysis, this novel flexivirus exhibited a close evolutionary relationship with members in the Deltaflexiviridae family (Figure 8B) and thus was designated as Ophiocordyceps sinensis deltaflexivirus 1 (OsDFV1).

Figure 8
www.frontiersin.org

Figure 8. Genome organization and phylogenetic analysis of OsDFV1. (A) Genome structures of OsDFV1. ORFs are shown as boxes. The sequence of the genome is incomplete. (B) Evolutionary relationships of the viruses in the Tymovirales. RdRp amino acid sequences from different families of Tymovirales were aligned to derive a phylogenetic tree using the maximum likelihood method with the LG + I + G substitution model. Narnaviridae is set as the outgroup. Bootstrap values (%) of 1,000 replicates that are greater than 50% are labeled on the branches. The virus OsDFV1 is marked with a red triangle.

3.11 Divergent viruses were detected in the data from Chinese cordyceps at different developmental stages

Samples of three distinct O. sinensis strains were analyzed across their various developmental stages. It is shown that distinct viral species and loads were observed in the data from the same strain at various points during the fungal lifecycle (Table 3 and Supplementary Table S5).

Table 3
www.frontiersin.org

Table 3. Summary of Chinese cordyceps mycovirus genomes at different developmental stages identified from transcriptomic datasets.

In PRJNA507459 (Supplementary Table S6), Chinese cordyceps was classified into six distinct developmental stages. The hyphae sample exhibited the highest viral count. Samples from young fruitbody, developed fruitbody and mature fruitbody, were detected with an equal number of viruses, hosted different viral species. In PRJNA673413 (Supplementary Table S7), Chinese cordyceps was divided into three developmental stages. The highest number of viruses was detected in the larva-fungus complex sample, suggesting a possible influence of the interaction between fungus and host on viral diversity or load. Interestingly, none of the 13 viruses mentioned above was identified in the fruitbody sample, highlighting the possibility of stage-specific viral resistance or clearance mechanisms in this strain. PRJNA625214 (Supplementary Table S8) and PRJNA600609 (Supplementary Table S9) represent the same strain with a total of seven distinct developmental stages. Viruses were detected in all the stages, with the highest virus counts being in the blastospores in stationary stage, containing five distinct species. Notably, OsOMV1 is consistently present in six developmental stages examined. Overall, the number of viruses carried did not favor any specific developmental stage.

4 Discussion

The study has isolated and sequenced dsRNA in O. sinensis strains. By combining the analysis of publicly available transcriptome data from the SRA, we provide a comprehensive report on the diversity of mycovirus within O. sinensis. We identified 13 distinct mycoviruses, obtaining the full genomes of eight through RACE, while assembling the remaining five from public data, albeit with incomplete sequences. Most of these viruses possess a (+) ssRNA genome, except for one dsRNA virus. These mycoviruses span five families (Partitiviridae, Mitoviridae, Narnaviridae, Botourmiaviridae, Deltaflexiviridae) as well as two unclassified lineages (Ormycovirus and Vivivirus). Notably, the Mitoviridae and Narnaviridae families were the most prevalent (Table 4), suggesting their significant ecological or evolutionary roles in O. sinensis.

Table 4
www.frontiersin.org

Table 4. Summary of mycoviruses in this research.

Mitoviridae viruses, typically localizing within mitochondria, provide stability and resistance against the host’s detoxification, and may evade host antiviral RNA silencing mechanisms (Shahi et al., 2019). Currently, the Mitoviridae family is divided into four distinct genera. Our phylogenetic analysis suggests OsMV1-4 all belong to the genus Unuamitovirus. The Narnaviridae family exhibits remarkable diversity in terms of genome structure and host range (Hillman and Cai, 2013). OsNV1–4, detected in this study, are typical of narnaviruses, possessing an unsegmented, linear (+) ssRNA genome without rORFs. The simplicity of these viral genome structures could indicate a streamlined replication strategy that may confer evolutionary advantages. A thorough search for ORFans led to the discovery of a unique group of highly differentiated RNA arboviruses, designated as ormycovirus. SbOMV1 is currently the only ormycovirus collected from a corresponding live biological sample. OsOMV1 was identified earlier using publicly accessible data (Forgia et al., 2022). We detected OsOMV1 in a dsRNA sample from the IOZ strain, offering compelling evidence for the existence of ormycovirus. Notably, this virus was abundantly present in the SRA samples analyzed, highlighting its potential significance in the Chinese cordyceps lifecycle.

Coinfection of mycoviruses is a common phenomenon. The likelihood of coinfection in nature is much greater than that of single mycovirus infections (Thapa and Roossinck, 2019). For instance, Sclerotinia sclerotiorum SX276 was infected by nine mycoviruses from different evolutionary lineages, including seven (+) ssRNA viruses, one dsRNA virus, and one negative-sense ssRNA virus. This coinfection resulted in an altered growth rate, abnormal colony morphology, and reduced pathogenicity, characteristics typical of a hypovirulent strain (Mu et al., 2021). B. bassiana isolates infected with Beauveria bassiana polymycovirus 1 (BbPmV-1) and Beauveria bassiana non-segmented virus 1 (BbNV-1) exhibited increased growth rate, biomass, and hypervirulence towards Galleria mellonella larvae (Kotta-Loizou and Coutts, 2017). According to our observation, O. sinensis infecting its host larvae differs from the infection patterns of entomopathogenic fungi, such as B. bassiana and Metarhizium robertsii. The fungus can evade host immune responses and sustain a prolonged infection (5 to 12 months) without inducing disease symptoms on the host body surface (Meng et al., 2021; Wu et al., 2022b). Among the IOZ strains, IOZ14 and IOZ29 harbored eight viruses, while the other strains contained seven. Except for SRR13286702 and SRR5446809, all other SRA samples displayed coinfections. Notably, data from the larva-fungus complex showed the highest viral count compared to data from other tissues in PRJNA673413. The prevalence of coinfection within O. sinensis suggests a complex viral ecology that may influence the fungal ability to infect and persist in the host over the long-term infection. According to our observations on the mycelia, no visible morphological differences among the six IOZ strains were noted, except that the colonies of strain IOZ31 were darker in color. Due to the diversity of mycoviruses present in the O. sinensis strains, it was challenging to determine their specific effects on the biological phenotypes of these strains. Isolating stains with a single mycovirus would help us unravel these complex mechanisms.

Chinese cordyceps has attracted significant attention due to its limited distribution range, mysterious lifecycle, ecological importance, and developmental biology. Based on the process of artificial cultivation of O. sinensis, the development of its fruiting body can be divided into different stages. During the analysis of the SRA samples, it was discovered that the data from the same strain showed various amounts and types of mycoviruses at different developmental stages (Table 3), suggesting a dynamic shift in viral composition during its development. However, the viral distribution did not show a clear preference. This might be because of their geographically disparate habits or the limited number of analyzed samples. OsOMV1 was identified in SRR10878121 sample but not in the SRR8258346 sample. Both samples were collected at the primordium stage. The virus could be consistently detected in all other developmental stage samples belonging to the same strain as SRR8258346. Viruses that ‘disappear’ within a strain at a certain stage were probably due to their low abundances and insufficient sequencing depth. Investigating the mycoviral abundance and distribution within more O. sinensis strains during the entire lifecycle may reveal virus-host dynamics and their ecological consequences.

5 Conclusion

In summary, this study not only maps the mycoviral landscape of O. sinensis but also opens new avenues for research into the ecological and evolutionary implications of mycoviruses in fungi. The insights gained here pave the way for studies exploring the intricate relationships between mycoviruses and their hosts, with potential applications in fungal biology, ecology, and biotechnology.

Data availability statement

The viral genomes characterized for this study can be found in the NCBI with the accession numbers listed in Table 4 for each viral genome segment. The IOZ strains transcriptome raw reads have been deposited in the SRA with accession number PRJNA1169795.

Author contributions

QK: Data curation, Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing. JZ: Investigation, Methodology, Supervision, Writing – review & editing. FC: Data curation, Formal analysis, Methodology, Writing – review & editing. CD: Resources, Writing – review & editing. QQ: Conceptualization, Supervision, Writing – review & editing. XL: Project administration, Writing – review & editing. HW: Investigation, Resources, Writing – review & editing. HZ: Conceptualization, Supervision, Writing – review & editing. QM: Conceptualization, Supervision, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was financially supported by Beijing Yun’an Bio-tech Co. Ltd.

Conflict of interest

JZ, QQ, XL, HW, HZ, and QM are employed by Beijing Yun’an Bio-tech Co. Ltd.

The remaining 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.

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.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2024.1493365/full#supplementary-material

Supplementary Table S1 | Primer pairs used to confirm viral sequences in IOZ strains.

Supplementary Table S2 | The terminal primer of the eight viruses in IOZ strains with RACE.

Supplementary Table S3 | Detailed information on contigs obtained from 13 different samples by analyzing the SRA database.

Supplementary Table S4 | Coverage (% genome length) of each mycovirus identified in the public transcriptome of different O. sinensis strains.

Supplementary Table S5 | Coverage (% genome length) of each mycovirus identified in the public transcriptome of Chinese cordyceps at different developmental stages.

Supplementary Table S6 | Detailed information on contigs obtained from 6 different samples by analyzing the PRJNA507459.

Supplementary Table S7 | Detailed information on contigs obtained from 3 different samples by analyzing the PRJNA673413.

Supplementary Table S8 | Detailed information on contigs obtained from 5 different samples by analyzing the PRJNA625214.

Supplementary Table S9 | Detailed information on contigs obtained from 2 different samples by analyzing the PRJNA600609.

Footnotes

References

Ayllón, M. A., Turina, M., Xie, J., Nerva, L., Marzano, S. L., Donaire, L., et al. (2020). Ictv virus taxonomy profile: botourmiaviridae. J. Gen. Virol. 101, 454–455. doi: 10.1099/jgv.0.001409

PubMed Abstract | Crossref Full Text | Google Scholar

Baral, B. (2017). Entomopathogenicity and biological attributes of himalayan treasured fungus Ophiocordyceps sinensis (yarsagumba). J. Fungi 3:4. doi: 10.3390/jof3010004

PubMed Abstract | Crossref Full Text | Google Scholar

Begeman, A., Babaian, A., and Lewis, S. C. (2023). Metatranscriptomic analysis uncovers prevalent viral orfs compatible with mitochondrial translation. Msystems 8:e0100222. doi: 10.1128/msystems.01002-22

PubMed Abstract | Crossref Full Text | Google Scholar

Capella-Gutiérrez, S., Silla-Martínez, J. M., and Gabaldón, T. (2009). Trimal: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 25, 1972–1973. doi: 10.1093/bioinformatics/btp348

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, S., Zhou, Y., Chen, Y., and Gu, J. (2018). Fastp: an ultra-fast all-in-one fastq preprocessor. Bioinformatics 34, i884–i890. doi: 10.1093/bioinformatics/bty560

PubMed Abstract | Crossref Full Text | Google Scholar

Chiapello, M., Rodriguez-Romero, J., Ayllon, M. A., and Turina, M. (2020). Analysis of the virome associated to grapevine downy mildew lesions reveals new mycovirus lineages. Virus Evol. 6:veaa0582. doi: 10.1093/ve/veaa058

PubMed Abstract | Crossref Full Text | Google Scholar

Danecek, P., Bonfield, J. K., Liddle, J., Marshall, J., Ohan, V., Pollard, M. O., et al. (2021). Twelve years of samtools and bcftools. Gigascience 10:giab0082. doi: 10.1093/gigascience/giab008

PubMed Abstract | Crossref Full Text | Google Scholar

Degola, F., Spadola, G., Forgia, M., Turina, M., Dramis, L., Chitarra, W., et al. (2021). Aspergillus goes viral: ecological insights from the geographical distribution of the mycovirome within an Aspergillus flavus population and its possible correlation with aflatoxin biosynthesis. J. Fungi 7:833. doi: 10.3390/jof7100833

PubMed Abstract | Crossref Full Text | Google Scholar

Edgar, R. C. (2022). Muscle5: high-accuracy alignment ensembles enable unbiased assessments of sequence homology and phylogeny. Nat. Commun. 13:6968. doi: 10.1038/s41467-022-34630-w

PubMed Abstract | Crossref Full Text | Google Scholar

Fan, B., and He, J. (2024). Mapping the Ophiocordyceps sinensis value chain: actors, profits and social institutions in south-West China. Oryx 58, 38–47. doi: 10.1017/S0030605323001588

Crossref Full Text | Google Scholar

Filippou, C., Garrido-Jurado, I., Meyling, N. V., Quesada-Moraga, E., Coutts, R. H. A., and Kotta-Loizou, I. (2018). Mycoviral population dynamics in Spanish isolates of the entomopathogenic fungus Beauveria bassiana. Viruses 10:665. doi: 10.3390/v10120665

PubMed Abstract | Crossref Full Text | Google Scholar

Forgia, M., Chiapello, M., Daghino, S., Pacifico, D., Crucitti, D., Oliva, D., et al. (2022). Three new clades of putative viral RNA dependent RNA polymerases with rare or unique catalytic triads discovered in libraries of orfans from powdery mildews and the yeast of oenological interest Starmerella bacillaris. Virus Evol. 8:veac0381. doi: 10.1093/ve/veac038

PubMed Abstract | Crossref Full Text | Google Scholar

Gilbert, K. B., Holcomb, E. E., Allscheid, R. L., and Carrington, J. C. (2019). Hiding in plain sight: new virus genomes discovered via a systematic analysis of fungal public transcriptomes. PLoS One 14:e0219207. doi: 10.1371/journal.pone.0219207

PubMed Abstract | Crossref Full Text | Google Scholar

Guo, J., Zhang, P., Wu, N., Liu, W., Liu, Y., Jin, H., et al. (2024). Transfection of entomopathogenic Metarhizium species with a mycovirus confers hypervirulence against two lepidopteran pests. Proc. Natl. Acad. Sci. 121:e2320572121. doi: 10.1073/pnas.2320572121

PubMed Abstract | Crossref Full Text | Google Scholar

Hillman, B. I., and Cai, G. (2013). The family narnaviridae: simplest of RNA viruses. Adv. Virus Res. 86, 149–176. doi: 10.1016/B978-0-12-394315-6.00006-4

Crossref Full Text | Google Scholar

Hollings, M. (1962). Viruses associated with a die-back disease of cultivated mushroom. Nature 196, 962–965. doi: 10.1038/196962a0

Crossref Full Text | Google Scholar

Hough, B., Steenkamp, E., Wingfield, B., and Read, D. (2023). Fungal viruses unveiled: a comprehensive review of mycoviruses. Viruses 15:1202. doi: 10.3390/v15051202

PubMed Abstract | Crossref Full Text | Google Scholar

Hwang, D. Y., Kim, S., Woo, S. D., Shin, T. Y., Coutts, R., and Kotta-Loizou, I. (2023). Incidence of putative RNA mycoviruses in entomopathogenic fungi in Korea. Arch. Virol. 168:145. doi: 10.1007/s00705-023-05765-y

PubMed Abstract | Crossref Full Text | Google Scholar

Kim, D., Paggi, J. M., Park, C., Bennett, C., and Salzberg, S. L. (2019). Graph-based genome alignment and genotyping with hisat2 and hisat-genotype. Nat. Biotechnol. 37, 907–915. doi: 10.1038/s41587-019-0201-4

PubMed Abstract | Crossref Full Text | Google Scholar

Kondo, H., Botella, L., and Suzuki, N. (2022). Mycovirus diversity and evolution revealed/inferred from recent studies. Annu. Rev. Phytopathol. 60, 307–336. doi: 10.1146/annurev-phyto-021621-122122

PubMed Abstract | Crossref Full Text | Google Scholar

Kotta-Loizou, I. (2021). Mycoviruses and their role in fungal pathogenesis. Curr. Opin. Microbiol. 63, 10–18. doi: 10.1016/j.mib.2021.05.007

Crossref Full Text | Google Scholar

Kotta-Loizou, I., and Coutts, R. H. (2017). Studies on the virome of the entomopathogenic fungus Beauveria bassiana reveal novel dsRNA elements and mild hypervirulence. PLoS Pathog. 13:e1006183. doi: 10.1371/journal.ppat.1006183

PubMed Abstract | Crossref Full Text | Google Scholar

Lambden, P. R., Cooke, S. J., Caul, O., and Clarke, I. N. (1992). Cloning of noncultivatable human rotavirus by single primer amplification. J. Virol. 66, 1817–1822. doi: 10.1128/jvi.66.3.1817-1822.1992

PubMed Abstract | Crossref Full Text | Google Scholar

Li, Y., Hsiang, T., Yang, R., Hu, X., Wang, K., Wang, W., et al. (2016). Comparison of different sequencing and assembly strategies for a repeat-rich fungal genome, Ophiocordyceps sinensis. J. Microbiol. Methods 128, 1–6. doi: 10.1016/j.mimet.2016.06.025

PubMed Abstract | Crossref Full Text | Google Scholar

Li, D., Luo, R., Liu, C., Leung, C., Ting, H., Sadakane, K., et al. (2016). Megahit v1.0: a fast and scalable metagenome assembler driven by advanced methodologies and community practices. Methods. 102, 3–11. doi: 10.1016/j.ymeth.2016.02.020

PubMed Abstract | Crossref Full Text | Google Scholar

Li, M., Meng, Q., Zhang, H., Shu, R., Zhao, Y., Wu, P., et al. (2020). Changes in transcriptomic and metabolomic profiles of morphotypes of Ophiocordyceps sinensis within the hemocoel of its host larvae, Thitarodes xiaojinensis. BMC Genomics 21:789. doi: 10.1186/s12864-020-07209-2

PubMed Abstract | Crossref Full Text | Google Scholar

Li, X., Wang, F., Liu, Q., Li, Q., Qian, Z., Zhang, X., et al. (2019). Developmental transcriptomics of Chinese cordyceps reveals gene regulatory network and expression profiles of sexual development-related genes. BMC Genomics 20:337. doi: 10.1186/s12864-019-5708-z

PubMed Abstract | Crossref Full Text | Google Scholar

Meng, Q., Wu, P. P., Li, M. M., Shu, R. H., Zhou, G. L., Zhang, J. H., et al. (2021). Distinct responses of Thitarodes xiaojinensis beta-1,3-glucan recognition protein-1 and immulectin-8 to Ophiocordyceps sinensis and Cordyceps militaris infection. J. Immunol. 207, 200–209. doi: 10.4049/jimmunol.2000447

PubMed Abstract | Crossref Full Text | Google Scholar

Minh, B. Q., Schmidt, H. A., Chernomor, O., Schrempf, D., Woodhams, M. D., von Haeseler, A., et al. (2020). Iq-tree 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol. Biol. Evol. 37, 1530–1534. doi: 10.1093/molbev/msaa015

Crossref Full Text | Google Scholar

Morris, T. J., and Dodds, J. A. (1979). Isolation and analysis of double-stranded-rna from virus-infected plant and fungal tissue. Phytopathology. 69, 854–858. doi: 10.1094/Phyto-69-854

Crossref Full Text | Google Scholar

Mu, F., Li, B., Cheng, S., Jia, J., Jiang, D., Fu, Y., et al. (2021). Nine viruses from eight lineages exhibiting new evolutionary modes that co-infect a hypovirulent phytopathogenic fungus. PLoS Pathog. 17:e1009823. doi: 10.1371/journal.ppat.1009823

PubMed Abstract | Crossref Full Text | Google Scholar

Pagnoni, S., Oufensou, S., Balmas, V., Bulgari, D., Gobbi, E., Forgia, M., et al. (2023). A collection of trichoderma isolates from natural environments in Sardinia reveals a complex virome that includes negative-sense fungal viruses with unprecedented genome organizations. Virus Evol 9:vead0422. doi: 10.1093/ve/vead042

PubMed Abstract | Crossref Full Text | Google Scholar

Qin, Q. L., Zhou, G. L., Zhang, H., Meng, Q., Zhang, J. H., Wang, H. T., et al. (2018). Obstacles and approaches in artificial cultivation of Chinese cordyceps. Mycology 9, 7–9. doi: 10.1080/21501203.2018.1442132

PubMed Abstract | Crossref Full Text | Google Scholar

Retallack, H., Popova, K. D., Laurie, M. T., Sunshine, S., and DeRisi, J. L. (2021). Persistence of ambigrammatic narnaviruses requires translation of the reverse open reading frame. J. Virol. 95:e0010921. doi: 10.1128/JVI.00109-21

PubMed Abstract | Crossref Full Text | Google Scholar

Sato, Y., and Suzuki, N. (2023). Continued mycovirus discovery expanding our understanding of virus lifestyles, symptom expression, and host defense. Curr. Opin. Microbiol. 75:102337. doi: 10.1016/j.mib.2023.102337

PubMed Abstract | Crossref Full Text | Google Scholar

Shahi, S., Eusebio-Cope, A., Kondo, H., Hillman, B. I., and Suzuki, N. (2019). Investigation of host range of and host defense against a mitochondrially replicating mitovirus. J. Virol. 93, e01503–186. doi: 10.1128/JVI.01503-18

PubMed Abstract | Crossref Full Text | Google Scholar

Shi, N., Zhu, Q., Yang, G., Wang, P., and Huang, B. (2024). Prevalence and species diversity of dsRNA mycoviruses from Beauveria bassiana strains in the China's guniujiang nature. Heliyon 10:e30186. doi: 10.1016/j.heliyon.2024.e30186

PubMed Abstract | Crossref Full Text | Google Scholar

Tamura, K., Stecher, G., Kumar, S., and Battistuzzi, F. U. (2021). Mega11: molecular evolutionary genetics analysis version 11. Mol. Biol. Evol. 38, 3022–3027. doi: 10.1093/molbev/msab120

PubMed Abstract | Crossref Full Text | Google Scholar

Tang, C. Y., Wang, J., Liu, X., Chen, J. B., Liang, J., Wang, T., et al. (2022). Medium optimization for high mycelial soluble protein content of Ophiocordyceps sinensis using response surface methodology. Front. Microbiol. 13:1055055. doi: 10.3389/fmicb.2022.1055055

PubMed Abstract | Crossref Full Text | Google Scholar

Thapa, V., and Roossinck, M. J. (2019). Determinants of coinfection in the mycoviruses. Front. Cell. Infect. Microbiol. 9:169. doi: 10.3389/fcimb.2019.00169

Crossref Full Text | Google Scholar

Tong, X., Zhang, H., Wang, F., Xue, Z., Cao, J., Peng, C., et al. (2020). Comparative transcriptome analysis revealed genes involved in the fruiting body development of Ophiocordyceps sinensis. PeerJ 8:e8379. doi: 10.7717/peerj.8379

PubMed Abstract | Crossref Full Text | Google Scholar

Villan Larios, D. C., Diaz Reyes, B. M., Pirovani, C. P., Loguercio, L. L., Santos, V. C., Góes-Neto, A., et al. (2023). Exploring the mycovirus universe: identification, diversity, and biotechnological applications. J. Fungi 9:361. doi: 10.3390/jof9030361

PubMed Abstract | Crossref Full Text | Google Scholar

Wu, P., Qin, Q., Zhang, J., Zhang, H., Li, X., Wang, H., et al. (2022a). The invasion process of the entomopathogenic fungus Ophiocordyceps sinensis into the larvae of ghost moths (Thitarodes xiaojinensis) using a GFP-labeled strain. Front. Microbiol. 13:974323. doi: 10.3389/fmicb.2022.974323

PubMed Abstract | Crossref Full Text | Google Scholar

Wu, P., Shu, R., Gao, X., Li, M., Zhang, J., Zhang, H., et al. (2022b). Immulectin-2 from the ghost moth, Thitarodes xiaojinensis (Lepidoptera: Hepialidae), modulates cellular and humoral responses against fungal infection. Dev. Comp. Immunol. 133:104429. doi: 10.1016/j.dci.2022.104429

PubMed Abstract | Crossref Full Text | Google Scholar

Xia, E., Yang, D., Jiang, J., Zhang, Q., Liu, Y., Liu, Y., et al. (2017). The caterpillar fungus, Ophiocordyceps sinensis, genome provides insights into highland adaptation of fungal pathogenicity. Sci. Rep. 7:1806. doi: 10.1038/s41598-017-01869-z

PubMed Abstract | Crossref Full Text | Google Scholar

Ye, T., Lu, Z., Li, H., Duan, J., Hai, D., Lin, Y., et al. (2023). Characterization of a fungal virus representing a novel genus in the family alphaflexiviridae. Viruses 15:3392. doi: 10.3390/v15020339

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, Y., Li, E., Wang, C., Li, Y., and Liu, X. (2012). Ophiocordyceps sinensis, the flagship fungus of China: terminology, life strategy and ecology. Mycology 3, 2–10. doi: 10.1080/21501203.2011.654354

Crossref Full Text | Google Scholar

Zhang, Y., Liang, X., Zhao, M., Qi, T., Guo, H., Zhao, J., et al. (2023). A novel ambigrammatic mycovirus, psv5, works hand in glove with wheat stripe rust fungus to facilitate infection. Plant Commun. 4:100505. doi: 10.1016/j.xplc.2022.100505

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, Y., Shi, N., Wang, P., Zhu, Q., Yang, G., and Huang, B. (2022). Molecular characterization of a novel alternavirus infecting the entomopathogenic fungus Cordyceps chanhua. Arch. Virol. 167, 1467–1470. doi: 10.1007/s00705-022-05446-2

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, H., Yue, P., Tong, X., Bai, J., Yang, J., and Guo, J. (2021). MRNA-seq and miRNA-seq profiling analyses reveal molecular mechanisms regulating induction of fruiting body in Ophiocordyceps sinensis. Sci. Rep. 11:12944. doi: 10.1038/s41598-021-91718-x

PubMed Abstract | Crossref Full Text | Google Scholar

Zhao, Y., Zhang, J., Meng, Q., Zhang, H., Zhou, G., Li, M., et al. (2020). Transcriptomic analysis of the orchestrated molecular mechanisms underlying fruiting body initiation in Chinese cordyceps. Gene 763:145061. doi: 10.1016/j.gene.2020.145061

Crossref Full Text | Google Scholar

Zhu, J. Z., Qiu, Z. L., Gao, B. D., Li, X. G., and Zhong, J. (2024). A novel partitivirus conferring hypovirulence by affecting vesicle transport in the fungus colletotrichum. MBio 15:e0253023. doi: 10.1128/mbio.02530-23

PubMed Abstract | Crossref Full Text | Google Scholar

Zhu, Q., Shi, N., Wang, P., Zhang, Y., Peng, F., Yang, G., et al. (2022a). A novel gammapartitivirus that causes changes in fungal development and multi-stress tolerance to important medicinal fungus Cordyceps chanhua. J. Fungi 8:1309. doi: 10.3390/jof8121309

PubMed Abstract | Crossref Full Text | Google Scholar

Zhu, Q., Shi, N., Zhang, Y., Peng, F., Yang, G., and Huang, B. (2022b). Complete genome sequence of a novel victorivirus infecting cicada flower (Cordyceps chanhua). Arch. Virol. 168:4. doi: 10.1007/s00705-022-05640-2

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: Ophiocordyceps sinensis, Mycovirus, transcriptome, public transcriptome mining, viral taxonomy

Citation: Kang Q, Zhang J, Chen F, Dong C, Qin Q, Li X, Wang H, Zhang H and Meng Q (2024) Unveiling mycoviral diversity in Ophiocordyceps sinensis through transcriptome analyses. Front. Microbiol. 15:1493365. doi: 10.3389/fmicb.2024.1493365

Received: 09 September 2024; Accepted: 05 November 2024;
Published: 25 November 2024.

Edited by:

Jiatao Xie, Huazhong Agricultural University, China

Reviewed by:

Ioly Kotta-Loizou, Imperial College London, United Kingdom
Shin-Yi Lee Marzano, United States Department of Agriculture, United States

Copyright © 2024 Kang, Zhang, Chen, Dong, Qin, Li, Wang, Zhang and Meng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Huan Zhang, emhhbmdodWFuQGlvei5hYy5jbg==; Qian Meng, bWVuZ3FpYW5AaW96LmFjLmNu

Disclaimer: 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.