New state-of-the-art techniques in sequencing offer valuable tools in both detection of mycobiota and in understanding of the molecular mechanisms of resistance against antifungal compounds and virulence. Introduction of new sequencing platform with enhanced capacity and a reduction in costs for sequence analysis provides a potential powerful tool in mycological diagnosis and research. In this review, we summarize the applications of next-generation sequencing techniques in mycology.

In the last few years, increasing numbers of viruses infecting fungi have been identified. In this study, we used an in silico approach for the analysis of deep RNA sequencing data in order to discover and characterize putative genomic ssRNA or dsRNA mycovirus sequences in Aspergillus fumigatus. RNA sequencing reads of A. fumigatus strains were mapped against the A. fumigatus Af293 reference genome. Unmapped reads were collected for de novo assembly. Contigs were analyzed by Blastx comparison with a mycovirus protein database. Assembled viral genomes were used as template for remapping of RNA sequencing reads. In total, deep RNA sequencing results from 11 A. fumigatus strains were analyzed for the presence of mycoviral genomic RNAs. In 9 out of 11 strains, putative mycoviral RNA genomes were identified. Three strains were infected with two different mycovirus species. Two strains were infected with Aspergillus fumigatus polymycovirus type-1 (AfuPmV-1). Four strains contained fully recovered genomic RNA of unknown narna-like viruses designated as Aspergillus fumigatus narnavirus-1 and Aspergillus fumigatus narnavirus-2 (AfuNV-1 and AfuNV-2). Both viruses showed 38% amino acid sequence identity to Beihai narna-like virus-21. Three strains contained partially recovered genomic RNA of an unknown narna-like virus. Two strains contained fully recovered genomic RNAs of an unknown partitivirus designated as Aspergillus fumigatus partitivirus-2 (AfuPV-2) which showed 50% amino acid sequence identity to Alternaria alternata partitivirus-1. Finally, one strain contained fully recovered genomic RNA of an unknown mitovirus designated as Aspergillus fumigatus mitovirus-1 (AfuMV-1) which showed 34% amino acid sequence identity to Sclerotina sclerotiorum mitovirus. In silico analysis of deep RNA sequencing results showed that a majority of the A. fumigatus strains used here were infected with mycoviruses. Four novel A. fumigatus RNA mycoviruses could be identified: two different Aspergillus fumigatus narna-like viruses, one Aspergillus fumigatus partitivirus, and one Aspergillus fumigatus mitovirus.

Aspergillus fumigatus is a saprophytic fungus that extensively produces conidia. These microscopic asexually reproductive structures are small enough to reach the lungs. Germination of conidia followed by hyphal growth inside human lungs is a key step in the establishment of infection in immunocompromised patients. RNA-Seq was used to analyze the transcriptome of dormant and germinating A. fumigatus conidia. Construction of a gene co-expression network revealed four gene clusters (modules) correlated with a growth phase (dormant, isotropic growth, polarized growth). Transcripts levels of genes encoding for secondary metabolites were high in dormant conidia. During isotropic growth, transcript levels of genes involved in cell wall modifications increased. Two modules encoding for growth and cell cycle/DNA processing were associated with polarized growth. In addition, the co-expression network was used to identify highly connected intermodular hub genes. These genes may have a pivotal role in the respective module and could therefore be compelling therapeutic targets. Generally, cell wall remodeling is an important process during isotropic and polarized growth, characterized by an increase of transcripts coding for hyphal growth and cell cycle/DNA processing when polarized growth is initiated.

Aspergillus flavus is a pathogenic fungal species that can cause pulmonary aspergillosis, and triazole compounds are used for the treatment of these infections. Prolonged exposure to azoles may select for compensatory mutations in the A. flavus genome, resulting in azole resistance. Here, we characterize a series of 11 isogenic A. flavus strains isolated from a patient with pulmonary aspergillosis. Over a period of three months, the initially azole-susceptible strain developed itraconazole and voriconazole resistance. Short tandem repeat analysis and whole-genome sequencing revealed the high genetic relatedness of all isolates, indicating an infection with one single isolate. In contrast, the isolates were macroscopically highly diverse, suggesting an adaptation to the environment due to (epi)genetic changes. The whole-genome sequencing of susceptible and azole-resistant strains showed a number of mutations that might be associated with azole resistance. The majority of resistant strains contain a Y119F mutation in the Cyp51A gene, which corresponds to the Y121F mutation found in A. fumigatus. One azole-resistant strain demonstrated a divergent set of mutations, including a V99A mutation in a major facilitator superfamily (MSF) multidrug transporter (AFLA 083950).

The opportunistic pathogen Aspergillus fumigatus is found on all continents and thrives in soil and agricultural environments. Its ability to readily adapt to novel environments and to produce billions of spores led to the spread of azole-resistant A. fumigatus across the globe, posing a threat to many immunocompromised patients, including critically ill patients with severe influenza or COVID-19. In our study, we sought to compare the adaptational response to azoles from A. fumigatus isolates that differ in azole susceptibility and genetic background. To gain more insight into how short-term adaptation to stressful azole compounds is managed through gene expression, we conducted an RNA-sequencing study on the response of A. fumigatus to itraconazole and the newest clinically approved azole, isavuconazole. We observed many similarities in ergosterol biosynthesis up-regulation across isolates, with the exception of the pan-azole-resistant isolate, which showed very little differential regulation in comparison to other isolates. Additionally, we found differential regulation of membrane efflux transporters, secondary metabolites, iron metabolism, and various stress response and cell signaling mechanisms.