Comprehensive profiling of the fission yeast transcription start site activity during stress and media response

Alternative TSS use is widespread in Cryptococcus fungi in response to environmental cues and regulated genome-wide by the transcription factor Tur1

Alternative transcription start site (TSS) usage regulation has been identified as a major means of gene expression regulation in metazoans. However, in fungi, its impact remains elusive as its study has thus far been restricted to model yeasts. Here, we first re-analyzed TSS-seq data to define genuine TSS clusters in 2 species of pathogenic Cryptococcus. We identified 2 types of TSS clusters associated with specific DNA sequence motifs. Our analysis also revealed that alternative TSS usage regulation in response to environmental cues is widespread in Cryptococcus, altering gene expression and protein targeting. Importantly, we performed a forward genetic screen to identify a unique transcription factor (TF) named Tur1, which regulates alternative TSS (altTSS) usage genome-wide when cells switch from exponential phase to stationary phase. ChiP-Seq and DamID-Seq analyses suggest that at some loci, the role of Tur1 might be direct. Tur1 has been previously shown to be essential for virulence in C. neoformans. We demonstrated here that a tur1Δ mutant strain is more sensitive to superoxide stress and phagocytosed more efficiently by macrophages than the wild-type (WT) strain.

mRNA cap-binding protein eIF4E1 is a novel regulator of Toxoplasma gondii latency

The protozoan parasite Toxoplasma gondii causes serious opportunistic disease due to its ability to persist in patients as latent tissue cysts. The molecular mechanisms coordinating conversion between proliferative parasites (tachyzoites) and latent cysts (bradyzoites) are not fully understood. We previously showed that phosphorylation of eIF2α accompanies bradyzoite formation, suggesting that this clinically relevant process involves regulation of mRNA translation. In this study, we investigated the composition and role of eIF4F multi-subunit complexes in translational control. Using CLIPseq, we find that the cap-binding subunit, eIF4E1, localizes to the 5'-end of all tachyzoite mRNAs, many of which show evidence of stemming from heterogeneous transcriptional start sites. We further show that eIF4E1 operates as the predominant cap-binding protein in two distinct eIF4F complexes. Using genetic and pharmacological approaches, we found that eIF4E1 deficiency triggers efficient spontaneous formation of bradyzoites without stress induction. Consistent with this result, we also show that stress-induced bradyzoites exhibit reduced eIF4E1 expression. Overall, our findings establish a novel role for eIF4F in translational control required for parasite latency and microbial persistence.

IMPORTANCE

Toxoplasma gondii is an opportunistic pathogen important to global human and animal health. There are currently no chemotherapies targeting the encysted form of the parasite. Consequently, a better understanding of the mechanisms controlling encystation is required. Here we show that the mRNA cap-binding protein, eIF4E1, regulates the encystation process. Encysted parasites reduce eIF4E1 levels, and depletion of eIF4E1 decreases the translation of ribosome-associated machinery and drives Toxoplasma encystation. Together, these data reveal a new layer of mRNA translational control that regulates parasite encystation and latency.

Factors governing the transcriptome changes and chronological lifespan of fission yeast during phosphate starvation

Starvation of Schizosaccharomyces pombe for inorganic phosphate elicits adaptive transcriptome changes in which mRNAs driving ribosome biogenesis, tRNA biogenesis, and translation are globally downregulated, while those for autophagy and phosphate mobilization are upregulated. Here, we interrogated three components of the starvation response: upregulated autophagy; the role of transcription factor Pho7 (an activator of the PHO regulon); and upregulated expression of ecl3, one of three paralogous genes (ecl1, ecl2, and ecl3) collectively implicated in cell survival during other nutrient stresses. Ablation of autophagy factor Atg1 resulted in early demise of phosphate-starved fission yeast, as did ablation of Pho7. Transcriptome profiling of phosphate-starved pho7Δ cells highlighted Pho7 as an activator of genes involved in phosphate acquisition and mobilization, not limited to the original three-gene PHO regulon, and additional starvation-induced genes (including ecl3) not connected to phosphate dynamics. Pho7-dependent gene induction during phosphate starvation tracked with the presence of Pho7 DNA-binding elements in the gene promoter regions. Fewer ribosome protein genes were downregulated in phosphate-starved pho7Δ cells versus WT, which might contribute to their shortened lifespan. An ecl3Δ mutant elicited no gene expression changes in phosphate-replete cells and had no impact on survival during phosphate starvation. By contrast, pan-ecl deletion (ecl123Δ) curtailed lifespan during chronic phosphate starvation. Phosphate-starved ecl123Δ cells experienced a more widespread downregulation of mRNAs encoding aminoacyl tRNA synthetases vis-à-vis WT or pho7Δ cells. Collectively, these results enhance our understanding of fission yeast phosphate homeostasis and survival during nutrient deprivation.

POMBOX: A Fission Yeast Cloning Toolkit for Molecular and Synthetic Biology

The fission yeast Schizosaccharomyces pombe is a popular model organism in molecular biology and cell physiology. With its ease of genetic manipulation and growth, supported by in-depth functional annotations in the PomBase database and genome-wide metabolic models,S. pombe is an attractive option for synthetic biology applications. However,S. pombe currently lacks modular tools for generating genetic circuits with more than 1 transcriptional unit. We developed a toolkit to address this gap. Adapted from the MoClo-YTK plasmid kit for Saccharomyces cerevisiae and using the same modular cloning grammar, our POMBOX toolkit is designed to facilitate fast, efficient, and modular construction of genetic circuits inS. pombe. It allows for interoperability when working with DNA sequences that are functional in bothS. cerevisiae and S. pombe (e.g., protein tags, antibiotic resistance cassettes, and coding sequences). Moreover, POMBOX enables the modular assembly of multigene pathways and increases the possible pathway length from 6 to 12 transcriptional units. We also adapted the stable integration vector homology arms to Golden Gate assembly and tested the genomic integration success rates depending on different sequence sizes, from 4 to 24 kb. We included 14 S. pombe promoters that we characterized using two fluorescent proteins, in both minimally defined (EMM2─Edinburgh minimal media) and complex (YES─yeast extract with supplements) media. Then, we examined the efficacy of 6 S. cerevisiae and 6 synthetic terminators in S. pombe. Finally, we used the POMBOX kit for a synthetic biology application in metabolic engineering and expressed plant enzymes in S. pombe to produce specialized metabolite precursors, namely, methylxanthine, amorpha-4,11-diene, and cinnamic acid from the purine, mevalonate, and aromatic amino acid pathways.

Functional Consequences of Shifting Transcript Boundaries in Glucose Starvation

Glucose is a major source of carbon and essential for the survival of many organisms, ranging from yeast to human. A sudden 60-fold reduction of glucose in exponentially growing fission yeast induces transcriptome-wide changes in gene expression. This regulation is multilayered, and the boundaries of transcripts are known to vary, with functional consequences at the protein level. By combining direct RNA sequencing with 5'-CAGE and short-read sequencing, we accurately defined the 5'- and 3'-ends of transcripts that are both poly(A) tailed and 5'-capped in glucose starvation, followed by proteome analysis. Our results confirm previous experimentally validated loci with alternative isoforms and reveal several transcriptome-wide patterns. First, we show that sense-antisense gene pairs are more strongly anticorrelated when a time lag is taken into account. Second, we show that the glucose starvation response initially elicits a shortening of 3'-UTRs and poly(A) tails, followed by a shortening of the 5'-UTRs at later time points. These result in domain gains and losses in proteins involved in the stress response. Finally, the relatively poor overlap both between differentially expressed genes (DEGs), differential transcript usage events (DTUs), and differentially detected proteins (DDPs) highlight the need for further study on post-transcriptional regulation mechanisms in glucose starvation.

Crosstalk between the tRNA methyltransferase Trm1 and RNA chaperone La influences eukaryotic tRNA maturation

tRNAs undergo an extensive maturation process involving posttranscriptional modifications often associated with tRNA structural stability and promoting the native fold. Impaired posttranscriptional modification has been linked to human disease, likely through defects in translation, mitochondrial function, and increased susceptibility to degradation by various tRNA decay pathways. More recently, evidence has emerged that bacterial tRNA modification enzymes can act as tRNA chaperones to guide tRNA folding in a manner independent from catalytic activity. Here, we provide evidence that the fission yeast tRNA methyltransferase Trm1, which dimethylates nuclear- and mitochondrial-encoded tRNAs at G26, can also promote tRNA functionality in the absence of catalysis. We show that WT and catalytic-dead Trm1 are active in an in vivo tRNA-mediated suppression assay and possess RNA strand annealing and dissociation activity in vitro, similar to previously characterized RNA chaperones. Trm1 and the RNA chaperone La have previously been proposed to function synergistically in promoting tRNA maturation, yet we surprisingly demonstrate that La binding to nascent pre-tRNAs decreases Trm1 tRNA dimethylation in vivo and in vitro. Collectively, these results support the hypothesis for tRNA modification enzymes that combine catalytic and noncatalytic activities to promote tRNA maturation, as well as expand our understanding of how La function can influence tRNA modification.

mRNA cap-binding protein eIF4E1 is a novel regulator of Toxoplasma gondii latency

The protozoan parasite Toxoplasma gondii causes serious opportunistic disease due to its ability to persist in patients as latent tissue cysts. The molecular mechanisms coordinating conversion between proliferative parasites (tachyzoites) and dormant cysts (bradyzoites) are not fully understood. We previously showed that phosphorylation of eIF2α accompanies bradyzoite formation, suggesting that this clinically relevant process involves regulation of mRNA translation. In this study, we investigated the composition and role of eIF4F multi-subunit complexes in translational control. Using CLIPseq, we find that the cap-binding subunit, eIF4E1, localizes to the 5'-end of all tachyzoite mRNAs, many of which show evidence of stemming from heterogenous transcriptional start sites. We further show that eIF4E1 operates as the predominant cap-binding protein in two distinct eIF4F complexes. Using genetic and pharmacological approaches, we found that eIF4E1 deficiency triggers efficient spontaneous formation of bradyzoites without stress induction. Consistent with this result, we also show that stress-induced bradyzoites exhibit reduced eIF4E1 expression. Overall, our findings establish a novel role for eIF4F in translational control required for parasite latency and microbial persistence.

The minimal intrinsic stochasticity of constitutively expressed eukaryotic genes is sub-Poissonian

Gene expression inherently gives rise to stochastic variation ("noise") in the production of gene products. Minimizing noise is crucial for ensuring reliable cellular functions. However, noise cannot be suppressed below a certain intrinsic limit. For constitutively expressed genes, this limit is typically assumed to be Poissonian noise, wherein the variance in mRNA numbers is equal to their mean. Here, we demonstrate that several cell division genes in fission yeast exhibit mRNA variances significantly below this limit. The reduced variance can be explained by a gene expression model incorporating multiple transcription and mRNA degradation steps. Notably, in this sub-Poissonian regime, distinct from Poissonian or super-Poissonian regimes, cytoplasmic noise is effectively suppressed through a higher mRNA export rate. Our findings redefine the lower limit of eukaryotic gene expression noise and uncover molecular requirements for achieving ultralow noise, which is expected to be important for vital cellular functions.

Cellular responses to long-term phosphate starvation of fission yeast: Maf1 determines fate choice between quiescence and death associated with aberrant tRNA biogenesis

Inorganic phosphate is an essential nutrient acquired by cells from their environment. Here, we characterize the adaptative responses of fission yeast to chronic phosphate starvation, during which cells enter a state of quiescence, initially fully reversible upon replenishing phosphate after 2 days but resulting in gradual loss of viability during 4 weeks of starvation. Time-resolved analyses of changes in mRNA levels revealed a coherent transcriptional program in which phosphate dynamics and autophagy were upregulated, while the machineries for rRNA synthesis and ribosome assembly, and for tRNA synthesis and maturation, were downregulated in tandem with global repression of genes encoding ribosomal proteins and translation factors. Consistent with the transcriptome changes, proteome analysis highlighted global depletion of 102 ribosomal proteins. Concomitant with this ribosomal protein deficit, 28S and 18S rRNAs became vulnerable to site-specific cleavages that generated temporally stable rRNA fragments. The finding that Maf1, a repressor of RNA polymerase III transcription, was upregulated during phosphate starvation prompted a hypothesis that its activity might prolong lifespan of the quiescent cells by limiting production of tRNAs. Indeed, we found that deletion of maf1 results in precocious death of phosphate-starved cells via a distinctive starvation-induced pathway associated with tRNA overproduction and dysfunctional tRNA biogenesis.

The minimal intrinsic stochasticity of constitutively expressed eukaryotic genes is sub-Poissonian

Stochastic variation in gene products ("noise") is an inescapable by-product of gene expression. Noise must be minimized to allow for the reliable execution of cellular functions. However, noise cannot be suppressed beyond an intrinsic lower limit. For constitutively expressed genes, this limit is believed to be Poissonian, meaning that the variance in mRNA numbers cannot be lower than their mean. Here, we show that several cell division genes in fission yeast have mRNA variances significantly below this limit, which cannot be explained by the classical gene expression model for low-noise genes. Our analysis reveals that multiple steps in both transcription and mRNA degradation are essential to explain this sub-Poissonian variance. The sub-Poissonian regime differs qualitatively from previously characterized noise regimes, a hallmark being that cytoplasmic noise is reduced when the mRNA export rate increases. Our study re-defines the lower limit of eukaryotic gene expression noise and identifies molecular requirements for ultra-low noise which are expected to support essential cell functions.

Comparative Research: Regulatory Mechanisms of Ribosomal Gene Transcription in Saccharomyces cerevisiae and Schizosaccharomyces pombe

Restricting ribosome biosynthesis and assembly in response to nutrient starvation is a universal phenomenon that enables cells to survive with limited intracellular resources. When cells experience starvation, nutrient signaling pathways, such as the target of rapamycin (TOR) and protein kinase A (PKA), become quiescent, leading to several transcription factors and histone modification enzymes cooperatively and rapidly repressing ribosomal genes. Fission yeast has factors for heterochromatin formation similar to mammalian cells, such as H3K9 methyltransferase and HP1 protein, which are absent in budding yeast. However, limited studies on heterochromatinization in ribosomal genes have been conducted on fission yeast. Herein, we shed light on and compare the regulatory mechanisms of ribosomal gene transcription in two species with the latest insights.

Alternative transcription start sites contribute to acute-stress-induced transcriptome response in human skeletal muscle

Background

More than half of human protein-coding genes have an alternative transcription start site (TSS). We aimed to investigate the contribution of alternative TSSs to the acute-stress-induced transcriptome response in human tissue (skeletal muscle) using the cap analysis of gene expression approach. TSSs were examined at baseline and during recovery after acute stress (a cycling exercise).

Results

We identified 44,680 CAGE TSS clusters (including 3764 first defined) belonging to 12,268 genes and annotated for the first time 290 TSSs belonging to 163 genes. The transcriptome dynamically changes during the first hours after acute stress; the change in the expression of 10% of genes was associated with the activation of alternative TSSs, indicating differential TSSs usage. The majority of the alternative TSSs do not increase proteome complexity suggesting that the function of thousands of alternative TSSs is associated with the fine regulation of mRNA isoform expression from a gene due to the transcription factor-specific activation of various alternative TSSs. We identified individual muscle promoter regions for each TSS using muscle open chromatin data (ATAC-seq and DNase-seq). Then, using the positional weight matrix approach we predicted time course activation of “classic” transcription factors involved in response of skeletal muscle to contractile activity, as well as diversity of less/un-investigated factors.

Conclusions

Transcriptome response induced by acute stress related to activation of the alternative TSSs indicates that differential TSSs usage is an essential mechanism of fine regulation of gene response to stress stimulus. A comprehensive resource of accurate TSSs and individual promoter regions for each TSS in muscle was created. This resource together with the positional weight matrix approach can be used to accurate prediction of TFs in any gene(s) of interest involved in the response to various stimuli, interventions or pathological conditions in human skeletal muscle.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40246-022-00399-8.

Proteasome-dependent truncation of the negative heterochromatin regulator Epe1 mediates antifungal resistance

Epe1 histone demethylase restricts H3K9-methylation-dependent heterochromatin, preventing it from spreading over, and silencing, gene-containing regions in fission yeast. External stress induces an adaptive response allowing heterochromatin island formation that confers resistance on surviving wild-type lineages. Here we investigate the mechanism by which Epe1 is regulated in response to stress. Exposure to caffeine or antifungals results in Epe1 ubiquitylation and proteasome-dependent removal of the N-terminal 150 residues from Epe1, generating truncated Epe1 (tEpe1) which accumulates in the cytoplasm. Constitutive tEpe1 expression increases H3K9 methylation over several chromosomal regions, reducing expression of underlying genes and enhancing resistance. Reciprocally, constitutive non-cleavable Epe1 expression decreases resistance. tEpe1-mediated resistance requires a functional JmjC demethylase domain. Moreover, caffeine-induced Epe1-to-tEpe1 cleavage is dependent on an intact cell integrity MAP kinase stress signaling pathway, mutations in which alter resistance. Thus, environmental changes elicit a mechanism that curtails the function of this key epigenetic modifier, allowing heterochromatin to reprogram gene expression, thereby bestowing resistance to some cells within a population. H3K9me-heterochromatin components are conserved in human and crop-plant fungal pathogens for which a limited number of antifungals exist. Our findings reveal how transient heterochromatin-dependent antifungal resistant epimutations develop and thus inform on how they might be countered.

Co-transcriptional RNA cleavage by Drosha homolog Pac1 triggers transcription termination in fission yeast

Transcription termination of protein-coding genes in eukaryotic cells usually relies on a tight coordination between the cleavage and polyadenylation of the pre-mRNA, and 5′-3′ degradation of the downstream nascent transcript. Here we investigated the contribution of the essential fission yeast endonuclease Pac1, a homolog of human Drosha that cleaves hairpin RNA structures, in triggering polyadenylation-independent transcription termination. Using ChIP-sequencing in Pac1-deficient cells, we found that Pac1 triggers transcription termination at snRNA and snoRNA genes as well as at specific protein-coding genes. Notably, we found that Pac1-dependent premature termination occurred at two genes encoding conserved transmembrane transporters whose expression were strongly repressed by Pac1. Analysis by genome editing indicated that a stem-loop structure in the nascent transcript directs Pac1-mediated cleavage and that the regions upstream and downstream of the Pac1 cleavage site in the targeted mRNAs were stabilized by mutation of nuclear 3′-5′ and 5′-3′ exonucleases, respectively. Our findings unveil a premature transcription termination pathway that uncouples co-transcriptional RNA cleavage from polyadenylation, triggering rapid nuclear RNA degradation.

Involvement of Pca1 in ROS-mediated apoptotic cell death induced by alpha-thujone in the fission yeast (Schizosaccharomyces pombe)

Alpha-thujone, widely used in beverages (1-5 mg/kg), is known to have cytotoxic effects, but the mode of action and the role of potential apoptotic proteins in yeast cell death should be unraveled. In this study, we used Schizosaccharomyces pombe, which is a promising unicellular model organism in mechanistic toxicology and cell biology, to investigate the involvement of pro-apoptotic factors in alpha-thujone-induced cell death. We showed alpha-thujone-induced ROS accumulation-dependent cytotoxicity and apoptosis. In addition, we used superoxide dismutase-deficient cells (sod1 and sod2 mutants) to understand the effect of oxidative stress. Alpha-thujone caused significant cytotoxicity and apoptotic cell death, particularly in sod mutants. Moreover, two potential apoptotic factors, pca1 and pnu1 (pombe caspase-1 and pombe nuc1) were investigated to understand which factor mediates alpha-thujone-induced cell death. Pca1-deficient cells showed increased survival rates and reduced apoptosis in comparison to parental cells after chemical treatment while pnu1 mutation did not cause any significant change and the response was found identical as of parental cells. Yeast responded to alpha-thujone in caspase-dependent manner which was very similar to that for acetic acid. In conclusion, alfa-thujone-induced apoptosis and accounting mechanisms, which were mediated by ROS and driven by Pca1, were clarified in the unicellular model, S. pombe.

The leucine-NH4+ uptake regulator Any1 limits growth as part of a general amino acid control response to loss of La protein by fission yeast

The sla1⁺ gene of Schizosachharoymces pombe encodes La protein which promotes proper processing of precursor-tRNAs. Deletion of sla1 (sla1Δ) leads to disrupted tRNA processing and sensitivity to target of rapamycin (TOR) inhibition. Consistent with this, media containing NH4⁺ inhibits leucine uptake and growth of sla1Δ cells. Here, transcriptome analysis reveals that genes upregulated in sla1Δ cells exhibit highly significant overalp with general amino acid control (GAAC) genes in relevant transcriptomes from other studies. Growth in NH4⁺ media leads to additional induced genes that are part of a core environmental stress response (CESR). The sla1Δ GAAC response adds to evidence linking tRNA homeostasis and broad signaling in S. pombe. We provide evidence that deletion of the Rrp6 subunit of the nuclear exosome selectively dampens a subset of GAAC genes in sla1Δ cells suggesting that nuclear surveillance-mediated signaling occurs in S. pombe. To study the NH4⁺-effects, we isolated sla1Δ spontaneous revertants (SSR) of the slow growth phenotype and found that GAAC gene expression and rapamycin hypersensitivity were also reversed. Genome sequencing identified a F32V substitution in Any1, a known negative regulator of NH4⁺-sensitive leucine uptake linked to TOR. We show that ³H-leucine uptake by SSR-any1-F32V cells in NH4⁺-media is more robust than by sla1Δ cells. Moreover, F32V may alter any1⁺ function in sla1Δ vs. sla1⁺ cells in a distinctive way. Thus deletion of La, a tRNA processing factor leads to a GAAC response involving reprogramming of amino acid metabolism, and isolation of the any1-F32V rescuing mutant provides an additional specific link.

Conferring thermotolerant phenotype to wild-type Yarrowia lipolytica improves cell growth and erythritol production

In microbial engineering, heat stress is an important environmental factor modulating cell growth, metabolic flux distribution and the synthesis of target products. Yarrowia lipolytica, as a GARS (generally recognized as safe) nonconventional yeast, has been widely used in the food industry, especially as the host of erythritol production. Biomanufacturing economics is limited by the high operational cost of cooling energy in large-scale fermentation. It is of great significance to select thermotolerant Y. lipolytica to reduce the cooling cost and elucidate the heat-resistant mechanism at molecular level. For this purpose, we performed adaptive evolution and obtained a thermotolerant strain named Y. lipolytica BBE-18. Transcriptome analysis allows us to identify four genes in thiamine metabolism pathway that are responsible for the complicated thermotolerant phenotype. The heat-resistant phenotype was validated with the model strain Y. lipolytica Po1f by overexpression of single and combined genes. Then, conferring the thermotolerant phenotype to the wild-type Y. lipolytica BBE-17 enable the strain to produce three-times more erythritol of the control strain with 3°C higher than optimal cultivation temperature. To our knowledge, this is the first report on engineering heat-resistant phenotype to improve the erythritol production in Y. lipolytica. However, due to the increase of culture temperature, a large amount of adenosine triphosphate is consumed to ensure the life activities of Y. lipolytica which limits the potential of cell synthetic products to a certain extent. Even so, this study provides a reference for Y. lipolytica to produce other products under high temperature.

Rapid and inexpensive preparation of genome-wide nucleosome footprints from model and non-model organisms

MNase-seq (micrococcal nuclease sequencing) is used to map nucleosome positions in eukaryotic genomes to study the relationship between chromatin structure and DNA-dependent processes. Current protocols require at least two days to isolate nucleosome-protected DNA fragments. We have developed a streamlined protocol for S. cerevisiae and other fungi which takes only three hours. Modified protocols were developed for wild fungi and mammalian cells. This method for rapidly producing sequencing-ready nucleosome footprints from several organisms makes MNase-seq faster and easier, with less chemical waste.

•

A fast way to prepare micrococcal nuclease nucleosome footprints for MNase-seq

•

Eliminates use of phenol and chloroform and reduces the amount of cells required

•

Adaptable for a variety of organisms

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.

The RNA exosome shapes the expression of key protein-coding genes

The ribonucleolytic exosome complex is central for nuclear RNA degradation, primarily targeting non-coding RNAs. Still, the nuclear exosome could have protein-coding (pc) gene-specific regulatory activities. By depleting an exosome core component, or components of exosome adaptor complexes, we identify ∼2900 transcription start sites (TSSs) from within pc genes that produce exosome-sensitive transcripts. At least 1000 of these overlap with annotated mRNA TSSs and a considerable portion of their transcripts share the annotated mRNA 3′ end. We identify two types of pc-genes, both employing a single, annotated TSS across cells, but the first type primarily produces full-length, exosome-sensitive transcripts, whereas the second primarily produces prematurely terminated transcripts. Genes within the former type often belong to immediate early response transcription factors, while genes within the latter are likely transcribed as a consequence of their proximity to upstream TSSs on the opposite strand. Conversely, when genes have multiple active TSSs, alternative TSSs that produce exosome-sensitive transcripts typically do not contribute substantially to overall gene expression, and most such transcripts are prematurely terminated. Our results display a complex landscape of sense transcription within pc-genes and imply a direct role for nuclear RNA turnover in the regulation of a subset of pc-genes.

Epigenetic gene silencing by heterochromatin primes fungal resistance

Heterochromatin that depends on histone H3 lysine 9 methylation (H3K9me) renders embedded genes transcriptionally silent^1-3. In the fission yeast Schizosaccharomyces pombe, H3K9me heterochromatin can be transmitted through cell division provided the counteracting demethylase Epe1 is absent^4,5. Heterochromatin heritability might allow wild-type cells under certain conditions to acquire epimutations, which could influence phenotype through unstable gene silencing rather than DNA change^6,7. Here we show that heterochromatin-dependent epimutants resistant to caffeine arise in fission yeast grown with threshold levels of caffeine. Isolates with unstable resistance have distinct heterochromatin islands with reduced expression of embedded genes, including some whose mutation confers caffeine resistance. Forced heterochromatin formation at implicated loci confirms that resistance results from heterochromatin-mediated silencing. Our analyses reveal that epigenetic processes promote phenotypic plasticity, letting wild-type cells adapt to unfavourable environments without genetic alteration. In some isolates, subsequent or coincident gene-amplification events augment resistance. Caffeine affects two anti-silencing factors: Epe1 is downregulated, reducing its chromatin association, and a shortened isoform of Mst2 histone acetyltransferase is expressed. Thus, heterochromatin-dependent epimutation provides a bet-hedging strategy allowing cells to adapt transiently to insults while remaining genetically wild type. Isolates with unstable caffeine resistance show cross-resistance to antifungal agents, suggesting that related heterochromatin-dependent processes may contribute to resistance of plant and human fungal pathogens to such agents.

Characterization of Arabidopsis thaliana Promoter Bidirectionality and Antisense RNAs by Inactivation of Nuclear RNA Decay Pathways

In animals, RNA polymerase II initiates transcription bidirectionally from gene promoters to produce pre-mRNAs on the forward strand and promoter upstream transcripts (PROMPTs) on the reverse strand. PROMPTs are degraded by the nuclear exosome. Previous studies based on nascent RNA approaches concluded that Arabidopsis (Arabidopsis thaliana) does not produce PROMPTs. Here, we used steady-state RNA sequencing in mutants defective in nuclear RNA decay including the exosome to reassess the existence of Arabidopsis PROMPTs. While they are rare, we identified ∼100 cases of exosome-sensitive PROMPTs in Arabidopsis. Such PROMPTs are sources of small interfering RNAs in exosome-deficient mutants, perhaps explaining why plants have evolved mechanisms to suppress PROMPTs. In addition, we found ∼200 long, unspliced and exosome-sensitive antisense RNAs that arise from transcription start sites within parts of the genome encoding 3'-untranslated regions on the sense strand. The previously characterized noncoding RNA that regulates expression of the key seed dormancy regulator, DELAY OF GERMINATION1, is a typical representative of this class of RNAs. Transcription factor genes are overrepresented among loci with exosome-sensitive antisense RNAs, suggesting a potential for widespread control of gene expression via this class of noncoding RNAs. Lastly, we assess the use of alternative promoters in Arabidopsis and compare the accuracy of existing TSS annotations.

Organismal benefits of transcription speed control at gene boundaries

RNA polymerase II (RNAPII) transcription is crucial for gene expression. RNAPII density peaks at gene boundaries, associating these key regions for gene expression control with limited RNAPII movement. The connections between RNAPII transcription speed and gene regulation in multicellular organisms are poorly understood. Here, we directly modulate RNAPII transcription speed by point mutations in the second largest subunit of RNAPII in Arabidopsis thaliana. A RNAPII mutation predicted to decelerate transcription is inviable, while accelerating RNAPII transcription confers phenotypes resembling auto-immunity. Nascent transcription profiling revealed that RNAPII complexes with accelerated transcription clear stalling sites at both gene ends, resulting in read-through transcription. The accelerated transcription mutant NRPB2-Y732F exhibits increased association with 5' splice site (5'SS) intermediates and enhanced splicing efficiency. Our findings highlight potential advantages of RNAPII stalling through local reduction in transcription speed to optimize gene expression for the development of multicellular organisms.

YeasTSS: an integrative web database of yeast transcription start sites

The transcription initiation landscape of eukaryotic genes is complex and highly dynamic. In eukaryotes, genes can generate multiple transcript variants that differ in 5' boundaries due to usages of alternative transcription start sites (TSSs), and the abundance of transcript isoforms are highly variable. Due to a large number and complexity of the TSSs, it is not feasible to depict details of transcript initiation landscape of all genes using text-format genome annotation files. Therefore, it is necessary to provide data visualization of TSSs to represent quantitative TSS maps and the core promoters (CPs). In addition, the selection and activity of TSSs are influenced by various factors, such as transcription factors, chromatin remodeling and histone modifications. Thus, integration and visualization of functional genomic data related to these features could provide a better understanding of the gene promoter architecture and regulatory mechanism of transcription initiation. Yeast species play important roles for the research and human society, yet no database provides visualization and integration of functional genomic data in yeast. Here, we generated quantitative TSS maps for 12 important yeast species, inferred their CPs and built a public database, YeasTSS (www.yeastss.org). YeasTSS was designed as a central portal for visualization and integration of the TSS maps, CPs and functional genomic data related to transcription initiation in yeast. YeasTSS is expected to benefit the research community and public education for improving genome annotation, studies of promoter structure, regulated control of transcription initiation and inferring gene regulatory network.

Diverse DNA Sequence Motifs Activate Meiotic Recombination Hotspots Through a Common Chromatin Remodeling Pathway

In meiosis, multiple different DNA sequence motifs help to position homologous recombination at hotspots in the genome. How do the seemingly disparate cis-acting regulatory modules each promote locally the activity of the basal recombination machinery? We defined molecular mechanisms of action for five different hotspot-activating DNA motifs (M26, CCAAT, Oligo-C, 4095, 4156) located independently at the same site within the ade6 locus of the fission yeast Schizosaccharomyces pombe Each motif promoted meiotic recombination (i.e., is active) within this context, and this activity required the respective binding proteins (transcription factors Atf1, Pcr1, Php2, Php3, Php5, Rst2). High-resolution analyses of chromatin structure by nucleosome scanning assays revealed that each motif triggers the displacement of nucleosomes surrounding the hotspot motif in meiosis. This chromatin remodeling required the respective sequence-specific binding proteins, was constitutive for two motifs, and was enhanced meiotically for three others. Hotspot activity of each motif strongly required the ATP-dependent chromatin remodeling enzyme Snf22 (Snf2/Swi2), with lesser dependence on Gcn5, Mst2, and Hrp3. These findings support a model in which most meiotic recombination hotspots are positioned by the binding of transcription factors to their respective DNA sites. The functional redundancy of multiple, sequence-specific protein-DNA complexes converges upon shared chromatin remodeling pathways that help provide the basal recombination machinery (Spo11/Rec12 complex) access to its DNA substrates within chromatin.

CAGEfightR: analysis of 5'-end data using R/Bioconductor

BACKGROUND

5'-end sequencing assays, and Cap Analysis of Gene Expression (CAGE) in particular, have been instrumental in studying transcriptional regulation. 5'-end methods provide genome-wide maps of transcription start sites (TSSs) with base pair resolution. Because active enhancers often feature bidirectional TSSs, such data can also be used to predict enhancer candidates. The current availability of mature and comprehensive computational tools for the analysis of 5'-end data is limited, preventing efficient analysis of new and existing 5'-end data.

RESULTS

We present CAGEfightR, a framework for analysis of CAGE and other 5'-end data implemented as an R/Bioconductor-package. CAGEfightR can import data from BigWig files and allows for fast and memory efficient prediction and analysis of TSSs and enhancers. Downstream analyses include quantification, normalization, annotation with transcript and gene models, TSS shape statistics, linking TSSs to enhancers via co-expression, identification of enhancer clusters, and genome-browser style visualization. While built to analyze CAGE data, we demonstrate the utility of CAGEfightR in analyzing nascent RNA 5'-data (PRO-Cap). CAGEfightR is implemented using standard Bioconductor classes, making it easy to learn, use and combine with other Bioconductor packages, for example popular differential expression tools such as limma, DESeq2 and edgeR.

CONCLUSIONS

CAGEfightR provides a single, scalable and easy-to-use framework for comprehensive downstream analysis of 5'-end data. CAGEfightR is designed to be interoperable with other Bioconductor packages, thereby unlocking hundreds of mature transcriptomic analysis tools for 5'-end data. CAGEfightR is freely available via Bioconductor: bioconductor.org/packages/CAGEfightR .

Studying fungal pathogens of humans and fungal infections: fungal diversity and diversity of approaches

Seminal work by Louis Pasteur revealed the contribution of fungi - yeasts and microsporidia to agroindustry and disease in animals, respectively. More than 150 years later, the impact of fungi on human health and beyond is an ever-increasing issue, although often underestimated. Recent studies estimate that fungal infections, especially those caused by Candida, Cryptococcus and Aspergillus species, kill more than one million people annually. Indeed, these neglected infections are in general very difficult to cure and the associated mortality remains very high even when antifungal treatments exist. The development of new antifungals and diagnostic tools that are both necessary to fight fungal diseases efficiently, requires greater insights in the biology of the fungal pathogens of humans in the context of the infection, on their epidemiology, and on their role in the human mycobiota. We also need a better understanding of the host immune responses to fungal pathogens as well as the genetic basis for the increased sensitivity of some individuals to fungal infections. Here, we highlight some recent progress made in these different areas of research, in particular based on work conducted in our own laboratories. These progresses should lay the ground for better management of fungal infections, as they provide opportunities for better diagnostic, vaccination, the development of classical antifungals but also strategies for targeting virulence factors or the host.

A step-by-step guide to analyzing CAGE data using R/Bioconductor

Cap Analysis of Gene Expression (CAGE) is one of the most popular 5'-end sequencing methods. In a single experiment, CAGE can be used to locate and quantify the expression of both Transcription Start Sites (TSSs) and enhancers. This is workflow is a case study on how to use the CAGEfightR package to orchestrate analysis of CAGE data within the Bioconductor project. This workflow starts from BigWig-files and covers both basic CAGE analyses such as identifying, quantifying and annotating TSSs and enhancers, advanced analysis such as finding interacting TSS-enhancer pairs and enhancer clusters, to differential expression analysis and alternative TSS usage. R-code, discussion and references are intertwined to help provide guidelines for future CAGE studies of the same kind.

Studying fungal pathogens of humans and fungal infections: fungal diversity and diversity of approaches

Seminal work by Louis Pasteur revealed the contribution of fungi-yeasts and microsporidia to agroindustry and disease in animals, respectively. More than 150 years later, the impact of fungi on human health and beyond is an ever-increasing issue, although often underestimated. Recent studies estimate that fungal infections, especially those caused by Candida, Cryptococcus and Aspergillus species, kill more than one million people annually. Indeed, these neglected infections are in general very difficult to cure and the associated mortality remains very high even when antifungal treatments exist. The development of new antifungals and diagnostic tools that are both necessary to fight fungal diseases efficiently, requires greater insights in the biology of the fungal pathogens of humans in the context of the infection, on their epidemiology, and on their role in the human mycobiota. We also need a better understanding of the host immune responses to fungal pathogens as well as the genetic basis for the increased sensitivity of some individuals to fungal infections. Here, we highlight some recent progress made in these different areas of research, in particular based on work conducted in our own laboratories. These progress should lay the ground for better management of fungal infections, as they provide opportunities for better diagnostic, vaccination, the development of classical antifungals but also strategies for targeting virulence factors or the host.

Transcription-driven chromatin repression of Intragenic transcription start sites

Progression of RNA polymerase II (RNAPII) transcription relies on the appropriately positioned activities of elongation factors. The resulting profile of factors and chromatin signatures along transcription units provides a “positional information system” for transcribing RNAPII. Here, we investigate a chromatin-based mechanism that suppresses intragenic initiation of RNAPII transcription. We demonstrate that RNAPII transcription across gene promoters represses their function in plants. This repression is characterized by reduced promoter-specific molecular signatures and increased molecular signatures associated with RNAPII elongation. The conserved FACT histone chaperone complex is required for this repression mechanism. Genome-wide Transcription Start Site (TSS) mapping reveals thousands of discrete intragenic TSS positions in fact mutants, including downstream promoters that initiate alternative transcript isoforms. We find that histone H3 lysine 4 mono-methylation (H3K4me1), an Arabidopsis RNAPII elongation signature, is enriched at FACT-repressed intragenic TSSs. Our analyses suggest that FACT is required to repress intragenic TSSs at positions that are in part characterized by elevated H3K4me1 levels. In sum, conserved and plant-specific chromatin features correlate with the co-transcriptional repression of intragenic TSSs. Our insights into TSS repression by RNAPII transcription promise to inform the regulation of alternative transcript isoforms and the characterization of gene regulation through the act of pervasive transcription across eukaryotic genomes.

Genes represent DNA elements that are transcribed into mRNA. However, the position where transcription actually starts can be dynamically regulated to expand the diversity of RNA isoforms produced from a single gene. Functionally, alternative Transcription Start Sites (TSSs) may generate protein isoforms with differing N-terminal regions and distinct cellular functions. In plants, light signaling regulates protein isoforms largely through regulated TSS selection, emphasizing the biological significance of this mechanism. Despite the importance of alternative TSS selection, little is known about the underlying molecular mechanisms. Here, we characterize for the first time how transcription initiation from an upstream promoter represses alternative downstream promoter activity in plants. This repression mechanism is associated with chromatin changes that are required to maintain precise gene expression control. Specific chromatin signatures are established during transcription via dynamic interactions between the transcription machinery and associated factors. The conserved histone chaperone complex FACT is one such factor involved in regulating the chromatin environment along genes during transcription. We find that mutant plants with reduced FACT activity specifically initiate transcription from thousands of intragenic positions, thus expanding RNA isoform diversity. Overall, our study reveals conserved and plant-specific chromatin features associated with the co-transcriptional repression of downstream intragenic TSSs. These findings promise to help inform the molecular mechanism underlying environmentally-triggered TSS regulation in plants.