Galaxy Dnpatterntools for Computational Analysis of Nucleosome Positioning Sequence Patterns.

Nucleosomes are basic units of DNA packing in eukaryotes. Their structure is well conserved from yeast to human and consists of the histone octamer core and 147 bp DNA wrapped around it. Nucleosomes are bound to a majority of the eukaryotic genomic DNA, including its regulatory regions. Hence, they also play a major role in gene regulation. For the latter, their precise positioning on DNA is essential. In the present paper, we describe Galaxy dnpatterntools-software package for nucleosome DNA sequence analysis and mapping. This software will be useful for computational biologists practitioners to conduct more profound studies of gene regulatory mechanisms.

Identification of the conserved long non-coding RNAs in myogenesis.

BACKGROUND: Our understanding of genome regulation is ever-evolving with the continuous discovery of new modes of gene regulation, and transcriptomic studies of mammalian genomes have revealed the presence of a considerable population of non-coding RNA molecules among the transcripts expressed. One such non-coding RNA molecule is long non-coding RNA (lncRNA). However, the function of lncRNAs in gene regulation is not well understood; moreover, finding conserved lncRNA across species is a challenging task. Therefore, we propose a novel approach to identify conserved lncRNAs and functionally annotate these molecules. RESULTS: In this study, we exploited existing myogenic transcriptome data and identified conserved lncRNAs in mice and humans. We identified the lncRNAs expressing differentially between the early and later stages of muscle development. Differential expression of these lncRNAs was confirmed experimentally in cultured mouse muscle C2C12 cells. We utilized the three-dimensional architecture of the genome and identified topologically associated domains for these lncRNAs. Additionally, we correlated the expression of genes in domains for functional annotation of these trans-lncRNAs in myogenesis. Using this approach, we identified conserved lncRNAs in myogenesis and functionally annotated them. CONCLUSIONS: With this novel approach, we identified the conserved lncRNAs in myogenesis in humans and mice and functionally annotated them. The method identified a large number of lncRNAs are involved in myogenesis. Further studies are required to investigate the reason for the conservation of the lncRNAs in human and mouse while their sequences are dissimilar. Our approach can be used to identify novel lncRNAs conserved in different species and functionally annotated them.

Nucleosome positioning sequence patterns as packing or regulatory.

Nucleosome positioning DNA sequence patterns (NPS)-usually distributions of particular dinucleotides or other sequence elements in nucleosomal DNA-at least partially determine chromatin structure and arrangements of nucleosomes that in turn affect gene expression. Statistically, NPS are defined as oscillations of the dinucleotide periodicity of about 10 base pairs (bp) which reflects the double helix period. We compared the nucleosomal DNA patterns in mouse, human and yeast organisms and observed few distinctive patterns that can be termed as packing and regulatory referring to distinctive modes of chromatin function. For the first time the NPS patterns in nucleus accumbens cells (NAC) in mouse brain were characterized and compared to the patterns in human CD4+ and apoptotic lymphocyte cells and well studied patterns in yeast. The NPS patterns in human CD4+ cells and mouse brain cells had very high positive correlation. However, there was no correlation between them and patterns in human apoptotic lymphocyte cells and yeast, but the latter two were highly correlated with each other. By their dinucleotide arrangements the analyzed NPS patterns classified into stable canonical WW/SS (W = A or T and S = C or G dinucleotide) and less stable RR/YY (R = A or G and Y = C or T dinucleotide) patterns and anti-patterns. In the anti-patterns positioning of the dinucleotides is flipped compared to those in the regular patterns. Stable canonical WW/SS patterns and anti-patterns are ubiquitously observed in many organisms and they had high resemblance between yeast and human apoptotic cells. Less stable RR/YY patterns had higher positive correlation between mouse and normal human cells. Our analysis and evidence from scientific literature lead to idea that various distinct patterns in nucleosomal DNA can be related to the two roles of the chromatin: packing (WW/SS) and regulatory (RR/YY and "anti").

Exploring a Drosophila Transcription Factor Interaction Network to Identify Cis-Regulatory Modules.

Multiple transcription factors (TFs) bind to specific sites in the genome and interact among themselves to form the cis-regulatory modules (CRMs). They are essential in modulating the expression of genes, and it is important to study this interplay to understand gene regulation. In the present study, we integrated experimentally identified TF binding sites collected from published studies with computationally predicted TF binding sites to identify Drosophila CRMs. Along with the detection of the previously known CRMs, this approach identified novel protein combinations. We determined high-occupancy target sites, where a large number of TFs bind. Investigating these sites revealed that Giant, Dichaete, and Knirp are highly enriched in these locations. A common TAG team motif was observed at these sites, which might play a role in recruiting other TFs. While comparing the binding sites at distal and proximal promoters, we found that certain regulatory TFs, such as Zelda, were highly enriched in enhancers. Our study has shown that, from the information available concerning the TF binding sites, the real CRMs could be predicted accurately and efficiently. Although we only may claim co-occurrence of these proteins in this study, it may actually point to their interaction (as known interaction proteins typically co-occur together). Such an integrative approach can, therefore, help us to provide a better understanding of the interplay among the factors, even though further experimental verification is required.

Snf2h Drives Chromatin Remodeling to Prime Upper Layer Cortical Neuron Development.

Alterations in the homeostasis of either cortical progenitor pool, namely the apically located radial glial (RG) cells or the basal intermediate progenitors (IPCs) can severely impair cortical neuron production. Such changes are reflected by microcephaly and are often associated with cognitive defects. Genes encoding epigenetic regulators are a frequent cause of intellectual disability and many have been shown to regulate progenitor cell growth, including our inactivation of the Smarca1 gene encoding Snf2l, which is one of two ISWI mammalian orthologs. Loss of the Snf2l protein resulted in dysregulation of Foxg1 and IPC proliferation leading to macrocephaly. Here we show that inactivation of the closely related Smarca5 gene encoding the Snf2h chromatin remodeler is necessary for embryonic IPC expansion and subsequent specification of callosal projection neurons. Telencephalon-specific Smarca5 cKO embryos have impaired cell cycle kinetics and increased cell death, resulting in fewer Tbr2+ and FoxG1+ IPCs by mid-neurogenesis. These deficits give rise to adult mice with a dramatic reduction in Satb2+ upper layer neurons, and partial agenesis of the corpus callosum. Mice survive into adulthood but molecularly display reduced expression of the clustered protocadherin genes that may further contribute to altered dendritic arborization and a hyperactive behavioral phenotype. Our studies provide novel insight into the developmental function of Snf2h-dependent chromatin remodeling processes during brain development.

Insights into interplay between rexinoid signaling and myogenic regulatory factor-associated chromatin state in myogenic differentiation.

While skeletal myogenesis is tightly coordinated by myogenic regulatory factors including MyoD and myogenin, chromatin modifications have emerged as vital mechanisms of myogenic regulation. We have previously established that bexarotene, a clinically approved agonist of retinoid X receptor (RXR), promotes the specification and differentiation of skeletal muscle lineage. Here, we examine the genome-wide impact of rexinoids on myogenic differentiation through integral RNA-seq and ChIP-seq analyses. We found that bexarotene promotes myoblast differentiation through the coordination of exit from the cell cycle and the activation of muscle-related genes. We uncovered a new mechanism of rexinoid action which is mediated by the nuclear receptor and largely reconciled through a direct regulation of MyoD gene expression. In addition, we determined a rexinoid-responsive residue-specific histone acetylation at a distinct chromatin state associated to MyoD and myogenin. Thus, we provide novel molecular insights into the interplay between RXR signaling and chromatin states pertinent to myogenic programs in early myoblast differentiation.

Structural dynamics of nucleosome mediated by acetylations at H3K56 and H3K115,122.

Post translational modifications have a profound role in the regulation of several biological processes such as transcription, replication, and DNA repair. Acetylation and phosphorylation form a major class of post translational modifications involved in nucleosomal regulation by modifying its structure. The effect of post translational modifications on nucleosome structure could be better explored when the molecular trajectories explaining the time dependent structural evolution over a period of time is examined at the atomic level. The present study attempts to highlight the importance of acetylation, especially at entry-exit (Lys56) and dyad (Lys115 and Lys122) regions in regulating the nucleosome accessibility and mobility using all atom simulations. It is evident from this study that acetylation at Lys56, Lys115, and Lys122 introduces local changes in the electrostatic nature of the lateral surface and thereby weakens the histone-DNA interactions. In addition, simulations also reveal significant changes in the dynamics of superhelical DNA. The acetylation at Lys56 promotes a high amplitude out-of-planar movement of entry-exit termini. Whereas, acetylation at Lys115 and Lys122 increases the flexibility of the superhelical DNA to facilitate the rolling of the superhelical DNA around the octameric histone. In essence, the present study highlights the role of acetylation at Lys56, Lys115, and Lys122 in transcriptional regulation by promoting high amplitude dynamics of superhelical DNA for a possible unwrapping as well as mobility of nucleosome.

Drosophila H2A and H2A.Z Nucleosome Sequences Reveal Different Nucleosome Positioning Sequence Patterns.

Nucleosomes are implicated in transcriptional regulation as well as in packing and stabilizing the DNA. Nucleosome positions affect the transcription by impeding or facilitating the binding of transcription factors. The DNA sequence, especially the periodic occurrences of dinucleotides, is a major factor that affects the nucleosome positioning. We analyzed the Drosophila DNA sequences bound by H2A and H2A.Z nucleosomes. Periodic patterns of dinucleotides (weak-weak/strong-strong or purine-purine/pyrimidine-pyrimidine) were identified as WW/SS and RR/YY nucleosome positioning sequence (NPS) patterns. The WW/SS NPS pattern of the H2A nucleosome has a 10-bp period of weak-weak/strong-strong (W = A or T; S = G or C) dinucleotides. The 10-bp periodicity, however, is disrupted in the middle of the sequence. At the dyad, the SS dinucleotide is preferred. On the other hand, the RR/YY NPS pattern has an 18-bp periodicity of purine-purine/pyrimidine-pyrimidine (R = A or G; Y = T or C) dinucleotides. The NPS patterns from H2A.Z nucleosomes differ from the NPS patterns from H2A nucleosomes. The RR/YY pattern of H2A.Z nucleosomes has major peaks shifted by 10 bp deviated from the H2A nucleosome pattern. The H2A and H2A.Z nucleosomes have different sequence preferences. The shifted peaks coincide with DNA regions interacting with the histone loops.

Voluntary Running Triggers VGF-Mediated Oligodendrogenesis to Prolong the Lifespan of Snf2h-Null Ataxic Mice.

Exercise has been argued to enhance cognitive function and slow progressive neurodegenerative disease. Although exercise promotes neurogenesis, oligodendrogenesis and adaptive myelination are also significant contributors to brain repair and brain health. Nonetheless, the molecular details underlying these effects remain poorly understood. Conditional ablation of the Snf2h gene impairs cerebellar development producing mice with poor motor function, progressive ataxia, and death between postnatal days 25-45. Here, we show that voluntary running induced an endogenous brain repair mechanism that resulted in a striking increase in hindbrain myelination and the long-term survival of Snf2h cKO mice. Further experiments identified the VGF growth factor as a major driver underlying this effect. VGF neuropeptides promote oligodendrogenesis in vitro, whereas Snf2h cKO mice treated with full-length VGF-encoding adenoviruses removed the requirement of exercise for survival. Together, these results suggest that VGF delivery could represent a therapeutic strategy for cerebellar ataxia and other pathologies of the CNS.

Nucleosome Presence at AML-1 Binding Sites Inversely Correlates with Ly49 Expression: Revelations from an Informatics Analysis of Nucleosomes and Immune Cell Transcription Factors.

Beyond its role in genomic organization and compaction, the nucleosome is believed to participate in the regulation of gene transcription. Here, we report a computational method to evaluate the nucleosome sensitivity for a transcription factor over a given stretch of the genome. Sensitive factors are predicted to be those with binding sites preferentially contained within nucleosome boundaries and lacking 10 bp periodicity. Based on these criteria, the Acute Myeloid Leukemia-1a (AML-1a) transcription factor, a regulator of immune gene expression, was identified as potentially sensitive to nucleosomal regulation within the mouse Ly49 gene family. This result was confirmed in RMA, a cell line with natural expression of Ly49, using MNase-Seq to generate a nucleosome map of chromosome 6, where the Ly49 gene family is located. Analysis of this map revealed a specific depletion of nucleosomes at AML-1a binding sites in the expressed Ly49A when compared to the other, silent Ly49 genes. Our data suggest that nucleosome-based regulation contributes to the expression of Ly49 genes, and we propose that this method of predicting nucleosome sensitivity could aid in dissecting the regulatory role of nucleosomes in general.

Apoptotic lymphocytes of H. sapiens lose nucleosomes in GC-rich promoters.

We analyzed two sets of human CD4+ nucleosomal DNA directly sequenced by Illumina (Solexa) high throughput sequencing method. The first set has ∼40 M sequences and was produced from the normal CD4+ T lymphocytes by micrococcal nuclease. The second set has ∼44 M sequences and was obtained from peripheral blood lymphocytes by apoptotic nucleases. The different nucleosome sets showed similar dinucleotide positioning AA/TT, GG/CC, and RR/YY (R is purine, Y--pyrimidine) patterns with periods of 10-10.4 bp. Peaks of GG/CC and AA/TT patterns were shifted by 5 bp from each other. Two types of promoters in H. sapiens: AT and GC-rich were identified. AT-rich promoters in apoptotic cell had +1 nucleosome shifts 50-60 bp downstream from those in normal lymphocytes. GC-rich promoters in apoptotic cells lost 80% of nucleosomes around transcription start sites as well as in total DNA. Nucleosome positioning was predicted by combination of {AA, TT}, {GG, CC}, {WW, SS} and {RR, YY} patterns. In our study we found that the combinations of {AA, TT} and {GG, CC} provide the best results and successfully mapped 33% of nucleosomes 147 bp long with precision ±15 bp (only 31/147 or 21% is expected).

Snf2h-mediated chromatin organization and histone H1 dynamics govern cerebellar morphogenesis and neural maturation.

Chromatin compaction mediates progenitor to post-mitotic cell transitions and modulates gene expression programs, yet the mechanisms are poorly defined. Snf2h and Snf2l are ATP-dependent chromatin remodelling proteins that assemble, reposition and space nucleosomes, and are robustly expressed in the brain. Here we show that mice conditionally inactivated for Snf2h in neural progenitors have reduced levels of histone H1 and H2A variants that compromise chromatin fluidity and transcriptional programs within the developing cerebellum. Disorganized chromatin limits Purkinje and granule neuron progenitor expansion, resulting in abnormal post-natal foliation, while deregulated transcriptional programs contribute to altered neural maturation, motor dysfunction and death. However, mice survive to young adulthood, in part from Snf2l compensation that restores Engrailed-1 expression. Similarly, Purkinje-specific Snf2h ablation affects chromatin ultrastructure and dendritic arborization, but alters cognitive skills rather than motor control. Our studies reveal that Snf2h controls chromatin organization and histone H1 dynamics for the establishment of gene expression programs underlying cerebellar morphogenesis and neural maturation.

Identification of cis-regulatory modules in promoters of human genes exploiting mutual positioning of transcription factors.

In higher organisms, gene regulation is controlled by the interplay of non-random combinations of multiple transcription factors (TFs). Although numerous attempts have been made to identify these combinations, important details, such as mutual positioning of the factors that have an important role in the TF interplay, are still missing. The goal of the present work is in silico mapping of some of such associating factors based on their mutual positioning, using computational screening. We have selected the process of myogenesis as a study case, and we focused on TF combinations involving master myogenic TF Myogenic differentiation (MyoD) with other factors situated at specific distances from it. The results of our work show that some muscle-specific factors occur together with MyoD within the range of ±100 bp in a large number of promoters. We confirm co-occurrence of the MyoD with muscle-specific factors as described in earlier studies. However, we have also found novel relationships of MyoD with other factors not specific for muscle. Additionally, we have observed that MyoD tends to associate with different factors in proximal and distal promoter areas. The major outcome of our study is establishing the genome-wide connection between biological interactions of TFs and close co-occurrence of their binding sites.

Optimized position weight matrices in prediction of novel putative binding sites for transcription factors in the Drosophila melanogaster genome.

Position weight matrices (PWMs) have become a tool of choice for the identification of transcription factor binding sites in DNA sequences. DNA-binding proteins often show degeneracy in their binding requirement and thus the overall binding specificity of many proteins is unknown and remains an active area of research. Although existing PWMs are more reliable predictors than consensus string matching, they generally result in a high number of false positive hits. Our previous study introduced a promising approach to PWM refinement in which known motifs are used to computationally mine putative binding sites directly from aligned promoter regions using composition of similar sites. In the present study, we extended this technique originally tested on single examples of transcription factors (TFs) and showed its capability to optimize PWM performance to predict new binding sites in the fruit fly genome. We propose refined PWMs in mono- and dinucleotide versions similarly computed for a large variety of transcription factors of Drosophila melanogaster. Along with the addition of many auxiliary sites the optimization includes variation of the PWM motif length, the binding sites location on the promoters and the PWM score threshold. To assess the predictive performance of the refined PWMs we compared them to conventional TRANSFAC and JASPAR sources. The results have been verified using performed tests and literature review. Overall, the refined PWMs containing putative sites derived from real promoter content processed using optimized parameters had better general accuracy than conventional PWMs.

Human milk metagenome: a functional capacity analysis.

BACKGROUND: Human milk contains a diverse population of bacteria that likely influences colonization of the infant gastrointestinal tract. Recent studies, however, have been limited to characterization of this microbial community by 16S rRNA analysis. In the present study, a metagenomic approach using Illumina sequencing of a pooled milk sample (ten donors) was employed to determine the genera of bacteria and the types of bacterial open reading frames in human milk that may influence bacterial establishment and stability in this primal food matrix. The human milk metagenome was also compared to that of breast-fed and formula-fed infants' feces (n = 5, each) and mothers' feces (n = 3) at the phylum level and at a functional level using open reading frame abundance. Additionally, immune-modulatory bacterial-DNA motifs were also searched for within human milk. RESULTS: The bacterial community in human milk contained over 360 prokaryotic genera, with sequences aligning predominantly to the phyla of Proteobacteria (65%) and Firmicutes (34%), and the genera of Pseudomonas (61.1%), Staphylococcus (33.4%) and Streptococcus (0.5%). From assembled human milk-derived contigs, 30,128 open reading frames were annotated and assigned to functional categories. When compared to the metagenome of infants' and mothers' feces, the human milk metagenome was less diverse at the phylum level, and contained more open reading frames associated with nitrogen metabolism, membrane transport and stress response (P < 0.05). The human milk metagenome also contained a similar occurrence of immune-modulatory DNA motifs to that of infants' and mothers' fecal metagenomes. CONCLUSIONS: Our results further expand the complexity of the human milk metagenome and enforce the benefits of human milk ingestion on the microbial colonization of the infant gut and immunity. Discovery of immune-modulatory motifs in the metagenome of human milk indicates more exhaustive analyses of the functionality of the human milk metagenome are warranted.

Dynamics of modeled oligonucleosomes and the role of histone variant proteins in nucleosome organization.

Elucidation of the structural dynamics of a nucleosome is of primary importance for understanding the molecular mechanisms that control the nucleosomal positioning. The presence of variant histone proteins in the nucleosome core raises the functional diversity of the nucleosomes in gene regulation and has the profound epigenetic consequences of great importance for understanding the fundamental issues like the assembly of variant nucleosomes, chromatin remodeling, histone posttranslational modifications, etc. Here, we report our observation of the dominant mechanisms of relaxation motions of the oligonucleosomes such as dimer, trimer, and tetramer (in the beads on a string model) with conventional core histones and role of variant histone H2A.Z in the chromatin dynamics using normal mode analysis. Analysis of the directionality of the global dynamics of the oligonucleosome reveals (i) the in-planar stretching as well as out-of-planar bending motions as the relaxation mechanisms of the oligonucleosome and (ii) the freedom of the individual nucleosome in expressing the combination of the above-mentioned motions as the global mode of dynamics. The highly dynamic N-termini of H3 and (H2A.Z-H2B) dimer evidence their participation in the transcriptionally active state. The key role of variant H2A.Z histone as a major source of vibrant motions via weaker intra- and intermolecular correlations is emphasized in this chapter.

Optimizing the GATA-3 position weight matrix to improve the identification of novel binding sites.

BACKGROUND: The identifying of binding sites for transcription factors is a key component of gene regulatory network analysis. This is often done using position-weight matrices (PWMs). Because of the importance of in silico mapping of tentative binding sites, we previously developed an approach for PWM optimization that substantially improves the accuracy of such mapping. RESULTS: The present work implements the optimization algorithm applied to the existing PWM for GATA-3 transcription factor and builds a new di-nucleotide PWM. The existing available PWM is based on experimental data adopted from Jaspar. The optimized PWM substantially improves the sensitivity and specificity of the TF mapping compared to the conventional applications. The refined PWM also facilitates in silico identification of novel binding sites that are supported by experimental data. We also describe uncommon positioning of binding motifs for several T-cell lineage specific factors in human promoters. CONCLUSION: Our proposed di-nucleotide PWM approach outperforms the conventional mono-nucleotide PWM approach with respect to GATA-3. Therefore our new di-nucleotide PWM provides new insight into plausible transcriptional regulatory interactions in human promoters.

Discovery, optimization and validation of an optimal DNA-binding sequence for the Six1 homeodomain transcription factor.

The Six1 transcription factor is a homeodomain protein involved in controlling gene expression during embryonic development. Six1 establishes gene expression profiles that enable skeletal myogenesis and nephrogenesis, among others. While several homeodomain factors have been extensively characterized with regards to their DNA-binding properties, relatively little is known of the properties of Six1. We have used the genomic binding profile of Six1 during the myogenic differentiation of myoblasts to obtain a better understanding of its preferences for recognizing certain DNA sequences. DNA sequence analyses on our genomic binding dataset, combined with biochemical characterization using binding assays, reveal that Six1 has a much broader DNA-binding sequence spectrum than had been previously determined. Moreover, using a position weight matrix optimization algorithm, we generated a highly sensitive and specific matrix that can be used to predict novel Six1-binding sites with highest accuracy. Furthermore, our results support the idea of a mode of DNA recognition by this factor where Six1 itself is sufficient for sequence discrimination, and where Six1 domains outside of its homeodomain contribute to binding site selection. Together, our results provide new light on the properties of this important transcription factor, and will enable more accurate modeling of Six1 function in bioinformatic studies.

Variety of genomic DNA patterns for nucleosome positioning.

Precise positioning of nucleosomes along DNA is important for a variety of gene regulatory processes. Among the factors directing nucleosome positioning, the DNA sequence is highly important. Two main classes of nucleosome positioning sequence (NPS) patterns have previously been described. In the first class, AA, TT, and other WW dinucleotides (where W is A or T) tend to occur together (in-phase) in the major groove of DNA closest to the histone octamer surface, while SS dinucleotides (where S is G or C) are predominantly positioned in the major groove facing outward. In the second class, AA and TT are structurally separated (AA backbone near the histone octamer, and TT backbone further away), but grouped with other RR (where R is purine A or G) and YY (where Y is pyrimidine C or T) dinucleotides. As a result, the RR/YY pattern includes counter-phase AA/TT distributions. We describe here anti-NPS patterns, which are inverse to the conventional NPS patterns: WW runs inverse to SS, and RR inverse to YY. Evidence for the biological relevance of anti-NPS patterns is presented.

A barrier nucleosome model for statistical positioning of nucleosomes throughout the yeast genome.

Most nucleosomes are well-organized at the 5' ends of S. cerevisiae genes where "-1" and "+1" nucleosomes bracket a nucleosome-free promoter region (NFR). How nucleosomal organization is specified by the genome is less clear. Here we establish and inter-relate rules governing genomic nucleosome organization by sequencing DNA from more than one million immunopurified S. cerevisiae nucleosomes (displayed at http://atlas.bx.psu.edu/). Evidence is presented that the organization of nucleosomes throughout genes is largely a consequence of statistical packing principles. The genomic sequence specifies the location of the -1 and +1 nucleosomes. The +1 nucleosome forms a barrier against which nucleosomes are packed, resulting in uniform positioning, which decays at farther distances from the barrier. We present evidence for a novel 3' NFR that is present at >95% of all genes. 3' NFRs may be important for transcription termination and anti-sense initiation. We present a high-resolution genome-wide map of TFIIB locations that implicates 3' NFRs in gene looping.