Genes with multiple chromosome locations

Hi, I am trying to annotate my list of gene IDs, some of which have multiple loci, e.g.:

   refseq_mrna hgnc_symbol   gene_biotype chromosome_name start_position end_position
7    NM_000076      CDKN1C protein_coding  HSCHR11_1_CTG7         115392       118091
8    NM_000076      CDKN1C protein_coding              11        2883213      2885775

What is the recommended practice for collapsing gene annotation with multiple entries of genes with alternative loci for a given gene?

gene annotation loci alternative

• 1,414 views • link

updated 2.4 years ago by Alex Reynolds

• written 2.4 years ago by Pac314 •

Where did this annotation originate from?

CDKN1C seems to be annotated only at one location. https://www.ncbi.nlm.nih.gov/gene/1028/ and https://www.genenames.org/data/gene-symbol-report/#!/hgnc_id/HGNC:1786

• 0 views • link

written 2.4 years ago by GenoMax

Thanks for your reply. I obtained this annotation using the biomaRt R library:

library(biomaRt)

ensembl <- useMart("ensembl", dataset = "hsapiens_gene_ensembl")

attrs <- c('refseq_mrna', 'hgnc_symbol', 'gene_biotype', 'chromosome_name', 
           'start_position', 'end_position')

annot_mrna <- getBM(attributes = attrs,
                     filters = 'refseq_mrna',
                     values = rownames(cts_mat),
                     mart = ensembl)

I have multiple gene IDs with different chromosome names like the above example.

• 0 views • link

written 2.4 years ago by Pac314 •

You might look for a canonical gene isoform. UCSC may have some useful advice: https://genome.ucsc.edu/FAQ/FAQgenes.html#singledownload

• 0 views • link

written 2.4 years ago by Alex Reynolds

0 answers

No answers yet.

More posts like this

Tutorial: Single-cell RNA-seq: Annotation: Celltype auto annotation with MetaTiME
written 2.6 years ago by Julia Ma •

MetaTiME learns data-driven, interpretable, and reproducible gene programs by integrating millions of single cells from hundreds of tumor scRNA-seq data. The idea is to learn …
Get SNPs with MAF 0.05 for all genes - attempt with BiomaRt
written 3.7 years ago by camerond

I have a working script which takes a list of genes and maps all the SNPs with MAF >= 0.05 to each gene using BiomaRt. …
Tutorial: Unsupervised clustering on gene expression data
written 4.8 years ago by Hamid Ghaedi

Clustering is a data mining method to identify unknown possible groups of items solely based on intrinsic features and no external variables. Basically, clustering includes …
gene ensemble ID- no symbol
written 5.4 years ago by Rob •

Hello friends, I used this code to convert ensemble gene ID to gene symbol: df <- read.csv("gene_Ens_ID.csv") library('biomaRt') mart <- useDataset("hsapiens_gene_ensembl", useMart("ensembl")) attributeNames <- c('ensembl_gene_id','external_gene_name','description','gene_biotype', …
Annotating DEXseq outputs
written 6.7 years ago by gayatri281 •

Hi All, I am trying to annotate to get gene symbols,from the outputs of DEXseq results that has Ensembl gene ids in the following pattern …
Duplicate and Different Annotations using biomaRt
written 6.8 years ago by JourneyToAbyss

For the affy HG-U133A_2, I used biomaRt to retrieve the annotation. That being said, I am finding multiple entities for the same probe ID. But …
How to avoid annotating multiple genes to a single CpG locus?
written 8.2 years ago by Mathias

I am downloading TCGA DNA methylation data using TCGAbiolinks. The rowData() function output looks a bit messy; all ensembl ids with which a cpg overlaps …
different positions for a gene
written 10.0 years ago by hamarillo •

Hello, I really don't understand why the same gene has a different number of transcripts, but moreover, a different start position. In this case for …
Retrieve All Genes Contained Within A Specific Chromosomal Region Using R And Biomart
written 13.4 years ago by secretjess

EDIT: Now fixed, using filters = c("chromosome_name","start","end","biotype") instead. Question still stands though if anyone knows? Original question: I've written some code which I think should …