Next Article in Journal
Genome-Scale Metabolic Model of Actinosynnema pretiosum ATCC 31280 and Its Application for Ansamitocin P-3 Production Improvement
Next Article in Special Issue
Golgi Reassembly and Stacking Protein (GRASP) Participates in Vesicle-Mediated RNA Export in Cryptococcus neoformans
Previous Article in Journal / Special Issue
Dual RNA-Seq Analysis of Trichophyton rubrum and HaCat Keratinocyte Co-Culture Highlights Important Genes for Fungal-Host Interaction
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 9
Issue 7
10.3390/genes9070363
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Comparative Analysis of Aspergillus fumigatus Strains: The Reference Genome as a Matter of Concern

by
Rocio Garcia-Rubio
1,†,
Sara Monzon
2,†,
Laura Alcazar-Fuoli
1,
Isabel Cuesta
2 and
Emilia Mellado
1,*
1
Mycology Reference Laboratory, National Centre for Microbiology, Instituto de Salud Carlos III (ISCIII), Majadahonda, 28220 Madrid, Spain
2
Bioinformatics Unit, Common Scientific Technical Units, ISCIII, Majadahonda, 28220 Madrid, Spain
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Genes 2018, 9(7), 363; https://0-doi-org.brum.beds.ac.uk/10.3390/genes9070363
Submission received: 28 June 2018 / Revised: 13 July 2018 / Accepted: 16 July 2018 / Published: 19 July 2018
(This article belongs to the Special Issue Fungal Pathogenesis in Humans: The Growing Threat)
Browse Figures
Versions Notes

Abstract

:
Aspergillus fumigatus is a ubiquitous saprophytic mold and a major pathogen in immunocompromised patients. The effectiveness of triazole compounds, the A. fumigatus first line treatment, is being threatened by a rapid and global emergence of azole resistance. Whole genome sequencing (WGS) has emerged as an invaluable tool for the analysis of genetic differences between A. fumigatus strains, their genetic background, and antifungal resistance development. Although WGS analyses can provide a valuable amount of novel information, there are some limitations that should be considered. These analyses, based on genome-wide comparative data and single nucleotide variant (SNV) calling, are dependent on the quality of sequencing, assembling, the variant calling criteria, as well as on the suitable selection of the reference genome, which must be genetically close to the genomes included in the analysis. In this study, 28 A. fumigatus genomes sequenced in-house and 73 available in public data bases have been analyzed. All genomes were distributed in four clusters and showed a variable number of SNVs depending on the genome used as reference (Af293 or A1163). Each reference genome belonged to a different cluster. The results highlighted the importance of choosing the most suitable A. fumigatus reference genome to avoid misleading conclusions.

1. Introduction

Aspergillus fumigatus is a saprophytic filamentous fungus and the principal causative agent of human aspergillosis [1,2]. There are many types of diseases caused by A. fumigatus and their symptoms vary according to the site of infection and host health condition. Predominantly, this species causes invasive infections, such as invasive pulmonary aspergillosis (IPA), with high mortality rates in immunocompromised or immunosuppressed patients [3,4,5]. Currently, the treatment options are limited to three classes of antifungal drugs (azoles, echinocandins and polyenes) and azoles are drugs of first line for treating Aspergillus-caused diseases. However, azole drug efficacy is threatened by the worldwide spreading of resistance [6].
Despite all the technological improvements developed in recent years, many basic aspects about the biology of this opportunistic pathogen remained largely unknown. Genomics and whole genome sequencing (WGS) have emerged as useful tools to greatly enhance knowledge and understanding of infectious diseases and clinical microbiology [7,8,9,10]. In this context, genome-wide sequencing using high throughput sequencing (HTS) together with alignment comparison analysis has been described as a good approach to identify single nucleotide variants (SNVs) or polymorphisms (SNPs) for the analysis of genetic differences between A. fumigatus strains [11]. The genetic diversity between isolates can be used to explain many different aspects such as the strain background lineage, strain specific virulence phenotypes, the potential factors involved in antifungal resistant development, and more commonly, determining the genetic strain relatedness in epidemiological surveillance and outbreak studies [12,13,14].
Although WGS studies provide a valuable amount of novel data, there are some considerations that should be taken into account. To date, most WGS analyses are based on identifying SNVs in a set of strains in comparison to a reference genome, and consequently the results will be completely dependent on the chosen reference. Therefore, the identified SNVs will have a different significance depending on how genetically close the reference genome will be from the strains included in the analysis. Classically, most advances have been focused on unravelling the genomics of two different A. fumigatus reference strains [15], Af293 (AAHF00000000.1) and A1163 (ABDB00000000.1), or their derivatives. First, the Af293 strain was fully sequenced in 2005 [16], and soon after, the A1163 strain was the second A. fumigatus genome sequenced by the J. Craig Venter Institute [17]. This strain was a derivative of another clinical isolate, the CBS144-89/CEA10. The first genome comparison between both strains showed that despite high synteny and identity in most of the regions of their genomes, there were some hundreds of genes unique from each strain and not present in the other [18]. This heterogeneity could potentially lead to confusion, so care must be taken in drawing conclusions on different aspects from a single genetic background [15]. A recent study carried out in our laboratory [19] showed that the Af293 reference genome belonged to a specific lineage of A. fumigatus strains which harbors five polymorphisms (F46Y, M172V, N248T, D255E, E427K) in cyp51A gene which codes for the azole target. Based on these phylogenetic studies from whole genome sequencing, the other classic reference genome, the A1163 strain, belonged to a different and distant phylogenetic lineage, which implies that both genomes are significantly different. Also, and apart from the relevance of choosing the most suitable reference genome, the quality of sequencing and assembling will be of great importance for getting reliable results.
In this work, we describe the in-house genome sequencing and comparative analysis of 28 A. fumigatus genomes, including the A. fumigatus genomes of six strains routinely used in Aspergillus basic research laboratories. All these genomes were distributed in four separate clusters, corresponding to very different lineages. In addition, 73 A. fumigatus genomes, including azole susceptible and resistant strains from a worldwide distribution, were included in the study. The whole sequence analysis of these 101 A. fumigatus genomes provided information that led to an increase in the knowledge of A. fumigatus genetic background and the development of azole resistance.

2. Materials and Methods

2.1. Aspergillus fumigatus Isolates

A total of 169 genomes were evaluated including in-house sequenced Aspergillus fumigatus isolates and genomes downloaded from NCBI Sequence Read Archive (SRA) public database (www.ncbi.nlm.nih.gov/sra). The genomic samples were selected following this criteria: (i) all genomes had to be sequenced with Illumina platform; and (ii) had to be paired-end. After quality control filtering (>85% mapping rate against Af293 reference genome and >95% of genome coverage at more than 10× depth of coverage), a total of 101 A. fumigatus strain genome sequences representing 8 different geographical locations were selected for further analysis. Among them, 28 A. fumigatus strains were from the Mycology Reference Laboratory Collection (Table 1). Most of them were clinical azole susceptible strains. Some of them had specific cyp51A modifications: 12 strains harbored three changes (F46Y, M172V, and E427K) and 3 strains had five (F46Y, M172V, N248T, D255E and E427K) [19]. Af293 was included in the latter group of five cyp51A modifications and sequenced again in-house (named as AF293). Another group of A. fumigatus azole susceptible strains were isolates that have been classically used as reference or control strains in most Aspergillus laboratories. This set of strains included the A1163 (also known as CEA10, CBS144.89, FGSC A1163 or AF10), AF237 (called CM237), ATCC204305, akuBKU80 (FGSC A1160) [20], and ATCC46645. Two azole resistant strains without modifications in cyp51A were also included and sequenced in this work (CM7510 and CM7555).
The remaining 73 strain genomes were obtained from public data bases and belonged to strains from the United Kingdom, the Netherlands, Portugal, Denmark, India, Japan and Canada [13,21,22] (NCBI SRA public database) (Table 2). Some of the isolates from The Netherlands and India came from clinical samples (n = 8) and others from environmental soil sources (n = 8). Strains from Japan, UK and Canada were all from clinical origin while the Portuguese strains had an unknown origin (Table 2).

2.2. DNA Extraction and Aspergillus fumigatus Identification

First, conidia from each selected strain were cultured in 3 mL of glucose-yeast extract-peptone (GYEP) broth (2% glucose, 0.3% yeast extract, 1% peptone; OXOID LTD, Basingstoke, Hampshire, England) and grown overnight at 37 °C, after which mycelium mats were harvested and DNA was extracted [23]. The A. fumigatus Spanish isolates included in this work were identified to the species level on the basis of PCR amplification and sequencing of the internal transcribed spacer (ITS) region and the partial amplification of the β-tubulin gene [24].

2.3. DNA Quality and Quantity Assessment

The DNA quality and quantity was assessed using a spectrophotometer (Thermo Scientific NanoDrop One Spectrophotometer, Waltham, MA, USA). Since genomic DNA (gDNA) must be clean and pure for preparing WGS libraries, when the 260/280 or 260/230 rates were lower than 1.8 or higher than 2.2, the DNA was rejected. The strain was extracted again from a new mycelium mat and its quality was assessed. Genomic DNA (gDNAs) were stored at −20 °C until further use.

2.4. Screening of cyp51A Changes: PCR Amplification and Sequence Analysis

The full coding sequences of cyp51A including its promoter, were amplified and sequenced using the PCR conditions described before [25]. To exclude the possibility that any change identified in the sequences was due to PCR-induced errors, each isolate was independently analyzed twice. Isolates were screened for the presence of tandem repeat insertions in the cyp51A promoter region, as well as for the presence of other cyp51A modifications. All the nucleotide sequences were analyzed using the DNASTAR Lasergene package (DNASTAR Inc., Madison, WI, USA).

2.5. DNA Library Preparation and Illumina Whole-Genome Sequencing

Genomic DNA was extracted as previously described. Genomic DNA were quantified afterwards using the QuantiFluor ® dsDNA System and the QuantiFluor ® ST Fluorometer (Promega, Madison, WI, USA) and their quality was determined with the Agilent 2100 Bioanalyzer (Agilent Technologies, Inc., Santa Clara, CA, USA). The preparation of fragmented gDNA libraries was performed using Nextera ® XT Library Prep Kit (Illumina Inc., San Diego, CA, USA), according to the manufacturer’s protocols. The mean fragment length of the libraries ranged from 800 to 1800 bp. Sequencing was conducted in paired-end 2 × 150 bp on a NextSeq 500 system, according to the manufacturer’s protocols (Illumina Inc., San Diego, CA, USA).

2.6. Whole-Genome Sequencing Alignment

The Illumina reads were trimmed using Trimmomatic (version 0.32) [26]. The sequencing adapters and sequences with low quality scores on 3′ ends (Phred score [Q], <20) were trimmed. Raw Illumina WGS reads were quality checked performing a quality control with FastQC (version 0.11.3; Babraham Institute). Data sets were analyzed against two different A. fumigatus reference genomes, the Af293 (GenBank accession number AAHF00000000.1) and the A1163 (GenBank accession number ABDB00000000.1) using WGS-outbreaker v1.0 (Instituto de Salud Carlos III, Madrid, Spain) (https://github.com/BU-ISCIII/WGS-Outbreaker) with default parameters. The pipeline comprised all steps needed for SNV analysis using whole genome sequencing data. Mapping against genome reference was performed with bwa mem (version 0.7.12-r1039) [27], duplicated reads removed using Picard (version 1.140) (http://broadinstitute.github.io/picard), and the bedtools coverage v2.26 program [28] was used to perform further quality controls. Hereafter, in order to identify genetic variations among strains, single nucleotide variant (SNV) detection (variant calling) and SNV matrix generation were performed using GATK version 3.8.0 [29] with best practices parameters. ENSEMBL variant effect predictor script (version 88) was used for variant annotation. The whole genome sequencing project has been deposited in NCBI SRA (project accession number SRP151231).

2.7. Phylogenetic Analysis and Single Nucleotide Variant Comparison

Final step of WGS-Outbreaker pipeline comprised Maximum-likelihood trees construction using RaxML software (version 8.2.9) [30] with GTRCAT model and 100 bootstrap replicates. Phylogenetic trees were visualized and annotated using ggtree R package [31]. SNV comparisons were performed using a custom R script, mapping all genomes to Af293 and A1163 reference genomes.

2.8. Analysis of Genetic Diversity

Genetic diversity was calculated using the Molecular Evolutionary Genetics Analysis software (MEGA, version 7) [32]. The A. fumigatus population was analyzed according to the formed groups (clusters and subclusters) from the previous phylogenic studies. In MEGA, evolutionary distances between sequences can be estimated by computing the number or proportion of nucleotide differences between sequences using the FASTA alignment against each reference genome [33]. Two different models have been used in this work, the p-distance model and another based on the number of differences; (i) the p-distance model calculates the distance as the proportion of different nucleotide sites compared. It is obtained by dividing the number of nucleotide differences (transitions and transversions) by the total number of nucleotides compared. As recommended [34], when the genetic distance was estimated, the complete-deletion option was used, since it normalizes the number of differences based on the number of valid sites compared, not taking into account the alignment gaps and missing data that sequences could contain. However, this model does not make any correction for multiple substitutions at the same site, substitution rate biases (for example, differences in the transitional and transversional rates), or differences in evolutionary rates among sites; (ii) The second model, based on the number of differences, estimates the genetic distances according to the number of sites at which the compared sequences differ (transitions and transversions). In this model, the complete-deletion option was also used.

2.9. Determining Modifications in Genes of Interest

In order to see if the population structure could be based on particular genomic modifications, some genes that have already been described as important in A. fumigatus biology were analyzed in depth in each of the A. fumigatus population clusters formed from the SNV comparisons. One of the genes was the cyp51A gene including its promoter (AFUA_4G06890 in Af293, and AFUB_063960 in A1163). This region was used to determine the azole resistant phenotype based on the resistance mechanism. Since both reference genomes, Af293 and A1163, have opposite mating types, the mating locus genes were included in the analysis: Mat1.1.1 (AFUB_042900), Mat1.2.1 (AFUA_3G06170), and Mat1.2.4 (AFUA_3G06160). The presence or absence of these loci was also tested in all genomes using srst2 software (version 0.1.8) [35]. Every gene of interest was analyzed independently in all genomes.

2.10. Visualization of Depth of Coverage

An exploratory approach using Circos (v 0.69.3) [36], a tool to represent visual data, was used to plot the depth of coverage of 20 whole-sequenced genomes selected to represent each cluster and subcluster in comparison to the Af293 reference genome. Mean coverage data in 10,000 length bins was calculated from bam files with bedtools (version 2.2.17) (http://bedtools.readthedocs.io/en/latest/).

2.11. Mating Type Analysis

A mapping approach was used in order to identify genes associated with each mating type. Since Af293 was Mat1.2 and A1163 was Mat1.1, all samples were mapped against both references using the workflow displayed in Figure S1: (i) samples that had been determined with srst2 as Mat1.2 and mapped against A1163 were selected, and unmapped reads were retained in this step; (ii) the unmapped reads were then mapped against Af293. Statistics of coverage and depth of coverage were calculated for all genes annotated in A1163 reference genome using the bam generated in this step; (iii) genes with at least 70% of coverage were selected; (iv) to discard genes that may be found due to the lack of presence in the A1163 assembly, the same procedure was performed for Mat1.1 samples; (v) finally, genes in common among all Mat1.2 samples minus one but not present in any Mat1.1 sample were exported as a result (Figure S1A). The reverse procedure was performed starting with Mat1.1 samples determining the genes only present in Mat1.1 samples (Figure S1B).

3. Results

3.1. Species Identification

The isolates were identified at the species level and also the azole resistance mechanism was analyzed by PCR amplification and sequencing of the cyp51A gene and its promoter as explained in Materials and Methods section. All the strains were identified as A. fumigatus sensu stricto. The azole susceptibility profile of most of the strains sequenced in the laboratory is described elsewhere [19]. For the remaining strains, the azole susceptibility profile was determined based on the known cyp51A modifications detected by WGS analysis. According to the azole resistance mechanism, the A. fumigatus strains included in this study harbored the following cyp51A modifications: (i) strains with cyp51A single point mutations (G54E/R, M220V/T/K, and G138C) or double point mutations (G448S together with H147Y) in the minority of cases (only one strain); or (ii) strains with tandem integrations and cyp51A modifications (TR34/L98H, TR34/L98H/S297T/F495I). Two azole resistant cyp51A wild type (WT) strains were also included. All these strain features (cyp51A modifications and azole susceptibility profile, as well as the geographical origin and mating type), were described in Table 1 and Table 2. Afterwards, the DNAs were further analyzed using WGS.

3.2. Whole-Genome Sequencing Analysis

All genome sequences mapped to a range of 85–99% genome coverage with at least 10x of depth of coverage and a mapping rate > 94% against Af293 reference genome. Similarly, sequences were mapped to a range of 84–99% genome coverage with at least 10× of depth of coverage against A1163 reference with a mapping rate > 88% (Table S1). The Af293 and the A1163 strains were resequenced in our laboratory (and named hereafter AF293 and CEA10) and had a total of 313 and 1009 SNVs, respectively, compared to their respective reference genomes.

3.3. Phylogenetic Analysis

Phylogenetic analysis using SNV data of Spanish genomes (Table 1) clearly showed that our strains were divided into four clusters independently of the genome used as reference (Figure 1, with Af293 as reference genome and Figure 2, with A1163 as reference genome). In this previous analysis, the different wild type reference strains used, akuBKU80 and CEA10 as derivative strains of A1163, AF293, CM237, ATCC204305 and ATCC46645 were located in separated clusters (Figure 1 and Figure 2). Those clusters were further confirmed when the remaining 73 A. fumigatus genomes available from other countries were included.
Therefore, based on phylogenetic analysis, the complete A. fumigatus population was divided into four well defined clusters (I, II, III, and IV). Two of them could be further subdivided into subclusters: cluster I split in 3 subclusters (I.1, I.2, I.3), and cluster II in 2 (II.1, II.2). Clusters III and IV remained undivided, because any subcluster was identified. Each subcluster was considered as an independent population group, and was used to compare itself against each reference genome (Table 3). The two reference genomes used were distributed in different clusters when compared against the other: A1163 in cluster I and Af293 in cluster III. The number of clusters and the tree topology was maintained regardless of the reference genome chosen.
Based on cyp51A modifications, azole susceptible cyp51A-WT strains together with azole resistant cyp51A single point mutation strains were grouped together in cluster I. In cluster II, there were azole susceptible and resistant strains, with both cyp51A single point mutations and TR34/L98H mechanisms. In cluster III, a set of strains with five cyp51A modifications (F46Y, M172V, N248T, D255E, E427K) grouped together. While in cluster IV strains with three cyp51A modifications (F46Y, M172V, E427K) also grouped together [19]. The Netherlands, Portugal and Canada genomes were distributed in two clusters. The United Kingdom genomes were distributed in three, while Japanese genomes (cyp51A-WT) were only in cluster I and Indian genomes (cyp51A-TR34) in cluster II. Except for the tightly clustered isolates from India in cluster II.2, there was no relationship between the geographical origin of the isolates and the cluster where they grouped.

3.4. Single Nucleotide Variant Analysis against Both References

To explore the genomic differences between the clusters formed from phylogenetic studies, total SNVs between populations were determined. There were a total of 93,609 and 71,844 SNV positions identified against the reference genomes, Af293 and A1163, respectively. No noticeable differences were found in the SNV ranges when all strains were mapped against Af293 (313–165,788) or A1163 (1043–163,367). However, the number of specific SNVs identified in each strain was notably different depending on the genome that was used as reference (Table 4). In fact, some particular SNVs previously detected in our laboratory by PCR (L98H, F46Y, M172V, N248T, D255E, E427K) were specifically searched for in both reference whole genomes. The L98H variant was found in all genomes that harbored this mutation since none of the reference genomes had this change in its own genome. However, the remaining 5 SNVs (F46Y, M172V, N248T, D255E, E427K) were only detected when the reference genome chosen was the A1163 genome, as Af293 specifically harbored these 5 same polymorphisms.
Although the tree topology was maintained regardless of the reference genome chosen, the number of SNVs within clusters changed significantly when comparing each reference against each cluster, showing that both reference genomes were genetically very different (Table 3). The number of SNVs of clusters II and IV were similar against both references, while remarkable differences in SNVs were present in clusters I (closer to A1163) and III (closer to Af293) (Table 3). Considering subclusters, the number of SNVs in subclusters I.1, I.2 and I.3 were almost half when the reference genome was the A1163, than the SNVs obtained using Af293 as reference genome. Another remarkable result was the small number of SNVs found in cluster III when compared to Af293 genome. The SNV differences found against both references were slighter but also considerable in cluster II, higher when compared to A1163. Cluster IV had the greatest number of SNVs, independent of the reference genome used.
Moreover, in order to explore any genomic differences between clusters and subclusters, some features of the SNVs found in each population were determined (Table 5). This table shows that all populations had a similar amount of SNVs, excluding cluster IV, which had much more variants, specifically missense and synonymous, compared to the other clusters.

3.5. Genetic Diversity

Evolutionary distances were estimated by computing the mean proportion or number of nucleotide differences between each subcluster and each reference genome using two different models provided by MEGA: the p-model and the number of differences model.
The pairwise distance matrices referred to the distance between each genome and the others in pairs and was calculated and expressed in proportion and number of differences (Tables S2–S5). The matrices obtained from the mean distance between groups in which each subcluster was compared to each other, were indicated also as the proportion and number of differences of all these comparisons (Tables S6 and S7). Cluster IV was identified as the most distant cluster, according to phylogenetic results, with 123,259, 118,271 and 114,815 differences in SNVs against cluster I, II and III, respectively. Moreover, while cluster III was more similar to cluster II with 36,031 differences, cluster I was closer to cluster II than to III with 38,197 and 50,040 differences, respectively. These data were calculated using Af293 genome as reference. Data using A1163 were not shown since no significant differences were found.
Homogeneity within a group or cluster could be inferred through the analysis of the mean distance within groups (Tables S8 and S9). Cluster IV was the most homogeneous since the number of differences was lower than in the remaining groups, followed by subcluster II.1, independently of the genome that was chosen as reference.

3.6. Genome-Wide Visualization of Depth of Coverage

The depth of coverage from a selected group of whole-sequenced genomes was plotted in a Circos image to make an exploratory analysis, and to visualize similarities and differences in genome structure (Figure 3). None of the 20 A. fumigatus genomes included in this study displayed chromosomal or copy number variations. However, small deletions were observed in multiple chromosomes: regions smaller than 100 kbp seemed to be deleted in chromosomes 3, 4, 5, 6, 7 and 8. Also, what seemed to be large-scale deletions (bigger than 300 kbp) were detected in the edges of chromosomes 1 and 7 in genomes of cluster I and II. It is noteworthy that there was a huge red region in chromosome 4 present in all genomes marked as an expected deletion, even in the two different Af293 genomes downloaded and sequenced in-house although the reference genome used was the same Af293 strain genome.

3.7. Genome Comparisons Based on Their Mating Type

There were no differences in the strain genetic background between genomes harboring the mating type 1 or 2. In fact, strain genomes of both mating types were present in all clusters. In order to detect if there was any differential presence or absence of certain genes based on their particular mating type, all genomes were divided according to their mating type and mapped to both reference genomes (Af293, Mat1.2 and A1163, Mat1.1). No differences in the presence of genes were found in Mat1.1 strain genomes other than those encoded by the Mat1.2 locus itself (Mat1.2.1 and Mat1.2.4). The same situation was found within strains with Mat1.2 genomes; there were no differences in genes other than the gene Mat1.1.1 itself encoded by the Mat1.1 locus.

4. Discussion

Different typing techniques have demonstrated remarkable genomic diversity among strains of A. fumigatus [13,37]. Aspergillus fumigatus is found in multiple environmental niches which contribute to its genotypic and phenotypic heterogeneity among isolates [38]. Therefore, understanding the association between genomic diversity and different biological and phenotypic features, such as the spectrum of pathogenicity and antifungal resistance, among others, is a challenge for multiple Aspergillus laboratories. In this context, WGS is emerging as a promising tool for strain characterization [18,39]. Since WGS considers the entire genetic material of every isolate, it provides the highest discriminatory power that any technique can reach, allowing for the differentiation of closely related strains [39]. In addition, WGS could emerge as a useful tool for gene detection and other sequence-based investigations [9]. This source of an unprecedented amount of novel and comparative data could aid in the development of a huge database which would contain information about the biological differences between A. fumigatus strains, such as their origin, diversity and population structure. Previous WGS studies with a suitable depth of coverage had shown variation in genome structure with small-scale insertions and deletions, and recombination events among A. fumigatus isolates [40]. However, due to the complexity of biological systems, some difficulties need to be overcome in order to correlate complex biological traits with genomic features.
To date, most of the WGS studies are based on genome-scale comparisons in which the number of SNVs of each strain are obtained comparing its genome to a previously selected reference genome. In this work, we have made all comparison independently using two A. fumigatus reference genomes, Af293 and A1163. The differential grouping of the A. fumigatus population in four clusters is driven by the different SNVs found. In fact, the number of SNVs in each cluster and subcluster differed depending on the reference genome used to determine the number of variants. Particularly, the small number of SNVs found in cluster III when it was compared to Af293 reference genome is due to the fact that the genomes included in this cluster were more similar to Af293 and had twice the SNVs obtained with A1163 as a reference. Because cluster III was the smallest, formed by only five genomes, the WGS analysis would be improved by the addition of more genomes when available. In addition, it is important to highlight that the specific five SNVs in cyp51A present in all genomes belonging to cluster III (F46Y, M172V, N248T, D255E, and E427K) would have not been detected when Af293 genome was used as reference since those changes are also found in its genome. Other important changes could have been missed during the analysis depending on which genome was used as a reference, highlighting the importance of choosing an appropriate reference genome for each case.
The highest number of SNVs was found in cluster IV, which means that was the furthest and most different population compared to both references (Table 3). In fact, this result was in agreement with those computed using MEGA. Cluster IV was the most distant (in number of differences as well as in proportion) from all the subclusters (Tables S6 and S7). When the differences within subclusters were analyzed (Tables S8 and S9), this cluster had the greatest homogeneity. In conclusion, cluster IV is the most different A. fumigatus population compared to the others but also the most homogeneously composed. The results obtained by MEGA to study the evolutionary distances among genomes were in line with the results obtained from WGS phylogenetic studies, since the evolutionary distances supported the formed clusters and subclusters (Figure 1 and Figure 2, Tables S6 and S7).
Regarding the Circos visual representation of the depth of coverage (Figure 3), it is noteworthy that the described large-scale deletions were detected in the edges of chromosomes 1 and 7, so these possible deletion events could be confined to telomeric or subtelomeric regions. Notably, all cluster III genomes did not have these deletions, as well as the reference genome used, which belonged to this same cluster. Remarkably, a huge deletion in chromosome 4 was present in all genomes. This might be due to a misassembly or to sequencing errors since two different Af293 genomes (downloaded and sequenced in-house) were included in the analysis and also had this deletion in spite of being the same strain that the genome used as reference. Similar results have been found by other authors although no comments were made about it, probably because it was difficult to realize without sequencing in-house the Af293 [13]. Areas with higher depth of coverage might be a consequence of possible duplication events, pseudogene presence or artefacts. However, since these regions appeared specifically in chromosomes 2, 4 and 6 in most of the genomes of cluster II, it would be worthy to explore them further. A specific mapping approach should be developed to analyze all these regions and to determine if there are real duplications and also if they represent a hallmark for cluster characterization. Further research is in progress to clarify which genes are in these differential regions and which are their biological processes or functions in which they are involved.
Experimental investigations on A. fumigatus biology and virulence have historically used only a few particular strains, of which Af293 and A1163 are the most used. More than in any other comparisons, the genetic background must be highly considered in the genomes used as references, since it can be responsible for the variability in many biological aspects between isolates and thus responsible for different phenotypes. The first genome comparison study between Af293 and A1163 was accomplished in 2007, showing that both genomes were highly syntenic in most regions, having core genes highly conserved and with levels of identity of 99.8% at the genome level [18]. These results were confirmed later by Fedorova et al. [17]. Specifically, the Af293 genome contained 143 [17] to 208 unique genes [18], while the A1163, 218 [17] to 320 [18]. In both cases, more than 60% of the strain-specific genes were related to cellular metabolism, secondary metabolism, signaling and transcriptional regulation [17,18]. These sequence divergences could explain the phenotypic differences described by other authors who reported a high heterogeneity among these two isolates regarding abiotic stimuli physiological responses, levels of growth, and virulence [15,38,41], as well as in immunogenic responses [42]. The variability in the total amount of unique genes in each genome could be due to sequencing errors or, less probably, to new mutations generated during the repeated subculture [13]. These results are in agreement with our study, in which the Af293 and the A1163 reference genomes belonged to different phylogenetic clusters which means that there are considerable differences between both genomes. Therefore, the appropriate election of the most suitable reference genome in each context will ultimately result in a more appropriate comparison and lead to solid conclusions.
Apart from the differences already discussed between Af293 and A1163, comparative analyses of the mating loci showed that both strains belonged to opposite mating types [40,43]. Mating type genes are not true idiomorphs, as they occupy adjacent positions on the chromosome. Mat1.1 strains have only one gene in this locus that encodes a transcriptional factor with an alpha box domain (AFUB_042900, Mat1.1.), while Mat1.2 strains contain a truncated copy of the high mobility group (HMG) box (Afu3g06170, Mat1.2.1) and another gene (Afu3g06160, Mat1.2.4) needed for heterokaryon formation as part of the mating process [44,45,46]. In this study we also checked if there were gene differences (presence or absence) other than those responsible for the mating type (Mat1.1, Mat1.2.1 and Mat1.2.4) among both groups of A. fumigatus genomes using an innovative approach (Figure S1). No gene absences were found in Mat1.1 strain genomes other than those encoded by the Mat1.2 locus itself (Mat1.2.1 and Mat1.2.4). Some previous works have reported differential gene presence between strains that coincidentally had different mating types [17,18]. However, the differences found must be due to their own different genetic background and not because of their opposite mating type.
One of the main concerns of invasive aspergillosis treatment is the global emergence of azole resistance. Regarding azole resistance mechanisms, Cyp51A protein is the target of these drugs and the cyp51A gene is a hotspot for mutations that confer azole resistance. The most common cyp51A modifications could be grouped into two categories: A. fumigatus strains that harbor cyp51A single point mutations (G54, G138, M220, or G448) [47,48,49,50], and isolates with specific point mutations in cyp51A gene together with various size tandem repeat (TR) integrations in the promoter of the gene (TR34/L98H, TR34/L98H/S297T/F495I, TR46/Y121F/T289A, or TR53) [24,51,52,53]. Moreover, all these cyp51A modifications have been described to evolve from two different azole resistance acquisition routes: the first set of strains come from the clinical setting as a consequence of the in-host drug adaptation after azole exposure in the patient [54], while the second set of isolates with TR insertions are hypothesized to develop from the azole exposure in the environmental setting [13,37,55,56,57,58]. Thanks to the high amount of isolates harboring different azole resistance mechanisms included in the study, a clear tendency in terms of grouping can be observed. In agreement with previous studies [13,37,55,56,58,59,60], all genomes included in our work that harbor a TR integration with single point mutations in cyp51A gene were grouped in cluster II, while isolates with punctual modifications in this gene were spread across clusters I and II. The same happened with cyp51A-WT strains that were distributed across those clusters. Therefore, our study supported the idea that isolates with punctual modifications in cyp51A (G54, G138, M220, or G448) and also cyp51A-WT strains have a greater genetic diversity than TR azole resistant isolates (TR34/L98H, TR34/L98H/S297T/F495I), which reinforces the previously suggested idea that TR resistance mechanisms have developed from a reduced set of clonally related strains with shorter genetic distances among them [59,60]. In this study, only genomes from strains which harbor the TR34/L98H resistance mechanism (with or without the S297T and F495I modifications) have been included. It would be most interesting to include isolates harboring TR46/Y121F/T289A, and TR53 azole resistance mechanisms. Furthermore, it is noteworthy that cluster III and IV were formed specifically by particular isolates that harbor 3 or 5 cyp51A modifications (F46Y, M172V, E427K and F46Y, M172V, N248T, D255E, E427K, respectively). This result was in agreement with previous studies developed in our group [19]. Further characterization of clado-specific genes would enhance our understanding of antifungal resistance mechanisms and help to unravel the environmental route resistance development.
Despite all these known azole resistance mechanisms, there is an increasing number of A. fumigatus azole resistant isolates for which the underlying mechanism remains unknown [61,62] or which are independent of cyp51A modifications [63]. The availability of A. fumigatus genomes to compare closely related strains but with differences in their azole susceptibility phenotype will increase the chance of finding genetic differences that could explain the phenotypic variations.
In conclusion, a great amount of novel and useful information can be derived from WGS studies. Here we particularly demonstrated that the selection of strains and reference genomes is crucial for comparative analysis at genomic and phenotypic level. Further compilation of A. fumigatus data generated by WGS studies could greatly enhance the understanding of molecular mechanisms involved in antifungal drug resistance and its development, as well as in many other biological functions that remain unknown. As the number of A. fumigatus genomes have increasingly been publically available, researchers will be able to increase the possibilities for data analysis which will ultimately allow them to come up with possible conclusions about an infinite number of hypotheses.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2073-4425/9/7/363/s1, Figure S1: Workflow for obtaining genes present in all Mat1.1 samples but not present in Mat1.2 (A); and on the contrary (B). Table S1: Mapping rate and coverage at 10× of 101 A. fumigatus strains against both reference genomes. Table S2: MEGA genetic distance estimations: pairwise distance matrix calculating the proportion of differences (p-model) and using Af293 as reference genome. Table S3: MEGA genetic distance estimations: pairwise distance matrix calculating the number of differences and using Af293 as reference genome. Table S4: MEGA genetic distance estimations: pairwise distance matrix calculating the proportion of differences (p-model) and using A1163 as reference genome. Table S5: MEGA genetic distance estimations: pairwise distance matrix calculating the number of differences and using A1163 as reference genome. Table S6: MEGA genetic distance estimations: mean distances within subclusters calculating the proportion and number of differences and using Af293 as reference genome. Table S7: MEGA genetic distance estimations: mean distances within subclusters calculating the proportion and number of differences and using A1163 as reference genome. Table S8: MEGA genetic distance estimations: matrix of the mean distances between subclusters calculating the proportion and number of differences and using Af293 as reference genome. Table S9: MEGA genetic distance estimations: matrix of the mean distances between subclusters calculating the proportion and number of differences and using A1163 as reference genome.

Author Contributions

Conceptualization, E.M.; Data curation, S.M. and I.C.; Formal analysis, R.G.-R. and S.M.; Funding acquisition, E.M.; Investigation, R.G.-R. and S.M.; Project administration, E.M.; Resources, E.M.; Software, S.M. and I.C.; Supervision, E.M.; Validation, I.C. and E.M.; Visualization, R.G.-R. and S.M.; Writing-original draft, R.G.-R. and E.M.; Writing-review & editing, R.G.-R., L.A.-F. and E.M.

Funding

This work has been supported by Fondo de Investigacion Sanitaria (FIS PI15/00019) and also by Plan Nacional de I+D+i 2013–2016 and Instituto de Salud Carlos III, Subdirección General de Redes y Centros de Investigación Cooperativa, Ministerio de Economía, Industria y Competitividad, Spanish Network for Research in Infectious Diseases (REIPI RD16/CIII/0004/0003), co-financed by European Development Regional Fund ERDF “A way to achieve Europe”, Operative program Intelligent Growth 2014–2020.

Conflicts of Interest

The authors declare no conflicts of interest

References

  1. Kontoyiannis, D.P.; Marr, K.A.; Park, B.J.; Alexander, B.D.; Anaissie, E.J.; Walsh, T.J.; Ito, J.; Andes, D.R.; Baddley, J.W.; Brown, J.M.; et al. Prospective surveillance for invasive fungal infections in hematopoietic stem cell transplant recipients, 2001–2006: Overview of the Transplant-Associated Infection Surveillance Network (TRANSNET) Database. Clin. Infect. Dis. 2010, 50, 1091–1100. [Google Scholar] [CrossRef] [PubMed]
  2. Pappas, P.G.; Alexander, B.D.; Andes, D.R.; Hadley, S.; Kauffman, C.A.; Freifeld, A.; Anaissie, E.J.; Brumble, L.M.; Herwaldt, L.; Ito, J.; et al. Invasive fungal infections among organ transplant recipients: Results of the Transplant-Associated Infection Surveillance Network (TRANSNET). Clin. Infect. Dis. 2010, 50, 1101–1111. [Google Scholar] [CrossRef] [PubMed]
  3. Denning, D.W. Invasive aspergillosis. Clin. Infect. Dis. 1998, 26, 781–805. [Google Scholar] [CrossRef] [PubMed]
  4. Latgé, J.P. Aspergillus fumigatus and aspergillosis. Clin. Microbiol. Rev. 1999, 12, 310–350. [Google Scholar] [PubMed]
  5. Kosmidis, C.; Denning, D.W. The clinical spectrum of pulmonary aspergillosis. Thorax 2015, 70, 270–277. [Google Scholar] [CrossRef] [PubMed]
  6. Verweij, P.E.; Chowdhary, A.; Melchers, W.J.G.; Meis, J.F. Azole Resistance in Aspergillus fumigatus: Can We Retain the Clinical Use of Mold-Active Antifungal Azoles? Clin. Infect. Dis. 2016, 62, 362–368. [Google Scholar] [CrossRef] [PubMed]
  7. Köser, C.U.; Ellington, M.J.; Peacock, S.J. Whole-genome sequencing to control antimicrobial resistance. Trends Genet. 2014, 30, 401–407. [Google Scholar] [CrossRef] [PubMed]
  8. Peacock, S. Bring microbial sequencing to hospitals. Nature 2014, 509, 557–559. [Google Scholar] [CrossRef] [PubMed]
  9. Kwong, J.C.; McCallum, N.; Sintchenko, V.; Howden, B.P. Whole genome sequencing in clinical and public health microbiology. Pathology 2015, 47, 199–210. [Google Scholar] [CrossRef] [PubMed]
  10. Zoll, J.; Snelders, E.; Verweij, P.E.; Melchers, W.J.G. Next-Generation Sequencing in the Mycology Lab. Curr. Fungal Infect. Rep. 2016, 10, 37–42. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Araujo, R. Towards the genotyping of fungi: Methods, benefits and challenges. Curr. Fungal Infect. Rep. 2014, 8, 203–210. [Google Scholar] [CrossRef]
  12. Hagiwara, D.; Takahashi, H.; Watanabe, A.; Takahashi-Nakaguchi, A.; Kawamoto, S.; Kamei, K.; Gonoi, T. Whole-genome comparison of Aspergillus fumigatus strains serially isolated from patients with aspergillosis. J. Clin. Microbiol. 2014, 52, 4202–4209. [Google Scholar] [CrossRef] [PubMed]
  13. Abdolrasouli, A.; Rhodes, J.; Beale, M.A.; Hagen, F.; Rogers, T.R.; Chowdhary, A.; Meis, J.F.; Armstrong-James, D.; Fisher, M.C. Genomic context of azole resistance mutations in Aspergillus fumigatus determined using whole-genome sequencing. MBio 2015, 6, e00536-15. [Google Scholar] [CrossRef] [PubMed]
  14. Wei, X.; Chen, P.; Gao, R.; Li, Y.; Zhang, A.; Liu, F.; Lu, L. Screening and characterization of a non-cyp51A mutation in an Aspergillus fumigatus cox10 strain conferring azole resistance. Antimicrob. Agents Chemother. 2017, 61, e02101-16. [Google Scholar] [CrossRef] [PubMed]
  15. Keller, N. Heterogeneity Confounds Establishment of “a” Model Microbial Strain. mBio 2017, 8, e00135-17. [Google Scholar] [CrossRef] [PubMed]
  16. Nierman, W.C.; Pain, A.; Anderson, M.J.; Wortman, J.R.; Kim, H.S.; Arroyo, J.; Berriman, M.; Abe, K.; Archer, D.B.; Bermejo, C.; et al. Genomic sequence of the pathogenic and allergenic filamentous fungus Aspergillus fumigatus. Nature 2005, 438, 1151–1156. [Google Scholar] [CrossRef] [PubMed]
  17. Fedorova, N.D.; Khaldi, N.; Joardar, V.S.; Maiti, R.; Amedeo, P.; Anderson, M.J.; Crabtree, J.; Silva, J.C.; Badger, J.H.; Albarraq, A.; et al. Genomic islands in the pathogenic filamentous fungus Aspergillus fumigatus. PLoS Genet. 2008, 4, e1000046. [Google Scholar] [CrossRef] [PubMed]
  18. Rokas, A.; Payne, G.; Fedorova, N.D.; Baker, S.E.; Machida, M.; Yu, J.; Georgianna, D.R.; Dean, R.A.; Bhatnagar, D.; Cleveland, T.E.; et al. What can comparative genomics tell us about species concepts in the genus Aspergillus? Stud. Mycol. 2007, 59, 11–17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Garcia-Rubio, R.; Alcazar-Fuoli, L.; Monteiro, M.C.; Monzon, S.; Cuesta, I.; Pelaez, T.; Mellado, E. Insight into the significance of Aspergillus fumigatus cyp51A polymorphisms. Antimicrob. Agents Chemother. 2018, 62, e00241-18. [Google Scholar] [CrossRef] [PubMed]
  20. Da Silva Ferreira, M.E.; Kress, M.R.V.Z.; Savoldi, M.; Goldman, M.H.S.; Härtl, A.; Heinekamp, T.; Brakhage, A.A.; Goldman, G.H. The akuBKU80 mutant deficient for nonhomologous end joining is a powerful tool for analyzing pathogenicity in Aspergillus fumigatus. Eukaryot. Cell 2006, 5, 207–211. [Google Scholar] [CrossRef] [PubMed]
  21. Howard, S.J.; Cerar, D.; Anderson, M.J.; Albarrag, A.; Fisher, M.C.; Pasqualotto, A.C.; Laverdiere, M.; Arendrup, M.C.; Perlin, D.S.; Denning, D.W. Frequency and evolution of Azole resistance in Aspergillus fumigatus associated with treatment failure. Emerg. Infect. Dis. 2009, 15, 1068–1076. [Google Scholar] [CrossRef] [PubMed]
  22. Takahashi-Nakaguchi, A.; Muraosa, Y.; Hagiwara, D.; Sakai, K.; Toyotome, T.; Watanabe, A.; Kawamoto, S.; Kamei, K.; Gonoi, T.; Takahashi, H. Genome sequence comparison of Aspergillus fumigatus strains isolated from patients with pulmonary aspergilloma and chronic necrotizing pulmonary aspergillosis. Med. Mycol. 2015, 53, 353–360. [Google Scholar] [CrossRef] [PubMed]
  23. Tang, C.M.; Cohen, J.; Holden, D.W. An Aspergillus fumigatus alkaline protease mutant constructed by gene disruption is deficient in extracellular elastase activity. Mol. Microbiol. 1992, 6, 1663–1671. [Google Scholar] [CrossRef] [PubMed]
  24. Alcazar-Fuoli, L.; Mellado, E.; Alastruey-Izquierdo, A.; Cuenca-Estrella, M.; Rodriguez-Tudela, J.L. Aspergillus section Fumigati: Antifungal susceptibility patterns and sequence-based identification. Antimicrob. Agents Chemother. 2008, 52, 1244–1251. [Google Scholar] [CrossRef] [PubMed]
  25. Mellado, E.; Garcia-Effron, G.; Alcazar-Fuoli, L.; Melchers, W.J.; Verweij, P.E.; Cuenca-Estrella, M.; Rodriguez-Tudela, J.L. A new Aspergillus fumigatus resistance mechanism conferring in vitro cross-resistance to azole antifungals involves a combination of cyp51A alterations. Antimicrob. Agents Chemother. 2007, 51, 1897–1904. [Google Scholar] [CrossRef] [PubMed]
  26. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [PubMed]
  27. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with Bowtie2. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef] [PubMed]
  28. Quinlan, A.R.; Hall, I.M. BEDTools: A flexible suite of utilities for comparing genomic features. Bioinformatics 2010, 26, 841–842. [Google Scholar] [CrossRef] [PubMed]
  29. McKenna, A.; Hanna, M.; Banks, E.; Sivachenko, A.; Cibulskis, K.; Kernytsky, A.; Garimella, K.; Altshuler, D.; Gabriel, S.; Daly, M.; et al. The Genome Analysis toolkit: A MapReduce framework for analyzing nextgeneration DNA sequencing data. Genome Res. 2010, 20, 1297–1303. [Google Scholar] [CrossRef] [PubMed]
  30. Stamatakis, A. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 2014, 30, 1312–1313. [Google Scholar] [CrossRef] [PubMed]
  31. Yu, G.; Smith, D.K.; Zhu, H.; Guan, Y.; Lam, T.T.-Y. ggtree: An R package for visualization and annotation of phylogenetic trees with their covariates and other associated data. Methods Ecol. Evol. 2017, 8, 28–36. [Google Scholar] [CrossRef]
  32. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis version 7.0. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [PubMed]
  33. Pearson, W.R.; Lipman, D.J. Improved tools for biological sequence comparison. Proc. Natl. Acad. Sci. USA 1988, 85, 2444–2448. [Google Scholar] [CrossRef] [PubMed]
  34. Nei, M.; Kumar, S. Molecular Evolution and Phylogenetics; Oxford University Press: New York, NY, USA, 2000; ISBN 9780195135855. [Google Scholar]
  35. Inouye, M.; Dashnow, H.; Raven, L.A.; Schultz, M.B.; Pope, B.J.; Tomita, T.; Zobel, J.; Holt, K.E. SRST2: Rapid genomic surveillance for public health and hospital microbiology labs. Genome Med. 2014, 6, 90. [Google Scholar] [CrossRef] [PubMed]
  36. Krzywinski, M.; Schein, J.; Birol, I.; Connors, J.; Gascoyne, R.; Horsman, D.; Jones, S.J.; Marra, M.A. Circos: An information aesthetic for comparative genomics. Genome Res. 2009, 19, 1639–1645. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Ashu, E.E.; Hagen, F.; Chowdhary, A.; Meis, J.F.; Xu, J. Global population genetic analysis of Aspergillus fumigatus. mSphere 2017, 2, e00019-17. [Google Scholar] [CrossRef] [PubMed]
  38. Kowalski, C.H.; Beattie, S.R.; Fuller, K.K.; McGurk, E.A.; Tang, Y.-W.; Hohl, T.M.; Obar, J.J.; Cramer, R.A., Jr. Heterogeneity among isolates reveals that fitness in low oxygen correlates with Aspergillus fumigatus virulence. mBio 2016, 7, e01515-16. [Google Scholar] [CrossRef] [PubMed]
  39. Popovich, K.J.; Snitkin, E.S. Whole Genome Sequencing—Implications for Infection Prevention and Outbreak Investigations. Curr. Infect. Dis. Rep. 2017, 19, 15. [Google Scholar] [CrossRef] [PubMed]
  40. Wortman, J.R.; Fedorova, N.; Crabtree, J.; Joardar, V.; Maiti, R.; Haas, B.J.; Amedeo, P.; Lee, E.; Angiuoli, S.V.; Jiang, B.; et al. Whole genome comparison of the A. fumigatus family. Med. Mycol. 2006, 44 (Suppl. 1), S3–S7. [Google Scholar] [CrossRef]
  41. Fuller, K.K.; Cramer, R.A.; Zegans, M.E.; Dunlap, J.C.; Loros, J.J. Aspergillus fumigatus photobiology illuminates the marked heterogeneity between isolates. mBio 2016, 7, e01517-16. [Google Scholar] [CrossRef] [PubMed]
  42. Rizzetto, L.; Giovannini, G.; Bromley, M.; Bowyer, P.; Romani, L.; Cavalieri, D. Strain dependent variation of immune responses to A. fumigatus: Definition of pathogenic species. PLoS ONE 2013, 8, e56651. [Google Scholar] [CrossRef] [PubMed]
  43. Losada, L.; Sugui, J.A.; Eckhaus, M.A.; Chang, Y.C.; Mounaud, S.; Figat, A.; Joardar, V.; Pakala, S.B.; Pakala, S.; Venepally, P. Genetic Analysis Using an Isogenic Mating Pair of Aspergillus fumigatus Identifies Azole Resistance Genes and Lack of MAT Locus’s Role in Virulence. PLoS Pathog. 2015, 11, e1004834. [Google Scholar] [CrossRef] [PubMed]
  44. Dyer, P.S.; O’Gorman, C.M. Sexual development and cryptic sexuality in fungi: Insights from Aspergillus species. FEMS Microbiol. Rev. 2012, 36, 165–192. [Google Scholar] [CrossRef] [PubMed]
  45. Dyer, P.S.; Inderbitzin, P.; Debuchy, R. Mating-type structure, function, regulation and evolution in the Pezizomycotina. In Growth, Differentiation and Sexuality. The Mycota (A Comprehensive Treatise on Fungi as Experimental Systems for Basic and Applied Research); Wendland, J., Ed.; Springer International Publishing: Basel, Switzerland, 2016; Volume I, pp. 351–385. ISBN 978-3-319-25844-7. [Google Scholar]
  46. Yu, Y.; Amich, J.; Will, C.; Eagle, C.E.; Dyer, P.S.; Krappmann, S. The novel Aspergillus fumigatus MAT1-2-4 mating-type gene is required for mating and cleistothecia formation. Fungal Genet. Biol. 2017, 108, 1–12. [Google Scholar] [CrossRef] [PubMed]
  47. Diaz-Guerra, T.M.; Mellado, E.; Cuenca-Estrella, M.; Rodriguez-Tudela, J.L. A point mutation in the 14alpha-sterol demethylase gene cyp51A contributes to itraconazole resistance in Aspergillus fumigatus. Antimicrob. Agents Chemother. 2003, 47, 1120–1124. [Google Scholar] [CrossRef] [PubMed]
  48. Mellado, E.; Garcia-Effron, G.; Alcazar-Fuoli, L.; Cuenca-Estrella, M.; Rodriguez-Tudela, J.L. Substitutions at methionine 220 in the 14-alpha sterol demethylase (Cyp51A) of Aspergillus fumigatus are responsible for resistance in vitro to azole antifungal drugs. Antimicrob. Agents Chemother. 2004, 48, 2747–2750. [Google Scholar] [CrossRef] [PubMed]
  49. Howard, S.J.; Webster, I.; Moore, C.B.; Gardiner, R.E.; Park, S.; Perlin, D.S.; Denning, D.W. Multi-azole resistance in Aspergillus fumigatus. Int. J. Antimicrob. Agents 2006, 28, 450–453. [Google Scholar] [CrossRef] [PubMed]
  50. Bellete, B.; Raberin, H.; Morel, J.; Flori, P.; Hafid, J.; Manhsung, R.T. Acquired resistance to voriconazole and itraconazole in a patient with pulmonary aspergilloma. Med. Mycol. 2010, 48, 197–200. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Hodiamont, C.J.; Dolman, K.M.; Ten Berge, I.J.; Melchers, W.J.; Verweij, P.E.; Pajkrt, D. Multiple-azole-resistant Aspergillus fumigatus osteomyelitis in a patient with chronic granulomatous disease successfully treated with long-term oral posaconazole and surgery. Med. Mycol. 2009, 47, 217–220. [Google Scholar] [CrossRef] [PubMed]
  52. Vermeulen, E.; Maertens, J.; Schoemans, H.; Lagrou, K. Azole-resistant Aspergillus fumigatus due to TR46/Y121F/T289A mutation emerging in Belgium, July 2012. Eur. Surveill. 2012, 17, 20326. [Google Scholar]
  53. Chen, Y.; Li, Z.; Han, X.; Tian, S.; Zhao, J.; Chen, F.; Su, X.; Zhao, J.; Zou, Z.; Gong, Y.; et al. Elevated MIC values to imidazole drugs among Aspergillus fumigatus isolates with TR34/L98H/S297T/F495I mutation. Antimicrob. Agents Chemother. 2018, 62, e01549-17. [Google Scholar] [CrossRef] [PubMed]
  54. Burgel, P.R.; Baixench, M.T.; Amsellem, M.; Audureau, E.; Chapron, J.; Kanaan, R.; Honoré, I.; Dupouy-Camet, J.; Dusser, D.; Klaassen, C.H.; et al. High prevalence of azole-resistant Aspergillus fumigatus in adults with cystic fibrosis exposed to itraconazole. Antimicrob. Agents Chemother. 2012, 56, 869–874. [Google Scholar] [CrossRef] [PubMed]
  55. Chowdhary, A.; Kathuria, S.; Xu, J.; Sharma, C.; Sundar, G.; Singh, P.K.; Gaur, S.N.; Hagen, F.; Klaassen, C.H.; Meis, J.F. Clonal expansion and emergence of environmental multiple-triazole-resistant Aspergillus fumigatus strains carrying the TR34/L98H mutations in the cyp51A gene in India. PLoS ONE 2012, 7, e52871. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Chang, H.; Ashu, E.; Sharma, C.; Kathuria, S.; Chowdhary, A.; Xu, J. Diversity and origins of Indian multi-triazole resistant strains of Aspergillus fumigatus. Mycoses 2016, 59, 450–466. [Google Scholar] [CrossRef] [PubMed]
  57. Garcia-Rubio, R.; Cuenca-Estrella, M.; Mellado, E. Triazole Resistance in Aspergillus Species: An Emerging Problem. Drugs 2017, 77, 599–613. [Google Scholar] [CrossRef] [PubMed]
  58. Wang, H.C.; Huang, J.C.; Lin, Y.H.; Chen, Y.H.; Hsieh, M.I.; Choi, P.C.; Lo, H.J.; Liu, W.L.; Hsu, C.S.; Shih, H.I.; et al. Prevalence, mechanisms and genetic relatedness of the human pathogenic fungus Aspergillus fumigatus exhibiting resistance to medical azoles in the environment of Taiwan. Environ. Microbiol. 2017, 20, 270–280. [Google Scholar] [CrossRef] [PubMed]
  59. Snelders, E.; Rijs, A.J.; Kema, G.H.; Melchers, W.J.; Verweij, P.E. Possible environmental origin of resistance of Aspergillus fumigatus to medical triazoles. Appl. Environ. Microbiol. 2009, 75, 4053–4057. [Google Scholar] [CrossRef] [PubMed]
  60. Camps, S.M.T.; Rijs, A.J.M.M.; Klaassen, C.H.W.; Meis, J.F.; O’Gorman, C.M.; Dyer, P.S.; Melchers, W.J.G.; Verweij, P.E. Molecular epidemiology of Aspergillus fumigatus isolates harboring the TR34/L98H azole resistance mechanism. J. Clin. Microbiol. 2012, 50, 2674–2680. [Google Scholar] [CrossRef] [PubMed]
  61. Chowdhary, A.; Sharma, C.; Hagen, F.; Meis, J.F. Exploring azole antifungal drug resistance in Aspergillus fumigatus with special reference to resistance mechanisms. Future Microbiol. 2014, 9, 697–711. [Google Scholar] [CrossRef] [PubMed]
  62. Moye-Rowley, W.S. Multiple mechanisms contribute to the development of clinically significant azole resistance in Aspergillus fumigatus. Front. Microbiol. 2015, 6, 70. [Google Scholar] [CrossRef] [PubMed]
  63. Fraczek, M.G.; Bromley, M.; Buied, A.; Moore, C.B.; Rajendran, R.; Rautemaa, R.; Ramage, G.; Denning, D.W.; Bowyer, P. The cdr1B efflux transporter is associated with non-cyp51 a mediated itraconazole resistance in Aspergillus fumigatus. J. Antimicrob. Chemother. 2013, 68, 1486–1496. [Google Scholar] [CrossRef] [PubMed]
Genes 09 00363 g001 550
Figure 1. Whole genome phylogenetic analysis of 101 Aspergillus fumigatus genomes included in the study. Dendrogram was performed using Af293 as reference genome.
Figure 1. Whole genome phylogenetic analysis of 101 Aspergillus fumigatus genomes included in the study. Dendrogram was performed using Af293 as reference genome.
Genes 09 00363 g001
Genes 09 00363 g002 550
Figure 2. Whole genome phylogenetic analysis of 101 Aspergillus fumigatus genomes included in the study. Dendrogram was performed using A1163 as reference genome.
Figure 2. Whole genome phylogenetic analysis of 101 Aspergillus fumigatus genomes included in the study. Dendrogram was performed using A1163 as reference genome.
Genes 09 00363 g002
Genes 09 00363 g003 550
Figure 3. Circos representation of whole-genome depth of coverage of 20 A. fumigatus genomes, averaged over 10,000-bp bins. Red marks showed null coverage, black marks a range of coverage between 1–500, and green marks showed a coverage higher than 500. Color legend shows the selected genomes of each cluster: I (blue), II (pink), III (orange) and IV (green).
Figure 3. Circos representation of whole-genome depth of coverage of 20 A. fumigatus genomes, averaged over 10,000-bp bins. Red marks showed null coverage, black marks a range of coverage between 1–500, and green marks showed a coverage higher than 500. Color legend shows the selected genomes of each cluster: I (blue), II (pink), III (orange) and IV (green).
Genes 09 00363 g003
Table
Table 1. Aspergillus fumigatus strains whole genome sequenced in our laboratory. The assigned lineage was based on the phylogenetic study.
Table 1. Aspergillus fumigatus strains whole genome sequenced in our laboratory. The assigned lineage was based on the phylogenetic study.
SamplesOrigincyp51A ModificationsAZL SCMating TypeSourceLineage
AF293SP5SNPsSM1.2CL3
akuBKU80SPWTSM1.1CL1.1
ATCC204305SPI242VSM1.1CL1.2
ATCC46645SPWTSM1.1CL2.1
CEA10SPWTSM1.1CL1.1
CM2141SPWTSM1.1CL2.2
CM237SPWTSM1.1CL1.1
CM2495SP3SNPsSM1.1CL4
CM2730SP3SNPsSM1.2CL4
CM2733SP3SNPsSM1.1CL4
CM3248SPN248KSM1.1CL1.1
CM3249SP3SNPsSM1.1CL4
CM3249bSP3SNPsSM1.1CL4
CM3262SP3SNPsSM1.1CL4
CM3720SP3SNPsSM1.1CL4
CM4602SP3SNPsSM1.2CL4
CM4946SP3SNPsSM1.2CL4
CM5419SPWTSM1.2CL2.2
CM5757SPWTSM1.2CL1.1
CM6126SPWTSM1.2CL1.1
CM6458SPWTSM1.2CL1.3
CM7510SPWTRM1.1CL1.1
CM7555SPWTRM1.2CL1.3
CM7560SP3SNPsSM1.1CL4
CM7570SP3SNPsSM1.1CL4
CM7632SP5SNPsSM1.1CL3
TP12SP5SNPsSM1.2CL3
TP32SP3SNPsSM1.1CL4
SP (Spain), SNPS (single nucleotide polymorphisms), 3SNPs cyp51A modifications (F46Y, M172V, E427K) and 5SNPs cyp51A modifications (F46Y, M172V, N248T, D255E, E427K), AZL susceptibility (azole susceptibility), S (susceptible), R (resistant), CL (clinical).
Table
Table 2. Aspergillus fumigatus genomes from public data bases included in this study. The assigned lineage is based on the phylogenetic study.
Table 2. Aspergillus fumigatus genomes from public data bases included in this study. The assigned lineage is based on the phylogenetic study.
SamplesOrigincyp51A ModificationsAZL SCMating TypeSourceLineageRef.
08-12-12-13NTTR34/L98H, S297T, F495IRM1.1CL2.2[13]
08-19-02-10NTTR34/L98HRM1.2ENV2.2[13]
08-19-02-30NTWTSM1.1ENV1.3[13]
08-19-02-46NTTR34/L98HRM1.1ENV2.2[13]
08-19-02-61NTTR34/L98HRM1.1ENV2.2[13]
08-31-08-91NTTR34/L98HRM1.1CL2.2[13]
08-36-03-25NTTR34/L98H, S297T, F495IRUnclearCL2.2[13]
09-7500806UKWTSM1.2CL1.3[13]
10-01-02-27NTTR34/L98HRM1.2CL2.2[13]
12-7504462UKWTSM1.2CL1.3[13]
12-7504652UKWTSM1.1CL2.1[13]
12-7505054UKWTSM1.1CL2.1[13]
12-7505220UKTR34/L98HRM1.1CL2.2[13]
12-7505446UKTR34/L98HRM1.2CL2.2[13]
Af293UK5SNPsSM1.2CL3[13]
AF41UKWTSM1.2CL1.1[13]
Af65UKWTSM1.2CL2.1[13]
AF72UKG54ERM1.2CL1.1[13]
AF90UKM220VRM1.1CL1.1[13]
Afu1042-09INTR34/L98HRM1.1CL2.2[13]
Afu124-E11INTR34/L98HRM1.1ENV2.2[13]
Afu166-E11INTR34/L98HRM1.1ENV2.2[13]
Afu218-E11INTR34/L98HRM1.1ENV2.2[13]
Afu257-E11INTR34/L98HRM1.1ENV2.2[13]
Afu343-P-11INTR34/L98HRUnclearCL2.2[13]
Afu591-12INTR34/L98HRM1.1CL2.2[13]
Afu942-09INTR34/L98HRM1.1CL2.2[13]
F12041UKG138CRM1.2CL2.1[21] *
F12219UKG54RRM1.2CL1.1[21] *
F12636UKG54ERM1.1CL2.1[21] *
F13535UKG138CRM1.2CL2.1[21] *
F13619UKH147Y, G448SRM1.1CL2.1[21] *
F13952UKG138CRM1.2CL2.1[21] *
F14403UKG54RRM1.1CL1.1[21] *
F14513GUKG138CRM1.2CL2.1[21] *
F14532UKM220TRM1.1CL2.1[21] *
F14946GUKWTSM1.2CL1.1[21] *
F15390UKM220TRM1.1CL2.1[21] *
F15927CN3SNPsSM1.1CL4[21] *
F16134DNM220KRM1.2CL1.1[21] *
F16216UKTR34/L98HRM1.2CL2.1[21] *
F16311UKWTSM1.1CL1.1[21] *
F17582UKWTSM1.2CL2.1[21] *
F17729UKWTSM1.2CL1.1[21] *
F17729WUKWTSM1.2CL1.1[21] *
F17764CNWTSM1.1CL1.1[21] *
F18085UKWTSM1.1CL1.1[21] *
F5211GUKWTSM1.1CL2.1[21] *
F7763UK5SNPsSM1.1CL3[21] *
IFM55369JPWTSM1.1CL1.2[22]
IFM58026JPN248KSM1.2CL1.1[22]
IFM58029JPWTSM1.2CL1.2[22]
IFM58401JPWTSM1.2CL1.1[22]
IFM59056JPWTSM1.1CL1.1[22]
IFM59073JPWTSM1.2CL1.2[22]
IFM59359JPWTSUnclearCL1.2[22]
IFM59361JPWTSM1.1CL1.3[22]
IFM59365JPWTSM1.2CL1.1[22]
IFM59777JPWTSM1.1CL1.2[22]
IFM60514JPWTSM1.2CL1.1[22]
IFM61118JPN248KSM1.1CL1.1[22]
IFM61407JPWTSUnclearCL1.1[22]
IFM61578JPP216LSM1.2CL1.1[22]
IFM61610JPWTSM1.2CL1.1[22]
IFM62516JPP329PSM1.1CL1.2[22]
SL143435PTN248KSM1.2UNK1.1*
SL143436PTWTSM1.2UNK1.1*
SL143437PTWTSM1.1UNK1.3*
SL143438PTN248KSM1.2UNK1.1*
SL143439PTWTSM1.1UNK1.3*
SL143440PT3SNPsSM1.1UNK4*
SL143441PTWTSM1.1UNK1.1*
SL146112PTN248KSM1.2UNK1.1*
NT (the Netherlands), UK (the United Kingdom), IN (India), CN (Canada), JP (Japan), PT (Portugal). SNPs (single nucleotide polymorphisms), 3SNPs cyp51A modifications (F46Y, M172V, E427K) and 5SNPs cyp51A modifications (F46Y, M172V, N248T, D255E, E427K), AZL SC (azole susceptibility), S (susceptible), R (resistant), CL (clinical), ENV (environmental), UNK (unknown), * NCBI SRA.
Table
Table 3. Mean of single nucleotide variants (SNVs). Variants within clusters and subclusters compared to both reference genomes.
Table 3. Mean of single nucleotide variants (SNVs). Variants within clusters and subclusters compared to both reference genomes.
Phylogenetic GroupsSNVs vs. A1163 Reference GenomeSNVs vs. Af293 Reference Genome
Subcluster I.134,54475,653
Subcluster I.240,12465,121
Subcluster I.346,81179,899
Subcluster II.181,47898,423
Subcluster II.280,99497,749
Cluster III82,23535,268
Cluster IV158,154159,742
Table
Table 4. Single nucleotide variant differences compared to both reference genomes. Variants of each A. fumigatus genome included in this study.
Table 4. Single nucleotide variant differences compared to both reference genomes. Variants of each A. fumigatus genome included in this study.
SamplesA1163Af293SamplesA1163Af293SamplesA1163Af293
08-12-12-1380,79191,956CM2495151,089148,524F1758279,59897,024
08-19-02-1075,40189,762CM2730157,505159,177F1772944,05479,106
08-19-02-3042,48787,831CM2733154,361160,494F17729W43,66378,416
08-19-02-4678,44996,243CM324822,83474,568F1776425,81270,013
08-19-02-6180,24496,416CM3249157,242159,900F1808551,24178,553
08-31-08-9188,97298,911CM3249b154,510157,780F5211G82,09199,726
08-36-03-2583,921102,006CM3262161,695160,926F776384,76455,332
09-750080645,28674,134CM3720156,600158,781IFM5536942,84964,003
10-01-02-2772,48793,299CM4602155,381161,767IFM5802622,18668,777
12-750446248,92775,301CM4946158,692161,082IFM5802936,52163,106
12-750465280,20393,323CM541989,35499,771IFM5840126,66963,489
12-750505486,952101,471CM575736,02774,420IFM5905630,84161,481
12-750522075,479100,229CM612655,87186,902IFM5907343,23068,775
12-750544686,194100,736CM645842,48271,670IFM5935939,30967,118
Af29383,486313CM751033,55979,283IFM5936142,79573,228
AF29378,045659CM755551,65777,679IFM5936521,83074,248
AF4134,87972,914CM7560158,657153,769IFM5977740,11765,300
Af6590,143103,614CM7570160,908162,292IFM6051439,62975,765
AF7235,27173,632CM763281,51925,490IFM6111821,90570,230
AF9040,26273,395F1204178,81293,712IFM6140732,02576,034
Afu1042-0980,34696,348F1221935,44073,872IFM6157826,20167,946
Afu124-E1181,11597,509F1263680,35998,501IFM6161041,69074,102
Afu166-E1180,94897,350F1353578,74193,777IFM6251641,22967,364
Afu218-E1180,57196,955F1361985,412107,985SL14343522,03781,031
Afu257-E1180,49997,008F1395278,46693,926SL14343652,60484,188
Afu343-P1173,797101,305F1440327,47674,454SL14343761,28991,394
Afu591-1280,78197,039F14513G78,64794,213SL14343822,67681,649
Afu942-0980,82797,149F1453279,931102,436SL14343939,56287,955
akuBKU80104378,152F14946G44,44279,719SL143440162,993165,788
ATCC 20430537,61460,180F1539079,743101,967SL14344138,49872,994
ATCC 4664584,994103,941F15927161,157164,076SL14611221,78281,137
CEA10100975,854F1613419,66077,718TP1283,36324,982
CM214188,704107,238F1621678,07390,736TP32163,367162,032
CM23748,40283,874F1631151,38778,630F1758279,59897,024
Table
Table 5. Types of SNVs and features. Variants found within each cluster compared to both genomes.
Table 5. Types of SNVs and features. Variants found within each cluster compared to both genomes.
VariantsA1163 Reference GenomeAf293 Reference Genome
IIIIIIIVIIIIIIIV
Frameshift variants256512477828446587133819
Inframe deletions14927826451124633094527
Inframe insertions19136929156132544791556
Intergenic variants640511,10712,819673013,42816,3943,7858792
Intron variants1568366834189066312043229608670
Missense variants434210,25110,00823,107876611,973256222,746
Protein altering81912181622421
Start lost91623572124657
Stop gained9523026429814121843288
Stop lost132732612533961
Stop retained102221521724552
Synonymous variants36759026900526,080818911,058236525,418

Share and Cite

MDPI and ACS Style

Garcia-Rubio, R.; Monzon, S.; Alcazar-Fuoli, L.; Cuesta, I.; Mellado, E. Genome-Wide Comparative Analysis of Aspergillus fumigatus Strains: The Reference Genome as a Matter of Concern. Genes 2018, 9, 363. https://0-doi-org.brum.beds.ac.uk/10.3390/genes9070363

AMA Style

Garcia-Rubio R, Monzon S, Alcazar-Fuoli L, Cuesta I, Mellado E. Genome-Wide Comparative Analysis of Aspergillus fumigatus Strains: The Reference Genome as a Matter of Concern. Genes. 2018; 9(7):363. https://0-doi-org.brum.beds.ac.uk/10.3390/genes9070363

Chicago/Turabian Style

Garcia-Rubio, Rocio, Sara Monzon, Laura Alcazar-Fuoli, Isabel Cuesta, and Emilia Mellado. 2018. "Genome-Wide Comparative Analysis of Aspergillus fumigatus Strains: The Reference Genome as a Matter of Concern" Genes 9, no. 7: 363. https://0-doi-org.brum.beds.ac.uk/10.3390/genes9070363

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop