Macrophages in the infarcted heart acquire a fibrogenic phenotype, expressing matricellular proteins, but do not undergo fibroblast conversion.

Although some studies have suggested that macrophages may secrete structural collagens, and convert to fibroblast-like cells, macrophage to fibroblast transdifferentiation in infarcted and remodeling hearts remains controversial. Our study uses linage tracing approaches and single cell transcriptomics to examine whether macrophages undergo fibroblast conversion, and to characterize the extracellular matrix expression profile of myeloid cells in myocardial infarction. To examine whether infarct macrophages undergo fibroblast conversion, we identified macrophage-derived progeny using the inducible CX3CR1CreER mice crossed with the PDGFRαEGFP reporter line for reliable fibroblast identification. The abundant fibroblasts that infiltrated the infarcted myocardium after 7 and 28 days of coronary occlusion were not derived from CX3CR1+ macrophages. Infarct macrophages retained myeloid cell characteristics and did not undergo conversion to myofibroblasts, endothelial or vascular mural cells. Single cell RNA-seq of CSF1R+ myeloid cells harvested from control and infarcted hearts showed no significant expression of fibroblast identity genes by myeloid cell clusters. Moreover, infarct macrophages did not express significant levels of genes encoding structural collagens. However, infarct macrophage and monocyte clusters were the predominant source of the fibrogenic growth factors Tgfb1 and Pdgfb, and of the matricellular proteins Spp1/Osteopontin, Thbs1/Thrombospondin-1, Emilin2, and Fn1/fibronectin, while expressing significant amounts of several other matrix genes, including Vcan/versican, Ecm1 and Sparc. ScRNA-seq data suggested similar patterns of matrix gene expression in human myocardial infarction. In conclusion, infarct macrophages do not undergo fibroblast or myofibroblast conversion and do not exhibit upregulation of structural collagens but may contribute to fibrotic remodeling by producing several fibrogenic matricellular proteins.

VISTA Emerges as a Promising Target against Immune Evasion Mechanisms in Medulloblastoma.

BACKGROUND: Relapsed medulloblastoma (MB) poses a significant therapeutic challenge due to its highly immunosuppressive tumor microenvironment. Immune checkpoint inhibitors (ICIs) have struggled to mitigate this challenge, largely due to low T-cell infiltration and minimal PD-L1 expression. Identifying the mechanisms driving low T-cell infiltration is crucial for developing more effective immunotherapies. METHODS: We utilize a syngeneic mouse model to investigate the tumor immune microenvironment of MB and compare our findings to transcriptomic and proteomic data from human MB. RESULTS: Flow cytometry reveals a notable presence of CD45hi/CD11bhi macrophage-like and CD45int/CD11bint microglia-like tumor-associated macrophages (TAMs), alongside regulatory T-cells (Tregs), expressing high levels of the inhibitory checkpoint molecule VISTA. Compared to sham control mice, the CD45hi/CD11bhi compartment significantly expands in tumor-bearing mice and exhibits a myeloid-specific signature composed of VISTA, CD80, PD-L1, CTLA-4, MHCII, CD40, and CD68. These findings are corroborated by proteomic and transcriptomic analyses of human MB samples. Immunohistochemistry highlights an abundance of VISTA-expressing myeloid cells clustering at the tumor-cerebellar border, while T-cells are scarce and express FOXP3. Additionally, tumor cells exhibit immunosuppressive properties, inhibiting CD4 T-cell proliferation in vitro. Identification of VISTA's binding partner, VSIG8, on tumor cells, and its correlation with increased VISTA expression in human transcriptomic analyses suggests a potential therapeutic target. CONCLUSIONS: This study underscores the multifaceted mechanisms of immune evasion in MB and highlights the therapeutic potential of targeting the VISTA-VSIG axis to enhance anti-tumor responses.

Molecular and network disruptions in neurodevelopment uncovered by single cell transcriptomics analysis of CHD8 heterozygous cerebral organoids.

More than 100 genes have been associated with significantly increased risks of autism spectrum disorders (ASD) with an estimate of ∼1000 genes that may contribute. The new challenge is to investigate the molecular and cellular functions of these genes during neural and brain development, and then even more challenging, to link the altered molecular and cellular phenotypes to the ASD clinical manifestations. In this study, we used single-cell RNA-seq analysis to study one of the top risk genes, CHD8, in cerebral organoids, which models early neural development. We identified 21 cell clusters in the organoid samples, representing non-neuronal cells, neural progenitors, and early differentiating neurons at the start of neural cell fate commitment. Comparisons of the cells with one copy of a CHD8 knockout allele, generated by CRISPR/Cas9 editing, and their isogenic controls uncovered thousands of differentially expressed genes, which were enriched with functions related to neural and brain development, cilium organization, and extracellular matrix organization. The affected genes were also enriched with genes and pathways previously implicated in ASD, but surprisingly not for schizophrenia and intellectual disability risk genes. The comparisons also uncovered cell composition changes, indicating potentially altered neural differential trajectories upon CHD8 reduction. Moreover, we found that cell-cell communications were affected in the CHD8 knockout organoids, including the interactions between neural and glial cells. Taken together, our results provide new data and information for understanding CHD8 functions in the early stages of neural lineage development and interaction.

Aging limits stemness and tumorigenesis in the lung by reprogramming iron homeostasis.

Aging is associated with a decline in the number and fitness of adult stem cells 1-4 . Aging-associated loss of stemness is posited to suppress tumorigenesis 5,6 , but this hypothesis has not been tested in vivo . Here, using physiologically aged autochthonous genetically engineered mouse models and primary cells 7,8 , we demonstrate aging suppresses lung cancer initiation and progression by degrading stemness of the alveolar cell of origin. This phenotype is underpinned by aging-associated induction of the transcription factor NUPR1 and its downstream target lipocalin-2 in the cell of origin in mice and humans, leading to a functional iron insufficiency in the aged cells. Genetic inactivation of the NUPR1-lipocalin-2 axis or iron supplementation rescue stemness and promote tumorigenic potential of aged alveolar cells. Conversely, targeting the NUPR1- lipocalin-2 axis is detrimental to young alveolar cells via induction of ferroptosis. We find that aging-associated DNA hypomethylation at specific enhancer sites associates with elevated NUPR1 expression, which is recapitulated in young alveolar cells by inhibition of DNA methylation. We uncover that aging drives a functional iron insufficiency, which leads to loss of stemness and tumorigenesis, but promotes resistance to ferroptosis. These findings have significant implications for the therapeutic modulation of cellular iron homeostasis in regenerative medicine and in cancer prevention. Furthermore, our findings are consistent with a model whereby most human cancers initiate in young individuals, revealing a critical window for such cancer prevention efforts.

TAD hierarchy restricts poised LTR activation and loss of TAD hierarchy promotes LTR co-option in cancer.

Transposable elements (TEs) are abundant in the human genome, and they provide the sources for genetic and functional diversity. The regulation of TEs expression and their functional consequences in physiological conditions and cancer development remain to be fully elucidated. Previous studies suggested TEs are repressed by DNA methylation and chromatin modifications. The effect of 3D chromatin topology on TE regulation remains elusive. Here, by integrating transcriptome and 3D genome architecture studies, we showed that haploinsufficient loss of NIPBL selectively activates alternative promoters at the long terminal repeats (LTRs) of the TE subclasses. This activation occurs through the reorganization of topologically associating domain (TAD) hierarchical structures and recruitment of proximal enhancers. These observations indicate that TAD hierarchy restricts transcriptional activation of LTRs that already possess open chromatin features. In cancer, perturbation of the hierarchical chromatin topology can lead to co-option of LTRs as functional alternative promoters in a context-dependent manner and drive aberrant transcriptional activation of novel oncogenes and other divergent transcripts. These data uncovered a new layer of regulatory mechanism of TE expression beyond DNA and chromatin modification in human genome. They also posit the TAD hierarchy dysregulation as a novel mechanism for alternative promoter-mediated oncogene activation and transcriptional diversity in cancer, which may be exploited therapeutically.

Comprehensive single cell transcriptomics analysis of murine osteosarcoma uncovers Skp2 function in metastasis, genomic instability and immune activation and reveals additional target pathways.

Osteosarcoma (OS) is the most common primary pediatric bone malignancy. One promising new therapeutic target is SKP2, encoding a substrate recognition factor of the SCF E3 ubiquitin ligase responsible for ubiquitination and proteasome degradation of substrate p27, thus driving cellular proliferation. We have shown previously that knockout of Skp2 in an immunocompetent transgenic mouse model of OS improved survival, drove apoptosis, and induced tumor inflammation. Here, we applied single-cell RNA-sequencing (scRNA-seq) to study primary OS tumors derived from Osx-Cre driven conditional knockout of Rb1 and Trp53. We showed that murine OS models recapitulate the tumor heterogeneity and microenvironment complexity observed in patient tumors. We further compared this model with OS models with functional disruption of Skp2: one with Skp2 knockout and the other with the Skp2-p27 interaction disrupted (resulting in p27 overexpression). We found reduction of T cell exhaustion and upregulation of interferon activation, along with evidence of replicative and endoplasmic reticulum-related stress in the Skp2 disruption models, and showed that interferon induction was correlated with improved survival in OS patients. Additionally, our scRNA-seq analysis uncovered decreased activities of metastasis-related gene signatures in the Skp2-disrupted OS, which we validated by observation of a strong reduction in lung metastasis in the Skp2 knockout mice. Finally, we report several potential mechanisms of escape from targeting Skp2 in OS, including upregulation of Myc targets, DNA copy number amplification and overexpression of alternative E3 ligase genes, and potential alternative lineage activation. These mechanistic insights into OS tumor biology and Skp2 function suggest novel targets for new, synergistic therapies, while the data and our comprehensive analysis may serve as a public resource for further big data-driven OS research.

Structure and transcription of integrated HPV DNA in vulvar carcinomas.

HPV infections are associated with a fraction of vulvar cancers. Through hybridization capture and DNA sequencing, HPV DNA was detected in five of thirteen vulvar cancers. HPV16 DNA was integrated into human DNA in three of the five. The insertions were in introns of human NCKAP1, C5orf67, and LRP1B. Integrations in NCKAP1 and C5orf67 were flanked by short direct repeats in the human DNA, consistent with HPV DNA insertions at sites of abortive, staggered, endonucleolytic incisions. The insertion in C5orf67 was present as a 36 kbp, human-HPV-hetero-catemeric DNA as either an extrachromosomal circle or a tandem repeat within the human genome. The human circularization/repeat junction was defined at single nucleotide resolution. The integrated viral DNA segments all retained an intact upstream regulatory region and the adjacent viral E6 and E7 oncogenes. RNA sequencing revealed that the only HPV genes consistently transcribed from the integrated viral DNAs were E7 and E6*I. The other two HPV DNA+ tumors had coinfections, but no evidence for integration. HPV-positive and HPV-negative vulvar cancers exhibited contrasting human, global gene expression patterns partially overlapping with previously observed differences between HPV-positive and HPV-negative cervical and oropharyngeal cancers. A substantial fraction of the differentially expressed genes involved immune system function. Thus, transcription and HPV DNA integration in vulvar cancers resemble those in other HPV-positive cancers. This study emphasizes the power of hybridization capture coupled with DNA and RNA sequencing to identify a broad spectrum of HPV types, determine human genome integration status of viral DNAs, and elucidate their structures.

Fine-tuning spatial-temporal dynamics and surface receptor expression support plasma cell-intrinsic longevity.

Durable serological memory following vaccination is critically dependent on the production and survival of long-lived plasma cells (LLPCs). Yet, the factors that control LLPC specification and survival remain poorly resolved. Using intravital two-photon imaging, we find that in contrast to most plasma cells (PCs) in the bone marrow (BM), LLPCs are uniquely sessile and organized into clusters that are dependent on APRIL, an important survival factor. Using deep, bulk RNA sequencing, and surface protein flow-based phenotyping, we find that LLPCs express a unique transcriptome and phenotype compared to bulk PCs, fine-tuning expression of key cell surface molecules, CD93, CD81, CXCR4, CD326, CD44, and CD48, important for adhesion and homing. Conditional deletion of Cxcr4 in PCs following immunization leads to rapid mobilization from the BM, reduced survival of antigen-specific PCs, and ultimately accelerated decay of antibody titer. In naïve mice, the endogenous LLPCs BCR repertoire exhibits reduced diversity, reduced somatic mutations, and increased public clones and IgM isotypes, particularly in young mice, suggesting LLPC specification is non-random. As mice age, the BM PC compartment becomes enriched in LLPCs, which may outcompete and limit entry of new PCs into the LLPC niche and pool.

scDAPP: a comprehensive single-cell transcriptomics analysis pipeline optimized for cross-group comparison.

Single-cell transcriptomics profiling has increasingly been used to evaluate cross-group differences in cell population and cell-type gene expression. This often leads to large datasets with complex experimental designs that need advanced comparative analysis. Concurrently, bioinformatics software and analytic approaches also become more diverse and constantly undergo improvement. Thus, there is an increased need for automated and standardized data processing and analysis pipelines, which should be efficient and flexible too. To address these, we develop the single-cell Differential Analysis and Processing Pipeline (scDAPP), a R-based workflow for comparative analysis of single cell (or nucleus) transcriptomic data between two or more groups and at the levels of single cells or "pseudobulking" samples. The pipeline automates many steps of pre-processing using data-learnt parameters, uses previously benchmarked software, and generates comprehensive intermediate data and final results that are valuable for both beginners and experts of scRNA-seq analysis. Moreover, the analytic reports, augmented by extensive data visualization, increase the transparency of computational analysis and parameter choices, while facilitate users to go seamlessly from raw data to biological interpretation. Availability and Implementation: scDAPP is freely available for non-commercial usage as an R package under the MIT license. Source code, documentation and sample data are available at the GitHub (https://github.com/bioinfoDZ/scDAPP).

TOP CAR with TMIGD2 as a safe and effective costimulatory domain in CAR cells treating human solid tumors.

Chimeric antigen receptor (CAR)-T cell therapy shows impressive efficacy treating hematologic malignancies but requires further optimization in solid tumors. Here, we developed a TMIGD2 optimized potent/persistent (TOP) CAR that incorporated the costimulatory domain of TMIGD2, a T and NK cell costimulator, and monoclonal antibodies targeting the IgV domain of B7-H3, an immune checkpoint expressed on solid tumors and tumor vasculature. Comparing second- and third-generation B7-H3 CARs containing TMIGD2, CD28, and/or 4-1BB costimulatory domains revealed superior antitumor responses in B7-H3.TMIGD2 and B7-H3.CD28.4-1BB CAR-T cells in vitro. Comparing these two constructs using in vivo orthotopic human cancer models demonstrated that B7-H3.TMIGD2 CAR-T cells had equivalent or superior antitumor activity, survival, expansion, and persistence. Mechanistically, B7-H3.TMIGD2 CAR-T cells maintained mitochondrial metabolism; produced less cytokines; and established fewer exhausted cells, more central memory cells, and a larger CD8/CD4 T cell ratio. These studies demonstrate that the TOP CAR with TMIGD2 costimulation offered distinct benefits from CD28.41BB costimulation and is effective against solid tumors.

KDM5-mediated transcriptional activation of ribosomal protein genes alters translation efficiency to regulate mitochondrial metabolism in neurons.

Genes encoding the KDM5 family of transcriptional regulators are disrupted in individuals with intellectual disability (ID). To understand the link between KDM5 and ID, we characterized five Drosophila strains harboring missense alleles analogous to those observed in patients. These alleles disrupted neuroanatomical development, cognition and other behaviors, and displayed a transcriptional signature characterized by the downregulation of many ribosomal protein genes. A similar transcriptional profile was observed in KDM5C knockout iPSC-induced human glutamatergic neurons, suggesting an evolutionarily conserved role for KDM5 proteins in regulating this class of gene. In Drosophila, reducing KDM5 changed neuronal ribosome composition, lowered the translation efficiency of mRNAs required for mitochondrial function, and altered mitochondrial metabolism. These data highlight the cellular consequences of altered KDM5-regulated transcriptional programs that could contribute to cognitive and behavioral phenotypes. Moreover, they suggest that KDM5 may be part of a broader network of proteins that influence cognition by regulating protein synthesis.

Metacell-based differential expression analysis identifies cell type specific temporal gene response programs in COVID-19 patient PBMCs.

By profiling gene expression in individual cells, single-cell RNA-sequencing (scRNA-seq) can resolve cellular heterogeneity and cell-type gene expression dynamics. Its application to time-series samples can identify temporal gene programs active in different cell types, for example, immune cells' responses to viral infection. However, current scRNA-seq analysis has limitations. One is the low number of genes detected per cell. The second is insufficient replicates (often 1-2) due to high experimental cost. The third lies in the data analysis-treating individual cells as independent measurements leads to inflated statistics. To address these, we explore a new computational framework, specifically whether "metacells" constructed to maintain cellular heterogeneity within individual cell types (or clusters) can be used as "replicates" for increasing statistical rigor. Toward this, we applied SEACells to a time-series scRNA-seq dataset from peripheral blood mononuclear cells (PBMCs) after SARS-CoV-2 infection to construct metacells, and used them in maSigPro for quadratic regression to find significantly differentially expressed genes (DEGs) over time, followed by clustering expression velocity trends. We showed that such metacells retained greater expression variances and produced more biologically meaningful DEGs compared to either metacells generated randomly or from simple pseudobulk methods. More specifically, this approach correctly identified the known ISG15 interferon response program in almost all PBMC cell types and many DEGs enriched in the previously defined SARS-CoV-2 infection response pathway. It also uncovered additional and more cell type-specific temporal gene expression programs. Overall, our results demonstrate that the metacell-pseudoreplicate strategy could potentially overcome the limitation of 1-2 replicates.

Lineage commitment pathways epigenetically oppose oncogenic Gαq/11-YAP signaling in dormant disseminated uveal melanoma.

The mechanisms driving late relapse in uveal melanoma (UM) patients remains a medical mystery and major challenge. Clinically it is inferred that UM disseminated cancer cells (DCCs) persist asymptomatic for years-to-decades mainly in the liver before they manifest as symptomatic metastasis. Here we reveal using Gαq/11 mut /BAP wt human uveal melanoma models and human UM metastatic samples, that the neural crest lineage commitment nuclear receptor NR2F1 is a key regulator of spontaneous UM DCC dormancy in the liver. Using a quiescence reporter, RNA-seq and multiplex imaging we revealed that rare dormant UM DCCs upregulate NR2F1 expression and genes related to neural crest programs while repressing gene related to cell cycle progression. Gain and loss of function assays showed that NR2F1 silences YAP1/TEAD1 transcription downstream of Gαq/11 signaling and that NR2F1 expression can also be repressed by YAP1. YAP1 expression is repressed by NR2F1 binding to its promoter and changing the histone H3 tail activation marks to repress YAP1 transcription. In vivo CRISPR KO of NR2F1 led dormant UM DCCs to awaken and initiate relentless liver metastatic growth. Cut&Run and bulk RNA sequencing further confirmed that NR2F1 epigenetically stimulates neuron axon guidance and neural lineage programs, and it globally represses gene expression linked to G-protein signaling to drive dormancy. Pharmacological inhibition of Gαq/11 mut signaling resulted in NR2F1 upregulation and robust UM growth arrest, which was also achieved using a novel NR2F1 agonist. Our work sheds light on the molecular underpinnings of UM dormancy revealing that transcriptional programs driven by NR2F1 epigenetically short-circuit Gαq/11 signaling to its downstream target YAP1. Highlights: Quiescent solitary uveal melanoma (UM) DCCs in the liver up- and down-regulate neural crest and cell cycle progression programs, respectively.NR2F1 drives solitary UM DCC dormancy by antagonizing the Gαq/11-YAP1 pathway; small molecule Gαq/11 inhibition restores NR2F1 expression and quiescence. NR2F1 short-circuits oncogenic YAP1 and G-protein signaling via a chromatin remodeling program. Loss of function of NR2F1 in dormant UM DCCs leads to aggressive liver metastasis.

Targeted inhibition of SCFSKP2 confers anti-tumor activities resulting in a survival benefit in osteosarcoma.

Osteosarcoma(OS) is a highly aggressive bone cancer for which treatment has remained essentially unchanged for decades. Although OS is characterized by extensive genomic heterogeneity and instability, RB1 and TP53 have been shown to be the most commonly inactivated tumor suppressors in OS. We previously generated a mouse model with a double knockout (DKO) of Rb1 and Trp53 within cells of the osteoblastic lineage, which largely recapitulates human OS with nearly complete penetrance. SKP2 is a repression target of pRb and serves as a substrate recruiting subunit of the SCFSKP2 complex. In addition, SKP2 plays a central role in regulating the cell cycle by ubiquitinating and promoting the degradation of p27. We previously reported the DKOAA transgenic model, which harbored a knock-in mutation in p27 that impaired its binding to SKP2. Here, we generated a novel p53-Rb1-SKP2 triple-knockout model (TKO) to examine SKP2 function and its potential as a therapeutic target in OS. First, we observed that OS tumorigenesis was significantly delayed in TKO mice and their overall survival was markedly improved. In addition, the loss of SKP2 also promoted an apoptotic microenvironment and reduced the stemness of DKO tumors. Furthermore, we found that small-molecule inhibitors of SKP2 exhibited anti-tumor activities in vivo and in OS organoids as well as synergistic effects when combined with a standard chemotherapeutic agent. Taken together, our results suggest that SKP2 inhibitors may reduce the stemness plasticity of OS and should be leveraged as next-generation adjuvants in this cancer.

TMIGD2 is an orchestrator and therapeutic target on human acute myeloid leukemia stem cells.

Acute myeloid leukemia (AML) is initiated and sustained by a hierarchy of leukemia stem cells (LSCs), and elimination of this cell population is required for curative therapies. Here we show that transmembrane and immunoglobulin domain containing 2 (TMIGD2), a recently discovered co-stimulatory immune receptor, is aberrantly expressed by human AML cells, and can be used to identify and enrich functional LSCs. We demonstrate that TMIGD2 is required for the development and maintenance of AML and self-renewal of LSCs but is not essential for normal hematopoiesis. Mechanistically, TMIGD2 promotes proliferation, blocks myeloid differentiation and increases cell-cycle of AML cells via an ERK1/2-p90RSK-CREB signaling axis. Targeting TMIGD2 signaling with anti-TMIGD2 monoclonal antibodies attenuates LSC self-renewal and reduces leukemia burden in AML patient-derived xenograft models but has negligible effect on normal hematopoietic stem/progenitor cells. Thus, our studies reveal the function of TMIGD2 in LSCs and provide a promising therapeutic strategy for AML.

Tenascin-C expressing touch dome keratinocytes exhibit characteristics of all epidermal lineages.

The touch dome (TD) keratinocytes are specialized epidermal cells that intimately associate with the light touch sensing Merkel cells (MCs). The TD keratinocytes function as a niche for the MCs and can induce de novo hair follicles upon stimulation; however, how the TD keratinocytes are maintained during homeostasis remains unclear. scRNA-seq identified a specific TD keratinocyte marker, Tenascin-C (TNC). Lineage tracing of Tnc-expressing TD keratinocytes revealed that these cells maintain themselves as an autonomous epidermal compartment and give rise to MCs upon injury. Molecular characterization uncovered that, while the transcriptional and chromatin landscape of the TD keratinocytes is remarkably similar to that of the interfollicular epidermal keratinocytes, it also shares certain molecular signatures with the hair follicle keratinocytes. Our study highlights that the TD keratinocytes in the adult skin have molecular characteristics of keratinocytes of diverse epidermal lineages.

Fine-tuning spatial-temporal dynamics and surface receptor expression support plasma cell-intrinsic longevity.

Durable serological memory following vaccination is critically dependent on the production and survival of long-lived plasma cells (LLPCs). Yet, the factors that control LLPC specification and survival remain poorly resolved. Using intra-vital two-photon imaging, we find that in contrast to most plasma cells in the bone marrow, LLPCs are uniquely sessile and organized into clusters that are dependent on April, an important survival factor. Using deep, bulk RNA sequencing, and surface protein flow-based phenotyping, we find that LLPCs express a unique transcriptome and proteome compared to bulk PCs, fine tuning expression of key cell surface molecules, CD93, CD81, CXCR4, CD326, CD44 and CD48, important for adhesion and homing, and phenotypically label LLPCs within mature PC pool. Conditional deletion of Cxcr4 in PCs following immunization leads to rapid mobilization from the BM, reduced survival of antigen-specific PCs, and ultimately accelerated decay of antibody titer. In naive mice, the endogenous LLPCs BCR repertoire exhibits reduced diversity, reduced somatic mutations, and increased public clones and IgM isotypes, particularly in young mice, suggesting LLPC specification is non-random. As mice age, the BM PC compartment becomes enriched in LLPCs, which may outcompete and limit entry of new PC into the LLPC niche and pool.

SKP2 Knockout in Rb1/p53-Deficient Mouse Models of Osteosarcoma Induces Immune Infiltration and Drives a Transcriptional Program with a Favorable Prognosis.

Osteosarcoma is an aggressive bone malignancy with a poor prognosis. One putative proto-oncogene in osteosarcoma is SKP2, encoding a substrate recognition factor of the SCF E3 ubiquitin ligase. We previously demonstrated that Skp2 knockout in murine osteosarcoma improved survival and delayed tumorigenesis. Here, we performed RNA sequencing (RNA-seq) on tumors from a transgenic osteosarcoma mouse model with conditional Trp53 and Rb1 knockouts in the osteoblast lineage ("DKO": Osx1-Cre;Rb1lox/lox;p53lox/lox) and a triple-knockout model with additional Skp2 germline knockout ("TKO": Osx1-Cre;Rb1lox/lox;p53lox/lox;Skp2-/-), followed by qPCR and immunohistochemistry validation. To investigate the clinical implications of our results, we analyzed a human osteosarcoma patient cohort ("NCI-TARGET OS") with RNA-seq and clinical data. We found large differences in gene expression after SKP2 knockout. Surprisingly, we observed increased expression of genes related to immune microenvironment infiltration in TKO tumors, especially the signature genes for macrophages and to a lesser extent, T cells, B cells, and vascular cells. We also uncovered a set of relevant transcription factors that may mediate these changes. In osteosarcoma patient cohorts, high expression of genes upregulated in TKO was correlated with favorable overall survival, which was largely explained by the macrophage gene signatures. This relationship was further supported by our finding that SKP2 expression was negatively correlated with macrophage infiltration in the NCI-TARGET osteosarcoma and the TCGA Sarcoma cohorts. Overall, our findings indicate that SKP2 may mediate immune exclusion from the osteosarcoma tumor microenvironment, suggesting that SKP2 modulation in osteosarcoma may induce antitumor immune activation.

Functional Contribution and Clinical Implication of Cancer-Associated Fibroblasts in Glioblastoma.

PURPOSE: The abundance and biological contribution of cancer-associated fibroblasts (CAF) in glioblastoma (GBM) are poorly understood. Here, we aim to uncover its molecular signature, cellular roles, and potential tumorigenesis implications. EXPERIMENTAL DESIGN: We first applied single-cell RNA sequencing (RNA-seq) and bioinformatics analysis to identify and characterize stromal cells with CAF transcriptomic features in human GBM tumors. Then, we performed functional enrichment analysis and in vitro assays to investigate their interactions with malignant GBM cells. RESULTS: We found that CAF abundance was low but significantly correlated with tumor grade, poor clinical outcome, and activation of extracellular matrix remodeling using three large cohorts containing bulk RNA-seq data and clinical information. Proteomic analysis of a GBM-derived CAF line and its secretome revealed fibronectin (FN1) as a critical candidate factor mediating CAF functions. This was validated using in vitro cellular models, which demonstrated that CAF-conditioned media and recombinant FN1 could facilitate the migration and invasion of GBM cells. In addition, we showed that CAFs were more abundant in the mesenchymal-like state (or subtype) than in other states of GBMs. Interestingly, cell lines resembling the proneural state responded to the CAF signaling better for the migratory and invasive phenotypes. CONCLUSIONS: Overall, this study characterized the molecular features and functional impacts of CAFs in GBM, alluding to novel cell interactions mediated by CAFs in the GBM microenvironment.