Heterogeneity in the transcriptional response of the human pathogen Aspergillus fumigatus to the antifungal agent caspofungin

An oxylipin signal confers protection against antifungal echinocandins in pathogenic aspergilli

Aspergillus fumigatus is the leading causative agent of life-threatening invasive aspergillosis in immunocompromised individuals. One antifungal class used to treat Aspergillus infections is the fungistatic echinocandins, semisynthetic drugs derived from naturally occurring fungal lipopeptides. By inhibiting beta-1,3-glucan synthesis, echinocandins cause both fungistatic stunting of hyphal growth and repeated fungicidal lysis of apical tip compartments. Here, we uncover an endogenous mechanism of echinocandin tolerance in A. fumigatus whereby the inducible oxylipin signal 5,8-diHODE confers protection against tip lysis via the transcription factor ZfpA. Treatment of A. fumigatus with echinocandins induces 5,8-diHODE synthesis by the fungal oxygenase PpoA in a ZfpA dependent manner resulting in a positive feedback loop. This protective 5,8-diHODE/ZfpA signaling relay is conserved among diverse isolates of A. fumigatus and in two other Aspergillus pathogens. Our findings reveal an oxylipin-directed growth program-possibly arisen through natural encounters with native echinocandin producing fungi-that enables echinocandin tolerance in pathogenic aspergilli.

Strain heterogeneity in a non-pathogenic fungus highlights factors contributing to virulence

Fungal pathogens exhibit extensive strain heterogeneity, including variation in virulence. Whether closely related non-pathogenic species also exhibit strain heterogeneity remains unknown. Here, we comprehensively characterized the pathogenic potentials (i.e., the ability to cause morbidity and mortality) of 16 diverse strains of Aspergillus fischeri, a non-pathogenic close relative of the major pathogen Aspergillus fumigatus. In vitro immune response assays and in vivo virulence assays using a mouse model of pulmonary aspergillosis showed that A. fischeri strains varied widely in their pathogenic potential. Furthermore, pangenome analyses suggest that A. fischeri genomic and phenotypic diversity is even greater. Genomic, transcriptomic, and metabolomic profiling identified several pathways and secondary metabolites associated with variation in virulence. Notably, strain virulence was associated with the simultaneous presence of the secondary metabolites hexadehydroastechrome and gliotoxin. We submit that examining the pathogenic potentials of non-pathogenic close relatives is key for understanding the origins of fungal pathogenicity.

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.

Transcription factors SltA and CrzA reversely regulate calcium homeostasis under calcium-limited conditions

IMPORTANCE

Calcium ions are ubiquitous intracellular signaling molecules for many signaling pathways regulating the fungal response to stress and antifungal drugs. The concentration of intracellular calcium is tightly regulated in its storage, release, and distribution. CrzA is the best-studied transcription factor that regulates this process under sufficient calcium or other external signals. However, CrzA was excluded from nuclei and then lost transcriptional activation under calcium-limited conditions. The regulators in the Ca2+ signaling pathway under calcium-limited conditions remain unclear. Here, we identified SltA as a key regulator in the Ca2+ signaling pathway under calcium-limited conditions, and the underlying mechanisms were further explored in Aspergillus fumigatus. These findings reveal a transcriptional control pathway that precisely regulates calcium homeostasis under calcium-limited conditions.

The Aspergillus fumigatus UPR is variably activated across nutrient and host environments and is critical for the establishment of corneal infection

The Aspergillus fumigatus unfolded protein response (UPR) is a two-component relay consisting of the ER-bound IreA protein, which splices and activates the mRNA of the transcription factor HacA. Spliced hacA accumulates under conditions of acute ER stress in vitro, and UPR null mutants are hypovirulent in a murine model of invasive pulmonary infection. In this report, we demonstrate that a hacA deletion mutant (ΔhacA) is furthermore avirulent in a model of fungal keratitis, a corneal infection, and an important cause of ocular morbidity and unilateral blindness worldwide. Interestingly, we demonstrate that A. fumigatus hacA is spliced in infected lung samples, but not in the cornea, suggesting the amount of ER stress experienced by the fungus varies upon the host niche. To better understand how the UPR contributes to fungal cell biology across a spectrum of ER-stress levels, we employed transcriptomics on the wild-type and ΔhacA strains in glucose minimal media (low stress), glucose minimal media with dithiothreitol (high stress), and gelatin minimal media as a proxy for the nutrient stress encountered in the cornea (mid-level stress). These data altogether reveal a unique HacA-dependent transcriptome under each condition, suggesting that HacA activity is finely-tuned and required for proper fungal adaptation in each environment. Taken together, our results indicate that the fungal UPR could serve as an important antifungal target in the setting of both invasive pulmonary and corneal infections.

Aspergillus fumigatus transcription factor ZfpA regulates hyphal development and alters susceptibility to antifungals and neutrophil killing during infection

Hyphal growth is essential for host colonization during Aspergillus infection. The transcription factor ZfpA regulates A. fumigatus hyphal development including branching, septation, and cell wall composition. However, how ZfpA affects fungal growth and susceptibility to host immunity during infection has not been investigated. Here, we use the larval zebrafish-Aspergillus infection model and primary human neutrophils to probe how ZfpA affects A. fumigatus pathogenesis and response to antifungal drugs in vivo. ZfpA deletion promotes fungal clearance and attenuates virulence in wild-type hosts and this virulence defect is abrogated in neutrophil-deficient zebrafish. ZfpA deletion also increases susceptibility to human neutrophils ex vivo while overexpression impairs fungal killing. Overexpression of ZfpA confers protection against the antifungal caspofungin by increasing chitin synthesis during hyphal development, while ZfpA deletion reduces cell wall chitin and increases caspofungin susceptibility in neutrophil-deficient zebrafish. These findings suggest a protective role for ZfpA activity in resistance to the innate immune response and antifungal treatment during A. fumigatus infection.

Aspergillus fumigatus transcription factor ZfpA regulates hyphal development and alters susceptibility to antifungals and neutrophil killing during infection

Hyphal growth is essential for host colonization during Aspergillus infection. The transcription factor ZfpA regulates A. fumigatus hyphal development including branching, septation, and cell wall composition. However, how ZfpA affects fungal growth and susceptibility to host immunity during infection has not been investigated. Here, we use the larval zebrafish- Aspergillus infection model and primary human neutrophils to probe how ZfpA affects A. fumigatus pathogenesis and response to antifungal drugs in vivo . ZfpA deletion promotes fungal clearance and attenuates virulence in wild-type hosts and this virulence defect is abrogated in neutrophil-deficient zebrafish. ZfpA deletion also increases susceptibility to human neutrophils ex vivo while overexpression impairs fungal killing. Overexpression of ZfpA confers protection against the antifungal caspofungin by increasing chitin synthesis during hyphal development, while ZfpA deletion reduces cell wall chitin and increases caspofungin susceptibility in neutrophil-deficient zebrafish. These findings suggest a protective role for ZfpA activity in resistance to the innate immune response and antifungal treatment during A. fumigatus infection.

Author Summary

Aspergillus fumigatus is a common environmental fungus that can infect immunocompromised people and cause a life-threatening disease called invasive aspergillosis. An important step during infection is the development of A. fumigatus filaments known as hyphae. A. fumigatus uses hyphae to acquire nutrients and invade host tissues, leading to tissue damage and disseminated infection. In this study we report that a regulator of gene transcription in A. fumigatus called ZfpA is important for hyphal growth during infection. We find that ZfpA activity protects the fungus from being killed by innate immune cells and decreases the efficacy of antifungal drugs during infection by regulating construction of the cell wall, an important protective layer for fungal pathogens. Our study introduces ZfpA as an important genetic regulator of stress tolerance during infection that protects A. fumigatus from the host immune response and antifungal drugs.

Chitin Biosynthesis in Aspergillus Species

The fungal cell wall (FCW) is a dynamic structure responsible for the maintenance of cellular homeostasis, and is essential for modulating the interaction of the fungus with its environment. It is composed of proteins, lipids, pigments and polysaccharides, including chitin. Chitin synthesis is catalyzed by chitin synthases (CS), and up to eight CS-encoding genes can be found in Aspergillus species. This review discusses in detail the chitin synthesis and regulation in Aspergillus species, and how manipulation of chitin synthesis pathways can modulate fungal growth, enzyme production, virulence and susceptibility to antifungal agents. More specifically, the metabolic steps involved in chitin biosynthesis are described with an emphasis on how the initiation of chitin biosynthesis remains unknown. A description of the classification, localization and transport of CS was also made. Chitin biosynthesis is shown to underlie a complex regulatory network, with extensive cross-talks existing between the different signaling pathways. Furthermore, pathways and recently identified regulators of chitin biosynthesis during the caspofungin paradoxical effect (CPE) are described. The effect of a chitin on the mammalian immune system is also discussed. Lastly, interference with chitin biosynthesis may also be beneficial for biotechnological applications. Even after more than 30 years of research, chitin biosynthesis remains a topic of current interest in mycology.

AfSec1 is a signal peptidase and removes signal peptides of 1,3-β-glucanosyltransferases in Aspergillus fumigatus

To identify the infection mechanism of Aspergillus fumigatus, which is an opportunistic fungal pathogen, we analyzed the expression profile of the whole genome of A. fumigatus during the infection of murine macrophages. A previously reported RNA-seq data analysis showed that many genes involved in cell wall synthesis were upregulated during the infection process. Interestingly, AfSec1 (3g12840), which encodes a putative signal peptidase, was upregulated dramatically, and its putative target protein Gel1, which encodes a 1,3-β-glucanosyltransferase, was also upregulated. Instead of the AfSec1 deletion strain, the AfSec1-ΔP strain was constructed, in which the promoter region of AfSec1 was deleted, and AfSec1 expression was not detected in the AfSec1-ΔP strain. The expression of AfSec1 was recovered by the introduction of the promoter region (the AfSec1-ΔP/P strain). The nonprocessed form of Gel1 was identified in the AfSec1-ΔP strain, which lacked the promoter, but mature forms of Gel1 were found in the wild-type and in AfSec1-ΔP/P, which was the promoter complementation strain. In the plate assay, the AfSec1-ΔP strain showed higher sensitivity against caspofungin than the wild-type. However, compared with the wild-type, the deletion strain showed no difference in the sensitivity to other antifungal drugs, such as amphotericin B and voriconazole, which inhibit different targets compared with caspofungin. The AfSec1-ΔP strain exhibited ∼20% lower levels of β-glucan in the cell wall than the wild-type. Finally, the virulence decreased when the promoter region of AfSec1 was deleted, as observed in the murine infection test and conidia-killing assay using human macrophages and neutrophils. These results suggest that AfSec1 exerts signal peptidase activity on its target Gel1 and has an important role in fungal pathogenesis.

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.

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.

The Caspofungin Paradoxical Effect is a Tolerant "Eagle Effect" in the Filamentous Fungal Pathogen Aspergillus fumigatus

Cell responses against antifungals other than resistance have rarely been studied in filamentous fungi, while terms such as tolerance and persistence are well-described for bacteria and increasingly examined in yeast-like organisms. Aspergillus fumigatus is a filamentous fungal pathogen that causes a disease named aspergillosis, for which caspofungin (CAS), a fungistatic drug, is used as a second-line therapy. Some A. fumigatus clinical isolates can survive and grow in CAS concentrations above the minimum effective concentration (MEC), a phenomenon known as "caspofungin paradoxical effect" (CPE). Here, we evaluated the CPE in 67 A. fumigatus clinical isolates by calculating recovery rate (RR) values, where isolates with an RR of ≥0.1 were considered CPE⁺ while isolates with an RR of <0.1 were classified as CPE^-. Conidia produced by three CPE⁺ clinical isolates, CEA17 (RR = 0.42), Af293 (0.59), and CM7555 (0.38), all showed the ability to grow in high levels of CAS, while all conidia produced by the CPE^- isolate IFM61407 (RR = 0.00) showed no evidence of paradoxical growth. Given the importance of the calcium/calcineurin/transcription factor-CrzA pathway in CPE regulation, we also demonstrated that all ΔcrzA^CEA17 (CPE⁺) conidia exhibited CPE while 100% of ΔcrzA^Af293 (CPE^-) did not exhibit CPE. Because all spores derived from an individual strain were phenotypically indistinct with respect to CPE, it is likely that CPE is a genetically encoded adaptive trait that should be considered an antifungal-tolerant phenotype. Because the RR parameter showed that the strength of the CPE was not uniform between strains, we propose that the mechanisms which govern this phenomenon are multifactorial. IMPORTANCE The "Eagle effect," initially described for bacterial species, which reflects the capacity of some strains to growth above the minimum inhibitory concentration (MIC) of specific antimicrobial agents, has been known for more than 70 years. However, its underlying mechanism of action in fungi is not fully understood and its connection with other phenomena such as tolerance or persistence is not clear yet. Here, based on the characterization of the "caspofungin paradoxical effect" in several Aspergillus fumigatus clinical isolates, we demonstrate that all conidia from A. fumigatus CPE⁺ strains are able to grow in high levels of the drug while all conidia produced by CPE^- strains show no evidence of paradoxical growth. This work fills a gap in the understanding of this multifactorial phenomenon by proposing that CPE in A. fumigatus should be considered a tolerant but not persistent phenotype.

Regulation of gliotoxin biosynthesis and protection in Aspergillus species

Aspergillus fumigatus causes a range of human and animal diseases collectively known as aspergillosis. A. fumigatus possesses and expresses a range of genetic determinants of virulence, which facilitate colonisation and disease progression, including the secretion of mycotoxins. Gliotoxin (GT) is the best studied A. fumigatus mycotoxin with a wide range of known toxic effects that impair human immune cell function. GT is also highly toxic to A. fumigatus and this fungus has evolved self-protection mechanisms that include (i) the GT efflux pump GliA, (ii) the GT neutralising enzyme GliT, and (iii) the negative regulation of GT biosynthesis by the bis-thiomethyltransferase GtmA. The transcription factor (TF) RglT is the main regulator of GliT and this GT protection mechanism also occurs in the non-GT producing fungus A. nidulans. However, the A. nidulans genome does not encode GtmA and GliA. This work aimed at analysing the transcriptional response to exogenous GT in A. fumigatus and A. nidulans, two distantly related Aspergillus species, and to identify additional components required for GT protection. RNA-sequencing shows a highly different transcriptional response to exogenous GT with the RglT-dependent regulon also significantly differing between A. fumigatus and A. nidulans. However, we were able to observe homologs whose expression pattern was similar in both species (43 RglT-independent and 11 RglT-dependent). Based on this approach, we identified a novel RglT-dependent methyltranferase, MtrA, involved in GT protection. Taking into consideration the occurrence of RglT-independent modulated genes, we screened an A. fumigatus deletion library of 484 transcription factors (TFs) for sensitivity to GT and identified 15 TFs important for GT self-protection. Of these, the TF KojR, which is essential for kojic acid biosynthesis in Aspergillus oryzae, was also essential for virulence and GT biosynthesis in A. fumigatus, and for GT protection in A. fumigatus, A. nidulans, and A. oryzae. KojR regulates rglT, gliT, gliJ expression and sulfur metabolism in Aspergillus species. Together, this study identified conserved components required for GT protection in Aspergillus species.

A. fumigatus secretes mycotoxins that are essential for its virulence and pathogenicity. Gliotoxin (GT) is a sulfur-containing mycotoxin, which is known to impair several aspects of the human immune response. GT is also toxic to different fungal species, which have evolved several GT protection strategies. To further decipher these responses, we used transcriptional profiling aiming to compare the response to GT in the GT producer A. fumigatus and the GT non-producer A. nidulans. This analysis allowed us to identify additional genes with a potential role in GT protection. We also identified 15 transcription factors (TFs) encoded in the A. fumigatus genome that are important for conferring resistance to exogenous gliotoxin. One of these TFs, KojR, which is essential for A. oryzae kojic acid production, is also important for virulence in A. fumigatus and GT protection in A. fumigatus, A. nidulans and A. oryzae. KojR regulates the expression of genes important for gliotoxin biosynthesis and protection, and sulfur metabolism. Together, this work identified conserved components required for gliotoxin protection in Aspergillus species.

Chromatin profiling reveals heterogeneity in clinical isolates of the human pathogen Aspergillus fumigatus

Invasive Pulmonary Aspergillosis, which is caused by the filamentous fungus Aspergillus fumigatus, is a life-threatening infection for immunosuppressed patients. Chromatin structure regulation is important for genome stability maintenance and has the potential to drive genome rearrangements and affect virulence and pathogenesis of pathogens. Here, we performed the first A. fumigatus global chromatin profiling of two histone modifications, H3K4me3 and H3K9me3, focusing on the two most investigated A. fumigatus clinical isolates, Af293 and CEA17. In eukaryotes, H3K4me3 is associated with active transcription, while H3K9me3 often marks silent genes, DNA repeats, and transposons. We found that H3K4me3 deposition is similar between the two isolates, while H3K9me3 is more variable and does not always represent transcriptional silencing. Our work uncovered striking differences in the number, locations, and expression of transposable elements between Af293 and CEA17, and the differences are correlated with H3K9me3 modifications and higher genomic variations among strains of Af293 background. Moreover, we further showed that the Af293 strains from different laboratories actually differ in their genome contents and found a frequently lost region in chromosome VIII. For one such Af293 variant, we identified the chromosomal changes and demonstrated their impacts on its secondary metabolites production, growth and virulence. Overall, our findings not only emphasize the influence of genome heterogeneity on A. fumigatus fitness, but also caution about unnoticed chromosomal variations among common laboratory strains.