Targeting of SUMOylation leads to cBAF complex stabilization and disruption of the SS18::SSX transcriptome in Synovial Sarcoma.

Synovial Sarcoma (SS) is driven by the SS18::SSX fusion oncoprotein. and is ultimately refractory to therapeutic approaches. SS18::SSX alters ATP-dependent chromatin remodeling BAF (mammalian SWI/SNF) complexes, leading to the degradation of canonical (cBAF) complex and amplified presence of an SS18::SSX-containing non-canonical BAF (ncBAF or GBAF) that drives an SS-specific transcription program and tumorigenesis. We demonstrate that SS18::SSX activates the SUMOylation program and SSs are sensitive to the small molecule SAE1/2 inhibitor, TAK-981. Mechanistically, TAK-981 de-SUMOylates the cBAF subunit SMARCE1, stabilizing and restoring cBAF on chromatin, shifting away from SS18::SSX-ncBAF-driven transcription, associated with DNA damage and cell death and resulting in tumor inhibition across both human and mouse SS tumor models. TAK-981 synergized with cytotoxic chemotherapy through increased DNA damage, leading to tumor regression. Targeting the SUMOylation pathway in SS restores cBAF complexes and blocks the SS18::SSX-ncBAF transcriptome, identifying a therapeutic vulnerability in SS, positioning the in-clinic TAK-981 to treat SS.

A novel multi-omics mendelian randomization method for gene set enrichment and its application to psychiatric disorders.

Genome-wide association studies (GWAS) of psychiatric disorders (PD) yield numerous loci with significant signals, but often do not implicate specific genes. Because GWAS risk loci are enriched in expression/protein/methylation quantitative loci (e/p/mQTL, hereafter xQTL), transcriptome/proteome/methylome-wide association studies (T/P/MWAS, hereafter XWAS) that integrate xQTL and GWAS information, can link GWAS signals to effects on specific genes. To further increase detection power, gene signals are aggregated within relevant gene sets (GS) by performing gene set enrichment (GSE) analyses. Often GSE methods test for enrichment of "signal" genes in curated GS while overlooking their linkage disequilibrium (LD) structure, allowing for the possibility of increased false positive rates. Moreover, no GSE tool uses xQTL information to perform mendelian randomization (MR) analysis. To make causal inference on association between PD and GS, we develop a novel MR GSE (MR-GSE) procedure. First, we generate a "synthetic" GWAS for each MSigDB GS by aggregating summary statistics for x-level (mRNA, protein or DNA methylation (DNAm) levels) from the largest xQTL studies available) of genes in a GS. Second, we use synthetic GS GWAS as exposure in a generalized summary-data-based-MR analysis of complex trait outcomes. We applied MR-GSE to GWAS of nine important PD. When applied to the underpowered opioid use disorder GWAS, none of the four analyses yielded any signals, which suggests a good control of false positive rates. For other PD, MR-GSE greatly increased the detection of GO terms signals (2,594) when compared to the commonly used (non-MR) GSE method (286). Some of the findings might be easier to adapt for treatment, e.g., our analyses suggest modest positive effects for supplementation with certain vitamins and/or omega-3 for schizophrenia, bipolar and major depression disorder patients. Similar to other MR methods, when applying MR-GSE researchers should be mindful of the confounding effects of horizontal pleiotropy on statistical inference.

Genome-wide analysis of hepatic DNA methylation reveals impact of epigenetic aging on xenobiotic metabolism and transport genes in an aged mouse model.

Hepatic xenobiotic metabolism and transport decline with age, while intact xenobiotic metabolism is associated with longevity. However, few studies have examined the genome-wide impact of epigenetic aging on these processes. We used reduced representation bisulfite sequencing (RRBS) to map DNA methylation changes in liver DNA from mice ages 4 and 24 months. We identified several thousand age-associated differentially methylated sites (a-DMS), many of which overlapped genes encoding Phase I and Phase II drug metabolizing enzymes, in addition to ABC and SLC classes of transporters. Notable genes harboring a-DMS were Cyp1a2, Cyp2d9, and Abcc2 that encode orthologs of the human drug metabolizing enzymes CYP1A2 and CYP2D6, and the multidrug resistance protein 2 (MRP2) transporter. Cyp2d9 hypermethylation with age was significantly associated with reduced gene expression, while Abcc2 expression was unchanged with age. Cyp1a2 lost methylation with age while, counterintuitively, its expression also reduced with age. We hypothesized that age-related dysregulation of the hepatic transcriptional machinery caused down-regulation of genes despite age-related hypomethylation. Bioinformatic analysis of hypomethylated a-DMS in our sample found them to be highly enriched for hepatic nuclear factor 4 alpha (HNF4α) binding sites. HNF4α promotes Cyp1a2 expression and is downregulated with age, which could explain the reduction in Cyp1a2 expression. Overall, our study supports the broad impact of epigenetic aging on xenobiotic metabolism and transport. Future work should evaluate the interplay between hepatic nuclear receptor function and epigenetic aging. These results may have implications for studies of longevity and healthy aging.

Machine and deep learning methods for predicting 3D genome organization.

Three-Dimensional (3D) chromatin interactions, such as enhancer-promoter interactions (EPIs), loops, Topologically Associating Domains (TADs), and A/B compartments play critical roles in a wide range of cellular processes by regulating gene expression. Recent development of chromatin conformation capture technologies has enabled genome-wide profiling of various 3D structures, even with single cells. However, current catalogs of 3D structures remain incomplete and unreliable due to differences in technology, tools, and low data resolution. Machine learning methods have emerged as an alternative to obtain missing 3D interactions and/or improve resolution. Such methods frequently use genome annotation data (ChIP-seq, DNAse-seq, etc.), DNA sequencing information (k-mers, Transcription Factor Binding Site (TFBS) motifs), and other genomic properties to learn the associations between genomic features and chromatin interactions. In this review, we discuss computational tools for predicting three types of 3D interactions (EPIs, chromatin interactions, TAD boundaries) and analyze their pros and cons. We also point out obstacles of computational prediction of 3D interactions and suggest future research directions.

Inducible deletion of Ezh2 in CD4+ T cells inhibits kidney T cell infiltration and prevents interstitial nephritis in MRL/ lpr lupus-prone mice.

Systemic lupus erythematosus is a remitting relapsing autoimmune disease characterized by autoantibody production and multi-organ involvement. T cell epigenetic dysregulation plays an important role in the pathogenesis of lupus. We have previously demonstrated upregulation of the key epigenetic regulator EZH2 in CD4+ T cells isolated from lupus patients. To further investigate the role of EZH2 in the pathogenesis of lupus, we generated a tamoxifen-inducible CD4+ T cell Ezh2 conditional knockout mouse on the MRL/ lpr lupus-prone background. We demonstrate that Ezh2 deletion abrogates lupus-like disease and prevents T cell differentiation. Single-cell analysis suggests impaired T cell function and activation of programed cell death pathways in EZH2-deficient mice. Ezh2 deletion in CD4+ T cells restricts TCR clonal repertoire and prevents kidney-infiltrating effector CD4+ T cell expansion and tubulointerstitial nephritis, which has been linked to end-stage renal disease in patients with lupus nephritis.

Coordinate transcriptional regulation of ErbB2/3 by C-terminal binding protein 2 signals sensitivity to ErbB2 inhibition in pancreatic adenocarcinoma.

There is a critical need to identify new therapeutic vulnerabilities in pancreatic ductal adenocarcinoma (PDAC). Transcriptional co-regulators C-terminal binding proteins (CtBP) 1 and 2 are highly overexpressed in human PDAC, and CRISPR-based homozygous deletion of Ctbp2 in a mouse PDAC cell line (CKP) dramatically decreased tumor growth, reduced metastasis, and prolonged survival in orthotopic mouse allografts. Transcriptomic profiling of tumors derived from CKP vs. Ctbp2-deleted CKP cells (CKP/KO) revealed significant downregulation of the EGFR-superfamily receptor Erbb3, the heterodimeric signaling partner for both EGFR and ErbB2. Compared with CKP cells, CKP/KO cells also demonstrated reduced Erbb2 expression and did not activate downstream Akt signaling after stimulation of Erbb3 by its ligand neuregulin-1. ErbB3 expression in human PDAC cell lines was similarly dependent on CtBP2 and depletion of ErbB3 in a human PDAC cell line severely attenuated growth, demonstrating the critical role of ErbB3 signaling in maintaining PDAC cell growth. Sensitivity to the ErbB2-targeted tyrosine kinase inhibitor lapatinib, but not the EGFR-targeted agent erlotinib, varied in proportion to the level of ErbB3 expression in mouse and human PDAC cells, suggesting that an ErBb2 inhibitor can effectively leverage CtBP2-driven transcriptional activation of physiologic ErbB2/3 expression and signaling in PDAC cells for therapeutic benefit.

Dysregulation of Type I Interferon (IFN-I) Signaling: A Potential Contributor to Racial Disparity in Hepatocellular Carcinoma (HCC).

African-American (AA)/Black hepatocellular carcinoma (HCC) patients have increased incidence and decreased survival rates compared to non-Hispanic (White) patients, the underlying molecular mechanism of which is not clear. Analysis of existing RNA-sequencing (RNA-seq) data in The Cancer Genome Atlas (TCGA) and in-house RNA-sequencing of 14 White and 18 AA/Black HCC patients revealed statistically significant activation of type I interferon (IFN-I) signaling pathway in AA/Black patients. A four-gene signature of IFN-stimulated genes (ISGs) showed increased expression in AA/Black HCC tumors versus White. HCC is a disease of chronic inflammation, and IFN-Is function as pro-inflammatory cytokines. We tested efficacy of ginger extract (GE), a dietary compound known for anti-inflammatory properties, on HCC cell lines derived from White (HepG2), AA/Black (Hep3B and O/20) and Asian (HuH-7) patients. GE exhibited a significantly lower IC50 on Hep3B and O/20 cells than on HepG2 and HuH-7 cells. The GE treatment inhibited the activation of downstream mediators of IFN-I signaling pathways and expression of ISGs in all four HCC cells. Our data suggest that ginger can potentially attenuate IFN-I-mediated signaling pathways in HCC, and cells from AA/Black HCC patients may be more sensitive to ginger. AA/Black HCC patients might benefit from a holistic diet containing ginger.

Yin Yang 1 controls cerebellar astrocyte maturation.

Diverse subpopulations of astrocytes tile different brain regions to accommodate local requirements of neurons and associated neuronal circuits. Nevertheless, molecular mechanisms governing astrocyte diversity remain mostly unknown. We explored the role of a zinc finger transcription factor Yin Yang 1 (YY1) that is expressed in astrocytes. We found that specific deletion of YY1 from astrocytes causes severe motor deficits in mice, induces Bergmann gliosis, and results in simultaneous loss of GFAP expression in velate and fibrous cerebellar astrocytes. Single cell RNA-seq analysis showed that YY1 exerts specific effects on gene expression in subpopulations of cerebellar astrocytes. We found that although YY1 is dispensable for the initial stages of astrocyte development, it regulates subtype-specific gene expression during astrocyte maturation. Moreover, YY1 is continuously needed to maintain mature astrocytes in the adult cerebellum. Our findings suggest that YY1 plays critical roles regulating cerebellar astrocyte maturation during development and maintaining a mature phenotype of astrocytes in the adult cerebellum.

MBD2a-NuRD binds to the methylated γ-globin gene promoter and uniquely forms a complex required for silencing of HbF expression.

During human development, there is a switch in the erythroid compartment at birth that results in silencing of expression of fetal hemoglobin (HbF). Reversal of this silencing has been shown to be effective in overcoming the pathophysiologic defect in sickle cell anemia. Among the many transcription factors and epigenetic effectors that are known to mediate HbF silencing, two of the most potent are BCL11A and MBD2-NuRD. In this report, we present direct evidence that MBD2-NuRD occupies the γ-globin gene promoter in adult erythroid cells and positions a nucleosome there that results in a closed chromatin conformation that prevents binding of the transcriptional activator, NF-Y. We show that the specific isoform, MBD2a, is required for the formation and stable occupancy of this repressor complex that includes BCL11A, MBD2a-NuRD, and the arginine methyltransferase, PRMT5. The methyl cytosine binding preference and the arginine-rich (GR) domain of MBD2a are required for high affinity binding to methylated γ-globin gene proximal promoter DNA sequences. Mutation of the methyl cytosine-binding domain (MBD) of MBD2 results in a variable but consistent loss of γ-globin gene silencing, in support of the importance of promoter methylation. The GR domain of MBD2a is also required for recruitment of PRMT5, which in turn results in placement of the repressive chromatin mark H3K8me2s at the promoter. These findings support a unified model that integrates the respective roles of BCL11A, MBD2a-NuRD, PRMT5, and DNA methylation in HbF silencing.

Stratification of Tamoxifen Synergistic Combinations for the Treatment of ER+ Breast Cancer.

Breast cancer alone accounts for the majority of cancer deaths among women, with the most commonly diagnosed subtype being estrogen receptor positive (ER+). Survival has greatly improved for patients with ER+ breast cancer, due in part to the development of antiestrogen compounds, such as tamoxifen. While treatment of the primary disease is often successful, as many as 30% of patients will experience recurrence and metastasis, mainly due to developed endocrine therapy resistance. In this study, we discovered two tamoxifen combination therapies, with simeprevir and VX-680, that reduce the tumor burden in animal models of ER+ breast cancer more than either compound or tamoxifen alone. Additionally, these tamoxifen combinations reduced the expression of HER2, a hallmark of tamoxifen treatment, which can facilitate acquisition of a treatment-resistant phenotype. These combinations could provide clinical benefit by potentiating tamoxifen treatment in ER+ breast cancer.

matchRanges: generating null hypothesis genomic ranges via covariate-matched sampling.

MOTIVATION: Deriving biological insights from genomic data commonly requires comparing attributes of selected genomic loci to a null set of loci. The selection of this null set is non-trivial, as it requires careful consideration of potential covariates, a problem that is exacerbated by the non-uniform distribution of genomic features including genes, enhancers, and transcription factor binding sites. Propensity score-based covariate matching methods allow the selection of null sets from a pool of possible items while controlling for multiple covariates; however, existing packages do not operate on genomic data classes and can be slow for large data sets making them difficult to integrate into genomic workflows. RESULTS: To address this, we developed matchRanges, a propensity score-based covariate matching method for the efficient and convenient generation of matched null ranges from a set of background ranges within the Bioconductor framework. AVAILABILITY AND IMPLEMENTATION: Package: https://bioconductor.org/packages/nullranges, Code: https://github.com/nullranges, Documentation: https://nullranges.github.io/nullranges.

excluderanges: exclusion sets for T2T-CHM13, GRCm39, and other genome assemblies.

SUMMARY: Exclusion regions are sections of reference genomes with abnormal pileups of short sequencing reads. Removing reads overlapping them improves biological signal, and these benefits are most pronounced in differential analysis settings. Several labs created exclusion region sets, available primarily through ENCODE and Github. However, the variety of exclusion sets creates uncertainty which sets to use. Furthermore, gap regions (e.g. centromeres, telomeres, short arms) create additional considerations in generating exclusion sets. We generated exclusion sets for the latest human T2T-CHM13 and mouse GRCm39 genomes and systematically assembled and annotated these and other sets in the excluderanges R/Bioconductor data package, also accessible via the BEDbase.org API. The package provides unified access to 82 GenomicRanges objects covering six organisms, multiple genome assemblies, and types of exclusion regions. For human hg38 genome assembly, we recommend hg38.Kundaje.GRCh38_unified_blacklist as the most well-curated and annotated, and sets generated by the Blacklist tool for other organisms. AVAILABILITY AND IMPLEMENTATION: https://bioconductor.org/packages/excluderanges/. Package website: https://dozmorovlab.github.io/excluderanges/.

bootRanges: flexible generation of null sets of genomic ranges for hypothesis testing.

MOTIVATION: Enrichment analysis is a widely utilized technique in genomic analysis that aims to determine if there is a statistically significant association between two sets of genomic features. To conduct this type of hypothesis testing, an appropriate null model is typically required. However, the null distribution that is commonly used can be overly simplistic and may result in inaccurate conclusions. RESULTS: bootRanges provides fast functions for generation of block bootstrapped genomic ranges representing the null hypothesis in enrichment analysis. As part of a modular workflow, bootRanges offers greater flexibility for computing various test statistics leveraging other Bioconductor packages. We show that shuffling or permutation schemes may result in overly narrow test statistic null distributions and over-estimation of statistical significance, while creating new range sets with a block bootstrap preserves local genomic correlation structure and generates more reliable null distributions. It can also be used in more complex analyses, such as accessing correlations between cis-regulatory elements (CREs) and genes across cell types or providing optimized thresholds, e.g. log fold change (logFC) from differential analysis. AVAILABILITY AND IMPLEMENTATION: bootRanges is freely available in the R/Bioconductor package nullranges hosted at https://bioconductor.org/packages/nullranges.

Rewiring of the 3D genome during acquisition of carboplatin resistance in a triple-negative breast cancer patient-derived xenograft.

Changes in the three-dimensional (3D) structure of the genome are an emerging hallmark of cancer. Cancer-associated copy number variants and single nucleotide polymorphisms promote rewiring of chromatin loops, disruption of topologically associating domains (TADs), active/inactive chromatin state switching, leading to oncogene expression and silencing of tumor suppressors. However, little is known about 3D changes during cancer progression to a chemotherapy-resistant state. We integrated chromatin conformation capture (Hi-C), RNA-seq, and whole-genome sequencing obtained from triple-negative breast cancer patient-derived xenograft primary tumors (UCD52) and carboplatin-resistant samples and found increased short-range (< 2 Mb) interactions, chromatin looping, formation of TAD, chromatin state switching into a more active state, and amplification of ATP-binding cassette transporters. Transcriptome changes suggested the role of long-noncoding RNAs in carboplatin resistance. Rewiring of the 3D genome was associated with TP53, TP63, BATF, FOS-JUN family of transcription factors and led to activation of aggressiveness-, metastasis- and other cancer-related pathways. Integrative analysis highlighted increased ribosome biogenesis and oxidative phosphorylation, suggesting the role of mitochondrial energy metabolism. Our results suggest that 3D genome remodeling may be a key mechanism underlying carboplatin resistance.

Discovering Synergistic Compounds with BYL-719 in PI3K Overactivated Basal-like PDXs.

Basal-like triple-negative breast cancer (TNBC) tumor cells are difficult to eliminate due to resistance mechanisms that promote survival. While this breast cancer subtype has low PIK3CA mutation rates when compared to estrogen receptor-positive (ER+) breast cancers, most basal-like TNBCs have an overactive PI3K pathway due to gene amplification or high gene expression. BYL-719 is a PIK3CA inhibitor that has been found to have low drug-drug interactions, which increases the likelihood that it could be useful for combinatorial therapy. Alpelisib (BYL-719) with fulvestrant was recently approved for treating ER+ breast cancer patients whose cancer had developed resistance to ER-targeting therapy. In these studies, a set of basal-like patient-derived xenograft (PDX) models was transcriptionally defined with bulk and single-cell RNA-sequencing and clinically actionable mutation profiles defined with Oncomine mutational profiling. This information was overlaid onto therapeutic drug screening results. BYL-719-based, synergistic two-drug combinations were identified with 20 different compounds, including everolimus, afatinib, and dronedarone, which were also found to be effective at minimizing tumor growth. These data support the use of these drug combinations towards cancers with activating PIK3CA mutations/gene amplifications or PTEN deficient/PI3K overactive pathways.

Ezh2 Knockout in B Cells Impairs Plasmablast Differentiation and Ameliorates Lupus-like Disease in MRL/lpr Mice.

OBJECTIVES: EZH2 regulates B cell development and differentiation. We previously demonstrated increased EZH2 expression in peripheral blood mononuclear cells from lupus patients. The goal of this study was to evaluate the role of EZH2 expression in B cells in the pathogenesis of lupus. METHODS: We generated an MRL/lpr mouse with floxed Ezh2, which was crossed with CD19-Cre mice to examine the effect of B cell EZH2 deficiency in MRL/lpr lupus-prone mice. Differentiation of B cells was assessed using flow cytometry. Single-cell RNA sequencing and single-cell B cell receptor sequencing were performed. In vitro B cell culture with an X-box binding protein 1 (XBP1) inhibitor was performed. EZH2 and XBP1 messenger RNA levels in CD19+ B cells isolated from lupus patients and healthy controls were analyzed. RESULTS: We show that Ezh2 deletion in B cells significantly decreased autoantibody production and improved glomerulonephritis. B cell development was altered in the bone marrow and spleen of EZH2-deficient mice. Differentiation of germinal center B cells and plasmablasts was impaired. Single-cell RNA sequencing showed that XBP1, a key transcription factor in B cell development, is down-regulated in the absence of EZH2. Inhibiting XBP1 in vitro impairs plasmablast development similar to EZH2 deficiency in mice. Single-cell B cell receptor RNA sequencing revealed defective immunoglobulin class-switch recombination in EZH2-deficient mice. In human lupus B cells, we observed a strong correlation between EZH2 and XBP1 messenger RNA expression levels. CONCLUSION: EZH2 overexpression in B cells contributes to disease pathogenesis in lupus.

Sphingosine Kinase 2 in Stromal Fibroblasts Creates a Hospitable Tumor Microenvironment in Breast Cancer.

Reciprocal interactions between breast cancer cells and the tumor microenvironment (TME) are important for cancer progression and metastasis. We report here that the deletion or inhibition of sphingosine kinase 2 (SphK2), which produces sphingosine-1-phosphate (S1P), markedly suppresses syngeneic breast tumor growth and lung metastasis in mice by creating a hostile microenvironment for tumor growth and invasion. SphK2 deficiency decreased S1P and concomitantly increased ceramides, including C16-ceramide, in stromal fibroblasts. Ceramide accumulation suppressed activation of cancer-associated fibroblasts (CAF) by upregulating stromal p53, which restrained production of tumor-promoting factors to reprogram the TME and to restrict breast cancer establishment. Ablation of p53 in SphK2-deficient fibroblasts reversed these effects, enabled CAF activation and promoted tumor growth and invasion. These data uncovered a novel role of SphK2 in regulating non-cell-autonomous functions of p53 in stromal fibroblasts and their transition to tumor-promoting CAFs, paving the way for the development of a strategy to target the TME and to enhance therapeutic efficacy. SIGNIFICANCE: Sphingosine kinase 2 (SphK2) facilitates the activation of stromal fibroblasts to tumor-promoting cancer-associated fibroblasts by suppressing host p53 activity, revealing SphK2 as a potential target to reprogram the TME.

Astrocyte Elevated Gene-1 Cys75 S-Palmitoylation by ZDHHC6 Regulates Its Biological Activity.

Nonalcoholic fatty liver disease is a major risk factor for hepatocellular carcinoma (HCC). Astrocyte elevated gene-1/Metadherin (AEG-1/MTDH) augments lipid accumulation (steatosis), inflammation, and tumorigenesis, thereby promoting the whole spectrum of this disease process. Targeting AEG-1 is a potential interventional strategy for nonalcoholic steatohepatitis (NASH) and HCC. Thus, proper understanding of the regulation of this molecule is essential. We found that AEG-1 is palmitoylated at residue cysteine 75 (Cys75). Mutation of Cys75 to serine (Ser) completely abolished AEG-1 palmitoylation. We identified ZDHHC6 as a palmitoyltransferase catalyzing the process in HEK293T cells. To obtain insight into how palmitoylation regulates AEG-1 function, we generated knock-in mice by CRISPR/Cas9 in which Cys75 of AEG-1 was mutated to Ser (AEG-1-C75S). No developmental or anatomical abnormality was observed between AEG-1-wild type (AEG-1-WT) and AEG-1-C75S littermates. However, global gene expression analysis by RNA-sequencing unraveled that signaling pathways and upstream regulators, which contribute to cell proliferation, motility, inflammation, angiogenesis, and lipid accumulation, were activated in AEG-1-C75S hepatocytes compared to AEG-1-WT. These findings suggest that AEG-1-C75S functions as dominant positive and that palmitoylation restricts oncogenic and NASH-promoting functions of AEG-1. We thus identify a previously unknown regulatory mechanism of AEG-1, which might help design new therapeutic strategies for NASH and HCC.

Melanoma differentiation associated gene-9/syndecan binding protein promotes hepatocellular carcinoma.

BACKGROUND AND AIMS: The oncogene Melanoma differentiation associated gene-9/syndecan binding protein (MDA-9/SDCBP) is overexpressed in many cancers, promoting aggressive, metastatic disease. However, the role of MDA-9 in regulating hepatocellular carcinoma (HCC) has not been well studied. APPROACH AND RESULTS: To unravel the function of MDA-9 in HCC, we generated and characterized a transgenic mouse with hepatocyte-specific overexpression of MDA-9 (Alb/MDA-9). Compared with wild-type (WT) littermates, Alb/MDA-9 mice demonstrated significantly higher incidence of N-nitrosodiethylamine/phenobarbital-induced HCC, with marked activation and infiltration of macrophages. RNA sequencing (RNA-seq) in naive WT and Alb/MDA-9 hepatocytes identified activation of signaling pathways associated with invasion, angiogenesis, and inflammation, especially NF-κB and integrin-linked kinase signaling pathways. In nonparenchymal cells purified from naive livers, single-cell RNA-seq showed activation of Kupffer cells and macrophages in Alb/MDA-9 mice versus WT mice. A robust increase in the expression of Secreted phosphoprotein 1 (Spp1/osteopontin) was observed upon overexpression of MDA-9. Inhibition of NF-κB pathway blocked MDA-9-induced Spp1 induction, and knock down of Spp1 resulted in inhibition of MDA-9-induced macrophage migration, as well as angiogenesis. CONCLUSIONS: Alb/MDA-9 is a mouse model with MDA-9 overexpression in any tissue type. Our findings unravel an HCC-promoting role of MDA-9 mediated by NF-κB and Spp1 and support the rationale of using MDA-9 inhibitors as a potential treatment for aggressive HCC.