Large-scale meta-analysis of cancer microarray data identifies common transcriptional profiles of neoplastic transformation and progression

Six drivers of aging identified among genes differentially expressed with age

Many studies have compared gene expression in young and old samples to gain insights on aging, the primary risk factor for most major chronic diseases. However, these studies only describe associations, failing to distinguish drivers of aging from compensatory geroprotective responses and incidental downstream effects. Here, we introduce a workflow to characterize the causal effects of differentially expressed genes on lifespan. First, we performed a meta-analysis of 25 gene expression datasets comprising samples of various tissues from healthy, untreated adult mammals (humans, dogs, and rodents) at two distinct ages. We ranked each gene according to the number of distinct datasets in which the gene was differentially expressed with age in a consistent direction. The top age-upregulated genes were TMEM176A, EFEMP1, CP, and HLA-A; the top age-downregulated genes were CA4, SIAH, SPARC, and UQCR10. Second, the effects of the top ranked genes on lifespan were measured by applying post-developmental RNA interference of the corresponding ortholog in the nematode C. elegans (two trials, with roughly 100 animals per genotype per trial). Out of 10 age-upregulated and 9 age-downregulated genes that were tested, two age-upregulated genes (csp-3/CASP1 and spch-2/RSRC1) and four age-downregulated genes (C42C1.8/DIRC2, ost-1/SPARC, fzy-1/CDC20, and cah-3/CA4) produced significant and reproducible lifespan extension. Notably, the data do not suggest that the direction of differential expression with age is predictive of the effect on lifespan. Our study provides novel insight into the relationship between differential gene expression and aging phenotypes, pilots an unbiased workflow that can be easily repeated and expanded, and pinpoints six genes with evolutionarily conserved, causal roles in the aging process for further study.

The roles and mechanisms of TGFB1 in acute myeloid leukemia chemoresistance

Relapsed or Refractory (R/R) Acute Myeloid Leukemia (AML) patients usually have very poor prognoses, and drug-resistance is one of the major limiting factors. In this study, we aimed to explore the functions of Transforming Growth Factor-β1 (TGFB1) in AML drug-resistance. First, TGFB1 levels in serum and bone marrow are higher in R/R patients compared with newly diagnosed patients, this phenomenon could be due to different sources of secreted TGFB1 according to immunohistochemistry of marrow biopsies. Similarly, TGFB1 expression in AML drug-resistant cell lines is higher than that in their parental cell lines, and blocking the TGFB signaling pathway by specific inhibitors decreased resistance to chemotherapeutic agents. On the other hand, exogenous TGFB1 can also promote AML parental cells senescence and chemotherapy resistance. Next, we found SOX4 level is upregulated in drug-resistant cells, and parental cells treated with exogenous TGFB1 induced upregulation of SOX4 levels. Interference of SOX4 expression by siRNA diminished the TGFB1-induced sensitivity to chemotherapeutic agents. Finally, we conduct metabolomic analysis and find Alanine, aspartate and glutamate metabolism pathway, and Glycerophospholipid metabolism pathway are decreased after inhibiting TGFB signaling pathway or interfering SOX4 expression. This study concludes that TGFB1 level in R/R AML patients and drug-resistant strains is significantly increased. Blocking the TGFB signaling pathway can enhance the chemosensitivity of drug-resistant cells by suppressing SOX4 expression and metabolic reprogramming.

Establishment and characterization of cisplatin-resistant cell lines from canine mammary gland tumors

The development of cisplatin resistance is one of the major causes of mammary cancer treatment failure, and is associated with changes in Sox4 gene expression. To investigate the characteristic changes that occur in canine mammary gland tumor (CMGT) cells following the development of acquired cisplatin resistance, along with the relationship between these changes and the Sox4 gene. We constructed cisplatin-resistant cell line, CHMpCIS, from the cell line CHMp, which was isolated from the primary lesion of a malignant CMGT. The biological characteristics of these cells were examined by Western blot analysis, Transwell assays, and mammosphere formation assays. Compared to CHMp cells, CHMpCIS cells exhibited elevated cisplatin resistance, apoptotic escape ability, enhanced epithelial-mesenchymal transition (EMT) and cancer stem cell (CSC) features, in addition to over-activation of the Wnt/β-catenin signaling pathway and increased Sox4 protein. In CMGT cases, CMGT tissues (CMGTT) expressed higher levels of Sox4 protein and mRNA compared to adjacent tissues (CAMGTT). We found that these changes were inhibited by silencing of Sox4 expression in CHMpCIS cells. Furthermore, activation of the Wnt/β-catenin signaling pathway increased Sox4 expression levels through a positive feedback loop. These results suggested that CHMpCIS cells circumvented the damage caused by cisplatin through altering the expression of the Sox4 gene and activating the Wnt/β-catenin pathway, thereby changing the cellular biological characteristics.

Spatial MS multiomics on clinical prostate cancer tissues

Mass spectrometry (MS) and MS imaging (MSI) are used extensively for both the spatial and bulk characterization of samples in lipidomics and proteomics workflows. These datasets are typically generated independently due to different requirements for sample preparation. However, modern omics technologies now provide higher sample throughput and deeper molecular coverage, which, in combination with more sophisticated bioinformatic and statistical pipelines, make generating multiomics data from a single sample a reality. In this workflow, we use spatial lipidomics data generated by matrix-assisted laser desorption/ionization MSI (MALDI-MSI) on prostate cancer (PCa) radical prostatectomy cores to guide the definition of tumor and benign tissue regions for laser capture microdissection (LCM) and bottom-up proteomics all on the same sample and using the same mass spectrometer. Accurate region of interest (ROI) mapping was facilitated by the SCiLS region mapper software and dissected regions were analyzed using a dia-PASEF workflow. A total of 5525 unique protein groups were identified from all dissected regions. Lysophosphatidylcholine acyltransferase 1 (LPCAT1), a lipid remodelling enzyme, was significantly enriched in the dissected regions of cancerous epithelium (CE) compared to benign epithelium (BE). The increased abundance of this protein was reflected in the lipidomics data with an increased ion intensity ratio for pairs of phosphatidylcholines (PC) and lysophosphatidylcholines (LPC) in CE compared to BE.

E2F1 Mediates SOX17 Deficiency-Induced Pulmonary Hypertension

BACKGROUND

Rare genetic variants and genetic variation at loci in an enhancer in SOX17 (SRY-box transcription factor 17) are identified in patients with idiopathic pulmonary arterial hypertension (PAH) and PAH with congenital heart disease. However, the exact role of genetic variants or mutations in SOX17 in PAH pathogenesis has not been reported.

METHODS

SOX17 expression was evaluated in the lungs and pulmonary endothelial cells (ECs) of patients with idiopathic PAH. Mice with Tie2Cre-mediated Sox17 knockdown and EC-specific Sox17 deletion were generated to determine the role of SOX17 deficiency in the pathogenesis of PAH. Human pulmonary ECs were cultured to understand the role of SOX17 deficiency. Single-cell RNA sequencing, RNA-sequencing analysis, and luciferase assay were performed to understand the underlying molecular mechanisms of SOX17 deficiency-induced PAH. E2F1 (E2F transcription factor 1) inhibitor HLM006474 was used in EC-specific Sox17 mice.

RESULTS

SOX17 expression was downregulated in the lung and pulmonary ECs from patients with idiopathic PAH. Mice with Tie2Cre-mediated Sox17 knockdown and EC-specific Sox17 deletion induced spontaneously mild pulmonary hypertension. Loss of endothelial Sox17 in EC exacerbated hypoxia-induced pulmonary hypertension in mice. Loss of SOX17 in lung ECs induced endothelial dysfunctions including upregulation of cell cycle programming, proliferative and antiapoptotic phenotypes, augmentation of paracrine effect on pulmonary arterial smooth muscle cells, impaired cellular junction, and BMP (bone morphogenetic protein) signaling. E2F1 signaling was shown to mediate the SOX17 deficiency-induced EC dysfunction. Pharmacological inhibition of E2F1 in Sox17 EC-deficient mice attenuated pulmonary hypertension development.

CONCLUSIONS

Our study demonstrated that endothelial SOX17 deficiency induces pulmonary hypertension through E2F1. Thus, targeting E2F1 signaling represents a promising approach in patients with PAH.

Prognostic significance of TM4SF1 and DDR1 expression in epithelial ovarian cancer

Transmembrane 4 L6 family member 1 (TM4SF1) and discoidin domain receptor 1 (DDR1) are expressed in numerous types of cancer, but their expression in epithelial ovarian cancer and the association between their expression and patient prognosis are unclear. The present study aimed to explore the expression of TM4SF1 and DDR1 and their relationship with prognosis in epithelial ovarian cancer. Firstly, the Oncomine and Gene Expression Profiling Interactive Analysis (GEPIA) platforms were used to compare the expression levels of TM4SF1 and DDR1 in ovarian cancer and normal ovarian tissue, and Kaplan-Meier plotter was used to analyze the association between gene expression and patient prognosis. The proteins interacting with TM4SF1 and DDR1 were analyzed using Search Tool for the Retrieval of Interacting Genes/Proteins (STRING), and enrichment analysis of Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathways was conducted for the interacting proteins. Furthermore, immunohistochemical staining was performed to detect the expression of TM4SF1 and DDR1 protein in epithelial ovarian cancer tissue and to analyze the association between expression and prognosis. The Oncomine and GEPIA analyses showed that the expression levels of TM4SF1 and DDR1 were significantly higher in epithelial ovarian cancer than in normal ovarian tissue, and the analysis of clinical samples revealed that TM4SF1 and DDR1 were coexpressed in some cases. STRING analysis indicated that the TM4SF1 and DDR1 proteins interact with each other. The overall survival and progression-free survival of patients whose epithelial ovarian cancer coexpressed TM4SF1 and DDR1 were significantly shorter than those of patients lacking TM4SF1 and DDR1 coexpression. Multivariate analysis indicated that TM4SF1 and DDR1 protein coexpression was an independent prognostic factor. In summary, TM4SF1 and DDR1 proteins were coexpressed in some epithelial ovarian cancer tissues and appear to be adverse prognostic factors for epithelial ovarian cancer. In addition, TM4SF1 and DDR1 may have an interactive or mutual regulatory mechanism.

Growth differentiation factor myostatin regulates epithelial-mesenchymal transition genes and enhances invasion by increasing serine protease inhibitors E1 and E2 in human trophoblast cells

Placental insufficiency disorders, including preeclampsia and intrauterine growth restriction, are major obstetric complications that can have devastating effects on both the mother and the fetus. These syndromes have underlying poor placental trophoblast cell invasion into uterine tissues. Placental invasion is controlled by many hormones and growth factors. Myostatin (MSTN) is a transforming growth factor-β superfamily member recognized for its important role in muscle growth control. MSTN has also been shown to be secreted and functioning in the placenta, and its serum and/or placental levels were found to be upregulated in preeclampsia and intrauterine growth restriction. Considering that the mechanistic role of MSTN in placentation remains poorly understood, we hypothesized that MSTN uses ALK4/5-SMAD2/3/4 signaling to increase human trophoblast invasion through a group of epithelial-mesenchymal transition genes including SERPINE2, PAI-1, and SOX4. mRNA sequencing of control and MSTN-treated primary human trophoblast cells (n = 5) yielded a total of 610 differentially expressed genes (false discovery rate <0.05) of which 380 genes were upregulated and 230 were downregulated. These differentially expressed genes were highly enriched in epithelial-mesenchymal transition genes, and a subset including SERPINE2, PAI-1, and SOX4 was investigated for its role in MSTN-induced trophoblast cell invasion. We found that MSTN induced upregulation of SERPINE2 via ALK4/5-SMAD2/3/4 signaling; however, SMAD2 was not involved in MSTN-induced PAI-1 upregulation. SOX4 was involved in MSTN-induced upregulation of SERPINE2, but not PAI-1. Collectively, this study discovers novel molecular mechanisms of MSTN-induced human trophoblast cell invasion and provides insight into the functional consequences of its dysregulation in placental insufficiency disorders.

SOX4-SMARCA4 complex promotes glycolysis-dependent TNBC cell growth through transcriptional regulation of Hexokinase 2

Tumor cells rely on increased glycolytic capacity to promote cell growth and progression. While glycolysis is known to be upregulated in the majority of triple negative (TNBC) or basal-like subtype breast cancers, the mechanism remains unclear. Here, we used integrative genomic analyses to identify a subset of basal-like tumors characterized by increased expression of the oncogenic transcription factor SOX4 and its co-factor the SWI/SNF ATPase SMARCA4. These tumors are defined by unique gene expression programs that correspond with increased tumor proliferation and activation of key metabolic pathways, including glycolysis. Mechanistically, we demonstrate that the SOX4-SMARCA4 complex mediates glycolysis through direct transcriptional regulation of Hexokinase 2 (HK2) and that aberrant HK2 expression and altered glycolytic capacity are required to mediate SOX4-SMARCA4-dependent cell growth. Collectively, we have defined the SOX4-SMARCA4-HK2 signaling axis in basal-like breast tumors and established that this axis promotes metabolic reprogramming which is required to maintain tumor cell growth.

UHRF1 inhibition epigenetically reprograms cancer stem cells to suppress the tumorigenic phenotype of hepatocellular carcinoma

Cancer stem cells (CSCs) contribute to tumor initiation, progression, and recurrence in many types of cancer, including hepatocellular carcinoma (HCC). Epigenetic reprogramming of CSCs has emerged as a promising strategy for inducing the transition from malignancy to benignity. Ubiquitin-like with PHD and ring finger domains 1 (UHRF1) is required for DNA methylation inheritance. Here, we investigated the role and mechanism of UHRF1 in regulating CSC properties and evaluated the impact of UHRF1 targeting on HCC. Hepatocyte-specific Uhrf1 knockout (Uhrf1^HKO) strongly suppressed tumor initiation and CSC self-renewal in both diethylnitrosamine (DEN)/CCl₄-induced and Myc-transgenic HCC mouse models. Ablation of UHRF1 in human HCC cell lines yielded consistent phenotypes. Integrated RNA-seq and whole genome bisulfite sequencing revealed widespread hypomethylation induced by UHRF1 silencing epigenetically reprogrammed cancer cells toward differentiation and tumor suppression. Mechanistically, UHRF1 deficiency upregulated CEBPA and subsequently inhibited GLI1 and Hedgehog signaling. Administration of hinokitiol, a potential UHRF1 inhibitor, significantly reduced tumor growth and CSC phenotypes in mice with Myc-driven HCC. Of pathophysiological significance, the expression levels of UHRF1, GLI1, and key axis proteins consistently increased in the livers of mice and patients with HCC. These findings highlight the regulatory mechanism of UHRF1 in liver CSCs and have important implications for the development of therapeutic strategies for HCC.

E2F1 Mediates SOX17 Deficiency-Induced Pulmonary Hypertension

Rationale

Rare genetic variants and genetic variation at loci in an enhancer in SRY-Box Transcription Factor 17 (SOX17) are identified in patients with idiopathic pulmonary arterial hypertension (PAH) and PAH with congenital heart disease. However, the exact role of genetic variants or mutation in SOX17 in PAH pathogenesis has not been reported.

Objectives

To investigate the role of SOX17 deficiency in pulmonary hypertension (PH) development.

Methods

Human lung tissue and endothelial cells (ECs) from IPAH patients were used to determine the expression of SOX17. Tie2Cre-mediated and EC-specific deletion of Sox17 mice were assessed for PH development. Single-cell RNA sequencing analysis, human lung ECs, and smooth muscle cell culture were performed to determine the role and mechanisms of SOX17 deficiency. A pharmacological approach was used in Sox17 deficiency mice for therapeutic implication.

Measurement and Main Results

SOX17 expression was downregulated in the lungs and pulmonary ECs of IPAH patients. Mice with Tie2Cre mediated Sox17 knockdown and EC-specific Sox17 deletion developed spontaneously mild PH. Loss of endothelial Sox17 in EC exacerbated hypoxia-induced PH in mice. Loss of SOX17 in lung ECs induced endothelial dysfunctions including upregulation of cell cycle programming, proliferative and anti-apoptotic phenotypes, augmentation of paracrine effect on pulmonary arterial smooth muscle cells, impaired cellular junction, and BMP signaling. E2F Transcription Factor 1 (E2F1) signaling was shown to mediate the SOX17 deficiency-induced EC dysfunction and PH development.

Conclusions

Our study demonstrated that endothelial SOX17 deficiency induces PH through E2F1 and targeting E2F1 signaling represents a promising approach in PAH patients.

Synthesizing genome regulation data with vote-counting

The increasing availability of high-throughput datasets allows amalgamating research information across a large body of genome regulation studies. Given the recent success of meta-analyses on transcriptional regulators, epigenetic marks, and enhancer:gene associations, we expect that such surveys will continue to provide novel and reproducible insights. However, meta-analyses are severely hampered by the diversity of available data, concurring protocols, an eclectic amount of bioinformatics tools, and myriads of conceivable parameter combinations. Such factors can easily bar life scientists from synthesizing omics data and substantially curb their interpretability. Despite statistical challenges of the method, we would like to emphasize the advantages of joining data from different sources through vote-counting and showcase examples that achieve a simple but highly intuitive data integration.

Population structure discovery in meta-analyzed microbial communities and inflammatory bowel disease using MMUPHin

Microbiome studies of inflammatory bowel diseases (IBD) have achieved a scale for meta-analysis of dysbioses among populations. To enable microbial community meta-analyses generally, we develop MMUPHin for normalization, statistical meta-analysis, and population structure discovery using microbial taxonomic and functional profiles. Applying it to ten IBD cohorts, we identify consistent associations, including novel taxa such as Acinetobacter and Turicibacter, and additional exposure and interaction effects. A single gradient of dysbiosis severity is favored over discrete types to summarize IBD microbiome population structure. These results provide a benchmark for characterization of IBD and a framework for meta-analysis of any microbial communities.

Signatures of Co-Deregulated Genes and Their Transcriptional Regulators in Lung Cancer

Despite the significant progress made towards comprehending the deregulated signatures in lung cancer, these vary from study to study. We reanalyzed 25 studies from the Gene Expression Omnibus (GEO) to detect and annotate co-deregulated signatures in lung cancer and in single-gene or single-drug perturbation experiments. We aimed to decipher the networks that these co-deregulated genes (co-DEGs) form along with their upstream regulators. Differential expression and upstream regulators were computed using Characteristic Direction and Systems Biology tools, including GEO2Enrichr and X2K. Co-deregulated gene expression profiles were further validated across different molecular and immune subtypes in lung adenocarcinoma (TCGA-LUAD) and lung adenocarcinoma (TCGA-LUSC) datasets, as well as using immunohistochemistry data from the Human Protein Atlas, before being subjected to subsequent GO and KEGG enrichment analysis. The functional alterations of the co-upregulated genes in lung cancer were mostly related to immune response regulating the cell surface signaling pathway, in contrast to the co-downregulated genes, which were related to S-nitrosylation. Networks of hub proteins across the co-DEGs consisted of overlapping TFs (SOX2, MYC, KAT2A) and kinases (MAPK14, CSNK2A1 and CDKs). Furthermore, using Connectivity Map we highlighted putative repurposing drugs, including valproic acid, betonicine and astemizole. Similarly, we analyzed the co-DEG signatures in single-gene and single-drug perturbation experiments in lung cancer cell lines. In summary, we identified critical co-DEGs in lung cancer providing an innovative framework for their potential use in developing personalized therapeutic strategies.

Connecting omics signatures and revealing biological mechanisms with iLINCS

There are only a few platforms that integrate multiple omics data types, bioinformatics tools, and interfaces for integrative analyses and visualization that do not require programming skills. Here we present iLINCS ( http://ilincs.org ), an integrative web-based platform for analysis of omics data and signatures of cellular perturbations. The platform facilitates mining and re-analysis of the large collection of omics datasets (>34,000), pre-computed signatures (>200,000), and their connections, as well as the analysis of user-submitted omics signatures of diseases and cellular perturbations. iLINCS analysis workflows integrate vast omics data resources and a range of analytics and interactive visualization tools into a comprehensive platform for analysis of omics signatures. iLINCS user-friendly interfaces enable execution of sophisticated analyses of omics signatures, mechanism of action analysis, and signature-driven drug repositioning. We illustrate the utility of iLINCS with three use cases involving analysis of cancer proteogenomic signatures, COVID 19 transcriptomic signatures and mTOR signaling.

Evolutionary Dynamics and Molecular Mechanisms of HORMA Domain Protein Signaling

Controlled assembly and disassembly of multi-protein complexes is central to cellular signaling. Proteins of the widespread and functionally diverse HORMA family nucleate assembly of signaling complexes by binding short peptide motifs through a distinctive safety-belt mechanism. HORMA proteins are now understood as key signaling proteins across kingdoms, serving as infection sensors in a bacterial immune system and playing central roles in eukaryotic cell cycle, genome stability, sexual reproduction, and cellular homeostasis pathways. Here, we describe how HORMA proteins' unique ability to adopt multiple conformational states underlies their functions in these diverse contexts. We also outline how a dedicated AAA+ ATPase regulator, Pch2/TRIP13, manipulates HORMA proteins' conformational states to activate or inactivate signaling in different cellular contexts. The emergence of Pch2/TRIP13 as a lynchpin for HORMA protein action in multiple genome-maintenance pathways accounts for its frequent misregulation in human cancers and highlights TRIP13 as a novel therapeutic target. Expected final online publication date for the Annual Review of Biochemistry, Volume 91 is June 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Transcription factor SOX4 facilitates BMP2-regulated gene expression during invasive trophoblast differentiation

The interplay between growth factors, signaling pathways and transcription factors during placental development is key to controlling trophoblast differentiation. Bone morphogenetic protein 2 (BMP2) has been implicated in trophoblast invasion and spiral artery remodeling during early placental development. However, the molecular mechanisms by which these are accomplished have not been fully elucidated, particularly for transcriptional regulation of key transcription factors. Here, we identified SOX4 as a direct target gene induced by BMP2 in first-trimester placental trophoblasts. Analysis of single-cell RNA-seq data from first-trimester placentas and decidua tissues revealed that SOX4 expression is mainly localized in extravillous trophoblast and decidual stromal cells. Moreover, gain- and loss-of-function approaches demonstrated that SOX4 exerts a pro-invasive role in human trophoblasts, and this effect contributes to BMP2-enhanced trophoblast invasion. Importantly, we found that SOX4 was required for BMP2-induced regulation of a subset of genes associated with cell migration and extracellular matrix organization. We also show that SOX4-dependent regulation of the BMP2 target SERPINE2 occurs via binding of SOX4 to regulatory elements such as enhancers, thereby promoting BMP2-induced trophoblast invasion. In conclusion, these findings uncover a novel mechanism involving SOX4 that shapes the BMP2-regulated transcriptional network during invasive trophoblast development.

Homeostatic control of nuclear-encoded mitochondrial gene expression by the histone variant H2A.Z is essential for neuronal survival

Histone variants are crucial regulators of chromatin structure and gene transcription, yet their functions within the brain remain largely unexplored. Here, we show that the H2A histone variant H2A.Z is essential for neuronal survival. Mice lacking H2A.Z in GABAergic neurons or Purkinje cells (PCs) present with a progressive cerebellar ataxia accompanied by widespread degeneration of PCs. Ablation of H2A.Z in other neuronal subtypes also triggers cell death. H2A.Z binds to the promoters of key nuclear-encoded mitochondrial genes to regulate their expression and promote organelle function. Bolstering mitochondrial activity genetically or by organelle transplant enhances the survival of H2A.Z-ablated neurons. Changes in bioenergetic status alter H2A.Z occupancy at the promoters of nuclear-encoded mitochondrial genes, an adaptive response essential for cell survival. Our results highlight that H2A.Z fulfills a key, conserved role in neuronal survival by acting as a transcriptional rheostat to regulate the expression of genes critical to mitochondrial function.

Meta-analysis of transcriptomic data reveals clusters of consistently deregulated gene and disease ontologies in Down syndrome

Trisomy of human chromosome 21 (HSA21) causes Down syndrome (DS). The trisomy does not simply result in the upregulation of HSA21--encoded genes but also leads to a genome-wide transcriptomic deregulation, which affect differently each tissue and cell type as a result of epigenetic mechanisms and protein-protein interactions. We performed a meta-analysis integrating the differential expression (DE) analyses of all publicly available transcriptomic datasets, both in human and mouse, comparing trisomic and euploid transcriptomes from different sources. We integrated all these data in a "DS network". We found that genome wide deregulation as a consequence of trisomy 21 is not arbitrary, but involves deregulation of specific molecular cascades in which both HSA21 genes and HSA21 interactors are more consistently deregulated compared to other genes. In fact, gene deregulation happens in "clusters", so that groups from 2 to 13 genes are found consistently deregulated. Most of these events of "co-deregulation" involve genes belonging to the same GO category, and genes associated with the same disease class. The most consistent changes are enriched in interferon related categories and neutrophil activation, reinforcing the concept that DS is an inflammatory disease. Our results also suggest that the impact of the trisomy might diverge in each tissue due to the different gene set deregulation, even though the triplicated genes are the same. Our original method to integrate transcriptomic data confirmed not only the importance of known genes, such as SOD1, but also detected new ones that could be extremely useful for generating or confirming hypotheses and supporting new putative therapeutic candidates. We created "metaDEA" an R package that uses our method to integrate every kind of transcriptomic data and therefore could be used with other complex disorders, such as cancer. We also created a user-friendly web application to query Ensembl gene IDs and retrieve all the information of their differential expression across the datasets.

TRIP13 modulates protein deubiquitination and accelerates tumor development and progression of B cell malignancies

Multiple myeloma (MM), a terminally differentiated B cell malignancy, remains difficult to cure. Understanding the molecular mechanisms underlying the progression of MM may identify therapeutic targets and lead to a fundamental shift in treatment of the disease. Deubiquitination, like ubiquitination, is a highly regulated process, implicated in almost every cellular process. Multiple deubiquitinating enzymes (DUBs) have been identified, but their regulation is poorly defined. Here, we determined that TRIP13 increases cellular deubiquitination. Overexpression of TRIP13 in mice and cultured cells resulted in excess cellular deubiquitination by enhancing the association of the DUB USP7 with its substrates. We show that TRIP13 is an oncogenic protein because it accelerates B cell tumor development in transgenic mice. TRIP13-induced resistance to proteasome inhibition can be overcome by a USP7 inhibitor in vitro and in vivo. These findings suggest that TRIP13 expression plays a critical role in B cell lymphoma and MM by regulating deubiquitination of critical oncogenic (NEK2) and tumor suppressor (PTEN, p53) proteins. High TRIP13 identifies a high-risk patient group amenable to adjuvant anti-USP7 therapy.

Feedback control of PLK1 by Apolo1 ensures accurate chromosome segregation

Stable transmission of genetic material during cell division requires accurate chromosome segregation. PLK1 dynamics at kinetochores control establishment of correct kinetochore-microtubule attachments and subsequent silencing of the spindle checkpoint. However, the regulatory mechanism responsible for PLK1 activity in prometaphase has not yet been affirmatively identified. Here we identify Apolo1, which tunes PLK1 activity for accurate kinetochore-microtubule attachments. Apolo1 localizes to kinetochores during early mitosis, and suppression of Apolo1 results in misaligned chromosomes. Using the fluorescence resonance energy transfer (FRET)-based PLK1 activity reporter, we found that Apolo1 sustains PLK1 kinase activity at kinetochores for accurate attachment during prometaphase. Apolo1 is a cognate substrate of PLK1, and the phosphorylation enables PP1γ to inactivate PLK1 by dephosphorylation. Mechanistically, Apolo1 constitutes a bridge between kinase and phosphatase, which governs PLK1 activity in prometaphase. These findings define a previously uncharacterized feedback loop by which Apolo1 provides fine-tuning for PLK1 to guide chromosome segregation in mitosis.

Prognostic and immunological value of ATP6AP1 in breast cancer: implications for SARS-CoV-2

Abnormal ATPase H+ Transporting Accessory Protein 1 (ATP6AP1) expression may promote carcinogenesis. We investigated the association of ATP6AP1 with breast cancer (BC) and COVID-19. The Oncomine, Gene Expression Profiling Interactive Analysis, Human Protein Atlas and Kaplan-Meier plotter databases were used to evaluate the expression and prognostic value of ATP6AP1 in BC. ATP6AP1 was upregulated in BC tissues, and higher ATP6AP1 expression was associated with poorer outcomes. Data from the Tumor Immune Estimation Resource, Tumor-Immune System Interaction Database and Kaplan-Meier plotter indicated that ATP6AP1 expression correlated with immune infiltration, and that its prognostic effects in BC depended on tumor-infiltrating immune cell subtype levels. Multiple databases were used to evaluate the association of ATP6AP1 with clinicopathological factors, assess the mutation and methylation of ATP6AP1, and analyze gene co-expression and enrichment. The ATP6AP1 promoter was hypomethylated in BC tissues and differentially methylated between different disease stages and subtypes. Data from the Gene Expression Omnibus indicated that ATP6AP1 levels in certain cell types were reduced after SARS-CoV-2 infections. Ultimately, higher ATP6AP1 expression was associated with a poorer prognosis and with higher or lower infiltration of particular immune cells in BC. BC patients may be particularly susceptible to SARS-CoV-2 infections, which may alter their prognoses.

Meta-Analysis of Microdissected Breast Tumors Reveals Genes Regulated in the Stroma but Hidden in Bulk Analysis

Transcriptome data provide a valuable resource for the study of cancer molecular mechanisms, but technical biases, sample heterogeneity, and small sample sizes result in poorly reproducible lists of regulated genes. Additionally, the presence of multiple cellular components contributing to cancer development complicates the interpretation of bulk transcriptomic profiles. To address these issues, we collected 48 microarray datasets derived from laser capture microdissected stroma or epithelium in breast tumors and performed a meta-analysis identifying robust lists of differentially expressed genes. This was used to create a database with carefully harmonized metadata that we make freely available to the research community. As predicted, combining the results of multiple datasets improved statistical power. Moreover, the separate analysis of stroma and epithelium allowed the identification of genes with different contributions in each compartment, which would not be detected by bulk analysis due to their distinct regulation in the two compartments. Our method can be profitably used to help in the discovery of biomarkers and the identification of functionally relevant genes in both the stroma and the epithelium. This database was made to be readily accessible through a user-friendly web interface.

A Prognostic DNA Damage Repair Genes Signature and Its Impact on Immune Cell Infiltration in Glioma

Objective

Glioma is the most frequent type of malignant cerebral tumors. DNA damage repair genes (DDRGs) play a crucial role in the development of cancer. In this study, we constructed a DDRGs signature and investigated the potential mechanisms involved in this disease.

Methods

RNA sequence data, microarray data, and corresponding clinical information of gliomas were downloaded from The Cancer Genome Atlas (TCGA), Chinese Glioma Genome Atlas (CGGA), and Gene Expression Omnibus (GEO). Subsequently, we identified candidate genes by differential analysis and Cox regression analysis. The least absolute shrinkage and selection operator Cox regression model was utilized to construct a DDRGs signature using TCGA training dataset. According to this signature, patients with glioma were divided into low- and high-risk groups. The predictive ability of the signature was validated by prognostic analysis, receiver operating characteristic curves, principal component analysis, and stratification analysis in TCGA testing and CGGA verification datasets. CIBERSORT and single-sample gene set enrichment analysis (ssGSEA) were used to evaluate the immune microenvironment of glioma. Moreover, we conducted GSEA to determine the functions and pathways in the low- and high-risk groups. Finally, a nomogram was constructed by combining the signature and other clinical features.

Results

A total of 1,431 samples of glioma (592 from TCGA, 686 from the CGGA, and 153 from the GEO) and 23 samples of normal brain tissue from the GEO were analyzed in this study. There were 51 prognostic differentially expressed DDRGs. Additionally, five DDRGs (CDK4、HMGB2、WEE1、SMC3 and GADD45G) were selected to construct a DDRGs signature for glioma, stratifying patients into low- and high-risk groups. The survival analysis showed that the DDRGs signature could differentiate the outcome of the low- and high-risk groups, showing that high-risk gliomas were associated with shorter overall survival. The immune microenvironment analysis revealed that more immunosuppressive cells, such as tumor associated macrophages and regulatory T cells, were recruited in the high-risk group. GSEA also showed that high-risk glioma was correlated with the immune and extracellular matrix pathways.

Conclusion

The five DDRGs signature and its impact on the infiltration of immunosuppressive cells could precisely predict the prognosis and provide guidance on the treatment of glioma.

Clinicopathological significance of SOX4 and epithelial-mesenchymal transition markers in patients with laryngeal squamous cell carcinoma

OBJECTIVE

Sex-determining region-Y-related high-mobility-group box 4 (SOX4) is associated with the metastasis and prognosis of many cancer types. However, studies on the role of SOX4 in laryngeal squamous cell carcinoma (LSCC) are few, and hence the mechanism is unclear. Epithelial-mesenchymal transition (EMT) allows neoplastic cells to gain the plasticity and motility required for tumor progression and metastasis. This study aimed to analyze the relationship between SOX4 and EMT, and their relationship with clinicopathological factors and related prognosis.

METHODS

Immunohistochemical staining was used to detect the positive expression of SOX4 protein, EMT-related transcription factor protein, and related marker protein in 127 LSCC tissue samples. At the same time, data on various parameters of clinical pathology and postoperative survival were collected.

RESULTS

The positive expression rate of SOX4 and Slug in LSCC was related to pathological differentiation, Lymph node metastasis (LNM), and pathological TNM of a tumor. The expression rates of ZEB1, E-cadherin, N-cadherin, and β-catenin in LSCC correlated with LNM and pTNM. The expression of SOX4, combined expression of SOX4 and ZEB1, and LNM were independent prognostic factors for the total survival time of patients with LSCC.

CONCLUSIONS

In summary, SOX4 was vital in the LSCC EMT process, which might be mediated by transcription factor ZEB1. SOX4 and ZEB1 might serve as potential biomarkers of metastasis and prognosis, as well as promising therapeutic targets of LSCC.

Investigation of the expression levels of CPEB4, APC, TRIP13, EIF2S3, EIF4A1, IFNg, PIK3CA and CTNNB1 genes in different stage colorectal tumors

Background/aim

The aim of the study is to assess expression levels of CPEB4, APC, TRIP13, EIF2S3, EIF4A1, IFNg, PIK3CA and CTNNB1 genes in tumors and peripheral bloods of colorectal cancer patients in stages I–IV.

Materials and methods

The mRNA levels of the genes were determined in tumor tissues and peripheral blood samples of 45 colorectal cancer patients and colon tissues and peripheral blood samples of 5 healthy individuals. Real-time polymerase chain reaction method was used for the analysis.

Results

The mRNA level of the CPEB4 gene was significantly downregulated in colorectal tumor tissues and was upregulated in the peripheral blood of colorectal cancer patients relative to the controls (P < 0.05). APC mRNA level was significantly downregulated in tissues and upregulated in the peripheral blood (P < 0.05). TRIP13 mRNA level was upregulated in peripheral blood and also significantly upregulated in colorectal tumor tissues (P < 0.05). EIF2S3 mRNA level was upregulated in tissues and also significantly upregulated in peripheral blood (P < 0.05). PIK3CA mRNA level was downregulated in tissues and upregulated in peripheral blood. EIF4A1 mRNA level was downregulated in tissues and significantly upregulated in peripheral blood (P < 0.05). CTNNB1 mRNA level was downregulated in tissues and upregulated in peripheral blood. IFNg mRNA level was upregulated in both colorectal cancer tumor tissues and peripheral blood.

Conclusion:

TRIP13 and CPEB4 mRNA up regulation in the peripheral blood of patients with colorectal cancer may be a potential target for early stage diagnosis. In addition to this evaluation, although there is not much study on EIF2S3 and EIF4A1 mRNA changes in cases with colorectal cancer, upregulation in peripheral blood draws attention in our study. These data will shed light on the new comprehensive studies.

SOX4 as biomarker in hepatitis B virus-associated hepatocellular carcinoma

Background: Hepatitis B virus infection is associated with liver disease, including cancers. In this study, we assessed the power of sex-determining region Y (SRY)-related high-mobility group (HMG)-box 4(SOX4) gene to predict the clinical course of hepatocellular carcinoma (HCC).

Methods: To evaluate the differential expression of SOX4 and its diagnostic and prognostic potential in HCC, we analyzed the GSE14520 dataset. Stratified analysis and joint-effect analysis were done using SOX4 and clinical factor. We then designed a nomogram for predicting the clinical course of HCC. Differential SOX4 expression and its correlation with tumor stage as well as its diagnostic and prognostic value were analyzed on the oncomine and GEPIA websites. Gene set enrichment analysis was explored as well as candidate gene ontology and metabolic pathways modulated by in SOX4 HCC.

Results: Our analysis revealed that the level of SOX4 was significantly upregulated in tumor issue (P <0.001). This observation was validated through oncomine dataset and MERAV analysis (all P <0.05). Diagnostic receiver operating characteristic (ROC) analysis of SOX4 suggested it has diagnostic potential in HCC (GSE14520 dataset: P <0.001, area under curve (AUC) = 0.782; Oncomine: (Wurmbach dataset) P = 0.002, AUC = 0.831 and (Mas dataset) P <0.001, AUC = 0.947). In addition, SOX4 exhibited high correlation with overall survival of HBV-associated HCC (adjusted P = 0.004, hazard ratio (HR) (95% confidence interval (CI)) = 2.055 (1.261-3.349) and recurrence-free survival (adjusted P = 0.008, HR (95% CI) = 1.721 (1.151-2.574). These observations which were verified by GEPIA analysis for overall survival (P = 0.007) and recurrence-free survival (P= 0.096). Gene enrichment analysis revealed that affected processes included lymphocyte differentiation, pancreatic endocrine pathways, and insulin signaling pathway. SOX4 prognostic value was evaluated using nomogram analysis for HCC 1, 3, and 5-year, survival.

Conclusion: Differential SOX4 expression presents an avenue of diagnosing and predicting clinical course of HCC. In HCC, SOX4 may affect TP53 metabolic processes, lymphocyte differentiation and the insulin signaling pathway.

Molecular interactions of miR-338 during tumor progression and metastasis

Background

Cancer, as one of the main causes of human deaths, is currently a significant global health challenge. Since the majority of cancer-related deaths are associated with late diagnosis, it is necessary to develop minimally invasive early detection markers to manage and reduce mortality rates. MicroRNAs (miRNAs), as highly conserved non-coding RNAs, target the specific mRNAs which are involved in regulation of various fundamental cellular processes such as cell proliferation, death, and signaling pathways. MiRNAs can also be regulated by long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs). They are highly stable in body fluids and have tumor-specific expression profiles, which suggest their suitability as efficient non-invasive diagnostic and prognostic tumor markers. Aberrant expression of miR-338 has been widely reported in different cancers. It regulates cell proliferation, migration, angiogenesis, and apoptosis in tumor cells.

Main body

In the present review, we have summarized all miR-338 interactions with other non-coding RNAs (ncRNAs) and associated signaling pathways to clarify the role of miR-338 during tumor progression.

Conclusions

It was concluded that miR-338 mainly functions as a tumor suppressor in different cancers. There were also significant associations between miR-338 and other ncRNAs in tumor cells. Moreover, miR-338 has a pivotal role during tumor progression using the regulation of WNT, MAPK, and PI3K/AKT signaling pathways. This review highlights miR-338 as a pivotal ncRNA in biology of tumor cells.

Solid tumours hijack the histone variant network

Cancer is a complex disease characterized by loss of cellular homeostasis through genetic and epigenetic alterations. Emerging evidence highlights a role for histone variants and their dedicated chaperones in cancer initiation and progression. Histone variants are involved in processes as diverse as maintenance of genome integrity, nuclear architecture and cell identity. On a molecular level, histone variants add a layer of complexity to the dynamic regulation of transcription, DNA replication and repair, and mitotic chromosome segregation. Because these functions are critical to ensure normal proliferation and maintenance of cellular fate, cancer cells are defined by their capacity to subvert them. Hijacking histone variants and their chaperones is emerging as a common means to disrupt homeostasis across a wide range of cancers, particularly solid tumours. Here we discuss histone variants and histone chaperones as tumour-promoting or tumour-suppressive players in the pathogenesis of cancer.

Clinicopathological significance of SOX4 and epithelial-mesenchymal transition markers in patients with laryngeal squamous cell carcinoma

BACKGROUND AND AIM

Sex-determining region-Y-related high-mobility-group box 4 (SOX4) is associated with the metastasis and prognosis of many cancer types. However, studies on the role of SOX4 in laryngeal squamous cell carcinoma (LSCC) are few, and hence the mechanism is unclear. Epithelial-mesenchymal transition (EMT) allows neoplastic cells to gain the plasticity and motility required for tumor progression and metastasis. This study aimed to analyze the relationship between SOX4 and EMT, and their relationship with clinicopathological factors and related prognosis.

METHODS

Immunohistochemical staining was used to detect the positive expression of SOX4 protein, EMT-related transcription factor protein, and related marker protein in 127 LSCC tissue samples. At the same time, data on various parameters of clinical pathology and postoperative survival were collected.

RESULTS

The positive expression rate of SOX4 and Slug in LSCC was related to pathological differentiation, lymphatic invasion, and pathological tumor node metastasis (TNM) of a tumor. The expression rates of ZEB1, Twist, E-cadherin, N-cadherin, and β-catenin in LSCC correlated with lymphatic invasion and pathological tumor node metastasis. The expression of SOX4, combined expression of SOX4 and ZEB1, and lymphatic invasion were independent prognostic factors for the total survival time of patients with LSCC.

CONCLUSIONS

In summary, SOX4 was vital in the LSCC EMT process, which might be mediated by transcription factor ZEB1. SOX4 and ZEB1 might serve as potential biomarkers of metastasis and prognosis, as well as promising therapeutic targets of LSCC.

Mass-spectrometry-based proteomic correlates of grade and stage reveal pathways and kinases associated with aggressive human cancers

Proteomic signatures associated with clinical measures of more aggressive cancers could yield molecular clues as to disease drivers. Here, utilizing the Clinical Proteomic Tumor Analysis Consortium (CPTAC) mass-spectrometry-based proteomics datasets, we defined differentially expressed proteins and mRNAs associated with higher grade or higher stage, for each of seven cancer types (breast, colon, lung adenocarcinoma, clear cell renal, ovarian, uterine, and pediatric glioma), representing 794 patients. Widespread differential patterns of total proteins and phosphoproteins involved some common patterns shared between different cancer types. More proteins were associated with higher grade than higher stage. Most proteomic signatures predicted patient survival in independent transcriptomic datasets. The proteomic grade signatures, in particular, involved DNA copy number alterations. Pathways of interest were enriched within the grade-associated proteins across multiple cancer types, including pathways of altered metabolism, Warburg-like effects, and translation factors. Proteomic grade correlations identified protein kinases having functional impact in vitro in uterine endometrial cancer cells, including MAP3K2, MASTL, and TTK. The protein-level grade and stage associations for all proteins profiled-along with corresponding information on phosphorylation, pathways, mRNA expression, and copy alterations-represent a resource for identifying new potential targets. Proteomic analyses are often concordant with corresponding transcriptomic analyses, but with notable exceptions.

CHRNA5 belongs to the secondary estrogen signaling network exhibiting prognostic significance in breast cancer

PURPOSE

Cholinergic signals can be important modulators of cellular signaling in cancer. We recently have shown that knockdown of nicotinic acetylcholine receptor subunit alpha 5, CHRNA5, diminishes the proliferative potential of breast cancer cells. However, modulation of CHRNA5 expression in the context of estrogen signaling and its prognostic implications in breast cancer remained unexplored.

METHODS

Meta-analyses of large breast cancer microarray cohorts were used to evaluate the association of CHRNA5 expression with estrogen (E2) treatment, estrogen receptor (ER) status and patient prognosis. The results were validated through RT-qPCR analyses of multiple E2 treated cell lines, CHRNA5 depleted MCF7 cells and across a breast cancer patient cDNA panel. We also calculated a predicted secondary (PS) score representing direct/indirect induction of gene expression by E2 based on a public dataset (GSE8597). Co-expression analysis was performed using a weighted gene co-expression network analysis (WGCNA) pipeline. Multiple other publicly available datasets such as CCLE, COSMIC and TCGA were also analyzed.

RESULTS

Herein we found that CHRNA5 expression was induced by E2 in a dose- and time-dependent manner in breast cancer cell lines. ER^- breast tumors exhibited higher CHRNA5 expression levels than ER⁺ tumors. Independent meta-analysis for survival outcome revealed that higher CHRNA5 expression was associated with a worse prognosis in untreated breast cancer patients. Furthermore, CHRNA5 and its co-expressed gene network emerged as secondarily induced targets of E2 stimulation. These targets were largely downregulated by exposure to CHRNA5 siRNA in MCF7 cells while the response of primary ESR1 targets was dependent on the direction of the PS-score. Moreover, primary and secondary target genes were uncoupled and clustered distinctly based on multiple public datasets.

CONCLUSION

Our findings strongly associate increased expression of CHRNA5 and its co-expression network with secondary E2 signaling and a worse prognosis in breast cancer.

Inhibition of Cathepsin S Restores TGF-β-induced Epithelial-to-mesenchymal Transition and Tight Junction Turnover in Glioblastoma Cells

Background: Invasive growth is one of the most typical features of aggressive types of malignant cancer, including glioblastoma. Lysosomal cysteine protease-cathepsin S (CTSS), has been reported to be involved in invasive growth and distant metastasis of cancer cells. However, the underlying mechanisms remained elusive. Methods: U87 and U251 human glioblastoma cell lines were applied in this study. Cell migration and invasion ability were measured by wound healing assay and transwell assay. Western blot was employed to detect the expression levels of proteins. Immunofluorescence assays of cells and tissues were used to visualize the localization and expression of proteins. The SPSS software was used for statistical analysis. Results: Our results showed that the high expression of CTSS was link with the grades of glioma tissues. The CTSS inhibitor-Z-FL-COCHO (ZFL), could attenuate TGF-β-induced invasive growth as proven by wound healing and transwell assays. Furthermore, inhibition of CTSS could reverse TGF-β-induced epithelial-to-mesenchymal transition (EMT) and restore TGF-β-triggered tight junction proteins turnover, thus decreasing glioblastoma cell mobility. We also observed that TGF-β could change the morphology of glioblastoma cells, redistribute intermediate-filament, vimentin, which was highly relevant to mesenchymal type cells and enhanced mobility. However, inhibition of CTSS could significantly restore this transformation. Our results proved that PI3K/AKT/mTOR pathway was significantly suppressed in the TGF-β+ZFL (CTSS inhibitor) groups, and AKT activator-SC79, could reverse the anti-invasion effect of CTSS, indicating an important role of PI3K/AKT/mTOR pathway in this process. Conclusion: Z-FL-COCHO (ZFL), a CTSS inhibitor, could reverse TGF-β-induced EMT and change of tight junction proteins via PI3K/AKT/mTOR pathway.

MiR-320b/RAD21 axis affects hepatocellular carcinoma radiosensitivity to ionizing radiation treatment through DNA damage repair signaling

Hepatocellular carcinoma (HCC) is one of the most common malignancies in the world and is associated with high mortality. Ionizing radiation (IR)-based therapy causes DNA damage, exerting a curative effect; however, DNA damage repair signaling pathways lead to HCC resistance to IR-based therapy. RAD21 is a component of the cohesion complex, crucial for chromosome segregation and DNA damage repair, while it is still unclear whether RAD21 is implicated in DNA damage and influences IR sensitivity in HCC. The current research explores the effect and upstream regulatory mechanism of RAD21 on IR sensitivity in HCC. In the present study, RAD21 mRNA and protein expression were increased within HCC tissue samples, particularly within IR-insensitive HCC tissues. The overexpression of RAD21 partially attenuated the roles of IR in HCC by promoting the viability and suppressing the apoptosis of HCC cells. RAD21 overexpression reduced the culture medium 8-hydroxy-2-deoxyguanosine concentration and decreased the protein levels of γH2AX and ATM, suggesting that RAD21 overexpression attenuated IR treatment-induced DNA damage to HCC cells. miR-320b targeted RAD21 3'-UTR to inhibit RAD21 expression. In HCC tissues, particularly in IR-insensitive HCC tissues, miR-320b expression was significantly downregulated. miR-320b inhibition also attenuated IR treatment-induced DNA damage to HCC cells; more importantly, RAD21 silencing significantly attenuated the effects of miR-320b inhibition on IR treatment-induced DNA damage, suggesting that miR-320b plays a role through targeting RAD21. In conclusion, an miR-320b/RAD21 axis modulating HCC sensitivity to IR treatment through acting on IR-induced DNA damage was demonstrated. The miR-320b/RAD21 axis could be a novel therapeutic target for further study of HCC sensitivity to IR treatment.

The molecular basis of gender disparities in smoking lung cancer patients

AIMS

Gender disparities exist in smoking-related lung cancer epidemiology, but the molecular basis has not been explored so far. We aimed at identifying genes with gender-bias expression pattern in smoking lung cancer patients for understanding the molecular basis of gender bias in smokers using meta-analysis of microarray gene expression data.

MATERIALS AND METHODS

Transcriptome of around 1100 samples from 13 studies were used in the meta-analysis to identify 'Lung Cancer genes specific to Female-Smokers' (LCFS) and 'Lung Cancer genes specific to Male-Smokers' (LCMS). The expression profiles of these genes were validated with an independent microarray report and TCGA-RNA-sequencing data. The molecular interactions, pathway, and other functional annotations were portrayed for the key genes identified.

KEY FINDINGS

We identified 1159 gender-biased genes in smoking lung cancer patients. Of these, 400 and 474 genes showed differential expression in cancerous compared to normal lung of women (LCFS) and men (LCMS), respectively. While many up-regulated LCFS were involved in 'immune responses' including T-cell activation, leukocyte cell-cell adhesion, the LCMS were mainly involved in 'positive regulation of gene expression', signaling pathways including RAS, VEGF, insulin-receptor signaling, and 'cell cycle'.

SIGNIFICANCE

The strategic-method identified genes, particularly, SNX20, GIMAP6, MTMR2, FAM171B, IDH1, MOBP, FBXO17, LPXN and WIPF1, which were consistently differentially expressed in at least 4 studies, and in agreement with RNA-Seq data. Exploring their functions could be beneficial to the gender-based diagnosis, prognosis, and treatment of lung cancer in smokers. The current meta-analysis supports existing knowledge of sexual-dimorphism of immune responses in cancer.

TRPM7 Induces Tumorigenesis and Stemness Through Notch Activation in Glioma

We have reported that transient receptor potential melastatin-related 7 (TRPM7) regulates glioma stem cells (GSC) growth and proliferation through Notch, STAT3-ALDH1, and CD133 signaling pathways. In this study, we determined the major contributor(s) to TRPM7 mediated glioma stemness by further deciphering each individual Notch signaling. We first determined whether TRPM7 is an oncotarget in glioblastoma multiforme (GBM) using the Oncomine database. Next, we determined whether TRPM7 silencing by siRNA TRPM7 (siTRPM7) induces cell growth arrest or apoptosis to reduce glioma cell proliferation using cell cycle analysis and annexin V staining assay. We then examined the correlations between the expression of TRPM7 and Notch signaling activity as well as the expression of GSC markers CD133 and ALDH1 in GBM by downregulating TRPM7 through siTRPM7 or upregulating TRPM7 through overexpression of human TRPM7 (M7-wt). To distinguish the different function of channel and kinase domain of TRPM7, we further determined how the α-kinase-dead mutants of TRPM7 (α-kinase domain deleted/M7-DK and K1648R point mutation/M7-KR) affect Notch activities and CD133 and ALDH1 expression. Lastly, we determined the changes in TRPM7-mediated regulation of glioma cell growth/proliferation, cell cycle, and apoptosis by targeting Notch1. The Oncomine data revealed a significant increase in TRPM7 mRNA expression in anaplastic astrocytoma, diffuse astrocytoma, and GBM patients compared to that in normal brain tissues. TRPM7 silencing reduced glioma cell growth by inhibiting cell entry into S and G2/M phases and promoting cell apoptosis. TRPM7 expression in GBM cells was found to be positively correlated with Notch1 signaling activity and CD133 and ALDH1 expression; briefly, downregulation of TRPM7 by siTRPM7 decreased Notch1 signaling whereas upregulation of TRPM7 increased Notch1 signaling. Interestingly, kinase-inactive mutants (M7-DK and M7-KR) resulted in reduced activation of Notch1 signaling and decreased expression of CD133 and ALDH1 compared to that of wtTRPM7. Finally, targeting Notch1 effectively suppressed TRPM7-induced growth and proliferation of glioma cells through cell G1/S arrest and apoptotic induction. TRPM7 is responsible for sustained Notch1 signaling activation, enhanced expression of GSC markers CD133 and ALDH1, and regulation of glioma stemness, which contributes to malignant glioma cell growth and invasion.

Loss of SFRP1 expression is a key progression event in gastrointestinal stromal tumor pathogenesis

The mechanism of high-grade transformation in gastrointestinal stromal tumors (GISTs) remains to be clarified. We aim to discover the key progression events by studying biphasic GISTs. The study group included 101 GISTs. Nineteen of these had been screened from 263 GISTs to represent the early stage of GIST high-grade transformation, characterized by juxtaposed low-grade and high-grade regions in the same tumor (so-called biphasic GISTs). Mutational analyses, fluorescence in situ hybridization (FISH), NanoString analyses, telomere analysis, and gene expression profiling were carried out, followed by in silico analyses, cell line study, and immunohistochemical validation. Using gene expression analysis, downregulation of SFRP1 was revealed to be the main event in GIST high-grade transformation (p = 0.013), accompanied by upregulation of EZH2. In silico analyses revealed that downregulation of SFRP1 was a common feature in GIST progression across several different series. Immunohistochemically, the expression of SFRP1 was validated to be significantly lower in high-grade GISTs (WHO risk group 3a or higher) than in low-grade GISTs (p < 0.001), and attenuation/loss of SFRP1 was associated with GIST tumor progression (p < 0.001). By NanoString and FISH analyses, chromosomal 9/9p loss was the only recurrent large-scale chromosome aberration in biphasic GISTs, with a correlation with SFRP1 downregulation. Subclones containing chromosome 9/9p loss could be appreciated in the low-grade parts of biphasic GISTs. TP53 mutation, RB1 loss, KIT/PDGFRA mutation, and alternative lengthening of telomeres did not play a significant role in GIST high-grade transformation. In conclusion, high-grade transformation of GISTs features SFRP1 downregulation and chromosome 9/9p loss.

MiR-212-3p mediates apoptosis and invasion of esophageal squamous cell carcinoma through inhibition of the Wnt/β-catenin signaling pathway by targeting SOX4

Background

Esophageal squamous cell carcinoma (ESCC) is the predominant subtype of esophageal cancer in East Asia, with approximately half of ESCC cases occurring in China. ESCC poses a serious threat to the quality of life of patients. MicroRNAs (miRNAs) are extremely important in the occurrence and development of ESCC. Current studies have shown that miR-212-3p is expressed at low levels in esophageal adenocarcinoma tumor tissues; however, its function and mechanism in ESCC have not been studied.

Methods

The expression levels of miR-212-3p and SOX4 were detected by quantitative polymerase chain reaction (qPCR) in ESCC tissues and adjacent tissues, and ESCC cell lines from 60 patients. The interaction of miR-212-3p and SOX4 was determined using a dual-luciferase reporter gene system. EC9706 Cells were transfected with miR-212-3p mimic, NC mimic, si-NC and si-SOX4 and miR-212-3p mimic + overexpressing-SOX4, and the effects on miR-212-3p and SOX4 expression were observed, respectively. An MTT assay was carried out to detect the proliferation ability of ESCC. The invasion ability and apoptosis level of the cells were determined by Transwell assay and flow cytometry, respectively. qPCR was used to detect expression of miR-212-3p and SOX4. Western blot was performed to observe the expression of SOX4, Wnt1, β-catenin, c-Myc, and Cyclin D1.

Results

miR-212-3P was observed to be down-regulated in ESCC tissues and cells, while SOX4 expression was up-regulated; the two were negatively correlated. The dual-luciferase reporter gene further confirmed that miR-212-3p targeted SOX4. miR-212-3p overexpression and interference with SOX4 significantly inhibited the proliferation and invasion of ESCC EC970 cells, and promoted apoptosis. Furthermore, the results of Western blot confirmed that miR-212-3p overexpression and interference with SOX4 down-regulated the expression of Wnt1, β-catenin, c-Myc, and Cyclin D1. Meanwhile, SOX4 overexpression reversed the effect of up-regulation of miR-212-3p on EC970 function.

Conclusions

miR-212-3p mediates the apoptosis and invasion of ESCC cells through inhibiting the Wnt/β-catenin signal pathway by targeting SOX4.

S100A2 promotes glycolysis and proliferation via GLUT1 regulation in colorectal cancer

The deregulation of S100A2 has been implicated in the pathogenesis of several types of cancers. However, the molecular mechanisms underlying the protumorigenic capacities of S100A2 have not been fully elucidated. Here, we demonstrated the molecular mechanisms underlying the roles of S100A2 in glycolysis reprogramming and proliferation of colorectal cancer (CRC) cells. The results indicated that S100A2 overexpression raises glucose metabolism and proliferation. Mechanistically, S100A2 activated the PI3K/AKT signaling pathway, upregulated GLUT1 expression, induced glycolytic reprogramming, and consequently increased proliferation. Clinical data showed significantly increased S100A2 levels in CRC tissues and the Oncomine database. In addition, analysis revealed a positive correlation between S100A2 and GLUT1 mRNA expression in CRC tissues. Together, these results demonstrate that the S100A2/GLUT1 axis can promote the progression of CRC by modulating glycolytic reprogramming. Our results further suggest that targeting S100A2 could present a promising therapeutic avenue for the prevention of colorectal cancer progression.

Differential Characteristics of HMGB2 Versus HMGB1 and their Perspectives in Ovary and Prostate Cancer

We have summarized common and differential functions of HMGB1 and HMGB2 proteins with reference to pathological processes, with a special focus on cancer. Currently, several "omic" approaches help us compare the relative expression of these 2 proteins in healthy and cancerous human specimens, as well as in a wide range of cancer-derived cell lines, or in fetal versus adult cells. Molecules that interfere with HMGB1 functions, though through different mechanisms, have been extensively tested as therapeutic agents in animal models in recent years, and their effects are summarized. The review concludes with a discussion on the perspectives of HMGB molecules as targets in prostate and ovarian cancers.

RNAAgeCalc: A multi-tissue transcriptional age calculator

Biological aging reflects decline in physiological functions and is an effective indicator of morbidity and mortality. Numerous epigenetic age calculators are available, however biological aging calculators based on transcription remain scarce. Here, we introduce RNAAgeCalc, a versatile across-tissue and tissue-specific transcriptional age calculator. By performing a meta-analysis of transcriptional age signature across multi-tissues using the GTEx database, we identify 1,616 common age-related genes, as well as tissue-specific age-related genes. Based on these genes, we develop new across-tissue and tissue-specific age predictors. We show that our transcriptional age calculator outperforms other prior age related gene signatures as indicated by the higher correlation with chronological age as well as lower median and median error. Our results also indicate that both racial and tissue differences are associated with transcriptional age. Furthermore, we demonstrate that the transcriptional age acceleration computed from our within-tissue predictor is significantly correlated with mutation burden, mortality risk and cancer stage in several types of cancer from the TCGA database, and offers complementary information to DNA methylation age. RNAAgeCalc is available at http://www.ams.sunysb.edu/~pfkuan/softwares.html#RNAAgeCalc, both as Bioconductor and Python packages, accompanied by a user-friendly interactive Shiny app.

Lysine-specific histone demethylase 1B (LSD2/KDM1B) represses p53 expression to promote proliferation and inhibit apoptosis in colorectal cancer through LSD2-mediated H3K4me2 demethylation

Epigenetic alterations have been reported to play critical roles in the development of colorectal cancer (CRC). However, the biological function of the lysine-specific histone demethylase 1B (LSD2/KDM1B) in CRC is not well understood. Therefore, we investigated the characteristics of LSD2 in CRC. We observed significant upregulation of LSD2 in CRC tissue compared to that in normal colorectal tissue. LSD2 promotes CRC cell proliferation and inhibits cell apoptosis through cell cycle regulation, promoting CRC progression both in vitro and in vivo. We found that LSD2 performs these functions by inhibiting the p53-p21-Rb pathway. Finally, we found that LSD2 directly binds to p53 and represses p53 expression via H3K4me2 demethylation at the p53 promoter. Our results revealed that LSD2 acts as an oncogene by binding and inhibiting p53 activity in CRC. Thus, LSD2 may be a new molecular target for CRC treatment.

SHISA3, an antagonist of the Wnt/β-catenin signaling, is epigenetically silenced and its ectopic expression suppresses growth in breast cancer

Breast cancer (BC) is the foremost cause of cancer related deaths in women globally. Currently there is a scarcity of reliable biomarkers for its early stage diagnosis and theranostics monitoring. Altered DNA methylation patterns leading to the silencing of tumor suppressor genes are considered as an important mechanism underlying tumor development and progression in various cancer types, including BC. Very recently, epigenetic silencing of SHISA3, an antagonist of β-catenin, has been reported in various types of tumor. However, the role of SHISA3 in BC has not been investigated yet. Therefore, we aimed at evaluating the contribution of SHISA3 in BC causation by analyzing its expression and methylation levels in BC cell lines (MDA-MB231, MCF-7 and BT-474) and in 103 paired BC tissue samples. The SHISA3 expression and methylation status was determined by qPCR and methylation specific PCR (MSP) respectively. The role of SHISA3 in BC tumorigenesis was evaluated by proliferation and migration assays after ectopic expression of SHISA3. The association between SHISA3 hypermethylation and clinicopathological parameters of BC patients was also studied. The downregulation of SHISA3 expression was found in three BC cell lines used and in all BC tissue samples. However, SHISA3 promoter region was hypermethylated in 61% (63/103) tumorous tissues in comparison to the 18% of their matched normal tissues. The 5-aza-2’-deoxycytidine treatment restored SHISA3 expression by reversing promoter hypermethylation in both MDA-MB231 and MCF-7 cells. Furthermore, ectopic expression of SHISA3 significantly reduced the proliferation and migration ability of these cells. Taken together, our findings for the first time reveal epigenetic silencing and tumor suppressing role of SHISA3 in BC. Henceforth, this study has identified SHISA3 as potentially powerful target for the development of new therapies against BC, as well as novel diagnostic and therapy response monitoring approaches.

METTL14-mediated N6-methyladenosine modification of SOX4 mRNA inhibits tumor metastasis in colorectal cancer

BACKGROUND

Colorectal cancer (CRC) is one of the leading causes of tumor-related death worldwide, and its main cause of death is distant metastasis. Methyltransferase-like 14(METTL14), a major RNA N6-adenosine methyltransferase, is involved in tumor progression via regulating RNA function. The goal of the study is to uncover the biological function and molecular mechanism of METTL14 in CRC.

METHODS

Quantitative real-time PCR (qRT-PCR), western blot and immunohistochemical (IHC) assays were employed to detect METTL14 and SOX4 in CRC cell lines and tissues. The biological functions of METTL14 were demonstrated using in vitro and in vivo experiments. Chromatin immunoprecipitation (ChIP), Transcrptomic RNA sequencing (RNA-Seq), m6A-RNA immunoprecipitation sequencing (MeRIP-Seq), RNA immunoprecipitation and luciferase reporter assays were used to explore the mechanism of METTL14 action.

RESULTS

METTL14 expression was significantly downregulated in CRC and decreased METTL14 was associated with poor overall survival (OS). Both the univariate and multivariate Cox regression analysis indicated that METTL14 was an independent prognostic factor in CRC. Moreover, lysine-specific histone demethylase 5C(KDM5C)-mediated demethylation of histone H3 lysine 4 tri-methylation(H3K4me3) in the promoter of METTL14 inhibited METTL14 transcription. Functionally, we verified that METTL14 inhibited CRC cells migration, invasion and metastasis through in vitro and in vivo assays, respectively. Furthermore, we identified SRY-related high-mobility-group box 4(SOX4) as a target of METTL14-mediated m6A modification. Knockdown of METTL14 markedly abolished SOX4 mRNA m6A modification and elevated SOX4 mRNA expression. We also revealed that METTL14-mediated SOX4 mRNA degradation relied on the YTHDF2-dependent pathway. Lastly, we demonstrated that METTL14 might inhibit CRC malignant process partly through SOX4-mediated EMT process and PI3K/Akt signals.

CONCLUSIONS

Decreased METTL14 facilitates tumor metastasis in CRC, suggesting that METTL14 might be a potential prognostic biomarker and effective therapeutic target for CRC.

metaRE R Package for Meta-Analysis of Transcriptome Data to Identify the cis-Regulatory Code behind the Transcriptional Reprogramming

At the molecular level, response to an external factor or an internal condition causes reprogramming of temporal and spatial transcription. When an organism undergoes physiological and/or morphological changes, several signaling pathways are activated simultaneously. Examples of such complex reactions are the response to temperature changes, dehydration, various biologically active substances, and others. A significant part of the regulatory ensemble in such complex reactions remains unidentified. We developed metaRE, an R package for the systematic search for cis-regulatory elements enriched in the promoters of the genes significantly changed their transcription in a complex reaction. metaRE mines multiple expression profiling datasets generated to test the same organism’s response and identifies simple and composite cis-regulatory elements systematically associated with differential expression of genes. Here, we showed metaRE performance for the identification of low-temperature-responsive cis-regulatory code in Arabidopsis thaliana and Danio rerio. MetaRE identified potential binding sites for known as well as unknown cold response regulators. A notable part of cis-elements was found in both searches discovering great conservation in low-temperature responses between plants and animals.

Progranulin/EphA2 axis: A novel oncogenic mechanism in bladder cancer

The growth factor progranulin plays a critical role in bladder cancer by modulating tumor cell motility and invasion. Progranulin regulates remodeling of the actin cytoskeleton by interacting with drebrin, an actin binding protein that regulates tumor growth. We previously discovered that progranulin depletion inhibits epithelial-to-mesenchymal transition and markedly reduces in vivo tumor growth. Moreover, progranulin depletion sensitizes urothelial cancer cells to cisplatin treatment, further substantiating a pro-survival function of progranulin. Until recently, the progranulin signaling receptor remained unidentified, precluding a full understanding of progranulin action in tumor cell biology. We recently identified EphA2, a member of a large family of receptor tyrosine-kinases, as the functional receptor for progranulin. However, it is not established whether EphA2 plays an oncogenic role in bladder cancer. Here we demonstrate that progranulin, and not ephrin-A1, the canonical ligand for EphA2, is the predominant EphA2 ligand in bladder cancer. Progranulin evoked Akt- and Erk1/2-mediated EphA2 phosphorylation at Ser897, which could drive bladder tumorigenesis. We discovered that EphA2 depletion severely blunted progranulin-dependent motility and anchorage-independent growth, and sensitized bladder cancer cells to cisplatin treatment. We further defined the mechanisms of progranulin/EphA2-dependent motility by identifying liprin-α1 as a novel progranulin-dependent EphA2 interacting protein and establishing its critical role in cell motility. The discovery of EphA2 as the functional signaling receptor for progranulin and the identification of novel downstream effectors offer a new avenue for understanding the underlying mechanism of progranulin action and may constitute novel clinical and therapeutic targets in bladder cancer.

Translating Cancer Molecular Variability into Personalized Information Using Bulk and Single Cell Approaches

Cancer research is striving toward new frontiers of assigning the correct personalized drug(s) to a given patient. However, extensive tumor heterogeneity poses a major obstacle. Tumors of the same type often respond differently to therapy, due to patient-specific molecular aberrations and/or untargeted tumor subpopulations. It is frequently not possible to determine a priori which patients will respond to a certain therapy or how an efficient patient-specific combined therapy should be designed. Large-scale datasets have been growing at an accelerated pace and various technologies and analytical tools for single cell and bulk level analyses are being developed to extract significant individualized signals from such heterogeneous data. However, personalized therapies that dramatically alter the course of the disease remain scarce, and most tumors still respond poorly to medical care. In this review, the basic concepts of bulk and single cell approaches are discussed, as well as their emerging role in individualized designs of drug therapies, including the advantages and limitations of their applications in personalized medicine.