E2F1 Drives Breast Cancer Metastasis by Regulating the Target Gene FGF13 and Altering Cell Migration

Up-regulation of NCAPG mediated by E2F1 facilitates the progression of osteosarcoma through the Wnt/β-catenin signaling pathway

Background

In recent years, there are few reports on non-SMC condensin I complex subunit G (NCAPG) in osteosarcoma. Our study aims to explore the biological role of NCAPG in osteosarcoma and its underlying molecular mechanism and to further clarify the reasons for the abnormal expression of NCAPG in osteosarcoma.

Methods

Here, we mined The Cancer Genome Atlas (TCGA) Program public database through bioinformatics methods, analyzed the differential expression of NCAPG in sarcoma tissue and normal tissue, and explored the relationship between NCAPG expression level and sarcoma tissue differentiation, including tumor recurrence, metastasis, and patient survival. Next, the transcription factors responsible for the abnormal expression of NCAPG in osteosarcoma tumors were predicted by multiple online website tools and verified via cellular experiments. Subsequently, loss of function and cell phenotype experiments were performed to confirm the effect of NCAPG on the malignant biological behavior of osteosarcoma cells. Mechanistically, by reviewing the literature, we found that NCAPG can affect the malignant progression of many solid tumors by regulating the Wnt/β-catenin signaling pathway. Therefore, we preliminarily investigated the potential effect of NCAPG on this pathway via western blot experiments in osteosarcoma.

Results

Increased expression of NCAPG was found in sarcoma compared to normal tissues, which was positively correlated with poor differentiation, metastasis, and poor prognosis. Combining the transcription factor prediction results, correlation analysis, and expression level in the TCGA public database with validation outcomes of in vitro cell assays, we found that E2F transcription factor 1 (E2F1) regulated the increased expression of NCAPG in osteosarcoma. The results of cell phenotype experiments showed that silencing NCAPG could inhibit the proliferation, migration, and invasion of osteosarcoma cells. The preliminary mechanistic investigation suggested that NCAPG may affect osteosarcoma progression through the Wnt/β-catenin pathway.

Conclusions

Our data reveal that E2F1 facilitates NCAPG expression in osteosarcoma by regulating the transcription of the NCAPG gene. Up-regulation of NCAPG promotes osteosarcoma progression via the Wnt/β-catenin signaling axis.

Decoding the ubiquitin language: Orchestrating transcription initiation and gene expression through chromatin remodelers and histones

Eukaryotic transcription is a finely orchestrated process and it is controlled by transcription factors as well as epigenetic regulators. Transcription factors and epigenetic regulators undergo different types of posttranslational modifications including ubiquitination to control transcription process. Ubiquitination, traditionally associated with protein degradation, has emerged as a crucial contributor to the regulation of chromatin structure through ubiquitination of histone and chromatin remodelers. Ubiquitination introduces new layers of intricacy to the regulation of transcription initiation through controlling the equilibrium between euchromatin and heterochromatin states. Nucleosome, the fundamental units of chromatin, spacing in euchromatin and heterochromatin states are regulated by histone modification and chromatin remodeling complexes. Chromatin remodeling complexes actively sculpt the chromatin architecture and thereby influence the transcriptional states of genes. Therefore, understanding the dynamic behavior of nucleosome spacing is critical as it impacts various cellular functions through controlling gene expression profiles. In this comprehensive review, we discussed the intricate interplay between ubiquitination and transcription initiation, and illuminated the underlying molecular mechanisms that occur in a variety of biological contexts. This exploration sheds light on the complex regulatory networks that govern eukaryotic transcription, providing important insights into the fine orchestration of gene expression and chromatin dynamics.

Blocking the E2F transcription factor 1/high-mobility group box 2 pathway enhances the intervention effects of α-santalol on the malignant behaviors of liver cancer cells

In view of the tumor-inhibiting effect of α-santalol in various cancers and the role of E2F transcription factor 1 (E2F1) as an important target for anticancer research, this study investigates the relation between α-santalol and E2F1, as well as the effect of α-santalol on liver cancer progression and the corresponding mechanism. Concretely, liver cancer cells were treated with different concentrations of α-santalol. The IC50 value of α-santalol was determined using Probit regression analysis. Then, transcription factors that are targeted by α-santalol and differentially expressed in liver cancer were screened out. The clinicopathological impact of E2F1 and its targets were evaluated and predicted. The expressions of E2F1 and high-mobility group box 2 (HMGB2) and their correlation in the liver cancer tissues were analyzed by bioinformatics. The effects of E2F1 and HMGB2 on the biological characteristics of liver cancer cells were examined through loss/gain-of-function and molecular assays. With the extension of treatment time, the inhibitory effects of 10 μmol/L and 20 μmol/L α-santalol on cancer cell survival rate were enhanced (P < 0.001). E2F1 and HMGB2 were highly expressed and positively correlated in liver cancer tissues (P < 0.05). High E2F1 expression was correlated with large tumors and high TNM stages (P < 0.05). E2F1 knockdown promoted the effects of α-santalol on dose-dependently inhibiting viability, colony formation, invasion and migration (P < 0.05). Moreover, E2F1 knockdown reduced the IC50 value and HMGB2 level, while HMGB2 overexpression produced opposite effects. HMGB2 overexpression and E2F1 knockdown mutually counteracted their effects on the IC50 value and on the viability and apoptosis of α-santalol-treated liver cancer cells (P < 0.01). Collectively, blocking the E2F1/HMGB2 pathway enhances the intervention effects of α-santalol on the proliferation, migration and invasion of liver cancer cells.

Molecular targets and therapeutic strategies for triple-negative breast cancer

Triple-negative breast cancer (TNBC) is known for its heterogeneous complexity and is often difficult to treat. TNBC lacks the expression of major hormonal receptors like estrogen receptor, progesterone receptor, and human epidermal growth factor receptor-2 and is further subdivided into androgen receptor (AR) positive and AR negative. In contrast, AR negative is also known as quadruple-negative breast cancer (QNBC). Compared to AR-positive TNBC, QNBC has a great scarcity of prognostic biomarkers and therapeutic targets. QNBC shows excessive cellular growth and proliferation of tumor cells due to increased expression of growth factors like EGF and various surface proteins. This study briefly reviews the limited data available as protein biomarkers that can be used as molecular targets in treating TNBC as well as QNBC. Targeted therapy and immune checkpoint inhibitors have recently changed cancer treatment. Many studies in medicinal chemistry continue to focus on the synthesis of novel compounds to discover new antiproliferative medicines capable of treating TNBC despite the abundance of treatments currently on the market. Drug repurposing is one of the therapeutic methods for TNBC that has been examined. Moreover, some additional micronutrients, nutraceuticals, and functional foods may be able to lower cancer risk or slow the spread of malignant diseases that have already been diagnosed with cancer. Finally, nanomedicines, or applications of nanotechnology in medicine, introduce nanoparticles with variable chemistry and architecture for the treatment of cancer. This review emphasizes the most recent research on nutraceuticals, medication repositioning, and novel therapeutic strategies for the treatment of TNBC.

Depleted hexokinase1 and lack of AMPKα activation favor OXPHOS-dependent energetics in retinoblastoma tumors

Lack of retinoblastoma (Rb) protein causes aggressive intraocular retinal tumors in children. Recently, Rb tumors have been shown to have a distinctly altered metabolic phenotype, such as reduced expression of glycolytic pathway proteins alongside altered pyruvate and fatty acid levels. In this study, we demonstrate that loss of hexokinase 1(HK1) in tumor cells rewires their metabolism allowing enhanced oxidative phosphorylation-dependent energy production. We show that rescuing HK1 or retinoblastoma protein 1 (RB1) in these Rb cells reduced cancer hallmarks such as proliferation, invasion, and spheroid formation and increased their sensitivity to chemotherapy drugs. Induction of HK1 was accompanied by a metabolic shift of the cells to glycolysis and a reduction in mitochondrial mass. Cytoplasmic HK1 bound Liver Kinase B1 and phosphorylated AMP-activated kinase-α (AMPKα Thr172), thereby reducing mitochondria-dependent energy production. We validated these findings in tumor samples from Rb patients compared to age-matched healthy retinae. HK1 or RB1 expression in Rb-/- cells led to a reduction in their respiratory capacity and glycolytic proton flux. HK1 overexpression reduced tumor burden in an intraocular tumor xenograft model. AMPKα activation by AICAR also enhanced the tumoricidal effects of the chemotherapeutic drug topotecan in vivo. Therefore, enhancing HK1 or AMPKα activity can reprogram cancer metabolism and sensitize Rb tumors to lower doses of existing treatments, a potential therapeutic modality for Rb.

Immuno-oncologic signature of malignant transformation in oral squamous cell carcinoma

OBJECTIVE

The purpose of this study is to identify the immuno-oncologic (IO) signature at the surgical tumor margin (TM) of oral squamous cell carcinoma (OSCC) that is involved in the process of malignant transformation.

STUDY DESIGN

Under institutional review board approval, TM of 73 OSCC were investigated using immunohistochemistry for the immune biomarker, programmed death ligand-1 (PD-L1). NanoString 770 IO-focused gene set was analyzed in 5 pairs of TM and invasive tumor (T). PD-L1 regulation in response to interferon-gamma (IFN-γ) was investigated in an oral potentially malignant cell line (OPMC).

RESULTS

Programmed death ligand-1 expression in the epithelial margin directly correlated with its expression in the underlying immune cells (P = .0082). Differential gene expression showed downregulation of PD-L1 and IFN-γ 6 gene signature in the TM relative to T pair.CD8 and macrophages were higher in TM. CNTFR, LYZ, C7, RORC, and FGF13 downregulation in T relative to TM. TDO2, ADAM12, MMP1, LAMC2, MB21D1, TYMP, OASL, COL5A1, exhausted_CD8, Tregs,and NK_CD56dim were upregulated in T relative to TM. Finally, IFN-γ induced upregulation of PD-L1 in the OPMC.

CONCLUSIONS

Our work suggests a role for IFN-γ in PD-L1 upregulation in OPMC and presents novel IO transcriptional signatures for frankly invasive OSCC relative to TM.

Architecture and topologies of gene regulatory networks associated with breast cancer, adjacent normal, and normal tissues

Most cancer studies employ adjacent normal tissues to tumors (ANTs) as controls, which are not completely normal and represent a pre-cancerous state. However, the regulatory landscape of ANTs compared to tumor and non-tumor-bearing normal tissues is largely unexplored. Among cancers, breast cancer is the most commonly diagnosed cancer and a leading cause of death in women worldwide, with a lack of sufficient treatment regimens for various reasons. Hence, we aimed to gain deeper insights into normal, pre-cancerous, and cancerous regulatory systems of breast tissues towards identifying ANT and subtype-specific candidate genes. For this, we constructed and analyzed eight gene regulatory networks (GRNs), including five subtypes (viz., Basal, Her2, Luminal A, Luminal B, and Normal-Like), one ANT, and two normal tissue networks. Whereas several topological properties of these GRNs enabled us to identify tumor-related features of ANT, escape velocity centrality (EVC+) identified 24 functionally significant common genes, including well-known genes such as E2F1, FOXA1, JUN, BRCA1, GATA3, ERBB2, and ERBB3 across all six tissues including subtypes and ANT. Similarly, the EVC+ also helped us to identify tissue-specific key genes (Basal: 18, Her2: 6, Luminal A: 5, Luminal B: 5, Normal-Like: 2, and ANT: 7). Additionally, differentially correlated switching gene pairs along with functional, pathway, and disease annotations highlighted the cancer-associated role of these genes. In a nutshell, the present study revealed ANT and subtype-specific regulatory features and key candidate genes, which can be explored further using in vitro and in vivo experiments for better and effective disease management at an early stage.

Integrative multi-omic sequencing reveals the MMTV-Myc mouse model mimics human breast cancer heterogeneity

BACKGROUND

Breast cancer is a complex and heterogeneous disease with distinct subtypes and molecular profiles corresponding to different clinical outcomes. Mouse models of breast cancer are widely used, but their relevance in capturing the heterogeneity of human disease is unclear. Previous studies have shown the heterogeneity at the gene expression level for the MMTV-Myc model, but have only speculated on the underlying genetics.

METHODS

Tumors from the microacinar, squamous, and EMT histological subtypes of the MMTV-Myc mouse model of breast cancer underwent whole genome sequencing. The genomic data obtained were then integrated with previously obtained matched sample gene expression data and extended to additional samples of each histological subtype, totaling 42 gene expression samples. High correlation was observed between genetic copy number events and resulting gene expression by both Spearman's rank correlation coefficient and the Kendall rank correlation coefficient. These same genetic events are conserved in humans and are indicative of poor overall survival by Kaplan-Meier analysis. A supervised machine learning algorithm trained on METABRIC gene expression data was used to predict the analogous human breast cancer intrinsic subtype from mouse gene expression data.

RESULTS

Herein, we examine three common histological subtypes of the MMTV-Myc model through whole genome sequencing and have integrated these results with gene expression data. Significantly, key genomic alterations driving cell signaling pathways were well conserved within histological subtypes. Genomic changes included frequent, co-occurring mutations in KIT and RARA in the microacinar histological subtype as well as SCRIB mutations in the EMT subtype. EMT tumors additionally displayed strong KRAS activation signatures downstream of genetic activating events primarily ascribed to KRAS activating mutations, but also FGFR2 amplification. Analogous genetic events in human breast cancer showed stark decreases in overall survival. In further analyzing transcriptional heterogeneity of the MMTV-Myc model, we report a supervised machine learning model that classifies MMTV-Myc histological subtypes and other mouse models as being representative of different human intrinsic breast cancer subtypes.

CONCLUSIONS

We conclude the well-established MMTV-Myc mouse model presents further opportunities for investigation of human breast cancer heterogeneity.

Methylation of the transcription factor E2F1 by SETD6 regulates SETD6 expression via a positive feedback mechanism

The protein lysine methyltransferase SET domain-containing protein 6 (SETD6) has been shown to influence different cellular activities and to be critically involved in the regulation of diverse developmental and pathological processes. However, the upstream signals that regulate the mRNA expression of SETD6 are not known. Bioinformatic analysis revealed that the SETD6 promoter has a binding site for the transcription factor E2F1. Using various experimental approaches, we show that E2F1 binds to the SETD6 promoter and regulates SETD6 mRNA expression. Our further observation that this phenomenon is SETD6 dependent suggested that SETD6 and E2F1 are linked. We next demonstrate that SETD6 monomethylates E2F1 specifically at K117 in vitro and in cells. Finally, we show that E2F1 methylation at K117 positively regulates the expression level of SETD6 mRNA. Depletion of SETD6 or overexpression of E2F1 K117R mutant, which cannot be methylated by SETD6, reverses the effect. Taken together, our data provide evidence for a positive feedback mechanism, which regulates the expression of SETD6 by E2F1 in a SETD6 methylation-dependent manner, and highlight the importance of protein lysine methyltransferases and lysine methylation signaling in the regulation of gene transcription.

Thyroid hormone enhances estrogen-mediated proliferation and cell cycle regulatory pathways in steroid receptor-positive breast Cancer

Estrogen receptor (ER) α expression and associated signaling is a major driver of over two-thirds of all breast cancers (BC). ER targeting strategies are typically used as a first-line therapy in patients with steroid receptor positive (SR+) disease. Secondary resistance to anti-estrogenic agents may occur with clonal expansion and disease progression. Mechanisms underlying hormone resistance are an expanding field of significant translational importance. Cross-talk with other nuclear hormones, receptors, and signaling pathways, including thyroid hormones (TH) and their receptors (THRs), have been shown to promote endocrine therapy resistance in some studies. We have shown that TH replacement therapy (THRT) was independently and significantly associated with higher rates of relapse and mortality in SR positive (+), node-negative (LN-) BC patients, whereas it showed no association with outcomes in SR negative (-) patients. LN-, SR+ patients receiving THRT and tamoxifen had the worst outcomes, suggesting a pro-carcinogenic interaction that significantly and independently shortened survival and increased mortality. Using in vivo and in vitro models, we previously showed hormonal cross-talk, altered gene signaling, target gene activation, and resistance to tamoxifen in the presence of TH. In this report, we show TH ± E2 ± tamoxifen inhibits cell cycle control signaling, reduces apoptosis, and enhances cell proliferation, tumor growth, tamoxifen resistance, and clonal expansion. Mechanistically these changes involve numerous genes and pathways, including critical cell cycle regulatory proteins and genes identified using various molecular methods. These studies facilitate a greater mechanistic understanding of the biological and molecular impact of TH on SR+ BC.

Cancer invasion and metastasis: Insights from murine pubertal mammary gland morphogenesis

Cancer invasion and metastasis accounts for the majority of cancer related mortality. A better understanding of the players that drive the aberrant invasion and migration of tumors cells will provide critical targets to inhibit metastasis. Postnatal pubertal mammary gland morphogenesis is characterized by highly proliferative, invasive, and migratory normal epithelial cells. Identifying the molecular regulators of pubertal gland development is a promising strategy since tumorigenesis and metastasis is postulated to be a consequence of aberrant reactivation of developmental stages. In this review, we summarize the pubertal morphogenesis regulators that are involved in cancer metastasis and revisit pubertal mammary gland transcriptome profiling to uncover both known and unknown metastasis genes. Our updated list of pubertal morphogenesis regulators shows that most are implicated in invasion and metastasis. This review highlights molecular linkages between development and metastasis and provides a guide for exploring novel metastatic drivers.

EZH2 regulates pancreatic cancer cells through E2F1, GLI1, CDK3, and Mcm4

Pancreatic cancer (PC) is one of the most common malignant tumors in digestive tract. To explore the role of epigenetic factor EZH2 in the malignant proliferation of PC, so as to provide effective medical help in PC. Sixty paraffin sections of PC were collected and the expression of EZH2 in PC tissues was detected by immunohistochemical assay. Three normal pancreas tissue samples were used as controls. The regulation of EZH2 gene on proliferation and migration of normal pancreatic cell and PC cell were determined by MTS, colony forming, Ki-67 antibody, scratch and Transwell assays. Through differential gene annotation and differential gene signaling pathway analysis, differentially expressed genes related to cell proliferation were selected and verified by RT-qPCR. EZH2 is mainly expressed in the nuclei of pancreatic tumor cells, but not in normal pancreatic cells. The results of cell function experiments showed that EZH2 overexpression could enhance the proliferation and migration ability of PC cell BXPC-3. Cell proliferation ability increased by 38% compared to the control group. EZH2 knockdown resulted in reduced proliferation and migration ability of cells. Compared with control, proliferation ability of cells reduced by 16%-40%. The results of bioinformatics analysis of transcriptome data and RT-qPCR demonstrated that EZH2 could regulate the expression of E2F1, GLI1, CDK3 and Mcm4 in normal and PC cells. The results revealed that EZH2 might regulate the proliferation of normal pancreatic cell and PC cell through E2F1, GLI1, CDK3 and Mcm4.

SOX2 expression in prostate cancer drives resistance to nuclear hormone receptor signaling inhibition through the WEE1/CDK1 signaling axis

The development of androgen receptor signaling inhibitor (ARSI) drug resistance in prostate cancer (PC) remains therapeutically challenging. Our group has described the role of sex determining region Y-box 2 (SOX2) overexpression in ARSI-resistant PC. Continuing this work, we report that NR3C1, the gene encoding glucocorticoid receptor (GR), is a novel SOX2 target in PC, positively regulating its expression. Similar to ARSI treatment, SOX2-positive PC cells are insensitive to GR signaling inhibition using a GR modulating therapy. To understand SOX2-mediated nuclear hormone receptor signaling inhibitor (NHRSI) insensitivity, we performed RNA-seq in SOX2-positive and -negative PC cells following NHRSI treatment. RNA-seq prioritized differentially regulated genes mediating the cell cycle, including G2 checkpoint WEE1 Kinase (WEE1) and cyclin-dependent kinase 1 (CDK1). Additionally, WEE1 and CDK1 were differentially expressed in PC patient tumors dichotomized by high vs low SOX2 gene expression. Importantly, pharmacological targeting of WEE1 (WEE1i) in combination with an ARSI or GR modulator re-sensitizes SOX2-positive PC cells to nuclear hormone receptor signaling inhibition in vitro, and WEE1i combined with ARSI significantly slowed tumor growth in vivo. Collectively, our data suggest SOX2 predicts NHRSI resistance, and simultaneously indicates the addition of WEE1i to improve therapeutic efficacy of NHRSIs in SOX2-positive PC.

Multimodal AI for prediction of distant metastasis in carcinoma patients

Metastasis of cancer is directly related to death in almost all cases, however a lot is yet to be understood about this process. Despite advancements in the available radiological investigation techniques, not all cases of Distant Metastasis (DM) are diagnosed at initial clinical presentation. Also, there are currently no standard biomarkers of metastasis. Early, accurate diagnosis of DM is however crucial for clinical decision making, and planning of appropriate management strategies. Previous works have achieved little success in attempts to predict DM from either clinical, genomic, radiology, or histopathology data. In this work we attempt a multimodal approach to predict the presence of DM in cancer patients by combining gene expression data, clinical data and histopathology images. We tested a novel combination of Random Forest (RF) algorithm with an optimization technique for gene selection, and investigated if gene expression pattern in the primary tissues of three cancer types (Bladder Carcinoma, Pancreatic Adenocarcinoma, and Head and Neck Squamous Carcinoma) with DM are similar or different. Gene expression biomarkers of DM identified by our proposed method outperformed Differentially Expressed Genes (DEGs) identified by the DESeq2 software package in the task of predicting presence or absence of DM. Genes involved in DM tend to be more cancer type specific rather than general across all cancers. Our results also indicate that multimodal data is more predictive of metastasis than either of the three unimodal data tested, and genomic data provides the highest contribution by a wide margin. The results re-emphasize the importance for availability of sufficient image data when a weakly supervised training technique is used. Code is made available at: https://github.com/rit-cui-lab/Multimodal-AI-for-Prediction-of-Distant-Metastasis-in-Carcinoma-Patients.

Epigenetic alterations impede epithelial-mesenchymal transition by modulating centrosome amplification and Myc/RAS axis in triple negative breast cancer cells

Alterations in centrosome proteins may result in centrosome abnormalities such as disorganized spindles and centrosome amplification, leading to aneuploidy and genomic instability. Centrosomes exhibit unique epigenetic properties in which structural or positional information is propagated through somatic lineage by non-genetic pathways. Excessive centrosome amplification in breast cancer is accompanied by efficient clustering and loss of E-cadherin, indicating an important adaptive mechanism of cancer. This study sought to elucidate the effect of epigenetic alterations on centrosome amplification, epithelial-mesenchymal transition (EMT) and apoptosis in triple negative human breast adenocarcinoma derived MDA-MB-231 cell line. The results obtained here show that siRNA mediated silencing of DNMT1 and specific inhibition of HDAC1 & HDAC2 by Tricostatin A (TSA) synergistically inhibit cell proliferation through modulation of centrosome proteins γ-tubulin, TUBGCP2 and pericentrin. In addition, induction of apoptosis was observed by downregulation of Bcl2, upregulation of Bax and activation of PARP cleavage. Inhibition of EMT was confirmed through upregulation of E-cadherin and downregulation of N-cadherin and vimentin. Similarly, downregulation of Myc, RAS and CDK2, which plays important roles in proliferation and survival, was observed. Nuclear protein analysis revealed downregulation in the nuclear translocation of E2F1, which regulates centrosome amplification and metastasis in breast cancer. In conclusion, this study confirmed the role of epigenetic regulators in centrosome amplification and suggests that inhibition of DNA methylation and histone deacetylation-mediated chromatin remodelling synergistically disrupt EMT through modulation of centrosome amplification and Myc/RAS axis to potentiate apoptosis and attenuate cell proliferation in triple negative breast cancer cells.

Deregulated E2F Activity as a Cancer-Cell Specific Therapeutic Tool

The transcription factor E2F, the principal target of the tumor suppressor pRB, plays crucial roles in cell proliferation and tumor suppression. In almost all cancers, pRB function is disabled, and E2F activity is enhanced. To specifically target cancer cells, trials have been undertaken to suppress enhanced E2F activity to restrain cell proliferation or selectively kill cancer cells, utilizing enhanced E2F activity. However, these approaches may also impact normal growing cells, since growth stimulation also inactivates pRB and enhances E2F activity. E2F activated upon the loss of pRB control (deregulated E2F) activates tumor suppressor genes, which are not activated by E2F induced by growth stimulation, inducing cellular senescence or apoptosis to protect cells from tumorigenesis. Deregulated E2F activity is tolerated in cancer cells due to inactivation of the ARF-p53 pathway, thus representing a feature unique to cancer cells. Deregulated E2F activity, which activates tumor suppressor genes, is distinct from enhanced E2F activity, which activates growth-related genes, in that deregulated E2F activity does not depend on the heterodimeric partner DP. Indeed, the ARF promoter, which is specifically activated by deregulated E2F, showed higher cancer-cell specific activity, compared to the E2F1 promoter, which is also activated by E2F induced by growth stimulation. Thus, deregulated E2F activity is an attractive potential therapeutic tool to specifically target cancer cells.

Identifying Associations between DCE-MRI Radiomic Features and Expression Heterogeneity of Hallmark Pathways in Breast Cancer: A Multi-Center Radiogenomic Study

BACKGROUND

To investigate the relationship between dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) radiomic features and the expression activity of hallmark pathways and to develop prediction models of pathway-level heterogeneity for breast cancer (BC) patients.

METHODS

Two radiogenomic cohorts were analyzed (n = 246). Tumor regions were segmented semiautomatically, and 174 imaging features were extracted. Gene set enrichment analysis (GSEA) and gene set variation analysis (GSVA) were performed to identify significant imaging-pathway associations. Random forest regression was used to predict pathway enrichment scores. Five-fold cross-validation and grid search were used to determine the optimal preprocessing operation and hyperparameters.

RESULTS

We identified 43 pathways, and 101 radiomic features were significantly related in the discovery cohort (p-value < 0.05). The imaging features of the tumor shape and mid-to-late post-contrast stages showed more transcriptional connections. Ten pathways relevant to functions such as cell cycle showed a high correlation with imaging in both cohorts. The prediction model for the mTORC1 signaling pathway achieved the best performance with the mean absolute errors (MAEs) of 27.29 and 28.61% in internal and external test sets, respectively.

CONCLUSIONS

The DCE-MRI features were associated with hallmark activities and may improve individualized medicine for BC by noninvasively predicting pathway-level heterogeneity.

Resistance to CDK4/6 inhibition: Mechanisms and strategies to overcome a therapeutic problem in the treatment of hormone receptor-positive metastatic breast cancer

Selective CDK4/6 inhibitors, such as palbociclib, ribociclib, and abemaciclib, have been approved in combination with hormone therapy for the treatment of patients with HR+, HER2-negative advanced or metastatic breast cancer (mBC). Despite their promising activity, approximately 10 % of patients have de novo resistance, while the rest of them will develop acquired resistance after 24-28 months when used as first-line therapy and after a shorter period when used as second-line therapy. Various mechanisms of resistance to CDK4/6 inhibitors have been described, including cell cycle-related mechanisms, such as RB loss, p16 amplification, CDK6 or CDK4 amplification, and cyclin E-CDK2 amplification. Other bypass mechanisms involve the activation of FGFR or PI3K/AKT/mTOR pathways. Identifying the different mechanisms by which resistance to CDK4/6 inhibitors occurs may help to design new treatment strategies to improve patient outcomes. This review presents the currently available knowledge on the mechanisms of resistance to CDK4/6 inhibitors, explores possible treatment strategies that could overcome this therapeutic problem, and summarizes relevant recent clinical trials.

E2F1-initiated transcription of PRSS22 promotes breast cancer metastasis by cleaving ANXA1 and activating FPR2/ERK signaling pathway

Breast cancer (BC) is the most common malignant tumor in women worldwide. Metastasis is the main cause of BC-related death. The specific mechanism underlying BC metastasis remains obscure. Recently, PRSS22 was discovered to be involved in tumor development, however, its detailed biological function and regulatory mechanism in BC are unclear. Here, we characterized that PRSS22 expression is upregulated in BC tissues compared with non-tumorous breast tissues. Dual luciferase assays, bioinformatics analyses and chromatin immunoprecipitation (ChIP) assays indicated that transcription factor E2F1 directly binds to the PRSS22 promoter region and activates its transcription. Functionally, upregulation of PRSS22 promoted invasion and metastasis of BC cells in vitro and in vivo, whereas knockdown of PRSS22 inhibited its function. Mechanistically, the combination of PRSS22 and ANXA1 protein in BC cells was first screened by protein mass spectrometry analysis, and then confirmed by co-immunoprecipitation (Co-IP) and western blot assays. Co-overexpression of PRSS22 and ANXA1 could promote BC cell migration and invasion. We further demonstrated that PRSS22 promotes the cleavage of ANXA1 and in turn generates an N-terminal peptide, which initiates the FPR2/ERK signaling axis to increase BC aggressiveness.

Integration of differential expression and network structure for 'omics data analysis

Differential expression (DE) analysis has been routinely used to identify molecular features that are statistically significantly different between distinct biological groups. In recent years, differential network (DN) analysis has emerged as a powerful approach to uncover molecular network structure changes from one biological condition to the other where the molecular features with larger topological changes are selected as biomarkers. Although a large number of DE and a few DN-based methods are available, they have been usually implemented independently. DE analysis ignores the relationship among molecular features while DN analysis does not account for the expression changes at individual level. Therefore, an integrative analysis approach that accounts for both DE and DN is required to identify disease associated key features. Although, a handful of methods have been proposed, there is no method that optimizes the combination of DE and DN. We propose a novel integrative analysis method, DNrank, to identify disease-associated molecular features that leverages the strengths of both DE and DN by calculating a weight using resampling based cross validation scheme within the algorithm. First, differential expression analysis of individual molecular features is carried out. Second, a differential network structure is constructed using the differential partial correlation analysis. Third, the molecular features are ranked in the order of their significances by integrating their DE measures and DN structure using the modified Google's PageRank algorithm. In the algorithm, the optimum combination of DE and DN analyses is achieved by evaluating the prediction performance of top-ranked features utilizing support vector machine classifier with Monte Carlo cross validation. The proposed method is illustrated using both simulated data and three real data sets. The results show that the proposed method has a better performance in identifying important molecular features with respect to predictive discrimination. Also, as compared to existing feature selection methods, the top-ranked features selected by our method had a higher stability in selection. DNrank allows the researchers to identify the disease-associated features by utilizing both expression and network topology changes between two groups.

Troponin C-1 Activated by E2F1 Accelerates Gastric Cancer Progression via Regulating TGF-β/Smad Signaling

BACKGROUND

Troponin C-1 (TNNC1) has been previously characterized as an oncogenic gene.

AIMS

This study aimed to reveal the roles of TNNC1 in gastric cancer and the potential underlying mechanisms.

METHODS

TNNC1 siRNAs and TNNC1 overexpression plasmid were used to alter its expression in AGS, MKN45, and HGC-27 cells. CCK-8 assay, colony formation, EdU assay, flow cytometry, transwell assay, and scratch test were conducted to measure the phenotype changes. In vivo effects of TNNC1 silence were confirmed by using a xenograft mouse model. Bioinformatics analysis was conducted to screen out the transcription factor and downstream signaling of TNNC1.

RESULTS

TNNC1 was highly expressed in gastric cancer tissues and cell lines, and its expression was associated with poor prognosis. TNNC1 silence suppressed the proliferation, migration, and invasion of AGS and MKN45 cells. However, TNNC1 silence induced apoptosis by mediating the cleavage of caspase-3 and caspase-9. Overexpression of TNNC1 in HGC-27 cells led to the contrary effects. The anti-tumor effects of TNNC1 silence were also confirmed in a xenograft animal model. E2F1 was validated as an upstream transcription factor of TNNC1. Effects of TNNC1 silence on AGS cell migration and invasion were attenuated by E2F1 overexpression. Besides, TGF-β/Smad was a downstream signaling pathway of TNNC1. The anti-tumor impacts of TNNC1 silence were weaken by SB431542 (a specific inhibitor of TGF-β signaling) while accelerated by TGF-β.

CONCLUSION

TNNC1 activated by E2F1 functioned as an oncogenic gene through regulating TGF-β/Smad signaling. TNNC1 was suggested as a potential molecular drug target of gastric cancer.

BluePrint breast cancer molecular subtyping recognizes single and dual subtype tumors with implications for therapeutic guidance

Purpose

BluePrint (BP) is an 80-gene molecular subtyping test that classifies early-stage breast cancer (EBC) into Basal, Luminal, and HER2 subtypes. In most cases, breast tumors have one dominant subtype, representative of a single activated pathway. However, some tumors show a statistically equal representation of more than one subtype, referred to as dual subtype. This study aims to identify and examine dual subtype tumors by BP to understand their biology and possible implications for treatment guidance.

Methods

The BP scores of over 15,000 tumor samples from EBC patients were analyzed, and the differences between the highest and the lowest scoring subtypes were calculated. Based upon the distribution of the differences between BP scores, a threshold was determined for each subtype to identify dual versus single subtypes.

Results

Approximately 97% of samples had one single activated BluePrint molecular subtype, whereas ~ 3% of samples were classified as BP dual subtype. The most frequently occurring dual subtypes were the Luminal-Basal-type and Luminal-HER2-type. Luminal-Basal-type displays a distinct biology from the Luminal single type and Basal single type. Burstein’s classification of the single and dual Basal samples showed that the Luminal-Basal-type is mostly classified as ‘luminal androgen receptor’ and ‘mesenchymal’ subtypes, supporting molecular evidence of AR activation in the Luminal-Basal-type tumors. Tumors classified as Luminal-HER2-type resemble features of both Luminal-single-type and HER2-single-type. However, patients with dual Luminal-HER2-type have a lower pathological complete response after receiving HER2-targeted therapies in addition to chemotherapy in comparison with patients with a HER2-single-type.

Conclusion

This study demonstrates that BP identifies tumors with two active functional pathways (dual subtype) with specific transcriptional characteristics and highlights the added value of distinguishing BP dual from single subtypes as evidenced by distinct treatment response rates.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10549-022-06698-x.

Co-inhibition of ATM and ROCK synergistically improves cell proliferation in replicative senescence by activating FOXM1 and E2F1

The multifaceted nature of senescent cell cycle arrest necessitates the targeting of multiple factors arresting or promoting the cell cycle. We report that co-inhibition of ATM and ROCK by KU-60019 and Y-27632, respectively, synergistically increases the proliferation of human diploid fibroblasts undergoing replicative senescence through activation of the transcription factors E2F1 and FOXM1. Time-course transcriptome analysis identified FOXM1 and E2F1 as crucial factors promoting proliferation. Co-inhibition of the kinases ATM and ROCK first promotes the G2/M transition via FOXM1 activation, leading to accumulation of cells undergoing the G1/S transition via E2F1 activation. The combination of both inhibitors increased this effect more significantly than either inhibitor alone, suggesting synergism. Our results demonstrate a FOXM1- and E2F1-mediated molecular pathway enhancing cell cycle progression in cells with proliferative potential under replicative senescence conditions, and treatment with the inhibitors can be tested for senomorphic effect in vivo.

Inhibitors to ATM and ROCK kinases increases the proliferation via transcriptional changes in human diploid fibroblasts undergoing replicative senescence in a synergistic manner

Role of germline variants in the metastasis of breast carcinomas

Most cancer-related deaths in breast cancer patients are associated with metastasis, a multistep, intricate process that requires the cooperation of tumour cells, tumour microenvironment and metastasis target tissues. It is accepted that metastasis does not depend on the tumour characteristics but the host’s genetic makeup. However, there has been limited success in determining the germline genetic variants that influence metastasis development, mainly because of the limitations of traditional genome-wide association studies to detect the relevant genetic polymorphisms underlying complex phenotypes. In this work, we leveraged the extreme discordant phenotypes approach and the epistasis networks to analyse the genotypes of 97 breast cancer patients. We found that the host’s genetic makeup facilitates metastases by the dysregulation of gene expression that can promote the dispersion of metastatic seeds and help establish the metastatic niche—providing a congenial soil for the metastatic seeds.

Neurospora crassa is a potential source of anti-cancer agents against breast cancer

Fungi are an excellent source of pharmaceuticals including anti-tumor agents. Neurospora crassa generates metabolites with diverse structural classes, however, its potential as an anti-tumor agent source has not been explored. The purpose of this study aimed to investigate the potential of Neurospora crassa mixture against breast cancer. The in vitro T-47D and MDA-MB-231 experiments showed that N. crassa mixture at the concentrations of both 1.7 and 0.85 µg/ml significantly inhibited tumor cell proliferation, migration and invasion, and 3D spheroid formation. However, the inhibition rates of MCF-10A ranged 10-20% at concentrations of 0.85 and 1.7 µg/ml. The mixture at the concentration of 0.85 µg/ml could significantly downregulate the expressions of transcription factors of E2F1 and E2F3, cancer stem cell-related genes of LIN28, HIWI, and CD133, and onco-lncRNA HOTAIR, and increase CASP3 activity in either T-47D or MDA-MD-231 breast cancer cell lines. In vivo breast cancer C3H mouse model results showed that N. crassa mixture significantly inhibited tumor growth. These findings suggest that N. crassa contains an antitumor component(s) against breast cancer invasiveness, which may inhibit the self-renewal and differentiation of breast cancer stem cells possibly by downregulating cancer stem cell-associated and/or transcription factor genes and oncogenes, and promoting apoptosis.

ENDORSE: a prognostic model for endocrine therapy in estrogen-receptor-positive breast cancers

Advanced and metastatic estrogen receptor-positive (ER+ ) breast cancers are often endocrine resistant. However, endocrine therapy remains the primary treatment for all advanced ER+ breast cancers. Treatment options that may benefit resistant cancers, such as add-on drugs that target resistance pathways or switching to chemotherapy, are only available after progression on endocrine therapy. Here we developed an endocrine therapy prognostic model for early and advanced ER+ breast cancers. The endocrine resistance (ENDORSE) model is composed of two components, each based on the empirical cumulative distribution function of ranked expression of gene signatures. These signatures include a feature set associated with long-term survival outcomes on endocrine therapy selected using lasso-regularized Cox regression and a pathway-based curated set of genes expressed in response to estrogen. We extensively validated ENDORSE in multiple ER+ clinical trial datasets and demonstrated superior and consistent performance of the model over clinical covariates, proliferation markers, and multiple published signatures. Finally, genomic and pathway analyses in patient data revealed possible mechanisms that may help develop rational stratification strategies for endocrine-resistant ER+ breast cancer patients.

Protrudin regulates FAK activation, endothelial cell migration and angiogenesis

During angiogenesis, endothelial cells form protrusive sprouts and migrate towards the angiogenic stimulus. In this study, we investigate the role of the endoplasmic reticulum (ER)-anchored protein, Protrudin, in endothelial cell protrusion, migration and angiogenesis. Our results demonstrate that Protrudin regulates angiogenic tube formation in primary endothelial cells, Human umbilical vein endothelial cells (HUVECs). Analysis of RNA sequencing data and its experimental validation revealed cell migration as a prominent cellular function affected in HUVECs subjected to Protrudin knockdown. Further, our results demonstrate that knockdown of Protrudin inhibits focal adhesion kinase (FAK) activation in HUVECs and human aortic endothelial cells (HAECs). This is associated with a loss of polarized phospho-FAK distribution upon Protrudin knockdown as compared to Protrudin expressing HUVECs. Reduction of Protrudin also results in a perinuclear accumulation of mTOR and a decrease in VEGF-mediated S6K activation. However, further experiments suggest that the observed inhibition of angiogenesis in Protrudin knockdown cells is not affected by mTOR disturbance. Therefore, our findings suggest that defects in FAK activation and its abnormal subcellular distribution upon Protrudin knockdown are associated with a detrimental effect on endothelial cell migration and angiogenesis. Furthermore, mice with global Protrudin deletion demonstrate reduced retinal vascular progression. To conclude, our results provide evidence for a novel key role of Protrudin in endothelial cell migration and angiogenesis.

Analysis of Breast Cancer Based on the Dysregulated Network

Breast cancer is a heterogeneous disease, and its development is closely associated with the underlying molecular regulatory network. In this paper, we propose a new way to measure the regulation strength between genes based on their expression values, and construct the dysregulated networks (DNs) for the four subtypes of breast cancer. Our results show that the key dysregulated networks (KDNs) are significantly enriched in critical breast cancer-related pathways and driver genes; closely related to drug targets; and have significant differences in survival analysis. Moreover, the key dysregulated genes could serve as potential driver genes, drug targets, and prognostic markers for each breast cancer subtype. Therefore, the KDN is expected to be an effective and novel way to understand the mechanisms of breast cancer.

Fibroblast growth factor receptor signalling dysregulation and targeting in breast cancer

Fibroblast Growth Factor Receptor (FGFR) signalling plays a critical role in breast embryonal development, tissue homeostasis, tumorigenesis and metastasis. FGFR, its numerous FGF ligands and signalling partners are often dysregulated in breast cancer progression and are one of the causes of resistance to treatment in breast cancer. Furthermore, FGFR signalling on epithelial cells is affected by signals from the breast microenvironment, therefore increasing the possibility of breast developmental abnormalities or cancer progression. Increasing our understanding of the multi-layered roles of the complex family of FGFRs, their ligands FGFs and their regulatory partners may offer novel treatment strategies for breast cancer patients, as a single agent or rational co-target, which will be explored in depth in this review.

Potential roles of hsa_circ_000839 and hsa_circ_0005986 in breast cancer

Background

Breast cancer (BC) is one of the leading causes of death among women around the world. Circular RNAs (circRNAs) are a newly discovered group of non‐coding RNAs that their roles are being investigated in BC and other cancer types. In this study, we evaluated the association of hsa_circ_0005986 and hsa_circ_000839 in tumor and adjacent normal tissues of BC patients with their clinicopathological characteristics.

Materials and methods

Total RNA was extracted from tumors and adjacent non‐tumor tissues by the Trizol isolation reagent, and cDNA was synthesized using First Strand cDNA Synthesis Kit (Thermo Scientific). The expression level of hsa_circ_0005986 and hsa_circ_000839 was quantified using RT‐qPCR. Online in silico tools were used for identifying potentially important competing endogenous RNA (ceRNA) networks of these two circRNAs.

Results

The expression level of hsa_circ_0005986 and hsa_circ_000839 was lower in the tumor as compared to adjacent tissues. The expression level of hsa_circ_0005986 in the patients who had used hair dye in the last 5 years was significantly lower. Moreover, a statistically significant negative correlation between body mass index (BMI) and hsa_circ_000839 expression was observed. In silico analysis of the ceRNA network of these circRNAs revealed mRNAs and miRNAs with crucial roles in BC.

Conclusion

Downregulation of hsa_circ_000839 and hsa_circ_0005986 in BC tumors suggests a tumor‐suppressive role for these circRNAs in BC, meriting the need for more experimentations to delineate the exact mechanism of their involvement in BC pathogenesis.

In-Silico Analysis of Differentially Expressed Genes and Their Regulating microRNA Involved in Lymph Node Metastasis in Invasive Breast Carcinoma

Axillary nodal metastasis is related to poor prognosis in breast cancer (BC). Key candidate genes in BC lymph node metastasis have been identified from Gene Expression Omnibus datasets and explored through functional enrichment database for annotation, visualization and integrated discovery (DAVID) , protein-protein interaction by Search Tool for the Retrieval of Interacting Genes and proteins (STRING), network visualization (Cytoscape), survival analysis (GEPIA, KM Plotter), and target prediction (miRNet). A total of 102 overlapping differentially expressed genes were found. In-silico survival and expression analyses revealed six candidate hub genes, Desmocollin 3 (DSC3), KRT5, KRT6B, KRT17, KRT81, and SERPINB5, to be significantly associated with nodal metastasis and overall survival, and 83 MicroRNA (miRNAs), which may be potential diagnostic markers and therapeutic targets in BC patients.

E2F1 transcriptionally regulates CCNA2 expression to promote triple negative breast cancer tumorigenicity

BACKGROUND

Triple-negative breast cancer (TNBC) is a highly malignant breast cancer subtype with a poor prognosis. The cell cycle regulator cyclin A2 (CCNA2) plays a role in tumor development. Herein, we explored the role of CCNA2 in TNBC.

METHODS

We analyzed CCNA2 expression in 15 pairs of TNBC and adjacent tissues and assessed the relationship between CCNA2 expression using the tissue microarray cohort. Furthermore, we used two TNBC cohort datasets to analyze the correlation between CCNA2 and E2F transcription factor 1 (E2F1) and a luciferase reporter to explore their association. Through rescue experiments, we analyzed the effects of E2F1 knockdown on CCNA2 expression and cellular behavior.

RESULTS

We found that CCNA2 expression in TNBC was significantly higher than that in adjacent tissues with similar observations in MDA-MB-231 and MDA-MB-468 cells. E2F1 was highly correlated with CCNA2 as observed through bioinformatics analysis (R= 0.80, P< 0.001) and through TNBC tissue verification analysis (R= 0.53, P< 0.001). We determined that E2F1 binds the +677 position within the CCNA2 promoter. Moreover, CCNA2 overexpression increased cell proliferation, invasion, and migration owing to E2F1 upregulation in TNBC.

CONCLUSION

Our data indicate that E2F1 promotes TNBC proliferation and invasion by upregulating CCNA2 expression. E2F1 and CCNA2 are potential candidates that may be targeted for effective TNBC treatment.

Exosomal miRNA-205 promotes breast cancer chemoresistance and tumorigenesis through E2F1

Breast cancer (BC) is a common malignant tumor in females. The challenge in treating BC is overcoming chemoresistance. Exosome-mediated transfer of miRNAs is a molecule-shuttle in intercellular communication. Thus, we aimed to investigate whether exosomal miRNA-205 could affect chemoresistance and tumorigenesis in recipient tumor cells and to elucidate the underlying mechanism in vivo and in vitro. Microarray and qRT-PCR assays demonstrated that miRNA-205 was upregulated in tamoxifen resistance MCF-7/TAMR-1 (M/T) cells and M/T cell-derived exosomes (M/T-Exo). The M/T-Exo was internalized by human BC cells (BCCs), causing increased expression of miRNA-205 in BCCs. Coculturing with M/T-Exo promoted tamoxifen resistance, proliferation, migration, and invasion while suppressed apoptosis in recipient BCCs, which were associated with activating the caspase pathway and phosphorylating Akt. Luciferase reporter assays showed that miRNA-205 directly targeted E2F Transcription Factor 1 (E2F1) in BCCs. Furthermore, knockdown of miRNA-205 or overexpression of E2F1 reversed the roles of M/T-Exo in BCCs. In vivo experiments showed that the intratumoral injection of M/T-Exo caused greater tamoxifen resistance and larger tumor size relative to mice treated with miRNA-205-knockdown or E2F1-overexpressing BCCs. Together, the results suggest that exosomal miRNA-205 may promote tamoxifen resistance and tumorigenesis in BC through targeting E2F1 in vivo and in vitro.

Comprehensive Analysis of Peritoneal Metastasis Sequencing Data to Identify LINC00924 as a Prognostic Biomarker in Gastric Cancer

Background

Gastric cancer peritoneal metastasis has high mortality. At present, there is no effective way to cure the patients diagnosed with gastric cancer peritoneal metastasis due to its indistinct molecular mechanism. Therefore, to understand the pathogenesis and help for further target therapy, we conduct comprehensive analysis of peritoneal metastasis by bioinformatics in gastric cancer.

Methods

Microarray sequencing was used to find differential mRNAs and long non-coding RNAs (lncRNAs) expression between primary foci and peritoneal metastases lesion in gastric cancer. RT-qPCR was used to verify the expression levels of lncRNAs in gastric cancer cells after co-culture with adipocytes. TCGA, Cytoscape, lnCAR, cBioPoratal and R packages (ggrisk, survival, survminer, timeROC, forestplot, immunedeconv, ggplot2, pheatmap and ggpubr) were applied in this work.

Results

Adipocytes co-culture model was used to mimic the peritoneal microenvironment and found that LINC01151 (NR_126348), FAM27B (NR_027422) and LINC00924 (NR_027133) were up-regulated in co-culture group. Increased LINC00924 expression was significantly associated with reduced overall survival and up-regulated percentage abundance of tumor-infiltrating CD8⁺ T, B, macrophage and NK immune cells; moreover, immune checkpoint blockers (ICBs) had a worse effect on the LINC00924 high expression group. Furthermore, through univariate and multivariate Cox regression analysis, we found that LINC00924-related PEX5L in CNC network was an independent prognostic factor in gastric cancer progression.

Conclusion

LINC00924 expression was associated with poor survival, immune infiltrations and worse response to ICBs. LINC00924 might be immunotherapeutic targets of advanced gastric cancer.

DSN1 is a prognostic biomarker and correlated with clinical characterize in breast cancer

DSN1 affects cell cycle progression and is associated with clinical-pathological features in colorectal and hepatocellular carcinomas. However, the biological function of DSN1 in breast cancer is still indistinct. In this study, we comprehensively analyzed the correlation between DSN1 expression in different molecular subtypes or stages of breast cancer and investigated the prognostic value of DSN1 in databases such as Oncomine, Cancer Cell Line Encyclopedia, UALCAN, Human Protein Atlas, Kaplan-Meier Plotter, OncoLnc, GEPIA. Moreover, we investigated the correlation of DSN1 with tumor-infiltrating immune cells in the different tumor microenvironments via Tumor Immune Estimation Resource database and explore DSN1 co-expression networks in breast cancer via LinkedOmics analysis, NetworkAnalyst database analysis. Finally, we also performed our immunohistochemical experiments to explore the expression of DSN1 in different stages or subtypes of breast cancer. The findings in this article shed light on the essential role of DSN1 in breast cancers as well as suggested that the upregulation of DSN1 expression was strongly associated with poor prognosis and decreased survival in breast cancer, and there were significant differences in its expression in different pathological subtypes and stages of breast cancer.

Gypenoside LI arrests the cell cycle of breast cancer in G0/G1 phase by down-regulating E2F1

ETHNOPHARMACOLOGICAL RELEVANCE

Gynostemma pentaphyllum (Thunb.) Makino, a traditional medicine in China, has been widely used for the treatment of various diseases. Gypenoside LI (Gyp LI) is a major constituent from steamed G. pentaphyllum. Previous studies have shown that gypnenoside LI possess inhibitory effect on the growth of many cancer cells. However, its pharmacological effect in breast cancer and the mechanism have not been reported yet.

AIM OF THE STUDY

To investigate the anti-breast cancer activity of gypenoside LI and underlying mechanisms of gypenoside LI in MDA-MB-231 and MCF-7 cells.

MATERIAL/METHODS

The cytotoxicity of gypenoside LI was determined by MTT, colony-formation and three-dimensional spheroid assay. The migration, cell apoptosis and the cell cycle were investigated through cell morphology observation, flow cytometry analysis and key proteins detection. The anticancer mechanisms of gypenoside LI were detected by RNA sequencing (RNA-seq) and Gene Set Enrichment Analysis (GSEA) transcriptome analysis.

RESULTS

Gypenoside LI inhibited cell proliferation, migration, induced cell apoptosis and cell cycle arrest. Gypenoside LI arrested cell cycle at G0/G1 phase by regulating E2F1. It also inhibited tumor proliferation by regulating the expression of ERCC6L. Interestingly, we found that E2F1 siRNA also down-regulated the expression of ERCC6L. Gypenoside LI showed potential anti-breast cancer cells activity, especially on triple-negative breast cancer cells.

CONCLUSIONS

These data indicate that gypenoside LI could inhibit human breast cancer cells through inhibiting proliferation and migration, inducing apoptosis, arresting cell cycle at G0/G1 phase by regulating E2F1. It could be used as potential multi-target chemopreventive agents for cancer.

WDNE: an integrative graphical model for inferring differential networks from multi-platform gene expression data with missing values

The mechanisms controlling biological process, such as the development of disease or cell differentiation, can be investigated by examining changes in the networks of gene dependencies between states in the process. High-throughput experimental methods, like microarray and RNA sequencing, have been widely used to gather gene expression data, which paves the way to infer gene dependencies based on computational methods. However, most differential network analysis methods are designed to deal with fully observed data, but missing values, such as the dropout events in single-cell RNA-sequencing data, are frequent. New methods are needed to take account of these missing values. Moreover, since the changes of gene dependencies may be driven by certain perturbed genes, considering the changes in gene expression levels may promote the identification of gene network rewiring. In this study, a novel weighted differential network estimation (WDNE) model is proposed to handle multi-platform gene expression data with missing values and take account of changes in gene expression levels. Simulation studies demonstrate that WDNE outperforms state-of-the-art differential network estimation methods. When applied WDNE to infer differential gene networks associated with drug resistance in ovarian tumors, cell differentiation and breast tumor heterogeneity, the hub genes in the estimated differential gene networks can provide important insights into the underlying mechanisms. Furthermore, a Matlab toolbox, differential network analysis toolbox, was developed to implement the WDNE model and visualize the estimated differential networks.

miRGTF-net: Integrative miRNA-gene-TF network analysis reveals key drivers of breast cancer recurrence

Analysis of regulatory networks is a powerful framework for identification and quantification of intracellular interactions. We introduce miRGTF-net, a novel tool for construction of miRNA-gene-TF networks. We consider multiple transcriptional and post-transcriptional interaction types, including regulation of gene and miRNA expression by transcription factors, gene silencing by miRNAs, and co-expression of host genes with their intronic miRNAs. The underlying algorithm uses information on experimentally validated interactions as well as integrative miRNA/mRNA expression profiles in a given set of samples. The latter ensures simultaneous tissue-specificity and biological validity of interactions. We applied miRGTF-net to paired miRNA/mRNA-sequencing data of breast cancer samples from The Cancer Genome Atlas (TCGA). Together with topological analysis of the constructed network we showed that considered players can form reliable prognostic gene signatures for ER-positive breast cancer. A number of signatures demonstrated remarkably high accuracy on transcriptomic data obtained by both microarrays and RNA sequencing from several independent patient cohorts. Furthermore, an essential part of prognostic genes were identified as direct targets of transcription factor E2F1. The putative interplay between estrogen receptor alpha and E2F1 was suggested as a potential recurrence factor in patients treated with tamoxifen. Source codes of miRGTF-net are available at GitHub (https://github.com/s-a-nersisyan/miRGTF-net).

Human-derived osteoblast-like cells and pericyte-like cells induce distinct metastatic phenotypes in primary breast cancer cells

Approximately 70% of advanced breast cancer patients will develop bone metastases, which accounts for ∼90% of cancer-related mortality. Breast cancer circulating tumor cells (CTCs) establish metastatic tumors in the bone after a close interaction with local bone marrow cells including pericytes and osteoblasts, both related to resident mesenchymal stem/stromal cells (BM-MSCs) progenitors. In vitro recapitulation of the critical cellular players of the bone microenvironment and infiltrating CTCs could provide new insights into their cross-talk during the metastatic cascade, helping in the development of novel therapeutic strategies. Human BM-MSCs were isolated and fractionated according to CD146 presence. CD146+ cells were utilized as pericyte-like cells (PLCs) given the high expression of the marker in perivascular cells, while CD146- cells were induced into an osteogenic phenotype generating osteoblast-like cells (OLCs). Transwell migration assays were performed to establish whether primary breast cancer cells (3384T) were attracted to OLC. Furthermore, proliferation of 3384T breast cancer cells was assessed in the presence of PLC- and OLC-derived conditioned media. Additionally, conditioned media cultures as well as transwell co-cultures of each OLCs and PLCs were performed with 3384T breast cancer cells for gene expression interrogation assessing their induced transcriptional changes with an emphasis on metastatic potential. PLC as well as their conditioned media increased motility and invasion potential of 3384T breast cancer cells, while OLC induced a dormant phenotype, downregulating invasiveness markers related with migration and proliferation. Altogether, these results indicate that PLC distinctively drive 3384T cancer cells to an invasive and migratory phenotype, while OLC induce a quiescence state, thus recapitulating the different phases of the in vivo bone metastatic process. These data show that phenotypic responses from metastasizing cancer cells are influenced by neighboring cells at the bone metastatic niche during the establishment of secondary metastatic tumors.

Triple negative aggressive phenotype controlled by miR-135b and miR-365: new theranostics candidates

Triple negative breast cancer (TNBC) accounts for about a fifth of all breast cancers and includes a diverse group of cancers. The heterogeneity of TNBC and the lack of target receptors on the cell surface make it difficult to develop specific therapeutic treatments. These aspects cause the high negative prognosis of patients with this type of tumor. The analysis of the molecular profiles of TNBC samples has allowed a better characterization of this tumor, supporting the search for new reliable diagnostic markers. To this end, we have developed a bioinformatic approach to integrate networks of genes differentially expressed in basal breast cancer compared to healthy tissues, with miRNAs able to regulate their expression. We studied the role of these miRNAs in TNBC subtype cell lines. We therefore identified two miRNAs, namely miR-135b and miR-365, with a central role in regulating the altered functional pathways in basal breast cancer. These two miRNAs are differentially expressed in human TNBC immunohistochemistry-selected tissues, and their modulation has been shown to play a role in the proliferation of tumor control and its migratory and invasive capacity in TNBC subtype cell lines. From the perspective of personalized medicine, we managed to modulate the expression of the two miRNAs in organotypic cultures, suggesting their possible use as diagnostic and therapeutic molecules. miR-135b and miR-365 have a key role in TNBC, controlling proliferation and invasion. Their detection could be helpful in TNBC diagnosis, while their modulation could become a new therapeutic tool for TNBC.

The molecular mechanisms of the long noncoding RNA SBF2-AS1 in regulating the proliferation of oesophageal squamous cell carcinoma

The long noncoding RNASBF2-AS1 can promote the occurrence and development of many kinds of tumours, but its role in oesophageal squamous cell carcinoma (ESCC) is unknown. We found that SBF2-AS1 was up-regulated in ESCC, and its expression was positively correlated with tumor size (P = 0.0001), but was not related to gender, age, TNM stage, histological grade, and lymphnode metastasis (P > 0.05). It was further found that the higher the expression of SBF2-AS1, the lower the survival rate. COX multivariate analysis showed that the expression of SBF2-AS1 was an independent prognostic factor. Functional experiments show that inhibition of SBF2-AS1 can inhibit the proliferation of ESCC through in vivo and in vitro, and overexpression of SBF2-AS1 can promote the proliferation of ESCC and inhibit its apoptosis. In mechanism, SBF2-AS1/miR-338-3P, miR-362-3P/E2F1 axis are involved in the regulation of ESCC growth. In general, SBF2-AS1 may be used as ceRNA to combine with miR-338-3P and miR-362-3P to up-regulate the expression ofE2F1, and ultimately play a role in promoting cancer. It may be used as a therapeutic target and a biomarker for prognosis.

Gene expression analysis of human prostate cell lines with and without tumor metastasis suppressor CD82

BACKGROUND

Tetraspanin CD82 is a tumor metastasis suppressor that is known to down regulate in various metastatic cancers. However, the exact mechanism by which CD82 prevents cancer metastasis is unclear. This study aims to identify genes that are regulated by CD82 in human prostate cell lines.

METHODS

We used whole human genome microarray to obtain gene expression profiles in a normal prostate epithelial cell line that expressed CD82 (PrEC-31) and a metastatic prostate cell line that does not express CD82 (PC3). Then, siRNA silencing was used to knock down CD82 expression in PrEC-31 while CD82 was re-expressed in PC3 to acquire differentially-expressed genes in the respective cell line.

RESULTS

Differentially-expressed genes with a P < 0.05 were identified in 3 data sets: PrEC-31 (+CD82) vs PrEC-31(-CD82), PC3-57 (+CD82) vs. PC3-5 V (-CD82), and PC3-29 (+CD82) vs. PC3-5 V (-CD82). Top 25 gene lists did not show overlap within the data sets, except (CALB1) the calcium binding protein calbindin 1 which was significantly up-regulated (2.8 log fold change) in PrEC-31 and PC3-29 cells that expressed CD82. Other most significantly up-regulated genes included serine peptidase inhibitor kazal type 1 (SPINK1) and polypeptide N-acetyl galactosaminyl transferase 14 (GALNT14) and most down-regulated genes included C-X-C motif chemokine ligand 14 (CXCL14), urotensin 2 (UTS2D), and fibroblast growth factor 13 (FGF13). Pathways related with cell proliferation and angiogenesis, migration and invasion, cell death, cell cycle, signal transduction, and metabolism were highly enriched in cells that lack CD82 expression. Expression of two mutually inclusive genes in top 100 gene lists of all data sets, runt-related transcription factor (RUNX3) and trefoil factor 3 (TFF3), could be validated with qRT-PCR.

CONCLUSION

Identification of genes and pathways regulated by CD82 in this study may provide additional insights into the role that CD82 plays in prostate tumor progression and metastasis, as well as identify potential targets for therapeutic intervention.

Quadruple-negative breast cancer: novel implications for a new disease

Based on the androgen receptor (AR) expression, triple-negative breast cancer (TNBC) can be subdivided into AR-positive TNBC and AR-negative TNBC, also known as quadruple-negative breast cancer (QNBC). QNBC characterization and treatment is fraught with many challenges. In QNBC, there is a greater paucity of prognostic biomarkers and therapeutic targets than AR-positive TNBC. Although the prognostic role of AR in TNBC remains controversial, many studies revealed that a lack of AR expression confers a more aggressive disease course. Literature characterizing QNBC tumor biology and uncovering novel biomarkers for improved management of the disease remains scarce. In this comprehensive review, we summarize the current QNBC landscape and propose avenues for future research, suggesting potential biomarkers and therapeutic strategies that warrant investigation.

Low E2F2 activity is associated with high genomic instability and PARPi resistance

The E2F family, classically known for a central role in cell cycle, has a number of emerging roles in cancer including angiogenesis, metabolic reprogramming, metastasis and DNA repair. E2F1 specifically has been shown to be a critical mediator of DNA repair; however, little is known about DNA repair and other E2F family members. Here we present an integrative bioinformatic and high throughput drug screening study to define the role of E2F2 in maintaining genomic integrity in breast cancer. We utilized in vitro E2F2 ChIP-chip and over expression data to identify transcriptional targets of E2F2. This data was integrated with gene expression from E2F2 knockout tumors in an MMTV-Neu background. Finally, this data was compared to human datasets to identify conserved roles of E2F2 in human breast cancer through the TCGA breast cancer, Cancer Cell Line Encyclopedia, and CancerRx datasets. Through these methods we predict that E2F2 transcriptionally regulates mediators of DNA repair. Our gene expression data supports this hypothesis and low E2F2 activity is associated with a highly unstable tumor. In human breast cancer E2F2, status was also correlated with a patient's response to PARP inhibition therapy. Taken together this manuscript defines a novel role of E2F2 in cancer progression beyond cell cycle and could impact patient treatment.

Cross-talk between myeloid-derived suppressor cells and Mucin1 in breast cancer vaccination: On the verge of a breakthrough

Although breast cancer is one of the leading troublesome cancers, the available therapeutic options have not fulfilled the desired outcomes. Immune-based therapy has gained special attention for breast cancer treatment. Although this approach is highly tolerable, its low response rate has rendered it as an undesirable approach. This review aims to describe the essential oncogenic pathways involved in breast cancer, elucidate the immunosuppression and oncogenic effect of Mucin1, and introduce myeloid-derived suppressor cells, which are the main culprits of anti-tumoral immune response attenuation. The various auto-inductive loops between Mucin1 and myeloid-derived suppressor cells are focal in the suppression of anti-tumoral immune responses in patients with breast cancer. These cross-talks between the Mucin1 and myeloid-derived suppressor cells can be the underlying causes of immunotherapy's impotence for patients with breast cancer. This approach can pave the road for the development of a potent vaccine for patients with breast cancer and is translated into clinical settings.

Analysis of Copy-Number Variations and Feline Mammary Carcinoma Survival

Feline mammary carcinomas (FMCs) are highly malignant. As the disease-free survival (DFS) and overall survival (OS) are short, prognostication is crucial. Copy-number variations (CNVs) analysis by next-generation sequencing serves to identify critical cancer-related genomic regions. Thirty-three female cats with FMCs were followed during two years after surgery. Tumours represented tubulopapillary and solid carcinomas encompassing six molecular subtypes. Regardless of the histopathological diagnosis, molecular subtypes showed important differences in survival. Luminal A tumours exhibited the highest DFS (p = 0.002) and cancer-specific OS (p = 0.001), and the lowest amount of CNVs (p = 0.0001). In contrast, basal-like triple-negative FMCs had the worst outcome (DFS, p < 0.0001; and OS, p < 0.00001) and were the most aberrant (p = 0.05). In the multivariate analysis, copy-number losses (CNLs) in chromosome B1 (1-23 Mb) harbouring several tumour-repressors (e.g. CSMD1, MTUS1, MSR1, DBC2, and TUSC3) negatively influenced DFS. Whereas, copy-number gains (CNGs) in B4 (1-29 Mb) and F2 (64-82.3 Mb) comprising epithelial to mesenchymal transition genes and metastasis-promoting transcription factors (e.g. GATA3, VIM, ZEB1, and MYC) negatively influenced DFS and cancer-specific OS. These data evidence an association between specific CNVs in chromosomes B1, B4 and F2, and poor prognosis in FMCs.