Previous Article | Next Article ![]()
Eukaryotic Cell, March 2009, p. 373-387, Vol. 8, No. 3
1535-9778/09/$08.00+0 doi:10.1128/EC.00387-08
Copyright © 2009, American Society for Microbiology. All Rights Reserved.
,
Aberdeen Fungal Group, Institute of Medical Sciences, University of Aberdeen, Aberdeen, United Kingdom
Received 12 December 2008/ Accepted 8 January 2009
|
|
|---|
mutant was mildly but significantly attenuated in terms of survival curves in the mouse model. The study has extended the range of properties known to differ between C. albicans clades and suggests a possible but minor role of phosphate metabolism in the virulence of the species. |
|
|---|
The most self-evident associations between the four major C. albicans clades and other properties are seen for geographical origins and for the presence of an intron in ribosomal DNA—the so-called ABC type (33, 34). While there is no absolute geographical association between any clade and a definable area where a clade is endemic, a high proportion of isolates within each clade usually originates from a broadly definable region (5, 38, 43, 48, 51). Moreover, there are statistically definable subsets of strains within the broader clade groupings that show more specific geographical origins (50). ABC typing of isolates shows a near-total predominance of type A isolates in clades 1 and 2 and a majority of type B isolates in clade 3 and of type C isolates in clade 4 (38).
Clade-specific associations extend to properties of potential relevance to the role of C. albicans as human commensal and pathogen. Isolates of C. albicans resistant to flucytosine are found predominantly within clade 1 (44, 51). These clade 1 isolates also show a common mechanism of flucytosine resistance, based on an amino acid change (R101C) in the FUR1 gene product, a ribosyltransferase enzyme (13, 23, 51). This mutation is not found in flucytosine-resistant isolates from other clades (13, 51). Amphotericin B resistance shows a less well pronounced but statistically significant association with C. albicans clades (4). Statistically significant differences between clades have been found for the numbers of tandem repeat sequences in alleles of the genes ALS3, ALS5, and ALS6, members of a gene family that encode C. albicans surface proteins that play a role in adhesion to host surfaces (41, 63).
The evidence clearly indicates that C. albicans major clades, delineated by MLST or Ca3 typing, comprise strains that differ in a number of properties not directly related to the particular gene sequences or tandem repeat sequences that are used to define their clade assignation. Each clade therefore comprises strains that have evolved in a manner different from strains in other clades. In this study we have prospectively tested a set of C. albicans isolates representing the four major clades for their virulence in the murine intravenous (i.v.) challenge model of C. albicans infection, for a range of phenotypic attributes, for tandem repeat alleles in several genes encoding glycosylphosphatidylinositol-anchored proteins, and for their gene expression profiles under two growth conditions.
|
|
|---|
/
) at the MAT locus as well as the more common a/
types (38), and a mix of types A, B, and C, where this was possible. The MLST data for these 43 isolates were reanalyzed in an P-distance dendrogram made by the unweighted-pair group method using arithmetic averages with phylogeny tested by 1,000 bootstraps in MEGA 3.1 software (25). For the isolates shown in Table 1, the 14 isolates in clade 1 and the 8 isolates in clade 3 coclustered in separate groups with 100% bootstrap values (Fig. 1). For the 10 isolates in clade 2 and the 11 isolates in clade 4, the bootstrap values for the cluster nodes were 90% and 86%, respectively (Fig. 1). |
View this table: [in a new window] |
TABLE 1. Details of the 43 C. albicans isolates included in this study
|
![]() View larger version (11K): [in a new window] |
FIG. 1. Dendrogram showing the similarity clusters (referred to as clades [shown to the right of illustration]) for the 43 C. albicans clinical isolates used in this study. The dendrogram was constructed by the unweighted-pair group method using arithmetic averages with MEGA 3.0 software as previously described (50). The numbers at nodal points indicate bootstrap values (as percentages) for 1,000 replications. The scale bar at the bottom of the figure shows nucleotide substitutions.
|
Deletion of PHO100 from C. albicans SC5314 and construction of reintegrant strains.
The CaPHO100 gene was disrupted by the standard "ura-blaster" protocol (ura, uracil) (15). The disruption cassette was produced by amplifying two regions of the gene by PCR; a 5' region, containing the start codon, was amplified with primer pair DMPHO100_1 and DMPHO100_2 (Table 2), and a 3' region, containing the stop codon, was amplified with primer pair DMPHO100_3 and DMPHO100_4 (Table 2). Restriction digests of the 5' region by SacI and BglII and the 3' region with SalI and SphI allowed the DNA fragments to be cloned into digested pMB7 (15). Digestion of the resulting plasmid with SacI and SphI released the disruption cassette. The disruption cassette consisted of the ura-blaster cassette flanked by complementary sequences of 340 bp at the 5' end and 460 bp at the 3' end of PHO100. PHO100 was sequentially disrupted by transformation of strain CAI-4 (15), with the Ura3 marker recycled by selection on SD medium (0.67% [wt/vol] yeast nitrogen base with ammonium sulfate without amino acids but with 2% glucose) containing 5-fluoroorotic acid (1 mg/ml) and uridine (50 µg/ml). To overcome any problems with URA3 expression from the PHO100 locus (6), the Ura-minus null mutant (DMY001; pho100
::hisG/pho100
::hisG) was transformed with StuI-digested CIp10 plasmid (35), which leads to URA3 expression at the RPS1 locus (DMY002; as DMY001 but RPS1/rps1
::CIp10). Correct disruption of the PHO100 alleles was analyzed by PCR (primer pair DMPHO100_1 and DMPHO100_5) and confirmed by Southern analysis. A reintegrant strain for use as a control was created by PCR amplifying the entire PHO100 ORF plus 885-bp promoter and 454-bp terminator sequences (primer pair DMORF_1 and DMORF_2), and cloning it into the pGEM-T vector (Promega Ltd., Southampton, United Kingdom). The insert was released by NotI digestion and cloned into the NotI site in CIp10. The resulting plasmid was digested with StuI and used to transform the Ura-minus pho100
null strain (DMY003; as DMY001 but RPS1/rps1
::CIp10-PHO100). The control strain for these experiments was NGY152 (6), which is strain CAI-4 transformed with StuI-digested CIp10.
|
View this table: [in a new window] |
TABLE 2. Oligonucleotide primers used for gene disruption and screening of C. albicans PHO100
|
The mice were observed for up to 28 days postchallenge and were humanely terminated when their body weight fell more than 20% below their initial weight or when the animals showed signs of serious illness. The dates of termination of animals were used to determine survival after challenge. At the time of death, both kidneys were dissected with aseptic precautions and homogenized together in 0.5-ml volumes of sterile water, and tissue burdens were measured by counting the viable cells on Sabouraud agar. The percent change in body weight of each mouse between the day of challenge and 3 days postchallenge was calculated relative to the weight of the mouse on the day of challenge.
The same model was used to test virulence of the pho100
mutant, but a single experiment was done with groups of six mice challenged i.v. with control and mutant strains.
Measurement of phenotypic properties. Several phenotypic properties were characterized for all 43 C. albicans isolates. For all phenotypic tests, cultures of each isolate were grown in 5 ml of NGY medium overnight at 30°C with constant rotation. Cells were harvested and washed with sterile water prior to use in phenotypic tests. Where required, cell numbers were determined by counting cells with a hemocytometer.
Biofilm formation was measured by resuspending cells in RPMI 1640 medium containing 10% (vol/vol) serum (RPMI 1640-10% serum) (1x RPMI 1640 medium, 2 mM L-glutamine, 1.65 mM morpholinepropanesulfonic acid [MOPS], 10% fetal calf serum) to a density of 2 x 106 cells/ml in triplicate wells of a 24-well tissue culture plate. Plates were incubated overnight at 37°C in 5% CO2 without agitation. Biofilm formation was gauged by reading the optical density at 600 nm (OD600) of the supernatant after agitation.
Germ tube formation and mean morphology index were determined for cells incubated for 3 h in fetal calf serum. Washed cells (3 x 106) were resuspended in 1 ml of 100% fetal calf serum (prewarmed to 37°C) and incubated in a 37°C water bath for 3 h. The cultures were fixed with formalin, and the morphology of cells was determined microscopically. For each strain, the morphology index was determined on three separate occasions.
For measurement of growth under proteinase-inducing conditions, overnight cultures in NGY medium were washed twice in sterile water and reinoculated into 5-ml volumes of 1.17% yeast carbon base (Difco)—1% glucose and 0.5% bovine serum albumin (Sigma) at an initial concentration of 104 yeasts/ml. Cultures were incubated for 72 h at 30°C with constant rotation. Turbidity ODs were measured versus a medium blank at 600 nm.
Acid phosphatase activity was determined in intact C. albicans cells. These cells were grown in phosphate-free NGY medium for 16 h, washed once in water, and resuspended in water to one-fifth of the original culture volume. Enzyme activity was measured by the following modification of the method of Lindhardt and Walter (29). The final reaction mixture contained 50 µl of 5.5 mM p-nitrophenyl phosphate (ICN Biomedicals Inc., United Kingdom) in 0.11 M citrate buffer (pH 5.0) and 5 µl cell suspension. The reaction in control wells was immediately stopped by the addition of 245 µl of 0.1 M NaOH. The plate was incubated for 15 min at 37°C, and then the reaction was stopped with the addition of NaOH to the remaining wells. The absorbance of the yellow p-nitrophenyl released was measured at 405 nm, and the amount of p-nitrophenol released was calculated from a standard curve. Experimental measurements were corrected by subtraction of control measurements. Specific activity was defined as nanomoles of p-nitrophenol/15 min/OD600 in the cell suspension.
Alcian blue binding to surface charged groups was measured for yeast cells grown in YPD (1% yeast extract, 2% mycological peptone, 2% glucose) at 30°C for 18 h with constant rotation according to the method of Hobson et al. (21). For each assay, 50 µl of cell suspension (
1.5 x 107 cells) was harvested, washed in water, and resuspended in 1 ml of 30 µg/ml alcian blue in 20 mM HCl. The cell suspensions were incubated for 10 min at room temperature, and then the cells spun down. A620 values were determined for the supernatants, allowing the amount of alcian blue bound by the cells in the assay to be determined from a standard curve. Specific alcian blue binding was expressed as micrograms of alcian blue bound per OD600 unit of cells.
The growth rates of cells in RPMI 1640 medium were determined by adding 10 µl of washed C. albicans cells to 5 ml RPMI 1640 with low MOPS (1x RPMI 1640, 2 mM L-glutamine, 1.65 mM MOPS [pH 6.5]). For each strain, 200 µl of cell suspension was dispensed into triplicate wells of 96-well flat-bottomed tissue culture plates. Growth curves were plotted at 40°C over 16 h by means of a Fluostar Optima plate reader (BMG Labtech Ltd., Aylesbury, United Kingdom). Growth rates were calculated from the exponential region of the growth curve and were recorded as means from the growth rates determined on three separate occasions.
To measure NaCl tolerance, 10 µl of washed C. albicans cells was added to 5 ml 2x RPMI 1640 with low MOPS (1x RPMI 1640, 2 mM L-glutamine, 1.65 mM MOPS [pH 6.5]), and 100-µl lots were dispensed into wells containing 100 µl 4 M NaCl, such that tolerance to 2 M NaCl was assayed. For each C. albicans isolate, triplicate wells were assayed on at least two occasions.
Adherence of fungi to BECs and catheter material. Buccal epithelial cells (BECs) were collected by scraping the gums and cheeks of healthy donors, washed, and resuspended (3 x 105 cells/ml) in physiological saline. Cells were monitored microscopically for prior colonization by yeast cells. Washed C. albicans cell suspensions were adjusted to 5 x 106 cells/ml in saline. Adherence assay solutions (600 µl) contained 6 x 104 BECs mixed with 1 x 106 C. albicans cells. The solutions were incubated at 30°C with constant rotation for 1 h. The number of C. albicans cells adherent to each BEC was determined microscopically, with at least 100 BECs examined for each assay. Assays were performed for each strain on at least three separate occasions.
Adherence to catheter material was measured by the method of Seidler et al. (47). Arrow-Howes Hands-Off multilumen central venous catheters (radiopaque polyurethane) were a kind gift from Nigel Webster, University of Aberdeen. Briefly, catheters were cut into 5-mm sections and weighed. Catheter pieces were preconditioned in 100% fetal calf serum at 37°C with agitation overnight. Catheter material was washed twice in phosphate-buffered saline and then immersed in 1 ml C. albicans cell suspension (2 x 106 cells/ml in RPMI 1640-10% serum), and the cell suspension containing catheter material was incubated for 2.5 h at 37°C with shaking to allow adherence. Sections were washed three times with phosphate-buffered saline and then added to 1 ml RPMI 1640-10% serum, and the catheter culture was incubated for 24 h at 37°C under 5% CO2. At this time, cell growth in catheter culture was estimated using Alamar Blue (AbD Serotec, Oxford, United Kingdom) by measuring the percent reduction Alamar Blue per gram catheter plastic.
Expression profiling. From the full panel of 43 C. albicans isolates, a subset of 20 was chosen for inclusion in expression profiling experiments. The subset comprised five isolates from each of the four clades, selected to represent as much diversity as possible with respect to mouse virulence, MAT homozygosity, and ABC type (Table 1).
For expression profiling, cells were grown either in NGY medium at 30°C or in RPMI 1640-10% serum at 37°C. For NGY-grown cells, an NGY preculture was used to inoculate 100 ml NGY at 30°C to an initial OD600 of approximately 0.15. Cells were grown with constant shaking at 200 rpm until they reached an OD600 of approximately 0.8. For cells grown in RPMI 1640-10% serum, 1 ml of an NGY preculture was harvested and the cells were washed and resuspended in 1 ml distilled H2O. This was used to inoculate 120 ml RPMI 1640-10% serum at 37°C in a tissue culture flask. Cultures were incubated at 37°C under 5% CO2 for 16 h without agitation. Cells were scraped off the flask surface and resuspended in medium. For both growth conditions, cells were harvested by centrifugation, resuspended in residual liquid, and flash-frozen in liquid nitrogen. RNA was extracted by the method of Hauser et al. (20). The quality of the RNA was analyzed by agarose gel electrophoresis, and RNA concentration was measured by spectrophotometry. For each isolate under each growth condition, three independent cultures were grown for RNA extraction.
For cells grown in both NGY medium and RPMI 1640-10% serum, SC5314 RNA samples were used to prepare a pooled control, which consisted of equal quantities of RNA from each of the three samples processed. The pooled control was used for comparison with each isolate, including SC5314. Labeled cDNAs were prepared (11), with 25 µg isolate RNA labeled with Cy3 and 25 µg pooled control labeled with Cy5. For each isolate and growth condition, three microarrays were carried out comparing the three independent RNAs extracted to the relevant pooled control. Microarray hybridizations, washes, and scanning and acquisition of data were performed as described previously (11).
The results of the expression profiling exercise were two sets of triplicate data for 20 C. albicans isolates, each grown in NGY and RPMI 1640 media and profiled with reference to pooled RNA from C. albicans SC5314 triplicate cultures grown under the same conditions. The data were normalized on a per spot and per chip intensity-dependent (Lowess) basis and on a per chip basis to the 50th percentile, and then the expression ratios were converted to base 2 logarithms. The Statistical Analysis of Microarrays (SAM) (54) software package (http://www-stat.stanford.edu/
tibs/SAM/) was used to determine significantly differentiated genes between groups of isolates. The analyses were done with expression profile data grouped in three strata according to overall virulence and in four strata by clade; multiclass SAM comparisons were done between strata. For each analysis, two approaches were used. In the first, the triplicate data were analyzed directly. In the second, a proprietary computer script was used to exclude the most extreme outlier among each set of triplicates, and SAM analysis was then applied to the most concordant pair from the three data for each gene and isolate. Data sets used in the SAM analysis are provided as tables in the supplemental material. Genes were considered to contribute significantly to differentiation between clades and virulence groups when their q value by SAM was <0.05. Only genes that emerged as significant multiclass differentiators from both the triplicate and the "best two of three" data sets were recorded as likely to represent differences in constitutive gene expression between the strata analyzed. Transcript profiling data can be accessed via ArrayExpress (http://www.ebi.ac.uk/microarray): identification numbers E-MEXP-1933 (RPMI 1640 data) and E-MEXP-1931 (NGY data).
ALS and HYR tandem repeat determinations. For approximation of the number of tandem repeats for ALS and HYR genes, genomic DNA was extracted from each of the 43 C. albicans isolates (22). PCR mixtures (50 µl) containing 100 ng genomic DNA, 0.4 µM forward and reverse primers, and 1x PCR Mastermix (Thermoscientific, Epsom, Surrey, United Kingdom) were set up for each isolate. Oligonucleotide primer sequences are listed in Table 3. PCR was performed as follows: one step of 1 min at 94°C, followed by 30 cycles (1 cycle consisting of 1 min at 94°C, 1 min at 50°C, and 3 min at 72°C), with a final step of 10 min at 72°C. PCR products were analyzed by agarose gel electrophoresis against appropriate molecular weight DNA markers (Promega UK, Southampton, United Kingdom).
|
View this table: [in a new window] |
TABLE 3. Oligonucleotide primers for analysis of cell wall gene repeat regions
|
For analysis of associations between tandem repeat data and clades, two approaches were used. In the first, results for each isolate were represented as the average number of repeats between the diploid alleles for each isolate; in the second, each allele was analyzed individually, so each isolate was represented in duplicate. The Kruskal-Wallis test was applied in both cases.
Differences in mouse survival were tested by Kaplan-Meier/log rank analyses, and differences in body weight changes and mouse kidney burdens were tested by the Mann-Whitney U test.
|
|
|---|
![]() View larger version (12K): [in a new window] |
FIG. 2. Scatterplot of mean change in body weight over 3 days postchallenge versus mean survival time for mice infected with 43 C. albicans isolates. The regression line is shown. The Pearson correlation coefficient r for the data was 0.79; Kendall tau (nonparametric correlation coefficient) was 0.65. For both measures of correlation, P was <0.001.
|
|
View this table: [in a new window] |
TABLE 4. Virulence of 43 C. albicans isolates in mice and other measured characteristics of the isolates tested in vitro
|
The three variables found to show statistically significant associations with the overall virulence ranking also emerged as significant factors in multiple regression analysis of individual virulence parameters. For MST as the dependent variable, acid phosphatase activity (P = 0.007) and growth at 40°C (P = 0.005) were the only two significantly associated variables in a two-step regression analysis. With weight change as the dependent variable, acid phosphatase activity (P = 0.003) and growth at 40°C (P = 0.003) again emerged as the most significant among the four originally most significant variables, although in this analysis, the clade (P = 0.025) also showed a possibly significant association. Finally, the analysis with kidney burden as the dependent variable showed only MAT heterozygosity (P = 0.002) as a significantly associated factor.
Associations between the clade and other properties of C. albicans isolates. Kruskal-Wallis analysis of the data (Table 4) for tests in vitro with the C. albicans clade as the grouping variable showed significant associations between the clade and acid phosphatase activity (P = 0.002), ABC type (P < 0.001), and growth in 2 M NaCl (P < 0.001). The percentage of BECs with C. albicans attached (but not the average numbers of C. albicans attached to BECs) showed a possibly significant association with clade (P = 0.033). Only 4 of the 43 isolates made poor or no biofilms on polystyrene under the conditions of the experiment: all were from clade 4, but this association was not significant in the Kruskal-Wallis analysis (P = 0.056). Table 5 summarizes the results of the tests done in vitro, stratified by clade. Data on MAT heterozygosity were not included, since the proportion of homozygous and heterozygous isolates depended on the choice of isolates. The data on ABC type were included as a positive control for validation of the statistical analysis. ABC types A, A, and B are already known to correlate strongly with clades 1, 2, and 3, respectively (38). The association of acid phosphatase activity with clades appeared largely as a generally lower level of activity among clade 2 isolates (Table 5). For clade 2 isolates, mean acid phosphatase specific activity was 249 ± 41 units compared with 735 ± 192, 603 ± 79, and 640 ± 113 units for isolates in clades 1, 3, and 4, respectively. The association of growth in 2 M NaCl with clades was largely the consequence of a much higher proportion of clade 1 isolates that grew at this salt concentration compared with the other clades (Table 5).
|
View this table: [in a new window] |
TABLE 5. Summary of results for strain properties and tests in vitro for isolates in each of the four main C. albicans clades
|
![]() View larger version (17K): [in a new window] |
FIG. 3. Alleles of midrepeat sequences for genes in the ALS and HYR families. For each of 43 C. albicans isolates, the number of sequence repeats (shown on the y axis) is indicated by a line joining the repeat numbers for the two alleles of each gene. A horizontal line thus indicates homozygous repeats of the numbers indicated in the y axis. A diagonal line indicates an isolate with heterozygous repeats, with the larger of the two allele repeats indicated by the left-hand position of the line. For isolates with identical allele pairs, the results have been jittered vertically to make each isolate line visible. The isolate data are grouped in columns representing results for isolates in clades 1, 2, 3, and 4. Note that some y axes have been compressed to emphasize interclade differences.
|
|
View this table: [in a new window] |
TABLE 6. Summary statistics of ALS and HYR tandem repeat length estimates for 43 C. albicans isolates
|
Expression profiling of 20 C. albicans isolates from clades 1 to 4. Five isolates selected from each of the four main C. albicans clades were grown to late-exponential phase in NGY and RPMI 1640 media. NGY is the medium used routinely in this laboratory for growth of challenge inocula in the mouse model of disseminated C. albicans infection; it supports the growth of C. albicans purely in the yeast form. It was chosen since the profile observed in this medium is that of cells immediately before they infect experimental animals. RPMI 1640 is a tissue culture medium that mimics the composition of human body fluids and is utilized in applications such as yeast susceptibility testing (37). C. albicans cells grow as mixtures of pseudohyphal and hyphal forms in this medium (7). For both media, expression profiles for the 20 isolates were normalized against the profile for pooled RNA from strain SC5314 grown under the same conditions. The quality of the expression profile data is illustrated in Fig. 4, which shows interpretation graphs for two randomly chosen isolates each in clades 2 and 3 relative to the reference isolate SC5314. Here the relative expression levels of all genes are depicted for 4 of the 20 profiled isolates and the reference isolate SC5314. Figure 4 shows that the replicate quality of the triplicate data was good, and the triplicate SC5314 profiles showed only minor expression variations from the SC5314 RNA pool used to normalize all the data. Figure 4 also illustrates the considerable interisolate diversity found in constitutive gene expression in the present experiments. The expression levels depicted in Fig. 4 suggested extensive interclade gene expression differences, but this level of differential expression was not sustained by analysis of expression profiles for all 20 isolates.
![]() View larger version (66K): [in a new window] |
FIG. 4. Interisolate differences in constitutive gene expression profiles. The figure depicts the GeneSpring interpretation graph for triplicate expression profiles of two C. albicans isolates each from clades 2 and 3, plus reference isolate SC5314, grown in NGY medium. Each line in the graph indicates the relative expression level for a single gene. The expression level of the gene in strain AM2003/0100 compared to SC5314 sets the color of the line, which is kept through all the other strains.
|
![]() View larger version (43K): [in a new window] |
FIG. 5. Heat maps indicating the average differences in RNA levels between genes found to be significant interclade differentiators (q < 0.05 by SAM analysis) for 20 C. albicans isolates grown in NGY medium (a) and in RPMI 1640 medium (b). The heat map colors are based on pale yellow indicating no upregulation relative to the pooled SC5314 RNA control. Colors ranging from purple (maximum of 20-fold) through red and orange indicate RNA levels greater than the level for the control, while green through pale blue (20-fold downregulation) colors indicate RNA levels below the level for the control. In both figures, the heat map for SC5314 (mean of triplicate data versus SC5314 RNA pool) is shown in the leftmost column.
|
SAM analysis of the expression profiling data stratified by mouse virulence phenotype (Fig. 6) determined significant multiclass differentiation for 21 ORFs (22 microarray spots) in NGY-grown cells. On the list of differentially expressed genes between the three groups (based on overall virulence ranking; Table 4), genes with reference numbers CA0636 (FDH4.3F), CA0639 (FDH13.3F), and CA3258 (IFC2) are no longer accepted in the current assembly (assembly 21) of the C. albicans genome database (CGD) (http://www.candidagenome.org/). Our finding of differential expression of these ORFs suggest that they may in fact be valid for inclusion in the database. Figure 6 lists the 18 significant and currently accepted genes (one significant gene was represented twice on the arrays, so that the list contains 19 entries) with multiclass differential expression based on the overall virulence of the 20 C. albicans isolates. For RPMI 1640-grown cells, only two genes (ORFs 19.6847 and 19.3075; Fig. 6) met the criteria for inclusion as differentially expressed between isolates in virulence groups 1, 2, and 3.
![]() View larger version (56K): [in a new window] |
FIG. 6. Heat map indicating the average differences in RNA levels between genes found to be significant multiclass differentiators between C. albicans isolates of high (group 1), intermediate (group 2) and low (group 3) virulence in the mouse model (see Table 4). The heat map was set up as explained in the legend to Fig. 5. CA ref., CA reference.
|
Characterization of a pho100
mutant in vitro and in vivo.
Because the expression profiles possibly associated with virulence contained three upregulated genes associated with phosphate uptake or metabolism, we elected to delete PHO100, which encodes a predicted periplasmic acid phosphatase enzyme (http://www.candidagenome.org). We reasoned that, in this surface location, the phosphatase may be directly involved in host-C. albicans interactions. Yeast cells of the pho100
mutant showed elevated acid phosphatase activity in vitro (specific activities of the mutant and pho100::rps1
::CIp10-PHO100 reintegrant strains were, respectively, 219% [±29%] and 114% [±37%]) relative to the control strain CAI-4 carrying plasmid CIp10.
In the mouse intravenous challenge model, the pho100
mutant showed modest but statistically significant attenuation of virulence in terms of time to euthanasia (Fig. 7) relative both to the CAI-4 control strain with CIp10 (P < 0.01) and to the pho100::PHO100 reintegrant (P < 0.05). There were no statistically significant differences in percent weight changes between day 0 and day 3 (means ± standard deviations [SDs] were –8.5% ± 4.8%, –13.9% ± 3.1%, and –10.1% ± 3.1%, respectively, for the mutant, CAI-4 control strain with CIp10, and reintegrant strain). Similarly, no significant differences were found between the three strains in terms of kidney fungal burdens at the time of demise. (The log10 CFU/g tissue data were as follows: 5.5 ± 1.6 for the null mutant, 6.0 ± 1.0 for the CAI-4 control strain with CIp10, and 5.8 ± 0.4 for the reintegrant.)
![]() View larger version (16K): [in a new window] |
FIG. 7. Time to euthanasia curves for groups of six mice infected i.v. with the control strain (NGY152; solid circles), pho100::pho100 mutant strain (DMY002; squares), and pho100::PHO100 reintegrant strain (DMY003; open circles).
|
|
|
|---|
The basis for clade-level differentiation of C. albicans strains by MLST is partly arbitrary (40) but reflects at least well-established geographical associations suggestive of an independent evolutionary history for isolates in clades 1 to 4 (40, 48, 51). Isolates from the four major clades should therefore exhibit differences detectable by fine-scale measurement but are unlikely to show gross phenotypic differences. In common with previous authors (41, 63), we found statistically significant differences in the numbers of midrepeat sequences for ALS3, ALS5, and ALS6, and we extended the list of genes with clade-associated variable midrepeat sequences to all the other members of the ALS gene family except ALS1 (Table 6) and to three midrepeat sequences in the genes HYR1 and HYR3 (Table 6). Since most of the products encoded by these gene families are surface proteins with known contributions to adherence of C. albicans to various surfaces (58-62) or to yeast-hypha morphogenesis (3, 8), the variations in midrepeat allele lengths between clades might be assumed to contribute to overall virulence of isolates. However, no associations were found between the various midrepeat alleles and strain virulence for mice. Since adherence and morphogenesis are phenomena mediated and regulated by many gene products (1, 9, 14, 19, 28, 42), it is likely that C. albicans possesses considerable functional redundancy with respect to these properties, so while differences may arise between adhesin structures within different strain clades, these are unlikely to translate to detectable effects on strain virulence.
Our expression profiling data generated lists of genes with statistically significant associations with strain clades. Microarray profiling of a complement of more than 6,000 genes with RNA extracted in triplicate from 20 isolates of a fungus might easily result in a set of apparently differential genes arising by chance. Even with appropriate statistical analysis, the sheer number of data points involved offers considerable opportunity for artifactual significant differences in gene expression; hence, we used a highly rigorous approach to data analysis. We compared the SAM outputs from analyses of all triplicate data for each isolate with the outputs for analyses done with the closest pair of results for each isolate-gene combination and accepted from the outputs only genes common to both analyses as potentially significant differentiators at the level of clade or virulence group. One hypothetical expectation for genes with expression differing between clades of strains is that at least some of such genes would encode enzymes in a common biochemical pathway. Our data provide no suggestion of pathway links between any of the genes emerging as differentially expressed between different clades (Fig. 5); the significant multiclass genes listed are all effectively singletons. We have reanalyzed our expression data with the software package GALGO (53), which combines Bayesian and multiple regression approaches to detect sets of genes among expression data that are differentially associated with grouping strata. The software generated a list of genes common with those in Fig. 5 but did not detect any gene combinations that related to strain clades (unpublished data). It is clear from Fig. 5 that expression profiling of none of the genes listed as significant in these experiments could be reliably used to differentiate new isolates at the clade level. However, it is unquestionable that individual isolates of C. albicans can differ markedly in their gene expression profiles when grown under identical conditions (Fig. 4). We cannot exclude the possibility that genes apparently expressed at lower levels in some strains than others may have DNA sequence differences from those of the oligonucleotides in the microarray. Resequencing of a well-selected set of isolates could address this possibility. Overall, we interpret our expression profiles in Fig. 5 as further indication of clade-associated phenomena that represent statistical rather than absolute interclade differences. Profiling for isolates grown under different conditions may, of course, reveal other gene sets that differ between clades.
Among our set of phenotypic properties, levels of acid phosphatase activity and ability to grow in 2 M NaCl were the only parameters significantly associated with strain clades. Average levels of acid phosphatase activity tended to be lower among clade 2 isolates than among isolates from other clades, and all but one of the clade 1 isolates grew in 2 M NaCl compared with 40% in clade 2 and fewer than 20% in clades 3 and 4. Acid phosphatase activity also emerged as a property associated with overall strain virulence (discussed below), but clades themselves were not significantly associated with overall virulence, so we conclude that acid phosphatase should be regarded as a weakly differential interclade property. The finding of generally higher salt tolerance among our clade 1 isolates is similarly difficult to interpret. Clade 1 isolates represent the most commonly encountered C. albicans strain types on a world-wide scale. There have been previous suggestions that clade 1 isolates may be slightly better adapted than others to colonize or superficially infect human hosts (17, 38, 40, 46), but we are not able to suggest how slightly enhanced salt tolerance might contribute to this situation.
Considered overall, the results of this study relating a diverse range of properties to C. albicans clades indicate that direct measurements of genetic diversity, such as determinations of midrepeat sequence lengths, are more successful than transcriptome data and phenotypic testing at revealing interclade differences. Clade-level phenotypic properties are apparent, notably the previously described association of clade 1 isolates with flucytosine resistance and the present finding of much greater salt tolerance among clade 1 isolates. However, none of these represents a characteristic unique to any clade, nor have there been any demonstrable inherent differences in the overall murine virulence of strains from different clades, even in this, the largest examination of interstrain differences in virulence in mice.
Our data relating properties of individual isolates with virulence in mice show more robust associations than the analyses of clade associations. Lower virulence in mice among MAT-homozygous isolates compared with their MAT-heterozygous counterparts has been demonstrated previously. In these tests small numbers of isolates were compared with laboratory-generated MAT-homozygous mutants alone and in experiments with competitive challenge designs. Survival was used as the virulence endpoint (24, 30). Other authors tested a larger panel of clinical isolates and found a lower virulence among MAT-homozygous isolates (57). Our data, based on virulence testing with 43 clinical isolates of C. albicans, fully confirm these previous findings, since no MAT-homozygous isolates were among the 18 most virulent we tested (Table 4). In common with Wu et al. (57), we were unable to confirm a prior suggestion that MAT
/
strains might be even less virulent than a/a types (24). Our results show a mixed distribution of
/
and a/a isolates among the ranks of isolates with intermediate and low virulence in mice (Table 4). The lower virulence of MAT-homozygous isolates is probably a consequence of homozygosity in other genes adjacent to the MAT locus rather than a specific MAT-related event (57).
Phosphate uptake and metabolism emerge from our data as properties associated with mouse virulence in C. albicans. This conclusion is based on two lines of evidence from this study. First, a statistically significant association was found between acid phosphatase activity in vitro and gross virulence of isolates (Table 4). None of the 14 most virulent isolates had acid phosphatase specific activities greater than 700 nmol/15 min/OD600, compared with 7/14 of the least virulent isolates. Second, expression of four genes encoding phosphate-specific activities was shown as significantly related to mouse virulence (Fig. 6): PHO84, PHO87, PHO100, and VTC4, which encode an enzyme involved in phosphate transport, an enzyme involved in phosphate uptake, a phosphomonoesterase, and a polyphosphate synthase, respectively. Acid phosphatase activity detectable in intact C. albicans cells was proposed many years ago as one of several hydrolytic enzymes with theoretical potential as a virulence factor (10). An inducible phosphomonoesterase, which can be released from intact cells by treatment with sulfhydryl reagents (18), was one of the first biological examples of a highly mannosylated protein with enzyme activity (39). In assembly 21 of the C. albicans genome database (http://www.candidagenome.org), five acid phosphatase genes are listed: PHO100, PHO111, PHO112, PHO113, and PHO114. Of these, only PHO100 is inducible and predicted to be located on the surface: we presume this encodes the inducible surface enzyme that has been previously purified and characterized (39) and listed as upregulated in expression profiles of biofilm-forming C. albicans (36). From the C. albicans genome database and our own literature searches, PHO100 has not emerged as a significantly differentially regulated gene in any other studies concerned with C. albicans and virulence. Deletion of the PHO100 gene in the SC5314 background generated a mutant with paradoxically elevated acid phosphatase activity. The mutant showed mild but significant attenuation in virulence in mice as measured by time to euthanasia, and animals infected with the mutant had smaller changes in body weight over the first 3 days postchallenge and the lowest mean kidney fungal burdens compared to the parental reference and PHO100 reintegrant strains, but these differences were not statistically significantly different. The increase in acid phosphatase activity in the pho100 mutant is compatible with the finding of higher acid phosphatase levels in the clinical isolates with the least virulence in mice (Table 4). We think that more-detailed investigation of phosphatase activity is warranted in the context of pathogenesis of C. albicans infections.
The finding of a significant association between gross virulence in mice and an isolate's growth rate in a tissue culture medium at 40°C but not at lower temperatures (Table 4) is compatible with the clinical observation of fever in patients with disseminated C. albicans infections. However, fever is not a characteristic of the mouse i.v. challenge model; body temperature tends to fall with progressive disease (49, 56). We are therefore unable to explain the significance of the high temperature growth rate association with virulence. Monitoring of body temperature to predict the onset of C. albicans terminal disease has been recommended in the mouse model (56). Our own data suggest that mouse body weight change over 3 days postchallenge is also a useful predictor of survival, since the two parameters showed a highly significant correlation. We are currently developing a virulence score based on day 3 weight changes and kidney burdens that will eliminate the need for prolonged survival experiments to determine the virulence of C. albicans strains.
|
|
|---|
Published ahead of print on 16 January 2009. ![]()
Supplemental material for this article may be found at http://ec.asm.org/. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»