Abstract Background SWI/SNF chromatin remodeling genes are required for normal acute responses to alcohol in C. elegans and are associated with alcohol use disorder in two human populations. In an effort to discover the downstream genes that are mediating this effect, we identified SWI/SNF-regulated genes in C. elegans. Results To identify SWI/SNF-regulated genes in adults, we compared mRNA expression in wild type and swsn-1(os22ts) worms under conditions that produce inactive swsn-1 in mature cells. To identify SWI/SNF-regulated genes in neurons, we compared gene expression in swsn-9(ok1354) null mutant worms that harbor a neuronal rescue or a control construct. RNA sequencing was performed to an average depth of 25 million reads per sample using 50-base, paired-end reads. We found that 6813 transcripts were significantly differentially expressed between swsn-1(os22ts) mutants and wild-type worms and 2412 transcripts were significantly differentially expressed between swsn-9(ok1354) mutants and swsn-9(ok1354) mutants with neuronal rescue. We examined the intersection between these two datasets and identified 603 genes that were differentially expressed in the same direction in both comparisons; we defined these as SWI/SNF-regulated genes in neurons and in adults. Among the differentially expressed genes was cbp-1, a C. elegans homolog of the mammalian CBP/p300 family of histone acetyltransferases. CBP has been implicated in the epigenetic regulation in response to alcohol in animal models and a polymorphism in the human CBP gene, CREBBP, has been associated with alcohol-related phenotypes. We found that cbp-1 is required for the development of acute functional tolerance to alcohol in C. elegans. Conclusions We identified 603 transcripts that were regulated by two different SWI/SNF complex subunits in adults and in neurons. The SWI/SNF-regulated genes were highly enriched for genes involved in membrane rafts, suggesting an important role for this membrane microdomain in the acute alcohol response. Among the differentially expressed genes was cbp-1; CBP-1 homologs have been implicated in alcohol responses across phyla and we found that C. elegans cbp-1 was required for the acute alcohol response in worms. Keywords: SWI/SNF chromatin remodeling, Ethanol, Alcohol, C. elegans, Transcriptome, cbp-1, Histone acetyltransferase, Differential gene expression Background Alcohol is a widely used and abused drug. The significant majority of adults in the United States uses alcohol, and the misuse of alcohol has serious societal impacts [[35]1]. Disturbingly, problems associated with alcohol misuse appear to be getting worse; for example, the rate of deaths associated with alcohol use has risen substantially in the last two decades, particularly among groups that have traditionally lagged in alcohol-related mortality [[36]2]. The liability to abuse alcohol is strongly influenced by an individual’s genetics [[37]3], although the specific genes involved are as yet incompletely understood. An individual’s level of response (LR) to alcohol manifests earlier than alcohol problems, and is a strong predictor of the likelihood to develop alcohol use disorder (AUD) [[38]4–[39]6]. LR to alcohol is a composite assessment that consists of both physiological measures, such as body sway, as well as subjective measures of intoxication. A lower sensitivity to the depressing effects of alcohol early in adulthood is correlated with a higher lifetime incidence of alcohol problems [[40]7–[41]9]. Identifying the genes that modulate the LR to alcohol is the goal of our work. We use the nematode Caenorhabditis elegans as a model for understanding the acute physiological response to alcohol (ethanol), and for identifying genes that are likely to affect LR in humans. We previously found that SWI/SNF chromatin remodeling complex genes were required for normal acute ethanol response behaviors in C. elegans [[42]10]. SWI/SNF are multi-protein complexes that modify chromatin structure by regulating the interaction of histones and DNA at particular genes, occluding or revealing transcription factor binding sites (reviewed in [[43]11]); these actions result in an epigenetic modulation of target gene expression. SWI/SNF complexes share common enzymatic core proteins, and different complexes have different accessory subunits that confer their sequence specificity (reviewed in [[44]12, [45]13]). In worms, most SWI/SNF complex members were involved in regulating acute ethanol response behaviors, but different SWI/SNF complexes affected different behaviors: BAF subunits were required for normal initial sensitivity to ethanol, while PBAF subunits were required for acute functional tolerance to ethanol [[46]10]. We found that allelic variation in several SWI/SNF complex members was associated with alcohol dependence in different adult human populations [[47]10, [48]14]. Allelic variation in different SWI/SNF complex members was also associated with antisocial behavior in adolescents, a phenotype that is strongly predictive of future alcohol use problems [[49]14]. Together, these data strongly argue that the SWI/SNF complex plays a functional role in human AUD, and support the investigation of the biological underpinnings of this role. Because SWI/SNF acts to regulate transcription, the downstream effectors through which it modulates the acute ethanol response are of significant interest. Here, we use transcriptional profiling in C. elegans to identify these downstream modulators. We took advantage of two genetic reagents, one that localizes expression of swsn-9 to neurons using tissue-specific rescue of a swsn-9 null allele, and one that temporally restricts function of swsn-1 to adults using a temperature-sensitive allele of swsn-1 [[50]15]. swsn-9 encodes a PBAF subunit, while swsn-1 encodes a core subunit that is predicted to be present in both BAF and PBAF; both genes are required for acute functional tolerance to ethanol [[51]10]. We found that 603 genes were regulated by the SWI/SNF complex in adults and in neurons. Among these genes was cbp-1, encoding a histone acetyltransferase. We found that animals carrying loss of function alleles of cbp-1 had defects in acute ethanol response behaviors, indicating that cbp-1 is an important SWI/SNF target for LR to ethanol. Results SWI/SNF-regulated genes in adults We previously showed that swsn-1 is required during adulthood for the development of acute functional tolerance (AFT) to ethanol [[52]10]. swsn-1(os22ts) is a temperature sensitive allele of swsn-1; os22ts is an amino acid substitution that causes the protein encoded by swsn-1(os22ts) to be functional at 15 °C, but destabilized so that its function is lost at 25 °C [[53]15]. For these experiments, we reared N2 and swsn-1(os22ts) at the permissive temperature (15 °C) from the time the eggs were laid through the fourth larval stage (L4), at which time we shifted the worms to the restrictive temperature (25 °C) and grew them to adulthood (Fig. [54]1a). In wild type worms, SWSN-1 is expressed both during development and in adulthood. In swsn-1(os22ts) mutant worms, this temperature shift protocol produces functional SWSN-1 during development and inactive SWSN-1 in adulthood. We found that worms lacking SWSN-1 during adulthood failed to develop AFT, while wild type worms reared under the same conditions developed AFT to ethanol [[55]10]. This result suggests that the SWI/SNF complex regulates gene expression in adults that is necessary for AFT in these animals. To identify genes that are regulated by SWSN-1 in adults, we performed RNA sequencing on wild type and swsn-1(os22ts) mutant worms that were subjected to the temperature shift protocol. Five biological replicates were generated on different days; mutant and wild type worms were reared in parallel, thus limiting the effect of environmental variation on our comparison. We isolated RNA from young adult worms approximately 6 h after the temperature shift, and at a time when they were beginning to lay eggs. RNA sequencing libraries were prepared and sequenced by the Genomic Services Laboratory at Hudson Alpha (Huntsville, AL). Fig. 1. [56]Fig. 1 [57]Open in a new tab SWI/SNF-regulated genes in adults. a A temperature-sensitive allele, swsn-1(os22ts), was used to demonstrate that SWI/SNF complexes are required in adult tissues for AFT [[58]10]. Diagram (top) shows the temperature shift regimen and table (bottom) summarizes alcohol response behaviors from this experiment. b Principal component analysis with gene expression plotted relative to the first two principal components (PC1 and PC2). Wild type and swsn-1(os22ts) replicates are most similar to replicates of the same genotype. c Volcano plot shows genes that are up-regulated (red) and down-regulated (blue) in swsn-1(os22ts) mutants. Dashed lines indicate the FDR and fold change cutoffs (FDR ≤ 0.05 and fold change ≥1) We assessed the correlation between biological replicates using principal component analysis and found that our biological replicates grouped together and the two genotypes cluster in distinct groups (Fig. [59]1b). The first two principal components accounted for 95% of the variance in the dataset, with principal component one (variation between sample types) accounting for 92% of the variance. We analyzed differential gene expression using DESeq2 [[60]16] and found that 6813 genes were differentially expressed between swsn-1(os22ts) mutant and wild type worms (FDR ≤ 0.05) (Additional file [61]1). A volcano plot shows the wide distribution of differentially expressed genes (DEGs) (Fig. [62]1c). Approximately equal numbers of genes had increased (3210) and decreased (3603) expression in swsn-1(os22ts) mutants when compared to wild type. This suggests that the SWI/SNF complex both activates and represses gene expression in adult tissues, which is consistent with previous observations of SWI/SNF regulatory activity [[63]17–[64]20]. SWI/SNF-regulated genes in neurons We previously showed that swsn-9(ok1354) mutant worms were unable to develop AFT to ethanol and that, if we expressed swsn-9::GFP in all neurons under the control of the ric-19 promoter, we could restore AFT to swsn-9(ok1354) mutants [[65]10]. This result indicates that the SWI/SNF complex regulates gene expression in neurons that is important for AFT. Our previous experiments used extra-chromosomal arrays to show that swsn-9 was required in neurons. Extra-chromosomal arrays are multicopy DNA concatamers that are transmitted with varying frequency through mitosis and meiosis [[66]21]. In order to have uniform transgene expression and normal diploid gene dosage, we generated single copy insertions of the neuronal expression construct Pric-19::swsn-9::GFP using the Mos-1 mediated single copy insertion (MosSCI) technique [[67]22]. The MosSCI targeting vector contains a selectable marker; therefore we also generated single copy insertions of the empty vector as a control (Fig. [68]2a). We crossed these insertions into swsn-9(ok1354) and tested them for rescue of the swsn-9(ok1354) AFT defect (Fig. [69]2b). Wild type worms develop AFT, while swsn-9(ok1354) worms fail to develop AFT. We found that single copy insertions of the empty vector did not rescue AFT, while single copy insertions containing Pric-19::swsn-9::GFP did rescue the AFT defect of swsn-9(ok1354) null mutants. Fig. 2. [70]Fig. 2 [71]Open in a new tab SWI/SNF-regulated genes in neurons. a Constructs for single-copy rescue of swsn-9(ok1354) mutants in neurons. The empty MosSCI vector (top) serves as a negative control and contains the unc-119+ selectable marker. The neuronal rescue plasmid (bottom) expresses SWSN-9::GFP under control of the ric-19 promoter; it also contains the unc-119+ selectable marker. b Animals were tested on 0 mM and 400 mM exogenous ethanol. Relative speeds were calculated as treated over untreated speeds. The development of AFT is indicated by a statistically significant recovery of speed between 10 and 30 min of exposure. All genotypes on the same graph were tested simultaneously on the same plates. Wild-type (WT) animals develop AFT, while swsn-9(ok1354) mutants fail to develop AFT (n = 6). A single copy insertion of Pric-19::SWSN-9::GFP (rescue) was able to rescue the AFT defect of swsn-9(ok1354) mutants, while a single copy insertion of the empty MosSCI vector (control) did not rescue the AFT defect. Error bars indicate S.E.M. Paired one-tailed Student’s t tests were used for statistical comparisons of speeds at 10 and 30 min. One-way ANOVA, with Tukey’s multiple comparison test, was used to compare initial sensitivity across the strains; no significant difference was observed. For indicated comparisons: ns, not significant; **p ≤ 0.01 c Principal component analysis with gene expression plotted relative to the first two principal components (PC1 and PC2). swsn-9(ok1354) mutant replicates containing the rescue or control constructs are most similar to replicates of the same genotype. d Volcano plot shows genes that are up-regulated (red) and down-regulated (blue) in swsn-9(ok1354) mutant neurons. Dashed lines indicate the FDR and fold change cutoffs (FDR ≤ 0.05 and fold change ≥1) In order to identify SWI/SNF-regulated neuronal genes, we performed RNA sequencing on swsn-9(ok1354) mutants harboring either the empty vector control or the neuronal swsn-9::GFP rescue construct. We collected five biological replicates on different days and rescue and control worms were reared in parallel. We isolated total RNA from young adult worms of each genotype. RNA sequencing libraries were prepared and sequenced by the Genomic Services Laboratory at Hudson Alpha (Huntsville, AL). We assessed the correlation between biological replicates using principal component analysis and found that our biological replicates grouped together (Fig. [72]2c). The first two principal components accounted for 75% of the variance in the dataset, with principal component one (variation between sample types) accounting for 57% of the variance. We analyzed differential gene expression using DESeq2 [[73]16] and found that 2412 genes were differentially expressed between swsn-9(ok1354) mutants with and without neuronal rescue (FDR ≤ 0.05) (Additional file [74]1). A volcano plot shows the distribution of DEGs (Fig. [75]2d). Consistent with our swsn-1(os22ts) data, we found that similar numbers of genes were down-regulated (1157) as were up-regulated (1255). SWI/SNF-regulated genes, in neurons, and in adults The intersection of these two datasets identifies genes that are differentially expressed in neurons and in adults, in SWI/SNF mutants relative to wild type. Because we isolated mRNA from whole adult animals and we expect that the relevant gene expression changes will be in a small number of cells, we have not used any fold-change cutoff for this analysis. We identified 603 genes that were differentially expressed in the same direction in both swsn-1(os22ts) and swsn-9(ok1354) mutants compared to wild type (FDR ≤ 0.05); 393 genes were up-regulated and 210 were down-regulated (Fig. [76]3a). These 603 genes are excellent candidates for mediators of the SWI/SNF effect on AFT. Fig. 3. [77]Fig. 3 [78]Open in a new tab SWI/SNF-regulated genes in neurons and adults. a Differential gene expression analysis identified 6813 genes with differential expression between swsn-1(os22ts) and wild type adult worms (purple) and 2412 genes between swsn-9(ok1354) mutants with and without neuronal rescue (blue). Of the 603 DEGs that were regulated in the same direction in both comparisons, 392 were up-regulated, and 210 were down-regulated. b PANTHER GO biological process and cellular compartment terms enriched in this gene set; there were no significant GO molecular function terms. GO terms are plotted against the fold enrichment relative to the expected number of genes in lists of these sizes GO term enrichment analysis We performed an analysis to determine if particular functional classes of genes are overrepresented in our DEGs. We found that gene ontology (GO) biological process terms associated with metabolic processes, such as “mitochondrial ATP synthesis coupled electron transport” (GO:0042775) and “carboxylic acid catabolic process” (GO:0046395), were overrepresented in the swsn-1(os22ts) DEGs (Additional file [79]2). By contrast, we found a wide range of overrepresented GO biological process terms for the swsn-9(ok1354) DEGs, ranging from “regulation of cell shape” (GO:0008360) to “translation initiation” (GO:0006413) to “mitochondrial gene expression” (GO:0140053) (Additional file [80]2). Interestingly, we found that the GO biological process term “fatty acid metabolic process” (GO:0006631) was overrepresented in both the swsn-1(os22ts) and swsn-9(ok1354) DEGs, with 2.4 fold enrichment in the swsn-9(ok1354) DEG set and 1.8 fold enrichment in the swsn-1(os22ts) DEGs. Finally, we performed gene ontology GO term overrepresentation tests for genes at the intersection of our two datasets (Additional file [81]2). We found 7.5 times as many genes with the GO cellular process term “membrane rafts” (GO:0045121) and 3.4 times as many genes with the GO biological process term “innate immune response” (GO:0045087) as would be expected for a gene list of this size (Fig. [82]3b). Comparison of swsn-1 and swsn-9 regulated genes SWI/SNF complexes are composed of multiple subunits that combine to form functionally distinct complexes [reviewed in [[83]12, [84]13]]. Our two RNA sequencing datasets utilized different SWI/SNF mutants: swsn-1 encodes a core subunit that is present in all SWI/SNF complexes, while swsn-9 encodes an accessory subunit that is specific to the PBAF complex in worms (Fig. [85]4a) [[86]10, [87]23]. We have an opportunity with these datasets to compare the effect of loss of different SWI/SNF subunits on gene expression. We examined the 2000 most variable genes in wild type, swsn-1(os22ts), and swsn-9(ok1354) using k-means clustering (Fig. [88]4b-d; Additional file [89]3). As expected, we identified clusters containing genes that were up- or down-regulated in swsn-1(os22ts) (clusters A, D, and E) and swsn-9(ok1354) mutants (clusters B and G). Note that clusters B and G contain genes that are regulated by swsn-9 and also genes that respond to temperature, since the wild type and swsn-1(os22ts) worms were shifted from 15 °C to 25 °C prior to gene expression analysis, whereas swsn-9(ok1354) mutants were not. We also identified clusters of genes that are regulated similarly in both swsn-1(os22ts) and swsn-9(ok1354) mutants (clusters C and F). These genes are regulated by the PBAF complex in adult neurons. Cluster C contains genes that are normally activated, while Cluster F contains genes that are normally repressed by PBAF. Fig. 4. [90]Fig. 4 [91]Open in a new tab Comparison of swsn-1 and swsn-9 gene expression. a Diagram of SWI/SNF complexes in C. elegans. Genetic studies support the presence of distinct BAF and PBAF complexes that form by the association of unique subunits with the core complex. The SWSN-1 (blue) and SWSN-9 (purple) proteins are indicated. b-d k-means clustering of the top 2000 most variable genes across three genotypes (five replicates per genotype). b t-SNE map shows distinct clustering for seven k-means clusters. c Description of the seven clusters. d Heatmap of clustered genes; genotypes are indicated on the bottom We asked if our PBAF target genes were regulated by PBAF in other systems. Chowdhury et al. identified 2464 genes that were differentially expressed in a renal cell carcinoma model upon re-expression of PBRM-1 [[92]24]. We identified human homologs of the genes in clusters C and F using the DRSC Integrative Ortholog Prediction Tool (DIOPT) [[93]25], retaining any genes with a DIOPT score greater than two. This resulted in lists of 336 and 312 unique human homologs of the genes in clusters C and F, respectively. We compared these lists of human genes with the DEGs found in the Chowdhury study and found that 11.2% of the cluster C and F homologs were regulated by PBRM-1 in cancer cell lines (Additional file [94]3). Human homologs of SWI/SNF-regulated genes are associated with alcohol-related phenotypes We are interested in finding SWI/SNF-regulated genes that are mediating the effects of the SWI/SNF complex on human alcohol-related phenotypes. We identified human homologs of the 603 worm genes identified at the intersection of our two datasets using DIOPT [[95]25]. This resulted in a list of 874 unique human genes. We generated a list of human genes for which there was a reported gene disease association in the literature using DisGenNet [[96]26]. We included disease terms related to alcohol or substance misuse, consumption, and withdrawal. Among the 874 human homologs of our SWI/SNF-regulated genes we found that 84 (9.6%) are linked to one or more alcohol or substance use phenotypes in the DisGeNET database (Additional file [97]4). Next, we used the human gene list to query PubMed for publications linking these genes to problematic alcohol use (search terms include “alcoholism”, “alcohol use disorder”, or “alcohol dependence”) or to alcohol more broadly (“alcohol” or “ethanol”). We found 418 genes with one or more alcohol publications and 224 genes with one or more alcohol phenotype publications (Additional file [98]4). We sought to identify epigenetic regulators among our SWI/SNF-regulated genes because epigenetic regulators, like SWI/SNF, often work together. We queried PubMed for any citations linking these genes to epigenetic regulation and found 518 genes with one or more PubMed citations containing the search term “epigenetic”. Finally, we generated a list of 195 genes for which there is literature evidence for both alcohol misuse phenotypes and epigenetic mechanisms of action (Additional file [99]4). These 195 genes were ranked based on the number of citations linking the genes to epigenetics; the top twenty genes are reported in Table [100]1. These are our top candidate genes for future study. Table 1. Top 20 human genes linked to alcohol problems and epigenetics Human gene PubMed Citations Worm Gene Expression Alcohol problems Epigenetics Gene Expression STAT3 87 363 sta-1 Up GSTP1 19 284 gst-16 Down gst-39 Up gst-38 Up gst-19 Down EP300 1 153 cbp-1 Down NOTCH1 7 149 glp-1 Down SLC6A4 193 123 snf-11 Up CREBBP 6 123 cbp-1 Down STAT1 26 113 sta-1 Up TLR4 202 102 tol-1 Down SFN 2 96 ftt-2 Up VDR 7 95 nhr-8 Up RHOA 7 62 rho-1 Down TLR2 42 58 tol-1 Down FRY 19 58 sax-2 Down ABCB1 6 58 pgp-5 Up AGO2 2 58 alg-2 Down alg-5 Down HNF4A 9 56 nhr-35 Up nhr-60 Down STAT4 4 45 sta-1 Up RAC1 6 42 ced-10 Down LMNA 2 41 ifd-2 Down ifd-1 Up KAT5 1 40 mys-4 Down [101]Open in a new tab The histone acetyltransferase CBP-1 is dose-dependently required for AFT Among our top candidate genes is the C. elegans gene encoding CBP, cbp-1, which is similarly homologous to two closely related human genes, CREBBP and EP300. Both CREBBP and EP300 appear in our top 20 genes related to human alcohol problems and epigenetics (Table [102]1), strongly supporting our interest in cpb-1. CBP homologs have well described roles in the nervous system response to alcohol in animal models [[103]27–[104]31], and CREBBP is associated with addiction in humans [[105]32]. We reasoned that if SWI/SNF is regulating ethanol response behaviors in part through its regulation of cbp-1, then disrupting function of cbp-1 itself should alter ethanol response behaviors. cbp-1 expression was decreased relative to wild-type in swsn-1(os22ts) and swsn-9(ok1354) mutants (Fig. [106]5b), both of which are defective in acute functional tolerance (AFT), so we predicted that a loss of cbp-1 should decrease or eliminate AFT. We therefore tested two different strong loss-of-function alleles of cbp-1 in our behavioral assays (Fig. [107]5a). Mutants carrying homozygous cbp-1 loss-of-function alleles die, so we tested cbp-1 heterozygotes. We found that cbp-1(ok1491) heterozygotes did not develop AFT (Fig. [108]5c). We found that cbp-1(bm2)/balancer heterozygotes had reduced AFT, whereas cbp-1(bm2)/+ heterozygyotes had no AFT to ethanol (Fig. [109]5d). This supports a model in which SWI/SNF acts, at least in part, through its regulation of cpb-1 in ethanol response behaviors. Fig. 5. [110]Fig. 5 [111]Open in a new tab cbp-1 is required for AFT. a A model of the cbp-1 gene indicating exon/intron structure and the location and extent of the deletion alleles used in this study. b Normalized cbp-1 expression in our RNA sequencing samples. c-d Animals were tested on 0 mM and 400 mM exogenous ethanol. Relative speeds were calculated as treated over untreated speeds (left Y-axes). AFT is quantified as the recovery of speed between 10 and 30 min of exposure (right Y-axes). All genotypes on the same graph were tested simultaneously on the same plates. Wild-type (N2) animals develop AFT; swsn-9(ok1354) mutant animals do not develop AFT. c Heterozygous cbp-1(ok1491) mutant animals over either a balancer (bal) or wild type (+) chromosome failed to develop AFT. d Heterozygous cbp-1(bm2) mutant animals maintained over a balancer (bal) chromosome failed to develop AFT, while heterozygous cbp-1(bm2) over a wild type chromosome did develop AFT; the recovery of speed was statistically different from wild-type for cbp-1(bm2)/balancer. Error bars indicate S.E.M. Paired one-tailed Student’s t tests were used for statistical comparisons of speeds at 10 and 30 min; One-way ANOVA, with Tukey’s multiple comparison test, was used to compare the initial sensitivity and recovery of speed across the strains. For indicated comparisons: ns, not significant; *p ≤ 0.05; **p ≤ 0.01; ***p ≤ 0.001 Discussion The SWI/SNF chromatin remodeling complex is required in C. elegans for the development of acute functional tolerance to ethanol [[112]10]. In this study, we identified the set of transcripts that are regulated by the SWI/SNF chromatin remodeling complex in adult neurons. To identify the transcripts under the regulatory control of swsn-9 in neurons, we generated a strain in which we provided swsn-9 expression only in neurons in an otherwise swsn-9(ok1354) mutant background (Fig. [113]1). By comparing the transcriptome of these animals to that of swsn-9(ok1354) mutant animals, we could find transcripts regulated by the SWI/SNF complex in neurons. We found that 2412 genes were differentially expressed in this sample, suggesting that these genes are under the control of a SWSN-9 containing SWI/SNF complex in neurons. To identify transcripts that are under the regulatory control of swsn-1 in adults, we used the special temperature sensitive allele swsn-1(os22ts) to provide SWSN-1 function throughout development, then eliminate SWSN-1 function in adults (Fig. [114]2). We compared the transcriptomes of animals with and without SWSN-1 function in adulthood to identifiy candidate AFT mediators. We found that 6813 genes are under regulatory control by swsn-1 in adults. These effects on gene expression may be direct or indirect; for example, some may be the result of SWI/SNF’s regulation of one or more transcriptional repressors or activators that themselves regulate expression of the genes. We observed that more than two times as many genes were regulated by swsn-1 than by swsn-9 in our study; this was consistent with our expectations because SWSN-9 is an accessory subunit that is predicted to be found only in PBAF complexes, while SWSN-1 is a core subunit that is predicted to be present in both BAF and PBAF complexes. Further, this difference in the number of regulated genes may also be due in part to the fact that our swsn-9(ok1354) dataset is confined to neuronally expressed genes, whereas the swsn-1(os22ts) dataset includes all tissues. Interestingly, in addition to different numbers of genes, we also found that the swsn-1 and swsn-9 datasets consisted of mostly non-overlapping sets of genes. If these two SWI/SNF complex subunits were acting together for most or all gene regulation by PBAF, then we would expect that the swsn-9 neuronally-regulated dataset would largely be a subset of the swsn-1-pan tissue regulated dataset. In fact, we observed the contrary; most of the genes in each dataset do not appear in the other dataset. There are several possible explanations for this lack of regulatory overlap. One possibility is that, in addition to acting together in some cases, these two subunits also participate in separate complexes that target different genes. This would imply that SWSN-9 can act with other as yet undetermined proteins to regulate transcription independent of its known function with SWSN-1. This would be a non-canonical function of SWSN-9, and while our data do not exclude this possibility, we do not favor this hypothesis. A second possibility is that swsn-1(os22ts) may be non-null at 25 °C, so that it does provide some function in our restrictive conditions. In such a case, we would see DEGs disrupted by the swsn-9(ok1354) null mutation that are not disrupted by swsn-1(os22ts) because of its residual activity, but those genes nevertheless do require the function of SWSN-1. We favor a third explanation: in our experimental paradigm, swsn-1(os22ts) function is only disrupted in the adult, and its function is more normal during development. In contrast, swsn-9(ok1354) is non-functional throughout development. We hypothesize that the swsn-9-specific DEGs are regulated by SWI/SNF before adulthood. We do not observe their dysregulation by the swsn-1(os22ts) allele because in this work it functions normally during development. Thus, the swsn-9-specific DEGs represent developmental gene regulation defects that persist into adulthood. Our datasets provide a rich resource with which to examine the functions of these SWI/SNF target genes. Although it remains possible that swsn-1 and swsn-9 regulate AFT through independent mechanisms, because we found that loss of either swsn-9 or swsn-1 function causes similar defects in AFT [[115]10], we favor the idea that they act together in the SWI/SNF complex for AFT. Therefore, our primary interest was in genes that are regulated by both of these subunits. We found that a total of 603 genes were regulated in the same direction by both swsn-9 and swsn-1. These are candidate AFT mediators. To begin to characterize these genes, we first looked at overrepresented classes of genes, with the goal of identifying biological processes that are likely to be involved in AFT. SWI/SNF regulates the membrane microenvironment We were particularly intrigued by our observation that the most overrepresented category in the intersection of the swsn-9 and swsn-1 regulated genes are those annotated to be involved in membrane rafts, which are a form of membrane microarchitecture characterized by particular lipid comoposition, thick bilayer and the presence of cholesterol. Membrane microachitecture can provide a platform for regulating the activity of many membrane proteins, both directly and indirectly. For example, membrane characteristics such as lipid makeup, bilayer thickness, and cholesterol content have all been shown to modulate the activity of the BK channel (reviewed in [[116]33, [117]34]), which is a well characterized ethanol target (reviewed in [[118]35, [119]36]). In addition, membrane microdomains such as lipid rafts can regulate the access of proteins to effector molecules, thereby affecting their activity [[120]37, [121]38]. Importantly, we have previously demonstrated that cholesterol, a key component of lipid rafts, is required for the development of AFT [[122]39]. In addition, we have shown that levels of the long chain ω-3 polyunsaturated fatty acid, eicosapentaneoic acid (EPA), tune the levels of AFT; EPA is a component of lipid rafts [[123]40, [124]41]. Taken together, these data point to a probable role for rafts in AFT. Our gene expression data implicate SWI/SNF complexes in the regulation of membrane raft components. These data have led us to propose a molecular model for the development of AFT in which the ethanol sensitivity of ethanol target proteins, such as BK channels, can be modulated by changing the membrane milleu in which they reside. SWI/SNF regulation of factors that modify lipid rafts may be required if one of the mechanisms of AFT is for ethanol target proteins to be moved into or out of lipid rafts to normalize their function if they are exposed to ethanol. AFT in our system is observable within 30 min, and such modulation could be very quick, since it would not rely on transcription or changes in protein abundance. If there are appropriate lipid raft components available, then the system could be poised to respond to ethanol quickly when it is encountered. SWI/SNF acts through cbp-1 to regulate AFT Generally, gene families have expanded in more complex organisms, so our 603 candidate C. elegans genes are homologous to 874 human genes. As a first analysis of our dataset, we asked if any of the homologs of the genes in our candidate dataset have been implicated in alcohol response behaviors in other studies, by searching PubMed for coincidence of the human gene name and alcohol response search terms such as ‘alcohol dependence’ and ‘alcohol use disorder,’ and then further narrowed our search by filtering these by a mention of epigenetic regulation. This narrowed our list to 195 genes (Additional file [125]4). We noted that among the top 20 of these genes (Table [126]1), two, CREBBP (encoding CBP) and EP300 (encoding p300), are homologs of the same C. elegans gene, cbp-1. Both CBP and p300 have been shown to be involved in ethanol effects in several systems [[127]28, [128]30, [129]42–[130]44]. cbp-1 is also interesting to us because CBP is known to act in concert with SWI/SNF complexes in epigenetic modulation of gene transcription [[131]45, [132]46]. Together, these make cbp-1 an excellent candidate mediator of the SWI/SNF effect on AFT. We tested strains heterozygous for two different loss of function mutations in cbp-1 and found that cbp-1/+ animals had defects in the development of AFT, demonstrating that this SWI/SNF target gene is important in this process. We favor a model in which SWI/SNF can modulate the transcription of cbp-1, and, CBP-1 and the SWI/SNF complex act in concert to regulate the transcription of effector genes that are required for the development of AFT. CBP-1 may also interact with other transcription factors to regulate the expression of effector genes. Intriguingly, both p300 and CBP have direct interactions with human CTBP family members in complex regulatory relationships [[133]47–[134]50]. CTBP1 and CTBP2 are transcriptional repressors that are homologous to C. elegans ctbp-1. We previously found that ctbp-1 is required for normal AFT in worms [[135]39], further suggesting a mechanism by which cbp-1 may play a role in AFT. Our ultimate goal is to identify the direct effectors of AFT, to further our molecular understanding of this important aspect of the physiological effects of alcohol. Conclusions This work describes the transcriptional targets of two subunits of the SWI/SNF chromatin remodeling complex in C. elegans. We generated one dataset of genes whose transcription is regulated by SWSN-9 in neurons. The second dataset is genes whose transcription requires SWSN-1 in adult tissues. Both of these genes are required for the development of acute functional tolerance, so we examined the intersection of these datasets, and identified 603 genes whose regulation is dependent on SWSN-1 in adults and SWSN-9 in neurons. Among these genes we expect to find the effectors of AFT. Pathway enrichment analysis identified genes involved in membrane rafts, suggesting that the regulation of membrane microarchitecture may play an important role in AFT. The candidate gene cbp-1 is homologous to mammalian genes that have been implicated in ethanol respose behaviors, and we found that cbp-1 function is required for normal AFT. This dataset provides a rich resource with which to examine the mechanisms underlying the development of acute functional tolerance to ethanol. Methods Strains C. elegans strains were cultured as described previously [[136]51, [137]52]. All strains were grown at 20 °C unless otherwise specified and were derived from the Bristol strain N2. Strains were obtained from the Caenorhabditis Genetics Center or were generated as described below. The lethal cbp-1 alleles were maintained over hT2[bli-4(e937) let-?(q782) qIs48], which serves as a balancer chromosome for LGI and LGIII. The following alleles were used in this study and are described in C. elegans II [[138]53], cited references, or this work: