The native cistrome and sequence motif families of the maize ear.

Elucidating the transcriptional regulatory networks that underlie growth and development requires robust ways to define the complete set of transcription factor (TF) binding sites. Although TF-binding sites are known to be generally located within accessible chromatin regions (ACRs), pinpointing these DNA regulatory elements globally remains challenging. Current approaches primarily identify binding sites for a single TF (e.g. ChIP-seq), or globally detect ACRs but lack the resolution to consistently define TF-binding sites (e.g. DNAse-seq, ATAC-seq). To address this challenge, we developed MNase-defined cistrome-Occupancy Analysis (MOA-seq), a high-resolution (< 30 bp), high-throughput, and genome-wide strategy to globally identify putative TF-binding sites within ACRs. We used MOA-seq on developing maize ears as a proof of concept, able to define a cistrome of 145,000 MOA footprints (MFs). While a substantial majority (76%) of the known ATAC-seq ACRs intersected with the MFs, only a minority of MFs overlapped with the ATAC peaks, indicating that the majority of MFs were novel and not detected by ATAC-seq. MFs were associated with promoters and significantly enriched for TF-binding and long-range chromatin interaction sites, including for the well-characterized FASCIATED EAR4, KNOTTED1, and TEOSINTE BRANCHED1. Importantly, the MOA-seq strategy improved the spatial resolution of TF-binding prediction and allowed us to identify 215 motif families collectively distributed over more than 100,000 non-overlapping, putatively-occupied binding sites across the genome. Our study presents a simple, efficient, and high-resolution approach to identify putative TF footprints and binding motifs genome-wide, to ultimately define a native cistrome atlas.

Topologically associating domains are stable units of replication-timing regulation.

Eukaryotic chromosomes replicate in a temporal order known as the replication-timing program. In mammals, replication timing is cell-type-specific with at least half the genome switching replication timing during development, primarily in units of 400-800 kilobases ('replication domains'), whose positions are preserved in different cell types, conserved between species, and appear to confine long-range effects of chromosome rearrangements. Early and late replication correlate, respectively, with open and closed three-dimensional chromatin compartments identified by high-resolution chromosome conformation capture (Hi-C), and, to a lesser extent, late replication correlates with lamina-associated domains (LADs). Recent Hi-C mapping has unveiled substructure within chromatin compartments called topologically associating domains (TADs) that are largely conserved in their positions between cell types and are similar in size to replication domains. However, TADs can be further sub-stratified into smaller domains, challenging the significance of structures at any particular scale. Moreover, attempts to reconcile TADs and LADs to replication-timing data have not revealed a common, underlying domain structure. Here we localize boundaries of replication domains to the early-replicating border of replication-timing transitions and map their positions in 18 human and 13 mouse cell types. We demonstrate that, collectively, replication domain boundaries share a near one-to-one correlation with TAD boundaries, whereas within a cell type, adjacent TADs that replicate at similar times obscure replication domain boundaries, largely accounting for the previously reported lack of alignment. Moreover, cell-type-specific replication timing of TADs partitions the genome into two large-scale sub-nuclear compartments revealing that replication-timing transitions are indistinguishable from late-replicating regions in chromatin composition and lamina association and accounting for the reduced correlation of replication timing to LADs and heterochromatin. Our results reconcile cell-type-specific sub-nuclear compartmentalization and replication timing with developmentally stable structural domains and offer a unified model for large-scale chromosome structure and function.

Regulated large-scale nucleosome density patterns and precise nucleosome positioning correlate with V(D)J recombination.

We show that the physical distribution of nucleosomes at antigen receptor loci is subject to regulated cell type-specific and lineage-specific positioning and correlates with the accessibility of these gene segments to recombination. At the Ig heavy chain locus (IgH), a nucleosome in pro-B cells is generally positioned over each IgH variable (VH) coding segment, directly adjacent to the recombination signal sequence (RSS), placing the RSS in a position accessible to the recombination activating gene (RAG) recombinase. These changes result in establishment of a specific chromatin organization at the RSS that facilitates accessibility of the genomic DNA for the RAG recombinase. In contrast, in mouse embryonic fibroblasts the coding segment is depleted of nucleosomes, which instead cover the RSS, thereby rendering it inaccessible. Pro-T cells exhibit a pattern intermediate between pro-B cells and mouse embryonic fibroblasts. We also find large-scale variations of nucleosome density over hundreds of kilobases, delineating chromosomal domains within IgH, in a cell type-dependent manner. These findings suggest that developmentally regulated changes in nucleosome location and occupancy, in addition to the known chromatin modifications, play a fundamental role in regulating V(D)J recombination. Nucleosome positioning-which has previously been observed to vary locally at individual enhancers and promoters-may be a more general mechanism by which cells can regulate the accessibility of the genome during development, at scales ranging from several hundred base pairs to many kilobases.

iSeg: an efficient algorithm for segmentation of genomic and epigenomic data.

BACKGROUND: Identification of functional elements of a genome often requires dividing a sequence of measurements along a genome into segments where adjacent segments have different properties, such as different mean values. Despite dozens of algorithms developed to address this problem in genomics research, methods with improved accuracy and speed are still needed to effectively tackle both existing and emerging genomic and epigenomic segmentation problems. RESULTS: We designed an efficient algorithm, called iSeg, for segmentation of genomic and epigenomic profiles. iSeg first utilizes dynamic programming to identify candidate segments and test for significance. It then uses a novel data structure based on two coupled balanced binary trees to detect overlapping significant segments and update them simultaneously during searching and refinement stages. Refinement and merging of significant segments are performed at the end to generate the final set of segments. By using an objective function based on the p-values of the segments, the algorithm can serve as a general computational framework to be combined with different assumptions on the distributions of the data. As a general segmentation method, it can segment different types of genomic and epigenomic data, such as DNA copy number variation, nucleosome occupancy, nuclease sensitivity, and differential nuclease sensitivity data. Using simple t-tests to compute p-values across multiple datasets of different types, we evaluate iSeg using both simulated and experimental datasets and show that it performs satisfactorily when compared with some other popular methods, which often employ more sophisticated statistical models. Implemented in C++, iSeg is also very computationally efficient, well suited for large numbers of input profiles and data with very long sequences. CONCLUSIONS: We have developed an efficient general-purpose segmentation tool and showed that it had comparable or more accurate results than many of the most popular segment-calling algorithms used in contemporary genomic data analysis. iSeg is capable of analyzing datasets that have both positive and negative values. Tunable parameters allow users to readily adjust the statistical stringency to best match the biological nature of individual datasets, including widely or sparsely mapped genomic datasets or those with non-normal distributions.

The spring-loaded genome: nucleosome redistributions are widespread, transient, and DNA-directed.

Nucleosome occupancy plays a key role in regulating access to eukaryotic genomes. Although various chromatin regulatory complexes are known to regulate nucleosome occupancy, the role of DNA sequence in this regulation remains unclear, particularly in mammals. To address this problem, we measured nucleosome distribution at high temporal resolution in human cells at hundreds of genes during the reactivation of Kaposi's sarcoma-associated herpesvirus (KSHV). We show that nucleosome redistribution peaks at 24 h post-KSHV reactivation and that the nucleosomal redistributions are widespread and transient. To clarify the role of DNA sequence in these nucleosomal redistributions, we compared the genes with altered nucleosome distribution to a sequence-based computer model and in vitro-assembled nucleosomes. We demonstrate that both the predicted model and the assembled nucleosome distributions are concordant with the majority of nucleosome redistributions at 24 h post-KSHV reactivation. We suggest a model in which loci are held in an unfavorable chromatin architecture and "spring" to a transient intermediate state directed by DNA sequence information. We propose that DNA sequence plays a more considerable role in the regulation of nucleosome positions than was previously appreciated. The surprising findings that nucleosome redistributions are widespread, transient, and DNA-directed shift the current perspective regarding regulation of nucleosome distribution in humans.

Differential nuclease sensitivity profiling of chromatin reveals biochemical footprints coupled to gene expression and functional DNA elements in maize.

The eukaryotic genome is organized into nucleosomes, the fundamental units of chromatin. The positions of nucleosomes on DNA regulate protein-DNA interactions and in turn influence DNA-templated events. Despite the increasing number of genome-wide maps of nucleosome position, how global changes in gene expression relate to changes in nucleosome position is poorly understood. We show that in nucleosome occupancy mapping experiments in maize (Zea mays), particular genomic regions are highly susceptible to variation introduced by differences in the extent to which chromatin is digested with micrococcal nuclease (MNase). We exploited this digestion-linked variation to identify protein footprints that are hypersensitive to MNase digestion, an approach we term differential nuclease sensitivity profiling (DNS-chip). Hypersensitive footprints were enriched at the 5' and 3' ends of genes, associated with gene expression levels, and significantly overlapped with conserved noncoding sequences and the binding sites of the transcription factor KNOTTED1. We also found that the tissue-specific regulation of gene expression was linked to tissue-specific hypersensitive footprints. These results reveal biochemical features of nucleosome organization that correlate with gene expression levels and colocalize with functional DNA elements. This approach to chromatin profiling should be broadly applicable to other species and should shed light on the relationships among chromatin organization, protein-DNA interactions, and genome regulation.

Label-Free Relative Quantitation of Isobaric and Isomeric Human Histone H2A and H2B Variants by Fourier Transform Ion Cyclotron Resonance Top-Down MS/MS.

Histone variants are known to play a central role in genome regulation and maintenance. However, many variants are inaccessible by antibody-based methods or bottom-up tandem mass spectrometry due to their highly similar sequences. For many, the only tractable approach is with intact protein top-down tandem mass spectrometry. Here, ultra-high-resolution FT-ICR MS and MS/MS yield quantitative relative abundances of all detected HeLa H2A and H2B isobaric and isomeric variants with a label-free approach. We extend the analysis to identify and relatively quantitate 16 proteoforms from 12 sequence variants of histone H2A and 10 proteoforms of histone H2B from three other cell lines: human embryonic stem cells (WA09), U937, and a prostate cancer cell line LaZ. The top-down MS/MS approach provides a path forward for more extensive elucidation of the biological role of many previously unstudied histone variants and post-translational modifications.

BAF250a Protein Regulates Nucleosome Occupancy and Histone Modifications in Priming Embryonic Stem Cell Differentiation.

The unique chromatin signature of ES cells is fundamental to the pluripotency and differentiation of ES cells. One key feature is the poised chromatin state of master developmental genes that are transcriptionally repressed in ES cells but ready to be activated in response to differentiation signals. Poised chromatin in ES cells contains both H3 Lys-4 trimethylation (H3K4me3) and H3 Lys-27 trimethylation (H3K27me3) methylation, indicating activating and repressing potential. However, the contribution of non-covalent chromatin structure to the poised state is not well understood. To address whether remodeling of nucleosomes is important to the poised state, we characterized the function of BAF250a, a key regulatory subunit of the ES cell ATP-dependent Brahma-associated factor (BAF) chromatin remodeling complex (esBAF). Acute deletion of BAF250a disrupted the differentiation potential of ES cells by altering the expression timing of key developmental genes and pluripotent genes. Our genome-wide nucleosome and histone modification analyses indicated that the disruption of gene expression timing was largely due to changes of chromatin structures at poised genes, particularly those key developmental genes mediated by BAF250a. Specifically, BAF250a deletion caused a nucleosome occupancy increase at H3K4me3- and/or H3K27me3-associated promoters. Moreover, H3K27me3 levels and the number of bivalent promoter genes were reduced in BAF250a KO ES cells. We revealed that BAF250a ablation led to elevated Brg1 but reduced Suz12 recruitment at nucleosome occupancy-increased regions, indicating an unexpected and complicated role of BAF250a in regulating esBAF and Polycomb repressive complex (PRC) activities. Together, our studies identified that BAF250a mediates esBAF and PRC functions to establish the poised chromatin configuration in ES cells, which is essential for the proper differentiation of ES cells.

Chromatin-interaction compartment switch at developmentally regulated chromosomal domains reveals an unusual principle of chromatin folding.

Several 400- to 800-kb murine chromosome domains switch from early to late replication during loss of pluripotency, accompanied by a stable form of gene silencing that is resistant to reprogramming. We found that, whereas enhanced nuclease accessibility correlated with early replication genome-wide, domains that switch replication timing during differentiation were exceptionally inaccessible even when early-replicating. Nonetheless, two domains studied in detail exhibited substantial changes in transcriptional activity and higher-order chromatin unfolding confined to the region of replication timing change. Chromosome conformation capture (4C) data revealed that in the unfolded state in embryonic stem cells, these domains interacted preferentially with the early-replicating chromatin compartment, rarely interacting even with flanking late-replicating domains, whereas after differentiation, these same domains preferentially associated with late-replicating chromatin, including flanking domains. In both configurations they retained local boundaries of self-interaction, supporting the replication domain model of replication-timing regulation. Our results reveal a principle of developmentally regulated, large-scale chromosome folding involving a subnuclear compartment switch of inaccessible chromatin. This unusual level of regulation may underlie resistance to reprogramming in replication-timing switch regions.

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.

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 initiation zone, and suggest that initiation zones may be regions of generally low nucleosome occupancy where flexible nucleosome positioning permits flexible pre-RC assembly sites.

Virus detection via programmable Type III-A CRISPR-Cas systems.

Among the currently available virus detection assays, those based on the programmable CRISPR-Cas enzymes have the advantage of rapid reporting and high sensitivity without the requirement of thermocyclers. Type III-A CRISPR-Cas system is a multi-component and multipronged immune effector, activated by viral RNA that previously has not been repurposed for disease detection owing in part to the complex enzyme reconstitution process and functionality. Here, we describe the construction and application of a virus detection method, based on an in vivo-reconstituted Type III-A CRISPR-Cas system. This system harnesses both RNA- and transcription-activated dual nucleic acid cleavage activities as well as internal signal amplification that allow virus detection with high sensitivity and at multiple settings. We demonstrate the use of the Type III-A system-based method in detection of SARS-CoV-2 that reached 2000 copies/µl sensitivity in amplification-free and 60 copies/µl sensitivity via isothermal amplification within 30 min and diagnosed SARS-CoV-2-infected patients in both settings. The high sensitivity, flexible reaction conditions, and the small molecular-driven amplification make the Type III-A system a potentially unique nucleic acid detection method with broad applications.

Multiple roles of H2A.Z in regulating promoter chromatin architecture in human cells.

Chromatin accessibility of a promoter is fundamental in regulating transcriptional activity. The histone variant H2A.Z has been shown to contribute to this regulation, but its role has remained poorly understood. Here, we prepare high-depth maps of the position and accessibility of H2A.Z-containing nucleosomes for all human Pol II promoters in epithelial, mesenchymal and isogenic cancer cell lines. We find that, in contrast to the prevailing model, many different types of active and inactive promoter structures are observed that differ in their nucleosome organization and sensitivity to MNase digestion. Key aspects of an active chromatin structure include positioned H2A.Z MNase resistant nucleosomes upstream or downstream of the TSS, and a MNase sensitive nucleosome at the TSS. Furthermore, the loss of H2A.Z leads to a dramatic increase in the accessibility of transcription factor binding sites. Collectively, these results suggest that H2A.Z has multiple and distinct roles in regulating gene expression dependent upon its location in a promoter.

Chromatin structure profile data from DNS-seq: Differential nuclease sensitivity mapping of four reference tissues of B73 maize (Zea mays L).

Presented here are data from Next-Generation Sequencing of differential micrococcal nuclease digestions of formaldehyde-crosslinked chromatin in selected tissues of maize (Zea mays) inbred line B73. Supplemental materials include a wet-bench protocol for making DNS-seq libraries, the DNS-seq data processing pipeline for producing genome browser tracks. This report also includes the peak-calling pipeline using the iSeg algorithm to segment positive and negative peaks from the DNS-seq difference profiles. The data repository for the sequence data is the NCBI SRA, BioProject Accession PRJNA445708.

Expression of the Brn-3b transcription factor correlates with expression of HSP-27 in breast cancer biopsies and is required for maximal activation of the HSP-27 promoter.

In breast cancer, overexpression of the small heat shock protein, HSP-27, is associated with increased anchorage-independent growth, increased invasiveness, and resistance to chemotherapeutic drugs and is associated with poor prognosis and reduced disease-free survival. Therefore, factors that increase the expression of HSP-27 in breast cancer are likely to affect the prognosis and outcome of treatment. In this study, we show a strong correlation between elevated levels of the Brn-3b POU transcription factor and high levels of HSP-27 protein in manipulated MCF-7 breast cancer cells as well as in human breast biopsies. Conversely, HSP-27 is decreased on loss of Brn-3b. In cotransfection assays, Brn-3b can strongly transactivate the HSP-27 promoter, supporting a role for direct regulation of HSP-27 expression. Brn-3b also cooperates with the estrogen receptor (ER) to facilitate maximal stimulation of the HSP-27 promoter, with significantly enhanced activity of this promoter observed on coexpression of Brn-3b and ER compared with either alone. RNA interference and site-directed mutagenesis support the requirement for the Brn-3b binding site on the HSP-27 promoter, which facilitates maximal transactivation either alone or on interaction with the ER. Chromatin immunoprecipitation provides evidence for association of Brn-3b with the HSP-27 promoter in the intact cell. Thus, Brn-3b can, directly and indirectly (via interaction with the ER), activate HSP-27 expression, and this may represent one mechanism by which Brn-3b mediates its effects in breast cancer cells.

Hierarchical regulation of the genome: global changes in nucleosome organization potentiate genome response.

Nucleosome occupancy is critically important in regulating access to the eukaryotic genome. Few studies in human cells have measured genome-wide nucleosome distributions at high temporal resolution during a response to a common stimulus. We measured nucleosome distributions at high temporal resolution following Kaposi's-sarcoma-associated herpesvirus (KSHV) reactivation using our newly developed mTSS-seq technology, which maps nucleosome distribution at the transcription start sites (TSS) of all human genes. Nucleosomes underwent widespread changes in organization 24 hours after KSHV reactivation and returned to their basal nucleosomal architecture 48 hours after KSHV reactivation. The widespread changes consisted of an indiscriminate remodeling event resulting in the loss of nucleosome rotational phasing signals. Additionally, one in six TSSs in the human genome possessed nucleosomes that are translationally remodeled. 72% of the loci with translationally remodeled nucleosomes have nucleosomes that moved to positions encoded by the underlying DNA sequence. Finally we demonstrated that these widespread alterations in nucleosomal architecture potentiated regulatory factor binding. These descriptions of nucleosomal architecture changes provide a new framework for understanding the role of chromatin in the genomic response, and have allowed us to propose a hierarchical model for chromatin-based regulation of genome response.

Comprehensive nucleosome mapping of the human genome in cancer progression.

Altered chromatin structure is a hallmark of cancer, and inappropriate regulation of chromatin structure may represent the origin of transformation. Important studies have mapped human nucleosome distributions genome wide, but the role of chromatin structure in cancer progression has not been addressed. We developed a MNase-Transcription Start Site Sequence Capture method (mTSS-seq) to map the nucleosome distribution at human transcription start sites genome-wide in primary human lung and colon adenocarcinoma tissue. Here, we confirm that nucleosome redistribution is an early, widespread event in lung (LAC) and colon (CRC) adenocarcinoma. These altered nucleosome architectures are consistent between LAC and CRC patient samples indicating that they may serve as important early adenocarcinoma markers. We demonstrate that the nucleosome alterations are driven by the underlying DNA sequence and potentiate transcription factor binding. We conclude that DNA-directed nucleosome redistributions are widespread early in cancer progression. We have proposed an entirely new hierarchical model for chromatin-mediated genome regulation.

Stimulation of the Drosophila immune system alters genome-wide nucleosome occupancy.

In eukaryotes, nucleosomes participate in all DNA-templated events by regulating access to the underlying DNA sequence. However, nucleosome dynamics during a genome response have not been well characterized [1,2]. We stimulated Drosophila S2 cells with heat-killed Gram-negative bacteria Salmonella typhimurium, and mapped genome-wide nucleosome occupancy at high temporal resolution by MNase-seq using Illumina HiSeq 2500. We show widespread nucleosome occupancy change in S2 cells during the immune response, with the significant nucleosomal loss occurring at 4 h after stimulation. Data have been deposited to the Gene Expression Omnibus (GEO) database repository with the dataset identifier GSE64507.

Nucleosome Repositioning: A Novel Mechanism for Nicotine- and Cocaine-Induced Epigenetic Changes.

Drugs of abuse modify behavior by altering gene expression in the brain. Gene expression can be regulated by changes in DNA methylation as well as by histone modifications, which alter chromatin structure, DNA compaction and DNA accessibility. In order to better understand the molecular mechanisms directing drug-induced changes in chromatin structure, we examined DNA-nucleosome interactions within promoter regions of 858 genes in human neuroblastoma cells (SH-SY5Y) exposed to nicotine or cocaine. Widespread, drug- and time-resolved repositioning of nucleosomes was identified at the transcription start site and promoter region of multiple genes. Nicotine and cocaine produced unique and shared changes in terms of the numbers and types of genes affected, as well as repositioning of nucleosomes at sites which could increase or decrease the probability of gene expression based on DNA accessibility. Half of the drug-induced nucleosome positions approximated a theoretical model of nucleosome occupancy based on physical and chemical characteristics of the DNA sequence, whereas the basal or drug naïve positions were generally DNA sequence independent. Thus we suggest that nucleosome repositioning represents an initial dynamic genome-wide alteration of the transcriptional landscape preceding more selective downstream transcriptional reprogramming, which ultimately characterizes the cell- and tissue-specific responses to drugs of abuse.

DNA-Encoded Chromatin Structural Intron Boundary Signals Identify Conserved Genes with Common Function.

The regulation of metazoan gene expression occurs in part by pre-mRNA splicing into mature RNAs. Signals affecting the efficiency and specificity with which introns are removed have not been completely elucidated. Splicing likely occurs cotranscriptionally, with chromatin structure playing a key regulatory role. We calculated DNA encoded nucleosome occupancy likelihood (NOL) scores at the boundaries between introns and exons across five metazoan species. We found that (i) NOL scores reveal a sequence-based feature at the introns on both sides of the intron-exon boundary; (ii) this feature is not part of any recognizable consensus sequence; (iii) this feature is conserved throughout metazoa; (iv) this feature is enriched in genes sharing similar functions: ATPase activity, ATP binding, helicase activity, and motor activity; (v) genes with these functions exhibit different genomic characteristics; (vi) in vivo nucleosome positioning data confirm ontological enrichment at this feature; and (vii) genes with this feature exhibit unique dinucleotide distributions at the intron-exon boundary. The NOL scores point toward a physical property of DNA that may play a role in the mechanism of pre-mRNA splicing. These results provide a foundation for identification of a new set of regulatory DNA elements involved in splicing regulation.