The Modifier of Transcription 1 (Mot1) ATPase and Spt16 Histone Chaperone Co-regulate Transcription through Preinitiation Complex Assembly and Nucleosome Organization*

Modifier of transcription 1 (Mot1) is a conserved and essential Swi2/Snf2 ATPase that can remove TATA-binding protein (TBP) from DNA using ATP hydrolysis and in so doing exerts global effects on transcription. Spt16 is also essential and functions globally in transcriptional regulation as a component of the facilitates chromatin transcription (FACT) histone chaperone complex. Here we demonstrate that Mot1 and Spt16 regulate a largely overlapping set of genes in Saccharomyces cerevisiae. As expected, Mot1 was found to control TBP levels at co-regulated promoters. In contrast, Spt16 did not affect TBP recruitment. On a global scale, Spt16 was required for Mot1 promoter localization, and Mot1 also affected Spt16 localization to genes. Interestingly, we found that Mot1 has an unanticipated role in establishing or maintaining the occupancy and positioning of nucleosomes at the 5′ ends of genes. Spt16 has a broad role in regulating chromatin organization in gene bodies, including those nucleosomes affected by Mot1. These results suggest that the large scale overlap in Mot1 and Spt16 function arises from a combination of both their unique and shared functions in transcription complex assembly and chromatin structure regulation.


Modifier of transcription 1 (Mot1) is a conserved and essential
Swi2/Snf2 ATPase that can remove TATA-binding protein (TBP) from DNA using ATP hydrolysis and in so doing exerts global effects on transcription. Spt16 is also essential and functions globally in transcriptional regulation as a component of the facilitates chromatin transcription (FACT) histone chaperone complex. Here we demonstrate that Mot1 and Spt16 regulate a largely overlapping set of genes in Saccharomyces cerevisiae. As expected, Mot1 was found to control TBP levels at co-regulated promoters. In contrast, Spt16 did not affect TBP recruitment. On a global scale, Spt16 was required for Mot1 promoter localization, and Mot1 also affected Spt16 localization to genes. Interestingly, we found that Mot1 has an unanticipated role in establishing or maintaining the occupancy and positioning of nucleosomes at the 5 ends of genes. Spt16 has a broad role in regulating chromatin organization in gene bodies, including those nucleosomes affected by Mot1. These results suggest that the large scale overlap in Mot1 and Spt16 function arises from a combination of both their unique and shared functions in transcription complex assembly and chromatin structure regulation.
Other studies have suggested that Mot1 may have a direct role in transcription activation, including facilitating transcription preinitiation complex (PIC) formation or remodeling chromatin (17)(18)(19)(20). Mot1 regulates over one-third of the protein-coding genome in budding yeast (17,(21)(22)(23), affecting both mRNA levels and RNA synthesis precision. Using a novel method to determine transcription length defects, we showed that the most common length defect in the temperature-sensitive mot1-42 strain is premature termination, linking Mot1 to productive transcription elongation as well as initiation (23). Chromatin immunoprecipitation (ChIP) experiments showed that Mot1 is present mainly at promoters (15,18,21), raising the question of how it can promote elongation.
Spt16 is a histone chaperone with functions in transcription, replication, and DNA repair (24 -28). It forms a heterodimer with Pob3, and together they loosely associate with the high mobility group protein Nhp6 to bind nucleosomes (29 -32). This complex of Spt16, Pob3, and Nhp6 is known as facilitates chromatin transcription (FACT) or yFACT in budding yeast. Spt16 can bind H2A/H2B dimers in vivo and in vitro and H3/H4 tetramers in vitro but with higher specificity for H2A/H2B (30,(33)(34)(35)(36)(37). Recently it was demonstrated that the middle domain of human Spt16 interacts with the H3/H4 tetramer, and the AID domain interacts with H2B (38). The best studied transcriptional role for Spt16 involves the catalysis of nucleosomal reorganization during transcriptional elongation (25, 39 -41). Spt16 has been suggested to travel with elongating RNA polymerase II (Pol II) (42,43), enabling it to remove or reorganize nucleosomes blocking transcription and replace them after Pol II passage (24,25,33,44). Consistent with this role, mutations in Spt16 led to exposure of adventitious promoters buried within the gene body, leading to a phenomenon known as cryptic transcription (45,46). There is some evidence that Spt16 can function during normal initiation as well by promoting transcription preinitiation complex formation (37,43,47). The functions of Mot1 and Spt16 in initiation as well as elongation prompted us to test whether they overlap functionally.

Mot1 and Spt16 Genetically and Physically Interact-Mot1
and Spt16 have been implicated in transcription initiation and elongation but through different mechanisms. However, a functional relationship between them was suggested by the identification of Spt16 as a Mot1-associated protein using mass spectrometry (48). To test for a functional relationship, we first compared growth of conditional single mutant strains with that of the mot1-42 spt16-197 double mutant. As shown in Fig. 1A, mot1-42 and spt16-197 strains were temperature-sensitive at 35°C. The double mutant displayed a synthetic growth defect at 30°C, indicating that MOT1 and SPT16 interact genetically. In reciprocal co-immunoprecipitation (co-IP) experiments, Spt16 and Mot1 co-associated in cell extracts (Fig. 1, B and C), substantiating the prior report of their physical association. Together, these results establish that the two factors interact both biochemically and genetically.
Genome-wide Analysis of RNA Expression Changes-To assess the global relationship between Mot1 and Spt16 in gene expression, we first compared mot1-42 genome-wide tiling array data from our previous study (23) with published spt16-197 genome-wide tiling array data (46). Initial analysis suggested overlap in genes regulated by Mot1 and Spt16 and significant correlation in expression changes between the two mutant strains (data not shown). The spt16-197 study used polyadenylated RNA and interrogated open reading frames (ORFs) with six probes per ORF for the main purpose of detecting cryptic initiation (46). In contrast, our mot1-42 study used total RNA and probes at 8-base pair (bp) intervals throughout the entire yeast genome to measure expression changes at higher resolution (23). For a more equivalent comparison, we repeated the spt16-197 genome-wide analysis using total RNA and probes at 8-bp intervals. log 2 median differential expression correlated highly between previous and new data sets as expected ( Fig. 2A; r ϭ 0.72). Our initial observation of the correlation between the mot1-42 data set and the previous spt16-197 data set was upheld by the new spt16-197 data set, showing that the expression changes in mot1-42 and spt16-197 correlate ( Fig. 2B; r ϭ 0.34).
We next identified genes for which changes in expression were significant. For each data set, upper and lower thresholds of deviation from a normal distribution were defined to classify genes as repressed (up-regulated in a mutant), activated (downregulated in a mutant), or unaffected (no substantial change) (Fig. 2, D, E, and F). Approximately 40% of genes classified as activated or repressed in one mutant were classified similarly in the other (Fig. 3A). The observed number of co-regulated genes (1,307 of 6,685 genes in the yeast genome) is much higher than chance expectation and indicates significant overlap in transcriptional regulation. Additionally, considering only genes significantly affected by Mot1 and Spt16, the initial correlation in Fig. 2B increases substantially (r ϭ 0.52). To determine the effects in the double mutant, we analyzed RNA from the mot1-42 spt16-197 double mutant obtained from a cross of the mot1-42 and spt16-197 strains. Because the spt16-197 strain was derived from a different background than the mot1-42 strain, we first computed two sets of differential expression by comparing the double mutant with each parental wild type (WT) strain. The two sets correlated highly, indicating that differences in strain background contribute relatively little to the differential gene expression changes (Fig. 2C). For this reason, subsequent analyses for the double mutant are in comparison with the mot1-42 parental WT strain. We observed that differential expression was similar between single mutants, and in the double mutant these effects were exacerbated (Fig.  3B). The results show that at co-regulated genes Mot1 and Spt16 act in combination to control gene expression.
To validate the differential effects inferred from the tiling array data and to select specific examples of co-regulated genes for further study, total RNA was isolated from each mutant, and gene-specific RNA levels were measured by RT-qPCR. Consistent with the tiling array data, RPS23B, TRX1, LEU9 and ACT1 all showed decreased expression in the three mutant strains, and most decreases were significant using a p value threshold of 0.05 (Fig. 3C). Three examples of genes identified by the global analysis as repressed by Mot1 and Spt16 were also confirmed by qPCR, and in most cases these were statistically significant (Fig. 3D).   . p values were determined by a one-tailed Fisher's exact test. B, comparison of expression in the single and double mutants. The heat map shows log 2 median differential expression for 6,685 genes with vertical lines as genes and colored as in the key. Genes regulated by both factors generally showed similar expression changes in the single mutants that were exacerbated in the double mutant. C, validation of gene expression effects at selected co-activated genes. The bar graph shows the relative levels of RNA for each gene relative to the level in the WT strain. The RNA levels for the Mot1 and Spt16 co-activated genes RPS23B, TRX1, LEU9, and ACT1 decreased compared with their levels in WT cells. PCF11 is an unaffected control gene. D, RNA levels for the Mot1 and Spt16 co-repressed genes INO1, GAD1, and GND2. In C and D, error bars show one standard deviation from at least three biological replicates, and asterisks denote p Ͻ 0.05 by Student's t test.

Mot1 and Spt16 Co-Regulation
Preferences for Promoter Elements in Mot1 and Spt16 Regulation-Mot1 preferentially represses genes containing the TATA consensus sequence (12,15,18,23,49) "TATAWAWR" (W ϭ A/T and R ϭ A/G) (50), so we tested whether Spt16 had a similar preference. A preference for regulation of TATA-containing genes by Spt16 would be consistent with the production of cryptic transcripts that initiate preferentially from sites containing a cryptic TATA element in spt16 cells (46). Most promoters in yeast historically classified as lacking a TATA element do contain a short degenerate TATA sequence that underlies the TBP binding location (51). Here, for simplicity, genes with TATA sequences in their promoters are termed "TATA-containing," and genes without high affinity TBP binding sites are referred to as "TATA-less." As expected for the mot1-42 data set, TATA-containing promoters were enriched in the repressed gene class, and TATA-less promoters were enriched in the activated gene class (Fig. 4A). Notably, spt16-197 and the double mutant strains had nearly identical enrichment patterns to the mot1-42 strain. In contrast, no such pattern was observed in a deletion strain for Set2, a histone methyltransferase that suppresses cryptic initiation (52), underscoring the significance of the preferential regulation for TATA or TATA-less promoters by both Mot1 and Spt16.
Mot1 and Spt16 Affect Recruitment of PIC Components-Mot1 and Spt16 have been shown to be involved in TBP recruitment to promoters (15,37,53), but their functional relationship to PIC assembly at co-regulated promoters has not been reported. Consistent with published ChIP data (17, 20 -22), in mot1-42 cells TBP levels increased at five of six Mot1-regulated promoters examined (Fig. 4B). TBP levels were unaffected in spt16-197 cells at all of these co-regulated promoters (Fig. 4B). In the mot1-42 spt16-197 double mutant, TBP levels increased at all co-regulated promoters and to a similar extent as in mot1-42 cells (Fig. 4B). These results suggest that Mot1 influences TBP recruitment at these promoters, whereas Spt16 does not. ChIP was also performed for TFIIB, which associates with TBP-DNA and forms a scaffold for association of Pol II (54). In contrast to TBP, TFIIB levels generally correlated with gene expression changes observed in each mutant ( Fig. 4C) (23). These results suggest that the increased TBP occupancy observed at Mot1-activated promoters in mot1-42 cells is indicative of transcriptionally incompetent complexes. The results also suggest that Mot1 and Spt16 facilitate formation of a relatively early promoter-bound intermediate that promotes TFIIB binding but that they do so through distinct biochemical mechanisms with respect to TBP.
Roles for Mot1 and Spt16 in RNA Synthesis Precision-Given the roles for Mot1 and Spt16 in transcription initiation, we next compared their roles in RNA synthesis precision. Using our method described previously (23), tiling array data were mined for RNA length defects defined by four categories: upstream initiation (UI), downstream or cryptic initiation (DI), premature termination (PT), and downstream termination (DT). Fig.  5 shows the relationship between the change in level and precision of gene expression in each mutant strain. Previous studies have shown that characteristic defects in mot1-42 and spt16-197 cells are PT and DI, respectively (23,46); however, we observed that spt16-197 RNA was highly enriched in all four defect categories relative to mot1-42 RNA, which indicates that Spt16 has a broader role in synthesis precision than reported previously. The mot1-42 spt16-197 double mutant RNA was also enriched in all four categories, and DI and PT were both prominent. The co-prominence of the main defect from each single mutant in the double mutant suggests distinct roles for Mot1 and Spt16 in RNA synthesis precision. Common to all three strains were significant enrichments in two of the four possible initiation/termination couplings (UI/DT and DI/PT), average differential expression near zero for DI and PT, and widely distributed and on average highly elevated differential expression for UI and DT. Despite non-equivalent roles for FIGURE 4. Mot1 and Spt16 have similar preferences for promoter elements but regulate PIC components differently. A, genes are grouped by promoter attribute and proportioned by category of differential expression: co-repressed (yellow), co-activated (blue), or unaffected (gray). The set2⌬ data (23) are a control that demonstrates a different relationship of TATA-containing and TATA-less genes to a regulatory factor. p values were calculated using a 2 test. p values associated with enrichment or depletion of each group of genes are indicated by the symbols (* and #) as defined in the table below the figure. B, to examine effects of Mot1 on TBP levels at co-regulated promoters, relative TBP levels in WT and mutant strains were examined at three co-activated and three co-repressed genes. TBP in mot1-42 cells increased at promoters regardless of promoter attribute, was unaffected in spt16-197 cells, and in the double mutant cells resembled mot1-42 cells. C, TFIIB recruitment was affected in mot1-42 and spt16-197 cells at co-regulated promoters. Changes in TFIIB levels correlated with changes in RNA in each mutant. In B and C, error bars are one standard deviation from at least three biological replicates, and asterisks indicate p Ͻ 0.05.
Mot1 and Spt16 in elongation, certain characteristics of each individual and coupled defect category are the same for both factors, suggesting a shared role that extends beyond the amount of RNA production to how RNA synthesis precision is determined.
Spt16 Affects Mot1 Localization Genome-wide-To better define the functional relationship between Mot1 and Spt16, we next determined the genome-wide localization profile of Mot1 by ChIP-Seq and compared the patterns in WT and spt16-197 cells, focusing on genes whose expression was regulated by both factors. Mot1 distributions at co-activated genes (Fig. 6A) and co-repressed genes (Fig. 6B) were plotted by computationally aligning genes to their transcription start site (TSS). On average, Mot1 ChIP signal was highest in the 5Ј regulatory region as expected (9,15,18,21). At both classes of genes, promoterlocalized Mot1 signals decreased in spt16-197 cells compared with WT with the largest effect occurring at co-activated gene promoters. Mot1 ChIP-Seq signal was also higher at co-activated genes than co-repressed genes (Fig. 6, black lines), consistent with published results showing that Mot1 signal in general scales with gene expression level (21,55). The global effect of Spt16 on the Mot1 profile supports a broad role for Spt16 in establishing the proper chromatin environment for Mot1 binding. Mutation of Spt16 also resulted in a decrease in Mot1 signal downstream of the TSS of co-activated genes (Fig. 6A, red line) and an increase at co-repressed genes (Fig. 6B, red line). Possible explanations for these differences are discussed below.
Mot1 Affects Spt16 Localization-For comparison with the Mot1 ChIP-Seq data just described, we performed ChIP-Seq for Spt16 in WT and mot1-42 cells. Spt16 profiles were plotted for Mot1-Spt16 co-activated genes (Fig. 7A) and co-repressed genes (Fig. 7B). Average Spt16 localization differed from Mot1 with enrichment throughout gene bodies and depletion from promoter regions. Like Mot1, Spt16 levels increased at co-activated genes and decreased at co-repressed genes in mot1-42 cells compared with WT (Fig. 7, A and B). The relative enrichment of Spt16 at co-activated compared with co-repressed genes in WT cells suggests that Spt16 directly activates its target genes, whereas Spt16-mediated repression may result from indirect effects. Spt16 enrichment in gene bodies is consistent with published work as expected for a factor with key roles in elongation (25,(41)(42)(43)(44)(45)(46). Global trends in Spt16 enrichment inferred from ChIP-Seq data were validated for seven genes by locus-specific ChIP and qPCR (Fig. 7C). Notably, Mot1 was required for full Spt16 recruitment to co-activated promoters, and Mot1 affected recruitment of Spt16 to two of three co-repressed promoters. From these results, it is unclear whether Mot1 plays a direct role in Spt16 recruitment or whether Spt16 occupancy is linked to gene activity as dictated by Pol II density.
Effects of Mot1 on RNA Polymerase II Levels-To better understand the relationship between Mot1 and Spt16 in transcription initiation and elongation, Pol II levels were measured by ChIP at a subset of co-activated and co-repressed genes. As shown in Fig. 8A, Pol II levels were affected in mot1-42 versus WT cells as expected based on the changes we observed in RNA levels at these same genes. Changes in Pol II levels were statistically significant and quantitatively in line with the relatively modest changes in gene expression (Fig. 3). A larger difference in Pol II was detectable at the GAL10 ORF during galactose induction (Fig. 8B), indicating that these modest changes are not due to insensitivity in the ChIP assay. Combined with the effects on RNA synthesis precision described above,  JULY 15, 2016 • VOLUME 291 • NUMBER 29

Mot1 and Spt16 Co-Regulation
although Mot1 and Spt16 are essential global transcriptional regulators, they appear to have primary roles in fine-tuning gene expression rather than being absolutely required for RNA synthesis. However, as cells have the capacity to buffer effects on key regulators of transcription by modifying the RNA decay rate (56 -58), it may be that the analysis of steady-state RNA reported here underestimates the magnitude of the effects of Mot1 and Spt16 on RNA synthesis.
Roles for Mot1 and Spt16 in Nucleosome Organization-Next, we used micrococcal nuclease (MNase) ChIP-Seq to define the chromatin landscape in WT and mutant cells. Mononucleosomes were obtained by MNase digestion followed by histone H3 immunoprecipitation, DNA isolation, and high throughput sequencing as described previously (see "Experimental Procedures" and Ref. 59). As shown in Fig. 9A (green line), in an average plot of all genes in WT cells, we observed a nucleosome-depleted promoter region flanked by peaks of nucleosomal density. As expected, nucleosomal peaks were most precisely defined proximal to the TSS and became progressively less well defined within the gene body and moving away from the TSS (60 -62). Nucleosome density and periodicity were diminished in spt16-197 cells, consistent with prior results implicating Spt16 in establishing chromatin structure on a global scale (32, 33, 39 -41, 45, 46). In addition, nucleo-  ChIP signal decreased at co-activated genes and increased at co-repressed genes, consistent with changes observed in A and B and RNA levels (Fig. 3, C and D). Error bars are one standard deviation from the mean from at least three biological replicates, and asterisks denote p Ͻ 0.05. nucleosomes, but other nucleosomes in the transcribed region were largely unaffected. As discussed below, the positions of the ϩ1 and ϩ2 nucleosomes appear shifted in mot1-42 cells com-pared with their average positions in WT cells, and promoter regions show somewhat more pronounced nucleosomal depletion in cells lacking Mot1. Nucleosome organization in mot1-42 spt16-197 cells was more severely perturbed than in either single mutant with an overall pattern similar to spt16-197 cells but with a greater loss of nucleosome density, particularly at the 5Ј ends of genes. Nucleosome density in promoter regions was restored to near WT levels in the double mutant compared with spt16-197 cells as well.
Similar trends for each mutant were observed at co-activated genes (Fig. 9B) compared with the all-gene average. Mot1 exerted effects on the nucleosomes proximal to the TSS, and there were more widespread effects on nucleosome structure at co-activated genes in spt16-197 and double mutant strains. At co-repressed genes (Fig. 9C), mutation of Mot1 led to diminished nucleosome occupancies at the 5Ј ends of genes, there was global disruption of nucleosome structure in spt16-197 cells, and the most severe effects were observed in mot1-42 spt16-197 cells. These effects are consistent with nucleosome disruption accompanying increased expression of these genes in the mutants. However, there were trends at co-repressed genes that were markedly different from patterns at other classes of genes as well. Comparison of the mot1-42 pattern in Fig. 9C with A and B shows that nucleosome depletion at the promoter is attributable primarily to repressed genes. For spt16-197, elevated histone H3 levels at the promoter specifically were not observed at repressed genes. For all strains, there is a less pronounced trough between the ϩ1 and ϩ2 peaks, which may be related to chromatin organization that occurs with up-regulation. Additionally, although mutation of Mot1 by itself mainly affected the ϩ1 and ϩ2 nucleosomes regardless of gene expression outcome, comparison of spt16-197 with the double mutant reveals a more expansive co-regulatory role for Mot1 in gene repression with effects on nucleosome occupancy that extend at least several hundred base pairs in both directions from the TSS.
To better understand the nature of nucleosome changes that gave rise to the different patterns in Fig. 9, A-C, a computational approach (60, 63) was applied to identify nucleosomal peaks that show statistically significant changes in occupancy versus positioning (also known as fuzziness; see Refs. 60, 64, and 65). Occupancy was defined by the height of the MNase ChIP-Seq peak, whereas fuzziness refers to an estimate of the width of the peak. The statistically significant (false discovery rate-corrected p values Ͻ0.05) differential nucleosome effects observed in each mutant strain compared with WT are summarized in Table 1. In general, those nucleosomes with significant changes in occupancy showed decreased occupancy in a given mutant compared with the WT strain. In addition, there were more changes in fuzziness than occupancy with most affected nucleosomes having increased fuzziness in the mutant compared with WT. As shown in Fig. 10A, the overlaps in nucleosomes showing differential occupancy and fuzziness were statistically significant with the greatest overlaps occurring between differentially fuzzy nucleosomes in spt16-197 and the double mutant and in agreement with their similar nucleosomal patterns shown in Fig. 9. In mot1-42 cells, the features of the genic nucleosomal array are still present, but the first two nucleosomal peaks are reduced and slightly shifted rightward into the gene body. The Ϫ1 and Ϫ2 nucleosomes immediately upstream of the promoter region were also reduced on average compared with WT. B, plot as in A but for all co-activated genes. C, plot as in A but for all co-repressed genes.

Mot1 and Spt16 Co-Regulation
Mot1 and Spt16 Affect the Position of the ϩ1 Nucleosome-In addition to affecting occupancy and fuzziness, it appeared from the gene average plots in Fig. 9 that loss of one or both of these factors shifted the median position of promoter-proximal nucleosomes with respect to the TSS. Such an effect has been reported previously for Spt16 (65). Nucleosomal peaks were canonically designated Ϫ1, ϩ1, ϩ2, etc. based on the distance of the nucleosome to the TSS. Nucleosomal dyad positions were tabulated for each differentially fuzzy nucleosome, and position shifts were calculated from the difference in dyad position for each nucleosome in the mutant compared with WT. Position shift patterns for spt16-197 cells were similar to published data (65) with the largest (ϩ10-bp) shift at the ϩ1 nucleosome. A ϩ10-bp median shift at the ϩ1 nucleosome was also observed in mot1-42 cells and the double mutant, which demonstrates that the two factors play non-redundant roles in positioning the ϩ1 nucleosome properly (Fig. 10B). Nucleosomes at ϩ2 and ϩ3 positions also showed significant position shifts in all three mutants as did ϩ4 and ϩ5 nucleosomes in spt16-197 cells, consistent with the loss of nucleosome organization in this mutant and the double mutant.

Discussion
Mot1 and Spt16 Co-Regulate Transcription-Here we report a functional overlap in transcriptional control between Mot1 and Spt16 that operates on a global scale with ϳ1,300 genes controlled by both factors. Co-regulated genes have similar differential expression patterns in the single mutant strains, and in the double mutant these effects tend to be exacerbated. The physical and genetic interactions between Mot1 and Spt16 further suggest that these two factors function together in vivo. This relationship is particularly interesting because they have distinct genome-wide localizations (Figs. 6 and 7). Spt16 is broadly distributed on the transcribed regions of genes but notably depleted in promoter regions where Mot1 is enriched. The localization of Spt16 is distinct from the localization of other elongation factors such as Spt6 and Bur1 (66). Mutation of Mot1 led to redistribution of Spt16 in accordance with changes in expression, most likely due to changes in Pol II levels and the requirement for Spt16 during elongation. It is interesting that mutation of Spt16 decreased Mot1 localization at coregulated promoters regardless of whether expression was positively or negatively affected. This suggests that although Spt16 is not localized to promoters it is nonetheless required to establish the chromatin environment responsible for the genomewide distribution of Mot1.
Although they are found at distinct locations in the genome, the physical interaction between Mot1 and Spt16 in extracts suggests how their activities may be coordinated. The physical interaction with Mot1 may increase the local concentration of Spt16 near gene start sites, allowing Spt16 to more efficiently act on promoter-proximal nucleosomes and/or to be handed off to Pol II during the early stages of elongation (25,42,43). Although Mot1 and Spt16 associate with one another in extracts, this may not be a direct interaction, and indeed, an interaction mediated by other regulatory factors could help explain how the activities of these factors are coordinated when  Fig. 9 established that the ϩ1 nucleosome was located from approximately Ϫ50 to ϩ100 bp with respect to the TSS. Nucleosomes downstream of the ϩ1 nucleosome were assigned to non-overlapping 150-bp bins based on the ϩ1 nucleosome position. The distance from the nucleosome dyad to the TSS was then determined for each nucleosome and in each strain. The plot shows the distribution of median differences in position with respect to the TSS for each nucleosome in the indicated mutant compared with WT. The ϩ1, ϩ2, and ϩ3 nucleosomes in mot1-42 (blue), spt16-197 (red), and mot1-42 spt16-197 (green) cells had significant shifts away from the TSS. The Ϫ1, ϩ4, and ϩ5 nucleosomes all shifted away from the TSS in spt16-197 and double mutant cells. Asterisks indicate p Ͻ 0.05 as determined by a Kolmogorov-Smirnov test. NDR, nucleosomedepleted region.

TABLE 1 Effects on nucleosomes in mutants compared with wild type
Mot1 and Spt16 Co-Regulation they are largely found at different chromosomal locations. Regardless of the nature of their physical association, such a mechanism would be consistent with previous observations of Spt16 localization within or near promoters and affecting nucleosome occupancy to properly initiate transcription (37, 66 -68). The synthetic genetic interaction observed between Mot1 and Spt16 is consistent with the observation that mutations in Mot1 are synthetically lethal with deletions of the FACT subunit genes NHP6A and NHP6B (69). Overall, these synthetic effects strongly suggest that Mot1 and the entire FACT complex function together to coordinate PIC assembly and nucleosome organization.
Although the results presented here demonstrate broad overlap in the set of genes controlled by Mot1 and the set controlled by Spt16, there are many genes apparently controlled by only one or the other as well (Fig. 3A). The transcriptome-wide overlap that we observe is limited by the achievable statistical power typical of such studies; allele specificity, growth conditions, and other technical issues may also have led to an underestimate of the functional overlap of Mot1 and Spt16. Conversely, we would expect that, given the complexity in the chromatin landscape and in gene regulatory mechanisms, not all affected genes require both factors for WT expression.
Regulation of TBP-A previous study (37) found that Spt16 was important for establishing the levels of TBP at the SER3 and ELP3 promoters, so it is possible that Spt16 activity impacts TBP directly at some promoters but not all. In our analysis, we did not detect any effect on TBP promoter localization in spt16-197 cells, but increased TBP levels were detected at co-regulated promoters in mot1-42 and mot1-42 spt16-197 cells. This is consistent with the global role of Mot1 in limiting TBP binding to chromatin that we reported previously (23). In contrast, both Mot1 and Spt16 apparently control TFIIB levels at promoters, which are in general better correlated with gene expression levels than TBP (23). Because mutation of Spt16 promotes cryptic initiation (45,46), the genome-wide redistribution of Mot1 in spt16-197 cells reported here is consistent with Mot1 tracking the presumed redistribution of TBP to DNA sites exposed as a result of the global perturbation of chromatin integrity that occurs when Spt16 activity is compromised.
Mot1 and Spt16 Have Roles in Elongation-Mutations in Mot1 and Spt16 gave rise to prominent but distinct RNA length defects, and the double mutant showed the combination of defects present in the single mutants. These results suggest that Mot1 and Spt16 have distinct roles in elongation. However, they do show similar associations between length precision and differential expression. In particular, both factors overlap in suppressing upstream initiation and downstream termination at co-repressed genes (Fig. 5). One conundrum is why Mot1 and Spt16 are essential when their effects on steady-state RNA levels, although affecting many genes, are relatively small in magnitude. As mentioned above, changes in transcription may be underestimated as a result of the cellular capacity to buffer RNA turnover rates in the mutants. However, the large numbers of RNA length defects observed in the mot1-42 and spt16-197 strains suggest that the essential roles of these factors may be more intimately related to ensuring proper start site and termination site selection rather than driving abundance of proper RNA species per se.
Mot1 and Spt16 Globally Affect Nucleosome Organization-Spt16 has been widely studied as a histone chaperone, but the role of Mot1 in nucleosomal organization had not been established. Mot1 was required for nucleosome remodeling at the induced GAL1 promoter (19), but the chromatin organization at the Mot1-activated URA1 promoter was unaffected by a mutation in Mot1 (20). A direct role for Mot1 in chromatin remodeling is attractive as most other members of the Swi2/ Snf2 enzyme family are ATP-dependent chromatin remodelers (70 -72). The results presented here implicate Mot1 in determining the occupancy and position of the ϩ1 and ϩ2 nucleosomes. It remains to be determined precisely how these nucleosomal changes arise, but as they were observed at genes regardless of the effect of Mot1 on expression, it is unlikely that these changes are entirely an indirect result of Mot1-mediated effects on transcription via interaction of Mot1 with TBP. The main roles for Mot1 in chromatin structure involve positioning and fuzziness of the ϩ1 and ϩ2 nucleosomes. These nucleosomes tended to be less precisely positioned (more fuzzy) in mot1-42 cells than WT, and the ϩ1 and ϩ2 nucleosomes had a median position shift of 10 bp away from the TSS. A shift was also apparent in spt16-197 cells and the double mutant, indicating that both factors contribute to establishing and/or maintaining this important aspect of the native chromatin state.

Experimental Procedures
Yeast Strains and Growth Conditions-Saccharomyces cerevisiae strains used in this study are listed in Table 2. The WT strain AY51 is referred to as MOT1-WT because it was used to generate the mot1-42 (AY87) strain; these strains were described previously (73,74). FY56 is referred to as SPT16-WT Co-immunoprecipitation-Yeast strains were grown as above to an OD of ϳ1.0, and cells were then collected and lysed in BA/150 lysis buffer (20 mM HEPES-KOH, pH 7.6, 2 mM EDTA, 2 mM EGTA, 0.25% Nonidet P-40, 150 mM potassium acetate, 5 mM DTT, and a protease inhibitor mixture tablet (one tablet/25 ml of buffer) (Roche Applied Science)). A total of 1 mg of protein was immunoprecipitated overnight at 4°C with either 15 l of ␣-myc (9E10, 1 mg/ml) or 2.5 l ␣-Spt16 (rabbit polyclonal antibody from Tim Formosa). Immunoprecipitated samples were incubated with Protein A-Sepharose beads (GE Healthcare) and washed with BA/250 lysis buffer (same as BA/150 lysis buffer but with 250 mM potassium acetate). Samples were separated on 4 -15% precast Bio-Rad SDS-polyacrylamide gels. After transfer to nitrocellulose membranes, samples were blocked in Tris-buffered saline with Tween (50 mM Tris HCl, pH 7.6, 150 mM NaCl, 0.05% Tween 20) plus 5% nonfat milk, then incubated with either ␣-myc (9E10; 1 mg/ml), ␣-HA (12CA5, Abcam), or ␣-Spt16 (rabbit polyclonal antibody) followed by incubation with secondary antibody, and developed using the ECL Prime kit (GE Healthcare).
RNA Isolation for Tiling Array Analysis and RT-PCR-Strains were grown as above in duplicate for tiling array analysis and in triplicate for RT-PCR. Total RNA was isolated, and samples were prepared for tiling array analysis as described previously (23). The University of Virginia Microarray Core Facility hybridized the samples to S. cerevisiae 1.0 Tiling Arrays (Affymetrix) and generated the raw data. After the total RNA was isolated for RT-qPCR analysis, the RNA was reverse transcribed with the iScript Select cDNA Synthesis kit (Bio-Rad) according to the manufacturer's instructions with the exception that twice the amount of starting RNA and reagents were used. Real time PCR was carried out as described previously (23).
Tiling Array Data Analysis-Integrated Genome Browsercompatible files were produced from the raw tiling array data using Tiling Analysis Software. Signal profiles were smoothed over a 50-bp bandwidth (101-bp sliding window), and differen-tial profiles were smoothed over a 250-bp bandwidth (501-bp sliding window). Additional Tiling Analysis Software parameters were as follows: log 2 signal scale, scale to target intensity of 100, two-sided test type, plus or minus 0.3 threshold, 100-bp maximum gap, and 50-bp minimum run. The Mot1 differential signal was generated by subtracting the mot1-42 signal (log 2 scale) by the MOT1-WT signal. The Spt16 differential signal was generated by subtracting the spt16-197 signal by the SPT16-WT signal. Two double mutant differentials were generated by subtracting the mot1-42 spt16-197 signal by each of the two WT signals. Two types of analysis were conducted with the differential profiles: unbiased and biased. The unbiased analysis considered genome-wide differential expression, whereas the biased analysis considered only the 6,685 annotated genes. For the biased analysis, median differential expression was calculated for each gene, and for the 309 intron-containing genes only the longest exon was considered. Because each mutant produces different magnitudes of changes in transcript levels, different thresholds were applied to each differential profile for classifying genes as underexpressed or overexpressed. These thresholds were determined by comparing the distributions of median differential expression with a normal distribution and identifying the levels above or below which median gene expression deviated from the normal distribution using quantile-quantile plots. Because both WT strains had relatively high correlations between the data sets (r ϭ 0.74), further analyses for the single and double mutants were conducted using comparisons with MOT1-WT.
Chromatin Immunoprecipitation-Yeast strains were grown as described above. After cultures reached an OD of ϳ1.0 and were heat-shocked, cells were collected, lysed, and sonicated as described previously (20). Immunoprecipitations were conducted using 1 mg of protein as described previously with either 2.5 l of ␣-TFIIB (Sua7 rabbit polyclonal), 5 l of ␣-TBP (58C9, Abcam), 5 l of ␣-myc (9E10; 1 mg/ml), or 5 l of ␣-Pol II (8WG16, Covance). After immunoprecipitation and cross-link reversal, DNA was purified using a Qiagen PCR purification kit according to the manufacturer's instructions. Purified DNA was then quantified by qPCR using the oligonucleotides (Life Technologies) listed in Table 3 and SYBR Green Master Mix (Bio-Rad). The mock values were subtracted from the IP values and then normalized to the input.
Analysis of RNA Length Precision-Categories of RNA length defects were assigned to genes using our previously published method (23). Intron-containing genes were excluded from the analysis. For the spt16-197 and mot1-42 spt16-197 differential RNA signals compared with WT, the signal cutoff was raised from 0.30 to 0.36 to account for overall greater magnitudes of differential expression and to reduce the likelihood of calling false positives. The other parameters were held constant.
ChIP-Seq-Immunoprecipitated DNA was prepared as above for ChIP. The resulting purified DNA was then processed using the TruSeq ChIP Library Prep kit from Illumina (IP-202-1012) according to the manufacturer's instructions. Briefly, DNA fragments were end-repaired and adenylated, adapter 6 was ligated onto the adenylated ends and gel-purified, and purified products were enriched and quantified via Qubit. The DNA samples were then analyzed using a Bioanalyzer to ensure that the samples mostly contained a single peak of mononucleosomal DNA. A total of 0.2 pmol of DNA was sequenced by the University of Virginia Sequencing Core using an Illumina MiSeq instrument, yielding 25.6 Ϯ 1.6 million 150-bp raw reads per sample. Data were obtained from two independent replicate samples for each strain.
MNase ChIP-Seq-MNase ChIP-Seq was performed as described previously with minor modifications (59). Strains were grown as described above. After collection, cells were lysed as described above for ChIP. MNase titrations were performed, and samples with the highest proportion of mononucleosomes and without subnucleosomal fragments were selected for each sample. Cross-links were reversed with heat, and protein was digested with proteinase K (Thermo Fisher). Samples were then immunoprecipitated with 10 g of ␣-H3 (ab1791, Abcam) overnight at 4°C, then incubated with Protein A-Sepharose, washed, and eluted as described previously. The immunoprecipitated material was prepared for sequencing as described above, and sequencing data were obtained using the Illumina MiSeq instrument in the University of Virginia Sequencing Core, yielding 27.9 Ϯ 5.7 million 150-bp raw reads per sample. Data were obtained from two independent replicate samples for each strain.
ChIP-Seq and MNase ChIP-Seq Data Analysis-Reads were aligned to the SacCer3 reference genome using Bowtie 2 2.2.6 with default mapping settings (75). More than 95% of the raw reads from each data set were mapped under these conditions. Unmapped reads and reads with mapping quality scores Ͻ30 were removed from the Bowtie 2-generated bam file using SAMtools 0.1.19 (76). This mapping quality threshold resulted in removal of 10% or fewer of the mapped reads. For viewing the data in Integrative Genomics Viewer (IGV) (77, 78), bam files were converted to BigWig by first sorting and indexing using SAMtools, then creating .bed and .cov files using BEDTools 2.18.2 (79), and finally converting the .cov files to BigWig using the bedGraphToBigWig converter (80). Gene average plots were produced using ngs.plot 2.4.7 (81). Nucleosome mapping and differential nucleosome analyses were performed using the dpos tool in the DANPOS 2.2.2 package with default settings that include -fold change normalization (multiplying each sample by a scale factor to account for differences in mapped read coverage) rather than quantile normalization as originally described (63). dpos yielded normalized wig files of the processed data. We observed that the total signal nonetheless varied among these normalized data sets by a small (ϳ5%) but not insignificant amount. For this reason, the dpos-generated wig data files were further normalized globally based on the total signal. The gene average plots shown in Fig. 9 were generated in R from the renormalized wig data by computing the mean signal at base pair resolution for all genes, activated genes, or repressed genes (as indicated) with respect to the TSS. Other statistical analyses and data visualization were performed using R v3.0.2. All of the genome-wide data sets from this study have been submitted to the NCBI Gene Expression Omnibus (www.ncbi.nlm.nih.gov/geo) under SuperSeries accession number GSE80235.   -TGC TAA GCA CTA CCG CAT TG-3Ј  RPS23B Promoter R 5Ј-GAA AGC GTG GAG ACA AGG AG-3Ј  TRX1 Promoter F  5Ј-CCA AAA CCC TGA AAC TGC AT-3Ј  TRX1 Promoter R  5Ј-ATT CGC TGG CAG TTT TGA AT-3Ј  ACT1 Promoter F  5Ј-CAG CTT TTA GAT TTT TCA CGC TTA-