Exploring the immune escape mechanisms in gastric cancer patients based on the deep AI algorithms and single-cell sequencing analysis.

Gastric cancer is a prevalent and deadly malignancy, and the response to immunotherapy varies among patients. This study aimed to develop a prognostic model for gastric cancer patients and investigate immune escape mechanisms using deep machine learning and single-cell sequencing analysis. Data from public databases were analysed, and a prediction model was constructed using 101 algorithms. The high-AIDPS group, characterized by increased AIDPS expression, exhibited worse survival, genomic variations and immune cell infiltration. These patients also showed immunotherapy tolerance. Treatment strategies targeting the high-AIDPS group identified three potential drugs. Additionally, distinct cluster groups and upregulated AIDPS-associated genes were observed in gastric adenocarcinoma cell lines. Inhibition of GHRL expression suppressed cancer cell activity, inhibited M2 polarization in macrophages and reduced invasiveness. Overall, AIDPS plays a critical role in gastric cancer prognosis, genomic variations, immune cell infiltration and immunotherapy response, and targeting GHRL expression holds promise for personalized treatment. These findings contribute to improved clinical management in gastric cancer.

In situ single-cell therapeutic response imaging facilitated by the TRIPODD fluorescence imaging platform.

Purpose: Small molecule drugs such as tyrosine kinase inhibitors (TKIs) targeting tumoral molecular dependencies have become standard of care for numerous cancer types. Notably, epidermal growth factor receptor (EGFR) TKIs (e.g., erlotinib, afatinib, osimertinib) are the current first-line treatment for non-small cell lung cancer (NSCLC) due to their improved therapeutic outcomes for EGFR mutated and overexpressing disease over traditional platinum-based chemotherapy. However, many NSCLC tumors develop resistance to EGFR TKI therapy causing disease progression. Currently, the relationship between in situ drug target availability (DTA), local protein expression and therapeutic response cannot be accurately assessed using existing analytical tools despite being crucial to understanding the mechanism of therapeutic efficacy. Procedure: We have previously reported development of our fluorescence imaging platform termed TRIPODD (Therapeutic Response Imaging through Proteomic and Optical Drug Distribution) that is capable of simultaneous quantification of single-cell DTA and protein expression with preserved spatial context within a tumor. TRIPODD combines two complementary fluorescence imaging techniques: intracellular paired agent imaging (iPAI) to measure DTA and cyclic immunofluorescence (cyCIF), which utilizes oligonucleotide conjugated antibodies (Ab-oligos) for spatial proteomic expression profiling on tissue samples. Herein, TRIPODD was modified and optimized to provide a downstream analysis of therapeutic response through single-cell DTA and proteomic response imaging. Results: We successfully performed sequential imaging of iPAI and cyCIF resulting in high dimensional imaging and biomarker assessment to quantify single-cell DTA and local protein expression on erlotinib treated NSCLC models. Pharmacodynamic and pharmacokinetic studies of the erlotinib iPAI probes revealed that administration of 2.5 mg/kg each of the targeted and untargeted probe 4 h prior to tumor collection enabled calculation of DTA values with high Pearson correlation to EGFR, the erlotinib molecular target, expression in the tumors. Analysis of single-cell biomarker expression revealed that a single erlotinib dose was insufficient to enact a measurable decrease in the EGFR signaling cascade protein expression, where only the DTA metric detected the presence of bound erlotinib. Conclusion: We demonstrated the capability of TRIPODD to evaluate therapeutic response imaging to erlotinib treatment as it relates to signaling inhibition, DTA, proliferation, and apoptosis with preserved spatial context.

Spatial Transcriptomics: A Powerful Tool in Disease Understanding and Drug Discovery.

Recent advancements in modern science have provided robust tools for drug discovery. The rapid development of transcriptome sequencing technologies has given rise to single-cell transcriptomics and single-nucleus transcriptomics, increasing the accuracy of sequencing and accelerating the drug discovery process. With the evolution of single-cell transcriptomics, spatial transcriptomics (ST) technology has emerged as a derivative approach. Spatial transcriptomics has emerged as a hot topic in the field of omics research in recent years; it not only provides information on gene expression levels but also offers spatial information on gene expression. This technology has shown tremendous potential in research on disease understanding and drug discovery. In this article, we introduce the analytical strategies of spatial transcriptomics and review its applications in novel target discovery and drug mechanism unravelling. Moreover, we discuss the current challenges and issues in this research field that need to be addressed. In conclusion, spatial transcriptomics offers a new perspective for drug discovery.

CD4+ T cell heterogeneity in gestational age and preeclampsia using single-cell RNA sequencing.

A balance between pro-inflammatory decidual CD4+ T cells and FOXP3+ regulatory T cells (FOXP3+ Tregs) is important for maintaining fetomaternal tolerance. Using single-cell RNA-sequencing and T cell receptor repertoire analysis, we determined that diversity and clonality of decidual CD4+ T cell subsets depend on gestational age. Th1/Th2 intermediate and Th1 subsets of CD4+ T cells were clonally expanded in both early and late gestation, whereas FOXP3+ Tregs were clonally expanded in late gestation. Th1/Th2 intermediate and FOXP3+ Treg subsets showed altered gene expression in preeclampsia (PE) compared to healthy late gestation. The Th1/Th2 intermediate subset exhibited elevated levels of cytotoxicity-related gene expression in PE. Moreover, increased Treg exhaustion was observed in the PE group, and FOXP3+ Treg subcluster analysis revealed that the effector Treg like subset drove the Treg exhaustion signatures in PE. The Th1/Th2 intermediate and effector Treg like subsets are possible inflammation-driving subsets in PE.

A comparative study on the dose-effect of low-dose radiation based on microdosimetric analysis and single-cell sequencing technology.

The biological mechanisms triggered by low-dose exposure still need to be explored in depth. In this study, the potential mechanisms of low-dose radiation when irradiating the BEAS-2B cell lines with a Cs-137 gamma-ray source were investigated through simulations and experiments. Monolayer cell population models were constructed for simulating and analyzing distributions of nucleus-specific energy within cell populations combined with the Monte Carlo method and microdosimetric analysis. Furthermore, the 10 × Genomics single-cell sequencing technology was employed to capture the heterogeneity of individual cell responses to low-dose radiation in the same irradiated sample. The numerical uncertainties can be found both in the specific energy distribution in microdosimetry and in differential gene expressions in radiation cytogenetics. Subsequently, the distribution of nucleus-specific energy was compared with the distribution of differential gene expressions to guide the selection of differential genes bioinformatics analysis. Dose inhomogeneity is pronounced at low doses, where an increase in dose corresponds to a decrease in the dispersion of cellular-specific energy distribution. Multiple screening of differential genes by microdosimetric features and statistical analysis indicate a number of potential pathways induced by low-dose exposure. It also provides a novel perspective on the selection of sensitive biomarkers that respond to low-dose radiation.

Label-free metabolic optical biomarkers track stem cell fate transition in real time.

Tracking stem cell fate transition is crucial for understanding their development and optimizing biomanufacturing. Destructive single-cell methods provide a pseudotemporal landscape of stem cell differentiation but cannot monitor stem cell fate in real time. We established a metabolic optical metric using label-free fluorescence lifetime imaging microscopy (FLIM), feature extraction and machine learning-assisted analysis, for real-time cell fate tracking. From a library of 205 metabolic optical biomarker (MOB) features, we identified 56 associated with hematopoietic stem cell (HSC) differentiation. These features collectively describe HSC fate transition and detect its bifurcate lineage choice. We further derived a MOB score measuring the "metabolic stemness" of single cells and distinguishing their division patterns. This score reveals a distinct role of asymmetric division in rescuing stem cells with compromised metabolic stemness and a unique mechanism of PI3K inhibition in promoting ex vivo HSC maintenance. MOB profiling is a powerful tool for tracking stem cell fate transition and improving their biomanufacturing from a single-cell perspective.

A new dawn for the study of cell type evolution.

Animal evolution is influenced by the emergence of new cell types, yet our understanding of this process remains elusive. This prompts the need for a broader exploration across diverse research organisms, facilitated by recent breakthroughs, such as gene editing tools and single-cell genomics. Essential to our understanding of cell type evolution is the accurate identification of homologous cells. We delve into the significance of considering developmental ontogeny and potential pitfalls when drawing conclusions about cell type homology. Additionally, we highlight recent discoveries in the study of cell type evolution through the application of single-cell transcriptomics and pinpoint areas ripe for further exploration.

Screening and validation of atherosclerosis PAN-apoptotic immune-related genes based on single-cell sequencing.

Background: Carotid atherosclerosis (CAS) is a complication of atherosclerosis (AS). PAN-optosome is an inflammatory programmed cell death pathway event regulated by the PAN-optosome complex. CAS's PAN-optosome-related genes (PORGs) have yet to be studied. Hence, screening the PAN-optosome-related diagnostic genes for treating CAS was vital. Methods: We introduced transcriptome data to screen out differentially expressed genes (DEGs) in CAS. Subsequently, WGCNA analysis was utilized to mine module genes about PANoptosis score. We performed differential expression analysis (CAS samples vs. standard samples) to obtain CAS-related differentially expressed genes at the single-cell level. Venn diagram was executed to identify PAN-optosome-related differential genes (POR-DEGs) associated with CAS. Further, LASSO regression and RF algorithm were implemented to were executed to build a diagnostic model. We additionally performed immune infiltration and gene set enrichment analysis (GSEA) based on diagnostic genes. We verified the accuracy of the model genes by single-cell nuclear sequencing and RT-qPCR validation of clinical samples, as well as in vitro cellular experiments. Results: We identified 785 DEGs associated with CAS. Then, 4296 module genes about PANoptosis score were obtained. We obtained the 7365 and 1631 CAS-related DEGs at the single-cell level, respectively. 67 POR-DEGs were retained Venn diagram. Subsequently, 4 PAN-optosome-related diagnostic genes (CNTN4, FILIP1, PHGDH, and TFPI2) were identified via machine learning. Cellular function tests on four genes showed that these genes have essential roles in maintaining arterial cell viability and resisting cellular senescence. Conclusion: We obtained four PANoptosis-related diagnostic genes (CNTN4, FILIP1, PHGDH, and TFPI2) associated with CAS, laying a theoretical foundation for treating CAS.

Diagnostic leukapheresis reveals distinct phenotypes of NSCLC circulating tumor cells.

BACKGROUND: Circulating tumor cells (CTCs) hold immense promise for unraveling tumor heterogeneity and understanding treatment resistance. However, conventional methods, especially in cancers like non-small cell lung cancer (NSCLC), often yield low CTC numbers, hindering comprehensive analyses. This study addresses this limitation by employing diagnostic leukapheresis (DLA) to cancer patients, enabling the screening of larger blood volumes. To leverage DLA's full potential, this study introduces a novel approach for CTC enrichment from DLAs. METHODS: DLA was applied to six advanced stage NSCLC patients. For an unbiased CTC enrichment, a two-step approach based on negative depletion of hematopoietic cells was used. Single-cell (sc) whole-transcriptome sequencing was performed, and CTCs were identified based on gene signatures and inferred copy number variations. RESULTS: Remarkably, this innovative approach led to the identification of unprecedented 3,363 CTC transcriptomes. The extensive heterogeneity among CTCs was unveiled, highlighting distinct phenotypes related to the epithelial-mesenchymal transition (EMT) axis, stemness, immune responsiveness, and metabolism. Comparison with sc transcriptomes from primary NSCLC cells revealed that CTCs encapsulate the heterogeneity of their primary counterparts while maintaining unique CTC-specific phenotypes. CONCLUSIONS: In conclusion, this study pioneers a transformative method for enriching CTCs from DLA, resulting in a substantial increase in CTC numbers. This allowed the creation of the first-ever single-cell whole transcriptome in-depth characterization of the heterogeneity of over 3,300 NSCLC-CTCs. The findings not only confirm the diagnostic value of CTCs in monitoring tumor heterogeneity but also propose a CTC-specific signature that can be exploited for targeted CTC-directed therapies in the future. This comprehensive approach signifies a major leap forward, positioning CTCs as a key player in advancing our understanding of cancer dynamics and paving the way for tailored therapeutic interventions.

Comparative single-cell multiplex immunophenotyping of therapy-naive patients with rheumatoid arthritis, systemic sclerosis, and systemic lupus erythematosus shed light on disease-specific composition of the peripheral immune system.

Introduction: Systemic autoimmune diseases (SADs) are a significant burden on the healthcare system. Understanding the complexity of the peripheral immunophenotype in SADs may facilitate the differential diagnosis and identification of potential therapeutic targets. Methods: Single-cell mass cytometric immunophenotyping was performed on peripheral blood mononuclear cells (PBMCs) from healthy controls (HCs) and therapy-naive patients with rheumatoid arthritis (RA), progressive systemic sclerosis (SSc), and systemic lupus erythematosus (SLE). Immunophenotyping was performed on 15,387,165 CD45+ live single cells from 52 participants (13 cases/group), using an antibody panel to detect 34 markers. Results: Using the t-SNE (t-distributed stochastic neighbor embedding) algorithm, the following 17 main immune cell types were determined: CD4+/CD57- T cells, CD4+/CD57+ T cells, CD8+/CD161- T cells, CD8+/CD161+/CD28+ T cells, CD8dim T cells, CD3+/CD4-/CD8- T cells, TCRγ/δ T cells, CD4+ NKT cells, CD8+ NKT cells, classic NK cells, CD56dim/CD98dim cells, B cells, plasmablasts, monocytes, CD11cdim/CD172dim cells, myeloid dendritic cells (mDCs), and plasmacytoid dendritic cells (pDCs). Seven of the 17 main cell types exhibited statistically significant frequencies in the investigated groups. The expression levels of the 34 markers in the main populations were compared between HCs and SADs. In summary, 59 scatter plots showed significant differences in the expression intensities between at least two groups. Next, each immune cell population was divided into subpopulations (metaclusters) using the FlowSOM (self-organizing map) algorithm. Finally, 121 metaclusters (MCs) of the 10 main immune cell populations were found to have significant differences to classify diseases. The single-cell T-cell heterogeneity represented 64MCs based on the expression of 34 markers, and the frequency of 23 MCs differed significantly between at least twoconditions. The CD3- non-T-cell compartment contained 57 MCs with 17 MCs differentiating at least two investigated groups. In summary, we are the first to demonstrate the complexity of the immunophenotype of 34 markers over 15 million single cells in HCs vs. therapy-naive patients with RA, SSc, and SLE. Disease specific population frequencies or expression patterns of peripheral immune cells provide a single-cell data resource to the scientific community.

The future of rapid and automated single-cell data analysis using reference mapping.

As the number of single-cell datasets continues to grow rapidly, workflows that map new data to well-curated reference atlases offer enormous promise for the biological community. In this perspective, we discuss key computational challenges and opportunities for single-cell reference-mapping algorithms. We discuss how mapping algorithms will enable the integration of diverse datasets across disease states, molecular modalities, genetic perturbations, and diverse species and will eventually replace manual and laborious unsupervised clustering pipelines.

Data normalization for addressing the challenges in the analysis of single-cell transcriptomic datasets.

BACKGROUND: Normalization is a critical step in the analysis of single-cell RNA-sequencing (scRNA-seq) datasets. Its main goal is to make gene counts comparable within and between cells. To do so, normalization methods must account for technical and biological variability. Numerous normalization methods have been developed addressing different sources of dispersion and making specific assumptions about the count data. MAIN BODY: The selection of a normalization method has a direct impact on downstream analysis, for example differential gene expression and cluster identification. Thus, the objective of this review is to guide the reader in making an informed decision on the most appropriate normalization method to use. To this aim, we first give an overview of the different single cell sequencing platforms and methods commonly used including isolation and library preparation protocols. Next, we discuss the inherent sources of variability of scRNA-seq datasets. We describe the categories of normalization methods and include examples of each. We also delineate imputation and batch-effect correction methods. Furthermore, we describe data-driven metrics commonly used to evaluate the performance of normalization methods. We also discuss common scRNA-seq methods and toolkits used for integrated data analysis. CONCLUSIONS: According to the correction performed, normalization methods can be broadly classified as within and between-sample algorithms. Moreover, with respect to the mathematical model used, normalization methods can further be classified into: global scaling methods, generalized linear models, mixed methods, and machine learning-based methods. Each of these methods depict pros and cons and make different statistical assumptions. However, there is no better performing normalization method. Instead, metrics such as silhouette width, K-nearest neighbor batch-effect test, or Highly Variable Genes are recommended to assess the performance of normalization methods.

Integrating machine learning algorithms and single-cell analysis to identify gut microbiota-related macrophage biomarkers in atherosclerotic plaques.

Objective: The relationship between macrophages and the gut microbiota in patients with atherosclerosis remains poorly defined, and effective biological markers are lacking. This study aims to elucidate the interplay between gut microbial communities and macrophages, and to identify biomarkers associated with the destabilization of atherosclerotic plaques. The goal is to enhance our understanding of the underlying molecular pathways and to pave new avenues for diagnostic approaches and therapeutic strategies in the disease. Methods: This study employed Weighted Gene Co-expression Network Analysis (WGCNA) and differential expression analysis on atherosclerosis datasets to identify macrophage-associated genes and quantify the correlation between these genes and gut microbiota gene sets. The Random Forest algorithm was utilized to pinpoint PLEK, IRF8, BTK, CCR1, and CD68 as gut microbiota-related macrophage genes, and a nomogram was constructed. Based on the top five genes, a Non-negative Matrix Factorization (NMF) algorithm was applied to construct gut microbiota-related macrophage clusters and analyze their potential biological alterations. Subsequent single-cell analyses were conducted to observe the expression patterns of the top five genes and the interactions between immune cells. Finally, the expression profiles of key molecules were validated using clinical samples from atherosclerosis patients. Results: Utilizing the Random Forest algorithm, we ultimately identified PLEK, IRF8, CD68, CCR1, and BTK as gut microbiota-associated macrophage genes that are upregulated in atherosclerotic plaques. A nomogram based on the expression of these five genes was constructed for use as an auxiliary tool in clinical diagnosis. Single-cell analysis confirmed the specific expression of gut microbiota-associated macrophage genes in macrophages. Clinical samples substantiated the high expression of PLEK in unstable atherosclerotic plaques. Conclusion: Gut microbiota-associated macrophage genes (PLEK, IRF8, CD68, CCR1, and BTK) may be implicated in the pathogenesis of atherosclerotic plaques and could serve as diagnostic markers to aid patients with atherosclerosis.

Single-cell transcriptome analysis reveals heterogeneity of neutrophils in non-small cell lung cancer.

BACKGROUND: Lung cancer stands out as a highly perilous malignant tumor with severe implications for human health. There has been a growing interest in neutrophils as a result of their role in promoting cancer in recent years. Thus, the present study aimed to investigate the heterogeneity of neutrophils in non-small cell lung cancer (NSCLC). METHODS: Single-cell RNA sequencing of tumor-associated neutrophils (TANs) and polymorphonuclear neutrophils sourced from the Gene Expression Omnibus database was analyzed. Moreover, cell-cell communication, differentiation trajectories and transcription factor analyses were performed. RESULTS: Neutrophils were found to be closely associated with macrophages. Four major types of TANs were identified: a transitional subcluster that migrated from blood to tumor microenvironment (TAN-0), an inflammatory subcluster (TAN-1), a subpopulation that displayed a distinctive transcriptional signature (TAN-2) and a final differentiation state that promoted tumor formation (TAN-3). Meanwhile, TAN-3 displayed a marked increase in glycolytic activity. Finally, transcription factors were analyzed to uncover distinct TAN cluster-specific regulons. CONCLUSIONS: The discovery of the dynamic characteristics of TANs in the present study is anticipated to contribute to yielding a better understanding of the tumor microenvironment and advancing the treatment of NSCLC.

EmptyDropsMultiome discriminates real cells from background in single-cell multiomics assays.

Multiomic droplet-based technologies allow different molecular modalities, such as chromatin accessibility and gene expression (scATAC-seq and scRNA-seq), to be probed in the same nucleus. We develop EmptyDropsMultiome, an approach that distinguishes true nuclei-containing droplets from background. Using simulations, we show that EmptyDropsMultiome has higher statistical power and accuracy than existing approaches, including CellRanger-arc and EmptyDrops. On real datasets, we observe that CellRanger-arc misses more than half of the nuclei identified by EmptyDropsMultiome and, moreover, is biased against certain cell types, some of which have a retrieval rate lower than 20%.

SpatialCells: automated profiling of tumor microenvironments with spatially resolved multiplexed single-cell data.

Cancer is a complex cellular ecosystem where malignant cells coexist and interact with immune, stromal and other cells within the tumor microenvironment (TME). Recent technological advancements in spatially resolved multiplexed imaging at single-cell resolution have led to the generation of large-scale and high-dimensional datasets from biological specimens. This underscores the necessity for automated methodologies that can effectively characterize molecular, cellular and spatial properties of TMEs for various malignancies. This study introduces SpatialCells, an open-source software package designed for region-based exploratory analysis and comprehensive characterization of TMEs using multiplexed single-cell data. The source code and tutorials are available at https://semenovlab.github.io/SpatialCells. SpatialCells efficiently streamlines the automated extraction of features from multiplexed single-cell data and can process samples containing millions of cells. Thus, SpatialCells facilitates subsequent association analyses and machine learning predictions, making it an essential tool in advancing our understanding of tumor growth, invasion and metastasis.

AnnoSpat annotates cell types and quantifies cellular arrangements from spatial proteomics.

Cellular composition and anatomical organization influence normal and aberrant organ functions. Emerging spatial single-cell proteomic assays such as Image Mass Cytometry (IMC) and Co-Detection by Indexing (CODEX) have facilitated the study of cellular composition and organization by enabling high-throughput measurement of cells and their localization directly in intact tissues. However, annotation of cell types and quantification of their relative localization in tissues remain challenging. To address these unmet needs for atlas-scale datasets like Human Pancreas Analysis Program (HPAP), we develop AnnoSpat (Annotator and Spatial Pattern Finder) that uses neural network and point process algorithms to automatically identify cell types and quantify cell-cell proximity relationships. Our study of data from IMC and CODEX shows the higher performance of AnnoSpat in rapid and accurate annotation of cell types compared to alternative approaches. Moreover, the application of AnnoSpat to type 1 diabetic, non-diabetic autoantibody-positive, and non-diabetic organ donor cohorts recapitulates known islet pathobiology and shows differential dynamics of pancreatic polypeptide (PP) cell abundance and CD8+ T cells infiltration in islets during type 1 diabetes progression.

Integrated analysis of single cell-RNA sequencing and Mendelian randomization identifies lactate dehydrogenase B as a target of melatonin in ischemic stroke.

AIMS: Despite the success of single-cell RNA sequencing in identifying cellular heterogeneity in ischemic stroke, clarifying the mechanisms underlying these associations of differently expressed genes remains challenging. Several studies that integrate gene expression and gene expression quantitative trait loci (eQTLs) with genome wide-association study (GWAS) data to determine their causal role have been proposed. METHODS: Here, we combined Mendelian randomization (MR) framework and single cell (sc) RNA sequencing to study how differently expressed genes (DEGs) mediating the effect of gene expression on ischemic stroke. The hub gene was further validated in the in vitro model. RESULTS: We identified 2339 DEGs in 10 cell clusters. Among these DEGs, 58 genes were associated with the risk of ischemic stroke. After external validation with eQTL dataset, lactate dehydrogenase B (LDHB) is identified to be positively associated with ischemic stroke. The expression of LDHB has also been validated in sc RNA-seq with dominant expression in microglia and astrocytes, and melatonin is able to reduce the LDHB expression and activity in vitro ischemic models. CONCLUSION: Our study identifies LDHB as a novel biomarker for ischemic stroke via combining the sc RNA-seq and MR analysis.

Single-cell sequencing of tumour infiltrating T cells efficiently identifies tumour-specific T cell receptors based on the T cell activation score.

Adoptively transferred T cell receptor-engineered T cells are a promising cancer treatment strategy, and the identification of tumour-specific TCRs is essential. Previous studies reported that tumour-reactive T cells and TCRs could be isolated based on the expression of activation markers. However, since T cells with different cell states could not respond uniformly to activation but show a heterogeneous expression profile of activation and effector molecules, isolation of tumour-reactive T cells based on single activation or effector molecules could result in the absence of tumour-reactive T cells; thus, combinations of multiple activation and effector molecules could improve the efficiency of isolating tumour-specific TCRs. We enrolled two patients with lung adenocarcinoma and obtained their tumour infiltrating lymphocytes (TILs) and autologous tumour cells (ATCs). TILs were cocultured with the corresponding ATCs for 12 h and subjected to single-cell RNA sequencing. First, we identified three TCRs with the highest expression levels of IFNG and TNFRSF9 mRNA for each patient, yet only the top one or two recognized the corresponding ATCs in each patient. Next, we defined the activation score based on normalized expression levels of IFNG, IL2, TNF, IL2RA, CD69, TNFRSF9, GZMB, GZMA, GZMK, and PRF1 mRNA for each T cell and then identified three TCRs with the highest activation score for each patient. We found that all three TCRs in each patient could specifically identify corresponding ATCs. In conclusion, we established an efficient approach to isolate tumour-reactive TCRs based on combinations of multiple activation and effector molecules through single-cell RNA sequencing.

Multiplexed bulk and single-cell RNA-seq hybrid enables cost-efficient disease modeling with chimeric organoids.

Disease modeling with isogenic Induced Pluripotent Stem Cell (iPSC)-differentiated organoids serves as a powerful technique for studying disease mechanisms. Multiplexed coculture is crucial to mitigate batch effects when studying the genetic effects of disease-causing variants in differentiated iPSCs or organoids, and demultiplexing at the single-cell level can be conveniently achieved by assessing natural genetic barcodes. Here, to enable cost-efficient time-series experimental designs via multiplexed bulk and single-cell RNA-seq of hybrids, we introduce a computational method in our Vireo Suite, Vireo-bulk, to effectively deconvolve pooled bulk RNA-seq data by genotype reference, and thereby quantify donor abundance over the course of differentiation and identify differentially expressed genes among donors. Furthermore, with multiplexed scRNA-seq and bulk RNA-seq, we demonstrate the usefulness and necessity of a pooled design to reveal donor iPSC line heterogeneity during macrophage cell differentiation and to model rare WT1 mutation-driven kidney disease with chimeric organoids. Our work provides an experimental and analytic pipeline for dissecting disease mechanisms with chimeric organoids.