Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Background set for GO enrichment analysis yvan.wenger Bioinformatics 1 01-10-2018 11:41 AM
Issue with BiSeq for DMR analysis shocker8786 Bioinformatics 13 10-12-2016 04:41 PM
SAAP-BS data analysis (RRBS) - BiSeq use and Coverage Histogram Maria_Maqueda Epigenetics 0 04-28-2015 07:41 AM
Which software/algorithm do you use to find DMR in chipseq methylation data? KYR Bioinformatics 0 10-10-2014 10:06 AM

Thread Tools
Old 03-05-2019, 12:18 PM   #1
Location: Sheffield, UK

Join Date: Feb 2012
Posts: 18
Default Background for GO term enrichment analysis of RRBS DMR data?


I am attempting to carry out GO term enrichment analysis on set of differentially methylated genomic regions (DMRs). The program I am attempting to do this in is DAVID.

My data is derived from reduced-representation bisulfite sequencing (RRBS) of liver tissue. Although I am working on a non-model species, there is relatively close reference genome available, and I have used this to annotate the DMRs where they overlap with known genes.

I therefore have a gene list that I would like to use in GO term enrichment analysis and explore other aspects of functional annotation. However, when carrying out these sorts of analyses in DAVID, a background list is required for statistical comparison. This usually involves taking a list of the genes known from that particular reference genome (e.g. 30,000 in humans) and then seeing if any particular gene category are over-represented in my gene list in comparison to the genomic background. With the RRBS data, I am only sequencing a subset of the available genome and genes in my gene list can only come from this subset, therefore to ask if genes categories are over-represented in my DMR data set by comparing against the entire number of genes known from that organism does not really make much sense to me.

One way of creating a background list might be to get a list of genes that overlap with any aligned methylated C in my data set (not just the significant ones) and then use this as the background to do the test. Would this be appropriate? One complication is that I am comparing three different pairwise treatments, so I guess I would need three independent background lists comprising methylated Cs present in at least one individual/replicate in both treatment groups being compared or is that nonsense?

More generally, does anyone have any advice on how to generate an appropriate background list for RRBS functional annotation analysis? I guess this might overlap with other reduced-representation sequencing techniques such as RAD-seq. Alternatively, does anyone know of any other methods/programs for carrying out GO term enrichment analysis that takes into account the biased sampling of the genome involved in RRBS?

Thanks in advance for any advice you can give,

pig_raffles is offline   Reply With Quote

annotation, epigenetics, go terms, methylation, rrbs

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 02:39 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO