Questions tagged [chip-seq]

The tag has no usage guidance.

Filter by
Sorted by
Tagged with
0 votes
0 answers
21 views

Error loading DiffBind sample sheet

I am trying to analyse a ChIP-Seq data set with DiffBind. These are the contents of my sample sheet ("diffbind_samples.csv"): ...
user avatar
3 votes
1 answer
50 views

Where do I find the primary evidence for regulatory relationships included in RegulonDB?

When I look up high-throughout (ChIP) evidence for binding of the ferric uptake regulator Fur to sites in the genome of E. coli, RegulonDB lists a whopping 4214 targets. (To repeat this search, go to ...
user avatar
1 vote
0 answers
117 views

Metagene analysis with deepTools; 5UTR-CDS-3UTR and Exon-Intron-Exon

I have some PCR-CLIP data. PAR-CLIP identifies targets for RNA binding proteins. I am interested to build metagene profile so that transcript wide binding pattern can be compared. What I am mostly ...
user avatar
  • 11
1 vote
0 answers
130 views

How to choose the proper way to analyze the differential binding peaks among different samples for ChIP-seq? DiffBind or CSAW?

I have used the TMM method in CSAW package to normalize the composition bias for my different ChIP-seq samples and also got the normalized bigwig files. Next step I want to analyze the differential ...
user avatar
0 votes
1 answer
49 views

Impact of merging ChIP-seq runs of the same sample on PCR duplicates identification?

I'd like to a follow up question to this question related to merging fastq files for ChIP-seq. Let's assume we have one sample library that is re-sequenced two or three times in order to achieve a ...
user avatar
  • 1
1 vote
2 answers
147 views

Problems about value log2(IP/Input) less than zero in ChIP-seq?

I have a question about the normalization for ChIP-seq. I used CPM to normalize my bam files of each IP and Input. Then I calculate the coverage of gene bodies for all genes on the genome. I have WT ...
user avatar
0 votes
1 answer
41 views

How to interpret track height in Chipseq?

I recently received some .bigwig files from a chipseq experiment. I loaded the files into IGV. The interpretation is rather intuitive. My question is the following: What exactly do the track heights ...
user avatar
0 votes
0 answers
20 views

How to identify specific motif for chip-seq data

For chip-seq data, I had analyzed the samples until motif discovery using both cluster profiler and Homer for each histone mark. I was wondering what is the efficient method to determine if each GO ...
user avatar
  • 95
1 vote
1 answer
93 views

chipseeker annotation issue

...
user avatar
  • 1,640
0 votes
1 answer
81 views

Chip-seq down-sampling problem

Control sample chip read stat ...
user avatar
  • 1,640
1 vote
1 answer
298 views

annotation using ChIPseeker package error

I have differential binding sites object obtain from diffBind (dba.report). I am using the ChIP Seeker package to annotate the peaks but keep getting the following error: Error in (function (classes, ...
user avatar
  • 95
0 votes
1 answer
24 views

Extracting concatenated metadata from a column in R

I've downloaded some ChIP-seq data from www.chip-atlas.org and I want to use this to conduct some enrichment tests with genomic regions from my experiment. The .bed file generated by that website is ...
user avatar
0 votes
1 answer
42 views

How can I add for several bed files the header : track type=narrowPeak name=“narrowPeak” preferably in python ,can handle with R

I want to create custom tracks from these files I can add the line : track type=narrowPeak name=“narrowPeak” manually by opening it with text editor: track type=narrowPeak name=“narrowPeak” but I ...
user avatar
0 votes
1 answer
157 views

How to visualize called narrowPeak files in UCSC Genome browser? [closed]

I have this file: and I get this error I googled and found this solution but it still doesn't work and I get this error Error File 'GSM2797523_FOXA2_IDR0.02_narrowPeak_try.bed' - line 1 of custom ...
user avatar
0 votes
2 answers
95 views

Recommended coverage depth for ChIP Seq (H3K9me3 )

According to ENCODE recommendations For broad-peak histone experiments, each replicate should have 45 million usable fragments According to my understanding, a usable fragment means a uniquely ...
user avatar
  • 167
0 votes
2 answers
474 views

Installing HiC-Pro issues

I am setting up a HiChIPseq pipeline, but installing software necessary for the HiC side of things has been unbecoming. Both my PI and I have tried to install Hi-C Pro onto our Mac Pro that we use as ...
user avatar
  • 177
0 votes
2 answers
37 views

Separating peaks of chip-seq with specific length

My data contain several chip-seq results. I have the peaks called by MACS2.I wanted to only look at those peaks that their size is e.g 500bp to 1000bp. How can I separate those peaks efficiently? I ...
user avatar
  • 95
2 votes
1 answer
26 views

How to calculate the number of peaks that are upstream/downstream of some other peaks

I have 3 histone marks,I have used Macs2 for peak calling and diffBind to analyze the peaks. I was wondering if you know any way to calculate the peak numbers of one specific histone mark that occur ...
user avatar
  • 95
1 vote
1 answer
42 views

DiffBind, diffferentially binding site

I have data for 3 histone marks (2 for silencing and 1 for activation) each mark has three replicates. when I run the diffBind package I have three contrast: ...
user avatar
  • 95
0 votes
3 answers
129 views

differential analysis of chip-seq data

I have several sets of chip-seq data. I called the peaks using Macs2. I am pretty new to the field and I will appreciate any help. I wanted to annotate the peaks and see which peaks are shared between ...
user avatar
  • 95
0 votes
2 answers
465 views

Hypergeometic test for the overlap of two set of genomic intervals (e.g. from Chip-seq data)

Suppose I have done two Chip-seq experiments. Now after the bias correction step such as MACS2, I get two set of genomic intervals (or peaks, as what is usually called) that corresponds to each of the ...
user avatar
  • 43
1 vote
1 answer
100 views

merge the matrix from computematrix of deeptools

I have several matrix from computematrix of deeptools. I need to merge two of them using "computeMatrixOperations cbind -m input1.mat.gz input2.mat.gz -o output.mat.gz" but I am running to error "/...
user avatar
  • 95
0 votes
2 answers
164 views

How to make Venn diagram for the Macs peak calling output of two data sets?

I have two outputs of macs2 peak calling for two of my data sets. I wanted to plot the Venn diagram to see how many peaks are shared. I mainly work on Unix/Linux. Do you know any way that I can have ...
user avatar
  • 95
1 vote
1 answer
172 views

plot correlation between several bam files

I have several mab files from medip-seq. they supposed to be replicates and I have to draw their correlation. First I used multiBamSummary from deeptools to merge all the bam files based on bi/bed. ...
user avatar
  • 95
2 votes
1 answer
63 views

what is the concept of N-channel geneChip?

I can understand two-color and one-color(Channel) gene chip; but there is another type of gene-Chip called "N-channel chips" ; I do not know how can several (> 2) samples be hybridize in one chip. I ...
user avatar
1 vote
1 answer
336 views

How to get the intersection of peaks after peak calling using MACS2?

I have following 5 narrow peak files after peak calling. K14_peaks.narrowPeak K15_peaks.narrowPeak K16_peaks.narrowPeak K3_peaks.narrowPeak K8_peaks.narrowPeak I ...
user avatar
0 votes
1 answer
52 views

Pol_II genome wide occupancy analysis

I would like to know how to perform this global Pol II occupancy which is shown here: Paper Figure 1G So in this figure they are analysing the Pol II occupancy in Control,PAF1 KD and LEO1. I have ...
user avatar
  • 1,640
1 vote
1 answer
132 views

When running ChIPQC bioconductor package after peak calling, should I use sorted BAM files or unsorted BAM files?

I have called peaks using MACS2 and now I'm trying to run ChIPQC bioconductor package. Can someone please tell me should I use sorted BAM files or unsorted BAM files? I have created sorted BAM files ...
user avatar
1 vote
1 answer
1k views

How to visualize called narrowPeak files in UCSC Genome browser or IGV?

I have called peaks using MACS2. Then I got a narrowPeak file like this. ...
user avatar
5 votes
1 answer
48 views

Find covariation of different DNA-binding protein binding between two conditions

I have ChIP-seq data of RNAPII and two other factors which we think follow RNAPII during transcription in two different conditions, and I'd like to show that the genes that lose RNAPII also lose the ...
user avatar
  • 280
0 votes
1 answer
455 views

Tag density plot form chip or atac seq data

I m quoting this paper , figure B where they write "B) Violin plots for the normalized ATAC-seq tag density for all peaks within the indicated OC/CO groups at the indicated time points. Data were ...
user avatar
  • 1,640
2 votes
2 answers
125 views

Is there any way to align ChIP-seq reads to telomeres?

I know that telomeres are highly repeated sequences, but is there any way to retain any reads that map to these regions (on HG38)? I recently managed to find some protein binding to centromeres, ...
user avatar
  • 280
2 votes
1 answer
121 views

RPGC normalisation creates artefacts at centromere

I am studying ChIP-Seq data in HeLa cells and I've started using the RPGC normalisation of deepTool's bamCoverage. MACS2 also uses this normalisation for its peak calling. I am seeing a large number ...
user avatar
  • 280
0 votes
1 answer
595 views

How to convert a .fastq chip-seq file into a wig/bigwig track and how to extract chip-seq peaks?

I have a reference .fasta file and a raw .fastq file with chip-seq data. I am trying to create a bigwig track from the and .fastq and .fasta ref file of the raw signal. Then I would like to do some ...
user avatar
  • 1,347
0 votes
1 answer
466 views

Error: package or namespace load failed for 'ChIPseeker'

How do I solve the problem of loading the bioconductor Chipseeker package? I have installed in the following way and I having an error while loading it and I do not know why. ...
user avatar
  • 497
1 vote
2 answers
557 views

Biological replicates on Chip-seq Transcription factor data

I have 2 biological replicates for chip-seq transcription factor data and I have applied macs2 peak calling for each replicate separately. How can i make use of the biological replicates to extract ...
user avatar
3 votes
1 answer
93 views

Occupancy of TFs with the target genes

The occupancy of SMARCD3 in the target genes listed below. I want to see average, normalized ChIP-seq signal at the promoter proximal region (1000bp upstream and downstream of the TSS). I have 4 ...
user avatar
  • 1,640
1 vote
0 answers
126 views

How can I convert Affymetrix microarray data to bedgraph format?

I've previously asked a different question about the ChIP-chip microarray data that I'm working on. I'm new to this so I'm sorry for asking here again. I have .CEL files and I'm using the Affy ...
user avatar
  • 11
1 vote
1 answer
101 views

How to convert ChiP-seq track in .bam file format into an array at fixed bin resolution?

Is there a out of the box to convert a .bam file to a .wig file at specific resolution/bin size? I was thinking to use igvtools, but I am not sure how it works: ...
user avatar
  • 1,347
1 vote
1 answer
108 views

Score predefined ChIP-seq peaks with MACS2 or equivalent

I have performed peak calling on a number of separate ChIP-seq experiments and would like to harmonise the peaks from each of the experiments in order to convert my data into a matrix for further ...
user avatar
3 votes
1 answer
888 views

Integrative analysis of omics studies using machine learning

I would like to use public omics datasets (ChIP-seq, RNA-seq, and ATAC-seq) from different studies to do an integrative analysis as follow: Normalise samples, within each type of omics, from ...
user avatar
  • 33
1 vote
0 answers
96 views

peakranger bcp segmentation fault

I'm running into a segmentation fault running peakranger's bcp on Arch Linux. Binaries were compiled from source just today. Any idea how I can prevent this? ...
user avatar
3 votes
2 answers
1k views

Input normalization in ChIP-seq

If I subtract input counts from ChIP counts (for every gene, since I have one peak per gene) I get negative values for most genes. This makes sense to me, because (as can be seen in the figure) input ...
user avatar
  • 547
3 votes
1 answer
738 views

determine if ChIP-seq peaks are broad or narrow

Is there a method to determine if the peaks are broad or narrow? ENCODE provides some guidelines: Although those cover common histone marks, there are many others. If you are using one of the ones ...
user avatar
  • 2,099
5 votes
1 answer
1k views

Macs2 peak calling?

I have paired end ChIP-seq data with 101 bp and 2 biological replicates for each one. I have done peak calling with macs2 but I have some questions about it. I also faced with an warning: WARNING @...
user avatar
  • 153
1 vote
1 answer
1k views

What are phantom peaks in ChIP-seq?

I've seen two uses of this term that seem to refer to completely different phenomena. Is my understanding correct? In Jain et al. "Active promoters give rise to false positive ‘Phantom Peaks’ in ChIP-...
user avatar
1 vote
0 answers
74 views

How long does it typically take to call peaks via `phantompeakqualtools`?

I'm running the script run_spp_nodups.R from phantompeakqualtools. The input is a pair of BED 3+3 ("tagalign") files that are ...
user avatar
3 votes
1 answer
651 views

Aligning ChIP-Seq reads to repeats for downstream peak analysis

This is a brief question regarding the above. I have previously used bowtie to map reads from paired-end ChIP-Seq sequencing, then used the positions for peak-calling. However, I'm trying to do ...
user avatar
  • 291
7 votes
1 answer
1k views

Where to download JASPAR TFBS motif bed file?

I am interested in determining if any transcription factor binding site motifs are enriched in some BED files from a DNA methylation experiment. I am looking for a database that has BED Files ...
user avatar
  • 367
4 votes
2 answers
778 views

Where are .motif files from homer knownResults?

I have been using homer's findMotifsGenome.pl, but with my new version (v4.9.1) of homer I don't get ...
user avatar
  • 3,541