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:AURKA-SOGA1

Fusion Gene and Fusion Protein Summary

check button Fusion gene summary
Fusion partner gene informationFusion gene name: AURKA-SOGA1
FusionPDB ID: 8490
FusionGDB2.0 ID: 8490
HgeneTgene
Gene symbol

AURKA

SOGA1

Gene ID

6790

140710

Gene nameaurora kinase Asuppressor of glucose, autophagy associated 1
SynonymsAIK|ARK1|AURA|BTAK|PPP1R47|STK15|STK6|STK7C20orf117|KIAA0889|SOGA
Cytomap

20q13.2

20q11.23

Type of geneprotein-codingprotein-coding
Descriptionaurora kinase Aaurora 2aurora/IPL1-like kinaseaurora/IPL1-related kinase 1breast tumor-amplified kinaseprotein phosphatase 1, regulatory subunit 47serine/threonine protein kinase 15serine/threonine-protein kinase 6serine/threonine-protein kinase aprotein SOGA1SOGA family member 1suppressor of glucose by autophagysuppressor of glucose from autophagysuppressor of glucose, autophagy-associated protein 1
Modification date2020032920200313
UniProtAcc

Q9NWT8

Main function of 5'-partner protein: FUNCTION: May act as a negative regulator of Aurora-A kinase, by down-regulation through proteasome-dependent degradation.

O94964

Main function of 5'-partner protein: FUNCTION: Regulates autophagy by playing a role in the reduction of glucose production in an adiponectin- and insulin-dependent manner. {ECO:0000250}.
Ensembl transtripts involved in fusion geneENST idsENST00000312783, ENST00000347343, 
ENST00000371356, ENST00000395907, 
ENST00000395909, ENST00000395911, 
ENST00000395913, ENST00000395914, 
ENST00000395915, 
ENST00000357779, 
ENST00000456801, ENST00000237536, 
ENST00000279034, 
Fusion gene scores for assessment (based on all fusion genes of FusionGDB 2.0)* DoF score6 X 8 X 3=14410 X 10 X 4=400
# samples 812
** MAII scorelog2(8/144*10)=-0.84799690655495
possibly effective Gene in Pan-Cancer Fusion Genes (peGinPCFGs).
DoF>8 and MAII<0
log2(12/400*10)=-1.73696559416621
possibly effective Gene in Pan-Cancer Fusion Genes (peGinPCFGs).
DoF>8 and MAII<0
Fusion gene context

PubMed: AURKA [Title/Abstract] AND SOGA1 [Title/Abstract] AND fusion [Title/Abstract]

Fusion neoantigen context

PubMed: AURKA [Title/Abstract] AND SOGA1 [Title/Abstract] AND neoantigen [Title/Abstract]

Most frequent breakpoint (based on all fusion genes of FusionGDB 2.0)AURKA(54956489)-SOGA1(35445872), # samples:2
Anticipated loss of major functional domain due to fusion event.AURKA-SOGA1 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.
AURKA-SOGA1 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.
AURKA-SOGA1 seems lost the major protein functional domain in Hgene partner, which is a epigenetic factor due to the frame-shifted ORF.
AURKA-SOGA1 seems lost the major protein functional domain in Hgene partner, which is a essential gene due to the frame-shifted ORF.
AURKA-SOGA1 seems lost the major protein functional domain in Hgene partner, which is a IUPHAR drug target due to the frame-shifted ORF.
AURKA-SOGA1 seems lost the major protein functional domain in Hgene partner, which is a kinase due to the frame-shifted ORF.
AURKA-SOGA1 seems lost the major protein functional domain in Tgene partner, which is a essential gene 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
HgeneAURKA

GO:0006468

protein phosphorylation

21600873|21820309

HgeneAURKA

GO:0009611

response to wounding

19435814

HgeneAURKA

GO:0032091

negative regulation of protein binding

21820309

HgeneAURKA

GO:0097421

liver regeneration

19435814



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/chr20:54956489/chr20:35445872)
- 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 AURKA (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 SOGA1 (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)
ENST00000371356AURKAchr2054956489-ENST00000237536SOGA1chr2035445872-137958371847511577
ENST00000371356AURKAchr2054956489-ENST00000279034SOGA1chr2035445872-38568371835301170
ENST00000312783AURKAchr2054956489-ENST00000237536SOGA1chr2035445872-1390694824348621539
ENST00000312783AURKAchr2054956489-ENST00000279034SOGA1chr2035445872-396794824336411132
ENST00000395913AURKAchr2054956489-ENST00000237536SOGA1chr2035445872-1385289415948081549
ENST00000395913AURKAchr2054956489-ENST00000279034SOGA1chr2035445872-391389415935871142
ENST00000371356AURKAchr2054956488-ENST00000237536SOGA1chr2035445872-137958371847511577
ENST00000371356AURKAchr2054956488-ENST00000279034SOGA1chr2035445872-38568371835301170
ENST00000312783AURKAchr2054956488-ENST00000237536SOGA1chr2035445872-1390694824348621539
ENST00000312783AURKAchr2054956488-ENST00000279034SOGA1chr2035445872-396794824336411132
ENST00000395913AURKAchr2054956488-ENST00000237536SOGA1chr2035445872-1385289415948081549
ENST00000395913AURKAchr2054956488-ENST00000279034SOGA1chr2035445872-391389415935871142

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
ENST00000371356ENST00000237536AURKAchr2054956489-SOGA1chr2035445872-0.0009886440.9990113
ENST00000371356ENST00000279034AURKAchr2054956489-SOGA1chr2035445872-0.0083886760.9916113
ENST00000312783ENST00000237536AURKAchr2054956489-SOGA1chr2035445872-0.001144840.9988551
ENST00000312783ENST00000279034AURKAchr2054956489-SOGA1chr2035445872-0.0098925890.9901074
ENST00000395913ENST00000237536AURKAchr2054956489-SOGA1chr2035445872-0.0010568320.99894315
ENST00000395913ENST00000279034AURKAchr2054956489-SOGA1chr2035445872-0.0090108210.99098915
ENST00000371356ENST00000237536AURKAchr2054956488-SOGA1chr2035445872-0.0009886440.9990113
ENST00000371356ENST00000279034AURKAchr2054956488-SOGA1chr2035445872-0.0083886760.9916113
ENST00000312783ENST00000237536AURKAchr2054956488-SOGA1chr2035445872-0.001144840.9988551
ENST00000312783ENST00000279034AURKAchr2054956488-SOGA1chr2035445872-0.0098925890.9901074
ENST00000395913ENST00000237536AURKAchr2054956488-SOGA1chr2035445872-0.0010568320.99894315
ENST00000395913ENST00000279034AURKAchr2054956488-SOGA1chr2035445872-0.0090108210.99098915

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 AURKA-SOGA1

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

Top

Potential FusionNeoAntigen Information of AURKA-SOGA1 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.
AURKA-SOGA1_54956488_35445872.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)
AURKA-SOGA1chr2054956488chr2035445872948HLA-B18:01DEQRTATH0.99850.8636614
AURKA-SOGA1chr2054956488chr2035445872948HLA-B14:01QRTATHSL0.99630.7688816
AURKA-SOGA1chr2054956488chr2035445872948HLA-B14:02QRTATHSL0.99630.7688816
AURKA-SOGA1chr2054956488chr2035445872948HLA-B15:10QRTATHSL0.99490.5936816
AURKA-SOGA1chr2054956488chr2035445872948HLA-B15:37QRTATHSL0.96480.6423816
AURKA-SOGA1chr2054956488chr2035445872948HLA-A30:08RTATHSLKK0.99870.7313918
AURKA-SOGA1chr2054956488chr2035445872948HLA-B14:01EQRTATHSL0.99720.7804716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B14:02EQRTATHSL0.99720.7804716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B39:01EQRTATHSL0.99220.8767716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B08:01EQRTATHSL0.98020.5556716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B39:13EQRTATHSL0.97270.8815716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B08:09EQRTATHSL0.95490.6158716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B15:18EQRTATHSL0.81030.5781716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B15:03EQRTATHSL0.78320.6171716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B15:02EQRTATHSL0.64650.8439716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B13:01EQRTATHSL0.50570.7939716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B18:01DEQRTATHSL0.89320.9113616
AURKA-SOGA1chr2054956488chr2035445872948HLA-B39:12QRTATHSL0.99810.816816
AURKA-SOGA1chr2054956488chr2035445872948HLA-B14:03QRTATHSL0.96210.762816
AURKA-SOGA1chr2054956488chr2035445872948HLA-B27:14QRTATHSLK0.99750.7633817
AURKA-SOGA1chr2054956488chr2035445872948HLA-B39:09EQRTATHSL0.99470.6775716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B39:12EQRTATHSL0.9890.8744716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B42:02EQRTATHSL0.98650.5521716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B39:05EQRTATHSL0.97560.8612716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B39:08EQRTATHSL0.91160.8716716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B14:03EQRTATHSL0.9070.8649716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B15:04EQRTATHSL0.8620.7526716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B15:07EQRTATHSL0.85370.5384716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B15:21EQRTATHSL0.64590.8099716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B51:07EQRTATHSL0.15670.5363716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B18:06DEQRTATH0.99870.8726614
AURKA-SOGA1chr2054956488chr2035445872948HLA-B18:05DEQRTATH0.99850.8636614
AURKA-SOGA1chr2054956488chr2035445872948HLA-B39:31QRTATHSL0.99840.8154816
AURKA-SOGA1chr2054956488chr2035445872948HLA-B18:03DEQRTATH0.99740.8559614
AURKA-SOGA1chr2054956488chr2035445872948HLA-A30:01RTATHSLKK0.99880.8772918
AURKA-SOGA1chr2054956488chr2035445872948HLA-B27:10QRTATHSLK0.99840.8621817
AURKA-SOGA1chr2054956488chr2035445872948HLA-B39:02EQRTATHSL0.9930.8801716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B39:31EQRTATHSL0.99240.8741716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B08:18EQRTATHSL0.98020.5556716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B08:12EQRTATHSL0.96920.6987716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B15:27EQRTATHSL0.9690.7591716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B39:11EQRTATHSL0.9280.8415716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B15:73EQRTATHSL0.84760.7471716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B15:09EQRTATHSL0.81090.5358716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B15:54EQRTATHSL0.81010.6697716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B15:53EQRTATHSL0.80640.6937716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B15:30EQRTATHSL0.78310.7284716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B48:02EQRTATHSL0.77220.8114716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B15:50EQRTATHSL0.73980.8459716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B15:12EQRTATHSL0.49590.7661716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B35:20EQRTATHSL0.46690.8592716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B18:04EQRTATHSL0.40430.8999716
AURKA-SOGA1chr2054956488chr2035445872948HLA-B18:08DEQRTATHSL0.8950.9047616
AURKA-SOGA1chr2054956488chr2035445872948HLA-B18:05DEQRTATHSL0.89320.9113616
AURKA-SOGA1chr2054956488chr2035445872948HLA-B18:06DEQRTATHSL0.88170.919616
AURKA-SOGA1chr2054956488chr2035445872948HLA-B18:03DEQRTATHSL0.87370.9048616
AURKA-SOGA1chr2054956488chr2035445872948HLA-B18:11DEQRTATHSL0.69740.8628616

Top

Potential FusionNeoAntigen Information of AURKA-SOGA1 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 AURKA-SOGA1

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
1077DEQRTATHSLKKRGAURKASOGA1chr2054956488chr2035445872948

Top

Filtering FusionNeoAntigens Through Checking the Interaction with HLAs in 3D of AURKA-SOGA1

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:023BVN1077DEQRTATHSLKKRG-7.15543-7.26883
HLA-B14:023BVN1077DEQRTATHSLKKRG-4.77435-5.80965
HLA-B52:013W391077DEQRTATHSLKKRG-6.80875-6.92215
HLA-B52:013W391077DEQRTATHSLKKRG-4.20386-5.23916
HLA-A11:014UQ21077DEQRTATHSLKKRG-7.5194-8.5547
HLA-A11:014UQ21077DEQRTATHSLKKRG-6.9601-7.0735
HLA-A24:025HGA1077DEQRTATHSLKKRG-7.52403-7.63743
HLA-A24:025HGA1077DEQRTATHSLKKRG-5.82433-6.85963
HLA-B27:056PYJ1077DEQRTATHSLKKRG-3.28285-4.31815
HLA-B44:053DX81077DEQRTATHSLKKRG-5.91172-6.94702
HLA-B44:053DX81077DEQRTATHSLKKRG-4.24346-4.35686

Top

Vaccine Design for the FusionNeoAntigens of AURKA-SOGA1

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
AURKA-SOGA1chr2054956488chr2035445872614DEQRTATHGATGAGCAGAGAACTGCTACTCAT
AURKA-SOGA1chr2054956488chr2035445872616DEQRTATHSLGATGAGCAGAGAACTGCTACTCATTCCTTG
AURKA-SOGA1chr2054956488chr2035445872716EQRTATHSLGAGCAGAGAACTGCTACTCATTCCTTG
AURKA-SOGA1chr2054956488chr2035445872816QRTATHSLCAGAGAACTGCTACTCATTCCTTG
AURKA-SOGA1chr2054956488chr2035445872817QRTATHSLKCAGAGAACTGCTACTCATTCCTTGAAG
AURKA-SOGA1chr2054956488chr2035445872918RTATHSLKKAGAACTGCTACTCATTCCTTGAAGAAA

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 AURKA-SOGA1

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
OVAURKA-SOGA1chr2054956488ENST00000312783chr2035445872ENST00000237536TCGA-61-1728

Top

Potential target of CAR-T therapy development for AURKA-SOGA1

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 AURKA-SOGA1

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

Top

Related Diseases to AURKA-SOGA1

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