# Visualizing reads and enriched regions

We can visualize the reads and enriched regions using the Partek _Genome Viewer_.

* Select **ChIP-Seq** from the spreadsheet tree
* Select **Chromosome View** from the _Visualization_ section of the _ChIP-Seq_ workflow

The _Chromosome View_ tab will open. The default tracks are the transcript tracks, the sequence visualization tracks, and the cytoband track (Figure 1).

![](https://1384254481-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FJVEESmJAPppJ3ijFq5aR%2Fuploads%2Fgit-blob-5de24fa9a48b4f335fcfb3202e5e1963016ee836%2F2017-08-10%2009_46_09-Partek%20Genomics%20Suite%20-%201%20\(ChIP-Seq\).png?alt=media)

Figure 1. Viewing ChIP-Seq data in Chromosome View

We can add additional tracks to view the results of our analysis.

* Select **New Track** from the left-hand side of the _Tracks_ tab
* Select **Add tracks from a list of spreadsheets**
* Select **Next >**
* Select the **p-value\_filtered.txt** and **Motifs\_instances.txt** tracks on the _Tracks_ panel
* De-select **Aligned Reads**
* Select **Create** (Figure 2)

![](https://1384254481-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FJVEESmJAPppJ3ijFq5aR%2Fuploads%2Fgit-blob-d8e8be945417d0e75bc32777b73599b7438b0bf3%2F2017-07-31%2010_45_38-Track%20Wizard.png?alt=media)

Figure 2. Adding spreadsheet tracks to Genome Viewer

This will display the enriched regions found in the samples and the locations of the motif instances from the _de novo_ motif discovery (Figure 3).

![](https://1384254481-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FJVEESmJAPppJ3ijFq5aR%2Fuploads%2Fgit-blob-99edf39fce526b24b67df80849794a78b61156eb%2F2017-08-10%2009_48_11-Partek%20Genomics%20Suite%20-%201_p-value_filtered_motifs_instances%20\(Motifs_instances.tx.png?alt=media)

Figure 3. Adding p-value and motif binding tracks

The two track display the detected regions at each location on the chromosome for the NRSF-enriched sample, _chip_, and aligned them to the _de novo_ discovered motif binding sites. If you have not gone through the steps for peak detection and motif discovery, these tracks will not be available.

To facilitate analysis, we will consolidate the tracks. To select a track, either left-click it in the _Tracks_ panel or activate Selection Mode by selecting the (![](https://1384254481-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FJVEESmJAPppJ3ijFq5aR%2Fuploads%2Fgit-blob-cd50f4a93daed4f289d4900fd0efb52dc26bfea0%2Fimage2017-7-31%2011_57_48.png?alt=media)) icon and left-click the track in the viewer. Selecting a track opens its configuration options in the track properties panel below the _Tracks_ panel.

* Select **RefSeq Transcripts - 2016-10-18 (hg18) (+)** on the _Tracks_ panel
* Set _Strand_ to **Both** using the drop-down menu in the Track section of the track properties panel
* Select **Apply**

The _RefSeq Transcripts_ track now shows transcripts from both strands. We can remove the dedicated (-) strand track.

* Select **RefSeq Transcripts - 2016-10-18 (hg18) (-)** on the _Tracks_ panel
* Select **Remove Track**

Several other tracks will not be used in our analysis; we can remove them as well.

* Select **Legend: Base Colors**, **Genome Sequences (hg18)**, and **Cytoband (hg18)** on the _Tracks_ panel by left clicking on each while holding the **Ctrl** key on your keyboard
* Select **Remove Track**

The height of many tracks can be adjusted by changing the _Track height_ setting in the _Profile_ tab.

* Select **Regions (1/p-value\_filtered (p-value filtered.txt))** on the _Tracks_ panel
* Set _Track height_ to **25** using the slider
* Select **Apply**
* Repeat for **Regions (1/p-value\_filtered/motifs/instances (Motifs\_instances.txt))**

This view clearly shows the _Bam Profile_ tracks for _chip_ and _mock_ samples alongside smaller track indicating motif binding sites, significant peaks, and RefSeq transcripts (Figure 4).

![](https://1384254481-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FJVEESmJAPppJ3ijFq5aR%2Fuploads%2Fgit-blob-8f61e1ae5b408c1c16911a4f390971a911fbdb72%2F2017-07-31%2011_19_08-Partek%20Genomics%20Suite%20-%201%20\(ChIP-Seq\).png?alt=media)

Figure 4. Genome Viewer tracks can be modified to facilitate data exploration and analysis

We can also see strand-specific information by color-coding forward and reverse strand read information on the _Bam Profile (chip)_ track.

* Select **Bam Profile (chip)** on the _Tracks_ panel
* Select **Alignments** in the _Display Color Options_ section of the _Style_ tab
* Select Strands under _Color by_
* Select **Apply**

A legend track will be added at the bottom of the viewer showing that forward strand reads are colored green and reverse strand reads are colored red. We can move this track to below _Bam Profile (chip)_.

* Left-click **Legend: strands** in the _Tracks_ panel and drag it up the list to below _Bam Profile (chip)_ (Figure 5)

![](https://1384254481-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FJVEESmJAPppJ3ijFq5aR%2Fuploads%2Fgit-blob-a8125d1123afbca0c2d31b74e8f35609bfe8adb2%2F2017-07-31%2011_37_58-Partek%20Genomics%20Suite%20-%201%20\(ChIP-Seq\).png?alt=media)

Figure 5. Move tracks by left-clicking and dragging them in the Tracks panel

_Chromosome View_ opens to a whole-chromosome view of chromosome 1. To analyze the data we can zoom in on the data (Figure 6).

![](https://1384254481-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FJVEESmJAPppJ3ijFq5aR%2Fuploads%2Fgit-blob-1259fb56326dbfa7daa58263413f7e7d0094e364%2F2017-08-10%2010_01_24-Partek%20Genomics%20Suite%20-%201%20\(ChIP-Seq\).png?alt=media)

Figure 6. Viewing the zoomed-in view of an enriched region showing two possible binding sites at the location

There are four methods for zooming in and out in _Chromosome View._ (1) Select (![](https://1384254481-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FJVEESmJAPppJ3ijFq5aR%2Fuploads%2Fgit-blob-2fa4103bc9ebbaea62a5b3bbb91f4164df2cf189%2Fnavigation%20mode%20icon-1.png?alt=media)) to activate Zoom/Navigate Mode. With this mode selected, left-clicking and drawing a box on a track will zoom to the region in the box. (2) Use the mouse scroll wheel (3) Select the zoom in (![](https://1384254481-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FJVEESmJAPppJ3ijFq5aR%2Fuploads%2Fgit-blob-e0464530c5ae59682616ce7d311a4ce349b95a47%2Fimage2017-7-31%2011_50_19.png?alt=media)) and zoom out (![](https://1384254481-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FJVEESmJAPppJ3ijFq5aR%2Fuploads%2Fgit-blob-d53ab4e1c846639fc6fcc4a207fdda1e279921d8%2Fimage2017-7-31%2011_50_55.png?alt=media)) on the zoom controls (![](https://1384254481-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FJVEESmJAPppJ3ijFq5aR%2Fuploads%2Fgit-blob-539fff7602b613a2560354633846670d03ea0b56%2F2017-07-31%2011_48_23-Partek%20Genomics%20Suite%20-%201%20\(ChIP-Seq\).png?alt=media)). Selecting (![](https://1384254481-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FJVEESmJAPppJ3ijFq5aR%2Fuploads%2Fgit-blob-4e80692f66acff9f236ad06c0cadd8624bf19f8b%2Fimage2017-7-31%2011_54_53.png?alt=media)) resets the zoom to whole chromosome. (4) Use the slider on the zoom controls (![](https://1384254481-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FJVEESmJAPppJ3ijFq5aR%2Fuploads%2Fgit-blob-539fff7602b613a2560354633846670d03ea0b56%2F2017-07-31%2011_48_23-Partek%20Genomics%20Suite%20-%201%20\(ChIP-Seq\)-1.png?alt=media)). (5) Zoom to a location or gene by typing the coordinates or gene name in the Navigation bar (![](https://1384254481-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FJVEESmJAPppJ3ijFq5aR%2Fuploads%2Fgit-blob-19db666bb456b4313311825555bf6f5f84add0a9%2F2017-07-31%2012_05_47-Partek%20Genomics%20Suite%20-%201%20\(ChIP-Seq\).png?alt=media)) and pressing Enter on your keyboard. (5) Starting on a spreadsheet, right-click a row with region information, and select **Browse to Location**.

We can practice this using a gene highlighted in the paper describing this data set, NEUROD1.

* Type **NEUROD1** in the navigation bar
* Select **Enter** on your keyboard to navigate to NEUROD1 (Figure 7)

![](https://1384254481-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FJVEESmJAPppJ3ijFq5aR%2Fuploads%2Fgit-blob-173906ec0136c3dd72527d56d9d24279859f9a27%2F2017-08-10%2010_09_50-Partek%20Genomics%20Suite%20-%201%20\(ChIP-Seq\).png?alt=media)

Figure 7. Zooming to the NEUROD1 gene

NEUROD1 contains a binding site for the NRSF motif. Notice that the enriched region for the NRSF transcription factor is within the NEUROD1 gene. As discussed in the Johnson et al. paper, NRSF is implicated in the repression of NEUROD1, but it was unknown exactly where the NRSF binding occurred. This data indicates that the binding site is within the NEUROD1 gene itself, as shown by the orange box in the Regions track. We can zoom in further to view the forward and reverse strand read histogram.

* Select (![](https://1384254481-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FJVEESmJAPppJ3ijFq5aR%2Fuploads%2Fgit-blob-2fa4103bc9ebbaea62a5b3bbb91f4164df2cf189%2Fnavigation%20mode%20icon-2.png?alt=media))
* Left-click and drag to draw a box around the predicted binding site (Figure 8)

![](https://1384254481-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FJVEESmJAPppJ3ijFq5aR%2Fuploads%2Fgit-blob-9f0cee01b3cd90bfc1e212806a02b681c9e026d6%2F2017-08-10%2010_07_04-Partek%20Genomics%20Suite%20-%201%20\(ChIP-Seq\).png?alt=media)

Figure 8. Using Zoom/Navigate Mode. Drawing a box around a region on any track will zoom all tracks to that region

Zoomed-in further, we can see the intersecting peaks from forward and reverse strand reads (Figure 9).

![](https://1384254481-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FJVEESmJAPppJ3ijFq5aR%2Fuploads%2Fgit-blob-95f3ab0a1f6bc4eb3f071a2eb7a9f3a0ae3e55a4%2F2017-08-10%2010_08_48-Partek%20Genomics%20Suite%20-%201%20\(ChIP-Seq\).png?alt=media)

Figure 9. Zoomed-in view of NEUROD1 binding site for NRSF

You can save the reads shown in the visible genome browser window by selecting (![](https://1384254481-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FJVEESmJAPppJ3ijFq5aR%2Fuploads%2Fgit-blob-cd50f4a93daed4f289d4900fd0efb52dc26bfea0%2Fimage2017-7-31%2011_57_48-1.png?alt=media)), right-clicking in the peak area on the _Bam Profile_ track, and selecting **Dump Displayed Reads to Spreadsheet** from the pop-up menu.

## Additional Assistance

If you need additional assistance, please visit [our support page](http://www.partek.com/support) to submit a help ticket or find phone numbers for regional support.
