ChIP-Seq workflow with analysis in ArrayStar

Lasergene Genomics’ ChIP-Seq analysis workflow enables you to locate the binding sites of DNA-associated proteins and determine how these proteins interact with the DNA to affect expression in nearby genes.

In this tutorial, you will use SeqMan NGen to create a ChIP-Seq assembly for two experimental replicates and one wild type control of E. coli transcription factor FlhD. Once the assembly is complete, you will import the data into ArrayStar’s ChIP-Seq workflow and filter the results for IP fragments that match in both replicates, and that meet a specified signal threshold. Finally, you will look for genes that the fragments intersect, and see whether they share any common characteristics.

Setting up and running the assembly in SeqMan NGen:

1. Download T2_ChIP-Seq.zip (2.4 GB) and extract the contents to any convenient location (e.g., your computer’s desktop). The folder contents consist of:
   
   
   • The E. coli reference sequence: Escherichia coli str. K-12 substr. MG1655.U00096.gbk
     
     
   • The sample sequences: six folders beginning with flhC, flhD and WT.

2. Launch SeqMan NGen and choose New Assembly.

3. Choose RNA-Seq / Transcriptomics on the left, then choose ChIP-Seq on the right.

4. In the Reference Sequence screen, click Add and add the reference sequence MG1655_U00096.3.gb. Press Next.

5. In the Input Sequences screen, uncheck the box next to Paired-end reads. Press the Add Folder button to add the folder starting with WT. Do the same for each of the next two sample folders. Click Next.

6. In the Set Up Experiments screen, check the box to the right of WT_rep1, indicating that it is the control. Press Next.

7. In the Assembly Options screen, check the box next to Maximum total reads and enter 5000000 (5 followed by 6 zeros) to reduce the assembly time. Click Next.

8. In the Analysis Options screen, leave the ChiP-Seq detection method at MACS and press Next.

9. In the Define Binding Proteins screen, use the Known binding site motif drop-down menu to select Transcription Factor Database. Use the Organism drop-down menu to select Bacteria. Press Select and choose FlhCD_CS / FlhCD, then press OK.
   
   
   
   Press Next.

10. In the Assembly Output screen, type “ChIP-seq” into the Project Name text box. Use the Browse button to specify a Project Folder for your assembly output files, then press Next.

11. In the Run Assembly Project screen, check whether the recommendation is to run the assembly locally or on the cloud. Press the corresponding link to begin assembly. Typical assembly times are 20 minutes (cloud) or 30-45 minutes (local).

12. Wait until being informed that assembly has finished, then click Next.

13. From the Assembly Summary screen, click Analyze peaks to open the results in ArrayStar.

Analyzing the results in ArrayStar:

1. In ArrayStar, click on the Fragment Table tab and ensure that approximately 40 IP fragments were found.

2. Similarly, click on the Peak Table tab and ensure that about 50-60 peaks were found.

3. To filter for fragments found in both replicate samples, choose Filter > Filter All and make selections to match the image below.

4. Press the Choose Signal Criteria button and make selections to match the following image. Press OK.

5. Press Search. The search should yield about 17 IP fragments.

6. Press the Select and Show Results in IP Fragment Table tool (), located above the result table.

7. In the Fragment Table, use the Quick-Filter tool to choose Show Only Selected IP Fragments.

8. Click on Add/Manage Columns. In the IP Fragment Info section, choose Intersecting Genes, then >Add Column>, then OK.

9. Look at the newly-added column and note that many of the fragments intersect genes beginning with “flg”, “fli” and “flh”.

These are the major genes for producing flagella in E. coli and are known to be regulated by the transcription factor encoded by the flhC and flhD genes. The che gene is involved in chemotaxis and are also known to be regulated by FlhCD.

This marks the end of this tutorial.