Loss of EGFR-ASAP1 signaling in metastatic and unresectable hepatoblastoma

Identification and experimental validation of druggable epigenetic targets in hepatoblastoma

BACKGROUND & AIMS

Hepatoblastoma (HB) is the most frequent childhood liver cancer. Patients with aggressive tumors have limited therapeutic options; therefore, a better understanding of HB pathogenesis is needed to improve treatment. HBs have a very low mutational burden; however, epigenetic alterations are increasingly recognized. We aimed to identify epigenetic regulators consistently dysregulated in HB and to evaluate the therapeutic efficacy of their targeting in clinically relevant models.

METHODS

We performed a comprehensive transcriptomic analysis of 180 epigenetic genes. Data from fetal, pediatric, adult, peritumoral (n = 72) and tumoral (n = 91) tissues were integrated. Selected epigenetic drugs were tested in HB cells. The most relevant epigenetic target identified was validated in primary HB cells, HB organoids, a patient-derived xenograft model, and a genetic mouse model. Transcriptomic, proteomic and metabolomic mechanistic analyses were performed.

RESULTS

Altered expression of genes regulating DNA methylation and histone modifications was consistently observed in association with molecular and clinical features of poor prognosis. The histone methyltransferase G9a was markedly upregulated in tumors with epigenetic and transcriptomic traits of increased malignancy. Pharmacological targeting of G9a significantly inhibited growth of HB cells, organoids and patient-derived xenografts. Development of HB induced by oncogenic forms of β-catenin and YAP1 was ablated in mice with hepatocyte-specific deletion of G9a. We observed that HBs undergo significant transcriptional rewiring in genes involved in amino acid metabolism and ribosomal biogenesis. G9a inhibition counteracted these pro-tumorigenic adaptations. Mechanistically, G9a targeting potently repressed the expression of c-MYC and ATF4, master regulators of HB metabolic reprogramming.

CONCLUSIONS

HBs display a profound dysregulation of the epigenetic machinery. Pharmacological targeting of key epigenetic effectors exposes metabolic vulnerabilities that can be leveraged to improve the treatment of these patients.

IMPACT AND IMPLICATIONS

In spite of recent advances in the management of hepatoblastoma (HB), treatment resistance and drug toxicity are still major concerns. This systematic study reveals the remarkable dysregulation in the expression of epigenetic genes in HB tissues. Through pharmacological and genetic experimental approaches, we demonstrate that the histone-lysine-methyltransferase G9a is an excellent drug target in HB, which can also be harnessed to enhance the efficacy of chemotherapy. Furthermore, our study highlights the profound pro-tumorigenic metabolic rewiring of HB cells orchestrated by G9a in coordination with the c-MYC oncogene. From a broader perspective, our findings suggest that anti-G9a therapies may also be effective in other c-MYC-dependent tumors.

ASAP1 Expression in Invasive Breast Cancer and Its Prognostic Role

Breast cancer is a major global health burden with high morbidity and mortality rates. Previous studies have reported that increased expression of ASAP1 is associated with poor prognosis in various types of cancer. This study was conducted on 452 breast cancer patients who underwent surgery at Hanyang University Hospital, Seoul, South Korea. Data on clinicopathological characteristics including molecular pathologic markers were collected. Immunohistochemical staining of ASAP1 expression level were used to classify patients into high and low groups. In total, 452 cases low ASAP1 expression group was associated with significantly worse recurrence-free survival (p = 0.029). In ER-positive cases (n = 280), the low ASAP1 expression group was associated with significantly worse overall survival (p = 0.039) and recurrence-free survival (p = 0.029). In multivariate cox analysis, low ASAP1 expression was an independent significant predictor of poor recurrence-free survival in the overall patient group (hazard ratio = 2.566, p = 0.002) and ER-positive cases (hazard ratio = 4.046, p = 0.002). In the analysis of the TCGA dataset, the low-expression group of ASAP1 protein demonstrated a significantly poorer progression-free survival (p = 0.005). This study reports that low ASAP1 expression was associated with worse recurrence-free survival in invasive breast cancer.

GepLiver: an integrative liver expression atlas spanning developmental stages and liver disease phases

Chronic liver diseases usually developed through stepwise pathological transitions under the persistent risk factors. The molecular changes during liver transitions are pivotal to improve liver diagnostics and therapeutics yet still remain elusive. Cumulative large-scale liver transcriptomic studies have been revealing molecular landscape of various liver conditions at bulk and single-cell resolution, however, neither single experiment nor databases enabled thorough investigations of transcriptomic dynamics along the progression of liver diseases. Here we establish GepLiver, a longitudinal and multidimensional liver expression atlas integrating expression profiles of 2469 human bulk tissues, 492 mouse samples, 409,775 single cells from 347 human samples and 27 liver cell lines spanning 16 liver phenotypes with uniformed processing and annotating methods. Using GepLiver, we have demonstrated dynamic changes of gene expression, cell abundance and crosstalk harboring meaningful biological associations. GepLiver can be applied to explore the evolving expression patterns and transcriptomic features for genes and cell types respectively among liver phenotypes, assisting the investigation of liver transcriptomic dynamics and informing biomarkers and targets for liver diseases.

Identify potential driver genes for PAX-FOXO1 fusion-negative rhabdomyosarcoma through frequent gene co-expression network mining

Background

Rhabdomyosarcoma (RMS) is a soft tissue sarcoma usually originated from skeletal muscle. Currently, RMS classification based on PAX-FOXO1 fusion is widely adopted. However, compared to relatively clear understanding of the tumorigenesis in the fusion-positive RMS, little is known for that in fusion-negative RMS (FN-RMS).

Methods

We explored the molecular mechanisms and the driver genes of FN-RMS through frequent gene co-expression network mining (fGCN), differential copy number (CN) and differential expression analyses on multiple RMS transcriptomic datasets.

Results

We obtained 50 fGCN modules, among which five are differentially expressed between different fusion status. A closer look showed 23% of Module 2 genes are concentrated on several cytobands of chromosome 8. Upstream regulators such as MYC, YAP1, TWIST1 were identified for the fGCN modules. Using in a separate dataset we confirmed that, comparing to FP-RMS, 59 Module 2 genes show consistent CN amplification and mRNA overexpression, among which 28 are on the identified chr8 cytobands. Such CN amplification and nearby MYC (also resides on one of the above cytobands) and other upstream regulators (YAP1, TWIST1) may work together to drive FN-RMS tumorigenesis and progression. Up to 43.1% downstream targets of Yap1 and 45.8% of the targets of Myc are differentially expressed in FN-RMS vs. normal comparisons, which also confirmed the driving force of these regulators.

Discussion

We discovered that copy number amplification of specific cytobands on chr8 and the upstream regulators MYC, YAP1 and TWIST1 work together to affect the downstream gene co-expression and promote FN-RMS tumorigenesis and progression. Our findings provide new insights for FN-RMS tumorigenesis and offer promising targets for precision therapy. Experimental investigation about the functions of identified potential drivers in FN-RMS are in progress.

Inhibition of Heat Shock Factor 1 Signaling Decreases Hepatoblastoma Growth via Induction of Apoptosis

Although rare compared with adult liver cancers, hepatoblastoma (HB) is the most common pediatric liver malignancy, and its incidence is increasing. Currently, the treatment includes surgical resection with or without chemotherapy, and in severe cases, liver transplantation in children. The effort to develop more targeted, HB-specific therapies has been stymied by the lack of fundamental knowledge about HB biology. Heat shock factor 1 (HSF1), a transcription factor, is a canonical inducer of heat shock proteins, which act as chaperone proteins to prevent or undo protein misfolding. Recent work has shown a role for HSF1 in cancer beyond the canonical heat shock response. The current study found increased HSF1 signaling in HB versus normal liver. It showed that less differentiated, more embryonic tumors had higher levels of HSF1 than more differentiated, more fetal-appearing tumors. Most strikingly, HSF1 expression levels correlated with mortality. This study used a mouse model of HB to test the effect of inhibiting HSF1 early in tumor development on cancer growth. HSF1 inhibition resulted in fewer and smaller tumors, suggesting HSF1 is needed for aggressive tumor growth. Moreover, HSF1 inhibition also increased apoptosis in tumor foci. These data suggest that HSF1 may be a viable pharmacologic target for HB treatment.

Loss of ASAP1 in the MMTV-PyMT model of luminal breast cancer activates AKT, accelerates tumorigenesis, and promotes metastasis

ASAP1 is a multi-domain adaptor protein that regulates cytoskeletal dynamics, receptor recycling and intracellular vesicle trafficking. Its expression is associated with poor prognosis in a variety of cancers, and can promote cell migration, invasion and metastasis. Although amplification and expression of ASAP1 has been associated with poor survival in breast cancer, we found that in the autochthonous MMTV-PyMT model of luminal breast cancer, ablation of ASAP1 resulted in an earlier onset of tumor initiation and increased metastasis. This was due to tumor cell-intrinsic effects of ASAP1 deletion, as ASAP1 deficiency in tumor, but not in stromal cells was sufficient to replicate the enhanced tumorigenicity and metastasis observed in the ASAP1-null MMTV-PyMT mice. Loss of ASAP1 in MMTV-PyMT mice had no effect on proliferation, apoptosis, angiogenesis or immune cell infiltration, but enhanced mammary gland hyperplasia and tumor cell invasion, indicating that ASAP1 can accelerate tumor initiation and promote dissemination. Mechanistically, these effects were associated with a potent activation of AKT. Importantly, lower ASAP1 levels correlated with poor prognosis and enhanced AKT activation in human ER+/luminal breast tumors, validating our findings in the MMTV-PyMT mouse model for this subtype of breast cancer. Taken together, our findings reveal that ASAP1 can have distinct functions in different tumor types and demonstrate a tumor suppressive activity for ASAP1 in luminal breast cancer.

New insights in gastrointestinal "pediatric" neoplasms in adult patients: pancreatoblastoma, hepatoblastoma and embryonal sarcoma of the liver. A practical approach by GIPPI-GIPAD Groups

Pediatric solid neoplasms are rare and very different from those observed in adults. The majority of them are referred to as embryonal because they arise as a result of alterations in the processes of organogenesis or normal growth and are characterized by proliferation of primitive cells, reproducing the corresponding tissue at various stages of embryonic development. This review will focus on embryonal gastrointestinal pediatric neoplasms in adult patients, including pancreatoblastoma, hepatoblastoma, and embryonal sarcoma of the liver. Although they are classically considered pediatric neoplasms, they may (rarely) occur in adult patients. Hepatoblastoma represents the most frequent liver neoplasm in the pediatric population, followed by hepatocellular carcinoma and embryonal sarcoma of the liver; while pancreatoblastoma is the most common malignant pancreatic tumor in childhood. Both in children and adults, the mainstay of treatment is complete surgical resection, either up front or following neoadjuvant chemotherapy. Unresectable and/or metastatic neoplasms may be amenable to complete delayed surgery after neoadjuvant chemotherapy. However, these neoplasms display a more aggressive behavior and overall poorer prognosis in adults than in children, probably because they are diagnosed in later stages of diseases.

A Summary of the Inaugural WHO Classification of Pediatric Tumors: Transitioning from the Optical into the Molecular Era

Pediatric tumors are uncommon, yet are the leading cause of cancer-related death in childhood. Tumor types, molecular characteristics, and pathogenesis are unique, often originating from a single genetic driver event. The specific diagnostic challenges of childhood tumors led to the development of the first World Health Organization (WHO) Classification of Pediatric Tumors. The classification is rooted in a multilayered approach, incorporating morphology, IHC, and molecular characteristics. The volume is organized according to organ sites and provides a single, state-of-the-art compendium of pediatric tumor types. A special emphasis was placed on "blastomas," which variably recapitulate the morphologic maturation of organs from which they originate. SIGNIFICANCE: In this review, we briefly summarize the main features and updates of each chapter of the inaugural WHO Classification of Pediatric Tumors, including its rapid transition from a mostly microscopic into a molecularly driven classification systematically taking recent discoveries in pediatric tumor genomics into account.

Human cell transformation by combined lineage conversion and oncogene expression

Cancer is the most complex genetic disease known, with mutations implicated in more than 250 genes. However, it is still elusive which specific mutations found in human patients lead to tumorigenesis. Here we show that a combination of oncogenes that is characteristic of liver cancer (CTNNB1, TERT, MYC) induces senescence in human fibroblasts and primary hepatocytes. However, reprogramming fibroblasts to a liver progenitor fate, induced hepatocytes (iHeps), makes them sensitive to transformation by the same oncogenes. The transformed iHeps are highly proliferative, tumorigenic in nude mice, and bear gene expression signatures of liver cancer. These results show that tumorigenesis is triggered by a combination of three elements: the set of driver mutations, the cellular lineage, and the state of differentiation of the cells along the lineage. Our results provide direct support for the role of cell identity as a key determinant in transformation and establish a paradigm for studying the dynamic role of oncogenic drivers in human tumorigenesis.

Upregulation of human endogenous retrovirus-K (HML-2) mRNAs in hepatoblastoma: Identification of potential new immunotherapeutic targets and biomarkers

PURPOSE

Hepatoblastoma is the most common liver malignancy in children. In order to advance therapy against hepatoblastoma, novel immunologic targets and biomarkers are needed. Our purpose in this investigation is to examine hepatoblastoma transcriptomes for the expression of a class of genomic elements known as Human Endogenous Retrovirus (HERVs). HERVs are abundant in the human genome and are biologically active elements that have been associated with multiple malignancies and proposed as immunologic targets in a subset of tumors. A sub-family of HERVs, HERV-K(HML-2) (HERV-K), have been shown to be tightly regulated in fetal development, making investigation of these elements in pediatric tumors paramount.

METHODS

We first created a HERVK-FASTA file utilizing 91 previously described HML-2 proviruses. We then concatenated the file onto the GRCh38.95 cDNA library from Ensembl. We used this reference database to evaluate existing RNA-seq data from 10 hepatoblastoma tumors and 3 normal liver controls (GEO accession ID: GSE8977575). Quantification and differential proviral expression analysis between hepatoblastoma and normal liver controls was performed using the pseudo-alignment program Salmon and DESeq2, respectively.

RESULTS

HERV-K mRNA was expressed in hepatoblastoma from multiple proviral loci. All expressed HERV-K proviral loci were upregulated in hepatoblastoma compared to normal liver controls. Five HERV-K proviruses (1q21.3, 3q27.2, 7q22.2, 12q24.33 and 17p13.1) were significantly differentially expressed (p-adjusted value <0.05, |log2 fold change| > 1.5) across conditions. The provirus at 17p13.1 had an approximately 300-fold increased expression in hepatoblastoma as compared to normal liver. This was in part due to the near absence of HERV-K mRNA at the 17p13.1 locus in fully differentiated liver samples.

CONCLUSIONS

Our investigation demonstrates that HERV-K is expressed from multiple loci in hepatoblastoma and that the expression is increased for several proviruses compared to normal liver controls. Our results suggest that HERV-K mRNA expression may be useful as a biomarker in hepatoblastoma, given the large differential expression profiles in hepatoblastoma, with very low mRNA levels in liver control samples.

Hsa_circ_0008934 promotes the proliferation and migration of osteosarcoma cells by targeting miR-145-5p to enhance E2F3 expression

OBJECTIVE

To investigate the role of hsa_circ_0008934 in osteosarcoma and the molecular mechanism involved in the regulation of the occurrence and development of osteosarcoma METHODS: Differentially expressed circRNAs in the osteosarcoma cell lines SaOS2 and MG63 and in the normal human osteoblast cell line hFOB1.19 were identified via next-generation RNA sequencing. The expression and circular morphology of hsa_circ_0008934 were analyzed via quantitative real-time polymerase chain reaction (qRT-PCR) and RT-PCR analysis, respectively. Proliferation, apoptosis, cell cycle progression, migration, and invasion of SaOS2 and MG63 cells with hsa_circ_0008934 silencing or overexpression were assessed using the MTS method, colony formation assay, flow cytometry, and the transwell system, respectively. The subcellular distribution of hsa_circ_0008934 was revealed via fluorescence in situ hybridization. The binding of hsa_circ_0008934 with microRNAs was confirmed using the dual-luciferase reporter assay. The oncogenic roles of hsa_circ_0008934 in osteosarcoma were determined using an in vivo tumorigenesis assay with nude mice. qRT-PCR, western blotting, TUNEL assay, and immunohistochemistry (IHC) were used to detect the tumorigenicity of hsa_circ_0008934 in osteosarcoma cells.

RESULTS

Many circRNAs were differentially expressed in SaOS2 and MG63 cells than in hFOB1.19 cells. Hsa_circ_0008934 expression was significantly elevated in SaOS2 and MG63 cells. Hsa_circ_0008934 silencing significantly reduced proliferation, enhanced apoptosis, blocked cell cycle progression, and impaired migration and invasion capacities of SaOS2 cells. Opposite cellular alterations were achieved by overexpressing hsa_circ_0008934 in MG63 cells. Hsa_circ_0008934 was mainly distributed in the cytosol and positively regulated E2F3 expression in osteosarcoma cells. In addition, it directly bound with miR-145-5p to repress E2F3 expression and enhanced the tumorigenesis of MG63 cells in nude mice. qRT-PCR revealed that the intracellular injection of hsa_circ_0008934 lentivirus resulted in hsa_circ_0008934 overexpression and miR-145-5p downregulation. Western blotting confirmed that E2F3 was upregulated. Moreover, the TUNEL assay showed that hsa_circ_0008934 overexpression inhibited the apoptosis of tumor cells. IHC detection revealed that the hsa_circ_0008934 overexpression could promote the expression of Ki67 and PCNA.

CONCLUSION

Elevated hsa_circ_0008934 expression promotes the proliferation and migration of osteosarcoma cells by sponging miR-145-5p to enhance E2F3 expression.

Human endogenous retrovirus-K mRNA expression and genomic alignment data in hepatoblastoma

Human Endogenous Retroviruses are a class of genomic elements that are the result of ancient retroviral infection of the human germline. Many are biologically active elements that have been implicated in multiple diseases including cancer. The most recent class to invade the human genome is the HERV-K(HML-2) (HERV-K) family. Approximately 90 HERV-K proviruses and many smaller elements have been identified to date in the human genome. Additional proviruses are continually being discovered with the rapid advancement of deep-sequencing and long-read sequencing technologies. HERV-K proviruses are poorly annotated in human transcriptome databases making their analysis in RNA-seq data difficult. To enable analysis, we compiled the sequences of 91 HERV-K proviruses identified in NCBI GenBank (ID JN675007-JN675097) and created a proviral alignment tool for visualizing RNA-seq reads aligned across individual proviruses. This allowed us to analyse publicly available RNA-seq data from 10 hepatoblastoma samples and 3 normal liver controls (GEO Accession ID: GSE89775). This data report includes the raw FASTA sequence files of the HERV-K proviruses from NCBI, a differential gene expression list between hepatoblastoma samples, and genomic alignment figures from 5 HERV-K proviruses identified as differentially expressed in the companion research article "Upregulation of Human Endogenous Retrovirus-K (HML-2) mRNAs in hepatoblastoma: Identification of potential new immunotherapeutic targets and biomarkers [1]. The data provided here are available for other research groups interested in evaluating individual HERV-K proviral expression using RNA-seq data. Furthermore, the data analysis is highly flexible and will accommodate the addition of other HERV-K proviruses.

Volasertib preclinical activity in high-risk hepatoblastoma

Relapsed and metastatic hepatoblastoma represents an unmet clinical need with limited chemotherapy treatment options. In a chemical screen, we identified volasertib as an agent with in vitro activity, inhibiting hepatoblastoma cell growth while sparing normal hepatocytes. Volasertib targets PLK1 and prevents the progression of mitosis, resulting in eventual cell death. PLK1 is overexpressed in hepatoblastoma biopsies relative to normal liver tissue. As a potential therapeutic strategy, we tested the combination of volasertib and the relapse-related hepatoblastoma chemotherapeutic irinotecan. We found both in vitro and in vivo efficacy of this combination, which may merit further preclinical investigation and exploration for a clinical trial concept.

Neurokinin-1 Receptor Antagonists against Hepatoblastoma

Hepatoblastoma (HB) is the most common malignant liver tumor that occurs during childhood. The prognosis of children with HB is favorable when a complete surgical resection of the tumor is possible, but for high-risk patients, the prognosis is much worse. New anti-HB strategies must be urgently developed. The undecapeptide substance P (SP) after binding to the neurokinin-1 receptor (NK-1R), regulates cancer cell proliferation, exerts an antiapoptotic effect, induces cell migration for invasion/metastasis, and triggers endothelial cell proliferation for neoangiogenesis. HB samples and cell lines overexpress NK-1R (the truncated form) and SP elicits HB cell proliferation. One of these strategies could be the use of non-peptide NK-1R antagonists. These antagonists exert, in a concentration-dependent manner, an antiproliferative action against HB cells (inhibit cell proliferation and induce the death of HB cells by apoptosis). NK-1R antagonists exerted a dual effect in HB: Decreased both tumor volume and angiogenic activity. Thus, the SP/NK-1R system is an important target in the HB treatment and NK-1R antagonists could act as specific drugs against HB cells. In this review, we update and discuss the use of NK-1R antagonists in the treatment of HB.

The presence of PIM3 increases hepatoblastoma tumorigenesis and tumor initiating cell phenotype and is associated with decreased patient survival

PURPOSE

Hepatoblastoma is the most common primary liver cancer of childhood and has few prognostic indicators. We have previously shown that Proviral Integration site for Moloney murine leukemia virus (PIM3) kinase decreased hepatoblastoma tumorigenicity. We sought to determine the effect of PIM3 overexpression on hepatoblastoma cells and whether expression of PIM3 correlated with patient/tumor characteristics or survival.

METHODS

The hepatoblastoma cell line, HuH6, and patient-derived xenograft, COA67, were utilized. Viability, proliferation, migration, sphere formation, and tumor growth in mice were assessed in PIM3-overexpressing cells. Immunohistochemistry was performed for PIM3 on patient samples. Correlation between stain score and clinical/pathologic characteristics was assessed.

RESULTS

PIM3 overexpression rescued the anti-proliferative effect observed with PIM3 knockdown. Sphere formation was increased in PIM3 overexpressing cells. Cells with PIM3 overexpression yielded larger tumors than those with empty vector. Seventy-four percent of samples expressed PIM3. There was no statistical difference in patient characteristics between subjects with strong versus weak PIM3 staining, but patients with strong PIM3 staining had decreased survival.

CONCLUSIONS

PIM3 expression plays a role in hepatoblastoma tumorigenesis. PIM3 was present in the majority of hepatoblastomas and higher PIM3 expression correlated with decreased survival. PIM3 warrants investigation as a therapeutic target and prognostic marker for hepatoblastoma.

Identification of differentially expressed genes, lncRNAs and miRNAs which are associated with tumor malignant phenotypes in hepatoblastoma patients

Hepatoblastoma (HB) is one of the most common hepatic malignancies in the pediatric population. HB are composed of a variety of tumors, which derived from different origins and had varying clinical outcomes. However, the unclear underlying mechanisms of HB limited exploring novel biomarkers and effective therapeutic targets. We searched microarray datasets on Gene Expression Omnibus (GEO) database and selected GSE75271 and GSE75283 datasets for comprehensive analysis. Weighted gene correlation network analysis (WGCNA) was employed to identify genes which were associated with tumor malignant phenotypes, including HB subtypes, Cairo classification and tumor stage. Coexpression analysis of identified genes was also performed and lncRNA-miRNA-mRNA network was finally conducted. Our results showed that a total of 22 lncRNAs, 13 miRNAs and 66 mRNAs were identified to be associated with tumor malignant phenotypes. Mechanistically, these molecules might promote the malignant phenotypes via regulating metabolic pathways. Among of them, 6 miRNAs (hsa-miR-106b, hsa-miR-130b, hsa-miR-19a, hsa-miR-19b, hsa-miR-20a and hsa-miR-301a), 8 lncRNAs (NR_102317, XR_245338, XR_428373, XR_924945, XR_929728, XR_931611, XR_935074 and XR_946696), and 6 mRNAs (EGFR, GAREM, INSIG1, KRT81, SAR1B and SDC1) were selected to conduct a lncRNA-miRNA-mRNA network. Taken together, our findings provide evidence for exploring molecular mechanisms of HB. Those identified malignant phenotype-associated molecules might be potential biomarkers and anti-cancer therapeutic targets in future.