Find DE genes between all TSCAN paths rooted from given cluster

Source: R/runTSCAN.R

This function finds all paths that root from a given cluster useCluster, and performs tests to identify significant features for each path, and are not significant and/or changing in the opposite direction in the other paths. Using a branching cluster (i.e. a node with degree > 2) may highlight features which are responsible for the branching event. MST has to be pre-calculated with runTSCAN.

runTSCANClusterDEAnalysis(
  inSCE,
  useCluster,
  useAssay = "logcounts",
  fdrThreshold = 0.05
)

Arguments

inSCE

Input SingleCellExperiment object.

useCluster

The cluster to be regarded as the root, has to existing in colData(inSCE)$TSCAN_clusters.

useAssay

Character. The name of the assay to use. This assay should contain log normalized counts. Default "logcounts".

fdrThreshold

Only out put DEGs with FDR value smaller than this value. Default 0.05.

Value

The input inSCE with results updated in metadata.

Author

Nida Pervaiz

Examples

data("mouseBrainSubsetSCE", package = "singleCellTK")
mouseBrainSubsetSCE <- runTSCAN(inSCE = mouseBrainSubsetSCE,
                                useReducedDim = "PCA_logcounts")
#> Fri May  3 07:19:38 2024 ... Running 'scran SNN clustering' with 'louvain' algorithm
#> Fri May  3 07:19:38 2024 ...   Identified 2 clusters
#> Fri May  3 07:19:38 2024 ... Running TSCAN to estimate pseudotime
#> Fri May  3 07:19:39 2024 ...   Clusters involved in path index 2 are: 1, 2
#> Fri May  3 07:19:39 2024 ...   Number of estimated paths is 1
mouseBrainSubsetSCE <- runTSCANClusterDEAnalysis(inSCE = mouseBrainSubsetSCE,
                                         useCluster = 1)
#> Fri May  3 07:19:39 2024 ... Finding DEG between TSCAN branches
#> Fri May  3 07:19:39 2024 ...   Clusters involved in path index 2 are: 1, 2