# Link ChIP-seq\`s peaks with nearby genes **URL:** \ **Category:** usegalaxy.eu support\ **Tags:** macs2, chip-seq, data-manipulation, epigenetics\ **Created:** [March 30, 2020, 9:15am UTC](https://help.galaxyproject.org/t/link-chip-seq-s-peaks-with-nearby-genes/3270 "2020-03-30T09:15:28Z")\ **Posts on this page:** 4\ **Page:** 1
**Author:** ![Manolis1](https://avatars.discourse-cdn.com/v4/letter/m/977dab/32.png) [@Manolis1](https://help.galaxyproject.org/u/Manolis1)\ **Post date:** [March 30, 2020, 9:15am UTC](https://help.galaxyproject.org/t/link-chip-seq-s-peaks-with-nearby-genes/3270/1 "2020-03-30T09:15:28Z")
Hello, I have followed successfully with my own ChIP-seq data the tutorial `From peaks to genes` . I haven´t got though what I was expecting to get at the end. I explain: At the end of the association of my peaks (from my ChIP-seq) with certain genes, I would ideally like to have a table including in seperate columns for each peak at least the info below: 1. Number of chromosome that the peak is located 2. Start of the peak 3. End of the peak 4. Position of the peak summit 5. The unique name of the peak (I have used MACS for the peak calling and the interval file that MACS gives assigns for each discovered peak a unique name in the 10th column (e.g.: MACS2.1.0\_in\_Galaxy\_peak\_1)) 6. ID/locus of the gene which is nearby (or associated) with the certain peak 7. Relative position of the peak to the TSS of the gene that is nearby. For example, is the peak 1.5 kb upstream from the TSS of the nearby gene, or perhaps 3 kb downstream from the TSS of the nearby gene? In short, where the peak is exactly located relative to the nearby gene? My question: is there any way to generate such a table in GALAXY? Many thanks in advance, Any help will be appreciated, Greetings, Manolis P.S.: I do get that if I apply the pipeline described in the tutorial `From peaks to genes` I can have at the end the overlap between peaks and related to them genes, but that’s it, nothing more. I cannot see how I could/or if I could at all, to get as a final output a table with all the columns mentioned above. ---
**Author:** ![jennaj](https://sea2.discourse-cdn.com/flex020/user_avatar/help.galaxyproject.org/jennaj/32/27_2.png) [@jennaj](https://help.galaxyproject.org/u/jennaj)\ **Post date:** [March 30, 2020, 8:01pm UTC](https://help.galaxyproject.org/t/link-chip-seq-s-peaks-with-nearby-genes/3270/2 "2020-03-30T20:01:20Z")
Hi @Manolis1 Yes. There isn’t a single tool that will create this type of table, but a combination of tools could be used, and a custom workflow created once you work out the processing. The tutorial you are following already covers several data manipulation tools as examples. The 101 also covers many data manipulation tools, including the different methods of joining data together based on common content. More tools that will be useful are in the tool groups “Bed” (bedtools), plus “Operate on Genomic Intervals” and “General Text Manipulation” (where most of the tools in the tutorial are grouped). Some experimentation on your part will be needed, not only with tool choices but with the reference annotation the peak/summit coordinates will be compared to. For example, to get the data for column 7, you will need the TSS coordinates for the gene to perform a distance calculation, and decide if you want to base that off the peak region (start or end) or the summit. In the final step, remove/rearrange columns of data to produce the final custom summary report (tool: `Cut`). The report will not be in `interval` format anymore but `tabular` – since it contains multiple coordinate regions – but you could always create properly formatted (and attribute labeled) `bed` or `interval` datasets from it as needed. Be aware that `bed` format is stricter than `interval` format. Derived `bed` datasets may include (in the 4th “name” column) the geneID and others the unique peak name, but not both, unless you merge the values together (without introducing whitespace). You may need to merge those two values together anyway – to use it as a common key to join intermediate datasets together. FAQs: - [Format help for Tabular/BED/Interval Datasets](https://galaxyproject.org/support/tabular/) - [Common datatypes explained](https://galaxyproject.org/learn/datatypes/) Hope that helps! ---
**Author:** ![Manolis1](https://avatars.discourse-cdn.com/v4/letter/m/977dab/32.png) [@Manolis1](https://help.galaxyproject.org/u/Manolis1)\ **Post date:** [April 1, 2020, 1:33pm UTC](https://help.galaxyproject.org/t/link-chip-seq-s-peaks-with-nearby-genes/3270/3 "2020-04-01T13:33:13Z")
Hi @jennaj Thank you very much for your response, the useful input and your suggestions. I am currently working on it… Just one more question, If you would like to answer: From where (or how) could I obtain the coordinates for TSS (or the ATG) for each gene? Greetings, Manolis ---
**Author:** ![jennaj](https://sea2.discourse-cdn.com/flex020/user_avatar/help.galaxyproject.org/jennaj/32/27_2.png) [@jennaj](https://help.galaxyproject.org/u/jennaj)\ **Post date:** [April 1, 2020, 4:35pm UTC](https://help.galaxyproject.org/t/link-chip-seq-s-peaks-with-nearby-genes/3270/4 "2020-04-01T16:35:38Z")
> [@Manolis1](#): > > From where (or how) could I obtain the coordinates for TSS (or the ATG) for each gene? The UCSC Table Browser is a good choice. This is where the tutorial extracts other bed data used. The TSS is the start coordinate of a transcript’s coding region. Please feel free to publish and post back a link to your updated workflow (or workflows!), it may help others. 🙂