We present the genome sequences of a new clinical isolate of the important human pathogen, Aspergillus fumigatus, A1163, and two closely related but rarely pathogenic species, Neosartorya fischeri NRRL181 and Aspergillus clavatus NRRL1. Comparative genomic analysis of A1163 with the recently sequenced A. fumigatus isolate Af293 has identified core, variable and up to 2% unique genes in each genome. While the core genes are 99.8% identical at the nucleotide level, identity for variable genes can be as low 40%. The most divergent loci appear to contain heterokaryon incompatibility (het) genes associated with fungal programmed cell death such as developmental regulator rosA. Cross-species comparison has revealed that 8.5%, 13.5% and 12.6%, respectively, of A. fumigatus, N. fischeri and A. clavatus genes are species-specific. These genes are significantly smaller in size than core genes, contain fewer exons and exhibit a subtelomeric bias. Most of them cluster together in 13 chromosomal islands, which are enriched for pseudogenes, transposons and other repetitive elements. At least 20% of A. fumigatus-specific genes appear to be functional and involved in carbohydrate and chitin catabolism, transport, detoxification, secondary metabolism and other functions that may facilitate the adaptation to heterogeneous environments such as soil or a mammalian host. Contrary to what was suggested previously, their origin cannot be attributed to horizontal gene transfer (HGT), but instead is likely to involve duplication, diversification and differential gene loss (DDL). The role of duplication in the origin of lineage-specific genes is further underlined by the discovery of genomic islands that seem to function as designated “gene dumps” and, perhaps, simultaneously, as “gene factories”.

BackgroundThe fungal genus Aspergillus is of critical importance to humankind. Species include those with industrial applications, important pathogens of humans, animals and crops, a source of potent carcinogenic contaminants of food, and an important genetic model. The genome sequences of eight aspergilli have already been explored to investigate aspects of fungal biology, raising questions about evolution and specialization within this genus.ResultsWe have generated genome sequences for ten novel, highly diverse Aspergillus species and compared these in detail to sister and more distant genera. Comparative studies of key aspects of fungal biology, including primary and secondary metabolism, stress response, biomass degradation, and signal transduction, revealed both conservation and diversity among the species. Observed genomic differences were validated with experimental studies. This revealed several highlights, such as the potential for sex in asexual species, organic acid production genes being a key feature of black aspergilli, alternative approaches for degrading plant biomass, and indications for the genetic basis of stress response. A genome-wide phylogenetic analysis demonstrated in detail the relationship of the newly genome sequenced species with other aspergilli.ConclusionsMany aspects of biological differences between fungal species cannot be explained by current knowledge obtained from genome sequences. The comparative genomics and experimental study, presented here, allows for the first time a genus-wide view of the biological diversity of the aspergilli and in many, but not all, cases linked genome differences to phenotype. Insights gained could be exploited for biotechnological and medical applications of fungi.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-017-1151-0) contains supplementary material, which is available to authorized users.

SummaryThe protein phosphatase calcineurin is an important mediator connecting calcium-dependent signalling to various cellular responses in multiple organisms. In fungi calcineurin acts largely through regulating Crz1p-like transcription factors. Here we characterize an Aspergillus fumigatus CRZ1 homologue, CrzA and demonstrate its mediation of cellular tolerance to increased concentrations of calcium and manganese. In addition to acute sensitivitiy to these ions, and decreased conidiation, the crzA null mutant suffers altered expression of calcium transporter mRNAs under high concentrations of calcium, and loss of virulence when compared with the corresponding complemented and wild-type strains. We use multiple expression analyses to probe the transcriptional basis of A. fumigatus calcium tolerance identifying several genes having calA and/or crzA dependent mRNA accumulation patterns. We also demonstrate that contrary to previous findings, the gene encoding the Aspergillus nidulans calcineurin subunit homologue, cnaA, is not essential and that the cnaA deletion mutant shares the morphological phenotypes observed in the corresponding A. fumigatus mutant, DcalA. Exploiting the A. nidulans model system, we have linked calcineurin activity with asexual developmental induction, finding that CrzA supports appropriate developmental induction in a calcineurin and brlA-dependent manner in both species.

Paracoccidioides is a fungal pathogen and the cause of paracoccidioidomycosis, a health-threatening human systemic mycosis endemic to Latin America. Infection by Paracoccidioides, a dimorphic fungus in the order Onygenales, is coupled with a thermally regulated transition from a soil-dwelling filamentous form to a yeast-like pathogenic form. To better understand the genetic basis of growth and pathogenicity in Paracoccidioides, we sequenced the genomes of two strains of Paracoccidioides brasiliensis (Pb03 and Pb18) and one strain of Paracoccidioides lutzii (Pb01). These genomes range in size from 29.1 Mb to 32.9 Mb and encode 7,610 to 8,130 genes. To enable genetic studies, we mapped 94% of the P. brasiliensis Pb18 assembly onto five chromosomes. We characterized gene family content across Onygenales and related fungi, and within Paracoccidioides we found expansions of the fungal-specific kinase family FunK1. Additionally, the Onygenales have lost many genes involved in carbohydrate metabolism and fewer genes involved in protein metabolism, resulting in a higher ratio of proteases to carbohydrate active enzymes in the Onygenales than their relatives. To determine if gene content correlated with growth on different substrates, we screened the non-pathogenic onygenale Uncinocarpus reesii, which has orthologs for 91% of Paracoccidioides metabolic genes, for growth on 190 carbon sources. U. reesii showed growth on a limited range of carbohydrates, primarily basic plant sugars and cell wall components; this suggests that Onygenales, including dimorphic fungi, can degrade cellulosic plant material in the soil. In addition, U. reesii grew on gelatin and a wide range of dipeptides and amino acids, indicating a preference for proteinaceous growth substrates over carbohydrates, which may enable these fungi to also degrade animal biomass. These capabilities for degrading plant and animal substrates suggest a duality in lifestyle that could enable pathogenic species of Onygenales to transfer from soil to animal hosts.

Here, we investigated which stress responses were influenced by the MpkC and SakA mitogen-activated protein kinases of the high-osmolarity glycerol (HOG) pathway in the fungal pathogen Aspergillus fumigatus. The ΔsakA and the double ΔmpkC ΔsakA mutants were more sensitive to osmotic and oxidative stresses, and to cell wall damaging agents. Both MpkC::GFP and SakA::GFP translocated to the nucleus upon osmotic stress and cell wall damage, with SakA::GFP showing a quicker response. The phosphorylation state of MpkA was determined post exposure to high concentrations of congo red and Sorbitol. In the wild-type strain, MpkA phosphorylation levels progressively increased in both treatments. In contrast, the ΔsakA mutant had reduced MpkA phosphorylation, and surprisingly, the double ΔmpkC ΔsakA had no detectable MpkA phosphorylation. A. fumigatus ΔsakA and ΔmpkC were virulent in mouse survival experiments, but they had a 40% reduction in fungal burden. In contrast, the ΔmpkC ΔsakA double mutant showed highly attenuated virulence, with approximately 50% mice surviving and a 75% reduction in fungal burden. We propose that both cell wall integrity (CWI) and HOG pathways collaborate, and that MpkC could act by modulating SakA activity upon exposure to several types of stresses and during CW biosynthesis.

Extracellular vesicles (EVs) has been considered an alternative process for intercellular communication. EVs release by filamentous fungi and the role of vesicular secretion during fungus-host cells interaction remain unknown. Here, we identified the secretion of EVs from the pathogenic filamentous fungus, Aspergillus fumigatus. Analysis of the structure of EVs demonstrated that A. fumigatus produces round shaped bilayer structures ranging from 100 to 200 nm size, containing ergosterol and a myriad of proteins involved in REDOX, cell wall remodeling and metabolic functions of the fungus. We demonstrated that macrophages can phagocytose A. fumigatus EVs. Phagocytic cells, stimulated with EVs, increased fungal clearance after A. fumigatus conidia challenge. EVs were also able to induce the production of TNF-α and CCL2 by macrophages and a synergistic effect was observed in the production of these mediators when the cells were challenged with the conidia. In bone marrow-derived neutrophils (BMDN) treated with EVs, there was enhancement of the production of TNF-α and IL-1β in response to conidia. Together, our results demonstrate, for the first time, that A. fumigatus produces EVs containing a diverse set of proteins involved in fungal physiology and virulence. Moreover, EVs are biologically active and stimulate production of inflammatory mediators and fungal clearance.

We investigated the evolution of resistance to the antifungal drug itraconazole in replicate populations of Aspergillus fumigatus that were founded from a strain with a genotype of sensitivity to a single drug and then propagated under uniform conditions. For each population, conidia were serially transferred 10 times to agar medium either with or without itraconazole. After 10 transfers in medium supplemented with itraconazole, 10 itraconazole-resistant mutant strains were isolated from two populations. These mutant strains had different growth rates and different levels of itraconazole resistance. Analysis of the ergosterol contents of these mutants showed that they accumulate ergosterol when they are grown in the presence of itraconazole. The replacement of the CYP51A gene of the wild-type strain changed the susceptibility pattern of this strain to one of itraconazole resistance only when CYP51A genes with N22D and M220I mutations were used as selectable marker genes. Real-time quantitative reverse transcription-PCR was used to assess the levels of expression of the Afumdr1, Afumdr2, Afumdr3, Afumdr4, AtrF transporter, CYP51A, and CYP51B genes in these mutant strains. Most mutants showed either constitutive high-level expression or induction upon exposure of Afumdr3, Afumdr4, and AtrF to itraconazole. Our results suggest that overexpression of drug efflux pumps and/or selection of drug target site mutations are at least partially responsible for itraconazole resistance and could be considered mechanisms for the emergence of clinical resistance to this drug.

Aspergillus fumigatus is an opportunistic fungal pathogen that causes invasive aspergillosis (IA), a life-threatening disease in immunocompromised humans. The echinocandin caspofungin, adopted as a second-line therapy in combating IA, is a β-1,3-glucan synthase inhibitor, which, when used in high concentrations, reverts the anticipated A. fumigatus growth inhibition, a phenomenon called the “caspofungin paradoxical effect” (CPE). The CPE has been widely associated with increased chitin content in the cell wall due to a compensatory upregulation of chitin synthase-encoding genes. Here, we demonstrate that the CPE is dependent on the cell wall integrity (CWI) mitogen-activated protein kinase MpkAMPK1 and its associated transcription factor (TF) RlmARLM1, which regulate chitin synthase gene expression in response to different concentrations of caspofungin. Furthermore, the calcium- and calcineurin-dependent TF CrzA binds to and regulates the expression of specific chitin synthase genes during the CPE. These results suggest that the regulation of cell wall biosynthetic genes occurs by several cellular signaling pathways. In addition, CrzA is also involved in cell wall organization in the absence of caspofungin. Differences in the CPE were also observed between two A. fumigatus clinical isolates, which led to the identification of a novel basic leucine zipper TF, termed ZipD. This TF functions in the calcium-calcineurin pathway and is involved in the regulation of cell wall biosynthesis genes. This study therefore unraveled additional mechanisms and novel factors governing the CPE response, which ultimately could aid in developing more effective antifungal therapies.