Best practices for single-cell analysis across modalities

Systematic data analysis pipeline for quantitative morphological cell phenotyping

Quantitative morphological phenotyping (QMP) is an image-based method used to capture morphological features at both the cellular and population level. Its interdisciplinary nature, spanning from data collection to result analysis and interpretation, can lead to uncertainties, particularly among those new to this actively growing field. High analytical specificity for a typical QMP is achieved through sophisticated approaches that can leverage subtle cellular morphological changes. Here, we outline a systematic workflow to refine the QMP methodology. For a practical review, we describe the main steps of a typical QMP; in each step, we discuss the available methods, their applications, advantages, and disadvantages, along with the R functions and packages for easy implementation. This review does not cover theoretical backgrounds, but provides several references for interested researchers. It aims to broaden the horizons for future phenome studies and demonstrate how to exploit years of endeavors to achieve more with less.

A scalable approach to topic modelling in single-cell data by approximate pseudobulk projection

Probabilistic topic modelling has become essential in many types of single-cell data analysis. Based on probabilistic topic assignments in each cell, we identify the latent representation of cellular states. A dictionary matrix, consisting of topic-specific gene frequency vectors, provides interpretable bases to be compared with known cell type-specific marker genes and other pathway annotations. However, fitting a topic model on a large number of cells would require heavy computational resources-specialized computing units, computing time and memory. Here, we present a scalable approximation method customized for single-cell RNA-seq data analysis, termed ASAP, short for Annotating a Single-cell data matrix by Approximate Pseudobulk estimation. Our approach is more accurate than existing methods but requires orders of magnitude less computing time, leaving much lower memory consumption. We also show that our approach is widely applicable for atlas-scale data analysis; our method seamlessly integrates single-cell and bulk data in joint analysis, not requiring additional preprocessing or feature selection steps.

Revealing gene function with statistical inference at single-cell resolution

Single-cell and spatial molecular profiling assays have shown large gains in sensitivity, resolution and throughput. Applying these technologies to specimens from human and model organisms promises to comprehensively catalogue cell types, reveal their lineage origins in development and discern their contributions to disease pathogenesis. Moreover, rapidly dropping costs have made well-controlled perturbation experiments and cohort studies widely accessible, illuminating mechanisms that give rise to phenotypes at the scale of the cell, the tissue and the whole organism. Interpreting the coming flood of single-cell data, much of which will be spatially resolved, will place a tremendous burden on existing computational pipelines. However, statistical concepts, models, tools and algorithms can be repurposed to solve problems now arising in genetic and molecular biology studies of development and disease. Here, I review how the questions that recent technological innovations promise to answer can be addressed by the major classes of statistical tools.

Benchmarking computational methods for single-cell chromatin data analysis

BACKGROUND

Single-cell chromatin accessibility assays, such as scATAC-seq, are increasingly employed in individual and joint multi-omic profiling of single cells. As the accumulation of scATAC-seq and multi-omics datasets continue, challenges in analyzing such sparse, noisy, and high-dimensional data become pressing. Specifically, one challenge relates to optimizing the processing of chromatin-level measurements and efficiently extracting information to discern cellular heterogeneity. This is of critical importance, since the identification of cell types is a fundamental step in current single-cell data analysis practices.

RESULTS

We benchmark 8 feature engineering pipelines derived from 5 recent methods to assess their ability to discover and discriminate cell types. By using 10 metrics calculated at the cell embedding, shared nearest neighbor graph, or partition levels, we evaluate the performance of each method at different data processing stages. This comprehensive approach allows us to thoroughly understand the strengths and weaknesses of each method and the influence of parameter selection.

CONCLUSIONS

Our analysis provides guidelines for choosing analysis methods for different datasets. Overall, feature aggregation, SnapATAC, and SnapATAC2 outperform latent semantic indexing-based methods. For datasets with complex cell-type structures, SnapATAC and SnapATAC2 are preferred. With large datasets, SnapATAC2 and ArchR are most scalable.

Designing SERS nanotags for profiling overexpressed surface markers on single cancer cells: A review

This review focuses on the chemical design and the use of Surface-Enhanced Raman Scattering (SERS)-active nanotags for measuring surface markers that can be overexpressed at the surface of single cancer cells. Indeed, providing analytical tools with true single-cell measurements capabilities is capital, especially since cancer research is increasingly leaning toward single-cell analysis, either to guide treatment decisions or to understand complex tumor behaviour including the single-cell heterogeneity and the appearance of treatment resistance. Over the past two decades, SERS nanotags have triggered significant interest in the scientific community owing their advantages over fluorescent tags, mainly because SERS nanotags resist photobleaching and exhibit sharper signal bands, which reduces possible spectral overlap and enables the discrimination between the SERS signals and the autofluorescence background from the sample itself. The extensive efforts invested in harnessing SERS nanotags for biomedical purposes, particularly in cancer research, highlight their potential as the next generation of optical labels for single-cell studies. The review unfolds in two main parts. The first part focuses on the structure of SERS nanotags, detailing their chemical composition and the role of each building block of the tags. The second part explores applications in measuring overexpressed surface markers on single-cells. The latter encompasses studies using single nanotags, multiplexed measurements, quantitative information extraction, monitoring treatment responses, and integrating phenotype measurements with SERS nanotags on single cells isolated from complex biological matrices. This comprehensive review anticipates SERS nanotags to persist as a pivotal technology in advancing single-cell analytical methods, particularly in the context of cancer research and personalized medicine.

Current and future directions in network biology

Summary

Network biology is an interdisciplinary field bridging computational and biological sciences that has proved pivotal in advancing the understanding of cellular functions and diseases across biological systems and scales. Although the field has been around for two decades, it remains nascent. It has witnessed rapid evolution, accompanied by emerging challenges. These stem from various factors, notably the growing complexity and volume of data together with the increased diversity of data types describing different tiers of biological organization. We discuss prevailing research directions in network biology, focusing on molecular/cellular networks but also on other biological network types such as biomedical knowledge graphs, patient similarity networks, brain networks, and social/contact networks relevant to disease spread. In more detail, we highlight areas of inference and comparison of biological networks, multimodal data integration and heterogeneous networks, higher-order network analysis, machine learning on networks, and network-based personalized medicine. Following the overview of recent breakthroughs across these five areas, we offer a perspective on future directions of network biology. Additionally, we discuss scientific communities, educational initiatives, and the importance of fostering diversity within the field. This article establishes a roadmap for an immediate and long-term vision for network biology.

Availability and implementation

Not applicable.

Multisensory gamma stimulation mitigates the effects of demyelination induced by cuprizone in male mice

Demyelination is a common pathological feature in a wide range of diseases, characterized by the loss of myelin sheath and myelin-supporting oligodendrocytes. These losses lead to impaired axonal function, increased vulnerability of axons to damage, and result in significant brain atrophy and neuro-axonal degeneration. Multiple pathomolecular processes contribute to neuroinflammation, oligodendrocyte cell death, and progressive neuronal dysfunction. In this study, we use the cuprizone mouse model of demyelination to investigate long-term non-invasive gamma entrainment using sensory stimulation as a potential therapeutic intervention for promoting myelination and reducing neuroinflammation in male mice. Here, we show that multisensory gamma stimulation mitigates demyelination, promotes oligodendrogenesis, preserves functional integrity and synaptic plasticity, attenuates oligodendrocyte ferroptosis-induced cell death, and reduces brain inflammation. Thus, the protective effects of multisensory gamma stimulation on myelin and anti-neuroinflammatory properties support its potential as a therapeutic approach for demyelinating disorders.

scTab: Scaling cross-tissue single-cell annotation models

Identifying cellular identities is a key use case in single-cell transcriptomics. While machine learning has been leveraged to automate cell annotation predictions for some time, there has been little progress in scaling neural networks to large data sets and in constructing models that generalize well across diverse tissues. Here, we propose scTab, an automated cell type prediction model specific to tabular data, and train it using a novel data augmentation scheme across a large corpus of single-cell RNA-seq observations (22.2 million cells). In this context, we show that cross-tissue annotation requires nonlinear models and that the performance of scTab scales both in terms of training dataset size and model size. Additionally, we show that the proposed data augmentation schema improves model generalization. In summary, we introduce a de novo cell type prediction model for single-cell RNA-seq data that can be trained across a large-scale collection of curated datasets and demonstrate the benefits of using deep learning methods in this paradigm.

Connectivity Network Feature Sharing in Single-Cell RNA Sequencing Data Identifies Rare Cells

Single-cell RNA sequencing is a valuable technique for identifying diverse cell subtypes. A key challenge in this process is that the detection of rare cells is often missed by conventional methods due to low abundance and subtle features of these cells. To overcome this, we developed SCLCNF (Local Connectivity Network Feature Sharing in Single-Cell RNA sequencing), a novel approach that identifies rare cells by analyzing features uniquely expressed in these cells. SCLCNF creates a cellular connectivity network, considering how each cell relates to its neighbors. This network helps to pinpoint coexpression patterns unique to rare cells, utilizing a rarity score to confirm their presence. Our method performs better in detecting rare cells than existing techniques, offering enhanced robustness. It has proven to be effective in human gastrula data sets for accurately pinpointing rare cells, and in sepsis data sets where it uncovers previously unidentified rare cell populations.

Functional specialisation and coordination of myonuclei

Myofibres serve as the functional unit for locomotion, with the sarcomere as fundamental subunit. Running the entire length of this structure are hundreds of myonuclei, located at the periphery of the myofibre, juxtaposed to the plasma membrane. Myonuclear specialisation and clustering at the centre and ends of the fibre are known to be essential for muscle contraction, yet the molecular basis of this regionalisation has remained unclear. While the 'myonuclear domain hypothesis' helped explain how myonuclei can independently govern large cytoplasmic territories, novel technologies have provided granularity on the diverse transcriptional programs running simultaneously within the syncytia and added a new perspective on how myonuclei communicate. Building upon this, we explore the critical cellular and molecular sources of transcriptional and functional heterogeneity within myofibres, discussing the impact of intrinsic and extrinsic factors on myonuclear programs. This knowledge provides new insights for understanding muscle development, repair, and disease, but also opens avenues for the development of novel and precise therapeutic approaches.

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.

Transformers in single-cell omics: a review and new perspectives

Recent efforts to construct reference maps of cellular phenotypes have expanded the volume and diversity of single-cell omics data, providing an unprecedented resource for studying cell properties. Despite the availability of rich datasets and their continued growth, current single-cell models are unable to fully capitalize on the information they contain. Transformers have become the architecture of choice for foundation models in other domains owing to their ability to generalize to heterogeneous, large-scale datasets. Thus, the question arises of whether transformers could set off a similar shift in the field of single-cell modeling. Here we first describe the transformer architecture and its single-cell adaptations and then present a comprehensive review of the existing applications of transformers in single-cell analysis and critically discuss their future potential for single-cell biology. By studying limitations and technical challenges, we aim to provide a structured outlook for future research directions at the intersection of machine learning and single-cell biology.

Pleiotropy, a feature or a bug? Toward co-ordinating plant growth, development, and environmental responses through engineering plant hormone signaling

The advent of gene editing technologies such as CRISPR has simplified co-ordinating trait development. However, identifying candidate genes remains a challenge due to complex gene networks and pathways. These networks exhibit pleiotropy, complicating the determination of specific gene and pathway functions. In this review, we explore how systems biology and single-cell sequencing technologies can aid in identifying candidate genes for co-ordinating specifics of plant growth and development within specific temporal and tissue contexts. Exploring sequence-function space of these candidate genes and pathway modules with synthetic biology allows us to test hypotheses and define genotype-phenotype relationships through reductionist approaches. Collectively, these techniques hold the potential to advance breeding and genetic engineering strategies while also addressing genetic diversity issues critical for adaptation and trait development.

Cell Electrokinetic Fingerprint: A Novel Approach Based on Optically Induced Dielectrophoresis (ODEP) for In-Flow Identification of Single Cells

A novel optically induced dielectrophoresis (ODEP) system that can operate under flow conditions is designed for automatic trapping of cells and subsequent induction of 2D multi-frequency cell trajectories. Like in a "ping-pong" match, two virtual electrode barriers operate in an alternate mode with varying frequencies of the input voltage. The so-derived cell motions are characterized via time-lapse microscopy, cell tracking, and state-of-the-art machine learning algorithms, like the wavelet scattering transform (WST). As a cell-electrokinetic fingerprint, the dynamic of variation of the cell displacements happening, over time, is quantified in response to different frequency values of the induced electric field. When tested on two biological scenarios in the cancer domain, the proposed approach discriminates cellular dielectric phenotypes obtained, respectively, at different early phases of drug-induced apoptosis in prostate cancer (PC3) cells and for differential expression of the lectine-like oxidized low-density lipoprotein receptor-1 (LOX-1) transcript levels in human colorectal adenocarcinoma (DLD-1) cells. The results demonstrate increased discrimination of the proposed system and pose an additional basis for making ODEP-based assays addressing cancer heterogeneity for precision medicine and pharmacological research.

Challenges of therapeutic applications and regenerative capacities of urine based stem cells in oral, and maxillofacial reconstruction

Urine-derived stem cells (USCs) have gained the attention of researchers in the biomedical field in the past few years . Regarding the several varieties of cells that have been used for this purpose, USCs have demonstrated mesenchymal stem cell-like properties, such as differentiation and immunomodulation. Furthermore, they could be differentiated into several lineages. This is very interesting for regenerative techniques based on cell therapy. This review will embark on describing their separation, and profiling. We will specifically describe the USCs characteristics, in addition to their differentiation potential. Then, we will introduce and explore the primary uses of USCs. These involve thier utilization as a platform to produce stem cells, however, we shall concentrate on the utilization of USCs for therapeutic, and regenerative orofacial applications, providing an in-depth evaluation of this purpose. The final portion will address the limitations and challenges of their implementation in regenerative dentistry.

scCross: a deep generative model for unifying single-cell multi-omics with seamless integration, cross-modal generation, and in silico exploration

Single-cell multi-omics data reveal complex cellular states, providing significant insights into cellular dynamics and disease. Yet, integration of multi-omics data presents challenges. Some modalities have not reached the robustness or clarity of established transcriptomics. Coupled with data scarcity for less established modalities and integration intricacies, these challenges limit our ability to maximize single-cell omics benefits. We introduce scCross, a tool leveraging variational autoencoders, generative adversarial networks, and the mutual nearest neighbors (MNN) technique for modality alignment. By enabling single-cell cross-modal data generation, multi-omics data simulation, and in silico cellular perturbations, scCross enhances the utility of single-cell multi-omics studies.

Unveiling hidden connections in omics data via pyPARAGON: an integrative hybrid approach for disease network construction

Network inference or reconstruction algorithms play an integral role in successfully analyzing and identifying causal relationships between omics hits for detecting dysregulated and altered signaling components in various contexts, encompassing disease states and drug perturbations. However, accurate representation of signaling networks and identification of context-specific interactions within sparse omics datasets in complex interactomes pose significant challenges in integrative approaches. To address these challenges, we present pyPARAGON (PAgeRAnk-flux on Graphlet-guided network for multi-Omic data integratioN), a novel tool that combines network propagation with graphlets. pyPARAGON enhances accuracy and minimizes the inclusion of nonspecific interactions in signaling networks by utilizing network rather than relying on pairwise connections among proteins. Through comprehensive evaluations on benchmark signaling pathways, we demonstrate that pyPARAGON outperforms state-of-the-art approaches in node propagation and edge inference. Furthermore, pyPARAGON exhibits promising performance in discovering cancer driver networks. Notably, we demonstrate its utility in network-based stratification of patient tumors by integrating phosphoproteomic data from 105 breast cancer tumors with the interactome and demonstrating tumor-specific signaling pathways. Overall, pyPARAGON is a novel tool for analyzing and integrating multi-omic data in the context of signaling networks. pyPARAGON is available at https://github.com/netlab-ku/pyPARAGON.

Single-cell analysis reveals transcriptional dynamics in healthy primary parathyroid tissue

Studies on human parathyroids are generally limited to hyperfunctioning glands owing to the difficulty in obtaining normal human tissue. We therefore obtained non-human primate (NHP) parathyroids to provide a suitable alternative for sequencing that would bear a close semblance to human organs. Single-cell RNA expression analysis of parathyroids from four healthy adult M. mulatta reveals a continuous trajectory of epithelial cell states. Pseudotime analysis based on transcriptomic signatures suggests a progression from GCM2 hi progenitors to mature parathyroid hormone (PTH)-expressing epithelial cells with increasing core mitochondrial transcript abundance along pseudotime. We sequenced, as a comparator, four histologically characterized hyperfunctioning human parathyroids with varying oxyphil and chief cell abundance and leveraged advanced computational techniques to highlight similarities and differences from non-human primate parathyroid expression dynamics. Predicted cell-cell communication analysis reveals abundant endothelial cell interactions in the parathyroid cell microenvironment in both human and NHP parathyroid glands. We show abundant RARRES2 transcripts in both human adenoma and normal primate parathyroid cells and use coimmunostaining to reveal high levels of RARRES2 protein (also known as chemerin) in PTH-expressing cells, which could indicate that RARRES2 plays an unrecognized role in parathyroid endocrine function. The data obtained are the first single-cell RNA transcriptome to characterize nondiseased parathyroid cell signatures and to show a transcriptomic progression of cell states within normal parathyroid glands, which can be used to better understand parathyroid cell biology.

CUT&Tag applied to zebrafish adult tail fins reveals a return of embryonic H3K4me3 patterns during regeneration

Regenerative potential is governed by a complex process of transcriptional reprogramming, involving chromatin reorganization and dynamics in transcription factor binding patterns throughout the genome. The degree to which chromatin and epigenetic changes contribute to this process remains only partially understood. Here we provide a modified CUT&Tag protocol suitable for improved characterization and interrogation of changes in chromatin modifications during adult fin regeneration in zebrafish. Our protocol generates data that recapitulates results from previously published ChIP-Seq methods, requires far fewer cells as input, and significantly improves signal to noise ratios. We deliver high-resolution enrichment maps for H3K4me3 of uninjured and regenerating fin tissues. During regeneration, we find that H3K4me3 levels increase over gene promoters which become transcriptionally active and genes which lose H3K4me3 become silenced. Interestingly, these reprogramming events recapitulate the H3K4me3 patterns observed in developing fin folds of 24-h old zebrafish embryos. Our results indicate that changes in genomic H3K4me3 patterns during fin regeneration occur in a manner consistent with reactivation of developmental programs, demonstrating CUT&Tag to be an effective tool for profiling chromatin landscapes in regenerating tissues.

Trophoblast Side-Population Markers are Dysregulated in Preeclampsia and Fetal Growth Restriction

Dysregulated progenitor cell populations may contribute to poor placental development and placental insufficiency pathogenesis. Side-population cells possess progenitor properties. Recent human trophoblast side-population isolation identified enrichment of 8 specific genes (CXCL8, ELL2, GATA6, HK2, HLA-DPB1, INTS6, SERPINE3 and UPP1) (Gamage et al. 2020, Stem Cell Rev Rep). We characterised these trophoblast side-population markers in human placenta and in placental insufficiency disorders: preeclampsia and fetal growth restriction (FGR). Trophoblast side-population markers localised to mononuclear trophoblasts lining the placental villous basement membrane in preterm control, preeclamptic and FGR placental sections (n = 3, panel of 3 markers/serial section). Analysis of single-cell transcriptomics of an organoid human trophoblast stem cell (hTSC) to extravillous trophoblast (EVT) differentiation model (Shannon et al. 2022, Development) identified that all side-population genes were enriched in mononuclear trophoblast and trophoblasts committed to differentiation under hTSC culture conditions. In vitro validation via 96 h time course hTSC differentiation to EVTs or syncytiotrophoblasts (n = 5) demonstrated ELL2 and HK2 increased with differentiation (p < 0.0024, p < 0.0039 respectively). CXCL8 and HLA-DPB1 were downregulated (p < 0.030, p < 0.011 respectively). GATA6 and INTS6 increased with EVT differentiation only, and UPP1 reduced with syncytialisation. SERPINE3 was undetectable. Trophoblast side-population marker mRNA was measured in human placentas (< 34-weeks' gestation; n = 78 preeclampsia, n = 30 FGR, and n = 18 gestation-matched controls). ELL2, HK2 and CXCL8 were elevated in preeclamptic (p = 0.0006, p < 0.0001, p = 0.0335 respectively) and FGR placentas (p = 0.0065, p < 0.0001, p = 0.0001 respectively) versus controls. Placental GATA6 was reduced in pregnancies with preeclampsia and FGR (p = 0.0014, p = 0.0146 respectively). Placental INTS6 was reduced with FGR only (p < 0.0001). This study identified the localisation of a unique trophoblast subset enriched for side-population markers. Aberrant expression of some side-population markers may indicate disruptions to unique trophoblast subtypes in placental insufficiency.

Immune disease dialogue of chemokine-based cell communications as revealed by single-cell RNA sequencing meta-analysis

Immune-mediated diseases are characterized by aberrant immune responses, posing significant challenges to global health. In both inflammatory and autoimmune diseases, dysregulated immune reactions mediated by tissue-residing immune and non-immune cells precipitate chronic inflammation and tissue damage that is amplified by peripheral immune cell extravasation into the tissue. Chemokine receptors are pivotal in orchestrating immune cell migration, yet deciphering the signaling code across cell types, diseases and tissues remains an open challenge. To delineate disease-specific cell-cell communications involved in immune cell migration, we conducted a meta-analysis of publicly available single-cell RNA sequencing (scRNA-seq) data across diverse immune diseases and tissues. Our comprehensive analysis spanned multiple immune disorders affecting major organs: atopic dermatitis and psoriasis (skin), chronic obstructive pulmonary disease and idiopathic pulmonary fibrosis (lung), ulcerative colitis (colon), IgA nephropathy and lupus nephritis (kidney). By interrogating ligand-receptor (L-R) interactions, alterations in cell proportions, and differential gene expression, we unveiled intricate disease-specific and common immune cell chemoattraction and extravasation patterns. Our findings delineate disease-specific L-R networks and shed light on shared immune responses across tissues and diseases. Insights gleaned from this analysis hold promise for the development of targeted therapeutics aimed at modulating immune cell migration to mitigate inflammation and tissue damage. This nuanced understanding of immune cell dynamics at the single-cell resolution opens avenues for precision medicine in immune disease management.

Computational estimation of clonal diversity in autoimmunity

Diversity is the cornerstone of the adaptive immune system, crucial for its effectiveness against constantly evolving pathogens that pose threats to higher vertebrates. Accurately measuring and interpreting this diversity presents challenges for immunologists, as changes in diversity and clonotype composition can tip the balance between protective immunity and autoimmunity. In this review, we present the current methods commonly used to measure diversity from single-cell T-cell receptor and B-cell receptor sequencing. We also discuss two case studies where single-cell sequencing and diversity estimations have led to breakthroughs in autoimmune disease discovery and therapeutic innovation, and reflect upon the necessity and importance of accurately defining and measuring lymphocyte diversity in these contexts.

Single-cell RNA-sequencing identifies unique cell-specific gene expression profiles in high-grade cardiac allograft vasculopathy

Background

Cardiac allograft vasculopathy (CAV), a diffuse thickening of the intima of the coronary arteries and microvasculature, is the leading cause of late graft failure and mortality after heart transplantation (HT). Diagnosis involves invasive coronary angiography, which carries substantial risk, and minimally-invasive approaches to CAV diagnosis are urgently needed. Using single-cell RNA-sequencing in peripheral blood mononuclear cells (PBMCs), we sought to identify cell-specific gene expression profiles in CAV.

Methods

Whole blood was collected from 22 HT recipients with angiographically-confirmed CAV and 18 HT recipients without CAV. PBMCs were isolated and subjected to single-cell RNA-sequencing using a 10X Genomics microfluidic platform. Downstream analyses focused on differential expression of genes, cell compositional changes, and T cell receptor repertoire analyses.

Results

Across 40 PBMC samples, we isolated 134,984 cells spanning 8 major clusters and 31 subclusters of cell types. Compositional analyses showed subtle, but significant increases in CD4+ T central memory cells, and CD14+ and CD16+ monocytes in high-grade CAV (CAV-2 and CAV-3) as compared to low-grade or absent CAV. After adjusting for age, gender, and prednisone use, 745 genes were differentially expressed in a cell-specific manner in high-grade CAV. Weighted gene co-expression network analyses showed enrichment for putative pathways involved in inflammation and angiogenesis. There were no significant differences in T cell clonality or diversity with increasing CAV severity.

Conclusions

Unbiased whole transcriptomic analyses at single-cell resolution identify unique, cell-specific gene expression patterns in CAV, suggesting the potential utility of peripheral gene expression biomarkers in diagnosing CAV.

Are We There Yet? Assessing the Readiness of Single-Cell Proteomics to Answer Biological Hypotheses

Single-cell analysis is an active area of research in many fields of biology. Measurements at single-cell resolution allow researchers to study diverse populations without losing biologically meaningful information to sample averages. Many technologies have been used to study single cells, including mass spectrometry-based single-cell proteomics (SCP). SCP has seen a lot of growth over the past couple of years through improvements in data acquisition and analysis, leading to greater proteomic depth. Because method development has been the main focus in SCP, biological applications have been sprinkled in only as proof-of-concept. However, SCP methods now provide significant coverage of the proteome and have been implemented in many laboratories. Thus, a primary question to address in our community is whether the current state of technology is ready for widespread adoption for biological inquiry. In this Perspective, we examine the potential for SCP in three thematic areas of biological investigation: cell annotation, developmental trajectories, and spatial mapping. We identify that the primary limitation of SCP is sample throughput. As proteome depth has been the primary target for method development to date, we advocate for a change in focus to facilitate measuring tens of thousands of single-cell proteomes to enable biological applications beyond proof-of-concept.

Recent Innovations and Technical Advances in High-Throughput Parallel Single-Cell Whole-Genome Sequencing Methods

Single-cell whole-genome sequencing (scWGS) detects cell heterogeneity at the aspect of genomic variations, which are inheritable and play an important role in life processes such as aging and cancer progression. The recent explosive development of high-throughput single-cell sequencing methods has enabled high-performance heterogeneity detection through a vast number of novel strategies. Despite the limitation on total cost, technical advances in high-throughput single-cell whole-genome sequencing methods are made for higher genome coverage, parallel throughput, and level of integration. This review highlights the technical advancements in high-throughput scWGS in the aspects of strategies design, data efficiency, parallel handling platforms, and their applications on human genome. The experimental innovations, remaining challenges, and perspectives are summarized and discussed.

Panpipes: a pipeline for multiomic single-cell and spatial transcriptomic data analysis

Single-cell multiomic analysis of the epigenome, transcriptome, and proteome allows for comprehensive characterization of the molecular circuitry that underpins cell identity and state. However, the holistic interpretation of such datasets presents a challenge given a paucity of approaches for systematic, joint evaluation of different modalities. Here, we present Panpipes, a set of computational workflows designed to automate multimodal single-cell and spatial transcriptomic analyses by incorporating widely-used Python-based tools to perform quality control, preprocessing, integration, clustering, and reference mapping at scale. Panpipes allows reliable and customizable analysis and evaluation of individual and integrated modalities, thereby empowering decision-making before downstream investigations.

GeneSurfer Enables Transcriptome-wide Exploration and Functional Annotation of Gene Co-expression Modules in 3D Spatial Transcriptomics Data

Gene co-expression provides crucial insights into biological functions, however, there is a lack of exploratory analysis tools for localized gene co-expression in large-scale datasets. We present GeneSurfer, an interactive interface designed to explore localized transcriptome-wide gene co-expression patterns in the 3D spatial domain. Key features of GeneSurfer include transcriptome-wide gene filtering and gene clustering based on spatial local co-expression within transcriptomically similar cells, multi-slice 3D rendering of average expression of gene clusters, and on-the-fly Gene Ontology term annotation of co-expressed gene sets. Additionally, GeneSurfer offers multiple linked views for investigating individual genes or gene co-expression in the spatial domain at each exploration stage. Demonstrating its utility with both spatial transcriptomics and single-cell RNA sequencing data from the Allen Brain Cell Atlas, GeneSurfer effectively identifies and annotates localized transcriptome-wide co-expression, providing biological insights and facilitating hypothesis generation and validation.

Stochastic Modeling of Biophysical Responses to Perturbation

Recent advances in high-throughput, multi-condition experiments allow for genome-wide investigation of how perturbations affect transcription and translation in the cell across multiple biological entities or modalities, from chromatin and mRNA information to protein production and spatial morphology. This presents an unprecedented opportunity to unravel how the processes of DNA and RNA regulation direct cell fate determination and disease response. Most methods designed for analyzing large-scale perturbation data focus on the observational outcomes, e.g., expression; however, many potential transcriptional mechanisms, such as transcriptional bursting or splicing dynamics, can underlie these complex and noisy observations. In this analysis, we demonstrate how a stochastic biophysical modeling approach to interpreting high-throughout perturbation data enables deeper investigation of the 'how' behind such molecular measurements. Our approach takes advantage of modalities already present in data produced with current technologies, such as nascent and mature mRNA measurements, to illuminate transcriptional dynamics induced by perturbation, predict kinetic behaviors in new perturbation settings, and uncover novel populations of cells with distinct kinetic responses to perturbation.

SingmiR: a single-cell miRNA alignment and analysis tool

Single-cell RNA sequencing (RNA-seq) has revolutionized our understanding of cell biology, developmental and pathophysiological molecular processes, paving the way toward novel diagnostic and therapeutic approaches. However, most of the gene regulatory processes on the single-cell level are still unknown, including post-transcriptional control conferred by microRNAs (miRNAs). Like the established single-cell gene expression analysis, advanced computational expertise is required to comprehensively process newly emerging single-cell miRNA-seq datasets. A web server providing a workflow tailored for single-cell miRNA-seq data with a self-explanatory interface is currently not available. Here, we present SingmiR, enabling the rapid (pre-)processing and quantification of human miRNAs from noncoding single-cell samples. It performs read trimming for different library preparation protocols, generates automated quality control reports and provides feature-normalized count files. Numerous standard and advanced analyses such as dimension reduction, clustered feature heatmaps, sample correlation heatmaps and differential expression statistics are implemented. We aim to speed up the prototyping pipeline for biologists developing single-cell miRNA-seq protocols on small to medium-sized datasets. SingmiR is freely available to all users without the need for a login at https://www.ccb.uni-saarland.de/singmir.

Trigeminal ganglion neurons are directly activated by influx of CSF solutes in a migraine model

Classical migraine patients experience aura, which is transient neurological deficits associated with cortical spreading depression (CSD), preceding headache attacks. It is not currently understood how a pathological event in cortex can affect peripheral sensory neurons. In this study, we show that cerebrospinal fluid (CSF) flows into the trigeminal ganglion, establishing nonsynaptic signaling between brain and trigeminal cells. After CSD, ~11% of the CSF proteome is altered, with up-regulation of proteins that directly activate receptors in the trigeminal ganglion. CSF collected from animals exposed to CSD activates trigeminal neurons in naïve mice in part by CSF-borne calcitonin gene-related peptide (CGRP). We identify a communication pathway between the central and peripheral nervous system that might explain the relationship between migrainous aura and headache.

Coding and Non-Coding Transcriptomic Landscape of Aortic Complications in Marfan Syndrome

Marfan syndrome (MFS) is a rare congenital disorder of the connective tissue, leading to thoracic aortic aneurysms (TAA) and dissection, among other complications. Currently, the most efficient strategy to prevent life-threatening dissection is preventive surgery. Periodic imaging applying complex techniques is required to monitor TAA progression and to guide the timing of surgical intervention. Thus, there is an acute demand for non-invasive biomarkers for diagnosis and prognosis, as well as for innovative therapeutic targets of MFS. Unraveling the intricate pathomolecular mechanisms underlying the syndrome is vital to address these needs. High-throughput platforms are particularly well-suited for this purpose, as they enable the integration of different datasets, such as transcriptomic and epigenetic profiles. In this narrative review, we summarize relevant studies investigating changes in both the coding and non-coding transcriptome and epigenome in MFS-induced TAA. The collective findings highlight the implicated pathways, such as TGF-β signaling, extracellular matrix structure, inflammation, and mitochondrial dysfunction. Potential candidates as biomarkers, such as miR-200c, as well as therapeutic targets emerged, like Tfam, associated with mitochondrial respiration, or miR-632, stimulating endothelial-to-mesenchymal transition. While these discoveries are promising, rigorous and extensive validation in large patient cohorts is indispensable to confirm their clinical relevance and therapeutic potential.

Multi-modal generative modeling for joint analysis of single-cell T cell receptor and gene expression data

Recent advances in single-cell immune profiling have enabled the simultaneous measurement of transcriptome and T cell receptor (TCR) sequences, offering great potential for studying immune responses at the cellular level. However, integrating these diverse modalities across datasets is challenging due to their unique data characteristics and technical variations. Here, to address this, we develop the multimodal generative model mvTCR to fuse modality-specific information across transcriptome and TCR into a shared representation. Our analysis demonstrates the added value of multimodal over unimodal approaches to capture antigen specificity. Notably, we use mvTCR to distinguish T cell subpopulations binding to SARS-CoV-2 antigens from bystander cells. Furthermore, when combined with reference mapping approaches, mvTCR can map newly generated datasets to extensive T cell references, facilitating knowledge transfer. In summary, we envision mvTCR to enable a scalable analysis of multimodal immune profiling data and advance our understanding of immune responses.

Single cell regulatory architecture of human pancreatic islets suggests sex differences in β cell function and the pathogenesis of type 2 diabetes

Type 2 and type 1 diabetes (T2D, T1D) exhibit sex differences in insulin secretion, the mechanisms of which are unknown. We examined sex differences in human pancreatic islets from 52 donors with and without T2D combining single cell RNA-seq (scRNA-seq), single nucleus assay for transposase-accessible chromatin sequencing (snATAC-seq), hormone secretion, and bioenergetics. In nondiabetic (ND) donors, sex differences in islet cells gene accessibility and expression predominantly involved sex chromosomes. Islets from T2D donors exhibited similar sex differences in sex chromosomes differentially expressed genes (DEGs), but also exhibited sex differences in autosomal genes. Comparing β cells from T2D vs. ND donors, gene enrichment of female β cells showed suppression in mitochondrial respiration, while male β cells exhibited suppressed insulin secretion. Thus, although sex differences in gene accessibility and expression of ND β cells predominantly affect sex chromosomes, the transition to T2D reveals sex differences in autosomes highlighting mitochondrial failure in females.

Opportunities and Challenges in Advancing Plant Research with Single-cell Omics

Plants possess diverse cell types and intricate regulatory mechanisms to adapt to the ever-changing environment of nature. Various strategies have been employed to study cell types and their developmental progressions, including single-cell sequencing methods which provide high-dimensional catalogs to address biological concerns. In recent years, single-cell sequencing technologies in transcriptomics, epigenomics, proteomics, metabolomics, and spatial transcriptomics have been increasingly used in plant science to reveal intricate biological relationships at the single-cell level. However, the application of single-cell technologies to plants is more limited due to the challenges posed by cell structure. This review outlines the advancements in single-cell omics technologies, their implications in plant systems, future research applications, and the challenges of single-cell omics in plant systems.

Applications of single-cell multi-omics in liver cancer

Primary liver cancer, more specifically hepatocellular carcinoma (HCC), remains a significant global health problem associated with increasing incidence and mortality. Clinical, biological, and molecular heterogeneity are well-known hallmarks of cancer and HCC is considered one of the most heterogeneous tumour types, displaying substantial inter-patient, intertumoural and intratumoural variability. This heterogeneity plays a pivotal role in hepatocarcinogenesis, metastasis, relapse and drug response or resistance. Unimodal single-cell sequencing techniques have already revolutionised our understanding of the different layers of molecular hierarchy in the tumour microenvironment of HCC. By highlighting the cellular heterogeneity and the intricate interactions among cancer, immune and stromal cells before and during treatment, these techniques have contributed to a deeper comprehension of tumour clonality, hematogenous spreading and the mechanisms of action of immune checkpoint inhibitors. However, major questions remain to be elucidated, with the identification of biomarkers predicting response or resistance to immunotherapy-based regimens representing an important unmet clinical need. Although the application of single-cell multi-omics in liver cancer research has been limited thus far, a revolution of individualised care for patients with HCC will only be possible by integrating various unimodal methods into multi-omics methodologies at the single-cell resolution. In this review, we will highlight the different established single-cell sequencing techniques and explore their biological and clinical impact on liver cancer research, while casting a glance at the future role of multi-omics in this dynamic and rapidly evolving field.

Building and analyzing metacells in single-cell genomics data

The advent of high-throughput single-cell genomics technologies has fundamentally transformed biological sciences. Currently, millions of cells from complex biological tissues can be phenotypically profiled across multiple modalities. The scaling of computational methods to analyze and visualize such data is a constant challenge, and tools need to be regularly updated, if not redesigned, to cope with ever-growing numbers of cells. Over the last few years, metacells have been introduced to reduce the size and complexity of single-cell genomics data while preserving biologically relevant information and improving interpretability. Here, we review recent studies that capitalize on the concept of metacells-and the many variants in nomenclature that have been used. We further outline how and when metacells should (or should not) be used to analyze single-cell genomics data and what should be considered when analyzing such data at the metacell level. To facilitate the exploration of metacells, we provide a comprehensive tutorial on the construction and analysis of metacells from single-cell RNA-seq data ( https://github.com/GfellerLab/MetacellAnalysisTutorial ) as well as a fully integrated pipeline to rapidly build, visualize and evaluate metacells with different methods ( https://github.com/GfellerLab/MetacellAnalysisToolkit ).

CellRank 2: unified fate mapping in multiview single-cell data

Single-cell RNA sequencing allows us to model cellular state dynamics and fate decisions using expression similarity or RNA velocity to reconstruct state-change trajectories; however, trajectory inference does not incorporate valuable time point information or utilize additional modalities, whereas methods that address these different data views cannot be combined or do not scale. Here we present CellRank 2, a versatile and scalable framework to study cellular fate using multiview single-cell data of up to millions of cells in a unified fashion. CellRank 2 consistently recovers terminal states and fate probabilities across data modalities in human hematopoiesis and endodermal development. Our framework also allows combining transitions within and across experimental time points, a feature we use to recover genes promoting medullary thymic epithelial cell formation during pharyngeal endoderm development. Moreover, we enable estimating cell-specific transcription and degradation rates from metabolic-labeling data, which we apply to an intestinal organoid system to delineate differentiation trajectories and pinpoint regulatory strategies.

Unveiling mitochondrial and ribosomal gene deregulation and tumor microenvironment dynamics in acute myeloid leukemia

Acute myeloid leukemia (AML) is a malignant clonal hematopoietic disease with a poor prognosis. Understanding the interaction between leukemic cells and the tumor microenvironment (TME) can help predict the prognosis of leukemia and guide its treatment. Re-analyzing the scRNA-seq data from the CSC and G20 cohorts, using a Python-based pipeline including machine-learning-based scVI-tools, recapitulated the distinct hierarchical structure within the samples of AML patients. Weighted correlation network analysis (WGCNA) was conducted to construct a weighted gene co-expression network and to identify gene modules primarily focusing on hematopoietic stem cells (HSCs), multipotent progenitors (MPPs), and natural killer (NK) cells. The analysis revealed significant deregulation in gene modules associated with aerobic respiration and ribosomal/cytoplasmic translation. Cell-cell communications were elucidated by the CellChat package, revealing an imbalance of activating and inhibitory immune signaling pathways. Interception of genes upregulated in leukemic HSCs & MPPs as well as in NKG2A-high NK cells was used to construct prognostic models. Normal Cox and artificial neural network models based on 10 genes were developed. The study reveals the deregulation of mitochondrial and ribosomal genes in AML patients and suggests the co-occurrence of stimulatory and inhibitory factors in the AML TME.

Inflammatory memory in psoriasis: From remission to recurrence

The routine use of targeted systemic immunomodulatory therapies has transformed outcomes for people with severe psoriasis, with skin clearance (clinical remission) rates up to 60% at 1 year of biologic treatment. However, psoriasis may recur following drug withdrawal, and as a result, patients tend to continue receiving costly treatment indefinitely. Here, we review research into the "inflammatory memory" in resolved psoriasis skin and the potential mechanisms leading to psoriasis recurrence following drug withdrawal. Research has implicated immune cells such as tissue resident memory T cells, Langerhans cells, and dermal dendritic cells, and there is growing interest in keratinocytes and fibroblasts. A better understanding of the interactions between these cell populations, enabled by single cell technologies, will help to elucidate the events underpinning the shift from remission to recurrence. This may inform the development of personalized strategies for sustaining remission while reducing long-term drug burden.

Deep 3D histology powered by tissue clearing, omics and AI

To comprehensively understand tissue and organism physiology and pathophysiology, it is essential to create complete three-dimensional (3D) cellular maps. These maps require structural data, such as the 3D configuration and positioning of tissues and cells, and molecular data on the constitution of each cell, spanning from the DNA sequence to protein expression. While single-cell transcriptomics is illuminating the cellular and molecular diversity across species and tissues, the 3D spatial context of these molecular data is often overlooked. Here, I discuss emerging 3D tissue histology techniques that add the missing third spatial dimension to biomedical research. Through innovations in tissue-clearing chemistry, labeling and volumetric imaging that enhance 3D reconstructions and their synergy with molecular techniques, these technologies will provide detailed blueprints of entire organs or organisms at the cellular level. Machine learning, especially deep learning, will be essential for extracting meaningful insights from the vast data. Further development of integrated structural, molecular and computational methods will unlock the full potential of next-generation 3D histology.

Multi-omics data integration using ratio-based quantitative profiling with Quartet reference materials

Characterization and integration of the genome, epigenome, transcriptome, proteome and metabolome of different datasets is difficult owing to a lack of ground truth. Here we develop and characterize suites of publicly available multi-omics reference materials of matched DNA, RNA, protein and metabolites derived from immortalized cell lines from a family quartet of parents and monozygotic twin daughters. These references provide built-in truth defined by relationships among the family members and the information flow from DNA to RNA to protein. We demonstrate how using a ratio-based profiling approach that scales the absolute feature values of a study sample relative to those of a concurrently measured common reference sample produces reproducible and comparable data suitable for integration across batches, labs, platforms and omics types. Our study identifies reference-free 'absolute' feature quantification as the root cause of irreproducibility in multi-omics measurement and data integration and establishes the advantages of ratio-based multi-omics profiling with common reference materials.

Protocol for quantifying stem-cell-derived cardiomyocyte maturity using transcriptomic entropy score

The inability to quantify cardiomyocyte (CM) maturation remains a significant barrier to evaluating the effects of ongoing efforts to produce adult-like CMs from pluripotent stem cells (PSCs). Here, we present a protocol to quantify stem-cell-derived CM maturity using a single-cell RNA sequencing-based metric "entropy score." We describe steps for generating an entropy score using customized R code. This tool can be used to quantify maturation levels of PSC-CMs and potentially other cell types. For complete details on the use and execution of this protocol, please refer to Kannan et al.1.

Size and Surface Properties of Functionalized Organosilica Particles Impact Cell-Particle Interactions Including Mitochondrial Activity

Understanding of the interactions between macrophages and multifunctional nanoparticles is important for development of novel macrophage-based immunotherapies. Here, we investigated the effects of fluorescent thiol-organosilica particle size and surface properties on cell-particle interactions, including mitochondrial activity, using the mouse macrophage cell line J774A.1. Three different sizes of thiol-organosilica particles (150, 400, and 680 nm in diameter) containing fluorescein (OS/F150, OS/F400, and OS/F680) and particles surface functionalized with polyethylenimine (PEI) (OS/F150PEI, OS/F400PEI, and OS/F680PEI) were prepared. Flow cytometric analysis, time-lapse imaging, and single-cell analysis of particle uptake and mitochondrial activity of J774A.1 cells demonstrated variations in uptake and kinetics depending on the particle size and surface as well as on each individual cell. Cells treated with OS/F150 and OS/F150PEI showed higher uptake and mitochondrial activity than those treated with other particles. The interaction between endosomes and mitochondria was observed using 3D fluorescent imaging and was characterized by the involvement of iron transport into mitochondria by iron-containing proteins adsorbed on the particle surface. Scanning electron microscopy of the cells treated with the particles revealed alterations in cell membrane morphology, depending on particle size and surface. We performed correlative light and electron microscopy combined with time-lapse and 3D imaging to develop an integrated correlation analysis of particle uptake, mitochondrial activity, and cell membrane morphology in single macrophages. These cell-specific characteristics of macrophages against functional particles and their evaluation methods are crucial for understanding the immunological functions of individual macrophages and developing novel immunotherapies.

Defining the Regulatory Logic of Breast Cancer Using Single-Cell Epigenetic and Transcriptome Profiling

Annotation of the cis-regulatory elements that drive transcriptional dysregulation in cancer cells is critical to improving our understanding of tumor biology. Herein, we present a compendium of matched chromatin accessibility (scATAC-seq) and transcriptome (scRNA-seq) profiles at single-cell resolution from human breast tumors and healthy mammary tissues processed immediately following surgical resection. We identify the most likely cell-of-origin for luminal breast tumors and basal breast tumors and then introduce a novel methodology that implements linear mixed-effects models to systematically quantify associations between regions of chromatin accessibility (i.e. regulatory elements) and gene expression in malignant cells versus normal mammary epithelial cells. These data unveil regulatory elements with that switch from silencers of gene expression in normal cells to enhancers of gene expression in cancer cells, leading to the upregulation of clinically relevant oncogenes. To translate the utility of this dataset into tractable models, we generated matched scATAC-seq and scRNA-seq profiles for breast cancer cell lines, revealing, for each subtype, a conserved oncogenic gene expression program between in vitro and in vivo cells. Together, this work highlights the importance of non-coding regulatory mechanisms that underlie oncogenic processes and the ability of single-cell multi-omics to define the regulatory logic of BC cells at single-cell resolution.

Identifying Reproducible Transcription Regulator Coexpression Patterns with Single Cell Transcriptomics

The proliferation of single cell transcriptomics has potentiated our ability to unveil patterns that reflect dynamic cellular processes, rather than cell type compositional effects that emerge from bulk tissue samples. In this study, we leverage a broad collection of single cell RNA-seq data to identify the gene partners whose expression is most coordinated with each human and mouse transcription regulator (TR). We assembled 120 human and 103 mouse scRNA-seq datasets from the literature (>28 million cells), constructing a single cell coexpression network for each. We aimed to understand the consistency of TR coexpression profiles across a broad sampling of biological contexts, rather than examine the preservation of context-specific signals. Our workflow therefore explicitly prioritizes the patterns that are most reproducible across cell types. Towards this goal, we characterize the similarity of each TR's coexpression within and across species. We create single cell coexpression rankings for each TR, demonstrating that this aggregated information recovers literature curated targets on par with ChIP-seq data. We then combine the coexpression and ChIP-seq information to identify candidate regulatory interactions supported across methods and species. Finally, we highlight interactions for the important neural TR ASCL1 to demonstrate how our compiled information can be adopted for community use.

Gene count normalization in single-cell imaging-based spatially resolved transcriptomics

BACKGROUND

Recent advances in imaging-based spatially resolved transcriptomics (im-SRT) technologies now enable high-throughput profiling of targeted genes and their locations in fixed tissues. Normalization of gene expression data is often needed to account for technical factors that may confound underlying biological signals.

RESULTS

Here, we investigate the potential impact of different gene count normalization methods with different targeted gene panels in the analysis and interpretation of im-SRT data. Using different simulated gene panels that overrepresent genes expressed in specific tissue regions or cell types, we demonstrate how normalization methods based on detected gene counts per cell differentially impact normalized gene expression magnitudes in a region- or cell type-specific manner. We show that these normalization-induced effects may reduce the reliability of downstream analyses including differential gene expression, gene fold change, and spatially variable gene analysis, introducing false positive and false negative results when compared to results obtained from gene panels that are more representative of the gene expression of the tissue's component cell types. These effects are not observed with normalization approaches that do not use detected gene counts for gene expression magnitude adjustment, such as with cell volume or cell area normalization.

CONCLUSIONS

We recommend using non-gene count-based normalization approaches when feasible and evaluating gene panel representativeness before using gene count-based normalization methods if necessary. Overall, we caution that the choice of normalization method and gene panel may impact the biological interpretation of the im-SRT data.

Olfactory neuroblastoma mimics molecular heterogeneity and lineage trajectories of small-cell lung cancer

The olfactory epithelium undergoes neuronal regeneration from basal stem cells and is susceptible to olfactory neuroblastoma (ONB), a rare tumor of unclear origins. Employing alterations in Rb1/Trp53/Myc (RPM), we establish a genetically engineered mouse model of high-grade metastatic ONB exhibiting a NEUROD1+ immature neuronal phenotype. We demonstrate that globose basal cells (GBCs) are a permissive cell of origin for ONB and that ONBs exhibit cell fate heterogeneity that mimics normal GBC developmental trajectories. ASCL1 loss in RPM ONB leads to emergence of non-neuronal histopathologies, including a POU2F3+ microvillar-like state. Similar to small-cell lung cancer (SCLC), mouse and human ONBs exhibit mutually exclusive NEUROD1 and POU2F3-like states, an immune-cold tumor microenvironment, intratumoral cell fate heterogeneity comprising neuronal and non-neuronal lineages, and cell fate plasticity-evidenced by barcode-based lineage tracing and single-cell transcriptomics. Collectively, our findings highlight conserved similarities between ONB and neuroendocrine tumors with significant implications for ONB classification and treatment.

SpotSweeper: spatially-aware quality control for spatial transcriptomics

Quality control (QC) is a crucial step to ensure the reliability and accuracy of the data obtained from RNA sequencing experiments, including spatially-resolved transcriptomics (SRT). Existing QC approaches for SRT that have been adopted from single-nucleus RNA sequencing (snRNA-seq) methods are confounded by spatial biology and are inappropriate for SRT data. In addition, no methods currently exist for identifying histological tissue artifacts unique to SRT. Here, we introduce SpotSweeper, spatially-aware QC methods for identifying local outliers and regional artifacts in SRT. SpotSweeper evaluates the quality of individual spots relative to their local neighborhood, thus minimizing bias due to biological heterogeneity, and uses multiscale methods to detect regional artifacts. Using SpotSweeper on publicly available data, we identified a consistent set of Visium barcodes/spots as systematically low quality and demonstrate that SpotSweeper accurately identifies two distinct types of regional artifacts, resulting in improved downstream clustering and marker gene detection for spatial domains.

Revisiting Cardiac Biology in the Era of Single Cell and Spatial Omics

Throughout our lifetime, each beat of the heart requires the coordinated action of multiple cardiac cell types. Understanding cardiac cell biology, its intricate microenvironments, and the mechanisms that govern their function in health and disease are crucial to designing novel therapeutical and behavioral interventions. Recent advances in single-cell and spatial omics technologies have significantly propelled this understanding, offering novel insights into the cellular diversity and function and the complex interactions of cardiac tissue. This review provides a comprehensive overview of the cellular landscape of the heart, bridging the gap between suspension-based and emerging in situ approaches, focusing on the experimental and computational challenges, comparative analyses of mouse and human cardiac systems, and the rising contextualization of cardiac cells within their niches. As we explore the heart at this unprecedented resolution, integrating insights from both mouse and human studies will pave the way for novel diagnostic tools and therapeutic interventions, ultimately improving outcomes for patients with cardiovascular diseases.