Expression tracks
Both tracks can be shown in a Table (
) and a Graphical (
) view.
The expression track table view has the following options (figure 31.11).
- The "Filter to selection" only displays pre-selected rows in the table.
- The "Create track from Selection" will create a Track using selected rows.
- The "Select Genes in Other Views" button finds and selects the currently selected genes and transcripts in all other open expression track table views.
- The "Copy Gene Names to Clipboard" button copies the currently selected gene names to the clipboard.
Figure 31.11: RNA-Seq results shown in a table view.
By creating a Track list, the graphical view can be shown together with the read mapping track and tracks from other samples:
File | New | Track List (
)
Select the mapping and expression tracks of the samples you wish to visualize together and select the annotation tracks used as reference for the RNA-Seq and click Finish.
Once the track list is shown, double-click the label of the expression track to show it in a table view.
Clicking a row in the table makes the track list view jump to that location, allowing for quick inspection of interesting parts of the RNA-Seq read mapping (see an example in figure 31.12).
Figure 31.12: RNA-Seq results shown in a split view with an expression track at the bottom and a track list with read mappings of two samples at the top.
Reads spanning two exons are shown with a dashed line between each end as shown in figure 31.12, and the thin solid line represents the connection between two reads in a pair.
When doing comparative analysis, double click on one of the Expression or Statistical Comparison tracks in a track list to get its table view. Then click on the "Select genes in other views" button in any other table or expression browser will cause the track list to zoom to the selected gene.
Expression tracks can also be used to annotate variants using the Annotate with Overlap Information tool. Select the variant track as input and annotate with the expression track. For variants inside genes or transcripts, information will be added about expression (counts, expression value etc) from the gene or transcript in the expression track. Read more about the annotation tool in Annotate with Overlap Information.
Gene-level expression
The gene-level expression track (GE) holds information about counts and expression values for each gene. It can be opened in a Table view () allowing sorting and filtering on all the information in the track (see figure 31.13 for an example subset of an expression track).
Figure 31.13: A subset of a result of an RNA-Seq analysis on the gene level. Not all columns are shown in this figure
Each row in the table corresponds to a gene (or reference sequence, if the One reference sequence per transcript option was used). The corresponding counts and other information is shown for each gene:
- Name. The name of the gene, or the name of the reference sequence if "one reference sequence per transcript" is used.
- Chromosome and region. The position of the gene on the genome.
- Expression value. This is based on the expression measure chosen as described in Calculating expression values from RNA-Seq.
- Gene length The length of the gene as annotated.
- TPM (Transcripts per million). This is computed as
, where the sum is over the RPKM values of all genes/transcripts (see http://bioinformatics.oxfordjournals.org/content/26/4/493.long).
- RPKM. This is the expression value measured in RPKM [Mortazavi et al., 2008]:
. See RPKM for a detailed definition.
- Unique gene reads. This is the number of reads that match uniquely to the gene or its transcripts.
- Total gene reads. This is all the reads that are mapped to this gene - both reads that map uniquely to the gene or its transcripts and reads that matched to more positions in the reference (but fewer than the 'Maximum number of hits for a read' parameter) which were assigned to this gene.
- Transcripts annotated. The number of transcripts annotated for the gene. Note that this is not based on the sequencing data - only on the annotations already on the reference sequence(s).
- Uniquely identified transcripts. The number of transcripts with at least one mapped read that matches only that transcript and no others. Note that if a gene has 4 detected transcripts, and 8 undetected transcripts, all 4+8=12 transcripts will have the value "Uniquely identified transcripts = 4".
- Exon length. The total length of all exons (not all transcripts).
- Exons. The total number of exons across all transcripts.
- Unique exon reads. The number of reads that match uniquely to the exons (including across exon-exon junctions).
- Total exon reads. The total number of reads assigned to an exon or an exon-exon junction of this gene. As for the 'Total gene reads' this includes both uniquely mapped reads and reads with multiple matches that were assigned to an exon of this gene.
- Ratio of unique to total (exon reads). The ratio of the unique reads to the total number of reads in the exons. This can be convenient for filtering the results to exclude the ones where you have low confidence because of a relatively high number of non-unique exon reads.
- Unique exon-exon reads. Reads that uniquely match across an exon-exon junction of the gene (as specified in figure 31.12). The read is only counted once even though it covers several exons.
- Total exon-exon reads. Reads that match across an exon-exon junction of the gene (as specified in figure 31.12). As for the 'Total gene reads' this includes both uniquely mapped reads and reads with multiple matches that were assigned to an exon-exon junction of this gene.
- Total intron reads. The total number of reads that map to an intron of the gene.
- Ratio of intron to total gene reads. This can be convenient to identify genes with poor or lacking transcript annotations. If one or more exons are missing from the annotations, there will be a relatively high number of reads mapping in the intron.
Transcript-level expression
If the "Genome annotated with genes and transcripts" option is selected in figure 31.4, a transcript-level expression track (TE) is also generated.
The track can be opened in a Table view () allowing sorting and filtering on all the information in the track. Each row in the table corresponds to an mRNA annotation in the mRNA track used as reference.
- Name. The name of the transcript, or the name of the reference sequence if "one reference sequence per transcript" is used.
- Chromosome and region. The position of the gene on the genome.
- Expression value. This is based on the expression measure chosen as described in Calculating expression values from RNA-Seq.
- TPM (Transcripts per million). This is computed as
, where the sum is over the RPKM values of all genes/transcripts (see http://bioinformatics.oxfordjournals.org/content/26/4/493.long).
- RPKM. This is the expression value measured in RPKM [Mortazavi et al., 2008]:
. See RPKM for a detailed definition.
- Relative RPKM. The RPKM for the transcript divided by the maximum of the RPKM values among all transcripts of the same gene. This value describes the relative expression of alternative transcripts for the gene.
- Gene name. The name of the corresponding gene.
- Transcript length. This is the length of the transcript.
- Exons. The total number of exons in the transcript.
- Transcript ID. The transcript ID is taken from the transcript_id note in the mRNA track annotations and can be used to differentiate between different transcripts of the same gene.
- Transcripts annotated. The number of transcripts based on the mRNA annotations on the reference. Note that this is not based on the sequencing data - only on the annotations already on the reference sequence(s).
- Uniquely identified transcripts. The number of transcripts with at least one mapped read that matches only that transcript and no others. Note that if a gene has 4 detected transcripts, and 8 undetected transcripts, all 4+8=12 transcripts will have the value "Uniquely identified transcripts = 4".
- Unique transcript reads. This is the number of reads in the mapping for the gene that are uniquely assignable to the transcript.
- Total transcript reads. Once the 'Unique transcript read's have been identified and their counts calculated for each transcript, the remaining (non-unique) transcript reads are assigned to one of the transcripts to which they match. The 'Total transcript reads' counts are the total number of reads that are assigned to the transcript once this assignment has been done. As for the assignment of reads among genes, the assignment of reads within a gene but among transcripts, is done by the EM estimation algorithm (EM estimation algorithm).
- Ratio of unique to total (transcript reads). The ratio of the unique reads to the total number of reads in the transcripts. This can be convenient for filtering the results to exclude the ones where you have low confidence because of a relatively high number of non-unique transcript reads.
Additional information to GE or TE tracks
Both GE and TE tables can offer additional information such as hyperlinks to various databases (e.g., ENSEMBL, HGNC (HUGO Gene Nomenclature Committee), RefSeq, GeneID, etc.) In cases where the mRNA track or the gene track provided have biotype information, a biotype column will be added to the table.