Questions tagged [chip-seq]
The chip-seq tag has no usage guidance.
56
questions
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"):
...
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 ...
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 ...
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 ...
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 ...
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 ...
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 ...
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 ...
1
vote
1
answer
93
views
0
votes
1
answer
81
views
Chip-seq down-sampling problem
Control sample chip read stat
...
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, ...
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 ...
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 ...
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 ...
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 ...
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 ...
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 ...
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 ...
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:
...
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 ...
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 ...
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 "/...
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 ...
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. ...
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 ...
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 ...
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 ...
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 ...
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.
...
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 ...
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 ...
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, ...
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 ...
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 ...
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.
...
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 ...
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 ...
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 ...
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:
...
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 ...
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 ...
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?
...
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 ...
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 ...
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 @...
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-...
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 ...
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 ...
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 ...
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 ...