Locating nearest genes in RRBS analysis
Entering edit mode
Mageiros • 0
Last seen 12 days ago

I am trying to analyze a data set of BS-seq (RRBS) data following the example at the section 4.7 of the edgeR user manual.

In the tutorial the function nearestTSS is used to annotate the nearest gene of each site.

Is there a way to retrieve information though for the two upstream and two downstream genes instead of the one that is closer?

I found the get2NearestFeature function (https://rdrr.io/bioc/hiAnnotator/man/get2NearestFeature.html) but I am not sure how to convert my DGElist object to a GRanges objects.

Note that I am using hg19 as a human reference strain.

Thank you very much in advance for your time and help.

Best Leonardos

hiAnnotator bsseq DifferentialMethylation edgeR RRBSdata • 69 views
Entering edit mode
Last seen 56 minutes ago
WEHI, Melbourne, Australia

You can form a GRanges object from the vector of chromosome names and the vector of locations, same as you input to nearestTSS. The help page you link to gives an example of how to do that. I'm not familiar with get2NearestFeature though so can't help you with it.

Regarding nearestTSS, please be aware that it only works with hg38, not with hg19, so you should not use it at all in your analysis.

Is there a reason for using hg19, since hg38 was released nearly a decade ago?


Login before adding your answer.

Traffic: 147 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6