RUNX2 Promotes Malignant Progression in Glioma

Runx2 Suppresses Astrocyte Activation and Astroglial Scar Formation After Spinal Cord Injury in Mice

After spinal cord injury, astrocytes undergo a reactive process and form an astroglial scar, which impedes the regeneration of axons. The role of Runx2 in promoting the transformation of astrocytes in the central nervous system is well-established. However, it remains unclear whether Runx2 also plays a role in the development of astroglial scar, and the precise underlying mechanism has yet to be identified. Recently, our study using cell culture and animal models has demonstrated that Runx2 actually suppresses astrocyte activation and the formation of astroglial scar following injury. The initial results demonstrated an increase in the expression of Runx2 in astrocytes following in vivo injury. Subsequently, the overexpression of Runx2 resulted in the inhibition of astrocyte activation, reduction in the total area of astroglial scar, and restoration of neural function after 14 days of injury. However, these effects were reversed by CADD522. These findings indicate that Runx2 could potentially serve as a therapeutic intervention for spinal cord injury (SCI). Furthermore, our findings suggest that the Nuclear-matrix-targeting signal (NMTS) of Runx2 is associated with its effect. In summary, the study's results propose that targeting Runx2 may be a promising treatment approach for reactive astrocytes and astroglial scar in the recovery of SCI.

Multi-omic integration of DNA methylation and gene expression data reveals molecular vulnerabilities in glioblastoma

Glioblastoma (GBM) is one of the most aggressive types of cancer and exhibits profound genetic and epigenetic heterogeneity, making the development of an effective treatment a major challenge. The recent incorporation of molecular features into the diagnosis of patients with GBM has led to an improved categorization into various tumour subtypes with different prognoses and disease management. In this work, we have exploited the benefits of genome-wide multi-omic approaches to identify potential molecular vulnerabilities existing in patients with GBM. Integration of gene expression and DNA methylation data from both bulk GBM and patient-derived GBM stem cell lines has revealed the presence of major sources of GBM variability, pinpointing subtype-specific tumour vulnerabilities amenable to pharmacological interventions. In this sense, inhibition of the AP-1, SMAD3 and RUNX1/RUNX2 pathways, in combination or not with the chemotherapeutic agent temozolomide, led to the subtype-specific impairment of tumour growth, particularly in the context of the aggressive, mesenchymal-like subtype. These results emphasize the involvement of these molecular pathways in the development of GBM and have potential implications for the development of personalized therapeutic approaches.

Cell-cycle dependent nuclear gene delivery enhances the effects of E-cadherin against tumor invasion and metastasis

Gene delivery is the process by which foreign DNA is transferred to host cells, released from intracellular vesicles, and transported to the nuclei for transcription. This process is frequently inefficient and difficult to control spatiotemporally. We developed a gene delivery strategy that uses ultrasound to directly deliver plasmid DNA into nuclei via gas vesicles (GVs)-based intracellular cavitation. pDNA-binding GVs can be taken up by cells and cause intracellular cavitation when exposed to acoustic irradiation and delivering their pDNA payloads into nuclei. Importantly, GVs can remain stable in the cytoplasm in the absence of acoustic irradiation, allowing for temporally controlled nuclear gene delivery. We were able to achieve spatiotemporal control of E-cadherin nuclear gene delivery in this manner, demonstrating its efficacy in tumor invasion and metastasis inhibition. Interestingly, we discovered that nuclear gene delivery of E-cadherin during the G2/M phase of the cell cycle in C6 tumor cells inhibited tumor invasion and metastasis more effectively than during the G1 and S phases. The gene delivery of E-cadherin at the G2/M phase resulted in significantly lower expression of Fam50a, which reduced Fam50a/Runx2 interaction and led to reduced transactivation of MMP13, an important factor for epithelial-mesenchymal transition, as observed in a molecular mechanism assay. Thus, using remote acoustic control of intracellular cavitation of pDNA-GVs, we developed a high spatiotemporally controllable gene delivery strategy and achieved stronger tumor invasion and metastasis inhibition effects by delivering the E-cadherin gene at the G2/M phase.

Long Non-coding RNA KTN1-AS1 Targets miR-505 to Promote Glioblastoma Progression

Glioblastoma (GBM) is a highly malignant cancer, the prognosis of which is pretty poor. Long non-coding RNAs (lncRNAs) are a class of non-coding RNAs, which play important roles in carcinogenesis process of many cancers including GBM. In this study, we want to clarify the expression, biological function, and molecular mechanism of lncRNA KTN1 antisense RNA 1 (KTN1-AS1) in GBM tumor progression. We found that KTN1-AS1 expression was upregulated in GBM tissues and cell lines. KTN1-AS1 played oncogenic roles to facilitate proliferation, migration, and invasion of GBM cells. Then, we revealed that miR-505 was a target of KTN1-AS1, and its expression was decreased in GBM. KTN1-AS1 contributed to GBM progression by mediating miR-505. Finally, we demonstrated that KTN1-AS1 upregulated some target oncogenes of miR-505 including ZEB2, HMGB1, and RUNX2 in GBM cells. All in all, we concluded that the highly expressed KTN1-AS1 in GBM played oncogenic roles to facilitate GBM progression by targeting miR-505.

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.

E2F transcription factor 1 elevates cyclin D1 expression by suppressing transcription of microRNA-107 to augment progression of glioma

BACKGROUND

Dysregulation of microRNAs has been frequently implicated in the progression of human diseases, including glioma. This study aims to explore the interaction between E2F transcription factor 1 (E2F1) and miR-107 in the progression of glioma.

METHODS

Expression of miR-107 in glioma tissues and cells was examined. Putative binding sites between E2F1 and the promoter region of miR-107, and between miR-107 and cyclin D1 (CCND1) mRNA were predicted via bioinformatic systems and validated via chromatin immunoprecipitation and luciferase reporter gene assays. Altered expression of miR-107, E2F1, and CCND1 was introduced in A172 and T98G cells to examine their roles in cell growth and the activity of the Wnt/β-catenin signaling. In vivo experiments were performed by injecting cells in nude mice.

RESULTS

miR-107 was poorly expressed, whereas E2F1 and CCND1 were highly expressed in glioma tissues and cells. E2F1 bound to the promoter region of miR-107 to induce transcriptional repression, and miR-107 directly bound to CCND1 mRNA to reduce its expression. Overexpression of miR-107 reduced proliferation, migration and invasion, and augmented apoptosis of glioma cells, and it reduced activity of the Wnt/β-catenin pathway. The anti-tumorigenic roles of miR-107 were blocked by E2F1 or CCND1 overexpression. Similar results were reproduced in vivo where miR-107 overexpression or E2F1 inhibition blocked tumor growth in nude mice.

CONCLUSION

This study suggested that E2F1 reduces miR-107 transcription to induce CCND1 upregulation, which leads to progression of glioma via Wnt/β-catenin signaling activation.

Analysis of Circulating miRNA Profile in Plasma Samples of Glioblastoma Patients

(1) Background: Glioblastoma multiforme (GBM) is among the most aggressive cancers with a poor prognosis. Treatment options are limited, clinicians lack efficient prognostic and predictive markers. Circulating miRNAs—besides being important regulators of cancer development—may have potential as diagnostic biomarkers of GBM. (2) Methods: In this study, profiling of 798 human miRNAs was performed on blood plasma samples from 6 healthy individuals and 6 patients with GBM, using a NanoString nCounter Analysis System. To validate our results, five miRNAs (hsa-miR-433-3p, hsa-miR-362-3p, hsa-miR-195-5p, hsa-miR-133a-3p, and hsa-miR-29a-3p) were randomly chosen for RT-qPCR detection. (3) Results: In all, 53 miRNAs were significantly differentially expressed in plasma samples of GBM patients when data were filtered for FC 1 and FDR 0.1. Target genes of the top 39 differentially expressed miRNAs were identified, and we carried out functional annotation and pathway enrichment analysis of target genes via GO and KEGG-based tools. General and cortex-specific protein–protein interaction networks were constructed from the target genes of top miRNAs to assess their functional connections. (4) Conclusions: We demonstrated that plasma microRNA profiles are promising diagnostic and prognostic molecular biomarkers that may find an actual application in the clinical practice of GBM, although more studies are needed to validate our results.

KRAS gene polymorphisms are associated with the risk of glioma: a two-center case-control study

BACKGROUND

Glioma, also known as neuroglioma, is the most common primary tumors of the central nervous system. Many previous studies have reported associations between RAS gene polymorphisms and multiple tumors. However, the role of RAS gene polymorphisms on glioma risk has not been investigated.

METHODS

We conducted a two-center case-control study to investigate whether the RAS gene polymorphisms predispose individuals to gliomas in 248 healthy controls and 191 glioma patients. RAS gene polymorphisms (rs12587 G>T, rs7973450 A>G, rs7312175 G>A in KRAS, rs2273267 A>T in NRAS) were genotyped by the TaqMan assay. The relationship between RAS gene functional single nucleotide polymorphisms (SNPs) and the risk of glioma was evaluated based on odds ratios (ORs) and 95% confidence intervals (CIs).

RESULTS

Individuals with KRAS rs7312175 GA genotype were more likely to develop glioma than those with GG genotype (adjusted OR =1.66, 95% CI: 1.05-2.64, P=0.030). However, the other three SNPs could not affect glioma risk. In stratified analysis of age, gender, subtypes, and clinical stages, rs7312175 GA carriers were more likely to develop glioma in the following subgroups: children less than 60 months, tumor derived from the astrocytic tumors, and clinical stages I.

CONCLUSIONS

The study showed that polymorphism rs7312175 GA in the KRAS gene was associated with increased glioma susceptibility. Further investigation is warranted to confirm these findings and to better elucidate the involved biological pathways.

RUNX2 interacts with BRG1 to target CD44 for promoting invasion and migration of colorectal cancer cells

Background

Cancer stem cells (CSCs) play an important role in tumor invasion and metastasis. CD44 is the most commonly used marker of CSCs, with the potential to act as a determinant against the invasion and migration of CSCs and as the key factor in epithelial-mesenchymal transition (EMT)-like changes that occur in colorectal cancer (CRC). Runt-related transcription factor-2 (RUNX2) is a mesenchymal stem marker for cancer that is involved in stem cell biology and tumorigenesis. However, whether RUNX2 is involved in CSC and in inducing EMT-like changes in CRC remains uncertain, warranting further investigation.

Methods

We evaluated the role of RUNX2 in the invasion and migration of CRC cells as a promoter of CD44-induced stem cell- and EMT-like modifications. For this purpose, western blotting was employed to analyze the expression of differential proteins in CRC cells. We conducted sphere formation, wound healing, and transwell assays to investigate the biological functions of RUNX2 in CRC cells. Cellular immunofluorescence and coimmunoprecipitation (co-IP) assays were performed to study the relationship between RUNX2 and BRG1. Real-time quantitative PCR (RT-qPCR) and immunohistochemistry (IHC) were performed to analyze the expressions of RUNX2, BRG1, and CD44 in the CRC tissues.

Results

We found that RUNX2 could markedly induce the CRC cell sphere-forming ability and EMT. Interestingly, the RUNX2-mediated EMT in CRC cell may be associated with the activation of CD44. Furthermore, RUNX2 was found to interact with BRG1 to promote the recruitment of RUNX2 to the CD44 promoter.

Conclusions

Our cumulative findings suggest that RUNX2 and BRG1 can form a compact complex to regulate the transcription and expression of CD44, which has possible involvement in the invasion and migration of CRC cells.

Clinical significance of transcription factor RUNX2 in lung adenocarcinoma and its latent transcriptional regulating mechanism

RUNX family transcription factor 2 (RUNX2) overexpression has been found in various human malignancies. However, the expression levels of RUNX2 mRNA and protein in lung adenocarcinoma (LUAD) were not investigated. This study aims to thoroughly analysis the expression level and potential mechanisms of RUNX2 mRNA in LUAD. We applied in-house immunohistochemistry, high-throughput RNA-sequencing, and gene microarrays to comprehensively investigate the expression level of RUNX2 in LUAD. A pool standard mean difference (SMD) and summary receiver operating characteristic curves (SROC) were calculated to assess the integrated expression value of RUNX2 in LUAD. The hazard ratios (HRs) were integrated to evaluate the overall prognostic effect of RUNX2 on the LUAD patients. The differentially expressed genes (DEGs) of LUAD, the potential target genes of RUNX2, and its co-expressed genes were overlapped to obtain a set of specific genes for GO and KEGG enrichment analyses. RUNX2 overexpression in LUAD was validated using a large number of cases (2 418 LUAD and 1 574 non-tumor lung samples). The pooled SMD was 0.85 (95 % CI: 0.64-1.05) and the area under the curve (AUC) of the SROC was 0.86 (95 %CI: 0.83-0.89). The integrated HR was 1.20 [1.04-1.38], indicating that increased expression of RUNX2 was an independent risk factor for the poor survival of the LUAD patients. RUNX2 and its transcriptionally regulates potential target genes may promote cell proliferation and drug resistance of LUAD by modulating the cell cycle and MAPK signaling pathways. RUNX2 can provide new research directions for targeted drug therapy and drug resistance for LUAD treatment.

Glucose-induced microRNA-218 suppresses the proliferation and promotes the apoptosis of human retinal pigment epithelium cells by targeting RUNX2

OBJECTIVE

MicroRNA-218 (miR-218) critical for preventing the progression of numerous diseases, including diseases of the retinal pigment epithelium (RPE). However, the mechanism by which miR-218 regulates the PRE in humans remains largely unknown. Our study investigated the effects of glucose-induced miR-218 expression on human RPE cells (ARPE-19), as well as its targeted regulatory effect.

METHODS

The levels of miR-218 and runt-related transcription factor 2 (RUNX2) expression were investigated by RT-qPCR or Western blot assays. Cell viability and apoptosis were assessed by CCK-8 assays, flow cytometry, and Hoechst staining. A luciferase reporter assay was performed to determine whether Runx2 is a target gene of miR-218.

RESULTS

Our results showed that glucose up-regulated miR-218 expression, suppressed proliferation, and induced the apoptosis of ARPE-19 cells. We verified that miR-218 could inhibit the proliferation and facilitate the apoptosis of ARPE-19 cells, while inhibition of miR-218 expression produced the opposite effects. In terms of mechanism, we demonstrated that RUNX2 was a direct target of miR-218. Functional experiments showed that Runx2 served as a miR-218 target to help inhibit the proliferation and induction of apoptosis in ARPE-19 cells.

CONCLUSION

Our findings suggest the miR-218/Runx2 axis as a potential target for treating diabetic retinopathy (DR).

Knowledge-primed neural networks enable biologically interpretable deep learning on single-cell sequencing data

BACKGROUND

Deep learning has emerged as a versatile approach for predicting complex biological phenomena. However, its utility for biological discovery has so far been limited, given that generic deep neural networks provide little insight into the biological mechanisms that underlie a successful prediction. Here we demonstrate deep learning on biological networks, where every node has a molecular equivalent, such as a protein or gene, and every edge has a mechanistic interpretation, such as a regulatory interaction along a signaling pathway.

RESULTS

With knowledge-primed neural networks (KPNNs), we exploit the ability of deep learning algorithms to assign meaningful weights in multi-layered networks, resulting in a widely applicable approach for interpretable deep learning. We present a learning method that enhances the interpretability of trained KPNNs by stabilizing node weights in the presence of redundancy, enhancing the quantitative interpretability of node weights, and controlling for uneven connectivity in biological networks. We validate KPNNs on simulated data with known ground truth and demonstrate their practical use and utility in five biological applications with single-cell RNA-seq data for cancer and immune cells.

CONCLUSIONS

We introduce KPNNs as a method that combines the predictive power of deep learning with the interpretability of biological networks. While demonstrated here on single-cell sequencing data, this method is broadly relevant to other research areas where prior domain knowledge can be represented as networks.

ZNF521 which is downregulated by miR-802 suppresses malignant progression of Hepatocellular Carcinoma through regulating Runx2 expression

Zinc finger protein 521 (ZNF521) plays an important role in the tumor development and process. However, its regulatory role in hepatocellular carcinoma (HCC) remains unclear. In this study, we demonstrated for the first time that ZNF521 mRNA and protein was down-regulated in HCC tissues and cell lines. Down-regulated ZNF521 expression was significantly associated with malignant prognostic features, including advanced TNM stage and large tumor size. For 5-year survival, ZNF521 served as a potential prognostic marker of HCC patients. Moreover, ZNF521 inhibited cell proliferation, colony formation and cell viability through Runx2 transcriptional inhibition and AKT phosphorylation pathway. Moreover, we demonstrated that ZNF521 expression was regulated by miR-802. In HCC tissues. MiR-802 has an inverse correlation with ZNF521 expression. In conclusion, we demonstrate for the first time that ZNF521 is down-regulated in HCC tissues and inhibits HCC growth through Runx2 transcriptional inhibition and AKT inactivation, which was regulated by miR-802, suggesting the potential therapeutic value for HCC.

The Thyroid Hormone Receptor-RUNX2 Axis: A Novel Tumor Suppressive Pathway in Breast Cancer

Metastatic breast cancer is refractory to conventional therapies and is an end-stage disease. RUNX2 is a transcription factor that becomes oncogenic when aberrantly expressed in multiple tumor types, including breast cancer, supporting tumor progression and metastases. Our previous work demonstrated that the thyroid hormone receptor beta (TRβ) inhibits RUNX2 expression and tumorigenic characteristics in thyroid cells. As TRβ is a tumor suppressor, we investigated the compelling question whether TRβ also regulates RUNX2 in breast cancer. The Cancer Genome Atlas indicates that TRβ expression is decreased in the most aggressive basal-like subtype of breast cancer. We established that modulated levels of TRβ results in corresponding changes in the high levels of RUNX2 expression in metastatic, basal-like breast cells. The MDA-MB-231 triple-negative breast cancer cell line exhibits low expression of TRβ and high levels of RUNX2. Increased expression of TRβ decreased RUNX2 levels. The thyroid hormone-mediated suppression of RUNX2 is TRβ specific as TRα overexpression failed to alter RUNX2 expression. Consistent with these findings, knockdown of TRβ in non-tumor MCF10A mammary epithelial-like cells results in an increase in RUNX2 and RUNX2 target genes. Mechanistically, TRβ directly interacts with the proximal promoter of RUNX2 through a thyroid hormone response element to reduce promoter activity. The TRβ suppression of the oncogene RUNX2 is a signaling pathway shared by thyroid and breast cancers. Our findings provide a novel mechanism for TRβ-mediated tumor suppression in breast cancers. This pathway may be common to many solid tumors and impact treatment for metastatic cancers.

Identification of Secreted Phosphoprotein 1 (SPP1) as a Prognostic Factor in Lower-Grade Gliomas

OBJECTIVE

Secreted phosphoprotein 1 (SPP1) is an important extracellular glycoprotein that is associated with immune regulation, tumorigenesis, and cell signaling. However, the prognostic value of SPP1 in patients with glioma has not yet been clarified, especially in lower-grade gliomas. The objective of this study is to evaluate the prognostic merit of SPP1 in lower-grade gliomas.

METHODS

The messenger RNA (mRNA) expression of SPP1 in about 1000 cancer cell lines was explored by using the data from the Cancer Cell Line Encyclopedia database. The Oncomine database was mined to evaluate the mRNA expression of SPP1 in lower-grade glioma, glioblastoma, and normal brain tissues. The correlation between SPP1 mRNA expression and overall survival of patients with glioma from The Cancer Genome Atlas database was analyzed.

RESULTS

SPP1 mRNA expression of glioma was ranked as the eighth highest of all cancer cell lines in the Cancer Cell Line Encyclopedia database. The data from the Oncomine database suggested that SPP1 expression was significantly high in glioblastoma compared with normal brain tissues but was not significantly high in lower-grade glioma compared with normal brain tissue. Analysis of the RNA-Seq data from The Cancer Genome Atlas database showed that the increased SPP1 mRNA expression in lower-grade glioma was significantly associated with poor survival outcomes in patients with lower-grade glioma. Multivariate Cox regression analysis showed that SPP1 might be considered as an independent prognostic factor in lower-grade gliomas.

CONCLUSIONS

The present study showed that SPP1 overexpression is related to worse overall survival in patients with lower-grade glioma. Moreover, SPP1 could be considered as an independent factor in lower-grade gliomas.

HDACs control RUNX2 expression in cancer cells through redundant and cell context-dependent mechanisms

Background

RUNX2 is a Runt-related transcription factor required during embryogenesis for skeletal development and morphogenesis of other organs including thyroid and breast gland. Consistent evidence indicates that RUNX2 expression is aberrantly reactivated in cancer and supports tumor progression. The mechanisms leading to RUNX2 expression in cancer has only recently began to emerge. Previously, we showed that suppressing the activity of the epigenetic regulators HDACs significantly represses RUNX2 expression highlighting a role for these enzymes in RUNX2 reactivation in cancer. However, the molecular mechanisms by which HDACs control RUNX2 are still largely unexplored. Here, to fill this gap, we investigated the role of different HDACs in RUNX2 expression regulation in breast and thyroid cancer, tumors that majorly rely on RUNX2 for their development and progression.

Methods

Proliferation assays and evaluation of RUNX2 mRNA levels by qRT-PCR were used to evaluate the effect of several HDACi and specific siRNAs on a panel of cancer cell lines. Moreover, ChIP and co-IP assays were performed to elucidate the molecular mechanism underneath the RUNX2 transcriptional regulation. Finally, RNA-sequencing unveiled a new subset of genes whose transcription is regulated by the complex RUNX2-HDAC6.

Results

In this study, we showed that Class I HDACs and in particular HDAC1 are required for RUNX2 efficient transcription in cancer. Furthermore, we found an additional and cell-specific function of HDAC6 in driving RUNX2 expression in thyroid cancer cells. In this model, HDAC6 likely stabilizes the assembly of the transcriptional complex, which includes HDAC1, on the RUNX2 P2 promoter potentiating its transcription. Since a functional interplay between RUNX2 and HDAC6 has been suggested, we used RNA-Seq profiling to consolidate this evidence in thyroid cancer and to extend the knowledge on this cooperation in a setting in which HDAC6 also controls RUNX2 expression.

Conclusions

Overall, our data provide new insights into the molecular mechanisms controlling RUNX2 in cancer and consolidate the rationale for the use of HDACi as potential pharmacological strategy to counteract the pro-oncogenic program controlled by RUNX2 in cancer cells.

Electronic supplementary material

The online version of this article (10.1186/s13046-019-1350-5) contains supplementary material, which is available to authorized users.