Single-cell trajectories reconstruction, exploration and mapping of omics data with STREAM

scRDEN: single-cell dynamic gene rank differential expression network and robust trajectory inference

The remarkable advancement of single-cell RNA sequencing (scRNA-seq) technology has empowered researchers to probe gene expression at the single-cell level with unprecedented precision. To gain a profound understanding of the heterogeneity inherent in cell fate determination, a central challenge lies in the comprehensive analysis of the dynamic regulatory alterations that underlie transcriptional differences and the accurate inference of the differentiation trajectory. Here, we propose the method scRDEN, a robust framework that infers important cell sub-populations and differential expression networks of multiple genes along the differentiation directions of each branch by converting the unstable gene expression values in cells into relatively stable gene-gene interactions (global features) and extracting the order of differential expression (network features), and further integrating the expression features of different dimension reduction methods. When applied to five published scRNA-seq datasets from human and mouse cell differentiation, scRDEN not only successfully captures the stable cell subpopulations with potential marker genes, measures the transcriptional differences of gene pairs to identify the rank differential expression network along the differentiation direction of each branch. In addition, in multiple gene rank differential expression networks, the rank expression directly related to transcription factors/marker genes shows a significant strengthening and weakening trend along with their expression changes, and the distribution of diversity and cluster coefficient show a non-monotonic change trend, including the cases of increasing first and then decreasing or decreasing first and then increasing. This may correspond to the mechanism of cells gradually differentiating into stable functions. It is particularly noteworthy that scRDEN method yielded exceptional results when applied to the large-scale, multi-branched, double-batch mouse dentate gyrus data. This outstanding performance provides novel and valuable insights into large-scale, multi-batch trajectory inference and the study of transcriptional mechanism regulation during the processes of differentiation and development.

Quantitative CT Scoring for Local COPD Severity

Chronic obstructive pulmonary disease (COPD) is complex, and its course is difficult to predict due to its diverse pathophysiology. Small airway disease (SAD), a key component of COPD and potential target for emerging therapeutics, may be reversible in mild COPD, but left unchecked, may worsen, leading to airway loss and emphysema. The dual nature of SAD complicates clinical management of COPD patients, necessitating more accurate monitoring methods. To meet this need, we developed elastic Parametric Response Mapping (ePRM), a tiered scoring system that classifies local lung volumes by the degree of PRM-derived SAD, normal, and emphysematous tissue. In individuals with or at risk for COPD, we demonstrate that chest CT ePRM can categorize local lung tissue into distinct tiers of disease severity that distinguish between tissue characterized by early reversible SAD and progressive destruction. This level of characterization is crucial to developing personalized treatment strategies for COPD.

Small molecule BLVRB redox inhibitor promotes megakaryocytopoiesis and stress thrombopoiesis in vivo

Biliverdin IXβ reductase (BLVRB) is an NADPH-dependent enzyme previously implicated in a redox-regulated mechanism of thrombopoiesis distinct from the thrombopoietin (TPO)/c-MPL axis. Here, we apply computational modeling to inform molecule design, followed by de novo syntheses and screening of unique small molecules retaining the capacity for selective BLVRB inhibition as a novel platelet-enhancing strategy. Two distinct classes of molecules are identified, and NMR spectroscopy and co-crystallization studies confirm binding modes within the BLVRB active site and ring stacking between the nicotinamide moiety of the NADP+ cofactor. A diazabicyclo derivative displaying minimal off-target promiscuity and excellent bioavailability characteristics promotes megakaryocyte speciation in biphenotypic (erythro/megakaryocyte) cellular models and synergizes with TPO-dependent megakaryocyte formation in hematopoietic stem cells. Upon oral delivery into mice, this inhibitor expands platelet recovery in stress thrombopoietic models with no adverse effects. In this work, we identify and validate a cellular redox inhibitor retaining the potential to selectively promote megakaryocytopoiesis and enhance stress-associated platelet formation in vivo distinct from TPO receptor agonists.

Evolutionary fingerprints of epithelial-to-mesenchymal transition

Mesenchymal plasticity has been extensively described in advanced epithelial cancers; however, its functional role in malignant progression is controversial1-5. The function of epithelial-to-mesenchymal transition (EMT) and cell plasticity in tumour heterogeneity and clonal evolution is poorly understood. Here we clarify the contribution of EMT to malignant progression in pancreatic cancer. We used somatic mosaic genome engineering technologies to trace and ablate malignant mesenchymal lineages along the EMT continuum. The experimental evidence clarifies the essential contribution of mesenchymal lineages to pancreatic cancer evolution. Spatial genomic analysis, single-cell transcriptomic and epigenomic profiling of EMT clarifies its contribution to the emergence of genomic instability, including events of chromothripsis. Genetic ablation of mesenchymal lineages robustly abolished these mutational processes and evolutionary patterns, as confirmed by cross-species analysis of pancreatic and other human solid tumours. Mechanistically, we identified that malignant cells with mesenchymal features display increased chromatin accessibility, particularly in the pericentromeric and centromeric regions, in turn resulting in delayed mitosis and catastrophic cell division. Thus, EMT favours the emergence of genomic-unstable, highly fit tumour cells, which strongly supports the concept of cell-state-restricted patterns of evolution, whereby cancer cell speciation is propagated to progeny within restricted functional compartments. Restraining the evolutionary routes through ablation of clones capable of mesenchymal plasticity, and extinction of the derived lineages, halts the malignant potential of one of the most aggressive forms of human cancer.

Temporal Exploration of Chronic Obstructive Pulmonary Disease Phenotypes: Insights from the COPDGene and SPIROMICS Cohorts

Rationale: Chronic obstructive pulmonary disease (COPD) exhibits considerable progression heterogeneity. We hypothesized that elastic principal graph analysis (EPGA) would identify distinct clinical phenotypes and their longitudinal relationships. Objectives: Our primary objective was to create a map of COPD phenotypes and their connectivity using EPGA. Secondarily, we used longitudinal and external data sets to test the validity and reproducibility of this map. Methods: Cross-sectional data from 8,972 tobacco-exposed COPDGene participants, with and without COPD, were used to train a model with EPGA, using thirty clinical, physiologic and CT features. 4,585 participants from COPDGene Phase 2 were used to test longitudinal trajectories. 2,652 participants from SPIROMICS tested external reproducibility. Measurements and Main Results: Our analysis used cross-sectional data to create an elastic principal tree, where time is associated with distance on the tree. Six clinically distinct tree segments were identified that differed by lung function, symptoms, and CT features: Subclinical (SC); Parenchymal Abnormality (PA); Chronic Bronchitis (CB); Emphysema Male (EM); Emphysema Female (EF); and Severe Airways (SA) disease. 5-year data from COPDGene mapped longitudinal changes onto the tree, and longitudinal trajectories demonstrated a net flow of patients from SC towards EM and EF, including trajectories through airway disease predominant phenotypes, CB and SA. Cross-sectional SPIROMICS data projected onto the tree showed clinically similar patient groupings. Conclusions: This novel analytic methodology provides an approach to defining longitudinal phenotypic trajectories using cross sectional data. These insights are clinically relevant and could facilitate precision therapy and future trials to modify disease progression. Clinical trial registered with www.clinicaltrials.gov (NCT00608764 and NCT01969344).

A feedback amplifier circuit with Notch and E2A orchestrates T-cell fate and suppresses the innate lymphoid cell lineages during thymic ontogeny

External signals from the thymic microenvironment and the activities of lineage-specific transcription factors (TFs) instruct T-cell versus innate lymphoid cell (ILC) fates. However, mechanistic insights into how factors such as Notch1-Delta-like-4 (Dll4) signaling and E-protein TFs collaborate to establish T-cell identity remain rudimentary. Using multiple in vivo approaches and single-cell multiome analysis, we identified a feedback amplifier circuit that specifies fetal and adult T-cell fates. In early T progenitors (ETPs) in the fetal thymus, Notch signaling minimally lowered E-protein antagonist Id2 levels, and high Id2 abundance favored the differentiation of ETPs into ILCs. Conversely, in the adult thymus, Notch signaling markedly decreased Id2 abundance in ETPs, substantially elevating E-protein DNA binding and in turn promoting the activation of a T-cell lineage-specific gene expression program linked with V(D)J gene recombination and T-cell receptor signaling. Our findings indicate that, in the fetal versus the adult thymus, a simple feedback amplifier circuit dictated by Notch-mediated signals and Id2 abundance enforces T-cell identity and suppresses ILC development.

Fuzzy-Based Identification of Transition Cells to Infer Cell Trajectory for Single-Cell Transcriptomics

With the continuous evolution of single-cell RNA sequencing technology, it has become feasible to reconstruct cell development processes using computational methods. Trajectory inference is a crucial downstream analytical task that provides valuable insights into understanding cell cycle and differentiation. During cell development, cells exhibit both stable and transition states, which makes it challenging to accurately identify these cells. To address this challenge, we propose a novel single-cell trajectory inference method using fuzzy clustering, named scFCTI. By introducing fuzzy clustering and quantifying cell uncertainty, scFCTI can identify transition cells within unstable cell states. Moreover, scFCTI can obtain refined cell classification by characterizing different cell stages, which gain more accurate single-cell trajectory reconstruction containing transition paths. To validate the effectiveness of scFCTI, we conduct experiments on five real datasets and four different structure simulation datasets, comparing them with several state-of-the-art trajectory inference methods. The results demonstrate that scFCTI outperforms these methods by successfully identifying unstable cell clusters and obtaining more accurate cell paths with transition states. Especially the experimental results demonstrate that scFCTI can reconstruct the cell trajectory more precisely.

Dynamic Modeling of Cell Cycle Arrest Through Integrated Single-Cell and Mathematical Modelling Approaches

Highly multiplexed imaging assays allow simultaneous quantification of multiple protein and phosphorylation markers, providing a static snapshots of cell types and states. Pseudo-time techniques can transform these static snapshots of unsynchronized cells into dynamic trajectories, enabling the study of dynamic processes such as development trajectories and the cell cycle. Such ordering also enables training of mathematical models on these data, but technical challenges have hitherto made it difficult to integrate multiple experimental conditions, limiting the predictive power and insights these models can generate. In this work, we propose data processing and model training approaches for integrating multiplexed, multi-condition immunofluorescence data with mathematical modelling. We devise training strategies that are applicable to datasets where cells exhibit oscillatory as well as arrested dynamics and use them to train a cell cycle model on a dataset of MCF-10A mammary epithelial exposed to cell-cycle arresting small molecules. We validate the model by investigating predicted growth factor sensitivities and responses to inhibitors of cells at different initial conditions. We anticipate that our framework will generalise to other highly multiplexed measurement techniques such as mass-cytometry, rendering larger bodies of data accessible to dynamic modelling and paving the way to deeper biological insights.

Single-cell delineation of the microbiota-gut-brain axis: Probiotic intervention in Chd8 haploinsufficient mice

Emerging research underscores the gut microbiome's impact on the nervous system via the microbiota-gut-brain axis, yet comprehensive insights remain limited. Using a CHD8-haploinsufficient model for autism spectrum disorder (ASD), we explored host-gut microbiota interactions by constructing a single-cell transcriptome atlas of brain and intestinal tissues in wild-type and mutant mice across three developmental stages. CHD8 haploinsufficiency caused delayed development of radial glial precursors and excitatory neural progenitors in the E14.5 brain, inflammation in the adult brain, immunodeficiency, and abnormal intestinal development. Selective CHD8 knockdown in intestinal epithelial cells generated Chd8ΔIEC mice, which exhibited normal sociability but impaired social novelty recognition. Probiotic intervention with Lactobacillus murinus selectively rescued social deficits in Chd8ΔIEC mice, with single-cell transcriptome analysis revealing underlying mechanisms. This study provides a detailed single-cell transcriptomic dataset of ASD-related neural and intestinal changes, advancing our understanding of the gut-brain axis and offering potential therapeutic strategies for ASD.

Single-cell omics: experimental workflow, data analyses and applications

Cells are the fundamental units of biological systems and exhibit unique development trajectories and molecular features. Our exploration of how the genomes orchestrate the formation and maintenance of each cell, and control the cellular phenotypes of various organismsis, is both captivating and intricate. Since the inception of the first single-cell RNA technology, technologies related to single-cell sequencing have experienced rapid advancements in recent years. These technologies have expanded horizontally to include single-cell genome, epigenome, proteome, and metabolome, while vertically, they have progressed to integrate multiple omics data and incorporate additional information such as spatial scRNA-seq and CRISPR screening. Single-cell omics represent a groundbreaking advancement in the biomedical field, offering profound insights into the understanding of complex diseases, including cancers. Here, we comprehensively summarize recent advances in single-cell omics technologies, with a specific focus on the methodology section. This overview aims to guide researchers in selecting appropriate methods for single-cell sequencing and related data analysis.

Cancer-specific innate and adaptive immune rewiring drives resistance to PD-1 blockade in classic Hodgkin lymphoma

Hodgkin Reed-Sternberg (HRS) cells of classic Hodgkin lymphoma (cHL), like many solid tumors, elicit ineffective immune responses. However, patients with cHL are highly responsive to PD-1 blockade, which largely depends on HRS cell-specific retention of MHC class II and implicates CD4+ T cells and additional MHC class I-independent immune effectors. Here, we utilize single-cell RNA sequencing and spatial analysis to define shared circulating and microenvironmental features of the immune response to PD-1 blockade in cHL. Compared with non-responders, responding patients have more circulating CD4+ naïve and central memory T cells and B cells, as well as more diverse CD4+ T cell and B cell receptor repertoires. Importantly, a population of circulating and tumor-infiltrating IL1β+ monocytes/macrophages is detectable in patients with cHL but not healthy donors, and a proinflammatory, tumor-promoting signature of these circulating IL1β+ monocytes is associated with resistance to PD-1 blockade in cHL. Altogether, our findings reveal extensive immune rewiring and complementary roles of CD4+ T cells, B cells and IL1β+ monocytes in the response to PD-1 blockade and suggest that these features can be captured with a peripheral blood test.

Single-cell and spatial transcriptomics reveal SPP1-CD44 signaling drives primary resistance to immune checkpoint inhibitors in RCC

BACKGROUND

Immune checkpoint inhibitors (ICIs) are a cornerstone therapy for advanced renal cell carcinoma (RCC). However, significant rates of primary resistance hinder their efficacy, and the underlying mechanisms remain poorly understood. This study aims to unravel the tumor-immune interactions and signaling pathways driving primary resistance to ICIs in RCC.

METHODS

We integrated single-cell RNA sequencing, spatial transcriptomics, and clinical sample analysis to investigate the tumor microenvironment and intercellular signaling. Advanced computational methods, including cell-cell communication networks, pseudotime trajectories, and gene set enrichment analysis (GSEA), were employed to uncover the underlying resistance mechanisms.

RESULTS

Compared to the sensitive group, the primary resistance group exhibited a significant increase in SPP1-CD44 signaling-mediated interactions between tumor cells and immune cells. These interactions disrupted antigen presentation in immune effector cells and suppressed key chemokine and cytokine pathways, thereby impairing effective immune responses. In contrast, the sensitive group showed more active antigen presentation and cytokine signaling, which facilitated stronger immune responses. Furthermore, the interaction between SPP1-secreting tumor cells and CD44-expressing exhausted CD8 + T cells activated the MAPK signaling pathway within CD8 + Tex cells, exacerbating T cell exhaustion and driving the development of ICI resistance in RCC.

CONCLUSION

Our findings reveal a potential mechanism by which SPP1-CD44 signaling mediates tumor-immune cell interactions leading to ICI resistance, providing a theoretical basis for targeting and disrupting this signaling to overcome primary resistance in RCC.

The frontier of precision medicine: application of single-cell multi-omics in preimplantation genetic diagnosis

The advent of single-cell multi-omics technologies has revolutionized the landscape of preimplantation genetic diagnosis (PGD), offering unprecedented insights into the genetic, transcriptomic, and proteomic profiles of individual cells in early-stage embryos. This breakthrough holds the promise of enhancing the accuracy, efficiency, and scope of PGD, thereby significantly improving outcomes in assisted reproductive technologies (ARTs) and genetic disease prevention. This review provides a comprehensive overview of the importance of PGD in the context of precision medicine and elucidates how single-cell multi-omics technologies have transformed this field. We begin with a brief history of PGD, highlighting its evolution and application in detecting genetic disorders and facilitating ART. Subsequently, we delve into the principles, methodologies, and applications of single-cell genomics, transcriptomics, and proteomics in PGD, emphasizing their role in improving diagnostic precision and efficiency. Furthermore, we review significant recent advances within this domain, including key experimental designs, findings, and their implications for PGD practices. The advantages and limitations of these studies are analyzed to assess their potential impact on the future development of PGD technologies. Looking forward, we discuss the emerging research directions and challenges, focusing on technological advancements, new application areas, and strategies to overcome existing limitations. In conclusion, this review underscores the pivotal role of single-cell multi-omics in PGD, highlighting its potential to drive the progress of precision medicine and personalized treatment strategies, thereby marking a new era in reproductive genetics and healthcare.

Inverse genetics tracing the differentiation pathway of human chondrocytes

OBJECTIVE

Mammalian somatic cells can be reprogrammed into induced pluripotent stem cells (iPSCs) via the forced expression of Yamanaka reprogramming factors. However, only a limited population of the cells that pass through a particular pathway can metamorphose into iPSCs, while the others do not. This study aimed to clarify the pathways that chondrocytes follow during the reprogramming process.

DESIGN

The fate of human articular chondrocytes under reprogramming was investigated through a time-coursed single-cell transcriptomic analysis, which we termed an inverse genetic approach. The iPS interference technique was also employed to verify that chondrocytes inversely return to pluripotency following the proper differentiation pathway.

RESULTS

We confirmed that human chondrocytes could be converted into cells with an iPSC phenotype. Moreover, it was clarified that a limited population that underwent the silencing of SOX9, a master gene for chondrogenesis, at a specific point during the proper transcriptome transition pathway, could eventually become iPSCs. Interestingly, the other cells, which failed to be reprogrammed, followed a distinct pathway toward cells with a surface zone chondrocyte phenotype. The critical involvement of cellular communication network factors (CCNs) in this process was indicated. The idea that chondrocytes, when reprogrammed into iPSCs, follow the differentiation pathway backward was supported by the successful iPS interference using SOX9.

CONCLUSIONS

This inverse genetic strategy may be useful for seeking candidates for the master genes for the differentiation of various somatic cells. The utility of CCNs in articular cartilage regeneration is also supported.

Spatially resolved epigenome sequencing via Tn5 transposition and deterministic DNA barcoding in tissue

Spatial epigenetic mapping of tissues enables the study of gene regulation programs and cellular functions with the dependency on their local tissue environment. Here we outline a complete procedure for two spatial epigenomic profiling methods: spatially resolved genome-wide profiling of histone modifications using in situ cleavage under targets and tagmentation (CUT&Tag) chemistry (spatial-CUT&Tag) and transposase-accessible chromatin sequencing (spatial-ATAC-sequencing) for chromatin accessibility. Both assays utilize in-tissue Tn5 transposition to recognize genomic DNA loci followed by microfluidic deterministic barcoding to incorporate spatial address codes. Furthermore, these two methods do not necessitate prior knowledge of the transcription or epigenetic markers for a given tissue or cell type but permit genome-wide unbiased profiling pixel-by-pixel at the 10 μm pixel size level and single-base resolution. To support the widespread adaptation of these methods, details are provided in five general steps: (1) sample preparation; (2) Tn5 transposition in spatial-ATAC-sequencing or antibody-controlled pA-Tn5 tagmentation in CUT&Tag; (3) library preparation; (4) next-generation sequencing; and (5) data analysis using our customed pipelines available at: https://github.com/dyxmvp/Spatial_ATAC-seq and https://github.com/dyxmvp/spatial-CUT-Tag . The whole procedure can be completed on four samples in 2-3 days. Familiarity with basic molecular biology and bioinformatics skills with access to a high-performance computing environment are required. A rudimentary understanding of pathology and specimen sectioning, as well as deterministic barcoding in tissue-specific skills (e.g., design of a multiparameter barcode panel and creation of microfluidic devices), are also advantageous. In this protocol, we mainly focus on spatial profiling of tissue region-specific epigenetic landscapes in mouse embryos and mouse brains using spatial-ATAC-sequencing and spatial-CUT&Tag, but these methods can be used for other species with no need for species-specific probe design.

Lack of SMARCB1 expression characterizes a subset of human and murine peripheral T-cell lymphomas

Peripheral T-cell lymphoma, not otherwise specified (PTCL-NOS) is a heterogeneous group of malignancies with poor outcome. Here, we identify a subgroup, PTCL-NOSSMARCB1-, which is characterized by the lack of the SMARCB1 protein and occurs more frequently in young patients. Human and murine PTCL-NOSSMARCB1- show similar DNA methylation profiles, with hypermethylation of T-cell-related genes and hypomethylation of genes involved in myeloid development. Single-cell analyses of human and murine tumors revealed a rich and complex network of interactions between tumor cells and an immunosuppressive and exhausted tumor microenvironment (TME). In a drug screen, we identified histone deacetylase inhibitors (HDACi) as a class of drugs effective against PTCL-NOSSmarcb1-. In vivo treatment of mouse tumors with SAHA, a pan-HDACi, triggered remodeling of the TME, promoting replenishment of lymphoid compartments and reversal of the exhaustion phenotype. These results provide a rationale for further exploration of HDACi combination therapies targeting PTCL-NOSSMARCB1- within the TME.

A primer on single-cell RNA-seq analysis using dendritic cells as a case study

Recent advances in single-cell (sc) transcriptomics have revolutionized our understanding of dendritic cells (DCs), pivotal players of the immune system. ScRNA-sequencing (scRNA-seq) has unraveled a previously unrecognized complexity and heterogeneity of DC subsets, shedding light on their ontogeny and specialized roles. However, navigating the rapid technological progress and computational methods can be daunting for researchers unfamiliar with the field. This review aims to provide immunologists with a comprehensive introduction to sc transcriptomic analysis, offering insights into recent developments in DC biology. Addressing common analytical queries, we guide readers through popular tools and methodologies, supplemented with references to benchmarks and tutorials for in-depth understanding. By examining findings from pioneering studies, we illustrate how computational techniques have expanded our knowledge of DC biology. Through this synthesis, we aim to equip researchers with the necessary tools and knowledge to navigate and leverage scRNA-seq for unraveling the intricacies of DC biology and advancing immunological research.

Immunotherapy drives mesenchymal tumor cell state shift and TME immune response in glioblastoma patients

BACKGROUND

Glioblastoma is a highly aggressive type of brain tumor for which there is no curative treatment available. Immunotherapies have shown limited responses in unselected patients, and there is an urgent need to identify mechanisms of treatment resistance to design novel therapy strategies.

METHODS

Here we investigated the phenotypic and transcriptional dynamics at single-cell resolution during nivolumab immune checkpoint treatment of glioblastoma patients.

RESULTS

We present the integrative paired single-cell RNA-seq analysis of 76 tumor samples from patients in a clinical trial of the PD-1 inhibitor nivolumab and untreated patients. We identify a distinct aggressive phenotypic signature in both tumor cells and the tumor microenvironment in response to nivolumab. Moreover, nivolumab-treatment was associated with an increased transition to mesenchymal stem-like tumor cells, and an increase in TAMs and exhausted and proliferative T cells. We verify and extend our findings in large external glioblastoma dataset (n = 298), develop a latent immune signature and find 18% of primary glioblastoma samples to be latent immune, associated with mesenchymal tumor cell state and TME immune response. Finally, we show that latent immune glioblastoma patients are associated with shorter overall survival following immune checkpoint treatment (P = .0041).

CONCLUSIONS

We find a resistance mechanism signature in one fifth of glioblastoma patients associated with a tumor-cell transition to a more aggressive mesenchymal-like state, increase in TAMs and proliferative and exhausted T cells in response to immunotherapy. These patients may instead benefit from neuro-oncology therapies targeting mesenchymal tumor cells.

Cellular proliferation biases clonal lineage tracing and trajectory inference

MOTIVATION

Lineage tracing and trajectory inference from single-cell RNA-sequencing data hold tremendous potential for uncovering the genetic programs driving development and disease. Single cell datasets are thought to provide an unbiased view on the diverse cellular architecture of tissues. Sampling bias, however, can skew single cell datasets away from the cellular composition they are meant to represent.

RESULTS

We demonstrate a novel form of sampling bias, caused by a statistical phenomenon related to repeated sampling from a growing, heterogeneous population. Relative growth rates of cells influence the probability that they will be sampled in clones observed across multiple time points. We support our probabilistic derivations with a simulation study and an analysis of a real time-course of T-cell development. We find that this bias can impact fate probability predictions, and we explore how to develop trajectory inference methods which are robust to this bias.

AVAILABILITY AND IMPLEMENTATION

Source code for the simulated datasets and to create the figures in this manuscript is freely available in python at https://github.com/rbonhamcarter/simulate-clones. A python implementation of the extension of the LineageOT method is freely available at https://github.com/rbonhamcarter/LineageOT/tree/multi-time-clones.

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

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

Limb connective tissue is organized in a continuum of promiscuous fibroblast identities during development

Connective tissue (CT), which includes tendon and muscle CT, plays critical roles in development, in particular as positional cue provider. Nonetheless, our understanding of fibroblast developmental programs is hampered because fibroblasts are highly heterogeneous and poorly characterized. Combining single-cell RNA-sequencing-based strategies including trajectory inference and in situ hybridization analyses, we address the diversity of fibroblasts and their developmental trajectories during chicken limb fetal development. We show that fibroblasts switch from a positional information to a lineage diversification program at the fetal period onset. Muscle CT and tendon are composed of several fibroblast populations that emerge asynchronously. Once the final muscle pattern is set, transcriptionally close populations are found in neighboring locations in limbs, prefiguring the adult fibroblast layers. We propose that the limb CT is organized in a continuum of promiscuous fibroblast identities, allowing for the robust and efficient connection of muscle to bone and skin.

Molecular heterogeneity of quiescent melanocyte stem cells revealed by single-cell RNA-sequencing

Melanocyte stem cells (McSCs) of the hair follicle are a rare cell population within the skin and are notably underrepresented in whole-skin, single-cell RNA sequencing (scRNA-seq) datasets. Using a cell enrichment strategy to isolate KIT+/CD45- cells from the telogen skin of adult female C57BL/6J mice, we evaluated the transcriptional landscape of quiescent McSCs (qMcSCs) at high resolution. Through this evaluation, we confirmed existing molecular signatures for qMcCS subpopulations (e.g., Kit+, Cd34+/-, Plp1+, Cd274+/-, Thy1+, Cdh3+/-) and identified novel qMcSC subpopulations, including two that differentially regulate their immune privilege status. Within qMcSC subpopulations, we also predicted melanocyte differentiation potential, neural crest potential, and quiescence depth. Taken together, the results demonstrate that the qMcSC population is heterogeneous and future studies focused on investigating changes in qMcSCs should consider changes in subpopulation composition.

An increment of diversity method for cell state trajectory inference of time-series scRNA-seq data

The increasing emergence of the time-series single-cell RNA sequencing (scRNA-seq) data, inferring developmental trajectory by connecting transcriptome similar cell states (i.e., cell types or clusters) has become a major challenge. Most existing computational methods are designed for individual cells and do not take into account the available time series information. We present IDTI based on the Increment of Diversity for Trajectory Inference, which combines time series information and the minimum increment of diversity method to infer cell state trajectory of time-series scRNA-seq data. We apply IDTI to simulated and three real diverse tissue development datasets, and compare it with six other commonly used trajectory inference methods in terms of topology similarity and branching accuracy. The results have shown that the IDTI method accurately constructs the cell state trajectory without the requirement of starting cells. In the performance test, we further demonstrate that IDTI has the advantages of high accuracy and strong robustness.

Modeling blood-brain barrier formation and cerebral cavernous malformations in human PSC-derived organoids

The human blood-brain barrier (hBBB) is a highly specialized structure that regulates passage across blood and central nervous system (CNS) compartments. Despite its critical physiological role, there are no reliable in vitro models that can mimic hBBB development and function. Here, we constructed hBBB assembloids from brain and blood vessel organoids derived from human pluripotent stem cells. We validated the acquisition of blood-brain barrier (BBB)-specific molecular, cellular, transcriptomic, and functional characteristics and uncovered an extensive neuro-vascular crosstalk with a spatial pattern within hBBB assembloids. When we used patient-derived hBBB assembloids to model cerebral cavernous malformations (CCMs), we found that these assembloids recapitulated the cavernoma anatomy and BBB breakdown observed in patients. Upon comparison of phenotypes and transcriptome between patient-derived hBBB assembloids and primary human cavernoma tissues, we uncovered CCM-related molecular and cellular alterations. Taken together, we report hBBB assembloids that mimic the core properties of the hBBB and identify a potentially underlying cause of CCMs.

Esm-1 mediates transcriptional polarization associated with diabetic kidney disease

Esm-1, endothelial cell-specific molecule-1, is a susceptibility gene for diabetic kidney disease (DKD) and is a secreted proteoglycan, with notable expression in kidney, which attenuates inflammation and albuminuria. However, little is known about Esm1 expression in mature tissues in the presence or absence of diabetes. We utilized publicly available single-cell RNA sequencing data to characterize Esm1 expression in 27,786 renal endothelial cells (RECs) obtained from three mouse and four human databases. We validated our findings using bulk transcriptome data from 20 healthy subjects and 41 patients with DKD and using RNAscope. In both mice and humans, Esm1 is expressed in a subset of all REC types and represents a minority of glomerular RECs. In patients, Esm1(+) cells exhibit conserved enrichment for blood vessel development genes. With diabetes, these cells are fewer in number and shift expression toward chemotaxis pathways. Esm1 correlates with a majority of genes within these pathways, delineating a glomerular transcriptional polarization reflected by the magnitude of Esm1 deficiency. Diabetes correlates with lower Esm1 expression and with changes in the functional characterization of Esm1(+) cells. Thus, Esm1 appears to be a marker for glomerular transcriptional polarization in DKD.NEW & NOTEWORTHY Esm-1 is primarily expressed in glomerular endothelium in humans. Cells expressing Esm1 exhibit a high degree of conservation in the enrichment of genes related to blood vessel development. In the context of diabetes, these cells are reduced in number and show a significant transcriptional shift toward the chemotaxis pathway. In diabetes, there is a transcriptional polarization in the glomerulus that is reflected by the degree of Esm1 deficiency.

Branching topology of the human embryo transcriptome revealed by Entropy Sort Feature Weighting

Analysis of single cell transcriptomics (scRNA-seq) data is typically performed after subsetting to highly variable genes (HVGs). Here, we show that Entropy Sorting provides an alternative mathematical framework for feature selection. On synthetic datasets, continuous Entropy Sort Feature Weighting (cESFW) outperforms HVG selection in distinguishing cell-state-specific genes. We apply cESFW to six merged scRNA-seq datasets spanning human early embryo development. Without smoothing or augmenting the raw counts matrices, cESFW generates a high-resolution embedding displaying coherent developmental progression from eight-cell to post-implantation stages and delineating 15 distinct cell states. The embedding highlights sequential lineage decisions during blastocyst development, while unsupervised clustering identifies branch point populations obscured in previous analyses. The first branching region, where morula cells become specified for inner cell mass or trophectoderm, includes cells previously asserted to lack a developmental trajectory. We quantify the relatedness of different pluripotent stem cell cultures to distinct embryo cell types and identify marker genes of naïve and primed pluripotency. Finally, by revealing genes with dynamic lineage-specific expression, we provide markers for staging progression from morula to blastocyst.

A review of the role of epigenetic studies for intramuscular fat deposition in beef cattle

Intramuscular fat (IMF) deposition profoundly influences meat quality and economic value in beef cattle production. Meanwhile, contemporary developments in epigenetics have opened new outlooks for understanding the molecular basics of IMF regulation, and it has become a key area of research for world scholars. Therefore, the aim of this paper was to provide insight and synthesis into the intricate relationship between epigenetic mechanisms and IMF deposition in beef cattle. The methodology involves a thorough analysis of existing literature, including pertinent books, academic journals, and online resources, to provide a comprehensive overview of the role of epigenetic studies in IMF deposition in beef cattle. This review summarizes the contemporary studies in epigenetic mechanisms in IMF regulation, high-resolution epigenomic mapping, single-cell epigenomics, multi-omics integration, epigenome editing approaches, longitudinal studies in cattle growth, environmental epigenetics, machine learning in epigenetics, ethical and regulatory considerations, and translation to industry practices from perspectives of IMF deposition in beef cattle. Moreover, this paper highlights DNA methylation, histone modifications, acetylation, phosphorylation, ubiquitylation, non-coding RNAs, DNA hydroxymethylation, epigenetic readers, writers, and erasers, chromatin immunoprecipitation followed by sequencing, whole genome bisulfite sequencing, epigenome-wide association studies, and their profound impact on the expression of crucial genes governing adipogenesis and lipid metabolism. Nutrition and stress also have significant influences on epigenetic modifications and IMF deposition. The key findings underscore the pivotal role of epigenetic studies in understanding and enhancing IMF deposition in beef cattle, with implications for precision livestock farming and ethical livestock management. In conclusion, this review highlights the crucial significance of epigenetic pathways and environmental factors in affecting IMF deposition in beef cattle, providing insightful information for improving the economics and meat quality of cattle production.

Distribution-Agnostic Deep Learning Enables Accurate Single-Cell Data Recovery and Transcriptional Regulation Interpretation

Single-cell RNA sequencing (scRNA-seq) is a robust method for studying gene expression at the single-cell level, but accurately quantifying genetic material is often hindered by limited mRNA capture, resulting in many missing expression values. Existing imputation methods rely on strict data assumptions, limiting their broader application, and lack reliable supervision, leading to biased signal recovery. To address these challenges, authors developed Bis, a distribution-agnostic deep learning model for accurately recovering missing sing-cell gene expression from multiple platforms. Bis is an optimal transport-based autoencoder model that can capture the intricate distribution of scRNA-seq data while addressing the characteristic sparsity by regularizing the cellular embedding space. Additionally, they propose a module using bulk RNA-seq data to guide reconstruction and ensure expression consistency. Experimental results show Bis outperforms other models across simulated and real datasets, showcasing superiority in various downstream analyses including batch effect removal, clustering, differential expression analysis, and trajectory inference. Moreover, Bis successfully restores gene expression levels in rare cell subsets in a tumor-matched peripheral blood dataset, revealing developmental characteristics of cytokine-induced natural killer cells within a head and neck squamous cell carcinoma microenvironment.

Deciphering cell types by integrating scATAC-seq data with genome sequences

The single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) technology provides insight into gene regulation and epigenetic heterogeneity at single-cell resolution, but cell annotation from scATAC-seq remains challenging due to high dimensionality and extreme sparsity within the data. Existing cell annotation methods mostly focus on the cell peak matrix without fully utilizing the underlying genomic sequence. Here we propose a method, SANGO, for accurate single-cell annotation by integrating genome sequences around the accessibility peaks within scATAC data. The genome sequences of peaks are encoded into low-dimensional embeddings, and then iteratively used to reconstruct the peak statistics of cells through a fully connected network. The learned weights are considered as regulatory modes to represent cells, and utilized to align the query cells and the annotated cells in the reference data through a graph transformer network for cell annotations. SANGO was demonstrated to consistently outperform competing methods on 55 paired scATAC-seq datasets across samples, platforms and tissues. SANGO was also shown to be able to detect unknown tumor cells through attention edge weights learned by the graph transformer. Moreover, from the annotated cells, we found cell-type-specific peaks that provide functional insights/biological signals through expression enrichment analysis, cis-regulatory chromatin interaction analysis and motif enrichment analysis.

Deciphering the spatial landscape and plasticity of immunosuppressive fibroblasts in breast cancer

Although heterogeneity of FAP+ Cancer-Associated Fibroblasts (CAF) has been described in breast cancer, their plasticity and spatial distribution remain poorly understood. Here, we analyze trajectory inference, deconvolute spatial transcriptomics at single-cell level and perform functional assays to generate a high-resolution integrated map of breast cancer (BC), with a focus on inflammatory and myofibroblastic (iCAF/myCAF) FAP+ CAF clusters. We identify 10 spatially-organized FAP+ CAF-related cellular niches, called EcoCellTypes, which are differentially localized within tumors. Consistent with their spatial organization, cancer cells drive the transition of detoxification-associated iCAF (Detox-iCAF) towards immunosuppressive extracellular matrix (ECM)-producing myCAF (ECM-myCAF) via a DPP4- and YAP-dependent mechanism. In turn, ECM-myCAF polarize TREM2+ macrophages, regulatory NK and T cells to induce immunosuppressive EcoCellTypes, while Detox-iCAF are associated with FOLR2+ macrophages in an immuno-protective EcoCellType. FAP+ CAF subpopulations accumulate differently according to the invasive BC status and predict invasive recurrence of ductal carcinoma in situ (DCIS), which could help in identifying low-risk DCIS patients eligible for therapeutic de-escalation.

Tumor Microenvironment Landscapes Supporting EGFR-mutant NSCLC Are Modulated at the Single-cell Interaction Level by Unesbulin Treatment

Lung cancer is the leading cause of cancer deaths. Lethal pulmonary adenocarcinomas (ADC) present with frequent mutations in the EGFR. Genetically engineered murine models of lung cancer expedited comprehension of the molecular mechanisms driving tumorigenesis and drug response. Here, we systematically analyzed the evolution of tumor heterogeneity in the context of dynamic interactions occurring with the intermingled tumor microenvironment (TME) by high-resolution transcriptomics. Our effort identified vulnerable tumor-specific epithelial cells, as well as their cross-talk with niche components (endothelial cells, fibroblasts, and tumor-infiltrating immune cells), whose symbiotic interface shapes tumor aggressiveness and is almost completely abolished by treatment with Unesbulin, a tubulin binding agent that reduces B cell-specific Moloney murine leukemia virus integration site 1 (BMI-1) activity. Simultaneous magnetic resonance imaging (MRI) analysis demonstrated decreased tumor growth, setting the stage for future investigations into the potential of novel therapeutic strategies for EGFR-mutant ADCs.

SIGNIFICANCE

Targeting the TME is an attractive strategy for treatment of solid tumors. Here we revealed how EGFR-mutant landscapes are affected at the single-cell resolution level during Unesbulin treatment. This novel drug, by targeting cancer cells and their interactions with crucial TME components, could be envisioned for future therapeutic advancements.

Statistical method scDEED for detecting dubious 2D single-cell embeddings and optimizing t-SNE and UMAP hyperparameters

Two-dimensional (2D) embedding methods are crucial for single-cell data visualization. Popular methods such as t-distributed stochastic neighbor embedding (t-SNE) and uniform manifold approximation and projection (UMAP) are commonly used for visualizing cell clusters; however, it is well known that t-SNE and UMAP's 2D embeddings might not reliably inform the similarities among cell clusters. Motivated by this challenge, we present a statistical method, scDEED, for detecting dubious cell embeddings output by a 2D-embedding method. By calculating a reliability score for every cell embedding based on the similarity between the cell's 2D-embedding neighbors and pre-embedding neighbors, scDEED identifies the cell embeddings with low reliability scores as dubious and those with high reliability scores as trustworthy. Moreover, by minimizing the number of dubious cell embeddings, scDEED provides intuitive guidance for optimizing the hyperparameters of an embedding method. We show the effectiveness of scDEED on multiple datasets for detecting dubious cell embeddings and optimizing the hyperparameters of t-SNE and UMAP.

Human archetypal pluripotent stem cells differentiate into trophoblast stem cells via endogenous BMP5/7 induction without transitioning through naive state

Primary human trophoblast stem cells (TSCs) and TSCs derived from human pluripotent stem cells (hPSCs) can potentially model placental processes in vitro. Yet, the pluripotent states and factors involved in the differentiation of hPSCs to TSCs remain poorly understood. In this study, we demonstrate that the primed pluripotent state can generate TSCs by activating pathways such as Epidermal Growth Factor (EGF) and Wingless-related integration site (WNT), and by suppressing tumor growth factor beta (TGFβ), histone deacetylases (HDAC), and Rho-associated protein kinase (ROCK) signaling pathways, all without the addition of exogenous Bone morphogenetic protein 4 (BMP4)-a condition we refer to as the TS condition. We characterized this process using temporal single-cell RNA sequencing to compare TS conditions with differentiation protocols involving BMP4 activation alone or BMP4 activation in conjunction with WNT inhibition. The TS condition consistently produced a stable, proliferative cell type that closely mimics first-trimester placental cytotrophoblasts, marked by the activation of endogenous retroviral genes and the absence of amnion expression. This was observed across multiple cell lines, including various primed induced pluripotent stem cell (iPSC) and embryonic stem cell (ESC) lines. Primed-derived TSCs can proliferate for over 30 passages and further specify into multinucleated syncytiotrophoblasts and extravillous trophoblast cells. Our research establishes that the differentiation of primed hPSCs to TSC under TS conditions triggers the induction of TMSB4X, BMP5/7, GATA3, and TFAP2A without progressing through a naive state. These findings propose that the primed hPSC state is part of a continuum of potency with the capacity to differentiate into TSCs through multiple routes.

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

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

Single-Nucleus Analysis Reveals Tumor Heterogeneity of Aldosterone-Producing Adenoma

BACKGROUND

Recent advances in omics techniques have allowed detailed genetic characterization of aldosterone-producing adenoma (APA). The pathogenesis of APA is characterized by tumorigenesis-associated aldosterone synthesis. The pathophysiological intricacies of APAs have not yet been elucidated at the level of individual cells. Therefore, a single-cell level analysis is speculated to be valuable in studying the differentiation process of APA.

METHODS

We conducted single-nucleus RNA sequencing of APAs with KCNJ5 mutation and nonfunctional adenomas obtained from 3 and 2 patients, respectively.

RESULTS

The single-nucleus RNA sequencing revealed the intratumoral heterogeneity of APA and identified cell populations consisting of a shared cluster of nonfunctional adenoma and APA. In addition, we extracted 2 cell fates in APA and obtained a cell population specialized in aldosterone synthesis. Genes related to ribosomes and neurodegenerative diseases were upregulated in 1 of these fates, whereas those related to the regulation of glycolysis were upregulated in the other fate. Furthermore, the total RNA reads in the nucleus were higher in hormonally activated clusters, indicating a marked activation of transcription per cell.

CONCLUSIONS

The single-nucleus RNA sequencing revealed intratumoral heterogeneity of APA with KCNJ5 mutation. The observation of 2 cell fates in KCNJ5-mutated APAs provides the postulation that a heterogeneous process of cellular differentiation was implicated in the pathophysiological mechanisms underlying APA tumors.

TDP-43-stratified single-cell proteomics of postmortem human spinal motor neurons reveals protein dynamics in amyotrophic lateral sclerosis

A limitation of conventional bulk-tissue proteome studies in amyotrophic lateral sclerosis (ALS) is the confounding of motor neuron (MN) signals by admixed non-MN proteins. Here, we leverage laser capture microdissection and nanoPOTS single-cell mass spectrometry-based proteomics to query changes in protein expression in single MNs from postmortem ALS and control tissues. In a follow-up analysis, we examine the impact of stratification of MNs based on cytoplasmic transactive response DNA-binding protein 43 (TDP-43)+ inclusion pathology on the profiles of 2,238 proteins. We report extensive overlap in differentially abundant proteins identified in ALS MNs with or without overt TDP-43 pathology, suggesting early and sustained dysregulation of cellular respiration, mRNA splicing, translation, and vesicular transport in ALS. Together, these data provide insights into proteome-level changes associated with TDP-43 proteinopathy and begin to demonstrate the utility of pathology-stratified trace sample proteomics for understanding single-cell protein dynamics in human neurologic diseases.

Lymph node and tumor-associated PD-L1+ macrophages antagonize dendritic cell vaccines by suppressing CD8+ T cells

Current immunotherapies provide limited benefits against T cell-depleted tumors, calling for therapeutic innovation. Using multi-omics integration of cancer patient data, we predict a type I interferon (IFN) responseHIGH state of dendritic cell (DC) vaccines, with efficacious clinical impact. However, preclinical DC vaccines recapitulating this state by combining immunogenic cancer cell death with induction of type I IFN responses fail to regress mouse tumors lacking T cell infiltrates. Here, in lymph nodes (LNs), instead of activating CD4+/CD8+ T cells, DCs stimulate immunosuppressive programmed death-ligand 1-positive (PD-L1+) LN-associated macrophages (LAMs). Moreover, DC vaccines also stimulate PD-L1+ tumor-associated macrophages (TAMs). This creates two anatomically distinct niches of PD-L1+ macrophages that suppress CD8+ T cells. Accordingly, a combination of PD-L1 blockade with DC vaccines achieves significant tumor regression by depleting PD-L1+ macrophages, suppressing myeloid inflammation, and de-inhibiting effector/stem-like memory T cells. Importantly, clinical DC vaccines also potentiate T cell-suppressive PD-L1+ TAMs in glioblastoma patients. We propose that a multimodal immunotherapy and vaccination regimen is mandatory to overcome T cell-depleted tumors.

Single-cell RNA sequencing of a new transgenic t(8;21) preleukemia mouse model reveals regulatory networks promoting leukemic transformation

T(8;21)(q22;q22), which generates the AML1-ETO fusion oncoprotein, is a common chromosomal abnormality in acute myeloid leukemia (AML) patients. Despite having favorable prognosis, 40% of patients will relapse, highlighting the need for innovative models and application of the newest technologies to study t(8;21) leukemogenesis. Currently, available AML1-ETO mouse models have limited utility for studying the pre-leukemic stage because AML1-ETO produces mild hematopoietic phenotypes and no leukemic transformation. Conversely, overexpression of a truncated variant, AML1-ETO9a (AE9a), promotes fully penetrant leukemia and is too potent for studying pre-leukemic changes. To overcome these limitations, we devised a germline-transmitted Rosa26 locus AE9a knock-in mouse model that moderately overexpressed AE9a and developed leukemia with long latency and low penetrance. We observed pre-leukemic alterations in AE9a mice, including skewing of progenitors towards granulocyte/monocyte lineages and replating of stem and progenitor cells. Next, we performed single-cell RNA sequencing to identify specific cell populations that contribute to these pre-leukemic phenotypes. We discovered a subset of common myeloid progenitors that have heightened granulocyte/monocyte bias in AE9a mice. We also observed dysregulation of key hematopoietic transcription factor target gene networks, blocking cellular differentiation. Finally, we identified Sox4 activation as a potential contributor to stem cell self-renewal during the pre-leukemic stage.

Genome-wide profiling of transcription factor activity in primary liver cancer using single-cell ATAC sequencing

Primary liver cancer (PLC) consists of two main histological subtypes; hepatocellular carcinoma (HCC) and intrahepatic cholangiocarcinoma (iCCA). The role of transcription factors (TFs) in malignant hepatobiliary lineage commitment between HCC and iCCA remains underexplored. Here, we present genome-wide profiling of transcription regulatory elements of 16 PLC patients using single-cell assay for transposase accessible chromatin sequencing. Single-cell open chromatin profiles reflect the compositional diversity of liver cancer, identifying both malignant and microenvironment component cells. TF motif enrichment levels of 31 TFs strongly discriminate HCC from iCCA tumors. These TFs are members of the nuclear/retinoid receptor, POU, or ETS motif families. POU factors are associated with prognostic features in iCCA. Overall, nuclear receptors, ETS and POU TF motif families delineate transcription regulation between HCC and iCCA tumors, which may be relevant to development and selection of PLC subtype-specific therapeutics.

A new and effective two-step clustering approach for single cell RNA sequencing data

BACKGROUND

The rapid devolvement of single cell RNA sequencing (scRNA-seq) technology leads to huge amounts of scRNA-seq data, which greatly advance the research of many biomedical fields involving tissue heterogeneity, pathogenesis of disease and drug resistance etc. One major task in scRNA-seq data analysis is to cluster cells in terms of their expression characteristics. Up to now, a number of methods have been proposed to infer cell clusters, yet there is still much space to improve their performance.

RESULTS

In this paper, we develop a new two-step clustering approach to effectively cluster scRNA-seq data, which is called TSC - the abbreviation of Two-Step Clustering. Particularly, by dividing all cells into two types: core cells (those possibly lying around the centers of clusters) and non-core cells (those locating in the boundary areas of clusters), we first clusters the core cells by hierarchical clustering (the first step) and then assigns the non-core cells to the corresponding nearest clusters (the second step). Extensive experiments on 12 real scRNA-seq datasets show that TSC outperforms the state of the art methods.

CONCLUSION

TSC is an effective clustering method due to its two-steps clustering strategy, and it is a useful tool for scRNA-seq data analysis.

Reconstructing human brown fat developmental trajectory in vitro

Brown adipocytes (BAs) represent a specialized cell type that is able to uncouple nutrient catabolism from ATP generation to dissipate energy as heat. In humans, the brown fat tissue is composed of discrete depots found throughout the neck and trunk region. BAs originate from a precursor common to skeletal muscle, but their developmental trajectory remains poorly understood. Here, we used single-cell RNA sequencing to characterize the development of interscapular brown fat in mice. Our analysis identified a transient stage of BA differentiation characterized by the expression of the transcription factor GATA6. We show that recapitulating the sequence of signaling cues identified in mice can lead to efficient differentiation of BAs in vitro from human pluripotent stem cells. These precursors can in turn be efficiently converted into functional BAs that can respond to signals mimicking adrenergic stimuli by increasing their metabolism, resulting in heat production.

Single-cell DNA methylation and 3D genome architecture in the human brain

Delineating the gene-regulatory programs underlying complex cell types is fundamental for understanding brain function in health and disease. Here, we comprehensively examined human brain cell epigenomes by probing DNA methylation and chromatin conformation at single-cell resolution in 517 thousand cells (399 thousand neurons and 118 thousand non-neurons) from 46 regions of three adult male brains. We identified 188 cell types and characterized their molecular signatures. Integrative analyses revealed concordant changes in DNA methylation, chromatin accessibility, chromatin organization, and gene expression across cell types, cortical areas, and basal ganglia structures. We further developed single-cell methylation barcodes that reliably predict brain cell types using the methylation status of select genomic sites. This multimodal epigenomic brain cell atlas provides new insights into the complexity of cell-type-specific gene regulation in adult human brains.

Toward the functional interpretation of somatic structural variations: bulk- and single-cell approaches

Structural variants (SVs) are genomic rearrangements that can take many different forms such as copy number alterations, inversions and translocations. During cell development and aging, somatic SVs accumulate in the genome with potentially neutral, deleterious or pathological effects. Generation of somatic SVs is a key mutational process in cancer development and progression. Despite their importance, the detection of somatic SVs is challenging, making them less studied than somatic single-nucleotide variants. In this review, we summarize recent advances in whole-genome sequencing (WGS)-based approaches for detecting somatic SVs at the tissue and single-cell levels and discuss their advantages and limitations. First, we describe the state-of-the-art computational algorithms for somatic SV calling using bulk WGS data and compare the performance of somatic SV detectors in the presence or absence of a matched-normal control. We then discuss the unique features of cutting-edge single-cell-based techniques for analyzing somatic SVs. The advantages and disadvantages of bulk and single-cell approaches are highlighted, along with a discussion of their sensitivity to copy-neutral SVs, usefulness for functional inferences and experimental and computational costs. Finally, computational approaches for linking somatic SVs to their functional readouts, such as those obtained from single-cell transcriptome and epigenome analyses, are illustrated, with a discussion of the promise of these approaches in health and diseases.

scDEED: a statistical method for detecting dubious 2D single-cell embeddings and optimizing t-SNE and UMAP hyperparameters

Two-dimensional (2D) embedding methods are crucial for single-cell data visualization. Popular methods such as t-SNE and UMAP are commonly used for visualizing cell clusters; however, it is well known that t-SNE and UMAP's 2D embedding might not reliably inform the similarities among cell clusters. Motivated by this challenge, we developed a statistical method, scDEED, for detecting dubious cell embeddings output by any 2D-embedding method. By calculating a reliability score for every cell embedding, scDEED identifies the cell embeddings with low reliability scores as dubious and those with high reliability scores as trustworthy. Moreover, by minimizing the number of dubious cell embeddings, scDEED provides intuitive guidance for optimizing the hyperparameters of an embedding method. Applied to multiple scRNA-seq datasets, scDEED demonstrates its effectiveness for detecting dubious cell embeddings and optimizing the hyperparameters of t-SNE and UMAP.

Dictys: dynamic gene regulatory network dissects developmental continuum with single-cell multiomics

Gene regulatory networks (GRNs) are key determinants of cell function and identity and are dynamically rewired during development and disease. Despite decades of advancement, challenges remain in GRN inference, including dynamic rewiring, causal inference, feedback loop modeling and context specificity. To address these challenges, we develop Dictys, a dynamic GRN inference and analysis method that leverages multiomic single-cell assays of chromatin accessibility and gene expression, context-specific transcription factor footprinting, stochastic process network and efficient probabilistic modeling of single-cell RNA-sequencing read counts. Dictys improves GRN reconstruction accuracy and reproducibility and enables the inference and comparative analysis of context-specific and dynamic GRNs across developmental contexts. Dictys' network analyses recover unique insights in human blood and mouse skin development with cell-type-specific and dynamic GRNs. Its dynamic network visualizations enable time-resolved discovery and investigation of developmental driver transcription factors and their regulated targets. Dictys is available as a free, open-source and user-friendly Python package.

LVPT: Lazy Velocity Pseudotime Inference Method

The emergence of RNA velocity has enriched our understanding of the dynamic transcriptional landscape within individual cells. In light of this breakthrough, we embarked on integrating RNA velocity with cellular pseudotime inference, aiming to improve the prediction of cell orders along biological trajectories beyond existing methods. Here, we developed LVPT, a novel method for pseudotime and trajectory inference. LVPT introduces a lazy probability to indicate the probability that the cell stays in the original state and calculates the transition matrix based on RNA velocity to provide the probability and direction of cell differentiation. LVPT shows better and comparable performance of pseudotime inference compared with other existing methods on both simulated datasets with different structures and real datasets. The validation results were consistent with prior knowledge, indicating that LVPT is an accurate and efficient method for pseudotime inference.

Transcriptomic forecasting with neural ordinary differential equations

Single-cell transcriptomics technologies can uncover changes in the molecular states that underlie cellular phenotypes. However, understanding the dynamic cellular processes requires extending from inferring trajectories from snapshots of cellular states to estimating temporal changes in cellular gene expression. To address this challenge, we have developed a neural ordinary differential-equation-based method, RNAForecaster, for predicting gene expression states in single cells for multiple future time steps in an embedding-independent manner. We demonstrate that RNAForecaster can accurately predict future expression states in simulated single-cell transcriptomic data with cellular tracking over time. We then show that by using metabolic labeling single-cell RNA sequencing (scRNA-seq) data from constitutively dividing cells, RNAForecaster accurately recapitulates many of the expected changes in gene expression during progression through the cell cycle over a 3-day period. Thus, RNAForecaster enables short-term estimation of future expression states in biological systems from high-throughput datasets with temporal information.

Omics data integration in computational biology viewed through the prism of machine learning paradigms

Important quantities of biological data can today be acquired to characterize cell types and states, from various sources and using a wide diversity of methods, providing scientists with more and more information to answer challenging biological questions. Unfortunately, working with this amount of data comes at the price of ever-increasing data complexity. This is caused by the multiplication of data types and batch effects, which hinders the joint usage of all available data within common analyses. Data integration describes a set of tasks geared towards embedding several datasets of different origins or modalities into a joint representation that can then be used to carry out downstream analyses. In the last decade, dozens of methods have been proposed to tackle the different facets of the data integration problem, relying on various paradigms. This review introduces the most common data types encountered in computational biology and provides systematic definitions of the data integration problems. We then present how machine learning innovations were leveraged to build effective data integration algorithms, that are widely used today by computational biologists. We discuss the current state of data integration and important pitfalls to consider when working with data integration tools. We eventually detail a set of challenges the field will have to overcome in the coming years.

Deciphering the Biology of Circulating Tumor Cells through Single-Cell RNA Sequencing: Implications for Precision Medicine in Cancer

Circulating tumor cells (CTCs) hold unique biological characteristics that directly involve them in hematogenous dissemination. Studying CTCs systematically is technically challenging due to their extreme rarity and heterogeneity and the lack of specific markers to specify metastasis-initiating CTCs. With cutting-edge technology, single-cell RNA sequencing (scRNA-seq) provides insights into the biology of metastatic processes driven by CTCs. Transcriptomics analysis of single CTCs can decipher tumor heterogeneity and phenotypic plasticity for exploring promising novel therapeutic targets. The integrated approach provides a perspective on the mechanisms underlying tumor development and interrogates CTCs interactions with other blood cell types, particularly those of the immune system. This review aims to comprehensively describe the current study on CTC transcriptomic analysis through scRNA-seq technology. We emphasize the workflow for scRNA-seq analysis of CTCs, including enrichment, single cell isolation, and bioinformatic tools applied for this purpose. Furthermore, we elucidated the translational knowledge from the transcriptomic profile of individual CTCs and the biology of cancer metastasis for developing effective therapeutics through targeting key pathways in CTCs.

simCAS: an embedding-based method for simulating single-cell chromatin accessibility sequencing data

MOTIVATION

Single-cell chromatin accessibility sequencing (scCAS) technology provides an epigenomic perspective to characterize gene regulatory mechanisms at single-cell resolution. With an increasing number of computational methods proposed for analyzing scCAS data, a powerful simulation framework is desirable for evaluation and validation of these methods. However, existing simulators generate synthetic data by sampling reads from real data or mimicking existing cell states, which is inadequate to provide credible ground-truth labels for method evaluation.

RESULTS

We present simCAS, an embedding-based simulator, for generating high-fidelity scCAS data from both cell- and peak-wise embeddings. We demonstrate simCAS outperforms existing simulators in resembling real data and show that simCAS can generate cells of different states with user-defined cell populations and differentiation trajectories. Additionally, simCAS can simulate data from different batches and encode user-specified interactions of chromatin regions in the synthetic data, which provides ground-truth labels more than cell states. We systematically demonstrate that simCAS facilitates the benchmarking of four core tasks in downstream analysis: cell clustering, trajectory inference, data integration, and cis-regulatory interaction inference. We anticipate simCAS will be a reliable and flexible simulator for evaluating the ongoing computational methods applied on scCAS data.

AVAILABILITY AND IMPLEMENTATION

simCAS is freely available at https://github.com/Chen-Li-17/simCAS.