Previous Article | Next Article ![]()
Journal of Bacteriology, April 2008, p. 2939-2946, Vol. 190, No. 8
0021-9193/08/$08.00+0 doi:10.1128/JB.01847-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
,
Laboratory of Bacterial Genomics, Department of Pathobiological Sciences, University of Wisconsin—Madison, 1656 Linden Drive, Madison, Wisconsin 53706
Received 23 November 2007/ Accepted 3 February 2008
|
|
|---|
|
|
|---|
To characterize the response of M. tuberculosis to the host environment, we previously utilized microarrays to compare the transcriptome of M. tuberculosis in a murine model to the transcriptome of M. tuberculosis growing in culture (31). A 34-kb region of the genome was identified as a region that is predominantly upregulated in mouse lungs compared to in vitro culture. This region was designated the in vivo-expressed genomic island (iVEGI) (31). Several genes in the iVEGI have been shown to be relevant to the pathogenesis of M. tuberculosis. For instance, the iVEGI genes mprA and mprB encode a two-component signal transduction system required for infection (40). Within the iVEGI, the Rv0967 gene was predicted to encode a transcriptional regulator. Using biochemical and structural experiments, it was shown that the Rv0967 gene encodes a copper-binding transcriptional repressor protein (14). Because the repressor regulates its own four-gene operon in a copper-dependent manner, this operon was designated cso (copper sensitive operon), and the Rv0967 transcriptional regulator was designated CsoR (cso repressor). The induction of cso genes during murine tuberculosis suggested that the copper concentration might fluctuate within the phagosome and that the ability to respond to this fluctuation could be important for mycobacterial virulence. In fact, a previous study of cultured macrophages showed that the level of intraphagosomal Cu increases from 25 to 500 µM after phagocytosis of M. tuberculosis, but no such increase occurs after phagocytosis of less virulent species of Mycobacterium (36).
To our knowledge, Cu homeostasis in M. tuberculosis has not been studied previously. However, bacterial mechanisms of copper homeostasis have been established in two model systems: the K-12 laboratory strain of Escherichia coli and Enterococcus hirae, a gram-positive, extracellular member of the intestinal flora (13, 26, 39). The mechanisms of copper homeostasis identified in these systems include copper transporters (e.g., CopA [20]), the detoxifying enzyme multicopper oxidase (CueO [7]), copper-sensing regulatory proteins (e.g., CueR and CopY [18, 28]), and copper chaperones (e.g., CopZ [5]). Based on previous studies, the most important aspect of copper homeostasis is considered to be copper exporters because of their ability to pump excess Cu outside the cell before it can damage intracellular components (32). When free inside the cell, copper can facilitate the generation of toxic reactive oxygen species, and Cu(I) has been shown to damage proteins, particularly the proteins containing thiol groups (19). Furthermore, excess Cu nonspecifically binds to proteins that require other cofactors (e.g., iron) for activity (32). In this study, we further characterized the cso genes and obtained evidence that CtpV is a copper export transporter. Further, we elucidated the whole-genome transcriptional response of M. tuberculosis to variable levels of Cu by identifying both novel and known metalloregulatory proteins triggered by Cu stimulation. Based on our results, we began to sketch the components utilized by M. tuberculosis to survive the fluctuation in Cu levels inside the phagosome.
|
|
|---|
RNA extraction and cDNA synthesis and labeling. RNA was extracted by resuspending mycobacterial cell pellets in 1 ml Trizol (Invitrogen, Carlsbad, CA) per 5-ml culture. The suspensions were placed in 2-ml screw-cap tubes containing 1 g zirconia/silica beads (Biospec Products, Inc., Bartlesville, OK). The tubes were pulsed twice for 45 s with a mini-Bead-Beater-8 (Biospec Products, Inc.). Following processing with the Bead-Beater, RNA extraction was performed as recommended by the manufacturer (Invitrogen), with modifications described previously (38). Extracted RNA was treated with DNase I (Ambion, Austin, TX) until no DNA was detected using PCR primers for the 16S rRNA gene (see Table S1 in the supplemental material). For qRT-PCR, cDNA was synthesized from 1 µg of total RNA using SuperScript III (Invitrogen) as directed by the manufacturer in the presence of 250 ng of mycobacterial genome-directed primers (30). For microarrays, double-stranded cDNA (ds-cDNA) was synthesized from 10 µg of total RNA using an Invitrogen SuperScript ds-cDNA synthesis kit as directed by the manufacturer in the presence of 250 ng genome-directed primers (30). The ds-cDNA was cleaned and labeled using the NimbleGen gene expression analysis protocol (NimbleGen Systems, Inc., Madison, WI). Briefly, ds-cDNA was incubated with 10 ng RNase A (Novagen) at 37°C for 10 min and cleaned using phenol-chloroform extraction, followed by ethanol precipitation. For Cy3 labeling, 1 µg ds-cDNA was incubated for 10 min at 98°C with 1 OD600 unit of Cy3-9mer Wobble primer (TriLink Biotechnologies, San Diego, CA). Then 8 mmol of deoxynucleoside triphosphates and 100 U of the Klenow fragment (New England Biolabs, Ipswich, MA) were added, which was followed by incubation at 37°C for 2 h. The reaction was stopped by adding 0.1 volume of 0.5 M EDTA, and the labeled ds-cDNA was cleaned by isopropanol precipitation.
DNA microarrays. Microarray chips were purchased from NimbleGen Systems, Inc., and they contained 19 replicates of each 60-mer oligonucleotide probe designed for the 3,989 open reading frames in the genome of M. tuberculosis H37Rv (3). Further, the whole genome was represented five times on each chip (i.e., five technical replicates/chip) so that there were a total of 95 probes per gene. Hybridization of 3 µg of ds-cDNA with NimbleGen hybridization buffer and Nimblegen hybridization component A was performed in hybridization chambers (TeleChem International, Inc., Sunnyvale, CA) overnight at 42°C. Following hybridization, washing was performed using NimbleGen wash solutions I, II, and III as recommended by the manufacturer. Slides were scanned using an Axon GenePix 4000B scanner (Molecular Devices Corporation, Sunnyvale, CA), and fluorescence intensity levels were determined using NimbleScan (NimbleGen) and normalized to a mean value of 1,000. Significantly changed genes were identified using a flexible empirical Bayes model (specifically, the LNN model in the EBArrays package in R [http://www.bioconductor.org]). A cutoff value of 0.50 for the probability of differential expression was used to determine significantly changed genes (12). To determine functional groups of genes enriched in the microarray data set, identified genes were assigned TIGR roles (http://cmr.tigr.org). Statistical enrichment of each role category within the microarray genes compared with the whole H37Rv genome was performed using a standard hypergeometric distribution function in Microsoft Excel. Complete microarray data are shown in Table S4 in the supplemental material and are available online at http://www.ahabs.wisc.edu/Faculty/talaat-a/lab/data.php.
qRT-PCR.
A SYBR green-based qRT-PCR protocol was used for confirmation of the microarray results. SYBR green qRT-PCR was performed using 100 ng cDNA as the template in a reaction with iTaq SYBR green Supermix with ROX (Bio-Rad Laboratories, Hercules, CA) in the presence of gene-specific primers (see Table S1 in the supplemental material) at a concentration of 200 nM. For experiments to determine metal ion specificity and induction of cso, a TaqMan-based qRT-PCR protocol was used. TaqMan qRT-PCR was performed using 100 ng cDNA as the template in a reaction with Platinum quantitative PCR SuperMix-UDG with ROX (Invitrogen) in the presence of gene-specific primers at a concentration of 200 nM and 6-carboxyfluorescein-labeled TaqMan probes (Biosearch Technologies, Novato, CA) (see Table S1 in the supplemental material). For both SYBR green and TaqMan qRT-PCR, the cycle conditions were 50°C for 2 min, 95°C for 3 min, and 40 cycles of 95°C for 15 s and 60°C for 30 s. Reactions were performed in triplicate with an AB7300 with fluorescence read at the 60°C step. The threshold cycle values were normalized to levels of 16S rRNA transcripts and then expressed as changes compared to the metal-free sample (
CT method) (30).
|
|
|---|
![]() View larger version (18K): [in a new window] |
FIG. 1. Growth of M. tuberculosis in the presence of different levels of Cu. Cultures were inoculated to obtain an OD600 of 0.1 and then grown for 14 days ( 14 doubling times) in the presence of 0, 5, 50, or 500 µM CuCl2. Colony counts were determined by plating on Middlebrook 7H10 medium supplemented with ADC. The results for a representative sample of two biological replicates are shown. The error bars indicate standard deviations of the mean colony counts.
|
![]() View larger version (12K): [in a new window] |
FIG. 2. Metal induction of genes in the cso operon. Cultures were grown to an OD600 of 0.6 in metal-free (Sauton's) medium and exposed to the metal ions indicated for 3 h prior to RNA extraction and TaqMan qRT-PCR analysis. Expression values were normalized to 16S rRNA values and are expressed as changes relative to the culture with no metal added. (A) qRT-PCR showing transcriptional induction of csoR and ctpV as the Cu concentration increases. (B) Induction of the cso operon in the presence of Cu and other divalent cations (Ag, Ni, and Zn). Concentrations were chosen based on the upper limit of toxicity for mycobacteria (14).
|
More in-depth sequence analysis revealed that, like other metal translocation P-type ATPases, CtpV is predicted to have eight transmembrane helices with three major cytoplasmic domains (11). Within these domains are conserved motifs important for ATPase activity and substrate translocation (Fig. 3). Helix 6 contains the residues thought to confer metal specificity in the family of metal translocation P-type ATPases, and CtpV contains the CPCALGLA motif found in most other Cu transporters (26). CtpV does not, however, contain the cytoplasm-facing Cu binding sites common in metal translocation P-type ATPases. These metal binding domains have been proposed to aid in the transfer of Cu ions from Cu chaperones to the transporter (10), but deletion of them in two separate human Cu transporters, the Wilson protein and the Menkes protein, did not affect Cu translocation activity (4, 35). The significance of the apparent loss of the metal binding domain motif in CtpV is unknown. The directionality of ion transport (import versus export) cannot be determined based on sequence data alone. However, the induction of ctpV with increasing concentrations of Cu and particularly its high level of induction with 500 µM Cu suggested that it may function as a Cu export pump.
![]() View larger version (23K): [in a new window] |
FIG. 3. Predicted structure and alignment of CtpV with experimentally characterized Cu transporters in E. coli (Ec) and E. hirae (Eh). The A, P, and N domains and the underlined motifs are common to P-type ATPases (34). The motifs in bold type are specific to metal transport P-type ATPases (27). Helix 6 motif is associated with metal specificity (33). MBD, metal binding domain; Mtb, M. tuberculosis.
|
![]() View larger version (16K): [in a new window] |
FIG. 4. Microarray analysis of the Cu-responsive transcriptome in M. tuberculosis. (A) Reproducibility of technical replicates in NimbleGen-based arrays following hybridization with Cu-treated samples. The correlation (r > 0.95) between signals of two copies of the genome on a chip hybridized with a sample containing no Cu is shown. (B) Correlation (r > 0.67) of signals between biological replicates. The results for two biological replicates for samples containing no Cu are shown. (C) Venn diagram showing common and unique sets of Cu-responsive genes. The overlap with a previously identified oxidative stress data set (24) is shown for a comparison to genes induced with both low and high levels of Cu or with only high levels of Cu. As expected, most of the overlap occurs at the toxic concentration of Cu (500 µM).
|
|
View this table: [in a new window] |
TABLE 1. Genes identified as genes that are significantly up- or downregulated in medium containing 50 or 500 µM CuCl2 compared with Cu-free medium
|
Transporters involved in copper response. By assigning the 30 genes in the Cu transcriptome to functional groups (http://cmr.tigr.org), we were able to determine which biological functions are statistically overrepresented in our list of Cu-responsive genes (see Table S3 in the supplemental material). Only the transport and binding protein category was identified as significantly enriched (P < 0.05) in response to Cu, and four putative transporters were identified as proteins that were significantly upregulated in the presence of 500 µM Cu. Only one of these transporters, CtpV, is predicted to specifically transport Cu, although two permease proteins (Rv0849 and Rv2963) that could potentially serve as nonspecific ion transporters were identified. Unexpectedly, none of the other predicted metal-translocating P-type ATPases were identified as proteins that were induced by Cu, including two ATPases (CtpA and CtpB) with the predicted Cu-specific helix 6 motif. The lack of induction of CtpA and CtpB was confirmed by qRT-PCR (data not shown).
Copper-responsive enzymes. Many enzymes were identified in the Cu transcriptome, and most of them were induced upon addition of Cu. With 70-fold induction after exposure to Cu, cadI is by far the gene that is induced most in the Cu-responsive transcriptome. Bioinformatic analyses identified CadI as a putative glyoxalase I metalloenzyme required for the conversion of methylglyoxal, a toxic by-product of glycolysis, to lactate. Because methylglyoxalase enzymes generally contain zinc, it is possible that Cu substitutes for the Zn, which would result in an inactive enzyme (25). Inactivity could then lead to upregulation via a feedback mechanism (e.g., a mechanism related to the accumulation of substrate). Thus, we speculate that the extreme induction of cadI with high levels of Cu reveals that cofactor substitution is a potential pathway for Cu toxicity in M. tuberculosis.
Other induced enzymes include proteins with functions likely related to oxidative stress (e.g., catalase-peroxidase KatG) and Cu-induced protein damage (e.g., cysteine synthase and thioredoxin). Only two enzymes were repressed in response to Cu, a putative cytochrome c oxidase subunit 4 (Rv2199c) and a putative secreted hydrolase (Rv0988). Because cytochrome c oxidase requires Cu for activity (16), Cu-stimulated regulation is expected. However, the mechanism for regulation of Rv2199c remains unknown. Surprisingly, the predicted multicopper oxidase enzyme (Rv0846c) was not identified as a Cu-responsive protein based on the statistical criteria used in our microarray study, despite the clear role of multicopper oxidases in both E. coli and Salmonella Cu responses (6, 21). Multicopper oxidase has been proposed to reduce Cu(I) to the less toxic form Cu(II), as well as to act as a storage mechanism for excess Cu. Further investigation of the microarray data indicated that Rv0846c transcripts were induced twofold in the presence of a high copper concentration, but induction did not meet the criteria used for statistical significance. However, the induction of Rv0846c with both 50 and 500 µM Cu was clearly demonstrated using qRT-PCR (see Fig. S1 in the supplemental material). Additionally, the protein chaperone DnaK was identified as a Cu-induced protein in the microarray study. In E. coli, DnaK is required to chaperone the multicopper oxidase CueR to the TatAB secretion system (8), which then exports it into the periplasm. The Cu-responsive induction of DnaK and the presence of a predicted Tat signal sequence in Rv0846 (1) suggest that a similar process may occur in M. tuberculosis. Because M. tuberculosis lacks periplasm, we hypothesize that the main role of Rv0846c could be Cu storage and extracellular export.
Copper-responsive transcriptional regulators. As expected, the microarray study identified the previously characterized Cu-responsive regulator CsoR as a protein that is upregulated in response to Cu. Three other transcriptional regulators were also identified, including FurA, an iron-binding regulator associated with oxidative stress (23), and two previously uncharacterized transcriptional regulators, Rv1994c and Rv2642. The Rv1994c protein is a member of the MerR family of transcriptional regulators. MerR regulators are typically metal binding proteins, and the MerR protein CueR is vital for Cu homeostasis in E. coli (9, 18). Similarly, Rv2642 is a member of the ArsR family of metalloregulatory proteins (2). Either of these proteins could complement the Cu-regulatory activities of CsoR, and these proteins are interesting candidates for future study.
In conclusion, preliminary evidence shows that the Cu level fluctuates in the natural environment of M. tuberculosis, and we showed that physiological Cu levels are toxic to the mycobacterial cell. Furthermore, M. tuberculosis transcriptionally regulates a set of genes to deal with the Cu stress both specifically (e.g., csoR, ctpV, and Rv0846) and nonspecifically via the induction of common stress response genes. The Cu-specific genes were found with both 50 and 500 µM Cu, while the stress response genes were found mainly when the 500 µM Cu data set was examined. Our data indicate that CtpV is the most important copper-regulated cation transporter. Together, these data begin to sketch a model for the interplay of Cu-activated genes within M. tuberculosis (Fig. 5). Currently, experiments are under way to elucidate the specific role(s) played by CtpV in Cu homeostasis and M. tuberculosis virulence, as well as to verify aspects of the proposed model. This analysis should further reveal the role of metal-regulated genes in bacterial pathogenesis in general.
![]() View larger version (22K): [in a new window] |
FIG. 5. Model for Cu response in M. tuberculosis. Based on microarray analysis combined with bioinformatics, four major pathways were identified. (i) For Cu export, the predicted Cu-translocating CtpV protein likely functions as a Cu export pump. Additionally, the Rv0849 and Rv2963 proteins were identified as possible nonspecific metal permeases involved in Cu response. (ii) For Cu regulation, CsoR is the Cu-specific regulator of CtpV and other unknown targets. Possible metalloregulators Rv1994c and Rv2642 are also Cu responsive and may contribute to Cu regulation. (iii) For Cu detoxification, the Cu-induced predicted multicopper oxidase Rv0846c may detoxify Cu(I) in the cell and/or may bind and export Cu ions. Its predicted chaperone (DnaK) and export system (Tat) also show Cu induction. (iv) There are also stress response mechanisms for various aspects of Cu-induced stress, including oxidative stress (e.g., FurA/KatG), protein stress (e.g., TrxB), and cofactor substitution (CadI).
|
This work was supported in part by grant NIH-R21AI066235, by Animal Formula funds provided to A.M.T., and by NIH training grant T32GM007215 to S.K.W.
Published ahead of print on 8 February 2008. ![]()
Supplemental material for this article may be found at http://jb.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»