Quantifying ChIP-seq data:A spiking method providing an internal reference for sample-to-sample normalization

Chromatin immunoprecipitation followed by deep sequencing (ChIP-seq) experiments are widely used to determine, within entire genomes, the occupancy sites of any protein of interest, including, for example, transcription factors, RNA polymerases, or histones with or without various modifications. In addition to allowing the determination of occupancy sites within one cell type and under one condition, this method allows, in principle, the establishment and comparison of occupancy maps in various cell types, tissues, and conditions. Such comparisons require, however, that samples be normalized. Widely used normalization methods that include a quantile normalization step perform well when factor occupancy varies at a subset of sites, but may miss uniform genome-wide increases or decreases in site occupancy. We describe a spike adjustment procedure (SAP) that, unlike commonly used normalization methods intervening at the analysis stage, entails an experimental step prior to immunoprecipitation. A constant, low amount from a single batch of chromatin of a foreign genome is added to the experimental chromatin. This "spike" chromatin then serves as an internal control to which the experimental signals can be adjusted. We show that the method improves similarity between replicates and reveals biological differences including global and largely uniform changes

The loss of circadian PAR bZip transcription factors results in epilepsy

DBP (albumin D-site-binding protein), HLF (hepatic leukemia factor), and TEF (thyrotroph embryonic factor) are the three members of the PAR bZip (proline and acidic amino acid-rich basic leucine zipper) transcription factor family. All three of these transcriptional regulatory proteins accumulate with robust circadian rhythms in tissues with high amplitudes of clock gene expression, such as the suprachiasmatic nucleus (SCN) and the liver. However, they are expressed at nearly invariable levels in most brain regions, in which clock gene expression only cycles with low amplitude. Here we show that mice deficient for all three PAR bZip proteins are highly susceptible to generalized spontaneous and audiogenic epilepsies that frequently are lethal. Transcriptome profiling revealed pyridoxal kinase (Pdxk) as a target gene of PAR bZip proteins in both liver and brain. Pyridoxal kinase converts vitamin B6 derivatives into pyridoxal phosphate (PLP), the coenzyme of many enzymes involved in amino acid and neurotransmitter metabolism. PAR bZip-deficient mice show decreased brain levels of PLP, serotonin, and dopamine, and such changes have previously been reported to cause epilepsies in other systems. Hence, the expression of some clock-controlled genes, such as Pdxk, may have to remain within narrow limits in the brain. This could explain why the circadian oscillator has evolved to generate only low-amplitude cycles in most brain regions

Letter from JFD to AEK, May 16, 1942

Daniel thanks Kober for sending him her sign-list.Classic

A multiplicity of factors contributes to selective RNA polymerase III occupancy of a subset of RNA polymerase III genes in mouse liver

The genomic loci occupied by RNA polymerase (RNAP) III have been characterized in human culture cells by genome-wide chromatin immunoprecipitations, followed by deep sequencing (ChIP-seq). These studies have shown that only ∼40% of the annotated 622 human tRNA genes and pseudogenes are occupied by RNAP-III, and that these genes are often in open chromatin regions rich in active RNAP-II transcription units. We have used ChIP-seq to characterize RNAP-III-occupied loci in a differentiated tissue, the mouse liver. Our studies define the mouse liver RNAP-III-occupied loci including a conserved mammalian interspersed repeat (MIR) as a potential regulator of an RNAP-III subunit-encoding gene. They reveal that synteny relationships can be established between a number of human and mouse RNAP-III genes, and that the expression levels of these genes are significantly linked. They establish that variations within the A and B promoter boxes, as well as the strength of the terminator sequence, can strongly affect RNAP-III occupancy of tRNA genes. They reveal correlations with various genomic features that explain the observed variation of 81% of tRNA scores. In mouse liver, loci represented in the NCBI37/mm9 genome assembly that are clearly occupied by RNAP-III comprise 50 Rn5s (5S RNA) genes, 14 known non-tRNA RNAP-III genes, nine Rn4.5s (4.5S RNA) genes, and 29 SINEs. Moreover, out of the 433 annotated tRNA genes, half are occupied by RNAP-III. Transfer RNA gene expression levels reflect both an underlying genomic organization conserved in dividing human culture cells and resting mouse liver cells, and the particular promoter and terminator strengths of individual genes