Genome Sequences of Two Strains of the Food Spoilage Mold Aspergillus fischeri

Genomic and Molecular Identification of Genes Contributing to the Caspofungin Paradoxical Effect in Aspergillus fumigatus

Aspergillus fumigatus is a deadly opportunistic fungal pathogen responsible for ~100,000 annual deaths. Azoles are the first line antifungal agent used against A. fumigatus, but azole resistance has rapidly evolved making treatment challenging. Caspofungin is an important second-line therapy against invasive pulmonary aspergillosis, a severe A. fumigatus infection. Caspofungin functions by inhibiting β-1,3-glucan synthesis, a primary and essential component of the fungal cell wall. A phenomenon termed the caspofungin paradoxical effect (CPE) has been observed in several fungal species where at higher concentrations of caspofungin, chitin replaces β-1,3-glucan, morphology returns to normal, and growth rate increases. CPE appears to occur in vivo, and it is therefore clinically important to better understand the genetic contributors to CPE. We applied genomewide association (GWA) analysis and molecular genetics to identify and validate candidate genes involved in CPE. We quantified CPE across 67 clinical isolates and conducted three independent GWA analyses to identify genetic variants associated with CPE. We identified 48 single nucleotide polymorphisms (SNPs) associated with CPE. We used a CRISPR/Cas9 approach to generate gene deletion mutants for seven genes harboring candidate SNPs. Two null mutants, ΔAfu3g13230 and ΔAfu4g07080 (dscP), resulted in reduced basal growth rate and a loss of CPE. We further characterized the dscP phosphatase-null mutant and observed a significant reduction in conidia production and extremely high sensitivity to caspofungin at both low and high concentrations. Collectively, our work reveals the contribution of Afu3g13230 and dscP in CPE and sheds new light on the complex genetic interactions governing this phenotype. IMPORTANCE This is one of the first studies to apply genomewide association (GWA) analysis to identify genes involved in an Aspergillus fumigatus phenotype. A. fumigatus is an opportunistic fungal pathogen that causes hundreds of thousands of infections and ~100,000 deaths each year, and antifungal resistance has rapidly evolved in this species. A phenomenon called the caspofungin paradoxical effect (CPE) occurs in some isolates, where high concentrations of the drug lead to increased growth rate. There is clinical relevance in understanding the genetic basis of this phenotype, since caspofungin concentrations could lead to unintended adverse clinical outcomes in certain cases. Using GWA analysis, we identified several interesting candidate polymorphisms and genes and then generated gene deletion mutants to determine whether these genes were important for CPE. Two of these mutant strains (ΔAfu3g13230 and ΔAfu4g07080/ΔdscP) displayed a loss of the CPE. This study sheds light on the genes involved in clinically important phenotype CPE.

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.

Comparative Genomics Reveals a Single Nucleotide Deletion in pksP That Results in White-Spore Phenotype in Natural Variants of Aspergillus fumigatus

Aspergillus fumigatus is a potentially deadly opportunistic human pathogen. A. fumigatus has evolved a variety of mechanisms to evade detection by the immune system. For example, the conidium surface is covered in a layer of 1,8-dihydroxynaphthalene (DHN) melanin which masks the antigen macrophages use for recognition. DHN melanin also protects conidia from ultraviolet radiation and gives A. fumigatus conidia their characteristic green-grayish color. Here, we conducted genomic analysis of two closely related white-spore natural variants of A. fumigatus in comparison to two closely related green-spore isolates to identify a genetic basis of the white-spore phenotype. Illumina whole-genome resequencing data of the four isolates was used to identify variants that were shared in the white-spore isolates and different from both the green-spore isolates and the Af293 reference genome (which is also a green-spore isolate). We identified 4,279 single nucleotide variants and 1,785 insertion/deletions fitting this pattern. Among these, we identified 64 variants predicted to be high impact, loss-of-function mutations. One of these variants is a single nucleotide deletion that results in a frameshift in pksP (Afu2g17600), the core biosynthetic gene in the DHN melanin encoding gene cluster. The frameshift mutation in the white-spore isolates leads to a truncated protein in which a phosphopantetheine attachment site (PP-binding domain) is interrupted and an additional PP-binding domain and a thioesterase domain are omitted. Growth rate analysis of white-spore and green-spore isolates at 37°C and 48°C revealed that white-spore isolates are thermosensitive. Growth rate of A. fumigatus Af293 and a pksP null mutant in the Af293 background suggests pksP is not directly involved in the thermosensitivity phenotype. Further, our study identified a mutation in a gene (Afu4g04740) associated with thermal sensitivity in yeasts which could also be responsible for the thermosensitivity of the white-spore mutants. Overall, we used comparative genomics to identify the mutation and protein alterations responsible for the white-spore phenotype of environmental isolates of A. fumigatus.

Genome Sequence of Aspergillus aculeatinus IC_8, Isolated from an Indoor Air Sample of an Urban Housing Complex in Abidjan, Ivory Coast

Aspergillus aculeatinus is an industrially important species of Aspergillus section Nigri capable of producing bioactive, antibiotic, and antitumor compounds. We sequenced the genome of a strain of A. aculeatinus that was isolated from the interior of a housing complex in Abidjan, Ivory Coast.

Genome-Wide Association for Itraconazole Sensitivity in Non-resistant Clinical Isolates of Aspergillus fumigatus

Aspergillus fumigatus is a potentially lethal opportunistic pathogen that infects over ~200,000 people and causes ~100,000 deaths per year globally. Treating A. fumigatus infections is particularly challenging because of the recent emergence of azole-resistance. The majority of studies focusing on the molecular mechanisms underlying azole resistance have examined azole-resistant isolates. However, isolates that are susceptible to azoles also display variation in their sensitivity, presenting a unique opportunity to identify genes contributing to azole sensitivity. Here, we used genome-wide association (GWA) analysis to identify loci involved in azole sensitivity by analyzing the association between 68,853 SNPs and itraconazole (ITCZ) minimum inhibitory concentration (MIC) in 76 clinical isolates of A. fumigatus from Japan. Population structure analysis suggests the presence of four distinct populations, with ITCZ MICs distributed relatively evenly across populations. We independently conducted GWA when treating ITCZ MIC as a quantitative trait and a binary trait, and identified two SNPs with strong associations in both analyses. These SNPs fell within the coding regions of Afu2g02220 and Afu2g02140. We functionally validated Afu2g02220 by knocking it out using a CRISPR/Cas9 approach, because orthologs of this gene are involved in sterol modification and ITCZ targets the ergosterol biosynthesis pathway. Knockout strains displayed no difference in growth compared to the parent strain in minimal media, yet a minor but consistent inhibition of growth in the presence of 0.15 μg/ml ITCZ. Our results suggest that GWA paired with efficient gene deletion is a powerful and unbiased strategy for identifying the genetic basis of complex traits in A. fumigatus.

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.