Search

Insects use seasonal diapause as an alternative strategy to endure adverse seasons. This developmental trajectory is induced by environmental cues like short‐day lengths in late summer and early fall, but how insects measure day length is unknown. The circadian clock has been implicated in regulating photoperiodic or seasonal responses in many insects, including the Northern house mosquito, Culex pipiens, which enters adult diapause. To investigate the potential control of diapause by circadian control, we employed ChIP‐sequencing to identify the downstream targets of a circadian transcription factor, PAR domain protein 1 (PDP1), that contribute to the hallmark features of diapause. We identified the nearest genes in a 10 kb region of the anticipated PDP1 binding sites, listed prospective targets and searched for PDP1‐specific binding sites. By examining the functional relevance to diapause‐specific behaviours and modifications such as metabolic pathways, lifespan extension, cell cycle regulation and stress tolerance, eight genes were selected as targets and validated using ChIP‐qPCR. In addition, qRT‐PCR demonstrated that the mRNA abundance of PDP1 targets increased in the heads of diapausing females during the middle of the scotophase (ZT17) compared with the early photophase (ZT1), in agreement with the peak and trough of PDP1 abundance. Thus, our investigation uncovered the mechanism by which PDP1 might generate a diapause phenotype in insects. [ABSTRACT FROM AUTHOR]

Summary: Rice blast, caused by Magnaporthe oryzae, significantly impacts grain yield, necessitating the identification of broad‐spectrum resistance genes and their functional mechanisms for disease‐resistant crop breeding. Here, we report that rice with knockdown OsHDAC1 gene expression displays enhanced broad‐spectrum blast resistance without effects on plant height and tiller numbers compared to wild‐type rice, while rice overexpressing OsHDAC1 is more susceptible to M. oryzae. We identify a novel blast resistance transcription factor, OsGRAS30, which genetically acts upstream of OsHDAC1 and interacts with OsHDAC1 to suppress its enzymatic activity. This inhibition increases the histone H3K27ac level, thereby boosting broad‐spectrum blast resistance. Integrating genome‐wide mapping of OsHDAC1 and H3K27ac targets with RNA sequencing analysis unveils how OsHDAC1 mediates the expression of OsSSI2, OsF3H, OsRLR1 and OsRGA5 to regulate blast resistance. Our findings reveal that the OsGRAS30–OsHDAC1 module is critical to rice blast control. Therefore, targeting either OsHDAC1 or OsGRAS30 offers a promising approach for enhancing crop blast resistance. [ABSTRACT FROM AUTHOR]

Background: Standard ChIP-seq and RNA-seq processing pipelines typically disregard sequencing reads whose origin is ambiguous ("multimappers"). This usual practice has potentially important consequences for the functional interpretation of the data: genomic elements belonging to clusters composed of highly similar members are left unexplored. Results: In particular, disregarding multimappers leads to the underrepresentation in epigenetic studies of recently active transposable elements, such as AluYa5, L1HS and SVAs. Furthermore, this common strategy also has implications for transcriptomic analysis: members of repetitive gene families, such the ones including major histocompatibility complex (MHC) class I and II genes, are under-quantified. Conclusion: Revealing inherent biases that permeate routine tasks such as functional enrichment analysis, our results underscore the urgency of broadly adopting multimapper-aware bioinformatic pipelines –currently restricted to specific contexts or communities– to ensure the reliability of genomic and transcriptomic studies. [ABSTRACT FROM AUTHOR]

Summary: Microglia, the major population of brain‐resident macrophages, are now recognized as a heterogeneous population comprising several cell subtypes with different (so far mostly supposed) functions in health and disease. A number of studies have performed molecular characterization of these different microglial activation states over the last years making use of "omics" technologies, that is transcriptomics, proteomics and, less frequently, epigenomics profiling. These approaches offer the possibility to identify disease mechanisms, discover novel diagnostic biomarkers, and develop new therapeutic strategies. Here, we focus on epigenetic profiling as a means to understand microglial immune responses beyond what other omics methods can offer, that is, revealing past and present molecular responses, gene regulatory networks and potential future response trajectories, and defining cell subtype‐specific disease relevance through mapping non‐coding genetic variants. We review the current knowledge in the field regarding epigenetic regulation of microglial identity and function, provide an exemplary analysis that demonstrates the advantages of performing joint transcriptomic and epigenomic profiling of single microglial cells and discuss how comprehensive epigenetic analyses may enhance our understanding of microglial pathophysiology. [ABSTRACT FROM AUTHOR]

To identify the inactive genes in cataract lenses and explore their function in lens epithelial cells (LECs). Lens epithelium samples obtained from both age-related cataract (ARC) patients and normal donors were subjected to two forms of histone H3 immunoprecipitation: H3K9ac and H3K27me3 chromatin immunoprecipitation (ChIP), followed by ChIP-seq. The intersection set of "active genes in normal controls" and "repressed genes in cataract lenses" was identified. To validate the role of a specific gene, ETV1, within this set, quantitative polymerase chain reaction (qPCR), western blot, and immunofluorescence were performed using clinical lens epithelium samples. Small interference RNA (siRNA) was utilized to reduce the mRNA level of ETV1 in cultured LECs. Following this, transwell assay and western blot was performed to examine the migration ability of the cells. Furthermore, RNA-seq analysis was conducted on both cell samples with ETV1 knockdown and control cells. Additionally, the expression level of ETV1 in LECs was examined using qPCR under H2O2 treatment. Six genes were identified in the intersection set of "active genes in normal controls" and "repressed genes in ARC lenses". Among these genes, ETV1 showed the most significant fold-change decrease in the cataract samples compared to the control samples. After ETV1 knockdown by siRNA in cultured LECs, reduced cell migration was observed, along with a decrease in the expression of β-Catenin and Vimentin, two specific genes associated with cell migration. In addition, under the oxidative stress induced by H2O2 treatment, the expression level of ETV1 in LECs significantly decreased. Based on the findings of this study, it can be concluded that ETV1 is significantly reduced in human ARC lenses. The repression of ETV1 in ARC lenses appears to contribute to the disrupted differentiation of lens epithelium, which is likely caused by the inhibition of both cell differentiation and migration processes. [ABSTRACT FROM AUTHOR]

Background: Cervical cancer is among the highest-ranking types of cancer worldwide, with human papillomavirus (HPV) as the agent driving the malignant process. One aspect of the infection's evolution is given by epigenetic modifications, mainly DNA methylation and chromatin alteration. These processes are guided by several chromatin remodeling complexes, including NuRD. The purpose of this study was to evaluate the genome-wide binding patterns of the NuRD complex components (MBD2 and MBD3) in the presence of active HPV16 E6 and E7 oncogenes and to determine the potential of identified genes through an experimental model to differentiate between cervical precursor lesions, with the aim of establishing their utility as biomarkers. Methods: The experimental model was built using the CaSki cell line and shRNA for E6 and E7 HPV16 silencing, ChIP-seq, qRT-PCR, and Western blot analyses. Selected genes' expression was also assessed in patients. Results: Several genes have been identified to exhibit altered transcriptional activity due to the influence of HPV16 E6/E7 viral oncogenes acting through the MBD2/MBD3 NuRD complex, linking them to viral infection and cervical oncogenesis. Conclusions: The impacted genes primarily play roles in governing gene transcription, mRNA processing, and regulation of translation. Understanding these mechanisms offers valuable insights into the process of HPV-induced oncogenesis. [ABSTRACT FROM AUTHOR]

Objective: Rare study of the non-coding and regulatory regions of the genome limits our ability to decode the mechanisms of fatty liver hemorrhage syndrome (FLHS) in chickens. Methods: Herein, we constructed the high-fat diet-induced FLHS chicken model to investigate the genome-wide active enhancers and transcriptome by H3K27ac target chromatin immunoprecipitation sequencing (ChIP-seq) and RNA sequencing (RNA-Seq) profiles of normal and FLHS liver tissues. Concurrently, an integrative analysis combining ChIP-seq with RNA-Seq and a comparative analysis with chicken FLHS, rat non-alcoholic fatty liver disease (NAFLD) and human NAFLD at the transcriptome level revealed the enhancer and super enhancer target genes and conservative genes involved in metabolic processes. Results: In total, 56 and 199 peak-genes were identified in upregulated peak-genes positively regulated by H3K27ac (Cor (peak-gene correlation) ≥0.5 and log2(FoldChange) ≥1) (PP) and downregulated peak-genes positively regulated by H3K27ac (Cor (peak-gene correlation) ≥0.5 and log2(FoldChange)≤-1) (PN), respectively; then we screened key regulatory targets mainly distributing in lipid metabolism (PCK1, APOA4, APOA1, INHBE) and apoptosis (KIT, NTRK2) together with MAPK and PPAR signaling pathway in FLHS. Intriguingly, PCK1 was also significantly covered in up-regulated super-enhancers (SEs), which further implied the vital role of PCK1 during the development of FLHS. Conclusion: Together, our studies have identified potential therapeutic biomarkers of PCK1 and elucidated novel insights into the pathogenesis of FLHS, especially for the epigenetic perspective. [ABSTRACT FROM AUTHOR]

In a previous report, we demonstrated that Cbx1, PurB and Sp3 inhibited cardiac muscle differentiation by increasing nucleosome density around cardiac muscle gene promoters. Since cardiac and skeletal muscle express many of the same proteins, we asked if Cbx1, PurB and Sp3 similarly regulated skeletal muscle differentiation. In a C2C12 model of skeletal muscle differentiation, Cbx1 and PurB knockdown increased myotube formation. In contrast, Sp3 knockdown inhibited myotube formation, suggesting that Sp3 played opposing roles in cardiac muscle and skeletal muscle differentiation. Consistent with this finding, Sp3 knockdown also inhibited various muscle-specific genes. The Cbx1, PurB and Sp3 proteins are believed to influence gene-expression in part by altering nucleosome position. Importantly, we developed a statistical approach to determine if changes in nucleosome positioning were significant and applied it to understanding the architecture of muscle-specific genes. Through this novel statistical approach, we found that during myogenic differentiation, skeletal muscle-specific genes undergo a set of unique nucleosome changes which differ significantly from those shown in commonly expressed muscle genes. While Sp3 binding was associated with nucleosome loss, there appeared no correlation with the aforementioned nucleosome changes. In summary, we have identified a novel role for Sp3 in skeletal muscle differentiation and through the application of quantifiable MNase-seq have discovered unique fingerprints of nucleosome changes for various classes of muscle genes during myogenic differentiation. [ABSTRACT FROM AUTHOR]

Background: Fatty liver hemorrhagic syndrome (FLHS) in the modern poultry industry is primarily caused by nutrition. Despite encouraging progress on FLHS, the mechanism through which nutrition influences susceptibility to FLHS is still lacking in terms of epigenetics. Results: In this study, we analyzed the genome-wide patterns of trimethylated lysine residue 27 of histone H3 (H3K27me3) enrichment by chromatin immunoprecipitation-sequencing (ChIP-seq), and examined its association with transcriptomes in healthy and FLHS hens. The study results indicated that H3K27me3 levels were increased in the FLHS hens on a genome-wide scale. Additionally, H3K27me3 was found to occupy the entire gene and the distant intergenic region, which may function as silencer-like regulatory elements. The analysis of transcription factor (TF) motifs in hypermethylated peaks has demonstrated that 23 TFs are involved in the regulation of liver metabolism and development. Transcriptomic analysis indicated that differentially expressed genes (DEGs) were enriched in fatty acid metabolism, amino acid, and carbohydrate metabolism. The hub gene identified from PPI network is fatty acid synthase (FASN). Combined ChIP-seq and transcriptome analysis revealed that the increased H3K27me3 and down-regulated genes have significant enrichment in the ECM-receptor interaction, tight junction, cell adhesion molecules, adherens junction, and TGF-beta signaling pathways. Conclusions: Overall, the trimethylation modification of H3K27 has been shown to have significant regulatory function in FLHS, mediating the expression of crucial genes associated with the ECM-receptor interaction pathway. This highlights the epigenetic mechanisms of H3K27me3 and provides insights into exploring core regulatory targets and nutritional regulation strategies in FLHS. [ABSTRACT FROM AUTHOR]

Background: There are two major genetic types of Epstein-Barr Virus (EBV): type 1 (EBV-1) and type 2 (EBV-2). EBV functions by manipulating gene expression in host B cells, using virus-encoded gene regulatory proteins including Epstein-Barr Nuclear Antigen 2 (EBNA2). While type 1 EBNA2 is known to interact with human transcription factors (hTFs) such as RBPJ, EBF1, and SPI1 (PU.1), type 2 EBNA2 shares only ~ 50% amino acid identity with type 1 and thus may have distinct binding partners, human genome binding locations, and functions. Results: In this study, we examined genome-wide EBNA2 binding in EBV-1 and EBV-2 transformed human B cells to identify shared and unique EBNA2 interactions with the human genome, revealing thousands of type-specific EBNA2 ChIP-seq peaks. Computational predictions based on hTF motifs and subsequent ChIP-seq experiments revealed that both type 1 and 2 EBNA2 co-occupy the genome with SPI1 and AP-1 (BATF and JUNB) hTFs. However, type 1 EBNA2 showed preferential co-occupancy with EBF1, and type 2 EBNA2 preferred RBPJ. These differences in hTF co-occupancy revealed possible mechanisms underlying type-specific gene expression of known EBNA2 human target genes: MYC (shared), CXCR7 (type 1 specific), and CD21 (type 2 specific). Both type 1 and 2 EBNA2 binding events were enriched at systemic lupus erythematosus (SLE) and multiple sclerosis (MS) risk loci, while primary biliary cholangitis (PBC) risk loci were specifically enriched for type 2 peaks. Conclusions: This study reveals extensive type-specific EBNA2 interactions with the human genome, possible differences in EBNA2 interaction partners, and a possible new role for type 2 EBNA2 in autoimmune disorders. Our results highlight the importance of considering EBV type in the control of human gene expression and disease-related investigations. [ABSTRACT FROM AUTHOR]

Background: Switchgrass (Panicum virgatum L.) is a warm-season perennial (C4) grass identified as an important biofuel crop in the United States. It is well adapted to the marginal environment where heat and moisture stresses predominantly affect crop growth. However, the underlying molecular mechanisms associated with heat and drought stress tolerance still need to be fully understood in switchgrass. The methylation of H3K4 is often associated with transcriptional activation of genes, including stress-responsive. Therefore, this study aimed to analyze genome-wide histone H3K4-tri-methylation in switchgrass under heat, drought, and combined stress. Results: In total, ~ 1.3 million H3K4me3 peaks were identified in this study using SICER. Among them, 7,342; 6,510; and 8,536 peaks responded under drought (DT), drought and heat (DTHT), and heat (HT) stresses, respectively. Most DT and DTHT peaks spanned 0 to + 2000 bases from the transcription start site [TSS]. By comparing differentially marked peaks with RNA-Seq data, we identified peaks associated with genes: 155 DT-responsive peaks with 118 DT-responsive genes, 121 DTHT-responsive peaks with 110 DTHT-responsive genes, and 175 HT-responsive peaks with 136 HT-responsive genes. We have identified various transcription factors involved in DT, DTHT, and HT stresses. Gene Ontology analysis using the AgriGO revealed that most genes belonged to biological processes. Most annotated peaks belonged to metabolite interconversion, RNA metabolism, transporter, protein modifying, defense/immunity, membrane traffic protein, transmembrane signal receptor, and transcriptional regulator protein families. Further, we identified significant peaks associated with TFs, hormones, signaling, fatty acid and carbohydrate metabolism, and secondary metabolites. qRT-PCR analysis revealed the relative expressions of six abiotic stress-responsive genes (transketolase, chromatin remodeling factor-CDH3, fatty-acid desaturase A, transmembrane protein 14C, beta-amylase 1, and integrase-type DNA binding protein genes) that were significantly (P < 0.05) marked during drought, heat, and combined stresses by comparing stress-induced against un-stressed and input controls. Conclusion: Our study provides a comprehensive and reproducible epigenomic analysis of drought, heat, and combined stress responses in switchgrass. Significant enrichment of H3K4me3 peaks downstream of the TSS of protein-coding genes was observed. In addition, the cost-effective experimental design, modified ChIP-Seq approach, and analyses presented here can serve as a prototype for other non-model plant species for conducting stress studies. [ABSTRACT FROM AUTHOR]

Background & Aims: Acinar-to-ductal metaplasia (ADM) serves as a precursor event in the development of pancreatic ductal adenocarcinoma (PDAC) upon constitutive environmental and genetical stress. While the role of ADM in PDAC progression has been established, the molecular mechanisms underlying human ADM remain elusive. We previously demonstrated the induction of ADM in human acinar cells through the transforming growth factor beta (TGFβ) signaling pathway. We aim to investigate the interaction between TGFβ and Hippo pathways in mediating ADM. Methods: RNA-sequencing was conducted on sorted normal primary human acinar, ductal, and AD (acinar cells that have undergone ADM) cells. ATAC-seq analysis was utilized to reveal the chromatin accessibility in these three cell types. ChIP-Seq of YAP1, SMAD4, and H3K27ac was performed to identify the gene targets of YAP1 and SMAD4. The role of YAP1/TAZ in ADM-driven cell proliferation, as well as in oncogenic KRAS driven proliferation, was assessed using sphere formation assay. Results: AD cells have a unique transcription profile, with upregulated genes in open chromatin states in acinar cells. YAP1 and SMAD4 co-occupy the loci of ADM-related genes, including PROM1, HES1, and MMP7, co-regulating biological functions such as cell adhesion, cell migration, and inflammation. Overexpression of YAP1/TAZ promoted acinar cell proliferation but still required the TGFβ pathway. YAP1/TAZ were also crucial for TGFβ-induced sphere formation and were necessary for KRAS-induced proliferation. Conclusions: Our study reveals the intricate transition between acinar and AD states in human pancreatic tissues. It unveils the complex interaction between the Hippo and TGF-β pathways during ADM, highlighting the pivotal role of YAP1/TAZ and SMAD4 in PDAC initiation. [ABSTRACT FROM AUTHOR]

Simple Summary: CIC-DUX4-rearranged sarcoma (CDS) is a rare and aggressive soft tissue tumor predominantly affecting young adults for which the only oncogenic driver is the CIC-DUX4 fusion protein, resulting from chromosomal rearrangements. The molecular mechanisms by which CIC-DUX4 drives CDS development remain largely unknown, hindering the development of targeted therapies. This study conducted genome-wide profiling of CIC-DUX4 DNA occupancy and associated chromatin states in CDS cell models and tumors. Our findings reveal CIC-DUX4 as a potent transcriptional activator at its binding sites, facilitated by direct interaction with the acetyltransferase p300. Significantly, inhibiting p300 impedes CDS tumor cell proliferation, suggesting a potential therapeutic avenue. This research enhances understanding of CIC-DUX4-mediated transcriptional regulation and proposes targeting p300 as a promising strategy for the clinical management of CDS. CIC-DUX4-rearranged sarcoma (CDS) is a rare and aggressive soft tissue tumor that occurs most frequently in young adults. The key oncogenic driver of this disease is the expression of the CIC-DUX4 fusion protein as a result of chromosomal rearrangements. CIC-DUX4 displays chromatin binding properties, and is therefore believed to function as an aberrant transcription factor. However, the chromatin remodeling events induced by CIC-DUX4 are not well understood, limiting our ability to identify new mechanism-based therapeutic strategies for these patients. Here, we generated a genome-wide profile of CIC-DUX4 DNA occupancy and associated chromatin states in human CDS cell models and primary tumors. Combining chromatin profiling, proximity ligation assays, as well as genetic and pharmacological perturbations, we show that CIC-DUX4 operates as a potent transcriptional activator at its binding sites. This property is in contrast with the repressive function of the wild-type CIC protein, and is mainly mediated through the direct interaction of CIC-DUX4 with the acetyltransferase p300. In keeping with this, we show p300 to be essential for CDS tumor cell proliferation; additionally, we find its pharmacological inhibition to significantly impact tumor growth in vitro and in vivo. Taken together, our study elucidates the mechanisms underpinning CIC-DUX4-mediated transcriptional regulation. [ABSTRACT FROM AUTHOR]

Summary: Gossypol and the related terpenoids are stored in the pigment gland to protect cotton plants from biotic stresses, but little is known about the synthetic sites of these metabolites.Here, we showed that GoPGF, a key gene regulating gland formation, was expressed in gland cells and roots. The chromatin immunoprecipitation sequencing (ChIP‐seq) analysis demonstrated that GoPGF targets GhJUB1 to regulate gland morphogenesis.RNA‐sequencing (RNA‐seq) showed high accumulation of gossypol biosynthetic genes in gland cells. Moreover, integrated analysis of the ChIP‐seq and RNA‐seq data revealed that GoPGF binds to the promoter of several gossypol biosynthetic genes. The cotton callus overexpressing GoPGF had dramatically increased the gossypol levels, indicating that GoPGF can directly activate the biosynthesis of gossypol. In addition, the gopgf mutant analysis revealed the existence of both GoPGF‐dependent and ‐independent regulation of gossypol production in cotton roots.Our study revealed that the pigment glands are synthetic sites of gossypol in aerial parts of cotton and that GoPGF plays a dual role in regulating gland morphogenesis and gossypol biosynthesis. The study provides new insights for exploring the complex relationship between glands and the metabolites they store in cotton and other plant species. [ABSTRACT FROM AUTHOR]

DNA double-strand breaks (DSBs) are among the most serious types of DNA damage, causing mutations and chromosomal rearrangements. In eukaryotes, DSBs are immediately repaired in coordination with chromatin remodeling for the deposition of DSB-related histone modifications and variants. To elucidate the details of DSB-dependent chromatin remodeling throughout the genome, artificial DSBs need to be reproducibly induced at various genomic loci. Recently, a comprehensive method for elucidating chromatin remodeling at multiple DSB loci via chemically induced expression of a restriction enzyme was developed in mammals. However, this DSB induction system is unsuitable for investigating chromatin remodeling during and after DSB repair, and such an approach has not been performed in plants. Here, we established a transgenic Arabidopsis plant harboring a restriction enzyme gene Sbf I driven by a heat-inducible promoter. Using this transgenic line, we performed chromatin immunoprecipitation followed by deep sequencing (ChIP-seq) of histones H4K16ac and H2A.Z and investigated the dynamics of these histone marks around the endogenous 623 Sbf I recognition sites. We also precisely quantified DSB efficiency at all cleavage sites using the DNA resequencing data obtained by the ChIP-seq procedure. From the results, Sbf I–induced DSBs were detected at 360 loci, which induced the transient deposition of H4K16ac and H2A.Z around these regions. Interestingly, we also observed the co-localization of H4K16ac and H2A.Z at some DSB loci. Overall, DSB-dependent chromatin remodeling was found to be highly conserved between plants and animals. These findings provide new insights into chromatin remodeling that occurs in response to DSBs in Arabidopsis. [ABSTRACT FROM AUTHOR]

Background: The specific recognition of a DNA locus by a given transcription factor is a widely studied issue. It is generally agreed that the recognition can be influenced not only by the binding motif but by the larger context of the binding site. In this work, we present a novel heuristic algorithm that can reconstruct the unique binding sites captured in a sequencing experiment without using the reference genome. Results: We present PAPerFly, the Partial Assembly-based Peak Finder, a tool for the binding site and binding context reconstruction from the sequencing data without any prior knowledge. This tool operates without the need to know the reference genome of the respective organism. We employ algorithmic approaches that are used during genome assembly. The proposed algorithm constructs a de Bruijn graph from the sequencing data. Based on this graph, sequences and their enrichment are reconstructed using a novel heuristic algorithm. The reconstructed sequences are aligned and the peaks in the sequence enrichment are identified. Our approach was tested by processing several ChIP-seq experiments available in the ENCODE database and comparing the results of Paperfly and standard methods. Conclusions: We show that PAPerFly, an algorithm tailored for experiment analysis without the reference genome, yields better results than an aggregation of ChIP-seq agnostic tools. Our tool is freely available at https://github.com/Caeph/paperfly/ or on Zenodo (https://doi.org/10.5281/zenodo.7116424). [ABSTRACT FROM AUTHOR]

Background: Estrogen receptor beta (ERβ, Esr2) plays a pivotal role in folliculogenesis and ovulation, yet its exact mechanism of action is mainly uncharacterized. Results: We here performed ERβ ChIP-sequencing of mouse ovaries followed by complementary RNA-sequencing of wild-type and ERβ knockout ovaries. By integrating the ERβ cistrome and transcriptome, we identified its direct target genes and enriched biological functions in the ovary. This demonstrated its strong impact on genes regulating organism development, cell migration, lipid metabolism, response to hypoxia, and response to estrogen. Cell-type deconvolution analysis of the bulk RNA-seq data revealed a decrease in luteal cells and an increased proportion of theca cells and a specific type of cumulus cells upon ERβ loss. Moreover, we identified a significant overlap with the gene regulatory network of liver receptor homolog 1 (LRH-1, Nr5a2) and showed that ERβ and LRH-1 extensively bound to the same chromatin locations in granulosa cells. Using ChIP-reChIP, we corroborated simultaneous ERβ and LRH-1 co-binding at the ERβ-repressed gene Greb1 but not at the ERβ-upregulated genes Cyp11a1 and Fkbp5. Transactivation assay experimentation further showed that ERβ and LRH-1 can inhibit their respective transcriptional activity at classical response elements. Conclusions: By characterizing the genome-wide endogenous ERβ chromatin binding, gene regulations, and extensive crosstalk between ERβ and LRH-1, along with experimental corroborations, our data offer genome-wide mechanistic underpinnings of ovarian physiology and fertility. [ABSTRACT FROM AUTHOR]

The global cancer burden remains high; thus, a better understanding of the molecular mechanisms driving carcinogenesis is needed to improve current prevention and treatment options. We previously detected the ZNF643/ZFP69B gene upregulated in multiple tumors, and we speculated it may play a role in tumor biology. To test this hypothesis, we employed TCGA-centered databases to correlate ZNF643 status with various clinicopathological parameters. We also performed RNA-seq analysis and in vitro studies assessing cancer cell phenotypes, and we searched for ZNF643-bound genomic loci. Our data indicated higher levels of ZNF643 in most analyzed tumors compared to normal samples, possibly due to copy number variations. ZNF643 mRNA correlated with diverse molecular and immune subtypes and clinicopathological features (tumor stage, grade, patient survival). RNA-seq analysis revealed that ZNF643 silencing triggers the deregulation of the genes implicated in various cancer-related processes, such as growth, adhesion, and immune system. Moreover, we observed that ZNF643 positively influences cell cycle, migration, and invasion. Finally, our ChIP-seq analysis indicated that the genes associated with ZNF643 binding are linked to adhesion and immune signaling. In conclusion, our data confirm the oncogenic properties of ZNF643 and pinpoint its impact on cell adhesion and immune processes. [ABSTRACT FROM AUTHOR]

Rhabdomyosarcoma (RMS) is a pediatric malignancy of the muscle with characteristics of cells blocked in differentiation. NOTCH1 is an oncogene that promotes self-renewal and blocks differentiation in the fusion negative-RMS sub-type. However, how NOTCH1 expression is transcriptionally maintained in tumors is unknown. Analyses of SNAI2 and CTCF chromatin binding and HiC analyses revealed a conserved SNAI2/CTCF overlapping peak downstream of the NOTCH1 locus marking a sub-topologically associating domain (TAD) boundary. Deletion of the SNAI2-CTCF peak showed that it is essential for NOTCH1 expression and viability of FN-RMS cells. Reintroducing constitutively activated NOTCH1-ΔE in cells with the SNAI2-CTCF peak deleted restored cell-viability. Ablation of SNAI2 using CRISPR/Cas9 reagents resulted in the loss of majority of RD and SMS-CTR FN-RMS cells. However, the few surviving clones that repopulate cultures have recovered NOTCH1. Cells that re-establish NOTCH1 expression after SNAI2 ablation are unable to differentiate robustly as SNAI2 shRNA knockdown cells; yet, SNAI2-ablated cells continued to be exquisitely sensitive to ionizing radiation. Thus, we have uncovered a novel mechanism by which SNAI2 and CTCF maintenance of a sub-TAD boundary promotes rather than represses NOTCH1 expression. Further, we demonstrate that SNAI2 suppression of apoptosis post-radiation is independent of SNAI2/NOTCH1 effects on self-renewal and differentiation. [ABSTRACT FROM AUTHOR]

The multifunctional adenovirus E1B-55K oncoprotein can induce cell transformation in conjunction with adenovirus E1A gene products. Previous data from transient expression studies and in vitro experiments suggest that these growth-promoting activities correlate with E1B-55K-mediated transcriptional repression of p53-targeted genes. Here, we analyzed genome-wide occupancies and transcriptional consequences of species C5 and A12 E1B-55Ks in transformed mammalian cells by combinatory ChIP and RNA-seq analyses. E1B-55K-mediated repression correlates with tethering of the viral oncoprotein to p53-dependent promoters via DNA-bound p53. Moreover, we found that E1B-55K also interacts with and represses transcription of numerous p53-independent genes through interactions with transcription factors that play central roles in cancer and stress signaling. Our results demonstrate that E1B-55K oncoproteins function as promiscuous transcriptional repressors of both p53-dependent and -independent genes and further support the model that manipulation of cellular transcription is central to adenovirus-induced cell transformation and oncogenesis. [ABSTRACT FROM AUTHOR]

Background: Transcription factors (TFs) exhibit heterogeneous DNA-binding specificities in individual cells and whole organisms under natural conditions, and de novo motif discovery usually provides multiple motifs, even from a single chromatin immunoprecipitation-sequencing (ChIP-seq) sample. Despite the accumulation of ChIP-seq data and ChIP-seq-derived motifs, the diversity of DNA-binding specificities across different TFs and cell types remains largely unexplored. Results: Here, we applied MOCCS2, our k-mer-based motif discovery method, to a collection of human TF ChIP-seq samples across diverse TFs and cell types, and systematically computed profiles of TF-binding specificity scores for all k-mers. After quality control, we compiled a set of TF-binding specificity score profiles for 2,976 high-quality ChIP-seq samples, comprising 473 TFs and 398 cell types. Using these high-quality samples, we confirmed that the k-mer-based TF-binding specificity profiles reflected TF- or TF-family dependent DNA-binding specificities. We then compared the binding specificity scores of ChIP-seq samples with the same TFs but with different cell type classes and found that half of the analyzed TFs exhibited differences in DNA-binding specificities across cell type classes. Additionally, we devised a method to detect differentially bound k-mers between two ChIP-seq samples and detected k-mers exhibiting statistically significant differences in binding specificity scores. Moreover, we demonstrated that differences in the binding specificity scores between k-mers on the reference and alternative alleles could be used to predict the effect of variants on TF binding, as validated by in vitro and in vivo assay datasets. Finally, we demonstrated that binding specificity score differences can be used to interpret disease-associated non-coding single-nucleotide polymorphisms (SNPs) as TF-affecting SNPs and provide candidates responsible for TFs and cell types. Conclusions: Our study provides a basis for investigating the regulation of gene expression in a TF-, TF family-, or cell-type-dependent manner. Furthermore, our differential analysis of binding-specificity scores highlights noncoding disease-associated variants in humans. [ABSTRACT FROM AUTHOR]

Hox genes encode Homeodomain-containing transcription factors, which specify segmental identities along the anterior–posterior axis. Functional changes in Hox genes have been directly implicated in the evolution of body plans across the metazoan lineage. The Hox protein Ultrabithorax (Ubx) is expressed and required in developing third thoracic (T3) segments in holometabolous insects studied so far, particularly, of the order Coleoptera, Lepidoptera and Diptera. Ubx function is key to specify differential development of the second (T2) and T3 thoracic segments in these insects. While Ubx is expressed in the third thoracic segment in developing larvae of Hymenopteran Apismellifera, the morphological differences between T2 and T3 are subtle. To identify evolutionary changes that are behind the differential function of Ubx in Drosophila and Apis, which are diverged for more than 350 million years, we performed comparative analyses of genome wide Ubx-binding sites between these two insects. Our studies reveal that a motif with a TAAAT core is a preferred binding site for Ubx in Drosophila, but not in Apis. Biochemical and transgenic assays suggest that in Drosophila, the TAAAT core sequence in the Ubx binding sites is required for Ubx-mediated regulation of two of its target genes studied here; CG13222, a gene that is normally upregulated by Ubx and vestigial (vg), whose expression is repressed by Ubx in T3. Interestingly, changing the TAAT site to a TAAAT site was sufficient to bring an otherwise unresponsive enhancer of the vg gene from Apis under the control of Ubx in a Drosophila transgenic assay. Taken together, our results suggest an evolutionary mechanism by which critical wing patterning genes might have come under the regulation of Ubx in the Dipteran lineage. [ABSTRACT FROM AUTHOR]

Breast cancer is the most prevalent cancer in female patients worldwide. Tissue factor pathway inhibitor 2 (TFPI-2) is identified as an important tumor suppressor in various cancers. Recent studies have shown that TFPI-2 translocates into the nucleus, where it modulates the transcription of the matrix metalloproteinase-2 (MMP-2) gene. However, its biological role and molecular mechanisms in the progression of breast cancer remain unclear. In this study, we identified 5125 differentially expressed genes (DEGs) from RNA sequencing (RNA-seq) in TFPI-2-overexpressing MDA231 cells compared with control cells. Gene ontology and Kyoto encyclopedia of genes and genomes (KEGG) analysis shown that cell cycle, cell differentiation, proteoglycans in cancer, and pathways associated with cancer were highly enriched in downregulated DEGs. Integration of the RNA-seq and ChIP-sequencing (ChIP-seq) data identified 73 genes directly controlled by TFPI-2 in MDA231 cells. Among them, melanocyte inducing transcription factor (MITF) gene expression was repressed by TFPI-2, which was further verified by a luciferase reporter assay and ChIP-quantitative PCR. Our study provides evidence of a novel role of TFPI-2 in human breast cancer involving targeting of the MITF. [ABSTRACT FROM AUTHOR]

Simple Summary: The protein TAZ plays an important role in the Hippo signaling pathway that regulates important biological processes such as cell proliferation, controlled cell death, cancer stem cell traits, tumorigenesis, and resistance to therapies. Specifically, TAZ functions as a transcription coactivator that influences the transcription of various genes. Despite the importance of TAZ, genome-wide occupancy sites for TAZ remain poorly defined. Here, a tetracycline (tet)-inducible gene expression system was used to turn on TAZ activation in mammary epithelial cells, and we characterized genome-wide binding sites at different TAZ activation time points. We found that most TAZ binds to distant enhancer and non-coding regions of genes at earlier TAZ activation time points and binding shifts to nearby promoter regions of genes at later TAZ activation time points. We also found that TAZ activation results in chromatin architecture alterations. These results could lead to the identification of new therapeutic targets for breast cancer. The transcriptional co-activator with PDZ binding motif (TAZ) is a key effector of the Hippo signaling pathway. We and others previously reported that high expression levels of TAZ are positively associated with decreased survival rates and shorter times to relapse in basal-like breast cancer (BLBC) patients. The oncogenic activity of TAZ involves the regulation of diverse signal transduction pathways that direct processes such as cell proliferation, migration, and resistance to apoptosis, albeit through poorly characterized gene expression programs. Here, using a tet-inducible system in mammary epithelial MCF10A cells, we have characterized the TAZ-regulated transcription program using RNA sequencing in a temporal and spatial manner. We further identified global TAZ binding sites at different TAZ activation time points by chromatin immunoprecipitation (ChIP) sequencing analysis. We found that the vast majority of TAZ was rapidly localized in enhancer regions at the early TAZ activation time point and then gradually spread to promoter regions. TAZ bound to enhancer regions following a switch in potential TEAD and FOSL2 transcription factor motifs. Furthermore, the ATAC sequencing analysis indicated that TAZ activation led to chromatin structural alterations. Together, our results have revealed the landscape of genome-wide TAZ binding sites and may lead to improvements in the current understanding of how TAZ regulates the gene expression program that contributes to the development of breast cancer. [ABSTRACT FROM AUTHOR]

Background: Despite well-documented effects on human health, the action modes of environmental pollutants are incompletely understood. Although transcriptome-based approaches are widely used to predict associations between chemicals and disorders, the molecular cues regulating pollutant-derived gene expression changes remain unclear. Therefore, we developed a data-mining approach, termed "DAR-ChIPEA," to identify transcription factors (TFs) playing pivotal roles in the action modes of pollutants. Methods: Large-scale public ChIP-Seq data (human, n = 15,155; mouse, n = 13,156) were used to predict TFs that are enriched in the pollutant-induced differentially accessible genomic regions (DARs) obtained from epigenome analyses (ATAC-Seq). The resultant pollutant–TF matrices were then cross-referenced to a repository of TF–disorder associations to account for pollutant modes of action. We subsequently evaluated the performance of the proposed method using a chemical perturbation data set to compare the outputs of the DAR-ChIPEA and our previously developed differentially expressed gene (DEG)-ChIPEA methods using pollutant-induced DEGs as input. We then adopted the proposed method to predict disease-associated mechanisms triggered by pollutants. Results: The proposed approach outperformed other methods using the area under the receiver operating characteristic curve score. The mean score of the proposed DAR-ChIPEA was significantly higher than that of our previously described DEG-ChIPEA (0.7287 vs. 0.7060; Q = 5.278 × 10–42; two-tailed Wilcoxon rank-sum test). The proposed approach further predicted TF-driven modes of action upon pollutant exposure, indicating that (1) TFs regulating Th1/2 cell homeostasis are integral in the pathophysiology of tributyltin-induced allergic disorders; (2) fine particulates (PM2.5) inhibit the binding of C/EBPs, Rela, and Spi1 to the genome, thereby perturbing normal blood cell differentiation and leading to immune dysfunction; and (3) lead induces fatty liver by disrupting the normal regulation of lipid metabolism by altering hepatic circadian rhythms. Conclusions: Highlighting genome-wide chromatin change upon pollutant exposure to elucidate the epigenetic landscape of pollutant responses outperformed our previously described method that focuses on gene-adjacent domains only. Our approach has the potential to reveal pivotal TFs that mediate deleterious effects of pollutants, thereby facilitating the development of strategies to mitigate damage from environmental pollution. [ABSTRACT FROM AUTHOR]

Background: Cells orchestrate histone biogenesis with strict temporal and quantitative control. To efficiently regulate histone biogenesis, the repetitive Drosophila melanogaster replication-dependent histone genes are arrayed and clustered at a single locus. Regulatory factors concentrate in a nuclear body known as the histone locus body (HLB), which forms around the locus. Historically, HLB factors are largely discovered by chance, and few are known to interact directly with DNA. It is therefore unclear how the histone genes are specifically targeted for unique and coordinated regulation. Results: To expand the list of known HLB factors, we performed a candidate-based screen by mapping 30 publicly available ChIP datasets of 27 unique factors to the Drosophila histone gene array. We identified novel transcription factor candidates, including the Drosophila Hox proteins Ultrabithorax (Ubx), Abdominal-A (Abd-A), and Abdominal-B (Abd-B), suggesting a new pathway for these factors in influencing body plan morphogenesis. Additionally, we identified six other factors that target the histone gene array: JIL-1, hormone-like receptor 78 (Hr78), the long isoform of female sterile homeotic (1) (fs(1)h) as well as the general transcription factors TBP associated factor 1 (TAF-1), Transcription Factor IIB (TFIIB), and Transcription Factor IIF (TFIIF). Conclusions: Our foundational screen provides several candidates for future studies into factors that may influence histone biogenesis. Further, our study emphasizes the powerful reservoir of publicly available datasets, which can be mined as a primary screening technique. [ABSTRACT FROM AUTHOR]

Much of today's molecular science revolves around next-generation sequencing. Frequently, the first step in analyzing such data is aligning sequencing reads to a reference genome. This step is often taken for granted, but any analysis downstream of the alignment will be affected by the aligner's ability to correctly map sequences. In most cases, for research into chromatin structure and nucleosome positioning, ATAC-seq, ChIP-seq, and MNase-seq experiments use short read lengths. How well aligners manage these reads is critical. Most aligner programs will output mapped reads and unmapped reads. However, from a biological point of view, reads will fall into one of three categories: correctly mapped, incorrectly mapped, and unmapped. While increased sequencing depth can often compensate for unmapped reads, incorrectly and correctly mapped reads appear algorithmically identical but can produce biologically significant alterations in the results. For this reason, we are benchmarking various alignment programs to determine their propensity to incorrectly map short reads. As short-read alignment is an important step in ATAC-seq, ChIP-seq, and MNase-seq experiments, caution should be taken in mapping reads to ensure that the most accurate conclusions can be made from the data generated. Our analysis is intended to help investigators new to the field pick the alignment program best suited for their experimental conditions. In general, the aligners we tested performed well. BWA, Bowtie2, and Chromap were all exceptionally accurate, and we recommend using them. Furthermore, we show that longer read lengths do in fact lead to more accurate mappings. [ABSTRACT FROM AUTHOR]

Background: The largest group of patients with breast cancer are estrogen receptor-positive (ER+) type. The estrogen receptor acts as a transcription factor and triggers cell proliferation and differentiation. Hence, investigating ER-DNA interaction genomic regions can help identify genes directly regulated by ER and understand the mechanism of ER action in cancer progression. Methods: In the present study, we employed a workflow to do a meta-analysis of ChIP-seq data of ER+ cell lines stimulated with 10 nM and 100 nM of E2. All publicly available data sets were re-analyzed with the same platform. Then, the known and unknown batch effects were removed. Finally, the meta-analysis was performed to obtain meta-differentially bound sites in estrogen-treated MCF7 cell lines compared to vehicles (as control). Also, the meta-analysis results were compared with the results of T47D cell lines for more precision. Enrichment analyses were also employed to find the functional importance of common meta-differentially bound sites and associated genes among both cell lines. Results: Remarkably, POU5F1B, ZNF662, ZNF442, KIN, ZNF410, and SGSM2 transcription factors were recognized in the meta-analysis but not in individual studies. Enrichment of the meta-differentially bound sites resulted in the candidacy of pathways not previously reported in breast cancer. PCGF2, HNF1B, and ZBED6 transcription factors were also predicted through the enrichment analysis of associated genes. In addition, comparing the meta-analysis results of both ChIP-seq and RNA-seq data showed that many transcription factors affected by ER were up-regulated. Conclusion: The meta-analysis of ChIP-seq data of estrogen-treated MCF7 cell line leads to the identification of new binding sites of ER that have not been previously reported. Also, enrichment of the meta-differentially bound sites and their associated genes revealed new terms and pathways involved in the development of breast cancer which should be examined in future in vitro and in vivo studies. [ABSTRACT FROM AUTHOR]

Ferroptosis is a regulated cell death induced by iron-dependent lipid peroxidation. The heme-responsive transcription factor BTB and CNC homology 1 (BACH1) promotes ferroptosis by repressing the transcription of genes involved in glutathione (GSH) synthesis and intracellular labile iron metabolism, which are key regulatory pathways in ferroptosis. We found that BACH1 re-expression in Bach1 −/− immortalized mouse embryonic fibroblasts (iMEFs) can induce ferroptosis upon 2-mercaptoethanol removal, without any ferroptosis inducers. In these iMEFs, GSH synthesis was reduced, and intracellular labile iron levels were increased upon BACH1 re-expression. We used this system to investigate whether the major ferroptosis regulators glutathione peroxidase 4 (Gpx4) and apoptosis-inducing factor mitochondria-associated 2 (Aifm2), the gene for ferroptosis suppressor protein 1, are target genes of BACH1. Neither Gpx4 nor Aifm2 was regulated by BACH1 in the iMEFs. However, we found that BACH1 represses AIFM2 transcription in human pancreatic cancer cells. These results suggest that the ferroptosis regulators targeted by BACH1 may vary across different cell types and animal species. Furthermore, we confirmed that the ferroptosis induced by BACH1 re-expression exhibited a propagating effect. BACH1 re-expression represents a new strategy for inducing ferroptosis after GPX4 or system Xc− suppression and is expected to contribute to future ferroptosis research. [ABSTRACT FROM AUTHOR]

Summary: Transcriptional corepressors of the Topless (TPL) family regulate plant hormone and immunity signaling. The lack of a genome‐wide profile of their chromatin associations limits understanding of the TPL family roles in transcriptional regulation.Chromatin immunoprecipitation with sequencing (ChIP‐Seq) was performed on Arabidopsis thaliana lines expressing GFP‐tagged Topless‐related 1 (TPR1‐GFP) with and without constitutive immunity via Enhanced Disease Susceptibility 1 (EDS1). RNA‐Seq profiling of the TPR1‐GFP lines and pathogen‐infected tpl/tpr mutants, combined with measuring immunity, growth, and physiological parameters was employed to investigate TPL/TPR roles in immunity and defense homeostasis.TPR1 was enriched at promoter regions of c. 1400 genes and c. 10% of the detected binding required EDS1 immunity signaling. In a tpr1 tpl tpr4 (t3) mutant, resistance to bacteria was slightly compromised, and defense‐related transcriptional reprogramming was weakly reduced or enhanced, respectively, at early (< 1 h) and late 24 h stages of bacterial infection. The t3 plants challenged with bacteria or pathogen‐associated molecular pattern nlp24 displayed photosystem II dysfunctions. Also, t3 plants were hypersensitive to phytocytokine pep1 at the level of root growth inhibition. Transgenic expression of TPR1 rescued these t3 physiological defects.We propose that TPR1 and TPL family proteins function in Arabidopsis to reduce detrimental effects associated with activated transcriptional immunity. [ABSTRACT FROM AUTHOR]

ELT-2 is the major transcription factor (TF) required for Caenorhabditis elegans intestinal development. ELT-2 expression initiates in embryos to promote development and then persists after hatching through the larval and adult stages. Though the sites of ELT-2 binding are characterized and the transcriptional changes that result from ELT-2 depletion are known, an intestine-specific transcriptome profile spanning developmental time has been missing. We generated this dataset by performing Fluorescence Activated Cell Sorting on intestine cells at distinct developmental stages. We analyzed this dataset in conjunction with previously conducted ELT-2 studies to evaluate the role of ELT-2 in directing the intestinal gene regulatory network through development. We found that only 33% of intestineenriched genes in the embryo were direct targets of ELT-2 but that number increased to 75% by the L3 stage. This suggests additional TFs promote intestinal transcription especially in the embryo. Furthermore, only half of ELT-2's direct target genes were dependent on ELT-2 for their proper expression levels, and an equal proportion of those responded to elt-2 depletion with over-expression as with under-expression. That is, ELT-2 can either activate or repress direct target genes. Additionally, we observed that ELT-2 repressed its own promoter, implicating new models for its autoregulation. Together, our results illustrate that ELT-2 impacts roughly 20-50% of intestine-specific genes, that ELT-2 both positively and negatively controls its direct targets, and that the current model of the intestinal regulatory network is incomplete as the factors responsible for directing the expression of many intestinal genes remain unknown. [ABSTRACT FROM AUTHOR]

Background: Genomics data is available to the scientific community after publication of research projects and can be investigated for a multitude of research questions. However, in many cases deposited data is only assessed and used for the initial publication, resulting in valuable resources not being exploited to their full depth. Main: A likely reason for this is that many wetlab-based researchers are not formally trained to apply bioinformatic tools and may therefore assume that they lack the necessary experience to do so themselves. In this article, we present a series of freely available, predominantly web-based platforms and bioinformatic tools that can be combined in analysis pipelines to interrogate different types of next-generation sequencing data. Additionally to the presented exemplary route, we also list a number of alternative tools that can be combined in a mix-and-match fashion. We place special emphasis on tools that can be followed and used correctly without extensive prior knowledge in programming. Such analysis pipelines can be applied to existing data downloaded from the public domain or be compared to the results of own experiments. Conclusion: Integrating transcription factor binding to chromatin (ChIP-seq) with transcriptional output (RNA-seq) and chromatin accessibility (ATAC-seq) can not only assist to form a deeper understanding of the molecular interactions underlying transcriptional regulation but will also help establishing new hypotheses and pre-testing them in silico. [ABSTRACT FROM AUTHOR]

Simple Summary: Skeletal muscle is the most abundant tissue in animals, accounting for 45–60% of the body weight of meat animals, and the growth and development of skeletal muscle has the most direct impact on lean meat percentage and muscle quality. Therefore, exploring the molecular regulation mechanism of skeletal muscle growth and development is an important basis for improving pig meat production performance, and it has also been a research focus in animal genetics in recent years. In this study, ChIP-seq and other methods were used to explore the transcriptional regulation mechanism of new transcription factor MAMSTR in pig skeletal muscle development, The results of this study can further reveal the regulatory mechanism of pig skeletal muscle growth and development, identify new genes related to pig skeletal muscle growth and development, and provide a better reference for pig-farming practices and breeding of high-quality meat pig breeds. Myocyte enhancer factor-2-activating motif and SAP domain-containing transcriptional regulator (MAMSTR) regulates its downstream through binding in its promoter regions. However, its molecular mechanism, particularly the DNA-binding sites, and coregulatory genes are quite unexplored. Therefore, to identify the genome-wide binding sites of the MAMSTR transcription factors and their coregulatory genes, chromatin immunoprecipitation sequencing was carried out. The results showed that MAMSTR was associated with 1506 peaks, which were annotated as 962 different genes. Most of these genes were involved in transcriptional regulation, metabolic pathways, and cell development and differentiation, such as AMPK signaling pathway, TGF-beta signaling pathway, transcription coactivator activity, transcription coactivator binding, adipocytokine signaling pathway, fat digestion and absorption, skeletal muscle fiber development, and skeletal muscle cell differentiation. Lastly, the expression levels and transcriptional activities of PID1, VTI1B, PRKAG1, ACSS2, and SLC28A3 were screened and verified via functional markers and analysis. Overall, this study has increased our understanding of the regulatory mechanism of MAMSTR during skeletal muscle fibroblast development and provided a reference for analyzing muscle development mechanisms. [ABSTRACT FROM AUTHOR]

Hexamerins, the proteins massively stored in the larval haemolymph of insects, are gradually used throughout metamorphosis as a source of raw material and energy for the development of adult tissues. Such behaviour defined hexamerins as storage proteins. Immunofluorescence experiments coupled with confocal microscopy show a hexamerin, HEX 70a, in the nucleus of the brain and fat body cells from honeybee workers, an unexpected localization for a storage protein. HEX 70a colocalizes with fibrillarin, a nucleolar‐specific protein and H3 histone, thus suggesting a potential role as a chromatin‐binding protein. This was investigated through chromatin immunoprecipitation and high‐throughput DNA sequencing (ChIP‐seq). The significant HEX 70a‐DNA binding sites were mainly localized at the intergenic, promoter and intronic regions. HEX 70a targeted DNA stretches mapped to the genomic regions encompassing genes with relevant functional attributes. Several HEX 70a targeted genes were associated with H3K27ac or/and H3K27me3, known as active and repressive histone marks. Brain and fat body tissues shared a fraction of the HEX 70 targeted genes, and tissue‐specific targets were also detected. The presence of overrepresented DNA motifs in the binding sites is consistent with specific HEX 70a‐chromatin association. In addition, a search for HEX 70a targets in RNA‐seq public libraries of fat bodies from nurses and foragers revealed differentially expressed targets displaying hex 70a‐correlated developmental expression, thus supporting a regulatory activity for HEX 70a. Our results support the premise that HEX 70a is a moonlighting protein that binds chromatin and has roles in the brain and fat body cell nuclei, apart from its canonical role as a storage protein. [ABSTRACT FROM AUTHOR]

Background: Changes in gene expression levels during brain development are thought to have played an important role in the evolution of human cognition. With the advent of high-throughput sequencing technologies, changes in brain developmental expression patterns, as well as human-specific brain gene expression, have been characterized. However, interpreting the origin of evolutionarily advanced cognition in human brains requires a deeper understanding of the regulation of gene expression, including the epigenomic context, along the primate genome. Here, we used chromatin immunoprecipitation sequencing (ChIP-seq) to measure the genome-wide profiles of histone H3 lysine 4 trimethylation (H3K4me3) and histone H3 lysine 27 acetylation (H3K27ac), both of which are associated with transcriptional activation in the prefrontal cortex of humans, chimpanzees, and rhesus macaques. Results: We found a discrete functional association, in which H3K4me3HP gain was significantly associated with myelination assembly and signaling transmission, while H3K4me3HP loss played a vital role in synaptic activity. Moreover, H3K27acHP gain was enriched in interneuron and oligodendrocyte markers, and H3K27acHP loss was enriched in CA1 pyramidal neuron markers. Using strand-specific RNA sequencing (ssRNA-seq), we first demonstrated that approximately 7 and 2% of human-specific expressed genes were epigenetically marked by H3K4me3HP and H3K27acHP, respectively, providing robust support for causal involvement of histones in gene expression. We also revealed the co-activation role of epigenetic modification and transcription factors in human-specific transcriptome evolution. Mechanistically, histone-modifying enzymes at least partially contribute to an epigenetic disturbance among primates, especially for the H3K27ac epigenomic marker. In line with this, peaks enriched in the macaque lineage were found to be driven by upregulated acetyl enzymes. Conclusions: Our results comprehensively elucidated a causal species-specific gene-histone-enzyme landscape in the prefrontal cortex and highlighted the regulatory interaction that drove transcriptional activation. [ABSTRACT FROM AUTHOR]

Hypoxia and a suppressive tumour microenvironment (TME) are both independent negative prognostic factors for muscle-invasive bladder cancer (MIBC) that contribute to treatment resistance. Hypoxia has been shown to induce an immune suppressive TME by recruiting myeloid cells that inhibit anti-tumour T cell responses. Recent transcriptomic analyses show hypoxia increases suppressive and anti-tumour immune signalling and infiltrates in bladder cancer. This study sought to investigate the relationship between hypoxia-inducible factor (HIF)-1 and -2, hypoxia, and immune signalling and infiltrates in MIBC. ChIP-seq was performed to identify HIF1α, HIF2α, and HIF1β binding in the genome of the MIBC cell line T24 cultured in 1% and 0.1% oxygen for 24 h. Microarray data from four MIBC cell lines (T24, J82, UMUC3, and HT1376) cultured under 1%, 0.2%, and 0.1% oxygen for 24 h were used. Differences in the immune contexture between high- and low-hypoxia tumours were investigated using in silico analyses of two bladder cancer cohorts (BCON and TCGA) filtered to only include MIBC cases. GO and GSEA were used with the R packages "limma" and "fgsea". Immune deconvolution was performed using ImSig and TIMER algorithms. RStudio was used for all analyses. Under hypoxia, HIF1α and HIF2α bound to ~11.5–13.5% and ~4.5–7.5% of immune-related genes, respectively (1–0.1% O2). HIF1α and HIF2α both bound to genes associated with T cell activation and differentiation signalling pathways. HIF1α and HIF2α had distinct roles in immune-related signalling. HIF1 was associated with interferon production specifically, whilst HIF2 was associated with generic cytokine signalling as well as humoral and toll-like receptor immune responses. Neutrophil and myeloid cell signalling was enriched under hypoxia, alongside hallmark pathways associated with Tregs and macrophages. High-hypoxia MIBC tumours had increased expression of both suppressive and anti-tumour immune gene signatures and were associated with increased immune infiltrates. Overall, hypoxia is associated with increased inflammation for both suppressive and anti-tumour-related immune signalling and immune infiltrates, as seen in vitro and in situ using MIBC patient tumours. [ABSTRACT FROM AUTHOR]

The analysis of super-enhancers (SEs) has recently attracted attention in elucidating the molecular mechanisms of cancer and other diseases. SEs are genomic structures that strongly induce gene expression and have been reported to contribute to the overexpression of oncogenes. Because the analysis of SEs and integrated analysis with other data are performed using large amounts of genome-wide data, artificial intelligence technology, with machine learning at its core, has recently begun to be utilized. In promoting precision medicine, it is important to consider information from SEs in addition to genomic data; therefore, machine learning technology is expected to be introduced appropriately in terms of building a robust analysis platform with a high generalization performance. In this review, we explain the history and principles of SE, and the results of SE analysis using state-of-the-art machine learning and integrated analysis with other data are presented to provide a comprehensive understanding of the current status of SE analysis in the field of medical biology. Additionally, we compared the accuracy between existing machine learning methods on the benchmark dataset and attempted to explore the kind of data preprocessing and integration work needed to make the existing algorithms work on the benchmark dataset. Furthermore, we discuss the issues and future directions of current SE analysis. [ABSTRACT FROM AUTHOR]

The chromosome immunoprecipitation sequencing technology（ChIP-seq）, which combines the chromatin immunoprecipitation technology（ChIP）with the second generation sequencing technology, is an important method for analyzing the epigenetic changes of the whole genome, and can quickly and effectively detect DNA and protein binding sites, transcription factor binding sites（TFBS）, histone post-translation modifications（hPTMs）, nucleosome localization and DNA methylation in the whole genome. However, for a long time, ChIP-seq needs a large number of cells to generate high-quality data sets, which limits the application of ChIP in some specific tissue and low-cell samples such as oocytes, early embryonic cells and other research fields. In recent years, based on ChIP, researchers have proposed a series of low-input ChIP-seq methods to reduce the initial sample amount and experimental cost and increase the sequencing quality, which has promoted the development of epigenomics. In this paper, we reviewed the principles of ChIP and the methodological development of ChIP-seq reducing low-input amount, compared several of the more important methods, and summarized the application of ChIP-seq with lowinput amount in epigenetic research. Finally, we prospected the application and development of low initial ChIP-seq technology, aiming to provide reference for the selection of different low-input ChIP-seq methods for low cell number samples. [ABSTRACT FROM AUTHOR]

In mammals, centromeres are epigenetically specified by the histone H3 variant CENP-A and are typically associated with satellite DNA. We previously described the first example of a natural satellite-free centromere on Equus caballus chromosome 11 (ECA11) and, subsequently, on several chromosomes in other species of the genus Equus. We discovered that these satellite-free neocentromeres arose recently during evolution through centromere repositioning and/or chromosomal fusion, after inactivation of the ancestral centromere, where, in many cases, blocks of satellite sequences were maintained. Here, we investigated by FISH the chromosomal distribution of satellite DNA families in Equus przewalskii (EPR), demonstrating a good degree of conservation of the localization of the major horse satellite families 37cen and 2PI with the domestic horse. Moreover, we demonstrated, by ChIP-seq, that 37cen is the satellite bound by CENP-A and that the centromere of EPR10, the ortholog of ECA11, is devoid of satellite sequences. Our results confirm that these two species are closely related and that the event of centromere repositioning which gave rise to EPR10/ECA11 centromeres occurred in the common ancestor, before the separation of the two horse lineages. [ABSTRACT FROM AUTHOR]

Background: Epigenomic profiling assays such as ChIP-seq have been widely used to map the genome-wide enrichment profiles of chromatin-associated proteins and posttranslational histone modifications. Sequencing depth is a key parameter in experimental design and quality control. However, due to variable sequencing depth requirements across experimental conditions, it can be challenging to determine optimal sequencing depth, particularly for projects involving multiple targets or cell types. Results: We developed the peaksat R package to provide target read depth estimates for epigenomic experiments based on the analysis of peak saturation curves. We applied peaksat to establish the distinctive read depth requirements for ChIP-seq studies of histone modifications in different cell lines. Using peaksat, we were able to estimate the target read depth required per library to obtain high-quality peak calls for downstream analysis. In addition, peaksat was applied to other sequence-enrichment methods including CUT&RUN and ATAC-seq. Conclusion: peaksat addresses a need for researchers to make informed decisions about whether their sequencing data has been generated to an adequate depth and subsequently sufficient meaningful peaks, and failing that, how many more reads would be required per library. peaksat is applicable to other sequence-based methods that include calling peaks in their analysis. [ABSTRACT FROM AUTHOR]

Grainyhead like 2 (GRHL2) is an essential transcription factor for development and function of epithelial tissues. It has dual roles in cancer by supporting tumor growth while suppressing epithelial to mesenchymal transitions (EMT). GRHL2 cooperates with androgen and estrogen receptors (ER) to regulate gene expression. We explore genome wide GRHL2 binding sites conserved in three ER⍺/GRHL2 positive luminal breast cancer cell lines by ChIP-Seq. Interaction with the ER⍺/FOXA1/GATA3 complex is observed, however, only for a minor fraction of conserved GRHL2 peaks. We determine genome wide transcriptional dynamics in response to loss of GRHL2 by nascent RNA Bru-seq using an MCF7 conditional knockout model. Integration of ChIP- and Bru-seq pinpoints candidate direct GRHL2 target genes in luminal breast cancer. Multiple connections between GRHL2 and proliferation are uncovered, including transcriptional activation of ETS and E2F transcription factors. Among EMT-related genes, direct regulation of CLDN4 is corroborated but several targets identified in other cells (including CDH1 and ZEB1) are ruled out by both ChIP- and Bru-seq as being directly controlled by GRHL2 in luminal breast cancer cells. Gene clusters correlating positively (including known GRHL2 targets such as ErbB3, CLDN4/7) or negatively (including TGFB1 and TGFBR2) with GRHL2 in the MCF7 knockout model, display similar correlation with GRHL2 in ER positive as well as ER negative breast cancer patients. Altogether, this study uncovers gene sets regulated directly or indirectly by GRHL2 in luminal breast cancer, identifies novel GRHL2-regulated genes, and points to distinct GRHL2 regulation of EMT in luminal breast cancer cells. 6_BaCm7NaLGTK4K5qg2XWL Video Abstract [ABSTRACT FROM AUTHOR]

Background: Elucidating the Transcription Factors (TFs) that drive the gene expression changes in a given experiment is a common question asked by researchers. The existing methods rely on the predicted Transcription Factor Binding Site (TFBS) to model the changes in the motif activity. Such methods only work for TFs that have a motif and assume the TF binding profile is the same in all cell types. Results: Given the wealth of the ChIP-seq data available for a wide range of the TFs in various cell types, we propose that gene expression modeling can be done using ChIP-seq "signatures" directly, effectively skipping the motif finding and TFBS prediction steps. We present xcore, an R package that allows TF activity modeling based on ChIP-seq signatures and the user's gene expression data. We also provide xcoredata a companion data package that provides a collection of preprocessed ChIP-seq signatures. We demonstrate that xcore leads to biologically relevant predictions using transforming growth factor beta induced epithelial-mesenchymal transition time-courses, rinderpest infection time-courses, and embryonic stem cells differentiated to cardiomyocytes time-course profiled with Cap Analysis Gene Expression. Conclusions: xcore provides a simple analytical framework for gene expression modeling using linear models that can be easily incorporated into differential expression analysis pipelines. Taking advantage of public ChIP-seq databases, xcore can identify meaningful molecular signatures and relevant ChIP-seq experiments. [ABSTRACT FROM AUTHOR]

Simple Summary: The sea cucumber Apostichopus japonicus is an ecologically significant echinoderm with considerable nutritional and therapeutic value. Color variation of A. japonicus is one of its most important cultivation traits, as color affects its bioactive compounds, taste, and market price. Microphthalmia-associated transcription factor (MITF) is one of the most critical genes in melanocyte development and melanin synthesis pathways. However, how MITF regulates body color and differentiation in sea cucumbers is poorly understood. Here, we analyzed the expression level and location of MITF in white, purple, and green sea cucumbers and identified the genes regulated by MITF using chromatin immunoprecipitation followed by sequencing. Our findings suggested that MITF contributed to the body color differentiation of green, purple, and white sea cucumbers A. japonicus through expression differences and regulation of downstream genes. These findings give researchers a foundation for further research into the processes that underlie the development of body color in sea cucumbers and provide light on how genes are regulated. The results will also provide support for sea cucumber breeding and a useful theoretical basis in order to contribute to the development of the sea cucumber industry. Melanin, which is a pigment produced in melanocytes, is an important contributor to sea cucumber body color. MITF is one of the most critical genes in melanocyte development and melanin synthesis pathways. However, how MITF regulates body color and differentiation in sea cucumbers is poorly understood. In this study, we analyzed the expression level and location of MITF in white, purple, and green sea cucumbers and identified the genes regulated by MITF using chromatin immunoprecipitation followed by sequencing. The mRNA and protein expression levels of MITF were all highest in purple morphs and lowest in white morphs. In situ hybridization indicated that MITF mRNA were mainly expressed in the epidermis. We also identified 984, 732, and 1191 peaks of MITF binding in green, purple, and white sea cucumbers, which were associated with 727, 557, and 887 genes, respectively. Our findings suggested that MITF contributed to the body color differentiation of green, purple, and white sea cucumbers through expression differences and regulation of downstream genes. These results provided a basis for future studies to determine the mechanisms underlying body color formation and provided insights into gene regulation in sea cucumbers. [ABSTRACT FROM AUTHOR]

The sensitivity of human immune system cells to gravity changes has been investigated in numerous studies. Human macrophages mediate innate and thus rapid immune defense on the one hand and activate T- and B-cell-based adaptive immune response on the other hand. In this process they finally act as immunoeffector cells, and are essential for tissue regeneration and remodeling. Recently, we demonstrated in the human Jurkat T cell line that genes are differentially regulated in cluster structures under altered gravity. In order to study an in vivo near system of immunologically relevant human cells under physically real microgravity, we performed parabolic flight experiments with primary human M1 macrophages under highly standardized conditions and performed chromatin immunoprecipitation DNA sequencing (ChIP-Seq) for whole-genome epigenetic detection of the DNA-binding loci of the main transcription complex RNA polymerase II and the transcription-associated epigenetic chromatin modification H3K4me3. We identified an overall downregulation of H3K4me3 binding loci in altered gravity, which were unequally distributed inter- and intrachromosomally throughout the genome. Three-quarters of all affected loci were located on the p arm of the chromosomes chr5, chr6, chr9, and chr19. The genomic distribution of the downregulated H3K4me3 loci corresponds to a substantial extent to immunoregulatory genes. In microgravity, analysis of RNA polymerase II binding showed increased binding to multiple loci at coding sequences but decreased binding to central noncoding regions. Detection of altered DNA binding of RNA polymerase II provided direct evidence that gravity changes can lead to altered transcription. Based on this study, we hypothesize that the rapid transcriptional response to changing gravitational forces is specifically encoded in the epigenetic organization of chromatin. [ABSTRACT FROM AUTHOR]

Emerging evidence illustrates that RhoC has divergent roles in cervical cancer progression where it controls epithelial to mesenchymal transition (EMT), migration, angiogenesis, invasion, tumor growth, and radiation response. Cancer stem cells (CSCs) are the primary cause of recurrence and metastasis and exhibit all of the above phenotypes. It, therefore, becomes imperative to understand if RhoC regulates CSCs in cervical cancer. In this study, cell lines and clinical specimen-based findings demonstrate that RhoC regulates tumor phenotypes such as clonogenicity and anoikis resistance. Accordingly, inhibition of RhoC abrogated these phenotypes. RNA-seq analysis revealed that RhoC over-expression resulted in up-regulation of 27% of the transcriptome. Further, the Infinium MethylationEPIC array showed that RhoC over-expressing cells had a demethylated genome. Studies divulged that RhoC via TET2 signaling regulated the demethylation of the genome. Further investigations comprising ChIP-seq, reporter assays, and mass spectrometry revealed that RhoC associates with WDR5 in the nucleus and regulates the expression of pluripotency genes such as Nanog. Interestingly, clinical specimen-based investigations revealed the existence of a subset of tumor cells marked by RhoC+/Nanog+ expression. Finally, combinatorial inhibition (in vitro) of RhoC and its partners (WDR5 and TET2) resulted in increased sensitization of clinical specimen-derived cells to radiation. These findings collectively reveal a novel role for nuclear RhoC in the epigenetic regulation of Nanog and identify RhoC as a regulator of CSCs. The study nominates RhoC and associated signaling pathways as therapeutic targets. [ABSTRACT FROM AUTHOR]

Objective: To assess the genetic and epigenetic effects promoted by Bisphenol A (BPA) exposure in adolescent males from the Spanish INMA-Granada birth cohort, and in human cells. Methods: DNA methylation was analysed using MEDIP. Repeat number variation in genomic DNA was evaluated, along with the analysis of H3K4me3 by using chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq). Analyses were performed with material extracted from whole blood of the adolescents, complemented by in vitro assessments of human (HeLa) cells exposed to 10 nM BPA, specifically, immunofluorescence evaluation of protein levels, gene expression analysis and ChIP‒qPCR analysis. Results: Adolescents in the high urinary BPA levels group presented a higher level of Satellite A (SATA) repetitive region copy numbers compared to those in the low BPA group and a tendency towards increase in telomere length. We also observed decreased DNA methylation at the promoters of the imprinted genes H19, KCNQ1, and IGF2; at LINE1 retroelements; and at the ARID2, EGFR and ESRRA and TERT genes. Genome-wide sequencing revealed increased H3K4me3 occupancy at the promoters of genes encoding histone acetyltransferases, telomeric DNA binding factors and DNA repair genes. Results were supported in HeLa cells exposed to 10 nM BPA in vitro. In accordance with the data obtained in blood samples, we observed higher H3K4me3 occupancy and lower DNA methylation at some specific targets in HeLa cells. In exposed cells, changes in the expression of genes encoding DNA repair factors (ATM, ARID2, TRP53) were observed, and increased expression of several genes encoding telomeric DNA binding factors (SMG7, TERT, TEN1, UPF1, ZBTB48) were also found. Furthermore, an increase in ESR1/ERa was observed in the nuclei of HeLa cells along with increased binding of ESR1 to KAT5, KMT2E and TERF2IP promoters and decreased ESR1 binding at the RARA promoter. The DNA damage marker p53/TP53 was also increased. Conclusion: In this pilot study, genome-wide analysis of histone trimethylation in adolescent males exposed to BPA revealed a global impact on the expression of genes encoding telomeric binding proteins and histone acetyltransferase factors with similar results in HeLa cells. Nevertheless, larger studies should confirm our findings. [ABSTRACT FROM AUTHOR]

Background: The cis-regulatory element became increasingly important for resistance breeding. There were many DNA variations identified by resequencing. To investigate the links between the DNA variations and cis-regulatory element was the fundamental work. DNA variations in cis-regulatory elements caused phenotype variations in general. Results: We used WGBS, ChIP-seq and RNA-seq technology to decipher the regulatory element landscape from eight hulless barley varieties under four kinds of abiotic stresses. We discovered 231,440 lowly methylated regions (LMRs) from the methylome data of eight varieties. The LMRs mainly distributed in the intergenic regions. A total of 97,909 enhancer-gene pairs were identified from the correlation analysis between methylation degree and expression level. A lot of enriched motifs were recognized from the tolerant-specific LMRs. The key transcription factors were screened out and the transcription factor regulatory network was inferred from the enhancer-gene pairs data for drought stress. The NAC transcription factor was predicted to target to TCP, bHLH, bZIP transcription factor genes. We concluded that the H3K27me3 modification regions overlapped with the LMRs more than the H3K4me3. The variation of single nucleotide polymorphism was more abundant in LMRs than the remain regions of the genome. Conclusions: Epigenetic regulation is an important mechanism for organisms to adapt to complex environments. Through the study of DNA methylation and histone modification, we found that many changes had taken place in enhancers and transcription factors in the abiotic stress of hulless barley. For example, transcription factors including NAC may play an important role. This enriched the molecular basis of highland barley stress response. [ABSTRACT FROM AUTHOR]

Background: Causal variants for complex traits, such as eQTL are often found in non-coding regions of the genome, where they are hypothesised to influence phenotypes by regulating gene expression. Many regulatory regions are marked by histone modifications, which can be assayed by chromatin immunoprecipitation followed by sequencing (ChIP-seq). Sequence reads from ChIP-seq form peaks at putative regulatory regions, which may reflect the amount of regulatory activity at this region. Therefore, eQTL which are also associated with differences in histone modifications are excellent candidate causal variants. Results: We assayed the histone modifications H3K4Me3, H3K4Me1 and H3K27ac and mRNA in the mammary gland of up to 400 animals. We identified QTL for peak height (histone QTL), exon expression (eeQTL), allele specific expression (aseQTL) and allele specific binding (asbQTL). By intersecting these results, we identify variants which may influence gene expression by altering regulatory regions of the genome, and may be causal variants for other traits. Lastly, we find that these variants are found in putative transcription factor binding sites, identifying a mechanism for the effect of many eQTL. Conclusions: We find that allele specific and traditional QTL analysis often identify the same genetic variants and provide evidence that many eQTL are regulatory variants which alter activity at regulatory regions of the bovine genome. Our work provides methodological and biological updates on how regulatory mechanisms interplay at multi-omics levels. [ABSTRACT FROM AUTHOR]

Background: Despite the advent of Chromatin Immunoprecipitation Sequencing (ChIP-seq) having revolutionised our understanding of the mammalian genome's regulatory landscape, many challenges remain. In particular, because of their repetitive nature, the sequencing reads derived from transposable elements (TEs) pose a real bioinformatics challenge, to the point that standard analysis pipelines typically ignore reads whose genomic origin cannot be unambiguously ascertained. Results: We show that discarding ambiguously mapping reads may lead to a systematic underestimation of the number of reads associated with young TE families/subfamilies. We also provide evidence suggesting that the strategy of randomly permuting the location of the read mappings (or the TEs) that is often used to compute the background for enrichment calculations at TE families/subfamilies can result in both false positive and negative enrichments. To address these problems, we present the Transposable Element Enrichment Estimator (T3E), a tool that makes use of ChIP-seq data to characterise the epigenetic profile of associated TE families/subfamilies. T3E weights the number of read mappings assigned to the individual TE copies of a family/subfamily by the overall number of genomic loci to which the corresponding reads map, and this is done at the single nucleotide level. In addition, T3E computes ChIP-seq enrichment relative to a background estimated based on the distribution of the read mappings in the input control DNA. We demonstrated the capabilities of T3E on 23 different ChIP-seq libraries. T3E identified enrichments that were consistent with previous studies. Furthermore, T3E detected context-specific enrichments that are likely to pinpoint unexplored TE families/subfamilies with individual TE copies that have been frequently exapted as cis-regulatory elements during the evolution of mammalian regulatory networks. Conclusions: T3E is a novel open-source computational tool (available for use at: https://github.com/michelleapaz/T3E) that overcomes some of the pitfalls associated with the analysis of ChIP-seq data arising from the repetitive mammalian genome and provides a framework to shed light on the epigenetics of entire TE families/subfamilies. [ABSTRACT FROM AUTHOR]

Knowledge of the transcription factor binding landscape (TFBL) is necessary to analyze gene regulatory networks for important agronomic traits. However, a low‐cost and high‐throughput in vivo chromatin profiling method is still lacking in plants. Here, we developed a transient and simplified cleavage under targets and tagmentation (tsCUT&Tag) that combines transient expression of transcription factor proteins in protoplasts with a simplified CUT&Tag without nucleus extraction. Our tsCUT&Tag method provided higher data quality and signal resolution with lower sequencing depth compared with traditional ChIP‐seq. Furthermore, we developed a strategy combining tsCUT&Tag with machine learning, which has great potential for profiling the TFBL across plant development. [ABSTRACT FROM AUTHOR]