An effective approach for identification of in vivo protein–DNA binding sites from paired-end ChIP-Seq data
Ontology highlight
ABSTRACT:
ORGANISM(S): Arabidopsis thaliana
PROVIDER: GSE16940 | GEO | 2010/06/04
SECONDARY ACCESSION(S): PRJNA117657
REPOSITORIES: GEO
Similar Datasets
Project description:ChIP-Seq, which combines chromatin immunoprecipitation (ChIP) with high-throughput massively parallel sequencing, is increasingly being used for identification of proteinM-bM-^@M-^SDNA interactions in-vivo in the genome. In general, current algorithms for ChIP-seq reads employ artificial estimation of the average length of DNA fragments for peak finding, leading to uncertain prediction of DNA-protein binding sites. Here, we present SIPeS (Site Identification from Paired-end Sequencing), a novel algorithm for precise identification of binding sites from short reads generated from paired-end Solexa ChIP-Seq technology. SIPeS uses a dynamic baseline directly via M-bM-^@M-^Xpiling upM-bM-^@M-^Y the corresponding fragments defined by the paired reads to efficiently find peaks corresponding to binding sites. The performance of SIPeS is demonstrated by analyzing the ChIP-Seq data of the Arabidopsis basic helix-loop-helix transcription factor ABORTED MICROSPORES (AMS). The robustness of SIPeS was demonstrated in higher sensitivity and spatial resolution in peak finding compared to three existing peak detection algorithms. Keywords: transcription factors (protein-DNA interactions) Examination of protein-DNA interactions in buds of Arabidopsis anther cell
2010-06-04 | E-GEOD-16940 | biostudies-arrayexpress
Project description:CUT&RUN is a powerful tool to study protein-DNA interactions in vivo. DNA fragments cleaved by the targeted micrococcal nuclease identify the footprints of DNA-binding proteins on the chromatin. We performed CUT&RUN on human lung carcinoma cell line A549 maintained in a multi-well cell culture plate to profile RNA polymerase II. Size-fractionation of aligned reads revealed that the collection of long (>270 bp) DNA fragments corresponded to the bimodal peak around the transcription start sites, as previously seen with chromatin immunoprecipitation. However, we found that shorter (<120 bp) fragments identify a well-defined peak localised at the transcription start sites. This distinct DNA footprint of shorter fragments, which constituted only about 5% of the total reads, indicates the unique positioning of the poised RNA polymerase II before transcription initiation, which has not been detected in the steady state of transcription by standard chromatin immunoprecipitation.
2020-10-20 | GSE155666 | GEO
Project description:ChIP-Seq is a technique used to analyse protein-DNA interactions. The protein-DNA complex is pulled down using a protein antibody, after which sequencing and analysis of the bound DNA fragments is performed. A key bioinformatics analysis step is “peak” calling - identifying regions of enrichment. Benchmarking studies have consistently shown that no optimal peak caller exists. Peak callers have distinct selectivity and specificity characteristics which are often not additive and seldom completely overlap in many scenarios. In the absence of a universal peak caller, we rationalized one ought to utilize multiple peak-callers to 1) gauge peak confidence as determined through detection by multiple algorithms, and 2) more thoroughly survey the protein-bound landscape by capturing peaks not detected by individual peak callers owing to algorithmic limitations and biases. We therefore developed an integrated ChIP-Seq Analysis Pipeline (ChIP AP) which performs all analysis steps from raw fastq files to final result, and utilizes four commonly used peak callers to more thoroughly and comprehensively analyse datasets. Results are integrated and presented in a single file enabling users to apply selectivity and sensitivity thresholds to select the consensus peak set, the union peak set, or any sub-set in-between to more confidently and comprehensively explore the protein bound landscape. (https://github.com/JSuryatenggara/ChIP-AP).
2022-01-03 | GSE172355 | GEO
Project description:ChIP-seq has become the method of choice for studying functional DNA-protein interactions on a genome wide scale. The method is based on co-immunoprecipitation of DNA binding proteins with formaldehyde cross-linked DNA, followed by deep sequencing of immunoprecipitated chromatin fragments, allowing for the identification of binding sites with high accuracy [1-5]. Traditionally, genome-wide tiling path microarrays were used for the detection of specifically immunoprecipitated DNA fragments, but current next-generation sequencers now generate sufficient data to assay multiple samples in a single sequencing run, making this method more time and cost effective compared to microarray-based approaches [2]. However, due to limitations in read length of next generation sequencing technologies (50-76bp), it is not possible to sequence the complete length of immunoprecipitated DNA fragments which can be up to 2kb long reflecting biology in case of big protein complexes. As a consequence only the ends of the immunoprecipitated DNA fragments are sequenced. Deconvolution of sequencing reads mapping to the positive and negative strand is required to identify the real DNA binding site [4-9]. This limitation is most obvious in case of large complex regions where multiple binding sites of various regulatory elements are clustered close to each other. Deconvolution of such regions and the exact identification of individual binding positions is very challenging [4]. Similar problems can occur when studying histone positioning in cases where the length of ChIP fragments is bigger than the average distance of 2 neighboring nucleosomes. In addition frequently only narrow size range of immunoprecipitated fragments is selected for sequencing and thus possible bias towards binding regions within the selected range can be expected by loosing larger fragments possible originated from protein complexes interacting with DNA. To circumvent these complications, we modified the procedure for the preparation of ChIP-seq samples by introducing an additional extensive fragmentation round after isolation of immunoprecipitated chromatin to generate 70-110 bp long DNA fragments. For testing, we choose Tcf4 protein which is a well-studied downstream element of the Wnt-pathway for which the genome wide binding site profile is known was already known from ChIP-on-CHIP experiments [10]. Using the modified approach, we were able to identify Tcf4 binding sites at near nucleotide resolution without computational deconvolution. Furthermore, high-resolution information on genome-wide binding site regions allowed for the identification of potential novel Tcf4 co-factors as well as target genes. 5 samples + 3 input samples
2010-11-01 | E-GEOD-18481 | biostudies-arrayexpress
Project description:We analyzed levels of 5-methyl cytosine nnnn CCCGGG target sites by sequential restriction digest by SmaI and XmaI enzymes, ligating Illumina adaptors to the restriction fragments and reading methylation-specific signatures at the ends of restriction fragments by paired ends Illumina high throughput sequencing. Digital restriction enzyme analysis of methylation (DREAM) was performed to determine the methylation profile of SW48 colon cancer cell line genomic DNA. Genomic DNA spiked in with unmethylated, partially methylated and fully methylated standards was sequentially cut at CCCGGG sites with the methylation-sensitive enzyme SmaI (blunt ends) and its methylation-tolerant neoschizomer XmaI (5'CCGG overhangs), creating different end sequences that represented methylation status of the CCCGGG sites. These end sequences were analyzed by Illumina high throughput sequencing. Methylation status at individual CCCGGG sites across the genome was determined by counting the methylated reads with the CCGGG signature and unmethylated reads with the GGG signature at the beginnings of the sequencing reads after alignment to the human genome.
2015-09-08 | E-GEOD-66296 | biostudies-arrayexpress
Project description:This data was generated by ENCODE. If you have questions about the data, contact the submitting laboratory directly (Yijun Ruan mailto:ruanyj@gis.a-star.edu.sg). If you have questions about the Genome Browser track associated with this data, contact ENCODE (mailto:genome@soe.ucsc.edu). This track was produced as part of the ENCODE Project. It shows the locations of protein factor mediated chromatin interactions determined by Chromatin Interaction Analysis with Paired-End Tag (ChIA-PET) data (Fullwood et al., 2010) extracted from five different human cancer cell lines (K562 (chronic myeloid leukemia), HCT116 (colorectal cancer), HeLa-S3 (cervical cancer), MCF-7 (breast cancer), and NB4 (promyelocytic)). A chromatin interaction is defined as the association of two regions of the genome that are far apart in terms of genomic distance, but are spatially proximate to each other in the 3-dimensional cellular nucleus. Additionally, ChIA-PET experiments generate transcription factor binding sites. A binding site is defined as a region of the genome that is highly enriched by specific Chromatin ImmunoPrecipitation (ChIP) against a transcription factor, which indicates that the transcription factor binds specifically to this region. The protein factors displayed in the track include estrogen receptor alpha (ERa), RNA polymerase II (RNAPII), and CCCTC binding factor (CTCF). For data usage terms and conditions, please refer to http://www.genome.gov/27528022 and http://www.genome.gov/Pages/Research/ENCODE/ENCODEDataReleasePolicyFinal2008.pdf Chromatin interaction analysis with paired-end tag sequencing (ChIA-PET) is a global de novo high-throughput method for characterizing the 3-dimensional structure of chromatin in the nucleus. In the ChIA-PET protocol, samples were cross-linked and fragmented, then subjected to chromatin immunoprecipitation. The DNA fragments that were brought together by the chromatin interactions were then proximity-ligated. During this proximity-ligation step, the half-linkers (created by the fragmentation) containing flanking MmeI sites (type IIS restriction enzymes) were first ligated to the DNA fragments and then ligated to each other to form full linkers. Full linkers bridge either two ends of a self-circularized fragment, or two ends of two different chromatin fragments. The material was then reverse cross-linked, purified and digested with MmeI. MmeI cuts 20 base pairs away from its recognition site. Tag-linker-tag (paired-end tag, PET) constructs were sequenced by ultra-high-throughput methods (Illumina or SOLiD paired-end sequencing). ChIA-PET reads were processed with the ChIA-PET Tool (Li et al., 2010) by the following steps: linker filtering, short reads mapping, PET classification, binding site identification, and interaction cluster identification. The high-confidence binding sites and chromatin interaction clusters were reported. Chromatin interactions identified by ChIA-PET have been validated by 3C, ChIP-3C, 4C and DNA-FISH (Fullwood et al., 2009).
2012-07-19 | E-GEOD-39495 | biostudies-arrayexpress
Project description:ChIP-seq has become the method of choice for studying functional DNA-protein interactions on a genome wide scale. The method is based on co-immunoprecipitation of DNA binding proteins with formaldehyde cross-linked DNA, followed by deep sequencing of immunoprecipitated chromatin fragments, allowing for the identification of binding sites with high accuracy [1-5]. Traditionally, genome-wide tiling path microarrays were used for the detection of specifically immunoprecipitated DNA fragments, but current next-generation sequencers now generate sufficient data to assay multiple samples in a single sequencing run, making this method more time and cost effective compared to microarray-based approaches [2]. However, due to limitations in read length of next generation sequencing technologies (50-76bp), it is not possible to sequence the complete length of immunoprecipitated DNA fragments which can be up to 2kb long reflecting biology in case of big protein complexes. As a consequence only the ends of the immunoprecipitated DNA fragments are sequenced. Deconvolution of sequencing reads mapping to the positive and negative strand is required to identify the real DNA binding site [4-9]. This limitation is most obvious in case of large complex regions where multiple binding sites of various regulatory elements are clustered close to each other. Deconvolution of such regions and the exact identification of individual binding positions is very challenging [4]. Similar problems can occur when studying histone positioning in cases where the length of ChIP fragments is bigger than the average distance of 2 neighboring nucleosomes. In addition frequently only narrow size range of immunoprecipitated fragments is selected for sequencing and thus possible bias towards binding regions within the selected range can be expected by loosing larger fragments possible originated from protein complexes interacting with DNA. To circumvent these complications, we modified the procedure for the preparation of ChIP-seq samples by introducing an additional extensive fragmentation round after isolation of immunoprecipitated chromatin to generate 70-110 bp long DNA fragments. For testing, we choose Tcf4 protein which is a well-studied downstream element of the Wnt-pathway for which the genome wide binding site profile is known was already known from ChIP-on-CHIP experiments [10]. Using the modified approach, we were able to identify Tcf4 binding sites at near nucleotide resolution without computational deconvolution. Furthermore, high-resolution information on genome-wide binding site regions allowed for the identification of potential novel Tcf4 co-factors as well as target genes.
2010-11-01 | GSE18481 | GEO
Project description:ATACseq analyses between iEG (ETO2-GLIS2) and CTRL fetal primary cells and myeloid or CD41 iEG cells lines. Briefly, cells were isolated from iEG mice and cultured for 24 hours in RPMI supplemented with 10% FBS, cytokines (mIL3, mIL6, mSCF, mTPO, mFLT3l) and 100 ng/ml Doxycyclin. After cell lysis, transposition and purification step, the transposed DNA fragments were amplified by polymerase chain reaction (PCR) between 12 and 18 depending on the number of cells at the beginning (50,000 to 6,000) using adapters from the Nextera index kit (illumina). PCR purification was performed using Agencourt AMPure XP magnetic beads (Beckman Coulter A63880) in order to remove large fragments and remaining primers. Library quality was assessed using an Agilent 2100 Bioanalyzer using a High Sensitivity DNA chip (Agilent Technologies 5067-4626). Libraries were sequenced using Novaseq-6000 sequencer (Illumina) (50bp paired-end reads). Quality control of reads was performed using FastQC 0.11.7 and multiQC 1.5. The reads were aligned to the reference genome mm10 with bwa (aln 0.7.17). After alignment, we removed reads mapping to the mitochondrial genome, PCR duplicate reads and reads with a mapping quality lower than 20 using samtools (v 1.9). Final read counts for all mouse datasets ranged from 37 to 128 million reads. Mapped reads were normalized to bins per million (BPM) and were converted to bigwig format using deeptools (v3.2.0). Peak calling, differential analysis, annotation and motif analysis was performed using macs2 (V 2.1.2), Diffbind R package (v 2.8.0 in R-3.5.1 with threshold log2(1.5)), and homer (v4.10.4, annotatePeak.pl and findMotifsGenome.pl) respectively.
2019-10-15 | E-MTAB-8375 | biostudies-arrayexpress
Project description:Pseudomonas aeruginosa, a human opportunistic pathogen, is a common cause of nosocomial infections. Its ability to survive under different conditions relies on a complex regulatory network engaging transcriptional regulators controlling metabolic pathways and capabilities to efficiently use the available resources. P. aeruginosa PA3973 encodes a putative TetR family transcriptional regulator, with a helix-turn-helix motif involved in DNA binding. We applied transcriptome profiling (RNA-seq), and genome-wide identification of binding sites using ChIP-seq to unravel the biological role of PA3973.
2022-11-22 | GSE211771 | GEO
Project description:Our proteomics data (PXD030253) indicate RNA polymerase II is co-localized with Aldh18a1. To confirm this observation using an orthogonal method we performed ChIP-seq experiments to identify the binding patterns of RNA polymerase II in general, and RNA polymerase II with S2-phosphorylated CTD. We also performed Cut&Run using Aldh18a1 antibody which is submitted separately. To do a ChIP-seq, mES cells (46C) were treated with Accutase to be detached. The cells were fixed by 1.5% formaldehyde in PBS for 15min. Chromatin was fragmented by sonication and sheared to < 500 bp fragments. The cell debris were removed by spinning. The sheared chromatin was treated with 10 ug antibodies (pan-Pol II Bethyl A304-405A, or S2P-Pol II Bethyl A304-407A) overnight at 4 °C. Dynabeads protein A beads (40ul) were added to the samples and the samples were rotated 2-3 hours to perform the immunoprecipitation. The beads were washed 5 times. Chromatin was eluted in a denaturing solution, and proteins were digested using proteinase K. DNA was purified by SPRI beads. The ChIP libraries were prepared using NEBNext Ultra II DNA library prep according to the manufacturer manual. The libraries were sequenced using 150nt reads on paired-end mode by HiSeq4000. Reads were trimmed, aligned to the mouse genome (mm10) using Bowtie2, and duplicated reads were removed with samtools. The ChIP quality was evaluated by cross-correlation using the “SPP'' tool as suggested by ENCODE ChIP-seq guidelines. Peak calling was performed using MACS2.
2022-06-01 | E-MTAB-11510 | biostudies-arrayexpress