FindPTC.A3SS.PosStrand {MARVEL}

Find premature terminal codon (PTC) for alternative 3' splice site (A3SS) located on the positive strand of the transcript

Description

Finds PTC(s) introduced by alternative exons into protein-coding transcripts.

Usage

FindPTC.A3SS.PosStrand(MarvelObject, tran_id, gene_id)

Arguments

MarvelObject

Marvel object. S3 object generated from CompareValues.PSI and ParseGTF function.

tran_id

Character string. Vector of tran_id to look for PTCs.

gene_id

Character string. Vector of gene_id corresponding to the tran_id argument.

Value

A data frame of transcripts containing splicing events meeting the psi.de.sig and psi.de.diff criteria are categorised based on the presence or absence of PTCs.

Examples

marvel.demo <- readRDS(system.file("extdata/data", "marvel.demo.rds", package="MARVEL"))

# Define relevant event type
results <- marvel.demo$DE$PSI$Table[["ad"]]
index.1 <- which(results$event_type=="A3SS")
index.2 <- grep(":+@", results$tran_id, fixed=TRUE)
index <- intersect(index.1, index.2)
results <- results[index, ]
tran_id <- results$tran_id[1]
gene_id <- results$gene_id[1]

# Find PTC
results <- FindPTC.A3SS.PosStrand(MarvelObject=marvel.demo,
                                tran_id=tran_id,
                                gene_id=gene_id
                                )

# Check output
head(results)