YeasTSS: an integrative web database of yeast transcription start sites

Tuning fungal promoters for the expression of eukaryotic proteins

Fungal systems, yeast as well as filamentous fungi, are effective platforms for producing recombinant eukaryotic proteins because of their efficient secretion, robust development features, and capacity for post-translational modification. However, to achieve optimum protein expression in fungal hosts, a precise regulation of gene expression levels is necessary. Promoters are critical cis-regulatory regions that drive gene expression. Therefore, understanding the structure and function of fungal promoters and the factors that influence their performance is an essential step in developing yeast and filamentous fungal platforms as hosts for the expression and secretion of eukaryotic proteins. However, literature on the characterization of filamentous fungal promoters is non-exhaustive. The present review attempts to provide a comprehensive account of available information and future applications of fungal promoters. The properties of promoters from different classes of fungi are discussed with respect to their general structure, the core and proximal components that constitute the fungal promoters, types of fungal promoters based on their functions etc. Furthermore, the utility of fungal promoters for applications in healthcare, biofuels, agriculture and biotechnology are also discussed. The comprehensive understanding of fungal promoters will help in developing tailored promoters, paving the way for the optimum production of economically important eukaryotic proteins in different host organisms.

Yeast Crf1p is an activator with different roles in regulation of target genes

Under stress conditions, ribosome biogenesis is downregulated. This process requires that expression of ribosomal RNA, ribosomal protein, and ribosome biogenesis genes be controlled in a coordinated fashion. The mechanistic Target of Rapamycin Complex 1 (mTORC1) participates in sensing unfavorable conditions to effect the requisite change in gene expression. In Saccharomyces cerevisiae, downregulation of ribosomal protein genes involves dissociation of the activator Ifh1p in a process that depends on Utp22p, a protein that also functions in pre-rRNA processing. Ifh1p has a paralog, Crf1p, which was implicated in communicating mTORC1 inhibition and hence was perceived as a repressor. We focus here on two ribosomal biogenesis genes, encoding Utp22p and the high mobility group protein Hmo1p, both of which are required for communication of mTORC1 inhibition to target genes. Crf1p functions as an activator on these genes as evidenced by reduced mRNA abundance and RNA polymerase II occupancy in a crf1Δ strain. Inhibition of mTORC1 has distinct effects on expression of HMO1 and UTP22; for example, on UTP22, but not on HMO1, the presence of Crf1p promotes the stable depletion of Ifh1p. Our data suggest that Crf1p functions as a weak activator, and that it may be required to prevent re-binding of Ifh1p to some gene promoters after mTORC1 inhibition in situations when Ifh1p is available. We propose that the inclusion of genes encoding proteins required for mTORC1-mediated downregulation of ribosomal protein genes in the same regulatory circuit as the ribosomal protein genes serves to optimize transcriptional responses during mTORC1 inhibition.

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.

Synthetic reversed sequences reveal default genomic states

Pervasive transcriptional activity is observed across diverse species. The genomes of extant organisms have undergone billions of years of evolution, making it unclear whether these genomic activities represent effects of selection or 'noise'1-4. Characterizing default genome states could help understand whether pervasive transcriptional activity has biological meaning. Here we addressed this question by introducing a synthetic 101-kb locus into the genomes of Saccharomyces cerevisiae and Mus musculus and characterizing genomic activity. The locus was designed by reversing but not complementing human HPRT1, including its flanking regions, thus retaining basic features of the natural sequence but ablating evolved coding or regulatory information. We observed widespread activity of both reversed and native HPRT1 loci in yeast, despite the lack of evolved yeast promoters. By contrast, the reversed locus displayed no activity at all in mouse embryonic stem cells, and instead exhibited repressive chromatin signatures. The repressive signature was alleviated in a locus variant lacking CpG dinucleotides; nevertheless, this variant was also transcriptionally inactive. These results show that synthetic genomic sequences that lack coding information are active in yeast, but inactive in mouse embryonic stem cells, consistent with a major difference in 'default genomic states' between these two divergent eukaryotic cell types, with implications for understanding pervasive transcription, horizontal transfer of genetic information and the birth of new genes.

CRISPRi-induced transcriptional regulation of IAH1 gene and its influence on volatile compounds profile in Kluyveromyces marxianus DU3

Mezcal is a traditional Mexican distilled beverage, known for its marked organoleptic profile, which is influenced by several factors, such as the fermentation process, where a wide variety of microorganisms are present. Kluyveromyces marxianus is one of the main yeasts isolated from mezcal fermentations and has been associated with ester synthesis, contributing to the flavors and aromas of the beverage. In this study, we employed CRISPR interference (CRISPRi) technology, using dCas9 fused to the Mxi1 repressor factor domain, to down-regulate the expression of the IAH1 gene, encoding for an isoamyl acetate-hydrolyzing esterase, in K. marxianus strain DU3. The constructed CRISPRi plasmid successfully targeted the IAH1 gene, allowing for specific gene expression modulation. Through gene expression analysis, we assessed the impact of IAH1 down-regulation on the metabolic profile of volatile compounds. We also measured the expression of other genes involved in volatile compound biosynthesis, including ATF1, EAT1, ADH1, and ZWF1 by RT-qPCR. Results demonstrated successful down-regulation of IAH1 expression in K. marxianus strain DU3 using the CRISPRi system. The modulation of IAH1 gene expression resulted in alterations in the production of volatile compounds, specifically ethyl acetate, which are important contributors to the beverage's aroma. Changes in the expression levels of other genes involved in ester biosynthesis, suggesting that the knockdown of IAH1 may generate intracellular alterations in the balance of these metabolites, triggering a regulatory response. The application of CRISPRi technology in K. marxianus opens the possibility of targeted modulation of gene expression, metabolic engineering strategies, and synthetic biology in this yeast strain.

Genome Editing, Transcriptional Regulation, and Forward Genetic Screening Using CRISPR-Cas12a Systems in Yarrowia lipolytica

Class II Type V endonucleases have increasingly been adapted to develop sophisticated and easily accessible synthetic biology tools for genome editing, transcriptional regulation, and functional genomic screening in a wide range of organisms. One such endonuclease, Cas12a, presents itself as an attractive alternative to Cas9-based systems. The ability to mature its own guide RNAs (gRNAs) from a single transcript has been leveraged for easy multiplexing, and its lack of requirement of a tracrRNA element, also allows for short gRNA expression cassettes. To extend these functionalities into the industrially relevant oleaginous yeast Yarrowia lipolytica, we developed a set of CRISPR-Cas12a vectors for easy multiplexed gene knockout, repression, and activation. We further extended the utility of this CRISPR-Cas12a system to functional genomic screening by constructing a genome-wide guide library targeting every gene with an eightfold coverage. Pooled CRISPR screens conducted with this library were used to profile Cas12a guide activities and develop a machine learning algorithm that could accurately predict highly efficient Cas12a gRNA. In this protocols chapter, we first present a method by which protein coding genes may be functionally disrupted via indel formation with CRISPR-Cas12a systems. Further, we describe how Cas12a fused to a transcriptional regulator can be used in conjunction with shortened gRNA to achieve transcriptional repression or activation. Finally, we describe the design, cloning, and validation of a genome-wide library as well as a protocol for the execution of a pooled CRISPR screen, to determine guide activity profiles in a genome-wide context in Y. lipolytica. The tools and strategies discussed here expand the list of available synthetic biology tools for facile genome engineering in this industrially important host.

ChEC-seq2: an improved Chromatin Endogenous Cleavage sequencing method and bioinformatic analysis pipeline for mapping in vivo protein-DNA interactions

Defining the in vivo DNA binding specificity of transcription factors (TFs) has relied nearly exclusively on chromatin immunoprecipitation (ChIP). While ChIP reveals TF binding patterns, its resolution is low. Higher resolution methods employing nucleases such as ChIP-exo, chromatin endogenous cleavage (ChEC-seq) and CUT&RUN resolve both TF occupancy and binding site protection. ChEC-seq, in which an endogenous TF is fused to micrococcal nuclease, requires neither fixation nor antibodies. However, the specificity of DNA cleavage during ChEC has been suggested to be lower than the specificity of the peaks identified by ChIP or ChIP-exo, perhaps reflecting non-specific binding of transcription factors to DNA. We have simplified the ChEC-seq protocol to minimize nuclease digestion while increasing the yield of cleaved DNA. ChEC-seq2 cleavage patterns were highly reproducible between replicates and with published ChEC-seq data. Combined with DoubleChEC, a new bioinformatic pipeline that removes non-specific cleavage sites, ChEC-seq2 identified high-confidence cleavage sites for three different yeast TFs that are strongly enriched for their known binding sites and adjacent to known target genes.

Molecular dynamics simulation-based trinucleotide and tetranucleotide level structural and energy characterization of the functional units of genomic DNA

Genomes of most organisms on earth are written in a universal language of life, made up of four units - adenine (A), thymine (T), guanine (G), and cytosine (C), and understanding the way they are put together has been a great challenge to date. Multiple efforts have been made to annotate this wonderfully engineered string of DNA using different methods but they lack a universal character. In this article, we have investigated the structural and energetic profiles of both prokaryotes and eukaryotes by considering two essential genomic sites, viz., the transcription start sites (TSS) and exon-intron boundaries. We have characterized these sites by mapping the structural and energy features of DNA obtained from molecular dynamics simulations, which considers all possible trinucleotide and tetranucleotide steps. For DNA, these physicochemical properties show distinct signatures at the TSS and intron-exon boundaries. Our results firmly convey the idea that DNA uses the same dialect for prokaryotes and eukaryotes and that it is worth going beyond sequence-level analyses to physicochemical space to determine the functional destiny of DNA sequences.

Development and applications of a CRISPR activation system for facile genetic overexpression in Candida albicans

For the fungal pathogen Candida albicans, genetic overexpression readily occurs via a diversity of genomic alterations, such as aneuploidy and gain-of-function mutations, with important consequences for host adaptation, virulence, and evolution of antifungal drug resistance. Given the important role of overexpression on C. albicans biology, it is critical to develop and harness tools that enable the analysis of genes expressed at high levels in the fungal cell. Here, we describe the development, optimization, and application of a novel, single-plasmid-based CRISPR activation (CRISPRa) platform for targeted genetic overexpression in C. albicans, which employs a guide RNA to target an activator complex to the promoter region of a gene of interest, thus driving transcriptional expression of that gene. Using this system, we demonstrate the ability of CRISPRa to drive high levels of gene expression in C. albicans, and we assess optimal guide RNA targeting for robust and constitutive overexpression. We further demonstrate the specificity of the system via RNA sequencing. We highlight the application of CRISPR activation to overexpress genes involved in pathogenesis and drug susceptibility, and contribute toward the identification of novel phenotypes. Consequently, this tool will facilitate a broad range of applications for the study of C. albicans genetic overexpression.

New insights into genome annotation in Podospora anserina through re-exploiting multiple RNA-seq data

BACKGROUND

Publicly available RNA-seq datasets are often underused although being helpful to improve functional annotation of eukaryotic genomes. This is especially true for filamentous fungi genomes which structure differs from most well annotated yeast genomes. Podospora anserina is a filamentous fungal model, which genome has been sequenced and annotated in 2008. Still, the current annotation lacks information about cis-regulatory elements, including promoters, transcription starting sites and terminators, which are instrumental to integrate epigenomic features into global gene regulation strategies.

RESULTS

Here we took advantage of 37 RNA-seq experiments that were obtained in contrasted developmental and physiological conditions, to complete the functional annotation of P. anserina genome. Out of the 10,800 previously annotated genes, 5'UTR and 3'UTR were defined for 7554, among which, 3328 showed differential transcriptional signal starts and/or transcriptional end sites. In addition, alternative splicing events were detected for 2350 genes, mostly due alternative 3'splice sites and 1732 novel transcriptionally active regions (nTARs) in unannotated regions were identified.

CONCLUSIONS

Our study provides a comprehensive genome-wide functional annotation of P. anserina genome, including chromatin features, cis-acting elements such as UTRs, alternative splicing events and transcription of non-coding regions. These new findings will likely improve our understanding of gene regulation strategies in compact genomes, such as those of filamentous fungi. Characterization of alternative transcripts and nTARs paves the way to the discovery of putative new genes, alternative peptides or regulatory non-coding RNAs.

Applicability of the heterologous yeast promoters for recombinant protein production in Pichia pastoris

Promoter choice is an important step in recombinant protein production, which directly determines the expression manner as constitutive or inducible and the expression level of the recombinant protein. This study aims to investigate the applicability of heterologous yeast promoters (Kluyveromyces marxianus TPI, Hansenula polymorpha PMA, Candida tropicalis ICL, and Saccharomyces cerevisiae CUP) in Pichia pastoris. The regulation mode of the CtICL and ScCUP promoters in P. pastoris was found to be inducible and that of the KmTPI and HpPMA was constitutive. The carbon sources in which the promoters exhibited the highest activity were determined as glycerol for PMA and TPI, glucose for CUP, and ethanol for ICL. The DNA region showing the highest activity was determined as 1000 bp for all promoters by promoter deletion analysis. Results from the study demonstrate the potential of inducible and constitutive heterologous promoters allowing expression under different conditions in the P. pastoris expression system and offers alternatives to frequently used promoters. KEY POINTS: • Heterologous promoters exhibited similar expression pattern in P. pastoris with its native host. • HpPMA has the highest promoter activity among the heterologous promoters tested. • Reporter gene expression with ScCUP is responsive to elevating Cu²⁺in P. pastoris.

TSSr: an R package for comprehensive analyses of TSS sequencing data

Transcription initiation is regulated in a highly organized fashion to ensure proper cellular functions. Accurate identification of transcription start sites (TSSs) and quantitative characterization of transcription initiation activities are fundamental steps for studies of regulated transcriptions and core promoter structures. Several high-throughput techniques have been developed to sequence the very 5'end of RNA transcripts (TSS sequencing) on the genome scale. Bioinformatics tools are essential for processing, analysis, and visualization of TSS sequencing data. Here, we present TSSr, an R package that provides rich functions for mapping TSS and characterizations of structures and activities of core promoters based on all types of TSS sequencing data. Specifically, TSSr implements several newly developed algorithms for accurately identifying TSSs from mapped sequencing reads and inference of core promoters, which are a prerequisite for subsequent functional analyses of TSS data. Furthermore, TSSr also enables users to export various types of TSS data that can be visualized by genome browser for inspection of promoter activities in association with other genomic features, and to generate publication-ready TSS graphs. These user-friendly features could greatly facilitate studies of transcription initiation based on TSS sequencing data. The source code and detailed documentations of TSSr can be freely accessed at https://github.com/Linlab-slu/TSSr.

Efficient multiplexed gene regulation in Saccharomyces cerevisiae using dCas12a

CRISPR Cas12a is an RNA-programmable endonuclease particularly suitable for gene regulation. This is due to its preference for T-rich PAMs that allows it to more easily target AT-rich promoter sequences, and built-in RNase activity which can process a single CRISPR RNA array encoding multiple spacers into individual guide RNAs (gRNAs), thereby simplifying multiplexed gene regulation. Here, we develop a flexible dCas12a-based CRISPRi system for Saccharomyces cerevisiae and systematically evaluate its design features. This includes the role of the NLS position, use of repression domains, and the position of the gRNA target. Our optimal system is comprised of dCas12a E925A with a single C-terminal NLS and a Mxi1 or a MIG1 repression domain, which enables up to 97% downregulation of a reporter gene. We also extend this system to allow for inducible regulation via an RNAP II-controlled promoter, demonstrate position-dependent effects in crRNA arrays, and use multiplexed regulation to stringently control a heterologous β-carotene pathway. Together these findings offer valuable insights into the design constraints of dCas12a-based CRISPRi and enable new avenues for flexible and efficient gene regulation in S. cerevisiae.

Flexible analysis of TSS mapping data and detection of TSS shifts with TSRexploreR

Heterogeneity in transcription initiation has important consequences for transcript stability and translation, and shifts in transcription start site (TSS) usage are prevalent in various developmental, metabolic, and disease contexts. Accordingly, numerous methods for global TSS profiling have been developed, including most recently Survey of TRanscription Initiation at Promoter Elements with high-throughput sequencing (STRIPE-seq), a method to profile transcription start sites (TSSs) on a genome-wide scale with significant cost and time savings compared to previous methods. In anticipation of more widespread adoption of STRIPE-seq and related methods for construction of promoter atlases and studies of differential gene expression, we built TSRexploreR, an R package for end-to-end analysis of TSS mapping data. TSRexploreR provides functions for TSS and transcription start region (TSR) detection, normalization, correlation, visualization, and differential TSS/TSR analyses. TSRexploreR is highly interoperable, accepting the data structures of TSS and TSR sets generated by several existing tools for processing and alignment of TSS mapping data, such as CAGEr for Cap Analysis of Gene Expression (CAGE) data. Lastly, TSRexploreR implements a novel approach for the detection of shifts in TSS distribution.

Cell-to-cell heterogeneity of phosphate gene expression in yeast is controlled by alternative transcription, 14-3-3 and Spl2

Dependent on phosphate availability the yeast Saccharomyces cerevisiae expresses either low or high affinity phosphate transporters. In the presence of phosphate yeast cells still express low levels of the high affinity phosphate transporter Pho84. The regulator Spl2 is expressed in approximately 90% of the cells, and is not expressed in the remaining cells. Here we report that deletion of RRP6, encoding an exonuclease degrading non-coding RNA, or BMH1, encoding the major 14-3-3 isoform, resulted in less cells expressing SPL2 and in increased levels of RNA transcribed from sequences upstream of the SPL2 coding region. SPL2 stimulates its own expression and that of PHO84 ensuing a positive feedback. Upon deletion of the region responsible for upstream SPL2 transcription almost all cells express SPL2. These results indicate that the cell-to-cell variation in PHO84 and SPL2 expression is dependent on a specific part of the SPL2 promoter and is controlled by Bmh1 and Spl2.

Regulation of the conserved 3'-5' exoribonuclease EXOSC10/Rrp6 during cell division, development and cancer

The conserved 3'-5' exoribonuclease EXOSC10/Rrp6 processes and degrades RNA, regulates gene expression and participates in DNA double-strand break repair and control of telomere maintenance via degradation of the telomerase RNA component. EXOSC10/Rrp6 is part of the multimeric nuclear RNA exosome and interacts with numerous proteins. Previous clinical, genetic, biochemical and genomic studies revealed the protein's essential functions in cell division and differentiation, its RNA substrates and its relevance to autoimmune disorders and oncology. However, little is known about the regulatory mechanisms that control the transcription, translation and stability of EXOSC10/Rrp6 during cell growth, development and disease and how these mechanisms evolved from yeast to human. Herein, we provide an overview of the RNA- and protein expression profiles of EXOSC10/Rrp6 during cell division, development and nutritional stress, and we summarize interaction networks and post-translational modifications across species. Additionally, we discuss how known and predicted protein interactions and post-translational modifications influence the stability of EXOSC10/Rrp6. Finally, we explore the idea that different EXOSC10/Rrp6 alleles, which potentially alter cellular protein levels or affect protein function, might influence human development and disease progression. In this review we interpret information from the literature together with genomic data from knowledgebases to inspire future work on the regulation of this essential protein's stability in normal and malignant cells.

Genome-scale reconstruction of Gcn4/ATF4 networks driving a growth program

Growth and starvation are considered opposite ends of a spectrum. To sustain growth, cells use coordinated gene expression programs and manage biomolecule supply in order to match the demands of metabolism and translation. Global growth programs complement increased ribosomal biogenesis with sufficient carbon metabolism, amino acid and nucleotide biosynthesis. How these resources are collectively managed is a fundamental question. The role of the Gcn4/ATF4 transcription factor has been best studied in contexts where cells encounter amino acid starvation. However, high Gcn4 activity has been observed in contexts of rapid cell proliferation, and the roles of Gcn4 in such growth contexts are unclear. Here, using a methionine-induced growth program in yeast, we show that Gcn4/ATF4 is the fulcrum that maintains metabolic supply in order to sustain translation outputs. By integrating matched transcriptome and ChIP-Seq analysis, we decipher genome-wide direct and indirect roles for Gcn4 in this growth program. Genes that enable metabolic precursor biosynthesis indispensably require Gcn4; contrastingly ribosomal genes are partly repressed by Gcn4. Gcn4 directly binds promoter-regions and transcribes a subset of metabolic genes, particularly driving lysine and arginine biosynthesis. Gcn4 also globally represses lysine and arginine enriched transcripts, which include genes encoding the translation machinery. The Gcn4 dependent lysine and arginine supply thereby maintains the synthesis of the translation machinery. This is required to maintain translation capacity. Gcn4 consequently enables metabolic-precursor supply to bolster protein synthesis, and drive a growth program. Thus, we illustrate how growth and starvation outcomes are both controlled using the same Gcn4 transcriptional outputs that function in distinct contexts.

The origin and evolution of a distinct mechanism of transcription initiation in yeasts

The molecular process of transcription by RNA Polymerase II is highly conserved among eukaryotes (“classic model”). A distinct way of locating transcription start sites (TSSs) has been identified in a budding yeast Saccharomyces cerevisiae (“scanning model”). Herein, we applied genomic approaches to elucidate the origin of the scanning model and its underlying genetic mechanisms. We first identified TSSs at single-nucleotide resolution for 12 yeast species using the nAnT-iCAGE technique, which significantly improved the annotations of these genomes by providing accurate 5′ boundaries for protein-coding genes. We then inferred the initiation mechanism of each species based on its TSS maps and genome sequences. We discovered that the scanning model likely originated after the split of Yarrowia lipolytica and the other budding yeasts. Species that use the scanning model showed an adenine-rich region immediately upstream of the TSS that might facilitate TSS selection. Both initiation mechanisms share a strong preference for pyrimidine–purine dinucleotides surrounding the TSS. Our results suggest that the purine is required to accurately recruit the first nucleotide, thereby increasing the chances of a messenger RNA of being capped during mRNA maturation, which is critical for efficient translation initiation during protein biosynthesis. Based on our findings, we propose a model for TSS selection in the scanning-model species, as well as a model for the stepwise process responsible for the origin and evolution of the scanning model.

Extensive 5'-surveillance guards against non-canonical NAD-caps of nuclear mRNAs in yeast

The ubiquitous redox coenzyme nicotinamide adenine dinucleotide (NAD) acts as a non-canonical cap structure on prokaryotic and eukaryotic ribonucleic acids. Here we find that in budding yeast, NAD-RNAs are abundant (>1400 species), short (<170 nt), and mostly correspond to mRNA 5′-ends. The modification percentage of transcripts is low (<5%). NAD incorporation occurs mainly during transcription initiation by RNA polymerase II, which uses distinct promoters with a YAAG core motif for this purpose. Most NAD-RNAs are 3′-truncated. At least three decapping enzymes, Rai1, Dxo1, and Npy1, guard against NAD-RNA at different cellular locations, targeting overlapping transcript populations. NAD-mRNAs are not translatable in vitro. Our work indicates that in budding yeast, most of the NAD incorporation into RNA seems to be disadvantageous to the cell, which has evolved a diverse surveillance machinery to prematurely terminate, decap and reject NAD-RNAs.

NAD (nicotinamide adenine dinucleotide) acts as a non-canonical RNA cap structure in bacteria and eukaryotes. Here the authors demonstrate the whole landscape of budding yeast NAD-RNAs which are subject to diverse surveillance pathways, suggesting that NAD caps in budding yeast are mostly dysfunctional.

Parallel Concerted Evolution of Ribosomal Protein Genes in Fungi and Its Adaptive Significance

Ribosomal protein (RP) genes encode structural components of ribosomes, the cellular machinery for protein synthesis. A single functional copy has been maintained in most of 78–80 RP families in animals due to evolutionary constraints imposed by gene dosage balance. Some fungal species have maintained duplicate copies in most RP families. The mechanisms by which the RP genes were duplicated and maintained and their functional significance are poorly understood. To address these questions, we identified all RP genes from 295 fungi and inferred the timing and nature of gene duplication events for all RP families. We found that massive duplications of RP genes have independently occurred by different mechanisms in three distantly related lineages: budding yeasts, fission yeasts, and Mucoromycota. The RP gene duplicates in budding yeasts and Mucoromycota were mainly created by whole genome duplication events. However, duplicate RP genes in fission yeasts were likely generated by retroposition, which is unexpected considering their dosage sensitivity. The sequences of most RP paralogs have been homogenized by repeated gene conversion in each species, demonstrating parallel concerted evolution, which might have facilitated the retention of their duplicates. Transcriptomic data suggest that the duplication and retention of RP genes increased their transcript abundance. Physiological data indicate that increased ribosome biogenesis allowed these organisms to rapidly consume sugars through fermentation while maintaining high growth rates, providing selective advantages to these species in sugar-rich environments.

Transcriptome visualization and data availability at the Saccharomyces Genome Database

The Saccharomyces Genome Database (SGD; www.yeastgenome.org) maintains the official annotation of all genes in the Saccharomyces cerevisiae reference genome and aims to elucidate the function of these genes and their products by integrating manually curated experimental data. Technological advances have allowed researchers to profile RNA expression and identify transcripts at high resolution. These data can be configured in web-based genome browser applications for display to the general public. Accordingly, SGD has incorporated published transcript isoform data in our instance of JBrowse, a genome visualization platform. This resource will help clarify S. cerevisiae biological processes by furthering studies of transcriptional regulation, untranslated regions, genome engineering, and expression quantification in S. cerevisiae.

Guide RNA Engineering Enables Dual Purpose CRISPR-Cpf1 for Simultaneous Gene Editing and Gene Regulation in Yarrowia lipolytica

Yarrowia lipolytica has fast become a biotechnologically significant yeast for its ability to accumulate lipids to high levels. While there exists a suite of synthetic biology tools for genetic engineering in this yeast, there is a need for multipurposed tools for rapid strain generation. Here, we describe a dual purpose CRISPR-Cpf1 system that is capable of simultaneous gene disruption and gene regulation. Truncating guide RNA spacer length to 16 nt inhibited nuclease activity but not binding to the target loci, enabling gene activation and repression with Cpf1-fused transcriptional regulators. Gene repression was demonstrated using a Cpf1-Mxi1 fusion achieving a 7-fold reduction in mRNA, while CRISPR-activation with Cpf1-VPR increased hrGFP expression by 10-fold. High efficiency disruptions were achieved with gRNAs 23-25 bp in length, and efficiency and repression levels were maintained with multiplexed expression of truncated and full-length gRNAs. The developed CRISPR-Cpf1 system should prove useful in metabolic engineering, genome wide screening, and functional genomics studies.

Modulating transcription through development of semi-synthetic yeast core promoters

Altering gene expression regulation by promoter engineering is a very effective way to fine-tune heterologous pathways in eukaryotic hosts. Typically, pathway building approaches in yeast still use a limited set of long, native promoters. With the today’s introduction of longer and more complex pathways, an expansion of this synthetic biology toolbox is necessary. In this study we elucidated the core promoter structure of the well-characterized yeast TEF1 promoter and determined the minimal length needed for sufficient protein expression. Furthermore, this minimal core promoter sequence was used for the creation of a promoter library covering different expression strengths. This resulted in a group of short, 69 bp promoters with an 8.0-fold expression range. One exemplar had a two and four times higher expression compared to the native CYC1 and ADH1 promoter, respectively. Additionally, as it was described that the protein expression range could be broadened by upstream activating sequences (UASs), we integrated earlier described single and multiple short, synthetic UASs in front of the strongest yeast core promoter. This approach resulted to further variation in protein expression and an overall promoter library spanning a 20-fold activity range and covering a length from 69 bp to maximally 129 bp. Furthermore, the robustness of this library was assessed on three alternative carbon sources besides glucose. As such, the suitability of short yeast core promoters for metabolic engineering applications on different media, either in an individual context or combined with UAS elements, was demonstrated.

A nucleotide resolution map of Top2-linked DNA breaks in the yeast and human genome

DNA topoisomerases are required to resolve DNA topological stress. Despite this essential role, abortive topoisomerase activity generates aberrant protein-linked DNA breaks, jeopardising genome stability. Here, to understand the genomic distribution and mechanisms underpinning topoisomerase-induced DNA breaks, we map Top2 DNA cleavage with strand-specific nucleotide resolution across the S. cerevisiae and human genomes-and use the meiotic Spo11 protein to validate the broad applicability of this method to explore the role of diverse topoisomerase family members. Our data characterises Mre11-dependent repair in yeast and defines two strikingly different fractions of Top2 activity in humans: tightly localised CTCF-proximal, and broadly distributed transcription-proximal, the latter correlated with gene length and expression. Moreover, single nucleotide accuracy reveals the influence primary DNA sequence has upon Top2 cleavage-distinguishing sites likely to form canonical DNA double-strand breaks (DSBs) from those predisposed to form strand-biased DNA single-strand breaks (SSBs) induced by etoposide (VP16) in vivo.

Pervasive and dynamic transcription initiation in Saccharomyces cerevisiae

Transcription initiation is finely regulated to ensure proper expression and function of genes. The regulated transcription initiation in response to various environmental stimuli in a classic model organism Saccharomyces cerevisiae has not been systematically investigated. In this study, we generated quantitative maps of transcription start sites (TSSs) at a single-nucleotide resolution for S. cerevisiae grown in nine different conditions using no-amplification nontagging Cap analysis of gene expression (nAnT-iCAGE) sequencing. We mapped ∼1 million well-supported TSSs, suggesting highly pervasive transcription initiation in the compact genome of the budding yeast. The comprehensive TSS maps allowed us to identify core promoters for ∼96% verified protein-coding genes. We corrected misannotation of translation start codon for 122 genes and suggested an alternative start codon for 57 genes. We found that 56% of yeast genes are controlled by multiple core promoters, and alternative core promoter usage by a gene is widespread in response to changing environments. Most core promoter shifts are coupled with altered gene expression, indicating that alternative core promoter usage might play an important role in controlling gene transcriptional activities. Based on their activities in responding to environmental cues, we divided core promoters into constitutive class (55%) and inducible class (45%). The two classes of core promoters display distinctive patterns in transcriptional abundance, chromatin structure, promoter shape, and sequence context. In summary, our study improved the annotation of the yeast genome and demonstrated a much more pervasive and dynamic nature of transcription initiation in yeast than previously recognized.