Visualization of single cell RNA-seq¶
In this part, we will introduce the tutorial of special plot of omicverse
.
import omicverse as ov
import scanpy as sc
#import scvelo as scv
ov.plot_set()
____ _ _ __ / __ \____ ___ (_)___| | / /__ _____________ / / / / __ `__ \/ / ___/ | / / _ \/ ___/ ___/ _ \ / /_/ / / / / / / / /__ | |/ / __/ / (__ ) __/ \____/_/ /_/ /_/_/\___/ |___/\___/_/ /____/\___/ Version: 1.5.10, Tutorials: https://omicverse.readthedocs.io/
We utilized single-cell RNA-seq data (GEO accession: GSE95753) obtained from the dentate gyrus of the hippocampus in mouse.
adata = ov.read('data/DentateGyrus/10X43_1.h5ad')
adata
AnnData object with n_obs × n_vars = 2930 × 13913 obs: 'clusters', 'age(days)', 'clusters_enlarged' uns: 'clusters_colors' obsm: 'X_umap' layers: 'ambiguous', 'spliced', 'unspliced'
Optimizing color mapping¶
Visualizing spatially resolved biological data with appropriate color mapping can significantly facilitate the exploration of underlying patterns and heterogeneity. Spaco (spatial colorization) provides a spatially constrained approach that generates discriminate color assignments for visualizing single-cell spatial data in various scenarios.
Jing Z, Zhu Q, Li L, Xie Y, Wu X, Fang Q, et al. Spaco: A comprehensive tool for coloring spatial data at single-cell resolution. Patterns. 2024;100915
Function: ov.pl.optim_palette
:
- adata: the datasets of scRNA-seq
- basis: he position on the plane should be set to
X_spatial
in spatial RNA-seq,X_umap
,X_tsne
,X_mde
in scRNA-seq and should not be set toX_pca
- colors: Specify the colour to be optimised, which should be for one of the columns in adata.obs, noting that it should have the colour first, and that we can use ov.pl.embedding to colour the cell types. If there is no colour then colour blind optimisation colour will be used.
- palette: You can also re-specify the colour bar that needs to be drawn, just set
palette=['#FFFFFF','#000000']
, we have preparedov.pl.red_color
,ov.pl.blue_color
,ov.pl.green_color
,ov.pl.orange_color
, by default.
optim_palette=ov.pl.optim_palette(adata,basis='X_umap',colors='clusters')
|-----> Calculating cluster distance graph... |-----------> Calculating cell neighborhood... |-----------> Filtering out neighborhood outliers... |-----------> Calculating cluster interlacement score... |-----------> Constructing cluster interlacement graph... |-----> Calculating color distance graph... |-----------> Calculating color perceptual distance... |-----------> Constructing color distance graph... |-----------> Difference of the most similar pair in the palette is 53.88 |-----> Optimizing color mapping...
import matplotlib.pyplot as plt
fig,ax=plt.subplots(figsize = (4,4))
ov.pl.embedding(adata,
basis='X_umap',
color='clusters',
frameon='small',
show=False,
palette=optim_palette,
ax=ax,)
plt.title('Cell Type of DentateGyrus',fontsize=15)
Stacked histogram of cell proportions¶
This is a graph that appears widely in various CNS-level journals, and is limited to the fact that scanpy
does not have a proper way of plotting it, and we provide ov.pl.cellproportion
for plotting it here.
Function: ov.pl.cellproportion
:
- adata: the datasets of scRNA-seq
- celltype_clusters: Specify the colour to plot, which should be for one of the columns in adata.obs, noting that it should have the colour first, and that we can use ov.pl.embedding to colour the cell types. If there is no colour then colour blind optimisation colour will be used.
- groupby: The group variable for the different groups of cell types we need to display, in this case we are displaying different ages, so we set it to
age(days)
- groupby_li: If there are too many groups, we can also select the ones we are interested in plotting, here we use groupby_li to plot the groups
- figsize: If we specify axes, then this variable can be left empty
- legend: Whether to show a legend
import matplotlib.pyplot as plt
fig,ax=plt.subplots(figsize = (1,4))
ov.pl.cellproportion(adata=adata,celltype_clusters='clusters',
groupby='age(days)',legend=True,ax=ax)
fig,ax=plt.subplots(figsize = (2,2))
ov.pl.cellproportion(adata=adata,celltype_clusters='age(days)',
groupby='clusters',groupby_li=['nIPC','Granule immature','Granule mature'],
legend=True,ax=ax)
A collection of some interesting embedded plot¶
Our first presentation is an embedding map with the number and proportion of cell types. This graph visualises the low-dimensional representation of cells in addition to the number of cell proportions, etc. It should be noted that the cell proportions plotted on the left side may be distorted when there are too many cell types, and we would be grateful if anyone would be interested in fixing this bug.
Function: ov.pl.embedding_celltype
:
- adata: the datasets of scRNA-seq
- figsize: Note that we don't usually provide the ax parameter for combinatorial graphs, this is due to the fact that combinatorial graphs are made up of multiple axes, so the figsize parameter is more important, here we set it to
figsize=(7,4)
. - basis: he position on the plane should be set to
X_spatial
in spatial RNA-seq,X_umap
,X_tsne
,X_mde
in scRNA-seq and should not be set toX_pca
- celltype_key: Specify the colour to be optimised, which should be for one of the columns in adata.obs, noting that it should have the colour first, and that we can use ov.pl.embedding to colour the cell types. If there is no colour then colour blind optimisation colour will be used.
- title: Note that the space entered in title is used to control the position.
- celltype_range: Since our number of cell types is different in each data, we want to have the flexibility to control where the cell scale plot is drawn, here we set it to
(1,10)
. You can also tweak the observations yourself to find the parameter that best suits your data. - embedding_range: As with the positional parameters of the cell types, they need to be adjusted several times on their own for optimal results.
ov.pl.embedding_celltype(adata,figsize=(7,4),basis='X_umap',
celltype_key='clusters',
title=' Cell type',
celltype_range=(1,10),
embedding_range=(4,10),)
Sometimes we want to be able to circle a certain type of cell that we are interested in, and here we use convex polygons to achieve this, while the shape of the convex polygons may be optimised in future versions.
Function: ov.pl.ConvexHull
:
- adata: the datasets of scRNA-seq
- basis: he position on the plane should be set to
X_spatial
in spatial RNA-seq,X_umap
,X_tsne
,X_mde
in scRNA-seq and should not be set toX_pca
- cluster_key: Specify the celltype to be optimised, which should be for one of the columns in adata.obs, noting that it should have the colour first, and that we can use ov.pl.embedding to colour the cell types. If there is no colour then colour blind optimisation colour will be used.
- hull_cluster: the target celltype to be circled.
import matplotlib.pyplot as plt
fig,ax=plt.subplots(figsize = (4,4))
ov.pl.embedding(adata,
basis='X_umap',
color=['clusters'],
frameon='small',
show=False,
ax=ax)
ov.pl.ConvexHull(adata,
basis='X_umap',
cluster_key='clusters',
hull_cluster='Granule mature',
ax=ax)
Besides, if you don't want to plot convexhull, you can plot the contour instead.
Function: ov.pl.contour
:
- adata: the datasets of scRNA-seq
- basis: he position on the plane should be set to
X_spatial
in spatial RNA-seq,X_umap
,X_tsne
,X_mde
in scRNA-seq and should not be set toX_pca
- groupby: Specify the celltype to be optimised, which should be for one of the columns in adata.obs, noting that it should have the colour first, and that we can use ov.pl.embedding to colour the cell types. If there is no colour then colour blind optimisation colour will be used.
- clusters: the target celltype to be circled.
- colors: the color of the contour
- linestyles: the linestyles of the contour
- **kwargs: more kwargs could be found from
plt.contour
import matplotlib.pyplot as plt
fig,ax=plt.subplots(figsize = (4,4))
ov.pl.embedding(adata,
basis='X_umap',
color=['clusters'],
frameon='small',
show=False,
ax=ax)
ov.pl.contour(ax=ax,adata=adata,groupby='clusters',clusters=['Granule immature','Granule mature'],
basis='X_umap',contour_threshold=0.1,colors='#000000',
linestyles='dashed',)
In scanpy's default embedding
plotting function, when we set legend=True, legend masking may occur. To solve this problem, we introduced ov.pl.embedding_adjust
in omicverse to automatically adjust the position of the legend.
Function: ov.pl.embedding_adjust
:
- adata: the datasets of scRNA-seq
- basis: he position on the plane should be set to
X_spatial
in spatial RNA-seq,X_umap
,X_tsne
,X_mde
in scRNA-seq and should not be set toX_pca
- groupby: Specify the celltype to be optimised, which should be for one of the columns in adata.obs, noting that it should have the colour first, and that we can use ov.pl.embedding to colour the cell types. If there is no colour then colour blind optimisation colour will be used.
- exclude: We can specify which cell types are not to be plotted, in this case we set it to
OL
- adjust_kwargs: We can manually specify the parameters of adjustText, the specific parameters see the documentation of adjustText, it should be noted that we have to use dict to specify the parameters here.
- text_kwargs: We can also specify the font colour manually by specifying the text_kwargs parameter
from matplotlib import patheffects
import matplotlib.pyplot as plt
fig, ax = plt.subplots(figsize=(4,4))
ov.pl.embedding(adata,
basis='X_umap',
color=['clusters'],
show=False, legend_loc=None, add_outline=False,
frameon='small',legend_fontoutline=2,ax=ax
)
ov.pl.embedding_adjust(
adata,
groupby='clusters',
exclude=("OL",),
basis='X_umap',
ax=ax,
adjust_kwargs=dict(arrowprops=dict(arrowstyle='-', color='black')),
text_kwargs=dict(fontsize=12 ,weight='bold',
path_effects=[patheffects.withStroke(linewidth=2, foreground='w')] ),
)
Sometimes we are interested in the distribution density of a certain class of cell types in a categorical variable, which is cumbersome to plot in the scanpy
implementation, so we have simplified the implementation in omicverse and ensured the same plotting.
Function: ov.pl.embedding_density
:
- adata: the datasets of scRNA-seq
- basis: he position on the plane should be set to
X_spatial
in spatial RNA-seq,X_umap
,X_tsne
,X_mde
in scRNA-seq and should not be set toX_pca
- groupby: Specify the celltype to be optimised, which should be for one of the columns in adata.obs, noting that it should have the colour first, and that we can use ov.pl.embedding to colour the cell types. If there is no colour then colour blind optimisation colour will be used.
- target_clusters: We can specify which cell types are to be plotted, in this case we set it to
Granule mature
- kwargs: other parameter can be found in
scanpy.pl.embedding
ov.pl.embedding_density(adata,
basis='X_umap',
groupby='clusters',
target_clusters='Granule mature',
frameon='small',
show=False,cmap='RdBu_r',alpha=0.8)
Bar graph with overlapping dots (Bar-dot) plot¶
In biological research, bardotplot plots are the most common class of graphs we use, but unfortunately, there is no direct implementation of plotting functions in either matplotlib, seaborn or scanpy. To compensate for this, we implement bardotplot plotting in omicverse and provide manual addition of p-values (it should be noted that manual addition refers to the manual addition of p-values for model fitting rather than making up p-values yourself).
ov.single.geneset_aucell(adata,
geneset_name='Sox',
geneset=['Sox17', 'Sox4', 'Sox7', 'Sox18', 'Sox5'])
ctxcore have been install version: 0.2.0
Function: ov.pl.embedding_density
:
- adata: the datasets of scRNA-seq
- groupby: Specify the celltype to be optimised, which should be for one of the columns in adata.obs, noting that it should have the colour first, and that we can use ov.pl.embedding to colour the cell types. If there is no colour then colour blind optimisation colour will be used.
- color: The size of the variable to be plotted, which can be a gene, stored in adata.var, or a cell value, stored in adata.obs.
- bar_kwargs: We provide the parameters of the barplot for input, see the matplotlib documentation for more details
- scatter_kwargs: We also provide the parameters of the scatter for input, see the matplotlib documentation for more details
Function: ov.pl.add_palue
:
- ax: the axes of bardotplot
- line_x1: The left side of the p-value line to be plotted
- line_x2: The right side of the p-value line to be plotted|
- line_y: The height of the p-value line to be plotted
- text_y: How much above the p-value line is plotted text
- text: the text of p-value, you can set
***
to insteadp<0.001
- fontsize: the fontsize of text
- fontcolor: the color of text
- horizontalalignment: the location of text
fig, ax = plt.subplots(figsize=(6,2))
ov.pl.bardotplot(adata,groupby='clusters',color='Sox_aucell',figsize=(6,2),
ax=ax,
ylabel='Expression',
bar_kwargs={'alpha':0.5,'linewidth':2,'width':0.6,'capsize':4},
scatter_kwargs={'alpha':0.8,'s':10,'marker':'o'})
ov.pl.add_palue(ax,line_x1=3,line_x2=4,line_y=0.1,
text_y=0.02,
text='$p={}$'.format(round(0.001,3)),
fontsize=11,fontcolor='#000000',
horizontalalignment='center',)
fig, ax = plt.subplots(figsize=(6,2))
ov.pl.bardotplot(adata,groupby='clusters',color='Sox17',figsize=(6,2),
ax=ax,
ylabel='Expression',xlabel='Cell Type',
bar_kwargs={'alpha':0.5,'linewidth':2,'width':0.6,'capsize':4},
scatter_kwargs={'alpha':0.8,'s':10,'marker':'o'})
ov.pl.add_palue(ax,line_x1=3,line_x2=4,line_y=2,
text_y=0.2,
text='$p={}$'.format(round(0.001,3)),
fontsize=11,fontcolor='#000000',
horizontalalignment='center',)
Boxplot with jitter points¶
A box plot, also known as a box-and-whisker plot, is a graphical representation used to display the distribution and summary statistics of a dataset. It provides a concise and visual way to understand the central tendency, spread, and potential outliers in the data.
Function: ov.pl.single_group_boxplot
:
- adata (AnnData object): The data object containing the information for plotting.
- groupby (str): The variable used for grouping the data
- color (str): The variable used for coloring the data points.
- type_color_dict (dict): A dictionary mapping group categories to specific colors.
- scatter_kwargs (dict): Additional keyword arguments for customizing the scatter plot.
- ax (matplotlib.axes.Axes): A pre-existing axes object for plotting (optional). (optional).(optional).
import pandas as pd
import seaborn as sns
#sns.set_style('white')
ov.pl.single_group_boxplot(adata,groupby='clusters',
color='Sox_aucell',
type_color_dict=dict(zip(pd.Categorical(adata.obs['clusters']).categories, adata.uns['clusters_colors'])),
x_ticks_plot=True,
figsize=(5,2),
kruskal_test=True,
ylabel='Sox_aucell',
legend_plot=False,
bbox_to_anchor=(1,1),
title='Expression',
scatter_kwargs={'alpha':0.8,'s':10,'marker':'o'},
point_number=15,
sort=False,
save=False,
)
plt.grid(False)
plt.xticks(rotation=90,fontsize=12)
(array([ 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13]), [Text(0, 0, 'Microglia'), Text(1, 0, 'OPC'), Text(2, 0, 'Granule immature'), Text(3, 0, 'GABA'), Text(4, 0, 'Astrocytes'), Text(5, 0, 'nIPC'), Text(6, 0, 'Cajal Retzius'), Text(7, 0, 'Endothelial'), Text(8, 0, 'Cck-Tox'), Text(9, 0, 'Mossy'), Text(10, 0, 'Granule mature'), Text(11, 0, 'OL'), Text(12, 0, 'Radial Glia-like'), Text(13, 0, 'Neuroblast')])
Complexheatmap¶
A complex heatmap, also known as a clustered heatmap, is a data visualization technique used to represent complex relationships and patterns in multivariate data. It combines several elements, including clustering, color mapping, and hierarchical organization, to provide a comprehensive view of data across multiple dimensions.
Function: ov.pl.single_group_boxplot
:
- adata (AnnData): Annotated data object containing single-cell RNA-seq data.
- groupby (str, optional): Grouping variable for the heatmap. Default is ''.
- figsize (tuple, optional): Figure size. Default is (6, 10).
- layer (str, optional): Data layer to use. Default is None.
- use_raw (bool, optional): Whether to use the raw data. Default is False.
- var_names (list or None, optional): List of genes to include in the heatmap. Default is None.
- gene_symbols (None, optional): Not used in the function.
- standard_scale (str, optional): Method for standardizing values. Options: 'obs', 'var', None. Default is None.
- col_color_bars (dict, optional): Dictionary mapping columns types to colors.
- col_color_labels (dict, optional): Dictionary mapping column labels to colors.
- left_color_bars (dict, optional): Dictionary mapping left types to colors.
- left_color_labels (dict, optional): Dictionary mapping left labels to colors.
- right_color_bars (dict, optional): Dictionary mapping right types to colors.
- right_color_labels (dict, optional): Dictionary mapping right labels to colors.
- marker_genes_dict (dict, optional): Dictionary mapping cell types to marker genes.
- index_name (str, optional): Name for the index column in the melted DataFrame. Default is ''.
- value_name (str, optional): Name for the value column in the melted DataFrame. Default is ''.
- cmap (str, optional): Colormap for the heatmap. Default is 'parula'.
- xlabel (str, optional): X-axis label. Default is ''.
- ylabel (str, optional): Y-axis label. Default is ''.
- label (str, optional): Label for the plot. Default is ''.
- save (bool, optional): Whether to save the plot. Default is False.
- save_pathway (str, optional): File path for saving the plot. Default is ''.
- legend_gap (int, optional): Gap between legend items. Default is 7.
- legend_hpad (int, optional): Horizontal space between the heatmap and legend, default is 2 [mm].
- show (bool, optional): Whether to display the plot. Default is False.
import pandas as pd
marker_genes_dict = {
'Sox':['Sox4', 'Sox7', 'Sox18', 'Sox5'],
}
color_dict = {'Sox':'#EFF3D8',}
gene_color_dict = {}
gene_color_dict_black = {}
for cell_type, genes in marker_genes_dict.items():
cell_type_color = color_dict.get(cell_type)
for gene in genes:
gene_color_dict[gene] = cell_type_color
gene_color_dict_black[gene] = '#000000'
cm = ov.pl.complexheatmap(adata,
groupby ='clusters',
figsize =(5,2),
layer = None,
use_raw = False,
standard_scale = 'var',
col_color_bars = dict(zip(pd.Categorical(adata.obs['clusters']).categories, adata.uns['clusters_colors'])),
col_color_labels = dict(zip(pd.Categorical(adata.obs['clusters']).categories, adata.uns['clusters_colors'])),
left_color_bars = color_dict,
left_color_labels = None,
right_color_bars = color_dict,
right_color_labels = gene_color_dict_black,
marker_genes_dict = marker_genes_dict,
cmap = 'coolwarm', #parula,jet
legend_gap = 15,
legend_hpad = 0,
left_add_text = True,
col_split_gap = 2,
row_split_gap = 1,
col_height = 6,
left_height = 4,
right_height = 6,
col_split = None,
row_cluster = False,
col_cluster = False,
value_name='Gene',
xlabel = "Expression of selected genes",
label = 'Gene Expression',
save = True,
show = False,
legend = False,
plot_legend = False,
#save_pathway = "complexheatmap.png",
)
PyComplexHeatmap have been install version: 1.7.3 Starting plotting.. Starting calculating row orders.. Reordering rows.. Starting calculating col orders.. Reordering cols.. Plotting matrix.. Starting plotting HeatmapAnnotations Starting plotting HeatmapAnnotations Starting plotting HeatmapAnnotations Collecting legends.. Collecting annotation legends.. Collecting annotation legends.. Collecting annotation legends..