skip to main content


Title: Integrative genomic analysis of mouse and human hepatocellular carcinoma

Cancer genomics has enabled the exhaustive molecular characterization of tumors and exposed hepatocellular carcinoma (HCC) as among the most complex cancers. This complexity is paralleled by dozens of mouse models that generate histologically similar tumors but have not been systematically validated at the molecular level. Accurate models of the molecular pathogenesis of HCC are essential for biomedical progress; therefore we compared genomic and transcriptomic profiles of four separate mouse models [MUP transgenic, TAK1-knockout, carcinogen-driven diethylnitrosamine (DEN), and Stelic Animal Model (STAM)] with those of 987 HCC patients with distinct etiologies. These four models differed substantially in their mutational load, mutational signatures, affected genes and pathways, and transcriptomes. STAM tumors were most molecularly similar to human HCC, with frequent mutations inCtnnb1, similar pathway alterations, and high transcriptomic similarity to high-grade, proliferative human tumors with poor prognosis. In contrast, TAK1 tumors better reflected the mutational signature of human HCC and were transcriptionally similar to low-grade human tumors. DEN tumors were least similar to human disease and almost universally carried theBrafV637E mutation, which is rarely found in human HCC. Immune analysis revealed that strain-specific MHC-I genotype can influence the molecular makeup of murine tumors. Thus, different mouse models of HCC recapitulate distinct aspects of HCC biology, and their use should be adapted to specific questions based on the molecular features provided here.

 
more » « less
NSF-PAR ID:
10076795
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
Proceedings of the National Academy of Sciences
Date Published:
Journal Name:
Proceedings of the National Academy of Sciences
Volume:
115
Issue:
42
ISSN:
0027-8424
Page Range / eLocation ID:
p. E9879-E9888
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. INTRODUCTION Balance between excitatory and inhibitory neuron (interneuron) populations in the cortex promotes normal brain function. Interneurons are primarily generated in the medial, caudal, and lateral ganglionic eminences (MGE, CGE, and LGE) of the ventral embryonic forebrain; these subregions give rise to distinct interneuron subpopulations. In rodents, the MGE generates cortical interneurons, the parvalbumin + (PV + ) and somatostatin + (SST + ) subtypes that connect with excitatory neurons to regulate their activity. Defects in interneuron production have been implicated in neurodevelopmental and psychiatric disorders including autism, epilepsy, and schizophrenia. RATIONALE How does the human MGE (hMGE) produce the number of interneurons required to populate the forebrain? The hMGE contains progenitor clusters distinct from what has been observed in the rodent MGE and other germinal zones of the human brain. This cytoarchitecture could be the key to understanding interneuron neurogenesis. We investigated the cellular and molecular properties of different compartments within the developing hMGE, from 14 gestational weeks (GW) to 39 GW (term), to study their contribution to the production of inhibitory interneurons. We developed a xenotransplantation assay to follow the migration and maturation of the human interneurons derived from this germinal region. RESULTS Within the hMGE, densely packed aggregates (nests) of doublecortin + (DCX + ) and LHX6 + cells were surrounded by nestin + progenitor cells and their processes. These DCX + cell–enriched nests (DENs) were observed in the hMGE but not in the adjacent LGE. We found that cells within DENs expressed molecular markers associated with young neurons, such as DCX, and polysialylated neural cell adhesion molecule (PSA-NCAM). A subpopulation also expressed Ki-67, a marker of proliferation; therefore, we refer to these cells as neuroblasts. A fraction of DCX + cells inside DENs expressed SOX2 and E2F1, transcription factors associated with progenitor and proliferative properties. More than 20% of DCX + cells in the hMGE were dividing, specifically within DENs. Proliferating neuroblasts in DENs persisted in the hMGE throughout prenatal human brain development. The division of DCX + cells was confirmed by transmission electron microscopy and time-lapse microscopy. Electron microscopy revealed adhesion contacts between cells within DENs, providing multiple sites to anchor DEN cells together. Neuroblasts within DENs express PCDH19, and nestin + progenitors surrounding DENs express PCDH10; these findings suggest a role for differential cell adhesion in DEN formation and maintenance. When transplanted into the neonatal mouse brain, dissociated hMGE cells reformed DENs containing proliferative DCX + cells, similar to DENs observed in the prenatal human brain. This suggests that DENs are generated by cell-autonomous mechanisms. In addition to forming DENs, transplanted hMGE-derived neuroblasts generated young neurons that migrated extensively into cortical and subcortical regions in the host mouse brain. One year after transplantation, these neuroblasts had differentiated into distinct γ-aminobutyric acid–expressing (GABAergic) interneuron subtypes, including SST + and PV + cells, that showed morphological and functional maturation. CONCLUSION The hMGE harbors DENs, where cells expressing early neuronal markers continue to divide and produce GABAergic interneurons. This MGE-specific arrangement of neuroblasts in the human brain is present until birth, supporting expanded neurogenesis for inhibitory neurons. Given the robust neurogenic output from this region, knowledge of the mechanisms underlying cortical interneuron production in the hMGE will provide insights into the cell types and developmental periods that are most vulnerable to genetic or environmental insults. Nests of DCX + cells in the ventral prenatal brain. Schematic of a coronal view of the embryonic human forebrain showing the medial ganglionic eminence (MGE, green), with nests of DCX + cells (DENs, green). Nestin + progenitor cells (blue) are present within the VZ and iSVZ and are intercalated in the oSVZ (where DENs reside). The initial segment of the oSVZ contains palisades of nestin + progenitors referred to as type I clusters (light blue cells) around DENs. In the outer part of the oSVZ, DENs transition to chains of migrating DCX + cells; surrounding nestin + progenitors are arranged into groups of cells referred to as type II clusters (white cells). In addition to proliferation of nestin + progenitors, cell division is present among DCX + cells within DENs, suggesting multiple progenitor states for the generation of MGE-derived interneurons in the human forebrain. ILLUSTRATION: NOEL SIRIVANSANTI 
    more » « less
  2. Abstract The microenvironment that surrounds pancreatic ductal adenocarcinoma (PDAC) is profoundly desmoplastic and immunosuppressive. Understanding triggers of immunosuppression during the process of pancreatic tumorigenesis would aid in establishing targets for effective prevention and therapy. Here, we interrogated differential molecular mechanisms dependent on cell of origin and subtype that promote immunosuppression during PDAC initiation and in established tumors. Transcriptomic analysis of cell-of-origin–dependent epithelial gene signatures revealed that Nt5e/CD73, a cell-surface enzyme required for extracellular adenosine generation, is one of the top 10% of genes overexpressed in murine tumors arising from the ductal pancreatic epithelium as opposed to those rising from acinar cells. These findings were confirmed by IHC and high-performance liquid chromatography. Analysis in human PDAC subtypes indicated that high Nt5e in murine ductal PDAC models overlaps with high NT5E in human PDAC squamous and basal subtypes, considered to have the highest immunosuppression and worst prognosis. Multiplex immunofluorescent analysis showed that activated CD8+ T cells in the PDAC tumor microenvironment express high levels of CD73, indicating an opportunity for immunotherapeutic targeting. Delivery of CD73 small-molecule inhibitors through various delivery routes reduced tumor development and growth in genetically engineered and syngeneic mouse models. In addition, the adenosine receptor Adora2b was a determinant of adenosine-mediated immunosuppression in PDAC. These findings highlight a molecular trigger of the immunosuppressive PDAC microenvironment elevated in the ductal cell of origin, linking biology with subtype classification, critical components for PDAC immunoprevention and personalized approaches for immunotherapeutic intervention. Significance: Ductal-derived pancreatic tumors have elevated epithelial and CD8+GZM+ T-cell CD73 expression that confers sensitivity to small-molecule inhibition of CD73 or Adora2b to promote CD8+ T-cell–mediated tumor regression. See related commentary by DelGiorno, p. 977 
    more » « less
  3. Hepatocellular carcinoma (HCC) is one of the most fatal cancers in the world. There is an urgent need to understand the molecular background of HCC to facilitate the identification of biomarkers and discover effective therapeutic targets. Published transcriptomic studies have reported a large number of genes that are individually significant for HCC. However, reliable biomarkers remain to be determined. In this study, built on max-linear competing risk factor models, we developed a machine learning analytical framework to analyze transcriptomic data to identify the most miniature set of differentially expressed genes (DEGs). By analyzing 9 public whole-transcriptome datasets (containing 1184 HCC samples and 672 nontumor controls), we identified 5 critical differentially expressed genes (DEGs) (ie, CCDC107, CXCL12, GIGYF1, GMNN, and IFFO1) between HCC and control samples. The classifiers built on these 5 DEGs reached nearly perfect performance in identification of HCC. The performance of the 5 DEGs was further validated in a US Caucasian cohort that we collected (containing 17 HCC with paired nontumor tissue). The conceptual advance of our work lies in modeling gene-gene interactions and correcting batch effect in the analytic framework. The classifiers built on the 5 DEGs demonstrated clear signature patterns for HCC. The results are interpretable, robust, and reproducible across diverse cohorts/populations with various disease etiologies, indicating the 5 DEGs are intrinsic variables that can describe the overall features of HCC at the genomic level. The analytical framework applied in this study may pave a new way for improving transcriptome profiling analysis of human cancers.

     
    more » « less
  4. Abstract

    Lineage plasticity has been proposed as a major source of intratumoral heterogeneity and therapeutic resistance. Here, by employing an inducible genetic engineered mouse model, we illustrate that lineage plasticity enables advanced Pancreatic Ductal Adenocarcinoma (PDAC) tumors to develop spontaneous relapse following elimination of the central oncogenic driver - Yap. Transcriptomic and immunohistochemistry analysis of a large panel of PDAC tumors reveals that within high-grade tumors, small niches of PDAC cells gradually evolve to re-activate pluripotent transcription factors (PTFs), which lessen their dependency on Yap. Comprehensive Cut&Tag analysis demonstrate that although acquisition of PTF expression is coupled with the process of epithelial-to-mesenchymal transition (EMT), PTFs form a core transcriptional regulatory circuitry (CRC) with Jun to overcome Yap dependency, which is distinct from the classic TGFb-induced EMT-TF network. A chemical-genetic screen and follow-up functional studies establish Brd4 as an epigenetic gatekeeper for the PTF-Jun CRC, and strong synergy between BET and Yap inhibitors in blocking PDAC growth.

     
    more » « less
  5. Abstract Motivation Somatic mutations result from processes related to DNA replication or environmental/lifestyle exposures. Knowing the activity of mutational processes in a tumor can inform personalized therapies, early detection, and understanding of tumorigenesis. Computational methods have revealed 30 validated signatures of mutational processes active in human cancers, where each signature is a pattern of single base substitutions. However, half of these signatures have no known etiology, and some similar signatures have distinct etiologies, making patterns of mutation signature activity hard to interpret. Existing mutation signature detection methods do not consider tumor-level clinical/demographic (e.g. smoking history) or molecular features (e.g. inactivations to DNA damage repair genes). Results To begin to address these challenges, we present the Tumor Covariate Signature Model (TCSM), the first method to directly model the effect of observed tumor-level covariates on mutation signatures. To this end, our model uses methods from Bayesian topic modeling to change the prior distribution on signature exposure conditioned on a tumor’s observed covariates. We also introduce methods for imputing covariates in held-out data and for evaluating the statistical significance of signature-covariate associations. On simulated and real data, we find that TCSM outperforms both non-negative matrix factorization and topic modeling-based approaches, particularly in recovering the ground truth exposure to similar signatures. We then use TCSM to discover five mutation signatures in breast cancer and predict homologous recombination repair deficiency in held-out tumors. We also discover four signatures in a combined melanoma and lung cancer cohort—using cancer type as a covariate—and provide statistical evidence to support earlier claims that three lung cancers from The Cancer Genome Atlas are misdiagnosed metastatic melanomas. Availability and implementation TCSM is implemented in Python 3 and available at https://github.com/lrgr/tcsm, along with a data workflow for reproducing the experiments in the paper. Supplementary information Supplementary data are available at Bioinformatics online. 
    more » « less