Skip to main content

Single-cell multi-omic datasets, in which multiple molecular modalities are profiled within the same cell, provide a unique opportunity to discover the interplay between cellular epigenomic and transcriptomic changes. To realize this potential, we developed MultiVelo, a mechanistic model of gene expression that extends the popular RNA velocity framework by incorporating epigenomic data.

MultiVelo uses a probabilistic latent variable model to estimate the switch time and rate parameters of gene regulation, providing a quantitative summary of the temporal relationship between epigenomic and transcriptomic changes. Fitting MultiVelo on single-cell multi-omic datasets revealed two distinct mechanisms of regulation by chromatin accessibility, quantified the degree of concordance or discordance between transcriptomic and epigenomic states within each cell, and inferred the lengths of time lags between transcriptomic and epigenomic changes.

Installation

Install through PyPI:

pip install multivelo

The package is also available on Bioconda. Install with:

conda install -c bioconda multivelo or mamba install -c bioconda multivelo

Documentation

We have a ReadTheDocs page.

Tutorial

New: we have added Jupyter notebooks showing how to reproduce the main figure panels, along with all required processed data files. These can be found under the Examples folder in this repository or on our ReadTheDocs page.

A tutorial showing how to run MultiVelo can be found here: (jupyter notebook)

The tutorial uses the embryonic E18 mouse brain from 10X Multiome as an example. CellRanger output files can be downloaded from 10X website. Crucially, the filtered feature barcode matrix folder, ATAC peak annotations TSV, and the feature linkage BEDPE file in the secondary analysis outputs folder will be needed in this demo.

You can download the processed data that we used for this analysis if you want to run the example yourself. Unspliced and spliced counts, as well as cell type annotations can be downloaded from the MultiVelo GitHub page. We provide the cell annotations for this dataset in “cell_annotations.tsv”. We also provide the nearest neighbor graph used to smooth chromatin accessibility values in the GitHub folder “seurat_wnn”, which contains a zip file of three files: “nn_cells.txt”, “nn_dist.txt”, and “nn_idx.txt”. Please unzip the archive after downloading. The R script used to generate these files can also be found in the same folder.

Citation

Li, C., Virgilio, M.C., Collins, K.L. & Welch J.D. Multi-omic single-cell velocity models epigenome–transcriptome interactions and improves cell fate prediction. Nat Biotechnol 41, 387-398 (2023).

Metadata

Release files for multivelo 0.1.5

For a detailed explanation of source distributions (sdists) and built distributions (wheels), please see the package formats documentation.

Built distribution (wheel)

Table of built distributions (wheels) for multivelo 0.1.5
File Interpreter ABI Platform
multivelo-0.1.5-py3-none-any.whl Python 3 none any Details

Release files / multivelo-0.1.5-py3-none-any.whl

Download URL multivelo-0.1.5-py3-none-any.whl
Size 527.5 kB
Tags Python 3
SHA-256 checksum
How to use checksums
22e4ad673c713819f7e742375badc852611176e55f07b90f9fd3cdd1059bcdda
BLAKE2b-256 checksum
How to use checksums
5a84aa290674fa57e87a3884d25b2d8ee9cde2182a7f89a605d01cb7517d3395
Upload date
Uploaded using Trusted Publishing?
What is trusted publishing?
No
Uploaded via twine/5.1.1 CPython/3.11.4

Release history Release notifications | RSS feed

This release

0.1.5 This release

1 release file

0.1.4

1 release file

0.1.3

2 release files

0.1.2

2 release files

0.1.1

2 release files

0.1.0

2 release files

Anthropic, PBC Visionary sponsor Bloomberg Visionary sponsor Hudson River Trading Visionary sponsor Meta Visionary sponsor NVIDIA Visionary sponsor Microsoft Sustainability sponsor Depot Continuous Integration AWS Cloud computing and Security Sponsor Datadog Monitoring Fastly CDN Google Download Analytics Sentry Error logging StatusPage Status page