Alevin efficiently estimates accurate gene abundances from dscRNA-seq data

recount3: summaries and queries for large-scale RNA-seq expression and splicing

We present recount3, a resource consisting of over 750,000 publicly available human and mouse RNA sequencing (RNA-seq) samples uniformly processed by our new Monorail analysis pipeline. To facilitate access to the data, we provide the recount3 and snapcount R/Bioconductor packages as well as complementary web resources. Using these tools, data can be downloaded as study-level summaries or queried for specific exon-exon junctions, genes, samples, or other features. Monorail can be used to process local and/or private data, allowing results to be directly compared to any study in recount3. Taken together, our tools help biologists maximize the utility of publicly available RNA-seq data, especially to improve their understanding of newly collected data. recount3 is available from http://rna.recount.bio .

Updates on Immunotherapy and Immune Landscape in Renal Clear Cell Carcinoma

Several clinicopathological features of clear cell renal cell carcinomas (ccRCC) contribute to make an "atypical" cancer, including resistance to chemotherapy, sensitivity to anti-angiogenesis therapy and ICIs despite a low mutational burden, and CD8+ T cell infiltration being the predictor for poor prognosis-normally CD8+ T cell infiltration is a good prognostic factor in cancer patients. These "atypical" features have brought researchers to investigate the molecular and immunological mechanisms that lead to the increased T cell infiltrates despite relatively low molecular burdens, as well as to decipher the immune landscape that leads to better response to ICIs. In the present study, we summarize the past and ongoing pivotal clinical trials of immunotherapies for ccRCC, emphasizing the potential molecular and cellular mechanisms that lead to the success or failure of ICI therapy. Single-cell analysis of ccRCC has provided a more thorough and detailed understanding of the tumor immune microenvironment and has facilitated the discovery of molecular biomarkers from the tumor-infiltrating immune cells. We herein will focus on the discussion of some major immune cells, including T cells and tumor-associated macrophages (TAM) in ccRCC. We will further provide some perspectives of using molecular and cellular biomarkers derived from these immune cell types to potentially improve the response rate to ICIs in ccRCC patients.

Building the mega single-cell transcriptome ocular meta-atlas

BACKGROUND

The development of highly scalable single-cell transcriptome technology has resulted in the creation of thousands of datasets, >30 in the retina alone. Analyzing the transcriptomes between different projects is highly desirable because this would allow for better assessment of which biological effects are consistent across independent studies. However it is difficult to compare and contrast data across different projects because there are substantial batch effects from computational processing, single-cell technology utilized, and the natural biological variation. While many single-cell transcriptome-specific batch correction methods purport to remove the technical noise, it is difficult to ascertain which method functions best.

RESULTS

We developed a lightweight R package (scPOP, single-cell Pick Optimal Parameters) that brings in batch integration methods and uses a simple heuristic to balance batch merging and cell type/cluster purity. We use this package along with a Snakefile-based workflow system to demonstrate how to optimally merge 766,615 cells from 33 retina datsets and 3 species to create a massive ocular single-cell transcriptome meta-atlas.

CONCLUSIONS

This provides a model for how to efficiently create meta-atlases for tissues and cells of interest.

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.

Single-cell RNA-seq unravels alterations of the human spermatogonial stem cell compartment in patients with impaired spermatogenesis

Despite the high incidence of male infertility, only 30% of infertile men receive a causative diagnosis. To explore the regulatory mechanisms governing human germ cell function in normal and impaired spermatogenesis (crypto), we performed single-cell RNA sequencing (>30,000 cells). We find major alterations in the crypto spermatogonial compartment with increased numbers of the most undifferentiated spermatogonia (PIWIL4⁺). We also observe a transcriptional switch within the spermatogonial compartment driven by increased and prolonged expression of the transcription factor EGR4. Intriguingly, the EGR4-regulated chromatin-associated transcriptional repressor UTF1 is downregulated at transcriptional and protein levels. This is associated with changes in spermatogonial chromatin structure and fewer A_dark spermatogonia, characterized by tightly compacted chromatin and serving as reserve stem cells. These findings suggest that crypto patients are disadvantaged, as fewer cells safeguard their germline’s genetic integrity. These identified spermatogonial regulators will be highly interesting targets to uncover genetic causes of male infertility.

Crypto(zoospermic) men show increased number of PIWIL4⁺/EGR4⁺ spermatogonia

Crypto undifferentiated spermatogonia over-activate the EGR4 regulatory network

The predicted EGR4 target UTF1 is downregulated in crypto spermatogonia

Crypto testes show reduced numbers of UTF1⁺ A_dark reserve spermatogonia

SC-JNMF: single-cell clustering integrating multiple quantification methods based on joint non-negative matrix factorization

Single-cell RNA-sequencing is a rapidly evolving technology that enables us to understand biological processes at unprecedented resolution. Single-cell expression analysis requires a complex data processing pipeline, and the pipeline is divided into two main parts: The quantification part, which converts the sequence information into gene-cell matrix data; the analysis part, which analyzes the matrix data using statistics and/or machine learning techniques. In the analysis part, unsupervised cell clustering plays an important role in identifying cell types and discovering cell diversity and subpopulations. Identified cell clusters are also used for subsequent analysis, such as finding differentially expressed genes and inferring cell trajectories. However, single-cell clustering using gene expression profiles shows different results depending on the quantification methods. Clustering results are greatly affected by the quantification method used in the upstream process. In other words, even if the original RNA-sequence data is the same, gene expression profiles processed by different quantification methods will produce different clusters. In this article, we propose a robust and highly accurate clustering method based on joint non-negative matrix factorization (joint-NMF) by utilizing the information from multiple gene expression profiles quantified using different methods from the same RNA-sequence data. Our joint-NMF can extract common factors among multiple gene expression profiles by applying each NMF under the constraint that one of the factorized matrices is shared among multiple NMFs. The joint-NMF determines more robust and accurate cell clustering results by leveraging multiple quantification methods compared to conventional clustering methods, which use only a single gene expression profile. Additionally, we showed the usefulness of discovering marker genes with the extracted features using our method.

Choice of pre-processing pipeline influences clustering quality of scRNA-seq datasets

BACKGROUND

Single-cell RNA sequencing (scRNA-seq) has quickly become one of the most dominant techniques in modern transcriptome assessment. In particular, 10X Genomics' Chromium system, with its high throughput approach, turn key and thorough user guide made this cutting-edge technique accessible to many laboratories using diverse animal models. However, standard pre-processing, including the alignment and cell filtering pipelines might not be ideal for every organism or tissue. Here we applied an alternative strategy, based on the pseudoaligner kallisto, on twenty-two publicly available single cell sequencing datasets from a wide range of tissues of eight organisms and compared the results with the standard 10X Genomics' Cell Ranger pipeline.

RESULTS

In most of the tested samples, kallisto produced higher sequencing read alignment rates and total gene detection rates in comparison to Cell Ranger. Although datasets processed with Cell Ranger had higher cell counts, outside of human and mouse datasets, these additional cells were routinely of low quality, containing low gene detection rates. Thorough downstream analysis of one kallisto processed dataset, obtained from the zebrafish pineal gland, revealed clearer clustering, allowing the identification of an additional photoreceptor cell type that previously went undetected. The finding of the new cluster suggests that the photoreceptive pineal gland is essentially a bi-chromatic tissue containing both green and red cone-like photoreceptors and implies that the alignment and pre-processing pipeline can affect the discovery of biologically-relevant cell types.

CONCLUSION

While Cell Ranger favors higher cell numbers, using kallisto results in datasets with higher median gene detection per cell. We could demonstrate that cell type identification was not hampered by the lower cell count, but in fact improved as a result of the high gene detection rate and the more stringent filtering. Depending on the acquired dataset, it can be beneficial to favor high quality cells and accept a lower cell count, leading to an improved classification of cell types.

Differential transcript usage analysis of bulk and single-cell RNA-seq data with DTUrtle

Motivation

Each year, the number of published bulk and single-cell RNA-seq datasets is growing exponentially. Studies analyzing such data are commonly looking at gene-level differences, while the collected RNA-seq data inherently represents reads of transcript isoform sequences. Utilizing transcriptomic quantifiers, RNA-seq reads can be attributed to specific isoforms, allowing for analysis of transcript-level differences. A differential transcript usage (DTU) analysis is testing for proportional differences in a gene’s transcript composition, and has been of rising interest for many research questions, such as analysis of differential splicing or cell-type identification.

Results

We present the R package DTUrtle, the first DTU analysis workflow for both bulk and single-cell RNA-seq datasets, and the first package to conduct a ‘classical’ DTU analysis in a single-cell context. DTUrtle extends established statistical frameworks, offers various result aggregation and visualization options and a novel detection probability score for tagged-end data. It has been successfully applied to bulk and single-cell RNA-seq data of human and mouse, confirming and extending key results. In addition, we present novel potential DTU applications like the identification of cell-type specific transcript isoforms as biomarkers.

Availability and implementation

The R package DTUrtle is available at https://github.com/TobiTekath/DTUrtle with extensive vignettes and documentation at https://tobitekath.github.io/DTUrtle/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Single-Cell RNAseq Analysis of lncRNAs

Mammalian genomes are pervasively transcribed and a small fraction of RNAs produced codify for proteins. The importance of noncoding RNAs for the maintenance of cell functions is well known (e.g., rRNAs, tRNAs), but only recently it was first demonstrated the involvement of microRNAs (miRNAs) in posttranscriptional regulation and then the activity of long noncoding RNAs (lncRNAs) in the regulation of miRNAs, DNA structure and protein function. LncRNAs have an expression more cell specific than other RNAs and basing on their subcellular localization exert different functions. In this book chapter we consider different protocols to evaluate the expression of lncRNAs at the single cell level using genome-wide approaches. We considered the skeletal muscle as example because the most abundant tissue in mammals involved in the regulation of metabolism and body movement. We firstly described how to isolate the smallest complete contractile system responsible for muscle metabolic and contractile traits (myofibers). We considered how to separate long and short RNAs to allow the sequencing of the full-length transcript using the SMART technique for the retrotranscription. Because of myofibers are multinucleated cells and because of it is better to perform single cell sequencing on fresh tissues we described the single-nucleus sequencing that can be applied to frozen tissues. The chapter concludes with a description of bioinformatics approaches to evaluate differential expression from single-cell or single-nucleus RNA sequencing.

The MURAL collection of prostate cancer patient-derived xenografts enables discovery through preclinical models of uro-oncology

Preclinical testing is a crucial step in evaluating cancer therapeutics. We aimed to establish a significant resource of patient-derived xenografts (PDXs) of prostate cancer for rapid and systematic evaluation of candidate therapies. The PDX collection comprises 59 tumors collected from 30 patients between 2012-2020, coinciding with availability of abiraterone and enzalutamide. The PDXs represent the clinico-pathological and genomic spectrum of prostate cancer, from treatment-naïve primary tumors to castration-resistant metastases. Inter- and intra-tumor heterogeneity in adenocarcinoma and neuroendocrine phenotypes is evident from bulk and single-cell RNA sequencing data. Organoids can be cultured from PDXs, providing further capabilities for preclinical studies. Using a 1 x 1 x 1 design, we rapidly identify tumors with exceptional responses to combination treatments. To govern the distribution of PDXs, we formed the Melbourne Urological Research Alliance (MURAL). This PDX collection is a substantial resource, expanding the capacity to test and prioritize effective treatments for prospective clinical trials in prostate cancer.

Genome-wide detection and classification of terpene synthase genes in Aquilaria agallochum

Agarwood, one of the precious woods in the globe, is produced by Aquilaria plant species during an upshot of wounding and infection. Produced as a defence response, the dark, fragrant resin gets secreted in the plant's duramen, which is impregnated with fragrant molecules with the due course. Agarwood has gained worldwide popularity due to its high aromatic oil, fragrance, and pharmaceutical value, which makes it highly solicited by numerous industries. Predominant chemical constituents of agarwood, sesquiterpenoids, and 2-(2-phenylethyl) chromones have been scrutinized to comprehend the scientific nature of the fragrant wood and develop novel products. However, the genes involved in the biosynthesis of these aromatic compounds are still not comprehensively studied in Aquilaria. In this study, publicly available genomic and transcriptomics data of Aquilaria agallochum were integrated to identify putative functional terpene synthase genes (TPSs). The in silico study enabled us to identify ninety-six TPSs, of which thirty-nine full-length genes were systematically classified into TPS-a, TPS-b, TPS-c, TPS-e, TPS-f, and TPS-g subfamilies based on their gene structure, conserve motif, and phylogenetic comparison with TPSs from other plant species. Analysis of the cis-regulatory elements present upstream of AaTPSs revealed their association with hormone, stress and light responses. In silico expression studies detected their up-regulation in stress induced tissue. This study provides a basic understanding of terpene synthase gene repertoire in Aquilaria agallochum and unlatches opportunities for the biochemical characterization and biotechnological exploration of these genes.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01040-z.

miQC: An adaptive probabilistic framework for quality control of single-cell RNA-sequencing data

Single-cell RNA-sequencing (scRNA-seq) has made it possible to profile gene expression in tissues at high resolution. An important preprocessing step prior to performing downstream analyses is to identify and remove cells with poor or degraded sample quality using quality control (QC) metrics. Two widely used QC metrics to identify a 'low-quality' cell are (i) if the cell includes a high proportion of reads that map to mitochondrial DNA (mtDNA) encoded genes and (ii) if a small number of genes are detected. Current best practices use these QC metrics independently with either arbitrary, uniform thresholds (e.g. 5%) or biological context-dependent (e.g. species) thresholds, and fail to jointly model these metrics in a data-driven manner. Current practices are often overly stringent and especially untenable on certain types of tissues, such as archived tumor tissues, or tissues associated with mitochondrial function, such as kidney tissue [1]. We propose a data-driven QC metric (miQC) that jointly models both the proportion of reads mapping to mtDNA genes and the number of detected genes with mixture models in a probabilistic framework to predict the low-quality cells in a given dataset. We demonstrate how our QC metric easily adapts to different types of single-cell datasets to remove low-quality cells while preserving high-quality cells that can be used for downstream analyses. Our software package is available at https://bioconductor.org/packages/miQC.

Prioritization of cell types responsive to biological perturbations in single-cell data with Augur

Advances in single-cell genomics now enable large-scale comparisons of cell states across two or more experimental conditions. Numerous statistical tools are available to identify individual genes, proteins or chromatin regions that differ between conditions, but many experiments require inferences at the level of cell types, as opposed to individual analytes. We developed Augur to prioritize the cell types within a complex tissue that are most responsive to an experimental perturbation. In this protocol, we outline the application of Augur to single-cell RNA-seq data, proceeding from a genes-by-cells count matrix to a list of cell types ranked on the basis of their separability following a perturbation. We provide detailed instructions to enable investigators with limited experience in computational biology to perform cell-type prioritization within their own datasets and visualize the results. Moreover, we demonstrate the application of Augur in several more specialized workflows, including the use of RNA velocity for acute perturbations, experimental designs with multiple conditions, differential prioritization between two comparisons, and single-cell transcriptome imaging data. For a dataset containing on the order of 20,000 genes and 20 cell types, this protocol typically takes 1-4 h to complete.

Pheniqs 2.0: accurate, high-performance Bayesian decoding and confidence estimation for combinatorial barcode indexing

Background

Systems biology increasingly relies on deep sequencing with combinatorial index tags to associate biological sequences with their sample, cell, or molecule of origin. Accurate data interpretation depends on the ability to classify sequences based on correct decoding of these combinatorial barcodes. The probability of correct decoding is influenced by both sequence quality and the number and arrangement of barcodes. The rising complexity of experimental designs calls for a probability model that accounts for both sequencing errors and random noise, generalizes to multiple combinatorial tags, and can handle any barcoding scheme. The needs for reproducibility and community benchmark standards demand a peer-reviewed tool that preserves decoding quality scores and provides tunable control over classification confidence that balances precision and recall. Moreover, continuous improvements in sequencing throughput require a fast, parallelized and scalable implementation.

Results and discussion

We developed a flexible, robustly engineered software that performs probabilistic decoding and supports arbitrarily complex barcoding designs. Pheniqs computes the full posterior decoding error probability of observed barcodes by consulting basecalling quality scores and prior distributions, and reports sequences and confidence scores in Sequence Alignment/Map (SAM) fields. The product of posteriors for multiple independent barcodes provides an overall confidence score for each read. Pheniqs achieves greater accuracy than minimum edit distance or simple maximum likelihood estimation, and it scales linearly with core count to enable the classification of > 11 billion reads in 1 h 15 m using < 50 megabytes of memory. Pheniqs has been in production use for seven years in our genomics core facility.

Conclusion

We introduce a computationally efficient software that implements both probabilistic and minimum distance decoders and show that decoding barcodes using posterior probabilities is more accurate than available methods. Pheniqs allows fine-tuning of decoding sensitivity using intuitive confidence thresholds and is extensible with alternative decoders and new error models. Any arbitrary arrangement of barcodes is easily configured, enabling computation of combinatorial confidence scores for any barcoding strategy. An optimized multithreaded implementation assures that Pheniqs is faster and scales better with complex barcode sets than existing tools. Support for POSIX streams and multiple sequencing formats enables easy integration with automated analysis pipelines.

Supplementary Information

The online version supplementary material available at 10.1186/s12859-021-04267-5.

Modular, efficient and constant-memory single-cell RNA-seq preprocessing

We describe a workflow for preprocessing of single-cell RNA-sequencing data that balances efficiency and accuracy. Our workflow is based on the kallisto and bustools programs, and is near optimal in speed with a constant memory requirement providing scalability for arbitrarily large datasets. The workflow is modular, and we demonstrate its flexibility by showing how it can be used for RNA velocity analyses.

BP4RNAseq: a babysitter package for retrospective and newly generated RNA-seq data analyses using both alignment-based and alignment-free quantification method

SUMMARY

Processing raw reads of RNA-sequencing (RNA-seq) data, no matter public or newly sequenced data, involves a lot of specialized tools and technical configurations that are often unfamiliar and time-consuming to learn for non-bioinformatics researchers. Here, we develop the R package BP4RNAseq, which integrates the state-of-art tools from both alignment-based and alignment-free quantification workflows. The BP4RNAseq package is a highly automated tool using an optimized pipeline to improve the sensitivity and accuracy of RNA-seq analyses. It can take only two non-technical parameters and output six formatted gene expression quantification at gene and transcript levels. The package applies to both retrospective and newly generated bulk RNA-seq data analyses and is also applicable for single-cell RNA-seq analyses. It, therefore, greatly facilitates the application of RNA-seq.

AVAILABILITY AND IMPLEMENTATION

The BP4RNAseq package for R and its documentation are freely available at https://github.com/sunshanwen/BP4RNAseq.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Single-Cell RNA Sequencing Analysis of Chicken Anterior Pituitary: A Bird's-Eye View on Vertebrate Pituitary

It is well-established that anterior pituitary contains multiple endocrine cell populations, and each of them can secrete one/two hormone(s) to regulate vital physiological processes of vertebrates. However, the gene expression profiles of each pituitary cell population remains poorly characterized in most vertebrate groups. Here we analyzed the transcriptome of each cell population in adult chicken anterior pituitaries using single-cell RNA sequencing technology. The results showed that: (1) four out of five known endocrine cell clusters have been identified and designated as the lactotrophs, thyrotrophs, corticotrophs, and gonadotrophs, respectively. Somatotrophs were not analyzed in the current study. Each cell cluster can express at least one known endocrine hormone, and novel marker genes (e.g., CD24 and HSPB1 in lactotrophs, NPBWR2 and NDRG1 in corticotrophs; DIO2 and SOUL in thyrotrophs, C5H11ORF96 and HPGDS in gonadotrophs) are identified. Interestingly, gonadotrophs were shown to abundantly express five peptide hormones: FSH, LH, GRP, CART and RLN3; (2) four non-endocrine/secretory cell types, including endothelial cells (expressing IGFBP7 and CFD) and folliculo-stellate cells (FS-cells, expressing S100A6 and S100A10), were identified in chicken anterior pituitaries. Among them, FS-cells can express many growth factors, peptides (e.g., WNT5A, HBEGF, Activins, VEGFC, NPY, and BMP4), and progenitor/stem cell-associated genes (e.g., Notch signaling components, CDH1), implying that the FS-cell cluster may act as a paracrine/autocrine signaling center and enrich pituitary progenitor/stem cells; (3) sexually dimorphic expression of many genes were identified in most cell clusters, including gonadotrophs and lactotrophs. Taken together, our data provides a bird's-eye view on the diverse aspects of anterior pituitaries, including cell composition, heterogeneity, cell-to-cell communication, and gene expression profiles, which facilitates our comprehensive understanding of vertebrate pituitary biology.

VIRTUS: a pipeline for comprehensive virus analysis from conventional RNA-seq data

Summary

The possibility that RNA transcripts from clinical samples contain plenty of virus RNAs has not been pursued actively so far. We here developed a new tool for analyzing virus-transcribed mRNAs, not virus copy numbers, in the data of bulk and single-cell RNA-sequencing of human cells. Our pipeline, named VIRTUS (VIRal Transcript Usage Sensor), was able to detect 762 viruses including herpesviruses, retroviruses and even SARS-CoV-2 (COVID-19), and quantify their transcripts in the sequence data. This tool thus enabled simultaneously detecting infected cells, the composition of multiple viruses within the cell, and the endogenous host-gene expression profile of the cell. This bioinformatics method would be instrumental in addressing the possible effects of covertly infecting viruses on certain diseases and developing new treatments to target such viruses.

Availability and implementation

: VIRTUS is implemented using Common Workflow Language and Docker under a CC-NC license. VIRTUS is freely available at https://github.com/yyoshiaki/VIRTUS.

Supplementary information

Supplementary data are available at Bioinformatics online.

Integrated analysis of multimodal single-cell data

The simultaneous measurement of multiple modalities represents an exciting frontier for single-cell genomics and necessitates computational methods that can define cellular states based on multimodal data. Here, we introduce "weighted-nearest neighbor" analysis, an unsupervised framework to learn the relative utility of each data type in each cell, enabling an integrative analysis of multiple modalities. We apply our procedure to a CITE-seq dataset of 211,000 human peripheral blood mononuclear cells (PBMCs) with panels extending to 228 antibodies to construct a multimodal reference atlas of the circulating immune system. Multimodal analysis substantially improves our ability to resolve cell states, allowing us to identify and validate previously unreported lymphoid subpopulations. Moreover, we demonstrate how to leverage this reference to rapidly map new datasets and to interpret immune responses to vaccination and coronavirus disease 2019 (COVID-19). Our approach represents a broadly applicable strategy to analyze single-cell multimodal datasets and to look beyond the transcriptome toward a unified and multimodal definition of cellular identity.

Statistical Modeling of High Dimensional Counts

Statistical modeling of count data from RNA sequencing (RNA-seq) experiments is important for proper interpretation of results. Here I will describe how count data can be modeled using count distributions, or alternatively analyzed using nonparametric methods. I will focus on basic routines for performing data input, scaling/normalization, visualization, and statistical testing to determine sets of features where the counts reflect differences in gene expression across samples. Finally, I discuss limitations and possible extensions to the models presented here.

UMIc: A Preprocessing Method for UMI Deduplication and Reads Correction

A recent refinement in high-throughput sequencing involves the incorporation of unique molecular identifiers (UMIs), which are random oligonucleotide barcodes, on the library preparation steps. A UMI adds a unique identity to different DNA/RNA input molecules through polymerase chain reaction (PCR) amplification, thus reducing bias of this step. Here, we propose an alignment free framework serving as a preprocessing step of fastq files, called UMIc, for deduplication and correction of reads building consensus sequences from each UMI. Our approach takes into account the frequency and the Phred quality of nucleotides and the distances between the UMIs and the actual sequences. We have tested the tool using different scenarios of UMI-tagged library data, having in mind the aspect of a wide application. UMIc is an open-source tool implemented in R and is freely available from https://github.com/BiodataAnalysisGroup/UMIc.

satuRn: Scalable analysis of differential transcript usage for bulk and single-cell RNA-sequencing applications

Alternative splicing produces multiple functional transcripts from a single gene. Dysregulation of splicing is known to be associated with disease and as a hallmark of cancer. Existing tools for differential transcript usage (DTU) analysis either lack in performance, cannot account for complex experimental designs or do not scale to massive single-cell transcriptome sequencing (scRNA-seq) datasets. We introduce satuRn, a fast and flexible quasi-binomial generalized linear modelling framework that is on par with the best performing DTU methods from the bulk RNA-seq realm, while providing good false discovery rate control, addressing complex experimental designs, and scaling to scRNA-seq applications.

scSNV: accurate dscRNA-seq SNV co-expression analysis using duplicate tag collapsing

Identifying single nucleotide variants has become common practice for droplet-based single-cell RNA-seq experiments; however, presently, a pipeline does not exist to maximize variant calling accuracy. Furthermore, molecular duplicates generated in these experiments have not been utilized to optimally detect variant co-expression. Herein, we introduce scSNV designed from the ground up to "collapse" molecular duplicates and accurately identify variants and their co-expression. We demonstrate that scSNV is fast, with a reduced false-positive variant call rate, and enables the co-detection of genetic variants and A>G RNA edits across twenty-two samples.

Cell-level metadata are indispensable for documenting single-cell sequencing datasets

Single-cell RNA sequencing (scRNA-seq) provides an unprecedented view of cellular diversity of biological systems. However, across the thousands of publications and datasets generated using this technology, we estimate that only a minority (<25%) of studies provide cell-level metadata information containing identified cell types and related findings of the published dataset. Metadata omission hinders reproduction, exploration, validation, and knowledge transfer and is a common problem across journals, data repositories, and publication dates. We encourage investigators, reviewers, journals, and data repositories to improve their standards and ensure proper documentation of these valuable datasets.

Most Gene Expression Omnibus (GEO) depositions of single-cell mRNA sequencing data do not include cell-level metadata generated by typical analysis pipelines; this Essay maintains that this omission greatly hinders reproduction, exploration, validation, and knowledge transfer and is a common problem across journals, data repositories, and publication dates.

Quantifying alternative polyadenylation in RNAseq data with LABRAT

Alternative polyadenylation (APA) generates transcript isoforms that differ in their 3' UTR content and may therefore be subject to different regulatory fates. Although the existence of APA has been known for decades, quantification of APA isoforms from high-throughput RNA sequencing data has been difficult. To facilitate the study of APA in large datasets, we developed an APA quantification technique called LABRAT (Lightweight Alignment-Based Reckoning of Alternative Three-prime ends). LABRAT leverages modern transcriptome quantification approaches to determine the relative abundances of APA isoforms. In this manuscript we describe how LABRAT produces its calculations, provide a step-by-step protocol for its use, and demonstrate its ability to quantify APA in single-cell RNAseq data.

Multi-cohort analysis of host immune response identifies conserved protective and detrimental modules associated with severity across viruses

Viral infections induce a conserved host response distinct from bacterial infections. We hypothesized that the conserved response is associated with disease severity and is distinct between patients with different outcomes. To test this, we integrated 4,780 blood transcriptome profiles from patients aged 0 to 90 years infected with one of 16 viruses, including SARS-CoV-2, Ebola, chikungunya, and influenza, across 34 cohorts from 18 countries, and single-cell RNA sequencing profiles of 702,970 immune cells from 289 samples across three cohorts. Severe viral infection was associated with increased hematopoiesis, myelopoiesis, and myeloid-derived suppressor cells. We identified protective and detrimental gene modules that defined distinct trajectories associated with mild versus severe outcomes. The interferon response was decoupled from the protective host response in patients with severe outcomes. These findings were consistent, irrespective of age and virus, and provide insights to accelerate the development of diagnostics and host-directed therapies to improve global pandemic preparedness.

Emerging Role of ODC1 in Neurodevelopmental Disorders and Brain Development

Ornithine decarboxylase 1 (ODC1 gene) has been linked through gain-of-function variants to a rare disease featuring developmental delay, alopecia, macrocephaly, and structural brain anomalies. ODC1 has been linked to additional diseases like cancer, with growing evidence for neurological contributions to schizophrenia, mood disorders, anxiety, epilepsy, learning, and suicidal behavior. The evidence of ODC1 connection to neural disorders highlights the need for a systematic analysis of ODC1 genotype-to-phenotype associations. An analysis of variants from ClinVar, Geno2MP, TOPMed, gnomAD, and COSMIC revealed an intellectual disability and seizure connected loss-of-function variant, ODC G84R (rs138359527, NC_000002.12:g.10444500C > T). The missense variant is found in ~1% of South Asian individuals and results in 2.5-fold decrease in enzyme function. Expression quantitative trait loci (eQTLs) reveal multiple functionally annotated, non-coding variants regulating ODC1 that associate with psychiatric/neurological phenotypes. Further dissection of RNA-Seq during fetal brain development and within cerebral organoids showed an association of ODC1 expression with cell proliferation of neural progenitor cells, suggesting gain-of-function variants with neural over-proliferation and loss-of-function variants with neural depletion. The linkage from the expression data of ODC1 in early neural progenitor proliferation to phenotypes of neurodevelopmental delay and to the connection of polyamine metabolites in brain function establish ODC1 as a bona fide neurodevelopmental disorder gene.

Integrated RNA Sequencing Reveals Epigenetic Impacts of Diesel Particulate Matter Exposure in Human Cerebral Organoids

Autism spectrum disorder (ASD) manifests early in childhood. While genetic variants increase risk for ASD, a growing body of literature has established that in utero chemical exposures also contribute to ASD risk. These chemicals include air-based pollutants like diesel particulate matter (DPM). A combination of single-cell and direct transcriptomics of DPM-exposed human-induced pluripotent stem cell-derived cerebral organoids revealed toxicogenomic effects of DPM exposure during fetal brain development. Direct transcriptomics, sequencing RNA bases via Nanopore, revealed that cerebral organoids contain extensive RNA modifications, with DPM-altering cytosine methylation in oxidative mitochondrial transcripts expressed in outer radial glia cells. Single-cell transcriptomics further confirmed an oxidative phosphorylation change in cell groups such as outer radial glia upon DPM exposure. This approach highlights how DPM exposure perturbs normal mitochondrial function and cellular respiration during early brain development, which may contribute to developmental disorders like ASD by altering neurodevelopment.

Dissecting the Role of BET Bromodomain Proteins BRD2 and BRD4 in Human NK Cell Function

Natural killer (NK) cells are innate lymphocytes that play a pivotal role in the immune surveillance and elimination of transformed or virally infected cells. Using a chemo-genetic approach, we identify BET bromodomain containing proteins BRD2 and BRD4 as central regulators of NK cell functions, including direct cytokine secretion, NK cell contact-dependent inflammatory cytokine secretion from monocytes as well as NK cell cytolytic functions. We show that both BRD2 and BRD4 control inflammatory cytokine production in NK cells isolated from healthy volunteers and from rheumatoid arthritis patients. In contrast, knockdown of BRD4 but not of BRD2 impairs NK cell cytolytic responses, suggesting BRD4 as critical regulator of NK cell mediated tumor cell elimination. This is supported by pharmacological targeting where the first-generation pan-BET bromodomain inhibitor JQ1(+) displays anti-inflammatory effects and inhibit tumor cell eradication, while the novel bivalent BET bromodomain inhibitor AZD5153, which shows differential activity towards BET family members, does not. Given the important role of both cytokine-mediated inflammatory microenvironment and cytolytic NK cell activities in immune-oncology therapies, our findings present a compelling argument for further clinical investigation.

EKLF/KLF1 expression defines a unique macrophage subset during mouse erythropoiesis

Erythroblastic islands are a specialized niche that contain a central macrophage surrounded by erythroid cells at various stages of maturation. However, identifying the precise genetic and transcriptional control mechanisms in the island macrophage remains difficult due to macrophage heterogeneity. Using unbiased global sequencing and directed genetic approaches focused on early mammalian development, we find that fetal liver macrophages exhibit a unique expression signature that differentiates them from erythroid and adult macrophage cells. The importance of erythroid Krüppel-like factor (EKLF)/KLF1 in this identity is shown by expression analyses in EKLF-/- and in EKLF-marked macrophage cells. Single-cell sequence analysis simplifies heterogeneity and identifies clusters of genes important for EKLF-dependent macrophage function and novel cell surface biomarkers. Remarkably, this singular set of macrophage island cells appears transiently during embryogenesis. Together, these studies provide a detailed perspective on the importance of EKLF in the establishment of the dynamic gene expression network within erythroblastic islands in the developing embryo and provide the means for their efficient isolation.

Network Modeling in Biology: Statistical Methods for Gene and Brain Networks

The rise of network data in many different domains has offered researchers new insight into the problem of modeling complex systems and propelled the development of numerous innovative statistical methodologies and computational tools. In this paper, we primarily focus on two types of biological networks, gene networks and brain networks, where statistical network modeling has found both fruitful and challenging applications. Unlike other network examples such as social networks where network edges can be directly observed, both gene and brain networks require careful estimation of edges using covariates as a first step. We provide a discussion on existing statistical and computational methods for edge esitimation and subsequent statistical inference problems in these two types of biological networks.

Compression of quantification uncertainty for scRNA-seq counts

Motivation

Quantification estimates of gene expression from single-cell RNA-seq (scRNA-seq) data have inherent uncertainty due to reads that map to multiple genes. Many existing scRNA-seq quantification pipelines ignore multi-mapping reads and therefore underestimate expected read counts for many genes. alevin accounts for multi-mapping reads and allows for the generation of ‘inferential replicates’, which reflect quantification uncertainty. Previous methods have shown improved performance when incorporating these replicates into statistical analyses, but storage and use of these replicates increases computation time and memory requirements.

Results

We demonstrate that storing only the mean and variance from a set of inferential replicates (‘compression’) is sufficient to capture gene-level quantification uncertainty, while reducing disk storage to as low as 9% of original storage, and memory usage when loading data to as low as 6%. Using these values, we generate ‘pseudo-inferential’ replicates from a negative binomial distribution and propose a general procedure for incorporating these replicates into a proposed statistical testing framework. When applying this procedure to trajectory-based differential expression analyses, we show false positives are reduced by more than a third for genes with high levels of quantification uncertainty. We additionally extend the Swish method to incorporate pseudo-inferential replicates and demonstrate improvements in computation time and memory usage without any loss in performance. Lastly, we show that discarding multi-mapping reads can result in significant underestimation of counts for functionally important genes in a real dataset.

Availability and implementation

makeInfReps and splitSwish are implemented in the R/Bioconductor fishpond package available at https://bioconductor.org/packages/fishpond. Analyses and simulated datasets can be found in the paper’s GitHub repo at https://github.com/skvanburen/scUncertaintyPaperCode.

Supplementary information

Supplementary data are available at Bioinformatics online.

Single-Cell RNA Sequencing to Disentangle the Blood System

The blood system is often represented as a tree-like structure with stem cells that give rise to mature blood cell types through a series of demarcated steps. Although this representation has served as a model of hierarchical tissue organization for decades, single-cell technologies are shedding new light on the abundance of cell type intermediates and the molecular mechanisms that ensure balanced replenishment of differentiated cells. In this Brief Review, we exemplify new insights into blood cell differentiation generated by single-cell RNA sequencing, summarize considerations for the application of this technology, and highlight innovations that are leading the way to understand hematopoiesis at the resolution of single cells. Graphic Abstract: A graphic abstract is available for this article.

Preprocessing choices affect RNA velocity results for droplet scRNA-seq data

Experimental single-cell approaches are becoming widely used for many purposes, including investigation of the dynamic behaviour of developing biological systems. Consequently, a large number of computational methods for extracting dynamic information from such data have been developed. One example is RNA velocity analysis, in which spliced and unspliced RNA abundances are jointly modeled in order to infer a 'direction of change' and thereby a future state for each cell in the gene expression space. Naturally, the accuracy and interpretability of the inferred RNA velocities depend crucially on the correctness of the estimated abundances. Here, we systematically compare five widely used quantification tools, in total yielding thirteen different quantification approaches, in terms of their estimates of spliced and unspliced RNA abundances in five experimental droplet scRNA-seq data sets. We show that there are substantial differences between the quantifications obtained from different tools, and identify typical genes for which such discrepancies are observed. We further show that these abundance differences propagate to the downstream analysis, and can have a large effect on estimated velocities as well as the biological interpretation. Our results highlight that abundance quantification is a crucial aspect of the RNA velocity analysis workflow, and that both the definition of the genomic features of interest and the quantification algorithm itself require careful consideration.

Decoding myofibroblast origins in human kidney fibrosis

Kidney fibrosis is the hallmark of chronic kidney disease progression; however, at present no antifibrotic therapies exist^1-3. The origin, functional heterogeneity and regulation of scar-forming cells that occur during human kidney fibrosis remain poorly understood^1,2,4. Here, using single-cell RNA sequencing, we profiled the transcriptomes of cells from the proximal and non-proximal tubules of healthy and fibrotic human kidneys to map the entire human kidney. This analysis enabled us to map all matrix-producing cells at high resolution, and to identify distinct subpopulations of pericytes and fibroblasts as the main cellular sources of scar-forming myofibroblasts during human kidney fibrosis. We used genetic fate-tracing, time-course single-cell RNA sequencing and ATAC-seq (assay for transposase-accessible chromatin using sequencing) experiments in mice, and spatial transcriptomics in human kidney fibrosis, to shed light on the cellular origins and differentiation of human kidney myofibroblasts and their precursors at high resolution. Finally, we used this strategy to detect potential therapeutic targets, and identified NKD2 as a myofibroblast-specific target in human kidney fibrosis.

Transcriptomic Bioinformatic Analyses of Atria Uncover Involvement of Pathways Related to Strain and Post-translational Modification of Collagen in Increased Atrial Fibrillation Vulnerability in Intensely Exercised Mice

Atrial Fibrillation (AF) is the most common supraventricular tachyarrhythmia that is typically associated with cardiovascular disease (CVD) and poor cardiovascular health. Paradoxically, endurance athletes are also at risk for AF. While it is well-established that persistent AF is associated with atrial fibrosis, hypertrophy and inflammation, intensely exercised mice showed similar adverse atrial changes and increased AF vulnerability, which required tumor necrosis factor (TNF) signaling, even though ventricular structure and function improved. To identify some of the molecular factors underlying the chamber-specific and TNF-dependent atrial changes induced by exercise, we performed transcriptome analyses of hearts from wild-type and TNF-knockout mice following exercise for 2 days, 2 or 6 weeks of exercise. Consistent with the central role of atrial stretch arising from elevated venous pressure in AF promotion, all 3 time points were associated with differential regulation of genes in atria linked to mechanosensing (focal adhesion kinase, integrins and cell-cell communications), extracellular matrix (ECM) and TNF pathways, with TNF appearing to play a permissive, rather than causal, role in gene changes. Importantly, mechanosensing/ECM genes were only enriched, along with tubulin- and hypertrophy-related genes after 2 days of exercise while being downregulated at 2 and 6 weeks, suggesting that early reactive strain-dependent remodeling with exercise yields to compensatory adjustments. Moreover, at the later time points, there was also downregulation of both collagen genes and genes involved in collagen turnover, a pattern mirroring aging-related fibrosis. By comparison, twofold fewer genes were differentially regulated in ventricles vs. atria, independently of TNF. Our findings reveal that exercise promotes TNF-dependent atrial transcriptome remodeling of ECM/mechanosensing pathways, consistent with increased preload and atrial stretch seen with exercise. We propose that similar preload-dependent mechanisms are responsible for atrial changes and AF in both CVD patients and athletes.

Goals and approaches for each processing step for single-cell RNA sequencing data

Single-cell RNA sequencing (scRNA-seq) has enabled researchers to study gene expression at the cellular level. However, due to the extremely low levels of transcripts in a single cell and technical losses during reverse transcription, gene expression at a single-cell resolution is usually noisy and highly dimensional; thus, statistical analyses of single-cell data are a challenge. Although many scRNA-seq data analysis tools are currently available, a gold standard pipeline is not available for all datasets. Therefore, a general understanding of bioinformatics and associated computational issues would facilitate the selection of appropriate tools for a given set of data. In this review, we provide an overview of the goals and most popular computational analysis tools for the quality control, normalization, imputation, feature selection and dimension reduction of scRNA-seq data.

A single-cell RNA-sequencing training and analysis suite using the Galaxy framework

Background

The vast ecosystem of single-cell RNA-sequencing tools has until recently been plagued by an excess of diverging analysis strategies, inconsistent file formats, and compatibility issues between different software suites. The uptake of 10x Genomics datasets has begun to calm this diversity, and the bioinformatics community leans once more towards the large computing requirements and the statistically driven methods needed to process and understand these ever-growing datasets.

Results

Here we outline several Galaxy workflows and learning resources for single-cell RNA-sequencing, with the aim of providing a comprehensive analysis environment paired with a thorough user learning experience that bridges the knowledge gap between the computational methods and the underlying cell biology. The Galaxy reproducible bioinformatics framework provides tools, workflows, and trainings that not only enable users to perform 1-click 10x preprocessing but also empower them to demultiplex raw sequencing from custom tagged and full-length sequencing protocols. The downstream analysis supports a range of high-quality interoperable suites separated into common stages of analysis: inspection, filtering, normalization, confounder removal, and clustering. The teaching resources cover concepts from computer science to cell biology. Access to all resources is provided at the singlecell.usegalaxy.eu portal.

Conclusions

The reproducible and training-oriented Galaxy framework provides a sustainable high-performance computing environment for users to run flexible analyses on both 10x and alternative platforms. The tutorials from the Galaxy Training Network along with the frequent training workshops hosted by the Galaxy community provide a means for users to learn, publish, and teach single-cell RNA-sequencing analysis.

Cumulus provides cloud-based data analysis for large-scale single-cell and single-nucleus RNA-seq

Massively parallel single-cell and single-nucleus RNA sequencing has opened the way to systematic tissue atlases in health and disease, but as the scale of data generation is growing, so is the need for computational pipelines for scaled analysis. Here we developed Cumulus-a cloud-based framework for analyzing large-scale single-cell and single-nucleus RNA sequencing datasets. Cumulus combines the power of cloud computing with improvements in algorithm and implementation to achieve high scalability, low cost, user-friendliness and integrated support for a comprehensive set of features. We benchmark Cumulus on the Human Cell Atlas Census of Immune Cells dataset of bone marrow cells and show that it substantially improves efficiency over conventional frameworks, while maintaining or improving the quality of results, enabling large-scale studies.

Computational Methods for Single-Cell RNA Sequencing

Single-cell RNA sequencing (scRNA-seq) has provided a high-dimensional catalog of millions of cells across species and diseases. These data have spurred the development of hundreds of computational tools to derive novel biological insights. Here, we outline the components of scRNA-seq analytical pipelines and the computational methods that underlie these steps. We describe available methods, highlight well-executed benchmarking studies, and identify opportunities for additional benchmarking studies and computational methods. As the biochemical approaches for single-cell omics advance, we propose coupled development of robust analytical pipelines suited for the challenges that new data present and principled selection of analytical methods that are suited for the biological questions to be addressed.

RainDrop: Rapid activation matrix computation for droplet-based single-cell RNA-seq reads

Background

Obtaining data from single-cell transcriptomic sequencing allows for the investigation of cell-specific gene expression patterns, which could not be addressed a few years ago. With the advancement of droplet-based protocols the number of studied cells continues to increase rapidly. This establishes the need for software tools for efficient processing of the produced large-scale datasets. We address this need by presenting RainDrop for fast gene-cell count matrix computation from single-cell RNA-seq data produced by 10x Genomics Chromium technology.

Results

RainDrop can process single-cell transcriptomic datasets consisting of 784 million reads sequenced from around 8.000 cells in less than 40 minutes on a standard workstation. It significantly outperforms the established Cell Ranger pipeline and the recently introduced Alevin tool in terms of runtime by a maximal (average) speedup of 30.4 (22.6) and 3.5 (2.4), respectively, while keeping high agreements of the generated results.

Conclusions

RainDrop is a software tool for highly efficient processing of large-scale droplet-based single-cell RNA-seq datasets on standard workstations written in C++. It is available at https://gitlab.rlp.net/stnieble/raindrop.

A Bayesian framework for inter-cellular information sharing improves dscRNA-seq quantification

Motivation

Droplet-based single-cell RNA-seq (dscRNA-seq) data are being generated at an unprecedented pace, and the accurate estimation of gene-level abundances for each cell is a crucial first step in most dscRNA-seq analyses. When pre-processing the raw dscRNA-seq data to generate a count matrix, care must be taken to account for the potentially large number of multi-mapping locations per read. The sparsity of dscRNA-seq data, and the strong 3’ sampling bias, makes it difficult to disambiguate cases where there is no uniquely mapping read to any of the candidate target genes.

Results

We introduce a Bayesian framework for information sharing across cells within a sample, or across multiple modalities of data using the same sample, to improve gene quantification estimates for dscRNA-seq data. We use an anchor-based approach to connect cells with similar gene-expression patterns, and learn informative, empirical priors which we provide to alevin’s gene multi-mapping resolution algorithm. This improves the quantification estimates for genes with no uniquely mapping reads (i.e. when there is no unique intra-cellular information). We show our new model improves the per cell gene-level estimates and provides a principled framework for information sharing across multiple modalities. We test our method on a combination of simulated and real datasets under various setups.

Availability and implementation

The information sharing model is included in alevin and is implemented in C++14. It is available as open-source software, under GPL v3, at https://github.com/COMBINE-lab/salmon as of version 1.1.0.

Reversing a model of Parkinson's disease with in situ converted nigral neurons

Parkinson's disease is characterized by loss of dopamine neurons in the substantia nigra¹. Similar to other major neurodegenerative disorders, there are no disease-modifying treatments for Parkinson's disease. While most treatment strategies aim to prevent neuronal loss or protect vulnerable neuronal circuits, a potential alternative is to replace lost neurons to reconstruct disrupted circuits². Here we report an efficient one-step conversion of isolated mouse and human astrocytes to functional neurons by depleting the RNA-binding protein PTB (also known as PTBP1). Applying this approach to the mouse brain, we demonstrate progressive conversion of astrocytes to new neurons that innervate into and repopulate endogenous neural circuits. Astrocytes from different brain regions are converted to different neuronal subtypes. Using a chemically induced model of Parkinson's disease in mouse, we show conversion of midbrain astrocytes to dopaminergic neurons, which provide axons to reconstruct the nigrostriatal circuit. Notably, re-innervation of striatum is accompanied by restoration of dopamine levels and rescue of motor deficits. A similar reversal of disease phenotype is also accomplished by converting astrocytes to neurons using antisense oligonucleotides to transiently suppress PTB. These findings identify a potentially powerful and clinically feasible approach to treating neurodegeneration by replacing lost neurons.

A scalable SCENIC workflow for single-cell gene regulatory network analysis

This protocol explains how to perform a fast SCENIC analysis alongside standard best practices steps on single-cell RNA-sequencing data using software containers and Nextflow pipelines. SCENIC reconstructs regulons (i.e., transcription factors and their target genes) assesses the activity of these discovered regulons in individual cells and uses these cellular activity patterns to find meaningful clusters of cells. Here we present an improved version of SCENIC with several advances. SCENIC has been refactored and reimplemented in Python (pySCENIC), resulting in a tenfold increase in speed, and has been packaged into containers for ease of use. It is now also possible to use epigenomic track databases, as well as motifs, to refine regulons. In this protocol, we explain the different steps of SCENIC: the workflow starts from the count matrix depicting the gene abundances for all cells and consists of three stages. First, coexpression modules are inferred using a regression per-target approach (GRNBoost2). Next, the indirect targets are pruned from these modules using cis-regulatory motif discovery (cisTarget). Lastly, the activity of these regulons is quantified via an enrichment score for the regulon's target genes (AUCell). Nonlinear projection methods can be used to display visual groupings of cells based on the cellular activity patterns of these regulons. The results can be exported as a loom file and visualized in the SCope web application. This protocol is illustrated on two use cases: a peripheral blood mononuclear cell data set and a panel of single-cell RNA-sequencing cancer experiments. For a data set of 10,000 genes and 50,000 cells, the pipeline runs in <2 h.

Fast analysis of scATAC-seq data using a predefined set of genomic regions

Background: Analysis of scATAC-seq data has been recently scaled to thousands of cells. While processing of other types of single cell data was boosted by the implementation of alignment-free techniques, pipelines available to process scATAC-seq data still require large computational resources. We propose here an approach based on pseudoalignment, which reduces the execution times and hardware needs at little cost for precision. Methods: Public data for 10k PBMC were downloaded from 10x Genomics web site. Reads were aligned to various references derived from DNase I Hypersensitive Sites (DHS) using kallisto and quantified with bustools. We compared our results with the ones publicly available derived by cellranger-atac. We subsequently tested our approach on scATAC-seq data for K562 cell line. Results: We found that kallisto does not introduce biases in quantification of known peaks; cells groups identified are consistent with the ones identified from standard method. We also found that cell identification is robust when analysis is performed using DHS-derived reference in place of de novo identification of ATAC peaks. Lastly, we found that our approach is suitable for reliable quantification of gene activity based on scATAC-seq signal, thus allows for efficient labelling of cell groups based on marker genes. Conclusions: Analysis of scATAC-seq data by means of kallisto produces results in line with standard pipelines while being considerably faster; using a set of known DHS sites as reference does not affect the ability to characterize the cell populations.

Tximeta: Reference sequence checksums for provenance identification in RNA-seq

Correct annotation metadata is critical for reproducible and accurate RNA-seq analysis. When files are shared publicly or among collaborators with incorrect or missing annotation metadata, it becomes difficult or impossible to reproduce bioinformatic analyses from raw data. It also makes it more difficult to locate the transcriptomic features, such as transcripts or genes, in their proper genomic context, which is necessary for overlapping expression data with other datasets. We provide a solution in the form of an R/Bioconductor package tximeta that performs numerous annotation and metadata gathering tasks automatically on behalf of users during the import of transcript quantification files. The correct reference transcriptome is identified via a hashed checksum stored in the quantification output, and key transcript databases are downloaded and cached locally. The computational paradigm of automatically adding annotation metadata based on reference sequence checksums can greatly facilitate genomic workflows, by helping to reduce overhead during bioinformatic analyses, preventing costly bioinformatic mistakes, and promoting computational reproducibility. The tximeta package is available at https://bioconductor.org/packages/tximeta.

Orchestrating single-cell analysis with Bioconductor

Recent technological advancements have enabled the profiling of a large number of genome-wide features in individual cells. However, single-cell data present unique challenges that require the development of specialized methods and software infrastructure to successfully derive biological insights. The Bioconductor project has rapidly grown to meet these demands, hosting community-developed open-source software distributed as R packages. Featuring state-of-the-art computational methods, standardized data infrastructure and interactive data visualization tools, we present an overview and online book (https://osca.bioconductor.org) of single-cell methods for prospective users.

Expression Atlas update: from tissues to single cells

Expression Atlas is EMBL-EBI's resource for gene and protein expression. It sources and compiles data on the abundance and localisation of RNA and proteins in various biological systems and contexts and provides open access to this data for the research community. With the increased availability of single cell RNA-Seq datasets in the public archives, we have now extended Expression Atlas with a new added-value service to display gene expression in single cells. Single Cell Expression Atlas was launched in 2018 and currently includes 123 single cell RNA-Seq studies from 12 species. The website can be searched by genes within or across species to reveal experiments, tissues and cell types where this gene is expressed or under which conditions it is a marker gene. Within each study, cells can be visualized using a pre-calculated t-SNE plot and can be coloured by different features or by cell clusters based on gene expression. Within each experiment, there are links to downloadable files, such as RNA quantification matrices, clustering results, reports on protocols and associated metadata, such as assigned cell types.

Algorithms for efficiently collapsing reads with Unique Molecular Identifiers

BACKGROUND

Unique Molecular Identifiers (UMI) are used in many experiments to find and remove PCR duplicates. There are many tools for solving the problem of deduplicating reads based on their finding reads with the same alignment coordinates and UMIs. However, many tools either cannot handle substitution errors, or require expensive pairwise UMI comparisons that do not efficiently scale to larger datasets.

RESULTS

We reformulate the problem of deduplicating UMIs in a manner that enables optimizations to be made, and more efficient data structures to be used. We implement our data structures and optimizations in a tool called UMICollapse, which is able to deduplicate over one million unique UMIs of length 9 at a single alignment position in around 26 s, using only a single thread and much less than 10 GB of memory.

CONCLUSIONS

We present a new formulation of the UMI deduplication problem, and show that it can be solved faster, with more sophisticated data structures.

Nonparametric expression analysis using inferential replicate counts

A primary challenge in the analysis of RNA-seq data is to identify differentially expressed genes or transcripts while controlling for technical biases. Ideally, a statistical testing procedure should incorporate the inherent uncertainty of the abundance estimates arising from the quantification step. Most popular methods for RNA-seq differential expression analysis fit a parametric model to the counts for each gene or transcript, and a subset of methods can incorporate uncertainty. Previous work has shown that nonparametric models for RNA-seq differential expression may have better control of the false discovery rate, and adapt well to new data types without requiring reformulation of a parametric model. Existing nonparametric models do not take into account inferential uncertainty, leading to an inflated false discovery rate, in particular at the transcript level. We propose a nonparametric model for differential expression analysis using inferential replicate counts, extending the existing SAMseq method to account for inferential uncertainty. We compare our method, Swish, with popular differential expression analysis methods. Swish has improved control of the false discovery rate, in particular for transcripts with high inferential uncertainty. We apply Swish to a single-cell RNA-seq dataset, assessing differential expression between sub-populations of cells, and compare its performance to the Wilcoxon test.