The Arabidopsis SWI/SNF protein BAF60 mediates seedling growth control by modulating DNA accessibility
- Teddy Jégu1, 5, 6,
- Alaguraj Veluchamy2,
- Juan S. Ramirez-Prado1, 2,
- Charley Rizzi-Paillet1,
- Magalie Perez1,
- Anaïs Lhomme1,
- David Latrasse1,
- Emeline Coleno1,
- Serge Vicaire3,
- Stéphanie Legras3,
- Bernard Jost3,
- Martin Rougée4,
- Fredy Barneche4,
- Catherine Bergounioux1,
- Martin Crespi1,
- Magdy M. Mahfouz2,
- Heribert Hirt1, 2,
- Cécile Raynaud1 and
- Moussa Benhamed1, 2Email author
© The Author(s). 2017
Received: 9 January 2017
Accepted: 26 May 2017
Published: 15 June 2017
Plant adaptive responses to changing environments involve complex molecular interplays between intrinsic and external signals. Whilst much is known on the signaling components mediating diurnal, light, and temperature controls on plant development, their influence on chromatin-based transcriptional controls remains poorly explored.
In this study we show that a SWI/SNF chromatin remodeler subunit, BAF60, represses seedling growth by modulating DNA accessibility of hypocotyl cell size regulatory genes. BAF60 binds nucleosome-free regions of multiple G box-containing genes, opposing in cis the promoting effect of the photomorphogenic and thermomorphogenic regulator Phytochrome Interacting Factor 4 (PIF4) on hypocotyl elongation. Furthermore, BAF60 expression level is regulated in response to light and daily rhythms.
These results unveil a short path between a chromatin remodeler and a signaling component to fine-tune plant morphogenesis in response to environmental conditions.
KeywordsSWI/SNF Chromatin Morphogenesis Phytochrome Interacting Factor G-box DNA accessibility PIF4
Plants take advantage of a complex series of sensory receptors to perceive and respond to environmental cues such as the intensity, wavelength, direction, and duration of ambient light signals. For example, the Pr/Pfr reversible forms of red (R) and far-red (FR) sensing phytochrome photoreceptors allow plants to discriminate between mid-day and dusk periods, or shade and plain light conditions, notably because of rapid changes in the R/FR ratio. UV-A/blue light cryptochrome photoreceptors also contribute to control plant adaptive responses such as photomorphogenesis, the developmental transition that seedlings undergo on their first exposure to light. Light signal transduction involves multiple central integrators and sequence-specific effectors [1, 2]. Among them, a subset of the basic helix-loop-helix (bHLH) transcription factors called phytochrome-interacting factors (PIFs) physically interact with both phytochrome and cryptochrome photoreceptors under different light conditions [3, 4].
PIF proteins (PIF1, PIF3, PIF4, PIF5, PIF6, and PIF7 in Arabidopsis) accumulate in the dark to promote etiolated development (skotomorphogenesis). Distinct PIF family members promote hypocotyl elongation by inducing the expression of cell size regulatory genes, notably via the modulation of auxin accumulation in cotyledons . Accordingly, pifq mutant plants, lacking a functional PIF1, PIF3, PIF4, and PIF5 quartet, display a short hypocotyl phenotype . Expression, stability, and activity of the PIFs are also tightly controlled by light and temperature during seedling and adult developmental stages via a combination of transcriptional and post-translational controls [6–9]. Among them, PIF4 represents a central hub for modulating various facets of plant morphogenesis, through the integration of multiple environmental signals, such as photomorphogenesis, the shade avoidance response (SAR), and thermomorphogenesis [10–12]. The expression and activity of PIF4 are also regulated by the circadian clock, possibly influencing time-dependent transcriptional responses to environmental cues .
The bZIP ELONGATED HYPOCOTYL 5 (HY5) transcription factor is another major regulator of thermo- and photomorphogenesis. HY5 regulates hypocotyl elongation in response to light conditions , suggesting that it acts downstream of phytochrome A (phyA), phyB, cryptochrome, and UVR8 photoreceptors [14–16]. The protein stability of PIF4 and HY5 is also regulated by the COP1-DET1 signal integrators acting downstream of these photoreceptors [17, 18]. They control a wide repertoire of target genes, largely by binding to consensus sequences such as G-boxes [14, 19]. In fact, it was recently proved that HY5 inhibits PIF4-directed hypocotyl elongation by competitive chromatin binding to common targets. Such result reveals two independent pathways of PIF4 regulation, one involving DET1/COP1 and other HY5 .
In addition to such sequence-specific factors, accumulating evidence indicates that chromatin modifications, which are crucial components of transcription regulation, participate in light-mediated gene expression [21, 22]. Since the first piece of evidence that increased acetylation of histone H3 and H4 of the Pisum sativum PetE gene promoter correlates with light-induced transcription , several functional studies have shown the influence of histone acetyltransferase and deacetylase activities that oppositely balance histone acetylation levels on Arabidopsis photomorphogenesis [24, 25]. Interestingly, in darkness the histone deacetylase HDA15 can associate with PIF3 and repress expression of photomorphogenic genes, providing a molecular link between a chromatin-modifying activity and light-responsive elements. Accordingly, PIF3 does not contribute to the HDA15 influence on hypocotyl elongation upon light exposure . Profiling of light-induced chromatin state changes along the Arabidopsis genome has further unraveled the wide extent of a chromatin-based program associated with light perception and photomorphogenesis [27–30], but the molecular links triggering light- and sequence-specific chromatin dynamics remain poorly understood.
Chromatin-based transcriptional regulation relies on the deposition and removal of multiple histone and DNA modifications, but also on ATP-dependent chromatin remodeling complexes (CRCs) that modulate histone–DNA contacts . These multi-protein complexes regulate access of genomic regions to transcription factors and to the transcriptional machinery by influencing the structure, dynamic incorporation/eviction of histone variants, or nucleosome positioning. A distinctive feature of CRCs is the presence of a central ATPase domain belonging to the SWI2/SNF2 family [32, 33]. In Arabidopsis, PICKLE, a member of the ATPases CHD3 subfamily (also referred to as PKL or GYMNOS/Suppressor of slr2 [SSL2]/CYTOKININ-HYPERSENSITIVE2)  regulates multiple plant developmental processes, including embryonic development, seed germination, root meristem activity, and photomorphogenesis [35–39]. PKL was found to physically interact with HY5 and to influence hypocotyl elongation by triggering an increase of histone H3 lysine 27 trimethylation (H3K27me3) of several cell elongation-related genes in response to light , a chromatin signature of Polycomb-repressed transcriptional states . However, how the respective activities of chromatin remodeling and light signaling integrate and contribute to control plant adaptive responses remains to be explored on a genome-wide scale.
Based on sequence similarity with metazoan DNA sequences, the Arabidopsis genome encompasses multiple CRCs . These include four SWI2/SNF2 ATPases (BRM, SPLAYED [SYD], MINU1/CHR12, and MINU2/CHR23), four SWI3 proteins (SWI3A to SWI3D), two ACTIN RELATED PROTEINS predicted to belong to SWI/SNF complexes (ARP4 and ARP7), a single protein termed BUSHY (BSH), and two SWI/SNF ASSOCIATED PROTEINS 73 (SWP73A/CHC2 and SWP73B/CHC1), also called BAF60 [42–47]. In addition to their primary activity modulating histone–DNA interactions, and the accessibility of sequence-specific binding proteins onto DNA , CRCs also impinge on various other chromatin-related processes such as histone acetylation, methylation, phosphorylation, and ubiquitination, altogether influencing local chromatin status . These complexes may also impact on the formation of local chromatin loops to modulate gene expression [47, 50–52]. Here, we report that the SWI/SNF chromatin-remodeling factor BAF60 negatively regulates hypocotyl elongation under light and high temperature conditions, and locally represses the expression of multiple cell size regulatory genes by influencing DNA accessibility. Hundreds of BAF60 target genes are shared with PIF4, raising the possibility that BAF60-containing CRCs represents a key factor governing the integration of environmental cues onto the epigenomic landscape.
BAF60 represses hypocotyl elongation
Considering that hypocotyl growth is mostly driven by cell elongation [55, 56] and that hypocotyl cell size strongly correlates with endoreduplication levels, we quantified by flow cytometry the nuclear DNA content of hypocotyls dissected from wild-type and RNAi seedlings. In agreement with their long hypocotyl phenotype, the two BAF60 RNAi lines displayed higher ploidy levels (Fig. 1c; Additional file 1: Figure S1g) and larger hypocotyl cells (Fig. 1d) compared to wild-type seedlings. Taken together, these observations show that BAF60 can repress hypocotyl elongation under light and high temperature conditions, possibly through the control of both cell elongation and endoreduplication.
Regulation of BAF60 gene expression by light and the circadian clock
To test the potential direct effect of light on BAF60 transcription, the BAF60 promoter was cloned upstream of a green fluorescent protein (GFP) coding sequence and stably introduced into wild-type plants (BAF60::WPP-GFP-BLRP). GFP transcript levels were quantified upon transferring light-grown plants to darkness, showing a significant decrease after 3 h (Fig. 2g), similar to the phenomenon observed when measuring the expression levels of the BAF60 locus through RT-qPCR. Comparison of GFP RNA levels between dark- and light-grown seedlings also showed more BAF60 promoter activity in the light condition (Additional file 1: Figure S2d). Finally, analysis of public genome-wide data allowed us to determine that the BAF60 promoter domain becomes more susceptible to DNase I footprinting during de-etiolation  (Additional file 1: Figure S2e), corroborating its enhanced transcriptional activity in response to light. Taken together, these results indicate that BAF60 expression is regulated at the transcriptional level, and possibly also at the post-transcriptional level, by light and the circadian clock.
BAF60 binds nucleosome-free regions of expressed genes
To understand the relationship between BAF60 and gene expression, we tested how its enrichment relates to mRNA levels. This revealed a positive correlation between the binding frequency of this protein and the number of transcripts of specific loci (Fig. 3e). We also analyzed the relationship between gene expression and DNA accessibility in a genome-wide fashion, finding that the most expressed genes display a highest level of accessibility in their 5′ and 3′ ends compared to less transcribed genes (Additional file 1: Figure S4). Finally, to determine whether BAF60 associates with specific chromatin contexts, we compared its genome-wide positioning with publically available epigenomic profiles of histone modifications. This showed that more than 50% of the BAF60 target genes are marked by histone 3 lysine 9 acetylation (H3K9ac), confirming that BAF60 frequently binds active chromatin regions (Fig. 3f). On the other hand, BAF60 binding co-localizes, but to a much lower extent, with the repressive marks H3K27me3 and H3K9me2, not only in the gene body but also in the promoter regions, upstream of the TSS (Additional file 1: Figure S5).
BAF60 targets G-box motifs and acts antagonistically to PIF4
To understand the impact of BAF60 on gene expression we performed a transcriptomic analysis of BAF60 RNAi lines and found that 1103 of the BAF60 targets displayed missregulation in both of the BAF60 RNAi lines, 404 downregulated and 699 upregulated (Fig. 4f). This result suggests that BAF60 regulates either positively or negatively the expression of several of its targets. A gene ontology analysis of the BAF60 target genes that are upregulated in both RNAi lines revealed that there is a significant enrichment of genes involved in the response to light stimuli and hormone-mediated signaling (Additional file 1: Figure S8).
Given the prominent role of PIF4 and the negative effect of BAF60 on hypocotyl elongation, we analyzed in detail the functional relationship of BAF60 with cell size regulatory genes that impact this process. These were found to be targeted by both PIF4 and BAF60 and upregulated in the BAF60 RNAi line and downregulated in the pifq mutant (IAA19, ST2a, XTR7, SDR, HFR1, and BEE1). Our ChIP-seq and ATAC-seq analyses indicated that both BAF60 and PIF4 target the promoter and/or the gene body of these genes around a G-box motif at loci with high DNA accessibility (Fig. 5c). Targeted ChIP-qPCR further confirmed that BAF60-CFP is enriched over these loci (Additional file 1: Figure S9). To confirm its influence on their expression, we quantified the transcripts of these genes upon knocking down BAF60. RT-qPCR analysis showed highly increased mRNA levels in the two BAF60 RNAi lines compared to wild-type seedlings for the six genes (Fig. 5d). Consistently, expression of these genes was reduced in overexpressing (OE) lines (Additional file 1: Figure S10a). Altogether, these observations indicate that BAF60 represses these genes in cis.
Several chromatin modifications such as histone methylation and acetylation have recently been shown to influence genome expression reprogramming during photo- and thermomorphogenesis [12, 21, 22]. Still, the role of chromatin remodeling in the environmental control of seedling morphogenesis is poorly understood. In this study, we found that the SWI/SNF CRC subunit BAF60 is required for efficient repression of hypocotyl elongation under light and high temperature conditions by repressing cell elongation and endoreduplication, acting on local chromatin status of cell size regulatory genes. Hypocotyl elongation relies only on cell growth and not on cell division; in hypocotyl cells, DNA replication leads to 4C and 8C ploidy levels in light- and dark-grown seedlings, respectively , a control of DNA content that is thought to be associated with cell size regulation [55, 56]. Here, we first observed that inhibition of endoreduplication by light in the hypocotyl requires BAF60. This adds to our former identification of BAF60 as a protein partner of KRP5, a cell cycle inhibitor that promotes endoreduplication and hypocotyl elongation . The functional significance of this interaction is still elusive, but KRP5 could hypothetically interact with BAF60 to reduce its activity under environmental conditions that enhance hypocotyl elongation, such as darkness and high temperature.
This work also unveiled that the expression of BAF60 itself is regulated at the transcriptional level. Greater accessibility of the BAF60 promoter during de-etiolation detected in DNase I footprinting datasets  is a good in vivo indication of enhanced transcriptional activity in response to light signals. Consistently, we confirmed that its promoter domain drives more transcription under light than under dark conditions using a pBAF60::WPP-GFP-BLRP transcriptional reporter.
Interestingly, this combination of regulatory events tightly control BAF60 protein level in response to both light and circadian signals. Circadian rhythms are set by diurnal fluctuations of light and temperature signals [61, 62]. In addition, clock-controlled transcriptional regulation is highly integrated with environmental signals such as light and temperature, notably converging onto a common repertoire of genes that usually exert their function in specific light conditions or daytime . Clock- or light-controlled changes in BAF60 abundance might be directly linked to its function. Accordingly, BAF60 downregulation strongly enhances hypocotyl elongation under photocycles. Given that BAF60 triggers the repression of IAA19, ST2a, XTR7, SDR, HFR1, and BEE1 genes, low BAF60 expression levels at night, or upon a shift to darkness, could therefore mediate direct control of these cell size regulatory genes when plants face these conditions.
Genome-wide profiling of nucleosome positioning has recently been performed in a number of eukaryotic organisms, showing that gene bodies have high nucleosome occupancy [64, 65]. In contrast, regulatory regions located in promoters or terminators tend to have low nucleosome occupancy and often contain a NFR. Recently, the INO80 and SWR1 chromatin remodeling complexes have been found to target NFR regions to regulate histone variant deposition . Using a combination of genome-wide methodologies (ChIP-seq, ATAC-seq, and MNase-seq), we found that BAF60 binds hundreds of NFRs, showing that, like in the yeast Saccharomyces cerevisiae, plant chromatin remodelers can target nucleosome-free regions. The position of the nucleosomes along a particular DNA domain and the presence of a NFR can have profound effects on its accessibility and on its transcriptional regulation by trans-acting factors . Nevertheless, these results do not imply that BAF60 acts on its targets as a transcriptional activator, but rather suggest that BAF60 regulates the expression of genes with a high expression level. As it was recently shown for the catalytic subunit BRM, BAF60 can function as both an activator or inhibitor of gene expression : genes that were upregulated in BAF60 RNAi lines and genes that were downregulated were both significantly enriched amongst BAF60 targets. In the case of cell size regulatory genes, our analyses by FAIRE-qPCR showed that DNA accessibility around NFRs was increased in BAF60 RNAi lines, which at the same time is correlated with increased gene expression. These findings are highly consistent with the local association of BAF60 protein with these domains in vivo and with its role as a negative regulator of their expression.
Other proteins might facilitate the recruitment of BAF60 to such NFRs. Indeed, the SWI/SNF chromatin remodeler of Arabidopsis, BRAHMA (BRM), was reported to be recruited by the H3K27 demethylase REF6 to specific sequences. The demethylase was found to bind its target motifs through their recognition by zinc-finger (ZnF) domains, facilitating BRM targeting . A similar mechanism may be occurring in the case of BAF60, and for future studies we will consider determining its DNA-binding protein interactors, which would possibly provide evidence on the existence of other molecules involved in the recruitment of BAF60 to its targets within NFRs.
Determining the distribution of BAF60 over the genome led us to identify its preferential enrichment over G-box sequence motifs. Interestingly, this motif is also recognized by PIF transcription factors . Consistently, we observed that BAF60 and PIF4 share 1477 target genes. Several PIF family members (PIF1, PIF3, PIF4, and PIF5) play central roles in the repression of hypocotyl elongation in darkness . Interestingly, PIF4 is also central for integrating light and high temperature signaling and for mediating their influence on seedling growth . In contrast to PIF4, BAF60 knockdown lines display a long hypocotyl phenotype under light and high temperature. This suggests that BAF60 and PIF4 antagonistically regulate a common set of genes. Furthermore, in our study, BAF60 and PIF4 were found to compete for their targets, a similar mechanism of PIF4 regulation to the one reported by Gangappa and Kumar . PIF4 activity appears to be regulated and antagonized by various factors, including BAF60 and HY5, which compete with PIF4 for binding to several targets, including cell size regulatory genes. Consistently, increased expression of IAA19, ST2a, XTR7, SDR, HFR1, and BEE1 genes in the BAF60 RNAi lines correlated with higher accessibility of their promoter domains. The observation that the 5′ regions of genes became hypersensitive to micrococcal nuclease upon gene activation in Drosophila was among the earliest demonstrations of this phenomenon [70–73]. The appearance of hypersensitive sites reflects a loss or a destabilization of nucleosomes at the promoters of transcribed genes . Several mechanisms are thought to act in concert to achieve this, and can involve SWI/SNF CRCs [74–76]. Hence, BAF60 could bind NFRs of genes involved in hypocotyl growth control to negatively regulate their transcription through the modulation of DNA accessibility.
Plant material and growth conditions
Arabidopsis thaliana seeds of both RNAi lines BAF60_1 (CS30982) and BAF60_2 (CS23961), in the Wassilewskija (Ws) background, were obtained from the Nottingham Arabidopsis Stock Centre (NASC). Plant lines overexpressing BAF60 35S::BAF60-CFP_1 and 35::BAF60-CFP_2 have been described previously . BAF60::WPP-GFP-BLRP, 35S::GFP plants were in the Columbia-0 background (Col0). Plants were grown in chambers at 20 °C on sterile half-strength MS medium and 0.8% agar under long days (16 h of light at 20 °C, 8 h of darkness at 18 °C; LD), short days (8 h of light at 20 °C, 16 h of darkness at 18 °C; SD) or darkness conditions (16 h of darkness at 20 °C, 8 h of darkness at 18 °C; D). Seeds were surface-sterilized by treatment with bayrochlore for 20 min, washed, and imbibed in sterile-water for 2–4 days at 4 °C to obtain homogeneous germination. Phenotyping under high temperature was performed by cultivating plants under continuous light at 23 °C during 3 days and then at 28 °C for 4 days. For hypocotyl growth assays, plates were scanned and hypocotyl length was measured using the Image J software (http://rsb.info.nih.gov/ij/).
Fourteen-day-old hypocotyls were chopped with a razor 587 blade in 1 mL of Galbraith buffer supplemented with 1% polyvinylpyrrolidone 10,000, 5 mM metabisulfite, and 5 mg/mL RNase from a stock solution at 50 units/mg. Propidium iodide was added to the filtered supernatants at a final concentration of 50 μg/ml. Endoreduplication levels of 5000–10,000 stained nuclei were determined using a Cyflow SL flow cytometer (Partec) with a 532 nm solid state laser (30 mW) excitation and emission collected after a 630/30 nm filter.
Cell size measurement
For cell size measurement, hypocotyls of 14-day-old plantlets were fixed in ethanol:acetic acid (3:1) and washed in 70% (v/v) ethanol during 20 min at room temperature. Plantlets were subsequently cleared with chloral hydrate (8 g of chloral hydrate [Sigma], 2 mL of 50% glycerol [w/v], and 1 mL of water) overnight. The day after, samples were mounted on slides in water under coverslips, and differential interference contrast macroscopy (AZ100; Nikon) was used to capture images with a Nikon RI1 video camera. Cell size measurement was performed with the ImageJ software.
RNA extraction and real-time quantitative PCR analysis
Total RNAs were extracted from 180 mg of seedlings with the RNeasy MiniPrep kit (Qiagen) for RT-qPCR or with the ZR Plant RNA MiniPrep kit (Zymo Research) for RNA-seq, according to the manufacturer’s instructions. First strand cDNA was synthesized from 2 μg of total RNAs using Improm-II reverse transcriptase (A3802, Promega) according to the manufacturer’s instructions. We mixed 1/25th of the synthesized cDNA with 500 nM of each primer and LightCycler® 480 Sybr Green I master mix (Roche Applied Science) for quantitative PCR analysis (qPCR). Products were amplified and fluorescent signals acquired with a LightCycler® 480 detection system. The specificity of amplification products was determined by melting curves. Ubq10 was used as an internal control for signal normalization. The Exor4 relative quantification software (Roche Applied Science) automatically calculates relative expression levels of the selected genes with algorithms based on the ΔΔCt method. Data were from duplicates of at least three biological replicates. All primer sequences used are given in Additional file 1: Table S1.
Formaldehyde-assisted isolation of regulatory elements assay
Formaldehyde-assisted isolation of regulatory elements (FAIRE) was performed as described by Ariel et al. . Two grams of 7-day-old seedlings were crosslinked in 1% (v/v) formaldehyde at room temperature for 15 min. Purified nuclei were resuspended in 500 μl, but only 50 μl were used (diluted to 500 μl) to continue with the protocol.
Chromatin immunoprecipitation assay
ChIP-seq assays were performed on 14-day-old BAF60-CFP overexpressing line 1 (OE BAF60-CFP_1), formerly described by Jegu and collaborators , using an anti-GFP antibody (Clontech 632592), and on 14-day-old wild type plant (Col0) using an anti-H3K9ac antibody (Millipore), an anti-H3K9me2 antibody (Abcam), and an anti-H3K27me3 antibody (Millipore). ChIP-qPCR assays were performed on 7-day-old BAF60-CFP overexpressing line 2 (OE BAF60-CFP_2) using an anti-GFP antibody (Clontech 632592) and an anti-IgG antibody (Millipore). The ChIP protocol was modified from Gendrel et al. . Briefly, after plant material fixation in 1% (v/v) formaldehyde, tissues were homogenized and nuclei isolated and lysed. Cross-linked chromatin was sonicated using a water bath Bioruptor UCD-200 (Diagenode, Liège, Belgium; 15 s on/15 s off pulses; 15 times). The complexes were immunoprecipitated with antibodies overnight at 4 °C with gentle shaking and incubated for 1 h at 4 °C with 50 μL of Protein AG UltraLink Resin (Thermo scientific). The beads were washed for 6 × 5 min in ChIP Dilution Buffer (SDS 0.01%, Triton X-100 1.1%, 1.2 mM EDTA, pH 8, 16.7 mM Tris-HCl, pH 8, and 167 mM NaCl) and twice in TE. ChIPed material was eluted by two 15-min incubations at room temperature with 250 μL elution buffer (SDS 1%, 0.1 M NaHCO3). Chromatin was reverse-crosslinked by adding 20 μL of NaCl 5 M and incubated overnight at 65 °C.
For qPCR, reverse-crosslinked DNA was then recovered using the IPure kit (Diagenode, Liège, Belgium) and analyzed by RT-qPCR. An aliquot of untreated sonicated chromatin was processed in parallel and used as the total input DNA control.
For ChIP-seq, reverse-crosslinked DNA was submitted to RNase and proteinase K digestion and extracted with phenol-chloroform. DNA was ethanol precipitated in the presence of 20 μg of glycogen and resuspended in 50 μL of H20 in a siliconized tube. IP or input DNA (10 ng) was used for ChIP-seq library construction using a NEBNext® Ultra DNA Library Prep Kit for Illumina®.
Microccocal nuclease sensitivity assay
Five grams of 14-day-old seedlings were ground and nuclei were isolated with 4 °C buffer (0.25 M sucrose, 10 mM Tris-HCl, 10 mM MgCl2, 1% Triton, 5 mM β-mercaptoethanol) containing proteinase inhibitor cocktail (Roche), filtered with a 63 μm filter, and incubated in MN buffer (20 mM Tris-HCl, 70 mM NaCl, 20 mM KCl, 5 mM MgCl2, 3 mM CaCl2). MNase (100 U; Invitrogen®) was added to initiate the kinetics. Adding EGTA and EDTA to a final concentration of 2 mM stopped the reaction. DNA was extracted with phenol/chloroform, precipitated with isopropanol, and resuspended in ultra-pure water. Purified DNA was run on a 1% agarose gel and the band corresponding to mononucleosome was excised and purified with a MinElute Gel Extraction Kit (Qiagen). Mononucleosome DNA (10 ng) was used for ChIP-seq library construction using a NEBNext® Ultra DNA Library Prep Kit for Illumina®.
Assay for transposase-accessible chromatin with high-throughput sequencing
Fourteen-day-old seedlings (100 mg) were ground and nuclei were isolated with 4 °C buffer (0.25 M sucrose, 10 mM Tris-HCl, 10 mM MgCl2, 1% Triton, 5 mM β-mercaptoethanol) containing proteinase inhibitor cocktail (Roche) and filtered with a 63 μm filter. Nuclei were resuspended in 1× TD buffer (Illumina FC-121-1030) and 2.5 μL of Tn5 transposes (Illumina FC-121-1030) were added. The transposition reaction was performed at 37 °C for 30 min, and DNA was purified using a Qiagen MinElute Kit. DNA libraries were amplified for a total of eight cycles as described by Buenrostro et al. .
Deep sequencing and bioinformatics
Single-end sequencing of ChIP samples was performed using Illumina GAIIx with a read length of 50 bp. Reads were quality controlled using FASTQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/). Trimmomatic was used for quality trimming. The reads were mapped onto the TAIR10 assembly using Bowtie  with mismatch permission of 1 bp. Unique mapping of reads was adopted. To identify regions that were significantly enriched, we used MACS2 . Visualization and analysis of genome-wide enrichment profiles were done with IGB. Peak annotations, such as proximity to genes and overlap on genomic features such as transposons and genes, were assigned using HOMER. NGSplot was used to profile the enrichment of this mark at transcriptional start sites (TSSs) and along the gene . To identify regions that were significantly enriched in histone modification data (H3K27me3), we used SICER  with parameters of W:200 (window length) and G:200 (gap size) for H3K9ac, and W:200, G:600 for H3K27me3 and H3K9me2. Nucleosome occupancy and shift between the wild type and BAF60 mutant were deduced from the MNase sequencing data using DANPOS . DiffReps was used to find the differential marking between two histone modifications . Spatial binding of the two peaks of BAF60 and PIF4 were done by position-wise comparison using a binning approach and plotted in hexplot.
BAF60 (AT5G14170), IAA19 (AT3G15540), ST2A (AT5G07010), XTR7 (AT4G14130), SDR (AT5G02540), HFR1 (AT1G02340), BEE1 (AT1G18400).
This work was funded by King Abdullah University of Science and Technology (KAUST) and Université Paris Sud. The IPS2 benefits from the support of the LabExSaclay Plant Sciences-SPS (ANR-10-LABX-0040-SPS).
Availability of data and materials
The datasets generated and/or analyzed during the current study are available in the Gene Expression Omnibus with accession GSE89346.
TJ performed the wet lab experiments and AV performed the bioinformatics analyses. JSRP, AL, CRP, DL, EC, MP, and MR prepared material. SV, BJ, and SL were in charge of the sequencing. TJ, CR, HH, JSRP, CB, MM, FB, and MB analyzed the data. TJ and MB designed the research. MB and JSRP wrote the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Galvão VC, Fankhauser C. Sensing the light environment in plants: photoreceptors and early signaling steps. Curr Opin Neurobiol. 2015;34:46–53.View ArticlePubMedGoogle Scholar
- Jiao Y, Lau OS, Deng XW. Light-regulated transcriptional networks in higher plants. Nat Rev Genet. 2007;8:217–30.View ArticlePubMedGoogle Scholar
- Ni M, Tepperman JM, Quail PH. PIF3, a phytochrome-interacting factor necessary for normal photoinduced signal transduction, is a novel basic helix-loop-helix protein. Cell. 1998;95:657–67.View ArticlePubMedGoogle Scholar
- Pedmale UV, Huang SSC, Zander M, Cole BJ, Hetzel J, Ljung K, et al. Cryptochromes interact directly with PIFs to control plant growth in limiting blue light. Cell. 2016;164:233–45. doi:10.1016/j.cell.2015.12.018.View ArticlePubMedGoogle Scholar
- Kohnen M, Schmid-Siegert E, Trevisan M, Allenbach Petrolati L, Sénéchal F, Müller-Moulé P, et al. Neighbor detection induces organ-specific transcriptomes, revealing patterns underlying hypocotyl-specific growth. Plant Cell. 2016;28:1–53. doi:10.1105/tpc.16.00463.
- Shin J, Kim K, Kang H, Zulfugarov IS, Bae G, Lee C-H, et al. Phytochromes promote seedling light responses by inhibiting four negatively-acting phytochrome-interacting factors. Proc Natl Acad Sci U S A. 2009;106:7660–5. doi:10.1073/pnas.0812219106.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen M, Chory J. Phytochrome signaling mechanisms and the control of plant development. Trends Cell Biol. 2011;21:664–71. doi:10.1016/j.tcb.2011.07.002.View ArticlePubMedPubMed CentralGoogle Scholar
- Al-Sady B, Ni W, Kircher S, Schäfer E, Quail PH. Photoactivated phytochrome induces rapid PIF3 phosphorylation prior to proteasome-mediated degradation. Mol Cell. 2006;23:439–46.View ArticlePubMedGoogle Scholar
- Leivar P, Tepperman JM, Monte E, Calderon RH, Liu TL, Quail PH. Definition of early transcriptional circuitry involved in light-induced reversal of PIF-imposed repression of photomorphogenesis in young Arabidopsis seedlings. Plant Cell. 2009;21:3535–53. doi:10.1105/tpc.109.070672.View ArticlePubMedPubMed CentralGoogle Scholar
- Wigge PA. Ambient temperature signalling in plants. Curr Opin Plant Biol. 2013;16:661–6. doi:10.1016/j.pbi.2013.08.004.View ArticlePubMedGoogle Scholar
- Leivar P, Monte E. PIFs: systems integrators in plant development. Plant Cell. 2014;26:56–78.View ArticlePubMedPubMed CentralGoogle Scholar
- Quint M, Delker C, Franklin KA, Wigge PA, Halliday KJ, van Zanten M. Molecular and genetic control of plant thermomorphogenesis. Nat Plants. 2016;2:15190. doi:10.1038/nplants.2015.190.View ArticlePubMedGoogle Scholar
- Nozue K, Maloof JN. Diurnal regulation of plant growth. Plant Cell Environ. 2006;29:396–408.View ArticlePubMedGoogle Scholar
- Oyama T, Shimura Y, Okada K. The Arabidopsis HY5gene encodes a bZIP protein that regulates stimulus-induced development of root and hypocotyl. Genes Dev. 1997;11:2983–95.Google Scholar
- Koornneef M, Rolff E, Spruit CJP. Genetic control of light-inhibited hypocotyl elongation in Arabidopsis thaliana (L.) Heynh. Z Pflanzenphysiol. 1980;100:147–60.View ArticleGoogle Scholar
- Ang LH, Chattopadhyay S, Wei N, Oyama T, Okada K, Batschauer a, et al. Molecular interaction between COP1 and HY5 defines a regulatory switch for light control of Arabidopsis development. Mol Cell. 1998;1:213–22.View ArticlePubMedGoogle Scholar
- Osterlund MT, Wei N, Deng XW. The roles of photoreceptor systems and the COP1- targeted destabilization of HY5 in light control of Arabidopsis seedling development. Plant Physiol. 2000;124:1520–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Delker C, Sonntag L, James G, Janitza P, Ibañez C, Ziermann H, et al. The DET1-COP1-HY5 pathway constitutes a multipurpose signaling module regulating plant photomorphogenesis and thermomorphogenesis. Cell Rep. 2014;9:1983–9.View ArticlePubMedGoogle Scholar
- Lee J, He K, Stolc V, Lee H, Figueroa P, Gao Y, et al. Analysis of transcription factor HY5 genomic binding sites revealed its hierarchical role in light regulation of development. Plant Cell Online. 2007;19:731–49. doi:10.1105/tpc.106.047688.View ArticleGoogle Scholar
- Gangappa SN, Kumar SV. DET1 and HY5 Control PIF4-mediated thermosensory elongation growth through distinct mechanisms. Cell Rep. 2017;18:344–51. doi:10.1016/j.celrep.2016.12.046.View ArticlePubMedPubMed CentralGoogle Scholar
- He G, Elling A, Deng XW. The epigenome and plant development. Annu Rev Plant Biol. 2011;62:411–35.View ArticlePubMedGoogle Scholar
- Barneche F, Malapeira J, Mas P. The impact of chromatin dynamics on plant light responses and circadian clock function. J Exp Bot. 2014;65:2895–913.View ArticlePubMedGoogle Scholar
- Chua YL, Watson L, Gray JC. The transcriptional enhancer of the pea plastocyanin gene associates with the nuclear matrix and regulates gene expression through histone acetylation. Plant Cell. 2003;15:1468–79. doi:10.1105/tpc.011825.(MARs).View ArticlePubMedPubMed CentralGoogle Scholar
- Bertrand C, Benhamed M, Li YF, Ayadi M, Lemonnier G, Renou JP, et al. Arabidopsis HAF2 gene encoding TATA-binding protein (TBP)-associated factor TAF1, is required to integrate light signals to regulate gene expression and growth. J Biol Chem. 2005;280:1465–73.View ArticlePubMedGoogle Scholar
- Benhamed M, Bertrand C, Servet C, Zhou D-X. Arabidopsis GCN5, HD1, and TAF1/HAF2 interact to regulate histone acetylation required for light-responsive gene expression. Plant Cell. 2006;18:2893–903.View ArticlePubMedPubMed CentralGoogle Scholar
- Liu M-J, Wu S-HS-H WJ-F, Lin W-D, Wu Y-C, Tsai T-Y, et al. Translational landscape of photomorphogenic Arabidopsis. Plant Cell. 2013;25:3699–710. doi:10.1105/tpc.113.114769.View ArticlePubMedPubMed CentralGoogle Scholar
- Charron J-BF, He H, Elling A, Deng XW. Dynamic landscapes of four histone modifications during deetiolation in Arabidopsis. Plant Cell. 2009;21:3732–48.View ArticlePubMedPubMed CentralGoogle Scholar
- Bourbousse C, Ahmed I, Roudier F, Zabulon G, Blondet E, Balzergue S, et al. Histone H2B monoubiquitination facilitates the rapid modulation of gene expression PLoS Genet. 2012;8:1–16.Google Scholar
- Sullivan AM, Arsovski AA, Lempe J, Bubb KL, Weirauch MT, Sabo PJ, et al. Mapping and dynamics of regulatory DNA and transcription factor networks in A. thaliana. Cell Rep. 2014;8:2015–30. doi:10.1016/j.celrep.2014.08.019.View ArticlePubMedGoogle Scholar
- Rutowicz K, Puzio M, Halibart-Puzio J, Lirski M, Kroteń MA, Kotliński M, et al. A specialized histone H1 variant is required for adaptive responses to complex abiotic stress and related DNA methylation in Arabidopsis. Plant Physiol. 2015;169:2080–101. doi:10.1104/pp.15.00493.PubMedPubMed CentralGoogle Scholar
- Clapier CR, Cairns BR. The biology of chromatin remodeling complexes. Annu Rev Biochem. 2009;78:273–304.View ArticlePubMedGoogle Scholar
- Bork P, Koonin EV. An expanding family of helicases within the ‘DEAD/H’ superfamily. Nucleic Acids Res. 1993;21:1992–3.View ArticleGoogle Scholar
- Eisen JA, Sweder KS, Hanawalt PC. Evolution of the SNF2 family of proteins: subfamilies with distinct sequences and functions. Nucleic Acids Res. 1995;23:2715–23.View ArticlePubMedPubMed CentralGoogle Scholar
- Ho KK, Zhang H, Golden BL, Ogas J. PICKLE is a CHD subfamily II ATP-dependent chromatin remodeling factor. Biochim Biophys Acta Gene Regul Mech. 2013;1829:199–210. doi:10.1016/j.bbagrm.2012.10.011.View ArticleGoogle Scholar
- Ogas J, Kaufmann S, Henderson J, Somerville C. PICKLE is a CHD3 chromatin-remodeling factor that regulates the transition from embryonic to vegetative development in Arabidopsis. Proc Natl Acad Sci U S A. 1999;96:13839–44.View ArticlePubMedPubMed CentralGoogle Scholar
- Fukaki H, Taniguchi N, Tasaka M. PICKLE is required for SOLITARY-ROOT/IAA14-mediated repression of ARF7 and ARF19 activity during Arabidopsis lateral root initiation. Plant J. 2006;48:380–9.View ArticlePubMedGoogle Scholar
- Perruc E, Kinoshita N, Lopez-Molina L. The role of chromatin-remodeling factor PKL in balancing osmotic stress responses during Arabidopsis seed germination. Plant J. 2007;52:927–36.View ArticlePubMedGoogle Scholar
- Aichinger E, Villar CBR, Di Mambro R, Sabatini S, Köhler C. The CHD3 chromatin remodeler PICKLE and polycomb group proteins antagonistically regulate meristem activity in the Arabidopsis root. Plant Cell. 2011;23:1047–60. doi:10.1105/tpc.111.083352.View ArticlePubMedPubMed CentralGoogle Scholar
- Jing Y, Zhang D, Wang X, Tang W, Wang W, Huai J, et al. Arabidopsis chromatin remodeling factor PICKLE interacts with transcription factor HY5 to regulate hypocotyl cell elongation. Plant Cell. 2013;25:242–56. doi:10.1105/tpc.112.105742.View ArticlePubMedPubMed CentralGoogle Scholar
- Derkacheva M, Hennig L. Variations on a theme: Polycomb group proteins in plants. J Exp Bot. 2014;65:2769–84.View ArticlePubMedGoogle Scholar
- Gentry M, Hennig L. Remodelling chromatin to shape development of plants. Exp Cell Res. 2014;321:40–6. doi:10.1016/j.yexcr.2013.11.010.View ArticlePubMedGoogle Scholar
- Meagher RB, Deal RB, Kandasamy MK, McKinney EC. Nuclear actin-related proteins as epigenetic regulators of development. Plant Physiol. 2005;139:1576–85.View ArticlePubMedPubMed CentralGoogle Scholar
- Jerzmanowski A. SWI/SNF chromatin remodeling and linker histones in plants. Biochim Biophys Acta Gene Struct Expr. 2007;1769:330–45.View ArticleGoogle Scholar
- Kwon CS, Wagner D. Unwinding chromatin for development and growth: a few genes at a time. Trends Genet. 2007;23:403–12.View ArticlePubMedGoogle Scholar
- Sang Y, Silva-Ortega CO, Wu S, Yamaguchi N, Wu MF, Pfluger J, et al. Mutations in two non-canonical Arabidopsis SWI2/SNF2 chromatin remodeling ATPases cause embryogenesis and stem cell maintenance defects. Plant J. 2012;72:1000–14.View ArticlePubMedPubMed CentralGoogle Scholar
- Vercruyssen L, Verkest A, Gonzalez N, Heyndrickx KS, Eeckhout D, Han S-K, et al. ANGUSTIFOLIA3 binds to SWI/snf chromatin remodeling complexes to regulate transcription during Arabidopsis leaf development. Plant Cell. 2014;26:210–29. doi:10.1105/tpc.113.115907.View ArticlePubMedPubMed CentralGoogle Scholar
- Jégu T, Latrasse D, Delarue M, Hirt H, Domenichini S, Ariel F, et al. The BAF60 subunit of the SWI/SNF chromatin-remodeling complex directly controls the formation of a gene loop at FLOWERING LOCUS C in Arabidopsis. Plant Cell. 2014;26:538–51. doi:10.1105/tpc.113.114454.View ArticlePubMedPubMed CentralGoogle Scholar
- Ho L, Crabtree GR. Chromatin remodelling during development. Nature. 2010;463:474–84. doi:10.1038/nature08911.View ArticlePubMedPubMed CentralGoogle Scholar
- Hargreaves DC, Crabtree GR. ATP-dependent chromatin remodeling: genetics, genomics and mechanisms. Cell Res. 2011;21:396–420. doi:10.1038/cr.2011.32.View ArticlePubMedPubMed CentralGoogle Scholar
- Euskirchen G, Auerbach RK, Snyder M. SWI/SNF chromatin-remodeling factors: Multiscale analyses and diverse functions. J Biol Chem. 2012;287:30897–905.View ArticlePubMedPubMed CentralGoogle Scholar
- Shi J, Whyte W a, Zepeda-mendoza CJ, Milazzo JP, Shen C, Roe J, et al. enhancer-mediated Myc regulation Role of SWI/SNF in acute leukemia maintenance and enhancer-mediated Myc regulation. Genes Dev. 2013;2648–62.Google Scholar
- Jegu T, Domenichini S, Blein T, Ariel F, Christ A, Kim SK, et al. A SWI/SNF chromatin remodelling protein controls cytokinin production through the regulation of chromatin architecture. PLoS One. 2015;10:1–18.View ArticleGoogle Scholar
- Sacharowski SP, Gratkowska DM, Sarnowska E, Kondrak P, Jancewicz I, Porri A, et al. SWP73 subunits of Arabidopsis SWI/SNF chromatin remodeling complexes play distinct roles in leaf and flower development. Plant Cell. 2015;27:1889–906. doi:10.1105/tpc.15.00233.View ArticlePubMedPubMed CentralGoogle Scholar
- Crane YM, Gelvin SB. RNAi-mediated gene silencing reveals involvement of Arabidopsis chromatin-related genes in Agrobacterium-mediated root transformation. Proc Natl Acad Sci U S A. 2007;104:15156–61.View ArticlePubMedPubMed CentralGoogle Scholar
- Gendreau E, Höfte H, Grandjean O, Brown S, Traas J. Phytochrome controls the number of endoreduplication cycles in the Arabidopsis thaliana hypocotyl. Plant J. 1998;13:221–30.View ArticlePubMedGoogle Scholar
- Sugimoto-Shirasu K, Roberts K. “Big it up”: Endoreduplication and cell-size control in plants. Curr Opin Plant Biol. 2003;6:544–53.View ArticlePubMedGoogle Scholar
- Buenrostro JD, Giresi PG, Zaba LC, Chang HY, Greenleaf WJ. Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position. Nat Methods. 2013;10:1213–8. doi:10.1038/nmeth.2688.View ArticlePubMedPubMed CentralGoogle Scholar
- Oh E, Zhu J-Y, Wang Z-Y. Interaction between BZR1 and PIF4 integrates brassinosteroid and environmental responses. Nat Cell Biol. 2012;14:802–9. doi:10.1038/ncb2545.View ArticlePubMedPubMed CentralGoogle Scholar
- Ariel F, Jegu T, Latrasse D, Romero-Barrios N, Christ A, Benhamed M, et al. Noncoding transcription by alternative RNA polymerases dynamically regulates an auxin-driven chromatin loop. Mol Cell. 2014;55:383–96.View ArticlePubMedGoogle Scholar
- Jegu T, Latrasse D, Delarue M, Mazubert C, Bourge M, Hudik E, et al. Multiple functions of Kip-related protein5 connect endoreduplication and cell elongation. Plant Physiol. 2013;161:1694–705. doi:10.1104/pp.112.212357.View ArticlePubMedPubMed CentralGoogle Scholar
- Doherty CJ, Kay SA. Circadian control of global gene expression patterns. Annu Rev Genet. 2010;44:419–44. doi:10.1146/annurev-genet-102209-163432.View ArticlePubMedPubMed CentralGoogle Scholar
- McClung CR. The genetics of plant clocks. 1st edition. Elsevier Inc.; 2011. doi:10.1016/B978-0-12-387690-4.00004-0.
- Rugnone ML, Faigón Soverna A, Sanchez SE, Schlaen RG, Hernando CE, Seymour DK, et al. LNK genes integrate light and clock signaling networks at the core of the Arabidopsis oscillator. Proc Natl Acad Sci U S A. 2013;110:12120–5. doi:10.1073/pnas.1302170110.View ArticlePubMedPubMed CentralGoogle Scholar
- Yuan G-C, Liu Y-J, Dion MF, Slack MD, Wu LF, Altschuler SJ, et al. Genome-scale identification of nucleosome positions in S. cerevisiae. Science. 2005;309:626–30. doi:10.1126/science.1112178.View ArticlePubMedGoogle Scholar
- Chodavarapu RK, Feng S, Bernatavichute YV, Chen P-Y, Stroud H, Yu Y, et al. Relationship between nucleosome positioning and DNA methylation. Nature. 2010;466:388–92.View ArticlePubMedPubMed CentralGoogle Scholar
- Yen K, Vinayachandran V, Franklin Pugh B. Genome-wide structural integration of SWR-C and INO80 chromatin remodelers at +1 nucleosomes. Cell. 2013;154:1246–56.View ArticlePubMedPubMed CentralGoogle Scholar
- Boeger H, Griesenbeck J, Strattan JS, Kornberg RD. Nucleosomes unfold completely at a transcriptionally active promoter. Mol Cell. 2003;11:1587–98.View ArticlePubMedGoogle Scholar
- Li C, Gu L, Gao L, Chen C, Wei C-Q, Qiu Q, Chien C-W, Wang S, Jiang L, Ai L-F, Chen C-Y, Yang S, Nguyen V, Qi Y, Snyder MP, Burlingame AL, Chen X, Cui Y. Concerted genomic targeting of H3K27 demethylase REF6 and chromatin-remodeling ATPase BRM in Arabidopsis. Nat Genet. 2016;48:687–93.View ArticlePubMedPubMed CentralGoogle Scholar
- Leivar P, Monte E, Oka Y, Liu T, Carle C, Castillon A, et al. Multiple phytochrome-interacting bHLH transcription factors repress premature seedling photomorphogenesis in darkness. Curr Biol. 2008;18:1815–23. doi:10.1016/j.cub.2008.10.058.View ArticlePubMedPubMed CentralGoogle Scholar
- Wu C, Bingham PM, Livak KJ, Holmgren R, Elgin SCR. The chromatin structure of specific genes: I. Evidence for higher order domains of defined DNA sequence. Cell. 1979;16:797–806.View ArticlePubMedGoogle Scholar
- Wu C. The 5′ ends of Drosophila heat shock genes in chromatin are hypersensitive to DNase I. Nature. 1980;286:854–60.View ArticlePubMedGoogle Scholar
- Keene MA, Corces V, Lowenhaupt K, Elgin SC. DNase I hypersensitive sites in Drosophila chromatin occur at the 5′ ends of regions of transcription. Proc Natl Acad Sci U S A. 1981;78:143–6. doi:10.1073/pnas.78.1.143.View ArticlePubMedPubMed CentralGoogle Scholar
- Heo JB, Sung S, Kim DH, Doyle MR, Sung S, Amasino RM, et al. Vernalization-mediated epigenetic silencing by a long intronic noncoding RNA. Science. 2011;331:76–9. doi:10.1126/science.1197349.View ArticlePubMedGoogle Scholar
- Tsukiyama T, Wu C. Purification and properties of an ATP-dependent nucleosome remodeling factor. Cell. 1995;83:1011–20.View ArticlePubMedGoogle Scholar
- Sudarsanam P, Winston F. The Swi/Snf family. Trends Genet. 2000;16:345–51. doi:10.1016/S0168-9525(00)02060-6.View ArticlePubMedGoogle Scholar
- Varga-Weisz P. ATP-dependent chromatin remodeling factors: nucleosome shufflers with many missions. Oncogene. 2001;20:3076–85.View ArticlePubMedGoogle Scholar
- Gendrel A-V, Lippman Z, Martienssen R, Colot V. Profiling histone modification patterns in plants using genomic tiling microarrays. Nat Methods. 2005;2:213–8.View ArticlePubMedGoogle Scholar
- Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, et al. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008;9:R137.View ArticlePubMedPubMed CentralGoogle Scholar
- Shen L, Shao N, Liu X, Nestler E. ngs.plot: Quick mining and visualization of next-generation sequencing data by integrating genomic databases. BMC Genomics. 2014;15:284.View ArticlePubMedPubMed CentralGoogle Scholar
- Zang C, Schones DE, Zeng C, Cui K, Zhao K, Peng W. A clustering approach for identification of enriched domains from histone modification ChIP-Seq data. Bioinformatics. 2009;25:1952–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen K, Xi Y, Pan X, Li Z, Kaestner K, Tyler J, et al. DANPOS: Dynamic analysis of nucleosome position and occupancy by sequencing. Genome Res. 2013;23:341–51.View ArticlePubMedPubMed CentralGoogle Scholar
- Shen L, Shao NY, Liu X, Maze I, Feng J, Nestler EJ. diffReps: Detecting differential chromatin modification sites from ChIP-seq data with biological replicates. PLoS One. 2013;8:1–13.Google Scholar