Different gene-expression profiles for the poorly differentiated carcinoma and the highly differentiated papillary adenocarcinoma in mammary glands support distinct metabolic pathways

A Bioinformatic Pipeline Places STAT5A as a miR-650 Target in Poorly Differentiated Aggressive Breast Cancer

Breast cancer (BRCA) is a leading cause of mortality among women. Tumors often acquire aggressive features through genomic aberrations affecting cellular programs, e.g., the epithelial to mesenchymal transition (EMT). EMT facilitates metastasis leading to poor prognosis. We previously observed a correlation between an amplification of miR-650 (Amp-650) and EMT features in BRCA samples isolated from Mexican patients. In this study, we explored the cBioportal database aiming to extend that observation and better understand the importance of Amp-650 for BRCA aggressiveness. We found that Amp-650 is more frequent in aggressive molecular subtypes of BRCA, as well as in high grade poorly differentiated tumors, which we confirmed in an external miRNA expression database. We performed differential expression analysis on samples harboring Amp-650, taking advantage of gene target prediction tools and tumor suppressor gene databases to mine several hundreds of differentially underexpressed genes. We observed STAT5A as a likely putative target gene for miR-650 in aggressive poorly differentiated BRCA. Samples with both Amp-650 and low expression of STAT5A had less overall survival than samples with either or none of the alterations. No target gene has been described for miR-650 in BRCA, thus, this bioinformatic study provides valuable information that should be corroborated experimentally.

H2AX Promoter Demethylation at Specific Sites Plays a Role in STAT5-Induced Tumorigenesis

Deregulated STAT5 activity in the mammary gland of transgenic mice results in parity-dependent latent tumorigenesis. The trigger for cell transformation was previously associated with hyperactivation of the H2AX proximal promoter in a small basal cell population during pregnancy. The current study focuses on the latent activation of tumor development. H2AX was highly expressed in carcinoma and adenocarcinoma as compared to the multiparous mammary gland, whereas pSTAT5 expression decreased in a tumor type-dependent manner. In contrast to the pregnant gland, no positive correlation between H2AX and pSTAT5 expression could be defined in carcinoma and adenocarcinoma. Using targeted methylation analysis, the methylation profile of the H2AX promoter was characterized in the intact gland and tumors. Average H2AX promoter methylation in the tumors was relatively high (~90%), but did not exceed that of the multiparous gland; 5mC methylation was higher in the differentiated tumors and negatively correlated with its oxidative product 5hmC and H2AX expression. Individual analysis of 25 H2AX promoter-methylation sites revealed two consecutive CpGs at positions -77 and - 54 that were actively demethylated in the multiparous gland, but not in their age-matched virgin counterpart. The different methylation profiles at these sites distinguished tumor types and may assume a prognostic role. In-silico and ChIP analyses revealed overlapping methylation-independent SP1-binding and methylation-dependent p53-binding to these sites. We propose that interference with SP1-assisted p53-binding to these sites abrogates H2AX's ability to arrest the cell cycle upon DNA damage, and contributes to triggering latent development of STAT5-induced tumors in estrapausal multiparous mice.

Luminal STAT5 mediates H2AX promoter activity in distinct population of basal mammary epithelial cells

Deregulated STAT5 activity in the mammary gland causes parity-dependent tumorigenesis. Epithelial cell cultures transfected with constitutively active STAT5 express higher levels of the histone H2AX than their non-transfected counterparts. Higher H2AX expression may be involved in tumorigenesis. Here, we aimed to link high STAT5 activity to H2AX–GFP expression by looking for distinct types of mammary cells that express these proteins. In vitro and in transgenic mice, only 0.2 and 0.02%, respectively, of the cells expressed the H2AX–GFP hybrid gene. Its expression correlated with that of the endogenous H2AX gene, suggesting that detectable H2AX–GFP expression marks high levels of H2AX transcript. Methylation of the H2AX promoter characterized non-GFP-expressing H2AX–GFP cells and was inversely correlated with promoter activity. Administration of 5-azacytidine increased H2AX promoter activity in an activated STAT5-dependent manner. In transgenic mice, H2AX–GFP expression peaked at pregnancy. The number of H2AX–GFP-expressing cells and GFP expression decreased in a Stat5a-null background and increased in mice expressing the hyperactivated STAT5. Importantly, H2AX–GFP activity was allocated to basal mammary cells lacking stem-cell properties, whereas STAT5 hyperactivity was detected in the adjacent luminal cells. Knockdown of RANKL by siRNA suggested its involvement in signaling between the two layers. These results suggest paracrine activation of H2AX via promoter demethylation in specific populations of basal mammary cells that is induced by a signal from neighboring luminal cells with hyper STAT5 activity. This pathway provides an alternative route for the luminally confined STAT5 to affect basal mammary cell activity.

Signal transducer and activator of transcription 5 as a key signaling pathway in normal mammary gland developmental biology and breast cancer

STAT5 consists of two proteins, STAT5A/B, that impact mammary cell differentiation, proliferation, and survival. In normal development, STAT5 expression and activity are regulated by prolactin signaling with JAK2/ELF5, EGF signaling networks that include c-Src, and growth hormone, insulin growth factor, estrogen, and progesterone signaling pathways. In cancer, erythropoietin signaling can also regulate STAT5. Activation levels are influenced by AKT, caveolin, PIKE-A, Pak1, c-Myb, Brk, beta-integrin, dystroglycan, other STATs, and STAT pathway molecules JAK1, Shp2, and SOCS. TGF-β and PTPN9 can downregulate prolactin- and EGF-mediated STAT5 activation, respectively. IGF, AKT, RANKL, cyclin D1, BCL6, and HSP90A lie downstream of STAT5.

Conditional repression of STAT5 expression during lactation reveals its exclusive roles in mammary gland morphology, milk-protein gene expression, and neonate growth

The role of Stat5 in maintaining adequate lactation was studied in Stat5a(-/-) mice expressing a conditionally suppressed transgenic STAT5 in their mammary glands. This system enables distinguishing STAT5's effects on lactation from its contribution to mammary development during gestation. Females were allowed to express STAT5 during their first pregnancy. After delivery, STAT5 levels were manipulated by doxycycline administration and withdrawal. In two lines of genetically modified mice, the absence of STAT5 expression during the first 10 days of lactation resulted in a decrease of 29% or 41% in newborn weight gain. The STAT5-dependent decrease in growth was recoverable, but not completely reversible, particularly when STAT5 expression was omitted for the first 4 days of lactation. Within the first 10 days of STAT5-omitted lactation, alveolar occupancy regressed by 50% compared to that measured at delivery. By Day 10, only 18% of the fat-pad area was involved in milk production. The alveolar regression caused by 4 days of STAT5 deficiency was reversible, but neonate growth remained delayed. STAT5 deficiency resulted in reduced estrogen receptor α and connexin 32 gene expression, accompanied by delayed induction of both anti- and pro-apoptotic Bcl-2 family members. An increase in Gata-3 expression may reflect an attempt to maintain alveolar progenitors. A decrease of 39% and 23% in WAP and α-lactalbumin expression, respectively, with no associated effects on β-casein, also resulted from lack of STAT5 expression in the first 10 days of lactation. This deficiency enhances the major effect of alveolar regression on delayed weight gain in newborns.

Distinct gene-expression profiles characterize mammary tumors developed in transgenic mice expressing constitutively active and C-terminally truncated variants of STAT5

Background

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.

Results

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.

Conclusion

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.

Approaching solid tumor heterogeneity on a cellular basis by tissue proteomics using laser capture microdissection and biological mass spectrometry

The purpose of this study was to examine solid tumor heterogeneity on a cellular basis using tissue proteomics that relies on a functional relationship between Laser Capture Microdissection (LCM) and biological mass spectrometry (MS). With the use of LCM, homogeneous regions of cells exhibiting uniform histology were isolated and captured from fresh frozen tissue specimens, which were obtained from a human lymph node containing breast carcinoma metastasis. Six specimens approximately 50,000 cell each (three from tumor proper and three from tumor stroma) were collected by LCM. Specimens were processed directly on LCM caps, using sonication in buffered methanol to lyse captured cells, solubilize, and digest extracted proteins. Prepared samples were analyzed by LC/MS/MS resulting in more than 500 unique protein identifications. Decoy database searching revealed a false-positive rate between 5 and 10%. Subcellular localization analysis for stromal cells revealed plasma membrane 14%, cytoplasm 39%, nucleus 11%, extracellular space 27%, and unknown 9%; and tumor cell results were 5%, 58%, 26%, 4%, and 7%, respectively. Western blot analysis confirmed specific linkage of validated proteins to underlying pathology and their potential role in solid tumor heterogeneity. With continued research and optimization of this method including analysis of additional clinical specimens, this approach may lead to an improved understanding of tumor heterogeneity, and serve as a platform for solid tumor biomarker discovery.