Citations

Abstract Epigenetic inheritance relies on the precise recycling of parental histones during DNA replication. The N-terminal region of Pol1, the catalytic subunit of DNA polymerase α, has previously been implicated in chaperoning parental histone recycling. In this study, we identify a single amino acid mutation in the C-terminal region of Saccharomyces cerevisiae Pol1 (pol1-F1396S) that destabilizes the protein and disrupts the transfer of parental histones to the lagging strand during DNA replication. As a primary consequence, pol1-F1396S cells exhibit decreased chromatin stability. Furthermore, we demonstrate that the cellular abundance—or dosage—of all three DNA polymerases (Pol α, Pol δ, and Pol ε) influences the distribution of parental histones between the leading and lagging strands, revealing a dosage-sensitive mechanism of epigenetic regulation during DNA replication.

Background and Aims: Transforming growth factor-beta 1 (TGFβ1) induces HSC activation into metastasis-promoting cancer-associated fibroblasts (CAFs), but how the process is fueled remains incompletely understood. We studied metabolic reprogramming induced by TGFβ1 in HSCs. Approaches and Results: Activation of cultured primary human HSCs was assessed by the expression of myofibroblast markers. Glucose transporter 1 (Glut1) of murine HSC was disrupted by Cre recombinase/LoxP sequence derived from bacteriophage P1 recombination (Cre/LoxP). Plasma membrane (PM) Glut1 and glycolysis were studied by biotinylation assay and the Angilent Seahorse XFe96 Analyzer. S.c. HSC/tumor co-implantation and portal vein injection of MC38 colorectal cancer cells into HSC-specific Glut1 knockout mice were performed to determine in vivo relevance. Transcriptome was obtained by RNA sequencing of HSCs and spatialomics with MC38 liver metastases. TGFβ1-induced CAF activation of HSCs was accompanied by elevation of PM Glut1, glucose uptake, and glycolysis. Targeting Glut1 or Src by short hairpin RNA, pharmacologic inhibition, or a Src SH3 domain deletion mutant abrogated TGFβ1-stimulated PM accumulation of Glut1, glycolysis, and CAF activation. Mechanistically, binding of the Src SH3 domain to SH3 domain–binding protein 5 led to a Src/SH3 domain–binding protein 5/Rab11/Glut1 complex that activated Rab11-dependent Glut1 PM transport under TGFβ1 stimulation. Deleting the Src SH3 domain or targeting Glut1 of HSCs by short hairpin RNA or Cre recombinase/LoxP sequence derived from bacteriophage P1 recombination suppressed CAF activation in mice and MC38 colorectal liver metastasis. Multi-omics revealed that Glut1 deficiency in HSCs/CAFs suppressed HSC expression of tumor-promoting factors and altered MC38 transcriptome, contributing to reduced MC38 liver metastases. Conclusion: The Src SH3 domain–facilitated metabolic reprogramming induced by TGFβ1 represents a target to inhibit CAF activation and the pro-metastatic liver microenvironment.

Abstract Background The development of Crohn’s disease [CD] involves immune cell signalling pathways regulated by epigenetic modifications. Aberrant DNA methylation has been identified in peripheral blood and bulk intestinal tissue from CD patients. However, the DNA methylome of disease-associated intestinal CD4+ lymphocytes has not been evaluated. Materials and Methods Genome-wide DNA methylation sequencing was performed from terminal ileum CD4+ cells from 21 CD patients and 12 age- and sex-matched controls. Data were analysed for differentially methylated CpGs [DMCs] and methylated regions [DMRs]. Integration was performed with RNA-sequencing data to evaluate the functional impact of DNA methylation changes on gene expression. DMRs were overlapped with regions of differentially open chromatin [by ATAC-seq] and CCCTC-binding factor [CTCF] binding sites [by ChIP-seq] between peripherally derived Th17 and Treg cells. Results CD4+ cells in CD patients had significantly increased DNA methylation compared to those from the controls. A total of 119 051 DMCs and 8113 DMRs were detected. While hypermethylated genes were mostly related to cell metabolism and homeostasis, hypomethylated genes were significantly enriched within the Th17 signalling pathway. The differentially enriched ATAC regions in Th17 cells [compared to Tregs] were hypomethylated in CD patients, suggesting heightened Th17 activity. There was significant overlap between hypomethylated DNA regions and CTCF-associated binding sites. Conclusions The methylome of CD patients shows an overall dominant hypermethylation yet hypomethylation is more concentrated in proinflammatory pathways, including Th17 differentiation. Hypomethylation of Th17-related genes associated with areas of open chromatin and CTCF binding sites constitutes a hallmark of CD-associated intestinal CD4+ cells.

Abstract Understanding the function of non-coding genomic sequence variants represents a challenge for biomedicine. Many diseases are products of gene-by-environment interactions with complex mechanisms. This study addresses these themes by mechanistic characterization of non-coding variants that influence gene expression only after drug or hormone exposure. Using glucocorticoid signaling as a model system, we integrated genomic, transcriptomic, and epigenomic approaches to unravel mechanisms by which variant function could be revealed by hormones or drugs. Specifically, we identified cis-regulatory elements and 3D interactions underlying ligand-dependent associations between variants and gene expression. One-quarter of the glucocorticoid-modulated variants that we identified had already been associated with clinical phenotypes. However, their affected genes were ‘unmasked’ only after glucocorticoid exposure and often with function relevant to the disease phenotypes. These diseases involved glucocorticoids as risk factors or therapeutic agents and included autoimmunity, metabolic and mood disorders, osteoporosis and cancer. For example, we identified a novel breast cancer risk gene, MAST4, with expression that was repressed by glucocorticoids in cells carrying the risk genotype, repression that correlated with MAST4 expression in breast cancer and treatment outcomes. These observations provide a mechanistic framework for understanding non-coding genetic variant-chemical environment interactions and their role in disease risk and drug response.

AbstractAlternative polyadenylation (APA) is a common genetic regulatory mechanism that generates distinct 3′ ends for RNA transcripts. Changes in APA have been associated with multiple biological processes and disease phenotypes. However, the role of hormones and their drug analogs in APA remains largely unknown. In this study, we investigated transcriptome‐wide the impact of glucocorticoids on APA in 30 human B‐lymphoblastoid cell lines. We found that glucocorticoids could regulate APA for a subset of genes, possibly by changing the expression of 142 RNA‐binding proteins, some with known APA‐regulating properties. Interestingly, genes with glucocorticoid‐mediated APA were enriched in viral translation‐related pathways, while genes with glucocorticoid‐mediated expression were enriched in interferon and interleukin pathways, suggesting that glucocorticoid‐mediated APA might result in functional consequences distinct from gene expression. For example, glucocorticoids, a pharmacotherapy for severe COVID‐19, were found to change the APA but not the expression of LY6E, an important antiviral inhibitor in coronavirus diseases. Glucocorticoid‐mediated APA was also cell‐type‐specific, suggesting an action of glucocorticoids that may be unique to immune regulation. We also observed evidence for genotype‐dependent glucocorticoid‐mediated APA (referred to as pharmacogenomic‐alterative polyadenylation quantitative trait loci), providing potential functional mechanisms for a series of common genetic variants that had previously been associated with immune disorders, but without a clear mechanism. In summary, this study reports a series of observations regarding the impact of glucocorticoids on APA, raising the possibility that this mechanism might have implications for both disease pathophysiology and drug therapy.

Matrix stiffness is a central regulator of fibroblast function. However, the transcriptional mechanisms linking matrix stiffness to changes in fibroblast phenotype are incompletely understood. Here, we evaluated the effect of matrix stiffness on genome-wide chromatin accessibility in freshly isolated lung fibroblasts using ATAC-seq. We found higher matrix stiffness profoundly increased global chromatin accessibility relative to lower matrix stiffness, and these alterations were in close genomic proximity to known profibrotic gene programs. Motif analysis of these regulated genomic loci identified ZNF416 as a putative mediator of fibroblast stiffness responses. Genome occupancy analysis using ChIP-seq confirmed that ZNF416 occupies a broad range of genes implicated in fibroblast activation and tissue fibrosis, with relatively little overlap in genomic occupancy with other mechanoresponsive and profibrotic transcriptional regulators. Using loss- and gain-of-function studies, we demonstrated that ZNF416 plays a critical role in fibroblast proliferation, extracellular matrix synthesis, and contractile function. Together, these observations identify ZNF416 as novel mechano-activated transcriptional regulator of fibroblast biology.

Abstract Adenosquamous cancer of the pancreas (ASCP) is a subtype of pancreatic cancer that has a worse prognosis and greater metastatic potential than the more common pancreatic ductal adenocarcinoma (PDAC) subtype. To distinguish the genomic landscape of ASCP and identify actionable targets for this lethal cancer, we applied DNA content flow cytometry to a series of 15 tumor samples including five patient-derived xenografts (PDX). We interrogated purified sorted tumor fractions from these samples with whole-genome copy-number variant (CNV), whole-exome sequencing, and Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) analyses. These identified a variety of somatic genomic lesions targeting chromatin regulators in ASCP genomes that were superimposed on well-characterized genomic lesions including mutations in TP53 (87%) and KRAS (73%), amplification of MYC (47%), and homozygous deletion of CDKN2A (40%) that are common in PDACs. Furthermore, a comparison of ATAC-seq profiles of three ASCP and three PDAC genomes using flow-sorted PDX models identified genes with accessible chromatin unique to the ASCP genomes, including the lysine methyltransferase SMYD2 and the pancreatic cancer stem cell regulator RORC in all three ASCPs, and a FGFR1-ERLIN2 fusion associated with focal CNVs in both genes in a single ASCP. Finally, we demonstrate significant activity of a pan FGFR inhibitor against organoids derived from the FGFR1-ERLIN2 fusion–positive ASCP PDX model. Our results suggest that the genomic and epigenomic landscape of ASCP provide new strategies for targeting this aggressive subtype of pancreatic cancer. Significance: These data provide a unique description of the ASCP genomic and epigenomic landscape and identify candidate therapeutic targets for this dismal cancer.

Abstract Triple Negative Breast Cancer (TNBC) affects approximately 15-20% of BC patients, yet accounts for a disproportionately higher rate of BC morbidity and mortality, in part due to lack of targeted therapies. We have shown that estrogen receptor beta (ERβ) is expressed in approximately 20% of TN breast tumors and that ligand-mediated activation of ERβ with estradiol (E2) or ERβ-selective agonists decreases tumor cell proliferation, invasion and migration in vitro and in vivo. Therefore, we aimed to elucidate the mechanisms by which ERβ elicits its anti-cancer effects in TNBC. RNAseq analysis of ERβ-expressing MDA-MB-231 cells demonstrated that ERβ significantly downregulates NFκB signaling in the presence of E2. ChIPseq for ERβ in these cells revealed that ERβ primarily associated with estrogen response elements, but 12% of all ERβ binding sites were located at NFκB consensus motifs. Using an NFκB reporter construct and qPCR, ERβ was shown to block TNFα-mediated induction of NFκB signaling and NFκB target gene expression. RNAseq analysis of MDA-MB-231-ERβ cells treated with TNFα or E2+TNFα revealed substantial global inhibition of TNFα regulated genes in the presence of E2. ChIPseq for NFκB demonstrated that ERβ significantly alters NFκB’s cistrome whereby it can both diminish NFκB binding and redistribute NFκB throughout the genome. ChIPseq also demonstrated that ligand-mediated activation of ERβ significantly diminished an activating histone mark (H3K27Ac) at many of these NFκB target genes while enhancing a repressive mark (H3K27Me3). The addition of H3K27Me3 at these loci was shown to occur through the recruitment of the histone methyltransferase, EZH2. Drug-mediated blockade of EZH2 activity reversed suppression of NFκB target gene expression by ERβ. Knockdown of NFκB or Mutation of ERβ’s DNA binding domain rendered ERβ incapable of associating with DNA, recruiting EZH2, methylating NFκB target gene loci, repressing NFκB target gene expression and inhibiting proliferation. Interestingly, ERβ was shown to elicit more potent anti-cancer effects in TNBC cells expression a constitutively active form of NFκB. These finding suggest that a primary mechanism by which ERβ functions as a tumor suppressor is through inhibition of NFκB pathway activity. Our studies have also revealed that ERβ functions as a molecular switch for EZH2 and repurposes it for tumor suppressive activities, as EZH2 has previously been reported to enrich NFκB signaling in TNBC. These findings could address the paradox that high EZH2 expression is associated with worse TNBC patient outcomes, while high H3K27Me3 expression is associated with improved patient outcomes. Currently, a Mayo Clinic Breast Cancer SPORE prospective phase II clinical trial is underway to investigate the efficacy of estradiol for the treatment of metastatic ERβ+TNBC and to further evaluate the cross-talk between ERβ, EZH2 and NFκB signaling.

Abstract Background: Triple Negative Breast Cancer (TNBC) affects approximately 15-20% of BC patients, yet accounts for a disproportionately higher rate of BC morbidity and mortality, in part due to lack of targeted therapies. Using well-validated antibodies, we have found that estrogen receptor beta (ERβ) protein is expressed in approximately 20% of TN breast tumors. We have demonstrated that ligand-mediated activation of ERβ with estradiol (E2) or ERβ-selective agonists decreases tumor cell proliferation, invasion, and migration in vitro, as well as primary tumor growth and metastatic spread in vivo. As such, we aimed to elucidate the mechanisms by which ERβ elicits its anti-cancer effects in TNBC. Methods: RNAseq and ChIPseq were used to identify ERβ and NFκB regulomes and cistromes in TNBC cells. The global effects of ERβ on H3K27Ac and H3K27Me3 were also assessed via ChIPseq. The influence of ERβ on NFκB signaling was determined using co-immunoprecipitation, luciferase assays transcriptional readouts and ChIPseq. The involvement of EZH2 in mediating ERβ’s suppression of NFκB signaling and cancer cell phenotypes was determined using pharmacological approaches. Results: Pathway analysis of ERβ-regulated genes identified NFκB and inflammatory signaling as being highly suppressed in response to E2 treatment. Numerous NFκB target genes were among the most down-regulated genes by ERβ. ChIPseq revealed that ERβ primarily associated with estrogen response elements (EREs), but was also highly enriched at NFκB binding sites following E2 treatment. However, ERβ had no impact on TNFα-mediated NFκB phosphorylation or nuclear localization, but was shown to physically associate with NFκB protein. Using an NFκB reporter construct and qPCR, ERβ was shown to block TNFα-mediated induction of NFκB signaling and NFκB target gene expression. Globally, RNAseq identified 373 genes to be significantly regulated by TNFα in TNBC cells, of which 107 were blocked in the presence of E2+TNFα. Interestingly, E2+TNFα treatment induced novel NFκB binding sites compared to TNFα alone as identified by ChIPseq for NFκB. ChIPseq also demonstrated that ligand-mediated activation of ERβ significantly diminished an activating histone mark (H3K27Ac) at many of these NFκB target genes while enhancing a repressive mark (H3K27Me3). These modifications were shown to occur through the recruitment of the histone methyltransferase, EZH2, to enhancer elements of these NFκB target genes. Drug-mediated blockade of EZH2 activity reversed suppression of NFκB target gene expression by ERβ. Finally, knockdown of NFκB expression inhibited the anti-cancer effects of ERβ while expression of a constitutively active form of NFκB enhanced ERβ’s inhibitory effects. Conclusions: Our data suggest that ERβ elicits its anti-cancer effects, in part, via formation of a novel co-repressor complex consisting of ERβ, NFκB, and EZH2. These data are in keeping with prior observations of the importance of NFκB signaling as it relates to TNBC cell proliferation and invasion, and that decreased expression of NFκB target genes is associated with improved outcomes in TNBC patients. Currently, a Mayo Breast SPORE prospective study is underway to investigate the efficacy of estradiol for the treatment of metastatic ERβ+TNBC and to further evaluate the cross-talk between ERβ and NFκB signaling. Citation Format: Kirsten GM Aspros, Adam W Nelson, Zhenqing Ye, Zhifu Sun, Igor Chernukhin, Jason S Carroll, James N Ingle, Matthew P Goetz, John R Hawse. Estrogen receptor beta repurposes EZH2 to inhibit oncogenic NFκB signaling in triple negative breast cancer [abstract]. In: Proceedings of the 2019 San Antonio Breast Cancer Symposium; 2019 Dec 10-14; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2020;80(4 Suppl):Abstract nr P6-05-02.

Abstract Motivation DNA methylation can be measured at the single CpG level using sodium bisulfite conversion of genomic DNA followed by sequencing or array hybridization. Many analytic tools have been developed, yet there is still a high demand for a comprehensive and multifaceted tool suite to analyze, annotate, QC and visualize the DNA methylation data. Results We developed the CpGtools package to analyze DNA methylation data generated from bisulfite sequencing or Illumina methylation arrays. The CpGtools package consists of three types of modules: (i) ‘CpG position modules’ focus on analyzing the genomic positions of CpGs, including associating other genomic and epigenomic features to a given list of CpGs and generating the DNA motif logo enriched in the genomic contexts of a given list of CpGs; (ii) ‘CpG signal modules’ are designed to analyze DNA methylation values, such as performing the PCA or t-SNE analyses, using Bayesian Gaussian mixture modeling to classify CpG sites into fully methylated, partially methylated and unmethylated groups, profiling the average DNA methylation level over user-specified genomics regions and generating the bean/violin plots and (iii) ‘differential CpG analysis modules’ focus on identifying differentially methylated CpGs between groups using different statistical methods including Fisher’s Exact Test, Student’s t-test, ANOVA, non-parametric tests, linear regression, logistic regression, beta-binomial regression and Bayesian estimation. Availability and implementation CpGtools is written in Python under the open-source GPL license. The source code and documentation are freely available at https://github.com/liguowang/cpgtools. Supplementary information Supplementary data are available at Bioinformatics online.