Single-cell transcriptional diversity is a hallmark of developmental potential

Single-cell and long-read sequencing to enhance modelling of splicing and cell-fate determination

Single-cell sequencing technologies have revolutionised the life sciences and biomedical research. Single-cell sequencing provides high-resolution data on cell heterogeneity, allowing high-fidelity cell type identification, and lineage tracking. Computational algorithms and mathematical models have been developed to make sense of the data, compensate for errors and simulate the biological processes, which has led to breakthroughs in our understanding of cell differentiation, cell-fate determination and tissue cell composition. The development of long-read (a.k.a. third-generation) sequencing technologies has produced powerful tools for investigating alternative splicing, isoform expression (at the RNA level), genome assembly and the detection of complex structural variants (at the DNA level). In this review, we provide an overview of the recent advancements in single-cell and long-read sequencing technologies, with a particular focus on the computational algorithms that help in correcting, analysing, and interpreting the resulting data. Additionally, we review some mathematical models that use single-cell and long-read sequencing data to study cell-fate determination and alternative splicing, respectively. Moreover, we highlight the emerging opportunities in modelling cell-fate determination that result from the combination of single-cell and long-read sequencing technologies.

Single-cell dissection of human hematopoietic reconstitution after allogeneic hematopoietic stem cell transplantation

Hematopoietic stem cell transplantation is an effective regenerative therapy for many malignant, inherited, or autoimmune diseases. However, our understanding of reconstituted hematopoiesis in transplant patients remains limited. Here, we uncover the reconstitution dynamics of human allogeneic hematopoietic stem and progenitor cells (HSPCs) at single-cell resolution after transplantation. Transplanted HSPCs underwent rapid and measurable changes during the first 30 days after transplantation, characterized by a strong proliferative response on the first day. Transcriptomic analysis of HSPCs enabled us to observe that immunoregulatory neutrophil progenitors expressing high levels of the S100A gene family were enriched in granulocyte colony-stimulating factor-mobilized peripheral blood stem cells. Transplant recipients who developed acute graft-versus-host disease (aGVHD) infused fewer S100A^high immunoregulatory neutrophil progenitors, immunophenotyped as Lin^-CD34⁺CD66b⁺CD177⁺, than those who did not develop aGVHD. Therefore, our study provides insights into the regenerative process of transplanted HSPCs in human patients and identifies a potential criterion for identifying patients at high risk for developing aGVHD early after transplant.

Bipotent transitional liver progenitor cells contribute to liver regeneration

Following severe liver injury, when hepatocyte-mediated regeneration is impaired, biliary epithelial cells (BECs) can transdifferentiate into functional hepatocytes. However, the subset of BECs with such facultative tissue stem cell potential, as well as the mechanisms enabling transdifferentiation, remains elusive. Here we identify a transitional liver progenitor cell (TLPC), which originates from BECs and differentiates into hepatocytes during regeneration from severe liver injury. By applying a dual genetic lineage tracing approach, we specifically labeled TLPCs and found that they are bipotent, as they either differentiate into hepatocytes or re-adopt BEC fate. Mechanistically, Notch and Wnt/β-catenin signaling orchestrate BEC-to-TLPC and TLPC-to-hepatocyte conversions, respectively. Together, our study provides functional and mechanistic insights into transdifferentiation-assisted liver regeneration.

Phase-specific signatures of wound fibroblasts and matrix patterns define cancer-associated fibroblast subtypes

Healing wounds and cancers present remarkable cellular and molecular parallels, but the specific roles of the healing phases are largely unknown. We developed a bioinformatics pipeline to identify genes and pathways that define distinct phases across the time-course of healing. Their comparison to cancer transcriptomes revealed that a resolution phase wound signature is associated with increased severity in skin cancer and enriches for extracellular matrix-related pathways. Comparisons of transcriptomes of early- and late-phase wound fibroblasts vs skin cancer-associated fibroblasts (CAFs) identified an "early wound" CAF subtype, which localizes to the inner tumor stroma and expresses collagen-related genes that are controlled by the RUNX2 transcription factor. A "late wound" CAF subtype localizes to the outer tumor stroma and expresses elastin-related genes. Matrix imaging of primary melanoma tissue microarrays validated these matrix signatures and identified collagen- vs elastin-rich niches within the tumor microenvironment, whose spatial organization predicts survival and recurrence. These results identify wound-regulated genes and matrix patterns with prognostic potential in skin cancer.

Dissecting insect cell heterogeneity during influenza VLP production using single-cell transcriptomics

The insect cell-baculovirus expression vector system (IC-BEVS) has been widely used to produce recombinant protein at high titers, including complex virus-like particles (VPLs). However, cell-to-cell variability upon infection is yet one of the least understood phenomena in virology, and little is known about its impact on production of therapeutic proteins. This study aimed at dissecting insect cell population heterogeneity during production of influenza VLPs in IC-BEVS using single-cell RNA-seq (scRNA-seq). High Five cell population was shown to be heterogeneous even before infection, with cell cycle being one of the factors contributing for this variation. In addition, infected insect cells were clustered according to the timing and level of baculovirus genes expression, with each cluster reporting similar influenza VLPs transgenes (i.e., hemagglutinin and M1) transcript counts. Trajectory analysis enabled to track infection progression throughout pseudotime. Specific pathways such as translation machinery, protein folding, sorting and degradation, endocytosis and energy metabolism were identified as being those which vary the most during insect cell infection and production of Influenza VLPs. Overall, this study lays the ground for the application of scRNA-seq in IC-BEVS processes to isolate relevant biological mechanisms during recombinant protein expression towards its further optimization.

The Heterogeneous Cellular States of Glioblastoma Stem Cells Revealed by Single Cell Analysis

Glioblastoma stem cells (GSCs) contributed to the progression, treatment resistance, and relapse of glioblastoma (GBM). However, current researches on GSCs were performed usually outside the human tumor microenvironment, ignoring the importance of the cellular states of primary GSCs. In this study, we leveraged single-cell transcriptome sequencing data of 6 independent GBM cohorts from public databases, and combined lineage and stemness features to identify primary GSCs. We dissected the cell states of GSCs and correlated them with the clinical outcomes of patients. As a result, we constructed a cellular hierarchy where GSCs resided at the center. In addition, we identified and characterized 2 different and recurrent GSCs subpopulations: proliferative GSCs (pGSCs) and quiescent GSCs (qGSCs). The pGSCs showed high cell cycle activity, indicating rapid cell division, while qGSCs showed a quiescent state. Then we traced the processes of tumor development by pseudo-time analysis and tumor phylogeny, and found that GSCs accumulated throughout the whole tumor development period. During the process, pGSCs mainly contributed to the early stage and qGSCs were enriched in the later stage. Finally, we constructed an 8-gene prognostic signature reflecting pGSCs activity and found that patients whose tumors were enriched for the pGSC signature had poor clinical outcomes. Our study highlights the primary GSCs heterogeneity and its correlation to tumor development and clinical outcomes, providing the potential targets for GBM treatment.

Diagnostic classification of childhood cancer using multiscale transcriptomics

The causes of pediatric cancers' distinctiveness compared to adult-onset tumors of the same type are not completely clear and not fully explained by their genomes. In this study, we used an optimized multilevel RNA clustering approach to derive molecular definitions for most childhood cancers. Applying this method to 13,313 transcriptomes, we constructed a pediatric cancer atlas to explore age-associated changes. Tumor entities were sometimes unexpectedly grouped due to common lineages, drivers or stemness profiles. Some established entities were divided into subgroups that predicted outcome better than current diagnostic approaches. These definitions account for inter-tumoral and intra-tumoral heterogeneity and have the potential of enabling reproducible, quantifiable diagnostics. As a whole, childhood tumors had more transcriptional diversity than adult tumors, maintaining greater expression flexibility. To apply these insights, we designed an ensemble convolutional neural network classifier. We show that this tool was able to match or clarify the diagnosis for 85% of childhood tumors in a prospective cohort. If further validated, this framework could be extended to derive molecular definitions for all cancer types.

Towards AI-driven longevity research: An overview

While in the past technology has mostly been utilized to store information about the structural configuration of proteins and molecules for research and medical purposes, Artificial Intelligence is nowadays able to learn from the existing data how to predict and model properties and interactions, revealing important knowledge about complex biological processes, such as aging. Modern technologies, moreover, can rely on a broader set of information, including those derived from the next-generation sequencing (e.g., proteomics, lipidomics, and other omics), to understand the interactions between human body and the external environment. This is especially relevant as external factors have been shown to have a key role in aging. As the field of computational systems biology keeps improving and new biomarkers of aging are being developed, artificial intelligence promises to become a major ally of aging research.

Identification of a physiologic vasculogenic fibroblast state to achieve tissue repair

Tissue injury to skin diminishes miR-200b in dermal fibroblasts. Fibroblasts are widely reported to directly reprogram into endothelial-like cells and we hypothesized that miR-200b inhibition may cause such changes. We transfected human dermal fibroblasts with anti-miR-200b oligonucleotide, then using single cell RNA sequencing, identified emergence of a vasculogenic subset with a distinct fibroblast transcriptome and demonstrated blood vessel forming function in vivo. Anti-miR-200b delivery to murine injury sites likewise enhanced tissue perfusion, wound closure, and vasculogenic fibroblast contribution to perfused vessels in a FLI1 dependent manner. Vasculogenic fibroblast subset emergence was blunted in delayed healing wounds of diabetic animals but, topical tissue nanotransfection of a single anti-miR-200b oligonucleotide was sufficient to restore FLI1 expression, vasculogenic fibroblast emergence, tissue perfusion, and wound healing. Augmenting a physiologic tissue injury adaptive response mechanism that produces a vasculogenic fibroblast state change opens new avenues for therapeutic tissue vascularization of ischemic wounds.

Bulk and single-cell characterisation of the immune heterogeneity of atherosclerosis identifies novel targets for immunotherapy

BACKGROUND

Immune cells that infiltrate lesions are important for atherosclerosis progression and immunotherapies. This study was aimed at gaining important new insights into the heterogeneity of these cells by integrating the sequencing results of multiple samples and using an enhanced single-cell sequencing workflow to overcome the limitations of a single study.

RESULTS

Integrative analyses identified 28 distinct subpopulations based on gene expression profiles. Further analysis demonstrated that these cells manifested high heterogeneity at the levels of tissue preferences, genetic perturbations, functional variations, immune dynamics, transcriptional regulators, metabolic changes, and communication patterns. Of the T cells, interferon-induced CD8⁺ T cells were involved in the progression of atherosclerosis. In contrast, proinflammatory CD4⁺ CD28^null T cells predicted a poor outcome in atherosclerosis. Notably, we identified two subpopulations of foamy macrophages that exhibit contrasting phenotypes. Among them, TREM2^- SPP1⁺ foamy macrophages were preferentially distributed in the hypoxic core of plaques. These glycolytic metabolism-enriched cells, with impaired cholesterol metabolism and robust pro-angiogenic capacity, were phenotypically regulated by CSF1 secreted by co-localised mast cells. Moreover, combined with deconvolution of the bulk datasets, we revealed that these dysfunctional cells had a higher proportion of ruptured and haemorrhagic lesions and were significantly associated with poor atherosclerosis prognoses.

CONCLUSIONS

We systematically explored atherosclerotic immune heterogeneity and identified cell populations underlying atherosclerosis progression and poor prognosis, which may be valuable for developing new and precise immunotherapies.

The early neutrophil-committed progenitors aberrantly differentiate into immunoregulatory monocytes during emergency myelopoiesis

Inflammatory stimuli cause a state of emergency myelopoiesis leading to neutrophil-like monocyte expansion. However, their function, the committed precursors, or growth factors remain elusive. In this study we find that Ym1⁺Ly6C^hi monocytes, an immunoregulatory entity of neutrophil-like monocytes, arise from progenitors of neutrophil 1 (proNeu1). Granulocyte-colony stimulating factor (G-CSF) favors the production of neutrophil-like monocytes through previously unknown CD81⁺CX3CR1^lo monocyte precursors. GFI1 promotes the differentiation of proNeu2 from proNeu1 at the cost of producing neutrophil-like monocytes. The human counterpart of neutrophil-like monocytes that also expands in response to G-CSF is found in CD14⁺CD16^- monocyte fraction. The human neutrophil-like monocytes are discriminated from CD14⁺CD16^- classical monocytes by CXCR1 expression and the capacity to suppress T cell proliferation. Collectively, our findings suggest that the aberrant expansion of neutrophil-like monocytes under inflammatory conditions is a process conserved between mouse and human, which may be beneficial for the resolution of inflammation.

A single-cell transcriptional atlas reveals resident progenitor cell niche functions in TMJ disc development and injury

The biological characteristics of the temporomandibular joint disc involve complex cellular network in cell identity and extracellular matrix composition to modulate jaw function. The lack of a detailed characterization of the network severely limits the development of targeted therapies for temporomandibular joint-related diseases. Here we profiled single-cell transcriptomes of disc cells from mice at different postnatal stages, finding that the fibroblast population could be divided into chondrogenic and non-chondrogenic clusters. We also find that the resident mural cell population is the source of disc progenitors, characterized by ubiquitously active expression of the NOTCH3 and THY1 pathways. Lineage tracing reveals that Myh11⁺ mural cells coordinate angiogenesis during disc injury but lost their progenitor characteristics and ultimately become Sfrp2⁺ non-chondrogenic fibroblasts instead of Chad⁺ chondrogenic fibroblasts. Overall, we reveal multiple insights into the coordinated development of disc cells and are the first to describe the resident mural cell progenitor during disc injury.

Epstein-Barr virus DNA seropositivity links distinct tumoral heterogeneity and immune landscape in nasopharyngeal carcinoma

Background

Epstein-Barr virus (EBV) DNA seronegative (Sero-) and seropositive (Sero+) nasopharyngeal carcinoma (NPC) are distinctly different disease subtypes. Patients with higher baseline EBV DNA titers seem to benefit less from anti-PD1 immunotherapy, but underlying mechanisms remain unclear. Tumor microenvironment (TME) characteristics could be the important factor affecting the efficacy of immunotherapy. Here, we illuminated the distinct multicellular ecosystems of EBV DNA Sero- and Sero+ NPCs from cellular compositional and functional perspectives at single-cell resolution.

Method

We performed single-cell RNA sequencing analyses of 28,423 cells from ten NPC samples and one non-tumor nasopharyngeal tissue. The markers, function, and dynamics of related cells were analyzed.

Results

We found that tumor cells from EBV DNA Sero+ samples exhibit low-differentiation potential, stronger stemness signature, and upregulated signaling pathways associated with cancer hallmarks than that of EBV DNA Sero- samples. Transcriptional heterogeneity and dynamics in T cells were associated with EBV DNA seropositivity status, indicating different immunoinhibitory mechanisms employed by malignant cells depending on EBV DNA seropositivity status. The low expression of classical immune checkpoints, early-triggered cytotoxic T-lymphocyte response, global activation of IFN-mediated signatures, and enhanced cell-cell interplays cooperatively tend to form a specific immune context in EBV DNA Sero+ NPC.

Conclusions

Collectively, we illuminated the distinct multicellular ecosystems of EBV DNA Sero- and Sero+ NPCs from single-cell perspective. Our study provides insights into the altered tumor microenvironment of NPC associated with EBV DNA seropositivity, which will help direct the development of rational immunotherapy strategies.

ZBTB12 is a molecular barrier to dedifferentiation in human pluripotent stem cells

Development is generally viewed as one-way traffic of cell state transition from primitive to developmentally advanced states. However, molecular mechanisms that ensure the unidirectional transition of cell fates remain largely unknown. Through exact transcription start site mapping, we report an evolutionarily conserved BTB domain-containing zinc finger protein, ZBTB12, as a molecular barrier for dedifferentiation of human pluripotent stem cells (hPSCs). Single-cell RNA sequencing reveals that ZBTB12 is essential for three germ layer differentiation by blocking hPSC dedifferentiation. Mechanistically, ZBTB12 fine-tunes the expression of human endogenous retrovirus H (HERVH), a primate-specific retrotransposon, and targets specific transcripts that utilize HERVH as a regulatory element. In particular, the downregulation of HERVH-overlapping long non-coding RNAs (lncRNAs) by ZBTB12 is necessary for a successful exit from a pluripotent state and lineage derivation. Overall, we identify ZBTB12 as a molecular barrier that safeguards the unidirectional transition of metastable stem cell fates toward developmentally advanced states.

Clonal cell states link Barrett's esophagus and esophageal adenocarcinoma

Barrett's esophagus is a common type of metaplasia and a precursor of esophageal adenocarcinoma. However, the cell states and lineage connections underlying the origin, maintenance, and progression of Barrett's esophagus have not been resolved in humans. To address this, we performed single-cell lineage tracing and transcriptional profiling of patient cells isolated from metaplastic and healthy tissue. Our analysis revealed discrete lineages in Barrett's esophagus, normal esophagus, and gastric cardia. Transitional basal progenitor cells of the gastroesophageal junction were unexpectedly related to both esophagus and gastric cardia cells. Barrett's esophagus was polyclonal, with lineages that contained all progenitor and differentiated cell types. In contrast, precancerous dysplastic foci were initiated by the expansion of a single molecularly aberrant Barrett's esophagus clone. Together, these findings provide a comprehensive view of the cell dynamics of Barrett's esophagus, linking cell states along the full disease trajectory, from its origin to cancer.

Asxl1 deletion disrupts MYC and RNA polymerase II function in granulocyte progenitors

Mutations in the gene Additional Sex-Combs Like 1 (ASXL1) are recurrent in myeloid malignancies as well as the pre-malignant condition clonal hematopoiesis, where they are universally associated with poor prognosis. However, the role of ASXL1 in myeloid lineage maturation is incompletely described. To define the role of ASXL1 in myelopoiesis, we employed single cell RNA sequencing and a murine model of hematopoietic-specific Asxl1 deletion. In granulocyte progenitors, Asxl1 deletion leads to hyperactivation of MYC and a quantitative decrease in neutrophil production. This loss of granulocyte production was not accompanied by significant changes in the landscape of covalent histone modifications. However, Asxl1 deletion results in a decrease in RNAPII promoter-proximal pausing in granulocyte progenitors, indicative of a global increase in productive transcription. These results suggest that ASXL1 inhibits productive transcription in granulocyte progenitors, identifying a new role for this epigenetic regulator in myeloid development.

The astroglial and stem cell functions of adult rat folliculostellate cells

The mammalian pituitary gland is a complex organ consisting of hormone-producing cells, anterior lobe folliculostellate cells (FSCs), posterior lobe pituicytes, vascular pericytes and endothelial cells, and Sox2-expressing stem cells. We present single-cell RNA sequencing and immunohistofluorescence analyses of pituitary cells of adult female rats with a focus on the transcriptomic profiles of nonhormonal cell types. Samples obtained from whole pituitaries and separated anterior and posterior lobe cells contained all expected pituitary resident cell types and lobe-specific vascular cell subpopulations. FSCs and pituicytes expressed S100B, ALDOC, EAAT1, ALDH1A1, and VIM genes and proteins, as well as other astroglial marker genes, some common and some cell type-specific. We also found that the SOX2 gene and protein were expressed in ~15% of pituitary cells, including FSCs, pituicytes, and a fraction of hormone-producing cells, arguing against its stem cell specificity. FSCs comprised two Sox2-expressing subclusters; FS1 contained more cells but lower genetic diversity, while FS2 contained proliferative cells, shared genes with hormone-producing cells, and expressed genes consistent with stem cell niche formation, regulation of cell proliferation and stem cell pluripotency, including the Hippo and Wnt pathways. FS1 cells were randomly distributed in the anterior and intermediate lobes, while FS2 cells were localized exclusively in the marginal zone between the anterior and intermediate lobes. These data indicate the identity of the FSCs as anterior pituitary-specific astroglia, with FS1 cells representing differentiated cells equipped for classical FSC roles and FS2 cells exhibiting additional stem cell-like features.

SUMO1 regulates post-infarct cardiac repair based on cellular heterogeneity

Small ubiquitin-related modifier (SUMOylation) is a dynamic post-translational modification that maintains cardiac function and can protect against a hypertrophic response to cardiac pressure overload. However, the function of SUMOylation after myocardial infarction (MI) and the molecular details of heart cell responses to SUMO1 deficiency have not been determined. In this study, we demonstrated that SUMO1 protein was inconsistently abundant in different cell types and heart regions after MI. However, SUMO1 knockout significantly exacerbated systolic dysfunction and infarct size after myocardial injury. Single-nucleus RNA sequencing revealed the differential role of SUMO1 in regulating heart cells. Among cardiomyocytes, SUMO1 deletion increased the Nppa ⁺ Nppb ⁺ Ankrd1 ⁺ cardiomyocyte subcluster proportion after MI. In addition, the conversion of fibroblasts to myofibroblasts subclusters was inhibited in SUMO1 knockout mice. Importantly, SUMO1 loss promoted proliferation of endothelial cell subsets with the ability to reconstitute neovascularization and expressed angiogenesis-related genes. Computational analysis of ligand/receptor interactions suggested putative pathways that mediate cardiomyocytes to endothelial cell communication in the myocardium. Mice preinjected with cardiomyocyte-specific AAV-SUMO1, but not the endothelial cell-specific form, and exhibited ameliorated cardiac remodeling following MI. Collectively, our results identified the role of SUMO1 in cardiomyocytes, fibroblasts, and endothelial cells after MI. These findings provide new insights into SUMO1 involvement in the pathogenesis of MI and reveal novel therapeutic targets.

Absolute scaling of single-cell transcriptomes identifies pervasive hypertranscription in adult stem and progenitor cells

Hypertranscription supports biosynthetically demanding cellular states through global transcriptome upregulation. Despite its potential widespread relevance, documented examples of hypertranscription remain few and limited to early development. Here, we demonstrate that absolute scaling of single-cell RNA-sequencing data enables the estimation of total transcript abundances per cell. We validate absolute scaling in known cases of developmental hypertranscription and apply it to adult cell types, revealing a remarkable dynamic range in transcriptional output. In adult organs, hypertranscription marks activated stem/progenitor cells with multilineage potential and is redeployed in conditions of tissue injury, where it precedes bursts of proliferation during regeneration. Our analyses identify a common set of molecular pathways associated with both adult and embryonic hypertranscription, including chromatin remodeling, DNA repair, ribosome biogenesis, and translation. These shared features across diverse cell contexts support hypertranscription as a general and dynamic cellular program that is pervasively employed during development, organ maintenance, and regeneration.

Single-cell RNA sequencing of solid pseudopapillary neoplasms of the pancreas in children

Solid pseudopapillary neoplasm (SPN) of the pancreas is a rare pancreatic tumor in children. Its origin remains elusive, along with its pathogenesis. Heterogeneity within SPN has not been previously described. In addition, low malignant but recurrent cases have occasionally been reported. To comprehensively unravel these profiles, single-cell RNA sequencing was performed using surgical specimens. We identified the cell types and suggested the origin of pancreatic endocrine progenitors. The Wnt/β-catenin pathway may be involved in tumorigenesis, while the epithelial-to-mesenchymal transition may be responsible for SPN recurrence. Furthermore, NOV, DCN were nominated as primary and S100A10, MGP as recurrent SPN marker genes, respectively. Our results provide insight into the pathogenesis of SPN.

Single-cell genomics in AML: extending the frontiers of AML research

The era of genomic medicine has allowed acute myeloid leukemia (AML) researchers to improve disease characterization, optimize risk-stratification systems, and develop new treatments. Although there has been significant progress, AML remains a lethal cancer because of its remarkably complex and plastic cellular architecture. This degree of heterogeneity continues to pose a major challenge, because it limits the ability to identify and therefore eradicate the cells responsible for leukemogenesis and treatment failure. In recent years, the field of single-cell genomics has led to unprecedented strides in the ability to characterize cellular heterogeneity, and it holds promise for the study of AML. In this review, we highlight advancements in single-cell technologies, outline important shortcomings in our understanding of AML biology and clinical management, and discuss how single-cell genomics can address these shortcomings as well as provide unique opportunities in basic and translational AML research.

Bioinformatic analysis of single-cell RNA sequencing dataset dissects cellular heterogeneity of triple-negative breast cancer in transcriptional profile, splicing event and crosstalk network

BACKGROUND

Triple-negative breast cancer (TNBC) is a subtype of breast cancer with high tumoral heterogeneity, while the detailed regulatory network is not well known.

METHODS

Via single-cell RNA-sequencing (scRNA-seq) data analysis, we comprehensively investigated the transcriptional profile of different subtypes of TNBC epithelial cells with gene regulatory network (GRN) and alternative splicing (AS) event analysis, as well as the crosstalk between epithelial and non-epithelial cells.

RESULTS

Of note, we found that luminal progenitor subtype exhibited the most complex GRN and splicing events. Besides, hnRNPs negatively regulates AS events in luminal progenitor subtype. In addition, we explored the cellular crosstalk among endothelial cells, stromal cells and immune cells in TNBC and discovered that NOTCH4 was a key receptor and prognostic marker in endothelial cells, which provide potential biomarker and target for TNBC intervention.

CONCLUSIONS

In summary, our study elaborates on the cellular heterogeneity of TNBC, revealing that NOTCH4 in endothelial cells was critical for TNBC intervention. This in-depth understanding of epithelial cell and non-epithelial cell network would provide theoretical basis for the development of new drugs targeting this sophisticated network in TNBC.

Paneth cells as the origin of intestinal cancer in the context of inflammation

Paneth cells (PCs), responsible for the secretion of antimicrobial peptides in the small intestine and for niche support to Lgr5+ crypt-base columnar stem cells (CBCs), have been shown to respond to inflammation by dedifferentiating into stem-like cells in order to sustain a regenerative response1,2. Therefore, PCs may represent the cells-of-origin of intestinal cancer in the context of inflammation. To test this hypothesis, we targeted Apc, Kras, and Tp53 mutations in Paneth cells by Cre-Lox technology and modelled inflammation by dextran sodium sulfate (DSS) administration. PC-specific loss of Apc resulted in multiple small intestinal tumors, whereas Kras or Tp53 mutations did not. Compound Apc and Kras mutations in PCs resulted in a striking increase in tumor multiplicity even in the absence of the inflammatory insult. By combining scRNAseq with lineage tracing to capture the conversion of PCs into bona fide tumor cells, we show that they progress through a "revival stem cell" (RSC) state characterized by high Clusterin (Clu) expression and Yap1 signaling, reminiscent of what has been previously observed upon irradiation of the mouse digestive tract3. Accordingly, comparison of PC- and Lgr5-derived murine intestinal tumors revealed differences related to Wnt signaling and inflammatory pathways which match the dichotomy of CBCs and injury-induced RSCs4 between human sporadic colon cancers and those arising in the context of inflammatory bowel diseases. Last, we show that western-style dietary habits, known to trigger a low-grade inflammation throughout the intestinal tract, underlie the analogous dedifferentiation of Paneth cells and their acquisition of stem-like features. Taken together, our results show that intestinal cancer arises in the context of inflammation through the dedifferentiation of committed secretory lineages such as Paneth cells and the activation of the revival stem cell state. As such, a true quiescent stem cell identity may be hidden in fully committed and postmitotic lineages which, upon inflammation, support the regenerative response by re-entering the cell cycle and dedifferentiating into RSCs. The chronic nature of the tissue insult in inflammatory bowel diseases and even in the context of western-style dietary habits is likely to result in the expansion of cell targets for tumor initiation and progression.

Multiplexed 3D atlas of state transitions and immune interaction in colorectal cancer

Advanced solid cancers are complex assemblies of tumor, immune, and stromal cells characterized by high intratumoral variation. We use highly multiplexed tissue imaging, 3D reconstruction, spatial statistics, and machine learning to identify cell types and states underlying morphological features of known diagnostic and prognostic significance in colorectal cancer. Quantitation of these features in high-plex marker space reveals recurrent transitions from one tumor morphology to the next, some of which are coincident with long-range gradients in the expression of oncogenes and epigenetic regulators. At the tumor invasive margin, where tumor, normal, and immune cells compete, T cell suppression involves multiple cell types and 3D imaging shows that seemingly localized 2D features such as tertiary lymphoid structures are commonly interconnected and have graded molecular properties. Thus, while cancer genetics emphasizes the importance of discrete changes in tumor state, whole-specimen imaging reveals large-scale morphological and molecular gradients analogous to those in developing tissues.

Characterization of immature ovarian teratomas through single-cell transcriptome

Introduction

Immature ovarian teratomas are a type of malignant germ cell tumor composed of complicated cell types and are characterized by pathological features of immature neuroectodermal tubules/rosettes. However, there is a lack of understanding of patient-derived immature ovarian teratomas (PDT) at the single cell level. Moreover, whether stem cell lines derived from immature teratomas (CDT) can be used as models for research on PDT remains to be elucidated.

Methods

Single-cell RNA sequencing (scRNA-seq) and subsequent bioinformatic analysis was performed on three patient-derived immature ovarian teratomas (PDT) samples to reveal the heterogeneity, evolution trajectory, and cell communication within the tumor microenvironment of PDT. Validations were conducted in additional seven samples through multiplex immunofluorescence.

Result

A total of qualified 22,153 cells were obtained and divided into 28 clusters, which can match to the scRNA-seq annotation of CDT as well as human fetal Cell Atlas, but with higher heterogeneity and more prolific cell-cell crosstalk. Radial glia cells (tagged by SOX2) and immature neuron (tagged by DCX) exhibited mutually exclusive expression and differentiated along distinct evolutionary trajectory from cycling neural progenitors. Proportions of these neuroectodermal cell subtypes may play important roles in PDT through contributing to the internal heterogeneity of PDTs. Moreover, the immune cells in PDTs were infiltrated rather than teratoma-derived, with more abundant macrophage in immature neuron than those in radial glia cells, and the infiltrated macrophage subtypes (i.e., M1 and M2) were significantly correlated to clinical grade. Overall, suppressed evolution process and transcriptome regulation in neuroectodermal cells, reduced cell-cell crosstalk, higher M1/M2 proportion ratio, and enhanced T cell effects in tumor microenvironment are enriched in patients with favorable prognosis.

Discussion

This study provides a comprehensive profile of PDT at the single cell level, shedding light on the heterogeneity and evolution of neuroectodermal cells within PDTs and the role of immune cells within the tumor microenvironment. Also, our findings highlight the potential usage of CDTs as a model for research on PDT.

Single-cell analysis reveals the chemotherapy-induced cellular reprogramming and novel therapeutic targets in relapsed/refractory acute myeloid leukemia

Chemoresistance and relapse are the leading cause of AML-related deaths. Utilizing single-cell RNA sequencing (scRNA-seq), we dissected the cellular states of bone marrow samples from primary refractory or short-term relapsed AML patients and defined the transcriptional intratumoral heterogeneity. We found that compared to proliferating stem/progenitor-like cells (PSPs), a subpopulation of quiescent stem-like cells (QSCs) were involved in the chemoresistance and poor outcomes of AML. By performing longitudinal scRNA-seq analyses, we demonstrated that PSPs were reprogrammed to obtain a QSC-like expression pattern during chemotherapy in refractory AML patients, characterized by the upregulation of CD52 and LGALS1 expression. Flow cytometric analysis further confirmed that the preexisting CD99⁺CD49d⁺CD52⁺Galectin-1⁺ (QSCs) cells at diagnosis were associated with chemoresistance, and these cells were further enriched in the residual AML cells of refractory patients. Interaction of CD52-SIGLEC10 between QSCs and monocytes may contribute to immune evading and poor outcomes. Furthermore, we identified that LGALS1 was a promising target for chemoresistant AML, and LGALS1 inhibitor could help eliminate QSCs and enhance the chemotherapy in patient-derived primary AML cells, cell lines, and AML xenograft models. Our results will facilitate a better understanding of the AML chemoresistance mechanism and the development of novel therapeutic strategies for relapsed/refractory AML patients.

Deciphering the sequential changes of monocytes/macrophages in the progression of IDD with longitudinal approach using single-cell transcriptome

Intervertebral disk degeneration (IDD) is a chronic inflammatory disease with intricate connections between immune infiltration and oxidative stress (OS). Complex cell niches exist in degenerative intervertebral disk (IVD) and interact with each other and regulate the disk homeostasis together. However, few studies have used longitudinal approach to describe the immune response of IDD progression. Here, we conducted conjoint analysis of bulk-RNA sequencing and single-cell sequencing, together with a series of techniques like weighted gene co-expression network analysis (WGCNA), immune infiltration analysis, and differential analysis, to systematically decipher the difference in OS-related functions of different cell populations within degenerative IVD tissues, and further depicted the longitudinal alterations of immune cells, especially monocytes/macrophages in the progression of IDD. The OS-related genes CYP1A1, MMP1, CCND1, and NQO1 are highly expressed and might be diagnostic biomarkers for the progression of IDD. Further landscape of IVD microenvironment showed distinct changes in cell proportions and characteristics at late degeneration compared to early degeneration of IDD. Monocytes/macrophages were classified into five distinct subpopulations with different roles. The trajectory lineage analysis revealed transcriptome alterations from effector monocytes/macrophages and regulatory macrophages to other subtypes during the evolution process and identified monocytes/macrophage subpopulations that had rapidly experienced the activation of inflammatory or anti-inflammatory responses. This study further proposed that personalized therapeutic strategies are needed to be formulated based on specific monocyte/macrophage subtypes and degenerative stages of IDD.

Identification of resident progenitors labeled with Top2a responsible for proximal tubular regeneration in ischemia reperfusion-induced acute kidney injury

BACKGROUND

Acute kidney injury is a common fatal disease with complex etiology and limited treatment methods. Proximal tubules (PTs) are the most vulnerable segment. Not only in injured kidneys but also in normal kidneys, shedding of PTs often happens. However, the source cells and mechanism of their regeneration remain unclear.

METHODS

ScRNA and snRNA sequencing data of acute injured or normal kidney were downloaded from GEO database to identify the candidate biomarker of progenitor of proximal tubules. SLICE algorithm and CytoTRACE analyses were employed to evaluate the stemness of progenitors. Then the repairing trajectory was constructed through pseudotime analyses. SCENIC algorithm was used to detect cell-type-specific regulon. With spatial transcriptome data, the location of progenitors was simulated. Neonatal/ adult/ aged mice and preconditioning AKI mice model and deconvolution of 2 RNA-seq data were employed for validation.

RESULTS

Through cluster identification, PT cluster expressed Top2a specifically was identified to increase significantly during AKI. With relatively strong stemness, the Top2a-labeled PT cluster tended to be the origin of the repairing trajectory. Moreover, the cluster was regulated by Pbx3-based regulon and possessed great segmental heterogeneity. Changes of Top2a between neonatal and aged mice and among AKI models validated the progenitor role of Top2a-labeled cluster.

CONCLUSIONS

Our study provided transcriptomic evidence that resident proximal tubular progenitors labeled with Top2a participated in regeneration. Considering the segmental heterogeneity, we find that there is a group of reserve progenitor cells in each tubular segment. When AKI occurs, the reserve progenitors of each tubular segment proliferate and replenish first, and PT-progenitors, a cluster with no obvious PT markers replenish each subpopulation of the reserve cells.

Functional and developmental heterogeneity of pituitary lactotropes in medaka

In fish, prolactin-producing cells (lactotropes) are located in the anterior part of the pituitary and play an essential role in osmoregulation. However, small satellite lactotrope clusters have been described in other parts of the pituitary in several species. The functional and developmental backgrounds of these satellite clusters are not known. We recently discovered two distinct prolactin-expressing cell types in Japanese medaka (Oryzias latipes), a euryhaline species, using single cell transcriptomics. In the present study, we characterize these two transcriptomically distinct lactotrope cell types and explore the hypothesis that they represent spatially distinct cell clusters, as found in other species. Single cell RNA sequencing shows that one of the two lactotrope cell types exhibits an expression profile similar to that of stem cell-like folliculo-stellate cell populations. Using in situ hybridization, we show that the medaka pituitary often develops additional small satellite lactotrope cell clusters, like in other teleost species. These satellite clusters arise early during development and grow in cell number throughout life regardless of the animal's sex. Surprisingly, our data do not show a correspondence between the stem cell-like lactotropes and these satellite lactotrope clusters. Instead, our data support a scenario in which the stem cell-like lactotropes are an intrinsic stage in the development of every spatially distinct lactotrope cluster. In addition, lactotrope activity in both spatially distinct lactotrope clusters decreases when environmental salinity increases, supporting their role in osmoregulation. However, this decrease appears weaker in the satellite lactotrope cell clusters, suggesting that these lactotropes are regulated differently.

Biology-inspired data-driven quality control for scientific discovery in single-cell transcriptomics

BACKGROUND

Quality control (QC) of cells, a critical first step in single-cell RNA sequencing data analysis, has largely relied on arbitrarily fixed data-agnostic thresholds applied to QC metrics such as gene complexity and fraction of reads mapping to mitochondrial genes. The few existing data-driven approaches perform QC at the level of samples or studies without accounting for biological variation.

RESULTS

We first demonstrate that QC metrics vary with both tissue and cell types across technologies, study conditions, and species. We then propose data-driven QC (ddqc), an unsupervised adaptive QC framework to perform flexible and data-driven QC at the level of cell types while retaining critical biological insights and improved power for downstream analysis. ddqc applies an adaptive threshold based on the median absolute deviation on four QC metrics (gene and UMI complexity, fraction of reads mapping to mitochondrial and ribosomal genes). ddqc retains over a third more cells when compared to conventional data-agnostic QC filters. Finally, we show that ddqc recovers biologically meaningful trends in gradation of gene complexity among cell types that can help answer questions of biological interest such as which cell types express the least and most number of transcripts overall, and ribosomal transcripts specifically.

CONCLUSIONS

ddqc retains cell types such as metabolically active parenchymal cells and specialized cells such as neutrophils which are often lost by conventional QC. Taken together, our work proposes a revised paradigm to quality filtering best practices-iterative QC, providing a data-driven QC framework compatible with observed biological diversity.

Single-Cell and CellChat Resolution Identifies Collecting Duct Cell Subsets and Their Communications with Adjacent Cells in PKD Kidneys

ADPKD is a genetic disorder with a molecular complexity that remains poorly understood. In this study, we sampled renal cells to construct a comprehensive and spatiotemporally resolved gene expression atlas in whole Pkd1 mutant polycystic mouse kidneys at single-cell resolution. We characterized cell diversity and identified novel collecting duct (CD) cell subtypes in cystic kidneys. We further found that CD cells appear to take different cell fate trajectories, and the first and the most important step might take place around day 14 in Pkd1 homozygous kidneys. After that day, increased numbers of CD cells showed highly proliferative and fibrotic characteristics, as detected in later-stage Pkd1 homozygous kidneys, both of which should contribute to cyst growth and renal fibrosis. With a newly developed modeling algorithm, called CellChat Explorer, we identify cell-to-cell communication networks mediated by the ligand receptor, such as MIF-CD44/CD74, in cystic kidneys, and confirm them via the expression patterns of ligands and receptors in four major cell types, which addresses the key question as to whether and how Pkd1 mutant renal epithelial cells affect their neighboring cells. The allele-specific gene expression profiles show that the secretion of cytokines by Pkd1 mutant epithelial cells may affect the gene expression profiles in recipient cells via epigenetic mechanisms, and vice versa. This study can be used to drive precision therapeutic targeting of ADPKD.

Pharmacological conversion of gut epithelial cells into insulin-producing cells lowers glycemia in diabetic animals

As a highly regenerative organ, the intestine is a promising source for cellular reprogramming for replacing lost pancreatic β cells in diabetes. Gut enterochromaffin cells can be converted to insulin-producing cells by forkhead box O1 (FoxO1) ablation, but their numbers are limited. In this study, we report that insulin-immunoreactive cells with Paneth/goblet cell features are present in human fetal intestine. Accordingly, lineage-tracing experiments show that, upon genetic or pharmacologic FoxO1 ablation, the Paneth/goblet lineage can also undergo conversion to the insulin lineage. We designed a screening platform in gut organoids to accurately quantitate β-like cell reprogramming and fine-tune a combination treatment to increase the efficiency of the conversion process in mice and human adult intestinal organoids. We identified a triple blockade of FOXO1, Notch, and TGF-β that, when tested in insulin-deficient streptozotocin (STZ) or NOD diabetic animals, resulted in near normalization of glucose levels, associated with the generation of intestinal insulin-producing cells. The findings illustrate a therapeutic approach for replacing insulin treatment in diabetes.

Cancer-selective metabolic vulnerabilities in MYC-amplified medulloblastoma

MYC-driven medulloblastoma (MB) is an aggressive pediatric brain tumor characterized by therapy resistance and disease recurrence. Here, we integrated data from unbiased genetic screening and metabolomic profiling to identify multiple cancer-selective metabolic vulnerabilities in MYC-driven MB tumor cells, which are amenable to therapeutic targeting. Among these targets, dihydroorotate dehydrogenase (DHODH), an enzyme that catalyzes de novo pyrimidine biosynthesis, emerged as a favorable candidate for therapeutic targeting. Mechanistically, DHODH inhibition acts on target, leading to uridine metabolite scarcity and hyperlipidemia, accompanied by reduced protein O-GlcNAcylation and c-Myc degradation. Pyrimidine starvation evokes a metabolic stress response that leads to cell-cycle arrest and apoptosis. We further show that an orally available small-molecule DHODH inhibitor demonstrates potent mono-therapeutic efficacy against patient-derived MB xenografts in vivo. The reprogramming of pyrimidine metabolism in MYC-driven medulloblastoma represents an unappreciated therapeutic strategy and a potential new class of treatments with stronger cancer selectivity and fewer neurotoxic sequelae.

Cancer immunotherapies transition endothelial cells into HEVs that generate TCF1+ T lymphocyte niches through a feed-forward loop

The lack of T cell infiltrates is a major obstacle to effective immunotherapy in cancer. Conversely, the formation of tumor-associated tertiary-lymphoid-like structures (TA-TLLSs), which are the local site of humoral and cellular immune responses against cancers, is associated with good prognosis, and they have recently been detected in immune checkpoint blockade (ICB)-responding patients. However, how these lymphoid aggregates develop remains poorly understood. By employing single-cell transcriptomics, endothelial fate mapping, and functional multiplex immune profiling, we demonstrate that antiangiogenic immune-modulating therapies evoke transdifferentiation of postcapillary venules into inflamed high-endothelial venules (HEVs) via lymphotoxin/lymphotoxin beta receptor (LT/LTβR) signaling. In turn, tumor HEVs boost intratumoral lymphocyte influx and foster permissive lymphocyte niches for PD1^- and PD1⁺TCF1⁺ CD8 T cell progenitors that differentiate into GrzB⁺PD1⁺ CD8 T effector cells. Tumor-HEVs require continuous CD8 and NK cell-derived signals revealing that tumor HEV maintenance is actively sculpted by the adaptive immune system through a feed-forward loop.

Dental niche cells directly contribute to tooth reconstitution and morphogenesis

Mammalian teeth develop from the inductive epithelial-mesenchymal interaction, an important mechanism shared by many organs. The cellular basis for such interaction remains elusive. Here, we generate a dual-fluorescence model to track and analyze dental cells from embryonic to postnatal stages, in which Pitx2⁺ epithelium and Msx1⁺ mesenchyme are sufficient for tooth reconstitution. Single-cell RNA sequencing and spatial mapping further revealed critical cellular dynamics during molar development, where tooth germs are organized by Msx1⁺Sdc1⁺ dental papilla and surrounding dental niche. Surprisingly, niche cells are more efficient in tooth reconstitution and can directly regenerate papilla cells through interaction with dental epithelium. Finally, from the dental niche, we identify a group of previously unappreciated migratory Msx1⁺ Sox9⁺ cells as the potential cell origin for dental papilla. Our results indicate that the dental niche cells directly contribute to tooth organogenesis and provide critical insights into the essential cell composition for tooth engineering.

Expressional and functional characteristics of checkpoint kinase 1 as a prognostic biomarker in hepatocellular carcinoma

Background

Hepatocellular carcinoma (HCC) is the most common pathological subtype of liver cancer and is the third leading cause of cancer death worldwide. Checkpoint kinase 1 (CHEK1), an essential serine/threonine kinase that regulates the cell cycle, is reported to be associated with carcinogenesis. However, the biological role and clinical significance of CHEK1 in HCC are still incompletely known.

Methods

In this research, CHEK1 messenger RNA (mRNA) levels in various liver hepatocellular carcinoma (LIHC) cohorts from the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) databases were evaluated. The Kaplan-Meier database was applied to identify the correlation between survival time and CHEK1 expression in patients with HCC. Gene set enrichment analysis (GSEA) was performed to explore the potential mechanism of CHEK1 in HCC, and NetworkAnalyst v. 3.0 (https://www.networkanalyst.ca/) was used to construct the regulatory networks of CHEK1 in HCC. Discriminant Regulon Expression Analysis (DoRothEA) was used to detect the activity of transcriptional factors (TFs) in gene-enriched cells (EC) with CHEK1 coexpression. In vitro experiments were conducted to investigate the effects of CHEK1 on the biological function of HCC cells.

Results

The CHEK1 mRNA level was overexpressed in HCC, and increased CHEK1 expression correlated with poor survival outcomes. The homo sapiens-microRNA-195 (hsa-miR-195) may have contributed to the upregulation of CHEK1 in HCC. GSEA and NetworkAnalyst v. 3.0 showed that CHEK1 played a crucial part in tumor proliferation of HCC and may be regulated by TF E2F1. DoRothEA showed increased transcriptional activity of E2F1 in gene-EC with CHEK1 coexpression. Moreover, experiments of cell function showed that the knockdown of CHEK1 weakened the aggressive behavior and proliferation of HCC cells. Overexpression of E2F1 increased the proliferation and invasion of HCC cells in vitro, while the silencing of CHEK1 dampened cell invasion induced by E2F1 overexpression.

Conclusions

These results identified the prognostic significance and expression characteristics of CHEK1 in HCC through bioinformatics analysis and experimental verification. This lays the foundation for further research on the diagnosis and treatment of HCC.

Cellular interactions in the pituitary stem cell niche

Stem cells in the anterior pituitary gland can give rise to all resident endocrine cells and are integral components for the appropriate development and subsequent maintenance of the organ. Located in discreet niches within the gland, stem cells are involved in bi-directional signalling with their surrounding neighbours, interactions which underpin pituitary gland homeostasis and response to organ challenge or physiological demand. In this review we highlight core signalling pathways that steer pituitary progenitors towards specific endocrine fate decisions throughout development. We further elaborate on those which are conserved in the stem cell niche postnatally, including WNT, YAP/TAZ and Notch signalling. Furthermore, we have collated a directory of single cell RNA sequencing studies carried out on pituitaries across multiple organisms, which have the potential to provide a vast database to study stem cell niche components in an unbiased manner. Reviewing published data, we highlight that stem cells are one of the main signalling hubs within the anterior pituitary. In future, coupling single cell sequencing approaches with genetic manipulation tools in vivo, will enable elucidation of how previously understudied signalling pathways function within the anterior pituitary stem cell niche.

Liver tumour immune microenvironment subtypes and neutrophil heterogeneity

The heterogeneity of the tumour immune microenvironment (TIME), organized by various immune and stromal cells, is a major contributing factor of tumour metastasis, relapse and drug resistance^1-3, but how different TIME subtypes are connected to the clinical relevance in liver cancer remains unclear. Here we performed single-cell RNA-sequencing (scRNA-seq) analysis of 189 samples collected from 124 patients and 8 mice with liver cancer. With more than 1 million cells analysed, we stratified patients into five TIME subtypes, including immune activation, immune suppression mediated by myeloid or stromal cells, immune exclusion and immune residence phenotypes. Different TIME subtypes were spatially organized and associated with chemokine networks and genomic features. Notably, tumour-associated neutrophil (TAN) populations enriched in the myeloid-cell-enriched subtype were associated with an unfavourable prognosis. Through in vitro induction of TANs and ex vivo analyses of patient TANs, we showed that CCL4⁺ TANs can recruit macrophages and that PD-L1⁺ TANs can suppress T cell cytotoxicity. Furthermore, scRNA-seq analysis of mouse neutrophil subsets revealed that they are largely conserved with those of humans. In vivo neutrophil depletion in mouse models attenuated tumour progression, confirming the pro-tumour phenotypes of TANs. With this detailed cellular heterogeneity landscape of liver cancer, our study illustrates diverse TIME subtypes, highlights immunosuppressive functions of TANs and sheds light on potential immunotherapies targeting TANs.

X-box binding protein 1 as a key modulator in “healing endothelial cells”, a novel EC phenotype promoting angiogenesis after MCAO

Background

Endothelial cells (ECs) play an important role in angiogenesis and vascular reconstruction in the pathophysiology of ischemic stroke. Previous investigations have provided a profound cerebral vascular atlas under physiological conditions, but have failed to identify new disease-related cell subtypes. We aimed to identify new EC subtypes and determine the key modulator genes.

Methods

Two datasets GSE174574 and GSE137482 were included in the study. Seurat was utilized as the standard quality-control pipeline. UCell was used to calculate single-cell scores to validate cellular identity. Monocle3 and CytoTRACE were utilized in aid of pseudo-time differentiation analysis. CellChat was utilized to infer the intercellular communication pathways. The angiogenesis ability of ECs was validated by MTS, Transwell, tube formation, flow cytometry, and immunofluorescence assays in vitro and in vivo. A synchrotron radiation-based propagation contrast imaging was introduced to comprehensively portray cerebral vasculature.

Results

We successfully identified a novel subtype of EC named “healing EC” that highly expressed pan-EC marker and pro-angiogenic genes but lowly expressed all the arteriovenous markers identified in the vascular single-cell atlas. Further analyses showed its high stemness to differentiate into other EC subtypes and potential to modulate inflammation and angiogenesis via excretion of signal molecules. We therefore identified X-box binding protein 1 (Xbp1) as a key modulator in the healing EC phenotype. In vitro and in vivo experiments confirmed its pro-angiogenic roles under both physiological and pathological conditions. Synchrotron radiation-based propagation contrast imaging further proved that Xbp1 could promote angiogenesis and recover normal vasculature conformation, especially in the corpus striatum and prefrontal cortex under middle cerebral artery occlusion (MCAO) condition.

Conclusions

Our study identified a novel disease-related EC subtype that showed high stemness to differentiate into other EC subtypes. The predicted molecule Xbp1 was thus confirmed as a key modulator that can promote angiogenesis and recover normal vasculature conformation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s11658-022-00399-5.

Estimating Metastatic Risk of Pancreatic Ductal Adenocarcinoma at Single-Cell Resolution

Pancreatic ductal adenocarcinoma (PDAC) is characterized by intra-tumoral heterogeneity, and patients are always diagnosed after metastasis. Thus, finding out how to effectively estimate metastatic risk underlying PDAC is necessary. In this study, we proposed scMetR to evaluate the metastatic risk of tumor cells based on single-cell RNA sequencing (scRNA-seq) data. First, we identified diverse cell types, including tumor cells and other cell types. Next, we grouped tumor cells into three sub-populations according to scMetR score, including metastasis-featuring tumor cells (MFTC), transitional metastatic tumor cells (TransMTC), and conventional tumor cells (ConvTC). We identified metastatic signature genes (MSGs) through comparing MFTC and ConvTC. Functional enrichment analysis showed that up-regulated MSGs were enriched in multiple metastasis-associated pathways. We also found that patients with high expression of up-regulated MSGs had worse prognosis. Spatial mapping of MFTC showed that they are preferentially located in the cancer and duct epithelium region, which was enriched with the ductal cells' associated inflammation. Further, we inferred cell-cell interactions, and observed that interactions of the ADGRE5 signaling pathway, which is associated with metastasis, were increased in MFTC compared to other tumor sub-populations. Finally, we predicted 12 candidate drugs that had the potential to reverse expression of MSGs. Taken together, we have proposed scMetR to estimate metastatic risk in PDAC patients at single-cell resolution which might facilitate the dissection of tumor heterogeneity.

Cell type diversity in a developing octopus brain

Octopuses are mollusks that have evolved intricate neural systems comparable with vertebrates in terms of cell number, complexity and size. The brain cell types that control their sophisticated behavioral repertoire are still unknown. Here, we profile the cell diversity of the paralarval Octopus vulgaris brain to build a cell type atlas that comprises mostly neural cells, but also multiple glial subtypes, endothelial cells and fibroblasts. We spatially map cell types to the vertical, subesophageal and optic lobes. Investigation of cell type conservation reveals a shared gene signature between glial cells of mouse, fly and octopus. Genes related to learning and memory are enriched in vertical lobe cells, which show molecular similarities with Kenyon cells in Drosophila. We construct a cell type taxonomy revealing transcriptionally related cell types, which tend to appear in the same brain region. Together, our data sheds light on cell type diversity and evolution in the octopus brain.

Density-based detection of cell transition states to construct disparate and bifurcating trajectories

Tree- and linear-shaped cell differentiation trajectories have been widely observed in developmental biologies and can be also inferred through computational methods from single-cell RNA-sequencing datasets. However, trajectories with complicated topologies such as loops, disparate lineages and bifurcating hierarchy remain difficult to infer accurately. Here, we introduce a density-based trajectory inference method capable of constructing diverse shapes of topological patterns including the most intriguing bifurcations. The novelty of our method is a step to exploit overlapping probability distributions to identify transition states of cells for determining connectability between cell clusters, and another step to infer a stable trajectory through a base-topology guided iterative fitting. Our method precisely re-constructed various benchmark reference trajectories. As a case study to demonstrate practical usefulness, our method was tested on single-cell RNA sequencing profiles of blood cells of SARS-CoV-2-infected patients. We not only re-discovered the linear trajectory bridging the transition from IgM plasmablast cells to developing neutrophils, and also found a previously-undiscovered lineage which can be rigorously supported by differentially expressed gene analysis.

Endocardium-to-coronary artery differentiation during heart development and regeneration involves sequential roles of Bmp2 and Cxcl12/Cxcr4

Endocardial cells lining the heart lumen are coronary vessel progenitors during embryogenesis. Re-igniting this developmental process in adults could regenerate blood vessels lost during cardiac injury, but this requires additional knowledge of molecular mechanisms. Here, we use mouse genetics and scRNA-seq to identify regulators of endocardial angiogenesis and precisely assess the role of CXCL12/CXCR4 signaling. Time-specific lineage tracing demonstrated that endocardial cells differentiated into coronary endothelial cells primarily at mid-gestation. A new mouse line reporting CXCR4 activity-along with cell-specific gene deletions-demonstrated it was specifically required for artery morphogenesis rather than angiogenesis. Integrating scRNA-seq data of endocardial-derived coronary vessels from mid- and late-gestation identified a Bmp2-expressing transitioning population specific to mid-gestation. Bmp2 stimulated endocardial angiogenesis in vitro and in injured neonatal mouse hearts. Our data demonstrate how understanding the molecular mechanisms underlying endocardial angiogenesis can identify new potential therapeutic targets promoting revascularization of the injured heart.

A multi-omics approach to visualize early neuronal differentiation from hESCs in 4D

Neuronal differentiation of pluripotent stem cells is an established method to study physiology, disease, and medication safety. However, the sequence of events in human neuronal differentiation and the ability of in vitro models to recapitulate early brain development are poorly understood. We developed a protocol optimized for the study of early human brain development and neuropharmacological applications. We comprehensively characterized gene expression and epigenetic profiles at four timepoints, because the cells differentiate from embryonic stem cells towards a heterogeneous population of progenitors, immature and mature neurons bearing telencephalic signatures. A multi-omics roadmap of neuronal differentiation, combined with searchable interactive gene analysis tools, allows for extensive exploration of early neuronal development and the effect of medications.

Complex interaction and heterogeneity among cancer stem cells in head and neck squamous cell carcinoma revealed by single-cell sequencing

BACKGROUND

Cancer stem cells (CSCs) have been characterized to be responsible for multidrug resistance, metastasis, recurrence, and immunosuppressive in head and neck squamous cell carcinoma (HNSCC). However, the diversity of CSCs remains to be investigated. In this study, we aimed to determine the heterogeneity of CSCs and its effect on the formation of tumor microenvironment (TME).

METHODS

We depicted the landscape of HNSCC transcriptome profile by single-cell RNA-sequencing analysis of 20 HNSCC tissues from public databases, to reveal the Cell components, trajectory changes, signaling network, malignancy status and functional enrichment of CSCs within tumors.

RESULTS

Immune checkpoint molecules CD276, LILRB2, CD47 were significantly upregulated in CSCs, enabling host antitumor response to be weakened or damaged. Notably, naive CSCs were divided to 2 different types of cells with different functions, exhibiting functional diversity. In addition, CSCs underwent self-renewal and tumor metastasis activity through WNT and ncWNT signaling. Among them, Regulon regulators (IRF1_394g, IRF7_160g, NFKB1_12g, NFKB2_33g and STAT1_356g) were activated in subgroups 2 and 3, suggesting their pivotal roles in the inflammatory response process in tumors. Among all CSCs, naive CSCs appear to be the most malignant resulting in a worse prognosis.

CONCLUSIONS

Our study reveals the major signal transduction and biological function of CSCs during HNSCC progression, highlighting the heterogeneity of CSCs and their underlying mechanisms in the formation of an immunosuppressive TME. Therefore, our study about heterogeneity of CSCs in HNSCC can bring new insights for the treatment of HNSCC.

Multiomic analysis reveals conservation of cancer-associated fibroblast phenotypes across species and tissue of origin

Cancer-associated fibroblasts (CAFs) are integral to the solid tumor microenvironment. CAFs were once thought to be a relatively uniform population of matrix-producing cells, but single-cell RNA sequencing has revealed diverse CAF phenotypes. Here, we further probed CAF heterogeneity with a comprehensive multiomics approach. Using paired, same-cell chromatin accessibility and transcriptome analysis, we provided an integrated analysis of CAF subpopulations over a complex spatial transcriptomic and proteomic landscape to identify three superclusters: steady state-like (SSL), mechanoresponsive (MR), and immunomodulatory (IM) CAFs. These superclusters are recapitulated across multiple tissue types and species. Selective disruption of underlying mechanical force or immune checkpoint inhibition therapy results in shifts in CAF subpopulation distributions and affected tumor growth. As such, the balance among CAF superclusters may have considerable translational implications. Collectively, this research expands our understanding of CAF biology, identifying regulatory pathways in CAF differentiation and elucidating therapeutic targets in a species- and tumor-agnostic manner.

Identification of a minority population of LMO2+ breast cancer cells that integrate into the vasculature and initiate metastasis

Metastasis is responsible for most breast cancer-related deaths; however, identifying the cellular determinants of metastasis has remained challenging. Here, we identified a minority population of immature THY1+/VEGFA+ tumor epithelial cells in human breast tumor biopsies that display angiogenic features and are marked by the expression of the oncogene, LMO2. Higher abundance of LMO2+ basal cells correlated with tumor endothelial content and predicted poor distant recurrence-free survival in patients. Using MMTV-PyMT/Lmo2CreERT2 mice, we demonstrated that Lmo2 lineage-traced cells integrate into the vasculature and have a higher propensity to metastasize. LMO2 knockdown in human breast tumors reduced lung metastasis by impairing intravasation, leading to a reduced frequency of circulating tumor cells. Mechanistically, we find that LMO2 binds to STAT3 and is required for STAT3 activation by tumor necrosis factor-α and interleukin-6. Collectively, our study identifies a population of metastasis-initiating cells with angiogenic features and establishes the LMO2-STAT3 signaling axis as a therapeutic target in breast cancer metastasis.

Human prefrontal cortex gene regulatory dynamics from gestation to adulthood at single-cell resolution

Human brain development is underpinned by cellular and molecular reconfigurations continuing into the third decade of life. To reveal cell dynamics orchestrating neural maturation, we profiled human prefrontal cortex gene expression and chromatin accessibility at single-cell resolution from gestation to adulthood. Integrative analyses define the dynamic trajectories of each cell type, revealing major gene expression reconfiguration at the prenatal-to-postnatal transition in all cell types followed by continuous reconfiguration into adulthood and identifying regulatory networks guiding cellular developmental programs, states, and functions. We uncover links between expression dynamics and developmental milestones, characterize the diverse timing of when cells acquire adult-like states, and identify molecular convergence from distinct developmental origins. We further reveal cellular dynamics and their regulators implicated in neurological disorders. Finally, using this reference, we benchmark cell identities and maturation states in organoid models. Together, this captures the dynamic regulatory landscape of human cortical development.