Otdfbt

What is wrong with my code? RGB potentiometer

How to make a language evolve quickly?

Is there a need for better software for writers?

Why do Thanos' punches not kill Captain America or at least cause vital wounds?

No such column 'DeveloperName' on entity 'RecordType' after Summer '19 release on sandbox

is it permitted to swallow spit on a fast day?

How to efficiently lower your karma

Extending Kan fibrations, without using minimal fibrations

Improving Sati-Sampajañña (situative wisdom)

Renting a house to a graduate student in my department

Does the 500 feet falling cap apply per fall, or per turn?

Best species to breed to intelligence

Why was the ancient one so hesitant to teach Dr Strange the art of sorcery

Series that evaluates to different values upon changing order of summation

Has magnetic core memory been used beyond the Moon?

Is it a Munchausen Number?

What does formal training in a field mean?

Company threw a surprise party for the CEO, 3 weeks later management says we have to pay for it, do I have to?

How are one-time password generators like Google Authenticator different from having two passwords?

Why is the Sun made of light elements only?

Can a surprised creature fall prone voluntarily on their turn?

Why are low spin tetrahedral complexes so rare?

Can 'sudo apt-get remove [write]' destroy my Ubuntu?

Why is it wrong to *implement* myself a known, published, widely believed to be secure crypto algorithm?

CDS length for each human gene

Is there a way to retrieve ENSEMBL IDs from a search query?Converting Ensembl Gene IDs to Entrez Gene IDs through biomartIdentifying relevant SNPs from a listCounting the number of paralogues for mouse genes gives me the wrong frequency in RA good tool for gene locus visualizationDownload proteomes from NCBI based only on binomial namesGet Gene Expression Matrix from GEOqueryFinding gene length using ensembl IDFinding gene name from human genome using SP1 transcrition factor binding site from Postion Weight MatrixRetrieve RNA sequencing data for human p53 colon cancer cell lines

1

$begingroup$

Does anyone know where and how could I download a list of all human genes and the length of the coding sequence for each gene? Is it possible to do this on the NCBI site, ensembl?

gene sequence-analysis ncbi ensembl

share|improve this question

edited May 1 at 7:51

Kamil S Jaron

3,082942

asked Apr 30 at 14:52

solimanelefantsolimanelefant

1083

$endgroup$

• 
  
  
  

  
  2
  

  
  
  
  
  
  

  
  $begingroup$
  Which coding sequence? I mean, do you just want whichever has been designated the 'canonical' transcript or do you want all possible isoforms?
  $endgroup$
  – terdon♦
  Apr 30 at 14:54
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  Hi terdon, thanks for the quick reply! Yes exactly, the canonical transcript is good enough!
  $endgroup$
  – solimanelefant
  Apr 30 at 14:56
  
  

  
  
  
  


• 
  
  
  

  
  1
  

  
  
  
  
  
  

  
  $begingroup$
  Michael G. suggests to take a look at relevant front-end, NCBI's eFetch. Which is supposedly perfect for what you need.
  $endgroup$
  – Kamil S Jaron
  May 1 at 7:53
  

  
  
  
  

add a comment | 

1

$begingroup$

Does anyone know where and how could I download a list of all human genes and the length of the coding sequence for each gene? Is it possible to do this on the NCBI site, ensembl?

gene sequence-analysis ncbi ensembl

share|improve this question

edited May 1 at 7:51

Kamil S Jaron

3,082942

asked Apr 30 at 14:52

solimanelefantsolimanelefant

1083

$endgroup$

• 
  
  
  

  
  2
  

  
  
  
  
  
  

  
  $begingroup$
  Which coding sequence? I mean, do you just want whichever has been designated the 'canonical' transcript or do you want all possible isoforms?
  $endgroup$
  – terdon♦
  Apr 30 at 14:54
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  Hi terdon, thanks for the quick reply! Yes exactly, the canonical transcript is good enough!
  $endgroup$
  – solimanelefant
  Apr 30 at 14:56
  
  

  
  
  
  


• 
  
  
  

  
  1
  

  
  
  
  
  
  

  
  $begingroup$
  Michael G. suggests to take a look at relevant front-end, NCBI's eFetch. Which is supposedly perfect for what you need.
  $endgroup$
  – Kamil S Jaron
  May 1 at 7:53
  

  
  
  
  

add a comment | 

$begingroup$

Does anyone know where and how could I download a list of all human genes and the length of the coding sequence for each gene? Is it possible to do this on the NCBI site, ensembl?

gene sequence-analysis ncbi ensembl

share|improve this question

edited May 1 at 7:51

Kamil S Jaron

3,082942

asked Apr 30 at 14:52

solimanelefantsolimanelefant

1083

$endgroup$

share|improve this question

edited May 1 at 7:51

Kamil S Jaron

3,082942

asked Apr 30 at 14:52

solimanelefantsolimanelefant

1083

share|improve this question

edited May 1 at 7:51

Kamil S Jaron

3,082942

asked Apr 30 at 14:52

solimanelefantsolimanelefant

1083

edited May 1 at 7:51

Kamil S Jaron

3,082942

edited May 1 at 7:51

Kamil S Jaron

3,082942

asked Apr 30 at 14:52

solimanelefantsolimanelefant

1083

asked Apr 30 at 14:52

solimanelefantsolimanelefant

1083

2 Answers
2

active

oldest

votes

4

$begingroup$

While I haven't found a way to limit the results to the canonical transcript only, you can get a list of genes, transcripts and their CDS lengths using Ensemble's BioMart. I have already set it up for you, you can see the results, and modify them, here (click on the "Results" link if you don't see them).

Essentially, you just need to go to BioMart, and

1. select "Ensembl Genes 96" (the number will change if the version changes) as the database and "uman Genes" as the dataset.




2. Click on "Filters", and set Gene type to coding and Transcript type to protein_coding.




3. From "Attributes", select whatever you want to see. The "CDS Length" is under "Structures".

share|improve this answer

edited Apr 30 at 22:23

answered Apr 30 at 15:48

terdon♦terdon

4,9902830

$endgroup$

add a comment | 

0

$begingroup$

Ensembl has an FTP site that allows you to select and download only the coding sequences from many different genomes.
https://useast.ensembl.org/info/data/ftp/index.html

To determine the length of those sequences, download the associated gtf or gff3 annotation file. The annotation file is tab delim. The fourth and fifth column represent genomic loci of the annotated region. Subtract the amount in the fourth column from the amount in the fifth column to yield the length of all the annotated features.
You can easily do this in R after loading the file in the environment using the Magrittr library. The following code will create a new column called gene_length with the associated gene lengths.

install.packages("magrittr")
# this only needs to be done once
library(magrittr)
# must be run each time the library is neaded
annotation.gtf <- read.table("path/to/annotation.gtf")
annotation.gtf$start <- annotation.gtf[,4]
annotation.gtf$end <- annotation.gtf[,5]
annotation.new-column.gtf <- annotation.gtf %>%
 mutate(gene_length=end-start)

share|improve this answer

edited May 1 at 13:03

Kamil S Jaron

3,082942

answered May 1 at 11:33

Drew J-HDrew J-H

12

$endgroup$

add a comment | 

Your Answer

StackExchange.ready(function()
var channelOptions =
tags: "".split(" "),
id: "676"
;
initTagRenderer("".split(" "), "".split(" "), channelOptions);

StackExchange.using("externalEditor", function()
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled)
StackExchange.using("snippets", function()
createEditor();
);

else
createEditor();

);

function createEditor()
StackExchange.prepareEditor(
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: false,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: null,
bindNavPrevention: true,
postfix: "",
imageUploader:
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
,
onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
);



);

draft saved

draft discarded

Sign up or log in

StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);

Sign up using Google

Sign up using Facebook

Sign up using Email and Password

Post as a guest

Name

Email

Required, but never shown

StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fbioinformatics.stackexchange.com%2fquestions%2f8552%2fcds-length-for-each-human-gene%23new-answer', 'question_page');

);

Post as a guest

Name

Email

Required, but never shown

4

$begingroup$

While I haven't found a way to limit the results to the canonical transcript only, you can get a list of genes, transcripts and their CDS lengths using Ensemble's BioMart. I have already set it up for you, you can see the results, and modify them, here (click on the "Results" link if you don't see them).

Essentially, you just need to go to BioMart, and

1. select "Ensembl Genes 96" (the number will change if the version changes) as the database and "uman Genes" as the dataset.




2. Click on "Filters", and set Gene type to coding and Transcript type to protein_coding.




3. From "Attributes", select whatever you want to see. The "CDS Length" is under "Structures".

share|improve this answer

edited Apr 30 at 22:23

answered Apr 30 at 15:48

terdon♦terdon

4,9902830

$endgroup$

add a comment | 

4

$begingroup$

While I haven't found a way to limit the results to the canonical transcript only, you can get a list of genes, transcripts and their CDS lengths using Ensemble's BioMart. I have already set it up for you, you can see the results, and modify them, here (click on the "Results" link if you don't see them).

Essentially, you just need to go to BioMart, and

1. select "Ensembl Genes 96" (the number will change if the version changes) as the database and "uman Genes" as the dataset.




2. Click on "Filters", and set Gene type to coding and Transcript type to protein_coding.




3. From "Attributes", select whatever you want to see. The "CDS Length" is under "Structures".

share|improve this answer

edited Apr 30 at 22:23

answered Apr 30 at 15:48

terdon♦terdon

4,9902830

$endgroup$

add a comment | 

$begingroup$

While I haven't found a way to limit the results to the canonical transcript only, you can get a list of genes, transcripts and their CDS lengths using Ensemble's BioMart. I have already set it up for you, you can see the results, and modify them, here (click on the "Results" link if you don't see them).

Essentially, you just need to go to BioMart, and

1. select "Ensembl Genes 96" (the number will change if the version changes) as the database and "uman Genes" as the dataset.




2. Click on "Filters", and set Gene type to coding and Transcript type to protein_coding.




3. From "Attributes", select whatever you want to see. The "CDS Length" is under "Structures".

share|improve this answer

edited Apr 30 at 22:23

answered Apr 30 at 15:48

terdon♦terdon

4,9902830

$endgroup$

share|improve this answer

edited Apr 30 at 22:23

answered Apr 30 at 15:48

terdon♦terdon

4,9902830

answered Apr 30 at 15:48

terdon♦terdon

4,9902830

answered Apr 30 at 15:48

terdon♦terdon

4,9902830

0

$begingroup$

Ensembl has an FTP site that allows you to select and download only the coding sequences from many different genomes.
https://useast.ensembl.org/info/data/ftp/index.html

To determine the length of those sequences, download the associated gtf or gff3 annotation file. The annotation file is tab delim. The fourth and fifth column represent genomic loci of the annotated region. Subtract the amount in the fourth column from the amount in the fifth column to yield the length of all the annotated features.
You can easily do this in R after loading the file in the environment using the Magrittr library. The following code will create a new column called gene_length with the associated gene lengths.

install.packages("magrittr")
# this only needs to be done once
library(magrittr)
# must be run each time the library is neaded
annotation.gtf <- read.table("path/to/annotation.gtf")
annotation.gtf$start <- annotation.gtf[,4]
annotation.gtf$end <- annotation.gtf[,5]
annotation.new-column.gtf <- annotation.gtf %>%
 mutate(gene_length=end-start)

share|improve this answer

edited May 1 at 13:03

Kamil S Jaron

3,082942

answered May 1 at 11:33

Drew J-HDrew J-H

12

$endgroup$

add a comment | 

0

$begingroup$

Ensembl has an FTP site that allows you to select and download only the coding sequences from many different genomes.
https://useast.ensembl.org/info/data/ftp/index.html

To determine the length of those sequences, download the associated gtf or gff3 annotation file. The annotation file is tab delim. The fourth and fifth column represent genomic loci of the annotated region. Subtract the amount in the fourth column from the amount in the fifth column to yield the length of all the annotated features.
You can easily do this in R after loading the file in the environment using the Magrittr library. The following code will create a new column called gene_length with the associated gene lengths.

install.packages("magrittr")
# this only needs to be done once
library(magrittr)
# must be run each time the library is neaded
annotation.gtf <- read.table("path/to/annotation.gtf")
annotation.gtf$start <- annotation.gtf[,4]
annotation.gtf$end <- annotation.gtf[,5]
annotation.new-column.gtf <- annotation.gtf %>%
 mutate(gene_length=end-start)

share|improve this answer

edited May 1 at 13:03

Kamil S Jaron

3,082942

answered May 1 at 11:33

Drew J-HDrew J-H

12

$endgroup$

add a comment | 

$begingroup$

Ensembl has an FTP site that allows you to select and download only the coding sequences from many different genomes.
https://useast.ensembl.org/info/data/ftp/index.html

To determine the length of those sequences, download the associated gtf or gff3 annotation file. The annotation file is tab delim. The fourth and fifth column represent genomic loci of the annotated region. Subtract the amount in the fourth column from the amount in the fifth column to yield the length of all the annotated features.
You can easily do this in R after loading the file in the environment using the Magrittr library. The following code will create a new column called gene_length with the associated gene lengths.

install.packages("magrittr")
# this only needs to be done once
library(magrittr)
# must be run each time the library is neaded
annotation.gtf <- read.table("path/to/annotation.gtf")
annotation.gtf$start <- annotation.gtf[,4]
annotation.gtf$end <- annotation.gtf[,5]
annotation.new-column.gtf <- annotation.gtf %>%
 mutate(gene_length=end-start)

share|improve this answer

edited May 1 at 13:03

Kamil S Jaron

3,082942

answered May 1 at 11:33

Drew J-HDrew J-H

12

$endgroup$

install.packages("magrittr")
# this only needs to be done once
library(magrittr)
# must be run each time the library is neaded
annotation.gtf <- read.table("path/to/annotation.gtf")
annotation.gtf$start <- annotation.gtf[,4]
annotation.gtf$end <- annotation.gtf[,5]
annotation.new-column.gtf <- annotation.gtf %>%
 mutate(gene_length=end-start)

share|improve this answer

edited May 1 at 13:03

Kamil S Jaron

3,082942

answered May 1 at 11:33

Drew J-HDrew J-H

12

edited May 1 at 13:03

Kamil S Jaron

3,082942

edited May 1 at 13:03

Kamil S Jaron

3,082942

answered May 1 at 11:33

Drew J-HDrew J-H

12

answered May 1 at 11:33

Drew J-HDrew J-H

12