converting transcripts to gene level
1
0
Entering edit mode
@abhishek-pratap-5083
Last seen 9.6 years ago
Hi All I have cufflinks generated transcripts file in gtf format. I wondering if I can use any R package to convert it to gene level gtf file. By this means all the transcripts are collapsed into genes and each gene contains all the exons seen in all the transcripts for that gene. Sample Data attached. In the attached GTF file I have 3 made up transcripts from 2 genes. What I?would want to get as a output is a gtf file containing two genes and all the exons associated with them. Thanks! -Abhi
convert convert • 2.9k views
ADD COMMENT
0
Entering edit mode
@michael-lawrence-3846
Last seen 2.3 years ago
United States
Hi, There's nothing that does exactly that. You can use rtracklayer::import(gtf, format = "gff2") to get the file into a RangedData or GRanges. And then after your manipulations write it out again with export(object, file, "gff2"). Michael On Wed, Feb 1, 2012 at 5:15 PM, Abhishek Pratap <apratap@lbl.gov> wrote: > Hi All > > I have cufflinks generated transcripts file in gtf format. I wondering > if I can use any R package to convert it to gene level gtf file. By > this means all the transcripts are collapsed into genes and each gene > contains all the exons seen in all the transcripts for that gene. > > > Sample Data attached. > > In the attached GTF file I have 3 made up transcripts from 2 genes. > What I would want to get as a output is a gtf file containing two > genes and all the exons associated with them. > > Thanks! > -Abhi > > _______________________________________________ > Bioconductor mailing list > Bioconductor@r-project.org > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: > http://news.gmane.org/gmane.science.biology.informatics.conductor > [[alternative HTML version deleted]]
ADD COMMENT
0
Entering edit mode
You may want to check out the cummeRbund package. It doesn't do the conversion you want in but offers other storage and manipulation options of Cufflinks data that may be of interest. It's written by the authors of TopHat and Cufflinks. Valerie On 02/02/2012 04:42 AM, Michael Lawrence wrote: > Hi, > > There's nothing that does exactly that. You can use > rtracklayer::import(gtf, format = "gff2") to get the file into a RangedData > or GRanges. And then after your manipulations write it out again with > export(object, file, "gff2"). > > Michael > > On Wed, Feb 1, 2012 at 5:15 PM, Abhishek Pratap<apratap at="" lbl.gov=""> wrote: > >> Hi All >> >> I have cufflinks generated transcripts file in gtf format. I wondering >> if I can use any R package to convert it to gene level gtf file. By >> this means all the transcripts are collapsed into genes and each gene >> contains all the exons seen in all the transcripts for that gene. >> >> >> Sample Data attached. >> >> In the attached GTF file I have 3 made up transcripts from 2 genes. >> What I would want to get as a output is a gtf file containing two >> genes and all the exons associated with them. >> >> Thanks! >> -Abhi >> >> _______________________________________________ >> Bioconductor mailing list >> Bioconductor at r-project.org >> https://stat.ethz.ch/mailman/listinfo/bioconductor >> Search the archives: >> http://news.gmane.org/gmane.science.biology.informatics.conductor >> > [[alternative HTML version deleted]] > > _______________________________________________ > Bioconductor mailing list > Bioconductor at r-project.org > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
ADD REPLY

Login before adding your answer.

Traffic: 803 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6