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

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?