Skip to main content

BRIE: Bayesian regression for isoform estimate

Project description

PyPI Docs Build Status

BRIE: Bayesian Regression for Isoform Estimate

About BRIE

Welcome to the new BRIE (>=2.0 or BRIE2), Bayesian Regression for Isoform Estimate, a scalable Bayesian method to accurately identify splicing phenotypes in single-cell RNA-seq experiments and quantify isoform proportions and their uncertainty.

BRIE2 supports the analysis of splicing processes at two molecular levels, either between alternative splicing isoforms or between unspliced and spliced RNAs. In either case, it returns cell-by-event or cell-by-gene matrices of PSI value and its 95% confidence interval (quantification) and the statistics for detecting DAS and DMG on each event or gene:

  1. Differential alternative splicing (DAS): This task is to quantify the proportions of alternative splicing isoforms and to detect DAS between groups of cells or along with a continuous covariate, e.g., pseudotime. BRIE2 is designed for two-isoform splicing events with a focus on exon skipping, but in principle also applicable for mutual exclusion, intron-retaining, alternative poly-A site, 3’ splice site and 5’ splice site.
  2. Differential momentum genes (DMG): This task is to quantify the proportions of unspliced and spliced RNAs in each gene and each cell. Similar to DAS, the DMG is a principled selection of genes that capture heterogeneity in transcriptional kinetics between cell groups, e.g., cell types, or continuous cell covariates, hence may enhance the RNA velocity analyses by focusing on dynamics informed genes.

Installation

BRIE2 is available through PyPI. To install, type the following command line, and add -U for upgrading:

pip install -U brie

Alternatively, you can install from this GitHub repository for the latest (often development) version with the following command line

pip install -U git+https://github.com/huangyh09/brie

In either case, add --user if you don’t have the write permission for your Python environment.

For more instructions, see the installation manual.

Manual and examples

In short, there are two steps for running BRIE2. First, obtain cell-by-gene or cell-by-event count matrices for each isoform. For the exon-skipping event, you can run brie-count, which will return count matrices and hdf5 file for AnnData. For spliced and unspliced matrices, we listed a few options in the manual.

Then you can use brie-quant to perform quantification of splicing ratio and detect differential alternative splicing or differential momentum genes.

Type command line brie-count -h and brie-quant -h to see the full arguments.

References

Project details


Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Source Distribution

brie-2.1.0.tar.gz (487.4 kB view hashes)

Uploaded source

Supported by

AWS AWS Cloud computing Datadog Datadog Monitoring Facebook / Instagram Facebook / Instagram PSF Sponsor Fastly Fastly CDN Google Google Object Storage and Download Analytics Huawei Huawei PSF Sponsor Microsoft Microsoft PSF Sponsor NVIDIA NVIDIA PSF Sponsor Pingdom Pingdom Monitoring Salesforce Salesforce PSF Sponsor Sentry Sentry Error logging StatusPage StatusPage Status page