E2F transcription factors and digestive system malignancies: How much do we know?

Characterizing the key genes of COVID-19 that regulate tumor immune microenvironment and prognosis in hepatocellular carcinoma

Hepatocellular carcinoma (HCC), a highly heterogeneous malignant tumor associated with a poor prognosis, is a common cause of cancer-related deaths worldwide, with a limited survival benefit for patients despite ongoing therapeutic breakthroughs. Coronavirus disease 2019 (COVID-19), a severe infectious disease caused by severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2), is a global pandemic and a serious threat to human health. The increased susceptibility to SARS-CoV-2 infection and a poor prognosis in patients with cancer necessitate the exploration of the potential link between the two. No studies have investigated the relationship of COVID-19 genes with the prognosis and tumor development in patients with HCC. We screened prognosis-related COVID-19 genes in HCC, performed molecular typing, developed a stable and reliable COVID-19 genes signature for predicting survival, characterized the immune microenvironment in HCC patients, and explored new molecular therapeutic targets. Datasets of HCC patients, including RNA sequencing data and clinical information, were obtained from The Cancer Genome Atlas (TCGA), International Cancer Genome Consortium (ICGC), and Gene Expression Omnibus (GEO) databases. Prognosis-related COVID-19 genes were identified by univariate Cox analysis. Molecular typing of HCC was performed using the consensus non-negative matrix factorization method (cNMF), followed by the analysis of survival, tumor microenvironment, and pathway enrichment for each subtype. Prognostic signatures were constructed using LASSO-Cox regression models, and receiver operating characteristic (ROC) curves were used to validate the predictive performance of the signature. The same approach was used for the test and external validation sets. Seven software packages were applied to determine the abundance of immune infiltration in HCC patients and investigate its relationship with the risk scores. Gene set enrichment analysis (GSEA) was used to explore the potential mechanisms by which the COVID-19 genes affect hepatocarcinogenesis and prognosis. Three types of machine learning methods were combined to identify the most critical genes in the signature and localize their expression at the single cell level. We identified 53 prognosis-related COVID-19 genes and classified HCC into two molecular subtypes (C1, C2) by using the NMF method. The prognosis of C2 was significantly better than that of C1, and the two subtypes differed remarkably in terms of the tumor immune microenvironment and biological functions. The 17 COVID-19 genes were screened using the LASSO regression method to develop a 17 COVID-19 genes signature, which demonstrated a good predictive performance for 1-, 2- and 3-year OS of patients with HCC. The risk score as an independent prognostic factor for HCC has better predictive accuracy than traditional clinical variables. Patients in the TCGA cohort were categorized by risk score into the high- and low-risk groups, with the high-risk group mainly enriched in the immune modulation-related pathways and the low-risk group mainly enriched in the metabolism-related pathways, suggesting that the COVID-19 genes may affect disease progression and prognosis by regulating the tumor immune microenvironment and metabolism in HCC. NOL10 was identified as the most critical gene in the signature and hypothesized to be a potential therapeutic target for HCC. Objectively, the COVID-19 genes signature developed in this study, as an independent prognostic factor in HCC patients, is closely associated with the prognosis and tumor immune microenvironment of HCC patients and indicates that they may regulate the development of HCC in multiple ways, providing us with new perspectives for understanding the molecular mechanisms of HCC and finding effective therapeutic targets.

LINC00941 Promotes Cell Malignant Behavior and Is One of Five Costimulatory Molecule-Related lncRNAs That Predict Prognosis in Renal Clear Cell Carcinoma

Background and Objectives: A significant role was played by costimulatory molecules in renal cancer. However, the lncRNAs regulating costimulatory molecules have not been fully investigated. Materials and Methods: Data from the next-sequence file and clinical data were downloaded from the Cancer Genome Atlas (TCGA) database. All analyses were conducted using the R and GraphPad Prism software. Results: A total of 1736 costimulatory molecule-related lncRNAs were determined under the threshold of |Cor| > 0.5 and p-value < 0.001. Furthermore, a prognosis prediction signature consisting of five lncRNAs: LINC00941, AC016773.1, AL162171.1, HOTAIRM1, and AL109741.1 was established with great prediction ability. By combining risk score and clinical parameters, a nomogram plot was constructed for better clinical practice. A biological enrichment analysis indicated that E2F targets, coagulation, IL6/JAK/STAT3 signaling, G2/M checkpoint, and allograft rejection pathways were activated in high-risk patients. Furthermore, a higher infiltration level of resting CD4+ T cell, M2 macrophage, and resting mast cells, while a lower CD8+ T cell infiltration was observed in high-risk patients. It is worthy of note that, low-risk patients might respond better to PD-1 checkpoint therapy. A correlation analysis of LINC00941 revealed that it was positively correlated with Th2 cells, Th1 cells, macrophages, and Treg cells, but negatively correlated with Th17 cells. A pathway enrichment analysis indicated that the pathways of the inflammatory response, G2M checkpoint, and IL6/JAK/STAT3 signaling were significantly activated in patients with high LINC00941 expression. In vitro experiments indicated that LINC00941 can enhance the malignant biological behaviors of renal cancer cells. Conclusions: Our study established a costimulatory molecule-related lncRNAs-based prognosis model with a great prediction prognosis. In addition, LINC00941 could enhance the malignant biological behaviors of renal cancer cells.

Comprehensive Analysis of the E2F Transcription Factor Family in Human Lung Adenocarcinoma

Background

E2F transcription factors (E2Fs), code a family of pivotal transcription factors, have been identified as key regulators in tumor tumorigenesis. However, the function of E2F family in human lung adenocarcinoma (LUAD) have not been fully elucidated.

Methods

Herein, The Cancer Genome Atlas (TCGA) databases, Kaplan-Meier plotter, cBioPortal and TIMER were used to analyze differential expression, prognostic value, genetic alteration and immune cell infiltration of E2Fs in LUAD patients.

Results

The expression levels of E2Fs (E2F1-8) were all significantly upregulated in LUAD tissues compared with normal lung tissues. All eight E2Fs had low rates of gene mutation in LUAD patients from cBioPortal databases. Survival analysis revealed that E2F2 (P=0.038; HR 1.36; 95% CI 1.02–1.81), E2F7 (P<0.001; HR 1.78; 95% CI 1.33–2.39) and E2F8 (P=0.03; HR 1.37; 95% CI 1.02–1.82) were significantly associated with poor prognosis. Multivariate cox regression analysis found that only E2F7 (P<0.001; HR 2.72; 95% CI 1.75–4.25) was an independent prognostic predictor in LUAD after adjusting common clinical parameters. The receiver operating characteristic (ROC) analysis also found that E2F7 had high diagnostic value for LUAD (AUC=0.901). Further analysis found that E2F7 was significantly associated with LUAD immune cell infiltration of B cell, T cell, neutrophil, and myeloid dendritic cell. E2F7 also have positive correlations with immune checkpoint genes including SIGLEC15, CD274, HAVCR2, PDCD1LG2, CTLA4, TIGIT, LAG3 and PDCD1 in LUAD.

Conclusion

Our findings showed various association of E2F7 in LUAD diagnostic and prognostic aspects, which suggested its potential in becoming a novel biomarker.

An integrated analysis of DNA promoter methylation, microRNA regulation, and gene expression in gastric adenocarcinoma

Background

Gastric adenocarcinoma (GAC), a common type of gastric cancer, poses a significant public health threat worldwide. This study aimed to determine the transcriptional regulatory mechanisms of GAC.

Methods

HTSeq-FPKM raw data were obtained from The Cancer Genome Atlas Stomach Adenocarcinoma data collection. Subsequently, the limma package in R was used to identify differentially expressed genes (DEGs). Differentially methylated genes (DMGs), DEGs, and differentially expressed microRNAs (miRNAs) in normal, and tumor tissues of the same patients were screened and compared using R software tools. A functional enrichment analysis was performed using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) for various DEGs, DMGs, promoter methylation, and miRNAs. DEG-specific methylation and transcription factors were analyzed using ENCODE ChIP-seq.

Results

DEGs were centrally modified by the histone trimethylation of lysine 27 on histone H3 (H3K27me3). Upstream transcription factors of DEGs were enriched in different ChIP-seq clusters, such as Forkhead Box M1, E2F Transcription Factor 4, and suppressor of zest 12. Integrated regulatory networks of DEGs, promoter methylation, and miRNAs were constructed. Two miRNAs (hsa-mir-1 and hsa-mir-133a) and four DEGs (A disintegrin and metalloproteinase domain 12, transcription factor AP-2 alpha, solute carrier family 5 member 7, and cadherin 19) separately played important roles in the integrated regulatory network. Therefore, these DEGs, DMGs, promoter methylation, and miRNAs may play an important role in GAC pathogenesis.

Conclusions

In summary, the present study results provide insights into the oncogenesis and progression of GAC, thus accelerating the development of novel targeted GAC therapies.

Identification and characterization of a 25-lncRNA prognostic signature for early recurrence in hepatocellular carcinoma

Background

Early recurrence is the major cause of poor prognosis in hepatocellular carcinoma (HCC). Long non-coding RNAs (lncRNAs) are deeply involved in HCC prognosis. In this study, we aimed to establish a prognostic lncRNA signature for HCC early recurrence.

Methods

The lncRNA expression profile and corresponding clinical data were retrieved from total 299 HCC patients in TCGA database. LncRNA candidates correlated to early recurrence were selected by differentially expressed gene (DEG), univariate Cox regression and least absolute shrinkage and selection operator (LASSO) regression analyses. A 25-lncRNA prognostic signature was constructed according to receiver operating characteristic curve (ROC). Kaplan-Meier and multivariate Cox regression analyses were used to evaluate the performance of this signature. ROC and nomogram were used to evaluate the integrated models based on this signature with other independent clinical risk factors. Gene set enrichment analysis (GSEA) was used to reveal enriched gene sets in the high-risk group. Tumor infiltrating lymphocytes (TILs) levels were analyzed with single sample Gene Set Enrichment Analysis (ssGSEA). Immune therapy response prediction was performed with TIDE and SubMap. Chemotherapeutic response prediction was conducted by using Genomics of Drug Sensitivity in Cancer (GDSC) pharmacogenomics database.

Results

Compared to low-risk group, patients in high-risk group showed reduced disease-free survival (DFS) in the training (p < 0.0001) and validation cohort (p = 0.0132). The 25-lncRNA signature, AFP, TNM and vascular invasion could serve as independent risk factors for HCC early recurrence. Among them, the 25-lncRNA signature had the best predictive performance, and combination of those four risk factors further improves the prognostic potential. Moreover, GSEA showed significant enrichment of “E2F TARGETS”, “G2M CHECKPOINT”, “MYC TARGETS V1” and “DNA REPAIR” pathways in the high-risk group. In addition, increased TILs were observed in the low-risk group compared to the high-risk group. The 25-lncRNA signature negatively associates with the levels of some types of antitumor immune cells. Immunotherapies and chemotherapies prediction revealed differential responses to PD-1 inhibitor and several chemotherapeutic drugs in the low- and high-risk group.

Conclusions

Our study proposed a 25-lncRNA prognostic signature for predicting HCC early recurrence, which may guide postoperative treatment and recurrence surveillance in HCC patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-08827-z.

Characteristic of molecular subtypes in lung adenocarcinoma based on m6A RNA methylation modification and immune microenvironment

BACKGROUND

Lung adenocarcinoma (LUAD) is a major subtype of lung cancer and closely associated with poor prognosis. N6-methyladenosine (m6A), one of the most predominant modifications in mRNAs, is found to participate in tumorigenesis. However, the potential function of m6A RNA methylation in the tumor immune microenvironment is still murky.

METHODS

The gene expression profile cohort and its corresponding clinical data of LUAD patients were downloaded from TCGA database and GEO database. Based on the expression of 21 m6A regulators, we identified two distinct subgroups by consensus clustering. The single-sample gene-set enrichment analysis (ssGSEA) algorithm was conducted to quantify the relative abundance of the fraction of 28 immune cell types. The prognostic model was constructed by Lasso Cox regression. Survival analysis and receiver operating characteristic (ROC) curves were used to evaluate the prognostic model.

RESULT

Consensus classification separated the patients into two clusters (clusters 1 and 2). Those patients in cluster 1 showed a better prognosis and were related to higher immune scores and more immune cell infiltration. Subsequently, 457 differentially expressed genes (DEGs) between the two clusters were identified, and then a seven-gene prognostic model was constricted. The survival analysis showed poor prognosis in patients with high-risk score. The ROC curve confirmed the predictive accuracy of this prognostic risk signature. Besides, further analysis indicated that there were significant differences between the high-risk and low-risk groups in stages, status, clustering subtypes, and immunoscore. Low-risk group was related to higher immune score, more immune cell infiltration, and lower clinical stages. Moreover, multivariate analysis revealed that this prognostic model might be a powerful prognostic predictor for LUAD. Ultimately, the efficacy of this prognostic model was successfully validated in several external cohorts (GSE30219, GSE50081 and GSE72094).

CONCLUSION

Our study provides a robust signature for predicting patients' prognosis, which might be helpful for therapeutic strategies discovery of LUAD.

ChIP-sequencing analysis of E2F transcription factor 2 reveals its role in various biological processes of rheumatoid arthritis synovial fibroblasts

The development and progression of rheumatoid arthritis (RA) are complex and the pathogenesis of this disease is not fully understood. E2F transcription factor 2 (E2F2) affects the development and progression of many diseases. To identify the mechanisms underlying the role of E2F2 in RA, chromatin immunoprecipitation was performed followed by sequencing (ChIP-seq) using the E2F2 antibody. Gene Ontology (GO) analysis of differentially expressed genes (DEGs) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment of captured downstream target genes and Metascape analysis of 22 protein molecules partly elucidated the mechanism by which E2F2 affects the progression of RA. Results indicated that E2F2 affects the metabolism of RASFs and the development of ribosome synthesis as well as the stress response. Results indicated that E2F2 can affect multiple biological processes involving RASFs and indicate a unique possibility of targeting E2F2 in the treatment of RA.

Identification of Spindle and Kinetochore-Associated Family Genes as Therapeutic Targets and Prognostic Biomarkers in Pancreas Ductal Adenocarcinoma Microenvironment

Aim

The role of spindle and kinetochore-associated (SKA) genes in tumorigenesis and cancer progression has been widely studied. However, so far, the oncogenic involvement of SKA family genes in pancreatic cancer and their prognostic potential remain unknown.

Methods

Here, we carried out a meta-analysis of the differential expression of SKA genes in normal and tumor tissue. Univariate and multivariate survival analyses were done to evaluate the correlation between SKA family gene expression and pancreas ductal adenocarcinoma (PDAC) prognosis. Joint-effect and stratified survival analysis as well as nomogram analysis were used to estimate the prognostic value of genes. The underlying regulatory and biological mechanisms were identified by Gene set enrichment analysis. Interaction between SKA prognosis-related genes and immune cell infiltration was assessed using the Tumor Immune Estimation Resource tool.

Results

We find that SKA1-3 are highly expressed in PDAC tissues relative to non-cancer tissues. Survival analysis revealed that high expression of SKA1 and SKA3 independently indicate poor prognosis but they are not associated with relapse-free survival. The prognostic value of SKA1 and SKA3 was further confirmed by the nomogram, joint-effect, and stratified survival analysis. Analysis of underlying mechanisms reveals that these genes influence cancer-related signaling pathways, kinases, miRNA, and E2F family genes. Notably, prognosis-related genes are inversely correlated with several immune cells infiltrating levels.

Conclusion

We find that SKA1 and SKA3 expression correlates with prognosis and immune cell infiltration in PDAC, highlighting their potential as pancreatic cancer prognostic biomarkers.

Haprolid Inhibits Tumor Growth of Hepatocellular Carcinoma through Rb/E2F and Akt/mTOR Inhibition

BACKGROUND

Hepatocellular carcinoma (HCC) represents a major health burden with limited curative treatment options. There is a substantial unmet need to develop innovative approaches to impact the progression of advanced HCC. Haprolid is a novel natural component isolated from myxobacteria. Haprolid has been reported as a potent selective cytotoxin against a panel of tumor cells in recent studies including HCC cells. The aims of this study are to evaluate the antitumor effect of haprolid in HCC and to understand its underlying molecular mechanisms.

METHODS

The efficacy of haprolid was evaluated in human HCC cell lines (Huh-7, Hep3B and HepG2) and xenograft tumors (NMRI-Foxn1^nu mice with injection of Hep3B cells). Cytotoxic activity of haprolid was determined by the WST-1 and crystal violet assay. Wound healing, transwell and tumorsphere assays were performed to investigate migration and invasion of HCC cells. Apoptosis and cell-cycle distribution were measured by flow cytometry. The effects of haprolid on the Rb/E2F and Akt/mTOR pathway were examined by immunoblotting and immunohistochemistry.

RESULTS

haprolid treatment significantly inhibited cell proliferation, migration and invasion in vitro. The epithelial-mesenchymal transition (EMT) was impaired by haprolid treatment and the expression level of N-cadherin, vimentin and Snail was downregulated. Moreover, growth of HCC cells in vitro was suppressed by inhibition of G1/S transition, and partially by induction of apoptosis. The drug induced downregulation of cell cycle regulatory proteins cyclin A, cyclin B and CDK2 and induced upregulation of p21 and p27. Further evidence showed that these effects of haprolid were associated with Rb/E2F downregulation and Akt/mTOR inhibition. Finally, in vivo nude mice experiments demonstrated significant inhibition of tumor growth upon haprolid treatment.

CONCLUSION

Our results show that haprolid inhibits the growth of HCC through dual inhibition of Rb/E2F and Akt/mTOR pathways. Therefore, haprolid might be considered as a new and promising candidate for the palliative therapy of HCC.

Long noncoding RNA SNHG6 promotes proliferation and angiogenesis of cholangiocarcinoma cells through sponging miR-101-3p and activation of E2F8

Cholangiocarcinoma (CCA) development is an extremely complex process with alterations occurring in numerous genes. SNHG6, a validated lncRNA, has been reported to regulate the expression of multiple tumor-related genes in hepatocellular carcinoma, colorectal cancer and breast cancer. Here, we elucidated the function and possible molecular mechanisms of SNHG6 in human CCA cells. Our results proved that the expression SNHG6 was upregulated in CCA tissues and cell lines. Ectopic expression of SNHG6 promoted cell proliferation, cell cycle progression, migration, and angiogenesis in CCA cells, whereas knockdown of SNHG6 repressed these cellular processes. Further mechanistic studies revealed that SNHG6 could compete with the transcription factor E2F8 to bind with miR-101-3p, thus affecting E2F8 expression. Taken together, these results provided a comprehensive analysis of the role of SNHG6 in CCA cells and offered important clues to understand the key roles of competing endogenous RNA (ceRNA) mechanisms in human cholangiocarcinoma.

Identification of microRNAs as novel biomarkers for esophageal squamous cell carcinoma: a study based on The Cancer Genome Atlas (TCGA) and bioinformatics

BACKGROUND

MicroRNAs (miRNAs) have played important roles in the regulation of gene expression in many cancers, but their roles in esophageal squamous cell carcinoma (ESCC) are still unclear. The aim of this study was to determine the potential ESCC-specific key miRNAs from a large sample dataset in The Cancer Genome Atlas (TCGA).

METHODS

Integrative bioinformatics analysis was used to identify key ESCC-specific miRNAs related to the ESCC patients' tumor histological grade and lymphatic metastasis from TCGA. Next, these key miRNA potential gene regulatory functions and relationships with ESCC patients' clinical characteristics and overall survival were analyzed. Finally, three key miRNAs were selected randomly and quantificational real-time polymerase chain reaction (qRT-PCR) was used to validate in 51 newly diagnosed ESCC patients' tissues samples (collected from Nov. 2017 to Feb. 2019, in Wuwei, China) whether the bioinformatics analyses results were reliable and valid. Two-tailed Student's t test, Pearson Chi-squared test and Kaplan-Meier survival analysis were used in this study.

RESULTS

Thirty-five ESCC-specific miRNAs from TCGA database were investigated (fold-change > 2.0, P < 0.05), and 28 participated in the miRNAs-mRNAs co-expression network construction, while 17 were related with ESCC patients' tumor histological grade, TNM stage, and lymphatic metastasis (P < 0.05). Meanwhile, six miRNAs (including miR-200b-3p, miR-31-5p, miR-15b-5p, miR-141-3p, miR-135b-5p, and miR-195-5p) were correlated with overall survival of ESCC patients (log-rank, P < 0.05). MiR-135b-5p, miR-15b-5p, and miR-195-5p were selected for verification of the expression levels in 51 ESCC patients' tissue samples by using qRT-PCR. We found that the fold-changes between qRT-PCR and TCGA were completely consistent. The results also suggested that miR-135b-5p, miR-15b-5p, and miR-195-5p were significantly correlated with tumor differentiation degrees (P < 0.05), miR-195-5p was significantly correlated with tumor TNM stage (P < 0.05), and miR-135b-5p was significantly correlated with lymph-node metastasis (P < 0.05). MiR-135b-5p, miR-15b-5p, and miR-195-5p expression levels, ESCC patient clinical features association analysis results and the aforementioned TCGA bioinformatics analyses were similar.

CONCLUSION

This study identified key ESCC-related miRNAs. The key miRNAs are worthy of further investigation as potential novel biomarkers for diagnosis, classification, and prognosis of ESCC.

A 5-Gene Prognostic Combination for Predicting Survival of Patients with Gastric Cancer

Background

The aim of the study was to identify a multigene prognostic factor in patients with gastric cancer (GC).

Material/Methods

Random survival forest (RSF) was performed to screen survival-related genes and develop a multigene combination based on the cumulative hazard function of each GC patient in TCGA-STAD and GSE15459. Kaplan-Meier curve and univariate and multivariable Cox proportional hazards regression model were applied to evaluate the prognostic performance of the 5-gene combination. C-index was used to compare the prognostic performance of the 5-gene combination and another 9-gene signature in GC. Gene set enrichment analysis (GSEA) was conducted.

Results

We obtained 19 survival-related genes through univariate Cox proportional hazards analysis in the training set, 5 of which were identified and were used to develop a 5-gene combination through RSF. Patients in the 5-gene combination low-risk group had better overall survival (OS) than those in the 5-gene combination high-risk group, and the 5-gene combination was demonstrated to be an independent prognostic factor in patients with GC. The 5-gene combination outperformed the 9-gene signature in predicting the OS of GC patients, and it might affect the prognosis of GC patients through E2F signaling, MYC signaling, and G2M checkpoint.

Conclusions

We introduce a 5-gene combination that can predict the survival of GC patients and might be an independent prognostic factor in GC.

Genetic Variability as a Regulator of TLR4 and NOD Signaling in Response to Bacterial Driven DNA Damage Response (DDR) and Inflammation: Focus on the Gastrointestinal (GI) Tract

The fundamental role of human Toll-like receptors (TLRs) and NOD-like receptors (NLRs), the two most studied pathogen recognition receptors (PRRs), is the protection against pathogens and excessive tissue injury. Recent evidence supports the association between TLR/NLR gene mutations and susceptibility to inflammatory, autoimmune, and malignant diseases. PRRs also interfere with several cellular processes, such as cell growth, apoptosis, cell proliferation, differentiation, autophagy, angiogenesis, cell motility and migration, and DNA repair mechanisms. We briefly review the impact of TLR4 and NOD1/NOD2 and their genetic variability in the process of inflammation, tumorigenesis and DNA repair, focusing in the gastrointestinal tract. We also review the available data on new therapeutic strategies utilizing TLR/NLR agonists and antagonists for cancer, allergic diseases, viral infections and vaccine development against both infectious diseases and cancer.

Analysis of cis-Regulatory Elements in Gene Co-expression Networks in Cancer

Analysis of gene co-expression networks is a powerful "data-driven" tool, invaluable for understanding cancer biology and mechanisms of tumor development. Yet, despite of completion of thousands of studies on cancer gene expression, there were few attempts to normalize and integrate co-expression data from scattered sources in a concise "meta-analysis" framework. Here we describe an integrated approach to cancer expression meta-analysis, which combines generation of "data-driven" co-expression networks with detailed statistical detection of promoter sequence motifs within the co-expression clusters. First, we applied Weighted Gene Co-Expression Network Analysis (WGCNA) workflow and Pearson's correlation to generate a comprehensive set of over 3000 co-expression clusters in 82 normalized microarray datasets from nine cancers of different origin. Next, we designed a genome-wide statistical approach to the detection of specific DNA sequence motifs based on similarities between the promoters of similarly expressed genes. The approach, realized as cisExpress software module, was specifically designed for analysis of very large data sets such as those generated by publicly accessible whole genome and transcriptome projects. cisExpress uses a task farming algorithm to exploit all available computational cores within a shared memory node.We discovered that although co-expression modules are populated with different sets of genes, they share distinct stable patterns of co-regulation based on promoter sequence analysis. The number of motifs per co-expression cluster varies widely in accordance with cancer tissue of origin, with the largest number in colon (68 motifs) and the lowest in ovary (18 motifs). The top scored motifs are typically shared between several tissues; they define sets of target genes responsible for certain functionality of cancerogenesis. Both the co-expression modules and a database of precalculated motifs are publically available and accessible for further studies.

Prostate cancer-associated mutation in SPOP impairs its ability to target Cdc20 for poly-ubiquitination and degradation

Recent studies revealed that mutations in SPOP (Speckle-type POZ protein) occur in up to 15% of patients with prostate cancer. However, the physiological role of SPOP in regulating prostate tumorigenesis remains elusive. Here, we identified the Cdc20 oncoprotein as a novel ubiquitin substrate of SPOP. As such, pharmacological inhibition of Cullin-based E3 ligases by MLN4924 could stabilize endogenous Cdc20 in cells. Furthermore, we found that Cullin 3, and, to a less extent, Cullin 1, specifically interacted with Cdc20. Depletion of Cullin 3, but not Cullin 1, could upregulate the abudance of Cdc20 largely via prolonging Cdc20 half-life. Moreover, SPOP, the adaptor protein of Cullin 3 family E3 ligase, specifically interacted with Cdc20, and promoted the poly-ubiquitination and subsequent degradation of Cdc20 in a degron-dependent manner. Importantly, prostate cancer-derived SPOP mutants failed to interact with Cdc20 to promote its degradation. As a result, SPOP-deficient prostate cancer cells with elevated Cdc20 expression became resistant to a pharmacological Cdc20 inhibitor. Therefore, our results revealed a novel role of SPOP in tumorigenesis in part by promoting the degradation of the Cdc20 oncoprotein.

Mammalian Metallothionein-2A and Oxidative Stress

Mammalian metallothionein-2A (MT2A) has received considerable attention in recent years due to its crucial pathophysiological role in anti-oxidant, anti-apoptosis, detoxification and anti-inflammation. For many years, most studies evaluating the effects of MT2A have focused on reactive oxygen species (ROS), as second messengers that lead to oxidative stress injury of cells and tissues. Recent studies have highlighted that oxidative stress could activate mitogen-activated protein kinases (MAPKs), and MT2A, as a mediator of MAPKs, to regulate the pathogenesis of various diseases. However, the molecule mechanism of MT2A remains elusive. A deeper understanding of the functional, biochemical and molecular characteristics of MT2A would be identified, in order to bring new opportunities for oxidative stress therapy.

Elevated E2F7 expression predicts poor prognosis in human patients with gliomas

E2F transcription factors have been studied extensively in a broad range of organisms as major regulators of cell cycle, apoptosis, and differentiation. The E2F family includes the atypical member E2F7, which has been rarely studied in gliomas. The aim of this study is to determine the expression status of E2F7 in gliomas, its relationship to clinicopathological features, and patients' outcome. The mRNA levels of E2F7 in the human brain and different grades of gliomas were analysed using datasets from the publically available Oncomine database. One of the most significant co-expression factors, CDK1, together with E2F7, was further validated by immunohistochemistry in 90 different grades of gliomas. Furthermore, univariate and multivariate analyses were performed to identify prognostic variables relative to patient and tumour characteristics and treatment modalities. E2F7 mRNA expression was found to be elevated in gliomas by Oncomine-database analysis. Immunohistochemistry showed an increase in E2F7 labelling index in high- versus low-grade gliomas (62.1±11.8% vs. 18.9±10.2%, p<0.0001). There was a positive correlation between E2F7 and CDK1 immunoreactivity (Spearman r=0.446, p=0.037). Clinicopathological evaluation suggested that E2F7 expression was associated with tumour grade (p<0.0001) and recurrence (p=0.025). In Cox multivariate analysis, pathological classification and recurrence were independent prognostic factors of gliomas, and E2F7 was significantly related to progression-free survival (p=0.011), but not overall survival (p=0.062). Our findings suggested that E2F7 might act as an independent prognostic factor of gliomas and might constitute a potential therapeutic target for this disease.

Exploring and exploiting the systemic effects of deregulated replication licensing

Maintenance and accurate propagation of the genetic material are key features for physiological development and wellbeing. The replication licensing machinery is crucial for replication precision as it ensures that replication takes place once per cell cycle. Thus, the expression status of the components comprising the replication licensing apparatus is tightly regulated to avoid re-replication; a form of replication stress that leads to genomic instability, a hallmark of cancer. In the present review we discuss the mechanistic basis of replication licensing deregulation, which leads to systemic effects, exemplified by its role in carcinogenesis and a variety of genetic syndromes. In addition, new insights demonstrate that above a particular threshold, the replication licensing factor Cdc6 acts as global transcriptional regulator, outlining new lines of exploration. The role of the putative replication licensing factor ChlR1/DDX11, mutated in the Warsaw Breakage Syndrome, in cancer is also considered. Finally, future perspectives focused on the potential therapeutic advantage by targeting replication licensing factors, and particularly Cdc6, are discussed.

Dickkopf-1-promoted vasculogenic mimicry in non-small cell lung cancer is associated with EMT and development of a cancer stem-like cell phenotype

To characterize the contributions of Dickkopf‐1 (DKK1) towards the induction of vasculogenic mimicry (VM) in non‐small cell lung cancer (NSCLC), we evaluated cohorts of primary tumours, performed in vitro functional studies and generated xenograft mouse models. Vasculogenic mimicry was observed in 28 of 205 NSCLC tumours, while DKK1 was detected in 133 cases. Notably, DKK1 was positively associated with VM. Statistical analysis showed that VM and DKK1 were both related to aggressive clinical course and thus were indicators of a poor prognosis. Moreover, expression of epithelial‐mesenchymal transition (EMT)‐related proteins (vimentin, Slug, and Twist), cancer stem‐like cell (CSC)‐related proteins (nestin and CD44), VM‐related proteins (MMP2, MMP9, and vascular endothelial‐cadherin), and β‐catenin‐nu were all elevated in VM‐positive and DKK1‐positive tumours, whereas the epithelial marker (E‐cadherin) was reduced in the VM‐positive and DKK1‐positive groups. Non‐small cell lung cancer cell lines with overexpressed or silenced DKK1 highlighted its role in the restoration of mesenchymal phenotypes and development of CSC characteristics. Moreover, DKK1 significantly promotes NSCLC tumour cells to migrate, invade and proliferate. In vivo animal studies demonstrated that DKK1 enhances the growth of transplanted human tumours cells, as well as increased VM formation, mesenthymal phenotypes and CSC properties. Our results suggest that DKK1 can promote VM formation via induction of the expression of EMT and CSC‐related proteins. As such, we feel that DKK1 may represent a novel target of NSCLC therapy.

An ANCCA/PRO2000-miR-520a-E2F2 regulatory loop as a driving force for the development of hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is one of the most common malignancies in Asia especially in China. We previously identified that ANCCA/PRO2000 as an important proliferation-associated protein predicted poor prognosis of patients with HCC. However, the molecular mechanisms of ANCCA/PRO2000 leading to hepatocarcinogenesis and progression are still obscure. In the present study, we found that ANCCA/PRO2000 overexpression in HCC specimens correlated with aggressive tumor behavior and poor survival. Furthermore, ANCCA/PRO2000 exerts strong oncogenic function in HCC and promotes cell proliferation by regulating E2F2 expression, a critical cell cycle regulator. Notably, miR-520a is an intermediate regulator between ANCCA/PRO2000 and E2F2. Mechanistically, ANCCA/PRO2000 not only interacts with E2F2 but also negatively regulates miR-520a that inhibits E2F2 to cooperatively promote in vitro and in vivo growth of HCC cells. Moreover, we demonstrated that ANCCA/PRO2000 enhances the migratory capacity of HCC cells partially by suppressing ERO1L and G3BP2 expression. Additional research identified that miR-372, as a prognostic factor for HCC, could directly target ANCCA/PRO2000. Our results suggest the ANCCA/PRO2000-miR-520a-E2F2 regulatory loop as a driving force for HCC development and ANCCA/PRO2000 as a potential therapeutic target for HCC.

Apoptosis or senescence? Which exit route do epithelial cells and fibroblasts preferentially follow?

Senescence and apoptosis constitute types of cellular responses that normally ensure homeostasis, when endogenous or exogenous signals occur. Their deregulation is often observed in various pathologies, such as age and non-age related diseases including cancer. Although epithelial cells and fibroblasts are capable to exert both functions, under a plethora of insults, the fact that they exhibit notable intrinsic differences in cell/tissue homeostasis properties, might be a crucial determinant of the mode of response to a certain stress signal. Sparse evidence in the literature reveals that in the same tissue/organ context and under the same conditions, the cell type seems to drive the differential counteraction between epithelia and fibroblasts. Based on the above notion we propose that, upon stress insults, human fibroblasts seem to predominantly respond via senescence, while epithelial cells prefer to exert apoptosis. We suggest that considering the tissue as a whole (epithelium and stroma) would benefit research into new therapeutic strategies for chronic diseases and cancer.

Effects of hyperthermia as a mitigation strategy in DNA damage-based cancer therapies

Utilization of thermal therapy (hyperthermia) is defined as the application of exogenous heat induction and represents a concept that is far from new as it goes back to ancient times when heat was used for treating various diseases, including malignancies. Such therapeutic strategy has gained even more popularity (over the last few decades) since various studies have shed light into understanding hyperthermia's underlying molecular mechanism(s) of action. In general, hyperthermia is applied as complementary (adjuvant) means in therapeutic protocols combining chemotherapy and/or irradiation both of which can induce irreversible cellular DNA damage. Furthermore, according to a number of in vitro, in vivo and clinical studies, hyperthermia has been shown to enhance the beneficial effects of DNA targeting therapeutic strategies by interfering with DNA repair response cascades. Therefore, the continuously growing evidence supporting hyperthermia's beneficial role in cancer treatment can also encourage its application as a DNA repair mitigation strategy. In this review article, we aim to provide detailed information on how hyperthermia acts on DNA damage and repair pathways and thus potentially contributing to various adjuvant therapeutic protocols relevant to more efficient cancer treatment strategies.