Hi everyone,

I'm wondering where could I find the number of protein-coding genes in each chromosome arm in human? I have searched the ensembl, but it only gives the the number of protein-coding genes in each chromosome instead if dividing with arms. Wondering if anyone could give any suggestions or if I could get the information from using biomart or something else?

Thank you!

Something like this

> library(biomaRt)
> mart <- useEnsembl("ensembl","hsapiens_gene_ensembl")
> z <- getBM(c("ensembl_gene_id", "band","chromosome_name"), c("chromosome_name","biotype"), list(c(1:22, "X","Y"), "protein_coding"), mart = mart)
> z$pq <- paste(z$chromosome_name, sapply(strsplit(z$band, ""), "[", 1), sep = "_")
> sapply(split(z$ensembl_gene_id, z$pq), length)
 1_p  1_q 10_p 10_q 11_p 11_q 12_p 12_q 13_q 14_q 15_q 16_p 16_q 17_p 17_q 18_p 
1078  983  153  577  449  868  291  746  322  616  601  475  382  346  840   66 
18_q 19_p 19_q  2_p  2_q 20_p 20_q 21_p 21_q 22_q  3_p  3_q  4_p  4_q  5_p  5_q 
 200  641  835  470  777  186  360    1  219  447  510  571  234  523  155  727 
 6_p  6_q  7_p  7_q  8_p  8_q  9_p  9_q  X_p  X_q  Y_p  Y_q 
 619  432  297  632  283  418  214  563  335  524   14   31