SLICE: determining cell differentiation and lineage based on single cell entropy

Single-nucleus transcriptomics reveal cardiac cell type-specific diversification in metabolic disease transgenic pigs

Cardiac damage is widely present in patients with metabolic diseases, but the exact pathophysiological mechanisms involved remain unclear. The porcine heart is an ideal material for cardiovascular research due to its similarities to the human heart. This study evaluated pathological features and performed single-nucleus RNA sequencing (snRNA-seq) on myocardial samples from both wild-type and metabolic disease-susceptible transgenic pigs (previously established). We found that transgenic pigs exhibited lipid metabolism disturbances and myocardial injury after a high-fat high-sucrose diet intervention. snRNA-seq reveals the cellular landscape of healthy and metabolically disturbed pig hearts and identifies the major cardiac cell populations affected by metabolic diseases. Within metabolic disorder hearts, metabolically active cardiomyocytes exhibited impaired function and reduced abundance. Moreover, massive numbers of reparative LYVE1+ macrophages were lost. Additionally, proinflammatory endothelial cells were activated with high expression of multiple proinflammatory cytokines. Our findings provide insights into the cellular mechanisms of metabolic disease-induced myocardial injury.

Deciphering the molecular landscape: evolutionary progression from gynecomastia to aggressive male breast cancer

BACKGROUND

Gynecomastia denotes the benign proliferation of glandular breast tissue and stands as a recognized risk factor for male breast cancer. Nonetheless, the underlying carcinogenic mechanisms orchestrating the progression from gynecomastia to cancer remain poorly understood.

METHODS

This study employed single-cell RNA sequencing (scRNA-seq) to meticulously dissect the cellular landscape of gynecomastia and unravel potential associations with male breast cancer at a single-cell resolution. Pseudotime and evolutionary analyses were executed to delineate the distinct features characterizing gynecomastia and male breast cancer. The TCGA database, along with cell-cell communication analysis and immunohistochemistry staining, was harnessed to validate differential gene expression, specifically focusing on CD13.

RESULT

From the copy number variation profiles and evolutionary tree, we inferred shared mutation characteristics (18p+ and 18q+) underpinning both conditions. The developmental trajectory unveiled an intriguing overlap between gynecomastia and malignant epithelial cells. Moreover, the differential gene CD13 emerged as a common denominator in both gynecomastia and male breast cancer when compared with normal mammary tissue. Cell-cell interaction analysis and communication dynamics within the tumor microenvironment spotlighted distinctions between CD13+ and CD13- subsets, with the former exhibiting elevated expression of FGFR1-FGF7.

CONCLUSIONS

Our investigation provides novel insights into the evolutionary progression from gynecomastia to male breast cancer, shedding light on the pivotal role of CD13 in driving this transition. The identification of CD13 as a potential therapeutic target suggests the feasibility of CD13-targeted interventions, specifically tailored for male breast cancer treatment.

Transcriptional programs mediating neuronal toxicity and altered glial-neuronal signaling in a Drosophila knock-in tauopathy model

Missense mutations in the gene encoding the microtubule-associated protein TAU (current and approved symbol is MAPT) cause autosomal dominant forms of frontotemporal dementia. Multiple models of frontotemporal dementia based on transgenic expression of human TAU in experimental model organisms, including Drosophila, have been described. These models replicate key features of the human disease but do not faithfully recreate the genetic context of the human disorder. Here we use CRISPR-Cas-mediated gene editing to model frontotemporal dementia caused by the TAU P301L mutation by creating the orthologous mutation, P251L, in the endogenous Drosophila tau gene. Flies heterozygous or homozygous for Tau P251L display age-dependent neurodegeneration, display metabolic defects, and accumulate DNA damage in affected neurons. To understand the molecular events promoting neuronal dysfunction and death in knock-in flies, we performed single-cell RNA sequencing on approximately 130,000 cells from brains of Tau P251L mutant and control flies. We found that expression of disease-associated mutant tau altered gene expression cell autonomously in all neuronal cell types identified. Gene expression was also altered in glial cells, suggestive of non-cell-autonomous regulation. Cell signaling pathways, including glial-neuronal signaling, were broadly dysregulated as were brain region and cell type-specific protein interaction networks and gene regulatory programs. In summary, we present here a genetic model of tauopathy that faithfully recapitulates the genetic context and phenotypic features of the human disease, and use the results of comprehensive single-cell sequencing analysis to outline pathways of neurotoxicity and highlight the potential role of non-cell-autonomous changes in glia.

Single-cell profiling of response to neoadjuvant chemo-immunotherapy in surgically resectable esophageal squamous cell carcinoma

BACKGROUND

The efficacy of neoadjuvant chemo-immunotherapy (NAT) in esophageal squamous cell carcinoma (ESCC) is challenged by the intricate interplay within the tumor microenvironment (TME). Unveiling the immune landscape of ESCC in the context of NAT could shed light on heterogeneity and optimize therapeutic strategies for patients.

METHODS

We analyzed single cells from 22 baseline and 24 post-NAT treatment samples of stage II/III ESCC patients to explore the association between the immune landscape and pathological response to neoadjuvant anti-PD-1 combination therapy, including pathological complete response (pCR), major pathological response (MPR), and incomplete pathological response (IPR).

RESULTS

Single-cell profiling identified 14 major cell subsets of cancer, immune, and stromal cells. Trajectory analysis unveiled an interesting link between cancer cell differentiation and pathological response to NAT. ESCC tumors enriched with less differentiated cancer cells exhibited a potentially favorable pathological response to NAT, while tumors enriched with clusters of more differentiated cancer cells may resist treatment. Deconvolution of transcriptomes in pre-treatment tumors identified gene signatures in response to NAT contributed by specific immune cell populations. Upregulated genes associated with better pathological responses in CD8 + effector T cells primarily involved interferon-gamma (IFNγ) signaling, neutrophil degranulation, and negative regulation of the T cell apoptotic process, whereas downregulated genes were dominated by those in the immune response-activating cell surface receptor signaling pathway. Natural killer cells in pre-treatment tumors from pCR patients showed a similar upregulation of gene expression in response to IFNγ but a downregulation of genes in the neutrophil-mediated immunity pathways. A decreased cellular contexture of regulatory T cells in ESCC TME indicated a potentially favorable pathological response to NAT. Cell-cell communication analysis revealed extensive interactions between CCL5 and its receptor CCR5 in various immune cells of baseline pCR tumors. Immune checkpoint interaction pairs, including CTLA4-CD86, TIGIT-PVR, LGALS9-HAVCR2, and TNFSF4-TNFRSF4, might serve as additional therapeutic targets for ICI therapy in ESCC.

CONCLUSIONS

This pioneering study unveiled an intriguing association between cancer cell differentiation and pathological response in esophageal cancer patients, revealing distinct subgroups of tumors for which neoadjuvant chemo-immunotherapy might be effective. We also delineated the immune landscape of ESCC tumors in the context of clinical response to NAT, which provides clinical insights for better understanding how patients respond to the treatment and further identifying novel therapeutic targets for ESCC patients in the future.

Mapping single-cell developmental potential in health and disease with interpretable deep learning

Single-cell RNA sequencing (scRNA-seq) has transformed our understanding of cell fate in developmental systems. However, identifying the molecular hallmarks of potency - the capacity of a cell to differentiate into other cell types - has remained challenging. Here, we introduce CytoTRACE 2, an interpretable deep learning framework for characterizing potency and differentiation states on an absolute scale from scRNA-seq data. Across 31 human and mouse scRNA-seq datasets encompassing 28 tissue types, CytoTRACE 2 outperformed existing methods for recovering experimentally determined potency levels and differentiation states covering the entire range of cellular ontogeny. Moreover, it reconstructed the temporal hierarchy of mouse embryogenesis across 62 timepoints; identified pan-tissue expression programs that discriminate major potency levels; and facilitated discovery of cellular phenotypes in cancer linked to survival and immunotherapy resistance. Our results illuminate a fundamental feature of cell biology and provide a broadly applicable platform for delineating single-cell differentiation landscapes in health and disease.

Transcriptional programs mediating neuronal toxicity and altered glial-neuronal signaling in a Drosophila knock-in tauopathy model

Missense mutations in the gene encoding the microtubule-associated protein tau cause autosomal dominant forms of frontotemporal dementia. Multiple models of frontotemporal dementia based on transgenic expression of human tau in experimental model organisms, including Drosophila, have been described. These models replicate key features of the human disease, but do not faithfully recreate the genetic context of the human disorder. Here we use CRISPR-Cas mediated gene editing to model frontotemporal dementia caused by the tau P301L mutation by creating the orthologous mutation, P251L, in the endogenous Drosophila tau gene. Flies heterozygous or homozygous for tau P251L display age-dependent neurodegeneration, metabolic defects and accumulate DNA damage in affected neurons. To understand the molecular events promoting neuronal dysfunction and death in knock-in flies we performed single-cell RNA sequencing on approximately 130,000 cells from brains of tau P251L mutant and control flies. We found that expression of disease-associated mutant tau altered gene expression cell autonomously in all neuronal cell types identified and non-cell autonomously in glial cells. Cell signaling pathways, including glial-neuronal signaling, were broadly dysregulated as were brain region and cell-type specific protein interaction networks and gene regulatory programs. In summary, we present here a genetic model of tauopathy, which faithfully recapitulates the genetic context and phenotypic features of the human disease and use the results of comprehensive single cell sequencing analysis to outline pathways of neurotoxicity and highlight the role of non-cell autonomous changes in glia.

Identification of female-enriched and disease-associated microglia (FDAMic) contributes to sexual dimorphism in late-onset Alzheimer's disease

BACKGROUND

Late-onset Alzheimer's disease (LOAD) is the most common form of dementia; it disproportionally affects women in terms of both incidence rates and severity of progression. The cellular and molecular mechanisms underlying this clinical phenomenon remain elusive and ill-defined.

METHODS

In-depth analyses were performed with multiple human LOAD single-nucleus transcriptome datasets to thoroughly characterize cell populations in the cerebral cortex. ROSMAP bulk human brain tissue transcriptome and DNA methylome datasets were also included for validation. Detailed assessments of microglial cell subpopulations and their relevance to sex-biased changes at the tissue level were performed. Clinical trait associations, cell evolutionary trajectories, and transcription regulon analyses were conducted.

RESULTS

The relative numbers of functionally defective microglia were aberrantly increased uniquely among affected females. Substratification of the microglia into different subtypes according to their transcriptomic signatures identified a group of female-enriched and disease-associated microglia (FDAMic), the numbers of which were positively associated with disease severity. Phenotypically, these cells exhibit transcriptomic signatures that support active proliferation, MHC class II autoantigen presentation and amyloid-β binding, but they are also likely defective in phagocytosis. FDAMic are likely evolved from female activated response microglia (ARMic) with an APOE4 background and compromised estrogen receptor (ER) signaling that is deemed to be active among most subtypes of microglia.

CONCLUSION

This study offered important insights at both the cellular and molecular levels into how ER signaling affects microglial heterogeneity and function. FDAMic are associated with more advanced pathologies and severe trends of cognitive decline. Their emergence could, at least in part, explain the phenomenon of greater penetrance of the APOE4 genotype found in females. The biases of FDAMic emergence toward female sex and APOE4 status may also explain why hormone replacement therapy is more effective in APOE4 carriers. The pathologic nature of FDAMic suggests that selective modulations of these cells may help to regain brain neuroimmune homeostasis, serving as a new target for future drug development.

StemnesScoRe: an R package to estimate the stemness of glioma cancer cells at single-cell resolution

Background/aim

Glioblastoma is the most heterogeneous and the most difficult-to-treat type of brain tumor and one of the deadliest among all cancers. The high plasticity of glioma cancer stem cells and the resistance they develop against multiple modalities of therapy, along with their high heterogeneity, are the main challenges faced during treatment of glioblastoma. Therefore, a better understanding of the stemness characteristics of glioblastoma cells is needed. With the development of various single-cell technologies and increasing applications of machine learning, indices based on transcriptomic and/or epigenomic data have been developed to quantitatively measure cellular states and stemness. In this study, we aimed to develop a glioma-specific stemness score model using scATAC-seq data for the first time.

Materials and methods

We first applied three powerful machine-learning algorithms, i.e. random forest, gradient boosting, and extreme gradient boosting, to glioblastoma scRNA-seq data to discover the most important genes associated with cellular states. We then identified promoter and enhancer regions associated with these genes. After downloading the scATAC-seq peaks and their read counts for each patient, we identified the overlapping regions between the single-cell peaks and the peaks of genes obtained through machine-learning algorithms. Then we calculated read counts that were mapped to these overlapping regions. We finally developed a model capable of estimating the stemness score for each glioma cell using overlapping regions and the importance of genes predictive of glioblastoma cellular states. We also created an R package, accessible to all researchers regardless of their coding proficiency.

Results

Our results showed that mesenchymal-like stem cells display higher stemness scores compared to neural-progenitor-, oligodendrocyte-progenitor-, and astrocyte-like cells.

Conclusion

scATAC-seq can be used to assess heterogeneity in glioblastoma and identify cells with high stemness characteristics. The package is publicly available at https://github.com/Necla/StemnesScoRe and includes documentation with implementation of a real-data experiment.

Deciphering the decisive factors driving fate bifurcations in somatic cell reprogramming

Single-cell studies have demonstrated that somatic cell reprogramming is a continuous process of cell fates transition. Only partial reprogramming intermediates can overcome the molecular bottlenecks to acquire pluripotency. To decipher the underlying decisive factors driving cell fate, we identified induced pluripotent stem cells or stromal-like cells (iPSCs/SLCs) and iPSCs or trophoblast-like cells (iPSCs/TLCs) fate bifurcations by reconstructing cellular trajectory. The mesenchymal-epithelial transition and the activation of pluripotency networks are the main molecular series in successful reprogramming. Correspondingly, intermediates diverge into SLCs accompanied by the inhibition of cell cycle genes and the activation of extracellular matrix genes, whereas the TLCs fate is characterized by the up-regulation of placenta development genes. Combining putative gene regulatory networks, seven (Taf7, Ezh2, Klf2, etc.) and three key factors (Cdc5l, Klf4, and Nanog) were individually identified as drivers of the successful reprogramming by triggering downstream pluripotent networks during iPSCs/SLCs and iPSCs/TLCs fate bifurcation. Conversely, 11 factors (Cebpb, Sox4, Junb, etc.) and four factors (Gata2, Jund, Ctnnb1, etc.) drive SLCs fate and TLCs fate, respectively. Our study sheds new light on the understanding of decisive factors driving cell fate, which is helpful for improving reprogramming efficiency through manipulating cell fates to avoid alternative fates.

Mapping the Single-cell Differentiation Landscape of Osteosarcoma

The genetic and intratumoral heterogeneity observed in human osteosarcomas (OS) poses challenges for drug development and the study of cell fate, plasticity, and differentiation, processes linked to tumor grade, cell metastasis, and survival. To pinpoint errors in OS differentiation, we transcriptionally profiled 31,527 cells from a tissue-engineered model that directs MSCs toward adipogenic and osteoblastic fates. Incorporating pre-existing chondrocyte data, we applied trajectory analysis and non-negative matrix factorization (NMF) to generate the first human mesenchymal differentiation atlas. This 'roadmap' served as a reference to delineate the cellular composition of morphologically complex OS tumors and quantify each cell's lineage commitment. Projecting these signatures onto a bulk RNA-seq OS dataset unveiled a correlation between a stem-like transcriptomic phenotype and poorer survival outcomes. Our study takes the critical first step in accurately quantifying OS differentiation and lineage, a prerequisite to better understanding global differentiation bottlenecks that might someday be targeted therapeutically.

Quantifying cancer cell plasticity with gene regulatory networks and single-cell dynamics

Phenotypic plasticity of cancer cells can lead to complex cell state dynamics during tumor progression and acquired resistance. Highly plastic stem-like states may be inherently drug-resistant. Moreover, cell state dynamics in response to therapy allow a tumor to evade treatment. In both scenarios, quantifying plasticity is essential for identifying high-plasticity states or elucidating transition paths between states. Currently, methods to quantify plasticity tend to focus on 1) quantification of quasi-potential based on the underlying gene regulatory network dynamics of the system; or 2) inference of cell potency based on trajectory inference or lineage tracing in single-cell dynamics. Here, we explore both of these approaches and associated computational tools. We then discuss implications of each approach to plasticity metrics, and relevance to cancer treatment strategies.

VarID2 quantifies gene expression noise dynamics and unveils functional heterogeneity of ageing hematopoietic stem cells

Variability of gene expression due to stochasticity of transcription or variation of extrinsic signals, termed biological noise, is a potential driving force of cellular differentiation. Utilizing single-cell RNA-sequencing, we develop VarID2 for the quantification of biological noise at single-cell resolution. VarID2 reveals enhanced nuclear versus cytoplasmic noise, and distinct regulatory modes stratified by correlation between noise, expression, and chromatin accessibility. Noise levels are minimal in murine hematopoietic stem cells (HSCs) and increase during differentiation and ageing. Differential noise identifies myeloid-biased Dlk1+ long-term HSCs in aged mice with enhanced quiescence and self-renewal capacity. VarID2 reveals noise dynamics invisible to conventional single-cell transcriptome analysis.

CD168+ macrophages promote hepatocellular carcinoma tumor stemness and progression through TOP2A/β-catenin/YAP1 axis

Liver cancer stem-like cells (LCSCs) are the main cause of heterogeneity and poor prognosis in hepatocellular carcinoma (HCC). In this study, we aimed to explore the origin of LCSCs and the role of the TOP2A/β-catenin/YAP1 axis in tumor stemness and progression. Using single-cell RNA-seq analysis, we identified TOP2A+CENPF+ LCSCs, which were mainly regulated by CD168+ M2-like macrophages. Furthermore, spatial location analysis and fluorescent staining confirmed that LCSCs were enriched at tumor margins, constituting the spatial heterogeneity of HCC. Mechanistically, TOP2A competitively binds to β-catenin, leading to disassociation of β-catenin from YAP1, promoting HCC stemness and overgrowth. Our study provides valuable insights into the spatial transcriptome heterogeneity of the HCC microenvironment and the critical role of TOP2A/β-catenin/YAP1 axis in HCC stemness and progression.

The Breast Cancer Single-Cell Atlas: Defining cellular heterogeneity within model cell lines and primary tumors to inform disease subtype, stemness, and treatment options

PURPOSE

Breast Cancer (BC) is the most diagnosed cancer in women; however, through significant research, relative survival rates have significantly improved. Despite progress, there remains a gap in our understanding of BC subtypes and personalized treatments. This manuscript characterized cellular heterogeneity in BC cell lines through scRNAseq to resolve variability in subtyping, disease modeling potential, and therapeutic targeting predictions.

METHODS

We generated a Breast Cancer Single-Cell Cell Line Atlas (BSCLA) to help inform future BC research. We sequenced over 36,195 cells composed of 13 cell lines spanning the spectrum of clinical BC subtypes and leveraged publicly available data comprising 39,214 cells from 26 primary tumors.

RESULTS

Unsupervised clustering identified 49 subpopulations within the cell line dataset. We resolve ambiguity in subtype annotation comparing expression of Estrogen Receptor, Progesterone Receptor, and Human Epidermal Growth Factor Receptor 2 genes. Gene correlations with disease subtype highlighted S100A7 and MUCL1 overexpression in HER2 + cells as possible cell motility and localization drivers. We also present genes driving populational drifts to generate novel gene vectors characterizing each subpopulation. A global Cancer Stem Cell (CSC) scoring vector was used to identify stemness potential for subpopulations and model multi-potency. Finally, we overlay the BSCLA dataset with FDA-approved targets to identify to predict the efficacy of subpopulation-specific therapies.

CONCLUSION

The BSCLA defines the heterogeneity within BC cell lines, enhancing our overall understanding of BC cellular diversity to guide future BC research, including model cell line selection, unintended sample source effects, stemness factors between cell lines, and cell type-specific treatment response.

Cell Features Reconstruction from Gene Association Network of Single Cell

Gene expression as an unstable form of cell characterization has been widely used for single-cell analyses. Although there are cell-specific networks (CSN) to explore stable gene associations within a single cell, the amount of information in CSN is huge and there is no method to measure the interaction level between genes. Therefore, this paper presents a two-level approach to reconstructing single-cell features, which transforms the original gene expression feature into the gene ontology feature and gene interaction feature. Specifically, we first squeeze all CSNs into a cell network feature matrix (CNFM) by fusing the global position and neighborhood influence of genes. Next, we propose a computational method of gene gravitation based on CNFM to quantify the extent of gene-gene interaction, and we can construct a gene gravitation network for single cells. Finally, we further design a novel index of gene gravitation entropy to quantitatively evaluate the level of single-cell differentiation. The experiments on eight different scRNA-seq datasets show the effectiveness and broad application prospects of our method.

The Single-Cell Landscape of Intratumoral Heterogeneity and The Immunosuppressive Microenvironment in Liver and Brain Metastases of Breast Cancer

Distant metastasis remains the major cause of morbidity for breast cancer. Individuals with liver or brain metastasis have an extremely poor prognosis and low response rates to anti-PD-1/L1 immune checkpoint therapy compared to those with metastasis at other sites. Therefore, it is urgent to investigate the underlying mechanism of anti-PD-1/L1 resistance and develop more effective immunotherapy strategies for these patients. Using single-cell RNA sequencing, a high-resolution map of the entire tumor ecosystem based on 44 473 cells from breast cancer liver and brain metastases is depicted. Identified by canonical markers and confirmed by multiplex immunofluorescent staining, the metastatic ecosystem features remarkable reprogramming of immunosuppressive cells such as FOXP3+ regulatory T cells, LAMP3+ tolerogenic dendritic cells, CCL18+ M2-like macrophages, RGS5+ cancer-associated fibroblasts, and LGALS1+ microglial cells. In addition, PD-1 and PD-L1/2 are barely expressed in CD8+ T cells and cancer/immune/stromal cells, respectively. Interactions of the immune checkpoint molecules LAG3-LGALS3 and TIGIT-NECTIN2 between CD8+ T cells and cancer/immune/stromal cells are found to play dominant roles in the immune escape. In summary, this study dissects the intratumoral heterogeneity and immunosuppressive microenvironment in liver and brain metastases of breast cancer for the first time, providing insights into the most appropriate immunotherapy strategies for these patients.

In vivo-like Culture of Monophagous Animal Organ using Dietary Components

Animals depend on other species to live, with monophagy being an extreme mode. Monophagous animals depend on their diet not only for nutritients but also for developmental and reproductive controls. Thus, dietary components may be useful in culturing tissues from monophagous animals. We hypothesized that a dedifferentiated tissue from the monophagous silkworm, Bombyx mori, would re-differentiate when cultured in a medium containing an extract of mulberry (Morus alba) leaves, the only food of B. mori. Over 40 fat-body transcriptomes were sequenced, and we concluded that it is possible to establish in vivo-like silkworm tissue cultures using their diet.

W, Chen H, Yu C, Yang L, Fu Y, Li J, Fu S, Sun Z, Fei L, Guo Q, Wang J, Xiao Y, Wang X, Zhang P, Ma L, Ge D, Xu S, Caballero-Pérez J, Cruz-Ramírez A, Zhou Y, Chen M, Fei JF, Han X, Guo G. Construction of the axolotl cell landscape using combinatorial hybridization sequencing at single-cell resolution

The Mexican axolotl (Ambystoma mexicanum) is a well-established tetrapod model for regeneration and developmental studies. Remarkably, neotenic axolotls may undergo metamorphosis, a process that triggers many dramatic changes in diverse organs, accompanied by gradually decline of their regeneration capacity and lifespan. However, the molecular regulation and cellular changes in neotenic and metamorphosed axolotls are still poorly investigated. Here, we develop a single-cell sequencing method based on combinatorial hybridization to generate a tissue-based transcriptomic landscape of the neotenic and metamorphosed axolotls. We perform gene expression profiling of over 1 million single cells across 19 tissues to construct the first adult axolotl cell landscape. Comparison of single-cell transcriptomes between the tissues of neotenic and metamorphosed axolotls reveal the heterogeneity of non-immune parenchymal cells in different tissues and established their regulatory network. Furthermore, we describe dynamic gene expression patterns during limb development in neotenic axolotls. This system-level single-cell analysis of molecular characteristics in neotenic and metamorphosed axolotls, serves as a resource to explore the molecular identity of the axolotl and facilitates better understanding of metamorphosis.

The Mexican axolotl is a well-established tetrapod model for regeneration and development. Here the authors report a scRNA-seq method to profile neotenic, metamorphic and limb development stages, highlighting unique perturbation patterns of cell type-related gene expression throughout metamorphosis.

Decoding the colorectal cancer ecosystem emphasizes the cooperative role of cancer cells, TAMs and CAFsin tumor progression

Background

Single-cell transcription data provided unprecedented molecular information, enabling us to directly encode the ecosystem of colorectal cancer (CRC). Characterization of the diversity of epithelial cells and how they cooperate with tumor microenvironment cells (TME) to endow CRC with aggressive characteristics at single-cell resolution is critical for the understanding of tumor progression mechanism.

Methods

In this study, we comprehensively analyzed the single-cell transcription data, bulk-RNA sequencing data and pathological tissue data. In detail, cellular heterogeneity of TME and epithelial cells were analyzed by unsupervised classification and consensus nonnegative matrix factorization analysis, respectively. Functional status of epithelial clusters was annotated by CancerSEA and its crosstalk with TME cells was investigated using CellPhoneDB and correlation analysis. Findings from single-cell transcription data were further validated in bulk-RNA sequencing data and pathological tissue data.

Results

A distinct cellular composition was observed between tumor and normal tissues, and tumors exhibited immunosuppressive phenotypes. Regarding epithelial cells, we identified one highly invasiveQuery cluster, C4, that correlated closely with tumor-associated macrophages (TAMs) and cancer-associated fibroblasts (CAFs). Further analysis emphasized the TAMs subclass TAM1 and CAFs subclass S5 are closely related with C4.

Conclusions

In summary, our study elaborates on the cellular heterogeneity of CRC, revealing that TAMs and CAFs were critical for crosstalk network epithelial cells and TME cells. This in-depth understanding of cancer cell-TME network provided theoretical basis for the development of new drugs targeting this sophisticated network in CRC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-022-03661-8.

Transcriptome Dedifferentiation Observed in Animal Primary Cultures is Essential to Plant Reprogramming

Tissue culture environment liberate cells from ordinary laws of multi-cellular organisms. This liberation enables cells several behaviors, such as growth, dedifferentiation, acquisition of pluripotency, immortalization and reprogramming. Each phenomenon is relating to each other and hardly to determine. Recently, dedifferentiation of animal cell was quantified as increasing liberality which is information entropy of transcriptome. The increasing liberality induced by tissue culture may reappear in plant cells too. Here we corroborated it. Measuring liberality during reprogramming of plant cells suggested that reprogramming is a combined phenomenon of dedifferentiation and re-differentiation.

Atlas of human dental pulp cells at multiple spatial and temporal levels based on single-cell sequencing analysis

The dental pulp plays a crucial role in the long-term maintenance of tooth function. The progress of endodontic treatment and pulp tissue regeneration engineering has made pulp-regeneration therapy promising in clinical practice. However, the mechanisms of pulp regeneration and the role of dental stem cells in development and regeneration have not been fully elucidated. Bridging the gaps between clinical operation and basic research is urgently needed. With the application of single-cell sequencing technology in dental research, the landscapes of human dental pulp cells have begun being outlined. However, the specific cellular heterogeneity of dental pulp cells, especially that of dental stem cells, at different spatial and temporal levels, is still unclear. In this study, we used single-cell RNA sequencing analysis of pulp samples at four different developmental stages and combined the findings with immunohistochemical staining to explore the development of dental pulp and stem cells. The results revealed temporal changes in the proportion of pulp cells during development. For example, mononuclear phagocytes accounted for a higher proportion in early samples. Odontoblasts identified by DMP1 had a higher expression of ion channel-related and neurodevelopment-related genes. Subpopulations were identified in fibroblasts, odontoblasts, and mesenchymal stem cells. We identified a subclass of odontoblasts that expresses DGKI and RRBP1 present in early developmental samples. A population of earlier mesenchymal stem cells expressed the SEPTIN gene, which may have greater proliferative and differentiation potential. Furthermore, dental pulp stem cells can differentiate into two directions: mineralization and myogenesis. In summary, the specific cellular heterogeneity of dental pulp cells was revealed at different spatial and temporal levels. These findings may shed light on the mechanism of tooth development. The gene expression profile of developing pulp cells may help to select cells for regenerative engineering and improve the success of dental pulp regeneration.

Phenotypic heterogeneity driven by plasticity of the intermediate EMT state governs disease progression and metastasis in breast cancer

The epithelial-to-mesenchymal transition (EMT) is frequently co-opted by cancer cells to enhance migratory and invasive cell traits. It is a key contributor to heterogeneity, chemoresistance, and metastasis in many carcinoma types, where the intermediate EMT state plays a critical tumor-initiating role. We isolate multiple distinct single-cell clones from the SUM149PT human breast cell line spanning the EMT spectrum having diverse migratory, tumor-initiating, and metastatic qualities, including three unique intermediates. Using a multiomics approach, we identify CBFβ as a key regulator of metastatic ability in the intermediate state. To quantify epithelial-mesenchymal heterogeneity within tumors, we develop an advanced multiplexed immunostaining approach using SUM149-derived orthotopic tumors and find that the EMT state and epithelial-mesenchymal heterogeneity are predictive of overall survival in a cohort of stage III breast cancer. Our model reveals previously unidentified insights into the complex EMT spectrum and its regulatory networks, as well as the contributions of epithelial-mesenchymal plasticity (EMP) in tumor heterogeneity in breast cancer.

Single-Cell RNA Sequencing Unravels Distinct Tumor Microenvironment of Different Components of Lung Adenocarcinoma Featured as Mixed Ground-Glass Opacity

Lung adenocarcinoma featured as mixed ground-glass opacity (mGGO) doubled its volume half of the time in comparison with that featured as pure ground-glass opacity (pGGO). The mechanisms underlying the heterogeneous appearance of mGGO remain elusive. In this study, we macro-dissected the solid (S) components and ground-glass (GG) components of mGGO and performed single-cell sequencing analyses of six paired components from three mGGO patients. A total of 19,391 single-cell profiles were taken into analysis, and the data of each patient were analyzed independently to obtain a common alteration. Cancer cells and macrophages were the dominant cell types in the S and GG components, respectively. Cancer cells in the S components, which showed relatively malignant phenotypes, were likely to originate from both the GG and S components and monitor the surrounding tumor microenvironment (TME) through an intricate cell interaction network. SPP1^hi macrophages were enriched in the S components and showed increased activity of chemoattraction, while macrophages in the GG components displayed an active antimicrobial process with a higher stress-induced state. In addition, the CD47–SIRPA axis was demonstrated to be critical in the maintenance of the GG components. Taken together, our study unraveled the alterations of cell components and transcriptomic features between different components in mGGOs.

Integrated single-cell analyses decode the developmental landscape of the human fetal spine

The spine has essential roles in supporting body weight, and passaging the neural elements between the body and the brain. In this study, we used integrated single-cell RNA sequencing and single-cell transposase-accessible chromatin sequencing analyses to reveal the cellular heterogeneity, lineage, and transcriptional regulatory network of the developing human spine. We found that EPYC + HAPLN1+ fibroblasts with stem cell characteristics could differentiate into chondrocytes by highly expressing the chondrogenic markers SOX9 and MATN4. Neurons could originate from neuroendocrine cells, and MEIS2 may be an essential transcription factor that promotes spinal neural progenitor cells to selectively differentiate into neurons during early gestation. Furthermore, the interaction of NRP2_SEMA3C and CD74_APP between macrophages and neurons may be essential for spinal cord development. Our integrated map provides a blueprint for understanding human spine development in the early and midgestational stages at single-cell resolution and offers a tool for investigating related diseases.

•

scRNA-seq and scATAC-seq analyses reveal the developmental landscape of the fetal spine

•

Chondrocytes may originate from EPYC + HAPLN1+ fibroblasts with stem cell characteristics

•

Neurons may originate from neuroendocrine cells with regulation by MEIS2

Optoplasmonic Modulation of Cell Metabolic State Promotes Rapid Cell Differentiation

Cell differentiation plays a vital role in mediating organ formation and tissue repair and regeneration. Although rapid and effective methods to stimulate cell differentiation for clinical purposes are highly desired, it remains a great challenge in the medical fields. Herein, a highly effective and conceptual optical method was developed based on a plasmonic chip platform (made of 2D AuNPs nanomembranes). through effective light-augmented plasmonic regulation of cellular bioenergetics (CBE) and an entropy effect at bionano interfaces, to promote rapid cell differentiation. Compared with traditional methods, the developed optoplasmonic method greatly shortens cell differentiation time from usually more than 10 days to only about 3 days. Upon the optoplasmonic treatment of cells, the conformational and vibration entropy changes of cell membranes were clearly revealed through theoretical simulation and fingerprint spectra of cell membranes. Meanwhile, during the treatment process, bioenergetics levels of cells were elevated with increasing mitochondrial membrane potential (Δψ_m), which accelerates cell differentiation and proliferation. The developed optoplasmonic method is highly efficient and easy to implement, provides a new perspective and avenue for cell differentiation and proliferation, and has potential application prospects in accelerating tissue repair and regeneration.

Distinct Epithelial Cell Profiles in Normal Versus Induced-Congenital Diaphragmatic Hernia Fetal Lungs

BACKGROUND

Recent studies identified a great diversity of cell types in precise number and position to create the architectural features of the lung that ventilation and respiration at birth depend on. With damaged respiratory function at birth, congenital diaphragmatic hernia (CDH) is one of the more severe causes of fetal lung hypoplasia with unspecified cellular dynamics.

OBJECTIVES

 To characterize the epithelial cell tissue in hypoplastic lungs, a careful analysis regarding pulmonary morphology and epithelial cell profile was conducted from pseudoglandular-to-saccular phases in normal versus nitrofen-induced CDH rat lungs.

DESIGN

Our analysis comprises three experimental groups, control, nitrofen (NF) and CDH, in which the relative expression levels (western blot) by group and developmental stage were analyzed in whole lung. Spatiotemporal distribution (immunohistochemistry) was revealed by pulmonary structure during normal and hypoplastic fetal lung development. Surfactant protein-C (SP-C), calcitonin gene-related peptide (CGRP), clara cell secretory protein (CCSP), and forkhead box J1 (FOXJ1) were the used molecular markers for alveolar epithelial cell type 2 (AEC2), pulmonary neuroendocrine, clara, and ciliated cell profiles, respectively.

RESULTS

Generally, we identified an aberrant expression of SP-C, CGRP, CCSP, and FOXJ1 in nitrofen-exposed lungs. For instance, the overexpression of FOXJ1 and CGRP in primordia of bronchiole defined the pseudoglandular stage in CDH lungs, whereas the increased expression of CGRP in bronchi; FOXJ1 and CGRP in terminal bronchiole; and SP-C in BADJ classified the canalicular and saccular stages in hypoplastic lungs. We also described higher expression levels in NF than CDH or control groups for both FOXJ1 in bronchi, terminal bronchiole and BADJ at canalicular stage, and SP-C in bronchi and terminal bronchiole at canalicular and saccular stages. Finally, we report an unexpected expression of FOXJ1 in BADJ at canalicular and saccular stages, whereas the multi cilia observed in bronchi were notably absent at embryonic day 21.5 in induced-CDH lungs.

CONCLUSION

The recognized alterations in the epithelial cell profile contribute to a better understanding of neonatal respiratory insufficiency in induced-CDH lungs and indicate a problem in the epithelial cell differentiation in hypoplastic lungs.

Computational Methods to Identify Cell-Fate Determinants, Identity Transcription Factors, and Niche-Induced Signaling Pathways for Stem Cell Research

The large-scale development of high-throughput sequencing technologies has not only allowed the generation of reliable omics data related to various regulatory layers but also the development of novel computational models in the field of stem cell research. These computational approaches have enabled the disentangling of a complex interplay between these interrelated layers of regulation by interpreting large quantities of biomedical data in a systematic way. In the context of stem cell research, network modeling of complex gene-gene interactions has been successfully used for understanding the mechanisms underlying stem cell differentiation and cellular conversion. Notably, it has proven helpful for predicting cell-fate determinants and signaling molecules controlling such processes. This chapter will provide an overview of various computational approaches that rely on single-cell and/or bulk RNA sequencing data for elucidating the molecular underpinnings of cell subpopulation identities, lineage specification, and the process of cell-fate decisions. Furthermore, we discuss how these computational methods provide the right framework for computational modeling of biological systems in order to address long-standing challenges in the stem cell field by guiding experimental efforts in stem cell research and regenerative medicine.

Cellular origins and lineage relationships of the intestinal epithelium

Knowledge of the development and hierarchical organization of tissues is key to understanding how they are perturbed in injury and disease, as well as how they may be therapeutically manipulated to restore homeostasis. The rapidly regenerating intestinal epithelium harbors diverse cell types and their lineage relationships have been studied using numerous approaches, from classical label-retaining and genetic lineage tracing methods to novel transcriptome-based annotations. Here, we describe the developmental trajectories that dictate differentiation and lineage specification in the intestinal epithelium. We focus on the most recent single-cell RNA-sequencing (scRNA-seq)-based strategies for understanding intestinal epithelial cell lineage relationships, underscoring how they have refined our view of the development of this tissue and highlighting their advantages and limitations. We emphasize how these technologies have been applied to understand the dynamics of intestinal epithelial cells in homeostatic and injury-induced regeneration models.

Transcriptomic entropy benchmarks stem cell-derived cardiomyocyte maturation against endogenous tissue at single cell level

The immaturity of pluripotent stem cell (PSC)-derived tissues has emerged as a universal problem for their biomedical applications. While efforts have been made to generate adult-like cells from PSCs, direct benchmarking of PSC-derived tissues against in vivo development has not been established. Thus, maturation status is often assessed on an ad-hoc basis. Single cell RNA-sequencing (scRNA-seq) offers a promising solution, though cross-study comparison is limited by dataset-specific batch effects. Here, we developed a novel approach to quantify PSC-derived cardiomyocyte (CM) maturation through transcriptomic entropy. Transcriptomic entropy is robust across datasets regardless of differences in isolation protocols, library preparation, and other potential batch effects. With this new model, we analyzed over 45 scRNA-seq datasets and over 52,000 CMs, and established a cross-study, cross-species CM maturation reference. This reference enabled us to directly compare PSC-CMs with the in vivo developmental trajectory and thereby to quantify PSC-CM maturation status. We further found that our entropy-based approach can be used for other cell types, including pancreatic beta cells and hepatocytes. Our study presents a biologically relevant and interpretable metric for quantifying PSC-derived tissue maturation, and is extensible to numerous tissue engineering contexts.

Ultra-fast scalable estimation of single-cell differentiation potency from scRNA-Seq data

Motivation

An important task in the analysis of single-cell RNA-Seq data is the estimation of differentiation potency, as this can help identify stem-or-multipotent cells in non-temporal studies or in tissues where differentiation hierarchies are not well established. A key challenge in the estimation of single-cell potency is the need for a fast and accurate algorithm, scalable to large scRNA-Seq studies profiling millions of cells.

Results

Here, we present a single-cell potency measure, called Correlation of Connectome and Transcriptome (CCAT), which can return accurate single-cell potency estimates of a million cells in minutes, a 100-fold improvement over current state-of-the-art methods. We benchmark CCAT against 8 other single-cell potency models and across 28 scRNA-Seq studies, encompassing over 2 million cells, demonstrating comparable accuracy than the current state-of-the-art, at a significantly reduced computational cost, and with increased robustness to dropouts.

Availability and implementation

CCAT is part of the SCENT R-package, freely available from https://github.com/aet21/SCENT.

Supplementary information

Supplementary data are available at Bioinformatics online.

Inferring transcriptomic cell states and transitions only from time series transcriptome data

Cellular stages of biological processes have been characterized using fluorescence-activated cell sorting and genetic perturbations, charting a limited landscape of cellular states. Time series transcriptome data can help define new cellular states at the molecular level since the analysis of transcriptional changes can provide information on cell states and transitions. However, existing methods for inferring cell states from transcriptome data use additional information such as prior knowledge on cell types or cell-type-specific markers to reduce the complexity of data. In this study, we present a novel time series clustering framework to infer TRAnscriptomic Cellular States (TRACS) only from time series transcriptome data by integrating Gaussian process regression, shape-based distance, and ranked pairs algorithm in a single computational framework. TRACS determines patterns that correspond to hidden cellular states by clustering gene expression data. TRACS was used to analyse single-cell and bulk RNA sequencing data and successfully generated cluster networks that reflected the characteristics of key stages of biological processes. Thus, TRACS has a potential to help reveal unknown cellular states and transitions at the molecular level using only time series transcriptome data. TRACS is implemented in Python and available at http://github.com/BML-cbnu/TRACS/ .

What are the applications of single-cell RNA sequencing in cancer research: a systematic review

Single-cell RNA sequencing (scRNA-seq) is a tool for studying gene expression at the single-cell level that has been widely used due to its unprecedented high resolution. In the present review, we outline the preparation process and sequencing platforms for the scRNA-seq analysis of solid tumor specimens and discuss the main steps and methods used during data analysis, including quality control, batch-effect correction, normalization, cell cycle phase assignment, clustering, cell trajectory and pseudo-time reconstruction, differential expression analysis and gene set enrichment analysis, as well as gene regulatory network inference. Traditional bulk RNA sequencing does not address the heterogeneity within and between tumors, and since the development of the first scRNA-seq technique, this approach has been widely used in cancer research to better understand cancer cell biology and pathogenetic mechanisms. ScRNA-seq has been of great significance for the development of targeted therapy and immunotherapy. In the second part of this review, we focus on the application of scRNA-seq in solid tumors, and summarize the findings and achievements in tumor research afforded by its use. ScRNA-seq holds promise for improving our understanding of the molecular characteristics of cancer, and potentially contributing to improved diagnosis, prognosis, and therapeutics.

The Cellular and Physiological Basis for Lung Repair and Regeneration: Past, Present, and Future

The respiratory system, which includes the trachea, airways, and distal alveoli, is a complex multi-cellular organ that intimately links with the cardiovascular system to accomplish gas exchange. In this review and as members of the NIH/NHLBI-supported Progenitor Cell Translational Consortium, we discuss key aspects of lung repair and regeneration. We focus on the cellular compositions within functional niches, cell-cell signaling in homeostatic health, the responses to injury, and new methods to study lung repair and regeneration. We also provide future directions for an improved understanding of the cell biology of the respiratory system, as well as new therapeutic avenues.

Mapping the regulatory landscape of auditory hair cells from single-cell multi-omics data

Auditory hair cells transduce sound to the brain and in mammals these cells reside together with supporting cells in the sensory epithelium of the cochlea, called the organ of Corti. To establish the organ's delicate function during development and differentiation, spatiotemporal gene expression is strictly controlled by chromatin accessibility and cell type-specific transcription factors, jointly representing the regulatory landscape. Bulk-sequencing technology and cellular heterogeneity obscured investigations on the interplay between transcription factors and chromatin accessibility in inner ear development. To study the formation of the regulatory landscape in hair cells, we collected single-cell chromatin accessibility profiles accompanied by single-cell RNA data from genetically labeled murine hair cells and supporting cells after birth. Using an integrative approach, we predicted cell type-specific activating and repressing functions of developmental transcription factors. Furthermore, by integrating gene expression and chromatin accessibility datasets, we reconstructed gene regulatory networks. Then, using a comparative approach, 20 hair cell-specific activators and repressors, including putative downstream target genes, were identified. Clustering of target genes resolved groups of related transcription factors and was utilized to infer their developmental functions. Finally, the heterogeneity in the single-cell data allowed us to spatially reconstruct transcriptional as well as chromatin accessibility trajectories, indicating that gradual changes in the chromatin accessibility landscape were lagging behind the transcriptional identity of hair cells along the organ's longitudinal axis. Overall, this study provides a strategy to spatially reconstruct the formation of a lineage specific regulatory landscape using a single-cell multi-omics approach.

NetMix: A Network-Structured Mixture Model for Reduced-Bias Estimation of Altered Subnetworks

A classic problem in computational biology is the identification of altered subnetworks: subnetworks of an interaction network that contain genes/proteins that are differentially expressed, highly mutated, or otherwise aberrant compared with other genes/proteins. Numerous methods have been developed to solve this problem under various assumptions, but the statistical properties of these methods are often unknown. For example, some widely used methods are reported to output very large subnetworks that are difficult to interpret biologically. In this work, we formulate the identification of altered subnetworks as the problem of estimating the parameters of a class of probability distributions that we call the Altered Subset Distribution (ASD). We derive a connection between a popular method, jActiveModules, and the maximum likelihood estimator (MLE) of the ASD. We show that the MLE is statistically biased, explaining the large subnetworks output by jActiveModules. Based on these insights, we introduce NetMix, an algorithm that uses Gaussian mixture models to obtain less biased estimates of the parameters of the ASD. We demonstrate that NetMix outperforms existing methods in identifying altered subnetworks on both simulated and real data, including the identification of differentially expressed genes from both microarray and RNA-seq experiments and the identification of cancer driver genes in somatic mutation data.

Inference and multiscale model of epithelial-to-mesenchymal transition via single-cell transcriptomic data

Rapid growth of single-cell transcriptomic data provides unprecedented opportunities for close scrutinizing of dynamical cellular processes. Through investigating epithelial-to-mesenchymal transition (EMT), we develop an integrative tool that combines unsupervised learning of single-cell transcriptomic data and multiscale mathematical modeling to analyze transitions during cell fate decision. Our approach allows identification of individual cells making transition between all cell states, and inference of genes that drive transitions. Multiscale extractions of single-cell scale outputs naturally reveal intermediate cell states (ICS) and ICS-regulated transition trajectories, producing emergent population-scale models to be explored for design principles. Testing on the newly designed single-cell gene regulatory network model and applying to twelve published single-cell EMT datasets in cancer and embryogenesis, we uncover the roles of ICS on adaptation, noise attenuation, and transition efficiency in EMT, and reveal their trade-off relations. Overall, our unsupervised learning method is applicable to general single-cell transcriptomic datasets, and our integrative approach at single-cell resolution may be adopted for other cell fate transition systems beyond EMT.

Identification of Differential Gene Groups From Single-Cell Transcriptomes Using Network Entropy

A complex tissue contains a variety of cells with distinct molecular signatures. Single-cell RNA sequencing has characterized the transcriptomes of different cell types and enables researchers to discover the underlying mechanisms of cellular heterogeneity. A critical task in single-cell transcriptome studies is to uncover transcriptional differences among specific cell types. However, the intercellular transcriptional variation is usually confounded with high level of technical noise, which masks the important biological signals. Here, we propose a new computational method DiffGE for differential analysis, adopting network entropy to measure the expression dynamics of gene groups among different cell types and to identify the highly differential gene groups. To evaluate the effectiveness of our proposed method, DiffGE is applied to three independent single-cell RNA-seq datasets and to identify the highly dynamic gene groups that exhibit distinctive expression patterns in different cell types. We compare the results of our method with those of three widely applied algorithms. Further, the gene function analysis indicates that these detected differential gene groups are significantly related to cellular regulation processes. The results demonstrate the power of our method in evaluating the transcriptional dynamics and identifying highly differential gene groups among different cell types.

The elephant in the lung: Integrating lineage-tracing, molecular markers, and single cell sequencing data to identify distinct fibroblast populations during lung development and regeneration

During lung development, the mesenchyme and epithelium are dependent on each other for instructive morphogenic cues that direct proliferation, cellular differentiation and organogenesis. Specification of epithelial and mesenchymal cell lineages occurs in parallel, forming cellular subtypes that guide the formation of both transitional developmental structures and the permanent architecture of the adult lung. While epithelial cell types and lineages have been relatively well-defined in recent years, the definition of mesenchymal cell types and lineage relationships has been more challenging. Transgenic mouse lines with permanent and inducible lineage tracers have been instrumental in identifying lineage relationships among epithelial progenitor cells and their differentiation into distinct airway and alveolar epithelial cells. Lineage tracing experiments with reporter mice used to identify fibroblast progenitors and their lineage trajectories have been limited by the number of cell specific genes and the developmental timepoint when the lineage trace was activated. In this review, we discuss major developmental mesenchymal lineages, focusing on time of origin, major cell type, and other lineage derivatives, as well as the transgenic tools used to find and define them. We describe lung fibroblasts using function, location, and molecular markers in order to compare and contrast cells with similar functions. The temporal and cell-type specific expression of fourteen "fibroblast lineage" genes were identified in single-cell RNA-sequencing data from LungMAP in the LGEA database. Using these lineage signature genes as guides, we clustered murine lung fibroblast populations from embryonic day 16.5 to postnatal day 28 (E16.5-PN28) and generated heatmaps to illustrate expression of transcription factors, signaling receptors and ligands in a temporal and population specific manner.

Human Cell Atlas and cell-type authentication for regenerative medicine

In modern biology, the correct identification of cell types is required for the developmental study of tissues and organs and the production of functional cells for cell therapies and disease modeling. For decades, cell types have been defined on the basis of morphological and physiological markers and, more recently, immunological markers and molecular properties. Recent advances in single-cell RNA sequencing have opened new doors for the characterization of cells at the individual and spatiotemporal levels on the basis of their RNA profiles, vastly transforming our understanding of cell types. The objective of this review is to survey the current progress in the field of cell-type identification, starting with the Human Cell Atlas project, which aims to sequence every cell in the human body, to molecular marker databases for individual cell types and other sources that address cell-type identification for regenerative medicine based on cell data guidelines.

Thermal Resonance and Cell Behavior

From a thermodynamic point of view, living cell life is no more than a cyclic process. It starts with the newly separated daughter cells and restarts when the next generations grow as free entities. During this cycle, the cell changes its entropy. In cancer, the growth control is damaged. In this paper, we analyze the role of the volume–area ratio in the cell in relation to the heat exchange between cell and its environment in order to point out its effect on cancer growth. The result holds to a possible control of the cancer growth based on the heat exchanged by the cancer toward its environment and the membrane potential variation, with the consequence of controlling the ions fluxes and the related biochemical reactions. This second law approach could represent a starting point for a possible future support for the anticancer therapies, in order to improve their effectiveness for the untreatable cancers.

The cytoskeletal regulator HEM1 governs B cell development and prevents autoimmunity

The WAVE regulatory complex (WRC) is crucial for assembly of the peripheral branched actin network constituting one of the main drivers of eukaryotic cell migration. Here, we uncover an essential role of the hematopoietic-specific WRC component HEM1 for immune cell development. Germline-encoded HEM1 deficiency underlies an inborn error of immunity with systemic autoimmunity, at cellular level marked by WRC destabilization, reduced filamentous actin, and failure to assemble lamellipodia. Hem1^-/- mice display systemic autoimmunity, phenocopying the human disease. In the absence of Hem1, B cells become deprived of extracellular stimuli necessary to maintain the strength of B cell receptor signaling at a level permissive for survival of non-autoreactive B cells. This shifts the balance of B cell fate choices toward autoreactive B cells and thus autoimmunity.

Glucocorticoid regulates mesenchymal cell differentiation required for perinatal lung morphogenesis and function

While antenatal glucocorticoids are widely used to enhance lung function in preterm infants, cellular and molecular mechanisms by which glucocorticoid receptor (GR) signaling influences lung maturation remain poorly understood. Deletion of the glucocorticoid receptor gene (Nr3c1) from fetal pulmonary mesenchymal cells phenocopied defects caused by global Nr3c1 deletion, while lung epithelial- or endothelial-specific Nr3c1 deletion did not impair lung function at birth. We integrated genome-wide gene expression profiling, ATAC-seq, and single cell RNA-seq data in mice in which GR was deleted or activated to identify the cellular and molecular mechanisms by which glucocorticoids control prenatal lung maturation. GR enhanced differentiation of a newly defined proliferative mesenchymal progenitor cell (PMP) into matrix fibroblasts (MFBs), in part by directly activating extracellular matrix-associated target genes, including Fn1, Col16a4, and Eln and by modulating VEGF, JAK-STAT, and WNT signaling. Loss of mesenchymal GR signaling blocked fibroblast progenitor differentiation into mature MFBs, which in turn increased proliferation of SOX9+ alveolar epithelial progenitor cells and inhibited differentiation of mature alveolar type II (AT₂) and AT₁ cells. GR signaling controls genes required for differentiation of a subset of proliferative mesenchymal progenitors into matrix fibroblasts, in turn, regulating signals controlling AT₂/AT₁ progenitor cell proliferation and differentiation and identifying cells and processes by which glucocorticoid signaling regulates fetal lung maturation.

Application of information theoretical approaches to assess diversity and similarity in single-cell transcriptomics

Single-cell transcriptomics offers a powerful way to reveal the heterogeneity of individual cells. To date, many information theoretical approaches have been proposed to assess diversity and similarity, and characterize the latent heterogeneity in transcriptome data. Diversity implies gene expression variations and can facilitate the identification of signature genes; while, similarity unravels co-expression patterns for cell type clustering. In this review, we summarized 16 measures of information theory used for evaluating diversity and similarity in single-cell transcriptomic data, provide references and shed light on selected theoretical properties when there is a need to select proper measurements in general cases. We further provide an R package assembling discussed approaches to improve the researchers own single-cell transcriptome study. At last, we prospected further applications of diversity and similarity measures in support of depicting heterogeneity in single-cell multi-omics data.

scRCMF: Identification of Cell Subpopulations and Transition States From Single-Cell Transcriptomes

Single cell technologies provide an unprecedented opportunity to explore the heterogeneity in a biological process at the level of single cells. One major challenge in analyzing single cell data is to identify cell subpopulations, stable cell states, and cells in transition between states. To elucidate the transition mechanisms in cell fate dynamics, it is highly desirable to quantitatively characterize cellular states and intermediate states. Here, we present scRCMF, an unsupervised method that identifies stable cell states and transition cells by adopting a nonlinear optimization model that infers the latent substructures from a gene-cell matrix. We incorporate a random coefficient matrix-based regularization into the standard nonnegative matrix decomposition model to improve the reliability and stability of estimating latent substructures. To quantify the transition capability of each cell, we propose two new measures: single-cell transition entropy (scEntropy) and transition probability (scTP). When applied to two simulated and three published scRNA-seq datasets, scRCMF not only successfully captures multiple subpopulations and transition processes in large-scale data, but also identifies transition states and some known marker genes associated with cell state transitions and subpopulations. Furthermore, the quantity scEntropy is found to be significantly higher for transition cells than other cellular states during the global differentiation, and the scTP predicts the "fate decisions" of transition cells within the transition. The present study provides new insights into transition events during differentiation and development.

Predicting cell lineages using autoencoders and optimal transport

Lineage tracing involves the identification of all ancestors and descendants of a given cell, and is an important tool for studying biological processes such as development and disease progression. However, in many settings, controlled time-course experiments are not feasible, for example when working with tissue samples from patients. Here we present ImageAEOT, a computational pipeline based on autoencoders and optimal transport for predicting the lineages of cells using time-labeled datasets from different stages of a cellular process. Given a single-cell image from one of the stages, ImageAEOT generates an artificial lineage of this cell based on the population characteristics of the other stages. These lineages can be used to connect subpopulations of cells through the different stages and identify image-based features and biomarkers underlying the biological process. To validate our method, we apply ImageAEOT to a benchmark task based on nuclear and chromatin images during the activation of fibroblasts by tumor cells in engineered 3D tissues. We further validate ImageAEOT on chromatin images of various breast cancer cell lines and human tissue samples, thereby linking alterations in chromatin condensation patterns to different stages of tumor progression. Our results demonstrate the promise of computational methods based on autoencoding and optimal transport principles for lineage tracing in settings where existing experimental strategies cannot be used.

Many key biological processes, such as development and disease progression, require analyzing lineages of cells backwards as well as forwards in time. However, current single-cell experiments tend to be destructive to cells, so that a single lineage can only be measured at one point in time. In this work, we introduce a computational framework for predicting the lineage of cells from a single snapshot in time based on measurements of other cells at other time points. The method generates these lineages by computing the most plausible path for a population of cells to transition from one time point to the next, assuming that a cell is more likely to transition to similar cells compared to dissimilar cells. We validate our computational method on imaging data of fibroblasts and cancer cells, though our method could also be applied to other modalities of single-cell data such as genomics and transcriptomics as well as multi-modal single-cell datasets.

Single-cell transcriptional diversity is a hallmark of developmental potential

Single-cell RNA sequencing (scRNA-seq) is a powerful approach for reconstructing cellular differentiation trajectories. However, inferring both the state and direction of differentiation is challenging. Here, we demonstrate a simple, yet robust, determinant of developmental potential-the number of expressed genes per cell-and leverage this measure of transcriptional diversity to develop a computational framework (CytoTRACE) for predicting differentiation states from scRNA-seq data. When applied to diverse tissue types and organisms, CytoTRACE outperformed previous methods and nearly 19,000 annotated gene sets for resolving 52 experimentally determined developmental trajectories. Additionally, it facilitated the identification of quiescent stem cells and revealed genes that contribute to breast tumorigenesis. This study thus establishes a key RNA-based feature of developmental potential and a platform for delineation of cellular hierarchies.

Construction of a human cell landscape at single-cell level

Single-cell analysis is a valuable tool for dissecting cellular heterogeneity in complex systems¹. However, a comprehensive single-cell atlas has not been achieved for humans. Here we use single-cell mRNA sequencing to determine the cell-type composition of all major human organs and construct a scheme for the human cell landscape (HCL). We have uncovered a single-cell hierarchy for many tissues that have not been well characterized. We established a 'single-cell HCL analysis' pipeline that helps to define human cell identity. Finally, we performed a single-cell comparative analysis of landscapes from human and mouse to identify conserved genetic networks. We found that stem and progenitor cells exhibit strong transcriptomic stochasticity, whereas differentiated cells are more distinct. Our results provide a useful resource for the study of human biology.

Deciphering Cell Fate Decision by Integrated Single-Cell Sequencing Analysis

Cellular differentiation is a common underlying feature of all multicellular organisms through which naïve cells progressively become fate restricted and develop into mature cells with specialized functions. A comprehensive understanding of the regulatory mechanisms of cell fate choices during de- velopment, regeneration, homeostasis, and disease is a central goal of mod- ern biology. Ongoing rapid advances in single-cell biology are enabling the exploration of cell fate specification at unprecedented resolution. Here, we review single-cell RNA sequencing and sequencing of other modalities as methods to elucidate the molecular underpinnings of lineage specification. We specifically discuss how the computational tools available to reconstruct lineage trajectories, quantify cell fate bias, and perform dimensionality re- duction for data visualization are providing new mechanistic insights into the process of cell fate decision. Studying cellular differentiation using single- cell genomic tools is paving the way for a detailed understanding of cellular behavior in health and disease.

Using single-cell RNA sequencing to unravel cell lineage relationships in the respiratory tract

The respiratory tract is lined by a pseudo-stratified epithelium from the nose to terminal bronchioles. This first line of defense of the lung against external stress includes five main cell types: basal, suprabasal, club, goblet and multiciliated cells, as well as rare cells such as ionocytes, neuroendocrine and tuft/brush cells. At homeostasis, this epithelium self-renews at low rate but is able of fast regeneration upon damage. Airway epithelial cell lineages during regeneration have been investigated in the mouse by genetic labeling, mainly after injuring the epithelium with noxious agents. From these approaches, basal cells have been identified as progenitors of club, goblet and multiciliated cells, but also of ionocytes and neuroendocrine cells. Single-cell RNA sequencing, coupled to lineage inference algorithms, has independently allowed the establishment of comprehensive pictures of cell lineage relationships in both mouse and human. In line with genetic tracing experiments in mouse trachea, studies using single-cell RNA sequencing (RNAseq) have shown that basal cells first differentiate into club cells, which in turn mature into goblet cells or differentiate into multiciliated cells. In the human airway epithelium, single-cell RNAseq has identified novel intermediate populations such as deuterosomal cells, ‘hybrid’ mucous-multiciliated cells and progenitors of rare cells. Novel differentiation dynamics, such as a transition from goblet to multiciliated cells have also been discovered. The future of cell lineage relationships in the respiratory tract now resides in the combination of genetic labeling approaches with single-cell RNAseq to establish, in a definitive manner, the hallmarks of cellular lineages in normal and pathological situations.

Cell lineage and communication network inference via optimization for single-cell transcriptomics

The use of single-cell transcriptomics has become a major approach to delineate cell subpopulations and the transitions between them. While various computational tools using different mathematical methods have been developed to infer clusters, marker genes, and cell lineage, none yet integrate these within a mathematical framework to perform multiple tasks coherently. Such coherence is critical for the inference of cell–cell communication, a major remaining challenge. Here, we present similarity matrix-based optimization for single-cell data analysis (SoptSC), in which unsupervised clustering, pseudotemporal ordering, lineage inference, and marker gene identification are inferred via a structured cell-to-cell similarity matrix. SoptSC then predicts cell–cell communication networks, enabling reconstruction of complex cell lineages that include feedback or feedforward interactions. Application of SoptSC to early embryonic development, epidermal regeneration, and hematopoiesis demonstrates robust identification of subpopulations, lineage relationships, and pseudotime, and prediction of pathway-specific cell communication patterns regulating processes of development and differentiation.