simpleaf: a simple, flexible, and scalable framework for single-cell data processing using alevin-fry

T-bet+ CXCR3+ B cells drive hyperreactive B-T cell interactions in multiple sclerosis

Multiple sclerosis (MS) is an autoimmune disease of the central nervous system (CNS). Self-peptide-dependent autoproliferation (AP) of B and T cells is a key mechanism in MS. Here, we show that pro-inflammatory B-T cell-enriched cell clusters (BTECs) form during AP and mirror features of a germinal center reaction. T-bet+CXCR3+ B cells are the main cell subset amplifying and sustaining their counterpart Th1 cells via interferon (IFN)-γ and are present in highly inflamed meningeal tissue. The underlying B cell activation signature is reflected by epigenetic modifications and receptor-ligand interactions with self-reactive T cells. AP+ CXCR3+ B cells show marked clonal evolution from memory to somatically hypermutated plasmablasts and upregulation of IFN-γ-related genes. Our data underscore a key role of T-bet+CXCR3+ B cells in the pathogenesis of MS in both the peripheral immune system and the CNS compartment, and thus they appear to be involved in both early relapsing-remitting disease and the chronic stage.

kallisto, bustools and kb-python for quantifying bulk, single-cell and single-nucleus RNA-seq

The term 'RNA-seq' refers to a collection of assays based on sequencing experiments that involve quantifying RNA species from bulk tissue, single cells or single nuclei. The kallisto, bustools and kb-python programs are free, open-source software tools for performing this analysis that together can produce gene expression quantification from raw sequencing reads. The quantifications can be individualized for multiple cells, multiple samples or both. Additionally, these tools allow gene expression values to be classified as originating from nascent RNA species or mature RNA species, making this workflow amenable to both cell-based and nucleus-based assays. This protocol describes in detail how to use kallisto and bustools in conjunction with a wrapper, kb-python, to preprocess RNA-seq data. Execution of this protocol requires basic familiarity with a command line environment. With this protocol, quantification of a moderately sized RNA-seq dataset can be completed within minutes.

Accurate quantification of nascent and mature RNAs from single-cell and single-nucleus RNA-seq

In single-cell and single-nucleus RNA sequencing (RNA-seq), the coexistence of nascent (unprocessed) and mature (processed) messenger RNA (mRNA) poses challenges in accurate read mapping and the interpretation of count matrices. The traditional transcriptome reference, defining the "region of interest" in bulk RNA-seq, restricts its focus to mature mRNA transcripts. This restriction leads to two problems: reads originating outside of the "region of interest" are prone to mismapping within this region, and additionally, such external reads cannot be matched to specific transcript targets. Expanding the "region of interest" to encompass both nascent and mature mRNA transcript targets provides a more comprehensive framework for RNA-seq analysis. Here, we introduce the concept of distinguishing flanking k-mers (DFKs) to improve mapping of sequencing reads. We have developed an algorithm to identify DFKs, which serve as a sophisticated "background filter", enhancing the accuracy of mRNA quantification. This dual strategy of an expanded region of interest coupled with the use of DFKs enhances the precision in quantifying both mature and nascent mRNA molecules, as well as in delineating reads of ambiguous status.

Optimizing single-cell RNA sequencing methods for human colon biopsies: droplet-based vs. picowell-based platforms

Background & Aims

Single-cell RNA sequencing (scRNA) has empowered many insights into gastrointestinal microenvironments. However, profiling human biopsies using droplet-based scRNA (D-scRNA) is challenging since it requires immediate processing to minimize epithelial cell damage. In contrast, picowell-based (P-scRNA) platforms permit short-term frozen storage before sequencing. We compared P- and D-scRNA platforms on cells derived from human colon biopsies.

Methods

Endoscopic rectosigmoid mucosal biopsies were obtained from two adults and conducted D-scRNA (10X Chromium) and P-scRNA (Honeycomb HIVE) in parallel using an individual's pool of single cells (> 10,000 cells/participant). Three experiments were performed to evaluate 1) P-scRNA with cells under specific storage conditions (immediately processed [fresh], vs. frozen at -20C vs. -80C [2 weeks]); 2) fresh P-scRNA versus fresh D-scRNA; and 3) P-scRNA stored at -80C with fresh D-scRNA.

Results

Significant recovery of loaded cells was achieved for fresh (80.9%) and -80C (48.5%) P-scRNA and D-scRNA (76.6%), but not -20C P-scRNA (3.7%). However, D-scRNA captures more typeable cells among recovered cells (71.5% vs. 15.8% Fresh and 18.4% -80C P-scRNA), and these cells exhibit higher gene coverage at the expense of higher mitochondrial read fractions across most cell types. Cells profiled using D-scRNA demonstrated more consistent gene expression profiles among the same cell type than those profiled using P-scRNA. Significant intra-cell-type differences were observed in profiled gene classes across platforms.

Conclusions

Our results highlight non-overlapping advantages of P-scRNA and D-scRNA and underscore the need for innovation to enable high-fidelity capture of colonic epithelial cells. The platform-specific variation highlights the challenges of maintaining rigor and reproducibility across studies that use different platforms.

Forseti: a mechanistic and predictive model of the splicing status of scRNA-seq reads

MOTIVATION

Short-read single-cell RNA-sequencing (scRNA-seq) has been used to study cellular heterogeneity, cellular fate, and transcriptional dynamics. Modeling splicing dynamics in scRNA-seq data is challenging, with inherent difficulty in even the seemingly straightforward task of elucidating the splicing status of the molecules from which sequenced fragments are drawn. This difficulty arises, in part, from the limited read length and positional biases, which substantially reduce the specificity of the sequenced fragments. As a result, the splicing status of many reads in scRNA-seq is ambiguous because of a lack of definitive evidence. We are therefore in need of methods that can recover the splicing status of ambiguous reads which, in turn, can lead to more accuracy and confidence in downstream analyses.

RESULTS

We develop Forseti, a predictive model to probabilistically assign a splicing status to scRNA-seq reads. Our model has two key components. First, we train a binding affinity model to assign a probability that a given transcriptomic site is used in fragment generation. Second, we fit a robust fragment length distribution model that generalizes well across datasets deriving from different species and tissue types. Forseti combines these two trained models to predict the splicing status of the molecule of origin of reads by scoring putative fragments that associate each alignment of sequenced reads with proximate potential priming sites. Using both simulated and experimental data, we show that our model can precisely predict the splicing status of many reads and identify the true gene origin of multi-gene mapped reads.

AVAILABILITY AND IMPLEMENTATION

Forseti and the code used for producing the results are available at https://github.com/COMBINE-lab/forseti under a BSD 3-clause license.

scCensus: Off-target scRNA-seq reads reveal meaningful biology

Single-cell RNA-sequencing (scRNA-seq) provides unprecedented insights into cellular heterogeneity. Although scRNA-seq reads from most prevalent and popular tagged-end protocols are expected to arise from the 3' end of polyadenylated RNAs, recent studies have shown that "off-target" reads can constitute a substantial portion of the read population. In this work, we introduced scCensus, a comprehensive analysis workflow for systematically evaluating and categorizing off-target reads in scRNA-seq. We applied scCensus to seven scRNA-seq datasets. Our analysis of intergenic reads shows that these off-target reads contain information about chromatin structure and can be used to identify similar cells across modalities. Our analysis of antisense reads suggests that these reads can be used to improve gene detection and capture interesting transcriptional activities like antisense transcription. Furthermore, using splice-aware quantification, we find that spliced and unspliced reads provide distinct information about cell clusters and biomarkers, suggesting the utility of integrating signals from reads with different splicing statuses. Overall, our results suggest that off-target scRNA-seq reads contain underappreciated information about various transcriptional activities. These observations about yet-unexploited information in existing scRNA-seq data will help guide and motivate the community to improve current algorithms and analysis methods, and to develop novel approaches that utilize off-target reads to extend the reach and accuracy of single-cell data analysis pipelines.

kallisto, bustools, and kb-python for quantifying bulk, single-cell, and single-nucleus RNA-seq

The term "RNA-seq" refers to a collection of assays based on sequencing experiments that involve quantifying RNA species from bulk tissue, from single cells, or from single nuclei. The kallisto, bustools, and kb-python programs are free, open-source software tools for performing this analysis that together can produce gene expression quantification from raw sequencing reads. The quantifications can be individualized for multiple cells, multiple samples, or both. Additionally, these tools allow gene expression values to be classified as originating from nascent RNA species or mature RNA species, making this workflow amenable to both cell-based and nucleus-based assays. This protocol describes in detail how to use kallisto and bustools in conjunction with a wrapper, kb-python, to preprocess RNA-seq data.