On, and orientation for the gene nearest to the ChIP-seq alignment. The medium thick dark line is the 59 utr of the gene and the thicker dark region is the first exon followed by a thin line with SPI1005 cost arrows which is intron 1; Conservation, the track of Phastcons for sequence similarity among placental mammals; ChIPseeqer peak, the black rectangular block is the location of the statistically-qualified peak of sequencing alignments called by the ChIPseeqer algorithm; HU Bcl-3, a representation of the.sam alignments for the Bcl-3 ChIPseq of the unloaded muscle; WB Bcl-3, a representation of the.sam alignments for the Bcl-3 ChIPseq of the weight bearing muscle; Input, 12926553 a representation of the.sam alignments for non-ChIP unloaded chromatin; NF-kB site, location of a NF-kB consensus site identified by JASPAR. doi:10.1371/journal.pone.0051478.gResults Characterization of Bcl-3 ChIP-seqSince transcriptional activation of the p50-Bcl-3 complex will not happen without Bcl-3 [11] we reasoned that its binding is the best for following the active NF-kB complex in unloaded muscle at a genome-wide level. Bcl-3 ChIP sequences from unloaded muscle were put through the peak SC 66 site finding algorithm of ChIPseeqer [15], which identifies peaks with increased Bcl-3 binding compared to weight bearing muscle. By using a low stringency peak height cutoff, 49,000 Bcl-3 peaks were found. These peaks were evenly distributed across the mouse genome (Figure S1). Using a web based tool called Nebula [20,21] a component of the Galaxy suite of programs, the distribution of Bcl-3 peaks from unloaded musclewere compared to random peaks found from the input fraction of chromatin (Figure 1). This showed that the Bcl-3 ChIP had succeeded in enriching many Bcl-3 binding sites (i.e., peaks) near the activation sites of transcription across the entire genome. We then took the sequence alignments from unloaded muscle Bcl-3 ChIP-seq and compared them to weight bearing muscle sequences in order to find peaks that were at least 2-fold increased using the peak finder in ChIPseeqer. By using this level of stringency for peak finding we obtained 2,817 Bcl-3 peaks in unloaded compared to weight bearing muscle. Phastcon analysis using the Cistrome/Galaxy program [21] was used to show that the peaks were located at phylogenetically conserved sites (Figure 2). Annotation of the parts of genes associated with unloading-induced peaks showed that they were mainly inA Bcl-3 Network Controls Muscle AtrophyTable 2. qPCR of selected proteolysis genes with increased Bcl-3 promoter binding.Gene Arih2 Ate1 Fbxo6 Itch Rlim Rnf13 Psmb7 Sod1 Trim63 Ubb Ubr1 Fold change, control vs. unloaded. doi:10.1371/journal.pone.0051478.tFold activation 2.1 1.5 1.8 1.4 1.6 1.5 1.9 1.8 2.0 1.8 1.Also found in the GO pathways and shown in Table 1 are genes that function 23388095 in the reduction of reactive oxygen species, including SOD1, and several phosphatases. The other GO terms having genes represented are those involved in regulating myogenesis, particularly in the Wnt pathway, and those in glucose metabolism, including glycogen phosphorylase and 7-phosphofructokinase, genes that liberate glucose and control its glycolytic metabolism respectively. Several of the genes, especially the E3 ligases found as Bcl-3 targets by ChIP-seq were subject to qPCR to verify gene activation during unloading and these data are shown in Table 2. The advantage of iPAGE is that it can find the most important functions of the overrepresented genes ha.On, and orientation for the gene nearest to the ChIP-seq alignment. The medium thick dark line is the 59 utr of the gene and the thicker dark region is the first exon followed by a thin line with arrows which is intron 1; Conservation, the track of Phastcons for sequence similarity among placental mammals; ChIPseeqer peak, the black rectangular block is the location of the statistically-qualified peak of sequencing alignments called by the ChIPseeqer algorithm; HU Bcl-3, a representation of the.sam alignments for the Bcl-3 ChIPseq of the unloaded muscle; WB Bcl-3, a representation of the.sam alignments for the Bcl-3 ChIPseq of the weight bearing muscle; Input, 12926553 a representation of the.sam alignments for non-ChIP unloaded chromatin; NF-kB site, location of a NF-kB consensus site identified by JASPAR. doi:10.1371/journal.pone.0051478.gResults Characterization of Bcl-3 ChIP-seqSince transcriptional activation of the p50-Bcl-3 complex will not happen without Bcl-3 [11] we reasoned that its binding is the best for following the active NF-kB complex in unloaded muscle at a genome-wide level. Bcl-3 ChIP sequences from unloaded muscle were put through the peak finding algorithm of ChIPseeqer [15], which identifies peaks with increased Bcl-3 binding compared to weight bearing muscle. By using a low stringency peak height cutoff, 49,000 Bcl-3 peaks were found. These peaks were evenly distributed across the mouse genome (Figure S1). Using a web based tool called Nebula [20,21] a component of the Galaxy suite of programs, the distribution of Bcl-3 peaks from unloaded musclewere compared to random peaks found from the input fraction of chromatin (Figure 1). This showed that the Bcl-3 ChIP had succeeded in enriching many Bcl-3 binding sites (i.e., peaks) near the activation sites of transcription across the entire genome. We then took the sequence alignments from unloaded muscle Bcl-3 ChIP-seq and compared them to weight bearing muscle sequences in order to find peaks that were at least 2-fold increased using the peak finder in ChIPseeqer. By using this level of stringency for peak finding we obtained 2,817 Bcl-3 peaks in unloaded compared to weight bearing muscle. Phastcon analysis using the Cistrome/Galaxy program [21] was used to show that the peaks were located at phylogenetically conserved sites (Figure 2). Annotation of the parts of genes associated with unloading-induced peaks showed that they were mainly inA Bcl-3 Network Controls Muscle AtrophyTable 2. qPCR of selected proteolysis genes with increased Bcl-3 promoter binding.Gene Arih2 Ate1 Fbxo6 Itch Rlim Rnf13 Psmb7 Sod1 Trim63 Ubb Ubr1 Fold change, control vs. unloaded. doi:10.1371/journal.pone.0051478.tFold activation 2.1 1.5 1.8 1.4 1.6 1.5 1.9 1.8 2.0 1.8 1.Also found in the GO pathways and shown in Table 1 are genes that function 23388095 in the reduction of reactive oxygen species, including SOD1, and several phosphatases. The other GO terms having genes represented are those involved in regulating myogenesis, particularly in the Wnt pathway, and those in glucose metabolism, including glycogen phosphorylase and 7-phosphofructokinase, genes that liberate glucose and control its glycolytic metabolism respectively. Several of the genes, especially the E3 ligases found as Bcl-3 targets by ChIP-seq were subject to qPCR to verify gene activation during unloading and these data are shown in Table 2. The advantage of iPAGE is that it can find the most important functions of the overrepresented genes ha.