The BASP1 transcriptional corepressor modifies chromatin through lipid-dependent and lipid-independent mechanisms

Summary The transcriptional corepressor BASP1 requires N-terminal myristoylation for its activity and functions through interactions with nuclear lipids. Here we determine the role of BASP1 lipidation in histone modification and the modulation of chromatin accessibility. We find that the removal of the active histone modifications H3K9ac and H3K4me3 by BASP1 requires the N-terminal myristoylation of BASP1. In contrast, the placement of the repressive histone modification, H3K27me3, by BASP1 does not require BASP1 lipidation. RNA-seq and ATAC-seq analysis finds that BASP1 regulates the activity of multiple transcription factors and induces extensive changes in chromatin accessibility. We find that ∼50% of BASP1 target genes show lipidation-dependent chromatin compaction and transcriptional repression. Our results suggest that BASP1 elicits both lipid-dependent and lipid-independent functions in histone modification and transcriptional repression. In accordance with this, we find that the tumor suppressor activity of BASP1 is also partially dependent on its myristoylation.

INTRODUCTION BASP1 was initially identified as a cytoplasmic signaling protein in neuronal cells (Mosevitsky, 2005) but has since been found to be widely expressed and localized to the nucleus in several cell types via a canonical nuclear localization sequence (Carpenter et al., 2004;Goodfellow et al., 2011;. A role for BASP1 in transcriptional regulation was first identified as a cofactor for the Wilms' tumor one protein WT1 (Carpenter et al., 2004). BASP1 binds to WT1 and converts it from a transcriptional activator to a repressor. BASP1 has subsequently been found to act as a repressive cofactor for other transcriptional regulators including MYC (Hartl et al., 2009), ERa (Marsh et al., 2017), and YY1 (Santiago et al., 2021), suggesting that it is broadly deployed in transcription control (Hartl and Schneider, 2019). Our understanding of how BASP1 regulates transcription is limited and has, to date, largely been studied as a transcriptional corepressor of WT1. BASP1 is N-terminally myristoylated and we demonstrated that this lipid motif is required for its function as a transcriptional repressor (Toska et al., 2012. BASP1 also contains a consensus cholesterol binding motif adjacent to the myristoyl moiety (Epand, 2008). Indeed, both phosphatidylinositol 4,5-bisphosphate (PIP 2 ) and cholesterol are recruited to promoter regions and are required for transcriptional repression by BASP1 (Toska et al., 2012Loats et al., 2021). The myristoylation of BASP1 is required for interaction with both PIP 2 and cholesterol (Toska et al., 2012;Loats et al., 2021). PIP 2 assists in the recruitment of HDAC1 to mediate the deacetylation of histone H3K9 and direct transcriptional repression (Toska et al., 2012).
Several studies have demonstrated that BASP1 acts as a tumor suppressor, for example in hepatocellular carcinoma (Tsunedomi et al., 2010), gastric cancer (Li et al., 2020), and breast cancer (Marsh et al., 2017). Consistent with its tumor suppressor activity, BASP1 has a role in maintaining the differentiated state and is required to maintain the functional differentiated state of taste receptor cells in mice (Gao et al., 2014(Gao et al., , 2019. BASP1 can also block the transformation of fibroblasts by v-myc (Hartl et al., 2009). Moreover, BASP1 interferes with the action of the Yamanaka proteins in the induction of pluripotent stem cells by repressing WT1 target genes (Blanchard et al., 2017). Taken together, current studies suggest that BASP1 plays a major role in the maintenance of the differentiated state through the modulation of the function of several transcriptional activator proteins.
The mechanism of action of BASP1 as a transcription cofactor is poorly understood, but its physiological activities suggest extensive roles in gene regulation. In this study, we set out to determine how BASP1 regulates transcription through histone modifications and the roles of the myristoyl motif of BASP1. We find that the lipidation of BASP1 plays a selective role in its action as a transcription cofactor and is required for specific histone modification leading to the modulation of the chromatin environment in order to elicit specific transcriptional effects. Our results suggest that BASP1 plays a widespread role in regulating the chromatin environment and transcription control.

RESULTS
BASP1 mediates the removal of active histone modifications and the placement of repressive H3K27me3 K562 cells express endogenous WT1 but do not express BASP1 (Goodfellow et al., 2011). As we have shown before, the introduction of BASP1 into K562 cells (B-K562 cells) leads to the transcriptional repression of the WT1 target genes AREG and VDR compared to the control K562 cells (V-K562 cells; Figure 1A; Goodfellow et al., 2011;Toska et al., 2012Loats et al., 2021). We next performed chromatin immunoprecipitation (ChIP) to detect WT1, BASP1, and a selection of specific histone modifications at the promoter regions of the WT1 target genes AREG and VDR ( Figure 1B). WT1 was present at the promoter region of the AREG and VDR genes in both control V-K562 cells and B-K562 cells, while BASP1 was present at the promoter region only in B-K562 cells. Analysis of H3K9ac and H3K4me3, indicators of transcriptionally active promoter regions, demonstrated that these histone modifications are removed when BASP1 is recruited to the promoter in B-K562 cells and the genes are transcriptionally repressed. These data are consistent with our previous studies of BASP1 function in histone modification (Toska et al., 2012;Loats et al., 2021). We also analyzed the H3K27me3 histone modification associated with transcriptional repression. We found that BASP1 induces the trimethylation of H3K27 at both the AREG and VDR promoter regions. Taken together, these results demonstrate that BASP1 functions by both the removal of active histone modifications (H3K9ac and H3K4me3) and the placement of the repressive histone mark H3K27me3.
We have previously used a conditional BASP1 mouse to demonstrate that BASP1 transcriptionally represses the WT1 target genes LEF1 and PTCH1 in taste receptor cells (Gao et al., 2014(Gao et al., , 2019. A floxed BASP1 mouse was crossed with a Krt8-Cre-ER mouse to delete BASP1 in the fully differentiated Krt8-expressing taste cells. We used this model to perform ChIP to analyze BASP1-dependent histone modifications at the LEF1 and PTCH1 promoters in taste receptor cells. Consistent with our previous studies, knockout of BASP1 expression in taste cells leads to the upregulation of LEF1 and PTCH1 expression ( Figure 1C). ChIP analysis confirmed that WT1 was present in the promoter region of the LEF1 and PTCH1 genes in both the control and BASP1 KO taste cells while BASP1 was only present in the promoter regions in the taste receptor cells from control mice ( Figure 1D). Analysis of H3K9ac, H3K4me3, and H3K27me3 demonstrated that loss of BASP1 leads to the accumulation of the active histone modifications H3K9ac and H3K4me3, and the removal of the H3K27me3 repressive modification. These findings, using an in-vivo model of BASP1 function, confirm and extend our analysis in K562 cells demonstrating that BASP1 directs the removal of the active histone modifications H3K9ac and H3K4me3 and the placement of repressive H3K27me3.

Myristoylated BASP1 is recruited to gene promoters
Our previous work demonstrated that the mutation BASP1-G2A, which prevents N-terminal myristoylation of BASP1, leads to loss of transcriptional repression function (Toska et al., 2012. Although total cellular BASP1 is stoichiometrically N-terminally myristoylated (Mosevitsky, 2005), it has not yet been demonstrated that this lipidated form of BASP1 is present within the nucleus and is directly recruited to gene promoter regions. We, therefore, used a click chemistry approach to analyze the myristoylation status of nuclear BASP1 in K562 cell line derivatives. V-K562 and B-K562 cells along with G-K562 cells (expressing BASP1-G2A) were incubated in lipid-free media supplemented with myristic acid alkyne that is utilized by cells as a substrate for N-terminal myristoylation (Wright et al., 2015). Nuclei were then prepared as we described before using a method that avoids contamination with either ER or Golgi complex (Loats et al., 2021). Using click chemistry, the alkyne-myristoyl moiety was crosslinked to azide-PEG3-biotin which was then detected using a streptavidin-linked fluorophore. BASP1 was simultaneously probed using immunocytochemistry. In control V-K562 cells, myristoyl was detected within the purified nuclei of K562 cells and was particularly evident at the nuclear membrane ( Figure 2A). This staining pattern is consistent with the known N-terminal myristoylation of lamins (Linde and Stick, 2010 iScience Article colocalization with myristoyl compared to wtBASP1 (Pearson quantitation of the proportion of BASP1 signal that overlaps with myristoyl is shown in Figure 2B).
We next used click chemistry combined with immunoprecipitation to determine the association between BASP1 and nuclear myristoyl. V-K562, B-K562, and G-K562 cells were cultured in the presence of myristic acid alkyne as above for 48 h and, following the click chemistry reaction with azide-PEG3-biotin, nuclear protein extracts were prepared then immunoprecipitation was performed with streptavidin-coated beads. wtBASP1, but not BASP1-G2A, co-immunoprecipitated with the alkyne-myristoyl ( Figure 2C). We have previously combined click chemistry with ChIP (Click-ChIP; Loats et al., 2021) and therefore used this technique A B C D Figure 1. BASP1 directs the removal of active histone modifications H3K9ac and H3K4me3 and the placement of repressive H3K27 trimethylation (A) Immunoblotting of extracts prepared from V-K562 and B-K562 cells to confirm BASP1 expression. b-tubulin immunoblotting was performed as a loading control. Molecular weight markers are shown at left (kDa). cDNA was prepared from V-K562 cells and B-K562 cells and expression of AREG and VDR was quantitated relative to GAPDH. Data are SD from the mean (SDM) for three independent experiments. * = p < 0.05 by students t-test. (B) V-K562 or B-K562 cells were subjected to chromatin immunoprecipitation (ChIP) with control (ctrl) antibodies or antibodies against WT1, BASP1 or the histone modifications indicated. Data are presented as fold-enrichment over a control genomic region and error bars are SDM of three independent experiments. * = p < 0.05 by students t-test.
(C) Krt8-BASP1-CRE mice were treated with tamoxifen for 8 days and then 7 days later taste buds were isolated and RNA prepared. cDNA was then used to monitor the expression of BASP1, LEF1, and PTCH1 compared to GAPDH. Error bars are SDM for three independent experiments. * = p < 0.05 by students t-test.
(D) Mice were treated as in part C but isolated taste cells were subjected to chromatin immunoprecipitation with the antibodies indicated. Error bars are SDM for three independent experiments. * = p < 0.05, ** = p < 0.01 and *** = p < 0.005 by students t-test.

OPEN ACCESS
iScience 25, 104796, August 19, 2022 3 iScience Article to determine if myristoyl is recruited to the promoter region of WT1 target genes. We found that myristoyl associates with the AREG and VDR promoters at a significantly increased level when wtBASP1 is also present at the promoter ( Figure 2D; compare V-K562 with B-K562), but at a reduced level in K562 cells that express BASP1-G2A (G-K562). Our previous study (confirmed in Figure 3B) demonstrated that BASP1 G2A is recruited to gene promoters at a level equivalent to wtBASP1 (Toska et al., 2012. These data demonstrate that the detection of myristoyl in the promoter region of WT1 target genes is dependent on BASP1 and that this association requires an intact G2 target site of myristoylation in BASP1.
Myristoylation of BASP1 is required for the removal of active chromatin modifications but not for the placement of repressive H3K27me3 Our data so far have demonstrated that myristoylated BASP1 is present at the promoter region of WT1 target genes. Our previous work (Toska et al., 2012, confirmed in Figures 3A and 3B, has shown A B C D Figure 2. Myristoylated BASP1 is present in the nucleus and is recruited to the promoter regions of WT1 target genes (A) The indicated cell lines were incubated in lipid-free media in the presence of 10 mg/mL myristic acid alkyne for 20 h. Nuclei were prepared and incubated with azide PEG3-biotin, then the click reaction was initiated using 2mM CuBF 4 . Immunohistochemistry was then performed with streptavidin-linked antibodies (Myristoyl) and BASP1 antibodies. Scale bar (blue) is 10mm. Control assays were performed that lacked azide PEG3-biotin (-azide control) or myristic acid alkyne (-alykyne control) and are shown at right. (B) Quantification of the colocalization of BASP1 and Myristoyl in B-K562 and G-K562 nuclei was analyzed using Pearson's correlation coefficient. *** = p < 0.005 following students t-test comparing Pearson's values from BK and GK nuclei over three independent experiments (n = 67).
(C) V-K562, B-K562, and G-K562 cells were incubated with ethanol (-alkyne control) or 10 mg/mL alkyne-Myristic acid as above and subjected to a click chemistry reaction. Nuclear extracts were prepared and precipitation was performed with streptavidin beads. Immunoblotting was then performed with either BASP1  iScience Article that BASP1-G2A lacks transcriptional corepressor activity despite being recruited along with WT1 to the promoter region of WT1 target genes. We demonstrated before that BASP1-G2A is defective in the deacetylation of H3K9 (Toska et al., 2012). We, therefore, tested if BASP1-G2A was able to mediate the demethylation of H3K4me3 or the placement of repressive H3K27me3. ChIP analysis confirmed that BASP1-G2A was defective in the removal of the active H3K9ac mark ( Figure 3C). We also found that BASP1-G2A was defective in the removal of H3K4me3. Thus, myristoylation of BASP1 is required for the removal of active histone modifications. In contrast, BASP1-G2A was still able to support the placement of repressive H3K27me3 at both the AREG and VDR promoters. Thus, the N-terminal myristoylation of BASP1 is required to mediate the removal of the active histone modifications, but is not required to place the repressive H3K27me3 modification. These findings suggest a dual mechanism of lipidation-dependent and lipidation-independent transcriptional repression by BASP1.

Myristoylation of BASP1 plays a role in the tumor suppressor activity of BASP1
BASP1 acts as a tumor suppressor in several cell types and slows the growth of K562 cells (Goodfellow et al., 2011;. We, therefore, tested the effects of wtBASP1 and BASP1-G2A on the formation of anchorage-independent colonies of K562 cells in soft agar assays ( Figure 3D). The number of V-K562, B-K562, and G-K562 colonies greater than 50 cells in size were counted over a three-week period ( Figure 3E). Expression of wtBASP1 led to a significant decrease in the colony-forming ability of K562 cells. Measurement of the average area of colonies revealed that, as well as forming fewer colonies, B-K562 cells form significantly smaller colonies ( Figure 3F). Expression of BASP1-G2A also led to a significant decrease in the colony-forming ability of K562 cells ( Figure 3E). However, the average area of G-K562 colonies did not significantly differ from that of V-K562 colonies ( Figure 3F). Thus, myristoylation of BASP1 is not required to suppress the formation of anchorage-independent colonies but contributes to inhibiting the growth of the established colonies. We conclude that the myristoylation of BASP1 plays a partial role in its tumor suppressor function in K562 cells.

Genome-wide analysis of the requirement for BASP1 myristoylation in transcriptional repression
Our results so far demonstrate that BASP1 directs histone modifications that are both dependent (H3K9 deacetylation and H3K4 demethylation) and independent (H3K27 methylation) on the N-terminal myristoylation of BASP1. In addition, myristoylation of BASP1 is partially required for its tumor suppressor activity. Taken together, these findings suggest that non-myristoylated BASP1 is likely to retain partial activity in transcriptional regulation that we have not observed here or before (Toska et al., 2012 in targeted gene analysis. We, therefore, determined the gene expression profiles of B-K562, G-K562, and control V-K562 cells using RNAseq. Comparing the data from V-K562 and B-K562 cells ( Figure 4A) supports our previous genome-wide analysis with 3064 genes regulated by BASP1 (2165 genes repressed and 899 genes activated; Table S1; Goodfellow et al., 2011). As before, the majority of the known WT1 target genes that show altered expression were transcriptionally repressed (70.3% of the 105 WT1 target genes that changed significantly (padj 0.05); Goodfellow et al., 2011;Figures 4B and 4C). Our data further show that BASP1-G2A is partially defective in transcriptional repression when compared to wtBASP1, repressing 58.8% of the genes that are repressed by wtBASP1 ( Figure 4A). This was reduced to 45.7% when considering only WT1 target genes ( Figures 4B and 4D). In fact, BASP1-G2A activated many of the target genes that are repressed by wtBASP1 ( Figure 4D).
Pathway enrichment analysis of genes repressed by wtBASP1 reveals several networks (canonical pathways and upstream regulators) that are enriched for significantly repressed differentially expressed genes in  iScience Article B-K562 cells, and further demonstrate their derepression in G-K562 cells (Figures 4E and 4F). We also identified other networks downregulated in B-K562, but not G-K562 cells, including those involved in the response to the H3K27 Demethylase Inhibitor GSK-J4, the selective estrogen receptor modulator Tamoxifen, the extra-cellular space, and pre-ribosome cellular processes ( Figure S1A). Together these data show the repressive effects of wtBASP1 are partially lost with the BASP1-G2A mutant derivative, specifically, the BASP1-dependent downregulation of transcription factors, upstream regulators, drug responses, and canonical pathways genome-wide ( Figures 4E, 4F, and S1A). BASP1 was also found to regulate the function of a number of transcription factors, in particular MYC, RELA, TFCP2, POU3F2, YY1, and CTCF ( Figure S1B; note that the expression level of the transcription factors was unchanged by wtBASP1 or BASP1 G2A compared to control cells). In the case of MYC, the effect of wtBASP1 was both activation and repression and a similar profile for MYC target genes was observed with K562 cells expressing BASP1-G2A. In contrast, wtBASP1 was predominantly repressive at RELA, TFCPP2, POU3F2, YY1, and CTCF target genes but BASP1-G2A elicited transcriptional activation at most genes that changed expression. Taken together the data in Figure 4 demonstrate that the myristoylation of BASP1 is required for its transcriptional repressor function at $50% of target genes. Furthermore, BASP1 regulates the expression of the target genes of both previously described BASP1-associated transcription factors (MYC, YY1, CTCF) and newly identified targets.

The effects of BASP1 and requirement for myristoylation on the chromatin landscape
Our targeted ChIP analysis of histone modifications (Figures 1 and 3) suggests that BASP1 likely makes extensive changes in the chromatin environment. Furthermore, the BASP1-G2A mutant derivative is defective in the removal of active chromatin modifications (H3K9ac and H3K4me3) but not the addition of repressive chromatin modification (H3K27me3). To investigate the effects of these BASP1-dependent chromatin remodeling events genome-wide, we conducted sequencing experiments to Assay for Transposase-Accessible Chromatin (ATAC-seq; Table S1). Comparing control V-K562 cells with B-K562 cells the data demonstrate that chromatin accessibility is substantially altered genome-wide by wtBASP1 with 79% of the peak changes in closed chromatin compared to 21% of the changes in open chromatin ( Figure 5A). These data are consistent with broad transcriptional repression by BASP1. 15% of the peak changes were within the upstream regions and transcription start sites of genes, but there were also general BASP1-dependent changes across the gene bodies. When we compared the ATAC-seq data with the RNA-seq data, we observed a strong correlation between the BASP1-dependent formation of closed chromatin within the gene proximity and transcriptional repression of the associated gene by BASP1 ( Figure 5B). Specifically, 514 genes that are transcriptionally repressed by BASP1 showed BASP1-dependent formation of closed chromatin peaks. Analysis of these genes by g:profiler highlighted several categories related to the physiological functions of BASP1 including developmental processes and cell differentiation ( Figure S2A). BASP1-G2A caused more peak changes in ATAC-Seq than wtBASP1 ( Figure S2B), although the ratio of changes in open chromatin versus changes in closed chromatin was increased ( Figure 4C). Direct comparison of wtBASP1 with BASP1-G2A showed that 66.3% of genes affected by wtBASP1 correlated with both decreased chromatin accessibility and increased transcriptional repression, while for BASP1-G2A, this only occurred for 45.5% of the genes ( Figure 4D). Indeed, BASP1-G2A was either completely ( Figure 4E, cluster 3) or partially ( Figure 4E cluster 4) defective in inducing closed regions of chromatin at most genomic regions modified by wtBASP1 (Non-relative data are shown in Figure S2C). Thus, although BASP1-G2A showed a defect in driving closed chromatin at many of the sites acted upon by wtBASP1, it did not show an overall lower activity in driving the formation of closed chromatin. This is consistent with BASP1-G2A acting at several sites that are not regulated by wtBASP1. Figure S3A shows a pathway analysis demonstrating that wtBASP1 and BASP1 G2A regulate many genes involved in distinct cellular processes. iScience Article Figure 5F shows the comparison of RNA-seq and ATAC-seq analysis across four WT1 target genes, AREG, VDR, EREG, and VEGFA. All four loci show a reduction of RNA-seq signal across the genes in B-K562 cells compared to V-K562 cells, which is coincident with a decrease in the chromatin accessibility around the transcription start sites of each gene (as measured by ATAC-Seq; Figure S3B shows an expanded view with surrounding genomic regions). In G-K562 cells, AREG, VDR, and EREG are not transcriptionally repressed (in-fact, BASP1-G2A activates the transcription of EREG), while VEGFA is still partially repressed by BASP1-G2A. Both EREG and AREG show that BASP1-G2A fails to induce closed ATAC peaks in the vicinity of the transcription start site which is coincident with the failure of BASP1-G2A to repress the transcription of these genes. BASP1-G2A produced an ATAC profile similar to wtBASP1 at the transcription start site region of the VDR gene but instead showed a reduction in ATAC peaks within the main body of the gene. Taken together, these data suggest that BASP1 transcriptionally represses WT1 target genes through the formation of closed chromatin. Furthermore, BASP1-G2A shows gene-specific effects in function, both in the regulation of chromatin compaction and in transcriptional repression.

DISCUSSION
BASP1 represses transcription through a novel mechanism that requires its lipidation to recruit PIP 2 and cholesterol to the promoter region of target genes. The localization of lipids, including PIP 2 and cholesterol, within the nucleus as well as their association with chromatin is well-documented, but their roles in nuclear processes are poorly understood (Garcia-Gil and Albi, 2017; Fernandes et al., 2018;Fiume et al., 2019;Barbosa and Siniossoglou, 2020;Gapa et al., 2022). Here we have presented evidence that myristoylated BASP1 is recruited to the promoter of WT1 target genes. We find that the N-terminal myristoylation of BASP1 controls its ability to modulate the chromatin environment and is needed for the full tumor suppressor activity of BASP1.
We have previously shown that BASP1interacts with PIP 2 to recruit HDAC1 to the gene promoter, which then deacetylates histone H3K9. The enzymes responsible for the BASP1-dependent demethylation of K3K4 and the trimethylation methylation of H3K27 have yet to be identified. It will be interesting to determine if BASP1 recruits the polycomb complex to trimethylate H3K27. BASP1-G2A retains the function to trimethylate H3K27 but is not able to mediate the deacetylation of H3K9 or demethylation of H3K4. This raises the possibility that BASP1-G2A leads to the formation of bivalent histone marks. Bivalent marks are associated with poised promoter regions that can either be activated or repressed (Voigt et al., 2013). The BASP1-G2A mutant derivative still causes significant changes in chromatin accessibility, suggesting that the H3K27 trimethylation function is sufficient, at least for a subset of target genes, to remodel chromatin and alter gene expression. Indeed, our RNAseq data demonstrate that BASP1-G2A retains transcriptional repressor activity at $46% of WT1 target genes and $59% overall. Genome-wide analysis of BASP1-and myristoylation-dependent histone modifications could be used to determine if BASP1-G2A leads to the formation of bivalent chromatin. Several previous studies have demonstrated the tumor suppressor activity of BASP1 (Hartl and Schneider, 2019). We found that BASP1-G2A is partially defective in tumor suppressor function which is consistent with the retention of chromatin modification and transcriptional regulatory function. It will be interesting to determine if BASP1 myristoylation plays a similar role in its tumor suppressor activity in other cell types.
The BASP1 G2A mutation prevents terminal myristoylation of BASP1 and abolishes its interaction with lipids. In that respect, the findings presented here raise the question of whether the lipid-binding function of wtBASP1 can be modulated to provide a regulatory mechanism for the transcriptional repression function of BASP1. Indeed, previous studies have demonstrated that the phosphorylation of Serine-6 within iScience Article BASP1 generates an unfavorable environment for the binding of lipids by BASP1 (Mosevitsky, 2005). We have previously reported that the phosphorylation of BASP1 disrupts its transcriptional repressor activity (Toska et al., 2012). This raises the possibility that the phosphorylation of BASP1 might regulate its ability to interact with different chromatin remodeling activities in the regulation of the chromatin landscape.
Our RNAseq data suggest that BASP1 can regulate the activity of several transcription factors. These include transcription factors previously identified as partners of BASP1, MYC (Hartl et al., 2009), YY1 (Santiago et al., 2021, and CTCF (Essafi et al., 2011). Several target genes of MYC and the associated MAZ protein were regulated by BASP1 which had both coactivator and corepressor function that was not dependent on the myristoylation of BASP1. This reduced dependence on BASP1 myristoylation might be owing to the proposed alternative mechanism whereby BASP1 interferes with the MYC-Calmodulin interaction. In contrast, YY1 and CTCF target genes showed BASP1-dependent repression that was substantially dependent on BASP1 myristoylation. RELA and TFAP2A/C target genes were also repressed by BASP1 which was dependent on BASP1 myristoylation. Taken together, these results suggest that BASP1 modulates the transcriptional activity of several transcription factors through both lipidation-dependent and lipidation-independent mechanisms.

Limitations of the study
In this study, we evaluated the effect of the disruption of the myristoyl motif of BASP1 on its function as a transcriptional coregulator. As the G2A mutation in BASP1 disrupts the interactions with both PIP 2 and cholesterol, the current study does not discriminate between the contribution of these two lipids to transcriptional regulation. Although we demonstrated that the BASP1-dependent effects on H3K9, H3K4, and H3K237 modification status occur in both K562 cells and in murine taste cells, the genome-wide study was limited to K562 cells. RNAseq and ATACseq in other cell types will be required to confirm the generality of these effects. Further studies to determine the genome-wide localization of BASP1 will be required to distinguish between direct and indirect effects of BASP1 or BASP1 G2A in the RNAseq and ATACseq datasets. Finally, the gene distance cut-off of 10kb that was applied to the analysis of potential regulatory elements is an accepted standard for such analyses, providing a suitable balance of sensitivity and specificity for the detection of cis-acting regulatory effects. However, we recognize that this approach excludes distant regulatory elements.

STAR+METHODS
Detailed methods are provided in the online version of this paper and include the following:

Materials availability
This study did not generate unique reagents.

Data and code availability
All data reported in this paper will be shared by the lead contact upon request. The RNA-seq and ATACseq data are available at the NCBI Gene Expression Omnibus GSE196525. This paper does not report original code. Any additional information required to reanalyze the data reported in this paper is available from the lead contact upon request.
Stably transfected K562 cell line derivatives selected as a pool (Goodfellow et al., 2011;Toska et al., 2012) were supplemented with 1 mg/mL G-418 (Sigma-Aldrich). All cell line derivatives were screened by immunofluorescence to ensure that BASP1 was expressed in >95% of the cells. Cells were incubated at 37 C in humidified 95% air and 5% CO 2 . Lipid-free media was prepared as above but using charcoalstripped fetal calf serum (Thermofisher).

Chromatin immunoprecipitation
K562 cells were harvested and resuspended to a density of 1 3 10 6 cells/mL in PBS. Crosslinking was initiated by the addition of formaldehyde to a final concentration of 1.42% (v/v) and incubation at room temperature for 15 min. Crosslinking was terminated by the addition of glycine to a final concentration of 125mM glycine and incubation at 4 C for 5 min. The cells were then harvested by centrifugation at 2000 x g for 5 min at 4 C and then washed with PBS. Cells were resuspended and lysed in 1mL IP buffer (150mM NaCl, 50mM Tris-HCl (pH 7.5), 5mM EDTA, 0.5% (v/v) NP-40 and 1% (v/v) Triton X-100) with protease inhibitor cocktail for 15 min on ice. After centrifugation at 2000 x g for 5 min at 4 C, the pellet resuspended in 1mL IP buffer. The chromatin was sheared by sonication using a QSonica Q500 at 60% amplitude and the lysate was cleared by centrifugation at 12000 x g for 10 min at 4 C. The samples were pre-cleared by iScience 25, 104796, August 19, 2022 iScience Article incubation for 1 h at 4 C on a rotator with 10mL Protein G magnetic beads (ThermoFisher). The indicated antibody was incubated in 600mL IP buffer containing 1mL 10 mg/mL acetylated BSA and 5mL Protein G magnetic beads for 4 h with rotation at 4 C. 200mL of the pre-cleared chromatin was added and incubation continued at 4 C overnight. A sample of the input chromatin was stored for later decrosslinking and processing.

qPCR analysis
Total RNA was prepared using the RNeasy kit (Qiagen) and then cDNA prepared using the Iscript cDNA synthesis kit (Bio-Rad). qPCR was performed in triplicate using the iQ SYBR green master mix (Bio-Rad) on a Bio-Rad Miniopticon system. ChIP analysis was performed using promoter-specific primers as indicated and fold-enrichment calculated against a control genomic region. Melt curve analysis was performed at the end of each run. Data were collected using the BioRad-CFX Manager software. The relative expression of each target gene is expressed relative to control GAPDH. Data are expressed as mean with SD(SDM). Data distribution and significance between different groups was analyzed in Excel or OriginPro 7.5 using Student's t test.
Immunofluorescence K562 cell line derivatives were collected by centrifugation at 1400 x g for 3 min and nuclei were isolated using the Nuclei EZ Prep nuclei isolation kit (Sigma-Aldrich). The nuclei were fixed in 4% (v/v) paraformaldehyde at room temperature for 15 min, then incubated with 50mM NH 4 Cl for 15 min with rotation. The fixed nuclei were then washed three times in PBS and then incubated in blocking buffer (2% (w/v) BSA, 0.25% (w/v) Gelatin, 0.2% (w/v) Glycine, 0.2% (v/v) Triton X-100 in PBS) for 1 h with rotation at room temperature. Primary antibody was added to each sample of nuclei in PBS containing 1% (w/v) BSA, 0.25% (w/v) gelatin and 0.2% (v/v) Triton X-100 and then incubated for 1 h with rotation at room temperature. The nuclei were then washed three times with washing buffer (0.2% (w/v) Gelatin in PBS). Fluorescent secondary antibody in PBS containing 1% (w/v) BSA, 0.25% (w/v) gelatin and 0.2% (v/v) Triton X-100 was then incubated with the nuclei for 45 min in with rotation in the dark. After washing three times the samples were incubated with DAP solution for 10 min and then resuspended in a minimum volume of DABCO mounting media (Sigma-Aldrich) an applied to poly-lysine coated slides. Nuclei were viewed using a Leica SP5-II AOBS confocal laser scanning microscope attached to a Leica DM I6000 inverted epifluorescence microscope with oil 63x lens. Images were processed using ImageJ or Volocity 6.3 software.

Click chemistry
K-562 cells were incubated in lipid free media containing 10 mg/mL myristic acid-alkyne (Avanti) for 16 h. Cells were then harvested and washed once in PBS and the nuclei isolated using the EZ Prep nuclei kit. The purified nuclei were then incubated with 50mM Azide-PEG3-biotin conjugate (Sigma-Aldrich) in Buffer A and added to samples. The click reaction was initiated via addition of 2mM CuBF4 in 2% (v/v) acetonitrile, and the reaction was left to proceed at 43 C for 30 min with rotation. The nuclei were then washed in 0.1M HEPES/KOH (pH 7.4) and immediately used for either preparation of nuclear extract for immunoprecipitation, chromatin immunoprecipitation or immunofluorescence.

Immunoprecipitation
Following completion of the click chemistry reaction, nuclear extracts were prepared by resuspending samples in one PCV of NE2 buffer (20 mM HEPES pH 8.0, 1.5 mM MgCl2, 25% (v/v) Glycerol, 420 mM NaCl, 0.2 mM EDTA, 1 mM DTT and 0.5 mM PMSF) and incubating on ice with regular stirring for 30 min. Nuclear debris was pelleted by 5 min in a microfuge at full speed. The cleared nuclear extracts were then subjected to immunoprecipitation using streptavidin magnetic beads (Invitrogen) and IP buffer (

Colony formation assay
The K562 cell line derivatives were seeded onto agar bases (0.7% (w/v) LMP Agarose, 10% (v/v) FBS, 1% (v/v) Pen Strep, 1% (v/v) L-glutamine in DMEM media supplemented with 0.15% (v/v) NaHCO3, 200 mM sodium pyruvate and 800 mM NaOH at a density of 2 3 105 cells per mL of agar. Agar dishes were fed with additional agar (2 mL) on days 7 and 14 post seeding. The colony formation efficiency (%CFE) was calculated 7, 14 and 21 days post seeding via counting colonies consisting of greater than 50 cells, in a minimum of ten 15.89 mm 2 fields of view per 60 mm 2 dish. Average colony area was measured using ImageJ software. All experiments were performed in triplicate.

Animals
The floxed BASP1 mice in a C57BL/6 background have been described before (Gao et al., 2019). We have not detected any sex differences related to the effects of BASP1 loss. However, we routinely kept records of both sex and age (2-5 months) to ensure no bias is introduced. Animals were cared for in compliance with the University at Buffalo Animal Care and Use Committee. BASP1fl+/À mice were mated with Krt8-Cre/ ERT2+ mice (The Jackson Laboratory) to obtain BASP1fl+/ÀKrt8-Cre/ERT2+ mice and then bred to obtain BASP1fl+/+Krt8-Cre/ERT2+ mice. BASP1flÀ/À; Krt8-Cre/ERT2+ or BASP1fl+/+; Krt8-Cre/ERT2À littermates were used as wild-type controls. Adult mice were gavaged with 100 mg/kg body weight of Tamoxifen dissolved in corn oil/ethanol (9:1 by volume) daily for 8 days 7 days later the mice were euthanized and taste receptor cells harvested from the circumvallate papillae as described before (Dutta Banik et al., 2018). The mice were euthanized by carbon dioxide followed by cervical dislocation and the tongues removed. A solution containing 0.7 mg/mL Collagenase B (Roche), 3 mg/mL Dispase II (Roche), and 1 mg/mL trypsin inhibitor (Sigma-Aldrich) was then injected beneath the lingual epithelium. The tongues were then incubated in in oxygenated Tyrode's solution for 15 min and then the epithelium peeled from the muscle and placed in Ca2+-free Tyrode's. The cells were used immediately for either RNA preparation or ChIP.

Sample preparation and sequencing
All samples were prepared in biological triplicates for V-K562, B-K562 and G-K562 cells harvested at 0.7-0.8 3 10 6 cells/mL (Loats et al., 2021;Toska et al., 2012). Cells taken for ATAC-seq were treated with DNase and incubated for 30 mins buffered with 25mM MgCl2 and 5mM CaCl2 prior to harvesting. Samples of 2 3 10 6 for RNA-seq and of 1 3 10 5 cells for ATAC-seq were then pelleted and snap-frozen in liquid nitrogen until library preparation. Total RNA was isolated from the cells using the RNeasy Mini Kit (Qiagen, cat 74104), RNA Integrity Number (RIN) was determined for each sample using Tape Station 4150 System (Agilent). Libraries were prepared according to manufacturer's instructions for TruSeq Stranded mRNA Library kit (Illumina, cat #20020594) and ATAC-Seq Kit (Active Motif, cat #53150). Libraries were sequenced on Illumina NextSeq 500 in 75-nt (RNA-seq) and 42-nt (ATAC-seq) experiments in paired-end mode.

RNA-seq and ATAC-seq data analysis
Reads were mapped to GRCh38 human reference genome using STAR (Dobin et al., 2013) (RNA-seq) and BWA-MEM (Li and Durbin, 2009) (ATAC-seq), duplicate reads were removed, uniquely mapping reads with %2 mismatches were kept for downstream analysis. RNA-seq paired reads were assigned to RefSeq transcripts using Subread and annotated with Entrez gene identifiers. ATAC peaks were called using MACS2, merged regions were defined by overlapping intervals between all samples, peak fragment densities for merged regions were used for comparative analysis and annotated with the closest gene within a 10kb window. DESeq2 was used for differential expression analysis to determine p-adjusted values (padj) and shrunken Log2 fold-change (FC) for RNA-seq genes and ATAC-seq merged regions. (Gaspar, 2018;Liao et al., 2013;Love et al., 2014) The 5000 most differentially expressed genes for BK vs VK at padj %0.05 were analyzed using iPathway to determine significantly differentially regulated networks. The 'Impact Analysis' was applied to underlying pathway topologies comprised of genes and their directional interactions. The probability of observing the number of DE genes in a given pathway that is R to that observed by random chance (pORA) and perturbation accumulation, based on the combined perturbation of all genes within the pathway as a function of normalised change in expression and weighting in relation to gene type, (pAcc), are the principal measures ll OPEN ACCESS iScience Article