gene coordinates from TxDb.Hsapiens.UCSC.hg19.knownGene do not match NCBI gene bank

I retrieved gene cordinates from TxDb.Hsapiens.UCSC.hg19.knownGene using:

gene_info = genes(TxDb.Hsapiens.UCSC.hg19.knownGene)

But I found some of the genes, for example, PTPN20 (id: 26095) is hugely different from NCBI gene bank. It is error information or I am doing something wrong?

subset(gene_info, gene_id=='26095')

GRanges object with 1 range and 1 metadata column:
        seqnames            ranges strand |     gene_id
           <Rle>         <IRanges>  <Rle> | <character>
  26095    chr10 46550123-48827924      - |       26095

The gene coordinates I got from gene bank: https://www.ncbi.nlm.nih.gov/gene/?term=26095

gene

• 2,496 views • link

updated 6.9 years ago by GenoMax

• written 6.9 years ago by yliueagle

It looks like this gene was annotated in two places in hg19 build (http://uswest.ensembl.org/Homo_sapiens/Gene/Summary?db=core;g=ENSG00000204179;r=10:46911396-47002488 look in other assemblies). This seems to have been resolved in hg38, where only one copy is annotated.

https://www.ncbi.nlm.nih.gov/gene/26095
https://www.genenames.org/data/gene-symbol-report/#!/hgnc_id/HGNC:23423

• 0 views • link

written 6.9 years ago by GenoMax

Adding on this, see hg19:

enter image description here

• 0 views • link

written 6.9 years ago by ATpoint

0 answers

No answers yet.

More posts like this

Adding gene to GRange objects with single base
written 2.4 years ago by QX •

Hi all, I have a Grange objects look like this: ``` GRanges object with 4598553 ranges and 2 metadata columns: seqnames ranges strand | ref …
transcriptome in Granges format
written 4.7 years ago by alexmondaini •

Hello everyone, I would like to have some opinions on whether my approach to get a transcriptome is correct, perhaps this is naive question. In …
Window Length to Genomic Coordinates
written 5.2 years ago by shivangi.agarwal800

Hello Guys, I am trying to convert window length into genomic coordinates using below mentioned commands, following a post **https://www.biostars.org/p/169171/** library(Homo.sapiens) human.genes = genes(TxDb.Hsapiens.UCSC.hg19.knownGene) human.genes …
[R] subsetByOverlaps (Granges list and Granges object) not giving specific hits
written 6.1 years ago by Alewa

Dear Biostars, I have GRangesList object of length 15 with metadata (of gene annotations & copy number), created small GRanges object to use to query …
GRanges: merge metadata from the same ranges
written 6.7 years ago by arseru •

Hi everyone, I'm having trouble trying to merge the metadata from 3 different GRanges object that fall in the same range. I'll give an example …
locateVariants missing some mutation site
written 8.8 years ago by noeD

Hello! I am using Variant Annotation package from Bioconductor to annotate my variant but it missed some variants in the output. For example, in `var` …
How to save GRanges object to csv file
written 8.9 years ago by akij

Currently this is what I do library(Homo.sapiens) l <- genes(TxDb.Hsapiens.UCSC.hg19.knownGene) > l GRanges object with 23056 ranges and 1 metadata column: seqnames ranges strand | …
GRanges intersect operation
written 10.1 years ago by bsmith030465

I am trying to find the intersection of two GRanges objects. The first object has the coordinates for genes (grgenes below), the second mimics my …
How to annotate one GenomicRanges object using another?
written 10.2 years ago by endrebak

I have two genomicranges, one containing the gencode GFF v24 library(GenomicFeatures) gen = makeTxDbFromGFF("gencode.gff") tss = promoters(gen, upstream = 5000, downstream = 5000) # > …
How to find all miRNA locations in on GRch37 ?
written 11.8 years ago by Aurelie MLB

Hello, I am trying to get the location of miRNAs on the human genome GRch37 from UCSC. I tried to use bioconductor to do this. …