Emerging roles of the EBF family of transcription factors in tumor suppression

Construction of a Prognostic Model for Hepatocellular Carcinoma Based on Macrophage Polarization-Related Genes

Background

The progression of hepatocellular carcinoma (HCC) is related to macrophage polarization (MP). Our aim was to identify genes associated with MP in HCC patients and develop a prognostic model based on these genes.

Results

We successfully developed a prognostic model consisting of six MP-related genes (SCN4A, EBF3, ADGRB2, HOXD9, CLEC1B, and MSC) to calculate the risk score for each patient. Patients were then classified into high- and low-risk groups based on their median risk score. The performance of the MP-related prognostic model was evaluated using Kaplan-Meier and ROC curves, which yielded favorable results. Additionally, the nomogram demonstrated good clinical effectiveness and displayed consistent survival predictions with actual observations. Gene Set Enrichment Analysis (GSEA) revealed enrichment of pathways related to KRAS signaling downregulation, the G2M checkpoint, and E2F targets in the high-risk group. Conversely, pathways associated with fatty acid metabolism, xenobiotic metabolism, bile acid metabolism, and adipogenesis were enriched in the low-risk group. The risk score positively correlated with the number of invasion-related genes. Immune checkpoint expression differed significantly between the two groups. Patients in the high-risk group exhibited increased sensitivity to mitomycin C, cisplatin, gemcitabine, rapamycin, and paclitaxel, while those in the low-risk group showed heightened sensitivity to doxorubicin. These findings suggest that the high-risk group may have more invasive HCC with greater susceptibility to specific drugs. IHC staining revealed higher expression levels of SCN4A in HCC tissues. Furthermore, experiments conducted on HepG2 cells demonstrated that supernatants from cells with reduced SCN4A expression promoted M2 macrophage polarization marker, CD163 in THP-1 cells. Reduced SCN4A expression induced HCC-related genes, while increased SCN4A expression reduced their expression in HepG2 cells.

Conclusion

The MP-related prognostic model comprising six MPRGs can effectively predict HCC prognosis, infer invasiveness, and guide drug therapy. SCN4A is identified as a suppressor gene in HCC.

Identification of prognostic stemness-related genes in kidney renal papillary cell carcinoma

BACKGROUND

Kidney renal papillary cell carcinoma (KIRP) is the second most prevalent malignant cancer originating from the renal epithelium. Nowadays, cancer stem cells and stemness-related genes (SRGs) are revealed to play important roles in the carcinogenesis and metastasis of various tumors. Consequently, we aim to investigate the underlying mechanisms of SRGs in KIRP.

METHODS

RNA-seq profiles of 141 KIRP samples were downloaded from the TCGA database, based on which we calculated the mRNA expression-based stemness index (mRNAsi). Next, we selected the differentially expressed genes (DEGs) between low- and high-mRNAsi groups. Then, we utilized weighted gene correlation network analysis (WGCNA) and univariate Cox analysis to identify prognostic SRGs. Afterwards, SRGs were included in the multivariate Cox regression analysis to establish a prognostic model. In addition, a regulatory network was constructed by Pearson correlation analysis, incorporating key genes, upstream transcription factors (TFs), and downstream signaling pathways. Finally, we used Connectivity map analysis to identify the potential inhibitors.

RESULTS

In total, 1124 genes were characterized as DEGs between low- and high-RNAsi groups. Based on six prognostic SRGs (CCKBR, GPR50, GDNF, SPOCK3, KC877982.1, and MYO15A), a prediction model was established with an area under curve of 0.861. Furthermore, among the TFs, genes, and signaling pathways that had significant correlations, the CBX2-ASPH-Notch signaling pathway was the most significantly correlated. Finally, resveratrol might be a potential inhibitor for KIRP.

CONCLUSIONS

We suggested that CBX2 could regulate ASPH through activation of the Notch signaling pathway, which might be correlated with the carcinogenesis, development, and unfavorable prognosis of KIRP.

A genome-wide association study identified PTPN2 as a population-specific susceptibility gene locus for primary biliary cholangitis

BACKGROUND AND AIMS

Previous genome-wide association studies (GWAS) have indicated the involvement of shared (population-nonspecific) and nonshared (population-specific) susceptibility genes in the pathogenesis of primary biliary cholangitis (PBC) among European and East-Asian populations. Although a meta-analysis of these distinct populations has recently identified more than 20 novel PBC susceptibility loci, analyses of population-specific genetic architecture are still needed for a more comprehensive search for genetic factors in PBC.

APPROACH AND RESULTS

Protein tyrosine phosphatase nonreceptor type 2 ( PTPN2) was identified as a novel PBC susceptibility gene locus through GWAS and subsequent genome-wide meta-analysis involving 2181 cases and 2699 controls from the Japanese population (GWAS-lead variant: rs8098858, p = 2.6 × 10 -8 ). In silico and in vitro functional analyses indicated that the risk allele of rs2292758, which is a primary functional variant, decreases PTPN2 expression by disrupting Sp1 binding to the PTPN2 promoter in T follicular helper cells and plasmacytoid dendritic cells. Infiltration of PTPN2-positive T-cells and plasmacytoid dendritic cells was confirmed in the portal area of the PBC liver by immunohistochemistry. Furthermore, transcriptomic analysis of PBC-liver samples indicated the presence of a compromised negative feedback loop in vivo between PTPN2 and IFNG in patients carrying the risk allele of rs2292758.

CONCLUSIONS

PTPN2 , a novel susceptibility gene for PBC in the Japanese population, may be involved in the pathogenesis of PBC through an insufficient negative feedback loop caused by the risk allele of rs2292758 in IFN-γ signaling. This suggests that PTPN2 could be a potential molecular target for PBC treatment.

Intratumoral IL-12 delivery via mesenchymal stem cells combined with PD-1 blockade leads to long-term antitumor immunity in a mouse glioblastoma model

BACKGROUND

Although PD-1 blockade is effective for treating several types of cancer, the efficacy of this agent in glioblastoma is largely limited. To overcome non-responders and the immunosuppressive tumor microenvironment, combinational immunotherapeutic strategies with anti-PD-1 need to be considered. Here, we developed IL-12-secreting mesenchymal stem cells (MSC_IL-12) with glioblastoma tropism and evaluated the therapeutic effects of anti-PD-1, MSC_IL-12, and their combination against glioblastoma.

METHODS

Therapeutic responses were evaluated using an immunocompetent mouse orthotopic model. Tumor-infiltrating lymphocytes (TILs) were analyzed using immunofluorescent imaging. Single-cell transcriptome was obtained from mouse brains after treatments.

RESULTS

Anti-PD-1 and MSC_IL-12 showed complete tumor remission in 25.0% (4/16) and 23.1% (3/13) of glioblastoma-implanted mice, respectively, and their combination yielded synergistic antitumor efficacy indicated by 50.0% (6/12) of complete tumor remission. Analyses of TILs revealed that anti-PD-1 increased CD8+ T cells, while MSC_IL-12 led to infiltration of CD4+ T cells and NK cells. Both therapies reduced frequencies of Tregs. All these aspects observed in each monotherapy group were superimposed in the combination group. Notably, no tumor growth was observed upon rechallenge in cured mice, indicating long-term immunity against glioblastoma provoked by the therapies. Single-cell RNA-seq data confirmed these results and revealed that the combined treatment led to immune-favorable tumor microenvironment-CD4+, CD8+ T cells, effector memory T cells, and activated microglia were increased, whereas exhausted T cells, Tregs, and M2 polarized microglia were reduced.

CONCLUSION

Anti-PD-1 and MSC_IL-12 monotherapies show long-term therapeutic responses, and their combination further enhances antitumor efficacy against glioblastoma via inducing immune-favorable tumor microenvironment.

The role of GATA family transcriptional factors in haematological malignancies: A review

GATA transcriptional factors are zinc finger DNA binding proteins that regulate transcription during development and cell differentiation. The 3 important GATA transcription factors GATA1, GATA2 and GATA3 play essential role in the development and maintenance of hematopoietic systems. GATA1 is required for the erythroid and Megakaryocytic commitment during hematopoiesis. GATA2 is crucial for the proliferation and survival of early hematopoietic cells, and is also involved in lineage specific transcriptional regulation as the dynamic partner of GATA1. GATA3 plays an essential role in T lymphoid cell development and immune regulation. As a result, mutations in gene encoding the GATA transcription factor or alteration in the protein expression level or their function have been linked to a variety of human haematological malignancies. This review presents a summary of recent understanding of how the disrupted biological function of GATA may contribute to hematologic diseases.

Targeted DNA Methylation Editing Using an All-in-One System Establishes Paradoxical Activation of EBF3

Cutaneous melanoma is rapidly on the rise globally, surpassing the growth rate of other cancers, with metastasis being the primary cause of death in melanoma patients. Consequently, understanding the mechanisms behind this metastatic process and exploring innovative treatments is of paramount importance. Recent research has shown promise in unravelling the role of epigenetic factors in melanoma progression to metastasis. While DNA hypermethylation at gene promoters typically suppresses gene expression, we have contributed to establishing the newly understood mechanism of paradoxical activation of genes via DNA methylation, where high methylation coincides with increased gene activity. This mechanism challenges the conventional paradigm that promoter methylation solely silences genes, suggesting that, for specific genes, it might actually activate them. Traditionally, altering DNA methylation in vitro has involved using global demethylating agents, which is insufficient for studying the mechanism and testing the direct consequence of gene methylation changes. To investigate promoter hypermethylation and its association with gene activation, we employed a novel approach utilising a CRISPR-SunTag All-in-one system. Here, we focused on editing the DNA methylation of a specific gene promoter segment (EBF3) in melanoma cells using the All-in-one system. Using bisulfite sequencing and qPCR with RNA-Seq, we successfully demonstrated highly effective methylation and demethylation of the EBF3 promoter, with subsequent gene expression changes, to establish and validate the paradoxical role of DNA methylation. Further, our study provides novel insights into the function of the EBF3 gene, which remains largely unknown. Overall, this study challenges the conventional view of methylation as solely a gene-silencing mechanism and demonstrates a potential function of EBF3 in IFN pathway signalling, potentially uncovering new insights into epigenetic drivers of malignancy and metastasis.

Inhibition of epigenetic and cell cycle-related targets in glioblastoma cell lines reveals that onametostat reduces proliferation and viability in both normoxic and hypoxic conditions

The choice of targeted therapies for treatment of glioblastoma patients is currently limited, and most glioblastoma patients die from the disease recurrence. Thus, systematic studies in simplified model systems are required to pinpoint the choice of targets for further exploration in clinical settings. Here, we report screening of 5 compounds targeting epigenetic writers or erasers and 6 compounds targeting cell cycle-regulating protein kinases against 3 glioblastoma cell lines following incubation under normoxic or hypoxic conditions. The viability/proliferation assay indicated that PRMT5 inhibitor onametostat was endowed with high potency under both normoxic and hypoxic conditions in cell lines that are strongly MGMT-positive (T98-G), weakly MGMT-positive (U-251 MG), or MGMT-negative (U-87 MG). In U-251 MG and U-87 MG cells, onametostat also affected the spheroid formation at concentrations lower than the currently used chemotherapeutic drug lomustine. In T98-G cell line, treatment with onametostat led to dramatic changes in the transcriptome profile by inducing the cell cycle arrest, suppressing RNA splicing, and down-regulating several major glioblastoma cell survival pathways. Further validation by immunostaining in three cell lines confirmed that onametostat affects cell cycle and causes reduction in nucleolar protein levels. In this way, inhibition of epigenetic targets might represent a viable strategy for glioblastoma treatment even in the case of decreased chemo- and radiation sensitivity, although further studies in clinically more relevant models are required.

MAPRE3 as an epigenetic target of EZH2 restricts ovarian cancer proliferation in vitro and in vivo

Ovarian cancer (OC) is a lethal gynecologic cancer and the common cause of death within women worldwide. The polycomb group protein enhancer of zeste homolog 2 (EZH2) is a histone methyltransferase highly expressed in various tumors, including OC. However, the mechanistic basis of EZH2 oncogenic activity in OC remain incompletely understood. Bioinformatics analysis showed that the expression of MAPRE3 was lower in OC tissues than in normal tissues, and was positively correlated with the overall survival. MAPRE3 overexpression decreased cell growth, inducing cell cycle arrest and apoptosis in OC cells, whereas MAPRE3 silencing promoted proliferation and accelerated cell cycle progression of OC cells. The in vivo study validated that overexpression of MAPRE3 impeded tumor formation and growth of OC xenografts in nude mice. In addition, knockdown of EZH2 in OC cells downregulated H3K27me3 expression and increased MAPRE3 expression. Inhibiting EZH2 in OC cells reduced the enrichment of H3K27me3 on the promoter of MAPRE3. Furthermore, MAPRE3 silencing significantly reversed changes in the expression of cell cycle and apoptosis-related markers and cell growth mediated by EZH2 knockdown in OC cells. MAPRE3 functions as a suppressor of OC and is epigenetic repressed by EZH2, suggesting a potential therapeutic strategy for OC by targeting EZH2/MAPRE3 axis.

EBF1, PAX5, and MYC: regulation on B cell development and association with hematologic neoplasms

During lymphocyte development, a diverse repertoire of lymphocyte antigen receptors is produced to battle against pathogens, which is the basis of adaptive immunity. The diversity of the lymphocyte antigen receptors arises primarily from recombination-activated gene (RAG) protein-mediated V(D)J rearrangement in early lymphocytes. Furthermore, transcription factors (TFs), such as early B cell factor 1 (EBF1), paired box gene 5 (PAX5), and proto-oncogene myelocytomatosis oncogene (MYC), play critical roles in regulating recombination and maintaining normal B cell development. Therefore, the aberrant expression of these TFs may lead to hematologic neoplasms.

Analysis of the genomic landscape of primary central nervous system lymphoma using whole-genome sequencing in Chinese patients

Primary central nervous system lymphoma (PCNSL) is an uncommon non-Hodgkin's lymphoma with poor prognosis. This study aimed to depict the genetic landscape of Chinese PCNSLs. Whole-genome sequencing was performed on 68 newly diagnosed Chinese PCNSL samples, whose genomic characteristics and clinicopathologic features were also analyzed. Structural variations were identified in all patients with a mean of 349, which did not significantly influence prognosis. Copy loss occurred in all samples, while gains were detected in 77.9% of the samples. The high level of copy number variations was significantly associated with poor progression-free survival (PFS) and overall survival (OS). A total of 263 genes mutated in coding regions were identified, including 6 newly discovered genes (ROBO2, KMT2C, CXCR4, MYOM2, BCLAF1, and NRXN3) detected in ⩾ 10% of the cases. CD79B mutation was significantly associated with lower PFS, TMSB4X mutation and high expression of TMSB4X protein was associated with lower OS. A prognostic risk scoring system was also established for PCNSL, which included Karnofsky performance status and six mutated genes (BRD4, EBF1, BTG1, CCND3, STAG2, and TMSB4X). Collectively, this study comprehensively reveals the genomic landscape of newly diagnosed Chinese PCNSLs, thereby enriching the present understanding of the genetic mechanisms of PCNSL.

Neurologic, Neuropsychologic, and Neuroradiologic Features ofEBF3-Related Syndrome

Background and Objectives

Heterozygous mutations or deletions of theEBF3gene are known to cause a syndrome characterized by intellectual disability, neurodevelopmental disorders, facial dysmorphisms, hypotonia, and ataxia; the latter is quite common despite in most patients brain MRI is reported to be normal. Despite the predominant neurologic involvement ofEBF3-related syndrome, a systematic definition of neurologic, cognitive/behavioral, and neuroradiologic features is lacking.

Methods

We report on 6 patients (2 females and 4 males, age range 2–12 years), of whom 4 carrying a heterozygous point mutation of theEBF3gene and 2 with 10q26 deletion encompassing the gene, diagnosed at Carlo Besta Neurologic Institute of Milan, Italy. Clinical evaluation was performed by a pediatric neurologist and pediatric dysmorphologist; ataxia severity was rated by Scale for the Assessment and Rating of Ataxia (SARA); brain MRIs were reviewed by expert neuroradiologists; general quotient levels were obtained through standardized Griffiths Mental Development Scales. Patients carrying a 10q26.3 deletion were diagnosed by array-CGH, whereasEBF3variants were detected by whole exome sequencing.

Results

Phenotype was consistent in all patients, but with wide variability in severity. Developmental milestones were invariably delayed and resulted in an extremely variable cognitive impairment. All patients showed ataxic signs, as confirmed by SARA scores, often associated with hypotonia. Brain MRI revealed in all children a cerebellar malformation with vermis hypoplasia and a peculiar foliation anomaly characterized by a radial disposition of cerebellar folia (dandelion sign). Neurophysiologic examinations were unremarkable.

Discussion

EBF3-related syndrome has been so far described as a neurodevelopmental condition with dysmorphic traits, with limited emphasis on the neurologic features; we highlight the predominant neurologic involvement of these patients, which can be explained at least in part by the underlying cerebellar malformation. We therefore propose thatEBF3-related syndrome should be classified and treated as a congenital, nonprogressive ataxia.

Identifying molecular and functional similarities and differences between human primary cardiac valve interstitial cells and ventricular fibroblasts

Introduction: Fibroblasts are mesenchymal cells that predominantly produce and maintain the extracellular matrix (ECM) and are critical mediators of injury response. In the heart, valve interstitial cells (VICs) are a population of fibroblasts responsible for maintaining the structure and function of heart valves. These cells are regionally distinct from myocardial fibroblasts, including left ventricular cardiac fibroblasts (LVCFBs), which are located in the myocardium in close vicinity to cardiomyocytes. Here, we hypothesize these subpopulations of fibroblasts are transcriptionally and functionally distinct. Methods: To compare these fibroblast subtypes, we collected patient-matched samples of human primary VICs and LVCFBs and performed bulk RNA sequencing, extracellular matrix profiling, and functional contraction and calcification assays. Results: Here, we identified combined expression of SUSD2 on a protein-level, and MEOX2, EBF2 and RHOU at a transcript-level to be differentially expressed in VICs compared to LVCFBs and demonstrated that expression of these genes can be used to distinguish between the two subpopulations. We found both VICs and LVCFBs expressed similar activation and contraction potential in vitro, but VICs showed an increase in ALP activity when activated and higher expression in matricellular proteins, including cartilage oligomeric protein and alpha 2-Heremans-Schmid glycoprotein, both of which are reported to be linked to calcification, compared to LVCFBs. Conclusion: These comparative transcriptomic, proteomic, and functional studies shed novel insight into the similarities and differences between valve interstitial cells and left ventricular cardiac fibroblasts and will aid in understanding region-specific cardiac pathologies, distinguishing between primary subpopulations of fibroblasts, and generating region-specific stem-cell derived cardiac fibroblasts.

Identification and verification of eight cancer-associated fibroblasts related genes as a prognostic signature for head and neck squamous cell carcinoma

Cancer-associated fibroblasts (CAFs) can exert their immunosuppressive effects by secreting various effectors that are involved in the regulation of tumor-infiltrating immune cells as well as other immune components in the tumor immune microenvironment (TIME), thereby promoting tumorigenesis, progression, metastasis, and drug resistance. Although a large number of studies suggest that CAFs play a key regulatory role in the development of head and neck squamous cell carcinoma (HNSCC), there are limited studies on the relevance of CAFs to the prognosis of HNSCC. In this study, we identified a prognostic signature containing eight CAF-related genes for HNSCC by univariate Cox analysis, lasso regression, stepwise regression, and multivariate Cox analysis. Our validation in primary cultures of CAFs from human HNSCC and four human HNSCC cell lines confirmed that these eight genes are indeed characteristic markers of CAFs. Immune cell infiltration differences analysis between high-risk and low-risk groups according to the eight CAF-related genes signature hinted at CAFs regulatory roles in the TIME, further revealing its potential role on prognosis. The signature of the eight CAF-related genes was validated in different independent validation cohorts and all showed that it was a valid marker for prognosis. The significantly higher overall survival (OS) in the low-risk group compared to the high-risk group was confirmed by Kaplan-Meier (K-M) analysis, suggesting that the signature of CAF-related genes can be used as a non-invasive predictive tool for HNSCC prognosis. The low-risk group had significantly higher levels of tumor-killing immune cell infiltration, as confirmed by CIBERSORT analysis, such as CD8⁺ T cells, follicular helper T cells, and Dendritic cells (DCs) in the low-risk group. In contrast, the level of infiltration of pro-tumor cells such as M0 macrophages and activated Mast cells (MCs) was lower. It is crucial to delve into the complex mechanisms between CAFs and immune cells to find potential regulatory targets and may provide new evidence for subsequently targeted immunotherapy. These results suggest that the signature of the eight CAF-related genes is a powerful indicator for the assessment of the TIME of HNSCC. It may provide a new and reliable potential indicator for clinicians to predict the prognosis of HNSCC, which may be used to guide treatment and clinical decision-making in HNSCC patients. Meanwhile, CAF-related genes are expected to become tumor biomarkers and effective targets for HNSCC.

ZNF521/EBF1 axis regulates AKR1B1 to promote the proliferation, migration, and invasion of gastric cancer cells

Although the incidence and death rates of gastric cancer (GC) are decreasing, approximately one million new cases and 800,000 GC-related deaths were reported worldwide in 2018. Currently, the oncogenesis of GC remains unclear, and the demand for novel treatment options are unmet. Here, we explored the role of aldo-keto reductase family 1 member B (AKR1B1) in the progression of GC. The proliferation, migration, and invasion of GC cells were evaluated by CCK-8 assay, wound healing assay, and transwell assay, respectively. The interaction between EBF transcription factor 1 (EBF1) and the promoter region of AKR1B1 was determined by luciferase reporter assay and chromatin immunoprecipitation (ChIP). Gene expression levels were measured by quantitative real-time polymerase chain reaction (qRT-PCR) and Western blotting assay. The expression of AKR1B1 was elevated in GC cells, resulting in increased cell proliferation, migration, and invasion. Meanwhile, EBF1 was a negative regulator of AKR1B1; its overexpression suppressed AKR1B1 expression and GC progression. Furthermore, knockdown of ZNF521 induced EBF1 expression, thus suppressing AKR1B1 expression and resulting in attenuated GC growth and invasiveness. Notably, knockdown of ZNF521 attenuated GC progression and was rescued by overexpression of AKR1B1. Our current study revealed a novel ZNF521/EBF1/AKR1B1 axis in GC and elaborated its important role in promoting GC progression, providing potential therapeutic targets for anti-GC treatments.

RETRACTED: hsa-miR-23a~27a~24-2 cluster members inhibit aggressiveness of breast cancer cells by commonly targeting NCOA1, NLK and RAP1B

This article has been retracted: please see Elsevier Policy on Article Withdrawal (http://www.elsevier.com/locate/withdrawalpolicy). This article has been retracted at the request of the Editor-in-Chief. The corresponding author notified the journal of three examples of image duplication within the published article (two in Figure 3D and one in Figure 4A), and requested a corrigendum. As per journal policy when considering corrigendum requests, the journal requested the authors to provide source data relating to these affected figures. The editorial team noticed 12 additional suspected image duplications within the supplied source data and the corresponding author was informed. Upon submission of revised source data, the editorial team noticed two new suspected image duplications. The editorial team have concerns about the provenance of the data and therefore the Editor-in-Chief decided to retract the article.

Early B cell factor 4 modulates FAS-mediated apoptosis and promotes cytotoxic function in human immune cells

Apoptosis is a genetically regulated program of cell death that plays a key role in immune disease processes. We identified EBF4, a little-studied member of the early B cell factor (EBF) family of transcription factors, in a whole-genome CRISPR screen for regulators of Fas/APO-1/CD95-mediated T cell death. Loss of EBF4 increases the half-life of the c-FLIP protein, and its presence in the Fas signaling complex impairs caspase-8 cleavage and apoptosis. Transcriptome analysis revealed that EBF4 regulates molecules such as TBX21, EOMES, granzyme, and perforin that are important for human natural killer (NK) and CD8⁺ T cell functions. Proximity-dependent biotin identification (Bio-ID) mass spectrometry analyses showed EBF4 binding to STAT3, STAT5, and MAP kinase 3 and a strong pathway relationship to interleukin-2 regulated genes, which are known to govern cytotoxicity pathways. Chromatin immunoprecipitation and DNA sequencing analysis defined a canonical EBF4 binding motif, 5'-CCCNNGG/AG-3', closely related to the EBF1 binding site; using a luciferase-based reporter, we found a dose-dependent transcriptional response of this motif to EBF4. We also conducted assay for transposase-accessible chromatin sequencing in EBF4-overexpressing cells and found increased chromatin accessibility upstream of granzyme and perforin and in topologically associated domains in human lymphocytes. Finally, we discovered that the EBF4 has basal expression in human but not mouse NK cells and CD8⁺ T cells and vanishes following activating stimulation. Together, our data reveal key features of a previously unknown transcriptional regulator of human cytotoxic immune function.

Identifying General Tumor and Specific Lung Cancer Biomarkers by Transcriptomic Analysis

The bioinformatic pipeline previously developed in our research laboratory is used to identify potential general and specific deregulated tumor genes and transcription factors related to the establishment and progression of tumoral diseases, now comparing lung cancer with other two types of cancer. Twenty microarray datasets were selected and analyzed separately to identify hub differentiated expressed genes and compared to identify all the deregulated genes and transcription factors in common between the three types of cancer and those unique to lung cancer. The winning DEGs analysis allowed to identify an important number of TFs deregulated in the majority of microarray datasets, which can become key biomarkers of general tumors and specific to lung cancer. A coexpression network was constructed for every dataset with all deregulated genes associated with lung cancer, according to DAVID’s tool enrichment analysis, and transcription factors capable of regulating them, according to oPOSSUM´s tool. Several genes and transcription factors are coexpressed in the networks, suggesting that they could be related to the establishment or progression of the tumoral pathology in any tissue and specifically in the lung. The comparison of the coexpression networks of lung cancer and other types of cancer allowed the identification of common connectivity patterns with deregulated genes and transcription factors correlated to important tumoral processes and signaling pathways that have not been studied yet to experimentally validate their role in lung cancer. The Kaplan–Meier estimator determined the association of thirteen deregulated top winning transcription factors with the survival of lung cancer patients. The coregulatory analysis identified two top winning transcription factors networks related to the regulatory control of gene expression in lung and breast cancer. Our transcriptomic analysis suggests that cancer has an important coregulatory network of transcription factors related to the acquisition of the hallmarks of cancer. Moreover, lung cancer has a group of genes and transcription factors unique to pulmonary tissue that are coexpressed during tumorigenesis and must be studied experimentally to fully understand their role in the pathogenesis within its very complex transcriptomic scenario. Therefore, the downstream bioinformatic analysis developed was able to identify a coregulatory metafirm of cancer in general and specific to lung cancer taking into account the great heterogeneity of the tumoral process at cellular and population levels.

Tracing the origins of glioblastoma by investigating the role of gliogenic and related neurogenic genes/signaling pathways in GBM development: a systematic review

Background

Glioblastoma is one of the most aggressive tumors. The etiology and the factors determining its onset are not yet entirely known. This study investigates the origins of GBM, and for this purpose, it focuses primarily on developmental gliogenic processes. It also focuses on the impact of the related neurogenic developmental processes in glioblastoma oncogenesis. It also addresses why glial cells are at more risk of tumor development compared to neurons.

Methods

Databases including PubMed, MEDLINE, and Google Scholar were searched for published articles without any date restrictions, involving glioblastoma, gliogenesis, neurogenesis, stemness, neural stem cells, gliogenic signaling and pathways, neurogenic signaling and pathways, and astrocytogenic genes.

Results

The origin of GBM is dependent on dysregulation in multiple genes and pathways that accumulatively converge the cells towards oncogenesis. There are multiple layers of steps in glioblastoma oncogenesis including the failure of cell fate-specific genes to keep the cells differentiated in their specific cell types such as p300, BMP, HOPX, and NRSF/REST. There are genes and signaling pathways that are involved in differentiation and also contribute to GBM such as FGFR3, JAK-STAT, and hey1. The genes that contribute to differentiation processes but also contribute to stemness in GBM include notch, Sox9, Sox4, c-myc gene overrides p300, and then GFAP, leading to upregulation of nestin, SHH, NF-κB, and others. GBM mutations pathologically impact the cell circuitry such as the interaction between Sox2 and JAK-STAT pathway, resulting in GBM development and progression.

Conclusion

Glioblastoma originates when the gene expression of key gliogenic genes and signaling pathways become dysregulated. This study identifies key gliogenic genes having the ability to control oncogenesis in glioblastoma cells, including p300, BMP, PAX6, HOPX, NRSF/REST, LIF, and TGF beta. It also identifies key neurogenic genes having the ability to control oncogenesis including PAX6, neurogenins including Ngn1, NeuroD1, NeuroD4, Numb, NKX6-1 Ebf, Myt1, and ASCL1. This study also postulates how aging contributes to the onset of glioblastoma by dysregulating the gene expression of NF-κB, REST/NRSF, ERK, AKT, EGFR, and others.

Exploration of Extracellular Vesicle miRNAs, Targeted mRNAs and Pathways in Prostate Cancer: Relation to Disease Status and Progression

Simple Summary

Prostate cancer lacks non-invasive specific biomarkers for aggressive disease. Urinary extracellular vesicles (uEV) could provide such markers; however, due to technical challenges, little is known regarding the pathogenesis pathways reflected in uEV. We performed a miRNA, target mRNA and pathway study focused on uEV, exploring the differences between cancer (1) status groups (Gleason score) and (2) progression groups. The uEV provided a surprisingly comprehensive presentation of differentially expressed miRNAs, target mRNAs and pathogenesis pathways. The miRNAs associated with prostate cancer status or progression were mostly unique, but still targeted overlapping sets of signalling, resistance, hormonal and immune pathways. Interestingly, mRNA targets of the key miRNAs (miR-892a, miR-223-3p, miR-146a-5p) were widely expressed in both uEV and plasma EV from PCa patients. The study thus suggests that uEV carry a vast presentation of PCa status and progression-linked RNAs that are worth further exploration in large personalized medicine trials.

Abstract

Background: Prostate cancer (PCa) lacks non-invasive specific biomarkers for aggressive disease. We studied the potential of urinary extracellular vesicles (uEV) as a liquid PCa biopsy by focusing on the micro RNA (miRNA) cargo, target messenger RNA (mRNA) and pathway analysis. Methods: We subjected uEV samples from 31 PCa patients (pre-prostatectomy) to miRNA sequencing and matched uEV and plasma EV (pEV) from three PCa patients to mRNA sequencing. EV quality control was performed by electron microscopy, Western blotting and particle and RNA analysis. We compared miRNA expression based on PCa status (Gleason Score) and progression (post-prostatectomy follow-up) and confirmed selected miRNAs by quantitative PCR. Expression of target mRNAs was mapped in matched EV. Results: Quality control showed typical small uEV, pEV, RNA and EV-protein marker enriched samples. Comparisons between PCa groups revealed mostly unique differentially expressed miRNAs. However, they targeted comprehensive and largely overlapping sets of cancer and progression-associated signalling, resistance, hormonal and immune pathways. Quantitative PCR confirmed changes in miR-892a (Gleason Score 7 vs. ≥8), miR-223-3p (progression vs. no progression) and miR-146a-5p (both comparisons). Their target mRNAs were expressed widely in PCa EV. Conclusions: PCa status and progression-linked RNAs in uEV are worth exploration in large personalized medicine trials.

A genome-wide association study of mammographic texture variation

Background

Breast parenchymal texture features, including grayscale variation (V), capture the patterns of texture variation on a mammogram and are associated with breast cancer risk, independent of mammographic density (MD). However, our knowledge on the genetic basis of these texture features is limited.

Methods

We conducted a genome-wide association study of V in 7040 European-ancestry women. V assessments were generated from digitized film mammograms. We used linear regression to test the single-nucleotide polymorphism (SNP)-phenotype associations adjusting for age, body mass index (BMI), MD phenotypes, and the top four genetic principal components. We further calculated genetic correlations and performed SNP-set tests of V with MD, breast cancer risk, and other breast cancer risk factors.

Results

We identified three genome-wide significant loci associated with V: rs138141444 (6q24.1) in ECT2L, rs79670367 (8q24.22) in LINC01591, and rs113174754 (12q22) near PGAM1P5. 6q24.1 and 8q24.22 have not previously been associated with MD phenotypes or breast cancer risk, while 12q22 is a known locus for both MD and breast cancer risk. Among known MD and breast cancer risk SNPs, we identified four variants that were associated with V at the Bonferroni-corrected thresholds accounting for the number of SNPs tested: rs335189 (5q23.2) in PRDM6, rs13256025 (8p21.2) in EBF2, rs11836164 (12p12.1) near SSPN, and rs17817449 (16q12.2) in FTO. We observed significant genetic correlations between V and mammographic dense area (r_g = 0.79, P = 5.91 × 10⁻⁵), percent density (r_g = 0.73, P = 1.00 × 10⁻⁴), and adult BMI (r_g =  − 0.36, P = 3.88 × 10⁻⁷). Additional significant relationships were observed for non-dense area (z =  − 4.14, P = 3.42 × 10⁻⁵), estrogen receptor-positive breast cancer (z = 3.41, P = 6.41 × 10⁻⁴), and childhood body fatness (z =  − 4.91, P = 9.05 × 10⁻⁷) from the SNP-set tests.

Conclusions

These findings provide new insights into the genetic basis of mammographic texture variation and their associations with MD, breast cancer risk, and other breast cancer risk factors.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13058-022-01570-8.

EBF3 reactivation by inhibiting the EGR1/EZH2/HDAC9 complex promotes metastasis via transcriptionally enhancing vimentin in nasopharyngeal carcinoma

Metastasis is the major reason for treatment failure and accounts for cancer-related death in patients with nasopharyngeal carcinoma. However, the genetic alterations and molecular mechanisms that cause nasopharyngeal carcinoma metastasis are elusive. Herein, we performed RNA sequencing in patients with or without metastasis, and found that the early B-cell factor 3 (EBF3) was significantly elevated in the samples with metastasis. Mechanistically, EBF3 promoted metastasis by directly combining with the promoter of Vimentin and transcriptionally upregulating it. In addition, EBF3 was epigenetically silenced by EGR1/EZH2/HDAC9 complexes via sustaining the high level of H3K27-Me3 at its promoter. Clinically, there was a positive correlation between EBF3 and Vimentin in nasopharyngeal carcinoma tissues. Moreover, high expression of EBF3 or Vimentin was correlated with poor overall survival, while the combination of high EBF3 and Vimentin expression was associated with more significant poor prognosis. Therefore, specific agents targeting EBF3 or stabilizing the EGR1/EZH2/HDAC9 complex could be novel therapeutic strategies for cancer metastasis.

EBF1 is expressed in pericytes and contributes to pericyte cell commitment

Early B-cell factor-1 (EBF1) is a transcription factor with an important role in cell lineage specification and commitment during the early stage of cell maturation. Originally described during B-cell maturation, EBF1 was subsequently identified as a crucial molecule for proper cell fate commitment of mesenchymal stem cells into adipocytes, osteoblasts and muscle cells. In vessels, EBF1 expression and function have never been documented. Our data indicate that EBF1 is highly expressed in peri-endothelial cells in both tumor vessels and in physiological conditions. Immunohistochemistry, quantitative reverse transcription polymerase chain reaction (RT-qPCR) and fluorescence-activated cell sorting (FACS) analysis suggest that EBF1-expressing peri-endothelial cells represent bona fide pericytes and selectively express well-recognized markers employed in the identification of the pericyte phenotype (SMA, PDGFRβ, CD146, NG2). This observation was also confirmed in vitro in human placenta-derived pericytes and in human brain vascular pericytes (HBVP). Of note, in accord with the key role of EBF1 in the cell lineage commitment of mesenchymal stem cells, EBF1-silenced HBVP cells showed a significant reduction in PDGFRβ and CD146, but not CD90, a marker mostly associated with a prominent mesenchymal phenotype. Moreover, the expression levels of VEGF, angiopoietin-1, NG2 and TGF-β, cytokines produced by pericytes during angiogenesis and linked to their differentiation and activation, were also significantly reduced. Overall, the data suggest a functional role of EBF1 in the cell fate commitment toward the pericyte phenotype.

VNTR polymorphism in the breakpoint region of ABL1 and susceptibility to bladder cancer

BACKGROUND

ABL1 is primarily known as a leukemia-related oncogene due to translocation, but about 2.2% of ABL1 mutations have been identified in bladder cancer, and high expression in solid cancer has also been detected.

METHODS

Here, we used the NCBI database, UCSC genome browser gateway and Tandem repeat finder program to investigate the structural characterization of the ABL1 breakpoint region and to identify the variable number of tandem repeats (VNTR). To investigate the relationship between ABL1-MS1 and bladder cancer, a case-controlled study was conducted in 207 controls and 197 bladder cancer patients. We also examined the level of transcription of the reporter gene driven by the ABL1 promoter to determine if the VNTR region affects gene expression.

RESULTS

In our study, one VNTR was identified in the breakpoint region, the intron 1 region of ABL1, and was named ABL1-MS1. In the control group, only two common alleles (TR13, TR15) were detected, but an additional two rare alleles (TR14, TR16) were detected in bladder cancer. A statistically significant association was identified between the rare ABL1-MS1 allele and bladder cancer risk: P = 0.013. Investigating the level of transcription of the reporter gene driven by the ABL1 promoter, VNTR showed inhibition of ABL1 expression in non-cancer cells 293 T, but not in bladder cancer cells. In addition, ABL1-MS1 was accurately passed on to offspring according to Mendelian inheritance through meiosis.

CONCLUSIONS

Therefore, the ABL1-MS1 region can affect ABL1 expression of bladder cancer. This study provides that ABL1-MS1 can be used as a DNA fingerprinting marker. In addition, rare allele detection can predict susceptibility to bladder cancer.

Neuroimaging Findings in Patients with EBF3 Mutations: Report of Two Cases

Early B cell factor 3 (EBF3) is a transcription factor involved in brain development. Heterozygous, loss-of-function mutations in EBF3 have been reported in an autosomal dominant neurodevelopmental syndrome characterized by hypotonia, ataxia, and developmental delay (sometimes described as "HADD"s). We report 2 unrelated cases with novel de novo EBF3 mutations: c.455G>T (p.Arg152Leu) and c.962dup (p.Tyr321*) to expand the genotype/phenotype correlations of this disorder; clinical, neuropsychological, and MRI studies were used to define the phenotype. IQ was in the normal range and diffusion tensor imaging revealed asymmetric alterations of the longitudinal fasciculus in both cases. Our results demonstrate that EBF3 mutations can underlie neurodevelopmental disorders without intellectual disability. Long tract abnormalities have not been previously recognized and suggest that they may be an unrecognized and characteristic feature in this syndrome.

Promoter hypermethylation of early B cell factor 1 (EBF1) is associated with cholangiocarcinoma progression

DNA hypermethylation in a promoter region causes gene silencing via epigenetic changes. We have previously reported that early B cell factor 1 (EBF1) was down-regulated in cholangiocarcinoma (CCA) tissues and related to tumor progression. Thus, we hypothesized that the DNA hypermethylation of EBF1 promoter would suppress EBF1 expression in CCA and induce its progression. In this study, the DNA methylation status of EBF1 and mRNA expression levels were analyzed in CCA and normal bile duct (NBD) tissues using a publicly available database of genome-wide association data. The results showed that the DNA methylation of EBF1 promoter region was significantly increased in CCA tissues compared with those of NBD. The degree of methylation was negatively correlated with EBF1 mRNA expression levels. Using methylation-specific PCR technique, the DNA methylation rates of EBF1 promoter region were investigated in CCA tissues (n=72). CCA patients with high methylation rates of EBF1 promoter region in the tumor tissues (54/72) had a poor prognosis. Higher methylation rates of EBF1 promoter region have shown in all CCA cell lines than that of an immortal cholangiocyte cell line (MMNK1). Upon treatment with the DNA methyltransferase inhibitor 5-Aza-dC, increased EBF1 expression levels and reduced DNA methylation rates were observed in CCA cells. Moreover, restoration of EBF1 expression in CCA cells led to inhibition of cell growth, migration and invasion. In addition, RNA sequencing analysis suggested that EBF1 is involved in suppression of numerous pathways in cancer. Taken together, DNA hypermethylation in the EBF1 promoter region suppresses EBF1 expression and induces CCA progression with aggressive clinical outcomes.

Genomic Space of MGMT in Human Glioma Revisited: Novel Motifs, Regulatory RNAs, NRF1, 2, and CTCF Involvement in Gene Expression

Background: The molecular regulation of increased MGMT expression in human brain tumors, the associated regulatory elements, and linkages of these to its epigenetic silencing are not understood. Because the heightened expression or non-expression of MGMT plays a pivotal role in glioma therapeutics, we applied bioinformatics and experimental tools to identify the regulatory elements in the MGMT and neighboring EBF3 gene loci. Results: Extensive genome database analyses showed that the MGMT genomic space was rich in and harbored many undescribed RNA regulatory sequences and recognition motifs. We extended the MGMT’s exon-1 promoter to 2019 bp to include five overlapping alternate promoters. Consensus sequences in the revised promoter for (a) the transcriptional factors CTCF, NRF1/NRF2, GAF, (b) the genetic switch MYC/MAX/MAD, and (c) two well-defined p53 response elements in MGMT intron-1, were identified. A putative protein-coding or non-coding RNA sequence was located in the extended 3′ UTR of the MGMT transcript. Eleven non-coding RNA loci coding for miRNAs, antisense RNA, and lncRNAs were identified in the MGMT-EBF3 region and six of these showed validated potential for curtailing the expression of both MGMT and EBF3 genes. ChIP analysis verified the binding site in MGMT promoter for CTCF which regulates the genomic methylation and chromatin looping. CTCF depletion by a pool of specific siRNA and shRNAs led to a significant attenuation of MGMT expression in human GBM cell lines. Computational analysis of the ChIP sequence data in ENCODE showed the presence of NRF1 in the MGMT promoter and this occurred only in MGMT-proficient cell lines. Further, an enforced NRF2 expression markedly augmented the MGMT mRNA and protein levels in glioma cells. Conclusions: We provide the first evidence for several new regulatory components in the MGMT gene locus which predict complex transcriptional and posttranscriptional controls with potential for new therapeutic avenues.

Genomic and epigenomic EBF1 alterations modulate TERT expression in gastric cancer

Transcriptional reactivation of telomerase catalytic subunit (TERT) is a frequent hallmark of cancer, occurring in 90% of human malignancies. However, specific mechanisms driving TERT reactivation remain obscure for many tumor types and in particular gastric cancer (GC), a leading cause of global cancer mortality. Here, through comprehensive genomic and epigenomic analysis of primary GCs and GC cell lines, we identified the transcription factor early B cell factor 1 (EBF1) as a TERT transcriptional repressor and inactivation of EBF1 function as a major cause of TERT upregulation. Abolishment of EBF1 function occurs through 3 distinct (epi)genomic mechanisms. First, EBF1 is epigenetically silenced via DNA methyltransferase, polycomb-repressive complex 2 (PRC2), and histone deacetylase activity in GCs. Second, recurrent, somatic, and heterozygous EBF1 DNA-binding domain mutations result in the production of dominant-negative EBF1 isoforms. Third, more rarely, genomic deletions and rearrangements proximal to the TERT promoter remobilize or abolish EBF1-binding sites, derepressing TERT and leading to high TERT expression. EBF1 is also functionally required for various malignant phenotypes in vitro and in vivo, highlighting its importance for GC development. These results indicate that multimodal genomic and epigenomic alterations underpin TERT reactivation in GC, converging on transcriptional repressors such as EBF1.

The mRNA landscape profiling reveals potential biomarkers associated with acute kidney injury AKI after kidney transplantation

Background

This study aims to identify potential biomarkers associated with acute kidney injury (AKI) post kidney transplantation.

Material and Methods

Two mRNA expression profiles from Gene Expression Omnibus repertory were downloaded, including 20 delayed graft function (DGF) and 68 immediate graft function (IGF) samples. Differentially expressed genes (DEGs) were identified between DGF and IGF group. The Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis of DEGs were performed. Then, a protein-protein interaction analysis was performed to extract hub genes. The key genes were searched by literature retrieval and cross-validated based on the training dataset. An external dataset was used to validate the expression levels of key genes. Receiver operating characteristic curve analyses were performed to evaluate diagnostic performance of key genes for AKI.

Results

A total of 330 DEGs were identified between DGF and IGF samples, including 179 up-regulated and 151 down-regulated genes. Of these, OLIG3, EBF3 and ETV1 were transcription factor genes. Moreover, LEP, EIF4A3, WDR3, MC4R, PPP2CB, DDX21 and GPT served as hub genes in PPI network. EBF3 was significantly up-regulated in validation GSE139061 dataset, which was consistently with our initial gene differential expression analysis. Finally, we found that LEP had a great diagnostic value for AKI (AUC = 0.740).

Conclusion

EBF3 may be associated with the development of AKI following kidney transplantation. Furthermore, LEP had a good diagnostic value for AKI. These findings provide deeper insights into the diagnosis and management of AKI post renal transplantation.

Biomarkers of response to ibrutinib plus nivolumab in relapsed diffuse large B-cell lymphoma, follicular lymphoma, or Richter's transformation

•

Biomarkers of response to ibrutinib + nivolumab were analyzed in diffuse large B-cell lymphoma (DLBCL), follicular lymphoma and Richter transformation.

•

DLBCL patients with elevated PD-L1 by immunohistochemistry tended to have better response and survival.

•

Whole exome sequencing identified gene mutations in alternate B-cell receptor pathways linked to response in DLBCL.

•

Enriched pathways by gene expression profiling were related to immune activation in responders and proliferation/replication in nonresponders.

•

This preliminary work may help to generate hypotheses on genetically defined subsets of patients most likely to benefit from ibrutinib + nivolumab.

We analyzed potential biomarkers of response to ibrutinib plus nivolumab in biopsies from patients with diffuse large B-cell lymphoma (DLBCL), follicular lymphoma (FL), and Richter's transformation (RT) from the LYM1002 phase I/IIa study, using programmed death ligand 1 (PD-L1) immunohistochemistry, whole exome sequencing (WES), and gene expression profiling (GEP). In DLBCL, PD-L1 elevation was more frequent in responders versus nonresponders (5/8 [62.5%] vs. 3/16 [18.8%]; p = 0.065; complete response 37.5% vs. 0%; p = 0.028). Overall response rates for patients with WES and GEP data, respectively, were: DLBCL (38.5% and 29.6%); FL (46.2% and 43.5%); RT (76.5% and 81.3%). In DLBCL, WES analyses demonstrated that mutations in RNF213 (40.0% vs. 6.2%; p = 0.055), KLHL14 (30.0% vs. 0%; p = 0.046), and LRP1B (30.0% vs. 6.2%; p = 0.264) were more frequent in responders. No responders had mutations in EBF1, ADAMTS20, AKAP9, TP53, MYD88, or TNFRSF14, while the frequency of these mutations in nonresponders ranged from 12.5% to 18.8%. In FL and RT, genes with different mutation frequencies in responders versus nonresponders were: BCL2 (75.0% vs. 28.6%; p = 0.047) and ROS1 (0% vs. 50.0%; p = 0.044), respectively. Per GEP, the most upregulated genes in responders were LEF1 and BTLA (overall), and CRTAM (germinal center B-cell–like DLBCL). Enriched pathways were related to immune activation in responders and resistance-associated proliferation/replication in nonresponders. This preliminary work may help to generate hypotheses regarding genetically defined subsets of DLBCL, FL, and RT patients most likely to benefit from ibrutinib plus nivolumab.

Methylome and transcriptome profiles in three yak tissues revealed that DNA methylation and the transcription factor ZGPAT co-regulate milk production

Background

Domestic yaks play an indispensable role in sustaining the livelihood of Tibetans and other ethnic groups on the Qinghai-Tibetan Plateau (QTP), by providing milk and meat. They have evolved numerous physiological adaptations to high-altitude environment, including strong blood oxygen transportation capabilities and high metabolism. The roles of DNA methylation and gene expression in milk production and high-altitudes adaptation need further exploration.

Results

We performed genome-wide DNA methylome and transcriptome analyses of breast, lung, and biceps brachii muscle tissues from yaks of different ages. We identified 432,350 differentially methylated regions (DMRs) across the age groups within each tissue. The post-mature breast tissue had considerably more differentially methylated regions (155,957) than that from the three younger age groups. Hypomethylated genes with high expression levels might regulate milk production by influencing protein processing in the endoplasmic reticulum. According to weighted gene correlation network analysis, the “hub” gene ZGPAT was highly expressed in the post-mature breast tissue, indicating that it potentially regulates the transcription of 280 genes that influence protein synthesis, processing, and secretion. The tissue network analysis indicated that high expression of HIF1A regulates energy metabolism in the lung.

Conclusions

This study provides a basis for understanding the epigenetic mechanisms underlying milk production in yaks, and the results offer insight to breeding programs aimed at improving milk production.

Supplementary information

Supplementary information accompanies this paper at 10.1186/s12864-020-07151-3.

Sox2-dependent maintenance of mouse oligodendroglioma involves the Sox2-mediated downregulation of Cdkn2b, Ebf1, Zfp423, and Hey2

Cancer stem cells (CSC) are essential for tumorigenesis. The transcription factor Sox2 is overexpressed in brain gliomas, and is essential to maintain CSC. In mouse high-grade glioma pHGG cells in culture, Sox2 deletion causes cell proliferation arrest and inability to reform tumors after transplantation in vivo; in Sox2-deleted cells, 134 genes are derepressed. To identify genes mediating Sox2 deletion effects, we overexpressed into pHGG cells nine among the most derepressed genes, and identified four genes, Ebf1, Hey2, Zfp423, and Cdkn2b, that strongly reduced cell proliferation in vitro and brain tumorigenesis in vivo. CRISPR/Cas9 mutagenesis of each gene, individually or in combination (Ebf1 + Cdkn2b), significantly antagonized the proliferation arrest caused by Sox2 deletion. The same genes also repressed clonogenicity in primary human glioblastoma-derived CSC-like lines. These experiments identify a network of critical tumor suppressive Sox2-targets whose inhibition by Sox2 is involved in glioma CSC maintenance, defining new potential therapeutic targets.

The Importance of Breast Adipose Tissue in Breast Cancer

Adipose tissue is a complex endocrine organ, with a role in obesity and cancer. Adipose tissue is generally linked to excessive body fat, and it is well known that the female breast is rich in adipose tissue. Hence, one can wonder: what is the role of adipose tissue in the breast and why is it required? Adipose tissue as an organ consists of adipocytes, an extracellular matrix (ECM) and immune cells, with a significant role in the dynamics of breast changes throughout the life span of a female breast from puberty, pregnancy, lactation and involution. In this review, we will discuss the importance of breast adipose tissue in breast development and its involvement in breast changes happening during pregnancy, lactation and involution. We will focus on understanding the biology of breast adipose tissue, with an overview on its involvement in the various steps of breast cancer development and progression. The interaction between the breast adipose tissue surrounding cancer cells and vice-versa modifies the tumor microenvironment in favor of cancer. Understanding this mutual interaction and the role of breast adipose tissue in the tumor microenvironment could potentially raise the possibility of overcoming breast adipose tissue mediated resistance to therapies and finding novel candidates to target breast cancer.

Comparative gene expression profiling between optic nerve and spinal cord injury in Xenopus laevis reveals a core set of genes inherent in successful regeneration of vertebrate central nervous system axons

BACKGROUND

The South African claw-toed frog, Xenopus laevis, is uniquely suited for studying differences between regenerative and non-regenerative responses to CNS injury within the same organism, because some CNS neurons (e.g., retinal ganglion cells after optic nerve crush (ONC)) regenerate axons throughout life, whereas others (e.g., hindbrain neurons after spinal cord injury (SCI)) lose this capacity as tadpoles metamorphose into frogs. Tissues from these CNS regions (frog ONC eye, tadpole SCI hindbrain, frog SCI hindbrain) were used in a three-way RNA-seq study of axotomized CNS axons to identify potential core gene expression programs for successful CNS axon regeneration.

RESULTS

Despite tissue-specific changes in expression dominating the injury responses of each tissue, injury-induced changes in gene expression were nonetheless shared between the two axon-regenerative CNS regions that were not shared with the non-regenerative region. These included similar temporal patterns of gene expression and over 300 injury-responsive genes. Many of these genes and their associated cellular functions had previously been associated with injury responses of multiple tissues, both neural and non-neural, from different species, thereby demonstrating deep phylogenetically conserved commonalities between successful CNS axon regeneration and tissue regeneration in general. Further analyses implicated the KEGG adipocytokine signaling pathway, which links leptin with metabolic and gene regulatory pathways, and a novel gene regulatory network with genes regulating chromatin accessibility at its core, as important hubs in the larger network of injury response genes involved in successful CNS axon regeneration.

CONCLUSIONS

This study identifies deep, phylogenetically conserved commonalities between CNS axon regeneration and other examples of successful tissue regeneration and provides new targets for studying the molecular underpinnings of successful CNS axon regeneration, as well as a guide for distinguishing pro-regenerative injury-induced changes in gene expression from detrimental ones in mammals.

Long intergenic noncoding RNA 00844 promotes apoptosis and represses proliferation of prostate cancer cells through upregulating GSTP1 by recruiting EBF1

Accumulating evidence have suggested the function of long noncoding RNAs as crucial players in the pathogenesis of prostate cancer (PC), a urologic tumor in male with poor prognosis. This study was designed to explore the functions of long intergenic noncoding RNA 00844 (LINC00844) in PC progression. The expression of LINC00844 and glutathione S-transferase P1-1 (GSTP1) was detected by reverse transcription quantitative polymerase chain reaction, followed by the identification of the relationship among LINC00844, GSTP1, and early B cell factor 1 (EBF1) by dual luciferase reporter gene assay, RNA immunoprecipitation assay, electrophoretic mobility shift assay, and chromatin immunoprecipitation assay. Using loss- and gain-of-function assays, the effects of LINC00844, GSTP1, and EBF1 on the biological characteristics of PC cells were assessed by cell counting kit-8 assay, 5-ethynyl-2'-deoxyuridine assay, and flow cytometry. Lastly, the results from in vitro experiments were verified in vivo by establishing a xenograft tumor model in nude mice. LINC00844 and GSTP1 both displayed low expression in PC tissues and cells. LINC00844 positively regulated the expression of GSTP1 via recruiting EBF1. Overexpression of LINC00844 reduced proliferation and elevated apoptosis of PC cells through recruiting EBF1, which subsequently upregulated GSTP1. In vivo experiments confirmed that LINC00844 or GSTP1 upregulation attenuated tumor growth. LINC00844 elevated GSTP1 expression by recruiting EBF1 to the promoter region of GSTP1, thereby suppressing PC progression. Hence, LINC00844 is a novel therapeutic target for the development of new treatment protocols for PC.

Phenotypic suppression caused by resonance with light-dark cycles indicates the presence of a 24-hours oscillator in yeast and suggests a new role of intrinsically disordered protein regions as internal mediators

The mutual interaction between environment and life is a main topic of biological sciences. An interesting aspect of this interaction is the existence of biological rhythms spanning all the levels of organisms from bacteria to humans. On the other hand, the existence of a coupling between external oscillatory stimuli and adaptation and evolution rate of biological systems is a still unexplored issue. Here we give the demonstration of a substantial increase of heritable phenotypic changes in yeast, an organism lacking a photoreception system, when growing at 12 h light/dark cycles, with respect to both stable dark (or light) or non-12 + 12 h cycling. The model system was a yeast strain lacking a gene whose product is at the crossroad of many different physiological regulations, so ruling out any simple explanation in terms of increase in reverse gene mutations. The abundance of intrinsically disordered protein regions (IDPRs) in both deleted gene product and in its vast ensemble of interactors supports the hypothesis that resonance with the environmental cycle might be mediated by intrinsic disorder-driven interactions of protein molecules. This result opens to the speculation of the effect of environment/biological resonance phenomena in evolution and of the role of protein intrinsically disordered regions as internal mediators.Communicated by Ramaswamy H. Sarma.

DNA methylation and cis-regulation of gene expression by prostate cancer risk SNPs

Genome-wide association studies have identified more than 100 SNPs that increase the risk of prostate cancer (PrCa). We identify and compare expression quantitative trait loci (eQTLs) and CpG methylation quantitative trait loci (meQTLs) among 147 established PrCa risk SNPs in primary prostate tumors (n = 355 from a Seattle-based study and n = 495 from The Cancer Genome Atlas, TCGA) and tumor-adjacent, histologically benign samples (n = 471 from a Mayo Clinic study). The role of DNA methylation in eQTL regulation of gene expression was investigated by data triangulation using several causal inference approaches, including a proposed adaptation of the Causal Inference Test (CIT) for causal direction. Comparing eQTLs between tumors and benign samples, we show that 98 of the 147 risk SNPs were identified as eQTLs in the tumor-adjacent benign samples, and almost all 34 eQTL identified in tumor sets were also eQTLs in the benign samples. Three lines of results support the causal role of DNA methylation. First, nearly 100 of the 147 risk SNPs were identified as meQTLs in one tumor set, and almost all eQTLs in tumors were meQTLs. Second, the loss of eQTLs in tumors relative to benign samples was associated with altered DNA methylation. Third, among risk SNPs identified as both eQTLs and meQTLs, mediation analyses suggest that over two-thirds have evidence of a causal role for DNA methylation, mostly mediating genetic influence on gene expression. In summary, we provide a comprehensive catalog of eQTLs, meQTLs and putative cancer genes for known PrCa risk SNPs. We observe that a substantial portion of germline eQTL regulatory mechanisms are maintained in the tumor development, despite somatic alterations in tumor genome. Finally, our mediation analyses illuminate the likely intermediary role of CpG methylation in eQTL regulation of gene expression.

Candidate Cancer Driver Mutations in Distal Regulatory Elements and Long-Range Chromatin Interaction Networks

A comprehensive catalog of cancer driver mutations is essential for understanding tumorigenesis and developing therapies. Exome-sequencing studies have mapped many protein-coding drivers, yet few non-coding drivers are known because genome-wide discovery is challenging. We developed a driver discovery method, ActiveDriverWGS, and analyzed 120,788 cis-regulatory modules (CRMs) across 1,844 whole tumor genomes from the ICGC-TCGA PCAWG project. We found 30 CRMs with enriched SNVs and indels (FDR < 0.05). These frequently mutated regulatory elements (FMREs) were ubiquitously active in human tissues, showed long-range chromatin interactions and mRNA abundance associations with target genes, and were enriched in motif-rewiring mutations and structural variants. Genomic deletion of one FMRE in human cells caused proliferative deficiencies and transcriptional deregulation of cancer genes CCNB1IP1, CDH1, and CDKN2B, validating observations in FMRE-mutated tumors. Pathway analysis revealed further sub-significant FMREs at cancer genes and processes, indicating an unexplored landscape of infrequent driver mutations in the non-coding genome.

DAXX in cancer: phenomena, processes, mechanisms and regulation

DAXX displays complex biological functions. Remarkably, DAXX overexpression is a common feature in diverse cancers, which correlates with tumorigenesis, disease progression and treatment resistance. Structurally, DAXX is modular with an N-terminal helical bundle, a docking site for many DAXX interactors (e.g. p53 and ATRX). DAXX's central region folds with the H3.3/H4 dimer, providing a H3.3-specific chaperoning function. DAXX has two functionally critical SUMO-interacting motifs. These modules are connected by disordered regions. DAXX's structural features provide a framework for deciphering how DAXX mechanistically imparts its functions and how its activity is regulated. DAXX modulates transcription through binding to transcription factors, epigenetic modifiers, and chromatin remodelers. DAXX's localization in the PML nuclear bodies also plays roles in transcriptional regulation. DAXX-regulated genes are likely important effectors of its biological functions. Deposition of H3.3 and its interactions with epigenetic modifiers are likely key events for DAXX to regulate transcription, DNA repair, and viral infection. Interactions between DAXX and its partners directly impact apoptosis and cell signaling. DAXX's activity is regulated by posttranslational modifications and ubiquitin-dependent degradation. Notably, the tumor suppressor SPOP promotes DAXX degradation in phase-separated droplets. We summarize here our current understanding of DAXX's complex functions with a focus on how it promotes oncogenesis.

Lack of zinc finger protein 521 upregulates dopamine β-hydroxylase expression in the mouse brain, leading to abnormal behavior

AIM

Previously, we reported that mice deficient in most of the Zfp521 coding region (Zfp521^Δ/Δ mice) displayed abnormal behaviors, including hyperlocomotion and lower anxiety. In this study, we aimed to elucidate the involvement and mechanisms of monoamine variation.

MAIN METHODS

First, we compared the levels of dopamine (DA), noradrenaline (NA), and serotonin in the brains of Zfp521^Δ/Δ and Zfp521^+/+ mice using enzyme-linked immunosorbent assay. Next, we elucidated the mechanisms using quantitative PCR and Western Blotting. Additionally, we administered inhibitory drug to the mice and performed behavioral tests.

KEY FINDINGS

Our results showed that the DA level decreased and the NA level increased in Zfp521^Δ/Δ mice. We found that ZFP521 suppresses the expression of dopamine β-hydroxylase (DBH), which converts DA into NA. We also demonstrated that paired homeodomain transcription factor 2 and early growth response protein-1, which are the transcription factors for Dbh, were involved in the upregulation of Dbh by ZFP521. The administration of nepicastat, a specific inhibitor of DBH, attenuated the abnormal behaviors of Zfp521^Δ/Δ mice.

SIGNIFICANCE

These results suggest that the lack of ZFP521 upregulates the expression of DBH, which leads to a decrease in the DA level and an increase in the NA level in the brain, resulting in abnormal behaviors.

miRNA profile obtained by next‑generation sequencing in metastatic breast cancer patients is able to predict the response to systemic treatments

Metastatic breast cancer (MBC) is a challenge for oncologists, and public efforts should focus on identifying additional molecular markers and therapeutic management to improve clinical outcomes. Among all diagnosed cases of breast cancer (BC; approximately 10%) involve metastatic disease; notably, approximately 40% of patients with early-stage BC develop metastasis within 5 years. The management of MBC consists of systemic therapy. Despite different treatment options, the 5-year survival rate is <20%, which may be due to a lack of response with de novo or acquired resistance. MicroRNAs (miRNAs or miRs) are promising biomarkers as they are readily detectable and have a broad spectrum and potential clinical applications. The aim of this study was to identify a miRNA profile for distinguishing patients with MBC who respond to systemic treatment. Patients with MBC were treated according to the National Comprehensive Cancer Network guidelines. We performed miRNA-Seq on 9 primary tumors using the Thermo Fisher Scientific Ion S5 system. To obtain global miRNA profiles, we carried out differentially expressed gene elimination strategy (DEGES) analysis between the responsive and non-responsive patients. The results identified a profile of 12 miRNAs associated with the response to systemic treatment. The data were validated in an independent cohort (TCGA database). Based on the results, the upregulation of miR-342-3p and miR-187-3p was associated with the response to systemic treatment, and with an increased progression-free survival (PFS) and overall survival (OS); by contrast, the downregulation of miR-301a-3p was associated with a higher PFS and OS. On the whole, the findings of this study indicate that these miRNAs may serve as biomarkers for the response to systemic treatment or the prognosis of patients with MBC. However, these data should be validated experimentally in other robust cohorts and using different specimens before implementing these miRNAs as biomarkers in clinical practice to benefit this group of patients.

EBF1-Mediated Upregulation of Ribosome Assembly Factor PNO1 Contributes to Cancer Progression by Negatively Regulating the p53 Signaling Pathway

The RNA-binding protein PNO1 is critical for ribosome biogenesis, but its potential role in cancer remains unknown. In this study, online data mining, cDNA, and tissue microarrays indicated that PNO1 expression was higher in colorectal cancer tissue than in noncancerous tissue, and its overexpression was associated with worse patient survival. Gain-of-function and loss-of-function studies demonstrated that PNO1 knockdown suppressed growth of colorectal cancer cells in vitro and in vivo, while PNO1 overexpression promoted colorectal cancer cell proliferation in vitro. In colorectal cancer cells expressing wild-type p53, PNO1 knockdown enhanced expression of p53 and its downstream gene p21, and reduced cell viability; these effects were prevented by p53 knockout and attenuated by the p53 inhibitor PFT-α. Moreover, PNO1 knockdown in HCT116 cells decreased levels of 18S rRNA, of 40S and 60S ribosomal subunits, and of the 80S ribosome. It also reduced global protein synthesis, increasing nuclear stress and inhibiting MDM2-mediated ubiquitination and p53 degradation. Overexpressing EBF1 suppressed PNO1 promoter activity and decreased PNO1 mRNA and protein, inhibiting cell proliferation and inducing cell apoptosis through the p53/p21 pathway. In colorectal cancer tissues, the expression of EBF1 correlated inversely with PNO1. Data mining of online breast and lung cancer databases showed increased PNO1 expression and association with poor patient survival; PNO1 knockdown reduced cell viability of cultured breast and lung cancer cells. Taken together, these findings indicate that PNO1 is overexpressed in colorectal cancer and correlates with poor patient survival, and that PNO1 exerts oncogenic effects, at least, in part, by altering ribosome biogenesis. SIGNIFICANCE: This study identifies the ribosome assembly factor PNO1 as a potential oncogene involved in tumor growth and progression of colorectal cancer.

MiR-23b-3p induces the proliferation and metastasis of esophageal squamous cell carcinomas cells through the inhibition of EBF3

MicroRNAs (miRNAs), some small non-coding RNAs that regulate gene expression at the posttranscriptional level, are always aberrantly expressed in carcinomas. In this study, we found that miR-23b-3p was remarkably up-regulated in human esophageal squamous cell carcinoma cells and tissues. Moreover, miR-23b-3p could induce the proliferation, invasion, and metastasis in vitro. EBF3 was identified as the direct downstream target gene of miR-23b-3p and ectogenic EBF3 could strongly inhibit the proliferation, invasion, and metastasis in vitro. Furthermore, it was found that miR-23b-3p could regulate epithelial-to-mesenchymal transition progress by blocking EBF3. Therefore, it was concluded that miR-23b-3p targeted EBF3 to accelerate the proliferation, invasion, and metastasis in ESCC.

Transcription Factor and lncRNA Regulatory Networks Identify Key Elements in Lung Adenocarcinoma

Lung cancer is the second most commonly diagnosed carcinoma and is the leading cause of cancer death. Although significant progress has been made towards its understanding and treatment, unraveling the complexities of lung cancer is still hampered by a lack of comprehensive knowledge on the mechanisms underlying the disease. High-throughput and multidimensional genomic data have shed new light on cancer biology. In this study, we developed a network-based approach integrating somatic mutations, the transcriptome, DNA methylation, and protein-DNA interactions to reveal the key regulators in lung adenocarcinoma (LUAD). By combining Bayesian network analysis with tissue-specific transcription factor (TF) and targeted gene interactions, we inferred 15 disease-related core regulatory networks in co-expression gene modules associated with LUAD. Through target gene set enrichment analysis, we identified a set of key TFs, including known cancer genes that potentially regulate the disease networks. These TFs were significantly enriched in multiple cancer-related pathways. Specifically, our results suggest that hepatitis viruses may contribute to lung carcinogenesis, highlighting the need for further investigations into the roles that viruses play in treating lung cancer. Additionally, 13 putative regulatory long non-coding RNAs (lncRNAs), including three that are known to be associated with lung cancer, and nine novel lncRNAs were revealed by our study. These lncRNAs and their target genes exhibited high interaction potentials and demonstrated significant expression correlations between normal lung and LUAD tissues. We further extended our study to include 16 solid-tissue tumor types and determined that the majority of these lncRNAs have putative regulatory roles in multiple cancers, with a few showing lung-cancer specific regulations. Our study provides a comprehensive investigation of transcription factor and lncRNA regulation in the context of LUAD regulatory networks and yields new insights into the regulatory mechanisms underlying LUAD. The novel key regulatory elements discovered by our research offer new targets for rational drug design and accompanying therapeutic strategies.

ZNF423: A New Player in Estrogen Receptor-Positive Breast Cancer

Preventive therapy can target hormone-responsive breast cancer (BC) by treatment with selective estrogen receptor modulators (SERMs) and reduce the incidence of BC. Genome-wide association studies have identified single nucleotide polymorphisms (SNPs) with relevant predictive values, SNPs in the ZNF423 gene were associated with decreased risk of BC during SERM therapy, and SNPs in the Cathepsin O gene with an increased risk. ZNF423, which was not previously associated with BC is a multifunctional transcription factor known to have a role in development, neurogenesis, and adipogenesis and is implicated in other types of cancer. ZNF423 is transcriptionally controlled by the homolog ZNF521, early B cell factor transcription factor, epigenetic silencing of the promoter by CpG island hyper-methylation, and also by ZNF423 itself in an auto-regulatory loop. In BC cells, ZNF423 expression is found to be induced by estrogen, dependent on the binding of the estrogen receptor and calmodulin-like 3 to SNPs in ZNP423 intronic sites in proximity to consensus estrogen response elements. ZNF423 has also been shown to play a mechanistic role by trans-activating the tumor suppressor BRCA1 and thus modulating the DNA damage response. Even though recent extensive trial studies did not classify these SNPs with the highest predictive values, for inclusion in polygenic SNP analysis, the mechanism unveiled in these studies has introduced ZNF423 as a factor important in the control of the estrogen response. Here, we aim at providing an overview of ZNF423 expression and functional role in human malignancies, with a specific focus on its implication in hormone-responsive BC.

Rare Variant, Gene-Based Association Study of Hereditary Melanoma Using Whole-Exome Sequencing

Background

Extraordinary progress has been made in our understanding of common variants in many diseases, including melanoma. Because the contribution of rare coding variants is not as well characterized, we performed an exome-wide, gene-based association study of familial cutaneous melanoma (CM) and ocular melanoma (OM).

Methods

Using 11 990 jointly processed individual DNA samples, whole-exome sequencing was performed, followed by large-scale joint variant calling using GATK (Genome Analysis ToolKit). PLINK/SEQ was used for statistical analysis of genetic variation. Four models were used to estimate the association among different types of variants. In vitro functional validation was performed using three human melanoma cell lines in 2D and 3D proliferation assays. In vivo tumor growth was assessed using xenografts of human melanoma A375 melanoma cells in nude mice (eight mice per group). All statistical tests were two-sided.

Results

Strong signals were detected for CDKN2A (Pmin = 6.16 × 10-8) in the CM cohort (n = 273) and BAP1 (Pmin = 3.83 × 10-6) in the OM (n = 99) cohort. Eleven genes that exhibited borderline association (P < 10-4) were independently validated using The Cancer Genome Atlas melanoma cohort (379 CM, 47 OM) and a matched set of 3563 European controls with CDKN2A (P = .009), BAP1 (P = .03), and EBF3 (P = 4.75 × 10-4), a candidate risk locus, all showing evidence of replication. EBF3 was then evaluated using germline data from a set of 132 familial melanoma cases and 4769 controls of UK origin (joint P = 1.37 × 10-5). Somatically, loss of EBF3 expression correlated with progression, poorer outcome, and high MITF tumors. Functionally, induction of EBF3 in melanoma cells reduced cell growth in vitro, retarded tumor formation in vivo, and reduced MITF levels.

Conclusions

The results of this large rare variant germline association study further define the mutational landscape of hereditary melanoma and implicate EBF3 as a possible CM predisposition gene.

Prolonged oxidative stress down-regulates Early B cell factor 1 with inhibition of its tumor suppressive function against cholangiocarcinoma genesis

Early B cell factor 1 (EBF1) is a transcription factor involved in the differentiation of several stem cell lineages and it is a negative regulator of estrogen receptors. EBF1 is down-regulated in many tumors, and is believed to play suppressive roles in cancer promotion and progression. However, the functional roles of EBF1 in carcinogenesis are unclear. Liver fluke-infection-associated cholangiocarcinoma (CCA) is an oxidative stress-driven cancer of bile duct epithelium. In this study, we investigated EBF1 expression in tissues from CCA patients, CCA cell lines (KKU-213, KKU-214 and KKU-156), cholangiocyte (MMNK1) and its oxidative stress-resistant (ox-MMNK1-L) cell lines. The formation of 8-oxo-7,8-dihydro-2'-deoxyguanosine (8-oxodG) was used as an oxidative stress marker. Our results revealed that EBF1 expression was suppressed in cancer cells compared with the individual normal bile duct cells at tumor adjacent areas of CCA tissues. CCA patients with low EBF1 expression and high formation of 8-oxodG were shown to correlate with poor survival. Moreover, EBF1 was suppressed in the oxidative stress-resistant cell line and all of CCA cell lines compared to the cholangiocyte cell line. This suggests that prolonged oxidative stress suppressed EBF1 expression and the reduced EBF1 level may facilitate CCA genesis. To elucidate the significance of EBF1 suppression in CCA genesis, EBF1 expression of the MMNK1 cell line was down-regulated by siRNA technique, and its effects on stem cell properties (CD133 and Oct3/4 expressions), tumorigenic properties (cell proliferation, wound healing and cell migration), estrogen responsive gene (TFF1), estrogen-stimulated wound healing, and cell migration were examined. The results showed that CD133, Oct3/4 and TFF1 expression levels, wound healing, and cell migration of EBF1 knockdown-MMNK1 cells were significantly increased. Also, cell migration of EBF1-knockdown cells was significantly enhanced after 17β-estradiol treatment. Our findings suggest that EBF1 down-regulation via oxidative stress induces stem cell properties, tumorigenic properties and estrogen responses of cholangiocytes leading to CCA genesis with aggressive clinical outcomes.

SOX5/6/21 Prevent Oncogene-Driven Transformation of Brain Stem Cells

Molecular mechanisms preventing self-renewing brain stem cells from oncogenic transformation are poorly defined. We show that the expression levels of SOX5, SOX6, and SOX21 (SOX5/6/21) transcription factors increase in stem cells of the subventricular zone (SVZ) upon oncogenic stress, whereas their expression in human glioma decreases during malignant progression. Elevated levels of SOX5/6/21 promoted SVZ cells to exit the cell cycle, whereas genetic ablation of SOX5/6/21 dramatically increased the capacity of these cells to form glioma-like tumors in an oncogene-driven mouse brain tumor model. Loss-of-function experiments revealed that SOX5/6/21 prevent detrimental hyperproliferation of oncogene expressing SVZ cells by facilitating an antiproliferative expression profile. Consistently, restoring high levels of SOX5/6/21 in human primary glioblastoma cells enabled expression of CDK inhibitors and decreased p53 protein turnover, which blocked their tumorigenic capacity through cellular senescence and apoptosis. Altogether, these results provide evidence that SOX5/6/21 play a central role in driving a tumor suppressor response in brain stem cells upon oncogenic insult. Cancer Res; 77(18); 4985-97. ©2017 AACR.

Genome-wide methylation sequencing of paired primary and metastatic cell lines identifies common DNA methylation changes and a role for EBF3 as a candidate epigenetic driver of melanoma metastasis

Epigenetic alterations are increasingly implicated in metastasis, whereas very few genetic mutations have been identified as authentic drivers of cancer metastasis. Yet, to date, few studies have identified metastasis-related epigenetic drivers, in part because a framework for identifying driver epigenetic changes in metastasis has not been established. Using reduced representation bisulfite sequencing (RRBS), we mapped genome-wide DNA methylation patterns in three cutaneous primary and metastatic melanoma cell line pairs to identify metastasis-related epigenetic drivers. Globally, metastatic melanoma cell lines were hypomethylated compared to the matched primary melanoma cell lines. Using whole genome RRBS we identified 75 shared (10 hyper- and 65 hypomethylated) differentially methylated fragments (DMFs), which were associated with 68 genes showing significant methylation differences. One gene, Early B Cell Factor 3 (EBF3), exhibited promoter hypermethylation in metastatic cell lines, and was validated with bisulfite sequencing and in two publicly available independent melanoma cohorts (n = 40 and 458 melanomas, respectively). We found that hypermethylation of the EBF3 promoter was associated with increased EBF3 mRNA levels in metastatic melanomas and subsequent inhibition of DNA methylation reduced EBF3 expression. RNAi-mediated knockdown of EBF3 mRNA levels decreased proliferation, migration and invasion in primary and metastatic melanoma cell lines. Overall, we have identified numerous epigenetic changes characterising metastatic melanoma cell lines, including EBF3-induced aggressive phenotypic behaviour with elevated EBF3 expression in metastatic melanoma, suggesting that EBF3 promoter hypermethylation may be a candidate epigenetic driver of metastasis.