This function create a set of synthetic transcripts from a provided annotation file in "gff3" or "gtf" format. As detailed in http://www.epigenesys.eu/en/protocols/bio-informatics/1283-guidelines-for-rna-seq-data-analysis, one major caveat of estimating gene expression using aligned RNA-Seq reads is that a single read, which originated from a single mRNA molecule, might sometimes align to several features (e.g. transcripts or genes) with alignments of equivalent quality. This, for example, might happen as a result of gene duplication and the presence of repetitive or common domains. To avoid counting unique mRNA fragments multiple times, the stringent approach is to keep only uniquely mapping reads - being aware of potential consequences. Not only can "multiple counting" arise from a biological reason, but also from technical artifacts, introduced mostly by poorly formatted gff3/gtf annotation files. To avoid this, it is best practice to adopt a conservative approach by collapsing all existing transcripts of a single gene locus into a "synthetic" transcript containing every exon of that gene. In the case of overlapping exons, the longest genomic interval is kept, i.e. an artificial exon is created. This process results in a flattened transcript - a gene structure with a one (gene) to one (transcript) relationship.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
## S4 method for signature 'AnnotParamCharacter' createSyntheticTranscripts( obj, features = c("mRNA", "miRNA", "tRNA", "transcript"), verbose = TRUE ) ## S4 method for signature 'character' createSyntheticTranscripts( obj, features = c("mRNA", "miRNA", "tRNA", "transcript"), verbose = TRUE, output = c("Genome_intervals", "GRanges"), input = c("gff3", "gtf") )
one or more of 'mRNA', 'miRNA', 'tRNA', 'transcript'
increase the verbosity (default TRUE)
the output type, one of 'Genome_intervals' or 'GRanges'
the type of input, one of 'gff3' or 'gtf'
createSyntheticTranscripts function implements this, taking
advantage of the hierarchical structure of the gff3/gtf file. Exon
features are related to their transcript (parent), which themselves derives
from their gene parents. Using this relationship, exons are combined per gene
into a flattened transcript structure. Note that this might not avoid multiple
counting if genes overlap on opposing strands. There, only strand specific
sequencing data has the power to disentangle these situations.
As gff3/gtf file can contain a large number of feature types, the
createSyntheticTranscripts currently only supports: mRNA,
miRNA, tRNA and transcript. Please contact me if you
need additional features to be considered. Note however, that I will only
add features that are part of the sequenceontology.org SOFA
Depending on the
character filename: depending on the selected
For the input:
For the output:
1 2 3 4 5 6 7 8 9 10
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.