FusionNeoAntigen Logo

Home

Download

Statistics

Examples

Help

Contact

Terms of Use

Center for Computational Systems Medicine
leaf

Fusion Gene and Fusion Protein Summary

leaf

Fusion Amino Acid Sequences (multiple BPs and multiple gene isoforms)

leaf

Fusion Protein Breakpoint Sequences - (for the Screening of the FusionNeoAntigens)

leaf

Potential FusionNeoAntigens in HLA I - (netMHCpan v4.1 + deepHLApan v1.1)

leaf

Potential FusionNeoAntigens in HLA II - (netMHCIIpan v4.1)

leaf

Fusion Breakpoint 14 AA Peptide Structure - (RoseTTAFold)

leaf

Filtering FusionNeoAntigens Through Checking the Interaction with HLAs in 3D - (Glide)

leaf

Vaccine Design for the FusionNeoAntigens (RNA/protein sequences)

leaf

Potential target of CAR-T therapy development

leaf

Information on the samples that have these potential fusion neoantigens

leaf

Fusion Protein Targeting Drugs - (Manual Curation)

leaf

Fusion Protein Related diseases - (Manual Curation)

Fusion Protein:RARA-IGFBP4

Fusion Gene and Fusion Protein Summary

check button Fusion gene summary
Fusion partner gene informationFusion gene name: RARA-IGFBP4
FusionPDB ID: 72291
FusionGDB2.0 ID: 72291
HgeneTgene
Gene symbol

RARA

IGFBP4

Gene ID

5914

3487

Gene nameretinoic acid receptor alphainsulin like growth factor binding protein 4
SynonymsNR1B1|RARBP-4|HT29-IGFBP|IBP4|IGFBP-4
Cytomap

17q21.2

17q21.2

Type of geneprotein-codingprotein-coding
Descriptionretinoic acid receptor alphaRAR-alphanuclear receptor subfamily 1 group B member 1nucleophosmin-retinoic acid receptor alpha fusion protein NPM-RAR long formretinoic acid nuclear receptor alpha variant 1retinoic acid nuclear receptor alpha variant 2insulin-like growth factor-binding protein 4IBP-4IGF-binding protein 4
Modification date2020032720200315
UniProtAcc

P10276

Main function of 5'-partner protein: FUNCTION: Receptor for retinoic acid (PubMed:19850744, PubMed:16417524, PubMed:20215566). Retinoic acid receptors bind as heterodimers to their target response elements in response to their ligands, all-trans or 9-cis retinoic acid, and regulate gene expression in various biological processes (PubMed:28167758). The RXR/RAR heterodimers bind to the retinoic acid response elements (RARE) composed of tandem 5'-AGGTCA-3' sites known as DR1-DR5 (PubMed:28167758, PubMed:19398580). In the absence of ligand, the RXR-RAR heterodimers associate with a multiprotein complex containing transcription corepressors that induce histone deacetylation, chromatin condensation and transcriptional suppression (PubMed:16417524). On ligand binding, the corepressors dissociate from the receptors and associate with the coactivators leading to transcriptional activation (PubMed:9267036, PubMed:19850744, PubMed:20215566). Formation of a complex with histone deacetylases might lead to inhibition of RARE DNA element binding and to transcriptional repression (PubMed:28167758). Transcriptional activation and RARE DNA element binding might be supported by the transcription factor KLF2 (PubMed:28167758). RARA plays an essential role in the regulation of retinoic acid-induced germ cell development during spermatogenesis (By similarity). Has a role in the survival of early spermatocytes at the beginning prophase of meiosis (By similarity). In Sertoli cells, may promote the survival and development of early meiotic prophase spermatocytes (By similarity). In concert with RARG, required for skeletal growth, matrix homeostasis and growth plate function (By similarity). Together with RXRA, positively regulates microRNA-10a expression, thereby inhibiting the GATA6/VCAM1 signaling response to pulsatile shear stress in vascular endothelial cells (PubMed:28167758). In association with HDAC3, HDAC5 and HDAC7 corepressors, plays a role in the repression of microRNA-10a and thereby promotes the inflammatory response (PubMed:28167758). {ECO:0000250|UniProtKB:P11416, ECO:0000269|PubMed:16417524, ECO:0000269|PubMed:19398580, ECO:0000269|PubMed:19850744, ECO:0000269|PubMed:20215566, ECO:0000269|PubMed:28167758, ECO:0000269|PubMed:9267036}.

P22692

Main function of 5'-partner protein: FUNCTION: IGF-binding proteins prolong the half-life of the IGFs and have been shown to either inhibit or stimulate the growth promoting effects of the IGFs on cell culture. They alter the interaction of IGFs with their cell surface receptors.
Ensembl transtripts involved in fusion geneENST idsENST00000254066, ENST00000394089, 
ENST00000425707, ENST00000394081, 
ENST00000394086, ENST00000420042, 
ENST00000269593, ENST00000542955, 
Fusion gene scores for assessment (based on all fusion genes of FusionGDB 2.0)* DoF score47 X 22 X 8=827213 X 12 X 5=780
# samples 7416
** MAII scorelog2(74/8272*10)=-3.48263900979862
possibly effective Gene in Pan-Cancer Fusion Genes (peGinPCFGs).
DoF>8 and MAII<0
log2(16/780*10)=-2.28540221886225
possibly effective Gene in Pan-Cancer Fusion Genes (peGinPCFGs).
DoF>8 and MAII<0
Fusion gene context

PubMed: RARA [Title/Abstract] AND IGFBP4 [Title/Abstract] AND fusion [Title/Abstract]

Fusion neoantigen context

PubMed: RARA [Title/Abstract] AND IGFBP4 [Title/Abstract] AND neoantigen [Title/Abstract]

Most frequent breakpoint (based on all fusion genes of FusionGDB 2.0)RARA(38487648)-IGFBP4(38609237), # samples:2
Anticipated loss of major functional domain due to fusion event.RARA-IGFBP4 seems lost the major protein functional domain in Hgene partner, which is a CGC by not retaining the major functional domain in the partially deleted in-frame ORF.
RARA-IGFBP4 seems lost the major protein functional domain in Hgene partner, which is a CGC by not retaining the major functional domain in the partially deleted in-frame ORF.
RARA-IGFBP4 seems lost the major protein functional domain in Hgene partner, which is a essential gene by not retaining the major functional domain in the partially deleted in-frame ORF.
RARA-IGFBP4 seems lost the major protein functional domain in Hgene partner, which is a essential gene by not retaining the major functional domain in the partially deleted in-frame ORF.
* DoF score (Degree of Frequency) = # partners X # break points X # cancer types
** MAII score (Major Active Isofusion Index) = log2(# samples/DoF score*10)

check button Gene ontology of each fusion partner gene with evidence of Inferred from Direct Assay (IDA) from Entrez
PartnerGeneGO IDGO termPubMed ID
HgeneRARA

GO:0007165

signal transduction

2825025

HgeneRARA

GO:0030853

negative regulation of granulocyte differentiation

19917671

HgeneRARA

GO:0032689

negative regulation of interferon-gamma production

18416830

HgeneRARA

GO:0032720

negative regulation of tumor necrosis factor production

18416830

HgeneRARA

GO:0032736

positive regulation of interleukin-13 production

18416830

HgeneRARA

GO:0032753

positive regulation of interleukin-4 production

18416830

HgeneRARA

GO:0032754

positive regulation of interleukin-5 production

18416830

HgeneRARA

GO:0045630

positive regulation of T-helper 2 cell differentiation

18416830

HgeneRARA

GO:0045892

negative regulation of transcription, DNA-templated

20080953

HgeneRARA

GO:0045893

positive regulation of transcription, DNA-templated

18845237|19850744|20080953

HgeneRARA

GO:0045944

positive regulation of transcription by RNA polymerase II

19850744|21131358

HgeneRARA

GO:0071300

cellular response to retinoic acid

19917671

HgeneRARA

GO:0071391

cellular response to estrogen stimulus

20080953



check button Four levels of functional features of fusion genes
Go to FGviewer search page for the most frequent breakpoint (https://ccsmweb.uth.edu/FGviewer/chr17:38487648/chr17:38609237)
- FGviewer provides the online visualization of the retention search of the protein functional features across DNA, RNA, protein, and pathological levels.
- How to search
1. Put your fusion gene symbol.
2. Press the tab key until there will be shown the breakpoint information filled.
4. Go down and press 'Search' tab twice.
4. Go down to have the hyperlink of the search result.
5. Click the hyperlink.
6. See the FGviewer result for your fusion gene.
FGviewer

check buttonRetention analysis results of each fusion partner protein across 39 protein features of UniProt such as six molecule processing features, 13 region features, four site features, six amino acid modification features, two natural variation features, five experimental info features, and 3 secondary structure features, are available here.

check buttonFusion gene breakpoints across RARA (5'-gene)
* Click on the image to open the UCSC genome browser with custom track showing this image in a new window.
all structure

check buttonFusion gene breakpoints across IGFBP4 (3'-gene)
* Click on the image to open the UCSC genome browser with custom track showing this image in a new window.
all structure


Top

Fusion Amino Acid Sequences


check buttonFusion information from ORFfinder translation from full-length transcript sequence from FusionPDB.
HenstTenstHgeneHchrHbpHstrandTgeneTchrTbpTstrandSeq length
(transcript)
BP loci
(transcript)
Predicted start
(transcript)
Predicted stop
(transcript)
Seq length
(amino acids)
ENST00000254066RARAchr1738487648+ENST00000542955IGFBP4chr1738609237+15616334401060206
ENST00000254066RARAchr1738487648+ENST00000269593IGFBP4chr1738609237+22096334401060206
ENST00000394089RARAchr1738487648+ENST00000542955IGFBP4chr1738609237+16367085151135206
ENST00000394089RARAchr1738487648+ENST00000269593IGFBP4chr1738609237+22847085151135206
ENST00000425707RARAchr1738487648+ENST00000542955IGFBP4chr1738609237+16327045111131206
ENST00000425707RARAchr1738487648+ENST00000269593IGFBP4chr1738609237+22807045111131206

check buttonDeepORF prediction of the coding potential based on the fusion transcript sequence of in-frame fusion genes. DeepORF is a coding potential classifier based on convolutional neural network by comparing the real Ribo-seq data. If the no-coding score < 0.5 and coding score > 0.5, then the in-frame fusion transcript is predicted as being likely translated.
HenstTenstHgeneHchrHbpHstrandTgeneTchrTbpTstrandNo-coding scoreCoding score
ENST00000254066ENST00000542955RARAchr1738487648+IGFBP4chr1738609237+0.45485890.5451411
ENST00000254066ENST00000269593RARAchr1738487648+IGFBP4chr1738609237+0.349277560.6507225
ENST00000394089ENST00000542955RARAchr1738487648+IGFBP4chr1738609237+0.415922370.5840776
ENST00000394089ENST00000269593RARAchr1738487648+IGFBP4chr1738609237+0.349632950.6503671
ENST00000425707ENST00000542955RARAchr1738487648+IGFBP4chr1738609237+0.387773160.6122268
ENST00000425707ENST00000269593RARAchr1738487648+IGFBP4chr1738609237+0.308664830.69133514

check button Predicted full-length fusion amino acid sequences. For individual full-length fusion transcript sequence from FusionPDB, we ran ORFfinder and chose the longest ORF among all the predicted ones.

Get the fusion protein sequences from here.

Fusion protein sequence information is available in the fasta format.
>FusionGDB ID_FusionGDB isoform ID_FGname_Hgene_Hchr_Hbp_Henst_Tgene_Tchr_Tbp_Tenst_length(fusion AA) seq_BP

Top

Fusion Protein Breakpoint Sequences for RARA-IGFBP4

check button +/-13 AA sequence from the breakpoints of the fusion protein sequences.
HgeneHchrHbpTgeneTchrTbpLength(fusion protein)BP in fusion proteinPeptide

Top

Potential FusionNeoAntigen Information of RARA-IGFBP4 in HLA I

check button Multiple sequence alignments of the potential FusionNeoAntigens per fusion breakpoints. If the MSA is empty, then it means that there were predicted fusion neoantigens in this fusion breakpoint, but those predicted fusion neoantigens were not across the breakpoint, which is not fusion-specific.

check button Potential FusionNeoAntigen Information
* We used NetMHCpan v4.1 (%rank<0.5) and deepHLApan v1.1 (immunogenic score>0.5)
Fusion geneHchrHbpTgeneTchrTbpHLA IFusionNeoAntigen peptideBinding scoreImmunogenic scoreNeoantigen start (at BP 13)Neoantigen end (at BP 13)

Top

Potential FusionNeoAntigen Information of RARA-IGFBP4 in HLA II

check button Multiple sequence alignments of the potential FusionNeoAntigens per fusion breakpoints. If the MSA is empty, then it means that there were predicted fusion neoantigens in this fusion breakpoint, but those predicted fusion neoantigens were not across the breakpoint, which is not fusion-specific.

check button Potential FusionNeoAntigen Information
* We used NetMHCIIpan v4.1 (%rank<0.5).
Fusion geneHchrHbpTgeneTchrTbpHLA IIFusionNeoAntigen peptideNeoantigen start (at BP 13)Neoantigen end (at BP 13)

Top

Fusion breakpoint peptide structures of RARA-IGFBP4

check button3D structures of the fusion breakpoint peptide of 14AA sequence that have potential fusion neoantigens
* The minimum length of the amino acid sequence in RoseTTAFold is 14AA. Here, we predicted the 14AA fusion protein breakpoint sequence not the fusion neoantigen peptide, which is shorter than 14 AA.

Top

Filtering FusionNeoAntigens Through Checking the Interaction with HLAs in 3D of RARA-IGFBP4

check buttonVirtual screening between 25 HLAs (from PDB) and FusionNeoAntigens
* We used Glide to predict the interaction between HLAs and neoantigens.
HLA allelePDB IDFile nameBPseqDocking scoreGlide score

Top

Vaccine Design for the FusionNeoAntigens of RARA-IGFBP4

check button mRNA and peptide sequences of FusionNeoAntigens that have potential interaction with HLA-Is.
Fusion geneHchrHbpTchrTbpStart in +/-13AAEnd in +/-13AAFusionNeoAntigen peptide sequenceFusionNeoAntigen RNA sequence

check button mRNA and peptide sequences of FusionNeoAntigens that have potential interaction with HLA-IIs.
Fusion geneHchrHbpTchrTbpStart in +/-13AAEnd in +/-13AAFusionNeoAntigen peptideFusionNEoAntigen RNA sequence

Top

Information of the samples that have these potential fusion neoantigens of RARA-IGFBP4

check button These samples were reported as having these fusion breakpoints. For individual breakpoints, we checked the open reading frames considering multiple gene isoforms and chose the in-frame fusion genes only. Then, we made fusion protein sequences and predicted the fusion neoantigens. These fusion-positive samples may have these potential fusion neoantigens.
Cancer typeFusion geneHchrHbpHenstTchrTbpTenstSample

Top

Potential target of CAR-T therapy development for RARA-IGFBP4

check button Predicted 3D structure. We used RoseTTAFold.

check buttonRetention analysis result of each fusion partner protein across 39 protein features of UniProt such as six molecule processing features, 13 region features, four site features, six amino acid modification features, two natural variation features, five experimental info features, and 3 secondary structure features. Here, to provide the retention of the transmembrane domain, we only show the protein feature retention information of those transmembrane features


* Minus value of BPloci means that the break point is located before the CDS.
- In-frame and retained 'Transmembrane'.
PartnerGeneHbpTbpENSTStrandBPexonTotalExonProtein feature loci*BPlociTotalLenProtein featureProtein feature note

check button Subcellular localization prediction of the transmembrane domain retained fusion proteins
* We used DeepLoc 1.0. The order of the X-axis of the barplot is as follows: Entry_ID, Localization, Type, Nucleus, Cytoplasm, Extracellular, Mitochondrion, Cell_membrane, Endoplasmic_reticulum, Plastid, Golgi.apparatus, Lysosome.Vacuole, Peroxisome. Y-axis is the output score of DeepLoc. Clicking the image will open a new tab with a large image.
HgeneHchrHbpHenstTgeneTchrTbpTenstDeepLoc result

Top

Related Drugs to RARA-IGFBP4

check button Drugs used for this fusion-positive patient.
(Manual curation of PubMed, 04-30-2022 + MyCancerGenome)
HgeneTgeneDrugSourcePMID

Top

Related Diseases to RARA-IGFBP4

check button Diseases that have this fusion gene.
(Manual curation of PubMed, 04-30-2022 + MyCancerGenome)
HgeneTgeneDiseaseSourcePMID

check button Diseases associated with fusion partners.
(DisGeNet 4.0)
PartnerGeneDisease IDDisease name# pubmedsSource
HgeneRARAC0023487Acute Promyelocytic Leukemia24CTD_human;ORPHANET
HgeneRARAC0036341Schizophrenia3PSYGENET
HgeneRARAC0006142Malignant neoplasm of breast1CTD_human
HgeneRARAC0009363Congenital ocular coloboma (disorder)1GENOMICS_ENGLAND
HgeneRARAC0010701Phyllodes Tumor1CTD_human
HgeneRARAC0085183Neoplasms, Second Primary1CTD_human
HgeneRARAC0086696Neoplasms, Therapy-Associated1CTD_human
HgeneRARAC0149940Sciatic Neuropathy1CTD_human
HgeneRARAC0154748Lesion of Sciatic Nerve1CTD_human
HgeneRARAC0206650Fibroadenoma1CTD_human
HgeneRARAC0242013Sciatic Neuritis1CTD_human
HgeneRARAC0525045Mood Disorders1PSYGENET
HgeneRARAC0600066Malignant Cystosarcoma Phyllodes1CTD_human
HgeneRARAC0678222Breast Carcinoma1CTD_human
HgeneRARAC0751924Neuralgia-Neuritis, Sciatic Nerve1CTD_human
HgeneRARAC0751925Sciatic Nerve Palsy1CTD_human
HgeneRARAC0877578Treatment related secondary malignancy1CTD_human
HgeneRARAC1257931Mammary Neoplasms, Human1CTD_human
HgeneRARAC1458155Mammary Neoplasms1CTD_human
HgeneRARAC2239176Liver carcinoma1CTD_human
HgeneRARAC4704874Mammary Carcinoma, Human1CTD_human