Kernel-based testing for single-cell differential analysis.

Single-cell technologies offer insights into molecular feature distributions, but comparing them poses challenges. We propose a kernel-testing framework for non-linear cell-wise distribution comparison, analyzing gene expression and epigenomic modifications. Our method allows feature-wise and global transcriptome/epigenome comparisons, revealing cell population heterogeneities. Using a classifier based on embedding variability, we identify transitions in cell states, overcoming limitations of traditional single-cell analysis. Applied to single-cell ChIP-Seq data, our approach identifies untreated breast cancer cells with an epigenomic profile resembling persister cells. This demonstrates the effectiveness of kernel testing in uncovering subtle population variations that might be missed by other methods.

Machine learning and deep learning to identifying subarachnoid haemorrhage macrophage-associated biomarkers by bulk and single-cell sequencing.

We investigated subarachnoid haemorrhage (SAH) macrophage subpopulations and identified relevant key genes for improving diagnostic and therapeutic strategies. SAH rat models were established, and brain tissue samples underwent single-cell transcriptome sequencing and bulk RNA-seq. Using single-cell data, distinct macrophage subpopulations, including a unique SAH subset, were identified. The hdWGCNA method revealed 160 key macrophage-related genes. Univariate analysis and lasso regression selected 10 genes for constructing a diagnostic model. Machine learning algorithms facilitated model development. Cellular infiltration was assessed using the MCPcounter algorithm, and a heatmap integrated cell abundance and gene expression. A 3 × 3 convolutional neural network created an additional diagnostic model, while molecular docking identified potential drugs. The diagnostic model based on the 10 selected genes achieved excellent performance, with an AUC of 1 in both training and validation datasets. The heatmap, combining cell abundance and gene expression, provided insights into SAH cellular composition. The convolutional neural network model exhibited a sensitivity and specificity of 1 in both datasets. Additionally, CD14, GPNMB, SPP1 and PRDX5 were specifically expressed in SAH-associated macrophages, highlighting its potential as a therapeutic target. Network pharmacology analysis identified some targeting drugs for SAH treatment. Our study characterised SAH macrophage subpopulations and identified key associated genes. We developed a robust diagnostic model and recognised CD14, GPNMB, SPP1 and PRDX5 as potential therapeutic targets. Further experiments and clinical investigations are needed to validate these findings and explore the clinical implications of targets in SAH treatment.

SC-Track: a robust cell-tracking algorithm for generating accurate single-cell lineages from diverse cell segmentations.

Computational analysis of fluorescent timelapse microscopy images at the single-cell level is a powerful approach to study cellular changes that dictate important cell fate decisions. Core to this approach is the need to generate reliable cell segmentations and classifications necessary for accurate quantitative analysis. Deep learning-based convolutional neural networks (CNNs) have emerged as a promising solution to these challenges. However, current CNNs are prone to produce noisy cell segmentations and classifications, which is a significant barrier to constructing accurate single-cell lineages. To address this, we developed a novel algorithm called Single Cell Track (SC-Track), which employs a hierarchical probabilistic cache cascade model based on biological observations of cell division and movement dynamics. Our results show that SC-Track performs better than a panel of publicly available cell trackers on a diverse set of cell segmentation types. This cell-tracking performance was achieved without any parameter adjustments, making SC-Track an excellent generalized algorithm that can maintain robust cell-tracking performance in varying cell segmentation qualities, cell morphological appearances and imaging conditions. Furthermore, SC-Track is equipped with a cell class correction function to improve the accuracy of cell classifications in multiclass cell segmentation time series. These features together make SC-Track a robust cell-tracking algorithm that works well with noisy cell instance segmentation and classification predictions from CNNs to generate accurate single-cell lineages and classifications.

Unravelling B cell heterogeneity: insights into flow cytometry-gated B cells from single-cell multi-omics data.

Introduction: B cells play a pivotal role in adaptive immunity which has been extensively characterised primarily via flow cytometry-based gating strategies. This study addresses the discrepancies between flow cytometry-defined B cell subsets and their high-confidence molecular signatures using single-cell multi-omics approaches. Methods: By analysing multi-omics single-cell data from healthy individuals and patients across diseases, we characterised the level and nature of cellular contamination within standard flow cytometric-based gating, resolved some of the ambiguities in the literature surrounding unconventional B cell subsets, and demonstrated the variable effects of flow cytometric-based gating cellular heterogeneity across diseases. Results: We showed that flow cytometric-defined B cell populations are heterogenous, and the composition varies significantly between disease states thus affecting the implications of functional studies performed on these populations. Importantly, this paper draws caution on findings about B cell selection and function of flow cytometric-sorted populations, and their roles in disease. As a solution, we developed a simple tool to identify additional markers that can be used to increase the purity of flow-cytometric gated immune cell populations based on multi-omics data (AlliGateR). Here, we demonstrate that additional non-linear CD20, CD21 and CD24 gating can increase the purity of both naïve and memory populations. Discussion: These findings underscore the need to reconsider B cell subset definitions within the literature and propose leveraging single-cell multi-omics data for refined characterisation. We show that single-cell multi-omics technologies represent a powerful tool to bridge the gap between surface marker-based annotations and the intricate molecular characteristics of B cell subsets.

Learning Consistency and Specificity of Cells From Single-Cell Multi-Omic Data.

Advancements in single-cell technologies concomitantly develop the epigenomic and transcriptomic profiles at the cell levels, providing opportunities to explore the potential biological mechanisms. Even though significant efforts have been dedicated to them, it remains challenging for the integration analysis of multi-omic data of single-cell because of the heterogeneity, complicated coupling and interpretability of data. To handle these issues, we propose a novel self-representation Learning-based Multi-omics data Integrative Clustering algorithm (sLMIC) for the integration of single-cell epigenomic profiles (DNA methylation or scATAC-seq) and transcriptomic (scRNA-seq), which the consistent and specific features of cells are explicitly extracted facilitating the cell clustering. Specifically, sLMIC constructs a graph for each type of single-cell data, thereby transforming omics data into multi-layer networks, which effectively removes heterogeneity of omic data. Then, sLMIC employs the low-rank and exclusivity constraints to separate the self-representation of cells into two parts, i.e., the shared and specific features, which explicitly characterize the consistency and diversity of omic data, providing an effective strategy to model the structure of cell types. Feature extraction and cell clustering are jointly formulated as an overall objective function, where latent features of data are obtained under the guidance of cell clustering. The extensive experimental results on 13 multi-omics datasets of single-cell from diverse organisms and tissues indicate that sLMIC observably exceeds the advanced algorithms regarding various measurements.

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.

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.

Decoding the principle of cell-fate determination for its reverse control.

Understanding and manipulating cell fate determination is pivotal in biology. Cell fate is determined by intricate and nonlinear interactions among molecules, making mathematical model-based quantitative analysis indispensable for its elucidation. Nevertheless, obtaining the essential dynamic experimental data for model development has been a significant obstacle. However, recent advancements in large-scale omics data technology are providing the necessary foundation for developing such models. Based on accumulated experimental evidence, we can postulate that cell fate is governed by a limited number of core regulatory circuits. Following this concept, we present a conceptual control framework that leverages single-cell RNA-seq data for dynamic molecular regulatory network modeling, aiming to identify and manipulate core regulatory circuits and their master regulators to drive desired cellular state transitions. We illustrate the proposed framework by applying it to the reversion of lung cancer cell states, although it is more broadly applicable to understanding and controlling a wide range of cell-fate determination processes.

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.

Integration of single-cell sequencing and bulk RNA-seq to identify and develop a prognostic signature related to colorectal cancer stem cells.

The prognosis for patients with colorectal cancer (CRC) remains worse than expected due to metastasis, recurrence, and resistance to chemotherapy. Colorectal cancer stem cells (CRCSCs) play a vital role in tumor metastasis, recurrence, and chemotherapy resistance. However, there are currently no prognostic markers based on CRCSCs-related genes available for clinical use. In this study, single-cell transcriptome sequencing was employed to distinguish cancer stem cells (CSCs) in the CRC microenvironment and analyze their properties at the single-cell level. Subsequently, data from TCGA and GEO databases were utilized to develop a prognostic risk model for CRCSCs-related genes and validate its diagnostic performance. Additionally, functional enrichment, immune response, and chemotherapeutic drug sensitivity of the relevant genes in the risk model were investigated. Lastly, the key gene RPS17 in the risk model was identified as a potential prognostic marker and therapeutic target for further comprehensive studies. Our findings provide new insights into the prognostic treatment of CRC and offer novel perspectives for a systematic and comprehensive understanding of CRC development.

Molecular mechanism of radiation tolerance in lung adenocarcinoma cells using single-cell RNA sequencing.

The efficacy of radiotherapy, a cornerstone in the treatment of lung adenocarcinoma (LUAD), is profoundly undermined by radiotolerance. This resistance not only poses a significant clinical challenge but also compromises patient survival rates. Therefore, it is important to explore this mechanism for the treatment of LUAD. Multiple public databases were used for single-cell RNA sequencing (scRNA-seq) data. We filtered, normalized and downscaled scRNA-seq data based on the Seurat package to obtain different cell subpopulations. Subsequently, the ssGSEA algorithm was used to assess the enrichment scores of the different cell subpopulations, and thus screen the cell subpopulations that are most relevant to radiotherapy tolerance based on the Pearson method. Finally, pseudotime analysis was performed, and a preliminary exploration of gene mutations in different cell subpopulations was performed. We identified HIST1H1D+ A549 and PIF1+ A549 as the cell subpopulations related to radiotolerance. The expression levels of cell cycle-related genes and pathway enrichment scores of these two cell subpopulations increased gradually with the extension of radiation treatment time. Finally, we found that the proportion of TP53 mutations in patients who had received radiotherapy was significantly higher than that in patients who had not received radiotherapy. We identified two cellular subpopulations associated with radiotherapy tolerance, which may shed light on the molecular mechanisms of radiotherapy tolerance in LUAD and provide new clinical perspectives.

Isolation of Nuclei from Human Intermuscular Adipose Tissue and Downstream Single-Nuclei RNA Sequencing.

Intermuscular adipose tissue (IMAT) is a relatively understudied adipose depot located between muscle fibers. IMAT content increases with age and BMI and is associated with metabolic and muscle degenerative diseases; however, an understanding of the biological properties of IMAT and its interplay with the surrounding muscle fibers is severely lacking. In recent years, single-cell and nuclei RNA sequencing have provided us with cell type-specific atlases of several human tissues. However, the cellular composition of human IMAT remains largely unexplored due to the inherent challenges of its accessibility from biopsy collection in humans. In addition to the limited amount of tissue collected, the processing of human IMAT is complicated due to its proximity to skeletal muscle tissue and fascia. The lipid-laden nature of the adipocytes makes it incompatible with single-cell isolation. Hence, single nuclei RNA sequencing is optimal for obtaining high-dimensional transcriptomics at single-cell resolution and provides the potential to uncover the biology of this depot, including the exact cellular composition of IMAT. Here, we present a detailed protocol for nuclei isolation and library preparation of frozen human IMAT for single nuclei RNA sequencing. This protocol allows for the profiling of thousands of nuclei using a droplet-based approach, thus providing the capacity to detect rare and low-abundant cell types.

Error modelled gene expression analysis (EMOGEA) provides a superior overview of time course RNA-seq measurements and low count gene expression.

Temporal RNA-sequencing (RNA-seq) studies of bulk samples provide an opportunity for improved understanding of gene regulation during dynamic phenomena such as development, tumor progression or response to an incremental dose of a pharmacotherapeutic. Moreover, single-cell RNA-seq (scRNA-seq) data implicitly exhibit temporal characteristics because gene expression values recapitulate dynamic processes such as cellular transitions. Unfortunately, temporal RNA-seq data continue to be analyzed by methods that ignore this ordinal structure and yield results that are often difficult to interpret. Here, we present Error Modelled Gene Expression Analysis (EMOGEA), a framework for analyzing RNA-seq data that incorporates measurement uncertainty, while introducing a special formulation for those acquired to monitor dynamic phenomena. This method is specifically suited for RNA-seq studies in which low-count transcripts with small-fold changes lead to significant biological effects. Such transcripts include genes involved in signaling and non-coding RNAs that inherently exhibit low levels of expression. Using simulation studies, we show that this framework down-weights samples that exhibit extreme responses such as batch effects allowing them to be modeled with the rest of the samples and maintain the degrees of freedom originally envisioned for a study. Using temporal experimental data, we demonstrate the framework by extracting a cascade of gene expression waves from a well-designed RNA-seq study of zebrafish embryogenesis and an scRNA-seq study of mouse pre-implantation and provide unique biological insights into the regulation of genes in each wave. For non-ordinal measurements, we show that EMOGEA has a much higher rate of true positive calls and a vanishingly small rate of false negative discoveries compared to common approaches. Finally, we provide two packages in Python and R that are self-contained and easy to use, including test data.

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.

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.

Multiregional transcriptomics identifies congruent consensus subtypes with prognostic value beyond tumor heterogeneity of colorectal cancer.

Intra-tumor heterogeneity compromises the clinical value of transcriptomic classifications of colorectal cancer. We investigated the prognostic effect of transcriptomic heterogeneity and the potential for classifications less vulnerable to heterogeneity in a single-hospital series of 1093 tumor samples from 692 patients, including multiregional samples from 98 primary tumors and 35 primary-metastasis sets. We show that intra-tumor heterogeneity of the consensus molecular subtypes (CMS) is frequent and has poor-prognostic associations independently of tumor microenvironment markers. Multiregional transcriptomics uncover cancer cell-intrinsic and low-heterogeneity signals that recapitulate the intrinsic CMSs proposed by single-cell sequencing. Further subclassification identifies congruent CMSs that explain a larger proportion of variation in patient survival than intra-tumor heterogeneity. Plasticity is indicated by discordant intrinsic phenotypes of matched primary and metastatic tumors. We conclude that multiregional sampling reconciles the prognostic power of tumor classifications from single-cell and bulk transcriptomics in the context of intra-tumor heterogeneity, and phenotypic plasticity challenges the reconciliation of primary and metastatic subtypes.

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.

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.

Integrating single-cell and spatial analysis reveals MUC1-mediated cellular crosstalk in mucinous colorectal adenocarcinoma.

BACKGROUND: Mucinous colorectal adenocarcinoma (MCA) is a distinct subtype of colorectal cancer (CRC) with the most aggressive pattern, but effective treatment of MCA remains a challenge due to its vague pathological characteristics. An in-depth understanding of transcriptional dynamics at the cellular level is critical for developing specialised MCA treatment strategies. METHODS: We integrated single-cell RNA sequencing and spatial transcriptomics data to systematically profile the MCA tumor microenvironment (TME), particularly the interactome of stromal and immune cells. In addition, a three-dimensional bioprinting technique, canonical ex vivo co-culture system, and immunofluorescence staining were further applied to validate the cellular communication networks within the TME. RESULTS: This study identified the crucial intercellular interactions that engaged in MCA pathogenesis. We found the increased infiltration of FGF7+/THBS1+ myofibroblasts in MCA tissues with decreased expression of genes associated with leukocyte-mediated immunity and T cell activation, suggesting a crucial role of these cells in regulating the immunosuppressive TME. In addition, MS4A4A+ macrophages that exhibit M2-phenotype were enriched in the tumoral niche and high expression of MS4A4A+ was associated with poor prognosis in the cohort data. The ligand-receptor-based intercellular communication analysis revealed the tight interaction of MUC1+ malignant cells and ZEB1+ endothelial cells, providing mechanistic information for MCA angiogenesis and molecular targets for subsequent translational applications. CONCLUSIONS: Our study provides novel insights into communications among tumour cells with stromal and immune cells that are significantly enriched in the TME during MCA progression, presenting potential prognostic biomarkers and therapeutic strategies for MCA. KEY POINTS: Tumour microenvironment profiling of MCA is developed. MUC1+ tumour cells interplay with FGF7+/THBS1+ myofibroblasts to promote MCA development. MS4A4A+ macrophages exhibit M2 phenotype in MCA. ZEB1+ endotheliocytes engage in EndMT process in MCA.

Spatiotemporal dynamics of the proton motive force on single bacterial cells.

Electrochemical gradients across biological membranes are vital for cellular bioenergetics. In bacteria, the proton motive force (PMF) drives essential processes like adenosine triphosphate production and motility. Traditionally viewed as temporally and spatially stable, recent research reveals a dynamic PMF behavior at both single-cell and community levels. Moreover, the observed lateral segregation of respiratory complexes could suggest a spatial heterogeneity of the PMF. Using a light-activated proton pump and detecting the activity of the bacterial flagellar motor, we perturb and probe the PMF of single cells. Spatially homogeneous PMF perturbations reveal millisecond-scale temporal dynamics and an asymmetrical capacitive response. Localized perturbations show a rapid lateral PMF homogenization, faster than proton diffusion, akin to the electrotonic potential spread observed in passive neurons, explained by cable theory. These observations imply a global coupling between PMF sources and consumers along the membrane, precluding sustained PMF spatial heterogeneity but allowing for rapid temporal changes.