PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
 
PLoS One. 2010; 5(9): e13013.
Published online 2010 September 24. doi:  10.1371/journal.pone.0013013
PMCID: PMC2945318

Tetracycline Regulator Expression Alters the Transcriptional Program of Mammalian Cells

Andrei L. Gartel, Editor

Abstract

Background

Tetracycline regulated ectopic gene expression is a widely used tool to study gene function. However, the tetracycline regulator (tetR) itself has been reported to cause certain phenotypic changes in mammalian cells. We, therefore, asked whether human myeloid U937 cells expressing the tetR in an autoregulated manner would exhibit alterations in gene expression upon removal of tetracycline.

Methodology/Principal Findings

Microarray analyses revealed that 172 and 774 unique genes were significantly differentially expressed by at least 2- or 1.5-fold, respectively, when tetR expressing U937 cells were maintained in media with or without the antibiotic.

Conclusions/Significance

These alterations in gene expression are likely to contribute to the phenotypic consequences of tetR expression. In addition, they need to be taken into consideration when using the tetR system for the identification of target genes of transcription factors or other genes of interest.

Introduction

Inducible ectopic expression is a widely used tool to study gene function. Its advantages compared to constitutive approaches are that early consequences of the presence of the investigated gene product can be studied, and that it facilitates the analysis also of growth inhibitory or cytotoxic gene products. Initial methods for inducible gene expression were based on mammalian promoters that could be activated, e.g., by metal ions. In the early 1990s, the adaptation of the bacterial tetracycline repressor/operator system for use in mammlian cells provided, for the first time, the possibility to regulate the expression of exogenous genes through a small molecule ligand that per se was expected not to affect the physiology of mammalian cells due to their lack of endogenous receptors for it [1]. In fact, mammalian cells need to be made responsive to tetracycline by introduction of one of several modified versions of the bacterial tetracycline regulator (tetR). The transcriptional activity of some of these derivatives is activated, and that of others repressed, by tetracycline or related compounds [1], [2]. In either case, the expression of genes cloned downstream of a promoter containing several copies of the bacterial tetracycline operator (tetO) can be regulated by removal or addition of tetracycline or tetracycline derivatives from/to the culture media. In further improvements of the system, tetR was placed under the control of tetO elements to create an autoregulatory loop, thus avoiding problems caused by toxicity of constitutive tetR expression and reducing background levels of the gene of interest in the repressed state [3], [4].

Tetracycline regulated gene expression has also been used to study the function of oncogenes involved in myeloid leukemia. A popular model for this purpose is the human cell line U937, which is able to differentiate into the monocytic lineage if stimulated with appropriate agents [5][16]. In the course of such studies, we and others have observed that control cells expressing only a tetR derivative, or containing a tetR derivative plus the tetO vector without a gene of interest, exhibited certain biological responses, e.g. reduced proliferation or differentiation in the absence of exogenous stimuli, to changes in the concentration of tetracycline [5], [6], [10], [12], [15].

Many leukemogenic oncogenes are transcription factors that are either aberrantly expressed or functionally altered as a consequence of chromosome rearrangements. U937 derivatives expressing such oncogenes in a tetracycline regulable manner have also been used in several microarray studies aimed at the identification of their downstream target genes [8], [14], [16]. Since, however, tetracycline removal also modifies the cellular behaviour of control U937 cells containing only the tetR but no gene of interest, we set out to ask whether these phenotypic alterations would be related to transcriptional changes. Indeed, microarray analyses revealed that several hundred genes were deregulated in tetR expressing U937 cells upon removal of tetracycline.

Results

Microarray analysis of U937 cells containing an autoregulated tetR

The cell line U937T_pUHD10S contains a plasmid coding for an autoregulated tetracycline repressor as well as the empty pUHD10S vector, which drives tet-regulated expression of its cDNA inserts. It was generated as a control for cell lines expressing a gene of interest from pUHD10S. Three replicate cultures of U937T_pUHD10S cells were transferred to media with or without tetracycline. 48 h later RNA was extracted and used for hybridization of Affymetrix U133 Plus 2.0 arrays. For statistical evaluation, different levels of stringency were applied so as to reflect stringency levels commonly found in the literature [8], [14], [16]. 1191 probe sets (774 unique genes) were deregulated at least 1.5-fold at p<0.05 upon removal of tetracycline from U937T_pUHD10S cells (Tables 1 and S1A). If the threshold for differential expression was increased to twofold and only genes with a false discovery rate (FDR) <10% after multiple hypothesis testing were considered, 90 probe sets (99 unique genes) resulted as deregulated (Tables 1 and S1B).

Table 1
Numbers of unique genes and probe sets deregulated in U937T_pUHD10S cells upon removal of tetracycline.

Confirmation of the regulation of selected genes by real time quantitative (RTQ-) RT-PCR

Of the genes deregulated in response to tetracycline withdrawal, CD36 (NM_000072; coding for the thrombospondin receptor) and ITGAL (NM_002209; coding for alpha L integrin) were selected for confirmation by RTQ-RT-PCR. Both genes form part of several gene ontology categories that were significantly enriched among genes induced by tetracycline removal (see below). RTQ-RT-PCR indeed corroborated the upregulation of these genes 48 h after removal of tetracycline from the culture media of U937T_pUHD10S cells. Both genes were also induced upon withdrawal of tetracycline from U937T cells, which contain the tetR, but not pUHD10S (Fig. 1A,B). A similar, albeit smaller and less reproducible, effect was already present 24 h after transfer of the cells to tetracycline free media (data not shown). To ask whether the changes in the expression of these genes were related to the induction of the autoregulated tetR, or represented an endogenous response of mammalian cells to altered tetracycline concentrations, native (i.e., untransfected) U937 cells were cultured in the presence of tetracycline for one week. Then they were washed and transferred to media with or without tetracycline in the same manner as U937T_pUHD10S cells. Neither CD36 nor ITGAL was induced upon transfer of native U937 cells to tetracycline free media (Fig. 1C,D).

Figure 1
RTQ-RT-PCR confirms the induction of CD36 and ITGAL upon tetracycline withdrawal from U937 cell derivatives expressing the tetR.

Gene ontology analysis

Gene ontology (GO) analysis revealed the terms "hemopoiesis", "myeloid cell differentiation", "leukocyte activation", "phagocytosis", "response to stress", "defense response", "immune system process", "inflammatory response", "innate immune response", "cell activation", "signal transduction", "chemotaxis", as well as several others, to be significantly enriched among the genes upregulated at least 1.5-fold at an FDR <10% in response to tetracycline withdrawal (Fig. 2, Table S2).

Figure 2
ClueGO analysis [22] of genes upregulated ≥1.5-fold at an FDR <10% in U937T_pUHD10S cells upon tetracycline withdrawal.

Discussion

Because we and others have found that tetR expressing cells show phenotypic changes in response to altered concentrations of antibiotic even if they do not contain a tetracycline regulated gene of interest [5], [6], [10], [12], [15], we asked whether these changes may correspond to altered gene expression patterns. Indeed, using microarray analysis we found that the levels of several hundred mRNAs were altered significantly after transfer of U937T_pUHD10S cells to tetracycline free media. The mechanism through which tetR affects gene expression is unclear and may or may not involve direct interactions with DNA. Irrespectively, it is plausible that the alterations in gene expression patterns in tetR expressing U937 cells contribute to the phenotypic changes observed in these cells after tetracycline withdrawal. Moreover, even though tetR certainly does not play a physiological role in human myeloid cells like U937, several of the gene ontology terms enriched among the genes induced upon tetracycline removal correspond to biological responses that might be expected as a consequence of the experimental expression of a physiologically or pathologically relevant gene. Similarly, altered gene expression patterns may well contribute to recently reported phenotypic consequences of tetR expression in transgenic mice, which were observed even in the absence of any exogenous gene of interest regulated by it [17].

Our results also have practical implications for studies aimed at the identification of target genes of transcription factors or other genes of interest. Because of the still considerable cost of transcriptome analyses, some authors restrict their experiments to the comparison of gene of interest-containing cells maintained in the presence or absence of tetracycline. Factors of 1.5- or 2-fold are often defined as thresholds for differential gene expression, and rigorous statistical analyses including corrections for multiple hypothesis testing are frequently not applied [8], [14], [16] (and according to our experience would in fact prevent further consideration of genes whose differential expression can be confirmed by RTQ-RT-PCR). However, even after multiple hypothesis testing, 99 unique genes were deregulated more than twofold in U937T_pUHD10S cells transferred to tetracycline free media. While this number may appear small at first sight, it is large enough to significantly confound the results of target gene analyses, especially if the investigated gene of interest has few true biological targets. Furthermore, the 774 genes that are deregulated ≥1.5-fold at p<0.05 may also lead to erroneous conclusions, and may do so even in analyses in which higher stringency levels are applied: a gene that would be regulated slightly less than twofold in control cells and slightly above this threshold in gene of interest expressing cells may still not be a relevant biological target of the investigated gene. Therefore, our findings need to be taken into consideration when performing gene expression experiments using the tetR system.

Materials and Methods

Cell lines and cell culture

Human U937 histiocytic lymphoma cells were cultured in a humidified incubator at 37°C and 5% CO2 in RPMI 1640 medium (Invitrogen) containing 10% FBS (Invitrogen). To study the effects of tetracycline removal, 1 µg/ml tetracycline (Sigma) was added to the culture media for one week, and cells were then processed as described below.

U937T cells have been derived from U937 cells by stable transfection with a construct driving tetracycline-regulable expression of the tetVP16 fusion protein [5]. They were cultured in a humidified incubator at 37°C and 5% CO2 in RPMI 1640 medium containing 10% FBS, 0.5 µg/ml puromycin (Sigma), and 1 µg/ml tetracycline. To obtain U937T_pUHD10S cells, U937T cells were electroporated (0.17 kV, 950 µF) with the vector pUHD10S, which contains seven copies of the tetracycline operator, thus facilitating tetracycline regulable expression of its cDNA inserts [18]. Transfectants were selected, cloned, and propagated in RPMI 1640 medium supplemented with 10% FBS, 0.5 µg/ml puromycin, 1 µg/ml tetracycline, and 500 µg/ml hygromycin (PAA). To remove tetracycline from the culture media, exponentially growing cells were washed 3 times with phosphate-buffered saline (PBS) and resuspended in growth media without tetracycline. Control cultures were washed in the same manner but resuspended in media with tetracycline.

Microarray analyses

RNA for microarray analyses was extracted from three independent replicate cultures of U937T_pUHD10S cells using the RNeasy Plus mini kit (Qiagen) according to the manufacturer's instructions. Quality control of the isolated RNA samples, labelling, and hybridization to U133 Plus 2.0 arrays (Affymetrix) were performed at the Center of Excellence for Fluorescent Bioanalytics (KFB) in Regensburg, Germany, which is an Affymetrix Service Provider and Core Facility. Sample preparation was carried out in accordance with the Affymetrix GeneChip Expression Analysis Technical Manual (one replicate) or the Affymetrix GeneChip IVT Express Kit Technical Manual (two replicates). Briefly, total RNA was reverse transcribed into double stranded cDNA, from which Biotin labeled cRNA was generated using the One-Cycle Target Labeling Kit (Affymetrix) or the 3′ IVT Express Labeling Kit (Affymetrix). The length of the cRNA products was assessed using an Agilent 2100 bioanalyzer (Agilent Technologies). Following fragmentation, cRNA products (15 µg) were hybridized to the array for 16 h at 45°C in a rotating chamber. Hybridized arrays were washed and stained in an Affymetrix Washing Station FS450 using Streptavidin-Phycoerythrin conjugate (Life Technologies) together with biotinylated anti-streptavidin antibody (Vector Laboratories), or, for the IVT protocol, preformulated Hyb, Wash, and Stain (HWS) components (Affymetrix). Fluorescent signals were measured with an Affymetrix GeneChip Scanner 3000. Primary data analysis was carried out using the Affymetrix MAS5 algorithm.

Real time quantitative reverse transcriptase polymerase chain reaction (RTQ-RT-PCR)

Total RNA for RTQ-RT-PCR was extracted using Trizol (Invitrogen) and reverse transcribed using random hexamer primers (Invitrogen) and M-MLV reverse transcriptase (Invitrogen) according to the manufacturer's instructions. RTQ-RT-PCR was carried out in an ABI Prism 7700 Sequence Detection System (Applied Biosystems) using the Mesa Green qPCR Master Mix Plus (Eurogentec) according to the manufacturers' instructions. RTQ-RT-PCR primers are shown in Table 2. All assays were carried out in triplicate. Expression values for the genes of interest relative to the housekeeping gene cyclophilinD and to a reference value were determined using the ΔΔCT method [19].

Table 2
Primers used for RTQ-RT-PCR.

Data evaluation and in silico analyses

Only probe sets with present calls in all six samples were considered for statistical analyses, which were performed on log2-transformed intensity values using the R/Bioconductor package limma employing a moderated t-test. Probe sets whose average expression differed at least 1.5- or 2-fold at p<0.05 between cells maintained in the presence or absence of tetracycline were considered as differentially expressed. Where indicated, multiple hypothesis testing according to the method of Benjamini and Hochberg [20] was performed, and probe sets with a false discovery rate (FDR) <10% were considered significantly deregulated. Differentially expressed unique genes were identified in the same manner, except that probe sets with no gene annotation or with gene names including the term “hypothetical” were omitted, and for probe sets sharing the same gene annotation only the most informative ones (i.e., those with the highest interquartile range) were kept. All analyses were performed using the statistical software environment R 2.10.1.

Over-represented gene ontology (GO) terms for up- or down-regulated genes were identified using DAVID [21] based on gene symbols. Up-regulated genes were functionally grouped into gene ontology networks using the Cytoscape plug-in ClueGO [22].

Supporting Information

Table S1

Genes deregulated at least 1.5-fold at p<0.05 (Table S1A) or at an FDR <10% (Table S1B) upon withdrawal of tetracycline from U937T_pUHD10S cells. FC, fold change; FDR, false discovery rate.

(0.19 MB XLS)

Table S2

Gene ontology terms significantly overrepresented among genes deregulated ≥1.5-fold at an FDR <10% upon withdrawal of tetracycline from U937T_pUHD10S cells. Count, number of genes upregulated or downregulated in respone to tetracycline withdrawal that are assigned to the named GO term; %, percentage of genes upregulated or downregulated in respone to tetracycline withdrawal that are assigned to the named GO term. FDR, false discovery rate.

(0.04 MB XLS)

Acknowledgments

U937T cells and the pUHD10S vector were kindly provided by Dr Gerard Grosveld, St Jude Children's Research Hospital, Memphis, TN, USA. Microarray hybridizations were performed at the Center of Excellence for Fluorescent Bioanalytics (KFB) in Regensburg, Germany. The authors wish to thank Dr Thomas Stempfl of the KFB for advice and analyses that went far beyond strict fulfillment of the business contract.

Footnotes

Competing Interests: The authors have declared that no competing interests exist.

Funding: This work was supported by Austrian Science Foundation (FWF), http://www.fwf.ac.at/, grant no P20920 to RW. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

References

1. Gossen M, Bujard H. Tight control of gene expression in mammalian cells by tetracycline-responsive promoters. Proc Natl Acad Sci U S A. 1992;89:5547–5551. [PubMed]
2. Gossen M, Freundlieb S, Bender G, Muller G, Hillen W, et al. Transcriptional activation by tetracyclines in mammalian cells. Science. 1995;268:1766–1769. [PubMed]
3. Shockett P, Difilippantonio M, Hellman N, Schatz D. A modified tetracycline-regulated system provides autoregulatory, inducible gene expression in cultured cells and transgenic mice. Proc Natl Acad Sci U S A. 1995;92:6522–6526. [PubMed]
4. Vignali D, Carson R, Chang B, Mittler R, Strominger J. The two membrane proximal domains of CD4 interact with the T cell receptor. J Exp Med. 1996;183:2097–2107. [PMC free article] [PubMed]
5. Boer J, Bonten-Surtel J, Grosveld G. Overexpression of the nucleoporin CAN/NUP214 induces growth arrest, nucleocytoplasmic transport defects, and apoptosis. Mol Cell Biol. 1998;18:1236–1247. [PMC free article] [PubMed]
6. Bromleigh V, Freedman L. p21 is a transcriptional target of HOXA10 in differentiating myelomonocytic cells. Genes Dev. 2000;14:2581–2586. [PubMed]
7. Caslini C, Shilatifard A, Yang L, Hess J. The amino terminus of the mixed lineage leukemia protein (MLL) promotes cell cycle arrest and monocytic differentiation. Proc Natl Acad Sci U S A. 2000;97:2797–2802. [PubMed]
8. McConnell M, Chevallier N, Berkofsky-Fessler W, Giltnane J, Malani R, et al. Growth Suppression by Acute Promyelocytic Leukemia-Associated Protein PLZF Is Mediated by Repression of c-myc Expression. Mol Cell Biol. 2003;23:9375–9388. [PMC free article] [PubMed]
9. Caslini C, Serna A, Rossi V, Introna M, Biondi A. Modulation of cell cycle by graded expression of MLL-AF4 fusion oncoprotein. Leukemia. 2004;18:1064–1071. [PubMed]
10. Hakansson P, Lassen C, Olofsson T, Baldetorp B, Karlsson A, et al. Establishment and phenotypic characterization of human U937 cells with inducible P210 BCR/ABL expression reveals upregulation of CEACAM1 (CD66a). Leukemia. 2004;18:538–547. [PubMed]
11. Hakansson P, Segal D, Lassen C, Gullberg U, Morse H, 3rd, et al. Identification of genes differentially regulated by the P210 BCR/ABL1 fusion oncogene using cDNA microarrays. Exp Hematol. 2004;32:476–482. [PubMed]
12. Hultquist A, Cetinkaya C, Wu S, Castell A, Erlandsson A, et al. Mad 1 inhibits cell growth and proliferation but does not promote differentiation or overall survival in human U-937 monoblasts. Mol Cancer Res. 2004;2:464–476. [PubMed]
13. Liang H, Samanta S, Nagarajan L. SSBP2, a candidate tumor suppressor gene, induces growth arrest and differentiation of myeloid leukemia cells. Oncogene. 2005;24:2625–2634. [PubMed]
14. Meester-Smoor M, Janssen M, Grosveld G, de Klein A, van IW, et al. MN1 affects expression of genes involved in hematopoiesis and can enhance as well as inhibit RAR/RXR-induced gene expression. Carcinogenesis. 2008;29:2025–2034. [PubMed]
15. Konrad T, Karger A, Hackl H, Schwarzinger I, Herbacek I, et al. Inducible expression of EVI1 in human myeloid cells causes phenotypes consistent with its role in myelodysplastic syndromes. J Leukoc Biol. 2009;86:813–822. [PMC free article] [PubMed]
16. Rice K, Hormaeche I, Doulatov S, Flatow J, Grimwade D, et al. Comprehensive genomic screens identify a role for PLZF-RARalpha as a positive regulator of cell proliferation via direct regulation of c-MYC. Blood. 2009;114:5499–5511. [PubMed]
17. Morimoto M, Kopan R. rtTA toxicity limits the usefulness of the SP-C-rtTA transgenic mouse. Dev Biol. 2009;325:171–178. [PMC free article] [PubMed]
18. Fornerod M, Boer J, van Baal S, Jaegle M, von Lindern M, et al. Relocation of the carboxyterminal part of CAN from the nuclear envelope to the nucleus as a result of leukemia-specific chromosome rearrangements. Oncogene. 1995;10:1739–1748. [PubMed]
19. Livak K, Schmittgen T. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25:402–408. [PubMed]
20. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Statist Soc B. 1995;57:289–300.
21. Huang da W, Sherman B, Lempicki R. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57. [PubMed]
22. Bindea G, Mlecnik B, Hackl H, Charoentong P, Tosolini M, et al. ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics. 2009;25:1091–1093. [PMC free article] [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science