Skip to content

colabobio/scsketch

Repository files navigation

scSketch project

scSketch is an interactive exploration tool of single-cell embeddings (UMAP, tSNE, etc.) for Python notebooks. It is based on the jupyter-scatter widget by Fritz Lekschas and it reimplements some features of the earlier SciViewer visualizer.

scSketch allows users to explore embeddings by selecting linear directions and cells of interest to see the genes/proteins that are changing the most among the selected cells and direction. It then allows users to see what biological pathways the genes/proteins are a part of (by connecting to Reactome database) and provides users with the diagram of biological pathways where users can see their gene/protein of interest and what other genes/proteins/molecules they interact with. The widget can handle millions of points and supports view linking.

scSketch directional analysis demo

Directional sketch → compute directional analysis → click a gene for pathway context.

scSketch differential expression demo

Select two groups → compute differential expression → browse results.

Usage

Quick Start

The easiest way to try scSketch is with the built-in demo (no installation required):

uvx scsketch demo

This single command will automatically install scSketch and all dependencies in an ephemeral environment, then launch the demo notebook. It requires uv, which is a fast Python package manager.

Alternatively, if you've cloned the repository, you can run the demo notebook directly with juv:

git clone https://github.com/colabobio/scsketch.git
cd scsketch
uvx juv run demo.ipynb

Then use in any notebook:

from scsketch import ScSketch

sketch = ScSketch(adata=...)
sketch.show()

See the demo notebook for more details

Use scSketch in your own notebook

scSketch requires Python 3.10 or later. Before installing, make sure your environment meets this requirement:

python --version   # should print Python 3.10.x or higher

If you are using conda and need to create or upgrade an environment:

# Create a new environment with the right Python version
conda create -n my-env python=3.10
conda activate my-env

Then install scSketch into the environment backing your Jupyter kernel:

pip install scsketch

Optional: Numba-accelerated kernels ([fast] extra)

Differential expression computations can be significantly accelerated by Numba. Numba is an optional dependency — scSketch works without it, but automatically uses it when available. To install with Numba:

pip install "scsketch[fast]"

Note: If you install scSketch into a conda environment, make sure you also have JupyterLab installed in that same environment so the kernel picks up the right packages:

pip install jupyterlab
jupyter lab

Then in a notebook:

import scanpy as sc
from scsketch import ScSketch

adata = sc.read_h5ad("my_data.h5ad")

# scSketch currently reads coordinates from `adata.obsm["X_umap"]`.
# If you have a different embedding (e.g. tSNE), you can copy it into `X_umap`:
# adata.obsm["X_umap"] = adata.obsm["X_tsne"]

sketch = ScSketch(
    adata=adata,
    metadata_cols=["louvain"],   # optional: columns in `adata.obs` for coloring
    color_by_default="louvain",  # optional: which metadata to color by initially
)

# If this isn't the last line in the cell, use: `from IPython.display import display; display(sketch.show())`
sketch.show()

Directional Search: keep brush selections roughly linear

Directional Search reduces your selection to a 1D “along-the-sketch” axis by projecting cells onto a single direction vector. If your brush selection is very curved, loops back, or spans multiple branches/blobs, that 1D projection can mix multiple directions of variation and produce hard-to-interpret results.

Practical tips:

  • Sketch along one clear gradient at a time (a selection closer to a straight line works best).
  • If the trajectory bends, split it into multiple shorter selections and compare results.

Gene IDs vs gene symbols

scSketch currently uses adata.var_names as the gene identifier for display/search. If your AnnData uses Ensembl IDs (e.g. ENSG...) in var_names, you will see Ensembl IDs in the UI.

If you have gene symbols in a column like adata.var["gene_symbols"], you can create a visualization-only copy that displays symbols:

adata_view = adata.copy()
adata_view.var["ensembl_id"] = adata_view.var_names
adata_view.var_names = adata_view.var["gene_symbols"].astype(str)
adata_view.var_names_make_unique()

sketch = ScSketch(adata=adata_view, metadata_cols=["louvain"], color_by_default="louvain")
sketch.show()

Using scSketch with uv / uvx

Here are the recommended ways to integrate scSketch when using uv / uvx.

Add to an existing uv project:

uv add scsketch          # base install
uv add "scsketch[fast]"  # with Numba-accelerated kernels (recommended)
uv run jupyter lab

Standalone notebook with juv:

juv runs notebooks in isolated environments defined by inline metadata — no pyproject.toml needed.

# Add scsketch to an existing notebook (writes inline dependency metadata)
uvx juv add my_notebook.ipynb "scsketch[fast]"

# Launch it in an auto-provisioned environment
uvx juv run my_notebook.ipynb

The notebook becomes fully self-contained and reproducible: anyone with juv can run it without any prior setup.

Try the built-in demo without installing anything:

uvx scsketch demo

This uses uvx to run scSketch ephemerally — nothing is permanently installed in your environment.

Running the original notebook with juv

To run the original inline notebook, first install juv and then call:

juv run demo.ipynb

Development

This project uses uv for development and dependency management.

Setup

  1. Clone the repository:

    git clone https://github.com/colabobio/scsketch.git
    cd scsketch
  2. Sync environment (installs dependencies):

    uv sync

Development Workflow

You can run commands inside the project's environment using uv run.

  • Launch Jupyter Lab for testing:

    uv run jupyter lab

    Open debug.ipynb to test changes.

  • Hot-reloading JS/CSS: The debug.ipynb notebook is pre-configured with ANYWIDGET_HMR=1. Any changes you save to files in src/scsketch/static/ will legally update the widget in your browser without reloading the page.

  • Linting (Ruff):

    uv run ruff check .
  • Testing: Run the unit test suite with pytest:

    uv run pytest tests/

Optional: Manual Activation

If you prefer to activate the environment in your shell:

source .venv/bin/activate
# Now you can use `jupyter`, `python`, `pip` directly
jupyter lab

Or with editable installs:

python -m venv .venv
source .venv/bin/activate
python -m pip install -e ".[dev,fast]"
jupyter lab demo.ipynb

Debugging with VS Code

To debug the Python side of your widgets step-by-step:

  1. Debugging Jupyter Notebooks:

    • Open debug.ipynb in VS Code.
    • Click the "Select Kernel" button at the top right and select your project environment (e.g., .venv or the one created by uv).
    • You can set breakpoints directly in the notebook cells or in the python files (src/scsketch/widgets/*.py).
    • To debug code in external files (like the widgets), use the "Debug Cell" option (often found in the dropdown menu next to the run button of a cell).
  2. Attaching to a running kernel: If you prefer using jupyter lab in your browser but want to debug Python code in VS Code:

    • Run jupyter lab (e.g., uv run jupyter lab).
    • Add this code to a cell at the beginning of your notebook:
      import debugpy
      debugpy.listen(5678)
      print("Waiting for debugger attach...")
      debugpy.wait_for_client()
      print("Debugger attached")
    • In VS Code, go to the Run and Debug view (Ctrl+Shift+D).
    • Select "Python: Attach to Local Process" from the dropdown and click the Play button.
    • VS Code will attach to your running kernel, and you can now use breakpoints in your local Python files.

Debugging Frontend (JavaScript)

Since the widgets run in the web browser (or VS Code's webview), you need to use browser developer tools to debug the JavaScript code (src/scsketch/static/*.js).

  1. Run the widget: Open debug.ipynb and run the cell that displays the widget.
  2. Open Developer Tools:
    • In Browser (Jupyter Lab): Right-click anywhere on the page > Inspect.
    • In VS Code: Open the command palette (Cmd+Shift+P) and run "Developer: Open Webview Developer Tools".
  3. Find your source:
    • Go to the Sources tab in the developer tools.
    • Use Cmd+P (Mac) or Ctrl+P (Windows/Linux) to search for your file (e.g., correlation_table.js).
    • Note: Because of how modules are loaded, the file path might look like localhost:xyz/.../correlation_table.js.
  4. Set Breakpoints: Click on the line number in the JS file to set a breakpoint.
  5. Trigger the code: Interact with the widget in the notebook. The debugger will pause on your breakpoint, allowing you to inspect variables and step through the code.

Publish a New Version

To bump the version use one of the following commands:

  1. uvx bump-my-version bump minor (e.g., v0.1.0 → v0.2.0)
  2. uvx bump-my-version bump patch (e.g., v0.1.0 → v0.1.1)
  3. uvx bump-my-version bump major (e.g., v0.1.0 → v1.0.0)

Afterward do git push --follow-tags. Github actions will handle the rest.

About

No description, website, or topics provided.

Resources

Stars

Watchers

Forks

Packages

 
 
 

Contributors