Distinct genes related to drug response identified in ER positive and ER negative breast cancer cell lines

Breast cancer patients have different responses to chemotherapeutic treatments. Genes associated with drug response can provide insight to understand the mechanisms of drug resistance, identify promising therapeutic opportunities, and facilitate personalized treatment. Estrogen receptor (ER) positive and ER negative breast cancer have distinct clinical behavior and molecular properties. However, to date, few studies have rigorously assessed drug response genes in them. In this study, our goal was to systematically identify genes associated with multidrug response in ER positive and ER negative breast cancer cell lines. We tested 27 human breast cell lines for response to seven chemotherapeutic agents (cyclophosphamide, docetaxel, doxorubicin, epirubicin, fluorouracil, gemcitabine, and paclitaxel). We integrated publicly available gene expression profiles of these cell lines with their in vitro drug response patterns, then applied meta-analysis to identify genes related to multidrug response in ER positive and ER negative cells separately. One hundred eighty-eight genes were identified as related to multidrug response in ER positive and 32 genes in ER negative breast cell lines. Of these, only three genes (DBI, TOP2A, and PMVK) were common to both cell types. TOP2A was positively associated with drug response, and DBI was negatively associated with drug response. Interestingly, PMVK was positively associated with drug response in ER positive cells and negatively in ER negative cells. Functional analysis showed that while cell cycle affects drug response in both ER positive and negative cells, most biological processes that are involved in drug response are distinct. A number of signaling pathways that are uniquely enriched in ER positive cells have complex cross talk with ER signaling, while in ER negative cells, enriched pathways are related to metabolic functions. Taken together, our analysis indicates that distinct mechanisms are involved in multidrug response in ER positive and ER negative breast cells.     

The Mutant p53-Targeting Compound APR-246 Induces ROS-Modulating Genes in Breast Cancer Cells

TP53 is the most frequently mutated gene in human cancer and thus an attractive target for novel cancer therapy. Several compounds that can reactive mutant p53 protein have been identified. APR-246 is currently being tested in a phase II clinical trial in high-grade serous ovarian cancer. We have used RNA-seq analysis to study the effects of APR-246 on gene expression in human breast cancer cell lines. Although the effect of APR-246 on gene expression was largely cell line dependent, six genes were upregulated across all three cell lines studied, i.e., TRIM16, SLC7A11, TXNRD1, SRXN1, LOC344887, and SLC7A11-AS1. We did not detect upregulation of canonical p53 target genes such as CDKN1A (p21), 14-3-3σ, BBC3 (PUMA), and PMAIP1 (NOXA) by RNA-seq, but these genes were induced according to analysis by qPCR. Gene ontology analysis showed that APR-246 induced changes in pathways such as response to oxidative stress, gene expression, cell proliferation, response to nitrosative stress, and the glutathione biosynthesis process. Our results are consistent with the dual action of APR-246, i.e., reactivation of mutant p53 and modulation of redox activity. SLC7A11, TRIM16, TXNRD1, and SRXN1 are potential new pharmacodynamic biomarkers for assessing the response to APR-246 in both preclinical and clinical studies.     

Molecular profile of breast versus ovarian cancer cells in response to treatment with the anticancer drugs cisplatin, carboplatin, doxorubicin, etoposide and taxol

We assessed changes in the apoptosis-related genes BCL2, BAX, BCL2L12, FAS and CASPASE-3 in OVCAR-3 human ovarian cancer cells and BT-20 human breast cancer cells to provide an insight into the molecular mechanisms involved in the response of these cells to treatment with anticancer drugs and to assess their value as potential biomarkers of chemotherapy response in breast and ovarian cancer. Cells were treated with different chemotherapeutic drugs (cisplatin, carboplatin, doxorubicin, etoposide and taxol) and assessed for changes in the expression of apoptosis-related genes at the mRNA level. Total RNA was extracted, reverse-transcribed into cDNA and amplified by PCR using gene-specific primers. GAPDH was used as a housekeeping gene. Cytotoxicity was assessed by MTT assay. Both cancer cell lines responded differentially at the molecular level to the drug treatments. OVCAR-3 cells showed more pronounced sensitivity and changes compared to BT-20 cells at the mRNA level for different apoptosis-related genes, leading to cell and cancer type dependence in conjunction with drug dependence.     

Down-regulation of HtrA1 activates the epithelial-mesenchymal transition and ATM DNA damage response pathways

Expression of the serine protease HtrA1 is decreased or abrogated in a variety of human primary cancers, and higher levels of HtrA1 expression are directly related to better response to chemotherapeutics. However, the precise mechanisms leading to HtrA1 down regulation during malignant transformation are unclear. To investigate HtrA1 gene regulation in breast cancer, we characterized expression in primary breast tissues and seven human breast epithelial cell lines, including two non-tumorigenic cell lines. In human breast tissues, HtrA1 expression was prominent in normal ductal glands. In DCIS and in invasive cancers, HtrA1 expression was greatly reduced or lost entirely. HtrA1 staining was also reduced in all of the human breast cancer cell lines, compared with the normal tissue and non-tumorigenic cell line controls. Loss of HtrA1 gene expression was attributable primarily to epigenetic silencing mechanisms, with different mechanisms operative in the various cell lines. To mechanistically examine the functional consequences of HtrA1 loss, we stably reduced and/or overexpressed HtrA1 in the non-tumorigenic MCF10A cell line. Reduction of HtrA1 levels resulted in the epithelial-to-mesenchymal transition with acquisition of mesenchymal phenotypic characteristics, including increased growth rate, migration, and invasion, as well as expression of mesenchymal biomarkers. A concomitant decrease in expression of epithelial biomarkers and all microRNA 200 family members was also observed. Moreover, reduction of HtrA1 expression resulted in activation of the ATM and DNA damage response, whereas overexpression of HtrA1 prevented this activation. Collectively, these results suggest that HtrA1 may function as a tumor suppressor by controlling the epithelial-to-mesenchymal transition, and may function in chemotherapeutic responsiveness by mediating DNA damage response pathways.     

Identification of circadian-related gene expression profiles in entrained breast cancer cell lines

Cancer cells have broken circadian clocks when compared to their normal tissue counterparts. Moreover, it has been shown in breast cancer that disruption of common circadian oscillations is associated with a more negative prognosis. Numerous studies, focused on canonical circadian genes in breast cancer cell lines, have suggested that there are no mRNA circadian-like oscillations. Nevertheless, cancer cell lines have not been extensively characterized and it is unknown to what extent the circadian oscillations are disrupted. We have chosen representative non-cancerous and cancerous breast cell lines (MCF-10A, MCF-7, ZR-75-30, MDA-MB-231 and HCC-1954) in order to determine the degree to which the circadian clock is damaged. We used serum shock to synchronize the circadian clocks in culture. Our aim was to initially observe the time course of gene expression using cDNA microarrays in the non-cancerous MCF-10A and the cancerous MCF-7 cells for screening and then to characterize specific genes in other cell lines. We used a cosine function to select highly correlated profiles. Some of the identified genes were validated by quantitative polymerase chain reaction (qPCR) and further evaluated in the other breast cancer cell lines. Interestingly, we observed that breast cancer and non-cancerous cultured cells are able to generate specific circadian expression profiles in response to the serum shock. The rhythmic genes, suggested via microarray and measured in each particular subtype, suggest that each breast cancer cell type responds differently to the circadian synchronization. Future results could identify circadian-like genes that are altered in breast cancer and non-cancerous cells, which can be used to propose novel treatments. Breast cell lines are potential models for in vitro studies of circadian clocks and clock-controlled pathways.     

Integrative genomics revealed RAI3 is a cell growth-promoting gene and a novel P53 transcriptional target

In this study, differential gene expression between normal human mammary epithelial cells and their malignant counterparts (eight well established breast cancer cell lines) was studied using Incyte GeneAlbum 1-6, which contains 65,873 cDNA clones representing 33,515 individual genes. 3,152 cDNAs showed a > or =3.0-fold expression level change in at least one of the human breast cancer cell lines as compared with normal human mammary epithelial cells. Integration of breast tumor gene expression data with the genes in the tumor suppressor p53 signaling pathway yielded 128 genes whose expression is altered in breast tumor cell lines and in response to p53 expression. A hierarchical cluster analysis of the 128 genes revealed that a significant portion of genes demonstrate an opposing expression pattern, i.e. p53-activated genes are down-regulated in the breast tumor lines, whereas p53-repressed genes are up-regulated. Most of these genes are involved in cell cycle regulation and/or apoptosis, consistent with the tumor suppressor function of p53. Follow-up studies on one gene, RAI3, suggested that p53 interacts with the promoter of RAI3 and repressed its expression at the onset of apoptosis. The expression of RAI3 is elevated in most tumor cell lines expressing mutant p53, whereas RAI3 mRNA is relatively repressed in the tumor cell lines expressing wild-type p53. Furthermore, ectopic expression of RAI3 in 293 cells promotes anchorage-independent growth and small interfering RNA-mediated depletion of RAI3 in AsPc-1 pancreatic tumor cells induces cell morphological change. Taken together, these data suggest a role for RAI3 in tumor growth and demonstrate the predictive power of integrative genomics.     

Gene Expression Profiling Reveals Epithelial Mesenchymal Transition (EMT) Genes Can Selectively Differentiate Eribulin Sensitive Breast Cancer Cells

Objectives

Eribulin mesylate is a synthetic macrocyclic ketone analog of the marine sponge natural product halichondrin B. Eribulin is a mechanistically unique inhibitor of microtubule dynamics. In this study, we investigated whether selective signal pathways were associated with eribulin activity compared to paclitaxel, which stabilizes microtubules, based on gene expression profiling of cell line panels of breast, endometrial, and ovarian cancer in vitro.

Results

We determined the sets of genes that were differentially altered between eribulin and paclitaxel treatment in breast, endometrial, and ovarian cancer cell line panels. Our unsupervised clustering analyses revealed that expression profiles of gene sets altered with treatments were correlated with the in vitro antiproliferative activities of the drugs. Several tubulin isotypes had significantly lower expression in cell lines treated with eribulin compared to paclitaxel. Pathway enrichment analyses of gene sets revealed that the common pathways altered between treatments in the 3 cancer panels were related to cytoskeleton remodeling and cell cycle regulation. The epithelial-mesenchymal transition (EMT) pathway was enriched in genes with significantly altered expression between the two drugs for breast and endometrial cancers, but not for ovarian cancer. Expression of genes from the EMT pathway correlated with eribulin sensitivity in breast cancer and with paclitaxel sensitivity in endometrial cancer. Alteration of expression profiles of EMT genes between sensitive and resistant cell lines allowed us to predict drug sensitivity for breast and endometrial cancers.

Conclusion

Gene expression analysis showed that gene sets that were altered between eribulin and paclitaxel correlated with drug in vitro antiproliferative activities in breast and endometrial cancer cell line panels. Among the panels, breast cancer provided the strongest differentiation between eribulin and paclitaxel sensitivities based on gene expression. In addition, EMT genes were predictive of eribulin sensitivity in the breast and endometrial cancer panels.     

Differential gene expression in normal and transformed human mammary epithelial cells in response to oxidative stress

Oxidative stress plays a key role in breast carcinogenesis. To investigate whether normal and malignant breast epithelial cells differ in their responses to oxidative stress, we examined the global gene expression profiles of three cell types, representing cancer progression from a normal to a malignant stage, under oxidative stress. Normal human mammary epithelial cells (HMECs), an immortalized cell line (HMLER-1), and a tumorigenic cell line (HMLER-5) were exposed to increased levels of reactive oxygen species (ROS) by treatment with glucose oxidase. Functional analysis of the metabolic pathways enriched with differentially expressed genes demonstrated that normal and malignant breast epithelial cells diverge substantially in their response to oxidative stress. Whereas normal cells exhibit the up-regulation of antioxidant mechanisms, cancer cells are unresponsive to the ROS insult. However, the gene expression response of normal HMECs under oxidative stress is comparable to that of the malignant cells under normal conditions, indicating that altered redox status is persistent in breast cancer cells, which makes them resistant to increased generation of ROS. We discuss some of the possible adaptation mechanisms of breast cancer cells under persistent oxidative stress that differentiate them from normal mammary epithelial cells as regards the response to acute oxidative stress.     

应用基因表达谱芯片从不同转移能力的 乳腺癌细胞中筛选转移相关基因

人乳腺癌细胞系 MCF-7 及其转移亚克隆 LM-MCF-7 为肿瘤转移分子机制的研究提供了细胞模型 . 应用基因芯片技术比较两种具有不同转移能力细胞系基因表达谱的差异,寻找乳腺癌转移相关基因 . 提取两种细胞总 RNA ,分别用 Cy5-dCTP 、 Cy3-dCTP 标记 LM-MCF-7 和 MCF-7 的 cDNA ,并与含有 21 329 个基因的芯片进行杂交并扫描,利用 GenePix Pro 4.0 图像分析软件处理数据判断基因是否在两个细胞中存在表达差异 . 经互换荧光标记物重复两次实验,共筛选出差异表达基因 67 个,其中 41 个在 LM-MCF-7 细胞中表达上调, 26 个在 LM-MCF-7 细胞中表达下调 . 应用实时定量 RT-PCR 对 7 个表达差异明显的基因进行了验证 . 生物信息学分析结果提示,上述发现的差异基因编码产物与细胞内信号转导、转录调节、应激反应、新陈代谢、发育、细胞运动、细胞凋亡和细胞粘连等功能有关 . 据文献报道,这些差异表达的基因中有 35 个与肿瘤有关,其中 9 个与乳腺癌转移有关, 6 个可能参与肿瘤浸润和转移过程 . 根据基因芯片检测的结果,从功能上对 LM-MCF-7 细胞和 MCF-7 细胞与细胞凋亡的关系进行了研究,发现具有高转移倾向的 LM-MCF-7 细胞与 MCF-7 细胞相比,抗凋亡能力较强 . 上述与肿瘤转移相关基因在肿瘤转移中的作用及其分子机理有待深入研究 .     

干扰小RNA序列非依赖性地影响胰腺癌细胞的基因表达

干扰小RNA(small interference RNA, siRNA)是基因敲减的常用工具,广泛用于基因沉默技术和基因功能研究,在临床疾病治疗等方面也有潜在的应用。一般认为,达到一定长度（比如大于27 bp）的双链RNA可以诱导干扰素反应,降低相关基因的表达。目前,siRNA对基因表达的非特异性作用尚不完全清楚。为研究siRNA干扰的非特异基因表达,本研究以胰腺癌细胞HPAC 和BxPC3 为模型,采用高通量测序技术对6种不同干扰小RNA处理及未作处理的HPAC和BxPC3细胞进行转录组测序分析,筛选出干扰小RNA处理后表达量共同下调的基因进行研究。通过生物信息学方法对表达下调基因的功能进行研究,并利用实时荧光定量PCR(qRT-PCR)技术对部分下调基因进行验证。结果表明,短片段双链小RNA能够显著改变细胞的基因表达,而这些基因表达谱的变化是有规律的,特定功能的基因优先发生变化。在表达下调的基因中,某些特定类型的基因变化非常显著,包括氨基酸代谢相关基因、Hedgehog信号途径基因和多巴胺受体D5基因等。这些结果表明,在使用siRNA时需要考虑其序列非依赖性地基因表达调控作用。     

Expression and prognostic significance of the autoimmune regulator gene in breast cancer cells

The autoimmune regulator gene (AIRE) plays a fundamental role in tolerance by promoting the expression of tissue-specific antigens in medullary thymic epithelial cells (mTECs). Recently, AIRE expression was detected also in human keratinocytes and in tumors originating in stratified epithelia. Here, we tested whether AIRE is expressed in cancer cells. We analyzed AIRE expression in cancer cases from The Cancer Genome Atlas (TCGA) RNA-seq dataset and we found association with better outcome. AIRE protein expression was verified by immunohistochemistry in a cohort of 39 human breast cancer specimens and its prognostic relevance was confirmed in microarray-based gene expression data set NKI-295 and KM-Plotter. Both in the RNA-seq and gene expression datasets analyzed, AIRE expression was an independent strong prognostic factor for relapse-free survival (RFS), particularly in estrogen receptor-positive tumors. Enrichment of translation-related pathways was observed in AIRE-expressing tumors by Ingenuity Pathway Analysis and a significant increase of cells in G1 phase and activation of caspase cascades was induced by AIRE transfection in breast cancer luminal cell lines, suggesting that AIRE-induced over-translation of proteins lead to cycle arrest and apoptosis. These data are the first to identify AIRE expression in breast cancer and an association with prognosis.