Single-cell atlases: shared and tissue-specific cell types across human organs

Identifying similar populations across independent single cell studies without data integration

Supervised and unsupervised methods have emerged to address the complexity of single cell data analysis in the context of large pools of independent studies. Here, we present ClusterFoldSimilarity (CFS), a novel statistical method design to quantify the similarity between cell groups across any number of independent datasets, without the need for data correction or integration. By bypassing these processes, CFS avoids the introduction of artifacts and loss of information, offering a simple, efficient, and scalable solution. This method match groups of cells that exhibit conserved phenotypes across datasets, including different tissues and species, and in a multimodal scenario, including single-cell RNA-Seq, ATAC-Seq, single-cell proteomics, or, more broadly, data exhibiting differential abundance effects among groups of cells. Additionally, CFS performs feature selection, obtaining cross-dataset markers of the similar phenotypes observed, providing an inherent interpretability of relationships between cell populations. To showcase the effectiveness of our methodology, we generated single-nuclei RNA-Seq data from the motor cortex and spinal cord of adult mice. By using CFS, we identified three distinct sub-populations of astrocytes conserved on both tissues. CFS includes various visualization methods for the interpretation of the similarity scores and similar cell populations.

The impact of breast radiotherapy on the tumor genome and immune ecosystem

Radiotherapy is a pillar of breast cancer treatment; however, it remains unclear how radiotherapy modulates the tumor microenvironment. We investigated this question in a cohort of 20 patients with estrogen-receptor positive (ER+) breast tumors who received neoadjuvant radiotherapy. Tumor biopsies were collected before and 7 days postradiation. Single-cell DNA sequencing (scDNA-seq) and scRNA-seq were conducted on 8 and 11 patients, respectively, at these two time points. The scRNA data showed increased infiltration of naive-like CD4 T cells and an early, activated CD8 T cell population following radiotherapy. Radiotherapy also eliminated existing cytotoxic T cells and resulted in myeloid cell increases. In tumor cells, the scDNA-seq data showed a high genomic selection of subclones in half of the patients with high ER expression, while the remaining number had low genomic selection and an interferon response. Collectively, these data provide insight into the impact of radiotherapy in ER+ breast cancer patients.

Unveiling cell-type-specific mode of evolution in comparative single-cell expression data

While methodology for determining the mode of evolution in coding sequences has been well established, evaluation of adaptation events in emerging types of phenotype data needs further development. Here we propose an analysis framework (expression variance decomposition, EVaDe) for comparative single-cell expression data based on phenotypic evolution theory. After decomposing the gene expression variance into separate components, we use two strategies to identify genes exhibiting large between-taxon expression divergence and small within-cell-type expression noise in certain cell types, attributing this pattern to putative adaptive evolution. In a dataset of primate prefrontal cortex, we find that such human-specific key genes enrich with neurodevelopment-related functions, while most other genes exhibit neutral evolution patterns. Specific neuron types are found to harbor more of these key genes than other cell types, thus likely to have experienced more extensive adaptation. Reassuringly, at molecular sequence level, the key genes are significantly associated with the rapidly evolving conserved non-coding elements. An additional case analysis comparing the naked mole-rat (NMR) with the mouse suggests that innate-immunity-related genes and cell types have undergone putative expression adaptation in NMR. Overall, the EVaDe framework may effectively probe adaptive evolution mode in single-cell expression data.

Integrating microbial GWAS and single-cell transcriptomics reveals associations between host cell populations and the gut microbiome

Microbial genome-wide association studies (GWAS) have uncovered numerous host genetic variants associated with gut microbiota. However, links between host genetics, the gut microbiome and specific cellular contexts remain unclear. Here we use a computational framework, scBPS (single-cell Bacteria Polygenic Score), to integrate existing microbial GWAS and single-cell RNA-sequencing profiles of 24 human organs, including the liver, pancreas, lung and intestine, to identify host tissues and cell types relevant to gut microbes. Analysing 207 microbial taxa and 254 host cell types, scBPS-inferred cellular enrichments confirmed known biology such as dominant communications between gut microbes and the digestive tissue module and liver epithelial cell compartment. scBPS also identified a robust association between Collinsella and the central-veinal hepatocyte subpopulation. We experimentally validated the causal effects of Collinsella on cholesterol metabolism in mice through single-nuclei RNA sequencing on liver tissue to identify relevant cell subpopulations. Mechanistically, oral gavage of Collinsella modulated cholesterol pathway gene expression in central-veinal hepatocytes. We further validated our approach using independent microbial GWAS data, alongside single-cell and bulk transcriptomic analyses, demonstrating its robustness and reproducibility. Together, scBPS enables a systematic mapping of the host-microbe crosstalk by linking cell populations to their interacting gut microbes.

Advances and applications in single-cell and spatial genomics

The applications of single-cell and spatial technologies in recent times have revolutionized the present understanding of cellular states and the cellular heterogeneity inherent in complex biological systems. These advancements offer unprecedented resolution in the examination of the functional genomics of individual cells and their spatial context within tissues. In this review, we have comprehensively discussed the historical development and recent progress in the field of single-cell and spatial genomics. We have reviewed the breakthroughs in single-cell multi-omics technologies, spatial genomics methods, and the computational strategies employed toward the analyses of single-cell atlas data. Furthermore, we have highlighted the advances made in constructing cellular atlases and their clinical applications, particularly in the context of disease. Finally, we have discussed the emerging trends, challenges, and opportunities in this rapidly evolving field.

Spatial genomics: Mapping the landscape of fibrosis

Organ fibrosis causes major morbidity and mortality worldwide. Treatments for fibrosis are limited, with organ transplantation being the only cure. Here, we review how various state-of-the-art spatial genomics approaches are being deployed to interrogate fibrosis across multiple organs, providing exciting insights into fibrotic disease pathogenesis. These include the detailed topographical annotation of pathogenic cell populations and states, detection of transcriptomic perturbations in morphologically normal tissue, characterization of fibrotic and homeostatic niches and their cellular constituents, and in situ interrogation of ligand-receptor interactions within these microenvironments. Together, these powerful readouts enable detailed analysis of fibrosis evolution across time and space.

Harnessing RNA base editing for diverse applications in RNA biology and RNA therapeutics

Recent advancements in molecular engineering have established RNA-based technologies as powerful tools for both fundamental research and translational applications. Among the various RNA-based technologies developed, RNA base editing has recently emerged as a groundbreaking advancement. It primarily involves the conversion of adenosine (A) to inosine (I) and cytidine (C) to uridine (U), which are mediated by ADAR and APOBEC enzymes, respectively. RNA base editing has been applied in both biological research and therapeutic contexts. It enables site-directed editing within target transcripts, offering reversible, dose-dependent effects, in contrast to the permanent or heritable changes associated with DNA base editing. Additionally, RNA editing-based profiling of RNA-binding protein (RBP) binding sites facilitates transcriptome-wide mapping of RBP-RNA interactions in specific tissues and at the single-cell level. Furthermore, RNA editing-based sensors have been utilized to express effector proteins in response to specific RNA species. As RNA base editing technologies continue to evolve, we anticipate that they will significantly drive advancements in RNA therapeutics, synthetic biology, and biological research.

Recurrent Composite Markers of Cell Types and States

Biological function is mediated by the hierarchical organization of cell types and states within tissue ecosystems. Identifying interpretable composite marker sets that both define and distinguish hierarchical cell identities is essential for decoding biological complexity, yet remains a major challenge. Here, we present RECOMBINE, an algorithm that identifies recurrent composite marker sets to define hierarchical cell identities. Validation using both simulated and biological datasets demonstrates that RECOMBINE achieves higher accuracy in identifying discriminative markers compared to existing approaches, including differential gene expression analysis. When applied to single-cell data and validated with spatial transcriptomics data from the mouse visual cortex, RECOMBINE identified key cell type markers and generated a robust gene panel for targeted spatial profiling. It also uncovered markers of CD8+; T cell states, including GZMK+;HAVCR2-; effector memory cells associated with anti-PD-1 therapy response, and revealed a rare intestinal subpopulation with composite markers in mice. Finally, using data from the Tabula Sapiens project, RECOMBINE identified composite marker sets across a broad range of human tissues. Together, these results highlight RECOMBINE as a robust, data-driven framework for optimized marker selection, enabling the discovery and validation of hierarchical cell identities across diverse tissue contexts.

Endometriosis and adenomyosis unveiled through single-cell glasses

Single-cell technologies are expanding our understanding of endometriosis and adenomyosis, which are sister disorders of the uterine endometrium that contain similar complements of lesion cell types but in different locations-outside and inside the uterus, respectively. Both diseases cause significant morbidity and impaired quality of life among those affected, and current therapies mitigate most of the symptoms although with highly variable efficacy, duration of effect, and frequent intolerable side effects. Thus, there is a pressing need for transformative approaches and to develop individualized therapies for the variety of presentations of endometriosis and adenomyosis symptoms and the heterogeneity of lesion types, both histologically and architecturally. Single-cell technologies are transforming the understanding of human physiology and pathophysiology in the reproductive system and beyond. This manuscript reviews the clinical characteristics of endometriosis and adenomyosis and the recent studies focused on eutopic endometrium and ectopic lesions at single-cell resolution, the myriad of cell types and subtypes, cell-cell communications, signaling pathways, applications for novel drug discovery and therapeutic approaches, and challenges and opportunities that accompany this type of research. Key take-home messages from the studies reviewed herein include the following: We conclude the review with an eye to the future-what Alice might see beyond the single-cell looking glass that connects endometrium and endometrial disorders with the trillions of cells of other tissues and organs in health and disease throughout the human body and the opportunities for novel diagnostic modalities and drug discovery for endometriosis, adenomyosis, and related uterine and inflammatory conditions.

High-throughput single cell -omics using semi-permeable capsules

Biological systems are inherently complex and heterogeneous. Deciphering this complexity increasingly relies on high-throughput analytical methods and tools that efficiently probe the cellular phenotype and genotype. While recent advancements have enabled various single-cell -omics assays, their broader applications are inherently limited by the challenge of efficiently conducting multi-step biochemical assays while retaining various biological analytes. Extending on our previous work (1) here we present a versatile technology based on semi-permeable capsules (SPCs), tailored for a variety of high-throughput nucleic acid assays, including digital PCR, genome sequencing, single-cell RNA-sequencing (scRNA-Seq) and FACS-based isolation of individual transcriptomes based on nucleic acid marker of interest. Being biocompatible, the SPCs support single-cell cultivation and clonal expansion over long periods of time - a fundamental limitation of droplet microfluidics systems. Using SPCs we perform scRNA-Seq on white blood cells from patients with hematopoietic disorders and demonstrate that capsule-based sequencing approach (CapSeq) offers superior transcript capture, even for the most challenging cell types. By applying CapSeq on acute myeloid leukemia (AML) samples, we uncover notable changes in transcriptomes of mature granulocytes and monocytes associated with blast and progenitor cell phenotypes. Accurate representation of the entirety of the cellular heterogeneity of clinical samples, driving new insights into the malfunctioning of the innate immune system, and ability to clonally expand individual cells over long periods of time, positions SPC technology as customizable, highly sensitive and broadly applicable tool for easy-to-use, scalable single-cell -omics applications.

Mosquito Cell Atlas: A single-nucleus transcriptomic atlas of the adult Aedes aegypti mosquito

The female mosquito's remarkable ability to hunt humans and transmit pathogens relies on her unique biology. Here, we present the Mosquito Cell Atlas (MCA), a comprehensive single-nucleus RNA sequencing dataset of more than 367,000 nuclei from 19 dissected tissues of adult female and male Aedes aegypti, providing cellular-level resolution of mosquito biology. We identify novel cell types and expand our understanding of sensory neuron organization of chemoreceptors to all sensory tissues. Our analysis uncovers male-specific cells and sexually dimorphic gene expression in the antenna and brain. In female mosquitoes, we find that glial cells in the brain, rather than neurons, undergo the most extensive transcriptional changes following blood feeding. Our findings provide insights into the cellular basis of mosquito behavior and sexual dimorphism. The MCA aims to serve as a resource for the vector biology community, enabling systematic investigation of cell-type specific expression across all mosquito tissues.

A versatile and efficient method to isolate nuclei from low-input cryopreserved tissues for single-nuclei transcriptomics

Clinical samples are vital for understanding diseases, but their scarcity requires refined research methods. Emerging single-cell technologies offer detailed views of tissue heterogeneity but need sufficient fully characterized tissues. We developed an optimized single-nuclei RNA sequencing (snRNA-seq) protocol to extract nuclei from just 15 mg of cryopreserved human tissue. Applied to four cancer tissues (brain, bladder, lung, prostate), it profiled 1550-7468 nuclei per tissue, revealing heterogeneity comparable to public single-cell atlases. This method enhances the use and sharing of rare, cryopreserved biospecimens, supporting research where sample quantity is limited and full tissue characterization is needed.

scCobra allows contrastive cell embedding learning with domain adaptation for single cell data integration and harmonization

The rapid advancement of single-cell technologies has created an urgent need for effective methods to integrate and harmonize single-cell data. Technical and biological variations across studies complicate data integration, while conventional tools often struggle with reliance on gene expression distribution assumptions and over-correction. Here, we present scCobra, a deep generative neural network designed to overcome these challenges through contrastive learning with domain adaptation. scCobra effectively mitigates batch effects, minimizes over-correction, and ensures biologically meaningful data integration without assuming specific gene expression distributions. It enables online label transfer across datasets with batch effects, allowing continuous integration of new data without retraining. Additionally, scCobra supports batch effect simulation, advanced multi-omic integration, and scalable processing of large datasets. By integrating and harmonizing datasets from similar studies, scCobra expands the available data for investigating specific biological problems, improving cross-study comparability, and revealing insights that may be obscured in isolated datasets.

Robust collection and processing for label-free single voxel proteomics

With advanced mass spectrometry (MS)-based proteomics, genome-scale proteome coverage can be achieved from bulk tissues. However, such bulk measurement lacks spatial resolution and obscures tissue heterogeneity, precluding proteome mapping of tissue microenvironment. Here we report an integrated wet collection of single microscale tissue voxels and Surfactant-assisted One-Pot voxel processing method termed wcSOP for robust label-free single voxel proteomics. wcSOP capitalizes on buffer droplet-assisted wet collection of single voxels dissected by LCM to the tube cap and SOP voxel processing in the same collection cap. This method enables reproducible, label-free quantification of approximately 900 and 4600 proteins for single voxels at 20 µm × 20 µm × 10 µm (~1 cell region) and 200 µm × 200 µm × 10 µm (~100 cell region) from fresh frozen human spleen tissue, respectively. It can reveal spatially resolved protein signatures and region-specific signaling pathways. Furthermore, wcSOP-MS is demonstrated to be broadly applicable for OCT-embedded and FFPE human archived tissues as well as for small-scale 2D proteome mapping of tissues at high spatial resolutions. wcSOP-MS may pave the way for routine robust single voxel proteomics and spatial proteomics.

inDrops-2: a flexible, versatile and cost-efficient droplet microfluidic approach for high-throughput scRNA-seq of fresh and preserved clinical samples

The expansion of single-cell analytical techniques has empowered the exploration of diverse biological questions at the individual cells. Droplet-based single-cell RNA sequencing (scRNA-seq) methods have been particularly widely used due to their high-throughput capabilities and small reaction volumes. While commercial systems have contributed to the widespread adoption of droplet-based scRNA-seq, their relatively high cost limits the ability to profile large numbers of cells and samples. Moreover, as the scale of single-cell sequencing continues to expand, accommodating diverse workflows and cost-effective multi-biospecimen profiling becomes more critical. Herein, we present inDrops-2, an open-source scRNA-seq technology designed to profile live or preserved cells with a sensitivity matching that of state-of-the-art commercial systems but at a 6-fold lower cost. We demonstrate the flexibility of inDrops-2, by implementing two prominent scRNA-seq protocols, based on exponential and linear amplification of barcoded-complementary DNA, and provide useful insights into the advantages and disadvantages inherent to each approach. We applied inDrops-2 to simultaneously profile multiple human lung carcinoma samples that had been subjected to cell preservation, long-term storage and multiplexing to obtain a multiregional cellular profile of the tumor microenvironment. The scalability, sensitivity and cost efficiency make inDrops-2 stand out among other droplet-based scRNA-seq methods, ideal for large-scale studies on rare cell molecular signatures.

Control of Cellular Differentiation Trajectories for Cancer Reversion

Cellular differentiation is controlled by intricate layers of gene regulation, involving the modulation of gene expression by various transcriptional regulators. Due to the complexity of gene regulation, identifying master regulators across the differentiation trajectory has been a longstanding challenge. To tackle this problem, a computational framework, single-cell Boolean network inference and control (BENEIN), is presented. Applying BENEIN to human large intestinal single-cell transcriptome data, MYB, HDAC2, and FOXA2 are identified as the master regulators whose inhibition induces enterocyte differentiation. It is found that simultaneous knockdown of these master regulators can revert colorectal cancer cells into normal-like enterocytes by synergistically inducing differentiation and suppressing malignancy, which is validated by in vitro and in vivo experiments.

Nuclei Isolation from Ocular Tissues of the Embryonic Chicken for Single-Nucleus Profiling

The generation of quality data from a single-nucleus profiling experiment requires nuclei to be isolated from tissues in a gentle and efficient manner. Nuclei isolation must be carefully optimized across tissue types to preserve nuclear architecture, prevent nucleic acid degradation, and remove unwanted contaminants. Here, we present an optimized workflow for generating a single-nucleus suspension from ocular tissues of the embryonic chicken that is compatible with various downstream workflows. The described protocol enables the rapid isolation of a high yield of aggregate-free nuclei from the embryonic chicken eye without compromising nucleic acid integrity, and the nuclei suspension is compatible with single-nucleus RNA and ATAC sequencing. We detail several stopping points, either via cryopreservation or fixation, to enhance workflow adaptability. Further, we provide a guide through multiple QC points and demonstrate proof-of-principle using two commercially available kits. Finally, we demonstrate that existing in silico genotyping methods can be adopted to computationally derive biological replicates from a single pool of chicken nuclei, greatly reducing the cost of biological replication and allowing researchers to consider sex as a variable during analysis. Together, this tutorial represents a cost-effective, simple, and effective approach to single-nucleus profiling of embryonic chicken eye tissues and is likely to be easily modified to be compatible with similar tissue types.

A CH4-Driven Ion Cloud-Stretched Approach Enables ICP-qMS for Multiplex Single-Cell Analysis

In the last 40 years, inductively coupled plasma quadrupole (q) mass spectrometry (ICP-qMS) has been recognized as one of the best tools for the quantification of multiple elements/isotopes and even the biomolecules they labeled in a homogeneous solution sample. However, it meets a tough challenge when acquiring multi-m/z signals from an intact single-cell dispersed in a cell suspension, since the single-cell ion cloud generated in ICP presents an intermittently transient event with a duration time of hundreds of microseconds while the dwell time plus settling time of the q is at the similar time scale when peak-hopping between different m/z. Herein, we report CH4 is able to stretch the single-cell ion cloud duration time to more than 7,000 μs in collision-reaction-cell (CRC), allowing multi-m/z signals acquisition by ICP-qMS. Quantification of single-cell's multiple phenotype protein markers can thus be achieved on ICP-(CH4-CRC)-qMS, not only revealing the heterogeneity between the single cells but also enabling an unambiguous cell-classification of their subtypes. CH4-driven ion cloud-stretched approach breaks through the long-standing bottleneck limited single-cell multiplex analysis on ICP-qMS, paving a path for more important applications of ICP-qMS in the fields related to single-cell analysis.

Single-cell transcriptomics reveal distinctive patterns of fibroblast activation in heart failure with preserved ejection fraction

Inflammation, fibrosis and metabolic stress critically promote heart failure with preserved ejection fraction (HFpEF). Exposure to high-fat diet and nitric oxide synthase inhibitor N[w]-nitro-l-arginine methyl ester (L-NAME) recapitulate features of HFpEF in mice. To identify disease-specific traits during adverse remodeling, we profiled interstitial cells in early murine HFpEF using single-cell RNAseq (scRNAseq). Diastolic dysfunction and perivascular fibrosis were accompanied by an activation of cardiac fibroblast and macrophage subsets. Integration of fibroblasts from HFpEF with two murine models for heart failure with reduced ejection fraction (HFrEF) identified a catalog of conserved fibroblast phenotypes across mouse models. Moreover, HFpEF-specific characteristics included induced metabolic, hypoxic and inflammatory transcription factors and pathways, including enhanced expression of Angiopoietin-like 4 (Angptl4) next to basement membrane compounds, such as collagen IV (Col4a1). Fibroblast activation was further dissected into transcriptional and compositional shifts and thereby highly responsive cell states for each HF model were identified. In contrast to HFrEF, where myofibroblast and matrifibrocyte activation were crucial features, we found that these cell states played a subsidiary role in early HFpEF. These disease-specific fibroblast signatures were corroborated in human myocardial bulk transcriptomes. Furthermore, we identified a potential cross-talk between macrophages and fibroblasts via SPP1 and TNFɑ with estimated fibroblast target genes including Col4a1 and Angptl4. Treatment with recombinant ANGPTL4 ameliorated the murine HFpEF phenotype and diastolic dysfunction by reducing collagen IV deposition from fibroblasts in vivo and in vitro. In line, ANGPTL4, was elevated in plasma samples of HFpEF patients and particularly high levels associated with a preserved global-longitudinal strain. Taken together, our study provides a comprehensive characterization of molecular fibroblast activation patterns in murine HFpEF, as well as the identification of Angiopoietin-like 4 as central mechanistic regulator with protective effects.

π-HuB: the proteomic navigator of the human body

The human body contains trillions of cells, classified into specific cell types, with diverse morphologies and functions. In addition, cells of the same type can assume different states within an individual's body during their lifetime. Understanding the complexities of the proteome in the context of a human organism and its many potential states is a necessary requirement to understanding human biology, but these complexities can neither be predicted from the genome, nor have they been systematically measurable with available technologies. Recent advances in proteomic technology and computational sciences now provide opportunities to investigate the intricate biology of the human body at unprecedented resolution and scale. Here we introduce a big-science endeavour called π-HuB (proteomic navigator of the human body). The aim of the π-HuB project is to (1) generate and harness multimodality proteomic datasets to enhance our understanding of human biology; (2) facilitate disease risk assessment and diagnosis; (3) uncover new drug targets; (4) optimize appropriate therapeutic strategies; and (5) enable intelligent healthcare, thereby ushering in a new era of proteomics-driven phronesis medicine. This ambitious mission will be implemented by an international collaborative force of multidisciplinary research teams worldwide across academic, industrial and government sectors.

Single cell variant to enhancer to gene map for coronary artery disease

Although genome wide association studies (GWAS) in large populations have identified hundreds of variants associated with common diseases such as coronary artery disease (CAD), most disease-associated variants lie within non-coding regions of the genome, rendering it difficult to determine the downstream causal gene and cell type. Here, we performed paired single nucleus gene expression and chromatin accessibility profiling from 44 human coronary arteries. To link disease variants to molecular traits, we developed a meta-map of 88 samples and discovered 11,182 single-cell chromatin accessibility quantitative trait loci (caQTLs). Heritability enrichment analysis and disease variant mapping demonstrated that smooth muscle cells (SMCs) harbor the greatest genetic risk for CAD. To capture the continuum of SMC cell states in disease, we used dynamic single cell caQTL modeling for the first time in tissue to uncover QTLs whose effects are modified by cell state and expand our insight into genetic regulation of heterogenous cell populations. Notably, we identified a variant in the COL4A1/COL4A2 CAD GWAS locus which becomes a caQTL as SMCs de-differentiate by changing a transcription factor binding site for EGR1/2. To unbiasedly prioritize functional candidate genes, we built a genome-wide single cell variant to enhancer to gene (scV2E2G) map for human CAD to link disease variants to causal genes in cell types. Using this approach, we found several hundred genes predicted to be linked to disease variants in different cell types. Next, we performed genome-wide Hi-C in 16 human coronary arteries to build tissue specific maps of chromatin conformation and link disease variants to integrated chromatin hubs and distal target genes. Using this approach, we show that rs4887091 within the ADAMTS7 CAD GWAS locus modulates function of a super chromatin interactome through a change in a CTCF binding site. Finally, we used CRISPR interference to validate a distal gene, AMOTL2, liked to a CAD GWAS locus. Collectively we provide a disease-agnostic framework to translate human genetic findings to identify pathologic cell states and genes driving disease, producing a comprehensive scV2E2G map with genetic and tissue level convergence for future mechanistic and therapeutic studies.

Molecular mechanisms of lipid droplets-mitochondria coupling in obesity and metabolic syndrome: insights and pharmacological implications

Abnormal lipid accumulation is a fundamental contributor to obesity and metabolic disorders. Lipid droplets (LDs) and mitochondria (MT) serve as organelle chaperones in lipid metabolism and energy balance. LDs play a crucial role in lipid storage and mobilization, working in conjunction with MT to regulate lipid metabolism within the liver, brown adipose tissue, and skeletal muscle, thereby maintaining metabolic homeostasis. The novelty of our review is the comprehensive description of LD and MT interaction mechanisms. We also focus on the current drugs that target this metabolism, which provide novel approaches for obesity and related metabolism disorder treatment.

Data-driven batch detection enhances single-cell omics data analysis

In single-cell omics studies, data are typically collected across multiple batches, resulting in batch effects: technical confounders that introduce noise and distort data distribution. Correcting these effects is challenging due to their unknown sources, nonlinear distortions, and the difficulty of accurately assigning data to batches that are optimal for integration methods.

Cross-tissue human fibroblast atlas reveals myofibroblast subtypes with distinct roles in immune modulation

Fibroblasts, known for their functional diversity, play crucial roles in inflammation and cancer. In this study, we conduct comprehensive single-cell RNA sequencing analyses on fibroblast cells from 517 human samples, spanning 11 tissue types and diverse pathological states. We identify distinct fibroblast subpopulations with universal and tissue-specific characteristics. Pathological conditions lead to significant shifts in fibroblast compositions, including the expansion of immune-modulating fibroblasts during inflammation and tissue-remodeling myofibroblasts in cancer. Within the myofibroblast category, we identify four transcriptionally distinct subpopulations originating from different developmental origins, with LRRC15+ myofibroblasts displaying terminally differentiated features. Both LRRC15+ and MMP1+ myofibroblasts demonstrate pro-tumor potential that contribute to the immune-excluded and immune-suppressive tumor microenvironments (TMEs), whereas PI16+ fibroblasts show potential anti-tumor functions in adjacent non-cancerous regions. Fibroblast-subtype compositions define patient subtypes with distinct clinical outcomes. This study advances our understanding of fibroblast biology and suggests potential therapeutic strategies for targeting specific fibroblast subsets in cancer treatment.

Universal and tissue-specific fibroblasts in chronic inflammation and cancer

In this issue of Cancer Cell, Gao et al. map fibroblast diversity across tumors and chronic inflammatory tissues. The authors uncover universal fibroblast subtypes such as LRRC15+ and MMP1+ myofibroblasts along with specialized tissue-specific subtypes. They reveal cellular roles of fibroblasts in immunosuppression through stromal niches and cell-cell interactions.

Advancements in Single-Cell RNA Sequencing and Spatial Transcriptomics for Central Nervous System Disease

The incidence of central nervous system (CNS) disease has persistently increased over the last several years. There is an urgent need for effective methods to improve the cure rates of CNS disease. However, the precise molecular basis underlying the development and progression of major CNS diseases remains elusive. A complete molecular map will contribute to research on CNS disease treatment strategies. Emerging technologies such as single-cell RNA sequencing (scRNA-seq) and Spatial Transcriptomics (ST) are potent tools for exploring the molecular complexity, cell heterogeneity, and functional specificity of CNS disease. scRNA-seq and ST can provide insights into the disease at cellular and spatial transcription levels. This review presents a survey of scRNA-seq and ST studies on CNS diseases, such as chronic neurodegenerative diseases, acute CNS injuries, and others. These studies offer novel perspectives in treating and diagnosing CNS diseases by discovering new cell types or subtypes associated with the disease, proposing new pathophysiological mechanisms, uncovering novel therapeutic targets, and identifying putative biomarkers.

Overcoming barriers to single-cell RNA sequencing adoption in low- and middle-income countries

The advent of single-cell resolution sequencing and spatial transcriptomics has enabled the delivery of cellular and molecular atlases of tissues and organs, providing new insights into tissue health and disease. However, if the full potential of these technologies is to be equitably realised, ancestrally inclusivity is paramount. Such a goal requires greater inclusion of both researchers and donors in low- and middle-income countries (LMICs). In this perspective, we describe the current landscape of ancestral inclusivity in genomic and single-cell transcriptomic studies. We discuss the collaborative efforts needed to scale the barriers to establishing, expanding, and adopting single-cell sequencing research in LMICs and to enable globally impactful outcomes of these technologies.

New frameworks for hematopoiesis derived from single-cell genomics

ABSTRACT

Recent advancements in single-cell genomics have enriched our understanding of hematopoiesis, providing intricate details about hematopoietic stem cell biology, differentiation, and lineage commitment. Technological advancements have highlighted extensive heterogeneity of cell populations and continuity of differentiation routes. Nevertheless, intermediate "attractor" states signify structure in stem and progenitor populations that link state transition dynamics to fate potential. We discuss how innovative model systems quantify lineage bias and how stress accelerates differentiation, thereby reducing fate plasticity compared with native hematopoiesis. We conclude by offering our perspective on the current model of hematopoiesis and discuss how a more precise understanding can translate to strategies that extend healthy hematopoiesis and prevent disease.

Spatial immunophenotyping of FFPE tissues by imaging mass cytometry

The immune compartment of a tissue is dynamic, changing to respond to infections, tumors, or therapeutic interventions. Within tissues, local microenvironments provide interaction partners and cytokines that can gear immune cells into distinct functional states. Thus, it is not just the immune composition of a tissue, but also the relative localization of immune cells that determines the outcome of a response. Conventional techniques like immunohistochemistry (IHC) have been used to describe infiltration of immune cells and their relative position within tissues. However, these technologies are limited on the number of targets that can be simultaneously imaged. Here, we describe a simple protocol using imaging mass cytometry (IMC) for immunophenotyping formalin-fixed, paraffin-embedded (FFPE) tissues. IMC has a 1-μm resolution and allows simultaneous detection of up to 40 targets, overcoming limitations of traditional methods. In this protocol, we detail the staining procedure, offer an example of a murine FFPE antibody panel for immunophenotyping, and additionally provide suggestions for initial image analysis. The herein presented workflow facilitates the characterization of immune niches and can be used to assess their alterations throughout immune responses or therapeutic interventions. With minimal alterations, this approach can be used on clinically relevant samples or animal models to investigate specific immune responses and better understand disease progression or treatment dynamics.

Building a learnable universal coordinate system for single-cell atlas with a joint-VAE model

A universal coordinate system that can ensemble the huge number of cells and capture their heterogeneities is of vital importance for constructing large-scale cell atlases as references for molecular and cellular studies. Studies have shown that cells exhibit multifaceted heterogeneities in their transcriptomic features at multiple resolutions. This nature of complexity makes it hard to design a fixed coordinate system through a combination of known features. It is desirable to build a learnable universal coordinate model that can capture major heterogeneities and serve as a controlled generative model for data augmentation. We developed UniCoord, a specially-tuned joint-VAE model to represent single-cell transcriptomic data in a lower-dimensional latent space with high interpretability. Each latent dimension can represent either discrete or continuous feature, and either supervised by prior knowledge or unsupervised. The latent dimensions can be easily reconfigured to generate pseudo transcriptomic profiles with desired properties. UniCoord can also be used as a pre-trained model to analyze new data with unseen cell types and thus can serve as a feasible framework for cell annotation and comparison. UniCoord provides a prototype for a learnable universal coordinate framework to enable better analysis and generation of cells with highly orchestrated functions and heterogeneities.

Scalable spatial single-cell transcriptomics and translatomics in 3D thick tissue blocks

Characterizing the transcriptional and translational gene expression patterns at the single-cell level within their three-dimensional (3D) tissue context is essential for revealing how genes shape tissue structure and function in health and disease. However, most existing spatial profiling techniques are limited to 5-20 μm thin tissue sections. Here, we developed Deep-STARmap and Deep-RIBOmap, which enable 3D in situ quantification of thousands of gene transcripts and their corresponding translation activities, respectively, within 200-μm thick tissue blocks. This is achieved through scalable probe synthesis, hydrogel embedding with efficient probe anchoring, and robust cDNA crosslinking. We first utilized Deep-STARmap in combination with multicolor fluorescent protein imaging for simultaneous molecular cell typing and 3D neuron morphology tracing in the mouse brain. We also demonstrate that 3D spatial profiling facilitates comprehensive and quantitative analysis of tumor-immune interactions in human skin cancer.

aKNNO: single-cell and spatial transcriptomics clustering with an optimized adaptive k-nearest neighbor graph

Typical clustering methods for single-cell and spatial transcriptomics struggle to identify rare cell types, while approaches tailored to detect rare cell types gain this ability at the cost of poorer performance for grouping abundant ones. Here, we develop aKNNO to simultaneously identify abundant and rare cell types based on an adaptive k-nearest neighbor graph with optimization. Benchmarking on 38 simulated and 20 single-cell and spatial transcriptomics datasets demonstrates that aKNNO identifies both abundant and rare cell types more accurately than general and specialized methods. Using only gene expression aKNNO maps abundant and rare cells more precisely compared to integrative approaches.

A systematic overview of single-cell transcriptomics databases, their use cases, and limitations

Rapid advancements in high-throughput single-cell RNA-seq (scRNA-seq) technologies and experimental protocols have led to the generation of vast amounts of transcriptomic data that populates several online databases and repositories. Here, we systematically examined large-scale scRNA-seq databases, categorizing them based on their scope and purpose such as general, tissue-specific databases, disease-specific databases, cancer-focused databases, and cell type-focused databases. Next, we discuss the technical and methodological challenges associated with curating large-scale scRNA-seq databases, along with current computational solutions. We argue that understanding scRNA-seq databases, including their limitations and assumptions, is crucial for effectively utilizing this data to make robust discoveries and identify novel biological insights. Such platforms can help bridge the gap between computational and wet lab scientists through user-friendly web-based interfaces needed for democratizing access to single-cell data. These platforms would facilitate interdisciplinary research, enabling researchers from various disciplines to collaborate effectively. This review underscores the importance of leveraging computational approaches to unravel the complexities of single-cell data and offers a promising direction for future research in the field.

Deciphering metabolic crosstalk in context: lessons from inflammatory diseases

Metabolism plays a crucial role in regulating the function of immune cells in both health and disease, with altered metabolism contributing to the pathogenesis of cancer and many inflammatory diseases. The local microenvironment has a profound impact on the metabolism of immune cells. Therefore, immunological and metabolic heterogeneity as well as the spatial organization of cells in tissues should be taken into account when studying immunometabolism. Here, we highlight challenges of investigating metabolic communication. Additionally, we review the capabilities and limitations of current technologies for studying metabolism in inflamed microenvironments, including single-cell omics techniques, flow cytometry-based methods (Met-Flow, single-cell energetic metabolism by profiling translation inhibition (SCENITH)), cytometry by time of flight (CyTOF), cellular indexing of transcriptomes and epitopes by sequencing (CITE-Seq), and mass spectrometry imaging. Considering the importance of metabolism in regulating immune cells in diseased states, we also discuss the applications of metabolomics in clinical research, as well as some hurdles to overcome to implement these techniques in standard clinical practice. Finally, we provide a flowchart to assist scientists in designing effective strategies to unravel immunometabolism in disease-relevant contexts.

Aging atlas reveals cell-type-specific effects of pro-longevity strategies

Organismal aging involves functional declines in both somatic and reproductive tissues. Multiple strategies have been discovered to extend lifespan across species. However, how age-related molecular changes differ among various tissues and how those lifespan-extending strategies slow tissue aging in distinct manners remain unclear. Here we generated the transcriptomic Cell Atlas of Worm Aging (CAWA, http://mengwanglab.org/atlas ) of wild-type and long-lived strains. We discovered cell-specific, age-related molecular and functional signatures across all somatic and germ cell types. We developed transcriptomic aging clocks for different tissues and quantitatively determined how three different pro-longevity strategies slow tissue aging distinctively. Furthermore, through genome-wide profiling of alternative polyadenylation (APA) events in different tissues, we discovered cell-type-specific APA changes during aging and revealed how these changes are differentially affected by the pro-longevity strategies. Together, this study offers fundamental molecular insights into both somatic and reproductive aging and provides a valuable resource for in-depth understanding of the diversity of pro-longevity mechanisms.

Rat Models of Hormone Receptor-Positive Breast Cancer

Hormone receptor-positive (HR+) breast cancer (BC) is the most common type of breast cancer among women worldwide, accounting for 70-80% of all invasive cases. Patients with HR+ BC are commonly treated with endocrine therapy, but intrinsic or acquired resistance is a frequent problem, making HR+ BC a focal point of intense research. Despite this, the malignancy still lacks adequate in vitro and in vivo models for the study of its initiation and progression as well as response and resistance to endocrine therapy. No mouse models that fully mimic the human disease are available, however rat mammary tumor models pose a promising alternative to overcome this limitation. Compared to mice, rats are more similar to humans in terms of mammary gland architecture, ductal origin of neoplastic lesions and hormone dependency status. Moreover, rats can develop spontaneous or induced mammary tumors that resemble human HR+ BC. To date, six different types of rat models of HR+ BC have been established. These include the spontaneous, carcinogen-induced, transplantation, hormone-induced, radiation-induced and genetically engineered rat mammary tumor models. Each model has distinct advantages, disadvantages and utility for studying HR+ BC. This review provides a comprehensive overview of all published models to date.

Machine learning integrative approaches to advance computational immunology

The study of immunology, traditionally reliant on proteomics to evaluate individual immune cells, has been revolutionized by single-cell RNA sequencing. Computational immunologists play a crucial role in analysing these datasets, moving beyond traditional protein marker identification to encompass a more detailed view of cellular phenotypes and their functional roles. Recent technological advancements allow the simultaneous measurements of multiple cellular components-transcriptome, proteome, chromatin, epigenetic modifications and metabolites-within single cells, including in spatial contexts within tissues. This has led to the generation of complex multiscale datasets that can include multimodal measurements from the same cells or a mix of paired and unpaired modalities. Modern machine learning (ML) techniques allow for the integration of multiple "omics" data without the need for extensive independent modelling of each modality. This review focuses on recent advancements in ML integrative approaches applied to immunological studies. We highlight the importance of these methods in creating a unified representation of multiscale data collections, particularly for single-cell and spatial profiling technologies. Finally, we discuss the challenges of these holistic approaches and how they will be instrumental in the development of a common coordinate framework for multiscale studies, thereby accelerating research and enabling discoveries in the computational immunology field.

A Molecular Voyage: Multiomics Insights into Circulating Tumor Cells

Circulating tumor cells (CTCs) play a pivotal role in metastasis, the leading cause of cancer-associated death. Recent improvements of CTC isolation tools, coupled with a steady development of multiomics technologies at single-cell resolution, have enabled an extensive exploration of CTC biology, unlocking insights into their molecular profiles. A detailed molecular portrait requires CTC interrogation across various levels encompassing genomic, epigenetic, transcriptomic, proteomic and metabolic features. Here, we review how state-of-the-art multiomics applied to CTCs are shedding light on how cancer spreads. Further, we highlight the potential implications of CTC profiling for clinical applications aimed at enhancing cancer diagnosis and treatment.

SIGNIFICANCE

Exploring the complexity of cancer progression through cutting-edge multiomics studies holds the promise of uncovering novel aspects of cancer biology and identifying therapeutic vulnerabilities to suppress metastasis.

Deep generative models in single-cell omics

Deep Generative Models (DGMs) are becoming instrumental for inferring probability distributions inherent to complex processes, such as most questions in biomedical research. For many years, there was a lack of mathematical methods that would allow this inference in the scarce data scenario of biomedical research. The advent of single-cell omics has finally made square the so-called "skinny matrix", allowing to apply mathematical methods already extensively used in other areas. Moreover, it is now possible to integrate data at different molecular levels in thousands or even millions of samples, thanks to the number of single-cell atlases being collaboratively generated. Additionally, DGMs have proven useful in other frequent tasks in single-cell analysis pipelines, from dimensionality reduction, cell type annotation to RNA velocity inference. In spite of its promise, DGMs need to be used with caution in biomedical research, paying special attention to its use to answer the right questions and the definition of appropriate error metrics and validation check points that confirm not only its correct use but also its relevance. All in all, DGMs provide an exciting tool that opens a bright future for the integrative analysis of single-cell -omics to understand health and disease.

Characterizing efficient feature selection for single-cell expression analysis

Unsupervised feature selection is a critical step for efficient and accurate analysis of single-cell RNA-seq data. Previous benchmarks used two different criteria to compare feature selection methods: (i) proportion of ground-truth marker genes included in the selected features and (ii) accuracy of cell clustering using ground-truth cell types. Here, we systematically compare the performance of 11 feature selection methods for both criteria. We first demonstrate the discordance between these criteria and suggest using the latter. We then compare the distribution of selected genes in their means between feature selection methods. We show that lowly expressed genes exhibit seriously high coefficients of variation and are mostly excluded by high-performance methods. In particular, high-deviation- and high-expression-based methods outperform the widely used in Seurat package in clustering cells and data visualization. We further show they also enable a clear separation of the same cell type from different tissues as well as accurate estimation of cell trajectories.

Stochasticity in genetics and gene regulation

Development from fertilized egg to functioning multi-cellular organism requires precision. There is no precision, and often no survival, without plasticity. Plasticity is conferred partly by stochastic variation, present inherently in all biological systems. Gene expression levels fluctuate ubiquitously through transcription, alternative splicing, translation and turnover. Small differences in gene expression are exploited to trigger early differentiation, conferring distinct function on selected individual cells and setting in motion regulatory interactions. Non-selected cells then acquire new functions along the spatio-temporal developmental trajectory. The differentiation process has many stochastic components. Meiotic segregation, mitochondrial partitioning, X-inactivation and the dynamic DNA binding of transcription factor assemblies-all exhibit randomness. Non-random X-inactivation generally signals deleterious X-linked mutations. Correct neural wiring, such as retina to brain, arises through repeated confirmatory activity of connections made randomly. In immune system development, both B-cell antibody generation and the emergence of balanced T-cell categories begin through stochastic trial and error followed by functional selection. Aberrant selection processes lead to immune dysfunction. DNA sequence variants also arise through stochastic events: some involving environmental fluctuation (radiation or presence of pollutants), or genetic repair system malfunction. The phenotypic outcome of mutations is also fluid. Mutations may be advantageous in some circumstances, deleterious in others. This article is part of a discussion meeting issue 'Causes and consequences of stochastic processes in development and disease'.

Single Cell Atlas: a single-cell multi-omics human cell encyclopedia

Single-cell sequencing datasets are key in biology and medicine for unraveling insights into heterogeneous cell populations with unprecedented resolution. Here, we construct a single-cell multi-omics map of human tissues through in-depth characterizations of datasets from five single-cell omics, spatial transcriptomics, and two bulk omics across 125 healthy adult and fetal tissues. We construct its complement web-based platform, the Single Cell Atlas (SCA, www.singlecellatlas.org ), to enable vast interactive data exploration of deep multi-omics signatures across human fetal and adult tissues. The atlas resources and database queries aspire to serve as a one-stop, comprehensive, and time-effective resource for various omics studies.

Prediction of protein-RNA interactions from single-cell transcriptomic data

Proteins are crucial in regulating every aspect of RNA life, yet understanding their interactions with coding and noncoding RNAs remains limited. Experimental studies are typically restricted to a small number of cell lines and a limited set of RNA-binding proteins (RBPs). Although computational methods based on physico-chemical principles can predict protein-RNA interactions accurately, they often lack the ability to consider cell-type-specific gene expression and the broader context of gene regulatory networks (GRNs). Here, we assess the performance of several GRN inference algorithms in predicting protein-RNA interactions from single-cell transcriptomic data, and propose a pipeline, called scRAPID (single-cell transcriptomic-based RnA Protein Interaction Detection), that integrates these methods with the catRAPID algorithm, which can identify direct physical interactions between RBPs and RNA molecules. Our approach demonstrates that RBP-RNA interactions can be predicted from single-cell transcriptomic data, with performances comparable or superior to those achieved for the well-established task of inferring transcription factor-target interactions. The incorporation of catRAPID significantly enhances the accuracy of identifying interactions, particularly with long noncoding RNAs, and enables the identification of hub RBPs and RNAs. Additionally, we show that interactions between RBPs can be detected based on their inferred RNA targets. The software is freely available at https://github.com/tartaglialabIIT/scRAPID.

scButterfly: a versatile single-cell cross-modality translation method via dual-aligned variational autoencoders

Recent advancements for simultaneously profiling multi-omics modalities within individual cells have enabled the interrogation of cellular heterogeneity and molecular hierarchy. However, technical limitations lead to highly noisy multi-modal data and substantial costs. Although computational methods have been proposed to translate single-cell data across modalities, broad applications of the methods still remain impeded by formidable challenges. Here, we propose scButterfly, a versatile single-cell cross-modality translation method based on dual-aligned variational autoencoders and data augmentation schemes. With comprehensive experiments on multiple datasets, we provide compelling evidence of scButterfly's superiority over baseline methods in preserving cellular heterogeneity while translating datasets of various contexts and in revealing cell type-specific biological insights. Besides, we demonstrate the extensive applications of scButterfly for integrative multi-omics analysis of single-modality data, data enhancement of poor-quality single-cell multi-omics, and automatic cell type annotation of scATAC-seq data. Moreover, scButterfly can be generalized to unpaired data training, perturbation-response analysis, and consecutive translation.

Technology-enabled great leap in deciphering plant genomes

Plant genomes provide essential and vital basic resources for studying many aspects of plant biology and applications (for example, breeding). From 2000 to 2020, 1,144 genomes of 782 plant species were sequenced. In the past three years (2021-2023), 2,373 genomes of 1,031 plant species, including 793 newly sequenced species, have been assembled, representing a great leap. The 2,373 newly assembled genomes, of which 63 are telomere-to-telomere assemblies and 921 have been generated in pan-genome projects, cover the major phylogenetic clades. Substantial advances in read length, throughput, accuracy and cost-effectiveness have notably simplified the achievement of high-quality assemblies. Moreover, the development of multiple software tools using different algorithms offers the opportunity to generate more complete and complex assemblies. A database named N3: plants, genomes, technologies has been developed to accommodate the metadata associated with the 3,517 genomes that have been sequenced from 1,575 plant species since 2000. We also provide an outlook for emerging opportunities in plant genome sequencing.

Cellular heterogeneity and dynamics of the human uterus in healthy premenopausal women

The human uterus is a complex and dynamic organ whose lining grows, remodels, and regenerates in every menstrual cycle or upon tissue damage. Here we applied single-cell RNA sequencing to profile more the 50,000 uterine cells from both the endometrium and myometrium of 5 healthy premenopausal individuals, and jointly analyzed the data with a previously published dataset from 15 subjects. The resulting normal uterus cell atlas contains more than 167K cells representing the lymphatic endothelium, blood endothelium, stromal, ciliated epithelium, unciliated epithelium, and immune cell populations. Focused analyses within each major cell type and comparisons with subtype labels from prior studies allowed us to document supporting evidence, resolve naming conflicts, and to propose a consensus annotation system of 39 subtypes. We release their gene expression centroids, differentially expressed genes, and mRNA patterns of literature-based markers as a shared community resource. We find many subtypes show dynamic changes over different phases of the cycle and identify multiple potential progenitor cells: compartment-wide progenitors for each major cell type, transitional cells that are upstream of other subtypes, and potential cross-lineage multipotent stromal progenitors that may be capable of replenishing the epithelial, stromal, and endothelial compartments. When compared to the healthy premenopausal samples, a postpartum and a postmenopausal uterus sample revealed substantially altered tissue composition, involving the rise or fall of stromal, endothelial, and immune cells. The cell taxonomy and molecular markers we report here are expected to inform studies of both basic biology of uterine function and its disorders.

SIGNIFICANCE

We present single-cell RNA sequencing data from seven individuals (five healthy pre-menopausal women, one post-menopausal woman, and one postpartum) and perform an integrated analysis of this data alongside 15 previously published scRNA-seq datasets. We identified 39 distinct cell subtypes across four major cell types in the uterus. By using RNA velocity analysis and centroid-centroid comparisons we identify multiple computationally predicted progenitor populations for each of the major cell compartments, as well as potential cross-compartment, multi-potent progenitors. While the function and interactions of these cell populations remain to be validated through future experiments, the markers and their "dual characteristics" that we describe will serve as a rich resource to the scientific community. Importantly, we address a significant challenge in the field: reconciling multiple uterine cell taxonomies being proposed. To achieve this, we focused on integrating historical and contemporary knowledge across multiple studies. By providing detailed evidence used for cell classification we lay the groundwork for establishing a stable, consensus cell atlas of the human uterus.

Principled and interpretable alignability testing and integration of single-cell data

Single-cell data integration can provide a comprehensive molecular view of cells, and many algorithms have been developed to remove unwanted technical or biological variations and integrate heterogeneous single-cell datasets. Despite their wide usage, existing methods suffer from several fundamental limitations. In particular, we lack a rigorous statistical test for whether two high-dimensional single-cell datasets are alignable (and therefore should even be aligned). Moreover, popular methods can substantially distort the data during alignment, making the aligned data and downstream analysis difficult to interpret. To overcome these limitations, we present a spectral manifold alignment and inference (SMAI) framework, which enables principled and interpretable alignability testing and structure-preserving integration of single-cell data with the same type of features. SMAI provides a statistical test to robustly assess the alignability between datasets to avoid misleading inference and is justified by high-dimensional statistical theory. On a diverse range of real and simulated benchmark datasets, it outperforms commonly used alignment methods. Moreover, we show that SMAI improves various downstream analyses such as identification of differentially expressed genes and imputation of single-cell spatial transcriptomics, providing further biological insights. SMAI's interpretability also enables quantification and a deeper understanding of the sources of technical confounders in single-cell data.

Uncertainty-aware single-cell annotation with a hierarchical reject option

MOTIVATION

Automatic cell type annotation methods assign cell type labels to new datasets by extracting relationships from a reference RNA-seq dataset. However, due to the limited resolution of gene expression features, there is always uncertainty present in the label assignment. To enhance the reliability and robustness of annotation, most machine learning methods address this uncertainty by providing a full reject option, i.e. when the predicted confidence score of a cell type label falls below a user-defined threshold, no label is assigned and no prediction is made. As a better alternative, some methods deploy hierarchical models and consider a so-called partial rejection by returning internal nodes of the hierarchy as label assignment. However, because a detailed experimental analysis of various rejection approaches is missing in the literature, there is currently no consensus on best practices.

RESULTS

We evaluate three annotation approaches (i) full rejection, (ii) partial rejection, and (iii) no rejection for both flat and hierarchical probabilistic classifiers. Our findings indicate that hierarchical classifiers are superior when rejection is applied, with partial rejection being the preferred rejection approach, as it preserves a significant amount of label information. For optimal rejection implementation, the rejection threshold should be determined through careful examination of a method's rejection behavior. Without rejection, flat and hierarchical annotation perform equally well, as long as the cell type hierarchy accurately captures transcriptomic relationships.

AVAILABILITY AND IMPLEMENTATION

Code is freely available at https://github.com/Latheuni/Hierarchical_reject and https://doi.org/10.5281/zenodo.10697468.

Tutorial on survival modeling with applications to omics data

MOTIVATION

Identification of genomic, molecular and clinical markers prognostic of patient survival is important for developing personalized disease prevention, diagnostic and treatment approaches. Modern omics technologies have made it possible to investigate the prognostic impact of markers at multiple molecular levels, including genomics, epigenomics, transcriptomics, proteomics and metabolomics, and how these potential risk factors complement clinical characterization of patient outcomes for survival prognosis. However, the massive sizes of the omics datasets, along with their correlation structures, pose challenges for studying relationships between the molecular information and patients' survival outcomes.

RESULTS

We present a general workflow for survival analysis that is applicable to high-dimensional omics data as inputs when identifying survival-associated features and validating survival models. In particular, we focus on the commonly used Cox-type penalized regressions and hierarchical Bayesian models for feature selection in survival analysis, which are especially useful for high-dimensional data, but the framework is applicable more generally.

AVAILABILITY AND IMPLEMENTATION

A step-by-step R tutorial using The Cancer Genome Atlas survival and omics data for the execution and evaluation of survival models has been made available at https://ocbe-uio.github.io/survomics.

Gene panel selection for targeted spatial transcriptomics

Targeted spatial transcriptomics hold particular promise in analyzing complex tissues. Most such methods, however, measure only a limited panel of transcripts, which need to be selected in advance to inform on the cell types or processes being studied. A limitation of existing gene selection methods is their reliance on scRNA-seq data, ignoring platform effects between technologies. Here we describe gpsFISH, a computational method performing gene selection through optimizing detection of known cell types. By modeling and adjusting for platform effects, gpsFISH outperforms other methods. Furthermore, gpsFISH can incorporate cell type hierarchies and custom gene preferences to accommodate diverse design requirements.