RNA-seq was performed on one fresh-frozen sample from P1 and one FFPE sample from P2. RNA was extracted using Promega Maxwell 16 MDX and RNA integrity was verified using the Agilent Bioanalyzer 2100 (Agilent Technologies). cDNA was synthesized from total RNA using Superscript III (Invitrogen). Sequencing was performed using Illumina HiSeq2500 in 2 × 75 bp paired-end mode, and generated an average of 45 million reads per sample. All reads were independently aligned with STAR_2.4.0f for sequence alignment against the human genome build hg19, downloaded via the UCSC genome browser (http://hgdownload.soe.ucsc.edu/goldenPath/hg19/bigZips/), and SAMTOOLS v0.1.19 for sorting and indexing reads (Li et al. 2009 (link)). For fusion analysis, we used STAR-fusion (STAR-Fusion_v0.5.1) and FusionCatcher (v0.99.3e). Fusions with significant support of junction reads and spanning pairs were selected and manually reviewed.