Your submission: WTC11_PacBioOnly on PacBio data
Background
Challenge 1 is evaluated according to four criteria:
- Broad GENCODE Annotation
- Subset of manually curated loci selected by GENCODE
- SIRV Lexogen Set 4
- Simulated data.
The LRGASP uses SQANTI categories to define evaluating features and metrics for Challenge 1.
LRGASP Challenge 1 Definitions:
- Full Splice Match (FSM): Transcripts matching a reference transcript at all splice junctions.
- Incomplete Splice Match (ISM): Transcripts matching consecutive, but not all, splice junctions of the reference transcripts.
- Novel in Catalog (NIC): Transcripts containing new combinations of already annotated splice junctions or novel splice junctions formed from already annotated donors and acceptors.
- Novel Not in Catalog (NNC): Transcripts using novel donors and/or acceptors.
- Reference Match (RM): FSM transcript with 5´ and 3´ends within 50 nts of the TSS/TTS annotation.
- 3´ polyA supported: Transcript with polyA motif support at the 3´end.
- 5´ CAGE supported: Transcript with CAGE support at the 5´end.
- 3´ reference supported: Transcript with 3´end within 50 nts from reference transcript or gene TTS.
- 5´ reference supported: Transcript with 5´end within 50 nts from reference transcript or gene TSS.
- Supported Reference Transcript Model (SRTM): FSM/ISM transcript with 5´ end within 50nts of the TSS or has CAGE support AND 3´ end within 50nts of the TTS or has polyA motif support
- Supported Novel Transcript Model (SNTM): NIC/NNC transcript with 5´ end within 50nts of the TSS or CAGE support AND 3´ end within 50nts of the TTS or polyA motif support AND Illumina read support at novel junctions
- Redundancy: # LR transcript models / reference model
- Intron retention (IR): level Number of IR within the NIC category
- Illumina Splice Junction (SJ): Support % SJ in transcript model with Illumina support
- Full Illumina Splice Junction Support: % transcripts in category with all SJ supported
- % Novel Junctions: # of new junctions / total # junctions
- % Non-canonical junctions: # of non-canonical junctions / total # junctions
- % Non-canonical transcripts % transcripts with at least one nc junction
- Intra-priming: Evidence of intra-priming (see SQANTI ref)
- RT-switching: Evidence of RT-switching (see SQANTI ref)
This document shows the performance of your pipeline for criteria 1. and 3. Critical data for evaluation according to 2. and 4. will be made available after the closure of the challenge, and therefore pre-evaluation reports cannot be provided. Note you’re your criterion 1 metrics reported here have been calculated using GENCODE human v38 and mouse M27 releases while the final evaluation will use human v39 and mouse M28 to be released after completing of the challenge.
Evaluation of detected transcripts for Challenge 1
Here, SIRVs are not taken into account for none of the metrics.
Global overview
| Number of genes detected |
18427 |
| Number of known genes detected |
13862 |
| Number of transcripts detected |
55724 |
| Number of transcripts associated to a known gene |
50411 |
| Number of unique SJ detected |
140927 |
|
| Novel SJ |
5434 |
0.04 |
| Non-canonical SJ |
418 |
0.00 |
|
Evaluation of FSM
| Number of isoforms |
31713 |
- |
| Reference Match |
12896 |
40.66 |
| 5’ reference supported (transcript) |
22998 |
72.52 |
| 3’ reference supported (transcript) |
17729 |
55.9 |
| 5’ reference supported (gene) |
27898 |
87.97 |
| 3’ reference supported (gene) |
24962 |
78.71 |
| 5’ CAGE supported |
27858 |
87.84 |
| 3’ polyA supported |
31060 |
97.94 |
| Supported Reference Transcript Model (SRTM) |
29665 |
93.54 |
| Intra-priming |
845 |
2.66 |
| Reference redundancy Level |
1.51 |
- |
|
Evaluation of ISM
| Number of isoforms |
3036 |
- |
| 5’ reference supported (transcript) |
1216 |
40.05 |
| 3’ reference supported (transcript) |
1372 |
45.19 |
| 5’ and 3’ reference supported (gene) |
1023 |
33.7 |
| 5’ reference supported (gene) |
1790 |
1790 |
| 3’ reference supported (gene) |
2045 |
67.36 |
| 5’ CAGE supported |
1985 |
65.38 |
| 3’ polyA supported |
2763 |
91.01 |
| Supported Reference Transcript Model (SRTM) |
1959 |
64.53 |
| Intra-priming |
231 |
7.61 |
| Reference redundancy Level |
1.3 |
- |
|
Evaluation NIC
| Number of isoforms |
12635 |
- |
| 5’ and 3’ reference supported (gene) |
9935 |
78.63 |
| 5’ reference supported (gene) |
11566 |
91.54 |
| 3’ reference supported (gene) |
10764 |
85.19 |
| 5’ CAGE supported |
11347 |
89.81 |
| 3’ polyA supported |
12318 |
97.49 |
| Supported Novel Transcript Model (SNTM) |
10853 |
85.9 |
| Intra-priming |
419 |
3.32 |
| Intron retention incidence |
3447 |
27.28 |
|
Evaluation NNC
| Number of isoforms |
3027 |
- |
| 5’ and 3’ reference supported (gene) |
1587 |
52.43 |
| 5’ reference supported (gene) |
2189 |
72.32 |
| 3’ reference supported (gene) |
2167 |
71.59 |
| 5’ CAGE supported |
2550 |
84.24 |
| 3’ polyA supported |
2962 |
97.85 |
| Supported Novel Transcript Model (SNTM) |
2318 |
76.58 |
| Non-canonical SJ incidence |
334 |
11.03 |
| Full Illumina SJ support |
2587 |
85.46 |
| Intra-priming |
109 |
3.6 |
| RT-switching incidence |
228 |
7.53 |
|
Evaluation of Spike-Ins (SIRVs)
The following metrics and definitions apply to SIRV transcripts:
- SIRV transcript: Transcript mapping to a SIRV sequence
- Reference SIRV (rSIRV): Lexogen SIRV model
- True Positive detections (TP): rSIRVs identified as RM
- Partial True Positive detections (PTP): rSIRVs identified as ISM or FSM_non_RM
- False Negative (FN): rSIRVs without FSM or ISM
- False Positive (FP): NIC + NNC + antisense + fusion SIRV_transcripts
- Sensitivity: TP/rSIRVs
- Precision: RM/ SIRV_transcripts
- Non_redundant Precision: TP/ SIRV_transcripts
- Positive Detection Rate: unique(TP+PTP)/rSIRVs
- False Discovery Rate: (FP + PTP)/SIRV_transcripts
- False Detection Rate: (FP)/SIRV_transcripts
- Redundancy: (FSM + ISM)/unique(TP+PTP)
ATENTION If in this chunk of the evaluation all the results are 0, please, check if the reference genome and/or transcriptome used for building your transcript-model contain information about spike-ins.
| SIRV transcripts |
85 |
| True Positive detections (TP) |
81 |
| SIRV transcripts associated to TP (Reference Match) |
81 |
| Partial True Positive detections (PTP) |
1 |
| SIRV transcripts associated to PTP |
2 |
| False Negative (FN) |
2 |
| False Positive (FP) |
2 |
| Sensitivity |
0.96 |
| Precision |
0.95 |
| Non Redundant Precision |
0.95 |
| Positive Detection Rate |
0.98 |
| False Discovery Rate |
0.04 |
| False Detection Rate |
0.02 |
| Redundancy |
1.01 |
|