PlotPSI.SE.MultiGroups {VALERIE}

Percent spliced-in (PSI) visualization for skipped exons (SE)

Description

PlotPSI.SE visualizes percent spliced-in (PSI) for each genomic coordinate for skipped exon (SE) event across groups of single cells.

Usage

PlotPSI.SE.MultiGroups(object, SampleInfo, ExonInfo,
  statistical.test = c("KW", "ANOVA"), multiple.testing = c("holm",
  "hochberg", "hommel", "bonferroni", "BH", "BY", "fdr", "none"), Plots,
  plot.width, plot.height)

Arguments

object

Object of class rehab generated using ComputePSI.

SampleInfo

Tab-delimited file describing the naming and grouping of the single cells. First column should contain the names of the binary alignment map (BAM) files. Second column indicates the grouping for each single cell, namely Group1, Group2, etc. Third column indicates the group names. Example file provided in extdata directory of the package.

ExonInfo

Tab-delimited file describing the alternative splicing events. First column contains the alternative splicing nomenclature as per BRIE (Huang et al, Genome Biology, 2019) or MISO (Katz et al, Nature Methods, 2010). Second column indicates the type of alternative splicing event, namely SE, MXE, RI, A5SS, and A3SS. Third column contains the gene name or any personal notation. Example file provided in extdata directory of the package.

statistical.test

Method for comparising PSI values at each genomic coordinates between more groups of single cells.

multiple.testing

Method for adjusting p-values for multiple comparisons.

Plots

Folder to output PSI plots.

plot.width

Width of outplot plots.

plot.height

Height of outplot plots.

Details

This function visualizes the percent spliced-in (PSI) at each genomic coordinate encompassing the alternative exon and its flanking constitutive exons for each single cell in the form of a heatmap. The PSI mean for the respective groups are also display in the form of a line graph to summarize the PSI distributions of the respective groups. Pair-wise comparison of PSI at each genomic coordinate is performed using either the parametric ANOVA or non-parameteric Kruskal-Wallis test. The p-values can be adjusted for multiple testing using the p.adjust function.

Value

For each alternative splicing event, a single plot consisting of three subplots arranged from top to bottom is returned. Bottom subplot is a line graph of PSI means at each genomic coordinate for the groups of single cells. Middle subplot is a line graph of p-values corresponding to the comparison of PSI values at each genomic coordinate between groups of single cells. Top subplot is a heatmap of PSI values at each genomic coordinate across all single cells. Location of plots as per specified in the Plots argument.

Author(s)

Sean Wen <[e-mail]>

Examples

PSI <- readRDS(system.file("extdata/PSI", "PSI_RED_Multi_Groups_small.rds", package="VALERIE"))
PlotPSI.SE.MultiGroups(PSI, SampleInfo=system.file("extdata/Sample_Info",
  "Sample_Info_RED_Multi_Groups.txt", package="VALERIE"),
  ExonInfo=system.file("extdata/Exon_Info", "Exon_Info_RED_small.txt", package="VALERIE"),
  statistical.test="KW",  multiple.testing="bonferroni",
  Plots=tempdir(), plot.width=5, plot.height=8)