Previous Article | Next Article ![]()
Eukaryotic Cell, June 2009, p. 806-820, Vol. 8, No. 6
1535-9778/09/$08.00+0 doi:10.1128/EC.00002-09
Copyright © 2009, American Society for Microbiology. All Rights Reserved.
,
Teresa T. Liu,2,3
Geneviève Boucher,1
P. David Rogers,2,3 and
Martine Raymond1,4*
Institute for Research in Immunology and Cancer, Université de Montréal, Montreal, Quebec, Canada H3T 1J4,1 Departments of Clinical Pharmacy, Pharmaceutical Sciences, Molecular Sciences, and Pediatrics, University of Tennessee Health Science Center, Memphis, Tennessee 38163,2 Children's Foundation Research Center of Memphis, Le Bonheur Children's Medical Center, Memphis, Tennessee 38103,3 Department of Biochemistry, Université de Montréal, Montreal, Quebec, Canada H3T 1J44
Received 2 January 2009/ Accepted 14 April 2009
|
|
|---|
0.01). Strikingly, Cap1p binding was detected not only at the promoter region of its target genes but also at their 3' ends and within their open reading frames, suggesting that Cap1p may associate with the transcriptional or chromatin remodeling machinery to exert its activity. Overrepresented functional groups of the Cap1p targets (P
0.02) included 11 genes involved in the OSR (CAP1, GLR1, TRX1, SOD1, CAT1, and others), 13 genes involved in response to drugs (PDR16, MDR1, FLU1, YCF1, FCR1, and others), 4 genes involved in phospholipid transport (PDR16, GIT1, RTA2, and orf19.932), and 3 genes involved in the regulation of nitrogen utilization (GST3, orf19.2693, and orf19.3121), suggesting that Cap1p has other cellular functions in addition to the OSR. Bioinformatic analyses of the bound sequences suggest that Cap1p recognizes the DNA motif 5'-MTKASTMA. Finally, transcriptome analyses showed that increased expression generally accompanies Cap1p binding at its targets, indicating that Cap1p functions as a transcriptional activator. |
|
|---|
Azoles, including both imidazoles (e.g., ketoconazole) and triazoles (e.g., fluconazole [FLC]), inhibit the function of the lanosterol demethylase enzyme Erg11p, a component of the ergosterol biosynthesis pathway, leading to methylsterol accumulation, sterol depletion, and consequently to growth arrest (1). This fungistatic property of azoles coupled to their repeated use in the clinic renders the surviving C. albicans cells prone to the selection of mutations conferring azole resistance. The clinical resistance of Candida spp. to azole is a challenging problem for clinicians. Azole resistance develops particularly in human immunodeficiency virus-infected patients with recurrent episodes of oropharyngeal or esophageal candidiases (18). The molecular mechanisms of clinical resistance to azole in C. albicans involve (i) mutations in the target of azoles, Erg11p, resulting in altered drug binding, and/or (ii) the constitutive overexpression of genes responsible for the drug resistance phenotype, including CDR1 and CDR2, which encode ATP binding cassette transporters; MDR1, which encodes a transporter of the major facilitator superfamily; PDR16, which codes for a phospholipid transferase; and/or ERG11 (1, 40, 46, 55, 56). Recent studies revealed the direct involvement in clinical azole resistance of gain-of-function mutations in genes encoding transcription factors of the fungus-specific zinc cluster family (42). It was shown that activating mutations in the transcription factor Tac1p (for transcriptional activator of CDR genes) leads to the constitutive overexpression of its target genes, CDR1, CDR2, and PDR16, in clinical isolates of C. albicans (15-17, 38, 70). Similarly, MDR1 constitutive overexpression is due to gain-of-function mutations in the zinc cluster transcription factor Mrr1p (for multidrug resistance regulator) (19, 47). Finally, an activating mutation in the transcription factor Upc2p (for uptake control) was shown to be responsible for clinical azole resistance and the upregulation of ERG11 and, to a lesser extent, MDR1 (20).
Genome-wide location and/or expression studies have shown that several targets of the transcription factors Tac1p, Upc2p, and Mrr1p have established or predicted roles in the oxidative stress response (OSR) (20, 38, 47, 71), suggesting that multidrug resistance and the OSR are interconnected processes in C. albicans. Studies of the OSR in Saccharomyces cerevisiae have shown that Yap1p, a basic region leucine zipper (bZIP) transcription factor homologous to mammalian activating protein 1 (AP-1), is a key regulator of this process (27, 48). Yap1p binds to Yap1 recognition elements (YRE, for Yap1 recognition elements) TTA(C/G)T(A/C)A, located in the promoter of its target genes (49), and controls the expression of genes encoding the majority of antioxidants and thiol-oxidoreductases, such as the glutathione reductase GLR1 and the thioredoxin reductase TRR1 (27, 48). Yap1p is also essential for the response to cadmium or drug exposures and can be activated by chemicals (e.g., diamide), antifungal agents (e.g., benomyl), ionizing radiation, or toxic endogenous cellular metabolites (e.g., methylglyoxal) (2, 27, 43, 45, 48, 49). Interestingly, Yap1p confers azole resistance in S. cerevisiae by activating the expression of FLR1, the functional homolog of C. albicans MDR1 (2). Yap1p is activated by a mechanism acting on its nuclear export (33, 34). In response to high levels of oxidants, Yap1p undergoes redox conformational changes caused by intramolecular bond formation between cysteine residues within the C-terminal cysteine-rich domain (CRD) of the protein (33, 34). This prevents the interaction of Yap1p with the nuclear exportin Crm1p, leading to nuclear retention and transcriptional activation (27, 35, 48, 67).
The C. albicans CAP1 gene encodes the functional homolog of Yap1p and has been isolated based on its ability to confer azole resistance when expressed in S. cerevisiae (2, 3, 68). The mechanisms whereby Cap1p exerts its function are reminiscent of S. cerevisiae Yap1p, as a truncation of the Cap1p CRD (CAP1-TR) or mutagenesis of the third cysteine residue of the Cap1p CRD (CAP1-C477A) results in enhanced resistance to toxic compounds, including azoles, the heavy metal cadmium, and the oxidative stress-inducing agent 4-nitroquinoline N-oxide (4-NQO), as well as Cap1p constitutive transcriptional activation and nuclear retention (3, 68). Interestingly, the Cap1p-TR protein was shown to constitutively activate MDR1 expression in azole-susceptible C. albicans cells, demonstrating that MDR1 is a direct or indirect target of Cap1p (3). However, deleting CAP1 in an azole-resistant strain overexpressing MDR1 did not decrease MDR1 RNA levels (3), indicating that another transcription factor, possibly Mrr1p, is responsible for the constitutive overexpression of MDR1 in that strain. In addition, Cap1p is involved in protecting C. albicans against the oxidative stress induced by neutrophils during the course of the immune response (13, 24). The OSR in C. albicans involves oxidant sensing and response to oxidative damage via two major pathways that appear to act distinctly, namely, the Cap1p pathway and the high osmolarity glycerol (HOG) mitogen-activated protein kinase pathway (through a mechanism involving Ssk1p) (13, 22). These pathways respond differently to the OSR in a concentration- and/or oxidant-dependent manner, reflecting a complex process. For instance, while CAP1 is required for growth on both low and high concentrations of H2O2, HOG1 is required for growth only on high concentrations of peroxide (22). Also, a CAP1-deficient strain appears to be more susceptible to cadmium but more resistant to menadione than a HOG1-deficient strain (4). Genome-wide expression and proteomic studies showed that Cap1p regulates the expression of many genes involved in the OSR as well as other metabolic pathways, including energy metabolism and substance transport (3, 36, 63, 64); however, it was not determined whether Cap1p regulates these genes directly. In this paper, we used genome-wide location and expression analyses to better characterize the Cap1p regulon as well as Cap1p function in C. albicans.
|
|
|---|
derivatives were grown in YPD broth (1% yeast extract, 2% peptone, and 1% dextrose) supplemented or not with the indicated drug. The Escherichia coli MC1061 strain was used for DNA cloning and maintenance of the plasmid constructs. |
View this table: [in a new window] |
TABLE 1. Strains used in this study
|
C. albicans transformations. C. albicans transformations were conducted as described in MacPherson et al. (41), using a modified standard lithium acetate procedure. The transformed cells were plated on SC-ura plates and incubated for 3 days at 30°C.
Antifungal drugs and susceptibility testing. Stock solutions of FLC (a gift from Pfizer) and 4-NQO (Sigma) were prepared at concentrations of 5 mg/ml and 500 µM in water or dimethyl sulfoxide (DMSO), respectively. Drug susceptibility testing was performed using spot assays. Cells were grown overnight on SC-ura-met-cys plates and resuspended in water to an optical density at 600 nm (OD600) of 0.1. Serial dilutions (10-fold) of each strain were spotted onto SC-ura-met-cys plates supplemented with 2 µg/ml of FLC and 1.5 µM of 4-NQO or the solvent alone (water or DMSO, respectively). The plates were incubated for 2 days at 30°C.
Total protein preparation and Western blotting. Total protein extracts were prepared as described for S. cerevisiae (54) from 2 OD units of two independent strains expressing CAP1-HA3 (SGY243-CaEXP-CAP1-HA clones A and B) or CAP1-CSE-HA3 (SGY243-CAP1-CSE-HA clones A and B) (Table 1) grown overnight in SC-ura-met-cys (PMET3-inducing conditions) or SC-ura+met+cys (PMET3-repressing conditions). Extracts were boiled for 1 min, and 25-µl extracts were separated from total extracts of 100 µl by electrophoresis on a sodium dodecyl sulfate-10% polyacrylamide gel. Proteins were transferred to a nitrocellulose membrane with a Trans-Blot SD semidry transfer apparatus (Bio-Rad, Hercules, CA), and the membrane was incubated with a mouse anti-HA monoclonal antibody (12CA5; Roche) at a dilution of 1:1,000, followed by incubation with rabbit anti-mouse immunoglobulin G antibodies coupled to alkaline phosphatase (Bio-Rad). The membrane was then developed with 5-bromo-4-chloro-3-indolylphosphate p-toluidine salt and nitroblue tetrazolium chloride substrates, as recommended by the manufacturer (Bio-Rad).
ChIP-on-chip (ChIP-chip) and data analysis. Three independent cultures (50-ml each) of strains SGY243-CaEXP-A (untagged; control strain) and SGY243-CAP1-HA-A or SGY243-CAP1-CSE-HA-A (tagged strains) (Table 1) were grown overnight in SC-ura+met+cys diluted to an OD600 of 0.005 in SC-ura-met-cys (to induce PMET3) and grown until the OD600 reached 1.0. The subsequent steps of DNA cross-linking, DNA shearing, chromatin immunoprecipitation (ChIP), DNA labeling with Cy dyes, hybridization to intergenic DNA microarrays, and data analysis were conducted exactly as described in Liu et al. (38). Both pools of labeled DNA from the tagged strain (SGY243-CAP1-HA-A or SGY243-CAP1-CSE-HA-A; Cy5-labeled) and the corresponding untagged control strain (SGY243-CaEXP-A; Cy3-labeled) were mixed and hybridized to a C. albicans whole-genome tiled oligonucleotide DNA microarray described elsewhere (61). Hybridization was performed as recommended by the manufacturer (NimbleGen Systems, Inc). Scanning of the slides (n = 3) was performed using a GenePix 4000B scanner (Molecular Devices). Scanned images were preprocessed using NimbleScan software (version 2.4; NimbleGen Systems, Inc). General feature format reports were created for the Cy5 (tagged strain) and Cy3 (untagged control strain) intensity signals from each independent replicate and were then imported into the Tilescope program (http://tilescope.gersteinlab.org:8080/mosaic/pipeline.html) (69). Quantile normalization was applied to the data (69). The parameters used were as follows: a window size of 400 bp, a maximum genomic distance of 60 bp, and a minimum length of 120 bp. The replicate data were combined, and peak finding (i.e., determining the Cap1p binding sites) was done using a pseudomedian signal threshold of at least twofold and a P value cutoff of 0.01 or less (69).
Q-PCR for confirmation of the ChIP-chip data.
Quantitative real-time PCR (Q-PCR) was performed with three independent SGY243-CaEXP-A and SGY243-CAP1-HA-A or SGY243-CAP1-CSE-HA-A ChIP samples prepared as described above. Quantification of the recovered DNA was performed using a Quant-iT PicoGreen double-stranded DNA assay kit (Molecular Probes-Invitrogen) as previously described (38). The DNA concentration ranged from 0.08 ng/µl to 0.60 ng/µl for the tagged strains and 0.61 ng/µl to 1.22 ng/µl for the untagged strains. Q-PCR assays were conducted using Universal ProbeLibrary (Roche Applied Science) or TaqMan (Integrated DNA Technologies) methodology (38). The different primers and probe combinations used for Q-PCR are listed in Table 2. Optimal specific primer sequences and probes for the CIP1, IFR1 (targets), FUR1 (control for statistical analyses by t test), and SPS4 (orf19.7568; reference for normalization) promoters were obtained using Universal ProbeLibrary Web-based ProbeFinder software (version 2.34; Roche Applied Sciences) as previously described (38). Design of the TaqMan probe and specific forward and reverse primers for the MDR1 target promoter have been described previously (71). Q-PCR mixtures, Q-PCR conditions, and data analyses were performed as described previously (71). Statistical significance was determined using Welsh's two-sample t test. The statistical significance threshold was set at
= 0.05.
|
View this table: [in a new window] |
TABLE 2. Primers and probe sequences used for Q-PCR binding assays
|
0.02. For motif discovery analyses, DNA sequences covered by the 189 or 117 peaks identified in Cap1p-HA3 or Cap1p-CSE-HA 3binding data, respectively, were extracted and used as input for motif discovery, using the SCOPE (Suite for Computational Identification of Promoter Elements) program (http://genie.dartmouth.edu/scope/) (10, 12). This program allows accurate determination of potential transcription factor binding sites in a set of sequences using three different motif discovery algorithms (10, 12). To search for the TTASTAA motif within Cap1p-bound sequences, the same sequences analyzed with the SCOPE program were used as input for the DNA pattern matching TTASTAA, using the pattern-matching tool from Regulatory Sequence Analysis Tools ([RSAT] http://rsat.ulb.ac.be/rsat/). As a control, up to 1.0 kb of promoter sequences upstream of the ATG translation start site of the 6,093 promoters of the C. albicans ORFs was retrieved from the RSAT database (http://rsat.ulb.ac.be/rsat/). To search for the MTKASTMA sequence (where M designates A or C, K designates G or T, and S designates C or G) within the promoter region of the genes modulated in CJD21/PMK-CAP1-CSE versus CJD/PMK-CAP1, up to 1.0 kb of the promoter sequence upstream of the ATG translation start site of each gene was retrieved from the assembly 21 genome sequence such that overlap with neighboring ORFs was prevented and was then used as input for the DNA pattern matching MTKASTMA, using the RSAT pattern-matching tool. RNA isolation. Strains were grown overnight in 10 ml YPD at 30°C. The next day, an aliquot of the overnight culture was used to inoculate 200 ml YPD broth to a starting OD600 of 0.2. This culture was grown for 3 h as before, cells were collected by centrifugation, and cell pellets were immediately frozen and stored at –80°C until RNA isolation. Benomyl-exposed cultures were treated with 25 µg/ml benomyl (Sigma-Aldrich, St. Louis, MO) for 30 min before cells were harvested. Three independently obtained sets of cell cultures were used. RNA was isolated from frozen cell pellets using the hot-phenol method (57). Briefly, cells were resuspended in 12 ml AE buffer (50 mM sodium acetate [pH 5.2] and 10 mM EDTA) at room temperature, after which 800 µl 25% sodium dodecyl sulfate and 12 ml acid phenol (Fisher Scientific, Houston, TX) were added. The cell lysate was then incubated for 10 min at 65°C with vigorous shaking each minute, cooled on ice for 5 min, and subjected to centrifugation for 15 min at 11,952 x g. The supernatants were transferred to new tubes containing 15 ml chloroform, mixed, and subjected to centrifugation at 200 x g for 10 min. RNA was precipitated from the resulting aqueous layer by mixing that portion in new tubes with 1 volume of isopropanol and 0.1 volume of 2 M sodium acetate [pH 5.0] and subjecting the mixture to centrifugation at 17,211 x g for 35 min at 4°C. The supernatants were removed, the pellet was resuspended in 10 ml 70% ethanol, and the RNA was collected by centrifugation at 17,211 x g for 20 min at 4°C. The supernatants were again removed, and the RNA was resuspended in 50 to 200 µl diethyl pyrocarbonate-treated water. The RNA was stored at –80°C until needed.
cRNA synthesis and microarray hybridization. Immediately prior to cDNA and subsequent cRNA syntheses, the purity and concentration of RNA samples were determined from A260/A280 readings and RNA integrity was determined by capillary electrophoresis, using an RNA 6000 Nano laboratory-on-a-chip kit and 2100 bioanalyzer (Agilent Technologies) per the manufacturer's instructions. First- and second-strand cDNA was synthesized from 15 µg total RNA, using a SuperScript double-stranded cDNA synthesis kit (Invitrogen) and an oligo-dT24-T7 primer (Proligo) according to the manufacturers' instructions. cRNA was synthesized and labeled with biotinylated UTP and CTP by in vitro transcription, using T7 promoter-coupled double-stranded cDNA as a template and a BioArray HighYield RNA transcript labeling kit (Enzo Diagnostics). Double-stranded cDNA synthesized from the previous steps was washed twice with 70% ethanol and suspended in 22 µl of RNase-free water. The cDNA was incubated as recommended with reaction buffer, biotin-labeled ribonucleotides, dithiothreitol, RNase inhibitor mix, and T7 RNA polymerase for 5 h at 37°C. The labeled cRNA was separated from unincorporated ribonucleotides with a Chroma Spin-100 column (Clontech) and was ethanol precipitated at –20°C overnight.
The cRNA pellet was suspended in 10 µl of RNase-free water, and 10 µg was fragmented at 95°C for 35 min in 200 mM Tris-acetate (pH 8.1), 500 mM potassium acetate, and 150 mM magnesium acetate. The fragmented cRNA was hybridized for 16 h at 45°C to either C. albicans Affymetrix GeneChip arrays (CAN04a530004N; manufactured by NimbleExpress) or to an Affymetrix custom expression array (CAN07a520619F; manufactured by Affymetrix) for C. albicans. Arrays were washed at 25°C with 6x SSPE (1x SSPE is 0.18 M NaCl, 10 mM NaH2PO4, and 1 mM EDTA [pH 7.7]) and 0.01% Tween 20 followed by a stringent wash at 50°C with 100 mM MES (morpholineethanesulfonic acid), 0.1 M NaCl, and 0.01% Tween 20. Affymetrix Fluidics Station 450 was used for hybridizations and washes according to standard EukGE-WS2v5 protocol. The arrays were then stained with phycoerythrin-conjugated streptavidin (Molecular Probes), and the fluorescence intensities were determined using a GCS 3000 high-resolution confocal laser scanner (Affymetrix). The scanned images were analyzed using software resident in the GeneChip operating system, version 2.0 (Affymetrix). Sample loading and variations in staining were standardized by scaling the average of the fluorescent intensities of all genes on an array to a constant target intensity of 250. The signal intensity for each gene was calculated as the average intensity difference, represented by
(PM – MM)/number of probe pairs, where PM and MM denote perfectly matched and mismatched probes, respectively.
Gene expression microarray data analysis. The scaled gene expression values from GeneChip operating system version 2.0 software were imported into GeneSpring 7.2 software (Agilent Technologies) for preprocessing and data analysis. Probe sets were deleted from subsequent analysis if they were called absent by the Affymetrix criterion and displayed an absolute value below 20 in all experiments. The expression value of each gene was normalized to the median expression of all genes in each chip as well as to the median expression for that gene across all chips in the study. A pairwise comparison of gene expression was performed for each matched experiment.
Q-PCR for expression data.
Real-time PCR was performed in follow-up experiments to validate the microarray results. First-strand cDNAs were synthesized from 2 µg of total RNA in a 21-µl volume of reaction mixture, using a SuperScript first-strand synthesis system for reverse transcription (RT)-PCR (Invitrogen, Carlsbad, CA) in accordance with the manufacturer's instructions. Q-PCRs were performed in triplicate using a 7000 sequence detection system (Applied Biosystems, Foster City, CA). Independent PCRs were performed using the same cDNA for both the gene of interest and the 18S rRNA, using SYBR green PCR master mix (Applied Biosystems). Gene-specific primers were designed for the gene of interest and the 18S rRNA, using Primer Express software (Applied Biosystems) and an oligo analysis and plotting tool (Qiagen, Valencia, CA), and are shown in Table 3. The PCR conditions consisted of AmpliTaq Gold activation at 95°C for 10 min, followed by 40 cycles of denaturation at 95°C for 15 s and annealing/extension at 60°C for 1 min. A dissociation curve was generated at the end of each PCR cycle to verify that a single product was amplified, using software provided with the 7000 sequence detection system. The change in fluorescence of SYBR green I dye in every cycle was monitored by the system software, and the threshold cycle (CT) above the background for each reaction was calculated. The CT value of 18S rRNA was subtracted from that of the gene of interest to obtain a
CT value. The
CT value of an arbitrary calibrator (e.g., untreated sample) was subtracted from the
CT value of each sample to obtain a 
CT value. The gene expression level relative to that of the calibrator was expressed as 2–
CT. Statistical analysis was performed using R software, version 2.5.0 (www.r-project.org). Relative changes were compared using a two-sample t test. The statistical significance threshold was set at
= 0.05.
|
View this table: [in a new window] |
TABLE 3. Primers used for Q-PCR expression analysis
|
|
|
|---|
![]() View larger version (33K): [in a new window] |
FIG. 1. Strategy for tagging Cap1p with a HA3 epitope and characterization of the tagged strains. (A) Schematic representation of the CAP1-tagging cassette. The wild-type or mutated version of the CAP1 gene [CAP1 (WT/CSE); light gray box], cloned as an in-frame fusion with a HA3 tag (dark gray box) in plasmid pCaEXP (9), is under the control of the MET3 promoter (PMET3; open arrow) and is followed by the C. albicans URA3 marker (black oval). The 5' and 3' fragments of the RP10 gene (5'RP10 and 3'RP10; black boxes) flank the cassette and allow targeted integration at the RP10 locus (9). (B) Western blot analysis of strains expressing HA3-tagged versions of the CAP1 gene (CAP1-HA3 or CAP1-CSE-HA3). Total proteins were extracted from two independent clones of the SGY243 transformants (A and B) grown in the absence (–) or presence (+) of 2.5 mM methionine (Met). Western blotting was performed using the anti-HA antibody 12CA5. Positions of the molecular mass standards are indicated on the left (kDa). (C) Drug resistance profiles of C. albicans strains expressing HA3-tagged CAP1 alleles. Two independent transformants (A and B) for each of the CAP1-HA3- or CAP1-CSE-HA3-expressing strains or the strain carrying the empty vector as the negative control (vector) were analyzed by spot assay for their ability to grow on SC-ura-met plates in the absence or presence of 2 µg/ml of FLC or 1.5 µM of 4-NQO. The plates were incubated for 2 days at 30°C.
|
0.01) (see Tables S1 and S2 in the supplemental material). We visualized the ChIP-chip data using a C. albicans genome browser representing the entire assembly 19 (see Materials and Methods). We found that a high proportion of Cap1p binding peaks clearly associated with ORFs, totaling 89 target genes for Cap1p-HA3 or Cap1p-CSE-HA3 (60 of the 89 genes were common to both proteins; 23 additional genes were specific to Cap1p-HA3, while 6 additional genes were specific to Cap1p-CSE-HA3) (Fig. 2; see also Tables S1 and S2 in the supplemental material). In some cases, more than one peak associated with one ORF, while in one occurrence one peak associated with two ORFs (orf19.3121 and orf19.3122) (see Tables S1 and S2 in the supplemental material). Interestingly, we also found that 101 Cap1p-HA 3binding peaks and 35-Cap1p-CSE-HA 3binding peaks did not clearly associate with defined ORFs, including peaks that were located in intergenic regions (see the bottom panel of Fig. 2 for an example), suggesting that these regions may encode unidentified ORFs or small RNAs. Strikingly, Cap1p binding was detected not only at the promoter region of its target genes but also at their 3' ends and within their ORFs (Fig. 2). This surprising binding profile could be explained by the association of Cap1p with chromatin-associated proteins and/or the transcriptional machinery (see Discussion). Finally, it is noteworthy that the Cap1p-CSE-HA3 signal was often increased relative to that of Cap1p-HA3 at common targets (Fig. 2; see also Tables S1 and S2 in the supplemental material).
![]() View larger version (18K): [in a new window] |
FIG. 2. Cap1p binding at selected C. albicans genomic regions. Plotted are the normalized log2-transformed signal intensities (lower graphs of each panel) of HA3-tagged wild-type (WT) or hyperactive (CSE) Cap1p binding (y-axis) versus the corresponding position of each signal (x-axis) in selected C. albicans genomic regions from assembly 19 (the corresponding contig 19 number is indicated at the top of each panel). Log2-transformed signal intensity values are indicated at the left of the y-axis. The y-axis intercept is the value 0 (i.e., a binding ratio of 1). The location of each selected region from the corresponding contig 19 is shown on the scaled upper axis of each panel in kilobases (K). The spacing is 1.0 kb between each major graduation and 0.2 kb between each minor graduation. The orientation of each ORF is depicted by the arrowed gray rectangle. Negative enrichment values in the FLU1, RTA2, and PDR16 panels may be due to background noise, preferential amplification, or normalization biases inherent to the ChIP-chip technology (6, 51).
|
|
View this table: [in a new window] |
TABLE 4. Overrepresented functional categories in Cap1p ChIP-chip data
|
![]() View larger version (13K): [in a new window] |
FIG. 3. Quantification of the in vivo enrichment of Cap1p-HA3 (Cap1-HA) and Cap1p-CSE-HA3 (Cap1-CSE-HA) binding at the MDR1, CIP1, IFR1, and FUR1 targets using Q-PCR. SGY243-CaEXP-A (control strain; untagged) and SGY243-CAP1-HA-A or SGY243-CAP1-CSE-HA-A (tagged strains) were submitted to ChIP (three biological replicates), and the recovered DNA samples were analyzed by Q-PCR, using Universal ProbeLibrary probes (Roche Diagnostics) for the CIP1, IFR1, FUR1, and SPS4 (reference for normalization) targets and a TaqMan probe for the MDR1 target. Relative enrichment values (n-fold) are presented in logarithmic scale. Error bars denote standard deviations. Asterisks denote statistical significance determined by using Welsh's two-sample t test (P 0.05).
|
![]() View larger version (22K): [in a new window] |
FIG. 4. Motif logo of a conserved sequence in Cap1p-HA3- or Cap1p-CSE-HA3-enriched DNA fragments. The 189 Cap1p-HA3- or 117 Cap1p-CSE-HA3-bound sequences were used as input for motif discovery, using the SCOPE program (http://genie.dartmouth.edu/scope/) (10, 12). The highest scoring motif MTKASTMA common to both conditions (Cap1p-HA3 and Cap1p-CSE-HA3) is shown.
|
2.0, (ii) their expression changed by at least 2.0-fold in each experiment, and (iii) the change (n-fold) was considered statistically significant by Student's t test. Based on these criteria, we found 51 differentially expressed genes (Table 5), including GLR1, GRE2, OYE32, CIP1, MDR1, orf19.2262, orf19.3537, orf19.5060, orf19.1167, orf19.1340, orf19.5517, and orf19.6757, all of which are involved in responses to oxidative stress or coding for oxidoreductases. Indeed, using the GO term Finder tool at the CGD, we found that "Response to oxidative stress" and "Oxidoreductase activity" were the most significantly overrepresented GO terms among the 51 modulated genes (P = 0.004 and P = 0.01, respectively). We performed real-time RT-PCR assays to validate the expression array results (see Materials and Methods) and confirmed the differential expression of MDR1 and GLR1 with relative expression changes (± standard errors) of 13.1 ± 0.3 and 1.9 ± 0.1, respectively (Fig. 5). |
View this table: [in a new window] |
TABLE 5. List of genes whose expression is modulated in CJD21/PMK-CAP1-CSE versus CJD21/PMK-CAP1
|
![]() View larger version (16K): [in a new window] |
FIG. 5. Quantitative real-time RT-PCR analysis of (A) MDR1 and (B) GLR1 genes differentially expressed in the microarray experiments. Bars in each graph indicate log-transformed relative changes in RNA expression of the samples indicated compared to their controls: CJD21/PMK-CAP1-CSE versus CJD21/PMK-CAP1, CJD21/PMK plus benomyl (Ben) versus CJD21/PMK plus DMSO and CJD21/PMK-CAP1 plus Ben versus CJD21/PMK-CAP1 plus DMSO. Asterisks denote statistical significance by the t test (P 0.05). Error bars denote standard errors.
|
Transcriptional response of Cap1p targets to benomyl treatment. Previous studies have shown that the genes responding to benomyl overlap significantly with those responding to H2O2 and thus to oxidative stress (30, 36), suggesting that benomyl is an inducer of Cap1p activity. To identify genes whose expression is induced by benomyl in a CAP1-dependent manner, we examined the gene expression profiles of the CAP1-expressing C. albicans strain CJD21/PMK-CAP1 treated with benomyl to those of CJD21/PMK-CAP1 treated with DMSO and those of the CAP1-deficient related strain CJD21/PMK also treated with benomyl to those of CJD21/PMK treated with DMSO (see Materials and Methods). When the gene expression profiles of benomyl-treated and diluent-treated CJD21/PMK-CAP1 were compared using our criteria (see Materials and Methods and Table S3 in the supplemental material), 432 genes were found to be differentially expressed. Similarly, when benomyl-treated and diluent-treated CAP1-deficient C. albicans gene expression profiles were compared using our criteria (see Materials and Methods and Table S4 in the supplemental material), 232 genes were found to be differentially expressed. By considering only those genes differentially expressed in the CAP1-expressing strain in response to benomyl compared to those differentially expressed in CJD21/PMK in response to benomyl, we found a set of 115 genes whose expression fit these criteria, qualifying them as genes induced by benomyl in a CAP1-dependent manner (see Table S5 in the supplemental material). We validated the expression microarray data by confirming the CAP1-dependent benomyl induction of MDR1 and GLR1, using real-time PCR (Fig. 5). This experiment showed that strain CJD21/PMK-CAP1 displayed a 517-fold induction of MDR1 expression in response to benomyl treatment, while the CAP1-deficient mutant CJD21/PMK displayed a 133-fold increase in MDR1 RNA upon treatment with benomyl (Fig. 5A), most likely due to the contribution of another transcriptional regulator (see Discussion). In addition, benomyl-treated CJD21/PMK-CAP1 displayed a 4.2-fold change in GLR1 expression, whereas CJD21/PMK treated with benomyl showed no significant change in GLR1 expression (Fig. 5B), demonstrating that Cap1p is the major regulator of GLR1 expression in response to benomyl.
Interestingly, among the 115 genes responding to benomyl in a CAP1-dependent manner, 31 were also identified as target genes of Cap1p-HA3 or Cap1p-CSE-HA3 in our ChIP-chip data, as illustrated in the Venn diagram representing the overlap between the ChIP-chip data set and the benomyl-induced CAP1-dependent gene data set (Fig. 6). These 31 genes include CAP1 itself, MDR1, EBP1, OYE32, OYE23, and GLR1 (Table 6), all of which are involved in the OSR. Taken together, these results indicate that benomyl induces Cap1p activity and suggest that it activates the OSR via Cap1p in C. albicans (see Discussion).
![]() View larger version (15K): [in a new window] |
FIG. 6. Venn diagram of the overlap between the benomyl-responsive CAP1-dependent gene data set (left circle) and the Cap1p-HA3 (right upper circle) and Cap1p-CSE-HA3 (right lower circle) ChIP-chip datasets. The numbers in the Venn diagram indicate the numbers of genes.
|
|
View this table: [in a new window] |
TABLE 6. List of Cap1p-bound targets whose expression is modulated upon benomyl treatment in a CAP1-dependent manner
|
|
|
|---|
A previous study was performed to compare genes coregulated with MDR1 in azole-resistant C. albicans strains with those induced in response to benomyl (30). Several genes that responded to benomyl treatment in a CAP1-dependent manner (see Table S5 in the supplemental material) were shown to be benomyl responsive by Karababa et al., including the stress genes orf19.3121, OYE23, EBP1, orf19.2262, OYE32, IFR2, orf19.251, orf19.5517, and TTR1 and the transporter genes MDR1, ARR3, and SNQ2 (30). In addition, genes of unknown function that were defined as CAP1 dependent upon benomyl treatment in our study were also found to be benomyl responsive by Karababa et al., for example, orf19.2285, orf19.6898, PRN1, orf19.7042, orf19.6586, YIM1, orf19.2043, and orf19.1162. Using our criteria (P
0.01; binding ratio,
2), we found a majority of these genes to be bound in vivo by Cap1p (except TTR1, SNQ2, PRN1, orf19.2043, and orf19.1162) (see Tables S1 and S2 in the supplemental material), indicating a direct transcriptional regulation by Cap1p. Although using a microarray with only partial coverage of the C. albicans genome, another previous study examined the C. albicans response to hydrogen peroxide, an inducer of oxidative stress (64). As expected, several genes responsive to benomyl treatment in a CAP1-dependent manner overlapped with those responsive to H2O2 treatment, including the stress genes EBP1, OYE23, orf19.2262, OYE32, orf19.5517, TRR1, orf19.251, IFR2, and GLR1. Another report by Wang et al. aimed at identifying genes differentially expressed under untreated conditions, also termed basal or stress-absent conditions, in a wild-type C. albicans strain versus those in a cap1
/cap1
mutant (63). This study identified 48 downregulated genes in the cap1 mutant relative to those in the wild-type strain, among which only EBP1, OYE23, OYE32, orf19.2262, ARR3, ZRT2, and orf19.868 were bound by Cap1p in the present study (see Tables S1 and S2 in the supplemental material), suggesting that many genes found by Wang et al. that were differentially expressed under untreated conditions are indirect Cap1p targets.
A striking finding was that Cap1p binds not only to the promoter region of its target genes but also within the ORF and the 3' region, including the transcriptional termination region (Fig. 2; see also Tables S1 and S2 in the supplemental material). Interestingly, binding of Cap1p-CSE-HA3 within the ORFs was more frequent than that of Cap1p-HA3 (Fig. 2; see also Tables S1 and S2 in the supplemental material), suggesting that activation of Cap1p enhances its propensity to occupy intragenic regions. To our knowledge, this is the first report of such a transcription factor binding profile in C. albicans. Previous studies reported similar associations of transcription factors to both promoter and coding regions of their targets, including the tumor suppressor p53 and the estrogen receptor (5, 11). Interestingly, binding of both transcription factors to intragenic regions overlapped with that of RNA polymerase II binding sites (5, 11). Moreover, it was shown that p53 physically associated with RNA polymerase II, with which it travels within target gene loci (5). The intragenic binding profile of Cap1p as well as its binding at the 3' end of target genes is consistent with a model in which Cap1p interacts with the transcriptional machinery and travels with the RNA polymerase II across the transcribed region of the target locus in a manner similar to that of p53. Interestingly, it was shown that Yap1p physically interacts with the general transcription factor TFIIA in S. cerevisiae and that TFIIA mutants exhibited an enhanced susceptibility to oxidative stress (32), highlighting a direct link between Yap1p and the RNA polymerase II complex. Another possibility could be that Cap1p interacts with chromatin-associated proteins.
Previous studies have implicated the zinc cluster transcription factors Mrr1p and Upc2p in regulating MDR1 gene expression in C. albicans (19, 20, 47, 71). While Mrr1p appears to act as a potent transcriptional activator of MDR1 in C. albicans azole-resistant clinical isolates overexpressing MDR1 (47), Upc2p appeared to act as a moderate activator or a repressor of MDR1, depending upon the activating signal (71). Consistently, a gain-of-function mutation in UPC2 from a C. albicans azole-resistant strain was shown to cause a moderate upregulation of MDR1 (20). We previously showed, using ChIP experiments, that Upc2p binds to the MDR1 promoter (71), whereas it is still unknown whether Mrr1p associates directly with the promoter region of MDR1. Based on previous studies and the present study, Cap1p appears to be a potent activator of MDR1 expression, as reflected by Northern blot and luciferase reporter analyses (3, 53) and by our expression microarray data (Table 5; see also Table S5 in the supplemental material). Another regulator that was shown to participate in the regulation of MDR1 expression is Mcm1p, which binds directly to an Mcm1p binding motif found in the promoter of MDR1 (52), a finding that was recently confirmed in vivo by genome-wide location analyses (37). Taken together, these observations suggest that a complex network of transcriptional regulators, including Cap1p, Mrr1p, Upc2p, and Mcm1p, is involved in the regulation of MDR1 in C. albicans. Studies are ongoing to determine whether Cap1p interacts physically with these other factors to regulate MDR1 expression.
Among the overrepresented functional categories of genes bound in vivo by Cap1p were "Response to drug," "Hyphal cell wall," "Phospholipid transport," and "Regulation of nitrogen utilization" (Table 3), suggesting other roles for Cap1p besides the OSR. Of particular interest was the grouping of PDR16, RTA2, MDR1, and FLU1 into the functional category "Response to drug." These genes were previously reported to participate in azole resistance in C. albicans (8, 29, 55, 66). While PDR16 and MDR1 are involved in clinical azole resistance (55, 66), FLU1 and RTA2 were shown to modulate cell susceptibility to azoles in C. albicans (8, 29), suggesting that Cap1p could confer azole resistance via PDR16, FLU1, RTA2, and/or MDR1 in C. albicans. Although no evidence of Cap1p gain-of-function mutations in azole-resistant clinical isolates of C. albicans has been reported to date, a recent study has reported the correlation of constitutive overexpression of FLU1 together with MDR1 with the development of clinical azole resistance (25). Whether upregulation of these two major facilitator-encoding genes is mediated by Cap1p remains to be determined. Another overrepresented functional category was "Hyphal cell wall," which included SSA2, EBP1, orf19.251, PDC11, and ADH1. As Cap1p appears to be important for protecting C. albicans against the oxidative stress induced by neutrophils during the course of the immune response (24), it might regulate processes contributing to virulence, including hyphal growth. Neutrophils induce oxidative stress through the release of nitric oxide anions; thus, it would not be surprising if Cap1p regulates processes involved in nitrogen metabolism, as reflected by the enrichment of the functional category "Regulation of nitrogen utilization," which included orf19.2693, orf19.3121, and GST3. Finally, "Phospholipid transport" was also among the overrepresented functional categories and might reflect a role for Cap1p in enhancing the dynamics of phospholipid metabolism, as these molecules are highly sensitive to oxidative stress.
Benomyl is an aneuploidogen, a toxic antimitotic drug also used as an antifungal agent that is thought to exert its effect by binding to tubulin and inhibiting tubulin assembly (44). Based on previous studies as well as the present study, it appears that benomyl induces the OSR in C. albicans. First, genes differentially expressed upon exposure to H2O2 overlap with those responsive to benomyl (30, 36). Second, we showed here that benomyl induces Cap1p activity, reflected by the upregulation of genes involved in oxidative stress in a CAP1-dependent manner. Indeed, the GO term "Response to oxidative stress" was among the most significantly enriched functional categories (P = 0.026) of the CAP1-dependent genes responding to benomyl (data not shown). A study with S. cerevisiae also showed that a subset of OSR genes was upregulated rapidly after the addition of benomyl (39). Interestingly, this transcriptional response involved Yap1p as the central regulator of the early response to benomyl treatment (39). Lucau-Danila et al. suggested that benomyl potentially acts at the level of Yap1p nuclear localization rather than at the level of DNA binding (39). Thus, it is possible that benomyl induces conformational changes within the Cap1p CRD by directly binding to the C-terminal domain of Cap1p, leading to nuclear retention of Cap1p and transcriptional activation. However, an effect on intracellular redox balance is not excluded, as the benomyl metabolite n-butylisocyanate, a cleavage product of benomyl, results in inhibitory effects on dehydrogenases or glutathione reductases (44), which in turn may activate the OSR via Cap1p.
This work was supported by research grants from the Canadian Institutes of Health Research (CIHR) to M.R. (grants MT-15679 and HOP-67260) and from the National Institutes of Health (R01 AI058145) to P.D.R. S.Z. is supported by a doctoral studentship from the Fonds de la Recherche en Santé du Québec (FRSQ). IRIC is supported by the Canadian Center of Excellence in Commercialization and Research, the Canadian Foundation for Innovation, and the Fonds de Recherche en Santé du Québec.
Published ahead of print on 24 April 2009. ![]()
Supplemental material for this article may be found at http://ec.asm.org/. ![]()
Present address: Génome Québec, Montréal, Quebec, Canada H3B 1S6. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2010 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»