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:ESR1-RFX4

Fusion Gene and Fusion Protein Summary

check button Fusion gene summary
Fusion partner gene informationFusion gene name: ESR1-RFX4
FusionPDB ID: 27550
FusionGDB2.0 ID: 27550
HgeneTgene
Gene symbol

ESR1

RFX4

Gene ID

2099

5992

Gene nameestrogen receptor 1regulatory factor X4
SynonymsER|ESR|ESRA|ESTRR|Era|NR3A1NYD-SP10
Cytomap

6q25.1-q25.2

12q23.3

Type of geneprotein-codingprotein-coding
Descriptionestrogen receptorE2 receptor alphaER-alphaestradiol receptorestrogen nuclear receptor alphaestrogen receptor alpha E1-E2-1-2estrogen receptor alpha E1-N2-E2-1-2nuclear receptor subfamily 3 group A member 1oestrogen receptor alphatranscription factor RFX4regulatory factor X, 4 (influences HLA class II expression)testis development protein NYD-SP10winged-helix transcription factor RFX4
Modification date2020032920200313
UniProtAcc

P03372

Main function of 5'-partner protein: FUNCTION: Nuclear hormone receptor. The steroid hormones and their receptors are involved in the regulation of eukaryotic gene expression and affect cellular proliferation and differentiation in target tissues. Ligand-dependent nuclear transactivation involves either direct homodimer binding to a palindromic estrogen response element (ERE) sequence or association with other DNA-binding transcription factors, such as AP-1/c-Jun, c-Fos, ATF-2, Sp1 and Sp3, to mediate ERE-independent signaling. Ligand binding induces a conformational change allowing subsequent or combinatorial association with multiprotein coactivator complexes through LXXLL motifs of their respective components. Mutual transrepression occurs between the estrogen receptor (ER) and NF-kappa-B in a cell-type specific manner. Decreases NF-kappa-B DNA-binding activity and inhibits NF-kappa-B-mediated transcription from the IL6 promoter and displace RELA/p65 and associated coregulators from the promoter. Recruited to the NF-kappa-B response element of the CCL2 and IL8 promoters and can displace CREBBP. Present with NF-kappa-B components RELA/p65 and NFKB1/p50 on ERE sequences. Can also act synergistically with NF-kappa-B to activate transcription involving respective recruitment adjacent response elements; the function involves CREBBP. Can activate the transcriptional activity of TFF1. Also mediates membrane-initiated estrogen signaling involving various kinase cascades. Isoform 3 is involved in activation of NOS3 and endothelial nitric oxide production. Isoforms lacking one or several functional domains are thought to modulate transcriptional activity by competitive ligand or DNA binding and/or heterodimerization with the full-length receptor. Essential for MTA1-mediated transcriptional regulation of BRCA1 and BCAS3. Isoform 3 can bind to ERE and inhibit isoform 1. {ECO:0000269|PubMed:10681512, ECO:0000269|PubMed:10816575, ECO:0000269|PubMed:11477071, ECO:0000269|PubMed:11682626, ECO:0000269|PubMed:14764652, ECO:0000269|PubMed:15078875, ECO:0000269|PubMed:15891768, ECO:0000269|PubMed:16043358, ECO:0000269|PubMed:16617102, ECO:0000269|PubMed:16684779, ECO:0000269|PubMed:17922032, ECO:0000269|PubMed:17932106, ECO:0000269|PubMed:18247370, ECO:0000269|PubMed:19350539, ECO:0000269|PubMed:20074560, ECO:0000269|PubMed:20705611, ECO:0000269|PubMed:21330404, ECO:0000269|PubMed:22083956, ECO:0000269|PubMed:7651415, ECO:0000269|PubMed:9328340}.
.
Ensembl transtripts involved in fusion geneENST idsENST00000206249, ENST00000338799, 
ENST00000427531, ENST00000440973, 
ENST00000443427, ENST00000456483, 
ENST00000406599, ENST00000482101, 
ENST00000229387, ENST00000357881, 
ENST00000392842, 
Fusion gene scores for assessment (based on all fusion genes of FusionGDB 2.0)* DoF score21 X 18 X 8=30247 X 7 X 5=245
# samples 3510
** MAII scorelog2(35/3024*10)=-3.11103131238874
possibly effective Gene in Pan-Cancer Fusion Genes (peGinPCFGs).
DoF>8 and MAII<0
log2(10/245*10)=-1.29278174922785
possibly effective Gene in Pan-Cancer Fusion Genes (peGinPCFGs).
DoF>8 and MAII<0
Fusion gene context

PubMed: ESR1 [Title/Abstract] AND RFX4 [Title/Abstract] AND fusion [Title/Abstract]

Fusion neoantigen context

PubMed: ESR1 [Title/Abstract] AND RFX4 [Title/Abstract] AND neoantigen [Title/Abstract]

Most frequent breakpoint (based on all fusion genes of FusionGDB 2.0)ESR1(152201906)-RFX4(107033109), # samples:1
Anticipated loss of major functional domain due to fusion event.ESR1-RFX4 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.
ESR1-RFX4 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.
ESR1-RFX4 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.
ESR1-RFX4 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
HgeneESR1

GO:0006366

transcription by RNA polymerase II

15831516

HgeneESR1

GO:0010629

negative regulation of gene expression

21695196

HgeneESR1

GO:0030520

intracellular estrogen receptor signaling pathway

9841876

HgeneESR1

GO:0032355

response to estradiol

15304487

HgeneESR1

GO:0043124

negative regulation of I-kappaB kinase/NF-kappaB signaling

7651415|16043358

HgeneESR1

GO:0043433

negative regulation of DNA-binding transcription factor activity

10816575

HgeneESR1

GO:0043627

response to estrogen

11581164

HgeneESR1

GO:0045893

positive regulation of transcription, DNA-templated

9841876|20074560

HgeneESR1

GO:0045899

positive regulation of RNA polymerase II transcriptional preinitiation complex assembly

9841876

HgeneESR1

GO:0045944

positive regulation of transcription by RNA polymerase II

11544182|12047722|15345745|15831516|18563714

HgeneESR1

GO:0051091

positive regulation of DNA-binding transcription factor activity

9328340|10681512

HgeneESR1

GO:0071392

cellular response to estradiol stimulus

15831516

TgeneRFX4

GO:0045944

positive regulation of transcription by RNA polymerase II

16893423



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/chr6:152201906/chr12:107033109)
- 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 ESR1 (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 RFX4 (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)
ENST00000440973ESR1chr6152201906+ENST00000392842RFX4chr12107033109+454111303703207945
ENST00000440973ESR1chr6152201906+ENST00000357881RFX4chr12107033109+339811303703207945
ENST00000338799ESR1chr6152201906+ENST00000392842RFX4chr12107033109+5631222014604297945
ENST00000338799ESR1chr6152201906+ENST00000357881RFX4chr12107033109+4488222014604297945
ENST00000456483ESR1chr6152201906+ENST00000392842RFX4chr12107033109+43899782183055945
ENST00000456483ESR1chr6152201906+ENST00000357881RFX4chr12107033109+32469782183055945
ENST00000443427ESR1chr6152201906+ENST00000392842RFX4chr12107033109+443210212613098945
ENST00000443427ESR1chr6152201906+ENST00000357881RFX4chr12107033109+328910212613098945
ENST00000206249ESR1chr6152201906+ENST00000392842RFX4chr12107033109+453311223623199945
ENST00000206249ESR1chr6152201906+ENST00000357881RFX4chr12107033109+339011223623199945
ENST00000427531ESR1chr6152201906+ENST00000392842RFX4chr12107033109+38844732052550781
ENST00000427531ESR1chr6152201906+ENST00000357881RFX4chr12107033109+27414732052550781

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
ENST00000440973ENST00000392842ESR1chr6152201906+RFX4chr12107033109+0.0005227310.9994773
ENST00000440973ENST00000357881ESR1chr6152201906+RFX4chr12107033109+0.0023977420.9976023
ENST00000338799ENST00000392842ESR1chr6152201906+RFX4chr12107033109+0.0011878430.9988122
ENST00000338799ENST00000357881ESR1chr6152201906+RFX4chr12107033109+0.0032978940.99670213
ENST00000456483ENST00000392842ESR1chr6152201906+RFX4chr12107033109+0.0004695670.99953043
ENST00000456483ENST00000357881ESR1chr6152201906+RFX4chr12107033109+0.002360650.9976394
ENST00000443427ENST00000392842ESR1chr6152201906+RFX4chr12107033109+0.0005022910.99949765
ENST00000443427ENST00000357881ESR1chr6152201906+RFX4chr12107033109+0.0024858770.9975141
ENST00000206249ENST00000392842ESR1chr6152201906+RFX4chr12107033109+0.0004543610.99954563
ENST00000206249ENST00000357881ESR1chr6152201906+RFX4chr12107033109+0.0020773040.9979227
ENST00000427531ENST00000392842ESR1chr6152201906+RFX4chr12107033109+0.0013579550.998642
ENST00000427531ENST00000357881ESR1chr6152201906+RFX4chr12107033109+0.0056825140.9943175

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 ESR1-RFX4

check button +/-13 AA sequence from the breakpoints of the fusion protein sequences.
HgeneHchrHbpTgeneTchrTbpLength(fusion protein)BP in fusion proteinPeptide
ESR1chr6152201906RFX4chr121070331091021253RLRKCYEVGMMKGDEEKENNRASKPH
ESR1chr6152201906RFX4chr121070331091122253RLRKCYEVGMMKGDEEKENNRASKPH
ESR1chr6152201906RFX4chr121070331091130253RLRKCYEVGMMKGDEEKENNRASKPH
ESR1chr6152201906RFX4chr121070331092220253RLRKCYEVGMMKGDEEKENNRASKPH
ESR1chr6152201906RFX4chr1210703310947389RLRKCYEVGMMKGDEEKENNRASKPH
ESR1chr6152201906RFX4chr12107033109978253RLRKCYEVGMMKGDEEKENNRASKPH

Top

Potential FusionNeoAntigen Information of ESR1-RFX4 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 ESR1-RFX4 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 ESR1-RFX4

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 ESR1-RFX4

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 ESR1-RFX4

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 ESR1-RFX4

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 ESR1-RFX4

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 ESR1-RFX4

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

Top

Related Diseases to ESR1-RFX4

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