This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Retrieving gene details

Hi, I have a list of genes and want to create a .tsv file containing gene name, respective chromosome, and positions of the left and right ends of the gene. I saw this post [link]on researchgate for the same task but I do not find any options to select datasets as mentioned in there, can anyone suggest any other alternative ways of achieving this goal.

The list looks like

TP53

BRCA1

The expected output is supposed to look like;

Gene chr start end

BRCA1 chr17 43044295 43170245

TP53 chr17 7661779 7687550

genes biomart
• 2,413 views • link
written by Mk • Mk

any options to select datasets as mentioned in there

What do you mean by that? Are you not finding the genome you are interested in? You may want to watch this tutorial on BioMart to see if that answers your question.

• 0 views • link
written by GenoMax GenoMax

Thank you genomax2, I was able to follow the tutorial you provided and find the required solutions.

• 0 views • link
written by Mk • Mk

2 answers

You could use UCSC Goldenpath; e.g. for hg19:

$ wget -qO- http://hgdownload.cse.ucsc.edu/goldenpath/hg19/database/knownGene.txt.gz | gunzip -c > knownGene.txt

You could use the kgXref table to map naming schemes via grep etc.:

$ wget -qO- http://hgdownload.cse.ucsc.edu/goldenpath/hg19/database/kgXref.txt.gz | gunzip -c > kgXref.txt

Or you could do the mapping via the web interface in the UCSC Genome Browser.

• 0 views • link
written by Alex Reynolds Alex Reynolds

You can use biomaRt to get information easily.

mart <- useMart("ensembl")
datasets <- listDatasets(mart)
mart <- useDataset("hsapiens_gene_ensembl",mart)

Use getBM() to filter characteristic you want.

filter_data <- getBM(
  attributes = c("external_gene_name", "chromosome_name", "transcript_start", "transcript_end"),
  mart = mart
)

And the result you get:

> head(filter_data)
  external_gene_name  chromosome_name transcript_start transcript_end
1          RNU6-280P CHR_HG2128_PATCH         67546651       67546754
2              Y_RNA CHR_HG2128_PATCH         67631019       67631127
3       RP11-222G7.2 CHR_HG2191_PATCH         74823667       74824187
4    Clostridiales-1 CHR_HG2022_PATCH         91356877       91357036
5      RP11-654C22.2  CHR_HG126_PATCH         72505075       72550889
6      RP11-315F22.1 CHR_HG2233_PATCH        239734956      239735538
• 0 views • link
written by zjhzwang zjhzwang

I do not know R programming, but thank you zjhzwang for the suggestion.

• 0 views • link
written by Mk • Mk

Log in to answer this question.

More posts like this