- Research article
- Open Access
Distinct gene-expression profiles characterize mammary tumors developed in transgenic mice expressing constitutively active and C-terminally truncated variants of STAT5
BMC Genomics volume 10, Article number: 231 (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
RNA was independently extracted from six tumors developed in mice carrying the BLG/STAT5Δ750 transgene and from seven tumors originating in mice carrying BLG/STAT5ca. Tumor phenotypes were either poorly differentiated carcinoma or papillary adenocarcinoma (Fig. 1 and described in ). Following microarry hybridization, we identified a set of 381 features (364 genes) exhibiting a significant (P < 0.05) twofold difference in expression level between tumors originated from mice carrying the two transgenic STAT5 variants [see Additional file 1, data sheet B]. Their grouping according to type of transgenic STAT5 variant was confirmed by principal component analysis (PCA, ) and unsupervised hierarchical clustering (Fig. 1).
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 .
Cyr61 can be also found among the differentially expressed genes regulating the formation of cellular protrusions and filopodia–cellular extensions needed for cellular interactions and movement . Most of the listed genes affecting this process, including the laminin α5 (Lama5)  and the VEGF receptor PVR , were expressed at higher levels in the STAT5Δ750-induced tumors. This suggests higher involvement of this variant in cell movement and tumor extension. The GDP/GTP exchange protein Fabrin 4 (FDG4,  and the ganglioside GD3-synthase (ST8SIA)1 which marks ER-negative breast cancer tumors  were the only genes in this context that were expressed at a higher level in the STAT5ca-induced tumors.
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
Unsupervised hierarchical clustering assembled the genes that were differentially expressed between the two sets of STAT5-affected tumors into four clusters (Fig. 2 and [see Additional file 1, data sheet D]). Cluster 1 contains genes that were more highly expressed in the STAT5ca-induced tumors than in their STAT5Δ750 counterparts, resulting from downregulation of the STAT5Δ750-induced genes compared to their expression in the mammary gland (represented by the value "1"). In contrast, genes allocated to clusters 2, 3 and 4 were more highly expressed in the STAT5Δ750-induced tumors. In cluster 2, the STAT5Δ750 affected genes were, in general, also more highly expressed than in the wild-type gland. In cluster 3, the STAT5Δ750-affected genes maintained mammary gland levels of expression while their STAT5ca counterparts were expressed at a lower level. Finally, in cluster 4, differential expression relative to the mammary gland was noted: higher expression of genes induced by the BLG/STAT5Δ750 transgene, lower expression of genes affected by its BLG/STAT5ca counterpart.
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
Specifically affected genes differed in their expression between an individual tumor set and its counterpart, and between that set and the intact mammary gland. To identify genes which were specifically affected by each of the STAT5 variants, three mammary gland samples were added to the analysis, and the level of expression of the 364 genes exhibiting differential expression among STAT5-induced tumors was re-evaluated (Fig. 3). This analysis identified 52 genes associated with the effect of the STAT5Δ750 transgene [(STAT5Δ750 ≠ STAT5ca = mammary gland (M.G.)], 50 genes that were specifically affected by STAT5ca (STAT5ca ≠ STAT5Δ750 = M.G.), and 62 genes with differential expression in the three groups (STAT5Δ750 ≠ STAT5ca ≠ M.G., and [see Additional file 1, data sheet F]). Unexpectedly, 94% of the STAT5ca-affected genes (47 genes) were downregulated compared to their expression in the STAT5Δ750-induced tumors and the wild-type mammary gland. A considerable number of these genes could be associated with anti-proliferative or tumor-suppressive processes, including apoptosis inducers Edg3 and Jip-2 [37, 38], differentiation inducers Msln, Il5ra and Ebf2 [39–41], Prp19 (also known as Pso4) which is involved in DNA repair , Ragef3 which blocks chemotaxis induced by angiogenic factors  and Tyrp, which is expressed in poorly metastatic breast cancer and whose downregulation augments metastasis .
A more equal deviation characterized the STAT5Δ750-affected genes. Only 38% of these genes were downregulated, whereas the rest were expressed at significantly higher levels compared to their expression in the STAT5ca-induced tumors or the wild-type mammary gland. In contrast to the STAT5ca-affected genes that were located in cluster 3, their STAT5Δ750 counterparts spread mainly among cluster 1 (20 genes, 38%) and cluster 2 (29 genes, 56%), with a residual presence in cluster 4 (3 genes, 6%). Upregulated and downregulated genes were annotated with different activities (Table 2). Upregulated genes affected angiogenesis, cell adhesion, cell-cell signaling, and progression through the cell cycle. Downregulated genes were involved in insulin receptor binding, DNA replication, chromatin remodeling and G-protein receptor pathways.
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
Our next step was to confirm the "transgenic signature" of the STAT5 variants in independent groups of tumors. From the list of 50 genes specifically affected by STAT5ca, six candidates were selected and their expression was analyzed by real-time-PCR in a separate set of tumors caused by the two STAT5 variants. This "test set" included seven tumors that had developed in BLG/STAT5ca-transgenic mice (four papillary adenocarcinomas, one squamous carcinoma and two poorly differentiated carcinomas) and six tumors that had developed in the BLG/STAT5Δ750 mice (one papillary adenocarcinoma, three micropapillary adenocarcinomas, one squamous carcinoma and one adenosquamous carcinoma). As demonstrated in Table 3, the differences between the expressions of genes affected by STAT5Δ750 vs. STAT5ca detected in the original set of tumors that was tested by the array analysis were comparable to those found in the "test set" which was subjected to real-time PCR.
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
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.
Bowman T, Garcia R, Turkson J, Jove R: STATs in oncogenesis. Oncogene. 2000, 19: 2474-2488. 10.1038/sj.onc.1203527.
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.
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.
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.
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.
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.
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.
Van Etten RA: Aberrant cytokine signaling in leukemia. Oncogene. 2007, 26: 6738-6749. 10.1038/sj.onc.1210758.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
The authors declare that they have no competing interests.
TE designed and performed the research and also helped in writing the manuscript. IB designed the research and wrote the manuscript.
Electronic supplementary material
Additional file 1: Supporting information. Primers list and bioinformatical analyses of differentially expressed genes in tumors caused by the STAT5ca and STAT5Δ750 variants. Data sheet A. List of primers used to amplify coding regions of the following listed genes. Data sheet B. List of features with a significant (P < 0.05), over twofold difference in expression between STAT5Δ750- and STAT5ca-induced tumors. Data sheet C. IPA classification of the differentially expressed genes in tumors induced by STAT5ca and the STAT5Δ750 transgenes. List of genes comprising the process annotation to biological function and/or disease is presented. Data sheet D. Clustering of genes with a significant (P < 0.05), over twofold difference in expression in tumors caused by the two STAT5 variants. Data sheet E. Clustering and IPA classification of the differentially expressed genes in tumors caused by STAT5ca and STAT5Δ750 variants. Data sheet F. List of specifically affected genes and those differentially expressed in each of the tumor sets relative to the mammary gland. Data sheet G. Set of 14 genes with expression levels that associate the effect of STAT5Δ750 with the resulting carcinoma phenotype and the effect of STAT5ca with the papillary adenocarcinoma phenotype. (XLS 180 KB)
Additional file 2: Correlations between the expression of STAT5ca and its specifically affected genes, and among genes specifically affected by the STAT5 variants. The table presents the correlation values between the expression levels of genes specifically affected by the STAT5ca and STAT5Δ750 variants in the tumors, and between the expression of STAT5ca and these genes. (DOC 159 KB)
About this article
Cite this article
Eilon, T., Barash, I. Distinct gene-expression profiles characterize mammary tumors developed in transgenic mice expressing constitutively active and C-terminally truncated variants of STAT5. BMC Genomics 10, 231 (2009). https://doi.org/10.1186/1471-2164-10-231
- Mammary Gland
- Mammary Tumor
- Papillary Adenocarcinoma
- Unsupervised Hierarchical Cluster
- Ingenuity Pathway Analysis Software