Original Research ARTICLE
Whole Genome Sequencing of Australian Candida glabrata Isolates Reveals Genetic Diversity and Novel Sequence Types
- 1Centre for Infectious Diseases and Microbiology-Public Health, Westmead Hospital, Sydney, NSW, Australia
- 2Westmead Clinical School, Faculty of Medicine and Health, The University of Sydney, Sydney, NSW, Australia
- 3Centre for Infectious Diseases and Microbiology, Westmead Institute for Medical Research, Westmead, NSW, Australia
- 4Marie Bashir Institute for Emerging Infectious Diseases and Biosecurity, The University of Sydney, Sydney, NSW, Australia
- 5Department of Infectious Diseases and Microbiology, New South Wales Health Pathology, Royal Prince Alfred Hospital, Faculty of Medicine and Health, The University of Sydney, Sydney, NSW, Australia
- 6Centre for Infectious Diseases and Microbiology Laboratory Services, ICPMR, New South Wales Health Pathology, Westmead Hospital, Sydney, NSW, Australia
- 7National Mycology Reference Centre, SA Pathology, Adelaide, SA, Australia
- 8Department of Microbiology and Infectious Diseases, Canberra Hospital & Health Services, Australian National University Medical School, Canberra, ACT, Australia
- 9Department of Microbiology and Infectious Diseases, St Vincent’s Hospital, Sydney, NSW, Australia
- 10Department of Infectious Diseases, Alfred Health and Monash University, Melbourne, VIC, Australia
- 11Department of Microbiology, PathWest Laboratory Medicine, Queen Elizabeth II Medical Centre, Perth, WA, Australia
- 12Department of Microbiology and Infectious Diseases, Royal North Shore Hospital, Sydney, NSW, Australia
- 13National Centre for Infections in Cancer, Peter MacCallum Cancer Centre, Melbourne, VIC, Australia
Candida glabrata is a pathogen with reduced susceptibility to azoles and echinocandins. Analysis by traditional multilocus sequence typing (MLST) has recognized an increasing number of sequence types (STs), which vary with geography. Little is known about STs of C. glabrata in Australia. Here, we utilized whole genome sequencing (WGS) to study the genetic diversity of 51 Australian C. glabrata isolates and sought associations between STs over two time periods (2002–2004, 2010–2017), and with susceptibility to fluconazole by principal component analysis (PCA). Antifungal susceptibility was determined using Sensititre YeastOneTM Y010 methodology and WGS performed on the NextSeq 500 platform (Illumina) with in silico MLST STs inferred by WGS data. Single nucleotide polymorphisms (SNPs) in genes linked to echinocandin, azole and 5-fluorocytosine resistance were analyzed. Of 51 isolates, WGS identified 18 distinct STs including four novel STs (ST123, ST124, ST126, and ST127). Four STs accounted for 49% of isolates (ST3, 15.7%; ST83, 13.7%; ST7, 9.8%; ST26, 9.8%). Split-tree network analysis resolved isolates to terminal branches; many of these comprised multiple isolates from disparate geographic settings but four branches contained Australian isolates only. ST3 isolates were common in Europe, United States and now Australia, whilst ST8 and ST19, relatively frequent in the United States, were rare/absent amongst our isolates. There was no association between ST distribution (genomic similarity) and the two time periods or with fluconazole susceptibility. WGS identified mutations in the FKS1 (S629P) and FKS2 (S663P) genes in three, and one, echinocandin-resistant isolate(s), respectively. Both mutations confer phenotypic drug resistance. Twenty-five percent (13/51) of isolates were fluconazole-resistant (MIC ≥ 64 μg/ml) of which 9 (18%) had non wild-type MICs to voriconazole and posaconazole. Multiple SNPs were present in genes linked to azole resistance such as CgPDR1 and CgCDR1, as well as several in MSH2; however, SNPs occurred in both azole-susceptible and azole-resistant isolates. Although no particular SNP in these genes was definitively associated with resistance, azole-resistant/non-wild type isolates had a propensity to harbor SNPs resulting in amino acid substitutions in Pdr1 beyond the first 250 amino acid positions. The presence of SNPs may be markers of STs. Our study shows the value of WGS for high-resolution sequence typing of C. glabrata, discovery of novel STs and potential to monitor trends in genetic diversity. WGS assessment for echinocandin resistance augments phenotypic susceptibility testing.
The opportunistic yeast Candida glabrata is the second most common cause of candidemia and invasive candidiasis (IC) in many countries (Arendrup et al., 2013; Giunea, 2014; Pfaller et al., 2014; Chapman et al., 2017). Its clinical importance as a species lies in its reduced susceptibility to azole antifungal agents; more recently, resistance to the echinocandins as well as resistance to both these drug classes (Pfaller et al., 2012; Wisplinghoff et al., 2014; Shields et al., 2015). This has prompted much investigation of the epidemiology and biological properties of C. glabrata infections (Vale-Silva and Sanglard, 2015).
Because the prevalence of C. glabrata candidiasis and drug resistance rates varies both between and within geographical region, local epidemiological data are essential to inform management (Pfaller et al., 2012; Giunea, 2014). The reasons for this variation are uncertain but likely include prior exposure to azoles, patient factors and geographic location-specific determinants. Fungus-specific factors such as genetic strain variation within species are also pertinent. Delineation of intraspecies variation is useful not only to elucidate the molecular epidemiology of C. glabrata infections but also to assess potential transmission routes, biological niches and population structure. Yet relatively little is known about the genomic variation between isolates from different regions or the clinical significance of such differences.
Genetic typing methods e.g., pulsed-field gel electrophoresis, microsatellite analysis and multilocus sequence typing (MLST) have been used to determine genetic relatedness of C. glabrata (Dodgson et al., 2003; Lin et al., 2007; Abbes et al., 2012). In particular, the use of a standardized 6-locus MLST system (Dodgson et al., 2003) has improved discrimination between isolates with good reproducibility and portability of data via internet-accessible databases. Major findings from MLST analyses highlight that despite description of a broad range of MLST sequence types (STs), C. glabrata appears to be highly clonal with infrequent emergence of novel STs, which may be restricted to various geographical regions (Lott et al., 2012; Hou et al., 2017; Amanloo et al., 2018). This observed clonality however, may be fluid with temporal shifts of the major C. glabrata subtypes documented over time in one study (Lott et al., 2010). More discriminatory methods for pathogen discrimination such as next generation sequencing (NGS) offer new insights into C. glabrata genetics including its molecular epidemiology and population dynamics. Global spread of previously isolated populations was inferred from genomic data in a recent study (Carrete et al., 2018). In addition, NGS has been utilized to elucidate mechanisms of drug resistance in this species from the diagnostic laboratory perspective (Singh-Babak et al., 2012; Biswas et al., 2017a).
In Australia, we observed a 1.7-fold increase in the proportion of Candida bloodstream infections caused by C. glabrata over a decade (2004–2006 vs. 2014–2015) (Chen et al., 2006; Chapman et al., 2017). Our laboratory is increasingly using whole genome sequencing (WGS) approaches, in line with international trends, in public health practice and investigations of nosocomial infections (Besser et al., 2018). Here we applied WGS to investigate the genetic diversity of Australian C. glabrata strains across more than a decade and sought associations between the frequency of sequence types and two time periods and with drug susceptibility to fluconazole.
Materials and Methods
All isolates were obtained from our culture collection spanning 10–20 years and represent previous surveillance isolates for which research ethics approval had been obtained. The present study was a laboratory-based epidemiological study. No identifiable patient data or medical records were accessed.
Isolates and Identification
Fifty-two C. glabrata (sensu stricto) isolates were studied. These comprised C. glabrata ATCC 90030 and 51 C. glabrata isolates from Australia obtained through the culture collection at the Clinical Mycology Laboratory, Westmead Hospital and the Molecular Mycology Research Laboratory, Westmead Millennium Institute for Medical Research, Sydney. With the exception of two isolates recovered from the same patient 3 weeks apart, all isolates represented single patient episodes of IC. The majority (>90%) of isolates were from the jurisdictions of New South Wales and Victoria. All isolates were re-confirmed as C. glabrata sensu stricto by matrix-assisted laser desorption ionization-time of flight technique (MALDI-TOF MS) supplemented by ITS sequencing as required (White et al., 1990).
Susceptibility to antifungal agents were determined using the Sensititre®YeastOneTM YO10 methodology (TREK Diagnostics, Cleveland, OH, United States) according to Clinical and Laboratory Standards Institute (CLSI) methodology (Clinical and Laboratory Standards Institute [CLSI], 2017b). Candida parapsilosis ATCC 22019 and Candida krusei ATCC 6258 were the quality control strains. MIC values were interpreted according to CLSI M60 guidelines for fluconazole and the echinocandins (Clinical and Laboratory Standards Institute [CLSI], 2017a); where there are no clinical breakpoints (CBPs) (for voriconazole, posaconazole and amphotericin B), species-specific epidemiological cut-off values (ECVs) defined isolates as wild-type (WT) or non-WT (Clinical and Laboratory Standards Institute [CLSI], 2018). There are neither CBPs or ECVs for 5-fluorocytosine.
DNA Extraction and Library Preparation for Whole Genome Sequencing
Candida glabrata ATCC 90030 and the 51 clinical isolates were subcultured on Sabouraud’s dextrose agar for 48 h at 35°C prior to testing to ensure purity. Genomic DNA was extracted using the Wizard®Genomic DNA Purification kit (Promega, Alexandria, NSW, Australia) and the concentration was quantified by Quant-iTTM PicoGreen®dsDNA Assay Kit (Life Technologies, Mulgrave, VIC, Australia). The Nextera XT kit (Illumina, San Diego, CA, United States) was used to construct genomic libraries. Tagmentation, PCR amplification and cleanup, library normalization and pooling, and sequencing on the NextSeq 500 platform (Illumina) were carried out with 2 X 150-bp paired-end chemistry as previously described (Biswas et al., 2017b).
Whole Genome Sequencing Data Analysis
Sequence reads were deposited in the NCBI Sequence Read Archive (SRA: project number PRJNA480138) and were mapped against the reference genome of C. glabrata CBS138 (GenBank Accession 4 No. GCA_0002545.21).
Obtained sequence reads were mapped to each chromosome independently (using C glabrata CBS 138 chromosomes A to M as the reference) employing Stampy v1.0.23 (Lunter and Goodson, 2011) with pre-BWA alignment. Analysis of the mitochondria was not included. Variants were called using FreeBayes v1.1.0-dirty and filtered for read depth (minimum 20), read base quality (minimum Phred score 30), mapping quality (minimum 30) and proportion of reads supporting the variant (>0.9). All indels were excluded from the mapping-based analysis. An aligned mapped file for each chromosome was constructed for all isolates using an in-house script. All probable recombination events were identified using Gubbins (Croucher et al., 2015) and subsequently masked prior to concatenating all chromosome sequences into a single SNP alignment.
To infer the phylogenetic relationship of the Australian isolates, the best-fitting substitution model (TVM+F+ASC+R2) was selected with the Bayesian Information Criterion using ModelFinder implemented in IQ-Tree v.1.6.2 (Nguyen et al., 2015; Kalyaanamoorthy et al., 2017). A maximum likelihood tree was then reconstructed using IQ-Tree using 1000 ultrafast bootstrap replicates (Hoang et al., 2018).
To place our data into global context, all 38 publicly available C. glabrata Illumina short read sequence data from seven countries other than Australia were downloaded and included in the analysis (Havlesrud and Gaustad, 2017; Carrete et al., 2018; Table 2). A network approach using SplitsTree4 (Huson and Bryant, 2006) was employed to examine the relationships between our isolates and the isolates from other countries.
Single Polynucleotide Polymorphism of Genes Associated With Antifungal Resistance
All SNPs in genes known for their role in drug resistance in C. glabrata (Sanguinetti et al., 2005; Berila et al., 2009; Arendrup and Perlin, 2014) were manually curated in CLC Genomics Workbench (CLC Bio version 7.0, Arrhus, Denmark) with only non-synonymous SNPs reported. Genes examined included FKS1, FKS2, FKS3 (for echinocandin resistance) FCY1, FCY2, CgFPS1, CgFPS2 (5-fluorocytosine resistance), ERG9, ERG11, CgCDR1, and CgPDR1 (azole resistance) and MSH2 (for multi-drug resistance). Only non-synonymous SNPs with a minimum read depth coverage of 20, defined as high-quality (hq), were included in the analysis.
MLST Sequence Types and Principal Component Analysis
In silico MLST sequence types (STs), inferred from whole genome sequence data (genome types) were obtained from assembled contigs using SPAdes v22.214.171.124 (Bankevich et al., 2012) and MLST software (Seeman, 2017, T2). All obtained STs were subsequently confirmed using a read based approach implemented through SRST2. Four novel C. glabrata MLST types (See Results) were submitted to the C. glabrata MLST database (Dr. Andrew Dodgson3 accessed September 11, 2018; and now Professor Oliver Bader, accessed November 5, 2018), and designated as ST123, ST124, ST126, and ST127. Clustering of isolates by genetic similarity according to time of isolation (2002–2004 vs. 2010–2017) and drug susceptibility (susceptible-dose-dependent (S-DD) or resistant (R) to fluconazole (Clinical and Laboratory Standards Institute [CLSI], 2017a)) was examined by Principal Component Analysis (PCA). PCA was based on a pairwise SNP distance matrix calculated under the K80 model with pairwise deletion, using the native stats R v3.4.1 package (R Core Team, 2013), and visualized with ggplot2 (Wickham, 2009).
Of 51 Australian clinical isolates, 49 (94%) were cultured from blood. Thirteen (25%) isolates were from the time period 2002–2004, and 38 (75%) from, 2010–2017 (Table 1).
TABLE 1. Candida glabrata isolates, Australia: year of isolation, body site, sequence types and in vitro susceptibility to five antifungal agents.
Table 1 summarizes the MIC values of the isolates against five antifungal agents. All clinical isolates tested had low MICs against 5-fluorocytosine (≤0.12 μg/ml) and WT MICs (<2 μg/ml) against amphotericin B. Four (7.8%) isolates (strains WM_18.26, WM_18.24, WM_18.64, and WM_18.63) were resistant to caspofungin (MIC range 2 to >8 μg/ml; Table 1), and cross resistant to micafungin and anidulafungin (results not shown) (Clinical and Laboratory Standards Institute [CLSI], 2017a). Thirteen isolates (Strains WM_03.419, WM_03.449, WM_04.113, WM_05.113, WM_18.33, WM_18.40, WM_18.27, WM_18.42, WM_18.49, WM_18.51, WM_18.53, WM_18.65, and WM_18.66) were classified as resistant to fluconazole with MICs of ≥ 64 ug/ml (Clinical and Laboratory Standards Institute [CLSI], 2017a); all also had non-WT MICs to voriconazole whilst posaconazole MICs ranged from 2 to >8 μg/ml (Table 1). However, six additional isolates (strains WM_18.44, WM_18.48, WM_18.56, WM_18.57, WM_18.59, and WM_18.62) had non-WT MICs to voriconazole and posaconazole although were susceptible-dose dependent to fluconazole, and yet 11 other isolates had non-WT MICs only to voriconazole (Table 1).
Overall, an average of 95% of sequencing reads were mapped to the C. glabrata reference genome with a median read depth coverage of 75-fold. After mapping to each reference chromosome independently, a number of recombination events were identified (range, n = 320–776, lowest for chromosome B and highest for chromosome I – data not shown).
In silico MLST and Global Phylogeny
By WGS, there were 18 distinct STs defined based on the alleles from six genetic loci (FKS, LEU2, NMT1, TRP1, UGP1, and URA3) among 52 isolates (Table 1 and Supplementary Table S1), including four new STs (referred to as ST123, ST124, ST126, and ST127) not previously recognized by the C. glabrata MLST database (https://pubmlst.org/cglabrata/ see Supplementary Table S1 for the allele numbers of the new STs). Strain ATCC 90030 typed as ST10.
Of the known STs, the commonest ST amongst the Australian isolates was ST3 (8/51, 15.7% of isolates) followed by ST83 (7/51 and 13.7%), ST7 and ST26 (each n = 5, 9.8%). Collectively, these four STs were responsible for almost half (n = 25; 49%) of the isolates. The most common new ST was ST123 (n = 4 isolates). Eight (17.6%) STs (ST6, ST18, ST36, ST45, ST59, ST124, ST126, and ST127) were represented by only a single isolate (Table 1). Despite a relatively small number of isolates sequenced, the number of ST types was considerable compared to previous studies suggesting a relatively high genetic diversity within Australian C. glabrata isolates. Three isolates cultured from body sites other than blood were of ST10, ST16 and ST26 (Table 1). The two isolates from the same patient (WM_18.63 and WM_18.64) were both ST8.
In general, the whole genome data clustered broadly within determined MLST types but with greater intra-cluster resolution. The split tree network analysis (Figure 1) resolved isolates to at least nine terminal branches containing four or more isolates; many of these comprise multiple isolates clustering together, and with isolates representing disparate geographic locations. However, four branches comprised of Australian isolates only (Figure 1). To place local isolates into a more global perspective, the STs of isolates reported from seven other countries are shown in Table 2. Certain STs were common to isolates from the regions studied herein e.g., ST3 (Belgium, France, Germany, United States, Australia), whilst others were either more restricted, or were more prevalent to one or two countries e.g., ST6 in Norway and France, ST8 in the United States and continental Europe. The two isolates from Taiwan showed new STs. Of note, isolates of ST19 were absent from Australia (vs. 4/12 US isolates; Table 2).
FIGURE 1. Unrooted network tree depicting the association between Australian Candida glabrata isolates and international isolates from seven countries based on their whole genome sequences. All clusters in the tree have been represented by different sequence types (STs) except Norway 5, Norway 6, Taiwan 1 and Taiwan 2 which have previously unassigned (new) STs. New sequence types (STs) from Australia are ST123, ST124, ST126 and ST127. Isolates representing a particular ST in branches, which contain multiple STs, are put in circles. The colors depict isolates from different countries: Black, Australia; Green, Belgium; Blue, France; Brown, Germany; Purple, Italy; Pink, Norway; Yellow, Taiwan; Red, United States. The Australian isolates have names starting with WM_ and the international isolates where named according to the country of origin, all followed by a numerical scheme. For isolates from same country in a cluster, the country name was followed by numerical identities of the isolates separated by commas. For example, in ST7 cluster, WM_18.33, 03.308,689, 04.242,387 (where 18, 03 and 04 are years of isolation followed by isolate number).
TABLE 2. Candida glabrata isolates from countries other than Australia with known multi-locus sequence types as a comparison with sequence types of Australian isolates.
Genomic Similarity According to Period of Isolation and Susceptibility to Fluconazole
On analysis of the sequenced genomes by PCA, there was no temporal association between the two periods of isolation (2002–2004 and 2010–2017) and genomic similarity (as represented by ST distribution), or between genomic similarity and fluconazole susceptibility as measured by MICs (fluconazole S-DD: n = 39 isolates vs. fluconazole-resistant: n = 12) (Figures 2A,B).
FIGURE 2. Principle component analysis (PCA) of Candida glabrata genomic SNP distances following masking of recombination and (A) phenotypic fluconazole susceptibility or (B) period of isolation.
The phylogenetic relationship of the 51 Australian isolates was also reconstructed with high bootstrap support (Figure 3). There was no association between drug susceptibility to fluconazole, voriconazole, posaconazole or caspofungin (results were similar for anidulafungin and micafungin) or resistance, and phylogenetic clustering. Rather the analysis illustrated that resistance or non-WT MICs emerged at several time points along the phylogeny.
FIGURE 3. Maximum likelihood phylogeny of Australian Candida glabrata isolates. Bootstrap support values of less than 100 are indicated at the corresponding nodes. Sequence types (STs) and minimum inhibitory concentrations (MICs) for fluconazole (FLU), voriconazole (VOR), posaconazole (POS) and caspofungin (CAS) of each isolate is indicated to the right of the tree. Isolates with MIC values that are classed as “resistant” or as “non-wild type” are colored in red.
Analysis for SNPS Related to Drug Resistance and Their Relationship to STs
Of two isolates with caspofungin MICs of ≥ 8 μg/ml, one (strain WM_18.24; ST16) contained the FKS1 mutation leading to the amino acid substitution Ser629Pro, whilst strain WM_18.26 (ST10) harbored the FKS2 mutation Ser663Pro. There were no other SNPs in any of FKS1, FKS2 or FKS3 in both strains. Isolates WM_18.63, and WM-18.64 (both ST8) recovered from the same patient harbored the FKS1 mutation Ser629Pro as well as a FKS2 mutation Glu784Gly. The Glu784Gly mutation was not present in any echinocandin susceptible isolates. Other SNPs were present in FKS1 (Gly14Ser) and FKS2 (Thr926Pro) but only in isolates of ST3. Several SNPs in FKS3 (Supplementary Table S2) were present in both echinocandin-susceptible and echinocandin-resistant isolates.
The presence of SNPs in genes linked to 5-fluorocytosine resistance e.g., CgFPS1, CgFPS2 and CgFCY1 and CgFCY2 broadly, varied with ST with no SNPs observed in isolates of ST6, ST22, ST55, ST59 and ST123 (Supplementary Table S2). Fluconazole-resistant isolates (Table 1) in general harbored mutations in CgPDR1 and to a lesser extent in CgCDR1, but overall, SNPs in these genes and in other efflux pump genes, CgFLR1 and CgSNQ2 (data not shown) were also present in azole-susceptible/WT isolates (Supplementary Table S2). However, although no particular SNP was definitively linked to the resistance phenotype, 5/13 fluconazole-resistant isolates harbored Pdr1 amino acid substitutions in the region beyond the first 200–250 amino acids with no fluconazole-susceptible isolate containing such changes. In addition, 2/6 isolates with non-WT MICs to voriconazole and posaconazole (strains WM_18.48 and WM_ 18.62; fluconazole MIC both 32 μg/ml) also harbored mutations resulting in amino acid substitutions outside this region. Conversely, substitutions within the first 250 amino acid positions in Pdr1 were common to both azole resistant and susceptible isolates (Supplementary Table S2). Eight fluconazole resistant isolates however, did not demonstrate a PDR1 mutation. Isolates exhibiting pan-azole resistance, or which had non-WT MICs also had mutation in CgCDR1 His58Tyr (6/10 isolates), but the last was also present in azole-susceptible isolates.
SNPs occurred in isolates of diverse ST. There were no SNPs in ERG11 and the few SNPs observed in ERG9 were predominantly in isolates of ST3 and ST26.
SNPs in the MSH2 gene were observed in 19 of 51 (37%) isolates with three main locations of mutations – Val239Leu (9 isolates), Glu456Asp (7 isolates) and Leu269Phe (3 isolates), with two isolates (strains WM_18.63 and WM_18.64) having two mutations at Val239Leu and Ala942Thr (Supplementary Table S2). The same MSH2 mutations were found in azole susceptible as well as azole-resistant/non-WT isolates. Overall, SNPs were identified in isolates of diverse STs including in isolates of the new STs, ST123 and ST127 (Supplementary Table S2). Whilst the mutation Glu456Asp was found in isolates of 5 different STs and that of Val239Leu in 4 different STs, the Leu269Phe mutation was found only in isolates of ST16. The combination of Val239Leu and Ala942Thr were only identified in ST8 isolates, both of which were pan-echinocandin resistant.
Understanding the genomic diversity of C. glabrata and its antifungal susceptibility patterns is key to optimal management of infections caused by this problematic pathogen. The few studies that have examined the genetic variation of large culture collections have employed traditional MLST and indicate a predominantly clonal population structure with infrequent recombination (Dodgson et al., 2005). Prevalence of circulating STs also showed geographical bias (Dodgson et al., 2003; Hou et al., 2017; Amanloo et al., 2018). Hence, genetic variation amongst isolates from one region cannot be generalized to another. Here, we determined for the first time using a WGS approach, the relative frequency of endemic STs among 51 Australian C. glabrata isolates from two time periods, and verified the applicability of WGS to determine STs, STs by WGS clustered isolates within similar “ST” groupings as in silico MLST with good intra-cluster resolution (Figure 1).
Sequence typing demonstrated relatively “large” genetic diversity amongst Australian C. glabrata isolates, with just under half of the isolates represented by only four STs. The remaining STs, not only represented Australian specific STs (Figures 1, 3) but suggested an overall diversity within this pathogen that is greater than previously appreciated. In a nationwide Chinese study (Hou et al., 2017) of 411 isolates, a “new” ST sequence type was encountered approximately every 11 isolates compared to our study, which observed a “new” ST every 5 isolates. These, and our data emphasize the regional differences, with 75.9% of Chinese isolates comprised of ST7 and ST3 compared to only 25.5% of Australian isolates. Another recent report noted a predominance of ST3 and ST7 (70% isolates) in Korea (Byun et al., 2018) whilst in Iran, three STs (ST59, ST74, and ST7) accounted for 50% of isolates further supporting the notion of low intraspecies diversity within C. glabrata (Amanloo et al., 2018). Isolates belonging to ST3 have been reported with relative high frequency from Europe and Asia and now, from Australia (Dodgson et al., 2003; this study). The presence of strains with the same ST on different continents demonstrates that clones may have arisen from the same ancestor and disseminated globally followed by local adaptation.
Whilst ST5 isolates were reportedly common in Europe (Dodgson et al., 2003), this ST was not found amongst our Australian isolates nor amongst those from a more recent study of European and US isolates (Tables 1, 2). Conversely, isolates of ST7 appear uncommon in Europe and the United States but are more prevalent in Japan, Korea and China (Dodgson et al., 2003; Hou et al., 2017; Byun et al., 2018; Carrete et al., 2018). Strains of ST8, ST18, and ST19 were the commonest types in the United States (Dodgson et al., 2003) whereas we identified only one ST18 isolate, no ST19 isolates and two ST8 isolates (from the same patient) in Australia. Broadly, there are more common STs between Australian and Asian isolates than between Australian and US/European isolates likely reflecting the geographical proximity between Asia and Australia. Further studies involving a larger number of C. glabrata isolates to test this hypothesis would be of interest. The observed geographic variation among STs highlights the importance of acquiring local data.
MLST analysis of US C. glabrata isolates collected over three time periods between 1992–2009 revealed a relatively small number of STs with little genetic differentiation (Lott et al., 2010). In the present study, the results indicate that there is no evidence of genomic similarity or ST distribution among the sequences isolated in the two timeframes studied (Figure 2A), suggesting that there have been no marked shifts in the present STs. We further found no association between WGS STs and susceptibility or /resistance to fluconazole (Figure 2B), consistent with that reported by traditional MLST studies (Dodgson et al., 2003; Amanloo et al., 2018). The phylogenetic tree (Figure 3) suggests that both azole and echinocandin resistance may arise at multiple time points, independent of strain clustering. The small numbers of isolates are acknowledged as a study limitation.
Strain typing is essential for epidemiological investigation. MLST has the advantage of providing easily comparable results via internet-accessible databases (Dodgson et al., 2003), but does not adequately capture the breadth of genetic diversity and is not readily available in diagnostic laboratories. The present study has illustrated the utility of WGS to delineate genome variability in C. glabrata and importantly offers both superior discriminatory power and convenience. The costs of the two techniques are nearly identical – approximately AUD 50/sample for MLST and AUD 80/sample for WGS. With decreasing footprint and technological advances, however, cost reduction for WGS is anticipated.
The echinocandin resistance rate of 7.8% (4/51 isolates) in the present study is influenced by sampling bias, being <2% across Australia (Chapman et al., 2017), lower than that in the United States (Alexander et al., 2013; Shields et al., 2015). The FKS1 mutation Ser629Pro (in one isolate) and the mutation FKS2 Ser663Pro (three isolates) identified are among the most common in C. glabrata strains with high-level resistance phenotypes (Garcia-Effron et al., 2009; Arendrup and Perlin, 2014). None of the isolates harbored other well-known mutations that confer echinocandin resistance e.g., R665G, R636S, and F659Y (Zimbeck et al., 2010; Shields et al., 2015). Conversely, the role of FKS2 Glu748Gly in isolates WM_18.63 and WM_18.64 in echinocandin resistance remains uncertain as this SNP has not been previously described. Interestingly, the SNP was absent in genomes of all ST8 isolates from other countries in that cluster (Italy2, France9, USA3, USA4, and USA11).
Approximately 25% of C. glabrata isolates in our study were fluconazole-resistant, comparable to that in the United States (20–30%) (Castanheira et al., 2014). Genome-wide sequencing revealed mutations in several multidrug resistance transporter genes (Supplementary Table S2) (e.g., CgPDR1 and CgCDR1) that are associated with resistance through activation of drug efflux pumps (Ferrari et al., 2009; Rodrigues et al., 2014). Although SNPs in these genes were found in both azole-resistant/non wild-type as well as azole-susceptible/wild-type isolates, as previously reported, the alterations were located within the first ∼250 amino acids of Pdr1 in both susceptible and resistant isolates. It is noteworthy that 5/13 fluconazole-resistant isolates had other amino acid substitutions located outside this region of Pdr1 (after the first 250 amino acid positions) (Supplementary Table S2) (Ferrari et al., 2009; Tsai et al., 2010; this study). The reasons for fluconazole resistance in the remaining 8/13 isolates lacking Pdr1 mutations are uncertain. One possibility is that these isolates are petite mutants due to the loss of functional mitochondria, common in C. glabrata upon fluconazole exposure (Brun et al., 2004; Ferrari et al., 2009). However, effects of specific SNPs in these genes and the potential of petite mutants leading to fluconazole resistance phenotype needs to be confirmed by functional and gene expression analyses, which was beyond the scope of the present study. Through gene deletion studies, it has been recently shown that when functionally active, all of CDR1, PDR1 and SNQ2 contribute to high level resistance to azoles (Whaley et al., 2018). The absence of ERG11 mutations predominant in other Candida species, such as azole-resistant C. albicans is also important (Morio et al., 2010).
Mutations in the DNA mismatch repair gene MSH2 are reported to be a genetic driver of multi-drug resistance and about 55% of C. glabrata isolates are expected to contain MSH2 gene mutations. Whilst over 35% of our isolates harbored a mutation in MSH2, these were present in both azole-susceptible and azole-resistant (or non-WT) isolates. Hence it is possible that MSH2 mutations are more a marker of “ST” of C. glabrata rather than an indicator of drug resistance, having been linked to isolates of ST16 (Healey et al., 2016; Delliére et al., 2016); similarly, in our study, the mutation Leu269Phe was only present in isolates of ST16 (Supplementary Table S2). However, the remaining MSH2 mutations were found in several STs including two of the novel STs identified herein. Of note, the mutation combination of Val239Leu and Ala942Thr were only identified in ST8.
Limitations of the present study include the relatively small numbers of isolates analyzed, which may have precluded the identification of associations between ST and period of isolation. In addition, the majority of isolates were from blood. However, Lott et al. demonstrated that bloodstream isolates of C. glabrata were genetically indistinguishable from those colonizing the host (Lott et al., 2012). By using genome-wide information in 33 strains, Carrete et al. inferred the population structure of C. glabrata where strains were clustered into highly divergent clades but with the structure suggesting recent global spread of previously isolated populations (Carrete et al., 2018). WGS with its superior discrimination is well placed to provide additional clues to evolutionary traits in this species.
In conclusion, we have shown the value of a WGS approach for high resolution sequence typing, discovery of novel STs of C. glabrata, and the potential to monitor trends in genetic diversity. We envisage useful contribution of our data including that of four novel STs to the global sequence repository. Our results suggest that azole, as well as echinocandin, resistance may arise at multiple time points independent of strain clustering or of STs. WGS assessment for echinocandin resistance has good potential to augment phenotypic susceptibility testing methods. Further study by WGS of C. glabrata STs and their evolution over time is warranted.
SC, VS, WM, CB, SV, VM, and CH conceived and designed the research. CB, VM, QW, EM, CH, SV, SC, VS, and WM performed the experimental work and analysis of genomic data. CB, CH, WM, SV, SC, SK, KK, DM, MS, IA, CM, and TS provided isolates and analyzed the data. SC, CB, VM, and SV wrote the manuscript. All authors read, provided scientific critique, edited, and approved the manuscript.
Conflict of Interest Statement
SV, MS, DM, CM, KK, TS, SK, CH, and IA are part of antifungal advisory boards and untied grants from MSD Australia and Gilead Sciences Inc.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors thank Ms. Krystyna Maszewska for assistance with the subculturing of Candida glabrata isolates from the 2002–2004 culture collection. The authors also thank the Sydney Informatics Hub of the University of Sydney for providing the high performance computing resources for comparative genomics.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.02946/full#supplementary-material
Abbes, S., Sellami, H., Sellami, A., Hadrich, I., Amouri, I., Mahfoudh, N., et al. (2012). Candida glabrata strain relatedness by new microsatellite markers. Eur. J. Clin. Microbiol. Infect. Dis. 31, 83–91. doi: 10.1007/s10096-011-1280-4
Alexander, B. D., Johnson, M. D., Pfeiffer, C. D., Jiménez-Ortigosa, C., Catania, J., Booker, R., et al. (2013). Increasing echinocandin resistance in Candida glabrata: clinical failure correlates with presence of FKS mutations and elevated minimum inhibitory concentrations. Clin. Infect. Dis. 56, 1724–1732. doi: 10.1093/cid/cit136
Amanloo, S., Shams-Ghahfarokhi, M., Ghahri, M., Rassaghi-Abyaneh, M. (2018). Genotyping of clinical isolates of Candida glabrata from Iran by multilocus sequence typing and determination f population structure and drug resistance profile. Med. Mycol. 56, 207–215. doi: 10.1093/mmy/myx030
Arendrup, M. C., Dzajic, E., Jensen, E., Johansen, H. K., Kjaeldgaard, P., Knudsen, J. D., et al. (2013). Epidemiological changes with potential implication for antifungal prescription recommendations for fungemia: data from a nationwide fungemia surveillance programme. Clin. Microbiol. Infect. 19, E343–E354. doi: 10.1111/1469-0691.12212
Bankevich, A., Nurk, S., Antipov, D., Alexey, A. G., Mikhail, D., Alexander, S. K., et al. (2012). SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 19, 455–477. doi: 10.1089/cmb.2012.0021
Berila, N., Borecka, S., Dzugasova, V., Bojnansky, J., and Subik, J. (2009). Mutations in the CgPDR1 and CgERG11 genes in azole-resistant Candida glabrata clinical isolates from Slovakia. Int. J. Antimicrob. Agents 33, 574–578. doi: 10.1016/j.ijantimicag.2008.11.011
Besser, J., Carleton, H. A., Gerner-Smidt, P., Lindsey, R. L., and Trees, E. (2018). Next-generation sequencing technologies and their application to the study and control of bacterial infections. Clin. Microbiol. Infect. 24, 335–341. doi: 10.1016/j.cmi.2017.10.013
Biswas, C., Chen, S. C. A., Halliday, C., Kennedy, K., Playford, E. G., Marriott, D. J., et al. (2017a). Identification of genetic markers of resistance to echinocandins, azoles and 5-fluorocytosine in Candida glabrata by next-generation sequencing: a feasibility study. Clin. Microbiol. Infect. 23, 676.e7–676.e10. doi: 10.1016/j.cmi.2017.03.014
Biswas, C., Chen, S. C. A., Halliday, C., Martinez, E., Rockett, R., Wang, Q., et al. (2017b). Whole genome sequencing of Candida glabrata of markers of antifungal drug resistance∗. J. Vis. Expt. 130:e56714. doi: 10.3791/56714
Brun, S., Berges, T., Poupard, P., Vauzelle-Moreau, C., Renier, G., Chabasse, D., et al. (2004). Mechanisms of azole resistance in petite mutants of Candida glabrata. Antimicrob. Agents Chemother. 48, 1788–1796. doi: 10.1128/AAC.48.5.1788-1796.2004
Byun, S. A., Won, E. J., Kim, M. N., Lee, W. G., Lee, K., Lee, H. S. et al. (2018). Multilocus sequence typing (MLST) genotypes of Candida glabrata bloodstream isolates in korea: association with antifungal resistance, mutations in mismatch repair gene (Msh2), and clinical outcomes. Front. Microbiol. 9:1523. doi: 10.3389/fmicb.2018.01523
Carrete, L., Ksiezopolska, E., Pegueroles, C., Gómez-Molero, E., Saus, E., Iraola-Guzmán, S., et al. (2018). Patterns of genomic variation in the opportunistic pathogen Candida glabrata suggest the existence of mating and a secondary association with humans. Curr. Biol. 28:15.e7–27.e7. doi: 10.1016/j.cub.2017.11.027
Castanheira, M., Messer, S. A., Jones, R. N., and Farrell, D. J., Pfaller, M. A. (2014). Activity of echinocandins and triazoles against a contemporary (2012) worldwide collection of yeast and moulds collected from invasive infections. Int. J. Antimicrob. Agents 44, 320–326. doi: 10.1016/j.ijantimicag.2014.06.007
Chapman, B., Slavin, M., Marriott, D., Halliday, C., Kidd, S., Arthur, I., et al. (2017). Changing epidemiology of candidemia in Australia. J. Antimicrob. Chemother. 72, 1103–1108. doi: 10.1093/jac/dkw422
Chen, S., Slavin, M., Nguyen, Q., Marriott, D., Playford, E. G., Ellis, D., et al. (2006). Active surveillance for candidaemia, Australia. Emerg. Infect. Dis. 12, 1508–1516. doi: 10.3201/eid1210.060389
Clinical and Laboratory Standards Institute [CLSI] (2017a). Performance Standards For Antifungal Susceptibility Testing Of Yeasts, M60, 1st Edn. Wayne, PA: Clinical and Laboratory Standards Institute document
Clinical and Laboratory Standards Institute [CLSI] (2017b). Reference Method For Broth Dilution Antifungal Susceptibility Testing Of Yeasts: Fourth Edition, M27, 4th Edn. Wayne, PA: Clinical and Laboratory Standards Institute document.
Clinical and Laboratory Standards Institute [CLSI] (2018). Epidemiological Cutoff Values for Antifungal Susceptibility Testing, M59, 2nd Edn. Wayne, PA: Clinical and Laboratory Standards Institute document
Croucher, N. J., Page, A. J., Connor, T. R., Delany, A. J., Keane, J. A., Bentley, S. D., et al. (2015). Rapid phylogenetic analysis of large samples of recombinant bacterial whole genome sequences using Gubbins. Nucleic Acids Res. 43:e15. doi: 10.1093/nar/gku1196
Delliére, S., Healey, K., Gits-Muselli, M., Carrara, B., Barbaro, A., Guigue, N., et al. (2016). Fluconazole and echinocandin resistance of Candida glabrata correlates better with antifungal drug exposure rather than with MSH2 mutator genotype in a french cohort of patients harboring low rates of resistance. Front. Microbiol. 7:2038. doi: 10.3389/fmicb.2016.02038
Dodgson, A. R., Pujol, C., Denning, D. W., Soll, D. R., and Fox, A. J. (2003). Multilocus sequence typing of Candida glabrata reveals geographically enriched clades. J. Clin. Microbiol. 41, 5709–5717. doi: 10.1128/JCM.41.12.5709-5717.2003
Ferrari, S., Ischer, F., Calabrese, D., Posteraro, B., Sanguinetti, M., Fadda, G., et al. (2009). Gain of function mutations in CgPDR1 of Candida glabrata not only mediate antifungal resistance but also enhance virulence. PLoS Pathog. 5:e100268. doi: 10.1371/journal.ppat.1000268
Garcia-Effron, G., Lee, S., Park, S., Cleary, J. D., and Perlin, D. S. (2009). Effect of Candida glabrata FKS1 and FKS2 mutations on echinocandin sensitivity and kinetics of 1,3-ß-d-glucan synthase: implications for the existing susceptibility breakpoint. Antimicrob. Agents Chemother. 53, 3690–3699. doi: 10.1128/AAC.00443-09
Hou, X., Xiao, M., Chen, S. C. A., Kong, F., Wang, H., Chu, Y. Z., et al. (2017). Molecular epidemiology and antifungal susceptibility of Candida glabrata in China (August 2009 to July 2014): a multicentre study. Front. Microbiol. 8:880. doi: 10.3389/fmicb.2017.00880
Kalyaanamoorthy, S., Minh, B. Q., Wong, T. K. F., von Haeseler, A., and Jermiin, L. S. (2017). ModelFinder: fast model selection for accurate phylogenetic estimates. Nat. Methods 14, 587–589. doi: 10.1038/nmeth.4285
Lin, C. Y., Chen, Y. C., Lo, H. J., Chen, K. W., Li, S. Y. (2007). Assessment of Candida glabrata strain relatedness by pulse-field gel electrophoresis and multilocus sequence typing. J. Clin. Microbiol. 45, 2452–2459. doi: 10.1128/JCM.00699-07
Lott, T. J., Frade, J. P., and Lockhart, S. R. (2010). MLST analysis reveals both clonality and recombination in Candida glabrata bloodstream isolates for US surveillance. Euk. Cell 9:6190625. doi: 10.1128/EC.00002-10
Lott, T. J., Frade, J. P., Lyon, G. M., Iqbal, N., Lockhart, S. R. (2012). Bloodstream and noninvasive isolates of Candida glabrata have similar population structures and fluconazole susceptibilities. Med. Mycol. 50, 136–142. doi: 10.3109/13693786.2011.592153
Morio, F., Loge, C., Besse, B., Hennequin, C., and Le Pape, P. (2010). Screening for amino acid substitutions in the Candida albicans Erg11 protein of azole-susceptible and azole-resistant clinical isolates: new substitutions and a review of the literature. Diagn. Microbiol. Infect. Dis. 66, 373–384. doi: 10.1016/j.diagmicrobio.2009.11.006
Nguyen, L. T., Schmidt, H. A., von Haeseler, A., and Minh, B. Q. (2015). IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 32, 268–274. doi: 10.1093/molbev/msu300
Pfaller, M. A., Andes, D. R., Diekema, D. J., Horn, D. L., Reboli, A. C., Rotsein, C., et al. (2014). Epidemiology and outcomes of invasive candidiasis due to non-albicans species of candida in 2,496 patients: data from the prospective antifungal therapy (PATH) registry 2004-2008. PLoS One 9:e101510. doi: 10.1371/journal.pone.0101510
Pfaller, M. A., Castanheira, M., Lockhart, S. R., Ahlquist, A. M., Messer, S. A., and Jones, R. (2012). Frequency of decreased susceptibility and resistance to echinocandins among fluconazole- resistant bloodstream isolates of Candida glabrata. J. Clin. Microbiol. 50, 1199–1203. doi: 10.1128/JCM.06112-11
Sanguinetti, M., Posteraro, B., Fiori, B., Rannom, S., Torelli, R., and Fadda, G. (2005). Mechanisms of azole resistance in clinical isolates of Candida glabrata collected during a hospital survey of antifungal resistance. Antimicrob. Agents Chemother. 49, 668–679. doi: 10.1128/AAC.49.2.668-679.2005
Seeman, T. (2017). MLST—Scan Contig Files Against PubMLST Typing Schemes. Available at: https://github.com/tseemann/mlst.
Shields, R. K., Nguyen, M. H., Press, E. G., Cumbie, R., Driscoll, E., Pasculle, A. W., et al. (2015). Rate of FKS mutations among consecutive Candida isolates causing bloodstream infection. Antimicrob. Agents Chemother. 59, 7465–7470. doi: 10.1128/AAC.01973-15
Singh-Babak, S., Babak, T., Diezmann, S., Hill, J. A., Xie, J. L., Chen, Y.-L., et al. (2012). Global analysis of the evolution and mechanism of echinocandin resistance in Candida glabrata. PLoS Pathog. 8:e1002718. doi: 10.1371/journal.ppat.1002718
Tsai, H.-F., Sammons, L. R., Zhang, X., Suffis, S. D., Su, Q., and Myeers, T. G. M., et al. (2010). Microarray and molecular analyses of the azole resistance mechanism in Candida glabrata oropharyngeal isolates. Antimicrob. Agents Chemother. 54, 3308–3317. doi: 10.1128/AAC.00535-10
Whaley, S. G., Zhang, Q., Candle, K. E., and Rogers, P. D. (2018). Relative contribution of the ABC transporters Cdr1, Pdh1, and Snq2 to azole resistance in Candida glabrata. Antimcirob. Agents Chemother. 62:e01070-18 doi: 10.1128/AAC.01070-18
White, T., Bruns, T., Lee, S., Taylor, J. (1990). “Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics,” in Pcr Protocols: A Guide To Methods And Applications, eds M. Innis, D. Gelfand, J. Sninsky, T. White (New York, NY: Academic Press), 315–322.
Wisplinghoff, H., Ebbers, J., Geurtz, L., Stefanik, D., Major, Y., Edmond, M. B., et al. (2014). Nosocomial bloodstream infections due to Candida spp. in the USA: species distribution, clinical features and antifungal susceptibilities. Int. J. Antimicrob. Agents 43, 78–81. doi: 10.1016/j.ijantimicag.2013.09.005
Zimbeck, A. J., Iqbal, N., Ahlquist, A. M., and Farley, M. M., Harrison, L. H., Chiller, T., et al. (2010). FKS mutations and elevated echinocandin MIC values among Candida glabrata isolates from U.S. population-based surveillance. Antimicrob. Agents Chemother. 54, 5042–5047. doi: 10.1128/AAC.00836-10
Keywords: whole genome sequencing, Candida glabrata, MLST, sequence type, Australia
Citation: Biswas C, Marcelino VR, Van Hal S, Halliday C, Martinez E, Wang Q, Kidd S, Kennedy K, Marriott D, Morrissey CO, Arthur I, Weeks K, Slavin MA, Sorrell TC, Sintchenko V, Meyer W and Chen SC-A (2018) Whole Genome Sequencing of Australian Candida glabrata Isolates Reveals Genetic Diversity and Novel Sequence Types. Front. Microbiol. 9:2946. doi: 10.3389/fmicb.2018.02946
Received: 12 September 2018; Accepted: 16 November 2018;
Published: 03 December 2018.
Edited by:Miguel Cacho Teixeira, Universidade de Lisboa, Portugal
Reviewed by:Erika Shor, Public Health Research Institute (PHRI), United States
Kelley R. Healey, William Paterson University, United States
Copyright © 2018 Biswas, Marcelino, Van Hal, Halliday, Martinez, Wang, Kidd, Kennedy, Marriott, Morrissey, Arthur, Weeks, Slavin, Sorrell, Sintchenko, Meyer and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Sharon C.-A. Chen, email@example.com
†These authors have contributed equally to this work
‡Equal senior authors