Fujun Qin

Recurrent chimeric fusion RNAs in non-cancer tissues and cells

Gene fusions and their products (RNA and protein) were once thought to be unique features to cancer. However, chimeric RNAs can also be found in normal cells. Here, we performed, curated and analyzed nearly 300 RNA-Seq libraries covering 30 different non-neoplastic human tissues and cells as well as 15 mouse tissues. A large number of fusion transcripts were found. Most fusions were detected only once, while 291 were seen in more than one sample. We focused on the recurrent fusions and performed RNA and protein level validations on a subset. We characterized these fusions based on various features of the fusions, and their parental genes. They tend to be expressed at higher levels relative to their parental genes than the non-recurrent ones. Over half of the recurrent fusions involve neighboring genes transcribing in the same direction. A few sequence motifs were found enriched close to the fusion junction sites. We performed functional analyses on a few widely expressed fusions, and found that silencing them resulted in dramatic reduction in normal cell growth and/or motility. Most chimeras use canonical splicing sites, thus are likely products of ‘intergenic splicing’. We also explored the implications of these non-pathological fusions in cancer and in evolution.

Comparative assessment of methods for the fusion transcripts detection from RNA-Seq data.

RNA-Seq made possible the global identification of fusion transcripts, i.e. “chimeric RNAs”. Even though various software packages have been developed to serve this purpose, they behave differently in different datasets provided by different developers. It is important for both users, and developers to have an unbiased assessment of the performance of existing fusion detection tools. Toward this goal, we compared the performance of 12 well-known fusion detection software packages. We evaluated the sensitivity, false discovery rate, computing time, and memory usage of these tools in four different datasets (positive, negative, mixed, and test). We conclude that some tools are better than others in terms of sensitivity, positive prediction value, time consumption and memory usage. We also observed small overlaps of the fusions detected by different tools in the real dataset (test dataset). This could be due to false discoveries by various tools, but could also be due to the reason that none of the tools are inclusive. We have found that the performance of the tools depends on the quality, read length, and number of reads of the RNA-Seq data. We recommend that users choose the proper tools for their purpose based on the properties of their RNA-Seq data.

Discovery of CTCF-sensitive Cis-spliced fusion RNAs between adjacent genes in human prostate cells.

Genes or their encoded products are not expected to mingle with each other unless in some disease situations. In cancer, a frequent mechanism that can produce gene fusions is chromosomal rearrangement. However, recent discoveries of RNA trans-splicing and cis-splicing between adjacent genes (cis-SAGe) support for other mechanisms in generating fusion RNAs. In our transcriptome analyses of 28 prostate normal and cancer samples, 30% fusion RNAs on average are the transcripts that contain exons belonging to same-strand neighboring genes. These fusion RNAs may be the products of cis-SAGe, which was previously thought to be rare. To validate this finding and to better understand the phenomenon, we used LNCaP, a prostate cell line as a model, and identified 16 additional cis-SAGe events by silencing transcription factor CTCF and paired-end RNA sequencing. About half of the fusions are expressed at a significant level compared to their parental genes. Silencing one of the in-frame fusions resulted in reduced cell motility. Most out-of-frame fusions are likely to function as non-coding RNAs. The majority of the 16 fusions are also detected in other prostate cell lines, as well as in the 14 clinical prostate normal and cancer pairs. By studying the features associated with these fusions, we developed a set of rules: 1) the parental genes are same-strand-neighboring genes; 2) the distance between the genes is within 30kb; 3) the 5′ genes are actively transcribing; and 4) the chimeras tend to have the second-to-last exon in the 5′ genes joined to the second exon in the 3′ genes. We then randomly selected 20 neighboring genes in the genome, and detected four fusion events using these rules in prostate cancer and non-cancerous cells. These results suggest that splicing between neighboring gene transcripts is a rather frequent phenomenon, and it is not a feature unique to cancer cells.

A Chimeric RNA Characteristic of Rhabdomyosarcoma in Normal Myogenesis Process

UNLABELLED Gene fusions and their chimeric products are common features of neoplasia. Given that many cancers arise by the dysregulated recapitulation of processes in normal development, we hypothesized that comparable chimeric gene products may exist in normal cells. Here, we show that a chimeric RNA, PAX3-FOXO1, identical to that found in alveolar rhabdomyosarcoma, is transiently present in cells undergoing differentiation from pluripotent cells into skeletal muscle. Unlike cells of rhabdomyosarcoma, these cells do not seem to harbor the t(2;13) chromosomal translocation. Importantly, both PAX3-FOXO1 RNA and protein could be detected in the samples of normal fetal muscle. Overexpression of the chimera led to continuous expression of MYOD and MYOG-two myogenic markers that are overexpressed in rhabdomyosarcoma cells. Our results are consistent with a developmental role of a specific chimeric RNA generated in normal cells without the corresponding chromosomal rearrangement at the DNA level seen in neoplastic cells presumably of the same lineage. SIGNIFICANCE A chimeric fusion RNA, PAX3-FOXO1, associated with alveolar rhabdomyosarcoma, is also present in normal non-cancer cells and tissues. Its transient expression nature and the absence of t(2;13) chromosomal translocation are consistent with a posttranscriptional mechanism. When constantly expressed, PAX3-FOXO1 interfered with the muscle differentiation process, which presumably contributes to tumorigenesis.

Recurrent fusion RNA DUS4L-BCAP29 in non-cancer human tissues and cells

Traditional gene fusions are involved in the development of various neoplasia. DUS4L-BCAP29, a chimeric fusion RNA, has been reported to be a cancer-fusion in prostate and gastric cancer, in addition to playing a tumorigenic role. Here, we showed that the DUS4L-BCAP29 fusion transcript exists in a variety of normal tissues. It is also present in non-cancer epithelial, as well as in fibroblast cell lines. Quantitatively, the fusion transcript has a comparable expression in non-cancerous, gastric and prostate cell lines and tissues as in the cancer cell lines and tissues. The loss-of-function approach as previously reported is not sufficient to prove the functionality of the fusion. On the other hand, the gain-of-function approach showed that overexpression of DUS4L-BCAP29 promotes cell growth and motility, even in non-cancer cells. Finally, we provide further evidence that the fusion transcript is a product of cis-splicing between adjacent genes. In summary, we believe that in contrast to traditional gene fusions, DUS4L-BCAP29 cannot be used as a cancer biomarker. Instead, it is a fusion transcript that exists in normal physiology and that its pro-growth effect is not unique to cancer cells.

Role of CTCF in Regulating SLC45A3-ELK4 Chimeric RNA.

The chimeric RNA, SLC45A3-ELK4, was found to be a product of cis-splicing between the two adjacent genes (cis-SAGe). Despite the biological and clinical significance of SLC45A3-ELK4, its generating mechanism has not been elucidated. It was shown in one cell line that the binding of transcription factor CTCF to the insulators located at or near the gene boundaries, inversely correlates with the level of the chimera. To investigate the mechanism of such cis-SAGe events, we sequenced potential regions that may play a role in such transcriptional read-through. We could not detect mutations at the transcription termination site, insulator sites, splicing sites, or within CTCF itself in LNCaP cells, thus suggesting a “soft-wired” mechanism in regulating the cis-SAGe event. To investigate the role CTCF plays in regulating the chimeric RNA expression, we compared the levels of CTCF binding to the insulators in different cell lines, as well as clinical samples. Surprisingly, we did not find an inverse correlation between CTCF level, or its bindings to the insulators and SLC45A3-ELK4 expression among different samples. However, in three prostate cancer cell lines, different environmental factors can cause the expression levels of the chimeric RNA to change, and these changes do inversely correlate with CTCF level, and/or its bindings to the insulators. We thus conclude that CTCF and its bindings to the insulators are not the primary reasons for differential SLC45A3-ELK4 expression in different cell lines, or clinical cases. However, they are the likely mechanism for the same cells to respond to different environmental cues, in order to regulate the expression of SLC45A3-ELK4 chimeric RNA. This response to different environmental cues is not general to other cis-SAGe events, as we only found one out of 16 newly identified chimeric RNAs showing a pattern similar to SLC45A3-ELK4.

Recurrent cis-SAGe chimeric RNA, D2HGDH-GAL3ST2, in prostate cancer

Neighboring genes transcribing in the same direction can form chimeric RNAs via cis-splicing (cis-SAGe). Previously, we reported 16 novel cis-SAGe chimeras in prostate cancer cell lines, and performed in silico validation on 14 pairs of normal and tumor samples from Chinese patients. However, whether these fusions exist in different populations, as well as their clinical implications, remains unclear. To investigate, we developed a bioinformatics pipeline using modified Spliced Transcripts Alignment to a Reference (STAR) to quantify these fusion RNAs simultaneously in silico. From RNA-Seq data of 100 paired normal and prostate cancer samples from TCGA, we find that most fusions are not specific to cancer. However, D2HGDH-GAL3ST2 is more frequently seen in cancer samples, and seems to be enriched in the African American group. Further validation with our own collection as well as from commercial sources did not detect this fusion RNA in 29 normal prostate samples, but in 19 of 93 prostate cancer samples. It is more frequently detected in late stage cancer, suggesting a role in cancer progression. Consistently, silencing this fusion resulted in dramatic reduction of cell proliferation rate and cell motility.

Fusion transcriptome profiling provides insights into alveolar rhabdomyosarcoma

Significance Here, we propose an approach to study connections between biological samples. By using binary input of fusion RNA expression, samples of same or similar tissue origin were clustered together. The concept was then put into use to gain insights for the pediatric alveolar rhabdomyosarcoma (ARMS). We found that the signature fusion RNA for ARMS, PAX3–FOXO1 (paired box gene 3 fused with forkhead box O1), and all the other chimeric RNAs expressed in ARMS cells are expressed at the same normal myogenic time point(s). Several chimeric RNAs were further confirmed to be specifically expressed in clinical rhabdomyosarcoma tumor cases. These results support the link between at least some ARMS tumors and the PAX3–FOXO1-expressing myogenic cell. Fusion RNA profiling is a tool to investigate the etiology of fusion-gene-associated cancers. 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.

SLC45A3-ELK4 functions as a long non-coding chimeric RNA

Gene fusions in cancer typically lead to the expression of a fusion protein or disrupt the expression of one of the parental genes. Here we report a new phenomenon whereby a fusion transcript functions as a long non-coding chimeric RNA (lnccRNA). This fusion RNA, SLC45A3-ELK4, generated by cis-splicing between neighboring genes, was found in prostate cancer. The fusion RNA encodes the same protein as ELK4. Intriguingly, we found that the fusion RNA level is less than 1% of wild type ELK4, unlikely to perturb the general pool of ELK4 protein. Nonetheless, when the fusion RNA, but not ELK4 is silenced, cell proliferation is inhibited in both androgen-dependent and castration-resistant prostate cancer cells. This growth arrest can be rescued by exogenous expression of the fusion and a mutant designed to prevent translation of the ELK4 protein. In the same setting, the mutant could also suppress CDKN1A and several other targets of SLC45A3-ELK4. In addition, similar to many long non-coding RNAs, the fusion RNA is enriched in the nuclear fraction. Altogether, these results indicate that SLC45A3-ELK4 regulates cancer cell proliferation by its transcript, not translated protein.

Fas-associated factor 1 inhibits tumor growth by suppressing Helicobacter pylori -induced activation of NF-κB signaling in human gastric carcinoma

Loss of Fas-associated factor 1 (FAF1) may act as a pro-survival signal in diseased cells, but whether this is true in gastric carcinoma remains unclear. Here we report that FAF1 was expressed at low levels in gastric carcinoma tissues and cell lines, and its expression correlated with larger tumors, higher histology grade, higher TNM stage, tumor infiltration, and lymph node metastasis. Univariate analysis and survival curve analysis identified low FAF1 expression as a predictor of poor prognosis. FAF1 overexpression in HGC-27 gastric cancer cells induced cell apoptosis and inhibited cell proliferation and growth. It also reduced colony formation in vitro and tumor growth in mice. We found that Helicobacter pylori, a risk factor for gastric cancer, down-regulated FAF1 expression via NF-κB signaling. Knock-down of IKKβ or p65 expression in gastric cancer cells reversed H. pylori-induced down-regulation of FAF1 expression and partially blocked H. pylori-induced secretion of inflammatory cytokines TNF-α and IL-8. Our results suggest that loss of FAF1 contributes to human gastric carcinogenesis by allowing H. pylori to activate NF-κB signaling.