cellxgene

An interactive, performant explorer for single cell transcriptomics data.

cellxgene is an open-source experiment in how to bring powerful tools from modern web development to visualize and explore large single-cell transcriptomics datasets. Started in the context of the Human Cell Atlas Consortium, cellxgene hopes to both enable scientists to explore their data and to equip developers with scalable, reusable patterns and frameworks for visualizing large scientific datasets.

Features

• Visualization at scale: built with WebGL, React & Redux to handle visualization of at least 1 million cells.

• Interactive exploration: select, cross-filter, and compare subsets of your data with performant indexing and data handling.

• Flexible API: the cellxgene client-server model is designed to support a range of existing analysis packages for backend computational tasks (eg scanpy), integrated with client-side visualization via a REST API.

Getting Started

Requirements

• OS: OSX, Windows, Linux -- the developers are currently testing on OSX and Windows (via WSL using Ubuntu). It should work on other platforms but if you are using something different and need help, please let us know.
• python 3.6
• python3 tkinter
• npm
• Google Chrome

Clone project

git clone https://github.com/chanzuckerberg/cellxgene.git

Install client

cd cellxgene
./bin/build-client

To use with virtual env for python (optional, but recommended)

ENV_NAME=cellxgene
python3 -m venv ${ENV_NAME}
source ${ENV_NAME}/bin/activate

Install server

 pip install -e .

Run (with demo data)

cellxgene launch --title PBMC3K example-dataset/pbmc3k.h5ad

Help

cellxgene --help

For help with the scanpy engine

cellxgene scanpy --help

Using your own data

Scanpy

To prepare your data you will need to format your data into AnnData format using scanpy and calculate PCA and nearest neighbors and save in h5ad format.

1. Load data into scanpy

   • Ensure that obs's index is the cell names: print(data.obs_names) should show your cell indices. If it shows gene names, you may need to just call data.transpose().

2. Calculate PCA

   sc.pp.pca(data) ## sc is scanpy.api

3. Calculate nearest neighbors (depending on layout algorithm)

   # For umap layout algorithm, you need to use the "umap" method for neighbors
   sc.pp.neighbors(data, method="umap", metric="euclidean", use_rep="X_pca")
   
   # For tsne layout algorithm, you can use either "umap" or "gauss"; we recommend "gauss"
   sc.pp.neighbors(data, method="gauss", metric="euclidean", use_rep="X_pca")
   

4. Save file

   # cellxgene requires file to be named data.h5ad
   data.write("data.h5ad")
   

Contributing

We warmly welcome contributions from the community. Please submit any bug reports and feature requests through github issues. Please submit any direct contributions via a branch + pull request.

Inspiration and collaboration

We’ve been inspired by several other related efforts in this space, including the UCSC Cell Browswer, Cytoscape, Xena, ASAP, Gene Pattern, & many others; we hope to explore collaborations where useful.

Help/Contact

Have questions, suggestions, or comments? You can contact us by joining CZI Science Slack and posting in the #cellxgene channel. Please submit any feature requests or bugs as an issue in github. We'd love to hear from you!

Reuse

This project was started with the sole goal of empowering the scientific community to explore and understand their data. As such, we whole-heartedly encourage other scientific tool builders to adopt the patterns, tools, and code from this project, and reach out to us with ideas or questions using Github Issues or Pull Requests. All code is freely available for reuse under the MIT license.

Acknowledgements

cellxgene is inspired by many innovative projects. We would like to specifically thank:

• Alex Wolf for the demo dataset.
• Mike Bostock and the crossfilter team for API inspiration.