Explore the vast range of titles available.

Multiplexed single-cell RNA-seq analysis of multiple samples using pooling is a promising experimental design, offering increased throughput while allowing to overcome batch variation. To reconstruct the sample identify of each cell, genetic variants that segregate between the samples in the pool have been proposed as natural barcode for cell demultiplexing. Existing demultiplexing strategies rely on availability of complete genotype data from the pooled samples, which limits the applicability of such methods, in particular when genetic variation is not the primary object of study. To address this, we here present Vireo, a computationally efficient Bayesian model to demultiplex single-cell data from pooled experimental designs. Uniquely, our model can be applied in settings when only partial or no genotype information is available. Using pools based on synthetic mixtures and results on real data, we demonstrate the robustness of Vireo and illustrate the utility of multiplexed experimental designs for common expression analyses.

Cell-cell communication is a key aspect of dissecting the complex cellular microenvironment. Existing single-cell and spatial transcriptomics-based methods primarily focus on identifying cell-type pairs for a specific interaction, while less attention has been paid to the prioritisation of interaction features or the identification of interaction spots in the spatial context. Here, we introduce SpatialDM, a statistical model and toolbox leveraging a bivariant Moran’s statistic to detect spatially co-expressed ligand and receptor pairs, their local interacting spots (single-spot resolution), and communication patterns. By deriving an analytical null distribution, this method is scalable to millions of spots and shows accurate and robust performance in various simulations. On multiple datasets including melanoma, Ventricular-Subventricular Zone, and intestine, SpatialDM reveals promising communication patterns and identifies differential interactions between conditions, hence enabling the discovery of context-specific cell cooperation and signalling. Spatial omics are increasingly being recognised to study cell-cell communications. Here, the authors present a bioinformatics toolbox for rapid identification of spatially co-expressed ligand-receptor and revealing cell-cell communication patterns.

The recent breakthrough of single-cell RNA velocity methods brings attractive promises to reveal directed trajectory on cell differentiation, states transition and response to perturbations. However, the existing RNA velocity methods are often found to return erroneous results, partly due to model violation or lack of temporal regularization. Here, we present UniTVelo, a statistical framework of RNA velocity that models the dynamics of spliced and unspliced RNAs via flexible transcription activities. Uniquely, it also supports the inference of a unified latent time across the transcriptome. With ten datasets, we demonstrate that UniTVelo returns the expected trajectory in different biological systems, including hematopoietic differentiation and those even with weak kinetics or complex branches. RNA velocity can detect the differentiation directionality by modelling sparse unspliced RNAs, but suffers from high estimation errors. Here, the authors develop a computational method called UniTVelo to reinforce the velocity estimation by introducing a unified time and a top-down model design.

RNA splicing is an important driver of heterogeneity in single cells through the expression of alternative transcripts and as a determinant of transcriptional kinetics. However, the intrinsic coverage limitations of scRNA-seq technologies make it challenging to associate specific splicing events to cell-level phenotypes. BRIE2 is a scalable computational method that resolves these issues by regressing single-cell transcriptomic data against cell-level features. We show that BRIE2 effectively identifies differential disease-associated alternative splicing events and allows a principled selection of genes that capture heterogeneity in transcriptional kinetics and improve RNA velocity analyses, enabling the identification of splicing phenotypes associated with biological changes.

Single-cell RNA-seq (scRNA-seq) provides a comprehensive measurement of stochasticity in transcription, but the limitations of the technology have prevented its application to dissect variability in RNA processing events such as splicing. Here, we present BRIE (Bayesian regression for isoform estimation), a Bayesian hierarchical model that resolves these problems by learning an informative prior distribution from sequence features. We show that BRIE yields reproducible estimates of exon inclusion ratios in single cells and provides an effective tool for differential isoform quantification between scRNA-seq data sets. BRIE, therefore, expands the scope of scRNA-seq experiments to probe the stochasticity of RNA processing.

The first-strand reads (often reads 1) of three-prime single-cell RNA-seq (3′ scRNA-seq) can contain informative cDNA for analysis of polyadenylation sites (PAS), but are often overlooked or trimmed. Here, we describe a computational method, scTail, to identify PAS using first-strand reads and quantify its expression leveraging second-strand reads, consequently enabling detection of alternative PAS usage. Compared with other methods, scTail detects PAS more precisely and retains high sensitivity. Furthermore, we demonstrated that combining scTail and BRIE2 can discover differential alternative PAS usage in various biological processes including cancers and time-series development, giving critical insight into PAS regulation.

Single-nucleotide-variant (SNV) clone assignment of high-covariance single-cell lineage tracing data remains a challenge due to hierarchical mutation structure and many missing signals. We develop SNPmanifold, a Python package that learns an SNV embedding manifold using a binomial variational autoencoder to give an efficient and interpretable cell-cell distance metric. We demonstrate that SNPmanifold is a suitable tool for analysis of complex, single-cell SNV mutation data, such as in the context of demultiplexing a large number of donors and somatic lineage tracing via mitochondrial SNV data and can reveal insights into single-cell clonality and lineages more accurately and comprehensively than existing methods.

Five-prime single-cell RNA-seq (scRNA-seq) has been widely employed to profile cellular transcriptomes, however, its power of analysing transcription start sites (TSS) has not been fully utilised. Here, we present a computational method suite, CamoTSS, to precisely identify TSS and quantify its expression by leveraging the cDNA on read 1, which enables effective detection of alternative TSS usage. With various experimental data sets, we have demonstrated that CamoTSS can accurately identify TSS and the detected alternative TSS usages showed strong specificity in different biological processes, including cell types across human organs, the development of human thymus, and cancer conditions. As evidenced in nasopharyngeal cancer, alternative TSS usage can also reveal regulatory patterns including systematic TSS dysregulations. Five-prime single-cell RNA-seq, especially the read 1, has precise capture of transcription start sites (TSS), but such information is often overlooked. Here, authors present a computational method suite, CamoTSS, to precisely identify TSS and quantify its expression, enabling effective detection of alternative TSS usage in different biological processes.

Mitochondrial mutations are increasingly recognised as informative endogenous genetic markers that can be used to reconstruct cellular clonal structure using single-cell RNA or DNA sequencing data. However, identifying informative mtDNA variants in noisy and sparse single-cell sequencing data is still challenging with few computation methods available. Here we present an open source computational tool MQuad that accurately calls clonally informative mtDNA variants in a population of single cells, and an analysis suite for complete clonality inference, based on single cell RNA, DNA or ATAC sequencing data. Through a variety of simulated and experimental single cell sequencing data, we showed that MQuad can identify mitochondrial variants with both high sensitivity and specificity, outperforming existing methods by a large extent. Furthermore, we demonstrate its wide applicability in different single cell sequencing protocols, particularly in complementing single-nucleotide and copy-number variations to extract finer clonal resolution. Mitochondrial variants are informative endogenous barcodes for clonal substructure. Here, the authors developed a computational method MQuad to effectively detect these clonal informed mtDNA variants from single-cell RNA, DNA or ATAC sequencing data.

Understanding human endometrial dynamics in the establishment of endometrial receptivity remains a challenge, which limits early diagnosis and treatment of endometrial-factor infertility. Here, we decode the endometrial dynamics of fertile women across the window of implantation and characterize the endometrial deficiency in women with recurrent implantation failure. A computational model capable of both temporal prediction and pattern discovery is used to analyze single-cell transcriptomic data from over 220,000 endometrial cells. The time-series atlas highlights a two-stage stromal decidualization process and a gradual transitional process of the luminal epithelial cells across the window of implantation. In addition, a time-varying gene set regulating epithelium receptivity is identified, based on which the recurrent implantation failure endometria are stratified into two classes of deficiencies. Further investigation uncovers a hyper-inflammatory microenvironment for the dysfunctional endometrial epithelial cells of recurrent implantation failure. The holistic characterization of the physiological and pathophysiological window of implantation and a computational tool trained on this temporal atlas provide a platform for future therapeutic developments. Time-series single-cell transcriptomic characterization of luteal-phase endometrium in fertile women with a variational autoencoder model uncovers the association of decreased epithelial receptivity and hyperinflammatory microenvironment in recurrent implantation failures.