A high-resolution protein architecture of the budding yeast genome

Revisiting the model for coactivator recruitment: Med15 can select its target sites independent of promoter-bound transcription factors

Activation domains (ADs) within transcription factors (TFs) induce gene expression by recruiting coactivators such as the Mediator complex. Coactivators lack DNA binding domains (DBDs) and are assumed to passively follow their recruiting TFs. This is supported by direct AD-coactivator interactions seen in vitro but has not yet been tested in living cells. To examine that, we targeted two Med15-recruiting ADs to a range of budding yeast promoters through fusion with different DBDs. The DBD-AD fusions localized to hundreds of genomic sites but recruited Med15 and induced transcription in only a subset of bound promoters, characterized by a fuzzy-nucleosome architecture. Direct DBD-Med15 fusions shifted DBD localization towards fuzzy-nucleosome promoters, including promoters devoid of the endogenous Mediator. We propose that Med15, and perhaps other coactivators, possess inherent promoter preference and thus actively contribute to the selection of TF-induced genes.

TFIIB-Termination Factor Interaction Affects Termination of Transcription on Genome-Wide Scale

Apart from its well-established role in the initiation of transcription, the general transcription factor TFIIB has been implicated in the termination step as well. The ubiquity of TFIIB involvement in termination as well as mechanistic details of its termination function, however, remain largely unexplored. Using GRO-seq analyses, we compared the terminator readthrough phenotype in the sua7-1 mutant (TFIIBsua7-1) and the isogenic wild type (TFIIBWT) strains. Approximately 74% of genes analyzed exhibited a 2-3-fold increase in readthrough of the poly(A)-termination signal in the TFIIBsua7-1 mutant compared to TFIIBWT cells. To understand the mechanistic basis of TFIIB's role in termination, we performed the mass spectrometry of TFIIB-affinity purified from chromatin and soluble cellular fractions-from TFIIBsua7-1 and TFIIBWT cells. TFIIB purified from the chromatin fraction of TFIIBWT cells exhibited significant enrichment of CF1A and Rat1 termination complexes. There was, however, a drastic decrease in TFIIB interaction with CF1A and Rat1 complexes in the TFIIBsua7-1 mutant. ChIP assays revealed about a 90% decline in the recruitment of termination factors in the TFIIBsua7-1 mutant compared to wild type cells. The overall conclusion of these results is that TFIIB affects the termination of transcription on a genome-wide scale, and the TFIIB-termination factor interaction plays a crucial role in the process.

Systematic dissection of sequence features affecting binding specificity of a pioneer factor reveals binding synergy between FOXA1 and AP-1

Despite the unique ability of pioneer factors (PFs) to target nucleosomal sites in closed chromatin, they only bind a small fraction of their genomic motifs. The underlying mechanism of this selectivity is not well understood. Here, we design a high-throughput assay called chromatin immunoprecipitation with integrated synthetic oligonucleotides (ChIP-ISO) to systematically dissect sequence features affecting the binding specificity of a classic PF, FOXA1, in human A549 cells. Combining ChIP-ISO with in vitro and neural network analyses, we find that (1) FOXA1 binding is strongly affected by co-binding transcription factors (TFs) AP-1 and CEBPB; (2) FOXA1 and AP-1 show binding cooperativity in vitro; (3) FOXA1's binding is determined more by local sequences than chromatin context, including eu-/heterochromatin; and (4) AP-1 is partially responsible for differential binding of FOXA1 in different cell types. Our study presents a framework for elucidating genetic rules underlying PF binding specificity and reveals a mechanism for context-specific regulation of its binding.

Transcriptomic balance and optimal growth are determined by cell size

Cell size and growth are intimately related across the evolutionary scale, but whether cell size is important to attain maximal growth or fitness is still an open question. We show that growth rate is a non-monotonic function of cell volume, with maximal values around the critical size of wild-type yeast cells. The transcriptome of yeast and mouse cells undergoes a relative inversion in response to cell size, which we associate theoretically and experimentally with the necessary genome-wide diversity in RNA polymerase II affinity for promoters. Although highly expressed genes impose strong negative effects on fitness when the DNA/mass ratio is reduced, transcriptomic alterations mimicking the relative inversion by cell size strongly restrain cell growth. In all, our data indicate that cells set the critical size to obtain a properly balanced transcriptome and, as a result, maximize growth and fitness during proliferation.

Genome-wide analysis of TFIIB's role in termination of transcription

This study provides evidence that the role of TFIIB extends beyond initiation to include the termination step of transcription. Using GRO-seq analyses, we compared terminator readthrough phenotype in sua7-1 mutant (TFIIB sua7-1 ) and the isogenic wild type (TFIIB WT ) strains. Approximately 74% of genes analyzed exhibited a 2-3-fold increase in readthrough of the poly(A)-termination signal in the TFIIB sua7-1 mutant compared to TFIIB WT cells. Mass spectrometry of affinity purified TFIIB from chromatin fraction found TFIIB exhibiting interaction with CF1A and Rat1 termination complexes in TFIIB WT cells. There was, however, a drastic decrease in TFIIB interaction with CF1A and Rat1 termination complexes in the TFIIB sua7-1 mutant. ChIP assays revealed about 90% decline in recruitment of termination factors in TFIIB sua7-1 mutant compared to wild type cells. The overall conclusion of these results is that TFIIB affects termination of transcription on a genome-wide scale, and TFIIB-termination factor interaction may play a crucial role in the process.

Chromatin endogenous cleavage provides a global view of RNA polymerase II transcription kinetics

Chromatin immunoprecipitation (ChIP-seq) is the most common approach to observe global binding of proteins to DNA in vivo. The occupancy of transcription factors (TFs) from ChIP-seq agrees well with an alternative method, chromatin endogenous cleavage (ChEC-seq2). However, ChIP-seq and ChEC-seq2 reveal strikingly different patterns of enrichment of yeast RNA polymerase II. We hypothesized that this reflects distinct populations of RNAPII, some of which are captured by ChIP-seq and some of which are captured by ChEC-seq2. RNAPII association with enhancers and promoters - predicted from biochemical studies - is detected well by ChEC-seq2 but not by ChIP-seq. Enhancer/promoter bound RNAPII correlates with transcription levels and matches predicted occupancy based on published rates of enhancer recruitment, preinitiation assembly, initiation, elongation and termination. The occupancy from ChEC-seq2 allowed us to develop a stochastic model for global kinetics of RNAPII transcription which captured both the ChEC-seq2 data and changes upon chemical-genetic perturbations to transcription. Finally, RNAPII ChEC-seq2 and kinetic modeling suggests that a mutation in the Gcn4 transcription factor that blocks interaction with the NPC destabilizes promoter-associated RNAPII without altering its recruitment to the enhancer.

Massively integrated coexpression analysis reveals transcriptional regulation, evolution and cellular implications of the yeast noncanonical translatome

BACKGROUND

Recent studies uncovered pervasive transcription and translation of thousands of noncanonical open reading frames (nORFs) outside of annotated genes. The contribution of nORFs to cellular phenotypes is difficult to infer using conventional approaches because nORFs tend to be short, of recent de novo origins, and lowly expressed. Here we develop a dedicated coexpression analysis framework that accounts for low expression to investigate the transcriptional regulation, evolution, and potential cellular roles of nORFs in Saccharomyces cerevisiae.

RESULTS

Our results reveal that nORFs tend to be preferentially coexpressed with genes involved in cellular transport or homeostasis but rarely with genes involved in RNA processing. Mechanistically, we discover that young de novo nORFs located downstream of conserved genes tend to leverage their neighbors' promoters through transcription readthrough, resulting in high coexpression and high expression levels. Transcriptional piggybacking also influences the coexpression profiles of young de novo nORFs located upstream of genes, but to a lesser extent and without detectable impact on expression levels. Transcriptional piggybacking influences, but does not determine, the transcription profiles of de novo nORFs emerging nearby genes. About 40% of nORFs are not strongly coexpressed with any gene but are transcriptionally regulated nonetheless and tend to form entirely new transcription modules. We offer a web browser interface ( https://carvunislab.csb.pitt.edu/shiny/coexpression/ ) to efficiently query, visualize, and download our coexpression inferences.

CONCLUSIONS

Our results suggest that nORF transcription is highly regulated. Our coexpression dataset serves as an unprecedented resource for unraveling how nORFs integrate into cellular networks, contribute to cellular phenotypes, and evolve.

nucMACC: An MNase-seq pipeline to identify structurally altered nucleosomes in the genome

Micrococcal nuclease sequencing is the state-of-the-art method for determining chromatin structure and nucleosome positioning. Data analysis is complex due to the AT-dependent sequence bias of the endonuclease and the requirement for high sequencing depth. Here, we present the nucleosome-based MNase accessibility (nucMACC) pipeline unveiling the regulatory chromatin landscape by measuring nucleosome accessibility and stability. The nucMACC pipeline represents a systematic and genome-wide approach for detecting unstable ("fragile") nucleosomes. We have characterized the regulatory nucleosome landscape in Drosophila melanogaster, Saccharomyces cerevisiae, and mammals. Two functionally distinct sets of promoters were characterized, one associated with an unstable nucleosome and the other being nucleosome depleted. We show that unstable nucleosomes present intermediate states of nucleosome remodeling, preparing inducible genes for transcriptional activation in response to stimuli or stress. The presence of unstable nucleosomes correlates with RNA polymerase II proximal pausing. The nucMACC pipeline offers unparalleled precision and depth in nucleosome research and is a valuable tool for future nucleosome studies.

Genome-wide analysis of TFIIB's role in termination of transcription

Apart from its well-established role in initiation of transcription, the general transcription factor TFIIB has been implicated in the termination step as well. The ubiquity of TFIIB involvement in termination as well as mechanistic details of its termination function, however, remains largely unexplored. To determine the prevalence of TFIIB's role in termination, we performed GRO-seq analyses in sua7-1 mutant (TFIIB sua7-1 ) and the isogenic wild type (TFIIB WT ) strains of yeast. Almost a three-fold increase in readthrough of the poly(A)-termination signal was observed in TFIIB sua7-1 mutant compared to the TFIIB WT cells. Of all genes analyzed in this study, nearly 74% genes exhibited a statistically significant increase in terminator readthrough in the mutant. To gain an understanding of the mechanistic basis of TFIIB involvement in termination, we performed mass spectrometry of TFIIB, affinity purified from chromatin and soluble cellular fractions, from TFIIB sua7-1 and TFIIB WT cells. TFIIB purified from the chromatin fraction of TFIIB WT cells exhibited significant enrichment of CF1A and Rat1 termination complexes. There was, however, a drastic decrease in TFIIB interaction with both CF1A and Rat1 termination complexes in TFIIB sua7-1 mutant. ChIP assay revealed that the recruitment of Pta1 subunit of CPF complex, Rna15 subunit of CF1 complex and Rat1 subunit of Rat1 complex registered nearly 90% decline in the mutant over wild type cells. The overall conclusion of these results is that TFIIB affects termination of transcription on a genome-wide scale, and TFIIB-termination factor interaction may play a crucial role in the process.

Nucleosomal DNA has topological memory

One elusive aspect of the chromosome architecture is how it constrains the DNA topology. Nucleosomes stabilise negative DNA supercoils by restraining a DNA linking number difference (∆Lk) of about -1.26. However, whether this capacity is uniform across the genome is unknown. Here, we calculate the ∆Lk restrained by over 4000 nucleosomes in yeast cells. To achieve this, we insert each nucleosome in a circular minichromosome and perform Topo-seq, a high-throughput procedure to inspect the topology of circular DNA libraries in one gel electrophoresis. We show that nucleosomes inherently restrain distinct ∆Lk values depending on their genomic origin. Nucleosome DNA topologies differ at gene bodies (∆Lk = -1.29), intergenic regions (∆Lk = -1.23), rDNA genes (∆Lk = -1.24) and telomeric regions (∆Lk = -1.07). Nucleosomes near the transcription start and termination sites also exhibit singular DNA topologies. Our findings demonstrate that nucleosome DNA topology is imprinted by its native chromatin context and persists when the nucleosome is relocated.

Genome-wide nucleosome and transcription factor responses to genetic perturbations reveal chromatin-mediated mechanisms of transcriptional regulation

Epigenetic mechanisms contribute to gene regulation by altering chromatin accessibility through changes in transcription factor (TF) and nucleosome occupancy throughout the genome. Despite numerous studies focusing on changes in gene expression, the intricate chromatin-mediated regulatory code remains largely unexplored on a comprehensive scale. We address this by employing a factor-agnostic, reverse-genetics approach that uses MNase-seq to capture genome-wide TF and nucleosome occupancies in response to the individual deletion of 201 transcriptional regulators in Saccharomyces cerevisiae, thereby assaying nearly one million mutant-gene interactions. We develop a principled approach to identify and quantify chromatin changes genome-wide, observing differences in TF and nucleosome occupancy that recapitulate well-established pathways identified by gene expression data. We also discover distinct chromatin signatures associated with the up- and downregulation of genes, and use these signatures to reveal regulatory mechanisms previously unexplored in expression-based studies. Finally, we demonstrate that chromatin features are predictive of transcriptional activity and leverage these features to reconstruct chromatin-based transcriptional regulatory networks. Overall, these results illustrate the power of an approach combining genetic perturbation with high-resolution epigenomic profiling; the latter enables a close examination of the interplay between TFs and nucleosomes genome-wide, providing a deeper, more mechanistic understanding of the complex relationship between chromatin organization and transcription.

ifDEEPre: large protein language-based deep learning enables interpretable and fast predictions of enzyme commission numbers

Accurate understanding of the biological functions of enzymes is vital for various tasks in both pathologies and industrial biotechnology. However, the existing methods are usually not fast enough and lack explanations on the prediction results, which severely limits their real-world applications. Following our previous work, DEEPre, we propose a new interpretable and fast version (ifDEEPre) by designing novel self-guided attention and incorporating biological knowledge learned via large protein language models to accurately predict the commission numbers of enzymes and confirm their functions. Novel self-guided attention is designed to optimize the unique contributions of representations, automatically detecting key protein motifs to provide meaningful interpretations. Representations learned from raw protein sequences are strictly screened to improve the running speed of the framework, 50 times faster than DEEPre while requiring 12.89 times smaller storage space. Large language modules are incorporated to learn physical properties from hundreds of millions of proteins, extending biological knowledge of the whole network. Extensive experiments indicate that ifDEEPre outperforms all the current methods, achieving more than 14.22% larger F1-score on the NEW dataset. Furthermore, the trained ifDEEPre models accurately capture multi-level protein biological patterns and infer evolutionary trends of enzymes by taking only raw sequences without label information. Meanwhile, ifDEEPre predicts the evolutionary relationships between different yeast sub-species, which are highly consistent with the ground truth. Case studies indicate that ifDEEPre can detect key amino acid motifs, which have important implications for designing novel enzymes. A web server running ifDEEPre is available at https://proj.cse.cuhk.edu.hk/aihlab/ifdeepre/ to provide convenient services to the public. Meanwhile, ifDEEPre is freely available on GitHub at https://github.com/ml4bio/ifDEEPre/.

DNA methylation-based high-resolution mapping of long-distance chromosomal interactions in nucleosome-depleted regions

3C-based methods have significantly advanced our understanding of 3D genome organization. However, it remains a formidable task to precisely capture long-range chromosomal interactions between individual loci, such as those between promoters and distal enhancers. Here, we present Methyltransferase Targeting-based chromosome Architecture Capture (MTAC), a method that maps the contacts between a target site (viewpoint) and the rest of the genome in budding yeast with high resolution and sensitivity. MTAC detects hundreds of intra- and inter-chromosomal interactions within nucleosome-depleted regions (NDRs) that cannot be captured by 4C, Hi-C, or Micro-C. By applying MTAC to various viewpoints, we find that (1) most long-distance chromosomal interactions detected by MTAC reflect tethering by the nuclear pore complexes (NPCs), (2) genes co-regulated by methionine assemble into inter-chromosomal clusters near NPCs upon activation, (3) mediated by condensin, the mating locus forms a highly specific interaction with the recombination enhancer (RE) in a mating-type specific manner, and (4) correlation of MTAC signals among NDRs reveal spatial mixing and segregation of the genome. Overall, these results demonstrate MTAC as a powerful tool to resolve fine-scale long-distance chromosomal interactions and provide insights into the 3D genome organization.

Exportin-1 functions as an adaptor for transcription factor-mediated docking of chromatin at the nuclear pore complex

Nuclear pore proteins (Nups) in yeast, flies and mammals physically interact with hundreds or thousands of chromosomal sites, which impacts transcriptional regulation. In budding yeast, transcription factors mediate interaction of Nups with enhancers of highly active genes. To define the molecular basis of this mechanism, we exploited a separation-of-function mutation in the Gcn4 transcription factor that blocks its interaction with the nuclear pore complex (NPC) without altering its DNA binding or activation domains. SILAC mass spectrometry revealed that this mutation reduces the interaction of Gcn4 with the highly conserved nuclear export factor Crm1/Xpo1. Crm1 both interacts with the same sites as Nups genome-wide and is required for Nup2 to interact with the yeast genome. In vivo, Crm1 undergoes extensive and stable interactions with the NPC. In vitro, Crm1 binds to Gcn4 and these proteins form a complex with the nuclear pore protein Nup2. Importantly, the interaction between Crm1 and Gcn4 does not require Ran-GTP, suggesting that it is not through the nuclear export sequence binding site. Finally, Crm1 stimulates DNA binding by Gcn4, supporting a model in which allosteric coupling between Crm1 binding and DNA binding permits docking of transcription factor-bound enhancers at the NPC.

General regulatory factors exert differential effects on nucleosome sliding activity of the ISW1a complex

BACKGROUND

Chromatin dynamics is deeply involved in processes that require access to DNA, such as transcriptional regulation. Among the factors involved in chromatin dynamics at gene regulatory regions are general regulatory factors (GRFs). These factors contribute to establishment and maintenance of nucleosome-depleted regions (NDRs). These regions are populated by nucleosomes through histone deposition and nucleosome sliding, the latter catalyzed by a number of ATP-dependent chromatin remodeling complexes, including ISW1a. It has been observed that GRFs can act as barriers against nucleosome sliding towards NDRs. However, the relative ability of the different GRFs to hinder sliding activity is currently unknown.

RESULTS

Considering this, we performed a comparative analysis for the main GRFs, with focus in their ability to modulate nucleosome sliding mediated by ISW1a. Among the GRFs tested in nucleosome remodeling assays, Rap1 was the only factor displaying the ability to hinder the activity of ISW1a. This effect requires location of the Rap1 cognate sequence on linker that becomes entry DNA in the nucleosome remodeling process. In addition, Rap1 was able to hinder nucleosome assembly in octamer transfer assays. Concurrently, Rap1 displayed the highest affinity for and longest dwell time from its target sequence, compared to the other GRFs tested. Consistently, through bioinformatics analyses of publicly available genome-wide data, we found that nucleosome occupancy and histone deposition in vivo are inversely correlated with the affinity of Rap1 for its target sequences in the genome.

CONCLUSIONS

Our findings point to DNA binding affinity, residence time and location at particular translational positions relative to the nucleosome core as the key features of GRFs underlying their roles played in nucleosome sliding and assembly.

Swd2/Cps35 determines H3K4 tri-methylation via interactions with Set1 and Rad6

BACKGROUND

Histone H3K4 tri-methylation (H3K4me3) catalyzed by Set1/COMPASS, is a prominent epigenetic mark found in promoter-proximal regions of actively transcribed genes. H3K4me3 relies on prior monoubiquitination at the histone H2B (H2Bub) by Rad6 and Bre1. Swd2/Cps35, a Set1/COMPASS component, has been proposed as a key player in facilitating H2Bub-dependent H3K4me3. However, a more comprehensive investigation regarding the relationship among Rad6, Swd2, and Set1 is required to further understand the mechanisms and functions of the H3K4 methylation.

RESULTS

We investigated the genome-wide occupancy patterns of Rad6, Swd2, and Set1 under various genetic conditions, aiming to clarify the roles of Set1 and Rad6 for occupancy of Swd2. Swd2 peaks appear on both the 5' region and 3' region of genes, which are overlapped with its tightly bound two complexes, Set1 and cleavage and polyadenylation factor (CPF), respectively. In the absence of Rad6/H2Bub, Set1 predominantly localized to the 5' region of genes, while Swd2 lost all the chromatin binding. However, in the absence of Set1, Swd2 occupancy near the 5' region was impaired and rather increased in the 3' region.

CONCLUSIONS

This study highlights that the catalytic activity of Rad6 is essential for all the ways of Swd2's binding to the transcribed genes and Set1 redistributes the Swd2 to the 5' region for accomplishments of H3K4me3 in the genome-wide level.

Transcript elongation by RNA polymerase II in plants: factors, regulation and impact on gene expression

Transcriptional elongation by RNA polymerase II (RNAPII) through chromatin is a dynamic and highly regulated step of eukaryotic gene expression. A combination of transcript elongation factors (TEFs) including modulators of RNAPII activity and histone chaperones facilitate efficient transcription on nucleosomal templates. Biochemical and genetic analyses, primarily performed in Arabidopsis, provided insight into the contribution of TEFs to establish gene expression patterns during plant growth and development. In addition to summarising the role of TEFs in plant gene expression, we emphasise in our review recent advances in the field. Thus, mechanisms are presented how aberrant intragenic transcript initiation is suppressed by repressing transcriptional start sites within coding sequences. We also discuss how transcriptional interference of ongoing transcription with neighbouring genes is prevented. Moreover, it appears that plants make no use of promoter-proximal RNAPII pausing in the way mammals do, but there are nucleosome-defined mechanism(s) that determine the efficiency of mRNA synthesis by RNAPII. Accordingly, a still growing number of processes related to plant growth, development and responses to changing environmental conditions prove to be regulated at the level of transcriptional elongation.

Genome-scale chromatin binding dynamics of RNA Polymerase II general transcription machinery components

A great deal of work has revealed, in structural detail, the components of the preinitiation complex (PIC) machinery required for initiation of mRNA gene transcription by RNA polymerase II (Pol II). However, less-well understood are the in vivo PIC assembly pathways and their kinetics, an understanding of which is vital for determining how rates of in vivo RNA synthesis are established. We used competition ChIP in budding yeast to obtain genome-scale estimates of the residence times for five general transcription factors (GTFs): TBP, TFIIA, TFIIB, TFIIE and TFIIF. While many GTF-chromatin interactions were short-lived ( < 1 min), there were numerous interactions with residence times in the range of several minutes. Sets of genes with a shared function also shared similar patterns of GTF kinetic behavior. TFIIE, a GTF that enters the PIC late in the assembly process, had residence times correlated with RNA synthesis rates. The datasets and results reported here provide kinetic information for most of the Pol II-driven genes in this organism, offering a rich resource for exploring the mechanistic relationships between PIC assembly, gene regulation, and transcription.

Multiple direct and indirect roles of Paf1C in elongation, splicing, and histone post-translational modifications

Paf1C is a highly conserved protein complex with critical functions during eukaryotic transcription. Previous studies have shown that Paf1C is multi-functional, controlling specific aspects of transcription, ranging from RNAPII processivity to histone modifications. However, it is unclear how specific Paf1C subunits directly impact transcription and coupled processes. We have compared conditional depletion to steady-state deletion for each Paf1C subunit to determine the direct and indirect contributions to gene expression in Saccharomyces cerevisiae. Using nascent transcript sequencing, RNAPII profiling, and modeling of transcription elongation dynamics, we have demonstrated direct effects of Paf1C subunits on RNAPII processivity and elongation rate and indirect effects on transcript splicing and repression of antisense transcripts. Further, our results suggest that the direct transcriptional effects of Paf1C cannot be readily assigned to any particular histone modification. This work comprehensively analyzes both the immediate and extended roles of each Paf1C subunit in transcription elongation and transcript regulation.

Yeast heterochromatin stably silences only weak regulatory elements by altering burst duration

Transcriptional silencing in Saccharomyces cerevisiae involves the generation of a chromatin state that stably represses transcription. Using multiple reporter assays, a diverse set of upstream activating sequence enhancers and core promoters were investigated for their susceptibility to silencing. We show that heterochromatin stably silences only weak and stress-induced regulatory elements but is unable to stably repress housekeeping gene regulatory elements, and the partial repression of these elements did not result in bistable expression states. Permutation analysis of enhancers and promoters indicates that both elements are targets of repression. Chromatin remodelers help specific regulatory elements to resist repression, most probably by altering nucleosome mobility and changing transcription burst duration. The strong enhancers/promoters can be repressed if silencer-bound Sir1 is increased. Together, our data suggest that the heterochromatic locus has been optimized to stably silence the weak mating-type gene regulatory elements but not strong housekeeping gene regulatory sequences.

Species-aware DNA language models capture regulatory elements and their evolution

BACKGROUND

The rise of large-scale multi-species genome sequencing projects promises to shed new light on how genomes encode gene regulatory instructions. To this end, new algorithms are needed that can leverage conservation to capture regulatory elements while accounting for their evolution.

RESULTS

Here, we introduce species-aware DNA language models, which we trained on more than 800 species spanning over 500 million years of evolution. Investigating their ability to predict masked nucleotides from context, we show that DNA language models distinguish transcription factor and RNA-binding protein motifs from background non-coding sequence. Owing to their flexibility, DNA language models capture conserved regulatory elements over much further evolutionary distances than sequence alignment would allow. Remarkably, DNA language models reconstruct motif instances bound in vivo better than unbound ones and account for the evolution of motif sequences and their positional constraints, showing that these models capture functional high-order sequence and evolutionary context. We further show that species-aware training yields improved sequence representations for endogenous and MPRA-based gene expression prediction, as well as motif discovery.

CONCLUSIONS

Collectively, these results demonstrate that species-aware DNA language models are a powerful, flexible, and scalable tool to integrate information from large compendia of highly diverged genomes.

The rise of single-cell transcriptomics in yeast

The field of single-cell omics has transformed our understanding of biological processes and is constantly advancing both experimentally and computationally. One of the most significant developments is the ability to measure the transcriptome of individual cells by single-cell RNA-seq (scRNA-seq), which was pioneered in higher eukaryotes. While yeast has served as a powerful model organism in which to test and develop transcriptomic technologies, the implementation of scRNA-seq has been significantly delayed in this organism, mainly because of technical constraints associated with its intrinsic characteristics, namely the presence of a cell wall, a small cell size and little amounts of RNA. In this review, we examine the current technologies for scRNA-seq in yeast and highlight their strengths and weaknesses. Additionally, we explore opportunities for developing novel technologies and the potential outcomes of implementing single-cell transcriptomics and extension to other modalities. Undoubtedly, scRNA-seq will be invaluable for both basic and applied yeast research, providing unique insights into fundamental biological processes.

Size fractionated NET-Seq reveals a conserved architecture of transcription units around yeast genes

Genomes from yeast to humans are subject to pervasive transcription. A single round of pervasive transcription is sufficient to alter local chromatin conformation, nucleosome dynamics and gene expression, but is hard to distinguish from background signals. Size fractionated native elongating transcript sequencing (sfNET-Seq) was developed to precisely map nascent transcripts independent of expression levels. RNAPII-associated nascent transcripts are fractionation into different size ranges before library construction. When anchored to the transcription start sites (TSS) of annotated genes, the combined pattern of the output metagenes gives the expected reference pattern. Bioinformatic pattern matching to the reference pattern identified 9542 transcription units in Saccharomyces cerevisiae, of which 47% are coding and 53% are noncoding. In total, 3113 (33%) are unannotated noncoding transcription units. Anchoring all transcription units to the TSS or polyadenylation site (PAS) of annotated genes reveals distinctive architectures of linked pairs of divergent transcripts approximately 200nt apart. The Reb1 transcription factor is enriched 30nt downstream of the PAS only when an upstream (TSS -60nt with respect to PAS) noncoding transcription unit co-occurs with a downstream (TSS +150nt) coding transcription unit and acts to limit levels of upstream antisense transcripts. The potential for extensive transcriptional interference is evident from low abundance unannotated transcription units with variable TSS (median -240nt) initiating within a 500nt window upstream of, and transcribing over, the promoters of protein-coding genes. This study confirms a highly interleaved yeast genome with different types of transcription units altering the chromatin landscape in distinctive ways, with the potential to exert extensive regulatory control.

Energy-driven genome regulation by ATP-dependent chromatin remodellers

The packaging of DNA into chromatin in eukaryotes regulates gene transcription, DNA replication and DNA repair. ATP-dependent chromatin remodelling enzymes (re)arrange nucleosomes at the first level of chromatin organization. Their Snf2-type motor ATPases alter histone-DNA interactions through a common DNA translocation mechanism. Whether remodeller activities mainly catalyse nucleosome dynamics or accurately co-determine nucleosome organization remained unclear. In this Review, we discuss the emerging mechanisms of chromatin remodelling: dynamic remodeller architectures and their interactions, the inner workings of the ATPase cycle, allosteric regulation and pathological dysregulation. Recent mechanistic insights argue for a decisive role of remodellers in the energy-driven self-organization of chromatin, which enables both stability and plasticity of genome regulation - for example, during development and stress. Different remodellers, such as members of the SWI/SNF, ISWI, CHD and INO80 families, process (epi)genetic information through specific mechanisms into distinct functional outputs. Combinatorial assembly of remodellers and their interplay with histone modifications, histone variants, DNA sequence or DNA-bound transcription factors regulate nucleosome mobilization or eviction or histone exchange. Such input-output relationships determine specific nucleosome positions and compositions with distinct DNA accessibilities and mediate differential genome regulation. Finally, remodeller genes are often mutated in diseases characterized by genome dysregulation, notably in cancer, and we discuss their physiological relevance.

Transcription as source of genetic heterogeneity in budding yeast

Transcription presents challenges to genome stability both directly, by altering genome topology and exposing single-stranded DNA to chemical insults and nucleases, and indirectly by introducing obstacles to the DNA replication machinery. Such obstacles include the RNA polymerase holoenzyme itself, DNA-bound regulatory factors, G-quadruplexes and RNA-DNA hybrid structures known as R-loops. Here, we review the detrimental impacts of transcription on genome stability in budding yeast, as well as the mitigating effects of transcription-coupled nucleotide excision repair and of systems that maintain DNA replication fork processivity and integrity. Interactions between DNA replication and transcription have particular potential to induce mutation and structural variation, but we conclude that such interactions must have only minor effects on DNA replication by the replisome with little if any direct mutagenic outcome. However, transcription can significantly impair the fidelity of replication fork rescue mechanisms, particularly Break Induced Replication, which is used to restart collapsed replication forks when other means fail. This leads to de novo mutations, structural variation and extrachromosomal circular DNA formation that contribute to genetic heterogeneity, but only under particular conditions and in particular genetic contexts, ensuring that the bulk of the genome remains extremely stable despite the seemingly frequent interactions between transcription and DNA replication.

Structural perspectives on transcription in chromatin

In eukaryotes, all genetic processes take place in the cell nucleus, where DNA is packaged as chromatin in 'beads-on-a-string' nucleosome arrays. RNA polymerase II (RNAPII) transcribes protein-coding and many non-coding genes in this chromatin environment. RNAPII elongates RNA while passing through multiple nucleosomes and maintaining the integrity of the chromatin structure. Recent structural studies have shed light on the detailed mechanisms of this process, including how transcribing RNAPII progresses through a nucleosome and reassembles it afterwards, and how transcription elongation factors, chromatin remodelers, and histone chaperones participate in these processes. Other studies have also illuminated the crucial role of nucleosomes in preinitiation complex assembly and transcription initiation. In this review we outline these advances and discuss future perspectives.

In vitro reconstitution of chromatin domains shows a role for nucleosome positioning in 3D genome organization

Eukaryotic genomes are organized into chromatin domains. The molecular mechanisms driving the formation of these domains are difficult to dissect in vivo and remain poorly understood. Here we reconstitute Saccharomyces cerevisiae chromatin in vitro and determine its 3D organization at subnucleosome resolution by micrococcal nuclease-based chromosome conformation capture and molecular dynamics simulations. We show that regularly spaced and phased nucleosome arrays form chromatin domains in vitro that resemble domains in vivo. This demonstrates that neither loop extrusion nor transcription is required for basic domain formation in yeast. In addition, we find that the boundaries of reconstituted domains correspond to nucleosome-free regions and that insulation strength scales with their width. Finally, we show that domain compaction depends on nucleosome linker length, with longer linkers forming more compact structures. Together, our results demonstrate that regular nucleosome positioning is important for the formation of chromatin domains and provide a proof-of-principle for bottom-up 3D genome studies.

Identification of 1600 replication origins in S. cerevisiae

There are approximately 500 known origins of replication in the yeast genome, and the process by which DNA replication initiates at these locations is well understood. In particular, these sites are made competent to initiate replication by loading of the Mcm replicative helicase prior to the start of S phase; thus, 'a site that binds Mcm in G1' might be considered to provide an operational definition of a replication origin. By fusing a subunit of Mcm to micrococcal nuclease, we previously showed that known origins are typically bound by a single Mcm double hexamer, loaded adjacent to the ARS consensus sequence (ACS). Here, we extend this analysis from known origins to the entire genome, identifying candidate Mcm binding sites whose signal intensity varies over at least three orders of magnitude. Published data quantifying single-stranded DNA (ssDNA) during S phase revealed replication initiation among the most abundant 1600 of these sites, with replication activity decreasing with Mcm abundance and disappearing at the limit of detection of ssDNA. Three other hallmarks of replication origins were apparent among the most abundant 5500 sites. Specifically, these sites: (1) appeared in intergenic nucleosome-free regions flanked on one or both sides by well-positioned nucleosomes; (2) were flanked by ACSs; and (3) exhibited a pattern of GC skew characteristic of replication initiation. We conclude that, if sites at which Mcm double hexamers are loaded can function as replication origins, then DNA replication origins are at least threefold more abundant than previously assumed, and we suggest that replication may occasionally initiate in essentially every intergenic region. These results shed light on recent reports that as many as 15% of replication events initiate outside of known origins, and this broader distribution of replication origins suggest that S phase in yeast may be less distinct from that in humans than widely assumed.

Binding profiles for 954 Drosophila and C. elegans transcription factors reveal tissue specific regulatory relationships

A catalog of transcription factor (TF) binding sites in the genome is critical for deciphering regulatory relationships. Here we present the culmination of the modERN (model organism Encyclopedia of Regulatory Networks) consortium that systematically assayed TF binding events in vivo in two major model organisms, Drosophila melanogaster (fly) and Caenorhabditis elegans (worm). We describe key features of these datasets, comprising 604 TFs identifying 3.6M sites in the fly and 350 TFs identifying 0.9 M sites in the worm. Applying a machine learning model to these data identifies sets of TFs with a prominent role in promoting target gene expression in specific cell types. TF binding data are available through the ENCODE Data Coordinating Center and at https://epic.gs.washington.edu/modERNresource, which provides access to processed and summary data, as well as widgets to probe cell type-specific TF-target relationships. These data are a rich resource that should fuel investigations into TF function during development.

A single fiber view of the nucleosome organization in eukaryotic chromatin

Eukaryotic cells are thought to arrange nucleosomes into extended arrays with evenly spaced nucleosomes phased at genomic landmarks. Here we tested to what extent this stereotypic organization describes the nucleosome organization in Saccharomyces cerevisiae using Fiber-Seq, a long-read sequencing technique that maps entire nucleosome arrays on individual chromatin fibers in a high throughput manner. With each fiber coming from a different cell, Fiber-Seq uncovers cell-to-cell heterogeneity. The long reads reveal the nucleosome architecture even over repetitive DNA such as the ribosomal DNA repeats. The absolute nucleosome occupancy, a parameter that is difficult to obtain with conventional sequencing approaches, is a direct readout of Fiber-Seq. We document substantial deviations from the stereotypical nucleosome organization with unexpectedly long linker DNAs between nucleosomes, gene bodies missing entire nucleosomes, cell-to-cell heterogeneity in nucleosome occupancy, heterogeneous phasing of arrays and irregular nucleosome spacing. Nucleosome array structures are indistinguishable throughout the gene body and with respect to the direction of transcription arguing against transcription promoting array formation. Acute nucleosome depletion destroyed most of the array organization indicating that nucleosome remodelers cannot efficiently pack nucleosomes under those conditions. Given that nucleosomes are cis-regulatory elements, the cell-to-cell heterogeneity uncovered by Fiber-Seq provides much needed information to understand chromatin structure and function.

Protein interaction network revealed by quantitative proteomic analysis links TFIIB to multiple aspects of the transcription cycle

Although TFIIB is widely regarded as an initiation factor, recent reports have implicated it in multiple aspects of eukaryotic transcription. To investigate the broader role of TFIIB in transcription, we performed quantitative proteomic analysis of yeast TFIIB. We purified two different populations of TFIIB; one from soluble cell lysate, which is not engaged in transcription, and the other from the chromatin fraction which yields the transcriptionally active form of the protein. TFIIB purified from the chromatin exhibits several interactions that explain its non-canonical roles in transcription. RNAPII, TFIIF and TFIIH were the only components of the preinitiation complex with a significant presence in chromatin TFIIB. A notable feature was enrichment of all subunits of CF1 and Rat1 3' end processing-termination complexes in chromatin-TFIIB preparation. Subunits of the CPF termination complex were also detected in both chromatin and soluble derived TFIIB preparations. These results may explain the presence of TFIIB at the 3' end of genes during transcription as well as its role in promoter-termination interaction.

Activity of the pleiotropic drug resistance transcription factors Pdr1p and Pdr3p is modulated by binding site flanking sequences

The transcription factors Pdr1p and Pdr3p regulate pleiotropic drug resistance (PDR) in Saccharomyces cerevisiae via the PDR responsive elements (PDREs) to modulate gene expression. However, the exact mechanisms underlying the differences in their regulons remain unclear. Employing genomic occupancy profiling (CUT&RUN), binding assays, and transcription studies, we characterized the differences in sequence specificity between transcription factors. Findings reveal distinct preferences for core PDRE sequences and the flanking sequences for both proteins. While flanking sequences moderately alter DNA binding affinity, they significantly impact Pdr1/3p transcriptional activity. Notably, both proteins demonstrated the ability to bind half sites, showing potential enhancement of transcription from adjacent PDREs. This insight sheds light on ways Pdr1/3p can differentially regulate PDR.

GenoPipe: identifying the genotype of origin within (epi)genomic datasets

Confidence in experimental results is critical for discovery. As the scale of data generation in genomics has grown exponentially, experimental error has likely kept pace despite the best efforts of many laboratories. Technical mistakes can and do occur at nearly every stage of a genomics assay (i.e. cell line contamination, reagent swapping, tube mislabelling, etc.) and are often difficult to identify post-execution. However, the DNA sequenced in genomic experiments contains certain markers (e.g. indels) encoded within and can often be ascertained forensically from experimental datasets. We developed the Genotype validation Pipeline (GenoPipe), a suite of heuristic tools that operate together directly on raw and aligned sequencing data from individual high-throughput sequencing experiments to characterize the underlying genome of the source material. We demonstrate how GenoPipe validates and rescues erroneously annotated experiments by identifying unique markers inherent to an organism's genome (i.e. epitope insertions, gene deletions and SNPs).

Single molecule MATAC-seq reveals key determinants of DNA replication origin efficiency

Stochastic origin activation gives rise to significant cell-to-cell variability in the pattern of genome replication. The molecular basis for heterogeneity in efficiency and timing of individual origins is a long-standing question. Here, we developed Methylation Accessibility of TArgeted Chromatin domain Sequencing (MATAC-Seq) to determine single-molecule chromatin accessibility of four specific genomic loci. MATAC-Seq relies on preferential modification of accessible DNA by methyltransferases combined with Nanopore-Sequencing for direct readout of methylated DNA-bases. Applying MATAC-Seq to selected early-efficient and late-inefficient yeast replication origins revealed large heterogeneity of chromatin states. Disruption of INO80 or ISW2 chromatin remodeling complexes leads to changes at individual nucleosomal positions that correlate with changes in their replication efficiency. We found a chromatin state with an accessible nucleosome-free region in combination with well-positioned +1 and +2 nucleosomes as a strong predictor for efficient origin activation. Thus, MATAC-Seq identifies the large spectrum of alternative chromatin states that co-exist on a given locus previously masked in population-based experiments and provides a mechanistic basis for origin activation heterogeneity during eukaryotic DNA replication. Consequently, our single-molecule chromatin accessibility assay will be ideal to define single-molecule heterogeneity across many fundamental biological processes such as transcription, replication, or DNA repair in vitro and ex vivo.

Locus-specific proteome decoding reveals Fpt1 as a chromatin-associated negative regulator of RNA polymerase III assembly

Transcription of tRNA genes by RNA polymerase III (RNAPIII) is tuned by signaling cascades. The emerging notion of differential tRNA gene regulation implies the existence of additional regulatory mechanisms. However, tRNA gene-specific regulators have not been described. Decoding the local chromatin proteome of a native tRNA gene in yeast revealed reprogramming of the RNAPIII transcription machinery upon nutrient perturbation. Among the dynamic proteins, we identified Fpt1, a protein of unknown function that uniquely occupied RNAPIII-regulated genes. Fpt1 binding at tRNA genes correlated with the efficiency of RNAPIII eviction upon nutrient perturbation and required the transcription factors TFIIIB and TFIIIC but not RNAPIII. In the absence of Fpt1, eviction of RNAPIII was reduced, and the shutdown of ribosome biogenesis genes was impaired upon nutrient perturbation. Our findings provide support for a chromatin-associated mechanism required for RNAPIII eviction from tRNA genes and tuning the physiological response to changing metabolic demands.

Nucleic acid sequence contributes to remodeler-mediated targeting of histone H2A.Z

The variant histone H2A.Z is inserted into nucleosomes immediately downstream of promoters and is important for transcription. The site-specific deposition of H2A.Z is catalyzed by SWR, a conserved chromatin remodeler with affinity for promoter-proximal nucleosome depleted regions (NDRs) and histone acetylation. By comparing the genomic distribution of H2A.Z in wild-type and SWR-deficient cells, we found that SWR is also responsible for depositing H2A.Z at thousands of non-canonical sites not directly linked to NDRs or histone acetylation. To understand the targeting mechanism of H2A.Z, we presented SWR with a library of nucleosomes isolated from yeast and characterized those preferred by SWR. We found that SWR prefers nucleosomes associated with intergenic over coding regions, especially when polyadenine tracks are present. Insertion of polyadenine sequences into recombinant nucleosomes near the H2A-H2B binding site stimulated the H2A.Z insertion activity of SWR. Therefore, the genome is encoded with information contributing to remodeler-mediated targeting of H2A.Z.

Different elongation factors distinctly modulate RNA polymerase II transcription in Arabidopsis

Various transcript elongation factors (TEFs) including modulators of RNA polymerase II (RNAPII) activity and histone chaperones tune the efficiency of transcription in the chromatin context. TEFs are involved in establishing gene expression patterns during growth and development in Arabidopsis, while little is known about the genomic distribution of the TEFs and the way they facilitate transcription. We have mapped the genome-wide occupancy of the elongation factors SPT4-SPT5, PAF1C and FACT, relative to that of elongating RNAPII phosphorylated at residues S2/S5 within the carboxyterminal domain. The distribution of SPT4-SPT5 along transcribed regions closely resembles that of RNAPII-S2P, while the occupancy of FACT and PAF1C is rather related to that of RNAPII-S5P. Under transcriptionally challenging heat stress conditions, mutant plants lacking the corresponding TEFs are differentially impaired in transcript synthesis. Strikingly, in plants deficient in PAF1C, defects in transcription across intron/exon borders are observed that are cumulative along transcribed regions. Upstream of transcriptional start sites, the presence of FACT correlates with nucleosomal occupancy. Under stress conditions FACT is particularly required for transcriptional upregulation and to promote RNAPII transcription through +1 nucleosomes. Thus, Arabidopsis TEFs are differently distributed along transcribed regions, and are distinctly required during transcript elongation especially upon transcriptional reprogramming.

Identification of 1600 replication origins in S. cerevisiae

There are approximately 500 known origins of replication in the yeast genome, and the process by which DNA replication initiates at these locations is well understood. In particular, these sites are made competent to initiate replication by loading of the Mcm replicative helicase prior to the start of S phase; thus, "a site to which MCM is bound in G1" might be considered to provide an operational definition of a replication origin. By fusing a subunit of Mcm to micrococcal nuclease, a technique referred to as "Chromatin Endogenous Cleavage", we previously showed that known origins are typically bound by a single Mcm double hexamer, loaded adjacent to the ARS consensus sequence (ACS). Here we extend this analysis from known origins to the entire genome, identifying candidate Mcm binding sites whose signal intensity varies over at least 3 orders of magnitude. Published data quantifying the production of ssDNA during S phase showed clear evidence of replication initiation among the most abundant 1600 of these sites, with replication activity decreasing in concert with Mcm abundance and disappearing at the limit of detection of ssDNA. Three other hallmarks of replication origins were apparent among the most abundant 5,500 sites. Specifically, these sites (1) appeared in intergenic nucleosome-free regions that were flanked on one or both sides by well-positioned nucleosomes; (2) were flanked by ACSs; and (3) exhibited a pattern of GC skew characteristic of replication initiation. Furthermore, the high resolution of this technique allowed us to demonstrate a strong bias for detecting Mcm double-hexamers downstream rather than upstream of the ACS, which is consistent with the directionality of Mcm loading by Orc that has been observed in vitro. We conclude that, if sites at which Mcm double-hexamers are loaded can function as replication origins, then DNA replication origins are at least 3-fold more abundant than previously assumed, and we suggest that replication may occasionally initiate in essentially every intergenic region. These results shed light on recent reports that as many as 15% of replication events initiate outside of known origins, and this broader distribution of replication origins suggest that S phase in yeast may be less distinct from that in humans than is widely assumed.

Systematic Dissection of Sequence Features Affecting the Binding Specificity of a Pioneer Factor Reveals Binding Synergy Between FOXA1 and AP-1

Despite the unique ability of pioneer transcription factors (PFs) to target nucleosomal sites in closed chromatin, they only bind a small fraction of their genomic motifs. The underlying mechanism of this selectivity is not well understood. Here, we design a high-throughput assay called ChIP-ISO to systematically dissect sequence features affecting the binding specificity of a classic PF, FOXA1. Combining ChIP-ISO with in vitro and neural network analyses, we find that 1) FOXA1 binding is strongly affected by co-binding TFs AP-1 and CEBPB, 2) FOXA1 and AP-1 show binding cooperativity in vitro, 3) FOXA1's binding is determined more by local sequences than chromatin context, including eu-/heterochromatin, and 4) AP-1 is partially responsible for differential binding of FOXA1 in different cell types. Our study presents a framework for elucidating genetic rules underlying PF binding specificity and reveals a mechanism for context-specific regulation of its binding.

Hsp90 shapes adaptation by controlling the fitness consequences of regulatory variation

The essential stress-responsive chaperone Hsp90 impacts development and adaptation from microbes to humans. Yet despite evidence of its role in evolution, pathogenesis, and oncogenic transformation, the molecular mechanisms by which Hsp90 alters the consequences of mutations remain vigorously debated. Here we exploit the power of nucleotide-resolution genetic mapping in Saccharomyces cerevisiae to uncover more than 1,000 natural variant-to-phenotype associations governed by this molecular chaperone. Strikingly, Hsp90 more frequently modified the phenotypic effects of cis-regulatory variation than variants that altered protein sequence. Moreover, these interactions made the largest contribution to Hsp90-dependent heredity. Nearly all interacting variants-both regulatory and protein-coding-fell within clients of Hsp90 or targets of its direct binding partners. Hsp90 activity affected mutations in evolutionarily young genes, segmental deletions, and heterozygotes, highlighting its influence on variation central to evolutionary novelty. Reconciling the diverse epistatic effects of this chaperone, synthetic transcriptional regulation and reconstructions of natural alleles by genome editing revealed a central role for Hsp90 in regulating the fundamental relationship between activity and phenotype. Our findings establish that non-coding variation is a core driver of Hsp90's influence on heredity, offering a mechanistic explanation for the chaperone's strong effects on evolution and development across species.

Puzzling out nuclear pore complex assembly

Nuclear pore complexes (NPCs) are sophisticated multiprotein assemblies embedded within the nuclear envelope and controlling the exchanges of molecules between the cytoplasm and the nucleus. In this review, we summarize the mechanisms by which these elaborate complexes are built from their subunits, the nucleoporins, based on our ever-growing knowledge of NPC structural organization and on the recent identification of additional features of this process. We present the constraints faced during the production of nucleoporins, their gathering into oligomeric complexes, and the formation of NPCs within nuclear envelopes, and review the cellular strategies at play, from co-translational assembly to the enrolment of a panel of cofactors. Remarkably, the study of NPCs can inform our perception of the biogenesis of multiprotein complexes in general - and vice versa.

Transcription-coupled nucleosome assembly

Eukaryotic transcription occurs on chromatin, where RNA polymerase II encounters nucleosomes during elongation. These nucleosomes must unravel for the DNA to enter the active site. However, in most transcribed genes, nucleosomes remain intact due to transcription-coupled chromatin assembly mechanisms. These mechanisms primarily involve the local reassembly of displaced nucleosomes to prevent (epi)genomic instability and the emergence of cryptic transcription. As a fail-safe mechanism, cells can assemble nucleosomes de novo, particularly in highly transcribed genes, but this may result in the loss of epigenetic information. This review examines transcription-coupled chromatin assembly, with an emphasis on studies in yeast and recent structural studies. These studies shed light on how elongation factors and histone chaperones coordinate to enable nucleosome recycling during transcription.

FUN-PROSE: A deep learning approach to predict condition-specific gene expression in fungi

mRNA levels of all genes in a genome is a critical piece of information defining the overall state of the cell in a given environmental condition. Being able to reconstruct such condition-specific expression in fungal genomes is particularly important to metabolically engineer these organisms to produce desired chemicals in industrially scalable conditions. Most previous deep learning approaches focused on predicting the average expression levels of a gene based on its promoter sequence, ignoring its variation across different conditions. Here we present FUN-PROSE-a deep learning model trained to predict differential expression of individual genes across various conditions using their promoter sequences and expression levels of all transcription factors. We train and test our model on three fungal species and get the correlation between predicted and observed condition-specific gene expression as high as 0.85. We then interpret our model to extract promoter sequence motifs responsible for variable expression of individual genes. We also carried out input feature importance analysis to connect individual transcription factors to their gene targets. A sizeable fraction of both sequence motifs and TF-gene interactions learned by our model agree with previously known biological information, while the rest corresponds to either novel biological facts or indirect correlations.

The Swi-Snf chromatin remodeling complex mediates gene repression through metabolic control

In eukaryotes, ATP-dependent chromatin remodelers regulate gene expression in response to nutritional and metabolic stimuli. However, altered transcription of metabolic genes may have significant indirect consequences which are currently poorly understood. In this study, we use genetic and molecular approaches to uncover a role for the remodeler Swi-Snf as a critical regulator of metabolism. We find that snfΔ mutants display a cysteine-deficient phenotype, despite growth in nutrient-rich media. This correlates with widespread perturbations in sulfur metabolic gene transcription, including global redistribution of the sulfur-sensing transcription factor Met4. Our findings show how a chromatin remodeler can have a significant impact on a whole metabolic pathway by directly regulating an important gene subset and demonstrate an emerging role for chromatin remodeling complexes as decisive factors in metabolic control.

Features of yeast RNA polymerase I with special consideration of the lobe binding subunits

Ribosomal RNAs (rRNAs) are structural components of ribosomes and represent the most abundant cellular RNA fraction. In the yeast Saccharomyces cerevisiae, they account for more than 60 % of the RNA content in a growing cell. The major amount of rRNA is synthesized by RNA polymerase I (Pol I). This enzyme transcribes exclusively the rRNA gene which is tandemly repeated in about 150 copies on chromosome XII. The high number of transcribed rRNA genes, the efficient recruitment of the transcription machinery and the dense packaging of elongating Pol I molecules on the gene ensure that enough rRNA is generated. Specific features of Pol I and of associated factors confer promoter selectivity and both elongation and termination competence. Many excellent reviews exist about the state of research about function and regulation of Pol I and how Pol I initiation complexes are assembled. In this report we focus on the Pol I specific lobe binding subunits which support efficient, error-free, and correctly terminated rRNA synthesis.

Synthesis of the ribosomal RNA precursor in human cells: mechanisms, factors and regulation

The ribosomal RNA precursor (pre-rRNA) comprises three of the four ribosomal RNAs and is synthesized by RNA polymerase (Pol) I. Here, we describe the mechanisms of Pol I transcription in human cells with a focus on recent insights gained from structure-function analyses. The comparison of Pol I-specific structural and functional features with those of other Pols and with the excessively studied yeast system distinguishes organism-specific from general traits. We explain the organization of the genomic rDNA loci in human cells, describe the Pol I transcription cycle regarding structural changes in the enzyme and the roles of human Pol I subunits, and depict human rDNA transcription factors and their function on a mechanistic level. We disentangle information gained by direct investigation from what had apparently been deduced from studies of the yeast enzymes. Finally, we provide information about how Pol I mutations may contribute to developmental diseases, and why Pol I is a target for new cancer treatment strategies, since increased rRNA synthesis was correlated with rapidly expanding cell populations.

Yeast Heterochromatin Only Stably Silences Weak Regulatory Elements by Altering Burst Duration

The interplay between nucleosomes and transcription factors leads to programs of gene expression. Transcriptional silencing involves the generation of a chromatin state that represses transcription and is faithfully propagated through DNA replication and cell division. Using multiple reporter assays, including directly visualizing transcription in single cells, we investigated a diverse set of UAS enhancers and core promoters for their susceptibility to heterochromatic gene silencing. These results show that heterochromatin only stably silences weak and stress induced regulatory elements but is unable to stably repress housekeeping gene regulatory elements and the partial repression did not result in bistable expression states. Permutation analysis of different UAS enhancers and core promoters indicate that both elements function together to determine the susceptibility of regulatory sequences to repression. Specific histone modifiers and chromatin remodellers function in an enhancer specific manner to aid these elements to resist repression suggesting that Sir proteins likely function in part by reducing nucleosome mobility. We also show that the strong housekeeping regulatory elements can be repressed if silencer bound Sir1 is increased, suggesting that Sir1 is a limiting component in silencing. Together, our data suggest that the heterochromatic locus has been optimized to stably silence the weak mating type gene regulatory elements but not strong housekeeping gene regulatory sequences which could help explain why these genes are often found at the boundaries of silenced domains.

The selfish yeast plasmid exploits a SWI/SNF-type chromatin remodeling complex for hitchhiking on chromosomes and ensuring high-fidelity propagation

Extra-chromosomal selfish DNA elements can evade the risk of being lost at every generation by behaving as chromosome appendages, thereby ensuring high fidelity segregation and stable persistence in host cell populations. The yeast 2-micron plasmid and episomes of the mammalian gammaherpes and papilloma viruses that tether to chromosomes and segregate by hitchhiking on them exemplify this strategy. We document for the first time the utilization of a SWI/SNF-type chromatin remodeling complex as a conduit for chromosome association by a selfish element. One principal mechanism for chromosome tethering by the 2-micron plasmid is the bridging interaction of the plasmid partitioning proteins (Rep1 and Rep2) with the yeast RSC2 complex and the plasmid partitioning locus STB. We substantiate this model by multiple lines of evidence derived from genomics, cell biology and interaction analyses. We describe a Rep-STB bypass system in which a plasmid engineered to non-covalently associate with the RSC complex mimics segregation by chromosome hitchhiking. Given the ubiquitous prevalence of SWI/SNF family chromatin remodeling complexes among eukaryotes, it is likely that the 2-micron plasmid paradigm or analogous ones will be encountered among other eukaryotic selfish elements.

The architecture of binding cooperativity between densely bound transcription factors

The binding of transcription factors (TFs) along genomes is restricted to a subset of sites containing their preferred motifs. TF-binding specificity is often attributed to the co-binding of interacting TFs; however, apart from specific examples, this model remains untested. Here, we define dependencies among budding yeast TFs that localize to overlapping promoters by profiling the genome-wide consequences of co-depleting multiple TFs. We describe unidirectional interactions, revealing Msn2 as a central factor allowing TF binding at its target promoters. By contrast, no case of mutual cooperation was observed. Particularly, Msn2 retained binding at its preferred promoters upon co-depletion of fourteen similarly bound TFs. Overall, the consequences of TF co-depletions were moderate, limited to a subset of promoters, and failed to explain the role of regions outside the DNA-binding domain in directing TF-binding preferences. Our results call for re-evaluating the role of cooperative interactions in directing TF-binding preferences.

Nucleosome retention by histone chaperones and remodelers occludes pervasive DNA-protein binding

DNA packaging within chromatin depends on histone chaperones and remodelers that form and position nucleosomes. Cells express multiple such chromatin regulators with overlapping in-vitro activities. Defining specific in-vivo activities requires monitoring histone dynamics during regulator depletion, which has been technically challenging. We have recently generated histone-exchange sensors in Saccharomyces cerevisiae, which we now use to define the contributions of 15 regulators to histone dynamics genome-wide. While replication-independent exchange in unperturbed cells maps to promoters, regulator depletions primarily affected gene bodies. Depletion of Spt6, Spt16 or Chd1 sharply increased nucleosome replacement sequentially at the beginning, middle or end of highly expressed gene bodies. They further triggered re-localization of chaperones to affected gene body regions, which compensated for nucleosome loss during transcription complex passage, but concurred with extensive TF binding in gene bodies. We provide a unified quantitative screen highlighting regulator roles in retaining nucleosome binding during transcription and preserving genomic packaging.