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:WDR33-BACE2

Fusion Gene and Fusion Protein Summary

check button Fusion gene summary
Fusion partner gene informationFusion gene name: WDR33-BACE2
FusionPDB ID: 98932
FusionGDB2.0 ID: 98932
HgeneTgene
Gene symbol

WDR33

BACE2

Gene ID

55339

25825

Gene nameWD repeat domain 33beta-secretase 2
SynonymsNET14|WDC146AEPLC|ALP56|ASP1|ASP21|BAE2|CDA13|CEAP1|DRAP
Cytomap

2q14.3

21q22.2-q22.3

Type of geneprotein-codingprotein-coding
Descriptionpre-mRNA 3' end processing protein WDR33WD repeat-containing protein 33WD repeat-containing protein WDC146WD repeat-containing protein of 146 kDabeta-secretase 256 kDa aspartic-like proteaseDown syndrome region aspartic proteaseSLCO3A1/BACE2 fusionaspartyl protease 1beta-site APP-cleaving enzyme 2beta-site amyloid beta A4 precursor protein-cleaving enzyme 2memapsin-1membrane-associated asp
Modification date2020031320200313
UniProtAcc.

Q9Y5Z0

Main function of 5'-partner protein: FUNCTION: Responsible for the proteolytic processing of the amyloid precursor protein (APP). Cleaves APP, between residues 690 and 691, leading to the generation and extracellular release of beta-cleaved soluble APP, and a corresponding cell-associated C-terminal fragment which is later released by gamma-secretase. It has also been shown that it can cleave APP between residues 671 and 672. Responsible also for the proteolytic processing of CLTRN in pancreatic beta cells (PubMed:21907142). {ECO:0000269|PubMed:10591213, ECO:0000269|PubMed:11083922, ECO:0000269|PubMed:11423558, ECO:0000269|PubMed:15857888, ECO:0000269|PubMed:16816112, ECO:0000269|PubMed:21907142}.
Ensembl transtripts involved in fusion geneENST idsENST00000322313, ENST00000393006, 
ENST00000409658, 
ENST00000466122, 
ENST00000328735, ENST00000347667, 
ENST00000330333, 
Fusion gene scores for assessment (based on all fusion genes of FusionGDB 2.0)* DoF score10 X 9 X 4=36014 X 9 X 8=1008
# samples 1014
** MAII scorelog2(10/360*10)=-1.84799690655495
possibly effective Gene in Pan-Cancer Fusion Genes (peGinPCFGs).
DoF>8 and MAII<0
log2(14/1008*10)=-2.84799690655495
possibly effective Gene in Pan-Cancer Fusion Genes (peGinPCFGs).
DoF>8 and MAII<0
Fusion gene context

PubMed: WDR33 [Title/Abstract] AND BACE2 [Title/Abstract] AND fusion [Title/Abstract]

Fusion neoantigen context

PubMed: WDR33 [Title/Abstract] AND BACE2 [Title/Abstract] AND neoantigen [Title/Abstract]

Most frequent breakpoint (based on all fusion genes of FusionGDB 2.0)WDR33(128520636)-BACE2(42647298), # samples:1
Anticipated loss of major functional domain due to fusion event.WDR33-BACE2 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.
WDR33-BACE2 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.
WDR33-BACE2 seems lost the major protein functional domain in Hgene partner, which is a essential gene due to the frame-shifted ORF.
WDR33-BACE2 seems lost the major protein functional domain in Tgene partner, which is a IUPHAR drug target due to the frame-shifted 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
TgeneBACE2

GO:0006509

membrane protein ectodomain proteolysis

10591213



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/chr2:128520636/chr21:42647298)
- 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 WDR33 (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 BACE2 (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)
ENST00000322313WDR33chr2128520636-ENST00000330333BACE2chr2142647298+21108831591136325
ENST00000393006WDR33chr2128520636-ENST00000330333BACE2chr2142647298+21349071831160325

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
ENST00000322313ENST00000330333WDR33chr2128520636-BACE2chr2142647298+0.0017758740.9982241
ENST00000393006ENST00000330333WDR33chr2128520636-BACE2chr2142647298+0.0017882210.99821174

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 WDR33-BACE2

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

Top

Potential FusionNeoAntigen Information of WDR33-BACE2 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.
WDR33-BACE2_128520636_42647298.msa

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)
WDR33-BACE2chr2128520636chr2142647298883HLA-B41:01REIAGAAV0.94810.95681220
WDR33-BACE2chr2128520636chr2142647298883HLA-B39:13REIAGAAV0.93430.97971220
WDR33-BACE2chr2128520636chr2142647298883HLA-B40:05REIAGAAV0.92670.51281220
WDR33-BACE2chr2128520636chr2142647298883HLA-B50:01REIAGAAV0.84020.95391220
WDR33-BACE2chr2128520636chr2142647298883HLA-B45:01EERILREIA0.99750.6248716
WDR33-BACE2chr2128520636chr2142647298883HLA-B45:01HEERILREI0.99020.5988615
WDR33-BACE2chr2128520636chr2142647298883HLA-B44:03HEERILREI0.9870.9438615
WDR33-BACE2chr2128520636chr2142647298883HLA-B45:01REIAGAAVS0.96380.9421221
WDR33-BACE2chr2128520636chr2142647298883HLA-B50:02REIAGAAVS0.9540.74831221
WDR33-BACE2chr2128520636chr2142647298883HLA-B41:01HEERILREI0.91490.7855615
WDR33-BACE2chr2128520636chr2142647298883HLA-B39:13HEERILREI0.80770.8704615
WDR33-BACE2chr2128520636chr2142647298883HLA-B41:01REIAGAAVS0.2840.97721221
WDR33-BACE2chr2128520636chr2142647298883HLA-B50:01REIAGAAVS0.12360.85781221
WDR33-BACE2chr2128520636chr2142647298883HLA-B40:06REIAGAAV0.99990.81221220
WDR33-BACE2chr2128520636chr2142647298883HLA-B40:06REIAGAAVS0.99130.78751221
WDR33-BACE2chr2128520636chr2142647298883HLA-B39:08HEERILREI0.92920.7229615
WDR33-BACE2chr2128520636chr2142647298883HLA-B73:01ERILREIAGA0.99690.8225818
WDR33-BACE2chr2128520636chr2142647298883HLA-B73:01ERILREIAGAA0.99910.8395819
WDR33-BACE2chr2128520636chr2142647298883HLA-B50:05REIAGAAV0.84020.95391220
WDR33-BACE2chr2128520636chr2142647298883HLA-B50:04REIAGAAV0.84020.95391220
WDR33-BACE2chr2128520636chr2142647298883HLA-B40:04HEERILREI0.99850.5515615
WDR33-BACE2chr2128520636chr2142647298883HLA-B44:13HEERILREI0.9870.9438615
WDR33-BACE2chr2128520636chr2142647298883HLA-B44:07HEERILREI0.9870.9438615
WDR33-BACE2chr2128520636chr2142647298883HLA-B44:26HEERILREI0.9870.9438615
WDR33-BACE2chr2128520636chr2142647298883HLA-B50:04REIAGAAVS0.12360.85781221
WDR33-BACE2chr2128520636chr2142647298883HLA-B50:05REIAGAAVS0.12360.85781221

Top

Potential FusionNeoAntigen Information of WDR33-BACE2 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.
WDR33-BACE2_128520636_42647298.msa

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)
WDR33-BACE2chr2128520636chr2142647298883DRB1-0102EERILREIAGAAVSE722
WDR33-BACE2chr2128520636chr2142647298883DRB1-0102ERILREIAGAAVSEI823
WDR33-BACE2chr2128520636chr2142647298883DRB1-0123EERILREIAGAAVSE722
WDR33-BACE2chr2128520636chr2142647298883DRB1-0123ERILREIAGAAVSEI823
WDR33-BACE2chr2128520636chr2142647298883DRB1-0901LREIAGAAVSEISGP1126
WDR33-BACE2chr2128520636chr2142647298883DRB1-0904LREIAGAAVSEISGP1126
WDR33-BACE2chr2128520636chr2142647298883DRB1-0906LREIAGAAVSEISGP1126
WDR33-BACE2chr2128520636chr2142647298883DRB1-0906ILREIAGAAVSEISG1025
WDR33-BACE2chr2128520636chr2142647298883DRB1-0909LREIAGAAVSEISGP1126
WDR33-BACE2chr2128520636chr2142647298883DRB1-1002EERILREIAGAAVSE722
WDR33-BACE2chr2128520636chr2142647298883DRB1-1410LREIAGAAVSEISGP1126
WDR33-BACE2chr2128520636chr2142647298883DRB1-1615EERILREIAGAAVSE722
WDR33-BACE2chr2128520636chr2142647298883DRB1-1615ERILREIAGAAVSEI823

Top

Fusion breakpoint peptide structures of WDR33-BACE2

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.
File nameBPseqHgeneTgeneHchrHbpTchrTbpAAlen
3299HEERILREIAGAAVWDR33BACE2chr2128520636chr2142647298883

Top

Filtering FusionNeoAntigens Through Checking the Interaction with HLAs in 3D of WDR33-BACE2

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
HLA-B14:023BVN3299HEERILREIAGAAV-7.15543-7.26883
HLA-B14:023BVN3299HEERILREIAGAAV-4.77435-5.80965
HLA-B52:013W393299HEERILREIAGAAV-6.80875-6.92215
HLA-B52:013W393299HEERILREIAGAAV-4.20386-5.23916
HLA-A11:014UQ23299HEERILREIAGAAV-7.5194-8.5547
HLA-A11:014UQ23299HEERILREIAGAAV-6.9601-7.0735
HLA-A24:025HGA3299HEERILREIAGAAV-7.52403-7.63743
HLA-A24:025HGA3299HEERILREIAGAAV-5.82433-6.85963
HLA-B27:056PYJ3299HEERILREIAGAAV-3.28285-4.31815
HLA-B44:053DX83299HEERILREIAGAAV-5.91172-6.94702
HLA-B44:053DX83299HEERILREIAGAAV-4.24346-4.35686

Top

Vaccine Design for the FusionNeoAntigens of WDR33-BACE2

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
WDR33-BACE2chr2128520636chr21426472981220REIAGAAVGAGAAATTGCAGGTGCTGCAGTGT
WDR33-BACE2chr2128520636chr21426472981221REIAGAAVSGAGAAATTGCAGGTGCTGCAGTGTCTG
WDR33-BACE2chr2128520636chr2142647298615HEERILREIATGAGGAAAGAATTCTCCGAGAAATTG
WDR33-BACE2chr2128520636chr2142647298716EERILREIAAGGAAAGAATTCTCCGAGAAATTGCAG
WDR33-BACE2chr2128520636chr2142647298818ERILREIAGAAAAGAATTCTCCGAGAAATTGCAGGTGCTG
WDR33-BACE2chr2128520636chr2142647298819ERILREIAGAAAAAGAATTCTCCGAGAAATTGCAGGTGCTGCAG

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

Top

Information of the samples that have these potential fusion neoantigens of WDR33-BACE2

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
STADWDR33-BACE2chr2128520636ENST00000322313chr2142647298ENST00000330333TCGA-CD-8525-01A

Top

Potential target of CAR-T therapy development for WDR33-BACE2

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
TgeneBACE2chr2:128520636chr21:42647298ENST0000032873568474_4940397.0TransmembraneHelical
TgeneBACE2chr2:128520636chr21:42647298ENST0000033033379474_4940519.0TransmembraneHelical
TgeneBACE2chr2:128520636chr21:42647298ENST0000034766768474_4940469.0TransmembraneHelical

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 WDR33-BACE2

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

Top

Related Diseases to WDR33-BACE2

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