Question: ChiPQC and paired-end data
0
gravatar for Ugo Borello
5.4 years ago by
Ugo Borello340
France
Ugo Borello340 wrote:
Dear Rory, I have few basic questions about ChiPQC. - Could I use ChiPQC with paired-end data? - I would like to use my mouse blacklist https://sites.google.com/site/anshulkundaje/projects/blacklists MOUSE (mm9): http://www.broadinstitute.org/~anshul/projects/mouse/blacklist/mm9-bla cklist .bed.gz I noticed that the formatting of this mm9-blacklist.bed is different from your blacklist_hg19.bed. For example your file has an header, not mine. Is it correct then to provide only a filename to my blacklist, when calling ChIPQC(), or is it better to import the bed file as Granges object first? - If I understand correctly, while constructing a ChIPQCexperiment object and setting consensus=TRUE, a consensus peak set is generated. This consensus peak set consists of the peaks overlapping at least two of the replicate samples, correct? Is it possible to access this consensus peak set? Thank you for your help and for this very useful tool, Ugo
chipqc • 885 views
ADD COMMENTlink modified 5.3 years ago by Carroll, Thomas20 • written 5.4 years ago by Ugo Borello340
Answer: ChiPQC and paired-end data
0
gravatar for Carroll, Thomas
5.3 years ago by
Carroll, Thomas20 wrote:
Hi Ugo, Sorry for the delay in reply. If you include the maintainers in the post then I will get a quicker notification. I usually supply blacklists as GRanges objects and this is what i have done when using mm9 blacklists myself form Anshul Kudaje's site. ChIPQC expects there to be a header when reading from a file. You could use paired-end data with ChIPQC but for now this would require you to filter your Bam file for only "first read in pair" in order to avoid double counting and biasing the cross-coverage scores plots. This should be pretty straightforward to do using samtools view. In the future we will include paired-end data analysis for ChIPQC. Not sure how Diffbind handles paired end data, perhaps Rory has an idea? To retrieve consensus sets you could take the DBA object from ChIPQCexperiment and work from there using the Diffbind package. dbaobject = ChIPQCexperiment at DBA Then simply follow DiffBind vignette to select a suitable consensus set. best, Tom Thomas Carroll Head of Bioinformatics MRC Clinical Sciences Centre, Faculty of Medicine, Imperial College London,Hammersmith Hospital Campus Du Cane Road,London,W12 0NN
ADD COMMENTlink written 5.3 years ago by Carroll, Thomas20
Thank you Thomas, > You could use paired-end data with ChIPQC but for now this would require you > to filter your Bam file for only "first read in pair" in order to avoid double > counting and biasing the cross-coverage scores plots. > This should be pretty straightforward to do using samtools view. To be sure I understand you mean: $ samtools view -f 64 > To retrieve consensus sets you could take the DBA object from ChIPQCexperiment > and work from there using the Diffbind package. > > dbaobject = ChIPQCexperiment at DBA Indeed, I finally got my consensus playing with the Diffbind package. I needed to understand how to link properly ChIPQC with Diffbind. Thank you again for your help. Ugo > From: "Carroll, Thomas" <thomas.carroll at="" csc.mrc.ac.uk=""> > Date: Mon, 28 Jul 2014 12:16:15 +0000 > To: "bioconductor at r-project.org" <bioconductor at="" r-project.org=""> > Cc: "ugo.borello at inserm.fr" <ugo.borello at="" inserm.fr=""> > Subject: Re: [BioC] ChiPQC and paired-end data > > Hi Ugo, > > Sorry for the delay in reply. If you include the maintainers in the post then > I will get a quicker notification. > > > I usually supply blacklists as GRanges objects and this is what i have done > when using mm9 blacklists myself form Anshul Kudaje's site. ChIPQC expects > there to be a header when reading from a file. > > You could use paired-end data with ChIPQC but for now this would require you > to filter your Bam file for only "first read in pair" in order to avoid double > counting and biasing the cross-coverage scores plots. > This should be pretty straightforward to do using samtools view. > > In the future we will include paired-end data analysis for ChIPQC. Not sure > how Diffbind handles paired end data, perhaps Rory has an idea? > > > > > To retrieve consensus sets you could take the DBA object from ChIPQCexperiment > and work from there using the Diffbind package. > > dbaobject = ChIPQCexperiment at DBA > > > Then simply follow DiffBind vignette to select a suitable consensus set. > > > best, > > Tom > > > Thomas Carroll > Head of Bioinformatics > MRC Clinical Sciences Centre, Faculty of Medicine, > Imperial College London,Hammersmith Hospital Campus > Du Cane Road,London,W12 0NN
ADD REPLYlink written 5.3 years ago by Ugo Borello340
Hi Ugo, >To be sure I understand you mean: >$ samtools view -f 64 That sounds right to get the first in pair. Cheers, tom Thomas Carroll Head of Bioinformatics MRC Clinical Sciences Centre, Faculty of Medicine, Imperial College London,Hammersmith Hospital Campus Du Cane Road,London,W12 0NN ________________________________________ From: Ugo Borello [ugo.borello@inserm.fr] Sent: Wednesday, July 30, 2014 3:05 PM To: Carroll, Thomas; bioconductor at r-project.org Subject: Re: [BioC] ChiPQC and paired-end data Thank you Thomas, > You could use paired-end data with ChIPQC but for now this would require you > to filter your Bam file for only "first read in pair" in order to avoid double > counting and biasing the cross-coverage scores plots. > This should be pretty straightforward to do using samtools view. To be sure I understand you mean: $ samtools view -f 64 > To retrieve consensus sets you could take the DBA object from ChIPQCexperiment > and work from there using the Diffbind package. > > dbaobject = ChIPQCexperiment at DBA Indeed, I finally got my consensus playing with the Diffbind package. I needed to understand how to link properly ChIPQC with Diffbind. Thank you again for your help. Ugo > From: "Carroll, Thomas" <thomas.carroll at="" csc.mrc.ac.uk=""> > Date: Mon, 28 Jul 2014 12:16:15 +0000 > To: "bioconductor at r-project.org" <bioconductor at="" r-project.org=""> > Cc: "ugo.borello at inserm.fr" <ugo.borello at="" inserm.fr=""> > Subject: Re: [BioC] ChiPQC and paired-end data > > Hi Ugo, > > Sorry for the delay in reply. If you include the maintainers in the post then > I will get a quicker notification. > > > I usually supply blacklists as GRanges objects and this is what i have done > when using mm9 blacklists myself form Anshul Kudaje's site. ChIPQC expects > there to be a header when reading from a file. > > You could use paired-end data with ChIPQC but for now this would require you > to filter your Bam file for only "first read in pair" in order to avoid double > counting and biasing the cross-coverage scores plots. > This should be pretty straightforward to do using samtools view. > > In the future we will include paired-end data analysis for ChIPQC. Not sure > how Diffbind handles paired end data, perhaps Rory has an idea? > > > > > To retrieve consensus sets you could take the DBA object from ChIPQCexperiment > and work from there using the Diffbind package. > > dbaobject = ChIPQCexperiment at DBA > > > Then simply follow DiffBind vignette to select a suitable consensus set. > > > best, > > Tom > > > Thomas Carroll > Head of Bioinformatics > MRC Clinical Sciences Centre, Faculty of Medicine, > Imperial College London,Hammersmith Hospital Campus > Du Cane Road,London,W12 0NN
ADD REPLYlink written 5.3 years ago by Carroll, Thomas20
Answer: ChiPQC and paired-end data
0
gravatar for Tom Carroll
5.3 years ago by
Tom Carroll10
Tom Carroll10 wrote:
Ugo Borello <ugo.borello at="" ...=""> writes: > > Dear Rory, > > I have few basic questions about ChiPQC. > > - Could I use ChiPQC with paired-end data? > > - I would like to use my mouse blacklist > https://sites.google.com/site/anshulkundaje/projects/blacklists > MOUSE (mm9): > http://www.broadinstitute.org/~anshul/projects/mouse/blacklist/mm9-b lacklist > .bed.gz > > I noticed that the formatting of this mm9-blacklist.bed is different from > your blacklist_hg19.bed. For example your file has an header, not mine. > Is it correct then to provide only a filename to my blacklist, when calling > ChIPQC(), or is it better to import the bed file as Granges object first? > > - If I understand correctly, while constructing a ChIPQCexperiment object > and setting consensus=TRUE, a consensus peak set is generated. This > consensus peak set consists of the peaks overlapping at least two of the > replicate samples, correct? > Is it possible to access this consensus peak set? > > Thank you for your help and for this very useful tool, > > Ugo > > _______________________________________________ > Bioconductor mailing list > Bioconductor at ... > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor > > Hi Ugo, Sorry for the delay in reply. Supplying a GRanges object should work well and typically what I do. ChIPQC expects there to be a header when reading from a file so your file should be good. You could use paired-end data with ChIPQC but for now this would require you to filter your Bam file for only "first read in pair" in order to avoid double counting and biasing the cross-coverage scores plots. This should be straightforward to do using something like samtools view. In the future I do want to include paired-end data analysis for ChIPQC to accommodate MNAse-seq. I'm not sure how Diffbind handles paired end data, perhaps Rory has an idea? To retrieve consensus sets you could take the DBA object from ChIPQCexperiment and work from there using the Diffbind package. dbaobject = ChIPQCexperimentObject at DBA best, Tom Thomas Carroll Head of Bioinformatics, MRC Clinical Sciences Centre, Faculty of Medicine, Imperial College London,Hammersmith Hospital Campus, Du Cane Road,London,W12 0NN
ADD COMMENTlink written 5.3 years ago by Tom Carroll10
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: 327 users visited in the last hour