MAnorm2 for quantitatively comparing groups of ChIP-seq samples

Suppression of GATA3 promotes epithelial-mesenchymal transition and simultaneous cellular senescence in human extravillous trophoblasts

The regulatory mechanism of the transcription factor GATA3 in the differentiation and maturation process of extravillous trophoblasts (EVT) in early pregnancy placenta, as well as its relevance to the occurrence of pregnancy disorders, remains poorly understood. This study leveraged single-cell RNA sequencing data from placental organoid models and placental tissue to explore the dynamic changes in GATA3 expression during EVT maturation. The expression pattern exhibited an initial upregulation followed by subsequent downregulation, with aberrant GATA3 localization observed in cases of recurrent miscarriage (RM). By identifying global targets regulated by GATA3 in primary placental EVT cells, JEG3, and HTR8/SVneo cell lines, this study offered insights into its regulatory mechanisms across different EVT cell models. Shared regulatory targets among these cell types and activation of trophoblast cell marker genes emphasized the importance of GATA3 in EVT differentiation and maturation. Knockdown of GATA3 in JEG3 cells led to repression of GATA3-induced epithelial-mesenchymal transition (EMT), as evidenced by changes in marker gene expression levels and enhanced migration ability. Additionally, interference with GATA3 accelerated cellular senescence, as indicated by reduced proliferation rates and increased activity levels for senescence-associated β-galactosidase enzyme, along with elevated expression levels for senescence-associated genes. This study provides comprehensive insights into the dual role of GATA3 in regulating EMT and cellular senescence during EVT differentiation, shedding light on the dynamic changes in GATA3 expression in normal and pathological placental conditions.

The transcriptome landscape of developing barley seeds

Cereal grains are an important source of food and feed. To provide comprehensive spatiotemporal information about biological processes in developing seeds of cultivated barley (Hordeum vulgare L. subsp. vulgare), we performed a transcriptomic study of the embryo, endosperm, and seed maternal tissues collected from grains 4-32 days after pollination. Weighted gene co-expression network and motif enrichment analyses identified specific groups of genes and transcription factors (TFs) potentially regulating barley seed tissue development. We defined a set of tissue-specific marker genes and families of TFs for functional studies of the pathways controlling barley grain development. Assessing selected groups of chromatin regulators revealed that epigenetic processes are highly dynamic and likely play a major role during barley endosperm development. The repressive H3K27me3 modification is globally reduced in endosperm tissues and at specific genes related to development and storage compounds. Altogether, this atlas uncovers the complexity of developmentally regulated gene expression in developing barley grains.

O-GlcNAcylation controls pro-fibrotic transcriptional regulatory signaling in myofibroblasts

Tissue injury causes activation of mesenchymal lineage cells into wound-repairing myofibroblasts (MFs), whose uncontrolled activity ultimately leads to fibrosis. Although this process is triggered by deep metabolic and transcriptional reprogramming, functional links between these two key events are not yet understood. Here, we report that the metabolic sensor post-translational modification O-linked β-D-N-acetylglucosaminylation (O-GlcNAcylation) is increased and required for myofibroblastic activation. Inhibition of protein O-GlcNAcylation impairs archetypal myofibloblast cellular activities including extracellular matrix gene expression and collagen secretion/deposition as defined in vitro and using ex vivo and in vivo murine liver injury models. Mechanistically, a multi-omics approach combining proteomic, epigenomic, and transcriptomic data mining revealed that O-GlcNAcylation controls the MF transcriptional program by targeting the transcription factors Basonuclin 2 (BNC2) and TEA domain transcription factor 4 (TEAD4) together with the Yes-associated protein 1 (YAP1) co-activator. Indeed, inhibition of protein O-GlcNAcylation impedes their stability leading to decreased functionality of the BNC2/TEAD4/YAP1 complex towards promoting activation of the MF transcriptional regulatory landscape. We found that this involves O-GlcNAcylation of BNC2 at Thr455 and Ser490 and of TEAD4 at Ser69 and Ser99. Altogether, this study unravels protein O-GlcNAcylation as a key determinant of myofibroblastic activation and identifies its inhibition as an avenue to intervene with fibrogenic processes.

Lineage regulators TFAP2C and NR5A2 function as bipotency activators in totipotent embryos

During the first lineage segregation, a mammalian totipotent embryo differentiates into the inner cell mass (ICM) and trophectoderm (TE). However, how transcription factors (TFs) regulate this earliest cell-fate decision in vivo remains elusive, with their regulomes primarily inferred from cultured cells. Here, we investigated the TF regulomes during the first lineage specification in early mouse embryos, spanning the pre-initiation, initiation, commitment, and maintenance phases. Unexpectedly, we found that TFAP2C, a trophoblast regulator, bound and activated both early TE and inner cell mass (ICM) genes at the totipotent (two- to eight-cell) stages ('bipotency activation'). Tfap2c deficiency caused downregulation of early ICM genes, including Nanog, Nr5a2, and Tdgf1, and early TE genes, including Tfeb and Itgb5, in eight-cell embryos. Transcription defects in both ICM and TE lineages were also found in blastocysts, accompanied by increased apoptosis and reduced cell numbers in ICMs. Upon trophoblast commitment, TFAP2C left early ICM genes but acquired binding to late TE genes in blastocysts, where it co-bound with CDX2, and later to extra-embryonic ectoderm (ExE) genes, where it cooperatively co-occupied with the former ICM regulator SOX2. Finally, 'bipotency activation' in totipotent embryos also applied to a pluripotency regulator NR5A2, which similarly bound and activated both ICM and TE lineage genes at the eight-cell stage. These data reveal a unique transcription circuity of totipotency underpinned by highly adaptable lineage regulators.

Accelerating Single-Cell Sequencing Data Analysis with SciDAP: A User-Friendly Approach

Single-cell (sc) RNA, ATAC and Multiome sequencing became powerful tools for uncovering biological and disease mechanisms. Unfortunately, manual analysis of sc data presents multiple challenges due to large data volumes and complexity of configuration parameters. This complexity, as well as not being able to reproduce a computational environment, affects the reproducibility of analysis results. The Scientific Data Analysis Platform (https://SciDAP.com) allows biologists without computational expertise to analyze sequencing-based data using portable and reproducible pipelines written in Common Workflow Language (CWL). Our suite of computational pipelines addresses the most common needs in scRNA-Seq, scATAC-Seq and scMultiome data analysis. When executed on SciDAP, it offers a user-friendly alternative to manual data processing, eliminating the need for coding expertise. In this protocol, we describe the use of SciDAP to analyze scMultiome data. Similar approaches can be used for analysis of scRNA-Seq, scATAC-Seq and scVDJ-Seq datasets.

Robust estimation of mean-variance relation

Accurate assessment of the mean-variance relation can benefit subsequent analysis in biomedical research. However, in most biomedical data, both the true mean and the true variance are unavailable. Instead, raw data are typically used to allow forming sample mean and sample variance in practice. In addition, different experimental conditions sometimes cause a slightly different mean-variance relation from the majority of the data in the same data set. To address these issues, we propose a semiparametric estimator, where we treat the uncertainty in the sample mean as a measurement error problem, the uncertainty in the sample variance as model error, and use a mixture model to account for different mean-variance relations. Asymptotic normality of the proposed method is established and its finite sample properties are demonstrated by simulation studies. The data application shows that the proposed method produces sensible results compared with methods either ignoring the uncertainty in the sample means or ignoring the potential different mean-variance relations.

Multifaceted SOX2-chromatin interaction underpins pluripotency progression in early embryos

Pioneer transcription factors (TFs), such as OCT4 and SOX2, play crucial roles in pluripotency regulation. However, the master TF-governed pluripotency regulatory circuitry was largely inferred from cultured cells. In this work, we investigated SOX2 binding from embryonic day 3.5 (E3.5) to E7.5 in the mouse. In E3.5 inner cell mass (ICM), SOX2 regulates the ICM-trophectoderm program but is dispensable for opening global enhancers. Instead, SOX2 occupies preaccessible enhancers in part opened by early-stage expressing TFs TFAP2C and NR5A2. SOX2 then widely redistributes when cells adopt naive and formative pluripotency by opening enhancers or poising them for rapid future activation. Hence, multifaceted pioneer TF-enhancer interaction underpins pluripotency progression in embryos, including a distinctive state in E3.5 ICM that bridges totipotency and pluripotency.

Maternal antibiotic exposure enhances ILC2 activation in neonates via downregulation of IFN1 signaling

Microbiota have an important function in shaping and priming neonatal immunity, although the cellular and molecular mechanisms underlying these effects remain obscure. Here we report that prenatal antibiotic exposure causes significant elevation of group 2 innate lymphoid cells (ILC2s) in neonatal lungs, in both cell numbers and functionality. Downregulation of type 1 interferon signaling in ILC2s due to diminished production of microbiota-derived butyrate represents the underlying mechanism. Mice lacking butyrate receptor GPR41 (Gpr41-/-) or type 1 interferon receptor IFNAR1 (Ifnar1-/-) recapitulate the phenotype of neonatal ILC2s upon maternal antibiotic exposure. Furthermore, prenatal antibiotic exposure induces epigenetic changes in ILC2s and has a long-lasting deteriorative effect on allergic airway inflammation in adult offspring. Prenatal supplementation of butyrate ameliorates airway inflammation in adult mice born to antibiotic-exposed dams. These observations demonstrate an essential role for the microbiota in the control of type 2 innate immunity at the neonatal stage, which suggests a therapeutic window for treating asthma in early life.

TET1 regulates responses to house dust mite by altering chromatin accessibility, DNA methylation, and gene expression in airway epithelial cells

Background

Previous studies have identified TET1 as a potential key regulator of genes linked to asthma. TET1 has been shown to transcriptionally respond to house dust mite extract, an allergen known to directly cause allergic asthma development, and regulate the expression of genes involved in asthma. How TET1 regulates expression of these genes, however, is unknown. TET1 is a DNA demethylase; therefore, most prior research on TET1-based gene regulation has focused on how TET1 affects methylation. However, TET1 can also interact directly with transcription factors and histone modifiers to regulate gene expression. Understanding how TET1 regulates expression to contribute to allergic responses and asthma development thus requires a comprehensive approach. To this end, we measured mRNA expression, DNA methylation, chromatin accessibility and histone modifications in control and TET1 knockdown human bronchial epithelial cells treated or untreated with house dust mite extract.

Results

Throughout our analyses, we detected strong similarities between the effects of TET1 knockdown alone and the effects of HDM treatment alone. One especially striking pattern was that both TET1 knockdown and HDM treatment generally led to decreased chromatin accessibility at largely the same genomic loci. Transcription factor enrichment analyses indicated that altered chromatin accessibility following the loss of TET1 may affect, or be affected by, CTCF and CEBP binding. TET1 loss also led to changes in DNA methylation, but these changes were generally in regions where accessibility was not changing.

Conclusions

TET1 regulates gene expression through different mechanisms (DNA methylation and chromatin accessibility) in different parts of the genome in the airway epithelial cells, which mediates inflammatory responses to allergen. Collectively, our data suggest novel molecular mechanisms through which TET1 regulates critical pathways following allergen challenges and contributes to the development of asthma.

NR5A2 connects zygotic genome activation to the first lineage segregation in totipotent embryos

Zygotic genome activation (ZGA) marks the beginning of the embryonic program for a totipotent embryo, which gives rise to the inner cell mass (ICM) where pluripotent epiblast arises, and extraembryonic trophectoderm. However, how ZGA is connected to the first lineage segregation in mammalian embryos remains elusive. Here, we investigated the role of nuclear receptor (NR) transcription factors (TFs), whose motifs are highly enriched and accessible from the 2-cell (2C) to 8-cell (8C) stages in mouse embryos. We found that NR5A2, an NR TF strongly induced upon ZGA, was required for this connection. Upon Nr5a2 knockdown or knockout, embryos developed beyond 2C normally with the zygotic genome largely activated. However, 4-8C-specific gene activation was substantially impaired and Nr5a2-deficient embryos subsequently arrested at the morula stage. Genome-wide chromatin binding analysis showed that NR5A2-bound cis-regulatory elements in both 2C and 8C embryos are strongly enriched for B1 elements where its binding motif is embedded. NR5A2 was not required for the global opening of its binding sites in 2C embryos but was essential to the opening of its 8C-specific binding sites. These 8C-specific, but not 2C-specific, binding sites are enriched near genes involved in blastocyst and stem cell regulation, and are often bound by master pluripotency TFs in blastocysts and embryonic stem cells (ESCs). Importantly, NR5A2 regulated key pluripotency genes Nanog and Pou5f1/Oct4, and primitive endoderm regulatory genes including Gata6 among many early ICM genes, as well as key trophectoderm regulatory genes including Tead4 and Gata3 at the 8C stage. By contrast, master pluripotency TFs NANOG, SOX2, and OCT4 targeted both early and late ICM genes in mouse ESCs. Taken together, these data identify NR5A2 as a key regulator in totipotent embryos that bridges ZGA to the first lineage segregation during mouse early development.

LHP1-mediated epigenetic buffering of subgenome diversity and defense responses confers genome plasticity and adaptability in allopolyploid wheat

Polyploidization is a major driver of genome diversification and environmental adaptation. However, the merger of different genomes may result in genomic conflicts, raising a major question regarding how genetic diversity is interpreted and regulated to enable environmental plasticity. By analyzing the genome-wide binding of 191 trans-factors in allopolyploid wheat, we identified like heterochromatin protein 1 (LHP1) as a master regulator of subgenome-diversified genes. Transcriptomic and epigenomic analyses of LHP1 mutants reveal its role in buffering the expression of subgenome-diversified defense genes by controlling H3K27me3 homeostasis. Stripe rust infection releases latent subgenomic variations by eliminating H3K27me3-related repression. The simultaneous inactivation of LHP1 homoeologs by CRISPR-Cas9 confers robust stripe rust resistance in wheat seedlings. The conditional repression of subgenome-diversified defenses ensures developmental plasticity to external changes, while also promoting neutral-to-non-neutral selection transitions and adaptive evolution. These findings establish an LHP1-mediated buffering system at the intersection of genotypes, environments, and phenotypes in polyploid wheat. Manipulating the epigenetic buffering capacity offers a tool to harness cryptic subgenomic variations for crop improvement.

Histone H3 lysine 27 trimethylation suppresses jasmonate biosynthesis and signaling to affect male fertility under high temperature in cotton

High-temperature (HT) stress causes male sterility in crops, thus decreasing yields. To explore the possible contribution of histone modifications to male fertility under HT conditions, we defined the histone methylation landscape for the marks histone H3 lysine 27 trimethylation (H3K27me3) and histone H3 lysine 4 trimethylation (H3K4me3) by chromatin immunoprecipitation sequencing (ChIP-seq) in two differing upland cotton (Gossypium hirsutum) varieties. We observed a global disruption in H3K4me3 and H3K27me3 modifications, especially H3K27me3, in cotton anthers subjected to HT. HT affected the bivalent H3K4me3-H3K27me3 modification more than either monovalent modification. We determined that removal of H3K27me3 at the promoters of jasmonate-related genes increased their expression, maintaining male fertility under HT in the HT-tolerant variety at the anther dehiscence stage. Modulating jasmonate homeostasis or signaling resulted in an anther indehiscence phenotype under HT. Chemical suppression of H3K27me3 deposition increased jasmonic acid contents and maintained male fertility under HT. In summary, our study provides new insights into the regulation of male fertility by histone modifications under HT and suggests a potential strategy for improving cotton HT tolerance.

Comprehensive dissection of meiotic DNA double-strand breaks and crossovers in cucumber

Meiotic recombination drives genetic diversity and crop genome optimization. In plant breeding, parents with favorable traits are crossed to create elite varieties. Different hybridizations produce diverse types of segment reshuffling between homologous chromosomes. However, little is known about the factors that cause hybrid-specific changes in crossovers (COs). Here, we constructed 2 F2 populations from crosses between a semiwild and 2 domesticated cucumber (Cucumis sativus) accessions and examined CO events. COs mainly occurred around genes and differed unevenly along chromosomes between the 2 hybrids. Fine-scale CO distributions were suppressed in regions of heterozygous structural variations (SVs) and were accelerated by high sequence polymorphism. C. sativus RADiation sensitive 51A (CsRAD51A) binding, histone H3 lysine 4 trimethylation (H3K4me3) modification, chromatin accessibility, and hypomethylation were positively associated with global CO landscapes and in local DNA double-strand break (DSB) hotspots and genes. The frequency and suppression of COs could be roughly predicted based on multiomic information. Differences in CO events between hybrids could be partially traced to distinct genetic and epigenetic features and were significantly associated with specific DSB hotspots and heterozygous SVs. Our findings identify the genomic and epigenetic features that contribute to CO formation and hybrid-specific divergence in cucumber and provide theoretical support for selecting parental combinations and manipulating recombination events at target genomic regions during plant breeding.

Analyzing histone ChIP-seq data with a bin-based probability of being signal

Histone ChIP-seq is one of the primary methods for charting the cellular epigenomic landscape, the components of which play a critical regulatory role in gene expression. Analyzing the activity of regulatory elements across datasets and cell types can be challenging due to shifting peak positions and normalization artifacts resulting from, for example, differing read depths, ChIP efficiencies, and target sizes. Moreover, broad regions of enrichment seen in repressive histone marks often evade detection by commonly used peak callers. Here, we present a simple and versatile method for identifying enriched regions in ChIP-seq data that relies on estimating a gamma distribution fit to non-overlapping 5kB genomic bins to establish a global background. We use this distribution to assign a probability of being signal (PBS) between zero and one to each 5 kB bin. This approach, while lower in resolution than typical peak-calling methods, provides a straightforward way to identify enriched regions and compare enrichments among multiple datasets, by transforming the data to values that are universally normalized and can be readily visualized and integrated with downstream analysis methods. We demonstrate applications of PBS for both broad and narrow histone marks, and provide several illustrations of biological insights which can be gleaned by integrating PBS scores with downstream data types.

Double DAP-seq uncovered synergistic DNA binding of interacting bZIP transcription factors

Many eukaryotic transcription factors (TF) form homodimer or heterodimer complexes to regulate gene expression. Dimerization of BASIC LEUCINE ZIPPER (bZIP) TFs are critical for their functions, but the molecular mechanism underlying the DNA binding and functional specificity of homo- versus heterodimers remains elusive. To address this gap, we present the double DNA Affinity Purification-sequencing (dDAP-seq) technique that maps heterodimer binding sites on endogenous genomic DNA. Using dDAP-seq we profile twenty pairs of C/S1 bZIP heterodimers and S1 homodimers in Arabidopsis and show that heterodimerization significantly expands the DNA binding preferences of these TFs. Analysis of dDAP-seq binding sites reveals the function of bZIP9 in abscisic acid response and the role of bZIP53 heterodimer-specific binding in seed maturation. The C/S1 heterodimers show distinct preferences for the ACGT elements recognized by plant bZIPs and motifs resembling the yeast GCN4 cis-elements. This study demonstrates the potential of dDAP-seq in deciphering the DNA binding specificities of interacting TFs that are key for combinatorial gene regulation.

cAMP/PKA signaling promotes AKT deactivation by reducing CIP2A expression, thereby facilitating decidualization

cAMP signaling is widely known to be indispensable for decidualization, but the details are not fully understood. Here, we show that cAMP signaling promotes AKT deactivation in endometrial stromal cells, which favors their decidualization. The deactivation of AKT is found to be a consequence of the reduced expression of several inhibitors of PP2A, the major phosphatase of AKT, with CIP2A being the most prominent. CIP2A reduction is obligatory for decidualization, as persistent CIP2A expression impairs chromatin remodeling and the expression of several decidualization markers (IGFBP1, PRL, MAOA, and IL-15). Furthermore, analyses of the responsiveness of the CIP2A promoter to cAMP signaling suggest the ETS family to be a bridge between cAMP signaling and CIP2A reduction. Our results provide novel insights into the role of cAMP signaling in decidualization and might benefit the development of novel therapies for decidualization deficiency, AKT-driven tumors, and the reverse, insulin resistance.

Multi-level interaction between HIF and AHR transcriptional pathways in kidney carcinoma

Hypoxia-inducible factor (HIF) and aryl hydrocarbon receptor (AHR) are members of the bHLH-PAS family of transcription factors that underpin cellular responses to oxygen and to endogenous and exogenous ligands, respectively, and have central roles in the pathogenesis of renal cancer. Composed of heterodimers, they share a common HIF-1β/ARNT subunit and similar DNA-binding motifs, raising the possibility of crosstalk between the two transcriptional pathways. Here, we identify both general and locus-specific mechanisms of interaction between HIF and AHR that act both antagonistically and cooperatively. Specifically, we observe competition for the common HIF-1β/ARNT subunit, in cis synergy for chromatin binding, and overlap in their transcriptional targets. Recently, both HIF and AHR inhibitors have been developed for the treatment of solid tumours. However, inhibition of one pathway may promote the oncogenic effects of the other. Therefore, our work raises important questions as to whether combination therapy targeting both of these pro-tumourigenic pathways might show greater efficacy than targeting each system independently.

Reorganization of three-dimensional chromatin architecture in Medicago truncatula under phosphorus deficiency

Emerging evidence reveals that the three-dimensional (3D) chromatin architecture plays a key regulatory role in various biological processes of plants. However, information on the 3D chromatin architecture of the legume model plant Medicago truncatula and its potential roles in the regulation of response to mineral nutrient deficiency are very limited. Using high-resolution chromosome conformation capture sequencing, we identified the 3D genome structure of M. truncatula in terms of A/B compartments, topologically associated domains (TADs) and chromatin loops. The gene density, expressional level, and active histone modification were higher in A compartments than in B compartments. Moreover, we analysed the 3D chromatin architecture reorganization in response to phosphorus (P) deficiency. The intra-chromosomal cis-interaction proportion was increased by P deficiency, and a total of 748 A/B compartment switch regions were detected. In these regions, density changes in H3K4me3 and H3K27ac modifications were associated with expression of P deficiency-responsive genes involved in root system architecture and hormonal responses. Furthermore, these genes enhanced P uptake and mobilization by increasing root surface area and strengthening signal transduction under P deficiency. These findings advance our understanding of the potential roles of 3D chromatin architecture in responses of plants in general, and in particular in M. truncatula, to P deficiency.

A chromosome-scale genome assembly of Dasypyrum villosum provides insights into its application as a broad-spectrum disease resistance resource for wheat improvement

Dasypyrum villosum is one of the most valuable gene resources in wheat improvement, especially for disease resistance. The mining of favorable genes from D. villosum is frustrated by the lack of a whole genome sequence. In this study, we generated a doubled-haploid line, 91C43^DH, using microspore culture and obtained a 4.05-GB high-quality, chromosome-scale genome assembly for D. villosum. The assembly contains39 727 high-confidence genes, and 85.31% of the sequences are repetitive. Two reciprocal translocation events were detected, and 7VS-4VL is a unique translocation in D. villosum. The prolamin seed storage protein-coding genes were found to be duplicated; in particular, the genes encoding low-molecular-weight glutenin at the Glu-V3 locus were significantly expanded. RNA sequencing (RNA-seq) analysis indicated that, after Blumeria graminearum f.sp tritici (Bgt) inoculation, there were more upregulated genes involved in the pattern-triggered immunity and effector-triggered immunity defense pathways in D. villosum than in Triticum urartu. MNase hypersensitive sequencing (MH-seq) identified two Bgt-inducible MH sites (MHSs), one in the promoter and one in the 3' terminal region of the powdery mildew resistance (Pm) gene NLR1-V. Each site had two subpeaks and they were termed MHS1 (MHS1.1/1.2) and MHS2 (MHS2.1/2.2). Bgt-inducible MHS2.2 was uniquely present in D. villosum, and MHS1.1 was more inducible in D. villosum than in wheat, suggesting that MHSs may be critical for regulation of NLR1-V expression and plant defense. In summary, this study provides a valuable genome resource for functional genomics studies and wheat-D. villosum introgression breeding. The identified regulatory mechanisms may also be exploited to develop new strategies for enhancing Pm resistance by optimizing gene expression in wheat.

The homeodomain of Oct4 is a dimeric binder of methylated CpG elements

Oct4 is essential to maintain pluripotency and has a pivotal role in establishing the germline. Its DNA-binding POU domain was recently found to bind motifs with methylated CpG elements normally associated with epigenetic silencing. However, the mode of binding and the consequences of this capability has remained unclear. Here, we show that Oct4 binds to a compact palindromic DNA element with a methylated CpG core (CpGpal) in alternative states of pluripotency and during cellular reprogramming towards induced pluripotent stem cells (iPSCs). During cellular reprogramming, typical Oct4 bound enhancers are uniformly demethylated, with the prominent exception of the CpGpal sites where DNA methylation is often maintained. We demonstrate that Oct4 cooperatively binds the CpGpal element as a homodimer, which contrasts with the ectoderm-expressed POU factor Brn2. Indeed, binding to CpGpal is Oct4-specific as other POU factors expressed in somatic cells avoid this element. Binding assays combined with structural analyses and molecular dynamic simulations show that dimeric Oct4-binding to CpGpal is driven by the POU-homeodomain whilst the POU-specific domain is detached from DNA. Collectively, we report that Oct4 exerts parts of its regulatory function in the context of methylated DNA through a DNA recognition mechanism that solely relies on its homeodomain.

Statistical Analysis in ChIP-seq-Related Applications

Chromatin immunoprecipitation sequencing (ChIP-seq) has been widely performed to identify protein binding information along the genome. The sequencing protocol is quite flexible and mature to measure different types of protein binding as long as sequencing parameters are properly tailored to accommodate protein features. Two distinct types of protein binding are point-source-like binding by transcription factors and diffused-distribution binding by histone modifications. Consequently, statistical approaches have been proposed to address ChIP-seq-related questions according to different protein features. In this chapter, we briefly summarize statistical principles, approaches, and tools that are widely implemented in modeling ChIP-seq data, from raw data quality control to final result reporting. We discuss the key solutions in addressing eight routine questions in ChIP-seq applications. We also include discussion on approaches fitting unique data features in different ChIP-seq types. We hope this chapter will serve as a brief guide, especially for ChIP-seq beginners, to provide them with a high-level overview to understand and design processing plans for their ChIP-seq experiments.

Dissecting Locus-Specific Chromatin Interactions by CRISPR CAPTURE

The spatiotemporal control of tissue-specific gene expression is coordinated by cis-regulatory elements (CREs) and associated trans-acting factors. Despite major advances in genome-wide annotation of candidate CREs, the in situ regulatory composition of the vast majority of CREs remain unknown. To address this challenge, we developed the CRISPR affinity purification in situ of regulatory elements (CAPTURE) toolbox that employs an in vivo biotinylated nuclease-deficient Cas9 (dCas9) protein and programmable single-guide RNAs (sgRNAs) to identify CRE-associated macromolecular complexes and chromatin looping. In this chapter, we provide a detailed protocol for implementing the latest iteration of the CRISPR-based CAPTURE methods to interrogate the molecular composition of locus-specific chromatin complexes and configuration in a mammalian genome.

β-hydroxybutyrate inhibits ferroptosis-mediated pancreatic damage in acute liver failure through the increase of H3K9bhb

Acute pancreatitis and hyperamylasemia are often seen in patients with acute liver failure (ALF). However, the underlying mechanisms remain elusive. This study describes pancreatic tissue damage and exocrine dysfunction in a mouse model of major-liver-resection-induced ALF. The analysis of 1,264 clinical cases of liver failure (LF) showed that the incidence of hyperamylasemia and hyperlipasemia in patients with LF is 5.5% and 20%, respectively. Metabolomic studies indicate that glutathione (GSH)-deficiency-caused ferroptosis contributes to pancreatic damage in mouse ALF. β-hydroxybutyrate (β-HB) is the only metabolite downregulated in the liver, serum, and pancreas. Our data suggest that β-HB protects pancreatic cells and tissues from GSH-deficiency-caused ferroptosis. β-HB administration in ALF mice restores the expression of ferroptosis-suppressor genes through histone H3 lysine 9 β-hydroxybutyrylation (H3K9bhb)-mediated chromatin opening. Our findings highlight β-HB as an endogenous metabolite regulating ferroptosis in the pancreas and extend our understanding of the pathophysiology of ALF-induced pancreatitis.

FoxA1 and FoxA2 control growth and cellular identity in NKX2-1-positive lung adenocarcinoma

Changes in cellular identity (also known as histologic transformation or lineage plasticity) can drive malignant progression and resistance to therapy in many cancers, including lung adenocarcinoma (LUAD). The lineage-specifying transcription factors FoxA1 and FoxA2 (FoxA1/2) control identity in NKX2-1/TTF1-negative LUAD. However, their role in NKX2-1-positive LUAD has not been systematically investigated. We find that Foxa1/2 knockout severely impairs tumorigenesis in KRAS-driven genetically engineered mouse models and human cell lines. Loss of FoxA1/2 leads to the collapse of a dual-identity state, marked by co-expression of pulmonary and gastrointestinal transcriptional programs, which has been implicated in LUAD progression. Mechanistically, FoxA1/2 loss leads to aberrant NKX2-1 activity and genomic localization, which in turn actively inhibits tumorigenesis and drives alternative cellular identity programs that are associated with non-proliferative states. This work demonstrates that FoxA1/2 expression is a lineage-specific vulnerability in NKX2-1-positive LUAD and identifies mechanisms of response and resistance to targeting FoxA1/2 in this disease.

DiffChIPL: a differential peak analysis method for high-throughput sequencing data with biological replicates based on limma

MOTIVATION

ChIP-seq detects protein-DNA interactions within chromatin, such as that of chromatin structural components and transcription machinery. ChIP-seq profiles are often noisy and variable across replicates, posing a challenge to the development of effective algorithms to accurately detect differential peaks. Methods have recently been designed for this purpose but sometimes yield conflicting results that are inconsistent with the underlying biology. Most existing algorithms perform well on limited datasets. To improve differential analysis of ChIP-seq, we present a novel Differential analysis method for ChIP-seq based on Limma (DiffChIPL).

RESULTS

DiffChIPL is adaptive to asymmetrical or symmetrical data and can accurately report global differences. We used simulated and real datasets for transcription factors (TFs) and histone modification marks to validate and benchmark our algorithm. DiffChIPL shows superior performance in sensitivity and false positive rate in different simulations and control datasets. DiffChIPL also performs well on real ChIP-seq, CUT&RUN, CUT&Tag and ATAC-seq datasets. DiffChIPL is an accurate and robust method, exhibiting better performance in differential analysis for a variety of applications including TF binding, histone modifications and chromatin accessibility.

AVAILABILITY AND IMPLEMENTATION

https://github.com/yancychy/DiffChIPL.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

HyperChIP: identification of hypervariable signals across ChIP-seq or ATAC-seq samples

Identifying genomic regions with hypervariable ChIP-seq or ATAC-seq signals across given samples is essential for large-scale epigenetic studies. In particular, the hypervariable regions across tumors from different patients indicate their heterogeneity and can contribute to revealing potential cancer subtypes and the associated epigenetic markers. We present HyperChIP as the first complete statistical tool for the task. HyperChIP uses scaled variances that account for the mean-variance dependence to rank genomic regions, and it increases the statistical power by diminishing the influence of true hypervariable regions on model fitting. A pan-cancer case study illustrates the practical utility of HyperChIP.

A balanced Oct4 interactome is crucial for maintaining pluripotency

Oct4 collaborates primarily with other transcriptional factors or coregulators to maintain pluripotency. However, how Oct4 exerts its function is still unclear. Here, we show that the Oct4 linker interface mediates competing yet balanced Oct4 protein interactions that are crucial for maintaining pluripotency. Oct4 linker mutant embryonic stem cells (ESCs) show decreased expression of self-renewal genes and increased expression of differentiation genes, resulting in impaired ESC self-renewal and early embryonic development. The linker mutation interrupts the balanced Oct4 interactome. In mutant ESCs, the interaction between Oct4 and Klf5 is decreased. In contrast, interactions between Oct4 and Cbx1, Ctr9, and Cdc73 are increased, disrupting the epigenetic state of ESCs. Control of the expression level of Klf5, Cbx1, or Cdc73 rebalances the Oct4 interactome and rescues the pluripotency of linker mutant ESCs, indicating that such factors interact with Oct4 competitively. Thus, we provide previously unidentified molecular insights into how Oct4 maintains pluripotency.

Enhancer-promoter interactions and transcription are largely maintained upon acute loss of CTCF, cohesin, WAPL or YY1

It remains unclear why acute depletion of CTCF (CCCTC-binding factor) and cohesin only marginally affects expression of most genes despite substantially perturbing three-dimensional (3D) genome folding at the level of domains and structural loops. To address this conundrum, we used high-resolution Micro-C and nascent transcript profiling in mouse embryonic stem cells. We find that enhancer-promoter (E-P) interactions are largely insensitive to acute (3-h) depletion of CTCF, cohesin or WAPL. YY1 has been proposed as a structural regulator of E-P loops, but acute YY1 depletion also had minimal effects on E-P loops, transcription and 3D genome folding. Strikingly, live-cell, single-molecule imaging revealed that cohesin depletion reduced transcription factor (TF) binding to chromatin. Thus, although CTCF, cohesin, WAPL or YY1 is not required for the short-term maintenance of most E-P interactions and gene expression, our results suggest that cohesin may facilitate TFs to search for and bind their targets more efficiently.

Characterizing batch effects and binding site-specific variability in ChIP-seq data

Multiple sources of variability can bias ChIP-seq data toward inferring transcription factor (TF) binding profiles. As ChIP-seq datasets increase in public repositories, it is now possible and necessary to account for complex sources of variability in ChIP-seq data analysis. We find that two types of variability, the batch effects by sequencing laboratories and differences between biological replicates, not associated with changes in condition or state, vary across genomic sites. This implies that observed differences between samples from different conditions or states, such as cell-type, must be assessed statistically, with an understanding of the distribution of obscuring noise. We present a statistical approach that characterizes both differences of interests and these source of variability through the parameters of a mixed effects model. We demonstrate the utility of our approach on a CTCF binding dataset composed of 211 samples representing 90 different cell-types measured across three different laboratories. The results revealed that sites exhibiting large variability were associated with sequence characteristics such as GC-content and low complexity. Finally, we identified TFs associated with high-variance CTCF sites using TF motifs documented in public databases, pointing the possibility of these being false positives if the sources of variability are not properly accounted for.

Differential ATAC-seq and ChIP-seq peak detection using ROTS

Changes in cellular chromatin states fine-tune transcriptional output and ultimately lead to phenotypic changes. Here we propose a novel application of our reproducibility-optimized test statistics (ROTS) to detect differential chromatin states (ATAC-seq) or differential chromatin modification states (ChIP-seq) between conditions. We compare the performance of ROTS to existing and widely used methods for ATAC-seq and ChIP-seq data using both synthetic and real datasets. Our results show that ROTS outperformed other commonly used methods when analyzing ATAC-seq data. ROTS also displayed the most accurate detection of small differences when modeling with synthetic data. We observed that two-step methods that require the use of a separate peak caller often more accurately called enrichment borders, whereas one-step methods without a separate peak calling step were more versatile in calling sub-peaks. The top ranked differential regions detected by the methods had marked correlation with transcriptional differences of the closest genes. Overall, our study provides evidence that ROTS is a useful addition to the available differential peak detection methods to study chromatin and performs especially well when applied to study differential chromatin states in ATAC-seq data.

A systematic dissection of the epigenomic heterogeneity of lung adenocarcinoma reveals two different subclasses with distinct prognosis and core regulatory networks

BACKGROUND

Lung adenocarcinoma (LUAD) is a highly malignant and heterogeneous tumor that involves various oncogenic genetic alterations. Epigenetic processes play important roles in lung cancer development. However, the variation in enhancer and super-enhancer landscapes of LUAD patients remains largely unknown. To provide an in-depth understanding of the epigenomic heterogeneity of LUAD, we investigate the H3K27ac histone modification profiles of tumors and adjacent normal lung tissues from 42 LUAD patients and explore the role of epigenetic alterations in LUAD progression.

RESULTS

A high intertumoral epigenetic heterogeneity is observed across the LUAD H3K27ac profiles. We quantitatively model the intertumoral variability of H3K27ac levels at proximal gene promoters and distal enhancers and propose a new epigenetic classification of LUAD patients. Our classification defines two LUAD subgroups which are highly related to histological subtypes. Group II patients have significantly worse prognosis than group I, which is further confirmed in the public TCGA-LUAD cohort. Differential RNA-seq analysis between group I and group II groups reveals that those genes upregulated in group II group tend to promote cell proliferation and induce cell de-differentiation. We construct the gene co-expression networks and identify group-specific core regulators. Most of these core regulators are linked with group-specific regulatory elements, such as super-enhancers. We further show that CLU is regulated by 3 group I-specific core regulators and works as a novel tumor suppressor in LUAD.

CONCLUSIONS

Our study systematically characterizes the epigenetic alterations during LUAD progression and provides a new classification model that is helpful for predicting patient prognosis.