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:VRK3-TGFB1

Fusion Gene and Fusion Protein Summary

check button Fusion gene summary
Fusion partner gene informationFusion gene name: VRK3-TGFB1
FusionPDB ID: 98549
FusionGDB2.0 ID: 98549
HgeneTgene
Gene symbol

VRK3

TGFB1

Gene ID

51231

7040

Gene nameVRK serine/threonine kinase 3transforming growth factor beta 1
Synonyms-CED|DPD1|IBDIMDE|LAP|TGF-beta1|TGFB|TGFbeta
Cytomap

19q13.33

19q13.2

Type of geneprotein-codingprotein-coding
Descriptioninactive serine/threonine-protein kinase VRK3serine/threonine-protein kinase VRK3serine/threonine-protein pseudokinase VRK3vaccinia related kinase 3transforming growth factor beta-1 proproteinTGF-beta-1latency-associated peptideprepro-transforming growth factor beta-1transforming growth factor beta1
Modification date2020031320200329
UniProtAcc.

TIAF1

Main function of 5'-partner protein: 115
Ensembl transtripts involved in fusion geneENST idsENST00000424804, ENST00000316763, 
ENST00000377011, ENST00000593919, 
ENST00000594948, ENST00000601341, 
ENST00000601912, ENST00000443401, 
ENST00000594092, ENST00000599538, 
ENST00000221930, 
Fusion gene scores for assessment (based on all fusion genes of FusionGDB 2.0)* DoF score7 X 6 X 4=1688 X 6 X 7=336
# samples 98
** MAII scorelog2(9/168*10)=-0.900464326449086
possibly effective Gene in Pan-Cancer Fusion Genes (peGinPCFGs).
DoF>8 and MAII<0
log2(8/336*10)=-2.0703893278914
possibly effective Gene in Pan-Cancer Fusion Genes (peGinPCFGs).
DoF>8 and MAII<0
Fusion gene context

PubMed: VRK3 [Title/Abstract] AND TGFB1 [Title/Abstract] AND fusion [Title/Abstract]

Fusion neoantigen context

PubMed: VRK3 [Title/Abstract] AND TGFB1 [Title/Abstract] AND neoantigen [Title/Abstract]

Most frequent breakpoint (based on all fusion genes of FusionGDB 2.0)VRK3(50504047)-TGFB1(41838186), # samples:2
Anticipated loss of major functional domain due to fusion event.VRK3-TGFB1 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.
VRK3-TGFB1 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.
VRK3-TGFB1 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.
VRK3-TGFB1 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.
VRK3-TGFB1 seems lost the major protein functional domain in Hgene partner, which is a IUPHAR drug target due to the frame-shifted ORF.
VRK3-TGFB1 seems lost the major protein functional domain in Hgene partner, which is a kinase due to the frame-shifted ORF.
VRK3-TGFB1 seems lost the major protein functional domain in Tgene partner, which is a tumor suppressor 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
TgeneTGFB1

GO:0001837

epithelial to mesenchymal transition

25893292|29529050

TgeneTGFB1

GO:0001933

negative regulation of protein phosphorylation

8053900

TgeneTGFB1

GO:0001934

positive regulation of protein phosphorylation

18625725

TgeneTGFB1

GO:0002062

chondrocyte differentiation

15040835

TgeneTGFB1

GO:0002244

hematopoietic progenitor cell differentiation

15451575

TgeneTGFB1

GO:0006611

protein export from nucleus

9770491|17438144|18588859

TgeneTGFB1

GO:0006754

ATP biosynthetic process

10513816

TgeneTGFB1

GO:0006796

phosphate-containing compound metabolic process

10513816

TgeneTGFB1

GO:0006954

inflammatory response

21147091

TgeneTGFB1

GO:0007050

cell cycle arrest

14555988

TgeneTGFB1

GO:0007093

mitotic cell cycle checkpoint

15334054

TgeneTGFB1

GO:0007173

epidermal growth factor receptor signaling pathway

18625725

TgeneTGFB1

GO:0007179

transforming growth factor beta receptor signaling pathway

9389648|11157754

TgeneTGFB1

GO:0007182

common-partner SMAD protein phosphorylation

20573232

TgeneTGFB1

GO:0007183

SMAD protein complex assembly

17438144

TgeneTGFB1

GO:0008284

positive regulation of cell proliferation

10513816|14633705

TgeneTGFB1

GO:0008285

negative regulation of cell proliferation

15334054

TgeneTGFB1

GO:0010628

positive regulation of gene expression

18625725|18832382|18941241|19913496|25322725|26634652|26687115|27162619|29167509

TgeneTGFB1

GO:0010629

negative regulation of gene expression

19913496|20067797|22269326|25163461|26634652|29167509|29529050

TgeneTGFB1

GO:0010718

positive regulation of epithelial to mesenchymal transition

17999987|18505915

TgeneTGFB1

GO:0010763

positive regulation of fibroblast migration

18555217

TgeneTGFB1

GO:0010800

positive regulation of peptidyl-threonine phosphorylation

18625725|19736306

TgeneTGFB1

GO:0010862

positive regulation of pathway-restricted SMAD protein phosphorylation

9389648|19736306|26634652

TgeneTGFB1

GO:0010936

negative regulation of macrophage cytokine production

20875417

TgeneTGFB1

GO:0016477

cell migration

25893292

TgeneTGFB1

GO:0017015

regulation of transforming growth factor beta receptor signaling pathway

15334054

TgeneTGFB1

GO:0022408

negative regulation of cell-cell adhesion

18593713

TgeneTGFB1

GO:0030214

hyaluronan catabolic process

17324121

TgeneTGFB1

GO:0030308

negative regulation of cell growth

15334054

TgeneTGFB1

GO:0030335

positive regulation of cell migration

19736306

TgeneTGFB1

GO:0031293

membrane protein intracellular domain proteolysis

25310401

TgeneTGFB1

GO:0031334

positive regulation of protein complex assembly

19366691

TgeneTGFB1

GO:0031663

lipopolysaccharide-mediated signaling pathway

21147091

TgeneTGFB1

GO:0032270

positive regulation of cellular protein metabolic process

15219857

TgeneTGFB1

GO:0032355

response to estradiol

18039789

TgeneTGFB1

GO:0032570

response to progesterone

18039789

TgeneTGFB1

GO:0032740

positive regulation of interleukin-17 production

18453574

TgeneTGFB1

GO:0032801

receptor catabolic process

17878231

TgeneTGFB1

GO:0032930

positive regulation of superoxide anion generation

22073128

TgeneTGFB1

GO:0032967

positive regulation of collagen biosynthetic process

19734317|22269326|25310401

TgeneTGFB1

GO:0033138

positive regulation of peptidyl-serine phosphorylation

18625725|19736306

TgeneTGFB1

GO:0035307

positive regulation of protein dephosphorylation

14555988

TgeneTGFB1

GO:0042307

positive regulation of protein import into nucleus

19366691

TgeneTGFB1

GO:0043117

positive regulation of vascular permeability

21168935

TgeneTGFB1

GO:0043406

positive regulation of MAP kinase activity

18625725

TgeneTGFB1

GO:0043536

positive regulation of blood vessel endothelial cell migration

18555217

TgeneTGFB1

GO:0043537

negative regulation of blood vessel endothelial cell migration

18555217

TgeneTGFB1

GO:0043552

positive regulation of phosphatidylinositol 3-kinase activity

18625725

TgeneTGFB1

GO:0045216

cell-cell junction organization

18505915

TgeneTGFB1

GO:0045599

negative regulation of fat cell differentiation

15040835

TgeneTGFB1

GO:0045662

negative regulation of myoblast differentiation

9770491

TgeneTGFB1

GO:0045786

negative regulation of cell cycle

11502704

TgeneTGFB1

GO:0045892

negative regulation of transcription, DNA-templated

15702480|18832382

TgeneTGFB1

GO:0045893

positive regulation of transcription, DNA-templated

9389648|14517293|15334054|16816361

TgeneTGFB1

GO:0045918

negative regulation of cytolysis

24586048

TgeneTGFB1

GO:0045930

negative regulation of mitotic cell cycle

14555988

TgeneTGFB1

GO:0045944

positive regulation of transcription by RNA polymerase II

18832382

TgeneTGFB1

GO:0048298

positive regulation of isotype switching to IgA isotypes

14988498

TgeneTGFB1

GO:0048642

negative regulation of skeletal muscle tissue development

9770491

TgeneTGFB1

GO:0050680

negative regulation of epithelial cell proliferation

9950587

TgeneTGFB1

GO:0050714

positive regulation of protein secretion

18505915

TgeneTGFB1

GO:0050731

positive regulation of peptidyl-tyrosine phosphorylation

21168935

TgeneTGFB1

GO:0050921

positive regulation of chemotaxis

18555217

TgeneTGFB1

GO:0051897

positive regulation of protein kinase B signaling

18625725

TgeneTGFB1

GO:0060389

pathway-restricted SMAD protein phosphorylation

11157754|17999987|18453574|25893292

TgeneTGFB1

GO:0060390

regulation of SMAD protein signal transduction

25893292

TgeneTGFB1

GO:0060391

positive regulation of SMAD protein signal transduction

9389648|19366691|29167509

TgeneTGFB1

GO:0070168

negative regulation of biomineral tissue development

26634652

TgeneTGFB1

GO:0070374

positive regulation of ERK1 and ERK2 cascade

25310401

TgeneTGFB1

GO:0070723

response to cholesterol

17878231

TgeneTGFB1

GO:0071407

cellular response to organic cyclic compound

21147091

TgeneTGFB1

GO:0071560

cellular response to transforming growth factor beta stimulus

19736306|22269326

TgeneTGFB1

GO:0085029

extracellular matrix assembly

19734317

TgeneTGFB1

GO:0090263

positive regulation of canonical Wnt signaling pathway

12893825|15040835

TgeneTGFB1

GO:0097191

extrinsic apoptotic signaling pathway

15334054

TgeneTGFB1

GO:1900126

negative regulation of hyaluronan biosynthetic process

17324121

TgeneTGFB1

GO:1900182

positive regulation of protein localization to nucleus

26634652

TgeneTGFB1

GO:1901666

positive regulation of NAD+ ADP-ribosyltransferase activity

22073128

TgeneTGFB1

GO:1902895

positive regulation of pri-miRNA transcription by RNA polymerase II

26311719|26493107

TgeneTGFB1

GO:1903077

negative regulation of protein localization to plasma membrane

21168935|24586048

TgeneTGFB1

GO:1903800

positive regulation of production of miRNAs involved in gene silencing by miRNA

18548003

TgeneTGFB1

GO:2000679

positive regulation of transcription regulatory region DNA binding

22073128

TgeneTGFB1

GO:2000727

positive regulation of cardiac muscle cell differentiation

25163461



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/chr19:50504047/chr19:41838186)
- 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 VRK3 (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 TGFB1 (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)
ENST00000599538VRK3chr1950504046-ENST00000221930TGFB1chr1941838186-231912776651318217
ENST00000443401VRK3chr1950504046-ENST00000221930TGFB1chr1941838186-159255088591167
ENST00000594092VRK3chr1950504046-ENST00000221930TGFB1chr1941838186-1779737125778217
ENST00000599538VRK3chr1950504047-ENST00000221930TGFB1chr1941838186-231912776651318217
ENST00000443401VRK3chr1950504047-ENST00000221930TGFB1chr1941838186-159255088591167
ENST00000594092VRK3chr1950504047-ENST00000221930TGFB1chr1941838186-1779737125778217
ENST00000599538VRK3chr1950504046-ENST00000221930TGFB1chr1941838185-231912776651318217
ENST00000443401VRK3chr1950504046-ENST00000221930TGFB1chr1941838185-159255088591167
ENST00000594092VRK3chr1950504046-ENST00000221930TGFB1chr1941838185-1779737125778217

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
ENST00000599538ENST00000221930VRK3chr1950504046-TGFB1chr1941838186-0.166067380.8339326
ENST00000443401ENST00000221930VRK3chr1950504046-TGFB1chr1941838186-0.22559210.7744079
ENST00000594092ENST00000221930VRK3chr1950504046-TGFB1chr1941838186-0.206010520.7939894
ENST00000599538ENST00000221930VRK3chr1950504047-TGFB1chr1941838186-0.166067380.8339326
ENST00000443401ENST00000221930VRK3chr1950504047-TGFB1chr1941838186-0.22559210.7744079
ENST00000594092ENST00000221930VRK3chr1950504047-TGFB1chr1941838186-0.206010520.7939894
ENST00000599538ENST00000221930VRK3chr1950504046-TGFB1chr1941838185-0.166067380.8339326
ENST00000443401ENST00000221930VRK3chr1950504046-TGFB1chr1941838185-0.22559210.7744079
ENST00000594092ENST00000221930VRK3chr1950504046-TGFB1chr1941838185-0.206010520.7939894

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 VRK3-TGFB1

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 VRK3-TGFB1 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 VRK3-TGFB1 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 VRK3-TGFB1

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 VRK3-TGFB1

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 VRK3-TGFB1

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 VRK3-TGFB1

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 VRK3-TGFB1

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 VRK3-TGFB1

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

Top

Related Diseases to VRK3-TGFB1

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