For the best experience,
browse Cellestial website on a computer.

Cellestial v0.53.1#

PyPI version License: Apache 2.0 Powered by Polars Built with Lets-Plot

The grammar of graphics for single-cell omics.

Installation#

pip install cellestial

uv add cellestial

poetry add cellestial

Examples#

Hover over the plot ‘geoms’ to see tooltips, or use toolbar above the plot for zooming and panning options.

Spatial#

Overlay categorical labels or gene expression on tissue coordinates.

Spatial categorical example
Show code
import cellestial as cl
import squidpy as sq
from lets_plot import *

data_spatial = cl.datasets.human_lymph_node()
data_hne = sq.datasets.visium_hne_adata()

s1 = cl.spatial(data_spatial, key="clusters")
s2 = cl.spatial(data_hne, key="leiden")
spatial_categorical = gggrid([s1, s2], ncol=2) + ggsize(1000, 400) + ggtb(size_zoomin=-1)
Spatial expression example
Show code
import cellestial as cl
import squidpy as sq
from lets_plot import *

data_spatial = cl.datasets.human_lymph_node()
data_hne = sq.datasets.visium_hne_adata()

s3 = cl.spatial(data_spatial, key="MS4A1") + scale_color_viridis(option="inferno")
s4 = cl.spatial(data_hne, key="Mef2c") + scale_color_viridis(option="inferno")
spatial_numeric = gggrid([s3, s4], ncol=2) + ggsize(1000, 400) + ggtb(size_zoomin=-1)

Dimensionality Reduction#

Cellestial offers dimensionality reduction plots such as UMAP and t-SNE, which accept a key for categorical and numeric columns.

Also, cellestial-specific layers, providing utilities, could be added to dimensionality reduction plots.

Dimensionality reduction example
Show code
import cellestial as cl
from lets_plot import *

data = cl.datasets.pbmc3k()

dims = gggrid(
    [
        cl.umap(data, key="cell_type_lvl1", axis_type="arrow", size=1.5, legend_ondata=True)
        + scale_color_hue(),
        cl.tsne(
            data,
            key="NEAT1",
            tooltips=["NEAT1", "cell_type_lvl1"],
            axis_type="arrow",
            size=2,
            color_high="#219B9D",
            color_low="#f6f6f6",
        ),
    ]
) + ggtb(size_zoomin=-1)
Dimensionality reduction layers example
Show code
import cellestial as cl
from lets_plot import *

data = cl.datasets.pbmc3k()
data_velocity = cl.datasets.pancreas()

outlined = (
    cl.umap(data, key="cell_type_lvl1", axis_type="arrow", size=1.5, legend_ondata=True)
    + scale_color_hue()
    + cl.cluster_outlines(groups=["Lymphocytes", "B Cells"])
)
streamed = (
    cl.umap(
        data_velocity,
        key="clusters_coarse",
        axis_type="arrow",
        size=4,
        alpha=0.4,
        legend_ondata=True,
        ondata_color="black",
    )
    + cl.stream()
)
layers = gggrid([outlined, streamed]) + ggtb(size_zoomin=-1)

Heatmaps#

Cellestial offers heatmap variants, with built-in dendrogram features.

Heatmap example
Show code
import cellestial as cl
from lets_plot import *

data = cl.datasets.pbmc3k()

markers = [
    "PSAP",
    "LYZ",
    "CST3",
    "CD79A",
    "CD79B",
    "IL7R",
    "CD3D",
    "CD3E",
    "CD4",
    "CD8A",
    "CD8B",
    "NKG7",
    "GNLY",
    "KLRD1",
    "HLA-DRA",
    "FCER1A",
]
htmp = (
    cl.heatmap(
        data,
        group_by="cell_type_lvl1",
        keys=markers,
        geom="raster",
        group_lines_size=0.5,
        group_lines_color="white",
        dendrogram=True,
        dendrogram_size=1,
        group_bars_labels=True,
        group_bars=True,
    )
    + scale_fill_viridis()
)
Matrixplot example
Show code
import cellestial as cl
from lets_plot import *

data = cl.datasets.pbmc3k()

markers = [
    "PSAP",
    "LYZ",
    "CST3",
    "CD79A",
    "CD79B",
    "IL7R",
    "CD3D",
    "CD3E",
    "CD4",
    "CD8A",
    "CD8B",
    "NKG7",
    "GNLY",
    "KLRD1",
    "HLA-DRA",
    "FCER1A",
]
mtrx = cl.matrixplot(
    data,
    group_by="leiden",
    keys=markers,
    group_lines_size=0.5,
    group_lines_color="white",
    dendrogram=True,
    dendrogram_size=1,
) + scale_fill_viridis(option="inferno", begin=0.1)
Stacked violin example
Show code
import cellestial as cl
from lets_plot import *

data = cl.datasets.pbmc3k()

markers = [
    "PSAP",
    "LYZ",
    "CST3",
    "CD79A",
    "CD79B",
    "IL7R",
    "CD3D",
    "CD3E",
    "CD4",
    "CD8A",
    "CD8B",
    "NKG7",
    "GNLY",
    "KLRD1",
    "HLA-DRA",
    "FCER1A",
]
stck = cl.stacked_violin(
    data,
    group_by="leiden_res_0.50",
    keys=markers,
    dendrogram=True,
    dendrogram_size=1,
)
Dotplot example
Show code
import cellestial as cl
from lets_plot import *

data = cl.datasets.pbmc3k()

markers = [
    "PSAP",
    "LYZ",
    "CST3",
    "CD79A",
    "CD79B",
    "IL7R",
    "CD3D",
    "CD3E",
    "CD4",
    "CD8A",
    "CD8B",
    "NKG7",
    "GNLY",
    "KLRD1",
    "HLA-DRA",
    "FCER1A",
]
dtplt = cl.dotplot(
    data,
    group_by="leiden_res_0.50",
    keys=markers,
    dendrogram=True,
    dendrogram_size=1,
)

Distribution#

Cellestial comes with distribution (boxplot and violin) plots with built-in statistical analysis.

Boxplot example
Show code
import cellestial as cl
from lets_plot import *

data = cl.datasets.pbmc3k()

bx = (
    cl.boxplot(
        data,
        key="CD3D",
        fill="cell_type_lvl1",
        threshold=0.1,
    )
    + scale_fill_hue()
)
Boxplot with brackets example
Show code
import cellestial as cl
from lets_plot import *

data = cl.datasets.pbmc3k()

bx = (
    cl.boxplot(
        data,
        key="CD3D",
        fill="cell_type_lvl1",
        threshold=0.1,
    )
    + scale_fill_hue()
)
bx_brckt = bx + cl.bracket(
    y_padding=0.2,
    label="pvalue",
    prefix_style="<",
    prefix="p",
    comparisons=[
        ("Lymphocytes", "Monocytes"),
        ("Monocytes", "Erythroid"),
        ("Monocytes", "B Cells"),
    ],
)
Violin example
Show code
import cellestial as cl
from lets_plot import *

data = cl.datasets.pbmc3k()

vln = (
    cl.violin(
        data,
        key="CD3D",
        fill="cell_type_lvl1",
        threshold=0.1,
    )
    + scale_fill_viridis()
)
Violin with brackets example
Show code
import cellestial as cl
from lets_plot import *

data = cl.datasets.pbmc3k()

vln = (
    cl.violin(
        data,
        key="CD3D",
        fill="cell_type_lvl1",
        threshold=0.1,
    )
    + scale_fill_viridis()
)
vln_brckt = vln + cl.bracket(
    y_padding=0.2,
    label="pvalue",
    prefix_style="<",
    prefix="p",
    comparisons=[
        ("Lymphocytes", "Monocytes"),
        ("Monocytes", "Erythroid"),
        ("Monocytes", "B Cells"),
    ],
)

Exploratory#

Cellestial ships with quick exploratory plots.

Ridge example
Show code
import cellestial as cl
from lets_plot import *

data = cl.datasets.pbmc3k()

rdg = (
    cl.ridge(
        data,
        key="B2M",
        alpha=0.6,
        group_by="cell_type_lvl1",
    )
    + scale_fill_hue()
)
Highest expressed genes example
Show code
import cellestial as cl
from lets_plot import *

data = cl.datasets.pbmc3k()

hexp = cl.highest_expressed_genes(data, n=20) + scale_fill_viridis()

Documentation#

About Lets-Plot#

Cellestial is built on top of a powerful Python library, Lets-Plot. It is the best Python implementation of ggplot2 with additional features such as tooltips and zooming and panning. Lets-Plot API