Identifying and removing the cell-cycle effect from single-cell RNA-Sequencing data

Mitigating Cell Cycle Effects in Multi-Omics Data: Solutions and Analytical Frameworks

Cell cycle structures vary significantly across cell types, which exhibit distinct phase compositions. Asynchronous DNA replication and dynamic cellular characteristics during the cell cycle result in considerable heterogeneity in DNA dosage, chromatin accessibility, methylation, and expression. Nonetheless, the consequences of cell cycle disruption in the interpretation of multi-omics data remain unclear. Here, we systematically assessed the influence of distinct cell phase structures on the interpretation of omics features in proliferating cells, and proposed solutions for each omics dataset. For copy number variation (CNV) calling, asynchronous replication timing (RT) interference induces false CNVs in cells with high S-phase ratio (SPR), which are significantly decreased following replication timing domain (RTD) correction. Similar noise is observed in the chromatin accessibility data. Moreover, for DNA methylation and transcriptomic analyses, cell cycle-sorted data outperformed direct comparison in elucidating the biological features of compared cells. Additionally, we established an integrated pipeline to identify differentially expressed genes (DEGs) after cell cycle phasing. Consequently, our study demonstrated extensive cell-cycle heterogeneity, warranting consideration in future studies involving cells with diverse cell-cycle structures. RTD correction or phase-specific comparison could reduce the influence of cell cycle composition on the analysis of the differences observed between stem and differentiated cells.

A single-cell resolved genotype-phenotype map using genome-wide genetic and environmental perturbations

Heterogeneity is inherent to living organisms and it determines cell fate and phenotypic variability. Despite its ubiquity, the underlying molecular mechanisms and the genetic basis linking genotype to-phenotype heterogeneity remain a central challenge. Here we construct a yeast knockout library with a clone and genotype RNA barcoding structure suitable for genome-scale analyses to generate a high-resolution single-cell yeast transcriptome atlas of 3500 mutants under control and stress conditions. We find that transcriptional heterogeneity reflects the coordinated expression of specific gene programs, generating a continuous of cell states that can be responsive to external insults. Cell state plasticity can be genetically modulated with mutants that act as state attractors and disruption of state homeostasis results in decreased adaptive fitness. Leveraging on intra-genetic variability, we establish that regulators of transcriptional heterogeneity are functionally diverse and influenced by the environment. Our multimodal perturbation-based single-cell Genotype-to-Transcriptome Atlas in yeast provides insights into organism-level responses.

Inhibition of ATP-citrate lyase by bempedoic acid protects against abdominal aortic aneurysm formation in mice

Abdominal aortic aneurysm (AAA) is a prevalent degenerative disease characterized by an exacerbated inflammation and destructive vascular remodeling. Unfortunately, effective pharmacological tools for the treatment of this disease remain a challenge. ATP-citrate lyase (ACLY), the primary enzyme responsible for acetyl-CoA biosynthesis, is a key regulator of inflammatory signaling in macrophages and lymphocytes. Here, we found increased levels of the active (phosphorylated) form of ACLY (p-ACLY) in the inflammatory infiltrate of AAA from patients and in aneurysmal lesions from angiotensin II (Ang II)-infused apolipoprotein E-deficient mice (ApoE-/-). Furthermore, plasma ACLY levels positively correlates with IL6 and IFNγ levels in patients with AAA, while inflammatory stimuli strongly upregulated ACLY expression in macrophages and Jurkat cells. The administration of the ACLY inhibitor bempedoic acid (BemA) protected against Ang II-induced AAA formation in ApoE-/- mice, limiting the progression of aortic dilatation and reducing mortality due to aortic rupture. BMS-303141, another ACLY inhibitor, also ameliorated AAA formation, although to a lesser extent. BemA attenuated vascular remodeling and the disorganization and rupture of elastic fibers induced by Ang II, as well as vascular inflammation, decreasing the recruitment of macrophages (CD68 +) and neutrophils (Ly-6G+) into the aortic wall. Moreover, BemA shifted splenic monocytes toward a functionally anti-inflammatory phenotype, and increased the percentage of CD4+CD69+ cells. Taken together, these results support the contribution of ACLY to AAA and point to BemA as a promising tool to be considered for future clinical trials addressing the management of this disease which is quite often associated with disorders of lipoprotein metabolism.

In vivo regulation of the monocyte phenotype by Mycobacterium marinum and the ESX-1 type VII secretion system

Pathogenic mycobacteria require the conserved ESX-1 type VII secretion system to cause disease. In a murine Mycobacterium marinum infection model we previously demonstrated that infiltrating monocytes and neutrophils represent the major bacteria-harbouring cell populations in infected tissue. In the current study we use this model, in combination with scRNA sequencing, to assess the impact of M. marinum infection on the transcriptional profile of infiltrating Ly6C⁺MHCII⁺ monocytes in vivo. Our findings demonstrate that infection of infiltrating monocytes with M. marinum alters their cytokine expression profile, induces glycolytic metabolism, hypoxia-mediated signaling, nitric oxide synthesis, tissue remodeling, and suppresses responsiveness to IFNγ. We further show that the transcriptional response of bystander monocytes is influenced by ESX-1-dependent mechanisms, including a reduced responsiveness to IFNγ. These findings suggest that mycobacterial infection has pleiotropic effects on monocyte phenotype, with potential implications in bacterial growth restriction and granuloma formation.

Ginkgolide B increases healthspan and lifespan of female mice

Various anti-aging interventions show promise in extending lifespan, but many are ineffective or even harmful to healthspan. Ginkgolide B (GB), derived from Ginkgo biloba, reduces aging-related morbidities such as osteoporosis, yet its effects on healthspan and longevity have not been fully understood. In this study, we found that continuous oral administration of GB to female mice beginning at 20 months of age extended median survival and median lifespan by 30% and 8.5%, respectively. GB treatment also decreased tumor incidence; enhanced muscle quality, physical performance and metabolism; and reduced systemic inflammation and senescence. Single-nucleus RNA sequencing of skeletal muscle tissue showed that GB ameliorated aging-associated changes in cell type composition, signaling pathways and intercellular communication. GB reduced aging-induced Runx1+ type 2B myonuclei through the upregulation of miR-27b-3p, which suppresses Runx1 expression. Using functional analyses, we found that Runx1 promoted senescence and cell death in muscle cells. Collectively, these findings suggest the translational potential of GB to extend healthspan and lifespan and to promote healthy aging.

Engineering mesenchymal stem cells for premature ovarian failure: overcoming challenges and innovating therapeutic strategies

Premature ovarian failure (POF) is a leading cause of infertility in women, causing significant psychological and physical distress. Current therapeutic options are limited, necessitating the exploration of new treatments. Mesenchymal stem cells (MSCs), known for their remarkable homing and regenerative properties, have emerged as a promising intervention for POF. However, their clinical efficacy has been inconsistent. This paper aims to address these challenges by examining the cellular heterogeneity within MSC populations, which is crucial for identifying and selecting specific functional subpopulations for clinical applications. Understanding this heterogeneity can enhance therapeutic efficacy and ensure treatment stability. Additionally, this review comprehensively examines the literature on the effectiveness, safety, and ethical considerations of MSCs for ovarian regeneration, with a focus on preclinical and clinical trials. We also discuss potential strategies involving genetically and tissue-engineered MSCs. By integrating insights from these studies, we propose new directions for the design of targeted MSC treatments for POF and related disorders, potentially improving outcomes, addressing safety concerns, and expanding therapeutic options while ensuring ethical compliance.

What Is a Plant Cell Type in the Age of Single-Cell Biology? It's Complicated

One of the fundamental questions in developmental biology is how a cell is specified to differentiate as a specialized cell type. Traditionally, plant cell types were defined based on their function, location, morphology, and lineage. Currently, in the age of single-cell biology, researchers typically attempt to assign plant cells to cell types by clustering them based on their transcriptomes. However, because cells are dynamic entities that progress through the cell cycle and respond to signals, the transcriptome also reflects the state of the cell at a particular moment in time, raising questions about how to define a cell type. We suggest that these complexities and dynamics of cell states are of interest and further consider the roles signaling, stochasticity, cell cycle, and mechanical forces play in plant cell fate specification. Once established, cell identity must also be maintained. With the wealth of single-cell data coming out, the field is poised to elucidate both the complexity and dynamics of cell states.

A Cell Cycle-Aware Network for Data Integration and Label Transferring of Single-Cell RNA-Seq and ATAC-Seq

In recent years, the integration of single-cell multi-omics data has provided a more comprehensive understanding of cell functions and internal regulatory mechanisms from a non-single omics perspective, but it still suffers many challenges, such as omics-variance, sparsity, cell heterogeneity, and confounding factors. As it is known, the cell cycle is regarded as a confounder when analyzing other factors in single-cell RNA-seq data, but it is not clear how it will work on the integrated single-cell multi-omics data. Here, a cell cycle-aware network (CCAN) is developed to remove cell cycle effects from the integrated single-cell multi-omics data while keeping the cell type-specific variations. This is the first computational model to study the cell-cycle effects in the integration of single-cell multi-omics data. Validations on several benchmark datasets show the outstanding performance of CCAN in a variety of downstream analyses and applications, including removing cell cycle effects and batch effects of scRNA-seq datasets from different protocols, integrating paired and unpaired scRNA-seq and scATAC-seq data, accurately transferring cell type labels from scRNA-seq to scATAC-seq data, and characterizing the differentiation process from hematopoietic stem cells to different lineages in the integration of differentiation data.

scHiCyclePred: a deep learning framework for predicting cell cycle phases from single-cell Hi-C data using multi-scale interaction information

The emergence of single-cell Hi-C (scHi-C) technology has provided unprecedented opportunities for investigating the intricate relationship between cell cycle phases and the three-dimensional (3D) structure of chromatin. However, accurately predicting cell cycle phases based on scHi-C data remains a formidable challenge. Here, we present scHiCyclePred, a prediction model that integrates multiple feature sets to leverage scHi-C data for predicting cell cycle phases. scHiCyclePred extracts 3D chromatin structure features by incorporating multi-scale interaction information. The comparative analysis illustrates that scHiCyclePred surpasses existing methods such as Nagano_method and CIRCLET across various metrics including accuracy (ACC), F1 score, Precision, Recall, and balanced accuracy (BACC). In addition, we evaluate scHiCyclePred against the previously published CIRCLET using the dataset of complex tissues (Liu_dataset). Experimental results reveal significant improvements with scHiCyclePred exhibiting improvements of 0.39, 0.52, 0.52, and 0.39 over the CIRCLET in terms of ACC, F1 score, Precision, and Recall metrics, respectively. Furthermore, we conduct analyses on three-dimensional chromatin dynamics and gene features during the cell cycle, providing a more comprehensive understanding of cell cycle dynamics through chromatin structure. scHiCyclePred not only offers insights into cell biology but also holds promise for catalyzing breakthroughs in disease research. Access scHiCyclePred on GitHub at https:// github.com/HaoWuLab-Bioinformatics/ scHiCyclePred .

Increased heterogeneity in expression of genes associated with cancer progression and drug resistance

Fluctuations in the number of regulatory molecules and differences in timings of molecular events can generate variation in gene expression among genetically identical cells in the same environmental condition. This variation, termed as expression noise, can create differences in metabolic state and cellular functions, leading to phenotypic heterogeneity. Expression noise and phenotypic heterogeneity have been recognized as important contributors to intra-tumor heterogeneity, and have been associated with cancer growth, progression, and therapy resistance. However, how expression noise changes with cancer progression in actual cancer patients has remained poorly explored. Such an analysis, through identification of genes with increasing expression noise, can provide valuable insights into generation of intra-tumor heterogeneity, and could have important implications for understanding immune-suppression, drug tolerance and therapy resistance. In this work, we performed a genome-wide identification of changes in gene expression noise with cancer progression using single-cell RNA-seq data of lung adenocarcinoma patients at different stages of cancer. We identified 37 genes in epithelial cells that showed an increasing noise trend with cancer progression, many of which were also associated with cancer growth, EMT and therapy resistance. We found that expression of several of these genes was positively associated with expression of mitochondrial genes, suggesting an important role of mitochondria in generation of heterogeneity. In addition, we uncovered substantial differences in sample-specific noise profiles which could have implications for personalized prognosis and treatment.

A Cell Cycle-aware Network for Data Integration and Label Transferring of Single-cell RNA-seq and ATAC-seq

In recent years, the integration of single-cell multi-omics data has provided a more comprehensive understanding of cell functions and internal regulatory mechanisms from a non-single omics perspective, but it still suffers many challenges, such as omics-variance, sparsity, cell heterogeneity and confounding factors. As we know, cell cycle is regarded as a confounder when analyzing other factors in single-cell RNA-seq data, but it's not clear how it will work on the integrated single-cell multi-omics data. Here, we developed a Cell Cycle-Aware Network (CCAN) to remove cell cycle effects from the integrated single-cell multi-omics data while keeping the cell type-specific variations. This is the first computational model to study the cell-cycle effects in the integration of single-cell multi-omics data. Validations on several benchmark datasets show the out-standing performance of CCAN in a variety of downstream analyses and applications, including removing cell cycle effects and batch effects of scRNA-seq datasets from different protocols, integrating paired and unpaired scRNA-seq and scATAC-seq data, accurately transferring cell type labels from scRNA-seq to scATAC-seq data, and characterizing the differentiation process from hematopoietic stem cells to different lineages in the integration of differentiation data.

SiFT: uncovering hidden biological processes by probabilistic filtering of single-cell data

Cellular populations simultaneously encode multiple biological attributes, including spatial configuration, temporal trajectories, and cell-cell interactions. Some of these signals may be overshadowed by others and harder to recover, despite the great progress made to computationally reconstruct biological processes from single-cell data. To address this, we present SiFT, a kernel-based projection method for filtering biological signals in single-cell data, thus uncovering underlying biological processes. SiFT applies to a wide range of tasks, from the removal of unwanted variation in the data to revealing hidden biological structures. We demonstrate how SiFT enhances the liver circadian signal by filtering spatial zonation, recovers regenerative cell subpopulations in spatially-resolved liver data, and exposes COVID-19 disease-related cells, pathways, and dynamics by filtering healthy reference signals. SiFT performs the correction at the gene expression level, can scale to large datasets, and compares favorably to state-of-the-art methods.

From G1 to M: a comparative study of methods for identifying cell cycle phases

Accurate identification of cell cycle phases in single-cell RNA-sequencing (scRNA-seq) data is crucial for biomedical research. Many methods have been developed to tackle this challenge, employing diverse approaches to predict cell cycle phases. In this review article, we delve into the standard processes in identifying cell cycle phases within scRNA-seq data and present several representative methods for comparison. To rigorously assess the accuracy of these methods, we propose an error function and employ multiple benchmarking datasets encompassing human and mouse data. Our evaluation results reveal a key finding: the fit between the reference data and the dataset being analyzed profoundly impacts the effectiveness of cell cycle phase identification methods. Therefore, researchers must carefully consider the compatibility between the reference data and their dataset to achieve optimal results. Furthermore, we explore the potential benefits of incorporating benchmarking data with multiple known cell cycle phases into the analysis. Merging such data with the target dataset shows promise in enhancing prediction accuracy. By shedding light on the accuracy and performance of cell cycle phase prediction methods across diverse datasets, this review aims to motivate and guide future methodological advancements. Our findings offer valuable insights for researchers seeking to improve their understanding of cellular dynamics through scRNA-seq analysis, ultimately fostering the development of more robust and widely applicable cell cycle identification methods.

A cell abundance analysis based on efficient PAM clustering for a better understanding of the dynamics of endometrial remodelling

BACKGROUND

Single-cell RNA sequencing (scRNA-seq) is a powerful tool for investigating cell abundance changes during tissue regeneration and remodeling processes. Differential cell abundance supports the initial clustering of all cells; then, the number of cells per cluster and sample are evaluated, and the dependence of these counts concerning the phenotypic covariates of the samples is studied. Analysis heavily depends on the clustering method. Partitioning Around Medoids (PAM or k-medoids) represents a well-established clustering procedure that leverages the downstream interpretation of clusters by pinpointing real individuals in the dataset as cluster centers (medoids) without reducing dimensions. Of note, PAM suffers from high computational costs and memory requirements.

RESULTS

This paper proposes a method for differential abundance analysis using PAM as a clustering method and negative binomial regression as a statistical model to relate covariates to cluster/cell counts. We used this approach to study the differential cell abundance of human endometrial cell types throughout the natural secretory phase of the menstrual cycle. We developed a new R package -scellpam-, that incorporates an efficient parallel C++ implementation of PAM, and applied this package in this study. We compared the PAM-BS clustering method with other methods and evaluated both the computational aspects of its implementation and the quality of the classifications obtained using distinct published datasets with known subpopulations that demonstrate promising results.

CONCLUSIONS

The implementation of PAM-BS, included in the scellpam package, exhibits robust performance in terms of speed and memory usage compared to other related methods. PAM allowed quick and robust clustering of sets of cells with a size ranging from 70,000 to 300,000 cells. https://cran.r-project.org/web/packages/scellpam/index.html . Finally, our approach provides important new insights into the transient subpopulations associated with the fertile time frame when applied to the study of changes in the human endometrium during the secretory phase of the menstrual cycle.

Computational single cell oncology: state of the art

Single cell computational analysis has emerged as a powerful tool in the field of oncology, enabling researchers to decipher the complex cellular heterogeneity that characterizes cancer. By leveraging computational algorithms and bioinformatics approaches, this methodology provides insights into the underlying genetic, epigenetic and transcriptomic variations among individual cancer cells. In this paper, we present a comprehensive overview of single cell computational analysis in oncology, discussing the key computational techniques employed for data processing, analysis, and interpretation. We explore the challenges associated with single cell data, including data quality control, normalization, dimensionality reduction, clustering, and trajectory inference. Furthermore, we highlight the applications of single cell computational analysis, including the identification of novel cell states, the characterization of tumor subtypes, the discovery of biomarkers, and the prediction of therapy response. Finally, we address the future directions and potential advancements in the field, including the development of machine learning and deep learning approaches for single cell analysis. Overall, this paper aims to provide a roadmap for researchers interested in leveraging computational methods to unlock the full potential of single cell analysis in understanding cancer biology with the goal of advancing precision oncology. For this purpose, we also include a notebook that instructs on how to apply the recommended tools in the Preprocessing and Quality Control section.

NIPMAP: niche-phenotype mapping of multiplex histology data by community ecology

Advances in multiplex histology allow surveying millions of cells, dozens of cell types, and up to thousands of phenotypes within the spatial context of tissue sections. This leads to a combinatorial challenge in (a) summarizing the cellular and phenotypic architecture of tissues and (b) identifying phenotypes with interesting spatial architecture. To address this, we combine ideas from community ecology and machine learning into niche-phenotype mapping (NIPMAP). NIPMAP takes advantage of geometric constraints on local cellular composition imposed by the niche structure of tissues in order to automatically segment tissue sections into niches and their interfaces. Projecting phenotypes on niches and their interfaces identifies previously-reported and previously-unreported spatially-driven phenotypes, concisely summarizes the phenotypic architecture of tissues, and reveals fundamental properties of tissue architecture. NIPMAP is applicable to both protein and RNA multiplex histology of healthy and diseased tissue. An open-source R/Python package implements NIPMAP.

scPrisma infers, filters and enhances topological signals in single-cell data using spectral template matching

Single-cell RNA sequencing has been instrumental in uncovering cellular spatiotemporal context. This task is challenging as cells simultaneously encode multiple, potentially cross-interfering, biological signals. Here we propose scPrisma, a spectral computational method that uses topological priors to decouple, enhance and filter different classes of biological processes in single-cell data, such as periodic and linear signals. We apply scPrisma to the analysis of the cell cycle in HeLa cells, circadian rhythm and spatial zonation in liver lobules, diurnal cycle in Chlamydomonas and circadian rhythm in the suprachiasmatic nucleus in the brain. scPrisma can be used to distinguish mixed cellular populations by specific characteristics such as cell type and uncover regulatory networks and cell-cell interactions specific to predefined biological signals, such as the circadian rhythm. We show scPrisma's flexibility in incorporating prior knowledge, inference of topologically informative genes and generalization to additional diverse templates and systems. scPrisma can be used as a stand-alone workflow for signal analysis and as a prior step for downstream single-cell analysis.

Microfluidic single-cell measurements of oxidative stress as a function of cell cycle position

Single-cell measurements routinely demonstrate high levels of variation between cells, but fewer studies provide insight into the analytical and biological sources of this variation. This is particularly true of chemical cytometry, in which individual cells are lysed and their contents separated, compared to more established single-cell measurements of the genome and transcriptome. To characterize population-level variation and its sources, we analyzed oxidative stress levels in 1278 individual Dictyostelium discoideum cells as a function of exogenous stress level and cell cycle position. Cells were exposed to varying levels of oxidative stress via singlet oxygen generation using the photosensitizer Rose Bengal. Single-cell data reproduced the dose-response observed in ensemble measurements by CE-LIF, superimposed with high levels of heterogeneity. Through experiments and data analysis, we explored possible biological sources of this heterogeneity. No trend was observed between population variation and oxidative stress level, but cell cycle position was a major contributor to heterogeneity in oxidative stress. Cells synchronized to the same stage of cell division were less heterogeneous than unsynchronized cells (RSD of 37-51% vs 93%), and mitotic cells had higher levels of reactive oxygen species than interphase cells. While past research has proposed changes in cell size during the cell cycle as a source of biological noise, the measurements presented here use an internal standard to normalize for effects of cell volume, suggesting a more complex contribution of cell cycle to heterogeneity of oxidative stress.

The circular RNA Ataxia Telangiectasia Mutated regulates oxidative stress in smooth muscle cells in expanding abdominal aortic aneurysms

An abdominal aortic aneurysm (AAA) is a pathological widening of the aortic wall characterized by loss of smooth muscle cells (SMCs), extracellular matrix degradation, and local inflammation. This condition is often asymptomatic until rupture occurs, leading to high morbidity and mortality rates. Diagnosis is mostly accidental and the only currently available treatment option remains surgical intervention. Circular RNAs (circRNAs) represent a novel class of regulatory non-coding RNAs that originate from backsplicing. Their highly stable loop structure, combined with a remarkable enrichment in body fluids, make circRNAs promising disease biomarkers. We investigated the contribution of circRNAs to AAA pathogenesis and their potential application to improve AAA diagnostics. Gene expression analysis revealed the presence of deregulated circular transcripts stemming from AAA-relevant gene loci. Among these, the circRNA to the Ataxia Telangiectasia Mutated gene (cATM) was upregulated in human AAA specimens, in AAA-derived SMCs, and serum samples collected from aneurysm patients. In primary aortic SMCs, cATM increased upon angiotensin II and doxorubicin stimulation, while its silencing triggered apoptosis. Higher cATM levels made AAA-derived SMCs less vulnerable to oxidative stress, compared with control SMCs. These data suggest that cATM contributes to elicit an adaptive oxidative-stress response in SMCs and provides a reliable AAA disease signature.

Targeting long non-coding RNA NUDT6 enhances smooth muscle cell survival and limits vascular disease progression

Long non-coding RNAs (lncRNAs) orchestrate various biological processes and regulate the development of cardiovascular diseases. Their potential therapeutic benefit to tackle disease progression has recently been extensively explored. Our study investigates the role of lncRNA Nudix Hydrolase 6 (NUDT6) and its antisense target fibroblast growth factor 2 (FGF2) in two vascular pathologies: abdominal aortic aneurysms (AAA) and carotid artery disease. Using tissue samples from both diseases, we detected a substantial increase of NUDT6, whereas FGF2 was downregulated. Targeting Nudt6 in vivo with antisense oligonucleotides in three murine and one porcine animal model of carotid artery disease and AAA limited disease progression. Restoration of FGF2 upon Nudt6 knockdown improved vessel wall morphology and fibrous cap stability. Overexpression of NUDT6 in vitro impaired smooth muscle cell (SMC) migration, while limiting their proliferation and augmenting apoptosis. By employing RNA pulldown followed by mass spectrometry as well as RNA immunoprecipitation, we identified Cysteine and Glycine Rich Protein 1 (CSRP1) as another direct NUDT6 interaction partner, regulating cell motility and SMC differentiation. Overall, the present study identifies NUDT6 as a well-conserved antisense transcript of FGF2. NUDT6 silencing triggers SMC survival and migration and could serve as a novel RNA-based therapeutic strategy in vascular diseases.

The small and large intestine contain related mesenchymal subsets that derive from embryonic Gli1+ precursors

The intestinal lamina propria contains a diverse network of fibroblasts that provide key support functions to cells within their local environment. Despite this, our understanding of the diversity, location and ontogeny of fibroblasts within and along the length of the intestine remains incomplete. Here we show that the small and large intestinal lamina propria contain similar fibroblast subsets that locate in specific anatomical niches. Nevertheless, we find that the transcriptional profile of similar fibroblast subsets differs markedly between the small intestine and colon suggesting region specific functions. We perform in vivo transplantation and lineage-tracing experiments to demonstrate that adult intestinal fibroblast subsets, smooth muscle cells and pericytes derive from Gli1-expressing precursors present in embryonic day 12.5 intestine. Trajectory analysis of single cell RNA-seq datasets of E12.5 and adult mesenchymal cells suggest that adult smooth muscle cells and fibroblasts derive from distinct embryonic intermediates and that adult fibroblast subsets develop in a linear trajectory from CD81⁺ fibroblasts. Finally, we provide evidence that colonic subepithelial PDGFRα^hi fibroblasts comprise several functionally distinct populations that originate from an Fgfr2-expressing fibroblast intermediate. Our results provide insights into intestinal stromal cell diversity, location, function, and ontogeny, with implications for intestinal development and homeostasis.

Variability of fluorescence intensity distribution measured by flow cytometry is influenced by cell size and cell cycle progression

The distribution of fluorescence signals measured with flow cytometry can be influenced by several factors, including qualitative and quantitative properties of the used fluorochromes, optical properties of the detection system, as well as the variability within the analyzed cell population itself. Most of the single cell samples prepared from in vitrocultures or clinical specimens contain a variable cell cycle component. Cell cycle, together with changes in the cell size, are two of the factors that alter the functional properties of analyzed cells and thus affect the interpretation of obtained results. Here, we describe the association between cell cycle status and cell size, and the variability in the distribution of fluorescence intensity as determined with flow cytometry, at population scale. We show that variability in the distribution of background and specific fluorescence signals is related to the cell cycle state of the selected population, with the 10% low fluorescence signal fraction enriched mainly in cells in their G0/G1 cell cycle phase, and the 10% high fraction containing cells mostly in the G2/M phase. Therefore we advise using caution and additional experimental validation when comparing populations defined by fractions at both ends of fluorescence signal distribution to avoid biases caused by the effect of cell cycle and cell size.

Effective Preprocessing of Single-Cell RNA-Seq for Unravelling Alzheimer's Disease Signatures

The development in the field of biomedical technology has brought significant progress in the diagnosis and prediction of many complex diseases. Part of this development is the single-cell RNA sequencing analysis, which allows the study of a complex disease in great depth at the cellular level. Such analyses can decipher the mechanisms that cause complex diseases, such as Alzheimer's disease (AD). However, the increasing depth in the collection of single-cell RNA sequencing data implies, in addition to greater challenges, the production of a large amount of information, which needs careful analysis. Toward this direction, we examine the approach to single-cell RNA sequencing data through the development of an exploratory data analysis methodology. For this purpose, a combination of various tools is presented for their effective and efficient processing. At the same time, reference is made to the relevant biological concepts, the goals and challenges of the studies, and the workflows of sequencing, preprocessing, and analysis of the data. Our framework is applied to Alzheimer's disease data providing evidence that such data are quite complex while the appropriate preprocess step can boost the machine learning processes for identifying AD signatures.

Engineered exosomes targeting MYC reverse the proneural-mesenchymal transition and extend survival of glioblastoma

Dysregulated Myc signaling is a key oncogenic pathway in glioblastoma multiforme (GBM). Yet, effective therapeutic targeting of Myc continues to be challenging. Here, we demonstrate that exosomes generated from human bone marrow mesenchymal stem cells (MSCs) engineered to encapsulate siRNAs targeting Myc (iExo-Myc) localize to orthotopic GBM tumors in mice. Treatment of late stage GBM tumors with iExo-Myc inhibits proliferation and angiogenesis, suppresses tumor growth, and extends survival. Transcriptional profiling of tumors reveals that the mesenchymal transition and estrogen receptor signaling pathways are impacted by Myc inhibition. Single nuclei RNA sequencing (snRNA-seq) shows that iExo-Myc treatment induces transcriptional repression of multiple growth factor and interleukin signaling pathways, triggering a mesenchymal to proneural transition and shifting the cellular landscape of the tumor. These data confirm that Myc is an effective anti-glioma target and that iExo-Myc offers a feasible, readily translational strategy to inhibit challenging oncogene targets for the treatment of brain tumors.

Insights into protein post-translational modification landscapes of individual human cells by trapped ion mobility time-of-flight mass spectrometry

Single cell proteomics is a powerful tool with potential for markedly enhancing understanding of cellular processes. Here we report the development and application of multiplexed single cell proteomics using trapped ion mobility time-of-flight mass spectrometry. When employing a carrier channel to improve peptide signal, this method allows over 40,000 tandem mass spectra to be acquired in 30 min. Using a KRASG12C model human-derived cell line, we demonstrate the quantification of over 1200 proteins per cell with high relative sequence coverage permitting the detection of multiple classes of post-translational modifications in single cells. When cells were treated with a KRASG12C covalent inhibitor, this approach revealed cell-to-cell variability in the impact of the drug, providing insight missed by traditional proteomics. We provide multiple resources necessary for the application of single cell proteomics to drug treatment studies including tools to reduce cell cycle linked proteomic effects from masking pharmacological phenotypes.

Brain milieu induces early microglial maturation through the BAX-Notch axis

Microglia are derived from primitive myeloid cells and gain their early identity in the embryonic brains. However, the mechanism by which the brain milieu confers microglial maturation signature remains elusive. Here, we demonstrate that the baxcq55 zebrafish and Baxtm1Sjk mouse embryos exhibit similarly defective early microglial maturation. BAX, a typical pro-apoptotic factor, is highly enriched in neuronal cells and regulates microglial maturation through both pro-apoptotic and non-apoptotic mechanisms. BAX regulates dlb via the CaMKII-CREB axis calcium-dependently in living neurons while ensuring the efficient Notch activation in the immigrated pre-microglia by apoptotic neurons. Notch signaling is conserved in supporting embryonic microglia maturation. Compromised microglial development occurred in the Cx3cr1Cre/+Rbpjfl/fl embryonic mice; however, microglia acquire their appropriate signature when incubated with DLL3 in vitro. Thus, our findings elucidate a BAX-CaMKII-CREB-Notch network triggered by the neuronal milieu in microglial development, which may provide innovative insights for targeting microglia in neuronal disorder treatment.

Decoding brain memory formation by single-cell RNA sequencing

To understand how distinct memories are formed and stored in the brain is an important and fundamental question in neuroscience and computational biology. A population of neurons, termed engram cells, represents the physiological manifestation of a specific memory trace and is characterized by dynamic changes in gene expression, which in turn alters the synaptic connectivity and excitability of these cells. Recent applications of single-cell RNA sequencing (scRNA-seq) and single-nucleus RNA sequencing (snRNA-seq) are promising approaches for delineating the dynamic expression profiles in these subsets of neurons, and thus understanding memory-specific genes, their combinatorial patterns and regulatory networks. The aim of this article is to review and discuss the experimental and computational procedures of sc/snRNA-seq, new studies of molecular mechanisms of memory aided by sc/snRNA-seq in human brain diseases and related mouse models, and computational challenges in understanding the regulatory mechanisms underlying long-term memory formation.

Comparison of cell state models derived from single-cell RNA sequencing data: graph versus multi-dimensional space

Single-cell sequencing technologies have revolutionized molecular and cellular biology and stimulated the development of computational tools to analyze the data generated from these technology platforms. However, despite the recent explosion of computational analysis tools, relatively few mathematical models have been developed to utilize these data. Here we compare and contrast two cell state geometries for building mathematical models of cell state-transitions with single-cell RNA-sequencing data with hematopoeisis as a model system; (i) by using partial differential equations on a graph representing intermediate cell states between known cell types, and (ii) by using the equations on a multi-dimensional continuous cell state-space. As an application of our approach, we demonstrate how the calibrated models may be used to mathematically perturb normal hematopoeisis to simulate, predict, and study the emergence of novel cell states during the pathogenesis of acute myeloid leukemia. We particularly focus on comparing the strength and weakness of the graph model and multi-dimensional model.

Single-cell landscape of immunocytes in patients with extrahepatic cholangiocarcinoma

Background

The intricate landscape of immunocytes in the tumor microenvironment (TME) is fundamental to immunotherapy but notably under-researched in extrahepatic cholangiocarcinoma (ECCA).

Methods

Single-cell RNA sequencing technology was conducted to make an in-depth analysis of immunocytes from matched tumor tissues, paratumor tissues and peripheral blood from ECCA patients. The potential cellular interactions between two cell populations were analyzed with software CellPhoneDB (v2.1.7).

Results

We obtained 13526 cells and characterized the transcriptomes and heterogeneity of different clusters and subclusters of immunocytes from ECCA, including CD4+ T cells, CD8+ T cells, B cells and myeloid immunocytes. We observed the rarely described immunocyte subclusters "intermediate" exhausted CD8+ T (CD8+ Tex) cells and “nonclassic” plasmacytes (CD27+ CD138+ CD38−). In addition, we identified potential immunotherapy targets, for example, ACP5, MAGEH1, TNFRSF9 and CCR8 for Tregs and MT1 for CD8+ Tex cells. We also found strong cellular interactions among Treg cells, M2 macrophages and CD8+ Tex cells through ligand–receptor analysis, implying that potential cellular cross-linkage promoted the immunosuppressive nature of the TME.

Conclusions

In a word, our study illuminated the components of the TME and revealed potential cellular interactions at the individual cellular level in ECCA, we aimed to provide a new perspective for further immunological studies and immunotherapy of ECCA.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-022-03424-5.

An activation to memory differentiation trajectory of tumor-infiltrating lymphocytes informs metastatic melanoma outcomes

There is a need for better classification and understanding of tumor-infiltrating lymphocytes (TILs). Here, we applied advanced functional genomics to interrogate 9,000 human tumors and multiple single-cell sequencing sets using benchmarked T cell states, comprehensive T cell differentiation trajectories, human and mouse vaccine responses, and other human TILs. Compared with other T cell states, enrichment of T memory/resident memory programs was observed across solid tumors. Trajectory analysis of single-cell melanoma CD8+ TILs also identified a high fraction of memory/resident memory-scoring TILs in anti-PD-1 responders, which expanded post therapy. In contrast, TILs scoring highly for early T cell activation, but not exhaustion, associated with non-response. Late/persistent, but not early activation signatures, prognosticate melanoma survival, and co-express with dendritic cell and IFN-γ response programs. These data identify an activation-like state associated to poor response and suggest successful memory conversion, above resuscitation of exhaustion, is an under-appreciated aspect of successful anti-tumoral immunity.

Enhancing biological signals and detection rates in single-cell RNA-seq experiments with cDNA library equalization

Considerable effort has been devoted to refining experimental protocols to reduce levels of technical variability and artifacts in single-cell RNA-sequencing data (scRNA-seq). We here present evidence that equalizing the concentration of cDNA libraries prior to pooling, a step not consistently performed in single-cell experiments, improves gene detection rates, enhances biological signals, and reduces technical artifacts in scRNA-seq data. To evaluate the effect of equalization on various protocols, we developed Scaffold, a simulation framework that models each step of an scRNA-seq experiment. Numerical experiments demonstrate that equalization reduces variation in sequencing depth and gene-specific expression variability. We then performed a set of experiments in vitro with and without the equalization step and found that equalization increases the number of genes that are detected in every cell by 17-31%, improves discovery of biologically relevant genes, and reduces nuisance signals associated with cell cycle. Further support is provided in an analysis of publicly available data.

CCPE: cell cycle pseudotime estimation for single cell RNA-seq data

Pseudotime analysis from scRNA-seq data enables to characterize the continuous progression of various biological processes, such as the cell cycle. Cell cycle plays an important role in cell fate decisions and differentiation and is often regarded as a confounder in scRNA-seq data analysis when analyzing the role of other factors. Therefore, accurate prediction of cell cycle pseudotime and identification of cell cycle stages are important steps for characterizing the development-related biological processes. Here, we develop CCPE, a novel cell cycle pseudotime estimation method to characterize cell cycle timing and identify cell cycle phases from scRNA-seq data. CCPE uses a discriminative helix to characterize the circular process of the cell cycle and estimates each cell's pseudotime along the cell cycle. We evaluated the performance of CCPE based on a variety of simulated and real scRNA-seq datasets. Our results indicate that CCPE is an effective method for cell cycle estimation and competitive in various applications compared with other existing methods. CCPE successfully identified cell cycle marker genes and is robust to dropout events in scRNA-seq data. Accurate prediction of the cell cycle using CCPE can also effectively facilitate the removal of cell cycle effects across cell types or conditions.

Single-Cell RNA Sequencing and Assay for Transposase-Accessible Chromatin Using Sequencing Reveals Cellular and Molecular Dynamics of Aortic Aging in Mice

OBJECTIVE

The impact of vascular aging on cardiovascular diseases has been extensively studied; however, little is known regarding the cellular and molecular mechanisms underlying age-related vascular aging in aortic cellular subpopulations. Approach and Results: Transcriptomes and transposase-accessible chromatin profiles from the aortas of 4-, 26-, and 86-week-old C57/BL6J mice were analyzed using single-cell RNA sequencing and assay for transposase-accessible chromatin sequencing. By integrating the heterogeneous transcriptome and chromatin accessibility data, we identified cell-specific TF (transcription factor) regulatory networks and open chromatin states. We also determined that aortic aging affects cell interactions, inflammation, cell type composition, dysregulation of transcriptional control, and chromatin accessibility. Endothelial cells 1 have higher gene set activity related to cellular senescence and aging than do endothelial cells 2. Moreover, construction of senescence trajectories shows that endothelial cell 1 and fibroblast senescence is associated with distinct TF open chromatin states and an mRNA expression model.

CONCLUSIONS

Our data provide a system-wide model for transcriptional and epigenetic regulation during aortic aging at single-cell resolution.

Circular RNAs in cell cycle regulation: Mechanisms to clinical significance

Circular RNAs (circRNAs), a type of non‐coding RNA, are single‐stranded circularized molecules characterized by high abundance, evolutionary conservation and cell development‐ and tissue‐specific expression. A large body of studies has found that circRNAs exert a wide variety of functions in diverse biological processes, including cell cycle. The cell cycle is controlled by the coordinated activation and deactivation of cell cycle regulators. CircRNAs exert mutifunctional roles by regulating gene expression via various mechanisms. However, the functional relevance of circRNAs and cell cycle regulation largely remains to be elucidated. Herein, we briefly describe the biogenesis and mechanistic models of circRNAs and summarize their functions and mechanisms in the regulation of critical cell cycle modulators, including cyclins, cyclin‐dependent kinases and cyclin‐dependent kinase inhibitors. Moreover, we highlight the participation of circRNAs in cell cycle‐related signalling pathways and the clinical value of circRNAs as promising biomarkers or therapeutic targets in diseases related to cell cycle disorder.

Single-Cell RNA Sequencing Reveals Endothelial Cell Transcriptome Heterogeneity Under Homeostatic Laminar Flow

Supplemental Digital Content is available in the text.

Endothelial cells (ECs) that form the innermost layer of all vessels exhibit heterogeneous cell behaviors and responses to pro-angiogenic signals that are critical for vascular sprouting and angiogenesis. Once vessels form, remodeling and blood flow lead to EC quiescence, and homogeneity in cell behaviors and signaling responses. These changes are important for the function of mature vessels, but whether and at what level ECs regulate overall expression heterogeneity during this transition is poorly understood. Here, we profiled EC transcriptomic heterogeneity, and expression heterogeneity of selected proteins, under homeostatic laminar flow.

Single-cell RNA sequencing of freshly isolated bovine milk cells and cultured primary mammary epithelial cells

Bovine mammary function at molecular level is often studied using mammary tissue or primary bovine mammary epithelial cells (pbMECs). However, bulk tissue and primary cells are heterogeneous with respect to cell populations, adding further transcriptional variation in addition to genetic background. Thus, understanding of the variation in gene expression profiles of cell populations and their effect on function are limited. To investigate the mononuclear cell composition in bovine milk, we analyzed a single-cell suspension from a milk sample. Additionally, we harvested cultured pbMECs to characterize gene expression in a homogeneous cell population. Using the Drop-seq technology, we generated single-cell RNA datasets of somatic milk cells and pbMECs. The final datasets after quality control filtering contained 7,119 and 10,549 cells, respectively. The pbMECs formed 14 indefinite clusters displaying intrapopulation heterogeneity, whereas the milk cells formed 14 more distinct clusters. Our datasets constitute a molecular cell atlas that provides a basis for future studies of milk cell composition and gene expression, and could serve as reference datasets for milk cell analysis.

Advances and Opportunities in Single-Cell Transcriptomics for Plant Research

Single-cell approaches are quickly changing our view on biological systems by increasing the spatiotemporal resolution of our analyses to the level of the individual cell. The field of plant biology has fully embraced single-cell transcriptomics and is rapidly expanding the portfolio of available technologies and applications. In this review, we give an overview of the main advances in plant single-cell transcriptomics over the past few years and provide the reader with an accessible guideline covering all steps, from sample preparation to data analysis. We end by offering a glimpse of how these technologies will shape and accelerate plant-specific research in the near future.

Machine Intelligence in Single-Cell Data Analysis: Advances and New Challenges

The rapid development of single-cell technologies allows for dissecting cellular heterogeneity at different omics layers with an unprecedented resolution. In-dep analysis of cellular heterogeneity will boost our understanding of complex biological systems or processes, including cancer, immune system and chronic diseases, thereby providing valuable insights for clinical and translational research. In this review, we will focus on the application of machine learning methods in single-cell multi-omics data analysis. We will start with the pre-processing of single-cell RNA sequencing (scRNA-seq) data, including data imputation, cross-platform batch effect removal, and cell cycle and cell-type identification. Next, we will introduce advanced data analysis tools and methods used for copy number variance estimate, single-cell pseudo-time trajectory analysis, phylogenetic tree inference, cell-cell interaction, regulatory network inference, and integrated analysis of scRNA-seq and spatial transcriptome data. Finally, we will present the latest analyzing challenges, such as multi-omics integration and integrated analysis of scRNA-seq data.

Cell identity specification in plants: lessons from flower development

Multicellular organisms display a fascinating complexity of cellular identities and patterns of diversification. The concept of 'cell type' aims to describe and categorize this complexity. In this review, we discuss the traditional concept of cell types and highlight the impact of single-cell technologies and spatial omics on the understanding of cellular differentiation in plants. We summarize and compare position-based and lineage-based mechanisms of cell identity specification using flower development as a model system. More than understanding ontogenetic origins of differentiated cells, an important question in plant science is to understand their position- and developmental stage-specific heterogeneity. Combinatorial action and crosstalk of external and internal signals is the key to cellular heterogeneity, often converging on transcription factors that orchestrate gene expression programs.

Quantifying the effect of experimental perturbations at single-cell resolution

Current methods for comparing single-cell RNA sequencing datasets collected in multiple conditions focus on discrete regions of the transcriptional state space, such as clusters of cells. Here we quantify the effects of perturbations at the single-cell level using a continuous measure of the effect of a perturbation across the transcriptomic space. We describe this space as a manifold and develop a relative likelihood estimate of observing each cell in each of the experimental conditions using graph signal processing. This likelihood estimate can be used to identify cell populations specifically affected by a perturbation. We also develop vertex frequency clustering to extract populations of affected cells at the level of granularity that matches the perturbation response. The accuracy of our algorithm at identifying clusters of cells that are enriched or depleted in each condition is, on average, 57% higher than the next-best-performing algorithm tested. Gene signatures derived from these clusters are more accurate than those of six alternative algorithms in ground truth comparisons.

Toward uncharted territory of cellular heterogeneity: advances and applications of single-cell RNA-seq

Among single-cell analysis technologies, single-cell RNA-seq (scRNA-seq) has been one of the front runners in technical inventions. Since its induction, scRNA-seq has been well received and undergone many fast-paced technical improvements in cDNA synthesis and amplification, processing and alignment of next generation sequencing reads, differentially expressed gene calling, cell clustering, subpopulation identification, and developmental trajectory prediction. scRNA-seq has been exponentially applied to study global transcriptional profiles in all cell types in humans and animal models, healthy or with diseases, including cancer. Accumulative novel subtypes and rare subpopulations have been discovered as potential underlying mechanisms of stochasticity, differentiation, proliferation, tumorigenesis, and aging. scRNA-seq has gradually revealed the uncharted territory of cellular heterogeneity in transcriptomes and developed novel therapeutic approaches for biomedical applications. This review of the advancement of scRNA-seq methods provides an exploratory guide of the quickly evolving technical landscape and insights of focused features and strengths in each prominent area of progress.

The transcriptome dynamics of single cells during the cell cycle

The cell cycle is among the most basic phenomena in biology. Despite advances in single-cell analysis, dynamics and topology of the cell cycle in high-dimensional gene expression space remain largely unknown. We developed a linear analysis of transcriptome data which reveals that cells move along a planar circular trajectory in transcriptome space during the cycle. Non-cycling gene expression adds a third dimension causing helical motion on a cylinder. We find in immortalized cell lines that cell cycle transcriptome dynamics occur largely independently from other cellular processes. We offer a simple method ("Revelio") to order unsynchronized cells in time. Precise removal of cell cycle effects from the data becomes a straightforward operation. The shape of the trajectory implies that each gene is upregulated only once during the cycle, and only two dynamic components represented by groups of genes drive transcriptome dynamics. It indicates that the cell cycle has evolved to minimize changes of transcriptional activity and the related regulatory effort. This design principle of the cell cycle may be of relevance to many other cellular differentiation processes.

An era of single-cell genomics consortia

The human body consists of 37 trillion single cells represented by over 50 organs that are stitched together to make us who we are, yet we still have very little understanding about the basic units of our body: what cell types and states make up our organs both compositionally and spatially. Previous efforts to profile a wide range of human cell types have been attempted by the FANTOM and GTEx consortia. Now, with the advancement in genomic technologies, profiling the human body at single-cell resolution is possible and will generate an unprecedented wealth of data that will accelerate basic and clinical research with tangible applications to future medicine. To date, several major organs have been profiled, but the challenges lie in ways to integrate single-cell genomics data in a meaningful way. In recent years, several consortia have begun to introduce harmonization and equity in data collection and analysis. Herein, we introduce existing and nascent single-cell genomics consortia, and present benefits to necessitate single-cell genomic consortia in a regional environment to achieve the universal human cell reference dataset.

Single-Cell RNA-Seq of the Developing Cardiac Outflow Tract Reveals Convergent Development of the Vascular Smooth Muscle Cells

Cardiac outflow tract (OFT) is a major hotspot for congenital heart diseases. A thorough understanding of the cellular diversity, transitions, and regulatory networks of normal OFT development is essential to decipher the etiology of OFT malformations. We performed single-cell transcriptomic sequencing of 55,611 mouse OFT cells from three developmental stages that generally correspond to the early, middle, and late stages of OFT remodeling and septation. Known cellular transitions, such as endothelial-to-mesenchymal transition, have been recapitulated. In particular, we identified convergent development of the vascular smooth muscle cell (VSMC) lineage where intermediate cell subpopulations were found to be involved in either myocardial-to-VSMC trans-differentiation or mesenchymal-to-VSMC transition. Finally, we uncovered transcriptional regulators potentially governing cellular transitions. Our study provides a single-cell reference map of cell states for normal OFT development and paves the way for further studies of the etiology of OFT malformations at the single-cell level.

Novel classifier orthologs of bovine and human oocytes matured in different melatonin environments

It has been demonstrated that melatonin influences the developmental competence of both in vivo and in vitro matured oocytes. It modulates oocyte-specific gene expression patterns among mammalian species. Due to differences among study systems, the identification of the classifier orthologs-the homologous genes related among mammals that could universally categorize oocytes matured in environments with varied melatonin levels is still limitedly studied. To gain insight into such orthologs, cross-species transcription profiling meta-analysis of in vitro matured bovine oocytes and in vivo matured human oocytes in low and high melatonin environments was demonstrated in the current study. RNA-Seq data of bovine and human oocytes were retrieved from the Sequence Read Archive database and pre-processed. The used datasets of bovine oocytes obtained from culturing in the absence of melatonin and human oocytes from old patients were regarded as oocytes in the low melatonin environment (Low). Datasets from bovine oocytes cultured in 10-9 M melatonin and human oocytes from young patients were considered as oocytes in the high melatonin environment (High). Candidate orthologs differentially expressed between Low and High melatonin environments were selected by a linear model, and were further verified by Zero-inflated regression analysis. Support Vector Machine (SVM) was applied to determine the potentials of the verified orthologs as classifiers of melatonin environments. According to the acquired results, linear model analysis identified 284 candidate orthologs differentially expressed between Low and High melatonin environments. Among them, only 15 candidate orthologs were verified by Zero-inflated regression analysis (FDR ≤ 0.05). Utilization of the verified orthologs as classifiers in SVM resulted in the precise classification of oocyte learning datasets according to their melatonin environments (Misclassification rates < 0.18, area under curves > 0.9). In conclusion, the cross-species RNA-Seq meta-analysis to identify novel classifier orthologs of matured oocytes under different melatonin environments was successfully demonstrated in this study-delivering candidate orthologs for future studies at biological levels. Such verified orthologs might provide valuable evidence about melatonin sufficiency in target oocytes-by which, the decision on melatonin supplementation could be implied.

Characterizing and inferring quantitative cell cycle phase in single-cell RNA-seq data analysis

Cellular heterogeneity in gene expression is driven by cellular processes, such as cell cycle and cell-type identity, and cellular environment such as spatial location. The cell cycle, in particular, is thought to be a key driver of cell-to-cell heterogeneity in gene expression, even in otherwise homogeneous cell populations. Recent advances in single-cell RNA-sequencing (scRNA-seq) facilitate detailed characterization of gene expression heterogeneity and can thus shed new light on the processes driving heterogeneity. Here, we combined fluorescence imaging with scRNA-seq to measure cell cycle phase and gene expression levels in human induced pluripotent stem cells (iPSCs). By using these data, we developed a novel approach to characterize cell cycle progression. Although standard methods assign cells to discrete cell cycle stages, our method goes beyond this and quantifies cell cycle progression on a continuum. We found that, on average, scRNA-seq data from only five genes predicted a cell's position on the cell cycle continuum to within 14% of the entire cycle and that using more genes did not improve this accuracy. Our data and predictor of cell cycle phase can directly help future studies to account for cell cycle-related heterogeneity in iPSCs. Our results and methods also provide a foundation for future work to characterize the effects of the cell cycle on expression heterogeneity in other cell types.

Latent periodic process inference from single-cell RNA-seq data

The development of a phenotype in a multicellular organism often involves multiple, simultaneously occurring biological processes. Advances in single-cell RNA-sequencing make it possible to infer latent developmental processes from the transcriptomic profiles of cells at various developmental stages. Accurate characterization is challenging however, particularly for periodic processes such as cell cycle. To address this, we develop Cyclum, an autoencoder approach identifying circular trajectories in the gene expression space. Cyclum substantially improves the accuracy and robustness of cell-cycle characterization beyond existing approaches. Applying Cyclum to removing cell-cycle effects substantially improves delineations of cell subpopulations, which is useful for establishing various cell atlases and studying tumor heterogeneity.

rCASC: reproducible classification analysis of single-cell sequencing data

Background

Single-cell RNA sequencing is essential for investigating cellular heterogeneity and highlighting cell subpopulation-specific signatures. Single-cell sequencing applications have spread from conventional RNA sequencing to epigenomics, e.g., ATAC-seq. Many related algorithms and tools have been developed, but few computational workflows provide analysis flexibility while also achieving functional (i.e., information about the data and the tools used are saved as metadata) and computational reproducibility (i.e., a real image of the computational environment used to generate the data is stored) through a user-friendly environment.

Findings

rCASC is a modular workflow providing an integrated analysis environment (from count generation to cell subpopulation identification) exploiting Docker containerization to achieve both functional and computational reproducibility in data analysis. Hence, rCASC provides preprocessing tools to remove low-quality cells and/or specific bias, e.g., cell cycle. Subpopulation discovery can instead be achieved using different clustering techniques based on different distance metrics. Cluster quality is then estimated through the new metric "cell stability score" (CSS), which describes the stability of a cell in a cluster as a consequence of a perturbation induced by removing a random set of cells from the cell population. CSS provides better cluster robustness information than the silhouette metric. Moreover, rCASC's tools can identify cluster-specific gene signatures.

Conclusions

rCASC is a modular workflow with new features that could help researchers define cell subpopulations and detect subpopulation-specific markers. It uses Docker for ease of installation and to achieve a computation-reproducible analysis. A Java GUI is provided to welcome users without computational skills in R.

Deciphering the transcriptomic landscape of tumor-infiltrating CD8 lymphocytes in B16 melanoma tumors with single-cell RNA-Seq

Recent studies have proposed that tumor-specific tumor-infiltrating CD8⁺ T lymphocytes (CD8 TIL) can be classified into two main groups: "exhausted" TILs, characterized by high expression of the inhibitory receptors PD-1 and TIM-3 and lack of transcription factor 1 (Tcf1); and "memory-like" TILs, with self-renewal capacity and co-expressing Tcf1 and PD-1. However, a comprehensive definition of the heterogeneity existing within CD8 TILs has yet to be clearly established. To investigate this heterogeneity at the transcriptomic level, we performed paired single-cell RNA and TCR sequencing of CD8 T cells infiltrating B16 murine melanoma tumors, including cells of known tumor specificity. Unsupervised clustering and gene-signature analysis revealed four distinct CD8 TIL states - exhausted, memory-like, naïve and effector memory-like (EM-like) - and predicted novel markers, including Ly6C for the EM-like cells, that were validated by flow cytometry. Tumor-specific PMEL T cells were predominantly found within the exhausted and memory-like states but also within the EM-like state. Further, T cell receptor sequencing revealed a large clonal expansion of exhausted, memory-like and EM-like cells with partial clonal relatedness between them. Finally, meta-analyses of public bulk and single-cell RNA-seq data suggested that anti-PD-1 treatment induces the expansion of EM-like cells. Our reference map of the transcriptomic landscape of murine CD8 TILs will help interpreting future bulk and single-cell transcriptomic studies and may guide the analysis of CD8IL subpopulations in response to therapeutic interventions.