Benchmarking Computational Doublet-Detection Methods for Single-Cell RNA Sequencing Data

Single-cell transcriptomic analysis of endometriosis

Endometriosis is a common condition in women that causes chronic pain and infertility and is associated with an elevated risk of ovarian cancer. We profiled transcriptomes of >370,000 individual cells from endometriomas (n = 8), endometriosis (n = 28), eutopic endometrium (n = 10), unaffected ovary (n = 4) and endometriosis-free peritoneum (n = 4), generating a cellular atlas of endometrial-type epithelial cells, stromal cells and microenvironmental cell populations across tissue sites. Cellular and molecular signatures of endometrial-type epithelium and stroma differed across tissue types, suggesting a role for cellular restructuring and transcriptional reprogramming in the disease. Epithelium, stroma and proximal mesothelial cells of endometriomas showed dysregulation of pro-inflammatory pathways and upregulation of complement proteins. Somatic ARID1A mutation in epithelial cells was associated with upregulation of pro-angiogenic and pro-lymphangiogenic factors and remodeling of the endothelial cell compartment, with enrichment of lymphatic endothelial cells. Finally, signatures of ciliated epithelial cells were enriched in ovarian cancers, reinforcing epidemiologic associations between these two diseases.

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 S³ (S³: 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 S³, along with general pointers on how to subsequently analyze the resultant single-cell RNA-seq data.

Vaeda computationally annotates doublets in single-cell RNA sequencing data

MOTIVATION

Single-cell RNA sequencing (scRNA-seq) continues to expand our knowledge by facilitating the study of transcriptional heterogeneity at the level of single cells. Despite this technology's utility and success in biomedical research, technical artifacts are present in scRNA-seq data. Doublets/multiplets are a type of artifact that occurs when two or more cells are tagged by the same barcode, and therefore they appear as a single cell. Because this introduces non-existent transcriptional profiles, doublets can bias and mislead downstream analysis. To address this limitation, computational methods to annotate and remove doublets form scRNA-seq datasets are needed.

RESULTS

We introduce vaeda (Variational Auto-Encoder for Doublet Annotation), a new approach for computational annotation of doublets in scRNA-seq data. Vaeda integrates a variational auto-encoder and Positive-Unlabeled learning to produce doublet scores and binary doublet calls. We apply vaeda, along with seven existing doublet annotation methods, to 16 benchmark datasets and find that vaeda performs competitively in terms of doublet scores and doublet calls. Notably, vaeda outperforms other python-based methods for doublet annotation. Altogether, vaeda is a robust and competitive method for scRNA-seq doublet annotation and may be of particular interest in the context of python-based workflows.

AVAILABILITY AND IMPLEMENTATION

Vaeda is available at https://github.com/kostkalab/vaeda, and the version used for the results we present here is archived at zenodo (https://doi.org/10.5281/zenodo.7199783).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Cellsnake: a user-friendly tool for single-cell RNA sequencing analysis

BACKGROUND

Single-cell RNA sequencing (scRNA-seq) provides high-resolution transcriptome data to understand the heterogeneity of cell populations at the single-cell level. The analysis of scRNA-seq data requires the utilization of numerous computational tools. However, nonexpert users usually experience installation issues, a lack of critical functionality or batch analysis modes, and the steep learning curves of existing pipelines.

RESULTS

We have developed cellsnake, a comprehensive, reproducible, and accessible single-cell data analysis workflow, to overcome these problems. Cellsnake offers advanced features for standard users and facilitates downstream analyses in both R and Python environments. It is also designed for easy integration into existing workflows, allowing for rapid analyses of multiple samples.

CONCLUSION

As an open-source tool, cellsnake is accessible through Bioconda, PyPi, Docker, and GitHub, making it a cost-effective and user-friendly option for researchers. By using cellsnake, researchers can streamline the analysis of scRNA-seq data and gain insights into the complex biology of single cells.

High-resolution single-cell atlas reveals diversity and plasticity of tissue-resident neutrophils in non-small cell lung cancer

Non-small cell lung cancer (NSCLC) is characterized by molecular heterogeneity with diverse immune cell infiltration patterns, which has been linked to therapy sensitivity and resistance. However, full understanding of how immune cell phenotypes vary across different patient subgroups is lacking. Here, we dissect the NSCLC tumor microenvironment at high resolution by integrating 1,283,972 single cells from 556 samples and 318 patients across 29 datasets, including our dataset capturing cells with low mRNA content. We stratify patients into immune-deserted, B cell, T cell, and myeloid cell subtypes. Using bulk samples with genomic and clinical information, we identify cellular components associated with tumor histology and genotypes. We then focus on the analysis of tissue-resident neutrophils (TRNs) and uncover distinct subpopulations that acquire new functional properties in the tissue microenvironment, providing evidence for the plasticity of TRNs. Finally, we show that a TRN-derived gene signature is associated with anti-programmed cell death ligand 1 (PD-L1) treatment failure.

Data analysis guidelines for single-cell RNA-seq in biomedical studies and clinical applications

The application of single-cell RNA sequencing (scRNA-seq) in biomedical research has advanced our understanding of the pathogenesis of disease and provided valuable insights into new diagnostic and therapeutic strategies. With the expansion of capacity for high-throughput scRNA-seq, including clinical samples, the analysis of these huge volumes of data has become a daunting prospect for researchers entering this field. Here, we review the workflow for typical scRNA-seq data analysis, covering raw data processing and quality control, basic data analysis applicable for almost all scRNA-seq data sets, and advanced data analysis that should be tailored to specific scientific questions. While summarizing the current methods for each analysis step, we also provide an online repository of software and wrapped-up scripts to support the implementation. Recommendations and caveats are pointed out for some specific analysis tasks and approaches. We hope this resource will be helpful to researchers engaging with scRNA-seq, in particular for emerging clinical applications.

Clustering single-cell multi-omics data with MoClust

MOTIVATION

Single-cell multi-omics sequencing techniques have rapidly developed in the past few years. Clustering analysis with single-cell multi-omics data may give us novel perspectives to dissect cellular heterogeneity. However, multi-omics data have the properties of inherited large dimension, high sparsity and existence of doublets. Moreover, representations of different omics from even the same cell follow diverse distributions. Without proper distribution alignment techniques, clustering methods will encounter less separable clusters easily affected by less informative omics data.

RESULTS

We developed MoClust, a novel joint clustering framework that can be applied to several types of single-cell multi-omics data. A selective automatic doublet detection module that can identify and filter out doublets is introduced in the pretraining stage to improve data quality. Omics-specific autoencoders are introduced to characterize the multi-omics data. A contrastive learning way of distribution alignment is adopted to adaptively fuse omics representations into an omics-invariant representation. This novel way of alignment boosts the compactness and separableness of clusters, while accurately weighting the contribution of each omics to the clustering object. Extensive experiments, over both simulated and real multi-omics datasets, demonstrated the powerful alignment, doublet detection and clustering ability features of MoClust.

AVAILABILITY AND IMPLEMENTATION

An implementation of MoClust is available from https://doi.org/10.5281/zenodo.7306504.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Challenges in translational machine learning

Machine learning (ML) algorithms are increasingly being used to help implement clinical decision support systems. In this new field, we define as "translational machine learning", joint efforts and strong communication between data scientists and clinicians help to span the gap between ML and its adoption in the clinic. These collaborations also improve interpretability and trust in translational ML methods and ultimately aim to result in generalizable and reproducible models. To help clinicians and bioinformaticians refine their translational ML pipelines, we review the steps from model building to the use of ML in the clinic. We discuss experimental setup, computational analysis, interpretability and reproducibility, and emphasize the challenges involved. We highly advise collaboration and data sharing between consortia and institutes to build multi-centric cohorts that facilitate ML methodologies that generalize across centers. In the end, we hope that this review provides a way to streamline translational ML and helps to tackle the challenges that come with it.

Comparative single-cell transcriptomes of dose and time dependent epithelial–mesenchymal spectrums

Epithelial–mesenchymal transition (EMT) is a cellular process involved in development and disease progression. Intermediate EMT states were observed in tumors and fibrotic tissues, but previous in vitro studies focused on time-dependent responses with single doses of signals; it was unclear whether single-cell transcriptomes support stable intermediates observed in diseases. Here, we performed single-cell RNA-sequencing with human mammary epithelial cells treated with multiple doses of TGF-β. We found that dose-dependent EMT harbors multiple intermediate states at nearly steady state. Comparisons of dose- and time-dependent EMT transcriptomes revealed that the dose-dependent data enable higher sensitivity to detect genes associated with EMT. We identified cell clusters unique to time-dependent EMT, reflecting cells en route to stable states. Combining dose- and time-dependent cell clusters gave rise to accurate prognosis for cancer patients. Our transcriptomic data and analyses uncover a stable EMT continuum at the single-cell resolution, and complementary information of two types of single-cell experiments.

Improving the diagnosis of thyroid cancer by machine learning and clinical data

Thyroid cancer is a common endocrine carcinoma that occurs in the thyroid gland. Much effort has been invested in improving its diagnosis, and thyroidectomy remains the primary treatment method. A successful operation without unnecessary side injuries relies on an accurate preoperative diagnosis. Current human assessment of thyroid nodule malignancy is prone to errors and may not guarantee an accurate preoperative diagnosis. This study proposed a machine learning framework to predict thyroid nodule malignancy based on our collected novel clinical dataset. The ten-fold cross-validation, bootstrap analysis, and permutation predictor importance were applied to estimate and interpret the model performance under uncertainty. The comparison between model prediction and expert assessment shows the advantage of our framework over human judgment in predicting thyroid nodule malignancy. Our method is accurate, interpretable, and thus useable as additional evidence in the preoperative diagnosis of thyroid cancer.

SAREV: A review on statistical analytics of single-cell RNA sequencing data

Due to the development of next-generation RNA sequencing (NGS) technologies, there has been tremendous progress in research involving determining the role of genomics, transcriptomics and epigenomics in complex biological systems. However, scientists have realized that information obtained using earlier technology, frequently called 'bulk RNA-seq' data, provides information averaged across all the cells present in a tissue. Relatively newly developed single cell (scRNA-seq) technology allows us to provide transcriptomic information at a single-cell resolution. Nevertheless, these high-resolution data have their own complex natures and demand novel statistical data analysis methods to provide effective and highly accurate results on complex biological systems. In this review, we cover many such recently developed statistical methods for researchers wanting to pursue scRNA-seq statistical and computational research as well as scientific research about these existing methods and free software tools available for their generated data. This review is certainly not exhaustive due to page limitations. We have tried to cover the popular methods starting from quality control to the downstream analysis of finding differentially expressed genes and concluding with a brief description of network analysis.

Statistical learning in preclinical drug proarrhythmic assessment

Torsades de pointes (TdP) is an irregular heart rhythm characterized by faster beat rates and potentially could lead to sudden cardiac death. Much effort has been invested in understanding the drug-induced TdP in preclinical studies. However, a comprehensive statistical learning framework that can accurately predict the drug-induced TdP risk from preclinical data is still lacking. We proposed ordinal logistic regression and ordinal random forest models to predict low-, intermediate-, and high-risk drugs based on datasets generated from two experimental protocols. Leave-one-drug-out cross-validation, stratified bootstrap, and permutation predictor importance were applied to estimate and interpret the model performance under uncertainty. The potential outlier drugs identified by our models are consistent with their descriptions in the literature. Our method is accurate, interpretable, and thus useable as supplemental evidence in the drug safety assessment.

From single-omics to interactomics: How can ligand-induced perturbations modulate single-cell phenotypes?

Cells suffer from perturbations by different stimuli, which, consequently, rise to individual alterations in their profile and function that may end up affecting the tissue as a whole. This is no different if we consider the effect of a therapeutic agent on a biological system. As cells are exposed to external ligands their profile can change at different single-omics levels. Detecting how these changes take place through different sequencing technologies is key to a better understanding of the effects of therapeutic agents. Single-cell RNA-sequencing stands out as one of the most common approaches for cell profiling and perturbation analysis. As a result, single-cell transcriptomics data can be integrated with other omics data sources, such as proteomics and epigenomics data, to clarify the perturbation effects and mechanism at the cell level. Appropriate computational tools are key to process and integrate the available information. This chapter focuses on the recent advances on ligand-induced perturbation and single-cell omics computational tools and algorithms, their current limitations, and how the deluge of data can be used to improve the current process of drug research and development.

Chord: an ensemble machine learning algorithm to identify doublets in single-cell RNA sequencing data

High-throughput single-cell RNA sequencing (scRNA-seq) is a popular method, but it is accompanied by doublet rate problems that disturb the downstream analysis. Several computational approaches have been developed to detect doublets. However, most of these methods may yield satisfactory performance in some datasets but lack stability in others; thus, it is difficult to regard a single method as the gold standard which can be applied to all types of scenarios. It is a difficult and time-consuming task for researchers to choose the most appropriate software. We here propose Chord which implements a machine learning algorithm that integrates multiple doublet detection methods to address these issues. Chord had higher accuracy and stability than the individual approaches on different datasets containing real and synthetic data. Moreover, Chord was designed with a modular architecture port, which has high flexibility and adaptability to the incorporation of any new tools. Chord is a general solution to the doublet detection problem.

For the unmet need to choose the suitable doublet detection method, an ensemble machine learning algorithm called Chord was developed, which integrates multiple methods and achieves higher accuracy and stability on different scRNA-seq datasets.

Doublet identification in single-cell sequencing data using scDblFinder

Doublets are prevalent in single-cell sequencing data and can lead to artifactual findings. A number of strategies have therefore been proposed to detect them. Building on the strengths of existing

approaches, we developed scDblFinder, a fast, flexible and accurate Bioconductor-based doublet detection method. Here we present the method, justify its design choices, demonstrate its performance on both single-cell RNA and accessibility (ATAC) sequencing data, and provide some observations on doublet formation, detection, and enrichment analysis. Even in complex datasets, scDblFinder can accurately identify most heterotypic doublets, and was already found by an independent benchmark to outcompete alternatives.

Reconstructing physical cell interaction networks from single-cell data using Neighbor-seq

Cell-cell interactions are the fundamental building blocks of tissue organization and multicellular life. We developed Neighbor-seq, a method to identify and annotate the architecture of direct cell–cell interactions and relevant ligand–receptor signaling from the undissociated cell fractions in massively parallel single cell sequencing data. Neighbor-seq accurately identifies microanatomical features of diverse tissue types such as the small intestinal epithelium, terminal respiratory tract, and splenic white pulp. It also captures the differing topologies of cancer-immune-stromal cell communications in pancreatic and skin tumors, which are consistent with the patterns observed in spatial transcriptomic data. Neighbor-seq is fast and scalable. It draws inferences from routine single-cell data and does not require prior knowledge about sample cell-types or multiplets. Neighbor-seq provides a framework to study the organ-level cellular interactome in health and disease, bridging the gap between single-cell and spatial transcriptomics.

Statistical and machine learning methods for spatially resolved transcriptomics data analysis

The recent advancement in spatial transcriptomics technology has enabled multiplexed profiling of cellular transcriptomes and spatial locations. As the capacity and efficiency of the experimental technologies continue to improve, there is an emerging need for the development of analytical approaches. Furthermore, with the continuous evolution of sequencing protocols, the underlying assumptions of current analytical methods need to be re-evaluated and adjusted to harness the increasing data complexity. To motivate and aid future model development, we herein review the recent development of statistical and machine learning methods in spatial transcriptomics, summarize useful resources, and highlight the challenges and opportunities ahead.

Single-Cell Antibody Sequencing in Atherosclerosis Research

The transcriptomic information obtained by single cell RNA sequencing (scRNA-seq) can be supplemented by information on the cell surface phenotype by using oligonucleotide-tagged monoclonal antibodies (scAb-Seq). This is of particular importance in immune cells, where the correlation between mRNA and cell surface expression is very weak. scAb-Seq is facilitated by the availability of commercial antibodies and antibody mixes. Now panels of up to 200 antibodies are available for human and mouse cells. Proteins are detected by antibodies conjugated to a tripartite DNA sequence that contains a primer for amplification and sequencing, a unique oligonucleotide that acts as an antibody barcode and a poly(dA) sequence, simultaneously detecting extension of antibody-specific DNA sequences and cDNAs in the same poly(dT)-primed reaction. For each cell, surface protein expression is captured and sequenced along with the cell's transcriptome. Here, we list the steps needed to produce antibody sequencing data from tissue or blood cells.

Guidelines for bioinformatics of single-cell sequencing data analysis in Alzheimer's disease: review, recommendation, implementation and application

Alzheimer's disease (AD) is the most common form of dementia, characterized by progressive cognitive impairment and neurodegeneration. Extensive clinical and genomic studies have revealed biomarkers, risk factors, pathways, and targets of AD in the past decade. However, the exact molecular basis of AD development and progression remains elusive. The emerging single-cell sequencing technology can potentially provide cell-level insights into the disease. Here we systematically review the state-of-the-art bioinformatics approaches to analyze single-cell sequencing data and their applications to AD in 14 major directions, including 1) quality control and normalization, 2) dimension reduction and feature extraction, 3) cell clustering analysis, 4) cell type inference and annotation, 5) differential expression, 6) trajectory inference, 7) copy number variation analysis, 8) integration of single-cell multi-omics, 9) epigenomic analysis, 10) gene network inference, 11) prioritization of cell subpopulations, 12) integrative analysis of human and mouse sc-RNA-seq data, 13) spatial transcriptomics, and 14) comparison of single cell AD mouse model studies and single cell human AD studies. We also address challenges in using human postmortem and mouse tissues and outline future developments in single cell sequencing data analysis. Importantly, we have implemented our recommended workflow for each major analytic direction and applied them to a large single nucleus RNA-sequencing (snRNA-seq) dataset in AD. Key analytic results are reported while the scripts and the data are shared with the research community through  GitHub. In summary, this comprehensive review provides insights into various approaches to analyze single cell sequencing data and offers specific guidelines for study design and a variety of analytic directions. The review and the accompanied software tools will serve as a valuable resource for studying cellular and molecular mechanisms of AD, other diseases, or biological systems at the single cell level.

Cellular and molecular atlas of the placenta from a COVID-19 pregnant woman infected at midgestation highlights the defective impacts on foetal health

Objectives

The impacts of the current COVID‐19 pandemic on maternal and foetal health are enormous and of serious concern. However, the influence of SARS‐CoV‐2 infection at early‐to‐mid gestation on maternal and foetal health remains unclear.

Materials and methods

Here, we report the follow‐up study of a pregnant woman of her whole infective course of SARS‐CoV‐2, from asymptomatic infection at gestational week 20 to mild and then severe illness state, and finally cured at Week 24. Following caesarean section due to incomplete uterine rupture at Week 28, histological examinations on the placenta and foetal tissues as well as single‐cell RNA sequencing (scRNA‐seq) for the placenta were performed.

Results

Compared with the gestational age‐matched control placentas, the placenta from this COVID‐19 case exhibited more syncytial knots and lowered expression of syncytiotrophoblast‐related genes. The scRNA‐seq analysis demonstrated impaired trophoblast differentiation, activation of antiviral and inflammatory CD8 T cells, as well as the tight association of increased inflammatory responses in the placenta with complement over‐activation in macrophages. In addition, levels of several inflammatory factors increased in the placenta and foetal blood.

Conclusion

These findings illustrate a systematic cellular and molecular signature of placental insufficiency and immune activation at the maternal–foetal interface that may be attributed to SARS‐CoV‐2 infection at the midgestation stage, which highly suggests the extensive care for maternal and foetal outcomes in pregnant women suffering from COVID‐19.

Epithelial GPR35 protects from Citrobacter rodentium infection by preserving goblet cells and mucosal barrier integrity

Goblet cells secrete mucin to create a protective mucus layer against invasive bacterial infection and are therefore essential for maintaining intestinal health. However, the molecular pathways that regulate goblet cell function remain largely unknown. Although GPR35 is highly expressed in colonic epithelial cells, its importance in promoting the epithelial barrier is unclear. In this study, we show that epithelial Gpr35 plays a critical role in goblet cell function. In mice, cell-type-specific deletion of Gpr35 in epithelial cells but not in macrophages results in goblet cell depletion and dysbiosis, rendering these animals more susceptible to Citrobacter rodentium infection. Mechanistically, scRNA-seq analysis indicates that signaling of epithelial Gpr35 is essential to maintain normal pyroptosis levels in goblet cells. Our work shows that the epithelial presence of Gpr35 is a critical element for the function of goblet cell-mediated symbiosis between host and microbiota.

Multi-omics integration in the age of million single-cell data

An explosion in single-cell technologies has revealed a previously underappreciated heterogeneity of cell types and novel cell-state associations with sex, disease, development and other processes. Starting with transcriptome analyses, single-cell techniques have extended to multi-omics approaches and now enable the simultaneous measurement of data modalities and spatial cellular context. Data are now available for millions of cells, for whole-genome measurements and for multiple modalities. Although analyses of such multimodal datasets have the potential to provide new insights into biological processes that cannot be inferred with a single mode of assay, the integration of very large, complex, multimodal data into biological models and mechanisms represents a considerable challenge. An understanding of the principles of data integration and visualization methods is required to determine what methods are best applied to a particular single-cell dataset. Each class of method has advantages and pitfalls in terms of its ability to achieve various biological goals, including cell-type classification, regulatory network modelling and biological process inference. In choosing a data integration strategy, consideration must be given to whether the multi-omics data are matched (that is, measured on the same cell) or unmatched (that is, measured on different cells) and, more importantly, the overall modelling and visualization goals of the integrated analysis.

huARdb: human Antigen Receptor database for interactive clonotype-transcriptome analysis at the single-cell level

T-cell receptors (TCRs) and B-cell receptors (BCRs) are critical in recognizing antigens and activating the adaptive immune response. Stochastic V(D)J recombination generates massive TCR/BCR repertoire diversity. Single-cell immune profiling with transcriptome analysis allows the high-throughput study of individual TCR/BCR clonotypes and functions under both normal and pathological settings. However, a comprehensive database linking these data is not yet readily available. Here, we present the human Antigen Receptor database (huARdb), a large-scale human single-cell immune profiling database that contains 444 794 high confidence T or B cells (hcT/B cells) with full-length TCR/BCR sequence and transcriptomes from 215 datasets. All datasets were processed in a uniform workflow, including sequence alignment, cell subtype prediction, unsupervised cell clustering, and clonotype definition. We also developed a multi-functional and user-friendly web interface that provides interactive visualization modules for biologists to analyze the transcriptome and TCR/BCR features at the single-cell level. HuARdb is freely available at https://huarc.net/database with functions for data querying, browsing, downloading, and depositing. In conclusion, huARdb is a comprehensive and multi-perspective atlas for human antigen receptors.

Single-Cell RNA-Seq of T Cells in B-ALL Patients Reveals an Exhausted Subset with Remarkable Heterogeneity

Characterization of functional T cell clusters is key to developing strategies for immunotherapy and predicting clinical responses in leukemia. Here, single-cell RNA sequencing is performed with T cells sorted from the peripheral blood of healthy individuals and patients with B cell-acute lymphoblastic leukemia (B-ALL). Unbiased bioinformatics analysis enabled the authors to identify 13 T cell clusters in the patients based on their molecular properties. All 11 major T cell subsets in healthy individuals are found in the patients with B-ALL, with the counterparts in the patients universally showing more activated characteristics. Two exhausted T cell populations, characterized by up-regulation of TIGIT, PDCD1, HLADRA, LAG3, and CTLA4 are specifically discovered in B-ALL patients. Of note, these exhausted T cells possess remarkable heterogeneity, and ten sub-clusters are further identified, which are characterized by different cell cycle phases, naïve states, and GNLY (coding granulysin) expression. Coupled with single-cell T cell receptor repertoire profiling, diverse originations of the exhausted T cells in B-ALL are suggested, and clonally expanded exhausted T cells are likely to originate from CD8+ effector memory/terminal effector cells. Together, these data provide for the first-time valuable insights for understanding exhausted T cell populations in leukemia.

Doublet identification in single-cell sequencing data using scDblFinder

Doublets are prevalent in single-cell sequencing data and can lead to artifactual findings. A number of strategies have therefore been proposed to detect them. Building on the strengths of existing approaches, we developed scDblFinder, a fast, flexible and accurate Bioconductor-based doublet detection method. Here we present the method, justify its design choices, demonstrate its performance on both single-cell RNA and accessibility (ATAC) sequencing data, and provide some observations on doublet formation, detection, and enrichment analysis. Even in complex datasets, scDblFinder can accurately identify most heterotypic doublets, and was already found by an independent benchmark to outcompete alternatives.

Genetic demultiplexing of pooled single-cell RNA-sequencing samples in cancer facilitates effective experimental design

BACKGROUND

Pooling cells from multiple biological samples prior to library preparation within the same single-cell RNA sequencing experiment provides several advantages, including lower library preparation costs and reduced unwanted technological variation, such as batch effects. Computational demultiplexing tools based on natural genetic variation between individuals provide a simple approach to demultiplex samples, which does not require complex additional experimental procedures. However, to our knowledge these tools have not been evaluated in cancer, where somatic variants, which could differ between cells from the same sample, may obscure the signal in natural genetic variation.

RESULTS

Here, we performed in silico benchmark evaluations by combining raw sequencing reads from multiple single-cell samples in high-grade serous ovarian cancer, which has a high copy number burden, and lung adenocarcinoma, which has a high tumor mutational burden. Our results confirm that genetic demultiplexing tools can be effectively deployed on cancer tissue using a pooled experimental design, although high proportions of ambient RNA from cell debris reduce performance.

CONCLUSIONS

This strategy provides significant cost savings through pooled library preparation. To facilitate similar analyses at the experimental design phase, we provide freely accessible code and a reproducible Snakemake workflow built around the best-performing tools found in our in silico benchmark evaluations, available at https://github.com/lmweber/snp-dmx-cancer.

AMULET: a novel read count-based method for effective multiplet detection from single nucleus ATAC-seq data

Detecting multiplets in single nucleus (sn)ATAC-seq data is challenging due to data sparsity and limited dynamic range. AMULET (ATAC-seq MULtiplet Estimation Tool) enumerates regions with greater than two uniquely aligned reads across the genome to effectively detect multiplets. We evaluate the method by generating snATAC-seq data in the human blood and pancreatic islet samples. AMULET has high precision, estimated via donor-based multiplexing, and high recall, estimated via simulated multiplets, compared to alternatives and identifies multiplets most effectively when a certain read depth of 25K median valid reads per nucleus is achieved.

Heterogeneity of immune cells in human atherosclerosis revealed by scRNA-Seq

Immune cells in atherosclerosis include T, B, natural killer (NK) and NKT cells, macrophages, monocytes, dendritic cells (DCs), neutrophils and mast cells. Advances in single cell RNA sequencing (sRNA-Seq) have refined our understanding of immune cell subsets. Four recent studies have used scRNA-Seq of immune cells in human atherosclerotic lesions and peripheral blood mononuclear cells (PBMCs), some including cell surface phenotypes revealed by oligonucleotide-tagged antibodies, which confirmed known and identified new immune cell subsets and identified genes significantly upregulated in PBMCs from HIV+ subjects with atherosclerosis compared to PBMCs from matched HIV+ subjects without atherosclerosis. The ability of scRNA-Seq to identify cell types is greatly augmented by adding cell surface phenotype using antibody sequencing. In this review we summarize the latest data obtained by scRNA-Seq on plaques and human PBMCs in human subjects with atherosclerosis.

Protocol for executing and benchmarking eight computational doublet-detection methods in single-cell RNA sequencing data analysis

The existence of doublets is a key confounder in single-cell RNA sequencing (scRNA-seq) data analysis. Computational techniques have been developed for detecting doublets from scRNA-seq data. We developed an R package DoubletCollection to integrate the installation and execution of eight doublet detection methods. DoubletCollection provides a unified interface to perform and visualize downstream analysis after doublet detection. Here, we present a protocol of using DoubletCollection to benchmark doublet-detection methods. This protocol can accommodate new doublet-detection methods in the fast-growing scRNA-seq field.

For details on the use and execution of this protocol, please refer to Xi and Li (2020).

•

Integrate the installation and execution of eight doublet-detection methods

•

An interface to perform and visualize downstream analysis after doublet detection

•

A collection of real and synthetic scRNA-seq datasets with doublet annotations

•

A user-friendly R package to perform doublet detection on scRNA-seq count matrices

doubletD: detecting doublets in single-cell DNA sequencing data

Motivation

While single-cell DNA sequencing (scDNA-seq) has enabled the study of intratumor heterogeneity at an unprecedented resolution, current technologies are error-prone and often result in doublets where two or more cells are mistaken for a single cell. Not only do doublets confound downstream analyses, but the increase in doublet rate is also a major bottleneck preventing higher throughput with current single-cell technologies. Although doublet detection and removal are standard practice in scRNA-seq data analysis, options for scDNA-seq data are limited. Current methods attempt to detect doublets while also performing complex downstream analyses tasks, leading to decreased efficiency and/or performance.

Results

We present doubletD, the first standalone method for detecting doublets in scDNA-seq data. Underlying our method is a simple maximum likelihood approach with a closed-form solution. We demonstrate the performance of doubletD on simulated data as well as real datasets, outperforming current methods for downstream analysis of scDNA-seq data that jointly infer doublets as well as standalone approaches for doublet detection in scRNA-seq data. Incorporating doubletD in scDNA-seq analysis pipelines will reduce complexity and lead to more accurate results.

Availability and implementation

https://github.com/elkebir-group/doubletD.

Supplementary information

Supplementary data are available at Bioinformatics online.

A yeast-optimized single-cell transcriptomics platform elucidates how mycophenolic acid and guanine alter global mRNA levels

Stochastic gene expression leads to inherent variability in expression outcomes even in isogenic single-celled organisms grown in the same environment. The Drop-Seq technology facilitates transcriptomic studies of individual mammalian cells, and it has had transformative effects on the characterization of cell identity and function based on single-cell transcript counts. However, application of this technology to organisms with different cell size and morphology characteristics has been challenging. Here we present yeastDrop-Seq, a yeast-optimized platform for quantifying the number of distinct mRNA molecules in a cell-specific manner in individual yeast cells. Using yeastDrop-Seq, we measured the transcriptomic impact of the lifespan-extending compound mycophenolic acid and its epistatic agent guanine. Each treatment condition had a distinct transcriptomic footprint on isogenic yeast cells as indicated by distinct clustering with clear separations among the different groups. The yeastDrop-Seq platform facilitates transcriptomic profiling of yeast cells for basic science and biotechnology applications.

Plasticity and heterogeneity of thermogenic adipose tissue

The perception of adipose tissue, both in the scientific community and in the general population, has changed dramatically in the past 20 years. While adipose tissue was thought for a long time to be a rather simple lipid storage entity, it is now recognized as a highly heterogeneous organ and a critical regulator of systemic metabolism, composed of many different subtypes of cells, with important endocrine functions. Additionally, adipose tissue is nowadays recognized to contribute to energy turnover, due to the presence of specialized thermogenic adipocytes, which can be found in many adipose depots. This review discusses the unprecedented insights that we have gained into the heterogeneity of thermogenic adipocytes and their respective precursors due to the technical developments in single-cell and nucleus technologies. These methodological advances have increased our understanding of how adipose tissue catabolic function is influenced by developmental and intercellular communication events.

scDesign2: a transparent simulator that generates high-fidelity single-cell gene expression count data with gene correlations captured

A pressing challenge in single-cell transcriptomics is to benchmark experimental protocols and computational methods. A solution is to use computational simulators, but existing simulators cannot simultaneously achieve three goals: preserving genes, capturing gene correlations, and generating any number of cells with varying sequencing depths. To fill this gap, we propose scDesign2, a transparent simulator that achieves all three goals and generates high-fidelity synthetic data for multiple single-cell gene expression count-based technologies. In particular, scDesign2 is advantageous in its transparent use of probabilistic models and its ability to capture gene correlations via copulas.

Anti-bias training for (sc)RNA-seq: experimental and computational approaches to improve precision

RNA-seq, including single cell RNA-seq (scRNA-seq), is plagued by insufficient sensitivity and lack of precision. As a result, the full potential of (sc)RNA-seq is limited. Major factors in this respect are the presence of global bias in most datasets, which affects detection and quantitation of RNA in a length-dependent fashion. In particular, scRNA-seq is affected by technical noise and a high rate of dropouts, where the vast majority of original transcripts is not converted into sequencing reads. We discuss these biases origins and implications, bioinformatics approaches to correct for them, and how biases can be exploited to infer characteristics of the sample preparation process, which in turn can be used to improve library preparation.