Systematic assessment of tissue dissociation and storage biases in single-cell and single-nucleus RNA-seq workflows

Spatial transcriptomics defines injury specific microenvironments and cellular interactions in kidney regeneration and disease

Kidney injury disrupts the intricate renal architecture and triggers limited regeneration, together with injury-invoked inflammation and fibrosis. Deciphering the molecular pathways and cellular interactions driving these processes is challenging due to the complex tissue structure. Here, we apply single cell spatial transcriptomics to examine ischemia-reperfusion injury in the mouse kidney. Spatial transcriptomics reveals injury-specific and spatially-dependent gene expression patterns in distinct cellular microenvironments within the kidney and predicts Clcf1-Crfl1 in a molecular interplay between persistently injured proximal tubule cells and their neighboring fibroblasts. Immune cell types play a critical role in organ repair. Spatial analysis identifies cellular microenvironments resembling early tertiary lymphoid structures and associated molecular pathways. Collectively, this study supports a focus on molecular interactions in cellular microenvironments to enhance understanding of injury, repair and disease.

Community assessment of methods to deconvolve cellular composition from bulk gene expression

We evaluate deconvolution methods, which infer levels of immune infiltration from bulk expression of tumor samples, through a community-wide DREAM Challenge. We assess six published and 22 community-contributed methods using in vitro and in silico transcriptional profiles of admixed cancer and healthy immune cells. Several published methods predict most cell types well, though they either were not trained to evaluate all functional CD8+ T cell states or do so with low accuracy. Several community-contributed methods address this gap, including a deep learning-based approach, whose strong performance establishes the applicability of this paradigm to deconvolution. Despite being developed largely using immune cells from healthy tissues, deconvolution methods predict levels of tumor-derived immune cells well. Our admixed and purified transcriptional profiles will be a valuable resource for developing deconvolution methods, including in response to common challenges we observe across methods, such as sensitive identification of functional CD4+ T cell states.

Profiling cell identity and tissue architecture with single-cell and spatial transcriptomics

Single-cell transcriptomics has broadened our understanding of cellular diversity and gene expression dynamics in healthy and diseased tissues. Recently, spatial transcriptomics has emerged as a tool to contextualize single cells in multicellular neighbourhoods and to identify spatially recurrent phenotypes, or ecotypes. These technologies have generated vast datasets with targeted-transcriptome and whole-transcriptome profiles of hundreds to millions of cells. Such data have provided new insights into developmental hierarchies, cellular plasticity and diverse tissue microenvironments, and spurred a burst of innovation in computational methods for single-cell analysis. In this Review, we discuss recent advancements, ongoing challenges and prospects in identifying and characterizing cell states and multicellular neighbourhoods. We discuss recent progress in sample processing, data integration, identification of subtle cell states, trajectory modelling, deconvolution and spatial analysis. Furthermore, we discuss the increasing application of deep learning, including foundation models, in analysing single-cell and spatial transcriptomics data. Finally, we discuss recent applications of these tools in the fields of stem cell biology, immunology, and tumour biology, and the future of single-cell and spatial transcriptomics in biological research and its translation to the clinic.

Single-nuclei sequencing of uterine serous carcinoma reveals racial differences in immune signaling

Significant racial disparities exist between Black and White patients with uterine serous carcinoma (USC). While the reasons for these disparities are unclear, several studies have demonstrated significantly different rates of driver mutations between racial groups, including TP53. However, limited research has investigated the transcriptional differences of tumors or the composition of the tumor microenvironment (TME) between these groups. Here, we report the single-nuclei RNA-sequencing profiles of primary USC tumors from diverse racial backgrounds. We find that there are significant differences between the tumors of Black and White patients. Tumors from Black patients exhibited higher expression of specific genes associated with aggressiveness, such as PAX8, and axon guidance and synaptic signaling pathways. We also demonstrated that T cell populations are reduced in the tumor tissue compared to matched benign, while anti-inflammatory macrophage populations are retained within the TME. Furthermore, we investigated the connection between PAX8 overexpression and immunosuppression in USC through regulation of several cytokines and chemokines. Notably, we show that PAX8 activity can influence macrophage gene expression and protein secretion. These studies provide a detailed understanding of the USC transcriptome and TME, and identify differences in tumor biology from patients of different racial backgrounds.

Novel Insights into Post-Myocardial Infarction Cardiac Remodeling through Algorithmic Detection of Cell-Type Composition Shifts

Background

Recent advances in single cell sequencing have led to an increased focus on the role of cell-type composition in phenotypic presentation and disease progression. Cell-type composition research in the heart is challenging due to large, frequently multinucleated cardiomyocytes that preclude most single cell approaches from obtaining accurate measurements of cell composition. Our in silico studies reveal that ignoring cell type composition when calculating differentially expressed genes (DEGs) can have significant consequences. For example, a relatively small change in cell abundance of only 10% can result in over 25% of DEGs being false positives.

Methods

We have implemented an algorithmic approach that uses snRNAseq datasets as a reference to accurately calculate cell type compositions from bulk RNAseq datasets through robust data cleaning, gene selection, and multi-sample cross-subject and cross-cell-type deconvolution. We applied our approach to cardiomyocyte-specific α1A adrenergic receptor (CM-α1A-AR) knockout mice. 8-12 week-old mice (either WT or CM-α1A-KO) were subjected to permanent left coronary artery (LCA) ligation or sham surgery (n=4 per group). Transcriptomes from the infarct border zones were collected 3 days later and analyzed using our algorithm to determine cell-type abundances, corrected differential expression calculations using DESeq2, and validated these findings using RNAscope.

Results

Uncorrected DEGs for the CM-α1A-KO X LCA interaction term featured many cell-type specific genes such as Timp4 (fibroblasts) and Aplnr (cardiomyocytes) and overall GO enrichment for terms pertaining to cardiomyocyte differentiation (P=3.1E-4). Using our algorithm, we observe a striking loss of cardiomyocytes and gain in fibroblasts in the α1A-KO + LCA mice that was not recapitulated in WT + LCA animals, although we did observe a similar increase in macrophage abundance in both conditions. This recapitulates prior results that showed a much more severe heart failure phenotype in CM-α1A-KO + LCA mice. Following correction for cell-type, our DEGs now highlight a novel set of genes enriched for GO terms such as cardiac contraction (P=3.7E-5) and actin filament organization (P=6.3E-5).

Conclusions

Our algorithm identifies and corrects for cell-type abundance in bulk RNAseq datasets opening new avenues for research on novel genes and pathways as well as an improved understanding of the role of cardiac cell types in cardiovascular disease.

Inflammatory breast cancer microenvironment repertoire based on DNA methylation data deconvolution reveals actionable targets to enhance the treatment efficacy

BACKGROUND

Although the clinical signs of inflammatory breast cancer (IBC) resemble acute inflammation, the role played by infiltrating immune and stromal cells in this aggressive disease is uncharted. The tumor microenvironment (TME) presents molecular alterations, such as epimutations, prior to morphological abnormalities. These changes affect the distribution and the intricate communication between the TME components related to cancer prognosis and therapy response. Herein, we explored the global DNA methylation profile of IBC and surrounding tissues to estimate the microenvironment cellular composition and identify epigenetically dysregulated markers.

METHODS

We used the HiTIMED algorithm to deconvolve the bulk DNA methylation data of 24 IBC and six surrounding non-tumoral tissues (SNT) (GSE238092) and determine their cellular composition. The prognostic relevance of cell types infiltrating IBC and their relationship with clinicopathological variables were investigated. CD34 (endothelial cell marker) and CD68 (macrophage marker) immunofluorescence staining was evaluated in an independent set of 17 IBC and 16 non-IBC samples.

RESULTS

We found lower infiltration of endothelial, stromal, memory B, dendritic, and natural killer cells in IBC than in SNT samples. Higher endothelial cell (EC) and stromal cell content were related to better overall survival. EC proportions positively correlated with memory B and memory CD8+ T infiltration in IBC. Immune and EC markers exhibited distinct DNA methylation profiles between IBC and SNT samples, revealing hypermethylated regions mapped to six genes (CD40, CD34, EMCN, HLA-G, PDPN, and TEK). We identified significantly higher CD34 and CD68 protein expression in IBC compared to non-IBC.

CONCLUSIONS

Our findings underscored cell subsets that distinguished patients with better survival and dysregulated markers potentially actionable through combinations of immunotherapy and epigenetic drugs.

Reciprocal regulation of T follicular helper cells and dendritic cells drives colitis development

The immunological mechanisms underlying chronic colitis are poorly understood. T follicular helper (TFH) cells are critical in helping B cells during germinal center reactions. In a T cell transfer colitis model, a lymphoid structure composed of mature dendritic cells (DCs) and TFH cells was found within T cell zones of colonic lymphoid follicles. TFH cells were required for mature DC accumulation, the formation of DC-T cell clusters and colitis development. Moreover, DCs promoted TFH cell differentiation, contributing to colitis development. A lineage-tracing analysis showed that, following migration to the lamina propria, TFH cells transdifferentiated into long-lived pathogenic TH1 cells, promoting colitis development. Our findings have therefore demonstrated the reciprocal regulation of TFH cells and DCs in colonic lymphoid follicles, which is critical in chronic colitis pathogenesis.

Concepts and new developments in droplet-based single cell multi-omics

Single cell sequencing technologies have become a fixture in the molecular profiling of cells due to their ease, flexibility, and commercial availability. In particular, partitioning individual cells inside oil droplets via microfluidic reactions enables transcriptomic or multi-omic measurements for thousands of cells in parallel. Complementing the multitude of biological discoveries from genomics analyses, the past decade has brought new capabilities from assay baselines to enable a deeper understanding of the complex data from single cell multi-omics. Here, we highlight four innovations that have improved the reliability and understanding of droplet microfluidic assays. We emphasize new developments that further orient principles of technology development and guidelines for the design, benchmarking, and implementation of new droplet-based methodologies.

Open-ST: High-resolution spatial transcriptomics in 3D

Spatial transcriptomics (ST) methods unlock molecular mechanisms underlying tissue development, homeostasis, or disease. However, there is a need for easy-to-use, high-resolution, cost-efficient, and 3D-scalable methods. Here, we report Open-ST, a sequencing-based, open-source experimental and computational resource to address these challenges and to study the molecular organization of tissues in 2D and 3D. In mouse brain, Open-ST captured transcripts at subcellular resolution and reconstructed cell types. In primary head-and-neck tumors and patient-matched healthy/metastatic lymph nodes, Open-ST captured the diversity of immune, stromal, and tumor populations in space, validated by imaging-based ST. Distinct cell states were organized around cell-cell communication hotspots in the tumor but not the metastasis. Strikingly, the 3D reconstruction and multimodal analysis of the metastatic lymph node revealed spatially contiguous structures not visible in 2D and potential biomarkers precisely at the 3D tumor/lymph node boundary. All protocols and software are available at https://rajewsky-lab.github.io/openst.

A spatiotemporal molecular atlas of mouse spinal cord injury identifies a distinct astrocyte subpopulation and therapeutic potential of IGFBP2

Spinal cord injury (SCI) triggers a cascade of intricate molecular and cellular changes that determine the outcome. In this study, we resolve the spatiotemporal organization of the injured mouse spinal cord and quantitatively assess in situ cell-cell communication following SCI. By analyzing existing single-cell RNA sequencing datasets alongside our spatial data, we delineate a subpopulation of Igfbp2-expressing astrocytes that migrate from the white matter (WM) to gray matter (GM) and become reactive upon SCI, termed Astro-GMii. Further, Igfbp2 upregulation promotes astrocyte migration, proliferation, and reactivity, and the secreted IGFBP2 protein fosters neurite outgrowth. Finally, we show that IGFBP2 significantly reduces neuronal loss and remarkably improves the functional recovery in a mouse model of SCI in vivo. Together, this study not only provides a comprehensive molecular atlas of SCI but also exemplifies how this rich resource can be applied to endow cells and genes with functional insight and therapeutic potential.

SingleCellGGM enables gene expression program identification from single-cell transcriptomes and facilitates universal cell label transfer

Gene co-expression analysis of single-cell transcriptomes, aiming to define functional relationships between genes, is challenging due to excessive dropout values. Here, we developed a single-cell graphical Gaussian model (SingleCellGGM) algorithm to conduct single-cell gene co-expression network analysis. When applied to mouse single-cell datasets, SingleCellGGM constructed networks from which gene co-expression modules with highly significant functional enrichment were identified. We considered the modules as gene expression programs (GEPs). These GEPs enable direct cell-type annotation of individual cells without cell clustering, and they are enriched with genes required for the functions of the corresponding cells, sometimes at levels greater than 10-fold. The GEPs are conserved across datasets and enable universal cell-type label transfer across different studies. We also proposed a dimension-reduction method through averaging by GEPs for single-cell analysis, enhancing the interpretability of results. Thus, SingleCellGGM offers a unique GEP-based perspective to analyze single-cell transcriptomes and reveals biological insights shared by different single-cell datasets.

SLAM/SAP signaling regulates discrete γδ T cell developmental checkpoints and shapes the innate-like γδ TCR repertoire

During thymic development, most γδ T cells acquire innate-like characteristics that are critical for their function in tumor surveillance, infectious disease, and tissue repair. The mechanisms, however, that regulate γδ T cell developmental programming remain unclear. Recently, we demonstrated that the SLAM-SAP signaling pathway regulates the development and function of multiple innate-like γδ T cell subsets. Here, we used a single-cell proteogenomics approach to identify SAP-dependent developmental checkpoints and to define the SAP-dependent γδ TCR repertoire. SAP deficiency resulted in both a significant loss of an immature Gzma + Blk + Etv5 + Tox2 + γδT17 precursor population, and a significant increase in Cd4 + Cd8+ Rorc + Ptcra + Rag1 + thymic γδ T cells. SAP-dependent diversion of embryonic day 17 thymic γδ T cell clonotypes into the αβ T cell developmental pathway was associated with a decreased frequency of mature clonotypes in neonatal thymus, and an altered γδ TCR repertoire in the periphery. Finally, we identify TRGV4/TRAV13-4(DV7)-expressing T cells as a novel, SAP-dependent Vγ4 γδT1 subset. Together, the data suggest that SAP-dependent γδ/αβ T cell lineage commitment regulates γδ T cell developmental programming and shapes the γδ TCR repertoire.

Translatome profiling reveals Itih4 as a novel smooth muscle cell-specific gene in atherosclerosis

AIMS

Vascular smooth muscle cells (SMCs) and their derivatives are key contributors to the development of atherosclerosis. However, studying changes in SMC gene expression in heterogeneous vascular tissues is challenging due to the technical limitations and high cost associated with current approaches. In this paper, we apply translating ribosome affinity purification sequencing to profile SMC-specific gene expression directly from tissue.

METHODS AND RESULTS

To facilitate SMC-specific translatome analysis, we generated SMCTRAP mice, a transgenic mouse line expressing enhanced green fluorescent protein (EGFP)-tagged ribosomal protein L10a (EGFP-L10a) under the control of the SMC-specific αSMA promoter. These mice were further crossed with the atherosclerosis model Ldlr-/-, ApoB100/100 to generate SMCTRAP-AS mice and used to profile atherosclerosis-associated SMCs in thoracic aorta samples of 15-month-old SMCTRAP and SMCTRAP-AS mice. Our analysis of SMCTRAP-AS mice showed that EGFP-L10a expression was localized to SMCs in various tissues, including the aortic wall and plaque. The TRAP fraction demonstrated high enrichment of known SMC-specific genes, confirming the specificity of our approach. We identified several genes, including Cemip, Lum, Mfge8, Spp1, and Serpina3, which are known to be involved in atherosclerosis-induced gene expression. Moreover, we identified several novel genes not previously linked to SMCs in atherosclerosis, such as Anxa4, Cd276, inter-alpha-trypsin inhibitor-4 (Itih4), Myof, Pcdh11x, Rab31, Serpinb6b, Slc35e4, Slc8a3, and Spink5. Among them, we confirmed the SMC-specific expression of Itih4 in atherosclerotic lesions using immunofluorescence staining of mouse aortic roots and spatial transcriptomics of human carotid arteries. Furthermore, our more detailed analysis of Itih4 showed its link to coronary artery disease through the colocalization of genome-wide association studies, splice quantitative trait loci (QTL), and protein QTL signals.

CONCLUSION

We generated a SMC-specific TRAP mouse line to study atherosclerosis and identified Itih4 as a novel SMC-expressed gene in atherosclerotic plaques, warranting further investigation of its putative function in extracellular matrix stability and genetic evidence of causality.

Quartet RNA reference materials improve the quality of transcriptomic data through ratio-based profiling

Certified RNA reference materials are indispensable for assessing the reliability of RNA sequencing to detect intrinsically small biological differences in clinical settings, such as molecular subtyping of diseases. As part of the Quartet Project for quality control and data integration of multi-omics profiling, we established four RNA reference materials derived from immortalized B-lymphoblastoid cell lines from four members of a monozygotic twin family. Additionally, we constructed ratio-based transcriptome-wide reference datasets between two samples, providing cross-platform and cross-laboratory 'ground truth'. Investigation of the intrinsically subtle biological differences among the Quartet samples enables sensitive assessment of cross-batch integration of transcriptomic measurements at the ratio level. The Quartet RNA reference materials, combined with the ratio-based reference datasets, can serve as unique resources for assessing and improving the quality of transcriptomic data in clinical and biological settings.

Deconstructing Intratumoral Heterogeneity through Multiomic and Multiscale Analysis of Serial Sections

Tumors may contain billions of cells, including distinct malignant clones and nonmalignant cell types. Clarifying the evolutionary histories, prevalence, and defining molecular features of these cells is essential for improving clinical outcomes, since intratumoral heterogeneity provides fuel for acquired resistance to targeted therapies. Here we present a statistically motivated strategy for deconstructing intratumoral heterogeneity through multiomic and multiscale analysis of serial tumor sections (MOMA). By combining deep sampling of IDH-mutant astrocytomas with integrative analysis of single-nucleotide variants, copy-number variants, and gene expression, we reconstruct and validate the phylogenies, spatial distributions, and transcriptional profiles of distinct malignant clones. By genotyping nuclei analyzed by single-nucleus RNA-seq for truncal mutations, we further show that commonly used algorithms for identifying cancer cells from single-cell transcriptomes may be inaccurate. We also demonstrate that correlating gene expression with tumor purity in bulk samples can reveal optimal markers of malignant cells and use this approach to identify a core set of genes that are consistently expressed by astrocytoma truncal clones, including AKR1C3, whose expression is associated with poor outcomes in several types of cancer. In summary, MOMA provides a robust and flexible strategy for precisely deconstructing intratumoral heterogeneity and clarifying the core molecular properties of distinct cellular populations in solid tumors.

Single-cell transcriptional profiling of clear cell renal cell carcinoma reveals a tumor-associated endothelial tip cell phenotype

Clear cell renal cell carcinoma (ccRCC) is the most prevalent form of renal cancer, accounting for over 75% of cases. The asymptomatic nature of the disease contributes to late-stage diagnoses and poor survival. Highly vascularized and immune infiltrated microenvironment are prominent features of ccRCC, yet the interplay between vasculature and immune cells, disease progression and response to therapy remains poorly understood. Using droplet-based single-cell RNA sequencing we profile 50,236 transcriptomes from paired tumor and healthy adjacent kidney tissues. Our analysis reveals significant heterogeneity and inter-patient variability of the tumor microenvironment. Notably, we discover a previously uncharacterized vasculature subpopulation associated with epithelial-mesenchymal transition. The cell-cell communication analysis reveals multiple modes of immunosuppressive interactions within the tumor microenvironment, including clinically relevant interactions between tumor vasculature and stromal cells with immune cells. The upregulation of the genes involved in these interactions is associated with worse survival in the TCGA KIRC cohort. Our findings demonstrate the role of tumor vasculature and stromal cell populations in shaping the ccRCC microenvironment and uncover a subpopulation of cells within the tumor vasculature that is associated with an angiogenic phenotype.

Single-nucleus RNA sequencing and lipidomics reveal characteristics of transcriptional and lipid composition in porcine longissimus dorsi muscle

BACKGROUND

Global per capita meat consumption continues to rise, especially pork. Meat quality is influenced by the content of intramuscular fat (IMF) as a key factor. The longissimus dorsi muscle of Dahe pigs (DHM, IMF: 7.98% ± 1.96%) and Dahe black pigs (DHBM, IMF: 3.30% ± 0.64%) was studied to explore cellular heterogeneity and differentially expressed genes (DEGs) associated with IMF deposition using single-nucleus RNA sequencing (snRNA-seq). The lipid composition was then analyzed using non-targeted lipidomics.

RESULTS

A total of seven cell subpopulations were identified, including myocytes, fibroblast/fibro/adipogenic progenitors (FAPs), satellite cells, endothelial cells, macrophages, pericytes, and adipocytes. Among them, FAPs and adipocytes were more focused because they could be associated with lipid deposition. 1623 DEGs in the FAPs subpopulation of DHBM were up-regulated compared with DHM, while 1535 were down-regulated. These DEGs enriched in the glycolysis/gluconeogenesis pathway. 109 DEGs were up-regulated and 806 were down-regulated in the adipocyte subpopulation of DHBM compared with DHM, which were mainly enriched in the PPAR signaling pathway and fatty acid (FA) biosynthesis. The expression level of PPARG, ABP4, LEP, and ACSL1 genes in DHM was higher than that in DHBM. Lipidomics reveals porcine lipid composition characteristics of muscle tissue. A total of 41 lipid classes and 2699 lipid species were identified in DHM and DHBM groups. The top ten relative peak areas of lipid classes in DHM and DHBM were triglyceride (TG), phosphatidylcholine (PC), phosphatidylethanolamine (PE), phosphatidylserine (PS), diglyceride (DG), cardiolipin (CL), ceramides (Cer), Simple Glc series (Hex1Cer), sphingomyelin (phSM), and phosphatidylinositol (PI). The relative peak areas of 35 lipid species in DHM were lower than DHBM, and 28 lipid species that were higher. There was a significant increase in the TG fatty acyl chains C6:0, C17:0, and C11:4, and a significant decrease in C16:0, C18:1, C18:2, and C22:4 in DHBM (p < 0.05).

CONCLUSIONS

C16:0 FA may downregulate the expression level of PPARG gene, which leads to the downregulation of fat metabolism-related genes such as ACSL, PLIN2, and FABP4 in DHBM compared with DHM. This may be the reason that the lipid deposition ability of Dahe pigs is stronger than that of Dahe black pigs, which need further investigation.

[Advancements in Single-cell RNA Sequencing Technology in the Study of the Tumor Microenvironment in Lung Cancer]

The immune microenvironment plays a key role in the development and progression of tumors. In recent years, with the rapid advancement of high-throughput sequencing technologies, researchers have gained a deeper understanding of the composition and function of immune cells in the tumor microenvironment. However, traditional bulk sequencing technologies are limited in resolving heterogeneity at the single-cell level, constraining a comprehensive understanding of the complexity of the tumor microenvironment. The advent of single-cell RNA sequencing technology has brought new opportunities to uncover the heterogeneity of the immune microenvironment in lung cancer. Currently, T-cell-centered immunotherapy in clinical settings is prone to side effects affecting prognosis, such as immunogenic drug resistance or immune-related pneumonia, with the key factor being changes in the interactions between immune cells and tumor cells in the tumor microenvironment. Single-cell RNA sequencing technology can reveal the origins and functions of different subgroups within the tumor microenvironment from perspectives such as intercellular interactions and pseudotime analysis, thereby discovering new cell subgroups or novel biomarkers, providing new avenues for uncovering resistance to immunotherapy and monitoring therapeutic efficacy. This review comprehensively discusses the newest research techniques and advancements in single-cell RNA sequencing technology for unveiling the heterogeneity of the tumor microenvironment after lung cancer immunotherapy, offering insights for enhancing the precision and personalization of immunotherapy.
.

Tackling neurodegeneration in vitro with omics: a path towards new targets and drugs

Drug discovery is a generally inefficient and capital-intensive process. For neurodegenerative diseases (NDDs), the development of novel therapeutics is particularly urgent considering the long list of late-stage drug candidate failures. Although our knowledge on the pathogenic mechanisms driving neurodegeneration is growing, additional efforts are required to achieve a better and ultimately complete understanding of the pathophysiological underpinnings of NDDs. Beyond the etiology of NDDs being heterogeneous and multifactorial, this process is further complicated by the fact that current experimental models only partially recapitulate the major phenotypes observed in humans. In such a scenario, multi-omic approaches have the potential to accelerate the identification of new or repurposed drugs against a multitude of the underlying mechanisms driving NDDs. One major advantage for the implementation of multi-omic approaches in the drug discovery process is that these overarching tools are able to disentangle disease states and model perturbations through the comprehensive characterization of distinct molecular layers (i.e., genome, transcriptome, proteome) up to a single-cell resolution. Because of recent advances increasing their affordability and scalability, the use of omics technologies to drive drug discovery is nascent, but rapidly expanding in the neuroscience field. Combined with increasingly advanced in vitro models, which particularly benefited from the introduction of human iPSCs, multi-omics are shaping a new paradigm in drug discovery for NDDs, from disease characterization to therapeutics prediction and experimental screening. In this review, we discuss examples, main advantages and open challenges in the use of multi-omic approaches for the in vitro discovery of targets and therapies against NDDs.

Single cell transcriptomic analysis of the canine duodenum in chronic inflammatory enteropathy and health

Chronic inflammatory enteropathy (CIE) is a common condition in dogs causing recurrent or persistent gastrointestinal clinical signs. Pathogenesis is thought to involve intestinal mucosal inflammatory infiltrates, but histopathological evaluation of intestinal biopsies from dogs with CIE fails to guide treatment, inform prognosis, or correlate with clinical remission. We employed single-cell RNA sequencing to catalog and compare the diversity of cells present in duodenal mucosal endoscopic biopsies from 3 healthy dogs and 4 dogs with CIE. Through characterization of 35,668 cells, we identified 31 transcriptomically distinct cell populations, including T cells, epithelial cells, and myeloid cells. Both healthy and CIE samples contributed to each cell population. T cells were broadly subdivided into GZMAhigh (putatively annotated as tissue resident) and IL7Rhigh (putatively annotated as non-resident) T cell categories, with evidence of a skewed proportion favoring an increase in the relative proportion of IL7Rhigh T cells in CIE dogs. Among the myeloid cells, neutrophils from CIE samples exhibited inflammatory (SOD2 and IL1A) gene expression signatures. Numerous differentially expressed genes were identified in epithelial cells, with gene set enrichment analysis suggesting enterocytes from CIE dogs may be undergoing stress responses and have altered metabolic properties. Overall, this work reveals the previously unappreciated cellular heterogeneity in canine duodenal mucosa and provides new insights into molecular mechanisms which may contribute to intestinal dysfunction in CIE. The cell type gene signatures developed through this study may also be used to better understand the subtleties of canine intestinal physiology in health and disease.

From whole bodies to single cells: A guide to transcriptomic approaches for ecology and evolutionary biology

RNA sequencing (RNAseq) methodology has experienced a burst of technological developments in the last decade, which has opened up opportunities for studying the mechanisms of adaptation to environmental factors at both the organismal and cellular level. Selecting the most suitable experimental approach for specific research questions and model systems can, however, be a challenge and researchers in ecology and evolution are commonly faced with the choice of whether to study gene expression variation in whole bodies, specific tissues, and/or single cells. A wide range of sometimes polarised opinions exists over which approach is best. Here, we highlight the advantages and disadvantages of each of these approaches to provide a guide to help researchers make informed decisions and maximise the power of their study. Using illustrative examples of various ecological and evolutionary research questions, we guide the readers through the different RNAseq approaches and help them identify the most suitable design for their own projects.

Simultaneous isolation of intact brain cells and cell-specific extracellular vesicles from cryopreserved Alzheimer's disease cortex

BACKGROUND

The neuronal and gliaI populations within the brain are tightly interwoven, making isolation and study of large populations of a single cell type from brain tissue a major technical challenge. Concurrently, cell-type specific extracellular vesicles (EVs) hold enormous diagnostic and therapeutic potential in neurodegenerative disorders including Alzheimer's disease (AD).

NEW METHOD

Postmortem AD cortical samples were thawed and gently dissociated. Following filtration, myelin and red blood cell removal, cell pellets were immunolabeled with fluorescent antibodies and analyzed by flow cytometry. The cell pellet supernatant was applied to a triple sucrose cushion for brain EV isolation.

RESULTS

Neuronal, astrocyte and microglial cell populations were identified. Cell integrity was demonstrated using calcein AM, which is retained by cells with esterase activity and an intact membrane. For some experiments cell pellets were fixed, permeabilized, and immunolabeled for cell-specific markers. Characterization of brain small EV fractions showed the expected size, depletion of EV negative markers, and enrichment in positive and cell-type specific markers.

COMPARISON WITH EXISTING METHODS AND CONCLUSIONS

We optimized and integrated established protocols, aiming to maximize information obtained from each human autopsy brain sample. The uniqueness of our method lies in its capability to isolate cells and EVs from a single cryopreserved brain sample. Our results not only demonstrate the feasibility of isolating specific brain cell subpopulations for RNA-seq but also validate these subpopulations at the protein level. The accelerated study of EVs from human samples is crucial for a better understanding of their contribution to neuron/glial crosstalk and disease progression.

Gene representation bias in spatial transcriptomics

For sequencing-based spatial transcriptomics data, the gene-spot count matrix is highly sparse. This feature is similar to scRNA-seq. The goal of this paper is to identify whether there exist genes that are frequently under-detected in Visium compared to bulk RNA-seq, and the underlying potential mechanism of under-detection in Visium. We collected paired Visium and bulk RNA-seq data for 28 human samples and 19 mouse samples, which covered diverse tissue sources. We compared the two data types and observed that there indeed exists a collection of genes frequently under-detected in Visium compared to bulk RNA-seq. We performed a motif search to examine the last 350 bp of the frequently under-detected genes, and we observed that the poly (T) motif was significantly enriched in genes identified from both human and mouse data, which matches with our previous finding about frequently under-detected genes in scRNA-seq. We hypothesized that the poly (T) motif may be able to form a hairpin structure with the poly (A) tails of their mRNA transcripts, making it difficult for their mRNA transcripts to be captured during Visium library preparation.

Single-cell transcriptomics reveals antigen-presenting capacity and therapeutic resistance potential of immunomodulatory endothelial cells in colorectal cancer

BACKGROUND

The heterogeneity of tumor endothelial cells (TECs) hinders the efficacy of antiangiogenic therapies (AATs). Only a small percentage of angiogenic TECs are considered effective targets for AATs. Immunomodulatory ECs (IMECs), as a newly focused functional subgroup of endothelial cells (ECs), are being evaluated for their ability to regulate tumor immune balance and influence existing AATs.

METHODS

Based on single-cell transcriptome data from colorectal cancer in a publicly available database, we conducted a wide array of bioinformatic approaches to study EC subsets that meet the IMECs definition. Our investigation encompassed the gene expression signatures of these subsets, cellular composition differences, cell-cell interactions.

RESULTS

Two subsets that meet the IMECs definition were found in tumors and para-cancerous tissues. Combined with the results of gene ontological analysis and interaction with CD4+ T cells, we found that IMECs can present MHC-II antigens to mature CD4+ T cells. There were differences in the level of interaction between IMECs and different types of mature CD4+ T cell subsets. In addition, IMEC subsets had different expression levels of angiogenesis related genes. The angiogenesis score of IMECs decreased after patients received immunotherapy. IMEC subsets do not depend on a single proangiogenic receptor and are involved in regulating angiogenesis, which may reduce the efficacy of AATs. The adverse effects of specific IMEC subsets on AATs were validated in the RNA-seq dataset of the bevacizumab treatment group.

CONCLUSION

Our study suggests the potential MHC-II antigen presentation capacity of IMECs and the enhanced angiogenesis characteristics within tumors. The function of IMECs in the vascular network may have a potentially adverse effect on AATs. Controlling the functional properties of IMECs may be a new angle for tumor therapy.

A novel approach to generate enzyme-free single cell suspensions from archived tissues for miRNA sequencing

Obtaining high-quality omics data at the single-cell level from archived human tissue samples is crucial for gaining insights into cellular heterogeneity and pushing the field of personalized medicine forward. In this technical brief we present a comprehensive methodological framework for the efficient enzyme-free preparation of tissue-derived single cell suspensions and their conversion into single-cell miRNA sequencing libraries. The resulting data from this study have the potential to deepen our understanding of miRNA expression at the single-cell level and its relevance in the context of the examined tissues. The workflow encompasses tissue collection, RNALater immersion, storage, thawing, TissueGrinder-mediated dissociation, miRNA lysis, library preparation, sequencing, and data analysis. Quality control measures ensure reliable miRNA data, with specific attention to sample quality. The UMAP analysis reveals tissue-specific cell clustering, while miRNA diversity reflects tissue variations. The presented workflow effectively processes preserved tissues, extending opportunities for retrospective analysis and biobank utilization.

Stress-free single-cell transcriptomic profiling and functional genomics of murine eosinophils

Eosinophils are a class of granulocytes with pleiotropic functions in homeostasis and various human diseases. Nevertheless, they are absent from conventional single-cell RNA sequencing atlases owing to technical difficulties preventing their transcriptomic interrogation. Consequently, eosinophil heterogeneity and the gene regulatory networks underpinning their diverse functions remain poorly understood. We have developed a stress-free protocol for single-cell RNA capture from murine tissue-resident eosinophils, which revealed distinct intestinal subsets and their roles in colitis. Here we describe in detail how to enrich eosinophils from multiple tissues of residence and how to capture high-quality single-cell transcriptomes by preventing transcript degradation. By combining magnetic eosinophil enrichment with microwell-based single-cell RNA capture (BD Rhapsody), our approach minimizes shear stress and processing time. Moreover, we report how to perform genome-wide CRISPR pooled genetic screening in ex vivo-conditioned bone marrow-derived eosinophils to functionally probe pathways required for their differentiation and intestinal maturation. These protocols can be performed by any researcher with basic skills in molecular biology and flow cytometry, and can be adapted to investigate other granulocytes, such as neutrophils and mast cells, thereby offering potential insights into their roles in both homeostasis and disease pathogenesis. Single-cell transcriptomics of eosinophils can be performed in 2-3 d, while functional genomics assays may require up to 1 month.

Multi-omics profiling of mouse polycystic kidney disease progression at a single cell resolution

Autosomal dominant polycystic kidney disease (ADPKD) is the most common hereditary kidney disease and causes significant morbidity, ultimately leading to end-stage kidney disease. PKD pathogenesis is characterized by complex and dynamic alterations in multiple cell types during disease progression, hampering a deeper understanding of disease mechanism and the development of therapeutic approaches. Here, we generate a single nucleus multimodal atlas of an orthologous mouse PKD model at early, mid and late timepoints, consisting of 125,434 single-nucleus transcriptomic and epigenetic multiomes. We catalogue differentially expressed genes and activated epigenetic regions in each cell type during PKD progression, characterizing cell-type-specific responses to Pkd1 deletion. We describe heterogeneous, atypical collecting duct cells as well as proximal tubular cells that constitute cyst epithelia in PKD. The transcriptional regulation of the cyst lining cell marker GPRC5A is conserved between mouse and human PKD cystic epithelia, suggesting shared gene regulatory pathways. Our single nucleus multiomic analysis of mouse PKD provides a foundation to understand the earliest changes molecular deregulation in a mouse model of PKD at a single-cell resolution.

CellsFromSpace: a fast, accurate, and reference-free tool to deconvolve and annotate spatially distributed omics data

Motivation

Spatial transcriptomics enables the analysis of cell crosstalk in healthy and diseased organs by capturing the transcriptomic profiles of millions of cells within their spatial contexts. However, spatial transcriptomics approaches also raise new computational challenges for the multidimensional data analysis associated with spatial coordinates.

Results

In this context, we introduce a novel analytical framework called CellsFromSpace based on independent component analysis (ICA), which allows users to analyze various commercially available technologies without relying on a single-cell reference dataset. The ICA approach deployed in CellsFromSpace decomposes spatial transcriptomics data into interpretable components associated with distinct cell types or activities. ICA also enables noise or artifact reduction and subset analysis of cell types of interest through component selection. We demonstrate the flexibility and performance of CellsFromSpace using real-world samples to demonstrate ICA's ability to successfully identify spatially distributed cells as well as rare diffuse cells, and quantitatively deconvolute datasets from the Visium, Slide-seq, MERSCOPE, and CosMX technologies. Comparative analysis with a current alternative reference-free deconvolution tool also highlights CellsFromSpace's speed, scalability and accuracy in processing complex, even multisample datasets. CellsFromSpace also offers a user-friendly graphical interface enabling non-bioinformaticians to annotate and interpret components based on spatial distribution and contributor genes, and perform full downstream analysis.

Availability and implementation

CellsFromSpace (CFS) is distributed as an R package available from github at https://github.com/gustaveroussy/CFS along with tutorials, examples, and detailed documentation.

Advances in long-read single-cell transcriptomics

Long-read single-cell transcriptomics (scRNA-Seq) is revolutionizing the way we profile heterogeneity in disease. Traditional short-read scRNA-Seq methods are limited in their ability to provide complete transcript coverage, resolve isoforms, and identify novel transcripts. The scRNA-Seq protocols developed for long-read sequencing platforms overcome these limitations by enabling the characterization of full-length transcripts. Long-read scRNA-Seq techniques initially suffered from comparatively poor accuracy compared to short read scRNA-Seq. However, with improvements in accuracy, accessibility, and cost efficiency, long-reads are gaining popularity in the field of scRNA-Seq. This review details the advances in long-read scRNA-Seq, with an emphasis on library preparation protocols and downstream bioinformatics analysis tools.

Evaluating performance and applications of sample-wise cell deconvolution methods on human brain transcriptomic data

Sample-wise deconvolution methods estimate cell-type proportions and gene expressions in bulk tissue samples, yet their performance and biological applications remain unexplored, particularly in human brain transcriptomic data. Here, nine deconvolution methods were evaluated with sample-matched data from bulk tissue RNA sequencing (RNA-seq), single-cell/nuclei (sc/sn) RNA-seq, and immunohistochemistry. A total of 1,130,767 nuclei per cells from 149 adult postmortem brains and 72 organoid samples were used. The results showed the best performance of dtangle for estimating cell proportions and bMIND for estimating sample-wise cell-type gene expressions. For eight brain cell types, 25,273 cell-type eQTLs were identified with deconvoluted expressions (decon-eQTLs). The results showed that decon-eQTLs explained more schizophrenia GWAS heritability than bulk tissue or single-cell eQTLs did alone. Differential gene expressions associated with Alzheimer's disease, schizophrenia, and brain development were also examined using the deconvoluted data. Our findings, which were replicated in bulk tissue and single-cell data, provided insights into the biological applications of deconvoluted data in multiple brain disorders.

Fourteen years of cellular deconvolution: methodology, applications, technical evaluation and outstanding challenges

Single-cell RNA sequencing (scRNA-Seq) is a recent technology that allows for the measurement of the expression of all genes in each individual cell contained in a sample. Information at the single-cell level has been shown to be extremely useful in many areas. However, performing single-cell experiments is expensive. Although cellular deconvolution cannot provide the same comprehensive information as single-cell experiments, it can extract cell-type information from bulk RNA data, and therefore it allows researchers to conduct studies at cell-type resolution from existing bulk datasets. For these reasons, a great effort has been made to develop such methods for cellular deconvolution. The large number of methods available, the requirement of coding skills, inadequate documentation, and lack of performance assessment all make it extremely difficult for life scientists to choose a suitable method for their experiment. This paper aims to fill this gap by providing a comprehensive review of 53 deconvolution methods regarding their methodology, applications, performance, and outstanding challenges. More importantly, the article presents a benchmarking of all these 53 methods using 283 cell types from 30 tissues of 63 individuals. We also provide an R package named DeconBenchmark that allows readers to execute and benchmark the reviewed methods (https://github.com/tinnlab/DeconBenchmark).

Spatial detection of mitochondrial DNA and RNA in tissues

Background

Mitochondrial health has gained attention in a number of diseases, both as an indicator of disease state and as a potential therapeutic target. The quality and amount of mitochondrial DNA (mtDNA) and RNA (mtRNA) can be important indicators of mitochondrial and cell health, but are difficult to measure in complex tissues.

Methods

mtDNA and mtRNA in zebrafish retina samples were fluorescently labeled using RNAscope™ in situ hybridization, then mitochondria were stained using immunohistochemistry. Pretreatment with RNase was used for validation. Confocal images were collected and analyzed, and relative amounts of mtDNA and mtRNA were reported. Findings regarding mtDNA were confirmed using qPCR.

Results

Signals from probes detecting mtDNA and mtRNA were localized to mitochondria, and were differentially sensitive to RNase. This labeling strategy allows for quantification of relative mtDNA and mtRNA levels in individual cells. As a demonstration of the method in a complex tissue, single photoreceptors in zebrafish retina were analyzed for mtDNA and mtRNA content. An increase in mtRNA but not mtDNA coincides with proliferation of mitochondria at night in cones. A similar trend was measured in rods.

Discussion

Mitochondrial gene expression is an important component of cell adaptations to disease, stress, or aging. This method enables the study of mtDNA and mtRNA in single cells of an intact, complex tissue. The protocol presented here uses commercially-available tools, and is adaptable to a range of species and tissue types.

Advanced sequencing-based high-throughput and long-read single-cell transcriptome analysis

Cells are the fundamental building blocks of living systems, exhibiting significant heterogeneity. The transcriptome connects the cellular genotype and phenotype, and profiling single-cell transcriptomes is critical for uncovering distinct cell types, states, and the interplay between cells in development, health, and disease. Nevertheless, single-cell transcriptome analysis faces daunting challenges due to the low abundance and diverse nature of RNAs in individual cells, as well as their heterogeneous expression. The advent and continuous advancements of next-generation sequencing (NGS) and third-generation sequencing (TGS) technologies have solved these problems and facilitated the high-throughput, sensitive, full-length, and rapid profiling of single-cell RNAs. In this review, we provide a broad introduction to current methodologies for single-cell transcriptome sequencing. First, state-of-the-art advancements in high-throughput and full-length single-cell RNA sequencing (scRNA-seq) platforms using NGS are reviewed. Next, TGS-based long-read scRNA-seq methods are summarized. Finally, a brief conclusion and perspectives for comprehensive single-cell transcriptome analysis are discussed.

Spatiotemporal analysis of gene expression in the human dentate gyrus reveals age-associated changes in cellular maturation and neuroinflammation

The dentate gyrus of the hippocampus is important for many cognitive functions, including learning, memory, and mood. Here, we investigated age-associated changes in transcriptome-wide spatial gene expression in the human dentate gyrus across the lifespan. Genes associated with neurogenesis and the extracellular matrix were enriched in infants, while gene markers of inhibitory neurons and cell proliferation showed increases and decreases in post-infancy, respectively. While we did not find evidence for neural proliferation post-infancy, we did identify molecular signatures supporting protracted maturation of granule cells. We also identified a wide-spread hippocampal aging signature and an age-associated increase in genes related to neuroinflammation. Our findings suggest major changes to the putative neurogenic niche after infancy and identify molecular foci of brain aging in glial and neuropil enriched tissue.

Intra-tumoral T cells in pediatric brain tumors display clonal expansion and effector properties

Brain tumors in children are a devastating disease in a high proportion of patients. Owing to inconsistent results in clinical trials in unstratified patients, the role of immunotherapy remains unclear. We performed an in-depth survey of the single-cell transcriptomes and clonal relationship of intra-tumoral T cells from children with brain tumors. Our results demonstrate that a large fraction of T cells in the tumor tissue are clonally expanded with the potential to recognize tumor antigens. Such clonally expanded T cells display enrichment of transcripts linked to effector function, tissue residency, immune checkpoints and signatures of neoantigen-specific T cells and immunotherapy response. We identify neoantigens in pediatric brain tumors and show that neoantigen-specific T cell gene signatures are linked to better survival outcomes. Notably, among the patients in our cohort, we observe substantial heterogeneity in the degree of clonal expansion and magnitude of T cell response. Our findings suggest that characterization of intra-tumoral T cell responses may enable selection of patients for immunotherapy, an approach that requires prospective validation in clinical trials.

Single-cell and single-nucleus RNA-sequencing from paired normal-adenocarcinoma lung samples provide both common and discordant biological insights

Whether single-cell RNA-sequencing (scRNA-seq) captures the same biological information as single-nucleus RNA-sequencing (snRNA-seq) remains uncertain and likely to be context-dependent. Herein, a head-to-head comparison was performed in matched normal-adenocarcinoma human lung samples to assess biological insights derived from scRNA-seq versus snRNA-seq and better understand the cellular transition that occurs from normal to tumoral tissue. Here, the transcriptome of 160,621 cells/nuclei was obtained. In non-tumor lung, cell type proportions varied widely between scRNA-seq and snRNA-seq with a predominance of immune cells in the former (81.5%) and epithelial cells (69.9%) in the later. Similar results were observed in adenocarcinomas, in addition to an overall increase in cell type heterogeneity and a greater prevalence of copy number variants in cells of epithelial origin, which suggests malignant assignment. The cell type transition that occurs from normal lung tissue to adenocarcinoma was not always concordant whether cells or nuclei were examined. As expected, large differential expression of the whole-cell and nuclear transcriptome was observed, but cell-type specific changes of paired normal and tumor lung samples revealed a set of common genes in the cells and nuclei involved in cancer-related pathways. In addition, we showed that the ligand-receptor interactome landscape of lung adenocarcinoma was largely different whether cells or nuclei were evaluated. Immune cell depletion in fresh specimens partly mitigated the difference in cell type composition observed between cells and nuclei. However, the extra manipulations affected cell viability and amplified the transcriptional signatures associated with stress responses. In conclusion, research applications focussing on mapping the immune landscape of lung adenocarcinoma benefit from scRNA-seq in fresh samples, whereas snRNA-seq of frozen samples provide a low-cost alternative to profile more epithelial and cancer cells, and yield cell type proportions that more closely match tissue content.

NKX2-2 based nuclei sorting on frozen human archival pancreas enables the enrichment of islet endocrine populations for single-nucleus RNA sequencing

BACKGROUND

Current approaches to profile the single-cell transcriptomics of human pancreatic endocrine cells almost exclusively rely on freshly isolated islets. However, human islets are limited in availability. Furthermore, the extensive processing steps during islet isolation and subsequent single cell dissolution might alter gene expressions. In this work, we report the development of a single-nucleus RNA sequencing (snRNA-seq) approach with targeted islet cell enrichment for endocrine-population focused transcriptomic profiling using frozen archival pancreatic tissues without islet isolation.

RESULTS

We cross-compared five nuclei isolation protocols and selected the citric acid method as the best strategy to isolate nuclei with high RNA integrity and low cytoplasmic contamination from frozen archival human pancreata. We innovated fluorescence-activated nuclei sorting based on the positive signal of NKX2-2 antibody to enrich nuclei of the endocrine population from the entire nuclei pool of the pancreas. Our sample preparation procedure generated high-quality single-nucleus gene-expression libraries while preserving the endocrine population diversity. In comparison with single-cell RNA sequencing (scRNA-seq) library generated with live cells from freshly isolated human islets, the snRNA-seq library displayed comparable endocrine cellular composition and cell type signature gene expression. However, between these two types of libraries, differential enrichments of transcripts belonging to different functional classes could be observed.

CONCLUSIONS

Our work fills a technological gap and helps to unleash frozen archival pancreatic tissues for molecular profiling targeting the endocrine population. This study opens doors to retrospective mappings of endocrine cell dynamics in pancreatic tissues of complex histopathology. We expect that our protocol is applicable to enrich nuclei for transcriptomics studies from various populations in different types of frozen archival tissues.

Immune checkpoint therapy responders display early clonal expansion of tumor infiltrating lymphocytes

Immune checkpoint therapy (ICT) causes durable tumour responses in a subgroup of patients, but it is not well known how T cell receptor beta (TCRβ) repertoire dynamics contribute to the therapeutic response. Using murine models that exclude variation in host genetics, environmental factors and tumour mutation burden, limiting variation between animals to naturally diverse TCRβ repertoires, we applied TCRseq, single cell RNAseq and flow cytometry to study TCRβ repertoire dynamics in ICT responders and non-responders. Increased oligoclonal expansion of TCRβ clonotypes was observed in responding tumours. Machine learning identified TCRβ CDR3 signatures unique to each tumour model, and signatures associated with ICT response at various timepoints before or during ICT. Clonally expanded CD8+ T cells in responding tumours post ICT displayed effector T cell gene signatures and phenotype. An early burst of clonal expansion during ICT is associated with response, and we report unique dynamics in TCRβ signatures associated with ICT response.

Donor regulatory T cells rapidly adapt to recipient tissues to control murine acute graft-versus-host disease

The adoptive transfer of regulatory T cells is a promising strategy to prevent graft-versus-host disease after allogeneic bone marrow transplantation. Here, we use a major histocompatibility complex-mismatched mouse model to follow the fate of in vitro expanded donor regulatory T cells upon migration to target organs. Employing comprehensive gene expression and repertoire profiling, we show that they retain their suppressive function and plasticity after transfer. Upon entering non-lymphoid tissues, donor regulatory T cells acquire organ-specific gene expression profiles resembling tissue-resident cells and activate hallmark suppressive and cytotoxic pathways, most evidently in the colon, when co-transplanted with graft-versus-host disease-inducing conventional T cells. Dominant T cell receptor clonotypes overlap between organs and across recipients and their relative abundance correlates with protection efficacy. Thus, this study reveals donor regulatory T cell selection and adaptation mechanisms in target organs and highlights protective features of Treg to guide the development of improved graft-versus-host disease prevention strategies.

Benchmark of cellular deconvolution methods using a multi-assay reference dataset from postmortem human prefrontal cortex

Background

Cellular deconvolution of bulk RNA-sequencing (RNA-seq) data using single cell or nuclei RNA-seq (sc/snRNA-seq) reference data is an important strategy for estimating cell type composition in heterogeneous tissues, such as human brain. Computational methods for deconvolution have been developed and benchmarked against simulated data, pseudobulked sc/snRNA-seq data, or immunohistochemistry reference data. A major limitation in developing improved deconvolution algorithms has been the lack of integrated datasets with orthogonal measurements of gene expression and estimates of cell type proportions on the same tissue sample. Deconvolution algorithm performance has not yet been evaluated across different RNA extraction methods (cytosolic, nuclear, or whole cell RNA), different library preparation types (mRNA enrichment vs. ribosomal RNA depletion), or with matched single cell reference datasets.

Results

A rich multi-assay dataset was generated in postmortem human dorsolateral prefrontal cortex (DLPFC) from 22 tissue blocks. Assays included spatially-resolved transcriptomics, snRNA-seq, bulk RNA-seq (across six library/extraction RNA-seq combinations), and RNAScope/Immunofluorescence (RNAScope/IF) for six broad cell types. The Mean Ratio method, implemented in the DeconvoBuddies R package, was developed for selecting cell type marker genes. Six computational deconvolution algorithms were evaluated in DLPFC and predicted cell type proportions were compared to orthogonal RNAScope/IF measurements.

Conclusions

Bisque and hspe were the most accurate methods, were robust to differences in RNA library types and extractions. This multi-assay dataset showed that cell size differences, marker genes differentially quantified across RNA libraries, and cell composition variability in reference snRNA-seq impact the accuracy of current deconvolution methods.

lute: estimating the cell composition of heterogeneous tissue with varying cell sizes using gene expression

Relative cell type fraction estimates in bulk RNA-sequencing data are important to control for cell composition differences across heterogenous tissue samples. Current computational tools estimate relative RNA abundances rather than cell type proportions in tissues with varying cell sizes, leading to biased estimates. We present lute, a computational tool to accurately deconvolute cell types with varying sizes. Our software wraps existing deconvolution algorithms in a standardized framework. Using simulated and real datasets, we demonstrate how lute adjusts for differences in cell sizes to improve the accuracy of cell composition. Software is available from https://bioconductor.org/packages/lute.

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.

Advancements in Single-Cell RNA Sequencing Research for Neurological Diseases

Neurological diseases are a major cause of the global burden of disease. Although the mechanisms of the occurrence and development of neurological diseases are not fully clear, most of them are associated with cells mediating neuroinflammation. Yet medications and other therapeutic options to improve treatment are still very limited. Single-cell RNA sequencing (scRNA-seq), as a delightfully potent breakthrough technology, not only identifies various cell types and response states but also uncovers cell-specific gene expression changes, gene regulatory networks, intercellular communication, and cellular movement trajectories, among others, in different cell types. In this review, we describe the technology of scRNA-seq in detail and discuss and summarize the application of scRNA-seq in exploring neurological diseases, elaborating the corresponding specific mechanisms of the diseases as well as providing a reliable basis for new therapeutic approaches. Finally, we affirm that scRNA-seq promotes the development of the neuroscience field and enables us to have a deeper cellular understanding of neurological diseases in the future, which provides strong support for the treatment of neurological diseases and the improvement of patients' prognosis.

Bioinformatics in urology - molecular characterization of pathophysiology and response to treatment

The application of bioinformatics has revolutionized the practice of medicine in the past 20 years. From early studies that uncovered subtypes of cancer to broad efforts spearheaded by the Cancer Genome Atlas initiative, the use of bioinformatics strategies to analyse high-dimensional data has provided unprecedented insights into the molecular basis of disease. In addition to the identification of disease subtypes - which enables risk stratification - informatics analysis has facilitated the identification of novel risk factors and drivers of disease, biomarkers of progression and treatment response, as well as possibilities for drug repurposing or repositioning; moreover, bioinformatics has guided research towards precision and personalized medicine. Implementation of specific computational approaches such as artificial intelligence, machine learning and molecular subtyping has yet to become widespread in urology clinical practice for reasons of cost, disruption of clinical workflow and need for prospective validation of informatics approaches in independent patient cohorts. Solving these challenges might accelerate routine integration of bioinformatics into clinical settings.

A single-cell atlas enables mapping of homeostatic cellular shifts in the adult human breast

Here we use single-cell RNA sequencing to compile a human breast cell atlas assembled from 55 donors that had undergone reduction mammoplasties or risk reduction mastectomies. From more than 800,000 cells we identified 41 cell subclusters across the epithelial, immune and stromal compartments. The contribution of these different clusters varied according to the natural history of the tissue. Age, parity and germline mutations, known to modulate the risk of developing breast cancer, affected the homeostatic cellular state of the breast in different ways. We found that immune cells from BRCA1 or BRCA2 carriers had a distinct gene expression signature indicative of potential immune exhaustion, which was validated by immunohistochemistry. This suggests that immune-escape mechanisms could manifest in non-cancerous tissues very early during tumor initiation. This atlas is a rich resource that can be used to inform novel approaches for early detection and prevention of breast cancer.

FixNCut: single-cell genomics through reversible tissue fixation and dissociation

The use of single-cell technologies for clinical applications requires disconnecting sampling from downstream processing steps. Early sample preservation can further increase robustness and reproducibility by avoiding artifacts introduced during specimen handling. We present FixNCut, a methodology for the reversible fixation of tissue followed by dissociation that overcomes current limitations. We applied FixNCut to human and mouse tissues to demonstrate the preservation of RNA integrity, sequencing library complexity, and cellular composition, while diminishing stress-related artifacts. Besides single-cell RNA sequencing, FixNCut is compatible with multiple single-cell and spatial technologies, making it a versatile tool for robust and flexible study designs.

CATD: a reproducible pipeline for selecting cell-type deconvolution methods across tissues

Motivation

Cell-type deconvolution methods aim to infer cell composition from bulk transcriptomic data. The proliferation of developed methods coupled with inconsistent results obtained in many cases, highlights the pressing need for guidance in the selection of appropriate methods. Additionally, the growing accessibility of single-cell RNA sequencing datasets, often accompanied by bulk expression from related samples enable the benchmark of existing methods.

Results

In this study, we conduct a comprehensive assessment of 31 methods, utilizing single-cell RNA-sequencing data from diverse human and mouse tissues. Employing various simulation scenarios, we reveal the efficacy of regression-based deconvolution methods, highlighting their sensitivity to reference choices. We investigate the impact of bulk-reference differences, incorporating variables such as sample, study and technology. We provide validation using a gold standard dataset from mononuclear cells and suggest a consensus prediction of proportions when ground truth is not available. We validated the consensus method on data from the stomach and studied its spillover effect. Importantly, we propose the use of the critical assessment of transcriptomic deconvolution (CATD) pipeline which encompasses functionalities for generating references and pseudo-bulks and running implemented deconvolution methods. CATD streamlines simultaneous deconvolution of numerous bulk samples, providing a practical solution for speeding up the evaluation of newly developed methods.

Availability and implementation

https://github.com/Papatheodorou-Group/CATD_snakemake.

Quality Control for Single Cell Analysis of High-plex Tissue Profiles using CyLinter

Tumors are complex assemblies of cellular and acellular structures patterned on spatial scales from microns to centimeters. Study of these assemblies has advanced dramatically with the introduction of high-plex spatial profiling. Image-based profiling methods reveal the intensities and spatial distributions of 20-100 proteins at subcellular resolution in 103-107 cells per specimen. Despite extensive work on methods for extracting single-cell data from these images, all tissue images contain artefacts such as folds, debris, antibody aggregates, optical aberrations and image processing errors that arise from imperfections in specimen preparation, data acquisition, image assembly, and feature extraction. We show that these artefacts dramatically impact single-cell data analysis, obscuring meaningful biological interpretation. We describe an interactive quality control software tool, CyLinter, that identifies and removes data associated with imaging artefacts. CyLinter greatly improves single-cell analysis, especially for archival specimens sectioned many years prior to data collection, such as those from clinical trials.

Recent advances in differential expression analysis for single-cell RNA-seq and spatially resolved transcriptomic studies

Differential expression (DE) analysis is a necessary step in the analysis of single-cell RNA sequencing (scRNA-seq) and spatially resolved transcriptomics (SRT) data. Unlike traditional bulk RNA-seq, DE analysis for scRNA-seq or SRT data has unique characteristics that may contribute to the difficulty of detecting DE genes. However, the plethora of DE tools that work with various assumptions makes it difficult to choose an appropriate one. Furthermore, a comprehensive review on detecting DE genes for scRNA-seq data or SRT data from multi-condition, multi-sample experimental designs is lacking. To bridge such a gap, here, we first focus on the challenges of DE detection, then highlight potential opportunities that facilitate further progress in scRNA-seq or SRT analysis, and finally provide insights and guidance in selecting appropriate DE tools or developing new computational DE methods.

Deconstructing intratumoral heterogeneity through multiomic and multiscale analysis of serial sections

Tumors may contain billions of cells including distinct malignant clones and nonmalignant cell types. Clarifying the evolutionary histories, prevalence, and defining molecular features of these cells is essential for improving clinical outcomes, since intratumoral heterogeneity provides fuel for acquired resistance to targeted therapies. Here we present a statistically motivated strategy for deconstructing intratumoral heterogeneity through multiomic and multiscale analysis of serial tumor sections (MOMA). By combining deep sampling of IDH-mutant astrocytomas with integrative analysis of single-nucleotide variants, copy-number variants, and gene expression, we reconstruct and validate the phylogenies, spatial distributions, and transcriptional profiles of distinct malignant clones. By genotyping nuclei analyzed by single-nucleus RNA-seq for truncal mutations, we further show that commonly used algorithms for identifying cancer cells from single-cell transcriptomes may be inaccurate. We also demonstrate that correlating gene expression with tumor purity in bulk samples can reveal optimal markers of malignant cells and use this approach to identify a core set of genes that is consistently expressed by astrocytoma truncal clones, including AKR1C3, whose expression is associated with poor outcomes in several types of cancer. In summary, MOMA provides a robust and flexible strategy for precisely deconstructing intratumoral heterogeneity and clarifying the core molecular properties of distinct cellular populations in solid tumors.