Evolving moldy murderers: Aspergillus section Fumigati as a model for studying the repeated evolution of fungal pathogenicity

Giant transposons promote strain heterogeneity in a major fungal pathogen

Fungal infections are difficult to prevent and treat in large part due to strain heterogeneity, which confounds diagnostic predictability. Yet, the genetic mechanisms driving strain-to-strain variation remain poorly understood. Here, we determined the extent to which Starships-giant transposons capable of mobilizing numerous fungal genes-generate genetic and phenotypic variability in the opportunistic human pathogen Aspergillus fumigatus. We analyzed 519 diverse strains, including 11 newly sequenced with long-read technology and multiple isolates of the same reference strain, to reveal 20 distinct Starships that are generating genomic heterogeneity over timescales relevant for experimental reproducibility. Starship-mobilized genes encode diverse functions, including known biofilm-related virulence factors and biosynthetic gene clusters, and many are differentially expressed during infection and antifungal exposure in a strain-specific manner. These findings support a new model of fungal evolution wherein Starships help generate variation in genome structure, gene content, and expression among fungal strains. Together, our results demonstrate that Starships are a previously hidden mechanism generating genotypic and, in turn, phenotypic heterogeneity in a major human fungal pathogen.IMPORTANCENo "one size fits all" option exists for treating fungal infections in large part due to genetic and phenotypic variability among strains. Accounting for strain heterogeneity is thus fundamental for developing efficacious treatments and strategies for safeguarding human health. Here, we report significant progress toward achieving this goal by uncovering a previously hidden mechanism generating heterogeneity in the human fungal pathogen Aspergillus fumigatus: giant transposons, called Starships, that span dozens of kilobases and mobilize fungal genes as cargo. By conducting a systematic investigation of these unusual transposons in a single fungal species, we demonstrate their contributions to population-level variation at the genome, pangenome, and transcriptome levels. The Starship compendium we develop will not only help predict variation introduced by these elements in laboratory experiments but will serve as a foundational resource for determining how Starships impact clinically relevant phenotypes, such as antifungal resistance and pathogenicity.

Giant transposons promote strain heterogeneity in a major fungal pathogen

Fungal infections are difficult to prevent and treat in large part due to strain heterogeneity. However, the genetic mechanisms driving pathogen variation remain poorly understood. Here, we determined the extent to which Starships-giant transposons capable of mobilizing numerous fungal genes-generate genetic and phenotypic variability in the human pathogen Aspergillus fumigatus. We analyzed 519 diverse strains, including 12 newly sequenced with long-read technology, to reveal 20 distinct Starships that are generating genomic heterogeneity over timescales potentially relevant for experimental reproducibility. Starship-mobilized genes encode diverse functions, including biofilm-related virulence factors and biosynthetic gene clusters, and many are differentially expressed during infection and antifungal exposure in a strain-specific manner. These findings support a new model of fungal evolution wherein Starships help generate variation in gene content and expression among fungal strains. Together, our results demonstrate that Starships are a previously hidden mechanism generating genotypic and, in turn, phenotypic heterogeneity in a major human fungal pathogen.

Aspergillus fumigatus conidial surface-associated proteome reveals factors for fungal evasion and host immunity modulation

Aspergillus fumigatus causes aspergillosis and relies on asexual spores (conidia) for initiating host infection. There is scarce information about A. fumigatus proteins involved in fungal evasion and host immunity modulation. Here we analysed the conidial surface proteome of A. fumigatus, two closely related non-pathogenic species, Aspergillus fischeri and Aspergillus oerlinghausenensis, as well as pathogenic Aspergillus lentulus, to identify such proteins. After identifying 62 proteins exclusively detected on the A. fumigatus conidial surface, we assessed null mutants for 42 genes encoding these proteins. Deletion of 33 of these genes altered susceptibility to macrophage, epithelial cells and cytokine production. Notably, a gene that encodes a putative glycosylasparaginase, modulating levels of the host proinflammatory cytokine IL-1β, is important for infection in an immunocompetent murine model of fungal disease. These results suggest that A. fumigatus conidial surface proteins are important for evasion and modulation of the immune response at the onset of fungal infection.

Genome-wide patterns of noncoding and protein-coding sequence variation in the major fungal pathogen Aspergillus fumigatus

Aspergillus fumigatus is a deadly fungal pathogen, responsible for >400,000 infections/year and high mortality rates. A. fumigatus strains exhibit variation in infection-relevant traits, including in their virulence. However, most A. fumigatus protein-coding genes, including those that modulate its virulence, are shared between A. fumigatus strains and closely related nonpathogenic relatives. We hypothesized that A. fumigatus genes exhibit substantial genetic variation in the noncoding regions immediately upstream to the start codons of genes, which could reflect differences in gene regulation between strains. To begin testing this hypothesis, we identified 5,812 single-copy orthologs across the genomes of 263 A. fumigatus strains. In general, A. fumigatus noncoding regions showed higher levels of sequence variation compared with their corresponding protein-coding regions. Focusing on 2,482 genes whose protein-coding sequence identity scores ranged between 75 and 99%, we identified 478 total genes with signatures of positive selection only in their noncoding regions and 65 total genes with signatures only in their protein-coding regions. Twenty-eight of the 478 noncoding regions and 5 of the 65 protein-coding regions under selection are associated with genes known to modulate A. fumigatus virulence. Noncoding region variation between A. fumigatus strains included single-nucleotide polymorphisms and insertions or deletions of at least a few nucleotides. These results show that noncoding regions of A. fumigatus genes harbor greater sequence variation than protein-coding regions, raising the hypothesis that this variation may contribute to A. fumigatus phenotypic heterogeneity.

Genomic factors shape carbon and nitrogen metabolic niche breadth across Saccharomycotina yeasts

Organisms exhibit extensive variation in ecological niche breadth, from very narrow (specialists) to very broad (generalists). Two general paradigms have been proposed to explain this variation: (i) trade-offs between performance efficiency and breadth and (ii) the joint influence of extrinsic (environmental) and intrinsic (genomic) factors. We assembled genomic, metabolic, and ecological data from nearly all known species of the ancient fungal subphylum Saccharomycotina (1154 yeast strains from 1051 species), grown in 24 different environmental conditions, to examine niche breadth evolution. We found that large differences in the breadth of carbon utilization traits between yeasts stem from intrinsic differences in genes encoding specific metabolic pathways, but we found limited evidence for trade-offs. These comprehensive data argue that intrinsic factors shape niche breadth variation in microbes.

Genome-wide patterns of non-coding sequence variation in the major fungal pathogen Aspergillus fumigatus

A.fumigatus is a deadly fungal pathogen, responsible for >400,000 infections/year and high mortality rates. A. fumigatus strains exhibit variation in infection-relevant traits, including in their virulence. However, most A. fumigatus protein-coding genes, including those that modulate its virulence, are shared between A. fumigatus strains and closely related non-pathogenic relatives. We hypothesized that A. fumigatus genes exhibit substantial genetic variation in the non-coding regions immediately upstream to the start codons of genes, which could reflect differences in gene regulation between strains. To begin testing this hypothesis, we identified 5,812 single-copy orthologs across the genomes of 263 A. fumigatus strains. A. fumigatus non-coding regions showed higher levels of sequence variation compared to their corresponding protein-coding regions. Specifically, we found that 1,274 non-coding regions exhibited <75% nucleotide sequence similarity (compared to 928 protein-coding regions) and 3,721 non-coding regions exhibited between 75% and 99% similarity (compared to 2,482 protein-coding regions) across strains. Only 817 non-coding regions exhibited ≥99% sequence similarity compared to 2,402 protein-coding regions. By examining 2,482 genes whose protein-coding sequence identity scores ranged between 75% and 99%, we identified 478 total genes with signatures of positive selection only in their non-coding regions and 65 total genes with signatures only in their protein-coding regions. 28 of the 478 non-coding regions and 5 of the 65 protein-coding regions under selection are associated with genes known to modulate A. fumigatus virulence. Non-coding region variation between A. fumigatus strains included single nucleotide polymorphisms and insertions or deletions of at least a few nucleotides. These results show that non-coding regions of A. fumigatus genes harbor greater sequence variation than protein-coding regions, raising the hypothesis that this variation may contribute to A. fumigatus phenotypic heterogeneity.

Modulatory immune responses in fungal infection associated with organ transplant - advancements, management, and challenges

Organ transplantation stands as a pivotal achievement in modern medicine, offering hope to individuals with end-stage organ diseases. Advancements in immunology led to improved organ transplant survival through the development of immunosuppressants, but this heightened susceptibility to fungal infections with nonspecific symptoms in recipients. This review aims to establish an intricate balance between immune responses and fungal infections in organ transplant recipients. It explores the fundamental immune mechanisms, recent advances in immune response dynamics, and strategies for immune modulation, encompassing responses to fungal infections, immunomodulatory approaches, diagnostics, treatment challenges, and management. Early diagnosis of fungal infections in transplant patients is emphasized with the understanding that innate immune responses could potentially reduce immunosuppression and promise efficient and safe immuno-modulating treatments. Advances in fungal research and genetic influences on immune-fungal interactions are underscored, as well as the potential of single-cell technologies integrated with machine learning for biomarker discovery. This review provides a snapshot of the complex interplay between immune responses and fungal infections in organ transplantation and underscores key research directions.

A phylogenetic approach to explore the Aspergillus fumigatus conidial surface-associated proteome and its role in pathogenesis

Aspergillus fumigatus, an important pulmonary fungal pathogen causing several diseases collectively called aspergillosis, relies on asexual spores (conidia) for initiating host infection. Here, we used a phylogenomic approach to compare proteins in the conidial surface of A. fumigatus, two closely related non-pathogenic species, Aspergillus fischeri and Aspergillus oerlinghausenensis, and the cryptic pathogen Aspergillus lentulus. After identifying 62 proteins uniquely expressed on the A. fumigatus conidial surface, we assessed null mutants for 42 genes encoding conidial proteins. Deletion of 33 of these genes altered susceptibility to macrophage killing, penetration and damage to epithelial cells, and cytokine production. Notably, a gene that encodes glycosylasparaginase, which modulates levels of the host pro-inflammatory cytokine IL-1β, is important for infection in an immunocompetent murine model of fungal disease. These results suggest that A. fumigatus conidial surface proteins and effectors are important for evasion and modulation of the immune response at the onset of fungal infection.

Genomic and ecological factors shaping specialism and generalism across an entire subphylum

Organisms exhibit extensive variation in ecological niche breadth, from very narrow (specialists) to very broad (generalists). Paradigms proposed to explain this variation either invoke trade-offs between performance efficiency and breadth or underlying intrinsic or extrinsic factors. We assembled genomic (1,154 yeast strains from 1,049 species), metabolic (quantitative measures of growth of 843 species in 24 conditions), and ecological (environmental ontology of 1,088 species) data from nearly all known species of the ancient fungal subphylum Saccharomycotina to examine niche breadth evolution. We found large interspecific differences in carbon breadth stem from intrinsic differences in genes encoding specific metabolic pathways but no evidence of trade-offs and a limited role of extrinsic ecological factors. These comprehensive data argue that intrinsic factors driving microbial niche breadth variation.

A phylogenetic approach to explore the Aspergillus fumigatus conidial surface-associated proteome and its role in pathogenesis

Aspergillus fumigatus , an important pulmonary fungal pathogen causing several diseases collectively called aspergillosis, relies on asexual spores or conidia for initiating host infection. Here, we used a phylogenomic approach to compare proteins in the conidial surface of A. fumigatus , two closely related non-pathogenic species, Aspergillus fischeri and Aspergillus oerlinghausenensis , and the cryptic pathogen Aspergillus lentulus . After identifying 62 proteins uniquely expressed on the A. fumigatus conidial surface, we deleted 42 genes encoding conidial proteins. We found deletion of 33 of these genes altered susceptibility to macrophage killing, penetration and damage to epithelial cells, and cytokine production. Notably, a gene that encodes glycosylasparaginase, which modulates levels of the host pro-inflammatory cytokine IL-1β, is important for infection in an immunocompetent murine model of fungal disease. These results suggest that A. fumigatus conidial surface proteins and effectors are important for evasion and modulation of the immune response at the onset of fungal infection.

Soil Aspergillus Species, Pathogenicity and Control Perspectives

Five Aspergillus sections have members that are established agricultural pests and producers of different metabolites, threatening global food safety. Most of these pathogenic Aspergillus species have been isolated from almost all major biomes. The soil remains the primary habitat for most of these cryptic fungi. This review explored some of the ecological attributes that have contributed immensely to the success of the pathogenicity of some members of the genus Aspergillus over time. Hence, the virulence factors of the genus Aspergillus, their ecology and others were reviewed. Furthermore, some biological control techniques were recommended. Pathogenic effects of Aspergillus species are entirely accidental; therefore, the virulence evolution prediction model in such species becomes a challenge, unlike their obligate parasite counterparts. In all, differences in virulence among organisms involved both conserved and species-specific genetic factors. If the impacts of climate change continue, new cryptic Aspergillus species will emerge and mycotoxin contamination risks will increase in all ecosystems, as these species can metabolically adjust to nutritional and biophysical challenges. As most of their gene clusters are silent, fungi continue to be a source of underexplored bioactive compounds. The World Soil Charter recognizes the relevance of soil biodiversity in supporting healthy soil functions. The question of how a balance may be struck between supporting healthy soil biodiversity and the control of toxic fungi species in the field to ensure food security is therefore pertinent. Numerous advanced strategies and biocontrol methods so far remain the most environmentally sustainable solution to the control of toxigenic fungi in the field.

Prolonged Subculturing of Aspergillus fumigatus on Galleria Extract Agar Results in Altered Virulence and Sensitivity to Antifungal Agents

Aspergillus fumigatus is an environmental saprophyte and opportunistic fungal pathogen of humans. The aim of the work presented here was to examine the effect of serially subculturing A. fumigatus on agar generated from Galleria mellonella larvae in order to characterize the alterations in the phenotypes that might occur. The passaged strains showed alterations in virulence, antifungal susceptibility, and in protein abundances that may indicate adaptation after 25 passages over 231 days on Galleria extract agar. Passaged strains demonstrated reduced virulence in G. mellonella larvae and increased tolerance to hemocyte-mediated killing, hydrogen peroxide, itraconazole, and amphotericin B. A label-free proteomic analysis of control and passaged A. fumigatus strains revealed a total of 3329 proteins, of which 1902 remained following filtration, and 32 proteins were statistically significant as well as differentially abundant. Proteins involved in the response to oxidative stress were altered in abundance in the passaged strain and included (S)-S-oxide reductase (+2.63-fold), developmental regulator FlbA (+2.27-fold), and histone H2A.Z (-1.82-fold). These results indicate that the prolonged subculturing of A. fumigatus on Galleria extract agar results in alterations in the susceptibility to antifungal agents and in the abundance of proteins associated with the oxidative stress response. The phenomenon may be a result of selection for survival in adverse conditions and highlight how A. fumigatus may adapt to tolerate the pulmonary immune response in cases of human infection.

A Simple CRISPR/Cas9 System for Efficiently Targeting Genes of Aspergillus Section Flavi Species, Aspergillus nidulans, Aspergillus fumigatus, Aspergillus terreus, and Aspergillus niger

For Aspergillus flavus, a pathogen of considerable economic and health concern, successful gene knockout work for more than a decade has relied nearly exclusively on using nonhomologous end-joining pathway (NHEJ)-deficient recipients via forced double-crossover recombination of homologous sequences. In this study, a simple CRISPR/Cas9 (clustered regularly interspaced short palindromic repeats/CRISPR-associated nuclease) genome editing system that gave extremely high (>95%) gene-targeting frequencies in A. flavus was developed. It contained a shortened Aspergillus nidulans AMA1 autonomously replicating sequence that maintained good transformation frequencies and Aspergillus oryzae ptrA as the selection marker for pyrithiamine resistance. Expression of the codon-optimized cas9 gene was driven by the A. nidulans gpdA promoter and trpC terminator. Expression of single guide RNA (sgRNA) cassettes was controlled by the A. flavus U6 promoter and terminator. The high transformation and gene-targeting frequencies of this system made generation of A. flavus gene knockouts with or without phenotypic changes effortless. Additionally, multiple-gene knockouts of A. flavus conidial pigment genes (olgA/copT/wA or olgA/yA/wA) were quickly generated by a sequential approach. Cotransforming sgRNA vectors targeting A. flavus kojA, yA, and wA gave 52%, 40%, and 8% of single-, double-, and triple-gene knockouts, respectively. The system was readily applicable to other section Flavi aspergilli (A. parasiticus, A. oryzae, A. sojae, A. nomius, A. bombycis, and A. pseudotamarii) with comparable transformation and gene-targeting efficiencies. Moreover, it gave satisfactory gene-targeting efficiencies (>90%) in A. nidulans (section Nidulantes), A. fumigatus (section Fumigati), A. terreus (section Terrei), and A. niger (section Nigri). It likely will have a broad application in aspergilli. IMPORTANCE CRISPR/Cas9 genome editing systems have been developed for many aspergilli. Reported gene-targeting efficiencies vary greatly and are dependent on delivery methods, repair mechanisms of induced double-stranded breaks, selection markers, and genetic backgrounds of transformation recipient strains. They are also mostly strain specific or species specific. This developed system is highly efficient and allows knocking out multiple genes in A. flavus efficiently either by sequential transformation or by cotransformation of individual sgRNA vectors if desired. It is readily applicable to section Flavi species and aspergilli in other sections ("section" is a taxonomic rank between genus and species). This cross-Aspergillus section system is for wild-type isolates and does not require homologous donor DNAs to be added, NHEJ-deficient strains to be created, or forced recycling of knockout recipients to be performed for multiple-gene targeting. Hence, it simplifies and expedites the gene-targeting process significantly.

The pan-genome of Aspergillus fumigatus provides a high-resolution view of its population structure revealing high levels of lineage-specific diversity driven by recombination

Aspergillus fumigatus is a deadly agent of human fungal disease where virulence heterogeneity is thought to be at least partially structured by genetic variation between strains. While population genomic analyses based on reference genome alignments offer valuable insights into how gene variants are distributed across populations, these approaches fail to capture intraspecific variation in genes absent from the reference genome. Pan-genomic analyses based on de novo assemblies offer a promising alternative to reference-based genomics with the potential to address the full genetic repertoire of a species. Here, we evaluate 260 genome sequences of A. fumigatus including 62 newly sequenced strains, using a combination of population genomics, phylogenomics, and pan-genomics. Our results offer a high-resolution assessment of population structure and recombination frequency, phylogenetically structured gene presence-absence variation, evidence for metabolic specificity, and the distribution of putative antifungal resistance genes. Although A. fumigatus disperses primarily via asexual conidia, we identified extraordinarily high levels of recombination with the lowest linkage disequilibrium decay value reported for any fungal species to date. We provide evidence for 3 primary populations of A. fumigatus, with recombination occurring only rarely between populations and often within them. These 3 populations are structured by both gene variation and distinct patterns of gene presence-absence with unique suites of accessory genes present exclusively in each clade. Accessory genes displayed functional enrichment for nitrogen and carbohydrate metabolism suggesting that populations may be stratified by environmental niche specialization. Similarly, the distribution of antifungal resistance genes and resistance alleles were often structured by phylogeny. Altogether, the pan-genome of A. fumigatus represents one of the largest fungal pan-genomes reported to date including many genes unrepresented in the Af293 reference genome. These results highlight the inadequacy of relying on a single-reference genome-based approach for evaluating intraspecific variation and the power of combined genomic approaches to elucidate population structure, genetic diversity, and putative ecological drivers of clinically relevant fungi.

Examination of Genome-Wide Ortholog Variation in Clinical and Environmental Isolates of the Fungal Pathogen Aspergillus fumigatus

Aspergillus fumigatus is both an environmental saprobe and an opportunistic human fungal pathogen. Knowledge of genomic variation across A. fumigatus isolates is essential for understanding the evolution of pathogenicity, virulence, and resistance to antifungal drugs. Here, we investigated 206 A. fumigatus isolates (133 clinical and 73 environmental isolates), aiming to identify genes with variable presence across isolates and test whether this variation was related to the clinical or environmental origin of isolates. The PanOrtho genome of A. fumigatus consists of 13,085 ortholog groups, of which 7,773 (59.4%) are shared by all isolates (core groups) and 5,312 (40.6%) vary in their gene presence across isolates (accessory groups plus singletons). Despite differences in the distribution of orthologs across all isolates, no significant differences were observed among clinical versus environmental isolates when phylogeny was accounted for. Orthologs that differ in their distribution across isolates tend to occur at low frequency and/or be restricted to specific isolates; thus, the degree of genomic conservation between orthologs of A. fumigatus is high. These results suggest that differences in the distribution of orthologs within A. fumigatus cannot be associated with the clinical or environmental origin of isolates. IMPORTANCE Aspergillus fumigatus is a cosmopolitan species of fungus responsible for thousands of cases of invasive disease annually. Clinical and environmental isolates of A. fumigatus exhibit extensive phenotypic differences, including differences related to virulence and antifungal drug resistance. A comprehensive survey of the genomic diversity present in A. fumigatus and its relationship to the clinical or environmental origin of isolates can contribute to the prediction of the mechanisms of evolution and infection of the species. Our results suggest that there is no significant variation in ortholog distribution between clinical and environmental isolates when accounting for evolutionary history. The work supports the hypothesis that environmental and clinical isolates of A. fumigatus do not differ in their gene contents.

Diversity of Parengyodontium spp. strains isolated from the cultural heritage environment: Phylogenetic diversity, phenotypical diversity, and occurrence

Parengyodontium album is a fungal species that frequently occurs in the cultural heritage environment. Although three subclades were initially described in the species, no study has sought to determine the occurrence of each subclade in the cultural heritage context. These subclades are easily distinguishable phylogenetically, but their morphological identification is more difficult. Eighteen strains isolated from different cultural sites and initially identified as P. album were studied phylogenetically, morphologically, and in terms of their susceptibility to econazole nitrate 0.2%, an antifungal product used as preservation treatment in cultural heritage domain. The phylogenetic study revealed that all studied strains belonged to P. album subclade 1 or P. torokii (P. album subclade 3) and none belonged to P. album subclade 2. The morphological study revealed the best characteristics to differentiate the three subclades/species, namely, the ability of the strains to grow at 32 C and 35 C on potato dextrose agar (PDA) medium and the shape of conidia. Finally, the strains displayed variable susceptibilities to econazole nitrate, with no apparent link to any particular subclade/species.

Extensive non-coding sequence divergence between the major human pathogen Aspergillus fumigatus and its relatives

Invasive aspergillosis is a deadly fungal disease; more than 400,000 patients are infected worldwide each year and the mortality rate can be as high as 50-95%. Of the ~450 species in the genus Aspergillus only a few are known to be clinically relevant, with the major pathogen Aspergillus fumigatus being responsible for ~50% of all invasive mold infections. Genomic comparisons between A. fumigatus and other Aspergillus species have historically focused on protein-coding regions. However, most A. fumigatus genes, including those that modulate its virulence, are also present in other pathogenic and non-pathogenic closely related species. Our hypothesis is that differential gene regulation - mediated through the non-coding regions upstream of genes' first codon - contributes to A. fumigatus pathogenicity. To begin testing this, we compared non-coding regions upstream of the first codon of single-copy orthologous genes from the two A. fumigatus reference strains Af293 and A1163 and eight closely related Aspergillus section Fumigati species. We found that these non-coding regions showed extensive sequence variation and lack of homology across species. By examining the evolutionary rates of both protein-coding and non-coding regions in a subset of orthologous genes with highly conserved non-coding regions across the phylogeny, we identified 418 genes, including 25 genes known to modulate A. fumigatus virulence, whose non-coding regions exhibit a different rate of evolution in A. fumigatus. Examination of sequence alignments of these non-coding regions revealed numerous instances of insertions, deletions, and other types of mutations of at least a few nucleotides in A. fumigatus compared to its close relatives. These results show that closely related Aspergillus species that vary greatly in their pathogenicity exhibit extensive non-coding sequence variation and identify numerous changes in non-coding regions of A. fumigatus genes known to contribute to virulence.

Evolution of the human pathogenic lifestyle in fungi

Fungal pathogens cause more than a billion human infections every year, resulting in more than 1.6 million deaths annually. Understanding the natural history and evolutionary ecology of fungi is helping us understand how disease-relevant traits have repeatedly evolved. Different types and mechanisms of genetic variation have contributed to the evolution of fungal pathogenicity and specific genetic differences distinguish pathogens from non-pathogens. Insights into the traits, genetic elements, and genetic and ecological mechanisms that contribute to the evolution of fungal pathogenicity are crucial for developing strategies to both predict emergence of fungal pathogens and develop drugs to combat them.

Fungal cell wall components modulate our immune system

Invasive fungal infections remain highly problematic for human health. Collectively, they account for more than 1 million deaths a year in addition to more than 100 million mucosal infections and 1 billion skin infections. To be able to make progress it is important to understand the pathobiology of fungal interactions with the immune system. Here, we highlight new advancements pointing out the pivotal role of fungal cell wall components (β-glucan, mannan, galactosaminogalactan and melanin) in modulating host immunity and discuss how these open new opportunities for the development of immunomodulatory strategies to combat deadly fungal infectious diseases.

Examining Signatures of Natural Selection in Antifungal Resistance Genes Across Aspergillus Fungi

Certain Aspergillus fungi cause aspergillosis, a set of diseases that typically affect immunocompromised individuals. Most cases of aspergillosis are caused by Aspergillus fumigatus, which infects millions of people annually. Some closely related so-called cryptic species, such as Aspergillus lentulus, can also cause aspergillosis, albeit at lower frequencies, and they are also clinically relevant. Few antifungal drugs are currently available for treating aspergillosis and there is increasing worldwide concern about the presence of antifungal drug resistance in Aspergillus species. Furthermore, isolates from both A. fumigatus and other Aspergillus pathogens exhibit substantial heterogeneity in their antifungal drug resistance profiles. To gain insights into the evolution of antifungal drug resistance genes in Aspergillus, we investigated signatures of positive selection in 41 genes known to be involved in drug resistance across 42 susceptible and resistant isolates from 12 Aspergillus section Fumigati species. Using codon-based site models of sequence evolution, we identified ten genes that contain 43 sites with signatures of ancient positive selection across our set of species. None of the sites that have experienced positive selection overlap with sites previously reported to be involved in drug resistance. These results identify sites that likely experienced ancient positive selection in Aspergillus genes involved in resistance to antifungal drugs and suggest that historical selective pressures on these genes likely differ from any current selective pressures imposed by antifungal drugs.

Genomic and Phenotypic Analysis of COVID-19-Associated Pulmonary Aspergillosis Isolates of Aspergillus fumigatus

The ongoing global pandemic caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for coronavirus disease 2019 (COVID-19), first described in Wuhan, China. A subset of COVID-19 patients has been reported to have acquired secondary infections by microbial pathogens, such as opportunistic fungal pathogens from the genus Aspergillus. To gain insight into COVID-19-associated pulmonary aspergillosis (CAPA), we analyzed the genomes and characterized the phenotypic profiles of four CAPA isolates of Aspergillus fumigatus obtained from patients treated in the area of North Rhine-Westphalia, Germany. By examining the mutational spectrum of single nucleotide polymorphisms, insertion-deletion polymorphisms, and copy number variants among 206 genes known to modulate A. fumigatus virulence, we found that CAPA isolate genomes do not exhibit significant differences from the genome of the Af293 reference strain. By examining a number of factors, including virulence in an invertebrate moth model, growth in the presence of osmotic, cell wall, and oxidative stressors, secondary metabolite biosynthesis, and the MIC of antifungal drugs, we found that CAPA isolates were generally, but not always, similar to A. fumigatus reference strains Af293 and CEA17. Notably, CAPA isolate D had more putative loss-of-function mutations in genes known to increase virulence when deleted. Moreover, CAPA isolate D was significantly more virulent than the other three CAPA isolates and the A. fumigatus reference strains Af293 and CEA17, but similarly virulent to two other clinical strains of A. fumigatus. These findings expand our understanding of the genomic and phenotypic characteristics of isolates that cause CAPA. IMPORTANCE The global pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the etiological agent of coronavirus disease 2019 (COVID-19), has already killed millions of people. COVID-19 patient outcome can be further complicated by secondary infections, such as COVID-19-associated pulmonary aspergillosis (CAPA). CAPA is caused by Aspergillus fungal pathogens, but there is little information about the genomic and phenotypic characteristics of CAPA isolates. We conducted genome sequencing and extensive phenotyping of four CAPA isolates of Aspergillus fumigatus from Germany. We found that CAPA isolates were often, but not always, similar to other reference strains of A. fumigatus across 206 genetic determinants of infection-relevant phenotypes, including virulence. For example, CAPA isolate D was more virulent than other CAPA isolates and reference strains in an invertebrate model of fungal disease, but similarly virulent to two other clinical strains. These results expand our understanding of COVID-19-associated pulmonary aspergillosis.

An evolutionary genomic approach reveals both conserved and species-specific genetic elements related to human disease in closely related Aspergillus fungi

Aspergillosis is an important opportunistic human disease caused by filamentous fungi in the genus Aspergillus. Roughly 70% of infections are caused by Aspergillus fumigatus, with the rest stemming from approximately a dozen other Aspergillus species. Several of these pathogens are closely related to A. fumigatus and belong in the same taxonomic section, section Fumigati. Pathogenic species are frequently most closely related to nonpathogenic ones, suggesting Aspergillus pathogenicity evolved multiple times independently. To understand the repeated evolution of Aspergillus pathogenicity, we performed comparative genomic analyses on 18 strains from 13 species, including 8 species in section Fumigati, which aimed to identify genes, both ones previously connected to virulence as well as ones never before implicated, whose evolution differs between pathogens and nonpathogens. We found that most genes were present in all species, including approximately half of those previously connected to virulence, but a few genes were section- or species-specific. Evolutionary rate analyses identified over 1700 genes whose evolutionary rate differed between pathogens and nonpathogens and dozens of genes whose rates differed between specific pathogens and the rest of the taxa. Functional testing of deletion mutants of 17 transcription factor-encoding genes whose evolution differed between pathogens and nonpathogens identified eight genes that affect either fungal survival in a model of phagocytic killing, host survival in an animal model of fungal disease, or both. These results suggest that the evolution of pathogenicity in Aspergillus involved both conserved and species-specific genetic elements, illustrating how an evolutionary genomic approach informs the study of fungal disease.

Interspecies Genomic Variation and Transcriptional Activeness of Secondary Metabolism-Related Genes in Aspergillus Section Fumigati

Filamentous fungi produce various bioactive compounds that are biosynthesized by sets of proteins encoded in biosynthesis gene clusters (BGCs). For an unknown reason, many BGCs are transcriptionally silent in laboratory conditions, which has hampered the discovery of novel fungal compounds. The transcriptional reactiveness of fungal secondary metabolism is not fully understood. To gain the comprehensive view, we conducted comparative genomic and transcriptomic analyses of nine closely-related species of Aspergillus section Fumigati (A. fumigatus, A. fumigatiaffinis, A. novofumigatus, A. thermomutatus, A. viridinutans, A. pseudoviridinutans, A. lentulus, A. udagawae, and Neosartorya fischeri). For expanding our knowledge, we newly sequenced genomes of A. viridinutans and A. pseudoviridinutans, and reassembled and reannotated the previously released genomes of A. lentulus and A. udagawae. Between 34 and 84 secondary metabolite (SM) backbone genes were identified in the genomes of these nine respective species, with 8.7-51.2% being unique to the species. A total of 247 SM backbone gene types were identified in the nine fungi. Ten BGCs are shared by all nine species. Transcriptomic analysis using A. fumigatus, A. lentulus, A. udagawae, A. viridinutans, and N. fischeri was conducted to compare expression levels of all SM backbone genes in four different culture conditions; 32-83% of SM backbone genes in these species were not expressed in the tested conditions, which reconfirmed that large part of fungal SM genes are hard to be expressed. The species-unique SM genes of the five species were expressed with lower frequency (18.8% in total) than the SM genes that are conserved in all five species (56%). These results suggest that the expression tendency of BGCs is correlated with their interspecies distribution pattern. Our findings increase understanding of the evolutionary processes associated with the regulation of fungal secondary metabolism.

Genomic and phenotypic analysis of COVID-19-associated pulmonary aspergillosis isolates of Aspergillus fumigatus

The ongoing global pandemic caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the coronavirus disease 2019 (COVID-19) first described from Wuhan, China. A subset of COVID-19 patients has been reported to have acquired secondary infections by microbial pathogens, such as fungal opportunistic pathogens from the genus Aspergillus . To gain insight into COVID-19 associated pulmonary aspergillosis (CAPA), we analyzed the genomes and characterized the phenotypic profiles of four CAPA isolates of Aspergillus fumigatus obtained from patients treated in the area of North Rhine-Westphalia, Germany. By examining the mutational spectrum of single nucleotide polymorphisms, insertion-deletion polymorphisms, and copy number variants among 206 genes known to modulate A. fumigatus virulence, we found that CAPA isolate genomes do not exhibit major differences from the genome of the Af293 reference strain. By examining virulence in an invertebrate moth model, growth in the presence of osmotic, cell wall, and oxidative stressors, and the minimum inhibitory concentration of antifungal drugs, we found that CAPA isolates were generally, but not always, similar to A. fumigatus reference strains Af293 and CEA17. Notably, CAPA isolate D had more putative loss of function mutations in genes known to increase virulence when deleted (e.g., in the FLEA gene, which encodes a lectin recognized by macrophages). Moreover, CAPA isolate D was significantly more virulent than the other three CAPA isolates and the A. fumigatus reference strains tested. These findings expand our understanding of the genomic and phenotypic characteristics of isolates that cause CAPA.

Genome-scale phylogeny and contrasting modes of genome evolution in the fungal phylum Ascomycota

Ascomycota, the largest and most well-studied phylum of fungi, contains three subphyla: Saccharomycotina (budding yeasts), Pezizomycotina (filamentous fungi), and Taphrinomycotina (fission yeasts). Despite its importance, we lack a comprehensive genome-scale phylogeny or understanding of the similarities and differences in the mode of genome evolution within this phylum. By examining 1107 genomes from Saccharomycotina (332), Pezizomycotina (761), and Taphrinomycotina (14) species, we inferred a robust genome-wide phylogeny that resolves several contentious relationships and estimated that the Ascomycota last common ancestor likely originated in the Ediacaran period. Comparisons of genomic properties revealed that Saccharomycotina and Pezizomycotina differ greatly in their genome properties and enabled inference of the direction of evolutionary change. The Saccharomycotina typically have smaller genomes, lower guanine-cytosine contents, lower numbers of genes, and higher rates of molecular sequence evolution compared with Pezizomycotina. These results provide a robust evolutionary framework for understanding the diversity and ecological lifestyles of the largest fungal phylum.

Comparative Genomics Analyses of Lifestyle Transitions at the Origin of an Invasive Fungal Pathogen in the Genus Cryphonectria

Forest and agroecosystems, as well as animal and human health, are threatened by emerging pathogens. Following decimation of chestnuts in the United States, the fungal pathogen Cryphonectria parasitica colonized Europe. After establishment, the pathogen population gave rise to a highly successful lineage that spread rapidly across the continent. Core to our understanding of what makes a successful pathogen is the genetic repertoire enabling the colonization and exploitation of host species. Here, we have assembled >100 genomes across two related genera to identify key genomic determinants leading to the emergence of chestnut blight. We found subtle yet highly specific changes in the transition from saprotrophy to latent pathogenicity mostly determined by enzymes involved in carbohydrate metabolism. Large-scale genomic analyses of genes underlying key nutrition modes can facilitate the detection of species with the potential to emerge as pathogens.

Emerging fungal pathogens are a threat to forest and agroecosystems, as well as animal and human health. How pathogens evolve from nonpathogenic ancestors is still poorly understood, making the prediction of future outbreaks challenging. Most pathogens have evolved lifestyle adaptations, which were enabled by specific changes in the gene content of the species. Hence, understanding transitions in the functions encoded by genomes gives valuable insight into the evolution of pathogenicity. Here, we studied lifestyle evolution in the genus Cryphonectria, including the prominent invasive pathogen Cryphonectria parasitica, the causal agent of chestnut blight on Castanea species. We assembled and compared the genomes of pathogenic and putatively nonpathogenic Cryphonectria species, as well as sister group pathogens in the family Cryphonectriaceae (Diaporthales, Ascomycetes), to investigate the evolution of genome size and gene content. We found a striking loss of genes associated with carbohydrate metabolism (CAZymes) in C. parasitica compared to other Cryphonectriaceae. Despite substantial CAZyme gene loss, experimental data suggest that C. parasitica has retained wood colonization abilities shared with other Cryphonectria species. Putative effectors substantially varied in number, cysteine content, and protein length among species. In contrast, secondary metabolite gene clusters show a high degree of conservation within the genus. Overall, our results underpin the recent lifestyle transition of C. parasitica toward a more pathogenic lifestyle. Our findings suggest that a CAZyme loss may have promoted pathogenicity of C. parasitica on Castanea species. Analyzing gene complements underlying key nutrition modes can facilitate the detection of species with the potential to emerge as pathogens.

IMPORTANCE Forest and agroecosystems, as well as animal and human health, are threatened by emerging pathogens. Following decimation of chestnuts in the United States, the fungal pathogen Cryphonectria parasitica colonized Europe. After establishment, the pathogen population gave rise to a highly successful lineage that spread rapidly across the continent. Core to our understanding of what makes a successful pathogen is the genetic repertoire enabling the colonization and exploitation of host species. Here, we have assembled >100 genomes across two related genera to identify key genomic determinants leading to the emergence of chestnut blight. We found subtle yet highly specific changes in the transition from saprotrophy to latent pathogenicity mostly determined by enzymes involved in carbohydrate metabolism. Large-scale genomic analyses of genes underlying key nutrition modes can facilitate the detection of species with the potential to emerge as pathogens.

Variation Among Biosynthetic Gene Clusters, Secondary Metabolite Profiles, and Cards of Virulence Across Aspergillus Species

Aspergillus fumigatus is a major human pathogen. In contrast, Aspergillus fischeri and the recently described Aspergillus oerlinghausenensis, the two species most closely related to A. fumigatus, are not known to be pathogenic. Some of the genetic determinants of virulence (or "cards of virulence") that A. fumigatus possesses are secondary metabolites that impair the host immune system, protect from host immune cell attacks, or acquire key nutrients. To examine whether secondary metabolism-associated cards of virulence vary between these species, we conducted extensive genomic and secondary metabolite profiling analyses of multiple A. fumigatus, one A. oerlinghausenensis, and multiple A. fischeri strains. We identified two cards of virulence (gliotoxin and fumitremorgin) shared by all three species and three cards of virulence (trypacidin, pseurotin, and fumagillin) that are variable. For example, we found that all species and strains examined biosynthesized gliotoxin, which is known to contribute to virulence, consistent with the conservation of the gliotoxin biosynthetic gene cluster (BGC) across genomes. For other secondary metabolites, such as fumitremorgin, a modulator of host biology, we found that all species produced the metabolite but that there was strain heterogeneity in its production within species. Finally, species differed in their biosynthesis of fumagillin and pseurotin, both contributors to host tissue damage during invasive aspergillosis. A. fumigatus biosynthesized fumagillin and pseurotin, while A. oerlinghausenensis biosynthesized fumagillin and A. fischeri biosynthesized neither. These biochemical differences were reflected in sequence divergence of the intertwined fumagillin/pseurotin BGCs across genomes. These results delineate the similarities and differences in secondary metabolism-associated cards of virulence between a major fungal pathogen and its nonpathogenic closest relatives, shedding light onto the genetic and phenotypic changes associated with the evolution of fungal pathogenicity.

Draft Genome Sequences of Four Aspergillus Section Fumigati Clinical Strains

Aspergillus fungi in section Fumigati include important human pathogens. Here, we sequenced the genomes of two strains of Aspergillus hiratsukae and two strains of Aspergillus felis. The average genome sizes are 29.5 Mb for A. hiratsukae and 31.8 Mb for A. felis.

Aspergillus fungi in section Fumigati include important human pathogens. Here, we sequenced the genomes of two strains of Aspergillus hiratsukae and two strains of Aspergillus felis. The average genome sizes are 29.5 Mb for A. hiratsukae and 31.8 Mb for A. felis.

Pathogenic Allodiploid Hybrids of Aspergillus Fungi

Interspecific hybridization substantially alters genotypes and phenotypes and can give rise to new lineages. Hybrid isolates that differ from their parental species in infection-relevant traits have been observed in several human-pathogenic yeasts and plant-pathogenic filamentous fungi but have yet to be found in human-pathogenic filamentous fungi. We discovered 6 clinical isolates from patients with aspergillosis originally identified as Aspergillus nidulans (section Nidulantes) that are actually allodiploid hybrids formed by the fusion of Aspergillus spinulosporus with an unknown close relative of Aspergillus quadrilineatus, both in section Nidulantes. Evolutionary genomic analyses revealed that these isolates belong to Aspergillus latus, an allodiploid hybrid species. Characterization of diverse infection-relevant traits further showed that A. latus hybrid isolates are genomically and phenotypically heterogeneous but also differ from A. nidulans, A. spinulosporus, and A. quadrilineatus. These results suggest that allodiploid hybridization contributes to the genomic and phenotypic diversity of filamentous fungal pathogens of humans.

Genomic and Phenotypic Heterogeneity of Clinical Isolates of the Human Pathogens Aspergillus fumigatus, Aspergillus lentulus, and Aspergillus fumigatiaffinis

Fungal pathogens are a global threat to human health. For example, fungi from the genus Aspergillus cause a spectrum of diseases collectively known as aspergillosis. Most of the >200,000 life-threatening aspergillosis infections per year worldwide are caused by Aspergillus fumigatus. Recently, molecular typing techniques have revealed that aspergillosis can also be caused by organisms that are phenotypically similar to A. fumigatus but genetically distinct, such as Aspergillus lentulus and Aspergillus fumigatiaffinis. Importantly, some of these so-called cryptic species are thought to exhibit different virulence and drug susceptibility profiles than A. fumigatus, however, our understanding of their biology and pathogenic potential has been stymied by the lack of genome sequences and phenotypic profiling of multiple clinical strains. To fill this gap, we phenotypically characterized the virulence and drug susceptibility of 15 clinical strains of A. fumigatus, A. lentulus, and A. fumigatiaffinis from Spain and sequenced their genomes. We found heterogeneity in drug susceptibility across species and strains. We further found heterogeneity in virulence within each species but no significant differences in the virulence profiles between the three species. Genes known to influence drug susceptibility (cyp51A and fks1) vary in paralog number and sequence among these species and strains and correlate with differences in drug susceptibility. Similarly, genes known to be important for virulence in A. fumigatus showed variability in number of paralogs across strains and across species. Characterization of the genomic similarities and differences of clinical strains of A. lentulus, A. fumigatiaffinis, and A. fumigatus that vary in disease-relevant traits will advance our understanding of the variance in pathogenicity between Aspergillus species and strains that are collectively responsible for the vast majority of aspergillosis infections in humans.