Question: Extracting gene structures from a GFF3 file for use with annotateWithGeneParts()
0
gravatar for alan.g.hudson
14 months ago by
alan.g.hudson0 wrote:

Hi,

I want to use genomation to annotate previously identified differentially methylated genomic regions and quantify the overall degree of overlap with different gene structures (promoters, TSSs, exons, introns, intergenic).

 

The annotation file for my specific reference genome is in GFF3 format. Whilst it seems relatively straight forward to use BED 12 format to extract the relevant genomic coordinates using readTranscriptFeatures() - I can't figure out how to do the same for a GFF3 format annotation file.

 

I have looked at GenomicFeatures and there it seems possible to extract some of the features listed above (promoters, exons, introns) but this would result in a mix of several GRanges and GRanges list objects and again I have no idea how to combine these objects into the single object of the right format to plug into annotateWithGeneParts()

 

Has any one successfully dealt with a similar situation and can give some advice? Alternatively, if anyone can suggest another annotation strategy that would also be gratefully received?

 

Cheers,

 

Alan

ADD COMMENTlink modified 14 months ago by altuna akalin20 • written 14 months ago by alan.g.hudson0
Answer: Extracting gene structures from a GFF3 file for use with annotateWithGeneParts()
1
gravatar for altuna akalin
14 months ago by
Germany/Berlin/MDC-BIMSB
altuna akalin20 wrote:
I already answered this exact question here: https://groups.google.com/forum/#!topic/genomation/5EFpaBBEWLQ you can do that if you can create a GRanges list that has the same structure with readTranscriptFeatures() output. It returns a GRangesList most likely or simple list with GRanges (I think it has TSS, exon, intron and promoter GRanges objects). You just have to replicate the order and the structure of that output then downstream stuff from genomation should work. On Tue, Jul 31, 2018 at 10:32 PM alan.g.hudson [bioc] < noreply@bioconductor.org> wrote: > Activity on a post you are following on support.bioconductor.org > > User alan.g.hudson <https: support.bioconductor.org="" u="" 16729=""/> wrote Question: > Extracting gene structures from a GFF3 file for use with > annotateWithGeneParts() <https: support.bioconductor.org="" p="" 111565=""/>: > > Hi, > > I want to use genomation to annotate previously identified differentially > methylated genomic regions and quantify the overall degree of overlap with > different gene structures (promoters, TSSs, exons, introns, intergenic). > > > > The annotation file for my specific reference genome is in GFF3 format. > Whilst it seems relatively straight forward to use BED 12 format to extract > the relevant genomic coordinates using readTranscriptFeatures() - I can't > figure out how to do the same for a GFF3 format annotation file. > > > > I have looked at GenomicFeatures and there it seems possible to extract > some of the features listed above (promoters, exons, introns) but this > would result in a mix of several GRanges and GRanges list objects and again > I have no idea how to combine these objects into the single object of the > right format to plug into annotateWithGeneParts() > > > > Has any one successfully dealt with a similar situation and can give some > advice? Alternatively, if anyone can suggest another annotation strategy > that would also be gratefully received? > > > > Cheers, > > > > Alan > ------------------------------ > > Post tags: genomicfeatures, genomation, genomicranges, methylation, > annotation > > You may reply via email or visit > Extracting gene structures from a GFF3 file for use with annotateWithGeneParts() >
ADD COMMENTlink written 14 months ago by altuna akalin20
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 16.09
Traffic: 388 users visited in the last hour