Search

Alternative polyadenylation generates numerous 3′ mRNA isoforms that can differ in their stability, structure, and function. These isoforms can be used to map mRNA stabilizing and destabilizing elements within 3′ untranslated regions (3′UTRs). Here, we examine how environmental conditions affect 3′ mRNA isoform turnover and structure in yeast cells on a transcriptome scale. Isoform stability broadly increases when cells grow more slowly, with relative half-lives of most isoforms being well correlated across multiple conditions. Surprisingly, dimethyl sulfate probing reveals that individual 3′ isoforms have similar structures across different conditions, in contrast to the extensive structural differences that can exist between closely related isoforms in an individual condition. Unexpectedly, most mRNA stabilizing and destabilizing elements function only in a single growth condition. The genes associated with some classes of condition-specific stability elements are enriched for different functional categories, suggesting that regulated mRNA stability might contribute to adaptation to different growth environments. Condition-specific stability elements do not result in corresponding condition-specific changes in steady-state mRNA isoform levels. This observation is consistent with a compensatory mechanism between polyadenylation and stability, and it suggests that condition-specific mRNA stability elements might largely reflect condition-specific regulation of mRNA 3′ end formation.

Alternative polyadenylation yields many mRNA isoforms whose 3’ termini occur disproportionately in clusters within 3’ untranslated regions. Previously, we showed that profiles of poly(A) site usage are regulated by the rate of transcriptional elongation by RNA polymerase (Pol) II (Geisberg et al., 2020). Pol II derivatives with slow elongation rates confer an upstream-shifted poly(A) profile, whereas fast Pol II strains confer a downstream-shifted poly(A) profile. Within yeast isoform clusters, these shifts occur steadily from one isoform to the next across nucleotide distances. In contrast, the shift between clusters – from the last isoform of one cluster to the first isoform of the next – is much less pronounced, even over large distances. GC content in a region 13–30 nt downstream from isoform clusters correlates with their sensitivity to Pol II elongation rate. In human cells, the upstream shift caused by a slow Pol II mutant also occurs continuously at single nucleotide resolution within clusters but not between them. Pol II occupancy increases just downstream of poly(A) sites, suggesting a linkage between reduced elongation rate and cluster formation. These observations suggest that (1) Pol II elongation speed affects the nucleotide-level dwell time allowing polyadenylation to occur, (2) poly(A) site clusters are linked to the local elongation rate, and hence do not arise simply by intrinsically imprecise cleavage and polyadenylation of the RNA substrate, (3) DNA sequence elements can affect Pol II elongation and poly(A) profiles, and (4) the cleavage/polyadenylation and Pol II elongation complexes are spatially, and perhaps physically, coupled so that polyadenylation occurs rapidly upon emergence of the nascent RNA from the Pol II elongation complex.

Alternative polyadenylation yields many mRNA isoforms whose 3’ termini occur disproportionately in clusters within 3’ UTRs. Previously, we showed that profiles of poly(A) site usage are regulated by the rate of transcriptional elongation by RNA polymerase (Pol) II (Geisberg et., 2020). Pol II derivatives with slow elongation rates confer an upstream-shifted poly(A) profile, whereas fast Pol II strains confer a downstream-shifted poly(A) profile. In yeast, upstream and downstream shifts within isoform clusters occur steadily at the nucleotide level. In contrast, changes from one isoform to the next are much smaller between clusters, even when the distances between them are relatively large. GC content in a region 13-30 nt downstream from isoform clusters is linked to Pol II elongation rate. In human cells, the upstream shift caused by a slow Pol II mutant also occurs continuously at the nucleotide level within clusters, but not between them. Pol II occupancy increases just downstream of the most speed-sensitive poly(A) sites, suggesting a linkage between reduced elongation rate and cluster formation. These observations suggest that 1) Pol II elongation speed affects the nucleotide-level dwell time allowing polyadenylation to occur, 2) poly(A) site clusters are linked to the local elongation rate and hence do not arise simply by intrinsically imprecise cleavage and polyadenylation of the RNA substrate, 3) DNA sequence elements can affect Pol II elongation and poly(A) profiles, and 4) the cleavage/polyadenylation and Pol II elongation complexes are spatially, and perhaps physically, coupled so that polyadenylation occurs rapidly upon emergence of the nascent RNA from the Pol II elongation complex.

The 3' ends of eukaryotic mRNAs are generated by cleavage of nascent transcripts followed by polyadenylation, which occurs at numerous sites within 3' untranslated regions (3' UTRs) but rarely within coding regions. An individual gene can yield many 3'-mRNA isoforms with distinct half-lives. We dissect the relative contributions of protein-coding sequences (open reading frames [ORFs]) and 3' UTRs to polyadenylation profiles in yeast. ORF-deleted derivatives often display strongly decreased mRNA levels, indicating that ORFs contribute to overall mRNA stability. Poly(A) profiles, and hence relative isoform half-lives, of most (9 of 10) ORF-deleted derivatives are very similar to their wild-type counterparts. Similarly, in-frame insertion of a large protein-coding fragment between the ORF and 3' UTR has minimal effect on the poly(A) profile in all 15 cases tested. Last, reciprocal ORF/3'-UTR chimeric genes indicate that the poly(A) profile is determined by the 3' UTR. Thus, 3' UTRs are self-contained modular entities sufficient to determine poly(A) profiles and relative 3'-isoform half-lives. In the one atypical instance, ORF deletion causes an upstream shift of poly(A) sites, likely because juxtaposition of an unusually high AT-rich stretch directs polyadenylation closely downstream. This suggests that long AT-rich stretches, which are not encountered until after coding regions, are important for restricting polyadenylation to 3' UTRs.

Significance Cells coordinate the rates of major biological processes in response to environmental conditions such as nutrient availability and temperature. Coordination mechanisms in which effects on one process are compensated for by effects on another process maintain an appropriate cellular balance. Here, we uncover a compensatory link between cleavage/polyadenylation, which generates poly(A) tails at mRNA 3′ ends, and mRNA decay. This compensatory link suggests that during cleavage/polyadenylation in the nucleus, mRNA isoforms are marked in a manner that persists upon translocation to the cytoplasm. This mark affects the activity of mRNA degradation machinery, thus influencing mRNA stability and allowing cells to maintain relatively similar levels of mRNA isoforms even when cleavage/polyadenylation rates change in response to environmental conditions., Cells have compensatory mechanisms to coordinate the rates of major biological processes, thereby permitting growth in a wide variety of conditions. Here, we uncover a compensatory link between cleavage/polyadenylation in the nucleus and messenger RNA (mRNA) turnover in the cytoplasm. On a global basis, same-gene 3′ mRNA isoforms with twofold or greater differences in half-lives have steady-state mRNA levels that differ by significantly less than a factor of 2. In addition, increased efficiency of cleavage/polyadenylation at a specific site is associated with reduced stability of the corresponding 3′ mRNA isoform. This inverse relationship between cleavage/polyadenylation and mRNA isoform half-life reduces the variability in the steady-state levels of mRNA isoforms, and it occurs in all four growth conditions tested. These observations suggest that during cleavage/polyadenylation in the nucleus, mRNA isoforms are marked in a manner that persists upon translocation to the cytoplasm and affects the activity of mRNA degradation machinery, thus influencing mRNA stability.

Yeast cells undergoing the diauxic response show a striking upstream shift in poly(A) site utilization, with increased use of ORF-proximal poly(A) sites resulting in shorter 3’ mRNA isoforms for most genes. This altered poly(A) pattern is extremely similar to that observed in cells containing Pol II derivatives with slow elongation rates. Conversely, cells containing derivatives with fast elongation rates show a subtle downstream shift in poly(A) sites. Polyadenylation patterns of many genes are sensitive to both fast and slow elongation rates, and a global shift of poly(A) utilization is strongly linked to increased purine content of sequences flanking poly(A) sites. Pol II processivity is impaired in diauxic cells, but strains with reduced processivity and normal Pol II elongation rates have normal polyadenylation profiles. Thus, Pol II elongation speed is important for poly(A) site selection and for regulating poly(A) patterns in response to environmental conditions.

Here we describe CLIP-READS, a technique that combines elements of CLIP (Cross-Linking and ImmunoPrecipitation) and 3’ READS (Jin et al., 2015) to provide a genome-wide map of mRNA 3’ isoform binding by a given messenger ribonucleoprotein (mRNP). In CLIP-READS, cells are grown to logarithmic phase and are irradiated with UV light (254nm) to form RNA-protein adducts. The protein-mRNA complexes are immunoprecipitated from cell extracts with an antibody specific to the protein of interest, after which the protein component is digested away with Pronase. Messenger RNAs are then subjected to 3’ region extraction and deep sequencing (3’ READS) (Jin et al., 2015). An input sample processed by 3’ READS in parallel allows for the relative quantification of isoform-specific binding by the mRNP of interest.

The DMS region extraction and deep sequencing (DREADS) procedure was designed to probe RNA structure in vivo and to link this structural information to specific 3' isoforms. Growing cells are treated with the alkylating agent dimethyl sulfate (DMS), which enters easily into cells and modifies RNA molecules at solvent-exposed A and C residues. RNA is isolated, and sequencing libraries are constructed in a manner that preserves the identities of individual mRNA isoforms arising from alternative cleavage/polyadenylation sites. During the cDNA synthesis step of library construction, the progress of reverse transcriptase (RT) is blocked when it encounters a DMS modification on the RNA, leading to disproportionate cDNA termination adjacent to DMS-modified positions. After paired-end deep sequencing, the downstream end of each sequenced fragment is mapped to a specific cleavage/poly(A) site representing an individual mRNA 3' isoform. The upstream mapped end of the sequenced fragment defines where the RT reaction stopped. Over the population of all sequenced fragments derived from a particular isoform, A and C positions that are overrepresented next to the upstream endpoints in the DMS sample (relative to a parallel untreated control) are inferred to have been DMS modified, and hence solvent exposed. This method thus allows in vivo structural information obtained using DMS to be linked to individual mRNA 3' isoforms. © 2019 by John Wiley & Sons, Inc.

Summary Alternative polyadenylation generates numerous 3′ mRNA isoforms that can vary in biological properties, such as stability and localization. We developed methods to obtain transcriptome-scale structural information and protein binding on individual 3′ mRNA isoforms in vivo. Strikingly, near-identical mRNA isoforms can possess dramatically different structures throughout the 3′ UTR. Analyses of identical mRNAs in different species or refolded in vitro indicate that structural differences in vivo are often due to trans-acting factors. The level of Pab1 binding to poly(A)-containing isoforms is surprisingly variable, and differences in Pab1 binding correlate with the extent of structural variation for closely spaced isoforms. A pattern encompassing single-strandedness near the 3′ terminus, double-strandedness of the poly(A) tail, and low Pab1 binding is associated with mRNA stability. Thus, individual 3′ mRNA isoforms can be remarkably different physical entities in vivo. Sequences responsible for isoform-specific structures, differential Pab1 binding, and mRNA stability are evolutionarily conserved, indicating biological function.

SummaryWe measured half-lives of 21,248 mRNA 3′ isoforms in yeast by rapidly depleting RNA polymerase II from the nucleus and performing direct RNA sequencing throughout the decay process. Interestingly, half-lives of mRNA isoforms from the same gene, including nearly identical isoforms, often vary widely. Based on clusters of isoforms with different half-lives, we identify hundreds of sequences conferring stabilization or destabilization upon mRNAs terminating downstream. One class of stabilizing element is a polyU sequence that can interact with poly(A) tails, inhibit the association of poly(A)-binding protein, and confer increased stability upon introduction into ectopic transcripts. More generally, destabilizing and stabilizing elements are linked to the propensity of the poly(A) tail to engage in double-stranded structures. Isoforms engineered to fold into 3′ stem-loop structures not involving the poly(A) tail exhibit even longer half-lives. We suggest that double-stranded structures at 3′ ends are a major determinant of mRNA stability.

In eukaryotes, RNA polymerase II-driven transcription and processing results in the formation of numerous mRNA 3' isoforms that for any given gene may differ from one another by as little as a single nucleotide. These 3' isoforms can vary in physical properties that may affect their function and stability. Here, we outline a systematic framework to measure individual mRNA 3' isoform half-lives on a genome-wide level in S. cerevisiae. Our approach utilizes the Anchor-Away system to sequester RNA polymerase II (Pol II) in the cytoplasm followed by direct single-molecule RNA sequencing to generate a highly detailed view of 3' isoform stability under most physiological conditions without many of the adverse effects associated with commonly used alternative approaches.

A cleaning protocol for the CryoXtract Automated Frozen Sample Aliquotter (AFSA) was tested for its effectiveness in eliminating sample-to-sample carryover during hands-free frozen sample aliquotting. The instrument's test platform was used to extract a statistically significant number of frozen cores from frozen bovine serum samples spiked with human DNA and from identical DNA-free frozen controls, processed in alternating succession, with a proprietary cleaning protocol performed in between samples. Measurements of DNA content by a highly sensitive real-time quantitative polymerase chain reaction assay revealed no traces of DNA in the frozen control cores, whereas they revealed the expected DNA amounts in all DNA-containing frozen cores. Analysis of high-copy DNA demonstrated that CryoXtract's AFSA's cleaning protocol effectively reduced DNA carryover by more than a million-fold.

In Saccharomyces cerevisiae, previous measurements of mRNA stabilities have been determined on a per-gene basis. We and others have recently shown that yeast genes give rise to a highly heterogeneous population of mRNAs thanks to extensive alternative 3' end formation. Typical genes can have fifty or more distinct mRNA isoforms with 3' endpoints differing by as little as one and as many as hundreds of nucleotides. In our recent paper [Geisberg et al. Cell (2014) 156: 812-824] we measured half-lives of individual mRNA isoforms in Saccharomyces cerevisiae by using the anchor away method for the rapid removal of Rpb1, the largest subunit of RNA Polymerase II, from the nucleus, followed by direct RNA sequencing of the cellular mRNA population over time. Combining these two methods allowed us to determine half-lives for more than 20,000 individual mRNA isoforms originating from nearly 5000 yeast genes. We discovered that different 3' mRNA isoforms arising from the same gene can have widely different stabilities, and that such half-life variability across mRNA isoforms from a single gene is highly prevalent in yeast cells. Determining half-lives for many different mRNA isoforms from the same genes allowed us to identify hundreds of RNA sequence elements involved in the stabilization and destabilization of individual isoforms. In many cases, the poly(A) tail is likely to participate in the formation of stability-enhancing secondary structures at mRNA 3' ends. Our results point to an important role for mRNA structure at 3' termini in governing transcript stability, likely by reducing the interaction of the mRNA with the degradation apparatus.

In eukaryotes, RNA polymerase II-driven transcription and processing results in the formation of numerous mRNA 3' isoforms that for any given gene may differ from one another by as little as a single nucleotide. These 3' isoforms can vary in physical properties that may affect their function and stability. Here, we outline a systematic framework to measure individual mRNA 3' isoform half-lives on a genome-wide level in S. cerevisiae. Our approach utilizes the Anchor-Away system to sequester RNA polymerase II (Pol II) in the cytoplasm followed by direct single-molecule RNA sequencing to generate a highly detailed view of 3' isoform stability under most physiological conditions without many of the adverse effects associated with commonly used alternative approaches.

Most eukaryotic genes are transcribed into mRNAs with alternative poly(A) sites. Emerging evidence suggests that mRNA isoforms with alternative poly(A) sites can perform critical regulatory functions in numerous biological processes. In recent years, a number of strategies utilizing high-throughput sequencing technologies have been developed to aid in the identification of genome-wide poly(A) sites. This unit describes a modified protocol for a recently published 3'READS (3' region extraction and deep sequencing) method that accurately identifies genome-wide poly(A) sites and that can be used to quantify the relative abundance of the resulting 3' mRNA isoforms. This approach minimizes nonspecific sequence reads due to internal priming and typically yields a high percentage of sequence reads that are ideally suited for accurate poly(A) identification.

In yeast, TFIID strongly associates with nearly all ribosomal protein (RP) promoters, but a TAF-independent form of TBP preferentially associates with other active promoters. RP promoters are regulated in response to growth stimuli, in most cases by a Rap1-containing activator. This Rap1-dependent activator is necessary and sufficient for TFIID recruitment, whereas other activators do not efficiently recruit TFIID. TAFs are recruited to RP promoters even when TBP and other general transcription factors are not associated, suggesting that TFIID recruitment involves a direct activator-TAF interaction. Most RP promoters lack canonical TATA elements, and they are preferentially activated by the Rap1-containing activator. These results demonstrate activator-specific recruitment of TFIID in vivo, and they suggest that TFIID recruitment is important for coordinate expression of RP genes.

Traditionally, methods for introducing specific new mutations at target loci in the yeast genome have involved the preparation of disruption or gene-replacement cassettes via multiple cloning steps. Sequences used for targeting these cassettes or integrating vectors are typically several hundred base pairs long. A variety of newer methods rely on the design of custom PCR oligonucleotides containing shorter sequence tails (∼50 nt) for targeting the locus of interest. These techniques obviate the need for cloning steps and allow construction of mutagenesis cassettes by PCR amplification. Such cassettes may be used for gene deletion, epitope tagging, or site-specific mutagenesis. The strategies differ in several ways, most notably with respect to whether they allow reuse of the selection marker and whether extra sequences are left behind near the target locus. This unit presents a summary of methods for targeted mutagenesis of Saccharomyces cerevisiae loci without cloning, including PCR-based allele replacement, delitto perfetto, and MIRAGE. Next, a protocol is provided for the delitto perfetto PCR- and oligonucleotide-based mutagenesis method, which offers particular advantages for generating several different mutant alleles of the same gene.

Nuclear proteins play key roles in the regulation of many important cellular processes. In Saccharomyces cerevisiae, many genes encoding nuclear proteins are essential. This unit describes a method termed Anchor Away that can be used to conditionally and rapidly deplete nuclear proteins of interest. It involves conditional export of the protein of interest out of the nucleus and its subsequent sequestration in the cytoplasm. This method can be used to simultaneously deplete multiple proteins from the nucleus.

Chromatin immunoprecipitation (ChIP) is a powerful and widely applied technique for detecting the association of individual proteins with specific genomic regions in vivo. Live cells are treated with formaldehyde to generate protein-protein and protein-DNA cross-links between molecules that are in close proximity on the chromatin template in vivo. DNA sequences that cross-link with a given protein are selectively enriched, and reversal of the formaldehyde cross-linking permits recovery and quantitative analysis of the immunoprecipitated DNA. As formaldehyde inactivates cellular enzymes essentially immediately upon addition to cells, ChIP provides snapshots of protein-protein and protein-DNA interactions at a particular time point, and hence is useful for kinetic analysis of events occurring on chromosomal sequences in vivo. In addition, ChIP can be combined with microarray technology to identify the location of specific proteins on a genome-wide basis. in this unit describes the ChIP procedure for Saccharomyces cerevisiae; describes the corresponding steps for mammalian cells.

Chromatin immunoprecipitation (ChIP) is a powerful and widely applied technique for detecting the association of individual proteins with specific genomic regions in vivo. Live cells are treated with formaldehyde to generate protein-protein and protein- DNA cross-links between molecules in close proximity on the chromatin template in vivo. DNA sequences that cross-link with a given protein are selectively enriched and reversal of the formaldehyde cross-link permits recovery and quantitative analysis of the immunoprecipitated DNA. As formaldehyde inactivates cellular enzymes essentially immediately upon addition to cells, ChIP provides snapshots of protein-protein and protein- DNA interactions at a particular time point, and hence is useful for kinetic analysis of events occurring on chromosomal sequences in vivo. In addition, ChIP can be combined with microarray technology to identify the location of specific proteins on a genome-wide basis. This unit describes the ChIP protocol for Saccharomyces cerevisiae; however, it is also applicable to other organisms.

Sequential Chromatin Immunoprecipitation (SeqChIP) is a powerful technique for analyzing the simultaneous association of two different proteins with genomic DNA sequences in vivo. Cellular Protein-DNA complexes are cross-linked with formaldehyde (unit Unavailable ), and are purified via two successive immunoprecipitations, with each immunoprecipitation targeting a different protein. Protein-DNA cross-links are then reversed and DNA sequences of interest are analyzed by quantitative PCR. At each genomic region, calculated SeqChIP co-occupancy values are compared to occupancy values of singly immunoprecipitated samples. The extent of enrichment brought about by the second immunoprecipitation relative to the singly immunoprecipitated sample is directly correlated with the degree of co-occupancy between the two proteins at the genomic location assayed. In principle, the technique is not limited to Saccharomyces cerevisiae. Cells from a wide variety of organisms can be used. Keywords: chromatin; immunoprecipitation; protein-DNA interactions; in vivo crosslinking

Sequential chromatin immunoprecipitation (SeqChIP) is a procedure in which formaldehyde-crosslinked, protein-DNA complexes from living cells are subjected to two sequential immunoprecipitations with antibodies of different specificity. SeqChIP has been used to address, in a qualitative manner, whether two proteins can simultaneously co-occupy a stretch of DNA in vivo. Here, we expand on our earlier work and describe theoretical and practical considerations for performing and interpreting SeqChIP experiments in a quantitative manner. We provide a detailed experimental procedure for designing and performing SeqChIP experiments as well as experimental examples of the three possible outcomes: full co-occupancy, no co-occupancy and partial co-occupancy. In some cases of partial co-occupancy, the order of immunoprecipitations in SeqChIP can strongly influence the outcome. We experimentally confirm a quantitative parameter that provides a measure of co-occupancy of two proteins on a given region of DNA and provide information on how to interpret the results of SeqChIP experiments. Our quantitative treatment of SeqChIP data substantially expands the usefulness of the technique for elucidating molecular mechanisms in vivo.

NC2 (Dr1-Drap1 or Bur6-Ydr1) has been characterized in vitro as a general negative regulator of RNA polymerase II (Pol II) transcription that interacts with TATA-binding protein (TBP) and inhibits its function. Here, we show that NC2 associates with promoters in vivo in a manner that correlates with transcriptional activity and with occupancy by basal transcription factors. NC2 rapidly associates with promoters in response to transcriptional activation, and it remains associated under conditions in which transcription is blocked after assembly of the Pol II preinitiation complex. NC2 positively and negatively affects approximately 17% of Saccharomyces cerevisiae genes in a pattern that resembles the response to general environmental stress. Relative to TBP, NC2 occupancy is high at promoters where NC2 is positively required for normal levels of transcription. Thus, NC2 is associated with the Pol II preinitiation complex, and it can play a direct and positive role at certain promoters in vivo.

Transcriptional activation by the adenovirus E1A 289R protein requires direct contacts with the TATA box-binding protein (TBP) and also displays a critical requirement for TBP-associated factors (TAFs) (T.G. Boyer and A. J. Berk, Genes Dev. 7:1810-1823, 1993; J. V. Geisberg, W. S. Lee, A. J. Berk, and R. P. Ricciardi, Proc. Natl. Acad. Sci. USA 91:2488-2492, 1994; W. S. Lee, C. C. Kao, G. O. Bryant, X. Liu, and A. J. Berk, Cell 67:365-376, 1991; and Q. Zhou, P. M. Lieberman, T. G. Boyer, and A. J. Berk, Genes Dev. 6:1964-1974, 1992). In this report, we demonstrate that the activation domain of E1A (CR3) specifically binds to two TAFs, human TAFII250 (hTAFII250) and Drosophila TAFII110 (dTAFII110). These interactions can take place both in vivo and in vitro and require the carboxy-terminal region of CR3; the zinc finger region of CR3, which binds TBP, is not needed to bind these TAFs. We mapped the E1A-binding sites on hTAFII250 to an internal region that contains a number of structural motifs, including an HMG box, a bromodomain, and direct repeats. This represents the first demonstration that hTAFII250 may serve as a target of a transcriptional activator. We also mapped the E1A binding on dTAFII110 to its C-terminal region. This is of significance since, by contrast, Sp1-mediated activation requires binding to the N-terminal domain of dTAFII110. Thus, distinct surfaces of dTAFII110 can serve as target sites for different activators. Our results indicate that E1A may activate transcription, in part, through direct contacts of the CR3 subdomains with selected components of the TFIID complex.