Intrinsic Sequence Preferences Research Articles

Computational prediction and characterization of cell-type-specific and shared binding sites.

Cell-type-specific gene expression is maintained in large part by transcription factors (TFs) selectively binding to distinct sets of sites in different cell types. Recent research works have provided evidence that such cell-type-specific binding is determined by TF's intrinsic sequence preferences, cooperative interactions with co-factors, cell-type-specific chromatin landscapes and 3D chromatin interactions. However, computational prediction and characterization of cell-type-specific and shared binding sites is rarely studied. In this article, we propose two computational approaches for predicting and characterizing cell-type-specific and shared binding sites by integrating multiple types of features, in which one is based on XGBoost and another is based on convolutional neural network (CNN). To validate the performance of our proposed approaches, ChIP-seq datasets of 10 binding factors were collected from the GM12878 (lymphoblastoid) and K562 (erythroleukemic) human hematopoietic cell lines, each of which was further categorized into cell-type-specific (GM12878- and K562-specific) and shared binding sites. Then, multiple types of features for these binding sites were integrated to train the XGBoost- and CNN-based models. Experimental results show that our proposed approaches significantly outperform other competing methods on three classification tasks. Moreover, we identified independent feature contributions for cell-type-specific and shared sites through SHAP values and explored the ability of the CNN-based model to predict cell-type-specific and shared binding sites by excluding or including DNase signals. Furthermore, we investigated the generalization ability of our proposed approaches to different binding factors in the same cellular environment. The source code is available at: https://github.com/turningpoint1988/CSSBS. Supplementary data are available at Bioinformatics online.

Translation initiation factor eIF4G1 preferentially binds yeast transcript leaders containing conserved oligo-uridine motifs.

Translational control of gene expression plays essential roles in cellular stress responses and organismal development by enabling rapid, selective, and localized control of protein production. Translational regulation depends on context-dependent differences in the protein output of mRNAs, but the key mRNA features that distinguish efficiently translated mRNAs are largely unknown. Here, we comprehensively determined the RNA-binding preferences of the eukaryotic initiation factor 4G (eIF4G) to assess whether this core translation initiation factor has intrinsic sequence preferences that may contribute to preferential translation of specific mRNAs. We identified a simple RNA sequence motif—oligo-uridine—that mediates high-affinity binding to eIF4G in vitro. Oligo(U) motifs occur naturally in the transcript leader (TL) of hundreds of yeast genes, and mRNAs with unstructured oligo(U) motifs were enriched in immunoprecipitations against eIF4G. Ribosome profiling following depletion of eIF4G in vivo showed preferentially reduced translation of mRNAs with long TLs, including those that contain oligo(U). Finally, TL oligo(U) elements are enriched in genes with regulatory roles and are conserved between yeast species, consistent with an important cellular function. Taken together, our results demonstrate RNA sequence preferences for a general initiation factor, which cells potentially exploit for translational control of specific mRNAs.

Nucleosome Free Regions in Yeast Promoters Result from Competitive Binding of Transcription Factors That Interact with Chromatin Modifiers

Because DNA packaging in nucleosomes modulates its accessibility to transcription factors (TFs), unraveling the causal determinants of nucleosome positioning is of great importance to understanding gene regulation. Although there is evidence that intrinsic sequence specificity contributes to nucleosome positioning, the extent to which other factors contribute to nucleosome positioning is currently highly debated. Here we obtained both in vivo and in vitro reference maps of positions that are either consistently covered or free of nucleosomes across multiple experimental data-sets in Saccharomyces cerevisiae. We then systematically quantified the contribution of TF binding to nucleosome positiong using a rigorous statistical mechanics model in which TFs compete with nucleosomes for binding DNA. Our results reconcile previous seemingly conflicting results on the determinants of nucleosome positioning and provide a quantitative explanation for the difference between in vivo and in vitro positioning. On a genome-wide scale, nucleosome positioning is dominated by the phasing of nucleosome arrays over gene bodies, and their positioning is mainly determined by the intrinsic sequence preferences of nucleosomes. In contrast, larger nucleosome free regions in promoters, which likely have a much more significant impact on gene expression, are determined mainly by TF binding. Interestingly, of the 158 yeast TFs included in our modeling, we find that only 10–20 significantly contribute to inducing nucleosome-free regions, and these TFs are highly enriched for having direct interations with chromatin remodelers. Together our results imply that nucleosome free regions in yeast promoters results from the binding of a specific class of TFs that recruit chromatin remodelers.

A Compendium of Nucleosome and Transcript Profiles Reveals Determinants of Chromatin Architecture and Transcription

Nucleosomes in all eukaryotes examined to date adopt a characteristic architecture within genes and play fundamental roles in regulating transcription, yet the identity and precise roles of many of the trans-acting factors responsible for the establishment and maintenance of this organization remain to be identified. We profiled a compendium of 50 yeast strains carrying conditional alleles or complete deletions of genes involved in transcriptional regulation, histone biology, and chromatin remodeling, as well as compounds that target transcription and histone deacetylases, to assess their respective roles in nucleosome positioning and transcription. We find that nucleosome patterning in genes is affected by many factors, including the CAF-1 complex, Spt10, and Spt21, in addition to previously reported remodeler ATPases and histone chaperones. Disruption of these factors or reductions in histone levels led genic nucleosomes to assume positions more consistent with their intrinsic sequence preferences, with pronounced and specific shifts of the +1 nucleosome relative to the transcription start site. These shifts of +1 nucleosomes appear to have functional consequences, as several affected genes in Ino80 mutants exhibited altered expression responses. Our parallel expression profiling compendium revealed extensive transcription changes in intergenic and antisense regions, most of which occur in regions with altered nucleosome occupancy and positioning. We show that the nucleosome-excluding transcription factors Reb1, Abf1, Tbf1, and Rsc3 suppress cryptic transcripts at their target promoters, while a combined analysis of nucleosome and expression profiles identified 36 novel transcripts that are normally repressed by Tup1/Cyc8. Our data confirm and extend the roles of chromatin remodelers and chaperones as major determinants of genic nucleosome positioning, and these data provide a valuable resource for future studies.

Global remodeling of nucleosome positions in C. elegans

BackgroundEukaryotic chromatin architecture is affected by intrinsic histone-DNA sequence preferences, steric exclusion between nucleosome particles, formation of higher-order structures, and in vivo activity of chromatin remodeling enzymes.ResultsTo disentangle sequence-dependent nucleosome positioning from the other factors, we have created two high-throughput maps of nucleosomes assembled in vitro on genomic DNA from the nematode worm Caenorhabditis elegans. A comparison of in vitro nucleosome positions with those observed in a mixed-stage, mixed-tissue population of C. elegans cells reveals that in vivo sequence preferences are modified on the genomic scale. Indeed, G/C dinucleotides are predicted to be most favorable for nucleosome formation in vitro but not in vivo. Nucleosome sequence read coverage in vivo is distinctly lower in chromosome arms than in central regions; the observed changes in apparent nucleosome sequence specificity, likely due to genome-wide chromatin remodeler activity, contribute to the formation of these megabase-scale chromatin domains. We also observe that the majority of well-positioned in vivo nucleosomes do not occupy thermodynamically favorable sequences observed in vitro. Finally, we find that exons are intrinsically more amenable to nucleosome formation compared to introns. Nucleosome occupancy of introns and exons consistently increases with G/C content in vitro but not in vivo, in agreement with our observation that G/C dinucleotide enrichment does not strongly promote in vivo nucleosome formation.ConclusionsOur findings highlight the importance of both sequence specificity and active nucleosome repositioning in creating large-scale chromatin domains, and the antagonistic roles of intrinsic sequence preferences and chromatin remodelers in C. elegans.Sequence read data has been deposited into Sequence Read Archive (http://www.ncbi.nlm.nih.gov/sra; accession number SRA050182). Additional data, software and computational predictions are available on the Nucleosome Explorer website (http://nucleosome.rutgers.edu).

71 Genomic regions flanking E-box binding sites influence DNA binding specificity of bHLH transcription factors through DNA shape

Transcriptional regulation of gene expression is enacted mainly through binding of transcription factors (TFs) to specific, short DNA sites in cis-regulatory regions of genes. Most TFs are members of protein families that share a common DNA-binding domain and thus recognize similar DNA-binding sequences. It is not well understood why paralogous TFs often bind different genomic target sites in vivo to effect different regulatory programs, despite apparently recognizing the same sequence motifs. Here, we designed custom protein-binding microarrays (PBMs) to analyze the DNA-binding specificities of two Saccharomyces cerevisiae basic helix-loop-helix (bHLH) proteins, Tye7 and Cbf1, as a model system. Our data reveal that E-box DNA-binding sequences (CAnnTG), when tested in the context of their native genomic flanking sequences, are bound differently by Cbf1 and Tye7. Computational models of the PBM data indicate that DNA sequence features located in the genomic sequences outside the E-box contribute to DNA-binding specificity in vitro. Our analyses suggest that these flanking regions affect DNA-binding specificity indirectly by influencing the three-dimensional structure of the E-box binding sites. Finally, we show that these subtle differences in intrinsic sequence preferences of Cbf1 and Tye7 in vitro help to explain their differential DNA-binding preferences in vivo. Our results provide further evidence that the local shape of DNA-binding sites may be an important feature in distinguishing the DNA-binding preferences among paralogous TFs and thus may play a widespread role in determining how transcriptional regulatory specificity within TF families is achieved.

Quantitative Test of the Barrier Nucleosome Model for Statistical Positioning of Nucleosomes Up- and Downstream of Transcription Start Sites

The positions of nucleosomes in eukaryotic genomes determine which parts of the DNA sequence are readily accessible for regulatory proteins and which are not. Genome-wide maps of nucleosome positions have revealed a salient pattern around transcription start sites, involving a nucleosome-free region (NFR) flanked by a pronounced periodic pattern in the average nucleosome density. While the periodic pattern clearly reflects well-positioned nucleosomes, the positioning mechanism is less clear. A recent experimental study by Mavrich et al. argued that the pattern observed in Saccharomyces cerevisiae is qualitatively consistent with a “barrier nucleosome model,” in which the oscillatory pattern is created by the statistical positioning mechanism of Kornberg and Stryer. On the other hand, there is clear evidence for intrinsic sequence preferences of nucleosomes, and it is unclear to what extent these sequence preferences affect the observed pattern. To test the barrier nucleosome model, we quantitatively analyze yeast nucleosome positioning data both up- and downstream from NFRs. Our analysis is based on the Tonks model of statistical physics which quantifies the interplay between the excluded-volume interaction of nucleosomes and their positional entropy. We find that although the typical patterns on the two sides of the NFR are different, they are both quantitatively described by the same physical model with the same parameters, but different boundary conditions. The inferred boundary conditions suggest that the first nucleosome downstream from the NFR (the +1 nucleosome) is typically directly positioned while the first nucleosome upstream is statistically positioned via a nucleosome-repelling DNA region. These boundary conditions, which can be locally encoded into the genome sequence, significantly shape the statistical distribution of nucleosomes over a range of up to ∼1,000 bp to each side.

Predicting nucleosome positions on the DNA: combining intrinsic sequence preferences and remodeler activities.

Nucleosome positions on the DNA are determined by the intrinsic affinities of histone proteins to a given DNA sequence and by the ATP-dependent activities of chromatin remodeling complexes that can translocate nucleosomes with respect to the DNA. Here, we report a theoretical approach that takes into account both contributions. In the theoretical analysis two types of experiments have been considered: in vitro experiments with a single reconstituted nucleosome and in vivo genome-scale mapping of nucleosome positions. The effect of chromatin remodelers was described by iteratively redistributing the nucleosomes according to certain rules until a new steady state was reached. Three major classes of remodeler activities were identified: (i) the establishment of a regular nucleosome spacing in the vicinity of a strong positioning signal acting as a boundary, (ii) the enrichment/depletion of nucleosomes through amplification of intrinsic DNA-sequence-encoded signals and (iii) the removal of nucleosomes from high-affinity binding sites. From an analysis of data for nucleosome positions in resting and activated human CD4+ T cells [Schones et al., Cell 132, p. 887] it was concluded that the redistribution of a nucleosome map to a new state is greatly facilitated if the remodeler complex translocates the nucleosome with a preferred directionality.

Two distinct modes of nucleosome modulation associated with different degrees of dependence of nucleosome positioning on the underlying DNA sequence

BackgroundThe nucleosome is the fundamental unit of eukaryotic genomes. Its positioning plays a central role in diverse cellular processes that rely on access to genomic DNA. Experimental evidence suggests that the genomic DNA sequence is one important determinant of nucleosome positioning. Yet it is less clear whether the role of the underlying DNA sequence in nucleosome positioning varies across different promoters. Whether different determinants of nucleosome positioning have characteristic influences on nucleosome modulation also remains to be elucidated.ResultsWe identified two typical promoter classes in yeast associated with high or low dependence of nucleosome positioning on the underlying DNA sequence, respectively. Importantly, the two classes have low or high intrinsic sequence preferences for nucleosomes, respectively. The two classes are further distinguished by multiple promoter features, including nucleosome occupancy, nucleosome fuzziness, H2A.Z occupancy, changes in nucleosome positions before and after transcriptional perturbation, and gene activity. Both classes have significantly high turnover rates of histone H3, but employ distinct modes of nucleosome modulation: The first class is characterized by hyperacetylation, whereas the second class is highly regulated by ATP-dependent chromatin remodelling.ConclusionOur results, coupled with the known features of nucleosome modulation, suggest that the two distinct modes of nucleosome modulation selectively employed by different genes are linked with the intrinsic sequence preferences for nucleosomes. The difference in modes of nucleosome modulation can account for the difference in the contribution of DNA sequence to nucleosome positioning between both promoter classes.

The DNA-encoded nucleosome organization of a eukaryotic genome

Nucleosome organization is critical for gene regulation. In living cells this organization is determined by multiple factors, including the action of chromatin remodellers, competition with site-specific DNA-binding proteins, and the DNA sequence preferences of the nucleosomes themselves. However, it has been difficult to estimate the relative importance of each of these mechanisms in vivo, because in vivo nucleosome maps reflect the combined action of all influencing factors. Here we determine the importance of nucleosome DNA sequence preferences experimentally by measuring the genome-wide occupancy of nucleosomes assembled on purified yeast genomic DNA. The resulting map, in which nucleosome occupancy is governed only by the intrinsic sequence preferences of nucleosomes, is similar to in vivo nucleosome maps generated in three different growth conditions. In vitro, nucleosome depletion is evident at many transcription factor binding sites and around gene start and end sites, indicating that nucleosome depletion at these sites in vivo is partly encoded in the genome. We confirm these results with a micrococcal nuclease-independent experiment that measures the relative affinity of nucleosomes for approximately 40,000 double-stranded 150-base-pair oligonucleotides. Using our in vitro data, we devise a computational model of nucleosome sequence preferences that is significantly correlated with in vivo nucleosome occupancy in Caenorhabditis elegans. Our results indicate that the intrinsic DNA sequence preferences of nucleosomes have a central role in determining the organization of nucleosomes in vivo.

Comparison of the Differential Context-dependence of DNA Deamination by APOBEC Enzymes: Correlation with Mutation Spectra in Vivo

To investigate the extent to which in vivo mutation spectra might reflect the intrinsic specificities of active mutators, genetic and biochemical assays were used to analyse the DNA target specificities of cytidine deaminases of the APOBEC family. The results reveal the critical importance of nucleotides immediately 5′ of the targeted C for the specificity of all three enzymes studied (AID, APOBEC1 and APOBEC3G). At position −1, APOBEC1 showed a marked preference for dT, AID for dA/dG and APOBEC3G a strong preference for dC. Furthermore, AID and APOBEC3G showed distinct dependence on the nucleotide at position −2 with dA/dT being favoured by AID and dC by APOBEC3G. Most if not all activity of the recombinant deaminases on free dC could be attributed to low-level contamination by host enzymes. The target preference of APOBEC3G supports it being a major but possibly not sole contributor to HIV hypermutation without making it a dominant contribution to general HIV sequence variation. The specificity of AID as deduced from the genetic assay (which relies on inactivation of sacB of Bacillus subtilis) agrees well with that deduced by Pham et al. using an in vitro assay although we postulate that major intrinsic mutational hotspots in immunoglobulin V genes in vivo might reflect favoured sites of AID action being generated by proximal DNA targets located on opposite DNA strands. The target specificity of AID also accords with the spectrum of mutations observed in B lymphoma-associated oncogenes. The possibility of deaminase involvement in non-lymphoid human tumours is hinted at by tissue-specific differences in the spectra of dC transitions in tumour-suppressor genes. Thus, the patterns of hypermutation in antibodies and retroviruses owe much to the intrinsic sequence preferences of the AID/APOBEC family of DNA deaminases: analogous biases might also contribute to the spectra of cancer-associated mutation.

A theoretical investigation on the sequence selective binding of daunomycin to double-stranded polynucleotides.

Theoretical computations are performed on the structural and energetical factors involved in the sequence selective binding of daunomycin (DNM) to six representative self-complementary double-stranded hexanucleotides: d(CGTACG)2,d(CGATCG)2,d(CITACI)2, d(TATATA)2, d(CGCGCG)2 and d(TACGTA)2. The conformational angles of the hexanucleotides are fixed in values found in the representative crystal structure of the d(CGTACG)2-DNM complex. The intermolecular DNM-hexanucleotide interaction energies and the conformational energy changes of DNM upon binding are computed and optimized in the framework of the SIBFA procedure, which uses empirical formulas based on ab initio SCF computations. Among the two regularly alternating hexanucleotides, d(TATATA)2 and d(CGCGCG)2, a stronger binding is predicted for the former, in agreement with experimental results obtained with poly(dA-dT).poly(dA-dT) and poly(dG-dC).poly(dG-dC). Altogether, however, among the six investigated sequences, the strongest complexes are computed for the mixed hexanucleotides d(CGATCG)2 and d(CGTACG)2, containing the intercalation site between two CG base pairs and an adjacent TA base pair. This situation may be related to the increased affinity of DNM for GC rich DNA's and to the situation in the crystal structure of the DNM-d(CGTACG)2 complex. Analysis of the intrinsic base sequence preferences expressed by the individual constituents of DNM, namely the daunosamine side chain, the chromophore ring and its two 9-hydroxyl and 9-acetoxy substituents, reveals that the overall sequence preference found is the result of a rather intricate interplay of intrinsic sequence preferences, in particular at the level of daunosamine and the 9-hydroxyl substituent.(ABSTRACT TRUNCATED AT 250 WORDS)