Protocol for the fast chromatin immunoprecipitation (ChIP) method

Integrating Arabidopsis and crop species gene discovery for crop improvement

Genome sequence assemblies form a durable and precise framework supporting nearly all areas of biological research, including evolutionary biology, taxonomy and conservation science, pathogen population diversity, crop domestication, and biochemistry. In the early days of plant genomics, resources were limited to a handful of tractable genomes, leading to a tension between focus on discovering mechanisms in experimental species such as Arabidopsis thaliana (Arabidopsis) and on trait analyses in crop species. This tension arose from challenges in translating knowledge of gene function across the large evolutionary distances between Arabidopsis and diverse crop species in the absence of comparative genome support. For some time, these clashing interests influenced funding priorities in plant science that limited both the acquisition of knowledge of mechanisms in Arabidopsis and the timely development of the capacity of crop science to incorporate emerging knowledge of genes and their mechanisms. In this review we show how advances in genomics analysis technologies are revealing a high degree of conservation of molecular mechanisms between evolutionarily distant plant species. This progress is bridging the model-species-to-crop barrier, resulting in ever-increasing unification of plant science that is now accelerating progress in understanding mechanisms underlying diverse traits in crops and improving their performance. We lay out some examples of important priorities and outcomes arising from these new opportunities.

Transcriptome analysis reveals the mechanism of tolerance to copper toxicity in the white rot fungus Trametes hirsuta AH28-2

Heavy metals, such as copper (Cu), are prevalent in the environment and pose a substantial threat to human health. White rot fungi, especially Trametes spp., display prominent Cu tolerance and removal capacity. However, how Trametes responds to environmental Cu stress remains poorly understood. Here, we found that Trametes hirsuta AH28-2 exhibits Cu removal efficiencies varying from 80.8 % at 1.25 mg/L to 57.6 % at 37 mg/L. Comparative transcriptome analysis identified 812, 1898, and 2110 differentially expressed genes (DEGs) at the Cu concentrations of 1.25, 12.5, and 25 mg/L, respectively. Some DEGs were associated with antioxidant defense systems, secondary metabolite biosynthesis (terpenoids and polyketides), transmembrane transport, and glutathione metabolism, potentially enhancing Cu tolerance. The activities of antioxidant enzymes such as superoxide dismutase, catalase, and laccase were increased under Cu stress. qRT-PCR confirmed the alterations in gene expression and demonstrated that glutathione S-transferases, catalases, cytochrome P450s, and laccases were involved in counteracting Cu-induced stress. Gene silencing experiments further confirmed the crucial roles of laccases in this process. Many transcription factors were enriched under Cu stress, including the Zn2Cys6 family transcription factor GME8421_g (TH8421), which was significantly upregulated at the Cu concentration of 12.5 mg/L. ChIP-seq identified five antioxidant enzyme-encoding genes as direct targets of TH8421, forming a regulatory network that protects against Cu stress. These findings offer insights into the molecular mechanisms driving Cu toxicity tolerance in Trametes fungi.

Antiepileptic drugs carbamazepine and valproic acid mediate transcriptional activation of CYP1A1 via aryl hydrocarbon receptor and regulation of estrogen metabolism

Cytochrome P450 1A1 (CYP1A1) actively catalyzes estrogen hydroxylation reactions and maintains the levels of neuroactive steroid estradiol. The widely prescribed first-line anti-epileptic drugs (AEDs) are considered to be a potent inducer of CYP1A1 and have also been observed to affect serum estradiol and calcium levels in patients with epilepsy. However, the ability of AEDs to interfere with CYP enzyme function and estrogen disposition is a relatively unexplored area. Here we investigate the effect of widely prescribed AEDs (carbamazepine and valproic acid) on CYP1A1 regulation and the levels of estradiol and calcium in cell supernatants of hepatocellular, HepG2, and neuronal, SH-SY5Y cells. We observed that both the AEDs significantly increased CYP1A1 expression and enzyme activity, which was accompanied by a decrease in estradiol and calcium levels in HepG2 cells. This induction of CYP1A1 mRNA and protein was fully prevented by aryl hydrocarbon receptor (AHR) knockdown and StemRegenin 1 (SR1) antagonism. Notably, the AEDs did not affect the AHR expression but regulated its nuclear translocation, potentially driving the transcriptional upregulation of CYP1A1. Furthermore, the knockdown of CYP1A1 in HepG2 cells elucidated a marked increase in estradiol and calcium levels. Later, this increase subsided upon AED exposure. Lastly, we observed a similar trend in estradiol and calcium alterations in SH-SY5Y cells on AED exposure, speculating the involvement of CYP1A1 induction via AEDs at neuronal sites. This work demonstrates that AEDs mediate the upregulation of CYP1A1 via an AHR-dependent mechanism and influence estrogen and calcium homeostasis.

Protocol for in vivo chromatin immunoprecipitation on purified chromatin isolated from mouse liver nuclei

Chromatin immunoprecipitation (ChIP) is used to investigate genome binding by transcription factors, but it can be problematic. We present a protocol to isolate fixed DNA-protein complexes from mouse liver prior to chromatin shearing. We describe steps for liver disaggregation and cross-linking, DNA-protein complex isolation, chromatin shearing, and quality control analysis as well as procedures for ChIP, DNA purification, and ChIP analysis. This protocol yields high-quality samples using commercial antibodies. For complete details on the use and execution of this protocol, please refer to Akl et al.1.

BRF110, an Orally Active Nurr1-RXRα-Selective Rexinoid, Enhances BDNF Expression without Elevating Triglycerides

We report the discovery of a Nurr1-RXRα heterodimer-selective rexinoid which emerged from the structural modification of aminopyrimidine XCT0135908. Although XCT0135908 demonstrated high selectivity for the Nurr1-RXRα heterodimer over other RXRα dimerization partners, its poor in vivo stability and limited brain penetration hindered its utility. Structure-activity relationship (SAR) studies alongside bioactivity evaluations of a diverse series of substituted pyrimidines led to BRF110, a brain-penetrant compound retaining the selective activation of the Nurr1-RXRα heterodimer. BRF110, as XCT0135908, protects dopaminergic cells against the Parkinson's disease-related toxin MPP+ and increases BDNF transcription in mice. Notably, BRF110, in contrast to the market-approved pan-RXR agonist bexarotene, did not elevate triglyceride levels, indicating that enhanced heterodimer selectivity can mitigate off-target in vivo side effects of rexinoids. These findings highlight the potential of heterodimer-selective scaffolds as a strategy for improving the therapeutic profile of rexinoids, addressing significant challenges in the clinical development of RXR-targeting molecules.

H3K56 acetylation regulates chromatin maturation following DNA replication

Following DNA replication, the newly reassembled chromatin is disorganized and must mature to its steady state to maintain both genome and epigenome integrity. However, the regulatory mechanisms governing this critical process remain poorly understood. Here, we show that histone H3K56 acetylation (H3K56ac), a mark on newly-synthesized H3, facilitates the remodeling of disorganized nucleosomes in nascent chromatin, and its removal at the subsequent G2/M phase of the cell cycle marks the completion of chromatin maturation. In vitro, H3K56ac enhances the activity of ISWI chromatin remodelers, including yeast ISW1 and its human equivalent SNF2h. In vivo, a deficiency of H3K56ac in nascent chromatin results in the formation of closely packed di-nucleosomes and/or tetra-nucleosomes. In contrast, abnormally high H3K56ac levels disrupt chromatin maturation, leading to genome instability. These findings establish a central role of H3K56ac in chromatin maturation and reveal a mechanism regulating this critical aspect of chromosome replication.

Analysis of DNA Methylation in Gliomas: Assessment of Preanalytical Variables

Precision oncology is driven by biomarkers. For glioblastoma multiforme (GBM), the most common malignant adult primary brain tumor, O6-methylguanine-DNA methyltransferase (MGMT) gene promoter methylation is an important prognostic and treatment clinical biomarker. Time-consuming preanalytical steps such as biospecimen storage, fixation, sampling, and processing are sources of data irreproducibility, and all these preanalytical variables are confounded by intratumor heterogeneity of MGMT promoter methylation. To assess the effect of preanalytical variables on GBM DNA methylation, tissue storage/sampling (CryoGrid), sample preparation multisonicator (PIXUL), and 5-methylcytosine DNA immunoprecipitation (Matrix-MeDIP-qPCR/seq) platforms were used. MGMT promoter methylation status assayed by MeDIP-qPCR was validated with methylation-specific polymerase chain reaction. MGMT promoter methylation levels in frozen and formalin-fixed paraffin-embedded sample pairs were not statistically different, confirming the reliability of formalin-fixed paraffin-embedded for MGMT promoter methylation analysis. Warm ex vivo ischemia (up to 4 hours at 37 °C) and 3 cycles of repeated sample thawing and freezing did not statistically impact 5-methylcytosine at MGMT promoter, exon, and enhancer regions, indicating the resistance of DNA methylation to common variations in sample processing conditions that might be encountered in research and clinical settings. Twenty-six percent to 34% of specimens exhibited intratumor heterogeneity in the MGMT DNA promoter methylation. These data demonstrate that variations in sample fixation, ischemia duration and temperature, and DNA methylation assay technique do not have a statistically significant impact on MGMT promoter methylation assessment. However, intratumor methylation heterogeneity underscores the value of multiple biopsies at different GBM geographic tumor sites in the evaluation of MGMT promoter methylation status. Matrix-MeDIP-seq analysis revealed that MGMT promoter methylation status clustered with other differentially methylated genomic loci (eg, HOXA and lncRNAs) that are resilient to variation in the above preanalytical conditions. These observations offer new opportunities to develop more granular data-based epigenetic GBM biomarkers. In this regard, the high-throughput CryoGrid-PIXUL-Matrix toolbox could be useful.

NKX3-2 Induces Ovarian Cancer Cell Migration by HDAC6-Mediated Repositioning of Lysosomes and Inhibition of Autophagy

Several soluble factors secreted by the stromal cells and cancer cells within the tumor microenvironment facilitate the progression and invasiveness of ovarian cancer. In ovarian cancer cells, lysophosphatidic acid (LPA) modulates the transcriptome profile and promotes cell invasiveness by the downregulation of autophagy. Here, we further elucidate this mechanism by focusing on the molecular and cellular events regulating autophagy. Transcriptomic and Western blotting analyses revealed NKX3-2, a transcriptional factor, to be among the genes hyperexpressed in LPA-stimulated ovarian cancer cells. Bioinformatic analyses revealed that in ovarian cancer patients, the expression of NKX3-2 positively correlates with genes involved in cell motility and migration, while it negatively correlates with macromolecular catabolic pathways. In various ovarian cancer cell lines, NKX3-2 silencing abrogated LPA-induced cell migration. Mechanistically, this effect is linked to the restoration of the HDAC6-mediated relocation of the lysosomes in the para-golgian area, and this results in an increase in autolysosome formation and the overall upregulation of autophagy. Silencing the expression of ATG7 or BECN1, two autophagy genes, rescued the migratory phenotype of the NKX3-2-silenced ovarian cancer cells. Taken together, these data reveal the mechanism by which the LPA-NKX3-2 axis promotes the invasiveness of ovarian cancer cells and supports the possibility of targeting NKX3-2 to reduce the migratory capacity of cancer cells in response to a permissive microenvironment.

PDNAPred: Interpretable prediction of protein-DNA binding sites based on pre-trained protein language models

Protein-DNA interactions play critical roles in various biological processes and are essential for drug discovery. However, traditional experimental methods are labor-intensive and unable to keep pace with the increasing volume of protein sequences, leading to a substantial number of proteins lacking DNA-binding annotations. Therefore, developing an efficient computational method to identify protein-DNA binding sites is crucial. Unfortunately, most existing computational methods rely on manually selected features or protein structure information, making these methods inapplicable to large-scale prediction tasks. In this study, we introduced PDNAPred, a sequence-based method that combines two pre-trained protein language models with a designed CNN-GRU network to identify DNA-binding sites. Additionally, to tackle the issue of imbalanced dataset samples, we employed focal loss. Our comprehensive experiments demonstrated that PDNAPred significantly improved the accuracy of DNA-binding site prediction, outperforming existing state-of-the-art sequence-based methods. Remarkably, PDNAPred also achieved results comparable to advanced structure-based methods. The designed CNN-GRU network enhances its capability to detect DNA-binding sites accurately. Furthermore, we validated the versatility of PDNAPred by training it on RNA-binding site datasets, showing its potential as a general framework for amino acid binding site prediction. Finally, we conducted model interpretability analysis to elucidate the reasons behind PDNAPred's outstanding performance.

Multi-scale DNA language model improves 6 mA binding sites prediction

DNA methylation at the N6 position of adenine (N6-methyladenine, 6 mA), which refers to the attachment of a methyl group to the N6 site of the adenine (A) of DNA, is an important epigenetic modification in prokaryotic and eukaryotic genomes. Accurately predicting the 6 mA binding sites can provide crucial insights into gene regulation, DNA repair, disease development and so on. Wet experiments are commonly used for analyzing 6 mA binding sites. However, they suffer from high cost and expensive time. Therefore, various deep learning methods have been widely used to predict 6 mA binding sites recently. In this study, we develop a framework based on multi-scale DNA language model named "iDNA6mA-MDL". "iDNA6mA-MDL" integrates multiple kmers and the nucleotide property and frequency method for feature embedding, which can capture a full range of DNA sequence context information. At the prediction stage, it also leverages DNABERT to compensate for the incomplete capture of global DNA information. Experiments show that our framework obtains average AUC of 0.981 on a classic 6 mA rice gene dataset, going beyond all existing advanced models under fivefold cross-validations. Moreover, "iDNA6mA-MDL" outperforms most of the popular state-of-the-art methods on another 11 6 mA datasets, demonstrating its effectiveness in 6 mA binding sites prediction.

Inhibition of heme-thiolate monooxygenase CYP1B1 prevents hepatic stellate cell activation and liver fibrosis by accumulating trehalose

Activation of extracellular matrix-producing hepatic stellate cells (HSCs) is a key event in liver fibrogenesis. We showed that the expression of the heme-thiolate monooxygenase cytochrome P450 1B1 (CYP1B1) was elevated in human and mouse fibrotic livers and activated HSCs. Systemic or HSC-specific ablation and pharmacological inhibition of CYP1B1 attenuated HSC activation and protected male but not female mice from thioacetamide (TAA)-, carbon tetrachloride (CCl4)-, or bile duct ligation (BDL)-induced liver fibrosis. Metabolomic analysis revealed an increase in the disaccharide trehalose in CYP1B1-deficient HSCs resulting from intestinal suppression of the trehalose-metabolizing enzyme trehalase, whose gene we found to be a target of RARα. Trehalose or its hydrolysis-resistant derivative lactotrehalose exhibited potent antifibrotic activity in vitro and in vivo by functioning as an HSC-specific autophagy inhibitor, which may account for the antifibrotic effect of CYP1B1 inhibition. Our study thus reveals an endobiotic function of CYP1B1 in liver fibrosis in males, mediated by liver-intestine cross-talk and trehalose. At the translational level, pharmacological inhibition of CYP1B1 or the use of trehalose/lactotrehalose may represent therapeutic strategies for liver fibrosis.

The role of hexokinases in epigenetic regulation: altered hexokinase expression and chromatin stability in yeast

BACKGROUND

Human hexokinase 2 (HK2) plays an important role in regulating Warburg effect, which metabolizes glucose to lactate acid even in the presence of ample oxygen and provides intermediate metabolites to support cancer cell proliferation and tumor growth. HK2 overexpression has been observed in various types of cancers and targeting HK2-driven Warburg effect has been suggested as a potential cancer therapeutic strategy. Given that epigenetic enzymes utilize metabolic intermediates as substrates or co-factors to carry out post-translational modification of histones and nucleic acids modifications in cells, we hypothesized that altering HK2 expression could impact the epigenome and, consequently, chromatin stability in yeast. To test this hypothesis, we established genetic models with different yeast hexokinase 2 (HXK2) expression in Saccharomyces cerevisiae yeast cells and investigated the effect of HXK2-dependent metabolism on parental nucleosome transfer, a key DNA replication-coupled epigenetic inheritance process, and chromatin stability.

RESULTS

By comparing the growth of mutant yeast cells carrying single deletion of hxk1Δ, hxk2Δ, or double-loss of hxk1Δ hxk2Δ to wild-type cells, we firstly confirmed that HXK2 is the dominant HXK in yeast cell growth. Surprisingly, manipulating HXK2 expression in yeast, whether through overexpression or deletion, had only a marginal impact on parental nucleosome assembly, but a noticeable trend with decrease chromatin instability. However, targeting yeast cells with 2-deoxy-D-glucose (2-DG), a clinical glycolysis inhibitor that has been proposed as an anti-cancer treatment, significantly increased chromatin instability.

CONCLUSION

Our findings suggest that in yeast cells lacking HXK2, alternative HXKs such as HXK1 or glucokinase 1 (GLK1) play a role in supporting glycolysis at a level that adequately maintains epigenomic stability. While our study demonstrated an increase in epigenetic instability with 2-DG treatment, the observed effect seemed to occur dependent on non-glycolytic function of Hxk2. Thus, additional research is needed to identify the molecular mechanism through which 2-DG influences chromatin stability.

A non-canonical role for a small nucleolar RNA in ribosome biogenesis and senescence

Cellular senescence is an irreversible state of cell-cycle arrest induced by various stresses, including aberrant oncogene activation, telomere shortening, and DNA damage. Through a genome-wide screen, we discovered a conserved small nucleolar RNA (snoRNA), SNORA13, that is required for multiple forms of senescence in human cells and mice. Although SNORA13 guides the pseudouridylation of a conserved nucleotide in the ribosomal decoding center, loss of this snoRNA minimally impacts translation. Instead, we found that SNORA13 negatively regulates ribosome biogenesis. Senescence-inducing stress perturbs ribosome biogenesis, resulting in the accumulation of free ribosomal proteins (RPs) that trigger p53 activation. SNORA13 interacts directly with RPL23, decreasing its incorporation into maturing 60S subunits and, consequently, increasing the pool of free RPs, thereby promoting p53-mediated senescence. Thus, SNORA13 regulates ribosome biogenesis and the p53 pathway through a non-canonical mechanism distinct from its role in guiding RNA modification. These findings expand our understanding of snoRNA functions and their roles in cellular signaling.

The cellular Notch1 protein promotes KSHV reactivation in an Rta-dependent manner

The cellular Notch signal transduction pathway is intimately associated with infections by Kaposi's sarcoma-associated herpesvirus (KSHV) and other gamma-herpesviruses. RBP-Jk, the cellular DNA binding component of the canonical Notch pathway, is the key Notch downstream effector protein in virus-infected and uninfected animal cells. Reactivation of KSHV from latency requires the viral lytic switch protein, Rta, to form complexes with RBP-Jk on numerous sites within the viral DNA. Constitutive Notch activity is essential for KSHV pathophysiology in models of Kaposi's sarcoma (KS) and Primary Effusion Lymphoma (PEL), and we demonstrate that Notch1 is also constitutively active in infected Vero cells. Although the KSHV genome contains >100 RBP-Jk DNA motifs, we show that none of the four isoforms of activated Notch can productively reactivate the virus from latency in a highly quantitative trans-complementing reporter virus system. Nevertheless, Notch contributed positively to reactivation because broad inhibition of Notch1-4 with gamma-secretase inhibitor (GSI) or expression of dominant negative mastermind-like1 (dnMAML1) coactivators severely reduced production of infectious KSHV from Vero cells. Reduction of KSHV production is associated with gene-specific reduction of viral transcription in both Vero and PEL cells. Specific inhibition of Notch1 by siRNA partially reduces the production of infectious KSHV, and NICD1 forms promoter-specific complexes with viral DNA during reactivation. We conclude that constitutive Notch activity is required for the robust production of infectious KSHV, and our results implicate activated Notch1 as a pro-viral member of a MAML1/RBP-Jk/DNA complex during viral reactivation.

IMPORTANCE

Kaposi's sarcoma-associated herpesvirus (KSHV) manipulates the host cell oncogenic Notch signaling pathway for viral reactivation from latency and cell pathogenesis. KSHV reactivation requires that the viral protein Rta functionally interacts with RBP-Jk, the DNA-binding component of the Notch pathway, and with promoter DNA to drive transcription of productive cycle genes. We show that the Notch pathway is constitutively active during KSHV reactivation and is essential for robust production of infectious virus progeny. Inhibiting Notch during reactivation reduces the expression of specific viral genes yet does not affect the growth of the host cells. Although Notch cannot reactivate KSHV alone, the requisite expression of Rta reveals a previously unappreciated role for Notch in reactivation. We propose that activated Notch cooperates with Rta in a promoter-specific manner that is partially programmed by Rta's ability to redistribute RBP-Jk DNA binding to the virus during reactivation.

Examining NF-κB genomic interactions by ChIP-seq and CUT&Tag

An understanding of the mechanisms and logic by which transcription factors coordinate gene regulation requires delineation of their genomic interactions at a genome-wide scale. Chromatin immunoprecipitation-sequencing (ChIP-seq) and more recent techniques, including CUT&Tag, typically reveal thousands of genomic interactions by transcription factors, but without insight into their functional roles. Due to cost and time considerations, optimization of ChIP experimental conditions is typically carried out only with representative interaction sites rather than through genome-wide analyses. Here, we describe insights gained from the titration of two chemical crosslinking reagents in genome-wide ChIP-seq experiments examining two members of the NF-κB family of transcription factors: RelA and c-Rel. We also describe a comparison of ChIP-seq and CUT&Tag. Our results highlight the large impact of ChIP-seq experimental conditions on the number of interactions detected, on the enrichment of consensus and non-consensus DNA motifs for the factor, and on the frequency with which the genomic interactions detected are located near potential target genes. We also found considerable consistency between ChIP-seq and CUT&Tag results, but with a substantial fraction of genomic interactions detected with only one of the two techniques. Together, the results demonstrate the dramatic impact of experimental conditions on the results obtained in a genome-wide analysis of transcription factor binding, highlighting the need for further scrutiny of the functional significance of these condition-dependent differences.

FGL1: a novel biomarker and target for non-small cell lung cancer, promoting tumor progression and metastasis through KDM4A/STAT3 transcription mechanism

Non-small cell lung cancer (NSCLC) is characterized by a high incidence rate and poor prognosis worldwide. A deeper insight into the pathogenesis of NSCLC and identification of novel therapeutic targets are essential to improve the prognosis of NSCLC. In this study, we revealed that fibrinogen-like protein 1 (FGL1) promotes proliferation, migration, and invasion of NSCLC cells. Mechanistically, we found that Stat3 acts as a transcription factor and can be recruited to the FGL1 promoter, enhancing FGL1 promoter activity. Lysine-specific demethylase 4A (KDM4A) interacts with Stat3 and facilitates the removal of methyl groups from H3K9me3, thereby enhancing Stat3-mediated transcription of FGL1. Furthermore, we observed that Stat3 and KDM4A promote NSCLC cell proliferation, migration, and invasion partly by upregulating FGL1 expression. Additionally, the expression of FGL1 was significantly higher in cancer tissues (n = 90) than in adjacent non-cancerous tissues (n = 90). Furthermore, patients with high FGL1 expression had a shorter overall survival (OS) compared to those with low FGL1 expression. We measured the expression levels of FGL1 on circulating tumor cells (CTCs) in 65 patients and found that patients with a dynamic decrease in FGL1 expression on CTCs exhibited a better therapeutic response. These findings suggest that the dynamic changes in FGL1 expression can serve as a potential biomarker for predicting treatment efficacy in NSCLC. Overall, this study revealed the significant role and regulatory mechanisms of FGL1 in the development of NSCLC, suggesting its potential as a therapeutic target for patients with NSCLC. Future studies should provide more personalized and effective treatment options for patients with NSCLC to improve clinical outcomes.

Human RNA Polymerase II Segregates from Genes and Nascent RNA and Transcribes in the Presence of DNA-Bound dCas9

RNA polymerase II (Pol II) dysfunction is frequently implied in human disease. Understanding its functional mechanism is essential for designing innovative therapeutic strategies. To visualize its supra-molecular interactions with genes and nascent RNA, we generated a human cell line carrying ~335 consecutive copies of a recombinant β-globin gene. Confocal microscopy showed that Pol II was not homogeneously concentrated around these identical gene copies. Moreover, Pol II signals partially overlapped with the genes and their nascent RNA, revealing extensive compartmentalization. Using a cell line carrying a single copy of the β-globin gene, we also tested if the binding of catalytically dead CRISPR-associated system 9 (dCas9) to different gene regions affected Pol II transcriptional activity. We assessed Pol II localization and nascent RNA levels using chromatin immunoprecipitation and droplet digital reverse transcription PCR, respectively. Some enrichment of transcriptionally paused Pol II accumulated in the promoter region was detected in a strand-specific way of gRNA binding, and there was no decrease in nascent RNA levels. Pol II preserved its transcriptional activity in the presence of DNA-bound dCas9. Our findings contribute further insight into the complex mechanism of mRNA transcription in human cells.

CX-5461 Preferentially Induces Top2α-Dependent DNA Breaks at Ribosomal DNA Loci

While genotoxic chemotherapeutic agents are among the most effective tools to combat cancer, they are often associated with severe adverse effects caused by indiscriminate DNA damage in non-tumor tissue as well as increased risk of secondary carcinogenesis. This study builds on our previous work demonstrating that the RNA Polymerase I (Pol I) transcription inhibitor CX-5461 elicits a non-canonical DNA damage response and our discovery of a critical role for Topoisomerase 2α (Top2α) in the initiation of Pol I-dependent transcription. Here, we identify Top2α as a mediator of CX-5461 response in the murine Eµ-Myc B lymphoma model whereby sensitivity to CX-5461 is dependent on cellular Top2α expression/activity. Most strikingly, and in contrast to canonical Top2α poisons, we found that the Top2α-dependent DNA damage induced by CX-5461 is preferentially localized at the ribosomal DNA (rDNA) promoter region, thereby highlighting CX-5461 as a loci-specific DNA damaging agent. This mechanism underpins the efficacy of CX-5461 against certain types of cancer and can be used to develop effective non-genotoxic anticancer drugs.

The PCNA-Pol δ complex couples lagging strand DNA synthesis to parental histone transfer for epigenetic inheritance

Inheritance of epigenetic information is critical for maintaining cell identity. The transfer of parental histone H3-H4 tetramers, the primary carrier of epigenetic modifications on histone proteins, represents a crucial yet poorly understood step in the inheritance of epigenetic information. Here, we show the lagging strand DNA polymerase, Pol δ, interacts directly with H3-H4 and that the interaction between Pol δ and the sliding clamp PCNA regulates parental histone transfer to lagging strands, most likely independent of their roles in DNA synthesis. When combined, mutations at Pol δ and Mcm2 that compromise parental histone transfer result in a greater reduction in nucleosome occupancy at nascent chromatin than mutations in either alone. Last, PCNA contributes to nucleosome positioning on nascent chromatin. On the basis of these results, we suggest that the PCNA-Pol δ complex couples lagging strand DNA synthesis to parental H3-H4 transfer, facilitating epigenetic inheritance.

Defective transfer of parental histone decreases frequency of homologous recombination by increasing free histone pools in budding yeast

Recycling of parental histones is an important step in epigenetic inheritance. During DNA replication, DNA polymerase epsilon subunit DPB3/DPB4 and DNA replication helicase subunit MCM2 are involved in the transfer of parental histones to the leading and lagging strands, respectively. Single Dpb3 deletion (dpb3Δ) or Mcm2 mutation (mcm2-3A), which each disrupts one parental histone transfer pathway, leads to the other's predominance. However, the biological impact of the two histone transfer pathways on chromatin structure and DNA repair remains elusive. In this study, we used budding yeast Saccharomyces cerevisiae to determine the genetic and epigenetic outcomes from disruption of parental histone H3-H4 tetramer transfer. We found that a dpb3Δ mcm2-3A double mutant did not exhibit the asymmetric parental histone patterns caused by a single dpb3Δ or mcm2-3A mutation, suggesting that the processes by which parental histones are transferred to the leading and lagging strands are independent. Surprisingly, the frequency of homologous recombination was significantly lower in dpb3Δ, mcm2-3A and dpb3Δ mcm2-3A mutants, likely due to the elevated levels of free histones detected in the mutant cells. Together, these findings indicate that proper transfer of parental histones during DNA replication is essential for maintaining chromatin structure and that lower homologous recombination activity due to parental histone transfer defects is detrimental to cells.

WSB1/2 target chromatin-bound lysine-methylated RelA for proteasomal degradation and NF-κB termination

Proteasome-mediated degradation of chromatin-bound NF-κB is critical in terminating the transcription of pro-inflammatory genes and can be triggered by Set9-mediated lysine methylation of the RelA subunit. However, the E3 ligase targeting methylated RelA remains unknown. Here, we find that two structurally similar substrate-recognizing components of Cullin-RING E3 ligases, WSB1 and WSB2, can recognize chromatin-bound methylated RelA for polyubiquitination and proteasomal degradation. We showed that WSB1/2 negatively regulated a subset of NF-κB target genes via associating with chromatin where they targeted methylated RelA for ubiquitination, facilitating the termination of NF-κB-dependent transcription. WSB1/2 specifically interacted with methylated lysines (K) 314 and 315 of RelA via their N-terminal WD-40 repeat (WDR) domains, thereby promoting ubiquitination of RelA. Computational modeling further revealed that a conserved aspartic acid (D) at position 158 within the WDR domain of WSB2 coordinates K314/K315 of RelA, with a higher affinity when either of the lysines is methylated. Mutation of D158 abolished WSB2's ability to bind to and promote ubiquitination of methylated RelA. Together, our study identifies a novel function and the underlying mechanism for WSB1/2 in degrading chromatin-bound methylated RelA and preventing sustained NF-κB activation, providing potential new targets for therapeutic intervention of NF-κB-mediated inflammatory diseases.

DUX4-induced HSATII transcription causes KDM2A/B-PRC1 nuclear foci and impairs DNA damage response

Polycomb repressive complexes regulate developmental gene programs, promote DNA damage repair, and mediate pericentromeric satellite repeat repression. Expression of pericentromeric satellite repeats has been implicated in several cancers and diseases, including facioscapulohumeral dystrophy (FSHD). Here, we show that DUX4-mediated transcription of HSATII regions causes nuclear foci formation of KDM2A/B-PRC1 complexes, resulting in a global loss of PRC1-mediated monoubiquitination of histone H2A. Loss of PRC1-ubiquitin signaling severely impacts DNA damage response. Our data implicate DUX4-activation of HSATII and sequestration of KDM2A/B-PRC1 complexes as a mechanism of regulating epigenetic and DNA repair pathways.

Analysis of gliomas DNA methylation: Assessment of pre-analytical variables

Precision oncology is driven by molecular biomarkers. For glioblastoma multiforme (GBM), the most common malignant adult primary brain tumor, O6-methylguanine-DNA methyltransferase ( MGMT ) gene DNA promoter methylation is an important prognostic and treatment clinical biomarker. Time consuming pre-analytical steps such as biospecimen storage before fixing, sampling, and processing are major sources of errors and batch effects, that are further confounded by intra-tumor heterogeneity of MGMT promoter methylation. To assess the effect of pre-analytical variables on GBM DNA methylation, tissue storage/sampling (CryoGrid), sample preparation multi-sonicator (PIXUL) and 5-methylcytosine (5mC) DNA immunoprecipitation (Matrix MeDIP-qPCR/seq) platforms were used. MGMT promoter CpG methylation was examined in 173 surgical samples from 90 individuals, 50 of these were used for intra-tumor heterogeneity studies. MGMT promoter methylation levels in paired frozen and formalin fixed paraffin embedded (FFPE) samples were very close, confirming suitability of FFPE for MGMT promoter methylation analysis in clinical settings. Matrix MeDIP-qPCR yielded similar results to methylation specific PCR (MS-PCR). Warm ex-vivo ischemia (37°C up to 4hrs) and 3 cycles of repeated sample thawing and freezing did not alter 5mC levels at MGMT promoter, exon and upstream enhancer regions, demonstrating the resistance of DNA methylation to the most common variations in sample processing conditions that might be encountered in research and clinical settings. 20-30% of specimens exhibited intratumor heterogeneity in the MGMT DNA promoter methylation. Collectively these data demonstrate that variations in sample fixation, ischemia duration and temperature, and DNA methylation assay technique do not have significant impact on assessment of MGMT promoter methylation status. However, intratumor methylation heterogeneity underscores the need for histologic verification and value of multiple biopsies at different GBM geographic tumor sites in assessment of MGMT promoter methylation. Matrix-MeDIP-seq analysis revealed that MGMT promoter methylation status clustered with other differentially methylated genomic loci (e.g. HOXA and lncRNAs), that are likewise resilient to variation in above post-resection pre-analytical conditions. These MGMT -associated global DNA methylation patterns offer new opportunities to validate more granular data-based epigenetic GBM clinical biomarkers where the CryoGrid-PIXUL-Matrix toolbox could prove to be useful.

The ARID1A-METTL3-m6A axis ensures effective RNase H1-mediated resolution of R-loops and genome stability

R-loops are three-stranded structures that can pose threats to genome stability. RNase H1 precisely recognizes R-loops to drive their resolution within the genome, but the underlying mechanism is unclear. Here, we report that ARID1A recognizes R-loops with high affinity in an ATM-dependent manner. ARID1A recruits METTL3 and METTL14 to the R-loop, leading to the m6A methylation of R-loop RNA. This m6A modification facilitates the recruitment of RNase H1 to the R-loop, driving its resolution and promoting DNA end resection at DSBs, thereby ensuring genome stability. Depletion of ARID1A, METTL3, or METTL14 leads to R-loop accumulation and reduced cell survival upon exposure to cytotoxic agents. Therefore, ARID1A, METTL3, and METTL14 function in a coordinated, temporal order at DSB sites to recruit RNase H1 and to ensure efficient R-loop resolution. Given the association of high ARID1A levels with resistance to genotoxic therapies in patients, these findings open avenues for exploring potential therapeutic strategies for cancers with ARID1A abnormalities.

The Role of Hexokinases in Epigenetic Regulation: Altered Hexokinase Expression and Chromatin Stability in Yeast

Background . Human hexokinase 2 ( HK2 ) plays an important role in regulating Warburg effect, which metabolizes glucose to lactate acid even in the presence of ample oxygen and provides intermediate metabolites to support cancer cell proliferation and tumor growth. HK2 overexpression has been observed in various types of cancers and targeting HK2 -driven Warburg effect has been suggested as a potential cancer therapeutic strategy. Given that epigenetic enzymes utilize metabolic intermediates as substrates or co-factors to carry out post-translational modification of DNA and histones in cells, we hypothesized that altering HK2 expression-mediated cellular glycolysis rates could impact the epigenome and, consequently, genome stability in yeast. To test this hypothesis, we established genetic models with different yeast hexokinase 2 ( HXK2) expression in Saccharomyces cerevisiae yeast cells and investigated the effect of HXK2 -dependent metabolism on parental nucleosome transfer, a key DNA replication-coupled epigenetic inheritance process, and chromatin stability. Results . By comparing the growth of mutant yeast cells carrying single deletion of hxk1Δ , hxk2Δ , or double-loss of hxk1Δ hxk2Δ to wild-type cells, we demonstrated that HXK2 is the dominant HXK in yeast cell growth. Surprisingly, manipulating HXK2 expression in yeast, whether through overexpression or deletion, had only a marginal impact on parental nucleosome assembly, but a noticeable trend with decrease chromatin instability. However, targeting yeast cells with 2-deoxy-D-glucose (2-DG), a HK2 inhibitor that has been proposed as an anti-cancer treatment, significantly increased chromatin instability. Conclusion . Our findings suggest that in yeast cells lacking HXK2 , alternative HXK s such as HXK1 or glucokinase 1 ( GLK1 ) play a role in supporting glycolysis at a level that adequately maintain epigenomic stability. While our study demonstrated an increase in epigenetic instability with 2-DG treatment, the observed effect seemed to occur independently of Hxk2-mediated glycolysis inhibition. Thus, additional research is needed to identify the molecular mechanism through which 2-DG influences chromatin stability.

Insulin Determines Transforming Growth Factor β Effects on Hepatocyte Nuclear Factor 4α Transcription in Hepatocytes

Loss of hepatocyte nuclear factor 4α (HNF4α) expression is frequently observed in end-stage liver disease and associated with loss of vital liver functions, thus increasing mortality. Loss of HNF4α expression is mediated by inflammatory cytokines, such as transforming growth factor (TGF)-β. However, details of how HNF4α is suppressed are largely unknown to date. Herein, TGF-β did not directly inhibit HNF4α but contributed to its transcriptional regulation by SMAD2/3 recruiting acetyltransferase CREB-binding protein/p300 to the HNF4α promoter. The recruitment of CREB-binding protein/p300 is indispensable for CCAAT/enhancer-binding protein α (C/EBPα) binding, another essential requirement for constitutive HNF4α expression in hepatocytes. Consistent with the in vitro observation, 67 of 98 patients with hepatic HNF4α expressed both phospho-SMAD2 and C/EBPα, whereas 22 patients without HNF4α expression lacked either phospho-SMAD2 or C/EBPα. In contrast to the observed induction of HNF4α, SMAD2/3 inhibited C/EBPα transcription. Long-term TGF-β incubation resulted in C/EBPα depletion, which abrogated HNF4α expression. Intriguingly, SMAD2/3 inhibitory binding to the C/EBPα promoter was abolished by insulin. Two-thirds of patients without C/EBPα lacked membrane glucose transporter type 2 expression in hepatocytes, indicating insulin resistance. Taken together, these data indicate that hepatic insulin sensitivity is essential for hepatic HNF4α expression in the condition of inflammation.

MYSM1 acts as a novel co-activator of ERα to confer antiestrogen resistance in breast cancer

Endocrine resistance is a crucial challenge in estrogen receptor alpha (ERα)-positive breast cancer (BCa). Aberrant alteration in modulation of E2/ERα signaling pathway has emerged as the putative contributor for endocrine resistance in BCa. Herein, we demonstrate that MYSM1 as a deubiquitinase participates in modulating ERα action via histone and non-histone deubiquitination. MYSM1 is involved in maintenance of ERα stability via ERα deubiquitination. MYSM1 regulates relevant histone modifications on cis regulatory elements of ERα-regulated genes, facilitating chromatin decondensation. MYSM1 is highly expressed in clinical BCa samples. MYSM1 depletion attenuates BCa-derived cell growth in xenograft models and increases the sensitivity of antiestrogen agents in BCa cells. A virtual screen shows that the small molecule Imatinib could potentially interact with catalytic MPN domain of MYSM1 to inhibit BCa cell growth via MYSM1-ERα axis. These findings clarify the molecular mechanism of MYSM1 as an epigenetic modifier in regulation of ERα action and provide a potential therapeutic target for endocrine resistance in BCa.

Wnt16 Promotes Vascular Smooth Muscle Contractile Phenotype and Function via Taz (Wwtr1) Activation in Male LDLR-/- Mice

Wnt16 is expressed in bone and arteries, and maintains bone mass in mice and humans, but its role in cardiovascular physiology is unknown. We show that Wnt16 protein accumulates in murine and human vascular smooth muscle (VSM). WNT16 genotypes that convey risk for bone frailty also convey risk for cardiovascular events in the Dallas Heart Study. Murine Wnt16 deficiency, which causes postnatal bone loss, also reduced systolic blood pressure. Electron microscopy demonstrated abnormal VSM mitochondrial morphology in Wnt16-null mice, with reductions in mitochondrial respiration. Following angiotensin-II (AngII) infusion, thoracic ascending aorta (TAA) dilatation was greater in Wnt16-/- vs Wnt16+/+ mice (LDLR-/- background). Acta2 (vascular smooth muscle alpha actin) deficiency has been shown to impair contractile phenotype and worsen TAA aneurysm with concomitant reductions in blood pressure. Wnt16 deficiency reduced expression of Acta2, SM22 (transgelin), and other contractile genes, and reduced VSM contraction induced by TGFβ. Acta2 and SM22 proteins were reduced in Wnt16-/- VSM as was Ankrd1, a prototypic contractile target of Yap1 and Taz activation via TEA domain (TEAD)-directed transcription. Wnt16-/- VSM exhibited reduced nuclear Taz and Yap1 protein accumulation. SiRNA targeting Wnt16 or Taz, but not Yap1, phenocopied Wnt16 deficiency, and Taz siRNA inhibited contractile gene upregulation by Wnt16. Wnt16 incubation stimulated mitochondrial respiration and contraction (reversed by verteporfin, a Yap/Taz inhibitor). SiRNA targeting Taz inhibitors Ccm2 and Lats1/2 mimicked Wnt16 treatment. Wnt16 stimulated Taz binding to Acta2 chromatin and H3K4me3 methylation. TEAD cognates in the Acta2 promoter conveyed transcriptional responses to Wnt16 and Taz. Wnt16 regulates cardiovascular physiology and VSM contractile phenotype, mediated via Taz signaling.

Biofortification of iron content by regulating a NAC transcription factor in maize

Iron (Fe) deficiency remains widespread among people in developing countries. To help solve this problem, breeders have been attempting to develop maize cultivars with high yields and high Fe concentrations in the kernels. We conducted a genome-wide association study and identified a gene, ZmNAC78 (NAM/ATAF/CUC DOMAIN TRANSCRIPTION FACTOR 78), that regulates Fe concentrations in maize kernels. We cultivated maize varieties with both high yield and high Fe concentrations in their kernels by using a molecular marker developed from a 42-base pair insertion or deletion (indel) in the promoter of ZmNAC78. ZmNAC78 expression is enriched in the basal endosperm transfer layer of kernels, and the ZmNAC78 protein directly regulates messenger RNA abundance of Fe transporters. Our results thus provide an approach to develop maize varieties with Fe-enriched kernels.

PARylated PDHE1α generates acetyl-CoA for local chromatin acetylation and DNA damage repair

Chromatin relaxation is a prerequisite for the DNA repair machinery to access double-strand breaks (DSBs). Local histones around the DSBs then undergo prompt changes in acetylation status, but how the large demands of acetyl-CoA are met is unclear. Here, we report that pyruvate dehydrogenase 1α (PDHE1α) catalyzes pyruvate metabolism to rapidly provide acetyl-CoA in response to DNA damage. We show that PDHE1α is quickly recruited to chromatin in a polyADP-ribosylation-dependent manner, which drives acetyl-CoA generation to support local chromatin acetylation around DSBs. This process increases the formation of relaxed chromatin to facilitate repair-factor loading, genome stability and cancer cell resistance to DNA-damaging treatments in vitro and in vivo. Indeed, we demonstrate that blocking polyADP-ribosylation-based PDHE1α chromatin recruitment attenuates chromatin relaxation and DSB repair efficiency, resulting in genome instability and restored radiosensitivity. These findings support a mechanism in which chromatin-associated PDHE1α locally generates acetyl-CoA to remodel the chromatin environment adjacent to DSBs and promote their repair.

ZEB1 promotes non-homologous end joining double-strand break repair

Repair of DSB induced by IR is primarily carried out by Non-Homologous End Joining (NHEJ), a pathway in which 53BP1 plays a key role. We have discovered that the EMT-inducing transcriptional repressor ZEB1 (i) interacts with 53BP1 and that this interaction occurs rapidly and is significantly amplified following exposure of cells to IR; (ii) is required for the localization of 53BP1 to a subset of double-stranded breaks, and for physiological DSB repair; (iii) co-localizes with 53BP1 at IR-induced foci (IRIF); (iv) promotes NHEJ and inhibits Homologous Recombination (HR); (v) depletion increases resection at DSBs and (vi) confers PARP inhibitor (PARPi) sensitivity on BRCA1-deficient cells. Lastly, ZEB1's effects on repair pathway choice, resection, and PARPi sensitivity all rely on its homeodomain. In contrast to the well-characterized therapeutic resistance of high ZEB1-expressing cancer cells, the novel ZEB1-53BP1-shieldin resection axis described here exposes a therapeutic vulnerability: ZEB1 levels in BRCA1-deficient tumors may serve as a predictive biomarker of response to PARPis.

CRISPR-ChIP reveals selective regulation of H3K79me2 by Menin in MLL leukemia

Chromatin regulation involves the selective recruitment of chromatin factors to facilitate DNA repair, replication and transcription. Here we demonstrate the utility of coupling unbiased functional genomics with chromatin immunoprecipitation (CRISPR-ChIP) to identify the factors associated with active chromatin modifications in mammalian cells. Specifically, an integrated reporter containing a cis-regulatory element of interest and a single guide RNA provide a chromatinized template for a direct readout for regulators of histone modifications associated with actively transcribed genes such as H3K4me3 and H3K79me2. With CRISPR-ChIP, we identify all the nonredundant COMPASS complex members required for H3K4me3 and demonstrate that RNA polymerase II is dispensable for the maintenance of H3K4me3. As H3K79me2 has a putative oncogenic function in leukemia cells driven by MLL translocations, using CRISPR-ChIP we reveal a functional partitioning of H3K79 methylation into two distinct regulatory units: an oncogenic DOT1L complex directed by the MLL fusion protein in a Menin-dependent manner and a separate endogenous DOT1L complex, where catalytic activity is directed by MLLT10. Overall, CRISPR-ChIP provides a powerful tool for the unbiased interrogation of the mechanisms underpinning chromatin regulation.

SMCHD1 has separable roles in chromatin architecture and gene silencing that could be targeted in disease

The interplay between 3D chromatin architecture and gene silencing is incompletely understood. Here, we report a novel point mutation in the non-canonical SMC protein SMCHD1 that enhances its silencing capacity at endogenous developmental targets. Moreover, it also results in enhanced silencing at the facioscapulohumeral muscular dystrophy associated macrosatellite-array, D4Z4, resulting in enhanced repression of DUX4 encoded by this repeat. Heightened SMCHD1 silencing perturbs developmental Hox gene activation, causing a homeotic transformation in mice. Paradoxically, the mutant SMCHD1 appears to enhance insulation against other epigenetic regulators, including PRC2 and CTCF, while depleting long range chromatin interactions akin to what is observed in the absence of SMCHD1. These data suggest that SMCHD1's role in long range chromatin interactions is not directly linked to gene silencing or insulating the chromatin, refining the model for how the different levels of SMCHD1-mediated chromatin regulation interact to bring about gene silencing in normal development and disease.

SnRK1 inhibits anthocyanin biosynthesis through both transcriptional regulation and direct phosphorylation and dissociation of the MYB/bHLH/TTG1 MBW complex

Plants have evolved an extensive specialized secondary metabolism. The colorful flavonoid anthocyanins, for example, not only stimulate flower pollination and seed dispersal, but also protect different tissues against high light, UV and oxidative stress. Their biosynthesis is highly regulated by environmental and developmental cues and induced by high sucrose levels. Expression of the biosynthetic enzymes involved is controlled by a transcriptional MBW complex, comprising (R2R3) MYB- and bHLH-type transcription factors and the WD40 repeat protein TTG1. Anthocyanin biosynthesis is not only useful, but also carbon- and energy-intensive and non-vital. Consistently, the SnRK1 protein kinase, a metabolic sensor activated in carbon- and energy-depleting stress conditions, represses anthocyanin biosynthesis. Here we show that Arabidopsis SnRK1 represses MBW complex activity both at the transcriptional and post-translational level. In addition to repressing expression of the key transcription factor MYB75/PAP1, SnRK1 activity triggers MBW complex dissociation, associated with loss of target promoter binding, MYB75 protein degradation and nuclear export of TTG1. We also provide evidence for direct interaction with and phosphorylation of multiple MBW complex proteins. These results indicate that repression of expensive anthocyanin biosynthesis is an important strategy to save energy and redirect carbon flow to more essential processes for survival in metabolic stress conditions.

Protein-DNA Binding Residues Prediction Using a Deep Learning Model With Hierarchical Feature Extraction

Biologically important effects occur when proteins bind to other substances, of which binding to DNA is a crucial one. Therefore, accurate identification of protein-DNA binding residues is important for further understanding of the protein-DNA interaction mechanism. Although wet-lab methods can accurately obtain the location of bound residues, it requires significant human, financial and time costs. There is thus an urgent need to develop efficient computational-based methods. Most current state-of-the-art methods are two-step approaches: the first step uses a sliding window technique to extract residue features; the second step uses each residue as an input to the model for prediction. This has a negative impact on the efficiency of prediction and ease of use. In this study, we propose a sequence-to-sequence (seq2seq) model that can input the entire protein sequence of variable length and use two modules, Transformer Encoder Block and Feature Extracting Block, for hierarchical feature extraction, where Transformer Encoder Block is used to extract global features, and then Feature Extracting Block is used to extract local features to further improve the recognition capability of the model. The comparison results on two benchmark datasets, namely PDNA-543 and PDNA-41, prove the effectiveness of our method in identifying protein-DNA binding residues.

Comparative genomic analyses reveal the genetic basis of the yellow-seed trait in Brassica napus

Yellow-seed trait is a desirable breeding characteristic of rapeseed (Brassica napus) that could greatly improve seed oil yield and quality. However, the underlying mechanisms controlling this phenotype in B. napus plants are difficult to discern because of their complexity. Here, we assemble high-quality genomes of yellow-seeded (GH06) and black-seeded (ZY821). Combining in-depth fine mapping of a quantitative trait locus (QTL) for seed color with other omics data reveal BnA09MYB47a, encoding an R2R3-MYB-type transcription factor, as the causal gene of a major QTL controlling the yellow-seed trait. Functional studies show that sequence variation of BnA09MYB47a underlies the functional divergence between the yellow- and black-seeded B. napus. The black-seed allele BnA09MYB47aZY821, but not the yellow-seed allele BnA09MYB47aGH06, promotes flavonoid biosynthesis by directly activating the expression of BnTT18. Our discovery suggests a possible approach to breeding B. napus for improved commercial value and facilitates flavonoid biosynthesis studies in Brassica crops.

The basic helix-loop-helix transcription factor gene, OsbHLH38, plays a key role in controlling rice salt tolerance

The plant hormone abscisic acid (ABA) is crucial for plant seed germination and abiotic stress tolerance. However, the association between ABA sensitivity and plant abiotic stress tolerance remains largely unknown. In this study, 436 rice accessions were assessed for their sensitivity to ABA during seed germination. The considerable diversity in ABA sensitivity among rice germplasm accessions was primarily reflected by the differentiation between the Xian (indica) and Geng (japonica) subspecies and between the upland-Geng and lowland-Geng ecotypes. The upland-Geng accessions were most sensitive to ABA. Genome-wide association analyses identified four major quantitative trait loci containing 21 candidate genes associated with ABA sensitivity of which a basic helix-loop-helix transcription factor gene, OsbHLH38, was the most important for ABA sensitivity. Comprehensive functional analyses using knockout and overexpression transgenic lines revealed that OsbHLH38 expression was responsive to multiple abiotic stresses. Overexpression of OsbHLH38 increased seedling salt tolerance, while knockout of OsbHLH38 increased sensitivity to salt stress. A salt-responsive transcription factor, OsDREB2A, interacted with OsbHLH38 and was directly regulated by OsbHLH38. Moreover, OsbHLH38 affected rice abiotic stress tolerance by mediating the expression of a large set of transporter genes of phytohormones, transcription factor genes, and many downstream genes with diverse functions, including photosynthesis, redox homeostasis, and abiotic stress responsiveness. These results demonstrated that OsbHLH38 is a key regulator in plant abiotic stress tolerance.

Integrated multiomic analysis and high-throughput screening reveal potential gene targets and synergetic drug combinations for osteosarcoma therapy

Although great advances have been made over the past decades, therapeutics for osteosarcoma are quite limited. We performed long-read RNA sequencing and tandem mass tag (TMT)-based quantitative proteome on osteosarcoma and the adjacent normal tissues, next-generation sequencing (NGS) on paired osteosarcoma samples before and after neoadjuvant chemotherapy (NACT), and high-throughput drug combination screen on osteosarcoma cell lines. Single-cell RNA sequencing data were analyzed to reveal the heterogeneity of potential therapeutic target genes. Additionally, we clarified the synergistic mechanisms of doxorubicin (DOX) and HDACs inhibitors for osteosarcoma treatment. Consequently, we identified 2535 osteosarcoma-specific genes and several alternative splicing (AS) events with osteosarcoma specificity and/or patient heterogeneity. Hundreds of potential therapeutic targets were identified among them, which showed the core regulatory roles in osteosarcoma. We also identified 215 inhibitory drugs and 236 synergistic drug combinations for osteosarcoma treatment. More interestingly, the multiomic analysis pointed out the pivotal role of HDAC1 and TOP2A in osteosarcoma. HDAC inhibitors synergized with DOX to suppress osteosarcoma both in vitro and in vivo. Mechanistically, HDAC inhibitors synergized with DOX by downregulating SP1 to transcriptionally modulate TOP2A expression. This study provided a comprehensive view of molecular features, therapeutic targets, and synergistic drug combinations for osteosarcoma.

Calcineurin regulates aldosterone production via dephosphorylation of NFATC4

The mineralocorticoid aldosterone, secreted by the adrenal zona glomerulosa (ZG), is critical for life, maintaining ion homeostasis and blood pressure. Therapeutic inhibition of protein phosphatase 3 (calcineurin, Cn) results in inappropriately low plasma aldosterone levels despite concomitant hyperkalemia and hyperreninemia. We tested the hypothesis that Cn participates in the signal transduction pathway regulating aldosterone synthesis. Inhibition of Cn with tacrolimus abolished the potassium-stimulated (K+-stimulated) expression of aldosterone synthase, encoded by CYP11B2, in the NCI-H295R human adrenocortical cell line as well as ex vivo in mouse and human adrenal tissue. ZG-specific deletion of the regulatory Cn subunit CnB1 diminished Cyp11b2 expression in vivo and disrupted K+-mediated aldosterone synthesis. Phosphoproteomics analysis identified nuclear factor of activated T cells, cytoplasmic 4 (NFATC4), as a target for Cn-mediated dephosphorylation. Deletion of NFATC4 impaired K+-dependent stimulation of CYP11B2 expression and aldosterone production while expression of a constitutively active form of NFATC4 increased expression of CYP11B2 in NCI-H295R cells. Chromatin immunoprecipitation revealed NFATC4 directly regulated CYP11B2 expression. Thus, Cn controls aldosterone production via the Cn/NFATC4 pathway. Inhibition of Cn/NFATC4 signaling may explain low plasma aldosterone levels and hyperkalemia in patients treated with tacrolimus, and the Cn/NFATC4 pathway may provide novel molecular targets to treat primary aldosteronism.

LAPTM4B-YAP loop feedback amplification enhances the stemness of hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is highly heterogeneous, and stemness signatures are frequently elevated in HCC tumor cells to generate heterogeneous subtypes via multidirectional differentiation. However, the mechanisms affecting the regulation of stemness in HCC remain unclear. In this study, we identified that lysosome-associated protein transmembrane-4β (LAPTM4B) was significantly overexpressed in stem-like tumor cell populations with multidirectional differentiation potential at the single cell level, and verified that LAPTM4B was closely related to stemness of HCC using in vitro and in vivo experiments. Mechanistically, elevated LAPTM4B suppresses Yes-associated protein (YAP) phosphorylation and ubiquitination degradation. In turn, stabilized YAP localizes to the nucleus and binds to cAMP responsive element binding protein-1 (CREB1), which promotes transcription of LAPTM4B. Overall, our findings suggest that LAPTM4B forms a positive feedback loop with YAP, which maintains the stemness of HCC tumor cells and leads to an unfavorable prognosis for HCC patients.

Functional and structural diversification of incomplete phosphotransferase system in cellulose-degrading clostridia

Carbohydrate utilization is critical to microbial survival. The phosphotransferase system (PTS) is a well-documented microbial system with a prominent role in carbohydrate metabolism, which can transport carbohydrates through forming a phosphorylation cascade and regulate metabolism by protein phosphorylation or interactions in model strains. However, those PTS-mediated regulated mechanisms have been underexplored in non-model prokaryotes. Here, we performed massive genome mining for PTS components in nearly 15,000 prokaryotic genomes from 4,293 species and revealed a high prevalence of incomplete PTSs in prokaryotes with no association to microbial phylogeny. Among these incomplete PTS carriers, a group of lignocellulose degrading clostridia was identified to have lost PTS sugar transporters and carry a substitution of the conserved histidine residue in the core PTS component, HPr (histidine-phosphorylatable phosphocarrier). Ruminiclostridium cellulolyticum was then selected as a representative to interrogate the function of incomplete PTS components in carbohydrate metabolism. Inactivation of the HPr homolog reduced rather than increased carbohydrate utilization as previously indicated. In addition to regulating distinct transcriptional profiles, PTS associated CcpA (Catabolite Control Protein A) homologs diverged from previously described CcpA with varied metabolic relevance and distinct DNA binding motifs. Furthermore, the DNA binding of CcpA homologs is independent of HPr homolog, which is determined by structural changes at the interface of CcpA homologs, rather than in HPr homolog. These data concordantly support functional and structural diversification of PTS components in metabolic regulation and bring novel understanding of regulatory mechanisms of incomplete PTSs in cellulose-degrading clostridia.

BAP18 facilitates CTCF-mediated chromatin accessible to regulate enhancer activity in breast cancer

The estrogen receptor alpha (ERα) signaling pathway is a crucial target for ERα-positive breast cancer therapeutic strategies. Co-regulators and other transcription factors cooperate for effective ERα-related enhancer activation. Recent studies demonstrate that the transcription factor CTCF is essential to participate in ERα/E2-induced enhancer transactivation. However, the mechanism of how CTCF is achieved remains unknown. Here, we provided evidence that BAP18 is required for CTCF recruitment on ERα-enriched enhancers, facilitating CTCF-mediated chromatin accessibility to promote enhancer RNAs transcription. Consistently, GRO-seq demonstrates that the enhancer activity is positively correlated with BAP18 enrichment. Furthermore, BAP18 interacts with SMARCA1/BPTF to accelerate the recruitment of CTCF to ERα-related enhancers. Interestingly, BAP18 is involved in chromatin accessibility within enhancer regions, thereby increasing enhancer transactivation and enhancer-promoter looping. BAP18 depletion increases the sensitivity of anti-estrogen and anti-enhancer treatment in MCF7 cells. Collectively, our study indicates that BAP18 coordinates with CTCF to enlarge the transactivation of ERα-related enhancers, providing a better understanding of BAP18/CTCF coupling chromatin remodeling and E-P looping in the regulation of enhancer transcription.

Human DUX4 and mouse Dux interact with STAT1 and broadly inhibit interferon-stimulated gene induction

DUX4 activates the first wave of zygotic gene expression in the early embryo. Mis-expression of DUX4 in skeletal muscle causes facioscapulohumeral dystrophy (FSHD), whereas expression in cancers suppresses IFNγ induction of major histocompatibility complex class I (MHC class I) and contributes to immune evasion. We show that the DUX4 protein interacts with STAT1 and broadly suppresses expression of IFNγ-stimulated genes by decreasing bound STAT1 and Pol-II recruitment. Transcriptional suppression of interferon-stimulated genes (ISGs) requires conserved (L)LxxL(L) motifs in the carboxyterminal region of DUX4 and phosphorylation of STAT1 Y701 enhances interaction with DUX4. Consistent with these findings, expression of endogenous DUX4 in FSHD muscle cells and the CIC-DUX4 fusion containing the DUX4 CTD in a sarcoma cell line inhibit IFNγ induction of ISGs. Mouse Dux similarly interacted with STAT1 and suppressed IFNγ induction of ISGs. These findings identify an evolved role of the DUXC family in modulating immune signaling pathways with implications for development, cancers, and FSHD.

Glucocorticoid activation of anti-inflammatory macrophages protects against insulin resistance

Insulin resistance (IR) during obesity is linked to adipose tissue macrophage (ATM)-driven inflammation of adipose tissue. Whether anti-inflammatory glucocorticoids (GCs) at physiological levels modulate IR is unclear. Here, we report that deletion of the GC receptor (GR) in myeloid cells, including macrophages in mice, aggravates obesity-related IR by enhancing adipose tissue inflammation due to decreased anti-inflammatory ATM leading to exaggerated adipose tissue lipolysis and severe hepatic steatosis. In contrast, GR deletion in Kupffer cells alone does not alter IR. Co-culture experiments show that the absence of GR in macrophages directly causes reduced phospho-AKT and glucose uptake in adipocytes, suggesting an important function of GR in ATM. GR-deficient macrophages are refractory to alternative ATM-inducing IL-4 signaling, due to reduced STAT6 chromatin loading and diminished anti-inflammatory enhancer activation. We demonstrate that GR has an important function in macrophages during obesity by limiting adipose tissue inflammation and lipolysis to promote insulin sensitivity.

Mycobacterium indicus pranii therapy suppresses systemic dissemination of tumor cells in B16F10 murine model of melanoma

Cancer associated morbidity is mostly attributed to the dissemination of tumor cells from their primary niche into the circulation known as "metastasis". Mycobacterium indicus pranii (MIP) an approved immunotherapeutic agent against lung cancer (NSCLC) has shown potent anti-tumor activity in prior studies. While evaluating anti-tumor activity of MIP in mouse model, MIP treated animals typically exhibited less metastatic lesions in their pulmonary compartment. To study the role of MIP in metastasis closely, B16F10 melanoma cells were implanted subcutaneously in the mice, and the dissemination of tumor cells from the solid tumor was evaluated over a period of time. When B16F10 melanoma cells were treated with MIP in vitro, downregulation of epithelial mesenchymal transition markers was observed in these cells, which in turn suppressed the invasion, migration and adhesion of tumor cells. Notably, MIP therapy was found to be effectively reducing the metastatic burden in murine model of melanoma. Molecular characterization of MIP treated tumor cells substantiated that MIP upregulates the PPARγ expression within the tumor cells, which attenuates the NFκB/p65 levels within the nucleus, resulting in the suppression of Mmp9 expression in tumor cells. Besides that, MIP also downregulated the surface expression of chemokine receptor CXCR4 in murine melanoma cells, where chromatin immunoprecipitation confirmed the impeded recruitment of p50 and c-Rel factors to the Cxcr4 promoter, resulting in its downregulation transcriptionally. Taken together, MIP suppressed the dissemination of tumor cells in vivo, by regulating the expression of MMP9 and CXCR4 on these cells.

Development of Chromatin Immunoprecipitation for the Analysis of Histone Modifications in Red Macroalga Neopyropia yezoensis (Rhodophyta)

Epigenetic regulation by histone modification can activate or repress transcription through changes in chromatin dynamics and regulates development and the response to environmental signals in both animals and plants. Chromatin immunoprecipitation (ChIP) is an indispensable tool to identify histones with specific post-translational modifications. The lack of a ChIP technique for macroalgae has hindered understanding of the role of histone modification in the expression of genes in this organism. In this study, a ChIP method with several modifications, based on existing protocols for plant cells, has been developed for the red macroalga, Neopyropia yezoensis, that consists of a heterogeneous alternation of macroscopic leaf-like gametophytes and microscopic filamentous sporophytes. ChIP method coupled with qPCR enables the identification of a histone mark in generation-specific genes from N. yezoensis. The results indicate that acetylation of histone H3 at lysine 9 in the 5' flanking and coding regions from generation-specific genes was maintained at relatively high levels, even in generation-repressed gene expression. The use of this ChIP method will contribute significantly to identify the epigenetic regulatory mechanisms through histone modifications that control a variety of biological processes in red macroalgae.

11-Cl-BBQ, a select modulator of AhR-regulated transcription, suppresses lung cancer cell growth via activation of p53 and p27Kip1

Aryl hydrocarbon receptor (AhR) is a ligand-activated transcription factor and functions as a tumour suppressor in different cancer models. In the present study, we report detailed characterization of 11-chloro-7H-benzimidazo[2,1-a]benzo[de]iso-quinolin-7-one (11-Cl-BBQ) as a select modulator of AhR-regulated transcription (SMAhRT) with anti-cancer actions. Treatment of lung cancer cells with 11-Cl-BBQ induced potent and sustained AhR-dependent anti-proliferative effects by promoting G1 phase cell cycle arrest. Investigation of 11-Cl-BBQ-induced transcription in H460 cells with or without the AhR expression by RNA-sequencing revealed activation of p53 signalling. In addition, 11-Cl-BBQ suppressed multiple pathways involved in DNA replication and increased expression of cyclin-dependent kinase inhibitors, including p27Kip1 , in an AhR-dependent manner. CRISPR/Cas9 knockout of individual genes revealed the requirement for both p53 and p27Kip1 for the AhR-mediated anti-proliferative effects. Our results identify 11-Cl-BBQ as a potential lung cancer therapeutic, highlight the feasibility of targeting AhR and provide important mechanistic insights into AhR-mediated-anticancer actions.

Endogenous bioluminescent reporters reveal a sustained increase in utrophin gene expression upon EZH2 and ERK1/2 inhibition

Duchenne muscular dystrophy (DMD) is an X-linked disorder caused by loss of function mutations in the dystrophin gene (Dmd), resulting in progressive muscle weakening. Here we modelled the longitudinal expression of endogenous Dmd, and its paralogue Utrn, in mice and in myoblasts by generating bespoke bioluminescent gene reporters. As utrophin can partially compensate for Dmd-deficiency, these reporters were used as tools to ask whether chromatin-modifying drugs can enhance Utrn expression in developing muscle. Myoblasts treated with different PRC2 inhibitors showed significant increases in Utrn transcripts and bioluminescent signals, and these responses were independently verified by conditional Ezh2 deletion. Inhibition of ERK1/2 signalling provoked an additional increase in Utrn expression that was also seen in Dmd-mutant cells, and maintained as myoblasts differentiate. These data reveal PRC2 and ERK1/2 to be negative regulators of Utrn expression and provide specialised molecular imaging tools to monitor utrophin expression as a therapeutic strategy for DMD.

MultiomicsTracks96: A high throughput PIXUL-Matrix-based toolbox to profile frozen and FFPE tissues multiomes

Background

The multiome is an integrated assembly of distinct classes of molecules and molecular properties, or "omes," measured in the same biospecimen. Freezing and formalin-fixed paraffin-embedding (FFPE) are two common ways to store tissues, and these practices have generated vast biospecimen repositories. However, these biospecimens have been underutilized for multi-omic analysis due to the low throughput of current analytical technologies that impede large-scale studies.

Methods

Tissue sampling, preparation, and downstream analysis were integrated into a 96-well format multi-omics workflow, MultiomicsTracks96. Frozen mouse organs were sampled using the CryoGrid system, and matched FFPE samples were processed using a microtome. The 96-well format sonicator, PIXUL, was adapted to extract DNA, RNA, chromatin, and protein from tissues. The 96-well format analytical platform, Matrix, was used for chromatin immunoprecipitation (ChIP), methylated DNA immunoprecipitation (MeDIP), methylated RNA immunoprecipitation (MeRIP), and RNA reverse transcription (RT) assays followed by qPCR and sequencing. LC-MS/MS was used for protein analysis. The Segway genome segmentation algorithm was used to identify functional genomic regions, and linear regressors based on the multi-omics data were trained to predict protein expression.

Results

MultiomicsTracks96 was used to generate 8-dimensional datasets including RNA-seq measurements of mRNA expression; MeRIP-seq measurements of m6A and m5C; ChIP-seq measurements of H3K27Ac, H3K4m3, and Pol II; MeDIP-seq measurements of 5mC; and LC-MS/MS measurements of proteins. We observed high correlation between data from matched frozen and FFPE organs. The Segway genome segmentation algorithm applied to epigenomic profiles (ChIP-seq: H3K27Ac, H3K4m3, Pol II; MeDIP-seq: 5mC) was able to recapitulate and predict organ-specific super-enhancers in both FFPE and frozen samples. Linear regression analysis showed that proteomic expression profiles can be more accurately predicted by the full suite of multi-omics data, compared to using epigenomic, transcriptomic, or epitranscriptomic measurements individually.

Conclusions

The MultiomicsTracks96 workflow is well suited for high dimensional multi-omics studies - for instance, multiorgan animal models of disease, drug toxicities, environmental exposure, and aging as well as large-scale clinical investigations involving the use of biospecimens from existing tissue repositories.

MicroRNA Profiling of Self-Renewing Human Neural Stem Cells Reveals Novel Sets of Differentially Expressed microRNAs During Neural Differentiation In Vitro

The involvement of microRNAs (miRNAs) in orchestrating self-renewal and differentiation of stem cells has been revealed in a number of recent studies. And while in human pluripotent stem cells, miRNAs have been directly linked to the core pluripotency network, including the cell cycle regulation and the maintenance of the self-renewing capacity, their role in the onset of differentiation in other contexts, such as determination of neural cell fate, remains poorly described. To bridge this gap, we used three model cell types to study miRNA expression patterns: human embryonic stem cells (hESCs), hESCs-derived self-renewing neural stem cells (NSCs), and differentiating NSCs. The comprehensive miRNA profiling presented here reveals novel sets of miRNAs differentially expressed during human neural cell fate determination in vitro. Furthermore, we report a miRNA expression profile of self-renewing human NSCs, which has been lacking to this date. Our data also indicates that miRNA clusters enriched in NSCs share the target-determining seed sequence with cell cycle regulatory miRNAs expressed in pluripotent hESCs. Lastly, our mechanistic experiments confirmed that cluster miR-17-92, one of the NSCs-enriched clusters, is directly transcriptionally regulated by transcription factor c-MYC.