Hi ATpoint, thank you for your suggestion and feedback on it! What I want to achieve is just basically check expression levels for my genes of interest and as they have been annotated with the latest update of HGNC, I thought I need to do the same with the TCGA gene names. My concern was that I might not be able to find my genes of interest just because TCGA gene names are still using some "previous" names. Going with dummy names is an interesting idea. I will give it a try! Thank you!
How do I update the gene names of TCGA data correctly?
Dear all,
I want to work with RNA-seq TCGA data and as I am working with a list with genes of interests that is annotated based on the latest update of HGNC (04.06.2024), I wanted to do the same with the TCGA gene names. However, when I do this (using the ensembl gene ID), there are roughly 3,800 genes that I cannot match. I also tried to match the names but there are even more genes that do not match.
I am still a beginner in bioinformatics and I would be greatful for any tips or suggestions on how to annotate/up-date the TCGA gene names!
Thank you!
Best,
Ivana
• 2,373 views
•
link
updated
by
Zhenyu Zhang
•
written
by
Ivana •
•
written
by
Ivana •
1 answer
Updates in gene annotations might result in deprecation of certain gene names, while others might be added. If you use existing TCGA quantifications (without starting from fastq files) then I would really just focus on gene IDs that currently have a match in the recent annotation you want to use, and give the others some dummy name, like missing_[0-9]+. There is no point trying to force-match them somehow. If they're gone in recent HGNC then they're gone. Only true "good" workaround would be to process TCGA from fastq files on, but that is access-restricted and tedious. You cannot expect old annotations to perfectly match recent ones. If that was the case then the new annotations would be pointless, no?
Or just use the existing annotations in the TCGA databases, without updating. Is it really critical to "update" here?
Log in to answer this question.
More posts like this
-
Extracting mutation status from MAF files - TCGA
written by yahn •
-
Interpretation of mean-variance trend in voom
written by Ivana •
-
hg18, hg19, hg38
written by GG •
-
Get Gene families from Mus musculus
written by JulianC •
-
Identify or update refseq transcript version on snpeff
written by NB
-
How to avoid annotating multiple genes to a single CpG locus?
written by Mathias
-
How to annotate a chloroplast genome to find enzymes of interest
written by katherinefleury1 •
-
Matrix index error when using dnet package to calculate semantic similarity of genes
written by dolevrahat •
-
What are the GRCh37/38 tRNA gene names
written by AlexAbdulkaderKheirallah
-
Annotation Of Probeset-Level Data Of Affymetrix Hugene-1_0-St-V1 Chip
written by munch •
What do you mean by "cannot match"? Can you give us an example?
This is a classic bioinformatics question, and there are no standard way to do so. You are balancing your mappings between FPs and FNs.
I normally ensemble all the following mappings
And then you can setup a rule. My rule is if a mapping is not unique, I will manually inspect it.