estimation of non synonymous over synonymous (Ka/Ks) SNPs over genome with RNAseq

piemmea · 03-15-2013, 01:17 PM

dear all
I performed RNAseq in closely related plant species and was asked to perform a Ka/Ks estimation of the transcriptome. The problem is, that a number of genes are not fully covered, so I was wondering just to consider just those who are 100% covered. However, in some papers it seems that people look at single codons, not whole CDS, so the question is: do I consider the first option or the second is correct as well? Clearly, if I would consider also genes with incomplete coverage, I could get higher numbers. In case, does anyone can suggest me a software?
Many thanks
Marco

upendra_35 · 05-25-2014, 10:32 PM

I don't know if it is too late for you but i recently did something in the lines that you mentioned above - using RNA-Seq data to calculate Ka/Ks ratio for two genomes. For this first i have reconstructed the genomes of these species by incorporating the SNPs and then used KaKs calculator (https://code.google.com/p/kaks-calculator/) to calculate Ka/Ks ratio. Let me know if you need further details. Good luck!

Thanks
Upendra

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
synonymous snps from vcf	bioinfun	Bioinformatics	15	09-19-2017 06:41 AM
determining if a SNP is synonymous or not	nlsullivan	Bioinformatics	6	10-13-2013 05:38 AM
Nucleotide to amino acid translation + synonymous/non-synonymous count	Darth Ghisus	Bioinformatics	2	09-04-2012 04:15 AM
how to tell SNPs synonymous v. nonsynonymous	Margot	Bioinformatics	31	09-29-2011 09:24 AM
Non synonymous SNP programs	colindaven	Bioinformatics	9	05-15-2010 10:50 PM

03-15-2013, 01:17 PM	#1
piemmea Junior Member Location: Italy Join Date: Mar 2013 Posts: 6	estimation of non synonymous over synonymous (Ka/Ks) SNPs over genome with RNAseq dear all I performed RNAseq in closely related plant species and was asked to perform a Ka/Ks estimation of the transcriptome. The problem is, that a number of genes are not fully covered, so I was wondering just to consider just those who are 100% covered. However, in some papers it seems that people look at single codons, not whole CDS, so the question is: do I consider the first option or the second is correct as well? Clearly, if I would consider also genes with incomplete coverage, I could get higher numbers. In case, does anyone can suggest me a software? Many thanks Marco

05-25-2014, 10:32 PM	#2
upendra_35 Senior Member Location: USA Join Date: Apr 2010 Posts: 102	I don't know if it is too late for you but i recently did something in the lines that you mentioned above - using RNA-Seq data to calculate Ka/Ks ratio for two genomes. For this first i have reconstructed the genomes of these species by incorporating the SNPs and then used KaKs calculator (https://code.google.com/p/kaks-calculator/) to calculate Ka/Ks ratio. Let me know if you need further details. Good luck! Thanks Upendra