Biofilms formed by Candida albicans bloodstream isolates display phenotypic and transcriptional heterogeneity that are associated with resistance and pathogenicity

1Bloodstream infections (BSI) caused by Candida species remain a frequent cause of morbidity and mortality, particularly within the immunocompromised population. Overall, Candida species have been identified as the most common fungal pathogen found in bloodstream infections in the United States, and are the fourth most common organism responsible for all BSI, and are the third most common within the intensive care unit (ICU). Candidaemia is often associated with the ability of Candida to adhere to and form biofilms on indwelling medical devices, such as central venous catheters (CVC) and prosthesis. Biofilms are a population of microorganisms attached to one another and/or a surface, surrounded by an extracellular matrix (ECM).

2A defining feature of biofilms is their resistance to antimicrobial therapy, with higher drug concentrations required to kill biofilms and their dispersed cells when compared to equivalent free-floating planktonic cells. Another feature of C. albicans biofilms is their enhanced pathogenicity. For example, cells detaching from biofilms have been shown to be more cytotoxic than their planktonic counterparts and significantly increase mortality within a murine model of infection. These observations have been demonstrated clinically, where a significant association was observed between C. albicans biofilm formation and mortality rates in candidaemia patients.

3Whilst there is growing evidence of the importance of Candida biofilms in clinical medicine, not all clinical isolates are able to form biofilms. There is therefore a fundamental gap in understanding exactly what drives biofilm formation and its clinical implications. Establishing methods to differentiate these isolates is challenging, as many studies rely on either metabolic assays or biomass, and these frequently use a variety of different substrates and media. Therefore, comparison between these studies is not possible, and further interpretation of the data to improve clinical management both for diagnostics and antifungal therapy is limited. The purpose of this study was therefore to investigate and characterise biofilm formation by clinical isolates of C. albicans using standard methodologies and subsequently analyse biofilm subsets phenotypically and transcriptionally. Here we report that C. albicans clinical isolates form biofilms that are heterogeneous, and this is associated with altered antifungal drug sensitivity and pathogenic potential.

4C. albicans bloodstream isolates displayed heterogeneity with respect to their biofilm biomass when grown in RPMI (Figure 1A). RPMI was shown to support the optimal growth of C. albicans over 24, 48 and 72 h (Additional file 1: Figure S1). Isolates were categorised as low biofilm formers (LBF) or high biofilm formers (HBF) if their biomass absorbance were less than the first quartile (Q1 OD570 = 0.565) or greater than the third quartile (Q3 OD570 = 1.682), respectively. Those isolates in between the first and third quartile (Q1-Q3) were defined as intermediate biofilm formers. When HBF were stained with crystal violet (cv), the extent of the biofilm formation was observed macroscopically, where the bottom of the well was clearly covered with cellular biomass (Figure 1A). In contrast, minimal staining was retained on isolates classed as LBF, as demonstrated by the well remaining almost colourless. We analysed a subset of isolates from the LBF and HBF group (n = 3) using dry weight measurements and confirmed our previous observations that biofilm biomass was significantly greater in isolates termed HBF (p = 0.0023) (Figure 1B). These differences are clearly evident when viewed under a SEM at low (Figure 1C [i, iii]) and high magnification (Figure 1C [ii, iv]). LBF isolates were characterised by a predominance of yeast cells and lack of hyphal cells (Figure 1C i, ii). In contrast, C. albicans HBF were highly filamentous with a multi-dimensional structure with very few yeast cells (Figure 1C iii, iv).

5The CSH of LBF and HBF isolates was quantified to determine whether it played a role in biofilm forming ability. Figure 2 illustrates that the hydrophobicity of an isolate significantly alters its ability to form a biofilm. For LBF the CSH increased by 32% and 31% in 4 h (p < 0.05) and 24 h (p < 0.0005) biofilms, respectively, when compared to planktonic cells. This trend was also observed in isolates with HBF where CSH increased by 50% in 4 h (p < 0.0001) biofilms and 81% in 24 h (p < 0.0001) biofilms, when compared with planktonic counterparts. Furthermore, a significant increase in CSH was found in isolates with HBF between early (4 h) and mature (24 h) phases of biofilm development, where hydrophobicity increased by 31% (p < 0.001). When the hydrophobicity of LBF and HBF was compared, CSH was significantly increased by 41% in HBF isolates at 24 h (p < 0.0001); however, no significant difference was observed between isolates with LBF and HBF in planktonic cells and 4 h biofilms.

6Sessile antifungal testing was performed on C. albicans isolates with LBF and HBF to determine if one group were more susceptible to VRZ or AMB treatment. VRZ was ineffective against all biofilms tested, showing no difference in activity against LBF and HBF (data not shown). However, a dose-dependent effect was evident in all isolates tested with AMB (Figure 3). Moreover, a significant difference was observed between LBF and HBF treated with 0.25 – 32 mg/L AMB (p < 0.05). LBF and HBF isolates both had a MIC50 of 0.25 mg/L AMB, yet isolates with LBF were significantly less viable than those with HBF at this concentration (p = 0.0307). In addition, LBF isolates achieved an ~80% kill at 4 mg/L, whereas HBF required 32 mg/L to reach the same kill. No significant differences were observed in the growth rates of either set of LBF and HBF isolates (data not shown).

7We next analysed the impact of the isolates ability to form biofilms based upon the severity of infection using a previously described G. mellonella model. The average rate of killing by three HBF, three LBF and a reference strain (SC5314) of C. albicans were calculated to plot a survival curve. Survival data showed a significant difference in the killing of larvae between HBF and LBF (p < 0.0001 [Figure 4A]). After 2 and 6 days, respectively, >50% and 100% larval death was recorded for HBF isolates, whereas larvae infected with LBF only achieved 20% killing after 7 days challenge. The reference strain SC5314 achieved 50% and 100% larval death by day 4 and 7, respectively. Similar kill rates to that of HBF were observed in the type strain however, when compared to LBF there was a significant difference in larval mortality (p = 0.0005).

8Host-pathogen interactions in this model were then investigated by microscopically observing the morphology of the infected larvae at 24, 48 and 72 h post-infection with C. albicans HBF and LBF (Figure 4B). At 24 h in both the LBF (Figure 4Bi) and HBF (Figure 4Biv), the nodule formation and melanin deposition were mainly observed under the cuticle and in the fat body, with mild to strong melanisation observed in the centre of the nodules, together with the presence of yeast cells and/or hyphae. The LBF nodules were smaller in dimension and dispersed mainly in the subcuticle area (Figure 4Bi), whereas the HBF nodules had a stronger melanisation with the tendency to converge in large aggregates, and were localised more deeply within the fat body (Figure 4Biv). At 48 h, the LBF were confined to the external part of the visceral organs, with a spot-like distribution (Figure 4Bii); whereas the HBF were found to display a pronounced filamentation all around the intestinal wall, with a PAS positive matrix visible surrounding the hyphae (Figure 4Bv). Furthermore at 72 h, there was a substantial invasion of both the gastrointestinal tract and the tracheal system with damaged gut epithelium, where yeast and hyphal cells both observed in the HBF infection (Figure 4Bvi). In contrast, a segmental invasion of the intestinal wall (Figure 4Biii) was observed with LBF infection and the progression of the infection was to a lesser extent than that by the HBF. Table 1 summarises the localisation and characterisation of the nodules with LBF and HBF infected larvae. Changes in the fat body morphology and composition including vacuolisation and haemocyte recruitment, were detected during the course of the infection and were more evident in the HBF.

9C. albicans clinical isolates defined as LBF and HBF were further assessed at a transcriptional level and the expression of genes related to biofilm formation was investigated. ACT1 was used as the housekeeping gene and was shown to be stably expressed throughout all biofilm conditions. Figure 5 illustrates the levels of gene expression of LBF (n = 10) versus HBF (n = 10) at both (A) 4 and (B) 24 h. Overall, the majority of the genes tested followed a trend of up-regulation in HBF compared to LBF. However, statistically significant differences were observed in the glycosylated mannoproteins MNN4 (p = 0.0313) and MNT2 (p = 0.0044) at 4 h, where expression was increased by ~2 fold. Furthermore, the resistance gene CDR1 was significantly increased in HBF by 4- and 6-fold at 4 h (p = 0.0113) and 24 h (p = 0.0239), respectively (Additional file 2: Table S1).

10Clustering the expression of 23 selected genes from 5 different functional groups in a heat map showed their relationship with one another and their variable expression in LBF and HBF over time (Figure 6). Here we found the adhesion genes ALS3 and HWP1 were closely related and highly expressed, particularly in HBF isolates at 4 and 24 h. Furthermore, genes from different functional groups were closely related to one another irrespective of whether LBF or HBF, such as the proteinase SAP5 and the adhesion genes ALS5 and EAP1. The remaining SAP genes were all closely related to one another, and interestingly the resistance gene MDR1 and the cell wall mannoprotein OCH1. Further analysis of SAP3 showed an increase in transcription within LBF at 24 h, despite no differences being observed at 4 h. In contrast, SAP5 expression was consistently high at 4 and 24 h within HBF.

11Analysis of Spearman rho coefficients found that out of 23 selected genes, 7 including those related to adhesion (ALS3, EAP1), filamentation (CPH1), hydrolytic enzymes (SAP5, PLB1) and resistance (CDR1 and ZAP1) showed a significant positive correlation (p < 0.05) with cv biomass data at 4 h (Table 2). Further analysis of the relationship between these seven genes and all the other genes tested presented various correlations (marked bold in Table 2). For example, PLB1 was significantly correlated with all other genes tested (94.11%) except ZAP1, followed by CPH1 (76.47%), SAP5 (76.47%), EAP1 (64.71%), CDR1 (41.17%), ALS3 (35.29%) and ZAP1 (11.76%). Correlation of individual genes with one another at the 4 h time point showed that PLB1, CPH1, MNN4 and HWP1 all correlated with 5 of the 7 key genes defined above. Notably, 24 h gene expression revealed very few significant correlations other than cv. In fact, a significant negative correlation was found between SAP3 and the biomass data (Rho = −0.465, p = 0.045). Furthermore, SAP3 was positively correlated with MNT2 (Rho = 0.468, p = 0.043) and SAP4 (Rho = 0.460, p = 0.048).

12C. albicans is an important pathogen. It is the fourth most common organism isolated from total bloodstream infections and continues to carry a high mortality. The presence of medical devices such as central venous catheters (CVC’s) are known to be important risk factors suggesting that biofilm formation is a key feature in the pathogenesis of candidaemia. The past decade has seen a significant leap in our knowledge and understanding of the biology of C. albicans biofilms, particularly with respect to the molecular basis of their development and homeostasis. However, in the clinical setting it is generally assumed that all C. albicans isolates have the capacity to form biofilms, but often with little regard to individual differences within the species when managing the infection. Here we demonstrate that C. albicans display heterogeneous biofilm characteristics, and these strain differences have important implications with respect to treatment and pathogenicity.

13Previous studies have reported very defined categories in their analysis of association with clinical outcomes, i.e. biofilm formers and non-biofilm formers. However, these important studies fail to take into account the heterogeneous nature of individual clinical isolates forming biofilms, which based on their metabolic XTT values can range from 0.125 to 1.358. When looking for clinical correlations with these phenotypes then important information can be missed, as the isolates at either end of the biofilm forming spectrum may lead to different clinical outcomes. We therefore aimed to evaluate and categorise C. albicans biofilms into distinct levels of biofilm formation to determine if their biological features were significantly different. Initially we categorised biofilms grown in RPMI using a biomass stain and followed this up with dry weight analysis, which differentiated clinical isolates into defined groupings. This approach was used in preference to metabolic assays due to the highly variable nature of XTT from strain to strain. Moreover, XTT’s mainstay usefulness is limited to antifungal drug testing of biofilms. Our classification, based initially on biomass, was supported by observations on a macro- and microscopic level where it was clear that numerous cells consisting of hyphae and yeasts were visible in HBF, whereas scant layers of yeast cells were observed for LBF. We also investigated CSH as an additional biofilm positive feature, as previous studies have also shown a link between biofilm biomass and CSH. This study confirms that CSH impacts different phases of biofilm development, which is in agreement with previous work where it was shown that cells dispersed from mature biofilms were more hydrophobic than those dispersed from earlier stages of biofilm development. Furthermore, it has been shown hydrophobic cells are more adherent, and therefore it is unsurprising that CSH was increased in HBF isolates. Based on our overall approach to biofilm categorisation we can be confident in the phenotypes selected for further detailed analysis. We do however concede there are caveats to defining levels of biofilm, and this requires further work and collaboration between groups to establish a standardised method.

14One of the key defining features of C. albicans biofilms is their insensitivity to sterol active antifungal agents. We examined azole treatment, which unsurprisingly demonstrated poor activity overall with no group differences, presumably through adaptive resistance mechanisms as previously described. Notably, AMB was less effective against HBF biofilms than LBF, which we hypothesise is due to the inability of the compound to permeate easily throughout the dense physical structure of the cells encased within ECM. We purposely excluded echinocandins from this study as these have been shown to be an effective anti-biofilm antifungals, therefore quantifying differences in activity against the two populations would be difficult. These observations may have implications to whether a patient responds to antifungal therapy, as Tumbarello and colleagues (2012) demonstrated that inadequate antifungal therapy (azoles) and the presence of an indwelling venous catheter were key predictors of patient mortality and hospital length of stay in patients infected with biofilm forming isolates. Guidelines have also suggested that removal of the catheter is an important factor in improving clinical outcomes, again supporting the notion that biofilm formation has a crucial role in clinical outcomes. Given the importance of these infections, efficient and appropriate treatment in candidaemia cases has been highlighted, as failure to treat quickly and effectively has profound consequences on mortality statistics.

15We decided to test the hypothesis that clinical isolates capable of forming robust biofilms were more pathogenic, which may be a reason for their apparent role in infections with increased mortality. Previous experimental work has shown that cells dispersed from biofilms are more cytotoxic and kill mice quicker than the equivalent planktonic cells. Using a G. mellonella model we showed that HBF isolates caused significantly greater mortality rates than LBF isolates, a finding supported elsewhere. In addition, another study investigated the virulence of C. albicans isolates with varying levels of biofilm formation and found that mice infected with a LBF had increased survival rates compared to isolates that were infected with HBF. Histological analysis of the infected larvae displayed similar cell morphology of yeast and filamentous hyphae as observed in SEM images of LBF and HBF, respectively. This is in agreement with a previous study that showed filamentation plays a role in killing G. mellonella larvae.

16Filamentous growth is a characteristic feature of C. albicans biofilm formation. Defective hyphal formation through deletion of EFG1 has been shown to lead to low levels of biofilm growth. Given our growing knowledge of key biofilm related genes we decided to investigate transcriptional changes to determine whether these are truly represented amongst clinical isolates, and therefore could be used as a more robust way to categorise biofilm formation and as potential diagnostic targets of HBF isolates. ACT1, the stably expressed housekeeping gene, as reported elsewhere, enabled these relative comparisons. Regulation of biofilm related genes were shown to influence an isolates biomass within clinical isolates, echoing work carried out by other groups. Cluster analysis of the selected biofilm related genes showed a good association with functional classes of genes, such as adhesins and proteinases, suggesting that both LBF and HBF had conserved pathways in the basic developmental phases of biofilm growth. However, individual gene expression profiles were inconclusive, showing very few clear independent significant differences, though gene expression proved interesting at 4 h. We investigated the overall biomass at 24 h and examined how 4 h gene expression related to this. Overall HWP1 was the most highly regulated at this time point in both LBF and HBF, as has been shown elsewhere, though no significant differences between the populations were observed. Seven other genes were however shown to have significant positive correlations with biomass. The most significant was CDR1, which was unsurprising as it has been shown to be transiently expressed in different biofilm studies, though does not correspond directly to antifungal resistance. PLB1 was showed to be significantly correlated with another 16 genes including biomass, though expression appeared constitutively low level within the biofilm, which is agreement with previous studies, and may have an accessory role in the degradation of host tissue alongside SAP’s. Of these, SAP5 was shown to be highly expressed in mature biofilms, and correlated with biomass and 13 other genes. We previously reported that SAP5 was associated with higher levels of expression in in vitro biofilms formed from denture stomatitis C. albicans isolates. In addition, Nailis and colleagues demonstrated its crucial role both in a reconstituted human epithelial model and within in vivo biofilms. Furthermore, the role of Sap5 in biofilm formation has recently been demonstrated in BSI, where its expression was significantly increased when compared to planktonic counterparts. Adhesins, such as ALS3 was also upregulated, which has previously been identified to be involved in biofilm formation, particularly at early stages of biofilm development (0–6 h), where C. albicans mutants lacking this gene produce sparse biofilms on catheter material in vitro. EAP1, though showing no clear independent association with biofilm formation per se, did show a clear correlation with biomass and 11 other genes. Its importance in biofilm formation has been reported previously. Of interest was the positive correlation with ZAP1 expression at 4 h, which is a negative regulator of matrix production. It did positively correlate with BCR1, the global regulator of biofilm formation, suggesting that the early interaction between their proteins may be important for downstream construction of the biofilm. Collectively the data highlighted the importance of looking at multiple genes at once opposed to single gene targets.

17Overall, we have categorised isolates based on biological properties relating to biofilm characteristics, and evaluated these in models of infection and treatment, where we have shown clear differences in virulence. In an attempt to create a molecular basis of categorising these strains we have used gene expression studies, and showed that individual gene expression analysis of the biofilm related genes to differentiate and categorise biofilm-forming isolates may be futile. Instead, we have shown that taking a defined panel of genes during early biofilm growth may be more informative. In particular, the panel of genes such as SAP5, HWP1, EAP1, PLB1 and CDR1 investigated in tandem could constitute an important step towards diagnostics of C. albicans biofilm formation, though the use of transcriptomics, such as RNA-Seq, may prove useful in identifying novel diagnostic targets. Further work is required to determine why some patients succumb to C. albicans biofilms whereas others do not, as the HBF isolates do have an increased pathogenic potential and are more difficult to manage with antifungal agents.

18Candida albicans SC5314, 3153A and a series of routine patient anonymised clinical bloodstream isolates (n = 42) collected under the approval of the NHS Scotland Caldicott Gaurdians from the Royal Hospital for Sick Children (Yorkhill Division), Glasgow, UK, as part of candidaemia epidemiology surveillance study. All clinical isolates obtained during this period were independently identified using Colorex Candida chromogenic plates (E&O Laboratories Ltd, Bonnybridge, UK) and were stored in Microbank® vials (Pro-Lab Diagnostics, Cheshire, UK) at −80°C until further use. These isolates were sub-cultured onto Sabouraud’s dextrose agar (SAB [Sigma-Aldrich, Dorset, UK]). Plates were incubated at 30°C for 48 h and maintained at 4°C. Isolates were propagated in yeast peptone dextrose (YPD) medium (Sigma-Aldrich, Dorset, UK), washed by centrifugation and resuspended in the appropriate media (Sigma-Aldrich) to the desired concentration, as described previously.

19All C. albicans clinical isolates (n = 42) were standardised to 1 × 106 cells/mL in RPMI-1640 and biofilms grown in flat-bottomed 96 well microtitre plates at 37°C for 24 h and biomass of each isolate assessed using the crystal violet (cv) assay as previously reported, and isolates were grouped based on their level of biomass distribution (OD570nm values). Isolates that fell below the 1st quartile (Q1) were classed as having low biofilm formation (LBF), strains with a biomass greater than the 3rd quartile (Q3) were deemed isolates with high biofilm formation (HBF), and those that lay in between were classified as intermediate biofilm formation (IBF Q2). C. albicans biomass was further assessed using dry weight measurements. Selected isolates with LBF and HBF were grown as biofilms in 12 well tissue culture plates for 24 h, as previously described, and the resulting biomass homogenised in 1 mL of PBS using a cell scraper (STARLAB, Milton Keynes, UK). This was then passed through a 0.22 μm filter disc (Satorius Stedim) using a vacuum and filters were dried at 40°C overnight before measuring each isolates dry weight. Uninoculated controls were used for background correction.

20Representative isolates were also grown within 12 well flat-bottomed tissue culture plates (Corning Incorporated, NY, USA) for 24 h before carefully washing with PBS, stained with cv and then digitally imaged (Canon IXUS 220 HS). For scanning electron microscopy (SEM) representative C. albicans clinical isolates defined as LBF and HBF were grown directly onto Thermanox™ coverslips (Nunc, Roskilde, Denmark), then processed and analysed as previously described.

21The cell surface hydrophobicity (CSH) was determined for selected C. albicans clinical isolates with LBF (n = 10) and HBF (n = 10). CSH was assessed using the microbial adhesion to hydrocarbon test, with a few modifications. Isolates were standardised to 1 × 106 cells/mL in RPMI-1640 and grown as biofilms for 4 and 24 h in 75 cm2 flasks (Nunc, Rochester, NY) at 37°C. These were then washed with PBS and the resultant biomass scraped off and homogenised in YPD. Cells were standardised (OD590nm 1.0) and cells transferred into a glass tube and overlaid with 1/5th volume of xylene. Contents were vortexed for 1 min and phases separated over 30 min. The aqueous phase was carefully removed and OD590nm measured. The percentage of hydrophobicity was calculated as ([OD590nm before xylene overlay - OD590nm after xylene overlay]/ OD590nm before xylene overlay) × 100%.

22Antifungal testing to determine minimum inhibitory concentration (MIC) of sessile cells was performed using voriconazole (VRZ) and amphotericin B (AMB) (Sigma-Aldrich, Gillingham, UK) on 24 h preformed biofilms, as previously described in flat-bottomed, 96 well microtitre plates. C. albicans LBF (n = 10) and HBF (n = 10) were tested in duplicate, on three separate occasions. Sessile minimum inhibitory concentrations (SMICs) were determined at 80% inhibition using an XTT (2,3-bis(2-methoxy-4-nitro-5-sulfo-phenyl)-2H-tetrazolium-5-caboxanilide) metabolic reduction assay.

23The pathogenicity of C. albicans isolates pre-defined as LBF (n = 3) and HBF (n = 3) were assessed using the G. mellonella killing assay, as described previously. This biological model has been shown previously to be useful in the study of fungal virulence. Sixth-instar G. mellonella larvae (Livefoods Direct Ltd, UK) were stored in the dark and used within 7 days of shipment. Ten random larvae with a bodyweight of between 200 to 300 mg were used for each group. Overnight YPD cultures of each isolate were washed and standardised to 1 × 107 cells/mL in PBS. Larvae were inoculated using a 50-μl Hamilton syringe with 26 g needle by injecting 10 μL aliquots (1 × 105 cells/larva) into the haemocoel, through the hindmost proleg. In addition, mock inoculated larvae pierced on the proleg with a sterile needle and a PBS inoculated control group were also included in each experiment. The infected larvae were placed in sterile petri dishes, incubated at 37°C and the number of dead larvae were scored daily. Larva was considered dead when it displayed no movement in response to touch together with a dark discolouration of the cuticle. Pathogenicity of LBF and HBF was assessed using a Kaplan-Meier plot with percentage survival monitored over 7 days.

24The morphology of the larvae infected with two isolates of C. albicans LBF and HBF was examined. Larvae were infected with the respective strains as described previously, and after post-infection (24, 48 and 72 h) larvae were fixed by a direct injection of formalin into the haemocoel and by formalin immersion at room temperature for 24 h. Paraffin embedded samples were then transversally sectioned into four-micron thickness using a microtome (microm HM 3335H, Thermoscientific). Sections were then stained with Periodic Acid-Schiff (PAS) to evaluate C. albicans infected cells. Whole larvae sections were examined for characterisation and localisation of nodules by light optical microscope visualisation (Leica microscope, model 020–519.502). Two larvae were processed for each isolate, carried out on three separate occasions.

25C. albicans clinical isolates exhibiting LBF (n = 10) and HBF (n = 10) were selected for the analysis of genes related to biofilm formation. Biofilms were grown in 24 well flat-bottomed plates for 4 and 24 h at 37°C, as described above. Following incubation, biofilms were washed with PBS, removed and homogenised using a bead beater, and RNA extracted using the TRIzol™ (Life Technologies, Paisley, UK) method as described previously by our group. Total RNA was DNase (Qiagen, Crawley, UK) treated and purified using an RNeasy MinElute clean up kit (Qiagen, Crawley, UK), as per manufacturer’s instructions. RNA was quantified and quality assessed using a NanoDrop spectrophotometer (ND-1000, ThermoScientific, Loughborough, UK). Next, cDNA was synthesised from 200 ng of total RNA using High Capacity RNA to cDNA kit (Life Technologies, Paisley, UK) in a MyCycler PCR machine (Bio-Rad Laboratories, Hertfordshire, UK), following manufacturers instructions.

26All primers utilised for this study for quantitative PCR (qPCR) were designed from their sequences obtained from the Candida Genome Database (CGD) website (http://www.candidagenome.org). The web-based primer design software program Primer3 (http://primer3.ut.ee/) was used. Primers were checked for specificity to C. albicans using the NIH-BLAST for short or exact nucleotide sequences (http://blast.ncbi.nlm.nih.gov/Blast.cgi). PCR amplification efficiencies of all designed primer sets were optimised prior to gene expression analysis, with efficiencies of 90-110% used in this study. Details of the oligonucleotides primers (Eurogentec, Southampton, UK) used in this study are listed in Table 3. 200 ng cDNA was used in a mastermix containing SYBR® GreenER™ (Life Technologies, Paisley, UK), UV-treated RNase-free water and forward/reverse primers (10 μM), following manufacturers’ instructions. Cycle conditions consisted of 2 min at 50°C, 10 min at 95°C and forty cycles of 15 s at 95°C and 60 s at 60°C. Each parameter (LBF n = 10, HBF n = 10 at 4 and 24 h) was analysed in duplicate using MxProP Quantitative PCR machine and MxProP 3000 software (Stratagene, Amsterdam, Netherlands) and controls consisted of reactions in which reverse transcriptase template were absent. Gene expression was calculated using the ΔCt method where the genes of interest were normalised to the housekeeping gene Act1.

27Differential expression of the selected genes from all isolates with LBF (n = 10) and HBF (n = 10) were assessed by clustering and heat map analysis using GenEx software (Exiqon, Vedbaek, Denmark). In brief, percentage expression data was pre-processed for log transformation and mean values calculated (n = 10 for both LBF and HBF) for each gene before heat map production. Each coloured cell in the heat map represents the variable expression of genes in LBF and HBF at 4 and 24 h time points. An increase in gene expression is represented by red and a down-regulation by green. Clustering techniques were used to show genes with similar expression patterns (co-regulated genes) in each set of isolates. The clustering was performed independently by average linkage and Euclidean distances used as a distance measure for both dimensions in the data.

28Graph production, data distribution and statistical analysis were performed using GraphPad Prism (version 4; La Jolla, CA, USA). After assessing whether data conformed to a normal distribution data were transformed where necessary and One Way Analysis of Variance (ANOVA) was used to investigate significant differences between independent groups. A Bonferroni post-test was used to determine statistically significant differences between groups. The G. mellonella survival curve was analysed using log rank test. Student t-tests were used to measure statistical differences between the two independent groups assessed in gene expression studies. Statistical significance was achieved if p < 0.05. IBM SPSS® (version 20) statistical analysis software was used for correlation analysis. Two-tailed Spearman rho correlation coefficient was determined separately for all 4 and 24 h selected genes expression versus 24 h biomass data. Genes that had a significant correlation with biomass were tested for correlations with the other genes as described above.

29BSI: Bloodstream infections; ICU: Intensive care unit; CVC: Central venous catheter; ECM: Extracellular matrix; AMB: Amphotericin b; LBF: Low biofilm former; HBF: High biofilm former; CSH: Cell surface hydrophobicity; VRZ: Voriconazole; RPMI: Roswell park memorial institute; YPD: Yeast peptone dextrose; FCS: Fetal calf serum; YNB: Yeast nitrogen base; PBS: Phosphate buffered saline; CV: Crystal violet; XTT: 2,3-bis(2-methoxy-4-nitro-5-sulfo-phenyl)-2H-tetrazolium-5-caboxanilide; SEM: Scanning electron microscopy; MIC: Minimum inhibitory concentration; qPCR: Quantitative polymerase chain reaction; OM: Original magnification; FB: Fat body; SC: Subcuticle; PI: Paraintestinal; PT: Paratracheal; SD: Standard deviation.

30None of the authors have any conflict of interest in publishing this work.

31LS and RR participated in the study design, carried out the experimental studies on biofilms, performed statistical analysis and were responsible for the manuscript. DFL participated in study design, assisted with statistical support and helped draft the manuscript. FP, MF, DT and EB participated in study design and undertook the experimental work related to the G. mellonella model. BLJ and CW contributed to study design, data analysis and supervised manuscript writing. KS participated in qPCR analysis and manuscript writing. CJN participated in study design, analysis and supervised manuscript writing. GR conceived the study, participated in study design, data analysis and was responsible for writing and submission of the final manuscript. All authors read and approved the manuscript.

32Infection and Immunity Research Group, Glasgow Dental School, School of Medicine, College of Medical, Veterinary and Life Sciences, University of Glasgow.