Predicting human nucleosome occupancy from primary sequence

Genome-wide prediction of nucleosome occupancy in maize reveals plant chromatin structural features at genes and other elements at multiple scales

The nucleosome is a fundamental structural and functional chromatin unit that affects nearly all DNA-templated events in eukaryotic genomes. It is also a biochemical substrate for higher order, cis-acting gene expression codes and the monomeric structural unit for chromatin packaging at multiple scales. To predict the nucleosome landscape of a model plant genome, we used a support vector machine computational algorithm trained on human chromatin to predict the nucleosome occupancy likelihood (NOL) across the maize (Zea mays) genome. Experimentally validated NOL plots provide a novel genomic annotation that highlights gene structures, repetitive elements, and chromosome-scale domains likely to reflect regional gene density. We established a new genome browser (http://www.genomaize.org) for viewing support vector machine-based NOL scores. This annotation provides sequence-based comprehensive coverage across the entire genome, including repetitive genomic regions typically excluded from experimental genomics data. We find that transposable elements often displayed family-specific NOL profiles that included distinct regions, especially near their termini, predicted to have strong affinities for nucleosomes. We examined transcription start site consensus NOL plots for maize gene sets and discovered that most maize genes display a typical +1 nucleosome positioning signal just downstream of the start site but not upstream. This overall lack of a -1 nucleosome positioning signal was also predicted by our method for Arabidopsis (Arabidopsis thaliana) genes and verified by additional analysis of previously published Arabidopsis MNase-Seq data, revealing a general feature of plant promoters. Our study advances plant chromatin research by defining the potential contribution of the DNA sequence to observed nucleosome positioning and provides an invariant baseline annotation against which other genomic data can be compared.

Changes in nucleosome position at transcriptional start sites of specific genes in Zea mays mediator of paramutation1 mutants

Nucleosomes facilitate compaction of DNA within the confines of the eukaryotic nucleus. This packaging of DNA and histone proteins must accommodate cellular processes, such as transcription and DNA replication. The repositioning of nucleosomes to facilitate cellular processes is likely regulated by several factors. In Zea mays, Mediator of paramutation1 (MOP1) has been demonstrated to be an epigenetic regulator of gene expression. Based on sequence orthology and mutant phenotypes, MOP1 is likely to function in an RNA-dependent pathway to mediate changes to chromatin. High-resolution microarrays were used to assay the distribution of nucleosomes across the transcription start sites (TSSs) of ~400 maize genes in wild type and mutant mop1-1 tissues. Analysis of nucleosome distribution in leaf, immature tassel and ear shoot tissues resulted in the identification of three genes showing consistent differences in nucleosome positioning and occupancy between wild type and mutant mop1-1. These specific changes in nucleosome distribution were located upstream as well as downstream of the TSS. No direct relationship between the specific changes in nucleosome distribution and transcription were observed through quantitative expression analysis in these tissues. In silico prediction suggests that nucleosome positioning is not dictated by intrinsic DNA sequence signals in the TSSs of two of the identified genes, suggesting a role for chromatin remodeling proteins in MOP1-mediated pathways. These results also indicate that MOP1 contributions to nucleosome position may be either separate from changes in gene expression, or cooperative with development and other levels of regulation in coordinating gene expression.

Whole-genome sequencing in autism identifies hot spots for de novo germline mutation

De novo mutation plays an important role in autism spectrum disorders (ASDs). Notably, pathogenic copy number variants (CNVs) are characterized by high mutation rates. We hypothesize that hypermutability is a property of ASD genes and may also include nucleotide-substitution hot spots. We investigated global patterns of germline mutation by whole-genome sequencing of monozygotic twins concordant for ASD and their parents. Mutation rates varied widely throughout the genome (by 100-fold) and could be explained by intrinsic characteristics of DNA sequence and chromatin structure. Dense clusters of mutations within individual genomes were attributable to compound mutation or gene conversion. Hypermutability was a characteristic of genes involved in ASD and other diseases. In addition, genes impacted by mutations in this study were associated with ASD in independent exome-sequencing data sets. Our findings suggest that regional hypermutation is a significant factor shaping patterns of genetic variation and disease risk in humans.

An analysis and prediction of nucleosome positioning based on information content

Nucleosome positioning plays a key role in the regulation of many biological processes. In this study, the statistical difference of information content was investigated in nucleosome and linker DNA regions across eukaryotic organisms. By analyzing the information redundancy, D k , in Saccharomyces cerevisiae, Drosophila melanogaster, and Caenorhabditis elegans genomes, the short-range dominance of nucleotide correlation in nucleosome and linker DNA regions was confirmed. Significant difference of the D k value between the nucleosome and linker DNA regions was also found. The underlying reason for many successful oligonucleotide-based predictions of nucleosome positioning in eukaryotic model organisms may be attributed to the short-range dominance of nucleotide correlation in the nucleosome and linker DNA regions. When applying power spectrum analysis to the nucleosome and linker DNA regions, some obvious differences in sequence periodic signals were observed. The parameter F k was introduced to describe particular base correlation. Furthermore, the support vector machine combining F k was used to classify nucleosome and linker DNA regions in Homo sapiens, Oryzias latipes, C. elegans, Candida albicans, and S. cerevisiae. Independent test demonstrated that a good performance can be achieved by using this algorithm. This result further revealed that base correlation information has an important role in nucleosome positioning.

The roles of the monomer length and nucleotide context of plant tandem repeats in nucleosome positioning

Similar to regularly spaced nucleosomes in chromatin, long tandem DNA arrays are composed of regularly alternating monomers that have almost identical primary DNA structures. Such a similarity in the structural organization makes these arrays especially interesting for studying the role of intrinsic DNA preferences in nucleosome positioning. We have studied the nucleosome formation potential of DNA tandem repeat families with different monomer lengths (ML). In total, 165 plant tandem repeat families from the PlantSat database (http://w3lamc.umbr.cas.cz/PlantSat/) were divided into two classes based on the number of nucleosome repeats in one DNA monomer. For predicting nucleosome formation potential, we developed the Phase method, which combines the advantages of multiple bioinformatics models. The Phase method was able to distinguish interfamily differences and intrafamily monomer variation and identify the influence of nucleotide context on nucleosome formation potential. Three main types of nucleosome arrangement in DNA tandem repeat arrays--regular, partially regular (partial), and flexible--were distinguished among a great variety of Phase profiles. The regular type, in which all nucleosomes of the monomer array are positioned in a context-dependent manner, is the most representative type of the class 1 families, with ML equal to or a multiple of the nucleosome repeat length (NRL). In the partially regular type, nucleotide context influences the positioning of only a subset of nucleosomes. The influence of the nucleotide context on nucleosome positioning has the least effect in the flexible type, which contains the greatest number of families (65). The majority of these families belong to class 2 and have nonmultiple ML to NRL ratios.

iNuc-PhysChem: a sequence-based predictor for identifying nucleosomes via physicochemical properties

Nucleosome positioning has important roles in key cellular processes. Although intensive efforts have been made in this area, the rules defining nucleosome positioning is still elusive and debated. In this study, we carried out a systematic comparison among the profiles of twelve DNA physicochemical features between the nucleosomal and linker sequences in the Saccharomyces cerevisiae genome. We found that nucleosomal sequences have some position-specific physicochemical features, which can be used for in-depth studying nucleosomes. Meanwhile, a new predictor, called iNuc-PhysChem, was developed for identification of nucleosomal sequences by incorporating these physicochemical properties into a 1788-D (dimensional) feature vector, which was further reduced to a 884-D vector via the IFS (incremental feature selection) procedure to optimize the feature set. It was observed by a cross-validation test on a benchmark dataset that the overall success rate achieved by iNuc-PhysChem was over 96% in identifying nucleosomal or linker sequences. As a web-server, iNuc-PhysChem is freely accessible to the public at http://lin.uestc.edu.cn/server/iNuc-PhysChem. For the convenience of the vast majority of experimental scientists, a step-by-step guide is provided on how to use the web-server to get the desired results without the need to follow the complicated mathematics that were presented just for the integrity in developing the predictor. Meanwhile, for those who prefer to run predictions in their own computers, the predictor's code can be easily downloaded from the web-server. It is anticipated that iNuc-PhysChem may become a useful high throughput tool for both basic research and drug design.

Translation efficiency in yeasts correlates with nucleosome formation in promoters

Elongation efficiency index (EEI) was suggested earlier to estimate gene expression efficiency by nucleotide context of coding sequence in unicellular organisms. We have analyzed association between EEI and nucleosome formation potential (NFP) in 5' regulatory regions upstream translation initiation site (TIS) from two yeast species. Theoretical estimations of NFP based on DNA sequence were obtained by Recon method. Experimental estimation of nucleosome occupancy was obtained by high-throughput sequencing data of nucleosomal DNA in Saccharomyces cerevisiae . For the sample of all genes correlation coefficient was calculated between two vectors: vector of NFP values for fixed position relative to TIS and vector of EEI values. Profiles of correlation coefficients of NFP and EEI were counted in (-600; +600) regions relative to TIS for gene sequences extracted from GenBank. We found regions of strong negative dependence between NFP and EEI for all genes as well as for 10% highly expressed genes in Schizosaccharomyces pombe (10% of EEI-highest genes). At the same time, we found positive dependence between NFP and EEI for all genes and for low expressed genes in S. cerevisiae (10% of EEI-lowest genes). The association between NFP and EEI could be explained by evolutionary selection of context characteristics of nucleotide sequences for gene expression optimization.

Linking genome to epigenome

Recent epigenomic studies have identified significant differences between developmental stages and cell types. While the importance of epigenetic regulation has been increasingly recognized, it remains unclear how the global epigenetic patterns are established and maintained. Here I review a number of recent studies with the emphasis on the role of the genomic sequence in shaping the epigenetic landscape. These studies strongly suggest that the sequence information is important not just for controlling target specificity but for orchestrating the diversity of epigenetic patterns among different cell types. The epigenome is maintained by the complex network of a large number of interactions. Integrative approaches are needed to gain insights into these networks.

Epigenetic changes in fetal hypothalamic energy regulating pathways are associated with maternal undernutrition and twinning

Undernutrition during pregnancy is implicated in the programming of offspring for the development of obesity and diabetes. We hypothesized that maternal programming causes epigenetic changes in fetal hypothalamic pathways regulating metabolism. This study used sheep to examine the effect of moderate maternal undernutrition (60 d before to 30 d after mating) and twinning to investigate changes in the key metabolic regulators proopiomelanocortin (POMC) and the glucocorticoid receptor (GR) in fetal hypothalami. Methylation of the fetal hypothalamic POMC promoter was reduced in underfed singleton, fed twin, and underfed twin groups (60, 73, and 63% decrease, respectively). This was associated with reduced DNA methyltransferase activity and altered histone methylation and acetylation. Methylation of the hypothalamic GR promoter was decreased in both twin groups and in maternally underfed singleton fetuses (52, 65, and 55% decrease, respectively). This correlated with changes in histone methylation and acetylation and increased GR mRNA expression in the maternally underfed singleton group. Alterations in GR were hypothalamic specific, with no changes in hippocampi. Unaltered levels of OCT4 promoter methylation indicated gene-specific effects. In conclusion, twinning and periconceptional undernutrition are associated with epigenetic changes in fetal hypothalamic POMC and GR genes, potentially resulting in altered energy balance regulation in the offspring.

Sequence-based classification using discriminatory motif feature selection

Most existing methods for sequence-based classification use exhaustive feature generation, employing, for example, all -mer patterns. The motivation behind such (enumerative) approaches is to minimize the potential for overlooking important features. However, there are shortcomings to this strategy. First, practical constraints limit the scope of exhaustive feature generation to patterns of length , such that potentially important, longer () predictors are not considered. Second, features so generated exhibit strong dependencies, which can complicate understanding of derived classification rules. Third, and most importantly, numerous irrelevant features are created. These concerns can compromise prediction and interpretation. While remedies have been proposed, they tend to be problem-specific and not broadly applicable. Here, we develop a generally applicable methodology, and an attendant software pipeline, that is predicated on discriminatory motif finding. In addition to the traditional training and validation partitions, our framework entails a third level of data partitioning, a discovery partition. A discriminatory motif finder is used on sequences and associated class labels in the discovery partition to yield a (small) set of features. These features are then used as inputs to a classifier in the training partition. Finally, performance assessment occurs on the validation partition. Important attributes of our approach are its modularity (any discriminatory motif finder and any classifier can be deployed) and its universality (all data, including sequences that are unaligned and/or of unequal length, can be accommodated). We illustrate our approach on two nucleosome occupancy datasets and a protein solubility dataset, previously analyzed using enumerative feature generation. Our method achieves excellent performance results, with and without optimization of classifier tuning parameters. A Python pipeline implementing the approach is available at http://www.epibiostat.ucsf.edu/biostat/sen/dmfs/.

Exploratory analysis of genomic segmentations with Segtools

Background

As genome-wide experiments and annotations become more prevalent, researchers increasingly require tools to help interpret data at this scale. Many functional genomics experiments involve partitioning the genome into labeled segments, such that segments sharing the same label exhibit one or more biochemical or functional traits. For example, a collection of ChlP-seq experiments yields a compendium of peaks, each labeled with one or more associated DNA-binding proteins. Similarly, manually or automatically generated annotations of functional genomic elements, including cis-regulatory modules and protein-coding or RNA genes, can also be summarized as genomic segmentations.

Results

We present a software toolkit called Segtools that simplifies and automates the exploration of genomic segmentations. The software operates as a series of interacting tools, each of which provides one mode of summarization. These various tools can be pipelined and summarized in a single HTML page. We describe the Segtools toolkit and demonstrate its use in interpreting a collection of human histone modification data sets and Plasmodium falciparum local chromatin structure data sets.

Conclusions

Segtools provides a convenient, powerful means of interpreting a genomic segmentation.

Evolutionary trends of GC/AT distribution patterns in promoters

Nucleotide distributions in genomes is known not to be random, showing the presence of specific motifs, long and short range correlations, periodicities, etc. Particularly, motifs are critical for the recognition by specific proteins affecting chromosome organization, transcription and DNA replication but little is known about the possible functional effects of nucleotide distributions on the conformational landscape of DNA, putatively leading to differential selective pressures throughout evolution. Promoter sequences have a fundamental role in the regulation of gene activity and a vast literature suggests that their conformational landscapes may be a critical factor in gene expression dynamics. On these grounds, with the aim of investigating the putative existence of phylogenetic patterns of promoter base distributions, we analyzed GC/AT ratios along the 1000 nucleotide sequences upstream of TSS in wide sets of promoters belonging to organisms ranging from bacteria to pluricellular eukaryotes. The data obtained showed very clear phylogenetic trends throughout evolution of promoter sequence base distributions. Particularly, in all cases either GC-rich or AT-rich monotone gradients were observed: the former being present in eukaryotes, the latter in bacteria along with strand biases. Moreover, within eukaryotes, GC-rich gradients increased in length from unicellular organisms to plants, to vertebrates and, within them, from ancestral to more recent species. Finally, results were thoroughly discussed with particular attention to the possible correlation between nucleotide distribution patterns, evolution, and the putative existence of differential selection pressures, deriving from structural and/or functional constraints, between and within prokaryotes and eukaryotes.

Predicting nucleosome positioning in genomes: physical and bioinformatic approaches

In eukaryotic genomes, nucleosomes are responsible for packaging DNA and controlling gene expression. For this reason, an increasing interest is arising on computational methods capable of predicting the nucleosome positioning along genomes. In this review we describe and compare bioinformatic and physical approaches adopted to predict nucleosome occupancy along genomes. Computational analyses attempt at decoding the experimental nucleosome maps of genomes in terms of certain dinucleotide step periodicity observed along DNA. Such investigations show that highly significant information about the occurrence of a nucleosome along DNA is intrinsic in certain features of the sequence suggesting that DNA of eukaryotic genomes encodes nucleosome organization. Besides the bioinformatic approaches, physical models were proposed based on the sequence dependent conformational features of the DNA chain, which govern the free energy needed to transform recurrent DNA tracts along the genome into the nucleosomal shape.

Formaldehyde-assisted isolation of regulatory elements

Formaldehyde-Assisted Isolation of Regulatory Elements (FAIRE) is based on locus-specific variations in the ability of protein components of chromatin to trap genomic DNA following formaldehyde treatment. This variation is mostly due to uneven nucleosome distribution since histones are the most abundant and highly crosslinkable components of chromatin. The method can identify and enrich for physically accessible DNA segments of the eukaryotic genome corresponding to known regulatory regions and regions that might have thus far unidentified structural role in the nuclear organization of chromatin. The enrichment patterns are cell type specific and thus might provide information about how transcriptional systems are organized and regulated in various tissues and how they might be disrupted in disease states. Analysis of a 268 kb region of chromosome 19 in human fibroblasts shown here demonstrates that while most DNA fragments detected by FAIRE correspond to sites of DNaseI hypersensitivity in active regions of chromatin, some are found in otherwise repressed chromatin domains and at other sites that are not found with other methods used to probe chromatin structure. Further exploration of FAIRE is warrented due to the simplicity of the protocol and recent advancements in massively parallel sequencing.

A human genome structural variation sequencing resource reveals insights into mutational mechanisms

Understanding the prevailing mutational mechanisms responsible for human genome structural variation requires uniformity in the discovery of allelic variants and precision in terms of breakpoint delineation. We develop a resource based on capillary end sequencing of 13.8 million fosmid clones from 17 human genomes and characterize the complete sequence of 1054 large structural variants corresponding to 589 deletions, 384 insertions, and 81 inversions. We analyze the 2081 breakpoint junctions and infer potential mechanism of origin. Three mechanisms account for the bulk of germline structural variation: microhomology-mediated processes involving short (2-20 bp) stretches of sequence (28%), nonallelic homologous recombination (22%), and L1 retrotransposition (19%). The high quality and long-range continuity of the sequence reveals more complex mutational mechanisms, including repeat-mediated inversions and gene conversion, that are most often missed by other methods, such as comparative genomic hybridization, single nucleotide polymorphism microarrays, and next-generation sequencing.

Pre-replication complex proteins assemble at regions of low nucleosome occupancy within the Chinese hamster dihydrofolate reductase initiation zone

Genome-scale mapping of pre-replication complex proteins has not been reported in mammalian cells. Poor enrichment of these proteins at specific sites may be due to dispersed binding, poor epitope availability or cell cycle stage-specific binding. Here, we have mapped sites of biotin-tagged ORC and MCM protein binding in G1-synchronized populations of Chinese hamster cells harboring amplified copies of the dihydrofolate reductase (DHFR) locus, using avidin-affinity purification of biotinylated chromatin followed by high-density microarray analysis across the DHFR locus. We have identified several sites of significant enrichment for both complexes distributed throughout the previously identified initiation zone. Analysis of the frequency of initiations across stretched DNA fibers from the DHFR locus confirmed a broad zone of de-localized initiation activity surrounding the sites of ORC and MCM enrichment. Mapping positions of mononucleosomal DNA empirically and computing nucleosome-positioning information in silico revealed that ORC and MCM map to regions of low measured and predicted nucleosome occupancy. Our results demonstrate that specific sites of ORC and MCM enrichment can be detected within a mammalian intitiation zone, and suggest that initiation zones may be regions of generally low nucleosome occupancy where flexible nucleosome positioning permits flexible pre-RC assembly sites.

Sequence-dependent histone variant positioning signatures

Background

Nucleosome, the fundamental unit of chromatin, is formed by wrapping nearly 147bp of DNA around an octamer of histone proteins. This histone core has many variants that are different from each other by their biochemical compositions as well as biological functions. Although the deposition of histone variants onto chromatin has been implicated in many important biological processes, such as transcription and replication, the mechanisms of how they are deposited on target sites are still obscure.

Results

By analyzing genomic sequences of nucleosomes bearing different histone variants from human, including H2A.Z, H3.3 and both (H3.3/H2A.Z, so-called double variant histones), we found that genomic sequence contributes in part to determining target sites for different histone variants. Moreover, dinucleotides CA/TG are remarkably important in distinguishing target sites of H2A.Z-only nucleosomes with those of H3.3-containing (both H3.3-only and double variant) nucleosomes.

Conclusions

There exists a DNA-related mechanism regulating the deposition of different histone variants onto chromatin and biological outcomes thereof. This provides additional insights into epigenetic regulatory mechanisms of many important cellular processes.

Contribution of histone sequence preferences to nucleosome organization: proposed definitions and methodology

We propose definitions and procedures for comparing nucleosome maps and discuss current agreement and disagreement on the effect of histone sequence preferences on nucleosome organization in vivo.

Evaluating genome-scale approaches to eukaryotic DNA replication

Mechanisms regulating where and when eukaryotic DNA replication initiates remain a mystery. Recently, genome-scale methods have been brought to bear on this problem. The identification of replication origins and their associated proteins in yeasts is a well-integrated investigative tool, but corresponding data sets from multicellular organisms are scarce. By contrast, standardized protocols for evaluating replication timing have generated informative data sets for most eukaryotic systems. Here, I summarize the genome-scale methods that are most frequently used to analyse replication in eukaryotes, the kinds of questions each method can address and the technical hurdles that must be overcome to gain a complete understanding of the nature of eukaryotic replication origins.

Nucleosome positioning: an essential component of the enhancer regulatory code?

The organization of transcription factor binding sites within enhancers plays a key role in the function of these elements. A recent study suggests that the intrinsic propensity to bind to nucleosomes can also play a decisive part in determining enhancer function.

Prediction of nucleosome positioning in genomes: limits and perspectives of physical and bioinformatic approaches

Nucleosomes, the fundamental repeating subunits of all eukaryotic chromatin, are responsible for packaging DNA into chromosomes inside the cell nucleus and controlling gene expression. While it has been well established that nucleosomes exhibit higher affinity for select DNA sequences, until recently it was unclear whether such preferences exerted a significant, genome-wide effect on nucleosome positioning in vivo. For this reason, an increasing interest is arising on a wide-ranging series of experimental and computational analyses capable of predicting the nucleosome positioning along genomes. Toward this goal, we propose a theoretical model for predicting nucleosome thermodynamic stability in terms of DNA sequence. Based on a statistical mechanical approach, the model allows the calculation of the sequence-dependent canonical ensemble free energy involved in nucleosome formation. The theoretical free energies were evaluated for 90 single nucleosome DNA tracts and successfully compared with those obtained with nucleosome competitive reconstitution. These results, obtained for single nucleosomes, could in principle allow the calculation of the intrinsic affinity of nucleosomes along DNA sequences virtually opening the possibility of predicting the nucleosome positioning along genomes on physical basis. The theoretical nucleosome distribution was compared and validated with that of yeast and human genome experimentally determined. The results interpret on a physical basis the experimental nucleosome positioning and are comparable with those obtained adopting models based on the identification of some recurrent sequence features obtained from the statistical analysis of a very large pool of nucleosomal DNA sequences provided by the positioning maps of genomes.

Learning a weighted sequence model of the nucleosome core and linker yields more accurate predictions in Saccharomyces cerevisiae and Homo sapiens

DNA in eukaryotes is packaged into a chromatin complex, the most basic element of which is the nucleosome. The precise positioning of the nucleosome cores allows for selective access to the DNA, and the mechanisms that control this positioning are important pieces of the gene expression puzzle. We describe a large-scale nucleosome pattern that jointly characterizes the nucleosome core and the adjacent linkers and is predominantly characterized by long-range oscillations in the mono, di- and tri-nucleotide content of the DNA sequence, and we show that this pattern can be used to predict nucleosome positions in both Homo sapiens and Saccharomyces cerevisiae more accurately than previously published methods. Surprisingly, in both H. sapiens and S. cerevisiae, the most informative individual features are the mono-nucleotide patterns, although the inclusion of di- and tri-nucleotide features results in improved performance. Our approach combines a much longer pattern than has been previously used to predict nucleosome positioning from sequence—301 base pairs, centered at the position to be scored—with a novel discriminative classification approach that selectively weights the contributions from each of the input features. The resulting scores are relatively insensitive to local AT-content and can be used to accurately discriminate putative dyad positions from adjacent linker regions without requiring an additional dynamic programming step and without the attendant edge effects and assumptions about linker length modeling and overall nucleosome density. Our approach produces the best dyad-linker classification results published to date in H. sapiens, and outperforms two recently published models on a large set of S. cerevisiae nucleosome positions. Our results suggest that in both genomes, a comparable and relatively small fraction of nucleosomes are well-positioned and that these positions are predictable based on sequence alone. We believe that the bulk of the remaining nucleosomes follow a statistical positioning model.

DNA in eukaryotes is packaged into a chromatin complex, the most basic element of which is the nucleosome. The precise positioning of the nucleosome cores allows for selective access to the DNA, and the mechanisms that control this positioning are important pieces of the gene expression puzzle. In this work, we describe a large-scale DNA sequence pattern that jointly characterizes the sequence preferences of the nucleosome core and the adjacent linkers. We show that this pattern can be used to predict nucleosome positions in both H. sapiens and S. cerevisiae more accurately than previously published methods. The model is most accurate in predicting the most stably positioned nucleosomes, and describes a sequence composition pattern that determines a locally optimal dyad (nucleosomal DNA mid-point) position. In contrast to some previous models, this model is not based primarily on excluding poly-A/T sequences, nor does the model prefer 10 bp periodicity. Our results suggest that local sequence composition is one of many factors that direct the positioning of nucleosomes, while dynamic processes such as transcriptional elongation and the actions of chromatin remodeling complexes also play a significant role in the overall chromatin landscape.

The genomic distribution and local context of coincident SNPs in human and chimpanzee

We have previously shown that there is an excess of sites that are polymorphic at orthologous positions in humans and chimpanzees and that this is most likely due to cryptic variation in the mutation rate. We showed that this might be a consequence of complex context effects since we found significant heterogeneity in triplet frequencies around coincident single nucleotide polymorphism (SNP) sites. Here, we show that the heterogeneity in triplet frequencies is not specifically associated with coincident SNPs but is instead driven by base composition bias around CpG dinucleotides. As a result, we suggest that cryptic variation in the mutation rate is truly cryptic, in the sense that the mutation rate does not appear to depend on any specific primary sequence context. Furthermore, we propose that the patterns around CpG dinucleotides are driven by the mutability of CpG dinucleotides in different DNA contexts. We also show that the genomic distribution of coincident SNPs is nonuniform and that there are some subtle differences between the distributions of single and coincident SNPs. Furthermore, we identify regions that contain high numbers of coincident SNPs and suggest that one in particular, a region containing the gene PRIM2, may be under balancing selection.

Contrasting chromatin organization of CpG islands and exons in the human genome

BACKGROUND

CpG islands and nucleosome-free regions are both found in promoters. However, their association has never been studied. On the other hand, DNA methylation is absent in promoters but is enriched in gene bodies. Intragenic nucleosomes and their modifications have been recently associated with RNA splicing. Because the function of intragenic DNA methylation remains unclear, I explored the possibility of its involvement in splicing regulation.

RESULTS

Here I show that CpG islands were associated not only with methylation-free promoters but also with nucleosome-free promoters. Nucleosome-free regions were observed only in promoters containing a CpG island. However, the DNA sequences of CpG islands predicted the opposite pattern, implying a limitation of sequence programs for the determination of nucleosome occupancy. In contrast to the methylation-and nucleosome-free states of CpG-island promoters, exons were densely methylated at CpGs and packaged into nucleosomes. Exon-enrichment of DNA methylation was specifically found in spliced exons and in exons with weak splice sites. The enrichment patterns were less pronounced in initial exons and in non-coding exons, potentially reflecting a lower need for their splicing. I also found that nucleosomes, DNA methylation, and H3K36me3 marked the exons of transcripts with low, medium, and high gene expression levels, respectively.

CONCLUSIONS

Human promoters containing a CpG island tend to remain nucleosome-free as well as methylation-free. In contrast, exons demonstrate a high degree of methylation and nucleosome occupancy. Exonic DNA methylation seems to function together with exonic nucleosomes and H3K36me3 for the proper splicing of transcripts with different expression levels.

Effects of Alu elements on global nucleosome positioning in the human genome

Background

Understanding the genome sequence-specific positioning of nucleosomes is essential to understand various cellular processes, such as transcriptional regulation and replication. As a typical example, the 10-bp periodicity of AA/TT and GC dinucleotides has been reported in several species, but it is still unclear whether this feature can be observed in the whole genomes of all eukaryotes.

Results

With Fourier analysis, we found that this is not the case: 84-bp and 167-bp periodicities are prevalent in primates. The 167-bp periodicity is intriguing because it is almost equal to the sum of the lengths of a nucleosomal unit and its linker region. After masking Alu elements, these periodicities were greatly diminished. Next, using two independent large-scale sets of nucleosome mapping data, we analyzed the distribution of nucleosomes in the vicinity of Alu elements and showed that (1) there are one or two fixed slot(s) for nucleosome positioning within the Alu element and (2) the positioning of neighboring nucleosomes seems to be in phase, more or less, with the presence of Alu elements. Furthermore, (3) these effects of Alu elements on nucleosome positioning are consistent with inactivation of promoter activity in Alu elements.

Conclusions

Our discoveries suggest that the principle governing nucleosome positioning differs greatly across species and that the Alu family is an important factor in primate genomes.

Characterization of transcription start sites of putative non-coding RNAs by multifaceted use of massively paralleled sequencer

On the basis of integrated transcriptome analysis, we show that not all transcriptional start site clusters (TSCs) in the intergenic regions (iTSCs) have the same properties; thus, it is possible to discriminate the iTSCs that are likely to have biological relevance from the other noise-level iTSCs. We used a total of 251 933 381 short-read sequence tags generated from various types of transcriptome analyses in order to characterize 6039 iTSCs, which have significant expression levels. We analyzed and found that 23% of these iTSCs were located in the proximal regions of the RefSeq genes. These RefSeq-linked iTSCs showed similar expression patterns with the neighboring RefSeq genes, had widely fluctuating transcription start sites and lacked ordered nucleosome positioning. These iTSCs seemed not to form independent transcriptional units, simply representing the by-products of the neighboring RefSeq genes, in spite of their significant expression levels. Similar features were also observed for the TSCs located in the antisense regions of the RefSeq genes. Furthermore, for the remaining iTSCs that were not associated with any RefSeq genes, we demonstrate that integrative interpretation of the transcriptome data provides essential information to specify their biological functions in the hypoxic responses of the cells.

The impact of intragenic CpG content on gene expression

The development of vaccine components or recombinant therapeutics critically depends on sustained expression of the corresponding transgene. This study aimed to determine the contribution of intragenic CpG content to expression efficiency in transiently and stably transfected mammalian cells. Based upon a humanized version of green fluorescent protein (GFP) containing 60 CpGs within its coding sequence, a CpG-depleted variant of the GFP reporter was established by carefully modulating the codon usage. Interestingly, GFP reporter activity and detectable protein amounts in stably transfected CHO and 293 cells were significantly decreased upon CpG depletion and independent from promoter usage (CMV, EF1α). The reduction in protein expression associated with CpG depletion was likewise observed for other unrelated reporter genes and was clearly reflected by a decline in mRNA copy numbers rather than translational efficiency. Moreover, decreased mRNA levels were neither due to nuclear export restrictions nor alternative splicing or mRNA instability. Rather, the intragenic CpG content influenced de novo transcriptional activity thus implying a common transcription-based mechanism of gene regulation via CpGs. Increased high CpG transcription correlated with changed nucleosomal positions in vitro albeit histone density at the two genes did not change in vivo as monitored by ChIP.

Schizosaccharomyces pombe genome-wide nucleosome mapping reveals positioning mechanisms distinct from those of Saccharomyces cerevisiae

Positioned nucleosomes limit the access of proteins to DNA and implement regulatory features encoded in eukaryotic genomes. Here we have generated the first genome-wide nucleosome positioning map for Schizosaccharomyces pombe and annotated transcription start and termination sites genome wide. Using this resource, we found surprising differences from the previously published nucleosome organization of the distantly related yeast Saccharomyces cerevisiae. DNA sequence guides nucleosome positioning differently: for example, poly(dA-dT) elements are not enriched in S. pombe nucleosome-depleted regions. Regular nucleosomal arrays emanate more asymmetrically-mainly codirectionally with transcription-from promoter nucleosome-depleted regions, but promoters harboring the histone variant H2A.Z also show regular arrays upstream of these regions. Regular nucleosome phasing in S. pombe has a very short repeat length of 154 base pairs and requires a remodeler, Mit1, that is conserved in humans but is not found in S. cerevisiae. Nucleosome positioning mechanisms are evidently not universal but evolutionarily plastic.

The organization of nucleosomes around splice sites

The occupancy of nucleosomes along chromosome is a key factor for gene regulation. However, except promoter regions, genome-wide properties and functions of nucleosome organization remain unclear in mammalian genomes. Using the computational model of Increment of Diversity with Quadratic Discriminant (IDQD) trained from the microarray data, the nucleosome occupancy score (NOScore) was defined and applied to splice junction regions of constitutive, cassette exon, alternative 3′ and 5′ splicing events in the human genome. We found an interesting relation between NOScore and RNA splicing: exon regions have higher NOScores compared with their flanking intron sequences in both constitutive and alternative splicing events, indicating the stronger nucleosome occupation potential of exon regions. In addition, NOScore valleys present at ∼25 bp upstream of the acceptor site in all splicing events. By defining folding diversity-to-energy ratio to describe RNA structural flexibility, we demonstrated that primary RNA transcripts from nucleosome occupancy regions are relatively rigid and those from nucleosome depleted regions are relatively flexible. The negative correlation between nucleosome occupation/depletion of DNA sequence and structural flexibility/rigidity of its primary transcript around splice junctions may provide clues to the deeper understanding of the unexpected role for nucleosome organization in the regulation of RNA splicing.

A statistical thermodynamic approach for predicting the sequence-dependent nucleosome positioning along genomes

Nucleosomes are the fundamental repeating unit of chromatin and constitute the structural building blocks of the eukaryotic genome. The distribution of nucleosomes along the genome is a significant aspect of chromatin structure and influences gene regulation through modulation of DNA accessibility. For this reason, an increasing interest is arising in models capable of predicting the nucleosome positioning along genomes. Toward this goal, we propose a theoretical model for predicting nucleosome thermodynamic stability in terms of DNA sequence. The model, based on a statistical mechanical approach allows the calculation of the canonical ensemble free energy involved in nucleosome formation. The theoretical free energies were evaluated for about one hundred nucleosome DNA tracts and successfully compared with those obtained in different laboratories with nucleosome competitive reconstitution (correlation coefficient equal to 0.92). We extended these results to the nucleosome positioning along genomes. To test our model, the theoretical nucleosome distribution was compared with that of yeast genome experimentally determined. The results are comparable with those obtained by different authors adopting models based on identifying some recurrent sequence features obtained from the statistical analysis of a very large pool of nucleosomal DNA sequences provided by the positioning maps of genomes.

Nucleosome landscape and control of transcription in the human malaria parasite

In eukaryotic cells, chromatin reorganizes within promoters of active genes to allow the transcription machinery and various transcription factors to access DNA. In this model, promoter-specific transcription factors bind DNA to initiate the production of mRNA in a tightly regulated manner. In the case of the human malaria parasite, Plasmodium falciparum, specific transcription factors are apparently underrepresented with regards to the size of the genome, and mechanisms underlying transcriptional regulation are controversial. Here, we investigate the modulation of DNA accessibility by chromatin remodeling during the parasite infection cycle. We have generated genome-wide maps of nucleosome occupancy across the parasite erythrocytic cycle using two complementary assays--the formaldehyde-assisted isolation of regulatory elements to extract protein-free DNA (FAIRE) and the MNase-mediated purification of mononucleosomes to extract histone-bound DNA (MAINE), both techniques being coupled to high-throughput sequencing. We show that chromatin architecture undergoes drastic upheavals throughout the parasite's cycle, contrasting with targeted chromatin reorganization usually observed in eukaryotes. Chromatin loosens after the invasion of the red blood cell and then repacks prior to the next cycle. Changes in nucleosome occupancy within promoter regions follow this genome-wide pattern, with a few exceptions such as the var genes involved in virulence and genes expressed at early stages of the cycle. We postulate that chromatin structure and nucleosome turnover control massive transcription during the erythrocytic cycle. Our results demonstrate that the processes driving gene expression in Plasmodium challenge the classical eukaryotic model of transcriptional regulation occurring mostly at the transcription initiation level.

Sequence-dependent DNA helical rise and nucleosome stability

BACKGROUND

Nucleosomes are the basic structural units of eukaryotic chromatin and play a key role in regulation of gene expression. After resolution of the nucleosome structure, the bipartite nature of this particle has revealed itself and has disclosed the presence, on the histone surface, of a symmetric distribution of positive charges, able to interact with their negative DNA phosphate counterpart.

RESULTS

We analyzed helical steps in known nucleosomal DNA sequences, observing a significant relationship between their symmetric distribution and nucleosome stability. Synthetic DNA sequences able to form stable nucleosomes were used to compare distances on the left and on the right side of the nucleosomal dyad axis, where DNA phosphates and charged residues of the (H3H4)2-tetramer interact. We observed a linear relationship between coincidence of distances and nucleosome stability, i. e., the more symmetric these distances the more stable the nucleosome.

CONCLUSION

Curves related to this symmetric distribution along the DNA sequence identify preferential sites for positioning of the dyad axis, which we termed palinstases. The comparison of our data with known nucleosome positions in archaeal and eukaryotic sequences shows many coincidences of location. Sequences that impair nucleosome formation and DNase I hypersensitive sites yield curves with a lower degree of symmetry. Analysis performed on DNA tracts of promoters close to the transcription start and termination sites identified peculiar patterns: in particular low affinity for nucleosome binding at the transcription start site and a high affinity exactly at the transcription termination site, suggesting a major role of nucleosomes in the termination of transcription.

Epigenetic silencing in Friedreich ataxia is associated with depletion of CTCF (CCCTC-binding factor) and antisense transcription

BACKGROUND

Over 15 inherited diseases are caused by expansion of triplet-repeats. Friedreich ataxia (FRDA) patients are homozygous for an expanded GAA triplet-repeat sequence in intron 1 of the FXN gene. The expanded GAA triplet-repeat results in deficiency of FXN gene transcription, which is reversed via administration of histone deacetylase inhibitors indicating that transcriptional silencing is at least partially due to an epigenetic abnormality.

METHODOLOGY/PRINCIPAL FINDINGS

We found a severe depletion of the chromatin insulator protein CTCF (CCCTC-binding factor) in the 5'UTR of the FXN gene in FRDA, and coincident heterochromatin formation involving the +1 nucleosome via enrichment of H3K9me3 and recruitment of heterochromatin protein 1. We identified FAST-1 (FXNAntisense Transcript - 1), a novel antisense transcript that overlaps the CTCF binding site in the 5'UTR, which was expressed at higher levels in FRDA. The reciprocal relationship of deficient FXN transcript and higher levels of FAST-1 seen in FRDA was reproduced in normal cells via knockdown of CTCF.

CONCLUSIONS/SIGNIFICANCE

CTCF depletion constitutes an epigenetic switch that results in increased antisense transcription, heterochromatin formation and transcriptional deficiency in FRDA. These findings provide a mechanistic basis for the transcriptional silencing of the FXN gene in FRDA, and broaden our understanding of disease pathogenesis in triplet-repeat diseases.

Stable transcription activities dependent on an orientation of Tam3 transposon insertions into Antirrhinum and yeast promoters occur only within chromatin

Transposon insertions occasionally occur in the promoter regions of plant genes, many of which are still capable of being transcribed. However, it remains unclear how transcription of such promoters is able to occur. Insertion of the Tam3 transposon into various genes of Antirrhinum majus can confer leaky phenotypes without its excision. These genes, named Tam3-permissible alleles, often contain Tam3 in their promoter regions. Two alleles at different anthocyanin biosynthesis loci, nivea(recurrensTam3) (niv(rec)) and pallida(recurrensTam3) (pal(rec)), both contain Tam3 at a similar position immediately upstream of the promoter TATA-box; however, these insertions had different phenotypic consequences. Under conditions where the inserted Tam3 is immobilized, the niv(rec) line produces pale red petals, whereas the pal(rec) line produces no pigment. These pigmentation patterns are correlated with the level of transcripts from the niv(rec) or pal(rec) alleles, and these transcriptional activities are independent of DNA methylation in their promoter regions. In niv(rec), Tam3 is inserted in an orientation that results in the 3' end of Tam3 adjacent to the 5' region of the gene coding sequence. In contrast, the pal(rec) allele contains a Tam3 insertion in the opposite orientation. Four of five different nonrelated genes that are also Tam3-permissible alleles and contain Tam3 within the promoter region share the same Tam3 orientation as niv(rec). The different transcriptional activities dependent on Tam3 orientation in the Antirrhinum promoters were consistent with expression of luciferase reporter constructs introduced into yeast chromosomes but not with transient expression of these constructs in Antirrhinum cells. These results suggest that for Tam3 to sustain stable transcriptional activity in various promoters it must be embedded in chromatin.

From DNA sequence to transcriptional behaviour: a quantitative approach

Complex transcriptional behaviours are encoded in the DNA sequences of gene regulatory regions. Advances in our understanding of these behaviours have been recently gained through quantitative models that describe how molecules such as transcription factors and nucleosomes interact with genomic sequences. An emerging view is that every regulatory sequence is associated with a unique binding affinity landscape for each molecule and, consequently, with a unique set of molecule-binding configurations and transcriptional outputs. We present a quantitative framework based on existing methods that unifies these ideas. This framework explains many experimental observations regarding the binding patterns of factors and nucleosomes and the dynamics of transcriptional activation. It can also be used to model more complex phenomena such as transcriptional noise and the evolution of transcriptional regulation.

Ride the wavelet: A multiscale analysis of genomic contexts flanking small insertions and deletions

Recent studies have revealed that insertions and deletions (indels) are more different in their formation than previously assumed. What remains enigmatic is how the local DNA sequence context contributes to these differences. To investigate the relative impact of various molecular mechanisms to indel formation, we analyzed sequence contexts of indels in the non protein- or RNA-coding, nonrepetitive (NCNR) portion of the human genome. We considered small (<or=30-bp) indels occurring in the human lineage since its divergence from chimpanzee and used wavelet techniques to study, simultaneously for multiple scales, the spatial patterns of short sequence motifs associated with indel mutagenesis. In particular, we focused on motifs associated with DNA polymerase activity, topoisomerase cleavage, double-strand breaks (DSBs), and their repair. We came to the following conclusions. First, many motifs are characterized by unique enrichment profiles in the vicinity of indels vs. indel-free portions of the genome, verifying the importance of sequence context in indel mutagenesis. Second, only limited similarity in motif frequency profiles is evident flanking insertions vs. deletions, confirming differences in their mutagenesis. Third, substantial similarity in frequency profiles exists between pairs of individual motifs flanking insertions (and separately deletions), suggesting "cooperation" among motifs, and thus molecular mechanisms, during indel formation. Fourth, the wavelet analyses demonstrate that all these patterns are highly dependent on scale (the size of an interval considered). Finally, our results depict a model of indel mutagenesis comprising both replication and recombination (via repair of paused replication forks and site-specific recombination).

Nucleosome positioning and gene regulation: advances through genomics

Knowing the precise locations of nucleosomes in a genome is key to understanding how genes are regulated. Recent 'next generation' ChIP-chip and ChIP-Seq technologies have accelerated our understanding of the basic principles of chromatin organization. Here we discuss what high-resolution genome-wide maps of nucleosome positions have taught us about how nucleosome positioning demarcates promoter regions and transcriptional start sites, and how the composition and structure of promoter nucleosomes facilitate or inhibit transcription. A detailed picture is starting to emerge of how diverse factors, including underlying DNA sequences and chromatin remodelling complexes, influence nucleosome positioning.