MYC-induced human acute myeloid leukemia requires a continuing IL-3/GM-CSF costimulus.

Hematopoietic clones with leukemogenic mutations arise in healthy people as they age, but progression to acute myeloid leukemia (AML) is rare. Recent evidence suggests that the microenvironment may play an important role in modulating human AML population dynamics. To investigate this concept further, we examined the combined and separate effects of an oncogene (c-MYC) and exposure to interleukin-3 (IL-3), granulocyte-macrophage colony-stimulating factor (GM-CSF), and stem cell factor (SCF) on the experimental genesis of a human AML in xenografted immunodeficient mice. Initial experiments showed that normal human CD34+ blood cells transduced with a lentiviral MYC vector and then transplanted into immunodeficient mice produced a hierarchically organized, rapidly fatal, and serially transplantable blast population, phenotypically and transcriptionally similar to human AML cells, but only in mice producing IL-3, GM-CSF, and SCF transgenically or in regular mice in which the cells were exposed to IL-3 or GM-CSF delivered using a cotransduction strategy. In their absence, the MYC+ human cells produced a normal repertoire of lymphoid and myeloid progeny in transplanted mice for many months, but, on transfer to secondary mice producing the human cytokines, the MYC+ cells rapidly generated AML. Indistinguishable diseases were also obtained efficiently from both primitive (CD34+CD38-) and late granulocyte-macrophage progenitor (GMP) cells. These findings underscore the critical role that these cytokines can play in activating a malignant state in normally differentiating human hematopoietic cells in which MYC expression has been deregulated. They also introduce a robust experimental model of human leukemogenesis to further elucidate key mechanisms involved and test strategies to suppress them.

Barcoding reveals complex clonal dynamics of de novo transformed human mammary cells.

Most human breast cancers have diversified genomically and biologically by the time they become clinically evident. Early events involved in their genesis and the cellular context in which these events occur have thus been difficult to characterize. Here we present the first formal evidence of the shared and independent ability of basal cells and luminal progenitors, isolated from normal human mammary tissue and transduced with a single oncogene (KRAS(G12D)), to produce serially transplantable, polyclonal, invasive ductal carcinomas within 8 weeks of being introduced either subrenally or subcutaneously into immunodeficient mice. DNA barcoding of the initial cells revealed a dramatic change in the numbers and sizes of clones generated from them within 2 weeks, and the first appearance of many 'new' clones in tumours passaged into secondary recipients. Both primary and secondary tumours were phenotypically heterogeneous and primary tumours were categorized transcriptionally as 'normal-like'. This system challenges previous concepts that carcinogenesis in normal human epithelia is necessarily a slow process requiring the acquisition of multiple driver mutations. It also presents the first description of initial events that accompany the genesis and evolution of malignant human mammary cell populations, thereby contributing new understanding of the rapidity with which heterogeneity in their properties can develop.

Integrative analysis of 111 reference human epigenomes.

The reference human genome sequence set the stage for studies of genetic variation and its association with human disease, but epigenomic studies lack a similar reference. To address this need, the NIH Roadmap Epigenomics Consortium generated the largest collection so far of human epigenomes for primary cells and tissues. Here we describe the integrative analysis of 111 reference human epigenomes generated as part of the programme, profiled for histone modification patterns, DNA accessibility, DNA methylation and RNA expression. We establish global maps of regulatory elements, define regulatory modules of coordinated activity, and their likely activators and repressors. We show that disease- and trait-associated genetic variants are enriched in tissue-specific epigenomic marks, revealing biologically relevant cell types for diverse human traits, and providing a resource for interpreting the molecular basis of human disease. Our results demonstrate the central role of epigenomic information for understanding gene regulation, cellular differentiation and human disease.

Hepatocyte Nuclear Factor 4-Alpha Is Essential for the Active Epigenetic State at Enhancers in Mouse Liver.

Cell-fate determination is influenced by interactions between master transcription factors (TFs) and cis-regulatory elements. Hepatocyte nuclear factor 4 alpha (HNF4A), a liver-enriched TF, acts as a master controller in specification of hepatic progenitor cells by regulating a network of TFs to control onset of hepatocyte cell fate. Using analysis of genome-wide histone modifications, DNA methylation, and hydroxymethylation in mouse hepatocytes, we show that HNF4A occupies active enhancers in hepatocytes and is essential for active histone and DNA signatures, especially acetylation of lysine 27 of histone 3 (H3K27ac) and 5-hydroxymethylcytosine (5hmC). In mice lacking HNF4A protein in hepatocytes, we observed a decrease in both H3K27ac and hydroxymethylation at regions bound by HNF4A. Mechanistically, HNF4A-associated hydroxymethylation (5hmC) requires its interaction with ten-eleven translocation methylcytosine dioxygenase 3 (TET3), a protein responsible for oxidation from 5mC to 5hmC. Furthermore, HNF4A regulates TET3 expression in liver by directly binding to an enhancer region. Conclusion: In conclusion, we identified that HNF4A is required for the active epigenetic state at enhancers that amplifies transcription of genes in hepatocytes.

The prolyl isomerase FKBP25 regulates microtubule polymerization impacting cell cycle progression and genomic stability.

FK506 binding proteins (FKBPs) catalyze the interconversion of cis-trans proline conformers in proteins. Importantly, FK506 drugs have anti-cancer and neuroprotective properties, but the effectors and mechanisms underpinning these properties are not well understood because the cellular function(s) of most FKBP proteins are unclear. FKBP25 is a nuclear prolyl isomerase that interacts directly with nucleic acids and is associated with several DNA/RNA binding proteins. Here, we show the catalytic FKBP domain binds microtubules (MTs) directly to promote their polymerization and stabilize the MT network. Furthermore, FKBP25 associates with the mitotic spindle and regulates entry into mitosis. This interaction is important for mitotic spindle dynamics, as we observe increased chromosome instability in FKBP25 knockdown cells. Finally, we provide evidence that FKBP25 association with chromatin is cell-cycle regulated by Protein Kinase C phosphorylation. This disrupts FKBP25-DNA contacts during mitosis while maintaining its interaction with the spindle apparatus. Collectively, these data support a model where FKBP25 association with chromatin and MTs is carefully choreographed to ensure faithful genome duplication. Additionally, they highlight that FKBP25 is a MT-associated FK506 receptor and potential therapeutic target in MT-associated diseases.

A transgenic mouse model demonstrating the oncogenic role of mutations in the polycomb-group gene EZH2 in lymphomagenesis.

The histone methyltransferase EZH2 is frequently mutated in germinal center-derived diffuse large B-cell lymphoma and follicular lymphoma. To further characterize these EZH2 mutations in lymphomagenesis, we generated a mouse line where EZH2(Y641F) is expressed from a lymphocyte-specific promoter. Spleen cells isolated from the transgenic mice displayed a global increase in trimethylated H3K27, but the mice did not show an increased tendency to develop lymphoma. As EZH2 mutations often coincide with other mutations in lymphoma, we combined the expression of EZH2(Y641F) by crossing these transgenic mice with Eµ-Myc transgenic mice. We observed a dramatic acceleration of lymphoma development in this combination model of Myc and EZH2(Y641F). The lymphomas show histologic features of high-grade disease with a shift toward a more mature B-cell phenotype, increased cycling and gene expression, and epigenetic changes involving important pathways in B-cell regulation and function. Furthermore, they initiate disease in secondary recipients. In summary, EZH2(Y641F) can collaborate with Myc to accelerate lymphomagenesis demonstrating a cooperative role of EZH2 mutations in oncogenesis. This murine lymphoma model provides a new tool to study global changes in the epigenome caused by this frequent mutation and a promising model system for testing novel treatments.

Identification and characterization of Hoxa9 binding sites in hematopoietic cells.

The clustered homeobox proteins play crucial roles in development, hematopoiesis, and leukemia, yet the targets they regulate and their mechanisms of action are poorly understood. Here, we identified the binding sites for Hoxa9 and the Hox cofactor Meis1 on a genome-wide level and profiled their associated epigenetic modifications and transcriptional targets. Hoxa9 and the Hox cofactor Meis1 cobind at hundreds of highly evolutionarily conserved sites, most of which are distant from transcription start sites. These sites show high levels of histone H3K4 monomethylation and CBP/P300 binding characteristic of enhancers. Furthermore, a subset of these sites shows enhancer activity in transient transfection assays. Many Hoxa9 and Meis1 binding sites are also bound by PU.1 and other lineage-restricted transcription factors previously implicated in establishment of myeloid enhancers. Conditional Hoxa9 activation is associated with CBP/P300 recruitment, histone acetylation, and transcriptional activation of a network of proto-oncogenes, including Erg, Flt3, Lmo2, Myb, and Sox4. Collectively, this work suggests that Hoxa9 regulates transcription by interacting with enhancers of genes important for hematopoiesis and leukemia.

Retrotransposon-induced heterochromatin spreading in the mouse revealed by insertional polymorphisms.

The "arms race" relationship between transposable elements (TEs) and their host has promoted a series of epigenetic silencing mechanisms directed against TEs. Retrotransposons, a class of TEs, are often located in repressed regions and are thought to induce heterochromatin formation and spreading. However, direct evidence for TE-induced local heterochromatin in mammals is surprisingly scarce. To examine this phenomenon, we chose two mouse embryonic stem (ES) cell lines that possess insertionally polymorphic retrotransposons (IAP, ETn/MusD, and LINE elements) at specific loci in one cell line but not the other. Employing ChIP-seq data for these cell lines, we show that IAP elements robustly induce H3K9me3 and H4K20me3 marks in flanking genomic DNA. In contrast, such heterochromatin is not induced by LINE copies and only by a minority of polymorphic ETn/MusD copies. DNA methylation is independent of the presence of IAP copies, since it is present in flanking regions of both full and empty sites. Finally, such spreading into genes appears to be rare, since the transcriptional start sites of very few genes are less than one Kb from an IAP. However, the B3galtl gene is subject to transcriptional silencing via IAP-induced heterochromatin. Hence, although rare, IAP-induced local heterochromatin spreading into nearby genes may influence expression and, in turn, host fitness.

Synovial Sarcoma Chromatin Dynamics Reveal a Continuum in SS18:SSX Reprograming.

Synovial sarcoma (SyS) is an aggressive soft-tissue malignancy characterized by a pathognomonic chromosomal translocation leading to the formation of the SS18::SSX fusion oncoprotein. SS18::SSX associates with mammalian BAF complexes suggesting deregulation of chromatin architecture as the oncogenic driver in this tumour type. To examine the epigenomic state of SyS we performed comprehensive multi-omics analysis on 52 primary pre-treatment human SyS tumours. Our analysis revealed a continuum of epigenomic states across the cohort at fusion target genes independent of rare somatic genetic lesions. We identify cell-of-origin signatures defined by enhancer states and reveal unexpected relationships between H2AK119Ub1 and active marks. The number of bivalent promoters, dually marked by the repressive H3K27me3 and activating H3K4me3 marks, has strong prognostic value and outperforms tumor grade in predicting patient outcome. Finally, we identify SyS defining epigenomic features including H3K4me3 expansion associated with striking promoter DNA hypomethylation in which SyS displays the lowest mean methylation level of any sarcoma subtype. We explore these distinctive features as potential vulnerabilities in SyS and identify H3K4me3 inhibition as a promising therapeutic strategy.

RUNX1 colludes with NOTCH1 to reprogram chromatin in T cell acute lymphoblastic leukemia.

Runt-related transcription factor 1 (RUNX1) is oncogenic in diverse types of leukemia and epithelial cancers where its expression is associated with poor prognosis. Current models suggest that RUNX1 cooperates with other oncogenic factors (e.g., NOTCH1, TAL1) to drive the expression of proto-oncogenes in T cell acute lymphoblastic leukemia (T-ALL) but the molecular mechanisms controlled by RUNX1 and its cooperation with other factors remain unclear. Integrative chromatin and transcriptional analysis following inhibition of RUNX1 and NOTCH1 revealed a surprisingly widespread role of RUNX1 in the establishment of global H3K27ac levels and that RUNX1 is required by NOTCH1 for cooperative transcription activation of key NOTCH1 target genes including MYC, DTX1, HES4, IL7R, and NOTCH3. Super-enhancers were preferentially sensitive to RUNX1 knockdown and RUNX1-dependent super-enhancers were disrupted following the treatment of a pan-BET inhibitor, I-BET151.

Meis1 establishes the pre-hemogenic endothelial state prior to Runx1 expression.

Hematopoietic stem and progenitor cells (HSPCs) originate from an endothelial-to-hematopoietic transition (EHT) during embryogenesis. Characterization of early hemogenic endothelial (HE) cells is required to understand what drives hemogenic specification and to accurately define cells capable of undergoing EHT. Using Cellular Indexing of Transcriptomes and Epitopes by Sequencing (CITE-seq), we define the early subpopulation of pre-HE cells based on both surface markers and transcriptomes. We identify the transcription factor Meis1 as an essential regulator of hemogenic cell specification in the embryo prior to Runx1 expression. Meis1 is expressed at the earliest stages of EHT and distinguishes pre-HE cells primed towards the hemogenic trajectory from the arterial endothelial cells that continue towards a vascular fate. Endothelial-specific deletion of Meis1 impairs the formation of functional Runx1-expressing HE which significantly impedes the emergence of pre-HSPC via EHT. Our findings implicate Meis1 in a critical fate-determining step for establishing EHT potential in endothelial cells.

A novel methylated cell-free DNA marker panel to monitor treatment response in metastatic prostate cancer.

Aim: This study examined circulating cell-free DNA (cfDNA) biomarkers associated with androgen treatment resistance in metastatic castration resistance prostate cancer (mCRPC). Materials & methods: We designed a panel of nine candidate cfDNA methylation markers using droplet digital PCR (Methyl-ddPCR) and assessed methylation levels in sequentially collected cfDNA samples from patients with mCRPC. Results: Increased cfDNA methylation in eight out of nine markers during androgen-targeted treatment correlated with a faster time to clinical progression. Cox proportional hazards modeling and logistic regression analysis further confirmed that higher cfDNA methylation during treatment was significantly associated with clinical progression. Conclusion: Overall, our findings have revealed a novel methylated cfDNA marker panel that could aid in the clinical management of metastatic prostate cancer.

CRIS: complete reconstruction of immunoglobulin V-D-J sequences from RNA-seq data.

MOTIVATION: B cells display remarkable diversity in producing B-cell receptors through recombination of immunoglobulin (Ig) V-D-J genes. Somatic hypermutation (SHM) of immunoglobulin heavy chain variable (IGHV) genes are used as a prognostic marker in B-cell malignancies. Clinically, IGHV mutation status is determined by targeted Sanger sequencing which is a resource-intensive and low-throughput procedure. Here, we describe a bioinformatic pipeline, CRIS (Complete Reconstruction of Immunoglobulin IGHV-D-J Sequences) that uses RNA sequencing (RNA-seq) datasets to reconstruct IGHV-D-J sequences and determine IGHV SHM status. RESULTS: CRIS extracts RNA-seq reads aligned to Ig gene loci, performs assembly of Ig transcripts and aligns the resulting contigs to reference Ig sequences to enumerate and classify SHMs in the IGHV gene sequence. CRIS improves on existing tools that infer the B-cell receptor repertoire from RNA-seq data using a portion IGHV gene segment by de novo assembly. We show that the SHM status identified by CRIS using the entire IGHV gene segment is highly concordant with clinical classification in three independent chronic lymphocytic leukemia patient cohorts. AVAILABILITY AND IMPLEMENTATION: The CRIS pipeline is available under the MIT License from https://github.com/Rashedul/CRIS. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics Advances online.

Distinct DNA methylation patterns associated with treatment resistance in metastatic castration resistant prostate cancer.

Androgens are a major driver of prostate cancer (PCa) and continue to be a critical treatment target for advanced disease, which includes castration therapy and antiandrogens. However, resistance to these therapies leading to metastatic castration-resistant prostate cancer (mCRPC), and the emergence of treatment-induced neuroendocrine disease (tNEPC) remains an ongoing challenge. Instability of the DNA methylome is well established as a major hallmark of PCa development and progression. Therefore, investigating the dynamics of the methylation changes going from the castration sensitive to the tNEPC state would provide insights into novel mechanisms of resistance. Using an established xenograft model of CRPC, genome-wide methylation analysis was performed on cell lines representing various stages of PCa progression. We confirmed extensive methylation changes with the development of CRPC and tNEPC using this model. This included key genes and pathways associated with cellular differentiation and neurodevelopment. Combined analysis of methylation and gene expression changes further highlighted genes that could potentially serve as therapeutic targets. Furthermore, tNEPC-related methylation signals from this model were detectable in circulating cell free DNA (cfDNA) from mCRPC patients undergoing androgen-targeting therapies and were associated with a faster time to clinical progression. These potential biomarkers could help with identifying patients with aggressive disease.

Human placental cytotrophoblast epigenome dynamics over gestation and alterations in placental disease.

The human placenta and its specialized cytotrophoblasts rapidly develop, have a compressed lifespan, govern pregnancy outcomes, and program the offspring's health. Understanding the molecular underpinnings of these behaviors informs development and disease. Profiling the extraembryonic epigenome and transcriptome during the 2nd and 3rd trimesters revealed H3K9 trimethylation overlapping deeply DNA hypomethylated domains with reduced gene expression and compartment-specific patterns that illuminated their functions. Cytotrophoblast DNA methylation increased, and several key histone modifications decreased across the genome as pregnancy advanced. Cytotrophoblasts from severe preeclampsia had substantially increased H3K27 acetylation globally and at genes that are normally downregulated at term but upregulated in this syndrome. In addition, some cases had an immature pattern of H3K27ac peaks, and others showed evidence of accelerated aging, suggesting subtype-specific alterations in severe preeclampsia. Thus, the cytotrophoblast epigenome dramatically reprograms during pregnancy, placental disease is associated with failures in this process, and H3K27 hyperacetylation is a feature of severe preeclampsia.

Dynamics of the cell-free DNA methylome of metastatic prostate cancer during androgen-targeting treatment.

Aim: We examined methylation changes in cell-free DNA (cfDNA) in metastatic castration-resistant prostate cancer (mCRPC) during treatment. Patients & methods: Genome-wide methylation analysis of sequentially collected cfDNA samples derived from mCRPC patients undergoing androgen-targeting therapy was performed. Results: Alterations in methylation states of genes previously implicated in prostate cancer progression were observed and patients that maintained methylation changes throughout therapy tended to have a longer time to clinical progression. Importantly, we also report that markers associated with a highly aggressive form of the disease, neuroendocrine-CRPC, were associated with a faster time to clinical progression. Conclusion: Our findings highlight the potential of monitoring the cfDNA methylome during therapy in mCRPC, which may serve as predictive markers of response to androgen-targeting agents.

The Pathognomonic FOXL2 C134W Mutation Alters DNA-Binding Specificity.

The somatic missense point mutation c.402C>G (p.C134W) in the FOXL2 transcription factor is pathognomonic for adult-type granulosa cell tumors (AGCT) and a diagnostic marker for this tumor type. However, the molecular consequences of this mutation and its contribution to the mechanisms of AGCT pathogenesis remain unclear. To explore these mechanisms, we engineered V5-FOXL2WT- and V5-FOXL2C134W-inducible isogenic cell lines and performed chromatin immunoprecipitation sequencing and transcriptome profiling. FOXL2C134W associated with the majority of the FOXL2 wild-type DNA elements as well as a large collection of unique elements genome wide. This model enabled confirmation of altered DNA-binding specificity for FOXL2C134W and identification of unique targets of FOXL2C134W including SLC35F2, whose expression increased sensitivity to YM155. Our results suggest FOXL2C134W drives AGCT by altering the binding affinity of FOXL2-containing complexes to engage an oncogenic transcriptional program. SIGNIFICANCE: A mechanistic understanding of FOXL2C134W-induced regulatory state alterations drives discovery of a rationally designed therapeutic strategy.

Epigenomic programming in early fetal brain development.

Aim: To provide a comprehensive understanding of gene regulatory networks in the developing human brain and a foundation for interpreting pathogenic deregulation. Materials & methods: We generated reference epigenomes and transcriptomes of dissected brain regions and primary neural progenitor cells (NPCs) derived from cortical and ganglionic eminence tissues of four normal human fetuses. Results: Integration of these data across developmental stages revealed a directional increase in active regulatory states, transcription factor activities and gene transcription with developmental stage. Consistent with differences in their biology, NPCs derived from cortical and ganglionic eminence regions contained common, region specific, and gestational week specific regulatory states. Conclusion: We provide a high-resolution regulatory network for NPCs from different brain regions as a comprehensive reference for future studies.

Histone Chaperone Paralogs Have Redundant, Cooperative, and Divergent Functions in Yeast.

Gene duplications increase organismal robustness by providing freedom for gene divergence or by increasing gene dosage. The yeast histone chaperones Fpr3 and Fpr4 are paralogs that can assemble nucleosomes in vitro; however, the genomic locations they target and their functional relationship is poorly understood. We refined the yeast synthetic genetic array approach to enable the functional dissection of gene paralogs. Applying this method to Fpr3 and Fpr4 uncovered redundant, cooperative, and divergent functions. While Fpr3 is uniquely involved in chromosome segregation, Fpr3 and Fpr4 cooperate to regulate genes involved in polyphosphate metabolism and ribosome biogenesis. We find that the TRAMP5 RNA exosome is critical for fitness in Δfpr3Δfpr4 yeast and leverage this information to identify an important role for Fpr4 at the 5' ends of protein coding genes. Additionally, Fpr4 and TRAMP5 negatively regulate RNAs from the nontranscribed spacers of ribosomal DNA. Yeast lacking Fpr3 and Fpr4 exhibit a genome instability phenotype at the ribosomal DNA, which implies that these histone chaperones regulate chromatin structure and DNA access at this location. Taken together. we provide genetic and transcriptomic evidence that Fpr3 and Fpr4 operate separately, cooperatively, and redundantly to regulate a variety of chromatin environments.