Disruption of a GATA2-TAL1-ERG regulatory circuit promotes erythroid transition in healthy and leukemic stem cells.

Changes in gene regulation and expression govern orderly transitions from hematopoietic stem cells to terminally differentiated blood cell types. These transitions are disrupted during leukemic transformation, but knowledge of the gene regulatory changes underpinning this process is elusive. We hypothesized that identifying core gene regulatory networks in healthy hematopoietic and leukemic cells could provide insights into network alterations that perturb cell state transitions. A heptad of transcription factors (LYL1, TAL1, LMO2, FLI1, ERG, GATA2, and RUNX1) bind key hematopoietic genes in human CD34+ hematopoietic stem and progenitor cells (HSPCs) and have prognostic significance in acute myeloid leukemia (AML). These factors also form a densely interconnected circuit by binding combinatorially at their own, and each other's, regulatory elements. However, their mutual regulation during normal hematopoiesis and in AML cells, and how perturbation of their expression levels influences cell fate decisions remains unclear. In this study, we integrated bulk and single-cell data and found that the fully connected heptad circuit identified in healthy HSPCs persists, with only minor alterations in AML, and that chromatin accessibility at key heptad regulatory elements was predictive of cell identity in both healthy progenitors and leukemic cells. The heptad factors GATA2, TAL1, and ERG formed an integrated subcircuit that regulates stem cell-to-erythroid transition in both healthy and leukemic cells. Components of this triad could be manipulated to facilitate erythroid transition providing a proof of concept that such regulatory circuits can be harnessed to promote specific cell-type transitions and overcome dysregulated hematopoiesis.

Shared roles for Scl and Lyl1 in murine platelet production and function.

The stem cell leukemia (Scl or Tal1) protein forms part of a multimeric transcription factor complex required for normal megakaryopoiesis. However, unlike other members of this complex such as Gata1, Fli1, and Runx1, mutations of Scl have not been observed as a cause of inherited thrombocytopenia. We postulated that functional redundancy with its closely related family member, lymphoblastic leukemia 1 (Lyl1) might explain this observation. To determine whether Lyl1 can substitute for Scl in megakaryopoiesis, we examined the platelet phenotype of mice lacking 1 or both factors in megakaryocytes. Conditional Scl knockout (KO) mice crossed with transgenic mice expressing Cre recombinase under the control of the mouse platelet factor 4 (Pf4) promoter generated megakaryocytes with markedly reduced but not absent Scl These Pf4Sclc-KO mice had mild thrombocytopenia and subtle defects in platelet aggregation. However, Pf4Sclc-KO mice generated on an Lyl1-null background (double knockout [DKO] mice) had severe macrothrombocytopenia, abnormal megakaryocyte morphology, defective pro-platelet formation, and markedly impaired platelet aggregation. DKO megakaryocytes, but not single-knockout megakaryocytes, had reduced expression of Gata1, Fli1, Nfe2, and many other genes that cause inherited thrombocytopenia. These gene expression changes were significantly associated with shared Scl and Lyl1 E-box binding sites that were also enriched for Gata1, Ets, and Runx1 motifs. Thus, Scl and Lyl1 share functional roles in platelet production by regulating expression of partner proteins including Gata1. We propose that this functional redundancy provides one explanation for the absence of Scl and Lyl1 mutations in inherited thrombocytopenia.

RAB27A promotes melanoma cell invasion and metastasis via regulation of pro-invasive exosomes.

Despite recent advances in targeted and immune-based therapies, advanced stage melanoma remains a clinical challenge with a poor prognosis. Understanding the genes and cellular processes that drive progression and metastasis is critical for identifying new therapeutic strategies. Here, we found that the GTPase RAB27A was overexpressed in a subset of melanomas, which correlated with poor patient survival. Loss of RAB27A expression in melanoma cell lines inhibited 3D spheroid invasion and cell motility in vitro, and spontaneous metastasis in vivo. The reduced invasion phenotype was rescued by RAB27A-replete exosomes, but not RAB27A-knockdown exosomes, indicating that RAB27A is responsible for the generation of pro-invasive exosomes. Furthermore, while RAB27A loss did not alter the number of exosomes secreted, it did change exosome size and altered the composition and abundance of exosomal proteins, some of which are known to regulate cancer cell movement. Our data suggest that RAB27A promotes the biogenesis of a distinct pro-invasive exosome population. These findings support RAB27A as a key cancer regulator, as well as a potential prognostic marker and therapeutic target in melanoma.

Ancestry informative markers (AIMs) for Korean and other East Asian and South East Asian populations.

Inference of ancestry from biological evidence can provide investigative information, especially for unknown DNA donors. Although tools for predicting ancestry have been developing, ancestry research focusing on populations relevant for South Korea is not common and markers are seldom chosen specifically to differentiate Koreans from other East Asian and South East Asian populations. Here, we report ancestry informative markers (AIMs) for distinguishing six East/South East Asian regional populations: China, Japan, Indonesia, Philippines, South Korea and Thailand. Individual genotypes from these six populations were available in PanSNPdb: The HUGO Pan-Asian SNP Database. To select AIMs, we calculated four population divergence metrics for each SNP: Nei's FST, Rosenberg's Informativeness (In), the average absolute allele frequency difference between populations (δFmean) and the maximum allele frequency difference between populations (δFmax). Based on these values, we selected 100 single nucleotide polymorphisms (SNPs) for distinguishing the six populations, 13 of which exhibited large allele frequency differences between Koreans and non-Koreans. To assess the performance of the AIMs, we performed principal coordinates analysis (PCoA) on the individuals from all six populations and inferred ancestral population clusters using the STRUCTURE program. In conclusion, we found that the selected AIMs can be applied to distinguish the six East/South East Asian groups and we suggest the markers in this study will be helpful to establish ancestry panels for Korea and neighbouring populations.

BloodChIP: a database of comparative genome-wide transcription factor binding profiles in human blood cells.

The BloodChIP database (http://www.med.unsw.edu.au/CRCWeb.nsf/page/BloodChIP) supports exploration and visualization of combinatorial transcription factor (TF) binding at a particular locus in human CD34-positive and other normal and leukaemic cells or retrieval of target gene sets for user-defined combinations of TFs across one or more cell types. Increasing numbers of genome-wide TF binding profiles are being added to public repositories, and this trend is likely to continue. For the power of these data sets to be fully harnessed by experimental scientists, there is a need for these data to be placed in context and easily accessible for downstream applications. To this end, we have built a user-friendly database that has at its core the genome-wide binding profiles of seven key haematopoietic TFs in human stem/progenitor cells. These binding profiles are compared with binding profiles in normal differentiated and leukaemic cells. We have integrated these TF binding profiles with chromatin marks and expression data in normal and leukaemic cell fractions. All queries can be exported into external sites to construct TF-gene and protein-protein networks and to evaluate the association of genes with cellular processes and tissue expression.

Mesoderm-derived PDGFRA+ cells regulate the emergence of hematopoietic stem cells in the dorsal aorta.

Mouse haematopoietic stem cells (HSCs) first emerge at embryonic day 10.5 (E10.5), on the ventral surface of the dorsal aorta, by endothelial-to-haematopoietic transition. We investigated whether mesenchymal stem cells, which provide an essential niche for long-term HSCs (LT-HSCs) in the bone marrow, reside in the aorta-gonad-mesonephros and contribute to the development of the dorsal aorta and endothelial-to-haematopoietic transition. Here we show that mesoderm-derived PDGFRA+ stromal cells (Mesp1der PSCs) contribute to the haemogenic endothelium of the dorsal aorta and populate the E10.5-E11.5 aorta-gonad-mesonephros but by E13.5 were replaced by neural-crest-derived PSCs (Wnt1der PSCs). Co-aggregating non-haemogenic endothelial cells with Mesp1der PSCs but not Wnt1der PSCs resulted in activation of a haematopoietic transcriptional programme in endothelial cells and generation of LT-HSCs. Dose-dependent inhibition of PDGFRA or BMP, WNT and NOTCH signalling interrupted this reprogramming event. Together, aorta-gonad-mesonephros Mesp1der PSCs could potentially be harnessed to manufacture LT-HSCs from endothelium.

Induction of muscle-regenerative multipotent stem cells from human adipocytes by PDGF-AB and 5-azacytidine.

Terminally differentiated murine osteocytes and adipocytes can be reprogrammed using platelet-derived growth factor-AB and 5-azacytidine into multipotent stem cells with stromal cell characteristics. We have now optimized culture conditions to reprogram human adipocytes into induced multipotent stem (iMS) cells and characterized their molecular and functional properties. Although the basal transcriptomes of adipocyte-derived iMS cells and adipose tissue-derived mesenchymal stem cells were similar, there were changes in histone modifications and CpG methylation at cis-regulatory regions consistent with an epigenetic landscape that was primed for tissue development and differentiation. In a non-specific tissue injury xenograft model, iMS cells contributed directly to muscle, bone, cartilage, and blood vessels, with no evidence of teratogenic potential. In a cardiotoxin muscle injury model, iMS cells contributed specifically to satellite cells and myofibers without ectopic tissue formation. Together, human adipocyte-derived iMS cells regenerate tissues in a context-dependent manner without ectopic or neoplastic growth.

A Menin-MLL Inhibitor Induces Specific Chromatin Changes and Eradicates Disease in Models of MLL-Rearranged Leukemia.

Inhibition of the Menin (MEN1) and MLL (MLL1, KMT2A) interaction is a potential therapeutic strategy for MLL-rearranged (MLL-r) leukemia. Structure-based design yielded the potent, highly selective, and orally bioavailable small-molecule inhibitor VTP50469. Cell lines carrying MLL rearrangements were selectively responsive to VTP50469. VTP50469 displaced Menin from protein complexes and inhibited chromatin occupancy of MLL at select genes. Loss of MLL binding led to changes in gene expression, differentiation, and apoptosis. Patient-derived xenograft (PDX) models derived from patients with either MLL-r acute myeloid leukemia or MLL-r acute lymphoblastic leukemia (ALL) showed dramatic reductions of leukemia burden when treated with VTP50469. Multiple mice engrafted with MLL-r ALL remained disease free for more than 1 year after treatment. These data support rapid translation of this approach to clinical trials.

Mitogenomics phylogenetic relationships of the current sloth's genera and species (Bradypodidae and Megalonychidae).

We sequenced the complete mitogenome of 39 sloths (19 Bradypus variegatus, 4 B. tridactylus, 1 B. pygmaeus, 1 B. torquatus, 4 Choloepus didactylus, and 10 C. hoffmanni). A Bayesian tree (BI) indicated a temporal split between Bradypus and Choloepus around 31 million years ago (MYA, Oligocene) and the other major splits within each genera during the Miocene and Pliocene. A haplotype network (MJN) estimated a lower temporal split between the sloth genera (around 23.5 MYA). Both methods detected the ancestor of B. torquatus as the first to diverge within Bradypus (21 for BI and 19 MJN), followed by that of the ancestor of B. tridactylus. The split of B. pygmaeus from the common ancestor with B. variegatus was around 12 MYA (BI) or 4.3 MYA (MJN). The splits among the previous populations of B. variegatus began around 8 MYA (BI) or 3.6 MYA (MJN). The trans-Andean population was the first to diverge from the remaining cis-Andean populations of B. variegatus. The genetic differentiation of the trans-Andean B. variegatus population relative to the cis-Andean B. variegatus is similar to that found for different species of sloths. The mitogenomic analysis resolved the differentiation of C. hoffmanni from the C. didactylus individuals of the Guiana Shield. However, one C. didactylus from the Colombian Amazon specimen was inside the C. hoffmanni clade. This could be the first example of possible natural hybridization in the Amazon of both Choloepus taxa or the existence of un-differentiable phenotypes of these two species in some Amazonian areas.

OncoCis: annotation of cis-regulatory mutations in cancer.

Whole genome sequencing has enabled the identification of thousands of somatic mutations within non-coding genomic regions of individual cancer samples. However, identification of mutations that potentially alter gene regulation remains a major challenge. Here we present OncoCis, a new method that enables identification of potential cis-regulatory mutations using cell type-specific genome and epigenome-wide datasets along with matching gene expression data. We demonstrate that the use of cell type-specific information and gene expression can significantly reduce the number of candidate cis-regulatory mutations compared with existing tools designed for the annotation of cis-regulatory SNPs. The OncoCis webserver is freely accessible at https://powcs.med.unsw.edu.au/OncoCis/.