Pan-cancer single-cell RNA-seq identifies recurring programs of cellular heterogeneity

Transcriptomic intratumor heterogeneity of breast cancer patient-derived organoids may reflect the unique biological features of the tumor of origin

BACKGROUND

The intratumor heterogeneity (ITH) of cancer cells plays an important role in breast cancer resistance and recurrence. To develop better therapeutic strategies, it is necessary to understand the molecular mechanisms underlying ITH and their functional significance. Patient-derived organoids (PDOs) have recently been utilized in cancer research. They can also be used to study ITH as cancer cell diversity is thought to be maintained within the organoid line. However, no reports investigated intratumor transcriptomic heterogeneity in organoids derived from patients with breast cancer. This study aimed to investigate transcriptomic ITH in breast cancer PDOs.

METHODS

We established PDO lines from ten patients with breast cancer and performed single-cell transcriptomic analysis. First, we clustered cancer cells for each PDO using the Seurat package. Then, we defined and compared the cluster-specific gene signature (ClustGS) corresponding to each cell cluster in each PDO.

RESULTS

Cancer cells were clustered into 3-6 cell populations with distinct cellular states in each PDO line. We identified 38 clusters with ClustGS in 10 PDO lines and used Jaccard similarity index to compare the similarity of these signatures. We found that 29 signatures could be categorized into 7 shared meta-ClustGSs, such as those related to the cell cycle or epithelial-mesenchymal transition, and 9 signatures were unique to single PDO lines. These unique cell populations appeared to represent the characteristics of the original tumors derived from patients.

CONCLUSIONS

We confirmed the existence of transcriptomic ITH in breast cancer PDOs. Some cellular states were commonly observed in multiple PDOs, whereas others were specific to single PDO lines. The combination of these shared and unique cellular states formed the ITH of each PDO.

Prediction for recurrent non-muscle invasive bladder cancer

Non-muscle invasive bladder cancer (NMIBC) has a high recurrence rate, which places a significant burden on both patients and the healthcare system. Hence, it holds significant importance to predict the recurrence risk following treatment for individuals diagnosed with non-muscle invasive bladder cancer (NMIBC). As new generation technologies continue to emerge, an increasing number of recurrence risk prediction tools are being developed and discovered. This article provides an overview of the primary recurrence risk prediction tools currently available, including the liquid biopsy, tissue biopsy, and risk prediction tables. Each of these tools is described in detail and illustrated with relevant examples. Furthermore, we conduct an analysis of the advantages and disadvantages of these tools. This article aims to enhance the reader's understanding of the current progress in recurrence prediction tools and encourage their practical utilization in the fields of precision medicine and public health.

Deciphering functional tumor states at single-cell resolution

Within a tumor, cancer cells exist in different states that are associated with distinct tumor functions, including proliferation, differentiation, invasion, metastasis, and resistance to anti-cancer therapy. The identification of the gene regulatory networks underpinning each state is essential for better understanding functional tumor heterogeneity and revealing tumor vulnerabilities. Here, we review the different studies identifying tumor states by single-cell sequencing approaches and the mechanisms that promote and sustain these functional states and regulate their transitions. We also describe how different tumor states are spatially distributed and interact with the specific stromal cells that compose the tumor microenvironment. Finally, we discuss how the understanding of tumor plasticity and transition states can be used to develop new strategies to improve cancer therapy.

Pancancer Analysis of the Prognostic and Immunotherapeutic Value of Progestin and AdipoQ Receptor 4

AdipoQ receptor 4 (PAQR4) belongs to the family of progestin and AdipoQ receptors. PAQR4 plays an oncogenic role in lung and breast cancer. However, systematic pancancer analyses of PAQR4 have not been performed. The purpose was to investigate the prognostic and immunological significance of PAQR4 across 31 tumor types. Data were obtained from the following sources: TCGA, GEO, UALCAN, TIMER, GEPIA2, KM plotter, and TISIDB databases. The results proved that PAQR4 expression was significantly elevatory in most cancer types. We then explored the utility of PAQR4 as a prognostic indicator across all cancers. Using Cox proportional risk regression models, it has been demonstrated that PAQR4 is an independent risk factor in. High PAQR4 expression was not associated with other prognostic indicators, including overall survival, disease-free interval, disease-specific survival, and progression-free period. Subsequently, we explored the immunological value of PAQR4 and found that PAQR4 expression significantly correlated with tumor mutational burden, microsatellite instability, neoantigen, and immune checkpoint genes in tumors. It also significantly negatively correlated with most tumors' ESTIMATE scores, indicating that PAQR4 can influence the cellular composition of the tumor microenvironment. Our findings suggest the immunotherapeutic potential of PAQR4 in tumors. Finally, we explored the role of PAQR4 in tumor drug resistance and found that PAQR4 expression affected the sensitivity to multiple chemotherapeutic agents. A significant role for PAQR4 in tumor immunity is evident in these studies, as well as its potential role in cancer diagnosis, prognosis, and treatment precision.

Data analysis guidelines for single-cell RNA-seq in biomedical studies and clinical applications

The application of single-cell RNA sequencing (scRNA-seq) in biomedical research has advanced our understanding of the pathogenesis of disease and provided valuable insights into new diagnostic and therapeutic strategies. With the expansion of capacity for high-throughput scRNA-seq, including clinical samples, the analysis of these huge volumes of data has become a daunting prospect for researchers entering this field. Here, we review the workflow for typical scRNA-seq data analysis, covering raw data processing and quality control, basic data analysis applicable for almost all scRNA-seq data sets, and advanced data analysis that should be tailored to specific scientific questions. While summarizing the current methods for each analysis step, we also provide an online repository of software and wrapped-up scripts to support the implementation. Recommendations and caveats are pointed out for some specific analysis tasks and approaches. We hope this resource will be helpful to researchers engaging with scRNA-seq, in particular for emerging clinical applications.

High-plex imaging of RNA and proteins at subcellular resolution in fixed tissue by spatial molecular imaging

Resolving the spatial distribution of RNA and protein in tissues at subcellular resolution is a challenge in the field of spatial biology. We describe spatial molecular imaging, a system that measures RNAs and proteins in intact biological samples at subcellular resolution by performing multiple cycles of nucleic acid hybridization of fluorescent molecular barcodes. We demonstrate that spatial molecular imaging has high sensitivity (one or two copies per cell) and very low error rate (0.0092 false calls per cell) and background (~0.04 counts per cell). The imaging system generates three-dimensional, super-resolution localization of analytes at ~2 million cells per sample. Cell segmentation is morphology based using antibodies, compatible with formalin-fixed, paraffin-embedded samples. We measured multiomic data (980 RNAs and 108 proteins) at subcellular resolution in formalin-fixed, paraffin-embedded tissues (nonsmall cell lung and breast cancer) and identified >18 distinct cell types, ten unique tumor microenvironments and 100 pairwise ligand-receptor interactions. Data on >800,000 single cells and ~260 million transcripts can be accessed at http://nanostring.com/CosMx-dataset .

K27M in canonical and noncanonical H3 variants occurs in distinct oligodendroglial cell lineages in brain midline gliomas

Canonical (H3.1/H3.2) and noncanonical (H3.3) histone 3 K27M-mutant gliomas have unique spatiotemporal distributions, partner alterations and molecular profiles. The contribution of the cell of origin to these differences has been challenging to uncouple from the oncogenic reprogramming induced by the mutation. Here, we perform an integrated analysis of 116 tumors, including single-cell transcriptome and chromatin accessibility, 3D chromatin architecture and epigenomic profiles, and show that K27M-mutant gliomas faithfully maintain chromatin configuration at developmental genes consistent with anatomically distinct oligodendrocyte precursor cells (OPCs). H3.3K27M thalamic gliomas map to prosomere 2-derived lineages. In turn, H3.1K27M ACVR1-mutant pontine gliomas uniformly mirror early ventral NKX6-1⁺/SHH-dependent brainstem OPCs, whereas H3.3K27M gliomas frequently resemble dorsal PAX3⁺/BMP-dependent progenitors. Our data suggest a context-specific vulnerability in H3.1K27M-mutant SHH-dependent ventral OPCs, which rely on acquisition of ACVR1 mutations to drive aberrant BMP signaling required for oncogenesis. The unifying action of K27M mutations is to restrict H3K27me3 at PRC2 landing sites, whereas other epigenetic changes are mainly contingent on the cell of origin chromatin state and cycling rate.

Cervical cancer heterogeneity: a constant battle against viruses and drugs

Cervical cancer is the first identified human papillomavirus (HPV) associated cancer and the most promising malignancy to be eliminated. However, the ever-changing virus subtypes and acquired multiple drug resistance continue to induce failure of tumor prevention and treatment. The exploration of cervical cancer heterogeneity is the crucial way to achieve effective prevention and precise treatment. Tumor heterogeneity exists in various aspects including the immune clearance of viruses, tumorigenesis, neoplasm recurrence, metastasis and drug resistance. Tumor development and drug resistance are often driven by potential gene amplification and deletion, not only somatic genomic alterations, but also copy number amplifications, histone modification and DNA methylation. Genomic rearrangements may occur by selection effects from chemotherapy or radiotherapy which exhibits genetic intra-tumor heterogeneity in advanced cervical cancers. The combined application of cervical cancer therapeutic vaccine and immune checkpoint inhibitors has become an effective strategy to address the heterogeneity of treatment. In this review, we will integrate classic and recently updated epidemiological data on vaccination rates, screening rates, incidence and mortality of cervical cancer patients worldwide aiming to understand the current situation of disease prevention and control and identify the direction of urgent efforts. Additionally, we will focus on the tumor environment to summarize the conditions of immune clearance and gene integration after different HPV infections and to explore the genomic factors of tumor heterogeneity. Finally, we will make a thorough inquiry into completed and ongoing phase III clinical trials in cervical cancer and summarize molecular mechanisms of drug resistance among chemotherapy, radiotherapy, biotherapy, and immunotherapy.

Opportunities and challenges in interpretable deep learning for drug sensitivity prediction of cancer cells

In precision oncology, therapy stratification is done based on the patients' tumor molecular profile. Modeling and prediction of the drug response for a given tumor molecular type will further improve therapeutic decision-making for cancer patients. Indeed, deep learning methods hold great potential for drug sensitivity prediction, but a major problem is that these models are black box algorithms and do not clarify the mechanisms of action. This puts a limitation on their clinical implementation. To address this concern, many recent studies attempt to overcome these issues by developing interpretable deep learning methods that facilitate the understanding of the logic behind the drug response prediction. In this review, we discuss strengths and limitations of recent approaches, and suggest future directions that could guide further improvement of interpretable deep learning in drug sensitivity prediction in cancer research.

In vivo growth of subclones derived from Lewis lung carcinoma is determined by the tumor microenvironment

Heterogeneity is a fundamental feature of human tumors and plays a major role in drug resistance and disease progression. In the present study, we selected single-cell-derived cell lines (SCDCLs) derived from Lewis lung carcinoma (LLC1) cells to investigate tumorigenesis and heterogeneity. SCDCLs were generated using limiting dilution. Five SCDCLs were subcutaneously injected into wild-type C57BL/6N mice; however, they displayed significant differences in tumor growth. Subclone SCC1 grew the fastest in vivo, whereas it grew slower in vitro. The growth pattern of SCC2 was the opposite to that of SCC1. Genetic differences in these two subclones showed marked differences in cell adhesion and proliferation. Pathway enrichment results indicate that signal transduction and immune system responses were the most significantly altered functional categories in SCC2 cells compared to those in SCC1 cells in vitro. The number and activation of CD3⁺ and CD8⁺ T cells and NK cells in the tumor tissue of tumor-bearing mice inoculated with SCC2 were significantly higher, whereas those of myeloid cells were significantly lower, than those in the SCC1 and LLC1 groups. Our results suggest that the in vivo growth of two subclones derived from LLC1 was determined by the tumor microenvironment rather than their intrinsic proliferative cell characteristics.

scSemiGAN: a single-cell semi-supervised annotation and dimensionality reduction framework based on generative adversarial network

MOTIVATION

Cell-type annotation plays a crucial role in single-cell RNA-seq (scRNA-seq) data analysis. As more and more well-annotated scRNA-seq reference data are publicly available, automatical label transference algorithms are gaining popularity over manual marker gene-based annotation methods. However, most existing methods fail to unify cell-type annotation with dimensionality reduction and are unable to generate deep latent representation from the perspective of data generation.

RESULTS

In this article, we propose scSemiGAN, a single-cell semi-supervised cell-type annotation and dimensionality reduction framework based on a generative adversarial network, to overcome these challenges, modeling scRNA-seq data from the aspect of data generation. Our proposed scSemiGAN is capable of performing deep latent representation learning and cell-type label prediction simultaneously. Through extensive comparison with four state-of-the-art annotation methods on diverse simulated and real scRNA-seq datasets, scSemiGAN achieves competitive or superior performance in multiple downstream tasks including cell-type annotation, latent representation visualization, confounding factor removal and enrichment analysis.

AVAILABILITY AND IMPLEMENTATION

The code and data of scSemiGAN are available on GitHub: https://github.com/rafa-nadal/scSemiGAN.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Oral cancer patients experience mechanical and chemical sensitivity at the site of the cancer

INTRODUCTION

Oral cancer patients suffer severe chronic and mechanically-induced pain at the site of the cancer. Our clinical experience is that oral cancer patients report new sensitivity to spicy foods. We hypothesized that in cancer patients, mechanical and chemical sensitivity would be greater when measured at the cancer site compared to a contralateral matched normal site.

METHODS

We determined mechanical pain thresholds (MPT) on the right and left sides of the tongue of 11 healthy subjects, and at the cancer and contralateral matched normal site in 11 oral cancer patients in response to von Frey filaments in the range of 0.008 to 300 g (normally not reported as painful). We evaluated chemical sensitivity in 13 healthy subjects and seven cancer patients, who rated spiciness/pain on a visual analog scale in response to exposure to six paper strips impregnated with capsaicin (0-10 mM).

RESULTS

Mechanical detection thresholds (MDT) were recorded for healthy subjects, but not MPTs. By contrast, MPTs were measured at the site of the cancer in oral cancer patients (7/11 patients). No MPTs were measured at the cancer patients' contralateral matched normal sites. Measured MPTs were correlated with patients' responses to the University of California Oral Cancer Pain Questionnaire. Capsaicin sensitivity at the site of the cancer was evident in cancer patients by a leftward shift of the cancer site capsaicin dose-response curve compared to that of the patient's contralateral matched normal site. We detected no difference in capsaicin sensitivity on the right and left sides of tongues of healthy subjects.

CONCLUSIONS

Mechanical and chemical sensitivity testing was well tolerated by the majority of oral cancer patients. Sensitivity is greater at the site of the cancer than at a contralateral matched normal site.

The role of transcription factors in the acquisition of the four latest proposed hallmarks of cancer and corresponding enabling characteristics

With the recent description of the molecular and cellular characteristics that enable acquisition of both core and new hallmarks of cancer, the consequences of transcription factor dysregulation in the hallmarks scheme has become increasingly evident. Dysregulation or mutation of transcription factors has long been recognized in the development of cancer where alterations in these key regulatory molecules can result in aberrant gene expression and consequential blockade of normal cellular differentiation. Here, we provide an up-to-date review of involvement of dysregulated transcription factor networks with the most recently reported cancer hallmarks and enabling characteristic properties. We present some illustrative examples of the impact of dysregulated transcription factors, specifically focusing on the characteristics of phenotypic plasticity, non-mutational epigenetic reprogramming, polymorphic microbiomes, and senescence. We also discuss how new insights into transcription factor dysregulation in cancer is contributing to addressing current therapeutic challenges.

Microfluidics for long-term single-cell time-lapse microscopy: Advances and applications

Cells are inherently dynamic, whether they are responding to environmental conditions or simply at equilibrium, with biomolecules constantly being made and destroyed. Due to their small volumes, the chemical reactions inside cells are stochastic, such that genetically identical cells display heterogeneous behaviors and gene expression profiles. Studying these dynamic processes is challenging, but the development of microfluidic methods enabling the tracking of individual prokaryotic cells with microscopy over long time periods under controlled growth conditions has led to many discoveries. This review focuses on the recent developments of one such microfluidic device nicknamed the mother machine. We overview the original device design, experimental setup, and challenges associated with this platform. We then describe recent methods for analyzing experiments using automated image segmentation and tracking. We further discuss modifications to the experimental setup that allow for time-varying environmental control, replicating batch culture conditions, cell screening based on their dynamic behaviors, and to accommodate a variety of microbial species. Finally, this review highlights the discoveries enabled by this technology in diverse fields, such as cell-size control, genetic mutations, cellular aging, and synthetic biology.

Transcriptional Heterogeneity of Cellular Senescence in Cancer

Cellular senescence plays a paradoxical role in tumorigenesis through the expression of diverse senescence-associated (SA) secretory phenotypes (SASPs). The heterogeneity of SA gene expression in cancer cells not only promotes cancer stemness but also protects these cells from chemotherapy. Despite the potential correlation between cancer and SA biomarkers, many transcriptional changes across distinct cell populations remain largely unknown. During the past decade, single-cell RNA sequencing (scRNA-seq) technologies have emerged as powerful experimental and analytical tools to dissect such diverse senescence-derived transcriptional changes. Here, we review the recent sequencing efforts that successfully characterized scRNA-seq data obtained from diverse cancer cells and elucidated the role of senescent cells in tumor malignancy. We further highlight the functional implications of SA genes expressed specifically in cancer and stromal cell populations in the tumor microenvironment. Translational research leveraging scRNA-seq profiling of SA genes will facilitate the identification of novel expression patterns underlying cancer susceptibility, providing new therapeutic opportunities in the era of precision medicine.

Single-cell sequencing and establishment of an 8-gene prognostic model for pancreatic cancer patients

Background

Single-cell sequencing (SCS) technologies enable analysis of gene structure and expression data at single-cell resolution. However, SCS analysis in pancreatic cancer remains largely unexplored.

Methods

We downloaded pancreatic cancer SCS data from different databases and applied appropriate dimensionality reduction algorithms. We identified 10 cell types and subsequently screened differentially expressed marker genes of these 10 cell types using FindAllMarkers analysis. Also, we evaluated the tumor immune microenvironment based on ESTIMATE and MCP-counter. Statistical enrichment was evaluated using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis. We used all candidate gene sets in KEGG database to perform gene set enrichment analysis. We used LASSO regression to reduce the number of genes in the pancreatic risk model by R package glmnet, followed by rtPCR to validate the expression of the signature genes in different pancreatic cancer cell lines.

Results

We identified 15 cell subpopulations by dimension reduction and data clustering. We divided the 15 subpopulations into 10 distinct cell types based on marker gene expression. Then, we performed functional enrichment analysis for the 352 marker genes in pancreatic cancer cells. Based on RNA expression data and prognostic information from TCGA and GEO datasets, we identified 42 prognosis-related genes, including 5 protective genes and 37 high-risk genes, which we used to identified two molecular subtypes. C1 subtype was associated with a better prognosis, whereas C2 subtype was associated with a worse prognosis. Moreover, chemokine and chemokine receptor genes were differentially expressed between C1 and C2 subtypes. Functional and pathway enrichment uncovered functional differences between C1 and C2 subtype. We identified eight genes that could serve as potential biomarkers for prognosis prediction in pancreatic cancer patients. These genes were used to establish an 8-gene pancreatic cancer prognostic model.

Conclusions

We established an 8-gene pancreatic cancer prognostic model. This model can meaningfully predict prognosis and treatment response in pancreatic cancer patients.

Gene expression based inference of cancer drug sensitivity

Inter and intra-tumoral heterogeneity are major stumbling blocks in the treatment of cancer and are responsible for imparting differential drug responses in cancer patients. Recently, the availability of high-throughput screening datasets has paved the way for machine learning based personalized therapy recommendations using the molecular profiles of cancer specimens. In this study, we introduce Precily, a predictive modeling approach to infer treatment response in cancers using gene expression data. In this context, we demonstrate the benefits of considering pathway activity estimates in tandem with drug descriptors as features. We apply Precily on single-cell and bulk RNA sequencing data associated with hundreds of cancer cell lines. We then assess the predictability of treatment outcomes using our in-house prostate cancer cell line and xenografts datasets exposed to differential treatment conditions. Further, we demonstrate the applicability of our approach on patient drug response data from The Cancer Genome Atlas and an independent clinical study describing the treatment journey of three melanoma patients. Our findings highlight the importance of chemo-transcriptomics approaches in cancer treatment selection.

Predicting treatment response in cancer remains a highly complex task. Here, the authors develop Precily, a deep neural network framework to predict treatment response in cancer by considering gene expression, pathway activity estimates and drug features, and test this method in multiple datasets and preclinical models.

Single-cell and WGCNA uncover a prognostic model and potential oncogenes in colorectal cancer

Background

Colorectal cancer (CRC) is one of the leading causes of cancer-related death worldwide. Single-cell transcriptome sequencing (scRNA-seq) can provide accurate gene expression data for individual cells. In this study, a new prognostic model was constructed by scRNA-seq and bulk transcriptome sequencing (bulk RNA-seq) data of CRC samples to develop a new understanding of CRC.

Methods

CRC scRNA-seq data were downloaded from the GSE161277 database, and CRC bulk RNA-seq data were downloaded from the TCGA and GSE17537 databases. The cells were clustered by the FindNeighbors and FindClusters functions in scRNA-seq data. CIBERSORTx was applied to detect the abundance of cell clusters in the bulk RNA-seq expression matrix. WGCNA was performed with the expression profiles to construct the gene coexpression networks of TCGA-CRC. Next, we used a tenfold cross test to construct the model and a nomogram to assess the independence of the model for clinical application. Finally, we examined the expression of the unreported model genes by qPCR and immunohistochemistry. A clone formation assay and orthotopic colorectal tumour model were applied to detect the regulatory roles of unreported model genes.

Results

A total of 43,851 cells were included after quality control, and 20 cell clusters were classified by the FindCluster () function. We found that the abundances of C1, C2, C4, C5, C15, C16 and C19 were high and the abundances of C7, C10, C11, C13, C14 and C17 were low in CRC tumour tissues. Meanwhile, the results of survival analysis showed that high abundances of C4, C11 and C13 and low abundances of C5 and C14 were associated with better survival. The WGCNA results showed that the red module was most related to the tumour and the C14 cluster, which contains 615 genes. Lasso Cox regression analysis revealed 8 genes (PBXIP1, MPMZ, SCARA3, INA, ILK, MPP2, L1CAM and FLNA), which were chosen to construct a risk model. In the model, the risk score features had the greatest impact on survival prediction, indicating that the 8-gene risk model can better predict prognosis. qPCR and immunohistochemistry analysis showed that the expression levels of MPZ, SCARA3, MPP2 and PBXIP1 were high in CRC tissues. The functional experiment results indicated that MPZ, SCARA3, MPP2 and PBXIP1 could promote the colony formation ability of CRC cells in vitro and tumorigenicity in vivo.

Conclusions

We constructed a risk model to predict the prognosis of CRC patients based on scRNA-seq and bulk RNA-seq data, which could be used for clinical application. We also identified 4 previously unreported model genes (MPZ, SCARA3, MPP2 and PBXIP1) as novel oncogenes in CRC. These results suggest that this model could potentially be used to evaluate the prognostic risk and provide potential therapeutic targets for CRC patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12575-022-00175-x.

Genetic and microenvironmental intra-tumor heterogeneity impacts colorectal cancer evolution and metastatic development

Colorectal cancer (CRC) is a highly diverse disease, where different genomic instability pathways shape genetic clonal diversity and tumor microenvironment. Although intra-tumor heterogeneity has been characterized in primary tumors, its origin and consequences in CRC outcome is not fully understood. Therefore, we assessed intra- and inter-tumor heterogeneity of a prospective cohort of 136 CRC samples. We demonstrate that CRC diversity is forged by asynchronous forms of molecular alterations, where mutational and chromosomal instability collectively boost CRC genetic and microenvironment intra-tumor heterogeneity. We were able to depict predictor signatures of cancer-related genes that can foresee heterogeneity levels across the different tumor consensus molecular subtypes (CMS) and primary tumor location. Finally, we show that high genetic and microenvironment heterogeneity are associated with lower metastatic potential, whereas late-emerging copy number variations favor metastasis development and polyclonal seeding. This study provides an exhaustive portrait of the interplay between genetic and microenvironment intra-tumor heterogeneity across CMS subtypes, depicting molecular events with predictive value of CRC progression and metastasis development.

A machine learning-based method for automatically identifying novel cells in annotating single-cell RNA-seq data

MOTIVATION

Single-cell RNA sequencing (scRNA-seq) has been widely used to decompose complex tissues into functionally distinct cell types. The first and usually the most important step of scRNA-seq data analysis is to accurately annotate the cell labels. In recent years, many supervised annotation methods have been developed and shown to be more convenient and accurate than unsupervised cell clustering. One challenge faced by all the supervised annotation methods is the identification of the novel cell type, which is defined as the cell type that is not present in the training data, only exists in the testing data. Existing methods usually label the cells simply based on the correlation coefficients or confidence scores, which sometimes results in an excessive number of unlabeled cells.

RESULTS

We developed a straightforward yet effective method combining autoencoder with iterative feature selection to automatically identify novel cells from scRNA-seq data. Our method trains an autoencoder with the labeled training data and applies the autoencoder to the testing data to obtain reconstruction errors. By iteratively selecting features that demonstrate a bi-modal pattern and reclustering the cells using the selected feature, our method can accurately identify novel cells that are not present in the training data. We further combined this approach with a support vector machine to provide a complete solution for annotating the full range of cell types. Extensive numerical experiments using five real scRNA-seq datasets demonstrated favorable performance of the proposed method over existing methods serving similar purposes.

AVAILABILITY AND IMPLEMENTATION

Our R software package CAMLU is publicly available through the Zenodo repository (https://doi.org/10.5281/zenodo.7054422) or GitHub repository (https://github.com/ziyili20/CAMLU).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Single cell RNA-sequencing: A powerful yet still challenging technology to study cellular heterogeneity

Almost all biomedical research to date has relied upon mean measurements from cell populations, however it is well established that what it is observed at this macroscopic level can be the result of many interactions of several different single cells. Thus, the observable macroscopic 'average' cannot outright be used as representative of the 'average cell'. Rather, it is the resulting emerging behaviour of the actions and interactions of many different cells. Single-cell RNA sequencing (scRNA-Seq) enables the comparison of the transcriptomes of individual cells. This provides high-resolution maps of the dynamic cellular programmes allowing us to answer fundamental biological questions on their function and evolution. It also allows to address medical questions such as the role of rare cell populations contributing to disease progression and therapeutic resistance. Furthermore, it provides an understanding of context-specific dependencies, namely the behaviour and function that a cell has in a specific context, which can be crucial to understand some complex diseases, such as diabetes, cardiovascular disease and cancer. Here, we provide an overview of scRNA-Seq, including a comparative review of emerging technologies and computational pipelines. We discuss the current and emerging applications and focus on tumour heterogeneity a clear example of how scRNA-Seq can provide new understanding of a complex disease. Additionally, we review the limitations and highlight the need of powerful computational pipelines and reproducible protocols for the broader acceptance of this technique in basic and clinical research.

Inferring gene regulation from stochastic transcriptional variation across single cells at steady state

Regulatory relationships between transcription factors (TFs) and their target genes lie at the heart of cellular identity and function; however, uncovering these relationships is often labor-intensive and requires perturbations. Here, we propose a principled framework to systematically infer gene regulation for all TFs simultaneously in cells at steady state by leveraging the intrinsic variation in the transcriptional abundance across single cells. Through modeling and simulations, we characterize how transcriptional bursts of a TF gene are propagated to its target genes, including the expected ranges of time delay and magnitude of maximum covariation. We distinguish these temporal trends from the time-invariant covariation arising from cell states, and we delineate the experimental and technical requirements for leveraging these small but meaningful cofluctuations in the presence of measurement noise. While current technology does not yet allow adequate power for definitively detecting regulatory relationships for all TFs simultaneously in cells at steady state, we investigate a small-scale dataset to inform future experimental design. This study supports the potential value of mapping regulatory connections through stochastic variation, and it motivates further technological development to achieve its full potential.

Glycosphingolipids are mediators of cancer plasticity through independent signaling pathways

The molecular repertoire promoting cancer cell plasticity is not fully elucidated. Here, we propose that glycosphingolipids (GSLs), specifically the globo and ganglio series, correlate and promote the transition between epithelial and mesenchymal cells. The epithelial character of ovarian cancer remains stable throughout disease progression, and spatial glycosphingolipidomics reveals elevated globosides in the tumor compartment compared with the ganglioside-rich stroma. CRISPR-Cas9 knockin mediated truncation of endogenous E-cadherin induces epithelial-to-mesenchymal transition (EMT) and decreases globosides. The transcriptomics analysis identifies the ganglioside-synthesizing enzyme ST8SIA1 to be consistently elevated in mesenchymal-like samples, predicting poor outcome. Subsequent deletion of ST8SIA1 induces epithelial cell features through mTOR^S2448 phosphorylation, whereas loss of globosides in ΔA4GALT cells, resulting in EMT, is accompanied by increased ERK^Y202/T204 and AKT^S124. The GSL composition dynamics corroborate cancer cell plasticity, and further evidence suggests that mesenchymal cells are maintained through ganglioside-dependent, calcium-mediated mechanisms.

Pivotal models and biomarkers related to the prognosis of breast cancer based on the immune cell interaction network

The effect of breast cancer heterogeneity on prognosis of patients is still unclear, especially the role of immune cells in prognosis of breast cancer. In this study, single cell transcriptome sequencing data of breast cancer were used to analyze the relationship between breast cancer heterogeneity and prognosis. In this study, 14 cell clusters were identified in two single-cell datasets (GSE75688 and G118389). Proportion analysis of immune cells showed that NK cells were significantly aggregated in triple negative breast cancer, and the proportion of macrophages was significantly increased in primary breast cancer, while B cells, T cells, and neutrophils may be involved in the metastasis of breast cancer. The results of ligand receptor interaction network revealed that macrophages and DC cells were the most frequently interacting cells with other cells in breast cancer. The results of WGCNA analysis suggested that the MEblue module is most relevant to the overall survival time of triple negative breast cancer. Twenty-four prognostic genes in the blue module were identified by univariate Cox regression analysis and KM survival analysis. Multivariate regression analysis combined with risk analysis was used to analyze 24 prognostic genes to construct a prognostic model. The verification result of our prognostic model showed that there were significant differences in the expression of PCDH12, SLIT3, ACVRL1, and DLL4 genes between the high-risk group and the low-risk group, which can be used as prognostic biomarkers.

Pan-Cancer Analysis of the Immunological Role of PDIA5: A Potential Target for Immunotherapy

The aberrant protein disulfide isomerase A5 (PDIA5) expression was relevant to the poor prognosis of patients with human cancers. However, its relationship with the epigenetic and genetic alterations and its effect on tumor immunity is still lacking. In the present study, we comprehensively analyzed the immune infiltration role of PDIA5 in human cancers based on large-scale bioinformatics analyses and in vitro experiments. Obvious DNA methylation and moderate alteration frequency of PDIA5 were observed in human cancers. The expression level of PDIA5 was significantly correlated with infiltrated immune cells, immune pathways, and other immune signatures. We found that cancer cells and macrophages exhibited high PDIA5 expression in human cancers using the single-cell RNA sequencing analysis. We also demonstrated the interaction between PDIA5 and immune cells in glioblastoma multiforme (GBM). Multiplex immunofluorescence staining showed the upregulated expression level of PDIA5 and the increased number of M2 macrophage markers-CD163 positive cells in pan-cancer samples. Notably, PDIA5 silencing resulted in upregulated expression of PD-L1 and SPP1 in U251 cells. Silencing of PDIA5 in hepG2 cells, U251 cells, and PC3 cells contributed to a decline in their ability of proliferation, clone formation, and invasion and inhibited the migration of cocultured M2 macrophages. Additionally, PDIA5 also displayed predictive value in the immunotherapy response of both murine and human cancer cohorts. Overall, our findings indicated that PDIA5 might be a potential target for immunotherapies in cancers.

Cancer cell states recur across tumor types and form specific interactions with the tumor microenvironment

Transcriptional heterogeneity among malignant cells of a tumor has been studied in individual cancer types and shown to be organized into cancer cell states; however, it remains unclear to what extent these states span tumor types, constituting general features of cancer. Here, we perform a pan-cancer single-cell RNA-sequencing analysis across 15 cancer types and identify a catalog of gene modules whose expression defines recurrent cancer cell states including 'stress', 'interferon response', 'epithelial-mesenchymal transition', 'metal response', 'basal' and 'ciliated'. Spatial transcriptomic analysis linked the interferon response in cancer cells to T cells and macrophages in the tumor microenvironment. Using mouse models, we further found that induction of the interferon response module varies by tumor location and is diminished upon elimination of lymphocytes. Our work provides a framework for studying how cancer cell states interact with the tumor microenvironment to form organized systems capable of immune evasion, drug resistance and metastasis.

scPCOR-seq enables co-profiling of chromatin occupancy and RNAs in single cells

Cell-to-cell variation in gene expression is a widespread phenomenon, which may play important roles in cellular differentiation, function, and disease development^1–9. Chromatin is implicated in contributing to the cellular heterogeneity in gene expression^10–16. Fully understanding the mechanisms of cellular heterogeneity requires simultaneous measurement of RNA and occupancy of histone modifications and transcription factors on chromatin due to their critical roles in transcriptional regulation^17,18. We generally term the occupancy of histone modifications and transcription factors as Chromatin occupancy. Here, we report a technique, termed scPCOR-seq (single-cell Profiling of Chromatin Occupancy and RNAs Sequencing), for simultaneously profiling genome-wide chromatin protein binding or histone modification marks and RNA expression in the same cell. We demonstrated that scPCOR-seq can profile either H3K4me3 or RNAPII and RNAs in a mixture of human H1, GM12878 and 293 T cells at a single-cell resolution and either H3K4me3, RNAPII, or RNA profile can correctly separate the cells. Application of scPCOR-seq to the in vitro differentiation of the erythrocyte precursor CD36 cells from human CD34 stem or progenitor cells revealed that H3K4me3 and RNA exhibit distinct properties in clustering cells during differentiation. Overall, our work provides a promising approach to understand the relationships among different omics layers.

scPCOR-seq is a single-cell sequencing technique that enables simultaneous profiling of genome-wide chromatin protein binding or histone modification marks and RNA expression in the same cell.

Epithelial-Mesenchymal Plasticity in Tumor Immune Evasion

Epithelial-mesenchymal transition (EMT) is a fundamental process that occurs during embryogenesis and tissue repair. However, EMT can be hijacked by malignant cells, where it may promote immune evasion and metastasis. Classically considered a dichotomous transition, EMT in cancer has recently been considered a plastic process whereby malignant cells display and interconvert among hybrid epithelial/mesenchymal (E/M) states. Epithelial-mesenchymal plasticity (EMP) and associated hybrid E/M states are divergent from classical EMT, with unique immunomodulatory effects. Here, we review recent insights into the EMP-immune cross-talk, highlighting possible mechanisms of immune evasion conferred by hybrid E/M states and roles of immune cells in EMP.

Computational estimation of quality and clinical relevance of cancer cell lines

Immortal cancer cell lines (CCLs) are the most widely used system for investigating cancer biology and for the preclinical development of oncology therapies. Pharmacogenomic and genome‐wide editing screenings have facilitated the discovery of clinically relevant gene–drug interactions and novel therapeutic targets via large panels of extensively characterised CCLs. However, tailoring pharmacological strategies in a precision medicine context requires bridging the existing gaps between tumours and in vitro models. Indeed, intrinsic limitations of CCLs such as misidentification, the absence of tumour microenvironment and genetic drift have highlighted the need to identify the most faithful CCLs for each primary tumour while addressing their heterogeneity, with the development of new models where necessary. Here, we discuss the most significant limitations of CCLs in representing patient features, and we review computational methods aiming at systematically evaluating the suitability of CCLs as tumour proxies and identifying the best patient representative in vitro models. Additionally, we provide an overview of the applications of these methods to more complex models and discuss future machine‐learning‐based directions that could resolve some of the arising discrepancies.

Single Cell RNA Sequencing Reveals the Pathogenesis of Aortic Dissection Caused by Hypertension and Marfan Syndrome

Aortic dissection (AD) is mainly caused by hypertension and Marfan syndrome. However, it is unclear whether the cellular components and functions are different between the two causes. A total of 11 aortic samples were collected for single-cell RNA analysis and 20 clusters were disclosed, including VSMCs, fibroblasts, endothelial cells, T cells, B cells, monocytes, macrophages, mast cells, and neutrophils components. There were differences in cell subclusters and function between hypertension and Marfan patients. The cells also had different differentiations. Cellchat identified cell ligand–receptor interactions that were associated with hypertension and Marfan-induced AD involving SMC, fibroblast, mo-macrophages, and T-cell subsets. This study revealed the heterogeneity of cellular components and gene changes in hypertension and Marfan-induced AD. Through functional analysis and the changes in intercellular communication, the possible mechanisms of different causes of AD were explained from a new perspective, so we can better understand the occurrence and development of diseases.

From single-omics to interactomics: How can ligand-induced perturbations modulate single-cell phenotypes?

Cells suffer from perturbations by different stimuli, which, consequently, rise to individual alterations in their profile and function that may end up affecting the tissue as a whole. This is no different if we consider the effect of a therapeutic agent on a biological system. As cells are exposed to external ligands their profile can change at different single-omics levels. Detecting how these changes take place through different sequencing technologies is key to a better understanding of the effects of therapeutic agents. Single-cell RNA-sequencing stands out as one of the most common approaches for cell profiling and perturbation analysis. As a result, single-cell transcriptomics data can be integrated with other omics data sources, such as proteomics and epigenomics data, to clarify the perturbation effects and mechanism at the cell level. Appropriate computational tools are key to process and integrate the available information. This chapter focuses on the recent advances on ligand-induced perturbation and single-cell omics computational tools and algorithms, their current limitations, and how the deluge of data can be used to improve the current process of drug research and development.

Droplet Microfluidics for High-Resolution Virology

Microfluidics has enabled a new era of cellular and molecular assays due to the small length scales, parallelization, and the modularity of various analysis and actuation functions. Droplet microfluidics, in particular, has been instrumental in providing new tools for biology with its ability to quickly and reproducibly generate drops that act as individual reactors. A notable beneficiary of this technology has been single-cell RNA sequencing, which has revealed new heterogeneities and interactions for the fundamental unit of life. However, viruses far surpass the diversity of cellular life, affect the dynamics of all ecosystems, and are a chronic source of global health crises. Despite their impact on the world, high-throughput and high-resolution viral profiling has been difficult, with conventional methods being limited to population-level averaging, large sample volumes, and few cultivable hosts. Consequently, most viruses have not been identified and studied. Droplet microfluidics holds the potential to address many of these limitations and offers new levels of sensitivity and throughput for virology. This Feature highlights recent efforts that have applied droplet microfluidics to the detection and study of viruses, including for diagnostics, virus-host interactions, and cell-independent virus assays. In combination with traditional virology methods, droplet microfluidics should prove a potent tool toward achieving a better understanding of the most abundant biological species on Earth.

Melanopsin (Opn4) is an oncogene in cutaneous melanoma

The search for new therapeutical targets for cutaneous melanoma and other cancers is an ongoing task. We expanded this knowledge by evaluating whether opsins, light- and thermo-sensing proteins, could display tumor-modulatory effects on melanoma cancer. Using different experimental approaches, we show that melanoma cell proliferation is slower in the absence of Opn4, compared to Opn4WT due to an impaired cell cycle progression and reduced melanocyte inducing transcription factor (Mitf) expression. In vivo tumor progression of Opn4KO cells is remarkably reduced due to slower proliferation, and higher immune system response in Opn4KO tumors. Using pharmacological assays, we demonstrate that guanylyl cyclase activity is impaired in Opn4KO cells. Evaluation of Tumor Cancer Genome Atlas (TCGA) database confirms our experimental data as reduced MITF and OPN4 expression in human melanoma correlates with slower cell cycle progression and presence of immune cells in the tumor microenvironment (TME). Proteomic analyses of tumor bulk show that the reduced growth of Opn4KO tumors is associated with reduced Mitf signaling, higher translation of G2/M proteins, and impaired guanylyl cyclase activity. Conversely, in Opn4WT tumors increased small GTPase and an immune-suppressive TME are found. Such evidence points to OPN4 as an oncogene in melanoma, which could be pharmacologically targeted.

Integrated analysis of single-cell RNA-seq dataset and bulk RNA-seq dataset constructs a prognostic model for predicting survival in human glioblastoma

BACKGROUND

Glioblastoma (GBM) is the most common primary malignant brain tumor in adults. For patients with GBM, the median overall survival (OS) is 14.6 months and the 5-year survival rate is 7.2%. It is imperative to develop a reliable model to predict the survival probability in new GBM patients. To date, most prognostic models for predicting survival in GBM were constructed based on bulk RNA-seq dataset, which failed to accurately reflect the difference between tumor cores and peripheral regions, and thus show low predictive capability. An effective prognostic model is desperately needed in clinical practice.

METHODS

We studied single-cell RNA-seq dataset and The Cancer Genome Atlas-glioblastoma multiforme (TCGA-GBM) dataset to identify differentially expressed genes (DEGs) that impact the OS of GBM patients. We then applied the least absolute shrinkage and selection operator (LASSO) Cox penalized regression analysis to determine the optimal genes to be included in our risk score prognostic model. Then, we used another dataset to test the accuracy of our risk score prognostic model.

RESULTS

We identified 2128 DEGs from the single-cell RNA-seq dataset and 6461 DEGs from the bulk RNA-seq dataset. In addition, 896 DEGs associated with the OS of GBM patients were obtained. Five of these genes (LITAF, MTHFD2, NRXN3, OSMR, and RUFY2) were selected to generate a risk score prognostic model. Using training and validation datasets, we found that patients in the low-risk group showed better OS than those in the high-risk group. We validated our risk score model with the training and validating datasets and demonstrated that it can effectively predict the OS of GBM patients.

CONCLUSION

We constructed a novel prognostic model to predict survival in GBM patients by integrating a scRNA-seq dataset and a bulk RNA-seq dataset. Our findings may advance the development of new therapeutic targets and improve clinical outcomes for GBM patients.

Stepwise-edited, human melanoma models reveal mutations' effect on tumor and microenvironment

Establishing causal relationships between genetic alterations of human cancers and specific phenotypes of malignancy remains a challenge. We sequentially introduced mutations into healthy human melanocytes in up to five genes spanning six commonly disrupted melanoma pathways, forming nine genetically distinct cellular models of melanoma. We connected mutant melanocyte genotypes to malignant cell expression programs in vitro and in vivo, replicative immortality, malignancy, rapid tumor growth, pigmentation, metastasis, and histopathology. Mutations in malignant cells also affected tumor microenvironment composition and cell states. Our melanoma models shared genotype-associated expression programs with patient melanomas, and a deep learning model showed that these models partially recapitulated genotype-associated histopathological features as well. Thus, a progressive series of genome-edited human cancer models can causally connect genotypes carrying multiple mutations to phenotype.

Single‐cell RNA sequencing reveals the multi‐cellular ecosystem in different radiological components of pulmonary part‐solid nodules

Background

Early‐stage lung adenocarcinoma that radiologically manifests as part‐solid nodules, consisting of both ground‐glass and solid components, has distinctive growth patterns and prognosis. The characteristics of the tumour microenvironment and transcriptional features of the malignant cells of different radiological phenotypes remain poorly understood.

Methods

Twelve treatment‐naive patients with radiological part‐solid nodules were enrolled. After frozen pathology was confirmed as lung adenocarcinoma, two regions (ground‐glass and solid) from each of the 12 part‐solid nodules and 5 normal lung tissues from 5 of the12 patients were subjected to single‐cell sequencing by 10x Genomics. We used Seurat v3.1.5 for data integration and analysis.

Results

We comprehensively dissected the multicellular ecosystem of the ground‐glass and solid components of part‐solid nodules at the single‐cell resolution. In tumours, these components had comparable proportions of malignant cells. However, the angiogenesis, epithelial‐to‐mesenchymal transition, KRAS, p53, and cell‐cycle signalling pathways were significantly up‐regulated in malignant cells within solid components compared to those within ground‐glass components. For the tumour microenvironment, the relative abundance of myeloid and NK cells tended to be higher in solid components than in ground‐glass components. Slight subtype composition differences existed between the ground‐glass and solid components. The T/NK cell subsets’ cytotoxic function and the macrophages’ pro‐inflammation function were suppressed in solid components. Moreover, pericytes in solid components had a stronger communication related to angiogenesis promotion with endothelial cells and tumour cells.

Conclusion

The cellular landscape of ground‐glass components is significantly different from that of normal tissue and similar to that of solid components. However, transcriptional differences exist in the vital signalling pathways of malignant and immune cells within these components.

Computational Analyses of YY1 and Its Target RKIP Reveal Their Diagnostic and Prognostic Roles in Lung Cancer

Simple Summary

Lung cancer (LC) is the tumor with the highest global mortality rate. Novel personalized therapies are currently being tested (e.g., targeted inhibitors, the immune-checkpoint inhibitors), but they cannot yet prevent the very frequent relapse and generalized metastases observed in a large population of LC patients. Currently, there is an urgent need for novel reliable biomarkers for the prognosis and diagnosis of LC. Through the systematic analysis of multiple deposited expression datasets, this report aims to explore the role of the Yin-Yang 1 (YY1) transcription factor and its target the Raf Kinase Inhibitory Protein (RKIP) in LC. The computational analysis suggested the predictive diagnostic and prognostic roles for both YY1 and RKIP stimulating further studies for proving their implication as novel biomarkers, as well as therapeutically druggable targets in LC.

Abstract

Lung cancer (LC) represents a global threat, being the tumor with the highest mortality rate. Despite the introduction of novel therapies (e.g., targeted inhibitors, immune-checkpoint inhibitors), relapses are still very frequent. Accordingly, there is an urgent need for reliable predictive biomarkers and therapeutically druggable targets. Yin-Yang 1 (YY1) is a transcription factor that may work either as an oncogene or a tumor suppressor, depending on the genotype and the phenotype of the tumor. The Raf Kinase Inhibitory Protein (RKIP), is a tumor suppressor and immune enhancer often found downregulated in the majority of the examined cancers. In the present report, the role of both YY1 and RKIP in LC is thoroughly explored through the analysis of several deposited RNA and protein expression datasets. The computational analyses revealed that YY1 negatively regulates RKIP expression in LC, as corroborated by the deposited YY1-ChIP-Seq experiments and validated by their robust negative correlation. Additionally, YY1 expression is significantly higher in LC samples compared to normal matching ones, whereas RKIP expression is lower in LC and high in normal matching tissues. These observed differences, unlike many current biomarkers, bear a diagnostic significance, as proven by the ROC analyses. Finally, the survival data support the notion that both YY1 and RKIP might represent strong prognostic biomarkers. Overall, the reported findings indicate that YY1 and RKIP expression levels may play a role in LC as potential biomarkers and therapeutic targets. However, further studies will be necessary to validate the in silico results.

Single-cell transcriptomics identifies Gadd45b as a regulator of herpesvirus-reactivating neurons

Single-cell RNA sequencing (scRNA-seq) is a powerful technique for dissecting the complexity of normal and diseased tissues, enabling characterization of cell diversity and heterogeneous phenotypic states in unprecedented detail. However, this technology has been underutilized for exploring the interactions between the host cell and viral pathogens in latently infected cells. Herein, we use scRNA-seq and single-molecule sensitivity fluorescent in situ hybridization (smFISH) technologies to investigate host single-cell transcriptome changes upon the reactivation of a human neurotropic virus, herpes simplex virus-1 (HSV-1). We identify the stress sensor growth arrest and DNA damage-inducible 45 beta (Gadd45b) as a critical antiviral host factor that regulates HSV-1 reactivation events in a subpopulation of latently infected primary neurons. We show that distinct subcellular localization of Gadd45b correlates with the viral late gene expression program, as well as the expression of the viral transcription factor, ICP4. We propose that a hallmark of a "successful" or "aborted" HSV-1 reactivation state in primary neurons is determined by a unique subcellular localization signature of the stress sensor Gadd45b.

Lung Cancer Computational Biology and Resources

Comprehensive clinical, pathological, and molecular data, when appropriately integrated with advanced computational approaches, are transforming the way we characterize and study lung cancer. Clinically, cancer registry and publicly available historical clinical trial data enable retrospective analyses to examine how socioeconomic factors, patient demographics, and cancer characteristics affect treatment and outcome. Pathologically, digital pathology and artificial intelligence are revolutionizing histopathological image analyses, not only with improved efficiency and accuracy, but also by extracting additional information for prognostication and tumor microenvironment characterization. Genetically and molecularly, individual patient tumors and preclinical models of lung cancer are profiled by various high-throughput platforms to characterize the molecular properties and functional liabilities. The resulting multi-omics data sets and their interrogation facilitate both basic research mechanistic studies and translation of the findings into the clinic. In this review, we provide a list of resources and tools potentially valuable for lung cancer basic and translational research. Importantly, we point out pitfalls and caveats when performing computational analyses of these data sets and provide a vision of future computational biology developments that will aid lung cancer translational research.

Modeling Progression of Single Cell Populations Through the Cell Cycle as a Sequence of Switches

Cell cycle is a biological process underlying the existence and propagation of life in time and space. It has been an object for mathematical modeling for long, with several alternative mechanistic modeling principles suggested, describing in more or less details the known molecular mechanisms. Recently, cell cycle has been investigated at single cell level in snapshots of unsynchronized cell populations, exploiting the new methods for transcriptomic and proteomic molecular profiling. This raises a need for simplified semi-phenomenological cell cycle models, in order to formalize the processes underlying the cell cycle, at a higher abstracted level. Here we suggest a modeling framework, recapitulating the most important properties of the cell cycle as a limit trajectory of a dynamical process characterized by several internal states with switches between them. In the simplest form, this leads to a limit cycle trajectory, composed by linear segments in logarithmic coordinates describing some extensive (depending on system size) cell properties. We prove a theorem connecting the effective embedding dimensionality of the cell cycle trajectory with the number of its linear segments. We also develop a simplified kinetic model with piecewise-constant kinetic rates describing the dynamics of lumps of genes involved in S-phase and G2/M phases. We show how the developed cell cycle models can be applied to analyze the available single cell datasets and simulate certain properties of the observed cell cycle trajectories. Based on our model, we can predict with good accuracy the cell line doubling time from the length of cell cycle trajectory.

scMRMA: single cell multiresolution marker-based annotation

Single-cell RNA sequencing has become a powerful tool for identifying and characterizing cellular heterogeneity. One essential step to understanding cellular heterogeneity is determining cell identities. The widely used strategy predicts identities by projecting cells or cell clusters unidirectionally against a reference to find the best match. Here, we develop a bidirectional method, scMRMA, where a hierarchical reference guides iterative clustering and deep annotation with enhanced resolutions. Taking full advantage of the reference, scMRMA greatly improves the annotation accuracy. scMRMA achieved better performance than existing methods in four benchmark datasets and successfully revealed the expansion of CD8 T cell populations in squamous cell carcinoma after anti-PD-1 treatment.

Deep learning tackles single-cell analysis-a survey of deep learning for scRNA-seq analysis

Since its selection as the method of the year in 2013, single-cell technologies have become mature enough to provide answers to complex research questions. With the growth of single-cell profiling technologies, there has also been a significant increase in data collected from single-cell profilings, resulting in computational challenges to process these massive and complicated datasets. To address these challenges, deep learning (DL) is positioned as a competitive alternative for single-cell analyses besides the traditional machine learning approaches. Here, we survey a total of 25 DL algorithms and their applicability for a specific step in the single cell RNA-seq processing pipeline. Specifically, we establish a unified mathematical representation of variational autoencoder, autoencoder, generative adversarial network and supervised DL models, compare the training strategies and loss functions for these models, and relate the loss functions of these models to specific objectives of the data processing step. Such a presentation will allow readers to choose suitable algorithms for their particular objective at each step in the pipeline. We envision that this survey will serve as an important information portal for learning the application of DL for scRNA-seq analysis and inspire innovative uses of DL to address a broader range of new challenges in emerging multi-omics and spatial single-cell sequencing.

Hallmarks of Cancer: New Dimensions

The hallmarks of cancer conceptualization is a heuristic tool for distilling the vast complexity of cancer phenotypes and genotypes into a provisional set of underlying principles. As knowledge of cancer mechanisms has progressed, other facets of the disease have emerged as potential refinements. Herein, the prospect is raised that phenotypic plasticity and disrupted differentiation is a discrete hallmark capability, and that nonmutational epigenetic reprogramming and polymorphic microbiomes both constitute distinctive enabling characteristics that facilitate the acquisition of hallmark capabilities. Additionally, senescent cells, of varying origins, may be added to the roster of functionally important cell types in the tumor microenvironment. SIGNIFICANCE: Cancer is daunting in the breadth and scope of its diversity, spanning genetics, cell and tissue biology, pathology, and response to therapy. Ever more powerful experimental and computational tools and technologies are providing an avalanche of "big data" about the myriad manifestations of the diseases that cancer encompasses. The integrative concept embodied in the hallmarks of cancer is helping to distill this complexity into an increasingly logical science, and the provisional new dimensions presented in this perspective may add value to that endeavor, to more fully understand mechanisms of cancer development and malignant progression, and apply that knowledge to cancer medicine.

Transcriptional census of epithelial-mesenchymal plasticity in cancer

Epithelial-mesenchymal plasticity (EMP) contributes to tumor progression, promoting therapy resistance and immune cell evasion. Definitive molecular features of this plasticity have largely remained elusive due to the limited scale of most studies. Leveraging single-cell RNA sequencing data from 266 tumors spanning eight different cancer types, we identify expression patterns associated with intratumoral EMP. Integrative analysis of these programs confirmed a high degree of diversity among tumors. These diverse programs are associated with combinations of various common regulatory mechanisms initiated from cues within the tumor microenvironment. We show that inferring regulatory features can inform effective therapeutics to restrict EMP.

Discovering dominant tumor immune archetypes in a pan-cancer census

Cancers display significant heterogeneity with respect to tissue of origin, driver mutations, and other features of the surrounding tissue. It is likely that individual tumors engage common patterns of the immune system-here "archetypes"-creating prototypical non-destructive tumor immune microenvironments (TMEs) and modulating tumor-targeting. To discover the dominant immune system archetypes, the University of California, San Francisco (UCSF) Immunoprofiler Initiative (IPI) processed 364 individual tumors across 12 cancer types using standardized protocols. Computational clustering of flow cytometry and transcriptomic data obtained from cell sub-compartments uncovered dominant patterns of immune composition across cancers. These archetypes were profound insofar as they also differentiated tumors based upon unique immune and tumor gene-expression patterns. They also partitioned well-established classifications of tumor biology. The IPI resource provides a template for understanding cancer immunity as a collection of dominant patterns of immune organization and provides a rational path forward to learn how to modulate these to improve therapy.

Single-cell transcriptomics of neuroblastoma identifies chemoresistance-associated genes and pathways

High-Risk neuroblastoma (NB) survival rate is still <50%, despite treatments being more and more aggressive. The biggest hurdle liable to cancer therapy failure is the drug resistance by tumor cells that is likely due to the intra-tumor heterogeneity (ITH). To investigate the link between ITH and therapy resistance in NB, we performed a single cell RNA sequencing (scRNAseq) of etoposide and cisplatin resistant NB and their parental cells. Our analysis showed a clear separation of resistant and parental cells for both conditions by identifying 8 distinct tumor clusters in etoposide-resistant/parental and 7 in cisplatin-resistant/parental cells. We discovered that drug resistance can affect NB cell identities; highlighting the bi-directional ability of adrenergic-to-mesenchymal transition of NB cells. The biological processes driving the identified resistant cell subpopulations revealed genes such as (BARD1, BRCA1, PARP1, HISTH1 axis, members of RPL family), suggesting a potential drug resistance due to the acquisition of DNA repair mechanisms and to the modification of the drug targets. Deconvolution analysis of bulk RNAseq data from 498 tumors with cell subpopulation signatures showed that the transcriptional heterogeneity of our cellular models reflected the ITH of NB tumors and allowed the identification of clusters associated with worse/better survival.

Our study demonstrates the distinct cell populations characterized by genes involved in different biological processes can have a role in NB drug treatment failure. These findings evidence the importance of ITH in NB drug resistance studies and the chance that scRNA-seq analysis offers in the identification of genes and pathways liable for drug resistance.