Question: How can I get Ensembl CDS / 3'UTR and 5'UTR sequences
1
gravatar for 2405592M
11 months ago by
2405592M90
2405592M90 wrote:

Hi all,

I've finished doing an RNA-seq analysis and I've done some differential expression on my data set. I now want to do some codon usage analysis between conditions for my up and down regulated genes. How can one get the 3'UTR / CDS / 5'UTR sequences with ensembl IDs?

Thanks in advance!

ADD COMMENTlink modified 11 months ago by Emily_Ensembl21k • written 11 months ago by 2405592M90
1
gravatar for Emily_Ensembl
11 months ago by
Emily_Ensembl21k
EMBL-EBI
Emily_Ensembl21k wrote:

BioMart! Filter by your list of Ensembl IDs and get the sequences as Attributes. There's a help video to get you started if you've never used BioMart before.

ADD COMMENTlink written 11 months ago by Emily_Ensembl21k

Hi Emily,

I'm trying to do this through R using the following:

library(biomaRt)
listMarts()

ensembl <- useMart("ensembl")

datasets <- listDatasets(ensembl)
head(datasets)

ensembl = useDataset("mmusculus_gene_ensembl", mart = ensembl)

attributeNames_codon_usage <- c("ensembl_gene_id", "5utr", "coding", "3utr")

ourFilterType <- "ensembl_gene_id"

filterValues <- rownames(Day4_CONvsCRE)

codon_usage_mm10_annot_Day4_CONvsCRE <- getBM(attributes = attributeNames_codon_usage,
                                              filters = ourFilterType,
                                              values = filterValues,
                                              mart = ensembl)

but I keep getting the following error:

Error in scan(file = file, what = what, sep = sep, quote = quote, dec = dec, : line 1 did not have 4 elements

ADD REPLYlink written 11 months ago by 2405592M90

You can't select multiple sequence types in one query, you would need to run them as three queries.

I would also recommend getting the transcript ID in your attributes, as you'll get a sequence for each one.

ADD REPLYlink modified 11 months ago • written 11 months ago by Emily_Ensembl21k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2105 users visited in the last hour