SGSeq (version 1.0.6)

convertToSGFeatures: Convert transcript features to splice graph features

Description

Convert transcript features, predicted from RNA-seq data or extracted from transcript annotation, to splice graph features.

Usage

convertToSGFeatures(x, coerce = FALSE)

Arguments

x

TxFeatures object

coerce

Logical indicating whether transcript features should be coerced to splice graph features without disjoining exons and omitting splice donor and acceptor sites

Value

An SGFeatures object

Details

Splice junctions are unaltered. Exons are disjoined into non-overlapping exon bins. Adjacent exon bins without a splice site at the shared boundary are merged. All exon bins are assigned type “E”.

Entries for splice donor and acceptor sites (positions immediately upstream and downstream of introns, respectively) are added.

In the returned SGFeatures object, column slots splice5p and splice3p indicate whether compatibility with an exon bin requires a fragment to be spliced at the 5' or 3' boundary, respectively. splice5p (splice3p) is TRUE if the first (last) position of the exon coincides with a splice acceptor (donor), and it is not adjacent to a neighboring exon bin.

Each feature is assigned a unique feature and gene identifier, stored in column slots featureID and geneID, respectively. The latter indicates features that belong to the same gene, represented by a connected component in the splice graph.

Examples

sgf <- convertToSGFeatures(txf)