ChIP-Seq+Data+Analysis+using+Galaxy+Spring+2015

= = flat =**Getting a Galaxy Account **= Galaxy is a web interface for running a lot of command line tools used for NGS analysis. It is great for beginners and helps your streamline your data analysis. You can create pipelines (workflows) and share with lab mates. You can also download other open source published workflows to do your own analysis. There is a free public galaxy at usegalaxy.org but today we setup our own private galaxy in the Amazon cloud to work through the RNA-Seq exercise. To access this galaxy first go to http://ec2-52-1-38-85.compute-1.amazonaws.com and signup for a temporary galaxy account. Note: This account is valid only for today's workshop and it is different for different galaxy instances you log in to.

=ChIP-Seq: FASTQ data and quality control =

We will ultimately be comparing the binding profiles of two transcription factors, Nanog and Pou5f1 (Oct4) from the [|HAIB TFBS ENCODE collection] that have been ChIP-sequenced with Illumina in h1-ESC cells. Do look up the transcription factor (TF) information if you are not familiar with them as it should give you an idea what to expect in terms of TF binding.

To keep things manageable and allow algorithms to finish within a few minutes we will work with a single sample (//"h1-hESC - Input replicate 1" //), using reads from 32.8 Mb of chromosome 12 (chr12:1,000,000-33,800,000). We will switch over to using all samples to compare biological replicates and look at differential binding by different transcription factors by retrieving them from the Shared Data folders as needed.

Exploring the FASTQ files
Open up Galaxy from your machine as before, start up a new blank history and give it a name. In the top menu, move your mouse cursor over Shared data and select Data libraries. This allows you to import files into Galaxy that have been shared with you.

Browse to the ChIP-Seq Example folder, expand the Unfiltered reads folder, and click on the arrow to the right of the filename H1hesc_Input_Rep1_chr12_unfiltered.fastq. In the context menu select Import this dataset into selected histories. On the next screen, select the history you are using. The file should now appear in your history

This is the kind of data you would expect to receive from your sequencing core facility: a large number of short nucleotide reads in [|FASTQ format] (PDF).


 * __ ** Take a look at the file contents using Galaxy's preview and file view functionality. What additional information does FASTQ contain over FASTA? How many reads are in your file? Pick the first sequence in your FASTQ file and note the identifier. What is the quality score of it's first and last nucleotide, respectively? ** __

<span style="color: #333333; font-family: 'Proxima Nova Bold','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 28px;">Quality Controls
<span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">The FASTQ file contains output reads from the sequencer that need to be mapped to a reference genome for us to understand where those reads came from on the sequenced genome. However, before we can delve into read mapping, we first need to make sure that our preliminary data is of sufficiently high quality. This involves several steps:


 * 1) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Obtaining summary quality statistics on the reads and review diagnostic graphs
 * 2) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Eliminate sequencing artifacts
 * 3) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Filter out genetic contaminants (primers, vectors, adaptors)
 * 4) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Filter out low-quality reads
 * 5) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Recalculate quality statistics and review diagnostic plots on filtered data

<span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">Iterate through steps 2-5 until the data is of sufficient quality before proceeding to mapping.

<span style="color: #333333; font-family: 'Proxima Nova Bold','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 21px;">Obtain Quality Statistics
<span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">While Galaxy has a number of different methods to assess sequence properties you will want to stick to using a standard QC tool such as [|FASTQC] in most cases.

<span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">Under the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">NGS: QC and manipulation tool find <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">FASTQC: Read QC and give it a try.


 * __ ** What does the quality plot tell you about your data? What additional information do you get with this tool? What can you say about the per-base GC content? Does this match the human genome? Why is there a shift in the GC distribution from the theoretical distribution? Finally, what is the meaning of the k-mer content tab? This publication from [|Schroeder et al] (PDF) has some background information ** __.

<span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">Compare your results with a [|pre-generated FastQC plot] obtained from the full FASTQ data set (i.e., using all reads).


 * __** Do you note any differences in the output or format? Why is that? **__

<span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">Some of the diagnostic plots from a [|Core conference call] might be useful if you are curious to see what different modes of failure frequently look like.

=<span style="color: #333333; font-family: 'Proxima Nova Bold','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 28px; text-decoration: none;">ChIP-Seq: Trimming and filtering =

<span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">Depending on your FASTQ results you may need to test for contamination (either by vector or other sourcers) and/or remove adapter sequences from your read prior to alignment. It is good to remove these sequences to increase your mapping efficiency. Note that this is not a problem for our sample data.

<span style="color: #333333; font-family: 'Proxima Nova Bold','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 28px;">Screen for contamination
<span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">Not all required tools are available through Galaxy just yet. For example, depending on the source of your data you probably want to screen your reads for contamination using tools such as [|ContEst] (PDF) or [|SeqTrim] (PDF). For now we will just explore how you can spot potential contamination from the FastQC results alone. > What is the most overrepresented sequence in the pre-generated FastQC results from the h1-hESC sample? <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">Compare your own FastQC results with a different [|pre-generated FastQC summary] that highlights a potential problem. > Study the report. Is there any vector contamination present? Is there any adapter contamination present? <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">As an emergency solution to remove contaminating sequences, Galaxy comes with a way to <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Clip adapter sequences. If you notice a particularly strong contamination with a vector or primer you can use your FASTQ file as the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Library to clip, keep the default <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Minimum sequence length , and enter the sequence you spotted in the FASTQC report as the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Adapter to clip. You probably want to keep both clipped and unclipped sequences in this case; keeping only reads that contained an adapter or barcode can be useful if you expect this to be a feature of all valid reads. For command-line use tools such as [|cutadapt] are more versatile and less labor-intensive to use. <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">Normally, you would take the filtered output and feed it back into Galaxy, but our dataset has very little contamination so we will move on with the original data set.

<span style="color: #333333; font-family: 'Proxima Nova Bold','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 28px;">Filter by quality
<span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">As we do not plan on using our data to identify genomic variation individual sequencing errors become less of a problem as long as the sequence can be still mapped to the reference genome. Still, you want to filter sequences where the overall quality is just too low. <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">After reviewing the quality diagnostics from your FASTQC report, choose an appropriate lower nucleotide quality that you think would be good for this data (for instance, a quality score below 20 means there is more than a 1/100 chance of an incorrect base call). Use this cutoff to filter out low-quality reads from your data with the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Filter By Quality tool. After you have filtered your data generate another <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">FASTQC quality report.
 * __** Do you have an improvement in read quality? What has changed and why? What percentage of reads did you retain at your cutoff? **__

<span style="color: #333333; font-family: 'Proxima Nova Bold','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 28px;">Read Trimming
<span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">After the general filtering which removes sequence reads of overall low quality it can be a good idea to trim the end of your reads, cutting away nucleotides at the end that are still of low-quality
 * __** Is this necessary in your case? **__

<span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">Use the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Trim sequences tool to trim if required. Make sure to name your finalized FASTQ file that has been filtered and improved. You will be using this file in the next section for mapping. Again, tools such as [|sickle] provide more flexibility by allowing for adaptive, window-based trimming, but are not yet part of the default Galaxy toolkit. <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">It might also be helpful to create a new history at this time and copy over your final sequence file. An easy way to do this is through the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">History menu and the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Copy datasets function which allows you to pick one or more items in your current history, copying them over to a new or existing one.

=<span style="color: #333333; font-family: 'Proxima Nova Bold','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 28px; text-decoration: none;">ChIP-Seq: Alignment and visualization =

<span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">If for some reason any of the previous steps did not work you can retrieve a quality-controlled, trimmed and filtered FASTQ file from the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Data Library in the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">ChIP-Seq Data library, subfolder <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Sequence and reference data (<span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">H1hesc_Input_Rep1_chr12_qualityfiltered.fastq ).

<span style="color: #333333; font-family: 'Proxima Nova Bold','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 28px;">Standard alignment
<span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">Next we are going to look at the steps we need to take once we have a clean, filtered FASTQ file that is ready for alignment. Use the filtered FASTQ file that you prepared in the previous section. The alignment process consists of the following steps:
 * <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Choose an appropriate reference genome to map your reads against
 * <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Choose an appropriate gene annotation model to guide your alignment
 * <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Perform the read alignment

<span style="color: #333333; font-family: 'Proxima Nova Bold','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 21px;">Load the filtered FASTQ file and Reference Genome
<span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">To avoid excessive runtimes during later steps we will not align the reads against the whole human genome, just <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">chr12:1,000,000-33,800,000. We have retrieved the sequence for this chromosome from UCSC GoldenPath (hg19) and added it to the "Data Library". This means we won't use the pre-built genome indices, but submit our own reference sequence. We are also using only a subset of reads from a single sample for the next few steps. <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">Using it generates a slight time overhead since the reference genome needs to be indexed prior to each alignment run, but this is acceptable for now.
 * 1) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Start up a new blank history and give it a name.
 * 2) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Either load your filtered FASTQ files from the previous tutorial, or download a filtered FASTQ file from the Shared Data Libraries in Galaxy.
 * 3) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Find the chromosome 12 genomic sequence in the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">ChIP-Seq Example folder as <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">hg19-chr12-reference_sequence.fa and import it into your new history.

<span style="color: #333333; font-family: 'Proxima Nova Bold','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 21px;">Perform Read Alignment
<span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">You will start by aligning your reads to chr12 using a standard aligner, <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Bowtie. If you'd like to learn more about its underlying algorithm, advantages and drawbacks a recent [|review from Bao et al] (PDF) might be of interest. It also references the original paper which walks you through the methods. <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">You can find the tool under the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">NGS:Mapping tool heading. Make sure you select the corresponding tool for <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Illumina. Do not use the built-in reference genomes, and instead, select <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Use one from the history and choose the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">chr12-reference_sequence.fa. Select <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Single-End under “library mate pairing”. In the Bowtie settings, choose <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Full parameter list. As you can see, there are many parameters that you can change depending on the type of data you are working with. Change the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Suppress all alignments for a read if more than n reportable alignments exist (-m) from <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">-1 to <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">1 to ensure that only uniquely mapping reads are reported. Click the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Execute button to launch the alignment. [//<span style="font-family: 'Proxima Nova Italic','Helvetica Neue',Arial,Helvetica,sans-serif;">Note: There's also a Bowtie Version 2 by now, but we'll stick to Bowtie 1 for the time being. //]
 * __** While the alignment is running, do [|look up the documentation] to see what the various tool options do. Take a look at the output file. Note it’s size. How long did it take to run? Now extrapolate to how long you would expect this tool to run when mapping to the entire genome (approximately). **__

<span style="color: #333333; font-family: 'Proxima Nova Bold','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 21px;">Filter unmapped reads and convert SAM files to BAM files
<span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">Most aligners produce output in [|"SAM" format] (Sequence Alignment/Map format by now, see [|publication from Heng Li] (PDF) for more details). > Explore the SAM format. What key information is contained? <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">The SAM formatted output contains all reads along with flags indicating whether they have been mapped or not, their quality values and the genomic coordinates for mapped reads. Filter out the unmapped reads by applying the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Filter SAM tool under <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">NGS:SAM Tools. Click on the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Add a new flag button. Select <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">The read is unmapped under <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Type and set the state for this flag to <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">No to discard unmapped reads. > How many lines are in this final file? This represents the number of aligned reads. Calculate the percentage of aligned reads for this experiment. <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">Another common preprocessing step for NGS data is to filter out duplicates. For ChIP-Seq data, these can be the result of a genuine signal or can be caused by PCR amplification artifacts, especially at high sequencing depth. Duplicates that result from ChIP enrichment form nice peaks surrounded by an otherwise more even distribution of reads. In general, the recommendation is to avoid directly filtering duplicates. MACS handles biases introduced by amplification and library preparation by removing duplicate reads if their count is higher than can be expected from the sequencing depth (binomial distribution with a p-value < 10e-05). <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">To speed up processing those need to get converted into their binary counterpart (BAM).
 * <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Apply the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">SAM-to-BAM tool under the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">NGS:SAM Tools heading to convert your output SAM file to the compressed BAM format, once again using the chromosome 12 reference sequence from your history as reference

<span style="color: #333333; font-family: 'Proxima Nova Bold','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 21px;">Explore alignment in a genome browser
<span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">Take a look at your alignment in IGV (Integrative Genomics Viewer). To do so, first expand the alignment result (the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">BAM file) in your history and click on the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">web current link beside <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">display with IGV. This will download IGV to your machine. It should start automatically after a few security warnings (if not, find the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">igv.jlnp file and click it yourself). <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">You should see a window that looks like this: <span style="color: #333333; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px; text-decoration: none;"> <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">Click in the highlighted window… <span style="color: #333333; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px; text-decoration: none;"> <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">… and enter the gene name <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">AKAP3. You should now be able to see the individual reads, you can zoom in further using the controls on the upper right and scroll around by clicking and dragging in the highlighted alignment track. Note that IGV automatically translated the gene symbol into the matching genomic coordinates for your genome build: <span style="color: #333333; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px; text-decoration: none;"> <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">Recall that this is the input file, i.e. the control. <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">ChIP-seq data is not distributed evenly across the genome and, even in the input sample, you will see clustering of reads around gene regions and especially around promoter areas, which is why we need the input/control sample. By comparing the number and distribution of reads between input and your sample, you essentially correct for this background read distribution. <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">Let's compare the input and the corresponding Nanog TF ChIP-Seq data side by side, and this will be immediately apparent. <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">Launch IGV on your local computer. Expand the alignment result (the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">BAM file) in your history by clicking on the file name and then click on the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">local link beside <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">display with IGV. This will automatically download the alignment file into the IGV window that is already open so you can compare the two alignment files.
 * __ ** Can you identify any potential regions of enrichment in this gene or others? ** __
 * <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Go to <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Shared Data, <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Data Libraries , open the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">ChIP-Seq Example folder, then open the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Mapped folder. Find the Nanog alignment data <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">H1hesc_Nanog_Rep1_chr12_aln.bam and import it into your current history.
 * __ ** Can you identify any potential Nanog binding sites based on the comparison of the input vs. Nanog ChIP results? ** __

=<span style="color: #333333; font-family: 'Proxima Nova Bold','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 28px; text-decoration: none;">ChIP-Seq: Peak calling =

<span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">As before, if any of the previous steps did not work out you can retrieve intermediate results from the "Snapshot" <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Data Library under "Shared Data". We have deposited the Bowtie-aligned reads in BAM format in the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Mapped library. <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">The goal of peak calling is to define regions with a high density of reads indicating where the studied factor bound. There are [|multiple algorithms to perform the peak-calling] step. In this workshop we will be using a tool called MACS. <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">ChIP signal strength is a continuum with far more weak sites than strong. As a result, the composition of the final peak list depends heavily on the specific parameter settings and algorithm as well as the quality of the experiment itself. Thresholds that are too relaxed lead to a high proportion of false positives for each replicate, but as you will see in the Handling Replicates section, subsequent analysis can remove false positives from a final joint peak determination.

<span style="color: #333333; font-family: 'Proxima Nova Bold','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 28px;">Peak calling using MACS
<span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">We will use MACS to call "Nanog" peaks in a ChIP (tag) sample versus a control (input) sample. Although MACS can be run without a control data set, an appropriate control is critical for robust analysis of any ChIP-seq experiment because DNA breakage during sonication is not uniform. Some regions of the genome, particularly those at open chromatin, are sheared more readily than others and these fragments can be over-represented in the sequencing step, giving rise to false positive sites of occupancy. The control data set for this study used DNA isolated from cells that were cross-linked and fragmented under the same conditions as the ChIP DNA but without the antibody. This is sometimes referred to as “reverse-cross-linked chromatin” or “Input DNA”. An alternative approach is a “mock IP” in which a control antibody IgG antibody is used that does not recognize any DNA binding protein (“IgG” control). Additional information can be obtained from the [|ENCODE/modENCODE ChIP-seq guidelines]. <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">While MACS is running, use the time to browse through its [|manual] or explore the [|MACS Protocol paper] (PDF) which describes how to set parameters for transcription factor binding and histone modification data sets, and how to interpret the output. A new version, [|MACS2], has been released but is not available in this Galaxy instance. <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">The algorithm will generate two files. The first output is a BED file with the coordinates of the peaks. Instructions for viewing the peaks in IGV follows below. <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">The second output is an HTML page that links to the results in various formats and lists the parameters and tasks that were performed. Although warning messages do not affect the success of a MACS run, the majority should still be carefully evaluated. For example, the warning message //<span style="font-family: 'Proxima Nova Italic','Helvetica Neue',Arial,Helvetica,sans-serif;">'Treatment tags and Control tags are uneven' //means that the FDR of the resulting peaks will be overestimated when the control sample has more reads and will be underestimated when the ChIP-seq sample is sequenced more deeply. The message //<span style="font-family: 'Proxima Nova Italic','Helvetica Neue',Arial,Helvetica,sans-serif;">'Fewer paired peaks X than 1,000' // means that MACS only identified X model peaks and may indicate potential data quality issues because 1,000 model peaks are needed to robustly estimate ChIP-DNA fragment size. The message //<span style="font-family: 'Proxima Nova Italic','Helvetica Neue',Arial,Helvetica,sans-serif;">'missing chromosome X data' // might suggest that the raw input file for that chromosome is incomplete. <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">Under <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Additional Files, click on the first file ending with <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">model.pdf. This file shows the peak model built by MACS. ChIP-Seq DNA fragments are sequenced from the 5’ end, resulting in alignment of tags that form two peaks, on each strand, flanking the location that was bound by the protein of interest. The red curve represents the distribution of reads on the positive strand at each base position, and the blue line models the reads on the negative strand. The black curve is the combined profile that is generated by shifting each distribution toward the center. The detected DNA fragment length is given by d and varies among different ChIP-seq libraries. >> **Note that the estimated fragment length (d) is 103 which is much smaller than the default band width parameter of 300 used by MACS. This accounts for the jagged peak model. If you run it again and modify the band width parameter to 100, you will see smooth peaks in the peak model diagram.** <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">The <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">peaks.xls file contains all of the information about the peaks identified by MACS, including key parameters, coordinates and enrichment statistics.
 * 1) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Copy the aligned control (input) BAM file and the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">H1hesc_Nanog_Rep1_chr12_aln.bam from the previous section to a new history.
 * 2) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Select <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">MACS from the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">NGS:Peak Calling menu.
 * 3) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Enter an <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Experiment Name (e.g., //<span style="font-family: 'Proxima Nova Italic','Helvetica Neue',Arial,Helvetica,sans-serif;">MACS_NanogRep1 //)
 * 4) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Select the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">H1hesc_NanogRep1_chr12_aln.bam file for the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">ChIP-Seq Tag File, and the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">H1hesc_InputRep1_chr12_aln.bam file for your <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">ChIP-Seq Control File.
 * 5) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Note the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Effective genome size, which is the size of the genome that is can be sequenced and mapped (minus telomeres, repeats, etc). Use the default since we are dealing with human data (2,700,000,000). For mouse, the value is 1,870,000,000. For less commonly used species, an estimate of 70-90% of the full genome size is reasonable.
 * 6) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Set the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Tag size to 36bp (based on examining the FastQC report).
 * 7) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Note the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Band width parameter. This value is used to scan the genome for model building. You can set this parameter as the expected sonication fragment size. We will leave it at the default value of 300.
 * 8) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Note the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">MFOLD parameter which is used to select genomic regions to build the experiment-specific peak model. A region is used for the model if it contains more than mfold reads than the region surrounding it. We will use the default value (32) but this value can be tweaked to obtain more or less peaks.
 * 9) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Check the box for <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Perform the new peak detection method (futurefdr) . This uses both the control and treatment data to calculate local bias and identify peak locations.
 * __ ** How many peaks (“regions”) were detected by MACS? ** __

<span style="color: #333333; font-family: 'Proxima Nova Bold','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 21px;">Visual exploration
<span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">We will use IGV again to study some of the results:
 * 1) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">View the resulting MACS <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">peaks: bed file in IGV. As Galaxy does not support automatic export of BED files to IGV, you will have to add it to IGV yourself. To do so, //<span style="font-family: 'Proxima Nova Italic','Helvetica Neue',Arial,Helvetica,sans-serif;">download // the file from Galaxy and take note of where you saved it. Now move over to IGV and click the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">File menu and then <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Load from file . Select the file you downloaded and open it in IGV to explore the MACS peak calls next to your Bowtie alignment.
 * 2) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Type <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Akap3 in the coordinate text box to zoom in on this region and observe the correlation between the MACS peak call and the alignment. Another interesting gene to look at is <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Sox5.

<span style="color: #333333; font-family: 'Proxima Nova Bold','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 28px;">Handling replicates
<span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">In your own analyses, you will have multiple conditions and (hopefully!) multiple replicates for those conditions. //<span style="font-family: 'Proxima Nova Italic','Helvetica Neue',Arial,Helvetica,sans-serif;">Technical // replicates are usually merged after FASTQC analysis to ensure that the sequencing quality is adequate. FastQ files can be concatenated at the read level or BAM files can be merged after alignment using SAMtools. <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">//<span style="font-family: 'Proxima Nova Italic','Helvetica Neue',Arial,Helvetica,sans-serif;">Biological // replicates should be analyzed individually first and then compared at the peak call level. Key questions to ask are whether replicates correlate well and if there is substantial peak overlap. The current state-of-the-art for analyzing ChIP-Seq replicates and setting thresholds is the [|Irreproducible Discovery Rate] (IDR), which compares a pair of ranked lists of peaks and creates a curve to quantitatively determine the extent to which the ranks are no longer consistent across the replicates. The method was developed and used extensively for the ENCODE and modENCODE projects and is part of their [|ChIP-seq guidelines]. The method is not yet available in the standard Galaxy installation, however. <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">In the meantime, there are simple things in Galaxy that you can do to analyze replicates and additional samples, such as intersecting and subtracting peaks to get common and unique regions in two samples.

<span style="color: #333333; font-family: 'Proxima Nova Bold','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 21px;">Bringing in replicates and another sample
<span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">Up to now we have been analyzing a single ChIP experiment (Nanog-Rep1-chr12 sample and Input-Rep1-chr12 control) to save time. Now we will bring three other ENCODE data sets into our analysis. We have generated the necessary files -- the second H1Hesc Nanog replicate and two replicates for H1Hesc Pou5f1 (Oct4) -- for you and placed them in the Shared <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Data Library under <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">ChIP-Seq Example ; they can be found in the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace;">MACS Outputs <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif;">folder. Open up this folder and import the following peak calls for each of the other 3 ENCODE samples to your history (making sure to specify the bed format): <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">To increase our confidence in the peak calls, we will combine our replicate files and take the intersection of the two. Click on <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Operate on Genomic Intervals and <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Concatenate the two Nanog replicates to combine them into one file. Next, use the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Cluster tool on the concatenated data to combine overlapping peak calls using the default parameters. <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">A simple heuristic used by the ENCODE project prior to IDR was to retain replicates for analysis if either 80% of the top 40% of identified targets in one replicate were among the targets in the second replicate; or 75% of target lists were in common between both replicates. <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">Note that the resulting peak calls are dominated by the //<span style="font-family: 'Proxima Nova Italic','Helvetica Neue',Arial,Helvetica,sans-serif;">“weakest” // replicate. You will probably have noticed that there is also an <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Intersect tool and may be wondering why we did not use it directly on the two replicates. If we had, only peaks in the first replicate that overlap with the second would be returned, i.e. overlapping regions would not be merged to form a larger region encompassing the coverage from both replicates.
 * <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">MACS_Nanog_Rep2.chr12.bed
 * <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">MACS_Pou5f1_Rep1.chr12.bed
 * <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">MACS_Pou5f1_Rep2.chr12.bed
 * __ ** How many peaks are retained? What fraction of the individual replicates are retained? ** __
 * __ ** Repeat these steps for the Pou5f1 replicates to prepare to look for differential enrichment of binding sites. ** __

<span style="color: #333333; font-family: 'Proxima Nova Bold','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 21px;">Differential binding
<span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">To identify sites that have differential binding, across two conditions, developmental stages, factors, etc., we will use Galaxy’s <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Subtract tool on the intersected peak calls for Nanog and Pou5f1. Select the clustered Pou5f1 file for the first input and the clustered Nanog file for the second input, and return <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">intervals with no overlap. The resulting dataset contains regions where Nanog binds but Pou5f1 does not. You can then switch the inputs to find regions where Pou5f1 binds but Nanog does not. <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">Note: You can also use MACS directly to find differential binding between two conditions by treating one of the samples as the control.


 * __ ** Find regions where both TFs bind using the approach above for intersecting two replicates. ** __

=<span style="color: #333333; font-family: 'Proxima Nova Bold','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 28px; text-decoration: none;">ChIP-Seq: Enrichment analysis =

<span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">Motif and functional enrichment analyses are unlikely to give reliable results given that you are only looking at 32.8 Mb of reads mapping to chr12. To obtain more signal, we will use the complete set of Nanog peak calls generated by ENCODE and look at binding sites for the //<span style="font-family: 'Proxima Nova Italic','Helvetica Neue',Arial,Helvetica,sans-serif;">whole // of chr12. This will also keep the sequence data size manageable for this tutorial.

<span style="color: #333333; font-family: 'Proxima Nova Bold','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 28px;">ChIP-Seq: Motif enrichment analysis
<span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">Start by importing <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Encode-hesc-Nanog.narrowPeak from the Galaxy “Data Library”. It is in [|//narrowPeak// format] and contains Nanog peak calls for the whole genome. The format is an extension of the BED format and was developed to describe peaks of signal enrichment for pooled, normalized data. <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">Convert this file to <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">BED format by editing the attributes and selecting the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Convert Format tab. Select <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Convert Genomic Intervals to BED and click <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Convert. <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">To perform motif enrichment, we need to retrieve the genomic DNA sequences for the set of peaks. To do this, you will also need the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">hg19-chr12-reference_sequence.fa we used previously for the alignments. Import this file into your current history from either a previous history or the Data Library. Recall that it was located in Sequence and reference data in the Shared Data library. <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">Note that only the chr12 regions are generated because we used the chr12 reference sequence. Download this file to your computer for use in the motif over-representation analysis. <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">To identify over-represented motifs, we will use [|DREME] from the MEME suite of sequence analysis tools. DREME is a motif discovery algorithm designed to find short, core DNA-binding motifs of eukaryotic transcription factors and is optimized to handle large ChIP-seq data sets. <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">You will be shown a status page describing the inputs and the selected parameters, as well as a link to the results at the top of the screen. Clicking the link displays an output page that is continually updated until the analysis has completed. DREME will also email you once the analysis is complete. This may take some time depending on the server load and the size of the file. While you wait, take a look at the expected result (obtained by clicking on the HTML output link once the analysis has completed). <span style="color: #333333; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px; text-decoration: none;"> <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">DREME’s HTML output provides a list of <span style="font-family: 'Proxima Nova Semi-Bold','Helvetica Neue',Arial,Helvetica,sans-serif;">Discovered Motifs displayed as sequence logos (in the forward and reverse complement (RC) orientations), along with an E-value for the significance of the result. Clicking on <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">More displays the number of times the motif was identified. <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">To determine if the identified motifs resemble the binding motifs of known transcription factors, we can submit them to [|Tomtom], which searches a database of known motifs to find potential matches and provides a statistical measure of motif-motif similarity. We can run the analysis individually for each motif prediction, or analyze them all at once. <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">You will be shown a status page describing the inputs and the selected parameters, as well as a link to the results at the top of the screen. Clicking the link displays an output page that is continually updated until the analysis has completed. Like DREME, Tomtom will also email you the results. <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">The HTML output for Tomtom shows a list of possible motif matches for each DREME motif prediction generated from your Nanog regions. Clicking on the match name shows you an alignment of the predicted motif versus the match.
 * 1) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Click on <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Fetch Sequences followed by <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Extract Genomic DNA.
 * 2) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Select the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Encode-Nanog.narrowPeak file.
 * 3) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Change the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Source for Genomic Data to History, and then select <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">hg19-chr12-reference_sequence.fasta.
 * 4) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Click <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Execute to generate a FASTA file containing the genomic DNA given by the coordinates in the file.
 * 1) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Visit the [|DREME website] and select the downloaded Nanog peaks as input to DREME
 * 2) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Enter your email address so that DREME can email you once the analysis is complete
 * 3) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Enter a job description so you will recognize which job has been emailed to you and then start the search.
 * 1) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Scroll to the bottom of the motif list and click <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Submit All . A dialog box will appear asking you to <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Select what you want to do or <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Select a program . Select <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Tomtom and click <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Send.
 * 2) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">This takes you to the input page. Tomtom allows you to select the database you wish to search against. We will use the default but keep in mind that there are other options when performing your own analyses.
 * 3) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Enter your email address and job description and start the search.
 * __ ** What are the possible identities for the motifs you discovered from the Nanog ChIP regions on chr12? Do these make sense given what you know about the Nanog transcription factor? Compare your results to results for the full ENCODE data set available from [|FactorBook] and take a look at a [|Nanog binding matrix](from [|Ferraris //et. al.//, 2013]). ** __

<span style="color: #333333; font-family: 'Proxima Nova Bold','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 28px;">Functional enrichment
<span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">A large number of methods and external tools are available to make sense of a gene list, most of which revolve around detecting functional enrichment in Gene Ontology, curated pathways or networks. A recent post on [|Getting Genetics Done] should get you started, and we will expand this section for future workshops. For now we will explore one approach that cover the basics. <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">In order to be able to identify functional enrichment we need to look at genome wide signals; instead of just looking at the ENCODE Nanog peak calls for chromosome 12 we will use the full set of calls. We will use [|GREAT] to perform a functional enrichment. GREAT takes a list of regions, associates them with nearby genes, and then analyzes the gene annotations to assign biological meaning to the data: <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">GREAT provides the output in HTML format organized by section. <span style="color: #333333; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px; text-decoration: none;"> <span style="color: #444444; font-family: 'Proxima Nova Regular',Helvetica,Arial,Helvetica,sans-serif; font-size: 16px;">The first section contains the <span style="font-family: 'Proxima Nova Semi-Bold','Helvetica Neue',Arial,Helvetica,sans-serif;">Job Description detailing the inputs and parameters. The <span style="font-family: 'Proxima Nova Semi-Bold','Helvetica Neue',Arial,Helvetica,sans-serif;">Region-Gene Association Graphs indicate the number of associated genes per region, and the distances of the regions to transcription start sites (TSS), providing insights into the genomic organization of Nanog ChIP regions.
 * 1) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Go back to Galaxy and download the BED file you created from the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Encode-Nanog.narrowPeak file for the motif analysis. Make sure it is the BED formatted file as this is a requirement for the next step.
 * 2) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">First select <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Human: GRCh37 (UCSC hg19, Feb/2009) for the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Species Assembly.
 * 3) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Choose the ENCODE-Nanog.bed file and use the <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Whole genome for <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Background regions.
 * 4) <span style="color: #444444; font-family: 'Proxima Nova Regular','Helvetica Neue',Helvetica,Arial,sans-serif; font-size: 16px;">Click <span style="background-color: #ffb446; color: rgba(0,0,0,0.74902); font-family: Menlo,Consolas,'DejaVu Sans Mono',Monaco,'Courier New',Courier,monospace; font-size: 0.8em;">Submit.
 * __ ** Explore the rest of the results. The output contains enrichment of GO terms, pathway / disease enrichments, miRNA and TF-targets. It is also worthwhile to click on the enriched annotation names to get a more detailed view of genes annotated with these terms. ** __

=** References **=

1) http://ngs.molgen.mpg.de/ngsuploads/Cornelius/ESGI/Chip.htm 2) http://jura.wi.mit.edu/bio/education/hot_topics/galaxy/GalaxyNov2012_ChIP-seq_toPost.pdf 3) https://sites.google.com/site/princetonhtseq/tutorials/chip-seq 4) http://scriptogr.am/ohofmann/chip-seq (this wiki was constructed from this)