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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Characterization of the human thyroid epigenome.

The thyroid gland, necessary for normal human growth and development, functions as an essential regulator of metabolism by the production and secretion of appropriate levels of thyroid hormone. However, assessment of abnormal thyroid function may be challenging suggesting a more fundamental understanding of normal function is needed. One way to characterize normal gland function is to study the epigenome and resulting transcriptome within its constituent cells. This study generates the first published reference epigenomes for human thyroid from four individuals using ChIP-seq and RNA-seq. We profiled six histone modifications (H3K4me1, H3K4me3, H3K27ac, H3K36me3, H3K9me3, H3K27me3), identified chromatin states using a hidden Markov model, produced a novel quantitative metric for model selection and established epigenomic maps of 19 chromatin states. We found that epigenetic features characterizing promoters and transcription elongation tend to be more consistent than regions characterizing enhancers or Polycomb-repressed regions and that epigenetically active genes consistent across all epigenomes tend to have higher expression than those not marked as epigenetically active in all epigenomes. We also identified a set of 18 genes epigenetically active and consistently expressed in the thyroid that are likely highly relevant to thyroid function. Altogether, these epigenomes represent a powerful resource to develop a deeper understanding of the underlying molecular biology of thyroid function and provide contextual information of thyroid and human epigenomic data for comparison and integration into future studies.

Vitamin C induces specific demethylation of H3K9me2 in mouse embryonic stem cells via Kdm3a/b.

BACKGROUND: Histone methylation patterns regulate gene expression and are highly dynamic during development. The erasure of histone methylation is carried out by histone demethylase enzymes. We had previously shown that vitamin C enhances the activity of Tet enzymes in embryonic stem (ES) cells, leading to DNA demethylation and activation of germline genes. RESULTS: We report here that vitamin C induces a remarkably specific demethylation of histone H3 lysine 9 dimethylation (H3K9me2) in naïve ES cells. Vitamin C treatment reduces global levels of H3K9me2, but not other histone methylation marks analyzed, as measured by western blot, immunofluorescence and mass spectrometry. Vitamin C leads to widespread loss of H3K9me2 at large chromosomal domains as well as gene promoters and repeat elements. Vitamin C-induced loss of H3K9me2 occurs rapidly within 24 h and is reversible. Importantly, we found that the histone demethylases Kdm3a and Kdm3b are required for vitamin C-induced demethylation of H3K9me2. Moreover, we show that vitamin C-induced Kdm3a/b-mediated H3K9me2 demethylation and Tet-mediated DNA demethylation are independent processes at specific loci. Lastly, we document Kdm3a/b are partially required for the upregulation of germline genes by vitamin C. CONCLUSIONS: These results reveal a specific role for vitamin C in histone demethylation in ES cells and document that DNA methylation and H3K9me2 cooperate to silence germline genes in pluripotent cells.

Analysis of Normal Human Mammary Epigenomes Reveals Cell-Specific Active Enhancer States and Associated Transcription Factor Networks.

The normal adult human mammary gland is a continuous bilayered epithelial system. Bipotent and myoepithelial progenitors are prominent and unique components of the outer (basal) layer. The inner (luminal) layer includes both luminal-restricted progenitors and a phenotypically separable fraction that lacks progenitor activity. We now report an epigenomic comparison of these three subsets with one another, with their associated stromal cells, and with three immortalized, non-tumorigenic human mammary cell lines. Each genome-wide analysis contains profiles for six histone marks, methylated DNA, and RNA transcripts. Analysis of these datasets shows that each cell type has unique features, primarily within genomic regulatory regions, and that the cell lines group together. Analyses of the promoter and enhancer profiles place the luminal progenitors in between the basal cells and the non-progenitor luminal subset. Integrative analysis reveals networks of subset-specific transcription factors.

Nucleosome Density ChIP-Seq Identifies Distinct Chromatin Modification Signatures Associated with MNase Accessibility.

Nucleosome position, density, and post-translational modification are widely accepted components of mechanisms regulating DNA transcription but still incompletely understood. We present a modified native ChIP-seq method combined with an analytical framework that allows MNase accessibility to be integrated with histone modification profiles. Application of this methodology to the primitive (CD34+) subset of normal human cord blood cells enabled genomic regions enriched in one versus two nucleosomes marked by histone 3 lysine 4 trimethylation (H3K4me3) and/or histone 3 lysine 27 trimethylation (H3K27me3) to be associated with their transcriptional and DNA methylation states. From this analysis, we defined four classes of promoter-specific profiles and demonstrated that a majority of bivalent marked promoters are heterogeneously marked at a single-cell level in this primitive cell type. Interestingly, extension of this approach to human embryonic stem cells revealed an altered relationship between chromatin modification state and nucleosome content at promoters, suggesting developmental stage-specific organization of histone methylation states.

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.

Precursor States of Brain Tumor Initiating Cell Lines Are Predictive of Survival in Xenografts and Associated with Glioblastoma Subtypes.

In glioblastoma multiforme (GBM), brain-tumor-initiating cells (BTICs) with cancer stem cell characteristics have been identified and proposed as primordial cells responsible for disease initiation, recurrence, and therapeutic resistance. However, the extent to which individual, patient-derived BTIC lines reflect the heterogeneity of GBM remains poorly understood. Here we applied a stem cell biology approach and compared self-renewal, marker expression, label retention, and asymmetric cell division in 20 BTIC lines. Through cluster analysis, we identified two subgroups of BTIC lines with distinct precursor states, stem- or progenitor-like, predictive of survival after xenograft. Moreover, stem and progenitor transcriptomic signatures were identified, which showed a strong association with the proneural and mesenchymal subtypes, respectively, in the TCGA cohort. This study proposes a different framework for the study and use of BTIC lines and provides precursor biology insights into GBM.