Profile and functional analysis of small RNAs derived from Aspergillus fumigatus infected with double-stranded RNA mycoviruses

Recent developments in Aspergillus fumigatus research: diversity, drugs, and disease

SUMMARYAdvances in modern medical therapies for many previously intractable human diseases have improved patient outcomes. However, successful disease treatment outcomes are often prevented due to invasive fungal infections caused by the environmental mold Aspergillus fumigatus. As contemporary antifungal therapies have not experienced the same robust advances as other medical therapies, defining mechanisms of A. fumigatus disease initiation and progression remains a critical research priority. To this end, the World Health Organization recently identified A. fumigatus as a research priority human fungal pathogen and the Centers for Disease Control has highlighted the emergence of triazole-resistant A. fumigatus isolates. The expansion in the diversity of host populations susceptible to aspergillosis and the complex and dynamic A. fumigatus genotypic and phenotypic diversity call for a reinvigorated assessment of aspergillosis pathobiological and drug-susceptibility mechanisms. Here, we summarize recent advancements in the field and discuss challenges in our understanding of A. fumigatus heterogeneity and its pathogenesis in diverse host populations.

Mycoviruses from Aspergillus fungi involved in fermentation of dried bonito

Fungi are exploited for fermentation of foods such as cheese, Japanese sake, and soy sauce. However, the diversity of viruses that infect fungi involved in food fermentation is poorly understood. Fermented dried bonito ("katsuobushi") is one of the most important processed marine products in Japan. Fungi involved in katsuobushi fermentation are called katsuobushi molds, and Aspergillus spp. have been reported to be dominant on the surface of katsuobushi during fermentation. Because various mycoviruses have been found in members of the genus Aspergillus, we hypothesized that katsuobushi molds are also infected with mycoviruses. Here, we describe seven novel mycoviruses belonging to six families (Chrysoviridae, Fusariviridae, Mitoviridae, Partitiviridae, Polymycoviridae, and Pseudototiviridae) from isolated katsuobushi molds (Aspergillus chevalieri and A. sulphureus) detected by fragmented and primer-ligated double-stranded RNA sequencing. Aspergillus chevalieri fusarivirus 1 has a unique bi-segmented genome, whereas other known fusariviruses have a single genomic segment. Phenotypic comparison between the parental A. chevalieri strain infected with Aspergillus chevalieri polymycovirus 1 (AchPmV1) and isogenic AchPmV1-free isolates indicated that AchPmV1 inhibits the early growth of the host. This study reveals the diversity of mycoviruses that infect katsuobushi molds, and provides insight into the effect of mycoviruses on fungi involved in fermentation.

Discovering the hidden function in fungal genomes

New molecular technologies have helped unveil previously unexplored facets of the genome beyond the canonical proteome, including microproteins and short ORFs, products of alternative splicing, regulatory non-coding RNAs, as well as transposable elements, cis-regulatory DNA, and other highly repetitive regions of DNA. In this Review, we highlight what is known about this 'hidden genome' within the fungal kingdom. Using well-established model systems as a contextual framework, we describe key elements of this hidden genome in diverse fungal species, and explore how these factors perform critical functions in regulating fungal metabolism, stress tolerance, and pathogenesis. Finally, we discuss new technologies that may be adapted to further characterize the hidden genome in fungi.

RNA interference of Aspergillus flavus in response to Aspergillus flavus partitivirus 1 infection

RNA interference (RNAi) is one of the important defense responses against viral infection, but its mechanism and impact remain unclear in mycovirus infections. In our study, reverse genetics and virus-derived small RNA sequencing were used to show the antiviral responses of RNAi components in Aspergillus flavus infected with Aspergillus flavus partitivirus 1 (AfPV1). qRT-PCR revealed that AfPV1 infection induced the expression of the RNAi components in A. flavus compared with noninfected A. flavus. Knock mutants of each RNAi component were generated, but the mutants did not exhibit any obvious phenotypic changes compared with the A. flavus parental strain. However, after AfPV1 inoculation, production of AfPV1 was significantly less than in the parental strain. Furthermore, sporulation was greater in each AfPV1-infected mutant compared with the AfPV1-infected parental A. flavus. We also investigated the sensitivity of virus-free and AfPV1-infected RNAi mutants and the parental strain to cell wall stress, osmotic stress, genotoxic stress, and oxidative stress. The mutants of DCLs and AGOs infected by AfPV1 displayed more changes than RDRP mutants in response to the first three stresses. Small RNA sequencing analysis suggested that AfPV1 infection reduced the number of unique reads of sRNA in A. flavus, although there were many vsiRNA derived from the AfPV1 genome. GO term and KEGG pathway analyses revealed that the functions of sRNA affected by AfPV1 infection were closely related to vacuole production. These results provide a better understanding of the functional role of RNAi in the impact of AfPV1 on the hypovirulence of A. flavus.

Experimental verification of strain-dependent relationship between mycovirus and its fungal host

Mycoviruses are viruses that infect fungi. Unlike mammalian infectious viruses, their life cycle does not generally have an extracellular stage, and a symbiosis-like relationship is maintained between virus and host fungi. Recently, mycoviruses have been reported to show effects on host fungi, altering biological properties such as growth rate, virulence, drug resistance, and metabolite production. In this study, we systematically elucidated the effects of viruses on host cells by comparing host phenotypes and transcriptomic responses in multiple sets of virus-infected and -eliminated Aspergillus flavus strains. The comparative study showed that mycoviruses affect several cellular activities at the molecular level in a virus- and host strain-dependent manner. The virus-swapping experiment revealed that difference with only three bases in the virus genome led to different host fungal response at the transcriptional level. Our results highlighted highly specific relationship between viruses and their host fungi.

The good, the bad and the cryptic: The multifaceted roles of mycoviruses and their potential applications for a sustainable agriculture

Mycoviruses are natural inhabitants of fungi and have been identified in almost all fungal taxonomic groups. Mycoviruses that infect phytopathogenic fungi are now becoming a hot research area due to their potential for the biocontrol of important plant pathogens. But, before considering a mycovirus for biocontrol, we should be fully aware of the effects it induces in a fungal host and its interactions with other viruses, fungal strains and even the host plants. Mycoviral infections are generally associated with different effects, ranging from hypovirulence to hypervirulence, but they can often be cryptic (latent infections). The cryptic lifestyle has been associated to many mycoviruses, but thanks to growing knowledge we are now aware that it is often associated to axenic conditions while the real effects can be observed only in nature. Other mycoviruses either promote (hypervirulence) or (hypovirulence) fungal pathogenicity by a strong impact on the fungal physiology or by blocking the production of toxins or effectors. Finally, indirect effects of mycoviral infections can also be provided to the plant that hosts the fungal isolate, highlighting not only their potential as direct biocontrol agents but also as priming agents for plant resilience to biotic and abiotic stresses. This review provides a broad overview of mycoviral interactions both with their hosts and with other mycoviruses, highlighting the most interesting examples. In contrast to what has been observed to date, we believe that the collective availability of these data will not only improve our understanding of mycoviruses, but also increase our confidence in considering them as alternative measures against fungal diseases to improve the sustainable production of food and feed commodities.

Disruption of the Aspergillus fumigatus RNA interference machinery alters the conidial transcriptome

The RNA interference (RNAi) pathway has evolved numerous functionalities in eukaryotes, with many on display in Kingdom Fungi. RNAi can regulate gene expression, facilitate drug resistance, or even be altogether lost to improve growth potential in some fungal pathogens. In the WHO fungal priority pathogen, Aspergillus fumigatus, the RNAi system is known to be intact and functional. To extend our limited understanding of A. fumigatus RNAi, we first investigated the genetic variation in RNAi-associated genes in a collection of 217 environmental and 83 clinical genomes, where we found that RNAi components are conserved even in clinical strains. Using endogenously expressed inverted-repeat transgenes complementary to a conditionally essential gene (pabA) or a nonessential gene (pksP), we determined that a subset of the RNAi componentry is active in inverted-repeat transgene silencing in conidia and mycelium. Analysis of mRNA-seq data from RNAi double-knockout strains linked the A. fumigatus dicer-like enzymes (DclA/B) and RNA-dependent RNA polymerases (RrpA/B) to regulation of conidial ribosome biogenesis genes; however, surprisingly few endogenous small RNAs were identified in conidia that could explain this broad change. Although RNAi was not clearly linked to growth or stress response defects in the RNAi knockouts, serial passaging of RNAi knockout strains for six generations resulted in lineages with diminished spore production over time, indicating that loss of RNAi can exert a fitness cost on the fungus. Cumulatively, A. fumigatus RNAi appears to play an active role in defense against double-stranded RNA species alongside a previously unappreciated housekeeping function in regulation of conidial ribosomal biogenesis genes.

Exploring the Mycovirus Universe: Identification, Diversity, and Biotechnological Applications

Viruses that infect fungi are known as mycoviruses and are characterized by the lack of an extracellular phase. In recent years, the advances on nucleic acids sequencing technologies have led to a considerable increase in the number of fungi-infecting viral species described in the literature, with a special interest in assessing potential applications as fungal biocontrol agents. In the present study, we performed a comprehensive review using Scopus, Web of Science, and PubMed databases to mine mycoviruses data to explore their molecular features and their use in biotechnology. Our results showed the existence of 267 mycovirus species, of which 189 are recognized by the International Committee on Taxonomy of Viruses (ICTV). The majority of the mycoviruses identified have a dsRNA genome (38.6%), whereas the Botourmiaviridae (ssRNA+) alone represents 14% of all mycoviruses diversity. Regarding fungal hosts, members from the Sclerotinicaeae appeared as the most common species described to be infected by mycoviruses, with 16 different viral families identified so far. It is noteworthy that such results are directly associated with the high number of studies and strategies used to investigate the presence of viruses in members of the Sclerotinicaeae family. The knowledge about replication strategy and possible impact on fungi biology is available for only a small fraction of the mycoviruses studied, which is the main limitation for considering these elements potential targets for biotechnological applications. Altogether, our investigation allowed us to summarize the general characteristics of mycoviruses and their hosts, the consequences, and the implications of this knowledge on mycovirus–fungi interactions, providing an important source of information for future studies.

Interference of Small RNAs in Fusarium graminearum through FgGMTV1 Infection

Small RNA (sRNA) plays a central role in RNA silencing in fungi. The genome of Fusarium graminearum gemytripvirus 1 (FgGMTV1) is comprised of three DNA segments: DNA-A, DNA-B, and DNA-C. DNA-A and DNA-B are associated with fungal growth and virulence reduction. To elucidate the role of RNA silencing during the interactions of fungi and viruses, the sRNA profiles of F. graminearum in association with FgGMTV1 were established, using an FgGMTV1-free library (S-S), a library for infection with the DNA-A and DNA-B segments (S-AB), and a library for infection with the DNA-A, DNA-B, and DNA-C segments (S-ABC). A large amount of virus-derived sRNA (vsiRNA) was detected in the S-AB and S-ABC libraries, accounting for 9.9% and 13.8% of the total sRNA, respectively, indicating that FgGMTV1 triggers host RNA silencing. The total numbers of sRNA reads differed among the three libraries, suggesting that FgGMTV1 infection interferes with host RNA silencing. In addition, the relative proportions of the different sRNA lengths were altered in the S-AB and S-ABC libraries. The genome distribution patterns of the mapping of vsiRNA to DNA-A and DNA-B in the S-AB and S-ABC libraries were also different. These results suggest the influence of DNA-C on host RNA silencing. Transcripts targeted by vsiRNAs were enriched in pathways that included flavin adenine dinucleotide binding, protein folding, and filamentous growth.

Similar Characteristics of siRNAs of Plant Viruses Which Replicate in Plant and Fungal Hosts

Simple Summary

RNA silencing in fungi was shown to confer antiviral defense against plant viruses. In this study, using high-throughput sequencing and bioinformatic analyses, we showed that small interfering RNAs (siRNAs) of cucumber mosaic virus and tobacco mosaic virus (TMV) which replicated in phytopathogenic fungi Rhizoctonia solani and Fusarium graminearum had similarities with viral siRNAs produced in plant hosts in regard to the size distributions, proportion of plus and minus senses, and nucleotide preference for the 5′ termini. Additionally, our results also determined that both F. graminearum DCL1 and DCL2 were involved in the production of TMV siRNAs. Thus, the fungal RNA silencing machineries have adaptive capabilities to recognize and process the genome of invading plant viruses.

Abstract

RNA silencing is a host innate antiviral mechanism which acts via the synthesis of viral-derived small interfering RNAs (vsiRNAs). We have previously reported the infection of phytopathogenic fungi by plant viruses such as cucumber mosaic virus (CMV) and tobacco mosaic virus (TMV). Furthermore, fungal RNA silencing was shown to suppress plant virus accumulation, but the characteristics of plant vsiRNAs associated with the antiviral response in this nonconventional host remain unknown. Using high-throughput sequencing, we characterized vsiRNA profiles in two plant RNA virus–fungal host pathosystems: CMV infection in phytopathogenic fungus Rhizoctonia solani and TMV infection in phytopathogenic fungus Fusarium graminearum. The relative abundances of CMV and TMV siRNAs in the respective fungal hosts were much lower than those in the respective experimental plant hosts, Nicotiana benthamiana and Nicotiana tabacum. However, CMV and TMV siRNAs in fungi had similar characteristics to those in plants, particularly in their size distributions, proportion of plus and minus senses, and nucleotide preference for the 5′ termini of vsiRNAs. The abundance of TMV siRNAs largely decreased in F. graminearum mutants with a deletion in either dicer-like 1 (dcl1) or dcl2 genes which encode key proteins for the production of siRNAs and antiviral responses. However, deletion of both dcl1 and dcl2 restored TMV siRNA accumulation in F. graminearum, indicating the production of dcl-independent siRNAs with no antiviral function in the absence of the dcl1 and dcl2 genes. Our results suggest that fungal RNA silencing recognizes and processes the invading plant RNA virus genome in a similar way as in plants.

Role of MicroRNA-Like RNAs in the Regulation of Spore Morphological Differences in the Entomopathogenic Fungus Metarhizium acridum

Metarhizium acridum is an important microbial pesticide. Conidia (CO) and blastospores (BS) are two types of spores that occur in different patterns in the M. acridum life cycle and exhibit significant differences in cell morphology, structure, and activity. It may suggest that the fungus has a complex gene regulation mechanism. While previous studies on the differences between CO and BS have mainly focused on cell structure and application, little is known regarding the differences between CO and BS in fungi on the transcriptome levels. MicroRNAs (miRNAs) are small noncoding RNAs crucial to gene regulation and cell function. Understanding the miRNA-like RNAs (milRNA) and mRNA expression profiles related to cell growth and cellular morphological changes would elucidate the roles of miRNAs in spore morphological differences. In this study, 4,646 differentially expressed genes (DEGs) were identified and mainly classified in the GO terms cell, cell part, biological process, and catalytic activity. The KEGG annotation suggested that they were enriched in amino acid biosynthesis, carbohydrate metabolism, ribosome, and oxidative phosphorylation and might be involved in cell activity and structure. There were 113 differentially expressed milRNAs (DEMs), targeting 493 DEGs. Target gene functional analysis revealed that the target genes were mainly enriched in RNA transport, purine metabolism, and the cell cycle. In addition, we identified essential genes from milRNA-mRNA pairs that might participate in cell budding growth and cell membrane and wall integrity, including adenosine deaminase, glycosyl hydrolase, and G-patch domain protein (dno-miR-328-3p), WD repeat-containing protein pop1 (age-miR-127), and GPI-anchored wall transfer protein (cgr-miR-598). MilRNAs might therefore play a crucial role in cell growth and cellular morphological changes as transcriptional and post-transcriptional regulators.

A Botybirnavirus Isolated from Alternaria tenuissima Confers Hypervirulence and Decreased Sensitivity of Its Host Fungus to Difenoconazole

Alternaria alternata botybirnavirus 1 (AaBRV1) was isolated from a strain of Alternaria alternata, causing watermelon leaf blight in our previous research. The effect of AaBRV1 on the phenotype of its host fungus, however, was not determined. In the present study, a novel strain of AaBRV1 was identified in A. tenuissima strain TJ-NH-51S-4, the causal agent of cotton Alternaria leaf spot, and designated as AaBRV1-AT1. A mycovirus AaBRV1-AT1-free strain TJ-NH-51S-4-VF was obtained by protoplast regeneration, which eliminated AaBRV1-AT1 from the mycovirus AaBRV1-AT1-infected strain TJ-NH-51S-4. Colony growth rate, spore production, and virulence of strain TJ-NH-51S-4 were greater than they were in TJ-NH-51S-4-VF, while the sensitivity of strain TJ-NH-51S-4 to difenoconazole, as measured by the EC50, was lower. AaBRV1-AT1 was capable of vertical transmission via asexual spores and horizontal transmission from strain TJ-NH-51S-4 to strain XJ-BZ-5-1hyg (another strain of A. tenuissima) through hyphal contact in pairing cultures. A total of 613 differentially expressed genes (DEGs) were identified in a comparative transcriptome analysis between TJ-NH-51S-4 and TJ-NH-51S-4-VF. Relative to strain TJ-NH-51S-4-VF, the number of up-regulated and down-regulated DEGs in strain TJ-NH-51S-4 was 286 and 327, respectively. Notably, the expression level of one DEG-encoding cytochrome P450 sterol 14α-demethylase and four DEGs encoding siderophore iron transporters were significantly up-regulated. To our knowledge, this is the first documentation of hypervirulence and reduced sensitivity to difenoconazole induced by AaBRV1-AT1 infection in A. tenuissima.

New Insights on the Integrated Management of Plant Diseases by RNA Strategies: Mycoviruses and RNA Interference

RNA-based strategies for plant disease management offer an attractive alternative to agrochemicals that negatively impact human and ecosystem health and lead to pathogen resistance. There has been recent interest in using mycoviruses for fungal disease control after it was discovered that some cause hypovirulence in fungal pathogens, which refers to a decline in the ability of a pathogen to cause disease. Cryphonectria parasitica, the causal agent of chestnut blight, has set an ideal model of management through the release of hypovirulent strains. However, mycovirus-based management of plant diseases is still restricted by limited approaches to search for viruses causing hypovirulence and the lack of protocols allowing effective and systemic virus infection in pathogens. RNA interference (RNAi), the eukaryotic cell system that recognizes RNA sequences and specifically degrades them, represents a promising. RNA-based disease management method. The natural occurrence of cross-kingdom RNAi provides a basis for host-induced gene silencing, while the ability of most pathogens to uptake exogenous small RNAs enables the use of spray-induced gene silencing techniques. This review describes the mechanisms behind and the potential of two RNA-based strategies, mycoviruses and RNAi, for plant disease management. Successful applications are discussed, as well as the research gaps and limitations that remain to be addressed.

ICTV Virus Taxonomy Profile: Polymycoviridae 2022

Members of the family Polymycoviridae are small viruses with multi-segmented and non-conventionally encapsidated double-stranded (ds) RNA genomes. Typically, polymycoviruses have four genomic segments, although some have up to eight. The genus Polymycovirus includes several species whose members infect fungi (ascomycetes and basidiomycetes), and oomycetes, altering host morphology, sporulation, growth and virulence. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the family Polymycoviridae, which is available at ictv.global/report/polymycoviridae.

Metagenomics Analysis of the Wheat Virome Identifies Novel Plant and Fungal-Associated Viral Sequences

Wheat viruses including wheat streak mosaic virus, Triticum mosaic virus, and barley yellow dwarf virus cost substantial losses in crop yields every year. Although there have been extensive studies conducted on these known wheat viruses, currently, there is limited knowledge about all components of the wheat (Triticum aestivum L.) virome. Here, we determined the composition of the wheat virome through total RNA deep sequencing of field-collected leaf samples. Sequences were de novo assembled after removing the host reads, and BLASTx searches were conducted. In addition to the documented wheat viruses, novel plant and fungal-associated viral sequences were identified. We obtained the full genome sequence of the first umbra-like associated RNA virus tentatively named wheat umbra-like virus in cereals. Moreover, a novel bi-segmented putative virus tentatively named wheat-associated vipovirus sharing low but significant similarity with both plant and fungal-associated viruses was identified. Additionally, a new putative fungal-associated tobamo-like virus and novel putative Mitovirus were discovered in wheat samples. The discovery and characterization of novel viral sequences associated with wheat is important to determine if these putative viruses may pose a threat to the wheat industry or have the potential to be used as new biological control agents for wheat pathogens either as wild-type or recombinant viruses.

Identification of microRNA-like RNAs from Trichoderma asperellum DQ-1 during its interaction with tomato roots using bioinformatic analysis and high-throughput sequencing

MicroRNA-like small RNAs (milRNAs) and their regulatory roles in the interaction between plant and fungus have recently aroused keen interest of plant pathologists. Trichoderma spp., one of the widespread biocontrol fungi, can promote plant growth and induce plant disease resistance. To investigate milRNAs potentially involved in the interaction between Trichoderma and tomato roots, a small RNA (sRNA) library expressed during the interaction of T. asperellum DQ-1 and tomato roots was constructed and sequenced using the Illumina HiSeqTM 2500 sequencing platform. From 13,464,142 sRNA reads, we identified 21 milRNA candidates that were similar to other known microRNAs in the miRBase database and 22 novel milRNA candidates that possessed a stable microRNA precursor hairpin structure. Among them, three milRNA candidates showed different expression level in the interaction according to the result of stem-loop RT-PCR indicating that these milRNAs may play a distinct regulatory role in the interaction between Trichoderma and tomato roots. The potential transboundary milRNAs from T. asperellum and their target genes in tomato were predicted by bioinformatics analysis. The results revealed that several interesting proteins involved in plant growth and development, disease resistance, seed maturation, and osmotic stress signal transduction might be regulated by the transboundary milRNAs. To our knowledge, this is the first report of milRNAs taking part in the process of interaction of T. asperellum and tomato roots and associated with plant promotion and disease resistance. The results might be useful to unravel the mechanism of interaction between Trichoderma and tomato.

Freeing Aspergillus fumigatus of Polymycovirus Infection Renders It More Resistant to Competition with Pseudomonas aeruginosa Due to Altered Iron-Acquiring Tactics

A virus-free (VF) A. fumigatus isolate has been shown to be resistant in competition with Pseudomonas as compared to the isogenic line infected with Aspergillus fumigatus polymycovirus 1 (AfuPmV-1), and this phenotype was apparently related to alterations in iron metabolism. Here we investigated further the mechanisms underpinning this phenotype. The extracellular siderophore profiles of five isogenic VF and virus-infected (VI) strains were sampled at 24, 31, 48, 54, and 72 h in submerged cultures and quantitatively examined by liquid chromatography and mass spectrometry. Intracellular profiles of conidia and cultures at the stationary growth phase were defined. VF A. fumigatus demonstrated the best fitness represented by the fastest onset of its exponential growth when grown on an iron-limited mineral medium. The exponential phase and transitional production phase of the extracellular triacetylfusarinine C (TafC) were achieved at 24 and 31 h, respectively, contrary to VI strains, which acted more slowly. As a result, the TafC reservoir was consumed sooner in the VF strain. Additionally, the VF strain had lower ferricrocin and higher hydroxyferricrocin content in the pellet during the stationary phase. All of these differences were significant (Kruskal–Wallis, p < 0.01). In our study, the siderophore reservoir of a VF strain was consumed sooner, improving the fitness of the VF strain in competition with P. aeruginosa.

Mycoviruses and their role in fungal pathogenesis

Nowadays, the focus of mycovirology research has expanded from plant pathogenic fungi and mycovirus mediated hypovirulence to include insect and human pathogenic fungi together with a range of mycovirus mediated phenotypes, such as hypervirulence, control of endophytic traits, regulation of metabolite production and drug resistance. In fungus-mycovirus-environmental interactions, the environment and both abiotic and biotic factors play crucial roles in whether and how mycovirus mediated phenotypes are manifest. Mycovirus infections result in alterations in the host transcriptome profile, via protein-protein interactions and triggering of antiviral RNA silencing in the fungus. These alterations, in combination with the environmental factors, may result in desirable phenotypic traits for the host, for us and in some cases for both.

Multiple Viral Infections Detected in Phytophthora condilina by Total and Small RNA Sequencing

Marine oomycetes have recently been shown to be concurrently infected by (−)ssRNA viruses of the order Bunyavirales. In this work, even higher virus variability was found in a single isolate of Phytophthora condilina, a recently described member of Phytophthora phylogenetic Clade 6a, which was isolated from brackish estuarine waters in southern Portugal. Using total and small RNA-seq the full RdRp of 13 different potential novel bunya-like viruses and two complete toti-like viruses were detected. All these viruses were successfully confirmed by reverse transcription polymerase chain reaction (RT-PCR) using total RNA as template, but complementarily one of the toti-like and five of the bunya-like viruses were confirmed when dsRNA was purified for RT-PCR. In our study, total RNA-seq was by far more efficient for de novo assembling of the virus sequencing but small RNA-seq showed higher read numbers for most viruses. Two main populations of small RNAs (21 nts and 25 nts-long) were identified, which were in accordance with other Phytophthora species. To the best of our knowledge, this is the first study using small RNA sequencing to identify viruses in Phytophthora spp.

Phenotypic and Molecular Biological Analysis of Polymycovirus AfuPmV-1M From Aspergillus fumigatus: Reduced Fungal Virulence in a Mouse Infection Model

The filamentous fungal pathogen Aspergillus fumigatus is one of the most common causal agents of invasive fungal infection in humans; the infection is associated with an alarmingly high mortality rate. In this study, we investigated whether a mycovirus, named AfuPmV-1M, can reduce the virulence of A. fumigatus in a mouse infection model. AfuPmV-1M has high sequence similarity to AfuPmV-1, one of the polymycovirus that is a capsidless four-segment double-stranded RNA (dsRNA) virus, previously isolated from the genome reference strain of A. fumigatus, Af293. However, we found the isolate had an additional fifth dsRNA segment, referred to as open reading frame 5 (ORF5), which has not been reported in AfuPmV-1. We then established isogenic lines of virus-infected and virus-free A. fumigatus strains. Mycovirus infection had apparent influences on fungal phenotypes, with the virus-infected strain producing a reduced mycelial mass and reduced conidial number in comparison with these features of the virus-free strain. Also, resting conidia of the infected strain showed reduced adherence to pulmonary epithelial cells and reduced tolerance to macrophage phagocytosis. In an immunosuppressed mouse infection model, the virus-infected strain showed reduced mortality in comparison with mortality due to the virus-free strain. RNA sequencing and high-performance liquid chromatography (HPLC) analysis showed that the virus suppressed the expression of genes for gliotoxin synthesis and its production at the mycelial stage. Conversely, the virus enhanced gene expression and biosynthesis of fumagillin. Viral RNA expression was enhanced during conidial maturation, conidial germination, and the mycelial stage. We presume that the RNA or translation products of the virus affected fungal phenotypes, including spore formation and toxin synthesis. To identify the mycovirus genes responsible for attenuation of fungal virulence, each viral ORF was ectopically expressed in the virus-free KU strain. We found that the expression of ORF2 and ORF5 reduced fungal virulence in the mouse model. In addition, ORF3 affected the stress tolerance of host A. fumigatus in culture. We hypothesize that the respective viral genes work cooperatively to suppress the pathogenicity of the fungal host.

Different Genetic Sources Contribute to the Small RNA Population in the Arbuscular Mycorrhizal Fungus Gigaspora margarita

RNA interference (RNAi) is a key regulatory pathway of gene expression in almost all eukaryotes. This mechanism relies on short non-coding RNA molecules (sRNAs) to recognize in a sequence-specific manner DNA or RNA targets leading to transcriptional or post-transcriptional gene silencing. To date, the fundamental role of sRNAs in the regulation of development, stress responses, defense against viruses and mobile elements, and cross-kingdom interactions has been extensively studied in a number of biological systems. However, the knowledge of the “RNAi world” in arbuscular mycorrhizal fungi (AMF) is still limited. AMF are obligate mutualistic endosymbionts of plants, able to provide several benefits to their partners, from improved mineral nutrition to stress tolerance. Here we described the RNAi-related genes of the AMF Gigaspora margarita and characterized, through sRNA sequencing, its complex small RNAome, considering the possible genetic sources and targets of the sRNAs. G. margarita indeed is a mosaic of different genomes since it hosts endobacteria, RNA viruses, and non-integrated DNA fragments corresponding to mitovirus sequences. Our findings show that G. margarita is equipped with a complete set of RNAi-related genes characterized by the expansion of the Argonaute-like (AGO-like) gene family that seems a common trait of AMF. With regards to sRNAs, we detected populations of sRNA reads mapping to nuclear, mitochondrial, and viral genomes that share similar features (25-nt long and 5′-end uracil read enrichments), and that clearly differ from sRNAs of endobacterial origin. Furthermore, the annotation of nuclear loci producing sRNAs suggests the occurrence of different sRNA-generating processes. In silico analyses indicate that the most abundant G. margarita sRNAs, including those of viral origin, could target transcripts in the host plant, through a hypothetical cross-kingdom RNAi.

Analysis of an Intrinsic Mycovirus Associated With Reduced Virulence of the Human Pathogenic Fungus Aspergillus fumigatus

Aspergillus fumigatus is an airborne fungal pathogen that causes severe infections with invasive growth in immunocompromised patients. Several mycoviruses have recently been isolated from A. fumigatus strains, but there are presently no reports of mycoviral-mediated reduction or elimination of fungal pathogenicity in vertebrate models. Here, we report the biological features of a novel mycovirus, A. fumigatus chrysovirus 41362 (AfuCV41362), isolated from the hypovirulent A. fumigatus strain IFM 41362. The AfuCV41362 genome is comprised of four dsRNAs, each of which contains a single ORF (ORF1-4). ORF1 encodes a protein with sequence similarity to RNA-dependent RNA polymerases of viruses in the family Chrysoviridae, while ORF3 encodes a putative capsid protein. Viral RNAs are expressed primarily during the germination stage, and RNA-seq analysis of virus-infected A. fumigatus at the germination stage suggested that the virus suppressed expression of several pathogenicity-associated host genes, including hypoxia adaptation and nitric oxide detoxification genes. In vitro functional analysis revealed that the virus-infected strain had reduced tolerance to environmental stressors. Virus-infected A. fumigatus strain IFM 41362 had reduced virulence in vivo compared to the virus-free strain in a mouse infection model. Furthermore, introduction of the mycovirus to a natively virus-free KU A. fumigatus strain induced virus-infected phenotypes. To identify mycovirus genes responsible for the reduced virulence of A. fumigatus, each viral ORF was ectopically expressed in the virus-free KU strain. Ectopic expression of the individual ORFs only nominally reduced virulence of the host fungus in a mouse infection model. However, we found that ORF3 and ORF4 reduced tolerance to environmental stresses in in vitro analysis. Based on these results, we suggest that the AfuCV41362 mycovirus ORF3 and ORF4 reduce fungal virulence by suppressing stress tolerance together with other viral genes, rather than alone.

Mycoviruses in Biological Control: From Basic Research to Field Implementation

Mycoviruses from plant pathogens can induce hypovirulence (reduced virulence) in their host fungi and have gained considerable attention as potential biocontrol tools. An increasing number of mycoviruses that induce fungal hypovirulence, from a wide variety of taxonomic groups, are currently being reported. Successful application of these viruses in disease management is greatly dependent on their ability to spread in the natural populations of the pathogen. Mycoviruses generally lack extracellular routes of transmission. Hyphal anastomosis is the main route of horizontal mycovirus transmission to other isolates, and conidia of vertical transmission to the progeny. Transmission efficiencies are influenced by both the fungal host and the infecting virus. Interestingly, artificial transfection methods have shown that potential biocontrol mycoviruses often have the ability to infect a variety of fungi. This expands their possible use to the control of pathogens others than those where they were identified. Mycovirus research is also focused on gaining insights into their complex molecular biology and the molecular bases of fungus-virus interactions. This knowledge could be exploited to manipulate the mycovirus and/or the host and generate combinations with enhanced properties in biological control. Finally, when exploring the use of mycoviruses in field conditions, the pathogen life style and the characteristics of the disease and crops affected will deeply impact the specific challenges to overcome, and the development of biocontrol formulations and delivery methods.

Integrated microRNA and mRNA analysis in the pathogenic filamentous fungus Trichophyton rubrum

BACKGROUND

Trichophyton rubrum (T. rubrum) is an important model organism of dermatophytes, which are the most common fungal pathogens worldwide. Despite the severity and prevalence of the infection caused by these pathogens, current therapies are not sufficient. MicroRNA (miRNA) is a class of small noncoding RNAs that are key factors in the regulation of gene expression. These miRNAs are reported to be highly conserved in different organisms and are involved in various essential cellular processes. In this study, we performed an integrated analysis of microRNA-like RNAs (milRNAs) and mRNAs between conidial and mycelial stages to investigate the roles of milRNAs in regulating the expression of target genes in T. rubrum.

RESULTS

A total of 158 conserved milRNAs and 12 novel milRNAs were identified in our study, corresponding to 5470 target genes, which were involved in various essential biological pathways. In addition, 137 target genes corresponding to 21 milRNAs were concurrent differentially expressed between the conidial and mycelial stages. Among these 137 target genes, 64 genes showed the opposite trend to their corresponding milRNAs in expression difference between the two stages, indicating possible negative regulation. Furthermore, 46% of differentially expressed target genes are involved in transcription, transcriptional and post-transcriptional regulation. Our results indicate that milRNAs might associate with other regulatory elements to control gene expression at both transcriptional and post-transcriptional level.

CONCLUSIONS

This study provides the first analysis of milRNA expression profile in T. rubrum as well as dermatophytes in general. The results revealed the roles of milRNAs in regulating gene expression between the two major growth stages of this fungus. Our study deepens our understanding of T. rubrum and will serve as a foundation for further investigations to combat this fungus.

De novo transcriptomic assembly and mRNA expression patterns of Botryosphaeria dothidea infection with mycoviruses chrysovirus 1 (BdCV1) and partitivirus 1 (BdPV1)

BACKGROUND

Pear ring rot, caused by Botryosphaeria species, is responsible for substantial economic losses by causing severe recession of pear tree growth in China. Mycovirus-mediated hypovirulence in plant pathogenic fungi is a crucial biological method to control fungal diseases.

METHODS

We conducted a large-scale and comprehensive transcriptome analysis to identify mRNA in B. dothidea in response to mycovirus. De novo sequencing technology from four constructed libraries of LW-C (Botryosphaeria dothidea chrysovirus 1, BdCV1), LW-P (Botryosphaeria dothidea partitivirus 1, BdPV1), LW-CP (LW-1 strain infection with BdCV1 and BdPV1), and Mock (free virus) was used to investigate and compare gene expression changes in B.dothidea strains infected with mycovirus.

RESULTS

In total, 30,058 Unigenes with an average length of 2128 bp were obtained from 4 libraries of B. dothidea strains. These were annotated to specify their classified function. We demonstrate that mRNAs of B. dothidea strains in response to mycovirus are differentially expressed. In total, 5598 genes were up-regulated and 3298 were down-regulated in the LW-CP group, 4468 were up-regulated and 4291 down-regulated in the LW-C group, and 2590 were up-regulated and 2325 down-regulated in the LW-P group. RT-qPCR was used to validate the expression of 9 selected genes. The B. dothidea transcriptome was more affected by BdCV1 infection than BdPV1. We conducted GO enrichment analysis to characterize gene functions regulated by B. dothidea with mycovirus infection. These involved metabolic process, cellular process, catalytic activity, transporter activity, signaling, and other biological pathways. KEGG function analysis demonstrated that the enriched differentially expressed genes are involved in metabolism, transcription, signal transduction, and ABC transport. mRNA is therefore involved in the interaction between fungi and mycovirus. In addition, changes in differential accumulation levels of cp and RdRp of BdCV1 and BdPV1 in B. dothidea strains were evaluated, revealing that the accumulation of BdCV1 and BdPV1 is related to the phenotype and virulence of B. dothidea strain LW-1.

CONCLUSIONS

The identification and analysis of mRNAs from B. dothidea was first reported at the transcriptome level. Our analysis provides further insight into the interaction of B. dothidea strains infection with chrysovirus 1 (BdCV1) and partitivirus 1 (BdPV1) at the transcriptome level.

Novel Partitivirus Enhances Virulence of and Causes Aberrant Gene Expression in Talaromyces marneffei

Talaromyces marneffei is the most important thermal dimorphic fungus causing systemic mycosis in Southeast Asia. We report the discovery of a novel partitivirus, Talaromyces marneffeipartitivirus-1 (TmPV1). TmPV1 was detected in 7 (12.7%) of 55 clinical T. marneffei isolates. Complete genome sequencing of the seven TmPV1 isolates revealed two double-stranded RNA (dsRNA) segments encoding RNA-dependent RNA polymerase (RdRp) and capsid protein, respectively. Phylogenetic analysis showed that TmPV1 occupied a distinct clade among the members of the genus Gammapartitivirus Transmission electron microscopy confirmed the presence of isometric, nonenveloped viral particles of 30 to 45 nm in diameter, compatible with partitiviruses, in TmPV1-infected T. marneffei Quantitative reverse transcription-PCR (qRT-PCR) demonstrated higher viral load of TmPV1 in the yeast phase than in the mycelial phase of T. marneffei Two virus-free isolates, PM1 and PM41, were successfully infected by purified TmPV1 using protoplast transfection. Mice challenged with TmPV1-infected T. marneffei isolates showed significantly shortened survival time (P < 0.0001) and higher fungal burden in organs than mice challenged with isogenic TmPV1-free isolates. Transcriptomic analysis showed that TmPV1 causes aberrant expression of various genes in T. marneffei, with upregulation of potential virulence factors and suppression of RNA interference (RNAi)-related genes. This is the first report of a mycovirus in a thermally dimorphic fungus. Further studies are required to ascertain the mechanism whereby TmPV1 enhances the virulence of T. marneffei in mice and the potential role of RNAi-related genes in antiviral defense in T. marneffeiIMPORTANCETalaromyces marneffei (formerly Penicillium marneffei) is the most important thermal dimorphic fungus in Southeast Asia, causing highly fatal systemic penicilliosis in HIV-infected and immunocompromised patients. We discovered a novel mycovirus, TmPV1, in seven clinical isolates of T. marneffei TmPV1 belongs to the genus Gammapartitivirus of the family Partitiviridae We showed that TmPV1 enhanced the virulence of T. marneffei in mice, with shortened survival time and higher fungal burden in the organs of mice challenged with TmPV1-infected T. marneffei isolates than in those of mice challenged with virus-free isogenic isolates. Transcriptomics analysis showed that TmPV1 altered the expression of genes involved in various cellular processes in T. marneffei, with upregulation of potential virulence factors and suppression of RNAi machinery which may be involved in antiviral defense. This is the first report of a mycovirus in a thermal dimorphic fungus. The present results offer insights into mycovirus-fungus interactions and pathogenesis of thermal dimorphic fungi.

Small RNA populations revealed by blocking rRNA fragments in Drosophila melanogaster reproductive tissues

RNA interference (RNAi) is a complex and highly conserved regulatory mechanism mediated via small RNAs (sRNAs). Recent technical advances in high throughput sequencing have enabled an increasingly detailed analysis of sRNA abundances and profiles in specific body parts and tissues. This enables investigations of the localized roles of microRNAs (miRNAs) and small interfering RNAs (siRNAs). However, variation in the proportions of non-coding RNAs in the samples being compared can hinder these analyses. Specific tissues may vary significantly in the proportions of fragments of longer non-coding RNAs (such as ribosomal RNA or transfer RNA) present, potentially reflecting tissue-specific differences in biological functions. For example, in Drosophila, some tissues contain a highly abundant 30nt rRNA fragment (the 2S rRNA) as well as abundant 5’ and 3’ terminal rRNA fragments. These can pose difficulties for the construction of sRNA libraries as they can swamp the sequencing space and obscure sRNA abundances. Here we addressed this problem and present a modified “rRNA blocking” protocol for the construction of high-definition (HD) adapter sRNA libraries, in D. melanogaster reproductive tissues. The results showed that 2S rRNAs targeted by blocking oligos were reduced from >80% to < 0.01% total reads. In addition, the use of multiple rRNA blocking oligos to bind the most abundant rRNA fragments allowed us to reveal the underlying sRNA populations at increased resolution. Side-by-side comparisons of sequencing libraries of blocked and non-blocked samples revealed that rRNA blocking did not change the miRNA populations present, but instead enhanced their abundances. We suggest that this rRNA blocking procedure offers the potential to improve the in-depth analysis of differentially expressed sRNAs within and across different tissues.

Mycoviruses in Aspergilli: A Comprehensive Review

Fungi, similar to all species, are susceptible to viral infection. Aspergillus is arguably the most well studied fungal genus because of its medical, ecological and economical significance. Mycoviruses were initially detected in Aspergillus species almost 50 years ago and the field continues to be active today with ground-breaking discoveries. The aim of the present review is to cover the scientific progress in all aspects of mycovirology as exemplified by Aspergillus-focused research. Initially an overview of the population studies illustrating the presence of mycoviruses in numerous important Aspergillus species, such as A. niger, A. flavus, and A. fumigatus with be presented. Moreover the intricacies of mycovirus transmission, both inter- and intra-species, will be discussed together with the methodologies used to investigate viral dispersion in a laboratory setting. Subsequently, the genomic features of all molecularly characterized mycoviruses to date will be analyzed in depth. These include members of established viral families, such as Partitiviridae, Chrysoviridae and Totiviridae, but also more recent, novel discoveries that led to the proposal of new viral families, such as Polymycoviridae, Alternaviridae and, in the context of the present review, Exartaviridae. Finally, the major issue of phenotypic effects of mycoviral infection on the host is addressed, including aflatoxin production in A. flavus, together with growth and virulence in A. fumigatus. Although the molecular mechanisms behind these phenomena are yet to be elucidated, recent studies suggest that by implication, RNA silencing may be involved.