ATAC-seq normalization method can significantly affect differential accessibility analysis and interpretation

Targeting the 3D genome by anthracyclines for chemotherapeutic effects

The chromatins are folded into three-dimensional (3D) structures inside cells, which coordinates the regulation of gene transcription by the non-coding regulatory elements. Aberrant chromatin 3D folding has been shown in many diseases, such as acute myeloid leukemia (AML), and may contribute to tumorigenesis. The anthracycline topoisomerase II inhibitors can induce histone eviction and DNA damage. We performed genome-wide high-resolution mapping of the chemotherapeutic effects of various clinically used anthracycline drugs. ATAC-seq was used to profile the histone eviction effects of different anthracyclines. TOP2A ChIP-seq was used to profile the potential DNA damage regions. Integrated analyses show that different anthracyclines have distinct target selectivity on epigenomic regions, based on their respective ATAC-seq and ChIP-seq profiles. We identified the underlying molecular mechanism that unique anthracycline variants selectively target chromatin looping anchors via disrupting CTCF binding, suggesting an additional potential therapeutic effect on the 3D genome. We further performed Hi-C experiments, and data from K562 cells treated with the selective anthracycline drugs indicate that the 3D chromatin organization is disrupted. Furthermore, AML patients receiving anthracycline drugs showed altered chromatin structures around potential looping anchors, which linked to distinct clinical outcomes. Our data indicate that anthracyclines are potent and selective epigenomic targeting drugs and can target the 3D genome for anticancer therapy, which could be used for personalized medicine to treat tumors with aberrant 3D chromatin structures.

Large-scale discovery of chromatin dysregulation induced by oncofusions and other protein-coding variants

Population-scale databases have expanded to millions of protein-coding variants, yet insight into their mechanistic consequences has lagged. Here we present PROD-ATAC, a high-throughput method for discovering the effects of protein-coding variants on chromatin regulation. A pooled variant library is expressed in a disease-agnostic cell line, and single-cell assay for transposase-accessible chromatin resolves each variant's effect on the chromatin landscape. Using PROD-ATAC, we characterized the effects of more than 100 oncofusions (cancer-causing chimeric proteins) and controls and revealed that chromatin remodeling is common to fusions spanning an enormous range of fusion frequencies. Furthermore, fusion-induced dysregulation can be context agnostic, as observed mechanisms often overlapped with cancer and cell-type-specific prior knowledge. We also showed that gain-of-function activity is common among oncofusions. This work begins to outline a global map of fusion-induced chromatin alterations. We suggest that there might be convergent mechanisms among disparate oncofusions and shared modes of dysregulation among fusions present in tumors at different frequencies. PROD-ATAC is generalizable to any set of protein-coding variants.

Dynamic PRC1-CBX8 stabilizes a porous structure of chromatin condensates

The compaction of chromatin is a prevalent paradigm in gene repression. Chromatin compaction is commonly thought to repress transcription by restricting chromatin accessibility. However, the spatial organisation and dynamics of chromatin compacted by gene-repressing factors are unknown. Using cryo-electron tomography, we solved the three-dimensional structure of chromatin condensed by the Polycomb Repressive Complex 1 (PRC1) in a complex with CBX8. PRC1-condensed chromatin is porous and stabilised through multivalent dynamic interactions of PRC1 with chromatin. Mechanistically, positively charged residues on the internally disordered regions (IDRs) of CBX8 mask negative charges on the DNA to stabilize the condensed state of chromatin. Within condensates, PRC1 remains dynamic while maintaining a static chromatin structure. In differentiated mouse embryonic stem cells, CBX8-bound chromatin remains accessible. These findings challenge the idea of rigidly compacted polycomb domains and instead provides a mechanistic framework for dynamic and accessible PRC1-chromatin condensates.

Molecular adaptations in response to exercise training are associated with tissue-specific transcriptomic and epigenomic signatures

Regular exercise has many physical and brain health benefits, yet the molecular mechanisms mediating exercise effects across tissues remain poorly understood. Here we analyzed 400 high-quality DNA methylation, ATAC-seq, and RNA-seq datasets from eight tissues from control and endurance exercise-trained (EET) rats. Integration of baseline datasets mapped the gene location dependence of epigenetic control features and identified differing regulatory landscapes in each tissue. The transcriptional responses to 8 weeks of EET showed little overlap across tissues and predominantly comprised tissue-type enriched genes. We identified sex differences in the transcriptomic and epigenomic changes induced by EET. However, the sex-biased gene responses were linked to shared signaling pathways. We found that many G protein-coupled receptor-encoding genes are regulated by EET, suggesting a role for these receptors in mediating the molecular adaptations to training across tissues. Our findings provide new insights into the mechanisms underlying EET-induced health benefits across organs.

A disease-associated gene desert directs macrophage inflammation through ETS2

Increasing rates of autoimmune and inflammatory disease present a burgeoning threat to human health1. This is compounded by the limited efficacy of available treatments1 and high failure rates during drug development2, highlighting an urgent need to better understand disease mechanisms. Here we show how functional genomics could address this challenge. By investigating an intergenic haplotype on chr21q22-which has been independently linked to inflammatory bowel disease, ankylosing spondylitis, primary sclerosing cholangitis and Takayasu's arteritis3-6-we identify that the causal gene, ETS2, is a central regulator of human inflammatory macrophages and delineate the shared disease mechanism that amplifies ETS2 expression. Genes regulated by ETS2 were prominently expressed in diseased tissues and more enriched for inflammatory bowel disease GWAS hits than most previously described pathways. Overexpressing ETS2 in resting macrophages reproduced the inflammatory state observed in chr21q22-associated diseases, with upregulation of multiple drug targets, including TNF and IL-23. Using a database of cellular signatures7, we identified drugs that might modulate this pathway and validated the potent anti-inflammatory activity of one class of small molecules in vitro and ex vivo. Together, this illustrates the power of functional genomics, applied directly in primary human cells, to identify immune-mediated disease mechanisms and potential therapeutic opportunities.

Chromatin remodeling protein BPTF regulates transcriptional stability in planarian stem cells

Trimethylation of histone H3 lysine 4 (H3K4me3) correlates strongly with gene expression in many different organisms, yet the question of whether it plays a causal role in transcriptional activity remains unresolved. Although H3K4me3 does not directly affect chromatin accessibility, it can indirectly affect genome accessibility by recruiting the ATP-dependent chromatin remodeling complex NuRF (Nucleosome Remodeling Factor). The largest subunit of NuRF, BPTF/NURF301, binds H3K4me3 specifically and recruits the NuRF complex to loci marked by this modification. Studies have shown that the strength and duration of BPTF binding likely also depends on additional chromatin features at these loci, such as lysine acetylation and variant histone proteins. However, the exact details of this recruitment mechanism vary between studies and have largely been tested in vitro. Here, we use stem cells isolated directly from live planarian animals to investigate the role of BPTF in regulating chromatin accessibility in vivo. We find that BPTF operates at gene promoters and is most effective at facilitating transcription at genes marked by Set1-dependent H3K4me3 peaks, which are significantly broader than those added by the lysine methyltransferase MLL1/2. Moreover, BPTF is essential for planarian stem cell biology and its loss of function phenotype mimics that of Set1 knockdown. Together, these data suggest that BPTF and H3K4me3 are important mediators of both transcription and in vivo stem cell function.

Evolutionary Dynamics of Chromatin Structure and Duplicate Gene Expression in Diploid and Allopolyploid Cotton

Polyploidy is a prominent mechanism of plant speciation and adaptation, yet the mechanistic understandings of duplicated gene regulation remain elusive. Chromatin structure dynamics are suggested to govern gene regulatory control. Here, we characterized genome-wide nucleosome organization and chromatin accessibility in allotetraploid cotton, Gossypium hirsutum (AADD, 2n = 4X = 52), relative to its two diploid parents (AA or DD genome) and their synthetic diploid hybrid (AD), using DNS-seq. The larger A-genome exhibited wider average nucleosome spacing in diploids, and this intergenomic difference diminished in the allopolyploid but not hybrid. Allopolyploidization also exhibited increased accessibility at promoters genome-wide and synchronized cis-regulatory motifs between subgenomes. A prominent cis-acting control was inferred for chromatin dynamics and demonstrated by transposable element removal from promoters. Linking accessibility to gene expression patterns, we found distinct regulatory effects for hybridization and later allopolyploid stages, including nuanced establishment of homoeolog expression bias and expression level dominance. Histone gene expression and nucleosome organization are coordinated through chromatin accessibility. Our study demonstrates the capability to track high-resolution chromatin structure dynamics and reveals their role in the evolution of cis-regulatory landscapes and duplicate gene expression in polyploids, illuminating regulatory ties to subgenomic asymmetry and dominance.

RNA expression changes driven by altered epigenetics status related to NASH etiology

Non-alcoholic fatty liver disease (NAFLD) is a growing health problem due to the increased obesity rates, among other factors. In its more severe stage (NASH), inflammation, hepatocellular ballooning and fibrosis are present in the liver, which can further evolve to total liver dysfunction or even hepatocarcinoma. As a metabolic disease, is associated to environmental factors such as diet and lifestyle conditions, which in turn can influence the epigenetic landscape of the cells, affecting to the gene expression profile and chromatin organization. In this study we performed ATAC-sequencing and RNA-sequencing to interrogate the chromatin status of liver biopsies in subjects with and without NASH and its effects on RNA transcription and NASH etiology. NASH subjects showed transcriptional downregulation for lipid and glucose metabolic pathways (e.g., ABC transporters, AMPK, FoxO or insulin pathways). A total of 229 genes were differentially enriched (ATAC and mRNA) in NASH, which were mainly related to lipid transport activity, nuclear receptor-binding, dicarboxylic acid transporter, and PPARA lipid regulation. Interpolation of ATAC data with known liver enhancer regions showed differential openness at 8 enhancers, some linked to genes involved in lipid metabolism, (i.e., FASN) and glucose homeostasis (i.e., GCGR). In conclusion, the chromatin landscape is altered in NASH patients compared to patients without this liver condition. This alteration might cause mRNA changes explaining, at least partially, the etiology and pathophysiology of the disease.

Genetic variation is a key determinant of chromatin accessibility and drives differences in the regulatory landscape of C57BL/6J and 129S1/SvImJ mice

Most common genetic variants associated with disease are located in non-coding regions of the genome. One mechanism by which they function is through altering transcription factor (TF) binding. In this study, we explore how genetic variation is connected to differences in the regulatory landscape of livers from C57BL/6J and 129S1/SvImJ mice fed either chow or a high-fat diet. To identify sites where regulatory variation affects TF binding and nearby gene expression, we employed an integrative analysis of H3K27ac ChIP-seq (active enhancers), ATAC-seq (chromatin accessibility) and RNA-seq (gene expression). We show that, across all these assays, the genetically driven (i.e. strain-specific) differences in the regulatory landscape are more pronounced than those modified by diet. Most notably, our analysis revealed that differentially accessible regions (DARs, N = 29635, FDR < 0.01 and fold change > 50%) are almost always strain-specific and enriched with genetic variation. Moreover, proximal DARs are highly correlated with differentially expressed genes. We also show that TF binding is affected by genetic variation, which we validate experimentally using ChIP-seq for TCF7L2 and CTCF. This study provides detailed insights into how non-coding genetic variation alters the gene regulatory landscape, and demonstrates how this can be used to study the regulatory variation influencing TF binding.

Chromatin Profiles Are Prognostic of Clinical Response to Bortezomib-Containing Chemotherapy in Pediatric Acute Myeloid Leukemia: Results from the COG AAML1031 Trial

The addition of the proteasome inhibitor bortezomib to standard chemotherapy did not improve survival in pediatric acute myeloid leukemia (AML) when all patients were analyzed as a group in the Children's Oncology Group phase 3 trial AAML1031 (NCT01371981). Proteasome inhibition influences the chromatin landscape and proteostasis, and we hypothesized that baseline proteomic analysis of histone- and chromatin-modifying enzymes (HMEs) would identify AML subgroups that benefitted from bortezomib addition. A proteomic profile of 483 patients treated with AAML1031 chemotherapy was generated using a reverse-phase protein array. A relatively high expression of 16 HME was associated with lower EFS and higher 3-year relapse risk after AML standard treatment compared to low expressions (52% vs. 29%, p = 0.005). The high-HME profile correlated with more transposase-accessible chromatin, as demonstrated via ATAC-sequencing, and the bortezomib addition improved the 3-year overall survival compared with standard therapy (62% vs. 75%, p = 0.033). These data suggest that there are pediatric AML populations that respond well to bortezomib-containing chemotherapy.

Omics-based construction of regulatory variants can be applied to help decipher pig liver-related traits

Genetic variants can influence complex traits by altering gene expression through changes to regulatory elements. However, the genetic variants that affect the activity of regulatory elements in pigs are largely unknown, and the extent to which these variants influence gene expression and contribute to the understanding of complex phenotypes remains unclear. Here, we annotate 90,991 high-quality regulatory elements using acetylation of histone H3 on lysine 27 (H3K27ac) ChIP-seq of 292 pig livers. Combined with genome resequencing and RNA-seq data, we identify 28,425 H3K27ac quantitative trait loci (acQTLs) and 12,250 expression quantitative trait loci (eQTLs). Through the allelic imbalance analysis, we validate two causative acQTL variants in independent datasets. We observe substantial sharing of genetic controls between gene expression and H3K27ac, particularly within promoters. We infer that 46% of H3K27ac exhibit a concomitant rather than causative relationship with gene expression. By integrating GWAS, eQTLs, acQTLs, and transcription factor binding prediction, we further demonstrate their application, through metabolites dulcitol, phosphatidylcholine (PC) (16:0/16:0) and published phenotypes, in identifying likely causal variants and genes, and discovering sub-threshold GWAS loci. We provide insight into the relationship between regulatory elements and gene expression, and the genetic foundation for dissecting the molecular mechanism of phenotypes.

Remodeling of the endothelial cell transcriptional program via paracrine and DNA-binding activities of MPO

Myeloperoxidase (MPO) is an enzyme that functions in host defense. MPO is released into the vascular lumen by neutrophils during inflammation and may adhere and subsequently penetrate endothelial cells (ECs) coating vascular walls. We show that MPO enters the nucleus of ECs and binds chromatin independently of its enzymatic activity. MPO drives chromatin decondensation at its binding sites and enhances condensation at neighboring regions. It binds loci relevant for endothelial-to-mesenchymal transition (EndMT) and affects the migratory potential of ECs. Finally, MPO interacts with the RNA-binding factor ILF3 thereby affecting its relative abundance between cytoplasm and nucleus. This interaction leads to change in stability of ILF3-bound transcripts. MPO-knockout mice exhibit reduced number of ECs at scar sites following myocardial infarction, indicating reduced neovascularization. In summary, we describe a non-enzymatic role for MPO in coordinating EndMT and controlling the fate of endothelial cells through direct chromatin binding and association with co-factors.

KDM6A/UTX promotes spermatogenic gene expression across generations and is not required for male fertility†

Paternal chromatin undergoes extensive structural and epigenetic changes during mammalian spermatogenesis, producing sperm with an epigenome optimized for the transition to embryogenesis. Lysine demethylase 6a (KDM6A, also called UTX) promotes gene activation in part via demethylation of H3K27me3, a developmentally important repressive modification abundant throughout the epigenome of spermatogenic cells and sperm. We previously demonstrated increased cancer risk in genetically wild-type mice derived from a paternal germ line lacking Kdm6a (Kdm6a cKO), indicating a role for KDM6A in regulating heritable epigenetic states. However, the regulatory function of KDM6A during spermatogenesis is not known. Here, we show that Kdm6a is transiently expressed in spermatogenesis, with RNA and protein expression largely limited to late spermatogonia and early meiotic prophase. Kdm6a cKO males do not have defects in fertility or the overall progression of spermatogenesis. However, hundreds of genes are deregulated upon loss of Kdm6a in spermatogenic cells, with a strong bias toward downregulation coinciding with the time when Kdm6a is expressed. Misregulated genes encode factors involved in chromatin organization and regulation of repetitive elements, and a subset of these genes was persistently deregulated in the male germ line across two generations of offspring of Kdm6a cKO males. Genome-wide epigenetic profiling revealed broadening of H3K27me3 peaks in differentiating spermatogonia of Kdm6a cKO mice, suggesting that KDM6A demarcates H3K27me3 domains in the male germ line. Our findings highlight KDM6A as a transcriptional activator in the mammalian male germ line that is dispensable for spermatogenesis but important for safeguarding gene regulatory state intergenerationally.

Lineage-specific intolerance to oncogenic drivers restricts histological transformation

Lung adenocarcinoma (LUAD) and small cell lung cancer (SCLC) are thought to originate from different epithelial cell types in the lung. Intriguingly, LUAD can histologically transform into SCLC after treatment with targeted therapies. In this study, we designed models to follow the conversion of LUAD to SCLC and found that the barrier to histological transformation converges on tolerance to Myc, which we implicate as a lineage-specific driver of the pulmonary neuroendocrine cell. Histological transformations are frequently accompanied by activation of the Akt pathway. Manipulating this pathway permitted tolerance to Myc as an oncogenic driver, producing rare, stem-like cells that transcriptionally resemble the pulmonary basal lineage. These findings suggest that histological transformation may require the plasticity inherent to the basal stem cell, enabling tolerance to previously incompatible oncogenic driver programs.

Positive Selection Drives cis-regulatory Evolution Across the Threespine Stickleback Y Chromosome

Allele-specific gene expression evolves rapidly on heteromorphic sex chromosomes. Over time, the accumulation of mutations on the Y chromosome leads to widespread loss of gametolog expression, relative to the X chromosome. It remains unclear if expression evolution on degrading Y chromosomes is primarily driven by mutations that accumulate through processes of selective interference, or if positive selection can also favor the down-regulation of coding regions on the Y chromosome that contain deleterious mutations. Identifying the relative rates of cis-regulatory sequence evolution across Y chromosomes has been challenging due to the limited number of reference assemblies. The threespine stickleback (Gasterosteus aculeatus) Y chromosome is an excellent model to identify how regulatory mutations accumulate on Y chromosomes due to its intermediate state of divergence from the X chromosome. A large number of Y-linked gametologs still exist across 3 differently aged evolutionary strata to test these hypotheses. We found that putative enhancer regions on the Y chromosome exhibited elevated substitution rates and decreased polymorphism when compared to nonfunctional sites, like intergenic regions and synonymous sites. This suggests that many cis-regulatory regions are under positive selection on the Y chromosome. This divergence was correlated with X-biased gametolog expression, indicating the loss of expression from the Y chromosome may be favored by selection. Our findings provide evidence that Y-linked cis-regulatory regions exhibit signs of positive selection quickly after the suppression of recombination and allow comparisons with recent theoretical models that suggest the rapid divergence of regulatory regions may be favored to mask deleterious mutations on the Y chromosome.

Endogenous retroviruses shape pluripotency specification in mouse embryos

The smooth and precise transition from totipotency to pluripotency is a key process in embryonic development, generating pluripotent stem cells capable of forming all cell types. While endogenous retroviruses (ERVs) are essential for early development, their precise roles in this transition remains mysterious. Using cutting-edge genetic and biochemical techniques in mice, we identify MERVL-gag, a retroviral protein, as a crucial modulator of pluripotent factors OCT4 and SOX2 during lineage specification. MERVL-gag tightly operates with URI, a prefoldin protein that concurs with pluripotency bias in mouse blastomeres, and which is indeed required for totipotency-to-pluripotency transition. Accordingly, URI loss promotes a stable totipotent-like state and embryo arrest at 2C stage. Mechanistically, URI binds and shields OCT4 and SOX2 from proteasome degradation, while MERVL-gag displaces URI from pluripotent factor interaction, causing their degradation. Our findings reveal the symbiotic coevolution of ERVs with their host cells to ensure the smooth and timely progression of early embryo development.

The Chromatin Remodeler CHD4 Sustains Ewing Sarcoma Cell Survival by Controlling Global Chromatin Architecture

Ewing sarcoma is an aggressive cancer with a defective response to DNA damage leading to an enhanced sensitivity to genotoxic agents. Mechanistically, Ewing sarcoma is driven by the fusion transcription factor EWS-FLI1, which reprograms the tumor cell epigenome. The nucleosome remodeling and deacetylase (NuRD) complex is an important regulator of chromatin function, controlling both gene expression and DNA damage repair, and has been associated with EWS-FLI1 activity. Here, a NuRD-focused CRISPR/Cas9 inactivation screen identified the helicase CHD4 as essential for Ewing sarcoma cell proliferation. CHD4 silencing induced tumor cell death by apoptosis and abolished colony formation. Although CHD4 and NuRD colocalized with EWS-FLI1 at enhancers and super-enhancers, CHD4 promoted Ewing sarcoma cell survival not by modulating EWS-FLI1 activity and its oncogenic gene expression program but by regulating chromatin structure. CHD4 depletion led to a global increase in DNA accessibility and induction of spontaneous DNA damage, resulting in an increased susceptibility to DNA-damaging agents. CHD4 loss delayed tumor growth in vivo, increased overall survival, and combination with PARP inhibition by olaparib treatment further suppressed tumor growth. Collectively, these findings highlight the NuRD subunit CHD4 as a therapeutic target in Ewing sarcoma that can potentiate the antitumor activity of genotoxic agents.

SIGNIFICANCE

CRISPR/Cas9 screening in Ewing sarcoma identifies a dependency on CHD4, which is crucial for the maintenance of chromatin architecture to suppress DNA damage and a promising therapeutic target for DNA damage repair-deficient malignancies.

Stem cell-like reprogramming is required for leukemia-initiating activity in B-ALL

B cell acute lymphoblastic leukemia (B-ALL) is a multistep disease characterized by the hierarchical acquisition of genetic alterations. However, the question of how a primary oncogene reprograms stem cell-like properties in committed B cells and leads to a preneoplastic population remains unclear. Here, we used the PAX5::ELN oncogenic model to demonstrate a causal link between the differentiation blockade, the self-renewal, and the emergence of preleukemic stem cells (pre-LSCs). We show that PAX5::ELN disrupts the differentiation of preleukemic cells by enforcing the IL7r/JAK-STAT pathway. This disruption is associated with the induction of rare and quiescent pre-LSCs that sustain the leukemia-initiating activity, as assessed using the H2B-GFP model. Integration of transcriptomic and chromatin accessibility data reveals that those quiescent pre-LSCs lose B cell identity and reactivate an immature molecular program, reminiscent of human B-ALL chemo-resistant cells. Finally, our transcriptional regulatory network reveals the transcription factor EGR1 as a strong candidate to control quiescence/resistance of PAX5::ELN pre-LSCs as well as of blasts from human B-ALL.

Epigenetic landscape analysis reveals the significance of early reduced chromatin accessibility in osteoclastogenesis

Recently improved techniques could provide snapshots of chromatin structure generated based on chromatin accessibility. Since chromatin accessibility determines transcriptional potential, it has been attempted in a variety of cell systems. However, there has been no genome-wide analysis of chromatin accessibility for the entire murine osteoclast (OC) differentiation process. We performed an Assay for Transposase-Accessible Chromatin (ATAC)-sequencing (seq) during RANKL-induced OC differentiation and found that global chromatin accessibility decreased, especially early in OC differentiation. The global histone H3K27Ac level, an active histone modification mark, was diminished during OC differentiation by western blot and histone extract experiments. Its genomic enrichment was also reduced based on publicly available H3K27Ac chromatin immunoprecipitation (ChIP)-seq data. ATAC-seq and H3K27Ac ChIP-seq data demonstrated that RANKL induced a less accessible chromatin state during OC differentiation. Restoration of reduced H3K27Ac, presumably representing accessible states upon acetate treatment, suppresses OC differentiation by provoking immune-related gene expression. Subsequential integrative analysis of ATAC-seq, RNA-seq after acetate treatment, and H3K27Ac ChIP-seq reveals that Irf8 and its downstream targets are the most vulnerable to chromatin accessibility changes and acetate supplementation. Taken together, our study generated chromatin accessibility maps during the whole OC differentiation and suggested perturbation of chromatin accessibility might be a potential therapeutic strategy for excessive OC diseases.

MESIA: multi-epigenome sample integration approach for precise peak calling

The assay for transposase-accessible chromatin with sequencing (ATAC-seq) is the most widely used method for measuring chromatin accessibility. Researchers have included multi-sample replication in ATAC-seq experimental designs. In epigenomic analysis, researchers should measure subtle changes in the peak by considering the read depth of individual samples. It is important to determine whether the peaks of each replication have an integrative meaning for the region of interest observed during multi-sample integration. We developed multi-epigenome sample integration approach for precise peak calling (MESIA), which integrates replication with high representativeness and reproducibility in multi-sample replication and determines the optimal peak. After identifying the reproducibility between all replications, our method integrated multiple samples determined as representative replicates. MESIA detected 6.06 times more peaks, and the value of the peaks was 1.32 times higher than the previously used method. MESIA is a shell-script-based open-source code that provides researchers involved in the epigenome with comprehensive insights.

Gene-regulation modules in nonalcoholic fatty liver disease revealed by single-nucleus ATAC-seq

We investigated the progression of nonalcoholic fatty liver disease from fatty liver to steatohepatitis using single-nucleus and bulk ATAC-seq on the livers of rats fed a high-fat diet (HFD). Rats fed HFD for 4 wk developed fatty liver, and those fed HFD for 8 wk further progressed to steatohepatitis. We observed an increase in the proportion of inflammatory macrophages, consistent with the pathological progression. Utilizing machine learning, we divided global gene regulation into modules, wherein transcription factors within a module could regulate genes within the same module, reaffirming known regulatory relationships between transcription factors and biological processes. We identified core genes-central to co-expression and protein-protein interaction-for the biological processes discovered. Notably, a large part of the core genes overlapped with genes previously implicated in nonalcoholic fatty liver disease. Single-nucleus ATAC-seq, combined with data-driven statistical analysis, offers insight into in vivo global gene regulation as a combination of modules and assists in identifying core genes of relevant biological processes.

Hand2 represses non-cardiac cell fates through chromatin remodeling at cis- regulatory elements

Developmental studies have revealed the importance of the transcription factor Hand2 in cardiac development. Hand2 promotes cardiac progenitor differentiation and epithelial maturation, while repressing other tissue types. The mechanisms underlying the promotion of cardiac fates are far better understood than those underlying the repression of alternative fates. Here, we assess Hand2-dependent changes in gene expression and chromatin remodeling in cardiac progenitors of zebrafish embryos. Cell-type specific transcriptome analysis shows a dual function for Hand2 in activation of cardiac differentiation genes and repression of pronephric pathways. We identify functional cis- regulatory elements whose chromatin accessibility are increased in hand2 mutant cells. These regulatory elements associate with non-cardiac gene expression, and drive reporter gene expression in tissues associated with Hand2-repressed genes. We find that functional Hand2 is sufficient to reduce non-cardiac reporter expression in cardiac lineages. Taken together, our data support a model of Hand2-dependent coordination of transcriptional programs, not only through transcriptional activation of cardiac and epithelial maturation genes, but also through repressive chromatin remodeling at the DNA regulatory elements of non-cardiac genes.

The ENCODE Imputation Challenge: a critical assessment of methods for cross-cell type imputation of epigenomic profiles

A promising alternative to comprehensively performing genomics experiments is to, instead, perform a subset of experiments and use computational methods to impute the remainder. However, identifying the best imputation methods and what measures meaningfully evaluate performance are open questions. We address these questions by comprehensively analyzing 23 methods from the ENCODE Imputation Challenge. We find that imputation evaluations are challenging and confounded by distributional shifts from differences in data collection and processing over time, the amount of available data, and redundancy among performance measures. Our analyses suggest simple steps for overcoming these issues and promising directions for more robust research.

Parallel recovery of chromatin accessibility and gene expression dynamics from frozen human regulatory T cells

Epigenetic features such as DNA accessibility dictate transcriptional regulation in a cell type- and cell state- specific manner, and mapping this in health vs. disease in clinically relevant material is opening the door to new mechanistic insights and new targets for therapy. Assay for Transposase Accessible Chromatin Sequencing (ATAC-seq) allows chromatin accessibility profiling from low cell input, making it tractable on rare cell populations, such as regulatory T (Treg) cells. However, little is known about the compatibility of the assay with cryopreserved rare cell populations. Here we demonstrate the robustness of an ATAC-seq protocol comparing primary Treg cells recovered from fresh or cryopreserved PBMC samples, in the steady state and in response to stimulation. We extend this method to explore the feasibility of conducting simultaneous quantitation of chromatin accessibility and transcriptome from a single aliquot of 50,000 cryopreserved Treg cells. Profiling of chromatin accessibility and gene expression in parallel within the same pool of cells controls for cellular heterogeneity and is particularly beneficial when constrained by limited input material. Overall, we observed a high correlation of accessibility patterns and transcription factor dynamics between fresh and cryopreserved samples. Furthermore, highly similar transcriptomic profiles were obtained from whole cells and from the supernatants recovered from ATAC-seq reactions. We highlight the feasibility of applying these techniques to profile the epigenomic landscape of cells recovered from cryopreservation biorepositories.

TGFβ1-Induced EMT in the MCF10A Mammary Epithelial Cell Line Model Is Executed Independently of SNAIL1 and ZEB1 but Relies on JUNB-Coordinated Transcriptional Regulation

Epithelial-mesenchymal transition (EMT) fosters cancer cell invasion and metastasis, the main cause of cancer-related mortality. Growing evidence that SNAIL and ZEB transcription factors, typically portrayed as master regulators of EMT, may be dispensable for this process, led us to re-investigate its mechanistic underpinnings. For this, we used an unbiased computational approach that integrated time-resolved analyses of chromatin structure and differential gene expression, to predict transcriptional regulators of TGFβ1-inducible EMT in the MCF10A mammary epithelial cell line model. Bioinformatic analyses indicated comparatively minor contributions of SNAIL proteins and ZEB1 to TGFβ1-induced EMT, whereas the AP-1 subunit JUNB was anticipated to have a much larger impact. CRISPR/Cas9-mediated loss-of-function studies confirmed that TGFβ1-induced EMT proceeded independently of SNAIL proteins and ZEB1. In contrast, JUNB was necessary and sufficient for EMT in MCF10A cells, but not in A549 lung cancer cells, indicating cell-type-specificity of JUNB EMT-regulatory capacity. Nonetheless, the JUNB-dependence of EMT-associated transcriptional reprogramming in MCF10A cells allowed to define a gene expression signature which was regulated by TGFβ1 in diverse cellular backgrounds, showed positively correlated expression with TGFβ signaling in multiple cancer transcriptomes, and was predictive of patient survival in several cancer types. Altogether, our findings provide novel mechanistic insights into the context-dependent control of TGFβ1-driven EMT and thereby may lead to improved diagnostic and therapeutic options.

Switching defective/sucrose non-fermenting chromatin remodeling complex coordinates meiotic gene activation via promoter remodeling and Meiosin activation in female germline

In mammals, primordial germ cells (PGCs) enter meiosis and differentiate into primary oocytes in embryonic ovaries. Previously, we demonstrated that meiotic gene induction and meiotic initiation were impaired in female germline cells of conditional knockout (CKO) mice lacking the Smarcb1 (Snf5) gene, which encodes a core subunit of the switching defective/sucrose non-fermenting (SWI/SNF) complex. In this study, we classified meiotic genes expressed at lower levels in Snf5 CKO females into two groups based on promoter accessibility. The promoters of 74% of these genes showed lower accessibility in mutant mice, whereas those of the remaining genes were opened without the SWI/SNF complex. Notably, the former genes included Meiosin, which encodes a transcriptional regulator essential for meiotic gene activation. The promoters of the former and the latter genes were mainly modified with H3K27me3/bivalent and H3K4me3 histone marks, respectively. A subset of the former genes was precociously activated in female PGCs deficient in polycomb repressive complexes (PRCs). Our results point to a mechanism through which the SWI/SNF complex coordinates meiotic gene activation via the remodeling of PRC-repressed genes, including Meiosin, in female germline cells.

Normalization benchmark of ATAC-seq datasets shows the importance of accounting for GC-content effects

The assay for transposase-accessible chromatin using sequencing (ATAC-seq) allows the study of epigenetic regulation of gene expression by assessing chromatin configuration for an entire genome. Despite its popularity, there have been limited studies investigating the analytical challenges related to ATAC-seq data, with most studies leveraging tools developed for bulk transcriptome sequencing. Here, we show that GC-content effects are omnipresent in ATAC-seq datasets. Since the GC-content effects are sample specific, they can bias downstream analyses such as clustering and differential accessibility analysis. We introduce a normalization method based on smooth-quantile normalization within GC-content bins and evaluate it together with 11 different normalization procedures on 8 public ATAC-seq datasets. Accounting for GC-content effects in the normalization is crucial for common downstream ATAC-seq data analyses, improving accuracy and interpretability. Through case studies, we show that exploratory data analysis is essential to guide the choice of an appropriate normalization method for a given dataset.

AP-1 Subunit JUNB Promotes Invasive Phenotypes in Endometriosis

Endometriosis is a disease defined by the presence of abnormal endometrium at ectopic sites, causing pain and infertility in 10% of women. Mutations in the chromatin remodeling protein ARID1A (AT-rich interactive domain-containing protein 1A) have been identified in endometriosis, particularly in the more severe deep infiltrating endometriosis and ovarian endometrioma subtypes. ARID1A has been shown to regulate chromatin at binding sites of the Activator Protein 1 (AP-1) transcription factor, and AP-1 expression has been shown in multiple endometriosis models. Here, we describe a role for AP-1 subunit JUNB in promoting invasive phenotypes in endometriosis. Through a series of knockdown experiments in the 12Z endometriosis cell line, we show that JUNB expression in endometriosis promotes the expression of epithelial-to-mesenchymal transition genes co-regulated by ARID1A including transcription factors SNAI1 and SNAI2, cell adhesion molecules ICAM1 and VCAM1, and extracellular matrix remodelers LOX and LOXL2. In highly invasive ARID1A-deficient endometriotic cells, co-knockdown of JUNB is sufficient to suppress invasion. These results suggest that AP-1 plays an important role in the progression of invasive endometriosis, and that therapeutic inhibition of AP-1 could prevent the occurrence of deep infiltrating endometriosis.

Lineage-specific, fast-evolving GATA-like gene regulates zygotic gene activation to promote endoderm specification and pattern formation in the Theridiidae spider

Background

The process of early development varies across the species-rich phylum Arthropoda. Owing to the limited research strategies for dissecting lineage-specific processes of development in arthropods, little is known about the variations in early arthropod development at molecular resolution. The Theridiidae spider, Parasteatoda tepidariorum, has its genome sequenced and could potentially contribute to dissecting early embryonic processes.

Results

We present genome-wide identification of candidate genes that exhibit locally restricted expression in germ disc forming stage embryos of P. tepidariorum, based on comparative transcriptomes of isolated cells from different regions of the embryo. A subsequent pilot screen by parental RNA interference identifies three genes required for body axis formation. One of them is a GATA-like gene that has been fast evolving after duplication and divergence from a canonical GATA family gene. This gene is designated fuchi nashi (fuchi) after its knockdown phenotypes, where the cell movement toward the formation of a germ disc was reversed. fuchi expression occurs in cells outside a forming germ disc and persists in the endoderm. Transcriptome and chromatin accessibility analyses of fuchi pRNAi embryos suggest that early fuchi activity regulates chromatin state and zygotic gene activation to promote endoderm specification and pattern formation. We also show that there are many uncharacterized genes regulated by fuchi.

Conclusions

Our genome-based research using an arthropod phylogenetically distant from Drosophila identifies a lineage-specific, fast-evolving gene with key developmental roles in one of the earliest, genome-wide regulatory events, and allows for molecular exploration of the developmental variations in early arthropod embryos.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01421-0.

Confined migration induces heterochromatin formation and alters chromatin accessibility

During migration, cells often squeeze through small constrictions, requiring extensive deformation. We hypothesized that nuclear deformation associated with such confined migration could alter chromatin organization and function. By studying cells migrating through microfluidic devices that mimic interstitial spaces in vivo, we found that confined migration results in increased H3K9me3 and H3K27me3 heterochromatin marks that persist for days. This "confined migration-induced heterochromatin" (CMiH) was distinct from heterochromatin formation during migration initiation. Confined migration decreased chromatin accessibility at intergenic regions near centromeres and telomeres, suggesting heterochromatin spreading from existing sites. Consistent with the overall decrease in accessibility, global transcription was decreased during confined migration. Intriguingly, we also identified increased accessibility at promoter regions of genes linked to chromatin silencing, tumor invasion, and DNA damage response. Inhibiting CMiH reduced migration speed, suggesting that CMiH promotes confined migration. Together, our findings indicate that confined migration induces chromatin changes that regulate cell migration and other functions.

PU.1-Dependent Enhancer Inhibition Separates Tet2-Deficient Hematopoiesis from Malignant Transformation

Cytosine hypermethylation in and around DNA-binding sites of master transcription factors, including PU.1, occurs in aging hematopoietic stem cells following acquired loss-of-function mutations of DNA methyl-cytosine dioxygenase ten-eleven translocation-2 (TET2), albeit functional relevance has been unclear. We show that Tet2-deficient mouse hematopoietic stem and progenitor cells undergo malignant transformation upon compromised gene regulation through heterozygous deletion of an upstream regulatory region (UREΔ/WT) of the PU.1 gene. Although compatible with multilineage blood formation at young age, Tet2-deficient PU.1 UREΔ/WT mice develop highly penetrant, transplantable acute myeloid leukemia (AML) during aging. Leukemic stem and progenitor cells show hypermethylation at putative PU.1-binding sites, fail to activate myeloid enhancers, and are hallmarked by a signature of genes with impaired expression shared with human AML. Our study demonstrates that Tet2 and PU.1 jointly suppress leukemogenesis and uncovers a methylation-sensitive PU.1-dependent gene network as a unifying molecular vulnerability associated with AML.

SIGNIFICANCE

We identify moderately impaired PU.1 mRNA expression as a biological modality predisposing Tet2-deficient hematopoietic stem and progenitor cells to malignant transformation. Our study furthermore uncovers a methylation-sensitive PU.1 gene network as a common feature of myeloid leukemia potentially allowing for the identification of patients at risk for malignant transformation. See related commentary by Schleicher and Pietras, p. 378. This article is highlighted in the In This Issue feature, p. 369.

Comprehensive assessment of differential ChIP-seq tools guides optimal algorithm selection

Background

The analysis of chromatin binding patterns of proteins in different biological states is a main application of chromatin immunoprecipitation followed by sequencing (ChIP-seq). A large number of algorithms and computational tools for quantitative comparison of ChIP-seq datasets exist, but their performance is strongly dependent on the parameters of the biological system under investigation. Thus, a systematic assessment of available computational tools for differential ChIP-seq analysis is required to guide the optimal selection of analysis tools based on the present biological scenario.

Results

We created standardized reference datasets by in silico simulation and sub-sampling of genuine ChIP-seq data to represent different biological scenarios and binding profiles. Using these data, we evaluated the performance of 33 computational tools and approaches for differential ChIP-seq analysis. Tool performance was strongly dependent on peak size and shape as well as on the scenario of biological regulation.

Conclusions

Our analysis provides unbiased guidelines for the optimized choice of software tools in differential ChIP-seq analysis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02686-y.

Inhibition of pyrimidine biosynthesis targets protein translation in acute myeloid leukemia

The mitochondrial enzyme dihydroorotate dehydrogenase (DHODH) catalyzes one of the rate‐limiting steps in de novo pyrimidine biosynthesis, a pathway that provides essential metabolic precursors for nucleic acids, glycoproteins, and phospholipids. DHODH inhibitors (DHODHi) are clinically used for autoimmune diseases and are emerging as a novel class of anticancer agents, especially in acute myeloid leukemia (AML) where pyrimidine starvation was recently shown to reverse the characteristic differentiation block in AML cells. Herein, we show that DHODH blockade rapidly shuts down protein translation in leukemic stem cells (LSCs) and has potent and selective activity against multiple AML subtypes. Moreover, we find that ablation of CDK5, a gene that is recurrently deleted in AML and related disorders, increases the sensitivity of AML cells to DHODHi. Our studies provide important molecular insights and identify a potential biomarker for an emerging strategy to target AML.

HNF4A modulates glucocorticoid action in the liver

The glucocorticoid receptor (GR) is a nuclear receptor critical to the regulation of energy metabolism and inflammation. The actions of GR are dependent on cell type and context. Here, we demonstrate the role of liver lineage-determining factor hepatocyte nuclear factor 4A (HNF4A) in defining liver specificity of GR action. In mouse liver, the HNF4A motif lies adjacent to the glucocorticoid response element (GRE) at GR binding sites within regions of open chromatin. In the absence of HNF4A, the liver GR cistrome is remodeled, with loss and gain of GR recruitment evident. Loss of chromatin accessibility at HNF4A-marked sites associates with loss of GR binding at weak GRE motifs. GR binding and chromatin accessibility are gained at sites characterized by strong GRE motifs, which show GR recruitment in non-liver tissues. The functional importance of these HNF4A-regulated GR sites is indicated by an altered transcriptional response to glucocorticoid treatment in the Hnf4a-null liver.

Functional reprogramming of monocytes in acute and convalescent severe COVID-19 patients

Severe COVID-19 disease is associated with dysregulation of the myeloid compartment during acute infection. Survivors frequently experience long-lasting sequelae, but little is known about the eventual persistence of this immune alteration. Herein, we evaluated TLR-induced cytokine responses in a cohort of mild to critical patients during acute or convalescent phases (n = 97). In the acute phase, we observed impaired cytokine production by monocytes in the patients with the most severe COVID-19. This capacity was globally restored in convalescent patients. However, we observed increased responsiveness to TLR1/2 ligation in patients who recovered from severe disease, indicating that these cells display distinct functional properties at the different stages of the disease. In patients with acute severe COVID-19, we identified a specific transcriptomic and epigenomic state in monocytes that can account for their functional refractoriness. The molecular profile of monocytes from recovering patients was distinct and characterized by increased chromatin accessibility at activating protein 1 (AP1) and MAF loci. These results demonstrate that severe COVID-19 infection has a profound impact on the differentiation status and function of circulating monocytes, during both the acute and the convalescent phases, in a completely distinct manner. This could have important implications for our understanding of short- and long-term COVID-19–related morbidity.

Cytotoxicity-Related Gene Expression and Chromatin Accessibility Define a Subset of CD4+ T Cells That Mark Progression to Type 1 Diabetes

Type 1 diabetes in children is heralded by a preclinical phase defined by circulating autoantibodies to pancreatic islet antigens. How islet autoimmunity is initiated and then progresses to clinical diabetes remains poorly understood. Only one study has reported gene expression in specific immune cells of children at risk associated with progression to islet autoimmunity. We analyzed gene expression with RNA sequencing in CD4+ and CD8+ T cells, natural killer (NK) cells, and B cells, and chromatin accessibility by assay for transposase-accessible chromatin sequencing (ATAC-seq) in CD4+ T cells, in five genetically at risk children with islet autoantibodies who progressed to diabetes over a median of 3 years ("progressors") compared with five children matched for sex, age, and HLA-DR who had not progressed ("nonprogressors"). In progressors, differentially expressed genes (DEGs) were largely confined to CD4+ T cells and enriched for cytotoxicity-related genes/pathways. Several top-ranked DEGs were validated in a semi-independent cohort of 13 progressors and 11 nonprogressors. Flow cytometry confirmed that progression was associated with expansion of CD4+ cells with a cytotoxic phenotype. By ATAC-seq, progression was associated with reconfiguration of regulatory chromatin regions in CD4+ cells, some linked to differentially expressed cytotoxicity-related genes. Our findings suggest that cytotoxic CD4+ T cells play a role in promoting progression to type 1 diabetes.

BMAL1 drives muscle repair through control of hypoxic NAD+ regeneration in satellite cells

The process of tissue regeneration occurs in a developmentally timed manner, yet the role of circadian timing is not understood. Here, we identify a role for the adult muscle stem cell (MuSC)-autonomous clock in the control of muscle regeneration following acute ischemic injury. We observed greater muscle repair capacity following injury during the active/wake period as compared with the inactive/rest period in mice, and loss of Bmal1 within MuSCs leads to impaired muscle regeneration. We demonstrate that Bmal1 loss in MuSCs leads to reduced activated MuSC number at day 3 postinjury, indicating a failure to properly expand the myogenic precursor pool. In cultured primary myoblasts, we observed that loss of Bmal1 impairs cell proliferation in hypoxia (a condition that occurs in the first 1-3 d following tissue injury in vivo), as well as subsequent myofiber differentiation. Loss of Bmal1 in both cultured myoblasts and in vivo activated MuSCs leads to reduced glycolysis and premature activation of prodifferentiation gene transcription and epigenetic remodeling. Finally, hypoxic cell proliferation and myofiber formation in Bmal1-deficient myoblasts are restored by increasing cytosolic NAD+ Together, we identify the MuSC clock as a pivotal regulator of oxygen-dependent myoblast cell fate and muscle repair through the control of the NAD+-driven response to injury.

Dynamic transcriptome analysis reveals signatures of paradoxical effect of vemurafenib on human dermal fibroblasts

BACKGROUND

Vemurafenib (PLX4032) is one of the most frequently used treatments for late-stage melanoma patients with the BRAFV600E mutation; however, acquired resistance to the drug poses as a major challenge. It remains to be determined whether off-target effects of vemurafenib on normal stroma components could reshape the tumor microenvironment in a way that contributes to cancer progression and drug resistance.

METHODS

By using temporally-resolved RNA- and ATAC-seq, we studied the early molecular changes induced by vemurafenib in human dermal fibroblast (HDF), a main stromal component in melanoma and other tumors with high prevalence of BRAFV600 mutations.

RESULTS

Transcriptomics analyses revealed a stepwise up-regulation of proliferation signatures, together with a down-regulation of autophagy and proteolytic processes. The gene expression changes in HDF strongly correlated in an inverse way with those in BRAFV600E mutant malignant melanoma (MaMel) cell lines, consistent with the observation of a paradoxical effect of vemurafenib, leading to hyperphosphorylation of MEK1/2 and ERK1/2. The transcriptional changes in HDF were not strongly determined by alterations in chromatin accessibility; rather, an already permissive chromatin landscape seemed to facilitate the early accessibility to MAPK/ERK-regulated transcription factor binding sites. Combinatorial treatment with the MEK inhibitor trametinib did not preclude the paradoxical activation of MAPK/ERK signaling in HDF. When administered together, vemurafenib partially compensated for the reduction of cell viability and proliferation induced by trametinib. These paradoxical changes were restrained by using the third generation BRAF inhibitor PLX8394, a so-called paradox breaker compound. However, the advantageous effects on HDF during combination therapies were also lost.

CONCLUSIONS

Vemurafenib induces paradoxical changes in HDF, enabled by a permissive chromatin landscape. These changes might provide an advantage during combination therapies, by compensating for the toxicity induced in stromal cells by less specific MAPK/ERK inhibitors. Our results highlight the relevance of evaluating the effects of the drugs on non-transformed stromal components, carefully considering the implications of their administration either as mono- or combination therapies. Video Abstract.

Mapping the evolution of T cell states during response and resistance to adoptive cellular therapy

To elucidate mechanisms by which T cells eliminate leukemia, we study donor lymphocyte infusion (DLI), an established immunotherapy for relapsed leukemia. We model T cell dynamics by integrating longitudinal, multimodal data from 94,517 bone marrow-derived single T cell transcriptomes in addition to chromatin accessibility and single T cell receptor sequencing from patients undergoing DLI. We find that responsive tumors are defined by enrichment of late-differentiated T cells before DLI and rapid, durable expansion of early differentiated T cells after treatment, highly similar to "terminal" and "precursor" exhausted subsets, respectively. Resistance, in contrast, is defined by heterogeneous T cell dysfunction. Surprisingly, early differentiated T cells in responders mainly originate from pre-existing and novel clonotypes recruited to the leukemic microenvironment, rather than the infusion. Our work provides a paradigm for analyzing longitudinal single-cell profiling of scenarios beyond adoptive cell therapy and introduces Symphony, a Bayesian approach to infer regulatory circuitry underlying T cell subsets, with broad relevance to exhaustion antagonists across cancers.

Genome accessibility dynamics in response to phosphate limitation is controlled by the PHR1 family of transcription factors in Arabidopsis

As phosphorus is one of the most limiting nutrients in many natural and agricultural ecosystems, plants have evolved strategies that cope with its scarcity. Genetic approaches have facilitated the identification of several molecular elements that regulate the phosphate (Pi) starvation response (PSR) of plants, including the master regulator of the transcriptional response to phosphate starvation PHOSPHATE STARVATION RESPONSE1 (PHR1). However, the chromatin modifications underlying the plant transcriptional response to phosphate scarcity remain largely unknown. Here, we present a detailed analysis of changes in chromatin accessibility during phosphate starvation in Arabidopsis thaliana root cells. Root cells undergo a genome-wide remodeling of chromatin accessibility in response to Pi starvation that is often associated with changes in the transcription of neighboring genes. Analysis of chromatin accessibility in the phr1 phl2 double mutant revealed that the transcription factors PHR1 and PHL2 play a key role in remodeling chromatin accessibility in response to Pi limitation. We also discovered that PHR1 and PHL2 play an important role in determining chromatin accessibility and the associated transcription of many genes under optimal Pi conditions, including genes involved in the PSR. We propose that a set of transcription factors directly activated by PHR1 in Pi-starved root cells trigger a second wave of epigenetic changes required for the transcriptional activation of the complete set of low-Pi-responsive genes.

CHIPIN: ChIP-seq inter-sample normalization based on signal invariance across transcriptionally constant genes

Background

Multiple studies rely on ChIP-seq experiments to assess the effect of gene modulation and drug treatments on protein binding and chromatin structure. However, most methods commonly used for the normalization of ChIP-seq binding intensity signals across conditions, e.g., the normalization to the same number of reads, either assume a constant signal-to-noise ratio across conditions or base the estimates of correction factors on genomic regions with intrinsically different signals between conditions. Inaccurate normalization of ChIP-seq signal may, in turn, lead to erroneous biological conclusions.

Results

We developed a new R package, CHIPIN, that allows normalizing ChIP-seq signals across different conditions/samples when spike-in information is not available, but gene expression data are at hand. Our normalization technique is based on the assumption that, on average, no differences in ChIP-seq signals should be observed in the regulatory regions of genes whose expression levels are constant across samples/conditions. In addition to normalizing ChIP-seq signals, CHIPIN provides as output a number of graphs and calculates statistics allowing the user to assess the efficiency of the normalization and qualify the specificity of the antibody used. In addition to ChIP-seq, CHIPIN can be used without restriction on open chromatin ATAC-seq or DNase hypersensitivity data. We validated the CHIPIN method on several ChIP-seq data sets and documented its superior performance in comparison to several commonly used normalization techniques.

Conclusions

The CHIPIN method provides a new way for ChIP-seq signal normalization across conditions when spike-in experiments are not available. The method is implemented in a user-friendly R package available on GitHub: https://github.com/BoevaLab/CHIPIN

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04320-3.

Simultaneous disruption of PRC2 and enhancer function underlies histone H3.3-K27M oncogenic activity in human hindbrain neural stem cells

Driver mutations in genes encoding histone H3 proteins resulting in p.Lys27Met substitutions (H3-K27M) are frequent in pediatric midline brain tumors. However, the precise mechanisms by which H3-K27M causes tumor initiation remain unclear. Here, we use human hindbrain neural stem cells to model the consequences of H3.3-K27M on the epigenomic landscape in a relevant developmental context. Genome-wide mapping of epitope-tagged histone H3.3 revealed that both the wild type and the K27M mutant incorporate abundantly at pre-existing active enhancers and promoters, and to a lesser extent at Polycomb repressive complex 2 (PRC2)-bound regions. At active enhancers, H3.3-K27M leads to focal H3K27ac loss, decreased chromatin accessibility and reduced transcriptional expression of nearby neurodevelopmental genes. In addition, H3.3-K27M deposition at a subset of PRC2 target genes leads to increased PRC2 and PRC1 binding and augmented transcriptional repression that can be partially reversed by PRC2 inhibitors. Our work suggests that, rather than imposing de novo transcriptional circuits, H3.3-K27M drives tumorigenesis by locking initiating cells in their pre-existing, immature epigenomic state, via disruption of PRC2 and enhancer functions.

Differential chromatin accessibility landscape of gain-of-function mutant p53 tumours

BACKGROUND

Mutations in TP53 not only affect its tumour suppressor activity but also exerts oncogenic gain-of-function activity. While the genome-wide mutant p53 binding sites have been identified in cancer cell lines, the chromatin accessibility landscape driven by mutant p53 in primary tumours is unknown. Here, we leveraged the chromatin accessibility data of primary tumours from The Cancer Genome Atlas (TCGA) to identify differentially accessible regions in mutant p53 tumours compared to wild-type p53 tumours, especially in breast and colon cancers.

RESULTS

We identified 1587 lost and 984 gained accessible chromatin regions in breast, and 1143 lost and 640 gained regions in colon cancers. However, only less than half of those regions in both cancer types contain sequence motifs for wild-type or mutant p53 binding. Whereas, the remaining showed enrichment for master transcriptional regulators, such as FOX-Family TFs and NF-kB in lost and SMAD and KLF TFs in gained regions of breast. In colon, ATF3 and FOS/JUN TFs were enriched in lost, and CDX family TFs and HNF4A in gained regions. By integrating the gene expression data, we identified known and novel target genes regulated by the mutant p53.

CONCLUSION

This study reveals the direct and indirect mechanisms by which gain-of-function mutant p53 targets the chromatin and subsequent gene expression patterns in a tumour-type specific manner. This furthers our understanding of the impact of mutant p53 in cancer development.

ATRX limits the accessibility of histone H3-occupied HSV genomes during lytic infection

Histones are rapidly loaded on the HSV genome upon entry into the nucleus of human fibroblasts, but the effects of histone loading on viral replication have not been fully defined. We showed recently that ATRX is dispensable for de novo deposition of H3 to HSV genomes after nuclear entry but restricted infection through maintenance of viral heterochromatin. To further investigate the roles that ATRX and other histone H3 chaperones play in restriction of HSV, we infected human fibroblasts that were systematically depleted of nuclear H3 chaperones. We found that the ATRX/DAXX complex is unique among nuclear H3 chaperones in its capacity to restrict ICP0-null HSV infection. Only depletion of ATRX significantly alleviated restriction of viral replication. Interestingly, no individual nuclear H3 chaperone was required for deposition of H3 onto input viral genomes, suggesting that during lytic infection, H3 deposition may occur through multiple pathways. ChIP-seq for total histone H3 in control and ATRX-KO cells infected with ICP0-null HSV showed that HSV DNA is loaded with high levels of histones across the entire viral genome. Despite high levels of H3, ATAC-seq analysis revealed that HSV DNA is highly accessible, especially in regions of high GC content, and is not organized largely into ordered nucleosomes during lytic infection. ATRX reduced accessibility of viral DNA to the activity of a TN5 transposase and enhanced accumulation of viral DNA fragment sizes associated with nucleosome-like structures. Together, these findings support a model in which ATRX restricts viral infection by altering the structure of histone H3-loaded viral chromatin that reduces viral DNA accessibility for transcription. High GC rich regions of the HSV genome, especially the S component inverted repeats of the HSV-1 genome, show increased accessibility, which may lead to increased ability to transcribe the IE genes encoded in these regions during initiation of infection.

Protocol for assaying chromatin accessibility using ATAC-seq in plants

Open or accessible regions of the genome are the primary positions of binding sites for transcription factors and chromatin regulators. Transposase-accessible chromatin sequencing (ATAC-seq) can probe chromatin accessibility in the intact nucleus. Here, we describe a protocol to generate ATAC-seq libraries from fresh Arabidopsis thaliana tissues and establish an easy-to-use bioinformatic analysis pipeline. Our method could be applied to other plants and other tissues and allows for the reliable detection of changes in chromatin accessibility throughout plant growth and development.

For complete details on the use and execution of this protocol, please refer to Wang et al. (2020).

•

A protocol to generate ATAC-seq libraries from fresh plant tissues

•

An easy-to-use bioinformatic analysis pipeline for ATAC-seq

•

Identification of differentially accessible peaks by ATAC-seq

Chromatin accessibility profiling methods

Chromatin accessibility, or the physical access to chromatinized DNA, is a widely studied characteristic of the eukaryotic genome. As active regulatory DNA elements are generally 'accessible', the genome-wide profiling of chromatin accessibility can be used to identify candidate regulatory genomic regions in a tissue or cell type. Multiple biochemical methods have been developed to profile chromatin accessibility, both in bulk and at the single-cell level. Depending on the method, enzymatic cleavage, transposition or DNA methyltransferases are used, followed by high-throughput sequencing, providing a view of genome-wide chromatin accessibility. In this Primer, we discuss these biochemical methods, as well as bioinformatics tools for analysing and interpreting the generated data, and insights into the key regulators underlying developmental, evolutionary and disease processes. We outline standards for data quality, reproducibility and deposition used by the genomics community. Although chromatin accessibility profiling is invaluable to study gene regulation, alone it provides only a partial view of this complex process. Orthogonal assays facilitate the interpretation of accessible regions with respect to enhancer-promoter proximity, functional transcription factor binding and regulatory function. We envision that technological improvements including single-molecule, multi-omics and spatial methods will bring further insight into the secrets of genome regulation.

SWI/SNF inactivation in the endometrial epithelium leads to loss of epithelial integrity

Although ARID1A mutations are a hallmark feature, mutations in other SWI/SNF (SWItch/Sucrose Non-Fermentable) chromatin remodeling subunits are also observed in endometrial neoplasms. Here, we interrogated the roles of Brahma/SWI2-related gene 1 (BRG1, SMARCA4), the SWI/SNF catalytic subunit, in the endometrial epithelium. BRG1 loss affects more than one-third of all active genes and highly overlaps with the ARID1A gene regulatory network. Chromatin immunoprecipitation studies revealed widespread subunit-specific differences in transcriptional regulation, as BRG1 promoter interactions are associated with gene activation, while ARID1A binding is associated with gene repression. However, we identified a physiologically relevant subset of BRG1 and ARID1A co-regulated epithelial identity genes. Mice were genetically engineered to inactivate BRG1 specifically in the endometrial epithelium. Endometrial glands were observed embedded in uterine myometrium, indicating adenomyosis-like phenotypes. Molecular similarities were observed between BRG1 and ARID1A mutant endometrial cells in vivo, including loss of epithelial cell adhesion and junction genes. Collectively, these studies illustrate overlapping contributions of multiple SWI/SNF subunit mutations in the translocation of endometrium to distal sites, with loss of cell integrity being a common feature in SWI/SNF mutant endometrial epithelia.