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).

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.

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.

Six3 and Six6 jointly regulate the identities and developmental trajectories of multipotent retinal progenitor cells in the mouse retina.

Formation, maintenance, and differentiation of tissue-specific progenitor cells are fundamental tasks during organogenesis. Retinal development is an excellent model for dissecting these processes; mechanisms of retinal differentiation can be harnessed for retinal regeneration toward curing blindness. Using single-cell RNA sequencing of embryonic mouse eye cups in which transcription factor Six3 was conditionally inactivated in peripheral retinas on top of germline deletion of its close paralog Six6 ("DKO"), we identified cell clusters and then inferred developmental trajectories in the integrated dataset. In control retinas, naïve retinal progenitor cells had two major trajectories leading to ciliary margin cells and retinal neurons, respectively. The ciliary margin trajectory was directly from naïve retinal progenitor cells at G1 phase, and the retinal neuron trajectory was through a neurogenic state marked by Atoh7 expression. Upon Six3 and Six6 dual deficiency, both naïve and neurogenic retinal progenitor cells were defective. Ciliary margin differentiation was enhanced, and multi-lineage retinal differentiation was disrupted. An ectopic neuronal trajectory lacking the Atoh7+ state led to ectopic neurons. Differential expression analysis not only confirmed previous phenotype studies but also identified novel candidate genes regulated by Six3/Six6 . Six3 and Six6 were jointly required for balancing the opposing gradients of the Fgf and Wnt signaling in the central-peripheral patterning of the eye cups. Taken together, we identify transcriptomes and developmental trajectories jointly regulated by Six3 and Six6, providing deeper insight into molecular mechanisms underlying early retinal differentiation.

SKP2 knockout in Rb1/p53 deficient mouse models of osteosarcoma induces immune infiltration and drives a transcriptional program with a favorable prognosis.

Purpose: Osteosarcoma (OS) is an aggressive bone malignancy with a poor prognosis. One putative proto-oncogene in OS is SKP2, encoding a substrate recognition factor of the SCF E3 ubiquitin ligase. We previously demonstrated that SKP2 knockout in murine OS improved survival and delayed tumorigenesis. Here we aim to define the SKP2 drives transcriptional program and its clinical implication in OS. Experimental Design: We performed RNA-sequencing (RNA-seq) on tumors from a transgenic OS 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-/-). We validated our RNA-seq findings using qPCR and immunohistochemistry. To investigate the clinical implications of our results, we analyzed a human OS patient cohort ("NCI-TARGET OS") with RNA-seq and clinical data. Results: We found large differences in gene expression after SKP2 knockout. Strikingly, we observed increased expression of genes related to immune microenvironment infiltration in TKO tumors. We observed significant increases in 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 the changes. In OS 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 OS and the TCGA Sarcoma cohort. Conclusion: Our findings indicate that SKP2 may mediate immune exclusion from the OS tumor microenvironment, suggesting that SKP2 modulation in OS may induce anti-tumor immune activation.

Single-cell transcriptomics uncovers a non-autonomous Tbx1-dependent genetic program controlling cardiac neural crest cell development.

Disruption of cardiac neural crest cells (CNCCs) results in congenital heart disease, yet we do not understand the cell fate dynamics as these cells differentiate to vascular smooth muscle cells. Here we performed single-cell RNA-sequencing of NCCs from the pharyngeal apparatus with the heart in control mouse embryos and when Tbx1, the gene for 22q11.2 deletion syndrome, is inactivated. We uncover three dynamic transitions of pharyngeal NCCs expressing Tbx2 and Tbx3 through differentiated CNCCs expressing cardiac transcription factors with smooth muscle genes. These transitions are altered non-autonomously by loss of Tbx1. Further, inactivation of Tbx2 and Tbx3 in early CNCCs results in aortic arch branching defects due to failed smooth muscle differentiation. Loss of Tbx1 interrupts mesoderm to CNCC cell-cell communication with upregulation and premature activation of BMP signaling and reduced MAPK signaling, as well as alteration of other signaling, and failed dynamic transitions of CNCCs leading to disruption of aortic arch artery formation and cardiac outflow tract septation.

Systematic Comparison of Pancreatic Ductal Adenocarcinoma Models Identifies a Conserved Highly Plastic Basal Cell State.

Intratumoral heterogeneity and cellular plasticity have emerged as hallmarks of cancer, including pancreatic ductal adenocarcinoma (PDAC). As PDAC portends a dire prognosis, a better understanding of the mechanisms underpinning cellular diversity in PDAC is crucial. Here, we investigated the cellular heterogeneity of PDAC cancer cells across a range of in vitro and in vivo growth conditions using single-cell genomics. Heterogeneity contracted significantly in two-dimensional and three-dimensional cell culture models but was restored upon orthotopic transplantation. Orthotopic transplants reproducibly acquired cell states identified in autochthonous PDAC tumors, including a basal state exhibiting coexpression and coaccessibility of epithelial and mesenchymal genes. Lineage tracing combined with single-cell transcriptomics revealed that basal cells display high plasticity in situ. This work defines the impact of cellular growth conditions on phenotypic diversity and uncovers a highly plastic cell state with the capacity to facilitate state transitions and promote intratumoral heterogeneity in PDAC. SIGNIFICANCE: This work provides important insights into how different model systems of pancreatic ductal adenocarcinoma mold the phenotypic space of cancer cells, highlighting the power of in vivo models.

The interaction of SKP2 with p27 enhances the progression and stemness of osteosarcoma.

Osteosarcoma is a highly aggressive malignancy for which treatment has remained essentially unchanged for years. Our previous studies found that the F-box protein SKP2 is overexpressed in osteosarcoma, acting as a proto-oncogene; p27Kip1 (p27) is an inhibitor of cyclin-dependent kinases and a downstream substrate of SKP2-mediated ubiquitination. Overexpression of SKP2 and underexpression of p27 are common characteristics of cancer cells. The SCFSKP2 E3 ligase ubiquitinates Thr187-phosphorylated p27 for proteasome degradation, which can be abolished by a Thr187Ala knock-in (p27T187A KI) mutation. RB1 and TP53 are two major tumor suppressors commonly coinactivated in osteosarcoma. We generated a mouse model with a double knockout (DKO) of Rb1 and Trp53 within cells of the osteoblastic lineage, which developed osteosarcoma with full penetrance. When p27T187A KI mice were crossed on to the DKO background, p27T187A protein was found to accumulate in osteosarcoma tumor tissues. Furthermore, p27T187A promoted apoptosis in DKO tumors, slowed disease progression, and significantly prolonged overall survival. RNA sequencing analysis also linked the SCFSKP2 -p27T187A axis to potentially reduced cancer stemness. Given that RB1 and TP53 loss or coinactivation is common in human osteosarcoma, our study suggests that inhibiting the SKP2-p27 axis may represent a desirable therapeutic strategy for this cancer.