Learn R Programming

Bisque

An R toolkit for accurate and efficient estimation of cell composition ('decomposition') from bulk expression data with single-cell information.

Bisque provides two modes of operation:

Reference-based decomposition

This method utilizes single-cell data to decompose bulk expression. We assume that both single-cell and bulk counts are measured from the same tissue. Specifically, the cell composition of the labeled single-cell data should match the expected physiological composition. While we don't explicitly require matched samples, we expect having samples with both single-cell and bulk expression measured will provide more accurate results.

Marker-based decomposition

This method utilizes marker genes alone to decompose bulk expression when a reference profile is not available. Single-cell data is not explicitly required but can be used to identify these marker genes. This method captures relative abundances of a cell type across individuals. Note that these abundances are not proportions, so they cannot be compared between different cell types.

Installation

The Bisque R package is available on CRAN

install.packages("BisqueRNA")

as well as Bioconda

conda install r-bisquerna

The package can also be installed from the GitHub repository

devtools::install_github("cozygene/bisque")

Getting Started

You can load Bisque as follows:

library(BisqueRNA)

The two modes of operation described above are called as follows:

res <- BisqueRNA::ReferenceBasedDecomposition(bulk.eset, sc.eset, markers)
res <- BisqueRNA::MarkerBasedDecomposition(bulk.eset, markers)

Each method returns a list of results with estimated cell proportions/abundances stored in res$bulk.props.

To see examples of these methods on simulated data, check out the vignette:

browseVignettes("BisqueRNA")

Copy Link

Version

Install

install.packages('BisqueRNA')

Monthly Downloads

451

Version

1.0.5

License

GPL-3

Maintainer

Brandon Jew

Last Published

May 23rd, 2021

Functions in BisqueRNA (1.0.5)

GetNumGenesWeighted

Get number of genes to use with weighted PCA
FilterUnexpressedGenes

Remove genes in Expression Set with zero expression in all samples
FilterZeroVarianceGenes

Remove genes in Expression Set with zero variance across samples
ReferenceBasedDecomposition

Performs reference-based decomposition of bulk expression using single-cell data
SemisupervisedTransformBulk

Transforms bulk expression of a gene using only single-cell data
SupervisedTransformBulk

Transforms bulk expression of a gene given overlapping data
SimulateData

Simulate data for decomposition illustration
CorTri

Correlate columns of data frame
MarkerBasedDecomposition

Performs marker-based decomposition of bulk expression using marker genes
GetUniqueMarkers

Get unique markers present in only 1 cell type
CalculateSCCellProportions

Calculate cell proportions based on single-cell data
CountsToCPM

Convert counts data in Expression Set to counts per million (CPM)
EstimatePCACellTypeProportions

Estimate cell type proportions using first PC of expression matrix
GetOverlappingSamples

Find overlapping samples in single-cell and bulk data
GetOverlappingGenes

Find overlapping genes in single-cell data, bulk data, and marker genes
GenerateSCReference

Generate reference profile for cell types identified in single-cell data
GetCTP

Return cell type proportions from bulk
GetNumGenes

Get number of genes to use with no weighted information
SeuratToExpressionSet

Converts Seurat object to Expression Set
SimulateBarcode

Simulate barcode for decomposition illustration