With the H4K12ac dataset, we obtained 5. 53 million total reads, of which 4. 04 million were selleck chem unique reads with an average coverage of 8. 7 reads per promoter and 123 reads per CDS. The higher sequence coverage of H4K5ac in control, 13. 3% more mapped reads compared to FC, may account for the larger number of genes identified in control with our exclusion criteria. The lower coverage in H4K12ac may also explain the smaller percentage of genes found to overlap with H4K5ac. Differential peak calling and data mining analysis Peak finding was performed using a Model based Ana lysis of ChIP Seq algorithm. To determine genes differentially enriched for H4K5ac in the respective groups, we ran MACS on fear conditioned against non fear conditioned control and vice versa.

H4K5ac peaks were identified in MACS with the following parame ters effective genome size 1. 87e 09, tag size 50, bandwidth 300, m fold 4, and P value cutoff 1. 00e 5. We also used the Statistical model for the Identification of chip Enriched Regions to call differentially acetylated peaks between groups. We used the following parameters for SICER redundancy threshold 1, window size 200, fragment size 150, effective genome fraction 0. 7, gap size 400, FDR 1. 00e 3, and filtered post analysis for genes with P value 1. 00e 5. We further compared results to the Genomatix NGS analyzer with Auto Claverie algorithm with the following parameters window size 100 and P value 0. 05, filtered post analysis for genes with P value 1. 00e 5. EpiChip analysis was performed according to standard protocols, except gene scoring was performed 5000 from the 5 start position.

H4K12ac ChIP Seq data, by CFC in young mice, was obtained from the public repository at Galaxy Central. Con trol ChIP Seq data for H4K12ac, for sample or experi mental condition, was not available. Gene ontology and pathway analysis To determine functional gene enrichment and inter action networks of genes differentially acetylated in fear conditioned compared to non fear conditioned controls, we used the genes identified in MACS for functional annotation. From the 241 differentially acet ylated regions identified in fear conditioned over con trol, 115 unique peaks were associated in the promoter or coding region of genes. From the 77 differentially acet ylated regions identified in control over fear conditioned, 42 unique peaks were associated with gene bodies.

We used The Database for Annotation, Visualization and Inte grated Discovery for the analysis of functionally enriched genes in our respective gene lists. Settings were set at a count threshold of 2 and EASE score of 0. 1, a more conservative test than Fishers Exact test. We also used Web based Gene Set Analysis Toolkit V2 for the analysis Drug_discovery of functionally enriched genes in our respective gene lists.