The SAGA histone acetyltransferase module targets SMC5/6 to specific genes

BACKGROUND

Structural Maintenance of Chromosomes (SMC) complexes are molecular machines driving chromatin organization at higher levels. In eukaryotes, three SMC complexes (cohesin, condensin and SMC5/6) play key roles in cohesion, condensation, replication, transcription and DNA repair. Their physical binding to DNA requires accessible chromatin.

RESULTS

We performed a genetic screen in fission yeast to identify novel factors required for SMC5/6 binding to DNA. We identified 79 genes of which histone acetyltransferases (HATs) were the most represented. Genetic and phenotypic analyses suggested a particularly strong functional relationship between the SMC5/6 and SAGA complexes. Furthermore, several SMC5/6 subunits physically interacted with SAGA HAT module components Gcn5 and Ada2. As Gcn5-dependent acetylation facilitates the accessibility of chromatin to DNA-repair proteins, we first analysed the formation of DNA-damage-induced SMC5/6 foci in the Δgcn5 mutant. The SMC5/6 foci formed normally in Δgcn5, suggesting SAGA-independent SMC5/6 localization to DNA-damaged sites. Next, we used Nse4-FLAG chromatin-immunoprecipitation (ChIP-seq) analysis in unchallenged cells to assess SMC5/6 distribution. A significant portion of SMC5/6 accumulated within gene regions in wild-type cells, which was reduced in Δgcn5 and Δada2 mutants. The drop in SMC5/6 levels was also observed in gcn5-E191Q acetyltransferase-dead mutant.

CONCLUSION

Our data show genetic and physical interactions between SMC5/6 and SAGA complexes. The ChIP-seq analysis suggests that SAGA HAT module targets SMC5/6 to specific gene regions and facilitates their accessibility for SMC5/6 loading.

H3.3 contributes to chromatin accessibility and transcription factor binding at promoter-proximal regulatory elements in embryonic stem cells

BACKGROUND

The histone variant H3.3 is enriched at active regulatory elements such as promoters and enhancers in mammalian genomes. These regions are highly accessible, creating an environment that is permissive to transcription factor binding and the recruitment of transcriptional coactivators that establish a unique chromatin post-translational landscape. How H3.3 contributes to the establishment and function of chromatin states at these regions is poorly understood.

RESULTS

We perform genomic analyses of features associated with active promoter chromatin in mouse embryonic stem cells (ESCs) and find evidence of subtle yet widespread promoter dysregulation in the absence of H3.3. Loss of H3.3 results in reduced chromatin accessibility and transcription factor (TF) binding at promoters of expressed genes in ESCs. Likewise, enrichment of the transcriptional coactivator p300 and downstream histone H3 acetylation at lysine 27 (H3K27ac) is reduced at promoters in the absence of H3.3, along with reduced enrichment of the acetyl lysine reader BRD4. Despite the observed chromatin dysregulation, H3.3 KO ESCs maintain transcription from ESC-specific genes. However, upon undirected differentiation, H3.3 KO cells retain footprinting of ESC-specific TF motifs and fail to generate footprints of lineage-specific TF motifs, in line with their diminished capacity to differentiate.

CONCLUSIONS

H3.3 facilitates DNA accessibility, transcription factor binding, and histone post-translational modification at active promoters. While H3.3 is not required for maintaining transcription in ESCs, it does promote de novo transcription factor binding which may contribute to the dysregulation of cellular differentiation in the absence of H3.3.

Epigenetic remodeling of the immune landscape in cancer: therapeutic hurdles and opportunities

The tumor immune microenvironment represents a sophisticated ecosystem where various immune cell subtypes communicate with cancer cells and stromal cells. The dynamic cellular composition and functional characteristics of the immune landscape along the trajectory of cancer development greatly impact the therapeutic efficacy and clinical outcome in patients receiving systemic antitumor therapy. Mounting evidence has suggested that epigenetic mechanisms are the underpinning of many aspects of antitumor immunity and facilitate immune state transitions during differentiation, activation, inhibition, or dysfunction. Thus, targeting epigenetic modifiers to remodel the immune microenvironment holds great potential as an integral part of anticancer regimens. In this review, we summarize the epigenetic profiles and key epigenetic modifiers in individual immune cell types that define the functional coordinates of tumor permissive and non-permissive immune landscapes. We discuss the immunomodulatory roles of current and prospective epigenetic therapeutic agents, which may open new opportunities in enhancing cancer immunotherapy or overcoming existing therapeutic challenges in the management of cancer.

Imaging Chromatin Accessibility by Assay of Transposase-Accessible Chromatin with Visualization

Chromatin accessibility is one of the fundamental structures regulating genome functions including transcription and DNA repair. Recent technological advantages to analyze chromatin accessibility begun to explore the dynamics of local chromatin structures. Here I describe protocols for Assay of Transposase-Accessible Chromatin with Visualization (ATAC-see), which allows us to analyze subnuclear localization of accessible chromatin and quantify accessible chromatin at single-cell level.

Integration of ATAC-Seq and RNA-Seq identifies key genes and pathways involved in the neuroprotection of S-adenosylmethionine against perioperative neurocognitive disorder

The gene-regulatory landscape is highly dynamic in healthy and diseased brains. DNA methylation is a well-known epigenetic modification that regulates gene expression, and our previous study demonstrated that S-adenosylmethionine (SAM), a methylome modulator, was a neuroprotectant against perioperative neurocognitive disorder (PND). However, the underlying mechanism remains to be elucidated. Here, we integrated an assay for transposase-accessible chromatin by sequencing (ATAC-seq) and RNA sequencing (RNA-seq) to identify the key genes and pathways involved in the neuroprotection of SAM against PND. Our RNA-seq data demonstrated that genes involved in biological processes such as Wnt signaling, inflammatory response, transcription and long-term potentiation likely mediate the neuroprotection of SAM. Our ATAC-seq data provided comprehensive maps of chromatin accessibility changes induced by laparotomy and laparotomy + SAM treatment, and functional annotation of the regions with high variations in chromatin accessibility highlighted the role of the Wnt signaling pathway in PND pathogenesis and SAM treatment. Further motif analysis identified key transcription factors (e.g., CTCF, TFDP1, TCFL5, KLF15, ZBTB14, TFAP2E) that may participate in the neuroprotection of SAM. In conclusion, the current study provides an epigenomic perspective to understand the pathogenesis of PND and its treatment by SAM.

Profiling Dynamic Changes in DNA Accessibility During Axon Regeneration After Optic Nerve Crush in Adult Zebrafish

A time-course series utilizing assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) can be used to detect changes in accessibility of DNA regulatory elements such as promoters and enhancers over the course of regeneration. This chapter describes methods for preparing ATAC-seq libraries from isolated zebrafish retinal ganglion cells (RGCs) following optic nerve crush at selected post-injury time points. These methods have been used for identifying dynamic changes in DNA accessibility that govern successful optic nerve regeneration in zebrafish. This method may be adapted to identify changes in DNA accessibility that accompany other types of insults to RGCs or to identify changes that occur over the course of development.

Genome-Wide Identification of Open Chromatin in Plants Using MH-Seq

Functional cis-regulatory elements (CREs) act as precise transcriptional switches for fine-tuning gene transcription. Identification of CREs is critical for understanding regulatory mechanisms of gene expression associated with various biological processes in eukaryotes. It is well known that CREs reside in open chromatin that exhibits hypersensitivity to enzyme cleavage and physical shearing. Currently, high-throughput methodologies, such as DNase-seq, ATAC-seq, and FAIRE-seq, have been widely applied in mapping open chromatin in various eukaryotic genomes. More recently, differential MNase (micrococcal nuclease) treatment has been successfully employed to map open chromatin in addition to profiling nucleosome landscape in both mammalian and plant species. We have developed a MNase hypersensitivity sequencing (MH-seq) technique in plants. The MH-seq procedure includes plant nuclei fixation and purification, differential treatments of purified nuclei with MNase, specific recovery of MNase-trimmed small DNA fragments within 20~100 bp in length, and MH-seq library construction followed by Illumina sequencing and data analysis. MH-seq has been successfully applied for global identification of open chromatin in both Arabidopsis thaliana and maize. It has been proven to be an attractive alternative for profiling open chromatin. Thus, MH-seq is expected to be valuable in probing chromatin accessibility on a genome-wide scale for other plants with sequenced genomes. Moreover, MHS data allow to implement footprinting assays to unveil binding sites of transcription factors.

Profiling Accessible Chromatin and Nucleosomes in the Mammalian Genome

Genomic DNA wraps around core histones to form nucleosomes, which provides steric constraints on how transcription factors (TFs) can interact with gene regulatory sequences. It is increasingly apparent that well-positioned, accessible nucleosomes are an inherent feature of active enhancers and can facilitate cooperative TF binding, referred to as nucleosome-mediated cooperativity. Thus, profiling chromatin and nucleosome properties (accessibility, positioning, and occupancy) on the genome is crucial to understand cell-type-specific gene regulation. Here we describe a simplified protocol to profile accessible nucleosomes in the mammalian genome using low-level and high-level micrococcal nuclease (MNase) digestion followed by genome-wide sequencing.

Assay for Transposase-Accessible Chromatin Using Sequencing of Freshly Isolated Muscle Stem Cells

Actively transcribed genes harbor cis-regulatory modules with comparatively low nucleosome occupancy and few high-order structures (="open chromatin"), whereas non-transcribed genes are characterized by high nucleosome density and extensive interactions between nucleosomes (="closed chromatin"), preventing transcription factor binding. Knowledge about chromatin accessibility is crucial to understand gene regulatory networks determining cellular decisions. Several techniques are available to map chromatin accessibility, among which the Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) is one of the most popular. ATAC-seq is based on a straightforward and robust protocol but requires adjustments for different cell types. Here, we describe an optimized protocol for ATAC-seq of freshly isolated murine muscle stem cells. We provide details for the isolation of MuSC, tagmentation, library amplification, double-sided SPRI bead cleanup, and library quality assessment and give recommendations for sequencing parameters and downstream analysis. The protocol should facilitate generation of high-quality data sets of chromatin accessibility in MuSCs, even for newcomers to the field.

Single-Cell Joint Profiling of Open Chromatin and Transcriptome by Paired-Seq

Simultaneous detection of chromatin accessibility and transcription from the same cells promises to greatly facilitate the dissection of cell-type-specific gene regulatory programs in complex tissues. Paired-seq enables joint analysis of open chromatin and nuclear transcriptome from up to a million cells in parallel. It achieves ultra-high-throughput single-cell multiomics with the use of a combinatorial barcoding strategy involving sequential ligation of multiplexed DNA barcodes to chromatin DNA fragments and reverse transcription products, followed by high-throughput DNA sequencing of the resulting DNA libraries and deconvolution of single-cell multiomic maps based on cell-specific barcodes.

Epigenomic Profiling of Dorsal Root Ganglia upon Regenerative and Non-regenerative Axonal Injury

RNA sequencing (RNA-seq), chromatin immunoprecipitation sequencing (ChIP-seq), and assay for transposase-accessible chromatin sequencing (ATAC-seq) are genome-wide techniques that provide information relative to gene expression, chromatin binding sites, and chromatin accessibility, respectively. Here we describe RNA-seq, H3K9ac, H3K27ac and H3K27me3 ChIP-seq, and ATAC-seq in dorsal root ganglia (DRG) after sciatic nerve or dorsal column axotomy, to characterize the transcriptional and epigenetic signatures of DRG upon regenerative vs non-regenerative axonal lesion.

Differential Grainy head binding correlates with variation in chromatin structure and gene expression in Drosophila melanogaster

Phenotypic evolution is often caused by variation in gene expression resulting from altered gene regulatory mechanisms. Genetic variation affecting chromatin remodeling has been identified as a potential source of variable gene expression; however, the roles of specific chromatin remodeling factors remain unclear. Here, we address this knowledge gap by examining the relationship between variation in gene expression, variation in chromatin structure, and variation in binding of the pioneer factor Grainy head between imaginal wing discs of two divergent strains of Drosophila melanogaster and their F₁ hybrid. We find that (1) variation in Grainy head binding is mostly due to sequence changes that act in cis but are located outside of the canonical Grainy head binding motif, (2) variation in Grainy head binding correlates with changes in chromatin accessibility, and (3) this variation in chromatin accessibility, coupled with variation in Grainy head binding, correlates with variation in gene expression in some cases but not others. Interactions among these three molecular layers is complex, but these results suggest that genetic variation affecting the binding of pioneer factors contributes to variation in chromatin remodeling and the evolution of gene expression.

Computational Integration of HSV-1 Multi-omics Data

Functional genomics techniques based on next-generation sequencing provide new avenues for studying host responses to viral infections at multiple levels, including transcriptional and translational processes and chromatin organization. This chapter provides an overview on the computational integration of multiple types of "omics" data on lytic herpes simplex virus 1 (HSV-1) infection. It summarizes methods developed and applied in two publications that combined 4sU-seq for studying de novo transcription, ribosome profiling for investigating active translation, RNA-seq of subcellular RNA fractions for determining subcellular location of transcripts, and ATAC-seq for profiling chromatin accessibility genome-wide. These studies revealed an unprecedented disruption of transcription termination in HSV-1 infection resulting in widespread read-through transcription beyond poly(A) sites for most but not all host genes. This impacts chromatin architecture by increasing chromatin accessibility selectively in downstream regions of affected genes. In this way, computational integration of multi-omics data identified novel and unsuspected mechanisms at play in lytic HSV-1 infection.

Integrated chromatin accessibility and DNA methylation analysis to reveal the critical epigenetic modification and regulatory mechanism in gonadal differentiation of the sequentially hermaphroditic fish, Monopterus albus

BACKGROUND

Monopterus albus is a hermaphroditic and economically farmed fish that undergoes sex reversal from ovary to testis via ovotestis during gonadal development. The epigenetic changes that are associated with gonadal development in this species remain unclear.

METHODS

We produced DNA methylome, transcriptome, and chromatin accessibility maps of the key stages of gonad development: ovary, ovotestis, and testis. The expression of the key candidate genes was detected using qRT-PCR and in situ hybridization and the methylation levels were analysed using bisulphite sequencing PCR. Promoter activity and regulation were assessed using dual-luciferase reporter assays.

RESULTS

Gonadal development exhibits highly dynamic transcriptomic, DNA methylation, and chromatin accessibility changes. We found that DNA methylation status, especially of the transcription start site, was significantly negatively correlated with gene expression while chromatin accessibility exhibited no correlation with gene expression during gonadal development. The epigenetic signatures revealed many novel regulatory elements and genes involved in sex reversal, which were validated. DNA methylation detection and site mutation of plastin-2 promoter, as a candidate gene, revealed that DNA methylation could impact the binding of transcription factor dmrt1 and foxl2 through methylation and demethylation to regulate plastin-2 expression during gonadal development.

CONCLUSIONS

These data provide novel insights into epigenetic modification and help elucidate the potential molecular mechanism by which dynamic modification of DNA methylation plays a crucial role in gonadal development.

Chromatin accessibility illuminates single-cell regulatory dynamics of rice root tips

BACKGROUND

Root development and function have central roles in plant adaptation to the environment. The modification of root traits has additionally been a major driver of crop performance since the green revolution; however, the molecular underpinnings and the regulatory programmes defining root development and response to environmental stress remain largely unknown. Single-cell reconstruction of gene regulatory programmes provides an important tool to understand the cellular phenotypic variation in complex tissues and their response to endogenous and environmental stimuli. While single-cell transcriptomes of several plant organs have been elucidated, the underlying chromatin landscapes associated with cell type-specific gene expression remain largely unexplored.

RESULTS

To comprehensively delineate chromatin accessibility during root development of an important crop, we applied single-cell ATAC-seq (scATAC-seq) to 46,758 cells from rice root tips under normal and heat stress conditions. Our data revealed cell type-specific accessibility variance across most of the major cell types and allowed us to identify sets of transcription factors which associate with accessible chromatin regions (ACRs). Using root hair differentiation as a model, we demonstrate that chromatin and gene expression dynamics during cell type differentiation correlate in pseudotime analyses. In addition to developmental trajectories, we describe chromatin responses to heat and identify cell type-specific accessibility changes to this key environmental stimulus.

CONCLUSIONS

We report chromatin landscapes during rice root development at single-cell resolution. Our work provides a framework for the integrative analysis of regulatory dynamics in this important crop organ at single-cell resolution.

Global changes in chromatin accessibility and transcription in growth hormone-secreting pituitary adenoma

PURPOSE

Growth hormone-secreting pituitary adenoma (GHPA) is an insidious disease with persistent hypersecretion of growth hormone and insulin-like growth factor 1, causing increased morbidity and mortality. Previous studies have investigated the transcription of GHPA. However, the gene regulatory landscape has not been fully characterized. The objective of our study was to unravel the changes in chromatin accessibility and transcription in GHPA.

METHODS

Six patients diagnosed with GHPA in the Department of Neurosurgery at Huashan Hospital were enrolled in our study. Primary pituitary adenoma tissues and adjacent normal pituitary specimens with no morphologic abnormalities from these six patients were obtained at surgery. RNA sequencing (RNA-seq) and assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) were applied to investigate the underlying relationship between gene expression and chromatin accessibility changes in GHPA.

RESULTS

Totally, 1528 differential expression genes (DEGs) were identified by transcriptomics analyses, including 725 up-regulated and 803 down-regulated. Further, we obtained 64 significantly DEGs including 10 DEGs were elevated and 54 DEGs were negligibly expressed in tumors tissues. The up-regulated DEGs were mainly involved in terms related to synapse formation, nervous system development and secretory pathway. In parallel, 3916 increased and 2895 decreased chromatin-accessible regions were mapped by ATAC-seq. Additionally, the chromatin accessible changes were frequently located adjacent to transcription factor CTCF and Rfx2 binding site.

CONCLUSIONS

Our results are the first to demonstrate the landscape of chromatin accessibility in GHPA, which may contribute to illustrate the underlying transcriptional regulation mechanism of this disease.

DNMT and EZH2 inhibitors synergize to activate therapeutic targets in hepatocellular carcinoma

The development of more effective targeted therapies for hepatocellular carcinoma (HCC) patients due to its aggressiveness is urgently needed. DNA methyltransferase inhibitors (DNMTis) represented the first clinical breakthrough to target aberrant cancer epigenomes. However, their clinical efficacies are still limited, in part due to an "epigenetic switch" in which a large group of genes that are demethylated by DNMTi treatment remain silenced by polycomb repressive complex 2 (PRC2) occupancy. EZH2 is the member of PRC2 that catalyzes the placement of H3K27me3 marks. EZH2 overexpression is correlated with poor HCC patient survival. We tested the combination of a DNMTi (5-aza-2'-deoxycytidine, DAC) and the EZH2 inhibitor (EZH2i) GSK126 in human HCC cell lines on drug sensitivity, DNA methylation, nucleosome accessibility, and gene expression profiles. Compared with single agent treatments, all HCC cell lines studied showed increased sensitivity after receiving both drugs concomitant with prolonged anti-proliferative changes and sustained reactivation of nascently-silenced genes. The increased number of up-regulated genes after combination treatment correlated with prolonged anti-proliferation effects and increased nucleosome accessibility. Combination treatments also activate demethylated promoters that are repressed by PRC2 occupancy. Furthermore, 13-31% of genes down-regulated by DNA methylation in primary HCC tumors were reactivated through this combination treatment scheme in vitro. Finally, the combination treatment also exacerbates anti-tumor immune responses, while most of these genes were downregulated in over 50% of primary HCC tumors. We have linked the anti-tumor effects of DAC and GSK126 combination treatments to detailed epigenetic alterations in HCC cells, identified potential therapeutic targets and provided a rationale for treatment efficacy for HCC patients.

METTL3 modulates chromatin and transcription dynamics during cell fate transition

Transcriptional programming plays a key role in determining the cell state. Timely reconfiguration of chromatin structure and attenuation of pluripotent genes are required for efficient embryonic stem cell (ESC) differentiation. Here, we identify METTL3, a core N⁶-methyladenosine (m⁶A) catalyzing enzyme, as a crucial modulator of dynamic transcription and chromatin accessibility upon ESC-derived cardiac differentiation. Genome-wide analysis of chromatin-associated RNAs revealed that depletion of METTL3 failed to dramatically attenuate the transcription of pluripotent genes, as well as activate nascent cardiomyocyte-specific transcripts upon differentiation. Consistently, ATAC-seq analysis showed that loss of METTL3 markedly attenuated the dynamic alteration of chromatin accessibility at both promoters and gene bodies, resulting in reduced sensitivity of ESC chromatin structure to cardiac differentiation signal. Furthermore, we found that METTL3 negatively regulated the histone modifications H3K4me3 and H3K36me3, which are involved in METTL3-modulated dynamic chromatin architecture during cell state transition. Unexpectedly, using chromatin-associated m⁶A sequencing, we found that nuclear m⁶A underwent a dramatic increase upon differentiation, which correlates with the decrease of chromatin accessibility. Collectively, our findings reveal that METTL3 and nuclear m⁶A epitranscriptome couple with chromatin state to ensure transcriptional regulation of cell fate transition.

Construction of a transposase accessible chromatin landscape reveals chromatin state of repeat elements and potential causal variant for complex traits in pigs

Background

A comprehensive landscape of chromatin states for multiple mammalian tissues is essential for elucidating the molecular mechanism underlying regulatory variants on complex traits. However, the genome-wide chromatin accessibility has been only reported in limited tissue types in pigs.

Results

Here we report a genome-wide landscape of chromatin accessibility of 20 tissues in two female pigs at ages of 6 months using ATAC-seq, and identified 557,273 merged peaks, which greatly expanded the pig regulatory element repository. We revealed tissue-specific regulatory elements which were associated with tissue-relevant biological functions. We identified both positive and negative significant correlations between the regulatory elements and gene transcripts, which showed distinct distributions in terms of their strength and distances from corresponding genes. We investigated the presence of transposable elements (TEs) in open chromatin regions across all tissues, these included identifications of porcine endogenous retroviruses (PERVs) exhibiting high accessibility in liver and homology of porcine specific virus sequences to universally accessible transposable elements. Furthermore, we prioritized a potential causal variant for polyunsaturated fatty acid in the muscle.

Conclusions

Our data provides a novel multi-tissues accessible chromatin landscape that serve as an important resource for interpreting regulatory sequences in tissue-specific and conserved biological functions, as well as regulatory variants of loci associated with complex traits in pigs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-022-00767-3.

Inhibition of Arid1a increases stem/progenitor cell-like properties of liver cancer

ARID1A, a key subunit of the SWI/SNF chromatin remodeling complex, exhibits recurrent mutations in various types of human cancers, including liver cancer. However, the function of ARID1A in the pathogenesis of liver cancer remains controversial. Here, we demonstrate that Arid1a knockout may result in states of different cell differentiation, as indicated by single-cell RNA sequencing (scRNA-seq) analysis. Bulk RNA-seq also revealed that Arid1a deficiency upregulated these genes related to cell stemness and differentiation, but downregulated genes related to the hepatic functions. Furthermore, we confirmed that deficiency of Arid1a increased the expression of hepatic stem/progenitor cell markers, such as Cd133 and Epcam, and enhanced the self-renewal ability of cells. Mechanistic studies revealed that Arid1a loss remodeled the chromatin accessibility of some genes related to liver functions. Thus, Arid1a deficiency might contribute to cancer development by increasing the number of stem/progenitor-like cells through dysregulating the expression of these genes related to cell stemness, differentiation and liver functions.

MSR1 characterized by chromatin accessibility mediates M2 macrophage polarization to promote gastric cancer progression

BACKGROUND

M2 macrophages play an important role in gastric cancer progression and metastasis, but the underlying tumor-promoting mechanisms are largely unknown.

METHODS

The TCIA database was used to identify the infiltration profile of macrophages. Integrated ATAC-seq, RNA-seq, and single-cell RNA sequencing (scRNA-seq) data from GC samples were used for the analysis. Using ATAC-seq profiles and RNA-seq datasets, combined with cox univariate survival analysis, we identified prognosis-related differentially expressed genes (DEGs) with chromatin accessibility, which were identified as hub genes. The CIBERSORTx algorithm was utilized to estimate the relative infiltration level of M2 macrophages, and Pearson correlation analysis was performed to investigate the relationship between hub genes and M2 macrophages. Multidimensional database validations were carried out to avoid biases. The expression level and function of hub genes in the clusters of macrophages were evaluated by using scRNA-seq data. The role of hub genes in the alternative activation of macrophages and gastric cancer malignant behaviors, as well as their potential regulatory mechanism in gastric cancer progression, were further explored.

RESULTS

17,334 genes were acquired with chromatin accessibility in promoter regions by ATAC-seq. 2,714 genes were identified with both chromatin accessibility and differential expression based on the gene expression profiles (RNA-seq). By performing Cox univariate survival analysis, 171 survival-related DEGs with chromatin accessibility were identified as hub genes. Through the CIBERSORTx algorithm and Pearson correlation analysis, the gene MSR1 most associated with M2 macrophages was screened out. According to the scRNA-seq analysis, MSR1 was highly expressed in the clusters of macrophages and may be involved in regulating M2 macrophage polarization. In vitro experiments confirmed that M2 macrophage polarization and its induced malignant behavior of gastric cancer cells were inhibited by knockdown of MSR1. Furthermore, MSR1 mediated M2 macrophage polarization by regulating arginine and proline metabolism, thereby activating the AMPK/mTOR pathway to promote gastric cancer progression.

CONCLUSION

We identified a gene-MSR1-characterized by chromatin accessibility, associated with poor prognosis in gastric cancer. This gene dictates the progression of gastric cancer by facilitating M2 macrophage polarization.

Fundamental and practical approaches for single-cell ATAC-seq analysis

Assays for transposase-accessible chromatin through high-throughput sequencing (ATAC-seq) are effective tools in the study of genome-wide chromatin accessibility landscapes. With the rapid development of single-cell technology, open chromatin regions that play essential roles in epigenetic regulation have been measured at the single-cell level using single-cell ATAC-seq approaches. The application of scATAC-seq has become as popular as that of scRNA-seq. However, owing to the nature of scATAC-seq data, which are sparse and noisy, processing the data requires different methodologies and empirical experience. This review presents a practical guide for processing scATAC-seq data, from quality evaluation to downstream analysis, for various applications. In addition to the epigenomic profiling from scATAC-seq, we also discuss recent studies in which the function of non-coding variants has been investigated based on cell type-specific cis-regulatory elements and how to use the by-product genetic information obtained from scATAC-seq to infer single-cell copy number variants and trace cell lineage. We anticipate that this review will assist researchers in designing and implementing scATAC-seq assays to facilitate research in diverse fields.

Polyvinyl alcohol/gelatin hydrogels regulate cell adhesion and chromatin accessibility

Cell adhesion has a critical influence on various processes such as cancer metastasis and wound healing. Many substrates have been used for studying cell adhesion and its related biological processes, it is still highly desirable to have a simply prepared and low-cost substrate suitable for regulating cell adhesion. In this study, we produced a series of polyvinyl alcohol/gelatin hydrogels with different gelatin concentrations via dry-annealing method. Our data showed that the protein adsorbing capability was enhanced and cell adhesion area and the ratio of non-spherical cells were increased with the increment of gelatin concentration. We also observed that varying cell adhesion conditions induced by polyvinyl alcohol /gelatin hydrogels resulted in expression level changes of genes involved in mechanotransduction from extracellular matrices (ECM) to the nucleus. In particular, we detected a widespread increase in chromatin accessibility under poor cell adhesion condition. This work provides a useful hydrogel system for regulating cell adhesion and opens up new possibilities for the design of biomaterials for cell adhesion study.

Relationship between ATOH1 and tumor microenvironment in colon adenocarcinoma patients with different microsatellite instability status

Background

Colon adenocarcinoma (COAD) is one of the major varieties of malignant tumors threatening human health today. Immune checkpoint inhibitors (ICIs) have recently begun to emerge as an effective option for the treatment of COAD patients, but not all patients can benefit from ICI treatment. Previous studies have suggested that ICIs boast significant clinical effects on patients with microsatellite instability-high (MSI-H), while conversely patients with microsatellite-stable/microsatellite instability-low (MSS/MSI-L) have shown limited response.

Methods

We used ATAC-seq, RNA-seq, and mutation data from The Cancer Genome Atlas Colon adenocarcinoma (TCGA-COAD) cohort to perform multi-omics differential analysis on COAD samples with different MSI statuses, then further screened genes by additionally combining these results with survival analysis. We analyzed the effects of the screened genes on the tumor microenvironment and immunogenicity of COAD patients, and subsequently determined their influence on the efficacy of ICIs in COAD patients using a series of predictive indexes.

Results

Twelve genes were screened in the TCGA-COAD cohort, and after the combined survival analysis, we identified ATOH1 as having significant effects. ATOH1 is characterized by high chromatin accessibility, high expression, and high mutation in COAD patients in the MSI-H group. COAD patients with high ATOH1 expression are associated with a better prognosis, unique immune microenvironment, and higher efficacy in ICI treatment. Enrichment analysis showed that COAD patients with high ATOH1 expression displayed significant upregulation in their humoral immunity and other related pathways.

Conclusions

We speculate that ATOH1 may influence the efficacy of ICIs therapy in patients with COAD by affecting the immune microenvironment and immunogenicity of the tumor.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-022-02651-6.

Nuclear receptor RORγ inverse agonists/antagonists display tissue- and gene-context selectivity through distinct activities in altering chromatin accessibility and master regulator SREBP2 occupancy

The nuclear receptor RORγ is a major driver of autoimmune diseases and certain types of cancer due to its aberrant function in T helper 17 (Th17) cell differentiation and tumor cholesterol metabolism, respectively. Compound screening using the classic receptor-coactivator interaction perturbation scheme led to identification of many small-molecule modulators of RORγ(t). We report here that inverse agonists/antagonists of RORγ such as VTP-43742 derivative VTP-23 and TAK828F, which can potently inhibit the inflammatory gene program in Th17 cells, unexpectedly lack high potency in inhibiting the growth of TNBC tumor cells. In contrast, antagonists such as XY018 and GSK805 that strongly suppress tumor cell growth and survival display only modest activities in reducing Th17-related cytokine expression. Unexpectedly, we found that VTP-23 significantly induces the cholesterol biosynthesis program in TNBC cells. Our further mechanistic analyses revealed that VTP-23 enhances the local chromatin accessibility, H3K27ac mark and the cholesterol master regulator SREBP2 recruitment at the RORγ binding sites, whereas XY018 exerts the opposite activities. Yet, they display similar inhibitory effects on circadian rhythm program. Similar distinctions and contrasting activities between TAK828F and SR2211 in their effects on local chromatin structure at Il17 genes were also observed. Together, our study shows for the first-time that structurally distinct RORγ antagonists possess different or even contrasting activities in tissue/cell-specific manner. Our findings also highlight that the activities at natural chromatin are key determinants of RORγ modulators' tissue selectivity.

DeepCAGE: Incorporating Transcription Factors in Genome-wide Prediction of Chromatin Accessibility

Although computational approaches have been complementing high-throughput biological experiments for the identification of functional regions in the human genome, it remains a great challenge to systematically decipher interactions between transcription factors (TFs) and regulatory elements to achieve interpretable annotations of chromatin accessibility across diverse cellular contexts. To solve this problem, we propose DeepCAGE, a deep learning framework that integrates sequence information and binding statuses of TFs, for the accurate prediction of chromatin accessible regions at a genome-wide scale in a variety of cell types. DeepCAGE takes advantage of a densely connected deep convolutional neural network architecture to automatically learn sequence signatures of known chromatin accessible regions and then incorporates such features with expression levels and binding activities of human core TFs to predict novel chromatin accessible regions. In a series of systematic comparisons with existing methods, DeepCAGE exhibits superior performance in not only the classification but also the regression of chromatin accessibility signals. In a detailed analysis of TF activities, DeepCAGE successfully extracts novel binding motifs and measures the contribution of a TF to the regulation with respect to a specific locus in a certain cell type. When applied to whole-genome sequencing data analysis, our method successfully prioritizes putative deleterious variants underlying a human complex trait and thus provides insights into the understanding of disease-associated genetic variants. DeepCAGE can be downloaded from https://github.com/kimmo1019/DeepCAGE.

Chromatin accessibility landscape of stromal subpopulations reveals distinct metabolic and inflammatory features of porcine subcutaneous and visceral adipose tissue

Background

Fat accumulation in visceral adipose tissue (VAT) confers increased risk for metabolic disorders of obesity, whereas accumulation of subcutaneous adipose tissue (SAT) is associated with lower risk and may be protective. Previous studies have shed light on the gene expression profile differences between SAT and VAT; however, the chromatin accessibility landscape differences and how the cis-regulatory elements govern gene expression changes between SAT and VAT are unknown.

Methods

Pig were used to characterize the differences in chromatin accessibility between the two adipose depots-derived stromal vascular fractions (SVFs) using DNase-sequencing (DNase-seq). Using integrated data from DNase-seq, H3K27ac ChIP-sequencing (ChIP-seq), and RNA-sequencing (RNA-seq), we investigated how the regulatory locus complexity regulated gene expression changes between SAT and VAT and the possible impact that these changes may have on the different biological functions of these two adipose depots.

Results

SVFs form SAT and VAT (S-SVF and V-SVF) have differential chromatin accessibility landscapes. The differential DNase I hypersensitive site (DHS)-associated genes, which indicate dynamic chromatin accessibility, were mainly involved in metabolic processes and inflammatory responses. Additionally, the Krüppel-like factor family of transcription factors were enriched in the differential DHSs. Furthermore, the chromatin accessibility data were highly associated with differential gene expression as indicated using H3K27ac ChIP-seq and RNA-seq data, supporting the validity of the differential gene expression determined using DNase-seq. Moreover, by combining epigenetic and transcriptomic data, we identified two candidate genes, NR1D1 and CRYM, could be crucial to regulate distinct metabolic and inflammatory characteristics between SAT and VAT. Together, these results uncovered differences in the transcription regulatory network and enriched the mechanistic understanding of the different biological functions between SAT and VAT.

Rapid changes in chromatin structure during dedifferentiation of primary hepatocytes in vitro

Primary hepatocytes are widely used in the pharmaceutical industry to screen drug candidates for hepatotoxicity, but hepatocytes quickly dedifferentiate and lose their mature metabolic function in culture. Attempts have been made to better recapitulate the in vivo liver environment in culture, but the full spectrum of signals required to maintain hepatocyte function ex vivo remains elusive. To elucidate molecular changes that accompany, and may contribute to dedifferentiation of hepatocytes ex vivo, we performed lineage tracing and comprehensive profiling of alterations in their gene expression profiles and chromatin landscape during culture. First, using genetically tagged hepatocytes we demonstrate that expression of the fetal gene alpha-fetoprotein in cultured hepatocytes comes from cells that previously expressed the mature gene albumin, and not from a population of albumin-negative precursor cells, proving mature hepatocytes undergo true dedifferentiation in culture. Next we studied the dedifferentiation process in detail through bulk RNA-sequencing of hepatocytes cultured over an extended period. We identified three distinct phases of dedifferentiation: an early phase, where mature hepatocyte genes are rapidly downregulated in a matter of hours; a middle phase, where fetal genes are activated; and a late phase, where initially rare contaminating non-parenchymal cells proliferate, taking over the culture. Lastly, to better understand the signaling events that result in the rapid downregulation of mature genes in hepatocytes, we examined changes in chromatin accessibility in these cells during the first 24 h of culture using Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq). We find that drastic and rapid changes in chromatin accessibility occur immediately upon the start of culture. Using binding motif analysis of the areas of open chromatin sharing similar temporal profiles, we identify several candidate transcription factors potentially involved in the dedifferentiation of primary hepatocytes in culture.

FindIT2: an R/Bioconductor package to identify influential transcription factor and targets based on multi-omics data

BACKGROUND

Transcription factors (TFs) play central roles in regulating gene expression. With the rapid growth in the use of high-throughput sequencing methods, there is a need to develop a comprehensive data processing and analyzing framework for inferring influential TFs based on ChIP-seq/ATAC-seq datasets.

RESULTS

Here, we introduce FindIT2 (Find Influential TFs and Targets), an R/Bioconductor package for annotating and processing high-throughput multi-omics data. FindIT2 supports a complete framework for annotating ChIP-seq/ATAC-seq peaks, identifying TF targets by the combination of ChIP-seq and RNA-seq datasets, and inferring influential TFs based on different types of data input. Moreover, benefited from the annotation framework based on Bioconductor, FindIT2 can be applied to any species with genomic annotations, which is particularly useful for the non-model species that are less well-studied.

CONCLUSION

FindIT2 provides a user-friendly and flexible framework to generate results at different levels according to the richness of the annotation information of user's species. FindIT2 is compatible with all the operating systems and is released under Artistic-2.0 License. The source code and documents are freely available through Bioconductor ( https://bioconductor.org/packages/devel/bioc/html/FindIT2.html ).

Extensive evaluation of ATAC-seq protocols for native or formaldehyde-fixed nuclei

BACKGROUND

The "Assay for Transposase Accessible Chromatin sequencing" (ATAC-seq) is an efficient and easy to implement protocol to measure chromatin accessibility that has been widely used in multiple applications studying gene regulation. While several modifications or variants of the protocol have been published since it was first described, there has not yet been an extensive evaluation of the effects of specific protocol choices head-to-head in a consistent experimental setting. In this study, we tested multiple protocol options for major ATAC-seq components (including three reaction buffers, two reaction temperatures, two enzyme sources, and the use of either native or fixed nuclei) in a well-characterized cell line. With all possible combinations of components, we created 24 experimental conditions with four replicates for each (a total of 96 samples). In addition, we tested the 12 native conditions in a primary sample type (mouse lung tissue) with two different input amounts. Through these extensive comparisons, we were able to observe the effect of different ATAC-seq conditions on data quality and to examine the utility and potential redundancy of various quality metrics.

RESULTS

In general, native samples yielded more peaks (particularly at loci not overlapping transcription start sites) than fixed samples, and the temperature at which the enzymatic reaction was carried out had a major impact on data quality metrics for both fixed and native nuclei. However, the effect of various conditions tested was not always consistent between the native and fixed samples. For example, the Nextera and Omni buffers were largely interchangeable across all other conditions, while the THS buffer resulted in markedly different profiles in native samples. In-house and commercial enzymes performed similarly.

CONCLUSIONS

We found that the relationship between commonly used measures of library quality differed across temperature and fixation, and so evaluating multiple metrics in assessing the quality of a sample is recommended. Notably, we also found that these choices can bias the functional class of elements profiled and so we recommend evaluating several formulations in any new experiments. Finally, we hope the ATAC-seq workflow formulated in this study on crosslinked samples will help to profile archival clinical specimens.

Computational modeling of chromatin accessibility identified important epigenomic regulators

Chromatin accessibility is essential for transcriptional activation of genomic regions. It is well established that transcription factors (TFs) and histone modifications (HMs) play critical roles in chromatin accessibility regulation. However, there is a lack of studies that quantify these relationships. Here we constructed a two-layer model to predict chromatin accessibility by integrating DNA sequence, TF binding, and HM signals. By applying the model to two human cell lines (GM12878 and HepG2), we found that DNA sequences had limited power for accessibility prediction, while both TF binding and HM signals predicted chromatin accessibility with high accuracy. According to the HM model, HM features determined chromatin accessibility in a cell line shared manner, with the prediction power attributing to five core HM types. Results from the TF model indicated that chromatin accessibility was determined by a subset of informative TFs including both cell line-specific and generic TFs. The combined model of both TF and HM signals did not further improve the prediction accuracy, indicating that they provide redundant information in terms of chromatin accessibility prediction. The TFs and HM models can also distinguish the chromatin accessibility of proximal versus distal transcription start sites with high accuracy.

Single-Molecule Multikilobase-Scale Profiling of Chromatin Accessibility Using m6A-SMAC-Seq and m6A-CpG-GpC-SMAC-Seq

A hallmark feature of active cis-regulatory elements (CREs) in eukaryotes is their nucleosomal depletion and, accordingly, higher accessibility to enzymatic treatment. This property has been the basis of a number of sequencing-based assays for genome-wide identification and tracking the activity of CREs across different biological conditions, such as DNAse-seq, ATAC-seq , NOMeseq, and others. However, the fragmentation of DNA inherent to many of these assays and the limited read length of short-read sequencing platforms have so far not allowed the simultaneous measurement of the chromatin accessibility state of CREs located distally from each other. The combination of labeling accessible DNA with DNA modifications and nanopore sequencing has made it possible to develop such assays. Here, we provide a detailed protocol for carrying out the SMAC-seq assay (Single-Molecule long-read Accessible Chromatin mapping sequencing), in its m6A-SMAC-seq and m6A-CpG-GpC-SMAC-seq variants, together with methods for data processing and analysis, and discuss key experimental and analytical considerations for working with SMAC-seq datasets.

Chromatin accessibility and microRNA expression in nephron progenitor cells during kidney development

Mammalian nephrons originate from a population of nephron progenitor cells, and changes in these cells' transcriptomes contribute to the cessation of nephrogenesis, an important determinant of nephron number. To characterize microRNA (miRNA) expression and identify putative cis-regulatory regions, we collected nephron progenitor cells from mouse kidneys at embryonic day 14.5 and postnatal day zero and assayed small RNA expression and transposase-accessible chromatin. We detect expression of 1104 miRNA (114 with expression changes), and 46,374 chromatin accessible regions (2103 with changes in accessibility). Genome-wide, our data highlight processes like cellular differentiation, cell migration, extracellular matrix interactions, and developmental signaling pathways. Furthermore, they identify new candidate cis-regulatory elements for Eya1 and Pax8, both genes with a role in nephron progenitor cell differentiation. Finally, we associate expression-changing miRNAs, including let-7-5p, miR-125b-5p, miR-181a-2-3p, and miR-9-3p, with candidate cis-regulatory elements and target genes. These analyses highlight new putative cis-regulatory loci for miRNA in nephron progenitors.

Delineating chromatin accessibility re-patterning at single cell level during early stage of direct cardiac reprogramming

Direct conversion of cardiac fibroblast into induced cardiomyocytes (iCMs) by forced expression of cardiac transcription factors, such as Mef2c, Gata4, and Tbx5 (MGT), holds great promise for regenerative medicine. The process of cardiac reprogramming consists of waves of transcriptome remodelling events. However, how this transcriptome remodelling is driven by the upstream chromatin landscape alteration is still unclear. In this study, we performed single-cell ATAC-seq (Assay for Transposase-Accessible Chromatin using sequencing) on early reprogramming iCMs given the known epigenetic changes as early as day 3. This approach unveiled networks of transcription factors (TFs) involved in the early shift of chromatin accessibility during cardiac reprogramming. Combining our analysis with functional assays, we identified Smad3 to be a bimodal TF in cardiac reprogramming, a barrier in the initiation of reprogramming and a facilitator during the intermediate stage of reprogramming. Moreover, integrative analysis of scATAC-seq with scRNA-seq data led to the identification of active TFs important for iCM conversion. Finally, we discovered a global rewiring of cis-regulatory interactions of cardiac genes along the reprogramming trajectory. Collectively, our scATAC-seq study and the integrative analysis with scRNA-seq data provided valuable resources to understand the epigenomic heterogeneity and its alteration in relation to transcription changes during early stage of cardiac reprogramming.

Analysis of Enhancers and Transcriptional Networks in Thermogenic Adipocytes

Transcription factor (TF) networks orchestrate the regulation of gene programs in mammalian cells, including white and brown adipocytes. In this protocol, we outline how genomics and transcriptomics data can be integrated to infer causal TFs of a given cellular response or cell type using "Integrated analysis of Motif Activity and Gene Expression changes of transcription factors" (IMAGE). Here, we show how key regulatory TFs controlling white and brown adipocyte gene programs can be predicted from chromatin accessibility and RNA-seq data. Furthermore, we demonstrate how information about target sites and target genes of the predicted key regulators can be integrated to propose testable hypotheses regarding the role and mechanisms of TFs.

Changes in chromatin accessibility landscape and histone H3 core acetylation during valproic acid-induced differentiation of embryonic stem cells

Directed differentiation of mouse embryonic stem cells (mESCs) or induced pluripotent stem cells (iPSCs) provides powerful models to dissect the molecular mechanisms leading to the formation of specific cell lineages. Treatment with histone deacetylase inhibitors can significantly enhance the efficiency of directed differentiation. However, the mechanisms are not well understood. Here, we use CUT&RUN in combination with ATAC-seq to determine changes in both histone modifications and genome-wide chromatin accessibility following valproic acid (VPA) exposure. VPA induced a significant increase in global histone H3 acetylation (H3K56ac), a core histone modification affecting nucleosome stability, as well as enrichment at loci associated with cytoskeletal organization and cellular morphogenesis. In addition, VPA altered the levels of linker histone H1 subtypes and the total histone H1/nucleosome ratio indicative of initial differentiation events. Notably, ATAC-seq analysis revealed changes in chromatin accessibility of genes involved in regulation of CDK serine/threonine kinase activity and DNA duplex unwinding. Importantly, changes in chromatin accessibility were evident at several key genomic loci, such as the pluripotency factor Lefty, cardiac muscle troponin Tnnt2, and the homeodomain factor Hopx, which play critical roles in cardiomyocyte differentiation. Massive parallel transcription factor (TF) footprinting also indicates an increased occupancy of TFs involved in differentiation toward mesoderm and endoderm lineages and a loss of footprints of POU5F1/SOX2 pluripotency factors following VPA treatment. Our results provide the first genome-wide analysis of the chromatin landscape following VPA-induced differentiation in mESCs and provide new mechanistic insight into the intricate molecular processes that govern departure from pluripotency and early lineage commitment.

Isolation of Plant Nuclei Compatible with Microfluidic Single-nucleus ATAC-sequencing

Gene expression depends on the binding of transcription factors with DNA regulatory sequences. The level of accessibility for these sequences varies between cells and cell types. Until recently, using the Tn5 assay for transposase-accessible chromatin for sequencing (ATAC-seq) technology allowed assessing the profiles of chromatin from an entire organ or, when coupled with the isolation of nuclei tagged in specific cell types (INTACT) method, from a cell-type. Recently, ATAC-seq experiments were conducted at the level of individual plant nuclei. Applying single nuclei ATAC-seq (sNucATAC-seq) technology to thousands of individual cells revealed more finely tuned profiles of chromatin accessibility. In this manuscript, we describe a method to isolate nuclei fom plant roots and green tissues, permeabilize the nuclear membrane using detergent to allow the penetration of the Tn5 transposase, and re-suspend them in a nuclei resuspension buffer compatible with the construction of sNucATAC-seq libraries using the 10× Genomic's Chromium technology. This protocol was successfully applied on Arabidopsis thaliana and Glycine max root nuclei.

Anticancer drug resistance: An update and perspective

Driver mutations promote initiation and progression of cancer. Pharmacological treatment can inhibit the action of the mutant protein; however, drug resistance almost invariably emerges. Multiple studies revealed that cancer drug resistance is based upon a plethora of distinct mechanisms. Drug resistance mutations can occur in the same protein or in different proteins; as well as in the same pathway or in parallel pathways, bypassing the intercepted signaling. The dilemma that the clinical oncologist is facing is that not all the genomic alterations as well as alterations in the tumor microenvironment that facilitate cancer cell proliferation are known, and neither are the alterations that are likely to promote metastasis. For example, the common KRasG12C driver mutation emerges in different cancers. Most occur in NSCLC, but some occur, albeit to a lower extent, in colorectal cancer and pancreatic ductal carcinoma. The responses to KRasG12C inhibitors are variable and fall into three categories, (i) new point mutations in KRas, or multiple copies of KRAS G12C which lead to higher expression level of the mutant protein; (ii) mutations in genes other than KRAS; (iii) original cancer transitioning to other cancer(s). Resistance to adagrasib, an experimental antitumor agent exerting its cytotoxic effect as a covalent inhibitor of the G12C KRas, indicated that half of the cases present multiple KRas mutations as well as allele amplification. Redundant or parallel pathways included MET amplification; emerging driver mutations in NRAS, BRAF, MAP2K1, and RET; gene fusion events in ALK, RET, BRAF, RAF1, and FGFR3; and loss-of-function mutations in NF1 and PTEN tumor suppressors. In the current review we discuss the molecular mechanisms underlying drug resistance while focusing on those emerging to common targeted cancer drivers. We also address questions of why cancers with a common driver mutation are unlikely to evolve a common drug resistance mechanism, and whether one can predict the likely mechanisms that the tumor cell may develop. These vastly important and tantalizing questions in drug discovery, and broadly in precision medicine, are the focus of our present review. We end with our perspective, which calls for target combinations to be selected and prioritized with the help of the emerging massive compute power which enables artificial intelligence, and the increased gathering of data to overcome its insatiable needs.

Chromatin accessibility analysis identifies GSTM1 as a prognostic marker in human glioblastoma patients

BACKGROUND

Glioblastoma (GBM) is a malignant human brain tumor that has an extremely poor prognosis. Classic mutations such as IDH (isocitrate dehydrogenase) mutations, EGFR (epidermal growth factor receptor) alternations, and MGMT (O6-methylguanine-methyltransferase) promoter hypermethylation have been used to stratify patients and provide prognostic significance. Epigenetic perturbations have been demonstrated in glioblastoma tumorigenesis. Despite the genetic markers used in the management of glioblastoma patients, new biomarkers that could predict patient survival independent of known biomarkers remain to be identified.

METHODS

ATAC-seq (assay for transposase accessible chromatin followed by sequencing) and RNA-seq have been used to profile chromatin accessible regions using glioblastoma patient samples with short-survival versus long-survival. Cell viability, cell cycle, and Western blot analysis were used to characterize the cellular phenotypes and identify signaling pathways.

RESULTS

Analysis of chromatin accessibility by ATAC-seq coupled with RNA-seq methods identified the GSTM1 (glutathione S-transferase mu-1) gene, which featured higher chromatin accessibility in GBM tumors with short survival. GSTM1 was confirmed to be a significant prognostic marker to predict survival using a different GBM patient cohort. Knockdown of GSTM1 decreased cell viability, caused cell cycle arrest, and decreased the phosphorylation levels of the NF-kB (nuclear factor kappa B) p65 subunit and STAT3 (signal transducer and activator of transcription 3) (pSer727).

CONCLUSIONS

This report demonstrates the use of ATAC-seq coupled with RNA-seq to identify GSTM1 as a prognostic marker of GBM patient survival. Activation of phosphorylation levels of NF-kB p65 and STAT3 (pSer727) by GSTM1 is shown. Analysis of chromatin accessibility in patient samples could generate an independent biomarker that can be used to predict patient survival.

Enhanced chromatin accessibility contributes to X chromosome dosage compensation in mammals

BACKGROUND

Precise gene dosage of the X chromosomes is critical for normal development and cellular function. In mice, XX female somatic cells show transcriptional X chromosome upregulation of their single active X chromosome, while the other X chromosome is inactive. Moreover, the inactive X chromosome is reactivated during development in the inner cell mass and in germ cells through X chromosome reactivation, which can be studied in vitro by reprogramming of somatic cells to pluripotency. How chromatin processes and gene regulatory networks evolved to regulate X chromosome dosage in the somatic state and during X chromosome reactivation remains unclear.

RESULTS

Using genome-wide approaches, allele-specific ATAC-seq and single-cell RNA-seq, in female embryonic fibroblasts and during reprogramming to pluripotency, we show that chromatin accessibility on the upregulated mammalian active X chromosome is increased compared to autosomes. We further show that increased accessibility on the active X chromosome is erased by reprogramming, accompanied by erasure of transcriptional X chromosome upregulation and the loss of increased transcriptional burst frequency. In addition, we characterize gene regulatory networks during reprogramming and X chromosome reactivation, revealing changes in regulatory states. Our data show that ZFP42/REX1, a pluripotency-associated gene that evolved specifically in placental mammals, targets multiple X-linked genes, suggesting an evolutionary link between ZFP42/REX1, X chromosome reactivation, and pluripotency.

CONCLUSIONS

Our data reveal the existence of intrinsic compensatory mechanisms that involve modulation of chromatin accessibility to counteract X-to-Autosome gene dosage imbalances caused by evolutionary or in vitro X chromosome loss and X chromosome inactivation in mammalian cells.

Potentials of single-cell genomics in deciphering cellular phenotypes

Single-cell genomics, particularly single-cell transcriptome profiling by RNA sequencing have transformed the possibilities to relate genes to functions, structures, and eventually phenotypes. We can now observe changes in each cell's transcriptome and among its neighborhoods, interrogate the sequence of transcriptional events, and assess their influence on subsequent events. This paradigm shift in biology enables us to infer causal relationships in these events with high accuracy. Here we review the latest single-cell studies in plants that uncover how cellular phenotypes emerge as a result of the transcriptome process such as waves of expression, trajectories of development and responses to the environment, and spatial information. With an eye on the advances made in animal and human studies, we further highlight some of the needed areas for future research and development, including computational methods.

CharPlant: A De Novo Open Chromatin Region Prediction Tool for Plant Genomes

Chromatin accessibility is a highly informative structural feature for understanding gene transcription regulation, because it indicates the degree to which nuclear macromolecules such as proteins and RNAs can access chromosomal DNA. Studies have shown that chromatin accessibility is highly dynamic during stress response, stimulus response, and developmental transition. Moreover, physical access to chromosomal DNA in eukaryotes is highly cell-specific. Therefore, current technologies such as DNase-seq, ATAC-seq, and FAIRE-seq reveal only a portion of the open chromatin regions (OCRs) present in a given species. Thus, the genome-wide distribution of OCRs remains unknown. In this study, we developed a bioinformatics tool called CharPlant for the de novo prediction of OCRs in plant genomes. To develop this tool, we constructed a three-layer convolutional neural network (CNN) and subsequently trained the CNN using DNase-seq and ATAC-seq datasets of four plant species. The model simultaneously learns the sequence motifs and regulatory logics, which are jointly used to determine DNA accessibility. All of these steps are integrated into CharPlant, which can be run using a simple command line. The results of data analysis using CharPlant in this study demonstrate its prediction power and computational efficiency. To our knowledge, CharPlant is the first de novo prediction tool that can identify potential OCRs in the whole genome. The source code of CharPlant and supporting files are freely available from https://github.com/Yin-Shen/CharPlant.

Profiling of chromatin accessibility identifies transcription factor binding sites across the genome of Aspergillus species

Background

The identification of open chromatin regions and transcription factor binding sites (TFBs) is an important step in understanding the regulation of gene expression in diverse species. ATAC-seq is a technique used for such purpose by providing high-resolution measurements of chromatin accessibility revealed through integration of Tn5 transposase. However, the existence of cell walls in filamentous fungi and associated difficulty in purifying nuclei have precluded the routine application of this technique, leading to a lack of experimentally determined and computationally inferred data on the identity of genome-wide cis-regulatory elements (CREs) and TFBs. In this study, we constructed an ATAC-seq platform suitable for filamentous fungi and generated ATAC-seq libraries of Aspergillus niger and Aspergillus oryzae grown under a variety of conditions.

Results

We applied the ATAC-seq assay for filamentous fungi to delineate the syntenic orthologue and differentially changed chromatin accessibility regions among different Aspergillus species, during different culture conditions, and among specific TF-deleted strains. The syntenic orthologues of accessible regions were responsible for the conservative functions across Aspergillus species, while regions differentially changed between culture conditions and TFs mutants drove differential gene expression programs. Importantly, we suggest criteria to determine TFBs through the analysis of unbalanced cleavage of distinct TF-bound DNA strands by Tn5 transposase. Based on this criterion, we constructed data libraries of the in vivo genomic footprint of A. niger under distinct conditions, and generated a database of novel transcription factor binding motifs through comparison of footprints in TF-deleted strains. Furthermore, we validated the novel TFBs in vivo through an artificial synthetic minimal promoter system.

Conclusions

We characterized the chromatin accessibility regions of filamentous fungi species, and identified a complete TFBs map by ATAC-seq, which provides valuable data for future analyses of transcriptional regulation in filamentous fungi.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01114-0.

Global gene expression and chromatin accessibility of the peripheral nervous system in animal models of persistent pain

BACKGROUND

Efforts to understand genetic variability involved in an individual's susceptibility to chronic pain support a role for upstream regulation by epigenetic mechanisms.

METHODS

To examine the transcriptomic and epigenetic basis of chronic pain that resides in the peripheral nervous system, we used RNA-seq and ATAC-seq of the rat dorsal root ganglion (DRG) to identify novel molecular pathways associated with pain hypersensitivity in two well-studied persistent pain models induced by chronic constriction injury (CCI) of the sciatic nerve and intra-plantar injection of complete Freund's adjuvant (CFA) in rats.

RESULTS

Our RNA-seq studies identify a variety of biological process related to synapse organization, membrane potential, transmembrane transport, and ion binding. Interestingly, genes that encode transcriptional regulators were disproportionately downregulated in both models. Our ATAC-seq data provide a comprehensive map of chromatin accessibility changes in the DRG. A total of 1123 regions showed changes in chromatin accessibility in one or both models when compared to the naïve and 31 shared differentially accessible regions (DAR)s. Functional annotation of the DARs identified disparate molecular functions enriched for each pain model which suggests that chromatin structure may be altered differently following sciatic nerve injury and hind paw inflammation. Motif analysis identified 17 DNA sequences known to bind transcription factors in the CCI DARs and 33 in the CFA DARs. Two motifs were significantly enriched in both models.

CONCLUSIONS

Our improved understanding of the changes in chromatin accessibility that occur in chronic pain states may identify regulatory genomic elements that play essential roles in modulating gene expression in the DRG.

Sex dependent glial-specific changes in the chromatin accessibility landscape in late-onset Alzheimer's disease brains

BACKGROUND

In the post-GWAS era, there is an unmet need to decode the underpinning genetic etiologies of late-onset Alzheimer's disease (LOAD) and translate the associations to causation.

METHODS

We conducted ATAC-seq profiling using NeuN sorted-nuclei from 40 frozen brain tissues to determine LOAD-specific changes in chromatin accessibility landscape in a cell-type specific manner.

RESULTS

We identified 211 LOAD-specific differential chromatin accessibility sites in neuronal-nuclei, four of which overlapped with LOAD-GWAS regions (±100 kb of SNP). While the non-neuronal nuclei did not show LOAD-specific differences, stratification by sex identified 842 LOAD-specific chromatin accessibility sites in females. Seven of these sex-dependent sites in the non-neuronal samples overlapped LOAD-GWAS regions including APOE. LOAD loci were functionally validated using single-nuclei RNA-seq datasets.

CONCLUSIONS

Using brain sorted-nuclei enabled the identification of sex-dependent cell type-specific LOAD alterations in chromatin structure. These findings enhance the interpretation of LOAD-GWAS discoveries, provide potential pathomechanisms, and suggest novel LOAD-loci.

Acute depletion of CTCF rewires genome-wide chromatin accessibility

BACKGROUND

The transcription factor CTCF appears indispensable in defining topologically associated domain boundaries and maintaining chromatin loop structures within these domains, supported by numerous functional studies. However, acute depletion of CTCF globally reduces chromatin interactions but does not significantly alter transcription.

RESULTS

Here, we systematically integrate multi-omics data including ATAC-seq, RNA-seq, WGBS, Hi-C, Cut&Run, and CRISPR-Cas9 survival dropout screens, and time-solved deep proteomic and phosphoproteomic analyses in cells carrying auxin-induced degron at endogenous CTCF locus. Acute CTCF protein degradation markedly rewires genome-wide chromatin accessibility. Increased accessible chromatin regions are frequently located adjacent to CTCF-binding sites at promoter regions and insulator sites associated with enhanced transcription of nearby genes. In addition, we use CTCF-associated multi-omics data to establish a combinatorial data analysis pipeline to discover CTCF co-regulatory partners. We successfully identify 40 candidates, including multiple established partners. Interestingly, many CTCF co-regulators that have alterations of their respective downstream gene expression do not show changes of their own expression levels across the multi-omics measurements upon acute CTCF loss, highlighting the strength of our system to discover hidden co-regulatory partners associated with CTCF-mediated transcription.

CONCLUSIONS

This study highlights that CTCF loss rewires genome-wide chromatin accessibility, which plays a critical role in transcriptional regulation.

Sequencing of methylase-accessible regions in integral circular extrachromosomal DNA reveals differences in chromatin structure

BACKGROUND

Although extrachromosomal DNA (ecDNA) has been intensively studied for several decades, the mechanisms underlying its tumorigenic effects have been revealed only recently. In most conventional sequencing studies, the high-throughput short-read sequencing largely ignores the epigenetic status of most ecDNA regions except for the junctional areas.

METHODS

Here, we developed a method of sequencing enzyme-accessible chromatin in circular DNA (CCDA-seq) based on the use of methylase to label open chromatin without fragmentation and exonuclease to enrich ecDNA sequencing depth, followed by long-read nanopore sequencing.

RESULTS

Using CCDA-seq, we observed significantly different patterns in nucleosome/regulator binding to ecDNA at a single-molecule resolution.

CONCLUSIONS

These results deepen the understanding of ecDNA regulatory mechanisms.

Early adaptive chromatin remodeling events precede pathologic phenotypes and are reinforced in the failing heart

The temporal nature of chromatin structural changes underpinning pathologic transcription are poorly understood. We measured chromatin accessibility and DNA methylation to study the contribution of chromatin remodeling at different stages of cardiac hypertrophy and failure. ATAC-seq and reduced representation bisulfite sequencing were performed in cardiac myocytes after transverse aortic constriction (TAC) or depletion of the chromatin structural protein CTCF. Early compensation to pressure overload showed changes in chromatin accessibility and DNA methylation preferentially localized to intergenic and intronic regions. Most methylation and accessibility changes observed in enhancers and promoters at the late phase (3 weeks after TAC) were established at an earlier time point (3 days after TAC), before heart failure manifests. Enhancers were paired with genes based on chromatin conformation capture data: while enhancer accessibility generally correlated with changes in gene expression, this feature, nor DNA methylation, was alone sufficient to predict transcription of all enhancer interacting genes. Enrichment of transcription factors and active histone marks at these regions suggests that enhancer activity coordinates with other epigenetic factors to determine gene transcription. In support of this hypothesis, ChIP-qPCR demonstrated increased enhancer and promoter occupancy of GATA4 and NKX2.5 at Itga9 and Nppa, respectively, concomitant with increased transcription of these genes in the diseased heart. Lastly, we demonstrate that accessibility and DNA methylation are imperfect predictors of chromatin structure at the scale of A/B compartmentalization-rather, accessibility, DNA methylation, transcription factors and other histone marks work within these domains to determine gene expression. These studies establish that chromatin reorganization during early compensation after pathologic stimuli is maintained into the later decompensatory phases of heart failure. The findings reveal the rules for how local chromatin features govern gene expression in the context of global genomic structure and identify chromatin remodeling events for therapeutic targeting in disease.

AIAP: A Quality Control and Integrative Analysis Package to Improve ATAC-seq Data Analysis

Assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) is a technique widely used to investigate genome-wide chromatin accessibility. The recently published Omni-ATAC-seq protocol substantially improves the signal/noise ratio and reduces the input cell number. High-quality data are critical to ensure accurate analysis. Several tools have been developed for assessing sequencing quality and insertion size distribution for ATAC-seq data; however, key quality control (QC) metrics have not yet been established to accurately determine the quality of ATAC-seq data. Here, we optimized the analysis strategy for ATAC-seq and defined a series of QC metrics for ATAC-seq data, including reads under peak ratio (RUPr), background (BG), promoter enrichment (ProEn), subsampling enrichment (SubEn), and other measurements. We incorporated these QC tests into our recently developed ATAC-seq Integrative Analysis Package (AIAP) to provide a complete ATAC-seq analysis system, including quality assurance, improved peak calling, and downstream differential analysis. We demonstrated a significant improvement of sensitivity (20%–60%) in both peak calling and differential analysis by processing paired-end ATAC-seq datasets using AIAP. AIAP is compiled into Docker/Singularity, and it can be executed by one command line to generate a comprehensive QC report. We used ENCODE ATAC-seq data to benchmark and generate QC recommendations, and developed qATACViewer for the user-friendly interaction with the QC report. The software, source code, and documentation of AIAP are freely available at https://github.com/Zhang-lab/ATAC-seq_QC_analysis.

Understanding the DNA double-strand break repair and its therapeutic implications

Repair of DNA double-strand breaks (DSBs) and its regulation are tightly integrated inside cells. Homologous recombination, nonhomologous end joining and microhomology mediated end joining are three major DSB repair pathways in mammalian cells. Targeting proteins associated with these repair pathways using small molecule inhibitors can prove effective in tumors, especially those with deregulated repair. Sensitization of cancer to current age therapy including radio and chemotherapy, using small molecule inhibitors is promising and warrant further development. Although several are under clinical trial, till date no repair inhibitor is approved for commercial use in cancer patients, with the exception of PARP inhibitors targeting single-strand break repair. Based on molecular profiling of repair proteins, better prognostic and therapeutic output can be achieved in patients. In the present review, we highlight the different mechanisms of DSB repair, chromatin dynamics to provide repair accessibility and modulation of inhibitors in association with molecular profiling and current gold standard treatment modalities for cancer.