Aging Human Hematopoietic Stem Cells Manifest Massive Epigenetic Reprogramming and Altered Gene Splicing of Key Hematopoietic Gene Sets

Blood ◽  
2016 ◽  
Vol 128 (22) ◽  
pp. 885-885
Author(s):  
Emmalee R Adelman ◽  
Jian Shi ◽  
Maria E Figueroa

Abstract Aging leads to impairment of hematopoietic stem cell (HSC) function with decreased self-renewal, imbalanced differentiation potential and an increased risk to develop myeloid malignancies. These malignancies are associated with epigenetic deregulation, which contributes to pathogenesis. Notably, studies in murine models have revealed epigenetic changes in aged HSC. However, it is unknown if this occurs in normal human HSC aging and whether it may contribute to HSC dysfunction. Therefore, we performed comprehensive epigenomic and transcriptional profiling in primary human HSC (Lin-, CD34+, CD38-) isolated from young (18-30 yo), mid (45-55 yo) and old (65-75 yo) healthy donors. Using a micro-ChIP-seq protocol we profiled H3K4me1, H3K4me3, H3K27me3 and H3K27ac in 4-7 donors per age group, as well as genome-wide DNA methylation (5mC), hydroxymethylation (5hmC) and RNA-seq. Analysis of enhancer-associated marks revealed that with age there is marked reduction in both H3K4me1 and H3K27ac (20,783 and 15,625 peaks lost, respectively; log10likelihood ratio >3). Gene ontology analysis of these lost peaks revealed their association with genes involved in hematopoiesis and, RNA splicing and chromatin organization, respectively (ChIPenrich, FDR<0.05). In addition, regions depleted in H3K4me1 are enriched for PU.1, FLI1, ETS, and CTCF binding sites (Homer, q<0.00001). We next asked if aging results in specific remodeling of poised (H3K4me1>H3K4me3, H3K27ac-) and active (H3K4me1>H3K4me3, H3K27ac+) enhancers. We found age-related loss of H3K4me1 enrichment at 10,696 poised enhancers, which are associated with hematopoiesis and T- and B-cell receptor signaling (FDR<0.05). We also identified 17,242 active enhancers in young HSC, 7,057 of which are depleted in old HSC. This loss of active enhancers targets genes associated with hematopoiesis, immune signaling and myeloid malignancies (FDR<0.05). Next we analyzed the impact of aging on promoter-associated marks, H3K4me3 and H3K27me3. Remarkably, while aging leads to loss of 22,689 H3K4me3 peaks, only 1,339 H3K27me3 peaks are lost. Loss of H3K4me3 targets genes involved in inflammatory response, development and WNT signaling (FDR<0.05). Given this uneven change in H3K4me3/H3K27me3 with aging, we hypothesized this may correlate with changes in bivalently marked promoters, which regulate key developmental genes. Out of 3,947 bivalent promoters in young HSC, 842 are lost in aged HSC. This loss of bivalency affects genes involved in WNT, Cadherin and Hedgehog signaling pathways (FDR<0.05). Next we analyzed changes in cytosine modifications. We observe widespread gain of 5hmC (n=14,554 differentially hydroxymethylated regions [DHMR]; FDR <0.005), with specific enrichment at introns and exons (p<2.2e-16), as well as enrichment for GATA and KLF binding sites (Homer, q<0.00001). These DHMR target genes involved in hematopoiesis, proteins regulated by alternative splicing, and pathways associated with cancer (FDR<0.05). In contrast, much more subtle changes are found in 5mC with HSC aging, with only 529 differentially methylated regions (q-value <0.05, meth.diff >20%). However, these subtle changes also target genes associated with cadherin and WNT signaling. Finally, RNA-seq analysis revealed that this age-associated epigenetic reprogramming is accompanied by an overall trend to gene downregulation. Amongst the genes most affected are the nuclear lamin gene LMNA (mutated in progeria syndrome), splicing factors SRSF7 and U2AF1 and, the transcription factors KLF3/6 and HIF1α (FDR <0.05, fold change >1.5). Notably, changes in expression also include significant differential exon usage, which may be mediated by DHMR at intron-exon boundaries: 575 genes show altered exon usage (FDR <0.05, fold change>1.5) including LMNA and the epigenetic modifiers BRD9, CITED2, KDM6A and SETD6. In summary, we have completed the first comprehensive epigenomic profiling of aging in human HSC. Our findings show massive epigenetic remodeling in aged HSC, consisting of loss of activating histone marks primarily targeting enhancers and bivalent promoters at genes involved in hematopoiesis and developmental pathways. Cytosine modifications show widespread changes in 5hmC, targeting intron-exon boundaries. Globally, this epigenetic reprogramming results in overall gene downregulation and altered splicing of genes important for HSC regulation. Disclosures No relevant conflicts of interest to declare.

Blood ◽  
2015 ◽  
Vol 126 (23) ◽  
pp. 302-302
Author(s):  
Jean-Baptiste Micol ◽  
Nicolas Duployez ◽  
Alessandro Pastore ◽  
Robert Williams ◽  
Eunhee Kim ◽  
...  

Abstract Mutations in Addition of Sex Combs Like 1 (ASXL1) are common in patients with myeloid leukemias. More recently, mutations in ASXL2, a paralog of ASXL1 with ~40% shared amino acid homology, have been discovered to occur specifically in patients with acute myeloid leukemia (AML) patients bearing the RUNX1-ETO (AML1-ETO; RUNX1-RUNX1T1) translocation and are amongst the most common mutations in RUNX1-ETO AML (mutated in 20-25% of patients). Although ASXL1 is critical for Polycomb Repressive Complex 2 function in myeloid hematopoietic cells and loss of Asxl1 recapitulates key aspects of myelodysplastic syndrome (MDS), the function of ASXL2 in normal or malignant hematopoiesis is unknown. We therefore set out to perform a functional comparison of ASXL1and ASXL2on hematopoiesis and transcription and determine the functional basis for frequent mutations in RUNX1-ETO AML. In vitro analyses of ASXL2 insertion/deletion mutations revealed that these mutations resulted in substantial reduction of ASXL2 protein expression, stability, and half-life. We therefore generated Asxl2 conditional knockout (cKO) mice to delineate the effect of ASXL2 loss on hematopoiesis. Competitive (Fig. 1A) and noncompetitive transplantation revealed that Asxl2 or compound Asxl1/2 loss resulted in cell-autonomous, rapid defects of hematopoietic stem cell function, self-renewal, and number with peripheral blood leukopenia and thrombocytopenia but without any obvious MDS features- phenotypes distinct from Asxl1 cKO mice. Mice with heterozygous deletion of Asxl2 demonstrated an intermediate phenotype between control and homozygous cKO mice indicating a gene dosage effect of Asxl2 loss. RNA sequencing (RNA-seq) of hematopoietic stem/progenitor cells from Asxl2- and Asxl1-deficient mice revealed twenty-fold greater differentially expressed genes in Asxl2 cKO mice relative to Asxl1 cKO mice. Interestingly, genes differentially expressed with Asxl2 loss significantly overlapped with direct transcriptional targets of RUNX1-ETO, findings not seen in Asxl1 cKO mice (Fig. 1B). Asxl2 target genes appeared to also be targets of RUNX1, a key gene repressed by RUNX1-ETO to promote leukemogenesis. Consistent with this, genome-wide analysis of Asxl2 binding sites through anti-Asxl2 ChIP-seq revealed that Asxl2 binding sites substantially overlap with those of Runx1. Overall, the above data suggest that Asxl2 may be a critical mediator of RUNX1-ETO mediated leukemogenesis by affecting the expression of RUNX1 and/or RUNX1-ETO target genes. RNA-seq of primary RUNX1-ETO AML patient samples revealed that ASXL2-mutant RUNX1-ETO patients form a distinct transcriptional subset of RUNX1-ETO AML (Fig. 1C) suggesting a specific role of ASXL2 in leukemogenesis. To functionally interrogate the role of ASXL2 loss in RUNX1-ETO mediated leukemogenesis we first utilized an in vitro model with RNAi-mediated depletion of ASXL1 or ASXL2 in the SKNO1 cell line (the only ASXL-wildtype human RUNX1-ETO cell line). RNA-seq revealed distinct target genes dysregulated by ASXL1 versus ASXL2 loss in these cells without any significant overlap. Anti-ASXL2, RUNX1, and RUNX1-ETO ChIPSeq in SKNO1 cells revealed significant co-occupancy of ASXL2 with RUNX1 and RUNX1-ETO binding sites. Moreover, analysis of histone modification ChIPSeq revealed an enrichment in intergenic and enhancer H3K4me1 abundance following ASXL2 loss in SKNO1 cells. Next, to understand the in vivo effects of Asxl2 loss in the context of RUNX1-ETO, we performed retroviral bone marrow (BM) transplantation assays using RUNX1-ETO9a in Asxl2 cKO mice. In contrast to the failure of hematopoietic stem cell function with Asxl2 deletion alone, mice reconstituted with BM cells expressing RUNX1-ETO9a in Asxl2-deficient background had a shortened leukemia-free survival compared to Asxl2 -wildtype control. Overall, these data reveal that ASXL2 is required for hematopoiesis and has differing biological and transcriptional functions from ASXL1. Moreover, this work identifies ASXL2 as a novel mediator of RUNX1-ETOtranscriptional function and provides a new model of penetrant RUNX1-ETO AML based on genetic events found in a substantial proportion of t(8;21) AML patients. Further interrogation of the enhancer alterations generated by ASXL2 loss in RUNX1-ETO AML may highlight new therapeutic approaches for this subset of AML. Figure 1. Figure 1. Disclosures No relevant conflicts of interest to declare.


Hypertension ◽  
2012 ◽  
Vol 60 (suppl_1) ◽  
Author(s):  
Catherina A Cuevas ◽  
Alexis A Gonzalez ◽  
Nivaldo C Inestrosa ◽  
Carlos P Vio ◽  
Minolfa C Prieto

The prorenin receptor (PRR) is upregulated in the kidney by high angiotensin II (Ang II) states such as those that occur with AngII-dependent hypertension and low salt diet. The PRR is an accessory protein of the vacuolar H-ATPase, which facilitates Wnt/β-catenin signaling. The Wnt/β-catenin pathway is involved in fibrosis processes. In the present study, we aimed to determine whether the stimulation of PRR in mouse collecting duct M-1 cells induces fibrotic genes independently of Ang II, and if this effect is mediated by activation of Wnt/β-catenin. Both Ang II (10 -7 M) and human recombinant prorenin (hRPr; 2,5 x 10 -8 M) treatments (8 and 16 hours) increased mRNA and protein levels of fibronectin and collagen I (1.5±0.08 and 1.5 ± 0.1 fold change, respectibely; p<0.05); however, the effects of hRPr were elicited earlier. Likewise, Ang II and hRPr stimulated the Wnt target genes, cyclin D1 and c-myc (cyclin D1: 2±0.2 for both; c-myc: 1.4 ± 0.03 and 1.2± 0.002 fold change for Ang II and hRPr, respectively; p<0.001). Ang II type 1 receptor (AT1R) blockade with candesartan (10 -7 M) completely prevented the Ang II-dependent stimulation but not the effects of hRPr on Wnt signaling genes. Upregulation of fibronectin and collagen I genes by Ang II or hRP at 16 h was prevented by Wnt signaling inhibition with Pyrvinium Pamoate (10 -7 M). The data indicate that in M-1 cells, activation of AT1R and PRR stimulate the synthesis of fibrotic genes via Wnt signaling by independent mechanisms.


2020 ◽  
Vol 48 (19) ◽  
pp. 10768-10784
Author(s):  
Ge Sun ◽  
Chunyu Wang ◽  
Shengli Wang ◽  
Hongmiao Sun ◽  
Kai Zeng ◽  
...  

Abstract Estrogen receptor alpha (ERα) signaling pathway is essential for ERα-positive breast cancer progression and endocrine therapy resistance. Bromodomain PHD Finger Transcription Factor (BPTF) associated protein of 18kDa (BAP18) has been recognized as a crucial H3K4me3 reader. However, the whole genomic occupation of BAP18 and its biological function in breast cancer is still elusive. Here, we found that higher expression of BAP18 in ERα-positive breast cancer is positively correlated with poor prognosis. ChIP-seq analysis further demonstrated that the half estrogen response elements (EREs) and the CCCTC binding factor (CTCF) binding sites are the significant enrichment sites found in estrogen-induced BAP18 binding sites. Also, we provide the evidence to demonstrate that BAP18 as a novel co-activator of ERα is required for the recruitment of COMPASS-like core subunits to the cis-regulatory element of ERα target genes in breast cancer cells. BAP18 is recruited to the promoter regions of estrogen-induced genes, accompanied with the enrichment of the lysine 4-trimethylated histone H3 tail (H3K4me3) in the presence of E2. Furthermore, BAP18 promotes cell growth and associates the sensitivity of antiestrogen in ERα-positive breast cancer. Our data suggest that BAP18 facilitates the association between ERα and COMPASS-like core subunits, which might be an essential epigenetic therapeutic target for breast cancer.


2021 ◽  
Vol 12 ◽  
Author(s):  
Diana García-Cortés ◽  
Enrique Hernández-Lemus ◽  
Jesús Espinal-Enríquez

Luminal A is the most common breast cancer molecular subtype in women worldwide. These tumors have characteristic yet heterogeneous alterations at the genomic and transcriptomic level. Gene co-expression networks (GCNs) have contributed to better characterize the cancerous phenotype. We have previously shown an imbalance in the proportion of intra-chromosomal (cis-) over inter-chromosomal (trans-) interactions when comparing cancer and healthy tissue GCNs. In particular, for breast cancer molecular subtypes (Luminal A included), the majority of high co-expression interactions connect gene-pairs in the same chromosome, a phenomenon that we have called loss of trans- co-expression. Despite this phenomenon has been described, the functional implication of this specific network topology has not been studied yet. To understand the biological role that communities of co-expressed genes may have, we constructed GCNs for healthy and Luminal A phenotypes. Network modules were obtained based on their connectivity patterns and they were classified according to their chromosomal homophily (proportion of cis-/trans- interactions). A functional overrepresentation analysis was performed on communities in both networks to observe the significantly enriched processes for each community. We also investigated possible mechanisms for which the loss of trans- co-expression emerges in cancer GCN. To this end we evaluated transcription factor binding sites, CTCF binding sites, differential gene expression and copy number alterations (CNAs) in the cancer GCN. We found that trans- communities in Luminal A present more significantly enriched categories than cis- ones. Processes, such as angiogenesis, cell proliferation, or cell adhesion were found in trans- modules. The differential expression analysis showed that FOXM1, CENPA, and CIITA transcription factors, exert a major regulatory role on their communities by regulating expression of their target genes in other chromosomes. Finally, identification of CNAs, displayed a high enrichment of deletion peaks in cis- communities. With this approach, we demonstrate that network topology determine, to at certain extent, the function in Luminal A breast cancer network. Furthermore, several mechanisms seem to be acting together to avoid trans- co-expression. Since this phenomenon has been observed in other cancer tissues, a remaining question is whether the loss of long distance co-expression is a novel hallmark of cancer.


2019 ◽  
Author(s):  
Martin Lewinski ◽  
Yannik Bramkamp ◽  
Tino Köster ◽  
Dorothee Staiger

AbstractBackgroundRNA-binding proteins interact with their target RNAs at specific sites. These binding sites can be determined genome-wide through individual nucleotide resolution crosslinking immunoprecipitation (iCLIP). Subsequently, the binding sites have to be visualized. So far, no visualization tool exists that is easily accessible but also supports restricted access so that data can be shared among collaborators.ResultsHere we present SEQing, a customizable interactive dashboard to visualize crosslink sites on target genes of RNA-binding proteins that have been obtained by iCLIP. Moreover, SEQing supports RNA-seq data that can be displayed in a diffrerent window tab. This allows, e.g. crossreferencing the iCLIP data with genes differentially expressed in mutants of the RBP and thus obtain some insights into a potential functional relevance of the binding sites. Additionally, detailed information on the target genes can be incorporated in another tab.ConclusionSEQing is written in Python3 and runs on Linux. The web-based access makes iCLIP data easily accessible, even with mobile devices. SEQing is customizable in many ways and has also the option to be secured by a password. The source code is available at https://github.com/malewins/SEQing.


2020 ◽  
Vol 16 (12) ◽  
pp. e1008488
Author(s):  
Mirko Ronzio ◽  
Andrea Bernardini ◽  
Giulio Pavesi ◽  
Roberto Mantovani ◽  
Diletta Dolfini

NF-Y is a trimeric Transcription Factor -TF- which binds with high selectivity to the conserved CCAAT element. Individual ChIP-seq analysis as well as ENCODE have progressively identified locations shared by other TFs. Here, we have analyzed data introduced by ENCODE over the last five years in K562, HeLa-S3 and GM12878, including several chromatin features, as well RNA-seq profiling of HeLa cells after NF-Y inactivation. We double the number of sequence-specific TFs and co-factors reported. We catalogue them in 4 classes based on co-association criteria, infer target genes categorizations, identify positional bias of binding sites and gene expression changes. Larger and novel co-associations emerge, specifically concerning subunits of repressive complexes as well as RNA-binding proteins. On the one hand, these data better define NF-Y association with single members of major classes of TFs, on the other, they suggest that it might have a wider role in the control of mRNA production.


Blood ◽  
2012 ◽  
Vol 120 (21) ◽  
pp. SCI-30-SCI-30
Author(s):  
Jay L. Hess ◽  
Cailin Collins ◽  
Joel Bronstein ◽  
Yuqing Sun ◽  
Surya Nagaraja

Abstract Abstract SCI-30 HOXA9 plays important roles in both development and hematopoiesis and is overexpressed in more than 50 percent of acute myeloid leukemias (AML). Nearly all cases of AML with mixed lineage leukemia (MLL) translocations show increased HOXA9 expression, as well as cases with mutation of the nucleophosmin gene NPM1, overexpression of CDX2, and fusions of NUP98. In most cases, upregulation of HOXA9 is accompanied by upregulation of its homeodomain-containing cofactor MEIS1, which directly interacts with HOXA9. While HOXA9 alone is sufficient for transformation of hematopoietic stem cells in culture, the addition of MEIS1 increases the transformation efficiency and results in rapidly fatal leukemias in transplanted animals. Despite the crucial role that HOXA9 plays in development, hematopoiesis, and leukemia, its transcriptional targets and mechanisms of action are poorly understood. We have used ChIP-seq to identify Hoxa9 and Meis1 binding sites on a genome-wide level in myeloblastic cells, profiled their associated epigenetic modifications, identified the target genes regulated by HOXA9 and identified HOXA9 interacting proteins. HOXA9 and MEIS1 cobind at hundreds of promoter distal, highly evolutionarily conserved sites showing high levels of histone H3K4 monomethylation and CBP/P300 binding. These include many proleukemogenic gene loci, such as Erg, Flt3, Myb, Lmo2, and Sox4. In addition, HOXA9 binding sites overlap a subset of enhancers previously implicated in myeloid differentiation and inflammation. HOXA9 binding at enhancers stabilizes association of MEIS1 and lineage-restricted transcription factors, including C/EBPα, PU.1, and STAT5A/B thereby promoting CBP/p300 recruitment, histone acetylation, and transcriptional activation. Current efforts are focused on using both biochemical and genetic approaches to assess the role of HOXA9 “enhanceosome” components C/EBPα, PU.1, and STAT5A/B in transcriptional regulation and leukemogenesis. Studies to date suggest that C/EBPα and PU.1 binding can occur in the absence of HOXA9/MEIS1, supporting a model in which these proteins act as pioneer transcription factors for establishment of poised, but not activated, HOXA9-regulated enhancers. Work is under way to assess the impact of high-level HOXA9 and MEIS1 on enhanceosome assembly and the role of recruitment of transcriptional coactivators involved in target gene up- or downregulation, including histone acetyltransferases and chromatin remodeling complexes. Collectively, our findings suggest that HOXA9-regulated enhancers are a fundamental mechanism of HOX-mediated transcription in normal development that is deregulated in leukemia. Disclosures: No relevant conflicts of interest to declare.


2020 ◽  
Vol 21 (8) ◽  
pp. 2800 ◽  
Author(s):  
Xi Wu ◽  
Yang Yang ◽  
Chaoyue Zhong ◽  
Yin Guo ◽  
Tengyu Wei ◽  
...  

Chromatin structure plays a pivotal role in maintaining the precise regulation of gene expression. Accessible chromatin regions act as the binding sites of transcription factors (TFs) and cis-elements. Therefore, information from these open regions will enhance our understanding of the relationship between TF binding, chromatin status and the regulation of gene expression. We employed an assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) and RNA-seq analyses in the gonads of protogynous hermaphroditic orange-spotted groupers during sex reversal to profile open chromatin regions and TF binding sites. We focused on several crucial TFs, including ZNF263, SPIB, and KLF9, and analyzed the networks of TF-target genes. We identified numerous transcripts exhibiting sex-preferred expression among their target genes, along with their associated open chromatin regions. We then investigated the expression patterns of sex-related genes as well as the mRNA localization of certain genes during sex reversal. We found a set of sex-related genes that—upon further study—might be identified as the sex-specific or cell-specific marker genes that trigger sex reversal. Moreover, we discovered the core genes (gnas, ccnb2, and cyp21a) of several pathways related to sex reversal that provide the guideposts for future study.


2016 ◽  
Author(s):  
Naoki Osato

AbstractBackgroundTranscriptional target genes show functional enrichment of genes. However, how many and how significantly transcriptional target genes include functional enrichments are still unclear. To address these issues, I predicted human transcriptional target genes using open chromatin regions, ChIP-seq data and DNA binding sequences of transcription factors in databases, and examined functional enrichment and gene expression level of putative transcriptional target genes.ResultsGene Ontology annotations showed four times larger numbers of functional enrichments in putative transcriptional target genes than gene expression information alone, independent of transcriptional target genes. To compare the number of functional enrichments of putative transcriptional target genes between cells or search conditions, I normalized the number of functional enrichment by calculating its ratios in the total number of transcriptional target genes. With this analysis, native putative transcriptional target genes showed the largest normalized number of functional enrichments, compared with target genes including 5 – 60% of randomly selected genes. The normalized number of functional enrichments was changed according to the criteria of enhancer-promoter interactions such as distance from transcriptional start sites and orientation of CTCF-binding sites. Forward-reverse orientation of CTCF-binding sites showed significantly higher normalized number of functional enrichments than the other orientations. Journal papers showed that the top five frequent functional enrichments were related to the cellular functions in the three cell types. The median expression level of transcriptional target genes changed according to the criteria of enhancer-promoter assignments (i.e. interactions) and was correlated with the changes of the normalized number of functional enrichments of transcriptional target genes.ConclusionsHuman putative transcriptional target genes showed significant functional enrichments. Functional enrichments were related to the cellular functions. The normalized number of functional enrichments of human putative transcriptional target genes changed according to the criteria of enhancer-promoter assignments and correlated with the median expression level of the target genes. These analyses and characters of human putative transcriptional target genes would be useful to examine the criteria of enhancer-promoter assignments and to predict the novel mechanisms and factors such as DNA binding proteins and DNA sequences of enhancer-promoter interactions.


eLife ◽  
2016 ◽  
Vol 5 ◽  
Author(s):  
Luke Isbel ◽  
Lexie Prokopuk ◽  
Haoyu Wu ◽  
Lucia Daxinger ◽  
Harald Oey ◽  
...  

We previously identified Wiz in a mouse screen for epigenetic modifiers. Due to its known association with G9a/GLP, Wiz is generally considered a transcriptional repressor. Here, we provide evidence that it may also function as a transcriptional activator. Wiz levels are high in the brain, but its function and direct targets are unknown. ChIP-seq was performed in adult cerebellum and Wiz peaks were found at promoters and transcription factor CTCF binding sites. RNA-seq in Wiz mutant mice identified genes differentially regulated in adult cerebellum and embryonic brain. In embryonic brain most decreased in expression and included clustered protocadherin genes. These also decreased in adult cerebellum and showed strong Wiz ChIP-seq enrichment. Because a precise pattern of protocadherin gene expression is required for neuronal development, behavioural tests were carried out on mutant mice, revealing an anxiety-like phenotype. This is the first evidence of a role for Wiz in neural function.


Sign in / Sign up

Export Citation Format

Share Document