Identification of microRNA-like RNAs in mycelial and yeast phases of the thermal dimorphic fungus Penicillium marneffei

Global discovery, expression pattern, and regulatory role of miRNA-like RNAs in Ascosphaera apis infecting the Asian honeybee larvae

Ascosphaera apis, a specialized fungal pathogen, causes lethal infection in honeybee larvae. miRNA-like small RNAs (milRNAs) are fungal small non-coding RNAs similar to miRNAs, which have been shown to regulate fungal hyphal growth, spore formation, and pathogenesis. Based on the transcriptome data, differentially expressed miRNA-like RNAs (DEmilRNAs) in A. apis infecting the Apis cerana cerana worker 4-, 5-, and 6-day-old larvae (Aa-4, Aa-5, and Aa-6) were screened and subjected to trend analysis, followed by target prediction and annotation as well as investigation of regulatory networks, with a focus on sub-networks relative to MAPK signaling pathway, glycerolipid metabolism, superoxide dismutase, and enzymes related to chitin synthesis and degradation. A total of 606 milRNAs, with a length distribution ranging from 18 nt to 25 nt, were identified. The first nucleotide of these milRNAs presented a bias toward U, and the bias patterns across bases of milRNAs were similar in the aforementioned three groups. There were 253 milRNAs, of which 68 up-and 54 down-regulated milRNAs shared by these groups. Additionally, the expression and sequences of three milRNAs were validated by stem-loop RT-PCR and Sanger sequencing. Trend analysis indicated that 79 DEmilRNAs were classified into three significant profiles (Profile4, Profile6, and Profile7). Target mRNAs of DEmilRNAs in these three significant profiles were engaged in 42 GO terms such as localization, antioxidant activity, and nucleoid. These targets were also involved in 120 KEGG pathways including lysine biosynthesis, pyruvate metabolism, and biosynthesis of antibiotics. Further investigation suggested that DEmilRNA-targeted mRNAs were associated with the MAPK signaling pathway, glycerolipid metabolism, superoxide dismutase, and enzymes related to chitin synthesis and degradation. Moreover, the binding relationships between aap-milR10516-x and ChsD as well as between aap-milR-2478-y and mkh1 were confirmed utilizing a combination of dual-luciferase reporter gene assay and RT-qPCR. Our data not only provide new insights into the A. apis proliferation and invasion, but also lay a basis for illustrating the DEmilRNA-modulated mechanisms underlying the A. apis infection.

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.

Diagnosis of Human Endemic Mycoses Caused by Thermally Dimorphic Fungi: From Classical to Molecular Methods

Human endemic mycoses are potentially fatal diseases caused by a diverse group of fungi that can alter their morphology in response to an increase in temperature. These thermally dimorphic fungi affect both healthy and immunocompromised hosts, causing a substantial health and economic burden. Despite this, the diagnosis of endemic mycoses is still a formidable challenge for several reasons, including similar symptomatology, limited utility of classical diagnostic methods, inaccessibility to reliable molecular approaches in most endemic areas, and a lack of clinical suspicion out of these regions. This review summarizes essential knowledge on thermally dimorphic fungi and the life-threatening diseases they cause. The principle, advantages and limitations of the methods traditionally used for their diagnosis are also described, along with the application status and future directions for the development of alternative diagnostic strategies, which could help to reduce the disease burden in endemic areas.

Unveiling microRNA-like small RNAs implicated in the initial infection of Fusarium oxysporum f. sp. cubense through small RNA sequencing

Banana Fusarium wilt (BFW), caused by Fusarium oxysporum f. sp. cubense (Foc), poses a major challenge to the worldwide banana industry. Fungal microRNA-like small RNAs (milRNAs) play crucial roles in regulating fungal growth, conidiation, development, and pathogenicity. However, the milRNAs and their functions in the pathogenesis of Foc remain poorly understood. In this study, we employed high-throughput sequencing and bioinformatics to profile Foc sRNAs during both pure culture and early infection stages. Our analysis identified six milRNAs exhibiting significantly upregulated expression at the initial Foc infection. Of these, milR106's biogenesis was found to be Dicer-dependent, whereas milR87, milR133, milR138, and milR148 were associated with Dicer and Argonaute proteins. Genetic manipulation and phenotype analysis confirmed that milR106 is crucial for Foc virulence by regulating conidiation, hydrogen peroxide sensitivity, and infective growth. Gene Ontology analysis of milRNA targets in the banana genome revealed enrichment in defence response to fungus and cellular response to hypoxia, implying the importance of these target genes in response to pathogen infection. In conclusion, our sRNA profiling of Foc identified several infection-induced milRNAs. The corresponding results provide valuable molecular targets for the development of an efficient strategy to control BFW.

MicroRNAs in Tetrahymena thermophila: An epigenetic regulatory mechanism in the response to cadmium stress

Among the epigenetic mechanisms based on non-coding RNA are microRNAs (miRNAs) that are involved in the post-transcriptional regulation of mRNAs. In many organisms, the expression of genes involved in the cellular response to biotic or abiotic stress depends on the regulation, generally inhibitory, performed by miRNAs. For the first time in the eukaryotic microorganism (ciliate-model) Tetrahymena thermophila, miRNAs involved in the post-transcriptional regulation of transcripts linked to the response to cadmium have been isolated and analyzed. Forty de novo miRNAs (we named tte-miRNAs) have been isolated from control and Cd-treated populations (1 or 24 h exposures). An exhaustive comparative analysis of the features of these mature tte-miRNAs and their precursor sequences (pre-tte-miRNAs) confirms that they are true miRNAs. In addition to the three types of miRNA isoforms previously described in other organisms, two new types are also described among the tte-miRNAs studied. A certain percentage of the pre-tte-miRNA sequences are in introns from genes with many introns, and have been defined as 5', 3'-tailed mirtrons. A qRT-PCR analysis of selected tte-miRNAs together with some of their targets has validated them. Cd is one of the most toxic metals for the cell, which must defend itself against its toxicity by various mechanisms, such as expulsion by membrane pumps, chelation by metallothioneins, among others. Like other toxic metals, Cd also causes a well-known series of cellular effects such as intense proteotoxicity. Many of the targets that are regulated by the tte-miRNAs are transcripts encoding proteins that fit into these defense mechanisms and toxic metal effects.

Genome-wide identification of Sclerotinia sclerotiorum small RNAs and their endogenous targets

BACKGROUND

Several phytopathogens produce small non-coding RNAs of approximately 18-30 nucleotides (nt) which post-transcriptionally regulate gene expression. Commonly called small RNAs (sRNAs), these small molecules were also reported to be present in the necrotrophic pathogen Sclerotinia sclerotiorum. S. sclerotiorum causes diseases in more than 400 plant species, including the important oilseed crop Brassica napus. sRNAs can further be classified as microRNAs (miRNAs) and short interfering RNAs (siRNAs). Certain miRNAs can activate loci that produce further sRNAs; these secondary sRNA-producing loci are called 'phased siRNA' (PHAS) loci and have only been described in plants. To date, very few studies have characterized sRNAs and their endogenous targets in S. sclerotiorum.

RESULTS

We used Illumina sequencing to characterize sRNAs from fungal mycelial mats of S. sclerotiorum spread over B. napus leaves. In total, eight sRNA libraries were prepared from in vitro, 12 h post-inoculation (HPI), and 24 HPI mycelial mat samples. Cluster analysis identified 354 abundant sRNA clusters with reads of more than 100 Reads Per Million (RPM). Differential expression analysis revealed upregulation of 34 and 57 loci at 12 and 24 HPI, respectively, in comparison to in vitro samples. Among these, 25 loci were commonly upregulated. Altogether, 343 endogenous targets were identified from the major RNAs of 25 loci. Almost 88% of these targets were annotated as repeat element genes, while the remaining targets were non-repeat element genes. Fungal degradome reads confirmed cleavage of two transposable elements by one upregulated sRNA. Altogether, 24 milRNA loci were predicted with both mature and milRNA* (star) sequences; these are both criteria associated previously with experimentally verified miRNAs. Degradome sequencing data confirmed the cleavage of 14 targets. These targets were related to repeat element genes, phosphate acetyltransferases, RNA-binding factor, and exchange factor. A PHAS gene prediction tool identified 26 possible phased interfering loci with 147 phasiRNAs from the S. sclerotiorum genome, suggesting this pathogen might produce sRNAs that function similarly to miRNAs in higher eukaryotes.

CONCLUSIONS

Our results provide new insights into sRNA populations and add a new resource for the study of sRNAs in S. sclerotiorum.

Intracellular metabolomics and microRNAomics unveil new insight into the regulatory network for potential biocontrol mechanism of stress-tolerant Tricho-fusants interacting with phytopathogen Sclerotium rolfsii Sacc

The present study employed microRNA (miRNA) sequencing and metabolome profiling of Trichoderma parental strains and fusants during normal growth and interaction with the phytopathogen Sclerotium rolfsii Sacc. In-vitro antagonism indicated that abiotic stress-tolerant Tricho-fusant FU21 was examined as a potent biocontroller with mycoparasitic action after 10 days. During interaction with the test pathogen, the most abundant uprising intracellular metabolite was recognized as l-proline, which corresponds to held-down l-alanine, associated with arginine and proline metabolism, biosynthesis of secondary metabolites, and nitrogen metabolism linked to predicted genes controlled by miRNAs viz., cel-miR-8210-3p, hsa-miR-3613-5p, and mml-miR-7174-3p. The miRNAs- mml-miR-320c and mmu-miR-6980-5p were found to be associated with phenylpropanoid biosynthesis, transcription factors, and signal transduction pathways, respectively, and were ascertained downregulated in potent FU21_IB compared with FU21_CB. The amino benzoate degradation and T cell receptor signaling pathways were regulated by miRNAs cel-miR-8210 and tca-miR-3824 as stress tolerance mechanisms of FU21. The intracellular metabolites l-proline, maleic acid, d-fructose, Myo-inositol, arabinitol, d-xylose, mannitol, and butane were significantly elevated as potential biocontrol and stress-tolerant constituents associated with miRNA regulatory pathways in potent FU21_IB. A network analysis between regulatory miRNA-predicted genes and intracellular metabolomics acknowledged possible biocontrol pathways/mechanisms in potent FU21_IB to restrain phytopathogen.

An Overlooked and Underrated Endemic Mycosis-Talaromycosis and the Pathogenic Fungus Talaromyces marneffei

Talaromycosis is an invasive mycosis endemic in tropical and subtropical Asia and is caused by the pathogenic fungus Talaromyces marneffei. Approximately 17,300 cases of T. marneffei infection are diagnosed annually, and the reported mortality rate is extremely high (~1/3). Despite the devastating impact of talaromycosis on immunocompromised individuals, particularly HIV-positive persons, and the increase in reported occurrences in HIV-uninfected persons, diagnostic and therapeutic approaches for talaromycosis have received far too little attention worldwide. In 2021, scientists living in countries where talaromycosis is endemic raised a global demand for it to be recognized as a neglected tropical disease. Therefore, T. marneffei and the infectious disease induced by this fungus must be treated with concern. T. marneffei is a thermally dimorphic saprophytic fungus with a complicated mycological growth process that may produce various cell types in its life cycle, including conidia, hyphae, and yeast, all of which are associated with its pathogenicity. However, understanding of the pathogenic mechanism of T. marneffei has been limited until recently. To achieve a holistic view of T. marneffei and talaromycosis, the current knowledge about talaromycosis and research breakthroughs regarding T. marneffei growth biology are discussed in this review, along with the interaction of the fungus with environmental stimuli and the host immune response to fungal infection. Importantly, the future research directions required for understanding this serious infection and its causative pathogenic fungus are also emphasized to identify solutions that will alleviate the suffering of susceptible individuals worldwide.

Fungal Secondary Metabolites and Small RNAs Enhance Pathogenicity during Plant-Fungal Pathogen Interactions

Fungal plant pathogens use proteinaceous effectors as well as newly identified secondary metabolites (SMs) and small non-coding RNA (sRNA) effectors to manipulate the host plant's defense system via diverse plant cell compartments, distinct organelles, and many host genes. However, most molecular studies of plant-fungal interactions have focused on secreted effector proteins without exploring the possibly equivalent functions performed by fungal (SMs) and sRNAs, which are collectively known as "non-proteinaceous effectors". Fungal SMs have been shown to be generated throughout the plant colonization process, particularly in the early biotrophic stages of infection. The fungal repertoire of non-proteinaceous effectors has been broadened by the discovery of fungal sRNAs that specifically target plant genes involved in resistance and defense responses. Many RNAs, particularly sRNAs involved in gene silencing, have been shown to transmit bidirectionally between fungal pathogens and their hosts. However, there are no clear functional approaches to study the role of these SM and sRNA effectors. Undoubtedly, fungal SM and sRNA effectors are now a treasured land to seek. Therefore, understanding the role of fungal SM and sRNA effectors may provide insights into the infection process and identification of the interacting host genes that are targeted by these effectors. This review discusses the role of fungal SMs and sRNAs during plant-fungal interactions. It will also focus on the translocation of sRNA effectors across kingdoms, the application of cross-kingdom RNA interference in managing plant diseases and the tools that can be used to predict and study these non-proteinaceous effectors.

Small RNA Analyses of a Ceratobasidium Isolate Infected with Three Endornaviruses

Isolates of three endornavirus species were identified co-infecting an unidentified species of Ceratobasidium, itself identified as a symbiont from within the roots of a wild plant of the terrestrial orchid Pterostylis vittata in Western Australia. Isogenic lines of the fungal isolate lacking all three mycoviruses were derived from the virus-infected isolate. To observe how presence of endornaviruses influenced gene expression in the fungal host, we sequenced fungus-derived small RNA species from the virus-infected and virus-free isogenic lines and compared them. The presence of mycoviruses influenced expression of small RNAs. Of the 3272 fungus-derived small RNA species identified, the expression of 9.1% (300 of 3272) of them were up-regulated, and 0.6% (18 of 3272) were down-regulated in the presence of the viruses. Fourteen novel micro-RNA-like RNAs (Cer-milRNAs) were predicted. Gene target prediction of the differentially expressed Cer-milRNAs was quite ambiguous; however, fungal genes involved in transcriptional regulation, catalysis, molecular binding, and metabolic activities such as gene expression, DNA metabolic processes and regulation activities were differentially expressed in the presence of the mycoviruses.

Comparative profiling of canonical and non-canonical small RNAs in the rice blast fungus, Magnaporthe oryzae

RNA interference (RNAi) is divided into canonical, Dicer-dependent and non-canonical, Dicer-independent pathways according to Dicer protein dependency. However, sRNAs processed in a Dicer-independent manner have not been reported in plant pathogenic fungi, including Magnaporthe oryzae. We comparatively profiled the Dicer-dependent and -independent sRNAs of M. oryzae. Dicer-dependent sRNAs were 19–24-nt in length, had low strand-specificity, and showed a preference for uracil at the 5′-end. By contrast, Dicer-independent sRNAs presented irregular patterns in length distribution, high strand-specificity, and a preference for cytosine at the penultimate position. Dicer-dependent sRNA loci were mainly associated with LTR-transposons, while Dicer-independent sRNAs were associated with protein-coding genes and transposons. We identified MoERI-1, a non-canonical RNAi component, and profiled the sRNA and mRNA transcriptomes of ΔMoeri-1 at the mycelia and conidiation stages, as the mutant showed increased conidiation. We found that genes involved in conidiation and cell cycle were upregulated by MoERI-1 deletion. Furthermore, a comparison between sRNA and mRNA transcriptome revealed that MoERI-1-dependent sRNAs mediate the regulation of gene expression. Overall, these results showed that M. oryzae has non-canonical RNAi pathways distinct to the Dicer-dependent manner and exploits MoERI-1-dependent sRNAs to regulate the conidiation process.

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 comparative study of microRNAs in different stages of Eimeria tenella

Apicomplexan parasites have divergent biogenesis machinery for small RNA generation. Analysis has shown that parasites in Plasmodium and Cryptosporidium as well as many species in Leishmania or Trypanosoma do not have a complete machinery in small RNA biogenesis. Recently, the miRNA-generating system of Toxoplasma has been identified as plant/fungal-like and its miRNAome has been elucidated. However, the microRNA (miRNA) expression profiles and their potential regulatory functions in different stages of Eimeria tenella remain largely unknown. In this study, we characterized the RNA silencing machinery of E. tenella and investigated the miRNA population distribution at different life stages by high-throughput sequencing. We characterized the expression of miRNAs in the unsporulated oocyst, sporulated oocyst and schizogony stages, obtaining a total of 392 miRNAs. We identified 58 differentially expressed miRNAs between USO (unsporulated oocysts) and SO (sporulated oocysts) that were significantly enriched for their potential target genes in the regulation of gene expression and chromatin binding, suggesting an epigenetic modulation of sporulating by these miRNAs. In comparing miRNA expression at endogenous and exogenous developmental stages, twenty-four miRNAs were identified differently expressed. Those were mainly associated with the regulation of genes with protein kinase activity, suggesting control of protein phosphorylation. This is the first study about the evolution of miRNA biogenesis system and miRNA control of gene expression in Eimeria species. Our data may lead to functional insights into of the regulation of gene expression during parasite life cycle in apicomplexan parasites.

RNA-Seq, Bioinformatic Identification of Potential MicroRNA-like Small RNAs in the Edible Mushroom Agaricus bisporus and Experimental Approach for Their Validation

Although genomes from many edible mushrooms are sequenced, studies on fungal micro RNAs (miRNAs) are scarce. Most of the bioinformatic tools are designed for plants or animals, but the processing and expression of fungal miRNAs share similarities and differences with both kingdoms. Moreover, since mushroom species such as Agaricus bisporus (A. bisporus, white button mushroom) are frequently consumed as food, controversial discussions are still evaluating whether their miRNAs might or might not be assimilated, perhaps within extracellular vesicles (i.e., exosomes). Therefore, the A. bisporus RNA-seq was studied in order to identify potential de novo miRNA-like small RNAs (milRNAs) that might allow their later detection in diet. Results pointed to 1 already known and 37 de novo milRNAs. Three milRNAs were selected for RT-qPCR experiments. Precursors and mature milRNAs were found in the edible parts (caps and stipes), validating the predictions carried out in silico. When their potential gene targets were investigated, results pointed that most were involved in primary and secondary metabolic regulation. However, when the human transcriptome is used as the target, the results suggest that they might interfere with important biological processes related with cancer, infection and neurodegenerative diseases.

Discovery of microRNA-like Small RNAs in Pathogenic Plant Fungus Verticillium nonalfalfae Using High-Throughput Sequencing and qPCR and RLM-RACE Validation

Verticillium nonalfalfae (V. nonalfalfae) is one of the most problematic hop (Humulus lupulus L.) pathogens, as the highly virulent fungal pathotypes cause severe annual yield losses due to infections of entire hop fields. In recent years, the RNA interference (RNAi) mechanism has become one of the main areas of focus in plant-fungal pathogen interaction studies and has been implicated as one of the major contributors to fungal pathogenicity. MicroRNA-like RNAs (milRNAs) have been identified in several important plant pathogenic fungi; however, to date, no milRNA has been reported in the V. nonalfalfae species. In the present study, using a high-throughput sequencing approach and extensive bioinformatics analysis, a total of 156 milRNA precursors were identified in the annotated V. nonalfalfae genome, and 27 of these milRNA precursors were selected as true milRNA candidates, with appropriate microRNA hairpin secondary structures. The stem-loop RT-qPCR assay was used for milRNA validation; a total of nine V. nonalfalfae milRNAs were detected, and their expression was confirmed. The milRNA expression patterns, determined by the absolute quantification approach, imply that milRNAs play an important role in the pathogenicity of highly virulent V. nonalfalfae pathotypes. Computational analysis predicted milRNA targets in the V. nonalfalfae genome and in the host hop transcriptome, and the activity of milRNA-mediated RNAi target cleavage was subsequently confirmed for two selected endogenous fungal target gene models using the 5' RLM-RACE approach.

Novel Antarctic yeast adapts to cold by switching energy metabolism and increasing small RNA synthesis

The novel extremophilic yeast Rhodotorula frigidialcoholis, formerly R. JG1b, was isolated from ice-cemented permafrost in University Valley (Antarctic), one of coldest and driest environments on Earth. Phenotypic and phylogenetic analyses classified R. frigidialcoholis as a novel species. To characterize its cold-adaptive strategies, we performed mRNA and sRNA transcriptomic analyses, phenotypic profiling, and assessed ethanol production at 0 and 23 °C. Downregulation of the ETC and citrate cycle genes, overexpression of fermentation and pentose phosphate pathways genes, growth without reduction of tetrazolium dye, and our discovery of ethanol production at 0 °C indicate that R. frigidialcoholis induces a metabolic switch from respiration to ethanol fermentation as adaptation in Antarctic permafrost. This is the first report of microbial ethanol fermentation utilized as the major energy pathway in response to cold and the coldest temperature reported for natural ethanol production. R. frigidialcoholis increased its diversity and abundance of sRNAs when grown at 0 versus 23 °C. This was consistent with increase in transcription of Dicer, a key protein for sRNA processing. Our results strongly imply that post-transcriptional regulation of gene expression and mRNA silencing may be a novel evolutionary fungal adaptation in the cryosphere.

Serum Exosomal microRNA Profiling in AIDS Complicated with Talaromyces marneffei Infection

Introduction

In order to find the early diagnostic markers of AIDS combined with TM infection, we detected and analyzed the serum exosomal miRNAs of AIDS patients with or without TM infection.

Materials and Methods

We profiled the expression levels of miRNAs via RNA sequencing in serum exosomes from the pooled samples of 17 AIDS patients combined without TM infection and 15 AIDS combined with TM infection patients. For external validation, we validated these results using real-time quantification polymerase chain reaction (qRT-PCR) in an independent cohort of 35 AIDS patients combined without TM infection and 33 AIDS combined with TM infection patients. Finally, bioinformatics was used to predict the target genes and pathways of meaningful miRNAs.

Results

A total of 131 serum exosomal miRNAs including 73 up-regulated and 59 down-regulated miRNAs were found to be differentially expressed (log2FC≥1 and FDR <0.01) in the two groups. A validation analysis revealed that three miRNAs (miR-192-5p, miR-194-5p and miR-1246) were upregulated in exosomes from AIDS combined with TM infection patients. ROC analyses showed that the AUC in combined diagnosis of the three miRNAs was 0.742, and the diagnostic sensitivity and specificity were 0.568 and 0.861, respectively. In the biological process analysis, all the 3 miRNAs were involved in transcription, DNA-templated and positive regulation of transcription from RNA polymerase II promoter. At the same time, the related pathways were involved in TGF-β signaling pathway, AMPK signaling pathway, Wnt signaling pathway, MAPK signaling pathway, cGMP-PKG and cAMP signaling pathway, etc.

Conclusion

miR-192-5p, miR-194-5p and miR-1246 in serum exosomes might be a potential biomarker for AIDS combined with TM infection patients, which may be involved in TGF-β signaling pathway, AMPK signaling pathway, Wnt signaling pathway, MAPK signaling pathway, cGMP-PKG and cAMP signaling pathway, etc. Further research is needed on the biological functions and mechanisms of these miRNAs.

A Novel DCL2-Dependent Micro-Like RNA Vm-PC-3p-92107_6 Affects Pathogenicity by Regulating the Expression of Vm-VPS10 in Valsa mali

Dicer proteins are mainly responsible for generating small RNAs (sRNAs), which are involved in gene silencing in most eukaryotes. In previous research, two DCL proteins in Valsa mali, the pathogenic fungus causing apple tree Valsa canker, were found associated with both the pathogenicity and generation of sRNAs. In this study, the differential expression of small interfering RNAs (siRNAs) and miRNA-like RNAs (milRNAs) was analyzed based on the deep sequencing of the wild type and Vm-DCL2 mutant, respectively. Overall, the generation of 40 siRNAs and 18 milRNAs was evidently associated with Vm-DCL2. The target genes of milRNAs were then identified using degradome sequencing; according to the prediction results, most candidate targets are related to pathogenicity. Further, expression of Vm-PC-3p-92107_6 was confirmed in the wild type but not in the Vm-DCL2 mutant. Moreover, the pathogenicity of Vm-PC-3p-92107_6 deletion mutants (ΔVm-PC-3p-92107_6) and the over-expression transformants (Vm-PC-3p-92107_6-OE) was significantly increased and decreased, respectively. Based on those degradome results, vacuolar protein sorting 10 (Vm-VPS10) was identified as the target of Vm-PC-3p-92107_6. Co-expression analysis in tobacco leaves further confirmed that Vm-PC-3p-92107_6 could suppress the expression of Vm-VPS10. Meanwhile, the expression levels of Vm-PC-3p-92107_6 and Vm-VPS10 displayed divergent trends in ΔVm-PC-3p-92107_6 and Vm-PC-3p-92107_6-OE, respectively. Perhaps most importantly, ΔVm-VPS10 featured a significant reduction in pathogenicity. Taken together, our results indicate that a DCL2-dependent milRNA Vm-PC-3p-92107_6 plays roles in pathogenicity by regulating the expression of Vm-VPS10. This study lays a foundation for the comprehensive analysis of pathogenic mechanisms of V. mali and deepens our understanding of the generation and function of fungal sRNA.

Role of Dicer-Dependent RNA Interference in Regulating Mycoparasitic Interactions

Dicer-like proteins (DCLs) play a vital role in RNA interference (RNAi), by cleaving RNA filament into small RNAs. Although DCL-mediated RNAi can regulate interspecific communication between pathogenic/mutualistic organisms and their hosts, its role in mycoparasitic interactions is yet to be investigated. In this study, we deleted dcl genes in the mycoparasitic fungus Clonostachys rosea and characterize the functions of DCL-dependent RNAi in mycoparasitism. Deletion of dcl2 resulted in a mutant with reduced secondary metabolite production, antagonism toward the plant-pathogenic fungus Botrytis cinerea, and reduced ability to control Fusarium foot rot disease on wheat, caused by Fusarium graminearum. Transcriptome sequencing of the in vitro interaction between the C. rosea Δdcl2 strain and B. cinerea or F. graminearum identified the downregulation of genes coding for transcription factors, membrane transporters, hydrolytic enzymes, and secondary metabolites biosynthesis enzymes putatively involved in antagonistic interactions, in comparison with the C. rosea wild-type interaction. A total of 61 putative novel microRNA-like RNAs (milRNAs) were identified in C. rosea, and 11 were downregulated in the Δdcl2 mutant. In addition to putative endogenous gene targets, these milRNAs were predicted to target B. cinerea and F. graminearum virulence factor genes, which showed an increased expression during interaction with the Δdcl2 mutant incapable of producing the targeting milRNAs. In summary, this study constitutes the first step in elucidating the role of RNAi in mycoparasitic interactions, with important implications for biological control of plant diseases, and poses the base for future studies focusing on the role of cross-species RNAi regulating mycoparasitic interactions. IMPORTANCE Small RNAs mediated RNA interference (RNAi) known to regulate several biological processes. Dicer-like endoribonucleases (DCLs) play a vital role in the RNAi pathway by generating sRNAs. In this study, we investigated a role of DCL-mediated RNAi in interference interactions between mycoparasitic fungus Clonostachys rosea and the two fungal pathogens Botrytis cinerea and Fusarium graminearum (here called mycohosts). We found that the dcl mutants were not able to produce 11 sRNAs predicted to finetune the regulatory network of genes known to be involved in production of hydrolytic enzymes, antifungal compounds, and membrane transporters needed for antagonistic action of C. rosea. We also found C. rosea sRNAs putatively targeting known virulence factors in the mycohosts, indicating RNAi-mediated cross-species communication. Our study expanded the understanding of underlying mechanisms of cross-species communication during interference interactions and poses a base for future works studying the role of DCL-based cross-species RNAi in fungal interactions.

Identification and physiological function of one microRNA (Po-MilR-1) in oyster mushroom Pleurotus ostreatus

MicroRNAs are essential regulators of gene expression and have been extensively studied in plants and animals; however, few reports have been published in mushrooms. Po-MilR-1 is a novel microRNA with a length of 22 bp in Pleurotus ostreatus. The secondary structures of five precursors and the target genes of Po-MilR-1 were predicted. Expression profile analysis showed Po-MilR-1 had specific expression in the primordium and fruiting body. To explore its physiological function, Po-MilR-1 was overexpressed in P. ostreatus. The transformants showed slow mycelium growth rate and abnormal pileus with irregular edge, which suggested Po-MilR-1 plays an important role in P. ostreatus development. Additionally, Po-MilR-1 and one of its target hydrophobin genes POH1 had opposite temporal expression profiles in the primordium and fruiting body, which revealed that Po-MilR-1 may perform its physiological function through the negative regulation of POH1. This study explored the development-related function of a mushroom microRNA and will provide a reference for other microRNAs.

Vm-milR37 contributes to pathogenicity by regulating glutathione peroxidase gene VmGP in Valsa mali

MicroRNAs play important roles in various biological processes by regulating their corresponding target genes. However, the function and regulatory mechanism of fungal microRNA-like RNAs (milRNAs) are still largely unknown. In this study, a milRNA (Vm-milR37) was isolated and identified from Valsa mali, which causes the most serious disease on the trunk of apple trees in China. Based on the results of deep sequencing and quantitative reverse transcription PCR, Vm-milR37 was found to be expressed in the mycelium, while it was not expressed during the V. mali infection process. Overexpression of Vm-milR37 did not affect vegetative growth, but significantly decreased pathogenicity. Based on degradome sequencing, the target of Vm-milR37 was identified as VmGP, a glutathione peroxidase. The expression of Vm-milR37 and VmGP showed a divergent trend in V. mali-apple interaction samples and Vm-milR37 overexpression transformants. The expression of VmGP could be suppressed significantly by Vm-milR37 when coexpressed in tobacco leaves. Deletion of VmGP showed significantly reduced pathogenicity compared with the wild type. VmGP deletion mutants showed more sensitivity to hydrogen peroxide. Apple leaves inoculated with Vm-milR37 overexpression transformants and VmGP deletion mutant displayed increased accumulation of reactive oxygen species compared with the wild type. Thus, Vm-milR37 plays a critical role in pathogenicity by regulating VmGP, which contributes to the oxidative stress response during V. mali infection. These results provide important evidence to define the roles of milRNAs and their corresponding target genes in pathogenicity.

Bioinformatic analysis of the pathogenic mechanism of talaromyces marneffei infection

BACKGROUND

Talaromyces marneffei (T marneffei), known as a significant pathogen in patients with AIDS in Southeast Asia, is a dimorphic fungus, which can cause deadly systematic infection in immunocompromised hosts. What is more, the dimorphic phase transition has been reported as a conspicuous process linked with virulence. Interestingly, the yeast form was found in infected individuals, representing the pathogenic phase. However, few researches were found to study the mechanism of dimorphic transition. Thus, a diverse insight into the dimorphic switch mechanism, is urgently needed and we are the first one to research the mechanism of dimorphism.

METHODS

Firstly, we investigated the microarray of T. marneffei in the Gene Expression Omnibus database (GEO) for differentially expressed genes (DEGs). Then Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.8 was employed to analyze the underlying enrichment and pathway in biological process of DEGs. Meanwhile, protein-protein interaction (PPI) network was constructed using STRING database. On the strength of the theory that similar amino acid sequences share similar structures, which play a decisive role on the function of protein, three dimensional structures of hub-genes were predicted to further investigate the likely function of hub-genes.

RESULTS

GSE51109 was elected as the eligible series for the purpose of our research, including GSM1238923 (GSM23), GSM1238924 (GSM24), and GSM1238925 (GSM25). PMAA_012920, PMAA_028730, PMAA_068140, PMAA_092900, PMAA_032350 were the most remarkable genes in all of the three PPI networks, thus, were viewed as hub-genes. With regard to the three-dimensional construction, except that there was no significant prediction structure of PMAA_092900 with the criterion seq identify > 30%, GMQE: 0-1, QMEAN4: -4-0, the parallel templates for four structures were Crystal structure of Saccharomyces cerevesiae mitochondrial NADP(+)-dependent isocitrate dehydrogenase in complex with isocitrate, Organellar two-pore channels (TPCs), Yeast Isocitrate Dehydrogenase (Apo Form) and Crystal Structure Of ATP-Dependent Phosphoenolpyruvate Carboxykinase From Thermus thermophilus HB8 in order.

CONCLUSION

The dimorphic transition of T. marneffei was viewed as a pathogenic factor and DEGs were observed. In-depth study of the function and pathway of DEGs revealed that PMAA_012920, PMAA_028730, PMAA_068140, PMAA_092900, PMAA_032350 were most likely acting as the hub-genes and were likely taking effect through regulating energy metabolism.

Characterization of microRNA-like RNAs associated with sclerotial development in Sclerotinia sclerotiorum

Sclerotinia sclerotiorum is a model necrotrophic pathogen causing great economic losses worldwide. Sclerotia are dormant structures that play significant biological and ecological roles in the life and disease cycles of S. sclerotiorum and other species of sclerotia-forming fungi. microRNA-like RNAs (milRNAs) as non-coding small RNAs play regulatory roles in fungal development and pathogenicity. Therefore, milRNAs associated with sclerotial development in S. sclerotiorum were investigated in this study. A total of 275 milRNAs with induced expression during sclerotia development were identified, in which 51 were differentially expressed. The target genes of all milRNAs were predicted. The putative functions of the targets regulated by milRNAs were annotated by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. The expression levels of six selected milRNAs that coordinated with their corresponding targets were validated by qRT-PCR. Among these six milRNAs, Ssc-milR-240 was potentially associated with sclerotial development by epigenetic regulation of its target histone acetyltransferase. This study will facilitate the better understanding of the milRNA regulation associated with sclerotial development in S. sclerotiorum and even other sclerotia-forming fungi. This work will provide novel insights into the molecular regulations of fungal morphogenesis and the candidate targets of milRNAs used for the sustainable management of plant diseases caused by S. sclerotiorum.

Computational prediction of microRNAs in Histoplasma capsulatum

MicroRNAs (miRNAs) are small and non-coding RNAs that regulate gene expression through post-transcriptional regulation. Although, the standard miRNA repository, MiRBase, lists more than 200 organisms having miRNA mediated regulation mechanism and thousands of miRNAs, there is not enough information about miRNAs of fungal species. Considering that there are various fungal pathogens causing disease phenotypes, it is important to search for miRNAs of those organisms. The leading cause of endemic mycosis in the USA is a fungal disease known as histoplasmosis, which is resulted by infection with a fungal intracellular parasite, Histoplasma capsulatum (H. capsulatum). In this work, genomes of H. capsulatum strains NAm1 and G217B were explored for potential miRNA like sequences and structures. Through a complex workflow involving miRNA detection and target prediction, several miRNA candidates of H. capsulatum and their possible targets in human were identified. The results presented here indicate that H. capsulatum might be one of the fungal pathogens having a miRNA based post-transcriptional gene regulation mechanism and it might have a miRNA mediated host - parasite interaction with human.

Plant-pathogen interactions: MicroRNA-mediated trans-kingdom gene regulation in fungi and their host plants

MicroRNAs (miRNAs) have been prevalently studied in plants, animals, and viruses. However, recent studies show evidences of miRNA-like RNAs (milRNAs) in fungi as well. It is known that after successful infection, pathogens hijack the host machinery and use it for their own growth and multiplication. Alternatively, resistant plants can overcome the pathogen attack by a variety of mechanisms. Based on this prior knowledge, we computationally predicted milRNAs from 13 fungi, and identified their targets in transcriptomes of the respective fungi as well as their host plants. The expressions of the milRNAs and targets were confirmed using qRT-PCR. We found that plant miRNAs targeted fungal virulence genes, while fungal milRNAs targeted plant resistance genes; corroborating miRNA-mediated trans-kingdom gene regulation and the roles of miRNAs in plant-pathogen interactions. Transgenic plants with miRNAs targeting fungal virulence genes, or anti-sense of fungal milRNAs, would be expected to be highly resistant to the fungal pathogens.

Silencing Dicer-Like Genes Reduces Virulence and sRNA Generation in Penicillium italicum, the Cause of Citrus Blue Mold

The Dicer protein is one of the most important components of RNAi machinery because it regulates the production of small RNAs (sRNAs) in eukaryotes. Here, Dicer1-like gene (Pit-DCL1) and Dicer2-like gene (Pit-DCL2) RNAi transformants were generated via pSilent-1 in Penicillium italicum (Pit), which is the causal agent of citrus blue mold. Neither transformant showed a change in mycelial growth or sporulation ability, but the pathogenicity of the Pit-DCL2 RNAi transformant to citrus fruits was severely impaired, compared to that of the Pit-DCL1 RNAi transformant and the wild type. We further developed a citrus wound-mediated RNAi approach with a double-stranded fragment of Pit-DCL2 generated in vitro, which achieved an efficiency in reducing Pi-Dcl2 expression and virulence that was similar to that of protoplast-mediated RNAi in P. italicum, suggesting that this approach is promising in the exogenous application of dsRNA to control pathogens on the surface of citrus fruits. In addition, sRNA sequencing revealed a total of 69.88 million potential sRNAs and 12 novel microRNA-like small RNAs (milRNAs), four of which have been predicated on target innate immunity or biotic stress-related genes in Valencia orange. These data suggest that both the Pit-DCL1 and Pit-DCL2 RNAi transformants severely disrupted the biogenesis of the potential milRNAs, which was further confirmed for some milRNAs by qRT-PCR or Northern blot analysis. These data suggest the sRNAs in P. italicum that may be involved in a molecular virulence mechanism termed cross-kingdom RNAi (ck-RNAi) by trafficking sRNA from P. italicum to citrus fruits.

milRNApredictor: Genome-free prediction of fungi milRNAs by incorporating k-mer scheme and distance-dependent pair potential

MicroRNA-like small RNAs (milRNAs) with length of 21-22 nucleotides are a type of small non-coding RNAs that are firstly found in Neurospora crassa in 2010. Identifying milRNAs of species without genomic information is a difficult problem. Here, knowledge-based energy features are developed to identify milRNAs by tactfully incorporating k-mer scheme and distance-dependent pair potential. Compared with k-mer scheme, features developed here can alleviate the inherent curse of dimensionality in k-scheme once k becomes large. In addition, milRNApredictor built on novel features performs comparably to k-mer scheme, and achieves sensitivity of 74.21%, and specificity of 75.72% based on 10-fold cross-validation. Furthermore, for novel miRNA prediction, there exists high overlap of results from milRNApredictor and state-of-the-art mirnovo. However, milRNApredictor is simpler to use with reduced requirements of input data and dependencies. Taken together, milRNApredictor can be used to de novo identify fungi milRNAs and other very short small RNAs of non-model organisms.

Exogenous RNAi mechanisms contribute to transcriptome adaptation by phased siRNA clusters in Paramecium

Extensive research has characterized distinct exogenous RNAi pathways interfering in gene expression during vegetative growth of the unicellular model ciliate Paramecium. However, role of RNAi in endogenous transcriptome regulation, and environmental adaptation is unknown. Here, we describe the first genome-wide profiling of endogenous sRNAs in context of different transcriptomic states (serotypes). We developed a pipeline to identify, and characterize 2602 siRNA producing clusters (SRCs). Our data show no evidence that SRCs produce miRNAs, and in contrast to other species, no preference for strand specificity of siRNAs. Interestingly, most SRCs overlap coding genes and a separate group show siRNA phasing along the entire open reading frame, suggesting that the mRNA transcript serves as a source for siRNAs. Integrative analysis of siRNA abundance and gene expression levels revealed surprisingly that mRNA and siRNA show negative as well as positive associations. Two RNA-dependent RNA Polymerase mutants, RDR1 and RDR2, show a drastic loss of siRNAs especially in phased SRCs accompanied with increased mRNA levels. Importantly, most SRCs depend on both RDRs, reminiscent to primary siRNAs in the RNAi against exogenous RNA, indicating mechanistic overlaps between exogenous and endogenous RNAi contributing to flexible transcriptome adaptation.

Comparison of the RNA Content of Extracellular Vesicles Derived from Paracoccidioides brasiliensis and Paracoccidioides lutzii

Paracoccidioides brasiliensis and P. lutzii cause human paracoccidioidomycosis. We have previously characterized the <200-nt RNA sub-populations contained in fungal extracellular vesicles (EVs) from P. brasiliensis Pb18 and other pathogenic fungi. We have presently used the RNA-seq strategy to compare the <200- and >200-nt RNA fractions contained in EVs isolated from culture supernatants of P. brasiliensis Pb18, Pb3, and P. lutzii Pb01. Shared mRNA sequences were related to protein modification, translation, and DNA metabolism/biogenesis, while those related to transport and oxidation-reduction were exclusive to Pb01. The presence of functional full-length mRNAs was validated by in vitro translation. Among small non-coding (nc)RNA, 15 were common to all samples; small nucleolar (sno)RNAs were enriched in P. brasiliensis EVs, whereas for P. lutzii there were similar proportions of snoRNA, rRNA, and tRNA. Putative exonic sRNAs were highly abundant in Pb18 EVs. We also found sRNA sequences bearing incomplete microRNA structures mapping to exons. RNA-seq data suggest that extracellular fractions containing Pb18 EVs can modulate the transcriptome of murine monocyte-derived dendritic cells in a transwell system. Considering that sRNA classes are involved in transcription/translation modulation, our general results may indicate that differences in virulence among fungal isolates can be related to their distinct EV-RNA content.

RNA interference core components identified and characterised in Verticillium nonalfalfae, a vascular wilt pathogenic plant fungi of hops

The conserved RNA interference mechanism (RNAi) in the fungal kingdom has become a focus of intense scientific investigation. The three catalytic core components, Dicer-like (DCL), Argonaute (AGO), and RNA-dependent RNA polymerase (RdRP), and their associated small interfering RNA molecules (siRNAs) have been identified and characterised in several fungal species. Recent studies have proposed that RNAi is a major contributor to the virulence of fungal pathogens as a result of so-called trans-kingdom RNA silencing. In the present study, we report on the existence of three core RNAi proteins in the pathogenic plant fungus Verticillium nonalfalfae, which is a soilborne plant pathogen that causes severe wilting disease in hops (Humulus lupulus L.). Two DCL proteins, two AGO proteins, and two RdRP proteins were identified, and their conserved RNAi domains were characterised. Our phylogeny results confirm the existing taxonomic relationships in the Ascomycete fungal phylum and show that the fungi of the Hypocreomycetidae subclass of the Sordariomycetes class have high amino acid sequence similarity. The expression analysis revealed a potential role of RNAi in the pathogenicity of the fungi, since all the RNAi genes were highly upregulated in the highly virulent isolate T2 and were also differentially expressed in the V. nonalfalfae-susceptible Celeia and V. nonalfalfae-resistant Wye Target cultivars.

Comprehensive analysis of full genome sequence and Bd-milRNA/target mRNAs to discover the mechanism of hypovirulence in Botryosphaeria dothidea strains on pear infection with BdCV1 and BdPV1

Pear ring rot disease, mainly caused by Botryosphaeria dothidea, is widespread in most pear and apple-growing regions. Mycoviruses are used for biocontrol, especially in fruit tree disease. BdCV1 (Botryosphaeria dothidea chrysovirus 1) and BdPV1 (Botryosphaeria dothidea partitivirus 1) influence the biological characteristics of B. dothidea strains. BdCV1 is a potential candidate for the control of fungal disease. Therefore, it is vital to explore interactions between B. dothidea and mycovirus to clarify the pathogenic mechanisms of B. dothidea and hypovirulence of B. dothidea in pear. A high-quality full-length genome sequence of the B. dothidea LW-Hubei isolate was obtained using Single Molecule Real-Time sequencing. It has high repeat sequence with 9.3% and DNA methylation existence in the genome. The 46.34 Mb genomes contained 14,091 predicted genes, which of 13,135 were annotated. B. dothidea was predicted to express 3833 secreted proteins. In bioinformatics analysis, 351 CAZy members, 552 transporters, 128 kinases, and 1096 proteins associated with plant-host interaction (PHI) were identified. RNA-silencing components including two endoribonuclease Dicer, four argonaute (Ago) and three RNA-dependent RNA polymerase (RdRp) molecules were identified and expressed in response to mycovirus infection. Horizontal transfer of the LW-C and LW-P strains indicated that BdCV1 induced host gene silencing in LW-C to suppress BdPV1 transmission. To investigate the role of RNA-silencing in B. dothidea defense, we constructed four small RNA libraries and sequenced B. dothidea micro-like RNAs (Bd-milRNAs) produced in response to BdCV1 and BdPV1 infection. Among these, 167 conserved and 68 candidate novel Bd-milRNAs were identified, of which 161 conserved and 20 novel Bd-milRNA were differentially expressed. WEGO analysis revealed involvement of the differentially expressed Bd-milRNA-targeted genes in metabolic process, catalytic activity, cell process and response to stress or stimulus. BdCV1 had a greater effect on the phenotype, virulence, conidiomata, vertical and horizontal transmission ability, and mycelia cellular structure biological characteristics of B. dothidea strains than BdPV1 and virus-free strains. The results obtained in this study indicate that mycovirus regulates biological processes in B. dothidea through the combined interaction of antiviral defense mediated by RNA-silencing and milRNA-mediated regulation of target gene mRNA expression.

Degradation of Fungal MicroRNAs Triggered by Short Tandem Target Mimics Is via the Small-RNA-Degrading Nuclease

MicroRNAs (miRNAs) have been recognized as sequence-specific regulators of the genome, transcriptome, and proteome in eukaryotes. However, the functions and working mechanisms of hundreds of fungal miRNA-like (miR-like) RNAs are obscure. Here, we report that a short tandem target mimic (STTM) triggered the degradation of several fungal miR-like RNAs in two different fungal species, Metarhizium robertsii and Aspergillus flavus, and that small-RNA-degrading nucleases (SDNs) were indispensable for such degradation. STTMs were most effective when the fungal polymerase II (Pol II) promoter was used for their expression, while the Pol III promoter was less effective. The length of the STTM spacer, approximately 48 to 96 nucleotides, and the number of miR-like RNA binding sites, from 2 to 4 copies, showed no significant difference in the degradation of miR-like RNAs. STTMs modulated the miR-like RNA expression levels in at least two different fungal species, which further impacted fungal asexual growth and sporulation. Further analysis showed that the degraded miR-like RNAs in STTM mutants led to the upregulation of potential target genes involved in fungal development and conidial production, which result in different phenotypes in these mutants. The STTM technology developed in this study is an effective and powerful tool for the functional dissection of fungal miR-like RNAs.IMPORTANCE The development and application of STTM technology to block miR-like RNAs in M. robertsii and A. flavus may allow for efficient generation of miR-like RNA mutants in various fungi, providing a powerful tool for functional genomics of small RNA molecules in fungi.

In silico characterization of microRNAs-like sequences in the genome of Paracoccidioides brasiliensis

Eukaryotic cells have different mechanisms of post-transcriptional regulation. Among these mechanisms, microRNAs promote regulation of targets by cleavage or degradation of the mRNA. Fungi of the Paracoccidioides complex are the etiological agents of the main systemic mycosis of Latin America. These fungi present a plasticity to adapt and survive in different conditions, and the presence of microRNAs-like molecules could be part of the mechanisms that provide such plasticity. MicroRNAs produced by the host influence the progression of this mycosis in the lungs besides regulating targets involved in apoptosis in macrophage, activation of T and B cells and the production of cytokines. Therefore, this work analyzed the presence of regions in the genome of this fungus with a potential to encode microRNAs-like molecules. Here we show by analysis of sequence similarity the presence of 18 regions, putatively coding for microRNAs-like molecules in the Paracoccidioides brasiliensis genome. We also described the conservation of dicer and argonaut proteins and the cognate transcripts induced in the yeast parasitic phase. This work represents a starting point for the analysis of the presence of those molecules in the morphological stages of the fungus and their role in fungal development.

milR4 and milR16 Mediated Fruiting Body Development in the Medicinal Fungus Cordyceps militaris

Cordyceps militaris readily performs sexual reproduction, thus providing a remarkably rich model for understanding the processes involved in sexual development. It could regulate expression of human genes by diet-derived miRNA-like RNAs (milRNAs). However, the study of miRNAs in C. militaris has been limited. In the present study, genes encoding Dicers, Argonautes, and RNA-dependent RNA polymerases were identified. Illumina deep sequencing was performed to characterize the milRNAs in C. militaris at asexual and sexual development stages. Total 38 milRNAs were identified and five milRNAs were validated by northern blot and qRT-PCR, out of which, 19 were specific for sexual development. Importantly, the fungi could not form fruiting bodies after disruption of milR4, while the perithecium was formed in advance after over-expression of milR4. Abnormal pale yellow fruiting body primordium, covered with abnormal primordium, was formed in the strain with miR16 disruption. Although no milR4 or milR16 target genes were identified, differential expression of many different genes involved in mycelium growth and sexual development (mating process, mating signaling, and fruiting body development) among these mutants were found. Overall, milRNAs play vital roles in sexual development in C. militaris.

Discovery and profiling of small RNAs from Puccinia triticina by deep sequencing and identification of their potential targets in wheat

Cross-kingdom RNAi is a well-documented phenomenon where sRNAs generated by host and pathogens may govern resistance or susceptible phenotypes during host-pathogen interaction. With the first example of the direct involvement of fungal generated sRNAs in virulence of plant pathogenic fungi Botrytis cinerea and recently from Puccinia striiformis f. sp. tritici, we attempted to identify sRNAs in Puccinia triticina (P. triticina). Four sRNA libraries were prepared and sequenced using Illumina sequencing technology and a total of ~ 1-1.28 million potential sRNAs and two microRNA-like small RNA (mil-RNAs) candidates were identified. Computational prediction of targets using a common set of sRNAs and P. triticina mil-RNAs (pt-mil-RNAs) within P. triticina and wheat revealed the majority of the targets as repetitive elements in P. triticina whereas in wheat, the target genes were identified to be involved in many biological processes including defense-related pathways. We found 9 receptor-like kinases (RLKs) and 14 target genes of each related to reactive oxygen species (ROS) pathway and transcription factors respectively, including significant numbers of target genes from various other categories. Expression analysis of twenty selected sRNAs, targeting host genes pertaining to ROS related, disease resistance, metabolic processes, transporter, apoptotic inhibitor, and transcription factors along with two pt-mil-RNAs by qRT-PCR showed distinct patterns of expression of the sRNAs in urediniospore-specific libraries. In this study, for the first time, we report identification of novel sRNAs identified in P. triticina including two pt-mil-RNAs that may play an important role in biotrophic growth and pathogenicity.

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.

Transcriptional and Small RNA Responses of the White Mold Fungus Sclerotinia sclerotiorum to Infection by a Virulence-Attenuating Hypovirus

Mycoviruses belonging to the family Hypoviridae cause persistent infection of many different host fungi. We previously determined that the white mold fungus, Sclerotiniasclerotiorum, infected with Sclerotinia sclerotiorum hypovirus 2-L (SsHV2-L) exhibits reduced virulence, delayed/reduced sclerotial formation, and enhanced production of aerial mycelia. To gain better insight into the cellular basis for these changes, we characterized changes in mRNA and small RNA (sRNA) accumulation in S.sclerotiorum to infection by SsHV2-L. A total of 958 mRNAs and 835 sRNA-producing loci were altered after infection by SsHV2-L, among which >100 mRNAs were predicted to encode proteins involved in the metabolism and trafficking of carbohydrates and lipids. Both S. sclerotiorum endogenous and virus-derived sRNAs were predominantly 22 nt in length suggesting one dicer-like enzyme cleaves both. Novel classes of endogenous small RNAs were predicted, including phasiRNAs and tRNA-derived small RNAs. Moreover, S. sclerotiorum phasiRNAs, which were derived from noncoding RNAs and have the potential to regulate mRNA abundance in trans, showed differential accumulation due to virus infection. tRNA fragments did not accumulate differentially after hypovirus infection. Hence, in-depth analysis showed that infection of S. sclerotiorum by a hypovirulence-inducing hypovirus produced selective, large-scale reprogramming of mRNA and sRNA production.

Cell Wall Synthesis, Development of Hyphae and Metabolic Pathways Are Processes Potentially Regulated by MicroRNAs Produced Between the Morphological Stages of Paracoccidioides brasiliensis

MicroRNAs are molecules involved in post-transcriptional gene regulation. In pathogenic fungi, microRNAs have been described at different morphological stages by regulating targets involved in processes such as morphogenesis and energy production. Members of the Paracoccidioides complex are the main etiological agents of a systemic mycosis in Latin America. Fungi of the Paracoccidioides complex present a wide range of plasticity to colonize different niches. In response to environmental changes these fungi undergo a morphological switch, remodel their cellular metabolism and modulate structural cell wall components. However, the underlying mechanisms regulating the gene expression is not well understood. By using high performance sequencing and bioinformatics analyses, this work characterizes microRNAs produced by Paracoccidioides brasiliensis. Here, we demonstrated that the transcript encoding proteins involved in microRNA biogenesis were differentially expressed in each morphological stage. In addition, 49 microRNAs were identified in cDNA libraries with 44 differentially regulated among the libraries. Sixteen microRNAs were differentially regulated in comparison to the mycelium in the mycelium-to-yeast transition phase. The yeast parasitic phase revealed a complete remodeling of the expression of these small RNAs. Analyses of targets of the induced microRNAs, from the different libraries, revealed that these molecules may potentially regulate in the cell wall, by repressing genes involved in the synthesis and degradation of glucans and chitin. Furthermore, mRNAs involved in cellular metabolism and development were predicted to be regulated by microRNAs. Therefore, this work describes a putative post transcriptional regulation, mediated by microRNAs in P. brasiliensis and its influence on the adaptive processes of thermal dimorphic fungus.

Discovery of microRNA-like RNAs during early fruiting body development in the model mushroom Coprinopsis cinerea

Coprinopsis cinerea is a model mushroom particularly suited for the study of fungal fruiting body development and the evolution of multicellularity in fungi. While microRNAs (miRNAs) have been extensively studied in animals and plants for their essential roles in post-transcriptional regulation of gene expression, miRNAs in fungi are less well characterized and their potential roles in controlling mushroom development remain unknown. To identify miRNA-like RNAs (milRNAs) in C. cinerea and explore their expression patterns during the early developmental transition of mushroom development, small RNA libraries of vegetative mycelium and primordium were generated and putative milRNA candidates were identified following the standards of miRNA prediction in animals and plants. Two out of 22 novel predicted milRNAs, cci-milR-12c and cci-milR-13e-5p, were validated by northern blot and stem-loop reverse transcription real-time PCR. Cci-milR-12c was differentially expressed whereas the expression levels of cci-milR-13e-5p were similar in the two developmental stages. Target prediction of the validated milRNAs resulted in genes associated with fruiting body development, including pheromone, hydrophobin, cytochrome P450, and protein kinase. Essential genes for miRNA biogenesis, including three coding for Dicer-like (DCL), one for Argonaute (AGO), one for AGO-like and one for quelling deficient-2 (QDE-2) proteins, were also identified in the C. cinerea genome. Phylogenetic analysis showed that the DCL and AGO proteins of C. cinerea were more closely related to those in other basidiomycetes and ascomycetes than to those in animals and plants. Taken together, our findings provided the first evidence for milRNAs in the model mushroom and their potential roles in regulating fruiting body development. New information on the evolutionary relationship of milRNA biogenesis proteins across kingdoms has also provided new insights for guiding further functional and evolutionary studies of miRNAs.

Small RNAs from cereal powdery mildew pathogens may target host plant genes

Small RNAs (sRNAs) play a key role in eukaryotic gene regulation, for example by gene silencing via RNA interference (RNAi). The biogenesis of sRNAs depends on proteins that are generally conserved in all eukaryotic lineages, yet some species that lack part or all the components of the mechanism exist. Here we explored the presence of the RNAi machinery and its expression as well as the occurrence of sRNA candidates and their putative endogenous as well as host targets in phytopathogenic powdery mildew fungi. We focused on the species Blumeria graminis, which occurs in various specialized forms (formae speciales) that each have a strictly limited host range. B. graminis f. sp. hordei and B. graminis f. sp. tritici, colonizing barley and wheat, respectively, have genomes that are characterized by extensive gene loss. Nonetheless, we find that the RNAi machinery appears to be largely complete and expressed during infection. sRNA sequencing data enabled the identification of putative sRNAs in both pathogens. While a considerable part of the sRNA candidates have predicted target sites in endogenous genes and transposable elements, a small proportion appears to have targets in planta, suggesting potential cross-kingdom RNA transfer between powdery mildew fungi and their respective plant hosts.

MADS-Box Transcription Factor MadsA Regulates Dimorphic Transition, Conidiation, and Germination of Talaromyces marneffei

The opportunistic human pathogen Talaromyces marneffei exhibits a temperature-dependent dimorphic transition, which is closely related with its pathogenicity. This species grows as multinucleate mycelia that produce infectious conidia at 25°C, while undergoes a dimorphic transition to generate uninucleate yeast form cells at 37°C. The mechanisms of phenotype switching are not fully understood. The transcription factor madsA gene is a member of the MADS-box gene family. Previously, it was found that overexpression of madsA gene resulted in mycelial growth instead of yeast form at 37°C. In the current study, the madsA deletion mutant (ΔmadsA) and complemented strain (CMA) were constructed by genetic manipulation. We compared the phenotypes, growth, conidiation, conidial germination and susceptibility to stresses (including osmotic and oxidative) of the ΔmadsA with the wild-type (WT) and CMA strains. The results showed that the ΔmadsA displayed a faster process of the yeast-to-mycelium transition than the WT and CMA. In addition, the deletion of madsA led to a delay in conidia production and conidial germination. The tolerance of ΔmadsA conidia to hydrogen peroxide was better than that of the WT and CMA strains. Then, RNA-seq was performed to identify differences in gene expression between the ΔmadsA mutant and WT strain during the yeast phase, mycelium phase, yeast-to-mycelium transition and mycelium-to-yeast transition, respectively. Gene ontology functional enrichment analyses indicated that some important processes such as transmembrane transport, oxidation-reduction process, protein catabolic process and response to oxidative stress were affected by the madsA deletion. Together, our results suggest that madsA functions as a global regulator involved in the conidiation and germination, especially in the dimorphic transition of T. marneffei. Its roles in the survival, pathogenicity and transmission of T. marneffei require further investigation.

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.

Taxonomy and evolution of Aspergillus, Penicillium and Talaromyces in the omics era - Past, present and future

Aspergillus, Penicillium and Talaromyces are diverse, phenotypically polythetic genera encompassing species important to the environment, economy, biotechnology and medicine, causing significant social impacts. Taxonomic studies on these fungi are essential since they could provide invaluable information on their evolutionary relationships and define criteria for species recognition. With the advancement of various biological, biochemical and computational technologies, different approaches have been adopted for the taxonomy of Aspergillus, Penicillium and Talaromyces; for example, from traditional morphotyping, phenotyping to chemotyping (e.g. lipotyping, proteotypingand metabolotyping) and then mitogenotyping and/or phylotyping. Since different taxonomic approaches focus on different sets of characters of the organisms, various classification and identification schemes would result. In view of this, the consolidated species concept, which takes into account different types of characters, is recently accepted for taxonomic purposes and, together with the lately implemented 'One Fungus - One Name' policy, is expected to bring a more stable taxonomy for Aspergillus, Penicillium and Talaromyces, which could facilitate their evolutionary studies. The most significant taxonomic change for the three genera was the transfer of Penicillium subgenus Biverticillium to Talaromyces (e.g. the medically important thermally dimorphic 'P. marneffei' endemic in Southeast Asia is now named T. marneffei), leaving both Penicillium and Talaromyces as monophyletic genera. Several distantly related Aspergillus-like fungi were also segregated from Aspergillus, making this genus, containing members of both sexual and asexual morphs, monophyletic as well. In the current omics era, application of various state-of-the-art omics technologies is likely to provide comprehensive information on the evolution of Aspergillus, Penicillium and Talaromyces and a stable taxonomy will hopefully be achieved.

Dicer-Like Proteins Regulate Sexual Development via the Biogenesis of Perithecium-Specific MicroRNAs in a Plant Pathogenic Fungus Fusarium graminearum

Ascospores act as the primary inoculum of Fusarium graminearum, which causes the destructive disease Fusarium head blight (FHB), or scab. MicroRNAs (miRNAs) have been reported in the F. graminearum vegetative stage, and Fgdcl2 is involved in microRNA-like RNA (milRNA) biogenesis but has no major impact on vegetative growth, abiotic stress or pathogenesis. In the present study, we found that ascospore discharge was decreased in the Fgdcl1 deletion mutant, and completely blocked in the double-deletion mutant of Fgdcl1 and Fgdcl2. Besides, more immature asci were observed in the double-deletion mutant. Interestingly, the up-regulated differentially expressed genes (DEGs) common to ΔFgdcl1 and ΔFgdcl1/2 were related to ion transmembrane transporter and membrane components. The combination of small RNA and transcriptome sequencing with bioinformatics analysis predicted 143 novel milRNAs in wild-type perithecia, and 138 of these milRNAs partly or absolutely depended on Fgdcl1, while only 5 novel milRNAs were still obtained in the Fgdcl1 and Fgdcl2 double-deletion mutant. Furthermore, 117 potential target genes were predicted. Overall, Fgdcl1 and Fgdcl2 genes were partly functionally redundant in ascospore discharge and perithecium-specific milRNA generation in F. graminearum, and these perithecium-specific milRNAs play potential roles in sexual development.

MicroRNA-like RNAs in plant pathogenic fungus Fusarium oxysporum f. sp. niveum are involved in toxin gene expression fine tuning

MicroRNA-like RNAs (milRNAs) are short non-coding regulatory sRNAs which play an important role in regulating gene expression at the post-transcriptional level by targeting mRNAs for degradation or inhibiting protein translation. To explore the presence of milRNAs in Fusarium oxysporum f. sp. niveum (Fon) and analyze their expression at different propagules, two categories of sRNAs were identified from Fon hyphae and microconidia using illumina sequencing. A total of 650,960 and 561,114 unique sRNAs were obtained from the hyphae and microconidia samples. With a previously constructed pipeline to search for microRNAs, 74 and 56 milRNA candidates were identified in hyphae and microconidia, respectively, based on the short hairpin structure analysis. Global expression analysis showed an extensively differential expression of sRNAs between the two propagules. Altogether, 78 significantly differently expressed milRNAs were identified in two libraries. Target prediction revealed two interesting genes involved in trichothecene production, necrosis and ethylene-inducing peptide 1 (NEP1) biosynthesis and in silico analysis indicated that they were down-regulated by Fon-miR7696a-3p and Fon-miR6108a. The expression levels of these two milRNAs were further validated by qRT-PCR and the results were consistent. The negative correlation of the expression levels between these two milRNAs and their potential target genes imply that they play a role in trichothecene and NEP1 biosynthesis. And this negative regulation for toxin-related gene expression is more specific in microconidia. The present study provides the first large-scale characterization of milRNAs in Fon and the comparison between hyphae and microconidia propagules gives an insight into how milRNAs are involve in toxin biosynthesis.

Small RNA Functions Are Required for Growth and Development of Magnaporthe oryzae

RNA interference (RNAi) is conserved in eukaryotic organisms, and it has been well studied in many animal and plant species and some fungal species, yet it is not well studied in fungal plant pathogens. In the rice blast fungus Magnaporthe oryzae, we examined small RNA (sRNA) and their biogenesis in the context of growth and pathogenicity. Through genetic and genomic analyses, we demonstrate that loss of a single gene encoding Dicer, RNA-dependent RNA polymerase, or Argonaute reduces sRNA levels. These three proteins are required for the biogenesis of sRNA-matching genome-wide regions (coding regions, repeats, and intergenic regions). The loss of one Argonaute reduced both sRNA and fungal virulence on barley leaves. Transcriptome analysis of multiple mutants revealed that sRNA play an important role in transcriptional regulation of repeats and intergenic regions in M. oryzae. Together, these data support that M. oryzae sRNA regulate developmental processes including, fungal growth and virulence.

Puccinia striiformis f. sp. tritici microRNA-like RNA 1 (Pst-milR1), an important pathogenicity factor of Pst, impairs wheat resistance to Pst by suppressing the wheat pathogenesis-related 2 gene

Small RNAs (sRNAs), an important type of pathogenicity factor, contribute to impairing host immune responses. However, little is known about sRNAs in Puccinia striiformis f. sp. tritici (Pst), one of the most destructive pathogens of wheat (Triticum aestivum L.). Here, we report a novel microRNA-like RNA (milRNA) from Pst termed microRNA-like RNA 1 (Pst-milR1), which suppresses wheat defenses during wheat-Pst interactions. We identified Pst-milR1 as a novel milRNA in Pst. Biological prediction and co-transformation showed that Pst-milR1 takes part in cross-kingdom RNA interference (RNAi) events by binding the wheat pathogenesis-related 2 (PR2) gene. Silencing of the Pst-milR1 precursor resulted in increased wheat resistance to the virulent Pst isolate CYR31. PR2 knock-down plants increased the susceptibility of wheat to the avirulent Pst isolate CYR23. This suggests that Pst-milR1 represses the plant immune response by suppressing the expression of PR2. Taking our findings together, we postulate that Pst-milR1 is an important pathogenicity factor in Pst, which acts as an effector to suppress host immunity. Our results provide significant new insights into the pathogenicity of the stripe rust pathogen.