- Research article
- Open Access
Chromatin accessibility profiling in Neurospora crassa reveals molecular features associated with accessible and inaccessible chromatin
BMC Genomics volume 22, Article number: 459 (2021)
Regulation of chromatin accessibility and transcription are tightly coordinated processes. Studies in yeast and higher eukaryotes have described accessible chromatin regions, but little work has been done in filamentous fungi.
Here we present a genome-scale characterization of accessible chromatin regions in Neurospora crassa, which revealed characteristic molecular features of accessible and inaccessible chromatin. We present experimental evidence of inaccessibility within heterochromatin regions in Neurospora, and we examine features of both accessible and inaccessible chromatin, including the presence of histone modifications, types of transcription, transcription factor binding, and relative nucleosome turnover rates. Chromatin accessibility is not strictly correlated with expression level. Accessible chromatin regions in the model filamentous fungus Neurospora are characterized the presence of H3K27 acetylation and commonly associated with pervasive non-coding transcription. Conversely, methylation of H3 lysine-36 catalyzed by ASH1 is correlated with inaccessible chromatin within promoter regions. Conclusions: In N. crassa, H3K27 acetylation is the most predictive histone modification for open chromatin. Conversely, our data show that H3K36 methylation is a key marker of inaccessible chromatin in gene-rich regions of the genome. Our data are consistent with an expanded role for H3K36 methylation in intergenic regions of filamentous fungi compared to the model yeasts, S. cerevisiae and S. pombe, which lack homologs of the ASH1 methyltransferase.
Functionally distinct chromatin environments are defined by characteristic molecular features including histone modifications, occupancy of DNA- and chromatin-binding proteins, and varying degrees of chromatin accessibility (for a review, see ). Euchromatin is traditionally defined as gene-rich, transcriptionally active and enriched for specific chromatin modifications including histone H3 acetylation (H3Kac), histone H3 lysine 4 methylation (H3K4me1/2/3), and histone H3 lysine 36 trimethylation (H3K36me3) [2,3,4,5,6]. Conversely, heterochromatin is transcriptionally repressed and marked by histone H3 lysine 9 methylation (H3K9me3) [6,7,8,9,10] or histone H3 lysine 27 methylation (H3K27me2/3) [6, 11,12,13,14,15,16]. Chromatin environments can also be defined by the presence of histone variants such as histone H2A.Z, which is found surrounding transcriptional start sites and within vertebrate and plant enhancers [17,18,19,20,21,22,23].
Euchromatic and heterochromatic environments display marked differences in the level of chromatin accessibility to DNA-binding proteins or other regulatory factors. The Assay for Transposase-Accessible Chromatin sequencing (ATAC-seq) method has proven to be useful tool for characterizing chromatin accessibility in eukaryotes . ATAC-seq has been applied in yeasts, animals, and plants to identify promoter and enhancer environments [24,25,26,27]; however, work in non-yeast fungi is limited. A variety of factors appear to contribute to accessible chromatin. For example, certain DNA sequences are thought to exclude nucleosomes and thus contribute to nucleosome depletion near the start of transcription . In addition, binding of pioneer transcription factors within promoter regions can drive the formation of an accessible chromatin environment and promote transcriptional activity [29,30,31,32]. In Neurospora, previous MNase experiments have shown that binding of the light-response master regulator white collar complex (WCC) subunit white collar 2 (wc-2) can change nucleosome organization at promoters of WCC-inducible genes . Additionally, application of ATAC-seq in the pathogenic basidiomycete Cryptococcus neoformans has demonstrated that transcription factor ZNF2 and the Switch/Sucrose non-fermentable (SWI/SNF) complex create chromatin accessibility at filamentation genes during filamentation . Another factor that may contribute to accessible chromatin is pervasive non-coding transcription [34,35,36]. Indeed, a recent study investigated histone modifications and small RNAs at cis-regulatory elements in cnidarians, fungi, and human cells, providing evidence that small RNAs or RNA polymerase activity are evolutionary conserved features of accessible chromatin . While this study characterized small RNA production in the filamentous fungus Neurospora, raising the possibility that transcription may contribute to accessible chromatin, the authors did not identify accessible chromatin regions in this fungus.
Other chromatin features are correlated with inaccessible chromatin. Constitutive heterochromatin domains marked by H3K9me3 and H3K27me3 are reported to be largely inaccessible in higher eukaryotes , but chromatin inaccessibility within these regions has not yet been directly measured in filamentous fungi. In Saccharomyces cerevisiae, methylation of H3K36 is mediated by a single histone methyltransferase, SET2, which associates with elongating RNA polymerase to methylate H3K36 in coding sequences [39,40,41,42]. Although H3K36me3 occurs in active genes, this modification promotes inaccessible chromatin via recruitment of chromatin remodelers and histone deacetylases that act to repress spurious intergenic transcription from cryptic promoters [43,44,45,46,47]. In certain filamentous fungi and in higher eukaryotes, H3K36 methylation is catalyzed by multiple lysine methyltransferase enzymes. Neurospora and the closely related Fusarium fujikuroi encode two H3K36 methyltransferases that have been genetically characterized (SET-2 and ASH1) [48,49,50]. Like yeast SET2, Neurospora SET-2 trimethylates H3K36 along gene bodies and is presumably targeted by elongating RNA polymerase . Neurospora ASH1, instead, deposits mono- and di-methylation on H3K36 residues across promoters and gene bodies of silent genes in a transcriptionally independent manner . Some regions of ASH1 catalyzed H3K36me2 co-localize with H3K27me2/3 and contribute to Polycomb Repressive Complex 2 (PRC2)-mediated repression .
Here, we report the development of an ATAC-seq protocol for filamentous fungi. Using ATAC-seq, we carried out a comprehensive analysis of chromatin features associated with accessible chromatin. We report a diversity of promotor structures in Neurospora and we show that histone acetylation and small RNA production are highly correlated with accessible chromatin, whereas repressive modifications including H3K9me3, H3K27me2/3 and ASH1-catalyzed H3K36 methylation are correlated with inaccessible chromatin in this fungus.
Characterization of accessible chromatin regions in Neurospora
To characterize accessible regions in Neurospora, we performed ATAC-seq on wild type mycelial cultures. Because this assay is sensitive to contamination by mitochondria, we first performed ATAC-seq in nuclei that were sorted by flow cytometry and nuclei that were not sorted. In both cases, sequence reads mapped to the Neurospora genome assembly at > 98.4% (Table S1). Fewer reads mapped to the mitochondrial genome in sorted nuclei (~ 15%) when compared to unsorted nuclei (~ 27%) (Fig. S1A; Table S1). We next examined the distribution of ATAC-seq reads from sorted versus unsorted nuclei. We first identified accessible chromatin regions (ACRs) using Genrich from wild type sorted and unsorted nuclei (using standard parameters: minimum area under the curve of 200 bp, significance threshold > 1, and p-value < 0.01). We identified 7132 ACRs unsorted cells and 5181 ACRs sorted cells from the wild-type strain. Of the ACRs we identified in sorted cells, 92% overlapped with ACRs called in unsorted cells by bedtools intersect (n = 4772; Table S2). We then plotted ATAC-seq enrichment in sorted and unsorted cells centered on the ACRs from unsorted nuclei. Comparison of these two experiments revealed that global patterns of chromatin accessibility in Neurospora is unchanged in sorted and unsorted nuclei (Fig. S1B). This is confirmed with an example genome browser image comparing ATAC-seq reads of sorted and unsorted nuclei (Fig. S1C). Based on these data, we concluded that sorting nuclei removed some mitochondrial contamination, but also reduced the sensitivity of the ATAC-seq assay. Thus, it is more efficient and cost-effective to carry out ATAC-seq analysis using unsorted nuclei followed by downstream removal of mitochondrial reads in silico. All remaining experiments were performed with unsorted nuclei.
We next wanted to determine what chromatin environments contain high enrichment of accessible chromatin. To do this, we first visually inspected ATAC-seq enrichment across the Neurospora genome using the Integrated Genomics Viewer (IGV) . We found that chromatin accessibility is highest in intergenic regions and limited or absent in gene bodies (Fig. 1A). Additionally, facultative and constitutive heterochromatin marked by H3K27me2/3 or H3K9me3, respectively, have low levels of chromatin accessibility (Fig. 1A).
To quantify the level of accessibility in heterochromatin and euchromatin, we used bedtools intersect to determine the number of ACRs that overlap with H3K27me2/3 or H3K9me3 (Table S3). H3K27me2/3 is enriched in 309 domains covering 6% of the genome . We identified 68/7132 ACRs that overlap with called H3K27me2/3 enriched regions. The average density of ACRs in these domains is 0.27 ACRs per 10 kb. We also wanted to quantify the number and density of ACRs in constitutive heterochromatin marked by H3K9me3. In total, 454 domains are enriched for H3K9me3 covering approximately 18% of the genome . 198 ACRs overlapped with H3K9me3 with an average peak density of 0.32 ACRs per 10 kb. The remaining 6866 ACRs are found in regions that lack H3K9me3 or H3K27me2/3, with an average density in euchromatin of 2.25 ACRs per 10 kb.
In general, accessible chromatin regions appeared to be restricted to promoter regions of Neurospora genes (Fig. 1A). To determine if accessibility was a global feature of promoters, we averaged ATAC-seq data across all promoters and plotted the distribution (Fig. 1B). This confirmed that ACRs are present in the 5′ regions of genes. We next wanted to ask if Neurospora promoters exhibit distinct patterns of chromatin accessibility. We did this by performing k-means clustering and plotting a heatmap of ATAC-seq read enrichment across 3000 bp (+/− 1500 bp) of DNA centered on the transcription start site (TSS) of each Neurospora gene (Fig. 1C). It was difficult to resolve organization of accessible regions at the single-gene level because the Neurospora genome is highly compact; thus, many intergenic regions contain a second promoter within the 3000 base pair window plotted for each gene (e.g. divergently transcribed genes share the same intergenic space; Fig. 1C, Left). To better ascertain the diversity of accessible chromatin structures in Neurospora gene regulatory regions, we examined the distribution of ATAC-seq reads only for promoters that could be unambiguously assigned to a single gene (i.e. gene promoters without another annotated transcriptional start site within 1500 base pairs in either direction, which is the average intergenic space in Neurospora (Table S4; see Materials and Methods) . We performed k-means clustering and plotted ATAC-seq read enrichment across the the TSSs of this gene set to group genes with similar patterns of chromatin accessibility within promoters and 5′ intergenic regions (Fig. 1C, Right; k = 6). This analysis reveals that the size and pattern of accessible chromatin varies for genes throughout the Neurospora genome. For example, Cluster 1 and 2 show broad, highly accessible domains, with Cluster 1 showing some extension into the gene body compared to Cluster 2. Cluster 3 shows a narrow region of chromatin accessibility at the TSS, whereas Cluster 4 has a similar enrichment pattern with a moderately lower level of enrichment. Cluster 5 enrichment is lower and broader than Clusters 1–4, and Cluster 6 shows general depletion of ATAC-seq signal.
Given the differences in ATAC-seq enrichment levels and patterns in Neurospora promoters, we wanted to know if there was a correlation between promoter structure and transcript level. We plotted the log2-normalized transcripts per million (TPM) of genes within each k-means cluster (Fig. 1D). In general, the average level of expression was higher in clusters with more ATAC-seq accessibility, whereas genes in Cluster 6 had lower transcript levels on average. Although we did observe differences in the average transcript levels between clusters, all clusters had genes with either very high or very low expression indicating that chromatin accessibility in Neurospora does not necessarily result in higher mRNA transcription (Fig. 1D). Instead, it is likely that accessible chromatin regions correspond to binding sites of both transcriptional activators and repressors.
We next asked if distinct promoter structures were enriched for any functional categories of genes. In particular, we were interested in genes with large, hyperaccessible chromatin regulatory regions. We used ACRs identified from unsorted nuclei, as described above (n = 7132), and we ranked them by size. We annotated ACRs greater than 2 kb (n = 84) to their closest genes (Table S2; Fig. 1E) and performed Gene Ontology (GO) analysis. This revealed that clusters of genes with large, hyperaccessible chromatin regions are significantly enriched for GO categories “sequence-specific DNA binding”, “transcription factor activity”, and “regulation of transcription” (Fig. 1F; Table S5). Closer inspection of these genes also reveals that they are involved in cellular processes including light response, conidial separation, and other master regulatory functions, including the carbon catabolite repressor cre-1, and a member of the light-response pioneer factor complex WCC, white collar 1 (wc-1).
Because gene annotation software determines annotated genes based on distance from the feature of interest, and these regions could annotate to more than one gene, we were concerned that the annotation software was giving us inaccurate results. We therefore completed the same analysis using Cluster 1 and Cluster 2 from Fig. 1C (Table S4). Our analyses revealed that Cluster 2 is enriched for genes which are annotated to have transcriptional roles (i.e. RNA Polymerase 2 (Pol2) activity, TF binding, and sequence-specific DNA binding), while Cluster 1 is more enriched for signal transduction activity (Fig. S2; Table S5). This supports our finding that genes with large, hyperaccessible chromatin regions are enriched for genes with regulatory functions.
ACRs are associated with high nucleosome turnover and transcription factor binding
Pioneer factors are DNA-binding proteins that can displace nucleosomes to create accessible chromatin regions. To determine if ACRs are associated with higher rates of nucleosome turnover, we examined rates of nucleosome exchange across the Neurospora genome using an inducible hH3-3xFLAG protein. We first introduced an H3-3xFLAG fusion protein (hH3:3xflag) under the control of a copper-inducible promoter, and we confirmed that the H3-3xFLAG protein is induced under copper limiting conditions and incorporated into chromatin (Fig. 2A). We then performed ChIP-seq following MNase digestion (MNase-ChIP) to identify regions with high H3-3xFLAG enrichment, reasoning that incorporation of newly synthesized H3-3xFLAG would preferentially occur in regions with high rates of nucleosome turnover (i.e. due to competition between TFs and nucleosomes). Conceptually similar approaches have been applied in other fungal systems and recently in N. crassa [54,55,56]. We plotted H3-3xFLAG enrichment after 8 h of induction and found high levels of H3 turnover correspond to ACRs (Figs. 2B, S3). We then plotted nucleosome turnover (i.e. H3-3xFLAG enrichment) across all ACRs, which revealed high nucleosome turnover at these regions (Fig. 2C). We confirmed significant overlap of ATAC-seq peaks and H3-3xFLAG peaks using a permutation test (p < 0.001).
We next asked if ATAC-seq could be used to identify pioneer transcription factors in Neurospora by asking if disruption of the White Collar Complex (WCC) would lead to altered patterns of accessible chromatin measured by ATAC-seq. Transcription factors have been demonstrated to provide a mechanism for creating chromatin accessibility at gene promoters, and Sancar and colleagues previously implicated WCC in maintenance of accessible chromatin using MNase assays in wild type and strains lacking a functional WCC. Thus, WCC is proposed to function as a pioneer transcription factor . We asked if ATAC-seq data might provide comparable or even more robust results than MNase assays by performing ATAC-seq assays in wild type and the wc-2 mutant. We plotted ATAC-seq reads across all WCC binding sites  for wild type and the wc-2 mutant. We observed a striking reduction of chromatin accessibility in the Δwc-2 mutant at these sites. These data confirm prior work indicating that the WCC functions as a pioneer transcription factor (Fig. 2D,E). To demonstrate that reduced chromatin accessibility at WC-2 binding sites in Δwc-2 is not an experimental artefact, we compared heatmaps of ATAC-seq data from wild type and Δwc-2 and found no global reduction in ATAC-seq read enrichment genome-wide (Fig. S4; Table S6). Taken together, these data suggest that ATAC-seq is a robust method for examining the relationship between transcription factors and chromatin accessibility in filamentous fungi.
Accessible chromatin is associated with non-coding transcriptional activity in N. crassa
We next wanted to identify additional chromatin features that correlate with accessible chromatin. Recent work identified small abortive transcripts (capped small RNAs; csRNAs) that map to Neurospora intergenic regions, raising the possibility that small RNAs or active transcription is a feature of accessible chromatin in filamentous fungi as reported for other eukaryotes [36, 37]. These small RNAs are short (< 60 bp); can be stable or unstable, as determined by read depth when compared to total RNA; and can be uni-directional or bi-directional. To determine if these abortive transcripts are features of accessible regions in Neurospora, we compared the previously published csRNA transcript data to our ATAC-seq data. We first asked if regions that generate csRNA transcripts are accessible. We created a heatmap centered on the TSS of csRNAs sorted by stability and directionality as defined above. When we compare these transcripts to chromatin accessibility data, we find that the TSSs of csRNAs are accessible (Fig. 3A). Similarly, a plot of csRNA transcripts across all ACRs reveals that csRNA transcripts occur adjacent to many of them, however not all ACRs contain high levels of csRNAs (Fig. 3B).
To quantify the extent of overlap between chromatin accessibility and csRNA transcripts, we used bedtools intersect to identify ACRs that contain csRNA transcripts (Table S7) and found that 73% of ACRs (n = 5213) overlap with at least one annotated csRNA transcript. Conversely, only 60% of csRNA transcripts (n = 9709) are found within accessible chromatin regions. Many of these csRNA transcripts appeared to flank accessible chromatin regions. This is a feature that is common among higher eukaryotic promoters and enhancers [58, 59]. We therefore examined enrichment of csRNAs within 100 bp of accessible chromatin in Neurospora using bedtools closest. We found that 83% of ACRs either contained or occurred within 100 bp of an annotated csRNA transcript (n = 5939). Conversely, approximately 9% of csRNAs were identified within these 100-bp regions (n = 1570). In total, 69% of csRNA transcripts either overlapped or were within 100 bp of an ACR (n = 11,393). Together, these data indicate that non-coding transcriptional activity is a widespread feature of accessible regions in Neurospora.
H3 acetylation is associated with ACRs and H3K36me3 marks inaccessible chromatin
To identify additional chromatin features that are associated with accessible regions, we plotted enrichment of ChIP-seq data for various histone modifications surrounding accessible regions (Fig. 4A). We performed k-means clustering to construct heatmaps using enrichment data for a suite of chromatin features, cs-RNA reads, and mRNA-seq reads to identify chromatin features that are associated with accessible and inaccessible chromatin regions. A heatmap of read density centered on ACRs reveals that virtually all regions of chromatin accessibility are adjacent to nucleosomes containing H3K27ac. The direction of productive transcription is evident by visualizing enrichment of mRNA-seq reads near the ACRs. The heatmap in Fig. 4A further supports the conclusion that accessible chromatin regions are not necessarily associated with active transcription, as shown by the low level of mRNA enrichment in Clusters 3 and 4. H3K4me was correlated with productive transcripts but was not associated with all ACRs.
H2A.Z is reported to promote chromatin accessibility in in embryonic stem cells . We compared previously published H2A.Z ChIP-seq data  to our ATAC-seq data and found that H2A.Z was enriched flanking most, but not all ACRs (e.g. see Cluster 3). We then carried out ATAC-seq experiments in wild type and ΔhH2Az, which revealed that loss of H2A.Z did not alter global patterns of chromatin accessibility in Neurospora (Fig. 4B). Thus, H2A.Z is not a major determinant of accessible chromatin in Neurospora.
Recent work has identified a novel role of ASH1-mediated H3K36me2 in gene repression in Neurospora . Namely, Bicocca et al. identified regions of the Neurospora genome that were both transcriptionally repressed and marked by ASH1-catalyzed H3K36me2 or both H3K27me2/3 and ASH1-catalyzed H3K36me2. We compared previously published enrichment data for ASH1-dependent H3K36me2 with our chromatin accessibility data and found a small number of ACRs that are marked by these repressive histone modifications (Cluster 4 of Fig. 4A). Small regions of chromatin accessibility have been identified in facultative heterochromatin in other organisms. In particular, Polycomb Response Elements (PREs) that function to recruit Polycomb complexes are nucleosome free and accessible, while surrounding regions are inaccessible [61, 62]. We asked if these ACRs were enriched for any specific DNA sequences, which would be consistent with a PRE-like function. We used Hypergeometric Optimization of Motif EnRichment (HOMER) to search for known and unknown DNA motifs enriched under ACRs within PRC2-target domains. Motif analysis did not identify unknown enriched sequences with high confidence (Fig. S5). Interestingly, the top two known motifs identified by HOMER contain a complete or partial identity to the Neurospora telomere repeat sequence, TTAGGG, which is sufficient to recruit PRC2 to establish a domain of H3K27me2/3 .
Next, we compared the pattern of chromatin accessibility, H3K27me2/3 and ASH1-catalyzed H3K36 di-methylation across the TSSs of genes with unambiguous 5′ regions, as in Fig. 1C. A heatmap centered on the TSS of distinct promoters shows that promoters marked by H3K27me2/3 and H3K36me2 or H3K36me2 alone are highly inaccessible (Fig. 4C). Together, these data demonstrate that H3K36 methylation within euchromatin domains is correlated with inaccessible chromatin in promoters.
We carried out a global analysis of chromatin accessibility in a filamentous fungus using ATAC-seq. Filamentous fungi are important in ecological, industrial, clinical, and evolutionary contexts, and Neurospora crassa is an important model for chromatin and epigenetic studies due to shared pathways with higher eukaryotes [9, 16]. Our analyses provide experimental evidence that facultative and constitutive heterochromatin regions marked by H3K27me2/3 and H3K9me3 are largely inaccessible in Neurospora. This result is interesting given that fungi lack Polycomb Repressive Complex 1 (PRC1), which is proposed to play a critical role in compacting chromatin within H3K27me3-enriched regions in higher eukaryotes [64, 65]. We also found that promoter chromatin accessibility is not tightly correlated with transcript levels. This suggests that accessible chromatin is occupied by both activators and repressors and therefore should not be used as a measure of active gene expression in filamentous fungi. Conversely, there is a strong correlation between accessible chromatin and H3 turnover. These data support a model in which accessible chromatin results from a competition between promoters and DNA-binding transcription factors. Indeed, we confirm using ATAC-seq that WCC is required to maintain accessible chromatin at WCC binding sites. It is likely that WCC is needed to drive formation of accessible chromatin and allow transcription factors to bind nearby sequences.
Our findings in combination with those of Duttke et al. demonstrate that abortive transcription is a key feature of accessible chromatin in Neurospora. The presence of non-coding transcription from distal accessible regions in Neurospora is reminiscent of enhancers in higher eukaryotes. Similar results were obtained in the cnidarian Nematostella [37, 66] supporting that model that non-coding transcription is an evolutionarily conserved mechanisms of eukaryotic gene regulation. In higher eukaryotes, short unstable RNAs typically flank enhancers and promoters [58, 59]. It is unclear if the small RNAs produced from accessible regions provide a specific function: (1) They could simply be results of promiscuous Pol II activity and themselves be non-functional; (2) These could be regulatory RNAs; or (3) The act of transcription itself may be important in maintaining accessible chromatin in filamentous fungi. Recently, a new model of promoter/enhancer characterization has been proposed in which genomic regulatory elements have enhancer and promoter potential . This model defines promoters and enhancers as regulatory elements at any accessible chromatin region bound by transcription factors, with the ability to recruit Pol II and initiate transcription at either edge, which can positively influence transcription initiation at other regulatory elements. Therefore, it is likely that these RNAs themselves may not be regulatory, but Pol II recruitment to regulatory elements leads to some level of transcription, and the occupancy of Pol II at these sites helps to maintain their accessibility.
We did find a subset of annotated csRNA loci that were not associated with ACRs; however, these may reflect differences in growth conditions used for our ATAC-seq experiments and the csRNA experiments carried out by Duttke et al. Experimental variations including growth conditions, single versus mixed tissue type, and light exposure can produce quite varied results in Neurospora expression profiles. It will be important, therefore, to perform more in-depth analysis at these putative cis-regulatory regions to determine the role and abundance of these types of RNAs. Importantly, the growth conditions of the csRNA-seq experiments performed by Duttke et al. vary from our ATAC-seq protocol: csRNA-seq cultures were grown for 2 days at constant light, while cultures used for ATAC-seq were grown overnight at constant light. Therefore, matched experiments will allow for more fine-tuned associations between csRNA production and chromatin accessibility.
The presence of H3K27me2/3 in a small number of accessible chromatin regions also presents an interesting finding, suggesting that PRE equivalents may exist in filamentous fungi. Recent work has identified telomere repeat sequences and the presence of H3K36me2 as prerequisites for H3K27me via PRC2 , and a recent study identified PRC2 accessory subunit (PAS) as a PRC2-interacting protein that recruits the complex to a subset of target domains . More work is needed to determine if these ACRs within PRC2-target domains are indeed playing a role in recruitment of PRC2 and establishment of H3K27me.
Our results support the findings reported in Bicocca et al. , and suggest that H3K36me may contribute to gene repression by maintaining an inaccessible chromatin structure within silent promoters. Indeed, H3K36 methylation may be a major determinant of inaccessible chromatin within active regions of the genome. Anti-sense transcription at the Neursopora frequency (frq) locus has been demonstrated to regulate circadian clock function via H3K36me3 [69, 70]. The anti-sense transcript of frq (qrf) is proposed to stall RNA Polymerase II at the frq locus, resulting in premature termination of transcription and establishment of H3K36me3 that precludes binding by the WCC . Our global data are consistent with a repressive role for H3K36me3 when found in regulatory regions. The recent work investigating ASH1-dependent H3K36me in Neurospora could provide a model for gene repression in other fungi, especially those that lack PRC2 and H3K27 methylation [71, 72]. F. fujikuroi strains lacking ASH1 mediated H3K36me2 exhibited increased expression of secondary metabolism biosynthetic clusters . Despite a lack of ASH1, H3K36me3 also performs a repressive role in S. cerevisiae and S. pombe, where it is responsible for preventing cryptic transcription and heterochromatic silencing via recruitment of nucleosome remodelers and histone deacetylases [43, 73, 74]. It is possible, therefore, that H3K36me plays a conserved role in gene repression in fungi, namely by contributing to an inaccessible local chromatin architecture within both gene bodies and promoters, as suggested by the work here.
Profiling accessible chromatin by ATAC-seq uncovered key features of accessible and inaccessible chromatin in the N. crassa genome. Accessible chromatin is found in promoters and 5′ regulatory regions of genes. Distinct patterns of accessibility exist, with large, highly accessible intergenic regions frequently associated with master regulatory genes. Accessible regions depend on transcription factor binding and are associated with high rates of nucleosome turnover. H3K27 acetylation is strongly correlated with accessible chromatin, whereas ASH-1-dependent H3K36 methylation, H3K27 methylation, and H3K9 methylation are correlated with inaccessible chromatin.
Culture growth conditions and strain construction
All liquid cultures were grown in Vogel’s minimal medium with 1.5% sucrose with shaking at 32 °C for 18 h. Strains were constructed using standard Neurospora protocols [75, 76]. All strains used in this study are listed in Table S8 and primers used for strain construction are listed in Table S9. To construct strain S564 containing a hH3-3xflag under control of the copper-regulated tcu-1 promoter , the hH3 coding sequence was amplified using primers, hH3 inF FP and hH3 inF RP, and then cloned using the In-Fusion HD Cloning HD kit (Clontech) into PacI site of a linearized version of pCCG-3xFLAG plasmid . An hH3-3xflag fragment was amplified using primers hH3 CDS FP and hH3 CDS + FLAG RP from the constructed vector (pCCG-C-hH3-3xFLAG). A fragment of the tcu-1 promoter was amplified using primers Ptcu-1 FP and Ptcu-1 RP + h3 CDS. To assemble Ptcu-1-hH3-3xflag, equal molar amounts of each DNA fragment was mixed and used as a template for PCR with primers, Ptcu-1 FP and 3xFLAG RP. The Ptcu-1-hH3-3xflag was cloned into the pMD19 TA cloning vector.
The csr-1 locus containing csr-1 5′ flank-CDS-csr-1 3′ flank was amplified by using primers, csr-1 5′ FLANK FP3 and csr-1 3′ FLANK RP3, and cloned into pMD19 TA cloning vector to create plasmid pMD19-CSR-1. To generate the CDS-less version of a linearized version of this vector, inverse PCR was performed using primers csr-1 5 flank RP + Ptcu-1 and csr-1 3 flank FP + FLAG. The Ptcu-1-hH3-3xflag fragment was cloned into this linearized vector using In-Fusion cloning. The construct was integrated into the csr-1 locus by transformation and selection on Cyclosporin A as described .
ATAC-seq was performed as described with minor modifications [24, 25]. Briefly, overnight mycelial cultures were lysed in lysis buffer (15 mM Tris pH 7.5; 2 mM EDTA; 0.5 mM spermine; 80 mM KCl; 20 mM NaCl; 15 mM (or 0.1% v/v) β-me; 0.3% TrixtonX-100) by placing mycelia in a petri dish and chopping mycelia with a razor blade. Nuclei were isolated by passing the lysate through miracloth followed by centrifugation at 1000 rpm for 5 min in a microcentrifuge. Nuclei were resuspended in 200 uL of lysis buffer and stored on ice for a maximum of 2 h. To perform transposon integration, nuclei were centrifuged at 1000 rpm for 5 min and washed in 50 uL 1X TAPS buffer (10 mM TAPS-NaOH, 5 mM MgCl2, 10% Dimethylformamide, pH 8.0). After centrifugation, Tn5 integration was performed in 30 uL of tagmentation buffer using modified Tn5 transposase pre-loaded with Illumina sequencing adaptors; samples were incubated for 30 min at 37 °C and fragmented DNA was isolated using the Qiagen MinElute kit (Cat#. 28,004). Sequencing libraries were constructed using Phusion polymerase (ThermoFisher Cat # F549N) and amplified with primers provided in . Libraries were cleaned using SeraPure beads (Fisher Cat. # 09–981-123) before submission for sequencing as described in [80, 81]. Libraries were sequenced on an Illumina NextSeq 500 instrument by Georgia Genomics Facility at University of Georgia.
ATAC-seq reads were trimmed and quality scored with TrimGalore! (https://github.com/FelixKrueger/TrimGalore). Trimmed reads were mapped to the NC12 genome (GCA_000182925.2_NC12) with BowTie2 using the –very-sensitive option and setting the maximum insert size set at 2000 bp .
Analysis of transcriptional start sites (TSSs) was performed on genes with unambiguous regulatory regions. Briefly, a list of promoters that could be unambiguously assigned to a single promoter was generated from a list of all TSSs in the Neurospora genome created with the GenomicRanges package in R [83, 84]. This list was compared to the genome annotation file using bedtools closest to identify the distances between TSSs . Only genes with TSSs that were at least 1500 base pairs apart were classified as genes with unambiguous promoters (n = 5200).
ACRs were called using Genrich (https://github.com/jsh58/Genrich) using the ATAC setting (−j). Reads were separated by insert size using R , and A less than 115 bp were called using MACS2 with a q-value cutoff of 0.01 . FunCat2  was used for genome ontology analysis. Heatmaps and metaplots were made using deepTools . All peak comparisons were performed with bedtools . ACRs were annotated using ChIPseeker . We used the peakPermTest function of the ChIPpeakAnno R package to test for significant overlap between ACRs and H3-3xFLAG turnover .
Paired-end RNA-seq data was acquired from Joint Genomes Institute Community Sequencing Project 54,089 and were mapped with HISAT2  using standard options. Counts were created with featureCounts  and analyzed in R . All heatmaps were made using deepTools, and reads were normalized in deepTools for counts per million (CPM) .
Chromatin immunoprecipitation (ChIP-seq)
ChIP-seq was performed as described in [52, 80, 81] with the following antibodies: H3K27me2/3 (Active Motif Cat# 39535, Lot# 1671401)2; H3K4me1 (Abcam Cat# ab9885); H3K4me2 (Active Motif Cat# 29679, Lot# 31713007); H3K4me3 (Abcam Cat# ab1012); H3K36me3 (Abcam Cat# ab9050, Lot# GR300388–1); H3K9me3 (Active Motif Cat# 39161, Lot# 14418003); H3K27ac; FLAG (Sigma Cat# F1804, Lot# SLBS3530V). Libraries were created and analyzed as described in . Libraries were sequenced by Georgia Genomics Facility at University of Georgia on a NextSeq 500 instrument. Motif analysis was performed with HOMER findMotif.pl using standard options . ChIP-seq peaks were called using MACS2, with a q-value cutoff of 0.01 .
MNase-seq for hH3-3xFLAG was performed in cells eight hours after copper induction. Cells were crosslinked in 1% formaldehyde and transferred to NPS buffer (10 mM Tris-HCl, pH 7.5; 50 mM NaCl; 25 mM MgCl; 1 mM CaCl) and sonicated to shear the chromatin. MNase was performed with 20 U/μl Takara MNase (Cat # 2910A) at 37 °C for 60 min. MNase digestion was quenched with 10 mM EDTA and 100 mM NaCl. Anti-flag ChIP was performed as described above. Previously published data for H3K36me2 in the ∆set-2 mutant, , hH2Az-GFP , and WC-2,  were re-mapped to the current Neurospora genome assembly and analyzed as described above.
csRNA-seq was downloaded from the NCBI Short Read archive and mapped to the NC12 genome (GCA_000182925.2_NC12) and analyzed using csRNA tools developed in HOMER as described in ..
Availability of data and materials
The new datasets supporting the conclusions of this article are available in the NCBI Gene Expression Omnibus (GEO) repository, accession number GSE154497. csRNA-seq was downloaded from the NCBI Short Reach Archive, accession #s: SRR9916743, SRR9916744, SRR91745, SRR9916746, SRR916747, SRR916748, SRR916749, SRR916750. Previously published ChIP-seq data for H3K36me2, hH2A.Z-GFP, and WC-2 were downloaded from the NCBI Short Reach Archive, accession #s: SRS3667122, SRS3667122, SRR1578072. RNA-seq data from wild type were obtained from the Joint Genome Institute Genome Portal (https://genome.jgi.doe.gov/portal/DetofExpression/DetofExpression.info.html).
Histone H3 acetylation
Histone H3 lysine 4 methylation
Histone H3 lysine 36 trimethylation
H3 lysine 9 methylation
Histone H3 lysine 27 methylation
Assay for Transposase-Accessible Chromatin sequencing
White collar complex
- wc-2 :
White collar 2
- wc-1 :
White collar 1
Polycomb Repressive Complex 2
Transcription start site
Accessible chromatin regions
Transcripts per million
RNA Polymerase 2
ChIP-seq following MNase digestion
Capped small RNAs
Polycomb Response Elements
Hypergeometric Optimization of Motif EnRichment
Polycomb Repressive Complex 1
PRC2 accessory subunit
- frq :
- qrf :
Anti-sense transcript of frq
Klemm SL, Shipony Z, Greenleaf WJ. Chromatin accessibility and the regulatory epigenome. Nat Rev Genet. 2019;20(4):207–20. https://doi.org/10.1038/s41576-018-0089-8.
Liu CL, Kaplan T, Kim M, Buratowski S, Schreiber SL, Friedman N, et al. Single-nucleosome mapping of histone modifications in S. cerevisiae. PLoS Biol. 2005;3(10):e328.
Pokholok DK, Harbison CT, Levine S, Cole M, Hannett NM, Lee TI, et al. Genome-wide map of nucleosome acetylation and methylation in yeast. Cell. 2005;122(4):517–27. https://doi.org/10.1016/j.cell.2005.06.026.
Mikkelsen TS, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, et al. Genome-wide maps of chromatin state in pluripotent and lineage-committed cells. Nature. 2007;448(7153):553–60. https://doi.org/10.1038/nature06008.
Rando OJ, Winston F. Chromatin and transcription in yeast. Genetics. 2012;190(2):351–87. https://doi.org/10.1534/genetics.111.132266.
Ho JW, Jung YL, Liu T, Alver BH, Lee S, Ikegami K, et al. Comparative analysis of metazoan chromatin organization. Nature. 2014;512(7515):449–52. https://doi.org/10.1038/nature13415.
Lippman Z, Gendrel AV, Black M, Vaughn MW, Dedhia N, McCombie WR, et al. Role of transposable elements in heterochromatin and epigenetic control. Nature. 2004;430(6998):471–6. https://doi.org/10.1038/nature02651.
Cam HP, Sugiyama T, Chen ES, Chen X, FitzGerald PC, Grewal SI. Comprehensive analysis of heterochromatin- and RNAi-mediated epigenetic control of the fission yeast genome. Nat Genet. 2005;37(8):809–19. https://doi.org/10.1038/ng1602.
Lewis ZA, Honda S, Khlafallah TK, Jeffress JK, Freitag M, Mohn F, et al. Relics of repeat-induced point mutation direct heterochromatin formation in Neurospora crassa. Genome Res. 2009;19(3):427–37. https://doi.org/10.1101/gr.086231.108.
Sasaki T, Lynch KL, Mueller CV, Friedman S, Freitag M, Lewis ZA. Heterochromatin controls gammaH2A localization in Neurospora crassa. Eukaryot Cell. 2014;13(8):990–1000. https://doi.org/10.1128/EC.00117-14.
Zhang X, Clarenz O, Cokus S, Bernatavichute YV, Pellegrini M, Goodrich J, et al. Whole-genome analysis of histone H3 lysine 27 trimethylation in Arabidopsis. PLoS Biol. 2007;5(5):e129. https://doi.org/10.1371/journal.pbio.0050129.
Jamieson K, Rountree MR, Lewis ZA, Stajich JE, Selker EU. Regional control of histone H3 lysine 27 methylation in Neurospora. Proc Natl Acad Sci U S A. 2013;110(15):6027–32. https://doi.org/10.1073/pnas.1303750110.
Connolly LR, Smith KM, Freitag M. The Fusarium graminearum histone H3 K27 methyltransferase KMT6 regulates development and expression of secondary metabolite gene clusters. PLoS Genet. 2013;9(10):e1003916. https://doi.org/10.1371/journal.pgen.1003916.
Basenko EY, Sasaki T, Ji L, Prybol CJ, Burckhardt RM, Schmitz RJ, et al. Genome-wide redistribution of H3K27me3 is linked to genotoxic stress and defective growth. Proc Natl Acad Sci U S A. 2015;112(46):E6339–48. https://doi.org/10.1073/pnas.1511377112.
Margueron R, Reinberg D. The Polycomb complex PRC2 and its mark in life. Nature. 2011;469(7330):343–9. https://doi.org/10.1038/nature09784.
Lewis ZA. Polycomb Group Systems in Fungi: new models for understanding Polycomb repressive complex 2. Trends Genet. 2017;33(3):220–31. https://doi.org/10.1016/j.tig.2017.01.006.
Guillemette B, Bataille AR, Gevry N, Adam M, Blanchette M, Robert F, et al. Variant histone H2A.Z is globally localized to the promoters of inactive yeast genes and regulates nucleosome positioning. PLoS Biol. 2005;3(12):e384.
Barski A, Cuddapah S, Cui K, Roh TY, Schones DE, Wang Z, et al. High-resolution profiling of histone methylations in the human genome. Cell. 2007;129(4):823–37. https://doi.org/10.1016/j.cell.2007.05.009.
Bargaje R, Alam MP, Patowary A, Sarkar M, Ali T, Gupta S, et al. Proximity of H2A.Z containing nucleosome to the transcription start site influences gene expression levels in the mammalian liver and brain. Nucleic Acids Res. 2012;40(18):8965–78. https://doi.org/10.1093/nar/gks665.
Weber CM, Ramachandran S, Henikoff S. Nucleosomes are context-specific, H2A.Z-modulated barriers to RNA polymerase. Mol Cell. 2014;53(5):819–30. https://doi.org/10.1016/j.molcel.2014.02.014.
Dai X, Bai Y, Zhao L, Dou X, Liu Y, Wang L, et al. H2A.Z represses gene expression by modulating promoter nucleosome structure and enhancer histone modifications in Arabidopsis. Mol Plant. 2017;10(10):1274–92. https://doi.org/10.1016/j.molp.2017.09.007.
Courtney AJ, Kamei M, Ferraro AR, Gai K, He Q, Honda S, et al. Normal patterns of histone H3K27 methylation require the histone variant H2A.Z in Neurospora crassa. Genetics. 2020;216(1):51–66. https://doi.org/10.1534/genetics.120.303442.
Lu Z, Marand AP, Ricci WA, Ethridge CL, Zhang X, Schmitz RJ. The prevalence, evolution and chromatin signatures of plant regulatory elements. Nat Plants. 2019;5(12):1250–9. https://doi.org/10.1038/s41477-019-0548-z.
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(12):1213–8. https://doi.org/10.1038/nmeth.2688.
Lu Z, Hofmeister BT, Vollmers C, DuBois RM, Schmitz RJ. Combining ATAC-seq with nuclei sorting for discovery of cis-regulatory regions in plant genomes. Nucleic Acids Res. 2017;45(6):e41. https://doi.org/10.1093/nar/gkw1179.
Satpathy AT, Saligrama N, Buenrostro JD, Wei Y, Wu B, Rubin AJ, et al. Transcript-indexed ATAC-seq for precision immune profiling. Nat Med. 2018;24(5):580–90. https://doi.org/10.1038/s41591-018-0008-8.
Lin J, Zhao Y, Ferraro AR, Yang E, Lewis ZA, Lin X. Transcription factor Znf2 coordinates with the chromatin remodeling SWI/SNF complex to regulate cryptococcal cellular differentiation. Commun Biol. 2019;2(1):412. https://doi.org/10.1038/s42003-019-0665-2.
Segal E, Fondufe-Mittendorf Y, Chen L, Thastrom A, Field Y, Moore IK, et al. A genomic code for nucleosome positioning. Nature. 2006;442(7104):772–8. https://doi.org/10.1038/nature04979.
Zaret KS, Carroll JS. Pioneer transcription factors: establishing competence for gene expression. Genes Dev. 2011;25(21):2227–41. https://doi.org/10.1101/gad.176826.111.
Magnani L, Eeckhoute J, Lupien M. Pioneer factors: directing transcriptional regulators within the chromatin environment. Trends Genet. 2011;27(11):465–74. https://doi.org/10.1016/j.tig.2011.07.002.
Soufi A, Garcia MF, Jaroszewicz A, Osman N, Pellegrini M, Zaret KS. Pioneer transcription factors target partial DNA motifs on nucleosomes to initiate reprogramming. Cell. 2015;161(3):555–68. https://doi.org/10.1016/j.cell.2015.03.017.
Meers MP, Janssens DH, Henikoff S. Pioneer factor-nucleosome binding events during differentiation are motif encoded. Mol Cell. 2019;75(3):562–75 e565. https://doi.org/10.1016/j.molcel.2019.05.025.
Sancar C, Ha N, Yilmaz R, Tesorero R, Fisher T, Brunner M, et al. Combinatorial control of light induced chromatin remodeling and gene activation in Neurospora. PLoS Genet. 2015;11(3):e1005105. https://doi.org/10.1371/journal.pgen.1005105.
Consortium EP, Birney E, Stamatoyannopoulos JA, Dutta A, Guigo R, Gingeras TR, et al. Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project. Nature. 2007;447(7146):799–816. https://doi.org/10.1038/nature05874.
Kapranov P, Cheng J, Dike S, Nix DA, Duttagupta R, Willingham AT, et al. RNA maps reveal new RNA classes and a possible function for pervasive transcription. Science. 2007;316(5830):1484–8. https://doi.org/10.1126/science.1138341.
Kim TK, Hemberg M, Gray JM, Costa AM, Bear DM, Wu J, et al. Widespread transcription at neuronal activity-regulated enhancers. Nature. 2010;465(7295):182–7. https://doi.org/10.1038/nature09033.
Duttke SH, Chang MW, Heinz S, Benner C. Identification and dynamic quantification of regulatory elements using total RNA. Genome Res. 2019;29(11):1836–46. https://doi.org/10.1101/gr.253492.119.
Spracklin G, Pradhan S. Protect-seq: genome-wide profiling of nuclease inaccessible domains reveals physical properties of chromatin. Nucleic Acids Res. 2020;48(3):e16. https://doi.org/10.1093/nar/gkz1150.
Krogan NJ, Kim M, Tong A, Golshani A, Cagney G, Canadien V, et al. Methylation of histone H3 by Set2 in Saccharomyces cerevisiae is linked to transcriptional elongation by RNA polymerase II. Mol Cell Biol. 2003;23(12):4207–18. https://doi.org/10.1128/MCB.23.12.4207-4218.2003.
Li B, Howe L, Anderson S, Yates JR 3rd, Workman JL. The Set2 histone methyltransferase functions through the phosphorylated carboxyl-terminal domain of RNA polymerase II. J Biol Chem. 2003;278(11):8897–903. https://doi.org/10.1074/jbc.M212134200.
Li J, Moazed D, Gygi SP. Association of the histone methyltransferase Set2 with RNA polymerase II plays a role in transcription elongation. J Biol Chem. 2002;277(51):49383–8. https://doi.org/10.1074/jbc.M209294200.
Schaft D, Roguev A, Kotovic KM, Shevchenko A, Sarov M, Shevchenko A, et al. The histone 3 lysine 36 methyltransferase, SET2, is involved in transcriptional elongation. Nucleic Acids Res. 2003;31(10):2475–82. https://doi.org/10.1093/nar/gkg372.
Smolle M, Venkatesh S, Gogol MM, Li H, Zhang Y, Florens L, et al. Chromatin remodelers Isw1 and Chd1 maintain chromatin structure during transcription by preventing histone exchange. Nat Struct Mol Biol. 2012;19(9):884–92. https://doi.org/10.1038/nsmb.2312.
Maltby VE, Martin BJ, Schulze JM, Johnson I, Hentrich T, Sharma A, et al. Histone H3 lysine 36 methylation targets the Isw1b remodeling complex to chromatin. Mol Cell Biol. 2012;32(17):3479–85. https://doi.org/10.1128/MCB.00389-12.
Venkatesh S, Smolle M, Li H, Gogol MM, Saint M, Kumar S, et al. Set2 methylation of histone H3 lysine 36 suppresses histone exchange on transcribed genes. Nature. 2012;489(7416):452–5. https://doi.org/10.1038/nature11326.
Keogh MC, Kurdistani SK, Morris SA, Ahn SH, Podolny V, Collins SR, et al. Cotranscriptional set2 methylation of histone H3 lysine 36 recruits a repressive Rpd3 complex. Cell. 2005;123(4):593–605. https://doi.org/10.1016/j.cell.2005.10.025.
Carrozza MJ, Li B, Florens L, Suganuma T, Swanson SK, Lee KK, et al. Histone H3 methylation by Set2 directs deacetylation of coding regions by Rpd3S to suppress spurious intragenic transcription. Cell. 2005;123(4):581–92. https://doi.org/10.1016/j.cell.2005.10.023.
Adhvaryu KK, Morris SA, Strahl BD, Selker EU. Methylation of histone H3 lysine 36 is required for normal development in Neurospora crassa. Eukaryot Cell. 2005;4(8):1455–64. https://doi.org/10.1128/EC.4.8.1455-1464.2005.
Janevska S, Baumann L, Sieber CMK, Munsterkotter M, Ulrich J, Kamper J, et al. Elucidation of the two H3K36me3 histone methyltransferases Set2 and Ash1 in Fusarium fujikuroi unravels their different chromosomal targets and a major impact of Ash1 on genome stability. Genetics. 2018;208(1):153–71. https://doi.org/10.1534/genetics.117.1119.
Bicocca VT, Ormsby T, Adhvaryu KK, Honda S, Selker EU. ASH1-catalyzed H3K36 methylation drives gene repression and marks H3K27me2/3-competent chromatin. Elife. 2018;7. https://doi.org/10.7554/eLife.41497.
Thorvaldsdottir H, Robinson JT, Mesirov JP. Integrative genomics viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinform. 2013;14(2):178–92. https://doi.org/10.1093/bib/bbs017.
Sasaki T, Lynch KL, Mueller CV, Friedman S, Freitag M. Heterochromatin controls γH2A localization in Neurospora crassa. Eukaryot Cell. 2014;13(8):990–1000. https://doi.org/10.1128/EC.00117-14.
Borkovich KA, Alex LA, Yarden O, Freitag M, Turner GE, Read ND, et al. Lessons from the genome sequence of Neurospora crassa: tracing the path from genomic blueprint to multicellular organism. Microbiol Mol Biol Rev. 2004;68(1):1–108. https://doi.org/10.1128/MMBR.68.1.1-108.2004.
Dion MF, Kaplan T, Kim M, Buratowski S, Friedman N, Rando OJ. Dynamics of replication-independent histone turnover in budding yeast. Science. 2007;315(5817):1405–8. https://doi.org/10.1126/science.1134053.
Choi ES, Shin JA, Kim HS, Jang YK. Dynamic regulation of replication independent deposition of histone H3 in fission yeast. Nucleic Acids Res. 2005;33(22):7102–10. https://doi.org/10.1093/nar/gki1011.
Storck WK, Abdulla SZ, Rountree MR, Bicocca VT, Selker EU. A light-inducible strain for genome-wide histone turnover profiling in Neurospora crassa. Genetics. 2020;215(3):569–78. https://doi.org/10.1534/genetics.120.303217.
Smith KM, Sancar G, Dekhang R, Sullivan CM, Li S, Tag AG, et al. Transcription factors in light and circadian clock signaling networks revealed by genomewide mapping of direct targets for neurospora white collar complex. Eukaryot Cell. 2010;9(10):1549–56. https://doi.org/10.1128/EC.00154-10.
Andersson R, Gebhard C, Miguel-Escalada I, Hoof I, Bornholdt J, Boyd M, et al. An atlas of active enhancers across human cell types and tissues. Nature. 2014;507(7493):455–61. https://doi.org/10.1038/nature12787.
Core LJ, Martins AL, Danko CG, Waters CT, Siepel A, Lis JT. Analysis of nascent RNA identifies a unified architecture of initiation regions at mammalian promoters and enhancers. Nat Genet. 2014;46(12):1311–20. https://doi.org/10.1038/ng.3142.
Hu G, Cui K, Northrup D, Liu C, Wang C, Tang Q, et al. H2A.Z facilitates access of active and repressive complexes to chromatin in embryonic stem cell self-renewal and differentiation. Cell Stem Cell. 2013;12(2):180–92. https://doi.org/10.1016/j.stem.2012.11.003.
Kahn TG, Schwartz YB, Dellino GI, Pirrotta V. Polycomb complexes and the propagation of the methylation mark at the Drosophila ubx gene. J Biol Chem. 2006;281(39):29064–75. https://doi.org/10.1074/jbc.M605430200.
Alhaj Abed J, Ghotbi E, Ye P, Frolov A, Benes J, Jones RS. De novo recruitment of Polycomb-group proteins in Drosophila embryos. Development. 2018;145(23):dev165027.
Jamieson K, McNaught KJ, Ormsby T, Leggett NA, Honda S, Selker EU. Telomere repeats induce domains of H3K27 methylation in Neurospora. Elife. 2018;7. https://doi.org/10.7554/eLife.31216.
Grau DJ, Chapman BA, Garlick JD, Borowsky M, Francis NJ, Kingston RE. Compaction of chromatin by diverse Polycomb group proteins requires localized regions of high charge. Genes Dev. 2011;25(20):2210–21. https://doi.org/10.1101/gad.17288211.
Plys AJ, Davis CP, Kim J, Rizki G, Keenen MM, Marr SK, et al. Phase separation of Polycomb-repressive complex 1 is governed by a charged disordered region of CBX2. Genes Dev. 2019;33(13–14):799–813. https://doi.org/10.1101/gad.326488.119.
Schwaiger M, Schonauer A, Rendeiro AF, Pribitzer C, Schauer A, Gilles AF, et al. Evolutionary conservation of the eumetazoan gene regulatory landscape. Genome Res. 2014;24(4):639–50. https://doi.org/10.1101/gr.162529.113.
Andersson R, Sandelin A. Determinants of enhancer and promoter activities of regulatory elements. Nat Rev Genet. 2020;21(2):71–87. https://doi.org/10.1038/s41576-019-0173-8.
McNaught KJ, Wiles ET, Selker EU. Identification of a PRC2 Accessory Subunit Required for Subtelomeric H3K27 Methylation in Neurospora crassa. Mol Cell Biol. 2020;40:e00003–20.
Sun G, Zhou Z, Liu X, Gai K, Liu Q, Cha J, et al. Suppression of WHITE COLLAR-independent frequency transcription by histone H3 lysine 36 methyltransferase SET-2 is necessary for clock function in Neurospora. J Biol Chem. 2016;291(21):11055–63. https://doi.org/10.1074/jbc.M115.711333.
Xue Z, Ye Q, Anson SR, Yang J, Xiao G, Kowbel D, et al. Transcriptional interference by antisense RNA is required for circadian clock function. Nature. 2014;514(7524):650–3. https://doi.org/10.1038/nature13671.
Reyes-Dominguez Y, Bok JW, Berger H, Shwab EK, Basheer A, Gallmetzer A, et al. Heterochromatic marks are associated with the repression of secondary metabolism clusters in Aspergillus nidulans. Mol Microbiol. 2010;76(6):1376–86. https://doi.org/10.1111/j.1365-2958.2010.07051.x.
Gacek-Matthews A, Berger H, Sasaki T, Wittstein K, Gruber C, Lewis ZA, et al. KdmB, a Jumonji histone H3 demethylase, regulates genome-wide H3K4 Trimethylation and is required for Normal induction of secondary metabolism in Aspergillus nidulans. PLoS Genet. 2016;12(8):e1006222. https://doi.org/10.1371/journal.pgen.1006222.
DeGennaro CM, Alver BH, Marguerat S, Stepanova E, Davis CP, Bahler J, et al. Spt6 regulates intragenic and antisense transcription, nucleosome positioning, and histone modifications genome-wide in fission yeast. Mol Cell Biol. 2013;33(24):4779–92. https://doi.org/10.1128/MCB.01068-13.
Suzuki S, Kato H, Suzuki Y, Chikashige Y, Hiraoka Y, Kimura H, et al. Histone H3K36 trimethylation is essential for multiple silencing mechanisms in fission yeast. Nucleic Acids Res. 2016;44(9):4147–62. https://doi.org/10.1093/nar/gkw008.
Colot HV, Park G, Turner GE, Ringelberg C, Crew CM, Litvinkova L, et al. A high-throughput gene knockout procedure for Neurospora reveals functions for multiple transcription factors. Proc Natl Acad Sci U S A. 2006;103(27):10352–7. https://doi.org/10.1073/pnas.0601456103.
Davis RH, de Serres FJ.  genetic and microbiological research techniques for Neurospora crassa. Methods Enzymol. 1970;17:79–143. https://doi.org/10.1016/0076-6879(71)17168-6.
Lamb TM, Vickery J, Bell-Pedersen D. Regulation of gene expression in Neurospora crassa with a copper responsive promoter. G3 (Bethesda). 2013;3(12):2273–80.
Honda S, Selker EU. Tools for fungal proteomics: multifunctional neurospora vectors for gene replacement, protein expression and protein purification. Genetics. 2009;182(1):11–23. https://doi.org/10.1534/genetics.108.098707.
Bardiya N, Shiu PK. Cyclosporin A-resistance based gene placement system for Neurospora crassa. Fungal Genet Biol. 2007;44(5):307–14. https://doi.org/10.1016/j.fgb.2006.12.011.
Seymour M, Ji L, Santos AM, Kamei M, Sasaki T, Basenko EY, et al. Histone H1 Limits DNA Methylation in Neurospora crassa. G3 (Bethesda). 2016;6(7):1879–89.
Ferraro AR, Lewis ZA. ChIP-Seq analysis in Neurospora crassa. Methods Mol Biol. 1775;2018:241–50.
Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9(4):357–9. https://doi.org/10.1038/nmeth.1923.
Lawrence M, Huber W, Pages H, Aboyoun P, Carlson M, Gentleman R, et al. Software for computing and annotating genomic ranges. PLoS Comput Biol. 2013;9(8):e1003118. https://doi.org/10.1371/journal.pcbi.1003118.
Team RC: R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2014. URL http://wwwR-projectorg/ 2014
Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26(6):841–2. https://doi.org/10.1093/bioinformatics/btq033.
Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, et al. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008;9(9):R137. https://doi.org/10.1186/gb-2008-9-9-r137.
Ruepp A, Zollner A, Maier D, Albermann K, Hani J, Mokrejs M, et al. The FunCat, a functional annotation scheme for systematic classification of proteins from whole genomes. Nucleic Acids Res. 2004;32(18):5539–45. https://doi.org/10.1093/nar/gkh894.
Ramirez F, Ryan DP, Gruning B, Bhardwaj V, Kilpert F, Richter AS, et al. DeepTools2: a next generation web server for deep-sequencing data analysis. Nucleic Acids Res. 2016;44(W1):W160–5. https://doi.org/10.1093/nar/gkw257.
Yu G, Wang LG, He QY. ChIPseeker: an R/Bioconductor package for ChIP peak annotation, comparison and visualization. Bioinformatics. 2015;31(14):2382–3. https://doi.org/10.1093/bioinformatics/btv145.
Zhu LJ, Gazin C, Lawson ND, Pages H, Lin SM, Lapointe DS, et al. ChIPpeakAnno: a Bioconductor package to annotate ChIP-seq and ChIP-chip data. BMC Bioinformatics. 2010;11(1):237. https://doi.org/10.1186/1471-2105-11-237.
Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12(4):357–60. https://doi.org/10.1038/nmeth.3317.
Liao Y, Smyth GK, Shi W. The subread aligner: fast, accurate and scalable read mapping by seed-and-vote. Nucleic Acids Res. 2013;41(10):e108. https://doi.org/10.1093/nar/gkt214.
Heinz S, Benner C, Spann N, Bertolino E, Lin YC, Laslo P, et al. Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol Cell. 2010;38(4):576–89. https://doi.org/10.1016/j.molcel.2010.05.004.
We thank members of the UGA Chromatin Group for constructive feedback on this project.
This work was funded by grants to ZAL from the National Institute of General Medical Sciences (R01GM132644) and the American Cancer Society (RSG-14-184-01-DMC). ARF and AJA were supported by NSF GRFP awards (DGE-1443117). RJS and ZL were supported by funds from the National Science Foundation (IOS-1856627) and by the UGA Office of Research.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Sorting nuclei does not fully remove mitochondrial contamination. (A) Bar chart of mapped reads per chromosome. Sorted and Unsorted nuclei both exhibit reads mapped to mitochondrial contigs (chromosomes beginning with KC or KI). (B) Heatmap of ATAC-seq enrichment surrounding ACRs called from unsorted nuclei. Patterns of enrichment do not change in sorted versus unsorted nuclei. (C) Representative browser shot of the left arm for LG I showing overlap of sorted and unsorted nuclei ATAC-seq experiments.
. Comparison of ATAC-seq experiments in wild type and Δwc-2. Heatmap showing ATAC-seq enrichment of all ACRs in wild type and Δwc-2. Heatmap is centered on wild type ACRs.
. Motif analysis results. Both known and de novo motif analyses of ACRs within H3K27me2/3 regions as determined by HOMER.
About this article
Cite this article
Ferraro, A.R., Ameri, A.J., Lu, Z. et al. Chromatin accessibility profiling in Neurospora crassa reveals molecular features associated with accessible and inaccessible chromatin. BMC Genomics 22, 459 (2021). https://doi.org/10.1186/s12864-021-07774-0
- Chromatin accessibility