ETV7 regulates breast cancer stem-like cell features by repressing IFN-response genes

Anoikis-related genes linked with patient outcome in pancreatic cancer

Anoikis is programmed cell death occurring upon cell detachment from the extracellular matrix. Cancer cells need to evade anoikis to be able to metastasize to distant sites. However, the molecular features and prognostic value of anoikis-related genes (ARGs) in pancreatic cancer remain unclear. In this study, we utilized transcriptome data from the TCGA and GSE102238 databases to identify 64 ARGs significantly associated with prognosis. We used the "ConsensusClusterPlus" R package to stratify patients into high and low-risk prognostic subgroups. The KEGG and GSEA analyses revealed that the clusters with poor prognosis were enriched for the ECM receptor interaction pathway, the TP53 signaling pathway, and the galactose metabolism pathway, and that the cell cycle pathway was upregulated. A prognostic model consisting of seven ARGs (SERPINE1, EGF, E2F1, MSLN, RAB27B, ETV7, MST1) was constructed using LASSO regression and when combined with clinicopathological parameters using Cox regression, a prognostic Nomogram was created, which demonstrated high prognostic utility. Among the biomarker candidates, we report ETV7 as a novel, independent prognostic marker in pancreatic cancer. ETV7 was highly expressed in KRAS and TP53 co-occurrent mutant TCGA patients, indicating that it may be regulated by the two major driver genes of pancreatic cancer. Therefore, targeting ETV7 could be a potential focus for future therapeutic studies.

Mitigation of Breast Cancer Cells' Invasiveness via Down Regulation of ETV7, Hippo, and PI3K/mTOR Pathways by Vitamin D3 Gold-Nanoparticles

Metastasis in breast cancer is the major cause of death in females (about 30%). Based on our earlier observation that Vitamin D3 downregulates mTOR, we hypothesized that Vitamin D3 conjugated to gold nanoparticles (VD3-GNPs) reduces breast cancer aggressiveness by downregulating the key cancer controller PI3K/AKT/mTOR. Western blots, migration/invasion assays, and other cell-based, biophysical, and bioinformatics studies are used to study breast cancer cell aggressiveness and nanoparticle characterization. Our VD3-GNP treatment of breast cancer cells (MCF-7 and MDA-MB-231) significantly reduces the aggressiveness (cancer cell migration and invasion rates > 45%) via the simultaneous downregulation of ETV7 and the Hippo pathway. Consistent with our hypothesis, we, indeed, found a downregulation of the PI3K/AKT/mTOR pathway. It is surprising that the extremely low dose of VD3 in the nano formulation (three orders of magnitude lower than in earlier studies) is quite effective in the alteration of cancer invasiveness and cell signaling pathways. Clearly, VD3-GNPs are a viable candidate for non-toxic, low-cost treatment for reducing breast cancer aggressiveness.

An atlas of causal and mechanistic drivers of interpatient heterogeneity in glioma

Grade IV glioma, formerly known as glioblastoma multiforme (GBM) is the most aggressive and lethal type of brain tumor, and its treatment remains challenging in part due to extensive interpatient heterogeneity in disease driving mechanisms and lack of prognostic and predictive biomarkers. Using mechanistic inference of node-edge relationship (MINER), we have analyzed multiomics profiles from 516 patients and constructed an atlas of causal and mechanistic drivers of interpatient heterogeneity in GBM (gbmMINER). The atlas has delineated how 30 driver mutations act in a combinatorial scheme to causally influence a network of regulators (306 transcription factors and 73 miRNAs) of 179 transcriptional "programs", influencing disease progression in patients across 23 disease states. Through extensive testing on independent patient cohorts, we share evidence that a machine learning model trained on activity profiles of programs within gbmMINER significantly augments risk stratification, identifying patients who are super-responders to standard of care and those that would benefit from 2 nd line treatments. In addition to providing mechanistic hypotheses regarding disease prognosis, the activity of programs containing targets of 2 nd line treatments accurately predicted efficacy of 28 drugs in killing glioma stem-like cells from 43 patients. Our findings demonstrate that interpatient heterogeneity manifests from differential activities of transcriptional programs, providing actionable strategies for mechanistically characterizing GBM from a systems perspective and developing better prognostic and predictive biomarkers for personalized medicine.

IFI35 limits antitumor immunity in triple-negative breast cancer via CCL2 secretion

Triple-negative breast cancer (TNBC) is the most aggressive subtype of breast cancer with poor prognosis due to the lack of therapeutic targets. Although immunotherapy brings survival benefits to patients diagnosed with TNBC, it remains limited and treatment resistance is widespread. Here we demonstrate that IFI35 is highly expressed in tumor tissues and can be induced by Interferon-γ in a time-dependent and concentration-dependent manner in breast cancer cells. In xenograft models, we reveal that IFI35 dramatically increases myeloid-derived suppressor cells infiltration in tumors, along with depletion and anergy of CD8+T cells. IFI35 ablation leads to prolonged survival of the mice. Mechanistically, RNA-sequencing reveals that IFI35 promotes CCL2 secretion, resulting in the remodeling of TNBC immune microenvironment. Ablation of IFI35 promotes the infiltration of effector CD8+T cells, and thereby sensitizes TNBC to anti-PD-1 immunotherapy. Our data suggest that IFI35 limits antitumor immunity and may be expected to become a new immunotherapy target in TNBC.

Epigenetic variation impacts individual differences in the transcriptional response to influenza infection

Humans display remarkable interindividual variation in their immune response to identical challenges. Yet, our understanding of the genetic and epigenetic factors contributing to such variation remains limited. Here we performed in-depth genetic, epigenetic and transcriptional profiling on primary macrophages derived from individuals of European and African ancestry before and after infection with influenza A virus. We show that baseline epigenetic profiles are strongly predictive of the transcriptional response to influenza A virus across individuals. Quantitative trait locus (QTL) mapping revealed highly coordinated genetic effects on gene regulation, with many cis-acting genetic variants impacting concomitantly gene expression and multiple epigenetic marks. These data reveal that ancestry-associated differences in the epigenetic landscape can be genetically controlled, even more than gene expression. Lastly, among QTL variants that colocalized with immune-disease loci, only 7% were gene expression QTL, while the remaining genetic variants impact epigenetic marks, stressing the importance of considering molecular phenotypes beyond gene expression in disease-focused studies.

ETV7 promotes colorectal cancer progression through upregulation of IFIT3

Members of the E26 transformation-specific (ETS) variant transcription factor family act as either tumor suppressors or oncogenic factors in numerous types of cancer. ETS variant transcription factor 7 (ETV7) participates in the development of malignant tumors, whereas its involvement in colorectal cancer (CRC) is less clear. In this study, The Cancer Genome Atlas (TCGA) and immunochemistry staining were applied to check the clinical relevance of ETV7 and interferon-induced protein with tetratricopeptide repeats 3 (IFIT3) in CRC patients. Overexpression and knockdown of ETV7 and IFIT3 were conducted by transfecting the cells with pCDNA3.1 plasmids and siRNAs, respectively. Western blotting was used to detect the protein expression of ETV7 in CRC cells. Cell Counting Kit-8, cell colony formation, and Transwell assays, as well as flow cytometry, were used to evaluate the proliferation, migration, cell cycle, and apoptosis of CRC cells. Furthermore, western blotting, RT-qPCR, and luciferase assay were used to explore the regulation of ETV7 on IFIT3. Rescue assay was used to investigate the significance of ETV7/IFIT3 axis on CRC progression. We found that ETV7 was upregulated in CRC tissues and cells. Overexpression of ETV7 stimulated the proliferation, migration, and cell cycle amplification, and reduced the apoptosis of CRC cells. Downregulation of ETV7 exerted the opposite effect on CRC cell progression. Moreover, we demonstrated that ETV7 stimulated the transcription activity, the mRNA and protein expression of IFIT3 in CRC cells. There was a positive correlation between ETV7 and IFIT3 in CRC patients. IFIT3 knockdown reversed the promotive effect exerted by overexpression of ETV7 on the amplification and migration of CRC cells. By contrast, overexpression of IFIT3 blocked the inhibitory effect of ETV7-targeting siRNA. In summary, ETV7 induces progression of CRC by activating the transcriptional expression of IFIT3. The EVT7/IFIT3 axis may be a novel target for CRC therapy.

Effect of 2-methoxyestradiol treatment on early- and late-stage breast cancer progression in a mouse model

The prevalence of breast cancer (BC) continues to increase and is the leading cause of cancer deaths in many countries. Numerous in vitro and in vivo studies have demonstrated that 2-methoxyestradiol (2-ME) has antiproliferative and antiangiogenic effects in BC, thereby inhibiting tumour growth and metastasis. We compared the effect of 2-ME in early- and late-stage BC using a transgenic mouse model-FVB/N-Tg(MMTV-PyVT)-of spontaneously development of aggressive mammary carcinoma with lung metastasis. Mice received 100 mg/kg 2-ME treatment immediately when palpable mammary tumours were identified (early-stage BC; Experimental group 1) and 28 days after palpable mammary tumours were detected (late-stage BC; Experimental group 2). 2-ME was administered via oral gavage three times a week for 28 days after initiation of treatment, whereas control mice received the vehicle containing 10% dimethyl sulfoxide and 90% sunflower oil for the same duration as the treatment group. Mammary tumours were measured weekly over the 28 days and at termination, blood, mammary and lung tissue were collected for analysis. Mice with a tumour volume threshold of 4000 mm3 were killed before the treatment regime was completed. 2-ME treatment of early-stage BC led to lower levels of mammary tumour necrosis, whereas tumour mass and volume were increased. Additionally, necrotic lesions and anti-inflammatory CD163-expressing cells were more frequent in pulmonary metastatic tumours in this group. In contrast, 2-ME treatment of late-stage BC inhibited tumour growth over the 28-day period and resulted in increased CD3+ cell number and tumour necrosis. Furthermore, 2-ME treatment slowed down pulmonary metastasis but did not increase survival of late-stage BC mice. Besides late-stage tumour necrosis, none of the other results were statistically significant. This study demonstrates that 2-ME treatment has an antitumour effect on late-stage BC, however, with no increase in survival rate, whereas the treatment failed to demonstrate any benefit in early-stage BC.

A novel CD8+ T cell-related gene signature for predicting the prognosis and immunotherapy efficacy in bladder cancer

OBJECTIVE

To identify CD8+ T cell-related molecular clusters and establish a novel gene signature for predicting the prognosis and efficacy of immunotherapy in bladder cancer (BCa).

METHODS

Transcriptome and clinical data of BCa samples were obtained from the Cancer Genome Atlas (TCGA) and GEO databases. The CD8+ T cell-related genes were screened through the CIBERSORT algorithm and correlation analysis. Consensus clustering analysis was utilized to identified CD8+ T cell-related molecular clusters. A novel CD8+ T cell-related prognostic model was developed using univariate Cox regression analysis and Lasso regression analysis. Internal and external validations were performed and the validity of the model was validated in a real-world cohort. Finally, preliminary experimental verifications were carried out to verify the biological functions of SH2D2A in bladder cancer.

RESULTS

A total of 52 CD8+ T cell-related prognostic genes were screened and two molecular clusters with notably diverse immune cell infiltration, prognosis and clinical features were developed. Then, a novel CD8+ T cell-related prognostic model was constructed. The patients with high-risk scores exhibited a significantly worse overall survival in training, test, whole TCGA and validating cohort. The AUC was 0.766, 0.725, 0.739 and 0.658 in the four cohorts sequentially. Subgroup analysis suggested that the novel prognostic model has a robust clinical application for selecting high-risk patients. Finally, we confirmed that patients in the low-risk group might benefit more from immunotherapy or chemotherapy, and validated the prognostic model in a real-world immunotherapy cohort. Preliminary experiment showed that SH2D2A was capable of attenuating proliferation, migration and invasion of BCa cells.

CONCLUSIONS

CD8+ T cell-related molecular clusters were successfully identified. Besides, a novel CD8+ T cell-related prognostic model with an excellent predictive performance in predicting survival rates and immunotherapy efficacy of BCa was developed.

RNAseq profiling of blood from patients with coronary artery disease: Signature of a T cell imbalance

Background

Cardiovascular disease had a global prevalence of 523 million cases and 18.6 million deaths in 2019. The current standard for diagnosing coronary artery disease (CAD) is coronary angiography either by invasive catheterization (ICA) or computed tomography (CTA). Prior studies employed single-molecule, amplification-independent RNA sequencing of whole blood to identify an RNA signature in patients with angiographically confirmed CAD. The present studies employed Illumina RNAseq and network co-expression analysis to identify systematic changes underlying CAD.

Methods

Whole blood RNA was depleted of ribosomal RNA (rRNA) and analyzed by Illumina total RNA sequencing (RNAseq) to identify transcripts associated with CAD in 177 patients presenting for elective invasive coronary catheterization. The resulting transcript counts were compared between groups to identify differentially expressed genes (DEGs) and to identify patterns of changes through whole genome co-expression network analysis (WGCNA).

Results

The correlation between Illumina amplified RNAseq and the prior SeqLL unamplified RNAseq was quite strong (r = 0.87), but there was only 9 % overlap in the DEGs identified. Consistent with the prior RNAseq, the majority (93 %) of DEGs were down-regulated ~1.7-fold in patients with moderate to severe CAD (>20 % stenosis). DEGs were predominantly related to T cells, consistent with known reductions in Tregs in CAD. Network analysis did not identify pre-existing modules with a strong association with CAD, but patterns of T cell dysregulation were evident. DEGs were enriched for transcripts associated with ciliary and synaptic transcripts, consistent with changes in the immune synapse of developing T cells.

Conclusions

These studies confirm and extend a novel mRNA signature of a Treg-like defect in CAD. The pattern of changes is consistent with stress-related changes in the maturation of T and Treg cells, possibly due to changes in the immune synapse.

ETV7 reduces inflammatory responses in breast cancer cells by repressing the TNFR1/NF-κB axis

The transcription factor ETV7 is an oncoprotein that is up-regulated in all breast cancer (BC) types. We have recently demonstrated that ETV7 promoted breast cancer progression by increasing cancer cell proliferation and stemness and was also involved in the development of chemo- and radio-resistance. However, the roles of ETV7 in breast cancer inflammation have yet to be studied. Gene ontology analysis previously performed on BC cells stably over-expressing ETV7 demonstrated that ETV7 was involved in the suppression of innate immune and inflammatory responses. To better decipher the involvement of ETV7 in these signaling pathways, in this study, we identified TNFRSF1A, encoding for the main receptor of TNF-α, TNFR1, as one of the genes down-regulated by ETV7. We demonstrated that ETV7 directly binds to the intron I of this gene, and we showed that the ETV7-mediated down-regulation of TNFRSF1A reduced the activation of NF-κB signaling. Furthermore, in this study, we unveiled a potential crosstalk between ETV7 and STAT3, another master regulator of inflammation. While it is known that STAT3 directly up-regulates the expression of TNFRSF1A, here we demonstrated that ETV7 reduces the ability of STAT3 to bind to the TNFRSF1A gene via a competitive mechanism, recruiting repressive chromatin remodelers, which results in the repression of its transcription. The inverse correlation between ETV7 and TNFRSF1A was confirmed also in different cohorts of BC patients. These results suggest that ETV7 can reduce the inflammatory responses in breast cancer through the down-regulation of TNFRSF1A.

ETS transcription factors: Multifaceted players from cancer progression to tumor immunity

The E26 transformation specific (ETS) family comprises 28 transcription factors, the majority of which are involved in tumor initiation and development. Serving as a group of functionally heterogeneous gene regulators, ETS factors possess a structurally conserved DNA-binding domain. As one of the most prominent families of transcription factors that control diverse cellular functions, ETS activation is modulated by multiple intracellular signaling pathways and post-translational modifications. Disturbances in ETS activity often lead to abnormal changes in oncogenicity, including cancer cell survival, growth, proliferation, metastasis, genetic instability, cell metabolism, and tumor immunity. This review systematically addresses the basics and advances in studying ETS factors, from their tumor relevance to clinical translational utility, with a particular focus on elucidating the role of ETS family in tumor immunity, aiming to decipher the vital role and clinical potential of regulation of ETS factors in the cancer field.

Inside the stemness engine: Mechanistic links between deregulated transcription factors and stemness in cancer

Cell identity is largely determined by its transcriptional profile. In tumour, deregulation of transcription factor expression and/or activity enables cancer cell to acquire a stem-like state characterised by capacity to self-renew, differentiate and form tumours in vivo. These stem-like cancer cells are highly metastatic and therapy resistant, thus warranting a more complete understanding of the molecular mechanisms downstream of the transcription factors that mediate the establishment of stemness state. Here, we review recent research findings that provide a mechanistic link between the commonly deregulated transcription factors and stemness in cancer. In particular, we describe the role of master transcription factors (SOX, OCT4, NANOG, KLF, BRACHYURY, SALL, HOX, FOX and RUNX), signalling-regulated transcription factors (SMAD, β-catenin, YAP, TAZ, AP-1, NOTCH, STAT, GLI, ETS and NF-κB) and unclassified transcription factors (c-MYC, HIF, EMT transcription factors and P53) across diverse tumour types, thereby yielding a comprehensive overview identifying shared downstream targets, highlighting unique mechanisms and discussing complexities.

Exploration the global single-cell ecological landscape of adenomyosis-related cell clusters by single-cell RNA sequencing

Background: Adenomyosis (AM) is a common benign uterine disease that threatens the normal life of patients. Cells associated with microenvironmental immune ecology are crucial in AM, although they are not as well understood at the cellular level.

Methods: Single-cell sequencing (scRNA-seq) data were used to construct an AM global single-cell map, to further identify relevant cell clusters and infer chromosomal copy number variation (CNV) in AM samples. The biological functions of cell clusters were explored and cellular evolutionary processes were inferred by enrichment analysis and pseudotime analysis. In addition, a gene regulatory network (GRN) analysis was constructed to explore the regulatory role of transcription factors in AM progression.

Results: We obtained the expression profiles of 42260 cells and identified 10 cell clusters. By comparing the differences in cell components between AM patients and controls, we found that significant abundance of endometrial cells (EC), epithelial cells (Ep), endothelial cells (En), and smooth muscle cells (SMC) in AM patients. Cell clusters with high CNV levels possessing tumour-like features existed in the ectopic endometrium samples. Moreover, the Ep clusters were significantly involved in leukocyte transendothelial cell migration and apoptosis, suggesting an association with cell apoptosis and migration. En clusters were mainly involved in pathways in cancer and apoptosis, indicating that En has certain malignant features.

Conclusion: This study identified cell clusters with immune-related features, investigated the changes in the immune ecology of the microenvironment of these cells during AM, and provided a new strategy for the treatment of AM.

A Three-Gene Signature Predicts Response to Selinexor in Multiple Myeloma

PURPOSE

Selinexor is the first selective inhibitor of nuclear export to be approved for the treatment of relapsed or refractory multiple myeloma (MM). Currently, there are no known genomic biomarkers or assays to help select MM patients at higher likelihood of response to selinexor. Here, we aimed to characterize the transcriptomic correlates of response to selinexor-based therapy.

METHODS

We performed RNA sequencing on CD138+ cells from the bone marrow of 100 patients with MM who participated in the BOSTON study, followed by differential gene expression and pathway analysis. Using the differentially expressed genes, we used cox proportional hazard models to identify a gene signature predictive of response to selinexor, followed by validation in external cohorts.

RESULTS

The three-gene signature predicts response to selinexor-based therapy in patients with MM in the BOSTON cohort. Then, we validated this gene signature in 64 patients from the STORM cohort of triple-class refractory MM and additionally in an external cohort of 35 patients treated in a real-world setting outside of clinical trials. We found that the signature tracks with both depth and duration of response, and it also validates in a different tumor type using a cohort of pretreatment tumors from patients with recurrent glioblastoma. Furthermore, the genes involved in the signature, WNT10A, DUSP1, and ETV7, reveal a potential mechanism through upregulated interferon-mediated apoptotic signaling that may prime tumors to respond to selinexor-based therapy.

CONCLUSION

In this study, we present a present a novel, three-gene expression signature that predicts selinexor response in MM. This signature has important clinical relevance as it could identify patients with cancer who are most likely to benefit from treatment with selinexor-based therapy.

Targeting Cancer Stem Cells through Epigenetic Modulation of Interferon Response

Cancer stem cells (CSCs) are a small subset of cancer cells and are thought to play a critical role in the initiation and maintenance of tumor mass. CSCs exhibit similar hallmarks to normal stem cells, such as self-renewal, differentiation, and homeostasis. In addition, CSCs are equipped with several features so as to evade anticancer mechanisms. Therefore, it is hard to eliminate CSCs by conventional anticancer therapeutics that are effective at clearing bulk cancer cells. Interferons are innate cytokines and are the key players in immune surveillance to respond to invaded pathogens. Interferons are also crucial for adaptive immunity for the killing of specific aliens including cancer cells. However, CSCs usually evolve to escape from interferon-mediated immune surveillance and to shape the niche as a “cold” tumor microenvironment (TME). These CSC characteristics are related to their unique epigenetic regulations that are different from those of normal and bulk cancer cells. In this review, we introduce the roles of epigenetic modifiers, focusing on LSD1, BMI1, G9a, and SETDB1, in contributing to CSC characteristics and discussing the interplay between CSCs and interferon response. We also discuss the emerging strategy for eradicating CSCs by targeting these epigenetic modifiers, which can elevate cytosolic nuclei acids, trigger interferon response, and reshape a “hot” TME for improving cancer immunotherapy. The key epigenetic and immune genes involved in this crosstalk can be used as biomarkers for precision oncology.

Altered Expression of Shorter p53 Family Isoforms Can Impact Melanoma Aggressiveness

Cutaneous melanoma is the most aggressive form of skin cancer. Despite the significant advances in the management of melanoma in recent decades, it still represents a challenge for clinicians. The TP53 gene, the guardian of the genome, which is altered in more than 50% of human cancers, is rarely mutated in melanoma. More recently, researchers started to appreciate the importance of shorter p53 isoforms as potential modifiers of the p53-dependent responses. We analyzed the expression of p53 and p73 isoforms both at the RNA and protein level in a panel of melanoma-derived cell lines with different TP53 and BRAF status, in normal conditions or upon treatment with common anti-cancer DNA damaging agents or targeted therapy. Using lentiviral vectors, we also generated stable clones of H1299 p53 null cells over-expressing the less characterized isoforms Δ160p53α, Δ160p53β, and Δ160p53γ. Further, we obtained two melanoma-derived cell lines resistant to BRAF inhibitor vemurafenib. We observed that melanoma cell lines expressed a wide array of p53 and p73 isoforms, with Δ160p53α as the most variable one. We demonstrated for the first time that Δ160p53α, and to a lesser extent Δ160p53β, can be recruited on chromatin, and that Δ160p53γ can localize in perinuclear foci; moreover, all Δ160p53 isoforms can stimulate proliferation and in vitro migration. Lastly, vemurafenib-resistant melanoma cells showed an altered expression of p53 and p73 isoforms, namely an increased expression of potentially pro-oncogenic Δ40p53β and a decrease in tumor-suppressive TAp73β. We therefore propose that p53 family isoforms can play a role in melanoma cells' aggressiveness.