Fusion transcriptome profiling provides insights into alveolar rhabdomyosarcoma

Gene fusions and fusion products were thought to be unique features of neoplasia. However, more and more studies have identified fusion RNAs in normal physiology. Through RNA sequencing of 27 human noncancer tissues, a large number of fusion RNAs were found. By analyzing fusion transcriptome, we observed close clusterings between samples of same or similar tissues, supporting the feasibility of using fusion RNA profiling to reveal connections between biological samples. To put the concept into use, we selected alveolar rhabdomyosarcoma (ARMS), a myogenic pediatric cancer whose exact cell of origin is not clear. PAX3-FOXO1 (paired box gene 3 fused with forkhead box O1) fusion RNA, which is considered a hallmark of ARMS, was recently found during normal muscle cell differentiation. We performed and analyzed RNA sequencing from various time points during myogenesis and uncovered many chimeric fusion RNAs. Interestingly, we found that the fusion RNA profile of RH30, an ARMS cell line, is most similar to the myogenesis time point when PAX3-FOXO1 is expressed. In contrast, full transcriptome clustering analysis failed to uncover this connection. Strikingly, all of the 18 chimeric RNAs in RH30 cells could be detected at the same myogenic time point(s). In addition, the seven chimeric RNAs that follow the exact transient expression pattern as PAX3-FOXO1 are specific to rhabdomyosarcoma cells. Further testing with clinical samples also confirmed their specificity to rhabdomyosarcoma. These results provide further support for the link between at least some ARMSs and the PAX3-FOXO1-expressing myogenic cells and demonstrate that fusion RNA profiling can be used to investigate the etiology of fusion-gene-associated cancers.