DNA unwinding in the CYC1 and DED1 yeast promoters

Deep learning of the regulatory grammar of yeast 5' untranslated regions from 500,000 random sequences

Our ability to predict protein expression from DNA sequence alone remains poor, reflecting our limited understanding of cis-regulatory grammar and hampering the design of engineered genes for synthetic biology applications. Here, we generate a model that predicts the protein expression of the 5′ untranslated region (UTR) of mRNAs in the yeast Saccharomyces cerevisiae. We constructed a library of half a million 50-nucleotide-long random 5′ UTRs and assayed their activity in a massively parallel growth selection experiment. The resulting data allow us to quantify the impact on protein expression of Kozak sequence composition, upstream open reading frames (uORFs), and secondary structure. We trained a convolutional neural network (CNN) on the random library and showed that it performs well at predicting the protein expression of both a held-out set of the random 5′ UTRs as well as native S. cerevisiae 5′ UTRs. The model additionally was used to computationally evolve highly active 5′ UTRs. We confirmed experimentally that the great majority of the evolved sequences led to higher protein expression rates than the starting sequences, demonstrating the predictive power of this model.

DNA tracts composed of only two bases concentrate in gene promoters

We have previously shown that long DNA tracts composed of only two of the bases ("binary DNA") are highly overrepresented in sequenced eukaryotic genomes. Here we examine gene promoter regions, by superposing all genes in a chromosome at their transcription start sites. We find that of the four motifs made of two bases, three are concentrated in gene promoters: Purine/pyrimidine tracts are highly overrepresented in the promoters of yeast chromosome IV, in Caenorhabditis elegans chromosome I, in Arabidopsis thaliana chromosome 2, and in human chromosomes 14, 21, and 22 (a subset). AT-rich tracts (W tracts) are enriched in the same chromosomes, as well as in Drosophila melanogaster chromosome 2R and in an archeon, M. jannaschii. A third motif, K.M tracts, shows some concentration in D. melanogaster promoters. A propensity of binary DNA to unwind is proposed to explain the high presence of the two-base motifs in gene promoters.

Coevolution of DNA-interacting proteins and genome "dialect"

Several species-specific characteristics of genome organization that are superimposed on its coding aspects were proposed earlier, including genome signature (GS), genome accent, and compositional spectrum (CS). These notions could be considered as representatives of genome dialect (GD). We measured within the Proteobacteria some GD representatives, the relative abundance of dinucleotides or GS, the profiles of occurrence of 10 nucleotide words (CS), and the profiles of occurrence of 20 nucleotide words, using a degenerate two-letter alphabet (purine-pyrimidine compositional spectra [PPCS]). Here, we show that the evolutionary distances between DNA repair and recombination orthologous enzymes (especially those of the nucleotide excision repair system) are highly correlated with PPCS and GS distances. Orthologous proteins involved in structural or metabolic processes (control group) have significantly lower correlations of their evolutionary distances with the PPCS and GS distances. We hypothesize that the high correlation of the evolutionary distances of the DNA repair orthologous enzymes with their GD is a result of the coevolution of the DNA repair enzymes' structures and GDs. Species GDs could be substantially influenced by the function of DNA polymerase I (the bacterial major DNA repair polymerase). This might cause the correlation of species GDs differentiation with evolutionary changes of species DNA polymerase I. Simultaneously, the structures of DNA repair-recombination enzymes might be evolutionarily sensitive and responsive to changes in the structure of their substrate-the DNA (including those that are represented by GD differentiation). We further discuss the rationale and mechanisms of the hypothesized coevolution. We suggest that stress might be an important cause of changes in the repair-recombination genes and the GD and the trigger of the aforementioned coevolution process. Other triggers might be massive horizontal gene transfer and ecological selection.

The over-representation of binary DNA tracts in seven sequenced chromosomes

BACKGROUND

DNA tracts composed of only two bases are possible in six combinations: A+G (purines, R), C+T (pyrimidines, Y), G+T (Keto, K), A+C (Imino, M), A+T (Weak, W) and G+C (Strong, S). It is long known that all-pyrimidine tracts, complemented by all-purines tracts ("R.Y tracts"), are excessively present in analyzed DNA. We have previously shown that R.Y tracts are in vast excess in yeast promoters, and brought evidence for their role in gene regulation. Here we report the systematic mapping of all six binary combinations on the level of complete sequenced chromosomes, as well as in their different subregions.

RESULTS

DNA tracts composed of the above binary base combinations have been mapped in seven sequenced chromosomes: Human chromosomes 21 and 22 (the major contigs); Drosophila melanogaster chr. 2R; Caenorhabditis elegans chr. I; Arabidopsis thaliana chr. II; Saccharomyces cerevisiae chr. IV and M. jannaschii. A huge over-representation, reaching million-folds, has been found for very long tracts of all binary motifs except S, in each of the seven organisms. Long R.Y tracts are the most excessive, except in D. melanogaster, where the K.M motif predominates. S (G, C rich) tracts are in excess mainly in CpG islands; the W motif predominates in bacteria. Many excessively long W tracts are nevertheless found also in the archeon and in the eukaryotes. The survey of complete chromosomes enables us, for the first time, to map systematically the intergenic regions. In human and other chromosomes we find the highest over-representation of the binary DNA tracts in the intergenic regions. These over-representations are only partly explainable by the presence of interspersed elements.

CONCLUSIONS

The over-representation of long DNA tracts composed of five of the above motifs is the largest deviation from randomness so far established for DNA, and this in a wide range of eukaryotic and archeal chromosomes. A propensity for ready DNA unwinding is proposed as the functional role, explaining the evolutionary conservation of the huge excesses observed.

TRACTS: A program to map oligopurine.oligopyrimidine and other binary DNA tracts

A program to map the locations and frequencies of DNA tracts composed of only two bases ('Binary DNA') is described. The program, TRACTS (URL http://bioportal.weizmann.ac.il/tracts/tracts.html and/or http://bip.weizmann.ac.il/miwbin/servers/tracts) is of interest because long tracts composed of only two bases are highly over-represented in most genomes. In eukaryotes, oligopurine.oligopyrimidine tracts ('R.Y tracts') are found in the highest excess. In prokaryotes, W tracts predominate (A,T 'rich'). A pre-program, ANEX, parses database annotation files of GenBank and EMBL, to produce a convenient one-line list of every gene (exon, intron) in a genome. The main unit lists and analyzes tracts of the three possible binary pairs (R.Y, K.M and S;W). As an example, the results of R.Y tract mapping of mammalian gene p53 is described.

Evidence implicating Ku antigen as a structural factor in RNA polymerase II-mediated transcription

Ku antigen is an abundant nuclear protein with multiple functions that depend mainly on Ku's prolific and highly verstatile interactions with DNA. We have shown previously that the direct binding of Ku in vitro to negative regulatory element 1 (NRE1), a transcriptional regulatory element in the long terminal repeat of mouse mammary tumour virus, correlates with the regulation of viral transcription by Ku. In this study, we have sought to explore the interaction of Ku with NRE1 in vivo in yeast one-hybrid experiments. Unexpectedly, we observed that human Ku70 carrying a transcriptional activation domain from the yeast Gal4 protein induced transcription of yeast reporter genes pleiotrophically, independent of NRE1, promoter, reporter gene and chromosomal location. Ku80 with the same activation domain had no effect on transcription when expressed alone, but reconstituted activation when co-expressed with native human Ku70. The requirements for transcriptional activation by Ku-Gal4 activation domain proteins correlated with previous descriptions of the requirements for DNA sequence-independent DNA binding by Ku, but were distinct from determinants for DNA-end binding by a truncated Ku heterodimer determined recently by crystallography. These results suggest a preferential targeting of Ku to transcriptionally active chromatin that indicate a possible function for Ku within the RNA polymerase II holoenzyme.