Advanced methods for gene network identification and noise decomposition from single-cell data.

Central to analyzing noisy gene expression systems is solving the Chemical Master Equation (CME), which characterizes the probability evolution of the reacting species' copy numbers. Solving CMEs for high-dimensional systems suffers from the curse of dimensionality. Here, we propose a computational method for improved scalability through a divide-and-conquer strategy that optimally decomposes the whole system into a leader system and several conditionally independent follower subsystems. The CME is solved by combining Monte Carlo estimation for the leader system with stochastic filtering procedures for the follower subsystems. We demonstrate this method with high-dimensional numerical examples and apply it to identify a yeast transcription system at the single-cell resolution, leveraging mRNA time-course experimental data. The identification results enable an accurate examination of the heterogeneity in rate parameters among isogenic cells. To validate this result, we develop a noise decomposition technique exploiting time-course data but requiring no supplementary components, e.g., dual-reporters.

Unravelling the metabolic landscape of cutaneous melanoma: Insights from single-cell sequencing analysis and machine learning for prognostic assessment of lactate metabolism.

This manuscript presents a comprehensive investigation into the role of lactate metabolism-related genes as potential prognostic markers in skin cutaneous melanoma (SKCM). Bulk-transcriptome data from The Cancer Genome Atlas (TCGA) and GSE19234, GSE22153, and GSE65904 cohorts from GEO database were processed and harmonized to mitigate batch effects. Lactate metabolism scores were assigned to individual cells using the 'AUCell' package. Weighted Co-expression Network Analysis (WGCNA) was employed to identify gene modules correlated with lactate metabolism. Machine learning algorithms were applied to construct a prognostic model, and its performance was evaluated in multiple cohorts. Immune correlation, mutation analysis, and enrichment analysis were conducted to further characterize the prognostic model's biological implications. Finally, the function of key gene NDUFS7 was verified by cell experiments. Machine learning resulted in an optimal prognostic model, demonstrating significant prognostic value across various cohorts. In the different cohorts, the high-risk group showed a poor prognosis. Immune analysis indicated differences in immune cell infiltration and checkpoint gene expression between risk groups. Mutation analysis identified genes with high mutation loads in SKCM. Enrichment analysis unveiled enriched pathways and biological processes in high-risk SKCM patients. NDUFS7 was found to be a hub gene in the protein-protein interaction network. After the expression of NDUFS7 was reduced by siRNA knockdown, CCK-8, colony formation, transwell and wound healing tests showed that the activity, proliferation and migration of A375 and WM115 cell lines were significantly decreased. This study offers insights into the prognostic significance of lactate metabolism-related genes in SKCM.

Optimal sequencing budget allocation for trajectory reconstruction of single cells.

BACKGROUND: Charting cellular trajectories over gene expression is key to understanding dynamic cellular processes and their underlying mechanisms. While advances in single-cell RNA-sequencing technologies and computational methods have pushed forward the recovery of such trajectories, trajectory inference remains a challenge due to the noisy, sparse, and high-dimensional nature of single-cell data. This challenge can be alleviated by increasing either the number of cells sampled along the trajectory (breadth) or the sequencing depth, i.e. the number of reads captured per cell (depth). Generally, these two factors are coupled due to an inherent breadth-depth tradeoff that arises when the sequencing budget is constrained due to financial or technical limitations. RESULTS: Here we study the optimal allocation of a fixed sequencing budget to optimize the recovery of trajectory attributes. Empirical results reveal that reconstruction accuracy of internal cell structure in expression space scales with the logarithm of either the breadth or depth of sequencing. We additionally observe a power law relationship between the optimal number of sampled cells and the corresponding sequencing budget. For linear trajectories, non-monotonicity in trajectory reconstruction across the breadth-depth tradeoff can impact downstream inference, such as expression pattern analysis along the trajectory. We demonstrate these results for five single-cell RNA-sequencing datasets encompassing differentiation of embryonic stem cells, pancreatic beta cells, hepatoblast and multipotent hematopoietic cells, as well as induced reprogramming of embryonic fibroblasts into neurons. By addressing the challenges of single-cell data, our study offers insights into maximizing the efficiency of cellular trajectory analysis through strategic allocation of sequencing resources.

Assessment of biological organ age using molecular pathology in pre-transplant kidney biopsies.

Organ shortage is a major challenge in kidney transplantation but the use of older donors, often with co-morbidities, is hampered by inconsistent outcomes. Methods of accurately stratifying marginal donor organs by clinical and histological assessment are lacking. To better understand organ variability, we profiled the transcriptomes of 271 kidneys from deceased donors at retrieval. Following correction for biopsy composition, we assessed molecular pathways that associated with delayed, and sub-optimal one-year graft function. Analysis of cortical biopsies identified an adaptive immune gene-rich module that significantly associated with increasing age and worse outcomes. Cellular deconvolution using human kidney reference single cell transcriptomes confirmed an increase in kidney-specific B and T cell signatures, as well as kidney macrophage, myofibroblast and fibroblast gene sets in this module. Surprisingly, innate immune pathway and neutrophil gene signature enrichment was associated with better outcomes. Thus, our work uncovers cellular molecular features of pathological organ ageing, identifiable at kidney retrieval, with translational potential.

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.

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.

Hidden Markov random field models for cell-type assignment of spatially resolved transcriptomics.

MOTIVATION: The recent development of spatially resolved transcriptomics (SRT) technologies has facilitated research on gene expression in the spatial context. Annotating cell types is one crucial step for downstream analysis. However, many existing algorithms use an unsupervised strategy to assign cell types for SRT data. They first conduct clustering analysis and then aggregate cluster-level expression based on the clustering results. This workflow fails to leverage the marker gene information efficiently. On the other hand, other cell annotation methods designed for single-cell RNA-seq data utilize the cell-type marker genes information but fail to use spatial information in SRT data. RESULTS: We introduce a statistical spatial transcriptomics cell assignment model, SPAN, to annotate clusters of cells or spots into known types in SRT data with prior knowledge of predefined marker genes and spatial information. The SPAN model annotates cells or spots from SRT data using predefined overexpressed marker genes and combines a mixture model with a hidden Markov random field to model the spatial dependency between neighboring spots. We demonstrate the effectiveness of SPAN against spatial and nonspatial clustering algorithms through extensive simulation and real data experiments. AVAILABILITY AND IMPLEMENTATION: https://github.com/ChengZ352/SPAN.

Profiling Transcriptional Heterogeneity with Seq-Well S3: A Low-Cost, Portable, High-Fidelity Platform for Massively Parallel Single-Cell RNA-Seq.

Seq-Well is a high-throughput, picowell-based single-cell RNA-seq technology that can be used to simultaneously profile the transcriptomes of thousands of cells (Gierahn et al. Nat Methods 14(4):395-398, 2017). Relative to its reverse-emulsion-droplet-based counterparts, Seq-Well addresses key cost, portability, and scalability limitations. Recently, we introduced an improved molecular biology for Seq-Well to enhance the information content that can be captured from individual cells using the platform. This update, which we call Seq-Well S3 (S3: Second-Strand Synthesis), incorporates a second-strand-synthesis step after reverse transcription to boost the detection of cellular transcripts normally missed when running the original Seq-Well protocol (Hughes et al. Immunity 53(4):878-894.e7, 2020). This chapter provides details and tips on how to perform Seq-Well S3, along with general pointers on how to subsequently analyze the resultant single-cell RNA-seq data.

A critical assessment of clustering algorithms to improve cell clustering and identification in single-cell transcriptome study.

Cell clustering is typically the initial step in single-cell RNA sequencing (scRNA-seq) analyses. The performance of clustering considerably impacts the validity and reproducibility of cell identification. A variety of clustering algorithms have been developed for scRNA-seq data. These algorithms generate cell label sets that assign each cell to a cluster. However, different algorithms usually yield different label sets, which can introduce variations in cell-type identification based on the generated label sets. Currently, the performance of these algorithms has not been systematically evaluated in single-cell transcriptome studies. Herein, we performed a critical assessment of seven state-of-the-art clustering algorithms including four deep learning-based clustering algorithms and commonly used methods Seurat, Cosine-based Tanimoto similarity-refined graph for community detection using Leiden's algorithm (CosTaL) and Single-cell consensus clustering (SC3). We used diverse evaluation indices based on 10 different scRNA-seq benchmarks to systematically evaluate their clustering performance. Our results show that CosTaL, Seurat, Deep Embedding for Single-cell Clustering (DESC) and SC3 consistently outperformed Single-Cell Clustering Assessment Framework and scDeepCluster based on nine effectiveness scores. Notably, CosTaL and DESC demonstrated superior performance in clustering specific cell types. The performance of the single-cell Variational Inference tools varied across different datasets, suggesting its sensitivity to certain dataset characteristics. Notably, DESC exhibited promising results for cell subtype identification and capturing cellular heterogeneity. In addition, SC3 requires more memory and exhibits slower computation speed compared to other algorithms for the same dataset. In sum, this study provides useful guidance for selecting appropriate clustering methods in scRNA-seq data analysis.

scVAEBGM: Clustering Analysis of Single-Cell ATAC-seq Data Using a Deep Generative Model.

A surge in research has occurred because of current developments in single-cell technologies. Above all, single-cell Assay for Transposase-Accessible Chromatin with high throughput sequencing (scATAC-seq) is a popular approach of analyzing chromatin accessibility differences at the level of single cell, either within or between groups. As a result, it is critical to examine cell heterogeneity at a previously unseen level and to identify both recognized and unknown cell types. However, with the ever-increasing number of cells engendered by technological development and the characteristics of the data, such as high noise, sparsity and dimension, challenges in distinguishing cell types have emerged. We propose scVAEBGM, which integrates a Variational Autoencoder (VAE) with a Bayesian Gaussian-mixture model (BGM) to process and analyze scATAC-seq data. This method combines and takes benefits of a Bayesian Gaussian mixture model to estimate the number of cell types without determining the cluster number in a beforehand. In other words, the size of the clusters is inferred from the data, thus avoiding biases introduced by subjective assessments when manually determining the size of the clusters. Additionally, the method is more robust to noise and can better represent single-cell data in lower dimensions. We also create a further clustering strategy. It is indicated by experiments that further clustering based on the already completed clustering can improve the clustering accuracy again. We test on six public datasets, and scVAEBGM outperforms various dimension reduction baselines. In downstream applications, scVAEBGM can reveal biological cell types.

Forest Fire Clustering for single-cell sequencing combines iterative label propagation with parallelized Monte Carlo simulations.

In the era of single-cell sequencing, there is a growing need to extract insights from data with clustering methods. Here, we introduce Forest Fire Clustering, an efficient and interpretable method for cell-type discovery from single-cell data. Forest Fire Clustering makes minimal prior assumptions and, different from current approaches, calculates a non-parametric posterior probability that each cell is assigned a cell-type label. These posterior distributions allow for the evaluation of a label confidence for each cell and enable the computation of "label entropies", highlighting transitions along developmental trajectories. Furthermore, we show that Forest Fire Clustering can make robust, inductive inferences in an online-learning context and can readily scale to millions of cells. Finally, we demonstrate that our method outperforms state-of-the-art clustering approaches on diverse benchmarks of simulated and experimental data. Overall, Forest Fire Clustering is a useful tool for rare cell type discovery in large-scale single-cell analysis.

A Markov random field model-based approach for differentially expressed gene detection from single-cell RNA-seq data.

The development of single-cell RNA-sequencing (scRNA-seq) technologies has offered insights into complex biological systems at the single-cell resolution. In particular, these techniques facilitate the identifications of genes showing cell-type-specific differential expressions (DE). In this paper, we introduce MARBLES, a novel statistical model for cross-condition DE gene detection from scRNA-seq data. MARBLES employs a Markov Random Field model to borrow information across similar cell types and utilizes cell-type-specific pseudobulk count to account for sample-level variability. Our simulation results showed that MARBLES is more powerful than existing methods to detect DE genes with an appropriate control of false positive rate. Applications of MARBLES to real data identified novel disease-related DE genes and biological pathways from both a single-cell lipopolysaccharide mouse dataset with 24 381 cells and 11 076 genes and a Parkinson's disease human data set with 76 212 cells and 15 891 genes. Overall, MARBLES is a powerful tool to identify cell-type-specific DE genes across conditions from scRNA-seq data.

Regulatory analysis of single cell multiome gene expression and chromatin accessibility data with scREG.

Technological development has enabled the profiling of gene expression and chromatin accessibility from the same cell. We develop scREG, a dimension reduction methodology, based on the concept of cis-regulatory potential, for single cell multiome data. This concept is further used for the construction of subpopulation-specific cis-regulatory networks. The capability of inferring useful regulatory network is demonstrated by the two-fold increment on network inference accuracy compared to the Pearson correlation-based method and the 27-fold enrichment of GWAS variants for inflammatory bowel disease in the cis-regulatory elements. The R package scREG provides comprehensive functions for single cell multiome data analysis.

MarkovHC: Markov hierarchical clustering for the topological structure of high-dimensional single-cell omics data with transition pathway and critical point detection.

Clustering cells and depicting the lineage relationship among cell subpopulations are fundamental tasks in single-cell omics studies. However, existing analytical methods face challenges in stratifying cells, tracking cellular trajectories, and identifying critical points of cell transitions. To overcome these, we proposed a novel Markov hierarchical clustering algorithm (MarkovHC), a topological clustering method that leverages the metastability of exponentially perturbed Markov chains for systematically reconstructing the cellular landscape. Briefly, MarkovHC starts with local connectivity and density derived from the input and outputs a hierarchical structure for the data. We firstly benchmarked MarkovHC on five simulated datasets and ten public single-cell datasets with known labels. Then, we used MarkovHC to investigate the multi-level architectures and transition processes during human embryo preimplantation development and gastric cancer procession. MarkovHC found heterogeneous cell states and sub-cell types in lineage-specific progenitor cells and revealed the most possible transition paths and critical points in the cellular processes. These results demonstrated MarkovHC's effectiveness in facilitating the stratification of cells, identification of cell populations, and characterization of cellular trajectories and critical points.

Molecular Fingerprinting by Single Cell Clone Analysis in Adverse Drug Reaction (ADR) Assessment.

Causality assessment for idiosyncratic ADRs mainly relies on epidemiology, signal detection and less often on proven or plausible mechanistic evidence of the drug at a cellular or organ level. Distinct clones of cells can exist within organs of individual patients, some conferring susceptibility to well-recognised Adverse Drug Reactions (ADRs). Recent advances in molecular biology have allowed the development of single-cell clonal techniques, including single-cell RNA sequencing (scRNA-seq) to molecularly fingerprint ADRs and distinguish between distinct clones of cells within organs in individuals, which may confer differing susceptibilities to ADRs. ScRNA- seq permits molecular fingerprinting of some serious ADRs, mainly in the skin, through the identification of Directly Expressed Genes (DEG) of interest within specific clones. Overexpressed DEGs provide an opportunity for targeted treatment strategies to be developed. scRN A-seq could be applied to a number of other ADRs involving tissues that can be biopsied/sampled (including skin, liver, kidney, blood, stem cells) as well as providing a molecular basis for rapid screening of potential therapeutic candidates, which may not otherwise be predictable from a class of toxicity/organ involvement. A framework for putative assessment for ADRs using scRNA-seq is proposed as well as speculating on potential regulatory implications for pharmacovigilance and drug development. Molecular fingerprinting of ADRs using scRNA-seq may allow better targeting for enhanced pharmacovigilance and risk minimisation measures for medicines with appropriate benefit-risk profiles, although cost-effectiveness and other factors, such as frequency/severity of individual ADRs and population differences, will still be relevant.

A Review of Advanced Impedance Biosensors with Microfluidic Chips for Single-Cell Analysis.

Electrical impedance biosensors combined with microfluidic devices can be used to analyze fundamental biological processes for high-throughput analysis at the single-cell scale. These specialized analytical tools can determine the effectiveness and toxicity of drugs with high sensitivity and demonstrate biological functions on a single-cell scale. Because the various parameters of the cells can be measured depending on methods of single-cell trapping, technological development ultimately determine the efficiency and performance of the sensors. Identifying the latest trends in single-cell trapping technologies afford opportunities such as new structural design and combination with other technologies. This will lead to more advanced applications towards improving measurement sensitivity to the desired target. In this review, we examined the basic principles of impedance sensors and their applications in various biological fields. In the next step, we introduced the latest trend of microfluidic chip technology for trapping single cells and summarized the important findings on the characteristics of single cells in impedance biosensor systems that successfully trapped single cells. This is expected to be used as a leading technology in cell biology, pathology, and pharmacological fields, promoting the further understanding of complex functions and mechanisms within individual cells with numerous data sampling and accurate analysis capabilities.

Joint Clustering of Single-Cell Sequencing and Fluorescence In Situ Hybridization Data for Reconstructing Clonal Heterogeneity in Cancers.

Aneuploidy and whole genome duplication (WGD) events are common features of cancers associated with poor outcomes, but the ways they influence trajectories of clonal evolution are poorly understood. Phylogenetic methods for reconstructing clonal evolution from genomic data have proven a powerful tool for understanding how clonal evolution occurs in the process of cancer progression, but extant methods so far have limited the ability to resolve tumor evolution via ploidy changes. This limitation exists in part because single-cell DNA-sequencing (scSeq), which has been crucial to developing detailed profiles of clonal evolution, has difficulty in resolving ploidy changes and WGD. Multiplex interphase fluorescence in situ hybridization (miFISH) provides a more unambiguous signal of single-cell ploidy changes but it is limited to profiling small numbers of single markers. Here, we develop a joint clustering method to combine these two data sources with the goal of better resolving ploidy changes in tumor evolution. We develop a probabilistic framework to maximize the probability of latent variables given the pre-clustered datasets, which we optimize via Markov chain Monte Carlo sampling combined with linear regression. We validate the method by using simulated data derived from a glioblastoma (GBM) case profiled by both scSeq and miFISH. We further apply the method to two GBM cases with scSeq and miFISH data by reconstructing a phylogenetic tree from the joint clustering results, demonstrating their synergistic value in understanding how focal copy number changes and WGD events can collectively contribute to tumor progression.

A Markov random field model for network-based differential expression analysis of single-cell RNA-seq data.

BACKGROUND: Recent development of single cell sequencing technologies has made it possible to identify genes with different expression (DE) levels at the cell type level between different groups of samples. In this article, we propose to borrow information through known biological networks to increase statistical power to identify differentially expressed genes (DEGs). RESULTS: We develop MRFscRNAseq, which is based on a Markov random field (MRF) model to appropriately accommodate gene network information as well as dependencies among cell types to identify cell-type specific DEGs. We implement an Expectation-Maximization (EM) algorithm with mean field-like approximation to estimate model parameters and a Gibbs sampler to infer DE status. Simulation study shows that our method has better power to detect cell-type specific DEGs than conventional methods while appropriately controlling type I error rate. The usefulness of our method is demonstrated through its application to study the pathogenesis and biological processes of idiopathic pulmonary fibrosis (IPF) using a single-cell RNA-sequencing (scRNA-seq) data set, which contains 18,150 protein-coding genes across 38 cell types on lung tissues from 32 IPF patients and 28 normal controls. CONCLUSIONS: The proposed MRF model is implemented in the R package MRFscRNAseq available on GitHub. By utilizing gene-gene and cell-cell networks, our method increases statistical power to detect differentially expressed genes from scRNA-seq data.