Search

Copyright of Molecular Biology Reports is the property of Springer Nature and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

Background: The relationship between epigenomics and plant specialised metabolism remains largely unexplored despite the fundamental importance of epigenomics in gene regulation and, potentially, yield of products of plant specialised metabolic pathways. The glandular trichomes of Cannabis sativa are an emerging model system that produce large quantities of cannabinoid and terpenoid specialised metabolites with known medicinal and commercial value. To address this lack of epigenomic data, we mapped H3K4 trimethylation, H3K56 acetylation, H3K27 trimethylation post-translational modifications and the histone variant H2A.Z, using chromatin immunoprecipitation, in C. sativa glandular trichomes, leaf, and stem tissues. Corresponding transcriptomic (RNA-seq) datasets were integrated, and tissue-specific analyses conducted to relate chromatin states to glandular trichome specific gene expression. Results: The promoters of cannabinoid and terpenoid biosynthetic genes, specialised metabolite transporter genes, defence related genes, and starch and sucrose metabolism were enriched specifically in trichomes for histone marks H3K4me3 and H3K56ac, consistent with active transcription. We identified putative trichome-specific enhancer elements by identifying intergenic regions of H3K56ac enrichment, a histone mark that maintains enhancer accessibility, then associated these to putative target genes using the tissue specific gene transcriptomic data. Bi-valent chromatin loci specific to glandular trichomes, marked with H3K4 trimethylation and H3K27 trimethylation, were associated with genes of MAPK signalling pathways and plant specialised metabolism pathways, supporting recent hypotheses that implicate bi-valent chromatin in plant defence. The histone variant H2A.Z was largely found in intergenic regions and enriched in chromatin that contained genes involved in DNA homeostasis. Conclusion: We report the first genome-wide histone post-translational modification maps for C. sativa glandular trichomes, and more broadly for glandular trichomes in plants. Our findings have implications in plant adaptation and stress responses and provide a basis for enhancer-mediated, targeted, gene transformation studies in plant glandular trichomes. [ABSTRACT FROM AUTHOR]

The present study aimed to investigate age-related changes in histone variant H3.3 and its role in the aging process of mouse tibialis anterior muscle. H3.3 level significantly increased with age and correlated with H3K27me3 level. Acute exercise successfully upregulated the target gene expression in 8-wk-old mice, whereas no upregulation was noted in 53-wk-old mice. H3K27me3 level was increased at these loci in response to acute exercise in 8-wk-old mice. However, in 53-wk-old mice, H3.3 and H3K27me3 levels were increased at rest and were not affected by acute exercise. Furthermore, forced H3.3 expression in the skeletal muscle of 8-wk-old mice led to a gradual improvement in motor function. The results suggest that age-related H3.3 accumulation induces the formation of repressive chromatin in the mouse tibialis anterior muscle. However, H3.3 accumulation also appears to play a positive role in enhancing skeletal muscle function. [ABSTRACT FROM AUTHOR]

Abstract Background Response to oxidative stress is universal in almost all organisms and the mitochondrial membrane protein, BbOhmm, negatively affects oxidative stress responses and virulence in the insect fungal pathogen, Beauveria bassiana. Nothing further, however, is known concerning how BbOhmm and this phenomenon is regulated. Results Three oxidative stress response regulating Zn2Cys6 transcription factors (BbOsrR1, 2, and 3) were identified and verified via chromatin immunoprecipitation (ChIP)-qPCR analysis as binding to the BbOhmm promoter region, with BbOsrR2 showing the strongest binding. Targeted gene knockout of BbOsrR1 or BbOsrR3 led to decreased BbOhmm expression and consequently increased tolerances to free radical generating compounds (H2O2 and menadione), whereas the ΔBbOsrR2 strain showed increased BbOhmm expression with concomitant decreased tolerances to these compounds. RNA and ChIP sequencing analysis revealed that BbOsrR1 directly regulated a wide range of antioxidation and transcription-associated genes, negatively affecting the expression of the BbClp1 cyclin and BbOsrR2. BbClp1 was shown to localize to the cell nucleus and negatively mediate oxidative stress responses. BbOsrR2 and BbOsrR3 were shown to feed into the Fus3-MAPK pathway in addition to regulating antioxidation and detoxification genes. Binding motifs for the three transcription factors were found to partially overlap in the promoter region of BbOhmm and other target genes. Whereas BbOsrR1 appeared to function independently, co-immunoprecipitation revealed complex formation between BbClp1, BbOsrR2, and BbOsrR3, with BbClp1 partially regulating BbOsrR2 phosphorylation. Conclusions These findings reveal a regulatory network mediated by BbOsrR1 and the formation of a BbClp1-BbOsrR2-BbOsrR3 complex that orchestrates fungal oxidative stress responses.

Background: Response to oxidative stress is universal in almost all organisms and the mitochondrial membrane protein, BbOhmm, negatively affects oxidative stress responses and virulence in the insect fungal pathogen, Beauveria bassiana. Nothing further, however, is known concerning how BbOhmm and this phenomenon is regulated. Results: Three oxidative stress response regulating Zn2Cys6 transcription factors (BbOsrR1, 2, and 3) were identified and verified via chromatin immunoprecipitation (ChIP)-qPCR analysis as binding to the BbOhmm promoter region, with BbOsrR2 showing the strongest binding. Targeted gene knockout of BbOsrR1 or BbOsrR3 led to decreased BbOhmm expression and consequently increased tolerances to free radical generating compounds (H2O2 and menadione), whereas the ΔBbOsrR2 strain showed increased BbOhmm expression with concomitant decreased tolerances to these compounds. RNA and ChIP sequencing analysis revealed that BbOsrR1 directly regulated a wide range of antioxidation and transcription-associated genes, negatively affecting the expression of the BbClp1 cyclin and BbOsrR2. BbClp1 was shown to localize to the cell nucleus and negatively mediate oxidative stress responses. BbOsrR2 and BbOsrR3 were shown to feed into the Fus3-MAPK pathway in addition to regulating antioxidation and detoxification genes. Binding motifs for the three transcription factors were found to partially overlap in the promoter region of BbOhmm and other target genes. Whereas BbOsrR1 appeared to function independently, co-immunoprecipitation revealed complex formation between BbClp1, BbOsrR2, and BbOsrR3, with BbClp1 partially regulating BbOsrR2 phosphorylation. Conclusions: These findings reveal a regulatory network mediated by BbOsrR1 and the formation of a BbClp1-BbOsrR2-BbOsrR3 complex that orchestrates fungal oxidative stress responses. [ABSTRACT FROM AUTHOR]

Activating transcription factor 5 (ATF5) is a transcription factor that belongs to the cAMP-response element-binding protein/ATF family and is essential for the differentiation and survival of sensory neurons in mouse olfactory organs. However, transcriptional target genes for ATF5 have yet to be identified. In the present study, chromatin immunoprecipitation-quantitative polymerase chain reaction (ChIP-qPCR) experiments were performed to verify ATF5 target genes in the main olfactory epithelium and vomeronasal organ in the postnatal pups. ChIP-qPCR was conducted using hemagglutinin (HA)-tagged ATF5 knock-in olfactory organs. The results obtained demonstrated that ATF5-HA fusion proteins bound to the CCAAT/enhancer-binding protein-ATF response element (CARE) site in the enhancer region of nescient helix-loop-helix 1 (Nhlh1), a transcription factor expressed in differentiating olfactory and vomeronasal sensory neurons. Nhlh1 mRNA expression was downregulated in ATF5-deficient (ATF5−/−) olfactory organs. The LIM/homeobox protein transcription factor Lhx2 co-localized with ATF5 in the nuclei of olfactory and vomeronasal sensory neurons and bound to the homeodomain site proximal to the CARE site in the Nhlh1 gene. The CARE region of the Nhlh1 gene was enriched by the active enhancer marker, acetyl-histone H3 (Lys27). The present study identified Nhlh1 as a novel target gene for ATF5 in murine olfactory organs. ATF5 may upregulate Nhlh1 expression in concert with Lhx2, thereby promoting the differentiation of olfactory and vomeronasal sensory neurons. [ABSTRACT FROM AUTHOR]

Summary: Transcription factor (TF) gene knockout or knockdown experiments provide comprehensive downstream effects on gene regulation. However, distinguishing primary direct effects from secondary effects remains challenging. To assess the direct effect of TF binding events, we present a protocol for establishing a doxycycline (Dox)-inducible CRISPRd system in human pluripotent stem cells (hPSCs). We describe the steps for establishing CRISPRd host hPSCs, designing and preparing single-guide RNA (sgRNA) expression lentivirus vectors, generating CRISPRd hPSCs transduced with sgRNAs, and analyzing CRISPRd TF-block effects by chromatin immunoprecipitation (ChIP)-qPCR.For complete details on the use and execution of this protocol, please refer to Matsui et al.1 : Publisher’s note: Undertaking any experimental protocol requires adherence to local institutional guidelines for laboratory safety and ethics.

Background & aimsThe stimulator of interferon genes (STING)/TANK-binding kinase 1 (TBK1) pathway is vital in mediating innate immune and inflammatory responses during oxidative/endoplasmic reticulum (ER) stress. However, it remains unknown whether macrophage thioredoxin-interacting protein (TXNIP) may regulate TBK1 function and cell death pathways during oxidative/ER stress.MethodsA mouse model of hepatic ischaemia/reperfusion injury (IRI), the primary hepatocytes, and bone marrow-derived macrophages were used in the myeloid-specific TXNIP knockout (TXNIPM-KO) and TXNIP-proficient (TXNIPFL/FL) mice.ResultsThe TXNIPM-KO mice were resistant to ischaemia/reperfusion (IR) stress-induced liver damage with reduced serum alanine aminotransferase (ALT)/aspartate aminotransferase (AST) levels, macrophage/neutrophil infiltration, and pro-inflammatory mediators compared with the TXNIPFL/FL controls. IR stress increased TXNIP, p-STING, and p-TBK1 expression in ischaemic livers. However, TXNIPM-KO inhibited STING, TBK1, interferon regulatory factor 3 (IRF3), and NF-κB activation with interferon-β (IFN-β) expression. Interestingly, TXNIPM-KO augmented nuclear factor (erythroid-derived 2)-like 2 (NRF2) activity, increased antioxidant gene expression, and reduced macrophage reactive oxygen species (ROS) production and hepatic apoptosis/necroptosis in IR-stressed livers. Mechanistically, macrophage TXNIP deficiency promoted cylindromatosis (CYLD), which colocalised and interacted with NADPH oxidase 4 (NOX4) to enhance NRF2 activity by deubiquitinating NOX4. Disruption of macrophage NRF2 or its target gene 2',5' oligoadenylate synthetase-like 1 (OASL1) enhanced Ras GTPase-activating protein-binding protein 1 (G3BP1) and TBK1-mediated inflammatory response. Notably, macrophage OASL1 deficiency induced hepatocyte apoptotic peptidase activating factor 1 (APAF1), cytochrome c, and caspase-9 activation, leading to increased caspase-3-initiated apoptosis and receptor-interacting serine/threonine-protein kinase 3 (RIPK3)-mediated necroptosis.ConclusionsMacrophage TXNIP deficiency enhances CYLD activity and activates the NRF2-OASL1 signalling, controlling IR stress-induced liver injury. The target gene OASL1 regulated by NRF2 is crucial for modulating STING-mediated TBK1 activation and Apaf1/cytochrome c/caspase-9-triggered apoptotic/necroptotic cell death pathway. Our findings underscore a novel role of macrophage TXNIP-mediated CYLD-NRF2-OASL1 axis in stress-induced liver inflammation and cell death, implying the potential therapeutic targets in liver inflammatory diseases.Lay summaryLiver inflammation and injury induced by ischaemia and reperfusion (the absence of blood flow to the liver tissue followed by the resupply of blood) is a significant cause of hepatic dysfunction and failure following liver transplantation, resection, and haemorrhagic shock. Herein, we uncover an underlying mechanism that contributes to liver inflammation and cell death in this setting and could be a therapeutic target in stress-induced liver inflammatory injury.

Abstract Background Bivalent chromatin is an exemplar of epigenetic plasticity. This co-occurrence of active-associated H3K4me3 and inactive-associated H3K27me3 histone modifications on opposite tails of the same nucleosome occurs predominantly at promoters that are poised for future transcriptional upregulation or terminal silencing. We know little of the dynamics, resolution, and regulation of this chromatin state outside of embryonic stem cells where it was first described. This is partly due to the technical challenges distinguishing bone-fide bivalent chromatin, where both marks are on the same nucleosome, from allelic or sample heterogeneity where there is a mix of H3K4me3-only and H3K27me3-only mononucleosomes. Results Here, we present a robust and sensitive method to accurately map bivalent chromatin genome-wide, along with controls, from as little as 2 million cells. We optimized and refined the sequential ChIP protocol which uses two sequential overnight immunoprecipitation reactions to robustly purify nucleosomes that are truly bivalent and contain both H3K4me3 and H3K27me3 modifications. Our method generates high quality genome-wide maps with strong peak enrichment and low background, which can be analyzed using standard bioinformatic packages. Using this method, we detect 8,789 bivalent regions in mouse embryonic stem cells corresponding to 3,918 predominantly CpG rich and developmentally regulated gene promoters. Furthermore, profiling Dppa2/4 knockout mouse embryonic stem cells, which lose both H3K4me3 and H3K27me3 at approximately 10% of bivalent promoters, demonstrated the ability of our method to capture bivalent chromatin dynamics. Conclusions Our optimized sequential reChIP method enables high-resolution genome-wide assessment of bivalent chromatin together with all required controls in as little as 2 million cells. We share a detailed protocol and guidelines that will enable bivalent chromatin landscapes to be generated in a range of cellular contexts, greatly enhancing our understanding of bivalent chromatin and epigenetic plasticity beyond embryonic stem cells.

Polycomb-group proteins play critical roles in gene silencing through the deposition of histone H3 lysine 27 trimethylation (H3K27me3) and chromatin compaction. This process is essential for embryonic stem cell (ESC) pluripotency, differentiation, and development. Polycomb repressive complex 2 (PRC2) can both read and write H3K27me3, enabling progressive spreading of H3K27me3 on the linear genome. Long-range Polycomb-associated DNA contacts have also been described, but their regulation and role in gene silencing remain unclear. Here, we apply H3K27me3 HiChIP, a protein-directed chromosome conformation method, and optical reconstruction of chromatin architecture to profile long-range Polycomb-associated DNA loops that span tens to hundreds of megabases across multiple topological associated domains in mouse ESCs and human induced pluripotent stem cells. We find that H3K27me3 loop anchors are enriched for Polycomb nucleation points and coincide with key developmental genes. Genetic deletion of H3K27me3 loop anchors results in disruption of spatial contact between distant loci and altered H3K27me3 in cis, both locally and megabases away on the same chromosome. In mouse embryos, loop anchor deletion leads to ectopic activation of the partner gene, suggesting that Polycomb-associated loops control gene silencing during development. Further, we find that alterations in PRC2 occupancy resulting from an RNA binding–deficient EZH2 mutant are accompanied by loss of Polycomb-associated DNA looping. Together, these results suggest PRC2 uses RNA binding to enhance long-range chromosome folding and H3K27me3 spreading. Developmental gene loci have unique roles in Polycomb spreading, emerging as important architectural elements of the epigenome.

CD8+ T cells are critical for the immune response to pathogens and tumors, and CD8+ T cell memory protects against repeat infections. In this study, we identify the activating transcription factor 7 interacting protein (ATF7ip) as a critical regulator of CD8+ T cell immune responses. Mice with a T cell-specific deletion of ATF7ip have a CD8+ T cell-intrinsic enhancement of Il7r expression and Il2 expression leading to enhanced effector and memory responses. Chromatin immunoprecipitation sequencing studies identified ATF7ip as a repressor of Il7r and Il2 gene expression through the deposition of the repressive histone mark H3K9me3 at the Il7r gene and Il2-Il21 intergenic region. Interestingly, ATF7ip targeted transposable elements for H3K9me3 deposition at both the IL7r locus and the Il2-Il21 intergenic region, indicating that ATF7ip silencing of transposable elements is important for regulating CD8+ T cell function. These results demonstrate a new epigenetic pathway by which IL-7R and IL-2 production are constrained in CD8+ T cells, and this may open up new avenues for modulating their production.

The transcriptional regulatory network in prokaryotes controls global gene expression mostly through transcription factors (TFs), which are DNA-binding proteins. Chromatin immunoprecipitation (ChIP) with DNA sequencing methods can identify TF binding sites across the genome, providing a bottom-up, mechanistic understanding of how gene expression is regulated. ChIP provides indispensable evidence toward the goal of acquiring a comprehensive understanding of cellular adaptation and regulation, including condition-specificity. ChIP-derived data's importance and labor-intensiveness motivate its broad dissemination and reuse, which is currently an unmet need in the prokaryotic domain. To fill this gap, we present proChIPdb (prochipdb.org), an information-rich, interactive web database. This website collects public ChIP-seq/-exo data across several prokaryotes and presents them in dashboards that include curated binding sites, nucleotide-resolution genome viewers, and summary plots such as motif enrichment sequence logos. Users can search for TFs of interest or their target genes, download all data, dashboards, and visuals, and follow external links to understand regulons through biological databases and the literature. This initial release of proChIPdb covers diverse organisms, including most major TFs of Escherichia coli, and can be expanded to support regulon discovery across the prokaryotic domain.

Regulation of gene expression requires the combinatorial binding of sequence-specific transcription factors (TFs) at promoters and enhancers. Prior studies showed that alterations in the spacing between TF binding sites can influence promoter and enhancer activity. However, the relative importance of TF spacing alterations resulting from naturally occurring insertions and deletions (InDels) has not been systematically analyzed. To address this question, we first characterized the genome-wide spacing relationships of 73 TFs in human K562 cells as determined by ChIP-seq (chromatin immunoprecipitation sequencing). We found a dominant pattern of a relaxed range of spacing between collaborative factors, including 45 TFs exclusively exhibiting relaxed spacing with their binding partners. Next, we exploited millions of InDels provided by genetically diverse mouse strains and human individuals to investigate the effects of altered spacing on TF binding and local histone acetylation. These analyses suggested that spacing alterations resulting from naturally occurring InDels are generally tolerated in comparison to genetic variants directly affecting TF binding sites. To experimentally validate this prediction, we introduced synthetic spacing alterations between PU.1 and C/EBPβ binding sites at six endogenous genomic loci in a macrophage cell line. Remarkably, collaborative binding of PU.1 and C/EBPβ at these locations tolerated changes in spacing ranging from 5 bp increase to >30 bp decrease. Collectively, these findings have implications for understanding mechanisms underlying enhancer selection and for the interpretation of non-coding genetic variation.

Accurate investigations of neural circuitry require specific genetic access to individual circuit elements, i.e., the myriad neuronal cell-types in the brain. However, native promoters cannot achieve this because while most genes are expressed in the brain, few are expressed in a single neuronal cell-type. We recently used enhancers, the subcomponents of the transcriptional apparatus which tell promoters when and where to express, combined with heterologous minimal promoters to increase specificity of transgene expression, an approach we call Enhancer-Driven Gene Expression (EDGE). As we discuss, EDGE is a marked improvement in specificity over native promoters, but still requires careful anatomical analysis to avoid off-target effects. In this study we present a more complete set of genomic markers from the mouse brain and characterize a novel EDGE viral vector capable of specifically driving expression in distinct subtypes of hippocampal neurons, even though it can express in other cell-types elsewhere. The advent of cell-type specific viral tools in wild-type animals provides a powerful strategy for neural circuit investigation and holds promise for studies using animal models for which transgenic tools are not available. [ABSTRACT FROM AUTHOR]

Background: Bivalent chromatin is an exemplar of epigenetic plasticity. This co-occurrence of active-associated H3K4me3 and inactive-associated H3K27me3 histone modifications on opposite tails of the same nucleosome occurs predominantly at promoters that are poised for future transcriptional upregulation or terminal silencing. We know little of the dynamics, resolution, and regulation of this chromatin state outside of embryonic stem cells where it was first described. This is partly due to the technical challenges distinguishing bone-fide bivalent chromatin, where both marks are on the same nucleosome, from allelic or sample heterogeneity where there is a mix of H3K4me3-only and H3K27me3-only mononucleosomes. Results: Here, we present a robust and sensitive method to accurately map bivalent chromatin genome-wide, along with controls, from as little as 2 million cells. We optimized and refined the sequential ChIP protocol which uses two sequential overnight immunoprecipitation reactions to robustly purify nucleosomes that are truly bivalent and contain both H3K4me3 and H3K27me3 modifications. Our method generates high quality genome-wide maps with strong peak enrichment and low background, which can be analyzed using standard bioinformatic packages. Using this method, we detect 8,789 bivalent regions in mouse embryonic stem cells corresponding to 3,918 predominantly CpG rich and developmentally regulated gene promoters. Furthermore, profiling Dppa2/4 knockout mouse embryonic stem cells, which lose both H3K4me3 and H3K27me3 at approximately 10% of bivalent promoters, demonstrated the ability of our method to capture bivalent chromatin dynamics. Conclusions: Our optimized sequential reChIP method enables high-resolution genome-wide assessment of bivalent chromatin together with all required controls in as little as 2 million cells. We share a detailed protocol and guidelines that will enable bivalent chromatin landscapes to be generated in a range of cellular contexts, greatly enhancing our understanding of bivalent chromatin and epigenetic plasticity beyond embryonic stem cells. [ABSTRACT FROM AUTHOR]

Decidualization of endometrial stroma is a key step in embryo implantation and its abnormality often leads to pregnancy failure. Stromal decidualization is a very complex process that is co‐regulated by estrogen, progesterone and many local factors. The signaling protein SHP2 encoded by PTPN11 is dynamically expressed in decidualized endometrial stroma and mediates and integrates various signals to govern the decidualization. In the present study, we investigate the mechanism of PTPN11 gene transcription. Estrogen, progesterone and cAMP co‐induced decidualization of human endometrial stromal cell in vitro, but only progesterone and cAMP induced SHP2 expression. Using the luciferase reporter, we refined a region from −229 bp to +1 bp in the PTPN11 gene promoter comprising the transcriptional core regions that respond to progesterone and cAMP. Progesterone receptor (PGR) and cAMP‐responsive element‐binding protein 1 (CREB1) were predicted to be transcription factors in this core region by bioinformatic methods. The direct binding of PGR and CREB1 on the PTPN11 promoter was confirmed by electrophoretic mobility and chromatin immunoprecipitation in vitro. Knockdown of PGR and CREB1 protein significantly inhibited the expression of SHP2 induced by medroxyprogesterone acetate and cAMP. These results demonstrate that transcription factors PGR and CREB1 bind to the PTPN11 promoter to regulate the expression of SHP2 in response to decidual signals. Our results explain the transcriptional expression mechanism of SHP2 during decidualization and promote the understanding of the mechanism of decidualization of stromal cells. [ABSTRACT FROM AUTHOR]

Sequence-specific transcription factors (TFs) recognize motifs of related nucleotide sequences at their DNA binding sites. Upon binding at these sites, TFs regulate critical molecular processes such as gene expression. It is widely assumed that a TF recognizes a single "canonical" motif, although recent studies have identified additional "non-canonical" motifs for some TFs. A comprehensive approach to identify non-canonical DNA binding motifs and the functional importance of those motifs' matches in the human genome is necessary for fully understanding the mechanisms of TF-regulated molecular processes in human cells. To address this need, we developed a statistical pipeline for in vitro HT-SELEX data that identifies and characterizes the distributions of non-canonical TF motifs in a stringent manner. Analyzing ~170 human TFs' HT-SELEX data, we found non-canonical motifs for 19 TFs (11%). These non-canonical motifs occur independently of the TFs' canonical motifs. Non-canonical motif occurrences in the human genome show similar evolutionary conservation to canonical motif occurrences, explain TF binding in locations without canonical motifs, and occur within gene promoters and epigenetically marked regulatory sequences in human cell lines and tissues. Our approach and collection of non-canonical motifs expand current understanding of functionally relevant DNA binding sites for human TFs.

Molecular-level understanding of plasma cell (PC) differentiation has been modeled using lipopolysaccharide (LPS) stimulation in vitro. However, this system does not involve the B-cell receptor (BCR)-a critical component of B cell biology. Here, we present a protocol for in vitro PC differentiation system dependent on BCR signaling that easily scales up for cell number-demanding applications, including protein complex purification. We describe how to set up this system and detail applications for endogenous complex purification of chromatin-associated proteins. For further details on the use and execution of this protocol, please refer to Sciammas et al. (2011) and Ochiai et al. (2018, 2020).

Chromatin immunoprecipitation followed by massively parallel, high throughput sequencing (ChIP-seq) is the method of choice for genome-wide identification of DNA segments bound by specific transcription factors or in chromatin with particular histone modifications. However, the quality of ChIP-seq datasets varies widely, with a substantial fraction being of intermediate to poor quality. Thus, it is important to discern and control the factors that contribute to variation in ChIP-seq. In this study, we focused on sonication, a user-controlled variable, to produce sheared chromatin. We systematically varied the amount of shearing of fixed chromatin from a mouse erythroid cell line, carefully measuring the distribution of resultant fragment lengths prior to ChIP-seq. This systematic study was complemented with a retrospective analysis of additional experiments. We found that the level of sonication had a pronounced impact on the quality of ChIP-seq signals. Over-sonication consistently reduced quality, while the impact of under-sonication differed among transcription factors, with no impact on sites bound by CTCF but frequently leading to the loss of sites occupied by TAL1 or bound by POL2. The bound sites not observed in low-quality datasets were inferred to be a mix of both direct and indirect binding. We leveraged these findings to produce a set of CTCF ChIP-seq datasets in rare, primary hematopoietic progenitor cells. Our observation that the amount of chromatin sonication is a key variable in success of ChIP-seq experiments indicates that monitoring the level of sonication can improve ChIP-seq quality and reproducibility and facilitate ChIP-seq in rare cell types.

BackgroundCoronary artery disease (CAD) is a multifactorial condition with both genetic and exogenous causes. The contribution of tissue-specific functional networks to the development of atherosclerosis remains largely unclear. The aim of this study was to identify and characterize central regulators and networks leading to atherosclerosis.MethodsBased on several hundred genes known to affect atherosclerosis risk in mouse (as demonstrated in knockout models) and human (as shown by genome-wide association studies), liver gene regulatory networks were modeled. The hierarchical order and regulatory directions of genes within the network were based on Bayesian prediction models, as well as experimental studies including chromatin immunoprecipitation DNA-sequencing, chromatin immunoprecipitation mass spectrometry, overexpression, small interfering RNA knockdown in mouse and human liver cells, and knockout mouse experiments. Bioinformatics and correlation analyses were used to clarify associations between central genes and CAD phenotypes in both human and mouse.ResultsThe transcription factor MAFF (MAF basic leucine zipper transcription factor F) interacted as a key driver of a liver network with 3 human genes at CAD genome-wide association studies loci and 11 atherosclerotic murine genes. Most importantly, expression levels of the low-density lipoprotein receptor (LDLR) gene correlated with MAFF in 600 CAD patients undergoing bypass surgery (STARNET [Stockholm-Tartu Atherosclerosis Reverse Network Engineering Task]) and a hybrid mouse diversity panel involving 105 different inbred mouse strains. Molecular mechanisms of MAFF were tested in noninflammatory conditions and showed positive correlation between MAFF and LDLR in vitro and in vivo. Interestingly, after lipopolysaccharide stimulation (inflammatory conditions), an inverse correlation between MAFF and LDLR in vitro and in vivo was observed. Chromatin immunoprecipitation mass spectrometry revealed that the human CAD genome-wide association studies candidate BACH1 (BTB domain and CNC homolog 1) assists MAFF in the presence of lipopolysaccharide stimulation with respective heterodimers binding at the MAF recognition element of the LDLR promoter to transcriptionally downregulate LDLR expression.ConclusionsThe transcription factor MAFF was identified as a novel central regulator of an atherosclerosis/CAD-relevant liver network. MAFF triggered context-specific expression of LDLR and other genes known to affect CAD risk. Our results suggest that MAFF is a missing link between inflammation, lipid and lipoprotein metabolism, and a possible treatment target.

Accurate investigations of neural circuitry require specific genetic access to individual circuit elements, i.e., the myriad neuronal cell-types in the brain. However, native promoters cannot achieve this because while most genes are expressed in the brain, few are expressed in a single neuronal cell-type. We recently used enhancers, the subcomponents of the transcriptional apparatus which tell promoters when and where to express, combined with heterologous minimal promoters to increase specificity of transgene expression, an approach we call Enhancer-Driven Gene Expression (EDGE). As we discuss, EDGE is a marked improvement in specificity over native promoters, but still requires careful anatomical analysis to avoid off-target effects. In this study we present a more complete set of genomic markers from the mouse brain and characterize a novel EDGE viral vector capable of specifically driving expression in distinct subtypes of hippocampal neurons, even though it can express in other cell-types elsewhere. The advent of cell-type specific viral tools in wild-type animals provides a powerful strategy for neural circuit investigation and holds promise for studies using animal models for which transgenic tools are not available.

Significant advancements in understanding disease mechanisms can occur through combined analysis of next-generation sequencing datasets generated using purified cell populations. Here, we detail our optimized protocol for purification of mouse hepatic macrophages (or other liver non-parenchymal populations) suitable for use in various next-generation sequencing protocols. An alternative framework is described for sorting pre-fixed hepatic nuclei populations. This strategy has the advantage of rapidly preserving the nuclei and can facilitate success with ChIP-seq for more challenging molecules. For complete details on the use and execution of these protocols, please refer to Muse et al. (2018), Sakai et al. (2019), and Seidman et al. (2020).

BackgroundCo-localized combinations of histone modifications ("chromatin states") have been shown to correlate with promoter and enhancer activity. Changes in chromatin states over multiple time points ("chromatin state trajectories") have previously been analyzed at promoter and enhancers separately. With the advent of time series Hi-C data it is now possible to connect promoters and enhancers and to analyze chromatin state trajectories at promoter-enhancer pairs.ResultsWe present TimelessFlex, a framework for investigating chromatin state trajectories at promoters and enhancers and at promoter-enhancer pairs based on Hi-C information. TimelessFlex extends our previous approach Timeless, a Bayesian network for clustering multiple histone modification data sets at promoter and enhancer feature regions. We utilize time series ATAC-seq data measuring open chromatin to define promoters and enhancer candidates. We developed an expectation-maximization algorithm to assign promoters and enhancers to each other based on Hi-C interactions and jointly cluster their feature regions into paired chromatin state trajectories. We find jointly clustered promoter-enhancer pairs showing the same activation patterns on both sides but with a stronger trend at the enhancer side. While the promoter side remains accessible across the time series, the enhancer side becomes dynamically more open towards the gene activation time point. Promoter cluster patterns show strong correlations with gene expression signals, whereas Hi-C signals get only slightly stronger towards activation. The code of the framework is available at https://github.com/henriettemiko/TimelessFlex .ConclusionsTimelessFlex clusters time series histone modifications at promoter-enhancer pairs based on Hi-C and it can identify distinct chromatin states at promoter and enhancer feature regions and their changes over time.

Trimethylated histone H3 lysine 27 (H3K27me3) regulates gene repression, cell-fate determination and differentiation. We report that a conserved bromo-adjacent homology (BAH) module of BAHCC1 (BAHCC1BAH) 'recognizes' H3K27me3 specifically and enforces silencing of H3K27me3-demarcated genes in mammalian cells. Biochemical, structural and integrated chromatin immunoprecipitation-sequencing-based analyses demonstrate that direct readout of H3K27me3 by BAHCC1 is achieved through a hydrophobic trimethyl-L-lysine-binding 'cage' formed by BAHCC1BAH, mediating colocalization of BAHCC1 and H3K27me3-marked genes. BAHCC1 is highly expressed in human acute leukemia and interacts with transcriptional corepressors. In leukemia, depletion of BAHCC1, or disruption of the BAHCC1BAH-H3K27me3 interaction, causes derepression of H3K27me3-targeted genes that are involved in tumor suppression and cell differentiation, leading to suppression of oncogenesis. In mice, introduction of a germline mutation at Bahcc1 to disrupt its H3K27me3 engagement causes partial postnatal lethality, supporting a role in development. This study identifies an H3K27me3-directed transduction pathway in mammals that relies on a conserved BAH 'reader'.

BACKGROUND:During transcription, numerous transcription factors (TFs) bind to targets in a highly coordinated manner to control the gene expression. Alterations in groups of TF-binding profiles (i.e. "co-binding changes") can affect the co-regulating associations between TFs (i.e. "rewiring the co-regulator network"). This, in turn, can potentially drive downstream expression changes, phenotypic variation, and even disease. However, quantification of co-regulatory network rewiring has not been comprehensively studied. RESULTS:To address this, we propose DiNeR, a computational method to directly construct a differential TF co-regulation network from paired disease-to-normal ChIP-seq data. Specifically, DiNeR uses a graphical model to capture the gained and lost edges in the co-regulation network. Then, it adopts a stability-based, sparsity-tuning criterion -- by sub-sampling the complete binding profiles to remove spurious edges -- to report only significant co-regulation alterations. Finally, DiNeR highlights hubs in the resultant differential network as key TFs associated with disease. We assembled genome-wide binding profiles of 104 TFs in the K562 and GM12878 cell lines, which loosely model the transition between normal and cancerous states in chronic myeloid leukemia (CML). In total, we identified 351 significantly altered TF co-regulation pairs. In particular, we found that the co-binding of the tumor suppressor BRCA1 and RNA polymerase II, a well-known transcriptional pair in healthy cells, was disrupted in tumors. Thus, DiNeR successfully extracted hub regulators and discovered well-known risk genes. CONCLUSIONS:Our method DiNeR makes it possible to quantify changes in co-regulatory networks and identify alterations to TF co-binding patterns, highlighting key disease regulators. Our method DiNeR makes it possible to quantify changes in co-regulatory networks and identify alterations to TF co-binding patterns, highlighting key disease regulators.

The organization of chromatin into higher order structures is essential for chromosome segregation, the repair of DNA-damage, and the regulation of gene expression. Using Micro-C XL to detect chromosomal interactions, we observed the pervasive presence of cohesin-dependent loops with defined positions throughout the genome of budding yeast, as seen in mammalian cells. In early S phase, cohesin stably binds to cohesin associated regions (CARs) genome-wide. Subsequently, positioned loops accumulate with CARs at the bases of the loops. Cohesin regulators Wpl1 and Pds5 alter the levels and distribution of cohesin at CARs, changing the pattern of positioned loops. From these observations, we propose that cohesin with loop extrusion activity is stopped by preexisting CAR-bound cohesins, generating positioned loops. The patterns of loops observed in a population of wild-type and mutant cells can be explained by this mechanism, coupled with a heterogeneous residency of cohesin at CARs in individual cells.

Rationale: Puerto Ricans have the highest childhood asthma prevalence in the United States (23.6%); however, the etiology is uncertain.Objectives: In this study, we sought to uncover the genetic architecture of lung function in Puerto Rican youth with and without asthma who were recruited from the island (n = 836).Methods: We used admixture-mapping and whole-genome sequencing data to discover genomic regions associated with lung function. Functional roles of the prioritized candidate SNPs were examined with chromatin immunoprecipitation sequencing, RNA sequencing, and expression quantitative trait loci data.Measurements and Main Results: We discovered a genomic region at 1q32 that was significantly associated with a 0.12-L decrease in the lung volume of exhaled air (95% confidence interval, -0.17 to -0.07; P = 6.62 × 10-8) with each allele of African ancestry. Within this region, two SNPs were expression quantitative trait loci of TMEM9 in nasal airway epithelial cells and MROH3P in esophagus mucosa. The minor alleles of these SNPs were associated with significantly decreased lung function and decreased TMEM9 gene expression. Another admixture-mapping peak was observed on chromosome 5q35.1, indicating that each Native American ancestry allele was associated with a 0.15-L increase in lung function (95% confidence interval, 0.08-0.21; P = 5.03 × 10-6). The region-based association tests identified four suggestive windows that harbored candidate rare variants associated with lung function.Conclusions: We identified common and rare genetic variants that may play a critical role in lung function among Puerto Rican youth. We independently validated an inflammatory pathway that could potentially be used to develop more targeted treatments and interventions for patients with asthma.

Marek's disease (MD) is a neoplastic disease of chickens caused by Marek's disease virus (MDV), a member of the subfamily Alphaherpesvirinae Like other alphaherpesviruses, MDV encodes a serine/threonine protein kinase, US3. The functions of US3 have been extensively studied in other alphaherpesviruses; however, the biological functions of MDV US3 and its substrates have not been studied in detail. In this study, we investigated potential cellular pathways that are regulated by MDV US3 and identified chicken CREB (chCREB) as a substrate of MDV US3. We show that wild-type MDV US3, but not kinase-dead US3 (US3-K220A), increases CREB phosphorylation, leading to recruitment of phospho-CREB (pCREB) to the promoter of the CREB-responsive gene and activation of CREB target gene expression. Using US3 deletion and US3 kinase-dead recombinant MDV, we identified US3-responsive MDV genes during infection and found that the majority of US3-responsive genes were located in the MDV repeat regions. Chromatin immunoprecipitation sequencing (ChIP-seq) studies determined that some US3-regulated genes colocalized with Meq (an MDV-encoded oncoprotein) recruitment sites. Chromatin immunoprecipitation-PCR (ChIP-PCR) further confirmed Meq binding to the ICP4/LAT region, which is also regulated by US3. Furthermore, biochemical studies demonstrated that MDV US3 interacts with Meq in transfected cells and MDV-infected chicken embryonic fibroblasts in a phosphorylation-dependent manner. Finally, in vitro kinase studies revealed that Meq is a US3 substrate. MDV US3 thus acts as an upstream kinase of the CREB signaling pathway to regulate the transcription function of the CREB/Meq heterodimer, which targets cellular and viral gene expression.IMPORTANCE MDV is a potent oncogenic herpesvirus that induces T-cell lymphoma in infected chickens. Marek's disease continues to have a significant economic impact on the poultry industry worldwide. US3 encoded by alphaherpesviruses is a multifunctional kinase involved in the regulation of various cellular pathways. Using an MDV genome quantitative reverse transcriptase PCR (qRT-PCR) array and chromatin immunoprecipitation, we elucidated the role of MDV US3 in viral and cellular gene regulation. Our results provide insights into how viral kinase regulates host cell signaling pathways to activate both viral and host gene expression. This is an important step toward understanding host-pathogen interaction through activation of signaling cascades.

In a substantial fraction of cancers TERT promoter (TERTp) mutations drive expression of the catalytic subunit of telomerase, contributing to their proliferative immortality. We conducted a pan-cancer analysis of cell lines and find a TERTp mutation expression signature dominated by epithelial-to-mesenchymal transition and MAPK signaling. These data indicate that TERTp mutants are likely to generate distinctive tumor microenvironments and intercellular interactions. Analysis of high-throughput screening tests of 546 small molecules on cell line growth indicated that TERTp mutants displayed heightened sensitivity to specific drugs, including RAS pathway inhibitors, and we found that inhibition of MEK1 and 2, key RAS/MAPK pathway effectors, inhibited TERT mRNA expression. Consistent with an enrichment of mesenchymal states in TERTp mutants, cell lines and some patient tumors displayed low expression of the central adherens junction protein E-cadherin, and we provide evidence that its expression in these cells is regulated by MEK1/2. Several mesenchymal transcription factors displayed elevated expression in TERTp mutants including ZEB1 and 2, TWIST1 and 2, and SNAI1. Of note, the developmental transcription factor SNAI2/SLUG was conspicuously elevated in a significant majority of TERTp-mutant cell lines, and knock-down experiments suggest that it promotes TERT expression. IMPLICATIONS: Cancers harboring TERT promoter mutations are often more lethal, but the basis for this higher mortality remains unknown. Our study identifies that TERTp mutants, as a class, associate with a distinct gene and protein expression signature likely to impact their biological and clinical behavior and provide new directions for investigating treatment approaches for these cancers.

Understanding the aberrant transcriptional landscape of neuroblastoma is necessary to provide insight to the underlying influences of the initiation, progression and persistence of this developmental cancer. Here, we present chromatin immunoprecipitation sequencing (ChIP-Seq) data for the oncogenic transcription factors, MYCN and MYC, as well as regulatory histone marks H3K4me1, H3K4me3, H3K27Ac, and H3K27me3 in ten commonly used human neuroblastoma-derived cell line models. In addition, for all of the profiled cell lines we provide ATAC-Seq as a measure of open chromatin. We validate specificity of global MYCN occupancy in MYCN amplified cell lines and functional redundancy of MYC occupancy in MYCN non-amplified cell lines. Finally, we show with H3K27Ac ChIP-Seq that these cell lines retain expression of key neuroblastoma super-enhancers (SE). We anticipate this dataset, coupled with available transcriptomic profiling on the same cell lines, will enable the discovery of novel gene regulatory mechanisms in neuroblastoma.

Abstract Gene expression in malaria parasites is subject to various layers of regulation, including histone post-translational modifications (PTMs). Gene regulatory mechanisms have been extensively studied during the main developmental stages of Plasmodium parasites inside erythrocytes, from the ring stage following invasion to the schizont stage leading up to egress. However, gene regulation in merozoites that mediate the transition from one host cell to the next is an understudied area of parasite biology. Here, we sought to characterize gene expression and the corresponding histone PTM landscape during this stage of the parasite lifecycle through RNA-seq and ChIP-seq on P. falciparum blood stage schizonts, merozoites, and rings, as well as P. berghei liver stage merozoites. In both hepatic and erythrocytic merozoites, we identified a subset of genes with a unique histone PTM profile characterized by a region of H3K4me3 depletion in their promoter. These genes were upregulated in hepatic and erythrocytic merozoites and rings, had roles in protein export, translation, and host cell remodeling, and shared a DNA motif. These results indicate that similar regulatory mechanisms may underlie merozoite formation in the liver and blood stages. We also observed that H3K4me2 was deposited in gene bodies of gene families encoding variant surface antigens in erythrocytic merozoites, which may facilitate switching of gene expression between different members of these families. Finally, H3K18me and H2K27me were uncoupled from gene expression and were enriched around the centromeres in erythrocytic schizonts and merozoites, suggesting potential roles in the maintenance of chromosomal organization during schizogony. Together, our results demonstrate that extensive changes in gene expression and histone landscape occur during the schizont-to-ring transition to facilitate productive erythrocyte infection. The dynamic remodeling of the transcriptional program in hepatic and erythrocytic merozoites makes this stage attractive as a target for novel anti-malarial drugs that may have activity against both the liver and blood stages.

Abstract Background Identification of the motifs bound by a transcription factor (TF) is important to reveal the function of TF. Previously, we built a transcription factor centered yeast one hybrid (TF-Centered Y1H) that could identify the motifs bound by a target TF. However, that method was difficult to comprehensively identify all the motifs bound by a TF. Results Here, we build an improved TF-Centered Y1H to comprehensively determine the motifs bound by a target TF. Recombination-mediated cloning in yeast was performed to construct a saturated prey library that contains 7 random base insertions. After TF-Centered Y1H screening, all the positive clones were pooled together to isolate pHIS2 vector. The insertion regions of pHIS2 were PCR amplified and the PCR product was subjected to high-throughput sequencing. The insertion sequences were then retrieved and analyzed using MEME program to identify the potential motifs bound by the TF. Using this technology, we studied the motifs bound by an ethylene-responsive factor (BpERF2) from birch. In total, 22 conserved motifs were identified, and most of them are novel cis-acting elements. Both the yeast one hybrid and electrophoretic mobility shift assay verified that the obtained motifs could be bound by BpERF2. In addition, chromatin immunoprecipitation (ChIP) study further suggested that the identified motifs can be bound by BpERF2 in cells of birch. These results together suggested that this technology is reliable and has biological significance. Conclusion This method will have wide application in DNA-protein interaction studies.

Abstract Chromatin immunoprecipitation (ChIP) is an antibody-based approach that is frequently utilized in chromatin biology and epigenetics. The challenge in experimental variability by unpredictable nature of usable input amounts from samples and undefined antibody titer in ChIP reaction still remains to be addressed. Here, we introduce a simple and quick method to quantify chromatin inputs and demonstrate its utility for normalizing antibody amounts to the optimal titer in individual ChIP reactions. For a proof of concept, we utilized ChIP-seq validated antibodies against the key enhancer mark, acetylation of histone H3 on lysine 27 (H3K27ac), in the experiments. The results indicate that the titration-based normalization of antibody amounts improves assay outcomes including the consistency among samples both within and across experiments for a broad range of input amounts.

Summary: Pioneer transcription factors (TFs) can directly establish higher-order chromatin interactions to instruct gene transcription. Here, we present a protocol for capturing TF-mediated 3D chromatin interactions using affinity tag-based bridge linker (BL)-Hi-chromatin immunoprecipitation (HiChIP). We describe steps for constructing FLAG-tagged TF, performing BL-HiChIP, and preparing the library. We then detail procedures for sequencing, data analysis, and quality control. This protocol has potential applications in 3D chromatin analysis centered on any specific TF in any type of cells without the need of optimal antibodies.For complete details on the use and execution of this protocol, please refer to Ren et al. (2022).1 : Publisher’s note: Undertaking any experimental protocol requires adherence to local institutional guidelines for laboratory safety and ethics.