novel pathways associated with bypassing cellular ... · downregulated in immortalized hpecs...

33
Novel Pathways Associated with Bypassing Cellular Senescence in Human Prostate Epithelial Cells Steven R. Schwarze 1 , Samuel E. DePrimo 2 *, Lisa M. Grabert 1 , Vivian X. Fu 1 , James D. Brooks 2 , and David F. Jarrard 1† 1 Department of Surgery, Division of Urology, University of Wisconsin Medical School, Molecular and Environmental Toxicology and the University of Wisconsin Comprehensive Cancer Center, Madison, WI 53972 and 2 Department of Urology, Stanford University, Stanford, CA 94305 *Current address: SUGEN, Inc., South San Francisco, CA 94080 Running title: “GENE EXPRESSION PROFILES IN SENESCENCE AND IMMORTALIZATION” This work was supported by: National Institutes of Health (CA76184-01) (D.F.J), the University of Wisconsin Comprehensive Cancer Center (D.F.J.), the Doris Duke Foundation (J.D.B.), and the Kovitz Foundation (J.D.B.). To whom correspondence should be addressed: David F. Jarrard, M. D. Department of Surgery, University of Wisconsin, 600 Highland Avenue, K6/530, Madison, WI 53792 Ph. 608-265-2225; Fax: 608-265-8133; E-mail: [email protected] The abbreviations used are: SA-ß-gal, senescence-associated ß-galactosidase; HPEC, human prostate epithelial cell; EST, expressed sequence tag; qRT-PCR, quantitative reverse transcriptase-PCR. Copyright 2002 by The American Society for Biochemistry and Molecular Biology, Inc. JBC Papers in Press. Published on February 8, 2002 as Manuscript M200373200 by guest on May 31, 2020 http://www.jbc.org/ Downloaded from

Upload: others

Post on 28-May-2020

8 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

Novel Pathways Associated with Bypassing Cellular Senescence in Human Prostate

Epithelial Cells

Steven R. Schwarze1, Samuel E. DePrimo2*, Lisa M. Grabert1, Vivian X. Fu1, James D. Brooks2,

and David F. Jarrard1†

1Department of Surgery, Division of Urology, University of Wisconsin Medical School,

Molecular and Environmental Toxicology and the University of Wisconsin Comprehensive

Cancer Center, Madison, WI 53972 and 2Department of Urology, Stanford University, Stanford,

CA 94305

*Current address: SUGEN, Inc., South San Francisco, CA 94080

Running title: “GENE EXPRESSION PROFILES IN SENESCENCE AND

IMMORTALIZATION”

This work was supported by: National Institutes of Health (CA76184-01) (D.F.J), the University

of Wisconsin Comprehensive Cancer Center (D.F.J.), the Doris Duke Foundation (J.D.B.), and

the Kovitz Foundation (J.D.B.).

†To whom correspondence should be addressed: David F. Jarrard, M. D. Department of Surgery,

University of Wisconsin, 600 Highland Avenue, K6/530, Madison, WI 53792

Ph. 608-265-2225; Fax: 608-265-8133; E-mail: [email protected]

The abbreviations used are: SA-ß-gal, senescence-associated ß-galactosidase; HPEC, human

prostate epithelial cell; EST, expressed sequence tag; qRT-PCR, quantitative reverse

transcriptase-PCR.

Copyright 2002 by The American Society for Biochemistry and Molecular Biology, Inc.

JBC Papers in Press. Published on February 8, 2002 as Manuscript M200373200 by guest on M

ay 31, 2020http://w

ww

.jbc.org/D

ownloaded from

Page 2: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

Gene expression profiles in senescence and immortalization

2

SUMMARY

Cellular senescence forms a barrier that inhibits the acquisition of an immortal phenotype, a

critical feature in tumorigenesis. The inactivation of multiple pathways that positively regulate

senescence are required for immortalization. To identify these pathways in an unbiased manner,

we performed DNA microarray analyses to assess the expression of 20,000 genes in human

prostate epithelial cells passaged to senescence. These gene expression patterns were then

compared to those of HPECs immortalized with the Human Papilloma Virus 16 E7 oncoprotein.

Senescent cells display gene expression patterns that reflect their non-proliferative, differentiated

phenotype and express secretory proteases and extra-cellular matrix components. A comparison

of genes transcriptionally upregulated in senescence to those whose expression is significantly

downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1,

and a member of the insulin-like growth factor axis, IGFBP-3. Expression of these genes is

found to be uniformly lost in human prostate cancer cell lines and xenografts and previously

their inactivation has been documented in tumor samples. Thus, these genes may function in

novel pathways that regulate senescence and are inactivated during immortalization. These

changes may be critical in not only allowing cells to bypass senescence in vitro, but in the

progression of prostate cancer in vivo.

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from

Page 3: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

Gene expression profiles in senescence and immortalization

3

INTRODUCTION

Normal primary cells proliferate in culture for a limited number of population doublings

prior to undergoing terminal growth arrest and acquiring a senescent phenotype (1). Senescent

cells are resistant to mitogen-induced proliferation, express senescence-associated ß-

galactosidase (SA-ß-gal), and assume a characteristic enlarged, flattened morphology (reviewed

in (2)). With progressive cell divisions, errors in DNA replication, oxidative metabolism and

environmental insults generate somatic mutations and place cells at risk for oncogenic

transformation. Accumulating data suggest that the terminal arrest associated with cellular

senescence represents a major barrier, analogous to apoptosis, that cells must circumvent to

become malignant (3). Research into the pathways that positively regulate senescence and ways

cells bypass senescence is therefore critical in understanding carcinogenesis.

In part, senescence is mediated by loss of telomere length that occurs with progressive

divisions of mortal cells (4). Cellular senescence also can be induced prematurely by a variety of

stimuli including DNA damage, perturbations of chromatin structure, and overexpression of

mitogenic signals including E2F1, oncogenic H-ras, Raf or MEK (5-8). Two genes implicated

in and critical for the induction and maintenance of senescence are pRb and p53 (3).

Overexpression of downstream members of these pathways including p16, the p53 effector

molecule p21, and p14ARF can elicit cell cycle arrest and senescence (5,9,10). These gain of

function experiments demonstrate senescence can be prematurely induced, however do not

directly address the complex pathways regulating senescence in unmanipulated, serially-

passaged, normal primary cells. These mechanisms may significantly differ. For example in

fibroblasts, the overexpression of E2F1 leads to senescence, however E2F1 expression is actually

repressed as unmanipulated cells senesce (11).

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from

Page 4: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

Gene expression profiles in senescence and immortalization

4

Inactivation of the p16/pRb and p53 pathways is necessary to bypass senescence,

however is insufficient for immortalization. Telomerase activation is also required, although the

overexpression of telomerase alone is not sufficient to confer immortalization in epithelial cells

(12). In addition, genetic analyses of immortalized prostate epithelial and other cell types reveal

a number of consistent genetic alterations are required including the amplification of

chromosome 20q and 8q (13,14). Based on these findings, and other microcell transfer

experiments, it has been proposed that alterations in multiple pathways are required for the

acquisition of the immortal phenotype (15). In a series of classical experiments the fusion of

different immortalized cell lines could restore senescence and these cell lines could be sorted into

four complementation groups. This work implies that there are a limited number of pathways

critical to senescence, that these pathways must be activated to initiate senescence, and that they

can override other transforming genetic events to yield a senescent phenotype (16).

Genes that have been implicated in senescence are also targeted in the development and

progression of human prostate and bladder cancer. The p16/pRb pathway is altered in 85% of

primary prostate cancers and predicts adverse clinical outcome (17). The tumor suppressor p53

is one of the most commonly mutated genes in human cancer and is altered frequently in

advanced prostate and bladder cancers (18,19). Several lines of evidence suggest that bypassing

senescence may also be important in cancer progression in vivo. Immortalized colon cancer cell

lines can only be established from large, pathologically advanced lesions (20). When non-

invasive papillary bladder neoplasms are grown in culture, they uniformly senesce and maintain

an intact p16/pRb pathway (19,21). In contrast, cells cultured from myoinvasive, aggressive

transitional cell bladder carcinomas do not senesce, have alterations of the p53 pathway and

inactivate either p16 or pRb. Furthermore, prostate tumors that grow as xenografts have been

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from

Page 5: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

Gene expression profiles in senescence and immortalization

5

generated solely from advanced, metastatic cancers (22). These results suggest that bypassing

senescence is important in cancer progression, and occurs late in the progression of cancers in

vivo.

To gain insight into novel pathways involved in bypassing senescence, we developed a

model in which transcript levels for over 20,000 unique sequences were tracked in replicate cultures

as human prostate epithelial cells senesced or were immortalized. A major advantage of this

approach is that it avoids exogenous gene expression for the induction of premature senescence that

can result in biased or aberrant results. Using this strategy, we have identified several genes, BRAK,

DOC1, and IGFBP-3, which are upregulated in senescence and inactivated during immortalization.

We find utilizing quantitative RT-PCR, that these genes are significantly downregulated in human

prostate cancer cell lines and xenografts. Interestingly, all of these genes have been previously

implicated in the tumorigenic phenotype (23-25). These data support the concept that pathways

activated at senescence and selectively inactivated with immortalization are important in

carcinogenesis and tumor progression. Furthermore, these results provide an important effort in

cataloging genes altered in senescence, a biologically relevant model of aging.

MATERIALS AND METHODS

Cell Culture—Prostate tissue was obtained under an approved IRB protocol from men (ages

44-66) undergoing cystoprostatectomy for bladder cancer at the University of Wisconsin Hospitals

and Clinics. Histology confirmed that no bladder or prostate cancer was present in the tissues

harvested for our studies. Prostate epithelial cultures were established as described previously (13).

Prostate tissues were minced with a scalpel and digested in a solution containing 500 U/ml

collagenase (Sigma) and plated on collagen-coated plates. Cells were maintained in Ham’s F-12

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from

Page 6: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

Gene expression profiles in senescence and immortalization

6

media (Life Technologies, Inc.) supplemented with regular insulin, 0.25 units/ml; hydrocortisone, 1

µg/ml; human transferrin, 5 µg/ml; dextrose, 2.7 mg/ml; non-essential amino acids, 0.1 mM;

penicillin, 100 units/ml; streptomycin,100 µg/ml; L-glutamine, 2 mM; cholera toxin, 10 ng/ml;

bovine pituitary extract, 25 µg/ml and 1% FBS (26). Cells were passaged when confluent after

incubation with trypsin-EDTA. Proliferating cells were harvested 2 days following the first passage.

Pre-senescent cells were harvested when cell growth ceased, 1-5% of all cells appeared

morphologically circular, flattened and enlarged, and SA-ß-gal expression was negative.

Terminally-senescent cells were harvested when >60% appeared morphologically senescent and

expressed high levels of SA-ß-gal.

Transformation of HPECs with HPV16 E6 and/or E7—Infections and characterization of

HPEC cell lines was carried out as described (13). We used a retrovirus construct carrying either the

HPV16 E6 and/or E7 gene(s) as well as a gentamicin resistance cassette (received from Dr. D.

Galloway, Seattle, WA) (27). Sub-confluent proliferating HPECs were infected with 103-105

infectious viral units at pre-passage (~5 X 105 cells per 100-mm dish) in 3 ml of 1% FBS, F-12+

containing 4 µg/ml polybrene (Sigma). After 6 hours, the virus containing media was exchanges and

infected cells were selected with 50 µg/ml of G418 (Life Technologies, Inc.) for a minimum of 7

days. Immortalized cell lines were screened for HPV16 E6 and/or E7 protein expression and loss of

p53 and pRB by western blot analysis. All immortalized lines were passaged in 1% FBS F-12+ for

well over 20 times to confirm their immortality.

Microarray Hybridizations and Data Analysis—Microarray manufacture and

hybridizations were carried out in accord with previously published methods available at (28).

For this set of experiments we used DNA microarrays with either 41,000 or 47,000 spots

representing over 20,000 unique human genes and expressed sequence tags (ESTs). For each

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from

Page 7: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

Gene expression profiles in senescence and immortalization

7

hybridization, two micrograms of poly(A) mRNA from each sample were reverse-transcribed

and labeled with fluorescence-tagged nucleotides (Cy3 for the reference sample, Cy5 for the

experimental sample). Samples were hybridized against a common reference pool of mRNA

derived from a panel of human cell lines, as described previously (29). Hybridizations were

carried out for 16-18 hours at 65°C and the arrays were washed. After drying, the microarrays

were scanned with a confocal laser GenePix microarray scanner (Axon Instruments) and were

analyzed with Genepix software. After visual inspection, spots of poor quality were flagged and

excluded from analysis. Data files containing fluorescence ratios were entered into the Stanford

Microarray Database, and compiled experiments were further analyzed with hierarchical

clustering software and visualized with Treeview software (30). Prior to cluster analysis, gene

entries were filtered to select those with >2.5-fold changes in signal intensity when comparing

proliferating cells to pre-senescent, terminally-senescent and immortalized cell lines on at least

1/3 of the array samples, with fluorescent intensity in each channel that was greater than 1.6

times the local background. Identities of spotted cDNAs are based on information available in

UniGene cluster Build 137 (June, 2001).

Statistical Analysis of Microarray Data—Genes with significant expression changes in

response to senescence and immortalization in normal prostate epithelial cells were identified

using the Significance Analysis of Microarrays (SAM) procedure (31). Changes in gene

expression for any single gene as measured in several array experiments provide a statistically

testable measure of robustness, regardless of the magnitude of change in expression. The SAM

procedure computes a two-sample T-statistic (e.g. proliferating vs. senescent HPECs) for the

normalized log ratios of gene expression levels for each gene. It thresholds the T-statistics to

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from

Page 8: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

Gene expression profiles in senescence and immortalization

8

provide a “significant” gene list and provides and estimate of the false discovery rate (the percent

of genes identified by chance alone) from randomly permuted data.

For the experiments described in our study, the raw expression ratio dataset was filtered,

using the program Cluster, for genes whose transcript levels differ from their median value by at

least 1.5-fold in androgen treated cells compared to controls in at least 2 experiments (with not

more than 30% of measurements discarded due to poor data quality for each entry). We selected

a dataset in which 0.5 is the median number of likely false positive genes (false detection rate of

0.088%) for further analysis. We selected this value because it provided statistical assurance that

most genes in the set were significantly altered, but was not so stringent as to exclude genes of

biological interest.

Quantitative RT-PCR analysis—Total RNA was isolated using the RNeasy RNA

isolation kit (Qiagen), DNAse treated and 1 µg was used to prepare cDNA. Quantitative RT-

PCR was performed by monitoring in real time the increase in fluorescence of the SYBR Green

dye as described using a iCycler detection system (Bio-Rad) (32,33). For comparison of

transcript levels between samples, a standard curve of cycle thresholds for several serial dilutions

of a cDNA sample was established. This value was then used to calculate the relative abundance

of each gene. These values were then normalized to the relative amounts of 18S cDNA. All

PCR reactions were performed in duplicate. Sequences of primers used for PCR analysis are

available upon request.

RESULTS AND DISCUSSION

Analysis of Senescence Pathways in Human Prostate Epithelial Cells—Genome-wide

RNA expression analysis of serially passaged primary human prostate epithelial cells (HPECs)

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from

Page 9: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

Gene expression profiles in senescence and immortalization

9

and HPV16 E7 immortalized HPECs was performed with DNA microarrays representing over

20,000 distinct human genes and ESTs. In order to identify genes induced at senescence and

inactivated at immortalization, HPEC cultures derived from three separate patients were

analyzed at three defined stages in their lifespan. Cells harvested in their first passage and

growing exponentially were defined as proliferating. Pre-senescent cells were harvested when

growth was halted, but the cells had not yet expressed the morphologic characteristics of

senescence and were not yet SA-ß-gal positive. We have previously shown that this is a distinct

point at which cells transiently elevate protein levels of selected cell cycle inhibitory genes, and

it appears critical in signaling the onset of senescence (34). At terminal-senescence, > 60% of

HPECs expressed SA-ß-gal and had acquired an enlarged, flattened morphology (13,34).

Three independent proliferating HPEC cultures were then immortalized utilizing HPV16

E7 retrovirus. We have previously used this molecular tool to dissect the genetic events involved

in bypassing senescence (13,35). The E7 oncoprotein acts, in part, by inactivating the pRb/p16

pathway, and via the binding of p21, the p53 pathway, thus recapitulating several genetic

alterations commonly observed in prostate tumors (13,18). HPV E7 immortalized cells are non-

tumorigenic (36) making them an applicable model with which to isolate genes associated with

the immortalization phenotype. However, because of some of the diverse activities of the HPV

oncoproteins, we note not all targets relevant to bypassing senescence will be identified utilizing

this approach. Poly(A)+ mRNA was isolated from each of these experimental samples, reverse

transcribed, labeled with Cy-5, and hybridized in comparison with a common reference pool of

cDNA (labeled with Cy-3) (29).

Determination of genes that significantly differed in expression levels between cellular

phenotypes was made utilizing the analytical tool Significance Analysis of Microarrays (SAM)

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from

Page 10: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

Gene expression profiles in senescence and immortalization

10

(31). Comparison of proliferating to terminally-senescent HPECs generated 262 unique entries

using SAM, of which 93 were induced at senescence (51 of known function) and 159 were

repressed (77 of known function). Complete details and transcript identities for all outputs can be

accessed in the supplemental data (Supplemental Fig. 1-2). SAM outputs comparing

proliferating to pre-senescent cells identified no genes with significant induction in expression in

pre-senescent cultures, but did find 101 genes whose expression levels were repressed. These

101 genes also displayed repressed transcript levels in terminally-senescent cultures. The results

generated by SAM were filtered to remove ESTs and genes of unknown function. The

remaining genes were then classified based on their cellular function and placed into one or more

related groups.

Consistent with a non-proliferative senescent state, 28 of the 77 known genes that were

repressed at senescence are positive growth regulators or transcription factors involved in promoting

cell cycle progression. An additional 28 repressed genes are directly involved in DNA replication,

repair or mitosis. Of greater interest to us were those genes induced as cells progressed to terminal-

senescence as these genes may lend insights into the induction of senescence. It has been postulated

that p16 and p53 function as tumor suppressors, in part, by inducing senescence pathways (37). We

have previously reported at the cdk inhibitor p16 protein level, a known marker of senescent cells, is

elevated at senescence (34). In this study, we find ~50-fold increase in p16 mRNA at senescence,

suggesting that p16 protein increases are due, at least in part, to increased transcription. Transcript

levels of other negative growth regulators were also increased including the p53-regulated targets

BTG2, p21, IGFBP-3 and mdm2 (Fig. 1A), suggesting elevated p53 activity at senescence.

Overexpression of BTG2 (38) , p21 (39), IGFBP-3 (25) and even mdm2 in normal diploid cells (40)

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from

Page 11: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

Gene expression profiles in senescence and immortalization

11

can result in G0/1 cell cycle arrest. These p53-inducible genes may reflect the participation of p53

in regulating the onset of senescence and lend further insight into senescence-associated p53 targets.

Senescent cells show permanent growth arrest and express markers associated with terminal

differentiation. Consistent with this phenotype, 15 of 51 genes induced in senescent cells are

associated with differentiation, 10 of which are commonly expressed in neuronal tissue (Fig. 1A).

These genes may be indicative of the terminally differentiated phenotype rather than the induction of

a neuronal phenotype as senescent HPECs do not resemble elongated, spindle-shaped cells

characteristic of neuronal cells in culture (13). Senescent HPECs also acquire a phenotype in which

extracellular matrix degrading proteases, growth factors and other inflammatory cytokines are

produced (41). 16 of the 51 genes induced at senescence code for extra-cellular matrix (ECM)

proteins and matrix proteases (Supplemental Fig. 1). This group includes transcripts for

chondroitan sulfated proteoglycan 2, matrix metalloproteinase 2 and micofibrillar-associated protein

2 and their respective increases (5.3, 4.9 and 4.3 fold by cDNA array, respectively) are among the

highest we observed. Several other transcripts encoding ECM proteins (CSPG3, COL1A1, and

fibronectin) and proteases (PAPPA and cathepsin F) are also induced (2.3 and 4.3 fold by cDNA

array, respectively). Senescent cells produce extracellular proteases that promote matrix and

basement membrane degradation, and these changes are thought to contribute to the aging process in

vivo (42). Our data suggest that senescence is associated with significant alterations of the ECM,

and provides some of the molecular underpinnings for these matrix changes.

The finding that the ras homolog gene family member B (RhoB) was induced at senescence

is novel. RhoB mRNA levels were induced 3.6-fold by cDNA array and confirmed with real-time

quantitative PCR (17-fold increase, Fig. 2A). It has been demonstrated that the overexpression of

oncogenic H-ras promotes senescence in primary, mortal epithelial cells and fibroblasts (6);

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from

Page 12: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

Gene expression profiles in senescence and immortalization

12

however we did not observe changes in H-ras expression during senescence in our HPEC expression

profiles. RhoB appears to have several functions in the cell including negative growth regulation,

the induction of cytoskeletal changes, and sensitizing cells to apoptosis (43,44). In aging mice,

RhoB expression was found to be induced in skeletal muscle and this alteration was attenuated with

caloric restriction, the only known mechanism for retarding the aging process (45). Collectively,

these data suggest that RhoB may be integral to both senescence and the aging process, potentially

contributing to both the growth arrest and the characteristic morphology observed in senescent cells.

To validate the alterations in gene expression observed with senescence and immortalization,

we performed real-time quantitative PCR (qRT-PCR) in additional HPEC cultures on 20 genes that

had high (> 3 fold) or low (< 3 fold) level expression changes (Fig. 2). The pattern of gene induction

or repression observed on the microarrays correlated with that seen by qRT-PCR consistent with

other validation experiments performed using this array system (46). However, microarrays tended

to underestimate changes in transcript levels when genes were induced or repressed over ~3 fold.

Many of the gene expression changes we observed are also observed in primary fibroblasts

undergoing senescence (Table 1). However, in contrast to senescent fibroblasts, HPECs show

increased transcript levels of the retinoic acid receptor α, interleukin 1α, interleukin 6, and IGF-I at

terminal-senescence. These differences in expression may be cell-type specific or may reflect

differences between epithelial cells and fibroblasts (29). Comparison of limited expression profiles

from fibroblasts to our data derived from HPECs are extraordinarily consistent suggesting that

pathways that regulate senescence are conserved among cells derived from different tissues.

Identification of pathways and genes altered following HPEC immortalization—To aid

in defining pathways important in bypassing senescence, we immortalized HPEC cell cultures

with the HPV16 E7 oncoprotein. Comparison of genes expressed by these cells to proliferating

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from

Page 13: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

Gene expression profiles in senescence and immortalization

13

non-immortalized cells revealed 229 significant alterations in expression of which 92 were

induced (40 of known function) and 137 repressed (70 of known function) (for a complete list

see Supplemental Fig. 3-4). In bypassing senescence, HPECs appear to lose many specialized

functions associated with proliferating, mortal prostate epithelial cells. Of 70 genes

downregulated upon immortalization, ~29 are associated with an epithelial cell function. For

example, HPV16 E7 immortalized HPECs lose expression of prostate epithelial cell specific

kallikreins 10 and 11, two secreted proteases potentially important in reproduction (47). Not

surprisingly, immortalized HPECs activate genes participating in pathways that promote

proliferation, including transcription factors and signaling proteins (Fig. 1B). In addition, genes

involved in DNA damage response were also upregulated during immortalization (Fig. 1B). It is

unclear whether induction of these repair genes reflects increased DNA surveillance, or an

enhanced DNA repair capacity.

Identification of genes associated with bypassing senescence in prostate cancer

cells—We hypothesized that a gene whose expression is induced during senescence (activation)

and repressed following immortalization (inactivation) could be targeted during prostate cancer

progression. Hierarchical clustering analysis was performed on genes whose expression changed

greater than 1.5-fold as cells progressed to senescence or following immortalization to identify a

group of genes with elevated transcript levels at senescence that returned to levels below or

comparable to normal proliferating cells upon immortalization (Fig. 3). To further identify genes

that may be inactivated during carcinogenesis, we focused on transcripts induced in senescent

HPECs and downregulated in the immortalized HPECs below levels found in HPECs

(determined using SAM). This comparison revealed that only BRAK, downregulated in ovarian

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from

Page 14: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

Gene expression profiles in senescence and immortalization

14

cancer (DOC1), and IGFBP-3, were significantly induced with senescence and repressed in

immortalized cells markedly below (>10-fold) levels in proliferating HPECs.

To extend these findings, we analyzed the expression of BRAK, DOC1, and IGFBP-3 levels

in the human prostate cancer cell lines LNCaP, Du145, TSU, PPC-1, PC-3 and DuPro using qRT-

PCR. In 6/6 tumorigenic, immortalized prostate cancer cell lines we found DOC-1 to be

downregulated roughly 23 fold, and BRAK and IGFBP-3 expression to be undetectable (Fig. 4). We

analyzed the expression levels of these three genes in two additional human prostate cancer

xenografts LAPC4 and LAPC9 to evaluate whether their expression is also lost in vivo. Similar to

the prostate cancer cells, expression of BRAK, DOC1 and IGFBP-3 were all reduced markedly in

these models compared to normal HPECs (Fig. 4). These prostate cancer xenografts have not been

passaged in vitro, suggesting that the decreased expression of BRAK, DOC1 and IGFBP-3 is not due

to artifacts introduced by growth in vitro (22). Our results indicate that genes and pathways

activated at senescence and repressed during HPV16 E7 immortalization are also downregulated in

human prostate cancer cell lines and xenografts and may have a role in human prostate cancer.

Previous reports suggest that BRAK, DOC1, and IGF family members have roles in the

progression of other tumor types in vivo. BRAK is a member of the CXC chemokine family that is

widely expressed in normal tissues, while commonly repressed in head and neck squamous cell

carcinomas (23). Its location on 5q harbors a putative susceptibility gene for prostate cancer based

on familial linkage studies (48). As its name implies, DOC1 was discovered as a transcript

downregulated in ovarian cancers (24). Its function remains unknown. The IGF axis includes

several ligands (IGF-I, IGF-II), cell-surface receptors (IGF-R) and binding proteins (IGFBPs) that

regulate cell growth and differentiation (reviewed in (49)). IGFBP-3 modulates interaction of the

IGFs to their receptors and also possesses IGF-independent functions (50). Decreased IGFBP-3

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from

Page 15: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

Gene expression profiles in senescence and immortalization

15

expression is associated with prostate cancer progression, demonstrating more frequent loss of

expression in advanced disease, in both human and mouse models (25,51,52). Our data suggest that

this may be due to bypassing senescence pathways. Another factor implicating IGFBP-3 in a

senescence pathway is that the HPV16 E7 oncoprotein targets IGFBP-3 for inactivation (53).

Further study of prostate cancer specimens will be necessary to evaluate the role of BRAK and

DOC1 specifically in prostate cancer. These findings underscore the potential importance of

inactivation of senescence in carcinogenesis and tumor progression.

In summary, using DNA microarray analysis, we have identified genes participating in

normal human prostate epithelial cell senescence. We suspect these pathways participate in the

induction of senescence in many other cell types and these genes are members of a limited number

of pathways that can be inactivated to bypass senescence during tumorigenesis (15). We also have

identified novel markers for the senescent phenotype that could assist in an analysis of the role of

senescence in normal and pathologic states in vivo (54). Future studies will explore whether

reactivation of these pathways in immortalized cells can also re-establish cellular senescence and

whether this may represent a novel therapeutic approach to prostate cancer.

Acknowledgements

We would like to thank Dr. Rob Reiter for the xenografts and Dr. John Svaren for the critical

reading of this manuscript.

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from

Page 16: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

Gene expression profiles in senescence and immortalization

16

REFERENCES

1. Hayflick, L. (1965) Exp Cell Res 25, 585-621

2. Campisi, J. (2000) In Vivo 14(1), 183-8

3. Campisi, J. (2001) Trends Cell Biol 11(11), S27-31.

4. Campisi, J., Kim, S., Lim, C. S., and Rubio, M. (2001) Exp Gerontol 36(10), 1619-37.

5. Dimri, G. P., Itahana, K., Acosta, M., and Campisi, J. (2000) Mol Cell Biol 20(1), 273-

85.

6. Serrano, M., Lin, A. W., McCurrach, M. E., Beach, D., and Lowe, S. W. (1997) Cell

88(5), 593-602

7. Lin, A. W., Barradas, M., Stone, J. C., van Aelst, L., Serrano, M., and Lowe, S. W.

(1998) Genes Dev 12(19), 3008-19

8. Zhu, J., Woods, D., McMahon, M., and Bishop, J. M. (1998) Genes Dev 12(19), 2997-

3007.

9. Vogt, M., Haggblom, C., Yeargin, J., Christiansen-Weber, T., and Haas, M. (1998) Cell

Growth Differ 9(2), 139-46

10. Brown, J. P., Wei, W., and Sedivy, J. M. (1997) Science 277(5327), 831-4

11. Dimri, G. P., Hara, E., and Campisi, J. (1994) J Biol Chem 269(23), 16180-6.

12. Kiyono, T., Foster, S. A., Koop, J. I., McDougall, J. K., Galloway, D. A., and

Klingelhutz, A. J. (1998) Nature 396(6706), 84-8.

13. Jarrard, D. F., Sarkar, S., Shi, Y., Yeager, T. R., Magrane, G., Kinoshita, H., Nassif, N.,

Meisner, L., Newton, M. A., Waldman, F. M., and Reznikoff, C. A. (1999) Cancer Res

59(12), 2957-64

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from

Page 17: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

Gene expression profiles in senescence and immortalization

17

14. El Gedaily, A., Bubendorf, L., Willi, N., Fu, W., Richter, J., Moch, H., Mihatsch, M. J.,

Sauter, G., and Gasser, T. C. (2001) Prostate 46(3), 184-90.

15. Sasaki, M., Honda, T., Yamada, H., Wake, N., Barrett, J. C., and Oshimura, M. (1994)

Cancer Res 54(23), 6090-3.

16. Pereira-Smith, O. M., and Smith, J. R. (1988) Proc Natl Acad Sci U S A 85(16), 6042-6.

17. Jarrard, D. F., Modder, J., Sandefur, C., Shi, Y., Heisey, D., Schwarze, S., and Friedl, A.

(submitted)

18. Heidenberg, H. B., Sesterhenn, I. A., Gaddipati, J. P., Weghorst, C. M., Buzard, G. S.,

Moul, J. W., and Srivastava, S. (1995) J Urol 154(2 Pt 1), 414-21.

19. Yeager, T. R., DeVries, S., Jarrard, D. F., Kao, C., Nakada, S. Y., Moon, T. D.,

Bruskewitz, R., Stadler, W. M., Meisner, L. F., Gilchrist, K. W., Newton, M. A.,

Waldman, F. M., and Reznikoff, C. A. (1998) Genes Dev 12(2), 163-74

20. Newbold, R. F., Cuthbert, A. P., Themis, M., Trott, D. A., Blair, A. L., and Li, W. (1993)

Toxicol Lett 67(1-3), 211-30.

21. Sarkar, S., Julicher, K. P., Burger, M. S., Della Valle, V., Larsen, C. J., Yeager, T. R.,

Grossman, T. B., Nickells, R. W., Protzel, C., Jarrard, D. F., and Reznikoff, C. A. (2000)

Cancer Res 60(14), 3862-71.

22. Klein, K. A., Reiter, R. E., Redula, J., Moradi, H., Zhu, X. L., Brothman, A. R., Lamb, D.

J., Marcelli, M., Belldegrun, A., Witte, O. N., and Sawyers, C. L. (1997) Nat Med 3(4),

402-8.

23. Frederick, M. J., Henderson, Y., Xu, X., Deavers, M. T., Sahin, A. A., Wu, H., Lewis, D.

E., El-Naggar, A. K., and Clayman, G. L. (2000) Am J Pathol 156(6), 1937-50.

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from

Page 18: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

Gene expression profiles in senescence and immortalization

18

24. Mok, S. C., Wong, K. K., Chan, R. K., Lau, C. C., Tsao, S. W., Knapp, R. C., and

Berkowitz, R. S. (1994) Gynecol Oncol 52(2), 247-52.

25. Grimberg, A., and Cohen, P. (2000) J Cell Physiol 183(1), 1-9.

26. Reznikoff, C. A., Loretz, L. J., Pesciotta, D. M., Oberley, T. D., and Ignjatovic, M. M.

(1987) J Cell Physiol 131(3), 285-301.

27. Halbert, C. L., Demers, G. W., and Galloway, D. A. (1991) J Virol 65(1), 473-8.

28. DeRisi, J., Penland, L., Brown, P. O., Bittner, M. L., Meltzer, P. S., Ray, M., Chen, Y.,

Su, Y. A., and Trent, J. M. (1996) Nat Genet 14(4), 457-60.

29. Ross, D. T., Scherf, U., Eisen, M. B., Perou, C. M., Rees, C., Spellman, P., Iyer, V.,

Jeffrey, S. S., Van de Rijn, M., Waltham, M., Pergamenschikov, A., Lee, J. C., Lashkari,

D., Shalon, D., Myers, T. G., Weinstein, J. N., Botstein, D., and Brown, P. O. (2000) Nat

Genet 24(3), 227-35.

30. Eisen, M. B., Spellman, P. T., Brown, P. O., and Botstein, D. (1998) Proc Natl Acad Sci

U S A 95(25), 14863-8.

31. Tusher, V. G., Tibshirani, R., and Chu, G. (2001) Proc Natl Acad Sci U S A 98(9), 5116-

21.

32. Wittwer, C. T., Herrmann, M. G., Moss, A. A., and Rasmussen, R. P. (1997)

Biotechniques 22(1), 130-1, 134-8

33. Morrison, T. B., Weis, J. J., and Wittwer, C. T. (1998) Biotechniques 24(6), 954-8, 960,

962

34. Schwarze, S. R., Shi, Y., and Jarrard, D. F. (2002) Oncogene 20, 8184-8192

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from

Page 19: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

Gene expression profiles in senescence and immortalization

19

35. Reznikoff, C. A., Belair, C., Savelieva, E., Zhai, Y., Pfeifer, K., Yeager, T., Thompson,

K. J., DeVries, S., Bindley, C., Newton, M. A., and et al. (1994) Genes Dev 8(18), 2227-

40.

36. Puthenveettil, J. A., Frederickson, S. M., and Reznikoff, C. A. (1996) Oncogene 13(6),

1123-31.

37. Reznikoff, C. A., Yeager, T. R., Belair, C. D., Savelieva, E., Puthenveettil, J. A., and

Stadler, W. M. (1996) Cancer Res 56(13), 2886-90

38. Guardavaccaro, D., Corrente, G., Covone, F., Micheli, L., D'Agnano, I., Starace, G.,

Caruso, M., and Tirone, F. (2000) Mol Cell Biol 20(5), 1797-815.

39. McConnell, B. B., Starborg, M., Brookes, S., and Peters, G. (1998) Curr Biol 8(6), 351-4

40. Brown, D. R., Thomas, C. A., and Deb, S. P. (1998) Embo J 17(9), 2513-25.

41. Krtolica, A., Parrinello, S., Lockett, S., Desprez, P. Y., and Campisi, J. (2001) Proc Natl

Acad Sci U S A 98(21), 12072-7.

42. Campisi, J. (1997) J Am Geriatr Soc 45(4), 482-8.

43. Fritz, G., and Kaina, B. (1997) J Biol Chem 272(49), 30637-44.

44. Liu, A., Cerniglia, G. J., Bernhard, E. J., and Prendergast, G. C. (2001) Proc Natl Acad

Sci U S A 98(11), 6192-7.

45. Lee, C. K., Klopp, R. G., Weindruch, R., and Prolla, T. A. (1999) Science 285(5432),

1390-3.

46. Iyer, V. R., Eisen, M. B., Ross, D. T., Schuler, G., Moore, T., Lee, J. C., Trent, J. M.,

Staudt, L. M., Hudson, J., Jr., Boguski, M. S., Lashkari, D., Shalon, D., Botstein, D., and

Brown, P. O. (1999) Science 283(5398), 83-7.

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from

Page 20: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

Gene expression profiles in senescence and immortalization

20

47. Saitoh, S., Kumamoto, Y., Shimamoto, K., and Iimura, O. (1987) Arch Androl 19(2),

133-47

48. Witte, J. S., Goddard, K. A., Conti, D. V., Elston, R. C., Lin, J., Suarez, B. K., Broman,

K. W., Burmester, J. K., Weber, J. L., and Catalona, W. J. (2000) Am J Hum Genet 67(1),

92-9.

49. O'Dell, S. D., and Day, I. N. (1998) Int J Biochem Cell Biol 30(7), 767-71.

50. Ferry, R. J., Jr., Cerri, R. W., and Cohen, P. (1999) Horm Res 51(2), 53-67

51. Kaplan, P. J., Mohan, S., Cohen, P., Foster, B. A., and Greenberg, N. M. (1999) Cancer

Res 59(9), 2203-9.

52. Hampel, O. Z., Kattan, M. W., Yang, G., Haidacher, S. J., Saleh, G. Y., Thompson, T. C.,

Wheeler, T. M., and Marcelli, M. (1998) J Urol 159(6), 2220-5.

53. Mannhardt, B., Weinzimer, S. A., Wagner, M., Fiedler, M., Cohen, P., Jansen-Durr, P.,

and Zwerschke, W. (2000) Mol Cell Biol 20(17), 6483-95.

54. Chang, B. D., Xuan, Y., Broude, E. V., Zhu, H., Schott, B., Fang, J., and Roninson, I. B.

(1999) Oncogene 18(34), 4808-18

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from

Page 21: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

Gene expression profiles in senescence and immortalization

21

Fig. 1. Identity and cellular role of specific categories of genes induced in terminally-senescent

HPECs (A) and repressed in immortalized HPECs (B). Transcripts of named genes were selected

from the SAM outputs and placed into functional categories. Included in (A) are genes induced at

senescence associated with differentiated cell or neuronal functions and those associated with growth

regulation. Genes repressed in immortalized cells that are involved in signaling, transcription and

DNA repair are listed in (B). The Genbank accession number and fold change from normal

proliferating HPEC proliferating values are listed for each entry (+ indicates gene induction and –

indicates repression).

Fig. 2. Quantitative RT-PCR validation of cDNA microarray results. Real-time PCR following the

incorporation of SYBR green into DNA was performed on cDNA derived from three independent

HPEC cultures at in the proliferating, pre-senescent, terminally-senescent stages and in HPV16 E7

immortalized HPECs. Relative expression levels were normalized to the 18S content, also

determined by RT-PCR, and plotted on a log scale as fold change compared to proliferating cells.

Expression levels from three separate microarray experiments were averaged and plotted with the

qRT-PCR data for comparison. Good correlation was observed for all transcript levels, although the

DNA microarrays tended to under estimate changes in transcript levels at higher magnitudes.

Fig. 3. Cluster analysis of transcripts from senescent and HPV16 E7 immortalized HPECs identified

by the SAM procedure. Experiments are grouped to show transcripts that are induced at pre- and

terminal-senescence and at or below levels of proliferating cells in the immortalized HPECs. Red

squares indicate expression levels higher than that of proliferating HPECs, green squares indicate

expression levels lower than that of proliferating HPECs, black squares indicate transcript levels

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from

Page 22: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

Gene expression profiles in senescence and immortalization

22

approximately equal to those in proliferating HPECs; gray squares indicate data of insufficient

quality. NP, normal prostate; HPVE7, HPV16 E7 immortalized HPECs; prol, proliferating; Pre-

sen., pre-senescence; T-sen., terminally-senescent.

Fig. 4. Expression of BRAK (A), DOC1 (B), and IGFBP-3 (C) is very low in immortalized human

prostate cancer cell lines. qRT-PCR was performed on 6 prostate cancer cell lines (LNCaP, Du145,

TSU, PPC-1, PC-3 and DuPro) and two xenografts (LAPC4 and LAPC9) and normalized to

proliferating HPEC values. Expression levels in all samples were standardized to 18S RNA levels.

Prol, proliferating; Pre-sen., pre-senescent; T-sen., terminally-senescent.

Supplemental Fig. 1. Identity and cellular role of specific categories of genes induced in terminally-

senescent HPECs. Transcripts of named genes were selected from the SAM outputs and placed into

functional categories. The Genbank accession number and fold change from normal proliferating

HPEC values are listed for each entry (+ indicates gene induction).

Supplemental Fig. 2. Identity and cellular role of specific categories of genes repressed in

terminally-senescent HPECs. Transcripts of named genes were selected from the SAM outputs and

placed into functional categories. The Genbank accession number and fold change from normal

proliferating HPEC values are listed for each entry (– indicates gene repression).

Supplemental Fig. 3. Identity and cellular role of specific categories of genes induced in

immortalized HPECs. Transcripts of named genes were selected from the SAM outputs and placed

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from

Page 23: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

Gene expression profiles in senescence and immortalization

23

into functional categories. The Genbank accession number and fold change from normal

proliferating HPEC values are listed for each entry (+ indicates gene induction).

Supplemental Fig. 4. Identity and cellular role of specific categories of genes repressed in

immortalized HPECs. Transcripts of named genes were selected from the SAM outputs and placed

into functional categories. The Genbank accession number and fold change from normal

proliferating HPEC values are listed for each entry (– indicates gene repression).

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from

Page 24: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

Gene expression profiles in senescence and immortalization

24

Table 1. Comparison of HPEC senescence-associated gene mRNA levels

to those reported in the literature (fibroblasts). Fold-change was calculated

by comparing relative transcript levels in proliferating and terminal-senescent

HPECs. Genes are grouped according to the following cellular functions:

(A) Growth regulatory; (B) Cytokine/growth factor; (C) Structural. Values

listed were determined by cDNA microarray. NC, no change.

Gene Reported change HPEC changeId-1 repressed repressed 1.25XId-2 repressed repressed 1.25Xc-fos repressed repressed 1.45XE2F1 repressed repressed 1.5XE2F4 repressed repressed 1.3XE2F5 repressed repressed 1.2Xp33ING induced induced 1.8Xhic-5 induced induced 1.25XDHFR repressed repressed 3.8Xthymidine kinase repressed repressed 9.5Xthymidylate synthase repressed repressed 2.6XDNA polymerase α repressed repressed 1.6Xcyclin A repressed repressed 5.4Xcyclin B repressed repressed 5.2Xcdk2 repressed repressed 2.5Xcdc2 repressed repressed 11.3Xp16INK4a induced induced 49Xp21WAF1/CIP1 induced induced 2.3Xp53 induced induced 2.0XRARα induced NCInterleukin 1α induced repressed 2XInterleukin 6 induced repressed 1.9XIGF-I repressed NCIGFBP-3 induced induced 2.6Xfibronectin 1 induced induced 2.9Xcornifin induced induced 2.4Xcollagenase (MMP2) induced induced 4.9Xstromelysin (MMP3) induced induced 1.6XPAI-1 induced induced 2.5XTIMP1 repressed repressed 4X

C

B

A

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from

Page 25: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

Figure 1

A

Differentiation/neuronalNM_017459 +4.3 microfibrillar-associated protein 2X07868 +4.4 insulin-like growth factor 2NM_006426 +3.5 dihydropyrimidinase-like 4NM_005461 +3.5 Kreisler maf-related leucine zipperNM_002391 +3.3 midkine (neurite growth-promoting)NM_013372 +2.9 melanoma adhesion moleculeAF183421 +2.8 RAB31, ras GTP exchangerNM_013372 +2.7 BMP antagonist 1AL136179 +2.7 SOX4NM_006763 +2.7 BTG2AF231124 +2.7 testicanNM__004112 +2.4 fibroblast growth factor 11NM_032088 +2.4 protocadherin gamma subfamily A, 8NM_004986 +2.3 chondroitan sulfate proteoglycan 3NM_018901 +2.2 protocadherin alpha 5

Growth regulator/tumor associationNM_000077 +49 p16 INK4a

NM_005672 +4.8 prostate stem cell antigenX07868 +4.4 insulin-like growth factor 2NM_014890 +3.2 downregulated in ovarian cancer 1NM_006763 +2.7 BTG2NM_002392 +2.4 MDM2NM_000389 +2.3 p21 WAF1/CIP1

NM_002430 +2.3 meningioma 1

AF144103 +5.3 BRAK

M35878 +2.3 IGFBP-3

B

DNA Repair

subunit 2 homolog

+7.8 proliferating cell nuclear antigenNM_002592+7.2 p53-inducible ribonucleotide small AB036063

+5.7 p53 regulated PA26 nuclear proteinNM_014454+4.5 mdm2NM_002392+4.2 damage-specific DNA binding proteinNM_000107

Signaling

NM_012324 +2.3 JNK-interacting protein 2 (JIP2)

+12 interleukin 1, betaNM_000576+7.2 TNF ligand, member 7NM_001252+6.5 epiregulinNM_001432+6.0 placental growth factorNM_002632+5.4 LRDDNM_018494+4.8 coagulation factor II (thrombin) receptorNM_001992+3.6 ras homolog, family member HNM_004310+3.4 inositol polyphosphate-5-phosphataseNM_005541

+2.4 c-jun amino terminal kinase 1 (JNK1)NM_002750+3.2NM_014288 integrin beta 3 binding protein

Transcription factors+7.6 forkhead box D1NM_004472+7.2 HOXB7NM_004502

S53374 +4.2 transcription factor 19

Gene expression profiles in senescence and immortalization

25

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from

Page 26: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

A

B

C

Figure 2 Gene expression profiles in senescence and immortalization

26

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from

Page 27: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

Figure 3

UteroglobinEVG1 proteinSRY (sex determining region Y)-box 4ESTshypothetical protein FLJ20059matrix metalloproteinase 2cDNA DKFZp56402364cathepsin FCdc42 effector protein 3sema domain, immunoglobulin domainprotocadherin gamma subfamily A, 8matrix Gla proteinhypothetical protein DKFZp434J0617hypothetical protein FLJ14213ESTsDKFZP564I1922 proteinESTsepoxide hyrolase 1ESTsS100 calcium-binding protein A9S100 calcium-binding protein A8STAT4ESTsESTs, similar to p53 regulated PA26-T2aldo-keto reductase family 1, member C2

downregulated in ovarian cancer 1small inducible cytokine subfamily, BRAKnicotinamide N-methyltransferasecytochrome P450, polypeptide

NP

47-P

rol

NP

104-

Pro

lN

P49

-Pro

l

NP

47-P

re-s

en.

NP

104-

Pre

-sen

.N

P49

-Pre

-sen

.

NP

47-T

-sen

.

NP

104-

T-se

n.N

P49

-T-s

en.

HP

VE

7-9

HP

VE

7-15

HP

VE

7-14

P311 proteinATP-binding cassette, sub-family A, 1ras homolog gene family, member B (rhoB)gap junction protein, alpha 1, 43kDmicrofibrillar-associated protein 2milk fat globule-EGF factor 8 proteinchondroitan sulfate proteoglycan 2

Gene expression profiles in senescence and immortalization

27

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from

Page 28: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

-100 -80 -60 -40 -20 1 20

RNA abundance(fold change from prol. HPEC)

HPEC-prol.

HPEC-Pre-sen

HPEC-T-sen.

Du145

DuPro

LNCaP

PC-3

PPC-1

Tsu

LAPC4

LAPC9

HPEC-immortal

BRAKDOC1IGFBP-3

Figure 4 Gene expression profiles in senescence and immortalization

28

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from

Page 29: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

Motility/MigrationNM_004040 +3.6 RhoBNM_004772 +2.6 P311AL136842 +2.4 Cdc42 effector protein 3

Stress/Injury/XenobioticNM_002964 +5.5 S100 calcium-binding protein A8NM_001354 +3.9 aldo-keto reductase family 1, C1NM_002965 +3.7 S100 calcium-binding protein A9NM_001992 +3.0 coagulation factor II (thrombin) receptorNM_005627 +2.9 serum/glucocorticoid regulated kinaseNM_005527 +2.9 heat shock 70kD protein-like 1NM_000120 +2.7 epoxide hydrolase 1, microsomalNM_006169 +2.6 nicotinamide N-methyltransferase

OtherNM_005518 +4.8 HMGCS2-ketone body synthesisNM_000165 +2.8 gap junction protein (connexin 43)U10991 +2.8 G2 protein-glucosidaseNM_005502 +2.7 ABC1-cholesterol metabolism NM_006763 +2.4 milk fat globule-egf factor 8 protein

Supplemental Figure 1

Differentiation/neuronalNM_017459 +4.3 microfibrillar-associated protein 2X07868 +4.4 insulin-like growth factor 2NM_006426 +3.5 dihydropyrimidinase-like 4NM_005461 +3.5 Kreisler maf-related leucine zipperNM_002391 +3.3 midkine (neurite growth-promoting)NM_013372 +2.9 melanoma adhesion moleculeAF183421 +2.8 RAB31, ras GTP exchangerNM_013372 +2.7 BMP antagonist 1AL136179 +2.7 SOX4NM_006763 +2.7 BTG2AF231124 +2.7 testicanNM__004112 +2.4 fibroblast growth factor 11NM_032088 +2.4 protocadherin gamma subfamily A, 8NM_004986 +2.3 chondroitan sulfate proteoglycan 3NM_018901 +2.2 protocadherin alpha 5

Growth regulator/tumor associationNM_000077 +49 p16 INK4a

NM_005672 +4.8 prostate stem cell antigenX07868 +4.4 insulin-like growth factor 2NM_014890 +3.2 downregulated in ovarian cancer 1NM_006763 +2.7 BTG2NM_002392 +2.4 MDM2NM_000389 +2.3 p21 WAF1/CIP1

NM_002430 +2.3 meningioma 1

AF144103 +5.3 BRAK

M35878 +2.3 IGFBP-3

Cytokine/growth factorNM_003357 +5.5 uteroglobin

NM_001874 +3.1 carboxypeptidase M

AF144103 +5.3 BRAKX07868 +4.4 insulin-like growth factor 2NM_002391 +3.3 midkine (neurite growth-promoting)

NM__004112 +2.4 fibroblast growth factor 11

Extracellular matrixU16306 +5.3 chondroitan sulfate proteoglycan 2NM_004530 +4.9 matrix metalloproteinase 2NM_017459 +4.3 microfibrillar-associated protein 2NM_002581 +3.5 pregnancy-associated plasma protein ANM_003793 +3.0 cathepsin FNM_000088 +3.0 collagen, type I, alpha 1X02761 +2.9 fibronectinNM_004986 +2.3 chondroitan sulfate proteoglycan 3

TranscriptionAF055376 +5.9 v-maf musculoaponeurotic fibrosarcomaNM_005461 +3.2 Kreisler maf-related leucine zipper NM_003151 +3.0 STAT4AF129290 +2.8 Cbp/p300 interacting transactivator, 2

NM_000168 +2.2 GLI3AL136179 +2.7 SOX4

Gene expression profiles in senescence and immortalization

29

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from

Page 30: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

DNA Replication-9.5 thymidine kinase 1NM_003258-9.2 ribonucleotide reductase M2NM_001034-6.8 uridine phosphrylaseNM_003364-5.6 CDC45NM_003504-4.7 topoisomerase II, alphaNM_001067

-4.2 DNA polymerase, epsilon 2NM_002692-3.8 dihydrofolate reductaseNM_000791-2.8 replication factor CNM_002915-2.6 primase, polypeptide 1NM_000946

-4.6NM_001254 Cdc6

DNA Repair-4.7 RAD51-interacting proteinNM_006479-4.2 translinNM_004622-3.3 RAD51NM_002875-3.0 BRCA1 associated RING domain 1NM_000465

Transcription factor-5.8 transcription factor 19S53374-4.4 v-myb oncogene homologNM_002466-4.3 forkhead Box M1NM_021953-3.2 pituitary tumor transforming 1NM_004219-3.1 prostate epithelium-specific EtsNM_012391

Apoptosis/Stress Response-5.9 BCL2-related protein A1NM_004049-5.5 transgluaminase 1NM_000359-4.0 catenin, alpha-like 1NM_003798-3.7 alcohol dehydrogenase 4NM_000670-2.8 biliverdin reductase BNM_000713-2.7 tumor rejection antigen 1NM_003299-2.6 CSR1 proteinNM_016240

Positive Growth regulator-11 cdc2NM_001786-11 PDZ-binding kinaseNM_018492-10 CDKN3NM_005192-7.9 polo-like kinaseNM_005030-7.9 HBP17NM_005130-7.3 G0/G1 switchNM_015714-7.2 claspinNM_022111-6.1 RHAMMNM_012484-5.4 cyclin A2NM_001237-4.8 serine/threonine kinase 15NM_003600-4.3 forkhead box M1NM_021953-4.4 v-myb oncogene homologNM_002466-4.2 pro-platelet basic proteinNM_002704-3.9 thyroid hormone receptor interactor 13NM_004237-3.7 anterior gradient 2 homologNM_006408-3.6 CDC28 protein kinase 2NM_001827-3.1 MAPK 13NM_002754-2.8 notch homolog 4NM_004557-3.2 pituitary tumor-transforming 1NM_004219-3.0 insulin-induced gene 1NM_005542-3.0 GRO3 oncogeneNM_002090-2.6 CDC28 protein kinase 1NM_001826-2.2 cdc25ANM_001789

Signaling-4.1 rho GDP dissociation inhibitor betaNM_001175-4.0 catenin, alpha-like 1NM_003798

-2.9 integrin, alpha 4NM_000885-2.9 small inducible cytokine A20NM_004591

-3.9NM_012242 dickkopf homolog 1

Mitosis-11 anillinNM_018685-9.1 NUF2RNM_031423-8.6 ubiquitin protein E2-CNM_007019-7.3 centromere protein FNM_005196-6.8 centromere protein ENM_001813-6.2 NIMA-related kinase 2NM_002497-5.9 BUB1NM_004336-5.9 kinesin-like 5NM_004856-4.2 ZW10 interactorNM_007057-3.2 MAD2-like 1NM_002358-3.0 insulin-induced gene 1NM_005542-2.9 survivinNM_001168-2.8 MCM6NM_005915-2.7 SMC4AL136877 Other

-8.5 follistatinNM_006350-6.9 ADAM8NM_001109-6.7 SLC14A1NM_015865-5.1 L2DTLNM_016448

-4.5 UAP1NM_003115

-4.1 EFEMP1NM_004105

-3.6 diazepam binding inhibitorNM_020548

-3.3 coagulation factor IIINM_001993

-3.1 TROAPNM_005480-3.2 MYT1NM_004203

-3.4NM_004631 apolipoprotein E receptor

-3.6NM_005573 lamin B1

-4.4AB032261 steoaroyl CoA desaturase

-4.2NM_000963 prostaglandin-endoperoxide synthase 2

Supplemental Figure 2 Gene expression profiles in senescence and immortalization

30

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from

Page 31: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

Supplemental Figure 3

Metabolism/Transport+5.1 solute carrier 26NM_000112+5.1 solute carrier 10NM_003049

+4.7 fucosidase, alphaNM_000147+4.6 L2DTL proteinNM_01448+4.0 aldehyde dehydrogenase 7NM_001182

+3.8 ferredoxin reductaseNM_004110+3.8 synoaptosomal-associated proteinAF278704+2.9 sorting nexin 5NM_014426

+4.0NM_000067 carbonic anhydrase II

+4.9NM_004117 FK506-binding protein

Proliferation+7.8 proliferating cell nuclear antigenNM_002592+4.3 Cdc7-like 1NM_003503

+3.2 S-phase kinase-associated protein 2NM_005983+3.1 dihyrofolate reductaseNM_000791

S53374 +4.2 transcription factor 19

Transcription factors+7.6 forkhead box D1NM_004472+7.2 HOXB7NM_004502

S53374 +4.2 transcription factor 19

DNA Repair

subunit 2 homolog

+7.8 proliferating cell nuclear antigenNM_002592+7.2 p53-inducible ribonucleotide small AB036063

+5.7 p53 regulated PA26 nuclear proteinNM_014454+4.5 mdm2NM_002392+4.2 damage-specific DNA binding proteinNM_000107

Growth factor/Cytokine+12 interleukin 1, betaNM_000576

+6.5 epiregulinNM_001432+6.0 placental growth factorNM_002632

Signaling

NM_012324 +2.3 JNK-interacting protein 2 (JIP2)

+12 interleukin 1, betaNM_000576+7.2 TNF ligand, member 7NM_001252+6.5 epiregulinNM_001432+6.0 placental growth factorNM_002632+5.4 LRDDNM_018494+4.8 coagulation factor II (thrombin) receptorNM_001992+3.6 ras homolog, family member HNM_004310+3.4 inositol polyphosphate-5-phosphatase (SHIP)NM_005541

+2.4 c-jun amino terminal kinase 1 (JNK1)NM_002750+3.2NM_014288 integrin beta 3 binding protein

Apoptosis+7.2 TNF ligand, member 7NM_001252+5.4 TNF receptor, member 10c, decoyNM_003841

Other+8.2 fibrillin 2NM_001999+6.3 IGF2 mRNA-binding protein 3NM_006547+4.9 protocadherin, alpha 5NM_018901+3.6 NCAMNM_005010+3.4 NPATNM_002519

+9.1NM_001874 carboxypeptidase M

Gene expression profiles in senescence and immortalization

31

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from

Page 32: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

Supplemental Figure 4

Signaling/growth factor/cytokine-27 small inducible cytokine A2NM_002982

-16 interleukin 11NM_000641-10 TGF§2NM_003238-11 pro-platelet basic proteinNM_002704-8.9 BRAKAF144103-7.1 interleukin 6NM_000600-7.0 IGFBP-3M35878-6.6 TNF ligand superfamily, 10NM_003810-5.6 PDGF§NM_002608-4.9 osteonectinNM_003118-4.8 IGFBP-5NM_000599-4.7 CEACAM6M18216-4.5 integrin, beta 8NM_002214-3.9 CEACAM5NM_004363-3.8 GRO1 oncogeneNM_001511-3.5 WNT5bNM_030775-3.4 interleukin 8NM_000584-3.2 neurotrophin-1NM_013246-3.1 PLAURNM_002659-3.1 guanylate binding protein 1NM_002053-2.8 CTGFNM_001901-2.8 protein kinase H11NM_014365-2.7 enhancer of filamentation 1NM_006404

-21 RAR responder 1NM_002888

ECM-component-58 matrix metalloproteinase 7NM_002423-30 collagen, type III, alpha 1NM_000090-15 collagen type IV, alpha 3NM_004369-10 clade B, member 3NM_006919-14 collagen type IV, alpha 1BC005159-11 hexabrachionNM_002160-7.6 kallikrein 10NM_002776-6.6 collagen type I, alpha 1NM_000088-5.3 tryptase, alphaAK027693-4.5 ADAMTS1NM_006988-4.2 collagen, type V, alpha 2NM_000393-3.4 kallikrein 11NM_006853-3.4 protease, serine, 23NM_007173-3.2 biglycanNM_001711-2.6 thrombospondin 1NM_003246

Cytoskeletal/adhesion-10 keratin 4X07695-12 keratin 13NM_002274-9.5 transgelinNM_003186 -4.9 myosin VBAB032945-4.5 integrin, beta 8NM_002214-4.4 keratin 6ANM_005555-4.2 myosin regulatory light chain 2AC074331-3.2 contactin 1NM_001843

Calcium regulation-6.3 S100 calcium-binding protein PNM_005980-5.6 calmodulin-like 3M36707-3.8 calcium channel, voltage-dependent, beta 2NM_000724-3.6 annexin A6NM_001155-2.9 tumor-associated calcium signal transducer 1 NM_002354

Other-8.1 anterior gradient 2NM_006408-7.0 odd Oz/tenAB032953-5.9 protein kinase inhibitor betaNM_032471-5.5 cerebellar degeneration related proteinNM_004065-4.8 downregulated in ovarian cancer 1NM_014890-4.2 sulfotransferase 2B, member 1NM_004605-3.8 netrin 4AF278532-3.3 solute carrier family 14NM_015865-3.0 transcription elongation factor BNM_003198-2.8 transcript release factorBC004295-2.7 purinergic receptorNM_005767

Stress/Xenobiotic-12 glutathione peroxidase 2NM_002083-6.8 cytochrome p450, subfamily 1NM_000104-6.8 nicotinamide N-methyltransferaseNM_006169-4.3 glutathione S-transferase theta 1NM_000853-3.0 aryl hydrocarbone receptorNM_001621

Gene expression profiles in senescence and immortalization

32

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from

Page 33: Novel Pathways Associated with Bypassing Cellular ... · downregulated in immortalized HPECs identified three genes: the chemokine BRAK, DOC1, and a member of the insulin-like growth

and David F. JarrardSteven R. Schwarze, Samuel E. DePrimo, Lisa M. Grabert, Vivian X. Fu, James D. Brooks

epithelial cellsNovel pathways associated with bypassing cellular senescence in human prostate

published online February 8, 2002J. Biol. Chem. 

  10.1074/jbc.M200373200Access the most updated version of this article at doi:

 Alerts:

  When a correction for this article is posted• 

When this article is cited• 

to choose from all of JBC's e-mail alertsClick here

by guest on May 31, 2020

http://ww

w.jbc.org/

Dow

nloaded from