Genome-wide mapping of nucleosome positions in Saccharomyces cerevisiae in response to different nitrogen conditions.

Well-organized chromatin is involved in a number of various transcriptional regulation and gene expression. We used genome-wide mapping of nucleosomes in response to different nitrogen conditions to determine both nucleosome profiles and gene expression events in Saccharomyces cerevisiae. Nitrogen conditions influence general nucleosome profiles and the expression of nitrogen catabolite repression (NCR) sensitive genes. The nucleosome occupancy of TATA-containing genes was higher compared to TATA-less genes. TATA-less genes in high or low nucleosome occupancy, showed a significant change in gene coding regions when shifting cells from glutamine to proline as the sole nitrogen resource. Furthermore, a correlation between the expression of nucleosome occupancy induced NCR sensitive genes or TATA containing genes in NCR sensitive genes, and nucleosome prediction were found when cells were cultured in proline or shifting from glutamine to proline as the sole nitrogen source compared to glutamine. These results also showed that variation of nucleosome occupancy accompany with chromatin-dependent transcription factor could influence the expression of a series of genes involved in the specific regulation of nitrogen utilization.

Nucleosomes form the basic repeating unit of eukaryotic chromatin, consisting of ~147 bp DNA wrapped around a histone octamer core123. Dynamic changes of genome-wide maps of nucleosome positions reflect a set of molecular processes, including DNA replication, gene expression and transcription regulation45. Upon changes in response to different environmental conditions, >20% genes switch their promoter states6. DNA-binding transcription factors and chromatin-remodeling enzymes are primarily responsible for the nucleosome architecture6. Meanwhile, different nucleosome profiles could be observed under the same culture conditions while cells’ growth conditions are distinct7. Precise nucleosome positions on the genome inhibit the binding of any other protein, such as a transcriptional activator or the general transcriptional machinery, to the target DNA region. The change of nucleosome positions is crucial for understanding the mechanism underlying the regulation of gene expression8.

In response to different nutritional conditions or perturbation of the environment, yeast cells tend to fine tune their transcriptional, translational and metabolic processes in order to maximize survival ability9. For different nitrogen sources, including preferred nitrogen sources glutamate and glutamine or non-preferred nitrogen sources, such as proline, Saccharomyces cerevisiae is able to use the nitrogen source that enables best growth by a mechanism called nitrogen catabolite repression (NCR)101112. The NCR genes are regulated primarily by an interplay of the four GATA family transcription factors transcriptional activators Gln3 and Gat1 as well as repressors Dal80 and Gzf31314. More than 392 genes are involved in the response to changes of nitrogen source type and quality and about 90 genes are regulated directly by the four GATA regulators91415. These NCR-sensitive genes are involved in some metabolic processes, including amino acid and allantoin metabolism as well as transport and transcription regulation.

Recent advances in high-throughput sequencing technology have dramatically expanded our understanding of genome-wide nucleosome organization61617. The nucleosome core can be isolated as a particle containing 147 bp DNA by the digestion of chromatin with micrococcal nuclease (MNase)1819. MNase-seq is considered the standard procedure for nucleosome profile prediction620. Nucleosomes overlapping the transcription start site (TSS) of gene-coding regions are defined as +1 nucleosomes, while a nucleosome immediately before a +1 nucleosome is defined as a −1 nucleosome, the so called +1 and −1 nucleosomes flank a considerably larger opening termed the nucleosome-free region (NFR)6. Evidence is presented that the NFRs at TSS is essential for successful transcription initiation2122. A 3′ NFR present in >95% of all genes might be important in transcription termination;23 however, the dynamic change of nucleosome positions in response to different nitrogen conditions has not been examined on a genome-wide scale in S. cerevisiae. To map the location of individual nucleosomes on a genomic scale in response to different nitrogen conditions, we used a MNase-seq method to sequence the ends of nucleosome-associated DNA. This study focused firstly on the regions around TSSs and stop codon sites, which are well-defined nucleosome distributions2324. Furthermore, the nucleosome profiles of the TATA-containing and TATA-less promoters were analyzed to determine the correlation between general gene distribution and nucleosome profiles. Finally, the expression of NCR-sensitive genes was calculated to investigate their relationship under different nitrogen conditions.

Results

Nucleosome profiles of the 5′ and 3′ ends of all genes under different nitrogen conditions

Knowledge of the detailed remodeling of nucleosome positioning across genomes and the mechanism underlying nucleosome profile changes is crucial for understanding gene regulation and expression. The nucleosome profiles of the 5′ and 3′ ends of a gene-coding region have been considered as a marker for gene transcription initiation and termination. Here, the whole genome data of nucleosome occupancy profiles were exhibited at the 5′ and 3′ ends of all genes (Fig. 1A,B). There is a notable NFR immediately upstream of the TSS and downstream of the stop codon; the high nucleosome occupancy region is immediately downstream of the TSS and upstream of the stop codon. Therefore, the nucleosome occupancy score differed under different nitrogen conditions: Gln, low; Pro, medium; and Gln-Pro, high.

Figure 1

Patterns of nucleosome occupancy of S. cerevisiae under different nitrogen conditions.

(A) Nucleosome occupancy profiles of all genes in the yeast genome from −1000 bp to +1000 bp with respect to the TSS. (B) Nucleosome occupancy profiles of all genes in the yeast genome from −1000 bp to +1000 bp with respect to the stop codon.

Nucleosome profiles of TATA-containing and TATA-less genes under different nitrogen conditions

Approximately 20% of yeast genes contain a TATA box. Based upon the presence or absence of a TATA box in the promoter region, the genes can be classified as TATA-containing or TATA-less genes25. Plots of the nucleosome profiles on their promoter revealed the differences of nucleosome occupancy on the TATA-containing and TATA-less genes under different nitrogen conditions (Fig. 2). There were no significant differences of the general nucleosome occupancy profiles between TATA-containing and TATA-less genes, and the general genes under different nitrogen conditions (Gln, low; Pro, medium; and Gln-Pro, high.). TATA-containing genes downstream of TSS had a higher nucleosome occupancy score compared to TATA-less genes; however, there was no significant difference upstream of TSS.

Figure 2

Nucleosome occupancy of TATA-containing and TATA-less genes in response to different nitrogen conditions.

Nucleosome occupancy for TATA-containing (973) and TATA-less (4382) promoters, aligned with respect to the TSS. (A) Nucleosome profiles of S. cerevisiae in M.Gln. (B) Nucleosome profiles of S. cerevisiae in M.Pro. (C) Nucleosome profiles of S. cerevisiae for M.Gln shift to M.Pro for 2 h.

Analysis of nucleosome profiles and general gene distribution

The nucleosome occupancy score of NFRs and genes coding regions were calculated to investigate their relationships under different nitrogen conditions. There was no significant relationship (R2 < 0.009) between the nucleosome occupancy amount of gene coding regions and NFRs (Fig. 3). Thus, nucleosome occupancy at NFRs and gene types were examined further. Nucleosome occupancy score in the 10% most strongly positioned nucleosome (High nucleosome occupancy) and the 10% most weakly positioned nucleosomes (low nucleosome occupancy) over the NFRs were used for correlation analysis of nucleosome occupancy at NFRs and the distribution of TATA-containing and TATA-less genes. Figure 4 shows the NFRs of TATA-less genes are frequently covered with high and low nucleosome occupancy. We compared the correlation between normalized nucleosome occupancy scores of gene coding regions and NFRs of TATA-containing and TATA-less genes with high and low nucleosome occupancy (Fig. 5). The nucleosome occupancy score at the coding regions of TATA-containing genes showed no significant change in response to different nitrogen conditions. However, the nucleosome occupancy at coding regions of TATA-less genes while shift from M.Gln to M.Pro was changed significantly with high and low nucleosome occupancy at NFRs compared to M.Gln.

Figure 3

Correlation between nucleosome occupancy at NFRs and gene coding regions.

There was no correlation between nucleosome occupancy at NFRs and gene coding regions. The nucleosome occupancy at NFRs were plotted against gene coding regions under three nitrogen conditions, and no correlation was observed with R2 < 0.009. (A) Correlation between nucleosome occupancy at NFRs and gene coding regions in M.Gln medium. (B) Correlation between nucleosome occupancy at NFRs and gene coding regions in M.Pro medium. (C) Correlation between nucleosome occupancy at NFRs and gene coding regions in response to the shift from glutamine to proline for 2 h.

Figure 4

Percentage of TATA-containing and TATA-less genes with high or low nucleosome occupancy.

Distribution of TATA-containing and TATA-less genes with high or low nucleosome occupancy. Gln, Pro and Gln-Pro represent the culture conditions in M.Gln, M.Pro and in response to the shift from M.Gln to M.Pro for 2 h, respectively.

Figure 5

Correlation between nucleosome occupancy at NFRs and gene coding regions for TATA-containing and TATA-less promoters.

(A) TATA-containing; (B) TATA-less promoters with high nucleosome occupancy; (C) TATA-containing; (D) TATA-less promoters with low nucleosome occupancy. Gln, Pro and Gln-Pro represent the culture conditions in M.Gln, M.Pro and in response to the shift from M.Gln to M.Pro for 2 h, respectively. P-values were calculated using the Wilcoxon rank sum test, P < 0.05 (*), < 0.01 (**) and <0.001 (***).

Analysis of the nucleosome occupancy and gene expression of nitrogen metabolic genes

In order to understand the correlation between nucleosome occupancy and gene expression under different nitrogen conditions, we compared the gene expression of 506 nitrogen metabolic genes that had reported previously15 and the nucleosome occupancy at NFR. As shown in Fig. 6, there are no significant correlation (R2 < 0.0022) between the nucleosome occupancy and nitrogen metabolic gene expression under different nitrogen conditions, although the correlation in M.Gln-Pro medium (Fig. 6C) showed a little higher than that in M.Gln medium (Fig. 6A) and in M.Pro (Fig. 6B).

Figure 6

Correlation between nucleosome occupancy at NFRs and the expression of nitrogen metabolic genes.

Correlation between the nucleosome occupancy at NFRs and the expression of 506 nitrogen metabolic genes were analyzed under different nitrogen conditions. (A) Correlation between the nucleosome occupancy and the expression of nitrogen metabolic genes in M.Gln medium; (B) Correlation between the nucleosome occupancy and the expression of nitrogen metabolic genes in M.Pro medium; (C) Correlation between the nucleosome occupancy and the expression of nitrogen metabolic genes in M.Gln-Pro medium.

Analysis of the nucleosome occupancy and gene expression of NCR-sensitive genes

In order to further examine the nucleosome profile prediction on NCR sensitive genes expression, more than 80 well-characterized NCR-sensitive genes were used to demonstrate their nucleosome occupancy and gene expression in response to different nitrogen conditions. Compared to M.Gln, the remodeling profiles at NFRs that were activated or repressed at least 2.0-fold (Log2ratio > 1 or Log2ratio < −1) and q-value < 1.0E-30 under nitrogen source change were analyzed (Supplementary Table S4). Based upon the nucleosome profiles prediction results, 24 and 26 genes were identified in M.Pro and M.Gln-Pro media, respectively by qRT-PCR. As shown in Fig. 7, when the nucleosome profiles were activated, most of these target genes were up-regulated in M.Pro (Fig. 7A) and M.Gln-Pro (Fig. 7C) media. While the nucleosome profiles were repressed, the gene expression of target genes showed no significant correlation related to nucleosome profiles change under nitrogen condition change, especially in M.Gln-Pro medium (Fig. 7B,D). Furthermore, compared to M.Gln, the nucleosome profiles of the TATA containing genes of NCR sensitive genes at NFRs were preferred to be activated in M.Pro and M.Gln-Pro; the nucleosome profiles of TATA-less genes of NCR sensitive genes at NFRs were preferred to be repressed (Supplementary Table S4). This phenomenon is especially obvious in M.Gln-Pro. When S cerevisiae was shifted from preferred nitrogen conditions to non-preferred nitrogen conditions, most TATA-containing genes of NCR sensitive genes were up-regulated according to nucleosome profile activation, while most TATA-less genes of NCR sensitive were also up-regulated even their nucleosome profiles were repressed.

Figure 7

Comparison of nucleosome profile prediction and NCR-sensitive genes.

The dark grey column represents the strength of gene expression detected by qRT-PCR. Each gene was examined for three times. (A) The expression of nucleosome profiles activated genes in M.Pro medium; (B) The expression of nucleosome profiles repressed genes in M.Pro medium; (C) The expression of nucleosome profiles activated genes in M.Gln-Pro medium; (D) The expression of nucleosome profiles repressed genes in M.Gln-Pro medium.

Discussion

The development of high-throughput sequencing methods has provided global nucleosome positions across the genome of various eukaryotic organisms and environmental conditions626. Much experimental work has focused on determining the relationship between genome-wide nucleosome positioning and the regulation of gene expression and how these changes of nucleosome positioning influence gene expression27. In this study, dynamic remodeling of nucleosomes in the S. cerevisiae S288C genome was mapped under different nitrogen conditions using a combination of micrococcal nuclease digestion, mononucleosome DNA isolation and the Illumina high-throughput sequencing technologies. Analysis of the response of nucleosome profiles to different nitrogen conditions showed the nucleosome profiles associated with different gene types and the expression of NCR-sensitive genes, especially for those TATA containing genes in NCR sensitive genes.

Dynamic remodeling of individual nucleosomes across the yeast genome changed in response to transcriptional perturbation27. The well-positioned nucleosome profiles around the TSS and 3′ end stop codon site were changed under different nitrogen conditions (Fig. 1). However, under each condition, TATA-containing genes showed higher nucleosome occupancy in +1 nucleosome regions compared to TATA-less genes (Fig. 2). Since TATA-containing genes associated with response to stress are highly regulated by nucleosomal and TATA-binding protein-targeted mechanisms28, it appears the TATA-containing genes were also affected easily under these nitrogen conditions. Meanwhile, compared to M.Gln, the nucleosome profile of TATA-containing genes was also easily changed while shifting from glutamine to proline for 2 h (Fig. 2A,C). It means that TATA-containing genes are likely to undergo nucleosome profiles changes29. By comparing the nucleosome occupancy score at gene coding regions of TATA-containing and TATA-less genes with high or low nucleosome occupancy, TATA-less genes displayed a high average nucleosome occupancy in coding regions in M.Gln-Pro medium (Fig. 5). This phenomenon is also similar with nitrogen depletion in fission yeast30.

Glutamine and proline were preferred and non-preferred nitrogen sources, respectively, in response to NCR regulation in S. cerevisiae S288C13. Upon comparing glutamine and proline, shifting from a preferred nitrogen source (glutamine) to a non-preferred nitrogen source (proline), the expression of a set of genes involving amino acid and nucleotide metabolism, protein biosynthesis and degradation, stress response, molecule transport and unknown biological processes are changed14. The NCR-sensitive genes are involved in a set of regulatory mechanisms, including general amino acid control (GAAC), and the unfolded protein response as well as the SPS (Ssy1-Ptr3-Ssy5)-sensing system1215. As shown in Fig. 6, the NCR-sensitive genes GAP1 and CAN1 were regulated by both NCR and GAAC, whereas MEP2, MEP3, DAL5, UGA4, PUT4 and DIP5 were regulated by NCR only12. The different regulatory mechanism might influence the correlation between gene expression and nucleosome prediction.

A shift from a preferred to a non-preferred nitrogen source results in a number of transcriptional changes9. A correlation between nucleosome profiles prediction and expression of NCR-sensitive genes was shown by comparing the NCR-sensitive gene expression with nucleosome occupancy changes and qRT-PCR (Fig. 7). Most of the nucleosome occupancy induced genes were associated with their qRT-PCR result, while the nucleosome occupancy repressed genes showed no correlation. Interestingly, shifting from glutamine to proline, most of the nucleosome occupancy induced NCR sensitive genes were TATA-containing genes, while most of the repressed genes were TATA-less genes (Supplement Table 4), Traditionally, a genome-wide remodeling of nucleosome profiles is related to two main types of regulation, ATP-dependent slip or eviction and post-modification of histone tails31. It is shown recently that the precise replacement of individual nucleosomes at promoters mechanistically regulates transcription also by modulating access of trans-acting factors to specific sites27. Hence, the expression of genes influenced by both nucleosome profile changes of transcriptional factor regulation is more easily distinguish from nucleosome prediction. In non-preferred nitrogen conditions, such as proline or shifting from glutamine to proline, two GATA transcriptional activators, Gln3p and Gat1p activate the NCR-sensitive genes for the utilization of non-preferred nitrogen source and cell growth912. Since the NCR-related proximal promoter motif, GAT1/GLN3/DAL80 occurs preferentially in TATA-less promoters and show a strong preference for transcriptional factor binding32, the GATA factor binding might disturb the nucleosome prediction of the nucleosome occupancy repressed genes or the TATA less genes in NCR sensitive genes.

Earlier studies of transcriptomics and proteomics analysis of NCR metabolism have been reported141533. This work presented here described the response of genome-wide dynamic nucleosome remodeling events to different nitrogen conditions. A correlation between NCR-sensitive genes and nucleosome profiles in response to different nitrogen condition changes was achieved. The expression of NCR-sensitive genes involving transcription regulation showed greater correlation with nucleosome prediction of nucleosome occupancy induced genes or the TATA containing genes in NCR sensitive genes. Nucleosome information involving histone modification and various types of gene regulation might result in the global regulation of gene expression2731. Furthermore, it raises the possibility the gene patterns classified by nucleosome profiles and the assignment of transcriptional regulation factors binding at promoters might be another sketch of the overall transcription model.

Materials and Methods

Strains, media and culture conditions

S. cerevisiae strain S288C was used in this study34. Yeast was grown on yeast nitrogen base (YNB) medium without ammonia, without amino acids and with 3% (w/v) glucose, supplemented with 1 g/L glutamine (M.Gln) or 1 g/L proline (M.Pro) as the sole nitrogen source. The shift from M.Gln to M.Pro (M.Gln-Pro) was achieved by recovering then filtering the cells grown on M.Gln and cultivating them in fresh M.Pro for 2 h14. YPD medium: 10 g/L yeast extract, 20 g/L peptone and 20 g/L glucose.

Overview of the major data

We used a MNase-seq method to isolate mononucleosome-associated DNA from S. cerevisiae S288C grown in media with different sole sources of nitrogen (M.Gln, M.Pro and M.Gln-Pro) and then the ends of the DNA fragments were sequenced. After aligning sequence reads to the reference genome, the genome positions of individual nucleosomes in the genome were estimated. The precise positions calculated by analysis corresponded to the known positions of these nucleosomes and NUCwave software was used for the mapped sequence reads3536. Nucleosomes overlapping the TSS of gene-coding regions were defined as +1 nucleosomes and a nucleosome immediately before a +1 nucleosome was defined as a −1 nucleosome. The space between them was defined as a nucleosome-free region (NFR).

Isolation of nucleosomal DNA

Yeast cultures (400 mL of strain S288C) were grown in YNB medium in 2 L flasks to optical density at 600 nm (OD600) of 0.6–0.8. Yeast cells were treated with 2% (v/v) formaldehyde for 30 min (with vigorous shaking every 10 min) and stopped with 125 mM glycine. Cells were collected, washed twice with sterile water and then suspended in 40 mL of buffer Z (1 M sorbitol, 50 mM Tris (pH 7.4), 10 mM β-mercaptoethanol) with Zymolyase®-20T (SEIKAGAKU BIOBUSINESS CORPORATION, Tokyo, Japan) at 0.25 mg/mL. Cells were spheroplasted at 28–30 °C for about 40 min, depending on how much time was necessary to remove the cell wall of each species completely. Spheroplasts were pelleted by centrifugation (3000g for 10 min) then suspended in 3.5 mL of NP buffer (1 M sorbitol, 50 mM NaCl, 10 mM Tris (pH 7.4), 5 mM MgCl2, 1 mM CaCl2 and 0.075% nonylphenoxypolyethoxylethanol (NP-40), with 1 mM β-mercaptoethanol and 500 μM spermidine freshly added). Samples (600 μL) of spheroplasts were added to six Eppendorf tubes containing a range of concentrations (typically 5–250 units per tube) of micrococcal nuclease (MNase) and incubated at 37 °C for 20 min. Reactions were stopped by the addition of 150 μL of stop buffer (5% (w/v) SDS, 50 mM EDTA), 5 μL of proteinase K (20 mg/mL) was added and reactions were incubated at 65 °C overnight. Samples were extracted with phenol/chloroform/isoamyl alcohol (25:24:1 by vol.), precipitated, with ethanol and treated with DNase-free RNase. Samples were examined by gel electrophoresis with ~90% mononucleosomal DNA, and mononucleosomal DNA was further purified by gel isolation and used for high-throughput sequencing with Illumina HiSeq 2500 (Illumina, San Diego, CA, USA)27.

MNase-seq data and analysis

MNase-seq data were aligned to S288C Saccharomyces cerevisiae Genome Assembly. The raw MNase-seq data and processed data files were deposited at the NCBI Gene Expression Omnibus (GEO) (http://www.ncbi.nlm.nih.gov/geo) under accession numbers GSE75705.

The genome wide nucleosome occupancy maps from sequence reads under different nitrogen conditions were generated, annotated and normalized by using a bioinformatics tool, NUCwave3536. The complete datasets could be found at the GEO with accession numbers GSE75705. The occupancy of a nucleosome was calculated as the logarithm of the estimated number of MNase-seq reads at the peak of a predicted nucleosome position. Furthermore, the average nucleosome occupancy score of each base pair over a range from −1000 bp to +1000 bp relative to TSS and the stop codon site were measured (Supplementary Table S1). The TATA-containing genes and TATA-less genes were classified according to previous description28. The average nucleosome occupancy score of TATA-containing genes and TATA-less genes at each base pair over a range from −1000 bp to +1000 bp relative to TSS were also measured (Supplementary Table S2). The nucleosome free region (NFR) of a gene was defined as the distance between the 5′ -coordinate of the +1 and the 3′ -coordinate of the −1 nucleosome. To further investigate nucleosome occupancy at each gene, their nucleosome occupancy levels were measured by calculating sequencing reads in NFR regions and gene coding regions (Supplementary Table S3). Gln_NFR, Pro_NFR and Gln-Pro are defined as the normalized reads at NFR of each gene; Gln_gene, Pro_gene and Gln-Pro_gene are defined as the normalized reads at gene coding regions of each gene. By comparing the normalized reads in the NFRs of each gene under different nitrogen conditions, we defined the values of Gln_NFR/Pro_NFR or Gln_NFR/Gln-Pro_NFR as the activation (>1) and repression (<1). A scaled difference chi-square test statistic was used to analysis the p-value between Gln_NFR/ Pro_NFR and Gln_NFR/Gln-Pro_NFR. There is an option to generate the statistically significant genes, which calculates an adjusted p-value using the Benjamini-Hochberg method, the adjusted p-value is q-value. While for TATA genes, TATA 0 represent TATA less genes; TATA 1 represent TATA containing genes.

RNA preparation and DNA synthesis

S. cerevisiae S288C was grown overnight at 30 °C in YPD medium. The cells were recovered, washed and diluted to OD600 = 0.1 then transferred to and cultured in M.Gln until OD600 = 0.6–0.8. Cells were pelleted by centrifugation, washed and transferred to M.Pro for 2 h, pelleted again, washed and stored at −80 °C. The procedures for RNA extraction and cDNA synthesis are described fully elsewhere37.

Quantitative real-time PCR (qRT-PCR)

Primers used for qRT-PCR were designed by Beacon Designer 7.0 software (Table 1). The qRT-PCR experiments used a SYBR® Premix Ex TaqTM kit (TaKaRa, Dalian, China). The PCR protocol was: incubation at 95 °C for 30 s then 40 cycles of amplification at 95 °C for 5 s, 55 °C for 20 s and, finally, cooling at 50 °C for 30 s. Reactions (in triplicate) took place in a LightCycler 480 II Real-time PCR instrument (Roche Applied Science, Mannheim, Germany) and mean values were used for further calculations. The fold change was determined by the 2–ΔΔ method normalized to the ACT1 gene38.

Table 1

Oligonucleotides used for qRT-PCR.

Gene	Name	Sequence (5′-3′)
YGK3	YGK3-F	TGTACGTCAGGGAAGGGAAA
	YGK3-R	GCATAAGGTCCTAGCCATTCAAT
DUR1,2	DUR1,2-F	ATGTATCTGGTGGTTCCT
	DUR1,2-R	GCCAATCAGGTTGTTCAA
AVT4	AVT4-F	GCTGGAGGCGACATAACG
	AVT4-R	AACCACTTGCTGATGACAATTTAG
ARG1	ARG1-F	TTCTTACGAGGCAGGTATCTTG
	ARG1-R	AAGTCTTGTTGTCGGTGTAGG
CAR1	CAR1-F	ACGGAATTAGAGCCCTCAATG
	CAR1-R	GGAATCTGTTCGCCTGGAC
DAL1	DAL1-F	AGACGGAGCCACCTACTT
	DAL1-R	CACTACTGAGCCTATAACACCTTC
GAP1	GAP1-F	GAAGCACCACTTGAAGAATAG
	GAP1-R	CCAGAGCCATAACCATAGC
CAN1	CAN1-F	ATGGAGGATGGCATAGGTGATGA
	CAN1-R	GCGTTGGTCAGAGGTGTGGATA
BAT2	BAT2-F	TCGTCCAGATATGAATATG
	BAT2-R	CAGGAACTAAGCACTTAT
YJR011C	YJR011C-F	CGCTTATTATTGGATAGTTGTGAGT
	YJR011C-R	CGACATTTAATCTGTTTAATTCCTGG
PEP4	PEP4-F	CCATACGATTACACGCTTGAAGT
	PEP4-R	CGCAAGAAGGCATCACCAA
UGA4	UGA4-F	ATCTGTGATGGGTGGTGG
	UGA4-R	AGCGGTAGGAATGGAACTTG
TMC1	TMC1-F	CCACGCTAAGATAAAGAAATCAGA
	TMC1-R	CCGCACTTGAACAGGTATC
RTG3	RTG3-F	GAGGAAGCAGAATCGCAATCG
	RTG3-R	ACAGGTGAGTCGTAAGCCATAG
PRB1	PRB1-F	AAGGTGGAGGAGAAGAAGATG
	PRB1-R	CAGTGGACACGAGAGGAG
APG14	APG14-F	GGTCCGTCTCGCTATATTGG
	APG14-R	ATGAGGTCCTGTGACTGTTG
DAL80	DAL80-F	GTGCTTAGTGATTCGTTGA
	DAL80-R	ACCGTGTTCATCTCTTCT
MIG2	MIG2-F	CAGGCAGTAAATGGAGGTA
	MIG2-R	AATTGGCATAGGAGAAGTTG
MOH1	MOH1-F	TGAGCTTCATCACCTACGGTTGTA
	MOH1-R	GTCGCCAGTCAACATTCTTCGT
ECM37	ECM37-F	GCCACCTGATTCCTCGTCAA
	ECM37-R	GCGAATGCCCTGTCTTACTCA
YLR257W	YLR257W-F	CTCCTTCTCTGCGTGCTTCCA
	YLR257W-R	CGTCGTCTGCCTCACCTGTAAT
GDH2	GDH2-F	CCATGCTCGTGTGAGGAACT
	GDH2-R	CGAACACCCAGAGAAGTCATACC
MEP3	MEP3-F	GGTGGTTGGTTGACGCATA
	MEP3-R	TCGGCTTCCTCAGTGACT
URK1	URK1-F	AGGTATAGGTGGTGCTTCAGGTT
	URK1-R	GGCTCTGGCTCTGTCCTCTG
GLN1	GLN1-F	CGCCAAGGAAGGTTACGGTTA
	GLN1-R	CTTCGTCATGTCAGCATTGTCAA
DIP5	DIP5-F	TCATCCACCTCTACATCACCTTCA
	DIP5-R	GAACCACCAATGGCGATCATAGAA
PUT4	PUT4-F	CTGAGATTCCGCAAGGCTATT
	PUT4-R	CCCGTTGGTAATGGTGATGAT
YHM1	YHM1-F	TGCCTCATACCGATAAGAAACAAT
	YHM1-R	TTGGTGTGGTTGGACATCAG
GZF3	GZF3-F	GCGAAGAAGAAGGGCATATTGG
	GZF3-R	GCTCCTGCTTGCGATTCAC
VID30	VID30-F	TGACGACGACGACGATGA
	VID30-R	TTCCTTCTTCTCCACCTTCCT
DAL5	DAL5-F	GGAAGACAGAGGAACAATTCACAA
	DAL5-R	GGTTCTCCAGCCTTTAATAGCATA
DAL4	DAL4-F	GGAGACCACTTACACCAGAG
	DAL4-R	CCGTCATATCAGCACCACAT
DAL7	DAL7-F	GGAGTACAGGTTGAAGGATGACAT
	DAL7-R	GAAGATGGAAGCACTAATCGGTTC
MEP2	MEP2-F	AAGTGGACTACAGTTGGT
	MEP2-R	AACAGCAAGGTTACATCC
CPS1	CPS1-F	TTATTATCATCGGCACCTTCTTC
	CPS1-R	GGACTTAATGGTTCAATCTTCTCA
NIT1	NIT1-F	TGAGAAGGAGATCAAGGAATC
	NIT1-R	TTCGGCAAGATACTTAGCATA
GAT1	GAT1-F	GACGATGACGATGATGATGAC
	GAT1-R	ACCAGAAGCATTATGTGAAGC
DSD1	DSD1-F	GAATGGTATGTAGCAAGGGTCTCT
	DSD1-R	ACAAGCGTGTTGAGGAAGGA
ALD4	ALD4-F	GAAGGTAGAGAGGACGATGT
	ALD4-R	ACCGTTCCAAGACCCATTA
PUT2	PUT2-F	AGAGTCAGCGGATGGTACTTGGAA
	PUT2-R	TAATTTGAAAGGGCAGCGGTTTGT
UGA1	UGA1-F	GCCGTACTATCGTTCAAGAGA
	UGA1-R	GCCAAGCAATGGTCATCCT
ARG4	ARG4-F	CGTCTCTTCCGTATGATTA
	ARG4-R	CCAATTCTGTCTCCGTTA
ACT1	ACT1-F	GGTTGCTGCTTTGGTTATTG
	ACT1-R	CCTTGGTGTCTTGGTCTAC

Additional Information

How to cite this article: Zhang, P. et al. Genome-wide mapping of nucleosome positions in Saccharomyces cerevisiae in response to different nitrogen conditions. Sci. Rep. 6, 33970; doi: 10.1038/srep33970 (2016).