Visualization of single cell RNA-seq¶

In this part, we will introduce the tutorial of special plot of omicverse.

We utilized single-cell RNA-seq data (GEO accession: GSE95753) obtained from the dentate gyrus of the hippocampus in mouse.

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 to X_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 prepared ov.pl.red_color,ov.pl.blue_color,ov.pl.green_color,ov.pl.orange_color, by default.

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

If you are interested in the changes in cell types in different groups, we recommend using a stacked area graph.

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 to X_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.

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 to X_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.

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 to X_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

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 to X_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

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 to X_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

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).

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 instead p<0.001
• fontsize: the fontsize of text
• fontcolor: the color of text
• horizontalalignment: the location of text

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).

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.

Marker gene plot¶

In single-cell analysis, a marker gene heatmap is a powerful visualization tool that helps researchers to understand the expression patterns of specific marker genes across different cell populations. Here we provide ov.pl.marker_heatmap for visualizing the patterns of marker genes.

We first preprocess the data and define the dictionary of cell type and marker gene. Please ensure that each gene in the dictionary appears only once (i.e. different cells cannot have the same marker gene, otherwise an error will be reported).

Function: ov.pl.marker_heatmap:

• adata: AnnData object Annotated data matrix.
• marker_genes_dict: dict A dictionary containing the marker genes for each cell type.
• groupby: str The key in adata.obs that will be used for grouping the cells.
• color_map: str The color map to use for the value of heatmap.
• use_raw: bool Whether to use the raw data of AnnDta object for plotting.
• standard_scale: str The standard scale for the heatmap.
• expression_cutoff: float The cutoff value for the expression of genes.
• bbox_to_anchor: tuple The position of the legend bbox (x, y) in axes coordinates.
• figsize: tuple The size of the plot figure in inches (width, height).
• spines: bool Whether to show the spines of the plot.
• fontsize: int The font size of the text in the plot.
• show_rownames: bool Whether to show the row names in the heatmap.
• show_colnames: bool Whether to show the column names in the heatmap.
• save_pathway: str The file path for saving the plot.
• ax: matplotlib.axes.Axes A pre-existing axes object for plotting (optional).