- 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. 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).
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. 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. 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. 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. 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. 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. 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. 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. 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.
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
Baral, B. (2017). Entomopathogenicity and biological attributes of himalayan treasured fungus Ophiocordyceps sinensis (yarsagumba). J. Fungi 3:4. doi: 10.3390/jof3010004
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
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
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
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
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
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
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
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
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
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
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
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
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
Hollings, M. (1962). Viruses associated with a die-back disease of cultivated mushroom. Nature 196, 962–965. doi: 10.1038/196962a0
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
Edited by:
Jiatao Xie, Huazhong Agricultural University, ChinaReviewed by:
Ioly Kotta-Loizou, Imperial College London, United KingdomShin-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