scholarly journals GPress: a framework for querying general feature format (GFF) files and expression files in a compressed form

2020 ◽  
Vol 36 (18) ◽  
pp. 4810-4812
Author(s):  
Qingxi Meng ◽  
Idoia Ochoa ◽  
Mikel Hernaez
Keyword(s):  
Single Cell ◽  
Data Streams ◽  
General Feature ◽  
Supplementary Information ◽  
Storage Space ◽  
Supplementary Data ◽  
Rna Seq ◽  
Sequencing Data ◽  
General Feature Format ◽  
Original File

Abstract Motivation Sequencing data are often summarized at different annotation levels for further analysis, generally using the general feature format (GFF) or its descendants, gene transfer format (GTF) and GFF3. Existing utilities for accessing these files, like gffutils and gffread, do not focus on reducing the storage space, significantly increasing it in some cases. We propose GPress, a framework for querying GFF files in a compressed form. GPress can also incorporate and compress expression files from both bulk and single-cell RNA-Seq experiments, supporting simultaneous queries on both the GFF and expression files. In brief, GPress applies transformations to the data which are then compressed with the general lossless compressor BSC. To support queries, GPress compresses the data in blocks and creates several index tables for fast retrieval. Results We tested GPress on several GFF files of different organisms, and showed that it achieves on average a 61% reduction in size with respect to gzip (the current de facto compressor for GFF files) while being able to retrieve all annotations for a given identifier or a range of coordinates in a few seconds (when run in a common laptop). In contrast, gffutils provides faster retrieval but doubles the size of the GFF files. When additionally linking an expression file, we show that GPress can reduce its size by more than 68% when compared to gzip (for both bulk and single-cell RNA-Seq experiments), while still retrieving the information within seconds. Finally, applying BSC to the data streams generated by GPress instead of to the original file shows a size reduction of more than 44% on average. Availability and implementation GPress is freely available at https://github.com/qm2/gpress. Supplementary information Supplementary data are available at Bioinformatics online.

scholarly journals DEsingle for detecting three types of differential expression in single-cell RNA-seq data

2017 ◽  
Author(s):  
Zhun Miao ◽  
Ke Deng ◽  
Xiaowo Wang ◽  
Xuegong Zhang
Keyword(s):  
Single Cell ◽  
Differential Expression ◽  
Negative Binomial ◽  
Single Cells ◽  
R Package ◽  
Supplementary Information ◽  
Binomial Model ◽  
Supplementary Data ◽  
Rna Seq ◽  
Real Zeros

AbstractSummaryThe excessive amount of zeros in single-cell RNA-seq data include “real” zeros due to the on-off nature of gene transcription in single cells and “dropout” zeros due to technical reasons. Existing differential expression (DE) analysis methods cannot distinguish these two types of zeros. We developed an R package DEsingle which employed Zero-Inflated Negative Binomial model to estimate the proportion of real and dropout zeros and to define and detect 3 types of DE genes in single-cell RNA-seq data with higher accuracy.Availability and ImplementationThe R package DEsingle is freely available at https://github.com/miaozhun/DEsingle and is under Bioconductor’s consideration [email protected] informationSupplementary data are available at bioRxiv online.

schex avoids overplotting for large single-cell RNA-sequencing datasets

2019 ◽  
Vol 36 (7) ◽  
pp. 2291-2292 ◽  
Author(s):  
Saskia Freytag ◽  
Ryan Lister
Keyword(s):  
Single Cell ◽  
Rna Sequencing ◽  
R Package ◽  
Supplementary Information ◽  
Supplementary Data ◽  
Sequencing Data ◽  
Single Cell Rna Sequencing

Abstract Summary Due to the scale and sparsity of single-cell RNA-sequencing data, traditional plots can obscure vital information. Our R package schex overcomes this by implementing hexagonal binning, which has the additional advantages of improving speed and reducing storage for resulting plots. Availability and implementation schex is freely available from Bioconductor via http://bioconductor.org/packages/release/bioc/html/schex.html and its development version can be accessed on GitHub via https://github.com/SaskiaFreytag/schex. Supplementary information Supplementary data are available at Bioinformatics online.

scholarly journals SPsimSeq: semi-parametric simulation of bulk and single cell RNA sequencing data

2019 ◽  
Author(s):  
Alemu Takele Assefa ◽  
Jo Vandesompele ◽  
Olivier Thas
Keyword(s):  
Gene Expression ◽  
Single Cell ◽  
Rna Sequencing ◽  
Empirical Distribution ◽  
Supplementary Information ◽  
Rna Seq ◽  
Sequencing Data ◽  
Actual Distribution ◽  
Wide Range ◽  
Single Cell Rna Sequencing

SummarySPsimSeq is a semi-parametric simulation method for bulk and single cell RNA sequencing data. It simulates data from a good estimate of the actual distribution of a given real RNA-seq dataset. In contrast to existing approaches that assume a particular data distribution, our method constructs an empirical distribution of gene expression data from a given source RNA-seq experiment to faithfully capture the data characteristics of real data. Importantly, our method can be used to simulate a wide range of scenarios, such as single or multiple biological groups, systematic variations (e.g. confounding batch effects), and different sample sizes. It can also be used to simulate different gene expression units resulting from different library preparation protocols, such as read counts or UMI counts.Availability and implementationThe R package and associated documentation is available from https://github.com/CenterForStatistics-UGent/SPsimSeq.Supplementary informationSupplementary data are available at bioRχiv online.

scholarly journals Per-sample standardization and asymmetric winsorization lead to accurate clustering of RNA-seq expression profiles

2021 ◽  
Author(s):  
Davide Risso ◽  
Stefano Maria Pagnotta
Keyword(s):  
Single Cell ◽  
Expression Profiles ◽  
Unsupervised Clustering ◽  
Supplementary Information ◽  
Supplementary Data ◽  
Rna Seq ◽  
Data Transformations ◽  
The Impact

Abstract Motivation Data transformations are an important step in the analysis of RNA-seq data. Nonetheless, the impact of transformation on the outcome of unsupervised clustering procedures is still unclear. Results Here, we present an Asymmetric Winsorization per Sample Transformation (AWST), which is robust to data perturbations and removes the need for selecting the most informative genes prior to sample clustering. Our procedure leads to robust and biologically meaningful clusters both in bulk and in single-cell applications. Availability The AWST method is available at https://github.com/drisso/awst. The code to reproduce the analyses is available at https://github.com/drisso/awst\_analysis. Supplementary information Supplementary data are available at Bioinformatics online.

scholarly journals scBatch: batch-effect correction of RNA-seq data through sample distance matrix adjustment

2020 ◽  
Vol 36 (10) ◽  
pp. 3115-3123 ◽  
Author(s):  
Teng Fei ◽  
Tianwei Yu
Keyword(s):  
Single Cell ◽  
Differential Expression Analysis ◽  
Distance Matrix ◽  
Real Data ◽  
R Package ◽  
Batch Effect ◽  
Supplementary Information ◽  
Rna Seq ◽  
Sequencing Data ◽  
Gene Differential Expression

Abstract Motivation Batch effect is a frequent challenge in deep sequencing data analysis that can lead to misleading conclusions. Existing methods do not correct batch effects satisfactorily, especially with single-cell RNA sequencing (RNA-seq) data. Results We present scBatch, a numerical algorithm for batch-effect correction on bulk and single-cell RNA-seq data with emphasis on improving both clustering and gene differential expression analysis. scBatch is not restricted by assumptions on the mechanism of batch-effect generation. As shown in simulations and real data analyses, scBatch outperforms benchmark batch-effect correction methods. Availability and implementation The R package is available at github.com/tengfei-emory/scBatch. The code to generate results and figures in this article is available at github.com/tengfei-emory/scBatch-paper-scripts. Supplementary information Supplementary data are available at Bioinformatics online.

scholarly journals NewWave: a scalable R/Bioconductor package for the dimensionality reduction and batch effect removal of single-cell RNA-seq data

2021 ◽  
Author(s):  
Federico Agostinis ◽  
Chiara Romualdi ◽  
Gabriele Sales ◽  
Davide Risso
Keyword(s):  
Dimensionality Reduction ◽  
Single Cell ◽  
R Package ◽  
Batch Effect ◽  
Supplementary Information ◽  
Bioconductor Package ◽  
Rna Seq ◽  
Sequencing Data ◽  
Bioconductor Project ◽  
Single Cell Rna Sequencing

Summary: We present NewWave, a scalable R/Bioconductor package for the dimensionality reduction and batch effect removal of single-cell RNA sequencing data. To achieve scalability, NewWave uses mini-batch optimization and can work with out-of-memory data, enabling users to analyze datasets with millions of cells. Availability and implementation: NewWave is implemented as an open-source R package available through the Bioconductor project at https://bioconductor.org/packages/NewWave/ Supplementary information: Supplementary data are available at Bioinformatics online.

scholarly journals “GPress: a framework for querying General Feature Format (GFF) files and feature expression files in a compressed form”

2019 ◽  
Author(s):  
Qingxi Meng ◽  
Idoia Ochoa ◽  
Mikel Hernaez
Keyword(s):  
Data Storage ◽  
General Feature ◽  
Storage Space ◽  
Sequencing Data ◽  
Unique Identifier ◽  
Storage And Retrieval ◽  
Data Formats ◽  
General Feature Format ◽  
Efficient Data ◽  
Feature Expression

1Abstract1.1MotivationSequencing data are often summarized at different annotation levels for further analysis. The general feature format (GFF) and its descendants, the gene transfer format (GTF) and GFF3, are the most commonly used data formats for genomic annotations. These files are extensively updated, queried and shared, and hence as the number of generated GFF files increases, efficient data storage and retrieval are becoming increasingly important. Existing GFF utilities for accessing these files, like gffutils and gffread, do not focus on reducing the storage space, significantly increasing it in some cases. Hence, we propose GPress, a framework for querying GFF files in a compressed form. In addition, GPress can also incorporate and compress feature expression files, supporting simultaneous queries on both files.1.2ResultsWe tested GPress on several GFF files of different organisms, and showed that it achieves on average a 98% reduction in size, while being able to retrieve all annotations for a given identifier or a range of coordinates in a few seconds. For example, on a Human GFF file, GPress can find all items with a unique identifier in 2.47 seconds and all items with coordinates within the range of 1,000 to 100,000 in 4.61 seconds. In contrast, gffutils provides faster retrieval but doubles the size of the GFF files. When additionally linking an expression file, we show that GPress can reduce the size of the expression file by more than 92%, while still retrieving the information within seconds. GPress is freely available at https://github.com/qm2/gpress.

scholarly journals The barcode, UMI, set format and BUStools

2019 ◽  
Vol 35 (21) ◽  
pp. 4472-4473 ◽  
Author(s):  
Páll Melsted ◽  
Vasilis Ntranos ◽  
Lior Pachter
Keyword(s):  
Single Cell ◽  
Supplementary Information ◽  
Supplementary Data ◽  
Rna Seq ◽  
Modular Technology ◽  
Rapid Quantification

Abstract Summary We introduce the Barcode-UMI-Set format (BUS) for representing pseudoalignments of reads from single-cell RNA-seq experiments. The format can be used with all single-cell RNA-seq technologies, and we show that BUS files can be efficiently generated. BUStools is a suite of tools for working with BUS files and facilitates rapid quantification and analysis of single-cell RNA-seq data. The BUS format therefore makes possible the development of modular, technology-specific and robust workflows for single-cell RNA-seq analysis. Availability and implementation http://BUStools.github.io/ and http://pachterlab.github.io/kallisto/singlecell.html. Supplementary information Supplementary data are available at Bioinformatics online.

scholarly journals An ultra-sensitive T-cell receptor detection method for TCR-Seq and RNA-Seq data

2020 ◽  
Vol 36 (15) ◽  
pp. 4255-4262
Author(s):  
Si-Yi Chen ◽  
Chun-Jie Liu ◽  
Qiong Zhang ◽  
An-Yuan Guo
Keyword(s):  
T Cell ◽  
Single Cell ◽  
High Performance ◽  
Cell Receptor ◽  
Computational Method ◽  
Read Length ◽  
Supplementary Information ◽  
De Bruijn Graph ◽  
Rna Seq ◽  
Sequencing Data

Abstract Motivation T-cell receptors (TCRs) function to recognize antigens and play vital roles in T-cell immunology. Surveying TCR repertoires by characterizing complementarity-determining region 3 (CDR3) is a key issue. Due to the high diversity of CDR3 and technological limitation, accurate characterization of CDR3 repertoires remains a great challenge. Results We propose a computational method named CATT for ultra-sensitive and precise TCR CDR3 sequences detection. CATT can be applied on TCR sequencing, RNA-Seq and single-cell TCR(RNA)-Seq data to characterize CDR3 repertoires. CATT integrated de Bruijn graph-based micro-assembly algorithm, data-driven error correction model and Bayesian inference algorithm, to self-adaptively and ultra-sensitively characterize CDR3 repertoires with high performance. Benchmark results of datasets from in silico and experimental data demonstrated that CATT showed superior recall and precision compared with existing tools, especially for data with short read length and small size and single-cell sequencing data. Thus, CATT will be a useful tool for TCR analysis in researches of cancer and immunology. Availability and implementation http://bioinfo.life.hust.edu.cn/CATT or https://github.com/GuoBioinfoLab/CATT. Supplementary information Supplementary data are available at Bioinformatics online.

scholarly journals snakePipes: facilitating flexible, scalable and integrative epigenomic analysis

2019 ◽  
Vol 35 (22) ◽  
pp. 4757-4759 ◽  
Author(s):  
Vivek Bhardwaj ◽  
Steffen Heyne ◽  
Katarzyna Sikora ◽  
Leily Rabbani ◽  
Michael Rauer ◽  
...  
Keyword(s):  
Single Cell ◽  
Source Code ◽  
Supplementary Information ◽  
Command Line ◽  
Supplementary Data ◽  
Rna Seq ◽  
Downstream Analysis ◽  
Scalable Analysis

Abstract Summary Due to the rapidly increasing scale and diversity of epigenomic data, modular and scalable analysis workflows are of wide interest. Here we present snakePipes, a workflow package for processing and downstream analysis of data from common epigenomic assays: ChIP-seq, RNA-seq, Bisulfite-seq, ATAC-seq, Hi-C and single-cell RNA-seq. snakePipes enables users to assemble variants of each workflow and to easily install and upgrade the underlying tools, via its simple command-line wrappers and yaml files. Availability and implementation snakePipes can be installed via conda: `conda install -c mpi-ie -c bioconda -c conda-forge snakePipes’. Source code (https://github.com/maxplanck-ie/snakepipes) and documentation (https://snakepipes.readthedocs.io/en/latest/) are available online. Supplementary information Supplementary data are available at Bioinformatics online.

Sign in / Sign up

Export Citation Format

Share Document