|Home | About | Journals | Submit | Contact Us | Français|
Like many bacteria, yeast species can form biofilms on several surfaces. Candida albicans colonizes the surfaces of catheters, prostheses, and epithelia, forming biofilms that are extremely resistant to antifungal drugs. We have used transcript profiling to investigate the specific properties of C. albicans biofilms. Biofilm and planktonic cultures produced under different conditions of nutrient flow, aerobiosis, or glucose concentration were compared by overall gene expression correlation. Correlation was much higher between biofilms than planktonic populations irrespective of the growth conditions, indicating that biofilm populations formed in different environments display very similar and specific transcript profiles. A first cluster of 325 differentially expressed genes was identified. In agreement with the overrepresentation of amino acid biosynthesis genes in this cluster, Gcn4p, a regulator of amino acid metabolism, was shown to be required for normal biofilm growth. To identify biofilm-related genes that are independent of mycelial development, we studied the transcriptome of biofilms produced by a wild-type, hypha-producing strain and a cph1/cph1 efg1/efg1 strain defective for hypha production. This analysis identified a cluster of 317 genes expressed independently of hypha formation, whereas 86 genes were dependent on mycelial development. Both sets revealed the activation of the sulfur-amino acid biosynthesis pathway as a feature of C. albicans biofilms.
The ability to grow in close association with surfaces is widespread among microorganisms. Despite the prevalence of so-called biofilm-like growth in most aqueous environments, it is only recently that the properties of surface-attached communities have been compared to those of free-living (planktonic) populations. Within biofilms, microorganisms lose their motility and produce an intercellular matrix of secreted materials and a degree of specialization is seen among individual cells (6, 32).
While bacteria have been the focus of many biofilm studies, little attention has been paid to fungal biofilms. Because many fungi are able to colonize nonaqueous environments, surface attachment is essential and free-living planktonic forms are usually restricted to dissemination processes. However, both forms are common in budding yeasts. The model yeast Saccharomyces cerevisiae adheres to plastic surfaces under low glucose concentrations and initiates biofilm formation (27). Pathogenic Candida species also establish well-developed biofilms. For example, Candida albicans colonizes several tissues such as vaginal and oral epithelia, developing a biofilm that, in immunocompromised patients, can disseminate into the bloodstream and cause fatal systemic infections (22). This species also forms biofilms on inert surfaces such as urinary or central venous catheters, dental prostheses, and other indwelling devices. These communities are often resistant to drug therapy and act as a source of reinfections (7).
C. albicans biofilm formation is initiated when planktonic cells adhere to a surface, presumably through cell wall-located adhesion molecules, and begin to aggregate into a microcolony. Candidate adhesion molecules are the proteins encoded by the ALS gene family, homologous to the mating-associated agglutinin of S. cerevisiae and involved in adhesion of C. albicans to host tissues (5, 15). Once a basal layer of cells is formed, an extracellular matrix of proteins and polysaccharides is produced to consolidate the early biofilm (5). C. albicans biofilms contain small amounts of extracellular material under static conditions, but copious quantities are secreted in dynamic models (12). The contribution of yeast or hyphal forms to the different layers varies with substrate, medium, or strain, but generally the inner layers are composed of yeast-like cells while the external ones include mostly hyphal forms (7).
Some properties of C. albicans biofilms differ markedly from those of planktonic populations. Their resistance to therapeutic drugs is remarkably high and remains high even when the biofilm is disrupted (16). Concomitantly increased expression of some drug efflux pumps (MDR1, CDR1, and CDR2) has been reported (24).
Even though some differential properties of biofilms have been observed and expression changes in a few genes have been reported, it is still unclear whether biofilm-like growth itself imposes a significant change in C. albicans physiology and whether this change is reflected in a particular pattern of gene expression. Genome-wide approaches are now available that allow this question to be fully investigated (2). Here we present for the first time evidence of a gene expression pattern specific to C. albicans biofilms, which is only marginally influenced by environmental parameters. We studied biofilm populations developed under different conditions of flow, oxygenation, and glucose concentration and compared their global transcript profiles to those of planktonic populations. Biofilm transcriptomes were remarkably similar under different environmental conditions and differed dramatically from those of planktonic populations. Hence, a detailed study of the expression patterns that accompany biofilm growth in C. albicans is presented.
The following strains of C. albicans were used: SC5314 (wild-type isolate), CAF2-1 (URA3/ura3::λimm434), and CAI4 (ura3::λimm434/ura3::λimm434) (9), GTC41 (CAI4 GCN4/gcn4::hisG-URA3-hisG), GTC43 (CAI4 gcn4::hisG-URA3-hisG/gcn4::hisG), and GTC65 (CAI4 gcn4::hisG/gcn4::hisG URA3 GCN4) (30); and CDB1 (CAI4 cph1::hisG/cph1::hisG efg1::hisG/efg1::hisG) (3).
The C. albicans strains were grown in 0.67% yeast nitrogen base (YNB; Difco) with 0.4 or 2% glucose. Strains CAF2-1, GTC41, GTC43, and GTC65 were cultured in YPD1/3 (0.33% yeast extract, 0.66% Bacto Peptone, 0.4% glucose) supplemented with histidine and arginine (20 mg/liter). Dilution of standard YPD was necessary to grow these strains in microfermenters. All media were supplemented with uridine (40 mg/liter) and adenine (20 mg/liter).
An inoculum was prepared from an early-stationary-phase culture grown in flasks at 30°C in an orbital shaker and diluted to an optical density at 600 nm (OD600) of 1. Biofilms b1, b2, b6, Bwt, and Bm were produced in microfermentors (11). These consist of a glass vessel with a 40-ml incubation chamber where two glass tubes are inserted to drive the entry of medium and air. Used medium is evacuated through a third tube. Medium flow is controlled by a recirculation pump (Ismatec) and pushed by the pressured air. Plastic slides (Thermanox; Nunc) glued to a glass spatula were immersed in the inoculum for 1 h at room temperature. After this adhesion period, the spatula was transferred to the chamber and incubated at 37°C (except for b6, which was incubated at 30°C), with the medium flow set to 0.6 ml/min and air supplied at 105 Pa; under these conditions, the growth of the planktonic phase is minimized and most of the cells remain on the spatula. Biofilm b3 was produced in wells of 96-well microtiter plates (Techno Plastic Products) as previously described (25). A 100-μl volume of the inoculum were incubated for 1 h in wells at room temperature; 200 μl of fresh medium was added, and the plates were incubated at 37°C for 48 h with shaking at 40 rpm. The planktonic phase was discarded, and the wells were washed with YNB before the biofilms were collected. Biofilm b4 was obtained using disks of Exacanal (MERCK-Eurolab) preadsorbed in calf serum for 24 h (5), rinsed, and incubated with the inoculum for 72 h at 37°C with shaking at 40 rpm. Disks were collected by filtration through gauze and washed with YNB, and cells were recovered. Planktonic cultures were produced in flasks as described for the biofilm preinoculum, except for p5, where cultures were incubated in microtiter plates. Cells collected from several microfermentors, microtiter plates, or flasks were pooled to obtain enough RNA for replicated transcript profiling experiments.
Biofilm formation in microtiter plates was quantified by the method of Ramage et al. (25). The XTT [2,3-bis(2-methoxy-4-nitro-5-sulfo-phenyl)-2H-tetrazolium-5-carboxanilide sodium salt] stock solution was prepared at 0.5 g/liter in phosphate-buffered saline PBS (pH 7.4), filter sterilized, and stored at −20°C. A colorimetric change resulting from XTT reduction was measured in a microtiter plate reader (Labsystem Multiscan RC) at 492 nm.
Total-RNA preparation, [α-33P]dCTP labeling, Northern blot hybridizations, and hybridization of macroarrays were carried out as previously described (19). Each RNA sample was independently labeled and hybridized a minimum of three times against swapped filters. Preliminary normalization, logarithmic transformation, and determination of correlation coefficients were performed with Excel 2000 (Microsoft). Data were loaded into GeneSpring version 5.0 (Silicon Genetics) and normalized in a per-chip (to the 50th percentile) and per-gene (to the median) normalization. Microarrays spotted in duplicate with probes for 5,907 C. albicans genes (information available at www.pasteur.fr/Galar_Fungail) were obtained from Eurogentec (Seraing, Belgium). For each of the four pair comparisons described in the text, RNA preparations were independently labeled with the two dyes and dyes were swapped in two replicate arrays. A total of 16 paired hybridizations were performed. Cy3/Cy5 labeling and hybridization were carried out as recommended by Eurogentec. Sixty pmoles of each dye were hybridized onto the array overnight at 42°C in a Corning chamber. Data were captured by GenePix Pro version 22.214.171.124. Each microarray was inspected to identify low-quality spots that were not included in the analysis. Low-quality hybridizations were excluded on the basis of the GenePix vital-statistic table. Median F635, F560, B635, and B560 of each spot were loaded into GeneSpring for calculation of signal and background intensities. Signals were log transformed and per-chip normalized by an intensity-dependent method (Lowess) applied to the print-tip region (35). Significance analysis was performed using GeneSpring version 5.0 as indicated above. Gene annotations were from http://genolist.pasteur.fr/CandidaDB. Expression data are available at http://www.pasteur.fr/recherche/unites/Galar_Fungail/Biofilm/biofilm.html.
Biofilm development in microfermentors was recorded with a Nikon Coolpix digital camera. Biofilms developed on Thermanox slides (wild type) or glass coverslips (cph1/cph1 efg1/efg1 strain) were examined by scanning electron microscopy (SEM) as described previously (23). Briefly, biofilm samples were fixed for 1 h in 0.07 M sodium cacodylate buffer (pH 7.3) containing 1.2% glutaraldehyde and 0.05% ruthenium red. The samples were then washed in the same buffer containing 0.05% ruthenium red and postfixed in 1% osmium tetroxide in cacodylate buffer, treated by the critical-point drying method, and observed on a Gemini DSM 982 scanning electron microscope. Transmission and scanning electron microscopy were performed by Brigitte Arbeille and Claude Lebos at the Laboratoire de Biologie Cellulaire et Microscopie Electronique, UFR Médecine, Tours, France.
A first objective of this study was to establish the extent to which C. albicans biofilm and planktonic populations differ at the transcriptome level. For this purpose, we performed multiple comparisons of biofilm and planktonic populations of C. albicans wild-type strain SC5314 obtained under different conditions of culture (Table (Table11).
Biofilm populations b1 and b2 were produced in a microfermentor (11). This model permits the recovery of pure biofilm cultures in milligram quantities, in a single and rapid way that does not interfere with RNA purification. Briefly, following adhesion of planktonic cells to plastic slides, the slides were incubated vertically at 37°C in a continuous flow of YNB-0.4% glucose (Fig. (Fig.1A).1A). After 48 h, the slides were covered by a ca. 3-mm-thick biofilm of fungal cells, developing a “cottony” surface with depressions and some filamentous projections (Fig. (Fig.1B).1B). Biofilm development continued for about 72 h. Therefore, the biofilms collected in samples b1 and b2 were representative of a growing population. Another biofilm (b6; Table Table1)1) was produced in microfermentors under different conditions of nutrient supply (YNB-2.0% glucose), time (72 h), and temperature (30°C).
In microfermentors, biofilms develop under a flow that supplies unlimited nutrients; air bubbling provides agitation and strongly aerobic conditions. Since flow parameters influence the secretion of matrix material (12), we also studied biofilms produced in static and more anaerobic environments. These were the surface of microtiter plates for sample b3 (Table (Table1)1) (25) and serum-treated catheters for sample b4 (Table (Table1)1) (5). The proportions of yeast and hyphal forms were variable in each model, with longer incubation times, higher temperature, and the inclusion of serum leading to increased hyphal content.
To rule out genes whose differential gene expression was due to differences in growth rate, we collected exponential (p2) and stationary-phase (p3) populations of planktonic cultures (0.4% glucose). Two further planktonic populations were produced in microtiter plates (p5) or in YNB-2% glucose (p4).
Transcript profiles were determined using macroarrays with probes for 2,002 C. albicans open reading frames that are a random representation of the whole genome (19). These arrays were the only tools available at the start of this study, but they also offer advantages for multiple transcriptome comparisons in that hybridizations are independently performed for each sample, thus limiting the bias due to the choice of pairs, as is the case in two-color experiments. To evaluate variations inherent in cell culture and RNA isolation, we first calculated the correlation between two replicate biofilms, b1 and b2. The Pearson coefficient for this comparison was 0.97, which is in agreement with previously reported correlations for two independent hybridizations of the same RNA preparation (19). This value also represents the limit for detecting similarity between two related populations. Correlation coefficients for the biofilm-versus-biofilm comparisons were greatest (0.97 to 0.80), followed by planktonic-versus-planktonic (0.85 to 0.54) and planktonic-versus-biofilm (0.81 to 0.54) comparisons (Fig. (Fig.2A).2A). Planktonic and biofilm cultures under the same glucose concentration (b6 versus p4) or produced in the same microtiter plate model (b3 versus p5) showed low correlation (0.63 and 0.65, respectively); however, two biofilms produced under such different conditions as in b2 and b4 were highly correlated (0.91). For most of the comparisons, differences in the growth conditions did not significantly affect the correlation between biofilms (6 of 10 comparisons produced correlation coefficients over 0.90), whereas they produced a noticeable dispersion between planktonic populations (none of them over 0.90). The relationships are illustrated in a dendrogram for clarity (Fig. (Fig.2B).2B). A cluster containing the biofilm populations was clearly separated from planktonic branches.
We applied significance analysis to detect genes with differential expression between biofilm and planktonic conditions. Normalized gene expression values for all the biofilm samples were averaged and compared to the average value for all the planktonic samples. A combination of nonparametric tests (P ≤ 0.02) was used for ranking the genes, and Benjamini's correction for the false discovery rate (FDR) was applied (31). Of 1,850 genes, 325 were found to be differentially expressed between biofilm and planktonic conditions (Table (Table2).2). The 325-gene cluster was divided by K-means into two subsets; the first contained 214 genes that were overexpressed in biofilms, and the second contained 111 genes that were underexpressed (see Supplementary Table A1 available at http://www.pasteur.fr/recherche/unites/Galar_Fungail/Biofilm/TableA1.xls). Differential expression of HXT62, ELF1, and SSB1 was confirmed by Northern hybridization (Fig. (Fig.3).3). To clarify the way in which regulated expression of these 325 genes contributes to the biofilm phenotype, C. albicans genes represented on the array were first assigned to 25 functional categories based on homology to S. cerevisiae genes (http://mips.gsf.de/proj/yeast/catalogues/funcat/index.html). The expected frequency of genes belonging to each category in a random group of 325 genes was compared to the observed frequency in the set of differentially expressed genes. The distribution of expected and observed frequencies differed significantly, as shown by a χ2 test. Rather than a nonspecific effect on all cell processes, the χ2 test suggests that some cellular functions are differentially expressed in biofilms (Table (Table2).2). In particular, the cluster is significantly enriched in genes involved in the following activities (categories are from Table Table22).
The largest difference between observed and expected frequencies was found in the protein synthesis category. The 34 genes included here were overexpressed. They encode 14 ribosomal proteins, several translation factors, aminoacyl-tRNA synthases, and protein turnover factors. This suggests that increased ribosome production and protein translation occur in biofilms.
Categories 1 and 3 are the metabolic categories most strongly represented in the overexpressed subset. Many of those genes are involved in the synthesis of aromatic (TRP4, TRP5, ARO3, and ARO4) or sulfur (MET3, ECM17/MET5 CYS4, CYS3, and SER33) amino acids. S. cerevisiae homologues of many genes (MET3, ILV3, SER33, HIS1, ARO3, ARO4, TRP4, TRP5, ADE8, and ADE12) in this category are subject to regulation by the transcriptional activator Gcn4p, a general regulator of amino acid metabolism (21). All these genes were overexpressed.
The overexpressed genes are involved in the biosynthesis of membrane-associated lipids such as ergosterol (ERG25 and ERG16), sphingolipids (SUR2 and FAT1), or mitochondrial phospholipids (PEL1). SEC14 encodes a phospholipid transfer protein that regulates several aspects of membrane traffic and is involved in the production of secretory vesicules from Golgi apparatus.
Many overexpressed genes encode general components of the transcriptional machinery, in particular RNA polymerase subunits. Other genes are involved in mRNA or rRNA processing and stability. Some specific transcriptional regulators were also overexpressed, e.g., MIG1. This gene plays a decisive role in glucose-mediated repression in S. cerevisiae, but its function in C. albicans is not fully equivalent (36).
Of the 13 differentially expressed genes classified here, 6 (ODC1, HCA4, ECM17, CIS2, PLC1, and SPT14) directly control cell wall biogenesis. ODC1 is a key regulator of polyamine biosynthesis and plays a pivotal role in cell wall organization, since polyamines compensate for anionic charges in this fungal structure. C. albicans ODC1 deletion leads to a defect in hyphal formation (13). Interestingly, expression of PTK2, encoding a protein kinase required for efficient polyamine uptake, was also significantly elevated in biofilms. IPF3518 (VPS16 homologue) and VPS35 are both included in the underexpressed cluster and encode vacuolar sorting proteins. In contrast, IPF13042, which encodes the C. albicans homologue of S. cerevisiae Sys3p involved at the vesicle-docking stage between endosome and distal Golgi, is overexpressed.
We analyzed the expression profiles of CDR, MDR, and ALS genes whose differential expression was previously demonstrated in biofilms (5, 24). CDR1, CDR2, CDR3, CDR4, and MDR1 (BMR1) were not differentially expressed in our study. Transcription of ALS genes has been found in a polymethylmethacrylate model using a probe that recognizes several members of the family (5). Our analysis identified ALS1 as the major overexpressed gene. ALS1 was dramatically overexpressed in all the biofilm populations (average B/P expression ratio = 19.04; see Supplementary Table A1 for detailed values), in agreement with the observation by Chandra et al. (5) and with the role of similar adhesins (FLO11) in S. cerevisiae biofilms (27). In contrast, ALS7 was underexpressed and expression of ALS10 and ALS5 was not significantly altered. Thus, each ALS family member makes independent contributions to the biofilm transcriptome.
Based on the fact that a high proportion of the 214-gene set are amino acid biosynthesis genes and homologues of genes that are Gcn4p targets in S. cerevisiae (14, 21), a role for general control of amino acid biosynthesis in biofilm development was envisaged. Like S. cerevisiae Gcn4p, C. albicans Gcn4p controls many amino acid starvation responses, some of which lead to hyphal development (30). We tested the ability of a C. albicans gcn4/gcn4 strain to form biofilms in a microfermentor. Since we observed decreased planktonic growth of the gcn4/gcn4 strain in minimal medium (data not shown), presumably because this medium imposes amino acid starvation, we assayed biofilm formation in the microfermentor under conditions in which the growth rate and morphology of the two strains are equivalent. In rich medium (YPD1/3), the total biomass produced by the two strains was similar (Fig. (Fig.4A)4A) and so were the proportions between yeast and hyphal forms (data not shown). However, the proportion of biofilm-growing cells was significantly lower in the mutant (Fig. (Fig.4B).4B). Quantification of biofilm growth in microtiter plates (Fig. (Fig.4C)4C) also confirmed the requirement for a functional Gcn4p for normal biofilm formation.
Our next objective was to study overlaps in the cellular functions required for hyphal and biofilm development. Dimorphism is not an absolute prerequisite for biofilm formation but can influence the structure of C. albicans biofilms (1). Since Efg1p and Cph1p are two major regulators of the yeast-to-hypha transition and their simultaneous inactivation prevents switching under most hypha-inducing conditions (8), we decided to use the cph1/cph1 efg1/efg1 mutant for a biofilm formation assay.
The C. albicans CPH1/CPH1 EFG1/EFG1 CAI4 and cph1/cph1 efg1/efg1 CDB1 strains were precontacted against plastic slides and incubated in microfermentors. Surprisingly, after 48 h the mutant had developed a thin, haze-like biofilm on the microfermentor glass wall and the glass spatula rather than the plastic slide itself. In contrast, the wild-type strain produced a thick, cottony biofilm which colonized the plastic surface (Fig. (Fig.1B1B and C, panels 1). When incubated directly in contact with the spatula, the wild-type strain was unable to colonize the glass surface.
SEM of the biofilms formed by both strains showed that wild-type biofilms are composed of microcolonies connected by “bridges” of hyphae (Fig. (Fig.1B,1B, panel 2) and contain blastospores, hyphae, and pseudohyphae in variable proportions (Fig. (Fig.1B,1B, panel 3). In contrast, cph1/cph1 efg1/efg1 biofilms are less reticulated, as a reflection of the lack of hyphal bridges, and contain mostly blastospores (Fig. (Fig.1C,1C, panels 2 and 3). Hyphae emerged occasionally as short and nonbranched tips, but most cells remained as slightly elongated yeasts, as observed for the mutant under hypha-inducing conditions (17). The structure of the first layer of cells appears highly organized: cells orientate very close to each other, in a parallel manner (Fig. (Fig.1C,1C, panel 3). All of these traits result in a quite “compact” but thin biofilm, whose resistance to detachment is similar to that of the wild type.
For both strains, budding of the new blastospores (Fig. (Fig.1B1B and C, panels 4) was frequently observed adjacent to the first bud (coaxial pattern), similar to that in planktonic cultures. Remarkably, however, in biofilms some hyphae also emerge coaxially while in planktonic cultures hyphal budding is bipolar or lateral. This often results in the production of bouquet-like structures (Fig. (Fig.1B,1B, panel 4) that might facilitate the three-dimensional development of the biofilm.
Similar to Ramage et al. (26), we have observed that the cph1/cph1 efg1/efg1 strain is unable to form a biofilm on a plastic surface. However, our results suggest that the CPH1 and EFG1 genes are not essential for biofilm formation since this can progress on a glass surface.
We decided to exploit the altered cph1/cph1 efg1/efg1 biofilm phenotype to extend our transcriptome analysis of biofilm-intrinsic expression patterns that are unrelated to hyphal development. Indeed, a minimum set of genes required for biofilm development should be expressed in the mutant biofilm despite the filamentation defect. RNAs were prepared from duplicate cultures of wild-type and mutant biofilms as in b1 (Bwt and Bm, respectively [Table [Table1])1]) and from duplicate planktonic cultures as in p3 (Pwt and Pm [Table [Table1]).1]). Four independent comparisons, Bwt/Pwt, Bm/Pm, Bwt/Bm, and Pwt/Pm, were performed in which two samples for each comparison were simultaneously hybridized on a microarray with duplicated probes for 5,907 C. albicans genes.
Significance analysis was performed using GeneSpring version 5.0 only on genes with a log ratio above 1.5 in at least one hybridization. This cutoff was introduced to control the stringency of the FDR correction, which is directly influenced by the number of genes analyzed (31). The number of differentially expressed genes is shown in Table Table33 (see also Supplementary Table A2, available at http://www.pasteur.fr/recherche/unites/Galar_Fungail/Biofilm/TableA2.xls). As expected, a considerable number of genes (748 and 856 genes) were affected in biofilm-versus-planktonic comparisons. Since the cph1/cph1 efg1/efg1 mutation has a remarkable impact on biofilm structure, expression of many genes (386 genes) was also altered in the biofilm-versus-biofilm comparison. Thirty-three of these genes were further affected by the mutation under planktonic conditions, indicating a global effect of the mutation rather than a biofilm-dependent effect.
We used Venn diagrams to delimit each group of differentially expressed genes. The intersection between the Bwt/Pwt and the Bm/Pm sets comprised 442 genes (Fig. (Fig.5A)5A) whose planktonic-versus-biofilm expression changed in both strains. However, some of these genes may also alter their expression in response to other conditions. We therefore excluded from the 442-gene set a group of 39 genes whose expression was also dependent on the cph1/cph1 efg1/efg1 mutation during planktonic growth (Fig. (Fig.5B).5B). The intersection of the remaining 403 genes with the set of genes changing between biofilms (Bwt/Bm, 386 genes) identified 86 shared genes and 317 unshared genes. The two subsets are proposed as potential biofilm-specific candidates since they are (i) differentially expressed in biofilms independently of whether these contain hyphae and (ii) not affected by the mutation under planktonic conditions. Besides, expression of the 317-gene set was not significantly altered when mycelial development in the biofilm was impaired. In contrast, the 86-gene set identified genes whose expression changed between biofilms that contained hyphae and those that did not. The pertinence of this approach was confirmed by the exclusion from the 317-gene subset of some of the biofilm-related genes mentioned above that have also been identified as hypha specific or Efg1p or Cph1p dependent. For instance, ALS1, a downstream effector of the Efg1p filamentation pathway (10), was not included in either the 317- or 86-gene set, since its expression is not induced in mutant biofilms; some other hypha-specific genes (20) are discussed below.
The intersection between the two B/P sets also recognized many biofilm-related genes identified by our first approach: 72 genes from the 325-gene set were included in the intersection, 49 of them were included in the 317-gene set, and 17 were included in the 86-gene set (Supplementary Table A2).
A χ2 test revealed significant alterations in the distribution of functional categories in the 317-gene set (Table (Table2).2). As observed for the 325-gene set, the most overrepresented categories are protein synthesis (category 12) and amino acid and nucleotide metabolism (categories 1 and 3). The group includes six genes (MET16, MET10, CYS4, CYS3, MET6, and SAM2) involved in sulfur amino acid biosynthesis and salvage pathways (Fig. (Fig.6).6). Carbohydrate metabolism and energy were overrepresented with respect to the 325-gene set, mostly due to the presence of several genes for the glycolytic pathway (PGI1, FBA1, TPI1, PGK1, GPM1, ENO1, and PDC1), 2 genes responsible for the conversion of the glycolytic intermediate dihydroxyacetone-phosphate to glycerol (GPD1 and RHR2), and the catabolic repressor MIG1. All these genes were overexpressed. Consistent with this, the underexpressed genes include known targets of glucose repression (SNF3, MLS1, and TPS3). The study of the representation of functional categories also underlined several genes related to phosphate transport (category 4; PHO84, PHO87, and PHO88).
Inspection of the 86-gene set revealed another 3 genes in the sulfur amino acid pathways (MET3, MET15, and SAH1), completing the biosynthesis and salvage route (Fig. (Fig.6).6). These genes, like those in the 317-gene set, were overexpressed in both B/P comparisons and their expression was higher in mutant than in wild-type biofilms. In contrast, the agglutinin genes ALS12 and ALS4 and the TUP1-repressed genes RBT2 and RBT5 were overexpressed in both B/P comparisons but lower in the cph1/cph1 efg1/efg1 biofilm than in the hypha-producing wild-type biofilm. Some of these genes respond to several types of hypha-switching stimuli under planktonic conditions (4). Of the 46 genes identified by Nantel et al. (20) as modulated during the yeast-to-hypha transition induced by serum and temperature, 6 [DDR48, IPF3912(CUP9), TYE7, IFU3(YLR63), IPF5915(YKR70), and HSP12] are included in the 86-gene set (Fig. (Fig.5).5). Since we did not expose B and P samples to serum, elevated temperature, or other environmental conditions known to induce hyphal switching, the inclusion of known hypha-modulated genes (20) in the 86-gene cluster and, more broadly, in the set of genes that show differential expression between wild-type and mutant biofilms (Fig. (Fig.5;5; Bwt-versus-Bm set) is likely to reflect a response to a switching stimulus that is specific for biofilms.
In most natural environments, microorganisms follow a sessile way of life in which they remain surface attached in close contact with each other. As a result, such populations exhibit a number of features that are very different from those of planktonic communities. Many forms of candidiasis develop on organ or implant surfaces. Consequently, C. albicans populations can follow a surface-attached, biofilm-like way of life that confers, among other characteristics, an increased resistance to antifungal drugs (7).
In this study we used transcript profiling to perform a global analysis of the biofilm-specific features of C. albicans. We produced biofilms in three laboratory models that involved very different environmental conditions. Transcript profile comparisons of these biofilms with planktonic populations led us to our initial observation that biofilm transcriptomes are highly correlated (Fig. (Fig.2A)2A) and that nutrient flow, aerobiosis, and other factors that strongly influence the transcriptome of planktonic populations do not have a noticeable effect on biofilm-growing communities. Hence, by comparison to planktonic populations, transcriptome invariance is proposed as a specific feature of biofilm populations. This observation is unlikely to reflect a slower growth of all biofilm populations than of planktonic populations, since both actively growing biofilms and stationary planktonic populations were included in the study.
The second observation is the identification of genes encoding biofilm-associated traits. As expected from the low correlation between biofilm and planktonic populations, the number of differentially expressed genes is considerable (Table (Table3).3). This is also true when both populations arise from a cph1/cph1 efg1/efg1 strain unable to develop hyphae (Table (Table3),3), indicating that the differential expression is not due only to differences in the hyphal content of planktonic and biofilm populations. We identified three significant sets of genes: the 325-gene set contains biofilm-related genes whose expression is independent of extrinsic conditions, the 317-gene set contains biofilm-related genes that are expressed independently of the intrinsic ability of the cell to form mycelia, and the 86-gene set identifies genes whose expression varies between hypha-containing and hypha-noncontaining biofilms. The classification of these genes into functional categories and the study of their representation in each set revealed that expression changes specifically affect some cellular functions (Fig. (Fig.3;3; Table Table22).
Genes involved in protein synthesis are the most highly represented in the three sets (Fig. (Fig.2;2; Table Table2),2), and they are exclusively overexpressed. One reason for this activation of the protein synthesis machinery could be that biofilms grow faster than planktonic populations. However, genes related to cell-cycle and DNA processing are underrepresented. Thus, the increase of protein synthesis is independent of cell division and unlikely to be caused by differences in growth rate. On the other hand, this excess protein might be used for the production of extracellular matrix, which is not present in planktonic cells. This is in keeping with the inclusion in the 325-gene set of several components of the secretion machinery, such as SEC14 or SYS3, and the dramatic increase in the expression of the extracellular agglutinin gene ALS1. Enhanced production and secretion of certain proteins with agglutination properties may facilitate the cohesion of cells within the biofilm. In addition to ALS1, several new genes that encode proteins with secretion signals and adhesin-like domains, such as IPF20161, IFP5185, or IPF20008, were found to be significantly overexpressed.
Among metabolic activities, amino acid biosynthetic pathways are strongly represented. We have demonstrated that GCN4 is required for efficient biofilm-like growth. In C. albicans, this transcription factor induces two different processes in response to amino acid starvation: the activation of amino acid biosynthesis and the triggering of morphogenesis (30). The latter is unlikely to be related to the biofilm-forming defect exhibited by the gcn4/gcn4 strain, since the proportion of yeast and hyphal cells was roughly the same in wild-type and gcn4/gcn4 biofilms (data not shown). Therefore, the effect of GCN4 on biofilm formation resembles a GCN-like response, although it was not observed under conditions that are known to impose amino acid starvation on planktonic populations. The situation might be different in biofilms that seem to have an increased need for protein synthesis, and amino acids might therefore become limiting faster. Under these conditions, the regulatory role of GCN4 could be decisive for biofilm progression.
A more specific role for the sulfur amino acid biosynthesis/salvage pathway is revealed by the 317- and 86-gene subsets, which identify overexpression of 9 of the 11 key genes in this route (Fig. (Fig.6).6). The route leads to the production of S-adenosylmethionine (SAM), with a separate branch for the production of cysteine (29). Since SAM is the precursor of polyamines and since ODC1 and PTK2 are also included in the 325-gene set, all of these changes could be associated with cell wall rearrangements that allow closer interaction of cells within the biofilm (13). Alternatively, biofilms may have a special requirement for methionine- and cysteine-rich proteins, an idea supported by the overexpression of genes in the cysteine biosynthesis branch. Finally, activation of the genes for SAM biosynthesis might be related to the production of a quorum-sensing molecule associated with biofilm formation. Indeed, the bacterial autoinducer AI-2 is produced from SAM, which has been proposed to serve as a “universal signal” for interspecies communication (34).
The activation of glycolytic flux and catabolic repression is underscored by the 317-gene set. In the wild-type-versus-cph1/cph1 efg1/efg1 experiment, biofilms were produced under a continuous flow of nutrients, in contrast to planktonic cultures (Table (Table1).1). Thus, activation of glycolysis could be merely a reflection of unlimited sugar availability in the microfermentor. However, activation of MIG1 in biofilms is revealed by the 325-gene set as well and is therefore independent of whether biofilms were produced under continuous flow. Overexpression of MIG1 in S. cerevisiae is associated with flocculation (28). Whether overexpression of MIG1 in C. albicans biofilms could promote flocculation and hence increase the aggregation of the cells in biofilms remains to be investigated.
Biofilm resistance to antifungals has been associated with increased expression of the MDR and CDR genes (24). Expression of these genes was studied previously in a single microtiter plate model, whereas we report a multiple-model comparison. Here, changes in the expression of these genes are not considered significant. However, CDR1 is included in the 86-gene set arising from the biofilm produced in the microfermenter model. Thus, expression of this gene could be significant under certain conditions of biofilm production, and it is also predicted to be influenced by the triggering of hyphal development concomitant with biofilm formation. Instead, we observed differential expression of ERG25 and ERG16 (in the 325-gene set) and ERG6 (in the 317-gene set). ERG16 mRNA levels correlate with azole resistance in clinical isolates C. albicans (33). An enrichment or redistribution of sterols in biofilm membranes could explain their resistance to azole-derived antifungal agents. This is consistent with the recent data of Mukherjee et al. (18) showing that expression of MDR and CDR genes in biofilms is phase specific and contributes to azole resistance only during the early phase of biofilm development whereas changes in sterol composition are involved in the resistance of the mature phase. Hence, our laboratory models for biofilm production and transcript profiling analysis converge with previous approaches to the identification of differentially expressed genes related to pathogenesis and underline the importance of the study of biofilm-like growth for a full understanding of virulence in C. albicans.
We thank J.-Y. Coppée, O. Sismeiro, M.-A. Dillies, A. Maitournam, and other members of the Institut Pasteur DNA microarray platform for their excellent assistance; F. Dromer and members of our respective laboratories for constant support; the Brown and Ernst laboratories for C. albicans strains; C. Lebos and B. Arbeille for SEM; and K. A. Charlton, C. Munro, and T. Pugsley for proofreading the manuscript. Sequence data from C. albicans were obtained from the Stanford Genome Technology Center (http://www.sequence.stanford.edu/group/candida).
Sequencing of C. albicans was accomplished with the support of the NIDR and the Burroughs Wellcome Fund. This work was supported by grants from Institut Pasteur (PTR 50), the Ministère de la Recherche et de la Technologie (P.R.F.M.M.I.P. ‘Réseau Infections Fongiques'), and the European Commission (QLK2-2000-00795; ‘Galar Fungail consortium’).