Contributions of transcription and mRNA decay to gene expression dynamics of fission yeast in response to oxidative stress

Modeling of mRNA deadenylation rates reveal a complex relationship between mRNA deadenylation and decay

Complete cytoplasmic polyadenosine tail (polyA-tail) deadenylation is thought to be essential for initiating mRNA decapping and subsequent degradation. To investigate this prevalent model, we conducted direct RNA sequencing of S. cerevisiae mRNAs derived from chase experiments under steady-state and stress condition. Subsequently, we developed a numerical model based on a modified gamma distribution function, which estimated the transcriptomic deadenylation rate at 10 A/min. A simplified independent method, based on the delineation of quantile polyA-tail values, showed a correlation between the decay and deadenylation rates of individual mRNAs, which appeared consistent within functional transcript groups and associated with codon optimality. Notably, these rates varied during the stress response. Detailed analysis of ribosomal protein-coding mRNAs (RPG mRNAs), constituting 40% of the transcriptome, singled out this transcript group. While deadenylation and decay of RPG mRNAs accelerated under heat stress, their degradation could proceed even when deadenylation was blocked, depending entirely on ongoing nuclear export. Our findings support the general primary function of deadenylation in dictating the onset of decapping, while also demonstrating complex relations between these processes.

Mutations in yeast Pcf11, a conserved protein essential for mRNA 3' end processing and transcription termination, elicit the Environmental Stress Response

The Pcf11 protein is an essential subunit of the large complex that cleaves and polyadenylates eukaryotic mRNA precursor. It has also been functionally linked to gene-looping, termination of RNA Polymerase II (Pol II) transcripts, and mRNA export. We have examined a poorly characterized but conserved domain (amino acids 142-225) of the Saccharomyces cerevisiae  Pcf11 and found that while it is not needed for mRNA 3' end processing or termination downstream of the poly(A) sites of protein-coding genes, its presence improves the interaction with Pol II and the use of transcription terminators near gene promoters. Analysis of genome-wide Pol II occupancy in cells with Pcf11 missing this region, as well as Pcf11 mutated in the Pol II CTD Interacting Domain, indicates that systematic changes in mRNA expression are mediated primarily at the level of transcription. Global expression analysis also shows that a general stress response, involving both activation and suppression of specific gene sets known to be regulated in response to a wide variety of stresses, is induced in the two pcf11 mutants, even though cells are grown in optimal conditions. The mutants also cause an unbalanced expression of cell wall-related genes that does not activate the Cell Wall Integrity pathway but is associated with strong caffeine sensitivity. Based on these findings, we propose that Pcf11 can modulate the expression level of specific functional groups of genes in ways that do not involve its well-characterized role in mRNA 3' end processing.

The ubiquitin conjugase Rad6 mediates ribosome pausing during oxidative stress

Oxidative stress causes K63-linked ubiquitination of ribosomes by the E2 ubiquitin conjugase Rad6. How Rad6-mediated ubiquitination of ribosomes affects translation, however, is unclear. We therefore perform Ribo-seq and Disome-seq in Saccharomyces cerevisiae and show that oxidative stress causes ribosome pausing at specific amino acid motifs, which also leads to ribosome collisions. However, these redox-pausing signatures are lost in the absence of Rad6 and do not depend on the ribosome-associated quality control (RQC) pathway. We also show that Rad6 is needed to inhibit overall translation in response to oxidative stress and that its deletion leads to increased expression of antioxidant genes. Finally, we observe that the lack of Rad6 leads to changes during translation that affect activation of the integrated stress response (ISR) pathway. Our results provide a high-resolution picture of the gene expression changes during oxidative stress and unravel an additional stress response pathway affecting translation elongation.

Codon Usage and mRNA Stability are Translational Determinants of Cellular Response to Canonical Ferroptosis Inducers

Ferroptosis is a non-apoptotic cell death mechanism characterized by the generation of lipid peroxides. While many effectors in the ferroptosis pathway have been mapped, its epitranscriptional regulation is not yet fully understood. Ferroptosis can be induced via system xCT inhibition (Class I) or GPX4 inhibition (Class II). Previous works have revealed important differences in cellular response to different ferroptosis inducers. Importantly, blocking mRNA transcription or translation appears to protect cells against Class I ferroptosis inducing agents but not Class II. In this work, we examined the impact of blocking transcription (via Actinomycin D) or translation (via Cycloheximide) on Erastin (Class I) or RSL3 (Class II) induced ferroptosis. Blocking transcription or translation protected cells against Erastin but was detrimental against RSL3. Cycloheximide led to increased levels of GSH alone or when co-treated with Erastin via the activation of the reverse transsulfuration pathway. RNA sequencing analysis revealed early activation of a strong alternative splice program before observed changes in transcription. mRNA stability analysis revealed divergent mRNA stability changes in cellular response to Erastin or RSL3. Importantly, codon optimality biases were drastically different in either condition. Our data also implicated translation repression and rate as an important determinant of the cellular response to ferroptosis inducers. Given that mRNA stability and codon usage can be influenced via the tRNA epitranscriptome, we evaluated the role of a tRNA modifying enzyme in ferroptosis stress response. Alkbh1, a tRNA demethylase, led to translation repression and increased the resistance to Erastin but made cells more sensitive to RSL3.

The impact of bZIP Atf1ortholog global regulators in fungi

Regulation of signal transduction pathways is crucial for the maintenance of cellular homeostasis and organismal development in fungi. Transcription factors are key elements of this regulatory network. The basic-region leucine zipper (bZIP) domain of the bZIP-type transcription factors is responsible for DNA binding while their leucine zipper structural motifs are suitable for dimerization with each other facilitiating the formation of homodimeric or heterodimeric bZIP proteins. This review highlights recent knowledge on the function of fungal orthologs of the Schizosaccharomyces pombe Atf1, Aspergillus nidulans AtfA, and Fusarium verticillioides FvAtfA, bZIP-type transcription factors with a special focus on pathogenic species. We demonstrate that fungal Atf1-AtfA-FvAtfA orthologs play an important role in vegetative growth, sexual and asexual development, stress response, secondary metabolite production, and virulence both in human pathogens, including Aspergillus fumigatus, Mucor circinelloides, Penicillium marneffei, and Cryptococcus neoformans and plant pathogens, like Fusarium ssp., Magnaporthe oryzae, Claviceps purpurea, Botrytis cinerea, and Verticillium dahliae. KEY POINTS: • Atf1 orthologs play crucial role in the growth and development of fungi. • Atf1 orthologs orchestrate environmental stress response of fungi. • Secondary metabolite production and virulence are coordinated by Atf1 orthologs.

Genome analysis of Plectus murrayi, a nematode from continental Antarctica

Plectus murrayi is one of the most common and locally abundant invertebrates of continental Antarctic ecosystems. Because it is readily cultured on artificial medium in the laboratory and highly tolerant to an extremely harsh environment, P. murrayi is emerging as a model organism for understanding the evolutionary origin and maintenance of adaptive responses to multiple environmental stressors, including freezing and desiccation. The de novo assembled genome of P. murrayi contains 225.741 million base pairs and a total of 14,689 predicted genes. Compared to Caenorhabditis elegans, the architectural components of P. murrayi are characterized by a lower number of protein-coding genes, fewer transposable elements, but more exons, than closely related taxa from less harsh environments. We compared the transcriptomes of lab-reared P. murrayi with wild-caught P. murrayi and found genes involved in growth and cellular processing were up-regulated in lab-cultured P. murrayi, while a few genes associated with cellular metabolism and freeze tolerance were expressed at relatively lower levels. Preliminary comparative genomic and transcriptomic analyses suggest that the observed constraints on P. murrayi genome architecture and functional gene expression, including genome decay and intron retention, may be an adaptive response to persisting in a biotically simplified, yet consistently physically harsh environment.

The bZIP transcription factor VdAtf1 regulates virulence by mediating nitrogen metabolism in Verticillium dahliae

The fungus Verticillium dahliae causes vascular wilt disease on hundreds of plant species. Homologs of the bZIP transcription factor Atf1 are required for virulence in most pathogenic fungi, but the molecular basis for their involvement is largely unknown. We performed targeted gene deletion, expression analysis, biochemistry and pathogenicity assays to demonstrate that VdAtf1 governs pathogenesis via the regulation of nitrosative resistance and nitrogen metabolism in V. dahliae. VdAtf1 controls pathogenesis via the regulation of nitric oxide (NO) resistance and inorganic nitrogen metabolism rather than oxidative resistance and is important for penetration peg formation in V. dahliae. VdAtf1 affects ammonium and nitrate assimilation in response to various nitrogen sources. VdAtf1 may be involved in regulating the expression of VdNut1. VdAtf1 responds to NO stress by strengthening the fungal cell wall, and by causing over-accumulation of methylglyoxal and glycerol, which in turn impacts NO detoxification. We also verified that the VdAtf1 ortholog in Fusarium graminearum mediates nitrogen metabolism, suggesting conservation of this function in related plant pathogenic fungi. Our findings revealed new functions of VdAtf1 in pathogenesis, response to nitrosative stress and nitrogen metabolism in V. dahliae. The results provide novel insights into the regulatory mechanisms of the transcription factor VdAtf1 in virulence.

Prenatal Glucocorticoid Exposure Results in Changes in Gene Transcription and DNA Methylation in the Female Juvenile Guinea Pig Hippocampus Across Three Generations

Synthetic glucocorticoids (sGC) are administered to women at risk for pre-term delivery, to mature the fetal lung and decrease neonatal morbidity. sGC also profoundly affect the fetal brain. The hippocampus expresses high levels of glucocorticoid (GR) and mineralocorticoid receptor (MR), and its development is affected by elevated fetal glucocorticoid levels. Antenatal sGC results in neuroendocrine and behavioral changes that persist in three generations of female guinea pig offspring of the paternal lineage. We hypothesized that antenatal sGC results in transgenerational changes in gene expression that correlate with changes in DNA methylation. We used RNASeq and capture probe bisulfite sequencing to investigate the transcriptomic and epigenomic effects of antenatal sGC exposure in the hippocampus of three generations of juvenile female offspring from the paternal lineage. Antenatal sGC exposure (F0 pregnancy) resulted in generation-specific changes in hippocampal gene transcription and DNA methylation. Significant changes in individual CpG methylation occurred in RNApol II binding regions of small non-coding RNA (snRNA) genes, which implicates alternative splicing as a mechanism involved in transgenerational transmission of the effects of antenatal sGC. This study provides novel perspectives on the mechanisms involved in transgenerational transmission and highlights the importance of human studies to determine the longer-term effects of antenatal sGC on hippocampal-related function.

Local Changes in Chromatin Accessibility and Transcriptional Networks Underlying the Nitrate Response in Arabidopsis Roots

Transcriptional regulation, determined by the chromatin structure and regulatory elements interacting at promoter regions, is a key step in plant responses to environmental cues. Nitrate (NO₃^-) is a nutrient signal that regulates the expression of hundreds of genes in Arabidopsis thaliana. Here, we integrate mRNA sequencing, genome-wide RNA polymerase II (RNPII), chromatin immunoprecipitation sequencing, and DNase sequencing datasets to establish the relationship between RNPII occupancy and chromatin accessibility in response to NO₃^- treatments in Arabidopsis roots. Genomic footprinting allowed us to identify in vivo regulatory elements controlling gene expression in response to NO₃^- treatments. NO₃^--modulated transcription factor (TF) footprints are important for a rapid increase in RNPII occupancy and transcript accumulation over time. We mapped key TF regulatory interactions and functionally validated the role of NAP, an NAC-domain containing TF, as a new regulatory factor in NO₃^- transport. Taken together, our study provides a comprehensive view of transcriptional networks in response to a nutrient signal in Arabidopsis roots.

Functions for fission yeast splicing factors SpSlu7 and SpPrp18 in alternative splice-site choice and stress-specific regulated splicing

Budding yeast spliceosomal factors ScSlu7 and ScPrp18 interact and mediate intron 3'ss choice during second step pre-mRNA splicing. The fission yeast genome with abundant multi-intronic transcripts, degenerate splice signals and SR proteins is an apt unicellular fungal model to deduce roles for core spliceosomal factors in alternative splice-site choice, intron retention and to study the cellular implications of regulated splicing. From our custom microarray data we deduce a stringent reproducible subset of S. pombe alternative events. We examined the role of factors SpSlu7 or SpPrp18 for these splice events and investigated the relationship to growth phase and stress. Wild-type log and stationary phase cells showed ats1+ exon 3 skipped and intron 3 retained transcripts. Interestingly the non-consensus 5'ss in ats1+ intron 3 caused SpSlu7 and SpPrp18 dependent intron retention. We validated the use of an alternative 5'ss in dtd1+ intron 1 and of an upstream alternative 3'ss in DUF3074 intron 1. The dtd1+ intron 1 non-canonical 5'ss yielded an alternative mRNA whose levels increased in stationary phase. Utilization of dtd1+ intron 1 sub-optimal 5' ss required functional SpPrp18 and SpSlu7 while compromise in SpSlu7 function alone hampered the selection of the DUF3074 intron 1 non canonical 3'ss. We analysed the relative abundance of these splice isoforms during mild thermal, oxidative and heavy metal stress and found stress-specific splice patterns for ats1+ and DUF3074 intron 1 some of which were SpSlu7 and SpPrp18 dependent. By studying ats1+ splice isoforms during compromised transcription elongation rates in wild-type, spslu7-2 and spprp18-5 mutant cells we found dynamic and intron context-specific effects in splice-site choice. Our work thus shows the combinatorial effects of splice site strength, core splicing factor functions and transcription elongation kinetics to dictate alternative splice patterns which in turn serve as an additional recourse of gene regulation in fission yeast.

The non-stop decay mRNA surveillance pathway is required for oxidative stress tolerance

Reactive oxygen species (ROS) are toxic by-products of normal aerobic metabolism. ROS can damage mRNAs and the translational apparatus resulting in translational defects and aberrant protein production. Three mRNA quality control systems monitor mRNAs for translational errors: nonsense-mediated decay, non-stop decay (NSD) and no-go decay (NGD) pathways. Here, we show that factors required for the recognition of NSD substrates and components of the SKI complex are required for oxidant tolerance. We found an overlapping requirement for Ski7, which bridges the interaction between the SKI complex and the exosome, and NGD components (Dom34/Hbs1) which have been shown to function in both NSD and NGD. We show that ski7 dom34 and ski7 hbs1 mutants are sensitive to hydrogen peroxide stress and accumulate an NSD substrate. We further show that NSD substrates are generated during ROS exposure as a result of aggregation of the Sup35 translation termination factor, which increases stop codon read-through allowing ribosomes to translate into the 3΄-end of mRNAs. Overexpression of Sup35 decreases stop codon read-through and rescues oxidant tolerance consistent with this model. Our data reveal an unanticipated requirement for the NSD pathway during oxidative stress conditions which prevents the production of aberrant proteins from NSD mRNAs.

Cis-regulatory elements explain most of the mRNA stability variation across genes in yeast

The stability of mRNA is one of the major determinants of gene expression. Although a wealth of sequence elements regulating mRNA stability has been described, their quantitative contributions to half-life are unknown. Here, we built a quantitative model for Saccharomyces cerevisiae based on functional mRNA sequence features that explains 59% of the half-life variation between genes and predicts half-life at a median relative error of 30%. The model revealed a new destabilizing 3' UTR motif, ATATTC, which we functionally validated. Codon usage proves to be the major determinant of mRNA stability. Nonetheless, single-nucleotide variations have the largest effect when occurring on 3' UTR motifs or upstream AUGs. Analyzing mRNA half-life data of 34 knockout strains showed that the effect of codon usage not only requires functional decapping and deadenylation, but also the 5'-to-3' exonuclease Xrn1, the nonsense-mediated decay genes, but not no-go decay. Altogether, this study quantitatively delineates the contributions of mRNA sequence features on stability in yeast, reveals their functional dependencies on degradation pathways, and allows accurate prediction of half-life from mRNA sequence.

Kinetic CRAC uncovers a role for Nab3 in determining gene expression profiles during stress

RNA-binding proteins play a key role in shaping gene expression profiles during stress, however, little is known about the dynamic nature of these interactions and how this influences the kinetics of gene expression. To address this, we developed kinetic cross-linking and analysis of cDNAs (χCRAC), an ultraviolet cross-linking method that enabled us to quantitatively measure the dynamics of protein-RNA interactions in vivo on a minute time-scale. Here, using χCRAC we measure the global RNA-binding dynamics of the yeast transcription termination factor Nab3 in response to glucose starvation. These measurements reveal rapid changes in protein-RNA interactions within 1 min following stress imposition. Changes in Nab3 binding are largely independent of alterations in transcription rate during the early stages of stress response, indicating orthogonal transcriptional control mechanisms. We also uncover a function for Nab3 in dampening expression of stress-responsive genes. χCRAC has the potential to greatly enhance our understanding of in vivo dynamics of protein-RNA interactions.Protein RNA interactions are dynamic and regulated in response to environmental changes. Here the authors describe 'kinetic CRAC', an approach that allows time resolved analyses of protein RNA interactions with minute time point resolution and apply it to gain insight into the function of the RNA-binding protein Nab3.

Measures of RNA metabolism rates: Toward a definition at the level of single bonds

We give an overview of experimental and computational methods to estimate RNA metabolism rates genome-wide. We then advocate a local definition of RNA metabolism rate at the level of individual phosphodiester bonds. Rates of formation and disappearance of individual bonds are unambiguously defined, in contrast to rates of complete transcripts. We show that over previous approaches, the recently developed transient transcriptome sequencing (TT-seq) protocol allows for estimation of metabolism rates of individual bonds with least positional bias.

The transcription factor ERG recruits CCR4-NOT to control mRNA decay and mitotic progression

Control of mRNA levels, a fundamental aspect in the regulation of gene expression, is achieved through a balance between mRNA synthesis and decay. E26-related gene (Erg) proteins are canonical transcription factors whose previously described functions are confined to the control of mRNA synthesis. Here, we report that ERG also regulates gene expression by affecting mRNA stability and identify the molecular mechanisms underlying this function in human cells. ERG is recruited to mRNAs via interaction with the RNA-binding protein RBPMS, and it promotes mRNA decay by binding CNOT2, a component of the CCR4-NOT deadenylation complex. Transcriptome-wide mRNA stability analysis revealed that ERG controls the degradation of a subset of mRNAs highly connected to Aurora signaling, whose decay during S phase is necessary for mitotic progression. Our data indicate that control of gene expression by mammalian transcription factors may follow a more complex scheme than previously anticipated, integrating mRNA synthesis and degradation.

Determinants of RNA metabolism in the Schizosaccharomyces pombe genome

To decrypt the regulatory code of the genome, sequence elements must be defined that determine the kinetics of RNA metabolism and thus gene expression. Here, we attempt such decryption in an eukaryotic model organism, the fission yeast S. pombe. We first derive an improved genome annotation that redefines borders of 36% of expressed mRNAs and adds 487 non‐coding RNAs (ncRNAs). We then combine RNA labeling in vivo with mathematical modeling to obtain rates of RNA synthesis and degradation for 5,484 expressed RNAs and splicing rates for 4,958 introns. We identify functional sequence elements in DNA and RNA that control RNA metabolic rates and quantify the contributions of individual nucleotides to RNA synthesis, splicing, and degradation. Our approach reveals distinct kinetics of mRNA and ncRNA metabolism, separates antisense regulation by transcription interference from RNA interference, and provides a general tool for studying the regulatory code of genomes.

Temporal endogenous gene expression profiles in response to lipid-mediated transfection

BACKGROUND

Design of efficient nonviral gene delivery systems is limited as a result of the rudimentary understanding of the specific molecules and processes that facilitate DNA transfer.

METHODS

Lipoplexes formed with Lipofectamine 2000 (LF2000) and plasmid-encoding green fluorescent protein (GFP) were delivered to the HEK 293T cell line. After treating cells with lipoplexes, HG-U133 Affymetrix microarrays were used to identify endogenous genes differentially expressed between treated and untreated cells (2 h exposure) or between flow-separated transfected cells (GFP+) and treated, untransfected cells (GFP-) at 8, 16 and 24 h after lipoplex treatment. Cell priming studies were conducted using pharmacologic agents to alter endogenous levels of the identified differentially expressed genes to determine effect on transfection levels.

RESULTS

Relative to untreated cells 2 h after lipoplex treatment, only downregulated genes were identified ≥ 30-fold: ALMS1, ITGB1, FCGR3A, DOCK10 and ZDDHC13. Subsequently, relative to GFP- cells, the GFP+ cell population showed at least a five-fold upregulation of RAP1A and PACSIN3 (8 h) or HSPA6 and RAP1A (16 and 24 h). Pharmacologic studies altering endogenous levels for ALMS1, FCGR3A, and DOCK10 (involved in filopodia protrusions), ITGB1 (integrin signaling), ZDDHC13 (membrane trafficking) and PACSIN3 (proteolytic shedding of membrane receptors) were able to increase or decrease transgene production.

CONCLUSIONS

RAP1A, PACSIN3 and HSPA6 may help lipoplex-treated cells overcome a transcriptional shutdown due to treatment with lipoplexes and provide new targets for investigating molecular mechanisms of transfection or for enhancing transfection through cell priming or engineering of the nonviral gene delivery system.

Temporal endogenous gene expression profiles in response to polymer-mediated transfection and profile comparison to lipid-mediated transfection

BACKGROUND

Design of efficient nonviral gene delivery systems is limited by the rudimentary understanding of specific molecules that facilitate transfection.

METHODS

Polyplexes using 25-kDa polyethylenimine (PEI) and plasmid-encoding green fluorescent protein (GFP) were delivered to HEK 293T cells. After treating cells with polyplexes, microarrays were used to identify endogenous genes differentially expressed between treated and untreated cells (2 h of exposure) or between flow-separated transfected cells (GFP+) and treated, untransfected cells (GFP-) at 8, 16 and 24 h after lipoplex treatment. Cell priming studies were conducted using pharmacologic agents to alter endogenous levels of the identified differentially expressed genes to determine effect on transfection levels. Differentially expressed genes in polyplex-mediated transfection were compared with those differentially expressed in lipoplex transfection to identify DNA carrier-dependent molecular factors.

RESULTS

Differentially expressed genes were RGS1, ARHGAP24, PDZD2, SNX24, GSN and IGF2BP1 after 2 h; RAP1A and ACTA1 after 8 h; RAP1A, WDR78 and ACTA1 after 16 h; and RAP1A, SCG5, ATF3, IREB2 and ACTA1 after 24 h. Pharmacologic studies altering endogenous levels for ARHGAP24, GSN, IGF2BP1, PDZD2 and RGS1 were able to increase or decrease transgene production. Comparing differentially expressed genes for polyplexes and lipoplexes, no common genes were identified at the 2-h time point, whereas, after the 8-h time point, RAP1A, ATF3 and HSPA6 were similarly expressed. SCG5 and PGAP1 were only upregulated in polyplex-transfected cells.

CONCLUSIONS

The identified genes and pharmacologic agents provide targets for improving transfection systems, although polyplex or lipoplex dependencies must be considered.

Regulation of mRNA Levels by Decay-Promoting Introns that Recruit the Exosome Specificity Factor Mmi1

In eukaryotic cells, inefficient splicing is surprisingly common and leads to the degradation of transcripts with retained introns. How pre-mRNAs are committed to nuclear decay is unknown. Here, we uncover a mechanism by which specific intron-containing transcripts are targeted for nuclear degradation in fission yeast. Sequence elements within these "decay-promoting" introns co-transcriptionally recruit the exosome specificity factor Mmi1, which induces degradation of the unspliced precursor and leads to a reduction in the levels of the spliced mRNA. This mechanism negatively regulates levels of the RNA helicase DDX5/Dbp2 to promote cell survival in response to stress. In contrast, fast removal of decay-promoting introns by co-transcriptional splicing precludes Mmi1 recruitment and relieves negative expression regulation. We propose that decay-promoting introns facilitate the regulation of gene expression. Based on the identification of multiple additional Mmi1 targets, including mRNAs, long non-coding RNAs, and sn/snoRNAs, we suggest a general role in RNA regulation for Mmi1 through transcript degradation.