Inference and effects of barcode multiplets in droplet-based single-cell assays

Droplet-based single-cell sequencing: Strategies and applications

Notable advancements in single-cell omics technologies have not only addressed longstanding challenges but also enabled unprecedented studies of cellular heterogeneity with unprecedented resolution and scale. These strides have led to groundbreaking insights into complex biological systems, paving the way for a more profound comprehension of human biology and diseases. The droplet microfluidic technology has become a crucial component in many single-cell sequencing workflows in terms of throughput, cost-effectiveness, and automation. Utilizing a microfluidic chip to encapsulate and profile individual cells within droplets has significantly improved single-cell research. Therefore, this review aims to comprehensively elaborate the droplet microfluidics-assisted omics methods from a single-cell perspective. The strategies for using droplet microfluidics in the realms of genomics, epigenomics, transcriptomics, and proteomics analyses are first introduced. On this basis, the focus then turns to the latest applications of this technology in different sequencing patterns, including mono- and multi-omics. Finally, the challenges and further perspectives of droplet-based single-cell sequencing in both foundational research and commercial applications are discussed.

High-throughput gene expression analysis with TempO-LINC sensitively resolves complex brain, lung and kidney heterogeneity at single-cell resolution

We report the development and performance of a novel genomics platform, TempO-LINC, for conducting high-throughput transcriptomic analysis on single cells and nuclei. TempO-LINC works by adding cell-identifying molecular barcodes onto highly selective and high-sensitivity gene expression probes within fixed cells, without having to first generate cDNA. Using an instrument-free combinatorial-indexing approach, all probes within the same fixed cell receive an identical barcode, enabling the reconstruction of single-cell gene expression profiles across as few as several hundred cells and up to 100,000+ cells per run. The TempO-LINC approach is easily scalable based on the number of barcodes and rounds of barcoding performed; however, for the experiments reported in this study, the assay utilized over 5.3 million unique barcodes. TempO-LINC has a robust protocol for fixing and banking cells and displays high-sensitivity gene detection from multiple diverse sample types. We show that TempO-LINC has an observed multiplet rate of less than 1.1% and a cell capture rate of ~50%. Although the assay can accurately profile the whole transcriptome (19,683 human or 21,400 mouse genes), it can be targeted to measure only actionable/informative genes and molecular pathways of interest - thereby reducing sequencing requirements. In this study, we applied TempO-LINC to profile the transcriptomes of 89,722 cells across multiple sample types, including nuclei from mouse lung, kidney and brain tissues. The data demonstrated the ability to identify and annotate at least 50 unique cell populations and positively correlate expression of cell type-specific molecular markers within them. TempO-LINC is a robust new single-cell technology that is ideal for large-scale applications/studies across thousands of samples with high data quality.

Concepts and new developments in droplet-based single cell multi-omics

Single cell sequencing technologies have become a fixture in the molecular profiling of cells due to their ease, flexibility, and commercial availability. In particular, partitioning individual cells inside oil droplets via microfluidic reactions enables transcriptomic or multi-omic measurements for thousands of cells in parallel. Complementing the multitude of biological discoveries from genomics analyses, the past decade has brought new capabilities from assay baselines to enable a deeper understanding of the complex data from single cell multi-omics. Here, we highlight four innovations that have improved the reliability and understanding of droplet microfluidic assays. We emphasize new developments that further orient principles of technology development and guidelines for the design, benchmarking, and implementation of new droplet-based methodologies.

Single-Cell Chromatin Accessibility Analysis Reveals Subgroup-Specific TF-NTR Regulatory Circuits in Medulloblastoma

Medulloblastoma (MB) stands as one of the prevalent malignant brain tumors among pediatric patients. Despite its prevalence, the intricate interplay between the regulatory program driving malignancy in MB cells and their interactions with the microenvironment remains insufficiently understood. Leveraging the capabilities of single-cell Assay for Transposase-Accessible Chromatin sequencing (scATAC-seq), the chromatin accessibility landscape is unveiled across 59,015 distinct MB cells. This expansive dataset encompasses cells belonging to discrete molecular subgroups, namely SHH, WNT, Group3, and Group4. Within these chromatin accessibility profiles, specific regulatory elements tied to individual subgroups are uncovered, shedding light on the distinct activities of transcription factors (TFs) that likely orchestrate the tumorigenesis process. Moreover, it is found that certain neurotransmitter receptors (NTRs) are subgroup-specific and can predict MB subgroup classification when combined with their associated transcription factors. Notably, targeting essential NTRs within tumors influences both the in vitro sphere-forming capability and the in vivo tumorigenic capacity of MB cells. These findings collectively provide fresh insights into comprehending the regulatory networks and cellular dynamics within MBs. Furthermore, the significance of the TF-NTR regulatory circuits is underscored as prospective biomarkers and viable therapeutic targets.

Systematic benchmarking of single-cell ATAC-sequencing protocols

Single-cell assay for transposase-accessible chromatin by sequencing (scATAC-seq) has emerged as a powerful tool for dissecting regulatory landscapes and cellular heterogeneity. However, an exploration of systemic biases among scATAC-seq technologies has remained absent. In this study, we benchmark the performance of eight scATAC-seq methods across 47 experiments using human peripheral blood mononuclear cells (PBMCs) as a reference sample and develop PUMATAC, a universal preprocessing pipeline, to handle the various sequencing data formats. Our analyses reveal significant differences in sequencing library complexity and tagmentation specificity, which impact cell-type annotation, genotype demultiplexing, peak calling, differential region accessibility and transcription factor motif enrichment. Our findings underscore the importance of sample extraction, method selection, data processing and total cost of experiments, offering valuable guidance for future research. Finally, our data and analysis pipeline encompasses 169,000 PBMC scATAC-seq profiles and a best practices code repository for scATAC-seq data analysis, which are freely available to extend this benchmarking effort to future protocols.

Tracking in situ checkpoint inhibitor-bound target T cells in patients with checkpoint-induced colitis

The success of checkpoint inhibitors (CPIs) for cancer has been tempered by immune-related adverse effects including colitis. CPI-induced colitis is hallmarked by expansion of resident mucosal IFNγ cytotoxic CD8+ T cells, but how these arise is unclear. Here, we track CPI-bound T cells in intestinal tissue using multimodal single-cell and subcellular spatial transcriptomics (ST). Target occupancy was increased in inflamed tissue, with drug-bound T cells located in distinct microdomains distinguished by specific intercellular signaling and transcriptional gradients. CPI-bound cells were largely CD4+ T cells, including enrichment in CPI-bound peripheral helper, follicular helper, and regulatory T cells. IFNγ CD8+ T cells emerged from both tissue-resident memory (TRM) and peripheral populations, displayed more restricted target occupancy profiles, and co-localized with damaged epithelial microdomains lacking effective regulatory cues. Our multimodal analysis identifies causal pathways and constitutes a resource to inform novel preventive strategies.

txci-ATAC-seq: a massive-scale single-cell technique to profile chromatin accessibility

We develop a large-scale single-cell ATAC-seq method by combining Tn5-based pre-indexing with 10× Genomics barcoding, enabling the indexing of up to 200,000 nuclei across multiple samples in a single reaction. We profile 449,953 nuclei across diverse tissues, including the human cortex, mouse brain, human lung, mouse lung, mouse liver, and lung tissue from a club cell secretory protein knockout (CC16-/-) model. Our study of CC16-/- nuclei uncovers previously underappreciated technical artifacts derived from remnant 129 mouse strain genetic material, which cause profound cell-type-specific changes in regulatory elements near many genes, thereby confounding the interpretation of this commonly referenced mouse model.

Comparison of Single Cell Transcriptome Sequencing Methods: Of Mice and Men

Single cell RNAseq has been a big leap in many areas of biology. Rather than investigating gene expression on a whole organism level, this technology enables scientists to get a detailed look at rare single cells or within their cell population of interest. The field is growing, and many new methods appear each year. We compared methods utilized in our core facility: Smart-seq3, PlexWell, FLASH-seq, VASA-seq, SORT-seq, 10X, Evercode, and HIVE. We characterized the equipment requirements for each method. We evaluated the performances of these methods based on detected features, transcriptome diversity, mitochondrial RNA abundance and multiplets, among others and benchmarked them against bulk RNA sequencing. Here, we show that bulk transcriptome detects more unique transcripts than any single cell method. While most methods are comparable in many regards, FLASH-seq and VASA-seq yielded the best metrics, e.g., in number of features. If no equipment for automation is available or many cells are desired, then HIVE or 10X yield good results. In general, more recently developed methods perform better. This also leads to the conclusion that older methods should be phased out, and that the development of single cell RNAseq methods is still progressing considerably.

BLR: a flexible pipeline for haplotype analysis of multiple linked-read technologies

Linked-read sequencing promises a one-method approach for genome-wide insights including single nucleotide variants (SNVs), structural variants, and haplotyping. We introduce Barcode Linked Reads (BLR), an open-source haplotyping pipeline capable of handling millions of barcodes and data from multiple linked-read technologies including DBS, 10× Genomics, TELL-seq and stLFR. Running BLR on DBS linked-reads yielded megabase-scale phasing with low (<0.2%) switch error rates. Of 13616 protein-coding genes phased in the GIAB benchmark set (v4.2.1), 98.6% matched the BLR phasing. In addition, large structural variants showed concordance with HPRC-HG002 reference assembly calls. Compared to diploid assembly with PacBio HiFi reads, BLR phasing was more continuous when considering switch errors. We further show that integrating long reads at low coverage (∼10×) can improve phasing contiguity and reduce switch errors in tandem repeats. When compared to Long Ranger on 10× Genomics data, BLR showed an increase in phase block N50 with low switch-error rates. For TELL-Seq and stLFR linked reads, BLR generated longer or similar phase block lengths and low switch error rates compared to results presented in the original publications. In conclusion, BLR provides a flexible workflow for comprehensive haplotype analysis of linked reads from multiple platforms.

scIBD: a self-supervised iterative-optimizing model for boosting the detection of heterotypic doublets in single-cell chromatin accessibility data

Application of the widely used droplet-based microfluidic technologies in single-cell sequencing often yields doublets, introducing bias to downstream analyses. Especially, doublet-detection methods for single-cell chromatin accessibility sequencing (scCAS) data have multiple assay-specific challenges. Therefore, we propose scIBD, a self-supervised iterative-optimizing model for boosting heterotypic doublet detection in scCAS data. scIBD introduces an adaptive strategy to simulate high-confident heterotypic doublets and self-supervise for doublet-detection in an iteratively optimizing manner. Comprehensive benchmarking on various simulated and real datasets demonstrates the outperformance and robustness of scIBD. Moreover, the downstream biological analyses suggest the efficacy of doublet-removal by scIBD.

Gene-regulation modules in nonalcoholic fatty liver disease revealed by single-nucleus ATAC-seq

We investigated the progression of nonalcoholic fatty liver disease from fatty liver to steatohepatitis using single-nucleus and bulk ATAC-seq on the livers of rats fed a high-fat diet (HFD). Rats fed HFD for 4 wk developed fatty liver, and those fed HFD for 8 wk further progressed to steatohepatitis. We observed an increase in the proportion of inflammatory macrophages, consistent with the pathological progression. Utilizing machine learning, we divided global gene regulation into modules, wherein transcription factors within a module could regulate genes within the same module, reaffirming known regulatory relationships between transcription factors and biological processes. We identified core genes-central to co-expression and protein-protein interaction-for the biological processes discovered. Notably, a large part of the core genes overlapped with genes previously implicated in nonalcoholic fatty liver disease. Single-nucleus ATAC-seq, combined with data-driven statistical analysis, offers insight into in vivo global gene regulation as a combination of modules and assists in identifying core genes of relevant biological processes.

Divergent single cell transcriptome and epigenome alterations in ALS and FTD patients with C9orf72 mutation

A repeat expansion in the C9orf72 (C9) gene is the most common genetic cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). Here we investigate single nucleus transcriptomics (snRNA-seq) and epigenomics (snATAC-seq) in postmortem motor and frontal cortices from C9-ALS, C9-FTD, and control donors. C9-ALS donors present pervasive alterations of gene expression with concordant changes in chromatin accessibility and histone modifications. The greatest alterations occur in upper and deep layer excitatory neurons, as well as in astrocytes. In neurons, the changes imply an increase in proteostasis, metabolism, and protein expression pathways, alongside a decrease in neuronal function. In astrocytes, the alterations suggest activation and structural remodeling. Conversely, C9-FTD donors have fewer high-quality neuronal nuclei in the frontal cortex and numerous gene expression changes in glial cells. These findings highlight a context-dependent molecular disruption in C9-ALS and C9-FTD, indicating unique effects across cell types, brain regions, and diseases.

DUPLETS: Deformability-Assisted Dual-Particle Encapsulation Via Electrically Activated Sorting

Co-encapsulation of bead carriers and biological cells in microfluidics has become a powerful technique for various biological assays in single-cell genomics and drug screening because of its distinct capability of single-cell confinement. However, current co-encapsulation approaches exist a trade-off between cell/bead pairing rate and probability of multiple cells in individual droplets, significantly limiting the effective throughput of single-paired cell-bead droplets production. Deformability-assisted dUal-Particle encapsuLation via Electrically acTivated Sorting (DUPLETS) system is reported to overcome this problem. The DUPLETS can differentiate the encapsulated content in individual droplets and sort out targeted droplets via a combined screening of mechanical and electrical characteristics of single droplets in label-free manners and with the highest effective throughput in comparison to current commercial platforms. The DUPLETS has been demonstrated to enrich single-paired cell-bead droplets to over 80% (above eightfold higher than current co-encapsulation techniques). It eliminates multicell droplets to 0.1% whereas up to ≈24% in 10× Chromium. It is believed that merging DUPLETS into the current co-encapsulation platforms can meaningfully elevate sample quality in terms of high purity of single-paired cell-bead droplets, low fraction of multicell droplets, and high cell viability, which can benefit a multitude of biological assay applications.

Opportunities and challenges in the application of single-cell and spatial transcriptomics in plants

Single-cell and spatial transcriptomics have diverted researchers' attention from the multicellular level to the single-cell level and spatial information. Single-cell transcriptomes provide insights into the transcriptome at the single-cell level, whereas spatial transcriptomes help preserve spatial information. Although these two omics technologies are helpful and mature, further research is needed to ensure their widespread applicability in plant studies. Reviewing recent research on plant single-cell or spatial transcriptomics, we compared the different experimental methods used in various plants. The limitations and challenges are clear for both single-cell and spatial transcriptomic analyses, such as the lack of applicability, spatial information, or high resolution. Subsequently, we put forth further applications, such as cross-species analysis of roots at the single-cell level and the idea that single-cell transcriptome analysis needs to be combined with other omics analyses to achieve superiority over individual omics analyses. Overall, the results of this review suggest that combining single-cell transcriptomics, spatial transcriptomics, and spatial element distribution can provide a promising research direction, particularly for plant research.

Molecular features driving cellular complexity of human brain evolution

Human-specific genomic changes contribute to the unique functionalities of the human brain1-5. The cellular heterogeneity of the human brain6,7 and the complex regulation of gene expression highlight the need to characterize human-specific molecular features at cellular resolution. Here we analysed single-nucleus RNA-sequencing and single-nucleus assay for transposase-accessible chromatin with sequencing datasets for human, chimpanzee and rhesus macaque brain tissue from posterior cingulate cortex. We show a human-specific increase of oligodendrocyte progenitor cells and a decrease of mature oligodendrocytes across cortical tissues. Human-specific regulatory changes were accelerated in oligodendrocyte progenitor cells, and we highlight key biological pathways that may be associated with the proportional changes. We also identify human-specific regulatory changes in neuronal subtypes, which reveal human-specific upregulation of FOXP2 in only two of the neuronal subtypes. We additionally identify hundreds of new human accelerated genomic regions associated with human-specific chromatin accessibility changes. Our data also reveal that FOS::JUN and FOX motifs are enriched in the human-specifically accessible chromatin regions of excitatory neuronal subtypes. Together, our results reveal several new mechanisms underlying the evolutionary innovation of human brain at cell-type resolution.

Integrating genetics with single-cell multiomic measurements across disease states identifies mechanisms of beta cell dysfunction in type 2 diabetes

Dysfunctional pancreatic islet beta cells are a hallmark of type 2 diabetes (T2D), but a comprehensive understanding of the underlying mechanisms, including gene dysregulation, is lacking. Here we integrate information from measurements of chromatin accessibility, gene expression and function in single beta cells with genetic association data to nominate disease-causal gene regulatory changes in T2D. Using machine learning on chromatin accessibility data from 34 nondiabetic, pre-T2D and T2D donors, we identify two transcriptionally and functionally distinct beta cell subtypes that undergo an abundance shift during T2D progression. Subtype-defining accessible chromatin is enriched for T2D risk variants, suggesting a causal contribution of subtype identity to T2D. Both beta cell subtypes exhibit activation of a stress-response transcriptional program and functional impairment in T2D, which is probably induced by the T2D-associated metabolic environment. Our findings demonstrate the power of multimodal single-cell measurements combined with machine learning for characterizing mechanisms of complex diseases.

Understanding cell fate acquisition in stem-cell-derived pancreatic islets using single-cell multiome-inferred regulomes

Pancreatic islet cells derived from human pluripotent stem cells hold great promise for modeling and treating diabetes. Differences between stem-cell-derived and primary islets remain, but molecular insights to inform improvements are limited. Here, we acquire single-cell transcriptomes and accessible chromatin profiles during in vitro islet differentiation and pancreas from childhood and adult donors for comparison. We delineate major cell types, define their regulomes, and describe spatiotemporal gene regulatory relationships between transcription factors. CDX2 emerged as a regulator of enterochromaffin-like cells, which we show resemble a transient, previously unrecognized, serotonin-producing pre-β cell population in fetal pancreas, arguing against a proposed non-pancreatic origin. Furthermore, we observe insufficient activation of signal-dependent transcriptional programs during in vitro β cell maturation and identify sex hormones as drivers of β cell proliferation in childhood. Altogether, our analysis provides a comprehensive understanding of cell fate acquisition in stem-cell-derived islets and a framework for manipulating cell identities and maturity.

Single-cell chromatin accessibility and transcriptome atlas of mouse embryos

Cis-regulatory elements regulate gene expression and lineage specification. However, the potential regulation of cis-elements on mammalian embryogenesis remains largely unexplored. To address this question, we perform single-cell assay for transposase-accessible chromatin using sequencing (ATAC-seq) and RNA-seq in embryonic day 7.5 (E7.5) and E13.5 mouse embryos. We construct the chromatin accessibility landscapes with cell spatial information in E7.5 embryos, showing the spatial patterns of cis-elements and the spatial distribution of potentially functional transcription factors (TFs). We further show that many germ-layer-specific cis-elements and TFs in E7.5 embryos are maintained in the cell types derived from the corresponding germ layers at later stages, suggesting that these cis-elements and TFs are important during cell differentiation. We also find a potential progenitor for Sertoli and granulosa cells in gonads. Interestingly, both Sertoli and granulosa cells exist in male gonads and female gonads during gonad development. Collectively, we provide a valuable resource to understand organogenesis in mammals.

Integration of single-cell multiomic measurements across disease states with genetics identifies mechanisms of beta cell dysfunction in type 2 diabetes

Altered function and gene regulation of pancreatic islet beta cells is a hallmark of type 2 diabetes (T2D), but a comprehensive understanding of mechanisms driving T2D is still missing. Here we integrate information from measurements of chromatin activity, gene expression and function in single beta cells with genetic association data to identify disease-causal gene regulatory changes in T2D. Using machine learning on chromatin accessibility data from 34 non-diabetic, pre-T2D and T2D donors, we robustly identify two transcriptionally and functionally distinct beta cell subtypes that undergo an abundance shift in T2D. Subtype-defining active chromatin is enriched for T2D risk variants, suggesting a causal contribution of subtype identity to T2D. Both subtypes exhibit activation of a stress-response transcriptional program and functional impairment in T2D, which is likely induced by the T2D-associated metabolic environment. Our findings demonstrate the power of multimodal single-cell measurements combined with machine learning for identifying mechanisms of complex diseases.

Highly sensitive single-cell chromatin accessibility assay and transcriptome coassay with METATAC

The thriving field of single-cell genomics allows researchers to dissect the complexity and heterogeneity of tissues at single-cell resolution at large scale, involving transcriptome and epigenome. However, single-cell chromatin accessibility profiling methods exhibit low sensitivity. Here, we increased accessible chromatin detection sensitivity in single cells with METATAC, a single-cell ATAC-seq technique, with the help of META amplification strategy and other biochemical modifications. METATAC reached the highest accessible chromatin region detection efficiency compared with existing techniques, allowing more accurate cis-regulatory element coaccessibility measurement and allele-specific chromatin accessibility analysis in complex tissue samples. In combination with a high-resolution single-cell RNA sequencing assay, we further developed a high-sensitivity joint single-cell ATAC–RNA strategy, which helps us to better resolve gene regulatory programs.

Recent advances in single-cell assay for transposase accessible chromatin using sequencing (scATAC-seq) and its coassays have transformed the field of single-cell epigenomics and transcriptomics. However, the low detection efficiency of current methods has limited our understanding of the true complexity of chromatin accessibility and its relationship with gene expression in single cells. Here, we report a high-sensitivity scATAC-seq method, termed multiplexed end-tagging amplification of transposase accessible chromatin (METATAC), which detects a large number of accessible sites per cell and is compatible with automation. Our high detectability and statistical framework allowed precise linking of enhancers to promoters without merging single cells. We systematically investigated allele-specific accessibility in the mouse cerebral cortex, revealing allele-specific accessibility of promotors of certain imprinted genes but biallelic accessibility of their enhancers. Finally, we combined METATAC with our high-sensitivity single-cell RNA sequencing (scRNA-seq) method, multiple annealing and looping based amplification cycles for digital transcriptomics (MALBAC-DT), to develop a joint ATAC–RNA assay, termed METATAC and MALBAC-DT coassay by sequencing (M2C-seq). M2C-seq achieved significant improvements for both ATAC and RNA compared with previous methods, with consistent performance across cell lines and early mouse embryos.

Practical bioinformatics pipelines for single-cell RNA-seq data analysis

Single-cell RNA sequencing (scRNA-seq) is a revolutionary tool to explore cells. With an increasing number of scRNA-seq data analysis tools that have been developed, it is challenging for users to choose and compare their performance. Here, we present an overview of the workflow for computational analysis of scRNA-seq data. We detail the steps of a typical scRNA-seq analysis, including experimental design, pre-processing and quality control, feature selection, dimensionality reduction, cell clustering and annotation, and downstream analysis including batch correction, trajectory inference and cell-cell communication. We provide guidelines according to our best practice. This review will be helpful for the experimentalists interested in analyzing their data, and will aid the users seeking to update their analysis pipelines.

Single-cell-resolved dynamics of chromatin architecture delineate cell and regulatory states in zebrafish embryos

DNA accessibility of cis-regulatory elements (CREs) dictates transcriptional activity and drives cell differentiation during development. While many genes regulating embryonic development have been identified, the underlying CRE dynamics controlling their expression remain largely uncharacterized. To address this, we produced a multimodal resource and genomic regulatory map for the zebrafish community, which integrates single-cell combinatorial indexing assay for transposase-accessible chromatin with high-throughput sequencing (sci-ATAC-seq) with bulk histone PTMs and Hi-C data to achieve a genome-wide classification of the regulatory architecture determining transcriptional activity in the 24-h post-fertilization (hpf) embryo. We characterized the genome-wide chromatin architecture at bulk and single-cell resolution, applying sci-ATAC-seq on whole 24-hpf stage zebrafish embryos, generating accessibility profiles for ∼23,000 single nuclei. We developed a genome segmentation method, ScregSeg (single-cell regulatory landscape segmentation), for defining regulatory programs, and candidate CREs, specific to one or more cell types. We integrated the ScregSeg output with bulk measurements for histone post-translational modifications and 3D genome organization and identified new regulatory principles between chromatin modalities prevalent during zebrafish development. Sci-ATAC-seq profiling of npas4l/cloche mutant embryos identified novel cellular roles for this hematovascular transcriptional master regulator and suggests an intricate mechanism regulating its expression. Our work defines regulatory architecture and principles in the zebrafish embryo and establishes a resource of cell-type-specific genome-wide regulatory annotations and candidate CREs, providing a valuable open resource for genomics, developmental, molecular, and computational biology.

Single-nucleus RNA sequencing of plant tissues using a nanowell-based system

Single-cell genomics provides unprecedented potential for research on plant development and environmental responses. Here, we introduce a generic procedure for plant nucleus isolation combined with nanowell-based library preparation. Our method enables the transcriptome analysis of thousands of individual plant nuclei. It serves as an alternative to the use of protoplast isolation, which is currently a standard methodology for plant single-cell genomics, although it can be challenging for some plant tissues. We show the applicability of our nucleus isolation method by using different plant materials from different species. The potential of our single-nucleus RNA sequencing method is shown through the characterization of transcriptomes of seedlings and developing flowers from Arabidopsis thaliana. We evaluated the transcriptome dynamics during the early stages of anther development, identified stage-specific activities of transcription factors regulating this process, and predicted potential target genes of these transcription factors. Our nucleus isolation procedure can be applied in different plant species and tissues, thus expanding the toolkit for plant single-cell genomics experiments.

Single-Cell RNA Sequencing (scRNA-seq) in Cardiac Tissue: Applications and Limitations

Cardiovascular diseases (CVDs) are a group of disorders of the blood vessels and heart, which are considered as the leading causes of death worldwide. The pathology of CVDs could be related to the functional abnormalities of multiple cell types in the heart. Single-cell RNA sequencing (scRNA-seq) technology is a powerful method for characterizing individual cells and elucidating the molecular mechanisms by providing a high resolution of transcriptomic changes at the single-cell level. Specifically, scRNA-seq has provided novel insights into CVDs by identifying rare cardiac cell types, inferring the trajectory tree, estimating RNA velocity, elucidating the cell-cell communication, and comparing healthy and pathological heart samples. In this review, we summarize the different scRNA-seq platforms and published single-cell datasets in the cardiovascular field, and describe the utilities and limitations of this technology. Lastly, we discuss the future perspective of the application of scRNA-seq technology into cardiovascular research.

Scalable, multimodal profiling of chromatin accessibility, gene expression and protein levels in single cells

Recent technological advances have enabled massively parallel chromatin profiling with scATAC-seq (single-cell assay for transposase accessible chromatin by sequencing). Here we present ATAC with select antigen profiling by sequencing (ASAP-seq), a tool to simultaneously profile accessible chromatin and protein levels. Our approach pairs sparse scATAC-seq data with robust detection of hundreds of cell surface and intracellular protein markers and optional capture of mitochondrial DNA for clonal tracking, capturing three distinct modalities in single cells. ASAP-seq uses a bridging approach that repurposes antibody:oligonucleotide conjugates designed for existing technologies that pair protein measurements with single-cell RNA sequencing. Together with DOGMA-seq, an adaptation of CITE-seq (cellular indexing of transcriptomes and epitopes by sequencing) for measuring gene activity across the central dogma of gene regulation, we demonstrate the utility of systematic multi-omic profiling by revealing coordinated and distinct changes in chromatin, RNA and surface proteins during native hematopoietic differentiation and peripheral blood mononuclear cell stimulation and as a combinatorial decoder and reporter of multiplexed perturbations in primary T cells.

Leveraging Novel Integrated Single-Cell Analyses to Define HIV-1 Latency Reversal

While suppressive antiretroviral therapy can effectively limit HIV-1 replication and evolution, it leaves behind a residual pool of integrated viral genomes that persist in a state of reversible nonproductive infection, referred to as the HIV-1 reservoir. HIV-1 infection models were established to investigate HIV-1 latency and its reversal; recent work began to probe the dynamics of HIV-1 latency reversal at single-cell resolution. Signals that establish HIV-1 latency and govern its reactivation are complex and may not be completely resolved at the cellular and regulatory levels by the aggregated measurements of bulk cellular-sequencing methods. High-throughput single-cell technologies that characterize and quantify changes to the epigenome, transcriptome, and proteome continue to rapidly evolve. Combinations of single-cell techniques, in conjunction with novel computational approaches to analyze these data, were developed and provide an opportunity to improve the resolution of the heterogeneity that may exist in HIV-1 reactivation. In this review, we summarize the published single-cell HIV-1 transcriptomic work and explore how cutting-edge advances in single-cell techniques and integrative data-analysis tools may be leveraged to define the mechanisms that control the reversal of HIV-1 latency.

Massively parallel single-cell mitochondrial DNA genotyping and chromatin profiling

Natural mitochondrial DNA (mtDNA) mutations enable the inference of clonal relationships among cells. mtDNA can be profiled along with measures of cell state, but has not yet been combined with the massively parallel approaches needed to tackle the complexity of human tissue. Here, we introduce a high-throughput, droplet-based mitochondrial single-cell assay for transposase-accessible chromatin with sequencing (scATAC-seq), a method that combines high-confidence mtDNA mutation calling in thousands of single cells with their concomitant high-quality accessible chromatin profile. This enables the inference of mtDNA heteroplasmy, clonal relationships, cell state and accessible chromatin variation in individual cells. We reveal single-cell variation in heteroplasmy of a pathologic mtDNA variant, which we associate with intra-individual chromatin variability and clonal evolution. We clonally trace thousands of cells from cancers, linking epigenomic variability to subclonal evolution, and infer cellular dynamics of differentiating hematopoietic cells in vitro and in vivo. Taken together, our approach enables the study of cellular population dynamics and clonal properties in vivo.

Clinicomolecular Identification of Conserved and Individualized Features of Granulomatous Uveitis

Objective

To identify molecular features that distinguish individuals with shared clinical features of granulomatous uveitis.

Design

Cross-sectional, observational study.

Participants

Four eyes from patients with active granulomatous uveitis.

Methods

We performed single-cell RNA-sequencing with antigen-receptor sequence analysis to obtain an unbiased gene expression survey of ocular immune cells and identify clonally expanded lymphocytes.

Main Outcomes Measures

For each inflamed eye, we measured the proportion of distinct immune cell types, the amount of B or T cell clonal expansion, and the transcriptional profile of T and B cells.

Results

Each individual had robust clonal expansion arising from a single T or B cell lineage, suggesting distinct, antigen-driven pathogenic processes in each patient. This variability in clonal expansion was mirrored by individual variability in CD4 T cell populations, whereas ocular CD8 T cells and B cells were more transcriptionally similar between patients. Finally, ocular B cells displayed evidence of class-switching and plasmablast differentiation within the ocular microenvironment, providing additional support for antigen-driven immune responses in granulomatous uveitis.

Conclusions

Collectively, our study identified both conserved and individualized features of granulomatous uveitis, illuminating parallel pathophysiologic mechanisms, and suggesting that future personalized therapeutic approaches may be warranted.

Tutorial: guidelines for the computational analysis of single-cell RNA sequencing data

Single-cell RNA sequencing (scRNA-seq) is a popular and powerful technology that allows you to profile the whole transcriptome of a large number of individual cells. However, the analysis of the large volumes of data generated from these experiments requires specialized statistical and computational methods. Here we present an overview of the computational workflow involved in processing scRNA-seq data. We discuss some of the most common tasks and the tools available for addressing central biological questions. In this article and our companion website ( https://scrnaseq-course.cog.sanger.ac.uk/website/index.html ), we provide guidelines regarding best practices for performing computational analyses. This tutorial provides a hands-on guide for experimentalists interested in analyzing their data as well as an overview for bioinformaticians seeking to develop new computational methods.

Beyond bulk single-chain sequencing: Getting at the whole receptor

Recent advancements in paired B-cell receptor sequencing technologies have accelerated the development of simpler, high-throughput pipelines for generating native antibody heavy and light chain pairs used to elucidate novel antibodies and provide insights into antibody response against pathogenic targets. These technologies involve single-cell isolation, using either single wells or emulsified droplets to maintain physical separation of individual cells, followed by sequencing. The development of novel single wells and emulsion-based workflows addresses key challenges by improving throughput of single-cell analyses, reducing method complexity, and integrating functional assays into existing workflows. Enabled by paired B-cell receptor sequencing, functional characterization of pathogen-specific antibodies reveals immunological insights beyond bulk sequencing.

Single-cell atlas of colonic CD8+ T cells in ulcerative colitis

Colonic antigen-experienced lymphocytes such as tissue-resident memory CD8+ T cells can respond rapidly to repeated antigen exposure. However, their cellular phenotypes and the mechanisms by which they drive immune regulation and inflammation remain unclear. Here we compiled an unbiased atlas of human colonic CD8+ T cells in health and ulcerative colitis (UC) using single-cell transcriptomics with T-cell receptor repertoire analysis and mass cytometry. We reveal extensive heterogeneity in CD8+ T-cell composition, including expanded effector and post-effector terminally differentiated CD8+ T cells. While UC-associated CD8+ effector T cells can trigger tissue destruction and produce tumor necrosis factor (TNF)-α, post-effector cells acquire innate signatures to adopt regulatory functions that may mitigate excessive inflammation. Thus, we identify colonic CD8+ T-cell phenotypes in health and UC, define their clonal relationships and characterize terminally differentiated dysfunctional UC CD8+ T cells expressing IL-26, which attenuate acute colitis in a humanized IL-26 transgenic mouse model.