The impact of nucleosome structure on CRISPR/Cas9 fidelity.

The clustered regularly interspaced short palindromic repeats (CRISPR) Cas system is a powerful tool that has the potential to become a therapeutic gene editor in the near future. Cas9 is the best studied CRISPR system and has been shown to have problems that restrict its use in therapeutic applications. Chromatin structure is a known impactor of Cas9 targeting and there is a gap in knowledge on Cas9's efficacy when targeting such locations. To quantify at a single base pair resolution how chromatin inhibits on-target gene editing relative to off-target editing of exposed mismatching targets, we developed the gene editor mismatch nucleosome inhibition assay (GEMiNI-seq). GEMiNI-seq utilizes a library of nucleosome sequences to examine all target locations throughout nucleosomes in a single assay. The results from GEMiNI-seq revealed that the location of the protospacer-adjacent motif (PAM) sequence on the nucleosome edge drives the ability for Cas9 to access its target sequence. In addition, Cas9 had a higher affinity for exposed mismatched targets than on-target sequences within a nucleosome. Overall, our results show how chromatin structure impacts the fidelity of Cas9 to potential targets and highlight how targeting sequences with exposed PAMs could limit off-target gene editing, with such considerations improving Cas9 efficacy and resolving current limitations.

Defining TP53 pioneering capabilities with competitive nucleosome binding assays.

Accurate gene expression requires the targeting of transcription factors (TFs) to regulatory sequences often occluded within nucleosomes. The ability to target a TF binding site (TFBS) within a nucleosome has been the defining characteristic for a special class of TFs known as pioneer factors. Recent studies suggest TP53 functions as a pioneer factor that can target its TFBS within nucleosomes, but it remains unclear how TP53 binds to nucleosomal DNA. To comprehensively examine TP53 nucleosome binding, we competitively bound TP53 to multiple in vitro-formed nucleosomes containing a high- or low-affinity TP53 TFBS located at differing translational and rotational positions within the nucleosome. Stable TP53-nucleosome complexes were isolated and quantified using next-generation sequencing. Our results demonstrate TP53 binding is limited to nucleosome edges with significant binding inhibition occurring within 50 bp of the nucleosome dyad. Binding site affinity only affects TP53 binding for TFBSs located at the same nucleosomal positions; otherwise, nucleosome position takes precedence. Furthermore, TP53 has strong nonspecific nucleosome binding facilitating its interaction with chromatin. Our in vitro findings were confirmed by examining TP53-induced binding in a cell line model, showing induced binding at nucleosome edges flanked by a nucleosome-free region. Overall, our results suggest that the pioneering capabilities of TP53 are driven by nonspecific nucleosome binding with specific binding at nucleosome edges.

Molecular dissection of the oncogenic role of ETS1 in the mesenchymal subtypes of head and neck squamous cell carcinoma.

Head and Neck Squamous Cell Carcinoma (HNSCC) is a heterogeneous disease of significant mortality and with limited treatment options. Recent genomic analysis of HNSCC tumors has identified several distinct molecular classes, of which the mesenchymal subtype is associated with Epithelial to Mesenchymal Transition (EMT) and shown to correlate with poor survival and drug resistance. Here, we utilize an integrated approach to characterize the molecular function of ETS1, an oncogenic transcription factor specifically enriched in Mesenchymal tumors. To identify the global ETS1 cistrome, we have performed integrated analysis of RNA-Seq, ChIP-Seq and epigenomic datasets in SCC25, a representative ETS1high mesenchymal HNSCC cell line. Our studies implicate ETS1 as a crucial regulator of broader oncogenic processes and specifically Mesenchymal phenotypes, such as EMT and cellular invasion. We found that ETS1 preferentially binds cancer specific regulator elements, in particular Super-Enhancers of key EMT genes, highlighting its role as a master regulator. Finally, we show evidence that ETS1 plays a crucial role in regulating the TGF-ß pathway in Mesenchymal cell lines and in leading-edge cells in primary HNSCC tumors that are endowed with partial-EMT features. Collectively our study highlights ETS1 as a key regulator of TGF-ß associated EMT and reveals new avenues for sub-type specific therapeutic intervention.

Pioneer factors and their in vitro identification methods.

Pioneer transcription factors are a special group of transcription factors that can interact with nucleosomal DNA and initiate regulatory events. Their binding to regulatory regions is the first event in gene activation and can occur in silent or heterochromatin regions. Several research groups have endeavored to define pioneer factors and study their binding characteristics using various techniques. In this review, we describe the in vitro methods used to define and characterize pioneer factors, paying particular attention to differences in methodologies and how these differences can affect results.

Evolutionary re-wiring of p63 and the epigenomic regulatory landscape in keratinocytes and its potential implications on species-specific gene expression and phenotypes.

Although epidermal keratinocyte development and differentiation proceeds in similar fashion between humans and mice, evolutionary pressures have also wrought significant species-specific physiological differences. These differences between species could arise in part, by the rewiring of regulatory network due to changes in the global targets of lineage-specific transcriptional master regulators such as p63. Here we have performed a systematic and comparative analysis of the p63 target gene network within the integrated framework of the transcriptomic and epigenomic landscape of mouse and human keratinocytes. We determined that there exists a core set of ∼1600 genomic regions distributed among enhancers and super-enhancers, which are conserved and occupied by p63 in keratinocytes from both species. Notably, these DNA segments are typified by consensus p63 binding motifs under purifying selection and are associated with genes involved in key keratinocyte and skin-centric biological processes. However, the majority of the p63-bound mouse target regions consist of either murine-specific DNA elements that are not alignable to the human genome or exhibit no p63 binding in the orthologous syntenic regions, typifying an occupancy lost subset. Our results suggest that these evolutionarily divergent regions have undergone significant turnover of p63 binding sites and are associated with an underlying inactive and inaccessible chromatin state, indicative of their selective functional activity in the transcriptional regulatory network in mouse but not human. Furthermore, we demonstrate that this selective targeting of genes by p63 correlates with subtle, but measurable transcriptional differences in mouse and human keratinocytes that converges on major metabolic processes, which often exhibit species-specific trends. Collectively our study offers possible molecular explanation for the observable phenotypic differences between the mouse and human skin and broadly informs on the prevailing principles that govern the tug-of-war between evolutionary forces of rigidity and plasticity over transcriptional regulatory programs.

Trypanosome RNA Editing Mediator Complex proteins have distinct functions in gRNA utilization.

Uridine insertion/deletion RNA editing is an essential process in kinetoplastid parasites whereby mitochondrial mRNAs are modified through the specific insertion and deletion of uridines to generate functional open reading frames, many of which encode components of the mitochondrial respiratory chain. The roles of numerous non-enzymatic editing factors have remained opaque given the limitations of conventional methods to interrogate the order and mechanism by which editing progresses and thus roles of individual proteins. Here, we examined whole populations of partially edited sequences using high throughput sequencing and a novel bioinformatic platform, the Trypanosome RNA Editing Alignment Tool (TREAT), to elucidate the roles of three proteins in the RNA Editing Mediator Complex (REMC). We determined that the factors examined function in the progression of editing through a gRNA; however, they have distinct roles and REMC is likely heterogeneous in composition. We provide the first evidence that editing can proceed through numerous paths within a single gRNA and that non-linear modifications are essential, generating commonly observed junction regions. Our data support a model in which RNA editing is executed via multiple paths that necessitate successive re-modification of junction regions facilitated, in part, by the REMC variant containing TbRGG2 and MRB8180.

Composition and diversity of the subgingival microbiome and its relationship with age in postmenopausal women: an epidemiologic investigation.

BACKGROUND: The extent to which the composition and diversity of the oral microbiome varies with age is not clearly understood. METHODS: The 16S rRNA gene of subgingival plaque in 1219 women, aged 53-81 years, was sequenced and its taxonomy annotated against the Human Oral Microbiome Database (v.14.5). Composition of the subgingival microbiome was described in terms of centered log(2)-ratio (CLR) transformed OTU values, relative abundance, and prevalence. Correlations between microbiota abundance and age were evelauted using Pearson Product Moment correlations. P-values were corrected for multiple testing using the Bonferroni method. RESULTS: Of the 267 species identified overall, Veillonella dispar was the most abundant bacteria when described by CLR OTU (mean 8.3) or relative abundance (mean 8.9%); whereas Streptococcus oralis, Veillonella dispar and Veillonella parvula were most prevalent (100%, all) when described as being present at any amount. Linear correlations between age and several CLR OTUs (Pearson r = - 0.18 to 0.18), of which 82 (31%) achieved statistical significance (P < 0.05). The correlations lost significance following Bonferroni correction. Twelve species that differed across age groups (each corrected P < 0.05); 5 (42%) were higher in women ages 50-59 compared to ≥70 (corrected P < 0.05), and 7 (48%) were higher in women 70 years and older. CONCLUSIONS: We identified associations between several bacterial species and age across the age range of postmenopausal women studied. Understanding the functions of these bacteria could identify intervention targets to enhance oral health in later life.

Suppressed hepatic bile acid signalling despite elevated production of primary and secondary bile acids in NAFLD.

OBJECTIVE: Bile acids are regulators of lipid and glucose metabolism, and modulate inflammation in the liver and other tissues. Primary bile acids such as cholic acid and chenodeoxycholic acid (CDCA) are produced in the liver, and converted into secondary bile acids such as deoxycholic acid (DCA) and lithocholic acid by gut microbiota. Here we investigated the possible roles of bile acids in non-alcoholic fatty liver disease (NAFLD) pathogenesis and the impact of the gut microbiome on bile acid signalling in NAFLD. DESIGN: Serum bile acid levels and fibroblast growth factor 19 (FGF19), liver gene expression profiles and gut microbiome compositions were determined in patients with NAFLD, high-fat diet-fed rats and their controls. RESULTS: Serum concentrations of primary and secondary bile acids were increased in patients with NAFLD. In per cent, the farnesoid X receptor (FXR) antagonistic DCA was increased, while the agonistic CDCA was decreased in NAFLD. Increased mRNA expression for cytochrome P450 7A1, Na+-taurocholate cotransporting polypeptide and paraoxonase 1, no change in mRNA expression for small heterodimer partner and bile salt export pump, and reduced serum FGF19 were evidence of impaired FXR and fibroblast growth factor receptor 4 (FGFR4)-mediated signalling in NAFLD. Taurine and glycine metabolising bacteria were increased in the gut of patients with NAFLD, reflecting increased secondary bile acid production. Similar changes in liver gene expression and the gut microbiome were observed in high-fat diet-fed rats. CONCLUSIONS: The serum bile acid profile, the hepatic gene expression pattern and the gut microbiome composition consistently support an elevated bile acid production in NAFLD. The increased proportion of FXR antagonistic bile acid explains, at least in part, the suppression of hepatic FXR-mediated and FGFR4-mediated signalling. Our study suggests that future NAFLD intervention may target the components of FXR signalling, including the bile acid converting gut microbiome.

Gut microbiome may contribute to insulin resistance and systemic inflammation in obese rodents: a meta-analysis.

A number of studies have associated obesity with altered gut microbiota, although results are discordant regarding compositional changes in the gut microbiota of obese animals. Herein we used a meta-analysis to obtain an unbiased evaluation of structural and functional changes of the gut microbiota in diet-induced obese rodents. The raw sequencing data of nine studies generated from high-fat diet (HFD)-induced obese rodent models were processed with QIIME to obtain gut microbiota compositions. Biological functions were predicted and annotated with KEGG pathways with PICRUSt. No significant difference was observed for alpha diversity and Bacteroidetes-to-Firmicutes ratio between obese and lean rodents. Bacteroidia, Clostridia, Bacilli, and Erysipelotrichi were dominant classes, but gut microbiota compositions varied among studies. Meta-analysis of the nine microbiome data sets identified 15 differential taxa and 57 differential pathways between obese and lean rodents. In obese rodents, increased abundance was observed for Dorea, Oscillospira, and Ruminococcus, known for fermenting polysaccharide into short chain fatty acids (SCFAs). Decreased Turicibacter and increased Lactococcus are consistent with elevated inflammation in the obese status. Differential functional pathways of the gut microbiome in obese rodents included enriched pyruvate metabolism, butanoate metabolism, propanoate metabolism, pentose phosphate pathway, fatty acid biosynthesis, and glycerolipid metabolism pathways. These pathways converge in the function of carbohydrate metabolism, SCFA metabolism, and biosynthesis of lipid. HFD-induced obesity results in structural and functional dysbiosis of gut microbiota. The altered gut microbiome may contribute to obesity development by promoting insulin resistance and systemic inflammation.

High-throughput sequencing of partially edited trypanosome mRNAs reveals barriers to editing progression and evidence for alternative editing.

Uridine insertion/deletion RNA editing in kinetoplastids entails the addition and deletion of uridine residues throughout the length of mitochondrial transcripts to generate translatable mRNAs. This complex process requires the coordinated use of several multiprotein complexes as well as the sequential use of noncoding template RNAs called guide RNAs. The majority of steady-state mitochondrial mRNAs are partially edited and often contain regions of mis-editing, termed junctions, whose role is unclear. Here, we report a novel method for sequencing entire populations of pre-edited partially edited, and fully edited RNAs and analyzing editing characteristics across populations using a new bioinformatics tool, the Trypanosome RNA Editing Alignment Tool (TREAT). Using TREAT, we examined populations of two transcripts, RPS12 and ND7-5', in wild-typeTrypanosoma brucei We provide evidence that the majority of partially edited sequences contain junctions, that intrinsic pause sites arise during the progression of editing, and that the mechanisms that mediate pausing in the generation of canonical fully edited sequences are distinct from those that mediate the ends of junction regions. Furthermore, we identify alternatively edited sequences that constitute plausible alternative open reading frames and identify substantial variability in the 5' UTRs of both canonical and alternatively edited sequences. This work is the first to use high-throughput sequencing to examine full-length sequences of whole populations of partially edited transcripts. Our method is highly applicable to current questions in the RNA editing field, including defining mechanisms of action for editing factors and identifying potential alternatively edited sequences.

A global analysis of the complex landscape of isoforms and regulatory networks of p63 in human cells and tissues.

BACKGROUND: The transcription factor p63 belongs to the p53/p63/p73 family and plays key functional roles during normal epithelial development and differentiation and in pathological states such as squamous cell carcinomas. The human TP63 gene, located on chromosome 3q28 is driven by two promoters that generate the full-length transactivating (TA) and N-terminal truncated (ΔN) isoforms. Furthermore alternative splicing at the C-terminus gives rise to additional α, ß, γ and likely several other minor variants. Teasing out the expression and biological function of each p63 variant has been both the focus of, and a cause for contention in the p63 field. RESULTS: Here we have taken advantage of a burgeoning RNA-Seq based genomic data-sets to examine the global expression profiles of p63 isoforms across commonly utilized human cell-lines and major tissues and organs. Consistent with earlier studies, we find ΔNp63 transcripts, primarily that of the ΔNp63α isoforms, to be expressed in most cells of epithelial origin such as those of skin and oral tissues, mammary glands and squamous cell carcinomas. In contrast, TAp63 is not expressed in the majority of normal cell-types and tissues; rather it is selectively expressed at moderate to high levels in a subset of Burkitt's and diffuse large B-cell lymphoma cell lines. We verify this differential expression pattern of p63 isoforms by Western blot analysis, using newly developed ΔN and TA specific antibodies. Furthermore using unsupervised clustering of human cell lines, tissues and organs, we show that ΔNp63 and TAp63 driven transcriptional networks involve very distinct sets of molecular players, which may underlie their different biological functions. CONCLUSIONS: In this study we report comprehensive and global expression profiles of p63 isoforms and their relationship to p53/p73 and other potential transcriptional co-regulators. We curate publicly available data generated in part by consortiums such as ENCODE, FANTOM and Human Protein Atlas to delineate the vastly different transcriptomic landscapes of ΔNp63 and TAp63. Our studies help not only in dispelling prevailing myths and controversies on p63 expression in commonly used human cell lines but also augur new isoform- and cell type-specific activities of p63.

Iron-responsive chromatin remodelling and MAPK signalling enhance adhesion in Candida albicans.

Recent cumulative data show that various transcription factors are recruited to the chromatin in an iron-responsive manner to affect diverse cellular functions in the pathogenic fungus Candida albicans. Here we identified groups of iron-responsive genes in C. albicans by chromatin remodelling analysis at gene promoters, using micrococcal nuclease (MNase) digestion followed by deep sequencing. Chromatin in the promoter regions of iron uptake and utilization genes showed repressed and active configuration, respectively, under iron-replete conditions. GO Term enrichment analysis of genes with differentially remodelled chromatin, in respective promoter locales, suggested that many genes involved in adhesion are also iron-responsive. C. albicans was observed to be more self-adherent (twofold increase) and formed higher biofilm mass (77% increase) in the presence of iron. Furthermore, we identified various known and novel adhesion-related genes with iron-dependent active chromatin profiles that are indicative of potential upregulation under iron-replete conditions. Transcription factor Cph1 that is activated upon Cek1 phosphorylation also showed an active chromatin profile under iron-replete conditions and cells showed iron-responsive Cek1 MAPK phosphorylation in the presence of iron. Thus, iron affects diverse biological functions by modulating chromatin profiles of large gene sets and by signalling through Cek1 MAPK in C. albicans.

A chromatin-mediated mechanism for specification of conditional transcription factor targets.

Organisms respond to changes in their environment, and many such responses are initiated at the level of gene transcription. Here, we provide evidence for a previously undiscovered mechanism for directing transcriptional regulators to new binding targets in response to an environmental change. We show that repressor-activator protein 1 (Rap1), a master regulator of yeast metabolism, binds to an expanded target set after glucose depletion despite decreasing protein levels and no evidence of posttranslational modification. Computational analysis predicts that proteins capable of recruiting the chromatin regulator Tup1 act to restrict the binding distribution of Rap1 in the presence of glucose. Deletion of the gene(s) encoding Tup1, recruiters of Tup1 or chromatin regulators recruited by Tup1 cause Rap1 to bind specifically and inappropriately to low-glucose targets. These data, combined with whole-genome measurements of nucleosome occupancy and Tup1 distribution, provide evidence for a mechanism of dynamic target specification that coordinates the genome-wide distribution of intermediate-affinity DNA sequence motifs with chromatin-mediated regulation of accessibility to those sites.

Role of chromatin and transcriptional co-regulators in mediating p63-genome interactions in keratinocytes.

BACKGROUND: The Transcription Factor (TF) p63 is a master regulator of epidermal development and differentiation as evident from the remarkable skin phenotype of p63 mouse knockouts. Furthermore, ectopic expression of p63 alone is sufficient to convert simple epithelium into stratified epithelial tissues in vivo and p63 is required for efficient transdifferentiation of fibroblasts into keratinocytes. However, little is known about the molecular mechanisms of p63 function, in particular how it selects its target sites in the genome. p63, which acts both as an activator and repressor of transcription, recognizes a canonical binding motif that occurs over 1 million times in the human genome. But, in human keratinocytes less than 12,000 of these sites are bound in vivo suggesting that underlying chromatin architecture and cooperating TFs mediate p63-genome interactions. RESULTS: We find that the chromatin architecture at p63-bound targets possess distinctive features and can be used to categorize p63 targets into proximal promoters (1%), enhancers (59%) and repressed or inactive (40%) regulatory elements. Our analysis shows that the chromatin modifications H3K4me1, H3K27me3, along with overall chromatin accessibility status can accurately predict bonafide p63-bound sites without a priori DNA sequence information. Interestingly, however there exists a qualitative correlation between the p63 binding motif and accessibility and H3K4me1 levels. Furthermore, we use a comprehensive in silico approach that leverages ENCODE data to identify several known TFs such as AP1, AP2 and novel TFs (RFX5 for e.g.) that can potentially cooperate with p63 to modulate its myriad biological functions in keratinocytes. CONCLUSIONS: Our analysis shows that p63 bound genomic locations in keratinocytes are accessible, marked by active histone modifications, and co-targeted by other developmentally important transcriptional regulators. Collectively, our results suggest that p63 might actively remodel and/or influence chromatin dynamics at its target sites and in the process dictate its own DNA binding and possibly that of adjacent TFs.

An integrative approach to understanding the combinatorial histone code at functional elements.

SUMMARY: The rapid advancement of genomic technology has revealed the enormous complexity and combinatorial nature of chromatin modifications. To facilitate interpretation of the combinatorial nature of chromatin, we have developed a novel method to integrate all chromatin datasets into distinct nucleosome types (nucleosome alphabet). We have applied this approach to Saccharomyces cerevisiae, generating a nucleosome alphabet, which forms chromatin motifs when mapped back to the genome. By applying novel chromatin alignment and global word search approaches, we have defined distinctive chromatin motifs for introns, origins of replication, tRNAs, antisense transcripts, double-strand-break hotspots and DNase hypersensitive sites, and can distinguish genes by expression level. We have also uncovered strong associations between transcription factor binding and specific types of nucleosomes. Our results demonstrate the uses and functionality of defining a chromatin alphabet and provide a unique and novel framework for exploring chromatin architecture. CONTACT: mjbuck@buffalo.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Chromatin architectures at fission yeast transcriptional promoters and replication origins.

We have used micrococcal nuclease (MNase) digestion followed by deep sequencing in order to obtain a higher resolution map than previously available of nucleosome positions in the fission yeast, Schizosaccharomyces pombe. Our data confirm an unusually short average nucleosome repeat length, ∼152 bp, in fission yeast and that transcriptional start sites (TSSs) are associated with nucleosome-depleted regions (NDRs), ordered nucleosome arrays downstream and less regularly spaced upstream nucleosomes. In addition, we found enrichments for associated function in four of eight groups of genes clustered according to chromatin configurations near TSSs. At replication origins, our data revealed asymmetric localization of pre-replication complex (pre-RC) proteins within large NDRs-a feature that is conserved in fission and budding yeast and is therefore likely to be conserved in other eukaryotic organisms.

Nucleosome-binding by TP53, TP63, and TP73 is determined by the composition, accessibility, and helical orientation of their binding sites.

The p53 family of transcription factors plays key roles in driving development and combating cancer by regulating gene expression. TP53, TP63, and TP73-the three members of the p53 family-regulate gene expression by binding to their DNA binding sites, many of which are situated within nucleosomes. To thoroughly examine the nucleosome-binding abilities of the p53 family, we used Pioneer-seq, a technique that assesses a transcription factor's binding affinity to its DNA binding sites at all possible positions within the nucleosome core particle. Using Pioneer-seq, we analyzed the binding affinity of TP53, TP63, and TP73 to 10 p53-family binding sites across the nucleosome core particle. We found that the affinity of TP53, TP63, and TP73 for nucleosomes was largely determined by the positioning of p53-family binding sites within nucleosomes; p53-family members bind strongly to the more accessible edges of nucleosomes but weakly to the less accessible centers of nucleosomes. We also found that the DNA-helical orientation of p53-family binding sites within nucleosomal DNA impacted the nucleosome-binding affinity of p53-family members. The composition of their binding sites also impacted each p53-family member's nucleosome-binding affinities only when the binding site was located in an accessible location. Taken together, our results show that the accessibility, composition, and helical orientation of p53-family binding sites collectively determine the nucleosome-binding affinities of TP53, TP63, and TP73. These findings help explain the rules underlying p53-family-nucleosome binding and thus provide requisite insight into how we may better control gene-expression changes involved in development and tumor suppression.

Defining Porphyromonas gingivalis strains associated with periodontal disease.

Porphyromonas gingivalis, a Gram-negative anaerobic bacterium commonly found in human subgingival plaque, is a major etiologic agent for periodontitis and has been associated with multiple systemic pathologies. Many P. gingivalis strains have been identified and different strains possess different virulence factors. Current oral microbiome approaches (16S or shotgun) have been unable to differentiate P. gingivalis strains. This study presents a new approach that aims to improve the accuracy of strain identification, using a detection method based on sequencing of the intergenic spacer region (ISR) which is variable between P. gingivalis strains. Our approach uses two-step PCR to amplify only the P. gingivalis ISR region. Samples are then sequenced with an Illumina sequencer and mapped to specific strains. Our approach was validated by examining subgingival plaque from 153 participants with and without periodontal disease. We identified the avirulent strain ATCC33277/381 as the most abundant strain across all sample types. The W83/W50 strain was significantly enriched in periodontitis, with 13% of participants harboring that strain. Overall, this approach can have significant implications not only for the diagnosis and treatment of periodontal disease but also for other diseases where P. gingivalis or its toxins have been implicated, such as Alzheimer's disease.

Discovering chromatin motifs using FAIRE sequencing and the human diploid genome.

BACKGROUND: Specific chromatin structures are associated with active or inactive gene transcription. The gene regulatory elements are intrinsically dynamic and alternate between inactive and active states through the recruitment of DNA binding proteins, such as chromatin-remodeling proteins. RESULTS: We developed a unique genome-wide method to discover DNA motifs associated with chromatin accessibility using formaldehyde-assisted isolation of regulatory elements with high-throughput sequencing (FAIRE-seq). We aligned the FAIRE-seq reads to the GM12878 diploid genome and subsequently identified differential chromatin-state regions (DCSRs) using heterozygous SNPs. The DCSR pairs represent the locations of imbalances of chromatin accessibility between alleles and are ideal to reveal chromatin motifs that may directly modulate chromatin accessibility. In this study, we used DNA 6-10mer sequences to interrogate all DCSRs, and subsequently discovered conserved chromatin motifs with significant changes in the occurrence frequency. To investigate their likely roles in biology, we studied the annotated protein associated with each of the top ten chromatin motifs genome-wide, in the intergenic regions and in genes, respectively. As a result, we found that most of these annotated motifs are associated with chromatin remodeling, reflecting their significance in biology. CONCLUSIONS: Our method is the first one using fully phased diploid genome and FAIRE-seq to discover motifs associated with chromatin accessibility. Our results were collected to construct the first chromatin motif database (CMD), providing the potential DNA motifs recognized by chromatin-remodeling proteins and is freely available at http://syslab.nchu.edu.tw/chromatin.

ArchTEx: accurate extraction and visualization of next-generation sequence data.

MOTIVATION: The extension of mapped sequence tags is a common step in the analysis of single-end next-generation sequencing (NGS) data from protein localization and chromatin studies. The optimal extension can vary depending on experimental and technical conditions. Improper extension of sequence tags can obscure or mislead the interpretation of NGS results. We present an algorithm, ArchTEx (Architectural Tag Extender), which identifies the optimal extension of sequence tags based on the maximum correlation between forward and reverse tags and extracts and visualizes sites of interest using the predicted extension. AVAILABILITY AND IMPLEMENTATION: ArchTEx requires Java 1.6 or newer. Source code and the compiled program are freely available at http://sourceforge.net/projects/archtex/ CONTACT: mjbuck@buffalo.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.