MARVEL (version 1.4.0)

FindPTC.RI.PosStrand: Find premature terminal codon (PTC) for retained-Intron (RI) located on the positive strand of the transcript

Description

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

Usage

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

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.

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.

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=="RI")
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.RI.PosStrand(MarvelObject=marvel.demo,
                                tran_id=tran_id,
                                gene_id=gene_id
                                )

# Check output
head(results)