A novel ATAC-seq approach reveals lineage-specific reinforcement of the open chromatin landscape via cooperation between BAF and p63
© Bao et al. 2015
Received: 18 August 2015
Accepted: 19 November 2015
Published: 18 December 2015
Open chromatin regions are correlated with active regulatory elements in development and are dysregulated in diseases. The BAF (SWI/SNF) complex is essential for development, and has been demonstrated to remodel reconstituted chromatin in vitro and to control the accessibility of a few individual regions in vivo. However, it remains unclear where and how BAF controls the open chromatin landscape to regulate developmental processes, such as human epidermal differentiation.
Using a novel “on-plate” ATAC-sequencing approach for profiling open chromatin landscapes with a low number of adherent cells, we demonstrate that the BAF complex is essential for maintaining 11.6 % of open chromatin regions in epidermal differentiation. These BAF-dependent open chromatin regions are highly cell-type-specific and are strongly enriched for binding sites for p63, a master epidermal transcription factor. The DNA sequences of p63 binding sites intrinsically favor nucleosome formation and are inaccessible in other cell types without p63 to prevent ectopic activation. In epidermal cells, BAF and p63 mutually recruit each other to maintain 14,853 open chromatin regions. We further demonstrate that BAF and p63 cooperatively position nucleosomes away from p63 binding sites and recruit transcriptional machinery to control tissue differentiation.
BAF displays high specificity in controlling the open chromatin landscape during epidermal differentiation by cooperating with the master transcription factor p63 to maintain lineage-specific open chromatin regions.
KeywordsBAF Differentiation Nucleosomes Open chromatin p63
Somatic tissue differentiation is an essential process that restricts cellular potential and enhances cell fate for specialized functions, such as the barrier function of epidermal tissue. Failure to engage tissue-specific differentiation underlies a spectrum of human diseases, including lipodystrophies , psoriasis , and cancers [3, 4]. It remains incompletely understood, however, how cell fate is determined and is restricted to the differentiation state.
An increasingly appreciated correlate of cell fate is the cell-type-specific open chromatin landscape , as defined by the specific regions of chromatin exhibiting high DNA accessibility in a given cell or tissue type. These open chromatin regions span an average of 300 base pairs , comprise 1–3 % of the total genome , and can be identified by DNase I hypersensitive sites sequencing (DNase-seq)  or formaldehyde-assisted isolation of regulatory elements sequencing (FAIRE-seq) . Comparison of open chromatin landscapes across different cell types has revealed distinct patterns. In contrast to open sites that are widely shared across diverse cell types, which can be near commonly active transcription start sites (TSSs) of housekeeping genes , many cell-type-specific open chromatin regions are intergenic and may correspond to tissue-specific regulatory regions.
The BAF chromatin remodeling complex, composed of one catalytic subunit (BRG1 or BRM) and 10 regulatory subunits, is essential for the differentiation of multiple somatic tissues such as muscle, fat, blood, and epidermis [10–15]. BAF is also the most frequently mutated chromatin-modifying complex in all human cancers [16, 17]. Using in vitro reconstituted chromatin and at only selected genomic loci in vivo, the BAF complex has been demonstrated to alter DNase I footprinting and to facilitate the DNA binding of a variety of transcription factors (TFs), including zinc-finger, helix-loop-helix, and rel-domain containing proteins [11, 18, 19]. However, it remains unclear how BAF controls the open chromatin landscape on a genome-wide scale, and whether BAF can systematically prioritize its interaction with specific TFs in a given developmental process such as tissue differentiation.
There are several technical barriers to understanding the control of open chromatin by BAF on a genome-wide scale. First, the BAF-dependent open chromatin regions cannot be simply represented by its chromatin immunoprecipitation sequencing (ChIP-seq) data, because the chromatin binding of the BAF complex can be mediated by its non-catalytic subunits in addition to the catalytic subunits, and because many BAF binding sites are also controlled redundantly by CHD and ISWI chromatin remodelers . Second, both DNase-seq and FAIRE-seq typically require tens of millions of cells [6, 7]. This amount of input material is more challenging in the context of BAF loss, which impairs the proliferation of many cell types . Lastly, the two BAF catalytic subunits BRG1 and BRM have been demonstrated to share redundant functions in somatic tissue development [12, 22]. However, limited genomic characterizations so far have only been attempted in the context of altered Brg1 function alone. By expressing a dominant-negative Brg1 mutant, only 1.9 % DNase hypersensitive sites (DHSs) were lost in the 3134 mouse mammary epithelial cell line . In mouse CD36+ cells, knocking down Brg1 alone also only modestly affected the nucleosome positioning at GATA1 sites without affecting GATA1 genomic binding . It remains unclear, however, how the open chromatin landscape could be altered with the loss of both BRG1 and BRM.
Here we address the role of the BAF complex in controlling the genome-wide open chromatin landscape using a model of human epidermal differentiation. We and others have demonstrated previously that loss of the BAF catalytic subunits BRG1 and BRM abolished epidermal tissue differentiation [12, 13], yet the mechanism underlying how BAF controls epidermal differentiation remains unclear. To characterize the open chromatin landscape in the context of BAF loss, we leveraged the recent development of the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq), which requires only 50,000 cells as input material . Because the original protocol was optimized for suspension cells, we further innovated a novel “on-plate ATAC-seq” method, which allows direct probing of open chromatin regions in adherent epidermal cells growing on a 96-well plate, bypassing the disruptive trypsinization process [24, 25]. Our high-quality data produced using this on-plate ATAC-seq method allowed exploration of the impact of the BAF complex on the open chromatin landscape during epidermal differentiation.
By characterizing the open chromatin landscape in the context of depleting both BAF catalytic subunits BRG1 and BRM, here we show that BAF is required to maintain 17,617 open chromatin sites (11.6 %) during epidermal cell differentiation, despite the fact that BAF binding is enriched in 50.3 % of open chromatin regions. BAF-dependent open chromatin regions are strikingly epidermal cell-type specific. They are also strongly enriched for binding sites of the epidermal master TF p63. p63 binding elements are enriched in ectodermal-derived tissue types . In keratinocytes, p63 binding is enriched in enhancers, and has a strong correlation with the histone modification H3K4me1 . However, it remains incompletely understood how p63 binding to DNA is facilitated in the context of chromatin. We show that both BAF and p63 co-regulate the expression of 236 epidermal differentiation genes, that they mutually reinforce each other’s binding, and that they are both required for sustaining open chromatin formation at p63 binding sites. Thus, these results demonstrate a new mode of cooperative action between BAF and p63 in controlling the keratinocyte-specific open chromatin landscape to enable epidermal differentiation. These data suggest that tissue terminal differentiation requires cooperation of highly selective BAF chromatin remodeling activity and a lineage-specific TF to sustain a tissue-specific open chromatin landscape.
BAF catalytic subunits are required to maintain open chromatin regions
The catalytic subunits of the BAF chromatin remodeling complex BRG1/BRM are essential for epidermal differentiation [12, 13]. We first reconfirmed this in primary human keratinocyte (NHEK) culture in vitro, where loss of BRG1/BRM strongly inhibited the induction of differentiation genes, such as KRT1, KRT10, and LCE3D (Additional file 1: Figure S1a). During the time course of 4-day keratinocyte differentiation, KRT1 expression was strongly induced in the control condition but not with BRG1/BRM loss (Additional file 1: Figure S1c). Differentiation day 4 was used as the time point for comparative genomic characterization in this study. RNA-seq analysis indicated that the gene set controlled by BAF is highly enriched with Gene Ontology (GO) terms such as “epidermis development” and “epidermal cell differentiation” (Additional file 1: Figure S1b; Additional file 2: Table S1). To interrogate the role of BAF in genome accessibility during differentiation, we used ATAC-seq to assess genome-wide accessibility to the Tn5 transposase, which acts in this context similarly to DNase I , in normal cells and in cells depleted of BRG1/BRM to remove functional catalytic subunits of the BAF complex (BAFi).
Because open chromatin regions comprise both cell-type-specific regions as well as common regions, we first compared both BAF-dependent and independent ATAC regions with the ENCODE DNase-seq open chromatin regions from other cell types. While the gained regions as well as the BAF-independent regions did not display strong cell-type specificity, the lost peaks with BAFi were significantly enriched with keratinocyte-specific open chromatin regions (p < 10−600; Fig. 1g, h). These data suggest that BAF is essential to maintain the cell-type-specific open chromatin landscape during epidermal differentiation.
BAF is required to maintain open chromatin regions that bind p63
To further characterize these BAF-dependent open chromatin regions, we performed de novo motif enrichment analysis on the 17,617 open chromatin regions that were lost with BAFi and the 7374 regions that gained ATAC-seq peaks with BAFi. The top motif enriched in BAF-dependent open chromatin corresponded to the p63 motif based on p63 ChIP-seq data (Fig. 2d–f) [30, 31]. p63 is the lineage-specifying TF essential for epidermal development as well as the TF required post-natally to maintain stratified epithelial tissue identity and terminal differentiation [32, 33]. Other p53 family member TFs, such as p53 and p73, were also expressed in epidermal cells; however, p63 was expressed at higher levels than these other two genes (Additional file 1: Figure S2b). p63 motif enrichment was striking (p = 1 × 10−1786), but only in genomic regions that lost openness with BAFi; it was not enriched at sites that gained accessibility with BAFi (Fig. 2d). Further analysis of ATAC-seq accessibility at individual p63 motif sites in p63 binding sites indicated that 21.1 % of these sites became inaccessible after BAF loss (Additional file 1: Figure S2c). Binding regions for 808 other TFs that are expressed at RPKM > 1, as measured by RNA-seq (Additional file 1: Figure S2d), including many that are functionally important in epidermal differentiation, such as HOPX1, CTNNB1, FOSL2, ZNF750, GRHL3, YAP1, RBPJ, OVOL1, MAF, MAFB, CEBPA, RARG, and KLF4, did not display such enrichment (Fig. 2d). ATAC-seq analysis further demonstrated a selective loss of accessibility at p63 motif regions with BAF depletion, but not at CTCF motif regions (Fig. 2g). In addition, ATAC-seq accessibility at KLF4 motif regions was not significantly altered at KLF4 motif regions (Additional file 1: Figure S2e). These data suggest that BAF maintains open chromatin regions with a strong selectivity for p63 binding in epidermal differentiation.
The dependence on BAF to maintain open chromatin regions containing p63 binding sites suggested that BAF and p63 may cooperate to control the expression of epidermal differentiation genes. To test this, we compared the target gene sets of BAF and p63 using RNA-seq. This demonstrated that BAF and p63-dependent target genes (fold change > 3 with depletion, FDR < 0.01) displayed highly significant overlap (p = 2.1 × 10−86), with 236 genes controlled jointly by both BAF and p63 (Fig. 2h; Additional file 1: Figure S2f, g; Additional file 2: Table S2). The GO terms for this shared target gene set were enriched with “epidermis development” (p = 4.8 × 10−5) and “epidermal cell differentiation” (p = 8.7 × 10−3) (Fig. 2i). Many of these shared genes were suppressed or induced in both BAFi and p63i conditions (Additional file 1: Figure S2e, f), suggesting that BAF and p63 may co-regulate epidermal tissue differentiation in a process that involves BAF-dependent open chromatin regions enriched for p63 binding.
BAF regulates nucleosome positioning and recruits transcriptional machinery
Additionally, in normal differentiation, there was high accessibility at the p63 motif sites; however, this was also decreased with BAF loss. To quantify the changes observed in these two conditions, we used the nucleosome-free fragments (<100 base pairs [bp]) and mononucleosome fragments (180–247 bp), as published previously for nucleosome positioning analysis . Consistent with the V-plot analysis, a decreased nucleosome-free signal was observed at p63 binding sites, indicating impaired accessibility at p63 motif regions and diminished nucleosome phasing adjacent to p63 motif with BAF loss. In addition, we found that the summits of the mononucleosome signals were shifted 20 bp towards the center of the p63 motif by BAF depletion (Fig. 3b), indicating an average of 40-bp nucleosome compaction at p63 motif regions in the absence of BAF. This altered nucleosome positioning with BAF loss was not universal, as no significant changes were observed at the CTCF motif (Additional file 1: Figure S3b-e). Therefore, nucleosome positioning analysis indicates that BAF is required to maintain genome accessibility at p63 binding motifs by positioning the nucleosomes away from the recognition sites.
We next hypothesized that BAF might help recruit transcriptional machinery and promote activation marks to sustain open chromatin at p63 binding regions. To test this, we performed ChIP-seq for the activating histone mark H3K27Ac, the histone acetyltransferase p300, RNA pol II, and the repressive histone mark H3K27me3, comparing control differentiating keratinocytes with BAF loss. This demonstrated that H3K27Ac, p300, and pol II levels were decreased with BAF loss, while H3K27me3 levels were increased (Fig. 3c–g). Therefore BAF is required for normal levels of p300, pol II, and H3K27Ac at p63 sites, where it is also required to oppose H3K27me3 and maintain open chromatin.
p63 is also required to maintain open chromatin regions
Because local DNA sequence can directly impact nucleosome positioning, we calculated the intrinsic nucleosome binding probability based on genomic sequence [36, 37]. The genomic DNA sequences containing p63 motif sites displayed high nucleosome binding probability, which peaked at p63 motif sites with a probability of ~80 % (Fig. 4e). This high probability for nucleosome binding was absent in the shuffled genomic regions (Additional file 1: Figure S3f). These data suggest that the DNA sequences of p63 binding sites strongly favor nucleosome formation, and that these regions are intrinsically inaccessible in the absence of DNA binding proteins and chromatin remodelers. These data are consistent with our observation that both BAF and p63 are required to maintain open chromatin at p63 sites and support the existence of strong intrinsic suppressive mechanisms at these sites to prevent ectopic epithelial gene activation in non-epithelial cell types.
BAF and p63 binding reinforce open chromatin
Because BAF loss caused nucleosome compaction and decreased accessibility at p63 binding sites, we hypothesized that BAF could be essential for normal p63 binding across the genome. To test this, we first compared ATAC-seq accessibility at single nucleotide resolution at p63 and CTCF motif sites, as a method to predict the binding of these TFs. We observed a decreased signal for p63, but not CTCF with BAF loss (Additional file 1: Figure S4a, b), suggesting that p63 DNA binding is dependent on BAF. p63 levels did not change significantly with BAF loss (Additional file 1: Figure S1c), confirming that impaired p63 binding in this setting was not due to decreased p63 protein levels. To validate these findings, p63 ChIP-seq was performed in differentiating human keratinocytes with and without BAF depletion. p63 ChIP-seq peaks in both control and BAF-loss conditions were called using MACS (version 1.4.2, p-value cutoff 1.00e − 05). Out of a total of 30,075 normal control p63 binding peaks, 11,930 fell below the MACS peak-calling threshold with BAF loss. Even among the 18,145 peaks still called by MACS in both conditions, the p63 binding signal was reduced on average to half with BAF loss (Additional file 1: Figure S4c–e). In the 14,853 BAF and p63 co-bound ATAC-seq regions, BAF loss impaired both ATAC accessibility and p63 binding in these regions (Fig. 5c, d). These data indicate that normal p63 genomic binding requires intact BAF.
Prior work suggested that the BAF complex can itself be targeted by TFs to chromatin [19, 38]. To test whether p63 facilitates BAF binding in differentiating epidermal cells, we performed BAF ChIP-seq with p63 depletion. p63 loss correlated with reduced BAF binding at the p63 co-binding sites, but not at CTCF binding sites (Fig. 5c, d; Additional file 1: Figure S4f), indicating that p63 functions selectively on BAF association at p63 binding sites but not universally on BAF association with other sites.
Taken together, our data demonstrate that BAF maintains the epidermal cell-specific open chromatin landscape to promote tissue differentiation. In epithelial tissue, these BAF-dependent open chromatin regions are specifically enriched with p63 sites, which are intrinsically inaccessible to transcriptional machinery in a manner that may ensure repression of epithelial-specific gene expression in non-epithelial tissues that lack p63 expression. The chromatin remodeling activity of BAF controls nucleosome position to increase the accessibility of p63 binding sites, facilitating p63 genomic binding. The binding of p63, in turn, also enhances the chromatin targeting of BAF. This cooperation between BAF and p63 further enables local association of p300 and pol II, sustains activating histone marks, and opposes repressive histone marks in concert with promoting open chromatin formation and differentiation gene expression.
Here we present data indicating that the BAF complex and the epidermal lineage-determining TF p63 cooperate to remodel the open chromatin landscape and condition p63 binding regions for gene activation. In this setting, BAF chromatin remodeling catalytic subunits cooperate with p63 to control nucleosome positioning and to promote epidermal tissue-specific gene expression during differentiation. In this context, BAF is required to displace nucleosomes at an average of approximately 40 bp laterally away from p63 binding sites. This BAF-dependent genome accessibility at p63 sites is associated with enabled p63 binding, activating histone marks, RNA polymerase II machinery recruitment, and target gene expression. Epidermal tissue differentiation is a dynamic process that involves the repression of progenitor genes and the induction of differentiation genes . Although BRG1/BRM and p63 are expressed in both progenitor and differentiation states, our data generation and analysis in this study only focused on differentiation day 4. It remains to be confirmed whether this cooperation also occurs in the epidermal progenitors to maintain the progenitor state.
The BAF complex has long been suggested to slide nucleosomes and to enhance TF binding to DNA, with early studies using in vitro DNase I cleavage experiments to probe changes on reconstituted chromatin templates in vitro [40, 41]. Using similar in vitro systems, it has also been indicated that transactivation domains, such as VP16, could facilitate BAF targeting to chromatin to stimulate its remodeling activity in vitro . Very little has been determined, however, regarding how the BAF complex remodels genome-wide open chromatin landscape to regulate gene expression in vivo . The present work indicates that BAF maintains DNA accessibility at ~11.6 % of open genomic regions in epidermal differentiation, with a strong selectivity towards lineage-specific regions rather than the common regions shared by multiple cell types. This suggests that the BAF complex does not function as a general chromatin remodeler, but displays strong selectivity in the context of developmental processes.
BAF has been demonstrated to associate with a variety of TFs from both in vitro and in vivo evidence [16, 25, 40, 41, 44, 45]. With the presence of hundreds of transcription regulators at a given developmental stage, it was unclear whether BAF displays any specificity in interacting selectively at binding sites of particular functionally active TFs to control open chromatin regions. Our data suggest that, during epidermal differentiation, the BAF complex exhibits strong specificity in cooperating with the lineage-determining TF p63 to activate gene expression by maintaining chromatin accessibility. Previous motif analysis indicated that a number of other TFs, such as RUNX1, RXRA, and TFAP2a, may function as co-activators for p63 . This may explain the variations we observed for the accessibility changes in individual p63 sites with BAF knockdown. The co-binding of other TFs near p63 binding sites may make these sites more resistant to decreased p63 binding caused by BAF loss. In addition, it has been shown in mouse skin that p63 regulates Satb1 expression to further control chromatin architecture . It would be of interest for future investigations to elucidate the crosstalk and the feedback regulation between chromatin modifiers and TFs to control gene expression in the context of tissue differentiation.
Our data also suggest that the BAF complex is required to maintain the accessibility at slightly more than 10 % of open chromatin regions genome-wide, at least in the context of epidermal differentiation. This raises the question of which factors maintain the open structure of the other ~90 %. One plausible possibility is that the other regions are controlled by the remodeling activity of other chromatin remodeling complexes, such as CHD, ISWI, and INO80 family members . It also remains to be determined how much of the other open chromatin regions are pre-determined by their DNA sequence to dis-favor nucleosome formation.
We observed two well-positioned nucleosomes adjacent to p63 binding sites in normal epithelial differentiation. The role of nucleosome positioning in gene expression, while incompletely understood, has been best studied at TSSs. At the TSSs for active genes, a total of eight (three upstream and five downstream of the TSS) well-positioned nucleosomes are detected; however, only one +1 nucleosome is detected for silenced genes . This suggests that well-positioned nucleosomes are positively correlated with gene expression. Much less is understood for intergenic DNA regions, including the distal regulatory regions observed here. Similar to what we have observed at p63 binding sites, well-positioned nucleosomes near TF binding sites have been reported at the Pu.1 sites in macrophages  and the androgen receptor sites in prostate cancer cells , although the basis for this positioning is unclear. Our data suggest that nucleosome phasing at regulatory regions requires the cooperation of a dominant TF, such as p63, and a chromatin remodeling complex, such as BAF, to overcome the intrinsic proclivity to nucleosome binding.
The BAF chromatin remodeling complex is a 2-MDa protein complex composed of 10 regulatory subunits and the interchangeable catalytic subunit BRG1 or BRM. The expression of BAF subunits differs in various cell types, and specific BAF subunits function to control specific gene expression . To study the role of the BAF chromatin remodeling complex, we focused on its catalytic subunits. Taken together with the cooperativity and selectivity of BAF with p63 in epithelial tissue, it would be interesting to investigate the specific subunit composition of the BAF complex that interacts with p63. Furthermore, the BAF complex is frequently mutated in human cancers, including epidermal squamous cell carcinoma of the skin (SCC) . Testing whether p63-BAF cooperativity in differentiation gene induction is compromised in SCC in a manner that inhibits tissue terminal differentiation and promotes cancer progression is also an intriguing future avenue of research. Given that p53 and p63 are structurally very similar to each other, it is also of future interest to investigate the degree of sustained cooperation between p53 and the BAF complex in both normal and cancer tissues.
We demonstrate that the BAF chromatin remodeling complex plays an essential role in controlling the open chromatin landscape in somatic tissue differentiation. Although BAF binds broadly across the open chromatin landscape, BAF displays strong preferences in controlling the lineage-specific regions by cooperating with the master TF p63 in epidermal differentiation. These lineage-specific regions are characterized by a high intrinsic affinity to nucleosomes, and are inaccessible in other lineages. Maintenance of accessibility at the lineage-specific regions requires cooperation between BAF and the master TF to position the nucleosomes away from the TF recognition sites and to further recruit the transcriptional machinery to these regions, enabling tissue-specific gene expression.
Epidermal cell culture
Primary human epidermal keratinocytes were isolated from fresh surgically discarded skin from multiple donors. Keratinocytes from a minimum of three donors were mixed and cultured in complete Keratinocyte-SFM (Life Technologies, Grand Island, New York, USA, #17005-142) and Medium 154 (Life Technologies #M-154-500). To induce terminal differentiation, keratinocytes were seeded in confluent condition with the addition of 1.2 mM Ca2+ in culture medium for 4 days. All sequencing experiments in this paper were carried out on differentiation day 4.
On-plate ATAC-seq for differentiated keratinocytes
The on-plate ATAC-seq was modified based on the original ATAC-seq method  with the following modifications. First, 50,000 keratinocytes were seeded on a 96-well plate (Breiner Bio-one, Monroe, North Carolina, USA, Cat# 675180). After inducing differentiation for 4 days, the cells were washed briefly with phosphate-buffered saline (PBS) two times, and permeabilized with lysis buffer (10 mM Tris-Cl pH 7.4, 10 mM NaCl, 3 mM MgCl2, 0.05 % NP-40) for 10 min. After this permeabilization step, the supernatant were removed, and replaced with 50 μL transposase and buffer mix (25 μL TD (2x reaction buffer from Nextera kit, Illumina, San Diego, California, USA), 5 μL TDE1(Tn5 Transposase from Nextera kit, Illumina), 20 μL nuclease-free H2O). After 30 min of incubation at 37 °C, a total volume of 250 μL Buffer PB (Qiagen, Hilden, Germany) was mixed with the reaction before DNA purification and polymerase chain reaction (PCR) enrichment. Each replicate of ATAC-seq data used in this study was generated individually using 50,000 keratinocytes growing in each well of a 96-well plate under differentiation conditions.
Gene transfer and knockdown
For siRNA knockdown of BRG1/BRM, 1 × 106 cells were electroporated with 1 nmol siRNA of each (ON-TARGETplus SMARTpool, Dharmacon, Lafayette, Colorado, USA) using Amaxa Human Keratinocyte Nucleofector Kit (Lonza, Basel, Switzerland, VPD-1002). To serve as control, 2 nmol total control siRNA were electroporated to 1 × 106 cells. For siRNA knockdown of p63, an siRNA targeting pan p63 was used as published previously (sequence: CGACAGUCUUGUACAAUUUUU, Dharmacon) .
Protein expression analysis
For immunoblot analysis, 20–50 μg of cell lysates were loaded per lane for sodium dodecyl sulfate polyacrylamide gel electrophoresis and transferred to polyvinylidene fluoride membranes. Primary antibodies were incubated at 4 °C overnight and secondary antibodies were incubated at room temperature for 1 h. The affinity purified J1 antiserum against BRG1/BRM was a generous gift of G Crabtree (Stanford University). Other antibodies used in this study include: anti-Krt1 (Covance, Princeton, New Jersey, USA), p63 (Santa Cruz Biotechnology, Santa Cruz, California, USA), p53 (Santa Cruz Biotechnology), and tubulin (Sigma, St. Louis, Missouri, USA).
qRT-PCR expression analysis
For qRT-PCR, total RNA was extracted using the RNeasy Plus kit (Qiagen) and subsequently subjected to reverse transcription using the SuperScript VILO cDNA synthesis kit (Invitrogen). qRT-PCR (Quantitative Reverse Transcription Polymerase Chain Reaction) analysis was performed using the Mx3000P instrument with the SYBR Green Master Mix (Fermentas, Burlington, Ontario, Canada). Samples were run in duplicate and normalized to levels of GAPDH mRNA or 18S ribosomal RNA for each reaction.
ChIP assays were performed essentially as previously described  with minor modifications. Human keratinocytes were cross-linked with 1 % formaldehyde at room temperature for 10 min. The chromatin was sonicated to achieve fragments with an average length of 200–500 bp. The sonicated chromatin was immunoprecipitated overnight at 4 °C with J1 antibody , Poll II (Millipore, Bedford, Massachusetts, USA), H3K27Ac (Active Motif, Calsbad, California, USA), H3K27me3 (Active Motif), H3K4me (Active Motif), or the same amount of IgG control (Cell Signaling). Following reverse cross-linking, the samples were treated with RNase and Protease K, and the DNA was purified using the Qiagen PCR Purification kit. For ChIP-seq, 10 ng of purified ChIP DNA were used to generate the sequencing library using NEBNext ChIP-Seq Library Prep Master Mix Set.
Analysis of cell-type specificity of chromatin accessibility
Filtered reads from DHS-seq were downloaded as BAM files from the ENCODE project (http://hgdownload.cse.ucsc.edu/goldenPath/hg19/encodeDCC/wgEncodeUwDnase/). A set of universal DHS “clusters” identified by ENCODE was download from the UCSC Genome Browser (wgEncodeRegDnaseClusteredV3.bed), and DHS reads from each cell type were counted in each universal DHS region. DHS read counts were transformed with the arcsin function, and a z-score was calculated for the signal of each cell type for each region. Box plots were made using the R boxplot() function, and heatmaps were made using the heatmap.2 function in the R gplots package.
Nucleosome probability analysis
Calculation of intrinsic nucleosome binding based on DNA sequence was performed using the Perl algorithm developed by Segal Lab . Sequences of 6000 bp flanking the p63 motif sites were used as input genomic sequences. To remove the boundary effect, a window size of 4000 bp flanking the motif sites was used to plot the output. The same number and length of randomly shuttled genomic regions generated by Bedtools were used as a negative control for the analysis.
Proximity ligation analysis
Primary human keratinocytes were cultured in differentiation condition on chamber slides, fixed by 4 % Formaldehyde Solution (Thermo Scientific, Rockford, Illinois, #28906) for 10 min at room temperature. Slides were blocked using blocking buffer (PBS with 0.3 % TritonX-100, 3 % normal horse serum) for 1 h at room temperature, followed by primary antibody incubation overnight at 4 °C. Slides were washed by PBS three times the next day, and were subsequently incubated with the PLA probes anti-mouse and anti-rabbit (Sigma). The Duolink In Situ Detection Reagents Orange (Sigma) was used for ligation and signal amplification.
Nucleosome-free and mononucleosome ATAC-seq analysis
We first isolated paired-end ATAC-seq fragments from two length ranges: fragments shorter than 100 bp, and fragment lengths between 180 and 247 bp. These ranges were used previously  to distinguish transposase insertion into nucleosome-free DNA and insertions flanking nucleosomes, respectively. A 1-kb window centered on each region of interest was divided into ten 100 equal-sized bins. For each region, the number of fragment centers (single base pair) overlapping each bin was counted for each of the two fragment size ranges. The average fragment count plotted in each bin was normalized to the average fragment count in the first five bins in order to set the background signal to one. In order to facilitate comparisons across conditions in which ATAC-seq libraries exhibited markedly different signal to noise ratios, CTCF binding sites were used as a stable reference feature to scale average profiles. The enrichment profile was computed as described above for CTCF binding sites, and then a simple linear transformation was applied to match CTCF profiles across conditions. The same transformation was applied to profiles derived from other features of interest. For profiles derived from low-depth ATAC-seq libraries, a spline smoothing function from R (smooth.spline package) was used with a smoothing parameter between 0.3 and 0.5.
V-plot analysis for nucleosome positioning
A volcano plot was constructed by plotting the normalized number of fragments. The length of the fragment was indicated by the y-axis; and the distance between the centers of the fragments to a p63 motif site was indicated by the x-axis. Volcano plots for both control and BAF loss were generated, enabling visual comparison of nucleosome positioned near p63 motif sites.
ATAC-seq mapping, peak calls, and differential peak calls
ATAC-seq paired-end reads were trimmed for Illumina adapter sequences using a custom script and mapped to hg19 using bowtie with parameters –S –X2000 –m1. Duplicate reads were discarded with samtools. Peaks were identified using macs2 with parameters callpeak –nomodel –shift -100 –extsize 200. Peak sets called in individual replicates were combined and individual peaks merged if overlapping within 300 bp to form a union peak set. Differentially accessible peaks from this union set were identified with edgeR by counting all read ends overlapping peaks in each condition. edgeR was run with default settings, a fold-change threshold of 2, and FDR < 0.01.
Single-nucleotide-resolution accessibility analysis
Because the Tn5 transposon binds as a dimer and inserts two adaptors separated by 9 bp , all reads aligning to the + strand were offset by +4 bp, and all reads aligning to the – strand were offset −5 bp. To predict the DNA binding of p63 and CTCF, we calculated the enrichment of the first base pair of each mapped and shifted ATAC-seq read in the motif sites within the TF binding sites from ChIP-seq analysis. p63 ChIP-seq in differentiated human keratinocytes  as well as the ENCODE CTCF ChIP-seq using human keratinocytes were used in this analysis.
Cross cell type DHS analysis
DHS profiles were downloaded as BigWig files from the ENCODE portal at the UCSC genome browser (http://hgdownload.cse.ucsc.edu/goldenPath/hg19/encodeDCC/wgEncodeUwDnase/). The Homer annotatePeaks tool (http://homer.salk.edu/homer/) was used to calculate the average DHS signal in 10-bp bins surrounding p63 or CTCF sites.
Metaplots of ChIP-seq signal (H3K27ac, H3K27me3, H3K4me1, RNA polymerase II, and BAF complex) were generated using the Homer annotatePeaks tool with reads counted in a 10-bp window and default library size normalization.
RNA-seq library generation and data analysis
RNA extraction was performed using RNeasy Plus (Qiagen) from day 4 differentiated human keratinocytes treated with Brg1/Brm siRNA, p63 siRNA, or Control siRNA in technical replicates. rRNA was removed from each RNA extraction using Ribo-Zero rRNA Removal Kit (Illumina) before proceeding with RNA-seq library construction (NEBNext Ultra RNA Library Prep Kit, NEB). Paired-end sequencing of 2 × 75 bp was performed by Stanford Functional Genomics Facility using a Nextseq Sequencer (Illumina). Sequencing reads were mapped to hg19 using TopHat, raw read counts for each gene were calculated using an in-house generated perl script, and differential expression was analyzed using DEseq with the cutoff of FDR < 0.01, fold change >3.
Our ATAC-seq, ChIP-seq, and RNA-seq data sets have been deposited with the Gene Expression Omnibus [GEO: GSE67382]. The sequencing depth of these sequencing data is included in Additional file 2: Table S3.
The DHS data, RNA-seq data for 15 cell types as well as the CTCF ChIP-seq data in NHEK used in this study were generated by the ENCODE project. p63 ChIP-seq data in normal keratinocytes and the p63 transcriptome profiling analysis were extracted from previously published data sets generated by our group . KLF4 motif sites within KLF4 ChIP-seq binding sites and RNA-seq data of differentiated human keratinocytes were also extracted from previously published studies from our group [55, 56].
This study was conducted under IRB-35324 (Stanford University). The keratinocytes used in this study were isolated from surgically discarded foreskin before January 2015, and were mixed from 3–5 de-identified donors. The tissue donation was under IRB exempt protocol at Stanford. We confirm that the experimental methods comply with the Helsinki Declaration.
Assay for transposase-accessible chromatin with high throughput sequencing
Brg/Brm associated factor
chromatin immunoprecipitation followed by high-throughput DNA sequencing
DNaseI hypersensitivity site
Encyclopedia of DNA Elements
formaldehyde-assisted isolation of regulatory elements followed by high-throughput sequencing
false discovery rate
normal human epidermal keratinocytes
polymerase chain reaction
proximity ligation analysis
squamous cell carcinoma
transcription start site
We thank A. Oro, G. Crabtree, V. Lopez-Pajares, A. Zehnder, and B. Barajas for pre-submission review; and R. Li, Y. Shen, and S. Tao for technical support. We are grateful for the suggestions from W. Greenleaf, B. Zarnegar, R. Flynn, and B. Bajaras. We thank L. Morcom and P. Bernstein for administrative assistance. This work is supported by the U.S. Department of Veterans Affairs Office of Research and Development, by NIH R01 AR45192 and AR49737 (PAK), by NIH F32 and NIH K99 Award (AR06123, AR06548) to XB, and by NIH P50 HG007735 and U01 HG007919 to HYC and PAK.
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.
- Agarwal AK, Garg A. Genetic disorders of adipose tissue development, differentiation, and death. Annu Rev Genomics Hum Genet. 2006;7:175–99.PubMedView ArticleGoogle Scholar
- Thelu J, Rossio P, Favier B. Notch signalling is linked to epidermal cell differentiation level in basal cell carcinoma, psoriasis and wound healing. BMC Dermatol. 2002;2:7.PubMed CentralPubMedView ArticleGoogle Scholar
- Reuter JA, Ortiz-Urda S, Kretz M, Garcia J, Scholl FA, Pasmooij AMG, et al. Modeling inducible human tissue neoplasia identifies an extracellular matrix interaction network involved in cancer progression. Cancer Cell. 2009;15:477–88.PubMed CentralPubMedView ArticleGoogle Scholar
- Ridky TW, Chow JM, Wong DJ, Khavari PA. Invasive three-dimensional organotypic neoplasia from multiple normal human epithelia. Nat Med. 2010;16:1450–5.PubMed CentralPubMedView ArticleGoogle Scholar
- Stergachis AB, Neph S, Reynolds A, Humbert R, Miller B, Paige SL, et al. Developmental fate and cellular maturity encoded in human regulatory DNA landscapes. Cell. 2013;154:888–903.PubMed CentralPubMedView ArticleGoogle Scholar
- Natarajan A, Yardimci GG, Sheffield NC, Crawford GE, Ohler U. Predicting cell-type-specific gene expression from regions of open chromatin. Genome Res. 2012;22:1711–22.PubMed CentralPubMedView ArticleGoogle Scholar
- Song L, Zhang Z, Grasfeder LL, Boyle AP, Giresi PG, Lee B-K, et al. Open chromatin defined by DNaseI and FAIRE identifies regulatory elements that shape cell-type identity. Genome Res. 2011;21:1757–67.PubMed CentralPubMedView ArticleGoogle Scholar
- Boyle AP, Davis S, Shulha HP, Meltzer P, Margulies EH, Weng Z, et al. High-resolution mapping and characterization of open chromatin across the genome. Cell. 2008;132:311–22.PubMed CentralPubMedView ArticleGoogle Scholar
- Giresi PG, Kim J, McDaniell RM, Iyer VR, Lieb JD. FAIRE (Formaldehyde-Assisted Isolation of Regulatory Elements) isolates active regulatory elements from human chromatin. Genome Res. 2007;17:877–85.PubMed CentralPubMedView ArticleGoogle Scholar
- Hang CT, Yang J, Han P, Cheng HL, Shang C, Ashley E, et al. Chromatin regulation by Brg1 underlies heart muscle development and disease. Nature. 2010; 466:62-67.
- De la Serna IL, Carlson KA, Imbalzano AN. Mammalian SWI/SNF complexes promote MyoD-mediated muscle differentiation. Nat Genet. 2001;27:187–90.PubMedView ArticleGoogle Scholar
- Indra AK, Dupé V, Bornert J-M, Messaddeq N, Yaniv M, Mark M, et al. Temporally controlled targeted somatic mutagenesis in embryonic surface ectoderm and fetal epidermal keratinocytes unveils two distinct developmental functions of BRG1 in limb morphogenesis and skin barrier formation. Development. 2005;132:4533–44.PubMedView ArticleGoogle Scholar
- Bao X, Tang J, Lopez-Pajares V, Tao S, Qu K, Crabtree GR, et al. ACTL6a enforces the epidermal progenitor state by suppressing SWI/SNF-dependent induction of KLF4. Cell Stem Cell. 2013;12:193–203.PubMed CentralPubMedView ArticleGoogle Scholar
- Hu G, Schones DE, Cui K, Ybarra R, Northrup D, Tang Q, et al. Regulation of nucleosome landscape and transcription factor targeting at tissue-specific enhancers by BRG1. Genome Res. 2011;21:1650–8.PubMed CentralPubMedView ArticleGoogle Scholar
- Gresh L, Bourachot B, Reimann A, Guigas B, Fiette L, Garbay S, et al. The SWI/SNF chromatin-remodeling complex subunit SNF5 is essential for hepatocyte differentiation. EMBO J. 2005;24:3313–24.PubMed CentralPubMedView ArticleGoogle Scholar
- Kadoch C, Hargreaves DC, Hodges C, Elias L, Ho L, Ranish J, et al. Proteomic and bioinformatic analysis of mammalian SWI/SNF complexes identifies extensive roles in human malignancy. Nat Genet. 2013;45:592–601.PubMed CentralPubMedView ArticleGoogle Scholar
- Wilson BG, Roberts CWM. SWI/SNF nucleosome remodellers and cancer. Nat Rev Cancer. 2011;11:481–92.PubMedView ArticleGoogle Scholar
- Utley RT. SWI/SNF Stimulates the formation of disparate activator-nucleosome complexes but is partially redundant with cooperative binding. J Biol Chem. 1997;272:12642–9.PubMedView ArticleGoogle Scholar
- De la Serna IL, Ohkawa Y, Berkes CA, Bergstrom DA, Dacwag CS, Tapscott SJ, et al. MyoD targets chromatin remodeling complexes to the myogenin locus prior to forming a stable DNA-bound complex. Mol Cell Biol. 2005;25:3997–4009.PubMed CentralPubMedView ArticleGoogle Scholar
- Morris SA, Baek S, Sung M-H, John S, Wiench M, Johnson TA, et al. Overlapping chromatin-remodeling systems collaborate genome wide at dynamic chromatin transitions. Nat Struct Mol Biol. 2014;21:73–81.PubMed CentralPubMedView ArticleGoogle Scholar
- Moshkin YM, Mohrmann L, van Ijcken WFJ, Verrijzer CP. Functional differentiation of SWI/SNF remodelers in transcription and cell cycle control. Mol Cell Biol. 2007;27:651–61.PubMed CentralPubMedView ArticleGoogle Scholar
- Willis MS, Homeister JW, Rosson GB, Annayev Y, Holley D, Holly SP, et al. Functional redundancy of SWI/SNF catalytic subunits in maintaining vascular endothelial cells in the adult heart. Circ Res. 2012;111:e111–22.PubMed CentralPubMedView ArticleGoogle 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.PubMed CentralPubMedView ArticleGoogle Scholar
- Huang H, Hsing H, Lai T, Chen Y, Lee T, Chan H, et al. Trypsin-induced proteome alteration during cell subculture in mammalian cells. J Biomed Sci. 2010;17:36.PubMed CentralPubMedView ArticleGoogle Scholar
- Banno T, Blumenberg M. Keratinocyte detachment-differentiation connection revisited, or Anoikis-Pityriasi Nexus Redux. PLoS One. 2014;9:1–12.View ArticleGoogle Scholar
- Sethi I, Sinha S, Buck MJ. Role of chromatin and transcriptional co-regulators in mediating p63-genome interactions in keratinocytes. BMC Genomics. 2014;15:1042.PubMed CentralPubMedView ArticleGoogle Scholar
- Ho L, Jothi R, Ronan JL, Cui K, Zhao K, Crabtree GR. An embryonic stem cell chromatin remodeling complex, esBAF, is an essential component of the core pluripotency transcriptional network. Proc Natl Acad Sci U S A. 2009;106:5187–91.PubMed CentralPubMedView ArticleGoogle Scholar
- Euskirchen GM, Auerbach RK, Davidov E, Gianoulis TA, Zhong G, Rozowsky J, et al. Diverse roles and interactions of the SWI/SNF chromatin remodeling complex revealed using global approaches. PLoS Genet. 2011;7:e1002008.PubMed CentralPubMedView ArticleGoogle Scholar
- Khavari PA, Peterson CL, Tamkun JW, Mendel DB, Crabtree GR. BRG1 contains a conserved domain of the SWI2/SNF2 family necessary for normal mitotic growth and transcription. Nature. 1993;366:170–4.PubMedView ArticleGoogle Scholar
- Zarnegar BJ, Webster DE, Lopez-Pajares V, Vander Stoep Hunt B, Qu K, Yan KJ, et al. Genomic profiling of a human organotypic model of AEC syndrome reveals ZNF750 as an essential downstream target of mutant TP63. Am J Hum Genet. 2012;91:435–43.PubMed CentralPubMedView ArticleGoogle Scholar
- Smeenk L, van Heeringen SJ, Koeppel M, van Driel MA, Bartels SJJ, Akkers RC, et al. Characterization of genome-wide p53-binding sites upon stress response. Nucleic Acids Res. 2008;36:3639–54.PubMed CentralPubMedView ArticleGoogle Scholar
- Truong AB, Kretz M, Ridky TW, Kimmel R, Khavari PA. p63 regulates proliferation and differentiation of developmentally mature keratinocytes. Genes Dev. 2006;20:3185–97.PubMed CentralPubMedView ArticleGoogle Scholar
- Mills AA, Zheng B, Wang XJ, Vogel H, Roop DR, Bradley A. p63 is a p53 homologue required for limb and epidermal morphogenesis. Nature. 1999;398:708–13.PubMedView ArticleGoogle Scholar
- Armstrong JA, Bieker JJ, Emerson BM. A SWI/SNF–related chromatin remodeling complex, E-RC1, is required for tissue-specific transcriptional regulation by EKLF in vitro. Cell. 1998;95:93–104.PubMedView ArticleGoogle Scholar
- Martens JA, Winston F. Recent advances in understanding chromatin remodeling by Swi/Snf complexes. Curr Opin Genet Dev. 2003;13:136–42.PubMedView ArticleGoogle Scholar
- Kaplan N, Moore IK, Fondufe-Mittendorf Y, Gossett AJ, Tillo D, Field Y, et al. The DNA-encoded nucleosome organization of a eukaryotic genome. Nature. 2009;458:362–6.PubMed CentralPubMedView ArticleGoogle Scholar
- Segal E, Fondufe-Mittendorf Y, Chen L, Thåström A, Field Y, Moore IK, et al. A genomic code for nucleosome positioning. Nature. 2006;442:772–8.PubMed CentralPubMedView ArticleGoogle Scholar
- Kowenz-Leutz E, Leutz A. A C/EBPβ isoform recruits the SWI/SNF complex to activate myeloid genes. Mol Cell. 1999;4:735–43.PubMedView ArticleGoogle Scholar
- Lopez-Pajares V, Qu K, Zhang J, Webster DE, Barajas BC, Siprashvili Z, et al. A LncRNA-MAF:MAFB transcription factor network regulates epidermal differentiation. Dev Cell. 2015;32:693–706.PubMedView ArticleGoogle Scholar
- Cote J, Peterson CL, Workman JL. Perturbation of nucleosome core structure by the SWI/SNF complex persists after its detachment, enhancing subsequent transcription factor binding. Proc Natl Acad Sci. 1998;95:4947–52.PubMed CentralPubMedView ArticleGoogle Scholar
- Kwon H, Imbalzano AN, Khavari PA, Kingston RE, Green MR. Nucleosome disruption and enhancement of activator binding by a human SW1/SNF complex. Nature. 1994;370:477–81.PubMedView ArticleGoogle Scholar
- Gutiérrez JL, Chandy M, Carrozza MJ, Workman JL. Activation domains drive nucleosome eviction by SWI/SNF. EMBO J. 2007;26:730–40.PubMed CentralPubMedView ArticleGoogle Scholar
- Hargreaves DC, Crabtree GR. ATP-dependent chromatin remodeling: genetics, genomics and mechanisms. Cell Res. 2011;21:396–420.PubMed CentralPubMedView ArticleGoogle Scholar
- Seo S, Richardson GA, Kroll KL. The SWI/SNF chromatin remodeling protein Brg1 is required for vertebrate neurogenesis and mediates transactivation of Ngn and NeuroD. Development. 2005;132:105–15.PubMedView ArticleGoogle Scholar
- Pedersen TA, Kowenz-Leutz E, Leutz A, Nerlov C. Cooperation between C/EBPalpha TBP/TFIIB and SWI/SNF recruiting domains is required for adipocyte differentiation. Genes Dev. 2001;15:3208–16.PubMed CentralPubMedView ArticleGoogle Scholar
- Kouwenhoven EN, Oti M, Niehues H, van Heeringen SJ, Schalkwijk J, Stunnenberg HG, et al. Transcription factor p63 bookmarks and regulates dynamic enhancers during epidermal differentiation. EMBO Rep. 2015;16:863–78.PubMed CentralPubMedView ArticleGoogle Scholar
- Fessing MY, Mardaryev AN, Gdula MR, Sharov AA, Sharova TY, Rapisarda V, et al. p63 regulates Satb1 to control tissue-specific chromatin remodeling during development of the epidermis. J Cell Biol. 2011;194:825–39.PubMed CentralPubMedView ArticleGoogle Scholar
- Clapier CR, Cairns BR. The biology of chromatin remodeling complexes. Annu Rev Biochem. 2009;78:273-304.
- Schones DE, Cui K, Cuddapah S, Roh T-Y, Barski A, Wang Z, et al. Dynamic regulation of nucleosome positioning in the human genome. Cell. 2008;132:887–98.PubMedView ArticleGoogle Scholar
- Barozzi I, Simonatto M, Bonifacio S, Yang L, Rohs R, Ghisletti S, et al. Coregulation of transcription factor binding and nucleosome occupancy through DNA features of mammalian enhancers. Mol Cell. 2014;54:844–57.PubMed CentralPubMedView ArticleGoogle Scholar
- He HH, Meyer CA, Shin H, Bailey ST, Wei G, Wang Q, et al. Nucleosome dynamics define transcriptional enhancers. Nat Genet. 2010;42:343–7.PubMed CentralPubMedView ArticleGoogle Scholar
- Wu JI, Lessard J, Crabtree GR. Understanding the words of chromatin regulation. Cell. 2009;136:200–6.PubMed CentralPubMedView ArticleGoogle Scholar
- Lee CS, Bhaduri A, Mah A, Johnson WL, Ungewickell A, Aros CJ, et al. Recurrent point mutations in the kinetochore gene KNSTRN in cutaneous squamous cell carcinoma. Nat Genet. 2014;46:1060–2.PubMed CentralPubMedView ArticleGoogle Scholar
- Adey A, Morrison HG, Asan, Xun X, Kitzman JO, Turner EH, et al. Rapid, low-input, low-bias construction of shotgun fragment libraries by high-density in vitro transposition. Genome Biol. 2010;11:R119.PubMed CentralPubMedView ArticleGoogle Scholar
- Boxer LD, Barajas B, Tao S, Zhang J, Khavari PA. ZNF750 interacts with KLF4 and RCOR1, KDM1A, and CTBP1/2 chromatin regulators to repress epidermal progenitor genes and induce differentiation genes. Genes Dev. 2014;28:2013–26.PubMed CentralPubMedView ArticleGoogle Scholar
- Kretz M, Webster DE, Flockhart RJ, Lee CS, Zehnder A, Lopez-Pajares V, et al. Suppression of progenitor differentiation requires the long noncoding RNA ANCR. Genes Dev. 2012;26:338–43.PubMed CentralPubMedView ArticleGoogle Scholar