Raymond A , Liu B , Liang H , Wei C , Guindani M , Lu Y , Liang S , St John LS , Molldrem J , Nagarajan L .

CA16672 NCI NIH HHS , P50-CA100632 NCI NIH HHS , P30 CA016672 NCI NIH HHS , P01 CA148600 NCI NIH HHS , P50 CA100632 NCI NIH HHS

	Figure 1. MIXL1 expression confers decreased sensitivity to doxorubicin in AML cells(A) Stable transfectants of U937 cells express MIXL1 at levels similar to those of endogenous MIXL1 in AML cell lines. MIXL1 was detected by probing 30 μg of whole cell lysates resolved on SDS-PAGE and transferred to PVDF membrane, with rabbit antibodies against N-terminal epitope of MIXL1 and with β-actin for a loading control [17]. (B) MIXL1 expression reduces sensitivity of U937 cells to doxorubicin. The cell lines were treated with 0–1.75 μM doxorubicin on day 0. Cell survival was measured at 24 hours by MTS assay as detailed in Materials and Methods. Absorbance of untreated cells was normalized to 1. Relative viability at varying concentrations of doxorubicin is denoted.
	Figure 2. Identification of direct MIXL1 transcriptional targets by ChIP-Sequencing(A) Venn diagram 1MIXL-Flag normalized to control-Flag and control-IgG (Set 1), 2MIXL-Flag normalized to control-Flag and control-IgG (Set 2), and 1MIXL-Flag and 2MIXL-Flag combined and normalized to control-Flag and control-IgG (Set 3). A total of 179 peaks shared by the three groups is denoted. (B) Pie chart depicting localization of MIXL1 in the human genome. Peaks were classified according to distance from the nearest transcribed gene using the following criteria: upstream was 5–25 kbp 5′ of the transcription start site, promoter was 0–5 kbp upstream of the transcription start site, body was between the transcription start site and end, TSE was 0–5 kbp downstream of the transcriptional end, downstream was 5–25 kbp downstream of the transcriptional end, and distant peaks were those not allocated to a gene. Note that the majority of peaks (64%) localized to gene promoters. (C) ChIP of five candidate peaks identified by ChIP-Se1. FLAG antibodies were used and ChIP-Seq (EIF1, c-REL, SLC39A13, SMYD5, and ZP3) showed specific MIXL1 binding to both 1MIXL and 2MIXL clones by ChIP normal mouse IgG served as control. Error bars represent standard deviation between triplicates. (D) The most common motif in the ChIP-seq peaks are Zinc-Finger binding sites. Motif1 and Motif2 were the two most statistically significant generated using the Multiple EM for Motif Elicitation [MEME] against the peak regions identified in the ChIP-seq analysis. Both motifs are C/G-heavy regions with similarity to known zinc finger motifs.
	Figure 3. MIXL1 up regulates c-REL expression to enhance anti apoptotic gene transcription(A) MIXL1-expressing clones show enhanced transcript levels for c-REL, BCL2A1, and BCL2L1. Quantitative RT-PCR results show the differences in c-REL, BCL2A1, and BCL2L1 expression levels between the U937 control, 1MIXL, and 2MIXL cells. Expression was normalized to 18S rRNA transcript levels. Error bars represent standard deviation between triplicates. *p < 0.05. (B) ChIP localizes endogenous MIXL1 to c-REL promoter in KG1 cells. Quantitative genomic PCR analysis shows specific enrichment of endogenous MIXL1 immunoprecipitated with either N-terminal or C-terminal MIXL1 antibodies on the c-REL promoter whereas an internal locus within the c-REL gene showed no MIXL1 occupancy. Error bars represent standard deviation between triplicates. (C) Knockdown of MIXL1 decreased while enforced expression of c-REL increased c-REL, BCL2A1, and BCL2L1 transcript levels. MIXL1 shRNA lentivirus and c-REL retrovirus were transduced into KG1 cells. RT-qPCR was performed in triplicate on RNAs isolated 48 hours after transduction. Expression was normalized to 18S rRNA levels, and error bars represent standard deviation between triplicates. (D) c-REL over-expression rescues MIXL1 knockdown–mediated growth arrest in KG1 cells. Growth was measured by MTS assay every 24 hours over a 4-day period in KG1 cells transduced with MIXL1 shRNA lentivirus and c-REL retrovirus. Absorbance was normalized to that of a non-transfected control sample. *p < 0.05.
	Figure 4. MIXL1 binds to the c-REL promoter(A) c-REL promoter peak region identified by ChIP-Seq, as generated by the University of California, Santa Cruz, genome browser is shown. The location and size of each promoter fragment used for the luciferase reporter assay is displayed underneath. (B) MIXL1 binds to a 550-bp region within the c-REL promoter. Regions of the DNA depicted in 4A were cloned into the reporter vector pBV-Luc luciferase, which were then transiently co-transfected into HEK293T cells with MIXL1, MIXL1 Homeobox-less, or empty expression vector. Equal amount of Renilla luciferase co-transfected with the reporter constructs allowed normalization. Luciferase activity of each combination was tested in triplicate after 48 hours. Error bars represent standard deviation between triplicates. (C) MZF1 binds to the same locus as MIXL1 on c-REL promoter. Quantitative PCR analysis of the identified c-REL promoter region and c-REL intron control region compared the abundance of each genomic locus immunoprecipitated by either IgG or MZF1 antibodies, normalized to a standard curve. Error bars represent standard deviation between triplicates. *p < 0.05. (D) Model depicting a potential MIXL1-TBX-MZF1 multiprotein complex activating c-REL transcription.
	Figure 5. High MIXL1 expression denotes a distinct subset of AML(A) Mutually exclusive expression of MIXL1 and HOXA9 in distinct FAB subsets. The RNA-seq data are publicly available from the TCGA website (https://tcga-data.nci.nih.gov/tcga/). 177 samples of Acute Myeloid Leukemia classified by leukemia French American British morphology code (FAB) and a total of 20319 genes with expression values in the RPKM format were included. The data were quantile normalized using the normalize Quantiles function from the limma package. The expression levels of MIXL1 and HOXA9 of the 177 samples from 8 FAB categories were plotted in the heatmap using the heatmap.2 function in the gplots package of R 3.1.1. Within each FAB category, the samples were ordered according to the expression value of the MIXL1 gene. (B) MIXL1 upregulation identifies a non-overlapping AML subset from those expressing CDX2, HOXA9, or HLX. TCGA AML patient dataset was queried for alterations in expression as determined by RNA-Seq across 166 AML cases through the cBioPortal database. Each column represents a case of AML. MIXL1 is amplified or upregulated in 13% of the total AML cases. Note the predominantly non-overlapping expression patterns of MIXL1, CDX2, HOXA9, and HLX. (C) Seventy-six percent (16/21) of MIXL1-expressing cases in TCGA AML dataset harbored somatic mutations common in AML (NPM1, FLT3, DNMT3A, IDH1, RUNX1 JAK3, and TP53). Each column represents a case. (D) Relapse-free survival of MIXL1-expressing cases is lower than that of non–MIXL1-expressing cases. TCGA AML cases were separated into two groups: MIXL1-expressing (increase in expression or amplified) and non–MIXL1-expressing. Relapse-free survival was then compared between the two groups using the Kaplan-Meier method.
	Figure 6. BMP4 induced MIXL1, an important survival axis and therapeutic target in AML(A) MIXL1 expression increased 2-fold in CD34+ HSPCs treated with BMP4. CD34+ HSPCs from three cord blood donors (unique donor number 47, 51, 60) were cultured with either 50 ng/ml BMP4 or 2 ng/ml TGF-β1 for 2 hours. MIXL1 transcript levels were quantified by RT-qPCR using 18S rRNA as normalization control. Error bars represent standard deviation between triplicates. *p < 0.05. (B) LDN-193189 at 3 μM was cytotoxic to OCI-AML2, ML3, KG1, and K562 cells but not U937 and HL60 cells. Each cell line was treated with vehicle or 3 μM LDN-193189 on day 0, and viability was measured every 24 hours by MTS assay in triplicate. (C) OCI-AML2, ML3, KG1, and K562 were sensitive to 200 nM LDN-193189, while the non–MIXL1-expressing lines U937 and HL60 were unaffected. Each cell line was treated in triplicate with 0–700 nM LDN-193189, with the drug or control medium replenished every 24 hours, for 4 days. Viability on day 4 was assayed by MTS assay. Absorbance was normalized to that of control samples treated with vehicle only.
	Figure 7. BMP4 induced MIXL1 an important survival axis and therapeutic target in AML(A) Canonical BMP signaling pathway. Upon BMP binding to type II receptor (BMPRII), type I receptor BMPRI is phosphorylated. Activated BMPRI phosphorylates SMAD1 or 5 which heterodimerizes with SMAD4 to up regulate MIXL1 expression. (B) BMPR1 mediated signal is hypothesized to induce endogenous MIXL1 expression in KG1, OCI-ML2, ML3 and K562 cell lines. These cells are sensitive to BMPR1 inhibitor LDN-193189. Enforced expression of MIXL1 in U937 cells 1MIXL1 is independent of BMP signaling and therefore insensitive to LDN-193189. In HL60 and U937 cells which are insensitive to LDN-193189, the pathway may be absent or modified.

Alfa, Murine and HIV-based retroviral vectors for in vitro and in vivo gene transfer. 2006, Pubmed

Alfa, Murine and HIV-based retroviral vectors for in vitro and in vivo gene transfer. 2006, Pubmed
Andreeff, HOX expression patterns identify a common signature for favorable AML. 2008, Pubmed
Argiropoulos, Hox genes in hematopoiesis and leukemogenesis. 2007, Pubmed
Ashburner, Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. 2000, Pubmed
Bailey, Fitting a mixture model by expectation maximization to discover motifs in biopolymers. 1994, Pubmed
Bailey, Methods and statistics for combining motif match scores. 1998, Pubmed
Bhatia, Bone morphogenetic proteins regulate the developmental program of human hematopoietic stem cells. 1999, Pubmed
Boer, Overexpression of the nucleoporin CAN/NUP214 induces growth arrest, nucleocytoplasmic transport defects, and apoptosis. 1998, Pubmed
Boergermann, Dorsomorphin and LDN-193189 inhibit BMP-mediated Smad, p38 and Akt signalling in C2C12 cells. 2010, Pubmed
Bogenberger, BCL-2 family proteins as 5-Azacytidine-sensitizing targets and determinants of response in myeloid malignancies. 2014, Pubmed
Bolstad, A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. 2003, Pubmed
Burnett, Therapeutic advances in acute myeloid leukemia. 2011, Pubmed
Caslini, The amino terminus of the mixed lineage leukemia protein (MLL) promotes cell cycle arrest and monocytic differentiation. 2000, Pubmed
Caudell, Expression of a CALM-AF10 fusion gene leads to Hoxa cluster overexpression and acute leukemia in transgenic mice. 2007, Pubmed
Cerami, The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. 2012, Pubmed
Chadee, Increased Ser-10 phosphorylation of histone H3 in mitogen-stimulated and oncogene-transformed mouse fibroblasts. 1999, Pubmed
Chen, The Rel/NF-kappaB family directly activates expression of the apoptosis inhibitor Bcl-x(L). 2000, Pubmed
Coenen, Pediatric acute myeloid leukemia with t(8;16)(p11;p13), a distinct clinical and biological entity: a collaborative study by the International-Berlin-Frankfurt-Munster AML-study group. 2013, Pubmed
Cuny, Structure-activity relationship study of bone morphogenetic protein (BMP) signaling inhibitors. 2008, Pubmed
Drakos, Differential expression of the human MIXL1 gene product in non-Hodgkin and Hodgkin lymphomas. 2007, Pubmed
Döhner, Diagnosis and management of acute myeloid leukemia in adults: recommendations from an international expert panel, on behalf of the European LeukemiaNet. 2010, Pubmed
Eklund, The role of HOX genes in malignant myeloid disease. 2007, Pubmed
Eklund, The role of Hox proteins in leukemogenesis: insights into key regulatory events in hematopoiesis. 2011, Pubmed
Feng, Specificity and versatility in tgf-beta signaling through Smads. 2005, Pubmed
Fields, Activin receptor antagonists for cancer-related anemia and bone disease. 2013, Pubmed
Fotinos, Bone morphogenetic protein-focused strategies to induce cytotoxicity in lung cancer cells. 2014, Pubmed
Fröhling, HOX gene regulation in acute myeloid leukemia: CDX marks the spot? 2007, Pubmed
Gao, Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. 2013, Pubmed
Garg, GATA4 mutations cause human congenital heart defects and reveal an interaction with TBX5. 2003, Pubmed , Xenbase
Gilmore, The c-Rel transcription factor and B-cell proliferation: a deal with the devil. 2004, Pubmed
Gilmore, The c-Rel Transcription Factor in Development and Disease. 2011, Pubmed
Glaser, Enforced expression of the homeobox gene Mixl1 impairs hematopoietic differentiation and results in acute myeloid leukemia. 2006, Pubmed , Xenbase
Golub, Molecular classification of cancer: class discovery and class prediction by gene expression monitoring. 1999, Pubmed
Gruber, An Inv(16)(p13.3q24.3)-encoded CBFA2T3-GLIS2 fusion protein defines an aggressive subtype of pediatric acute megakaryoblastic leukemia. 2012, Pubmed
Gu, Phosphotyrosine profiling identifies the KG-1 cell line as a model for the study of FGFR1 fusions in acute myeloid leukemia. 2006, Pubmed
Guo, A human Mix-like homeobox gene MIXL shows functional similarity to Xenopus Mix.1. 2002, Pubmed , Xenbase
Haferlach, Another piece of the AML puzzle. 2013, Pubmed
Hromas, Hematopoietic transcriptional regulation by the myeloid zinc finger gene, MZF-1. 1996, Pubmed
Hu, Stable expression of a dominant negative mutant of CCAAT binding factor/NF-Y in mouse fibroblast cells resulting in retardation of cell growth and inhibition of transcription of various cellular genes. 2000, Pubmed
Huang, Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. 2009, Pubmed
Huang, Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. 2009, Pubmed
Huang, Extracting biological meaning from large gene lists with DAVID. 2009, Pubmed
Hughes, Clonal architecture of secondary acute myeloid leukemia defined by single-cell sequencing. 2014, Pubmed
Hutton, Bone morphogenetic protein 4 contributes to the maintenance of primitive cord blood hematopoietic progenitors in an ex vivo stroma-noncontact co-culture system. 2006, Pubmed
Iancu-Rubin, Stromal cell-mediated inhibition of erythropoiesis can be attenuated by Sotatercept (ACE-011), an activin receptor type II ligand trap. 2013, Pubmed
Kandilci, Reintroduction of CEBPA in MN1-overexpressing hematopoietic cells prevents their hyperproliferation and restores myeloid differentiation. 2009, Pubmed
Kawahara, H2.0-like homeobox regulates early hematopoiesis and promotes acute myeloid leukemia. 2012, Pubmed
Khurana, A novel role of BMP4 in adult hematopoietic stem and progenitor cell homing via Smad independent regulation of integrin-α4 expression. 2013, Pubmed
Kroon, Hoxa9 transforms primary bone marrow cells through specific collaboration with Meis1a but not Pbx1b. 1998, Pubmed
Laperrousaz, Primitive CML cell expansion relies on abnormal levels of BMPs provided by the niche and on BMPRIb overexpression. 2013, Pubmed
Lawrence, Frequent co-expression of the HOXA9 and MEIS1 homeobox genes in human myeloid leukemias. 1999, Pubmed
Ley, A pilot study of high-throughput, sequence-based mutational profiling of primary human acute myeloid leukemia cell genomes. 2003, Pubmed
Ley, DNA sequencing of a cytogenetically normal acute myeloid leukaemia genome. 2008, Pubmed
Ley, Genomic and epigenomic landscapes of adult de novo acute myeloid leukemia. 2013, Pubmed
Li, Model-based analysis of oligonucleotide arrays: expression index computation and outlier detection. 2001, Pubmed
Liang, SSBP2, a candidate tumor suppressor gene, induces growth arrest and differentiation of myeloid leukemia cells. 2005, Pubmed , Xenbase
Link, Identification of a novel TP53 cancer susceptibility mutation through whole-genome sequencing of a patient with therapy-related AML. 2011, Pubmed
Mardis, Recurring mutations found by sequencing an acute myeloid leukemia genome. 2009, Pubmed
Mead, Cloning of Mix-related homeodomain proteins using fast retrieval of gel shift activities, (FROGS), a technique for the isolation of DNA-binding proteins. 1998, Pubmed , Xenbase
Mead, BMP-4-responsive regulation of dorsal-ventral patterning by the homeobox protein Mix.1. 1996, Pubmed , Xenbase
Metcalf, The preleukemic state of mice reconstituted with Mixl1-transduced marrow cells. 2007, Pubmed
Mizutani, Cell type-specific target selection by combinatorial binding of Smad2/3 proteins and hepatocyte nuclear factor 4alpha in HepG2 cells. 2011, Pubmed
Mullighan, Pediatric acute myeloid leukemia with NPM1 mutations is characterized by a gene expression profile with dysregulated HOX gene expression distinct from MLL-rearranged leukemias. 2007, Pubmed
Métais, BCL2A1a over-expression in murine hematopoietic stem and progenitor cells decreases apoptosis and results in hematopoietic transformation. 2012, Pubmed
Nakamura, Cooperative activation of Hoxa and Pbx1-related genes in murine myeloid leukaemias. 1996, Pubmed
Pereira, Pdgfrα and Flk1 are direct target genes of Mixl1 in differentiating embryonic stem cells. 2012, Pubmed
Pereira, Brachyury and related Tbx proteins interact with the Mixl1 homeodomain protein and negatively regulate Mixl1 transcriptional activity. 2011, Pubmed
Ravasi, An atlas of combinatorial transcriptional regulation in mouse and man. 2010, Pubmed
Rawat, Overexpression of CDX2 perturbs HOX gene expression in murine progenitors depending on its N-terminal domain and is closely correlated with deregulated HOX gene expression in human acute myeloid leukemia. 2008, Pubmed
Rosa, Mix.1, a homeobox mRNA inducible by mesoderm inducers, is expressed mostly in the presumptive endodermal cells of Xenopus embryos. 1989, Pubmed , Xenbase
Sanvitale, A new class of small molecule inhibitor of BMP signaling. 2013, Pubmed
Scholl, The homeobox gene CDX2 is aberrantly expressed in most cases of acute myeloid leukemia and promotes leukemogenesis. 2007, Pubmed
Sergeeva, An anti-PR1/HLA-A2 T-cell receptor-like antibody mediates complement-dependent cytotoxicity against acute myeloid leukemia progenitor cells. 2011, Pubmed
Silliman, Alternative splicing in wild-type AF10 and CALM cDNAs and in AF10-CALM and CALM-AF10 fusion cDNAs produced by the t(10;11)(p13-14;q14-q21) suggests a potential role for truncated AF10 polypeptides. 1998, Pubmed
Steinbicker, Inhibition of bone morphogenetic protein signaling attenuates anemia associated with inflammation. 2011, Pubmed
Trompouki, Lineage regulators direct BMP and Wnt pathways to cell-specific programs during differentiation and regeneration. 2011, Pubmed
Vassiliou, Mutant nucleophosmin and cooperating pathways drive leukemia initiation and progression in mice. 2011, Pubmed
Vogt, The specificities of small molecule inhibitors of the TGFß and BMP pathways. 2011, Pubmed
Walter, Clonal architecture of secondary acute myeloid leukemia. 2012, Pubmed
Wang, Persistent transactivation by meis1 replaces hox function in myeloid leukemogenesis models: evidence for co-occupancy of meis1-pbx and hox-pbx complexes on promoters of leukemia-associated genes. 2006, Pubmed
Willey, Acceleration of mesoderm development and expansion of hematopoietic progenitors in differentiating ES cells by the mouse Mix-like homeodomain transcription factor. 2006, Pubmed
Wilson, Cooperative dimerization of paired class homeo domains on DNA. 1993, Pubmed
Wilson, Conservation and diversification in homeodomain-DNA interactions: a comparative genetic analysis. 1996, Pubmed
Xi, A poised chromatin platform for TGF-β access to master regulators. 2011, Pubmed
Yamazaki, Nonmyelinating Schwann cells maintain hematopoietic stem cell hibernation in the bone marrow niche. 2011, Pubmed
Yu, BMP type I receptor inhibition reduces heterotopic [corrected] ossification. 2008, Pubmed
Zhang, Transcriptional activation by the Mixl1 homeodomain protein in differentiating mouse embryonic stem cells. 2009, Pubmed
Zhang, Model-based analysis of ChIP-Seq (MACS). 2008, Pubmed
Zhao, Noncanonical NF-κB signaling regulates hematopoietic stem cell self-renewal and microenvironment interactions. 2012, Pubmed