Previous Article | Next Article ![]()
Molecular and Cellular Biology, October 2008, p. 5951-5964, Vol. 28, No. 19
0270-7306/08/$08.00+0 doi:10.1128/MCB.00305-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Maria F. Pino,
Luo Jia Tang, and
Jennifer A. Pietenpol*
Department of Biochemistry, Center in Molecular Toxicology, Vanderbilt-Ingram Cancer Center, Vanderbilt University School of Medicine, Nashville, Tennessee 37232
Received 22 February 2008/ Returned for modification 17 May 2008/ Accepted 27 July 2008
|
|
|---|
|
|
|---|
The integral role of p53 in tumor suppression has prompted many laboratories to perform extensive analyses of signaling pathways downstream of the p53 family of sequence-specific DNA binding transcription factors (p53 and its homologs p63 and p73). It is estimated that p53 has 1,600 binding sites in the human genome, only 22% of which are near promoters (19, 104). A more recent study has identified
5,800 binding sites for p63 (106). By using integrative genomic tools, hundreds of novel target genes have been identified for all three family members (37, 44, 77, 100, 104, 106). Similar analyses with many transcription factors have led to an explosion of genomic binding site and target gene expression data (42). These data sets hold great potential for much more than the characterization of downstream pathways; in particular we predict that they can be used to define the signaling pathways that reside upstream of transcription factors of interest.
Despite the ability of p73 to regulate many p53 family target genes, little is known about the specific pathways that modulate p73 during development, tumorigenesis, and tumor therapy (67). Unlike p53, which is mutated in more than 50% of human cancers, p73 is not mutated during tumorigenesis but instead can be overexpressed (16, 28, 29, 110, 111). There has been much interest in modulating p73, due to its high expression level in p53-deficient tumors and its ability to activate p53 target genes, leading to apoptosis of tumor cells (2, 9, 22, 36, 59, 68). Given the above, drug-inducible pathways upstream of p73 are of therapeutic interest.
We used gene expression signatures downstream of p73 to identify novel upstream regulators. A p73 gene signature was created using a combination of genomic tools and was queried against a database of drug-related profiles known as the Connectivity Map (48, 61, 105). Pattern-matching software (61) was used to identify potential p73-activating drugs. A link between p73 and mammalian Target of Rapamycin (mTOR), a kinase important in energy homeostasis and tumorigenesis, was identified and validated, demonstrating the utility of this approach to identify critical signaling nodes upstream of transcription factors.
|
|
|---|
Rapamycin (Calbiochem), metformin (Sigma), and pyrvinium (United States Pharmacopeial Convention) were used as described. For experiments involving rapamycin, cells were plated at 3 x 105 to 4 x 105 cells per 10-cm2 dish. (HMECs were plated at 5 x 105 cells per 10-cm2 dish.) After cells had completely attached, medium was changed 12 h prior to addition of drug (to avoid experimental variation due to the effect of medium replacement on mTOR [20]), or drug was added to serum-free medium as described. Medium without antibiotics was used for treatment.
For cell growth experiments, MDA-MB-231 cells were plated in triplicate and treated with short hairpin RNA (shRNA)-expressing lentivirus, as described below. After 2 days, cells were treated with 20 nM rapamycin or vehicle control and total cell numbers were determined at the indicated times.
The status of p63 in breast cancer cell lines was confirmed using adenoviral shRNA kindly provided by L. Ellisen (18).
Cell transfection/infection and shRNA. The following sequences were used for small interfering RNA: p73-1, 5'-GCAATAATCTCTCGCAGTA-3'; p73-2, 5'-GAGACGAGGACACGTACTA-3'; green fluorescent protein (GFP), 5'-GAAGGTGATACCCTTGTTA-3'; mTOR (FRAP1), 5'-GCATTTACTGCTGCCTCCTAT-3'; p73β, 5'-TCAAGGAGGAGTTCACGGA-3'; and TAp73, 5'-GAACCAGACAGCACCTACT-3'. 293T cells were transfected using Fugene 6 (Roche, Indianapolis, IN). For knockdown of p73, the pSicoR lentivirus system was used (99). Production of virus and transduction were performed as described previously (82). For knockdown of mTOR (FRAP1), the pGIPZ system was used according to the manufacturer's protocol (OpenBiosystems, Huntsville, AL).
For microarray and ChIP assays, MDA-MB-231 cells were infected with adenovirus expressing hemagglutinin-TAp73β (pAdEasy-1:HA-TAp73β) or with a control adenovirus, and the cells were harvested after 80% transduction efficiency was reached, as monitored by GFP fluorescence. The recombinant adenoviruses were generated as previously described (43), pAdEasy was kindly provided by B. Vogelstein (Johns Hopkins University). Overexpression of multiple p73 isoforms was performed using the pcDNA3 backbone (kindly provided by C. Backendorf and G. Melino [24, 73]), and transfection was performed using Lipofectamine (Invitrogen, Carlsbad, CA) as previously described (78).
Protein lysate preparation and Western analysis.
Cells were trypsinized and lysed as previously described (92). Western analysis was done as previously described (35) with the following antibodies: p73 monoclonal antibodies IMG-246, IMG-259, and IMG-313 (Imgenex, San Diego, CA); p73 monoclonal antibody cocktail Ab-4 (Neomarkers, Fremont, CA); mdm2 monoclonal antibody SMP14, β-actin polyclonal antibody I-19, mTOR polyclonal antibody N-19 (
-FRAP), p63 monoclonal antibody 4A4, and p53 monoclonal antibody DO-1 (Santa Cruz Biotechnology, Santa Cruz, CA); glyceraldehyde-3-phosphate dehydrogenase (GAPDH) monoclonal antibody MAB374 (Chemicon, Temecula, CA); p21 monoclonal antibody Ab-1 (Calbiochem, San Diego, CA); phospho-4EBP1 Thr37/46 polyclonal antibody, phospho-S6 Ser235/236 polyclonal antibody 2F9, and total S6 monoclonal antibody 54D2 (Cell Signaling Technology, Danvers, MA); MAP1LC3B antibody (Abgent, San Diego, CA); and p73 antibody (Bethyl Laboratories, Montgomery, TX). p73 was immunoprecipitated for ChIP with Ab-4 or p73 antibody using conditions previously described (92). A Fluor-S Max MultiImager (Bio-Rad, Hercules, CA) was used to quantify Western signals.
Flow cytometry. Flow cytometry was performed as described elsewhere (6, 12).
RNA isolation, microarray experiments, and statistical analyses. Microarray experiments were performed in duplicate as follows: H1299 cells were infected with adenoviruses expressing GFP or TAp73β for 5 h, and RNA was isolated using the Aurum total RNA minikit (Bio-Rad, Richmond, CA) and submitted to the Vanderbilt-Ingram Cancer Center Microarray Shared Resource for quality control. The RNA was processed, and the microarray was hybridized as previously described (7). Microarray data analyses were performed using the ArrayAssist software platform (Stratagene, La Jolla, CA). Algorithms similar to those described elsewhere (7) were used to create a list of probes with changes in gene expression for p73-overexpressing samples versus GFP controls. The following software programs were used for statistical analyses, gene annotations, and determination of categorical enrichment as indicated: ArrayAssist (Stratagene), WebGestalt (Bioinformatics Resource Center at Vanderbilt University) (114), Ingenuity Pathway Analysis (Ingenuity Systems, Redwood City, CA), NCBI DAVID, and the Connectivity Map (61). KEGG and gene ontology analyses was accessed through WebGestalt, using statistical tests coupled to the WebGestalt interface (114).
Quantitative reverse transcription-PCR. Total RNA was purified and reverse transcribed, and quantitative real-time PCR was performed as previously described (7); primer sequences are available upon request.
ChIP and ChIPSeq. Formaldehyde cross-linking and chromatin preparation and immunoprecipitation (ChIP) were carried out as described previously (44, 90, 92). For ChIPSeq and semiquantitative ChIP experiments, cells were cross-linked and submitted to GenPathway, Inc. (San Diego, CA), according to their FactorPath protocol.
Microarray and ChIPSeq data set accession numbers. All gene expression data from microarray experiments and genomic location information from the ChIPSeq experiment have been deposited in the National Center for Biotechnology Information Gene Expression Omnibus (GEO; http://www.ncbi.nlm.nih.gov/geo/) and are accessible through GEO series accession numbers GSE11626 and GSE11672.
|
|
|---|
) had been overexpressed (1). The lists of genes both induced and repressed by these transcription factors were queried against the Connectivity Map resource (61). In brief, this resource consists of pattern-matching software that compares an input gene signature to a database of signatures from 164 small-molecule bioactive compounds (dubbed perturbagens), 85 of which are classified as pharmaceutical drugs. A connectivity score from +1 to –1 is assigned based on the degree of similarity or dissimilarity between the two signatures (61). Thus, a drug with a high connectivity score has a gene signature very similar to the query signature and might be hypothesized to act on a pathway in parallel with the transcription factor that generated the query signature. This score allows the user to choose perturbagens irrespective of P value, which is of particular relevance because perturbagens are profiled at different doses, in different cell lines, or for a different number of experimental replicates (61). Both the average connectivity score and the maximum connectivity score (from the best instance of treatment, dose, and cell line) are informative; for example, some cell lines might not express the target of interest, and some doses might not be effective, bringing down the average.
In this manner we evaluated the gene signatures of p53 and PPAR
that we had created to test the feasibility of our approach. Many of the well-studied chemical agents that activate p53 are not included in the Connectivity Map. Nevertheless, when we analyzed the p53 gene signature, two known activators of p53 were identified: nocodazole, a microtubule inhibitor (23), and tioguanine, a chemotherapy drug known to induce p53-mediated autophagy (95, 113). These agents were ranked 6th and 18th, respectively, out of all 85 pharmaceutical perturbagens by average connectivity score (data not shown). Analysis of the PPAR
signature resulted in the identification of the PPAR
activator troglitazone, used to treat diabetes mellitus type 2. Troglitazone was ranked 22nd out of 85 compounds (data not shown). Given these results, we considered the 30 highest-ranking perturbagens to be of likely relevance in terms of modulating a transcription factor of interest.
Having ascertained the feasibility of the in silico arm of our approach, we sought to identify novel signaling upstream of the transcription factor p73. To establish a collection of p73-regulated genes, Affymetrix GeneChips were used to quantify transcript levels after ectopic p73 expression in the H1299 lung carcinoma cell line. This cell line does not have readily detectable expression of any of the p53 family members (65), making it ideal for an analysis of p73 without the confounding effect of its homologs. Since p73 can be expressed as multiple isoforms, we ectopically expressed an isoform designated TAp73β that is capable of strong transactivation (24). p73 transcriptional activity was evidenced by induction of known target genes p21 and mdm2 (Fig. 1A to C), and 11 out of the top 20 regulated genes as shown by microarray are known p53 family targets (Table 1).
![]() View larger version (28K): [in a new window] |
FIG. 1. Generation of a multitiered p73 signature. H1299 cells were transduced with TAp73β- or GFP-expressing adenoviruses. (A) Protein lysates were harvested after transduction, and p73, GAPDH, and downstream targets mdm2 and p21 were analyzed by Western blotting. (B and C) p73 regulates known target genes when expressed in H1299 cells. (B) Total RNA was purified and reverse transcribed, and quantitative real-time PCR was performed with primers for p21 and mdm2. The samples were normalized to GAPDH, and the results are presented as increases over values for GFP control. Error bars represent standard deviations from three experiments. (C) For ChIP analysis, p73 was immunoprecipitated (IP) from formaldehyde-cross-linked H1299 cells transduced with adenoviral p73 or a GFP control. Associated DNA fragments were PCR amplified using primers flanking the p53 family response elements in p21 and mdm2. Nonspecific binding was assayed by immunoprecipitation (IP) of p73 from non-cross-linked lysates or cross-linked lysates with an isotype-matched antibody (–, specific immunoprecipitation). (D) Schematic showing the number of genes that increase or decrease upon p73 overexpression relative to GFP control by microarray alone, the number of genes identified by ChIPSeq analysis alone, and the number of genes that are present in both data sets. (E and F) DNA fragments were created from ChIP as for panel D, and analysis of p73 binding at genomic regions near the indicated genes was performed by semiquantitative PCR for genes showing high levels of binding (E) and lower levels of binding (F). "Neg." represents a negative-control region. Error bars represent standard deviations from three experiments.
|
|
View this table: [in a new window] |
TABLE 1. Top 20 genes identified by microarray analysis
|
Perturbagens identified through the Connectivity Map increase p73 levels. Multiple tiers of information were used to create a rank-ordered list of target genes for further analysis, including (i) presence in microarray data set, (ii) change in expression level with p73 overexpression, (iii) presence in ChIP data set, (iv) presence of a cluster of ChIP tags at gene locus, (v) number of tags per cluster, (vi) cluster length, and (vii) functional annotation by a variety of methods that are described in greater detail below. The microarray data set was used to query the Connectivity Map because it contains a large number of genes analyzed in the proper format. However, we later used all tiers of information to choose target genes for follow-up analysis.
To begin our analysis, genes whose transcript levels increased or decreased twofold after ectopic p73 expression were analyzed using the Connectivity Map. Because we were ultimately interested in connecting p73 to pathways rather than just drugs per se, we focused on the "pharmaceutical" subset of 85 perturbagens, which we found to be better affiliated with known molecular targets and signaling pathways. Five of the top 30 perturbagens predicted to induce p73 activity were either direct or indirect inhibitors of mTOR signaling (Table 2). This includes sirolimus (known as rapamycin), a drug that binds to FKBP12 to form a complex that inhibits mTOR (40), as well as metformin, a drug widely used to treat diabetes mellitus type 2 that has been shown to activate AMP kinase and inhibit mTOR in cell culture (112). We also identified phenformin, a drug that is in the same class as metformin but is no longer used therapeutically; pyrvinium, an inhibitor of the Akt kinase that is an upstream activator of mTOR (32); and dexamethasone, known to inhibit mTOR in muscle cells (103).
|
View this table: [in a new window] |
TABLE 2. Top 30 pharmaceutical perturbagens identified through the Connectivity Map that induce a p73 signature
|
![]() View larger version (23K): [in a new window] |
FIG. 2. Enrichment of genes by function and signaling pathway in the p73 gene signature. (A and B) Enrichment of major biological processes among genes regulated by p73. Gene ontology enrichment is shown for sets of genes that are both present in the ChIP data set and increased (A) or decreased (B) twofold over GFP values with p73 overexpression in H1299 cells by microarray analysis. Processes with P values by hypergeometric test of less than 0.01 and containing two or more genes as annotated by WebGestalt are graphed. (C) Analysis of KEGG signaling pathways enriched among all genes that were upregulated twofold or more in p73-overexpressing H1299 cells by microarray. Enrichment is shown as the number of observed genes in the data set compared to the expected number of genes as calculated using the WebGestalt software.
|
and TAp73β, that differ only in the presence or absence of a C-terminal sterile alpha motif of unknown function (26). Only the TAp73β isoform is elevated after treatment with rapamycin and metformin in Rh30 cells (Fig. 3B). Pyrvinium pyruvate inhibits mTOR through the upstream kinase Akt (32). In a breast cancer cell line expressing high levels of Akt (MDA-MB-468) (109), we also observed an increase in p73 protein levels 12 h after treatment with pyrvinium (Fig. 3C). Thus, drugs that inhibit mTOR by multiple mechanisms were found to elevate p73 levels.
![]() View larger version (35K): [in a new window] |
FIG. 3. Western analysis of perturbagen effect on p73. (A) p73 levels were increased in MDA-MB-231 cells treated with rapamycin (rap) (left panel) or metformin (met) (right panel) for 24 h. (B) Rh30 cells treated with rapamycin (left panel) or metformin (right panel) for 24 h. (C) MDA-MB-468 cells treated with pyrvinium (pyr) for 36 h. (D) Serum starvation enhances rapamycin-induced regulation of p73. Left panel: 20 nM rapamycin was added to MDA-MB-231 cells 12 h after replacement of medium containing 10% serum. Right panel: 20 nM rapamycin was added to MDA-MB-231 cells in fresh serum-free medium. For all panels, "C" is vehicle control. Protein lysates were analyzed by Western blotting for p73, p4EBP1, pS6, total S6, and actin as indicated. Panels are representative of at least three independent experiments.
|
![]() View larger version (22K): [in a new window] |
FIG. 4. Changes in p73 protein levels do not correspond to changes in p73 RNA levels. MDA-MB-231 cells (A) and Rh30 cells (B) were treated with rapamycin for 24 h in 10% serum and analyzed by Western blotting for p73, pS6, total S6, and actin. Total RNA was purified 24 h after treatment and reverse transcribed, and quantitative real-time PCR was performed with primers for TAp73. The samples were normalized to GAPDH, and the results are presented as increase over vehicle control (C). Error bars represent standard deviations from three experiments. Densitometry was performed on p73 Western signals, followed by normalization to actin. The increases in TAp73β protein levels over vehicle control were 2.2-fold in panel A and 4.3-fold in panel B.
|
![]() View larger version (29K): [in a new window] |
FIG. 5. General cell cycle inhibition does not increase p73 levels. MDA-MB-231 and Rh30 cells were treated with 100 nM hydroxyurea (Hu) or 500 µM mimosine (Mim) or the appropriate vehicle control (C) for 24 h. Protein lysates were harvested and analyzed for p73 and actin by Western blotting, and cells were analyzed by flow cytometry to assess cell cycle profile.
|
![]() View larger version (64K): [in a new window] |
FIG. 6. Differential regulation of p53 family members by rapamycin. HMECs were treated with vehicle control (C) or rapamcyin for 12 h. Protein lysates were analyzed by Western blotting for p53, p63, p73, actin, pS6, and total S6. Results are representative of three independent experiments.
|
![]() View larger version (23K): [in a new window] |
FIG. 7. mTOR regulates p73 levels and activity. (A) mTOR knockdown induces p73. MDA-MB-231 cells were transduced with lentivirus engineered to express shRNA against the FRAP1 subunit of mTOR or with control lentivirus. Protein lysates were harvested 3 days after transduction, and mTOR, p73, p4EBP1, and actin were analyzed by Western blotting to demonstrate knockdown of mTOR levels and activity and induction of p73 levels. Western blots are representative of at least three independent experiments. (B) MDA-MB-231 cells were transduced with lentivirus engineered to express shRNA against GFP or p73. Protein lysates were harvested 5 days after transduction, and reduction of p73 levels was confirmed by Western blotting. (C) p73 activity is induced in MDA-MB-231 cells treated with 20 nM rapamycin, with or without concurrent serum starvation, and the result was verified using p73 RNAi. Cells were serum starved by preincubation in serum-free medium overnight before treatment with rapamycin. p73 RNAi was performed by transducing cells with lentivirus engineered to express shRNA against either GFP or p73 72 h before treatment. Total RNA was purified 48 h after treatment and reverse transcribed, and quantitative real-time PCR was performed with primers for the indicated genes. The samples were normalized to GAPDH, and the results are presented as increase over vehicle control values for an average of three experiments. Samples that exhibited a 30% or greater increase relative to control are indicated in red. Twelve of 17 genes exhibited a p73-dependent increase in RNA levels after rapamycin treatment and serum starvation. (D) Change in RNA as in panel C for INSR, TSC1, and XDH shows rapamycin/serum starvation-induced changes that are p73 dependent. Error bars represent standard deviations for three experiments. (E) Semiquantitative ChIP was performed to assess levels of p73 binding to genomic regions in INSR, TSC1, and XDH promoters or introns in Rh30 cells treated with vehicle or 40 nM rapamycin for 24 h. (F) MDA-MB-231 cells were transduced with shRNA lentivirus as for panel B and treated with 20 nM rapamycin or vehicle. At the indicated times cells from treated or control cultures were counted, and changes in cell growth rates due to p73 and rapamycin are shown. Error bars represent standard deviations from three experiments.
|
Three target genes further demonstrate our selection methods. The insulin receptor (INSR) was identified by both microarray analysis and ChIP and was selected based on known cellular function. Similarly, tuberous sclerosis 1 (TSC1) was identified by ChIP and is a known component of the mTOR signaling pathway (34). Finally, xanthine dehydrogenase (XDH) is an example of a target gene chosen both for its ontology and for its presence as a cluster of tags in the ChIPSeq data set (data not shown). XDH regulates cellular metabolism and cellular response to reactive oxidative stress (3), a process that can by regulated by mTOR during hypoxia (27).
XDH, TSC1, and INSR are three examples of a larger subset of genes whose regulation by mTOR is p73 dependent at the transcriptional level. An increase in all three genes after treatment with rapamycin in the absence of serum is abrogated by p73 knockdown, as depicted in Fig. 7C and in greater detail in Fig. 7D. In addition, semiquantitative ChIP in Rh30 cells revealed an increase in p73 binding to promoter or intronic sequences in all three genes in response to rapamycin (Fig. 7E). These genes, which were selected using multiple criteria, thus serve as readout of an mTOR-p73 signaling axis.
Also consistent with mTOR being an upstream modifier of p73 signaling, knockdown of p73 prevented a decrease in cell growth rate in MDA-MB-231 cells after treatment with rapamycin (Fig. 7F). These data suggest that cellular processes affecting growth rate are modulated by mTOR and p73.
To assess more-specific functions of p73 that might be regulated by mTOR, we focused on autophagy, a form of degradative cell death in response to energy starvation recently shown to be induced by p73 (21). Because our gene signature was generated using the TAp73β isoform and only TAp73β is induced by rapamycin, we tested the effect of selective knockdown of this isoform on autophagy by using two different RNA interference (RNAi) constructs targeting either the N-terminal TA domain or the C-terminal β domain. In both cases, knockdown abrogated baseline autophagy as measured by detection of the cleaved forms of MAP1LC3 (LC-I and LC-II) by Western blotting (Fig. 8A). MAP1LC3-II is the only protein known to be associated with the completed autophagolysosome and is considered a marker of autophagy (83).
![]() View larger version (24K): [in a new window] |
FIG. 8. Analysis of p73-regulated genes in profiling studies of starvation and starvation-induced autophagy. (A) p73β knockdown decreases levels of autophagy markers. MDA-MB-231 cells were transduced with lentivirus engineered to express shRNA against TAp73 isoforms or against p73β isoforms. Protein lysates were prepared, and p73, actin, and the autophagy markers LC3-I and LC3-II were detected by Western blotting. (B and C) Genes from the p73 signature were assessed using publicly available data sets in which T98G glioblastoma cells were grown asynchronously or serum starved for 3 days before RNA harvest and microarray analysis (15) (B) or in which Awells B-lymphoblastoid cells were serum starved for 6 h or 24 h, inducing autophagy (25) (C). Known p53, p63, and/or p73 target genes indicated in orange are DDIT4 (31), DDB2 (53), DFNA5 (66), CDKN1C (10), GADD45A (49), JAG1 (88), and SESN2 (13). Cell lines are arranged in columns, grouped by treatment as indicated. Genes (annotated from Affymetrix probes) are in rows that are ordered based on hierarchical clustering results (data not shown). Color range shown indicates baseline transformed expression level on a log scale. Genes that are upregulated during autophagy or serum starvation that were selected for additional analysis are in blue. (D) Changes in RNA as in Fig. 4 for KLHL24 and LOC153222 indicate increases in RNA levels after rapamycin treatment in a p73-dependent manner in MDA-MB-231 cells.
|
|
|
|---|
Determinants of our approach include overexpression of the protein of interest in a cell line that does not express high levels of the protein. This should be easily applied to transcription factors that are not essential for cell cycle progression or proliferation. In addition, we predict that modifications of our approach such as using a knockdown instead of an overexpression technique would be successful. One of the striking findings from the completion of the Connectivity Map was the extent to which connections could be made across different cell lines and tissues (60, 61). Likewise, we did not encounter any problems due to the cell lines that we chose for analysis, although this might be a significant issue for a modified knockdown approach. In our study, p73 was regulated by mTOR inhibitors in primary cells as well as multiple cancer cell lines.
The link between mTOR and p73 sheds new light on p73 biology, further demonstrating that p73 can be regulated by upstream pathways outside of canonical p53 signaling (55, 58, 64, 85, 96, 97, 102). For example, it has been shown that Akt can regulate p73 and p73-dependent apoptosis (2, 8). It would be intriguing to determine the role of mTOR (and additional members of the mTOR signaling pathway such as Akt) in canonical p53 family functions such as apoptosis and cell cycle arrest. Because Akt has been shown to regulate an activator of p73, YAP, directly (8), this mechanism may function in parallel with mTOR to affect p73-dependent transcriptional regulation.
Other functions downstream of mTOR may also be mediated by p73. Recent studies show that metformin induces apoptosis specifically in p53-null tumor cells (14). Given our results showing that metformin can increase p73 levels, additional studies are ongoing to test the role of p73 in metformin-induced apoptosis and tumor toxicity, in the presence or absence of p53.
It is the isoform TAp73β that is induced by rapamycin and metformin, with apparent coordinate downregulation of TAp73
(Fig. 3B and 4B). Few proteins that differentially modify the alpha versus beta isoforms of p73 have been identified (50, 72). We identified a putative E3 ubiquitin ligase in our data sets, TRIM32. The RNA levels of TRIM32 increase 142% upon serum-free treatment with rapamycin (Fig. 7C and data not shown). The increase is decreased to 37% after treatment with concomitant p73 RNAi-mediated knockdown (Fig. 7C and data not shown). Perhaps TRIM32 or another regulatory protein differentially modifies TAp73
and TAp73β isoforms in response to signaling from the mTOR pathway.
We demonstrated that endogenous regulation of caspase-independent cell death, known as autophagy, is a specific function of the TAp73β isoform (Fig. 8A). Although the mechanism that links p73 to autophagy is unknown (21), the p73 signature presented in this work contains candidate target genes, including unknown factors such as KLHL24 and LOC153222, that are increased in autophagy by microarray profiling (Fig. 8). KLHL24 contains a BTB/POZ domain that can function in transcriptional repression, and LOC153222 contains a basic leucine zipper domain suggesting transcription factor activity (4). These factors may regulate a transcriptional program inhibited by mTOR but induced by p73 during autophagy.
Our study suggests a larger role for mTOR in regulating the p53 family as a whole. Recent studies have placed p53 upstream of mTOR through regulation of PTEN, TSC1, and IGFBP3 (33), as well as downstream of mTOR in vivo in hamartomas (62). In general the mTOR signaling pathway integrates multiple inputs and can involve alteration of and feedback regulation by Akt and the mTORC2 complex (86, 87, 91, 101). It would be interesting to determine the signaling mechanisms by which p73 connects the p53 family to this pathway. Our data suggest that p73 is regulated by mTOR at the posttranscriptional level, either by altering protein stability or by other mechanisms. For example, recent studies suggest that p53 family transcripts may have internal ribosomal entry sites that would allow for cap-independent translation (80, 89, 107). Cap-independent translation may play a role in mTOR regulation of the p53 family.
As has been suggested by Golub, Lamb, and colleagues, our study demonstrates the potential impact of an expanded Connectivity Map as a community-wide resource that would include many more drugs and perturbagens than were included in the first build (60, 61). Here we show utility for linking the pathways represented by the Connectivity Map to the study of transcription factors in particular, but the power of this approach is dependent on the number of drugs and drug-inducible pathways that are testable through this resource. In our study we chose to assess the top 30 compounds predicted to induce a p73 gene signature. This was based on our pilot study using transcription factors with known activators. In contrast, few compounds were identified that would be predicted to repress p73 (i.e., compounds with a negative connectivity score), and we did not do a formal analysis of these drugs. However, one compound with a negative connectivity score was estradiol (data not shown), a known activator of the mTOR pathway (108). Preliminary data suggest that estradiol can decrease p73 levels (data not shown). Therefore, we predict that using our approach to identify compounds that inhibit transcription factors could also elucidate signaling connections. An expanded Connectivity Map that included more drugs that impact the same pathway at multiple levels would allow for more rigorous statistical testing and thus increase an investigator's ability to look for enrichment of pathways and agents among the top hits.
Already, genomic technologies have created a wealth of information downstream of transcription factors. We propose that this information can be used not only to characterize downstream signaling pathways but also to map upstream drug-inducible pathways. Using a genomic-based screening procedure and cell culture-based validation techniques, we identified a novel link between p73 and mTOR, an important kinase in energy homeostasis and tumorigenesis. Because p73 isoforms are overexpressed in many tumors, they may provide effective targets for cancer treatment (9). Recent studies suggest that p73 may be present in an inactive form in select head and neck squamous cell carcinomas and breast tumors and that activation of p73 would lead to tumor cell apoptosis (28, 63, 81). Thus, our results have implications for ongoing and future clinical trials examining the efficacy of mTOR inhibitors in tumors expressing p73.
This work was supported by the National Institutes of Health grants CA70856 and CA105436 (J. A. Pietenpol), ES00267 and CA68485 (core services), US Army Medical Research and Materiel Command grant W81XWH-04-1-0304 (J. M. Rosenbluth), and GM07347 (MSTP training).
Published ahead of print on 4 August 2008. ![]()
These authors contributed equally to this work. ![]()
|
|
|---|
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»