Comments (9)
Hello,
This is because some of the datasets were processed using Ensemble v87. We updated the version later on.
Thanks,
Naim
from grein.
Hi team,
Thanks for your reply. Unfortunately, I also checked ensembl v87 and didn't find these genes named "COX1","ND5","COX3","CYTB","COX2" and more, which show up in GSE58375_GeneLevel_Raw_data.csv. I was thinking that if Ensembl v91 gtf is applied to processed data recently so I checked newly processed data (a month ago) such as GSE97615 and GSE115348 and the situation is the same. Could you please check it again? Thank you.
Best,
Zheng
from grein.
Hi Zheng,
This is actually v86 (sorry for the typo!). Our back-end Bioconductor version of Ensembl DB was 86. That's why the discrepancy. We are in the process of upgrading to v91.
Best,
Naim
from grein.
Hi team,
Thanks for your reply. unfortunately I checked ensembl v86 and the situation remains. Could you please check your gtf version again? Thanks.
Best,
Zheng
from grein.
Hi Zheng,
So these are alternative gene symbols available in the AnnotationDbi package. Please check with the corresponding Ensemble id of the genes. I just checked and they are all there. Please let me know if you cannot find them.
from grein.
Hi team,
Thanks for your reply. I actually downloaded ensembl v86(along with v91 and v87) human gif from ensembl and found these genes don't belong to any part of ensembl v86/v91/v87 gtf file. And why did you use AnnotationDbi packages? To convert ENSG to gene symbol? If so, why don't you just extract ENSG-gene_symbol pairs and convert ENSG to gene symbol based on corresponding gtf file? Even though, I still can't find any trace in ensembl v91/v87/v86 gtf by searching the corresponding ENSG ID. It's disturbing that we cannot know which gtf file you actually use, whether the gene names are correct(It has huge effects on downstream analysis), and thus, the credibility of this pipeline should be concerned.
Best,
Zheng
from grein.
As I said, we are using Ensembl version 86. I used "Homo_sapiens.GRCh38.86.gtf" file and double-checked with the ensemble ids. For example, search with ENSG00000198804 for COX1 (alternative symbol MT-CO1). Even if you simply google with the ensembl id, you will see the symbols.
Here, a snapshot of v86 for ensg id ENSG00000198804:
seqnames start end width strand source type score phase gene_id gene_version gene_name gene_source gene_biotype
2574686 MT 5904 7445 1542 + insdc gene NA NA ENSG00000198804 2 MT-CO1 insdc protein_coding
2574687 MT 5904 7445 1542 + ensembl transcript NA NA ENSG00000198804 2 MT-CO1 insdc protein_coding
2574688 MT 5904 7445 1542 + ensembl exon NA NA ENSG00000198804 2 MT-CO1 insdc protein_coding
Since, it is an automated pipeline, AnnotationDbi is the easiest/quickest way for annotation conversion rather than reading the huge gtf file every time.
from grein.
from grein.
As I said, we are using Ensembl version 86. I used "Homo_sapiens.GRCh38.86.gtf" file and double-checked with the ensemble ids. For example, search with ENSG00000198804 for COX1 (alternative symbol MT-CO1). Even if you simply google with the ensembl id, you will see the symbols.
Here, a snapshot of v86 for ensg id ENSG00000198804:seqnames start end width strand source type score phase gene_id gene_version gene_name gene_source gene_biotype
2574686 MT 5904 7445 1542 + insdc gene NA NA ENSG00000198804 2 MT-CO1 insdc protein_coding
2574687 MT 5904 7445 1542 + ensembl transcript NA NA ENSG00000198804 2 MT-CO1 insdc protein_coding
2574688 MT 5904 7445 1542 + ensembl exon NA NA ENSG00000198804 2 MT-CO1 insdc protein_codingSince, it is an automated pipeline, AnnotationDbi is the easiest/quickest way for annotation conversion rather than reading the huge gtf file every time.
Thank you for your reply. your answer really helps me. Last question. When do you consider mapping old version alignment RNA-Seq to an uniform genome version? And can we know which dataset is mapped to ensembl v86 or ensembl v91?
Zheng
from grein.
Related Issues (20)
- website down? HOT 3
- Analysis not done HOT 2
- All counts are zero HOT 1
- Gene number different from different dataset HOT 2
- Combine multiple dataset from GREIN HOT 2
- Duplicated gene symbol HOT 2
- GEO_data_processing.R get error
- GEO_data_processing error
- Site not loading HOT 2
- SIte not loading (bis) HOT 1
- new datasets HOT 3
- ALL SAMPLES ARE NOT INCLUDED FOR ANALYSIS HOT 2
- Count table at gene level HOT 7
- Hi. I downloaded raw counts for GSE83577, but the data corresponds to normalized counts. Is this a problem just for this dataset of a platform bug? thanks
- Any possibility for mass download or API for dataset HOT 1
- Incomplete Dataset and New Datasets HOT 4
- Some datasets are not loading completely
- Fast download of Count Matrix? HOT 1
- API for download and TPM calculation HOT 1
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from grein.