Distinct gene-expression profiles characterize mammary tumors developed in transgenic mice expressing constitutively active and C-terminally truncated variants of STAT5
© Eilon and Barash; licensee BioMed Central Ltd. 2009
Received: 04 December 2008
Accepted: 18 May 2009
Published: 18 May 2009
Stat5 is a latent transcription factor that regulates essential growth and survival functions in normal cells. Constitutive activity of Stat5 and the involvement of its C-terminally truncated variant have been implicated in blood cell malignancies and mammary or breast cancer. To distinguish the individual contributions of the Stat5 variants to mammary tumorigenesis, global gene-expression profiling was performed on transgenic STAT5-induced tumors.
We identified 364 genes exhibiting differential expression in mammary tumors developed in transgenic mice expressing constitutively active STAT5 (STAT5ca) vs. its C-terminally truncated variant (STAT5Δ750). These genes mediate established Stat5 effects on cellular processes such as proliferation and cell death, as well as yet-unrelated homeostatic features, e.g. carbohydrate metabolism. A set of 14 genes linked STAT5Δ750 expression to the poorly differentiated carcinoma phenotype and STAT5ca to the highly differentiated papillary adenocarcinoma.
Specifically affected genes exhibited differential expression in an individual tumor set vs. its counterpart and the intact mammary gland: 50 genes were specifically affected by STAT5ca, and 94% of these were downregulated, the latter involved in suppression of tumor suppressors and proliferation antagonistics. This substantial downregulation distinguishes the STAT5ca-induced tumorigenic consequences from the relatively equal effect of the STAT5Δ750 on gene expression, which included significant elevation in the expression of oncogenes and growth mediators.
STAT5Δ750 mRNA expression was below detection levels in the tumors and the amount of STAT5ca transcript was not correlated with the expression of its specifically affected genes. Interestingly, we identified several groups of three to eight genes affected by a particular STAT5 variant with significant correlated expression at distinct locations in the clustergram.
The different gene-expression profiles in mammary tumors caused by the STAT5Δ750 and STAT5ca variants, corroborated by the absence of a direct link to transgenic STAT5 expression, imply distinct metabolic consequences for their oncogenic role which probably initiate early in tumor development. Tumorigenesis may involve induction of growth factor and oncogenes by STAT5Δ750 or suppression of tumor suppressors and growth antagonists by STAT5ca. The list of genes specifically affected by the STAT5 variants may provide a basis for the development of a marker set for their distinct oncogenic role.
Stat5 is a latent transcription factor that is activated upon binding of specific cytokines to their cognate membranal receptors. Activation of Stat5 involves phosphorylation of specific tyrosine residues by Janus kinase 2 (Jak2), and the translocation of dimerized Stat5 molecules into the nucleus to bind target sites (TTCC(A>T)GGAA) in individual gene promoters. Consequently, induction of a wide variety of signaling events ensues, observed mainly in the hematopoietic system and mammary gland. Stat5's transactivation domain (TAD) has been mapped to amino acids 750 to 772 in the C-terminal domain of the molecule . This domain interacts with co-factors and is essential for transcription activation . The TAD encompasses the main diversity between STAT5a, which is prevalent in the mammary gland, and Stat5b which is more abundant in the liver. Several naturally occurring C-terminally truncated variants of Stat5 have been identified. These variants are generated by alternative splicing or proteolytic cleavage, most commonly by nuclear serine proteases [3–5].
In the blood, naive T cells in the peripheral blood mononuclear cell fraction exclusively express the C-terminally truncated Stat5 . Upon activation by mitogenic stimuli, the truncated Stat5a and Stat5b are replaced by the full-length Stat5, implying that the truncated proteins have distinct functions. Naturally truncated forms of Stat5 have been implicated in blood-cell cancer: 94% of patients with relapsed leukemia expressed this Stat5 variant, suggesting that it controls progression of the disease . Constitutive activation of Stat5 has also been linked to a variety of blood-derived malignancies such as BCR-ABL-induced chronic myeloid leukemia (CML ), acute myeloid leukemia (AML) and acute lymphoid leukemia (ALL ).
In the mammary gland, Stat5's expression and activity are induced during pregnancy and lactation and decrease upon involution. During these stages, Stat5 controls epithelial cell proliferation, final differentiation, lactogenesis, cell survival and tissue remodeling. These cellular processes play a major role in the structural and functional adaptation of the gland to the specific stages of the female reproductive cycle (reviewed in ). Overexpression of STAT5 in transgenic mice induced proliferation of mammary epithelial cells during pregnancy, increased β-casein synthesis upon lactation, and delayed involution . In contrast, the C-terminally truncated variant was unable to induce β-casein/luciferase activity upon prolactin stimulation in vitro. Its expression in the mammary glands of transgenic animals resulted in reduced rates of cell proliferation at pregnancy and increased apoptosis during involution. Morphological signs of milk secretion upon parturition were delayed .
Stat5 has also been associated with breast cancer. Stat5a nuclear localization was observed in 76% of breast cancer specimens and a positive correlation was established between its nuclear localization and the level of histological differentiation of the tumors . Inactivation of Stat5a in transgenic mice expressing transforming growth factor (TGF) α or the SV40 T antigen delayed hyperplasia  and mammary cancer progression. More recently, a direct effect for Stat5 on mammary tumorigenesis was established [15, 16]. Overexpression and forced activation of STAT5 caused parity-dependent development, most frequently of differentiated tumors, in transgenic post-estropausal female mice. Surprisingly, comparable rates of tumors (~8%) with similar latency periods were monitored in mice expressing the C-terminally truncated STAT5 protein, though the amount of poorly differentiated tumors in these mice was higher.
Both the constitutively active and C-terminally truncated variants of Stat5 are potent oncogenes. In this study, we sought to determine the distinctness of their effects and the exclusive contribution of each variant to the oncogenic profiles of gene expression.
Profiling global gene expression in breast cancer has improved our understanding of the clinical diversity of this disease, allowing a better classification of its subtypes and definition of their response to drug treatment. Attempts to predict survival rates have also been reported [17–20]. Transgenic mouse models have been used to study the signatures of specific genes involved in initiation and maintenance of the disease. Genetic analysis of these mouse models could be highly relevant to human cancer [21–23]. The data also suggest that genes involved in the same pathway generate tumors with similar expression profiles, which are distinct from the profiles of tumors arising from other transgenic pathways .
Gene-expression profiles were compared in tumors caused by the constitutively active STAT5 (STAT5ca) and its C-terminally truncated variant (STAT5Δ750). A set of 364 differentially expressed genes was identified. Analysis and classification of these genes suggest that defined differences, possibly triggered at early stages of tumor development, characterize these Stat5 variants.
Mouse mammary-tumor samples
Mammary tumors were derived from transgenic mice carrying one of two Stat5 variants on a FVB/N background: (i) constitutively activated STAT5, termed STAT5ca, comprising sequences from three genes: amino acids 1–750 from ovine Stat5, which is homologous to mouse Stat5a, 677–847 from human Stat6, and 757–1129 from mouse Jak2 and (ii) a deleted construct, STAT5Δ750, prepared by introducing a stop codon at the respective site of the native Stat5 DNA sequence, thus eliminating the expression of its TAD. These constructs were inserted into the β-lactoglobulin (BLG) multiple-cloning site for mammary-gland-specific expression . Upon identification, tumors were excised and snap-frozen for RNA isolation and validation studies. The pathological analysis of the tumors was performed by Dr. Robert Cardiff (University of California, Davis) as previously described . All animals used in this study were treated humanely. Study protocols were in compliance with the regulations of the Israeli Ministry of Health and local institutional policies (approval no. IL- 39-03).
RNA extraction and microarray hybridization
The protocols for RNA extraction from tumors and mammary glands and microarray hybridization were as previously described . Briefly, RNA was extracted from individual tumors or mammary glands with TRIZOL and reverse-transcribed. Equal amounts of complementary (c) RNA from each tumor were hybridized to an Affymetrix GeneChip® Mouse Genome 430A 2.0 array (Affymetrix, Santa Clara, CA), which includes approximately 14,000 annotated genes from the mouse genome. Hybridization and signal quantitation were performed according to Affymetrix's protocol by the Biological Services of the Weizmann Institute of Science (Rehovot, Israel). Total RNA (15 μg) was reversed-transcribed using a T7-oligo(dT) promoter-primer in the first-strand DNA-synthesis reaction. Following RNase H-mediated second-strand cDNA synthesis, the double-stranded cDNA was purified and used as a template for the subsequent in-vitro transcription reaction. This reaction was carried out in the presence of T7-RNA polymerase and a biotinylated nucleotide analogue/ribonucleotide mix for cRNA amplification and biotin labeling. The biotinylated cRNA targets were then cleaned up, fragmented, and hybridized to the GeneChip expression array. The chip was reacted with streptavidin-phycoerythrin and then with biotinylated anti-streptavidin antibody (Vector Laboratories, Burlingame, CA). Arrays were scanned by GeneArray scanner G2500A (Hewlett Packard, Palo Alto, CA), visually inspected for hybridization imperfections and analyzed using Affymetrix Microarray Suite software version 5.0 by scaling to an average intensity of 250. The raw data have been deposited in the public repository "Gene Expression Omnibus" (GEO), accession no. GSE15119.
Statistical, hierarchical clustering and functional annotation analyses
The data were analyzed with GeneSpring (Silicon Genetics, Redwood City, CA) using the MAS5 algorithm . Gene-expression data were normalized "per chip" and "per gene". For "per chip" normalization, all expression data on a chip were normalized to the 50th percentile of the measurements taken from all values on that chip. For "per gene" normalization, each gene's measurement in the selected samples was divided by the median of the gene's measurements in the respective control group, according to the type of comparison being made. When profiles were compared between tumors induced by each of the two transgenic Stat5 variants, the expression of a given gene was normalized to the median of the expression level in the wild-type mammary gland samples. This allowed relating the level of expression of a given gene to that in the mammary gland . When the analysis also included values obtained from the mammary glands, the expression of a given gene was normalized to the median of the expression levels of all genes from all samples. The normalized data were log-transformed and the differences in gene expression (based on the individual values obtained from each tumor) were calculated using one-way statistical analysis of variance (ANOVA). The statistical analysis, which discriminates between the effects of transgenic STAT5 variants, was cross-validated by the K-nearest-neighbor algorithm using the leave-one-out methodology . The differences in gene expression between the tumors and the mammary gland were determined by post-hoc Student-Newman- Keuls analysis. Note that the different expression patterns of selected genes in the array had been previously confirmed by semi-quantitative PCR analyses of selected genes, and hybridization to the respective probes .
Genes exhibiting a significant (P ≤ 0.05) twofold difference in expression between transgenic variants were categorized into "biological functions and/or diseases" using Ingenuity Pathways Analysis (IPA) software (Ingenuity Systems, Mountain View, CA). The probability of each term being identified by random chance was calculated using Fisher's exact test http://www.ingenuity.com/. Unsupervised hierarchical clustering (GeneSpring) organized these genes according to their similarity or dissimilarity in expression profiles, placing the cases with similar expression profiles together as neighboring rows in the clustergram.
Multivariate correlations were determined between the expressions of transgenic STAT5 and specifically affected genes, and among the expressions of the affected genes in the individual tumors using JMP statistical software (Version 7.0, SAS Institute, Inc., Cary, NC).
Quantitative real-time PCR analyses were performed in an ABI Prism 7700 (Applied Biosystems, Foster City, CA) in a 20-μl reaction volume containing 4 μl cDNA (diluted 1:100), 10 μl SYBR Green PCR Master Mix (Applied Biosystems) and 10 μM primers. The thermal-cycle conditions consisted of 2 min at 50°C, 2 min at 95°C, and 40 cycles of 15 s at 95°C and 1 min at 60°C. The primers were designed so that the PCR would yield a single product without any primer dimerization, and the product was verified using a dissociation protocol [see Additional file 1, data sheet A]. The primers were designed across exon-exon junctions to ensure that there was no genomic DNA contamination in the cDNA samples. The amplification curves for the selected genes (70–100 bp) were parallel.
Establishment of distinct expression profiles for the effects of STAT5ca and STAT5Δ750 in mammary tumors
IPA software was used to annotate the differentially expressed genes into processes involved in cellular metabolism and cancer (Table 1). In both categories, the highest number of differentially expressed genes in the STAT5ca- vs. STAT5Δ750-induced tumors was associated with two opposing processes: cell proliferation and cell death [see Additional file 1, data sheet C]. The involvement of some of the genes affecting cancer cell death–Adm, Cflar, Cyr61, Ddit4, Itgb1, Mapk1, Mapk8, Tgfβ2, Tpm1, Vegfα and Wasf1, was demonstrated in breast cancer cell lines [see Additional file 1, data sheet C]. The activities of Itgb1 (coding for integrin β1), TGFβ2 and Vegfα are mediated, at least in part, by Cyr61 (CCN1, [28–30])–a secreted matrix protein involved in the clinical progression of breast cancer to an invasive phenotype .
Cellular metabolism and cancer processes affected by genes that are differentially expressed in tumors caused by the different transgenic STAT5 variants.
No. of Molecules
formation of inositol phosphate
production of carbohydrate
transport of carbohydrate
utilization of carbohydrate
Cardiovascular System Development and Function
angiogenesis of cells
cardiovascular process of cornea
G1 phase of eukaryotic cells
G2 phase of tumor cell lines
interphase of eukaryotic cells
length of telomeres
cell death of cell lines
cell death of tumor cell lines
Cellular Assembly and Organization
formation of cellular protrusions
formation of filopodia
formation of plasma membrane projections
growth of plasma membrane projections
Cellular Growth and Proliferation
formation of eukaryotic cells
growth of cell lines
proliferation of eukaryotic cells
cell movement of eukaryotic cells
cell movement of tumor cell lines
homing of eukaryotic cells
migration of tumor cell lines
DNA Replication, Recombination and Repair
synthesis of DNA
Hematological System Development and Function
Nervous System Development and Function
growth of neurites
contraction of tissue
carcinoma in situ
carcinoma in situ
apoptosis of tumor cell lines
cell death of tumor cell lines
cell death of breast cancer cell lines
cell movement of tumor cell lines
growth of tumor
growth of tumor cell lines
migration of tumor cell lines
proliferation of tumor cell lines
transformation of cells
The differential expression of genes mediating carbohydrate metabolism indicates involvement of the STAT5 variants in cellular homeostasis as well. Most of the genes involved in the transport and utilization of carbohydrates were more highly expressed in tumors caused by the BLG/STAT5Δ750 transgene, for example, Aqp7 and Aqp9 which code for proteins operating as glycerol channels , or glucokinase (Gck) which encodes a protein that catalyzes the conversion of glucose to glucose-6-phosphate, thus maintaining glucose homeostasis .
Hierarchical clustering of genes affected by transgenic STAT5 variants
Annotation of the 76 genes in cluster 1 assembled them into non-malignant cellular processes, including amino acid modification, hematopoiesis and neurological disorders [see Additional file 1, data sheet E]. In contrast, all other clusters which included genes expressed at higher levels in the STAT5Δ750-induced tumors contained the "cancer" category which assembled genes into specific tumorigenic processes. In this category, genes mediating "invasion", such as Adm, Itg3 and Igf2, converged to cluster 2. In contrast, most genes involved in tumor growth, transformation and differentiation were expressed at lower levels and were therefore allocated to cluster 3.
Transgenic STAT5 variants exert distinct effects on a limited set of specifically affected genes
Classification, annotation and expression levels of genes, specifically affected by the STAT5 variants.
Specifically affected by BLG/STAT5ca
Ras small GTPase
Rasd1 Rasd2 GNAO1
Tnni1, Cnn1, Rsnl2
Organization of microtubules
Msln, ptn, Rapgef3, Adra2c
Csf2ra, il5ra, Foxji, H2-K1, EDG3
Elavl3, Wnt8a, Gfra2, Phgdh
DNA damage response/cell cycle
Prp19 (DDR), Magea2 (inhibits p53 activation)
JIP-2, Rapgef3, Mapk8
Mapk1, JIP-2, Rapgef3, Mapk8, Gfra2, Adra2c, Rasd1, Rasd2, GNAO1, Chn1, EDG3, Wnt8a, JIP-2
Protein phosphates inhibitor activity
Regulation of nucleotide and nucleic acid metabolism
Bacterial cell wall degradation
Regulating hemoglobin oxygen affinity
Expressed in the early embryo
Tyrp1, 1700065I17Rik, Serpinb6c, 2010003K11Rik, 4931412G03Rik (Trpd52l3)
Specifically affected by BLG/STAT5Δ750
Adm, Mrpplf4, Cyr61
Cspg2, Itgβ1, Itgβ3, Scgb1a1, Pcdhb5, Rin1, Lana5, Cyr61, Cav2
Itgβ1, Foxk1, Rnf134, Peg3, Mkrn3
Progression through cell cycle
Peg3, Kif5c, Pcdhb5, Ivns1abp, Mkrn3
Foxk1, Rnf134, Peg3
Transcription factor activity
GTPase activator activity
Tec, Trhr, Oprl, V1ra9, 5930418K15Rik
G-protein-coupled receptor protein signaling pathway
Adm, Wnt7b, Pcdhb5
Insulin receptor binding
Translation initiation factor activity
Translation repressor activity
Estradiol 17-beta-dehydrogenase activity
H2-Eb1, Tec, Igbp1b, Nudcd1
Nudcd1, Mkrn3, Msmb, Cav2, Gnmt, Cryl1
Transferase activity, transferring glycosyl groups
Cytoskeleton organization and biogenesis
Isoprenoid biosynthetic process
Phosphotransferase activity, alcohol group as acceptor
Ribosome biogenesis and assembly
PPAR signaling pathway/lipid metabolic process
Breaks down the protein building block (amino acid) methionine
C80171, 1810030O07Rik, Es2el, Pter
The expression of STAT5ca in individual tumors does not correlate with the expression of its specifically affected genes
Does the deviation in gene expression between the two sets of tumors result from an early transgenic STAT5 effect, or does continuously deregulated STAT5 expression in the tumors govern the gene-expression profile?
STAT5ca and STATΔ750 expression levels were analyzed by real-time PCR in the respective tumors. In contrast to measurable STAT5ca transcripts, STAT5Δ750 mRNA did not reach detection levels. There was no significant correlation between expression of STAT5ca in the tumors and that of any of its affected genes [see additional file 2, bottom row of upper section]. Interestingly, within the genes affected specifically by STAT5ca, a few groups were identified which were linked by significant correlations between the expressions of their individual members. The genes with the highest number of internal correlations are presented [see Additional file 2]. The largest such group included Mapk8 (top row), Wnt8b, Ptn, Pramel6, Klra5, Edg3, Ebf2 and Tyrp1. Significant gene-expression correlations were also detected among the STAT5Δ750-affected tumors, the largest group including Wnt7b, Scgb1a1, Itga3, Mrpplf4 and Cspg2. Analyses of these gene groups by DAVID  and IPA software could not assemble the individual members into a single metabolic pathway or identify a mutual effector.
Distinct distribution into clusters was observed when the specifically affected genes were allocated. With the list of highly correlated genes as a core, the STAT5ca-affected genes were assembled into a small region in the upper section of cluster 3 (Fig. 2), confirming a lower level of expression compared to both the STAT5Δ750-induced tumors and the wild-type mammary gland. In contrast, genes affected by STAT5Δ750 were distinctly located in clusters 2 and 1. Their location in cluster 2 implied a higher level of expression compared to those specifically affected by the STAT5ca transgene or the wild-type mammary gland. Their presence in cluster 1 as well reflected the more even distribution between upregulated and downregulated genes.
Reproducing transgenic STAT5 effects in a "test set" of tumors
Validation of the differences in gene expression determined by the array analysis using real-time PCR of gene expression in a test set of tumors.
Tumor set 1
(Microarray analysis, arbitrary units)
Tumor set 2 (test set)
(Real-Time PCR analysis, 2-Δct × 1,000)
Fold change (BLG/STAT5Δ750/BLG/STAT5ca)
Fold change (BLG/STAT5Δ750/BLG/STAT5ca)
1.59 ± 0.21
0.43 ± 0.15
0.10 ± 0.02
0.04 ± 0.01
1.30 ± 0.19
0.51 ± 0.12
0.64 ± 0.189
0.43 ± .05
1.41 ± 0.31
0.44 ± 0.14
1.41 ± 0.380
0.37 ± 0.21
1.69 ± 0.48
0.46 ± 0.25
1.01 ± 0.219
0.48 ± 0.07
3.49 ± 0.97
0.86 ± 0.19
0.35 ± 0.118
0.25 ± 0.02
1.37 ± 0.21
0.41 ± 0.17
2.22 ± 1.159
0.47 ± 0.18
The constitutively active Stat5 and its C-terminally truncated variant have been implicated in cancer in laboratory animals and humans. Here we established a model system, based on mice expressing the individual transgenic STAT5 variants, to study the distinct variants' oncogenic role in the mammary gland.
Stat5 mediates proliferation in the mammary gland of pregnant females , and different effects of its two variants have been demonstrated in transgenic mice during the reproductive cycle: STAT5ca-induced proliferation during pregnancy and delayed apoptosis and tissue remodeling during involution . The opposite has been shown for the STAT5Δ750 variant [11, 12]. This diversity results from the structural differences between the STAT5 variants which involves either forced activation or lack of TAD-mediated interactions with the proteins, such as CBP/p300, that recruit acetylases . Both STAT5 variants seem to have a persistent effect, probably on chromatin structure and accessibility [16, 48, 49], which results in different profiles of gene expression in the developing tumors.
Differential expression of 364 genes was found between mammary tumors developed in BLG/STAT5ca- and BLG/STAT5Δ750-transgenic mice. These genes were involved in a set of cellular activities, many of which could still be associated with the two main processes mediated by Stat5 in the intact tissue: cell proliferation and cell death. The individual genes that established the specific mark of the STAT5 variants in the tumors were identified by gene-array analysis and their differential expression in the population was confirmed in a distinct "test set" of tumors using a different detection method, real-time PCR. Of special interest was the gene encoding the proto-oncogene Met, which was expressed in the STAT5Δ750-induced tumors at a 22-fold higher level than in their STAT5ca-induced counterparts. The product of Met is the hepatocyte growth factor receptor which encodes tyrosine-kinase activity. The ligand-activated cytoplasmic domain of the c-Met receptor induces growth motility, morphogenesis and angiogenesis . In breast cancer, c-Met overexpression is associated with tumor progression (reviewed in ) and has an independent predictive value for poor survival, even in early-stage patients with negative lymph nodes . Expression levels of Met were upregulated relative to intact tissue in the STAT5Δ750-induced tumors and downregulated in tumors developed in mice carrying the BLG/STAT5ca transgene. This differential expression suggests a more aggressive downstream cascade in the former. Other genes with exceptionally high expression levels in the STAT5Δ750-induced tumors, which were not affected by the STAT5ca variant, were proliferin and Igf2. Proliferin is a member of the prolactin family which is involved in progenitor cell expansion along the luminal and myoepithelial lineage  and Igf2 plays a pivotal role in fetal and cancer development by signaling via the IGF-I and insulin receptors, and activating the estrogen-signaling cascade . Unfortunately, studies on these genes [55–57] do not include or base additional lists of genes with altered expression profiles that might be compared to our data and aid in delineating the pathway(s) involved in the tumorigenic effect of STAT5Δ750.
A group of genes regulating carbohydrate metabolism and transport could also be discerned. These genes were not associated with Stat5 effects in the intact mammary gland but were differentially expressed in the two sets of tumors. Their role in the tumors could be related to altered levels of metabolism (i.e. higher metabolic rate in the STAT5Δ750-induced tumors) rather than cancer growth per se.
Expression of an additional set of 14 genes linked the effects of the BLG/STAT5Δ750 and the BLG/STAT5ca transgenes to the resulting phenotypes of poorly differentiated carcinoma or highly differentiated papillary adenocarcinoma, respectively (r = 0.97) [see Additional file 1, data sheet G and ]. These genes cover a wide variety of cellular functions: calcium sensitivity of the myofibrils (troponin I and tropomyosin 1), interaction between the cell and the extracellular matrix (endomucin), normal adipose tissue development (lipin 3), lipid metabolism, cellular growth and apoptosis (caveolin 2), mRNA metabolism (5'-3' exoribonuclease 1), mitochondrial fatty acid oxidation (acyl coenzyme A dehydrogenase), as well as kinase and protease activities (TAU tubulin kinase 1 and corin, respectively). Taken together, their diverse expression is most likely involved in determining tumor phenotype, i.e. associating higher proportions of the poorly differentiated carcinomas or the highly differentiated papillary adenocarcinomas with the expression of STAT5Δ750 or STAT5ca, respectively .
Overall, 94% of the genes specifically affected by STAT5ca were downregulated relative to their expression in the host tissue. This contrasts with the more equal specific effect of the STAT5Δ750 variant and provides additional evidence for the distinct routes via which the two STAT5 variants impose their mark on tumor growth and maintenance. The high number of tumor and growth suppressors defined among the STAT5ca-downregulated genes vs. the potent oncogenes (Met, Igf2) that were induced by the STAT5Δ750 variant may serve to further distinguish the routes via which they initiate and maintain tumorigenesis.
Substantial downregulation of gene expression has been demonstrated in breast cancers with bone marrow (BM) micrometastasis , and during the molecular transition from organ-confined to metastatic prostate cancer . Apparently, transcription repression is important for the metastatic process in these tissues. When compared with the list of genes downregulated by STAT5ca-, no overlap was observed for those associated with BM micrometastasis, and only four genes (Ptn, Cflar, Cnn1 and Mpdz) shared the list of downregulated genes mediating prostate metastasis. Thus, the resultant tumor characteristics are probably determined by a combination of the phenomenon of gene downregulation per se and the properties of the specifically affected genes (in this case by STAT5ca).
Our attempt to characterize the different roles of the STAT5 variants in mammary cancer development did not produce any evidence for the effect of transgenic STAT5 expression in the tumors. Neither distinct metabolic pathways nor a central mediator(s) were located. This suggests domination of an earlier STAT5 effect on the resulting gene-expression profiles. Several groups of genes, each composed of three to eight members, with significant internal correlations were identified among the sets specifically affected by STAT5ca and STAT5Δ750. This correlation was not associated with physical linkage. Within a single group, genes with correlated expression might control a distinct range of cellular functions such as angiogenesis (Ptn and Tyrp1 [60, 61]), apoptosis (Cflar, ) and morphogenesis (Ebf2, ). Their comparable levels of expression, defined by their location in a confined region within the cluster, suggest that although these genes are mapped to several loci, they may colocalize to a shared transcription site. The concept of several loci being targeted to a shared transcription site where they generate "transcription factories" with similar levels of expression was proposed a few years ago and is reviewed in [64, 65]. The relevance of this biological system to the specific STAT5 effects presented in this study remains to be determined. Regardless of the detailed mechanism involved in the effect of STAT5Δ750 and STAT5ca on tumor growth and maintenance, this study establishes the feasibility of identifying and distinguishing mammary tumors according the variants' signature on specific gene-expression profiles. To the best of our knowledge, this signature is specific to the effect of the STAT5 variants. However, its uniqueness will only be confirmed by the contribution of further gene-expression profiles that are specific for the effects of other oncogenes.
The C-terminally truncated form of Stat5 and its constitutively active variant retain oncogenic potency that is conveyed via distinct pathways and probably initiated during early stages of tumor development. The detailed results of this study may have clinical implications regarding the decision of whether, and when, to use putative anti-Stat5 therapy .
signal transducer and activator of transcription 5
constitutively active STAT5
This study was supported by a grant from the Israel Science Foundation, Israel Academy of Sciences, contract number: 706/04 to IB.
Contribution No. 536/08 from the ARO, The Volcani Center, Bet-Dagan, Israel.
- Moriggl R, Gouilleux-Gruart V, Jahne R, Berchtold S, Gartmann C, Liu X, Hennighausen L, Sotiropoulos A, Groner B, Gouilleux F: Deletion of the carboxyl-terminal transactivation domain of MGF-Stat5 results in sustained DNA binding and a dominant negative phenotype. Mol Cell Biol. 1996, 16: 5691-5700.PubMed CentralView ArticlePubMed
- Bowman T, Garcia R, Turkson J, Jove R: STATs in oncogenesis. Oncogene. 2000, 19: 2474-2488. 10.1038/sj.onc.1203527.View ArticlePubMed
- Hendry L, John S: Regulation of STAT signalling by proteolytic processing. Eur J Biochem. 2004, 271: 4613-4620. 10.1111/j.1432-1033.2004.04424.x.View ArticlePubMed
- Azam M, Erdjument-Bromage H, Kreider BL, Xia M, Quelle F, Basu R, Saris C, Tempst P, Ihle JN, Schindler C: Interleukin-3 signals through multiple isoforms of Stat5. Embo J. 1995, 14: 1402-1411.PubMed CentralPubMed
- Meyer J, Jucker M, Ostertag W, Stocking C: Carboxyl-truncated STAT5beta is generated by a nucleus-associated serine protease in early hematopoietic progenitors. Blood. 1998, 91: 1901-1908.PubMed
- Mitchell TJ, Whittaker SJ, John S: Dysregulated expression of COOH-terminally truncated Stat5 and loss of IL2-inducible Stat5-dependent gene expression in Sezary Syndrome. Cancer Res. 2003, 63: 9048-9054.PubMed
- Xia Z, Sait SN, Baer MR, Barcos M, Donohue KA, Lawrence D, Ford LA, Block AM, Baumann H, Wetzler M: Truncated STAT proteins are prevalent at relapse of acute myeloid leukemia. Leuk Res. 2001, 25: 473-482. 10.1016/S0145-2126(00)00158-2.View ArticlePubMed
- Ilaria RL, Van Etten RA: P210 and P190(BCR/ABL) induce the tyrosine phosphorylation and DNA binding activity of multiple specific STAT family members. J Biol Chem. 1996, 271: 31704-31710. 10.1074/jbc.271.49.31704.View ArticlePubMed
- Van Etten RA: Aberrant cytokine signaling in leukemia. Oncogene. 2007, 26: 6738-6749. 10.1038/sj.onc.1210758.PubMed CentralView ArticlePubMed
- Wagner KU, Rui H: Jak2/Stat5 signaling in mammogenesis, breast cancer initiation and progression. J Mammary Gland Biol Neoplasia. 2008, 13: 93-103. 10.1007/s10911-008-9062-z.View ArticlePubMed
- Iavnilovitch E, Groner B, Barash I: Overexpression and forced activation of stat5 in mammary gland of transgenic mice promotes cellular proliferation, enhances differentiation, and delays postlactational apoptosis. Mol Cancer Res. 2002, 1: 32-47.PubMed
- Iavnilovitch E, Eilon T, Groner B, Barash I: Expression of a carboxy terminally truncated Stat5 with no transactivation domain in the mammary glands of transgenic mice inhibits cell proliferation during pregnancy, delays onset of milk secretion, and induces apoptosis upon involution. Mol Reprod Dev. 2006, 73: 841-849. 10.1002/mrd.20479.View ArticlePubMed
- Cotarla I, Ren S, Zhang Y, Gehan E, Singh B, Furth PA: Stat5a is tyrosine phosphorylated and nuclear localized in a high proportion of human breast cancers. Int J Cancer. 2004, 108: 665-671. 10.1002/ijc.11619.View ArticlePubMed
- Humphreys RC, Hennighausen L: Signal transducer and activator of transcription 5a influences mammary epithelial cell survival and tumorigenesis. Cell Growth Differ. 1999, 10: 685-694.PubMed
- Iavnilovitch E, Cardiff RD, Groner B, Barash I: Deregulation of Stat5 expression and activation causes mammary tumors in transgenic mice. Int J Cancer. 2004, 112: 607-619. 10.1002/ijc.20484.View ArticlePubMed
- Eilon T, Groner B, Barash I: Tumors caused by overexpression and forced activation of Stat5 in mammary epithelial cells of transgenic mice are parity-dependent and developed in aged, postestropausal females. Int J Cancer. 2007, 121: 1892-1902. 10.1002/ijc.22954.View ArticlePubMed
- Perou CM, Sorlie T, Eisen MB, Rijn van de M, Jeffrey SS, Rees CA, Pollack JR, Ross DT, Johnsen H, Akslen LA, Fluge O, Pergamenschikov A, Williams C, Zhu SX, Lønning PE, Børresen-Dale AL, Brown PO, Botstein D: Molecular portraits of human breast tumors. Nature. 2000, 406: 747-752. 10.1038/35021093.View ArticlePubMed
- van 't Veer LJ, Dai H, Vijver van de MJ, He YD, Hart AA, Mao M, Peterse HL, Kooy van der K, Marton MJ, Witteveen AT, Schreiber GJ, Kerkhoven RM, Roberts C, Linsley PS, Bernards R, Friend SH: Gene expression profiling predicts clinical outcome of breast cancer. Nature. 2002, 415: 530-536. 10.1038/415530a.View ArticlePubMed
- Ein-Dor L, Kela I, Getz G, Givol D, Domany E: Outcome signature genes in breast cancer: is there a unique set. Bioinformatics. 2005, 21: 171-178. 10.1093/bioinformatics/bth469.View ArticlePubMed
- Dressman HK, Hans C, Bild A, Olson JA, Rosen E, Marcom PK, Liotcheva VB, Jones EL, Vujaskovic Z, Marks J, Dewhirst MW, West M, Nevins JR, Blackwell K: Gene expression profiles of multiple breast cancer phenotypes and response to neoadjuvant chemotherapy. Clin Cancer Res. 2006, 12: 819-826. 10.1158/1078-0432.CCR-05-1447.View ArticlePubMed
- Ye Y, Qiu TH, Kavanaugh C, Green JE: Molecular mechanisms of breast cancer progression: lessons from mouse mammary cancer models and gene expression profiling. Breast Dis. 2004, 19: 69-82.PubMed
- Qiu TH, Chandramouli GV, Hunter KW, Alkharouf NW, Green JE, Liu ET: Global expression profiling identifies signatures of tumor virulence in MMTV-PyMT-transgenic mice: correlation to human disease. Cancer Res. 2004, 64: 5973-5981. 10.1158/0008-5472.CAN-04-0242.View ArticlePubMed
- Eilon T, Barash I: Different gene-expression profiles for the poorly differentiated carcinoma and the highly differentiated papillary adenocarcinoma in mammary glands support distinct metabolic pathways. BMC Cancer. 2008, 8: 270-10.1186/1471-2407-8-270.PubMed CentralView ArticlePubMed
- Desai KV, Xiao N, Wang W, Gangi L, Greene J, Powell JI, Dickson R, Furth P, Hunter K, Kucherlapati R, Simon R, Liu ET, Green JE: Initiating oncogenic event determines gene-expression patterns of human breast cancer models. Proc Natl Acad Sci USA. 2002, 99: 6967-6972. 10.1073/pnas.102172399.PubMed CentralView ArticlePubMed
- Lauren PD: Algorithm to model gene expression on Affymetrix chips without the use of MM cells. IEEE Trans Nanobioscience. 2003, 2: 163-170. 10.1109/TNB.2003.817020.View ArticlePubMed
- Miller LD, Long PM, Wong L, Mukherjee S, McShane LM, Liu ET: Optimal gene expression analysis by microarrays. Cancer Cell. 2002, 2: 353-361. 10.1016/S1535-6108(02)00181-2.View ArticlePubMed
- Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci USA. 1998, 95: 14863-14868. 10.1073/pnas.95.25.14863.PubMed CentralView ArticlePubMed
- Yoshida Y, Togi K, Matsumae H, Nakashima Y, Kojima Y, Yamamoto H, Ono K, Nakamura T, Kita T, Tanaka M: CCN1 protects cardiac myocytes from oxidative stress via beta1 integrin-Akt pathway. Biochem Biophys Res Commun. 2007, 355: 611-618. 10.1016/j.bbrc.2007.01.195.View ArticlePubMed
- Gressner OA, Lahme B, Siluschek M, Rehbein K, Herrmann J, Weiskirchen R, Gressner AM: Activation of TGF-beta within cultured hepatocytes and in liver injury leads to intracrine signaling with expression of connective tissue growth factor. J Cell Mol Med. 2008, 12: 2717-2730. 10.1111/j.1582-4934.2008.00260.x.PubMed CentralView ArticlePubMed
- Yu Y, Gao Y, Wang H, Huang L, Qin J, Guo R, Song M, Yu S, Chen J, Cui B, Gao P: The matrix protein CCN1 (CYR61) promotes proliferation, migration and tube formation of endothelial progenitor cells. Exp Cell Res. 2008, 314: 3198-3208. 10.1016/j.yexcr.2008.08.001.View ArticlePubMed
- O'Kelly J, Chung A, Lemp N, Chumakova K, Yin D, Wang HJ, Said J, Gui D, Miller CW, Karlan BY, Koeffler HP: Functional domains of CCN1 (Cyr61) regulate breast cancer progression. Int J Oncol. 2008, 33: 59-67.PubMed
- Egea M, Meton I, Cordoba M, Fernandez F, Baanante IV: Role of Sp1 and SREBP-1a in the insulin-mediated regulation of glucokinase transcription in the liver of gilthead sea bream (Sparus aurata). Gen Comp Endocrinol. 2008, 155: 359-367. 10.1016/j.ygcen.2007.06.018.View ArticlePubMed
- Fukumoto S, Miner JH, Ida H, Fukumoto E, Yuasa K, Miyazaki H, Hoffman MP, Yamada Y: Laminin alpha5 is required for dental epithelium growth and polarity and the development of tooth bud and shape. J Biol Chem. 2006, 281: 5008-5016. 10.1074/jbc.M509295200.View ArticlePubMed
- Ikeda W, Nakanishi H, Takekuni K, Itoh S, Takai Y: Identification of splicing variants of Frabin with partly different functions and tissue distribution. Biochem Biophys Res Commun. 2001, 286: 1066-1072. 10.1006/bbrc.2001.5481.View ArticlePubMed
- Ruckhaberle E, Rody A, Engels K, Gaetje R, von Minckwitz G, Schiffmann S, Grosch S, Geisslinger G, Holtrich U, Karn T, Kaufmann M: Microarray analysis of altered sphingolipid metabolism reveals prognostic significance of sphingosine kinase 1 in breast cancer. Breast Cancer Res Treat. 2008, 112: 41-52. 10.1007/s10549-007-9836-9.View ArticlePubMed
- Catalan V, Gomez-Ambrosi J, Pastor C, Rotellar F, Silva C, Rodriguez A, Gil MJ, Cienfuegos JA, Salvador J, Vendrell J, Frühbeck G: Influence of Morbid Obesity and Insulin Resistance on Gene Expression Levels of AQP7 in Visceral Adipose Tissue and AQP9 in Liver. Obes Surg. 2008, 18: 695-701. 10.1007/s11695-008-9453-7.View ArticlePubMed
- Shin JH, Choi GS, Kang WH, Myung KB: Sphingosine 1-phosphate triggers apoptotic signal for B16 melanoma cells via ERK and caspase activation. J Korean Med Sci. 2007, 22: 298-304. 10.3346/jkms.2007.22.2.298.PubMed CentralView ArticlePubMed
- Kim IJ, Lee KW, Park BY, Lee JK, Park J, Choi IY, Eom SJ, Chang TS, Kim MJ, Yeom YI, Chang SK, Lee YD, Choi EJ, Han PL: Molecular cloning of multiple splicing variants of JIP-1 preferentially expressed in brain. J Neurochem. 1999, 72: 1335-1343. 10.1046/j.1471-4159.1999.721335.x.View ArticlePubMed
- Uehara N, Unami A, Kiyozuka Y, Shikata N, Oishi Y, Tsubura A: Parous mammary glands exhibit distinct alterations in gene expression and proliferation responsiveness to carcinogenic stimuli in Lewis rats. Oncol Rep. 2006, 15: 903-911.PubMed
- Emslie D, D'Costa K, Hasbold J, Metcalf D, Takatsu K, Hodgkin PO, Corcoran LM: Oct2 enhances antibody-secreting cell differentiation through regulation of IL-5 receptor alpha chain expression on activated B cells. J Exp Med. 2008, 205: 409-421. 10.1084/jem.20072049.PubMed CentralView ArticlePubMed
- Jimenez MA, Akerblad P, Sigvardsson M, Rosen ED: Critical role for Ebf1 and Ebf2 in the adipogenic transcriptional cascade. Mol Cell Biol. 2007, 27: 743-757. 10.1128/MCB.01557-06.PubMed CentralView ArticlePubMed
- Mahajan KN, Mitchell BS: Role of human Pso4 in mammalian DNA repair and association with terminal deoxynucleotidyl transferase. Proc Natl Acad Sci USA. 2003, 100: 10746-10751. 10.1073/pnas.1631060100.PubMed CentralView ArticlePubMed
- Hong J, Doebele RC, Lingen MW, Quilliam LA, Tang WJ, Rosner MR: Anthrax edema toxin inhibits endothelial cell chemotaxis via Epac and Rap1. J Biol Chem. 2007, 282: 19781-19787. 10.1074/jbc.M700128200.View ArticlePubMed
- Wang-Rodriguez J, Urquidi V, Rivard A, Goodison S: Elevated osteopontin and thrombospondin expression identifies malignant human breast carcinoma but is not indicative of metastatic status. Breast Cancer Res. 2003, 5: R136-143. 10.1186/bcr620.PubMed CentralView ArticlePubMed
- Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4: P3-10.1186/gb-2003-4-5-p3.View ArticlePubMed
- Cui Y, Riedlinger G, Miyoshi K, Tang W, Li C, Deng CX, Robinson GW, Hennighausen L: Inactivation of Stat5 in mouse mammary epithelium during pregnancy reveals distinct functions in cell proliferation, survival, and differentiation. Mol Cell Biol. 2004, 24: 8037-8047. 10.1128/MCB.24.18.8037-8047.2004.PubMed CentralView ArticlePubMed
- Pfitzner E, Jahne R, Wissler M, Stoecklin E, Groner B: p300/CREB-binding protein enhances the prolactin-mediated transcriptional induction through direct interaction with the transactivation domain of Stat5, but does not participate in the Stat5-mediated suppression of the glucocorticoid response. Mol Endocrinol. 1998, 12: 1582-1593. 10.1210/me.12.10.1582.View ArticlePubMed
- Rascle A, Johnston JA, Amati B: Deacetylase activity is required for recruitment of the basal transcription machinery and transactivation by STAT5. Mol Cell Biol. 2003, 23: 4162-4173. 10.1128/MCB.23.12.4162-4173.2003.PubMed CentralView ArticlePubMed
- Rascle A, Lees E: Chromatin acetylation and remodeling at the Cis promoter during STAT5-induced transcription. Nucleic Acids Res. 2003, 31: 6882-6890. 10.1093/nar/gkg907.PubMed CentralView ArticlePubMed
- Weidner KM, Sachs M, Birchmeier W: The Met receptor tyrosine kinase transduces motility, proliferation, and morphogenic signals of scatter factor/hepatocyte growth factor in epithelial cells. J Cell Biol. 1993, 121: 145-154. 10.1083/jcb.121.1.145.View ArticlePubMed
- Lindemann K, Resau J, Nahrig J, Kort E, Leeser B, Annecke K, Welk A, Schafer J, Woude Vande GF, Lengyel E, Harbeck N: Differential expression of c-Met, its ligand HGF/SF and HER2/neu in DCIS and adjacent normal breast tissue. Histopathology. 2007, 51: 54-62. 10.1111/j.1365-2559.2007.02732.x.View ArticlePubMed
- Ghoussoub RA, Dillon DA, D'Aquila T, Rimm EB, Fearon ER, Rimm DL: Expression of c-met is a strong independent prognostic factor in breast carcinoma. Cancer. 1998, 82: 1513-1520. 10.1002/(SICI)1097-0142(19980415)82:8<1513::AID-CNCR13>3.0.CO;2-7.View ArticlePubMed
- Wang XY, Yin Y, Yuan H, Sakamaki T, Okano H, Glazer RI: Musashi1 modulates mammary progenitor cell expansion through proliferin-mediated activation of the Wnt and Notch pathways. Mol Cell Biol. 2008, 28: 3589-3599. 10.1128/MCB.00040-08.PubMed CentralView ArticlePubMed
- Singh SK, Moretta D, Almaguel F, Wall NR, De Leon M, De Leon D: Differential effect of proIGF-II and IGF-II on resveratrol induced cell death by regulating survivin cellular localization and mitochondrial depolarization in breast cancer cells. Growth Factors. 2007, 25: 363-372. 10.1080/08977190801886905.PubMed CentralView ArticlePubMed
- Cardiff RD, Wellings SR: The comparative pathology of human and mouse mammary glands. J Mammary Gland Biol Neoplasia. 1999, 4: 105-122. 10.1023/A:1018712905244.View ArticlePubMed
- Gallego MI, Bierie B, Hennighausen L: Targeted expression of HGF/SF in mouse mammary epithelium leads to metastatic adenosquamous carcinomas through the activation of multiple signal transduction pathways. Oncogene. 2003, 22: 8498-8508. 10.1038/sj.onc.1207063.View ArticlePubMed
- Welm AL, Kim S, Welm BE, Bishop JM: MET and MYC cooperate in mammary tumorigenesis. Proc Natl Acad Sci USA. 2005, 102: 4324-4329. 10.1073/pnas.0500470102.PubMed CentralView ArticlePubMed
- Woelfle U, Cloos J, Sauter G, Riethdorf L, Janicke F, van Diest P, Brakenhoff R, Pantel K: Molecular signature associated with bone marrow micrometastasis in human breast cancer. Cancer Res. 2003, 63: 5679-5684.PubMed
- Varambally S, Dhanasekaran SM, Zhou M, Barrette TR, Kumar-Sinha C, Sanda MG, Ghosh D, Pienta KJ, Sewalt RG, Otte AP, Rubin MA, Chinnaiyan AM: The polycomb group protein EZH2 is involved in progression of prostate cancer. Nature. 2002, 419: 624-629. 10.1038/nature01075.View ArticlePubMed
- Chang Y, Zuka M, Perez-Pinera P, Astudillo A, Mortimer J, Berenson JR, Deuel TF: Secretion of pleiotrophin stimulates breast cancer progression through remodeling of the tumor microenvironment. Proc Natl Acad Sci USA. 2007, 104: 10888-10893. 10.1073/pnas.0704366104.PubMed CentralView ArticlePubMed
- Rousseau B, Larrieu-Lahargue F, Javerzat S, Guilhem-Ducleon F, Beermann F, Bikfalvi A: The tyrp1-Tag/tyrp1-FGFR1-DN bigenic mouse: a model for selective inhibition of tumor development, angiogenesis, and invasion into the neural tissue by blockade of fibroblast growth factor receptor activity. Cancer Res. 2004, 64: 2490-2495. 10.1158/0008-5472.CAN-03-3623.View ArticlePubMed
- Yeh WC, Itie A, Elia AJ, Ng M, Shu HB, Wakeham A, Mirtsos C, Suzuki N, Bonnard M, Goeddel DV, Mak TW: Requirement for Casper (c-FLIP) in regulation of death receptor-induced apoptosis and embryonic development. Immunity. 2000, 12: 633-642. 10.1016/S1074-7613(00)80214-9.View ArticlePubMed
- Zhao LY, Niu Y, Santiago A, Liu J, Albert SH, Robertson KD, Liao D: An EBF3-mediated transcriptional program that induces cell cycle arrest and apoptosis. Cancer Res. 2006, 66: 9445-9452. 10.1158/0008-5472.CAN-06-1713.View ArticlePubMed
- Marenduzzo D, Faro-Trindade I, Cook PR: What are the molecular ties that maintain genomic loops?. Trends Genet. 2007, 23: 126-133. 10.1016/j.tig.2007.01.007.View ArticlePubMed
- Trinkle-Mulcahy L, Lamond AI: Nuclear functions in space and time: gene expression in a dynamic, constrained environment. FEBS Lett. 2008, 582: 1960-1970. 10.1016/j.febslet.2008.04.029.View ArticlePubMed
- Desrivieres S, Kunz C, Barash I, Vafaizadeh V, Borghouts C, Groner B: The biological functions of the versatile transcription factors STAT3 and STAT5 and new strategies for their targeted inhibition. J Mammary Gland Biol Neoplasia. 2006, 11: 75-87. 10.1007/s10911-006-9014-4.View ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.