# Quality control of transcription start site selection by nonsense-mediated-mRNA decay

Christophe Malabat, Frank Feuerbach, Laurence Ma, Cosmin Saveanu, Alain

Jacquier

► To cite this version:

Christophe Malabat, Frank Feuerbach, Laurence Ma, Cosmin Saveanu, Alain Jacquier. Quality control of transcription start site selection by nonsense-mediated-mRNA decay. eLife, 2015, 4, 10.7554/eLife.06722. pasteur-01404014

HAL Id: pasteur-01404014

https://pasteur.hal.science/pasteur-01404014

Submitted on 6 Dec 2016

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Quality control of transcription start site selection by nonsense-mediated-mRNA decay

Christophe Malabat¹†, Frank Feuerbach¹†, Laurence Ma², Cosmin Saveanu¹, Alain Jacquier¹*

¹Institut Pasteur, UMR3525, Génétique des Interactions Macromoléculaires, Centre National de la Recherche Scientifique, Paris, France; ²Plate-Forme Génomique, Institut Pasteur, Paris, France

**Abstract** Nonsense-mediated mRNA decay (NMD) is a translation-dependent RNA quality-control pathway targeting transcripts such as messenger RNAs harboring premature stop-codons or short upstream open reading frame (uORFs). Our transcription start sites (TSSs) analysis of *Saccharomyces cerevisiae* cells deficient for RNA degradation pathways revealed that about half of the pervasive transcripts are degraded by NMD, which provides a fail-safe mechanism to remove spurious transcripts that escaped degradation in the nucleus. Moreover, we found that the low specificity of RNA polymerase II TSSs selection generates, for 47% of the expressed genes, NMD-sensitive transcript isoforms carrying uORFs or starting downstream of the ATG START codon. Despite the low abundance of this last category of isoforms, their presence seems to constrain genomic sequences, as suggested by the significant bias against in-frame ATGs specifically found at the beginning of the corresponding genes and reflected by a depletion of methionines in the N-terminus of the encoded proteins.

DOI: 10.7554/eLife.06722.001

## Introduction

Recent advances in sequencing technologies led to the detection of a wealth of new RNA transcripts and revealed that eukaryotic genomes are pervasively transcribed. In human cells, roughly 75% of the genome gives rise to RNA transcripts of various length but only an estimated 2% corresponds to protein-coding messenger RNAs (mRNAs) (*Djebali et al., 2012*). Even in a compact genome such as the budding yeast *Saccharomyces cerevisiae* (*S. cerevisiae*), hundreds of such pervasive non-coding RNAs (ncRNAs) were identified in addition to the stable ncRNAs such as transfer RNAs (tRNAs), small nucleolar RNAs (sn(o)RNAs), and ribosomal RNA (rRNAs) (*Neil et al., 2009*; *Xu et al., 2009*; *van Dijk et al., 2011*; *Geisler et al., 2012*).

Yeast ncRNAs are transcribed from nucleosome-free regions (NFRs) present throughout the genome, most often at gene promoters and terminators, and part of these RNAs represent by-products of divergent transcription initiation (*Neil et al., 2009*; *Xu et al., 2009*). Different mechanisms of RNA quality control prevent their accumulation in wild-type cells. In the nucleus, transcription from bidirectional promoters is restricted by early termination-coupled RNA degradation pathways (*Arigo et al., 2006*; *Thiebaut et al., 2006*; *Almada et al., 2013*; *Ntini et al., 2013*). In the yeast *S. cerevisiae*, degradation-coupled transcription termination of pervasive transcripts relies on the Nrd1-Nab3-Sen1 (NNS) complex, which interacts with the carboxy-terminal domain (CTD) of the RNA polymerase II (RNAPII) phosphorylated on serine 5 and triggers termination upon recognition of short sequences on the nascent RNA (*Arigo et al., 2006*; *Thiebaut et al., 2006*; *Gudipati et al., 2008*; *Vasiljeva et al., 2008*; *Schulz et al., 2013*). Early termination of these transcripts (named CUTs for

<footnote>*For correspondence: alain.jacquier@pasteur.fr</footnote>

<footnote>†These authors contributed equally to this work</footnote>

<footnote>**Competing interests:** The authors declare that no competing interests exist.</footnote>

<footnote>**Funding:** See page 21</footnote>

<footnote>**Received:** 27 January 2015</footnote>

<footnote>**Accepted:** 22 April 2015</footnote>

<footnote>**Published:** 23 April 2015</footnote>

<footnote>**Reviewing editor:** Torben Heick Jensen, Aarhus University, Denmark</footnote>

<footnote>cc Copyright Malabat et al. This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.</footnote>
**eLife digest** Eukaryotes such as animals, plants and fungi store their DNA within the nucleus of each of their cells. Genes within this DNA contain the instructions needed to make molecules of RNA; some of which can leave the nucleus and be decoded to build proteins. However, not all of the DNA that is copied into RNA actually codes for proteins. Instead, some RNA molecules are important parts of the cell's protein-making machinery in their own right, and others help to regulate the expression of genes as RNAs or proteins.

Nevertheless, many non-coding RNAs don't have such clear roles. Often these RNAs—which are called 'pervasive transcripts'—are quickly destroyed within the nucleus, but it is likely that some molecules will escape this quality-control mechanism. If the cell's protein-making machinery decodes these RNAs, it could lead to the production of faulty or harmful proteins. Recent research suggested that another quality-control mechanism, which typically eradicates incorrectly processed protein-coding RNAs, could also destroy unneeded or harmful pervasive transcripts. But it was not clear how common it was for this process—called 'nonsense-mediated decay'—to be used for this purpose.

Now Malabat, Feuerbach et al. have engineered yeast cells that lacked either the genes required to carry out nonsense-mediated decay or the ability to destroy RNA molecules in the nucleus. Experiments with these yeast cells revealed that about half of all pervasive transcripts can be destroyed via nonsense-mediated decay; this suggests that this mechanism serves as a fail-safe to prevent the build-up of these potentially harmful molecules.

Malabat, Feuerbach et al. also revealed that the enzyme complex that copies gene sequences to make RNA molecules will often also copy some extra DNA sequence from before the start of the gene. On the other hand, it is also common for this enzyme complex to miss the start of the gene and produce an RNA molecule that lacks some of the instructions needed to build the correct protein. Further experiments showed that in yeast these two kinds of incorrectly made protein-coding RNAs could both be identified and destroyed by nonsense-mediated decay as well. The next challenge will be to see to what extent these phenomena are conserved in other eukaryotes.

DOI: 10.7554/eLife.06722.002

cryptic unstable transcripts) is coupled with the recruitment of the Trf4-Air2-Mtr4 (TRAMP) complex and the Rrp6-containing nuclear RNA-exosome for rapid degradation (Wyers et al., 2005; Davis and Ares, 2006; Tudek et al., 2014). However, a significant proportion of the pervasive transcripts can escape this early nuclear quality-control step and be exported to the cytoplasm where they are targeted for degradation by the cytoplasmic 5'-3' exonuclease Xrn1 (XUT; van Dijk et al., 2011). Some pervasive transcripts seem to be immune enough to both surveillance pathways to be detectable in wild-type cells and are called SUTs (stable unannotated transcripts; Xu et al., 2009). The distinction between these different classes of transcripts is not very stringent. Their behavior can be similar, as illustrated by the fact that in the absence of Xrn1p the average levels of SUTs and CUTs increase by 7.9-fold and 3.6-fold respectively (van Dijk et al., 2011). What makes pervasive transcripts (XUTs in the first instance) highly sensitive to Xrn1p-dependent cytoplasmic degradation remains to be determined but may be linked to the presence of small spurious open reading frames (ORF) in these transcripts that will make them substrates for the nonsense-mediated mRNA decay (NMD) as has been recently demonstrated for a new class of unannotated transcripts (Smith et al., 2014).

Originally described more than two decades ago as a quality-control pathway targeting for degradation mRNAs containing a premature termination codon (Leeds et al., 1991), NMD is a translation-coupled quality-control pathway affecting cytoplasmic transcripts with restricted coding capacities (for a review see Kervestin and Jacobson, 2012). NMD targets include mRNAs harboring 'upstream open reading frames' (uORFs) (He et al., 2003; Arribere and Gilbert, 2013) and unspliced or incorrectly spliced transcripts exported to the cytoplasm (Sayani et al., 2008; Kawashima et al., 2014) accounting for roughly 10-15% of mRNAs in yeast as well as in human cells (Kervestin and Jacobson, 2012 and references therein). This number may be an underestimate since recent studies underscored the structural heterogeneity of most mRNAs at both their 5' and 3'-ends (Ozsolak et al., 2010; Arribere and Gilbert, 2013; Pelechano et al., 2013; Waern and Snyder, 2013).

Transcript heterogeneity at the 5'-end is particularly relevant for NMD targeting since the use of different transcription start sites (TSSs) will generate, for a given gene, mRNAs with different 5'-UTRs
that may include uORFs. Although NMD targets were previously analysed genome-wide using NMD-deficient strains of S. cerevisiae, the techniques used to perform these analyses (He et al., 2003; Kawashima et al., 2014; Smith et al., 2014) did not allow monitoring such individual transcript isoforms and were thus susceptible to have missed a large number of transcripts degraded by this pathway. In order to address this question in a global and systematic manner, we used a modified 5'-RACE approach (Hashimoto et al., 2009; Arribere and Gilbert, 2013) to perform a genome-wide analysis of TSSs in wild-type S. cerevisiae cells as well as in cells deficient for nuclear and NMD RNA quality-control pathways. Minor transcript isoforms targeted for degradation by NMD were identified for almost half of yeast protein-coding genes, underscoring the low specificity of TSS selection by RNAPII. In particular, our study revealed for the majority of protein-coding genes the use of TSSs downstream the ATG start codons. This phenomenon has the potential to generate N-terminally truncated proteins if the first ATG encountered by ribosomes translating this 5' truncated transcript is in the right reading frame. Such an undesirable outcome seems, however, to be counteracted by the significant depletion of in-frame relative to out-of-frame ATGs at the beginning of protein-coding genes in yeast. This bias increases the probability of such transcripts to have very short coding regions and be efficiently degraded by NMD.

Our analysis also showed that NMD restricts the accumulation of cryptic transcripts initiating inside transcribed ORFs as a consequence of altered chromatin structure and provides a fail-safe control mechanism for the removal of pervasive transcripts that escaped degradation by the nuclear quality-control pathway, reminiscent of what has been previously shown for some unspliced pre-mRNAs (Sayani and Chanfreau, 2012). This includes not only a large fraction of the XUTs and SUTs but also a fraction of the CUTs and previously unannotated transcripts, which accumulate to substantial levels only when both pathways are inactive.

## Results

### TSS sequencing technique

To identify TSSs genome-wide and with high specificity we used a modified genomic 5'-RACE approach (Hashimoto et al., 2009; Arribere and Gilbert, 2013) that we called TSS sequencing, which involves a biotin purification step and allows the selective enrichment of the 5'-ends of capped transcripts (see Figure 1A and 'Materials and methods'). We evaluated the 'false-discovery' rate of the method for the identification of TSSs (that is, the proportion of sequencing reads not actually mapping to the 5'-end of capped RNAs) by comparing libraries made using samples treated or not with the tobacco acid pyrophosphatase (TAP), which is required for efficient and specific ligation of the biotinylated primer to the 5'-end of capped-RNAs. To extend this analysis to transcripts that are unstable in wild-type cells, we prepared libraries using RNAs extracted from *upf1Δrrp6Δ* double mutant strains. To minimize polymerase chain reaction (PCR) amplification biases and provide an internal control, the S. cerevisiae poly(A)⁺ RNAs, treated or not with TAP, were mixed with an equal amount of TAP-treated poly(A)⁺ RNAs from *Schizosaccharomyces pombe* just prior to the ligation step. After normalization using the S. pombe sequencing reads and removal of the ribosomal DNA ones, these experiments generated 2,844,877 reads when TAP was used compared with 72,435 reads when TAP was omitted. The ratio between these two numbers provides an upper limit for the false-discovery rate of TSS identification of 2.5%. Analysis of the cumulative 5'-end read counts per nucleotide around the start codon for all protein-coding genes showed that, genome-wide, 77% (2,190,211) of the reads obtained for TAP-treated samples could be mapped within a 200 nucleotide region upstream of the start codons with a maximum at around 30 nucleotides upstream from ATGs, while only 22% mapped to the same region when this treatment was omitted (Figure 1B). However, even in the latter case, the number of mapped reads also peaked at around 30 nucleotides upstream of the ATGs, suggesting that a substantial proportion of these sequences likely correspond to genuine TSSs even though they were generated in the absence of TAP treatment. These data suggest that the real false-discovery rate for TSS is substantially less than 2.5%.

### Repeatability and reproducibility of TSS sequencing

To estimate the repeatability of the experiment, two independent biological replicates were generated in parallel from wild-type cells and sequenced on different lanes of an Illumina HiSeq 2500 sequencer. The number of reads mapping to the same genomic position were highly correlated
<figure><img src="image_3.png" /></figure>

**Figure 1.** Transcription start site sequencing. (A) Schematic view of the methodology used to produce the transcription start site (TSS) sequence-tag libraries (RNA molecules are in blue, DNA molecules in red). (B) All protein-coding genes were aligned on the A of their annotated ATG start codon and the distribution of the TSSs read counts was computed for each position in a window from -200 to +200 nucleotides for samples treated (blue curve) or not (red curve) with tobacco acid pyrophosphatase (TAP). The insert within the figure shows a zoomed view of the -100 to +50 nucleotide region. (C) Correlation between two biologically independent replicates (replicate 1: library L5p_03.WT, replicate 2: library L5p_04.WT; see **Table 1**). The Pearson's correlation value (ρ) between the read counts of the 174,151 TSSs identified in the two data sets is indicated.

DOI: 10.7554/eLife.06722.003

The following figure supplements are available for figure 1:

**Figure supplement 1.** Reproducibility of transcription start site (TSS) sequencing.

DOI: 10.7554/eLife.06722.004

**Figure supplement 2.** Transcription start site (TSS) consensus sequences.

DOI: 10.7554/eLife.06722.005

(Figure 1C; Pearson's correlation coefficient (ρ) = 0.89, see also Table 1). When we compared data sets for the wild-type strain generated months apart in our laboratory (reproducibility test), we observed lower correlation coefficients (ρ = 0.71) likely reflecting biological variation inherent to samples prepared at different moments (Figure 1—figure supplement 1A). However, the correlation between data generated within our lab is still higher than the correlation between our data and results previously published by *Pelechano et al.* (2013) (ρ = 0.58; Figure 1—figure supplement 1B) or Arribere and Gilbert (ρ = 0.48; *Arribere and Gilbert, 2013*) as well as between these two sets of data (ρ = 0.55).

## TSS consensus sequences

The complementary DNA (cDNA) sequences were aligned with the genomic sequence (allowing one mismatch in the seed sequence; see 'Materials and methods'). Analysis of the TSSs and their surrounding sequences using the Web-LOGO algorithm (Crooks et al., 2004) identified a consensus
<table><caption>**Table 1.** Libraries generated and analysed in this study</caption><thead><tr><th>Library</th><th>Genotype</th><th>Total reads count</th><th>Unique reads count</th><th>Unique reads count mapped on S. cerevisiae genome</th><th>Unique reads count mapped on S. pombe genome</th></tr></thead><tbody><tr><td colspan="6">TSS-sequencing</td></tr><tr><td>L5p_01</td><td>WT</td><td>8,650,655</td><td>5,862,245</td><td>5,631,307</td><td>-</td></tr><tr><td>L5p_01</td><td>upf1Δ</td><td>9,624,626</td><td>7,161,842</td><td>6,733,121</td><td>-</td></tr><tr><td>L5p_01</td><td>rrp6Δ</td><td>10,949,693</td><td>7,640,701</td><td>7,155,101</td><td>-</td></tr><tr><td>L5p_01</td><td>upf1Δrrp6Δ</td><td>10,839,311</td><td>7,681,910</td><td>7,231,052</td><td>-</td></tr><tr><td>L5p_02</td><td>WT</td><td>10,139,398</td><td>4,207,854</td><td>3,949,042</td><td>-</td></tr><tr><td>L5p_02</td><td>upf1Δ</td><td>10,607,408</td><td>4,325,398</td><td>4,067,154</td><td>-</td></tr><tr><td>L5p_02</td><td>rrp6Δ</td><td>22,861,396</td><td>10,269,822</td><td>9,711,133</td><td>-</td></tr><tr><td>L5p_02</td><td>upf1Δrrp6Δ</td><td>13,803,020</td><td>6,526,805</td><td>6,207,357</td><td>-</td></tr><tr><td>L5p_03</td><td>WT</td><td>11,525,631</td><td>4,670,586</td><td>4,343,264</td><td>-</td></tr><tr><td>L5p_03</td><td>upf1Δ</td><td>9,832,590</td><td>3,338,679</td><td>3,122,942</td><td>-</td></tr><tr><td>L5p_03</td><td>rrp6Δ</td><td>11,811,879</td><td>5,228,409</td><td>4,905,881</td><td>-</td></tr><tr><td>L5p_03</td><td>upf1Δrrp6Δ</td><td>17,157,244</td><td>7,897,030</td><td>7,430,231</td><td>-</td></tr><tr><td>L5p_04</td><td>WT</td><td>13,282,069</td><td>10,019,742</td><td>6,032,695</td><td>-</td></tr><tr><td>L5p_04</td><td>upf1Δ</td><td>14,162,741</td><td>10,872,408</td><td>6,724,910</td><td>-</td></tr><tr><td>L5p_04</td><td>set2Δ</td><td>14,714,187</td><td>11,048,398</td><td>6,591,345</td><td>-</td></tr><tr><td>L5p_04</td><td>upf1Δset2Δ</td><td>16,958,167</td><td>12,121,894</td><td>8,081,113</td><td>-</td></tr><tr><td>L5p_05</td><td>WT</td><td>11,270,824</td><td>4,446,988</td><td>4,172,618</td><td>-</td></tr><tr><td>L5p_05</td><td>upf1Δ</td><td>12,093,631</td><td>4,599,307</td><td>4,323,962</td><td>-</td></tr><tr><td>L5p_05</td><td>set2Δ</td><td>18,047,134</td><td>7,132,585</td><td>6,724,933</td><td>-</td></tr><tr><td>L5p_05</td><td>upf1Δset2Δ</td><td>12,719,564</td><td>5,637,310</td><td>5,333,480</td><td>-</td></tr><tr><td>L5p_06</td><td>WT</td><td>12,253,253</td><td>4,800,799</td><td>4,463,020</td><td>-</td></tr><tr><td>L5p_06</td><td>upf1Δ</td><td>10,481,413</td><td>3,402,317</td><td>3,181,752</td><td>-</td></tr><tr><td>L5p_06</td><td>set2Δ</td><td>12,179,190</td><td>4,448,762</td><td>4,167,876</td><td>-</td></tr><tr><td>L5p_06</td><td>upf1Δset2Δ</td><td>14,269,228</td><td>5,584,019</td><td>5,227,982</td><td>-</td></tr><tr><td>L5p_07</td><td>upf1Δrrp6Δ + TAP</td><td>8,890,286</td><td>2,973,612</td><td>1,621,712</td><td>1,134,673</td></tr><tr><td>L5p_07</td><td>upf1Δrrp6Δ - TAP</td><td>13,836,172</td><td>3,214,470</td><td>149,092</td><td>2,700,018</td></tr><tr><td>L5p_08</td><td>upf1Δrrp6Δ + TAP</td><td>9,689,188</td><td>2,768,418</td><td>1,341,975</td><td>1,209,108</td></tr><tr><td>L5p_08</td><td>upf1Δrrp6Δ - TAP</td><td>9,765,672</td><td>2,390,794</td><td>82,685</td><td>2,074,096</td></tr><tr><td>L5p_09</td><td>upf1Δset2Δ + TAP</td><td>11,885,938</td><td>3,793,668</td><td>1,976,161</td><td>1,555,818</td></tr><tr><td>L5p_09</td><td>upf1Δset2Δ - TAP</td><td>11,105,585</td><td>2,936,749</td><td>117,202</td><td>2,552,264</td></tr><tr><td>L5p_10</td><td>upf1Δset2Δ + TAP</td><td>11,665,986</td><td>3,797,708</td><td>1,343,091</td><td>2,147,100</td></tr><tr><td>L5p_10</td><td>upf1Δset2Δ - TAP</td><td>11,000,429</td><td>2,476,135</td><td>59,450</td><td>2,145,439</td></tr><tr><td colspan="6">RNAseq</td></tr><tr><td>LT_01</td><td>WT</td><td>8,104,047</td><td>7,257,423</td><td>6,634,522</td><td>1,761,046</td></tr><tr><td>LT_01</td><td>upf1Δ</td><td>11,137,269</td><td>10,129,315</td><td>9,257,158</td><td>2,440,195</td></tr><tr><td>LT_01</td><td>xrn1Δ</td><td>11,619,211</td><td>10,631,126</td><td>9,737,924</td><td>2,310,274</td></tr><tr><td>LT_01</td><td>upf1Δxrn1Δ</td><td>7,947,627</td><td>7,299,151</td><td>6,645,353</td><td>1,678,281</td></tr><tr><td>LT_02</td><td>WT</td><td>34,611,003</td><td>27,825,118</td><td>22,658,503</td><td>3,635,624</td></tr><tr><td>LT_02</td><td>upf1Δ</td><td>29,379,233</td><td>24,421,717</td><td>19,815,913</td><td>2,966,963</td></tr><tr><td>LT_02</td><td>xrn1Δ</td><td>26,816,267</td><td>22,686,420</td><td>18,245,251</td><td>2,493,471</td></tr><tr><td>LT_02</td><td>upf1Δxrn1Δ</td><td>25,466,016</td><td>21,800,532</td><td>16,996,728</td><td>2,436,568</td></tr></tbody></table>

DOI: 10.7554/eLife.06722.023
sequence around TSSs (Figure 1—figure supplement 2A) similar to the previously reported one derived from a smaller data set (Zhang and Dietrich, 2005). In particular, we observed a very strong bias to start at a purine (88% of mapped TSSs), usually following a pyrimidine (76% of the mapped TSSs), and the enrichment for an A at position -8 relative to the TSS ($A(N)_6PyPu$ consensus). Surprisingly, 58% of TSS reads starting with a pyrimidine when aligned on the genome (12% of the mapped TSSs) show a mismatch at their first nucleotide, most of the time an A instead of the encoded pyrimidine (Figure 1—figure supplement 2B). Moreover, in 32% of these cases the surrounding genomic sequences exhibited a specific consensus, $A(N)_6PyAAA$ (where the underlined base is the mapped TSS; Figure 1—figure supplement 2C). These observations suggest that, in these cases, transcription actually initiates on the A following the pyrimidine and that an additional A is added at the 5'-end of the transcript, possibly by a back-tracking mechanism as proposed in the model described in Figure 1—figure supplement 2D. This may also apply to other TSSs mapped on a pyrimidine and showing a mismatched first nucleotide. It thus appears from this observation that, even though 12% of the TSSs mapped on a pyrimidine, transcription initiation actually occurred on a pyrimidine in less than 5% of cases (42% of non-mismatched cDNAs out of 12% of cDNAs aligned on a pyrimidine).

## The nuclear and cytoplasmic RNA degradation pathways cooperate in the removal of pervasive transcripts

To assess the relative contributions of the cytoplasmic NMD and nuclear RNA control pathways in shaping the yeast transcriptome, the genes encoding Upf1, an RNA helicase essential for NMD, or Rrp6, a nuclear exosome catalytic subunit, were deleted (He et al., 2003; Wyers et al., 2005). Since, for a given gene, several closely spaced TSSs can be identified (Pelechano et al., 2013), we used a peak-calling procedure to define, in three biological replicates, TSS clusters corresponding to transcript isoforms with closely spaced 5'-ends (TSSCs; see 'Materials and methods'). This analysis allowed us to identify 17,812 TSSCs, among which 5927 could be assigned to 5'-ends of mRNAs corresponding to 5231 ORFs (O category in Supplementary file 1), as previously defined (Pelechano et al., 2013). Among the remaining TSSCs, 502 and 3644 were assigned to stable ncRNAs or repeated sequences (F category in Supplementary file 1) and previously described pervasive transcripts respectively (C, X and S categories for CUTs, XUTs and SUTs in Supplementary file 1; Neil et al., 2009; Xu et al., 2009; van Dijk et al., 2011).

In contrast to mRNA-TSSCs, the majority of which were unaffected in the absence of Upf1, Rrp6 or both (Figure 2A) and, as expected (Wyers et al., 2005), CUT-TSSCs were strongly stabilized in the absence of Rrp6 (Figure 2B; 79% stabilized significantly, as determined using the moderated estimation of fold change and dispersion function of DESeq2—Love et al., 2014; see 'Materials and methods'). In contrast, only 23% of the CUT-TSSCs increased significantly in the absence of Upf1. The opposite situation was observed for XUTs and SUTs, which were more sensitive to the cytoplasmic NMD quality-control pathway (28% vs 52% significantly stabilized in the absence of Rrp6 or Upf1, respectively; Figure 2C,D and Supplementary file 1). Collectively, 39% and 54% of the pervasive transcripts (CUTs, SUTs and XUTs) were significantly sensitive to the absence of Upf1 in a wild-type or *rrp6Δ* background, respectively, indicating that NMD targets about half of the pervasive transcripts. The absence of both Rrp6 and Upf1 had an additive effect on the stabilization of pervasive transcripts. In the double mutant, 76% and 92% of the XUTs/SUTs and CUTs, respectively, were significantly stabilized, suggesting that NMD provides a fail-safe control mechanism for pervasive transcripts that escaped degradation in the nucleus, reminiscent of what has been previously shown for some unspliced pre-mRNAs (Sayani and Chanfreau, 2012).

## XUTs and SUTs carry short spurious ORFs that target them for degradation by NMD

The accumulation of XUTs in the absence of Upf1 suggests that these transcripts, originally described as Xrn1-sensitive (van Dijk et al., 2011), may be primarily targeted for degradation by NMD. We tested this hypothesis by performing northern blot hybridization for three XUTs (Figure 3A) and genome-wide transcriptome analyses of wild-type, *upf1Δ*, *xrn1Δ*, and *upf1Δxrn1Δ* cells (Figure 3B–D). Since the absence of Xrn1 leads to the accumulation of decapped RNAs (Hsu and Stevens, 1993), we could not use the TSS sequencing methodology applied for the other mutants and we thus used, for this particular experiment, a classical RNA sequencing approach for transcript quantification
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Differential effect of *UPF1* and/or *RRP6* deletion on mRNAs and pervasive transcripts. Frequency distribution of the ratios of transcription start site clusters (TSSCs) read counts in *upf1Δ* (red), *rrp6Δ* (blue), or *upf1Δrrp6Δ* (purple) compared with wild-type (WT) for mRNAs (A), cryptic unstable transcripts (CUTs) (B), Xrn1-sensitive transcripts (XUTs) (C) and stable unannotated transcripts (SUTs) (D). The dashed vertical lines mark a twofold increase in TSSC read counts in the various mutants relative to the wild-type. The number of identified TSSCs and of features to which they were assigned is indicated for each transcript class. Note that CUTs, XUTs and SUTs constituting overlapping transcript populations, when a TSSC was assigned to a pervasive transcript annotated in more than one of these classes, we arbitrarily associated the corresponding TSSC in priority to CUTs, then to XUTs and finally to SUTs.</figcaption></figure>

DOI: 10.7554/eLife.06722.006

(see 'Materials and methods'). Moreover, since the absence of Xrn1 affects the overall cellular mRNA content (Sun et al., 2013), an aliquot of a S. pombe culture was added to the cell pellet before RNA extraction to provide an independent internal control for normalization of the results. After having verified that the overall quantifications with 'RNAseq' gave results similar to those obtained by TSS sequencing for the wild-type and *upf1Δ* strains (Figure 3B and Supplementary file 2), we used the former technique to analyse the genome-wide effect of the *XRN1* deletion. Consistent with previous reports (van Dijk et al., 2011; Sun et al., 2013), deletion of *XRN1* resulted in a global increase of mRNAs, XUTs and SUTs compared with the wild-type (Figure 3C,D). In contrast to mRNAs, the extent to which XUTs and SUTs were stabilized in *xrn1Δ* was very similar to the one observed in *upf1Δ* cells (Figure 3C,D). Furthermore, deleting *XRN1* in an *upf1Δ* background had almost no additional stabilizing effect on XUTs and SUTs, in contrast to mRNAs. This epistatic relationship observed between *xrn1Δ* and *upf1Δ* for the stabilization of XUTs and SUTs is consistent with these transcripts being primarily targeted by NMD, and with Xrn1 acting as a downstream effector of this pathway.

In contradiction to the fact that XUTs and SUTs have been designated as 'non-coding', the strong effect NMD inactivation had on these transcripts indicates that they must, at some point, be translated. Indeed, all XUTs and SUTs carry spurious ORFs (Figure 3—figure supplement 1A) and some of them have been identified in ribosomal profiling experiments and shown to encode short peptides (Ingolia et al., 2009; Smith et al., 2014). However, the high sensitivity to NMD of XUTs and SUTs compared with mRNAs suggests the presence of specific features. Assuming that the first
<figure><img src="image_5.png" /><figcaption>**Figure 3.** XUTs are primarily targeted for degradation by the nonsense-mediated mRNA decay (NMD). (A) Northern blot analysis of three different Xrn1-sensitive transcripts (XUTs) in wild-type (WT), *upf1Δ*, *xrn1Δ*, and *upf1Δxrn1Δ* cells. ACT1 was used as a loading control. (B) Frequency distribution of the ratios of transcription start site clusters (TSSCs) read counts between *upf1Δ* and wild-type cells obtained with two different methods for library preparation and normalization procedures. 'TSSCs' refers to data obtained using the protocol developed to identify TSSs and 'RNAseq' to the protocol used to obtain the whole transcriptome (see 'Materials and methods'). Only reads corresponding to annotated mRNAs, XUTs and SUTs were included in the analysis (see **Supplementary file 2**). (C) and (D) Frequency distribution of the ratios of read counts for *upf1Δ* (red), *xrn1Δ* (blue) or *upf1Δxrn1Δ* (purple) compared with wild-type for mRNAs and XUTs and SUTs respectively. The vertical dashed lines mark a twofold increase in read counts.</figcaption></figure>

DOI: 10.7554/eLife.06722.007

The following figure supplement is available for figure 3:

**Figure supplement 1.** The presence of short open reading frames (ORFs) and long 3'-UTRs is a hallmark of natural nonsense-mediated mRNA decay (NMD) substrates.

DOI: 10.7554/eLife.06722.008

encountered ORF is translated, XUTs and SUTs carry, on average, much shorter ORFs and longer 3'-UTRs than mRNAs (Figure 3—figure supplement 1A). The presence of both a long 3'-UTR and a short ORF increases NMD efficiency in budding yeast (Decourty et al., 2014) and thus explain why these transcripts are efficiently targeted for degradation by this pathway. Furthermore, the few XUTs and SUTs found not to be sensitive to NMD have on average significantly shorter 3'-UTRs (Figure 3—figure supplement 1B), which is in agreement with the recently published observation that, in ribosomal profiling experiments, the length of RNA downstream of the ribosome protected region is significantly longer for NMD-sensitive compared with NMD-insensitive unannotated RNA transcripts (Smith et al., 2014).

## Identification of additional pervasive transcripts

Amongst the 17,812 TSSCs identified in this study, 7739 could not be assigned to previously annotated transcripts (ORFs, stable ncRNAs or pervasive transcripts) and were only detected in mutants cells (Figure 4). These newly identified transcripts, expressed at low levels (Figure 4—figure
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Novel transcripts revealed upon deletion of UPF1 and/or RRP6. (A–C) Frequency distribution of the ratios of transcription start site clusters (TSSCs) read counts in *upf1Δ* (red), *rrp6Δ* (blue) or *upf1Δrrp6Δ* (purple) compared with wild-type for transcripts initiating within intergenic regions—'intergenic TSSCs' (A), from within an mRNA transcribed region but antisense to the mRNA—A-TSSCs (B) or within an mRNA transcribed region, but in the sense orientation with respect to the mRNA—B-TSSCs (C). (D) Schematic representation of the various classes of TSSCs described above. The small blue and orange vertical bars represent individual TSSs within TSSCs (dashed lines). (E) Northern blot analysis of poly(A)⁺ RNA from wild-type, *upf1Δ*, *rrp6Δ* and *upf1Δrrp6Δ*. The category to which the transcripts belong is indicated on the right. ACT1 was used as a loading control.</figcaption></figure>

DOI: 10.7554/eLife.06722.009

The following figure supplements are available for figure 4:

**Figure supplement 1.** Distribution of read counts for different classes of transcripts.

DOI: 10.7554/eLife.06722.010

**Figure supplement 2.** Transcription start site (TSS) consensus sequences for the different classes of transcripts.

DOI: 10.7554/eLife.06722.011

**supplement 1)**, originated from intergenic regions (3011 'intergenic-TSSCs' originating from 1843 intergenic regions) as well as from within mRNA transcribed regions either in sense (2978 'B-TSSCs' within 1889 mRNAs) or in antisense orientation (1750 'A-TSSCs' antisense to 1186 mRNAs; **Supplementary file 1**). While 'B-TSSCs' were sensitive to *upf1Δ* and almost not affected by the
deletion of *RRP6* (Figure 4C and Figure 4—figure supplement 1), 'intergenic' and 'A-TSSCs' were affected by the two mutations (Figure 4A,B and Figure 4—figure supplement 1). Furthermore, the absence of both *Rrp6* and *Upf1* had an additive effect on the accumulation of these last two classes, as is the case for previously identified pervasive transcripts. Combining these two mutations might sometimes even have a synergistic and not only an additive effect, since some of these transcripts were readily detectable only in the *upf1Δrrp6Δ* double mutant (e.g., the transcript found antisense to *MAL12/32* in Figure 4E). Analysis of individual TSSs signatures for these three classes of previously unannotated transcripts using the Web-LOGO algorithm (Crooks et al., 2004) yielded a consensus sequence almost indistinguishable from the one obtained for individual TSSs assigned to known mRNAs or previously identified pervasive transcripts (Figure 4—figure supplement 2), suggesting that the underlying DNA sequence plays an important role in TSS selection by the scanning polymerase.

## Co-transcriptional histone modifications and NMD cooperate to restrict the accumulation of internally initiated transcripts

Even though transcripts initiated inside ORFs are expected to be targeted for degradation by NMD due to the presence of short spurious ORFs, the high number of internally initiated transcripts ('A and B-TSSCs') identified in cells lacking *UPF1* was surprising. The synthesis of such transcripts is normally repressed within transcribed regions unless the proper chromatin structure cannot be re-established in the wake of RNAPII (Smolle and Workman, 2013). In particular, histone methylation by Set2 was shown to be a key determinant of this repression. We therefore analysed the impact of *SET2* deletion on 'A' and 'B' TSSCs identified in the *upf1Δ* strain. While deletion of *SET2* had little effect on the internal TSSCs identified in the single *upf1Δ* mutant, it revealed new Set2-sensitive 'A' and 'B' TSSCs (Figure 5 and Supplementary file 3). As expected and in contrast to A-TSSCs and B-TSSCs, deletion of *SET2* had no global effect on ORF TSSCs and only a marginal effect on CUTs, XUTs and SUTs or 'intergenic' TSSCs (Figure 5—figure supplement 1). In agreement with the observed enrichment in Set2-catalyzed H3K36 methylation towards the 3' end of ORFs (Pokholok et al., 2005), the 'A' and 'B' TSSCs more sensitive to Set2 were located further away from the mRNA 5'-ends (Figure 6A–C). Analysis of the sequence surrounding individual TSSs from these TSSCs identified a pattern almost identical to the one found for ORF-TSSs (Figure 6—figure supplement 1), confirming that they corresponded to bona fide transcription initiation events. However, unlike the TSSs associated with other features identified in this study, the 'B' TSSs were not associated with a strong NFR

<figure><img src="image_7.png" /><figcaption>**Figure 5.** Effect of the absence of *SET2* on transcription start sites (TSSs) identified inside open reading frames (ORFs). (A) Frequency distribution of the ratios of TSS clusters (TSSCs) read counts in *upf1Δset2Δ* compared with *upf1Δ* cells for the intragenic A-TSSCs. The comparison was performed for all the TSSCs (blue line) and for the ones identified in a single *upf1Δ* mutant (red line). (B) As in (A) but for B-TSSCs.</figcaption></figure>

DOI: 10.7554/eLife.06722.012

The following figure supplement is available for figure 5:

**Figure supplement 1.** Deletion of *SET2* specifically increases the expression level of intragenic transcription start site clusters (TSSCs).

DOI: 10.7554/eLife.06722.013
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Transcription start site (TSSs) identified inside open reading frames (ORFs) show a differential sensitivity to the absence of *SET2* according to their position along the mRNA. (A) Visualization of the TSS reads at the *YIL136w* and *YMR114c* loci in the wild-type (WT; black), *upf1Δ* (red), *set2Δ* (green), and *upf1Δset2Δ* (orange) cells. The blue arrows represent ORFs and the horizontal red bar the position of the probes used for the Northern blots displayed on the right. Arrowheads in the right panel indicate the position of the full-length and internally initiated (B1, B2, and B) transcripts. (B) and (C) Frequency distribution of A-TSSCs and B-TSSCs read counts respectively in *upf1Δset2Δ* vs *upf1Δ* according to the distance from their associated mRNA TSS. Blue and red lines are for TSSCs sensitive and insensitive to the deletion of *SET2* respectively. (D) Frequency distribution of read counts for the A- (blue) and B- (red) TSSCs in *upf1Δset2Δ* compared to *set2Δ* cells.
DOI: 10.7554/eLife.06722.014</figcaption></figure>

The following figure supplements are available for figure 6:

**Figure supplement 1.** Consensus sequences around the transcription start sites (TSSs) for (A) A-TSSCs and (B) B-TSSCs identified in the *upf1Δset2Δ* mutant generated using the Web-LOGO algorithm (Crooks et al., 2004).
DOI: 10.7554/eLife.06722.015

**Figure supplement 2.** Transcription start sites (TSSs) are associated within nucleosome-free regions.
DOI: 10.7554/eLife.06722.016
(Figure 6—figure supplement 2). Instead, analysing the nucleosome density around these TSSs revealed the presence of a weak NFR surrounded by two well-positioned nucleosomes. This particular pattern might explain the high sensitivity of 'B' TSSs to mutations affecting chromatin structure, such as the deletion of the Set2 histone methyl-transferase (Smolle and Workman, 2013).

Importantly, 55% of the A-TSSCs and 40% of the B-TSSCs, whether repressed by Set2 or not, were sensitive to Upf1 (i.e., increased significantly in *upf1Δset2Δ* compared with set2Δ cells; Figure 6D and Supplementary file 3), revealing NMD to be an important quality-control mechanism to eliminate internally initiated transcripts.

## Almost half of the coding genes produced transcript isoforms sensitive to NMD

In contrast to their strong impact on the steady state levels of a variety of pervasive transcripts, including previously unannotated ones, inactivation of nuclear and/or cytoplasmic RNA quality-control pathways had a relatively minor global effect on the steady state levels of coding transcripts (Figure 2).

Consistent with previously reported data (He et al., 2003; Wyers et al., 2005), the effect on mRNAs was larger upon deletion of *UPF1* (with 17% of the mRNAs-TSSCs showing a significant increase in the absence of Upf1) than upon deletion of *RRP6* and deletion of both genes simultaneously did not show any additive effect (Figure 2A). Yet, due to the mRNA 5'-end heterogeneity, analysing individual TSSs gave a different picture. Indeed, transcript isoforms in which TSSs where followed by uORFs were, globally, stabilized in the absence of *UPF1*, while those in which the annotated start codon was the first ATG downstream the TSS were mostly unaffected (Figure 7A, solid lines; Figure 7B). A large fraction of these NMD-sensitive transcripts corresponded to minor isoforms, explaining why NMD had only a weak effect on the overall mRNA-TSSC levels (Figure 2A). Yet, for 1129 out of the 5231 active genes (22%), a fraction of their transcript isoforms carried at least one uORF and was significantly sensitive to NMD. Isoforms carrying more than one uORF appeared even more sensitive to NMD (Figure 7A), suggesting that the first AUGs of uORFs containing transcripts are not always efficiently used for translation initiation, likely because they are not in a favorable context (Arribere and Gilbert, 2013).

Unexpectedly, another important category of transcript isoforms found to be sensitive to NMD corresponded to TSSs mapping downstream to the annotated ORF ATG (here called iTSSs; Supplementary file 4). Only observed after TAP treatment (see insert in Figure 1B), these iTSSs must correspond to genuine capped-RNAs. We considered genes as having iTSSs (3327 genes; 64% of the expressed genes) if they contained at least four iTSS reads, which accounted for 99% of all iTTS reads. A fraction of transcript isoforms initiated at iTSSs followed by an out-of-frame ATG and significantly stabilized in the absence of Upf1 could be identified for 1821 out of the 5231 expressed genes (35%, Figure 7A,B). In contrast to the weak overall effects observed for TSSs located directly upstream of the start codon, the mean expression level of iTSSs significantly increased in the *upf1Δ* strain relative to the wild-type (Figure 7C, and Figure 7—figure supplement 1). Note that RNAPII transcription initiating preferentially at purines, the A and G of the annotated start codons are quite frequently used as sites of transcription initiation and the first nucleotide of codons are preferential sites of internal initiation because these positions are enriched in purines (Figure 7—figure supplement 1B,C; Mackiewicz et al., 1999). The consensus sequence for the iTSSs was not different from the general consensus sequence for RNAPII transcription initiation (Figure 7—figure supplement 2).

Altogether, transcript isoforms significantly stabilized in NMD-deficient cells, either because they carry a uORF or because they initiated at an out-of-frame iTSS, were found in 2437 genes; that is, 47% of the 5231 expressed genes.

In transcripts initiated at iTSSs, the first encountered ATGs are most of the time out-of-frame relative to the main ORF (Figure 8A). The first ORFs of the corresponding transcripts are thus short and followed by a long 3'-UTR, making them excellent NMD substrates. This bias towards out-of-frame ATGs was specifically observed at the beginning of genes for which iTSSs were identified (Figure 8B,C). This correlated with a significantly lower frequency of methionine in the N-terminal part of the corresponding proteins (Figure 8D), even though methionines were found to be globally underrepresented at the beginning of all yeast proteins. The frequent use of iTSSs by RNAPII thus does not generally result in the production of truncated proteins and tends to generate transcripts that are sensitive to NMD.
<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 7.** Deletion of *UPF1* reveals numerous minor mRNA-associated transcription start sites (TSSs). (A) Genes were aligned by their start codon and the log₂ of ratios of TSS reads for *upf1Δ* vs wild-type was plotted for TSSs upstream or downstream (iTSSs) the annotated ATG start codons, as depicted in the right panel. The main open reading frames (ORFs) are represented by large blue arrows and upstream ORFs (uORFs) or small internal out-of-frame ORFs by small orange arrows. The thin blue arrows indicate the TSSs. (B) The curve represents the probability at each nucleotide position that the distributions of reads corresponding to the red and blue curves shown in **A** are the same. The dashed red line marks the 0.05 p-value. (C) Genes were aligned by their annotated start codons (A of the ATG at position +1) and the cumulative TSS read counts per nucleotide (smoothed over 11 nucleotides) was plotted for the wild-type (black) and *upf1Δ* (red) cells. Inset: Magnification of the +1 to +50 region. The p-value < 2.2 x 10⁻¹⁶ is the probability (ANOVA test) that the distributions of the values, per nucleotide, for the red and black curves are the same within the +1 to +50 nucleotides interval.</figcaption></figure>

DOI: 10.7554/eLife.06722.017

The following figure supplements are available for figure 7:

**Figure supplement 1.** Initiation of iTSS at purines reflects a bias in codon composition.

DOI: 10.7554/eLife.06722.018

**Figure supplement 2.** (iTSSs) consensus sequences.

DOI: 10.7554/eLife.06722.019

Discussion

Genome-wide identification of TSSs

We describe here the use of a modified 5'-RACE technique to identify the 5'-ends of capped RNA, which map RNAPII TSSs. The method, which includes a streptavidin-biotin purification step, is highly specific and can be applied to any eukaryotic organism. Our results provide the most comprehensive genome-wide identification of TSSs in budding yeast, not only for mRNAs but also for a wide range of pervasive transcripts, including previously non-annotated ones. This was made possible by combining this highly specific TSS determination technique with the use of mutants affecting both the nuclear exosome and the cytoplasmic NMD pathway.

Detailed analyses of the sequencing reads and the corresponding upstream genomic sequences extended and validated the previously identified sequence consensus surrounding TSSs in S.
<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /></figure>

**Figure 8.** Distribution of in-frame vs out-of-frame ATG codons and methionines for genes with or without transcription start sites downstream of annotated ATGs (iTSSs). (A) Cumulative total number of transcription start sites (TSSs) per nucleotide in *upf1Δ* cells for TSSs directly followed by the annotated ATGs (blue line) or a upstream open reading frame (uORF) (red line) and for TSSs downstream the annotated ATGs (iTSSs) and followed by an out-of-frame ATG (red dashed line) or by an in-frame ATG (blue dashed line). (B) Proportion of in-phase ATGs (+1 frame) following an annotated start codon, binned over nine nucleotides, for protein-coding genes with (red; 3327 genes) or without (blue; 1904 genes) iTSS reads in *upf1Δ*, as defined in the text. The dashed line indicates the expected value for a random distribution. Genes with iTSSs reads are significantly depleted of in-frame ATGs relative to all ATGs in the first 100 nucleotides when compared with genes without iTSSs reads (p-value = 3.13 10⁻⁵; ANOVA). (C) Proportion (log₂) along the ORFs of in-frame codons, normalized over the regions downstream the first 300 nucleotides of genes, for all codons (binned over nine nucleotides) for the two sets of genes defined in **B**. The ATG codon is in red. (D) Frequency of methionines per amino acid position along the ORFs for the two sets of genes defined in **B**. The p-values (ANOVA method) for the difference in methionine composition over the regions between 0–100 and 100–1000 nucleotides are indicated on the figure.

DOI: 10.7554/eLife.06722.020
cerevisiae (A(N)₆PyPu; Figure 1—figure supplement 2A and Figure 4—figure supplement 2, and Zhang and Dietrich 2005). Importantly, this consensus sequence was observed for the TSSs of all the different classes of transcripts identified in this study. Note that 1,775,474 sites in the genome conform to this consensus but only 85,714 (4.8%) were used as TSSs.

We also identified cases in which transcription started on what appears to be a 'slippery' sequence (PyAAA), which resulted in the incorporation of an additional non-encoded A at the 5'-end of the transcripts. The simplest explanation for this observation is that the very short nascent RNA formed after two or three nucleotide incorporation is able to shift back by one nucleotide, probably together with the RNAPII, before transcription resumes (Figure 1—figure supplement 2D), suggesting that the ternary complex formed by the transcribing polymerase, the RNA moiety and chromatin is rather labile at this early stage of transcription.

## Role of nonsense-mediated RNA decay in removal of pervasive and cryptic transcripts

Our analyses of TSSs in NMD-deficient cells revealed the prevalent role of this cytoplasmic RNA quality-control pathway in eliminating pervasive transcripts. We found that 52% of the SUTs and XUTs were significantly stabilized in *upf1Δ* cells when using the TSSCs from the TSS-sequencing experiments (see *Supplementary file 1*). When using the RNAseq approach, this proportion was even higher (70%; see *Supplementary file 2*), possibly because of the overall higher sequence read counts for a given transcript, which provides more power to statistical analyses. Since NMD is a translation-coupled RNA degradation pathway, the stabilizing effect seen upon *UPF1* deletion on pervasive transcripts indicated that once these transcripts reached the cytoplasm they associate with the translation machinery. Indeed, once they reached the cytoplasm the pervasive transcripts (being capped and poly-adenylated) cannot be distinguished from normal mRNAs and associate with the translation machinery. However, since they usually have short ORFs followed by relatively long 3' UTR regions (Figure 3—figure supplement 1), a hallmark of NMD substrates in yeast (*Decourty et al., 2014*), they will be targeted for degradation by this RNA quality-control pathway. This is supported by the identification of pervasive transcripts (mainly XUTs and SUTs) in ribosomal profiling experiments performed in wild-type cells and the recent discovery of a new class of unannotated transcripts (uRNAs) associated with polyribosomes and encoding short peptides (*Ingolia et al., 2009; Smith et al., 2014*).

The association of so-called 'non-coding' transcripts with the translation machinery and their ensuing high sensitivity to NMD has also been observed in multicellular organisms. In mouse embryonic stem cells, many long non-coding RNAs (lncRNAs) were found to be associated with translating ribosomes and an estimated 17.4% of lncRNAs were found to be upregulated in the absence of *UPF1*, compared with only 4% for protein-coding genes (*Ingolia et al., 2014; Smith et al., 2014*).

Our results also showed that XUTs, originally described as Xrn1-sensitive transcripts (van Dijk et al., 2011), are in fact primarily targeted by NMD by virtue of their poor coding capacities, with the Xrn1 exonuclease acting as a downstream effector of this pathway. Thus, we propose to redefine the various classes of pervasive transcripts identified in yeast into transcripts primarily degraded in the cytoplasm (SUTs and XUTs) and transcripts more sensitive to the nuclear exosome (CUTs) even though each of these degradation pathways similarly affect some transcripts and can show a synergic effect (see below).

Our analysis also revealed the prevalent role of NMD in the removal of cryptic transcripts initiating from intragenic promoters either in sense (here called B) or in antisense orientation (here called A). The strong impact of NMD inactivation on these transcripts indicated that they are exported to the cytoplasm and associate with the translation machinery, as is the case for the pervasive transcripts (see above). Even though we cannot rule out an indirect effect of the *UPF1* deletion on chromatin structure, the high number of intragenic transcripts (whether sense or antisense) identified in the *upf1Δ* single mutant suggests that the inhibition of transcription initiation by transcription-coupled chromatin modifications is not fully efficient and thus that some of these transcripts are produced as part of the normal transcription cycle in wild-type cells, with NMD playing an important role to eliminate them.

## Cooperation between nuclear and cytoplasmic RNA quality-control pathways

Deletions of *UPF1* and *RRP6* had an additive effect on the accumulation of CUTs, SUTs and XUTs. However, a number of pervasive transcripts, in particular 'A' and 'intergenic' ones, were stabilized to
substantial levels only when both the nuclear and the cytoplasmic RNA quality-control pathways were compromised, suggesting that they can act synergistically.

Two mechanisms of transcription termination have been reported in yeast. One depends on the cleavage and poly-adenylation complex (CPF-CFI/II), which generates the poly-adenylated mRNAs that get exported and translated in the cytoplasm, and another, which involves the NNS complex and is shared by CUTs and sn(o)RNAs precursors (see 'Introduction'). However, the demarcation between the two modes of termination is far from being strict as some terminators can often be recognized by both pathways depending on their distance from the TSS (Porrua et al., 2012). In the early phase of transcription elongation the RNAPII CTD repeats are mainly phosphorylated at Ser5, which favor recruitment of the NNS complex; while transcription proceeds, Ser2 gets phosphorylated at the expense of Ser5, promoting the recruitment of the CPF-CFI/II complex (Ahn et al., 2004; Kim et al., 2004). For some pervasive transcripts, these two termination pathways may compete and generate transcripts terminated by the NNS pathways and degraded by the nuclear exosome, as well as transcripts terminated by the CPF-CFI/II pathway and exported to the cytoplasm where they are targeted for degradation by the NMD. The produced RNAs would thus accumulate to detectable levels only when both RNA quality-control pathways are inactivated. Since the nuclear exosome has recently been shown to act with the NNS complex to promote early transcription termination at specific targets in S. cerevisiae (Fox et al., 2015), we cannot rule out that in wild-type or *upf1Δ* cells transcription will normally be terminated by the NNS complex and the synthesized transcripts degraded by the nuclear exosome, while in *rrp6Δ* cells transcription will proceed until the polymerase encounters the next CPF-CFI/II termination signal giving rise to longer transcripts exported to the cytoplasm and targeted for degradation by the NMD.

## NMD removes numerous transcripts arising from the low specificity of transcription initiation

NMD was previously shown to target a few hundred uORF-containing mRNAs (He et al., 2003; Guan et al., 2006; Johansson et al., 2007; Arribere and Gilbert, 2013). However, the use of NMD-deficient cells allowed us to identify 22% of the active genes for which a fraction of their transcript isoforms carried at least one uORF and were significantly sensitive to NMD. Moreover, 35% of the expressed genes generated transcript isoforms initiating at iTSSs and significantly stabilized in the absence of Upf1 (see Figure 7A,B). Some transcripts initiated at iTSSs were previously described in wild-type cells and shown to allow the synthesis of N-terminal variants of proteins exhibiting differential stabilities or localizations (Wu and Tzagoloff, 1987; Gammie et al., 1999; Arribere and Gilbert, 2013; Pelechano et al., 2013) but their number was vastly underestimated, probably because of their sensitivity to NMD (see Figure 7C). These two observations reveal an important role for NMD in getting rid of numerous undesired transcripts, the majority of which likely resulted from the low specificity of TSS selection by RNAPII. Yet, we cannot rule out that some of these transcripts have a biological function. For example, the use of alternative TSSs giving rise to transcripts with very different sensitivity to NMD could be used for regulatory purposes, in a way similar to that described for some genes of the nucleotide biosynthetic pathway (Kuehner and Brow, 2008; Thiebaut et al., 2008).

Altogether, transcript isoforms carrying a uORF or starting at an iTSS and targeted by NMD constituted the major transcript isoforms for only 446 mRNAs and minor transcript isoforms for almost half of the expressed genes. Therefore, although not quantitatively impacting the overall mRNA levels substantially (Figures 2A, 7C), NMD affects qualitatively a large fraction of the mRNA transcription units. NMD also has a strong impact on the accumulation of pervasive and cryptic transcripts and thus appears as a major player shaping the yeast transcriptome. In addition, we observed that transcripts initiated at iTSSs are significantly depleted in ATGs in the +1 frame, precluding the synthesis of N-terminally truncated proteins and ensuring their efficient NMD degradation. It thus suggests that this phenomenon might be important enough to have imposed an evolutionary constraint on the methionine content of the N-terminus of the yeast proteins.

## Materials and methods

### Yeast strains and culture

All the strains are a derivative of BY4741 and were obtained directly from the Euroscarf deletion collection (http://web.uni-frankfurt.de/fb15/mikro/euroscarf/) or generated by crossing with a can1Δ
derivative of BY4741 (LMA1057, see Table 2). Cells were grown to mid-exponential phase in YPD-rich medium at 30°C in a microturbidostat as previously described (Decourty et al., 2008), harvested by centrifugation and the pellet was frozen in liquid nitrogen.

## RNA extraction and analysis

Total RNA was extracted using the hot acid phenol protocol (Collart and Oliviero, 2001). Poly(A)⁺-RNA were obtained by two successive rounds of purification using oligo (dT)₂₅ magnetic beads (New England Biolabs, Ipswich, MA) following the manufacturer's protocol. Northern blots were carried out on poly(A)⁺-RNA as described in Neil et al. (2009) using ³²P-labeled riboprobes except for ACT1 for which a ³²P-labeled oligonucleotide was used.

## Library preparation

Approximately 500 ng of poly(A)⁺-RNA was mixed with 10 units of Antarctic phosphatase (New England Biolabs) in a final volume of 50 µl. After 1 hr at 37°C, the reaction was treated with phenol/chloroform and ethanol precipitated. The RNA pellet was resuspended in 44 µl of water and 10 units (1 µl) of TAP (Epicentre, Madison, WI) and 5 µl of 10× TAP buffer was added. The reaction was incubated 1 hr at 37°C followed by phenol/chloroform extraction and ethanol precipitation. The RNA pellet was resuspended in 5 µl of water. Next, the RNAs were ligated overnight at 16°C with 50 pmoles of the biotinylated oligonucleotide 3041 (Table 3) in a 20 µl reaction containing 10 units (1 µl) of T4 RNA ligase I (New England Biolabs) and ATP at a final concentration of 1 mM. RNAs were subsequently fragmented by incubation for 10 min at 70°C after addition of 5 µl of a 50 mM ZnCl₂, 50 mM Tris-HCl pH7.4 solution. The reaction was stopped by the addition of 1 µl EDTA 0.5 M and biotinylated RNA were purified using streptavidin magnetic beads according to the manufacturer's protocol (Dynabeads, MyOne streptavidin C1, Life Technologies, Carlsbad, CA). After washing, the beads were resuspended in 20 µl of water and the bound RNAs eluted by incubation for 5 min at 90°C. This fraction is enriched for 5'-ends of capped RNA molecules. Note that the supernatant of the first step of the purification procedure containing RNA fragments corresponding to the body and the 3'-end of the gene not attached to the biotinylated oligonucleotide can be recovered and used to prepare independent libraries. This RNA population can be further fractionated using oligo (dT)₂₅ magnetic beads to enrich for 3'-end of RNA molecules. The ~18.5 µl eluate from the streptavidin beads was mixed with 50 pmoles of oligonucleotide 3038 (see Table 3), heat denatured for 5 min at 70°C and slowly cooled down to 30°C in a Biorad iCycler PCR machine. Once the temperature had reached 30°C, 6 µl of 5× RT buffer, 1.5 µl of a 10 mM dNTPs solution, 1.5 µl of RNasin (Promega), 180 ng of actinomycin D, 300 units of RevertAid reverse transcriptase (Thermo Scientific, Waltham, MA), and water qsp 30 µl were added. The reaction was incubated 10 min at 30°C, followed by 40 min at

<table><caption>**Table 2.** Yeast strains used in this study</caption><thead><tr><td>Strain</td><td>Genotype</td><td>Reference</td></tr></thead><tbody><tr><td>BY4741</td><td>Mat a, his3Δ1, ura3Δ0, leu2Δ0, met15Δ0</td><td>(Brachmann et al., 1998)</td></tr><tr><td>LMA1057/3401</td><td>BY4741 can1Δ</td><td>This study</td></tr><tr><td>LMA1774/2759</td><td>BY4741 can1Δ, upf1Δ::HIS3MX6</td><td>This study</td></tr><tr><td>LMA1676/3405</td><td>BY4741 can1Δ, rrp6Δ::hphMX6</td><td>This study</td></tr><tr><td>LMA1772</td><td>BY4741 can1Δ, upf1Δ::KANMX6, rrp6Δ::HPHMX6</td><td>This study</td></tr><tr><td>LMA1790</td><td>BY4741 can1Δ, upf1Δ::KANMX6</td><td>This study</td></tr><tr><td>LMA2758</td><td>BY4741 can1Δ</td><td>This study</td></tr><tr><td>LMA2760</td><td>BY4741 can1Δ, xrn1Δ::KANMX6</td><td>This study</td></tr><tr><td>LMA2762</td><td>BY4741 can1Δ, xrn1Δ::KANMX6, upf1Δ::HIS3MX6</td><td>This study</td></tr><tr><td>LMA2921/3403</td><td>BY4741 can1Δ, set2Δ::KANMX6</td><td>This study</td></tr><tr><td>LMA2922</td><td>BY4741 can1Δ, set2Δ::KANMX6, upf1Δ::HIS3MX6</td><td>This study</td></tr><tr><td>LMA3409</td><td>BY4741 can1Δ, upf1Δ::HIS3MX6, rrp6Δ::HPHMX6</td><td>This study</td></tr></tbody></table>

DOI: 10.7554/eLife.06722.021
<table><caption>**Table 3.** Oligonucleotides used in this study</caption><thead><tr><th>Name</th><th>Sequence 5'-3'</th><th></th></tr></thead><tbody><tr><td>ACT1-1407-</td><td>ACACTTGTGGTGAACGATAGATGG</td><td>P<sup>32</sup> labelled probe</td></tr><tr><td>YMR114c-839+</td><td>ATCGAGGTGTAAAGGGTG</td><td>Synthesis of probe</td></tr><tr><td>T7-YMR114C-1068-*</td><td><b>TAATACGACTCACTATAGGG</b>CCTCTGGAGTCTTTCTGG</td><td>Synthesis of probe</td></tr><tr><td>YIL136w-1013+</td><td>ACTGGTGGTCTGGATGG</td><td>Synthesis of probe</td></tr><tr><td>T7-YIL136w(+)+115-*</td><td><b>TAATACGACTCACTATAGGG</b>TGCCACTAATTTACTCCG</td><td>Synthesis of probe</td></tr><tr><td>NEL025c-35+</td><td>AACAAATGCCAAGTCGGGAC</td><td>Synthesis of probe</td></tr><tr><td>T7-NEL025c-263-*</td><td><b>TAATACGACTCACTATAGGG</b>AAACGTTTGGTAAGAACTC</td><td>Synthesis of probe</td></tr><tr><td>SUT093_fwd</td><td>GAGTCCAGCGTCTCTACAC</td><td>Synthesis of probe</td></tr><tr><td>T7-SUT093_rev*</td><td><b>TAATACGACTCACTATAGGG</b>GACTTAATTGTCGTTGCTAGGAC</td><td>Synthesis of probe</td></tr><tr><td>SUT338_fwd</td><td>GAAAGACCGAAGGTGAAGAG</td><td>Synthesis of probe</td></tr><tr><td>T7-SUT338_rev*</td><td><b>TAATACGACTCACTATAGGG</b>TGGTACAGCCCTGTGTTCC</td><td>Synthesis of probe</td></tr><tr><td>SUT779_fwd</td><td>AACGAGGGAACTAGCCAG</td><td>Synthesis of probe</td></tr><tr><td>T7-SUT779_rev*</td><td><b>TAATACGACTCACTATAGGG</b>CTCTTCATCATCTGTGGAG</td><td>Synthesis of probe</td></tr><tr><td>TPO2(+)131-</td><td>GTATGTAGAAATGTCCGACG</td><td>Synthesis of probe</td></tr><tr><td>T7-TPO2-1798+*</td><td><b>TAATACGACTCACTATAGGG</b>GTAAGGGCTTGAGAC</td><td>Synthesis of probe</td></tr><tr><td>MAL12/32-1723-</td><td>GATTCTACCTTCCCATGG</td><td>Synthesis of probe</td></tr><tr><td>T7-MAL12/32-1161+*</td><td><b>TAATACGACTCACTATAGGG</b>TCAAGGTCAGGAGATAGG</td><td>Synthesis of probe</td></tr><tr><td>XUT3F5-fwd</td><td>AGGAAAAATGGGACTACAG</td><td>Synthesis of probe</td></tr><tr><td>T7-XUT3F5-rev*</td><td><b>TAATACGACTCACTATAGGG</b>TGTAAAAGGGCACAGTC</td><td>Synthesis of probe</td></tr><tr><td>3041†</td><td>5BioTEG/CTTTCCCTACACGACGCTCTTCCGATCT<br><b>NNNNCGCGCrGrGrNrN</b></td><td>Ligation with TAP<br>treated RNA</td></tr><tr><td>3118†</td><td>5BioTEG/CTTTCCCTACACGACGCTCTTCCGATCT<br><b>NNNNGCCGrCrGrNrN</b></td><td>Ligation with<br>fragmented RNA</td></tr><tr><td>3038†</td><td>GTTCAGACGTGTGCTCTTCCGATCTNNNNNN</td><td>Reverse transcription</td></tr></tbody></table>

42°C, 10 min at 55°C, 10 min at 60°C and 15 min at 75°C. RNAs were then degraded by incubation for 10 min at 75°C after the addition of 3 µl of 1 N NaOH. The reaction was quenched by the addition of 3 µl of 1 N HCl and precipitated by the addition of 3 volumes of 100% ethanol and 0.1 volume of 3 M sodium acetate pH 5.2. The precipitated cDNA was subjected to 5 cycles of PCR amplification in 20 µl with Phusion DNA polymerase (Thermo Scientific) using Illumina (San Diego, CA) multiplexing PCR primer 1.0 and 2.0 followed by an additional 6 to 8 cycles of amplification using Illumina multiplexing PCR primer 1.0 and one PCR primer index. The reaction was purified with Agencourt AmPure XP beads (Beckman Coulter, Indianapolis, IN) following manufacturer's instructions at a 1.8 x concentration and eluted in 20 µl water. Libraries were quantified using Qubit (Life Technologies) and sequenced on an Illumina HiSeq 2000 or 2500 with 50 or 100 bases single-end reads. To validate our method for the identification of TSSs, we compared libraries made with or without treatment of the RNAs with TAP. 500 ng of poly(A)<sup>+</sup>-RNA obtained by purification with oligo (dT)<sub>25</sub> magnetic beads starting from S. cerevisiae or S. pombe (used as an internal reference for normalization) were dephosphorylated using Antarctic phosphatase as described. Then, S. pombe RNAs and one half of the RNAs extracted from S. cerevisiae were treated with TAP. The second half of S. cerevisiae RNAs was mock treated. Before ligation with the biotinylated oligonucleotide 3041, each half of S. cerevisiae RNAs was mixed with an equivalent quantity of calf intestinal alkaline phosphatase (CIP)/ TAP-treated RNAs from S. pombe. The subsequent steps of the library preparation were identical to the ones described above.

<footnote>*The sequence in bold face corresponds to the T7 promoter sequence.</footnote>

<footnote>†The sequence in bold face corresponds to the tag used to identify the 5' end of the cDNAs. r stands for ribonucleotide.</footnote>

<footnote>DOI: 10.7554/eLife.06722.022</footnote>
For analysis of the *xrn1*1Δ mutant compared with the *upf1*1Δ mutant, since the absence of Xrn1 leads to the accumulation of decapped RNAs (*Hsu and Stevens, 1993*), we had to modify the protocol used for library preparation. Instead of the CIP/TAP treatment before ligation with the biotinylated oligonucleotide 3118, RNA was fragmented with ZnCl₂ and subsequently phosphorylated with T4 polynucleotide kinase (Thermo scientific). The subsequent steps of library preparation were identical to the ones described above for the 5'-end. Moreover, since the absence of Xrn1 affects the overall mRNA content (*Sun et al., 2013*), an aliquot of a *S. pombe* culture was added to the cell pellet before RNA extraction for library preparation to provide an independent internal control for normalization.

# Data analyses

## Illumina reads treatments

Duplicated reads, identified using the random sequence tags within the ligated oligonucleotides (oligonucleotides 3041 and 3118 in Table 3) were first filtered out. Then reads corresponding to the 5'-ends of cDNA fragments were extracted using the tag present in the added oligonucleotide (see Table 3) for all the unique reads. After removal of the tag, the resulting reads were mapped using bowtie (version 2.2.3 with the following parameters: -N 1 -p 1 -no-unal -D 15 -R 2 -L 22 -I S,1,1.15) and a compilation of *S. cerevisiae* genome (S288C reference sequence, Release 64 obtained from the Saccharomyces Genome Database (SGD) [http://www.yeastgenome.org/]) and *S. pombe* genome (ASM294 reference sequence, v2.19 obtained from PomBase [http://www.pombase.org/]) as reference genomes.

## Mapped reads processing

For libraries L5p_01 to L5p_10 (see Table 1) used for TSS sequencing, the 5'-end positions of the resulting mapped reads were used as TSS positions and extracted to wig files. For samples of the libraries LT_01 and LT_02, used for RNAseq, reads corresponding to the whole transcripts and full read coverage were extracted to wig files.

## Peak calling

Libraries L5p_01 to L5p_03 and L5p_04 to L5p_06 were regrouped into clusters (TSSCs) using the peak calling method described in *Neil et al. (2009)*. Parameters were optimized by maximizing the number of ORFs assigned to only one TSSC while minimizing the number of ORFs not assigned to any TSSC. As parameters, we finally chose a threshold of four reads in at least one of the samples used for the clustering and a maximum distance between two consecutive TSSs within a cluster of 50 nucleotides. For library L5p_01 to L5p_03, TSSCs were generated with the code TSSC_upf1_rrp6 (*Supplementary file 1*). For library L5p_04 to L5p_06, TSSCs were generated with the code TSSC_upf1_set2 (*Supplementary file 3*).

## Differential expression

TSSCs and transcript differential expression were calculated using DESeq2 (*Love et al., 2014*) within the SARTools pipeline (https://github.com/PF2-pasteur-fr/SARTools). For TSSCs, the mean TSS read counts in the region 50 to 10 nucleotides upstream of the ATG of each non-dubious ORF were used as an internal standard for size factor determination. For transcripts, the *S. pombe* transcripts read counts were used as the internal standard for size factor determination.

## Identification of TSSCs

TSSCs were assigned to annotations and sorted in a defined class according to their relative positions to the linked features. TSSCs overlapping the 5'-UTR mRNA sequences in the sense orientation were assigned to mRNAs (O-TSSCs in *Supplementary files 1, 3*). TSSCs mapped within the ORF or the 3'-UTR sequences and in the sense orientation were called B-TSSCs. TSSCs overlapping an mRNA sequence but in the antisense orientation were called A-TSSCs. C, S, X, and F TSSCs (*Supplementary files 1, 3*) correspond to TSSCs overlapping in the sense orientation with the 5'-end or found within CUTs, SUTs, XUTs and stable ncRNAs and transposable elements, respectively. All other TSSCs were classified as intergenic (I in *Supplementary files 1, 3*). CUTs, XUTs, and SUTs constituting overlapping transcript populations, when a TSSC was assigned to a pervasive transcript annotated in more than one of these classes we arbitrarily associated the corresponding TSSC in priority to CUTs, then to XUTs and finally to SUTs. mRNA coordinates were extracted from data of *Pelechano et al. (2013)*. ORFs, tRNAs, rRNAs, small nuclear RNAs (snRNAs), and sn(o)RNAs coordinates were retrieved from the SGD (http://www.yeastgenome.org/). CUT and SUT coordinates were retrieved from *Xu et al. (2009)*, while XUT
coordinates were from van Dijk et al. (2011). We also used the CUTs described in Neil et al. (2009).
ORF ATG coordinates were corrected for 150 of them identified as misannotated by Park et al. (2014).

Filtering of TSSCs

False positive TSSCs were filtered out using two criteria: the signal to noise ratio (SNR) and TSS accuracy.

SNR was defined as the log2 of the ratios of read counts of TAP samples vs no-TAP samples. For TSSC_upf1_rrp6 we used libraries L5p_07 and L5p_08 and for TSSC_upf1_set2 we used libraries L5p_09 and L5p_10. We choose a SNR limit of 0.55.

TSS accuracy was defined as the proportion of TSSs within a TSSC mapped on a purine preceded by a pyrimidine vs the total count of TSSs mapped within the TSSC. We choose a minimum TSS accuracy of 80%.

Distribution of individual TSSs around start codons

Individual TSS counts and read counts were used to determine the repartition of the TSSs around the ATGs of annotated ORFs (–200 to +100 nucleotides around ORFs ATGs). We used the ORF ATG positions retrieved from the SGD and corrected as described above in the 'Identification of TSSCs' section.

Identification of ORFs following individual TSSs

For each TSS identified in the region from –200 to +100 nucleotides around ORF ATGs, we looked for the two subsequent ATGs. When the TSS was upstream of the annotated ATG start codon, we determined whether or not it was the first encountered one. When the TSS was within the ORF, we determined whether the ATGs were in the same phase as the natural ORF or not.

Determination of ORF and 3'-UTR sizes for the different features

For mRNAs, ORF start, ORF end and 3'-end coordinates obtained, as described above in the 'Identification of TSSCs' section, were used to calculate the sizes of the ORF and the 3'-UTR. For XUTs and SUTs, ORF start was determined as the coordinate of the first ATG after the annotated start obtained in Xu et al. (2009) and van Dijk et al. (2011). ORF ends were determined as the coordinate of the first in phase stop codon encountered after the previously selected ATG start codon. The region between the end of this potential ORF and the end of the XUT or SUT was considered to be the 3'UTR.

Nucleotide use

We determined the proportion of A, T, G, and C at each position in the region spanning –200 to +100 nucleotides around the ATGs of ORFs according to the annotation found in the SGD and corrected as described above.

Codons and amino acids use

We determined the number of ATGs within the first 1000 nucleotides after the ORF start codon based on the information available in the SGD. To avoid the overrepresented value of the start ATG codon, the first three nucleotides were not included in the analysis. The ratios between phased ATG and total ATG counts were computed in a window of nine nucleotides. The same method was applied for all other codons. To estimate the differential use of a codon between the start and the body of the ORFs, we normalized the value for each codon taking into account the mean ratio for this particular codon in the region comprised between 500 and 1000 nucleotides after the start codon.

We used the protein sequences retrieved from SGD and omitted the first methionine of the proteins.

Transcripts expression levels

We used the coverage of LT_01 and LT_02 libraries to calculate transcript expression levels of mRNAs, SUTs and XUTs. The coordinates of the different types of transcripts were obtained as described above. We calculated the expression levels for the four samples of both libraries and then used DESeq2 for normalization, estimation of fold change and mean signal calculation (Supplementary file 2).

Mean NFR density

The mean NFR density has been calculated using a nucleosome definition in Kaplan et al. (2009). We collected the nucleosome signal in regions from –1000 to +1000 nucleotides around each TSS, aligned on TSS position and calculated means for each position.

Accession number

The data reported here have been deposited in NCBI GEO under the accession number GSE64139.
# Acknowledgements

We thank D Libri, B Séraphin and M Fromont-Racine for discussions and for critical reading of the manuscript.

# Additional information

## Funding

<table><thead><tr><td>Funder</td><td>Grant reference</td><td>Author</td></tr></thead><tbody><tr><td>Institut Pasteur</td><td></td><td>Christophe Malabat</td></tr><tr><td>Centre National de la Recherche Scientifique</td><td></td><td>Alain Jacquier</td></tr><tr><td>Agence Nationale de la Recherche</td><td>ANR-2011-BSV6-011-02</td><td>Alain Jacquier</td></tr><tr><td>Agence Nationale de la Recherche</td><td>ANR-14-CE-10-0014-01</td><td>Alain Jacquier</td></tr></tbody></table>

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

## Author contributions

CM, CS, Analysis and interpretation of data, Drafting or revising the article; FF, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article; LM, Acquisition of data, Analysis and interpretation of data; AJ, Conception and design, Analysis and interpretation of data, Drafting or revising the article

# Additional files

## Supplementary files

* Supplementary file 1. TSS clusters identified in WT, *upf1Δ*, *rrp6Δ* and *upf1Δrrp6Δ* cells.
DOI: 10.7554/eLife.06722.024

* Supplementary file 2. Reads count for whole transcriptome analysis in wild type, *xrn1Δ*, *upf1Δ* and *upf1Δxrn1Δ* cells.
DOI: 10.7554/eLife.06722.025

* Supplementary file 3. TSS clusters identified in WT, *upf1Δ*, *set2Δ* and *set2Δupf1Δ* cells.
DOI: 10.7554/eLife.06722.026

* Supplementary file 4. List of genes with identified iTSSs and/or uORFs.
DOI: 10.7554/eLife.06722.027

## Major datasets

The following dataset was generated:

<table><thead><tr><td>Author(s)</td><td>Year</td><td>Dataset title</td><td>Dataset ID and/or URL</td><td>Database, license, and accessibility information</td></tr></thead><tbody><tr><td>Feuerbach F, Malabat C, Ma L, Saveanu C, Jacquier A</td><td>2015</td><td>Quality control of transcription start site selection by nonsense-mediated-mRNA Decay</td><td>http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE64139</td><td>Publicly available at NCBI Gene Expression Omnibus (GSE64139).</td></tr></tbody></table>

The following previously published datasets were used:

<table><thead><tr><td>Author(s)</td><td>Year</td><td>Dataset title</td><td>Dataset ID and/or URL</td><td>Database, license, and accessibility information</td></tr></thead><tbody><tr><td>Xu Z, Wei W, Gagneur J, Perocchi F, Clauder-Münster S, Camblong J, Guffanti E, Stutz F, Huber W, Steinmetz LM</td><td>2009</td><td>Transcription profiling of wild type yeast grown with ethanol, glucose and galactose and the deletion mutant of Rrp6 to identify transcription start and end positions</td><td>http://www.ebi.ac.uk/arrayexpress/experiments/E-TABM-590/</td><td>Publicly available at ArrayExpress (E-TABM-590).</td></tr></tbody></table>
<table><thead><tr><th>Author(s)</th><th>Year</th><th>Dataset title</th><th>Dataset ID and/or URL</th><th>Database, license, and accessibility information</th></tr></thead><tbody><tr><td>Arribere JA, Gilbert WV</td><td>2013</td><td>Roles for Transcript Leaders in Translation and mRNA Decay Revealed by Transcript Leader Sequencing</td><td>http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE39074</td><td>Publicly available at NCBI Gene Expression Omnibus (GSE39074).</td></tr></tbody></table>

## References

* Ahn SH, Kim M, Buratowski S. 2004. Phosphorylation of serine 2 within the RNA polymerase II C-terminal domain couples transcription and 3' end processing. *Molecular Cell* **13**:67-76. doi: 10.1016/S1097-2765(03)00492-1.

* Almada AE, Wu X, Kriz AJ, Burge CB, Sharp PA. 2013. Promoter directionality is controlled by U1 snRNP and polyadenylation signals. *Nature* **499**:360-363. doi: 10.1038/nature12349.

* Arigo JT, Eyler DE, Carroll KL, Corden JL. 2006. Termination of cryptic unstable transcripts is directed by yeast RNA-binding proteins Nrd1 and Nab3. *Molecular Cell* **23**:841-851. doi: 10.1016/j.molcel.2006.07.024.

* Arribere JA, Gilbert WV. 2013. Roles for transcript leaders in translation and mRNA decay revealed by transcript leader sequencing. *Genome Research* **23**:977-987. doi: 10.1101/gr.150342.112.

* Brachmann CB, Davies A, Cost GJ, Caputo E, Li J, Hieter P, Boeke JD. 1998. Designer deletion strains derived from *Saccharomyces cerevisiae* S288C: a useful set of strains and plasmids for PCR-mediated gene disruption and other applications. *Yeast* **14**:115-132. doi: 10.1002/(SICI)1097-0061(19980130)14:2<115::AID-YEA204>3.0.CO;2-2.

* Collart MA, Oliviero S. 2001. Preparation of yeast RNA. In: Ausubel FM, Brent R, Kingston RE, Moore DD, Seidman JG, Smith JA, Struhl K, editors. *Current protocols in molecular biology*. Hoboken, NJ, USA: John Wiley & Sons.

* Crooks GE, Hon G, Chandonia JM, Brenner SE. 2004. WebLogo: a sequence logo generator. *Genome Research* **14**:1188-1190. doi: 10.1101/gr.849004.

* Davis CA, Ares M Jr. 2006. Accumulation of unstable promoter-associated transcripts upon loss of the nuclear exosome subunit Rrp6p in *Saccharomyces cerevisiae*. *Proceedings of the National Academy of Sciences of USA* **103**:3262-3267. doi: 10.1073/pnas.0507783103.

* Decourty L, Doyen A, Malabat C, Frachon E, Rispal D, Séraphin B, Feuerbach F, Jacquier A, Saveanu C. 2014. Long open reading frame transcripts escape nonsense-mediated mRNA decay in yeast. *Cell Reports* **6**:593-598. doi: 10.1016/j.celrep.2014.01.025.

* Decourty L, Saveanu C, Zemam K, Hantraye F, Frachon E, Rousselle J-C, Fromont-Racine M, Jacquier A. 2008. Linking functionally related genes by sensitive and quantitative characterization of genetic interaction profiles. *Proceedings of the National Academy of Sciences of USA* **105**:5821-5826. doi: 10.1073/pnas.0710533105.

* Djebali S, Davis CA, Merkel A, Dobin A, Lassmann T, Mortazavi A, Tanzer A, Lagarde J, Lin W, Schlesinger F, Xue C, Marinov GK, Khatun J, Williams BA, Zaleski C, Rozowsky J, Röder M, Kokocinski F, Abdelhamid RF, Alioto T, Antoshechkin I, Baer MT, Bar NS, Batut P, Bell K, Bell I, Chakrabortty S, Chen X, Chrast J, Curado J, Derrien T, Drenkow J, Dumais E, Dumais J, Duttagupta R, Falconnet E, Fastuca M, Fejes-Toth K, Ferreira P, Foissac S, Fullwood MJ, Gao H, Gonzalez D, Gordon A, Gunawardena H, Howald C, Jha S, Johnson R, Kapranov P, King B, Kingswood C, Luo OJ, Park E, Persaud K, Preall JB, Ribeca P, Risk B, Robyr D, Sammeth M, Schaffer L, See L-H, Shahab A, Skancke J, Suzuki AM, Takahashi H, Tilgner H, Trout D, Walters N, Wang H, Wrobel J, Yu Y, Ruan X, Hayashizaki Y, Harrow J, Gerstein M, Hubbard T, Reymond A, Antonarakis SE, Hannon G, Giddings MC, Ruan Y, Wold B, Carninci P, Guigó R, Gingeras TR. 2012. Landscape of transcription in human cells. *Nature* **489**:101-108. doi: 10.1038/nature11233.

* Fox MJ, Gao H, Smith-Kinnaman WR, Liu Y, Mosley AL. 2015. The exosome component Rrp6 is required for RNA polymerase II termination at specific targets of the Nrd1-Nab3 pathway. *PLOS Genetics* **10**:e1004999. doi: 10.1371/journal.pgen.1004999.

* Gammie AE, Stewart BG, Scott CF, Rose MD. 1999. The two forms of karyogamy transcription factor Kar4p are regulated by differential initiation of transcription, translation, and protein turnover. *Molecular and Cellular Biology* **19**:817-825.

* Geisler S, Lojek L, Khalil AM, Baker KE, Coller J. 2012. Decapping of long noncoding RNAs regulates inducible genes. *Molecular Cell* **45**:279-291. doi: 10.1016/j.molcel.2011.11.025.

* Guan Q, Zheng W, Tang S, Liu X, Zinkel RA, Tsui K-W, Yandell BS, Culbertson MR. 2006. Impact of nonsense-mediated mRNA decay on the global expression profile of budding yeast. *PLOS Genetics* **2**:e203. doi: 10.1371/journal.pgen.0020203.

* Gudipati RK, Villa T, Boulay J, Libri D. 2008. Phosphorylation of the RNA polymerase II C-terminal domain dictates transcription termination choice. *Nature Structural & Molecular Biology* **15**:786-794. doi: 10.1038/nsmb.1460.

* Hashimoto S, Qu W, Ahsan B, Ogoshi K, Sasaki A, Nakatani Y, Lee Y, Ogawa M, Ametani A, Suzuki Y, Sugano S, Lee CC, Nutter RC, Morishita S, Matsushima K. 2009. High-resolution analysis of the 5'-end transcriptome using a next generation DNA sequencer. *PLOS ONE* **4**:e4108. doi: 10.1371/journal.pone.0004108.

* He F, Li X, Spatrick P, Casillo R, Dong S, Jacobson A. 2003. Genome-wide analysis of mRNAs regulated by the nonsense-mediated and 5' to 3' mRNA decay pathways in yeast. *Molecular Cell* **12**:1439-1452. doi: 10.1016/S1097-2765(03)00446-5.

* Hsu CL, Stevens A. 1993. Yeast cells lacking 5'->3' exoribonuclease 1 contain mRNA species that are poly(A) deficient and partially lack the 5' cap structure. *Molecular and Cellular Biology* **13**:4826-4835.
Ingolia NT, Brar GA, Stern-Ginossar N, Harris MS, Talhouarne GJS, Jackson SE, Wills MR, Weissman JS. 2014. Ribosome profiling reveals pervasive translation outside of annotated protein-coding genes. *Cell Reports* **8**: 1365-1379. doi: 10.1016/j.celrep.2014.07.045.

Ingolia NT, Ghaemmaghami S, Newman JRS, Weissman JS. 2009. Genome-wide analysis in vivo of translation with nucleotide resolution using ribosome profiling. *Science* **324**:218-223. doi: 10.1126/science.1168978.

Johansson MJO, He F, Spatrick P, Li C, Jacobson A. 2007. Association of yeast Upf1p with direct substrates of the NMD pathway. *Proceedings of the National Academy of Sciences of USA* **104**:20872-20877. doi: 10.1073/pnas.0709257105.

Kaplan N, Moore IK, Fondufe-Mittendorf Y, Gossett AJ, Tillo D, Field Y, LeProust EM, Hughes TR, Lieb JD, Widom J, Segal E. 2009. The DNA-encoded nucleosome organization of a eukaryotic genome. *Nature* **458**:362-366. doi: 10.1038/nature07667.

Kawashima T, Douglass S, Gabunilas J, Pellegrini M, Chanfreau GF. 2014. Widespread use of non-productive alternative splice sites in *Saccharomyces cerevisiae*. *PLOS Genetics* **10**:e1004249. doi: 10.1371/journal.pgen.1004249.

Kervestin S, Jacobson A. 2012. NMD: a multifaceted response to premature translational termination. *Nature Reviews Molecular Cell Biology* **13**:700-712. doi: 10.1038/nrm3454.

Kim M, Ahn S-H, Krogan NJ, Greenblatt JF, Buratowski S. 2004. Transitions in RNA polymerase II elongation complexes at the 3' ends of genes. *The EMBO Journal* **23**:354-364. doi: 10.1038/sj.emboj.7600053.

Kuehner JN, Brow DA. 2008. Regulation of a eukaryotic gene by GTP-dependent start site selection and transcription attenuation. *Molecular Cell* **31**:201-211. doi: 10.1016/j.molcel.2008.05.018.

Leeds P, Peltz SW, Jacobson A, Culbertson MR. 1991. The product of the yeast UPF1 gene is required for rapid turnover of mRNAs containing a premature translational termination codon. *Genes & Development* **5**:2303-2314. doi: 10.1101/gad.5.12a.2303.

Love MI, Huber W, Anders S. 2014. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. *Genome Biology* **15**:550. doi: 10.1186/s13059-014-0550-8.

Mackiewicz P, Kowalczuk M, Gierlik A, Dudek MR, Cebrat S. 1999. Origin and properties of non-coding ORFs in the yeast genome. *Nucleic Acids Research* **27**:3503-3509. doi: 10.1093/nar/27.17.3503.

Neil H, Malabat C, d'Aubenton-Carafa Y, Xu Z, Steinmetz LM, Jacquier A. 2009. Widespread bidirectional promoters are the major source of cryptic transcripts in yeast. *Nature* **457**:1038-1042. doi: 10.1038/nature07747.

Ntini E, Järvelin AI, Bornholdt J, Chen Y, Boyd M, Jørgensen M, Andersson R, Hoof I, Schein A, Andersen PR, Andersen PK, Preker P, Valen E, Zhao X, Pelechano V, Steinmetz LM, Sandelin A, Jensen TH. 2013. Polyadenylation site-induced decay of upstream transcripts enforces promoter directionality. *Nature Structural & Molecular Biology* **20**:923-928. doi: 10.1038/nsmb.2640.

Ozsolak F, Kapranov P, Foissac S, Kim SW, Fishilevich E, Monaghan AP, John B, Milos PM. 2010. Comprehensive polyadenylation site maps in yeast and human reveal pervasive alternative polyadenylation. *Cell* **143**:1018-1029. doi: 10.1016/j.cell.2010.11.020.

Park D, Morris AR, Battenhouse A, Iyer VR. 2014. Simultaneous mapping of transcript ends at single-nucleotide resolution and identification of widespread promoter-associated non-coding RNA governed by TATA elements. *Nucleic Acids Research* **42**:3736-3749. doi: 10.1093/nar/gkt1366.

Pelechano V, Wei W, Steinmetz LM. 2013. Extensive transcriptional heterogeneity revealed by isoform profiling. *Nature* **497**:127-131. doi: 10.1038/nature12121.

Pokholok DK, Harbison CT, Levine S, Cole M, Hannett NM, Lee TI, Bell GW, Walker K, Rolfe PA, Herbolsheimer E, Zeitlinger J, Lewitter F, Gifford DK, Young RA. 2005. Genome-wide map of nucleosome acetylation and methylation in yeast. *Cell* **122**:517-527. doi: 10.1016/j.cell.2005.06.026.

Porrua O, Hobor F, Boulay J, Kubicek K, D'Aubenton-Carafa Y, Gudipati RK, Stefl R, Libri D. 2012. In vivo SELEX reveals novel sequence and structural determinants of Nrd1-Nab3-Sen1-dependent transcription termination. *The EMBO Journal* **31**:3935-3948. doi: 10.1038/emboj.2012.237.

Sayani S, Chanfreau GF. 2012. Sequential RNA degradation pathways provide a fail-safe mechanism to limit the accumulation of unspliced transcripts in *Saccharomyces cerevisiae*. *RNA* **18**:1563-1572. doi: 10.1261/rna.033779.112.

Sayani S, Janis M, Lee CY, Toesca I, Chanfreau GF. 2008. Widespread impact of nonsense-mediated mRNA decay on the yeast intronome. *Molecular Cell* **31**:360-370. doi: 10.1016/j.molcel.2008.07.005.

Schulz D, Schwalb B, Kiesel A, Baejen C, Torkler P, Gagneur J, Soeding J, Cramer P. 2013. Transcriptome surveillance by selective termination of noncoding RNA synthesis. *Cell* **155**:1075-1087. doi: 10.1016/j.cell.2013.10.024.

Smith JE, Alvarez-Dominguez JR, Kline N, Huynh NJ, Geisler S, Hu W, Coller J, Baker KE. 2014. Translation of small open reading frames within unannotated RNA transcripts in *Saccharomyces cerevisiae*. *Cell Reports* **7**: 1858-1866. doi: 10.1016/j.celrep.2014.05.023.

Smolle M, Workman JL. 2013. Transcription-associated histone modifications and cryptic transcription. *Biochimica Et Biophysica Acta* **1829**:84-97. doi: 10.1016/j.bbagrm.2012.08.008.

Sun M, Schwalb B, Pirkl N, Maier KC, Schenk A, Failmezger H, Tresch A, Cramer P. 2013. Global analysis of eukaryotic mRNA degradation reveals xrn1-dependent buffering of transcript levels. *Molecular Cell* **52**:52-62. doi: 10.1016/j.molcel.2013.09.010.

Thiebaut M, Colin J, Neil H, Jacquier A, Séraphin B, Lacroute F, Libri D. 2008. Futile cycle of transcription initiation and termination modulates the response to nucleotide shortage in *S. cerevisiae*. *Molecular Cell* **31**:671-682. doi: 10.1016/j.molcel.2008.08.010.

Thiebaut M, Kisseleva-Romanova E, Rougemaille M, Boulay J, Libri D. 2006. Transcription termination and nuclear degradation of cryptic unstable transcripts: a role for the nrd1-nab3 pathway in genome surveillance. *Molecular Cell* **23**:853-864. doi: 10.1016/j.molcel.2006.07.029.
Tudek A, Porrua O, Kabzinski T, Lidschreiber M, Kubicek K, Fortova A, Lacroute F, Vanacova S, Cramer P, Stefl R, Libri D. 2014. Molecular basis for coordinating transcription termination with noncoding RNA degradation. *Molecular Cell* **55**:467-481. doi: 10.1016/j.molcel.2014.05.031.

van Dijk EL, Chen CL, d' Aubenton-Carafa Y, Gourvennec S, Kwapisz M, Roche V, Bertrand C, Silvain M, Legoix-Né P, Loeillet S, Nicolas A, Thermes C, Morillon A. 2011. XUTs are a class of Xrn1-sensitive antisense regulatory non-coding RNA in yeast. *Nature* **475**:114-117. doi: 10.1038/nature10118.

Vasiljeva L, Kim M, Mutschler H, Buratowski S, Meinhart A. 2008. The Nrd1-Nab3-Sen1 termination complex interacts with the Ser5-phosphorylated RNA polymerase II C-terminal domain. *Nature Structural & Molecular Biology* **15**:795-804. doi: 10.1038/nsmb.1468.

Waern K, Snyder M. 2013. Extensive transcript diversity and novel upstream open reading frame regulation in yeast. *G3* **3**:343-352. doi: 10.1534/g3.112.003640.

Wu M, Tzagoloff A. 1987. Mitochondrial and cytoplasmic fumarases in *Saccharomyces cerevisiae* are encoded by a single nuclear gene FUM1. *The Journal of Biological Chemistry* **262**:12275-12282.

Wyers F, Rougemaille M, Badis G, Rousselle J-C, Dufour ME, Boulay J, Régnault B, Devaux F, Namane A, Séraphin B, Libri D, Jacquier A. 2005. Cryptic pol II transcripts are degraded by a nuclear quality control pathway involving a new poly(A) polymerase. *Cell* **121**:725-737. doi: 10.1016/j.cell.2005.04.030.

Xu Z, Wei W, Gagneur J, Perocchi F, Clauder-Münster S, Camblong J, Guffanti E, Stutz F, Huber W, Steinmetz LM. 2009. Bidirectional promoters generate pervasive transcription in yeast. *Nature* **457**:1033-1037. doi: 10.1038/nature07728.

Zhang Z, Dietrich FS. 2005. Mapping of transcription start sites in *Saccharomyces cerevisiae* using 5' SAGE. *Nucleic Acids Research* **33**:2838-2851. doi: 10.1093/nar/gki583.